Protein kinases are enzymes that act as molecular switches for signaling pathways in cells. Most kinases have an “on”-state that promotes signaling, and an “off”-state that does not, and these two states have different shapes or conformations. A number of diseases are associated with kinases not functioning properly, so kinases are important targets for researchers working in drug development.
In general, when a protein kinase is switched on, ATP (a molecule that provides cells with energy) binds to the active site of the kinase, and the kinase facilitates the transfer of a phosphate group from ATP to a protein. This means that it is possible to inhibit a kinase by getting a small molecule to bind to the active site instead of ATP. Some of these inhibitors can bind to just one state, but others can bind to both the on-state and the off-state.
To develop better drugs it is important to understand how kinase conformations relate to signaling and how inhibitor molecules select for these conformations. Now, in eLife, Natalie Ahn (University of Colorado, Boulder) and colleagues – including Jake Anderson as first author – report the results of experiments on a kinase called ERK2 that is involved in many signaling pathways in cells (Anderson et al., 2023). ERK2 is of particular interest because it is a key player in the MAPK signaling pathway that is often misregulated in human cancers (Roskoski, 2019; Lavoie et al., 2020).
In many kinases the on-state and the off-state have very different confirmations, but that it not the case for ERK2 (Zhang et al., 1994; Canagarajah et al., 1997). Indeed, the difference is so subtle that it is difficult to distinguish between the two states with X-ray crystallography, a technique that is widely used to determine the static structures of biomolecules. Fortunately, the dynamics of the two states are different, so Anderson et al. were able to use two established techniques for the study of dynamics – nuclear magnetic resonance and hydrogen-deuterium exchange mass spectrometry – to examine the dynamics of ERK2. They found that the the on-state corresponds to at least two conformational substates (named R and L) that interconvert on a timescale of milliseconds.
First, the researchers – who are based the University of Colorado in Boulder, the University of Colorado Anschutz Medical Center and Genentech – tested the impact of three inhibitors on the ERK2 kinase. They found ERK2 is predominantly in the R-state after it binds to an inhibitor called VTX11e, and similarly when it binds to BVD523. However, when ERK2 binds to an inhibitor called GDC0994, it remains in both the R-state and the L-state. What chemical motifs in VTX11e and BVD523 select for the ERK2 R-state? In an efort to answer this question the researchers tested a series of 17 ERK2 inhibitors that were similar to GDC0994 in various ways (see Figure 3 of Anderson et al., 2023). For three of these inhibitors the ERK2 kinase remained in both the R-state and the L-state, but in thirteen cases it selected the R-state (similar to the behavior observed for VTX11e and BVD523.) The remaining inhibitor did not fall neatly into either category.
To better understand the mechanism underlying such selection, Anderson et al. determined the three-dimensional structure of several inhibitors bound to ERK2. They observed that the large bulky chemical group on the right side of GDC0994 may have prevented conformational selection of the R-state, as it did for other inhibitors. The three inhibitors for which the ERK2 kinase remained in both the R-state and the L-state also had similar bulky chemical groups on their right sides, whereas the thirteen inhibitors that selected the R-state did not (Figure 1A). This observation prompted the researchers to ask what would happen if the same bulky group was added to the right side of an inhibitor that would normally select the R-state of ERK2. Luckily, an inhibitor called ATG017 has such a structure, and Anderson et al. found that it was able to bind to both the R-state and the L-state of ERK2, just as GDC0994 does (Figure 1B). This is consistent with previous work on ATG017 (also known as AZD0364) by AstraZeneca (Ward et al., 2019; Flemington et al., 2021).
(A) BVD523, VTX11e and 13 of the 17 ERK2 inhibitors studied by Anderson et al. select for the R-state when binding to the ERK2 kinase. GCD0994, ATG017 and 3 of the 17 inhibitors do not select for …
The work of Anderson et al. also sheds light on how the spatial configuration of atoms in ERK2, the inhibitors and even the surrounding solvent contribute to promoting these conformational preferences. For example, when an R-state-selective inhibitor binds to ERK2, movements in the upper lobe of the kinase bring the residues required for catalytic activity closer together. Furthermore, ERK2 is regulated in an allosteric manner (that is, the regulatory site and the active site are different), and Anderson et al. show that the residues involved in allosteric regulation re-arrange their positions in order to relay an allosteric signal from the regulatory site to the active site. Additionally, certain R-state-selective inhibitors influenced the flexibility of the activation loop of ERK2, which prompted Anderson et al. to suggest that these inhibitors might also influence how ERK2 is regulated by other proteins.
In summary, Anderson et al. show how small changes in ERK2 inhibitors might lead to larger changes in ERK2 dynamics. Specifically, their work provides a roadmap for how to engineer and study compounds that preferentially bind the ERK2 kinase in the R-state. Moreover, it will be interesting to explore if comparable phenomena occur with other similar kinases.
© 2023, Outhwaite and Seeliger
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
The spike protein is essential to the SARS-CoV-2 virus life cycle, facilitating virus entry and mediating viral-host membrane fusion. The spike contains a fatty acid (FA) binding site between every two neighbouring receptor-binding domains. This site is coupled to key regions in the protein, but the impact of glycans on these allosteric effects has not been investigated. Using dynamical nonequilibrium molecular dynamics (D-NEMD) simulations, we explore the allosteric effects of the FA site in the fully glycosylated spike of the SARS-CoV-2 ancestral variant. Our results identify the allosteric networks connecting the FA site to functionally important regions in the protein, including the receptor-binding motif, an antigenic supersite in the N-terminal domain, the fusion peptide region, and another allosteric site known to bind heme and biliverdin. The networks identified here highlight the complexity of the allosteric modulation in this protein and reveal a striking and unexpected link between different allosteric sites. Comparison of the FA site connections from D-NEMD in the glycosylated and non-glycosylated spike revealed that glycans do not qualitatively change the internal allosteric pathways but can facilitate the transmission of the structural changes within and between subunits.
Deep Mutational Scanning (DMS) is an emerging method to systematically test the functional consequences of thousands of sequence changes to a protein target in a single experiment. Because of its utility in interpreting both human variant effects and protein structure-function relationships, it holds substantial promise to improve drug discovery and clinical development. However, applications in this domain require improved experimental and analytical methods. To address this need, we report novel DMS methods to precisely and quantitatively interrogate disease-relevant mechanisms, protein-ligand interactions, and assess predicted response to drug treatment. Using these methods, we performed a DMS of the melanocortin-4 receptor (MC4R), a G-protein-coupled receptor (GPCR) implicated in obesity and an active target of drug development efforts. We assessed the effects of >6600 single amino acid substitutions on MC4R’s function across 18 distinct experimental conditions, resulting in >20 million unique measurements. From this, we identified variants that have unique effects on MC4R-mediated Gαs- and Gαq-signaling pathways, which could be used to design drugs that selectively bias MC4R’s activity. We also identified pathogenic variants that are likely amenable to a corrector therapy. Finally, we functionally characterized structural relationships that distinguish the binding of peptide versus small molecule ligands, which could guide compound optimization. Collectively, these results demonstrate that DMS is a powerful method to empower drug discovery and development.