Computational Study of the Phosphoryl Transfer Catalyzed by a Cyclin-Dependent Kinase

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.

Structural insights and influence of V599 mutations on the overall dynamics of BRAF protein against its kinase domains

Mutations in the BRAF gene are well known for their oncogenic effects. Point mutations in V599 are particularly oncogenic and are considered important for therapeutic purposes. Along with wild type, other V599 mutated BRAF variants viz. V599E, V599D and V599R are reported and crystals of the former two with inhibitor (BAY43-9006) are further detailed. Both wild-type and mutated BRAF forms show similar interaction patterns with BAY43-9006, but the 599th residue did not show any involvement in the interactions. Upon BAY43-9006 binding, kinase domains of both forms were found adopting essentially identical conformations. However, BAY43-9006 shows a varied activity profile in the case of the wild and V599E variant of the BRAF protein. Furthermore, MMGBSA binding energy results for all four BRAF variants, further revealed the importance of the 599th residue. In-depth analysis viz. molecular dynamics, residue correlation studies and residue interaction network (RIN) analyses were conducted, providing a deep insight into the 599th residue and its impact on the overall dynamics of BRAF protein. Our findings reveal that the mutated residue at the 599th position not only changed the BAY43-9006-BRAF binding behaviour but also produced a massive impact on the overall dynamic behaviour of the protein. The insights obtained herein could be of great relevance for designing new BRAF inhibitors aimed at getting ideal activity against all BRAF forms.

A Transient and Flexible Cation-π Interaction Promotes Hydrolysis of Nucleic Acids in DNA and RNA Nucleases

Metal-dependent DNA and RNA nucleases are enzymes that cleave nucleic acids with great efficiency and precision. These enzyme-mediated hydrolytic reactions are fundamental for the replication, repair, and storage of genetic information within the cell. Here, extensive classical and quantum-based free-energy molecular simulations show that a cation-π interaction is transiently formed in situ at the metal core of Bacteriophage-λ Exonuclease (Exo-λ), during catalysis. This noncovalent interaction (Lys131-Tyr154) triggers nucleophile activation for nucleotide excision. Then, our simulations also show the oscillatory dynamics and swinging of the newly formed cation-π dyad, whose conformational change may favor proton release from the cationic Lys131 to the bulk solution, thus restoring the precatalytic protonation state in Exo-λ. Altogether, we report on the novel mechanistic character of cation-π interactions for catalysis. Structural and bioinformatic analyses support that flexible orientation and transient formation of mobile cation-π interactions may represent a common catalytic strategy to promote nucleic acid hydrolysis in DNA and RNA nucleases.

Are peptides a solution for the treatment of hyperactivated JAK3 pathways?

While the inactivation mutations that eliminate JAK3 function lead to the immunological disorders such as severe combined immunodeficiency, activation mutations, causing constitutive JAK3 signaling, are known to trigger various types of cancer or are responsible for autoimmune diseases, such as rheumatoid arthritis, psoriasis, or inflammatory bowel diseases. Treatment of hyperactivated JAK3 is still an obstacle, due to different sensibility of mutation types to conventional drugs and unwanted side effects, because these drugs are not absolutely specific for JAK3, thus inhibiting other members of the JAK family, too. Lack of information, in which way sole inhibition of JAK3 is necessary for elimination of the disease, calls for the development of isoform-specific JAK3 inhibitors. Beside this strategy, up to date peptides are a rising alternative as chemo- or immunotherapeutics, but still sparsely represented in drug development and clinical trials. Beyond a possible direct inhibition function, crossing the cancer cell membrane and interfering in disease-causing pathways or triggering apoptosis, peptides could be used in future as adjunct remedies to potentialize traditional therapy and preserve non-affected cells. To discuss such feasible topics, this review deals with the knowledge about the structure-function of JAK3 and the actual state-of-the-art of isoform-specific inhibitor development, as well as the function of currently approved drugs or those currently being tested in clinical trials. Furthermore, several strategies for the application of peptide-based drugs for cancer therapy and the physicochemical and structural relations to peptide efficacy are discussed, and an overview of peptide sequences, which were qualified for clinical trials, is given.

Quantum Mechanical/Molecular Mechanical Analysis of the Catalytic Mechanism of Phosphoserine Phosphatase

Phosphoserine phosphatase (PSP), a member of the haloacid dehalogenase (HAD) superfamily that comprises the vast majority of phosphotransferases, is likely a steady-state regulator of the level of d-serine in the brain. The proposed catalytic cycle of PSP consists of a two-step mechanism: formation of a phospho-enzyme intermediate by phosphate transfer to Asp11 and its subsequent hydrolysis. Our combined quantum mechanical/molecular mechanical (QM/MM) calculations of the reaction pathways favour a dissociative mechanism of nucleophilic substitution via a trigonal-planar metaphosphate-like configuration for both steps, associated with proton transfer to the leaving group or from the nucleophile. This proton transfer is facilitated by active site residue Asp13 that acts as both a general base and a general acid. Free energy calculation on the reaction pathways further support the structural role of the enzymatic environment and the active site architecture. The choice of a proper reaction coordinate along which to bias the free energy calculations can be guided by a projection of the canonical reaction coordinate obtained from a chain-of-state optimisation onto important internal coordinates.

Reduced state transition barrier of CDK6 from open to closed state induced by Thr177 phosphorylation and its implication in binding modes of inhibitors

BACKGROUND

CDK6 is considered as a highly validated anticancer drug target due to its essential role in regulating cell cycle progression at G1 restriction point. Activation of CDK6 requires the phosphorylation of Thr177 on A-loop, but the structural insights of the activation mechanism remain unclear.

METHODS

Herein, all-atoms molecular dynamics (MD) simulations were used to study the effects of Thr177 phosphorylation on the dynamic structure of CDK6-Vcyclin complex.

RESULTS

MD results indicated that the free energy barrier of the transition from open to closed state decreased ~47.2% after Thr177 phosphorylation. Key steps along the state transition process were obtained from a cluster analysis. Binding preference of ten different inhibitors to open or closed state were also investigated through molecular docking along with MD simulations methods.

CONCLUSIONS

Our results indicated that Thr177 phosphorylation increased the flexibility around the ATP-binding pocket. The transition of the ATP-binding pocket between open and closed states should be considered for understanding the binding of CDK6 inhibitors.

GENERAL SIGNIFICANCE

This work could deepen the understanding of CDKs activation mechanism, and provide useful information for the discovery of new CDKs inhibitors with high affinity and specificity.

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.

Unraveling the Catalytic Pathway of Metalloenzyme Farnesyltransferase through QM/MM Computation

The protein farnesyltransferase (FTase) is a Zn(2+)-metalloenzyme that catalyzes the farnesylation reaction, i.e., the transfer of the 15-carbon atom farnesyl group from farnesyl diphosphate (FPP) to a specific cysteine of protein substrates. Oncogenic Ras proteins, which are among the FTase substrates, are observed in about 20-30% of human cancer cells. Thus, FTase represents a target for anticancer drug design. Herein, we present a classical force-field-based and quantum mechanics/molecular mechanics (QM/MM) computational study of the FTase reaction mechanism. Our findings offer a detailed picture of the FTase catalytic pathway, describing structural features and the energetics of its saddle points. A moderate dissociation of the diphosphate group from the FPP is observed during the nucleophilic attack of the zinc-bound thiolate. At the transition state, a resonance structure is observed, which indicates the formation of a metastable carbocation. However, no stable intermediate is found along the reaction pathway. Thus, the reaction occurs via an associative mechanism with dissociative character, in agreement with the mechanism proposed by Fierke et al. ( Biochemistry 2000, 39, 2593-2602 and Biochemistry 2003, 42, 9741-9748 ). Moreover, a fluorine-substituted FPP analogue (CF3-FPP) is used to investigate the inhibitory effect of fluorine, which in turn provides additional agreement with experimental data.

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.

Quantum mechanical modeling: a tool for the understanding of enzyme reactions

Most enzyme reactions involve formation and cleavage of covalent bonds, while electrostatic effects, as well as dynamics of the active site and surrounding protein regions, may also be crucial. Accordingly, special computational methods are needed to provide an adequate description, which combine quantum mechanics for the reactive region with molecular mechanics and molecular dynamics describing the environment and dynamic effects, respectively. In this review we intend to give an overview to non-specialists on various enzyme models as well as established computational methods and describe applications to some specific cases. For the treatment of various enzyme mechanisms, special approaches are often needed to obtain results, which adequately refer to experimental data. As a result of the spectacular progress in the last two decades, most enzyme reactions can be quite precisely treated by various computational methods.

Why nature really chose phosphate

Phosphoryl transfer plays key roles in signaling, energy transduction, protein synthesis, and maintaining the integrity of the genetic material. On the surface, it would appear to be a simple nucleophile displacement reaction. However, this simplicity is deceptive, as, even in aqueous solution, the low-lying d-orbitals on the phosphorus atom allow for eight distinct mechanistic possibilities, before even introducing the complexities of the enzyme catalyzed reactions. To further complicate matters, while powerful, traditional experimental techniques such as the use of linear free-energy relationships (LFER) or measuring isotope effects cannot make unique distinctions between different potential mechanisms. A quarter of a century has passed since Westheimer wrote his seminal review, 'Why Nature Chose Phosphate' (Science 235 (1987), 1173), and a lot has changed in the field since then. The present review revisits this biologically crucial issue, exploring both relevant enzymatic systems as well as the corresponding chemistry in aqueous solution, and demonstrating that the only way key questions in this field are likely to be resolved is through careful theoretical studies (which of course should be able to reproduce all relevant experimental data). Finally, we demonstrate that the reason that nature really chose phosphate is due to interplay between two counteracting effects: on the one hand, phosphates are negatively charged and the resulting charge-charge repulsion with the attacking nucleophile contributes to the very high barrier for hydrolysis, making phosphate esters among the most inert compounds known. However, biology is not only about reducing the barrier to unfavorable chemical reactions. That is, the same charge-charge repulsion that makes phosphate ester hydrolysis so unfavorable also makes it possible to regulate, by exploiting the electrostatics. This means that phosphate ester hydrolysis can not only be turned on, but also be turned off, by fine tuning the electrostatic environment and the present review demonstrates numerous examples where this is the case. Without this capacity for regulation, it would be impossible to have for instance a signaling or metabolic cascade, where the action of each participant is determined by the fine-tuned activity of the previous piece in the production line. This makes phosphate esters the ideal compounds to facilitate life as we know it.

The increasing role of QM/MM in drug discovery

Since its first appearance in 1976, the quantum mechanics/molecular mechanics (QM/MM) approach has mostly been used to study the chemical reactions of enzymes, which are frequently the target of drug discovery programs. In principle, a detailed understanding of the enzymatic mechanism should help researchers to design a potent enzyme inhibitor or new drug. However, QM/MM has not yet had a widespread impact on structure-based drug design. This is mostly due to its high computational cost. We expect this to change with the recent and extraordinary increases in computational power, and with the availability of more efficient algorithms for QM/MM calculations. Here, we report on some representative examples of QM/MM studies, including our own research, of pharmaceutically relevant enzymes, such as ribonuclease H and fatty acid amide hydrolase (FAAH). We aim to show how QM/MM has traditionally been used to study enzymatic catalysis. In this regard, we discuss its potential to become a routinely used drug design tool. To support this, we also discuss selected computational studies where QM/MM insights have been helpful in improving the potency of covalent inhibitors of FAAH.

Cyclin-dependent kinases: bridging their structure and function through computations

Cyclin-dependent kinases (CDKs) are one of the most promising target families for drug discovery for several diseases, such as cancer and neurodegenerative disorders. Over the years, structural insights on CDKs have demonstrated high protein plasticity, with several cases where two or more structures of the same protein adopt different conformations. This has generated a great deal of interest in understanding the relationship between CDK structure and function. Here, we highlight how computer simulations have recently contributed in characterizing some key rare and transient events in CDKs, such as the reaction transition state and activation loop movement. Although not yet fully defined, we can now portray the enzymatic mechanism and plasticity of CDKs at high spatial and temporal resolution. These theoretical studies bridge with experiments and highlight structural determinants that could help in designing specific CDK inhibitors.

Insights into the phosphoryl transfer mechanism of cyclin-dependent protein kinases from ab initio QM/MM free-energy studies

Phosphorylation reactions catalyzed by kinases and phosphatases play an indispensible role in cellular signaling, and their malfunctioning is implicated in many diseases. A better understanding of the catalytic mechanism will help design novel and effective mechanism-based inhibitors of these enzymes. In this work, ab initio quantum mechanical/molecular mechanical studies are reported for the phosphoryl transfer reaction catalyzed by a cyclin-dependent kinase, CDK2. Our results suggest that an active-site Asp residue, rather than ATP as previously proposed, serves as the general base to activate the Ser nucleophile. The corresponding transition state features a dissociative, metaphosphate-like structure, stabilized by the Mg(2+) ion and several hydrogen bonds. The calculated free-energy barrier is consistent with experimental values. Implications of our results in this and other protein kinases are discussed.

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

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

[Modeling the transition state of enzyme-catalyzed phosphoryl transfer reaction using QM/MM method]

Reversible phosphorylation of proteins is a post-translational modification that regulates diverse biological processes. The molecular mechanism underlying phosphoryl transfer catalyzed by enzymes, in particular the nature of transition state (TS), remains a subject of active debate. Structural evidence supports an associative TS, whereas physical organic studies point to a dissociative character. In this article, we briefly introduce our recent effort using the hybrid quantum mechanics/molecular mechanics (QM/MM) simulations to resolve the controversy. We perform QM/MM simulations for the reversible phosphorylation of phosphoserine phosphatase (PSP), which belongs to one of the largest phosphotransferase families characterized to data. Both phosphorylation and dephosphorylation reactions are investigated based on the two-dimensional energy surfaces along phosphoryl and proton transfer coordinates. The resultant structures of the active site at TS in both reactions have compact geometries but a less electron density of the phosphoryl group. This suggests that the TS of PSP has a geometrically associative yet electronically dissociative character and strongly depends on proton transfer being coupled with phosphoryl transfer. Structure and literature database searches on phosphotransferases suggest that such a hybrid TS is consistent with many structures and physical organic studies and likely holds for most enzymes catalyzing phosphoryl transfer.

Molecular modeling and molecular dynamics simulation studies of the GSK3β/ATP/substrate complex: understanding the unique P+4 primed phosphorylation specificity for GSK3β substrates

Substrate specificity of protein kinases is of fundamental importance for the integrity and fidelity of signaling pathways. Glycogen synthase kinase 3β (GSK3β) has a unique substrate specificity that prefers phosphorylation of its substrates at the P+4 serine before it can further phosphorylate the substrate at the P0 serine in the canonical motif SXXXS(p), where S(p) is the primed phosphorylation site. The detailed phosphorylation mechanism, however, is not clearly understood. In this study, a three-dimensional (3D) model of the ternary complex of GSK3β, ATP, and the phosphorylated glycogen synthase (pGS), termed GSK3β/ATP/pGS, is constructed using a hierarchical approach and by integrating molecular modeling and molecular dynamics (MD) simulations. Based on the 3D model, the substrate primed phosphorylation mechanism is investigated via two 12 ns comparative MD simulations of the GSK3β/ATP/pGS and GSK3β/ATP/GS systems, which differ in the phosphate group bound to the P+4 serine of GS. In agreement with structural analysis, computed binding free energies reveal that the binding of pGS to GSK3β is favored in the prephosphorylated state compared with the GS native state. More importantly, comparison with the system simulated without primed phosphorylation in the GSK3β/ATP/GS complex shows that for an optimal phosphorylation reaction to occur, the pGS priming phosphate in the GSK3β/ATP/pGS system optimizes the proper orientation of the GSK3β N- and C-terminal domains and clamps the P0 serine of pGS in the appropriate configuration for interaction with the ATP γ-phosphate within the catalytic groove.

Quantum mechanics/molecular mechanics investigation of the mechanism of phosphate transfer in human uridine-cytidine kinase 2

The mechanisms of enzyme-catalyzed phosphate transfer and hydrolysis reactions are of great interest due to their importance and abundance in biochemistry. The reaction may proceed in a stepwise fashion, with either a pentavalent phosphorane or a metaphosphate anion intermediate, or by a concerted SN2 mechanism. Despite much theoretical work focused on a few key enzymes, a consensus for the mechanism has not been reached, and examples of all three possibilities have been demonstrated. We have investigated the mechanism of human uridine-cytidine kinase 2 (UCK2, EC 2.7.1.48), which catalyzes the transfer of a phosphate group from ATP to the ribose 5'-hydroxyl of cytidine and uridine. UCK2 is normally expressed in human placenta, but is overexpressed in certain cancer cells, where it is responsible for activating a class of antitumor prodrugs. The UCK2 mechanism was investigated by generating a 2D potential energy surface as a function of the P-O bonds forming and breaking, with energies calculated using a quantum mechanics/molecular mechanics potential (B3LYP/6-31G(d):AMBER). The mechanism of phosphate transfer is shown to be concerted, and is accompanied by concerted proton transfer from the 5'-hydroxyl to a conserved active site aspartic acid that serves as a catalytic base. The calculated barrier for this reaction is 15.1 kcal/mol, in relatively good agreement with the experimental barrier of 17.5 kcal/mol. The interactions of the enzyme active site with the reactant, transition state, and product are examined for their implications on the design of anticancer prodrugs or positron emission tomography (PET) reporter probes for this enzyme.

How mitogen-activated protein kinases recognize and phosphorylate their targets: A QM/MM study

Mitogen-activated protein kinase (MAPK) signaling pathways play an essential role in the transduction of environmental stimuli to the nucleus, thereby regulating a variety of cellular processes, including cell proliferation, differentiation, and programmed cell death. The components of the MAPK extracellular activated protein kinase (ERK) cascade represent attractive targets for cancer therapy, as their aberrant activation is a frequent event among highly prevalent human cancers. To understand how MAPKs recognize and phosphorylate their targets is key to unravel their function. However, these events are still poorly understood because of the lack of complex structures of MAPKs with their bound targets in the active site. Here we have modeled the interaction of ERK with a target peptide and analyzed the specificity toward Ser/Thr-Pro motifs. By using a quantum mechanics/molecular mechanics (QM/MM) approach, we propose a mechanism for the phosphoryl transfer catalyzed by ERK that offers new insights into MAPK function. Our results suggest that (1) the proline residue has a role in both specificity and phospho transfer efficiency, (2) the reaction occurs in one step, with ERK2 Asp(147) acting as the catalytic base, (3) a conserved Lys in the kinase superfamily that is usually mutated to check kinase activity strongly stabilizes the transition state, and (4) the reaction mechanism is similar with either one or two Mg(2+) ions in the active site. Taken together, our results provide a detailed description of the molecular events involved in the phosphorylation reaction catalyzed by MAPK and contribute to the general understanding of kinase activity.

QM/MM methods for biomolecular systems

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

Cell cycle kinases predicted from conserved biophysical properties

Machine-learning techniques can classify functionally related proteins where homology-transfer as well as sequence and structure motifs fail. Here, we present a method that aimed at complementing homology-transfer in the identification of cell cycle control kinases from sequence alone. First, we identified functionally significant residues in cell cycle proteins through their high sequence conservation and biophysical properties. We then incorporated these residues and their features into support vector machines (SVM) to identify new kinases and more specifically to differentiate cell cycle kinases from other kinases and other proteins. As expected, the most informative residues tend to be highly conserved and tend to localize in the ATP binding regions of the kinases. Another observation confirmed that ATP binding regions are typically not found on the surface but in partially buried sites, and that this fact is correctly captured by accessibility predictions. Using these highly conserved, semi-buried residues and their biophysical properties, we could distinguish cell cycle S/T kinases from other kinase families at levels around 70-80% accuracy and 62-81% coverage. An application to the entire human proteome predicted at least 97 human proteins with limited previous annotations to be candidates for cell cycle kinases.

Phosphodiester cleavage in ribonuclease H occurs via an associative two-metal-aided catalytic mechanism

Ribonuclease H (RNase H) belongs to the nucleotidyl-transferase (NT) superfamily and hydrolyzes the phosphodiester linkages that form the backbone of the RNA strand in RNA x DNA hybrids. This enzyme is implicated in replication initiation and DNA topology restoration and represents a very promising target for anti-HIV drug design. Structural information has been provided by high-resolution crystal structures of the complex RNase H/RNA x DNA from Bacillus halodurans (Bh), which reveals that two metal ions are required for formation of a catalytic active complex. Here, we use classical force field-based and quantum mechanics/molecular mechanics calculations for modeling the nucleotidyl transfer reaction in RNase H, clarifying the role of the metal ions and the nature of the nucleophile (water versus hydroxide ion). During the catalysis, the two metal ions act cooperatively, facilitating nucleophile formation and stabilizing both transition state and leaving group. Importantly, the two Mg(2+) metals also support the formation of a meta-stable phosphorane intermediate along the reaction, which resembles the phosphorane intermediate structure obtained only in the debated beta-phosphoglucomutase crystal (Lahiri, S. D.; et al. Science 2003, 299 (5615), 2067-2071). The nucleophile formation (i.e., water deprotonation) can be achieved in situ, after migration of one proton from the water to the scissile phosphate in the transition state. This proton transfer is actually mediated by solvation water molecules. Due to the highly conserved nature of the enzymatic bimetal motif, these results might also be relevant for structurally similar enzymes belonging to the NT superfamily.

The associative nature of adenylyl transfer catalyzed by T4 DNA ligase

DNA ligase seals nicks in dsDNA using chemical energy of the phosphoanhydride bond in ATP or NAD(+) and assistance of a divalent metal cofactor Mg(2+). Molecular details of ligase catalysis are essential for understanding the mechanism of metal-promoted phosphoryl transfer reactions in the living cell responsible for a wide range of processes, e.g., DNA replication and transcription, signaling and differentiation, energy coupling and metabolism. Here we report a single-turnover (31)P solid-state NMR study of adenylyl transfer catalyzed by DNA ligase from bacteriophage T4. Formation of a high-energy covalent ligase-nucleotide complex is triggered in situ by the photo release of caged Mg(2+), and sequentially formed intermediates are monitored by NMR. Analyses of reaction kinetics and chemical-shift changes indicate that the pentacoordinated phosphorane intermediate builds up to 35% of the total reacting species after 4-5 h of reaction. This is direct experimental evidence of the associative nature of adenylyl transfer catalyzed by DNA ligase. NMR spectroscopy in rotating solids is introduced as an analytical tool for recording molecular movies of reaction processes. Presented work pioneers a promising direction in structural studies of biochemical transformations.

Regulatory phosphorylation of cyclin-dependent kinase 2: insights from molecular dynamics simulations

The structures of fully active cyclin-dependent kinase-2 (CDK2) complexed with ATP and peptide substrate, CDK2 after the catalytic reaction, and CDK2 inhibited by phosphorylation at Thr14/Tyr15 were studied using molecular dynamics (MD) simulations. The structural details of the CDK2 catalytic site and CDK2 substrate binding box were described. Comparison of MD simulations of inhibited complexes of CDK2 was used to help understand the role of inhibitory phosphorylation at Thr14/Tyr15. Phosphorylation at Thr14/Tyr15 causes ATP misalignment for the phosphate-group transfer, changes in the Mg(2+) coordination sphere, and changes in the H-bond network formed by CDK2 catalytic residues (Asp127, Lys129, Asn132). The inhibitory phosphorylation causes the G-loop to shift from the ATP binding site, which leads to opening of the CDK2 substrate binding box, thus probably weakening substrate binding. All these effects explain the decrease in kinase activity observed after inhibitory phosphorylation at Thr14/Tyr15 in the G-loop. Interaction of the peptide substrate, and the phosphorylated peptide product, with CDK2 was also studied and compared. These results broaden hypotheses drawn from our previous MD studies as to why a basic residue (Arg/Lys) is preferred at the P(+2) substrate position.