RNA catalysis

Role and Perspective of Molecular Simulation-Based Investigation of RNA-Ligand Interaction: From Small Molecules and Peptides to Photoswitchable RNA Binding

Aberrant RNA–protein complexes are formed in a variety of diseases. Identifying the ligands that interfere with their formation is a valuable therapeutic strategy. Molecular simulation, validated against experimental data, has recently emerged as a powerful tool to predict both the pose and energetics of such ligands. Thus, the use of molecular simulation may provide insight into aberrant molecular interactions in diseases and, from a drug design perspective, may allow for the employment of less wet lab resources than traditional in vitro compound screening approaches. With regard to basic research questions, molecular simulation can support the understanding of the exact molecular interaction and binding mode. Here, we focus on examples targeting RNA–protein complexes in neurodegenerative diseases and viral infections. These examples illustrate that the strategy is rather general and could be applied to different pharmacologically relevant approaches. We close this study by outlining one of these approaches, namely the light-controllable association of small molecules with RNA, as an emerging approach in RNA-targeting therapy.

Guidelines for SHAPE Reagent Choice and Detection Strategy for RNA Structure Probing Studies

Chemical probing is an important tool for characterizing the complex folded structures of RNA molecules, many of which play key cellular roles. Electrophilic SHAPE reagents create adducts at the 2'-hydroxyl position on the RNA backbone of flexible ribonucleotides with relatively little dependence on nucleotide identity. Strategies for adduct detection such as mutational profiling (MaP) allow accurate, automated calculation of relative adduct frequencies for each nucleotide in a given RNA or group of RNAs. A number of alternative reagents and adduct detection strategies have been proposed, especially for use in living cells. Here we evaluate five SHAPE reagents: three previously well-validated reagents 1M7 (1-methyl-7-nitroisatoic anhydride), 1M6 (1-methyl-6-nitroisatoic anhydride), and NMIA ( N-methylisatoic anhydride), one more recently proposed NAI (2-methylnicotinic acid imidazolide), and one novel reagent 5NIA (5-nitroisatoic anhydride). We clarify the importance of carefully designed software in reading out SHAPE experiments using massively parallel sequencing approaches. We examine SHAPE modification in living cells in diverse cell lines, compare MaP and reverse transcription-truncation as SHAPE adduct detection strategies, make recommendations for SHAPE reagent choice, and outline areas for future development.

Divalent Metal Ion Activation of a Guanine General Base in the Hammerhead Ribozyme: Insights from Molecular Simulations

The hammerhead ribozyme is a well-studied nucleolytic ribozyme that catalyzes the self-cleavage of the RNA phosphodiester backbone. Despite experimental and theoretical efforts, key questions remain about details of the mechanism with regard to the activation of the nucleophile by the putative general base guanine (G12). Straightforward interpretation of the measured activity-pH data implies the pK_a value of the N1 position in the G12 nucleobase is significantly shifted by the ribozyme environment. Recent crystallographic and biochemical work has identified pH-dependent divalent metal ion binding at the N7/O6 position of G12, leading to the hypothesis that this binding mode could induce a pK_a shift of G12 toward neutrality. We present computational results that support this hypothesis and provide a model that unifies the interpretation of available structural and biochemical data and paints a detailed mechanistic picture of the general base step of the reaction. Experimentally testable predictions are made for mutational and rescue effects on G12, which will give further insights into the catalytic mechanism. These results contribute to our growing knowledge of the potential roles of divalent metal ions in RNA catalysis.

Specific Reaction Parametrization of the AM1/d Hamiltonian for Phosphoryl Transfer Reactions: H, O, and P Atoms

A semiempirical AM1/d Hamiltonian is developed to model phosphoryl transfer reactions catalyzed by enzymes and ribozymes for use in linear-scaling calculations and combined quantum mechanical/molecular mechanical simulations. The model, designated AM1/d-PhoT, is parametrized for H, O, and P atoms to reproduce high-level density-functional results from a recently constructed database of quantum calculations for RNA catalysis ( http://theory.chem.umn.edu/Database/QCRNA ), including geometries and relative energies of minima, transition states and reactive intermediates, dipole moments, proton affinities, and other relevant properties. The model is tested in the gas phase and in solution using a QM/MM potential. The results indicate that the method provides significantly higher accuracy than MNDO/d, AM1, and PM3 methods and, for the transphosphorylation reactions, is in close agreement with the density-functional calculations at the B3LYP/6-311++G(3df,2p) level with a reduction in computational cost of 3-4 orders of magnitude. The model is expected to have considerable impact on the application of semiempirical QM/MM methods to transphosphorylation reactions in solution, enzymes, and ribozymes and to ultimately facilitate the design of improved next-generation multiscale quantum models.

Bridging the gap between theory and experiment to derive a detailed understanding of hammerhead ribozyme catalysis

Herein we summarize our progress toward the understanding of hammerhead ribozyme (HHR) catalysis through a multiscale simulation strategy. Simulation results collectively paint a picture of HHR catalysis: HHR first folds to form an electronegative active site pocket to recruit a threshold occupation of cationic charges, either a Mg(2+) ion or multiple monovalent cations. Catalytically active conformations that have good in-line fitness are supported by specific metal ion coordination patterns that involve either a bridging Mg(2+) ion or multiple Na(+) ions, one of which is also in a bridging coordination pattern. In the case of a single Mg(2+) ion bound in the active site, the Mg(2+) ion undergoes a migration that is coupled with deprotonation of the nucleophile (C17:O2'). As the reaction proceeds, the Mg(2+) ion stabilizes the accumulating charge of the leaving group and significantly increases the general acid ability of G8:O2'. Further computational mutagenesis simulations suggest that the disruptions due to mutations may severely impact HHR catalysis at different stages of the reaction. Catalytic mechanisms supported by the simulation results are consistent with available structural and biochemical experiments, and together they advance our understanding of HHR catalysis.

Alkaline hydrolysis of ethylene phosphate: an ab initio study by supermolecule model and polarizable continuum approach

The alkaline hydrolysis reaction of ethylene phosphate (EP) has been investigated using a supermolecule model, in which several explicit water molecules are included. The structures and single-point energies for all of the stationary points are calculated in the gas phase and in solution at the B3LYP/6-31++G(df,p) and MP2/6-311++G(df,2p) levels. The effect of water bulk solvent is introduced by the polarizable continuum model (PCM). Water attack and hydroxide attack pathways are taken into account for the alkaline hydrolysis of EP. An associative mechanism is observed for both of the two pathways with a kinetically insignificant intermediate. The water attack pathway involves a water molecule attacking and a proton transfer from the attacking water to the hydroxide in the first step, followed by an endocyclic bond cleavage to the leaving group. While in the first step of the hydroxide attack pathway the nucleophile is the hydroxide anion. The calculated barriers in aqueous solution for the water attack and hydroxide attack pathways are all about 22 kcal/mol. The excellent agreement between the calculated and observed values demonstrates that both of the two pathways are possible for the alkaline hydrolysis of EP.

Computational mutagenesis studies of hammerhead ribozyme catalysis

Computational studies of the mutational effects at the C3, G8, and G5 positions of the hammerhead ribozyme (HHR) are reported, based on a series of twenty-four 100-ns molecular dynamics simulations of the native and mutated HHR in the reactant state and in an activated precursor state (G8:2'OH deprotonated). Invoking the assumptions that G12 acts as the general base while the 2'OH of G8 acts as a general acid, the simulations are able to explain the origins of experimentally observed mutational effects, including several that are not easily inferred from the crystal structure. Simulations suggest that the Watson-Crick base-pairing between G8 and C3, the hydrogen bond network between C17 and G5, and the base stacking interactions between G8 and C1.1, collectively, are key to maintaining an active site structure conducive for catalytic activity. Mutation-induced disruption of any of these interactions will adversely affect activity. The simulation results predict that the C3U/G8D double mutant, where D is 2,6-diaminopurine, will have a rescue effect relative to the corresponding single mutations. Two general conclusions about the simulations emerge from this work. First, mutation simulations may require 30 ns or more to suitably relax such that the mutational effects become apparent. Second, in some cases, it is necessary to look beyond the reactant state in order to interpret mutational effects in terms of catalytically active structure. The present simulation results lead to better understanding of the origin of experimental mutational effects and provide insight into the key conserved features necessary to maintain the integrity of the active site architecture.

Role of Mg2+ in hammerhead ribozyme catalysis from molecular simulation

Molecular dynamics simulations have been performed to investigate the role of Mg2+ in the full-length hammerhead ribozyme cleavage reaction. In particular, the aim of this work is to characterize the binding mode and conformational events that give rise to catalytically active conformations and stabilization of the transition state. Toward this end, a series of eight 12 ns molecular dynamics simulations have been performed with different divalent metal binding occupations for the reactant, early and late transition state using recently developed force field parameters for metal ions and reactive intermediates in RNA catalysis. In addition, hybrid QM/MM calculations of the early and late transition state were performed to study the proton-transfer step in general acid catalysis that is facilitated by the catalytic Mg2+ ion. The simulations suggest that Mg2+ is profoundly involved in the hammerhead ribozyme mechanism both at structural and catalytic levels. Binding of Mg2+ in the active site plays a key structural role in the stabilization of stem I and II and to facilitate formation of near attack conformations and interactions between the nucleophile and G12, the implicated general base catalyst. In the transition state, Mg2+ binds in a bridging position where it stabilizes the accumulated charge of the leaving group while interacting with the 2'OH of G8, the implicated general acid catalyst. The QM/MM simulations provide support that, in the late transition state, the 2'OH of G8 can transfer a proton to the leaving group while directly coordinating the bridging Mg2+ ion. The present study provides evidence for the role of Mg2+ in hammerhead ribozyme catalysis. The proposed simulation model reconciles the interpretation of available experimental structural and biochemical data, and provides a starting point for more detailed investigation of the chemical reaction path with combined QM/MM methods.

Structure-function relationships in RNA and RNP enzymes: recent advances

The structural biology of ribozymes and ribonucleoprotein (RNP) enzymes is now sufficiently advanced that a true dialogue between structural and functional studies is possible. In this review, we consider three important systems in which an integration of structural and biochemical data has recently led to major advances in mechanistic understanding. In the hammerhead ribozyme, application-driven biochemical studies led to the discovery of a key structural interaction that had been omitted from previously-studied constructs. A new crystal structure of the resulting, tertiary-stabilized hammerhead has resolved a remarkable number of longstanding paradoxes in the structure-function relationship of this ribozyme. In the Group I intron ribozyme, a flurry of high-resolution structures has largely confirmed, but in some cases refined or challenged, a detailed model of a metalloenzyme active site that had previously been derived by meticulous quantitative metal ion rescue experiments. Finally, for the peptidyl transferase center of the ribosome, recent biochemical and chemical results motivated by the pioneering crystal structures have suggested a picture of a catalytic mechanism dominated by proximity and orientation effects and substrate-assisted catalysis. These results refocus attention on catalysis as a property of the integrated RNP machinery as a whole, as opposed to a narrow concern with the RNA functional groups in immediate contact with the reactive center.

Density functional study of the in-line mechanism of methanolysis of cyclic phosphate and thiophosphate esters in solution: insight into thio effects in RNA transesterification

Density functional calculations of thio effects on the in-line mechanism of methanolysis of ethylene phosphate (a reverse reaction model for RNA phosphate transesterification) are presented. A total of 12 reaction mechanisms are examined using the B3LYP functional with large basis sets, and the effects of solvation were treated using the PCM, CPCM, and SM5 solvation models. Single thio substitutions at all of the distinct phosphoryl oxygen positions (2', 3', 5', pro-R) and a double thio substitution at the nonbridging (pro-R/pro-S) positions were considered. Profiles for each reaction were calculated in the dianionic and monoanionic/monoprotic states, corresponding to reaction models under alkaline and nonalkaline conditions, respectively. These models provide insight into the mechanisms of RNA transesterification thio effects and serve as a set of high-level quantum data that can be used in the design of new semiempirical quantum models for hybrid quantum mechanical/molecular mechanical simulations and linear-scaling electronic structure calculations.

CHARMM force field parameters for simulation of reactive intermediates in native and thio-substituted ribozymes

Force field parameters specifically optimized for residues important in the study of RNA catalysis are derived from density-functional calculations, in a fashion consistent with the CHARMM27 all-atom empirical force field. Parameters are presented for residues that model reactive RNA intermediates and transition state analogs, thio-substituted phosphates and phosphoranes, and bound Mg(2+) and di-metal bridge complexes. Target data was generated via density-functional calculations at the B3LYP/6-311++G(3df,2p)// B3LYP/6-31++G(d,p) level. Partial atomic charges were initially derived from CHelpG electrostatic potential fitting and subsequently adjusted to be consistent with the CHARMM27 charges. Lennard-Jones parameters were determined to reproduce interaction energies with water molecules. Bond, angle, and torsion parameters were derived from the density-functional calculations and renormalized to maintain compatibility with the existing CHARMM27 parameters for standard residues. The extension of the CHARMM27 force field parameters for the nonstandard biological residues presented here will have considerable use in simulations of ribozymes, including the study of freeze-trapped catalytic intermediates, metal ion binding and occupation, and thio effects.

Insight into the role of Mg in hammerhead ribozyme catalysis from X-ray crystallography and molecular dynamics simulation

Results of a series of 12 ns molecular dynamics (MD) simulations of the reactant state (with and without a Mg(2+) ion), early and late transition state mimics are presented based on a recently reported crystal structure of a full-length hammerhead RNA. The simulation results support a catalytically active conformation with a Mg(2+) ion bridging the A9 and scissile phosphates. In the reactant state, the Mg(2+) spends significant time closely associated with the 2'OH of G8, but remains fairly distant from the leaving group O(5') position. In the early TS mimic simulation, where the nucleophilic O(2') and leaving group O(5') are equidistant from the phosphorus, the Mg(2+) ion remains tightly coordinated to the 2'OH of G8, but is positioned closer to the O(5') leaving group, stabilizing the accumulating charge. In the late TS mimic simulation, the coordination around the bridging Mg(2+) ion undergoes a transition whereby the coordination with the 2'OH of G8 is replace by the leaving group O(5') that has developed significant charge. At the same time, the 2'OH of G8 forms a hydrogen bond with the leaving group O(5') and is positioned to act as a general acid catalyst. This work represents the first reported simulations of the full-length hammerhead structure and TS mimics, and provides direct evidence for the possible role of a bridging Mg(2+) ion in catalysis that is consistent with both crystallographic and biochemical data.

QCRNA 1.0: a database of quantum calculations for RNA catalysis

This work outlines a new on-line database of quantum calculations for RNA catalysis (QCRNA) available via the worldwide web at http://theory.chem.umn.edu/QCRNA. The database contains high-level density functional calculations for a large range of molecules, complexes and chemical mechanisms important to phosphoryl transfer reactions and RNA catalysis. Calculations are performed using a strict, consistent protocol such that a wealth of cross-comparisons can be made to elucidate meaningful trends in biological phosphate reactivity. Currently, around 2000 molecules have been collected in varying charge states in the gas phase and in solution. Solvation was treated with both the PCM and COSMO continuum solvation models. The data can be used to study important trends in reactivity of biological phosphates, or used as benchmark data for the design of new semiempirical quantum models for hybrid quantum mechanical/molecular mechanical simulations.

The chemical nature of the 2'-substituent in the pentose-sugar dictates the pseudoaromatic character of the nucleobase (pKa) in DNA/RNA

We here show that the pKa (error limit: 0.01 to 0.03 pKa unit) of a nucleobase in a nucleotide can be modulated by the chemical nature of the 2'-substituent at the sugar moiety. This has been evidenced by the measurement of nucleobase pKa in 47 different model nucleoside 3',5'-bis- and 3'-mono-ethylphosphates. The fact that the electronic character of each of the 2'-substituents (Fig. 1) alters the chemical shift of the H2' sugar proton, and also alters the pKa of the nucleobase in the nucleotides has been evidenced by a correlation plot of pKa of N3 of pyrimidine (T/C/U) or pKa of N7 of 9-guaninyl with the corresponding deltaH2' chemical shifts at the neutral pH, which shows linear correlation with high Pearson's correlation coefficients (R = 0.85-0.97). That this modulation of the pKa of the nucleobase by a 2'-substituent is a through-bond as well as through-space effect has been proven by ab initio determined pKa estimation. Interestingly, experimental pKas of nucleobases from NMR titration and the calculated pKas (by ab initio calculations utilizing closed shell HF 6-31G** basis set) are linearly correlated with R = 0.98. It has also been observed that the difference of ground and protonated/de-protonated HOMO orbital energies (DeltaHOMO, a.u.) for the nucleobases (A/G/C/T/U) are well correlated with their pK(a)s in different 2'-substituted 3',5'-bis-ethylphosphate analogs suggesting that only the orbital energy of HOMO can be successfully used to predict the modulation of the chemical reactivity of the nucleobase by the 2'-substituent. It has also been demonstrated that pKa values of nucleobases in 3',5'-bis-ethylphosphates (Table 1) are well correlated with the change in dipole moment for the respective nucleobases after protonation or de-protonation. This work thus unambiguously shows that alteration of the thermodynamic stability (Tm) of the donor-acceptor complexes [ref. 20], as found with various 2'-modified duplexes in the antisense, siRNA or in triplexes by many workers in the field, is a result of alteration of the pseudoaromatic character of the nucleobases engineered by alteration of the chemical nature of the 2'-substitution.

Pseudorotation Barriers of Biological Oxyphosphoranes: A Challenge for Simulations of Ribozyme Catalysis

Pseudorotation reactions of biologically relevant oxyphosphoranes were studied by using density functional and continuum solvation methods. A series of 16 pseudorotation reactions involving acyclic and cyclic oxyphosphoranes in neutral and monoanionic (singly deprotonated) forms were studied, in addition to pseudorotation of PF5. The effect of solvent was treated by using three different solvation models for comparison. The barriers to pseudorotation ranged from 1.5 to 8.1 kcal mol(-1) and were influenced systematically by charge state, apicophilicity of ligands, intramolecular hydrogen bonding, cyclic structure and solvation. Barriers to pseudorotation for monoanionic phosphoranes occur with the anionic oxo ligand as the pivotal atom, and are generally lower than for neutral phosphoranes. The OCH3 groups were observed to be more apicophilic than OH groups, and hence pseudorotations that involve axial OCH3/equatorial OH exchange had higher reaction and activation free energy values. Solvent generally lowered barriers relative to the gas-phase reactions. These results, together with isotope 18O exchange experiments, support the assertion that dianionic phosphoranes are not sufficiently long-lived to undergo pseudorotation. Comparison of the density functional results with those from several semiempirical quantum models highlight a challenge for new-generation hybrid quantum mechanical/molecular mechanical potentials for non-enzymatic and enzymatic phosphoryl transfer reactions: the reliable modeling of pseudorotation processes.

Hybrid QM/MM study of thio effects in transphosphorylation reactions: the role of solvation

Transphosphorylation thio effects in solution are studied using hybrid QM/MM calculations with a d-orbital semiempirical Hamiltonian. Activated dynamics simulations were performed for a 3' ribose-phosphate model in an explicit 20 A sphere of TIP3P water surrounded by a solvent boundary potential, and free energy analysis was performed using the weighted histogram analysis method. Single thio-substitutions at all of the phosphoryl oxygen positions and a double thio-substitution at the nonbridging positions were considered. The reaction free energy profiles are compared with available experimental data, and the role of solvation on the barrier heights and reaction coordinate is discussed. These results provide an important step in the characterization of thio effects in reactions of biological phosphates that may aid in the interpretation of kinetic data and ultimately help to unravel the catalytic mechanisms of ribozymes.

Existence of efficient divalent metal ion-catalyzed and inefficient divalent metal ion-independent channels in reactions catalyzed by a hammerhead ribozyme

The hammerhead ribozyme is generally accepted as a well characterized metalloenzyme. However, the precise nature of the interactions of the RNA with metal ions remains to be fully defined. Examination of metal ion-catalyzed hammerhead reactions at limited concentrations of metal ions is useful for evaluation of the role of metal ions, as demonstrated in this study. At concentrations of Mn2+ ions from 0.3 to 3 mM, addition of the ribozyme to the reaction mixture under single-turnover conditions enhances the reaction with the product reaching a fixed maximum level. Further addition of the ribozyme inhibits the reaction, demonstrating that a certain number of divalent metal ions is required for proper folding and also for catalysis. At extremely high concentrations, monovalent ions, such as Na+ ions, can also serve as cofactors in hammerhead ribozyme-catalyzed reactions. However, the catalytic efficiency of monovalent ions is extremely low and, thus, high concentrations are required. Furthermore, addition of monovalent ions to divalent metal ion-catalyzed hammerhead reactions inhibits the divalent metal ion-catalyzed reactions, suggesting that the more desirable divalent metal ion-ribozyme complexes are converted to less desirable monovalent metal ion-ribozyme complexes via removal of divalent metal ions, which serve as a structural support in the ribozyme complex. Even though two channels appear to exist, namely an efficient divalent metal ion-catalyzed channel and an inefficient monovalent metal ion-catalyzed channel, it is clear that, under physiological conditions, hammerhead ribozymes are metalloenzymes that act via the significantly more efficient divalent metal ion-dependent channel. Moreover, the observed kinetic data are consistent with Lilley's and DeRose's two-phase folding model that was based on ground state structure analyses.

Recent advances in the elucidation of the mechanisms of action of ribozymes

The cleavage of RNA can be accelerated by a number of factors. These factors include an acidic group (Lewis acid) or a basic group that aids in the deprotonation of the attacking nucleophile, in effect enhancing the nucleophilicity of the nucleophile; an acidic group that can neutralize and stabilize the leaving group; and any environment that can stabilize the pentavalent species that is either a transition state or a short-lived intermediate. The catalytic properties of ribozymes are due to factors that are derived from the complicated and specific structure of the ribozyme-substrate complex. It was postulated initially that nature had adopted a rather narrowly defined mechanism for the cleavage of RNA. However, recent findings have clearly demonstrated the diversity of the mechanisms of ribozyme-catalyzed reactions. Such mechanisms include the metal-independent cleavage that occurs in reactions catalyzed by hairpin ribozymes and the general double-metal-ion mechanism of catalysis in reactions catalyzed by the Tetrahymena group I ribozyme. Furthermore, the architecture of the complex between the substrate and the hepatitis delta virus ribozyme allows perturbation of the pK(a) of ring nitrogens of cytosine and adenine. The resultant perturbed ring nitrogens appear to be directly involved in acid/base catalysis. Moreover, while high concentrations of monovalent metal ions or polyamines can facilitate cleavage by hammerhead ribozymes, divalent metal ions are the most effective acid/base catalysts under physiological conditions.

RNA conformation and folding studied with fluorescence resonance energy transfer

Fluorescence resonance energy transfer (FRET) results from nonradiative coupling of two fluorophores and reports on distances in the range 10-100 A. It is therefore a suitable probe to determine distances in RNA molecules and define their global structure, to follow kinetics of RNA conformational changes during folding in real time, to monitor ion binding, or to analyze conformational equilibria and assess the thermodynamic stability of tertiary structure conformers. Along with the basic principles of steady-state and time-resolved fluorescence resonance energy transfer measurements, approaches to investigate RNA conformational transitions and folding are described and illustrated with selected examples. The versatility of FRET-based techniques has recently been demonstrated by implementations of FRET in high-throughput screening of potential drugs as well as studies of energy transfer that monitor RNA conformational changes on the single-molecule level.

Protein-RNA interactions: a structural analysis

A detailed computational analysis of 32 protein-RNA complexes is presented. A number of physical and chemical properties of the intermolecular interfaces are calculated and compared with those observed in protein-double-stranded DNA and protein-single-stranded DNA complexes. The interface properties of the protein-RNA complexes reveal the diverse nature of the binding sites. van der Waals contacts played a more prevalent role than hydrogen bond contacts, and preferential binding to guanine and uracil was observed. The positively charged residue, arginine, and the single aromatic residues, phenylalanine and tyrosine, all played key roles in the RNA binding sites. A comparison between protein-RNA and protein-DNA complexes showed that whilst base and backbone contacts (both hydrogen bonding and van der Waals) were observed with equal frequency in the protein-RNA complexes, backbone contacts were more dominant in the protein-DNA complexes. Although similar modes of secondary structure interactions have been observed in RNA and DNA binding proteins, the current analysis emphasises the differences that exist between the two types of nucleic acid binding protein at the atomic contact level.

Multiple conformational states of the hammerhead ribozyme, broad time range of relaxation and topology of dynamics

The dynamics of a hammerhead ribozyme was analyzed by measurements of fluorescence-detected temperature jump relaxation. The ribozyme was substituted at different positions by 2-aminopurine (2-AP) as fluorescence indicator; these substitutions do not inhibit catalysis. The general shape of relaxation curves reported from different positions of the ribozyme is very similar: a fast decrease of fluorescence, mainly due to physical quenching, is followed by a slower increase of fluorescence due to conformational relaxation. In most cases at least three relaxation time constants in the time range from a few microseconds to approximately 200 ms are required for fitting. Although the relaxation at different positions of the ribozyme is similar in general, suggesting a global type of ribozyme dynamics, a close examination reveals differences, indicating an individual local response. For example, 2-AP in a tetraloop reports mainly the local loop dynamics known from isolated loops, whereas 2-AP located at the core, e.g. at the cleavage site or its vicinity, also reports relatively large amplitudes of slower components of the ribozyme dynamics. A variant with an A-->G substitution in domain II, resulting in an inactive form, leads to the appearance of a particularly slow relaxation process (tau approximately 200 ms). Addition of Mg(2+) ions induces a reduction of amplitudes and in most cases a general increase of time constants. Differences between the hammerhead variants are clearly demonstrated by subtraction of relaxation curves recorded under corresponding conditions. The changes induced in the relaxation response by Mg(2+) are very similar to those induced by Ca(2+). The relaxation data do not provide any evidence for formation of Mg(2+)-inner sphere complexes in hammerhead ribozymes, because a Mg(2+)-specific relaxation effect was not visible. However, a Mg(2+)-specific effect was found for a dodeca-riboadenylate substituted with 2-AP, showing that the fluorescence of 2-AP is able to indicate inner sphere complexation. Amplitudes and time constants show that the equilibrium constant of inner sphere complexation is 1.2, corresponding to 55% inner sphere state of the Mg(2+) complexes; the rate constant 6.6 x 10(3) s(-1) for inner sphere complexation is relatively low and shows the existence of some barrier(s) on the way to inner sphere complexes.

Composing life

Textbooks often assert that life began with specialized complex molecules, such as RNA, that are capable of making their own copies. This scenario has serious difficulties, but an alternative has remained elusive. Recent research and computer simulations have suggested that the first steps toward life may not have involved biopolymers. Rather, non-covalent protocellular assemblies, generated by catalyzed recruitment of diverse amphiphilic and hydrophobic compounds, could have constituted the first systems capable of information storage, inheritance and selection. A complex chain of evolutionary events, yet to be deciphered, could then have led to the common ancestors of today's free-living cells, and to the appearance of DNA, RNA and protein enzymes.

Construction of new ribozymes requiring short regulator oligonucleotides as a cofactor

A hairpin loop and an oligonucleotide bound to the loop form one-half of the pseudoknot structure. We have designed an allosteric hammerhead ribozyme, which is activated by the introduction of this motif by using a short complementary oligonucleotide as a cofactor. Stem II of the hammerhead ribozyme was substituted with a non-self-complementary loop sequence (loop II) to abolish the cleavage activity. The new ribozyme had almost no cleavage activity of the target RNA. However, it exhibited the cleavage activity in the presence of a cofactor oligoribonucleotide, which is complementary to loop II of the ribozyme. The activity is assumed to be derived from the formation of a pseudo-stem structure between the cofactor oligonucleotide and loop II. The structure including the loop may be similar to the pseudo-half-knot structure. The activation efficiencies of the cofactor oligonucleotides were decreased as the lengths of the oligonucleotides increased, and the ribozyme with a longer loop II was more active than that with a short loop II. Oligoribonucleotides with 3'-dangling purine bases served as efficient cofactors of the ribozyme, and a 2'-O-methyloligonucleotide enhanced the cleavage activity of the ribozyme most efficiently, by as much as about 750-fold as compared with that in the absence of the oligonucleotide. Cofactor oligonucleotides with a cytidine base at the 3'-end also activated a ribozyme with the G10.1.G11.1 mutation, which eliminates the cleavage activity in the wild-type. The binding sites of the oligonucleotide were identified by photo-crosslinking experiments and were found to be the predicted sites in the loop. This is the first report of a design aimed at positively controlling the activity of ribozymes by employing a structural motif. This method can be applied to control the activities of other functional RNAs with hairpin loops.

Allosteric nucleic acid catalysts

Endowing nucleic acid catalysts with allosteric properties provides new prospects for RNA and DNA as controllable therapeutic agents or as sensors of their cognate effector compounds. The ability to engineer RNA catalysts that are allosterically regulated by effector binding has been propelled by the union of modular rational design principles and powerful combinatorial strategies.

Structure and function of the hairpin ribozyme

The hairpin ribozyme belongs to the family of small catalytic RNAs that cleave RNA substrates in a reversible reaction that generates 2',3'-cyclic phosphate and 5'-hydroxyl termini. The hairpin catalytic motif was discovered in the negative strand of the tobacco ringspot virus satellite RNA, where hairpin ribozyme-mediated self-cleavage and ligation reactions participate in processing RNA replication intermediates. The self-cleaving hairpin, hammerhead, hepatitis delta and Neurospora VS RNAs each adopt unique structures and exploit distinct kinetic and catalytic mechanisms despite catalyzing the same chemical reactions. Mechanistic studies of hairpin ribozyme reactions provided early evidence that, like protein enzymes, RNA enzymes are able to exploit a variety of catalytic strategies. In contrast to the hammerhead and Tetrahymena ribozyme reactions, hairpin-mediated cleavage and ligation proceed through a catalytic mechanism that does not require direct coordination of metal cations to phosphate or water oxygens. The hairpin ribozyme is a better ligase than it is a nuclease while the hammerhead reaction favors cleavage over ligation of bound products by nearly 200-fold. Recent structure-function studies have begun to yield insights into the molecular bases of these unique features of the hairpin ribozyme.

Characterization of the Tetrahymena ribozyme folding pathway using the kinetic footprinting reagent peroxynitrous acid

Large RNAs fold into complex structures which determine their biological activities. A full understanding of both RNA structure and dynamics will include the description of the pathways by which these structures are formed. Kinetic footprinting [Sclavi, B., et al. (1997) J. Mol. Biol. 266, 144-159] has been shown to be a powerful method for the study of dynamic processes involving RNA. Here we describe the use of a readily available reagent, peroxynitrous acid, as a kinetic footprinting tool for the study of RNA folding. Hydroxyl radicals generated from this reagent were used to footprint the Tetrahymena ribozyme during its magnesium-dependent folding-in agreement with synchroton X-ray footprinting [Sclavi, B., et al. (1998) Science 279, 1940-1943] and oligonucleotide/hybridization cleavage experiments [Zarrinkar, P. P., and Williamson, J. R. (1994) Science 265, 918-924], this work suggests an ordered, hierarchical folding pathway for the ribozyme. Several slow steps in the folding pathway were observed in the peroxynitrous acid footprinting, but none of these corresponded to the rate-determining step of folding. This suggests that the formation of the global, protected structure is followed by one or more slow local rearrangements to yield the final active structure. These studies illustrate the utility of peroxynitrous acid as a reagent for the elucidation of RNA folding pathways and the study of RNA dynamics.