J. Javier Ruiz-Pernía

Unraveling the role of protein dynamics in dihydrofolate reductase catalysis

Significance The role of protein dynamics in enzyme catalysis remains a topic of considerable debate. Here, we use a combination of experimental and computational methods to identify the origins of the observed changes in reactivity on isotopic substitution of dihydrofolate reductase from Escherichia coli. Isotopic substitution causes differences in environmental coupling to the hydride transfer step and protein dynamics have therefore a small but measurable effect on the chemical reaction rate. Protein dynamics have controversially been proposed to be at the heart of enzyme catalysis, but identification and analysis of dynamical effects in enzyme-catalyzed reactions have proved very challenging. Here, we tackle this question by comparing an enzyme with its heavy (15N, 13C, 2H substituted) counterpart, providing a subtle probe of dynamics. The crucial hydride transfer step of the reaction (the chemical step) occurs more slowly in the heavy enzyme. A combination of experimental results, quantum mechanics/molecular mechanics simulations, and theoretical analyses identify the origins of the observed differences in reactivity. The generally slightly slower reaction in the heavy enzyme reflects differences in environmental coupling to the hydride transfer step. Importantly, the barrier and contribution of quantum tunneling are not affected, indicating no significant role for “promoting motions” in driving tunneling or modulating the barrier. The chemical step is slower in the heavy enzyme because protein motions coupled to the reaction coordinate are slower. The fact that the heavy enzyme is only slightly less active than its light counterpart shows that protein dynamics have a small, but measurable, effect on the chemical reaction rate.

Increased Dynamic Effects in a Catalytically Compromised Variant of Escherichia coli Dihydrofolate Reductase

Isotopic substitution (15N, 13C, 2H) of a catalytically compromised variant of Escherichia coli dihydrofolate reductase, EcDHFR-N23PP/S148A, has been used to investigate the effect of these mutations on catalysis. The reduction of the rate constant of the chemical step in the EcDHFR-N23PP/S148A catalyzed reaction is essentially a consequence of an increase of the quasi-classical free energy barrier and to a minor extent of an increased number of recrossing trajectories on the transition state dividing surface. Since the variant enzyme is less well set up to catalyze the reaction, a higher degree of active site reorganization is needed to reach the TS. Although millisecond active site motions are lost in the variant, there is greater flexibility on the femtosecond time scale. The “dynamic knockout” EcDHFR-N23PP/S148A is therefore a “dynamic knock-in” at the level of the chemical step, and the increased dynamic coupling to the chemical coordinate is in fact detrimental to catalysis. This finding is most likely applicable not just to hydrogen transfer in EcDHFR but also to other enzymatic systems.

Studying the role of protein dynamics in an SN2 enzyme reaction using free-energy surfaces and solvent coordinates

Conformational changes are known to be able to drive an enzyme through its catalytic cycle, allowing, for example, substrate binding or product release. However, the influence of protein motions on the chemical step is a controversial issue. One proposal is that the simple equilibrium fluctuations incorporated into transition-state theory are insufficient to account for the catalytic effect of enzymes and that protein motions should be treated dynamically. Here, we propose the use of free-energy surfaces, obtained as a function of both a chemical coordinate and an environmental coordinate, as an efficient way to elucidate the role of protein structure and motions during the reaction. We show that the structure of the protein provides an adequate environment for the progress of the reaction, although a certain degree of flexibility is needed to attain the full catalytic effect. However, these motions do not introduce significant dynamical corrections to the rate constant and can be described as equilibrium fluctuations.

Heavy enzymes — experimental and computational insights in enzyme dynamics

The role of protein motions in the chemical step of enzyme-catalyzed reactions is the subject of an open debate in the scientific literature. The systematic use of isotopically substituted enzymes has been revealed as a useful tool to quantify the role of these motions. According to the Born-Oppenheimer approximation, changing the mass of the protein does not change the forces acting on the system but alters the frequencies of the protein motions, which in turn can affect the rate constant. Experimental and theoretical studies carried out in this field are presented in this article and discussed in the framework of Transition State Theory.

Toward an Automatic Determination of Enzymatic Reaction Mechanisms and Their Activation Free Energies

We present a combination of the string method and a path collective variable for the exploration of the free energy surface associated to a chemical reaction in condensed environments. The on-the-fly string method is employed to find the minimum free energy paths on a multidimensional free energy surface defined in terms of interatomic distances, which is a convenient selection to study bond forming/breaking processes. Once the paths have been determined, a reaction coordinate is defined as a measure of the advance of the system along these paths. This reaction coordinate can be then used to trace the reaction Potential of Mean Force from which the activation free energy can be obtained. This combination of methodologies has been here applied to the study, by means of Quantum Mechanics/Molecular Mechanics simulations, of the reaction catalyzed by guanidinoacetate methyltransferase. This enzyme catalyzes the methylation of guanidinoacetate by S-adenosyl-l-methionine, a reaction that involves a methyl transfer and a proton transfer and for which different reaction mechanisms have been proposed.

Hydrolysis of phosphotriesters: a theoretical analysis of the enzymatic and solution mechanisms.

A theoretical study on the alkaline hydrolysis of paraoxon, one of the most popular organophosphorus pesticides, in aqueous solution and in the active site of Pseudomonas diminuta phosphotriesterase (PTE) is presented. Simulations by means of hybrid quantum mechanics/molecular mechanics (QM/MM) potentials show that the hydrolysis of paraoxon takes place through an A(N)D(N) or associative mechanism both in solution and in the active site of PTE. The results correctly reproduce the magnitude of the activation free energies and can be used to rationalize the observed kinetic isotope effects (KIEs) for the hydrolysis of paraoxon in both media. Enzymatic hydrolysis of O,O-diethyl p-chlorophenyl phosphate, a phosphotriester having a leaving group with higher pK(a) than paraoxon, was also simulated. Hydrolysis of this phosphotriester by PTE follows a A(N)+D(N) mechanism with a pentacoordinate intermediate. Moreover, the leaving group of this new substrate coordinates to one of the zinc ions of the bimetallic active site in order to stabilize the large negative charge developed on the oxygen atom of the leaving group when the P-O bond is broken in the products state. To accommodate this new ligand in the coordination shell, carbamylated Lys169 must be displaced from one zinc ion to the other, which in turn affects the acidity of Asp301, a residue originally bound to the second zinc ion. This ability to displace some of the ligands of the coordination shell of the zinc centers would explain the promiscuity of this enzyme, which is capable of catalyzing hydrolysis of different substrate by means of different mechanisms.

Chemical Ligation and Isotope Labeling to Locate Dynamic Effects during Catalysis by Dihydrofolate Reductase

Abstract Chemical ligation has been used to alter motions in specific regions of dihydrofolate reductase from E. coli and to investigate the effects of localized motional changes on enzyme catalysis. Two isotopic hybrids were prepared; one with the mobile N‐terminal segment containing heavy isotopes (2H, 13C, 15N) and the remainder of the protein with natural isotopic abundance, and the other one with only the C‐terminal segment isotopically labeled. Kinetic investigations indicated that isotopic substitution of the N‐terminal segment affected only a physical step of catalysis, whereas the enzyme chemistry was affected by protein motions from the C‐terminal segment. QM/MM studies support the idea that dynamic effects on catalysis mostly originate from the C‐terminal segment. The use of isotope hybrids provides insights into the microscopic mechanism of dynamic coupling, which is difficult to obtain with other studies, and helps define the dynamic networks of intramolecular interactions central to enzyme catalysis.

Analysis of the Decarboxylation Step in Mammalian Histidine Decarboxylase A COMPUTATIONAL STUDY

We report a hybrid quantum mechanics/molecular mechanics theoretical study on the reaction mechanism of mammalian histidine decarboxylase that allows us to obtain valuable insights on the structure of the cofactor-substrate adduct (external aldimine) in the active site of rat histidine decarboxylase. By means of molecular dynamics simulations, we traced the potential of mean force corresponding to the decarboxylation reaction of the adduct both in the active site of the enzyme and in aqueous solution. By comparing this process in both media, we have identified the key electrostatic interactions that explain the lowering of the free energy barrier in the enzyme. Our analysis also offers a validation of Dunathans hypothesis (Dunathan, H. (1966) Proc. Natl. Acad. Sci. U. S. A. 55, 712–716) regarding the role of stereoelectronic effects in the enzymatic decarboxylation process.

Protein isotope effects in dihydrofolate reductase from Geobacillus stearothermophilus show entropic-enthalpic compensatory effects on the rate constant.

Catalysis by dihydrofolate reductase from the moderately thermophilic bacterium Geobacillus stearothermophilus (BsDHFR) was investigated by isotope substitution of the enzyme. The enzyme kinetic isotope effect for hydride transfer was close to unity at physiological temperatures but increased with decreasing temperatures to a value of 1.65 at 5 °C. This behavior is opposite to that observed for DHFR from Escherichia coli (EcDHFR), where the enzyme kinetic isotope effect increased slightly with increasing temperature. These experimental results were reproduced in the framework of variational transition-state theory that includes a dynamical recrossing coefficient that varies with the mass of the protein. Our simulations indicate that BsDHFR has greater flexibility than EcDHFR on the ps-ns time scale, which affects the coupling of the environmental motions of the protein to the chemical coordinate and consequently to the recrossing trajectories on the reaction barrier. The intensity of the dynamic coupling in DHFRs is influenced by compensatory temperature-dependent factors, namely the enthalpic barrier needed to achieve an ideal transition-state configuration with minimal nonproductive trajectories and the protein disorder that disrupts the electrostatic preorganization required to stabilize the transition state. Together with our previous studies of other DHFRs, the results presented here provide a general explanation why protein dynamic effects vary between enzymes. Our theoretical treatment demonstrates that these effects can be satisfactorily reproduced by including a transmission coefficient in the rate constant calculation, whose dependence on temperature is affected by the protein flexibility.

A QM/MM Exploration of the Potential Energy Surface of Pyruvate to Lactate Transformation Catalyzed by LDH. Improving the Accuracy of Semiempirical Descriptions.

We present a QM/MM study of the potential energy surface of the pyruvate to lactate transformation catalyzed by L-lactate dehydrogenase. The transformation involves a hydride and a proton transfer which are followed by means of the corresponding antisymmetric combination of the distances from the hydrogen atom to the donor and the acceptor atoms. To discriminate among the possible reaction mechanisms we have considered different improvements of the AM1/MM description:  reoptimization of the van der Waals parameters and inclusion of corrections to the QM energy associated with both transfer coordinates. The QM subsystem has been also enlarged to include charge-transfer effects from the substrate to some specific residues. In our best treatment, the transformation is described as a concerted process through a single transition structure in which the hydride transfer is more advanced than the proton transfer. From the methodological point of view, the correction schemes tested here improve the quality of the semiempirical potential energy surface although they also present deficiencies attributed to consideration of the proton and hydride transfer corrections as separable ones.