Crystal structure of a hairpin ribozyme-inhibitor complex with implications for catalysis

Variational Electrostatic Projection (VEP) Methods for Efficient Modeling of the Macromolecular Electrostatic and Solvation Environment in Activated Dynamics Simulations

New methods for the calculation of electrostatic interactions between the active dynamical region and surrounding external solvated macromolecular environment in hybrid quantum mechanical/molecular mechanical (QM/MM) simulations are presented. The variational electrostatic projection (VEP) method, and related variational reverse-mapping procedure (VEP-RVM) utilize an expansion in Gaussian surface elements for representation of electrostatic interactions. The use of a discretized surface that surrounds the active dynamical region greatly reduces the number of interactions with the particles of the external environment. The methods are tested on two catalytic RNA systems: the hammerhead and the hairpin ribozymes. It is shown that with surface elements numbering from 302 to 1202 points the direct VEP and VEP-RVM methods are able to obtain relative force errors in the range of 0.5-0.05% and 0.09-0.0001%, respectively, using a 4.0 A projection buffer. These results are encouraging and provide an essential step in the development of new variational macromolecular solvent-boundary methods for QM/MM calculations of enzyme reactions.

The pre-tRNA nucleotide base and 2'-hydroxyl at N(-1) contribute to fidelity in tRNA processing by RNase P

Fidelity in tRNA processing by the RNase P RNA from Escherichia coli depends, in part, on interactions with the nucleobase and 2' hydroxyl group of N(-1), the nucleotide immediately upstream of the site of RNA strand cleavage. Here, we report a series of biochemical and structure-function studies designed to address how these interactions contribute to cleavage site selection. We find that simultaneous disruption of cleavage site nucleobase and 2' hydroxyl interactions results in parallel reactions leading to correct cleavage and mis-cleavage one nucleotide upstream (5') of the correct site. Changes in Mg(2+) concentration and pH can influence the fraction of product that is incorrectly processed, with pH effects attributable to differences in the rate-limiting steps for the correct and mis-cleavage reaction pathways. Additionally, we provide evidence that interactions with the 2' hydroxyl group adjacent to the reactive phosphate group also contribute to catalysis at the mis-cleavage site. Finally, disruption of the adjacent 2'-hydroxyl contact has a greater effect on catalysis when pairing between the ribozyme and N(-1) is also disrupted, and the effects of simultaneously disrupting these contacts on binding are also non-additive. One implication of these results is that mis-cleavage will result from any combination of active site modifications that decrease the rate of correct cleavage beyond a certain threshold. Indeed, we find that inhibition of correct cleavage and corresponding mis-cleavage also results from disruption of any combination of active site contacts including metal ion interactions and conserved pairing interactions with the 3' RCCA sequence. Such redundancy in interactions needed for maintaining fidelity may reflect the necessity for multiple substrate recognition in vivo. These studies provide a framework for interpreting effects of substrate modifications on RNase P cleavage fidelity and provide evidence for interactions with the nucleobase and 2' hydroxyl group adjacent to the reactive phosphate group in the transition state.

Crystal structure of a group I intron splicing intermediate

A recently reported crystal structure of an intact bacterial group I self-splicing intron in complex with both its exons provided the first molecular view into the mechanism of RNA splicing. This intron structure, which was trapped in the state prior to the exon ligation reaction, also reveals the architecture of a complex RNA fold. The majority of the intron is contained within three internally stacked, but sequence discontinuous, helical domains. Here the tertiary hydrogen bonding and stacking interactions between the domains, and the single-stranded joiner segments that bridge between them, are fully described. Features of the structure include: (1) A pseudoknot belt that circumscribes the molecule at its longitudinal midpoint; (2) two tetraloop-tetraloop receptor motifs at the peripheral edges of the structure; (3) an extensive minor groove triplex between the paired and joiner segments, P6-J6/6a and P3-J3/4, which provides the major interaction interface between the intron's two primary domains (P4-P6 and P3-P9.0); (4) a six-nucleotide J8/7 single stranded element that adopts a mu-shaped structure and twists through the active site, making critical contacts to all three helical domains; and (5) an extensive base stacking architecture that realizes 90% of all possible stacking interactions. The intron structure was validated by hydroxyl radical footprinting, where strong correlation was observed between experimental and predicted solvent accessibility. Models of the pre-first and pre-second steps of intron splicing are proposed with full-sized tRNA exons. They suggest that the tRNA undergoes substantial angular motion relative to the intron between the two steps of splicing.

U2-U6 RNA folding reveals a group II intron-like domain and a four-helix junction

Intron removal in nuclear precursor mRNA is catalyzed through two transesterification reactions by a multi-megaDalton ribonucleoprotein machine called the spliceosome. A complex between U2 and U6 small nuclear RNAs is a core component of the spliceosome. Here we present an NMR structural analysis of a protein-free U2-U6 complex from Saccharomyces cerevisiae. The observed folding of the U2-U6 complex is a four-helix junction, in which the catalytically important AGC triad base-pairs only within U6 and not with U2. The base-pairing of the AGC triad extends the U6 intramolecular stem-loop (U6 ISL), and the NMR structure of this extended U6 ISL reveals structural similarities with domain 5 of group II self-splicing introns. The observed conformation of the four-helix junction could be relevant to the first, but not the second, step of splicing and may help to position the U6 ISL adjacent to the 5' splice site.

Principles of RNA Compaction: Insights from the Equilibrium Folding Pathway of the P4-P6 RNA Domain in Monovalent Cations

Counterions are required for RNA folding, and divalent metal ions such as Mg(2+) are often critical. To dissect the role of counterions, we have compared global and local folding of wild-type and mutant variants of P4-P6 RNA derived from the Tetrahymena group I ribozyme in monovalent and in divalent metal ions. A remarkably simple picture of the folding thermodynamics emerges. The equilibrium folding pathway in monovalent ions displays two phases. In the first phase, RNA molecules that are initially in an extended conformation enforced by charge-charge repulsion are relaxed by electrostatic screening to a state with increased flexibility but without formation of long-range tertiary contacts. At higher concentrations of monovalent ions, a state that is nearly identical to the native folded state in the presence of Mg(2+) is formed, with tertiary contacts that involve base and backbone interactions but without the subset of interactions that involve specific divalent metal ion-binding sites. The folding model derived from these and previous results provides a robust framework for understanding the equilibrium and kinetic folding of RNA.

Use of both direct and indirect 13C tags for probing nitrogen interactions in hairpin ribozyme models by 15N NMR

We have used the synthesis and 15N NMR study of separate loop A and loop B domains of the hairpin ribozyme to demonstrate that multiple 15N atoms can be incorporated into an RNA strand and be unambiguously distinguished through a combination of direct and indirect tagging by 13C atoms. Absence of 15N chemical shift changes shows that the G8N1 in loop A does not become deprotonated up to pH 8, and that the G21N7 of loop B does not bind to Mg2+.

Observation of internal cleavage and ligation reactions of a ribozyme

We have used single-molecule spectroscopy to untangle conformational dynamics and internal chemistry in the hairpin ribozyme. The active site of the ribozyme is stably formed by docking two internal loops, but upon cleavage undocking is accelerated by two orders of magnitude. The markedly different kinetic properties allow us to differentiate cleaved and ligated forms, and thereby observe multiple cycles of internal cleavage and ligation of a ribozyme in a uniquely direct way. The position of the internal equilibrium is biased toward ligation, but the cleaved ribozyme undergoes several undocking events before ligation, during which products may dissociate. Formation of the stably docked active site, rapid undocking after cleavage, and a strong bias toward ligation should combine to generate a stable circular template for the synthesis of the viral (+) strand and thus ensure a productive replication cycle.

The solution structure of the VS ribozyme active site loop reveals a dynamic "hot-spot"

The VS ribozyme is the largest ribozyme in its class and is also the least structurally characterized thus far. The current working model of the VS ribozyme locates the active site in stem-loop VI. The solution structure of this active site loop was determined using high resolution NMR spectroscopy. The structure reveals that the ground-state conformation of the active site differs significantly from that determined previously from chemical structure probing and mutational analysis of the ribozyme in its active conformation, which contains several looped out bases. In contrast, the base-pairing scheme found for the isolated loop contains three mismatched base-pairs: an A+-C, a G-U wobble, and a sheared G-A base-pair and no looped out bases. Dynamics observed within the active site loop provide insight into the mechanism by which the RNA can rearrange its secondary structure into an "activated" conformation prior to cleavage. These findings lend support to the idea that RNA secondary structure is more fluid than once believed and that a better understanding of structure and dynamic features of ribozymes is required to unravel the intricacies of their catalytic abilities.

Representation, searching and discovery of patterns of bases in complex RNA structures

We describe a graph theoretic method designed to perform efficient searches for substructural patterns in nucleic acid structural coordinate databases using a simplified vectorial representation. Two vectors represent each nucleic acid base and the relative positions of bases with respect to one another are described in terms of distances between the defined start and end points of the vectors on each base. These points comprise the nodes and the distances the edges of a graph, and a pattern search can then be performed using a subgraph isomorphism algorithm. The minimal representation was designed to facilitate searches for complex patterns but was first tested on simple, well-characterised arrangements of bases such as base pairs and GNRA-tetraloop receptor interactions. The method performed very well for these interaction types. A survey of side-by-side base interactions, of which the adenosine platform is the best known example, also locates examples of similar base rearrangements that we consider to be important in structural regulation. A number of examples were found, with GU platforms being particularly prevalent. A GC platform in the RNA of the Thermus thermophilus small ribosomal subunit is in an analogous position to an adenosine platform in other species. An unusual GG platform is also observed close to one of the substrate binding sites in Haloarcula marismortui large ribosomal subunit RNA.

Role of an active site guanine in hairpin ribozyme catalysis probed by exogenous nucleobase rescue

The hairpin ribozyme is a small catalytic RNA with reversible phosphodiester cleavage activity. Biochemical and structural studies exclude a requirement for divalent metal cation cofactors and implicate one active site nucleobase in particular, G8, in the catalytic mechanism. Our previous work demonstrated that the cleavage activity that is lost when G8 is replaced by an abasic residue is restored when certain nucleobases are provided in solution. The specificity and pH dependence of exogenous nucleobase rescue were consistent with several models of the rescue mechanism, including general acid base catalysis, electrostatic stabilization of negative charge in the transition state or a requirement for protonation to facilitate exogenous nucleobase binding. Detailed analyses of exogenous nucleobase rescue for both cleavage and ligation reactions now allow us to refine models of the rescue mechanism. Activity increased with increasing pH for both unmodified ribozyme reactions and unrescued reactions of abasic variants lacking G8. This similarity in pH dependence argues against a role for G8 as a general base catalyst, because G8 deprotonation could not be responsible for the pH-dependent transition in the abasic variant. Exogenous nucleobase rescue of both cleavage and ligation activity increased with decreasing pH, arguing against a role for rescuing nucleobases in general acid catalysis, because a nucleobase that contributes general acid catalysis in the cleavage pathway should provide general base catalysis in ligation. Analysis of the concentration dependence of cytosine rescue at high and low pH demonstrated that protonation promotes catalysis within the nucleobase-bound ribozyme complex but does not stabilize nucleobase binding in the ground state. These results support an electrostatic stabilization mechanism in which exogenous nucleobase binding counters negative charge that develops in the transition state.

The synthesis and properties of oligoribonucleotide-spermine conjugates

Polyamines stabilise nucleic acids against chemical and enzymatic degradation, facilitate the formation of secondary and tertiary structures and enhance cellular uptake. Therefore methods for the syntheses of polyamine-nucleic acid conjugates are of interest. A route for the syntheses of RNA-spermine conjugates has been developed. The polyamine was introduced to the C-5 position of uridine via an ethyl tether and the molecule elaborated into a synthon suitable for oligoribonucleotide assembly. The resultant oligomers were components of the hairpin ribozyme. Characterisation of the spermine-conjugated catalytic RNA revealed that attachment of the polyamine was well tolerated in three of four positions, namely U41, U37 and U34, suggesting that conjugation to C-5 brings about minimal structural perturbation.

Single-molecule enzymology of RNA: essential functional groups impact catalysis from a distance

The hairpin ribozyme is a minimalist paradigm for studying RNA folding and function. In this enzyme, two domains dock by induced fit to form a catalytic core that mediates a specific backbone cleavage reaction. Here, we have fully dissected its reversible reaction pathway, which comprises two structural transitions (docking/undocking) and a chemistry step (cleavage/ligation), by applying a combination of single-molecule fluorescence resonance energy transfer (FRET) assays, ensemble cleavage assays, and kinetic simulations. This has allowed us to quantify the effects that modifications of essential functional groups remote from the site of catalysis have on the individual rate constants. We find that all ribozyme variants show similar fractionations into effectively noninterchanging molecule subpopulations of distinct undocking rate constants. This leads to heterogeneous cleavage activity as commonly observed for RNA enzymes. A modification at the domain junction additionally leads to heterogeneous docking. Surprisingly, most modifications not only affect docking/undocking but also significantly impact the internal chemistry rate constants over a substantial distance from the site of catalysis. We propose that a network of coupled molecular motions connects distant parts of the RNA with its reaction site, which suggests a previously undescribed analogy between RNA and protein enzymes. Our findings also have broad implications for applications such as the action of drugs and ligands distal to the active site or the engineering of allostery into RNA.

Sensing complex regulatory networks by conformationally controlled hairpin ribozymes

The hairpin ribozyme catalyses RNA cleavage by a mechanism utilizing its conformational flexibility during the docking of two independently folded internal loop domains A and B. Based on this mechanism, we designed hairpin ribozyme variants that can be induced or repressed by external effector oligonucleotides influencing the docking process. We incorporated a third domain C to assimilate alternate stable RNA motifs such as a pseudo-half-knot or an internal stem-loop structure. Small sequence changes in domain C allowed targeted switching of ribozyme activity: the same effector oligonucleotide can either serve as an inducer or repressor. The ribozymes were applied to trp leader mRNA, the RNA sequence tightly bound by l-tryptophan-activated trp-RNA-binding attenuation protein (TRAP). When domain C is complementary to this mRNA, ribozyme activity can be altered by annealing trp leader mRNA, then specifically reverted by its TRAP/tryptophan-mediated sequestration. This approach allows to precisely sense the activity status of a protein controlled by its metabolite molecule.

Crystal structure of a self-splicing group I intron with both exons

The discovery of the RNA self-splicing group I intron provided the first demonstration that not all enzymes are proteins. Here we report the X-ray crystal structure (3.1-A resolution) of a complete group I bacterial intron in complex with both the 5'- and the 3'-exons. This complex corresponds to the splicing intermediate before the exon ligation step. It reveals how the intron uses structurally unprecedented RNA motifs to select the 5'- and 3'-splice sites. The 5'-exon's 3'-OH is positioned for inline nucleophilic attack on the conformationally constrained scissile phosphate at the intron-3'-exon junction. Six phosphates from three disparate RNA strands converge to coordinate two metal ions that are asymmetrically positioned on opposing sides of the reactive phosphate. This structure represents the first splicing complex to include a complete intron, both exons and an organized active site occupied with metal ions.

Enzymatic synthesis of phosphoroselenoate DNA using thymidine 5'-(alpha-P-seleno)triphosphate and DNA polymerase for X-ray crystallography via MAD

We report here the first study of enzymatic synthesis of two phosphoroselenoate (PSe) DNAs using the two alpha-Se-TTP diastereomers (Sp and Rp) and DNA polymerase. The experimental results indicate that Klenow equally recognizes the two individual diastereomers at the same level as natural TTP. The incorporations of the PSe groups at the expected sites have been confirmed by the digestion resistance to exonuclease III, and the different patterns of the digestion resistance of DNA I and II indicate the configurational differences of the PSe centers (Sp or Rp). Unlike chemical synthesis, which is limited to short DNAs and where the separation of the PSe DNA diastereomers is necessary, this enzymatic method can be used to prepare longer DNAs without diastereomer separation. This quantitative enzymatic approach is particular valuable for the synthesis of longer DNAs with multiple PSe groups in large scale for their X-ray crystal structure determination by the MAD phasing technique.

Folding of the natural hammerhead ribozyme is enhanced by interaction of auxiliary elements

It has been shown that the activity of the hammerhead ribozyme at microM magnesium ion concentrations is markedly increased by the inclusion of loops in helices I and II. We have studied the effect of such loops on the magnesium ion-induced folding of the ribozyme, using fluorescence resonance energy transfer. We find that with the loops in place, folding into the active conformation occurs in a single step, in the microM range of magnesium ion concentration. Disruption of the loop-loop interaction leads to a reversion to two-step folding, with the second stage requiring mM concentrations of magnesium ion. Sodium ions also promote the folding of the natural form of the ribozyme at high concentrations, but the folding occurs as a two-stage process. The loops clearly act as important auxiliary elements in the function of the ribozyme, permitting folding to occur efficiently under physiological conditions.

A conformational switch controls hepatitis delta virus ribozyme catalysis

Ribozymes enhance chemical reaction rates using many of the same catalytic strategies as protein enzymes. In the hepatitis delta virus (HDV) ribozyme, site-specific self-cleavage of the viral RNA phosphodiester backbone requires both divalent cations and a cytidine nucleotide. General acid-base catalysis, substrate destabilization and global and local conformational changes have all been proposed to contribute to the ribozyme catalytic mechanism. Here we report ten crystal structures of the HDV ribozyme in its pre-cleaved state, showing that cytidine is positioned to activate the 2'-OH nucleophile in the precursor structure. This observation supports its proposed role as a general base in the reaction mechanism. Comparison of crystal structures of the ribozyme in the pre- and post-cleavage states reveals a significant conformational change in the RNA after cleavage and that a catalytically critical divalent metal ion from the active site is ejected. The HDV ribozyme has remarkable chemical similarity to protein ribonucleases and to zymogens for which conformational dynamics are integral to biological activity. This finding implies that RNA structural rearrangements control the reactivity of ribozymes and ribonucleoprotein enzymes.

Chemical syntheses of inhibitory substrates of the RNA-RNA ligation reaction catalyzed by the hairpin ribozyme

The chemical syntheses of RNA oligomers containing modifications on the 5'-carbon of the 5'-terminal nucleoside for crystallographic and mechanistic studies of the hairpin ribozyme are reported. Phosphoramidites 4 and 8 were prepared and used in solid phase syntheses of RNA oligomers containing the sequence 5'-N'UCCUCUCC, where N' indicates either 5'-chloro-5'-deoxyguanosine or 5'-amino-5'-deoxyguanosine, respectively. A ribozyme ligation assay with the 5'-chloro- and 5'-amino-modified RNA oligomers demonstrated their inhibition of the hairpin-catalyzed RNA-RNA ligation reaction.

The Chirality of a Four-Way Helical Junction in RNA

We have used a novel electrophoretic approach to determine the chirality of a four-way helical junction in RNA. From our experiments we conclude that the handedness of the helical junction in the complete ribozyme is opposite to that found in the free junction and is therefore constrained by the interaction between the unpaired loops of the ribozyme.

Selenium derivatization and crystallization of DNA and RNA oligonucleotides for X-ray crystallography using multiple anomalous dispersion

We report here the solid phase synthesis of RNA and DNA oligonucleotides containing the 2'-selenium functionality for X-ray crystallography using multiwavelength anomalous dispersion. We have synthesized the novel 2'-methylseleno cytidine phosphoramidite and improved the accessibility of the 2'-methylseleno uridine phosphoramidite for the synthesis of many selenium-derivatized DNAs and RNAs in large scales. The yields of coupling these Se-nucleoside phosphoramidites into DNA or RNA oligonucleotides were over 99% when 5-(benzylmercapto)-1H-tetrazole was used as the coupling reagent. The UV melting study of A-form dsDNAs indicated that the 2'-selenium derivatization had no effect on the stability of the duplexes with the 3'-endo sugar pucker. Thus, the stems of functional RNA molecules with the same 3'-endo sugar pucker appear to be the ideal sites for the selenium derivatization with 2'-Se-C and 2'-Se-U. Crystallization of the selenium-derivatized oligonucleotides is also reported here. The results demonstrate that this 2'-selenium functionality is suitable for RNA and A-form DNA derivatization in X-ray crystallography.

Conformational flexibility of four-way junctions in RNA

Helical junctions are common architectural features in RNA. They are particularly important in autonomously folding molecules, as exemplified by the hairpin ribozyme. We have used single-molecule fluorescence spectroscopy to study the dynamic properties of the perfect (4H) four-way helical junction derived from the hairpin ribozyme. In the presence of Mg(2+), the junction samples parallel and antiparallel conformations and both stacking conformers, with a bias towards one antiparallel stacking conformer. There is continual interconversion between the forms, such that there are several transitions per second under physiological conditions. Our data suggest that interconversion proceeds via an open intermediate with reduced cation binding in which coaxial stacking between helices is disrupted. The rate of interconversion becomes slower at higher Mg(2+) concentrations, yet the activation barrier decreases under these conditions, indicating that entropic effects are important. Transitions also occur in the presence of Na(+) only; however, the coaxial stacking appears incomplete under these conditions. The polymorphic and dynamic character of the four-way RNA junction provides a source of structural diversity, from which particular conformations required for biological function might be stabilised by additional RNA interactions or protein binding.

Modifications and deletions of helices within the hairpin ribozyme-substrate complex: an active ribozyme lacking helix 1

Within the hairpin ribozyme, structural elements required for formation of the active tertiary structure are localized in two independently folding domains, each consisting of an internal loop flanked by helical elements. Here, we present results of a systematic examination of the relationship between the structure of the helical elements and the ability of the RNA to form the catalytically active tertiary structure. Deletions and mutational analyses indicate that helix 1 (H1) in domain A can be entirely eliminated, while segments of helices 2, 3, and 4 can also be deleted. From these results, we derive a new active minimal ribozyme that contains three helical elements, an internal loop, and a terminal loop. A three-dimensional model of this truncated ribozyme was generated using MC-SYM, and confirms that the catalytic core of the minimized construct can adopt a tertiary structure that is very similar to that of the nontruncated version. A new strategy is described to study the functional importance of various residues and chemical groups and to identify specific interdomain interactions. This approach uses two physically separated and truncated domains derived from the minimal motif.

The kink-turn motif in RNA is dimorphic, and metal ion-dependent

The kink-turn (K-turn) is a new motif in RNA structure that was identified by examination of the crystal structures of the ribosome. We examined the structural and dynamic properties of this element in free solution. The K-turn RNA exists in a dynamic equilibrium between a tightly kinked conformation and a more open structure similar to a simple bulge bend. The highly kinked form is stabilized by the noncooperative binding of metal ions, but a significant population of the less-kinked form is present even in the presence of relatively high concentrations of divalent metal ions. The conformation of the tightly kinked population is in excellent agreement with that of the K-turn structures observed in the ribosome by crystallography. The end-to-end FRET efficiency of this species agrees closely with that of the ribosomal K-turn, and the direction of the bend measured by comparative gel electrophoresis also corresponds very well. These results show that the tightly kinked conformation of the K-turn requires stabilization by other factors, possibly by protein binding, for example. The K-turn is therefore unlikely to be of itself a primary organizing feature in RNA.

The Varkud satellite ribozyme

The VS ribozyme is the largest nucleolytic ribozyme, for which there is no crystal structure to date. The ribozyme consists of five helical sections, organized by two three-way junctions. The global structure has been determined by solution methods, particularly FRET. The substrate stem-loop binds into a cleft formed between two helices, while making a loop-loop contact with another section of the ribozyme. The scissile phosphate makes a close contact with an internal loop (the A730 loop), the probable active site of the ribozyme. This loop contains a particularly critical nucleotide A756. Most changes to this nucleotide lead to three-orders of magnitude slower cleavage, and the Watson-Crick edge is especially important. NAIM experiments indicate that a protonated base is required at this position for the ligation reaction. A756 is thus a strong candidate for nucleobase participation in the catalytic chemistry.

Sequence-Specific Detection of MicroRNAs by Signal-Amplifying Ribozymes

The rational and straightforward design of hairpin ribozymes that can be sequence-specifically induced by external oligonucleotides is described. Due to intrinsic signal amplification, their sensitivity is at least an order of magnitude increased compared to standard molecular beacons. We applied this system to the detection of microRNAs, a recently discovered class of small endogenous RNA molecules that are involved in gene regulation. We show that the cognate microRNA can reliably and sensitively be detected at low concentrations in a mix of other microRNA sequences. These probes may be useful in applications that require direct detection of minute amounts of small DNAs or RNAs.

The origins of RNA catalysis in ribozymes

The discovery of RNA catalysis provided a paradigm shift in biology, insight into the evolution of life on the planet and a challenge to understand its mechanistic origins. RNA has limited catalytic resources that must be used to maximal effect. Consequently, RNA catalysis tends to be multifactorial, with several processes contributing to an overall significant enhancement of reaction rate. These include general acid-base catalysis, electrostatic effects, and substrate orientation and proximity. The main players are the RNA nucleobases and bound metal ions. Although most ribozymes carry out phosphoryl transfer, the same considerations appear to apply to peptidyl transfer in the ribosome.

The hairpin ribozyme

The hairpin ribozyme is a naturally occurring RNA that catalyzes sequence-specific cleavage and ligation of RNA. It has been the subject of extensive biochemical and structural studies, perhaps the most detailed for any catalytic RNA to date. Comparison of the structures of its constituent domains free and fully assembled demonstrates that the RNA undergoes extensive structural rearrangement. This rearrangement results in a distortion of the substrate RNA that primes it for cleavage. This ribozyme is known to achieve catalysis employing exclusively RNA functional groups. Metal ions or other catalytic cofactors are not used. Current experimental evidence points to a combination of at least four mechanistic strategies by this RNA: (1) precise substrate orientation, (2) preferential transition state binding, (3) electrostatic catalysis, and (4) general acid base catalysis.

In vitro selection of second site revertants analysis of the hairpin ribozyme active site

We have used in vitro genetics to evaluate the function and interactions of the conserved base G8 in the hairpin ribozyme catalytic RNA. Second site revertant selection for a G8X mutant, where X is any of the other three natural nucleobases, yielded a family of second site suppressors of the G8U mutant, but not of G8C or G8A, indicating that only G and U can be tolerated at position 8 of the ribozyme. This result is consistent with recent observations that point to the functional importance of G8 N-1 in the chemistry of catalysis by this ribozyme reaction. Suppression of the G8U mutation was observed when changes were made directly across loop A from the mutated base at substrate position +2 or positions +2 and +3 in combination. The same changes made in the context of the natural G8 sequence resulted in a very large drop in activity. Thus, the G8U mutation results in a change in specificity of the ribozyme from 5'-N / GUC-3' to 5'-N / GCU-3'. The results presented imply that G8 interacts directly with U+2 during catalysis. We propose that this interaction favors the correct positioning of the catalytic determinants of G8. The implications for the folding of the ribozyme and the catalytic mechanism are discussed.

A helical twist-induced conformational switch activates cleavage in the hammerhead ribozyme

We have captured the structure of a pre-catalytic conformational intermediate of the hammerhead ribozyme using a phosphodiester tether formed between I and Stem II. This phosphodiester tether appears to mimic interactions in the wild-type hammerhead RNA that enable switching between nuclease and ligase activities, both of which are required in the replicative cycles of the satellite RNA viruses from which the hammerhead ribozyme is derived. The structure of this conformational intermediate reveals how the attacking nucleophile is positioned prior to cleavage, and demonstrates how restricting the ability of Stem I to rotate about its helical axis, via interactions with Stem II, can inhibit cleavage. Analogous covalent crosslinking experiments have demonstrated that imposing such restrictions on interhelical movement can change the hammerhead ribozyme from a nuclease to a ligase. Taken together, these results permit us to suggest that switching between ligase and nuclease activity is determined by the helical orientation of Stem I relative to Stem II.

Zinc finger-dependent HIV-1 nucleocapsid protein-TAR RNA interactions

In the minus-strand transfer step of HIV-1 reverse transcription, the nucleocapsid protein (NC) promotes annealing of the 3' 'R' (repeat) region of the RNA genome to its complementary sequence located in the newly synthesized minus-strand strong-stop DNA. The R region contains the highly stable transactivation response (TAR) RNA hairpin. To gain insights into the molecular details of TAR RNA-NC interactions, we carried out hydroxyl radical footprinting, as well as gel-shift and fluorescence anisotropy binding assays using wild-type and mutant forms of NC. Our results support the conclusion that NC variants with mutations in their zinc finger domains have dramatically altered TAR RNA binding interactions relative to wild-type NC. These data demonstrate that a specific zinc finger architecture is required for optimal TAR RNA binding, and help to explain the requirement for the zinc finger motifs of NC in its role as a nucleic acid chaperone in minus-strand transfer.

A four-way junction accelerates hairpin ribozyme folding via a discrete intermediate

The natural form of the hairpin ribozyme comprises two major structural elements: a four-way RNA junction and two internal loops carried by adjacent arms of the junction. The ribozyme folds into its active conformation by an intimate association between the loops, and the efficiency of this process is greatly enhanced by the presence of the junction. We have used single-molecule spectroscopy to show that the natural form fluctuates among three distinct states: the folded state and two additional, rapidly interconverting states (proximal and distal) that are inherited from the junction. The proximal state juxtaposes the two loop elements, thereby increasing the probability of their interaction and thus accelerating folding by nearly three orders of magnitude and allowing the ribozyme to fold rapidly in physiological conditions. Therefore, the hairpin ribozyme exploits the dynamics of the junction to facilitate the formation of the active site from its other elements. Dynamic interplay between structural elements, as we demonstrate for the hairpin ribozyme, may be a general theme for other functional RNA molecules.

Single-molecule transition-state analysis of RNA folding

How RNA molecules fold into functional structures is a problem of great significance given the expanding list of essential cellular RNA enzymes and the increasing number of applications of RNA in biotechnology and medicine. A critical step toward solving the RNA folding problem is the characterization of the associated transition states. This is a challenging task in part because the rugged energy landscape of RNA often leads to the coexistence of multiple distinct structural transitions. Here, we exploit single-molecule fluorescence spectroscopy to follow in real time the equilibrium transitions between conformational states of a model RNA enzyme, the hairpin ribozyme. We clearly distinguish structural transitions between effectively noninterchanging sets of unfolded and folded states and characterize key factors defining the transition state of an elementary folding reaction where the hairpin ribozyme's two helical domains dock to make several tertiary contacts. Our single-molecule experiments in conjunction with site-specific mutations and metal ion titrations show that the two RNA domains are in a contact or close-to-contact configuration in the transition state even though the native tertiary contacts are at most partially formed. Such a compact transition state without well formed tertiary contacts may be a general property of elementary RNA folding reactions.

Unique tertiary and neighbor interactions determine conservation patterns of Cis Watson-Crick A/G base-pairs

X-ray, phylogenetic and quantum chemical analysis of molecular interactions and conservation patterns of cis Watson-Crick (W.C.) A/G base-pairs in 16S rRNA, 23S rRNA and other molecules was carried out. In these base-pairs, the A and G nucleotides interact with their W.C. edges with glycosidic bonds oriented cis relative to each other. The base-pair is stabilised by two hydrogen bonds, the C1'-C1' distance is enlarged and the G(N2) amino group is left unpaired. Quantum chemical calculations show that, in the absence of other interactions, the unpaired amino group is substantially non-planar due to its partial sp(3) pyramidalization, while the whole base-pair is internally propeller twisted and very flexible. The unique molecular properties of the cis W.C. A/G base-pairs make them distinct from other base-pairs. They occur mostly at the ends of canonical helices, where they serve as interfaces between the helix and other motifs. The cis W.C. A/G base-pairs play crucial roles in natural RNA structures with salient sequence conservation patterns. The key contribution to conservation is provided by the unpaired G(N2) amino group that is involved in a wide range of tertiary and neighbor contacts in the crystal structures. Many of them are oriented out of the plane of the guanine base and utilize the partial sp(3) pyramidalization of the G(N2). There is a lack of A/G to G/A covariation, which, except for the G(N2) position, would be entirely isosteric. On the contrary, there is a rather frequent occurrence of G/A to G/U covariation, as the G/U wobble base-pair has an unpaired amino group in the same position as the cis W.C. G/A base-pair. The cis W.C. A/G base-pairs are not conserved when there is no tertiary or neighbor interaction. Obtaining the proper picture of the interactions and phylogenetic patterns of the cis W.C. A/G base-pairs requires a detailed analysis of the relation between the molecular structures and the energetics of interactions at a level of single H-bonds and contacts.

Primordial genetics: phenotype of the ribocyte

The idea that the ancestors of modern cells were RNA cells (ribocytes) can be investigated by asking whether all essential cellular functions might be performed by RNAs. This requires isolating suitable molecules by selection-amplification when the predicted molecules are presently extinct. In fact, RNAs with many properties required during a period in which RNA was the major macromolecular agent in cells (an RNA world) have been selected in modern experiments. There is, accordingly, reason to inquire how such a ribocyte might appear, based on the properties of the RNAs that composed it. Combining the intrinsic qualities of RNA with the fundamental characteristics of selection from randomized sequence pools, one predicts ribocytes with a cell cycle measured (roughly) in weeks. Such cells likely had a rapidly varying genome, composed of many small genetic and catalytic elements made of tens of ribonucleotides. There are substantial arguments that, at the mid-RNA era, a subset of these nucleotides are reproducibly available and resemble the modern four. Such cells are predicted to evolve rapidly. Instead of modifying preexisting genes, ribocytes frequently draw new functions from an internal pool containing zeptomoles (<1 attomole) of predominantly inactive random sequences.

RNAsoft: A suite of RNA secondary structure prediction and design software tools

DNA and RNA strands are employed in novel ways in the construction of nanostructures, as molecular tags in libraries of polymers and in therapeutics. New software tools for prediction and design of molecular structure will be needed in these applications. The RNAsoft suite of programs provides tools for predicting the secondary structure of a pair of DNA or RNA molecules, testing that combinatorial tag sets of DNA and RNA molecules have no unwanted secondary structure and designing RNA strands that fold to a given input secondary structure. The tools are based on standard thermodynamic models of RNA secondary structure formation. RNAsoft can be found online at http://www.RNAsoft.ca.

A crystallographic study of the binding of 13 metal ions to two related RNA duplexes

Metal ions, and magnesium in particular, are known to be involved in RNA folding by stabilizing secondary and tertiary structures, and, as cofactors, in RNA enzymatic activity. We have conducted a systematic crystallographic analysis of cation binding to the duplex form of the HIV-1 RNA dimerization initiation site for the subtype-A and -B natural sequences. Eleven ions (K+, Pb2+, Mn2+, Ba2+, Ca2+, Cd2+, Sr2+, Zn2+, Co2+, Au3+ and Pt4+) and two hexammines [Co (NH3)6]3+ and [Ru (NH3)6]3+ were found to bind to the DIS duplex structure. Although the two sequences are very similar, strong differences were found in their cation binding properties. Divalent cations bind almost exclusively, as Mg2+, at 'Hoogsteen' sites of guanine residues, with a cation-dependent affinity for each site. Notably, a given cation can have very different affinities for a priori equivalent sites within the same molecule. Surprisingly, none of the two hexammines used were able to efficiently replace hexahydrated magnesium. Instead, [Co (NH3)4]3+ was seen bound by inner-sphere coordination to the RNA. This raises some questions about the practical use of [Co (NH3)6]3+ as a [Mg (H2O)6]2+ mimetic. Also very unexpected was the binding of the small Au3+ cation exactly between the Watson-Crick sites of a G-C base pair after an obligatory deprotonation of N1 of the guanine base. This extensive study of metal ion binding using X-ray crystallography significantly enriches our knowledge on the binding of middleweight or heavy metal ions to RNA, particularly compared with magnesium.

Advances in the Development of Ribozymes and Antisense Oligodeoxynucleotides as Antiviral Agents for Human Papillomaviruses

Urogenital human papillomavirus (HPV) infections are the most common viral sexually transmitted disease in women. On a worldwide basis cervical cancer is the second most prevalent cancer of women. Although HPV infection is not sufficient to induce cancer, the causal relation between high-risk HPV infection and cervical cancer is well established. Over 99% of cervical cancers are positive for high-risk HPV. Therefore, there is a need for newer approaches to treat HPV infection. Two novel approaches for inactivating gene expression involve ribozymes and oligonucleotides. Methods for identification of target genes involved in neoplastic transformation and tumour growth have been established, and these will lead to therapeutic approaches without any damage to normal cellular RNA molecules, which is often associated with conventional therapeutics. Ribozymes and oligonucleotides represent rational antiviral approaches for inhibiting the growth of cervical lesions and carcinomas by interfering with E6/E7 RNA production. The E6 and E7 genes of high-risk HPVs cooperate to immortalize primary epithelial cells and because they are found in cervical cancer are considered the hallmark of cervical cancer. The use and modification of ribozymes and antisense oligodeoxynucleotides can inhibit the growth of HPV-16 and HPV-18 immortalized cells, and tumour cells by eliminating E6/E7 transcript. Hammerhead and hairpin ribozymes have been widely studied because of their potential use for gene therapy and their place as therapeutic tools for cervical cancer is being evaluated. Although antiviral ribozymes and anti-sense molecules have been effective as in vitro or in vivo inhibitors of high-risk HPV-positive cells, none is currently in clinical trial. There are, however, a number of other antisense therapies in Phase I–III clinical trial for several oncogenes.

Biology and polymer physics at the single-molecule level

The ability to look at individual molecules has given us new insights into molecular processes. Examples of our recent work are given to illustrate how behaviour that may otherwise be hidden from view can be clearly seen in single-molecule experiments.

Hepatitis C virus internal ribosome entry site-mediated translation is stimulated by specific interaction of independent regions of human La autoantigen

The human La autoantigen has been shown to interact with the internal ribosome entry site (IRES) of hepatitis C virus (HCV) in vitro. Using a yeast three-hybrid system, we demonstrated that, in addition to full-length La protein, both N- and C-terminal halves were able to interact with HCV IRES in vivo. The exogenous addition of purified full-length and truncated La proteins in rabbit reticulocyte lysate showed dose-dependent stimulation of HCV IRES-mediated translation. However, an additive effect was achieved adding the terminal halves together in the reaction, suggesting that both might play critical roles in achieving full stimulatory activity of the full-length La protein. Using computational analysis, three-dimensional structures of the RNA recognition motifs (RRM) of the La protein were independently modeled. Of the three putative RRMs, RRM2 was predicted to have a good binding pocket for the interaction with the HCV IRES around the GCAC motif near the initiator AUG and RRM3 binds perhaps in a different location. This observation was further investigated by the filter-binding and toe-printing assays. The results presented here strongly suggest that both the N- and C-terminal halves can interact independently with the HCV IRES and are involved in stimulating internal initiation of translation.

Ribozymes: recent advances in the development of RNA tools

The discovery 20 years ago that some RNA molecules, called ribozymes, are able to catalyze chemical reactions was a breakthrough in biology. Over the last two decades numerous natural RNA motifs endowed with catalytic activity have been described. They all fit within a few well-defined types that respond to a specific RNA structure. The prototype catalytic domain of each one has been engineered to generate trans-acting ribozymes that catalyze the site-specific cleavage of other RNA molecules. On the 20th anniversary of ribozyme discovery we briefly summarize the main features of the different natural catalytic RNAs. We also describe progress towards developing strategies to ensure an efficient ribozyme-based technology, dedicating special attention to the ones aimed to achieve a new generation of therapeutic agents.

In vitro selection of adenine-dependent hairpin ribozymes

Adenine-dependent hairpin ribozymes were isolated by in vitro selection from a degenerated hairpin ribozyme population. Two new adenine-dependent ribozymes catalyze their own reversible cleavage in the presence of free adenine. Both aptamers have Mg(2+) requirements for adenine-assisted cleavage similar to the wild-type hairpin ribozyme. Cleavage kinetics studies in the presence of various other small molecules were compared. The data suggest that adenine does not induce RNA self-cleavage in the same manner for both aptamers. In addition, investigations of pH effects on catalytic rates show that both adenine-dependent aptamers are more active in basic conditions, suggesting that they use new acid/base catalytic strategies in which adenine could be involved directly. The discovery of hairpin ribozymes dependent on adenine for their reversible self-cleavage presents considerable biochemical and evolutionary interests because we show that RNA is able to use exogenous reactive molecules to enhance its own catalytic activity. Such a mechanism may have been a means by which the ribozymes of the RNA world enlarged their chemical repertoire.

Structural basis for a lethal mutation in U6 RNA

U6 RNA is essential for nuclear pre-mRNA splicing and has been implicated directly in catalysis of intron removal. The U80G mutation at the essential magnesium binding site of the U6 3' intramolecular stem-loop region (ISL) is lethal in yeast. To further understand the structure and function of the U6 ISL, we have investigated the structural basis for the lethal U80G mutation by NMR and optical spectroscopy. The NMR structure reveals that the U80G mutation causes a structural rearrangement within the ISL resulting in the formation of a new Watson-Crick base pair (C67 x G80), and disrupts a protonated C67 x A79 wobble pair that forms in the wild-type structure. Despite the structural change, the accessibility of the metal binding site is unperturbed, and cadmium titration produces similar phosphorus chemical shift changes for both the U80G mutant and wild-type RNAs. The thermodynamic stability of the U80G mutant is significantly increased (Delta Delta G(fold) = -3.6 +/- 1.9 kcal/mol), consistent with formation of the Watson-Crick pair. Our structural and thermodynamic data, in combination with previous genetic data, suggest that the lethal basis for the U80G mutation is stem-loop hyperstabilization. This hyperstabilization may prevent the U6 ISL melting and rearrangement necessary for association with U4.

Catalytic roles for proton transfer and protonation in ribozymes

Utilization of proton transfer in catalysis, which is well known in the mechanisms of protein enzymes, has been described only relatively recently for RNA enzymes. In this article, we present a current understanding of proton transfer by nucleic acids. Rate enhancement and specificity conferred by general acid-base catalysis are discussed. We also present possibilities for electrostatic catalysis from general acids and bases as well as cationic base pairs. The microenvironments of a large RNA provide the possibility of histidine-like pK(a)s for proton transfer, as well as lysine- and arginine-like pK(a)s for electrostatic catalysis. Discussion on proton transfer focuses on the hepatitis delta virus (HDV) and hairpin ribozymes, with select examples drawn from the protein literature. Discussion on electrostatic catalysis also draws on these two ribozymes, and a postulate for electrostatic catalysis by a cationic base pair in the mechanism of peptidyl transfer in the ribosome is presented. We also provide a perspective on possibilities for phosphoryl transfer mechanisms involving phosphorane intermediates and unusual tautomeric forms of the bases. Lastly, a distinction is made between ground state and "transition state" pK(a)s. We favor a model in which changes in pH lead to changes in the distribution of reactive and nonreactive ionizations of the ribozyme molecules in the ground state, and therefore suggest that "pK(a) changes in the transition state" do not provide an acceptable explanation for observed pH-rate profiles.

Outersphere and innersphere coordinated metal ions in an aminoacyl-tRNA synthetase ribozyme

Metal ions are essential cofactors for various ribozymes. Here we dissect the roles of metal ions in an aminoacyl-tRNA synthetase-like ribozyme (ARS ribozyme), which was evolved in vitro. This ribozyme can charge phenylalanine on tRNA in cis, where it is covalently attached to the 5'-end of tRNA (i.e. a form of precursor tRNA), as well as in trans, where it can act as a catalyst. The presence of magnesium ion is essential for this ribozyme to exhibit full catalytic activity. Metal-dependent kinetics, as well as structural mappings using Tb3+ in competition with Mg2+ or Co(NH3)6(3+), identified two potential metal-binding sites which are embedded near the tRNA-binding site. The high affinity metal-binding site can be filled with either Mg2+ or Co(NH3)6(3+) and thus the activity relies on a metal ion that is fully coordinated with water or ammonium ions. This site also overlaps with the amino acid-binding site, suggesting that the metal ion plays a role in constituting the catalytic core. The weak metal-binding site is occupied only by a metal ion(s) that can form innersphere contacts with ligands in the ribozyme and, hence, Mg2+ can enhance ribozyme activity, but Co(NH3)6(3+) cannot. The experiments described in this work establish the roles of metal ions that have distinct coordination properties in the ARS ribozyme.

Catalytic strategies of the hepatitis delta virus ribozymes

The hepatitis delta virus (HDV) ribozymes are self-cleaving RNA sequences critical to the replication of a small RNA genome. A recently determined crystal structure together with biochemical and biophysical studies provides new insight into the possible catalytic mechanism of these ribozymes. The HDV ribozymes are examples of naturally occurring small ribozymes that catalyze cleavage of the RNA backbone with a rate enhancement of 10(6)- to 10(7)-fold over the uncatalyzed rate. To achieve this level of rate enhancement, the HDV ribozymes have been proposed to employ several catalytic strategies that include the use of metal ions, intrinsic binding energy, and a novel example of general acid-base catalysis with a cytosine side chain acting as a proton donor or acceptor.

A mechanistic framework for co-transcriptional folding of the HDV genomic ribozyme in the presence of downstream sequence

Hepatitis delta virus (HDV) is a circular pathogenic RNA that uses self-cleavage by closely related 84 nt genomic and antigenomic ribozymes to facilitate the replication of its genome. Downstream of each ribozyme is a stretch of nucleotides termed the attenuator that functions to base-pair with and unfold the ribozyme into a rod-like fold. The competing rates of RNA synthesis, ribozyme folding and cleavage, and rod folding are therefore likely to affect the efficiency of co-transcriptional self-cleavage. In these studies, co-transcriptional folding of the genomic ribozyme was assayed in vitro by monitoring co-transcriptional self-cleavage of transcripts having variable lengths of sequence downstream of the ribozyme. Co-transcriptional cleavage data were simulated successfully only with kinetic models in which cleavage-inactive channels were populated during transcription. Partitioning to and escape from these channels was influenced, in part, by whether the available attenuator sequence could form structures with the ribozyme, and by the stability of such structures. Surprisingly, only 23 nt of attenuator were needed for strong inactivation of cleavage. Self-cleavage of certain 3'-virus-containing sequences could be restored, partially, by renaturation; however, self-cleavage of transcripts with a full-length attenuator could not be restored efficiently by renaturation in vitro. This suggests that in the presence of the attenuator, the cleavage-active ribozyme fold is not the thermodynamically most stable species. In accordance with this model, the efficiency of self-cleavage of the ribozyme followed by a full-length attenuator was increased by decreasing the rate of transcription. These results suggest that, in the absence of additional factors, efficient co-transcriptional cleavage of the full-length genomic HDV RNA may require cleavage to occur prior to synthesis of the attenuator.

Transition state stabilization by a catalytic RNA

The hairpin ribozyme catalyzes sequence-specific cleavage of RNA through transesterification of the scissile phosphate. Vanadate has previously been used as a transition state mimic of protein enzymes that catalyze the same reaction. Comparison of the 2.2 angstrom resolution structure of a vanadate-hairpin ribozyme complex with structures of precursor and product complexes reveals a rigid active site that makes more hydrogen bonds to the transition state than to the precursor or product. Because of the paucity of RNA functional groups capable of general acid-base or electrostatic catalysis, transition state stabilization is likely to be an important catalytic strategy for ribozymes.