A transition path ensemble study reveals a linchpin role for Mg(2+) during rate-limiting ADP release from protein kinase A

Crystal Structures Reveal Hidden Domain Mechanics in Protein Kinase A (PKA)

Cyclic-AMP-dependent protein kinase A (PKA) is a critical enzyme involved in various signaling pathways that plays a crucial role in regulating cellular processes including metabolism, gene transcription, cell proliferation, and differentiation. In this study, the mechanisms of allostery in PKA were investigated by analyzing the vast repertoire of crystal structures available in the RCSB database. From existing structures of murine and human PKA, we elucidated the conformational ensembles and protein dynamics that are altered in a ligand-dependent manner. Distance metrics to analyze conformations of the G-loop were proposed to delineate different states of PKA and were compared to existing structural metrics. Furthermore, ligand-dependent flexibility was investigated through normalized B'-factors to better understand the inherent dynamics in PKA. The presented study provides a contemporary approach to traditional methods in engaging the use of crystal structures for understanding protein dynamics. Importantly, our studies provide a deeper understanding into the conformational ensemble of PKA as the enzyme progresses through its catalytic cycle. These studies provide insights into kinase regulation that can be applied to both PKA individually and protein kinases as a class.

Edmond Fischer's kinase legacy: History of the protein kinase inhibitor and protein kinase A

Although Fischer's extraordinary career came to focus mostly on the protein phosphatases, after his co-discovery of Phosphorylase Kinase with Ed Krebs he was clearly intrigued not only by cAMP-dependent protein kinase (PKA), but also by the heat-stable, high-affinity protein kinase inhibitor (PKI). PKI is an intrinsically disordered protein that contains at its N-terminus a pseudo-substrate motif that binds synergistically and with high-affinity to the PKA catalytic (C) subunit. The sequencing and characterization of this inhibitor peptide (IP20) were validated by the structure of the PKA C-subunit solved first as a binary complex with IP20 and then as a ternary complex with ATP and two magnesium ions. A second motif, nuclear export signal (NES), was later discovered in PKI. Both motifs correspond to amphipathic helices that convey high-affinity binding. The dynamic features of full-length PKI, recently captured by NMR, confirmed that the IP20 motif becomes dynamically and sequentially ordered only in the presence of the C-subunit. The type I PKA regulatory (R) subunits also contain a pseudo-substrate ATPMg2-dependent high-affinity inhibitor sequence. PKI and PKA, especially the Cβ subunit, are highly expressed in the brain, and PKI expression is also cell cycle-dependent. In addition, PKI is now linked to several cancers. The full biological importance of PKI and PKA signaling in the brain, and their importance in cancer thus remains to be elucidated.

The "violin model": Looking at community networks for dynamic allostery

Allosteric regulation of proteins continues to be an engaging research topic for the scientific community. Models describing allosteric communication have evolved from focusing on conformation-based descriptors of protein structural changes to appreciating the role of internal protein dynamics as a mediator of allostery. Here, we explain a "violin model" for allostery as a contemporary method for approaching the Cooper-Dryden model based on redistribution of protein thermal fluctuations. Based on graph theory, the violin model makes use of community network analysis to functionally cluster correlated protein motions obtained from molecular dynamics simulations. This Review provides the theory and workflow of the methodology and explains the application of violin model to unravel the workings of protein kinase A.

Molecular-dynamics simulation methods for macromolecular crystallography

It is investigated whether molecular-dynamics (MD) simulations can be used to enhance macromolecular crystallography (MX) studies. Historically, protein crystal structures have been described using a single set of atomic coordinates. Because conformational variation is important for protein function, researchers now often build models that contain multiple structures. Methods for building such models can fail, however, in regions where the crystallographic density is difficult to interpret, for example at the protein-solvent interface. To address this limitation, a set of MD-MX methods that combine MD simulations of protein crystals with conventional modeling and refinement tools have been developed. In an application to a cyclic adenosine monophosphate-dependent protein kinase at room temperature, the procedure improved the interpretation of ambiguous density, yielding an alternative water model and a revised protein model including multiple conformations. The revised model provides mechanistic insights into the catalytic and regulatory interactions of the enzyme. The same methods may be used in other MX studies to seek mechanistic insights.

Inflammatory macrophages exacerbate neutrophil-driven joint damage through ADP/P2Y1 signaling in rheumatoid arthritis

Rheumatoid arthritis (RA) is a chronic autoimmune disease that primarily affects the joints and is associated with excessive immune cell infiltration. However, the complex interactions between the immune cell populations in the RA synovium remain unknown. Here, we demonstrate that inflammatory macrophages in the synovium exacerbate neutrophil-driven joint damage in RA through ADP/P2Y₁ signaling. We show that extracellular ADP (eADP) and its receptors are obviously increased in synovial tissues of RA patients as well as collagen-induced arthritis (CIA) mice, and eADP enhances neutrophil infiltration into joints through macrophages producing the chemokine CXCL2, aggravating disease development. Accordingly, the arthritis mouse model had more neutrophils in inflamed joints following ADP injection, whereas P2Y₁ deficiency and pharmacologic inhibition restored arthritis severity to basal levels, suggesting a dominant role of ADP/P2Y₁ signaling in RA pathology. Moreover, cellular activity of ADP/P2Y₁-mediated CXCL2 production was dependent on the G_αq/Ca²⁺-NF-κB/NFAT pathway in macrophages. Overall, this study reveals a non-redundant role of eADP as a trigger in the pathogenesis of RA through neutrophil recruitment and disrupted tissue homeostasis and function.

Drugging the Undruggable: How Isoquinolines and PKA Initiated the Era of Designed Protein Kinase Inhibitor Therapeutics

In 1984, Japanese researchers led by the biochemist Hiroyoshi Hidaka described the first synthetic protein kinase inhibitors based on an isoquinoline sulfonamide structure (Hidaka et al. Biochemistry, 1984 Oct 9; 23(21): 5036-41. doi: 10.1021/bi00316a032). These led to the first protein kinase inhibitor approved for medical use (fasudil), an inhibitor of the AGC subfamily Rho kinase. With potencies strong enough to compete against endogenous ATP, the isoquinoline compounds established the druggability of the ATP binding site. Crystal structures of their protein kinase complexes, including with cAMP-dependent protein kinase (PKA), showed interactions that, on the one hand, could mimic ATP but, on the other hand, could be optimized for high potency binding, kinase selectivity, and diversification away from adenosine. They also showed the flexibility of the glycine-rich loop, and PKA became a major prototype for crystallographic and nuclear magnetic resonance (NMR) studies of protein kinase mechanism and dynamic activity control. Since fasudil, more than 70 kinase inhibitors have been approved for clinical use, involving efforts that progressively have introduced new paradigms of data-driven drug discovery. Publicly available data alone comprise over 5000 protein kinase crystal structures and hundreds of thousands of binding data. Now, new methods, including artificial intelligence techniques and expansion of protein kinase targeting approaches, together with the expiration of patent protection for optimized inhibitor scaffolds, promise even greater advances in drug discovery. Looking back to the time of the first isoquinoline hinge binders brings the current state-of-the-art into stark contrast. Appropriately for this Perspective article, many of the milestone papers during this time were published in Biochemistry (now ACS Biochemistry).

From structure to the dynamic regulation of a molecular switch: A journey over 3 decades

It is difficult to imagine where the signaling community would be today without the Protein Data Bank. This visionary resource, established in the 1970s, has been an essential partner for sharing information between academics and industry for over 3 decades. We describe here the history of our journey with the protein kinases using cAMP-dependent protein kinase as a prototype. We summarize what we have learned since the first structure, published in 1991, why our journey is still ongoing, and why it has been essential to share our structural information. For regulation of kinase activity, we focus on the cAMP-binding protein kinase regulatory subunits. By exploring full-length macromolecular complexes, we discovered not only allostery but also an essential motif originally attributed to crystal packing. Massive genomic data on disease mutations allows us to now revisit crystal packing as a treasure chest of possible protein:protein interfaces where the biological significance and disease relevance can be validated. It provides a new window into exploring dynamic intrinsically disordered regions that previously were deleted, ignored, or attributed to crystal packing. Merging of crystallography with cryo-electron microscopy, cryo-electron tomography, NMR, and millisecond molecular dynamics simulations is opening a new world for the signaling community where those structure coordinates, deposited in the Protein Data Bank, are just a starting point!

Druggable Transient Pockets in Protein Kinases

Drug discovery using small molecule inhibitors is reaching a stalemate due to low selectivity, adverse off-target effects and inevitable failures in clinical trials. Conventional chemical screening methods may miss potent small molecules because of their use of simple but outdated kits composed of recombinant enzyme proteins. Non-canonical inhibitors targeting a hidden pocket in a protein have received considerable research attention. Kii and colleagues identified an inhibitor targeting a transient pocket in the kinase DYRK1A during its folding process and termed it FINDY. FINDY exhibits a unique inhibitory profile; that is, FINDY does not inhibit the fully folded form of DYRK1A, indicating that the FINDY-binding pocket is hidden in the folded form. This intriguing pocket opens during the folding process and then closes upon completion of folding. In this review, we discuss previously established kinase inhibitors and their inhibitory mechanisms in comparison with FINDY. We also compare the inhibitory mechanisms with the growing concept of “cryptic inhibitor-binding sites.” These sites are buried on the inhibitor-unbound surface but become apparent when the inhibitor is bound. In addition, an alternative method based on cell-free protein synthesis of protein kinases may allow the discovery of small molecules that occupy these mysterious binding sites. Transitional folding intermediates would become alternative targets in drug discovery, enabling the efficient development of potent kinase inhibitors.

Molecular Insights into the Trapping Effect of Ca2+ in Protein Kinase A: A Molecular Dynamics Study

Protein kinase A has become a model system for the study of kinases, and therefore, a comprehensive understanding of the underlying molecular mechanisms in its catalytic cycle is of crucial importance. One of the aspects that has received recent attention is the role that metal cofactors play in the catalytic cycle. Although Mg²⁺ is the well-known physiological ion used by protein kinases, Ca²⁺ ions can also assist the phosphoryl transfer reaction but with lower catalytic activities. This inhibitory effect has been attributed to the ability of Ca²⁺ to trap the reaction products at the active site, and it has been proposed as a possible regulatory mechanism of the enzyme. Thus, in order to get a clearer understanding of these molecular events, computational simulations in the product state of PKA, in the presence of Mg²⁺ and Ca²⁺ ions, were performed through molecular dynamics (MD). Different protonation states of the active site were considered in order to model the different mechanistic pathways that have been proposed. Our results show that different protonation states of the phosphorylated serine residue at the peptide substrate (pSer21), as well as the protonation state of residue Asp166, can have a marked influence on the flexibility of regions surrounding the active site. This is the case of the glycine-rich loop, a structural motif that is directly involved in the release of the products from the PKA active site. MD simulations were capable to reproduce the crystallographic conformations but also showed other conformations not previously reported in the crystal structures that may be involved in enhancing the affinity of pSP20 to PKA in the presence of Ca²⁺. Hydrogen bonding interactions at the PKA-pSP20 interface were influenced whether by the protonation state of the active site or by the metal cofactor used by the enzyme. Altogether, our results provide molecular aspects into the inhibitory mechanism of Ca²⁺ in PKA and suggest which is the most probable protonation state of the phosphorylated product at the active site.

Two PKA RIα holoenzyme states define ATP as an isoform-specific orthosteric inhibitor that competes with the allosteric activator, cAMP

Protein kinase A (PKA) holoenzyme, comprised of a cAMP-binding regulatory (R)-subunit dimer and 2 catalytic (C)-subunits, is the master switch for cAMP-mediated signaling. Of the 4 R-subunits (RIα, RIβ, RIIα, RIIβ), RIα is most essential for regulating PKA activity in cells. Our 2 RIα2C2 holoenzyme states, which show different conformations with and without ATP, reveal how ATP/Mg2+ functions as a negative orthosteric modulator. Biochemical studies demonstrate how the removal of ATP primes the holoenzyme for cAMP-mediated activation. The opposing competition between ATP/cAMP is unique to RIα. In RIIβ, ATP serves as a substrate and facilitates cAMP-activation. The isoform-specific RI-holoenzyme dimer interface mediated by N3A-N3A' motifs defines multidomain cross-talk and an allosteric network that creates competing roles for ATP and cAMP. Comparisons to the RIIβ holoenzyme demonstrate isoform-specific holoenzyme interfaces and highlights distinct allosteric mechanisms for activation in addition to the structural diversity of the isoforms.

Zooming in on protons: Neutron structure of protein kinase A trapped in a product complex

The question vis-à-vis the chemistry of phosphoryl group transfer catalyzed by protein kinases remains a major challenge. The neutron diffraction structure of the catalytic subunit of cAMP-dependent protein kinase (PKA-C) provides a more complete chemical portrait of key proton interactions at the active site. By using a high-affinity protein kinase substrate (PKS) peptide, we captured the reaction products, dephosphorylated nucleotide [adenosine diphosphate (ADP)] and phosphorylated PKS (pPKS), bound at the active site. In the complex, the phosphoryl group of the peptide is protonated, whereas the carboxyl group of the catalytic Asp166 is not. Our structure, including conserved waters, shows how the peptide links the distal parts of the cleft together, creating a network that engages the entire molecule. By comparing slow-exchanging backbone amides to those determined by the NMR analysis of PKA-C with ADP and inhibitor peptide (PKI), we identified exchangeable amides that likely distinguish catalytic and inhibited states.

Divalent metal ions control activity and inhibition of protein kinases

Protein kinases are key enzymes in the regulation of eukaryotic signal transduction. As metalloenzymes they employ divalent cations for catalysis and regulation. We used the catalytic (C) subunit of cAMP-dependent protein kinase (PKA) as a model protein to investigate the role of a variety of physiologically or pathophysiologically relevant divalent metal ions in distinct steps within the catalytic cycle. It is established that divalent metal ions play a crucial role in co-binding of nucleotides and also assist in catalysis. Our studies reveal that besides the physiologically highly relevant magnesium, metals like zinc and manganese can assist in phosphoryl transfer, however, only a few support efficient substrate turnover (turnover catalysis). Those trace metals allow for substrate binding and phosphotransfer but hamper product release. We further established the unique role of magnesium as the common biologically relevant divalent metal ideally enabling (co-) substrate binding and orientation. Magnesium allows stable substrate binding and, on the other hand accelerates product release, thus being extremely efficient in turnover catalysis. We extended our studies to non-catalytic functions of protein kinases looking at pseudokinases, a subfamily of protein kinases inherently lacking critical residues for catalysis. Recently, pseudokinases have been linked to human diseases. Some pseudokinases are still capable of binding metal ions, yet have lost the ability to transfer the phosphoryl group from ATP to a given substrate. Here metal ions stabilize an active like, though catalytically unproductive conformation and are therefore crucial to maintain non-catalytic function. Finally, we demonstrate for the canonical kinase PKA that the trace metal manganese alone can stabilize protein kinases in an active like conformation allowing them to bind substrates even in the absence of nucleotides.

The Molecular Basis for Specificity at the Level of the Protein Kinase a Catalytic Subunit

Assembly of multi enzyme complexes at subcellular localizations by anchoring- and scaffolding proteins represents a pivotal mechanism for achieving spatiotemporal regulation of cellular signaling after hormone receptor targeting [for review, see (1)]. In the 3' 5'-cyclic adenosine monophosphate (cAMP) dependent protein kinase (PKA) signaling pathway it is generally accepted that specificity is secured at several levels. This includes at the first level stimulation of receptors coupled to heterotrimeric G proteins which through stimulation of adenylyl cyclase (AC) forms the second messenger cAMP. Cyclic AMP has several receptors including PKA. PKA is a tetrameric holoenzyme consisting of a regulatory (R) subunit dimer and two catalytic (C) subunits. The R subunit is the receptor for cAMP and compartmentalizes cAMP signals through binding to cell and tissue-specifically expressed A kinase anchoring proteins (AKAPs). The current dogma tells that in the presence of cAMP, PKA dissociates into an R subunit dimer and two C subunits which are free to phosphorylate relevant substrates in the cytosol and nucleus. The release of the C subunit has raised the question how specificity of the cAMP and PKA signaling pathway is maintained when the C subunit no longer is attached to the R subunit-AKAP complex. An increasing body of evidence points toward a regulatory role of the cAMP and PKA signaling pathway by targeting the C subunits to various C subunit binding proteins in the cytosol and nucleus. Moreover, recent identification of isoform specific amino acid sequences, motifs and three dimensional structures have together provided new insight into how PKA at the level of the C subunit may act in a highly isoform-specific fashion. Here we discuss recent understanding of specificity of the cAMP and PKA signaling pathway based on C subunit subcellular targeting as well as evolution of the C subunit structure that may contribute to the dynamic regulation of C subunit activity.

Understanding how cAMP-dependent protein kinase can catalyze phosphoryl transfer in the presence of Ca2+and Sr2+: a QM/MM study

Theoretical results demonstrate for the first time at the molecular level that the overall PKAc-catalyzed phosphoryl-transfer reaction is plausible with Ca²⁺and Sr²⁺, alkaline earth metal ions other than Mg²⁺, which is in good agreement with experiments.

Protein Kinase A Catalytic Subunit Primed for Action: Time-Lapse Crystallography of Michaelis Complex Formation

The catalytic subunit of the cyclic AMP-dependent protein kinase A (PKAc) catalyzes the transfer of the γ-phosphate of bound Mg2ATP to a serine or threonine residue of a protein substrate. Here, time-lapse X-ray crystallography was used to capture a series of complexes of PKAc with an oligopeptide substrate and unreacted Mg2ATP, including the Michaelis complex, that reveal important geometric rearrangements in and near the active site preceding the phosphoryl transfer reaction. Contrary to the prevailing view, Mg(2+) binds first to the M1 site as a complex with ATP and is followed by Mg(2+) binding to the M2 site. Concurrently, the target serine hydroxyl of the peptide substrate rotates away from the active site toward the bulk solvent, which breaks the hydrogen bond with D166. Lastly, the serine hydroxyl of the substrate rotates back toward D166 to form the Michaelis complex with the active site primed for phosphoryl transfer.

Divalent Metal Ions Mg²⁺ and Ca²⁺ Have Distinct Effects on Protein Kinase A Activity and Regulation

cAMP-dependent protein kinase (PKA) is regulated primarily in response to physiological signals while nucleotides and metals may provide fine-tuning. PKA can use different metal ions for phosphoryl transfer, yet some, like Ca(2+), do not support steady-state catalysis. Fluorescence Polarization (FP) and Surface Plasmon Resonance (SPR) were used to study inhibitor and substrate interactions with PKA. The data illustrate how metals can act differentially as a result of their inherent coordination properties. We found that Ca(2+), in contrast to Mg(2+), does not induce high-affinity binding of PKA to pseudosubstrate inhibitors. However, Ca(2+) works in a single turnover mode to allow for phosphoryl-transfer. Using a novel SPR approach, we were able to directly monitor the interaction of PKA with a substrate in the presence of Mg(2+)ATP. This allows us to depict the entire kinase reaction including complex formation as well as release of the phosphorylated substrate. In contrast to Mg(2+), Ca(2+) apparently slows down the enzymatic reaction. A focus on individual reaction steps revealed that Ca(2+) is not as efficient as Mg(2+) in stabilizing the enzyme:substrate complex. The opposite holds true for product dissociation where Mg(2+) easily releases the phospho-substrate while Ca(2+) traps both reaction products at the active site. This explains the low steady-state activity in the presence of Ca(2+). Furthermore, Ca(2+) is able to modulate kinase activity as well as inhibitor binding even in the presence of Mg(2+). We therefore hypothesize that the physiological metal ions Mg(2+) and Ca(2+) both play a role in kinase activity and regulation. Since PKA is localized close to calcium channels and may render PKA activity susceptible to Ca(2+), our data provide a possible mechanism for novel crosstalk between cAMP and calcium signaling.

Bridging scales through multiscale modeling: a case study on protein kinase A

The goal of multiscale modeling in biology is to use structurally based physico-chemical models to integrate across temporal and spatial scales of biology and thereby improve mechanistic understanding of, for example, how a single mutation can alter organism-scale phenotypes. This approach may also inform therapeutic strategies or identify candidate drug targets that might otherwise have been overlooked. However, in many cases, it remains unclear how best to synthesize information obtained from various scales and analysis approaches, such as atomistic molecular models, Markov state models (MSM), subcellular network models, and whole cell models. In this paper, we use protein kinase A (PKA) activation as a case study to explore how computational methods that model different physical scales can complement each other and integrate into an improved multiscale representation of the biological mechanisms. Using measured crystal structures, we show how molecular dynamics (MD) simulations coupled with atomic-scale MSMs can provide conformations for Brownian dynamics (BD) simulations to feed transitional states and kinetic parameters into protein-scale MSMs. We discuss how milestoning can give reaction probabilities and forward-rate constants of cAMP association events by seamlessly integrating MD and BD simulation scales. These rate constants coupled with MSMs provide a robust representation of the free energy landscape, enabling access to kinetic, and thermodynamic parameters unavailable from current experimental data. These approaches have helped to illuminate the cooperative nature of PKA activation in response to distinct cAMP binding events. Collectively, this approach exemplifies a general strategy for multiscale model development that is applicable to a wide range of biological problems.

Phosphoryl Transfer Reaction Snapshots in Crystals: INSIGHTS INTO THE MECHANISM OF PROTEIN KINASE A CATALYTIC SUBUNIT

To study the catalytic mechanism of phosphorylation catalyzed by cAMP-dependent protein kinase (PKA) a structure of the enzyme-substrate complex representing the Michaelis complex is of specific interest as it can shed light on the structure of the transition state. However, all previous crystal structures of the Michaelis complex mimics of the PKA catalytic subunit (PKAc) were obtained with either peptide inhibitors or ATP analogs. Here we utilized Ca(2+) ions and sulfur in place of the nucleophilic oxygen in a 20-residue pseudo-substrate peptide (CP20) and ATP to produce a close mimic of the Michaelis complex. In the ternary reactant complex, the thiol group of Cys-21 of the peptide is facing Asp-166 and the sulfur atom is positioned for an in-line phosphoryl transfer. Replacement of Ca(2+) cations with Mg(2+) ions resulted in a complex with trapped products of ATP hydrolysis: phosphate ion and ADP. The present structural results in combination with the previously reported structures of the transition state mimic and phosphorylated product complexes complete the snapshots of the phosphoryl transfer reaction by PKAc, providing us with the most thorough picture of the catalytic mechanism to date.

The energy landscape of adenylate kinase during catalysis

Kinases perform phosphoryl-transfer reactions in milliseconds; without enzymes, these reactions would take about 8,000 years under physiological conditions. Despite extensive studies, a comprehensive understanding of kinase energy landscapes, including both chemical and conformational steps, is lacking. Here we scrutinize the microscopic steps in the catalytic cycle of adenylate kinase, through a combination of NMR measurements during catalysis, pre-steady-state kinetics, molecular-dynamics simulations and crystallography of active complexes. We find that the Mg(2+) cofactor activates two distinct molecular events: phosphoryl transfer (>10(5)-fold) and lid opening (10(3)-fold). In contrast, mutation of an essential active site arginine decelerates phosphoryl transfer 10(3)-fold without substantially affecting lid opening. Our results highlight the importance of the entire energy landscape in catalysis and suggest that adenylate kinases have evolved to activate key processes simultaneously by precise placement of a single, charged and very abundant cofactor in a preorganized active site.

A QM/MM study of Kemptide phosphorylation catalyzed by protein kinase A. The role of Asp166 as a general acid/base catalyst

In this work a theoretical study of the γ-phosphoryl group transfer from ATP to Ser17 of the synthetic substrate Kemptide (LRRASLG) in protein kinase A (PKA) has been carried out with a solvated model of the PKA-Mg2ATP-Kemptide system based on the X-ray crystallographic structure. We have used high levels (B3LYP/MM and MP2/MM) of theory to determine the overall reaction paths of the so-called concerted loose mechanism trying to clarify some aspects of that mechanism still under debate. Our calculations demonstrate for the first time in a complete model of the ternary system the viability of the final step of the catalytic mechanism in which the protonation of the phosphokemptide product by Asp166 takes place. Asp166 is a base catalyst that abstracts the HγSer17 of Kemptide thus facilitating the phosphoryl transfer, but it also acts as an acid catalyst by donating the proton just accepted from Ser17 to the O2γATP atom of the phosphoryl group.

Dynamic architecture of a protein kinase

Protein kinases are dynamically regulated signaling proteins that act as switches in the cell by phosphorylating target proteins. To establish a framework for analyzing linkages between structure, function, dynamics, and allostery in protein kinases, we carried out multiple microsecond-scale molecular-dynamics simulations of protein kinase A (PKA), an exemplar active kinase. We identified residue-residue correlated motions based on the concept of mutual information and used the Girvan-Newman method to partition PKA into structurally contiguous "communities." Most of these communities included 40-60 residues and were associated with a particular protein kinase function or a regulatory mechanism, and well-known motifs based on sequence and secondary structure were often split into different communities. The observed community maps were sensitive to the presence of different ligands and provide a new framework for interpreting long-distance allosteric coupling. Communication between different communities was also in agreement with the previously defined architecture of the protein kinase core based on the "hydrophobic spine" network. This finding gives us confidence in suggesting that community analyses can be used for other protein kinases and will provide an efficient tool for structural biologists. The communities also allow us to think about allosteric consequences of mutations that are linked to disease.

A QM/MM study of the associative mechanism for the phosphorylation reaction catalyzed by protein kinase A and its D166A mutant

Here we analyze in detail the possible catalytic role of the associative mechanism in the γ-phosphoryl transfer reaction in the catalytic subunit of the mammalian cyclic AMP-dependent protein kinase (PKA) enzyme and its D166A mutant. We have used a complete solvated model of the ATP-Mg2-Kemptide/PKA system and good levels of theory (B3LYP/MM and MP2/MM) to determine several potential energy paths from different MD snapshots, and we present a deep analysis of the interaction distances and energies between ligands, metals and enzyme residues. We have also tested the electrostatic stabilization of the transition state structures localized herein with the charge balance hypothesis. Overall, the results obtained in this work reopen the discussion about the plausibility of the associative reaction pathway and highlight the proposed role of the catalytic triad Asp166-Lys168-Thr201.

Molecular features of product release for the PKA catalytic cycle

Although ADP release is the rate limiting step in product turnover by protein kinase A, the steps and motions involved in this process are not well resolved. Here we report the apo and ADP bound structures of the myristylated catalytic subunit of PKA at 2.9 and 3.5 Å resolution, respectively. The ADP bound structure adopts a conformation that does not conform to the previously characterized open, closed, or intermediate states. In the ADP bound structure, the C-terminal tail and Gly-rich loop are more closed than in the open state adopted in the apo structure but are also much more open than the intermediate or closed conformations. Furthermore, ADP binds at the active site with only one magnesium ion, termed Mg2 from previous structures. These structures thus support a model where ADP release proceeds through release of the substrate and Mg1 followed by lifting of the Gly-rich loop and disengagement of the C-terminal tail. Coupling of these two structural elements with the release of the first metal ion fills in a key step in the catalytic cycle that has been missing and supports an ensemble of correlated conformational states that mediate the full catalytic cycle for a protein kinase.

Theoretical studies of cyclic adenosine monophosphate dependent protein kinase: native enzyme and ground-state and transition-state analogues

The mechanisms of phosphoryl transfer enzymes have garnered considerable attention. Cyclic AMP-dependent protein kinase (cAPK) catalyzes the transfer of the γ phosphoryl group of ATP to the serine hydroxyl group of a peptide chain. Metal-containing fluoro species have been used as transition-state and ground-state analogues in a variety of phosphoryl transfer enzymes and have shed light on the nature of the requirements in the active site to catalyze phosphoryl transfer. For cAPK, we present computational studies of the mechanism of phosphoryl transfer and the structure and (19)F NMR spectra of various ground- (BeF3(-)) and transition-state (MgF3(-), AlF4(-), and AlF3(0)) analogues. With native substrate, the phosphoryl transfer proceeds through a five-coordinate phosphorane transition state, i.e., there is not a five-coordinate phosphorane intermediate. Comparisons of simulated and experimental (19)F NMR spectra show cAPK prefers a monoanionic analogue (MgF3(-) or AlF4(-)) over a neutral analogue (AlF3), supporting the charge balance hypothesis.

Unique kinase catalytic mechanism of AceK with a single magnesium ion

Isocitrate dehydrogenase kinase/phosphatase (AceK) is the founding member of the protein phosphorylation system in prokaryotes. Based on the novel and unique structural characteristics of AceK recently uncovered, we sought to understand its kinase reaction mechanism, along with other features involved in the phosphotransfer process. Herein we report density functional theory QM calculations of the mechanism of the phosphotransfer reaction catalysed by AceK. The transition states located by the QM calculations indicate that the phosphorylation reaction, catalysed by AceK, follows a dissociative mechanism with Asp457 serving as the catalytic base to accept the proton delivered by the substrate. Our results also revealed that AceK prefers a single Mg²⁺-containing active site in the phosphotransfer reaction. The catalytic roles of conserved residues in the active site are discussed.

Evolutionary paths of the cAMP-dependent protein kinase (PKA) catalytic subunits

3',5'-cyclic adenosine monophosphate (cAMP) dependent protein kinase or protein kinase A (PKA) has served as a prototype for the large family of protein kinases that are crucially important for signal transduction in eukaryotic cells. The PKA catalytic subunits Cα and Cβ, encoded by the two genes PRKACA and PRKACB, respectively, are among the best understood and characterized human kinases. Here we have studied the evolution of this gene family in chordates, arthropods, mollusks and other animals employing probabilistic methods and show that Cα and Cβ arose by duplication of an ancestral PKA catalytic subunit in a common ancestor of vertebrates. The two genes have subsequently been duplicated in teleost fishes. The evolution of the PRKACG retroposon in simians was also investigated. Although the degree of sequence conservation in the PKA Cα/Cβ kinase family is exceptionally high, a small set of signature residues defining Cα and Cβ subfamilies were identified. These conserved residues might be important for functions that are unique to the Cα or Cβ clades. This study also provides a good example of a seemingly simple phylogenetic problem which, due to a very high degree of sequence conservation and corresponding weak phylogenetic signals, combined with problematic nonphylogenetic signals, is nontrivial for state-of-the-art probabilistic phylogenetic methods.

Phosphoryl transfer by protein kinase A is captured in a crystal lattice

The catalytic (C) subunit of cAMP-dependent protein kinase (PKA) is a serine/threonine kinase responsible for most of the effects of cAMP signaling, and PKA serves as a prototype for the entire kinase family. Despite multiple studies of PKA, the steps involved in phosphoryl transfer, the roles of the catalytically essential magnesium ions, and the processes that govern the rate-limiting step of ADP release are unresolved. Here we identified conditions that yielded slow phosphoryl transfer of the γ-phosphate from the generally nonhydrolyzable analog of ATP, adenosine-5'-(β,γ-imido)triphosphate (AMP-PNP), onto a substrate peptide within protein crystals. By trapping both products in the crystal lattice, we now have a complete resolution profile of all the catalytic steps. One crystal structure refined to 1.55 Å resolution shows two states of the protein with 55% displaying intact AMP-PNP and an unphosphorylated substrate and 45% displaying transfer of the γ-phosphate of AMP-PNP onto the substrate peptide yielding AMP-PN and a phosphorylated substrate. Another structure refined to 2.15 Å resolution displays complete phosphoryl transfer to the substrate. These structures, in addition to trapping both products in the crystal lattice, implicate one magnesium ion, previously termed Mg2, as the more stably bound ion. Following phosphoryl transfer, Mg2 recruits a water molecule to retain an octahedral coordination geometry suggesting the strong binding character of this magnesium ion, and Mg2 remains in the active site following complete phosphoryl transfer while Mg1 is expelled. Loss of Mg1 may thus be an important part of the rate-limiting step of ADP release.

Simulating GTP:Mg and GDP:Mg with a simple force field: a structural and thermodynamic analysis

Di- and tri-phosphate nucleotides are essential cofactors for many proteins, usually in an Mg(2+) -bound form. Proteins like GTPases often detect the difference between NDP and NTP and respond by changing conformations. To study such complexes, simple, fixed charge force fields have been used, which allow long simulations and precise free energy calculations. The preference for NTP or NDP binding depends on many factors, including ligand structure and Mg(2+) coordination and the changes they undergo upon binding. Here, we use a simple force field to examine two Mg(2+) coordination modes for the unbound GDP and GTP: direct, or "Inner Sphere" (IS) coordination by one or more phosphate oxygens and indirect, "Outer Sphere" (OS) coordination involving one or more bridging waters. We compare GTP: and GDP:Mg binding with OS and IS coordination; combining the results with experimental data then indicates that GTP prefers the latter. We also examine different kinds of IS coordination and their sensitivity to a key force field parameter: the optimal Mg:oxygen van der Waals distance Rmin . Increasing Rmin improves the Mg:oxygen distances, the GTP: and GDP:Mg binding affinities, and the fraction of GTP:Mg with β + γ phosphate coordination, but does not improve or change the GTP/GDP affinity difference, which remains much larger than experiment. It has no effect on the free energy of GDP binding to a GTPase.

Price to be paid for two-metal catalysis: magnesium ions that accelerate chemistry unavoidably limit product release from a protein kinase

Incorporation of divalent metal ions into an active site is a fundamental catalytic tool used by diverse enzymes. Divalent cations are used by protein kinases to both stabilize ATP binding and accelerate chemistry. Kinetic analysis establishes that Cyclin-dependent kinase 2 (CDK2) requires simultaneous binding of two Mg(2+) ions for catalysis of phosphoryl transfer. This tool, however, comes with a price: the rate-acceleration effects are opposed by an unavoidable rate-limiting consequence of the use of two Mg(2+) ions by CDK2. The essential metal ions stabilize ADP product binding and limit the overall rate of the reaction. We demonstrate that product release is rate limiting for activated CDK2 and evaluate the effects of the two catalytically essential Mg(2+) ions on the stability of the ADP product within the active site. We present two new crystal structures of CDK2 bound to ADP showing how the phosphate groups can be coordinated by either one or two Mg(2+) ions, with the occupancy of one site in a weaker equilibrium. Molecular dynamics simulations indicate that ADP phosphate mobility is more restricted when ADP is coordinated by two Mg(2+) ions compared to one. The structural similarity between the rigid ADP·2Mg product and the cooperatively assembled transition state provides a mechanistic rational for the rate-limiting ADP release that is observed. We demonstrate that although the simultaneous binding of two Mg(2+) ions is essential for efficient phosphoryl transfer, the presence of both Mg(2+) ions in the active site also cooperatively increases ADP affinity and opposes its release. Evolution of protein kinases must have involved careful tuning of the affinity for the second Mg(2+) ion in order to balance the needs to stabilize the chemical transition state and allow timely product release. The link between Mg(2+) site affinity and activity presents a chemical handle that may be used by regulatory factors as well as explain some mutational effects.

An ATP-site on-off switch that restricts phosphatase accessibility of Akt

The protein serine-threonine kinase Akt undergoes a substantial conformational change upon activation, which is induced by the phosphorylation of two critical regulatory residues, threonine 308 and serine 473. Paradoxically, treating cells with adenosine 5'-triphosphate (ATP)-competitive inhibitors of Akt results in increased phosphorylation of both residues. We show that binding of ATP-competitive inhibitors stabilized a conformation in which both phosphorylated sites were inaccessible to phosphatases. ATP binding also produced this protection of the phosphorylated sites, whereas interaction with its hydrolysis product adenosine 5'-diphosphate (ADP) or allosteric Akt inhibitors resulted in increased accessibility of these phosphorylated residues. ATP-competitive inhibitors mimicked ATP by targeting active Akt. Forms of Akt activated by an oncogenic mutation or myristoylation were more potently inhibited by the ATP-competitive inhibitors than was wild-type Akt. These data support a new model of kinase regulation, wherein nucleotides modulate an on-off switch in Akt through conformational changes, which is disrupted by ATP-competitive inhibitors.

Identification and characterization of novel mutations in the human gene encoding the catalytic subunit Calpha of protein kinase A (PKA)

The genes PRKACA and PRKACB encode the principal catalytic (C) subunits of protein kinase A (PKA) Cα and Cβ, respectively. Cα is expressed in all eukaryotic tissues examined and studies of Cα knockout mice demonstrate a crucial role for Cα in normal physiology. We have sequenced exon 2 through 10 of PRKACA from the genome of 498 Norwegian donors and extracted information about PRKACA mutations from public databases. We identified four interesting nonsynonymous point mutations, Arg45Gln, Ser109Pro, Gly186Val, and Ser263Cys, in the Cα1 splice variant of the kinase. Cα variants harboring the different amino acid mutations were analyzed for kinase activity and regulatory (R) subunit binding. Whereas mutation of residues 45 and 263 did not alter catalytic activity or R subunit binding, mutation of Ser(109) significantly reduced kinase activity while R subunit binding was unaltered. Mutation of Cα Gly(186) completely abrogated kinase activity and PKA type I but not type II holoenzyme formation. Gly(186) is located in the highly conserved DFG motif of Cα and mutation of this residue to Val was predicted to result in loss of binding of ATP and Mg(2+), which may explain the kinetic inactivity. We hypothesize that individuals born with mutations of Ser(109) or Gly(186) may be faced with abnormal development and possibly severe disease.

Multi-Scale Continuum Modeling of Biological Processes: From Molecular Electro-Diffusion to Sub-Cellular Signaling Transduction

This article provides a brief review of multi-scale modeling at the molecular to cellular scale, with new results for heart muscle cells. A finite element-based simulation package (SMOL) was used to investigate the signaling transduction at molecular and sub-cellular scales (http://mccammon.ucsd.edu/smol/, http://FETK.org) by numerical solution of time-dependent Smoluchowski equations and a reaction-diffusion system. At the molecular scale, SMOL has yielded experimentally-validated estimates of the diffusion-limited association rates for the binding of acetylcholine to mouse acetylcholinesterase using crystallographic structural data. The predicted rate constants exhibit increasingly delayed steady-state times with increasing ionic strength and demonstrate the role of an enzyme's electrostatic potential in influencing ligand binding. At the sub-cellular scale, an extension of SMOL solves a non-linear, reaction-diffusion system describing Ca²⁺ ligand buffering and diffusion in experimentally-derived rodent ventricular myocyte geometries. Results reveal the important role for mobile and stationary Ca²⁺ buffers, including Ca²⁺ indicator dye. We found that the alterations in Ca²⁺-binding and dissociation rates of troponin C (TnC) and total TnC concentration modulate subcellular Ca²⁺ signals. Model predicts that reduced off-rate in whole troponin complex (TnC, TnI, TnT) versus reconstructed thin filaments (Tn, Tm, actin) alters cytosolic Ca²⁺ dynamics under control conditions or in disease-linked TnC mutations. The ultimate goal of these studies is to develop scalable methods and theories for integration of molecular-scale information into simulations of cellular-scale systems.

Conformational selection through electrostatics: Free energy simulations of GTP and GDP binding to archaeal initiation factor 2

Archaeal Initiation Factor 2 is a GTPase involved in protein biosynthesis. In its GTP-bound, "ON" conformation, it binds an initiator tRNA and carries it to the ribosome. In its GDP-bound, "OFF" conformation, it dissociates from tRNA. To understand the specific binding of GTP and GDP and their dependence on the conformational state, molecular dynamics free energy simulations were performed. The ON state specificity was predicted to be weak, with a GTP/GDP binding free energy difference of -1 kcal/mol, favoring GTP. The OFF state specificity is larger, 4 kcal/mol, favoring GDP. The overall effects result from a competition among many interactions in several complexes. To interpret them, we use a simpler, dielectric continuum model. Several effects are robust with respect to the model details. Both nucleotides have a net negative charge, so that removing them from solvent into the binding pocket carries a desolvation penalty, which is large for the ON state, and strongly disfavors GTP binding compared to GDP. Short-range interactions between the additional GTP phosphate group and ionized sidechains in the binding pocket offset most, but not all of the desolvation penalty; more distant groups also contribute significantly, and the switch 1 loop only slightly. The desolvation penalty is lower for the more open, wetter OFF state, and the GTP/GDP difference much smaller. Short-range interactions in the binding pocket and with more distant groups again make a significant contribution. Overall, the simulations help explain how conformational selection is achieved with a single phosphate group.

Autoregulation of kinase dephosphorylation by ATP binding in AGC protein kinases

AGC kinases, including the three Akt (protein kinase B) isoforms, protein kinase A (PKA) and all protein kinase C (PKC) isoforms, require activation loop phosphorylation (threonine 308 in Akt1) as well as phosphorylation of a C-terminal residue (serine 473 in Akt1) for catalytic activity and phosphorylation of downstream targets. Conversely, phosphatases reverse these phosphorylations. Virtually all cellular processes are affected by AGC kinases, a circumstance that has led to intense scrutiny of the molecular mechanisms that regulate phosphorylation of these kinases. Here, we review a new layer of control of phosphorylation in Akt, PKA and PKC pointing to ATP binding pocket occupancy as a means to decelerate dephosphorylation of these and, potentially, other kinases. This additional level of kinase regulation opens the door to search for new functional motifs for the rational design of non- ATP-competitive kinase inhibitors that discriminate within and between protein kinase families.

Conformational selection by the aIF2 GTPase: a molecular dynamics study of functional pathways

Archaeal initiation factor 2 (aIF2) is a GTPase involved in protein biosynthesis. In its GTP-bound, "ON" conformation, it binds an initiator tRNA and carries it to the ribosome. In its GDP-bound, "OFF" conformation, it dissociates from tRNA. To improve our understanding of the role of each conformational state in the aIF2 "life cycle", we start from the state immediately after GTP hydrolysis, ON:GDP:P(i) (where P(i) is inorganic phosphate), and consider the possible next steps on the pathway to the OFF:GDP product. The first possibility is P(i) dissociation, leading to ON:GDP, which could then relax into OFF:GDP. We use molecular dynamics simulations to compute the P(i) dissociation free energy and show that dissociation is highly favorable. The second possibility is conformational relaxation into the OFF state before P(i) dissociation, to form OFF:GDP:P(i). We estimate the corresponding free energy approximately, 2 ± 3.5 kcal/mol, so that this is an uphill or weakly downhill process. A third possibility is relaxation into another conformation, neither ON nor OFF. Indeed, a third, "MIXED" conformation was seen recently in a crystal structure of the aIF2:GDP:P(i) complex. For this conformational state, P(i) dissociation is weakly unfavorable, in contrast to the ON and OFF states. From this, we will deduce that if the MIXED:GDP complex is not too unstable, the ON:GDP:P(i) → MIXED:GDP:P(i) transformation is a downhill process, which can occur spontaneously. This suggests that the MIXED state could be a functional intermediate.

Resistance of Akt kinases to dephosphorylation through ATP-dependent conformational plasticity

Phosphorylation of a threonine residue (T308 in Akt1) in the activation loop of Akt kinases is a prerequisite for deregulated Akt activity frequently observed in neoplasia. Akt phosphorylation in vivo is balanced by the opposite activities of kinases and phosphatases. Here we describe that targeting Akt kinase to the cell membrane markedly reduced sensitivity of phosphorylated Akt to dephosphorylation by protein phosphatase 2A. This effect was amplified by occupancy of the ATP binding pocket by either ATP or ATP-competitive inhibitors. Mutational analysis revealed that R273 in Akt1 and the corresponding R274 in Akt2 are essential for shielding T308 in the activation loop against dephosphorylation. Thus, occupancy of the nucleotide binding pocket of Akt kinases enables intramolecular interactions that restrict phosphatase access and sustain Akt phosphorylation. This mechanism provides an explanation for the "paradoxical" Akt hyperphosphorylation induced by ATP-competitive inhibitor, A-443654. The lack of phosphatase resistance further contributes insight into the mechanism by which the human Akt2 R274H missense mutation may cause autosomal-dominant diabetes mellitus.

Briefly bound to activate: transient binding of a second catalytic magnesium activates the structure and dynamics of CDK2 kinase for catalysis

We have determined high-resolution crystal structures of a CDK2/Cyclin A transition state complex bound to ADP, substrate peptide, and MgF(3)(-). Compared to previous structures of active CDK2, the catalytic subunit of the kinase adopts a more closed conformation around the active site and now allows observation of a second Mg(2+) ion in the active site. Coupled with a strong [Mg(2+)] effect on in vitro kinase activity, the structures suggest that the transient binding of the second Mg(2+) ion is necessary to achieve maximum rate enhancement of the chemical reaction, and Mg(2+) concentration could represent an important regulator of CDK2 activity in vivo. Molecular dynamics simulations illustrate how the simultaneous binding of substrate peptide, ATP, and two Mg(2+) ions is able to induce a more rigid and closed organization of the active site that functions to orient the phosphates, stabilize the buildup of negative charge, and shield the subsequently activated γ-phosphate from solvent.

A conserved non-canonical motif in the pseudoactive site of the ROP5 pseudokinase domain mediates its effect on Toxoplasma virulence

The ROP5 family is a closely related set of polymorphic pseudokinases that are critical to the ability of Toxoplasma to cause disease. Polymorphisms in ROP5 also make it a major determinant of strain-specific differences in virulence. ROP5 possesses all of the major kinase motifs required for catalysis except for a substitution at the catalytic Asp. We show that this substitution in the catalytic loop of ROP5 is part of a motif conserved in other pseudokinases of both Toxoplasma and human origin, and that this motif is required for the full activity in vivo of ROP5. This suggests evolutionary selection at this site for a biochemical function, rather than simple drift away from catalysis. We present the crystal structures of a virulent isoform of ROP5 both in its ATP-bound and -unbound states and have demonstrated that despite maintaining the canonical ATP-binding motifs, ROP5 binds ATP in a distorted conformation mediated by unusual magnesium coordination sites that would not be predicted from the primary sequence. In addition, we have mapped the polymorphisms spread throughout the primary sequence of ROP5 to two major surfaces, including the activation segment of ROP5. This suggests that the pseudoactive site of this class of pseudokinases may have evolved to use the canonical ATP-binding motifs for non-catalytic signaling through allostery.

Free energy simulations of a GTPase: GTP and GDP binding to archaeal initiation factor 2

Archaeal initiation factor 2 (aIF2) is a protein involved in the initiation of protein biosynthesis. In its GTP-bound, "ON" conformation, aIF2 binds an initiator tRNA and carries it to the ribosome. In its GDP-bound, "OFF" conformation, it dissociates from tRNA. To understand the specific binding of GTP and GDP and its dependence on the ON or OFF conformational state of aIF2, molecular dynamics free energy simulations (MDFE) are a tool of choice. However, the validity of the computed free energies depends on the simulation model, including the force field and the boundary conditions, and on the extent of conformational sampling in the simulations. aIF2 and other GTPases present specific difficulties; in particular, the nucleotide ligand coordinates a divalent Mg(2+) ion, which can polarize the electronic distribution of its environment. Thus, a force field with an explicit treatment of electronic polarizability could be necessary, rather than a simpler, fixed charge force field. Here, we begin by comparing a fixed charge force field to quantum chemical calculations and experiment for Mg(2+):phosphate binding in solution, with the force field giving large errors. Next, we consider GTP and GDP bound to aIF2 and we compare two fixed charge force fields to the recent, polarizable, AMOEBA force field, extended here in a simple, approximate manner to include GTP. We focus on a quantity that approximates the free energy to change GTP into GDP. Despite the errors seen for Mg(2+):phosphate binding in solution, we observe a substantial cancellation of errors when we compare the free energy change in the protein to that in solution, or when we compare the protein ON and OFF states. Finally, we have used the fixed charge force field to perform MDFE simulations and alchemically transform GTP into GDP in the protein and in solution. With a total of about 200 ns of molecular dynamics, we obtain good convergence and a reasonable statistical uncertainty, comparable to the force field uncertainty, and somewhat lower than the predicted GTP/GDP binding free energy differences. The sign and magnitudes of the differences can thus be interpreted at a semiquantitative level, and are found to be consistent with the experimental binding preferences of ON- and OFF-aIF2.

Dynamically committed, uncommitted, and quenched states encoded in protein kinase A revealed by NMR spectroscopy

Protein kinase A (PKA) is a ubiquitous phosphoryl transferase that mediates hundreds of cell signaling events. During turnover, its catalytic subunit (PKA-C) interconverts between three major conformational states (open, intermediate, and closed) that are dynamically and allosterically activated by nucleotide binding. We show that the structural transitions between these conformational states are minimal and allosteric dynamics encode the motions from one state to the next. NMR and molecular dynamics simulations define the energy landscape of PKA-C, with the substrate allowing the enzyme to adopt a broad distribution of conformations (dynamically committed state) and the inhibitors (high magnesium and pseudosubstrate) locking it into discrete minima (dynamically quenched state), thereby reducing the motions that allow turnover. These results unveil the role of internal dynamics in both kinase function and regulation.