Dynamics connect substrate recognition to catalysis in protein kinase A

Unbiased clustering of residues undergoing synchronous motions in proteins using NMR spin relaxation data

Biological macromolecules are dynamic entities that transition between various conformational states, often playing a vital role in biological functions. Their inherent flexibility spans a broad range of timescales. Motions occurring within the microsecond to millisecond range are especially important, as they are integral to processes such as enzyme catalysis, folding, ligand binding, and allostery. NMR Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion measurements are the preferred method for characterizing macromolecular motions at atomic resolution. However, it is still uncertain whether the functional motions of multiple residues in macromolecules need to be coordinated and/or synchronized within the protein matrix in order to perform the desired function. Here, we illustrate an unbiased method to analyze NMR relaxation dispersion and identify dynamic clusters of residues that fluctuate on similar timescales within proteins. The method requires relaxation dispersion data for backbone amides or side-chain methyl groups, which are globally fitted using the Bloch-McConnell equations for each pair of residues. The goodness of the pairwise fitting serves as a metric to construct two-dimensional synchronous dynamics (SyncDyn) maps, allowing us to identify residue clusters whose dynamics are influenced by ligand binding. We applied our method to the catalytic subunit of the cAMP-dependent protein kinase A (PKAC) and the T17A mutant of ribonuclease A (RNAse A). The SyncDyn maps for PKA-C showed distinct clusters of residues located in critical allosteric sites. Nucleotide binding activates the movement of residues at the interface between the two lobes and also those distal to the active site. In the case of RNAse A, the SyncDyn maps show that residues fluctuating with the same time scale are interspersed in both lobes of the enzyme. Overall, our approach eliminates arbitrary manual selection of residues for dynamic clustering and objectively identifies all possible residue pairs that fluctuate synchronously, i.e. on the same timescale.

Structural insights of AKT and its activation mechanism for drug development

AKT1, a serine/threonine kinase, is pivotal in signaling and regulating cell survival, proliferation, and metabolism. This review focuses on the structural insights and the essential features required for its active conformation. AKT belongs to the AGC kinase group and has three isoforms: AKT1, AKT2, and AKT3. AKT has three functional regions: PH domain, kinase domain, and hydrophobic motif. AKT1 activation involves intricate conformational changes, including transitions in the αC-in, DFG-in, G-loop, activation loop, and PH domain out, S-spine and R-spine formation, as well as phosphorylation at Thr 308 and Ser 473, which enable AKT1 to adopt active conformation. The analysis highlights the limitations of the AlphaFold-predicted AKT1 structure, which lacks key elements of the active state, including ATP, magnesium ion coordination, phosphatidylinositol-(1,3,4,5)-tetraphosphate, substrate peptide, and phosphorylation at Thr 308 and Ser 473. This study underscores the necessity of these features for stabilizing the kinase domain and facilitating efficient substrate phosphorylation. By consolidating structural insights and activation mechanisms, this review aims to inform the development of computational models and targeted therapeutics for AKT1 activators in diseases such as hepatic ischemia-reperfusion injury, cerebral ischemia, acute hepatic failure, subarachnoid hemorrhage, and alzheimer's disease.

Allostery

Allostery describes the ability of biological macromolecules to transmit signals spatially through the molecule from an allosteric site – a site that is distinct from orthosteric binding sites of primary, endogenous ligands – to the functional or active site. This review starts with a historical overview and a description of the classical example of allostery – hemoglobin – and other well-known examples (aspartate transcarbamoylase, Lac repressor, kinases, G-protein-coupled receptors, adenosine triphosphate synthase, and chaperonin). We then discuss fringe examples of allostery, including intrinsically disordered proteins and inter-enzyme allostery, and the influence of dynamics, entropy, and conformational ensembles and landscapes on allosteric mechanisms, to capture the essence of the field. Thereafter, we give an overview over central methods for investigating molecular mechanisms, covering experimental techniques as well as simulations and artificial intelligence (AI)-based methods. We conclude with a review of allostery-based drug discovery, with its challenges and opportunities: with the recent advent of AI-based methods, allosteric compounds are set to revolutionize drug discovery and medical treatments.

Crystallographic, kinetic, and calorimetric investigation of PKA interactions with L-type calcium channels and Rad GTPase

β-Adrenergic signaling activates cAMP-dependent PKA, which regulates the activity of L-type voltage-gated calcium channels such as CaV1.2. Several PKA target sites in the C-terminal tail of CaV1.2 have been identified, and their phosphorylation has been suggested to increase currents in specific tissues or heterologous expression systems. However, augmentation of CaV1.2 currents in the heart is instead mediated by phosphorylation of Rad, a small GTPase that can inhibit CaV1.2. It is unclear how each of the proposed target sites in CaV1.2 and Rad rank toward their recognition by PKA, which could reveal a preferential phosphorylation. Here, we used quantitative assays on three CaV1.2 and four Rad sites. Isothermal titration calorimetry and enzyme kinetics show that there are two Tiers of targets, with CaV1.2 residue Ser1981 and Rad residues Ser25 and Ser272 forming tier one substrates for PKA. These share a common feature with two Arginine residues at specific positions that can anchor the peptide into the substrate binding cleft of PKA. In contrast, PKA shows minimal activity for the other, tier two substrates, characterized by low kcat values and undetectable binding via isothermal titration calorimetry. The existence of two tiers suggests that PKA regulation of the CaV1.2 complex may occur in a graded fashion. We report crystal structures of the PKA catalytic subunit with and without a CaV1.2 and test the importance of several anchoring residues via mutagenesis. Different target sites utilize different anchors, highlighting the plasticity of PKAc to recognize substrates.

Tripartite interactions of PKA catalytic subunit and C-terminal domains of cardiac Ca2+ channel may modulate its β-adrenergic regulation

BACKGROUND

The β-adrenergic augmentation of cardiac contraction, by increasing the conductivity of L-type voltage-gated CaV1.2 channels, is of great physiological and pathophysiological importance. Stimulation of β-adrenergic receptors (βAR) activates protein kinase A (PKA) through separation of regulatory (PKAR) from catalytic (PKAC) subunits. Free PKAC phosphorylates the inhibitory protein Rad, leading to increased Ca2+ influx. In cardiomyocytes, the core subunit of CaV1.2, CaV1.2α1, exists in two forms: full-length or truncated (lacking the distal C-terminus (dCT)). Signaling efficiency is believed to emanate from protein interactions within multimolecular complexes, such as anchoring PKA (via PKAR) to CaV1.2α1 by A-kinase anchoring proteins (AKAPs). However, AKAPs are inessential for βAR regulation of CaV1.2 in heterologous models, and their role in cardiomyocytes also remains unclear.

RESULTS

We show that PKAC interacts with CaV1.2α1 in heart and a heterologous model, independently of Rad, PKAR, or AKAPs. Studies with peptide array assays and purified recombinant proteins demonstrate direct binding of PKAC to two domains in CaV1.2α1-CT: the proximal and distal C-terminal regulatory domains (PCRD and DCRD), which also interact with each other. Data indicate both partial competition and possible simultaneous interaction of PCRD and DCRD with PKAC. The βAR regulation of CaV1.2α1 lacking dCT (which harbors DCRD) was preserved, but subtly altered, in a heterologous model, the Xenopus oocyte.

CONCLUSIONS

We discover direct interactions between PKAC and two domains in CaV1.2α1. We propose that these tripartite interactions, if present in vivo, may participate in organizing the multimolecular signaling complex and fine-tuning the βAR effect in cardiomyocytes.

Allosteric Changes in the Conformational Landscape of Src Kinase upon Substrate Binding

Precise regulation of protein kinase activity is crucial in cell functions, and its loss is implicated in various diseases. The kinase activity is regulated by interconverting active and inactive states in the conformational landscape. However, how protein kinases switch conformations in response to different signals such as the binding at distinct sites remains incompletely understood. Here, we predict the binding mode for the peptide substrate to Src tyrosine kinase using enhanced conformational sampling simulations (totaling 24 μs) and then investigate changes in the conformational landscape upon substrate binding by conducting unbiased molecular dynamics simulations (totaling 50 μs) initiated from the apo and substrate-bound forms. Unexpectedly, the peptide substrate binding significantly facilitates the transitions from active to inactive conformations in which the αC helix is directed outward, the regulatory spine is broken, and the ATP-binding domain is perturbed. We also explore an underlying residue-contact network responsible for the allosteric conformational changes. Our results are in accord with the recent experiments reporting the negative cooperativity between the peptide substrate and ATP binding to tyrosine kinases and will contribute to advancing our understanding of the regulation mechanisms for kinase activity.

Exploring the conformational landscape of protein kinases

Protein kinases are dynamic enzymes that display complex regulatory mechanisms. Although they possess a structurally conserved catalytic domain, significant conformational dynamics are evident both within a single kinase and across different kinases in the kinome. Here, we highlight methods for exploring this conformational space and its dynamics using kinase domains from ABL1 (Abelson kinase), PKA (protein kinase A), AurA (Aurora A), and PYK2 (proline-rich tyrosine kinase 2) as examples. Such experimental approaches combined with AI-driven methods, such as AlphaFold, will yield discoveries about kinase regulation, the catalytic process, substrate specificity, the effect of disease-associated mutations, as well as new opportunities for structure-based drug design.

Parsing Dynamics of Protein Backbone NH and Side-Chain Methyl Groups using Molecular Dynamics Simulations

Experimental NMR spectroscopy and theoretical molecular dynamics (MD) simulations provide complementary insights into protein conformational dynamics and hence into biological function. The present work describes an extensive set of backbone NH and side-chain methyl group generalized order parameters for the Escherichia coli ribonuclease HI (RNH) enzyme derived from 2-μs microsecond MD simulations using the OPLS4 and AMBER-FF19SB force fields. The simulated generalized order parameters are compared with values derived from NMR 15N and 13CH2D spin relaxation measurements. The squares of the generalized order parameters, S2 for the N-H bond vector and Saxis2 for the methyl group symmetry axis, characterize the equilibrium distribution of vector orientations in a molecular frame of reference. Optimal agreement between simulated and experimental results was obtained by averaging S2 or Saxis2 calculated by dividing the simulated trajectories into 50 ns blocks (∼five times the rotational diffusion correlation time for RNH). With this procedure, the median absolute deviations (MAD) between experimental and simulated values of S2 and Saxis2 are 0.030 (NH) and 0.061 (CH3) for OPLS4 and 0.041 (NH) and 0.078 (CH3) for AMBER-FF19SB. The MAD between OPLS4 and AMBER-FF19SB are 0.021 (NH) and 0.072 (CH3). The generalized order parameters for the methyl group symmetry axis can be decomposed into contributions from backbone fluctuations, between-rotamer dihedral angle transitions, and within-rotamer dihedral angle fluctuations. Analysis of the simulation trajectories shows that (i) backbone and side chain conformational fluctuations exhibit little correlation and that (ii) fluctuations within rotamers are limited and highly uniform with values that depend on the number of dihedral angles considered. Low values of Saxis2, indicative of enhanced side-chain flexibility, result from between-rotamer transitions that can be enhanced by increased local backbone flexibility.

The αC-β4 loop controls the allosteric cooperativity between nucleotide and substrate in the catalytic subunit of protein kinase A

Allosteric cooperativity between ATP and substrates is a prominent characteristic of the cAMP-dependent catalytic subunit of protein kinase A (PKA-C). This long-range synergistic action is involved in substrate recognition and fidelity, and it may also regulate PKA's association with regulatory subunits and other binding partners. To date, a complete understanding of this intramolecular mechanism is still lacking. Here, we integrated NMR(Nuclear Magnetic Resonance)-restrained molecular dynamics simulations and a Markov State Model to characterize the free energy landscape and conformational transitions of 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. We validated the second excited state by analyzing the F100A mutant of PKA-C, assessing its structural response to ATP and substrate binding. While PKA-CF100A preserves its catalytic efficiency with Kemptide, this mutation rearranges the αC-β4 loop conformation, interrupting the coupling of the two lobes and abolishing the allosteric binding cooperativity. The highly conserved αC-β4 loop emerges as a pivotal element to control the synergistic binding of nucleotide and substrate, explaining how mutations or insertions near or within this motif affect the function and drug sensitivity in homologous kinases.

Genetic Encoding of Phosphorylated Amino Acids into Proteins

Reversible phosphorylation is a fundamental mechanism for controlling protein function. Despite the critical roles phosphorylated proteins play in physiology and disease, our ability to study individual phospho-proteoforms has been hindered by a lack of versatile methods to efficiently generate homogeneous proteins with site-specific phosphoamino acids or with functional mimics that are resistant to phosphatases. Genetic code expansion (GCE) is emerging as a transformative approach to tackle this challenge, allowing direct incorporation of phosphoamino acids into proteins during translation in response to amber stop codons. This genetic programming of phospho-protein synthesis eliminates the reliance on kinase-based or chemical semisynthesis approaches, making it broadly applicable to diverse phospho-proteoforms. In this comprehensive review, we provide a brief introduction to GCE and trace the development of existing GCE technologies for installing phosphoserine, phosphothreonine, phosphotyrosine, and their mimics, discussing both their advantages as well as their limitations. While some of the technologies are still early in their development, others are already robust enough to greatly expand the range of biologically relevant questions that can be addressed. We highlight new discoveries enabled by these GCE approaches, provide practical considerations for the application of technologies by non-GCE experts, and also identify avenues ripe for further development.

Missense variants in phospholamban and cardiac myosin binding protein identified in patients with a family history and clinical diagnosis of dilated cardiomyopathy

As the genetic landscape of cardiomyopathies continues to expand, the identification of missense variants in disease-associated genes frequently leads to a classification of variant of uncertain significance (VUS). For the proper reclassification of such variants, functional characterization is an important contributor to the proper assessment of pathogenic potential. Several missense variants in the calcium transport regulatory protein phospholamban have been associated with dilated cardiomyopathy. However, >40 missense variants in this transmembrane peptide are currently known and most remain classified as VUS with little clinical information. Similarly, missense variants in cardiac myosin binding protein have been associated with hypertrophic cardiomyopathy. However, hundreds of variants are known and many have low penetrance and are often found in control populations. Herein, we focused on novel missense variants in phospholamban, an Ala15-Thr variant found in a 4-year-old female and a Pro21-Thr variant found in a 60-year-old female, both with a family history and clinical diagnosis of dilated cardiomyopathy. The patients also harbored a Val896-Met variant in cardiac myosin binding protein. The phospholamban variants caused defects in the function, phosphorylation, and dephosphorylation of this calcium transport regulatory peptide, and we classified these variants as potentially pathogenic. The variant in cardiac myosin binding protein alters the structure of the protein. While this variant has been classified as benign, it has the potential to be a low-risk susceptibility variant because of the structural change in cardiac myosin binding protein. Our studies provide new biochemical evidence for missense variants previously classified as benign or VUS.

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

Revealing eEF-2 kinase: recent structural insights into function

The α-kinase eukaryotic elongation factor 2 kinase (eEF-2K) regulates translational elongation by phosphorylating its ribosome-associated substrate, the GTPase eEF-2. eEF-2K is activated by calmodulin (CaM) through a distinctive mechanism unlike that in other CaM-dependent kinases (CAMK). We describe recent structural insights into this unique activation process and examine the effects of specific regulatory signals on this mechanism. We also highlight key unanswered questions to guide future structure-function studies. These include structural mechanisms which enable eEF-2K to interact with upstream/downstream partners and facilitate its integration of diverse inputs, including Ca2+ transients, phosphorylation mediated by energy/nutrient-sensing pathways, pH changes, and metabolites. Answering these questions is key to establishing how eEF-2K harmonizes translation with cellular requirements within the boundaries of its molecular landscape.

Novel LRR-ROC Motif That Links the N- and C-terminal Domains in LRRK2 Undergoes an Order-Disorder Transition Upon Activation

Mutations in LRRK2, a large multi-domain protein kinase, create risk factors for Parkinson's Disease (PD). LRRK2 has seven well-folded domains that include three N-terminal scaffold domains (NtDs) and four C-terminal domains (CtDs). In full-length inactive LRRK2 there is an additional well-folded motif, the LRR-ROC Linker, that lies between the NtDs and the CtDs. This motif, which is stabilized by hydrophobic residues in the LRR and ROC/COR-A domains, is anchored to the C-Lobe of the kinase domain. The LRR-ROC Linker becomes disordered when the NtDs are unleashed from the CtDs following activation by Rab29 or by various PD mutations. A key residue within the LRR-ROC Linker, W1295, sterically blocks access of substrate proteins. The W1295A mutant blocks cis-autophosphorylation of S1292 and reduces phosphorylation of heterologous Rab substrates. GaMD simulations show that the LRR-Linker motif, P + 1 loop and the inhibitory helix in the DYGψ motif are very stable. Finally, in full-length inactive LRRK2 ATP is bound to the kinase domain and GDP:Mg to the GTPase/ROC domain. The fundamentally different mechanisms for binding nucleotide (G-Loop vs P-Loop) are captured by these GaMD simulations. In this model, where ATP binds with low affinity (μM range) to N-Lobe capping residues, the known auto-phosphorylation sites are located in the space that is sampled by the flexible phosphates thus providing a potential mechanism for cis-autophosphorylation.

An NMR portrait of functional and dysfunctional allosteric cooperativity in cAMP-dependent protein kinase A

The cAMP-dependent protein kinase A (PKA) is the archetypical eukaryotic kinase. The catalytic subunit (PKA-C) structure is highly conserved among the AGC-kinase family. PKA-C is a bilobal enzyme with a dynamic N-lobe, harbouring the Adenosine-5'-triphosphate (ATP) binding site and a more rigid helical C-lobe. The substrate-binding groove resides at the interface of the two lobes. A distinct feature of PKA-C is the positive binding cooperativity between nucleotide and substrate. Several PKA-C mutations lead to the development of adenocarcinomas, myxomas, and other rare forms of liver tumours. Nuclear magnetic resonance (NMR) spectroscopy shows that these mutations disrupt the allosteric communication between the two lobes, causing a drastic decrease in binding cooperativity. The loss of cooperativity correlates with changes in substrate fidelity and reduced kinase affinity for the endogenous protein kinase inhibitor (PKI). The similarity between PKI and the inhibitory sequence of the kinase regulatory subunits suggests that the overall mechanism of regulation of the kinase may be disrupted. We surmise that a reduced or obliterated cooperativity may constitute a common trait for both orthosteric and allosteric mutations of PKA-C that may lead to dysregulation and disease.

Activation loop phosphorylation tunes conformational dynamics underlying Pyk2 tyrosine kinase activation

Pyk2 is a multidomain non-receptor tyrosine kinase that undergoes a multistage activation mechanism. Activation is instigated by conformational rearrangements relieving autoinhibitory FERM domain interactions. The kinase autophosphorylates a central linker residue to recruit Src kinase. Pyk2 and Src mutually phosphorylate activation loops to confer full activation. While the mechanisms of autoinhibition are established, the conformational dynamics associated with autophosphorylation and Src recruitment remain unclear. We employ hydrogen/deuterium exchange mass spectrometry and kinase activity profiling to map the conformational dynamics associated with substrate binding and Src-mediated activation loop phosphorylation. Nucleotide engagement stabilizes the autoinhibitory interface, while phosphorylation deprotects both FERM and kinase regulatory surfaces. Phosphorylation organizes active site motifs linking catalytic loop with activation segment. Dynamics of the activation segment anchor propagate to EF/G helices to prevent reversion of the autoinhibitory FERM interaction. We employ targeted mutagenesis to dissect how phosphorylation-induced conformational rearrangements elevate kinase activity above the basal autophosphorylation rate.

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.

Binding and Functional Folding (BFF): A Physiological Framework for Studying Biomolecular Interactions and Allostery

EF-hand Ca2+-binding proteins (CBPs), such as S100 proteins (S100s) and calmodulin (CaM), are signaling proteins that undergo conformational changes upon increasing intracellular Ca2+. Upon binding Ca2+, S100 proteins and CaM interact with protein targets and induce important biological responses. The Ca2+-binding affinity of CaM and most S100s in the absence of target is weak (CaKD > 1 μM). However, upon effector protein binding, the Ca2+ affinity of these proteins increases via heterotropic allostery (CaKD < 1 μM). Because of the high number and micromolar concentrations of EF-hand CBPs in a cell, at any given time, allostery is required physiologically, allowing for (i) proper Ca2+ homeostasis and (ii) strict maintenance of Ca2+-signaling within a narrow dynamic range of free Ca2+ ion concentrations, [Ca2+]free. In this review, mechanisms of allostery are coalesced into an empirical "binding and functional folding (BFF)" physiological framework. At the molecular level, folding (F), binding and folding (BF), and BFF events include all atoms in the biomolecular complex under study. The BFF framework is introduced with two straightforward BFF types for proteins (type 1, concerted; type 2, stepwise) and considers how homologous and nonhomologous amino acid residues of CBPs and their effector protein(s) evolved to provide allosteric tightening of Ca2+ and simultaneously determine how specific and relatively promiscuous CBP-target complexes form as both are needed for proper cellular function.

MptpA Kinetics Enhanced by Allosteric Control of an Active Conformation

Understanding allostery in the Mycobacterium tuberculosis low molecular weight protein tyrosine phosphatase (MptpA) is a subject of great interest since MptpA is one of two protein tyrosine phosphatases (PTPs) from the pathogenic organism Mycobacterium tuberculosis expressed during host cell infection. Here, we combine computational modeling with solution NMR spectroscopy and we find that Q75 is an allosteric site. Removal of the polar side chain of Q75 by mutation to leucine results in a cascade of events that reposition the acid loop over the active site and relocates the catalytic aspartic acid (D126) at an optimal position for proton donation to the leaving aryl group of the substrate and for subsequent hydrolysis of the thiophosphoryl intermediate. The computational analysis is consistent with kinetic data, and NMR spectroscopy, showing that the Q75L mutant exhibits enhanced reaction kinetics with similar substrate binding affinity. We anticipate that our findings will motivate further studies on the possibility that MptpA remains passivated during the chronic state of infection and increases its activity as part of the pathogenic life cycle of M. tuberculosis possibly via allosteric means.

M, Walker C, Kim J, Camilloni C, De Simone A, Vendruscolo M, Bernlohr DA, Taylor SS, Veglia G. ATP-competitive inhibitors modulate the substrate binding cooperativity of a kinase by altering its conformational entropy

ATP-competitive inhibitors are currently the largest class of clinically approved drugs for protein kinases. By targeting the ATP-binding pocket, these compounds block the catalytic activity, preventing substrate phosphorylation. A problem with these drugs, however, is that inhibited kinases may still recognize and bind downstream substrates, acting as scaffolds or binding hubs for signaling partners. Here, using protein kinase A as a model system, we show that chemically different ATP-competitive inhibitors modulate the substrate binding cooperativity by tuning the conformational entropy of the kinase and shifting the populations of its conformationally excited states. Since we found that binding cooperativity and conformational entropy of the enzyme are correlated, we propose a new paradigm for the discovery of ATP-competitive inhibitors, which is based on their ability to modulate the allosteric coupling between nucleotide and substrate-binding sites.

Raf Kinase Inhibitory Protein regulates the cAMP-dependent protein kinase signaling pathway through a positive feedback loop

Raf Kinase Inhibitory Protein (RKIP) maintains cellular robustness and prevents the progression of diseases such as cancer and heart disease by regulating key kinase cascades including MAP kinase and protein kinase A (PKA). Phosphorylation of RKIP at S153 by Protein Kinase C (PKC) triggers a switch from inhibition of Raf to inhibition of the G protein coupled receptor kinase 2 (GRK2), enhancing signaling by the β-adrenergic receptor (β-AR) that activates PKA. Here we report that PKA-phosphorylated RKIP promotes β-AR-activated PKA signaling. Using biochemical, genetic, and biophysical approaches, we show that PKA phosphorylates RKIP at S51, increasing S153 phosphorylation by PKC and thereby triggering feedback activation of PKA. The S51V mutation blocks the ability of RKIP to activate PKA in prostate cancer cells and to induce contraction in primary cardiac myocytes in response to the β-AR activator isoproterenol, illustrating the functional importance of this positive feedback circuit. As previously shown for other kinases, phosphorylation of RKIP at S51 by PKA is enhanced upon RKIP destabilization by the P74L mutation. These results suggest that PKA phosphorylation at S51 may lead to allosteric changes associated with a higher-energy RKIP state that potentiates phosphorylation of RKIP at other key sites. This allosteric regulatory mechanism may have therapeutic potential for regulating PKA signaling in disease states.

Mapping the conformational energy landscape of Abl kinase using ClyA nanopore tweezers

Protein kinases play central roles in cellular regulation by catalyzing the phosphorylation of target proteins. Kinases have inherent structural flexibility allowing them to switch between active and inactive states. Quantitative characterization of kinase conformational dynamics is challenging. Here, we use nanopore tweezers to assess the conformational dynamics of Abl kinase domain, which is shown to interconvert between two major conformational states where one conformation comprises three sub-states. Analysis of kinase-substrate and kinase-inhibitor interactions uncovers the functional roles of relevant states and enables the elucidation of the mechanism underlying the catalytic deficiency of an inactive Abl mutant G321V. Furthermore, we obtain the energy landscape of Abl kinase by quantifying the population and transition rates of the conformational states. These results extend the view on the dynamic nature of Abl kinase and suggest nanopore tweezers can be used as an efficient tool for other members of the human kinome.

Quantitative characterization of kinase conformational dynamics remains challenging. Here, the authors show that protein nanopore tweezers allow analyzing the conformational energy landscape and ligand binding of the Abl kinase domain.

Quantitative Analysis of Dynamic Allostery

Dynamic allostery refers to one important class of allosteric regulation that does not involve noticeable conformational changes upon effector binding. In recent years, many "quasi"-dynamic allosteric proteins have been found to only experience subtle conformational changes during allosteric regulation. However, as enthalpic and entropic contributions are coupled to each other and even tiny conformational changes could bring in noticeable free energy changes, a quantitative description is essential to understand the contribution of pure dynamic allostery. Here, by developing a unified anisotropic elastic network model (uANM) considering both side-chain information and ligand heavy atoms, we quantitatively estimated the contribution of pure dynamic allostery in a dataset of known allosteric proteins by excluding the conformational changes upon ligand binding. We found that the contribution of pure dynamic allostery is generally small (much weaker than previously expected) and robustly exhibits an allosteric activation effect, which exponentially decays with the distance between the substrate and the allosteric ligand. We further constructed toy models to study the determinant factors of dynamic allostery in monomeric and oligomeric proteins using the uANM. Analysis of the toy models revealed that a short distance, a small angle between the two ligands, strong protein-ligand interactions, and weak protein internal interactions lead to strong dynamic allostery. Our study provides a quantitative estimation of pure dynamic allostery and facilitates the understanding of dynamic-allostery-controlled biological processes and the design of allosteric drugs and proteins.

Structures of PKA-phospholamban complexes reveal a mechanism of familial dilated cardiomyopathy

Several mutations identified in phospholamban (PLN) have been linked to familial dilated cardiomyopathy (DCM) and heart failure, yet the underlying molecular mechanism remains controversial. PLN interacts with sarco/endoplasmic reticulum Ca²⁺-ATPase (SERCA) and regulates calcium uptake, which is modulated by the protein kinase A (PKA)-dependent phosphorylation of PLN during the fight-or-flight response. Here, we present the crystal structures of the catalytic domain of mouse PKA in complex with wild-type and DCM-mutant PLNs. Our structures, combined with the results from other biophysical and biochemical assays, reveal a common disease mechanism: the mutations in PLN reduce its phosphorylation level by changing its conformation and weakening its interactions with PKA. In addition, we demonstrate that another more ubiquitous SERCA-regulatory peptide, called another-regulin (ALN), shares a similar mechanism mediated by PKA in regulating SERCA activity.

"ADPKD-omics": determinants of cyclic AMP levels in renal epithelial cells

The regulation of cyclic adenosine monophosphate (cAMP) levels in kidney epithelial cells is important in at least 2 groups of disorders, namely water balance disorders and autosomal dominant polycystic kidney disease. Focusing on the latter, we review genes that code for proteins that are determinants of cAMP levels in cells. We identify which of these determinants are expressed in the 14 kidney tubule segments using recently published RNA-sequencing and protein mass spectrometry data ("autosomal dominant polycystic kidney disease-omics"). This includes G protein-coupled receptors, adenylyl cyclases, cyclic nucleotide phosphodiesterases, cAMP transporters, cAMP-binding proteins, regulator of G protein-signaling proteins, G protein-coupled receptor kinases, arrestins, calcium transporters, and calcium-binding proteins. In addition, compartmentalized cAMP signaling in the primary cilium is discussed, and a specialized database of the proteome of the primary cilium of cultured "IMCD3" cells is provided as an online resource (https://esbl.nhlbi.nih.gov/Databases/CiliumProteome/). Overall, this article provides a general resource in the form of a curated list of proteins likely to play roles in determination of cAMP levels in kidney epithelial cells and, therefore, likely to be determinants of progression of autosomal dominant polycystic kidney disease.

Physics of biomolecular recognition and conformational dynamics

Biomolecular recognition usually leads to the formation of binding complexes, often accompanied by large-scale conformational changes. This process is fundamental to biological functions at the molecular and cellular levels. Uncovering the physical mechanisms of biomolecular recognition and quantifying the key biomolecular interactions are vital to understand these functions. The recently developed energy landscape theory has been successful in quantifying recognition processes and revealing the underlying mechanisms. Recent studies have shown that in addition to affinity, specificity is also crucial for biomolecular recognition. The proposed physical concept of intrinsic specificity based on the underlying energy landscape theory provides a practical way to quantify the specificity. Optimization of affinity and specificity can be adopted as a principle to guide the evolution and design of molecular recognition. This approach can also be used in practice for drug discovery using multidimensional screening to identify lead compounds. The energy landscape topography of molecular recognition is important for revealing the underlying flexible binding or binding-folding mechanisms. In this review, we first introduce the energy landscape theory for molecular recognition and then address four critical issues related to biomolecular recognition and conformational dynamics: (1) specificity quantification of molecular recognition; (2) evolution and design in molecular recognition; (3) flexible molecular recognition; (4) chromosome structural dynamics. The results described here and the discussions of the insights gained from the energy landscape topography can provide valuable guidance for further computational and experimental investigations of biomolecular recognition and conformational dynamics.

The structural and functional workings of KEOPS

KEOPS (Kinase, Endopeptidase and Other Proteins of Small size) is a five-subunit protein complex that is highly conserved in eukaryotes and archaea and is essential for the fitness of cells and for animal development. In humans, mutations in KEOPS genes underlie Galloway-Mowat syndrome, which manifests in severe microcephaly and renal dysfunction that lead to childhood death. The Kae1 subunit of KEOPS catalyzes the universal and essential tRNA modification N6-threonylcarbamoyl adenosine (t6A), while the auxiliary subunits Cgi121, the kinase/ATPase Bud32, Pcc1 and Gon7 play a supporting role. Kae1 orthologs are also present in bacteria and mitochondria but function in distinct complexes with proteins that are not related in structure or function to the auxiliary subunits of KEOPS. Over the past 15 years since its discovery, extensive study in the KEOPS field has provided many answers towards understanding the roles that KEOPS plays in cells and in human disease and how KEOPS carries out these functions. In this review, we provide an overview into recent advances in the study of KEOPS and illuminate exciting future directions.

Dynamic equilibria in protein kinases

Structural changes involved in protein kinase activation and ligand binding have been determined from a wealth of X-ray crystallographic evidence. Recent solution studies using NMR, EPR, HX-MS, and fluorescence techniques have deepened this understanding by highlighting the underlying energetics and dynamics of multistate conformational ensembles. This new research is showing how activation mechanisms and ligand binding alter the internal motions of kinases and enable allosteric coupling between distal regulatory regions and the active site.

Everything you ever wanted to know about PKA regulation and its involvement in mammalian sperm capacitation

The 3', 5'-cyclic adenosine monophosphate (cAMP) dependent protein kinase (PKA) is a tetrameric holoenzyme comprising a set of two regulatory subunits (PKA-R) and two catalytic (PKA-C) subunits. The PKA-R subunits act as sensors of cAMP and allow PKA-C activity. One of the first signaling events observed during mammalian sperm capacitation is PKA activation. Thus, understanding how PKA activity is restricted in space and time is crucial to decipher the critical steps of sperm capacitation. It is widely accepted that PKA specificity depends on several levels of regulation. Anchoring proteins play a pivotal role in achieving proper localization signaling, subcellular targeting and cAMP microdomains. These multi-factorial regulation steps are necessary for a precise spatio-temporal activation of PKA. Here we discuss recent understanding of regulatory mechanisms of PKA in mammalian sperm, such as post-translational modifications, in the context of its role as the master orchestrator of molecular events conducive to capacitation.

Revealing the Unbinding Kinetics and Mechanism of Type I and Type II Protein Kinase Inhibitors by Local-Scaled Molecular Dynamics Simulations

Protein kinase inhibitors disrupt phosphorylation of the target kinases, which are an important class of drug for treating cancer and other diseases. Conventional structure-based design methods (such as molecular docking) focus on the static binding mode of the kinase inhibitor with its target. However, dissociation kinetic properties of a drug molecule are found to correlate with its residence time in vivo and thus have drawn the attention of drug designers in recent years. In this study, we have applied the local-scaled molecular dynamics (MD) simulation enabled in GROMACS software to explore the unbinding mechanism of a total of 41 type I and type II kinase inhibitors. Our simulation considered multiple starting configurations as well as possible protonation states of kinase inhibitors. Based on our local-scaled MD results, we discovered that the integrals of the favorable binding energy during dissociation correlated well (R² = 0.64) with the experimental dissociation rate constants of those kinase inhibitors on the entire data set. Given its accuracy and technical advantage, this method may serve as a practical option for estimating this important property in reality. Our simulation also provided a reasonable explanation of the dynamic properties of kinase and its inhibitor as well as the role of relevant water molecules in dissociation.

Inter-lobe Motions Allosterically Regulate the Structure and Function of EGFR Kinase

Protein kinases play important roles in cellular signaling and have been one of the best-studied drug targets. The kinase domain of epidermal growth factor receptor (EGFR) is a receptor tyrosine kinase that has been extensively studied for cancer drug discovery and for understanding the unique activation mechanism by dimerization. Here, we analyzed all available 206 crystal structures of the EGFR kinase and the dynamics observed in molecular simulations to identify how these structures are determined. It was found that the arrangement between the N- and C-terminal lobes plays a key role in regulating the kinase structure by sensitively responding to the intermolecular interactions, or the crystal environment. A whole variety of crystal forms in the database is thus reflected in the broad distribution of the inter-lobe arrangement. The configuration of the catalytically important motifs as well as the bound ATP is closely coupled with the inter-lobe motion. When the intermolecular interactions are those of the activating asymmetric dimer, EGFR kinase takes the open-lobe arrangement that constructs the catalytically active configuration.

Multi-state recognition pathway of the intrinsically disordered protein kinase inhibitor by protein kinase A

In the nucleus, the spatiotemporal regulation of the catalytic subunit of cAMP-dependent protein kinase A (PKA-C) is orchestrated by an intrinsically disordered protein kinase inhibitor, PKI, which recruits the CRM1/RanGTP nuclear exporting complex. How the PKA-C/PKI complex assembles and recognizes CRM1/RanGTP is not well understood. Using NMR, SAXS, fluorescence, metadynamics, and Markov model analysis, we determined the multi-state recognition pathway for PKI. After a fast binding step in which PKA-C selects PKI's most competent conformations, PKI folds upon binding through a slow conformational rearrangement within the enzyme's binding pocket. The high-affinity and pseudo-substrate regions of PKI become more structured and the transient interactions with the kinase augment the helical content of the nuclear export sequence, which is then poised to recruit the CRM1/RanGTP complex for nuclear translocation. The multistate binding mechanism featured by PKA-C/PKI complex represents a paradigm on how disordered, ancillary proteins (or protein domains) are able to operate multiple functions such as inhibiting the kinase while recruiting other regulatory proteins for nuclear export.

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.

Molecular Basis of the Mechanisms Controlling MASTL

The human MASTL (Microtubule-associated serine/threonine kinase-like) gene encodes an essential protein in the cell cycle. MASTL is a key factor preventing early dephosphorylation of M-phase targets of Cdk1/CycB. Little is known about the mechanism of MASTL activation and regulation. MASTL contains a non-conserved insertion of 550 residues within its activation loop, splitting the kinase domain, and making it unique. Here, we show that this non-conserved middle region (NCMR) of the protein is crucial for target specificity and activity. We performed a phosphoproteomic assay with different MASTL constructs identifying key phosphorylation sites for its activation and determining whether they arise from autophosphorylation or exogenous kinases, thus generating an activation model. Hydrogen/deuterium exchange data complements this analysis revealing that the C-lobe in full-length MASTL forms a stable structure, whereas the N-lobe is dynamic and the NCMR and C-tail contain few localized regions with higher-order structure. Our results indicate that truncated versions of MASTL conserving a cryptic C-Lobe in the NCMR, display catalytic activity and different targets, thus establishing a possible link with truncated mutations observed in cancer-related databases.

NMR and computational methods for molecular resolution of allosteric pathways in enzyme complexes

Allostery is a ubiquitous biological mechanism in which a distant binding site is coupled to and drastically alters the function of a catalytic site in a protein. Allostery provides a high level of spatial and temporal control of the integrity and activity of biomolecular assembles composed of proteins, nucleic acids, or small molecules. Understanding the physical forces that drive allosteric coupling is critical to harnessing this process for use in bioengineering, de novo protein design, and drug discovery. Current microscopic models of allostery highlight the importance of energetics, structural rearrangements, and conformational fluctuations, and in this review, we discuss the synergistic use of solution NMR spectroscopy and computational methods to probe these phenomena in allosteric systems, particularly protein-nucleic acid complexes. This combination of experimental and theoretical techniques facilitates an unparalleled detection of subtle changes to structural and dynamic equilibria in biomolecules with atomic resolution, and we provide a detailed discussion of specialized NMR experiments as well as the complementary methods that provide valuable insight into allosteric pathways in silico. Lastly, we highlight two case studies to demonstrate the adaptability of this approach to enzymes of varying size and mechanistic complexity.

A Dynamic Switch in Inactive p38γ Leads to an Excited State on the Pathway to an Active Kinase

The inactive state of mitogen-activated protein kinases (MAPKs) adopts an open conformation while the active state exists in a compact form stabilized by phosphorylation. In the active state, eukaryotic kinases undergo breathing motions related to substrate binding and product release that have not previously been detected in the inactive state. However, docking interactions of partner proteins with inactive MAPK kinases exhibit allostery in binding of activating kinases. Interactions at a site distant from the activation loop are coupled to the configuration of the activation loop, suggesting that the inactive state may also undergo concerted dynamics. X-ray crystallographic studies of nonphosphorylated, inactive p38γ reveal differences in domain orientations and active site structure in the two molecules in the asymmetric unit. One molecule resembles an inactive kinase with an open active site. The second molecule has a rotation of the N-lobe that leads to partial compaction of the active site, resulting in a conformation that is intermediate between the inactive open state and the fully closed state of the activated kinase. Although the compact state of apo p38γ displays several of the features of the activated enzyme, it remains catalytically inert. In solution, the kinase fluctuates on a millisecond time scale between the open ground state and a weakly populated excited state that is similar in structure to the compact state observed in the crystal. The nuclear magnetic resonance and crystal structure data imply that interconversion between the open and compact states involves a molecular switch associated with the DFG loop.

Characterization of the dynamics and the conformational entropy in the binding between TAZ1 and CTAD-HIF-1α

The interaction between the C-terminal transactivation domain of HIF-1α (CTAD-HIF-1α) and the transcriptional adapter zinc binding 1 (TAZ1) domain of CREB binding protein participate in the initiation of gene transcription during hypoxia. Unbound CTAD-HIF-1α is disordered but undergoes a disorder-to-order transition upon binding to TAZ1. We have here performed NMR side chain and backbone relaxation studies on TAZ1 and side chain relaxation measurements on CTAD-HIF-1α in order to investigate the role of picosecond to nanosecond dynamics. We find that the internal motions are significantly affected upon binding, both on the side chain and the backbone level. The dynamic response corresponds to a conformational entropy change that contributes substantially to the binding thermodynamics for both binding partners. Furthermore, the conformational entropy change for the well-folded TAZ1 varies upon binding to different IDP targets. We further identify a cluster consisting of side chains in bound TAZ1 and CTAD-HIF-1α that experience extensive dynamics and are part of the binding region that involves the N-terminal end of the LPQL motif in CTAD-HIF-1α; a feature that might have an important role in the termination of the hypoxic response.

Universality of fold-encoded localized vibrations in enzymes

Enzymes speed up biochemical reactions at the core of life by as much as 15 orders of magnitude. Yet, despite considerable advances, the fine dynamical determinants at the microscopic level of their catalytic proficiency are still elusive. In this work, we use a powerful mathematical approach to show that rate-promoting vibrations in the picosecond range, specifically encoded in the 3D protein structure, are localized vibrations optimally coupled to the chemical reaction coordinates at the active site. Remarkably, our theory also exposes an hithertho unknown deep connection between the unique localization fingerprint and a distinct partition of the 3D fold into independent, foldspanning subdomains that govern long-range communication. The universality of these features is demonstrated on a pool of more than 900 enzyme structures, comprising a total of more than 10,000 experimentally annotated catalytic sites. Our theory provides a unified microscopic rationale for the subtle structure-dynamics-function link in proteins.

Cushing's syndrome driver mutation disrupts protein kinase A allosteric network, altering both regulation and substrate specificity

Genetic alterations in the PRKACA gene coding for the catalytic α subunit of the cAMP-dependent protein kinase A (PKA-C) are linked to cortisol-secreting adrenocortical adenomas, resulting in Cushing's syndrome. Among those, a single mutation (L205R) has been found in up to 67% of patients. Because the x-ray structures of the wild-type and mutant kinases are essentially identical, the mechanism explaining aberrant function of this mutant remains under active debate. Using NMR spectroscopy, thermodynamics, kinetic assays, and molecular dynamics simulations, we found that this single mutation causes global changes in the enzyme, disrupting the intramolecular allosteric network and eliciting losses in nucleotide/pseudo-substrate binding cooperativity. Remarkably, by rewiring its internal allosteric network, PKA-CL205R is able to bind and phosphorylate non-canonical substrates, explaining its changes in substrate specificity. Both the lack of regulation and change in substrate specificity reveal the complex role of this mutated kinase in the formation of cortisol-secreting adrenocortical adenomas.

Dynamic allostery-based molecular workings of kinase:peptide complexes

A dense interplay between structure and dynamics underlies the working of proteins, especially enzymes. Protein kinases are molecular switches that are optimized for their regulation rather than catalytic turnover rates. Using long-simulations dynamic allostery analysis, this study describes an exploration of the dynamic kinase:peptide complex. We have used protein kinase A (PKA) as a model system as a generic prototype of the protein kinase superfamily of signaling enzymes. Our results explain the role of dynamic coupling of active-site residues that must work in coherence to provide for a successful activation or inhibition response from the kinase. Amino acid networks-based community analysis allows us to ponder the conformational entropy of the kinase:nucleotide:peptide ternary complex. We use a combination of 7 peptides that include 3 types of PKA-binding partners: Substrates, products, and inhibitors. The substrate peptides provide for dynamic insights into the enzyme:substrate complex, while the product phospho-peptide allows for accessing modes of enzyme:product release. Mapping of allosteric communities onto the PKA structure allows us to locate the more unvarying and flexible dynamic regions of the kinase. These distributions, when correlated with the structural elements of the kinase core, allow for a detailed exploration of key dynamics-based signatures that could affect peptide recognition and binding at the kinase active site. These studies provide a unique dynamic allostery-based perspective to kinase:peptide complexes that have previously been explored only in a structural or thermodynamic context.

AKAP6 and phospholamban colocalize and interact in HEK-293T cells and primary murine cardiomyocytes

Phospholamban (PLN) is an important Ca²⁺ modulator at the sarcoplasmic reticulum (SR) of striated muscles. It physically interacts and inhibits sarcoplasmic reticulum Ca²⁺ATPase (SERCA2) function, whereas a protein kinase A (PKA)‐dependent phosphorylation at its serine 16 reverses the inhibition. The underlying mechanism of this post‐translational modification, however, remains not fully understood. Using publicly available databases, we identified A‐kinase anchoring protein 6 (AKAP6) as a candidate that might play some roles in PLN phosphorylation. Immunofluorescence showed colocalization between GFP‐AKAP6 and PLN in transfected HEK‐293T cells and cultured mouse neonatal cardiomyocytes (CMNCs). Co‐immunoprecipitation confirmed the functional interaction between AKAP6 and PLN in HEK‐293T and isolated adult rat cardiomyocytes in response to isoproterenol stimulation. Functionally, AKAP6 promoted Ca²⁺ uptake activity of SERCA1 in cotransfected HEK‐293T cells despite the presence of PLN. These results were further confirmed in adult rat cardiomyocytes. Immunofluorescence showed colocalization of both proteins around the perinuclear region, while protein–protein interaction was corroborated by immunoprecipitation of the nucleus‐enriched fraction of rat hearts. Our findings suggest AKAP6 as a novel interacting partner to PLN in HEK‐293T and murine cardiomyocytes.

The Cardiac Ryanodine Receptor Phosphorylation Hotspot Embraces PKA in a Phosphorylation-Dependent Manner

Ryanodine receptors (RyRs) are intracellular Ca²⁺ release channels controlling essential cellular functions. RyRs are targeted by cyclic AMP (cAMP)-dependent protein kinase A (PKA), a controversial regulation implicated in disorders ranging from heart failure to Alzheimer's. Using crystal structures, we show that the phosphorylation hotspot domain of RyR2 embraces the PKA catalytic subunit, with an extensive interface not seen in PKA complexes with peptides. We trapped an intermediary open-form PKA bound to the RyR2 domain and an ATP analog, showing that PKA can engage substrates in an open form. Phosphomimetics or prior phosphorylation at nearby sites in RyR2 either enhance or reduce the activity of PKA. Finally, we show that a phosphomimetic at S2813, a well-known target site for calmodulin-dependent kinase II, induces the formation of an alpha helix in the phosphorylation domain, resulting in increased interactions and PKA activity. This shows that the different phosphorylation sites in RyR2 are not independent.

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.

Globally correlated conformational entropy underlies positive and negative cooperativity in a kinase's enzymatic cycle

Enzymes accelerate the rate of chemical transformations by reducing the activation barriers of uncatalyzed reactions. For signaling enzymes, substrate recognition, binding, and product release are often rate-determining steps in which enthalpy-entropy compensation plays a crucial role. While the nature of enthalpic interactions can be inferred from structural data, the molecular origin and role of entropy in enzyme catalysis remains poorly understood. Using thermocalorimetry, NMR, and MD simulations, we studied the conformational landscape of the catalytic subunit of cAMP-dependent protein kinase A, a ubiquitous phosphoryl transferase involved in a myriad of cellular processes. Along the enzymatic cycle, the kinase exhibits positive and negative cooperativity for substrate and nucleotide binding and product release. We found that globally coordinated changes of conformational entropy activated by ligand binding, together with synchronous and asynchronous breathing motions of the enzyme, underlie allosteric cooperativity along the kinase's cycle.

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.

Effects of the Arg9Cys and Arg25Cys mutations on phospholamban's conformational equilibrium in membrane bilayers

Approximately, 70% of the Ca2+ ion transport into the sarcoplasmic reticulum is catalyzed by the sarcoplasmic reticulum Ca2+-ATPase (SERCA), whose activity is endogenously regulated by phospholamban (PLN). PLN comprises a TM inhibitory region and a cytoplasmic regulatory region that harbors a consensus sequence for cAMP-dependent protein kinase (PKA). The inhibitory region binds the ATPase, reducing its apparent Ca2+ binding affinity. β-adrenergic stimulation activates PKA, which phosphorylates PLN at Ser 16, reversing its inhibitory function. Mutations and post-translational modifications of PLN may lead to dilated cardiomyopathy (DCM) and heart failure. PLN's cytoplasmic region interconverts between a membrane-associated T state and a membrane-detached R state. The importance of these structural transitions on SERCA regulation is emerging, but the effects of natural occurring mutations and their relevance to the progression of heart disease are unclear. Here we use solid-state NMR spectroscopy to investigate the structural dynamics of two lethal PLN mutations, R9C and R25C, which lead to DCM. We found that the R25C mutant enhances the dynamics of PLN and shifts the conformational equilibrium toward the R state confirmation, whereas the R9C mutant drives the amphipathic cytoplasmic domain toward the membrane-associate state, enriching the T state population. The changes in membrane interactions caused by these mutations may explain the aberrant regulation of SERCA.

Conformational Plasticity in Tyrosine Kinase Inhibitor-Kinase Interactions Revealed with Fluorescence Spectroscopy and Theoretical Calculations

To understand drug-protein dynamics, it is necessary to account for drug molecular flexibility and binding site plasticity. Herein, we exploit fluorescence from a tyrosine kinase inhibitor, AG1478, as a reporter of its conformation and binding site environment when complexed with its cognate kinase. Water-soluble kinases, aminoglycoside phosphotransferase APH(3')-Ia and mitogen-activated protein kinase 14 (MAPK14), were chosen for this study. On the basis of our prior work, the AG1478 conformation (planar or twisted) was inferred from the fluorescence excitation spectrum and the polarity of the AG1478-binding site was deduced from the fluorescence emission spectrum, while red-edge excitation shift (REES) probed the heterogeneity of the binding site (protein conformation and hydration) distributions in the protein conformational ensemble. In the AG1478-APH(3')-Ia complex, both twisted (or partially twisted) and planar AG1478 conformations were evidenced from emission wavelength-dependent excitation spectra. The binding site environment provided by APH(3')-Ia was moderately polar (λ_max = 480 nm) with evidence for considerable heterogeneity (REES = 34 nm). In contrast, in the AG1478-MAPK14 complex, AG1478 was in a predominantly planar conformation with a lower degree of conformational heterogeneity. The binding site environment provided by the MAPK14 protein was of relatively low polarity (λ_max = 430 nm) with a smaller degree of heterogeneity (REES = 11 nm). The results are compared with the available X-ray data and discussed in the context of our current understanding of tyrosine kinase inhibitor conformation and protein conformational ensembles.

Dynamic activation and regulation of the mitogen-activated protein kinase p38

Mitogen-activated protein kinases, which include p38, are essential for cell differentiation and autophagy. The current model for p38 activation involves activation-loop phosphorylation with subsequent substrate binding leading to substrate phosphorylation. Despite extensive efforts, the molecular mechanism of activation remains unclear. Here, using NMR spectroscopy, we show how the modulation of protein dynamics across timescales activates p38. We find that activation-loop phosphorylation does not change the average conformation of p38; rather it quenches the loop ps-ns dynamics. We then show that substrate binding to nonphosphorylated and phosphorylated p38 results in uniform µs-ms backbone dynamics at catalytically essential regions and across the entire molecule, respectively. Together, these results show that phosphorylation and substrate binding flatten the energy landscape of the protein, making essential elements of allostery and activation dynamically accessible. The high degree of structural conservation among ser/thr kinases suggests that elements of this mechanism may be conserved across the kinase family.

Water-mediated conformational preselection mechanism in substrate binding cooperativity to protein kinase A

Substrate binding cooperativity in protein kinase A (PKA) seems to involve allosteric coupling between the two binding sites. It received significant attention, but its molecular basis still remains not entirely clear. Based on long molecular dynamics of PKA and its complexes, we characterized an allosteric pathway that links ATP binding to the redistribution of states adopted by a protein substrate positioning segment in favor of those that warrant correct binding. We demonstrate that the cooperativity mechanism critically depends on the presence of water in two distinct, buried hydration sites. One holds just a single water molecule, which acts as a switchable hydrogen bond bridge along the allosteric pathway. The second, filled with partially disordered solvent, is essential for providing a smooth free energy landscape underlying conformational transitions of the peptide binding region. Our findings remain in agreement with experimental data, also concerning the cooperativity abolishing effect of the Y204A mutation, and indicate a plausible molecular mechanism contributing to experimentally observed binding cooperativity of the two substrates.