Molecular Insights into the Regulation of 3-Phosphoinositide-Dependent Protein Kinase 1: Modeling the Interaction between the Kinase and the Pleckstrin Homology Domains.

The 3-phosphoinositide-dependent protein kinase 1 (PDK1) K465E mutant kinase can still activate protein kinase B (PKB) at the membrane in a phosphatidylinositol-3,4,5-trisphosphate (PIP3, PtdIns(3,4,5)P3) independent manner. To understand this new PDK1 regulatory mechanism, docking and molecular dynamics calculations were performed for the first time to simulate the wild-type kinase domain-pleckstrin homology (PH) domain complex with PH-in and PH-out conformations. These simulations were then compared to the PH-in model of the KD-PH(mutant K465E) PDK1 complex. Additionally, three KD-PH complexes were simulated, including a substrate analogue bound to a hydrophobic pocket (denominated the PIF-pocket) substrate-docking site. We find that only the PH-out conformation, with the PH domain well-oriented to interact with the cellular membrane, is active for wild-type PDK1. In contrast, the active conformation of the PDK1 K465E mutant is PH-in, being ATP-stable at the active site while the PIF-pocket is more accessible to the peptide substrate. We corroborate that both the docking-site binding and the catalytic activity are in fact enhanced in knock-in mouse samples expressing the PDK1 K465E protein, enabling the phosphorylation of PKB in the absence of PIP3 binding.

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

Recent experimental results have challenged conventional views on the role metals play in the chemistry of protein kinases because it has been shown that (cAMP)-dependent protein kinase (PKA) is active in the presence of other divalent alkaline earth metal cations besides physiological Mg2+ ions. This has raised the important possibility that Ca2+ may also be a physiological cofactor of protein kinases. In this work, QM/MM calculations, at the DFT and MP2 levels for the QM part, on complete solvated models of PKAc-M2ATP-substrate ternary complexes, with PKAc as the catalytic subunit of PKA, M denoting Ca2+ or Sr2+ and substrate denoting SP20 or Kemptide, have been carried out for the overall phosphoryl transfer reaction. In accordance with the experimental data, our theoretical results show for the first time at the molecular level how the overall PKAc-catalyzed phosphorylation of SP20, via a dissociative mechanism, is plausible with Ca2+ and Sr2+. The viability of the catalytic reaction with Kemptide and Ca2+ is also verified here. The energy barrier of the rate-limiting phosphoryl-transfer step does not depend on different coordination environments of the alkaline earth metal cations whereas the proton-transfer step region is metal dependent making the global chemical process more exoergic on going from Mg2+ to Sr2+. This trend is in agreement with the less effective release of the phosphorylated product observed experimentally in the presence of Ca2+versus Mg2+, and would explain also the lower activity of PKAc with Ca2+, since phospho-substrate and ADP releases are rate limiting for catalytic turnover. For the same reason, we predict an even lower activity of PKAc with Sr2+. Moreover, the active sites of the in silico reactant and product complexes and the available X-ray crystallographic structures show good agreement.

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.

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.

Theoretical study of the mechanism of the hydride transfer between ferredoxin-NADP+ reductase and NADP+: the role of Tyr303.

During photosynthesis, ferredoxin-NADP(+) reductase (FNR) catalyzes the electron transfer from ferredoxin to NADP(+) via its FAD cofactor. The final hydride transfer event between FNR and the nucleotide is a reversible process. Two different transient charge-transfer complexes form prior to and upon hydride transfer, FNR(rd)-NADP(+) and FNR(ox)-NADPH, regardless of the hydride transfer direction. Experimental structures of the FNR(ox):NADP(+) interaction have suggested a series of conformational rearrangements that might contribute to attaining the catalytically competent complex, but to date, no direct experimental information about the structure of this complex is available. Recently, a molecular dynamics (MD) theoretical approach was used to provide a putative organization of the active site that might represent a structure close to the transient catalytically competent interaction of Anabaena FNR with its coenzyme, NADP(+). Using this structure, we performed fully microscopic simulations of the hydride transfer processes between Anabaena FNR(rd)/FNR(ox) and NADP(+)/H, accounting also for the solvation. A dual-level QM/MM hybrid approach was used to describe the potential energy surface of the whole system. MD calculations using the finite-temperature string method combined with the WHAM method provided the potential of mean force for the hydride transfer processes. The results confirmed that the structural model of the reactants evolves to a catalytically competent transition state through very similar free energy barriers for both the forward and reverse reactions, in good agreement with the experimental hydride transfer rate constants reported for this system. This theoretical approach additionally provides subtle structural details of the mechanism in wild-type FNR and provides an explanation why Tyr303 makes possible the photosynthetic reaction, a process that cannot occur when this Tyr is replaced by a Ser.

Theoretical analysis of the catalytic mechanism of Helicobacter pylori glutamate racemase.

One of the most challenging open key questions behind the stereoinversion of D-glutamate and L-glutamate catalyzed by glutamate racemases is how those enzymes manage to generate the thermodynamically unfavorable reverse protonation state of the catalytic residue cysteine required for the proton abstraction from the α-carbon of glutamate. In this paper, we have used molecular dynamics (MD) simulations with a molecular mechanics force field along with QM/MM calculations starting from the crystal structure and from different MD snapshots to study the enantiomeric conversion of D-glutamate to L-glutamate catalyzed by the Helicobacter pylori glutamate racemase. Our results show that structural fluctuations of the enzyme-substrate complex, represented by the different snapshots, lead to reaction paths with different features and fates. The whole reaction, when it occurs, involves four successive proton transfers in three or four different steps. In the first step, Asp7 assists the deprotonation of D-glutamate by participating in general base catalysis with neutral Cys70 thiol. An analogous mechanism was previously found by some of us for the case of Bacillus subtilis glutamate racemase. This fact explains why that aspartate belongs to the group of strictly conserved residues.

Influence of the enzyme phosphorylation state and the substrate on PKA enzyme dynamics.

cAMP-dependent protein kinase (PKA) is one of the simplest and best understood members of the protein kinase family. In a previous study, we have theoretically studied the complex between PKA and the heptapeptide substrate Kemptide by classical molecular dynamics. On the basis of the results obtained for Kemptide, the aim of the present work is to explore how the different conditions, such as phosphorylation state, substrate, and mutations of key residues affect the enzyme dynamics. We have built different models of the complex; particularly we have focused our attention on two crystallographic structures which main difference consists in their phosphorylation state. The first one has the residue Thr197 modified into a phospho-threonine (pThr197); the second one, in addition to the same Thr197, has also the residue Ser338 modified into a phospho-serine (pSer338). In addition, we have analyzed the effect of the choice of the substrate by building a model of the PKA-SP20 Michaelis complex. Finally, we have theoretically studied the effect of the mutation of the highly conserved residue Asp166 that, experimentally, leads to a decrease of the reaction rate. The results of this study give insight into the dynamical states of the enzyme and their relationship with different elements of the model, which correspond to different natural or human guided situations of the active biological system.

Temperature dependence of the kinetic isotope effects in thymidylate synthase. A theoretical study.

In recent years, the temperature dependence of primary kinetic isotope effects (KIE) has been used as indicator for the physical nature of enzyme-catalyzed H-transfer reactions. An interactive study where experimental data and calculations examine the same chemical transformation is a critical means to interpret more properly temperature dependence of KIEs. Here, the rate-limiting step of the thymidylate synthase-catalyzed reaction has been studied by means of hybrid quantum mechanics/molecular mechanics (QM/MM) simulations in the theoretical framework of the ensemble-averaged variational transition-state theory with multidimensional tunneling (EA-VTST/MT) combined with Grote-Hynes theory. The KIEs were calculated across the same temperature range examined experimentally, revealing a temperature independent behavior, in agreement with experimental findings. The calculations show that the H-transfer proceeds with ∼91% by tunneling in the case of protium and ∼80% when the transferred protium is replaced by tritium. Dynamic recrossing coefficients are almost invariant with temperature and in all cases far from unity, showing significant coupling between protein motions and the reaction coordinate. In particular, the relative movement of a conserved arginine (Arg166 in Escherichia coli ) promotes the departure of a conserved cysteine (Cys146 in E. coli ) from the dUMP by polarizing the thioether bond thus facilitating this bond breaking that takes place concomitantly with the hydride transfer. These promoting vibrations of the enzyme, which represent some of the dimensions of the real reaction coordinate, would limit the search through configurational space to efficiently find those decreasing both barrier height and width, thereby enhancing the probability of H-transfer by either tunneling (through barrier) or classical (over-the-barrier) mechanisms. In other words, the thermal fluctuations that are coupled to the reaction coordinate, together with transition-state geometries and tunneling, are the same in different bath temperatures (within the limited experimental range examined). All these terms contribute to the observed temperature independent KIEs in thymidylate synthase.

A QM/MM study of the phosphoryl transfer to the Kemptide substrate catalyzed by protein kinase A. The effect of the phosphorylation state of the protein on the mechanism.

We present here a theoretical study of the phosphoryl transfer catalytic mechanism of protein kinase A, which is the best known member of the large protein kinase family. We have built different theoretical models of the complete PKA-Mg(2)-ATP-substrate system to explore the two most accepted reaction pathways, using for the first time in a reaction mechanism theoretical study, the heptapeptide substrate Kemptide, which is relevant for its high efficiency and small size. The effect of the protein configuration, as modeled by two different X-ray structures with different phosphorylation states and degrees of flexibility, has been analyzed. The results indicate that the environmental conditions can influence the availability of the pathways and thus the choice of the mechanism to be followed. In addition, the roles of the two active site conserved residues, Asp166 and Lys168, have been analyzed for each reaction mechanism.

Mechanism of the hydride transfer between Anabaena Tyr303Ser FNR(rd)/FNR(ox) and NADP+/H. A combined pre-steady-state kinetic/ensemble-averaged transition-state theory with multidimensional tunneling study.

The flavoenzyme ferredoxin-NADP(+) reductase (FNR) catalyzes the production of NADPH during photosynthesis. The hydride-transfer reactions between the Anabaena mutant Tyr303Ser FNR(rd)/FNR(ox) and NADP(+)/H have been studied both experimentally and theoretically. Stopped-flow pre-steady-state kinetic measurements have shown that, in contrast to that observed for WT FNR, the physiological hydride transfer from Tyr303Ser FNR(rd) to NADP(+) does not occur. Conversely, the reverse reaction does take place with a rate constant just slightly slower than for WT FNR. This latter process shows temperature-dependent rates, but essentially temperature independent kinetic isotope effects, suggesting the reaction takes place following the vibration-driven tunneling model. In turn, ensemble-averaged variational transition-state theory with multidimensional tunneling calculations provide reaction rate constant values and kinetic isotope effects that agree with the experimental results, the experimental and the theoretical values for the reverse process being noticeably similar. The reaction mechanism behind these hydride transfers has been analyzed. The formation of a close contact ionic pair FADH(-):NADP(+) surrounded by the polar environment of the enzyme in the reactant complex of the mutant might be the cause of the huge difference between the direct and the reverse reaction.

How the substrate D-glutamate drives the catalytic action of Bacillus subtilis glutamate racemase.

Molecular Dynamics simulations with a Molecular Mechanics force field and a quite complete exploration of the QM/MM potential energy surfaces have been performed to study the D-glutamate --> L-glutamate reaction catalyzed by Bacillus subtilis glutamate racemase. The results show that the whole process involves four successive proton transfers that occur in three different steps. The Michaelis complex is already prepared to make the first proton transfer (from Cys74 to Asp10) possible. The second step involves two proton transfers (from the alpha-carbon to Cys74, and from Cys185 to the alpha-carbon), which occurs in a concerted way, although highly asynchronic. Finally, in the third step, the nascent deprotonated Cys185 is protonated by His187. The positively charged ammonium group of the substrate plays a very important key role in the reaction. It accompanies each proton transfer in a concerted and coupled way, but moving itself in the opposite direction from Asp10 to His187. Thus, the catalytic action of Bacillus subtilis glutamate racemase is driven by its own substrate of the reaction, D-glutamate.

Critical role of substrate conformational change in the proton transfer process catalyzed by 4-oxalocrotonate tautomerase.

4-Oxalocrotonate tautomerase enzyme (4-OT) catalyzes the isomerization of 2-oxo-4-hexenedioate to 2-oxo-3-hexenedioate. The chemical process involves two proton transfers, one from a carbon of the substrate to the nitrogen of Pro1 and another from this nitrogen atom to a different carbon of the substrate. In this paper the isomerization has been studied using the combined quantum mechanical and molecular mechanical method with a dual-level treatment of the quantum subsystem employing the MPW1BK density functional as the higher level. Exploration of the potential energy surface shows that the process is stepwise, with a stable intermediate state corresponding to the deprotonated substrate and a protonated proline. The rate constant of the overall process has been evaluated using ensemble-averaged variational transition state theory, including the quantized vibrational motion of a primary zone of active-site atoms and a transmission coefficient based on an ensemble of optimized reaction coordinates to account for recrossing trajectories and optimized multidimensional tunneling. The two proton-transfer steps have similar free energy barriers, but the transition state associated with the first proton transfer is found to be higher in energy. The calculations show that reaction progress is coupled to a conformational change of the substrate, so it is important that the simulation allows this flexibility. The coupled conformational change is promoted by changes in the electron distribution of the substrate that take place as the proton transfers occur.

Comparative study of the prereactive protein kinase A Michaelis complex with kemptide substrate.

In the present work we have modeled the Michaelis complex of the cyclic-Adenosine Monophosphate Dependent (cAMD) Protein Kinase A (PKA) with Mg(2)ATP and the heptapeptide substrate Kemptide by classical molecular dynamics. The chosen synthetic substrate is relevant for its high efficiency and small size, and it has not been used in previous theoretical studies. The structural analysis of the data generated along the 6 ns simulation indicates that the modeled substrate-enzyme complex mimics the substrate binding pattern known for PKA. The values of the average prereactive distances obtained from the simulation do not exclude any of the two limiting situations proposed as mechanisms in the literature for the phosphorylation reaction (dissociative and associative) because the system oscillates between configurations compatible with each of them. Furthermore, the results obtained for the average interaction distances between active site residues concord in suggesting the plausibility of an alternative third reaction mechanism.

New insights into the reaction mechanism catalyzed by the glutamate racemase enzyme: pH titration curves and classical molecular dynamics simulations.

The mechanism of the reactions catalyzed by the pyridoxal-phosphate-independent amino acid racemases and epimerases faces the difficult task of deprotonating a relatively low acidicity proton, the amino acid's alpha-hydrogen, with a relatively poor base, a cysteine. In this work, we propose a mechanism for one of these enzymes, glutamate racemase (MurI), about which many controversies exist, and the roles that its active site residues may play. The titration curves and the pK1/2 values of all of the ionizable residues for different structures leading from reactants to products have been analyzed. From these results a concerted mechanism has been proposed in which the Cys70 residue would deprotonate the alpha-hydrogen of the substrate while, at the same time, being deprotonated by the Asp7 residue. To study the consistency of this mechanism classical molecular dynamics (MD) simulations have been carried out along with pK1/2 calculations on the MD-generated structures.

A theoretical analysis of rate constants and kinetic isotope effects corresponding to different reactant valleys in lactate dehydrogenase.

In some enzymatic systems large conformational changes are coupled to the chemical step, in such a way that dispersion of rate constants can be observed in single-molecule experiments, each corresponding to the reaction from a different reactant valley. Under this perspective here we present a computational study of pyruvate to lactate transformation catalyzed by lactate dehydrogenase. The reaction consists of a hydride transfer and a proton transfer that seem to take place concertedly. The degree of asynchronicity and the energy barrier depend on the particular starting reactant valley. In order to estimate rate constants we used a free energy perturbation technique adapted to follow the intrinsic reaction coordinate for several possible reaction paths. Tunneling effects are also obtained with a slightly modified version of the ensemble-averaged variational transition state theory with multidimensional tunneling contributions. According to our results the closure of the active site by means of a flexible loop can lead to the formation of different reactant complexes displaying different features in the disposition of some key residues (such as Arg109), interactions with the substrate and number of water molecules in the active site. The chemical step of the reaction takes place with a different reaction rate from each of these complexes. Finally, primary kinetic isotope effects for replacement of the transferring hydrogen of the cofactor with a deuteride are in good agreement with experimental observations, thus validating our methodology.

Enzyme dynamics and tunneling enhanced by compression in the hydrogen abstraction catalyzed by soybean lipoxygenase-1.

A fully microscopical simulation of the rate-limiting hydrogen abstraction catalyzed by soybean lipoxygenase-1 (SLO-1) has been carried out. This enzyme exhibits the largest, and weakly temperature dependent, experimental H/D kinetic isotope effect (KIE) reported for a biological system. The theoretical model used here includes the complete enzyme with a solvation shell of water molecules, the Fe(III)-OH- cofactor, and the linoleic acid substrate. We have used a hybrid QM(PM3/d-SRP)/MM method to describe the potential energy surface of the whole system, and the ensemble-averaged variational transition-state theory with multidimensional tunneling (EA-VTST/MT) to calculate the rate constant and the primary KIE. The computational results show that the compression of the wild-type active site enzyme results in the huge contribution of tunneling (99%) to the rate of the hydrogen abstraction. Importantly, the active site becomes more flexible in the Ile553Ala mutant reactant complex simulation (for which a markedly temperature dependent KIE has been experimentally determined), thus justifying the proposed key role of the gating promoting mode in the reaction catalyzed by SLO-1. Finally, the results indicate that the calculated KIE for the wild-type enzyme has an important dependence on the barrier width.

On the ionization state of the substrate in the active site of glutamate racemase. A QM/MM study about the importance of being zwitterionic.

Computer simulations on a QM/MM potential energy surface have been carried out to gain insights into the catalytic mechanism of glutamate racemase (MurI). Understanding such a mechanism is a challenging task from the chemical point of view because it involves the deprotonation of a low acidic proton by a relatively weak base to give a carbanionic intermediate. First, we have examined the dependency of the kinetics and thermodynamics of the racemization process catalyzed by MurI on the ionization state of the substrate (glutamate) main chain. Second, we have employed an energy decomposition procedure to study the medium effect on the enzyme-substrate electrostatic and polarization interactions along the reaction. Importantly, the present theoretical results quantitatively support the mechanistic proposal by Rios et al. [J. Am. Chem. Soc. 2000, 122, 9373-9385] for the PLP-independent amino acid racemases.

Path Integral Simulations of Proton Transfer Reactions in Aqueous Solution Using Combined QM/MM Potentials.

A bisection sampling method was implemented in path integral simulations of chemical reactions in solution in the framework of the quantized classical path approach. In the present study, we employ a combined quantum mechanical and molecular mechanical (QM/MM) potential to describe the potential energy surface and the path integral method to incorporate nuclear quantum effects. We examine the convergence of the bisection method for two proton-transfer reactions in aqueous solution at room temperature. The first reaction involves the symmetrical proton transfer between an ammonium ion and an ammonia molecule. The second reaction is the ionization of nitroethane by an acetate ion. To account for nuclear quantum mechanical corrections, it is sufficient to quantize the transferring light atom in the ammonium ion-ammonia reaction, while it is necessary to also quantize the donor and acceptor atoms in the nitroethane-acetate ion reaction. Kinetic isotope effects have been computed for isotopic substitution of the transferring proton by a deuteron in the nitroethane-acetate reaction. In all computations, it is important to employ a sufficient number of polymer beads along with a large number of configurations to achieve convergence in these simulations.

Nonperfect synchronization of reaction center rehybridization in the transition state of the hydride transfer catalyzed by dihydrofolate reductase.

It has been suggested that the magnitudes of secondary kinetic isotope effects (2 degrees KIEs) of enzyme-catalyzed reactions are an indicator of the extent of reaction-center rehybridization at the transition state. A 2 degrees KIE value close to the corresponding secondary equilibrium isotope effects (2 degrees EIE) is conventionally interpreted as indicating a late transition state that resembles the final product. The reliability of using this criterion to infer the structure of the transition state is examined by carrying out a theoretical investigation of the hybridization states of the hydride donor and acceptor in the Escherichia coli dihydrofolate reductase (ecDHFR)-catalyzed reaction for which a 2 degrees KIE close to the 2 degrees EIE was reported. Our results show that the donor carbon at the hydride transfer transition state resembles the reactant state more than the product state, whereas the acceptor carbon is more productlike, which is a symptom of transition state imbalance. The conclusion that the isotopically substituted carbon is reactant-like disagrees with the conclusion that would have been derived from the criterion of 2 degrees KIEs and 2 degrees EIEs, but the breakdown of the correlation with the equilibrium isotope effect can be explained by considering the effect of tunneling.

A Molecular Dynamics Simulation of the Binding Modes of d-Glutamate and d-Glutamine to Glutamate Racemase.

Classical molecular dynamics simulations of the d-Gln/Aquifex pyrophilus MurI and d-Glu/Aquifex pyrophilus MurI complexes have been carried out. Since the active site of the enzyme contains many charged and polar residues, several binding modes are possible. Thus, three very different stable conformations of the substrate analogue d-Gln have been found, and at least three binding modes are possible for the substrate d-Glu. These qualitative results give an explanation for the apparent disagreement between the d-Gln bound MurI X-ray crystal structure and the expected position and orientation of the substrate d-Glu in order to make it possible the assumed Cα deprotonation (by Cys70)/reprotonation (by Cys178) racemization mechanism.