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-8719), but the reaction mechanisms are noticeably different. Our simulations reproduce the experimentally determined catalytic effect of ribosome that can be explained by the different behavior of the two environments. While the solvent reorganizes during the chemical process involving an entropic penalty, the ribosome is preorganized in the formation of the Michaelis complex and does not suffer important changes along the reaction, dampening the charge redistribution of the chemical system.

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 different conformers of this species, four different paths were found: three of them taking place through concerted mechanisms of four-, six- and eight-membered ring transition states, and a stepwise mechanism through a neutral intermediate. Comparison of the free energy profiles indicates that the concerted mechanisms would be kinetically favored, with free energy barriers in very good agreement with experimental data. Calculations of kinetic isotope effects, when including the solute interactions with the first solvation shell, show that the 8-membered ring TS renders values in better agreement with available experimental data. Quantitative discrepancies can be attributed to different employed models in experiments and calculations.

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.

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 substrate implies the existence of a huge number of possible conformers, some of them implying formation of intramolecular hydrogen bond interaction that can determine the energetics of the conformers defining the different states along the reaction profile. The results show PM3 provides results that are in qualitative agreement with DFT calculations. The free energy profiles show a step-wise mechanism that is kinetically determined by the nucleophilic attack of a water molecule activated by the proton transfer to the carboxylate group of the substrate (the first step). However, since the main role of the ß-lactamase would be reducing the free energy barrier of the first step, and keeping in mind the barrier obtained from second intermediate to products, population of this second intermediate could be significant and consequently experimentally detected in ß-lactamases, as shown in the literature.

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 salicylate production from isochorismate, that appears to be pericyclic in nature. Findings also explain the role of the magnesium cation and the pH dependence activity experimentally observed in MbtI. Mg(2+) would be polarizing and pre-organizing the substrate and active site, as well as shifting the pK(a) values of key active site residues.

Hybrid schemes based on quantum mechanics/molecular mechanics simulations goals to success, problems, and perspectives.

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 science that allows one to reach the molecular design of tailor-made enzyme catalysts. This is even more so when one refers to enzyme catalysts, for which the detailed structure and composition are known and can be manipulated to introduce well-defined residues which can be implicated in the chemical rearrangements taking place in the active site. The methods and techniques of theoretical and computational chemistry are becoming more and more important in both understanding the fundamental biological roles of enzymes and facilitating their utilization in biotechnology. Improvement of the catalytic function of enzymes is important from scientific and industrial viewpoints, and to put this fact in the actual perspective as well as the potentialities, we recommend the very recent report of Sanderson [Sanderson, K. (2011). Chemistry: enzyme expertise. Nature 471, 397.]. Great fundamental advances have been made toward the ab initio design of enzyme catalysts based on molecular modeling. This has been based on the molecular mechanistic knowledge of the reactions to be catalyzed, together with the development of advanced synthesis and characterization techniques. The corresponding molecular mechanism can be studied by means of powerful quantum chemical calculations. The catalytic active site can be optimized to improve the transition state analogues (TSA) and to enhance the catalytic activity, even improve the active site to favor a desired direction of some promiscuous enzymes. In this chapter, we give a brief introduction, the state of the art, and future prospects and implications of enzyme design. Current computational tools to assist experimentalists for the design and engineering of proteins with desired catalytic properties are described. The interplay between enzyme design, molecular simulations, and experiments will be presented to emphasize the interdisciplinary nature of this research field. This text highlights the recent advances and examples selected from our laboratory are shown, of how the applications of these tools are a first attempt to de novo design of protein active sites. Identification of neutral/advantageous/deleterious mutation platforms can be exploited to penetrate some of Nature's closely guarded secrets of chemical reactivity. In this chapter, we give a brief introduction, the state of the art, and future prospects and implications of enzyme design. The first part describes briefly how the molecular modeling is carried out. Then, we discuss the requirements of hybrid quantum mechanical/molecular mechanics molecular dynamics (QM/MM MD) simulations, analyzing what are the basis of these theoretical methodologies, how we can use them with a view to its application in the study of enzyme catalysis, and what are the best methodologies for assessing its catalytic potential. In the second part, we focus on some selected examples, taking as a common guide the chorismate to prephenate rearrangement, studying the corresponding molecular mechanism in vacuo, in solution and in an enzyme environment. In addition, examples involving catalytic antibodies (CAs) and promiscuous enzymes will be presented. Finally, a special emphasis is made to provide some hints about the logical evolution that can be anticipated in this research field. Moreover, it helps in understanding the open directions in this area of knowledge and highlights the importance of computational approaches in discovering specific drugs and the impact on the rational design of tailor-made enzymes.

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. According to our simulations the specialization in the hydrolysis of phosphomonoesters or phosphodiesters, developed in different members of the superfamily, is a consequence of the interactions established between the protein and the oxygen atoms of the phosphate group and, in particular, with the oxygen atom that bears the additional alkyl group when the substrate is a diester. A water molecule, belonging to the coordination shell of the Mg(2+) ion, and residue Lys328 seem to play decisive roles stabilizing a phosphomonoester substrate, but the latter contributes to increase the energy barrier for the hydrolysis of phosphodiesters. Then, mutations affecting the nature or positioning of Lys328 lead to an increased diesterase activity in AP. Finally, the capacity of this enzymatic family to catalyze the reaction of phosphoesters having different leaving groups, or substrate promiscuity, is explained by the ability of the enzyme to stabilize different charge distributions in the leaving group using different interactions involving either one of the zinc centers or residues placed on the outer side of the catalytic site.

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 dissociative character. The reason for this change is rationalized in terms of the very different nature of the electrostatic interactions established in each of the environments: while the aqueous solution reduces the repulsion between the negatively charged reacting fragments, assisting their approach, the NPP active site stabilizes the charge distribution of dissociative TS structures, allowing the reaction to proceed with a significantly reduced free energy cost. Interestingly, the NPP active site is able to accommodate different substrates, and it seems that the nature of the TSs depends on their electronic characteristics. So, in the case of the MpNPP(-) substrate, the nitro group establishes hydrogen-bond interactions with water molecules and residues found in the outer part of the catalytic site, while the leaving group oxygen atom does not coordinate directly with any of the zinc atoms of the active site. If methyl phenyl phosphate is used as substrate, then the charge on the leaving group is supported to larger extent by the oxygen atom and the phenolate anion can be then coordinated to one of the two zinc atoms present in the active site.

Mechanism and plasticity of isochorismate pyruvate lyase: a computational study.

The isochorismate pyruvate lyase (IPL) from Pseudomonas aeruginosa, designated as PchB, catalyzes the transformation of isochorismate into pyruvate and salicylate, but it also catalyzes the rearrangement of chorismate into prephenate, suggesting that both reactions may proceed by a pericyclic mechanism. In this work, molecular dynamics simulations employing hybrid quantum mechanics/molecular mechanics methods have been carried out to get a detailed knowledge of the reaction mechanism of PchB. The results provide a theoretical rate constant enhancement by comparison with the reaction in solution, in agreement with the experimental data, and confirm the pericyclic nature of the reaction mechanism. The robustness of this promiscuous enzyme has been checked by considering the impact of Ala37Ile mutation, previously proposed by us to improve the secondary chorismate mutase (CM) activity. The effect of this mutation, which was shown to increase the rate constant for the CM activity by a factor of 10(3), also increases the IPL catalytic efficiency, although only by a factor of 6.

Theoretical modeling of the reaction mechanism of phosphate monoester hydrolysis in alkaline phosphatase.

The reaction mechanism of phosphate monoester hydrolysis in alkaline phosphatase is analyzed by means of hybrid QM/MM simulations. A recently developed semiempirical Hamiltonian, AM1/d-PhoT, which takes into account the d orbitals on the phosphorus atom, has been employed. The reaction mechanism obtained is either associative or dissociative, depending on the size of the QM subsystem. The results are rationalized on the basis of the degree of charge transfer from the reacting fragments to the two zinc ions present in the active site, which has been observed to be dependent on whether or not metal atoms and their coordination spheres are included in the QM region. The description obtained using the largest QM region agrees with the picture obtained from experimental data.

Influence of ionization on the conformational preferences of peptide models. Ramachandran surfaces of N-formyl-glycine amide and N-formyl-alanine amide radical cations.

Ramachandran maps of neutral and ionized HCO-Gly-NH2 and HCO-Ala-NH2 peptide models have been built at the B3LYP/6-31++G(d,p) level of calculation. Direct optimizations using B3LYP and the recently developed MPWB1K functional have also been carried out, as well as single-point calculations at the CCSD(T) level of theory with the 6-311++G(2df,2p) basis set. Results indicate that for both peptide models ionization can cause drastic changes in the shape of the PES in such a way that highly disallowed regions in neutral PES become low-energy regions in the radical cation surface. The structures localized in such regions, epsilonL+* and epsilonD+* are highly stabilized due to the formation of 2-centre-3-electron interactions between the two carbonyl oxygens. Inclusion of solvent effects by the conductor-like polarizable continuum model (CPCM) shows that the solute-solvent interaction energy plays an important role in determining the stability order.

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.

Cu2+/+ cation coordination to adenine--thymine base pair. Effects on intermolecular proton-transfer processes.

Intermolecular proton-transfer processes in the Watson & Crick adenine-thymine Cu+ and Cu2+ cationized base pairs have been studied using the density functional theory (DFT) methods. Cationized systems subject to study are those resulting from cation coordination to the main basic sites of the base pair, N7 and N3 of adenine and O2 of thymine. For Cu+ coordinated to N7 or N3 of adenine, only the double proton-transferred product is found to be stable, similarly to the neutral system. However, when Cu+ interacts with thymine, through the O2 carbonyl atom, the single proton transfer from thymine to adenine becomes thermodynamically spontaneous, and thus rare forms of the DNA bases may spontaneously appear. For Cu2+ cation, important effects on proton-transfer processes appear due to oxidation of the base pair, which stabilizes the different single proton-transfer products. Results for hydrated systems show that the presence of the water molecules interacting with the metal cation (and their mode of coordination) can strongly influence the ability of Cu2+ to induce oxidation on the base pair.

Attractive strain: the disadvantages of rigid multiple H-bond donors and acceptors. A theoretical analysis of the hydrogen-bonding interactions in complexes of tetraazaanthracenedione with pyridylureas.

We report density functional theory calculations at the B3LYP/D95(d,p) level on the hydrogen-bonding complexes of tetraazaanthracenedione, 1, with N-(pyridin-2-yl)urea, 2H, or N-(6-aminopyridin-2-yl)urea, 2N. The interaction energy of the 1-2H complex exceeds that of 1-2N, despite the fact that 1-2N contains a strong N-H...O interaction in place of a weak C-H...O interaction in 1-2H. We show that the 1-2N interaction is weaker than the sum of the four normal individual H-bonding interactions because the steric constraints of the complex prevent the H-bonding donors and acceptors from optimally approaching each other to form the two central H-bonds. This steric phenomenon, which we call attractive strain, is likely present to some extent in most H-bonding systems that contain more than two H-bonds between rigid monomers. Attractive strain is unusually important in 1-2N. Attractive strain can be conceived of as an enthalpic cost for the entropic benefits of freezing the dihedral angles of the multiple H-bond donors and acceptors by designing rigid systems.

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, together with internal energy minimizations, has allowed free-energy profiles to be obtained for both CAs. The profiles show a smaller barrier for the matured form, which is in agreement with the experimental observation. Free-energy profiles were obtained with this methodology, thereby avoiding the much more demanding two-dimensional calculations of the energy surfaces that are normally required to study this kind of reaction. Structural analysis and energy evaluations of substrate-protein interactions have been performed from averaged structures, which allows understanding of how the single mutations carried out during the maturation process can be responsible for the observed fourfold enhancement of the catalytic rate constant. The conclusion is that the mutation effect in this studied germline CA produces a complex indirect effect through coupled movements of the backbone of the protein and the substrate.

Influence of the Side Chain in the Structure and Fragmentation of Amino Acids Radical Cations.

The conformational properties of ionized amino acids (Gly, Ala, Ser, Cys, Asp, Gln, Phe, Tyr, and His) have been theoretically analyzed using the hybrid B3LYP and the hybrid-meta MPWB1K functionals as well as with the post-Hartree Fock CCSD(T) level of theory. As a general trend, ionization is mainly localized at the -NH2 group, which becomes more planar and acidic, the intramolecular hydrogen bond in which -NH2 acts as proton donor being strengthened upon ionization. For this reason, the so-called conformer IV(+) becomes the most stable for nonaromatic amino acid radical cations. Aromatic amino acids do not follow this trend because ionization takes place mainly at the side chain. For these amino acids for which ionization of the side chain prevails over the -NH2 group, structures III(+) and II(+) become competitive. The Cα-X fragmentations of the ionized systems have also been studied. Among the different decompositions considered, the one that leads to the loss of COOH(•) is the most favorable one. Nevertheless, for aromatic amino acids fragmentations leading to R(•) or R(+) start being competitive. In fact, for His and Tyr, results indicate that the fragmentation leading to R(+) is the most favorable process.

Stereoselectivity behavior of the AZ28 antibody catalyzed oxy-Cope rearrangement.

Catalytic antibodies are very interesting not only because of the rate enhancement of the reactions that they catalyze but also because of the selectivities they can achieve that are sometimes not present in natural enzyme processes. We have selected the study of the stereoselectivity of the matured AZ28 that catalyzes an oxy-Cope rearrangement. For this particular case, the presence of a chiral center in the substrate provokes the existence of two different enantiomers, R and S. Furthermore, it is also possible to locate two different orientations for the hydroxyl group in the central ring of the substrate in the transition state, equatorial and axial, rendering two different conformers. In this paper we present the free energy profiles obtained for different substrate isomers in the cavity created by the matured catalytic antibody. Our simulations have reproduced the stereoselectivity of the matured AZ28, differentiating between the axial or equatorial orientations and preferentially stabilizing the S forms, at a qualitative level. Finally, the inclusion of the substrate-CA interactions in a flexible molecular model has allowed us to observe the different pattern of interactions that are related to different interaction energies, which seem to be crucial in the stereoselectivity behavior of the catalytic antibody.

On the nature of the transition state in catechol O-methyltransferase. A complementary study based on molecular dynamics and potential energy surface explorations.

The way in which enzymes influence the rate of chemical processes is still a question of debate. The protein promotes the catalysis of biochemical processes by lowering the free energy barrier in comparison with the reference uncatalyzed reaction in solution. In this article we are reporting static and dynamic aspects of the enzyme catalysis in a bimolecular reaction, namely a methyl transfer from S-adenosylmethionine to the hydroxylate oxygen of a substituted catechol catalyzed by catechol O-methyltransferase. From QM/MM optimizations, we will first analyze the participation of the environment on the transition vector. The study of molecular dynamics trajectories will allow us to estimate the transmission coefficient from a previously localized transition state as the maximum in the potential of mean force profile. The analysis of the reactive and nonreactive trajectories in the enzyme environment and in solution will also allow studying the geometrical and electronic changes, with special attention to the chemical system movements and the coupling with the environment. The main result, coming from both analyses, is the approximation of the magnesium cation to the nucleophilic and the hydroxyl group of the catecholate as a result of a general movement of the protein, stabilizing in this way the transition state. Consequently, the free energy barrier of the enzyme reaction is dramatically decreased with respect to the reaction in solution.