Mireia Garcia-Viloca

Hydride transfer catalyzed by xylose isomerase: Mechanism and quantum effects

We have applied molecular dynamics umbrella‐sampling simulation and ensemble‐averaged variational transition state theory with multidimensional tunneling (EA‐VTST/MT) to calculate the reaction rate of xylose‐to‐ xylulose isomerization catalyzed by xylose isomerase in the presence of two Mg2+ ions. The calculations include determination of the free energy of activation profile and ensemble averaging in the transmission coefficient. The potential energy function is approximated by a combined QM/MM/SVB method involving PM3 for the quantum mechanical (QM) subsystem, CHARMM22 and TIP3P for the molecular mechanical (MM) environment, and a simple valence bond (SVB) local function of two bond distances for the hydride transfer reaction. The simulation confirms the essential features of a mechanism postulated on the basis of kinetics and X‐ray data by Whitlow et al. (Whitlow, M.; Howard, A. J.; Finzel, B. C.; Poulos, T. L.; Winborne, E.; Gilliland, G. L. Proteins 1991, 9, 153) and Ringe, Petsko, and coworkers (Labie, A.; Allen, K.‐N.; Petsko, G. A.; Ringe, D. Biochemistry 1994, 33, 5469). This mechanism involves a rate‐determining 1,2‐hydride shift with prior and post proton transfers. Inclusion of quantum mechanical vibrational energy is important for computing the free energy of activation, and quantum mechanical tunneling effects are essential for computing kinetic isotope effects (KIEs). It is found that 85% of the reaction proceeds by tunneling and 15% by overbarrier events. The computed KIE for the ratio of hydride to deuteride transfer is in good agreement with the experimental results. The molecular dynamics simulations reveal that proton and hydride transfer reactions are assisted by breathing motions of the mobile Mg2+ ion in the active site, providing evidence for concerted motion of Mg2+ during the hydride transfer step.

Inclusion of Quantum Mechanical Vibrational Energy in Reactive Potentials of Mean Force

Classical molecular dynamics and Monte Carlo simulations typically exclude quantum effects on the vibrations of reactants and transition states, and this may lead to significant errors in the computed potential of mean force. To correct this deficiency, a simple approximate procedure is proposed for the inclusion of quantum-mechanical vibrational energy in the computation of reactive potentials of mean force in condensed phases. The method is illustrated by a hydrogen atom transfer and a proton transfer reaction in water, in particular, the 1,5-sigmatropic shift in malonaldehyde and the intermolecular proton shift between ammonium ion and ammonia in an encounter complex. In both cases, quantum-mechanical vibrational energy makes significant contributions by reducing the free energy of activation by 2 to 3 kcal/mol. This finding has important implications in developing empirical potential functions for the study of enzyme reactions, and it is essential to quantize vibrational energy in the computed potenti...

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.

Importance of substrate and cofactor polarization in the active site of dihydrofolate reductase

By using a combined quantum-mechanical and molecular-mechanical potential in molecular dynamics simulations, we have investigated the effects of the enzyme electric field of dihydrofolate reductase on the electronic polarization of its 5-protonated dihydrofolate substrate at various stages of the catalyzed hydride transfer reaction. Energy decomposition of the total electrostatic interaction energy between the ligands and the enzyme shows that the polarization effect is 4% of the total electrostatic interaction energy, and, significantly, it accounts for 9kcal/mol of transition state stabilization relative to the reactant state. Therefore it is essential to take account of substrate polarization for quantitative interpretation of enzymatic function and for calculation of binding free energies of inhibitors to a protein. Atomic polarizations are calculated as the differences in the average atomic charges on the atoms in gas phase and in molecular simulations of the enzyme; this analysis shows that the glutamate tail and the pterin ring are the highly polarized regions of the substrate. Electron density difference plots of the reactant and product complexes at instantaneous configurations in the enzyme active center confirm the inferences made on the basis of partial atomic charges.

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.

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.

Mechanism of the Hydride Transfer between Anabaena Tyr303Ser FNRrd/FNRox 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.

Sensitivity of molecular dynamics simulations to the choice of the X-ray structure used to model an enzymatic reaction

A subject of great practical importance that has not received much attention is the question of the sensitivity of molecular dynamics simulations to the initial X‐ray structure used to set up the calculation. We have found two cases in which seemingly similar structures lead to quite different results, and in this article we present a detailed analysis of these cases. The first case is acyl‐CoA dehydrogenase, and the chief difference of the two structures is attributed to a slight shift in a backbone carbonyl that causes a key residue (the proton‐abstracting base) to be in a bad conformation for reaction. The second case is xylose isomerase, and the chief difference of the two structures appears to be the ligand sphere of a Mg2+ metal cofactor that plays an active role in catalysis.

Solvent effects, reaction coordinates, and reorganization energies on nucleophilic substitution reactions in aqueous solution

Publisher Summary This chapter reviews a combined quantum mechanical and molecular mechanical (QM/MM) method that couples molecular orbital and valence bond theories, called the MOVB method, to determine the free energy reaction profiles, or potentials of mean force (PMF), for chemical reactions in solution. The combined QM-MOVB/MM method is applied to three types of nucleophilic substitution reactions and the question of the dependence of the computed PMF on the choice of the reaction coordinate and its implications on the computed solvent reorganization energy have been addressed. A mixed MOVB model has been described for describing the potential energy surface of reactive systems, and results from applications to S N 2 reactions in aqueous solution have been presented. The MOVB model is based on a block-localized wave function method to define diabatic electronic state functions. A configuration interaction Hamiltonian has been constructed using these diabatic VB states as basis functions. The computed geometrical and energetic results for these systems are in accord with the experimental and theoretical studies. These studies show that the MOVB model can be adequately used as a mapping potential to derive solvent reaction coordinates for condensed phase processes.

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.