Quantum mechanics/molecular mechanics investigation of the chemical reaction in Dpo4 reveals water-dependent pathways and requirements for active site reorganization

The nucleotidyl-transfer reaction coupled with the conformational transitions in DNA polymerases is critical for maintaining the fidelity and efficiency of DNA synthesis. We examine here the possible reaction pathways of a Y-family DNA polymerase, Sulfolobus solfataricus DNA polymerase IV (Dpo4), for the correct insertion of dCTP opposite 8-oxoguanine using the quantum mechanics/molecular mechanics (QM/MM) approach, both from a chemistry-competent state and a crystal closed state. The latter examination is important for understanding pre-chemistry barriers to interpret the entire enzyme mechanism, since the crystal closed state is not an ideal state for initiating the chemical reaction. The most favorable reaction path involves initial deprotonation of O3'H via two bridging water molecules to O1A, overcoming an overall potential energy barrier of approximately 20.0 kcal/mol. The proton on O1A-P(alpha) then migrates to the gamma-phosphate oxygen of the incoming nucleotide as O3' attacks P(alpha), and the P(alpha)-O3A bond breaks. The other possible pathway in which the O3'H proton is transferred directly to O1A on P(alpha) has an overall energy barrier of 25.0 kcal/mol. In both reaction paths, the rate-limiting step is the initial deprotonation, and the trigonal-bipyramidal configuration for P(alpha) occurs during the concerted bond formation (O3'-P(alpha)) and breaking (P(alpha)-O3A), indicating the associative nature of the chemical reaction. In contrast, the Dpo4/DNA complex with an imperfect active-site geometry corresponding to the crystal state must overcome a much higher activation energy barrier (29.0 kcal/mol) to achieve a tightly organized site due to hindered O3'H deprotonation stemming from larger distances and distorted conformation of the proton acceptors. This significant difference demonstrates that the pre-chemistry reorganization in Dpo4 costs approximately 4.0 to 9.0 kcal/mol depending on the primer terminus environment. Compared to the higher fidelity DNA polymerase beta from the X-family, Dpo4 has a higher chemical reaction barrier (20.0 vs 15.0 kcal/mol) due to the more solvent-exposed active site.     

A quantum mechanical investigation of possible mechanisms for the nucleotidyl transfer reaction catalyzed by DNA polymerase beta

Several quantum mechanical (QM) and hybrid quantum/molecular mechanical (QM/MM) studies have been employed recently to analyze the nucleotidyl transfer reaction in DNA polymerase beta (pol beta). Our examination reveals strong dependence of the reported mechanism on the initial molecular model. Thus, we explore here several model systems by QM methods to investigate pol beta's possible pathway variations. Although our most favorable pathway involves a direct proton transfer from O3'(primer) to O2alpha(Palpha), we also discuss other initial proton-transfer steps--to an adjacent water, to triphosphate, or to aspartic units--and the stabilizing effect of crystallographic water molecules in the active site. Our favored reaction route has an energetically undemanding initial step of less than 1.0 kcal/mol (at the B3LYP/6-31G(d,p) level), and involves a slight rearrangement in the geometry of the active site. This is followed by two major steps: (1) direct proton transfer from O3'(primer) to O2alpha(Palpha) leading to the formation of a pentavalent, trigonal bipyramidal Palpha center, via an associative mechanism, at a cost of about 28 kcal/mol, and (2) breakage of the triphosphate unit (exothermic process, approximately 22 kcal/mol) that results in the full transfer of the nucleotide to the DNA and the formation of pyrophosphate. These energy values are expected to be lower in the physical system when full protein effects are incorporated. We also discuss variations from this dominant pathway, and their impact on the overall repair process. Our calculated barrier for the chemical reaction clearly indicates that chemistry is rate-limiting overall for correct nucleotide insertion in pol beta, in accord with other studies. Protonation studies on relevant intermediates suggest that, although protonation at a single aspartic residue may occur, the addition of a second proton to the system significantly disturbs the active site. We conclude that the active site rearrangement step necessary to attain a reaction-competent geometry is essential and closely related to the "pre-chemistry" avenue described recently as a key step in the overall kinetic cycle of DNA polymerases. Thus, our work emphasizes the many possible ways for DNA polymerase beta's chemical reaction to occur, determined by the active site environment and initial models.     

Computer simulation of the chemical catalysis of DNA polymerases: discriminating between alternative nucleotide insertion mechanisms for T7 DNA polymerase

Understanding the chemical step in the catalytic reaction of DNA polymerases is essential for elucidating the molecular basis of the fidelity of DNA replication. The present work evaluates the free energy surface for the nucleotide transfer reaction of T7 polymerase by free energy perturbation/empirical valence bond (FEP/EVB) calculations. A key aspect of the enzyme simulation is a comparison of enzymatic free energy profiles with the corresponding reference reactions in water using the same computational methodology, thereby enabling a quantitative estimate for the free energy of the nucleotide insertion reaction. The reaction is driven by the FEP/EVB methodology between valence bond structures representing the reactant, pentacovalent intermediate, and the product states. This pathway corresponds to three microscopic chemical steps, deprotonation of the attacking group, a nucleophilic attack on the P(alpha) atom of the dNTP substrate, and departure of the leaving group. Three different mechanisms for the first microscopic step, the generation of the RO(-) nucleophile from the 3'-OH hydroxyl of the primer, are examined: (i) proton transfer to the bulk solvent, (ii) proton transfer to one of the ionic oxygens of the P(alpha) phosphate group, and (iii) proton transfer to the ionized Asp654 residue. The most favorable reaction mechanism in T7 pol is predicted to involve the proton transfer to Asp654. This finding sheds light on the long standing issue of the actual role of conserved aspartates. The structural preorganization that helps to catalyze the reaction is also considered and analyzed. The overall calculated mechanism consists of three subsequent steps with a similar activation free energy of about 12 kcal/mol. The similarity of the activation barriers of the three microscopic chemical steps indicates that the T7 polymerase may select against the incorrect dNTP substrate by raising any of these barriers. The relative height of these barriers comparing right and wrong dNTP substrates should therefore be a primary focus of future computational studies of the fidelity of DNA polymerases.     

DNA polymerase beta catalysis: are different mechanisms possible?

DNA polymerases are crucial constituents of the complex cellular machinery for replicating and repairing DNA. Discerning mechanistic pathways of DNA polymerase on the atomic level is important for revealing the origin of fidelity discrimination. Mammalian DNA polymerase beta (pol beta), a small (39 kDa) member of the X-family, represents an excellent model system to investigate polymerase mechanisms. Here, we explore several feasible low-energy pathways of the nucleotide transfer reaction of pol beta for correct (according to Watson-Crick hydrogen bonding) G:C basepairing versus the incorrect G:G case within a consistent theoretical framework. We use mixed quantum mechanics/molecular mechanics (QM/MM) techniques in a constrained energy minimization protocol to effectively model not only the reactive core but also the influence of the rest of the enzymatic environment and explicit solvent on the reaction. The postulated pathways involve initial proton abstraction from the terminal DNA primer O3'H group, nucleophilic attack that extends the DNA primer chain, and elimination of pyrophosphate. In particular, we analyze several possible routes for the initial deprotonation step: (i) direct transfer to a phosphate oxygen O(Palpha) of the incoming nucleotide, (ii) direct transfer to an active site Asp group, and (iii) transfer to explicit water molecules. We find that the most probable initial step corresponds to step (iii), involving initial deprotonation to water, which is followed by proton migration to active site Asp residues, and finally to the leaving pyrophosphate group, with an activation energy of about 15 kcal/mol. We argue that initial deprotonation steps (i) and (ii) are less likely as they are at least 7 and 11 kcal/mol, respectively, higher in energy. Overall, the rate-determining step for both the correct and the incorrect nucleotide cases is the initial deprotonation in concert with nucleophilic attack at the phosphate center; however, the activation energy we obtain for the mismatched G:G case is 5 kcal/mol higher than that of the matched G:C complex, due to active site structural distortions. Taken together, our results support other reported mechanisms and help define a framework for interpreting nucleotide specificity differences across polymerase families, in terms of the concept of active site preorganization or the so-called "pre-chemistry avenue".     

Biochemical Characterization of Translesion Synthesis by Sulfolobus acidocaldarius DNA Polymerases

To study the DNA synthesis mechanism of Sulfolobus acidocaldarius, a thermophilic species from Crenarchaeota, two DNA polymerases of B family(polB1 and polB3), and one DNA polymerase of Y family(polIV) were recombinantly expressed, purified and biochemically characterized. Both DNA polymerases polB1(Saci_1537) and polB3(Saci_0074) possessed DNA polymerase and 3' to 5' exonuclease activities; however, both the activities of B3 were very inefficient in vitro. The polIV(Saci_0554) was a polymerase, not an exonuclease. The activities of all the three DNA polymerases were dependent on divalent metal ions Mn²⁺ and Mg²⁺. They showed the highest activity at pH values ranging from 8.0 to 9.5. Their activities were inhibited by KCl with high concentration. The optimal reaction temperatures for the three DNA polymerases were between 60 and 70℃. Deaminated bases dU and dI on DNA template strongly hindered primer extension by the two DNA polymerases of B family, not by the DNA polymerase of Y family. DNA polymerase of Y Family bypassed the two AP site analogues dSpacer and propane on template more easily than DNA polymerases of B family. Our results suggest that the three DNA polymerases coordinate to fulfill various DNA synthesis in Sulfolobus acidocaldarius cell.     

Fundamental reaction pathway and free energy profile for inhibition of proteasome by Epoxomicin

First-principles quantum mechanical/molecular mechanical free energy calculations have been performed to provide the first detailed computational study on the possible mechanisms for reaction of proteasome with a representative peptide inhibitor, Epoxomicin (EPX). The calculated results reveal that the most favorable reaction pathway consists of five steps. The first is a proton transfer process, activating Thr1-O(γ) directly by Thr1-N(z) to form a zwitterionic intermediate. The next step is nucleophilic attack on the carbonyl carbon of EPX by the negatively charged Thr1-O(γ) atom, followed by a proton transfer from Thr1-N(z) to the carbonyl oxygen of EPX (third step). Then, Thr1-N(z) attacks on the carbon of the epoxide group of EPX, accompanied by the epoxide ring-opening (S(N)2 nucleophilic substitution) such that a zwitterionic morpholino ring is formed between residue Thr1 and EPX. Finally, the product of morpholino ring is generated via another proton transfer. Noteworthy, Thr1-O(γ) can be activated directly by Thr1-N(z) to form the zwitterionic intermediate (with a free energy barrier of only 9.9 kcal/mol), and water cannot assist the rate-determining step, which is remarkably different from the previous perception that a water molecule should mediate the activation process. The fourth reaction step has the highest free energy barrier (23.6 kcal/mol) which is reasonably close to the activation free energy (～21-22 kcal/mol) derived from experimental kinetic data. The obtained novel mechanistic insights should be valuable for not only future rational design of more efficient proteasome inhibitors but also understanding the general reaction mechanism of proteasome with a peptide or protein.     

Combined QM(DFT)/MM molecular dynamics simulations of the deamination of cytosine by yeast cytosine deaminase (yCD)

Extensive combined quantum mechanical (B3LYP/6‐31G*) and molecular mechanical (QM/MM) molecular dynamics simulations have been performed to elucidate the hydrolytic deamination mechanism of cytosine to uracil catalyzed by the yeast cytosine deaminase (yCD). Though cytosine has no direct binding to the zinc center, it reacts with the water molecule coordinated to zinc, and the adjacent conserved Glu64 serves as a general acid/base to shuttle protons from water to cytosine. The overall reaction consists of several proton‐transfer processes and nucleophilic attacks. A tetrahedral intermediate adduct of cytosine and water binding to zinc is identified and similar to the crystal structure of yCD with the inhibitor 2‐pyrimidinone. The rate‐determining step with the barrier of 18.0 kcal/mol in the whole catalytic cycle occurs in the process of uracil departure where the proton transfer from water to Glu64 and nucleophilic attack of the resulting hydroxide anion to C2 of the uracil ring occurs synchronously. © 2016 Wiley Periodicals, Inc.     

A theoretical study of the mechanism of the nucleotidyl transfer reaction catalyzed by yeast RNA polymerase II

The mechanism of the nucleotidyl transfer reaction catalyzed by yeast RNA polymerase II has been investigated using molecular mechanics and quantum mechanics methods.Molecular dynamics(MD) simulations were carried out using the TIP3 water model and generalized solvent boundary potential(GSBP) by CHARMM based on the X-ray crystal structure.Two models of the ternary elongation complex were constructed based on CHARMM MD calculations.All the species including reactants,transition states,intermediates,and products were optimized using the DFT-PBE method coupled with the basis set DZVP and the auxiliary basis set GEN-A2.Three pathways were explored using the DFT method.The most favorable reaction pathway involves indirect proton migration from the RNA primer 3’-OH to the oxygen atom of-phosphate via a solvent water molecule,proton rotation from the oxygen atom of-phosphate to the-phosphate side,the RNA primer 3’-O nucleophilic attack on the-phosphorus atom,and P-O bond breakage.The corresponding reaction potential profile was obtained.The rate limiting step,with a barrier height of 21.5 kcal/mol,is the RNA primer 3’-O nucleophilic attack,rather than the commonly considered proton transfer process.A high-resolution crystal structure including crystallographic water molecules is required for further studies.     

Ab initio QM/MM study shows there is no general acid in the reaction catalyzed by 4-oxalocrotonate tautomerase

The mechanism for the reaction catalyzed by the 4-oxalocrotonate tautomerase (4-OT) enzyme has been studied using a quantum mechanical/molecular mechanical (QM/MM) method developed in our laboratory. Total free energy barriers were obtained for the two steps involved in this reaction. In the first step, Pro-1 acts as a general base to abstract a proton from the third carbon of the substrate, 2-oxo-4-hexenedioate, creating a negative charge on the oxygen at C-2 of this substrate. In the second step, the same hydrogen abstracted by the N-terminal Pro-1 is shuttled back to the fifth carbon of the substrate to form the product, 2-oxo-3-hexenedioate. The calculated total free energy barriers are 14.54 and 16.45 kcal/mol for the first and second steps, respectively. Our calculations clearly show that there is no general acid in the reaction. Arg-39' ', which is hydrogen bonded to the carboxylate group of the substrate, and an ordered water, which moves closer to the site of the charge formed in the transition state and intermediate, play the main role in transition state/intermediate stabilization without acting as general acids in the reaction.     

Hydration effects on the reaction with an open-shell transition state: QM/MM-ER study for the dehydration reaction of alcohol in hot water

The reaction mechanism for the dehydration of 1,4-butanediol in hot water has been investigated by means of the hybrid quantum mechanical/molecular mechanical approach combined with the theory of energy representation (QM/MM-ER). We have assumed that the proton transfers along the hydrogen bonds of the water molecules catalyze the reaction, where the transition state (TS) forms a singlet biradical electronic structure. It has been revealed by the simulation that the biradical electronic state at the TS changes to zwitterionic structure in solution due to the hydration of the polar solvent. Such the electronic structure change gives rise to the substantial stabilization of the TS in hot water. As a result, the water-catalytic path becomes more favorable in aqueous solution than another possible path that proceeds without proton transfers as opposed to the reaction mechanism in the gas phase. Furthermore, the activation free energy computed by the present method is in excellent agreement with the experimental result.     

Encoding phenotype in bacteria with an alternative genetic set

An unnatural base-pair architecture with base pairs 2.4 ? larger than the natural DNA-based genetic system (xDNA) is evaluated for its ability to function like DNA, encoding amino acids in the context of living cells. xDNA bases are structurally analogous to natural bases but homologated by the width of a benzene ring, increasing their sizes and resulting in a duplex that is wider than native B-DNA. Plasmids encoding green fluorescent protein were constructed to contain single and multiple xDNA bases (as many as eight) in both strands and were transformed into Escherichia coli. Although they yielded fewer colonies than the natural control plasmid, in all cases in which a modified plasmid (containing one, two, three, or four consecutive size-expanded base pairs) was used, the correct codon bases were substituted, yielding green colonies. All four xDNA bases (xA, xC, xG, and xT) were found to encode the correct partners in the replicated plasmid DNA, both alone and in longer segments of xDNA. Controls with mutant cell lines having repair functions deleted were found to express the gene correctly, ruling out repair of xDNA and confirming polymerase reading of the unnatural bases. Preliminary experiments with polymerase deletion mutants suggested combined roles of replicative and lesion-bypass polymerases in inserting correct bases opposite xDNA bases and in bypassing the xDNA segments. These experiments demonstrate a biologically functioning synthetic genetic set with larger-than-natural architecture.     

Total synthesis of dehydroaltenusin

The first total synthesis of dehydroaltenusin, a natural enzyme inhibitor, is described. The key step involves Suzuki-coupling reaction of an aryl triflate prepared from 2,4,6-trihydroxybenzoic acid with a catechol-derived boronic acid or boronic ester. The synthetic product was evaluated as a potent inhibitor against eukaryotic DNA polymerase α and other DNA polymerases.     

Computational analysis of the proton translocation from Asp96 to schiff base in bacteriorhodopsin

The potential energy change during the M --> N process in bacteriorhodopsin has been evaluated by ab initio quantum chemical and advanced quantum chemical calculations following molecular dynamics (MD) simulations. Many previous experimental studies have suggested that the proton transfer from Asp96 to the Schiff base occurs under the following two conditions: (1) the hydrogen bond between Thr46 and Asp96 breaks and Thr46 is detached from Asp96 and (2) a stable chain of four water molecules spans an area from Asp96 --> Schiff base. In this work, we successfully reproduced the proton-transfer process occurring under these two conditions by molecular dynamics and quantum chemical calculations. The quantum chemical computation revealed that the proton transfer from Asp96 to Shiff base occurs in two-step reactions via an intermediate in which an H(3)O(+) appears around Ala215. The activation energy for the proton transfer in the first reaction was calculated to be 9.7 kcal/mol, which enables fast and efficient proton pump action. Further QM/MM (quantum mechanical/molecular mechanical) and FMO (fragment molecular orbital) calculations revealed that the potential energy change during the proton transfer is tightly regulated by the composition and the geometry of the surrounding amino acid residues of bacteriorhodopsin. Here, we report in detail the Asp96 --> Schiff base proton translocation mechanism of bacteriorhodopsin. Additionally, we discuss the effectiveness of combining quantum chemical calculations with truncated cluster models followed by advanced quantum chemical calculations applied to a whole protein to elucidate its reaction mechanism.     

DNA polymerase template interactions probed by degenerate isosteric nucleobase analogs

The development of novel artificial nucleobases and detailed X-ray crystal structures for primer/template/DNA polymerase complexes provide opportunities to assess DNA-protein interactions that dictate specificity. Recent results have shown that base pair shape recognition in the context of DNA polymerase must be considered a significant component. The isosteric azole carboxamide nucleobases (compounds 1-5; ) differ only in the number and placement of nitrogen atoms within a common shape and therefore present unique electronic distributions that are shown to dictate the selectivity of template-directed nucleotide incorporation by DNA polymerases. The results demonstrate how nucleoside triphosphate substrate selection by DNA polymerase is a complex phenomenon involving electrostatic interactions in addition to hydrogen bonding and shape recognition. These azole nucleobase analogs offer unique molecular tools for probing nonbonded interactions dictating substrate selection and fidelity of DNA polymerases.     

Arrest of Rolling Circle Amplification by Protein‐Binding DNA Aptamers

Certain DNA polymerases, such as ?29 DNA polymerase, can isothermally copy the sequence of a circular template round by round in a process known as rolling circle amplification (RCA), which results in super‐long single‐stranded (ss) DNA molecules made of tandem repeats. The power of RCA reflects the high processivity and the strand‐displacement ability of these polymerases. In this work, the ability of ?29DNAP to carry out RCA over circular templates containing a protein‐binding DNA aptamer sequence was investigated. It was found that protein–aptamer interactions can prevent this DNA polymerase from reading through the aptameric domain. This finding indicates that protein‐binding DNA aptamers can form highly stable complexes with their targets in solution. This novel observation was exploited by translating RCA arrest into a simple and convenient colorimetric assay for the detection of specific protein targets, which continues to showcase the versatility of aptamers as molecular recognition elements for biosensing applications.     

Nucleotide-mimetic synthetic ligands for DNA-recognizing enzymes One-step purification of Pfu DNA polymerase

The commercial availability of DNA polymerases has revolutionized molecular biotechnology and certain sectors of the bio-industry. Therefore, the development of affinity adsorbents for purification of DNA polymerases is of academic interest and practical importance. In the present study we describe the design, synthesis and evaluation of a combinatorial library of novel affinity ligands for the purification of DNA polymerases (Pols). Pyrococcus furiosus DNA polymerase (Pfu Pol) was employed as a proof-of-principle example. Affinity ligand design was based on mimicking the natural interactions between deoxynucleoside-triphosphates (dNTPs) and the B-motif, a conserved structural moiety found in Pol-I and Pol-II family of enzymes. Solid-phase 'structure-guided' combinatorial chemistry was used to construct a library of 26 variants of the B-motif-binding 'lead' ligand X-Trz-Y (X is a purine derivative and Y is an aliphatic/aromatic sulphonate or phosphonate derivative) using 1,3,5-triazine (Trz) as the scaffold for assembly. The 'lead' ligand showed complementarity against a Lys and a Tyr residue of the polymerase B-motif. The ligand library was screened for its ability to bind and purify Pfu Pol from Escherichia coli extract. One immobilized ligand (oABSAd), bearing 9-aminoethyladenine (AEAd) and sulfanilic acid (oABS) linked on the triazine scaffold, displayed the highest purifying ability and binding capacity (0,55 mg Pfu Pol/g wet gel). Adsorption equilibrium studies with this affinity ligand and Pfu Pol determined a dissociation constant (K(D)) of 83 nM for the respective complex. The oABSAd affinity adsorbent was exploited in the development of a facile Pfu Pol purification protocol, affording homogeneous enzyme (>99% purity) in a single chromatography step. Quality control tests showed that Pfu Pol purified on the B-motif-complementing ligand is free of nucleic acids and contaminating nuclease activities, therefore, suitable for experimental use.     

"Proton holes" in long-range proton transfer reactions in solution and enzymes: A theoretical analysis

Proton transfers are fundamental to chemical processes in solution and biological systems. Often, the well-known Grotthuss mechanism is assumed where a series of sequential "proton hops" initiates from the donor and combines to produce the net transfer of a positive charge over a long distance. Although direct experimental evidence for the sequential proton hopping has been obtained recently, alternative mechanisms may be possible in complex molecular systems. To understand these events, all accessible protonation states of the mediating groups should be considered. This is exemplified by transfers through water where the individual water molecules can exist in three protonation states (water, hydronium, and hydroxide); as a result, an alternative to the Grotthuss mechanism for a proton transfer through water is to generate a hydroxide by first protonating the acceptor and then transfer the hydroxide toward the donor through water. The latter mechanism can be most generally described as the transfer of a "proton hole" from the acceptor to the donor where the "hole" characterizes the deprotonated state of any mediating molecule. This pathway is distinct and is rarely considered in the discussion of proton-transfer processes. Using a calibrated quantum mechanical/molecular mechanical (QM/MM) model and an effective sampling technique, we study proton transfers in two solution systems and in Carbonic Anhydrase II. Although the relative weight of the "proton hole" and Grotthuss mechanisms in a specific system is difficult to determine precisely using any computational approach, the current study establishes an energetics motivated framework that hinges on the donor/acceptor pKa values and electrostatics due to the environment to argue that the "proton hole" transfer is likely as important as the classical Grotthuss mechanism for proton transport in many complex molecular systems.     

Theoretical thermodynamic study of the adenine–thymine tautomeric equilibrium: Electronic structure calculations and steered molecular dynamic simulations

Free energy of the tautomeric equilibrium A‐T ? A*‐T* between the canonical and noncanonical DNA base equilibrium in aqueous solution was theoretically determined by applying electronic structure methods (at the M06‐2X‐PCM/6‐311++G(d,p) level) and steered molecular dynamic simulations. Concerted and stepwise mechanisms were considered for the double proton transfer in an effort to explain the anomalous behavior of this system where an unfavorable process without a transition state can be observed depending on the level of calculation used. Of the different mechanisms used in the simulations, the stepwise mechanism, in which the first step implies the transference of a proton from thymine to adenine, and a second step with the transference of a different proton from adenine to thymine, was the only one that showed two transition states and a reaction intermediate. However, a concerted and stepwise mechanism has similar kinetic and thermodynamic behavior, with similar reaction and activation energies. Simple proton transfer was more favorable for the transference of the hydrogen from the adenine to the thymine. The inclusion of an aqueous medium in this study only slightly modified these energies, but the barrier energy was higher when the solvent was described as a discrete medium. Transition states and intermediate structures were analyzed at molecular dynamic level.     

Catalytic mechanism of RNA backbone cleavage by ribonuclease H from quantum mechanics/molecular mechanics simulations

We use quantum mechanics/molecular mechanics simulations to study the cleavage of the ribonucleic acid (RNA) backbone catalyzed by ribonuclease H. This protein is a prototypical member of a large family of enzymes that use two-metal catalysis to process nucleic acids. By combining Hamiltonian replica exchange with a finite-temperature string method, we calculate the free energy surface underlying the RNA-cleavage reaction and characterize its mechanism. We find that the reaction proceeds in two steps. In a first step, catalyzed primarily by magnesium ion A and its ligands, a water molecule attacks the scissile phosphate. Consistent with thiol-substitution experiments, a water proton is transferred to the downstream phosphate group. The transient phosphorane formed as a result of this nucleophilic attack decays by breaking the bond between the phosphate and the ribose oxygen. In the resulting intermediate, the dissociated but unprotonated leaving group forms an alkoxide coordinated to magnesium ion B. In a second step, the reaction is completed by protonation of the leaving group, with a neutral Asp132 as a likely proton donor. The overall reaction barrier of ～15 kcal mol(-1), encountered in the first step, together with the cost of protonating Asp132, is consistent with the slow measured rate of ～1-100/min. The two-step mechanism is also consistent with the bell-shaped pH dependence of the reaction rate. The nonmonotonic relative motion of the magnesium ions along the reaction pathway agrees with X-ray crystal structures. Proton-transfer reactions and changes in the metal ion coordination emerge as central factors in the RNA-cleavage reaction.