Molecular dynamics suggest multifunctionality of an adenine imino group in acid-base catalysis of the hairpin ribozyme

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.

Thermodynamics and kinetics of RNA tertiary structure formation in the junctionless hairpin ribozyme

The hairpin ribozyme consists of two RNA internal loops that interact to form the catalytically active structure. This docking transition is a rare example of intermolecular formation of RNA tertiary structure without coupling to helix annealing. We have used temperature-dependent surface plasmon resonance (SPR) to characterize the thermodynamics and kinetics of RNA tertiary structure formation for the junctionless form of the ribozyme, in which loops A and B reside on separate molecules. We find docking to be strongly enthalpy-driven and to be accompanied by substantial activation barriers for association and dissociation, consistent with the structural reorganization of both internal loops upon complex formation. Comparisons with the parallel analysis of a ribozyme variant carrying a 2'-O-methyl modification at the self-cleavage site and with published data in other systems reveal a surprising diversity of thermodynamic signatures, emphasizing the delicate balance of contributions to the free energy of formation of RNA tertiary structure.

Intrinsic Base-Pair Rearrangement in the Hairpin Ribozyme Directs RNA Conformational Sampling and Tertiary Interface Formation

Dynamic fluctuations in RNA structure enable conformational changes that are required for catalysis and recognition. In the hairpin ribozyme, the catalytically active structure is formed as an intricate tertiary interface between two RNA internal loops. Substantial alterations in the structure of each loop are observed upon interface formation, or docking. The very slow on-rate for this relatively tight interaction has led us to hypothesize a double conformational capture mechanism for RNA-RNA recognition. We used extensive molecular dynamics simulations to assess conformational sampling in the undocked form of the loop domain containing the scissile phosphate (loop A). We observed several major accessible conformations with distinctive patterns of hydrogen bonding and base stacking interactions in the active-site internal loop. Several important conformational features characteristic of the docked state were observed in well-populated substates, consistent with the kinetic sampling of docking-competent states by isolated loop A. Our observations suggest a hybrid or multistage binding mechanism, in which initial conformational selection of a docking-competent state is followed by induced-fit adjustment to an in-line, chemically reactive state only after formation of the initial complex with loop B.

Insight into G-quadruplex-hemin DNAzyme/RNAzyme: adjacent adenine as the intramolecular species for remarkable enhancement of enzymatic activity

G-quadruplex (G4) with stacked G-tetrads structure is able to bind hemin (iron (III)-protoporphyrin IX) to form a unique type of DNAzyme/RNAzyme with peroxidase-mimicking activity, which has been widely employed in multidisciplinary fields. However, its further applications are hampered by its relatively weak activity compared with protein enzymes. Herein, we report a unique intramolecular enhancement effect of the adjacent adenine (EnEAA) at 3' end of G4 core sequences that significantly improves the activity of G4 DNAzymes. Through detailed investigations of the EnEAA, the added 3' adenine was proved to accelerate the compound I formation in catalytic cycle and thus improve the G4 DNAzyme activity. EnEAA was found to be highly dependent on the unprotonated state of the N1 of adenine, substantiating that adenine might function as a general acid-base catalyst. Further adenine analogs analysis supported that both N1 and exocyclic 6-amino groups in adenine played key role in the catalysis. Moreover, we proved that EnEAA was generally applicable for various parallel G-quadruplex structures and even G4 RNAzyme. Our studies implied that adenine might act analogously as the distal histidine in protein peroxidases, which shed light on the fundamental understanding and rational design of G4 DNAzyme/RNAzyme catalysts with enhanced functions.

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.

Wobble pairs of the HDV ribozyme play specific roles in stabilization of active site dynamics

The hepatitis delta virus (HDV) is the only known human pathogen whose genome contains a catalytic RNA motif (ribozyme). The overall architecture of the HDV ribozyme is that of a double-nested pseudoknot, with two GU pairs flanking the active site. Although extensive studies have shown that mutation of either wobble results in decreased catalytic activity, little work has focused on linking these mutations to specific structural effects on catalytic fitness. Here we use molecular dynamics simulations based on an activated structure to probe the active site dynamics as a result of wobble pair mutations. In both wild-type and mutant ribozymes, the in-line fitness of the active site (as a measure of catalytic proficiency) strongly depends on the presence of a C75(N3H3+)N1(O5') hydrogen bond, which positions C75 as the general acid for the reaction. Our mutational analyses show that each GU wobble supports catalytically fit conformations in distinct ways; the reverse G25U20 wobble promotes high in-line fitness, high occupancy of the C75(N3H3+)G1(O5') general-acid hydrogen bond and stabilization of the G1U37 wobble, while the G1U37 wobble acts more locally by stabilizing high in-line fitness and the C75(N3H3+)G1(O5') hydrogen bond. We also find that stable type I A-minor and P1.1 hydrogen bonding above and below the active site, respectively, prevent local structural disorder from spreading and disrupting global conformation. Taken together, our results define specific, often redundant architectural roles for several structural motifs of the HDV ribozyme active site, expanding the known roles of these motifs within all HDV-like ribozymes and other structured RNAs.

Dynamics of riboswitches: Molecular simulations

Riboswitch RNAs play key roles in bacterial metabolism and represent a promising new class of antibiotic targets for treatment of infectious disease. While many studies of riboswitches have been performed, the exact mechanism of riboswitch operation is still not fully understood at the atomistic level of detail. Molecular dynamics simulations are useful for interpreting existing experimental data and producing predictions for new experiments. Here, a wide range of computational studies on riboswitches is reviewed. By elucidating the key principles of riboswitch operation, computation may aid in the effort to design more specific antibiotics with affinities greater than those of the native ligand. Such a detailed understanding may be required to improve efficacy and reduce side effects. These studies are laying the groundwork for understanding the action mechanism of new compounds that inhibit riboswitch activity. Future directions such as magnesium effects, large-scale conformational changes, expression platforms and co-transcriptional folding are also discussed. This article is part of a Special Issue entitled: Riboswitches.

Comparison of ab Initio, DFT, and Semiempirical QM/MM Approaches for Description of Catalytic Mechanism of Hairpin Ribozyme

We have analyzed the capability of state-of-the-art multiscale computational approaches to provide atomic-resolution electronic structure insights into possible catalytic scenarios of the hairpin ribozyme by evaluating potential and free energy surfaces of the reactions by various hybrid QM/MM methods. The hairpin ribozyme is a unique catalytic RNA that achieves rate acceleration similar to other small self-cleaving ribozymes but without direct metal ion participation. Guanine 8 (G8) and adenine 38 (A38) have been identified as the catalytically essential nucleobases. However, their exact catalytic roles are still being investigated. In line with the available experimental data, we considered two reaction scenarios involving protonated A38H(+) as a general acid which is further assisted by either canonical G8 or deprotonated G8(-) forms. We used the spin-component scaled Møller-Plesset (SCS-MP2) method at the complete basis set limit as the reference method. The semiempirical AM1/d-PhoT and SCC-DFTBPR methods provided acceptable activation barriers with respect to the SCS-MP2 data but predicted significantly different reaction pathways. DFT functionals (BLYP and MPW1K) yielded the same reaction pathway as the SCS-MP2 method. The activation barriers were slightly underestimated by the GGA BLYP functional, although with accuracy comparable to the semiempirical methods. The SCS-MP2 method and hybrid MPW1K functional gave activation barriers that were closest to those derived from experimentally measured rate constants.

Principles of ion recognition in RNA: insights from the group II intron structures

Metal ions promote both RNA folding and catalysis, thus being essential in stabilizing the structure and determining the function of large RNA molecules, including group II introns. The latter are self-splicing metalloribozymes, containing a heteronuclear four-metal-ion center within the active site. In addition to these catalytic ions, group II introns bind many other structural ions, including delocalized ions that bind the RNA diffusively and well-ordered ions that bind the RNA tightly with high occupancy. The latter ions, which can be studied by biophysical methods, have not yet been analyzed systematically. Here, we compare crystal structures of the group IIC intron from Oceanobacillus iheyensis and classify numerous site-bound ions, which are primarily localized in the intron core and near long-range tertiary contacts. Certain ion-binding sites resemble motifs observed in known RNA structures, while others are idiosyncratic to the group II intron. Particularly interesting are (1) ions proximal to the active site, which may participate in splicing together with the catalytic four-metal-ion center, (2) organic ions that bind regions predicted to interact with intron-encoded proteins, and (3) unusual monovalent ions bound to GU wobble pairs, GA mismatches, the S-turn, the tetraloop-receptor, and the T-loop. Our analysis extends the general principles by which ions participate in RNA structural organization and it will aid in the determination and interpretation of future RNA structures.

Automated modeling of RNA 3D structure

This chapter gives an overview over the current methods for automated modeling of RNA structures, with emphasis on template-based methods. The currently used approaches to RNA modeling are presented with a side view on the protein world, where many similar ideas have been used. Two main programs for automated template-based modeling are presented: ModeRNA assembling structures from fragments and MacroMoleculeBuilder performing a simulation to satisfy spatial restraints. Both approaches have in common that they require an alignment of the target sequence to a known RNA structure that is used as a modeling template. As a way to find promising template structures and to align the target and template sequences, we propose a pipeline combining the ParAlign and Infernal programs on RNA family data from Rfam. We also briefly summarize template-free methods for RNA 3D structure prediction. Typically, RNA structures generated by automated modeling methods require local or global optimization. Thus, we also discuss methods that can be used for local or global refinement of RNA structures.

RNA structure and dynamics: a base pairing perspective

RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson-Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson-Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.

Towards Accurate Prediction of Protonation Equilibrium of Nucleic Acids

The role of protonated nucleotides in modulating the pH-dependent properties of nucleic acids is one of the emerging frontiers in the field of nucleic acid biology. The recent development of a constant pH molecular dynamics simulation (CPHMD^MSλD) framework for simulating nucleic acids has provided a tool for realistic simulations of pH-dependent dynamics. We enhanced the CPHMD^MSλD framework with pH-based replica exchange (pH-REX), which significantly improves the sampling of both titration and spatial coordinates. The results from our pK_a calculations for the GAAA tetraloop, which was predicted with lower accuracy previously due to sampling challenges, demonstrates that pH-REX reduces the average unsigned error (AUE) to 0.7 pK_a units, and the error of the most poorly predicted residue A17 was drastically reduced from 2.9 to 1.2 pK_a unit. Lastly, we show that pH-REX CPHMD^MSλD simulations can be used to identify the dominant conformation of nucleic acid structures in alternate pH environments. This work suggests that pH-REX CPHMD^MSλD simulations provide a practical tool for predicting nucleic acid protonation equilibrium from first-principles, and offering structural and mechanistic insight into the study of pH-dependent properties of nucleic acids.

pH-dependent dynamics of complex RNA macromolecules

The role of pH-dependent protonation equilibrium in modulating RNA dynamics and function is one of the key unanswered questions in RNA biology. Molecular dynamics (MD) simulations can provide insight into the mechanistic roles of protonated nucleotides, but it is only capable of modeling fixed protonation states and requires prior knowledge of the key residue's protonation state. Recently, we developed a framework for constant pH molecular dynamics simulations (CPHMDMSλD) of nucleic acids, where the nucleotides' protonation states are modeled as dynamic variables that are coupled to the structural dynamics of the RNA. In the present study, we demonstrate the application of CPHMDMSλD to the lead-dependent ribozyme; establishing the validity of this approach for modeling complex RNA structures. We show that CPHMDMSλD accurately predicts the direction of the pKa shifts and reproduces experimentally-measured microscopic pKa values with an average unsigned error of 1.3 pKa units. The effects of coupled titration states in RNA structures are modeled, and the importance of conformation sampling is highlighted. The general accuracy of CPHMDMSλD simulations in reproducing pH-dependent observables reported in this work demonstrates that constant pH simulations provides a powerful tool to investigate pH-dependent processes in nucleic acids.

Lewis acid catalysis of phosphoryl transfer from a copper(II)-NTP complex in a kinase ribozyme

The chemical strategies used by ribozymes to enhance reaction rates are revealed in part from their metal ion and pH requirements. We find that kinase ribozyme K28(1-77)C, in contrast with previously characterized kinase ribozymes, requires Cu²⁺ for optimal catalysis of thiophosphoryl transfer from GTPγS. Phosphoryl transfer from GTP is greatly reduced in the absence of Cu²⁺, indicating a specific catalytic role independent of any potential interactions with the GTPγS thiophosphoryl group. In-line probing and ATPγS competition both argue against direct Cu²⁺ binding by RNA; rather, these data establish that Cu²⁺ enters the active site within a Cu²⁺•GTPγS or Cu²⁺•GTP chelation complex, and that Cu²⁺•nucleobase interactions further enforce Cu²⁺ selectivity and position the metal ion for Lewis acid catalysis. Replacing Mg²⁺ with [Co(NH₃)₆]³⁺ significantly reduced product yield, but not k_obs, indicating that the role of inner-sphere Mg²⁺ coordination is structural rather than catalytic. Replacing Mg²⁺ with alkaline earths of increasing ionic radii (Ca²⁺, Sr²⁺ and Ba²⁺) gave lower yields and approximately linear rates of product accumulation. Finally, we observe that reaction rates increased with pH in log-linear fashion with an apparent pKa = 8.0 ± 0.1, indicating deprotonation in the rate-limiting step.

The DNA and RNA sugar-phosphate backbone emerges as the key player. An overview of quantum-chemical, structural biology and simulation studies

Knowledge of geometrical and physico-chemical properties of the sugar-phosphate backbone substantially contributes to the comprehension of the structural dynamics, function and evolution of nucleic acids. We provide a side by side overview of structural biology/bioinformatics, quantum chemical and molecular mechanical/simulation studies of the nucleic acids backbone. We highlight main features, advantages and limitations of these techniques, with a special emphasis given to their synergy. The present status of the research is then illustrated by selected examples which include classification of DNA and RNA backbone families, benchmark structure-energy quantum chemical calculations, parameterization of the dihedral space of simulation force fields, incorporation of arsenate into DNA, sugar-phosphate backbone self-cleavage in small RNA enzymes, and intricate geometries of the backbone in recurrent RNA building blocks. Although not apparent from the current literature showing limited overlaps between the QM, simulation and bioinformatics studies of the nucleic acids backbone, there in fact should be a major cooperative interaction between these three approaches in studies of the sugar-phosphate backbone.

General acid-base catalysis mediated by nucleobases in the hairpin ribozyme

The catalytic mechanism by which the hairpin ribozyme accelerates cleavage or ligation of the phosphodiester backbone of RNA has been incompletely understood. There is experimental evidence for an important role for an adenine (A38) and a guanine (G8), and it has been proposed that these act in general acid-base catalysis. In this work we show that a large reduction in cleavage rate on substitution of A38 by purine (A38P) can be reversed by replacement of the 5'-oxygen atom at the scissile phosphate by sulfur (5'-PS), which is a much better leaving group. This is consistent with A38 acting as the general acid in the unmodified ribozyme. The rate of cleavage of the 5'-PS substrate by the A38P ribozyme increases with pH log-linearly, indicative of a requirement for a deprotonated base with a relatively high pK(a). On substitution of G8 by diaminopurine, the 5'-PS substrate cleavage rate at first increases with pH and then remains at a plateau, exhibiting an apparent pK(a) consistent with this nucleotide acting in general base catalysis. Alternative explanations for the pH dependence of hairpin ribozyme reactivity are discussed, from which we conclude that general acid-base catalysis by A38 and G8 is the simplest and most probable explanation consistent with all the experimental data.

Constant pH Molecular Dynamics Simulations of Nucleic Acids in Explicit Solvent

The nucleosides of adenine and cytosine have pKa values of 3.50 and 4.08, respectively, and are assumed to be unprotonated under physiological conditions. However, evidence from recent NMR and X-Ray crystallography studies has revealed the prevalence of protonated adenine and cytosine in RNA macromolecules. Such nucleotides with elevated pKa values may play a role in stabilizing RNA structure and participate in the mechanism of ribozyme catalysis. With the work presented here, we establish the framework and demonstrate the first constant pH MD simulations (CPHMD) for nucleic acids in explicit solvent in which the protonation state is coupled to the dynamical evolution of the RNA system via λ-dynamics. We adopt the new functional form λ(Nexp) for λ that was recently developed for Multi-Site λ-Dynamics (MSλD) and demonstrate good sampling characteristics in which rapid and frequent transitions between the protonated and unprotonated states at pH = pKa are achieved. Our calculated pKa values of simple nucleotides are in a good agreement with experimentally measured values, with a mean absolute error of 0.24 pKa units. This work demonstrates that CPHMD can be used as a powerful tool to investigate pH-dependent biological properties of RNA macromolecules.

QM/MM studies of hairpin ribozyme self-cleavage suggest the feasibility of multiple competing reaction mechanisms

The hairpin ribozyme is a prominent member of small ribozymes since it does not require metal ions to achieve catalysis. Guanine 8 (G8) and adenine 38 (A38) have been identified as key participants in self-cleavage and -ligation. We have carried out hybrid quantum-mechanical/molecular mechanical (QM/MM) calculations to evaluate the energy along several putative reaction pathways. The error of our DFT description of the QM region was tested and shown to be ~1 kcal/mol. We find that self-cleavage of the hairpin ribozyme may follow several competing microscopic reaction mechanisms, all with calculated activation barriers in good agreement with those from experiment (20-21 kcal/mol). The initial nucleophilic attack of the A-1(2'-OH) group on the scissile phosphate is predicted to be rate-limiting in all these mechanisms. An unprotonated G8(-) (together with A38H(+)) yields a feasible activation barrier (20.4 kcal/mol). Proton transfer to a nonbridging phosphate oxygen also leads to feasible reaction pathways. Finally, our calculations consider thio-substitutions of one or both nonbridging oxygens of the scissile phosphate and predict that they have only a negligible effect on the reaction barrier, as observed experimentally.

Noncanonical hydrogen bonding in nucleic acids. Benchmark evaluation of key base-phosphate interactions in folded RNA molecules using quantum-chemical calculations and molecular dynamics simulations

RNA molecules are stabilized by a wide range of noncanonical interactions that are not present in DNA. Among them, the recently classified base-phosphate (BPh) interactions belong to the most important ones. Twelve percent of nucleotides in the ribosomal crystal structures are involved in BPh interactions. BPh interactions are highly conserved and provide major constraints on RNA sequence evolution. Here we provide assessment of the energetics of BPh interactions using MP2 computations extrapolated to the complete basis set of atomic orbitals and corrected for higher-order electron correlation effects. The reference computations are compared with DFT-D and DFT-D3 approaches, the SAPT method, and the molecular mechanics force field. The computations, besides providing the basic benchmark for the BPh interactions, allow some refinements of the original classification, including identification of some potential doubly bonded BPh patterns. The reference computations are followed by analysis of some larger RNA fragments that consider the context of the BPh interactions. The computations demonstrate the complexity of interaction patterns utilizing the BPh interactions in real RNA structures. The BPh interactions are often involved in intricate interaction networks. We studied BPh interactions of protonated adenine that can contribute to catalysis of hairpin ribozyme, the key BPh interaction in the S-turn motif of the sarcin-ricin loop, which may predetermine the S-turn topology and complex BPh patterns from the glmS riboswitch. Finally, the structural stability of BPh interactions in explicit solvent molecular dynamics simulations is assessed. The simulations well preserve key BPh interactions and allow dissection of structurally/functionally important water-meditated BPh bridges, which could not be considered in earlier bioinformatics classification of BPh interactions.

The pH dependence of hairpin ribozyme catalysis reflects ionization of an active site adenine

Understanding how self-cleaving ribozymes mediate catalysis is crucial in light of compelling evidence that human and bacterial gene expression can be regulated through RNA self-cleavage. The hairpin ribozyme catalyzes reversible phosphodiester bond cleavage through a mechanism that does not require divalent metal cations. Previous structural and biochemical evidence implicated the amidine group of an active site adenosine, A38, in a pH-dependent step in catalysis. We developed a way to determine microscopic pK(a) values in active ribozymes based on the pH-dependent fluorescence of 8-azaadenosine (8azaA). We compared the microscopic pK(a) for ionization of 8azaA at position 38 with the apparent pK(a) for the self-cleavage reaction in a fully functional hairpin ribozyme with a unique 8azaA at position 38. Microscopic and apparent pK(a) values were virtually the same, evidence that A38 protonation accounts for the decrease in catalytic activity with decreasing pH. These results implicate the neutral unprotonated form of A38 in a transition state that involves formation of the 5'-oxygen-phosphorus bond.

Metal ions: supporting actors in the playbook of small ribozymes

Since the 1980s, several small RNA motifs capable of chemical catalysis have been discovered. These small ribozymes, composed of between approximately 40 and 200 nucleotides, have been found to play vital roles in the replication of subviral and viral pathogens, as well as in gene regulation in prokaryotes, and have recently been discovered in noncoding eukaryotic RNAs. All of the known natural small ribozymes - the hairpin, hammerhead, hepatitis delta virus, Varkud satellite, and glmS ribozymes--catalyze the same self-cleavage reaction as RNase A, resulting in two products, one bearing a 2'-3' cyclic phosphate and the other a 5'-hydroxyl group. Although originally thought to be obligate metalloenzymes like the group I and II self-splicing introns, the small ribozymes are now known to support catalysis in a wide variety of cations that appear to be only indirectly involved in catalysis. Nevertheless, under physiologic conditions, metal ions are essential for the proper folding and function of the small ribozymes, the most effective of these being magnesium. Metal ions contribute to catalysis in the small ribozymes primarily by stabilizing the catalytically active conformation, but in some cases also by activating RNA functional groups for catalysis, directly participating in catalytic acid-base chemistry, and perhaps by neutralizing the developing negative charge of the transition state. Although interactions between the small ribozymes and cations are relatively nonspecific, ribozyme activity is quite sensitive to the types and concentrations of metal ions present in solution, suggesting a close evolutionary relationship between cellular metal ion homeostasis and cation requirements of catalytic RNAs, and perhaps RNA in general.

Deciphering the role of glucosamine-6-phosphate in the riboswitch action of glmS ribozyme

The GlmS ribozyme is believed to exploit a general acid-base catalytic mechanism in the presence of glucosamine-6-phosphate (GlcN6P) to accelerate self-cleavage by approximately six orders of magnitude. The general acid and general base are not known, and the role of the GlcN6P cofactor is even less well understood. The amine group of GlcN6P has the ability to either accept or donate a proton and could therefore potentially act as an acid or a base. In order to decipher the role of GlcN6P in the self-cleavage of glmS, we have determined the preferred protonation state of the amine group in the wild-type and an inactive G40A mutant using molecular dynamics simulations and free energy calculations. Here we show that, upon binding of GlcN6P to wild-type glmS, the pK(a) of the amine moiety is altered by the active site environment, decreasing by about 2.2 from a solution pK(a) of about 8.2. On the other hand, we show that the pK(a) of the amine group slightly increases to about 8.4 upon binding to the G40A inactive mutant of glmS. These results suggest that GlcN6P acts as a general acid in the self-cleavage of glmS. Upon binding to glmS, GlcN6P can easily release a proton to the 5'-oxygen of G1 during self-cleavage of the backbone phosphodiester bond. However, in the G40A inactive mutant of glmS, the results suggest that the ability of GlcN6P to easily release its proton is diminished, in addition to the possible lack of G40 as an effective base.

Structural dynamics of the box C/D RNA kink-turn and its complex with proteins: the role of the A-minor 0 interaction, long-residency water bridges, and structural ion-binding sites revealed by molecular simulations

Kink-turns (K-turns) are recurrent elbow-like RNA motifs that participate in protein-assisted RNA folding and contribute to RNA dynamics. We carried out a set of molecular dynamics (MD) simulations using parm99 and parmbsc0 force fields to investigate structural dynamics of the box C/D RNA and its complexes with two proteins: native archaeal L7ae protein and human 15.5 kDa protein, originally bound to very similar structure of U4 snRNA. The box C/D RNA forms K-turn with A-minor 0 tertiary interaction between its canonical (C) and noncanonical (NC) stems. The local K-turn architecture is thus different from the previously studied ribosomal K-turns 38 and 42 having A-minor I interaction. The simulations reveal visible structural dynamics of this tertiary interaction involving altogether six substates which substantially contribute to the elbow-like flexibility of the K-turn. The interaction can even temporarily shift to the A-minor I type pattern; however, this is associated with distortion of the G/A base pair in the NC-stem of the K-turn. The simulations show reduction of the K-turn flexibility upon protein binding. The protein interacts with the apex of the K-turn and with the NC-stem. The protein-RNA interface includes long-residency hydration sites. We have also found long-residency hydration sites and major ion-binding sites associated with the K-turn itself. The overall topology of the K-turn remains stable in all simulations. However, in simulations of free K-turn, we observed instability of the key C16(O2')-A7(N1) H-bond, which is a signature interaction of K-turns and which was visibly more stable in simulations of K-turns possessing A-minor I interaction. It may reflect either some imbalance of the force field or it may be a correct indication of early stages of unfolding since this K-turn requires protein binding for its stabilization. Interestingly, the 16(O2')-7(N1) H- bond is usually not fully lost since it is replaced by a water bridge with a tightly bound water, which is adenine-specific similarly as the original interaction. The 16(O2')-7(N1) H-bond is stabilized by protein binding. The stabilizing effect is more visible with the human 15.5 kDa protein, which is attributed to valine to arginine substitution in the binding site. The behavior of the A-minor interaction is force-field-dependent because the parmbsc0 force field attenuates the A-minor fluctuations compared to parm99 simulations. Behavior of other regions of the box C/D RNA is not sensitive to the force field choice. Simulation with net-neutralizing Na(+) and 0.2 M excess salt conditions appear in all aspects equivalent. The simulations show loss of a hairpin tetraloop, which is not part of the K-turn. This was attributed to force field limitations.

Protonation states of the key active site residues and structural dynamics of the glmS riboswitch as revealed by molecular dynamics

The glmS catalytic riboswitch is part of the 5'-untranslated region of mRNAs encoding glucosamine-6-phosphate (GlcN6P) synthetase (glmS) in numerous gram-positive bacteria. Binding of the cofactor GlcN6P induces site-specific self-cleavage of the RNA. However, the detailed reaction mechanism as well as the protonation state of the glmS reactive form still remains elusive. To probe the dominant protonation states of key active site residues, we carried out explicit solvent molecular dynamic simulations involving various protonation states of three crucial active site moieties observed in the available crystal structures: (i) guanine G40 (following the Thermoanaerobacter tengcongensis numbering), (ii) the GlcN6P amino/ammonium group, and (iii) the GlcN6P phosphate moiety. We found that a deprotonated G40(-) seems incompatible with the observed glmS active site architecture. Our data suggest that the canonical form of G40 plays a structural role by stabilizing an in-line attack conformation of the cleavage site A-1(2'-OH) nucleophile, rather than a more direct chemical role. In addition, we observe weakened cofactor binding upon protonation of the GlcN6P phosphate moiety, which explains the experimentally observed increase in K(m) with decreasing pH. Finally, we discuss a possible role of cofactor binding and its interaction with the G65 and G1 purines in structural stabilization of the A-1(2'-OH) in-line attack conformation. On the basis of the identified dominant protonation state of the reaction precursor, we propose a hypothesis of the self-cleavage mechanism in which A-1(2'-OH) is activated as a nucleophile by the G1(pro-R(p)) nonbridging oxygen of the scissile phosphate, whereas the ammonium group of GlcN6P acts as the general acid protonating the G1(O5') leaving group.

Extensive molecular dynamics simulations showing that canonical G8 and protonated A38H+ forms are most consistent with crystal structures of hairpin ribozyme

The hairpin ribozyme is a prominent member of the group of small catalytic RNAs (RNA enzymes or ribozymes) because it does not require metal ions to achieve catalysis. Biochemical and structural data have implicated guanine 8 (G8) and adenine 38 (A38) as catalytic participants in cleavage and ligation catalyzed by the hairpin ribozyme, yet their exact role in catalysis remains disputed. To gain insight into dynamics in the active site of a minimal self-cleaving hairpin ribozyme, we have performed extensive classical, explicit-solvent molecular dynamics (MD) simulations on time scales of 50-150 ns. Starting from the available X-ray crystal structures, we investigated the structural impact of the protonation states of G8 and A38, and the inactivating A-1(2'-methoxy) substitution employed in crystallography. Our simulations reveal that a canonical G8 agrees well with the crystal structures while a deprotonated G8 profoundly distorts the active site. Thus MD simulations do not support a straightforward participation of the deprotonated G8 in catalysis. By comparison, the G8 enol tautomer is structurally well tolerated, causing only local rearrangements in the active site. Furthermore, a protonated A38H(+) is more consistent with the crystallography data than a canonical A38. The simulations thus support the notion that A38H(+) is the dominant form in the crystals, grown at pH 6. In most simulations, the canonical A38 departs from the scissile phosphate and substantially perturbs the structures of the active site and S-turn. Yet, we occasionally also observe formation of a stable A-1(2'-OH)...A38(N1) hydrogen bond, which documents the ability of the ribozyme to form this hydrogen bond, consistent with a potential role of A38 as general base catalyst. The presence of this hydrogen bond is, however, incompatible with the expected in-line attack angle necessary for self-cleavage, requiring a rapid transition of the deprotonated 2'-oxyanion to a position more favorable for in-line attack after proton transfer from A-1(2'-OH) to A38(N1). The simulations revealed a potential force field artifact, occasional but irreversible formation of "ladder-like", underwound A-RNA structure in one of the external helices. Although it does not affect the catalytic center of the hairpin ribozyme, further studies are under way to better assess possible influence of such force field behavior on long RNA simulations.

Molecular dynamics and quantum mechanics of RNA: conformational and chemical change we can believe in

Structure and dynamics are both critical to RNA’s vital functions in biology. Numerous techniques can elucidate the structural dynamics of RNA, but computational approaches based on experimental data arguably hold the promise of providing the most detail. In this Account, we highlight areas wherein molecular dynamics (MD) and quantum mechanical (QM) techniques are applied to RNA, particularly in relation to complementary experimental studies.

We have expanded on atomic-resolution crystal structures of RNAs in functionally relevant states by applying explicit solvent MD simulations to explore their dynamics and conformational changes on the submicrosecond time scale. MD relies on simplified atomistic, pairwise additive interaction potentials (force fields). Because of limited sampling, due to the finite accessible simulation time scale and the approximated force field, high-quality starting structures are required.

Despite their imperfection, we find that currently available force fields empower MD to provide meaningful and predictive information on RNA dynamics around a crystallographically defined energy minimum. The performance of force fields can be estimated by precise QM calculations on small model systems. Such calculations agree reasonably well with the Cornell et al. AMBER force field, particularly for stacking and hydrogen-bonding interactions. A final verification of any force field is accomplished by simulations of complex nucleic acid structures.

The performance of the Cornell et al. AMBER force field generally corresponds well with and augments experimental data, but one notable exception could be the capping loops of double-helical stems. In addition, the performance of pairwise additive force fields is obviously unsatisfactory for inclusion of divalent cations, because their interactions lead to major polarization and charge-transfer effects neglected by the force field. Neglect of polarization also limits, albeit to a lesser extent, the description accuracy of other contributions, such as interactions with monovalent ions, conformational flexibility of the anionic sugar−phosphate backbone, hydrogen bonding, and solute polarization by solvent. Still, despite limitations, MD simulations are a valid tool for analyzing the structural dynamics of existing experimental structures. Careful analysis of MD simulations can identify problematic aspects of an experimental RNA structure, unveil structural characteristics masked by experimental constraints, reveal functionally significant stochastic fluctuations, evaluate the structural role of base ionization, and predict structurally and potentially functionally important details of the solvent behavior, including the presence of tightly bound water molecules. Moreover, combining classical MD simulations with QM calculations in hybrid QM/MM approaches helps in the assessment of the plausibility of chemical mechanisms of catalytic RNAs (ribozymes).

In contrast, the reliable prediction of structure from sequence information is beyond the applicability of MD tools. The ultimate utility of computational studies in understanding RNA function thus requires that the results are neither blindly accepted nor flatly rejected, but rather considered in the context of all available experimental data, with great care given to assessing limitations through the available starting structures, force field approximations, and sampling limitations. The examples given in this Account showcase how the judicious use of basic MD simulations has already served as a powerful tool to help evaluate the role of structural dynamics in biological function of RNA.

Direct Raman measurement of an elevated base pKa in the active site of a small ribozyme in a precatalytic conformation

Catalytic RNA molecules can achieve rate acceleration by shifting base pK(a) values toward neutrality. Prior evidence has suggested that base A38 of the hairpin ribozyme plays an important role in phosphoryl transfer, possibly functioning as a general acid, or by orienting a specific water molecule for proton transfer. To address the role of A38, we used Raman spectroscopy to measure directly the pK(a) of the N1-imino moiety in the context of hairpin ribozyme crystals representative of a "precatalytic" conformation. The results revealed that the pK(a) of A38 is shifted to 5.46 +/- 0.05 relative to 3.68 +/- 0.06 derived from a reference solution of the nucleotide AMP. The elevated pK(a) correlates well with the first titration point of the macroscopic pH-rate profile of the hairpin ribozyme in solution and strongly supports A38 as a general acid catalyst in bond scission. The results confirm that A38 is protonated before the transition state, which would promote phosphorane development. Overall, the results establish a cogent structure-function paradigm that expands our understanding of how RNA structure can enhance nucleobase reactivity to catalyze biological reactions.

Scope and mechanism of interstrand cross-link formation by the C4'-oxidized abasic site

The C4'-oxidized abasic site (C4-AP) is a commonly formed DNA lesion, which generates two types of interstrand cross-links (ICLs). The kinetically favored cross-link consists of two full length strands and forms reversibly and exclusively with dA. Cross-link formation is attributed to condensation of C4-AP with the N6-amino group of dA. Formation of the thermodynamic ICL involves cleavage of the strand containing C4-AP on the 3'-side of the lesion. The ratios and yields of the ICLs are highly dependent upon the local sequence. Product analysis of enzyme-digested material reveals that the ICL with dA is a cyclic adduct. Formation of the thermodynamically favored cross-link is catalyzed by the surrounding DNA sequence and occurs favorably with dC and dA but not with dG or dT. Mechanistic studies indicate that beta-elimination from C4-AP is the rate-limiting step in the formation of the thermodynamic ICL and that the local DNA environment determines the rate constant for this reaction. The efficiency of ICL formation, the stability of the thermodynamic products, and their possible formation in cells (Regelus, P.; et al. Proc. Natl. Acad. Sci. U.S.A. 2007, 104, 14032) suggest that these lesions will be deleterious to the biological system in which they are produced.

Theoretical studies of RNA catalysis: hybrid QM/MM methods and their comparison with MD and QM

Hybrid QM/MM methods combine the rigor of quantum mechanical (QM) calculations with the low computational cost of empirical molecular mechanical (MM) treatment allowing to capture dynamic properties to probe critical atomistic details of enzyme reactions. Catalysis by RNA enzymes (ribozymes) has only recently begun to be addressed with QM/MM approaches and is thus still a field under development. This review surveys methodology as well as recent advances in QM/MM applications to RNA mechanisms, including those of the HDV, hairpin, and hammerhead ribozymes, as well as the ribosome. We compare and correlate QM/MM results with those from QM and/or molecular dynamics (MD) simulations, and discuss scope and limitations with a critical eye on current shortcomings in available methodologies and computer resources. We thus hope to foster mutual appreciation and facilitate collaboration between experimentalists and theorists to jointly advance our understanding of RNA catalysis at an atomistic level.

Exploring ribozyme conformational changes with X-ray crystallography

Relating three-dimensional fold to function is a central challenge in RNA structural biology. Toward this goal, X-ray crystallography has long been considered the "gold standard" for structure determinations at atomic resolution, although NMR spectroscopy has become a powerhouse in this arena as well. In the area of dynamics, NMR remains the dominant technique to probe the magnitude and timescales of molecular motion. Although the latter area remains largely unassailable by conventional crystallographic methods, inroads have been made on proteins using Laue radiation on timescales of ms to ns. Proposed 'fourth generation' radiation sources, such as free-electron X-ray lasers, promise ps- to fs-timescale resolution, and credible evidence is emerging that supports the feasibility of single molecule imaging. At present however, the preponderance of RNA structural information has been derived from timescale and motion insensitive crystallographic techniques. Importantly, developments in computing, automation and high-flux synchrotron sources have propelled the rapidity of 'conventional' RNA crystal structure determinations to timeframes of hours once a suitable set of phases is obtained. With a sufficient number of crystal structures, it is possible to create a structural ensemble that can provide insight into global and local molecular motion characteristics that are relevant to biological function. Here we describe techniques to explore conformational changes in the hairpin ribozyme, a representative non-protein-coding RNA catalyst. The approaches discussed include: (i) construct choice and design using prior knowledge to improve X-ray diffraction; (ii) recognition of long-range conformational changes and (iii) use of single-base or single-atom changes to create ensembles. The methods are broadly applicable to other RNA systems.

Identification of an imino group indispensable for cleavage by a small ribozyme

The hairpin ribozyme is a small, noncoding RNA (ncRNA) that catalyzes a site-specific phosphodiester bond cleavage reaction. Prior biochemical and structural analyses pinpointed the amidine moiety of base Ade38 as a key functional group in catalysis, but base changes designed to probe function resulted in localized misfolding of the active site. To define the requirements for chemical activity using a conservative modification, we synthesized and incorporated N1-deazaadenosine into the full-length ribozyme construct. This single-atom variant severely impairs activity, although the active-site fold remains intact in the accompanying crystal structures. The results demonstrate the essentiality of the imino moiety as well as the importance of its interaction with the substrate in the precatalytic and transition-state conformations. This work demonstrates the efficacy of single-atom approaches in the analysis of ncRNA structure-function relationships.

Classification and energetics of the base-phosphate interactions in RNA

Structured RNA molecules form complex 3D architectures stabilized by multiple interactions involving the nucleotide base, sugar and phosphate moieties. A significant percentage of the bases in structured RNA molecules in the Protein Data Bank (PDB) hydrogen-bond with phosphates of other nucleotides. By extracting and superimposing base-phosphate (BPh) interactions from a reduced-redundancy subset of 3D structures from the PDB, we identified recurrent phosphate-binding sites on the RNA bases. Quantum chemical calculations were carried out on model systems representing each BPh interaction. The calculations show that the centers of each cluster obtained from the structure superpositions correspond to energy minima on the potential energy hypersurface. The calculations also show that the most stable phosphate-binding sites occur on the Watson-Crick edge of guanine and the Hoogsteen edge of cytosine. We modified the 'Find RNA 3D' (FR3D) software suite to automatically find and classify BPh interactions. Comparison of the 3D structures of the 16S and 23S rRNAs of Escherichia coli and Thermus thermophilus revealed that most BPh interactions are phylogenetically conserved and they occur primarily in hairpin, internal or junction loops or as part of tertiary interactions. Bases that form BPh interactions, which are conserved in the rRNA 3D structures are also conserved in homologous rRNA sequence alignments.