Abstract
Allosteric cooperativity between ATP and substrates is a prominent characteristic of the cAMP-dependent catalytic subunit of protein kinase A (PKA). Not only this long-range synergistic action is involved in substrate recognition and fidelity, but it is also likely to regulate PKA association with regulatory subunits and other binding partners. To date, a complete understanding of the molecular determinants for this intramolecular mechanism is still lacking.
Here, we integrated NMR-restrained molecular dynamics simulations and a Markov State Model to characterize the free energy landscape and conformational transitions of the catalytic subunit of protein kinase A (PKA-C). We found that the apoenzyme populates a broad free energy basin featuring a conformational ensemble of the active state of PKA-C (ground state) and other basins with lower populations (excited states). The first excited state corresponds to a previously characterized inactive state of PKA-C with the αC helix swinging outward. The second excited state displays a disrupted hydrophobic packing around the regulatory (R) spine, with a flipped configuration of the F100 and F102 residues at the αC-β4 loop. To experimentally validate the second excited state, we mutated F100 into alanine (F100A) and used NMR spectroscopy to characterize the structural response of the kinase to ATP and substrate binding. While the catalytic efficiency of PKA-CF100A with a canonical peptide substrate remains unaltered, this mutation rearranges the αC-β4 loop conformation, interrupting the structural coupling of the two lobes and abolishing the allosteric binding cooperativity of the enzyme. The highly conserved αC-β4 loop emerges as a pivotal element able to control the synergistic binding between nucleotide and substrate. These results may explain how mutations or insertions near or within this motif affect the function and drug sensitivity in other homologous kinases.
eLife assessment
This important study provides an example of integrating computational and experimental approaches that lead to new insights into the energy landscape of a model kinase. Compelling use of molecular dynamics simulations and NMR spectroscopy provide a conformational description of active and excited states of the kinase; one of which has not been captured in previously solved crystal structures. Overall, this comprehensive study expands our understanding of the architecture and allosteric features of the conserved bilobal kinase domain structure.
Significance of findings
important: Findings that have theoretical or practical implications beyond a single subfield
- landmark
- fundamental
- important
- valuable
- useful
Strength of evidence
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
Introduction
Eukaryotic protein kinases (EPKs) are plastic enzymes of paramount importance in signaling processes, catalyzing phosphoryl transfer reactions, or acting as scaffolds for other enzymes and/or binding partners 1,2. Of all kinases, the catalytic subunit of PKA (PKA-C) was the first to be structurally characterized by X-ray crystallography.3,4 In its inhibited state, PKA-C assembles into a heterotetrametric holoenzyme comprising two catalytic (C) and two regulatory (R) subunits.5 The canonical activation mechanism of PKA involves the binding of two cAMP molecules that disassemble the holoenzyme, unleashing active PKA-C monomers, which target signaling partners.6 In 2017, however, Scott and coworkers proposed an alternative activation mechanism in which the holoenzyme does not disassemble under physiological conditions; rather, it is clustered in signaling islands, localized near the substrates by A-kinase anchoring proteins (AKAPs), facilitating targeted signaling interactions.7 To date, the activation mechanism of PKA is still under active investigation since PKA-C appears to be fully liberated under pathological conditions.8–10
X-ray crystallography studies revealed that PKA-C is a bilobal enzyme, with a dynamic N-terminal lobe (N-lobe) comprising four β-sheets and an αC helix, while the C-terminal lobe (C-lobe) is more rigid and mostly composed of α-helices (Figure 1A).3,4,11 The N-lobe harbors the nucleotide-binding site, whereas the substrate binding cleft lays at the interface between the N- and C-lobe. The three-dimensional structure of PKA-C features a highly conserved hydrophobic core decorated with catalytically important motifs, i.e., the Gly-rich loop, DFG-loop, activation loop, and magnesium and peptide positioning loops.12 In the catalytically active state, these motifs are all poised for phosphoryl transfer,13–15 a condition necessary but not sufficient to define an active kinase. More recent studies revealed a critical role of the enzyme’s hydrophobic core, which is crossed by a catalytic (C) spine and a regulatory (R) spine surrounded by shell residues (Figure 1B).16,17 The C spine comprises an array of hydrophobic residues and is assembled upon binding ATP, whereas the R spine is engaged when the activation loop is phosphorylated at Thr197, which positions the αC-helix in an active configuration.18
A distinct property of PKA-C is the binding cooperativity between ATP and substrate.19 In the catalytic cycle, the kinase binds ATP and unphosphorylated substrates with positive binding cooperativity, whereas a negative binding cooperativity between ADP and phosphorylated substrate characterizes the exit complex.20 The biological importance of these phenomena has been emphasized by our recent studies on disease-driven mutations,21–23 which all feature disrupted cooperativity between nucleotides and protein kinase inhibitor (PKI) or canonical substrates.21–23 Since the recognition sequence of substrates is highly homologous to that of the R subunits,5 a loss of cooperativity may affect not only substrate binding fidelity but also regulatory processes. Thus far, the molecular determinants for the binding cooperativity between nucleotide and sub-strate and its role in PKA signalosome remain elusive.
Aiming to define the allosteric mechanism for PKA-C cooperativity, we combined NMR-restrained replica-averaged metadynamics (RAM)24,25 and Markov State Model (MSM)26 and charted the conformational landscape and dynamic transitions of the main forms of PKA-C. We found that the apo kinase occupies three distinct basins: (i) a most populated ground state (GS) with constitutively active conformations competent for catalysis, (ii) a first high free energy basin representative of typical inactive states with a dislodged configuration of the αC helix (ES1), and (iii) a second high free energy basin with a disrupted hydrophobic array of residues at the core of the enzyme (ES2). Notably, the equilibrium between the most populated ground state and the other low-populated states agrees with previous NMR Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion and chemical exchange saturation transfer (CEST) measurements.27 To compound the existence of the second inactive basin, we mutated F100 at the αC-β4 loop into Ala (PKA-CF100A) and characterized its ligand binding thermodynamics and structural response by isothermal titration calorimetry (ITC) and NMR spectroscopy, respectively. We found that PKA-CF100A phosphorylates a canonical peptide with a catalytic efficiency similar to the wild-type kinase (PKA-CWT); however, this mutation abolishes the binding cooperativity between nucleotide and substrate. The NMR experiments reveal that F100A perturbs the hydrophobic packing around the αC-β4 loop and interrupts the allosteric communication between the two lobes of the enzyme. Overall, these results emphasize the pivotal role of the αC-β4 loop in kinase function and may explain why single-site mutations or insertion mutations stabilizing this motif in homologous kinases result in oncogenes and confer differential drug sensitivity.28–30
Results
The free energy landscape of PKA-C charted by NMR-restrained replica-averaged metadynamics (RAM) simulations
To characterize the experimentally accessible conformational landscape of PKA-C in the µs-to-ms timescale, we performed NMR-restrained metadynamics simulations within the RAM framework.24,25 Using four replicas restrained with backbone chemical shifts (CS), we simulated the apo, binary (PKA-C/ATP), and ternary (PKA-C/ATP/PKI5-24) complexes of PKA-C. Incorporating the CS restrains in the simulations ensures a close agreement with the conformational space explored by the kinase under our NMR experimental conditions,31,32 whereas the enhanced sampling through metadynamics simulations boosts the conformational plasticity of the enzyme along different degrees of freedom (i.e., collective variables, CVs. See Methods) (Figure 2 - figure supplement 1).33 In our case, back-calculations of the CS values with Sparta+ 34 show that the restrained simulations improved the agreement between theoretical and experimental CS values of ∼0.2 ppm for the amide N atoms and ∼0.1 ppm for the remainder backbone atoms (Figure 2 - figure supplement 2). Also, the bias-exchange metadynamics allow each replica to span a significantly broader conformational space relative to classical MD simulations (Figure 2 - figure supplement 3). The metadynamics simulations were equilibrated and converged in 300 ns (Figure 2A, Figure 2 - figure supplement 4). The full free energy landscape was then reconstructed by sampling an extra 100 ns of production phase with reduced bias along each CV. We found that the population of the PKA-C conformers is modulated by ATP and pseudosubstrate binding within the NMR detection limit of sparsely populated states (∼ 0.5% or ΔG < 3.2 kcal/mol).
Specifically, the apo PKA-C populates preferentially a ground state and five readily accessible low-populated excited states (Figure 2B, Figure 2 – supplementary table 1). The binary complex populates the same ground state and a broad higher energy basin (Figure 2C, Figure 2 – supplementary table 1). Finally, the ternary complex occupies a narrower dominant ground state (Figure 2D, Figure 2 – supplementary table 1). Overall, these calculations revealed that the apo kinase explores multiple minima along each CV. The binding of ATP shifts the population of the ensemble into an additional minimum, featuring structures that are committed to substrate binding. Finally, the conformational heterogeneity of the ensemble is significantly reduced upon formation of the ternary complex, which represents a catalytically committed state13 (Figure 2A). The free energy landscape obtained from these RAM simulations is consistent with the qualitative picture previously inferred from our NMR spin relaxation experiments and hydrogen/deuterium exchange studies,35–37 while providing a detailed structural characterization of the ground and sparsely populated conformationally excited states.
MSM reveals the conformational transitions of PKA-C from ground to high free energy (excited) states
To refine the free energy landscape and delineate the kinetics of conformational transitions of the kinase, we performed additional unbiased sampling to build Markov State Models (MSM).38,39 MSMs are commonly used to describe the dynamic transitions of various metastable states for macromolecules in terms of their populations and transition rates. MSMs are typically created by combining thousands of short unbiased simulations.38,39 Following this strategy, we performed several short simulations (10 – 20 ns) using thousands of the low free energy conformations (ΔG < 3.2 kcal/mol) chosen randomly from the three forms of PKA-C as starting structures. The conformational ensembles were clustered into microstates and seeded to start a second round of adaptive sampling (see Methods). This iterative process was repeated three times to assure convergence and yielded 100 μs trajectories for both the apo and binary forms, whereas 60 μs trajectories were collected for the less dynamic ternary complex. Once we reached a sufficient sampling, we built an MSM including L95, V104, L106, M118, M120, Y164, and F185 to investigate the dynamic transitions of the hydrophobic R spine and shell residues (Figure 3 – figure supplement 1). These residues are ideal reporters of the dynamic processes governing the activation/deactivation of the kinase.40 To compare the free energy landscape of different complexes, we projected the conformational ensembles of three forms and existing crystal structures of PKA-C onto the first two time-lagged independent components (tICs) of the apo form, which were obtained by a time-lagged independent component analysis (tICA). These tICs represent the directions of the slowest motion of the kinase and visualize the conformational transitions of V104, L95, and F185 (Figure 3 – figure supplement 1). The MSM shows that the hydrophobic core of the apo PKA-C accesses three major basins. The broadest basin represents the ground state (GS) and corresponds to the conformations of PKA-C captured in several crystal structures (Figure 3A). The other two states (conformationally excited states) are less populated. The first excited state (ES1) features a disrupted hydrophobic packing of L95, V104, and F185. This conformation matches an inactive state with an outward orientation of the αC-helix typical of the inhibited states of PKA-C (Figure 3 – Figure supplement 2). The second state (ES2) displays a flipped configuration of the V104 side chain and disrupted hydrophobic interactions responsible for anchoring the αC-β4 loop to the C-lobe, with F100 and F102 adopting a gauche+ configuration and forming a stable π-π stacking. (Figure 3A). In contrast, the active GS ensemble features a trans configuration for the F100 and F102 side chains that stabilize the π-π stacking with Y156 at the αE-helix and cation-π interactions with R308 at the C-tail (Figure 3A). To our knowledge, the ES2 state was never observed in the available crystal structures. Upon ATP binding, the confor-mational space spanned by the kinase becomes narrower, and the conformers populate mostly the GS, with a small fraction in the ES1 state (Figure 3B). This is consistent with the role of the nucleotide as an allosteric effector that enhances the affinity of the enzyme for the substrate.41 In the ternary form (ATP and PKI-bound), PKA-C populates only the GS consistent with the competent conformation observed in the ternary structure (i.e., 1ATP). In this case, the αC-β4 loop of the enzyme is locked in a well-defined, active configuration (Figure 3A).
Moreover, in the apo form, the αC-β4 loop is quite dynamic due to transient hydrophobic interactions involving F100 and V104, as well as W222 and the APE motif (A206 and P207) (Figure 4A). The binding of both nucleotide and PKI rigidifies the residues near F100 and V104, e.g., V103 and F185 (Figure 4A) and makes more persistent several electrostatic interactions essential for catalysis (D166-N171, K168-T201, and Y204-E230), which are transient in the apo PKA-C. Finally, we used kinetic Monte Carlo sampling to characterize the slow transition between different states.42 In the GS state of the apo PKA-C, F100 and F102 adopt trans configurations that stabilize the interactions with the αE and αJ helices, and, together with the nucleotide, they lock the αC-β4 loop in an active state (Figure 3A). The GS to ES1 transition features the disruption of the K72-E91 salt bridge, typical of the inactive kinase. In contrast, the GS to ES2 transition involves a 120° flip of the F100 aromatic group that interacts with V104, a conformation found only in the catalytically uncommitted apoenzyme (Figure 3A). Interestingly, the GS to ES2 transition involves a concerted disruption of the D166-N171 and K168-T201 electrostatic interactions and the destabilization of the packing between W222 and the APE motif (A206 and P207) required for substrate recognition (Figure 3C). Overall, these calculations show that conformational transitions from the active GS to ES1 and ES2 represent two independent pathways toward inactive states of the kinase.
Direct correspondence between the conformationally excited states identified by MD simulations and NMR data
NMR CPMG relaxation dispersion and CEST experiments performed on the apo PKA-C revealed the presence of conformationally excited states for several residues embedded into the hydrophobic core of the enzyme.27 The chemical shifts of one of these states follow the previously observed open-to-closed transitions of the enzyme,41 whereas another excited state observable in the CEST experiments does not follow the same trend. We surmised that this new state may represent an alternate inactive state.27 Since ES2 features a disrupted hydrophobic packing of the αC-β4 loop and several of the methyl groups displaying the excited state in the CEST experiments are near the αC-β4 loop, we compared the CS differences of the methyl group from CEST experiments27 (ΔωExp) with those from the MD simulations (ΔωPred). We analyzed the methyl groups of L103, V104, I150, L172, and I180, sampling more than 500 snapshots as representative conformations of the kinase in the ES2 and GS states (Figure 3 – figure supplement 3). We then computed the distribution of the methyl 13C CS using ShiftX2.43 For Val104 and Ile150, these calculations yielded ΔωPred values of 0.87 ± 0.02 and 0.85 ± 0.02 ppm between the ES2 and GS conformations, respectively, which are in good agreement with ΔωExp (1.10 ± 0.06 and 0.94 ± 0.09 ppm) obtained from the CEST experiments or fitting the CPMG dispersion curves (Figure 5A-B). The ΔωPred values and the directions of the CS changes obtained for the remainder three sites are also in good agreement with the CEST profiles (Figure 5 – figure supplement 1). Figure 4C shows a linear fitting of ΔωPred and ΔωExp with a slope of 0.86 and R2 of 0.82. The latter supports the hypothesis that the excited state observed by NMR may correspond to the simulated structural ensemble of the ES2 basin. Finally, MSM estimates a population of ES2 of 6 ± 2%, which is consistent with the 5 ± 1% population found by fitting the NMR data.40 Altogether, the combination of NMR and MD simulations support the presence of an alternate inactive state that features disrupted hydrophobic interactions near the αC-β4 loop that perturb the anchoring with the αE helix, and in turn, the structural couplings between the two lobes of the kinase.
The F100A mutation disrupts the allosteric network of the kinase
The above analysis suggests that the destabilization of the αC-β4 loop may promote the GS to ES2 transition of the kinase. To validate this hypothesis, we generated the F100A mutant seeking to increase the flexibility of the αC-β4 loop. Starting from the coordinates of the X-ray structure of the ternary complex (PDB ID: 4WB5), we simulated the F100A dynamics in an explicit water environment. First, we performed a short equilibration using classical MD simulations. During 1 μs of MD simulations, the αC-β4 loop of the F100A mutant undergoes a significant motion as manifested by increased values of the backbone root mean square deviation (rmsd) and a conformational change (flipping) of the F102 side chain (Figure 6A). This region in PKA-CWT adopts a stable β-turn, with a persistent hydrogen bond between the backbone oxygen of F100 and the amide hydrogen of L103. In contrast, the hydrogen bond in F100A is formed more frequently between A100 and F102, resulting in an average γ-turn conformation (Figure 6B). Such local rearrangement not only disrupts the hydrogen bond between N99 and Y156, altering the interactions between the αC-β4 loop and the αE helix but also destabilizes the cation-π interaction between F102 and R308 at the C-tail (Figure 6B – supplementary figure 1)
The structural changes of the αC-β4 loop caused by the F100A mutation propagate through the protein core across the R spine, C spines, and shell residues (Figure 7A), and alter the response of the kinase to nucleotide binding. In the PKA-CWT, ATP binding shifts the conformational ensemble of the kinase toward an intermediate state competent to substrate binding, as represented by the changes in the population densities of the hydrophobic core residues as a function of the rmsd (Figure 7B). For F100A, the population of the C spine follows the same trend of the WT kinase, with the shell residues already populating the intermediate state. In contrast, the R spine residues fail to adopt a competent state (Figure 7B). The latter is due to the perturbation of hydrophobic packing of L95 and L106 in the R spine, and V104 in the shell near the mutation site. Using the lowest principal components, we also analyzed the global dynamic response of the kinase to ATP binding. Not only does F100A change the breathing mode of the two lobes (PC1), but it also alters the shearing motion (PC2) of the binary complex, emphasizing the importance of this allosteric site for the inter-lobe communication (Figure 7C-E).
To analyze the internal correlated motions of the kinase, we further analyzed the MD trajectories of both wild-type and mutant enzyme using mutual information (Figure 8).44 This method identifies correlated changes of backbone and side chain rotamers in proteins and detects clusters of residues that are responsible for intramolecular communication between active sites.44 For PKA-CWT, we observed numerous strong correlations that are clustered within each lobe and across the enzyme, connecting key motifs, such as the Gly-rich loop, αC-β4 loop, activation loop, and C-terminal tail (Figure 8A). These correlation patterns are reminiscent of the chemical shift correlations experimentally found for PKA-CWT,20 which constitute the hallmark of the dynamically-committed state.35,36 In contrast, the mutual information matrix for the F100A mutant displays weaker and interspersed correlations across the enzyme. New correlations also are present across the helical C-lobe suggesting a rewiring of the internal allosteric communication (Figure 8B). A possible explanation is that the increased motion caused by the elimination of the F100 aromatic side chain at the αC-β4 loop propagates throughout the entire hydrophobic core. This structural reorganization causes the F100A kinase to adopt a dynamically uncommitted state.35 Taken together, the simulations of F100A suggest that perturbations of the aC-b4 loop increase the local flexibility and break the structural connection between the two lobes, disrupting the cor-related breathing/shearing motions highlighted in previous MD simulation studies45,46. Although the single F100A mutation does not drive the kinase into a completely inactive state (ES2), it is sufficient to abolish the structural couplings between the two lobes.
Effects of the F100A mutation on the catalytic efficiency and binding thermodynamics of PKA-C
Based on the MD simulations of the F100A mutant, we hypothesized that the disruption of the structural coupling between the N-lobe and C-lobe would affect the binding cooperativity between nucleotide (ATPγN) and pseudosubstrate inhibitor (PKI5-24). To test this hypothesis, we expressed, purified, and evaluated the catalytic efficiency for the WT and F100A kinases by carrying out steady-state coupled enzyme assays47 using the standard substrate Kemptide.48 F100A showed only a slight increase in KM and Vmax compared to PKA-CWT, resulting in a small reduction of the catalytic efficiency (kcat/KM = 0.50 ± 0.04 and 0.41 ± 0.08 for PKA-CWT and PKA-CF100A, respectively; Supplementary table 2). We then performed ITC49 to obtain ΔG, ΔH, -TΔS, and Kd, and determine the cooperativity coefficient (σ) for ATPγN and PKI5-24 binding. We first analyzed the binding of ATPγN to the apo PKA-CF100A and, subsequently, the binding of PKI5-24 to the ATPγN-saturated PKA-CF100A (Supplementary Table 3). We found that PKA-CWT and PKACF100A have similar binding affinities for ATPγN (Kd = 83 ± 8 μM and 73 ± 2 μM, respectively). However, in the apo form, F100A showed a 3-fold higher binding affinity for the pseudosubstrate relative to PKA-CWT (Kd = 5 ± 1 μM and 17 ± 2 μM, respectively - Supplementary Table 3). Upon saturation with ATPγN, PKA-CF100A displayed a 12-fold reduction in binding affinity for PKI5-24, resulting in a σ of ∼3, a value significantly lower than the WT enzyme (σ greater than 100). These data support the predictions of MD simulations and mutual information that ATP binding to F100A does not promote a conformational state fully competent with substrate binding.
NMR mapping of nucleotide/pseudosubstrate binding response
To elucidate the atomic details of the disrupted structural coupling and binding cooperativity for PKA-CF100A, we used solution NMR spectroscopy mapping its response to the nucleotide (ATPγN) and pseudosubstrate (PKI5-24) binding. Specifically, we monitored the 1H and 15N chemical shift perturbations (CSPs, Δδ) of the amide fingerprints for the binary (PKA-CF100A/ATPγN) and ternary (PKA-CF100A/ATPγN/PKI5-24) complexes and compared them with the corresponding complexes of PKA-CWT (Figure 9 and Figure 9 – figure supplement 1). The interaction with ATPγN causes a dramatic broadening of several amide resonances throughout the structure of PKA-CF100A, suggesting the presence of an intermediate conformational exchange as previously observed.50 For the remainder of the amide peaks, both PKA-CWT and PKA-CF100A exhibit similar CSP profiles upon binding ATPγN (Figure 9A-B), although the extent of the CS changes is significantly attenuated for F100A (Figure 9 – figure supplement 1B). This reduction is apparent for residues in the N-lobe (β2-β3 region), around the mutation site (β4), and at the C-terminal tail. A similar pattern emerges upon binding PKI5-24 to the ATPγN-saturated PKA-CF100A, where a substantial decrease in CSP is observed for residues of the C-lobe localized near the motifs critical for substrate-binding (i.e., αE, αF, and αG helices, Figure 9C-D).
To determine the global response to ligand binding for WT and F100A, we analyzed the chemical shifts of the amide fingerprints of the two proteins using the COordiNated ChemIcal Shifts bEhavior (CONCISE).51 CONCISE describes the overall changes of the protein fingerprint resonances by providing the probability density (population) of each state along the conformational equilibrium for binding phenomena that follow linear CS trajectories.51 For PKA-CWT, both nucleotide and pseudosubstrate binding shift the overall population of the amides from an open to an intermediate and a fully closed state (Figure 9 – figure supplement 2). A similar trend is observed for PKA-CF100A, though the probability densities are broader, indicating that the amide resonances follow a less coordinated response.51 Also, the maximum of the probability density for the ternary complex is shifted toward the left, indicating that the mutant adopts a more open structure than the corresponding wild-type enzyme. Overall, the shapes and the positions of the probability distributions of the resonances suggests that several residues do not respond in a coordinated manner to the nucleotide binding, and PKI5-24 shifts the conformation of the kinase toward a partially closed state, which explains the loss in binding cooperativity as previously observed.21-
Discussion
Our structural and dynamic studies suggest that the binding cooperativity of PKA-C originates from the allosteric coupling between the nucleotide-binding pocket and the interfacial region between the two lobes, which harbors the substrate binding cleft.21–23 The latter has been emphasized for other kinases such as Src and ERK2.56,57 The biological significance of this intramolecular communication has emerged from our recent studies on pathological mutations situated in the activation loop of PKA-C and linked to Cushing’s syndrome.21–23 These mutations drastically reduce substrate binding affinity, and more importantly, disrupt the communication between the two ligand binding pockets of the kinase.21–23 Interestingly, the E31V mutation distal from the active site and also related to the Cushing’s syndrome, displays disfunction similar to the other orthosteric mutations, suggesting that it is possible to modulate substrate recognition allosterically.23 Indeed, these mutations are connected to allosteric nodes that, once perturbed, radiate their effects towards the periphery of the enzyme interrupting the coordinated dynamic coupling between the two lobes.23 Interestingly, these mutations do not prevent Kemptide phosphorylation; rather, they cause a loss of substrate fidelity with consequent aberrant phosphorylation of downstream substrates.10,58–62 Additionally, recent thermodynamic and NMR analyses of PKA-C binding to different nucleotides and inhibitors demonstrated that it is possible to control substrate binding affinity by changing the chemistry of the ligand at the ATP binding pocket.22 Altogether, these studies put forward the idea that substrate phosphorylation and cooperativity between ATP and substrates may be controlled independently.
By integrating our NMR data with RAM simulations and MSM, we comprehensively mapped the free energy landscape of PKA-C in various forms. We found that the active kinase unleashed from the regulatory subunits occupies a broad energy basin (GS) that corresponds to the conformation of the ternary structure of PKA-C with ATP and pseudosubstrate (PKI5-24), exemplifying a catalytic competent state poised for phosphoryl transfer.15 We also identified two orthogonal conformationally excited states, ES1 and ES2. While ES1 corresponds to canonical inactive kinase conformations, the ES2 state was never observed in the crystallized structures. Our previous CEST NMR measurements suggested the presence of an additional sparsely populated state that, at that time, we were unable to characterize structurally.27 The chemical shifts of this state did not follow the opening and closing of the kinase active cleft and were assigned to a possible alternate inactivation pathway.27 These new simulations and MSM show that this sparsely populated state may be attributed to the transition from GS to ES2, which features a disruption of hydrophobic packing, with a conformational rearrangement for the αC-β4 loop that causes a partial disruption of the hydrophobic R spine. These structural changes rewire the allosteric coupling between the two lobes, as shown by mutual information analysis. A single mutation (F100A), suggested by our simulations, promotes the flip of the αC-β4 loop and reproduces the hypothesized structural un-coupling between the two lobes of PKA-C. We experimentally tested the effects of the F100A mutation and found that it prevents the enzyme from adopting a conformation competent for substrate binding, resulting in a drastic reduction of the cooperativity between ATP and substrate. These NMR data further support our working model, showing that if the inter-lobe communication is interrupted, the binding response of the kinase is attenuated. Indeed, our investigation further emphasizes the integral role of the kinase hydrophobic core and demonstrates that alterations in the spines and shell residues may lead to a dysfunctional kinase by either preventing phosphorylation (see V104G and I150A mutations)40 or by disrupting the binding cooperativity as for F100A.
The αC-β4 loop is a regulatory element present in all EPKs and its importance has been stressed in bioinformatics studies63 and supported by computational work.55 F100 and F102 in the αC-β4 loop are at intersections of various structural communities and constitute a critical hydrophobic node that anchors the N- to the C-lobe. Also, studies on EGFR and ErbB2 kinases led to the hypothesis that the αC-β4 loop may act as a molecular brake 64,65 or an autoinhibitory switch.66 Therefore, it is not surprising that activating mutations and in-frame insertions in the αC-β4 loop are frequently found in kinase-related cancers and somatic oncogenic mutations such as P101S, P101L, and L103F.67 Also, an elegant study by Kannan and coworkers emphasized the role of the αC-β4 loop in dimerization and aberrant activation of EGFR by insertion mutations.28 Moreover, in a recent paper Zhang and co-workers found a similar behavior for the HER2 exon 20 insertions.30 These latter two studies propose that in-frame insertions alter the conformational landscape of the kinase rigidifying the structure around the αC-β4 loop and restricting the kinase conformational ensemble in a constitutively active state. Remarkably, in the case of EGRF insertion mutations, the activation of the kinase occurs in a length-dependent manner. 28 Additionally, these human mutations display a gradual response to drugs, which can be exploited for designing selective inhibitors against EGFR pathological mutants.68
The data presented here show that it is possible to abolish the binding cooperativity of a kinase by turning the dial in the opposite direction, i.e., increasing the flexibility of the αC-β4 loop and disconnecting the allosteric network between the N- and C-lobes. The identification of this new, partially inactivating pathway provides further understanding of how to control the dynamics and function of kinases.
Material and methods
Replica-averaged metadynamics (RAM) simulations
I. System setup
As a starting structure, we used the coordinates of the crystal structure of PKA-CWT (PDB ID: 1ATP) and added the missing residues 1-14 at the N terminus. The protonation state of histidine residues was set as previously reported.69 The kinase was solvated in a rhombic dodecahedron solvent box using the three-point charge TIP3P model70 for water, which extended approximately 10 Å away from the surface of the proteins. Counter ions (K+ and Cl-) were added to ensure electrostatic neutrality to a final ionic concentration of ∼ 150 mM. All protein covalent bonds were constrained using the LINCS algorithm71 and long-range electrostatic interactions were simulated using the particle-mesh Ewald method with a real-space cut-off of 10 Å.72 The simulations of apo, binary (one Mg2+ ion and one ATP), and ternary (two Mg2+ ions, one ATP, and PKI5-24 ) forms were performed simultaneously using GROMACS 4.673 with CHARMM36a1 force field.74 For F100A, the corresponding residues were mutated through the mutagenesis wizard of Pymol.75
II. Standard MD simulations
Each system was minimized using the steepest descent algorithm to remove the geometric distortions, and then gradually heated to 300 K at a constant volume over 1 ns using harmonic restraints with a force constant of 1000 kJ/(mol *Å2) on heavy atoms for both the kinase and nucleotide. The restraints were gradually released over the following 12 ns of simulations at constant pressure (1 atm) and temperature (300 K). The systems were equili-brated for an additional 20 ns without positional restraints. A Parrinello-Rahman barostat76 was used to keep the pressure constant, while a V-rescale thermostat 77 with a time step of 2 fs was used to keep the temperature constant. Each system was simulated for 1.05 µs, with snapshots recorded every 20 ps.
III. Replica exchange (REX) simulations
Following standard MD simulations, parallel replica exchange (REX) simulations were set up on the apo, binary, and ternary forms of PKA-C. Four replicas were used for each REX simulation, and the initial structures were randomly chosen from the µs-scale unbiased simulations. Chemical shifts of PKA-C for N, CA, CO, CB, and HN from NMR experiments were imposed as restraints based on the following penalty function:
where α is the force constant, k runs over all residues of the protein, l denotes the different backbone atoms, and m runs over the four replicas. was computed using CamShift, a module of ALMOST-2.1.78 The force constant was gradually increased from 0 (unbiased) to 20 (maximum restraints for production) over 50 ns. All other settings were identical to the unbiased simulations. REX simulations were carried out with GROMACS 4.6,73 with the replica exchange controlled by the PLUMED 2.1.1 module.79 Approximately, 100 ns of REX simulations were further carried out for each replica in the three PKA-C forms.
IV. Replica-averaged metadynamics (RAM) simulations
RAM simulations were started from the final structures of REX simulations. The CS restraints were imposed in the same manner as for the REX simulations. Four CVs are chosen to increase the conformational plasticity around the catalytic cores (detailed in Figure S1): (CVI) the ψ angles of the backbone of all the loops that are not in contact with ATP (Back-far), (CVII) the ϕ angles of the backbone of all the loops that are in contact with ATP (Back-close), (CVIII) the χ1 angles of side-chains of all the loops that are in contact with ATP (Side-close), (CVIV) the radius of gyration calculated over the rigid part (i.e., residues 50-300) of the protein (rgss). Gaussian deposition rate was performed with an initial rate of 0.125 kJ/mol/ps, where the σ values were set to 0.5, 0.2, 0.2, and 0.01 nm for the four CVs, respectively. The RAM simulations were also carried out with GROMACS 4.6 in conjunction with PLUMED 2.1 and ALMOST 2.1, and continued for ∼400 ns for each replica with exchange trails every 1 ps.
V. Reconstruction of free energy surface (FES) from the RAM simulations
After about 300 ns, the sampling along the first 3 CVs reached convergence with fluctuations of bias within 1 kcal/mol, allowing a reliable reconstruction of the corresponding FES. The production run was continued for another 100 ns to sample enough conformations. These conformations were first clustered into microstates using the regular spatial method (cut-off radius of 0.13), and then the free energy of each state was reweighted according to the deposited potential along each CV, which can be obtained from the analysis module of METAGUI.80 To visualize the distribution of these microstates and their relative energy differences, we further performed principal component analysis to project the microstates represented in the 3-dimensional CV space into 2-dimensional space spanned by PC1 and PC2. Then we plotted these microstates in a 3-dimensional space spanned by PC1, PC2, and free energy differences ΔG.
VI. Independent validation of chemical shifts with SPARTA+
During the REX and RAM simulations, the chemical shifts were computed via CamShift in ALMOST 2.1.78 As an independent validation for the efficacy of the bias, we further calibrated the chemical shifts of 2000 snapshots with MDTraj81 and SPARTA+.34
VII. Adaptive sampling
The first round of adaptive sampling started from the snapshots of low-energy microstates obtained from the previous step, i.e., 1200 structures for the apo form, 400 structures for the binary form, and 200 structures for the ternary form. The initial velocities were randomly generated to satisfy the Maxwell distribution at 300K. For the apo form, a 10 ns simulation was performed for each run, whereas for the binary, each simulation lasted 30 ns, resulting in a total of 12 μs trajectories for both the apo and binary forms. Three rounds of adaptive sampling were started from the 400 microstates, derived via K-mean clustering of all snapshots of previous ensembles, to obtain a converged free energy landscape. Therefore, a total of 100 µs trajectories and 100,0000 snapshots (100 ps per frame) were collected for both the apo and binary form after three rounds of adaptive sampling, and a total of 60 µs trajectories were collected for the ternary form.
VIII. Markov state model (MSM) and time-lagged independent component analysis (tICA)
The Cartesian coordinates of key hydrophobic residues, including R-spine residues, L95, L106, Y164, and F185, and the shell residues, V104, M118, and M120, were chosen as the metrics to characterize the conformational transition of the hydrophobic core of PKA-C. Specifically, each snapshot was first aligned to the same reference structure by superimposing the αE (residues 140-160) and αF helices (residues 217-233) and represented by the deviation of the Cartesian coordinates of key residues. The representation in this metric space was further reduced to 10-dimension vectors using time-lagged independent component analysis (tICA)82 at a lag time of 1 ns. All snapshots were clustered into 400 microstates with K-mean clustering. An MSM was built upon the transition counts between these microstates.
IX. Kinetic Monte Carlo trajectory of PKA-C in different forms
Long trajectories were generated using a kinetic Monte Carlo method based on the MSM transition probability matrix of the three forms of PKA-C. Specifically, the discrete jumps between the 100 microstates were sampled for 60 μs. Random conformations were chosen for each state from all the available snapshots. Subsequently, the time series of various order parameters were analyzed and plotted.
Protein expression and purification
The recombinant human Cα subunit of cAMP-dependent protein kinase with the Phe to Ala mutation in position 100 (PKA-CF100A) was generated from the human PKA-Cα wild-type using Quik-Change Lightning mutagenesis kit (Agilent genomics). The key resource table lists the PCR primers used to modify the pET-28a expression vector encoding for the wild-type human PKA-Cα gene (PRKACA – uniprot P17612) 21–23. The unlabeled and uniformly 15N-labeled PKA-CF100A mutant was expressed and purified following the same protocols used for the wild-type protein.41 Briefly, transformed E. coli BL21 (DE3) pLyss cells (Agilent) were cultured overnight at 30 °C in Luria-Bertani (LB) medium. The next morning, the cells were transferred to fresh LB medium for the overexpression of the unlabeled protein or to M9 minimal medium supplied with 15NH4Cl (Cambridge Isotope Laboratories Inc.) as the only nitrogen source for the labeled protein overexpression. In both cases, protein overexpression was induced with 0.4 mM of β-D-thiogalactopyranoside (IPTG) and carried out for 16 hours at 20 °C. The cells were harvested by centrifugation and resuspended in 50 mM Tris-HCl, 30 mM KH2PO4, 100 mM NaCl, 5 mM 2-mercaptoethanol, 0.15 mg/mL lysozyme, 200 μM ATP, DNaseI, 1 tablet of cOmplete ULTRA EDTA-free protease inhibitors (Roche Applied Science) (pH 8.0) and lysed using French press at 1,000 psi. The cell lysate was cleared by centrifugation (60,000 × g, 4 °C, 45 min), and the supernatant was batch-bound with Ni2+-NTA agarose affinity resin (Thermo Scientific). The his-tagged PKA-CF100A was eluted with 50 mM Tris-HCl, 30 mM KH2PO4, 100 mM NaCl, 5 mM 2-mercaptoethanol, 0.5 mM phenylmethylsulfonyl fluoride (PMSF) (pH 8.0) supplied with 200 mM of imidazole. The tail of poly-His was cleaved using a stoichiometric amount of recombinant tobacco etch virus (TEV) protease in 20 mM KH2PO4, 25 mM KCl, 5 mM 2-mercaptoethanol, 0.1 mM PMSF (pH 6.5), over-night at 4 °C. The different phosphorylation states of PKA-CF100A were separated using a cation exchange column (HiTrap Q-SP, GE Healthcare Life Sciences) using a linear gradient of KCl in 20 mM KH2PO4 at pH6.5.83 The purified protein isoforms were then stored in phosphate buffer containing 10 mM dithiothreitol (DTT), 10 mM MgCl2, and 1 mM NaN3 at 4 °C. The protein purity was assessed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS–PAGE).
Peptide synthesis
The Kemptide (LRRASLG) and PKI5–24 (TTYADFIASGRTGRRNAIHD) peptides were synthesized using a CEM Liberty Blue microwave synthesizer using standard Fmoc chemistry. All peptides were cleaved with Reagent K (82.5% TFA, 5% phenol, 5% thioanisole, 2.5% ethanedithiol, and 5% water) for 3 h and purified using a semipreparative Supelco C18 reverse-phase HPLC column at 3 mL/min. The purified peptides were concentrated, lyophilized, and stored at −20 °C. Molecular weight and quantity were verified by MALDI-TOF and/or amino-acid analysis (Texas A&M University).
Isotherm titration calorimetry (ITC) measurements
PKA-CF100A was dialyzed into 20 mM MOPS, 90 mM KCl, 10 mM DTT, 10 mM MgCl2, and 1 mM NaN3 (pH 6.5) and concentrated using conical spin concentrator (10 kDa membrane cut-off, Millipore) to a solution at 80-100 μM, as confirmed by A280 = 55,475 M−1 cm−1. Approximately 300 μL of protein was used for each experiment, with 50 μL of 2 mM ATPγN and/or 1 mM PKI5–24 in the titrant syringe. All measurements were performed at 300 K in triplicates with a low-volume NanoITC (TA Instruments). The binding was assumed to be 1:1, and curves were analyzed with the NanoAnalyze software (TA Instruments) using the Wiseman isotherm 49
where d[MX] is the change in total complex relative to the change in total protein concentration, d[Xtot] is dependent on r (the ratio of Kd relative to the total protein concentration), and Rm (the ratio between total ligand and total protein concentration). The heat of dilution of the ligand into the buffer was considered for all experiments and subtracted.
The free energy of binding was determined from:
where R is the universal gas constant and T is the temperature at measurement (300 K). The entropic contribution to binding was calculated using:
The degree of cooperativity (σ) was calculated as:
where Kd apo is the dissociation constant of PKI5–24 binding to the apo-enzyme, and Kd Nucleotide is the corresponding dissociation constant for PKI5–24 binding to the nucleotide-bound kinase.
Enzyme assays
Steady-state coupled enzyme activity assays using Kemptide as substrate were performed under saturating ATP concentrations and spectrophotometrically at 298 K, as reported by Cook et al.47 The values of Vmax and KM were obtained from a nonlinear fit of the initial velocities to the Michaelis-Menten equation.
NMR spectroscopy
NMR measurements were performed on a Bruker Avance NEO spectrometer operating at a 1H Larmor frequency of 600 MHz equipped with a cryogenic probe or on a Bruker Avance III 850 MHz spectrometer equipped with a TCI cryoprobe. The NMR experiments were recorded at 300K in 20 mM KH2PO4 (pH 6.5), 90 mM KCl, 10 mM MgCl2, 10 mM DTT, 1 mM NaN3, 5% D2O, and 0.1% 4-benzene sulfonyl fluoride hydrochloride (AEBSF, Pefabl–c - Sigma-Aldrich). Concentrations for samples were 0.15 mM of uniformly 15N-labeled PKA-CF100A, as determined by A280 measurements, 12 mM ATPγN or ADP was added for the nucleotide-bound form, and 0.3 mM PKI5-24 for the ternary complex. [1H, 15N]-WADE-TROSY-HSQC pulse sequence84 was used to record the amide fingerprint spectra of PKA-CF100A in the apo, nucleotide-bound (ADP- or ATPγN-bound – binary form), and ternary complex (PKACF100A/ATPγN/PKI5-24). All [1H, 15N]-WADE-TROSY-HSQC experiments were acquired with 2048 (proton) and128 (nitrogen) complex points, processed using NMRPipe85 and visualized using NMRFAM-SPARKY 86 and POKY.87 Combined chemical shift perturbations (CSP) were calculated using 1H and 15N chemical shifts according to:
in which Δδ is the CSP; ΔδH and ΔδN are the differences of 1H and 15N chemical shifts, respectively, between the first and last point of the titration; and 0.154 is the scaling factor for nitrogen.88 COordiNated ChemIcal Shift bEhavior (CONCISE). The normal distributions reported in the CONCISE plot were calculated using principal component analysis for residues whose chemical shifts responded linearly to ligand binding.51 In this work, we use the 1H and 15N chemical shifts derived from the [1H, 15N]-WADE-TROSY-HSQC experiments for the apo, ADP, ATPγN, and ATPγN/PKI5–24 bound forms of PKA-C.
CHEmical Shift Covariance Analysis (CHESCA)
This analysis was used to identify and characterize allosteric networks of residues showing concerted responses to nucleotide and pseudosubstrate binding. To identify inter-residue correlations, four states were used: apo, ATPγN-bound, ADP-bound, and ATPγN/PKI5–24. The identification of inter-residue correlations by CHESCA relies on agglomerative clustering (AC) and singular value decomposition (SVD).53 Pairwise correlations of chemical shift changes were calculated to identify networks. When plotted on a correlation matrix, the 2D correlations allow the identification of local and long-range correlated residues. For each residue, the maximum change in the chemical shift was calculated in both the 1H (x) and 15N (y) dimensions (Δδx,y). The residues included in the CHESCA analysis were those satisfying Δδx,y > ½ ΔνxA,yA+½ ΔνxB,yB, where A and B correspond to two different forms analyzed, and Δν denotes the linewidth. Correlation scores were used to quantify the CHESCA correlation of a single residue or a group of residues with another group. Correlation scores were evaluated for single residues or for the entire protein using the following expression:
where Rij is the correlation coefficient and i and j denote (a) a single residue and the remainder residues of the protein, respectively, or (b) both represent all the assigned residues in the entire protein. For all the analyses a Rcutoff of 0.98 was used. Community CHESCA analysis utilizes a similar approach to map correlations between functional communities within the kinase. Each community represents a group of residues associated with a function or regulatory mechanism.55 To represent community-based CHESCA analysis, we utilized Rcutoff > 0.8. For instance, to represent a chemical shift correlation between two community (X and Y) with nA and nB number of assigned residues, respectively, the correlation score between A and B is defined as
where Rij is the correlation coefficient between residue i (belonging to community A) and residue j (belonging to community B). Rcutoff is the correlation value cutoff. RA,B assumes values from 0 (no correlation between residues in A and B) to 1 (all residues in A have correlations > cut-off with all residues in B).
Acknowledgements
This work was supported by the National Institute of Health GM 100310 and HL 144130 to GV. The authors would like to acknowledge the Minnesota Supercomputing Institute for MD calculations. YW would like to thank the Guangdong Pearl River Talent Program (2021QN02Y618) and the National Natural Science Foundation of China (22007069, 92269102) for part of the MD analysis carried out at the Shenzhen Bay Laboratory Supercomputing Centre.
Competing interest
The authors declare no financial and non-financial competing interests.
Data Availability
All the data generated or analyzed in this study are included in the manuscript and supporting files. The NMR chemical shifts and MD trajectories are deposited in the Data Repository for the University of Minnesota, https://hdl.handle.net/11299/261043.
References
- 1.The Protein Kinase Complement of the Human GenomeScience 298:1912–1934
- 2.The Conformational Plasticity of Protein KinasesCell 109:275–282
- 3.Crystal structure of the catalytic subunit of cyclic adenosine monophosphate-dependent protein kinaseScience 253:407–14
- 4.Structure of a peptide inhibitor bound to the catalytic subunit of cyclic adenosine monophosphate-dependent protein kinaseScience 253:414–20
- 5.Assembly of allosteric macromolecular switches: lessons from PKANat Rev Mol Cell Biol 13:646–58
- 6.An adenosine 3’,5’-monophosphate-dependant protein kinase from rabbit skeletal muscleJ Biol Chem 243:3763–5
- 7.Local protein kinase A action proceeds through intact holoenzymesScience 356:1288–1293
- 8.Liberated PKA Catalytic Subunits Associate with the Membrane via Myristoylation to Preferentially Phosphorylate Membrane SubstratesCell Rep 19:617–629
- 9.KA regulation of neuronal function requires the dissociation of catalytic subunits from regulatory subunitsbioRxiv
- 10.Mislocalization of protein kinase A drives pathology in Cushing’s syndromeCell Rep 40
- 11.Crystallization studies of cAMP-dependent protein kinase. Cocrystals of the catalytic subunit with a 20 amino acid residue peptide inhibitor and MgATP diffract to 3.0 A resolutionJ Mol Biol 220:217–20
- 12.Dynamics of cAMP-dependent protein kinaseChem Rev 101:2243–70
- 13.Phosphoryl Transfer by Protein Kinase A Is Captured in a Crystal LatticeJournal of the American Chemical Society 135:4788–4798
- 14.Phosphoryl Transfer Reaction Snapshots in Crystals: INSIGHTS INTO THE MECHANISM OF PROTEIN KINASE A CATALYTIC SUBUNITJ Biol Chem 290:15538–15548
- 15.Zooming in on protons: Neutron structure of protein kinase A trapped in a product complexSci Adv 5
- 16.Surface comparison of active and inactive protein kinases identifies a conserved activation mechanismProc Natl Acad Sci U S A 103:17783–8
- 17.Conserved spatial patterns across the protein kinase familyBiochim Biophys Acta 1784:238–43
- 18.Kinase regulation by hydrophobic spine assembly in cancerMol Cell Biol 35:264–76
- 19.Synergistic binding of nucleotides and inhibitors to cAMP-dependent protein kinase examined by acrylodan fluorescence spectroscopyJ Biol Chem 272:1507–13
- 20.Globally correlated conformational entropy underlies positive and negative cooperativity in a kinase’s enzymatic cycleNat Commun 10
- 21.Defective internal allosteric network imparts dysfunctional ATP/substrate-binding cooperativity in oncogenic chimera of protein kinase ACommun Biol 4
- 22.Cushing’s syndrome driver mutation disrupts protein kinase A allosteric network, altering both regulation and substrate specificitySci Adv 5
- 23.Is Disrupted Nucleotide-Substrate Cooperativity a Common Trait for Cushing’s Syndrome Driving Mutations of Protein Kinase A?J Mol Biol 433
- 24.Assessment of the use of NMR chemical shifts as replica-averaged structural restraints in molecular dynamics simulations to characterize the dynamics of proteinsJ Phys Chem B 117:1838–43
- 25.Statistical mechanics of the denatured state of a protein using replica-averaged metadynamicsJ Am Chem Soc 136:8982–91
- 26.Markov State Models: From an Art to a ScienceJ Am Chem Soc 140:2386–2396
- 27.ATP-competitive inhibitors modulate the substrate binding cooperativity of a kinase by altering its conformational entropySci Adv 8
- 28.Altered conformational landscape and dimerization dependency underpins the activation of EGFR by alpha C-beta 4 loop insertion mutationsProceedings of the National Academy of Sciences of the United States of America 115
- 29.Analogous regulatory sites within the αC-β4 loop regions of ZAP-70 tyrosine kinase and AGC kinasesBiochimica Et Biophysica Acta-Proteins and Proteomics 1784:27–32
- 30.Conformational Landscapes of HER2 Exon 20 Insertions Explain Their Sensitivity to Kinase Inhibitors in Lung AdenocarcinomaJ Thorac Oncol 15:962–972
- 31.Using NMR chemical shifts as structural restraints in molecular dynamics simulations of proteinsStructure 18:923–33
- 32.Enhancing Biomolecular Simulations with Hybrid Potentials Incorporating NMR DataJ Chem Theory Comput 18:7733–7750
- 33.A bias-exchange approach to protein foldingJ Phys Chem B 111:4553–9
- 34.SPARTA+: a modest improvement in empirical NMR chemical shift prediction by means of an artificial neural networkJ Biomol NMR 48:13–22
- 35.Dynamically committed, uncommitted, and quenched states encoded in protein kinase A revealed by NMR spectroscopyProc Natl Acad Sci U S A 108:6969–74
- 36.Dynamics connect substrate recognition to catalysis in protein kinase ANat Chem Biol 6:821–8
- 37.Mapping the Hydrogen Bond Networks in the Catalytic Subunit of Protein Kinase A Using H/D Fractionation FactorsBiochemistry 54:4042–4049
- 38.Cloud-based simulations on Google Exacycle reveal ligand modulation of GPCR activation pathwaysNat Chem 6:15–21
- 39.Complete protein-protein association kinetics in atomic detail revealed by molecular dynamics simulations and Markov modellingNature Chemistry 9:1005–1011
- 40.A dynamic hydrophobic core orchestrates allostery in protein kinasesSci Adv 3
- 41.Allosteric cooperativity in protein kinase AProc Natl Acad Sci U S A 105:506–11
- 42.Activation pathway of Src kinase reveals intermediate states as targets for drug designNat Commun 5
- 43.SHIFTX2: significantly improved protein chemical shift predictionJournal of Biomolecular Nmr 50:43–57
- 44.Comparing Conformational Ensembles Using the Kullback-Leibler Divergence ExpansionJ Chem Theory Comput 8:2115–2126
- 45.600 ps molecular dynamics reveals stable substructures and flexible hinge points in cAMP dependent protein kinaseBiopolymers 50:513–24
- 46.Release of ADP from the catalytic subunit of protein kinase A: a molecular dynamics simulation studyProtein Sci 14:159–68
- 47.Adenosine cyclic 3’,5’-monophosphate dependent protein kinase: kinetic mechanism for the bovine skeletal muscle catalytic subunitBiochemistry 21:5794–9
- 48.Role of Multiple Basic Residues in Determining Substrate-Specificity of Cyclic Amp-Dependent Protein-KinaseJournal of Biological Chemistry 252:4888–4894
- 49.Rapid measurement of binding constants and heats of binding using a new titration calorimeterAnalytical Biochemistry 179:131–137
- 50.Synchronous Opening and Closing Motions Are Essential for cAMP-Dependent Protein Kinase A SignalingStructure 22:1735–43
- 51.NMR mapping of protein conformational landscapes using coordinated behavior of chemical shifts upon ligand bindingPhys Chem Chem Phys 16:6508–18
- 52.Signaling through dynamic linkers as revealed by PKAProc Natl Acad Sci U S A 110:14231–6
- 53.Mapping allostery through the covariance analysis of NMR chemical shiftsProc Natl Acad Sci U S A 108:6133–8
- 54.A conserved helix motif complements the protein kinase coreProc Natl Acad Sci U S A 90:10618–22
- 55.Dynamic architecture of a protein kinaseProc Natl Acad Sci U S A 111:E4623–31
- 56.A dynamically coupled allosteric network underlies binding cooperativity in Src kinaseNat Commun 6
- 57.Activation Loop Dynamics Are Coupled to Core Motions in Extracellular Signal-Regulated Kinase-2Biochemistry 59:2698–2706
- 58.Constitutive activation of PKA catalytic subunit in adrenal Cushing’s syndromeN Engl J Med 370:1019–28
- 59.Mechanisms of Aberrant PKA Activation by Cα Subunit MutationsHorm Metab Res 49:307–314
- 60.PKA catalytic subunit mutations in adrenocortical Cushing’s adenoma impair association with the regulatory subunitNature Communications 5
- 61.Novel somatic mutations in the catalytic subunit of the protein kinase A as a cause of adrenal Cushing’s syndrome: a European multicentric studyJ Clin Endocrinol Metab 99:E2093–100
- 62.Cushing’s syndrome mutant PKA(L)(205R) exhibits altered substrate specificityFEBS Lett 591:459–467
- 63.Did protein kinase regulatory mechanisms evolve through elaboration of a simple structural component?Journal of Molecular Biology 351:956–972
- 64.A molecular brake in the kinase hinge region regulates the activity of receptor tyrosine kinasesMol Cell 27:717–30
- 65.Structural and dynamic insights into the energetics of activation loop rearrangement in FGFR1 kinaseNature Communications
- 66.Mutational activation of ErbB2 reveals a new protein kinase autoinhibition mechanismJ Biol Chem 283:1588–1596
- 67.ProKinO: a unified resource for mining the cancer kinomeHum Mutat 36:175–86
- 68.Emerging roles of the alphaC-beta4 loop in protein kinase structure, function, evolution, and diseaseIUBMB Life 72:1189–1202
- 69.Conformational equilibrium of N-myristoylated cAMP-dependent protein kinase A by molecular dynamics simulationsBiochemistry 51:10186–96
- 70.Comparison of Simple Potential Functions for Simulating Liquid WaterJournal of Chemical Physics 79:926–935
- 71.LINCS: A linear constraint solver for molecular simulationsJournal of Computational Chemistry 18:1463–1472
- 72.Particle Mesh Ewald - an N.Log(N) Method for Ewald Sums in Large SystemsJournal of Chemical Physics 98:10089–10092
- 73.GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular SimulationJ Chem Theory Comput 4:435–47
- 74.Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone phi, psi and side-chain chi(1) and chi(2) dihedral anglesJ Chem Theory Comput 8:3257–3273
- 75.The PyMOL molecular graphics system, version 0.99 rc6. DeLano Scientific: San CarlosCA
- 76.Crystal-Structure and Pair Potentials - a Molecular-Dynamics StudyPhysical Review Letters 45:1196–1199
- 77.Canonical sampling through velocity rescalingJ Chem Phys 126
- 78.Fast and accurate predictions of protein NMR chemical shifts from interatomic distancesJ Am Chem Soc 131:13894–5
- 79.PLUMED: A portable plugin for free-energy calculations with molecular dynamicsComputer Physics Communications 180:1961–1972
- 80.METAGUI. A VMD interface for analyzing metadynamics and molecular dynamics simulationsComputer Physics Communications 183:203–211
- 81.MDTraj: A Modern Open Library for the Analysis of Molecular Dynamics TrajectoriesBiophys J 109:1528–32
- 82.Improvements in Markov State Model Construction Reveal Many Non-Native Interactions in the Folding of NTL9J Chem Theory Comput 9:2000–2009
- 83.Identification of phosphorylation sites in the recombinant catalytic subunit of cAMP-dependent protein kinaseJ Biol Chem 268:18626–32
- 84.Design and applications of water irradiation devoid RF pulses for ultra-high field biomolecular NMR spectroscopyPhysical Chemistry Chemical Physics 24:18477–18481
- 85.NMRPipe: a multidimensional spectral processing system based on UNIX pipesJ Biomol NMR 6:277–93
- 86.NMRFAM-SPARKY: enhanced software for biomolecular NMR spectroscopyBioinformatics 31:1325–1327
- 87.POKY software tools encapsulating assignment strategies for solution and solid-state protein NMR dataJ Struct Biol X 6
- 88.Using chemical shift perturbation to characterise ligand bindingProg Nucl Magn Reson Spectrosc 73:1–16
Article and author information
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Copyright
© 2023, Olivieri 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
- views
- 621
- downloads
- 56
- citations
- 2
Views, downloads and citations are aggregated across all versions of this paper published by eLife.