Conformation selection by ATP-competitive inhibitors and allosteric communication in ERK2
eLife assessment
This fundamental study provides compelling evidence to explain how chemical variations within a set of kinase inhibitors drive the selection of specific Erk2 conformations. Conformational selection plays a critical role in targeting medically relevant kinases such as Erk2 and the findings reported here open new avenues for designing small molecule inhibitors that block the active site while also steering the population of the enzyme into active or inactive conformations. Since protein dynamics and conformational ensembles are essential for enzyme function, this work will be of broad interest to those working in drug development, signal transduction, and enzymology.
https://doi.org/10.7554/eLife.91507.3.sa0Fundamental: Findings that substantially advance our understanding of major research questions
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Compelling: Evidence that features methods, data and analyses more rigorous than the current state-of-the-art
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
Activation of the extracellular signal-regulated kinase-2 (ERK2) by phosphorylation has been shown to involve changes in protein dynamics, as determined by hydrogen-deuterium exchange mass spectrometry (HDX-MS) and NMR relaxation dispersion measurements. These can be described by a global exchange between two conformational states of the active kinase, named ‘L’ and ‘R,’ where R is associated with a catalytically productive ATP-binding mode. An ATP-competitive ERK1/2 inhibitor, Vertex-11e, has properties of conformation selection for the R-state, revealing movements of the activation loop that are allosterically coupled to the kinase active site. However, the features of inhibitors important for R-state selection are unknown. Here, we survey a panel of ATP-competitive ERK inhibitors using HDX-MS and NMR and identify 14 new molecules with properties of R-state selection. They reveal effects propagated to distal regions in the P+1 and helix αF segments surrounding the activation loop, as well as helix αL16. Crystal structures of inhibitor complexes with ERK2 reveal systematic shifts in the Gly loop and helix αC, mediated by a Tyr-Tyr ring stacking interaction and the conserved Lys-Glu salt bridge. The findings suggest a model for the R-state involving small movements in the N-lobe that promote compactness within the kinase active site and alter mobility surrounding the activation loop. Such properties of conformation selection might be exploited to modulate the protein docking interface used by ERK substrates and effectors.
Introduction
The extracellular signal-regulated protein kinases, ERK1/2, are key components of the MAP kinase signaling pathway controlling cell proliferation, survival, morphology, and migration (Lavoie et al., 2020; Roskoski, 2012; Plotnikov et al., 2011). Due to the prevalence of oncogenic mutations in Ras, RAF, and MEK, components of this pathway have been targeted for the development of successful anti-cancer therapies (Halle and Johnson, 2021; Tangella et al., 2021; Canon et al., 2019; Ryan et al., 2020). Although inhibitors targeting ERK1/2 have not advanced to clinical use, clinical trials are ongoing for several, including ulixertinib/BVD523 (Biomed Valley Diagnostics), ravoxertinib/GDC0994 (Genentech), and tizakertib/ATG017/AZD0364 (Antengene), for the treatment of cancers with activating mutations in KRAS, BRAF, NRAS, MEK1/2, and ERK1/2 (Grierson et al., 2023; Goodwin et al., 2023; Wu et al., 2021; Varga et al., 2020; Kong et al., 2023; Chen et al., 2021). Importantly, reversible ATP-competitive inhibitors targeting ERK1/2 have been shown to be effective in cancer cells and tumors with acquired resistance to BRAF and MEK inhibitors (Ryan et al., 2015; Morris et al., 2013; Hatzivassiliou et al., 2012; Roskoski, 2019; Pan et al., 2022). Therefore, future anticancer drug development would benefit from a more complete understanding of the properties of ERK inhibitors.
Many protein kinases undergo relatively large conformational changes when they switch from inactive to active states. Such changes lead to the repositioning of conserved active site residues, including a Lys-Glu salt bridge in the N-terminal lobe between strand β3 and helix αC that hydrogen bonds with nucleotide phosphate, an Asp residue in the DFG motif that coordinates Mg2+, and an Asp in the HRD motif that catalyzes base-assisted phosphoryl transfer to the substrate hydroxyl acceptor (Taylor et al., 2019; Amatya et al., 2019). By contrast, ERK2 is unusual in that only small conformational changes occur within the active site between the inactive, unphosphorylated (0P) and the active, dual phosphorylated (2P) forms of the kinase (Figure 1). X-ray structures show that dual phosphorylation of ERK2 at T183 and Y185 remodels the activation loop by forming phosphate salt bridges with several Arg residues (Zhang et al., 1994; Canagarajah et al., 1997). This accommodates the S/T-P recognition motif common to ERK substrates and exposes a pocket for binding a hydrophobic docking motif found in substrates, scaffolds, and other effectors (Lee et al., 2004). But within the ATP binding site, the conserved residues that form the K52-E69 salt bridge, the DFG metal binding site (D165), and the HRD catalytic base (D147) are largely overlapping (Figure 1C). Thus, the changes in the active site that accompany ERK2 activation are small and difficult to discern crystallographically.
Solution measurements revealed that the activation of ERK2 is accompanied by changes in protein dynamics. NMR Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion measurements showed that dual phosphorylation results in a global exchange behavior that can be modeled by two conformational states, named ‘R’ and ‘L,’ which exchange on a millisecond timescale (Xiao et al., 2014; Xiao et al., 2015; Iverson et al., 2020). Thus, 2P-ERK2 can be characterized by the interconversion of at least two conformers. Whereas the L-state resembles that seen with 0P-ERK2, the R-state appears only in 2P-ERK2, where it exchanges with an R:L population ratio of 80:20 at 25 °C. Crystal structures of nucleotide-bound kinase show significant differences between 0P- and 2P-ERK2. In 2P-ERK2, the ATP analog, AMP-PNP, is positioned close to the catalytic base and in-line for nucleophilic attack on the Pγ phosphate (Pegram et al., 2019; Lechtenberg et al., 2017). By contrast, ATP bound to 0P-ERK2 is markedly distorted, moving Pγ away from the base by more than 6 Å (Zhang et al., 2012; Smorodinsky-Atias et al., 2016). Thus, despite the structural similarity between the apoenzyme forms, 2P-ERK2 promotes ATP binding interactions that are more productive for catalytic function. The unique appearance of the R-state in 2P-ERK2 has been proposed to confer the ability to form a Michaelis complex poised for phosphoryl transfer, while L may represent a state that promotes ADP product release (Pegram et al., 2019).
Interestingly, ATP competitive inhibitors of ERK show properties of conformation selection, where ligand binding traps conformational states in 2P-ERK2 that normally undergo reversible exchange in the apoenzyme. Binding of the inhibitor, Vertex-11e (VTX11e), to 2P-ERK2 shifts the L↔R equilibrium completely to the R-state, while the inhibitor SCH772984 shifts the equilibrium completely to the L-state (Rudolph et al., 2015; Pegram et al., 2019). Conformation selection is absent in the inhibitor, GDC0994, where R and L retain a population ratio similar to that seen in the apoenzyme (Pegram et al., 2019). This selective recognition could explain why the binding of VTX11e is stronger to 2P-ERK2 than 0P-ERK2 (Rudolph et al., 2015). In contrast, GDC0994 and SCH772984 bind both active and inactive kinase forms with equal affinity (Pegram et al., 2019; Ward et al., 2019).
HDX-MS also reveals differences between VTX11e, GDC0994, and SCH772984 in their ability to induce perturbations at the P+1 segment (Pegram et al., 2019). P+1 is distal to the kinase active site, and together with the activation loop forms part of the interface for substrate and effector binding. This suggests that R↔L conformational shifts within the active site can propagate out to the activation loop and surrounding regions. NMR-CPMG measurements confirm this, by showing that residues in the activation loop participate in the same global exchange motion as residues surrounding the active site (Iverson et al., 2020). Therefore, solution measurements of conformation selective inhibitors reveal allosteric coupling between the active site and distal regions around the activation loop.
A consequence of this behavior is that conformation-selective inhibitors can alter the interactions of ERK2 with other proteins. For example, the binding of VTX11e partially inhibits the rate of dephosphorylation at the activation loop by the MAP kinase phosphatase, MKP3/DUSP6, while SCH772984 enhances this rate (Pegram et al., 2019). The features of VTX11e controlling its property of R-state selection are unknown. An attempt to address this compared crystal structures of complexes formed between 2P-ERK2 and VTX11e (100% R), GDC0994 (80% R), or SCHCPD336, a close analog of SCH772984 (0% R) (Pegram et al., 2019; Chaikuad et al., 2014; Blake et al., 2016). But the results were inconclusive. Binding of SCHCPD336 disrupted Gly loop and helix αC interactions leading to a disordered activation loop, suggesting that the L-state might communicate N-lobe conformational changes to the activation loop dynamics. In contrast, VTX11e and GDC0994 formed similar close contacts with active site residues with no major distortions in the ATP pocket. Given the limited set of compounds with small structural differences, it was difficult to identify features of binding that could explain why VTX11e is selective for the R-state but GDC0994 is not.
Here, we survey a panel of small molecules with varying chemistries, using HDX-MS and NMR to classify their properties of conformation selection. Of the 19 inhibitors examined, 14 behave in a similar manner to VTX11e, with conformational selection for the R-state. These include 13 novel inhibitors and ulixertinib/BVD523, a promising ERK inhibitor with blood-brain barrier permeability (Sullivan et al., 2018; Yu et al., 2022; Sigaud et al., 2023). Four inhibitors resemble GDC0994 in allowing R↔L exchange with populations comparable to apoenzyme, and one allows exchange but partially shifts the equilibrium toward the R-state. NMR chemical shifts reveal perturbations unique to R-state inhibitors located in the P+1 segment and helix αL16, and remote from the inner-sphere region of ligand contact. Co-crystal structures of R-state selective inhibitors complexed with ERK2 show small but systematic conformational shifts in the Gly loop, helix αC, and helix αL16, in which the Gly loop moves closer to the DFG segment, while helicesαC and αL16 shift outwards. This movement rotates the N- and C-lobes into a more closed configuration and moves the K52-E69 salt bridge closer to the conserved DFG motif. Such small but systematic conformational shifts suggest a working model for R-state selection, involving movements of the N-lobe that favor a more compact active site with allosteric coupling to the activation loop.
Results
Solution measurements reveal R-state conformation selection
NMR experiments were used to examine conformation selection by inhibitor binding. Previously, NMR-CPMG measurements of ERK2, selectively [methyl-13C,1H]-labeled on Ile, Leu, and Val (ILV) residues, showed that phosphorylation introduced global motions involving multiple residues in the N- and C-lobes surrounding the active site, as well as in the activation loop (Xiao et al., 2014; Iverson et al., 2020). These could be modeled by a global two-state exchange process (R↔L) with an exchange rate constant kex = 300 s–1. As shown previously (Xiao et al., 2014; Iverson et al., 2020), the R:L populations were 80:20 at 25 °C and 50:50 at 5 °C for the 2P-ERK2 apoenzyme, and 0:100 for the 0P-ERK2 apoenzyme at both 25°C and 5°C.
NMR spectroscopy was used to test the effect of inhibitors on conformational exchange in 2P-ERK2. 2D (13C,1H)-HMQC spectra were collected on complexes of ERK2 with BVD523 and VTX11e, which share chemical features including a central amido-linked pyrrole scaffold, and GDC0994, which has a distinct central pyridone scaffold (Figure 2A; Figure 2—figure supplement 1; Figure 2—figure supplement 2). Each inhibitor was added to ERK2 in saturating amounts (1:1.2 ERK2:inhibitor), as illustrated for VTX11e and GDC0994 (Figure 2—figure supplement 3; Figure 2—figure supplement 4). Figure 2B–D illustrates methyl probes with slow exchange (Δω>>kex), where the relative populations of R and L separated by their chemical shift differences can be observed in HMQC NMR spectra. Saturation binding of BVD523 shifted the R↔L equilibrium to 100% R at both 25°C and 5°C, matching the effects of VTX11e (Figure 2B–D). In contrast, GDC0994 had little effect on the R and L populations relative to apoenzyme (Figure 2B–D). The results show that BVD523, like VTX11e, exhibits properties of R-state conformational selection. This behavior might explain why BVD523 preferentially binds the activated kinase, with ~10 fold higher affinity for 2P-ERK2 than 0P-ERK2 (Germann et al., 2017).
HDX-MS experiments were conducted as a complementary approach to compare the effects of BVD523 to VTX11e and GDC0994 on conformational selection (Supplementary file 1). Binding of all three inhibitors reduced deuteration at conserved regions in the active site, including the Gly-rich loop, helix αC, the hinge, helix αE, and helix αL16 to a comparable degree between 0P- and 2P-ERK2 (Figure 3A–C, light green; Figure 3—figure supplement 1; Figure 3—figure supplement 2; Figure 3—figure supplement 3A–E and G). All three inhibitors also reduced HDX in the β7-β8 loop, although only in 2P-ERK2 (Figure 3A–C, dark green; Figure 3—figure supplement 1; Figure 3—figure supplement 3F). Such patterns reflect steric protection from solvent due to the occupancy of ligands in the nucleotide-binding site. In contrast, the effects of inhibitors differed in the DFG motif and adjacent β9 strand, where binding of BVD523 and VTX11e led to greater HDX protection compared to GDC0994 (Figure 3A–C, blue; Figure 3D; Figure 3—figure supplement 4). BVD523 and VTX11e also induced HDX protection greater than GDC0994 in regions farther from the active site, including the P+1 segment and helix αF (Figure 3A–C, blue; Figure 3E and F). These effects at P+1 and helix αF reveal allosteric coupling between the active site and distal regions surrounding the activation loop. Thus, solution measurements reported conformation selection for the R-state by BVD523 and VTX11e, but not GDC0994. The similarities in chemical structure between BVD523 and VTX11e suggested that features in their central scaffold and/or right-side regions might contribute to their shared properties of R-state selection.
Conformation selection among small molecule inhibitors of ERK1/2
To expand the structural diversity of studied compounds, we surveyed a panel of 17 small molecule inhibitors of ERK2 (Figure 4A) that had been produced in the course of chemical development of GDC0994 (Blake et al., 2016). All of these were ATP-competitive, and all inhibited ERK1/2 with subnanomolar Ki. The compounds were chosen for variations in their left-side region, which interacts with the kinase hinge and solvent channel; their central scaffold, which contacts hydrophobic residues from strands β7-β8 and the Gly loop; and their right-side region, which interacts with the Gly loop. The left-side varied between pyrazole and tetrahydropyran groups, the central scaffold alternated pyridone and triazolopyrazine ring structures, and the right-side accommodated many groups with different sizes and chemistries (Figure 4B). Residue contacts with each region are illustrated for VTX11e and GDC0994 in Figure 4C.
HDX-MS was used to examine the effect of each inhibitor on the conformational mobility of ERK2. Thirteen inhibitors (#2, #3, #5-12, #14, #16, #17) produced HDX patterns resembling BVD523 and VTX11e, with strong protection of the DFG motif in both 0P- and 2P-ERK2, and protection of the P+1 segment and helix αF in 2P-ERK2 (Figure 5A-C; Figure 5—figure supplement 1). In contrast, three inhibitors (#4, #13, #15) showed HDX patterns resembling GDC0994, indicating lesser protection of the DFG motif and no protection in the regions surrounding the activation loop (Figure 5A-C; Figure 5—figure supplement 2). One inhibitor (#1) showed strong protection of DFG, but weak protection of the activation loop regions compared to BVD523/VTX11e (Figure 5A-C). The effect of each inhibitor on deuterium uptake was quantified for these peptide segments by summing the differences in centroid mass between apo- and ligand-bound enzymes over all time points. This yielded a difference area measurement (dAUC), which was then ranked for each segment. Figure 6 shows rankings plotted to compare the DFG and P+1 peptides. The results showed good agreement between the behavior of each segment, with obvious separation in both dimensions between thirteen inhibitors clustering with BVD523 and VTX11e, three clustering with GDC0994, and one in between.
The HDX-MS results indicated that many compounds that are structurally distinct from BVD523 and VTX11e nevertheless share properties of R-state selection when bound to 2P-ERK2, while others resemble GDC0994 in maintaining conformational exchange. This was confirmed for a subset of seven inhibitors by HMQC NMR, chosen to represent variations in left-side, central scaffold, and right-side substituents. Methyl probes reporting R↔L exchange in 2P-ERK2 showed that four inhibitors that clustered with BVD523 and VTX11e by HDX also shifted populations to 100% R by NMR, while two inhibitors that clustered with GDC0994 maintained R:L population ratios of approximately 80:20 (Figure 7A–C; Figure 7—figure supplement 1; Figure 7—figure supplement 2). Inhibitor #1, which appeared between the two main clusters in Figure 6, shifted to the R-state, but only to a partial degree that still allowed exchange (Figure 7A–C). Thus, the HDX and NMR solution measurements are in good agreement and revealed R-state selection in 2P-ERK2 among a majority of the 17 compounds surveyed, even though they are chemically distinct from VTX11e and BVD523.
NMR chemical shifts report inhibitor chemistry and long-distance allostery
Next, we examined NMR chemical shifts to identify [methyl-13C,1H]-ILV probes responsive to binding the different inhibitors (Figure 8—figure supplement 1; Supplementary file 2), as distinct from the methyl probes above reporting R↔L exchange. Chemical shift perturbations from apoenzyme were examined for correlations with different inhibitor features (Figure 8). All inhibitor complexes led to peak disappearances attributed to broadening (indicated by red asterisks below the x-axes in Figure 8—figure supplement 1). These occurred at residues expected to be in direct contact with ligands, based on kinase-inhibitor structures (Roskoski, 2016). They mapped to strand β7 (L154) below the inhibitor binding site, and the αC-β4 loop (I82) which contacts the inhibitor from the back pocket (Figure 8, black; see Figure 4C). One residue forming Gly loop contacts (V37), as well as second-shell residues in strand β5 (I101) and the N-terminal end (V12) appeared broadened with all inhibitors except VTX11e and BVD523 (Figure 8).
Four methyl probes showed chemical shifts that could be correlated with the identify of the left-side substituent, which in ERK2:GDC0994 and ERK2:VTX11e are located near the hinge and interact with the solvent-exposed front pocket (Figure 4C; Pegram et al., 2019; Chaikuad et al., 2014; Blake et al., 2016). Figure 8 shows the locations of these residue probes in green. Here, chemical shift perturbations separated the inhibitors #4, #6, #8, and #15, which share tetrahydropyran at the left-side position, from #1, #5, #16, and GDC0994, which share pyrazole. Three methyl probes (L105, L155, L161) were located near residues that directly contact the inhibitor left-side, and one (L26) was located in strand β1, nearby a hydrophobic side chain (I29) that contacts the pyrazole ring in ERK2:GDC0994 (Pegram et al., 2019; Blake et al., 2016). Thus, several peaks displayed chemical shifts or broadening consistent with direct binding interactions.
In contrast, two probes showed chemical shift behaviors that correlated with R-state conformation selection (Figure 8, red). These residues were located away from the active site, within the P+1 segment (I196) and the C-terminal helix αL16 (I345). Here, R-state selective inhibitors (VTX11e, BVD523, #5, #6, #8, #16) induced greater chemical shift perturbations, while inhibitors that retained conformational exchange (GDC0994, #4, #15) were closer to apoenzyme. The location of I196 was consistent with the HDX patterns showing greater protection in the P+1 loop by each of the R-state selective inhibitors. The location of I345 in helix αL16 suggested an association of the R-state with N-lobe perturbations proximal to helix αC, which interacts with αL16 (Figure 1). Thus, NMR identified changes occurring remotely from the ligand binding site that were sensitive to shifts towards the R-state, as reported by HDX-MS. The results helped confirm that the binding of the new R-state inhibitors to the ATP pocket couples to changes in distal regions, including those around the activation loop.
R-state selection properties correlate with conformational shifts in the N-lobe
We asked if structural features accompanying R-state selection might be observable by X-ray crystallography. Structures of 2P-ERK2 complexed with inhibitors #8 and #16 were solved (Supplementary file 3) and compared to previous structures of 2P-ERK2 complexed with VTX11e and GDC0994 (Pegram et al., 2019). Both inhibitors #8 and #16 formed contacts with active site residues that matched those contacting VTX11e and GDC0994 (Figure 9A). They included hydrogen bonds between the left-side aminopyrimidine and the hinge backbone (M106) and tetrahydropyran with the front pocket (K112); direct or water-mediated hydrogen bonds between the central triazolopyridine scaffold and the gatekeeper (Q103), Lys-Glu salt bridge (K52), and DFG motif (D165 backbone amide); and van der Waals interactions between the right-side chlorophenyl or isopentane groups and the Gly loop.
While conformational differences between each of the inhibitor complexes were small, systematic shifts nevertheless appeared in the N-lobe, involving the Gly loop, helix αC, and helix αL16. In the co-crystal structure with VTX11e, these elements shifted outwards (away from the ligand and towards the activation loop) relative to GDC0994 (Figure 9B). This difference could be attributed to a residue contact unique to VTX11e, whose right-side 3-chlorobenzyl substituent forms a Cl-π contact with Y34 on the Gly loop (Figure 9B left panel, see Figure 4C). Y34 in turn forms π-π ring-stacking interactions with Y62 on helix αC, suggesting cooperative ternary anion-π-π interactions (Lucas et al., 2016). These contacts, together with the K52-E69 salt bridge, likely couple movements of the β1-β2-β3 sheet to helix αC, which in turn promotes movement of the adjoining helix αL16. This results in the repositioning of helix αC and αL16 to accommodate the bulky chlorobenzyl group on VTX11e and moves the Gly loop to a more closed conformation. In contrast, the corresponding 3-fluoro 4-chlorobenzyl substituent on GDC0994 is positioned higher in the binding site and tilted upwards into the Gly loop relative to VTX11e (Figure 9B, right panel). This holds the Gly loop in a more open conformation and shifts the position of helix αC and αL16 inwards (towards the ligand). Thus, the variations between right-side substituents and their interactions with the Y34-Y62 stack underlie small conformational shifts in the N-lobe that affect domain opening and closure.
Like VTX11e, complexes of 2P-ERK2 with inhibitors #8 and #16 (PDBID: 8U8K, 8U8J) showed outward shifts of the Y34-Y62 stack, the Gly loop, and helices αC and αL16, relative to GDC0994 (Figure 9C and D, left panel). Unlike VTX11e, the 2-chloropyridine and isopentane substituents on inhibitors #8 and #16 respectively were located more than 6 Å from Y34. Both contacted the backbone of strand β1 at angles that allowed but were unlikely to force Gly loop closure (Figure 9C and D, right panels). Inspection of their triazolopyridine central scaffold showed interactions that were distinct from the pyridone scaffold in GDC0994. In GDC0994, the pyridone carbonyl oxygen formed hydrogen bonds with both K52 and a bound water (W1) coordinated to the gatekeeper side chain (Q103) and the DFG main chain amide (D165) (Figure 9E). In contrast, the two nitrogens in the triazole ring of inhibitors #8 and #16 formed separate hydrogen bonds to K52 and W1, respectively. This separation had the effect of shifting the K52-E69 salt bridge position relative to GDC0994. As a result, the Gly loop shifts to a more closed state, allowing the pocket below the Gly loop to accommodate two water molecules that form OH-π interactions with Y34, and favoring outward movement of helicesαC and αL16 (Figure 9C and D). In this way, the R-state inhibitors VTX11e, #8, and #16 share small but systematic movements that lead to Gly loop closure and movements of helixαC and αL16 in the direction of the activation loop. Differences were also seen in the hydrogen bond network around the conserved DFG region. Here, a second bound water (W2) formed a bridge between the catalytic residues K52 and D165 in complexes with VTX11e, #8, and #10, but was absent in GDC0994 (Figure 9B and D). This revealed a relayed network of hydrogen bond interactions linking the DFG motif to the K52-E69 salt bridge that appeared common to the three R-state inhibitors. Conceivably, this might reflect new bonding interactions that favor greater compactness of the active site.
Taken together, the crystal structures showed only small conformational variations among these closely related ATP-competitive inhibitors. Nevertheless, certain features were common to each of the R-state inhibitors, involving Gly loop closure and outward movements of helix αC and αL16. These appeared to be favored either by differential contacts between the inhibitor and Y34, or different hydrogen bond configurations between the central scaffold and K52. In each case, coupling to helix αC was enabled by Y34-Y62 and K52-E69 residue interactions. In contrast, the bulky right-side substituent on GDC0994 that prevented Gly loop closure and movements of helix αC/αL16 correlated with retention of conformational exchange.
The results so far revealed correspondence in R-state selection between the amido-linked pyrrole scaffold of VTX11e and BVD523 and the triazolopyridine scaffold of the new inhibitors. We asked what might happen, if the central scaffold were chosen to more closely resemble VTX11e, while choosing a right-side substituent similar to GDC0994. Tizakertib/ATG017/AZD0364 is an ATP-competitive ERK inhibitor with a tetrahydropyrrolodiazepenone core that resembles the amidopyrrole group in VTX11e and BVD523, and a right-side 3,4-fluorobenzyl group that resembles the 3-fluoro 4-chlorobenzyl group in GDC0994 (Figure 10A; Flemington et al., 2021; Ward et al., 2019). Interestingly, the X-ray structure of ATG017 complexed with 0P-ERK2 (Figure 10B, green; PDBID:6SLG, Ward et al., 2019) showed an open Gly loop, relative to complexes of BVD523 and VTX11e with 0P-ERK2 (shown for 0P-ERK2:BVD523, Figure 10B, light blue; PDBID:6GDQ, Heightman et al., 2018). This was due to the rotation of the 3,4-fluorobenzyl ring in ATG017 away from Y34 to form closer interactions with backbone atoms in the β1 and β2 strands. Unlike inhibitors #8 and #16, bound waters between the 3,4-fluorobenzyl ring and Y34 were absent, due to the open Gly loop conformation. As a result, Y34 shifts inwards along with helices αC and αL16, and the side chain of Y64 appears disordered (Figure 10B). The inward direction of this conformational shift suggested that ATG017 might share exchange properties with GDC0994 and the other pyridone-based inhibitors, despite their differences in central scaffold chemistry.
To examine this possibility, ATG017 was complexed with 2P-ERK2 and examined by NMR and HDX-MS. The NMR HMQC spectra revealed a small population shift to the R state, but with overall retention of R↔L exchange properties, similar to apoenzyme (Figure 10C). Likewise, HDX-MS showed little protection around the P+1 and helix αF segments (Figure 10D), with a pattern that was clearly distinct from the larger HDX protection induced by VTX11e or BVD523 (Figure 3E and F). Thus, ATG017 showed properties of conformational exchange and inward shifts of helices αC and αL16 that were shared with GDC0994 and inhibitors #4 and #15. This could be attributed to its right-side difluorobenzyl substituent that held the Gly loop in an open state. Interestingly, a bound water (W2) bridging K52 and D165 was present in the ATG017 complex (Figure 10B), and stronger HDX protection around the DFG segment and β9 strand was observed (Figure 10D; Figure 10—figure supplement 1). This implies that the HDX protection pattern around the DFG motif is functionally linked to the water bridge (W2) linking K52 and D165. The results show that while changes around the DFG region may contribute to R-state selection, the movements of the Gly loop and helices αC and αL16 appear more essential for conformation selection and allosteric communication.
Discussion
ERK inhibitors are promising therapeutic targets for cancers harboring oncogenic BRAF or RAS with acquired resistance to RAF and MEK inhibitors (Ryan et al., 2015; Morris et al., 2013; Hatzivassiliou et al., 2012; Roskoski, 2019; Pan et al., 2022). A feature of ERK2 activation that is somewhat unique among protein kinases is its global exchange behavior, involving millisecond interconversion between at least two states in its active, dual phosphorylated form. In this study, we present new ATP-competitive inhibitors of ERK1/2 that, despite their chemical divergence from VTX11e and BVD523, nevertheless share properties of R-state selection. They display sensitivity to the conformational differences between the R and L-states and the ability to induce long-distance perturbations at the P+1 segment and helix αL16. Crystallographic structures show correlations between the property of R-state selection and small but systematic shifts in the N-lobe that move helices αC and αL16 in the direction of the activation loop. This appears to be facilitated by movements of strands β1-β2-β3 coupled to helix αC. The findings reveal the importance of contacts between strand β3 and helix αC via the positioning of the K52-E69 salt bridge in ERK2, as well as the unique π−π contacts between Tyr residues in the Gly loop and helix αC, in guiding conformation selection for the R-state. By contrast, the binding of SCH772984 disrupts β3-αC contacts in ERK2 (Chaikuad et al., 2014), which may explain why this inhibitor favors the L-state (Pegram et al., 2019).
Unlike many other kinases, activation of ERK leads to relatively small structural changes within the active site. In 0P-ERK2, the N and C lobes are held apart by the backbone atoms of G167 in the DFG motif, which interdigitates with R65 and Q64 in helix αC (Zhang et al., 1994). This domain separation has been proposed to account for the low activity of unphosphorylated ERK2. Following phosphorylation, residue R65 rotates along with helix αC in the direction of the activation loop and pT183, partially relieving the barrier formed by G167. The resulting domain rotation between N- and C-lobes is only 5 degrees, but it brings K52 closer to D165 by 1.5 Å. The small magnitude of these changes are in part due to structural constraints unique to ERK2. Here, residues in the C-terminal L16 segment extending from the C-lobe to the N-lobe (e.g. Y316, D319, E320, F327, and M331) form close contacts with residues in helices αC and αE and the activation loop (Figure 11A). They culminate in helix αL16, where residues I345 and F346 fill a hydrophobic pocket formed by residues in helix αC and the β4-β5 loop (Figure 11B). This placement of helix αL16 aligns structurally with hydrophobic motif (HM) sequences at the C-termini of PKA and members of the AGC kinase subfamily (Taylor et al., 2022; Kannan et al., 2007; Baffi and Newton, 2022), as well as cyclins complexed with CDK (Tatum and Endicott, 2020). Thus helix αL16 coincides with regions in other kinases also involved in long-distance allostery.
Many of the key [methyl-13C,1H]-ILV probes that report global exchange behavior in ERK2 are part of the extensive network of hydrophobic residues connecting these substructures. Thus, the side chain of I345 fills a deep hydrophobic pocket formed by T66, L67, and I70 in helix αC, F57 in the β3-αC loop, V99 in strand β5, and M96 in loop β4-β5 (Figure 11B). This network also includes contacts formed by αL16 residue F346 with R89, I93, and M96 in the β4-β5 loop, and F352 with L74, I84, and I87 in αC and β4. Furthermore, side chains from F327 and M331 in L16 fill a hydrophobic pocket formed by I72, R75, and F76 in helix αC, V143 in the αE-β6 loop, and V171 in the activation loop (Figure 11A). In turn, I72 and F76 form hydrophobic contacts with L73 in helix αC, I138 in helix αE, and F166 in the DFG motif. Thus, αL16 and L16 contact the N-lobe via a hydrophobic network that merges with the R-spine (L73, I84, F166). Such extensive packing interactions may help to constrain ERK2, preventing the backbone rotational movements needed to allow a DFG-out configuration (Vijayan et al., 2015), which is never observed in the wild-type kinase. Conceivably, the conformations involved in R↔L global exchange could reflect systematic movements within this network, which propagate to produce the chemical shift perturbations observed in I345. Although such conformational shifts are small, the coupling between movements of the Gly loop and helices αC and αL16 are complementary to the N-lobe movements seen upon ERK2 phosphorylation. It is remarkable that inhibitors that share similar binding interactions and active site contacts are sensitive to such small conformational differences.
This study shows how HDX-MS and NMR are sensitive to slight perturbations within the active site that lead to significant long-distance responses at the activation loop and P+1 segment. Large structural rearrangements around the activation loop were not obvious in the X-ray structures of 2P-ERK2 complexed with GDC0994, VTX11e, or inhibitors #8 and #16. However, evidence that these measurements truly report conformational changes around pT183 and pY185 was supported by the ability of VTX11e to inhibit their rate of dephosphorylation by MAP kinase phosphatase (MKP3/DUSP6, Pegram et al., 2019). This demonstrates the protection of the activation loop by VTX11e binding. Conceivably, the conformation of the activation loop in solution may be more variable than in crystal structures of 2P-ERK2, in part due to lattice restraints. This is suggested by the extensive lattice contacts with the activation loop and P+1 segment in complexes formed by VTX11e, GDC0994, and inhibitors #8 and #16 (Figure 11—figure supplement 1). In fact, crystal contacts are positioned in a manner that might influence the activation loop in most X-ray structures of ERK2 (Pegram et al., 2023). This highlights the need for solution measurements to more fully understand allosteric regulation in ERK2.
Recently, we conducted extended molecular dynamics simulations (225 μs) starting from the X-ray structure of 2P-ERK2 (PDBID: 2ERK) (Pegram et al., 2023). The results showed unexpected flexibility of the activation loop, which deviated from the X-ray fold in all trajectories and was able to adopt multiple long-lived conformational states with RMSF ≤1.2 Å and lasting for more than 5 μs. Each state appeared only after 1–5 μs, perhaps due to biases introduced by lattice contacts in the starting model. Importantly, the states varied significantly in the contacts formed between the activation loop, L16, and the P+1 segment. These variable activation loop conformations were associated with differential dynamics within the active site, supporting allosteric coupling between these regions. Interestingly, in one of the novel states, the Arg salt bridges formed between pY185 and the P+1 segment were broken, allowing the phosphorylated side chain to move away and form a new salt bridge with helix αC. This reorganization resembled that observed in a co-crystal structure of 2P-ERK2 and the adaptor protein, PEA15 (Mace et al., 2013). The results suggest that distinct conformational states of the activation loop appear thermally accessible and likely to exist in solution. We speculate that they may function to accommodate variable binding modes for substrates and effectors. To the extent this is true, the coupling between the active site and activation loop in ERK2 might allow conformation selective inhibitors to control the binding of substrates and regulatory effectors. Similarly, allosteric properties in PKA allow different ATP-competitive inhibitors to vary the cooperativity of substrate binding (Olivieri et al., 2022).
Indeed, remote conformational changes induced by conformation selective inhibitors have been shown to alter noncatalytic scaffolding functions of protein kinases via interactions with regulatory enzymes (Agudo-Ibáñez et al., 2023; Leroux and Biondi, 2020; Fang et al., 2020; Sonti et al., 2018; Joseph et al., 2022). In fact, an intriguing recent study demonstrates that ubiquitination and degradation of 2P-ERK2 can be enhanced by BVD523, suggesting that conformation selective inhibitors may act as ‘kinase degraders’ that preferentially target the activated kinase (Balmanno et al., 2023). Other noncatalytic functions of ERK2 include the activation of topoisomerase IIα and poly-ADP ribose polymerase via direct binding of the phosphorylated kinase (Shapiro et al., 1999; Cohen-Armon et al., 2007). Preferential binding of DNA to 2P-ERK2 has also been shown to function in ERK1/2 binding to chromatin DNA, with a proposed regulation of transcription through phosphorylation-independent mechanisms (Hu et al., 2009; Lawrence et al., 2008; McReynolds et al., 2016; Tee et al., 2014). Thus, the ability of R-state inhibitors to control solution conformations of the activation loop may influence emerging noncatalytic functions of ERK.
In summary, the case of ERK2 illustrates how very small conformational changes in protein kinases can lead to considerable effects on activity and inhibition. Correspondingly, small variations between inhibitors which otherwise share similar binding interactions and active site contacts, differentially control long-distance allosteric effects that induce conformational changes at the activation loop. Such differences may be important for cellular responses, by altering the conformation of ERK at a key interface for substate and effector recognition. This highlights the importance of understanding the properties of conformation selection in ERK, and how these unique features of ERK regulation impact inhibitor design. The elements of protein architecture that enable small active site movements to yield large effects in distal regions remain an important and exciting feature of MAP kinases to explore in future investigations.
Materials and methods
Reagents and purified proteins
Request a detailed protocolVTX11e, GDC0994, and ATG017 were obtained from Selleck Chemicals (Cat# S7709, S7554, S8708; Houston, TX). Dimethylsulfoxide (DMSO), dithiothreitol (DTT), formic acid, and D2O (99.9%, Cat#756882) were obtained from Sigma-Aldrich (Burlington, MA). Ni-NTA agarose was obtained from Qiagen (Cat#30210, Germantown, MD). Inhibitors #1-#17 were synthesized in-house at Genentech, Inc, following methods described in WIPO Patents WO2015103137 A1 20150709 and WO2013130976 A1 20130906.
NMR spectroscopy
Request a detailed protocol[methyl-13C,1H]-Ile, Leu, and Val-labeled rat His6-ERK2 was prepared as described (Xiao et al., 2014; Xiao et al., 2015) in buffer containing 50 mM Tris-d11 pH 7.4 (D2O), 150 mM NaCl, 5 mM MgSO4, 0.1 mM EDTA, 1 mM D,L-1,4-DTT-d10, 100%(v/v) D2O, and 2.5% (v/v) glycerol-d8. Stock solutions (25 mM) of all inhibitors (VTX11e, BVD523, GDC0994, ATG017, Inhibitors #1, #4, #5, #6, #8, #15, #16) were prepared in DMSO-d6, and added to ERK2 (150 μM) to form complexes with 100% binding stoichiometry ([ERK2]:[inhibitor]=1.0:1.2). The DMSO-d6 concentration in the final protein sample was ~1% (v/v).
Two-dimensional (2D) (13C,1H) heteronuclear multiple quantum coherence (HMQC) spectra of 0P- and 2P-ERK2 were collected on a Varian 900 MHz NMR spectrometer. Data were collected at 25°C and 5°C each for a total time of 10 hr. Each spectrum was acquired with nonuniform sampling (NUS) with 160 complex points (900 MHz) in the t1 (13C) dimension, corresponding to 29.1 ms at 900 MHz, and 1024 complex points in the acquisition period. WURST40 13C decoupling was applied during the 85 ms acquisition period, and a 1.5 s delay period was used between each scan. The spectral processing was performed with the software package NMRPipe (Delaglio et al., 1995). Time domain data in the 1H dimension were apodized by a cosine window function and zero-filled prior to Fourier transformation. The indirect dimension (13C) was apodized by a cosine window function and zero-filled prior to Fourier transformation. NUS spectral reconstruction utilized the SMILE plugin (Ying et al., 2017). Spectral visualization and analysis were achieved using CCPNMR Analysis software (Vranken et al., 2005). NMR spectra of 0P-ERK2 and 2P-ERK2 apoenzymes were collected in two or more experiments on separate enzyme preparations with nearly identical results. Spectra of ERK2 complexed with inhibitors were each collected in one experiment.
X-ray crystallography
Request a detailed protocolX-ray structure images were created with Pymol 2.3.4 software (Schrödinger LLC, New York, NY). Co-crystals of human His6-2P-ERK2 in complex with Inhibitors #8 and #16 were grown at 20 °C in vapor diffusion hanging drops, mixing 0.5 µL protein +0.5 µL well solution (32–36% (w/v) PEG3350, 0.1 M Tris pH 8.8–9.0). Crystals were cryoprotected in well solution +10%(v/v) glycerol prior to data collection. X-ray diffraction data were collected on a Rigaku FR-E generator (Rigaku Americas Corporation, The Woodlands, TX) equipped with a Pilatus 2 M detector (Dectris USA, Philadelphia, PA). Data were processed using the Mosflm and Aimless programs in the CCP4 suite (Winn et al., 2011). Structures of ERK2 complexed with inhibitors were each collected in one experiment.
Kinase inhibition assays
Request a detailed protocolKinetic parameters (Ki) were determined from aggregate kinase assays of human 2P-ERK2 performed with the fluorescent Omnia peptide substrate, S/T17 (Cat. KNZ1171C, Invitrogen, Carlsbad, CA). Omnia assays contained 50 mM HEPES, pH 7.3, 10 mM MgCl2, 1 mM dithiothreitol, 0.005% Triton-X100, 5 nM ERK2, 6.25μΜ S/T17 peptide substrate, and 25 μΜ ATP (corresponding to the observed Km) in a total reaction volume of 25 μL. Assays were run at ambient temperature in 384-well plates, collecting time points every 50 s for 30 min on an Envision plate reader (PerkinElmer, Waltham, MA) with excitation 340 nm and emission 495 nm. Phosphorylation rates were normalized to controls (no inhibitor) and plotted against inhibitor concentration to obtain IC50 values using a four-parameter fit.
HDX-MS
Request a detailed protocolWild-type rat 0P-ERK2 and 2P-ERK2 were prepared as described, after bacterial expression from a plasmid constructed by Marcelo Sousa and Sandra Metzner, U. Colorado, Boulder (Pegram et al., 2019). ERK2 protein (17 μM) was incubated with 20 μM inhibitor (VTX11e, BVD523, GDC0994, ATG017, Inhibitors #1-#17) for 30 min in HDX buffer (50 mM potassium phosphate pH 7.2, 100 mM NaCl, 5 mM DTT, 10 mM MgCl2 and 0.5% (v/v) DMSO from concentrated ligand solutions). Deuterium uptake was initiated by adding 54 μL D2O (containing an equal concentration of ligand, DMSO, and MgCl2) to 6 μL ERK2 (complexed with or without ligand), to reach a final buffer concentration of 5 mM potassium phosphate pH 7.2, 10 mM NaCl, 0.5 mM DTT, 10 mM MgCl2, and 0.5% DMSO. Reactions were incubated at 25.0 °C for varying times between 30 s and 180 min (0.5, 1, 3, 8, 16, 30, 60, 90, 180 min). Exchange was quenched with 48 μL 100 mM potassium phosphate pH 2.2. Additionally, an in-exchange control was performed by the addition of quenching solution prior to initiating deuterium uptake, with no incubation period. All reactions were immediately injected in 100 μL volumes onto a Waters nanoAcquity HDX Manager UPLC system for proteolysis on an immobilized pepsin column (Enzymate, Waters Corp, Milford, MA) at 10 °C. Proteolysis, as well as peptide desalting (Vanguard Acquity UPLC BEH C18), were carried out in isocratic Solvent A (0.1%(v/v) formic acid in H2O) with flowrate 100 μL/min. Peptide separations (Waters Acquity UPLC BEH C18, 1.7 μm, 1.0 × 100 mm) were carried out at 3 °C using a 12 min linear gradient from 8–85% Solvent B (0.1% formic acid in acetonitrile) at 40 μL/min.
Peptides were analyzed by ESI-MS/MS using a Waters Synapt G2 HDMS Q-TOF mass spectrometer in positive ion mode. Continuum data were collected over 50–2000 m/z with 0.23 s scan time. Lock mass correction was achieved using Glu-Fibrinogen peptide (MH2+2 = 785.8426 m/z), with resolution = 20,000 at m/z 956. Undeuterated pepsin-cleaved peptides from 0P- or 2P-ERK2 were identified by MSe sequencing using PLGS 3.0 (Waters). Searches of rat ERK2 allowed nonspecific digestion and variable modifications of oxidized Met and phosphorylated Ser/Thr/Tyr, with parameter settings of intensity threshold = 750 counts, lock mass window = 0.4 Da, and low energy threshold = 135 counts. Peptides were accepted after filtering for high confidence, which required a maximum MH+ error of 15 ppm, minimum sequence length 3, maximum sequence length 25, minimum product ions 3, and identification in at least 2 of 4 replicate runs. Deuterium uptake was calculated for high-confidence peptides in each hydrogen exchange dataset using DynamX 3.0 (Waters). Raw files were processed using a 0.35 Da lock mass window and a low energy threshold of 100. For every experiment, automated ion assignments for deuterated peptides were manually inspected and validated. The weighted average mass of isotope distributions was calculated for each verified peptide using DynamX, and referenced to the weighted average mass of its unlabeled form. The recorded uptake from the in-exchange control was subtracted from the uptake values of all timepoints (0.5, 1, 3, 8, 16, 30, 60, 90, 180 min) to yield corrected uptake measurements that were used in plots of time courses and dAUC calculations. Supplementary file 1 contains deuterium uptake values without this correction. HDX time courses of 0P-ERK2 and 2P-ERK2 apoenzymes were collected in two or more experiments on separate enzyme preparations with nearly identical results. Time courses of ERK2 complexed with inhibitors were each collected in one experiment.
Dataset availability
Request a detailed protocolCrystallography data files have been deposited in the RSCB PDB under accession codes 8U8J and 8U8K. NMR datafiles have been deposited in BMRbig under accession code BMRbig96. Data files generated from HDX experiments have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository (Perez-Riverol et al., 2022), with the dataset identifier PXD048311.
Materials availability
Request a detailed protocolProprietary ERK inhibitors are available upon request from Dr. John Moffat, Genentech, Inc Recipients will be asked to sign a materials transfer agreement and will be restricted from distributing these compounds or use them in any applications related to for-profit research. Plasmids for the expression of ERK2 and constitutively active MKK1-G7B (Xiao et al., 2014; Xiao et al., 2015; Pegram et al., 2019) are available upon request from Dr. Natalie Ahn, University of Colorado - Boulder.
Data availability
Crystallography data have been deposited in the PDB under accession codes 8U8J and 8U8K. Data generated from HDX experiments, including Supplementary Dataset S1, coverage maps, and representative raw spectra have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository Perez-Riverol et al., 2022, with the dataset identifier PXD048311. Data generated from NMR experiments have been deposited in BMRbig under accession code BMRbig96.
-
PRIDEID PXD048311. HDX-MS Analysis of Conformation Selection by ATP-competitive Inhibitors and Allosteric Communication in ERK2.
-
BMRbigID BMRbig96. NMR Analysis of Chemical Shift Perturbations by ATP-competitive Inhibitors complexed with ERK2.
-
RCSB Protein Data BankID 8U8J. Co-crystal structure of phosphorylated ERK2 in complex with ERK1/2 inhibitor #16.
-
RCSB Protein Data BankID 8U8K. Co-crystal structure of phosphorylated ERK2 in complex with ERK1/2 inhibitor #8.
References
-
Dynamic regulatory features of the protein tyrosine kinasesBiochemical Society Transactions 47:1101–1116.https://doi.org/10.1042/BST20180590
-
mTOR regulation of AGC kinases: new twist to an old tailMolecular Pharmacology 101:213–218.https://doi.org/10.1124/molpharm.121.000310
-
Discovery of (S)-1-(1-(4-chloro-3-fluorophenyl)-2-hydroxyethyl)-4-(2-((1-methyl-1H-pyrazol-5-yl)amino)pyrimidin-4-yl)pyridin-2(1H)-one (GDC-0994), an extracellular signal-regulated kinase 1/2 (ERK1/2) inhibitor in early clinical developmentJournal of Medicinal Chemistry 59:5650–5660.https://doi.org/10.1021/acs.jmedchem.6b00389
-
A unique inhibitor binding site in ERK1/2 is associated with slow binding kineticsNature Chemical Biology 10:853–860.https://doi.org/10.1038/nchembio.1629
-
NMRPipe: a multidimensional spectral processing system based on UNIX pipesJournal of Biomolecular NMR 6:277–293.https://doi.org/10.1007/BF00197809
-
Defining and targeting BRAF mutations in solid tumorsCurrent Treatment Options in Oncology 22:30.https://doi.org/10.1007/s11864-021-00827-2
-
ERK inhibition overcomes acquired resistance to MEK inhibitorsMolecular Cancer Therapeutics 11:1143–1154.https://doi.org/10.1158/1535-7163.MCT-11-1010
-
The conformational state of the BTK substrate PLCγ contributes to ibrutinib resistanceJournal of Molecular Biology 434:167422.https://doi.org/10.1016/j.jmb.2021.167422
-
ERK signalling: a master regulator of cell behaviour, life and fateNature Reviews. Molecular Cell Biology 21:607–632.https://doi.org/10.1038/s41580-020-0255-7
-
Structure-guided strategy for the development of potent bivalent ERK inhibitorsACS Medicinal Chemistry Letters 8:726–731.https://doi.org/10.1021/acsmedchemlett.7b00127
-
Renaissance of allostery to disrupt protein kinase interactionsTrends in Biochemical Sciences 45:27–41.https://doi.org/10.1016/j.tibs.2019.09.007
-
Development of small molecule extracellular signal-regulated kinases (ERKs) inhibitors for cancer therapyActa Pharmaceutica Sinica B 12:2171–2192.https://doi.org/10.1016/j.apsb.2021.12.022
-
Activation loop plasticity and active site coupling in the MAP kinase, ERK2Journal of Molecular Biology 435:168309.https://doi.org/10.1016/j.jmb.2023.168309
-
The PRIDE database resources in 2022: A hub for mass spectrometry-based proteomics evidencesNucleic Acids Research 50:D543–D552.https://doi.org/10.1093/nar/gkab1038
-
The MAPK cascades: signaling components, nuclear roles and mechanisms of nuclear translocationBiochimica et Biophysica Acta 1813:1619–1633.https://doi.org/10.1016/j.bbamcr.2010.12.012
-
ERK1/2 MAP kinases: structure, function, and regulationPharmacological Research 66:105–143.https://doi.org/10.1016/j.phrs.2012.04.005
-
Targeting ERK1/2 protein-serine/threonine kinases in human cancersPharmacological Research 142:151–168.https://doi.org/10.1016/j.phrs.2019.01.039
-
Targeting RAS -mutant cancers: Is ERK the key?Trends in Cancer 1:183–198.https://doi.org/10.1016/j.trecan.2015.10.001
-
Vertical pathway inhibition overcomes adaptive feedback resistance to KRASG12C inhibitionClinical Cancer Research 26:1633–1643.https://doi.org/10.1158/1078-0432.CCR-19-3523
-
Extracellular signal-regulated kinase activates topoisomerase IIα through a mechanism independent of phosphorylationMolecular and Cellular Biology 19:3551–3560.https://doi.org/10.1128/MCB.19.5.3551
-
ATP site ligands determine the assembly state of the abelson kinase regulatory core via the activation loop conformationJournal of the American Chemical Society 140:1863–1869.https://doi.org/10.1021/jacs.7b12430
-
Resistance mechanisms to targeted therapy in BRAF-mutant melanoma - A mini reviewBiochimica et Biophysica Acta. General Subjects 1865:129736.https://doi.org/10.1016/j.bbagen.2020.129736
-
Chatterboxes: the structural and functional diversity of cyclinsSeminars in Cell & Developmental Biology 107:4–20.https://doi.org/10.1016/j.semcdb.2020.04.021
-
The tails of protein kinase AMolecular Pharmacology 101:219–225.https://doi.org/10.1124/molpharm.121.000315
-
Conformational analysis of the DFG-out kinase motif and biochemical profiling of structurally validated type II inhibitorsJournal of Medicinal Chemistry 58:466–479.https://doi.org/10.1021/jm501603h
-
Overview of the CCP4 suite and current developmentsActa Crystallographica D 67:2352.https://doi.org/10.1107/S0907444910045749
-
Structure of extracellular signal-regulated kinase 2 in complex with ATP and ADPActa Crystallographica. Section F, Structural Biology and Crystallization Communications 68:1434–1439.https://doi.org/10.1107/S1744309112042972
Article and author information
Author details
Funding
National Institutes of Health (R35GM136392)
- Natalie G Ahn
National Institutes of Health (R01GM114594)
- Natalie G Ahn
National Institutes of Health (S10RR026641)
- Natalie G Ahn
National Institutes of Health (S10OD025267)
- Natalie G Ahn
National Institutes of Health (S10OD025020)
- David N Jones
National Institutes of Health (P30CA046934)
- David N Jones
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We are indebted to Drs. Marcelo Sousa, Elan Eisenmesser, and Johannes Rudolph for guidance on structural analyses and many valuable discussions, and Josh Anderson for insightful discussions about data analysis and presentation. We also thank Drs. Thomas Lee, Chris Ebmeier, and Rachel Mehaffey for mass spectrometric analyses of purified proteins, Chang Liu for preliminary HDX studies, and Daniel Lee for Python script coding for HDX datasets. Special thanks to the Array Structural Biology group for their valuable contributions to X-ray structure determination, and to Dr. Yao Xiao for NMR inhibitor binding titrations in this study. This work was supported by NIH Research Awards R35GM136392 and R01GM114594 (NGA), NIH Instrumentation Awards S10RR026641 and S10OD025267 (mass spectrometry) and S10OD025020 (NMR), and Cancer Center Support Grant P30CA046934 (CU Anschutz NMR Shared Resource Facility). Work on the development and synthesis of ERK inhibitors and determination of X-ray structures was funded by Genentech, Inc.
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.91507. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2023, Anderson et al.
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.
Metrics
-
- 913
- views
-
- 129
- downloads
-
- 7
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
- Cell Biology
Our latest Special Issue brings together research on the allosteric regulation of kinase activity.
-
- Biochemistry and Chemical Biology
The conformational ensemble and function of intrinsically disordered proteins (IDPs) are sensitive to their solution environment. The inherent malleability of disordered proteins, combined with the exposure of their residues, accounts for this sensitivity. One context in which IDPs play important roles that are concomitant with massive changes to the intracellular environment is during desiccation (extreme drying). The ability of organisms to survive desiccation has long been linked to the accumulation of high levels of cosolutes such as trehalose or sucrose as well as the enrichment of IDPs, such as late embryogenesis abundant (LEA) proteins or cytoplasmic abundant heat-soluble (CAHS) proteins. Despite knowing that IDPs play important roles and are co-enriched alongside endogenous, species-specific cosolutes during desiccation, little is known mechanistically about how IDP-cosolute interactions influence desiccation tolerance. Here, we test the notion that the protective function of desiccation-related IDPs is enhanced through conformational changes induced by endogenous cosolutes. We find that desiccation-related IDPs derived from four different organisms spanning two LEA protein families and the CAHS protein family synergize best with endogenous cosolutes during drying to promote desiccation protection. Yet the structural parameters of protective IDPs do not correlate with synergy for either CAHS or LEA proteins. We further demonstrate that for CAHS, but not LEA proteins, synergy is related to self-assembly and the formation of a gel. Our results suggest that functional synergy between IDPs and endogenous cosolutes is a convergent desiccation protection strategy seen among different IDP families and organisms, yet the mechanisms underlying this synergy differ between IDP families.