The ONIOM molecular dynamics method for biochemical applications: Cytidine deaminase

A new step in kinetic proofreading due to misacylated-tRNA during ribosomal peptide bond formation

The translational accuracy in protein synthesis is contributed to by several mechanisms in the ribosome, generally called kinetic proofreading. This process in the ribosome inhibits the non-cognate codon-anticodon interaction. However, it is not sufficient for fidelity of protein synthesis since a wrong amino acid can easily be added to the growing polypeptide chain if a tRNA while cognate to the mRNA, carries a non-cognate amino acid. Therefore, additional to the kinetic proofreading, there must be some hitherto unknown characteristic in misacylated-tRNAs to stop the process of protein synthesis if such misacylated-tRNA is accommodated in the ribosomal A-site. In order to understand this characteristic, we have performed computational quantum chemistry analysis on five different tRNA molecules, each one attached to five different amino acids with one being cognate to the tRNA and the other four non-cognate. This study shows the importance of aminoacyl-tRNA binding energy in ensuring fidelity of protein synthesis.

Covalent simulations of covalent/irreversible enzyme inhibition in drug discovery: a reliable technical protocol

Aim: Irreversible covalent inhibition of biological targets in disease pathogenesis is an emerging field in drug design. Computational techniques have assumed a critical role in understanding covalent enzyme inhibition. However, a gap currently exists with regards to the reliability and reproducibility of currently available protocols available in literature and open scientific forums. Methodology/results: Appropriate ligand and protein target are selected, docked covalently or noncovalently using respective docking tools. Both components are subjected to premolecular dynamic preparations. This was followed by parameterization of the ligand, protein and covalent complex, respectively. The production runs were initiated and the resulting trajectories are saved and analyzed. Conclusion: This protocol is reliable and reproducible, hence would advance the development of irreversible covalent inhibitors toward disease treatment.

Product release mechanism and the complete enzyme catalysis cycle in yeast cytosine deaminase (yCD): A computational study

Yeast cytosine deaminase (yCD) is critical in gene-directed enzyme prodrug therapy as it catalyzes the hydrolytic deamination of cytosine. The product (uracil) release process is considered as rate-limiting in the whole enzymatic catalysis and includes the cleavage of the uracil-metal bond and the delivery of free uracil out of the reactive site. Herein extensive combined random acceleration molecular dynamics (RAMD) and molecular dynamics (MD) simulations coupled with the umbrella sampling technique have been performed to study the product transport mechanism. Five channels have been identified, and the thermodynamic and dynamic characterizations for the two most favorable channels have been determined and analyzed. The free energy barrier for the most beneficial pathway is about 13kcal/mol and mainly results from the cleavage of hydrogen bonds between the ligand uracil and surrounding residues Asn51, Glu64, and Asp155. The conjugated rings of Phe114 and Trp152 play gating and guiding roles in the product delivery via π⋯π van der Waals interactions with the product. Finally, the full cycle of the enzymatic catalysis has been determined, making the whole process computationally more precise.

Quantum Mechanical and Molecular Dynamics Studies of the Reaction Mechanism of the Nucleophilic Substitution at the Si Atom

The mechanism of the nucleophilic substitution at the Si atom, SiH3Cl + Cl*(-) → SiH3Cl* + Cl(-), is examined by both quantum mechanical (QM) and molecular dynamics (MD) methods. This reaction proceeds by two steps with the inversion or retention of the configuration passing through an intermediate with the trigonal bipyramid (TBP) structure, although the conventional SN2 reaction at the C atom proceeds by one step with the inversion of the configuration passing through a transition state with the TBP structure. We followed by the QM calculations all the possible paths of the substitution reaction that undergo the TBP intermediates with the cis and trans forms produced by the frontside and backside attacks of Cl(-). As a result, it was thought that TBPcis1 produced with a high probability is readily transformed to the energetically more stable TBPtrans. This fact was also shown by the MD simulations. In order to obtain more information concerning the trajectory of Cl(-) on the dissociation from TBPtrans, which we cannot clarify on the basis of the energy profile determined by the QM method, the MD simulations with and without the water solvent were conducted and analyzed in detail. The QM-MD simulations without the water solvent revealed that the dissociation of Cl(-) from TBPtrans occurs without passing through TBPcis1'. The ONIOM-MD simulations with the water solvent further suggested that the thermal fluctuation of the water solvent significantly affects the oscillation of the kinetic and potential energies of the substrate to facilitate the isomerization of the TBP intermediate from the cis form to the trans form and the subsequent dissociation of Cl(-) from TBPtrans.

Enzyme catalysis by entropy without Circe effect

Entropic effects have often been invoked to explain the extraordinary catalytic power of enzymes. In particular, the hypothesis that enzymes can use part of the substrate-binding free energy to reduce the entropic penalty associated with the subsequent chemical transformation has been very influential. The enzymatic reaction of cytidine deaminase appears to be a distinct example. Here, substrate binding is associated with a significant entropy loss that closely matches the activation entropy penalty for the uncatalyzed reaction in water, whereas the activation entropy for the rate-limiting catalytic step in the enzyme is close to zero. Herein, we report extensive computer simulations of the cytidine deaminase reaction and its temperature dependence. The energetics of the catalytic reaction is first evaluated by density functional theory calculations. These results are then used to parametrize an empirical valence bond description of the reaction, which allows efficient sampling by molecular dynamics simulations and computation of Arrhenius plots. The thermodynamic activation parameters calculated by this approach are in excellent agreement with experimental data and indeed show an activation entropy close to zero for the rate-limiting transition state. However, the origin of this effect is a change of reaction mechanism compared the uncatalyzed reaction. The enzyme operates by hydroxide ion attack, which is intrinsically associated with a favorable activation entropy. Hence, this has little to do with utilization of binding free energy to pay the entropic penalty but rather reflects how a preorganized active site can stabilize a reaction path that is not operational in solution.

ONIOM and ab-initio calculations on the mechanism of uncatalyzed peptide bond formation

Finding a proper transition structure for the peptide bond formation process can lead one to a better understanding of the role of ribosome in catalyzing this reaction. Using computer simulations, we performed the potential energy surface scan on the ester bond dissociation of P-site aminoacyl-tRNA and the peptide bond formation of P-site and A-site amino acids. The full fragments of initiator tRNA(i)(met) and elongator tRNA(phe) are attached to both cognate and non-cognate amino acids as the P-site substrate. The A-site amino acid for all four calculations is methionine. We used ONIOM calculations to reduce the computational cost. Our study illustrates the reduced rate of peptide bond formation for misacylated tRNA(i)(met) in the absence of ribosomal bases. The misacylated elongator tRNA(phe), however, did not show any difference in its PES compared with that for the phe-tRNA(phe). This demonstrates the structural specification of initiator tRNA(i)(met) for the amino acids side chain.

The role of initiator tRNAimet in fidelity of initiation of protein synthesis

The proper arrangement of amino acids in a protein determines its proper function, which is vital for the cellular metabolism. This indicates that the process of peptide bond formation requires high fidelity. One of the most important processes for this fidelity is kinetic proofreading. As biochemical experiments suggest that kinetic proofreading plays a major role in ensuring the fidelity of protein synthesis, it is not certain whether or not a misacylated tRNA would be corrected by kinetic proofreading during the peptide bond formation. Using 2-layered ONIOM (QM/MM) computational calculations, we studied the behavior of misacylated tRNAs and compared the results with these for cognate aminoacyl-tRNAs during the process of peptide bond formation to investigate the effect of nonnative amino acids on tRNAs. The difference between the behavior of initiator tRNA(i) (met) compared to the one for the elongator tRNAs indicates that only the initiator tRNA(i) (met) specifies the amino acid side chain.

Systematic fragmentation method and the effective fragment potential: an efficient method for capturing molecular energies

The systematic fragmentation method fragments a large molecular system into smaller pieces, in such a way as to greatly reduce the computational cost while retaining nearly the accuracy of the parent ab initio electronic structure method. In order to attain the desired (sub-kcal/mol) accuracy, one must properly account for the nonbonded interactions between the separated fragments. Since, for a large molecular species, there can be a great many fragments and therefore a great many nonbonded interactions, computations of the nonbonded interactions can be very time-consuming. The present work explores the efficacy of employing the effective fragment potential (EFP) method to obtain the nonbonded interactions since the EFP method has been shown previously to capture nonbonded interactions with an accuracy that is often comparable to that of second-order perturbation theory. It is demonstrated that for nonbonded interactions that are not high on the repulsive wall (generally >2.7 A), the EFP method appears to be a viable approach for evaluating the nonbonded interactions. The efficacy of the EFP method for this purpose is illustrated by comparing the method to ab initio methods for small water clusters, the ZOVGAS molecule, retinal, and the alpha-helix. Using SFM with EFP for nonbonded interactions yields an error of 0.2 kcal/mol for the retinal cis-trans isomerization and a mean error of 1.0 kcal/mol for the isomerization energies of five small (120-170 atoms) alpha-helices.

Ab Initio ONIOM−Molecular Dynamics (MD) Study on the Deamination Reaction by Cytidine Deaminase

We applied the ONIOM-molecular dynamics (MD) method to the hydrolytic deamination of cytidine by cytidine deaminase, which is an essential step of the activation process of the anticancer drug inside the human body. The direct MD simulations were performed for the realistic model of cytidine deaminase by calculating the energy and its gradient by the ab initio ONIOM method on the fly. The ONIOM-MD calculations including the thermal motion show that the neighboring amino acid residue is an important factor of the environmental effects and significantly affects not only the geometry and energy of the substrate trapped in the pocket of the active site but also the elementary step of the catalytic reaction. We successfully simulate the second half of the catalytic cycle, which has been considered to involve the rate-determining step, and reveal that the rate-determining step is the release of the NH3 molecule.

An insight into the environmental effects of the pocket of the active site of the enzyme.Ab initio ONIOM-molecular dynamics (MD) study on cytosine deaminase

We applied the ONIOM-molecular dynamics (MD) method to cytosine deaminase to examine the environmental effects of the amino acid residues in the pocket of the active site on the substrate taking account of their thermal motion. The ab initio ONIOM-MD simulations show that the substrate uracil is strongly perturbed by the amino acid residue Ile33, which sandwiches the uracil with His62, through the steric contact due to the thermal motion. As a result, the magnitude of the thermal oscillation of the potential energy and structure of the substrate uracil significantly increases.