How does the cAMP-dependent protein kinase catalyze the phosphorylation reaction: an ab initio QM/MM study

SERKs serve as co-receptors for SYR1 to trigger systemin-mediated defense responses in tomato

Systemin, the first peptide hormone identified in plants, was initially isolated from tomato (Solanum lycopersicum) leaves. Systemin mediates local and systemic wound-induced defense responses in plants, conferring resistance to necrotrophic fungi and herbivorous insects. Systemin is recognized by the leucine-rich-repeat receptor-like kinase (LRR-RLK) receptor SYSTEMIN RECEPTOR1 (SYR1), but how the systemin recognition signal is transduced to intracellular signaling pathways to trigger defense responses is poorly understood. Here, we demonstrate that SERK family LRR-RLKs function as co-receptors for SYR1 to mediate systemin signal transduction in tomato. By using chemical genetic approaches coupled with engineered receptors, we revealed that the association of the cytoplasmic kinase domains of SYR1 with SERKs leads to their mutual trans-phosphorylation and the activation of SYR1, which in turn induces a wide range of defense responses. Systemin stimulates the association between SYR1 and all tomato SERKs (SlSERK1, SlSERK3A, and SlSERK3B). The resulting SYR1-SlSERK heteromeric complexes trigger the phosphorylation of TOMATO PROTEIN KINASE 1B (TPK1b), a receptor-like cytoplasmic kinase that positively regulates systemin responses. Additionally, upon association with SYR1, SlSERKs are cleaved by the Pseudomonas syringae effector HopB1, further supporting the finding that SlSERKs are activated by systemin-bound SYR1. Finally, genetic analysis using Slserk mutants showed that SlSERKs are essential for systemin-mediated defense responses. Collectively, these findings demonstrate that the systemin-mediated association of SYR1 and SlSERKs activates defense responses against herbivorous insects.

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.

Specificity and Reactivity of Mycobacterium tuberculosis Serine/Threonine Kinases PknG and PknB

Mycobacterium tuberculosis (Mtb), the causative agent of Tuberculosis, has 11 eukaryotic-like serine/threonine protein kinases, which play essential roles in cell growth, signal transduction, and pathogenesis. Protein kinase G (PknG) regulates the carbon and nitrogen metabolism by phosphorylation of the glycogen accumulation regulator (GarA) protein at Thr21. Protein kinase B (PknB) is involved in cell wall synthesis and cell shape, as well as phosphorylates GarA but at Thr22. While PknG seems to be constitutively activated and recognition of GarA requires phosphorylation in its unstructured tail, PknB activation is triggered by phosphorylation of its activation loop, which allows binding of the forkhead-associated domain of GarA. In the present work, we used molecular dynamics and quantum-mechanics/molecular mechanics simulations of the catalytically competent complex and kinase activity assays to understand PknG/PknB specificity and reactivity toward GarA. Two hydrophobic residues in GarA, Val24 and Phe25, seem essential for PknG binding and allow specificity for Thr21 phosphorylation. On the other hand, phosphorylated residues in PknB bind Arg26 in GarA and regulate its specificity for Thr22. We also provide a detailed analysis of the free energy profile for the phospho-transfer reaction and show why PknG has a constitutively active conformation not requiring priming phosphorylation in contrast to PknB. Our results provide new insights into these two key enzymes relevant for Mtb and the mechanisms of serine/threonine phosphorylation in bacteria.

Activation loop phosphorylaton of a non-RD receptor kinase initiates plant innate immune signaling

Receptor kinases (RKs) are fundamental for extracellular sensing and regulate development and stress responses across kingdoms. In plants, leucine-rich repeat receptor kinases (LRR-RKs) are primarily peptide receptors that regulate responses to myriad internal and external stimuli. Phosphorylation of LRR-RK cytoplasmic domains is among the earliest responses following ligand perception, and reciprocal transphosphorylation between a receptor and its coreceptor is thought to activate the receptor complex. Originally proposed based on characterization of the brassinosteroid receptor, the prevalence of complex activation via reciprocal transphosphorylation across the plant RK family has not been tested. Using the LRR-RK ELONGATION FACTOR TU RECEPTOR (EFR) as a model, we set out to understand the steps critical for activating RK complexes. While the EFR cytoplasmic domain is an active protein kinase in vitro and is phosphorylated in a ligand-dependent manner in vivo, catalytically deficient EFR variants are functional in antibacterial immunity. These results reveal a noncatalytic role for EFR in triggering immune signaling and indicate that reciprocal transphoshorylation is not a ubiquitous requirement for LRR-RK complex activation. Rather, our analysis of EFR along with a detailed survey of the literature suggests a distinction between LRR-RKs with RD- versus non-RD protein kinase domains. Based on newly identified phosphorylation sites that regulate the activation state of the EFR complex in vivo, we propose that LRR-RK complexes containing a non-RD protein kinase may be regulated by phosphorylation-dependent conformational changes of the ligand-binding receptor, which could initiate signaling either allosterically or through driving the dissociation of negative regulators of the complex.

Contrasting Effects of Inhibitors Li+ and Be2+ on Catalytic Cycle of Glycogen Synthase Kinase-3β

Ionic lithium shows rare effectiveness for treating bipolar disorder and is a potential drug for neurodegenerative diseases. Unfortunately, lithium suffers from significant drawbacks, mainly a narrow therapeutic window. Among the targets of lithium, glycogen synthase kinase 3β (GSK-3β) may be responsible for its therapeutic effects. The development of alternative, selective inhibitors of this kinase could prevent lithium side effects, but such efforts have met little success so far. An atomistic understanding of Li⁺ inhibition and the GSK-3β phosphorylation reaction would therefore facilitate the development of new drugs. In this study, we use extensive sampling of catalytic states with our mixed quantum-classical dynamics method QM/DMD and binding affinities from a competitive metal affinity (CMA) approach to expand the atomistic picture of Li⁺ GSK-3β inhibition. We compare Li⁺ action with Be²⁺ and find our results in agreement with in vitro kinetics studies. Ultimately, our simulations show that Li⁺ inhibition is driven by decreasing the phosphorylation reaction rate, rather than reducing catalytic turnover through tight binding to different GSK-3β states like Be²⁺ inhibition. The effect of these metals derive from electrostatic differences and especially their smaller atomic radii compared to the native Mg²⁺ and thus provide insight for the development of GSK-3β inhibitors based on other paradigms.

Structures of a deAMPylation complex rationalise the switch between antagonistic catalytic activities of FICD

The endoplasmic reticulum (ER) Hsp70 chaperone BiP is regulated by AMPylation, a reversible inactivating post-translational modification. Both BiP AMPylation and deAMPylation are catalysed by a single ER-localised enzyme, FICD. Here we present crystallographic and solution structures of a deAMPylation Michaelis complex formed between mammalian AMPylated BiP and FICD. The latter, via its tetratricopeptide repeat domain, binds a surface that is specific to ATP-state Hsp70 chaperones, explaining the exquisite selectivity of FICD for BiP’s ATP-bound conformation both when AMPylating and deAMPylating Thr518. The eukaryotic deAMPylation mechanism thus revealed, rationalises the role of the conserved Fic domain Glu234 as a gatekeeper residue that both inhibits AMPylation and facilitates hydrolytic deAMPylation catalysed by dimeric FICD. These findings point to a monomerisation-induced increase in Glu234 flexibility as the basis of an oligomeric state-dependent switch between FICD’s antagonistic activities, despite a similar mode of engagement of its two substrates — unmodified and AMPylated BiP.

The ER chaperone BiP is regulated by FICD-mediated AMPylation and deAMPylation. Here, the authors characterise the structure of mammalian AMPylated BiP bound to FICD, by X-ray crystallography and neutron scattering, providing insights into the mechanism of BiP AMPylation and deAMPylation.

ATP-Magnesium Coordination: Protein Structure-Based Force Field Evaluation and Corrections

In the numerous molecular recognition and catalytic processes across biochemistry involving adenosine triphosphate (ATP), the common bioactive form is its magnesium chelate, ATP·Mg2+. In aqueous solution, two chelation geometries predominate, distinguished by bidentate and tridentate Mg2+-phosphate coordination. These are approximately isoenergetic but separated by a high energy barrier. Force field-based atomistic simulation studies of this complex require an accurate representation of its structure and energetics. Here we focused on the energetics of ATP·Mg2+ coordination. Applying an enhanced sampling scheme to circumvent prohibitively slow sampling of transitions between coordination modes, we observed striking contradictions between Amber and CHARMM force field descriptions, most prominently in opposing predictions of the favored coordination mode. Through further configurational free energy calculations, conducted against a diverse set of ATP·Mg2+-protein complex structures to supplement otherwise limited experimental data, we quantified systematic biases for each force field. The force field calculations were strongly predictive of experimentally observed coordination modes, enabling additive corrections to the coordination free energy that deliver close agreement with experiment. We reassessed the applicability of the thus corrected force field descriptions of ATP·Mg2+ for biomolecular simulation and observed that, while the CHARMM parameters display an erroneous preference for overextended triphosphate configurations that will affect many common biomolecular simulation applications involving ATP, the force field energy landscapes broadly agree with experimental measurements of solution geometry and the distribution of ATP·Mg2+ structures found in the Protein Data Bank. Our force field evaluation and correction approach, based on maximizing consistency with the large and heterogeneous collection of structural information encoded in the PDB, should be broadly applicable to many other systems.

Conformational dynamics modulate the catalytic activity of the molecular chaperone Hsp90

The heat shock protein 90 (Hsp90) is a molecular chaperone that employs the free energy of ATP hydrolysis to control the folding and activation of several client proteins in the eukaryotic cell. To elucidate how the local ATPase reaction in the active site couples to the global conformational dynamics of Hsp90, we integrate here large-scale molecular simulations with biophysical experiments. We show that the conformational switching of conserved ion pairs between the N-terminal domain, harbouring the active site, and the middle domain strongly modulates the catalytic barrier of the ATP-hydrolysis reaction by electrostatic forces. Our combined findings provide a mechanistic model for the coupling between catalysis and protein dynamics in Hsp90, and show how long-range coupling effects can modulate enzymatic activity.

The chaperone Hsp90 uses the free energy from ATP hydrolysis to control the folding of client proteins in eukaryotic cells. Here the authors provide mechanistic insights into how its catalytic activity is coupled to conformational changes by combining large-scale molecular simulations with NMR, FRET and SAXS experiments.

Kinase Activation by Small Conformational Changes

Protein kinases (PKs) are allosteric enzymes that play an essential role in signal transduction by regulating a variety of key cellular processes. Most PKs suffer conformational rearrangements upon phosphorylation that strongly enhance the catalytic activity. Generally, it involves the movement of the phosphorylated loop toward the active site and the rotation of the whole C-terminal lobe. However, not all kinases undergo such a large configurational change: The MAPK extracellular signal-regulated protein kinases ERK1 and ERK2 achieve a 50 000 fold increase in kinase activity with only a small motion of the C-terminal region. In the present work, we used a combination of molecular simulation tools to characterize the conformational landscape of ERK2 in the active (phosphorylated) and inactive (unphosphorylated) states in solution in agreement with NMR experiments. We show that the chemical reaction barrier is strongly dependent on ATP conformation and that the "active" low-barrier configuration is subtly regulated by phosphorylation, which stabilizes a key salt bridge between the conserved Lys52 and Glu69 belonging to helix-C and promotes binding of a second Mg ion. Our study highlights that the on-off switch embedded in the kinase fold can be regulated by small, medium, and large conformational changes.

Mechanistic insights into the phosphoryl transfer reaction in cyclin-dependent kinase 2: A QM/MM study

Cyclin-dependent kinase 2 (CDK2) is an important member of the CDK family exerting its most important function in the regulation of the cell cycle. It catalyzes the transfer of the gamma phosphate group from an ATP (adenosine triphosphate) molecule to a Serine/Threonine residue of a peptide substrate. Due to the importance of this enzyme, and protein kinases in general, a detailed understanding of the reaction mechanism is desired. Thus, in this work the phosphoryl transfer reaction catalyzed by CDK2 was revisited and studied by means of hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. Our results suggest that the base-assisted mechanism is preferred over the substrate-assisted pathway when one Mg²⁺ is present in the active site, in agreement with a previous theoretical study. The base-assisted mechanism resulted to be dissociative, with a potential energy barrier of 14.3 kcal/mol, very close to the experimental derived value. An interesting feature of the mechanism is the proton transfer from Lys129 to the phosphoryl group at the second transition state, event that could be helping in neutralizing the charge on the phosphoryl group upon the absence of a second Mg²⁺ ion. Furthermore, important insights into the mechanisms in terms of bond order and charge analysis were provided. These descriptors helped to characterize the synchronicity of bond forming and breaking events, and to characterize charge transfer effects. Local interactions at the active site are key to modulate the charge distribution on the phosphoryl group and therefore alter its reactivity.

Two symmetric arginine residues play distinct roles in Thermus thermophilus Argonaute DNA guide strand-mediated DNA target cleavage

Bacterium Thermus thermophilus Argonaute (Ago; TtAgo) is a prokaryotic Ago (pAgo) that acts as the host defense against the uptake and propagation of foreign DNA by catalyzing the DNA cleavage reaction. The TtAgo active site consists of a plugged-in glutamate finger with two arginine residues (R545 and R486) located symmetrically around it. An interesting challenge is to understand how they can collaboratively facilitate enzymatic catalysis. In Kluyveromyces polysporus Ago, a eukaryotic Ago, the evolutionarily symmetrical residues are arginine and histidine, both of which function to stabilize the plugged-in catalytic tetrad conformation. Surprisingly, our simulation results indicated that, in TtAgo, only R545 is involved in the cleavage reaction by serving as a critical structural anchor to stabilize the catalytic tetrad Asp-Glu-Asp-Asp that is completed by the insertion of the glutamate finger, whereas R486 is not involved in target cleavage. The TtAgo-mediated target DNA cleavage occurs in a substrate-assisted mechanism, in which the pro-Rp (Rp, a tetrahedral phosphorus center with "R-type" chirality) oxygen of scissile phosphate acts as a general base to activate the nucleophilic water. Our unexpected theoretical findings on distinct roles played by R545 and R486 in TtAgo catalysis have been validated by single-point site-mutagenesis experiments, wherein the target cleavage is abolished for all mutants of R545. In sharp contrast, the cleavage activity is maintained for all mutants of R486. Our work provides mechanistic insights on the catalytic specificity of Ago proteins and could facilitate the design of new gene-editing tools in the long term.

Study of the affinity between the protein kinase PKA and homoarginine-containing peptides derived from kemptide: Free energy perturbation (FEP) calculations

Protein kinases (PKs) discriminate between closely related sequences that contain serine, threonine, and/or tyrosine residues. Such specificity is defined by the amino acid sequence surrounding the phosphorylatable residue, so that it is possible to identify an optimal recognition motif (ORM) for each PK. The ORM for the protein kinase A (PKA), a well-known member of the PK family, is the sequence RRX(S/T)X, where arginines at the -3 and -2 positions play a key role with respect to the primed phosphorylation site. In this work, differential affinities of PKA for the peptide substrate Kemptide (LRRASLG) and mutants that substitute the arginine residues by the unnatural peptide homoarginine were evaluated through molecular dynamics (MD) and free energy perturbation (FEP) calculations. The FEP study for the homoarginine mutants required previous elaboration of a CHARMM "arginine to homoarginine" (R2B) hybrid topology file which is available in this manuscript as Supporting Information. Mutants substituting the arginine residues by alanine, lysine, and histidine were also considered in the comparison by using the same protocol. FEP calculations allowed estimating the free energy changes from the free PKA to PKA-substrate complex (ΔΔG_E→ES ) when Kemptide structure was mutated. Both ΔΔG_S→ES values for homoarginine mutants were predicted with a difference below 1 kcal/mol. In addition, FEP correctly predicted that all the studied mutations decrease the catalytic efficiency of Kemptide for PKA. © 2018 Wiley Periodicals, Inc.

Mechanistic Insights on Human Phosphoglucomutase Revealed by Transition Path Sampling and Molecular Dynamics Calculations

Human α-phosphoglucomutase 1 (α-PGM) catalyzes the isomerization of glucose-1-phosphate into glucose-6-phosphate (G6P) through two sequential phosphoryl transfer steps with a glucose-1,6-bisphosphate (G16P) intermediate. Given that the release of G6P in the gluconeogenesis raises the glucose output levels, α-PGM represents a tempting pharmacological target for type 2 diabetes. Here, we provide the first theoretical study of the catalytic mechanism of human α-PGM. We performed transition-path sampling simulations to unveil the atomic details of the two catalytic chemical steps, which could be key for developing transition state (TS) analogue molecules with inhibitory properties. Our calculations revealed that both steps proceed through a concerted S_N 2-like mechanism, with a loose metaphosphate-like TS. Even though experimental data suggests that the two steps are identical, we observed noticeable differences: 1) the transition state ensemble has a well-defined TS region and a late TS for the second step, and 2) larger coordinated protein motions are required to reach the TS of the second step. We have identified key residues (Arg23, Ser117, His118, Lys389), and the Mg²⁺ ion that contribute in different ways to the reaction coordinate. Accelerated molecular dynamics simulations suggest that the G16P intermediate may reorient without leaving the enzymatic binding pocket, through significant conformational rearrangements of the G16P and of specific loop regions of the human α-PGM.

Fragment molecular orbital study of the cAMP-dependent protein kinase catalyzed phosphoryl transfer: a comparison with the differential transition state stabilization method

The importance of key residues to the activity of the cAMP-dependent protein kinase catalyzed phosphoryl transfer and to the stabilization of the transition state of the reaction has been investigated by means of the fragment molecular orbital (FMO) method. To evaluate the accuracy of the method and its capability of fragmenting covalent bonds, we have compared stabilization energies due to the interactions between individual residues and the reaction center to results obtained with the differential transition state stabilization method (Szarek, et al., J. Phys. Chem. B, 2008, 112, 11819-11826) and observe, despite a size difference in the fragment describing the reaction center, near-quantitative agreement. We have also computed deletion energies to investigate the effect of virtual deletion of key residues on the activation energy. These results are consistent with the stabilization energies and yield additional information as they clearly capture the effect of secondary interactions, i.e. interactions in the second coordination layer of the reaction center. We find that using FMO to calculate deletion energies is a powerful and time efficient approach to analyze the importance of key residues to the activity of an enzyme catalyzed reaction.

Theoretical Study of the Phosphoryl Transfer Reaction from ATP to Dha Catalyzed by DhaK from Escherichia coli

Protein kinases, representing one of the largest protein families involved in almost all aspects of cell life, have become one of the most important targets for the development of new drugs to be used in, for instance, cancer treatments. In this article an exhaustive theoretical study of the phosphoryl transfer reaction from adenosine triphosphate (ATP) to dihydroxyacetone (Dha) catalyzed by DhaK from Escherichia coli (E. coli) is reported. Two different mechanisms, previously proposed for the phosphoryl transfer from ATP to the hydroxyl side chain of specific serine, threonine, or tyrosine residues, have been explored based on the generation of free energy surfaces (FES) computed with hybrid QM/MM potentials. The results suggest that the substrate-assisted phosphoryl and proton-transfer mechanism is kinetically more favorable than the mechanism where an aspartate would be activating the Dha. Although the details of the mechanisms appear to be dramatically dependent on the level of theory employed in the calculations (PM3/MM, B3LYP:PM3/MM, or B3LYP/MM), the transition states (TSs) for the phosphoryl transfer step appear to be described as a concerted step with different degrees of synchronicity in the breaking and forming bonds process in both explored mechanisms. Residues of the active site belonging to different subunits of the protein, such as Gly78B, Thr79A, Ser80A, Arg178B, and one Mg²⁺ cation, would be stabilizing the transferred phosphate in the TS. Asp109A would have a structural role by posing the Dha and other residues of the active site in the proper orientation. The information derived from our calculations not only reveals the role of the enzyme and the particular residues of its active site, but it can assist in the rational design of new more specific inhibitors.

The crystal structure of PknI from Mycobacterium tuberculosis shows an inactive, pseudokinase-like conformation

Eukaryotic-like Ser/Thr protein kinases (ePKs) have been identified in many bacterial species, where they are known to mediate signalling mechanisms that share several features with their eukaryotic counterparts. In Mycobacterium tuberculosis, PknI is one of the 11 predicted ePKs and it has been related to bacterial virulence. In order to better understand the molecular basis of its role in mycobacterial signalling, we solved the crystal structure of the PknI cytoplasmic domain. We found that even though PknI possesses most conserved elements characteristic of Hanks-type kinases, it is degraded in several motifs that are essential for the ePKs catalytic activity. Most notably, PknI presents a remarkably short activation segment lacking a peptide-substrate binding site. Consistent with this observation and similar to earlier findings for eukaryotic pseudokinases, no kinase activity was detected for the catalytic domain of PknI, against different substrates and in various experimental conditions. Based on these results, we conclude that PknI may rely on unconventional mechanism(s) for kinase activity and/or it could play alternative role(s) in mycobacterial signalling.

DATABASE

Atomic coordinates and structure factors for the catalytic domain of M. tuberculosis PknI are in the Protein Data Bank under the accession codes 5M06 (wild-type PknI + ADP), 5M07 (PknI_C20A), 5M08 (PknI_C20A_R136A) and 5M09 (PknI_C20A_R136N).

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.

The influence of active site conformations on the hydride transfer step of the thymidylate synthase reaction mechanism

The hydride transfer from C6 of tetrahydrofolate to the reaction's exocyclic methylene-dUMP intermediate is the rate limiting step in thymidylate synthase (TSase) catalysis. This step has been studied by means of QM/MM molecular dynamics simulations to generate the corresponding free energy surfaces. The use of two different initial X-ray structures has allowed exploring different conformational spaces and the existence of chemical paths with not only different reactivities but also different reaction mechanisms. The results confirm that this chemical conversion takes place preferentially via a concerted mechanism where the hydride transfer is conjugated to thiol-elimination from the product. The findings also confirm the labile character of the substrate-enzyme covalent bond established between the C6 of the nucleotide substrate and a conserved cysteine residue. The calculations also reproduce and rationalize a normal H/T 2° kinetic isotope effect measured for that step. From a computational point of view, the results demonstrate that the use of an incomplete number of coordinates to describe the real reaction coordinate can render biased results.

Phosphoryl transfer reaction catalyzed by membrane diacylglycerol kinase: a theoretical mechanism study

Diacylglycerol kinase is an integral membrane protein which catalyzes phosphoryl transfer from ATP to diacylglycerol. As the smallest kinase known, it shares no sequence homology with conventional kinases and possesses a distinct trimer structure. Thus far, its catalytic mechanism remains elusive. Using molecular dynamics and quantum mechanics calculations, we investigated the co-factor and the substrate binding and phosphoryl transfer mechanism. Based on the analysis of density functional theory calculations, we reveal that the phosphorylation reaction of diacylglycerol kinase features the same phosphoryl transfer mechanism as other kinases, despite its unique structural properties. Our results further show that the active site is relatively open and able to accommodate ligands in multiple orientations, suggesting that the optimization of binding orientations and conformational changes would occur prior to actual phosphoryl transfer.

Dual Role of the Active Site Residues of Thermus thermophilus 3-Isopropylmalate Dehydrogenase: Chemical Catalysis and Domain Closure

The key active site residues K185, Y139, D217, D241, D245, and N102 of Thermus thermophilus 3-isopropylmalate dehydrogenase (Tt-IPMDH) have been replaced, one by one, with Ala. A drastic decrease in the kcat value (0.06% compared to that of the wild-type enzyme) has been observed for the K185A and D241A mutants. Similarly, the catalytic interactions (Km values) of these two mutants with the substrate IPM are weakened by more than 1 order of magnitude. The other mutants retained some (1-13%) of the catalytic activity of the wild-type enzyme and do not exhibit appreciable changes in the substrate Km values. The pH dependence of the wild-type enzyme activity (pK = 7.4) is shifted toward higher values for mutants K185A and D241A (pK values of 8.4 and 8.5, respectively). For the other mutants, smaller changes have been observed. Consequently, K185 and D241 may constitute a proton relay system that can assist in the abstraction of a proton from the OH group of IPM during catalysis. Molecular dynamics simulations provide strong support for the neutral character of K185 in the resting state of the enzyme, which implies that K185 abstracts the proton from the substrate and D241 assists the process via electrostatic interactions with K185. Quantum mechanics/molecular mechanics calculations revealed a significant increase in the activation energy of the hydride transfer of the redox step for both D217A and D241A mutants. Crystal structure analysis of the molecular contacts of the investigated residues in the enzyme-substrate complex revealed their additional importance (in particular that of K185, D217, and D241) in stabilizing the domain-closed active conformation. In accordance with this, small-angle X-ray scattering measurements indicated the complete absence of domain closure in the cases of D217A and D241A mutants, while only partial domain closure could be detected for the other mutants. This suggests that the same residues that are important for catalysis are also essential for inducing domain closure.

Identification and classification of small molecule kinases: insights into substrate recognition and specificity

BACKGROUND

Many prokaryotic kinases that phosphorylate small molecule substrates, such as antibiotics, lipids and sugars, are evolutionarily related to Eukaryotic Protein Kinases (EPKs). These Eukaryotic-Like Kinases (ELKs) share the same overall structural fold as EPKs, but differ in their modes of regulation, substrate recognition and specificity-the sequence and structural determinants of which are poorly understood.

RESULTS

To better understand the basis for ELK specificity, we applied a Bayesian classification procedure designed to identify sequence determinants responsible for functional divergence. This reveals that a large and diverse family of aminoglycoside kinases, characterized members of which are involved in antibiotic resistance, fall into major sub-groups based on differences in putative substrate recognition motifs. Aminoglycoside kinase substrate specificity follows simple rules of alternating hydroxyl and amino groups that is strongly correlated with variations at the DFG + 1 position.

CONCLUSIONS

Substrate specificity determining features in small molecule kinases are mostly confined to the catalytic core and can be identified based on quantitative sequence and crystal structure comparisons.

Reaction Pathway for Cocaine Hydrolase-Catalyzed Hydrolysis of (+)-Cocaine

A recently designed and discovered cocaine hydrolase (CocH), engineered from human butyrylcholinesterase (BChE), has been proven promising as a novel enzyme therapy for treatment of cocaine overdose and addiction because it is highly efficient in catalyzing hydrolysis of naturally occurring (-)-cocaine. It has been known that the CocH also has a high catalytic efficiency against (+)-cocaine, a synthetic enantiomer of cocaine. Reaction pathway and the corresponding free energy profile for the CocH-catalyzed hydrolysis of (+)-cocaine have been determined, in the present study, by performing first-principles pseudobond quantum mechanical/molecular mechanical (QM/MM)-free energy (FE) calculations. Acordingt to the QM/MM-FE results, the catalytic hydrolysis process is initiated by the nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester via hydroxyl oxygen of S198 side chain, and the second reaction step (i.e. dissociation of benzoyl ester) is rate-determining. This finding for CocH-catalyzed hydrolysis of (+)-cocaine is remarkably different from that for the (+)-cocaine hydrolysis catalyzed by bacterial cocaine esterase in which the first reaction step of the deacylation is associated with the highest free energy barrier (~17.9 kcal/mol). The overall free energy barrier (~16.0 kcal/mol) calculated for the acylation stage of CocH-catalyzed hydrolysis of (+)-cocaine is in good agreement with the experimental free energy barrier of ~14.5 kcal/mol derivated from the experimental kinetic data.

Role of Protein Dynamics in Allosteric Control of the Catalytic Phosphoryl Transfer of Insulin Receptor Kinase

The catalytic and allosteric mechanisms of insulin receptor kinase (IRK) are investigated by a combination of ab initio and semiempirical quantum mechanical and molecular mechanical (QM/MM) methods and classical molecular dynamics (MD) simulations. The simulations reveal that the catalytic reaction proceeds in two steps, starting with the transfer of a proton from substrate Tyr to the catalytic Asp1132, followed by the phosphoryl transfer from ATP to substrate Tyr. The enhancement of the catalytic rate of IRK upon phosphorylations in the enzyme's activation loop is found to occur mainly via changes to the free energy landscape of the proton transfer step, favoring the proton transfer in the fully phosphorylated enzyme. In contrast, the effects of the phosphorylations on the phosphoryl transfer are smaller. Equilibrium MD simulations show that IRK phosphorylations affect the protein dynamics of the enzyme before the proton transfer to Asp1132 with only a minor effect after the proton transfer. This finding is consistent with the large change in the proton transfer free energy and the smaller change in the free energy barrier of phosphoryl transfer found by QM/MM simulations. Taken together, the present results provide details on how IRK phosphorylation exerts allosteric control of the catalytic activity via modifications of protein dynamics and free energy landscape of catalytic reaction. The results also highlight the importance of protein dynamics in connecting protein allostery and catalysis to control catalytic activity of enzymes.

Characterization of the Catalytic and Nucleotide Binding Properties of the α-Kinase Domain of Dictyostelium Myosin-II Heavy Chain Kinase A

The α-kinases are a widely expressed family of serine/threonine protein kinases that exhibit no sequence identity with conventional eukaryotic protein kinases. In this report, we provide new information on the catalytic properties of the α-kinase domain of Dictyostelium myosin-II heavy chain kinase-A (termed A-CAT). Crystallization of A-CAT in the presence of MgATP yielded structures with AMP or adenosine in the catalytic cleft together with a phosphorylated Asp-766 residue. The results show that the β- and α-phosphoryl groups are transferred either directly or indirectly to the catalytically essential Asp-766. Biochemical assays confirmed that A-CAT hydrolyzed ATP, ADP, and AMP with kcat values of 1.9, 0.6, and 0.32 min(-1), respectively, and showed that A-CAT can use ADP to phosphorylate peptides and proteins. Binding assays using fluorescent 2'/3'-O-(N-methylanthraniloyl) analogs of ATP and ADP yielded Kd values for ATP, ADP, AMP, and adenosine of 20 ± 3, 60 ± 20, 160 ± 60, and 45 ± 15 μM, respectively. Site-directed mutagenesis showed that Glu-713, Leu-716, and Lys-645, all of which interact with the adenine base, were critical for nucleotide binding. Mutation of the highly conserved Gln-758, which chelates a nucleotide-associated Mg(2+) ion, eliminated catalytic activity, whereas loss of the highly conserved Lys-722 and Arg-592 decreased kcat values for kinase and ATPase activities by 3-6-fold. Mutation of Asp-663 impaired kinase activity to a much greater extent than ATPase, indicating a specific role in peptide substrate binding, whereas mutation of Gln-768 doubled ATPase activity, suggesting that it may act to exclude water from the active site.

Catalytic mechanisms and regulation of protein kinases

Protein kinases transfer a phosphoryl group from ATP onto target proteins and play a critical role in signal transduction and other cellular processes. Here, we review the kinase kinetic and chemical mechanisms and their application in understanding kinase structure and function. Aberrant kinase activity has been implicated in many human diseases, in particular cancer. We highlight applications of technologies and concepts derived from kinase mechanistic studies that have helped illuminate how kinases are regulated and contribute to pathophysiology.

Determining the Functions of HIV-1 Tat and a Second Magnesium Ion in the CDK9/Cyclin T1 Complex: A Molecular Dynamics Simulation Study

The current paradigm of cyclin-dependent kinase (CDK) regulation based on the well-established CDK2 has been recently expanded. The determination of CDK9 crystal structures suggests the requirement of an additional regulatory protein, such as human immunodeficiency virus type 1 (HIV-1) Tat, to exert its physiological functions. In most kinases, the exact number and roles of the cofactor metal ions remain unappreciated, and the repertoire has thus gained increasing attention recently. Here, molecular dynamics (MD) simulations were implemented on CDK9 to explore the functional roles of HIV-1 Tat and the second Mg²⁺ ion at site 1 (Mg₁²⁺). The simulations unveiled that binding of HIV-1 Tat to CDK9 not only stabilized hydrogen bonds (H-bonds) between ATP and hinge residues Asp104 and Cys106, as well as between ATP and invariant Lys48, but also facilitated the salt bridge network pertaining to the phosphorylated Thr186 at the activation loop. By contrast, these H-bonds cannot be formed in CDK9 owing to the absence of HIV-1 Tat. MD simulations further revealed that the Mg₁²⁺ ion, coupled with the Mg₂²⁺ ion, anchored to the triphosphate moiety of ATP in its catalytic competent conformation. This observation indicates the requirement of the Mg₁²⁺ ion for CDK9 to realize its function. Overall, the introduction of HIV-1 Tat and Mg₁²⁺ ion resulted in the active site architectural characteristics of phosphorylated CDK9. These data highlighted the functional roles of HIV-1 Tat and Mg₁²⁺ ion in the regulation of CDK9 activity, which contributes an important complementary understanding of CDK molecular underpinnings.

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.

Study of the affinity between the protein kinase PKA and peptide substrates derived from kemptide using molecular dynamics simulations and MM/GBSA

We have carried out a protocol in computational biochemistry including molecular dynamics (MD) simulations and MM/GBSA free energy calculations on the complex between the protein kinase A (PKA) and the specific peptide substrate Kemptide (LRRASLG). We made the same calculations on other PKA complexes that contain Kemptide derivatives (with mutations of the arginines, and with deletions of N and C-terminal amino acids). We predicted shifts in the free energy changes from the free PKA to PKA-substrate complex (ΔΔG_E→ES) when Kemptide structure is modified (we consider that the calculated shifts correlate with the experimental shifts of the free energy changes from the free PKA to the transition states (ΔΔG_E→TS) determined by the catalytic efficiency (k_cat/K_M) changes). Our results demonstrate that it is possible to predict the kinetic properties of protein kinases using simple computational biochemistry methods. As an additional benefit, these methods give detailed molecular information that permit the analysis of the atomic forces that contribute to the affinity between protein kinases and their substrates.

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.

ErbB/HER protein-tyrosine kinases: Structures and small molecule inhibitors

The epidermal growth factor receptor (EGFR) family consists of four members that belong to the ErbB lineage of proteins (ErbB1-4). These receptors consist of an extracellular domain, a single hydrophobic transmembrane segment, and an intracellular portion with a juxtamembrane segment, a protein kinase domain, and a carboxyterminal tail. The ErbB proteins function as homo and heterodimers. Growth factor binding to EGFR induces a large conformational change in the extracellular domain. Two ligand-EGFR complexes unite to form a back-to-back dimer in which the ligands are on opposite sides of the aggregate. Following ligand binding, EGFR intracellular kinase domains form an asymmetric dimer. The carboxyterminal lobe of the activator kinase of the dimer interacts with the amino-terminal lobe of the receiver kinase thereby leading to its allosteric stimulation. Several malignancies are associated with the mutation or increased expression of members of the ErbB family including lung, breast, stomach, colorectal, head and neck, and pancreatic carcinomas. Gefitinib, erlotinib, and afatinib are orally effective protein-kinase targeted quinazoline derivatives that are used in the treatment of ERBB1-mutant lung cancer and lapatinib is an orally effective quinazoline derivative used in the treatment of ErbB2-overexpressing breast cancer. Moreover, monoclonal antibodies that target the extracellular domain of ErbB2 are used for the treatment of ErbB2-positive breast cancer and monoclonal antibodies that target ErbB1 and are used for the treatment of colorectal cancer. Cancers treated with these targeted drugs eventually become resistant to them, and a current goal of research is to develop drugs that are effective against drug-resistant tumors.

Metal-free cAMP-dependent protein kinase can catalyze phosphoryl transfer

X-ray structures of several ternary product complexes of the catalytic subunit of cAMP-dependent protein kinase (PKAc) have been determined with no bound metal ions and with Na(+) or K(+) coordinated at two metal-binding sites. The metal-free PKAc and the enzyme with alkali metals were able to facilitate the phosphoryl transfer reaction. In all studied complexes, the ATP and the substrate peptide (SP20) were modified into the products ADP and the phosphorylated peptide. The products of the phosphotransfer reaction were also found when ATP-γS, a nonhydrolyzable ATP analogue, reacted with SP20 in the PKAc active site containing no metals. Single turnover enzyme kinetics measurements utilizing (32)P-labeled ATP confirmed the phosphotransferase activity of the enzyme in the absence of metal ions and in the presence of alkali metals. In addition, the structure of the apo-PKAc binary complex with SP20 suggests that the sequence of binding events may become ordered in a metal-free environment, with SP20 binding first to prime the enzyme for subsequent ATP binding. Comparison of these structures reveals conformational and hydrogen bonding changes that might be important for the mechanism of catalysis.

Computational delineation of tyrosyl-substrate recognition and catalytic landscapes by the epidermal growth factor receptor tyrosine kinase domain

The epidermal growth factor receptor (EGFR) is a receptor tyrosine kinase (RTK), which catalyzes protein phosphorylation reactions by transferring the γ-phosphoryl group from an ATP molecule to the hydroxyl group of tyrosine residues in protein substrates. EGFR is an important drug target in the treatment of cancers and a better understanding of the receptor function is critical to discern cancer mechanisms. We employ a suite of molecular simulation methods to explore the mechanism of substrate recognition and to delineate the catalytic landscape of the phosphoryl transfer reaction. Based on our results, we propose that a highly conserved region corresponding to Val852-Pro853-Ile854-Lys855-Trp856 in the EGFR tyrosine kinase domain (TKD) is essential for substrate binding. We also provide a possible explanation for the established experimental observation that protein tyrosine kinases (including EGFR) select substrates with a glutamic acid at the P - 1 position and a large hydrophobic amino acid at the P + 1 position. Furthermore, our mixed quantum mechanics/molecular mechanics (QM/MM) simulations show that the EGFR protein kinase favors the dissociative mechanism, although an alternative channel through the formation of an associative transition state is also possible. Our simulations establish some key molecular rules in the operation for substrate-recognition and for phosphoryl transfer in the EGFR TKD.

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.

Molecular dynamics simulation studies on the positive cooperativity of the Kemptide substrate with protein kinase A induced by the ATP ligand

The positive cooperativity of the Kemptide substrate or the ATP molecule with the PKA catalytic subunit has been studied by dynamics simulations and free energy calculations on a series of binary and ternary models. The results revealed that the first ATP binding to the PKA catalytic subunit is energetically favorable for the successive Kemptide binding, confirming the positive cooperativity. The key residues Thr51, Glu170, and Phe187 in PKA contributing to the positive cooperativity have been found. The binding of ATP to PKA induces the positive cooperativity through one direct allosteric communication network in PKA from the ATP binding sites in the catalytic loop of the large lobe to the Kemptide binding sites in the activation segment of the large lobe, two indirect ones from those in the glycine-rich loop and the β3 strand of the small lobe, and from those in the catalytic loop to those in the activation segment via the αF helix media. The Tyr204Ala mutation in the activation segment of PKA causes both the decoupling of the cooperativity and the disruption of the corresponding allosteric network through the αF helix media.

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.

Insights into the phosphoryl transfer catalyzed by cAMP-dependent protein kinase: an X-ray crystallographic study of complexes with various metals and peptide substrate SP20

X-ray structures of several ternary substrate and product complexes of the catalytic subunit of cAMP-dependent protein kinase (PKAc) have been determined with different bound metal ions. In the PKAc complexes, Mg²⁺, Ca²⁺, Sr²⁺, and Ba²⁺ metal ions could bind to the active site and facilitate the phosphoryl transfer reaction. ATP and a substrate peptide (SP20) were modified, and the reaction products ADP and the phosphorylated peptide were found trapped in the enzyme active site. Finally, we determined the structure of a pseudo-Michaelis complex containing Mg²⁺, nonhydrolyzable AMP-PCP (β,γ-methyleneadenosine 5′-triphosphate) and SP20. The product structures together with the pseudo-Michaelis complex provide snapshots of different stages of the phosphorylation reaction. Comparison of these structures reveals conformational, coordination, and hydrogen bonding changes that might occur during the reaction and shed new light on its mechanism, roles of metals, and active site residues.

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.

How calcium inhibits the magnesium-dependent kinase gsk3β: a molecular simulation study

Glycogen synthase kinase 3β (GSK3β) is a ubiquitous serine/threonine kinase that plays a pivotal role in many biological processes. GSK3β catalyzes the transfer of γ-phosphate of ATP to the unique substrate Ser/Thr residues with the assistance of two natural activating cofactors Mg(2+). Interestingly, the biological observation reveals that a non-native Ca(2+) ion can inhibit the GSK3β catalytic activity. Here, the inhibitory mechanism of GSK3β by the displacement of native Mg(2+) at site 1 by Ca(2+) was investigated by means of 80 ns comparative molecular dynamics (MD) simulations of the GSK3β···Mg(2+)-2/ATP/Mg(2+) -1 and GSK3β···Mg(2+)-2/ATP/Ca(2+)-1 systems. MD simulation results revealed that using the AMBER point charge model force field for Mg(2+) was more appropriate in the reproduction of the active site architectural characteristics of GSK3β than using the magnesium-cationic dummy atom model force field. Compared with the native Mg(2+) bound system, the misalignment of the critical triphosphate moiety of ATP, the erroneous coordination environments around the Mg(2+) ion at site 2, and the rupture of the key hydrogen bond between the invariant Lys85 and the ATP O(β2) atom in the Ca(2+) substituted system were observed in the MD simulation due to the Ca(2+) ion in active site in order to achieve its preferred sevenfold coordination geometry, which adequately abolish the enzymatic activity. The obtained results are valuable in understanding the possible mechanism by why Ca(2+) inhibits the GSK3β activity and also provide insights into the mechanism of Ca(2+) inhibition in other structurally related protein kinases.