Inverse thio effects in the hepatitis delta virus ribozyme reveal that the reaction pathway is controlled by metal ion charge density

RNA nanotechnology on the horizon: Self-assembly, chemical modifications, and functional applications

RNA nanotechnology harnesses the unique chemical and structural properties of RNA to build nanoassemblies and supramolecular structures with dynamic and functional capabilities. This review focuses on design and assembly approaches to building RNA structures, the RNA chemical modifications used to enhance stability and functionality, and modern-day applications in therapeutics, biosensing, and bioimaging.

Simulations predict preferred Mg2+ coordination in a nonenzymatic primer-extension reaction center

The mechanism by which genetic information was copied prior to the evolution of ribozymes is of great interest because of its importance to the origin of life. The most effective known process for the nonenzymatic copying of an RNA template is primer extension by a two-step pathway in which 2-aminoimidazole-activated nucleotides first react with each other to form an imidazolium-bridged intermediate that subsequently reacts with the primer. Reaction kinetics, structure-activity relationships, and X-ray crystallography have provided insight into the overall reaction mechanism, but many puzzles remain. In particular, high concentrations of Mg2+ are required for efficient primer extension, but the mechanism by which Mg2+ accelerates primer extension remains unknown. By analogy with the mechanism of DNA and RNA polymerases, a role for Mg2+ in facilitating the deprotonation of the primer 3'-hydroxyl is often assumed, but no catalytic metal ion is seen in crystal structures of the primer-extension complex. To explore the potential effects of Mg2+ binding in the reaction center, we performed atomistic molecular dynamics simulations of a series of modeled complexes in which a Mg2+ ion was placed in the reaction center with inner-sphere coordination with different sets of functional groups. Our simulations suggest that coordination of a Mg2+ ion with both O3' of the terminal primer nucleotide and the pro-Sp nonbridging oxygen of the reactive phosphate of an imidazolium-bridged dinucleotide would help to pre-organize the structure of the primer/template substrate complex to favor the primer-extension reaction. Our results suggest that the catalytic metal ion may play an important role in overcoming electrostatic repulsion between a deprotonated O3' and the reactive phosphate of the bridged dinucleotide and lead to testable predictions of the mode of Mg2+ binding that is most relevant to catalysis of primer extension.

Geometric deep learning for the prediction of magnesium-binding sites in RNA structures

Magnesium ions (Mg2+) are essential for the folding, functional expression, and structural stability of RNA molecules. However, predicting Mg2+-binding sites in RNA molecules based solely on RNA structures is still challenging. The molecular surface, characterized by a continuous shape with geometric and chemical properties, is important for RNA modelling and carries essential information for understanding the interactions between RNAs and Mg2+ ions. Here, we propose an approach named RNA-magnesium ion surface interaction fingerprinting (RMSIF), a geometric deep learning-based conceptual framework to predict magnesium ion binding sites in RNA structures. To evaluate the performance of RMSIF, we systematically enumerated decoy Mg2+ ions across a full-space grid within the range of 2 to 10 Å from the RNA molecule and made predictions accordingly. Visualization techniques were used to validate the prediction results and calculate success rates. Comparative assessments against state-of-the-art methods like MetalionRNA, MgNet, and Metal3DRNA revealed that RMSIF achieved superior success rates and accuracy in predicting Mg2+-binding sites. Additionally, in terms of the spatial distribution of Mg2+ ions within the RNA structures, a majority were situated in the deep grooves, while a minority occupied the shallow grooves. Collectively, the conceptual framework developed in this study holds promise for advancing insights into drug design, RNA co-transcriptional folding, and structure prediction.

A ribozyme that uses lanthanides as cofactor

To explore how an early, RNA-based life form could have functioned, in vitro selection experiments have been used to develop catalytic RNAs (ribozymes) with relevant functions. We previously identified ribozymes that use the prebiotically plausible energy source cyclic trimetaphosphate (cTmp) to convert their 5'-hydroxyl group to a 5'-triphosphate. While these ribozymes were developed in the presence of Mg2+, we tested here whether lanthanides could also serve as catalytic cofactors because lanthanides are ideal catalytic cations for this reaction. After an in vitro selection in the presence of Yb3+, several active sequences were isolated, and the most active RNA was analyzed in more detail. This ribozyme required lanthanides for activity, with highest activity at a 10:1 molar ratio of cTmp : Yb3+. Only the four heaviest lanthanides gave detectable signals, indicating a high sensitivity of ribozyme catalysis to the lanthanide ion radius. Potassium and Magnesium did not facilitate catalysis alone but they increased the lanthanide-mediated kOBS by at least 100-fold, with both K+ and Mg2+ modulating the ribozyme's secondary structure. Together, these findings show that RNA is able to use the unique properties of lanthanides as catalytic cofactor. The results are discussed in the context of early life forms.

Dissociative Transition State in Hepatitis Delta Virus Ribozyme Catalysis

Ribonucleases and small nucleolytic ribozymes are both able to catalyze RNA strand cleavage through 2'-O-transphosphorylation, provoking the question of whether protein and RNA enzymes facilitate mechanisms that pass through the same or distinct transition states. Here, we report the primary and secondary 18O kinetic isotope effects for hepatitis delta virus ribozyme catalysis that reveal a dissociative, metaphosphate-like transition state in stark contrast to the late, associative transition states observed for reactions catalyzed by specific base, Zn2+ ions, or ribonuclease A. This new information provides evidence for a discrete ribozyme active site design that modulates the RNA cleavage pathway to pass through an altered transition state.

Sensing Metal Ions with Phosphorothioate-Modified DNAzymes

Phosphorothioate (PS) modification refers to replacing one of the nonbridging oxygen atoms in nucleic acids with sulfur. PS modifications can be easily introduced during solid-phase DNA synthesis. It has been extensively used in ribozyme and DNAzyme research to achieve a bioinorganic understanding of metal binding, bioanalytical applications of metal detection, and chemical biology of DNA modification. It allows for the access of new chemistry, not available to natural DNA. Since each PS modification is accompanied by the production of a chiral phosphorus center, a key technical challenge is to separate the two diastereomers called R_p and S_p. In this chapter, we describe our methods of HPLC-based separation followed by ligation to generate a long and fluorescently modified DNAzyme substrate. Subsequently, the use of the modified substrate for activity assay to understand metal binding and for metal ion detection is also described.

Synthesis of Oligoribonucleotides Containing a 2'-Amino-5'-S-phosphorothiolate Linkage

Oligoribonucleotides containing a photocaged 2′-amino-5′-S-phophorothiolate linkage have potential applications as therapeutic agents and biological probes to investigate the RNA structure and function. We envisioned that oligoribonucleotides containing a 2′-amino-5′-S-phosphorothiolate linkage could provide an approach to identify the general base within catalytic RNAs by chemogenetic suppression. To enable preliminary tests of this idea, we developed synthetic approaches to a dinucleotide, trinucleotide, and oligoribonucleotide containing a photocaged 2′-amino-5′-S-phosphorothiolate linkage. We incorporated the photocaged 2′-amino-5′-S-phosphorothiolate linkage into an oligoribonucleotide substrate for the hepatitis delta virus (HDV) ribozyme and investigated the pH dependence of its cleavage following UV irradiation both in the presence and absence of the ribozyme. The substrate exhibited a pH-rate profile characteristic of the modified linkage but reacted slower when bound to the ribozyme. Cleavage inhibition by the HDV ribozyme could reflect a non-productive ground-state interaction with the modified substrate’s nucleophilic 2′-NH₂ or a poor fit of the modified transition state at the ribozyme’s active site.

Single-pass transcription by T7 RNA polymerase

RNA molecules can be conveniently synthesized in vitro by the T7 RNA polymerase (T7 RNAP). In some experiments, such as cotranscriptional biochemical analyses, continuous synthesis of RNA is not desired. Here, we propose a method for a single-pass transcription that yields a single transcript per template DNA molecule using the T7 RNAP system. We hypothesized that stalling the polymerase downstream from the promoter region and subsequent cleavage of the promoter by a restriction enzyme (to prevent promoter binding by another polymerase) would allow synchronized production of a single transcript per template. The single-pass transcription was verified in two different scenarios: a short self-cleaving ribozyme and a long mRNA. The results show that a controlled single-pass transcription using T7 RNAP allows precise measurement of cotranscriptional ribozyme activity, and this approach will facilitate the study of other kinetic events.

In vitro Selection of Chemically Modified DNAzymes

DNAzymes are in vitro selected DNA oligonucleotides with catalytic activities. RNA cleavage is one of the most extensively studied DNAzyme reactions. To expand the chemical functionality of DNA, various chemical modifications have been made during and after selection. In this review, we summarize examples of RNA-cleaving DNAzymes and focus on those modifications introduced during in vitro selection. By incorporating various modified nucleotides via polymerase chain reaction (PCR) or primer extension, a few DNAzymes were obtained that can be specifically activated by metal ions such as Zn2+ and Hg2+. In addition, some modifications were introduced to mimic RNase A that can cleave RNA substrates in the absence of divalent metal ions. In addition, single modifications at the fixed regions of DNA libraries, especially at the cleavage junctions, have been tested, and examples of DNAzymes with phosphorothioate and histidine-glycine modified tertiary amine were successfully obtained specific for Cu2+, Cd2+, Zn2+, and Ni2+. Labeling fluorophore/quencher pair right next to the cleavage junction was also used to obtain signaling DNAzymes for detecting various metal ions and cells. Furthermore, we reviewed work on the cleavage of 2'-5' linked RNA and L-RNA substrates. Finally, applications of these modified DNAzymes as biosensors, RNases, and biochemical probes are briefly described with a few future research opportunities outlined at the end.

Sensitivity of a classic DNAzyme for Pb2+ modulated by cations, anions and buffers

Both cations and anions in salt strongly affect the activity of a classic Pb²⁺ specific DNAzyme, which in turn can affect the sensitivity of related biosensors.

From general base to general acid catalysis in a sodium-specific DNAzyme by a guanine-to-adenine mutation

Recently, a few Na+-specific RNA-cleaving DNAzymes were reported, where nucleobases are likely to play critical roles in catalysis. The NaA43 and NaH1 DNAzymes share the same 16-nt Na+-binding motif, but differ in one or two nucleotides in a small catalytic loop. Nevertheless, they display an opposite pH-dependency, implicating distinct catalytic mechanisms. In this work, rational mutation studies locate a catalytic adenine residue, A22, in NaH1, while previous studies found a guanine (G23) to be important for the catalysis of NaA43. Mutation with pKa-perturbed analogs, such as 2-aminopurine (∼3.8), 2,6-diaminopurine (∼5.6) and hypoxanthine (∼8.7) affected the overall reaction rate. Therefore, we propose that the N1 position of G23 (pKa ∼6.6) in NaA43 functions as a general base, while that of A22 (pKa ∼6.3) in NaH1 as a general acid. Further experiments with base analogs and a phosphorothioate-modified substrate suggest that the exocyclic amine in A22 and both of the non-bridging oxygens at the scissile phosphate are important for catalysis for NaH1. This is an interesting example where single point mutations can change the mechanism of cleavage from general base to general acid, and it can also explain this Na+-dependent DNAzyme scaffold being sensitive to a broad range of metal ions and molecules.

Molecular simulations of the pistol ribozyme: unifying the interpretation of experimental data and establishing functional links with the hammerhead ribozyme

The pistol ribozyme (Psr) is among the most recently discovered RNA enzymes and has been the subject of experiments aimed at elucidating the mechanism. Recent biochemical studies have revealed exciting clues about catalytic interactions in the active site not apparent from available crystallographic data. The present work unifies the interpretation of the existing body of structural and functional data on Psr by providing a dynamical model for the catalytically active state in solution from molecular simulation. Our results suggest that a catalytic Mg2+ ion makes inner-sphere contact with G33:N7 and outer-sphere coordination to the pro-RP of the scissile phosphate, promoting electrostatic stabilization of the dianionic transition state and neutralization of the developing charge of the leaving group through a metal-coordinated water molecule that is made more acidic by a hydrogen bond donated from the 2'OH of P32. This model is consistent with experimental activity-pH and mutagenesis data, including sensitivity to G33(7cG) and phosphorothioate substitution/metal ion rescue. The model suggests several experimentally testable predictions, including the response of cleavage activity to mutations at G42 and P32 positions in the ribozyme, and thio substitutions of the substrate in the presence of different divalent metal ions. Further, the model identifies striking similarities of Psr to the hammerhead ribozyme (HHr), including similar global fold, organization of secondary structure around an active site three-way junction, catalytic metal ion binding mode, and guanine general base. However, the specific binding mode and role of the Mg2+ ion, as well as a conserved 2'-OH in the active site, are interrelated but subtly different between the ribozymes.

Design of highly active double-pseudoknotted ribozymes: a combined computational and experimental study

Design of RNA sequences that adopt functional folds establishes principles of RNA folding and applications in biotechnology. Inverse folding for RNAs, which allows computational design of sequences that adopt specific structures, can be utilized for unveiling RNA functions and developing genetic tools in synthetic biology. Although many algorithms for inverse RNA folding have been developed, the pseudoknot, which plays a key role in folding of ribozymes and riboswitches, is not addressed in most algorithms. For the few algorithms that attempt to predict pseudoknot-containing ribozymes, self-cleavage activity has not been tested. Herein, we design double-pseudoknot HDV ribozymes using an inverse RNA folding algorithm and test their kinetic mechanisms experimentally. More than 90% of the positively designed ribozymes possess self-cleaving activity, whereas more than 70% of negative control ribozymes, which are predicted to fold to the necessary structure but with low fidelity, do not possess it. Kinetic and mutation analyses reveal that these RNAs cleave site-specifically and with the same mechanism as the WT ribozyme. Most ribozymes react just 50- to 80-fold slower than the WT ribozyme, and this rate can be improved to near WT by modification of a junction. Thus, fast-cleaving functional ribozymes with multiple pseudoknots can be designed computationally.

Silicon Nitride: A Bioceramic with a Gift

In the closing decades of the 20th century, silicon nitride (Si₃N₄) was extensively developed for high-temperature gas turbine applications. Technologists attempted to take advantage of its superior thermal and mechanical properties to improve engine reliability and fuel economy. Yet, this promise was never realized in spite of the worldwide research, which was conducted at that time. Notwithstanding this disappointment, its use in medical applications in the early 21st century has been an unexpected gift. While retaining all of its engineered mechanical properties, it is now recognized for its peculiar surface chemistry. When immersed in an aqueous environment, the slow elution of silicon and nitrogen from its surface enhances healing of soft and osseous tissue, inhibits bacterial proliferation, and eradicates viruses. These benefits permit it to be used in a wide array of different disciplines inside and outside of the human body including orthopedics, dentistry, virology, agronomy, and environmental remediation. Given the global public health threat posed by mutating viruses and bacteria, silicon nitride offers a valid and straightforward alternative approach to fighting these pathogens. However, there is a conundrum behind these recent discoveries: How can this unique bioceramic be both friendly to mammalian cells while concurrently lysing invasive pathogens? This unparalleled characteristic can be explained by the pH-dependent kinetics of two ammonia species-NH₄⁺ and NH₃-both of which are leached from the wet Si₃N₄ surface.

Phosphorothioate DNA Mediated Sequence-Insensitive Etching and Ripening of Silver Nanoparticles

Many DNA-functionalized nanomaterials and biosensors have been reported, but most have ignored the influence of DNA on the stability of nanoparticles. We observed that cytosine-rich DNA oligonucleotides can etch silver nanoparticles (AgNPs). In this work, we showed that phosphorothioate (PS)-modified DNA (PS-DNA) can etch AgNPs independently of DNA sequence, suggesting that the thio-modifications are playing the major role in etching. Compared to unmodified DNA (e.g., poly-cytosine DNA), the concentration of required PS DNA decreases sharply, and the reaction rate increases. Furthermore, etching by PS-DNA occurs quite independent of pH, which is also different from unmodified DNA. The PS-DNA mediated etching could also be controlled well by varying DNA length and conformation, and the number and location of PS modifications. With a higher activity of PS-DNA, the process of etching, ripening, and further etching was taken place sequentially. The etching ability is inhibited by forming duplex DNA and thus etching can be used to measure the concentration of complementary DNA.

Evidence That Nucleophile Deprotonation Exceeds Bond Formation in the HDV Ribozyme Transition State

Steric constraints imposed by the active sites of protein and RNA enzymes pose major challenges to the investigation of structure-function relationships within these systems. As a strategy to circumvent such constraints in the HDV ribozyme, we have synthesized phosphoramidites from propanediol derivatives and incorporated them at the 5'-termini of RNA and DNA oligonucleotides to generate a series of novel substrates with nucleophiles perturbed electronically through geminal fluorination. In nonenzymatic, hydroxide-catalyzed intramolecular transphosphorylation of the DNA substrates, pH-rate profiles revealed that fluorine substitution reduces the maximal rate and the kinetic p K_a, consistent with the expected electron-withdrawing effect. In HDV ribozyme reactions, we observed that the RNA substrates undergo transphosphorylation relatively efficiently, suggesting that the conformational constraints imposed by a ribofuranose ring are not strictly required for ribozyme catalysis. In contrast to the nonenzymatic reactions, however, substrate fluorination modestly increases the ribozyme reaction rate, consistent with a mechanism in which (1) the 2'-hydroxyl nucleophile exists predominantly in its neutral, protonated form in the ground state and (2) the 2'-hydroxyl bears some negative charge in the rate-determining step, consistent with a transition state in which the extent of 2'-OH deprotonation exceeds the extent of P-O bond formation.

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.

Prebiotic synthesis of nucleic acids and their building blocks at the atomic level - merging models and mechanisms from advanced computations and experiments

The origin of life on Earth is one of the most fascinating questions of contemporary science. Extensive research in the past decades furnished diverse experimental proposals for the emergence of first informational polymers that could form the basis of the early terrestrial life. Side by side with the experiments, the fast development of modern computational chemistry methods during the last 20 years facilitated the use of in silico modelling tools to complement the experiments. Modern computations can provide unique atomic-level insights into the structural and electronic aspects as well as the energetics of key prebiotic chemical reactions. Many of these insights are not directly obtainable from the experimental techniques and the computations are thus becoming indispensable for proper interpretation of many experiments and for qualified predictions. This review illustrates the synergy between experiment and theory in the origin of life research focusing on the prebiotic synthesis of various nucleic acid building blocks and on the self-assembly of nucleotides leading to the first functional oligonucleotides.

Exploring RNA structure and dynamics through enhanced sampling simulations

RNA function is intimately related to its structural dynamics. Molecular dynamics simulations are useful for exploring biomolecular flexibility but are severely limited by the accessible timescale. Enhanced sampling methods allow this timescale to be effectively extended in order to probe biologically relevant conformational changes and chemical reactions. Here, we review the role of enhanced sampling techniques in the study of RNA systems. We discuss the challenges and promises associated with the application of these methods to force-field validation, exploration of conformational landscapes and ion/ligand-RNA interactions, as well as catalytic pathways. Important technical aspects of these methods, such as the choice of the biased collective variables and the analysis of multi-replica simulations, are examined in detail. Finally, a perspective on the role of these methods in the characterization of RNA dynamics is provided.

Elucidation of Catalytic Strategies of Small Nucleolytic Ribozymes From Comparative Analysis of Active Sites

A number of small, self-cleaving ribozyme classes have been identified including the hammerhead, hairpin, hepatitis delta virus (HDV), Varkud satellite (VS), glmS, twister, hatchet, pistol, and twister sister ribozymes. Within the active sites of these ribozymes, myriad functional groups contribute to catalysis. There has been extensive structure-function analysis of individual ribozymes, but the extent to which catalytic devices are shared across different ribozyme classes is unclear. As such, emergent catalytic principles for ribozymes may await discovery. Identification of conserved catalytic devices can deepen our understanding of RNA catalysis specifically and of enzymic catalysis generally. To probe similarities and differences amongst ribozyme classes, active sites from more than 80 high-resolution crystal structures of self-cleaving ribozymes were compared computationally. We identify commonalities amongst ribozyme classes pertaining to four classic catalytic devices: deprotonation of the 2'OH nucleophile (γ), neutralization of the non-bridging oxygens of the scissile phosphate (β), neutralization of the O5' leaving group (δ), and in-line nucleophilic attack (α). In addition, we uncover conservation of two catalytic devices, each of which centers on the activation of the 2'OH nucleophile by a guanine: one to acidify the 2'OH by hydrogen bond donation to it (γ') and one to acidify the 2'OH by releasing it from non-productive interactions by competitive hydrogen bonding (γ''). Our findings reveal that the amidine functionalities of G, A, and C are especially important for these strategies, and help explain absence of U at ribozyme active sites. The identified γ' and γ'' catalytic strategies help unify the catalytic strategies shared amongst catalytic RNAs and may be important for large ribozymes, as well as protein enzymes that act on nucleic acids.

Quantum mechanical force fields for condensed phase molecular simulations

Molecular simulations are powerful tools for providing atomic-level details into complex chemical and physical processes that occur in the condensed phase. For strongly interacting systems where quantum many-body effects are known to play an important role, density-functional methods are often used to provide the model with the potential energy used to drive dynamics. These methods, however, suffer from two major drawbacks. First, they are often too computationally intensive to practically apply to large systems over long time scales, limiting their scope of application. Second, there remain challenges for these models to obtain the necessary level of accuracy for weak non-bonded interactions to obtain quantitative accuracy for a wide range of condensed phase properties. Quantum mechanical force fields (QMFFs) provide a potential solution to both of these limitations. In this review, we address recent advances in the development of QMFFs for condensed phase simulations. In particular, we examine the development of QMFF models using both approximate and ab initio density-functional models, the treatment of short-ranged non-bonded and long-ranged electrostatic interactions, and stability issues in molecular dynamics calculations. Example calculations are provided for crystalline systems, liquid water, and ionic liquids. We conclude with a perspective for emerging challenges and future research directions.

Two Completely Different Mechanisms for Highly Specific Na+ Recognition by DNAzymes

Our view of the interaction between Na⁺ and nucleic acids was changed by a few recently discovered Na⁺ -specific RNA-cleaving DNAzymes. In addition to nonspecific electrostatic interactions, highly specific recognition is also possible. Herein, two such DNAzymes, named EtNa and Ce13d, are compared to elucidate their mechanisms of Na⁺ binding. Mutation studies indicate that they have different sequence requirements. Phosphorothioate (PS) substitution at the scissile phosphate drops the activity of EtNa 140-fold, and it cannot be rescued by thiophilic Cd²⁺ or Mn²⁺ , whereas the activity of PS-modified Ce13d can be rescued. Na⁺ -dependent activity assays indicate that two Na⁺ ions bind cooperatively in EtNa, and each Na⁺ likely interacts with a nonbridging oxygen atom in the scissile phosphate, whereas Ce13d binds only one Na⁺ ion in a well-defined Na⁺ aptamer, and this Na⁺ ion does not directly interact with the scissile phosphate. Both DNAzymes display a normal pH-rate profile, with a single deprotonation reaction required for catalysis. For EtNa, Na⁺ fails to protect the conserved nucleotides from dimethyl sulfate attack, and no specific Na⁺ binding is detected by 2-aminopurine fluorescence, both of which are different from those observed for Ce13d. This work suggests that EtNa binds Na⁺ mainly through its scissile phosphate without significant involvement of the nucleotides in the enzyme strand, whereas Ce13d has a well-defined aptamer for Na⁺ binding. Therefore, DNA has at least two distinct ways to achieve highly selective Na⁺ binding.

Folding of the silver aptamer in a DNAzyme probed by 2-aminopurine fluorescence

The RNA-cleaving Ag10c DNAzyme was recently isolated via in vitro selection and it can bind two Ag⁺ ions for activity. The Ag10c contains a well-defined Ag⁺ binding aptamer as indicated by DMS footprinting. Since aptamer binding is often accompanied with conformational changes, we herein used 2-aminopurine (2AP) to probe its folding in the presence of Ag⁺. The Ag10c was respectively labeled with 2AP at three different positions, both in the substrate strand and in the enzyme strand, one at a time. Ag⁺-induced folding was observed at the substrate cleavage junction and the A9 position of the enzyme strand, consistent with aptamer binding. The measured K_d at the A9 position was 18 μM Ag⁺ with a Hill coefficient of 2.17, similar to those obtained from the previous cleavage activity based assays. However, labeling a 2AP at the A2 position inhibited the activity and folding. Compared to other metal ions, Ag⁺ has a unique sigmoidal folding profile indicative of multiple silver binding cooperatively. This suggests that Ag⁺ can induce a local folding in the enzyme loop and this folding is important for activity. This study provides important biophysical insights into this new DNAzyme, suggesting the possibility of designing folding-based biosensors for Ag⁺.

Thermal Stability of RNA Structures with Bulky Cations in Mixed Aqueous Solutions

Bulky cations are used to develop nucleic-acid-based technologies for medical and technological applications in which nucleic acids function under nonaqueous conditions. In this study, the thermal stability of RNA structures was measured in the presence of various bulky cations in aqueous mixtures with organic solvents or polymer additives. The stability of oligonucleotide, transfer RNA, and polynucleotide structures was decreased in the presence of salts of tetrabutylammonium and tetrapentylammonium ions, and the stability and salt concentration dependences were dependent on cation sizes. The degree to which stability was dependent on salt concentration was correlated with reciprocals of the dielectric constants of mixed solutions, regardless of interactions between the cosolutes and RNA. Our results show that organic solvents affect the strength of electrostatic interactions between RNA and cations. Analysis of ion binding to RNA indicated greater enhancement of cation binding to RNA single strands than to duplexes in media with low dielectric constants. Furthermore, background bulky ions changed the dependence of RNA duplex stability on the concentration of metal ion salts. These unique properties of large tetraalkylammonium ions are useful for controlling the stability of RNA structures and its sensitivity to metal ion salts.

Understanding in-line probing experiments by modeling cleavage of nonreactive RNA nucleotides

Ribonucleic acid (RNA) is involved in many regulatory and catalytic processes in the cell. The function of any RNA molecule is intimately related with its structure. In-line probing experiments provide valuable structural data sets for a variety of RNAs and are used to characterize conformational changes in riboswitches. However, the structural determinants that lead to differential reactivities in unpaired nucleotides have not been investigated yet. In this work, we used a combination of theoretical approaches, i.e., classical molecular dynamics simulations, multiscale quantum mechanical/molecular mechanical calculations, and enhanced sampling techniques in order to compute and interpret the differential reactivity of individual residues in several RNA motifs, including members of the most important GNRA and UNCG tetraloop families. Simulations on the multinanosecond timescale are required to converge the related free-energy landscapes. The results for uGAAAg and cUUCGg tetraloops and double helices are compared with available data from in-line probing experiments and show that the introduced technique is able to distinguish between nucleotides of the uGAAAg tetraloop based on their structural predispositions toward phosphodiester backbone cleavage. For the cUUCGg tetraloop, more advanced ab initio calculations would be required. This study is the first attempt to computationally classify chemical probing experiments and paves the way for an identification of tertiary structures based on the measured reactivity of nonreactive nucleotides.

The GlcN6P cofactor plays multiple catalytic roles in the glmS ribozyme

RNA enzymes (ribozymes) have remarkably diverse biological roles despite having limited chemical diversity. Protein enzymes enhance their reactivity through recruitment of cofactors; likewise, the naturally occurring glmS ribozyme uses the glucosamine-6-phosphate (GlcN6P) organic cofactor for phosphodiester bond cleavage. Prior structural and biochemical studies have implicated GlcN6P as the general acid. Here we describe new catalytic roles of GlcN6P through experiments and calculations. Large stereospecific normal thio effects and a lack of metal-ion rescue in the holoribozyme indicate that nucleobases and the cofactor play direct chemical roles and align the active site for self-cleavage. Large stereospecific inverse thio effects in the aporibozyme suggest that the GlcN6P cofactor disrupts an inhibitory interaction of the nucleophile. Strong metal-ion rescue in the aporibozyme reveals that this cofactor also provides electrostatic stabilization. Ribozyme organic cofactors thus perform myriad catalytic roles, thereby allowing RNA to compensate for its limited functional diversity.

A Silver-Specific DNAzyme with a New Silver Aptamer and Salt-Promoted Activity

Most RNA-cleaving DNAzymes require a metal ion to interact with the scissile phosphate for activity. Therefore, few unmodified DNAzymes work with thiophilic metals because of their low affinity for phosphate. Recently, an Ag⁺-specific Ag10c DNAzyme was reported via in vitro selection. Herein, Ag10c is characterized to rationalize the role of the strongly thiophilic Ag⁺. Systematic mutation studies indicate that Ag10c is a highly conserved DNAzyme and its Ag⁺ binding is unrelated to C-Ag⁺-C interaction. Its activity is enhanced by increasing Na⁺ concentrations in buffer. At the same metal concentration, activity decreases in the following order: Li⁺ > Na⁺ > K⁺. Ag10c binds one Na⁺ ion and two Ag⁺ ions for catalysis. The pH-rate profile has a slope of ∼1, indicating a single deprotonation step. Phosphorothioate substitution at the scissile phosphate suggests that Na⁺ interacts with the pro-R_p oxygen of the phosphate, and dimethyl sulfate footprinting indicates that the DNAzyme loop is a silver aptamer binding two Ag⁺ ions. Therefore, Ag⁺ exerts its function allosterically, while the scissile phosphate interacts with Na⁺, Li⁺, Na⁺, or Mg²⁺. This work suggests the possibility of isolating thiophilic metal aptamers based on DNAzyme selection, and it also demonstrates a new Ag⁺ aptamer.

Probing fast ribozyme reactions under biological conditions with rapid quench-flow kinetics

Reaction kinetics on the millisecond timescale pervade the protein and RNA fields. To study such reactions, investigators often perturb the system with abiological solution conditions or substrates in order to slow the rate to timescales accessible by hand mixing; however, such perturbations can change the rate-limiting step and obscure key folding and chemical steps that are found under biological conditions. Mechanical methods for collecting data on the millisecond timescale, which allow these perturbations to be avoided, have been developed over the last few decades. These methods are relatively simple and can be conducted on affordable and commercially available instruments. Here, we focus on using the rapid quench-flow technique to study the fast reaction kinetics of RNA enzymes, or ribozymes, which often react on the millisecond timescale under biological conditions. Rapid quench of ribozymes is completely parallel to the familiar hand-mixing approach, including the use of radiolabeled RNAs and fractionation of reactions on polyacrylamide gels. We provide tips on addressing and preventing common problems that can arise with the rapid-quench technique. Guidance is also offered on ensuring the ribozyme is properly folded and fast-reacting. We hope that this article will facilitate the broader use of rapid-quench instrumentation to study fast-reacting ribozymes under biological reaction conditions.

Ambient-Potential Composite Ewald Method for ab Initio Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulation

A new approach for performing Particle Mesh Ewald in ab initio quantum mechanical/molecular mechanical (QM/MM) simulations with extended atomic orbital basis sets is presented. The new approach, the Ambient-Potential Composite Ewald (CEw) method, does not perform the QM/MM interaction with Mulliken charges nor electrostatically fit charges. Instead the nuclei and electron density interact directly with the MM environment, but in a manner that avoids the use of dense Fourier transform grids. By performing the electrostatics with the underlying QM density, the CEw method avoids self-consistent field instabilities that have been encountered with simple charge mapping procedures. Potential of mean force (PMF) profiles of the p-nitrophenyl phosphate dissociation reaction in explicit solvent are computed from PBE0/6-31G* QM/MM molecular dynamics simulations with various electrostatic protocols. The CEw profiles are shown to be stable with respect to real-space Ewald cutoff, whereas the PMFs computed from truncated and switched electrostatics produce artifacts. PBE0/6-311G**, AM1/d-PhoT, and DFTB2 QM/MM simulations are performed to generate two-dimensional PMF profiles of the phosphoryl transesterification reactions with ethoxide and phenoxide leaving groups. The semiempirical models incorrectly produce a concerted ethoxide mechanism, whereas PBE0 correctly produces a stepwise mechanism. The ab initio reaction barriers agree more closely to experiment than the semiempirical models. The failure of Mulliken-charge QM/MM-Ewald is analyzed.

VR-SCOSMO: A smooth conductor-like screening model with charge-dependent radii for modeling chemical reactions

To better represent the solvation effects observed along reaction pathways, and of ionic species in general, a charge-dependent variable-radii smooth conductor-like screening model (VR-SCOSMO) is developed. This model is implemented and parameterized with a third order density-functional tight binding quantum model, DFTB3/3OB-OPhyd, a quantum method which was developed for organic and biological compounds, utilizing a specific parameterization for phosphate hydrolysis reactions. Unlike most other applications with the DFTB3/3OB model, an auxiliary set of atomic multipoles is constructed from the underlying DFTB3 density matrix which is used to interact the solute with the solvent response surface. The resulting method is variational, produces smooth energies, and has analytic gradients. As a baseline, a conventional SCOSMO model with fixed radii is also parameterized. The SCOSMO and VR-SCOSMO models shown have comparable accuracy in reproducing neutral-molecule absolute solvation free energies; however, the VR-SCOSMO model is shown to reduce the mean unsigned errors (MUEs) of ionic compounds by half (about 2-3 kcal/mol). The VR-SCOSMO model presents similar accuracy as a charge-dependent Poisson-Boltzmann model introduced by Hou et al. [J. Chem. Theory Comput. 6, 2303 (2010)]. VR-SCOSMO is then used to examine the hydrolysis of trimethylphosphate and seven other phosphoryl transesterification reactions with different leaving groups. Two-dimensional energy landscapes are constructed for these reactions and calculated barriers are compared to those obtained from ab initio polarizable continuum calculations and experiment. Results of the VR-SCOSMO model are in good agreement in both cases, capturing the rate-limiting reaction barrier and the nature of the transition state.

In Vitro Selection of a DNAzyme Cooperatively Binding Two Lanthanide Ions for RNA Cleavage

Trivalent lanthanide ions (Ln(3+)) were recently employed to select RNA-cleaving DNAzymes, and three new DNAzymes have been reported so far. In this work, dysprosium (Dy(3+)) was used with a library containing 50 random nucleotides. After six rounds of in vitro selection, a new DNAzyme named Dy10a was obtained and characterized. Dy10a has a bulged hairpin structure cleaving a RNA/DNA chimeric substrate. Dy10a is highly active in the presence of the five Ln(3+) ions in the middle of the lanthanide series (Sm(3+), Eu(3+), Gd(3+), Tb(3+), and Dy(3+)), while its activity descends on the two sides. The cleavage rate reaches 0.6 min(-1) at pH 6 with just 200 nM Sm(3+), which is the fastest among all known Ln(3+)-dependent enzymes. Dy10a binds two Ln(3+) ions cooperatively. When a phosphorothioate (PS) modification is introduced at the cleavage junction, the activity decreases by >2500-fold for both the Rp and Sp diastereomers, and thiophilic Cd(2+) cannot rescue the activity. The pH-rate profile has a slope of 0.37 between pH 4.2 and 5.2, and the slope was even lower at higher pH. On the basis of these data, a model of metal binding is proposed. Finally, a catalytic beacon sensor that can detect Ho(3+) down to 1.7 nM is constructed.

An Ultrasensitive Light-up Cu(2+) Biosensor Using a New DNAzyme Cleaving a Phosphorothioate-Modified Substrate

Cu(2+) is a very important metal ion in biology, environmental science, and industry. Developing biosensors for Cu(2+) is a key topic in analytical chemistry. DNAzyme-based sensors are highly attractive for their excellent sensitivity, stability, and programmability. In the past decade, a few Cu(2+) biosensors were reported using DNAzymes with DNA cleavage or DNA ligation activity. However, they require unstable ascorbate or imidazole activation. So far, no RNA-cleaving DNAzymes specific for Cu(2+) are known. In this work, a phosphorothioate (PS) RNA-containing library was used for in vitro selection, and a few new Cu(2+)-specific RNA-cleaving DNAzymes were isolated. Among them, a DNAzyme named PSCu10 was studied further. It has only eight nucleotides in the enzyme loop with a cleavage rate of 0.1 min(-1) in the presence of 1 μM Cu(2+) at pH 6.0 (its optimal pH). Between the two diastereomers of the PS RNA chiral center, the R(p) isomer is 37 times more active than the S(p) one. Among the other divalent metal ions, only Hg(2+) can cleave the substrate due to its extremely high thiophilicity. A catalytic beacon sensor was designed with a detection limit of 1.6 nM Cu(2+) and extremely high selectivity. PSCu10 is specific for Cu(2+), and it has no cleavage in the presence of ascorbate, which reduces Cu(2+) to Cu(+).

A Two-Metal-Ion-Mediated Conformational Switching Pathway for HDV Ribozyme Activation

RNA enzymes serve as a potentially powerful platform from which to design catalysts and engineer new biotechnology. A fundamental understanding of these systems provides insight to guide design. The hepatitis delta virus ribozyme (HDVr) is a small, self-cleaving RNA motif widely distributed in nature, that has served as a paradigm for understanding basic principles of RNA catalysis. Nevertheless, questions remain regarding the precise roles of divalent metal ions and key nucleotides in catalysis. In an effort to establish a reaction mechanism model consistent with available experimental data, we utilize molecular dynamics simulations to explore different conformations and metal ion binding modes along the HDVr reaction path. Building upon recent crystallographic data, our results provide a dynamic model of the HDVr reaction mechanism involving a conformational switch between multiple non-canonical G25:U20 base pair conformations in the active site. These local nucleobase dynamics play an important role in catalysis by modulating the metal binding environments of two Mg²⁺ ions that support catalysis at different steps of the reaction pathway. The first ion plays a structural role by inducing a base pair flip necessary to obtain the catalytic fold in which C75 moves towards to the scissile phosphate in the active site. Ejection of this ion then permits a second ion to bind elsewhere in the active site and facilitate nucleophile activation. The simulations collectively describe a mechanistic scenario that is consistent with currently available experimental data from crystallography, phosphorothioate substitutions, and chemical probing studies. Avenues for further experimental verification are suggested.

Sodium and Potassium Interactions with Nucleic Acids

Metal ions are essential cofactors for the structure and functions of nucleic acids. Yet, the early discovery in the 70s of the crucial role of Mg(2+) in stabilizing tRNA structures has occulted for a long time the importance of monovalent cations. Renewed interest in these ions was brought in the late 90s by the discovery of specific potassium metal ions in the core of a group I intron. Their importance in nucleic acid folding and catalytic activity is now well established. However, detection of K(+) and Na(+) ions is notoriously problematic and the question about their specificity is recurrent. Here we review the different methods that can be used to detect K(+) and Na(+) ions in nucleic acid structures such as X-ray crystallography, nuclear magnetic resonance or molecular dynamics simulations. We also discuss specific versus non-specific binding to different structures through various examples.

Chemistry and Biology of Self-Cleaving Ribozymes

Self-cleaving ribozymes were discovered 30 years ago, but their biological distribution and catalytic mechanisms are only beginning to be defined. Each ribozyme family is defined by a distinct structure, with unique active sites accelerating the same transesterification reaction across the families. Biochemical studies show that general acid-base catalysis is the most common mechanism of self-cleavage, but metal ions and metabolites can be used as cofactors. Ribozymes have been discovered in highly diverse genomic contexts throughout nature, from viroids to vertebrates. Their biological roles include self-scission during rolling-circle replication of RNA genomes, co-transcriptional processing of retrotransposons, and metabolite-dependent gene expression regulation in bacteria. Other examples, including highly conserved mammalian ribozymes, suggest that many new biological roles are yet to be discovered.

Assessment of metal-assisted nucleophile activation in the hepatitis delta virus ribozyme from molecular simulation and 3D-RISM

The hepatitis delta virus ribozyme is an efficient catalyst of RNA 2'-O-transphosphorylation and has emerged as a key experimental system for identifying and characterizing fundamental features of RNA catalysis. Recent structural and biochemical data have led to a proposed mechanistic model whereby an active site Mg(2+) ion facilitates deprotonation of the O2' nucleophile, and a protonated cytosine residue (C75) acts as an acid to donate a proton to the O5' leaving group as noted in a previous study. This model assumes that the active site Mg(2+) ion forms an inner-sphere coordination with the O2' nucleophile and a nonbridging oxygen of the scissile phosphate. These contacts, however, are not fully resolved in the crystal structure, and biochemical data are not able to unambiguously exclude other mechanistic models. In order to explore the feasibility of this model, we exhaustively mapped the free energy surfaces with different active site ion occupancies via quantum mechanical/molecular mechanical (QM/MM) simulations. We further incorporate a three-dimensional reference interaction site model for the solvated ion atmosphere that allows these calculations to consider not only the rate associated with the chemical steps, but also the probability of observing the system in the presumed active state with the Mg(2+) ion bound. The QM/MM results predict that a pathway involving metal-assisted nucleophile activation is feasible based on the rate-controlling transition state barrier departing from the presumed metal-bound active state. However, QM/MM results for a similar pathway in the absence of Mg(2+) are not consistent with experimental data, suggesting that a structural model in which the crystallographically determined Mg(2+) is simply replaced with Na(+) is likely incorrect. It should be emphasized, however, that these results hinge upon the assumption of the validity of the presumed Mg(2+)-bound starting state, which has not yet been definitively verified experimentally, nor explored in depth computationally. Thus, further experimental and theoretical study is needed such that a consensus view of the catalytic mechanism emerges.