Atomistic details of the associative phosphodiester cleavage in human ribonuclease H

Mechanism of Nucleic Acid Phosphodiester Bond Cleavage by Human Endonuclease V: MD and QM/MM Calculations Reveal a Versatile Metal Dependence

Human endonuclease V (EndoV) catalytically removes deaminated nucleobases by cleaving the phosphodiester bond as part of RNA metabolism. Despite being implicated in several diseases (cancers, cardiovascular diseases, and neurological disorders) and potentially being a useful tool in biotechnology, details of the human EndoV catalytic pathway remain unclear due to limited experimental information beyond a crystal structure of the apoenzyme and select mutational data. Since a mechanistic understanding is critical for further deciphering the central roles and expanding applications of human EndoV in medicine and biotechnology, molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations were used to unveil the atomistic details of the catalytic pathway. Due to controversies surrounding the number of metals required for nuclease activity, enzyme-substrate models with different numbers of active site metals and various metal-substrate binding configurations were built based on structural data for other nucleases. Subsequent MD simulations revealed the structure and stability of the human EndoV-substrate complex for a range of active site metal binding architectures. Four unique pathways were then characterized using QM/MM that vary in metal number (one versus two) and modes of substrate coordination [direct versus indirect (water-mediated)], with several mechanisms being fully consistent with experimental structural, kinetic, and mutational data for related nucleases, including members of the EndoV family. Beyond uncovering key roles for several active site amino acids (D240 and K155), our calculations highlight that while one metal is essential for human EndoV activity, the enzyme can benefit from using two metals due to the presence of two suitable metal binding sites. By directly comparing one- versus two-metal-mediated P-O bond cleavage reactions within the confines of the same active site, our work brings a fresh perspective to the "number of metals" controversy.

Is Metal Stabilization of the Leaving Group Required or Can Lysine Facilitate Phosphodiester Bond Cleavage in Nucleic Acids? A Computational Study of EndoV

Endonuclease V (EndoV) is a single-metal-dependent enzyme that repairs deaminated DNA nucleobases in cells by cleaving the phosphodiester bond, and this enzyme has proven to be a powerful tool in biotechnology and medicine. The catalytic mechanism used by EndoV must be understood to design new disease detection and therapeutic solutions and further exploit the enzyme in interdisciplinary applications. This study has used a mixed molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) approach to compare eight distinct catalytic pathways and provides the first proposed mechanism for bacterial EndoV. The calculations demonstrate that mechanisms involving either direct or indirect metal coordination to the leaving group of the substrate previously proposed for other nucleases are unlikely for EndoV, regardless of the general base (histidine, aspartate, and substrate phosphate moiety). Instead, distinct catalytic pathways are characterized for EndoV that involve K139 stabilizing the leaving group, a metal-coordinated water stabilizing the transition structure, and either H214 or a substrate phosphate group activating the water nucleophile. In silico K139A and H214A mutational results support the newly proposed roles of these residues. Although this is a previously unseen combination of general base, general acid, and metal-binding architecture for a one-metal-dependent endonuclease, our proposed catalytic mechanisms are fully consistent with experimental kinetic, structural, and mutational data. In addition to substantiating a growing body of literature, suggesting that one metal is enough to catalyze P-O bond cleavage in nucleic acids, this new fundamental understanding of the catalytic function will promote the exploration of new and improved applications of EndoV.

Controlled sulfur-based engineering confers mouldability to phosphorothioate antisense oligonucleotides

Phosphorothioates (PS) have proven their effectiveness in the area of therapeutic oligonucleotides with applications spanning from cancer treatment to neurodegenerative disorders. Initially, PS substitution was introduced for the antisense oligonucleotides (PS ASOs) because it confers an increased nuclease resistance meanwhile ameliorates cellular uptake and in-vivo bioavailability. Thus, PS oligonucleotides have been elevated to a fundamental asset in the realm of gene silencing therapeutic methodologies. But, despite their wide use, little is known on the possibly different structural changes PS-substitutions may provoke in DNA·RNA hybrids. Additionally, scarce information and significant controversy exists on the role of phosphorothioate chirality in modulating PS properties. Here, through comprehensive computational investigations and experimental measurements, we shed light on the impact of PS chirality in DNA-based antisense oligonucleotides; how the different phosphorothioate diastereomers impact DNA topology, stability and flexibility to ultimately disclose pro-Sp S and pro-Rp S roles at the catalytic core of DNA Exonuclease and Human Ribonuclease H; two major obstacles in ASOs-based therapies. Altogether, our results provide full-atom and mechanistic insights on the structural aberrations PS-substitutions provoke and explain the origin of nuclease resistance PS-linkages confer to DNA·RNA hybrids; crucial information to improve current ASOs-based therapies.

Stimuli-responsive nanoformulations for CRISPR-Cas9 genome editing

The CRISPR-Cas9 technology has changed the landscape of genome editing and has demonstrated extraordinary potential for treating otherwise incurable diseases. Engineering strategies to enable efficient intracellular delivery of CRISPR-Cas9 components has been a central theme for broadening the impact of the CRISPR-Cas9 technology. Various non-viral delivery systems for CRISPR-Cas9 have been investigated given their favorable safety profiles over viral systems. Many recent efforts have been focused on the development of stimuli-responsive non-viral CRISPR-Cas9 delivery systems, with the goal of achieving efficient and precise genome editing. Stimuli-responsive nanoplatforms are capable of sensing and responding to particular triggers, such as innate biological cues and external stimuli, for controlled CRISPR-Cas9 genome editing. In this Review, we overview the recent advances in stimuli-responsive nanoformulations for CRISPR-Cas9 delivery, highlight the rationale of stimuli and formulation designs, and summarize their biomedical applications.

Study of the role of Mg2+ in dsRNA processing mechanism by bacterial RNase III through QM/MM simulations

The ribonuclease III (RNase III) cleaves dsRNA in specific positions generating mature RNAs. RNase III enzymes play important roles in RNA processing, post-transcriptional gene expression, and defense against viral infection. The enzyme's active site contains Mg²⁺ ions bound by a network of acidic residues and water molecules, but there is a lack of information about their specific roles. In this work, multiple steered molecular dynamics simulations at QM/MM level were performed to explore the hydrolysis reaction carried out by the enzyme. Free energy profiles modifying the features of the active site are obtained and the role of Mg²⁺ ions, the solvent molecules and the residues of the active site are discussed in detail. Our results show that Mg²⁺ ions carry out different roles in the hydrolysis process positioning the substrate for the attack from a coordinated nucleophile and activating it to perform hydrolysis reaction, cleaving the dsRNA backbone in a S_N2 substitution. In addition, water molecules present in the active site lower the energy barrier of the process. RNase III hydrolyzes dsRNA to generate mature RNAs. For this purpose, its active site contains Mg²⁺ which has an important role during the reaction. Results show that the Mg²⁺ activates the solvent molecule that produces the nucleophilic attack and the surrounding waters contribute significantly to the hydrolysis process.

Structure and Dynamics of the CRISPR-Cas9 Catalytic Complex

CRISPR-Cas9 is a bacterial immune system with exciting applications for genome editing. In spite of extensive experimental characterization, the active site chemistry of the RuvC domain-which performs DNA cleavages-has remained elusive. Its knowledge is key for structure-based engineering aimed at improving DNA cleavages. Here, we deliver an in-depth characterization by using quantum-classical (QM/MM) molecular dynamics (MD) simulations and a Gaussian accelerated MD method, coupled with bioinformatics analysis. We disclose a two-metal aided architecture in the RuvC active site, which is poised to operate DNA cleavages, in analogy with other DNA/RNA processing enzymes. The conformational dynamics of the RuvC domain further reveals that an "arginine finger" stably contacts the scissile phosphate, with the function of stabilizing the active complex. Remarkably, the formation of a catalytically competent state of the RuvC domain is only observed upon the conformational activation of the other nuclease domain of CRISPR-Cas9-i.e., the HNH domain-such allowing concerted cleavages of double stranded DNA. This structure is in agreement with the available experimental data and remarkably differs from previous models based on classical mechanics, demonstrating also that only quantum mechanical simulations can accurately describe the metal-aided active site in CRISPR-Cas9. This fully catalytic structure-in which both the HNH and RuvC domains are prone to perform DNA cleavages-constitutes a stepping-stone for understanding DNA cleavage and specificity. It calls for novel experimental verifications and offers the structural foundations for engineering efforts aimed at improving the genome editing capability of CRISPR-Cas9.

Two symmetric arginine residues play distinct roles in Thermus thermophilus Argonaute DNA guide strand-mediated DNA target cleavage

Bacterium Thermus thermophilus Argonaute (Ago; TtAgo) is a prokaryotic Ago (pAgo) that acts as the host defense against the uptake and propagation of foreign DNA by catalyzing the DNA cleavage reaction. The TtAgo active site consists of a plugged-in glutamate finger with two arginine residues (R545 and R486) located symmetrically around it. An interesting challenge is to understand how they can collaboratively facilitate enzymatic catalysis. In Kluyveromyces polysporus Ago, a eukaryotic Ago, the evolutionarily symmetrical residues are arginine and histidine, both of which function to stabilize the plugged-in catalytic tetrad conformation. Surprisingly, our simulation results indicated that, in TtAgo, only R545 is involved in the cleavage reaction by serving as a critical structural anchor to stabilize the catalytic tetrad Asp-Glu-Asp-Asp that is completed by the insertion of the glutamate finger, whereas R486 is not involved in target cleavage. The TtAgo-mediated target DNA cleavage occurs in a substrate-assisted mechanism, in which the pro-Rp (Rp, a tetrahedral phosphorus center with "R-type" chirality) oxygen of scissile phosphate acts as a general base to activate the nucleophilic water. Our unexpected theoretical findings on distinct roles played by R545 and R486 in TtAgo catalysis have been validated by single-point site-mutagenesis experiments, wherein the target cleavage is abolished for all mutants of R545. In sharp contrast, the cleavage activity is maintained for all mutants of R486. Our work provides mechanistic insights on the catalytic specificity of Ago proteins and could facilitate the design of new gene-editing tools in the long term.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Catalytic metal ions and enzymatic processing of DNA and RNA

CONSPECTUS: Two-metal-ion-dependent nucleases cleave the phosphodiester bonds of nucleic acids via the two-metal-ion (2M) mechanism. Several high-resolution X-ray structures portraying the two-metal-aided catalytic site, together with mutagenesis and kinetics studies, have demonstrated a functional role of the ions for catalysis in numerous metallonucleases. Overall, the experimental data confirm the general mechanistic hypothesis for 2M-aided phosphoryl transfer originally reported by Steitz and Steitz ( Proc. Natl. Acad. Sci. U.S.A. 1993 , 90 ( 14 ), 6498 - 6502 ). This seminal paper proposed that one metal ion favors the formation of the nucleophile, while the nearby second metal ion facilitates leaving group departure during RNA hydrolysis. Both metals were suggested to stabilize the enzymatic transition state. Nevertheless, static X-ray structures alone cannot exhaustively unravel how the two ions execute their functional role along the enzymatic reaction during processing of DNA or RNA strands when moving from reactants to products, passing through metastable intermediates and high-energy transition states. In this Account, we discuss the role of multiscale molecular simulations in further disclosing mechanistic insights of 2M-aided catalysis for two prototypical enzymatic targets for drug discovery, namely, ribonuclease H (RNase H) and type II topoisomerase (topoII). In both examples, first-principles molecular simulations, integrated with structural data, emphasize a cooperative motion of the bimetal motif during catalysis. The coordinated motion of both ions is crucial for maintaining a flexible metal-centered structural architecture exquisitely tailored to accommodate the DNA or RNA sugar-phosphate backbone during phosphodiester bond cleavage. Furthermore, our analysis of RNase H and the N-terminal domain (PAN) of influenza polymerase shows that classical molecular dynamics simulations coupled with enhanced sampling techniques have contributed to describe the modulatory effect of metal ion concentration and metal uptake on the 2M mechanism and efficiency. These aspects all point to the emerging and intriguing role of additional adjacent ions potentially involved in the modulation of phosphoryl transfer reactions and enzymatic turnover in 2M-catalysis, as recently observed experimentally in polymerase η and homing endonuclease I-DmoI. These computational results, integrated with experimental findings, describe and reinforce the nascent concept of a functional and cooperative dynamics of the catalytic metal ions during the 2M-dependent enzymatic processing of DNA and RNA. Encouraged by the insights provided by computational approaches, we foresee further experiments that will feature the functional and joint dynamics of the catalytic metal ions for nucleic acid processing. This could impact the de novo design of artificial metallonucleases and the rational design of potent metal-chelating inhibitors of pharmaceutically relevant enzymes.

Visualizing phosphodiester-bond hydrolysis by an endonuclease

The enzymatic hydrolysis of DNA phosphodiester bonds has been widely studied, but the chemical reaction has not yet been observed. Here we follow the generation of a DNA double-strand break (DSB) by the Desulfurococcus mobilis homing endonuclease I-DmoI, trapping sequential stages of a two-metal-ion cleavage mechanism. We captured intermediates of the different catalytic steps, and this allowed us to watch the reaction by 'freezing' multiple states. We observed the successive entry of two metals involved in the reaction and the arrival of a third cation in a central position of the active site. This third metal ion has a crucial role, triggering the consecutive hydrolysis of the targeted phosphodiester bonds in the DNA strands and leaving its position once the DSB is generated. The multiple structures show the orchestrated conformational changes in the protein residues, nucleotides and metals during catalysis.

Evidence from molecular dynamics simulations of conformational preorganization in the ribonuclease H active site

Ribonuclease H1 (RNase H) enzymes are well-conserved endonucleases that are present in all domains of life and are particularly important in the life cycle of retroviruses as domains within reverse transcriptase. Despite extensive study, especially of the E. coli homolog, the interaction of the highly negatively charged active site with catalytically required magnesium ions remains poorly understood. In this work, we describe molecular dynamics simulations of the E. coli homolog in complex with magnesium ions, as well as simulations of other homologs in their apo states. Collectively, these results suggest that the active site is highly rigid in the apo state of all homologs studied and is conformationally preorganized to favor the binding of a magnesium ion. Notably, representatives of bacterial, eukaryotic, and retroviral RNases H all exhibit similar active-site rigidity, suggesting that this dynamic feature is only subtly modulated by amino acid sequence and is primarily imposed by the distinctive RNase H protein fold.

Calcium inhibition of ribonuclease H1 two-metal ion catalysis

Most phosphate-processing enzymes require Mg(2+) as a cofactor to catalyze nucleotide cleavage and transfer reactions. Ca(2+) ions inhibit many of these enzymatic activities, despite Ca(2+) and Mg(2+) having comparable binding affinities and overall biological abundances. Here we study the molecular details of the calcium inhibition mechanism for phosphodiester cleavage, an essential reaction in the metabolism of nucleic acids and nucleotides, by comparing Ca(2+)- and Mg(2+) catalyzed reactions. We study the functional roles of the specific metal ion sites A and B in enabling the catalytic cleavage of an RNA/DNA hybrid substrate by B. halodurans ribonuclease (RNase) H1 using hybrid quantum-mechanics/molecular mechanics (QM/MM) free energy calculations. We find that Ca(2+) substitution of either of the two active-site Mg(2+) ions substantially increases the height of the reaction barrier and thereby abolishes the catalytic activity. Remarkably, Ca(2+) at the A site is inactive also in Mg(2+)-optimized active-site structures along the reaction path, whereas Mg(2+) substitution recovers activity in Ca(2+)-optimized structures. Geometric changes resulting from Ca(2+) substitution at metal ion site A may thus be a secondary factor in the loss of catalytic activity. By contrast, at metal ion site B geometry plays a more important role, with only a partial recovery of activity after Mg(2+) substitution in Ca(2+)-optimized structures. Ca(2+)-substitution also leads to a change in mechanism, with deprotonation of the water nucleophile requiring a closer approach to the scissile phosphate, which in turn increases the barrier. As a result, Ca(2+) is less efficient in activating the water. As a likely cause for the different reactivities of Mg(2+) and Ca(2+) ions in site A, we identify differences in charge transfer to the ions and the associated decrease in the pKa of the oxygen nucleophile attacking the phosphate group.

Novel complex MAD phasing and RNase H structural insights using selenium oligonucleotides

The crystal structures of protein-nucleic acid complexes are commonly determined using selenium-derivatized proteins via MAD or SAD phasing. Here, the first protein-nucleic acid complex structure determined using selenium-derivatized nucleic acids is reported. The RNase H-RNA/DNA complex is used as an example to demonstrate the proof of principle. The high-resolution crystal structure indicates that this selenium replacement results in a local subtle unwinding of the RNA/DNA substrate duplex, thereby shifting the RNA scissile phosphate closer to the transition state of the enzyme-catalyzed reaction. It was also observed that the scissile phosphate forms a hydrogen bond to the water nucleophile and helps to position the water molecule in the structure. Consistently, it was discovered that the substitution of a single O atom by a Se atom in a guide DNA sequence can largely accelerate RNase H catalysis. These structural and catalytic studies shed new light on the guide-dependent RNA cleavage.

Structures, dynamics, and stabilities of fully modified locked nucleic acid (β-D-LNA and α-L-LNA) duplexes in comparison to pure DNA and RNA duplexes

Locked nucleic acid (LNA) is a chemical modification which introduces a -O-CH2- linkage in the furanose sugar of nucleic acids and blocks its conformation in a particular state. Two types of modifications, namely, 2'-O,4'-C-methylene-β-D-ribofuranose (β-D-LNA) and 2'-O,4'-C-methylene-α-L-ribofuranose (α-L-LNA), have been shown to yield RNA and DNA duplex-like structures, respectively. LNA modifications lead to increased melting temperatures of DNA and RNA duplexes, and have been suggested as potential therapeutic agents in antisense therapy. In this study, molecular dynamics (MD) simulations were performed on fully modified LNA duplexes and pure DNA and RNA duplexes sharing a similar sequence to investigate their structure, stabilities, and solvation properties. Both LNA duplexes undergo unwinding of the helical structure compared to the pure DNA and RNA duplexes. Though the α-LNA substituent has been proposed to mimic deoxyribose sugar in its conformational properties, the fully modified duplex was found to exhibit unique structural and dynamic properties with respect to the other three nucleic acid structures. Free energy calculations accurately capture the enhanced stabilization of the LNA duplex structures compared to DNA and RNA molecules as observed in experiments. π-stacking interaction between bases from complementary strands is shown to be one of the contributors to enhanced stabilization upon LNA substitution. A combination of two factors, namely, nature of the -O-CH2- linkage in the LNAs vs their absence in the pure duplexes and similar conformations of the sugar rings in DNA and α-LNA vs the other two, is suggested to contribute to the stark differences among the four duplexes studied here in terms of their structural, dynamic, and energetic properties.

Evidence for a dual functional role of a conserved histidine in RNA·DNA heteroduplex cleavage by human RNase H1

Ribonuclease H1 is a conserved enzyme that cleaves the RNA strand of RNA·DNA heteroduplexes and has important functions in the nuclear and mitochondrial compartments. The therapeutic action of antisense oligodeoxynucleotides involves the recruitment of RNase H1 to cleave disease-relevant RNA targets. Recombinant human (Hs) RNase H1 was purified from a bacterial expression host, and conditions were identified that provided optimal oligonucleotide-directed RNA cleavage in vitro. Hs-RNase H1 exhibits optimal catalytic activity in pH 7.5 HEPES buffer and a salt (KCl) concentration of ~ 100-150 mm. Mg(2+) best supports Hs-RNase H1 with an optimal concentration of 10 mm, but at higher concentrations inhibits enzyme activity. Mn(2+) and Co(2+) also support catalytic activity, while Ni(2+) and Zn(2+) exhibit only modest activities as cofactors. The optimized assay was used to show that an antisense oligonucleotide, added in substoichiometric amounts to initiate RNA cleavage, supports up to 30 rounds of reaction in 30 min. Mutation to alanine of the conserved histidine at position 264 causes an ~ 100-fold decrease in k(cat) under multiple-turnover conditions, but does not alter K(m) . Under single-turnover conditions, the H264A mutant exhibits a 12-fold higher exponential time constant for substrate cleavage. The defective activity of the H264A mutant is not rescued in either assay condition by higher Mg(2+) concentrations. These data implicate the H264 side chain in phosphodiester hydrolysis as well as in product release, and are consistent with a proposed model in which the H264 side chain interacts with a divalent metal ion to support catalysis.

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 processes. We stress on the advantages and problems of the different computational strategies and their applicability in different cases.

Theoretical investigation of the enzymatic phosphoryl transfer of β-phosphoglucomutase: revisiting both steps of the catalytic cycle

Enzyme catalyzed phosphate transfer is a part of almost all metabolic processes. Such reactions are of central importance for the energy balance in all organisms and play important roles in cellular control at all levels. Mutases transfer a phosphoryl group while nucleases cleave the phosphodiester linkages between two nucleotides. The subject of our present study is the Lactococcus lactis β-phosphoglucomutase (β-PGM), which effectively catalyzes the interconversion of β-D-glucose-1-phosphate (β-G1P) to β-D-glucose-6-phosphate (β-G6P) and vice versa via stabile intermediate β-D-glucose-1,6-(bis)phosphate (β-G1,6diP) in the presence of Mg(2+). In this paper we revisited the reaction mechanism of the phosphoryl transfer starting from the bisphosphate β-G1,6diP in both directions (toward β-G1P and β-G6P) combining docking techniques and QM/MM theoretical method at the DFT/PBE0 level of theory. In addition we performed NEB (nudged elastic band) and free energy calculations to optimize the path and to identify the transition states and the energies involved in the catalytic cycle. Our calculations reveal that both steps proceed via dissociative pentacoordinated phosphorane, which is not a stabile intermediate but rather a transition state. In addition to the Mg(2+) ion, Ser114 and Lys145 also play important roles in stabilizing the large negative charge on the phosphate through strong coordination with the phosphate oxygens and guiding the phosphate group throughout the catalytic process. The calculated energy barrier of the reaction for the β-G1P to β-G1,6diP step is only slightly higher than for the β-G1,6diP to β-G6P step (16.10 kcal mol(-1) versus 15.10 kcal mol(-1)) and is in excellent agreement with experimental findings (14.65 kcal mol(-1)).

Density functional calculations on the effect of sulfur substitution for 2'-hydroxypropyl-p-nitrophenyl phosphate: C-O vs. P-O bond cleavage

Density functional theory calculations have been used to investigate the intra-molecular attack of 2'-hydroxypropyl-p-nitrophenyl phosphate (HPpNP) and its analogous compound 2-thiouridyl-p-nitrophenyl phosphate (s-2'pNP). Bulk solvent effect has been tested at the geometry optimization level with the polarized continuum model. It is found that the P-path involving the intra-molecular attack at the phosphorus atom and C-path involving the attack at the beta carbon atom proceed through the S(N)2-type mechanism for HPpNP and s-2'pNP. The calculated results indicate that the P-path with the free energy barrier of about 11 kcal/mol is more accessible than the C-path for the intra-molecular attack of HPpNP, which favors the formation of the five-membered phosphate diester. While for s-2'pNP, the C-path with the free energy barrier of about 21 kcal/mol proceeds more favorably than the P-path. The calculated energy barriers of the favorable pathways for HPpNP and s-2'pNP are both in agreement with the experimental results.

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.

Nucleotide docking: prediction of reactant state complexes for ribonuclease enzymes

Ribonuclease enzymes (RNases) play key roles in the maturation and metabolism of all RNA molecules. Computational simulations of the processes involved can help to elucidate the underlying enzymatic mechanism and is often employed in a synergistic approach together with biochemical experiments. Theoretical calculations require atomistic details regarding the starting geometries of the molecules involved, which, in the absence of crystallographic data, can only be achieved from computational docking studies. Fortunately, docking algorithms have improved tremendously in recent years, so that reliable structures of enzyme-ligand complexes can now be successfully obtained from computation. However, most docking programs are not particularly optimized for nucleotide docking. In order to assist our studies on the cleavage of RNA by the two most important ribonuclease enzymes, RNase A and RNase H, we evaluated four docking tools-MOE2009, Glide 5.5, QXP-Flo+0802, and Autodock 4.0-for their ability to simulate complexes between these enzymes and RNA oligomers. To validate our results, we analyzed the docking results with respect to the known key interactions between the protein and the nucleotide. In addition, we compared the predicted complexes with X-ray structures of the mutated enzyme as well as with structures obtained from previous calculations. In this manner, we were able to prepare the desired reaction state complex so that it could be used as the starting structure for further DFT/B3LYP QM/MM reaction mechanism studies.