Exploring chemical reactivity of complex systems with path-based coordinates: role of the distance metric.

Path-based reaction coordinates constitute a valuable tool for free-energy calculations in complex processes. When a reference path is defined by means of collective variables, a nonconstant distance metric that incorporates the nonorthonormality of these variables should be taken into account. In this work, we show that, accounting for the correct metric tensor, these kind of variables can provide iso-hypersurfaces that coincide with the iso-committor surfaces and that activation free energies equal the value that would be obtained if the committor function itself were used as reaction coordinate. The advantages of the incorporation of the variable metric tensor are illustrated with the an…

Catalytic Reaction Mechanism in Native and Mutant Catechol- O-methyltransferase from the Adaptive String Method and Mean Reaction Force Analysis.

Catechol- O-methyltransferase is an enzyme that catalyzes the methylation reaction of dopamine by S-adenosylmethionine, increasing the reaction rate by almost 16 orders of magnitude compared to the reaction in aqueous solution. Here, we combine the recently introduced adaptive string method and the mean reaction force method, in combination with the structural and electronic descriptors to characterize the reaction mechanism. The catalytic effect of the enzyme is addressed by the comparison of the reaction in the human wild-type enzyme, in the less effective Y68A mutant, and in aqueous solution. The influence of these different environments at different stages of the chemical process and th…

Elucidating the catalytic reaction mechanism of orotate phosphoribosyltransferase by means of X-ray crystallography and computational simulations

15 p.-8 fig.-2 tab.1 graph. abst.+ 8 fig. supl.-1 tab. supl.

Quantifying the limits of transition state theory in enzymatic catalysis

Significance Transition state theory (TST) is the most popular theory to calculate the rates of enzymatic reactions. However, in some cases TST could fail due to the violation of the nonrecrossing hypothesis at the transition state. In the present work we show that even for one of the most controversial enzymatic reactions—the hydride transfer catalyzed by dihydrofolate reductase—the error associated to TST represents only a minor correction to the reaction rate. Moreover, this error is actually larger for the reaction in solution than in the enzymatic active site. Based on this finding and on previous studies we propose an “enzymatic shielding” hypothesis which encompasses various aspects …

Publisher’s Note: “Transition state ensemble optimization for reactions of arbitrary complexity” [J. Chem. Phys. 143, 134111 (2015)]

Transition state ensemble optimization for reactions of arbitrary complexity.

In the present work, we use Variational Transition State Theory (VTST) to develop a practical method for transition state ensemble optimization by looking for an optimal hyperplanar dividing surface in a space of meaningful trial collective variables. These might be interatomic distances, angles, electrostatic potentials, etc. Restrained molecular dynamics simulations are used to obtain on-the-fly estimates of ensemble averages that guide the variations of the hyperplane maximizing the transmission coefficient. A central result of our work is an expression that quantitatively estimates the importance of the coordinates used for the localization of the transition state ensemble. Starting fro…

Dehydrochlorination of Hexachlorocyclohexanes Catalyzed by the LinA Dehydrohalogenase. A QM/MM Study

The elucidation of the catalytic role of LinA dehydrohalogenase in the degradation processes of hexachlorocyclohexane (HCH) isomers is extremely important to further studies on the bioremediation of HCH polluted areas. Herein, QM/MM free energy simulations are employed to provide the details of the dehydrochlorination reaction of two HCH isomers (γ and β). In particular, the role of the protonation state of one of the catalytic residues-His73-is explored. Based on our calculations, two distinct minimum free energy pathways (concerted and stepwise) were found for γ-HCH and β-HCH. The choice of the reaction channel for the dehydrochlorination reactions of γ- and β-HCH was shown to depend on t…

Linking Electrostatic Effects and Protein Motions in Enzymatic Catalysis. A Theoretical Analysis of Catechol O-Methyltransferase

The role of protein motions in enzymatic catalysis is the subject of a hot scientific debate. We here propose the use of an explicit solvent coordinate to analyze the impact of environmental motions during the reaction process. The example analyzed here is the reaction catalyzed by catechol O-methyltransferase, a methyl transfer reaction from S-adenosylmethionine (SAM) to the nucleophilic oxygen atom of catecholate. This reaction proceeds from a charged reactant to a neutral product, and then a large electrostatic coupling with the environment could be expected. By means of a two-dimensional free energy surface, we show that a large fraction of the environmental motions needed to attain the…

Reaction coordinates and transition states in enzymatic catalysis

Enzymatic reactions are complex chemical processes taking place in complex dynamic environments. Theoretical characterization of these reactions requires the determination of the reaction coordinate and the transition state ensemble. This is not an easy task because many degrees of freedom may be involved in principle. We present recent efforts to find good enzymatic reaction coordinates and the implications of these findings in the interpretation of enzymatic efficiency. In particular, we analyze different strategies based on the use of minimum free energy paths and direct localization of the dividing surface on multidimensional free energy surfaces. Another strategy is based on the genera…

Modeling caspase-1 inhibition: Implications for catalytic mechanism and drug design.

Abstract The metabolic product of caspase-1, IL-1β, is an important mediator in inflammation and pyroptosis cell death process. Alzheimer's disease, septic shock and rheumatoid arthritis are IL-1β mediated diseases, making the caspase-1 an interesting target of pharmacological value. Many inhibitors have been developed until now, most of them are peptidomimetic with improved potency. In the present study, all-atom molecular dynamics simulations and the MM/GBSA method were employed to reproduce and interpret the results obtained by in vitro experiments for a series of inhibitors. The analysis shows that the tautomeric state of the catalytic His237 impact significantly the performance of the …

Thermal Isomerization Mechanism in Dronpa and Its Mutants.

The photoswitching speed of the reversibly switchable fluorescent proteins (RSFPs) from the family of green fluorescent proteins (GFPs) changes upon mutation which is of direct importance for various high-resolution techniques. Dronpa is one of the most used RSFPs. Its point mutants rsFastLime (Dronpa V157G) and rsKame (Dronpa V157L) exhibit a striking difference in their photoswitching speed. Here the QM/MM on-the-fly string method is used in order to explore the details of the thermal isomerization mechanism. The four principal ways in which isomerization may occur have been scrutinized for each of the three proteins. It has been shown that thermal isomerization occurs via a one-bond-flip…

A Collective Coordinate to Obtain Free Energy Profiles for Complex Reactions in Condensed Phases.

Exploration of chemical reactions in complex explicit environments has become an affordable task with the use of hybrid quantum mechanics/molecular mechanics potentials which allow calculating free energy profiles of chemical reactions under the influence of the surroundings. Tracing these free energy profiles requires the selection of a reaction coordinate, which can be cumbersome for those processes involving more than a single chemical event in a concerted step. We here propose a collective coordinate to be used in the calculation of free energy profiles for complex reactions in condensed phases. This coordinate is based in the definition of the advance along a path introduced by Brandua…

Studying the phosphoryl transfer mechanism of the E. coli phosphofructokinase-2: from X-ray structure to quantum mechanics/molecular mechanics simulations

Phosphofructokinases (Pfks) catalyze the ATP-dependent phosphorylation of fructose-6-phosphate (F6P) and they are regulated in a wide variety of organisms. Although numerous aspects of the kinetics and regulation have been characterized for Pfks, the knowledge about the mechanism of the phosphoryl transfer reaction and the transition state lags behind. In this work, we describe the X-ray crystal structure of the homodimeric Pfk-2 from E. coli, which contains products in one site and reactants in the other, as well as an additional ATP molecule in the inhibitory allosteric site adjacent to the reactants. This complex was previously predicted when studying the kinetic mechanism of ATP inhibit…

How a Second Mg2+ Ion Affects the Phosphoryl-Transfer Mechanism in a Protein Kinase: A Computational Study

Mg2+ ions are essential for the proper functioning of protein kinases, and their roles in kinase activity have been studied for years. However, recent investigations have shed light on how these me...

A first peek into sub-picosecond dynamics of spin energy levels in magnetic biomolecules

We estimate the time- and temperature-evolution of spin energy levels in a metallopeptide by combining molecular dynamics with crystal field analysis. Fluctuations of tens of cm-1 for spin energy levels at fs times gradually average out at longer times. We confirm that local vibrations are key in spin dynamics.

Studying the phosphoryl transfer mechanism of the

Phosphofructokinases catalyze the ATP-dependent phosphorylation of fructose-6-phosphate and they are highly regulated.

Unraveling the Reaction Mechanism of Enzymatic C5-Cytosine Methylation of DNA. A Combined Molecular Dynamics and QM/MM Study of Wild Type and Gln119 Variant

M.HhaI is a DNA methyltransferase from Haemophilus hemolyticus that catalyzes the transfer of a methyl group from S-adenosyl-l-methionine (SAM) to the C5 position of a cytosine. This enzyme is a paradigmatic model for C5 DNA methyltransferases due to its major homology to mammalian enzymes and to the availability of high-resolution structures of the DNA–enzyme complex. In spite of the number of experimental and theoretical analyses carried out for this system, many mechanistic details remain unraveled. We have used full atomistic classical molecular dynamics simulations to explore the protein–SAM–DNA ternary complex, where the target cytosine base is flipped out into the active site for bot…

Dynamics and reactivity in Thermus aquaticus N6-adenine methyltransferase.

M.TaqI is a DNA methyltransferase from Thermus aquaticus that catalyzes the transfer of a methyl group from S-adenosyl-l-methionine to the N6 position of an adenine, a process described only in prokaryotes. We have used full atomistic classical molecular dynamics simulations to explore the protein–SAM–DNA ternary complex where the target adenine is flipped out into the active site. Key protein–DNA interactions established by the target adenine in the active site are described in detail. The relaxed structure was used for a combined quantum mechanics/molecular mechanics exploration of the reaction mechanism using the string method. According to our free energy calculations the reaction takes…

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. Th…

Adaptive Finite Temperature String Method in Collective Variables.

Here we present a modified version of the on-the-fly string method for the localization of the minimum free energy path in a space of arbitrary collective variables. In the proposed approach the shape of the biasing potential is controlled by only two force constants, defining the width of the potential along the string and orthogonal to it. The force constants and the distribution of the string nodes are optimized during the simulation, improving the convergence. The optimized parameters can be used for umbrella sampling with a path CV along the converged string as the reaction coordinate. We test the new method with three fundamentally different processes: chloride attack to chloromethane…

How a Second Mg2+ Ion Affects the Phosphoryl Transfer Mechanism in a Protein Kinase: A Computational Study

<div>In this contribution, the phosphoryl transfer reaction in CDK2 has been studied in detail considering the presence of an additional Mg2+ ion in the active site. For this purpose, QM/MM (quantum mechanics/molecular mechanics) free energy calculations with the adaptive string method were performed, which showed that indeed the system containing two Mg2+ ions exhibits a lower activation free energy, corroborating the experimental observations.</div>