A hydrogen-bonding triad stabilizes the chemical transition state of a group I ribozyme

Protocells and RNA Self-Replication

The general notion of an "RNA world" is that, in the early development of life on the Earth, genetic continuity was assured by the replication of RNA, and RNA molecules were the chief agents of catalytic function. Assuming that all of the components of RNA were available in some prebiotic locale, these components could have assembled into activated nucleotides that condensed to form RNA polymers, setting the stage for the chemical replication of polynucleotides through RNA-templated RNA polymerization. If a sufficient diversity of RNAs could be copied with reasonable rate and fidelity, then Darwinian evolution would begin with RNAs that facilitated their own reproduction enjoying a selective advantage. The concept of a "protocell" refers to a compartment where replication of the primitive genetic material took place and where primitive catalysts gave rise to products that accumulated locally for the benefit of the replicating cellular entity. Replication of both the protocell and its encapsulated genetic material would have enabled natural selection to operate based on the differential fitness of competing cellular entities, ultimately giving rise to modern cellular life.

Phosphodiester models for cleavage of nucleic acids

Nucleic acids that store and transfer biological information are polymeric diesters of phosphoric acid. Cleavage of the phosphodiester linkages by protein enzymes, nucleases, is one of the underlying biological processes. The remarkable catalytic efficiency of nucleases, together with the ability of ribonucleic acids to serve sometimes as nucleases, has made the cleavage of phosphodiesters a subject of intensive mechanistic studies. In addition to studies of nucleases by pH-rate dependency, X-ray crystallography, amino acid/nucleotide substitution and computational approaches, experimental and theoretical studies with small molecular model compounds still play a role. With small molecules, the importance of various elementary processes, such as proton transfer and metal ion binding, for stabilization of transition states may be elucidated and systematic variation of the basicity of the entering or departing nucleophile enables determination of the position of the transition state on the reaction coordinate. Such data is important on analyzing enzyme mechanisms based on synergistic participation of several catalytic entities. Many nucleases are metalloenzymes and small molecular models offer an excellent tool to construct models for their catalytic centers. The present review tends to be an up to date summary of what has been achieved by mechanistic studies with small molecular phosphodiesters.

The origins of the RNA world

The general notion of an "RNA World" is that, in the early development of life on the Earth, genetic continuity was assured by the replication of RNA and genetically encoded proteins were not involved as catalysts. There is now strong evidence indicating that an RNA World did indeed exist before DNA- and protein-based life. However, arguments regarding whether life on Earth began with RNA are more tenuous. It might be imagined that all of the components of RNA were available in some prebiotic pool, and that these components assembled into replicating, evolving polynucleotides without the prior existence of any evolved macromolecules. A thorough consideration of this "RNA-first" view of the origin of life must reconcile concerns regarding the intractable mixtures that are obtained in experiments designed to simulate the chemistry of the primitive Earth. Perhaps these concerns will eventually be resolved, and recent experimental findings provide some reason for optimism. However, the problem of the origin of the RNA World is far from being solved, and it is fruitful to consider the alternative possibility that RNA was preceded by some other replicating, evolving molecule, just as DNA and proteins were preceded by RNA.

Testing the nearest neighbor model for canonical RNA base pairs: revision of GU parameters

Thermodynamic parameters for GU pairs are important for predicting the secondary structures of RNA and for finding genomic sequences that code for structured RNA. Optical melting curves were measured for 29 RNA duplexes with GU pairs to improve nearest neighbor parameters for predicting stabilities of helixes. The updated model eliminates a prior penalty assumed for terminal GU pairs. Six additional duplexes with the 5'GG/3'UU motif were added to the single representation in the previous database. This revises the ΔG°(37) for the 5'GG/3'UU motif from an unfavorable 0.5 kcal/mol to a favorable -0.2 kcal/mol. Similarly, the ΔG°(37) for the 5'UG/3'GU motif changes from 0.3 to -0.6 kcal/mol. The correlation coefficients between predicted and experimental ΔG°(37), ΔH°, and ΔS° for the expanded database are 0.95, 0.89, and 0.87, respectively. The results should improve predictions of RNA secondary structure.

Analysis of helicase-RNA interactions using nucleotide analog interference mapping

Nucleotide analog interference mapping (NAIM) is a combinatorial approach that probes individual atoms and functional groups in an RNA molecule and identifies those that are important for a specific biochemical function. Here, we show how NAIM can be adapted to reveal functionally important atoms and groups on RNA substrates of helicases. We explain how NAIM can be used to investigate translocation and unwinding mechanisms of helicases and discuss the advantages and limitations of this powerful chemogenetic approach.

Impact of steric constraints on the product distribution of phosphate-branched oligonucleotide models of the large ribozymes

To assess the extent to which steric constraints may influence the product distribution of the reactions of the large ribozymes, phosphate-branched oligonucleotides of varying length and sequence have been synthesized and their alkaline hydrolysis studied over a wide temperature range. At low temperatures, the branching trinucleoside-3',3',5'-monophosphate moiety is hydrolyzed almost exclusively by P-O3' fission. At higher temperatures, P-O5' fission competes, accounting at most for 22% of the overall reaction. The results suggest that steric constraints imposed by the secondary structure of the reaction site may significantly contribute to the observed regioselectivity of the transesterification reactions catalyzed by the large ribozymes.

Tightening of active site interactions en route to the transition state revealed by single-atom substitution in the guanosine-binding site of the Tetrahymena group I ribozyme

Protein enzymes establish intricate networks of interactions to bind and position substrates and catalytic groups within active sites, enabling stabilization of the chemical transition state. Crystal structures of several RNA enzymes also suggest extensive interaction networks, despite RNA's structural limitations, but there is little information on the functional and the energetic properties of these inferred networks. We used double mutant cycles and presteady-state kinetic analyses to probe the putative interaction between the exocyclic amino group of the guanosine nucleophile and the N7 atom of residue G264 of the Tetrahymena group I ribozyme. As expected, the results supported the presence of this interaction, but remarkably, the energetic penalty for introducing a CH group at the 7-position of residue G264 accumulates as the reaction proceeds toward the chemical transition state to a total of 6.2 kcal/mol. Functional tests of neighboring interactions revealed that the presence of the CH group compromises multiple contacts within the interaction network that encompass the reactive elements, apparently forcing the nucleophile to bind and attack from an altered, suboptimal orientation. The energetic consequences of this indirect disruption of neighboring interactions as the reaction proceeds demonstrate that linkage between binding interactions and catalysis hinges critically on the precise structural integrity of a network of interacting groups.

A rearrangement of the guanosine-binding site establishes an extended network of functional interactions in the Tetrahymena group I ribozyme active site

Protein enzymes appear to use extensive packing and hydrogen bonding interactions to precisely position catalytic groups within active sites. Because of their inherent backbone flexibility and limited side chain repertoire, RNA enzymes face additional challenges relative to proteins in precisely positioning substrates and catalytic groups. Here, we use the group I ribozyme to probe the existence, establishment, and functional consequences of an extended network of interactions in an RNA active site. The group I ribozyme catalyzes a site-specific attack of guanosine on an oligonucleotide substrate. We previously determined that the hydrogen bond between the exocyclic amino group of guanosine and the 2'-hydroxyl group at position A261 of the Tetrahymena group I ribozyme contributes to overall catalysis. We now use functional data, aided by double mutant cycles, to probe this hydrogen bond in the individual reaction steps of the catalytic cycle. Our results indicate that this hydrogen bond is not formed upon guanosine binding to the ribozyme but instead forms at a later stage of the catalytic cycle. Formation of this hydrogen bond is correlated with other structural rearrangements in the ribozyme's active site that are promoted by docking of the oligonucleotide substrate into the ribozyme's active site, and disruption of this interaction has deleterious consequences for the chemical transformation within the ternary complex. These results, combined with earlier results, provide insight into the nature of the multiple conformational steps used by the Tetrahymena group I ribozyme to achieve its active structure and reveal an intricate, extended network of interactions that is used to establish catalytic interactions within this RNA's active site.

Phosphorane intermediate vs. leaving group stabilization by intramolecular hydrogen bonding in the cleavage of trinucleoside monophosphates: implications for understanding catalysis by the large ribozymes

Hydrolysis of 2',3'-O-methyleneadenosin-5'-yl 5'-O-methyluridin-2'-yl 5'-O-methyl-2'-trifluoroacetamido-2'-deoxyuridin-3'-yl phosphate (1b) has been followed by HPLC over a wide pH range to study the effects of potential hydrogen bonding interactions of the 2'-trifluoroacetamido function on the rate and product distribution of the reaction. At pH < 2, decomposition of 1b (and its 3',3',5'-isomer 1a) is first-order in hydronium-ion concentration and cleavage of the P-O3' bond of the 2'-trifluoroacetamido-modified nucleoside is slightly favored over cleavage of the P-O5' bond. Between pH 2 and 4, the overall hydrolysis is pH-independent and the P-O3' and P-O5' bonds are cleaved at comparable rates. At pH 5, the reaction becomes first-order in hydroxide-ion concentration, with P-O3' bond cleavage predominating. At 10 mmol L(-1) aqueous sodium hydroxide, no P-O5' bond cleavage is observed. Compared to the 2'-OH counterpart , a modest rate enhancement is observed over the entire pH range studied. The absence of P-O5' fission under alkaline conditions suggests hydrogen bond stabilization of the departing 3'-oxyanion by the neighboring 2'-trifluoroacetamido function.

2'-amino-modified ribonucleotides as probes for local interactions within RNA

The 2'-hydroxyl group plays an integral role in RNA structure and catalysis. This ubiquitous component of the RNA backbone can participate in multiple interactions essential for RNA function, such as hydrogen bonding and metal ion coordination, but the multifunctional nature of the 2'-hydroxyl renders identification of these interactions a significant challenge. By virtue of their versatile physicochemical properties, such as distinct metal coordination preferences, hydrogen bonding properties, and ability to be protonated, 2'-amino-2'-deoxyribonucleotides can serve as tools for probing local interactions involving 2'-hydroxyl groups within RNA. The 2'-amino group can also serve as a chemoselective site for covalent modification, permitting the introduction of probes for investigation of RNA structure and dynamics. In this chapter, we describe the use of 2'-aminonucleotides for investigation of local interactions within RNA, focusing on interactions involving 2'-hydroxyl groups required for RNA structure, function, and catalysis.

Nucleotide analog interference mapping

Nucleotide analog interference mapping (NAIM) is a powerful chemogenetic technique that rapidly identifies chemical groups essential for RNA function. Using a series of phosphorothioate-tagged nucleotide analogs, each carrying different modifications of nucleobase or backbone functionalities, it is possible to simultaneously, yet individually, assess the contribution of particular functional groups to an RNA's activity at every position within the molecule. In contrast to traditional mutagenesis, which modifies RNA on the nucleobase level, the smallest mutable unit in a NAIM analysis is a single atom, providing a detailed description of interactions at critical nucleotides. Because the method introduces modified nucleotides by in vitro transcription, NAIM offers a straightforward and efficient approach to study any RNA that has a selectable function, and it can be applied to RNAs of nearly any length.

Thermodynamic characterization of naturally occurring RNA single mismatches with G-U nearest neighbors

Due to their prevalence and roles in biological systems, single mismatches adjacent to G-U pairs are important RNA structural elements. Since there are only limited experimental values for the stability of single mismatches adjacent to G-U pairs, current algorithms using free energy minimization to predict RNA secondary structure from sequence assign predicted thermodynamic values to these types of single mismatches. Here, thermodynamic data are reported for frequently occurring single mismatches adjacent to at least one G-U pair. This experimental data can be used in place of predicted thermodynamic values in algorithms that predict secondary structure from sequence using free energy minimization. When predicting the thermodynamic contributions of previously unmeasured single mismatches, most algorithms apply the same thermodynamic penalty for an A-U pair adjacent to a single mismatch and a G-U pair adjacent to a single mismatch. A recent study, however, suggests that the penalty for a G-U pair adjacent to a tandem mismatch should be 1.2 +/- 0.1 kcal/mol, and the penalty for an A-U pair adjacent to a tandem mismatch should be 0.5 +/- 0.2 kcal/mol [Christiansen, M. E. and Znosko, B. M. (2008) Biochemistry 47, 4329-4336]. Therefore, the data reported here are combined with the existing thermodynamic dataset of single mismatches, and nearest neighbor parameters are derived for an A-U pair adjacent to a single mismatch (1.1 +/- 0.1 kcal/mol) and a G-U pair adjacent to a single mismatch (1.4 +/- 0.1 kcal/mol).

Probing RNA structure and function by nucleotide analog interference mapping

Nucleotide analog interference mapping (NAIM) can be used to simultaneously, yet individually, identify structurally or catalytically important functional groups within an RNA molecule. Phosphorothioate-tagged nucleotides and nucleotide analogs are randomly incorporated into an RNA of interest by in vitro transcription. The phosphorothioate tag marks the site of substitution and identifies sites at which the modification affects the structure or function of the RNA molecule. This technique has been expanded to include identification of hydrogen bonding pairs (NAIS), ionizable functional groups, metal ion ligands, and the energetics of protein binding (QNAIM). The analogs, techniques, and data analysis used in NAIM are described here.

The 2'-hydroxyl group of the guanosine nucleophile donates a functionally important hydrogen bond in the tetrahymena ribozyme reaction

In the first step of self-splicing, group I introns utilize an exogenous guanosine nucleophile to attack the 5'-splice site. Removal of the 2'-hydroxyl of this guanosine results in a 10 (6)-fold loss in activity, indicating that this functional group plays a critical role in catalysis. Biochemical and structural data have shown that this hydroxyl group provides a ligand for one of the catalytic metal ions at the active site. However, whether this hydroxyl group also engages in hydrogen-bonding interactions remains unclear, as attempts to elaborate its function further usually disrupt the interactions with the catalytic metal ion. To address the possibility that this 2'-hydroxyl contributes to catalysis by donating a hydrogen bond, we have used an atomic mutation cycle to probe the functional importance of the guanosine 2'-hydroxyl hydrogen atom. This analysis indicates that, beyond its role as a ligand for a catalytic metal ion, the guanosine 2'-hydroxyl group donates a hydrogen bond in both the ground state and the transition state, thereby contributing to cofactor recognition and catalysis by the intron. Our findings continue an emerging theme in group I intron catalysis: the oxygen atoms at the reaction center form multidentate interactions that function as a cooperative network. The ability to delineate such networks represents a key step in dissecting the complex relationship between RNA structure and catalysis.

Site-specific fluorescent labeling of large RNAs with pyrene

Pyrene is a useful chromophore for monitoring the tertiary structure and folding of large RNAs. This unit describes the general preparation of a large RNA (>80 nucleotides in length) that has been site-specifically modified with pyrene at the 2'-position of an individual internal nucleotide. A protocol is provided for derivatizing a 2'-amino-RNA oligonucleotide with a suitably activated pyrene reagent. This pyrene-labeled oligonucleotide is then assembled with other RNA(s) either by covalent ligation or by noncovalent hybridization to form a full-length structured RNA, which may then be studied by equilibrium and stopped-flow fluorescence spectroscopy.

Thio effects on the departure of the 3'-linked ribonucleoside from diribonucleoside 3',3'-phosphorodithioate diesters and triribonucleoside 3',3',5'-phosphoromonothioate triesters: implications for ribozyme catalysis

To provide a solid chemical basis for the mechanistic interpretations of the thio effects observed for large ribozymes, the cleavage of triribonucleoside 3',3',5'-phosphoromonothioate triesters and diribonucleoside 3',3'-phosphorodithioate diesters has been studied. To elucidate the role of the neighboring hydroxy group of the departing 3'-linked nucleoside, hydrolysis of 2',3'-O-methyleneadenosin-5'-yl bis[5'-O-methyluridin-3'-yl] phosphoromonothioate (1 a) has been compared to the hydrolysis of 2',3'-O-methyleneadenosin-5'-yl 5'-O-methyluridin-3'-yl 2',5'-di-O-methyluridin-3'-yl phosphoromonothioate (1 b) and the hydrolysis of bis[uridin-3'-yl] phosphorodithioate (2 a) to the hydrolysis of uridin-3'-yl 2',5'-di-O-methyluridin-3'-yl phosphorodithioate (2 b). The reactions have been followed by RP HPLC over a wide pH range. The phosphoromonothioate triesters 1 a,b undergo two competing reactions: the starting material is cleaved to a mixture of 3',3'- and 3',5'-diesters, and isomerized to 2',3',5'- and 2',2',5'-triesters. With phosphorodithioate diesters 2 a,b, hydroxide-ion-catalyzed cleavage of the P--O3' bond is the only reaction detected at pH >6, but under more acidic conditions desulfurization starts to compete with the cleavage. The 3',3'-diesters do not undergo isomerization. The hydroxide-ion-catalyzed cleavage reaction with both 1 a and 2 a is 27 times as fast as that compared with their 2'-O-methylated counterparts 1 b and 2 b. The hydroxide-ion-catalyzed isomerization of the 3',3',5'-triester to 2',3',5'- and 2',2',5'-triesters with 1 a is 11 times as fast as that compared with 1 b. These accelerations have been accounted for by stabilization of the anionic phosphorane intermediate by hydrogen bonding with the 2'-hydroxy function. Thio substitution of the nonbridging oxygens has an almost negligible influence on the cleavage of 3',3'-diesters 2 a,b, but the hydrolysis of phosphoromonothioate triesters 1 a,b exhibits a sizable thio effect, k(PO)/k(PS)=19. The effects of metal ions on the rate of the cleavage of diesters and triesters have been studied and discussed in terms of the suggested hydrogen-bond stabilization of the thiophosphorane intermediates derived from 1 a and 2 a.

Atomic level architecture of group I introns revealed

Twenty-two years after their discovery as ribozymes, the self-splicing group I introns are finally disclosing their architecture at the atomic level. The crystal structures of three group I introns solved at moderately high resolution (3.1-3.8A) reveal a remarkably conserved catalytic core bound to the metal ions required for activity. The structure of the core is stabilized by an intron-specific set of long-range interactions that involves peripheral elements. Group I intron structures thus provide much awaited and extremely valuable snapshots of how these ribozymes coordinate substrate binding and catalysis.

General deoxyribozyme-catalyzed synthesis of native 3'-5' RNA linkages

An elusive goal for nucleic acid enzymology has been deoxyribozymes that ligate RNA rapidly, sequence-generally, with formation of native 3'-5' linkages, and in preparatively useful yield. Using in vitro selection, we have identified Mg2+- and Zn2+-dependent deoxyribozymes that simultaneously fulfill all four of these criteria. The new deoxyribozymes operate under practical incubation conditions and have modest RNA substrate sequence requirements, specifically D downward arrowRA for 9DB1 and A downward arrowR for 7DE5 (D = A, G, or U; R = A or G). These requirements are comparable to those of deoxyribozymes such as 10-23 and 8-17, which are already widely used as biochemical tools for RNA cleavage. We anticipate that the 9DB1 and 7DE5 deoxyribozymes will find immediate practical application for RNA ligation.

Structural Evidence for a Two-Metal-Ion Mechanism of Group I Intron Splicing

We report the 3.4 angstrom crystal structure of a catalytically active group I intron splicing intermediate containing the complete intron, both exons, the scissile phosphate, and all of the functional groups implicated in catalytic metal ion coordination, including the 2'-OH of the terminal guanosine. This structure suggests that, like protein phosphoryltransferases, an RNA phosphoryltransferase can use a two-metal-ion mechanism. Two Mg2+ ions are positioned 3.9 angstroms apart and are directly coordinated by all six of the biochemically predicted ligands. The evolutionary convergence of RNA and protein active sites on the same inorganic architecture highlights the intrinsic chemical capacity of the two-metal-ion catalytic mechanism for phosphoryl transfer.

Functional identification of catalytic metal ion binding sites within RNA

The viability of living systems depends inextricably on enzymes that catalyze phosphoryl transfer reactions. For many enzymes in this class, including several ribozymes, divalent metal ions serve as obligate cofactors. Understanding how metal ions mediate catalysis requires elucidation of metal ion interactions with both the enzyme and the substrate(s). In the Tetrahymena group I intron, previous work using atomic mutagenesis and quantitative analysis of metal ion rescue behavior identified three metal ions (M_A, M_B, and M_C) that make five interactions with the ribozyme substrates in the reaction's transition state. Here, we combine substrate atomic mutagenesis with site-specific phosphorothioate substitutions in the ribozyme backbone to develop a powerful, general strategy for defining the ligands of catalytic metal ions within RNA. In applying this strategy to the Tetrahymena group I intron, we have identified the pro-S_P phosphoryl oxygen at nucleotide C262 as a ribozyme ligand for M_C. Our findings establish a direct connection between the ribozyme core and the functionally defined model of the chemical transition state, thereby extending the known set of transition-state interactions and providing information critical for the application of the recent group I intron crystallographic structures to the understanding of catalysis.

A combination of substrate atomic mutagenesis with site-specific substitutions in the ribozyme backbone allow the ligands of catalytic metal ions to be identified.

Hoechst 33258 selectively inhibits group I intron self-splicing by affecting RNA folding

Fungal pathogens are increasing in prevalence due to an increase in resistant strains and the number of immunocompromised humans. Candida albicans is one of these pathogens, and approximately 40% of strains contain a group I self-splicing intron, which is a potential RNA drug target, in their large subunit rRNA precursor. Here, we report that Hoechst 33258 and derivatives thereof are selective inhibitors of C. albicans group I intron self-splicing with an IC50 of 17 microM in 2 mM Mg2+. Chemical probing of the intron in the presence of Hoechst 33258 reveals that the folding of several nucleotides in the P4/P6 region of the intron is affected. A nucleotide near the J4/5 region is protected from chemical modification in the presence of Hoechst 33258 and several nearby are more reactive; this suggests that this region is the molecule's binding site. These results expand the available information on small-molecule targeting of RNA and suggest that the RNA-targeting scaffold provided by Hoechst may prove valuable in designing compounds that inhibit the functions of RNA.

Phosphodiester cleavage of guanylyl-(3',3')-(2'-amino-2'-deoxyuridine): rate acceleration by the 2'-amino function

Hydrolytic reactions of the structural analogue of guanylyl-(3',3')-uridine, guanylyl-(3',3')-(2'-amino-2'-deoxyuridine), having one of the 2'-hydroxyl groups replaced with an amino function, have been followed by RP HPLC in the pH range 0-13 at 90 degrees C. The results are compared to those obtained earlier with guanylyl-(3',3')-uridine, guanylyl-(3',3')-(2',5'-di-O-methyluridine), and uridylyl-(3',5')-uridine. Under basic conditions (pH > 8), the hydroxide ion-catalyzed cleavage of the P-O3' bond (first-order in [OH(-)]) yields a mixture of 2'-amino-2'-deoxyuridine and guanosine 2',3'-cyclic phosphate which is hydrolyzed to guanosine 2'- and 3'-phosphates. Under these conditions, guanylyl-(3',3')-(2'-amino-2'-deoxyuridine) is 10 times less reactive than guanylyl-(3',3')-uridine. Under acidic and neutral conditions (pH 3-8), where the pH-rate profile for the cleavage consists of two pH-independent regions (from pH 3 to pH 4 and from 6 to 8), guanylyl-(3',3')-(2'-amino-2'-deoxyuridine) is considerably reactive. For example, in the latter pH range, guanylyl-(3',3')-(2'-amino-2'-deoxyuridine) is more than 2 orders of magnitude more labile than guanylyl-(3',3')-(2',5'-di-O-methyluridine), while in the former pH range the reactivity difference is 1 order of magnitude. Under very acidic conditions (pH < 3), the isomerization giving guanylyl-(2',3')-(2'-amino-2'-deoxyuridine) and depurination yielding guanine (both first-order in [H(+)]) compete with the cleavage. The Zn(2+)-promoted cleavage ([Zn(2+)] = 5 mmol L(-)(1)) is 15 times faster than the uncatalyzed reaction at pH 5.6. The mechanisms of the reactions of guanylyl-(3',3')-(2'-amino-2'-deoxyuridine) are discussed, particularly focusing on the possible stabilization of phosphorane intermediate and/or transition state via an intramolecular hydrogen bonding by the 2'-amino group.

Highly regio- and enantioselective thermal [2 + 2] cycloaddition of coumarin in a crystalline inclusion complex under high vacuum

Coumarin and thiocoumarin, which are usually thermally unreactive for dimerization, were dimerized in their inclusion crystalline complexes under thermal and high vacuum conditions.

Acidity of secondary hydroxyls in ATP and adenosine analogues and the question of a 2',3'-hydrogen bond in ribonucleosides

As sugar-modified nucleosides, nucleotides, and oligonucleotides are much used in pharmacology and enzymatic studies, knowledge on the relative pKa values of the secondary hydroxyls can be of great help in design and interpretations. To obtain this, we have determined the acidity constants for ATP, dATP, 2'-F-dATP, 2'-NH2-ATP, and 2'-O-methyl ATP in aqueous solution. The influence of the relative acidities seem to be mainly from inductive effects since a good correlation between pKa values and group electronegativity is found. There is no clear energetic contribution from a much suggested H-bond between the 2'-OH and 3'-oxyanion in adenosine. To clarify if this kind of H-bond would make a more-prominent energetic contribution in less-polar solvents, we also determined the acidity of secondary hydroxyls in adenosine, 2'-O-methyladenosine, and 3'-O-methyladenosine in water, methanol, and DMSO. The relative differences in pKa values were, however, quite similar in all solvents, suggesting that no major energetic contribution is made by an intramolecular H-bond in adenosine.

Crystal structure of a phage Twort group I ribozyme–product complex

Group I introns are catalytic RNAs capable of orchestrating two sequential phosphotransesterification reactions that result in self-splicing. To understand how the group I intron active site facilitates catalysis, we have solved the structure of an active ribozyme derived from the orf142-I2 intron from phage Twort bound to a four-nucleotide product RNA at a resolution of 3.6 A. In addition to the three conserved domains characteristic of all group I introns, the Twort ribozyme has peripheral insertions characteristic of phage introns. These elements form a ring that completely envelops the active site, where a snug pocket for guanosine is formed by a series of stacked base triples. The structure of the active site reveals three potential binding sites for catalytic metals, and invokes a role for the 2' hydroxyl of the guanosine substrate in organization of the active site for catalysis.

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.

Substrate-assisted catalysis of peptide bond formation by the ribosome

The ribosome accelerates the rate of peptide bond formation by at least 10(7)-fold, but the catalytic mechanism remains controversial. Here we report evidence that a functional group on one of the tRNA substrates plays an essential catalytic role in the reaction. Substitution of the P-site tRNA A76 2' OH with 2' H or 2' F results in at least a 10(6)-fold reduction in the rate of peptide bond formation, but does not affect binding of the modified substrates. Such substrate-assisted catalysis is relatively uncommon among modern protein enzymes, but it is a property predicted to be essential for the evolution of enzymatic function. These results suggest that substrate assistance has been retained as a catalytic strategy during the evolution of the prebiotic peptidyl transferase center into the modern ribosome.

An Atomic Mutation Cycle for Exploring RNA's 2‘-Hydroxyl Group

The 2'-hydroxyl group fulfills numerous structural and functional roles in RNA, including those of hydrogen bond donor and acceptor. While loss of function upon 2'-deoxynucleotide substitution establishes the importance of specific 2'-hydroxyl groups within RNA, this approach provides no information about how these hydroxyl groups impart their functional contribution. We use an atomic mutation cycle to evaluate the functional importance of the 2'-hydroxyl group's hydrogen atom. Using the Tetrahymena ribozyme reaction, we challenge the cycle to expose the catalytic contribution of the cleavage site 2'-hydroxyl group and its associated hydrogen bond network. The results establish the viability of this cycle as an approach to reveal 2'-hydroxyl groups that donate functionally significant hydrogen bonds.

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.

Catalysis by RNase P RNA: unique features and unprecedented active site plasticity

Metal ions are essential cofactors for precursor tRNA (ptRNA) processing by bacterial RNase P. The ribose 2'-OH at nucleotide (nt) -1 of ptRNAs is known to contribute to positioning of catalytic Me2+. To investigate the catalytic process, we used ptRNAs with single 2'-deoxy (2'-H), 2'-amino (2'-N), or 2'-fluoro (2'-F) modifications at the cleavage site (nt -1). 2' modifications had small (2.4-7.7-fold) effects on ptRNA binding to E. coli RNase P RNA in the ground state, decreasing substrate affinity in the order 2'-OH > 2'-F > 2'-N > 2'-H. Effects on the rate of the chemical step (about 10-fold for 2'-F, almost 150-fold for 2'-H and 2'-N) were much stronger, and, except for the 2'-N modification, resembled strikingly those observed in the Tetrahymena ribozyme-catalyzed reaction at corresponding position. Mn2+ rescued cleavage of the 2'-N but also the 2'-H-modified ptRNA, arguing against a direct metal ion coordination at this location. Miscleavage between nt -1 and -2 was observed for the 2'-N-ptRNA at low pH (further influenced by the base identities at nt -1 and +73), suggesting repulsion of a catalytic metal ion due to protonation of the amino group. Effects caused by the 2'-N modification at nt -1 of the substrate allowed us to substantiate a mechanistic difference in phosphodiester hydrolysis catalyzed by Escherichia coli RNase P RNA and the Tetrahymena ribozyme: a metal ion binds next to the 2' substituent at nt -1 in the reaction catalyzed by RNase P RNA, but not at the corresponding location in the Tetrahymena ribozyme reaction.

Artificial modules for enhancing rate constants of a Group I intron ribozyme without a P4-P6 core element

In this paper we report newly selected artificial modules that enhance the kcat values comparable with or higher than those of the wild-type ribozyme with broad substrate specificity. The elements required for the catalysis of Group I intron ribozymes are concentrated in the P3-P7 domain of their core region, which consists of two conserved helical domains, P4-P6 and P3-P7. Previously, we reported the in vitro selection of artificial modules residing at the peripheral region of a mutant Group I ribozyme lacking P4-P6. We found that derivatives of the ribozyme containing the modules performed the reversal of the first step of the self-splicing reaction efficiently by using their affinity to the substrate RNA, although their kcat values and substrate specificity were uninfluenced and limited, respectively. The results show that it is possible to add a variety of new domains at the peripheral region that play a role comparable with that of the conserved P4-P6 domain.

Optimization and generality of a small deoxyribozyme that ligates RNA

In vitro evolution was previously used to identify a small deoxyribozyme, 7Q10, that ligates RNA with formation of a 2'-5' phosphodiester linkage from a 2',3'-cyclic phosphate and a 5'-hydroxyl group. Ligation occurs in a convenient "binding arms" format analogous to that of the well-known 10-23 and 8-17 RNA-cleaving deoxyribozymes. Here, we report the optimization and generality of 7Q10 as a 2'-5' RNA ligase. By comprehensive mutagenesis of its 16-nucleotide enzyme region, the parent 7Q10 sequence is shown to be optimal for RNA ligation yield, although several mutations are capable of increasing the ligation rate approximately fivefold at the expense of yield. The 7Q10 deoxyribozyme ligates any RNA substrates that form the sequence motif UA GR (arrowhead=ligation site and R=purine), providing at least 30% yield of ligated RNA in approximately 1-2 hours at 37 degrees C and pH 9.0. Comparable yields are obtained in approximately 12-24 hours at pH 7.5, which may be more suitable for larger RNAs that are more sensitive to non-specific degradation. For RNA substrates that form the related ligation junction UA GY (Y=pyrimidine), somewhat lower yields are obtained, but significant ligation activity is still observed. These data establish that 7Q10 is a generally applicable RNA ligase. A plot of log(k(obs)) versus pH from pH 6.9 to 9.0 has a slope of just under 1, suggesting that a single deprotonation occurs during the rate-determining reaction step. The compact 7Q10 deoxyribozyme has both practical utility and the potential for increasing our structural and mechanistic understanding of how nucleic acids can mediate chemical reactions.

Deoxyribozymes with 2'-5' RNA ligase activity

In vitro selection was used to identify deoxyribozymes that ligate two RNA substrates. In the ligation reaction, a 2'-5' RNA phosphodiester linkage is created from a 2',3'-cyclic phosphate and a 5'-hydroxyl group. The new Mg(2+)-dependent deoxyribozymes provide 50-60% yield of ligated RNA in overnight incubations at pH 7.5 and 37 degrees C, and they afford 40-50% yield in 1 h at pH 9.0 and 37 degrees C. Various RNA substrate sequences may be joined by simple Watson-Crick covaration of the DNA binding arms that interact with the two RNA substrates. The current deoxyribozymes have some RNA substrate sequence requirements at the nucleotides immediately surrounding the ligation junction (either UAUA GGAA or UAUN GGAA, where the arrow denotes the ligation site and N equals any nucleotide). One of the new deoxyribozymes was used to prepare by ligation the Tetrahymena group I intron RNA P4-P6 domain, a representative structured RNA. Nondenaturing gel electrophoresis revealed that a 2'-5' linkage between nucleotides A233 and G234 of P4-P6 does not disrupt its Mg(2+)-dependent folding (DeltaDeltaG degrees ' < 0.2 kcal/mol). This demonstrates that a 2'-5' linkage does not necessarily interfere with structure in a folded RNA. Therefore, these non-native linkages may be acceptable in modified RNAs when structure/function relationships are investigated. Deoxyribozymes that ligate RNA should be particularly useful for preparing site-specifically modified RNAs for studies of RNA structure, folding, and catalysis.

Catalysis of amide synthesis by RNA phosphodiester and hydroxyl groups

The functional groups found among the RNA bases and in the phosphoribose backbone represent a limited repertoire from which to construct a ribozyme active site. This work investigates the possibility that simple RNA phosphodiester and hydroxyl functional groups could catalyze amide bond synthesis. Reaction of amine groups with activated esters would be catalyzed by a group that stabilizes the partial positive charge on the amine nucleophile in the transition state. 2'-Amine substitutions adjacent to 3'-phosphodiester or 3'-hydroxyl groups react efficiently with activated esters to form 2'-amide and peptide products. In contrast, analogs in which the 3'-phosphodiester is replaced by an uncharged phosphotriester or is constrained in a distal conformation react at least 100-fold more slowly. Similarly, a nucleoside in which the 3'-hydroxyl group is constrained trans to the 2'-amine is also unreactive. Catalysis of synthetic reactions by RNA phosphodiester and ribose hydroxyl groups is likely to be even greater in the context of a preorganized and solvent-excluding catalytic center. One such group is the 2'-hydroxyl of the ribosome-bound P-site adenosine substrate, which is close to the amine nucleophile in the peptidyl synthesis reaction. Given ubiquitous 2'-OH groups in RNA, there exists a decisive advantage for RNA over DNA in catalyzing reactions of biological significance.

Monitoring intermediate folding states of the td group I intron in vivo

Group I introns consist of two major structural domains, the P4-P6 and P3-P9 domains, which assemble through interactions with peripheral extensions to fold into an active ribozyme. To assess group I intron folding in vivo, we probed the structure of td wild-type and mutant introns using dimethyl sulfate. The results suggest that the majority of the intron population is in the native state in accordance with the current structural model, which was refined to include two novel tertiary contacts. The importance of the loop E motif in the P7.1-P7.2 extension in assisting ribozyme folding was deduced from modeling and mutational analyses. Destabilization of stem P6 results in a deficiency in tertiary structure formation in both major domains, while weakening of stem P7 only interferes with folding of the P3-P9 domain. The different impact of mutations on the tertiary structure suggests that they interfere with folding at different stages. These results provide a first insight into the structure of folding intermediates and suggest a putative order of events in a hierarchical folding pathway in vivo.

Modular engineering of a Group I intron ribozyme

All Group I intron ribozymes contain a conserved core region consisting of two helical domains, P4-P6 and P3-P7. Recent studies have demonstrated that the elements required for catalysis are concentrated in the P3-P7 domain. We carried out in vitro selection experiments by using three newly constructed libraries on a variant of the T4 td Group I ribozyme containing only a P3-P7 domain in its core. Selected variants with new peripheral elements at L7.1, L8 or L9 after nine cycles efficiently catalyzed the reversal reaction of the first step of self-splicing. The variants from this selection contained a short sequence complementary to the substrate RNA without exception. The most active variant, which was 3-fold more active than the parental wild-type ribozyme, was developed from the second selection by employing a clone from the first selection. The results show that the P3-P7 domain can stand as an independent catalytic module to which a variety of new domains for enhancing the activity of the ribozyme can be added.

Leaving group stabilization by metal ion coordination and hydrogen bond donation is an evolutionarily conserved feature of group I introns

To understand the behavior of group I introns on a biologically fundamental level, we must distinguish those traits that arise as the products of natural selection (selected traits) from those that arise as the products of neutral drift (non-selected traits). In practice, this distinction relies on comparing the similarities and differences among widely divergent introns to identify conserved traits. Here we address whether the strategies used by the eukaryotic group I intron from the Tetrahymena ciliate to stabilize the leaving group during splicing are maintained in the group I intron from the widely divergent Azoarcus bacterium. A substrate analogue containing a 3'-phosphorothiolate linkage, in which a sulfur atom replaces the bridging 3'-oxygen atom of the scissile phosphate, reacts 20-fold slower in the Azoarcus reaction than the corresponding unmodified substrate in the presence of Mg(II) as the only divalent cation. However, Mn(II) relieves this negative effect such that the 3'-S-P bond cleaves 21-fold faster than does the 3'O-P bond. Other thiophilic divalent metal ions such as Co(II), Cd(II), and Zn(II) similarly support cleavage of the S-P bond. These results indicate that a metal ion directly coordinates to the leaving group in the transition state of the Azoarcus ribozyme reaction. Additionally, the 3'-sulfur substitution eliminates the approximately 10(3)-fold contribution of the adjacent 2'-OH to transition state stabilization. Considering that sulfur accepts hydrogen bonds weakly compared to oxygen, this result suggests that the 2'-OH contributes to catalysis by donating a hydrogen bond to the 3'-oxygen leaving group in the transition state, presumably acting in conjunction with the metal ion to stabilize the developing negative charge. These same catalytic strategies of metal ion coordination and hydrogen bond donation operate in the Tetrahymena ribozyme reaction, suggesting that these features of catalysis have been conserved during evolution and thus extend to all group I introns. The two ribozymes also exhibit quantitative differences in their response to 3'-sulfur substitution. The Azoarcus ribozyme binds and cleaves the phosphorothiolate substrate more efficiently relative to the natural substrate than the Tetrahymena ribozyme under the same conditions, suggesting that the Azoarcus ribozyme better accommodates the phosphorothiolate at the active site both in the ground state and in the transition state. These differences may reflect either a less tightly knit Azoarcus structure and/or spatial deviations between backbone atoms in the two ribozymes that arise during divergent evolution, analogous to the well-documented relationship between protein sequence and structure.

Self-splicing of the Tetrahymena group I ribozyme without conserved base-triples

BACKGROUND

Group I introns share a conserved core region consisting of two domains, P8-P3-P7 and P4-P6, joined by four base-triples. We showed previously that the T4 td intron can perform phosphoester transfer reactions at two splice sites in the absence of both P4-P6 and the conserved base-triples, whereas it is barely able to perform the intact splicing reaction due to the difficulty of conducting the sequential reactions.

RESULTS

Based on previous findings, we constructed a bimolecular ribozyme lacking a large portion of P4-P6 and the base-triples from the Tetrahymena intron, on the assumption that the long-range interactions of the peripheral regions in the two RNAs can compensate for the deteriorated core. The bimolecular ribozyme performed the intact splicing reaction.

CONCLUSION

The present analysis indicates that the base-triples are nonessential, but that L4 and the distal part of P4 in P4-P6 are important for conducting the splicing reaction. The reconstituted self-splicing ribozyme provides an amenable system for analysing the role(s) of elements in the core region in the self-splicing reaction mechanism.

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.

Efficient construction of long DNA duplexes with internal non-Watson-Crick motifs and modifications

We have developed a semi-synthetic approach for preparing long stretches of DNA (>100 bp) containing internal chemical modifications and/or non-Watson-Crick structural motifs which relies on splint-free, cell-free DNA ligations and recycling of side-products by non-PCR thermal cycling. A double-stranded DNA PCR fragment containing a polylinker in its middle is digested with two restriction enzymes and a small insert ( approximately 20 bp) containing the modification or non-Watson-Crick motif of interest is introduced into the middle. Incorrect products are recycled to starting materials by digestion with appropriate restriction enzymes, while the correct product is resistant to digestion since it does not contain these restriction sites. This semi-synthetic approach offers several advantages over DNA splint-mediated ligations, including fewer steps, substantially higher yields ( approximately 60% overall yield) and ease of use. This method has numerous potential applications, including the introduction of modifications such as fluorophores and cross-linking agents into DNA, controlling the shape of DNA on a large scale and the study of non-sequence-specific nucleic acid-protein interactions.

A phosphoramidate substrate analog is a competitive inhibitor of the Tetrahymena group I ribozyme

BACKGROUND

Phosphoramidate oligonucleotide analogs containing N3'-P5' linkages share many structural properties with natural nucleic acids and can be recognized by some RNA-binding proteins. Therefore, if the N-P bond is resistant to nucleolytic cleavage, these analogs may be effective substrate analog inhibitors of certain enzymes that hydrolyze RNA. We have explored the ability of the Tetrahymena group I intron ribozyme to bind and cleave DNA and RNA phosphoramidate analogs.

RESULTS

The Tetrahymena group I ribozyme efficiently binds to phosphoramidate oligonucleotides but is unable to cleave the N3'-P5' bond. Although it adopts an A-form helical structure, the deoxyribo-phosphoramidate analog, like DNA, does not dock efficiently into the ribozyme catalytic core. In contrast, the ribo-phosphoramidate analog docks similarly to the native RNA substrate, and behaves as a competitive inhibitor of the group I intron 5' splicing reaction.

CONCLUSIONS

Ribo-N3'-P5' phosphoramidate oligonucleotides are useful tools for structural and functional studies of ribozymes as well as protein-RNA interactions.

A chemical phylogeny of group I introns based upon interference mapping of a bacterial ribozyme

Despite its small size, the 205 nt group I intron from Azoarcus tRNA(Ile) is an exceptionally stable self-splicing RNA. This IC3 class intron retains the conserved secondary structural elements common to group I ribozymes, but lacks several peripheral helices. These features make it an ideal system to establish the conserved chemical basis of group I intron activity. We collected nucleotide analog interference mapping (NAIM) data of the Azoarcus intron using 14 analogs that modified the phosphate backbone, the ribose sugar, or the purine base functional groups. In conjunction with a complete interference set collected on the Tetrahymena group I intron (IC1 class), these data define a "chemical phylogeny" of functional groups that are important for the activity of both introns and that may be common chemical features of group I intron catalysts. The data identify the functional moieties most likely to play a conserved role as ligands for catalytic metal ions, the substrate helix, and the guanosine cofactor. These include backbone functional groups whose nucleotide identity is not conserved, and hence are difficult to identify by standard phylogenetic sequence comparisons. The data suggest that both introns utilize an equivalent set of long range tertiary interactions for 5'-splice site selection between the P1 substrate helix and its receptor in the J4/5 asymmetric bulge, as well as an equivalent set of 2'-OH groups for P1 helix docking into most of the single stranded segment J8/7. However, the Azoarcus intron appears to make an alternative set of interactions at the base of the P1 helix and at the 5'-end of the J8/7. Extensive differences were observed within the intron peripheral domains, particularly in P2 and P8 where the Azoarcus data strongly support the proposed formation of a tetraloop-tetraloop receptor interaction. This chemical phylogeny for group I intron catalysis helps to refine structural models of the RNA active site and identifies functional groups that should be carefully investigated for their role in transition state stabilization.