Understanding the effect of magnesium ion concentration on the catalytic activity of ribonuclease H through computation: does a third metal binding site modulate endonuclease catalysis?

Third Metal Ion Dictates the Catalytic Activity of the Two-Metal-Ion Pre-Ribosomal RNA-Processing Machinery

Nucleic acid processing enzymes use a two-Mg2+-ion motif to promote the formation and cleavage of phosphodiester bonds. Yet, recent evidence demonstrates the presence of spatially conserved second-shell cations surrounding the catalytic architecture of proteinaceous and RNA-dependent enzymes. The RNase mitochondrial RNA processing (MRP) complex, which cleaves the ribosomal RNA (rRNA) precursor at the A3 cleavage site to yield mature 5'-end of 5.8S rRNA, hosts in the catalytic core one atypically-located Mg2+ ion, in addition to the ions forming the canonical catalytic motif. Here, we employ biased quantum classical molecular dynamics simulations of RNase MRP to discover that the third Mg2+ ion inhibits the catalytic process. Instead, its displacement in favour of a second-shell monovalent K+ ion propels phosphodiester bond cleavage by enabling the formation of a specific hydrogen bonding network that mediates the essential proton transfer step. This study points to a direct involvement of a transient K+ ion in the catalytic cleavage of the phosphodiester bond and implicates cation trafficking as a general mechanism in nucleic acid processing enzymes and ribozymes.

Optimizing Messenger RNA Analysis Using Ultra-Wide Pore Size Exclusion Chromatography Columns

Biopharmaceutical products, in particular messenger ribonucleic acid (mRNA), have the potential to dramatically improve the quality of life for patients suffering from respiratory and infectious diseases, rare genetic disorders, and cancer. However, the quality and safety of such products are particularly critical for patients and require close scrutiny. Key product-related impurities, such as fragments and aggregates, among others, can significantly reduce the efficacy of mRNA therapies. In the present work, the possibilities offered by size exclusion chromatography (SEC) for the characterization of mRNA samples were explored using state-of-the-art ultra-wide pore columns with average pore diameters of 1000 and 2500 Å. Our investigation shows that a column with 1000 Å pores proved to be optimal for the analysis of mRNA products, whatever the size between 500 and 5000 nucleotides (nt). We also studied the influence of mobile phase composition and found that the addition of 10 mM magnesium chloride (MgCl2) can be beneficial in improving the resolution and recovery of large size variants for some mRNA samples. We demonstrate that caution should be exercised when increasing column length or decreasing the flow rate. While these adjustments slightly improve resolution, they also lead to an apparent increase in the amount of low-molecular-weight species (LMWS) and monomer peak tailing, which can be attributed to the prolonged residence time inside the column. Finally, our optimal SEC method has been successfully applied to a wide range of mRNA products, ranging from 1000 to 4500 nt in length, as well as mRNA from different suppliers and stressed/unstressed samples.

Impact of hierarchical water dipole orderings on the dynamics of aqueous salt solutions

Ions exhibit highly ion-specific complex behaviours when solvated in water, which remains a mystery despite the fundamental importance of ion solvation in nature, science, and technology. Here we explain these ion-specific properties by the ion-induced hierarchical dipolar, translational, and bond-orientational orderings of ion hydration shell under the competition between ion-water electrostatic interactions and inter-water hydrogen bonding. We first characterise this competition by a new length λHB(q), explaining the ion-specific effects on solution dynamics. Then, by continuously tuning ion size and charge, we find that the bond-orientational order of the ion hydration shell highly develops for specific ion size and charge combinations. This ordering drastically stabilises the hydration shell; its degree changes the water residence time around ions by 11 orders of magnitude for main-group ions. These findings are fundamental to ionic processes in aqueous solutions, providing a physical principle for electrolyte design and application.

Metal-binding and folding thermodynamics of Escherichia coli ribonuclease HI related to its catalytic function

Escherichia coli ribonuclease HI (RNH) hydrolyzes the RNA strands of RNA/DNA hybrids in the presence of Mg²⁺ at the highest level, relative to other metal ions. The Mg²⁺ binding affinity was 8.39 × 10³ M^-1, which was lower than those of other metal ions. The low-affinity binder can express the maximum catalytic activity of RNH. The stability of RNH increased with increasing metal ion concentration, except for Zn²⁺. The thermodynamic origin for enhancing the stability of RNH with Mg²⁺ was more favorable entropy compared to those with other metal ions, indicating that Mg²⁺ binding changes the RNH structure while maintaining flexibility. Upon H124A mutation, the metal ion binding affinities decreased for Mn²⁺ and Zn²⁺ to a relatively large extent. The present thermodynamic analyses provide information on the structural dynamics of RNH with metal ion exchangeable binding, which can reasonably explain the metal-ion-dependent catalytic activity.

The metal dependence of single-metal mediated phosphodiester bond cleavage: a QM/MM study of a multifaceted human enzyme

Nucleases catalyze the cleavage of phosphodiester bonds in nucleic acids using a range of metal cofactors. Although it is well accepted that many nucleases rely on two metal ions, the one-metal mediated pathway is debated. Furthermore, one-metal mediated nucleases maintain activity in the presence of many different metals, but the underlying reasons for this broad metal specificity are unknown. The human apurinic/apyrimidinic endonuclease (APE1), which plays a key role in DNA repair, transcription regulation, and gene expression, is a prototypical example of a one-metal dependent nuclease. Although Mg²⁺ is the native metal cofactor, APE1 remains catalytically active in the presence of several metals, with the rate decreasing as Mg²⁺ > Mn²⁺ > Ni²⁺ > Zn²⁺, while Ca²⁺ completely abolished the activity. The present work uses quantum mechanics-molecular mechanics techniques to map APE1-facilitated phosphodiester bond hydrolysis in the presence of these metals. The structural differences in stationary points along the reaction pathway shed light on the interplay between several factors that allow APE1 to remain catalytically active for various metals, with the trend in the barrier heights correlating with the experimentally reported APE1 catalytic activity. In contrast, Ca²⁺ significantly changes the metal coordination and active site geometry, and thus completely inhibits catalysis. Our work thereby provides support for the controversial single-metal mediated phosphodiester bond cleavage and clarifies uncertainties regarding the role of the metal and metal identity in this important reaction. This information is key for future medicinal and biotechnological applications including disease diagnosis and treatment, and protein engineering.

The catalytic mechanism, metal dependence, substrate specificity, and biodiversity of ribonuclease H

Ribonucleoside monophosphates are inevitably misincorporated into the DNA genome inside cells, and they need to be excised to avoid chromosome instability. Ribonucleases H (RNases H) are enzymes that specifically hydrolyze the RNA strand of RNA/DNA hybrids or the RNA moiety from DNA containing a stretch of RNA, they therefore are required for DNA integrity. Extensive studies have drawn a mostly clear picture of the mechanisms of RNase H catalysis, but some questions are still lacking definitive answers. This review summarizes three alternative models of RNase H catalysis. The two-metal model is prevalent, but a three-metal model suggests the involvement of a third cation in catalysis. Apparently, the mechanisms underlying metal-dependent hydrolyzation are more complicated than initially thought. We also discuss the metal choices of RNases H and analyze how chemically similar cations function differently. Substrate and cleavage-site specificities vary among RNases H, and this is explicated in detail. An intriguing phenomenon is that organisms have diverse RNase H combinations, which may provide important hints to how rnh genes were transferred during evolution. Whether RNase H is essential for cellular growth, a key question in the study of in vivo functions, is also discussed. This article may aid in understanding the mechanisms underlying RNase H and in developing potentially promising applications of it.

DNA-gold nanoprobe-based integrated biosensing technology for non-invasive liquid biopsy of serum miRNA: A new frontier in prostate cancer diagnosis

The low specificity of prostate-specific antigen contributes to overdiagnosis and ov ertreatment of prostate cancer (PCa) patients. Hence, there is an urgent need for inclusive diagnostic platforms that could improve the diagnostic accuracy of PCa. Dysregulated miRNAs are closely associated with the progression and recurrence and have emerged as promising diagnostic and prognostic biomarkers for PCa. Nevertheless, simple, rapid, and ultrasensitive quantification of serum miRNAs is highly challenging. This study designed, synthesized, and demonstrated the practicability of DNA-linked gold nanoprobes (DNA-AuNPs) for the single-step quantification of miR-21/miR-141/miR-375. In preclinical study, the assay differented PCa Pten conditional knockout (Pten^cKO) mice compared to their age-matched Pten wild-type (Pten^WT) control mice. In human sera, receiver operating characteristic (ROC) curve-based correlation analyses revealed clear discrimination between PCa patients from normal healthy controls using training and validation sets. Overall, we established integrated nano-biosensing technology for the PCR-free, non-invasive liquid biopsies of multiple miRNAs for PCa diagnosis.

Sulfur modification in natural RNA and therapeutic oligonucleotides

Sulfur modifications have been discovered on both DNA and RNA. Sulfur substitution of oxygen atoms at nucleobase or backbone locations in the nucleic acid framework led to a wide variety of sulfur-modified nucleosides and nucleotides. While the discovery, regulation and functions of DNA phosphorothioate (PS) modification, where one of the non-bridging oxygen atoms is replaced by sulfur on the DNA backbone, are important topics, this review focuses on the sulfur modification in natural cellular RNAs and therapeutic nucleic acids. The sulfur modifications on RNAs exhibit diversity in terms of modification location and cellular function, but the various sulfur modifications share common biosynthetic strategies across RNA species, cell types and domains of life. The first section reviews the post-transcriptional sulfur modifications on nucleobases with an emphasis on thiouridine on tRNA and phosphorothioate modification on RNA backbones, as well as the functions of the sulfur modifications on different species of cellular RNAs. The second section reviews the biosynthesis of different types of sulfur modifications and summarizes the general strategy for the biosynthesis of sulfur-containing RNA residues. One of the main goals of investigating sulfur modifications is to aid the genomic drug development pipeline and enhance our understandings of the rapidly growing nucleic acid-based gene therapies. The last section of the review focuses on the current drug development strategies employing sulfur substitution of oxygen atoms in therapeutic RNAs.

Ruthenium(II) Complex-Based G-quadruplex DNA Selective Luminescent 'Light-up' Probe for RNase H Activity Detection

A bis-heteroleptic ruthenium(II) complex, 1[PF₆ ]₂ of benzothiazole amide substituted 2,2'-bipyridine ligand (bmbbipy) has been synthesized for the selective detection of G-quadruplex (GQ) DNA and luminescence-assay-based RNase H activity monitoring. Compound 1[PF₆ ]₂ exhibited aggregation-caused quenching (ACQ) in water. Aggregate formation was supported by DLS, UV-vis, and ¹ H NMR spectroscopy results, and the morphology of aggregated particles was witnessed by SEM and TEM. 1[PF₆ ]₂ acted as an efficient GQ DNA-selective luminescent light-up probe over single-stranded and double-stranded DNA. The competency of 1[PF₆ ]₂ for selective GQ structure detection was established by PL and CD spectroscopy. For 1[PF₆ ]₂ , the PL light-up is exclusively due to the rigidification of the benzothiazole amide side arm in the presence of GQ-DNA. The interaction between the probe and GQ-DNA was analyzed by molecular docking analysis. The GQ structure detection capability of 1[PF₆ ]₂ was further applied in the luminescent 'off-on' RNase H activity detection. The assay utilized an RNA:DNA hybrid, obtained from 22AG2-RNA and 22AG2-DNA sequences. RNase H solely hydrolyzed the RNA of the RNA:DNA duplex and released G-rich 22AG2-DNA, which was detected via the PL enhancement of 1[PF₆ ]₂ . The selectivity of RNase H activity detection over various other restriction enzymes was also demonstrated.

Controlled Trafficking of Multiple and Diverse Cations Prompts Nucleic Acid Hydrolysis

Recent in crystallo reaction intermediates have detailed how nucleic acid hydrolysis occurs in the RNA ribonuclease H1 (RNase H1), a fundamental metalloenzyme involved in maintaining the human genome. At odds with the previous characterization, these in crystallo structures unexpectedly captured multiple metal ions (K⁺ and Mg²⁺) transiently bound in the vicinity of the two-metal-ion active site. Using multi-microsecond all-atom molecular dynamics and free-energy simulations, we investigated the functional implications of the dynamic exchange of multiple K⁺ and Mg²⁺ ions at the RNase H1 reaction center. We found that such ions are timely positioned at non-overlapping locations near the active site, at different stages of catalysis, being crucial for both reactants' alignment and leaving group departure. We also found that this cation trafficking is tightly regulated by variations of the solution's ionic strength and is aided by two conserved second-shell residues, E188 and K196, suggesting a mechanism for the cations' recruitment during catalysis. These results indicate that controlled trafficking of multi-cation dynamics, opportunely prompted by second-shell residues, is functionally essential to the complex enzymatic machinery of the RNase H1. These findings revise the current knowledge on the RNase H1 catalysis and open new catalytic possibilities for other similar metalloenzymes including, but not limited to, CRISPR-Cas9, group II intron ribozyme and the human spliceosome.

Structural basis for Ca2+-dependent catalysis of a cutinase-like enzyme and its engineering: application to enzymatic PET depolymerization

A cutinase-like enzyme from Saccharomonospora viridis AHK190, Cut190, can depolymerize polyethylene terephthalate (PET). As high activity at approximately 70°C is required for PET depolymerization, structure-based protein engineering of Cut190 was carried out. Crystal structure information of the Cut190 mutants was used for protein engineering and for evaluating the molecular basis of activity and thermal stability. A variety of biophysical methods were employed to unveil the mechanisms underlying the unique features of Cut190, which included the regulation of its activity and thermal stability by Ca²⁺. Ca²⁺ association and dissociation can change the enzyme conformation to regulate catalytic activity. Weak metal-ion binding would be required for the naïve conformational change of Cut190, while maintaining its fluctuation, to “switch” the enzyme on and off. The activity of Cut190 is regulated by the weak Ca²⁺ binding to the specific site, Site 1, while thermal stability is mainly regulated by binding to another Site 2, where a disulfide bond could be introduced to increase the stability. Recent results on the structure-activity relationship of engineered Cut190 are reviewed, including the application for PET depolymerization by enzymes.

Active-Site Models of Streptococcus pyogenes Cas9 in DNA Cleavage State

CRISPR-Cas9 is a powerful tool for target genome editing in living cells. Significant advances have been made to understand how this system cleaves target DNA. HNH is a nuclease domain, which shares structural similarity with the HNH endonuclease characterzied by a beta-beta-alpha-metal fold. Therefore, based on one- and two-metal-ion mechanisms, homology modeling and molecular dynamics (MD) simulation are suitable tools for building an atomic model of Cas9 in the DNA cleavage state. Here, by modeling and MD, we presented an atomic model of SpCas9-sgRNA-DNA complex with the cleavage state. This model shows that the HNH and RuvC conformations resemble their DNA cleavage state where the active-sites in the complex coordinate with DNA, Mg²⁺ ions, and water. Among them, residues D10, E762, H983, and D986 locate at the first shell of the RuvC active-site and interact with the ions directly, residues H982 or/and H985 are general (Lewis) bases, and the coordinated water is located at the positions for nucleophilic attack of the scissile phosphate. Meanwhile, this catalytic model led us to engineer a new SpCas9 variant (SpCas9-H982A + H983D) with reduced off-target effects. Thus, our study provided new mechanistic insights into the CRISPR-Cas9 system in the DNA cleavage state and offered useful guidance for engineering new CRISPR-Cas9 editing systems with improved specificity.

The Protection Role of Magnesium Ions on Coupled Transcription and Translation in Lyophilized Cell-Free System

Cell-free protein synthesis (CFPS) is a promising platform for protein engineering and synthetic biology. The storage of a CFPS system usually involves lyophilization, during which preventing the conformational damage of involved enzymes is critical to the activity. Herein, we report the protection role of magnesium ions on coupled transcription and translation in a lyophilized cell-free system. Mg²⁺ prevents the inactivation of the CFPS system from direct colyophilization of enzymes and substrates (nucleotides, and amino acids), and furthermore activates the CFPS system. We propose two-metal-ion regulation of Mg²⁺: Mg²⁺ (I) acts as an allosteric role for enzymes to prevent the conformational damage of enzymes from direct binding with substrates during lyophilization which locks up inactive enzyme-substrate complex; Mg²⁺ (II) consequently binds to enzymes to activate the CFPS system. Our work provides important implications for maximizing protein yields by using a cell-free system in protein engineering and understanding the functions of Mg²⁺ in biological systems.

Recruiting Mechanism and Functional Role of a Third Metal Ion in the Enzymatic Activity of 5' Structure-Specific Nucleases

Enzymes of the 5′ structure-specific nuclease family are crucial for DNA repair, replication, and recombination. One such enzyme is the human exonuclease 1 (hExo1) metalloenzyme, which cleaves DNA strands, acting primarily as a processive 5′-3′ exonuclease and secondarily as a 5′-flap endonuclease. Recently, in crystallo reaction intermediates have elucidated how hExo1 exerts hydrolysis of DNA phosphodiester bonds. These hExo1 structures show a third metal ion intermittently bound close to the two-metal-ion active site, to which recessed ends or 5′-flap substrates bind. Evidence of this third ion has been observed in several nucleic-acid-processing metalloenzymes. However, there is still debate over what triggers the (un)binding of this transient third ion during catalysis and whether this ion has a catalytic function. Using extended molecular dynamics and enhanced sampling free-energy simulations, we observed that the carboxyl side chain of Glu89 (located along the arch motif in hExo1) flips frequently from the reactant state to the product state. The conformational flipping of Glu89 allows one metal ion to be recruited from the bulk and promptly positioned near the catalytic center. This is in line with the structural evidence. Additionally, our simulations show that the third metal ion assists the departure, through the mobile arch, of the nucleotide monophosphate product from the catalytic site. Structural comparisons of nuclease enzymes suggest that this Glu(Asp)-mediated mechanism for third ion recruitment and nucleic acid hydrolysis may be shared by other 5′ structure-specific nucleases.

A Transient and Flexible Cation-π Interaction Promotes Hydrolysis of Nucleic Acids in DNA and RNA Nucleases

Metal-dependent DNA and RNA nucleases are enzymes that cleave nucleic acids with great efficiency and precision. These enzyme-mediated hydrolytic reactions are fundamental for the replication, repair, and storage of genetic information within the cell. Here, extensive classical and quantum-based free-energy molecular simulations show that a cation-π interaction is transiently formed in situ at the metal core of Bacteriophage-λ Exonuclease (Exo-λ), during catalysis. This noncovalent interaction (Lys131-Tyr154) triggers nucleophile activation for nucleotide excision. Then, our simulations also show the oscillatory dynamics and swinging of the newly formed cation-π dyad, whose conformational change may favor proton release from the cationic Lys131 to the bulk solution, thus restoring the precatalytic protonation state in Exo-λ. Altogether, we report on the novel mechanistic character of cation-π interactions for catalysis. Structural and bioinformatic analyses support that flexible orientation and transient formation of mobile cation-π interactions may represent a common catalytic strategy to promote nucleic acid hydrolysis in DNA and RNA nucleases.

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.

Second-Shell Basic Residues Expand the Two-Metal-Ion Architecture of DNA and RNA Processing Enzymes

Synthesis and scission of phosphodiester bonds in DNA and RNA regulate vital processes within the cell. Enzymes that catalyze these reactions operate mostly via the recognized two-metal-ion mechanism. Our analysis reveals that basic amino acids and monovalent cations occupy structurally conserved positions nearby the active site of many two-metal-ion enzymes for which high-resolution (<3 Å) structures are known, including DNA and RNA polymerases, nucleases such as Cas9, and splicing ribozymes. Integrating multiple-sequence and structural alignments with molecular dynamics simulations, electrostatic potential maps, and mutational data, we found that these elements always interact with the substrates, suggesting that they may play an active role for catalysis, in addition to their electrostatic contribution. We discuss possible mechanistic implications of this expanded two-metal-ion architecture, including inferences on medium-resolution cryoelectron microscopy structures. Ultimately, our analysis may inspire future experiments and strategies for enzyme engineering or drug design to modulate nucleic acid processing.

•

Basic residues in the active site of two-metal-ion enzymes are structurally conserved

•

These residues are also conserved in evolution

•

Mutagenesis suggests these residues may exert an effect on DNA- and RNA processing

•

Our work offers insights into CRISPR/Cas9, spliceosome, and DNA/RNA polymerases

Viral genome packaging terminase cleaves DNA using the canonical RuvC-like two-metal catalysis mechanism

Bacteriophages and large dsDNA viruses encode sophisticated machinery to translocate their DNA into a preformed empty capsid. An essential part of this machine, the large terminase protein, processes viral DNA into constituent units utilizing its nuclease activity. Crystal structures of the large terminase nuclease from the thermophilic bacteriophage G20c show that it is most similar to the RuvC family of the RNase H-like endonucleases. Like RuvC proteins, the nuclease requires either Mn²⁺, Mg²⁺ or Co²⁺ ions for activity, but is inactive with Zn²⁺ and Ca²⁺. High resolution crystal structures of complexes with different metals reveal that in the absence of DNA, only one catalytic metal ion is accommodated in the active site. Binding of the second metal ion may be facilitated by conformational variability, which enables the two catalytic aspartic acids to be brought closer to each other. Structural comparison indicates that in common with the RuvC family, the location of the two catalytic metals differs from other members of the RNase H family. In contrast to a recently proposed mechanism, the available data do not support binding of the two metals at an ultra-short interatomic distance. Thus we postulate that viral terminases cleave DNA by the canonical RuvC-like mechanism.

Development of Site-Specific Mg(2+)-RNA Force Field Parameters: A Dream or Reality? Guidelines from Combined Molecular Dynamics and Quantum Mechanics Simulations

The vital contribution of Mg²⁺ ions to RNA biology is challenging to dissect at the experimental level. This calls for the integrative support of atomistic simulations, which at the classical level are plagued by limited accuracy. Indeed, force fields intrinsically neglect nontrivial electronic effects that Mg²⁺ exerts on its surrounding ligands in varying RNA coordination environments. Here, we present a combined computational study based on classical molecular dynamics (MD) and Density Functional Theory (DFT) calculations, aimed at characterizing (i) the performance of five Mg²⁺ force field (FF) models in RNA systems and (ii) how charge transfer and polarization affect the binding of Mg²⁺ ions in different coordination motifs. As a result, a total of ∼2.5 μs MD simulations (100/200 ns for each run) for two prototypical Mg²⁺-dependent ribozymes showed remarkable differences in terms of populations of inner-sphere coordination site types. Most importantly, complementary DFT calculations unveiled that differences in charge transfer and polarization among recurrent Mg²⁺-RNA coordination motifs are surprisingly small. In particular, the charge of the Mg²⁺ ions substantially remains constant through different coordination sites, suggesting that the common philosophy of developing site-specific Mg²⁺ ion parameters is not in line with the physical origin of the Mg²⁺-RNA MD simulations inaccuracies. Overall, this study constitutes a guideline for an adept use of current Mg²⁺ models and provides novel insights for the rational development of next-generation Mg²⁺ FFs to be employed for atomistic simulations of RNA.

A Self-Activated Mechanism for Nucleic Acid Polymerization Catalyzed by DNA/RNA Polymerases

The enzymatic polymerization of DNA and RNA is the basis for genetic inheritance for all living organisms. It is catalyzed by the DNA/RNA polymerase (Pol) superfamily. Here, bioinformatics analysis reveals that the incoming nucleotide substrate always forms an H-bond between its 3'-OH and β-phosphate moieties upon formation of the Michaelis complex. This previously unrecognized H-bond implies a novel self-activated mechanism (SAM), which synergistically connects the in situ nucleophile formation with subsequent nucleotide addition and, importantly, nucleic acid translocation. Thus, SAM allows an elegant and efficient closed-loop sequence of chemical and physical steps for Pol catalysis. This is markedly different from previous mechanistic hypotheses. Our proposed mechanism is corroborated via ab initio QM/MM simulations on a specific Pol, the human DNA polymerase-η, an enzyme involved in repairing damaged DNA. The structural conservation of DNA and RNA Pols supports the possible extension of SAM to Pol enzymes from the three domains of life.

A Computational Model for Predicting RNase H Domain of Retrovirus

RNase H (RNH) is a pivotal domain in retrovirus to cleave the DNA-RNA hybrid for continuing retroviral replication. The crucial role indicates that RNH is a promising drug target for therapeutic intervention. However, annotated RNHs in UniProtKB database have still been insufficient for a good understanding of their statistical characteristics so far. In this work, a computational RNH model was proposed to annotate new putative RNHs (np-RNHs) in the retroviruses. It basically predicts RNH domains through recognizing their start and end sites separately with SVM method. The classification accuracy rates are 100%, 99.01% and 97.52% respectively corresponding to jack-knife, 10-fold cross-validation and 5-fold cross-validation test. Subsequently, this model discovered 14,033 np-RNHs after scanning sequences without RNH annotations. All these predicted np-RNHs and annotated RNHs were employed to analyze the length, hydrophobicity and evolutionary relationship of RNH domains. They are all related to retroviral genera, which validates the classification of retroviruses to a certain degree. In the end, a software tool was designed for the application of our prediction model. The software together with datasets involved in this paper can be available for free download at https://sourceforge.net/projects/rhtool/files/?source=navbar.

Cooperative motion of a key positively charged residue and metal ions for DNA replication catalyzed by human DNA Polymerase-η

Trans-lesion synthesis polymerases, like DNA Polymerase-η (Pol-η), are essential for cell survival. Pol-η bypasses ultraviolet-induced DNA damages via a two-metal-ion mechanism that assures DNA strand elongation, with formation of the leaving group pyrophosphate (PPi). Recent structural and kinetics studies have shown that Pol-η function depends on the highly flexible and conserved Arg61 and, intriguingly, on a transient third ion resolved at the catalytic site, as lately observed in other nucleic acid-processing metalloenzymes. How these conserved structural features facilitate DNA replication, however, is still poorly understood. Through extended molecular dynamics and free energy simulations, we unravel a highly cooperative and dynamic mechanism for DNA elongation and repair, which is here described by an equilibrium ensemble of structures that connect the reactants to the products in Pol-η catalysis. We reveal that specific conformations of Arg61 help facilitate the recruitment of the incoming base and favor the proper formation of a pre-reactive complex in Pol-η for efficient DNA editing. Also, we show that a third transient metal ion, which acts concertedly with Arg61, serves as an exit shuttle for the leaving PPi. Finally, we discuss how this effective and cooperative mechanism for DNA repair may be shared by other DNA-repairing polymerases.

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.

Synthesis and structures of soluble magnesium and zinc carboxylates containing intramolecular NH⋯O hydrogen bonds in nonpolar solvents

Magnesium and zinc carboxylates containing intramolecular NH⋯O hydrogen bonds showed a fast trans–cis isomerization in nonpolar solvents and were converted into anionic tris(carboxylate)s by the addition of an equimolar ligand.

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.

Ultradeep pyrosequencing and molecular modeling identify key structural features of hepatitis B virus RNase H, a putative target for antiviral intervention

Last-generation nucleoside/nucleotide analogues are potent against hepatitis B virus (HBV) and have a high barrier to resistance. However, delayed responses have been observed in patients previously exposed to other drugs of the same class, long-term resistance is possible, and cure of infection cannot be achieved with these therapies, emphasizing the need for alternative therapeutic approaches. The HBV RNase H represents an interesting target because its enzyme activity is essential to the HBV life cycle. The goal of our study was to characterize the structure of the HBV RNase H by computing a 3-dimensional molecular model derived from E. coli RNase H and analyzing 2,326 sequences of all HBV genotypes available in public databases and 958,000 sequences generated by means of ultradeep pyrosequencing of sequences from a homogenous population of 73 treatment-naive patients infected with HBV genotype D. Our data revealed that (i) the putative 4th catalytic residue displays unexpected variability that could be explained by the overlap of the HBx gene and has no apparent impact on HBV replicative capacity and that (ii) the C-helix-containing basic protrusion, which is required to guide the RNA/DNA heteroduplex into the catalytic site, is highly conserved and bears unique structural properties that can be used to target HBV-specific RNase H inhibitors without cross-species activity. The model shows substantial differences from other known RNases H and paves the way for functional and structural studies as a prerequisite to the development of new inhibitors of the HBV cell cycle specifically targeting RNase H activity.

Conserved structural chemistry for incision activity in structurally non-homologous apurinic/apyrimidinic endonuclease APE1 and endonuclease IV DNA repair enzymes

Non-coding apurinic/apyrimidinic (AP) sites in DNA form spontaneously and as DNA base excision repair intermediates are the most common toxic and mutagenic in vivo DNA lesion. For repair, AP sites must be processed by 5' AP endonucleases in initial stages of base repair. Human APE1 and bacterial Nfo represent the two conserved 5' AP endonuclease families in the biosphere; they both recognize AP sites and incise the phosphodiester backbone 5' to the lesion, yet they lack similar structures and metal ion requirements. Here, we determined and analyzed crystal structures of a 2.4 Å resolution APE1-DNA product complex with Mg(2+) and a 0.92 Å Nfo with three metal ions. Structural and biochemical comparisons of these two evolutionarily distinct enzymes characterize key APE1 catalytic residues that are potentially functionally similar to Nfo active site components, as further tested and supported by computational analyses. We observe a magnesium-water cluster in the APE1 active site, with only Glu-96 forming the direct protein coordination to the Mg(2+). Despite differences in structure and metal requirements of APE1 and Nfo, comparison of their active site structures surprisingly reveals strong geometric conservation of the catalytic reaction, with APE1 catalytic side chains positioned analogously to Nfo metal positions, suggesting surprising functional equivalence between Nfo metal ions and APE1 residues. The finding that APE1 residues are positioned to substitute for Nfo metal ions is supported by the impact of mutations on activity. Collectively, the results illuminate the activities of residues, metal ions, and active site features for abasic site endonucleases.

Designing chemically modified oligonucleotides for targeted gene silencing

Oligonucleotides (ONs), and their chemically modified mimics, are now routinely used in the laboratory as a means to control the expression of fundamentally interesting or therapeutically relevant genes. ONs are also under active investigation in the clinic, with many expressing cautious optimism that at least some ON-based therapies will succeed in the coming years. In this review, we will discuss several classes of ONs used for controlling gene expression, with an emphasis on antisense ONs (AONs), small interfering RNAs (siRNAs), and microRNA-targeting ONs (anti-miRNAs). This review provides a current and detailed account of ON chemical modification strategies for the optimization of biological activity and therapeutic application, while clarifying the biological pathways, chemical properties, benefits, and limitations of oligonucleotide analogs used in nucleic acids research.

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.

Metal and ligand binding to the HIV-RNase H active site are remotely monitored by Ile556

HIV-1 reverse transcriptase (RT) contains a C-terminal ribonuclease H (RH) domain on its p66 subunit that can be expressed as a stable, although inactive protein. Recent studies of several RH enzymes demonstrate that substrate binding plays a major role in the creation of the active site. In the absence of substrate, the C-terminal helix E of the RT RNase H domain is dynamic, characterized by severe exchange broadening of its backbone amide resonances, so that the solution characterization of this region of the protein has been limited. Nuclear magnetic resonance studies of 13C-labeled RH as a function of experimental conditions reveal that the δ1 methyl resonance of Ile556, located in a short, random coil segment following helix E, experiences a large 13C shift corresponding to a conformational change of Ile556 that results from packing of helix E against the central β-sheet. This shift provides a useful basis for monitoring the effects of various ligands on active site formation. Additionally, we report that the RNase H complexes formed with one or both divalent ions can be individually observed and characterized using diamagnetic Zn2+ as a substitute for Mg2+. Ordering of helix E results specifically from the interaction with the lower affinity binding to the A divalent ion site.

Probe Position-Dependent Counterion Dynamics in DNA: Comparison of Time-Resolved Stokes Shift of Groove-Bound to Base-Stacked Probes in the Presence of Different Monovalent Counterions

Time-resolved fluorescence Stokes shifts (TRFSS) of 4',6-diamidino-2-phenylindole (DAPI) inside the minor groove of DNA are measured in the presence of three different monovalent counterions: sodium (Na(+)), rubidium (Rb(+)), and tetrabutylammonium (TBA(+)). Fluorescence up-conversion and time-correlated single photon counting are combined to obtain the time-resolved emission spectra (TRES) of DAPI in DNA from 100 fs to 10 ns. Time-resolved Stokes shift data suggest that groove-bound DAPI can not sense the counterion dynamics because the ions are displaced by DAPI far from the probe-site. However, when these results are compared to the earlier base-stacked coumarin data, the same ions are found to affect the nanosecond dynamics significantly. This suggests that the ions come close to the probe-site, such that they can affect the dynamics when measured by base-stacked coumarin. These results support previous molecular dynamics (MD) simulation data of groove-bound and base-stacked probes inside DNA.

The structural role of Mg2+ ions in a class I RNA polymerase ribozyme: a molecular simulation study

According to the RNA world hypothesis, self-replicating ribozymes, storing the genetic information and being able to perform catalysis, were the constituents of the first living organisms. In particular, RNA polymerase ribozymes, similar to current proteinaceous enzymatic polymerases, may have been able to promote the synthesis of RNA strands in a primitive world. Polymerase catalysis is usually assisted by Mg(2+) ions, but it is not always trivial to find out experimentally the number of Mg(2+) ions placed in the active site as well as the identity and the number of their coordination ligands. Here, we addressed this issue in an artificial class I ligase ribozyme. On the basis of a recently solved crystal structure, we constructed computational models of reactant and product states of this ribozyme, considering monometallic and bimetallic species. Our models were relaxed by force field based molecular dynamics (MD) simulations and mixed quantum-classical (QM/MM) MD. The structural and dynamical properties of these models were consistent with experimental data and were validated by a comparison with the catalytic sites of proteinaceous DNA and RNA polymerases. Consistently with enzymatic polymerases, our results suggest that class I RNA ligases most probably contain two magnesium ions in the active site and they may, therefore, catalyze the junction of two RNA strands via "a two Mg(2+) ions" mechanism.

The increasing role of QM/MM in drug discovery

Since its first appearance in 1976, the quantum mechanics/molecular mechanics (QM/MM) approach has mostly been used to study the chemical reactions of enzymes, which are frequently the target of drug discovery programs. In principle, a detailed understanding of the enzymatic mechanism should help researchers to design a potent enzyme inhibitor or new drug. However, QM/MM has not yet had a widespread impact on structure-based drug design. This is mostly due to its high computational cost. We expect this to change with the recent and extraordinary increases in computational power, and with the availability of more efficient algorithms for QM/MM calculations. Here, we report on some representative examples of QM/MM studies, including our own research, of pharmaceutically relevant enzymes, such as ribonuclease H and fatty acid amide hydrolase (FAAH). We aim to show how QM/MM has traditionally been used to study enzymatic catalysis. In this regard, we discuss its potential to become a routinely used drug design tool. To support this, we also discuss selected computational studies where QM/MM insights have been helpful in improving the potency of covalent inhibitors of FAAH.

Insight into the phosphodiesterase mechanism from combined QM/MM free energy simulations

Molecular dynamics simulations employing a combined quantum mechanical and molecular mechanical potential have been carried out to elucidate the reaction mechanism of the hydrolysis of a cyclic nucleotide cAMP substrate by phosphodiesterase 4B (PDE4B). PDE4B is a member of the PDE superfamily of enzymes that play crucial roles in cellular signal transduction. We have determined a two-dimensional potential of mean force (PMF) for the coupled phosphoryl bond cleavage and proton transfer through a general acid catalysis mechanism in PDE4B. The results indicate that the ring-opening process takes place through an S(N)2 reaction mechanism, followed by a proton transfer to stabilize the leaving group. The computed free energy of activation for the PDE4B-catalyzed cAMP hydrolysis is about 13 kcal·mol(-1) and an overall reaction free energy is about -17 kcal·mol(-1), both in accord with experimental results. In comparison with the uncatalyzed reaction in water, the enzyme PDE4B provides a strong stabilization of the transition state, lowering the free energy barrier by 14 kcal·mol(-1). We found that the proton transfer from the general acid residue His234 to the O3' oxyanion of the ribosyl leaving group lags behind the nucleophilic attack, resulting in a shallow minimum on the free energy surface. A key contributing factor to transition state stabilization is the elongation of the distance between the divalent metal ions Zn(2+) and Mg(2+) in the active site as the reaction proceeds from the Michaelis complex to the transition state.