Effects of Mg2+ and the 2' OH of guanosine on steps required for substrate binding and reactivity with the Tetrahymena ribozyme reveal several local folding transitions

Coordination of two sequential ester-transfer reactions: exogenous guanosine binding promotes the subsequent omegaG binding to a group I intron

Self-splicing of group I introns is accomplished by two sequential ester-transfer reactions mediated by sequential binding of two different guanosine ligands, but it is yet unclear how the binding is coordinated at a single G-binding site. Using a three-piece trans-splicing system derived from the Candida intron, we studied the effect of the prior GTP binding on the later ωG binding by assaying the ribozyme activity in the second reaction. We showed that adding GTP simultaneously with and prior to the esterified ωG in a substrate strongly accelerated the second reaction, suggesting that the early binding of GTP facilitates the subsequent binding of ωG. GTP-mediated facilitation requires C2 amino and C6 carbonyl groups on the Watson–Crick edge of the base but not the phosphate or sugar groups, suggesting that the base triple interactions between GTP and the binding site are important for the subsequent ωG binding. Strikingly, GTP binding loosens a few local structures of the ribozyme including that adjacent to the base triple, providing structural basis for a rapid exchange of ωG for bound GTP.

Predicting RNA folding thermodynamics with a reduced chain representation model

Based on the virtual bond representation for the nucleotide backbone, we develop a reduced conformational model for RNA. We use the experimentally measured atomic coordinates to model the helices and use the self-avoiding walks in a diamond lattice to model the loop conformations. The atomic coordinates of the helices and the lattice representation for the loops are matched at the loop-helix junction, where steric viability is accounted for. Unlike the previous simplified lattice-based models, the present virtual bond model can account for the atomic details of realistic three-dimensional RNA structures. Based on the model, we develop a statistical mechanical theory for RNA folding energy landscapes and folding thermodynamics. Tests against experiments show that the theory can give much more improved predictions for the native structures, the thermal denaturation curves, and the equilibrium folding/unfolding pathways than the previous models. The application of the model to the P5abc region of Tetrahymena group I ribozyme reveals the misfolded intermediates as well as the native-like intermediates in the equilibrium folding process. Moreover, based on the free energy landscape analysis for each and every loop mutation, the model predicts five lethal mutations that can completely alter the free energy landscape and the folding stability of the molecule.

Exploration of the transition state for tertiary structure formation between an RNA helix and a large structured RNA

Docking of the P1 duplex into the pre-folded core of the Tetrahymena group I ribozyme exemplifies the formation of tertiary interactions in the context of a complex, structured RNA. We have applied Phi-analysis to P1 docking, which compares the effects of modifications on the rate constant for docking (k(dock)) with the effects on the docking equilibrium (K(dock)). To accomplish this we used a single molecule fluorescence resonance energy transfer assay that allows direct determination of the rate constants for formation of thermodynamically favorable, as well as unfavorable, states. Modification of the eight groups of the P1 duplex that make tertiary interactions with the core and changes in solution conditions decrease K(dock) up to 500-fold, whereas k(dock) changes by </=2-fold. The absence of effects on k(dock), both from atomic modifications and global perturbations, strongly suggests that the transition state for docking is early and does not closely resemble the docked state. These results, the slow rate of docking of 3s(-1), and the observation that a modification that is expected to increase the degrees of freedom between the P1 duplex and the ribozyme core accelerates docking, suggest a model in which a kinetic trap(s) slows docking substantially. Nonetheless, urea does not increase k(dock), suggesting that there is little change in the exposed surface area between the trapped, undocked state and the transition state. The findings highlight that urea and temperature dependencies can be inadequate to diagnose the presence of kinetic traps in a folding process. The results described here, combined with previous work, provide an in-depth view of an RNA tertiary structure formation event and suggest that large, highly structured RNAs may have local regions that are misordered.

RNA folding energy landscapes

Using a statistical mechanical treatment, we study RNA folding energy landscapes. We first validate the theory by showing that, for the RNA molecules we tested having only secondary structures, this treatment (i) predicts about the same native structures as the Zuker method, and (ii) qualitatively predicts the melting curve peaks and shoulders seen in experiments. We then predict thermodynamic folding intermediates. For one hairpin sequence, unfolding is a simple unzipping process. But for another sequence, unfolding is more complex. It involves multiple stable intermediates and a rezipping into a completely non-native conformation before unfolding. The principle that emerges, for which there is growing experimental support, is that although protein folding tends to involve highly cooperative two-state thermodynamic transitions, without detectable intermediates, the folding of RNA secondary structures may involve rugged landscapes, often with more complex intermediate states.

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

BACKGROUND

The group I intron is an RNA enzyme capable of efficiently catalyzing phosphoryl-transfer reactions. Functional groups that stabilize the chemical transition state of the cleavage reaction have been identified, but they are all located within either the 5'-exon (P1) helix or the guanosine cofactor, which are the substrates of the reaction. Functional groups within the ribozyme active site are also expected to assist in transition-state stabilization, and their role must be explored to understand the chemical basis of group I intron catalysis.

RESULTS

Using nucleotide analog interference mapping and site-specific functional group substitution experiments, we demonstrate that the 2'-OH at A207, a highly conserved nucleotide in the ribozyme active site, specifically stabilizes the chemical transition state by approximately 2 kcal mol-1. The A207 2'-OH only makes its contribution when the U(-1) 2'-OH immediately adjacent to the scissile phosphate is present, suggesting that the 2'-OHs of A207 and U(-1) interact during the chemical step.

CONCLUSIONS

These data support a model in which the 3'-oxyanion leaving group of the transesterification reaction is stabilized by a hydrogen-bonding triad consisting of the 2'-OH groups of U(-1) and A207 and the exocyclic amine of G22. Because all three nucleotides occur within highly conserved non-canonical base pairings, this stabilization mechanism is likely to occur throughout group I introns. Although this mechanism utilizes functional groups distinctive of RNA enzymes, it is analogous to the transition states of some protein enzymes that perform similar phosphoryl-transfer reactions.