Charged nucleobases and their potential for RNA catalysis

Insights into the A-C Mismatch Conformational Ensemble in Duplex DNA and its Role in Genetic Processes through a Structure-based Review

Knowing the conformational ensembles formed by mismatches is crucial for understanding how they are generated and repaired and how they contribute to genomic instability. Here, we review structural and energetic studies of the A-C mismatch in duplex DNA and use the information to identify critical conformational states in its ensemble and their significance in genetic processes. In the 1970s, Topal and Fresco proposed the A-C wobble stabilized by two hydrogen bonds, one requiring protonation of adenine-N1. Subsequent NMR and X-ray crystallography studies showed that the protonated A-C wobble was in dynamic equilibrium with a neutral inverted wobble. The mismatch was shown to destabilize duplex DNA in a sequence- and pH-dependent manner by 2.4-3.8 kcal/mol and to have an apparent pKa ranging between 7.2 and 7.7. The A-C mismatch conformational repertoire expanded as structures were determined for damaged and protein-bound DNA. These structures included Watson-Crick-like conformations forming through tautomerization of the bases that drive replication errors, the reverse wobble forming through rotation of the entire nucleotide proposed to increase the fidelity of DNA replication, and the Hoogsteen base-pair forming through the flipping of the adenine base which explained the unusual specificity of DNA polymerases that bypass DNA damage. Thus, the A-C mismatch ensemble encompasses various conformational states that can be selectively stabilized in response to environmental changes such as pH shifts, intermolecular interactions, and chemical modifications, and these adaptations facilitate critical biological processes. This review also highlights the utility of existing 3D structures to build ensemble models for nucleic acid motifs.

Radiation-Induced Transfer of Charge, Atoms, and Energy within Isolated Biomolecular Systems

In biological tissues, ionizing radiation interacts with a variety of molecules and the consequences include cell killing and the modification of mechanical properties. Applications of biological radiation action are for instance radiotherapy, sterilization, or the tailoring of biomaterial properties. During the first femtoseconds to milliseconds after the initial radiation action, biomolecular systems typically respond by transfer of charge, atoms, or energy. In the condensed phase, it is usually very difficult to distinguish direct effects from indirect effects. A straightforward solution for this problem is the use of gas-phase techniques, for instance from the field of mass spectrometry. In this review, we survey mainly experimental but also theoretical work, focusing on radiation-induced intra- and inter-molecular transfer of charge, atoms, and energy within biomolecular systems in the gas phase. Building blocks of DNA, proteins, and saccharides, but also antibiotics are considered. The emergence of general processes as well as their timescales and mechanisms are highlighted.

In vivo-like nearest neighbor parameters improve prediction of fractional RNA base-pairing in cells

We conducted a thermodynamic analysis of RNA stability in Eco80 artificial cytoplasm, which mimics in vivo conditions, and compared it to transcriptome-wide probing of mRNA. Eco80 contains 80% of Escherichia coli metabolites, with biological concentrations of metal ions, including 2 mM free Mg2+ and 29 mM metabolite-chelated Mg2+. Fluorescence-detected binding isotherms (FDBI) were used to conduct a thermodynamic analysis of 24 RNA helices and found that these helices, which have an average stability of -12.3 kcal/mol, are less stable by ΔΔGo37 ∼1 kcal/mol. The FDBI data was used to determine a set of Watson-Crick free energy nearest neighbor parameters (NNPs), which revealed that Eco80 reduces the stability of three NNPs. This information was used to adjust the NN model using the RNAstructure package. The in vivo-like adjustments have minimal effects on the prediction of RNA secondary structures determined in vitro and in silico, but markedly improve prediction of fractional RNA base pairing in E. coli, as benchmarked with our in vivo DMS and EDC RNA chemical probing data. In summary, our thermodynamic and chemical probing analyses of RNA helices indicate that RNA secondary structures are less stable in cells than in artificially stable in vitro buffer conditions.

An RNA-Cleaving DNAzyme That Requires an Organic Solvent to Function

Engineering functional nucleic acids that are active under unusual conditions will not only reveal their hidden abilities but also lay the groundwork for pursuing them for unique applications. Although many DNAzymes have been derived to catalyze diverse chemical reactions in aqueous solutions, no prior study has been set up to purposely derive DNAzymes that require an organic solvent to function. Herein, we utilized in vitro selection to isolate RNA-cleaving DNAzymes from a random-sequence DNA pool that were "compelled" to accept 35 % dimethyl sulfoxide (DMSO) as a cosolvent, via counter selection in a purely aqueous solution followed by positive selection in the same solution containing 35 % DMSO. This experiment led to the discovery of a new DNAzyme that requires 35 % DMSO for its catalytic activity and exhibits drastically reduced activity without DMSO. This DNAzyme also requires divalent metal ions for catalysis, and its activity is enhanced by monovalent ions. A minimized, more efficient DNAzyme was also derived. This work demonstrates that highly functional, organic solvent-dependent DNAzymes can be isolated from random-sequence DNA libraries via forced in vitro selection, thus expanding the capability and potential utility of catalytic DNA.

Characterizing Protonation-Coupled Conformational Ensembles in RNA via pH-Differential Mutational Profiling with DMS Probing

RNA molecules undergo conformational transitions in response to cellular and environmental stimuli. Site-specific protonation, a fundamental chemical property, can alter the conformational landscape of RNA to regulate their functions. However, characterizing protonation-coupled RNA conformational ensembles on a large scale remains challenging. Here, we present pH-differential mutational profiling (PD-MaP) with dimethyl sulfate probing for high-throughput detection of protonation-coupled conformational ensembles in RNA. We demonstrated this approach on microRNA-21 precursor (pre-miR-21) and recapitulated a previously discovered A+-G-coupled conformational ensemble. Additionally, we identified a secondary protonation event involving an A+-C mismatch. We validated the occurrence of both protonation-coupled ensembles in pre-miR-21 using NMR relaxation dispersion spectroscopy. Furthermore, the application of PD-MaP on a library of well-annotated human primary microRNAs uncovered widespread protonation-coupled conformational ensembles, suggesting their potentially broad functions in biology.

In vitro evolution of ribonucleases from expanded genetic alphabets

The ability of nucleic acids to catalyze reactions (as well as store and transmit information) is important for both basic and applied science, the first in the context of molecular evolution and the origin of life and the second for biomedical applications. However, the catalytic power of standard nucleic acids (NAs) assembled from just four nucleotide building blocks is limited when compared with that of proteins. Here, we assess the evolutionary potential of libraries of nucleic acids with six nucleotide building blocks as reservoirs for catalysis. We compare the outcomes of in vitro selection experiments toward RNA-cleavage activity of two nucleic acid libraries: one built from the standard four independently replicable nucleotides and the other from six, with the two added nucleotides coming from an artificially expanded genetic information system (AEGIS). Results from comparative experiments suggest that DNA libraries with increased chemical diversity, higher information density, and larger searchable sequence spaces are one order of magnitude richer reservoirs of molecules that catalyze the cleavage of a phosphodiester bond in RNA than DNA libraries built from a standard four-nucleotide alphabet. Evolved AEGISzymes with nitro-carrying nucleobase Z appear to exploit a general acid–base catalytic mechanism to cleave that bond, analogous to the mechanism of the ribonuclease A family of protein enzymes and heavily modified DNAzymes. The AEGISzyme described here represents a new type of catalysts evolved from libraries built from expanded genetic alphabets.

Scaling Catalytic Contributions of Small Self-Cleaving Ribozymes

Nucleolytic ribozymes utilize general acid-base catalysis to perform phosphodiester cleavage. In most ribozyme classes, a conserved active site guanosine is positioned to act as general base, thereby activating the 2'-OH group to attack the scissile phosphate (γ-catalysis). Here, we present an atomic mutagenesis study for the pistol ribozyme class. Strikingly, "general base knockout" by replacement of the guanine N1 atom by carbon results in only 2.7-fold decreased rate. Therefore, the common view that γ-catalysis critically depends on the N1 moiety becomes challenged. For pistol ribozymes we found that γ-catalysis is subordinate in overall catalysis, made up by two other catalytic factors (α and δ). Our approach allows scaling of the different catalytic contributions (α, β, γ, δ) with unprecedented precision and paves the way for a thorough mechanistic understanding of nucleolytic ribozymes with active site guanines.

Scaling Catalytic Contributions of Small Self-Cleaving Ribozymes

Nucleolytic ribozymes utilize general acid-base catalysis to perform phosphodiester cleavage. In most ribozyme classes, a conserved active site guanosine is positioned to act as general base, thereby activating the 2'-OH group to attack the scissile phosphate (γ-catalysis). Here, we present an atomic mutagenesis study for the pistol ribozyme class. Strikingly, "general base knockout" by replacement of the guanine N1 atom by carbon results in only 2.7-fold decreased rate. Therefore, the common view that γ-catalysis critically depends on the N1 moiety becomes challenged. For pistol ribozymes we found that γ-catalysis is subordinate in overall catalysis, made up by two other catalytic factors (α and δ). Our approach allows scaling of the different catalytic contributions (α, β, γ, δ) with unprecedented precision and paves the way for a thorough mechanistic understanding of nucleolytic ribozymes with active site guanines.

Ozonation of organic compounds in water and wastewater: A critical review

Ozonation has been applied in water treatment for more than a century, first for disinfection, later for oxidation of inorganic and organic pollutants. In recent years, ozone has been increasingly applied for enhanced municipal wastewater treatment for ecosystem protection and for potable water reuse. These applications triggered significant research efforts on the abatement efficiency of organic contaminants and the ensuing formation of transformation products. This endeavor was accompanied by developments in analytical and computational chemistry, which allowed to improve the mechanistic understanding of ozone reactions. This critical review assesses the challenges of ozonation of impaired water qualities such as wastewaters and provides an up-to-date compilation of the recent kinetic and mechanistic findings of ozone reactions with dissolved organic matter, various functional groups (olefins, aromatic compounds, heterocyclic compounds, aliphatic nitrogen-containing compounds, sulfur-containing compounds, hydrocarbons, carbanions, β-diketones) and antibiotic resistance genes.

Structure and mechanism of the methyltransferase ribozyme MTR1

RNA-catalysed RNA methylation was recently shown to be part of the catalytic repertoire of ribozymes. The methyltransferase ribozyme MTR1 catalyses the site-specific synthesis of 1-methyladenosine (m¹A) in RNA, using O⁶-methylguanine (m⁶G) as methyl group donor. Here we report the crystal structure of MTR1 at a resolution of 2.8 Å, which reveals a guanine binding site reminiscent of natural guanine riboswitches. The structure represents the postcatalytic state of a split ribozyme in complex with the m¹A-containing RNA product and the demethylated cofactor guanine. The structural data suggest the mechanistic involvement of a protonated cytidine in the methyl transfer reaction. A synergistic effect of two 2'-O-methylated ribose residues in the active site results in accelerated methyl group transfer. Supported by these results, it seems plausible that modified nucleotides may have enhanced early RNA catalysis and that metabolite-binding riboswitches may resemble inactivated ribozymes that have lost their catalytic activity during evolution.

Not making the cut: Techniques to prevent RNA cleavage in structural studies of RNase-RNA complexes

RNases are varied in the RNA structures and sequences they target for cleavage and are an important type of enzyme in cells. Despite the numerous examples of RNases known, and of those with determined three-dimensional structures, relatively few examples exist with the RNase bound to intact cognate RNA substrate prior to cleavage. To better understand RNase structure and sequence specificity for RNA targets, in vitro methods used to assemble these enzyme complexes trapped in a pre-cleaved state have been developed for a number of different RNases. We have surveyed the Protein Data Bank for such structures and in this review detail methodologies that have successfully been used and relate them to the corresponding structures. We also offer ideas and suggestions for future method development. Many strategies within this review can be used in combination with X-ray crystallography, as well as cryo-EM, and other structure-solving techniques. Our hope is that this review will be used as a guide to resolve future yet-to-be-determined RNase-substrate complex structures.

Understanding the adsorption of plasmid DNA and RNA molecules onto arginine-agarose chromatographic resin

Background

The production of nucleic acids (plasmid DNA or mRNA) in response to the development of new advanced vaccine platforms has greatly increased recently, mostly resulting from the pandemic situation. Due to the intended pharmaceutical use, nucleic acids preparations must fulfill all the required specifications in terms of purity and quality. Chromatography is a standard operation used to isolate these molecules from impurities, playing a central role in the manufacturing processes. However, the mechanism of nucleic acid adsorption in chromatographic resins is poorly understood, often leading to low adsorption capacities and a lack of specificity.

Methods and results

Here we investigated the adsorption of plasmid DNA and RNA molecules onto arginine-agarose, a resin with potential for large-scale application. Equilibrium batch studies were performed through pre-purified samples, using arginine-based ligands by varying the adsorption conditions in the pH value range from 6.0 to 9.0. Langmuir and Freundlich isotherm models were used to describe the adsorption equilibrium. The best fit for both nucleic acids was achieved using the Freundlich model. The correct choice of pH showed critical for controlling the efficacy of arginine-nucleic acid interaction, due to its influence on the nucleic acid structures. This type of analysis is necessary for the improvement of the selectivity and binding capacities of the resins used for plasmid DNA or mRNA purification.

Conclusions

The results presented here indicate that adsorption conditions can be tuned to enhance separation between pDNA and RNA, an important feature in the purification of nucleic acids for vaccine production.

A Combined Experimental and Computational Study of Halogen and Hydrogen Bonding in Molecular Salts of 5-Bromocytosine

Although natural or artificial modified pyrimidine nucleobases represent important molecules with valuable properties as constituents of DNA and RNA, no systematic analyses of the structural aspects of bromo derivatives of cytosine have appeared so far in the literature. In view of the biochemical and pharmaceutical relevance of these compounds, six different crystals containing proton-transfer derivatives of 5-bromocytosine are prepared and analyzed in the solid-state by single crystal X-ray diffraction. All six compounds are organic salts, with proton transfer occurring to the N_imino atom of the pyridine ring. Experimental results are then complemented with Hirshfeld surface analysis to quantitively evaluate the contribution of different intermolecular interactions in the crystal packing. Furthermore, theoretical calculations, based on different arrangements of molecules extracted from the crystal structure determinations, are carried out to analyze the formation mechanism of halogen bonds (XBs) in these compounds and provide insights into the nature and strength of the observed interactions. The results show that the supramolecular architectures of the six molecular salts involve extensive classical intermolecular hydrogen bonds. However, in all but one proton-transfer adducts, weak to moderate XBs are revealed by C-Br^…O short contacts between the bromine atom in the fifth position, which acts as XB donor (electron acceptor). Moreover, the lone pair electrons of the oxygen atom of adjacent pyrimidine nucleobases and/or counterions or water molecules, which acts as XB acceptor (electron donor).

Prebiotically-relevant low polyion multivalency can improve functionality of membraneless compartments

Multivalent polyions can undergo complex coacervation, producing membraneless compartments that accumulate ribozymes and enhance catalysis, and offering a mechanism for functional prebiotic compartmentalization in the origins of life. Here, we evaluate the impact of lower, more prebiotically-relevant, polyion multivalency on the functional performance of coacervates as compartments. Positively and negatively charged homopeptides with 1–100 residues and adenosine mono-, di-, and triphosphate nucleotides are used as model polyions. Polycation/polyanion pairs are tested for coacervation, and resulting membraneless compartments are analyzed for salt resistance, ability to provide a distinct internal microenvironment (apparent local pH, RNA partitioning), and effect on RNA structure formation. We find that coacervates formed by phase separation of the shorter polyions more effectively generated distinct pH microenvironments, accumulated RNA, and preserved duplexes than those formed by longer polyions. Hence, coacervates formed by reduced multivalency polyions are not only viable as functional compartments for prebiotic chemistries, they can outperform higher molecular weight analogues.

Short cationic peptides and nucleotides can form complex coacervates, but the influence of reduced multivalency on coacervate functionality was not investigated. Here, the authors report that coacervates formed from short polyions generate distinct pH microenvironments, accumulate RNA and preserve nucleic acid duplexes more efficiently than their longer counterparts.

Microhydration of protonated biomolecular building blocks: protonated pyrimidine

Protonation and hydration of biomolecules govern their structure, conformation, and function. Herein, we explore the microhydration structure in mass-selected protonated pyrimidine-water clusters (H⁺Pym-W_n, n = 1-4) by a combination of infrared photodissociation spectroscopy (IRPD) between 2450 and 3900 cm^-1 and density functional theory (DFT) calculations at the dispersion-corrected B3LYP-D3/aug-cc-pVTZ level. We further present the IR spectrum of H⁺Pym-N₂ to evaluate the effect of solvent polarity on the intrinsic molecular parameters of H⁺Pym. Our combined spectroscopic and computational approach unequivocally shows that protonation of Pym occurs at one of the two equivalent basic ring N atoms and that the ligands in H⁺Pym-L (L = N₂ or W) preferentially form linear H-bonds to the resulting acidic NH group. Successive addition of water ligands results in the formation of a H-bonded solvent network which increasingly weakens the NH group. Despite substantial activation of the N-H bond upon microhydration, no intracluster proton transfer occurs up to n = 4 because of the balance of relative proton affinities of Pym and W_n and the involved solvation energies. Comparison to neutral Pym-W_n clusters reveals the drastic effects of protonation on microhydration with respect to both structure and interaction strength.

An IRMPD spectroscopic and computational study of protonated guanine-containing mismatched base pairs in the gas phase

Infrared multiple photon dissociation spectroscopy has been used to probe the structures of the three protonated base-pair mismatches containing 9-ethylguanine (9eG) in the gas phase. Some of these protonated base-pairs have been identified in RNA.

The effect of adenine protonation on RNA phosphodiester backbone bond cleavage elucidated by deaza-nucleobase modifications and mass spectrometry

The catalytic strategies of small self-cleaving ribozymes often involve interactions between nucleobases and the ribonucleic acid (RNA) backbone. Here we show that multiply protonated, gaseous RNA has an intrinsic preference for the formation of ionic hydrogen bonds between adenine protonated at N3 and the phosphodiester backbone moiety on its 5'-side that facilitates preferential phosphodiester backbone bond cleavage upon vibrational excitation by low-energy collisionally activated dissociation. Removal of the basic N3 site by deaza-modification of adenine was found to abrogate preferential phosphodiester backbone bond cleavage. No such effects were observed for N1 or N7 of adenine. Importantly, we found that the pH of the solution used for generation of the multiply protonated, gaseous RNA ions by electrospray ionization affects phosphodiester backbone bond cleavage next to adenine, which implies that the protonation patterns in solution are at least in part preserved during and after transfer into the gas phase. Our study suggests that interactions between protonated adenine and phosphodiester moieties of RNA may play a more important mechanistic role in biological processes than considered until now.

Evidence for Hidden Involvement of N3-Protonated Guanine in RNA Structure and Function

Charged nucleobases have been found to occur in several known RNA molecules and are considered essential for their structure and function. The mechanism of their involvement is however not yet fully understood. Revelation of the role of N7-protonated guanine, in modulating the geometry and stability of noncanonical base pairs formed through its unprotonated edges [Watson-Crick (WC) and sugar], has triggered the need to evaluate the feasibility of similar roles of other protonated nucleobases [Halder et al., Phys Chem Chem Phys, 2015, 17, 26249]. In this context, N3 protonation of guanine makes an interesting case as its influence on the charge distribution of the WC edge is similar to that of N7 protonation, though its thermodynamic cost of protonation is significantly higher. In this work, we have carried out structural bioinformatics analyses and quantum mechanics-based calculations to show that N3 protonation of guanine may take place in a cellular environment, at least in the G:C W:W Trans and G:G W:H Cis base pairs. Our results provide a reasonable starting point for future investigations in order to address the larger mechanistic question.

In vivo RNA structural probing of uracil and guanine base-pairing by 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide (EDC)

Many biological functions performed by RNAs arise from their in vivo structures. The structure of the same RNA can differ in vitro and in vivo owing in part to the influence of molecules ranging from protons to secondary metabolites to proteins. Chemical reagents that modify the Watson-Crick (WC) face of unprotected RNA bases report on the absence of base-pairing and so are of value to determining structures adopted by RNAs. Reagents have thus been sought that can report on the native RNA structures that prevail in living cells. Dimethyl sulfate (DMS) and glyoxal penetrate cell membranes and inform on RNA secondary structure in vivo through modification of adenine (A), cytosine (C), and guanine (G) bases. Uracil (U) bases, however, have thus far eluded characterization in vivo. Herein, we show that the water-soluble carbodiimide 1-ethyl-3-(3-dimethylaminopropyl)carbodiimide (EDC) is capable of modifying the WC face of U and G in vivo, favoring the former nucleobase by a factor of ∼1.5, and doing so in the eukaryote rice, as well as in the Gram-negative bacterium Escherichia coli While both EDC and glyoxal target Gs, EDC reacts with Gs in their typical neutral state, while glyoxal requires Gs to populate the rare anionic state. EDC may thus be more generally useful; however, comparison of the reactivity of EDC and glyoxal may allow the identification of Gs with perturbed pK_as in vivo and genome-wide. Overall, use of EDC with DMS allows in vivo probing of the base-pairing status of all four RNA bases.

Temperature-Induced Replacement of Phosphate Proton with Metal Ion Captured in Neutron Structures of A-DNA

Nucleic acids can fold into well-defined 3D structures that help determine their function. Knowing precise nucleic acid structures can also be used for the design of nucleic acid-based therapeutics. However, locations of hydrogen atoms, which are key players of nucleic acid function, are normally not determined with X-ray crystallography. Accurate determination of hydrogen atom positions can provide indispensable information on protonation states, hydrogen bonding, and water architecture in nucleic acids. Here, we used neutron crystallography in combination with X-ray diffraction to obtain joint X-ray/neutron structures at both room and cryo temperatures of a self-complementary A-DNA oligonucleotide d[GTGG(CSe)CAC]2 containing 2'-SeCH3 modification on Cyt5 (CSe) at pH 5.6. We directly observed protonation of a backbone phosphate oxygen of Ade7 at room temperature. The proton is replaced with hydrated Mg2+ upon cooling the crystal to 100 K, indicating that metal binding is favored at low temperature, whereas proton binding is dominant at room temperature.

Energetic, Topological and Electric Field Analyses of Cation-Cation Nucleic Acid Interactions in Watson-Crick Disposition

A theoretical study of the effect of the diprotonation on the nucleic acid bases (A : U, A : T and G : C) in Watson-Crick conformation has been carried out by means of DFT computational methods in vacuum. In addition, the corresponding neutral and monoprotonated binary complexes have been considered. Most of the diprotonated species studied are stable, even though the binding energy is positive due to the overall repulsive electrostatic term. Local electrostatic attractive forces in the regions of hydrogen bonds (HBs) are responsible for equilibrium geometries, as shown by the electric field lines connecting the electrophilic and nucleophilic sites involved in the HB interactions. Secondary electrostatic effects also affect the assembling of the nucleic acid complexes in either neutral or cationic form. In particular, the electric field lines flowing from electrophilic sites in one base to nucleophilic sites in the other reinforce the linking between them. Hence, when the nucleophilic site concerns the free lone pair of the heteroatom involved in the HB interaction as acceptor, the HB distance shortens. However, if the free lone pair of the HB acceptor interacts with an electrophilic site in the same molecule, the HB distance elongates, weakening the HB interaction. The topological analysis of the electron density distribution in HB regions indicates that neutral, monoprotonated and diprotonated complexes show no differences in the nature of their HB's.

Adenine protonation enables cyclic-di-GMP binding to cyclic-GAMP sensing riboswitches

In certain structural or functional contexts, RNA structures can contain protonated nucleotides. However, a direct role for stably protonated nucleotides in ligand binding and ligand recognition has not yet been demonstrated unambiguously. Previous X-ray structures of c-GAMP binding riboswitch aptamer domains in complex with their near-cognate ligand c-di-GMP suggest that an adenine of the riboswitch either forms two hydrogen bonds to a G nucleotide of the ligand in the unusual enol tautomeric form or that the adenine in its N1 protonated form binds the G nucleotide of the ligand in its canonical keto tautomeric state. By using NMR spectroscopy we demonstrate that the c-GAMP riboswitches bind c-di-GMP using a stably protonated adenine in the ligand binding pocket. Thereby, we provide novel insights into the putative biological functions of protonated nucleotides in RNA, which in this case influence the ligand selectivity in a riboswitch.

Protonated nucleobases are not fully ionized in their chloride salt crystals and form metastable base pairs further stabilized by the surrounding anions

This paper presents experimental charge-density studies of cytosinium chloride, adeninium chloride hemihydrate and guaninium dichloride crystals based on ultra-high-resolution X-ray diffraction data and extensive theoretical calculations. The results confirm that the cohesive energies of the studied systems are dominated by contributions from intermolecular electrostatic interactions, as expected for ionic crystals. Electrostatic interaction energies (Ees) usually constitute 95% of the total interaction energy. The Ees energies in this study were several times larger in absolute value when compared, for example, with dimers of neutral nucleobases. However, they were not as large as some theoretical calculations have predicted. This was because the molecules appeared not to be fully ionized in the studied crystals. Apart from charge transfer from chlorine to the protonated nucleobases, small but visible charge redistribution within the nucleobase cations was observed. Some dimers of singly protonated bases in the studied crystals, namely a cytosinium-cytosinium trans sugar/sugar edge pair and an adeninium-adeninium trans Hoogsteen/Hoogsteen edge pair, exhibited attractive interactions (negative values of Ees) or unusually low repulsion despite identical molecular charges. The pairs are metastable as a result of strong hydrogen bonding between bases which overcompensates the overall cation-cation repulsion, the latter being weakened due to charge transfer and molecular charge-density polarization.

Single-Molecule Mechanical Folding and Unfolding of RNA Hairpins: Effects of Single A-U to A·C Pair Substitutions and Single Proton Binding and Implications for mRNA Structure-Induced -1 Ribosomal Frameshifting

A wobble A·C pair can be protonated at near physiological pH to form a more stable wobble A⁺·C pair. Here, we constructed an RNA hairpin (rHP) and three mutants with one A-U base pair substituted with an A·C mismatch on the top (near the loop, U22C), middle (U25C), and bottom (U29C) positions of the stem, respectively. Our results on single-molecule mechanical (un)folding using optical tweezers reveal the destabilization effect of A-U to A·C pair substitution and protonation-dependent enhancement of mechanical stability facilitated through an increased folding rate, or decreased unfolding rate, or both. Our data show that protonation may occur rapidly upon the formation of an apparent mechanical folding transition state. Furthermore, we measured the bulk -1 ribosomal frameshifting efficiencies of the hairpins by a cell-free translation assay. For the mRNA hairpins studied, -1 frameshifting efficiency correlates with mechanical unfolding force at equilibrium and folding rate at around 15 pN. U29C has a frameshifting efficiency similar to that of rHP (∼2%). Accordingly, the bottom 2-4 base pairs of U29C may not form under a stretching force at pH 7.3, which is consistent with the fact that the bottom base pairs of the hairpins may be disrupted by ribosome at the slippery site. U22C and U25C have a similar frameshifting efficiency (∼1%), indicating that both unfolding and folding rates of an mRNA hairpin in a crowded environment may affect frameshifting. Our data indicate that mechanical (un)folding of RNA hairpins may mimic how mRNAs unfold and fold in the presence of translating ribosomes.

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.

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.

The structure of an E. coli tRNAfMet A1-U72 variant shows an unusual conformation of the A1-U72 base pair

Translation initiation in eukaryotes and archaea involves a methionylated initiator tRNA delivered to the ribosome in a ternary complex with e/aIF2 and GTP. Eukaryotic and archaeal initiator tRNAs contain a highly conserved A₁-U₇₂ base pair at the top of the acceptor stem. The importance of this base pair to discriminate initiator tRNAs from elongator tRNAs has been established previously using genetics and biochemistry. However, no structural data illustrating how the A₁-U₇₂ base pair participates in the accurate selection of the initiator tRNAs by the translation initiation systems are available. Here, we describe the crystal structure of a mutant E. coli initiator tRNA_f^MetA₁-U₇₂, aminoacylated with methionine, in which the C₁:A₇₂ mismatch at the end of the tRNA acceptor stem has been changed to an A₁-U₇₂ base pair. Sequence alignments show that the mutant E. coli tRNA is a good mimic of archaeal initiator tRNAs. The crystal structure, determined at 2.8 Å resolution, shows that the A₁-U₇₂ pair adopts an unusual arrangement. A₁ is in a syn conformation and forms a single H-bond interaction with U₇₂ This interaction requires protonation of the N1 atom of A₁ Moreover, the 5' phosphoryl group folds back into the major groove of the acceptor stem and interacts with the N7 atom of G₂ A possible role of this unusual geometry of the A₁-U₇₂ pair in the recognition of the initiator tRNA by its partners during eukaryotic and archaeal translation initiation is discussed.

How to understand atomistic molecular dynamics simulations of RNA and protein-RNA complexes?

We provide a critical assessment of explicit-solvent atomistic molecular dynamics (MD) simulations of RNA and protein/RNA complexes, written primarily for non-specialists with an emphasis to explain the limitations of MD. MD simulations can be likened to hypothetical single-molecule experiments starting from single atomistic conformations and investigating genuine thermal sampling of the biomolecules. The main advantage of MD is the unlimited temporal and spatial resolution of positions of all atoms in the simulated systems. Fundamental limitations are the short physical time-scale of simulations, which can be partially alleviated by enhanced-sampling techniques, and the highly approximate atomistic force fields describing the simulated molecules. The applicability and present limitations of MD are demonstrated on studies of tetranucleotides, tetraloops, ribozymes, riboswitches and protein/RNA complexes. Wisely applied simulations respecting the approximations of the model can successfully complement structural and biochemical experiments. WIREs RNA 2017, 8:e1405. doi: 10.1002/wrna.1405 For further resources related to this article, please visit the WIREs website.

Catalytic DNA: Scope, Applications, and Biochemistry of Deoxyribozymes

The discovery of natural RNA enzymes (ribozymes) prompted the pursuit of artificial DNA enzymes (deoxyribozymes) by in vitro selection methods. A key motivation is the conceptual and practical advantages of DNA relative to proteins and RNA. Early studies focused on RNA-cleaving deoxyribozymes, and more recent experiments have expanded the breadth of catalytic DNA to many other reactions. Including modified nucleotides has the potential to widen the scope of DNA enzymes even further. Practical applications of deoxyribozymes include their use as sensors for metal ions and small molecules. Structural studies of deoxyribozymes are only now beginning; mechanistic experiments will surely follow. Following the first report 21 years ago, the field of deoxyribozymes has promise for both fundamental and applied advances in chemistry, biology, and other disciplines.

Anions in Nucleic Acid Crystallography

Nucleic acid crystallization buffers contain a large variety of chemicals fitting specific needs. Among them, anions are often solely considered for pH-regulating purposes and as cationic co-salts while their ability to directly bind to nucleic acid structures is rarely taken into account. Here we review current knowledge related to the use of anions in crystallization buffers along with data on their biological prevalence. Chloride ions are frequently identified in crystal structures but display low cytosolic concentrations. Hence, they are thought to be distant from nucleic acid structures in the cell. Sulfate ions are also frequently identified in crystal structures but their localization in the cell remains elusive. Nevertheless, the characterization of the binding properties of these ions is essential for better interpreting the solvent structure in crystals and consequently, avoiding mislabeling of electron densities. Furthermore, understanding the binding properties of these anions should help to get clues related to their potential effects in crowded cellular environments.

A DNAzyme requiring two different metal ions at two distinct sites

Most previously reported RNA-cleaving DNAzymes require only a single divalent metal ion for catalysis. We recently reported a general trivalent lanthanide-dependent DNAzyme named Ce13d. This work shows that Ce13d requires both Na⁺ and a trivalent lanthanide (e.g. Ce³⁺), simultaneously. This discovery is facilitated by the sequence similarity between Ce13d and a recently reported Na⁺-specific DNAzyme, NaA43. The Ce13d cleavage rate linearly depends on the concentration of both metal ions. Sensitized Tb³⁺ luminescence and DMS footprinting experiments indicate that the guanines in the enzyme loop are important for Na⁺-binding. The Na⁺ dissociation constants of Ce13d measured from the cleavage activity assay, Tb³⁺ luminescence and DMS footprinting are 24.6, 16.3 and 47 mM, respectively. Mutation studies indicate that the role of Ce³⁺ might be replaced by G₂₃ in NaA43. Ce³⁺ functions by stabilizing the transition state phosphorane, thus promoting cleavage. G₂₃ competes favorably with low concentration Ce³⁺ (below 1 μM). The G₂₃-to-hypoxanthine mutation suggests the N1 position of the guanine as a hydrogen bond donor. Together, Ce13d has two distinct metal binding sites, each fulfilling a different role. DNAzymes can be quite sophisticated in utilizing metal ions for catalysis and molecular recognition, similar to protein metalloenzymes.

How do metal ions direct ribozyme folding?

Ribozymes, which carry out phosphoryl-transfer reactions, often require Mg(2+) ions for catalytic activity. The correct folding of the active site and ribozyme tertiary structure is also regulated by metal ions in a manner that is not fully understood. Here we employ coarse-grained molecular simulations to show that individual structural elements of the group I ribozyme from the bacterium Azoarcus form spontaneously in the unfolded ribozyme even at very low Mg(2+) concentrations, and are transiently stabilized by the coordination of Mg(2+) ions to specific nucleotides. However, competition for scarce Mg(2+) and topological constraints that arise from chain connectivity prevent the complete folding of the ribozyme. A much higher Mg(2+) concentration is required for complete folding of the ribozyme and stabilization of the active site. When Mg(2+) is replaced by Ca(2+) the ribozyme folds, but the active site remains unstable. Our results suggest that group I ribozymes utilize the same interactions with specific metal ligands for both structural stability and chemical activity.

The effects of tautomerization and protonation on the adenine-cytosine mismatches: a density functional theory study

In the present work, we demonstrate the results of a theoretical study concerned with the question how tautomerization and protonation of adenine affect the various properties of adenine-cytosine mismatches. The calculations, in gas phase and in water, are performed at B3LYP/6-311++G(d,p) level. In gas phase, it is observed that any tautomeric form of investigated mismatches is more stabilized when adenine is protonated. As for the neutral mismatches, the mismatches containing amino form of cytosine and imino form of protonated adenine are more stable. The role of aromaticity on the stability of tautomeric forms of mismatches is investigated by NICS(1)ZZ index. The stability of mispairs decreases by going from gas phase to water. It can be explained using dipole moment parameter. The influence of hydrogen bonds on the stability of mismatches is examined by atoms in molecules and natural bond orbital analyses. In addition to geometrical parameters and binding energies, the study of the topological properties of electron charge density aids in better understanding of these mispairs.

Protonation-Dependent Base Flipping at Neutral pH in the Catalytic Triad of a Self-Splicing Bacterial Group II Intron

NMR spectroscopy has revealed pH-dependent structural changes in the highly conserved catalytic domain 5 of a bacterial group II intron. Two adenines with pK(a) values close to neutral pH were identified in the catalytic triad and the bulge. Protonation of the adenine opposite to the catalytic triad is stabilized within a G(syn)-AH(+) (anti) base pair. The pH-dependent anti-to-syn flipping of this G in the catalytic triad modulates the known interaction with the linker region between domains 2 and 3 (J23) and simultaneously the binding of the catalytic Mg(2+) ion to its backbone. Hence, this here identified shifted pK(a) value controls the conformational change between the two steps of splicing.

8-Azapurines as isosteric purine fluorescent probes for nucleic acid and enzymatic research

The 8-azapurines, and their 7-deaza and 9-deaza congeners, represent a unique class of isosteric (isomorphic) analogues of the natural purines, frequently capable of substituting for the latter in many biochemical processes. Particularly interesting is their propensity to exhibit pH-dependent room-temperature fluorescence in aqueous medium, and in non-polar media. We herein review the physico-chemical properties of this class of compounds, with particular emphasis on the fluorescence emission properties of their neutral and/or ionic species, which has led to their widespread use as fluorescent probes in enzymology, including enzymes involved in purine metabolism, agonists/antagonists of adenosine receptors, mechanisms of catalytic RNAs, RNA editing, etc. They are also exceptionally useful fluorescent probes for analytical and clinical applications in crude cell homogenates.

Feasibility of occurrence of different types of protonated base pairs in RNA: a quantum chemical study

Protonated nucleobases have significant roles in facilitating catalytic functions of RNA, and in stabilizing different structural motifs. Reported pKa values of nucleobase protonation suggest that the population of neutral nucleobases is 10(3)-10(4) times higher than that of protonated nucleobases under physiological conditions (pH ∼ 7.4). Therefore, a molecular level understanding of various putative roles of protonated nucleobases cannot be achieved without addressing the question of how their occurrence propensities and stabilities are related to the free energy costs associated with the process of protonation under physiological conditions. With water as the proton donor, we use advanced QM methods to evaluate the site specific protonation propensities of nucleobases in terms of their associated free energy changes (ΔGprot). Quantitative follow up on the energetics of base pair formation and database search for evaluating their occurrence frequencies, reveal a lack of correlation between base pair stability and occurrence propensities on the one hand, and ease of protonation on the other. For example, although N7 protonated adenine (ΔGprot = 40.0 kcal mol(-1)) is found to participate in stable base pairing, base pairs involving N7 protonated guanine (ΔGprot = 36.8 kcal mol(-1)), on geometry optimization, converge to a minima where guanine transfers its extra proton to its partner base. Such observations, along with examples of weak base pairs involving N3 protonation of cytosine (ΔGprot = 37.0 kcal mol(-1)) are rationalized by analysing the protonation induced charge redistributions which are found to significantly influence, both positively and negatively, the hydrogen bonding potentials of different functional sites of individual nucleobases. Protonation induced charge redistribution is also found to strongly influence (i) the aromatic character of the rings of the participating bases and (ii) hydrogen bonding potential of the free edges of the protonated base pair. Comprehensive analysis of a non-redundant RNA crystal structure dataset further reveals that, while availability of stabilization possibilities determine the feasibility of occurrence of protonated bases, their occurrence context and specific functional roles are important factors determining their occurrence propensities.

Determination of hepatitis delta virus ribozyme N(-1) nucleobase and functional group specificity using internal competition kinetics

Biological catalysis involves interactions distant from the site of chemistry that can position the substrate for reaction. Catalysis of RNA 2'-O-transphosphorylation by the hepatitis delta virus (HDV) ribozyme is sensitive to the identity of the N(-1) nucleotide flanking the reactive phosphoryl group. However, the interactions that affect the conformation of this position, and in turn the 2'O nucleophile, are unclear. Here, we describe the application of multiple substrate internal competition kinetic analyses to understand how the N(-1) nucleobase contributes to HDV catalysis and test the utility of this approach for RNA structure-function studies. Internal competition reactions containing all four substrate sequence variants at the N(-1) position in reactions using ribozyme active site mutations at A77 and A78 were used to test a proposed base-pairing interaction. Mutants A78U, A78G, and A79G retain significant catalytic activity but do not alter the specificity for the N(-1) nucleobase. Effects of nucleobase analog substitutions at N(-1) indicate that U is preferred due to the ability to donate an H-bond in the Watson-Crick face and avoid minor groove steric clash. The results provide information essential for evaluating models of the HDV active site and illustrate multiple substrate kinetic analyses as a practical approach for characterizing structure-function relationships in RNA reactions.

Integration of kinetic isotope effect analyses to elucidate ribonuclease mechanism

The well-studied mechanism of ribonuclease A is believed to involve concerted general acid-base catalysis by two histidine residues, His12 and His119. The basic features of this mechanism are often cited to explain rate enhancement by both protein and RNA enzymes that catalyze RNA 2'-O-transphosphorylation. Recent kinetic isotope effect analyses and computational studies are providing a more chemically detailed description of the mechanism of RNase A and the rate limiting transition state. Overall, the results support an asynchronous mechanism for both solution and ribonuclease catalyzed reactions in which breakdown of a transient dianoinic phosphorane intermediate by 5'OP bond cleavage is rate limiting. Relative to non-enzymatic reactions catalyzed by specific base, a smaller KIE on the 5'O leaving group and a less negative βLG are observed for RNase A catalysis. Quantum mechanical calculations consistent with these data support a model in which electrostatic and H-bonding interactions with the non-bridging oxygens and proton transfer from His119 render departure of the 5'O less advanced and stabilize charge buildup in the transition state. Both experiment and computation indicate advanced 2'OP bond formation in the rate limiting transition state. However, this feature makes it difficult to resolve the chemical steps involved in 2'O activation. Thus, modeling the transition state for RNase A catalysis underscores those elements of its chemical mechanism that are well resolved, as well as highlighting those where ambiguity remains. This article is part of a Special Issue entitled: Enzyme Transition States from Theory and Experiment.

Multiscale methods for computational RNA enzymology

RNA catalysis is of fundamental importance to biology and yet remains ill-understood due to its complex nature. The multidimensional "problem space" of RNA catalysis includes both local and global conformational rearrangements, changes in the ion atmosphere around nucleic acids and metal ion binding, dependence on potentially correlated protonation states of key residues, and bond breaking/forming in the chemical steps of the reaction. The goal of this chapter is to summarize and apply multiscale modeling methods in an effort to target the different parts of the RNA catalysis problem space while also addressing the limitations and pitfalls of these methods. Classical molecular dynamics simulations, reference interaction site model calculations, constant pH molecular dynamics (CpHMD) simulations, Hamiltonian replica exchange molecular dynamics, and quantum mechanical/molecular mechanical simulations will be discussed in the context of the study of RNA backbone cleavage transesterification. This reaction is catalyzed by both RNA and protein enzymes, and here we examine the different mechanistic strategies taken by the hepatitis delta virus ribozyme and RNase A.

The role of N7 protonation of guanine in determining the structure, stability and function of RNA base pairs

Ab initio computations and bioinformatics studies reveal that stabilization of some important RNA structural motifs might involve N7 protonation of guanine.

Role of tautomerism in RNA biochemistry

Heterocyclic nucleic acid bases and their analogs can adopt multiple tautomeric forms due to the presence of multiple solvent-exchangeable protons. In DNA, spontaneous formation of minor tautomers has been speculated to contribute to mutagenic mispairings during DNA replication, whereas in RNA, minor tautomeric forms have been proposed to enhance the structural and functional diversity of RNA enzymes and aptamers. This review summarizes the role of tautomerism in RNA biochemistry, specifically focusing on the role of tautomerism in catalysis of small self-cleaving ribozymes and recognition of ligand analogs by riboswitches. Considering that the presence of multiple tautomers of nucleic acid bases is a rare occurrence, and that tautomers typically interconvert on a fast time scale, methods for studying rapid tautomerism in the context of nucleic acids under biologically relevant aqueous conditions are also discussed.

Molecular dynamic simulations of protein/RNA complexes: CRISPR/Csy4 endoribonuclease

BACKGROUND

Many prokaryotic genomes comprise Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) offering defense against foreign nucleic acids. These immune systems are conditioned by the production of small CRISPR-derived RNAs matured from long RNA precursors. This often requires a Csy4 endoribonuclease cleaving the RNA 3'-end.

METHODS

We report extended explicit solvent molecular dynamic (MD) simulations of Csy4/RNA complex in precursor and product states, based on X-ray structures of product and inactivated precursor (55 simulations; ~3.7μs in total).

RESULTS

The simulations identify double-protonated His29 and deprotonated terminal phosphate as the likely dominant protonation states consistent with the product structure. We revealed potential substates consistent with Ser148 and His29 acting as the general base and acid, respectively. The Ser148 could be straightforwardly deprotonated through solvent and could without further structural rearrangements deprotonate the nucleophile, contrasting similar studies investigating the general base role of nucleobases in ribozymes. We could not locate geometries consistent with His29 acting as general base. However, we caution that the X-ray structures do not always capture the catalytically active geometries and then the reactive structures may be unreachable by the simulation technique.

CONCLUSIONS

We identified potential catalytic arrangement of the Csy4/RNA complex but we also report limitations of the simulation technique. Even for the dominant protonation state we could not achieve full agreement between the simulations and the structural data.

GENERAL SIGNIFICANCE

Potential catalytic arrangement of the Csy4/RNA complex is found. Further, we provide unique insights into limitations of simulations of protein/RNA complexes, namely, the influence of the starting experimental structures and force field limitations. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Building a stable RNA U-turn with a protonated cytidine

The U-turn is a classical three-dimensional RNA folding motif first identified in the anticodon and T-loops of tRNAs. It also occurs frequently as a building block in other functional RNA structures in many different sequence and structural contexts. U-turns induce sharp changes in the direction of the RNA backbone and often conform to the 3-nt consensus sequence 5'-UNR-3' (N = any nucleotide, R = purine). The canonical U-turn motif is stabilized by a hydrogen bond between the N3 imino group of the U residue and the 3' phosphate group of the R residue as well as a hydrogen bond between the 2'-hydroxyl group of the uridine and the N7 nitrogen of the R residue. Here, we demonstrate that a protonated cytidine can functionally and structurally replace the uridine at the first position of the canonical U-turn motif in the apical loop of the neomycin riboswitch. Using NMR spectroscopy, we directly show that the N3 imino group of the protonated cytidine forms a hydrogen bond with the backbone phosphate 3' from the third nucleotide of the U-turn analogously to the imino group of the uridine in the canonical motif. In addition, we compare the stability of the hydrogen bonds in the mutant U-turn motif to the wild type and describe the NMR signature of the C+-phosphate interaction. Our results have implications for the prediction of RNA structural motifs and suggest simple approaches for the experimental identification of hydrogen bonds between protonated C-imino groups and the phosphate backbone.

Altered (transition) states: mechanisms of solution and enzyme catalyzed RNA 2'-O-transphosphorylation

Although there have been great strides in defining the mechanisms of RNA strand cleavage by 2'-O-transphosphorylation, long-standing questions remain. How do different catalytic modes such as acid/base and metal ion catalysis influence transition state charge distribution? Does the large rate enhancement characteristic of biological catalysis result in different transition states relative to solution reactions? Answering these questions is important for understanding biological catalysis in general, and revealing principles for designing small molecule inhibitors. Recent application of linear free energy relationships and kinetic isotope effects together with multi-scale computational simulations are providing tentative answers to these questions for this fundamentally important class of phosphoryl transfer reactions.