Seeking the Source of Catalytic Efficiency of Lindane Dehydrochlorinase, LinA.

Herein we present the results of an in-depth simulation study of LinA and its two variants. In our analysis, we combined the exploration of protein conformational dynamics with and without bound substrates (hexachlorocyclohexane (HCH) isomers) performed using molecular dynamics simulation followed by the extraction of the most frequently visited conformations and their characteristics with a detailed description of the interactions taking place in the active site between the respective HCH molecule and the first shell residues by using symmetry-adapted perturbation theory (SAPT) calculations. A detailed investigation of the conformational space of LinA substates has been accompanied by desc…

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…

Theoretical insights in enzyme catalysis

In this tutorial review we show how the methods and techniques of computational chemistry have been applied to the understanding of the physical basis of the rate enhancement of chemical reactions by enzymes. This is to answer the question: Why is the activation free energy in enzyme catalysed reactions smaller than the activation free energy observed in solution? Two important points of view are presented: Transition State (TS) theories and Michaelis Complex (MC) theories. After reviewing some of the most popular computational methods employed, we analyse two particular enzymatic reactions: the conversion of chorismate to prephenate catalysed by Bacillus subtilis chorismate mutase, and a m…

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…

Minimization of dynamic effects in the evolution of dihydrofolate reductase

Protein isotope labeling is a powerful technique to probe functionally important motions in enzyme catalysis and can be applied to investigate the conformational dynamics of proteins.

Molecular surface area and hydrophobic effect.

Calculation of binding energy using BLYP/MM for the HIV-1 integrase complexed with the S-1360 and two analogues.

Abstract Integrase (IN) is one of the three human immunodeficiency virus type 1 (HIV-1) enzymes essential for effective viral replication. S-1360 is a potent and selective inhibitor of HIV-1 IN. In this work, we have carried out molecular dynamics (MD) simulations using a hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) approach, to determine the protein–ligand interaction energy for S-1360 and two analogues. Analysis of the MD trajectories reveals that the strongest protein–inhibitor interactions, observed in the three studied complexes, are established with Lys-159 residue and Mg 2+ cation. Calculations of binding energy using BLYP/MM level of theory reveal that there is a direct rela…

Reversibility and Diffusion in Mandelythiamin Decarboxylation. Searching Dynamical Effects in Decarboxylation Reactions

Decarboxylation of mandelylthiamin in aqueous solution is analyzed by means of quantum mechanics/molecular mechanics simulations including solvent effects. The free energy profile for the decarboxylation reaction was traced, assuming equilibrium solvation, while reaction trajectories allowed us to incorporate nonequilibrium effects due to the solvent degrees of freedom as well as to evaluate the rate of the diffusion process in competition with the backward reaction. Our calculations that reproduce the experimental rate constant show that decarboxylation takes place with a non-negligible free energy barrier for the backward reaction and that diffusion of carbon dioxide is very fast compared…

Peptide Bond Formation Mechanism Catalyzed by Ribosome

In this paper we present a study of the peptide bond formation reaction catalyzed by ribosome. Different mechanistic proposals have been explored by means of Free Energy Perturbation methods within hybrid QM/MM potentials, where the chemical system has been described by the M06-2X functional and the environment by means of the AMBER force field. According to our results, the most favorable mechanism in the ribosome would proceed through an eight-membered ring transition state, involving a proton shuttle mechanism through the hydroxyl group of the sugar and a water molecule. This transition state is similar to that described for the reaction in solution (J. Am. Chem. Soc. 2013, 135, 8708–871…

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 (QM/MM) 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 subs…

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 …

Water-Assisted Alkaline Hydrolysis of Monobactams: A Theoretical Study

A theoretical study of the water-assisted alkaline hydrolysis of 2-azetidinone, 3-formylamino-2-azetidinone and 3-formylamino-2-azetidine-1-sulfonate ion is carried out at the B3LYP/6-31+G* level. The effect of bulk solvent is taken into account using the PCM solvation model while specific solvent effects are represented by the inclusion of an ancillary water molecule along the reaction profile. The calculated free energy barriers in solution are in reasonable agreement with experimental values. The observed substituent effects due to the presence of the 3-formylamino and the SO(3) groups attached to the beta-lactam ring are crucial factors determining the hydrolysis of monobactam antibioti…

The Catalytic Mechanism of Carboxylesterases: A Computational Study

The catalytic mechanism of carboxylesterases (CEs, EC 3.1.1.1) is explored by computational means. CEs hydrolyze ester, amide, and carbamate bonds found in xenobiotics and endobiotics. They can also perform transesterification, a reaction important, for instance, in cholesterol homeostasis. The catalytic mechanisms with three different substrates (ester, thioester, and amide) have been established at the M06-2X/6-311++G**//B3LYP/6-31G* level of theory. It was found that the reactions proceed through a mechanism involving four steps instead of two as is generally proposed: (i) nucleophilic attack of serine to the substrate, forming the first tetrahedral intermediate, (ii) formation of the ac…

Theoretical Modeling on the Reaction Mechanism of p-Nitrophenylmethylphosphate Alkaline Hydrolysis and its Kinetic Isotope Effects

We have studied the alkaline hydrolysis of p-nitrophenylmethylphosphate (p-NPmP) in aqueous solution by means of polarizable continuum models and by hybrid quantum-mechanical/molecular-mechanical (QM/MM) methods. The theoretical predictions of kinetic isotope effects (KIEs) are in very good agreement with the experimental data, confirming a concerted asynchronous molecular mechanism. In addition, comparison of high level DFT theory with semiempirical AM1/d Hamiltonian has allowed checking the reliability of the later to be used in modeling very large molecular models containing phosphorus atoms.

The mechanism of formamide hydrolysis in water from ab initio calculations and simulations.

The neutral hydrolysis of formamide in water is a suitable reference to quantify the efficiency of proteolytic enzymes. However, experimental data for this reaction has only very recently been obtained and the kinetic constant determined experimentally is significantly higher than that predicted by previous theoretical estimations. In this work, we have investigated in detail the possible mechanisms of this reaction. Several solvent models have been considered that represent a considerable improvement on those used in previous studies. Density functional and ab initio calculations have been carried out on a system which explicitly includes the first solvation shell of the formamide molecule…

Molecular dynamics simulation in aqueous solution of N -methylazetidinone as a model of β -lactam antibiotics

In this article, we analyze the results of a molecular dynamics simulation in aqueous solution of the N-methylazetidinone molecule, often used to model β-lactam antibiotics. The radial distribution functions (RDFs) corresponding to the most interesting atoms, in terms of reactivity, are presented. We focus our study on the effect of a polar environment on the molecule. The solvent structure around the system is compared to the structure of β-lactam-water complexes, as obtained in a previous study of reaction mechanisms for the neutral and alkaline hydrolyses of N-methylazetidinone. Two types of complexes have been considered which are related to different hydrolysis mechanisms having simila…

Promiscuity in alkaline phosphatase superfamily. Unraveling evolution through molecular simulations.

We here present a theoretical study of the alkaline hydrolysis of a phosphodiester (methyl p-nitrophenyl phosphate or MpNPP) in the active site of Escherichia coli alkaline phosphatase (AP), a monoesterase that also presents promiscuous activity as a diesterase. The analysis of our simulations, carried out by means of molecular dynamics (MD) simulations with hybrid quantum mechanics/molecular mechanics (QM/MM) potentials, shows that the reaction takes place through a D(N)A(N) or dissociative mechanism, the same mechanism employed by AP in the hydrolysis of monoesters. The promiscuous activity observed in this superfamily can be then explained on the basis of a conserved reaction mechanism. …

Hypoxia-Selective Dissociation Mechanism of a Nitroimidazole Nucleoside in a DNA Environment

Photodynamic therapy is a promising approach to treat a variety of superficial tumors and other diseases. One of its major limitations arises from its dependence on molecular oxygen, which decreases the efficiency of the therapy in hypoxia conditions commonly developed by solid tumors. The present contribution reveals the molecular mechanism of a modified thymine bearing a nitroimidazole substituent, a photosensitizer able to produce highly harmful interstrand cross-links in the DNA double strand after irradiation selectively in absence of oxygen. The mechanism is resolved at a fully atomistic and electronic level relying on quantum mechanics (CASPT2, coupled-cluster, DFT, and TD-DFT method…

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.

Computational design of biological catalysts

The purpose of this tutorial review is to illustrate the way to design new and powerful catalysts. The first possibility to get a biological catalyst for a given chemical process is to use existing enzymes that catalyze related reactions. The second possibility is the use of immune systems that recognize stable molecules resembling the transition structure of the target reaction. We finally show how computational techniques are able to provide an enormous quantity of information, providing clues to guide the development of new biological catalysts

The catalytic mechanism of glyceraldehyde 3-phosphate dehydrogenase from Trypanosoma cruzi elucidated via the QM/MM approach

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) has been identified as a key enzyme involved in glycolysis processes for energy production in the Trypanosoma cruzi parasite. This enzyme catalyses the oxidative phosphorylation of glyceraldehyde 3-phosphate (G3P) in the presence of inorganic phosphate (Pi) and nicotinamide adenosine dinucleotide (NAD+). The catalytic mechanism used by GAPDH has been intensively investigated. However, the individual roles of Pi and the C3 phosphate of G3P (Ps) sites, as well as some residues such as His194 in the catalytic mechanism, remain unclear. In this study, we have employed Molecular Dynamics (MD) simulations within hybrid quantum mechanical/molecular …

Preorganization and reorganization as related factors in enzyme catalysis: the chorismate mutase case.

In this paper a deeper insight into the chorismate-to prephenate-rearrangement, catalyzed by Bacillus subtilis chorismate mutase, is provided by means of a combination of statistical quantum mechanics/molecular mechanics simulation methods and hybrid potential energy surface exploration techniques. The main aim of this work is to present an estimation of the preorganization and reorganization terms of the enzyme catalytic rate enhancement. To analyze the first of these, we have studied different conformational equilibria of chorismate in aqueous solution and in the enzyme active site. Our conclusion is that chorismate mutase preferentially binds the reactive conformer of the substrate--that…

Molecular surface calculations on organic compounds

Abstract The molecular surface area was calculated for 82 hydrocarbons, esters, ethers, alcohols and ketones including linear, branched and cyclic ones. The relationship between the aqueous solubility and the molecular surface area was determined for each family of compounds and for all the families together. The results show that solubility is mainly determined by molecular surface area and that the influence of functional group is the same for all the oxygen containing compounds studied.

Multiscale Simulations of SARS-CoV-2 3CL Protease Inhibition with Aldehyde Derivatives. Role of Protein and Inhibitor Conformational Dynamics in the Reaction Mechanism

<p>We here investigate the mechanism of SARS-CoV-2 3CL protease inhibition by one of the most promising families of inhibitors, those containing an aldehyde group as warhead. These compounds are covalent inhibitors that inactivate the protease forming a stable hemithioacetal complex. Inhibitor 11a is a potent inhibitor that has been already tested in vitro and in animals. Using a combination of classical and QM/MM simulations we determined the binding mode of the inhibitor into the active site and the preferred rotameric state of the catalytic histidine. In the noncovalent complex the aldehyde group is accommodated into the oxyanion hole formed by the NH main chain groups of residues …

Regioselectivity of the OH Radical Addition to Uracil in Nucleic Acids. A Theoretical Approach Based on QM/MM Simulations.

Oxidation of nucleic acids is ubiquitous in living beings under metabolic impairments and/or exposed to external agents such as radiation, pollutants, or drugs, playing a central role in the development of many diseases mediated by DNA/RNA degeneration. Great efforts have been devoted to unveil the molecular mechanisms behind the OH radical additions to the double bonds of nucleobases; however, the specific role of the biological environment remains relatively unexplored. The present contribution tackles the study of the OH radical addition to uracil from the gas phase to a full RNA macromolecule by means of quantum-chemistry methods combined with molecular dynamics simulations. It is shown…

Do dynamic effects play a significant role in enzymatic catalysis? A theoretical analysis of formate dehydrogenase.

A theoretical study of the protein dynamic effects on the hydride transfer between the formate anion and nicotinamide adenine dinucleotide (NAD + ), catalyzed by formate dehydrogenase (FDH), is presented in this paper. The analysis of free downhill molecular dynamic trajectories, performed in the enzyme and compared with the reaction in aqueous solution, has allowed the study of the dynamic coupling between the reacting fragments and the protein or the solvent water molecules, as well as an estimation of the dynamic effect contribution to the catalytic effect from calculation of the transmission coefficient in the enzyme and in solution. The obtained transmission coefficients for the enzyme…

Unraveling the SARS-CoV-2 Main Protease Mechanism Using Multiscale DFT/MM Methods

<p>We present a detailed theoretical analysis of the reaction mechanism of proteolysis catalyzed by the main protease of SARS-CoV-2. Using multiscale simulation methods, we have characterized the interactions stablished by a peptidic substrate in the active site and then we have explored the free energy landscape associated to the acylation and de-acylation steps of the proteolysis reaction, characterizing the transition states of the process. Our mechanistic proposals can explain most of the experimental observations made on the highly similar ortholog protease of SARS-CoV. We point out to some key interactions that may facilitate the acylation process and thus can be crucial in the …

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

Substrate promiscuity in DNA methyltransferase M.PvuII. A mechanistic insight

M.PvuII is a DNA methyltransferase from the bacterium Proteus vulgaris that catalyzes methylation of cytosine at the N4 position. This enzyme also displays promiscuous activity catalyzing methylation of adenine at the N6 position. In this work we use QM/MM methods to investigate the reaction mechanism of this promiscuous activity. We found that N6 methylation in M.PvuII takes place by means of a stepwise mechanism in which deprotonation of the exocyclic amino group is followed by the methyl transfer. Deprotonation involves two residues of the active site, Ser53 and Asp96, while methylation takes place directly from the AdoMet cofactor to the target nitrogen atom. The same reaction mechanism…

Theoretical study of the temperature dependence of dynamic effects in thymidylate synthase.

A theoretical study of the temperature dependence of dynamic effects in the rate limiting step of the reaction catalyzed by thymidylate synthase is presented in this paper. From hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) optimizations of transition state structures within a fully flexible molecular model, free downhill molecular dynamics trajectories have been performed at four different temperatures. The analysis of the reactive and non-reactive trajectories in the enzyme environment has allowed us to study the geometric and electronic coupling between the substrate, the cofactor and the protein. The results show how the contribution of dynamic effects to the rate enhancement mea…

The solvent-excluding surface as a descriptor of ionic channels: Gramicidin-A

Abstract We have set out a methodology for checking the performance of the methods used to compute the Solvent-Excluding Surface. The method consists of computing the area of the Solvent-Excluding Surface of a model of channel, with a fixed pore size, for several values of the solvent radius. The graphical representation of the value of the area versus the radius of the solvent sphere shows a sharp change just at the radius of the pore. With this model we may analyze the ability of each method to describe small changes of the surface. We made the study with GEPOL93, older versions of GEPOL and MSDOT. The study is applied to a natural protein channel, as is Gramicidin-A, showing that this ty…

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…

Enzymatic effects on reactant and transition states. The case of chalcone isomerase.

Chalcone isomerase catalyzes the transformation of chalcone to naringerin as a part of flavonoid biosynthetic pathways. The global reaction takes place through a conformational change of the substrate followed by chemical reaction, being thus an excellent example to analyze current theories about enzyme catalysis. We here present a detailed theoretical study of the enzymatic action on the conformational pre-equilibria and on the chemical steps for two different substrates of this enzyme. Free-energy profiles are obtained in terms of potentials of mean force using hybrid quantum mechanics/molecular mechanics potentials. The role of the enzyme becomes clear when compared to the counterpart eq…

A molecular dynamics study on the role of the protonation state in the biosynthesis of R-PAC by AHAS

Abstract The effect of the protonation state of the hydroxyl-ethylthiamin diphosphate intermediate, HEThDP, on the enzyme-substrate interactions and their consequences on the biosynthesis of R-phenylacetylcarbinol, R-PAC, by the acetohydroxy acid synthase, AHAS, is addressed by molecular dynamics simulations. It is found that the form of HEThDP, which favors the formation of R-PAC, is that having the 4-aminopyrimidine ring with the N1′ atom protonated and the N4′ atom as aminopyrimidinium ion. Under this form both active sites of AHAS have the ability to perform the catalysis, unlike that observed for the other possible protonation states of N1′ and N4′ atoms.

Binding and Reactivity of a Nitrile Oral Inhibitor of SARS-CoV-2 Main Protease Revealed by Computational Simulations

We present a detailed analysis of the binding mode and reactivity of the novel oral inhibitor PF-07321332 developed against SARS-CoV-2 3CL protease. Classical and QM/MM Molecular Dynamics simulations are used to quantify the contributions to the binding free energy and the reaction mechanism for covalent inhibition. The small size of the nitrile warhead conferes additional advantadges to this inhibitor.

A QM/MM study of the reaction mechanism for the 3′-processing step catalyzed by HIV-1 integrase

Integrase (IN) is one of the three human immunodeficiency virus type 1 enzymes (HIV-1) essential for effective viral replication. This viral enzyme is involved in the integration of HIV DNA into host chromosomal DNA. In this work we have carried out molecular dynamics simulations using a hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) approach to study the reaction mechanism for the 3′-processing step of DNA integration using a model substrate. The results obtained by QM(AM1)/MM and QM(PM3)/MM simulations have been improved by single-point corrections using an ab initio method to describe the quantum subsystem. The results obtained within this computational model can be used to obtain …

Singlet Oxygen Attack on Guanine: Reactivity and Structural Signature within the B-DNA Helix

International audience; Oxidatively generated DNA lesions are numerous and versatile, and have been the subject of intensive research since the discovery of 8-oxoguanine in 1984. Even for this prototypical lesion, the precise mechanism of formation remains elusive due to the inherent difficulties in characterizing high-energy intermediates. We have probed the stability of the guanine endoperoxide in B-DNA as a key intermediate and determined a unique activation free energy of around 6 kcal mol−1 for the formation of the first C−O covalent bond upon the attack of singlet molecular oxygen (1O2) on the central guanine of a solvated 13 base-pair poly(dG-dC), described by means of quantum mechan…

DFT Study on the Interaction of Tris(benzene-1,2-dithiolato)molybdenum Complex with Water. A Hydrolysis Mechanism Involving a Feasible Seven-Coordinate Aquomolybdenum Intermediate

In the present work, the reactivity of the tris(benzene-1,2-dithiolato)molybdenum complex ([Mo(bdt)3]) toward water is studied by means of the density functional theory (DFT). DFT calculations were performed using the M06, B3P86, and B3PW91 hybrid functionals for comparison purposes. The M06 method was employed to elucidate the reaction pathway, relative stability of the intermediate products, nature of the Mo–S bond cleavage, and electronic structure of the involved molybdenum species. This functional was also used to study the transference of electrons from the molybdenum center toward the ligands. The reaction pathway confirms that [Mo(bdt)3] undergoes hydrolysis, yielding dihydroxo-bis(…

Free energy profiles for two ubiquitous damaging agents: methylation and hydroxylation of guanine in B-DNA

International audience; DNA methylation and hydroxylation are two ubiquitous reactions in DNA damage induction, yet insights are scarce concerning the free energy of activation within B-DNA. We resort to multiscale simulations to investigate the attack of a hydroxyl radical and of the primary diazonium onto a guanine embedded in a solvated dodecamer. Reaction free energy profiles characterize two strongly exergonic processes, yet allow unprecedented quantification of the barrier towards this damage reaction, not higher than 6 kcal mol−1 and sometimes inexistent, and of the exergonicities. In the case of the [G(C8)-OH]˙ intermediate, we challenge the functional dependence of such simulations…

QM/MM Determination of Kinetic Isotope Effects for COMT-Catalyzed Methyl Transfer Does Not Support Compression Hypothesis

Secondary alpha-D3 kinetic isotope effects calculated by the hybrid AM1/TIP3P/CHARMM method for the reaction of S-adenosylmethionine with catecholate anion in aqueous solution and catalyzed by rat liver catechol O-methyltransferase at 298 K are 0.94 and 0.85, respectively, in good accord with experiment. The large inverse effect for the enzymatic reaction is not due to compression but arises from significant increases in the stretching and bending force constants involving the isotopically substituted atoms of the transferring methyl group as between the reactant complex and the transition structure, larger than for the reaction in water.

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…

A quantum mechanics/molecular mechanics study of the protein-ligand interaction for inhibitors of HIV-1 integrase.

Human immunodeficiency virus type-1 integrase (HIV-1 IN) is an essential enzyme for effective viral replication. Diketo acids such as L-731,988 and S-1360 are potent and selective inhibitors of HIV-1 IN. In this study, we used molecular dynamics simulations, within the hybrid quantum mechanics/molecular mechanics (QM/MM) approach, to determine the protein-ligand interaction energy between HIV-1 IN and L-731,988 and 10 of its derivatives and analogues. This hybrid methodology has the advantage that it includes quantum effects such as ligand polarisation upon binding, which can be very important when highly polarisable groups are embedded in anisotropic environments, as for example in metal-c…

Amino Acid Chemistry in Solution:  Structural Study and Vibrational Dynamics of Glutamine in Solution. An ab Initio Reaction Field Model

The self-consistent reaction field (SCRF) theory was used to study structural and vibrational features of the amino acid L-glutamine in solution. Raman and infrared spectra of this molecule in solutions of H 2O and D2O were recorded and measured. The bands were firstly assigned on the basis of the isotopic shifts. An ab initio quadratic force field at the 6-31 +G* level was achieved. The calculation simulated a polar solvent by placing the molecule in an ellipsoidal cavity surrounded by a continuum dielectric. The theoretical results, in terms of structural parameters, vibrational frequencies and descriptions, and infrared intensities, were in satisfactory agreement with the experimental da…

QM/MM calculations of kinetic isotope effects in the chorismate mutase active site.

Kinetic isotope effects have been computed for the Claisen rearrangement of chorismate to prephenate in aqueous solution and in the active site of chorismate mutase from B. subtilus. These included primary 13C and 18O and secondary 3H effects for substitutions at the bond-making and bond-breaking positions. The initial structures of the putative stationary points on the potential energy surface, required for the calculations of isotope effects using the CAMVIB/CAMISO programs, have been selected from hybrid QM/MM molecular dynamical simulations using the DYNAMO program. Refinement of the reactant complex and transition-state structures has been carried out by means of AM1/CHARMM24/TIP3P cal…

FUNDAMENTAL PRINCIPLES GOVERNING SOLVENTS USE

Mechanistic insights into the phosphoryl transfer reaction in cyclin-dependent kinase 2: a QM/MM study

AbstractCyclin-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 show that the base-assisted mechanism is preferred over the substrat…

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…

Application of Grote-Hynes theory to the reaction catalyzed by thymidylate synthase.

A theoretical study of dynamic effects on the rate-limiting step of the thymidylate synthase catalyzed reaction has been carried out by means of Grote-Hynes theory, successfully predicting the values of the recrossing effects for a chemical reaction that involves the transfer of a classical light particle. The transmission coefficients, obtained at 278, 293, 303, and 313 K, are almost invariant and in all cases far from unity, revealing a significant coupling of the environment motions and the reaction coordinate. Nevertheless, their energetic contribution to the activation free energy represents less than 0.50 kcal/mol for each of the four tested temperatures. Calculation of the transmissi…

Do zwitterionic species exist in the non-enzymatic peptide bond formation?

The use of proper computational methods and models has allowed answering the controversial question of whether zwitterionic species exist in the mechanism of peptide bond synthesis in aqueous solution. In fact, the different conformations of zwitterionic species open the door to different mechanistic paths.

A Computational Study of the Protein-Ligand Interactions in CDK2 Inhibitors: Using Quantum Mechanics/Molecular Mechanics Interaction Energy as a Predictor of the Biological Activity

ABSTRACT: We report a combined quantum mechanics/molecular mechanics (QM/MM) study to determine the protein-ligand interaction energy between CDK2 (cyclin-dependent kinase 2) and five inhibitors with the N2 -substituted 6-cyclohexylmethoxypurine scaffold. The computational results in this work show that the QM/MM interaction energy is strongly correlated to the biological activity and can be used as a predictor, at least within a family of substrates. A detailed analysis of the protein-ligand structures obtained from molecular dynamics simulations shows specific interactions within the active site that, in some cases, have not been reported before to our knowledge. The computed interaction …

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…

The transition state and cognate concepts

Abstract This review aims firstly to clarify the meanings of key terms and concepts associated with the idea of the transition state, as developed by theoreticians and applied by experimentalist, and secondly to provide an update to the meaning and significance of the transition state in an era when computational simulation, in which complexity is being increasingly incorporated, is commonly employed as a means by which to bridge the realms of theory and experiment. The relationship between the transition state and the potential-energy surface for an elementary reaction is explored, with discussion of the following terms: saddle point, minimum-energy reaction path, reaction coordinate, acti…

Hybrid Quantum Mechanics/Molecular Mechanics Simulations with Two-Dimensional Interpolated Corrections:  Application to Enzymatic Processes

Hybrid quantum mechanics/molecular mechanics (QM/MM) techniques are widely used to study chemical reactions in large systems. Because of the computational cost associated with the high dimensionality of these systems, the quantum description is usually restricted to low-level methods, such as semiempirical Hamiltonians. In some cases, the description obtained at this computational level is quite poor and corrections must be considered. We here propose a simple but efficient way to include higher-level corrections to be used in potential energy surface explorations and in the calculation of potentials of mean force. We evaluate a correction energy term as the difference between a high-level …

Mechanism of sulfur transfer across protein-protein interfaces: The cysteine desulfurase model system

CsdA cysteine desulfurase (the sulfur donor) and the CsdE sulfur acceptor are involved in biological sulfur trafficking and in iron-sulfur cluster assembly in the model bacterium Escherichia coli. CsdA and CsdE form a stable complex through a polar interface that includes CsdA Cys328 and CsdE Cys61, the two residues known to be involved in the sulfur transfer reaction. Although mechanisms for the transfer of a sulfur moiety across protein-protein interfaces have been proposed based on the IscS-IscU and IscS-TusA structures, the flexibility of the catalytic cysteine loops involved has precluded a high resolution view of the active-site geometry and chemical environment for sulfur transfer. H…

Internal rotation of 1,2-dichloroethane in haloalkane dehalogenase. A test case for analyzing electrostatic effects in enzymes

1,2-Dichoroethane (DCE) is a prototypical molecule for studying electrostatic solvent effects on molecular conformation as far as rotation around the carbon−carbon bond notably changes the electric...

Molecular Mechanism of Inhibition of DNA Methylation by Zebularine

In this work, we have analyzed the molecular mechanism of inhibition of a C5-DNA methyltransferase by zebularine using classical and QM/MM simulations. We found that the reaction proceeds with the addition of an unprotonated cysteine to the C6 position of the ring followed by methyl transfer to the C5 position. However, while the first step is reversible and presents a moderate free-energy barrier, the second step presents a large free-energy barrier, preventing the formation of the methylated complex. This mechanistic proposal agrees with recent experimental observations that point to the formation of a reversible covalent complex between DNA containing zebularine and methyltransferases. T…

GEPOL: An improved description of molecular surfaces. III. A new algorithm for the computation of a solvent-excluding surface

To understand and calculate the interactions of a solute with a solvent, a good method of computing the molecular surface is needed. Three kinds of surfaces may be used: the van der Waals Surface, the Accessible Surface, and the Molecular Surface. The latter is redefined in this article as the Solvent-Excluding Surface. The new algorithm for computing the Solvent-Excluding Surface included in the GEPOL93 program is described. GEPOL93 follows the same concept as former versions of GEPOL but with a full new algorithm. Thus, it computes the Solvent-Excluding Surface by filling the spaces not accessible to the solvent with a set of new spheres. The computation is controlled by three parameters:…

Exploring Chemical Reactivity in Enzyme Catalyzed Processes Using QM/MM Methods: An Application to Dihydrofolate Reductase

Enzymes are the catalysts used by living organisms to accelerate chemical processes under physiological conditions. In this chapter, we illustrate the current view about the origin of their extraordinary rate enhancement based on molecular simulations and, in particular, on methods based on the combination of Quantum Mechanics and Molecular Mechanics potentials which provide a solution to treat the chemical reactivity of these large and complex molecular systems. Computational studies on Dihydrofolate Reductase have been selected as a conductor wire to present the evolution and difficulties to model chemical reactivity in enzymes. The results discussed here show that experimental observatio…

Modeling methods for studying post-translational and transcriptional modifying enzymes.

Biological catalysis is a complex chemical process that involves not only electronic reorganization in the substrate but also the reorganization of the catalyst. This complexity is even larger in the case of post-transcriptional and post-translational modifications which may involve the interaction between two biomacromolecules. However, the development over the past decades of new computational methods and strategies is offering a detailed molecular picture of the catalytic event and a deep understanding of the mechanisms of chemical reactions in biological environments. Here we review the efforts made in the last years to model catalysis in post-transcriptional and post-translational proc…

A quantum mechanics/molecular mechanics study of the acylation reaction of TEM1 β-lactamase and penicillanate

The acylation step in β-lactamase catalyzed hydrolysis of β-lactams has been explored by means of a quantum mechanics/molecular mechanics approach (AM1/CHARMM). The TEM1 enzyme, a class A β-lactamase, and the penicillanate constitute the system employed in our study. The entire molecular system is divided into a quantum and a classical region: the quantum part is composed by the substrate, the serine Ser70 and the essential moieties of key active site residues, Lys73, Ser130 and Glu166, as well as a water molecule present in the active site region, while the classical part is formed by the remaining residues and structural waters of the enzyme. In particular, the sequence of steps proposed …

A simple model for barrier frequencies for enzymatic reactions.

We present a simple model to rationalize the effects of environment on the reaction barrier frequencies derived from free energy profiles. These frequencies are relevant in deviations of a rate constant from its transition state theory value and in determining which environmental dynamics participate in the reaction. In particular, this simple model can be used to understand the changes in the reaction barrier frequencies of an enzymatic catalyzed reaction and the corresponding uncatalyzed process in aqueous solution, a change which has implications for dynamical environmental effects on the enzymatic reaction. Two possible cases are analyzed, in which the polarity (charge separation/locali…

Protein Conformational Landscapes and Catalysis. Influence of Active Site Conformations in the Reaction Catalyzed by L-Lactate Dehydrogenase

In the past decade, L-Lactate Dehydrogenase (LDH) has become an extremely useful marker in both clinical diagnosis and in monitoring the course of many human diseases. It has been assumed since the 1980s that the full catalytic process of LDH starts with the binding of the cofactor and the substrate followed by the enclosure of the active site by a mobile loop of the protein before the reaction takes place. In this paper, we show that the chemical step of the LDH-catalyzed reaction can proceed within the open loop conformation, and the different reactivity of the different protein conformations would be in agreement with the broad range of rate constants measured in single-molecule spectrom…

GEPOL: An improved description of molecular surfaces II. Computing the molecular area and volume

The algorithm used by the program GEPOL for a finer description of molecular surface (for a fast calculation of molecular area and volume and for an efficient selection of sampling points) is presented in detail. Different types of surfaces such as van der Waals and Richard's molecular surfaces can be computed. As we described in the first article (J.L. Pascual-Ahuir and E. Silla, J. Comp. Chem., 11, 1047(1990)), GEPOL begins by building a set of spherical surfaces which fill the space which is not solvent accessible. In this second article, a triangular tessellation approach to select the parts of these spherical surfaces which form the molecular surface is described. By using a data coded…

Conformational equilibrium of chorismate. A QM/MM theoretical study combining statistical simulations and geometry optimisations in gas phase and in aqueous solution

We report a theoretical study on the conformational equilibrium of chorismate that precedes its rearrangement to prephenate, an important enzyme-catalyzed reaction. In first place we show that the usual classification of chorismate conformers based on the relative position of the hydroxyl and ether bridge, pseudo-diaxial and pseudo-diequatorial, is not the only relevant factor from the point of view of the a posteriori rearrangement. Here we also analyse another complementary geometrical classification based on the interatomic distance between the carbon atoms to be bounded. Using the umbrella sampling approach and this distance as distinguished internal reaction coordinate, the gas phase A…

A Microscopic Description of SARS-CoV-2 Main Protease Inhibition with Michael Acceptors. Strategies for Improving Inhibitors Design

The irreversible inhibition of the main protease of SARS-CoV-2 by a Michael acceptor known as N3 has been investigated using multiscale methods. The noncovalent enzyme–inhibitor complex was simulated using classical molecular dynamics techniques and the pose of the inhibitor in the active site was compared to that of the natural substrate, a peptide containing the Gln–Ser scissile bond. The formation of the covalent enzyme–inhibitor complex was then simulated using hybrid QM/MM free energy methods. After binding, the reaction mechanism was found to be composed of two steps: (i) the activation of the catalytic dyad (Cys145 and His41) to form an ion pair and (ii) a Michael addition where the …

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 motion…

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 …

Decision letter: Structure and mechanism of a phage-encoded SAM lyase revises catalytic function of enzyme family

Origin of Enzymatic Kinetic Isotope Effects in Human Purine Nucleoside Phosphorylase

Here we report a study of the effect of heavy isotope labeling on the reaction catalyzed by human purine nucleoside phosphorylase (hPNP) to elucidate the origin of its catalytic effect and of the enzymatic kinetic isotope effect (EKIE). Using quantum mechanical and molecular mechanical (QM/MM) molecular dynamics (MD) simulations, we study the mechanism of the hPNP enzyme and the dynamic effects by means of the calculation of the recrossing transmission coefficient. A free energy surface (FES), as a function of both a chemical and an environmental coordinate, is obtained to show the role of the environment on the chemical reaction. Analysis of reactive and nonreactive trajectories allows us …

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 …

Improving the QM/MM Description of Chemical Processes:  A Dual Level Strategy To Explore the Potential Energy Surface in Very Large Systems.

Potential energy surfaces are fundamental tools for the analysis of reaction mechanisms. The accuracy of these surfaces for reactions in very large systems is often limited by the size of the system even if hybrid quantum mechanics/molecular mechanics (QM/MM) strategies are employed. The large number of degrees of freedom of the system requires hundreds or even thousands of optimization steps to reach convergence. Reactions in condensed media (such as enzymes or solutions) are thus usually restricted to be analyzed using low level quantum mechanical methods, thus introducing a source of error in the description of the QM region. In this paper, an alternative method is proposed, coupled to t…

Computing kinetic isotope effects for chorismate mutase with high accuracy. A new DFT/MM strategy.

A novel procedure has been applied to compute experimentally unobserved intrinsic kinetic isotope effects upon the rearrangement of chorismate to prephenate catalyzed by B. subtilis chorismate mutase. In this modified QM/MM approach, the "low-level" QM description of the quantum region is corrected during the optimization procedure by means of a "high-level" calculation in vacuo, keeping the QM-MM interaction contribution at a quantum "low-level". This allows computation of energies, gradients, and Hessians including the polarization of the QM subsystem and its interaction with the MM environment, both terms calculated using the low-level method at a reasonable computational cost. New infor…

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…

Hybrid Schemes Based on Quantum Mechanics/Molecular Mechanics Simulations

The development of characterization techniques, advanced synthesis methods, as well as molecular modeling has transformed the study of systems in a well-established research field. The current research challenges in biocatalysis and biotransformation evolve around enzyme discovery, design, and optimization. How can we find or create enzymes that catalyze important synthetic reactions, even reactions that may not exist in nature? What is the source of enzyme catalytic power? To answer these and other related questions, the standard strategies have evolved from trial-and-error methodologies based on chemical knowledge, accumulated experience, and common sense into a clearly multidisciplinary …

Computational simulations on the binding and reactivity of a nitrile inhibitor of the SARS-CoV-2 main protease.

We present a detailed computational analysis of the binding mode and reactivity of the novel oral inhibitor PF-07321332 developed against the SARS-CoV-2 3CL protease. Alchemical free energy calculations suggest that positions P3 and P4 could be susceptible to improvement in order to get a larger binding strength. QM/MM simulations unveil the reaction mechanism for covalent inhibition, showing that the nitrile warhead facilitates the recruitment of a water molecule for the proton transfer step.

Dependence of enzyme reaction mechanism on protonation state of titratable residues and QM level description: lactate dehydrogenase

We have studied the dependence of the chemical reaction mechanism of L-lactate dehydrogenase (LDH) on the protonation state of titratable residues and on the level of the quantum mechanical (QM) description by means of hybrid quantum-mechanical/molecular-mechanical (QM/MM) methods; this methodology has allowed clarification of the timing of the hydride transfer and proton transfer components that hitherto had not been possible to state definitively. Ferrer Castillo, Silvia, Silvia.Ferrer@uv.es, Silla Santos, Estanislao, Estanislao.Silla@uv.es ; Tuñon Garcia de Vicuña, Ignacio Nilo, Ignacio.Tunon@uv.es

Towards a Rational Design of Antibody Catalysts through Computational Chemistry

Can hydrophobic interactions be correctly reproduced by the continuum models?

The ability of the continuum models to describe hydrophobic interactions is investigated. In this work we have studied the interactions between two methane molecules in aqueous solution by means of a continuum model. The resulting potential of mean force is in good agreement with those obtained using Monte Carlo and molecular dynamics techniques. The three energy contributions appearing in the continuum energy partition (electrostatic, dispersion−repulsion, and cavitation) have been analyzed. The cavitation free energy plays the most important role of the three, determining the existence of an energy barrier between the contact minimum and the separated methane monomers. This barrier, which…

Theoretical site-directed mutagenesis: Asp168Ala mutant of lactate dehydrogenase

Molecular simulations based on the use of hybrid quantum mechanics/molecular mechanics methods are able to provide detailed information about the complex enzymatic reactions and the consequences of specific mutations on the activity of the enzyme. In this work, the reduction of pyruvate to lactate catalysed by wild-type and Asp168Ala mutant lactate dehydrogenase (LDH) has been studied by means of simulations using a very flexible molecular model consisting of the full tetramer of the enzyme, together with the cofactor NADH, the substrate and solvent water molecules. Our results indicate that the Asp168Ala mutation provokes a shift in the p K a value of Glu199 that becomes unprotonated at n…

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…

Computational strategies for the design of new enzymatic functions

In this contribution, recent developments in the design of biocatalysts are reviewed with particular emphasis in the de novo strategy. Studies based on three different reactions, Kemp elimination, Diels–Alder and Retro-Aldolase, are used to illustrate different success achieved during the last years. Finally, a section is devoted to the particular case of designed metalloenzymes. As a general conclusion, the interplay between new and more sophisticated engineering protocols and computational methods, based on molecular dynamics simulations with Quantum Mechanics/Molecular Mechanics potentials and fully flexible models, seems to constitute the bed rock for present and future successful desig…

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

Are there dynamical effects in enzyme catalysis? Some thoughts concerning the enzymatic chemical step.

Highlights • The possible role of enzymatic reaction dynamical effects is examined. • Solution reactions usefully inform the issue of dynamical effects in enzymes. • Division into regions containing and away from the transition state is important. • Motions in passage to/from the transition state need not lead to dynamical effects. • Transition State Theory is usually a reasonable description of enzyme kinetics.

A quantum mechanics-molecular mechanics study of dissociative electron transfer : The methylchloride radical anion in aqueous solution

The dissociative electron transfer reaction CH3Cl+e−→CH3•+Cl− in aqueous solution is studied by using a QM/MM method. In this work the quantum subsystem (a methylchloride molecule plus an electron) is described using density functional theory while the solvent (300 water molecules) is described using the TIP3P classical potential. By means of molecular dynamics simulations and the thermodynamic integration technique we obtained the potential of mean force (PMF) for the carbon–chlorine bond dissociation of the neutral and radical anion species. Combining these two free energy curves we found a quadratic dependence of the activation free energy on the reaction free energy in agreement with Ma…

Unraveling the SARS-CoV-2 Main Protease Mechanism Using Multiscale Methods

We present a detailed theoretical analysis of the reaction mechanism of proteolysis catalyzed by the main protease of SARS-CoV-2. Using multiscale simulation methods, we have characterized the interactions established by a peptidic substrate in the active site, and then we have explored the free energy landscape associated with the acylation and deacylation steps of the proteolysis reaction, characterizing the transition states of the process. Our mechanistic proposals can explain most of the experimental observations made on the highly similar ortholog protease of SARS-CoV. We point to some key interactions that may facilitate the acylation process and thus can be crucial in the design of …

Are Heme-Dependent Enzymes Always Using a Redox Mechanism? A Theoretical Study of the Kemp Elimination Catalyzed by a Promiscuous Aldoxime Dehydratase

The design of biocatalysts is a goal to improve the rate, selectivity and environmental friendship of chemical processes in biotechnology. In this regard, the use of computational techniques has provided valuable assistance in the design of enzymes with remarkable catalytic activity. In this paper, hybrid QM/MM simulations have allowed getting an insight into the mechanism of a promiscuous aldoxime dehydratase (OxdA) for the Kemp elimination. We first demonstrate that, based on the use of linear response approximation (LRA) methods, the lowest energy electronic state of the benzisoxazole placed in the active sit of OxdA corresponds to a singlet state, being the triplet and the quintet state…

Translocation of enzymes into a mesoporous MOF for enhanced catalytic activity under extreme conditions

Translocation of protease into mesoporous MIL-101-NH2 results in enhanced catalytic activity, excellent recyclability and tolerance to competing enzymes.

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.

Heavy enzymes and the rational redesign of protein catalysts

Abstract An unsolved mystery in biology concerns the link between enzyme catalysis and protein motions. Comparison between isotopically labelled “heavy” dihydrofolate reductases and their natural‐abundance counterparts has suggested that the coupling of protein motions to the chemistry of the catalysed reaction is minimised in the case of hydride transfer. In alcohol dehydrogenases, unnatural, bulky substrates that induce additional electrostatic rearrangements of the active site enhance coupled motions. This finding could provide a new route to engineering enzymes with altered substrate specificity, because amino acid residues responsible for dynamic coupling with a given substrate present…

ChemInform Abstract: Computational Design of Biological Catalysts

The purpose of this tutorial review is to illustrate the way to design new and powerful catalysts. The first possibility to get a biological catalyst for a given chemical process is to use existing enzymes that catalyze related reactions. The second possibility is the use of immune systems that recognize stable molecules resembling the transition structure of the target reaction. We finally show how computational techniques are able to provide an enormous quantity of information, providing clues to guide the development of new biological catalysts.

Understanding the different activities of highly promiscuous MbtI by computational methods

Salicylate synthase from Mycobacterium tuberculosis, MbtI, is a highly promiscuous Mg(2+) dependent enzyme with up to four distinct activities detected in vitro: isochorismate synthase (IS), isochorismate pyruvate lyase (IPL), salicylate synthase (SS) and chorismate mutase (CM). In this paper, Molecular Dynamic (MD) simulations employing hybrid quantum mechanics/molecular mechanics (QM/MM) potentials have been carried out to get a detailed knowledge of the IS and the IPL activities at the molecular level. According to our simulations, the architecture of the MbtI active site allows catalyzing the two reactions: the isochorismate formation, by means of a stepwise mechanism, and the salicylat…

Multiscale Simulations of SARS-CoV-2 3CL Protease Inhibition with Aldehyde Derivatives. Role of Protein and Inhibitor Conformational Changes in the Reaction Mechanism

We here investigate the mechanism of SARS-CoV-2 3CL protease inhibition by one of the most promising families of inhibitors, those containing an aldehyde group as a warhead. These compounds are covalent inhibitors that inactivate the protease, forming a stable hemithioacetal complex. Inhibitor 11a is a potent inhibitor that has been already tested in vitro and in animals. Using a combination of classical and QM/MM simulations, we determined the binding mode of the inhibitor into the active site and the preferred rotameric state of the catalytic histidine. In the noncovalent complex, the aldehyde group is accommodated into the oxyanion hole formed by the NH main-chain groups of residues 143 …

Theoretical study of phosphodiester hydrolysis in nucleotide pyrophosphatase/phosphodiesterase. Environmental effects on the reaction mechanism.

We here present a theoretical study of the alkaline hydrolysis of methyl p-nitrophenyl phosphate (MpNPP(-)) in aqueous solution and in the active site of nucleotide pyrophosphatase/phosphodiesterase (NPP). The analysis of our simulations, carried out by means of hybrid quantum mechanics/molecular mechanics (QM/MM) methods, shows that the reaction takes place through different reaction mechanisms depending on the environment. Thus, while in aqueous solution the reaction occurs by means of an A(N)D(N) mechanism, the enzymatic process takes place through a D(N)A(N) mechanism. In the first case, we found associative transition-state (TS) structures, while in the enzyme TS structures have dissoc…

Protein Flexibility and Preorganization in the Design of Enzymes. The Kemp Elimination Catalyzed by HG3.17

A recently designed enzyme, HG3.17, obtained by directed evolution, has shown a catalytic activity close to natural enzymes. Hybrid QM/MM molecular dynamics simulations for the Kemp elimination in this new enzyme have provided a deep insight into the origin of its catalytic efficiency. In this case, we have first demonstrated the presence of different conformations with significantly different reactivity. The larger reactivity is related with a better electrostatic preorganization of the environment that creates a more favorable electrostatic potential for the reaction to proceed. In HG3.17, efforts to improve the catalytic properties must be focused in possible mutations increasing the pre…

Exploration of the Activation Mechanism of the Epigenetic Regulator MLL3: A QM/MM Study

The mixed lineage leukemia 3 or MLL3 is the enzyme in charge of the writing of an epigenetic mark through the methylation of lysine 4 from the N-terminal domain of histone 3 and its deregulation has been related to several cancer lines. An interesting feature of this enzyme comes from its regulation mechanism, which involves its binding to an activating dimer before it can be catalytically functional. Once the trimer is formed, the reaction mechanism proceeds through the deprotonation of the lysine followed by the methyl-transfer reaction. Here we present a detailed exploration of the activation mechanism through a QM/MM approach focusing on both steps of the reaction, aiming to provide new…

Loss of Hyperconjugative Effects Drives Hydride Transfer during Dihydrofolate Reductase Catalysis

Hydride transfer is widespread in nature and has an essential role in applied research. However, the mechanisms of how this transformation occurs in living organisms remain a matter of vigorous debate. Here, we examined dihydrofolate reductase (DHFR), an enzyme that catalyzes hydride from C4′ of NADPH to C6 of 7,8-dihydrofolate (H2F). Despite many investigations of the mechanism of this reaction, the contribution of polarization of the π-bond of H2F in driving hydride transfer remains unclear. H2F was stereospecifically labeled with deuterium β to the reacting center, and β-deuterium kinetic isotope effects were measured. Our experimental results combined with analysis derived from QM/MM si…

A Novel Strategy to Study Electrostatic Effects in Chemical Reactions: Differences between the Role of Solvent and the Active Site of Chalcone Isomerase in a Michael Addition

The electrostatic behavior of active site residues in enzyme catalysis is quite different from that of water molecules in solution. To highlight the electrostatic differences between both environments, we propose a QM/MM strategy to study the role of the environment in chemical reactions. The novelty of the present communication is that free energy surfaces are generated by means of two distinguished reaction coordinates: a solute coordinate and the electrostatic potential created by the environment. This is applied to analyze the origin of catalysis in the transformation of a chalcone into a flavanone, a Michael addition that requires the desolvation of the nucleophile.

Why Are Some Enzymes Dimers? Flexibility and Catalysis in Thermotoga maritima Dihydrofolate Reductase

[Image: see text] Dihydrofolate reductase from Thermotoga maritima (TmDFHFR) is a dimeric thermophilic enzyme that catalyzes the hydride transfer from the cofactor NADPH to dihydrofolate less efficiently than other DHFR enzymes, such as the mesophilic analogue Escherichia coli DHFR (EcDHFR). Using QM/MM potentials, we show that the reduced catalytic efficiency of TmDHFR is most likely due to differences in the amino acid sequence that stabilize the M20 loop in an open conformation, which prevents the formation of some interactions in the transition state and increases the number of water molecules in the active site. However, dimerization provides two advantages to the thermophilic enzyme: …

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…

On the Nature of the Enzyme–Substrate Complex and the Reaction Mechanism in Human Arginase I. A Combined Molecular Dynamics and QM/MM Study

We present here a detailed theoretical analysis of L-arginine hydrolysis catalyzed by Human Arginase I (HARGI). Our study combines classical molecular dynamic simulations of different model for the...

Analysis of the Decarboxylation Step in Mammalian Histidine Decarboxylase

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 va…

Molecular Mechanism of the site-specific self-cleavage of the RNA phosphodiester backbone by a Twister Ribozyme

Published as part of the special collection of articles derived from the 10th Congress on Electronic Structure: Principles and Applications (ESPA-2016). The catalytic activity of some classes of natural RNA, named as ribozymes, has been discovered just in the past decades. In this paper, the cleavage of the RNA phosphodiester backbone has been studied in aqueous solution and in a twister ribozyme from Oryza sativa. The free energy profiles associated with a baseline substrate-assisted mechanism for the reaction in the enzyme and in solution were computed by means of free energy perturbation methods within hybrid QM/MM potentials, describing the chemical system by the M06-2× functional and t…

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…

Theoretical Study of Catalytic Efficiency of a Diels–Alderase Catalytic Antibody: An Indirect Effect Produced During the Maturation Process

The Diels–Alder reaction is one of the most important and versatile transformations available to organic chemists for the construction of complex natural products, therapeutics agents, and synthetic materials. Given the lack of efficient enzymes capable of catalyzing this kind of reaction, it is of interest to ask whether a biological catalyst could be designed from an antibody-combining site. In the present work, a theoretical study of the different behavior of a germline catalytic antibody (CA) and its matured form, 39 A-11, that catalyze a Diels–Alder reaction has been carried out. A free-energy perturbation technique based on a hybrid quantum-mechanics/molecular-mechanics scheme, togeth…

Theoretical Studies of the Self Cleavage Pistol Ribozyme Mechanism

AbstractRibozymes are huge complex biological catalysts composed of a combination of RNA and proteins. Nevertheless, there is a reduced number of small ribozymes, the self-cleavage ribozymes, that are formed just by RNA and, apparently, they existed in cells of primitive biological systems. Unveiling the details of these “fossils” enzymes can contribute not only to the understanding of the origins of life but also to the development of new simplified artificial enzymes. A computational study of the reactivity of the pistol ribozyme carried out by means of classical MD simulations and QM/MM hybrid calculations is herein presented to clarify its catalytic mechanism. Analysis of the geometries…

Amino acid chemistry in solution: structural properties and vibrational dynamics of serine using density functional theory and a continuum solvent model

A structural and vibrational study of amino acid serine in aqueous solution has been carried out using Fourier transform spectroscopies and quantum mechanical calculations. FT-IR and FT-Raman spectra of serine in H2O and D2O solutions were recorded and a general assignment of the observed bands was proposed on the basis of a zwitterionic structure for serine. Main criteria were the observed wavenumber shifts upon deuteration and previous assignments for other amino acids. A quadratic force field was computed using ab initio methodology at the 6-31+G** level and the hybrid functional B3LYP. The solvent effect was simulated by placing the serine molecule into an ellipsoidal cavity surrounded …

Modeling ?-lactam interactions in aqueous solution through combined quantum mechanics-molecular mechanics methods

Transition structure selectivity in enzyme catalysis: a QM/MM study of chorismate mutase

Two different transition structures (TSs) have been located and characterized for the chorismate conversion to prephenate in Bacillus subtilis chorismate mutase by means of hybrid quantum-mechanical/molecular-mechanical (QM/MM) calculations. GRACE software, combined with an AM1/CHARMM24/TIP3P potential, has been used involving full gradient relaxation of the position of ca. 3300 atoms. These TSs have been connected with their respective reactants and products by the intrinsic reaction coordinate (IRC) procedure carried out in the presence of the protein environment, thus obtaining for the first time a realistic enzymatic reaction path for this reaction. Similar QM/MM computational schemes h…

Unraveling the role of protein dynamics in dihydrofolate reductase catalysis

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 ((15)N, (13)C, (2)H 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 e…

A theoretical study of solvent effects on the conformational equilibria of neutral glycine in aqueous solution

Abstract In this work conformational equilibrium of neutral glycine in solution is systematically investigated by using DFT and MP2 methods combined with solvent continuum models. A systematic exploration of the potential energy surface and full geometry optimizations for several conformers have been carried out in the gas phase and aqueous solution at the MP2/6-31+G** and B3LYP/6-31+G** levels. Zero-point and thermal contributions to the free energy have been obtained at the B3LYP level. Both theoretical levels lead to very similar results, in geometrical and energetic terms, both in the gas phase and in solution. Solvent effects play an important role on the conformational equilibria of n…

Computational simulation of the lifetime of the methoxymethyl cation in water. A simple model for a glycosyl cation: when is an intermediate an intermediate?

A two-dimensional free-energy surface is constructed for transfer of the methoxymethyl cation between two water molecules. These atoms are treated quantum mechanically within a box of >1000 classical solvent water molecules, and the molecular dynamics of the whole system is considered at 300 K. This provides a simple model for glycosyl transfer in water. The best surface obtained (MPWB1K/6-31+G(d,p) corrected AMI/TIP3P) contains a shallow free-energy well corresponding to an oxacarbenium ion intermediate in a stepwise mechanism. Molecular dynamics analysis at three temperatures leads to a classical estimate of the lifetime of the methoxymethyl cation in water; when quantum corrections fo…

A Quantum Mechanic/Molecular Mechanic Study of the Wild-Type and N155S Mutant HIV-1 Integrase Complexed with Diketo Acid

Integrase (IN) is one of the three human immunodeficiency virus type 1 (HIV-1) enzymes essential for effective viral replication. Recently, mutation studies have been reported that have shown that a certain degree of viral resistance to diketo acids (DKAs) appears when some amino acid residues of the IN active site are mutated. Mutations represent a fascinating experimental challenge, and we invite theoretical simulations for the disclosure of still unexplored features of enzyme reactions. The aim of this work is to understand the molecular mechanisms of HIV-1 IN drug resistance, which will be useful for designing anti-HIV inhibitors with unique resistance profiles. In this study, we use mo…

Vibrational Dynamics of Histamine Monocation in Solution: An Experimental (FT-IR, FT-Raman) and Theoretical (SCRF-DFT) Study

spectra of histamine monocation in solution, based on the isotopic shifts and a previous vibrational study in solid state. Force field and normal coordinate calculations were computed to support these assignments. The ab initio force constants were transformed into a set of locally symmetrized internal coordinates and subsequently scaled to the experimental frequencies by using one specific and two generic scaling factors. The comparison in terms of vibrational frequencies and normal coordinate descriptions supported most of the proposed assignments. The theoretical infrared spectra for the two isotopomers on the basis of the ab initio intensities also showed a good correlation with the exp…

A DFT study on the chiral synthesis of R-phenylacetyl carbinol within the quantum chemical cluster approach

Abstract The reaction pathway leading to R-phenylacetyl carbinol within the quantum chemical cluster approach is addressed by means of density functional theory (DFT) calculations. The study includes calculation of Fukui functions, activation free energies, and potential energy surface scans, both in gas and solution phase. The protonation states of the nitrogen atoms of the pyrimidine moiety are determined. The reaction appears to be slightly exergonic (ΔG 0  = −5.6 and −4.0 kcal/mol for gas and solution phase, respectively) following a concerted synchronous mechanism having activation free energy barriers of 16.2 and 13.3 kcal/mol, in gas phase and solution phase, respectively.

Comparative computational analysis of different active site conformations and substrates in a chalcone isomerase catalyzed reaction.

Chalcone isomerase catalyzes the transformation of chalcones to flavanones. We present a computational study of the rate-limiting chemical step, an intramolecular Michael addition of a 2'-oxyanion to the alpha,beta-double bound. By using quantum mechanical/molecular mechanical hybrid methods we traced the free-energy profiles associated with the reaction of two different substrates (chalcone and 6'-deoxychalcone) in two different conformations of the active site that are described in the different crystallographic structures available. We have obtained significant differences (about 4 kcal/mol) in the free-energy barriers calculated for the two active sites. According to our results, the ac…

Insights into the inhibited form of the redox-sensitive SufE-like sulfur acceptor CsdE

17 p.-8 fig.

Intramolecular proton transfer of serine in aqueous solution. Mechanism and energetics

Serine amino acid in aqueous solution is theoretically studied at the B3PW91/6-31+G** level using a dielectric continuum solvent model. Neutral and zwitterionic structures in the gas phase and in solution are described and the proton-transfer mechanism is discussed. A neutral conformation in which the carboxyl hydrogen atom is already oriented toward the amino group seems to be the absolute energy minimum in the gas phase and the most stable neutral form in solution. The absolute energy minimum in solution is a zwitterionic form. The energy barrier for proton transfer is predicted to be very small, in particular when zero-point-energy contributions are added. Our calculations allow the dyna…

Continuum-uniform approach calculations of the solubility of hydrocarbons in water

Abstract The ransfer free energies from gas phase to water for some hydrocarbons are calculated by means of a continuum-uniform model of the solvent. For the calculation of the cavitation energy a model based on the surface tension is proposed. The calculated values are compared with the experimental free energies obtained with and without a corrective factor that accounts for the difference in the solute—solvent sizes. Good agreement between the theoretical free energies and the corrected experimental data is obtained. Our calculations seem to show that the hydrophobic effect is directly related to the molecular surface area.

Enzyme molecular mechanism as a starting point to design new inhibitors: a theoretical study of O-GlcNAcase.

O-Glycoprotein 2-acetamino-2-deoxy-β-d-glucopyranosidase (O-GlcNAcase) hydrolyzes O-linked 2-acetamido-2-deoxy-β-d-glucopyranoside (O-GlcNAc) residues from post-translationally modified serine/threonine residues of nucleocytoplasmic protein. The chemical process involves substrate-assisted catalysis, where two aspartate residues have been identified as the two key catalytic residues of O-GlcNAcase. In this report, the first step of the catalytic mechanism used by O-GlcNAcase involving substrate-assisted catalysis has been studied using a hybrid quantum mechanical/molecular mechanical (QM/MM) Molecular Dynamics (MD) calculations. The free energy profile shows that the formation of the oxazol…

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 sho…

FT-Raman and QM/MM study of the interaction between histamine and DNA

The interaction between histamine and highly polymerized calf-thymus DNA has been investigated using FT-Raman spectroscopy and the hybrid QM/MM (quantum mechanics/molecular mechanics) methodology. Raman spectra of solutions containing histamine and calf-thymus DNA, at different molar ratios, were recorded. Solutions were prepared at physiological settings of pH and ionic strength, using both natural and heavy water as the solvent. The analysis of the spectral changes on the DNA Raman spectra when adding different concentrations of histamine allowed us to identify the reactive sites of DNA and histamine, which were used to built two minor groove and one intercalated binding models. They were…

Digitalis purpurea P5 beta R2, encoding steroid 5 beta-reductase, is a novel defense-related gene involved in cardenolide biosynthesis.

The stereospecific 5 beta-reduction of progesterone is a required step for cardiac glycoside biosynthesis in foxglove plants. Recently, we have isolated the gene P5 beta R, and here we investigate the function and regulation of P5 beta R2, a new progesterone 5 beta-reductase gene from Digitalis purpurea. P5 beta R2 cDNA was isolated from a D. purpurea cDNA library and further characterized at the biochemical, structural and physiological levels. Like P5 beta R, P5 beta R2 catalyzes the 5 beta-reduction of the Delta(4) double bond of several steroids and is present in all plant organs. Under stress conditions or on treatment with chemical elicitors, P5 beta R expression does not vary, wherea…

Insights on the origin of catalysis on glycine N-methyltransferase from computational modeling.

The origin of enzyme catalysis remains a question of debate despite much intense study. We report a QM/MM theoretical study of the SN2 methyl transfer reaction catalyzed by a glycine N-methyltransferase (GNMT) and three mutants to test whether recent experimental observations of rate-constant reductions and variations in inverse secondary α-3H kinetic isotope effects (KIEs) should be attributed to changes in the methyl donor−acceptor distance (DAD): is catalysis due to a compression effect? Semiempirical (AM1) and DFT (M06-2X) methods were used to describe the QM subset of atoms, while OPLS-AA and TIP3P classical force fields were used for the protein and water molecules, respectively. The …

Computational study on hydrolysis of cefotaxime in gas phase and in aqueous solution

We are presenting a theoretical study of the hydrolysis of a β-lactam antibiotic in gas phase and in aqueous solution by means of hybrid quantum mechanics/molecular mechanics potentials. After exploring the potential energy surfaces at semiempirical and density functional theory (DFT) level, potentials of mean force have been computed for the reaction in solution with hybrid PM3/TIP3P calculations and corrections with the B3LYP and M06-2X functionals. Inclusion of the full molecule of the antibiotic, Cefotaxime, in the gas phase molecular model has been demonstrated to be crucial since its carboxylate group can activate a nucleophilic water molecule. Moreover, the flexibility of the substra…

Mechanistic study of the biosynthesis of R-phenylcarbinol by acetohydroxyacid synthase enzyme using hybrid quantum mechanics/molecular mechanics simulations

Abstract The biosynthesis of R-phenylacetylcarbinol (R-PAC) by the acetohydroxy acid synthase, (AHAS) is addressed by molecular dynamics simulations (MD), hybrid quantum mechanics/molecular mechanics (QM/MM), and QM/MM free energy calculations. The results show the reaction starts with the nucleophilic attack of the C2α atom of the HEThDP intermediate on the Cβ atom of the carbonyl group of benzaldehyde substrate via the formation of a transition state (TS1) with the HEThDP intermediate under 4′-aminopyrimidium (APH+) form. The calculated activation free energy for this step is 17.4 kcal mol−1 at 27 °C. From this point, the reaction continues with the abstraction of Hβ atom of the HEThDP in…

Coupling of the guanosine glycosidic bond conformation and the ribonucleotide cleavage reaction: implications for barnase catalysis.

To examine the possible relationship of guanine-dependent GpA conformations with ribonucleotide cleavage, two potential of mean force (PMF) calculations were performed in aqueous solution. In the first calculation, the guanosine glycosidic (Gχ) angle was used as the reaction coordinate, and computations were performed on two GpA ionic species: protonated (neutral) or deprotonated (negatively charged) guanosine ribose O2 ′. Similar energetic profiles featuring two minima corresponding to the anti and syn Gχ regions were obtained for both ionic forms. For both forms the anti conformation was more stable than the syn, and barriers of ∼4 kcal/mol were obtained for the anti → syn transition. Str…

Convergence of Theory and Experiment on the Role of Preorganization, Quantum Tunneling, and Enzyme Motions into Flavoenzyme-Catalyzed Hydride Transfer

Hydride transfer is one of the most common reactions catalyzed by enzymatic systems, and it has become an object of study because of possible significant quantum tunneling effects. In the present work, we provide a combination of theoretical QM/MM simulations and experimental measurements of the rate constants and kinetic isotopic effects (KIEs) for the hydride transfer reaction catalyzed by morphinone reductase, MR. Quantum mechanical tunneling coefficients, computed in the framework of variational transition-state theory, play a significant role in this reaction, reaching values of 23.8 ± 5.5 for the lightest isotopologue—one of the largest values reported for enzymatic systems. This pred…

Targeting the JAK/STAT Pathway: A Combined Ligand- and Target-Based Approach

Janus kinases (JAKs) are a family of proinflammatory enzymes able to mediate the immune responses and the inflammatory cascade by modulating multiple cytokine expressions as well as various growth factors. In the present study, the inhibition of the JAK-signal transducer and activator of transcription (STAT) signaling pathway is explored as a potential strategy for treating autoimmune and inflammatory disorders. A computationally driven approach aimed at identifying novel JAK inhibitors based on molecular topology, docking, and molecular dynamics simulations was carried out. For the best candidates selected, the inhibitory activity against JAK2 was evaluated in vitro. Two hit compounds with…

On the tautomerization process of glycine in aqueous solution

The experimental activation energy for the tautomerization of glycine zwitterion neutral form has been reported to be 14.6 kcalrmol. It has been generally assumed that this energy barrier is needed for proton transfer to occur. However, previous theoretical results do not support this interpretation. In the present work, we examine this question using density functional calculations, extended basis sets and a polarizable continuum solvent model. Our results suggest that the limiting step for the tautomerization process corresponds basically to H-atom reorientation in the -COOH group. This could be a general feature in the tautomerization of amino acids. q 2000 Elsevier Science B.V. All righ…

Catalysis in glycine N-methyltransferase: testing the electrostatic stabilization and compression hypothesis.

Glycine N-methyltransferase (GNMT) is an S-adenosyl-l-methionine dependent enzyme that catalyzes glycine transformation to sarcosine. Here, we present a hybrid quantum mechanics/molecular mechanics (QM/MM) computational study of the reaction compared to the counterpart process in water. The process takes place through an SN2 mechanism in both media with a transition state in which the transferring methyl group is placed in between the donor (SAM) and the acceptor (the amine group of glycine). Comparative analysis of structural, electrostatic, and electronic characteristics of the in-solution and enzymatic transition states allows us to get a deeper insight into the origins of the enzyme's c…

Correlation effects in proton transfer reactions in solution

Abstract The effects of correlation energy on proton transfer reactions in solution for [H 2 OHOH 2 ] + and [NH 3 HH 2 O] + systems have been studied. Solvent effects have been represented by means of a continuum model. The proton transfer energy profiles for fixed proton donor-proton acceptor distances have been obtained in the gas phase and in solution, both at the HF/6-311G ∗∗ and MP2/6-311G ∗∗ //HF levels of theory. Differences between the correlation energies calculated in the gas phase and in solution are negligible, showing that solvent effects can be correctly described for these proton transfer processes at the HF level.

Theoretical Study of the Catalytic Mechanism of DNA-(N4-Cytosine)-Methyltransferase from the Bacterium Proteus vulgaris

In this paper the reaction mechanism for methylation of cytosine at the exocyclic N4 position catalyzed by M.PvuII has been explored by means of hybrid quantum mechanics/molecular mechanics (QM/MM) methods. A reaction model was prepared by placing a single cytosine base in the active site of the enzyme. In this model the exocyclic amino group of the base establishes hydrogen bond interactions with the hydroxyl oxygen atom of Ser53 and the carbonyl oxygen atom of Pro54. The reaction mechanism involves a direct methyl transfer from AdoMet to the N4 atom and a proton transfer from this atom to Ser53, which in turn transfers a proton to Asp96. Different timings for the proton transfers and meth…

Role of Protein Flexibility in Enzymatic Catalysis:  Quantum Mechanical−Molecular Mechanical Study of the Deacylation Reaction in Class A β-Lactamases

We present a theoretical study of a mechanism for the hydrolysis of the acyl-enzyme complex formed by a class A beta-lactamase (TEM1) and an antibiotic (penicillanate), as a part of the process of antibiotic's inactivation by this type of enzymes. In the presented mechanism the carboxylate group of a particular residue (Glu166) activates a water molecule, accepting one of its protons, and afterward transfers this proton directly to the acylated serine residue (Ser70). In our study we employed a quantum mechanics (AM1)-molecular mechanics partition scheme (QM/MM) where all the atoms of the system were allowed to relax. For this purpose we used the GRACE procedure in which part of the system …

Dynamic Effects on Reaction Rates in a Michael Addition Catalyzed by Chalcone Isomerase. Beyond the Frozen Environment Approach

We present a detailed microscopic study of the dynamics of the Michael addition reaction leading from 6'-deoxychalcone to the corresponding flavanone. The reaction dynamics are analyzed for both the uncatalyzed reaction in aqueous solution and the reaction catalyzed by Chalcone Isomerase. By means of rare event simulations of trajectories started at the transition state, we have computed the transmission coefficients, obtaining 0.76 +/- 0.04 and 0.87 +/- 0.03, in water and in the enzyme, respectively. According to these simulations, the Michael addition can be seen as a formation of a new intramolecular carbon-oxygen bond accompanied by a charge transfer essentially taking place from the nu…

Role of Solvent on Nonenzymatic Peptide Bond Formation Mechanisms and Kinetic Isotope Effects

Based on the hypothesis that similar mechanisms are involved in the peptide bond formation in aqueous solution and in the ribosome, the aminolysis of esters in aqueous solution has been the subject of numerous studies as the reference reaction for the catalyzed process. The mechanisms proposed in the literature have been explored in the present paper by hybrid QM/MM molecular dynamics simulations. The free energy profiles have been computed with the QM region of the system described at semiempirical AM1 level and by DFT within the M06-2X functional. According to the results, the formation of adduct zwitterion species is a preliminary step required for all possible mechanisms. Then, from dif…

Computational Modeling of Biological Systems: The LDH Story

Lactate dehydrogenases, LDH, catalyzed reaction has been used in this chapter as a conductor wire to present the evolution and difficulties on computing methods to model chemical reactions in enzymes, since the early calculations based at semiempirical level carried out in gas phase to the recent sophisticated simulations based on hybrid Quantum Mechanical/Molecular Mechanics Dynamics (QM/MM MD) schemes. LDH catalyzes the reversible transformation of pyruvate into lactate. The chemical step consists in a hydride and a proton transfer from the cofactor (NADH) and a protonated histidine (His195), respectively. This fact has generated a lot of controversy about the timing of both transfers in …

Revealing the Origin of the Efficiency of the De Novo Designed Kemp Eliminase HG-3.17 by Comparison with the Former Developed HG-3

The design of new biocatalysts is a goal in biotechnology to improve the rate, selectivity and environmental impact of industrial chemical processes. In this regard, the use of computational techniques has provided valuable assistance in the design of new enzymes with remarkable catalytic activity. In this paper, hybrid QM/MM molecular dynamics simulations have allowed insights to be gained on the origin of the limited efficiency of a computationally designed enzyme for the Kemp elimination; the HG-3. Comparison of results derived from this enzyme with those of a more evolved protein containing additional point mutations, HG-3.17, rendered important information that should be taken into acc…

Enzyme Promiscuity in Enolase Superfamily. Theoretical Study of o-Succinylbenzoate Synthase Using QM/MM Methods

The promiscuous activity of the enzyme o-succinylbenzoate synthase (OSBS) from the actinobacteria Amycolatopsis is investigated by means of QM/MM methods, using both density functional theory and semiempirical Hamiltonians. This enzyme catalyzes not only the dehydration of 2-succinyl-6R-hydroxy-2,4-cyclohexadiene-1R-carboxylate but also catalyzes racemization of different acylamino acids, with N-succinyl-R-phenylglycine being the best substrate. We investigated the molecular mechanisms for both reactions exploring the potential energy surface. Then, molecular dynamics simulations were performed to obtain the free energy profiles and the averaged interaction energies of enzymatic residues wi…

Aminoacid zwitterions in solution : Geometric, energetic, and vibrational analysis using density functional theory-continuum model calculations

Glycine and alanine aminoacids chemistry in solution is explored using a hybrid three parameters density functional (B3PW91) together with a continuum model. Geometries, energies, and vibrational spectra of glycine and alanine zwitterions are studied at the B3PW91/6-31+G∗∗ level and the results compared with those obtained at the HF and MP2/6-31+G∗∗ levels. Solvents effects are incorporated by means of an ellipsoidal cavity model with a multipolar expansion (up to sixth order) of the solute’s electrostatic potential. Our results confirm the validity of the B3PW91 functional for studying aminoacid chemistry in solution. Taking into account the more favorable scaling behavior of density funct…

Theoretical rotational constants of MeCnN species

Abstract By means of SCF HF “ab initio” calculations with STO-3G and 6-31G basis sets, the geometric parameters of methylcyanopolyynes (MeCnN n=3, 5, 7 and 9) have been obtained. B0=0.3748 GHz for MeC7N and B0=0.2708 GHz for MeC9N, with a STO-3G basis set, were obtained. Both species, unstable in the laboratory, are still undetected in the interstellar medium, although their existence is very probable.

Structural and vibrational study of the tautomerism of histamine free-base in solution.

Infrared and Raman spectroscopy in H(2)O and D(2)O and quantum Density Functional calculations were used to determine the structure of histamine free-base in aqueous solution. A quantum mechanical study of the tautomeric equilibrium of histamine free-base in solution was performed at the 6-311G level. Electronic correlation energies were included by using the hybrid functional B3LYP. The solvent was simulated as a continuum characterized by a dielectric constant, and the quantum system (solute) was placed in an ellipsoidal cavity. Solute-solvent electrostatic interaction was calculated by means a multipolar moment expansion introduced in the Hamiltonian. Four relevant histamine conformation…

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…

Inhibition Mechanism of SARS‐CoV‐2 Main Protease with Ketone‐Based Inhibitors Unveiled by Multiscale Simulations: Insights for Improved Designs**

Abstract We present the results of classical and QM/MM simulations for the inhibition of SARS‐CoV‐2 3CL protease by a hydroxymethylketone inhibitor, PF‐00835231. In the noncovalent complex the carbonyl oxygen atom of the warhead is placed in the oxyanion hole formed by residues 143 to 145, while P1–P3 groups are accommodated in the active site with interactions similar to those observed for the peptide substrate. According to alchemical free energy calculations, the P1′ hydroxymethyl group also contributes to the binding free energy. Covalent inhibition of the enzyme is triggered by the proton transfer from Cys145 to His41. This step is followed by the nucleophilic attack of the Sγ atom on …

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…

Molecular dynamics simulations of elementary chemical processes in liquid water using combined density functional and molecular mechanics potentials. II. Charge separation processes

A new approach to carry out molecular dynamics simulations of chemical reactions in solution using combined density functional theory/molecular mechanics potentials is presented. We focus our attention on the analysis of reactive trajectories, dynamic solvent effects and transmission coefficient rather than on the evaluation of free energy which is another important topic that will be examined elsewhere. In a previous paper we have described the generalities of this hybrid molecular dynamics method and it has been employed to investigate low energy barrier proton transfer process in water. The study of processes with activation energies larger than a few kT requires the use of specific tech…

Computer-Aided Rational Design of Catalytic Antibodies: The 1F7 Case.

Molecular dynamics simulations of elementary chemical processes in liquid water using combined density functional and molecular mechanics potentials. I. Proton transfer in strongly H-bonded complexes

The first molecular dynamics (MD) simulation of a chemical process in solution with an ab initio description of the reactant species and a classical representation of the solvent is presented. We study the dynamics of proton (deuterium) transfer in strongly hydrogen-bonded systems characterized by an energy surface presenting a double well separated by a low activation barrier. We have chosen the hydroxyl-water complex in liquid water to analyze the coupling between the reactive system and the environment. The proton is transferred from one well to the other with a frequency close to 1 ps−1 which is comparable to the low-frequency band associated to hindered translations, diffusional transl…

Dynamics of the excited-state hydrogen transfer in a (dG)·(dC) homopolymer: intrinsic photostability of DNA

Multiscale molecular dynamics simulations reveal out-of-plane distortions that favour DNA photostability. A novel photostability mechanism involving four proton transfers and triggered by a nearby Na+ ion is also unveiled.

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>

Does glycosyl transfer involve an oxacarbenium intermediate? Computational simulation of the lifetime of the methoxymethyl cation in water

2D free-energy surfaces for transfer of the methoxymethyl cation between two water molecules are constructed from molecular dynamics (MD) simulations in which these atoms are treated quantum-mechanically within a box of 1030 classical solvent water molecules at 300 K. This provides a simple model for glycosyl transfer in water. The AM1/TIP3P surfaces with 2D-spline corrections at either MPWB1K/6-31+G(d,p) or MP2/6-31+G(d,p) contain a shallow free-energy well corresponding to an oxacarbenium ion intermediate in a DN*AN mechanism. MD analysis at three temperatures leads to a classical estimate of the lifetime of the methoxymethyl cation in water; when quantum corrections for vibrational zero-…

Química Física avanzada (2009/2010)

La asignatura Química Física Avanzada es una asignatura troncal de carácter anual. En el plan de estudios actualmente en vigor consta de un total de 9.0 créditos ECTS.