General acid-base catalysis in the mechanism of a hepatitis delta virus ribozyme

pKa shifting in double-stranded RNA is highly dependent upon nearest neighbors and bulge positioning

Shifting of pKa's in RNA is important for many biological processes; however, the driving forces responsible for shifting are not well understood. Herein, we determine how structural environments surrounding protonated bases affect pKa shifting in double-stranded RNA (dsRNA). Using (31)P NMR, we determined the pKa of the adenine in an A(+)·C base pair in various sequence and structural environments. We found a significant dependence of pKa on the base pairing strength of nearest neighbors and the location of a nearby bulge. Increasing nearest neighbor base pairing strength shifted the pKa of the adenine in an A(+)·C base pair higher by an additional 1.6 pKa units, from 6.5 to 8.1, which is well above neutrality. The addition of a bulge two base pairs away from a protonated A(+)·C base pair shifted the pKa by only ~0.5 units less than a perfectly base paired hairpin; however, positioning the bulge just one base pair away from the A(+)·C base pair prohibited formation of the protonated base pair as well as several flanking base pairs. Comparison of data collected at 25 °C and 100 mM KCl to biological temperature and Mg(2+) concentration revealed only slight pKa changes, suggesting that similar sequence contexts in biological systems have the potential to be protonated at biological pH. We present a general model to aid in the determination of the roles protonated bases may play in various dsRNA-mediated processes including ADAR editing, miRNA processing, programmed ribosomal frameshifting, and general acid-base catalysis in ribozymes.

Thio effects and an unconventional metal ion rescue in the genomic hepatitis delta virus ribozyme

Metal ion and nucleobase catalysis are important for ribozyme mechanism, but the extent to which they cooperate is unclear. A crystal structure of the hepatitis delta virus (HDV) ribozyme suggested that the pro-RP oxygen at the scissile phosphate directly coordinates a catalytic Mg(2+) ion and is within hydrogen bonding distance of the amine of the general acid C75. Prior studies of the genomic HDV ribozyme, however, showed neither a thio effect nor metal ion rescue using Mn(2+). Here, we combine experiment and theory to explore phosphorothioate substitutions at the scissile phosphate. We report significant thio effects at the scissile phosphate and metal ion rescue with Cd(2+). Reaction profiles with an SP-phosphorothioate substitution are indistinguishable from those of the unmodified substrate in the presence of Mg(2+) or Cd(2+), supporting the idea that the pro-SP oxygen does not coordinate metal ions. The RP-phosphorothioate substitution, however, exhibits biphasic kinetics, with the fast-reacting phase displaying a thio effect of up to 5-fold and the slow-reacting phase displaying a thio effect of ~1000-fold. Moreover, the fast- and slow-reacting phases give metal ion rescues in Cd(2+) of up to 10- and 330-fold, respectively. The metal ion rescues are unconventional in that they arise from Cd(2+) inhibiting the oxo substrate but not the RP substrate. This metal ion rescue suggests a direct interaction of the catalytic metal ion with the pro-RP oxygen, in line with experiments with the antigenomic HDV ribozyme. Experiments without divalent ions, with a double mutant that interferes with Mg(2+) binding, or with C75 deleted suggest that the pro-RP oxygen plays at most a redundant role in positioning C75. Quantum mechanical/molecular mechanical (QM/MM) studies indicate that the metal ion contributes to catalysis by interacting with both the pro-RP oxygen and the nucleophilic 2'-hydroxyl, supporting the experimental findings.

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.

Identification of the catalytic Mg²⁺ ion in the hepatitis delta virus ribozyme

The hepatitis delta virus ribozyme catalyzes an RNA cleavage reaction using a catalytic nucleobase and a divalent metal ion. The catalytic base, C75, serves as a general acid and has a pK(a) shifted toward neutrality. Less is known about the role of metal ions in the mechanism. A recent crystal structure of the precleavage ribozyme identified a Mg²⁺ ion that interacts through its partial hydration sphere with the G25·U20 reverse wobble. In addition, this Mg²⁺ ion is in position to directly coordinate the nucleophile, the 2'-hydroxyl of U(-1), suggesting it can serve as a Lewis acid to facilitate deprotonation of the 2'-hydroxyl. To test the role of the active site Mg²⁺ ion, we replaced the G25·U20 reverse wobble with an isosteric A25·C20 reverse wobble. This change was found to significantly reduce the negative potential at the active site, as supported by electrostatics calculations, suggesting that active site Mg²⁺ binding could be adversely affected by the mutation. The kinetic analysis and molecular dynamics of the A25·C20 double mutant suggest that this variant stably folds into an active structure. However, pH-rate profiles of the double mutant in the presence of Mg²⁺ are inverted relative to the profiles for the wild-type ribozyme, suggesting that the A25·C20 double mutant has lost the active site metal ion. Overall, these studies support a model in which the partially hydrated Mg²⁺ positioned at the G25·U20 reverse wobble is catalytic and could serve as a Lewis acid, a Brønsted base, or both to facilitate deprotonation of the nucleophile.

Role of pK(a) of nucleobases in the origins of chemical evolution

The formation of canonical base pairs through Watson-Crick hydrogen bonding sits at the heart of the genetic apparatus. The specificity of the base pairing of adenine with thymine/uracil and guanine with cytosine preserves accurate information for the biochemical blueprint and replicates the instructions necessary for carrying out biological function. The chemical evolution question of how these five canonical nucleobases were selected over various other possibilities remains intriguing. Since these and alternative nucleobases would have been available for chemical evolution, the reasons for the emergence of this system appear to be primarily functional. While investigating the base-pairing properties of structural nucleic acid analogs, we encountered a relationship between the pK(a) of a series of nonstandard (and canonical) nucleobases and the pH of the aqueous medium. This relationship appeared to correspond with the propensity of these molecules to self-assemble via Watson-Crick-type base-pairing interactions. A simple correlation of the "magnitude of the difference between the pK(a) and pH" (pK(a)-pH correlation) enables a general prediction of which types of heterocyclic recognition elements form hydrogen-bonded base pairs in aqueous media. Using the pK(a)-pH relationship, we can rationalize why nature chose the canonical nucleobases in terms of hydrophobic and hydrophilic interactions, and further extrapolate its significance within the context of chemical evolution. The connection between the physicochemical properties of bioorganic compounds and the interactions with their aqueous environment directly affects structure and function, at both a molecular and a supramolecular level. A general structure-function pattern emerges in biomolecules and biopolymers in aqueous media near neutral pH. A pK(a) - pH < 2 generally prompts catalytic functions, central to metabolism, but a difference in pK(a) - pH > 2 seems to result in the emergence of structure, central to replication. While this general trend is observed throughout extant biology, it could have also been an important factor in chemical evolution.

The glmS ribozyme cofactor is a general acid-base catalyst

The glmS ribozyme is the first natural self-cleaving ribozyme known to require a cofactor. The d-glucosamine-6-phosphate (GlcN6P) cofactor has been proposed to serve as a general acid, but its role in the catalytic mechanism has not been established conclusively. We surveyed GlcN6P-like molecules for their ability to support self-cleavage of the glmS ribozyme and found a strong correlation between the pH dependence of the cleavage reaction and the intrinsic acidity of the cofactors. For cofactors with low binding affinities, the contribution to rate enhancement was proportional to their intrinsic acidity. This linear free-energy relationship between cofactor efficiency and acid dissociation constants is consistent with a mechanism in which the cofactors participate directly in the reaction as general acid-base catalysts. A high value for the Brønsted coefficient (β ~ 0.7) indicates that a significant amount of proton transfer has already occurred in the transition state. The glmS ribozyme is the first self-cleaving RNA to use an exogenous acid-base catalyst.

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.

Insights into aquatic toxicities of the antibiotics oxytetracycline and ciprofloxacin in the presence of metal: complexation versus mixture

Co-contamination of ligand-like antibiotics (e.g., tetracyclines and quinolones) and heavy metals prevails in the environment, and thus the complexation between them is involved in environmental risks of antibiotics. To understand toxicological significance of the complex, effects of metal coordination on antibiotics' toxicity were investigated. The complexation of two antibiotics, oxytetracycline and ciprofloxacin, with three heavy metals, copper, zinc, and cadmium, was verified by spectroscopic techniques. The antibiotics bound metals via multiple coordination sites and rendered a mixture of various complexation speciations. Toxicity analysis indicated that metal coordination did modify the toxicity of the antibiotics and that antibiotic, metal, and their complex acted primarily as concentration addition. Comparison of EC(50) values revealed that the complex commonly was highest toxic and predominately correlated in toxicity to the mixture. Finally, environmental scenario analysis demonstrated that ignoring complexation would improperly classify environmental risks of the antibiotics.

Mechanisms of RNA catalysis

Ribozymes are RNA molecules that act as chemical catalysts. In contemporary cells, most known ribozymes carry out phosphoryl transfer reactions. The nucleolytic ribozymes comprise a class of five structurally-distinct species that bring about site-specific cleavage by nucleophilic attack of the 2'-O on the adjacent 3'-P to form a cyclic 2',3'-phosphate. In general, they will also catalyse the reverse reaction. As a class, all these ribozymes appear to use general acid-base catalysis to accelerate these reactions by about a million-fold. In the Varkud satellite ribozyme, we have shown that the cleavage reaction is catalysed by guanine and adenine nucleobases acting as general base and acid, respectively. The hairpin ribozyme most probably uses a closely similar mechanism. Guanine nucleobases appear to be a common choice of general base, but the general acid is more variable. By contrast, the larger ribozymes such as the self-splicing introns and RNase P act as metalloenzymes.

Mechanistic analysis of the hepatitis delta virus (HDV) ribozyme: methods for RNA preparation, structure mapping, solvent isotope effects, and co-transcriptional cleavage

Small ribozymes such as the hairpin, hammerhead, VS, glm S, and hepatitis delta virus (HDV) are self-cleaving RNAs that are typically characterized by kinetics and structural methods. Working with these RNAs requires attention to numerous experimental details. In this chapter we focus on four different experimental aspects of ribozyme studies: preparing the RNA, mapping its structure with reverse transcription and end-labeled techniques, solvent isotope experiments, and co-transcriptional cleavage assays. Although the focus of these methods is the HDV ribozyme, the methods should be applicable to other ribozymes.

Charged nucleobases and their potential for RNA catalysis

Catalysis in living cells is carried out by both proteins and RNA. Protein enzymes have been known for over 200 years, but RNA enzymes, or "ribozymes", were discovered only 30 years ago. Developing insight into RNA enzyme mechanisms is invaluable for better understanding both extant biological catalysis as well as the primitive catalysis envisioned in an early RNA-catalyzed life. Natural ribozymes include large RNAs such as the group I and II introns; small RNAs such as the hepatitis delta virus and the hairpin, hammerhead, VS, and glmS ribozymes; and the RNA portion of the ribosome and spliceosome. RNA enzymes use many of the same catalytic strategies as protein enzymes, but do so with much simpler side chains. Among these strategies are metal ion, general acid-base, and electrostatic catalysis. In this Account, we examine evidence for participation of charged nucleobases in RNA catalysis. Our overall approach is to integrate direct measurements on catalytic RNAs with thermodynamic studies on oligonucleotide model systems. The charged amino acids make critical contributions to the mechanisms of nearly all protein enzymes. Ionized nucleobases should be critical for RNA catalysis as well. Indeed, charged nucleobases have been implicated in RNA catalysis as general acid-bases and oxyanion holes. We provide an overview of ribozyme studies involving nucleobase catalysis and the complications involved in developing these mechanisms. We also consider driving forces for perturbation of the pK(a) values of the bases. Mechanisms for pK(a) values shifting toward neutrality involve electrostatic stabilization and the addition of hydrogen bonding. Both mechanisms couple protonation with RNA folding, which we treat with a thermodynamic formalism and conceptual models. Furthermore, ribozyme reaction mechanisms can be multichannel, which demonstrates the versatility of ribozymes but makes analysis of experimental data challenging. We examine advances in measuring and analyzing perturbed pK(a) values in RNA. Raman crystallography and fluorescence spectroscopy have been especially important for pK(a) measurement. These methods reveal pK(a) values for the nucleobases A or C equal to or greater than neutrality, conferring potential histidine- and lysine/arginine-like behavior on them. Structural support for ionization of the nucleobases also exists: an analysis of RNA structures in the databases conducted herein suggests that charging of the bases is neither especially uncommon nor difficult to achieve under cellular conditions. Our major conclusions are that cationic and anionic charge states of the nucleobases occur in RNA enzymes and that these states make important catalytic contributions to ribozyme activity. We conclude by considering outstanding questions and possible experimental and theoretical approaches for further advances.

The molecular interactions that stabilize RNA tertiary structure: RNA motifs, patterns, and networks

RNA molecules adopt specific three-dimensional structures critical to their function. Many essential metabolic processes, including protein synthesis and RNA splicing, are carried out by RNA molecules with elaborate tertiary structures (e.g. 3QIQ, right). Indeed, the ribosome and self-splicing introns are complex RNA machines. But even the coding regions in messenger RNAs and viral RNAs are flanked by highly structured untranslated regions, which provide regulatory information necessary for gene expression. RNA tertiary structure is defined as the three-dimensional arrangement of RNA building blocks, which include helical duplexes, triple-stranded structures, and other components that are held together through connections collectively termed RNA tertiary interactions. The structural diversity of these interactions is now a subject of intense investigation, involving the techniques of NMR, X-ray crystallography, chemical genetics, and phylogenetic analysis. At the same time, many investigators are using biophysical techniques to elucidate the driving forces for tertiary structure formation and the mechanisms for its stabilization. RNA tertiary folding is promoted by maximization of base stacking, much like the hydrophobic effect that drives protein folding. RNA folding also requires electrostatic stabilization, both through charge screening and site binding of metals, and it is enhanced by desolvation of the phosphate backbone. In this Account, we provide an overview of the features that specify and stabilize RNA tertiary structure. A major determinant for overall tertiary RNA architecture is local conformation in secondary-structure junctions, which are regions from which two or more duplexes project. At junctions and other structures, such as pseudoknots and kissing loops, adjacent helices stack on one another, and these coaxial stacks play a major role in dictating the overall architectural form of an RNA molecule. In addition to RNA junction topology, a second determinant for RNA tertiary structure is the formation of sequence-specific interactions. Networks of triple helices, tetraloop-receptor interactions, and other sequence-specific contacts establish the framework for the overall tertiary fold. The third determinant of tertiary structure is the formation of stabilizing stacking and backbone interactions, and many are not sequence specific. For example, ribose zippers allow 2'-hydroxyl groups on different RNA strands to form networks of interdigitated hydrogen bonds, serving to seal strands together and thereby stabilize adjacent substructures. These motifs often require monovalent and divalent cations, which can interact diffusely or through chelation to specific RNA functional groups. As we learn more about the components of RNA tertiary structure, we will be able to predict the structures of RNA molecules from their sequences, thereby obtaining key information about biological function. Understanding and predicting RNA structure is particularly important given the recent discovery that although most of our genome is transcribed into RNA molecules, few of them have a known function. The prevalence of RNA viruses and pathogens with RNA genomes makes RNA drug discovery an active area of research. Finally, knowledge of RNA structure will facilitate the engineering of supramolecular RNA structures, which can be used as nanomechanical components for new materials. But all of this promise depends on a better understanding of the RNA parts list, and how the pieces fit together.

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.

Quantum Mechanical/Molecular Mechanical Study of the HDV Ribozyme: Impact of the Catalytic Metal Ion on the Mechanism

A recent crystal structure of the precleaved HDV ribozyme along with biochemical data support a mechanism for phosphodiester bond self-cleavage in which C75 acts as a general acid and bound Mg(2+) ion acts as a Lewis acid. Herein this precleaved crystal structure is used as the basis for quantum mechanical/molecular mechanical calculations. These calculations indicate that the self-cleavage reaction is concerted with a phosphorane-like transition state when a divalent ion, Mg(2+) or Ca(2+), is bound at the catalytic site but is sequential with a phosphorane intermediate when a monovalent ion, such as Na(+), is at this site. Electrostatic potential calculations suggest that the divalent metal ion at the catalytic site lowers the pK(a) of C75, leading to the concerted mechanism in which the proton is partially transferred to the leaving group in the phosphorane-like transition state. These observations are consistent with experimental data, including pK(a) measurements, reaction kinetics, and proton inventories with divalent and monovalent ions.

Characterization of the Structure and Dynamics of the HDV Ribozyme at Different Stages Along the Reaction Path

The structure and dynamics of the hepatitis delta virus ribozyme (HDVr) are studies using molecular dynamics simulations at several stages along its catalytic reaction path, including reactant, activated precursor, transition state mimic and product states, departing from an initial structure based on the C75U mutant crystal structure (PDB: 1VC7). Results of five 350 ns molecular dynamics simulations reveal a spontaneous rotation of U-1 that leads to an in-line conformation and support the role of protonated C75 as the general acid in the transition state. Our results provide rationale for the interpretation of several important experimental results, and make experimentally testable predictions regarding the roles of key active site residues that are not obvious from any available crystal structures.

Two distinct catalytic strategies in the hepatitis δ virus ribozyme cleavage reaction

The hepatitis delta virus (HDV) ribozyme and related RNAs are widely dispersed in nature. This RNA is a small nucleolytic ribozyme that self-cleaves to generate products with a 2',3'-cyclic phosphate and a free 5'-hydroxyl. Although small ribozymes are dependent on divalent metal ions under biologically relevant buffer conditions, they function in the absence of divalent metal ions at high ionic strengths. This characteristic suggests that a functional group within the covalent structure of small ribozymes is facilitating catalysis. Structural and mechanistic analyses have demonstrated that the HDV ribozyme active site contains a cytosine with a perturbed pK(a) that serves as a general acid to protonate the leaving group. The reaction of the HDV ribozyme in monovalent cations alone never approaches the velocity of the Mg(2+)-dependent reaction, and there is significant biochemical evidence that a Mg(2+) ion participates directly in catalysis. A recent crystal structure of the HDV ribozyme revealed that there is a metal binding pocket in the HDV ribozyme active site. Modeling of the cleavage site into the structure suggested that this metal ion can interact directly with the scissile phosphate and the nucleophile. In this manner, the Mg(2+) ion can serve as a Lewis acid, facilitating deprotonation of the nucleophile and stabilizing the conformation of the cleavage site for in-line attack of the nucleophile at the scissile phosphate. This catalytic strategy had previously been observed only in much larger ribozymes. Thus, in contrast to most large and small ribozymes, the HDV ribozyme uses two distinct catalytic strategies in its cleavage reaction.

The glmS ribozyme tunes the catalytically critical pK(a) of its coenzyme glucosamine-6-phosphate

The glmS ribozyme riboswitch is the first known natural catalytic RNA that employs a small-molecule cofactor. Binding of glucosamine-6-phosphate (GlcN6P) uncovers the latent self-cleavage activity of the RNA, which adopts a catalytically competent conformation that is nonetheless inactive in the absence of GlcN6P. Structural and analogue studies suggest that the amine of GlcN6P functions as a general acid-base catalyst, while its phosphate is important for binding affinity. However, the solution pK(a) of the amine is 8.06 ± 0.05, which is not optimal for proton transfer. Here we used Raman crystallography directly to determine the pK(a)'s of GlcN6P bound to the glmS ribozyme. Binding to the RNA lowers the pK(a) of the amine of GlcN6P to 7.26 ± 0.09 and raises the pK(a) of its phosphate to 6.35 ± 0.09. Remarkably, the pK(a)'s of these two functional groups are unchanged from their values for free GlcN6P (8.06 ± 0.05 and 5.98 ± 0.05, respectively) when GlcN6P binds to the catalytically inactive but structurally unperturbed G40A mutant of the ribozyme, thus implicating the ribozyme active site guanine in pK(a) tuning. This is the first demonstration that a ribozyme can tune the pK(a) of a small-molecule ligand. Moreover, the anionic glmS ribozyme in effect stabilizes the neutral amine of GlcN6P by lowering its pK(a). This is unprecedented and illustrates the chemical sophistication of ribozyme active sites.

Modeling the RNA 2'OH activation: possible roles of metal ion and nucleobase as catalysts in self-cleaving ribozymes

The RNA 2'OH activation as taking place in the first chemical step of self-cleaving ribozymes is studied theoretically by DFT and MP2 methods using a continuum solvation model (CPCM). The reaction of proton transfer is studied in the presence of two kinds of catalysts: a fully hydrated metal ion (Mg(2+)) or partially hydrated nucleobase (guanine), taken separately or together leading to three different modes of activation. The metal ion is either directly bound (inner-sphere) or indirectly bound (outer-sphere) to the 2'OH group and a hydroxide ion acts as a general or specific base; the nucleobase is taken in anionic or in neutral enol-tautomeric forms playing itself the role of general base. The presence of a close metal ion (outer-sphere) lowers the pK(a) value of the 2'OH group by several log units in both metal-ion and nuleobase catalysis. The direct metal coordination to the 2'OH group (inner-sphere) further stabilizes the developing negative charge on the nucleophile. The switching from the inner-sphere to the outer-sphere coordination appears to be driven by the energy cost for reorganizing the first coordination shell rather than by the electrostatic repulsion between the ligands. The metal-ion catalysis is more effective with a specific base in the dianionic mechanism. On the other hand, the nucleobase catalysis is more effective in the monoanionic mechanism and in the presence of a metal ion acting as a cofactor through nonspecific electrostatic interactions. The results establish a baseline to study the possible roles of metal and nucleobase catalysts and their environment in more realistic models for self-cleaving ribozymes.

Identification of an hepatitis delta virus-like ribozyme at the mRNA 5'-end of the L1Tc retrotransposon from Trypanosoma cruzi

L1Tc is a non-LTR LINE element from Trypanosoma cruzi that encodes its transposition machinery and bears an internal promoter. Herewith, we report the identification of an in vitro active hepatitis delta virus-like ribozyme located in the first 77 nt at the 5′-end of the L1Tc mRNA (L1TcRz). The data presented show that L1TcRz has a co-transcriptional function. Using gel-purified uncleaved RNA transcripts, the data presented indicate that the kinetics of the self-cleaving, in a magnesium-dependent reaction, fits to a two-phase decay curve. The cleavage point identified by primer extension takes place at +1 position of the element. The hydroxyl nature of the 5′-end of the 3′-fragment generated by the cleavage activity of L1TcRz was confirmed. Since we have previously described that the 77-nt long fragment located at the 5′-end of L1Tc has promoter activity, the existence of a ribozyme in L1Tc makes this element to be the first described non-LTR retroelement that has an internal promoter–ribozyme dual function. The L1Tc nucleotides located downstream of the ribozyme catalytic motif appear to inhibit its activity. This inhibition may be influenced by the existence of a specific L1Tc RNA conformation that is recognized by RNase P.

Mechanistic strategies in the HDV ribozyme: chelated and diffuse metal ion interactions and active site protonation

The crystal structure of the precleaved form of the hepatitis delta virus (HDV) ribozyme reveals two G•U wobbles near the active site: a rare reverse G•U wobble involving a syn G base, and a standard G•U wobble at the cleavage site. The catalytic mechanism for this ribozyme has been proposed to involve a Mg(2+) ion bound to the reverse G•U wobble, as well as a protonated C75 base. We carried out molecular dynamics simulations to analyze metal ion interaction with the reverse and standard G•U wobbles and to investigate the impact of C75 protonation on the structure and motions of the ribozyme. We identified two types of Mg(2+) ions associated with the ribozyme, chelated and diffuse, at the reverse and standard G•U wobbles, respectively, which appear to contribute to catalysis and stability, respectively. These two metal ion sites exhibit relatively independent behavior. Protonation of C75 was observed to locally organize the active site in a manner that facilitates the catalytic mechanism, in which C75(+) acts as a general acid and Mg(2+) as a Lewis acid. The simulations also indicated that the overall structure and thermal motions of the ribozyme are not significantly influenced by the catalytic Mg(2+) interaction or C75 protonation. This analysis suggests that the reaction pathway of the ribozyme is dominated by small local motions at the active site rather than large-scale global conformational changes. These results are consistent with a wealth of experimental data.

Base ionization and ligand binding: how small ribozymes and riboswitches gain a foothold in a protein world

Genome sequencing has produced thousands of nonprotein coding (nc)RNA sequences including new ribozymes and riboswitches. Such RNAs are notable for their extraordinary functionality, which entails exquisite folding that culminates in biocatalytic or ligand-binding capabilities. Here we discuss advances in relating ncRNA form to function with an emphasis on base pK(a) shifting by the hairpin and hepatitis delta virus ribozymes. We then describe ligand binding by the two smallest riboswitches, which target preQ(1) and S-adenosyl-(l)-homocysteine, followed by an analysis of a second-messenger riboswitch that binds cyclic-di-GMP. Each riboswitch is then compared to a protein that binds the same ligand to contrast binding properties. The results showcase the breadth of functionality attainable from ncRNAs, as well as molecular features notable for antibacterial design.

The N-terminal nucleophile serine of cephalosporin acylase executes the second autoproteolytic cleavage and acylpeptide hydrolysis

Cephalosporin acylase (CA) precursor is translated as a single polypeptide chain and folds into a self-activating pre-protein. Activation requires two peptide bond cleavages that excise an internal spacer to form the mature αβ heterodimer. Using Q-TOF LC-MS, we located the second cleavage site between Glu(159) and Gly(160), and detected the corresponding 10-aa spacer (160)GDPPDLADQG(169) of CA mutants. The site of the second cleavage depended on Glu(159): moving Glu into the spacer or removing 5-10 residues from the spacer sequence resulted in shorter spacers with the cleavage at the carboxylic side of Glu. The mutant E159D was cleaved more slowly than the wild-type, as were mutants G160A and G160L. This allowed kinetic measurements showing that the second cleavage reaction was a first-order, intra-molecular process. Glutaryl-7-aminocephalosporanic acid is the classic substrate of CA, in which the N-terminal Ser(170) of the β-subunit, is the nucleophile. Glu and Asp resemble glutaryl, suggesting that CA might also remove N-terminal Glu or Asp from peptides. This was indeed the case, suggesting that the N-terminal nucleophile also performed the second proteolytic cleavage. We also found that CA is an acylpeptide hydrolase rather than a previously expected acylamino acid acylase. It only exhibited exopeptidase activity for the hydrolysis of an externally added peptide, supporting the intra-molecular interaction. We propose that the final CA activation is an intra-molecular process performed by an N-terminal nucleophile, during which large conformational changes in the α-subunit C-terminal region are required to bridge the gap between Glu(159) and Ser(170).

Metal binding motif in the active site of the HDV ribozyme binds divalent and monovalent ions

The hepatitis delta virus (HDV) ribozyme uses both metal ion and nucleobase catalysis in its cleavage mechanism. A reverse G·U wobble was observed in a recent crystal structure of the precleaved state. This unusual base pair positions a Mg(2+) ion to participate in catalysis. Herein, we used molecular dynamics (MD) and X-ray crystallography to characterize the conformation and metal binding characteristics of this base pair in product and precleaved forms. Beginning with a crystal structure of the product form, we observed formation of the reverse G·U wobble during MD trajectories. We also demonstrated that this base pair is compatible with the diffraction data for the product-bound state. During MD trajectories of the product form, Na(+) ions interacted with the reverse G·U wobble in the RNA active site, and a Mg(2+) ion, introduced in certain trajectories, remained bound at this site. Beginning with a crystal structure of the precleaved form, the reverse G·U wobble with bound Mg(2+) remained intact during MD simulations. When we removed Mg(2+) from the starting precleaved structure, Na(+) ions interacted with the reverse G·U wobble. In support of the computational results, we observed competition between Na(+) and Mg(2+) in the precleaved ribozyme crystallographically. Nonlinear Poisson-Boltzmann calculations revealed a negatively charged patch near the reverse G·U wobble. This anionic pocket likely serves to bind metal ions and to help shift the pK(a) of the catalytic nucleobase, C75. Thus, the reverse G·U wobble motif serves to organize two catalytic elements, a metal ion and catalytic nucleobase, within the active site of the HDV ribozyme.

Catalysis by the nucleolytic ribozymes

The nucleolytic ribozymes use general acid-base catalysis to contribute significantly to their rate enhancement. The VS (Varkud satellite) ribozyme uses a guanine and an adenine nucleobase as general base and acid respectively in the cleavage reaction. The hairpin ribozyme is probably closely similar, while the remaining nucleolytic ribozymes provide some interesting contrasts.

Molecular recognition of 6'-N-5-hexynoate kanamycin A and RNA 1x1 internal loops containing CA mismatches

In our previous study to identify the RNA internal loops that bind an aminoglycoside derivative, we determined that 6'-N-5-hexynoate kanamycin A prefers to bind 1x1 nucleotide internal loops containing C·A mismatches. In this present study, the molecular recognition between a variety of RNAs that are mutated around the C·A loop and the ligand was investigated. Studies show that both loop nucleotides and loop closing pairs affect binding affinity. Most interestingly, it was shown that there is a correlation between the thermodynamic stability of the C·A internal loops and ligand affinity. Specifically, C·A loops that had relatively high or low stability bound the ligand most weakly whereas loops with intermediate stability bound the ligand most tightly. In contrast, there is no correlation between the likelihood that a loop forms a C-A(+) pair at lower pH and ligand affinity. It was also found that a 1x1 nucleotide C·A loop that bound to the ligand with the highest affinity is identical to the consensus site in RNAs that are edited by adenosine deaminases acting on RNA type 2 (ADAR2). These studies provide a detailed investigation of factors affecting small molecule recognition of internal loops containing C·A mismatches, which are present in a variety of RNAs that cause disease.

Do the hairpin and VS ribozymes share a common catalytic mechanism based on general acid-base catalysis? A critical assessment of available experimental data

The active centers of the hairpin and VS ribozymes are both generated by the interaction of two internal loops, and both ribozymes use guanine and adenine nucleobases to accelerate cleavage and ligation reactions. The centers are topologically equivalent and the relative positioning of key elements the same. There is good evidence that the cleavage reaction of the VS ribozyme is catalyzed by the guanine (G638) acting as general base and the adenine (A756) as general acid. We now critically evaluate the experimental mechanistic evidence for the hairpin ribozyme. We conclude that all the available data are fully consistent with a major contribution to catalysis by general acid-base catalysis involving the adenine (A38) and guanine (G8). It appears that the two ribozymes are mechanistically equivalent.

Protonation of base pairs in RNA: context analysis and quantum chemical investigations of their geometries and stabilities

Base pairs involving protonated nucleobases play important roles in mediating global macromolecular conformational changes and in facilitation of catalysis in a variety of functional RNA molecules. Here we present our attempts at understanding the role of such base pairs by detecting possible protonated base pairs in the available RNA crystal structures using BPFind software, in their specific structural contexts, and by the characterization of their geometries, interaction energies, and stabilities using advanced quantum chemical computations. We report occurrences of 18 distinct protonated base pair combinations from a representative data set of RNA crystal structures and propose a theoretical model for one putative base pair combination. Optimization of base pair geometries was carried out at the B3LYP/cc-pVTZ level, and the BSSE corrected interaction energies were calculated at the MP2/aug-cc-pVDZ level of theory. The geometries for each of the base pairs were characterized in terms of H-bonding patterns observed, rmsd values observed on optimization, and base pair geometrical parameters. In addition, the intermolecular interaction in these complexes was also analyzed using Morokuma energy decomposition. The gas phase interaction energies of the base pairs range from -24 to -49 kcal/mol and reveal the dominance of Hartree-Fock component of interaction energy constituting 73% to 98% of the total interaction energy values. On the basis of our combined bioinformatics and quantum chemical analysis of different protonated base pairs, we suggest resolution of structural ambiguities and correlate their geometric and energetic features with their structural and functional roles. In addition, we also examine the suitability of specific base pairs as key elements in molecular switches and as nucleators for higher order structures such as base triplets and quartets.

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.

The glmS ribozyme: use of a small molecule coenzyme by a gene-regulatory RNA

The glmS ribozyme is the first known example of a natural ribozyme that has evolved to require binding of an exogenous small molecule for activity. In Gram-positive bacteria, this RNA domain is part of the messenger RNA (mRNA) encoding the essential enzyme that synthesizes glucosamine-6-phosphate (GlcN6P). When present at physiologic concentration, this small molecule binds to the glmS ribozyme and uncovers a latent self-cleavage activity that ultimately leads to degradation of the mRNA. Biochemical and structural studies reveal that the RNA adopts a rigid fold stabilized by three pseudoknots and the packing of a peripheral domain against the ribozyme core. GlcN6P binding to this pre-organized RNA does not induce conformational changes; rather, the small molecule functions as a coenzyme, providing a catalytically essential amine group to the active site. The ribozyme is not a passive player, however. Active site functional groups are essential for catalysis, even in the presence of GlcN6P. In addition to being a superb experimental system with which to analyze how RNA catalysts can exploit small molecule coenzymes to broaden their chemical versatility, the presence of the glmS ribozyme in numerous pathogenic bacteria make this RNA an attractive target for the development of new antibiotics and antibacterial strategies.

Small self-cleaving ribozymes

The hammerhead, hairpin, hepatitis delta virus (HDV), Varkud Satellite (VS), and glmS ribozymes catalyze sequence-specific intramolecular cleavage of RNA. They range between 50 and 150 nucleotides in length, and are known as the "small self-cleaving ribozymes." Except for the glmS ribozyme that functions as a riboswitch in Gram-positive bacteria, they were originally discovered as domains of satellite RNAs. However, recent studies show that several of them are broadly distributed in genomes of organisms from many phyla. Each of these ribozymes has a unique overall architecture and active site organization. Crystal structures have revealed how RNA active sites can bind preferentially to the transition state of a reaction, whereas mechanistic studies have shown that nucleobases can efficiently perform general acid-base and electrostatic catalysis. This versatility explains the abundance of ribozymes in contemporary organisms and also supports a role for catalytic RNAs early in evolution.

Investigation of the catalytic mechanism of a synthetic DNAzyme with protein-like functionality: an RNaseA mimic?

The protein enzyme ribonuclease A (RNaseA) cleaves RNA with catalytic perfection, although with little sequence specificity, by a divalent metal ion (M(2+))-independent mechanism in which a pair of imidazoles provides general acid and base catalysis, while a cationic amine provides electrostatic stabilization of the transition state. Synthetic imitation of this remarkable organo-catalyst ("RNaseA mimicry") has been a longstanding goal in biomimetic chemistry. The 9(25)-11 DNAzyme contains synthetically modified nucleotides presenting both imidazole and cationic amine side chains, and catalyzes RNA cleavage with turnover in the absence of M(2+) similarly to RNaseA. Nevertheless, the catalytic roles, if any, of the "protein-like" functional groups have not been defined, and hence the question remains whether 9(25)-11 engages any of these functionalities to mimic aspects of the mechanism of RNaseA. To address this question, we report a mechanistic investigation of 9(25)-11 catalysis wherein we have employed a variety of experiments, such as DNAzyme functional group deletion, mechanism-based affinity labeling, and bridging and nonbridging phosphorothioate substitution of the scissile phosphate. Several striking parallels exist between the results presented here for 9(25)-11 and the results of analogous experiments applied previously to RNaseA. Specifically, our results implicate two particular imidazoles in general acid and base catalysis and suggest that a specific cationic amine stabilizes the transition state via diastereoselective interaction with the scissile phosphate. Overall, 9(25)-11 appears to meet the minimal criteria of an RNaseA mimic; this demonstrates how added synthetic functionality can expand the mechanistic repertoire available to a synthetic DNA-based catalyst.

Long-distance communication in the HDV ribozyme: insights from molecular dynamics and experiments

The hepatitis delta virus ribozyme is a small, self-cleaving RNA with a compact tertiary structure and buried active site that is important in the life cycle of the virus. The ribozyme's function in nature is to cleave an internal phosphodiester bond and linearize concatemers during rolling circle replication. Crystal structures of the ribozyme have been solved in both pre-cleaved and post-cleaved (product) forms and reveal an intricate network of interactions that conspire to catalyze bond cleavage. In addition, extensive biochemical studies have been performed to work out a mechanism for bond cleavage in which C75 and a magnesium ion catalyze the reaction by general acid-base chemistry. One issue that has remained unclear in this ribozyme and in other ribozymes is the nature of long-distance communication between peripheral regions of the RNA and the buried active site. We performed molecular dynamics simulations on the hepatitis delta virus ribozyme in the product form and assessed communication between a distal structural portion of the ribozyme-the protonated C41 base triple-and the active site containing the critical C75. We varied the ionization state of C41 in both the wild type and a C41 double mutant variant and determined the impact on the active site. In all four cases, effects at the active site observed in the simulations agree with experimental studies on ribozyme activity. Overall, these studies indicate that small functional RNAs have the potential to communicate interactions over long distances and that wild-type RNAs may have evolved ways to prevent such interactions from interfering with catalysis.

Cation-dependent cleavage of the duplex form of the subtype-B HIV-1 RNA dimerization initiation site

The crystal structure of subtype-B HIV-1 genomic RNA Dimerization Initiation Site duplex revealed chain cleavage at a specific position resulting in 3'-phosphate and 5'-hydroxyl termini. A crystallographic analysis showed that Ba(2+), Mn(2+), Co(2+) and Zn(2+) bind specifically on a guanine base close to the cleaved position. The crystal structures also point to a necessary conformational change to induce an 'in-line' geometry at the cleavage site. In solution, divalent cations increased the rate of cleavage with pH/pKa compensation, indicating that a cation-bound hydroxide anion is responsible for the cleavage. We propose a 'Trojan horse' mechanism, possibly of general interest, wherein a doubly charged cation hosted near the cleavage site as a 'harmless' species is further transformed in situ into an 'aggressive' species carrying a hydroxide anion.

The human HDV-like CPEB3 ribozyme is intrinsically fast-reacting

Self-cleaving RNAs have recently been identified in mammalian genomes. A small ribozyme related in structure to the hepatitis delta virus (HDV) ribozyme occurs in a number of mammals, including chimpanzees and humans, within an intron of the CPEB3 gene. The catalytic mechanisms for the CPEB3 and HDV ribozymes appear to be similar, generating cleavage products with 5'-hydroxyl and 2',3'-cyclic phosphate termini; nonetheless, the cleavage rate reported for the CPEB3 ribozyme is more than 6000-fold slower than for the fastest HDV ribozyme. Herein, we use full-length RNA and cotranscriptional self-cleavage assays to compare reaction rates among human CPEB3, chimp CPEB3, and HDV ribozymes. Our data reveal that a single base change of the upstream flanking sequence, which sequesters an intrinsically weak P1.1 pairing in a misfold, increases the rate of the wild-type human CPEB3 ribozyme by approximately 250-fold; thus, the human ribozyme is intrinsically fast-reacting. Secondary structure determination and native gel analyses reveal that the cleaved population of the CPEB3 ribozyme has a single, secondary structure that closely resembles the HDV ribozyme. In contrast, the precleavage population of the CPEB3 ribozyme appears to have a more diverse secondary structure, possibly reflecting misfolding with the upstream sequence and dynamics intrinsic to the ribozyme. Prior identification of expressed sequence tags (ESTs) in human cells indicated that cleavage activity of the human ribozyme is tissue-specific. It is therefore possible that cellular factors interact with regions upstream of the CPEB3 ribozyme to unmask its high intrinsic reactivity.

Nucleobase-mediated general acid-base catalysis in the Varkud satellite ribozyme

Existing evidence suggests that the Varkud satellite (VS) ribozyme accelerates the cleavage of a specific phosphodiester bond using general acid-base catalysis. The key functionalities are the nucleobases of adenine 756 in helix VI of the ribozyme, and guanine 638 in the substrate stem loop. This results in a bell-shaped dependence of reaction rate on pH, corresponding to groups with pK(a) = 5.2 and 8.4. However, it is not possible from those data to determine which nucleobase is the acid, and which the base. We have therefore made substrates in which the 5' oxygen of the scissile phosphate is replaced by sulfur. This labilizes the leaving group, removing the requirement for general acid catalysis. This substitution restores full activity to the highly impaired A756G ribozyme, consistent with general acid catalysis by A756 in the unmodified ribozyme. The pH dependence of the cleavage of the phosphorothiolate-modified substrates is consistent with general base catalysis by nucleobase at position 638. We conclude that cleavage of the substrate by the VS ribozyme is catalyzed by deprotonation of the 2'-O nucleophile by G638 and protonation of the 5'-O leaving group by A756.

Driving forces for nucleic acid pK(a) shifting in an A(+).C wobble: effects of helix position, temperature, and ionic strength

Secondary structure plays critical roles in nucleic acid function. Mismatches in DNA can lead to mutation and disease, and some mismatches involve a protonated base. Among protonated mismatches, A(+).C wobble pairs form near physiological pH and have relatively minor effects on helix geometry, making them especially important in biology. Herein, we investigate effects of helix position, temperature, and ionic strength on pK(a) shifting in A(+).C wobble pairs in DNA. We observe that pK(a) shifting is favored by internal A(+).C wobbles, which have low cooperativities of folding and make large contributions to stability, and disfavored by external A(+).C wobbles, which have high folding cooperativities but make small contributions to stability. An inverse relationship between pK(a) shifting and temperature is also found, which supports a model in which protonation is enthalpically favored overall and entropically correlated with cooperativity of folding. We also observe greater pK(a) shifts as the ionic strength decreases, consistent with anticooperativity between proton binding and counterion-condensed monovalent cation. Under the most favorable temperature and ionic strength conditions tested, a pK(a) of 8.0 is observed for the A(+).C wobble pair, which represents an especially large shift ( approximately 4.5 pK(a) units) from the unperturbed pK(a) value of adenosine. This study suggests that protonated A(+).C wobble pairs exist in DNA under biologically relevant conditions, where they can drive conformational changes and affect replication and transcription fidelity.

Competition between Co(NH(3)(6)3+ and inner sphere Mg2+ ions in the HDV ribozyme

Divalent cations play critical structural and functional roles in many RNAs. While the hepatitis delta virus (HDV) ribozyme can undergo self-cleavage in the presence of molar concentrations of monovalent cations, divalent cations such as Mg(2+) are required for efficient catalysis under physiological conditions. Moreover, the cleavage reaction can be inhibited with Co(NH(3))(6)(3+), an analogue of Mg(H(2)O)(6)(2+). Here, the binding of Mg(2+) and Co(NH(3))(6)(3+) to the HDV ribozyme is studied by Raman microscopic analysis of crystals. Raman difference spectra acquired at different metal ion conditions reveal changes in the ribozyme. When Mg(2+) alone is introduced to the ribozyme, inner sphere coordination of Mg(H(2)O)(x)(2+) (x </= 5) to nonbridging PO(2)(-) oxygen and changes in base stretches and phosphodiester group conformation are observed. In addition, binding of Mg(2+) induces deprotonation of a cytosine assigned to the general acid C75, consistent with solution studies. When Co(NH(3))(6)(3+) alone is introduced, deprotonation of C75 is again observed, as are distinctive changes in base vibrational ring modes and phosphodiester backbone conformation. In contrast to Mg(2+) binding, Co(NH(3))(6)(3+) binding does not perturb PO(2)(-) group vibrations, consistent with its ability to make only outer sphere contacts. Surprisingly, competitive binding studies reveal that Co(NH(3))(6)(3+) ions displace some inner sphere-coordinated magnesium species, including ions coordinated to PO(2)(-) groups or the N7 of a guanine, likely G1 at the active site. These observations contrast with the tenet that Co(NH(3))(6)(3+) ions displace only outer sphere magnesium ions. Overall, our data support two classes of inner sphere Mg(2+)-PO(2)(-) binding sites: sites that Co(NH(3))(6)(3+) can displace and others it cannot.

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.

Modulating RNA structure and catalysis: lessons from small cleaving ribozymes

RNA is a key molecule in life, and comprehending its structure/function relationships is a crucial step towards a more complete understanding of molecular biology. Even though most of the information required for their correct folding is contained in their primary sequences, we are as yet unable to accurately predict both the folding pathways and active tertiary structures of RNA species. Ribozymes are interesting molecules to study when addressing these questions because any modifications in their structures are often reflected in their catalytic properties. The recent progress in the study of the structures, the folding pathways and the modulation of the small ribozymes derived from natural, self-cleaving, RNA motifs have significantly contributed to today's knowledge in the field.

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.

Raman crystallography of RNA

Raman crystallography is the application of Raman spectroscopy to single crystals. This technique has been applied to a variety of protein molecules where it has provided unique information about biopolymer folding, substrate binding, and catalysis. Here, we describe the application of Raman crystallography to functional RNA molecules. RNA represents unique opportunities and challenges for Raman crystallography. One issue that confounds studies of RNA is its tendency to adopt multiple non-functional folds. Raman crystallography has the advantage that it isolates a single state of the RNA within the crystal and can evaluate its fold, metal ion binding properties (ligand identity, stoichiometry, and affinity), proton binding properties (identity, stoichiometry, and affinity), and catalytic potential. In particular, base-specific stretches can be identified and then associated with the binding of metal ions and protons. Because measurements are carried out in the hanging drop at ambient, rather than cryo, conditions and because RNA crystals tend to be approximately 70% solvent, RNA dynamics and conformational changes become experimentally accessible. This review focuses on experimental setup and procedures, acquisition and interpretation of Raman data, and determination of physicochemical properties of the RNA. Raman crystallographic and solution biochemical experiments on the HDV RNA enzyme are summarized and found to be in excellent agreement. Remarkably, characterization of the crystalline state has proven to help rather than hinder functional characterization of functional RNA, most likely because the tendency of RNA to fold heterogeneously is limited in a crystalline environment. Future applications of Raman crystallography to RNA are briefly discussed.