Cooperative tertiary interaction network guides RNA folding

Noncoding RNAs form unique 3D structures, which perform many regulatory functions. To understand how RNAs fold uniquely despite a small number of tertiary interaction motifs, we mutated the major tertiary interactions in a group I ribozyme by single-base substitutions. The resulting perturbations to the folding energy landscape were measured using SAXS, ribozyme activity, hydroxyl radical footprinting, and native PAGE. Double- and triple-mutant cycles show that most tertiary interactions have?a small effect on the stability of the native state. Instead, the formation of core and peripheral structural motifs is cooperatively linked in near-native folding intermediates, and this cooperativity depends on the native helix orientation. The emergence of a cooperative interaction network at an early stage of folding suppresses nonnative structures and guides the search for the native state. We suggest that cooperativity in noncoding RNAs arose from natural selection of architectures conducive to forming?a unique, stable fold.     

Entropy-driven folding of an RNA helical junction: an isothermal titration calorimetric analysis of the hammerhead ribozyme

Helical junctions are extremely common motifs in naturally occurring RNAs, but little is known about the thermodynamics that drive their folding. Studies of junction folding face several challenges: non-two-state folding behavior, superposition of secondary and tertiary structural energetics, and drastically opposing enthalpic and entropic contributions to folding. Here we describe a thermodynamic dissection of the folding of the hammerhead ribozyme, a three-way RNA helical junction, by using isothermal titration calorimetry of bimolecular RNA constructs. By using this method, we show that tertiary folding of the hammerhead core occurs with a highly unfavorable enthalpy change, and is therefore entropically driven. Furthermore, the enthalpies and heat capacities of core folding are the same whether supported by monovalent or divalent ions. These properties appear to be general to the core sequence of bimolecular hammerhead constructs. We present a model for the ion-induced folding of the hammerhead core that is similar to those advanced for the folding of much larger RNAs, involving ion-induced collapse to a structured, non-native state accompanied by rearrangement of core residues to produce the native fold. In agreement with previous enzymological and structural studies, our thermodynamic data suggest that the hammerhead structure is stabilized in vitro predominantly by diffusely bound ions. Our approach addresses several significant challenges that accompany the study of junction folding, and should prove useful in defining the thermodynamic determinants of stability in these important RNA motifs.     

A simple molecular model for thermophilic adaptation of functional nucleic acids

RNA molecules have numerous functions including catalysis and small molecule recognition, which typically arise from a tertiary structure. There is increasing interest in mechanisms for the thermostability of functional RNA molecules. Sosnick, Pan, and co-workers introduced the notion of "functional stability" as the free energy of the tertiary (functional) state relative to the next most stable (nonfunctional) state. We investigated the extent to which secondary structure stability influences the functional stability of nucleic acids. Intramolecularly folding DNA triplexes containing alternating T*AT and C+*GC base triples were used as a three-state model for the folding of nucleic acids with functional tertiary structures. A four-base-pair tunable region was included adjacent to the triplex-forming portion of the helix to allow secondary structure strength to be modulated. The degree of folding cooperativity was controlled by pH, with high cooperativity maintained by lower pH (5.5), and no cooperativity by higher pH (7.0). We find a linear relationship between functional free energy and the free energy of the secondary structure element adjacent to tertiary interactions, but only when folding is cooperative. We translate the definition of functional stability into equations and perform simulations of the thermodynamic data, which lend support to this model. The ability to increase the melting temperature of tertiary structure by strengthening base-pairing interactions separate from tertiary interactions provides a simple means for evolving thermostability in functional RNAs.     

Structural basis for altering the stability of homologous RNAs from a mesophilic and a thermophilic bacterium

Tertiary RNA structures from thermophilic bacteria generally are more stable than their mesophilic homologs. To understand the structural basis of the increase in stability, we investigated equilibrium folding of the specificity domain (S-domain) of RNase P RNA from a mesophilic (Escherichia coli) and a thermophilic (Thermus thermophilus) bacterium. Equilibrium folding of both S-domains is described by a minimal, three-state folding scheme, U-to-I-to-N. In the I-to-N transition of the thermophilic S-domain, more structure forms and protections are stronger against T1 nuclease and hydroxyl radical reactions. Phylogenetic comparison in the context of the native structure reveals that among 39 nucleotide differences between these S-domains, 12 likely contribute to higher stability. These residues participate in extensive networks of hydrogen bonding, stacking, and metal ion coordination throughout the molecule. The thermophilic S-domain achieves higher stability by mutating strategic base pairs to G-C, decreasing surface accessibility of the native state, and increasing the amount of structure formation in the native folding transition. An E. coli S-domain mutant containing these 12 nt has the same stability and folding cooperativity as the T. thermophilus S-domain. E. coli S-domain mutants containing a subset of 4 or 6 nt have the same stability as the T. thermophilus S-domain but the same folding cooperativity as the E. coli S-domain. These results show that increasing stability can be accomplished by mutations within a local structure, but increasing folding cooperativity needs concerted changes among multiple structural units.     

Exploring local and non-local interactions for protein stability by structural motif engineering

In order to probe the relative contribution of local and non-local interactions to the thermodynamic stability of proteins, we have devised an experimental approach based on a combination of motif engineering and sequence shuffling. Candidate chain segments in an immunoglobulin V(L) domain were identified whose conformation is proposed to be dominated by non-local interactions. Locally interacting structural motifs of a different conformation were then constructed as replacements, by introducing motif consensus sequences. We find that all nine replacements we constructed systematically reduce the folding cooperativity. By comparing this destabilising effect with the folding transitions of shuffled sequences for three of these motifs, we estimate the contribution of local, native interactions to the free energy of folding. Our results suggest that local and non-local interactions contribute to stability by an approximately equal amount, but that local interactions stabilise by increasing the resistance to denaturation while non-local interactions increase folding cooperativity. The systematic loss of stability by sequence shuffling in these host-guest experiments suggests that the designed interactions indeed are present in the native state, thus consensus sequence engineering may be a useful tool in structure design, but non-local interactions must be taken into account for global stability engineering. Statistical approaches are powerful tools for engineering protein structure and stability, but an analysis based on local sequence propensities alone does not adequately represent the balance of sequence and context in protein structures.     

In vitro selection of RNAs with increased tertiary structure stability.

An in vitro selection system was devised to select RNAs based on their tertiary structural stability, independent of RNA activity. Selection studies were conducted on the P4-P6 domain from the Tetrahymena thermophila group I intron, an autonomous self-folding unit that contains several important tertiary folding motifs including the tetraloop receptor and the A-rich bulge. Partially randomized P4-P6 molecules were selected based on their ability to fold into compact structures using native gel electrophoresis in the presence of decreasing concentrations of MgCl2. After 10 rounds of the selection process, a number of sequence alterations were identified that stabilized the P4-P6 RNA. One of these, a single base deletion of C209 within the P4 helix, significantly stabilized the P4-P6 molecule and would not have been identified by an activity-based selection because of its essential role for ribozyme function. Additionally, the sequence analysis provided evidence that stabilization of secondary structure may contribute to overall tertiary stability for RNAs. This system for probing RNA structure irrespective of RNA activity allows analysis of RNA structure/function relationships by identifying nucleotides or motifs important for folding and then comparing them with RNA sequences required for function.     

Stability and cooperativity of individual tertiary contacts in RNA revealed through chemical denaturation

For proteins, understanding tertiary interactions involved in local versus global unfolding has become increasingly important for understanding the nature of the native state ensemble, the mechanisms of unfolding, and the stability of both the native and intermediate states in folding. In this work we have addressed related questions with respect to RNA structure by combining chemical denaturation and hydroxyl radical footprinting methods. We have determined unfolding isotherms for each of 26 discrete sites of protection located throughout the Tetrahymena thermophila group I ribozyme. The cooperativity of folding, m-value, and the free energy, DeltaG degrees N-U, associated with formation of each tertiary contact was determined by analysis of the isotherms. The DeltaG degrees N-U values measured in this study vary from 1.7 +/- 0.2 to 7. 6 +/- 1.2 kcal mol-1. Thus, the stability of these discrete tertiary contacts vary by almost 104. In addition, an intradomain contact and three interdomain contacts show high cooperativity (m-values of 1.1 +/- 0.2 to 1.7 +/- 0.3 kcal mol-1 M-1) indicating that these contacts exhibit global cooperatively in their folding behavior. This new approach to examining RNA stability provides an exciting comparison to our understanding of protein structure and folding mechanisms.     

Influence of specific mutations on the thermal stability of the td group I intron in vitro and on its splicing efficiency in vivo: a comparative study.

Group I introns constitute excellent systems for analyzing the relationship between RNA tertiary folding and catalysis. Within a hierarchical framework interpretation of RNA folding, secondary structure motifs subtend RNA three-dimensional (3D) architecture. Thus, mutations in two-dimensional motifs are expected to have effects different from those disrupting 3D contacts. Using UV spectroscopy, we have studied the influence of nucleotide substitutions, in both secondary and tertiary structure elements, on the thermal stability of the tertiary folding of the bacteriophage T4 td group I intron. Further, we present a quantitative analysis of the relationship between the splicing efficiency in vivo and the stability of the intron structure as monitored by UV melting curves. We conclude that the stability of the tertiary structure of a group I intron as measured by UV melting is generally a good indication of its ability to splice in vivo.     

High-Pressure NMR and SAXS Reveals How Capping Modulates Folding Cooperativity of the pp32 Leucine-rich Repeat Protein

Many repeat proteins contain capping motifs, which serve to shield the hydrophobic core from solvent and maintain structural integrity. While the role of capping motifs in enhancing the stability and structural integrity of repeat proteins is well documented, their contribution to folding cooperativity is not. Here we examined the role of capping motifs in defining the folding cooperativity of the leucine-rich repeat protein, pp32, by monitoring the pressure- and urea-induced unfolding of an N-terminal capping motif (N-cap) deletion mutant, pp32-?N-cap, and a C-terminal capping motif destabilization mutant pp32-Y131F/D146L, using residue-specific NMR and small-angle X-ray scattering. Destabilization of the C-terminal capping motif resulted in higher cooperativity for the unfolding transition compared to wild-type pp32, as these mutations render the stability of the C-terminus similar to that of the rest of the protein. In contrast, deletion of the N-cap led to strong deviation from two-state unfolding. In both urea- and pressure-induced unfolding, residues in repeats 1–3 of pp32-ΔN-cap lost their native structure first, while the C-terminal half was more stable. The residue-specific free energy changes in all regions of pp32-ΔN-cap were larger in urea compared to high pressure, indicating a less cooperative destabilization by pressure. Moreover, in contrast to complete structural disruption of pp32-ΔN-cap at high urea concentration, its pressure unfolded state remained compact. The contrasting effects of the capping motifs on folding cooperativity arise from the differential local stabilities of pp32, whereas the contrasting effects of pressure and urea on the pp32-ΔN-cap variant arise from their distinct mechanisms of action.     

Effect of an alternative disulfide bond on the structure, stability, and folding of human lysozyme

Human lysozyme has four disulfide bonds, one of which, Cys65-Cys81, is included in a long loop of the beta-domain. A cysteine-scanning mutagenesis in which the position of Cys65 was shifted within a continuous segment from positions 61 to 67, with fixed Cys81, has previously shown that only the mutant W64CC65A, which has a nonnative Cys64-Cys81 disulfide, can be correctly folded and secreted by yeast. Here, using the W64CC65A mutant, we investigated the effects of an alternative disulfide bond on the structure, stability, and folding of human lysozyme using circular dichroism (CD) and fluorescence spectroscopy combined with a stopped-flow technique. Although the mutant is expected to have a different main-chain structure from that of the wild-type protein around the loop region, far- and near-UV CD spectra show that the native state of the mutant has tightly packed side chains and secondary structure similar to that of the wild-type. Guanidine hydrochloride-induced equilibrium unfolding transition of the mutant is reversible, showing high stability and cooperativity of folding. In the kinetic folding reaction, both proteins accumulate a similar burst-phase intermediate having pronounced secondary structure within the dead time of the measurement and fold into the native structure by means of a similar folding mechanism. Both the kinetic refolding and unfolding reactions of the mutant protein are faster than those of the wild-type, but the increase in the unfolding rate is larger than that of the refolding rate. The Gibbs' free-energy diagrams obtained from the kinetic analysis suggest that the structure around the loop region in the beta-domain of human lysozyme is formed after the transition state of folding, and thus, the effect of the alternative disulfide bond on the structure, stability, and folding of human lysozyme appears mainly in the native state.     

Evidence that folding of an RNA tetraloop hairpin is less cooperative than its DNA counterpart

Hairpin secondary structural elements play important roles in the folding and function of RNA and DNA molecules. Previous work from our lab on small DNA hairpin loop motifs, d(cGNAg) and d(cGNABg) (where B is C, G, or T), showed that folding is highly cooperative and obeys indirect coupling, consistent with a concerted transition. Herein, we investigate folding of the related, exceptionally stable RNA hairpin motif, r(cGNRAg) (where R is A or G). Previous NMR characterization identified a complex network of seven hydrogen bonds in this loop. We inserted three carbon (C3) spacers throughout the loop and found coupling between G1 of the loop and the CG closing base pair, similar to that found in DNA. These data support a GNRA motif being expandable at any position but before the G. Thermodynamic measurements of nucleotide-analogue-substituted oligonucleotides revealed pairwise-coupling free energies ranging from weak to strong. When coupling free energies were remeasured in the background of changes at a third site, they remained essentially unchanged even though all of the sites were coupled to each other. This type of coupling, referred to as "direct", is peculiar to the RNA loop. The data suggest that, for small stable loops, folding of RNA obeys a model with nearest-neighbor interactions, while folding of DNA follows a more concerted process in which the stabilizing interactions are linked through a conformational change. The lesser cooperativity in RNA loops may provide a more robust loop that can withstand mutations without a severe loss in stability. These differences may enhance the ability of RNA to evolve.     

Folding cooperativity in RNA and DNA is dependent on position in the helix

Secondary structural motifs play essential roles in the folding and function of RNA and DNA molecules. Previous work from our lab compared the folding of small DNA and RNA hairpin loops containing a sheared GA pair [Moody, E. M., Feerar, J. C., and Bevilacqua, P. C. (2004) Biochemistry 43, 7992-7998]. We found that the small DNA hairpins fold in a highly cooperative manner with indirect coupling, while their RNA counterparts fold in a much less cooperative fashion and display direct coupling. Herein, we extend this study to the double-stranded helix. We carried out double mutant cycles on base pairs having identical nearest-neighbor contexts but located in either external or internal helical registers. In the external register, both RNA and DNA exhibit extensive folding cooperativity between the penultimate and terminal base pair, which is independent of mismatch identity. In contrast, DNA exhibits virtually no folding cooperativity in the center of the helix, while RNA maintains substantial coupling, which is dependent on mismatch identity. Two models account for these non-nearest-neighbor effects: one involves the unfavorable entropy of helix initiation common to DNA and RNA, and the other involves steric and electrostatic strain peculiar to RNA. These data show that RNA can display cooperativity less than, greater than, or equal to that of DNA depending on context and position.     

A collapsed state functions to self-chaperone RNA folding into a native ribonucleoprotein complex

Most large RNAs achieve their active, native structures only as complexes with one or more cofactor proteins. By varying the Mg(2+) concentration, the catalytic core of the bI5 group I intron RNA can be manipulated into one of three states, expanded, collapsed or native, or into balanced equilibria between these states. Under near-physiological conditions, the bI5 RNA folds rapidly to a collapsed but non-native state. Hydroxyl radical footprinting demonstrates that assembly with the CBP2 protein cofactor chases the RNA from the collapsed state to the native state. In contrast, CBP2 also binds to the RNA in the expanded state to form many non-native interactions. This structural picture is reinforced by functional splicing experiments showing that RNA in an expanded state forms a non-productive, kinetically trapped complex with CBP2. Thus, rapid folding to the collapsed state functions to self-chaperone bI5 RNA folding by preventing premature interaction with its protein cofactor. This productive, self-chaperoning role for RNA collapsed states may be especially important to avert misassembly of large multi-component RNA-protein machines in the cell.     

Concerted kinetic folding of a multidomain ribozyme with a disrupted loop-receptor interaction

The free energy landscape for the folding of large, multidomain RNAs is rugged, and kinetically trapped, misfolded intermediates are a hallmark of RNA folding reactions. Here, we examine the role of a native loop-receptor interaction in determining the ruggedness of the energy landscape for folding of the Tetrahymena ribozyme. We demonstrate a progressive smoothing of the energy landscape for ribozyme folding as the strength of the loop-receptor interaction is reduced. Remarkably, with the most severe mutation, global folding is more rapid than for the wild-type ribozyme and proceeds in a concerted fashion without the accumulation of long-lived kinetic intermediates. The results demonstrate that a complex interplay between native tertiary interactions, divalent ion concentration, and non-native secondary structure determines the ruggedness of the energy landscape. Furthermore, the results suggest that kinetic folding transitions involving large regions of highly structured RNAs can proceed in a concerted fashion, in the absence of significant stable, preorganized tertiary structure.     

Effects of magnesium ions on the stabilization of RNA oligomers of defined structures

Optical melting was used to determine the stabilities of 11 small RNA oligomers of defined secondary structure as a function of magnesium ion concentration. The oligomers included helices composed of Watson-Crick base pairs, GA tandem base pairs, GU tandem base pairs, and loop E motifs (both eubacterial and eukaryotic). The effect of magnesium ion concentration on stability was interpreted in terms of two simple models. The first assumes an uptake of metal ion upon duplex formation. The second assumes nonspecific electrostatic attraction of metal ions to the RNA oligomer. For all oligomers, except the eubacterial loop E, the data could best be interpreted as nonspecific binding of metal ions to the RNAs. The effect of magnesium ions on the stability of the eubacterial loop E was distinct from that seen with the other oligomers in two ways. First, the extent of stabilization by magnesium ions (as measured by either change in melting temperature or free energy) was three times greater than that observed for the other helical oligomers. Second, the presence of magnesium ions produces a doubling of the enthalpy for the melting transition. These results indicate that magnesium ion stabilizes the eubacterial loop E sequence by chelating the RNA specifically. Further, these results on a rather small system shed light on the large enthalpy changes observed upon thermal unfolding of large RNAs like group I introns. It is suggested that parts of those large enthalpy changes observed in the folding of RNAs may be assigned to variations in the hydration states and types of coordinating atoms in some specifically bound magnesium ions and to an increase in the observed cooperativity of the folding transition due to the binding of those magnesium ions coupling the two stems together. Brownian dynamic simulations, carried out to visualize the metal ion binding sites, reveal rather delocalized ionic densities in all oligomers, except for the eubacterial loop E, in which precisely located ion densities were previously calculated.     

Probing possible downhill folding: native contact topology likely places a significant constraint on the folding cooperativity of proteins with approximately 40 residues

Experiments point to appreciable variations in folding cooperativity among natural proteins with approximately 40 residues, indicating that the behaviors of these proteins are valuable for delineating the contributing factors to cooperative folding. To explore the role of native topology in a protein's propensity to fold cooperatively and how native topology might constrain the degree of cooperativity achievable by a given set of physical interactions, we compared folding/unfolding kinetics simulated using three classes of native-centric C^α chain models with different interaction schemes. The approach was applied to two homologous 45-residue fragments from the peripheral subunit-binding domain family and a 39-residue fragment of the N-terminal domain of ribosomal protein L9. Free-energy profiles as functions of native contact number were computed to assess the heights of thermodynamic barriers to folding. In addition, chevron plots of folding/unfolding rates were constructed as functions of native stability to facilitate comparison with available experimental data. Although common Gō-like models with pairwise Lennard-Jones-type interactions generally fold less cooperatively than real proteins, the rank ordering of cooperativity predicted by these models is consistent with experiment for the proteins investigated, showing increasing folding cooperativity with increasing nonlocality of a protein's native contacts. Models that account for water-expulsion (desolvation) barriers and models with many-body (nonadditive) interactions generally entail higher degrees of folding cooperativity indicated by more linear model chevron plots, but the rank ordering of cooperativity remains unchanged. A robust, experimentally valid rank ordering of model folding cooperativity independent of the multiple native-centric interaction schemes tested here argues that native topology places significant constraints on how cooperatively a protein can fold.     

RNA tertiary interactions mediate native collapse of a bacterial group I ribozyme

Large RNAs collapse into compact intermediates in the presence of counterions before folding to the native state. We previously found that collapse of a bacterial group I ribozyme correlates with the formation of helices within the ribozyme core, but occurs at Mg2+ concentrations too low to support stable tertiary structure and catalytic activity. Here, using small-angle X-ray scattering, we show that Mg2+-induced collapse is a cooperative folding transition that can be fit by a two-state model. The Mg2+ dependence of collapse is similar to the Mg2+ dependence of helix assembly measured by partial ribonuclease T1 digestion and of an unfolding transition measured by UV hypochromicity. The correspondence between multiple probes of RNA structure further supports a two-state model. A mutation that disrupts tertiary contacts between the L9 tetraloop and its helical receptor destabilized the compact state by 0.8 kcal/mol, while mutations in the central triplex were less destabilizing. These results show that native tertiary interactions stabilize the compact folding intermediates under conditions in which the RNA backbone remains accessible to solvent.     

Mg2+-dependent compaction and folding of yeast tRNAPhe and the catalytic domain of the B. subtilis RNase P RNA determined by small-angle X-ray scattering

We apply synchrotron-based small-angle X-ray scattering to investigate the relationship between compaction, metal binding, and structure formation of two RNAs at 37 degrees C: the 76 nucleotide yeast tRNA(Phe) and the 255 nucleotide catalytic domain of the Bacillus subtilis RNase P RNA. For both RNAs, this method provides direct evidence for the population of a distinct folding intermediate. The relative compaction between the intermediate and the native state does not correlate with the size of the RNA but does correlate well with the amount of surface burial as quantified previously by the urea-dependent m-value. The total compaction process can be described in two major stages. Starting from a completely unfolded state (4-8 M urea, no Mg(2+)), the major amount of compaction occurs upon the dilution of the denaturant and the addition of micromolar amounts of Mg(2+) to form the intermediate. The native state forms in a single transition from the intermediate state upon cooperative binding of three to four Mg(2+) ions. The characterization of this intermediate by small-angle X-ray scattering lends strong support for the cooperative Mg(2+)-binding model to describe the stability of a tertiary RNA.     

Molecular crowders and cosolutes promote folding cooperativity of RNA under physiological ionic conditions

Folding mechanisms of functional RNAs under idealized in vitro conditions of dilute solution and high ionic strength have been well studied. Comparatively little is known, however, about mechanisms for folding of RNA in vivo where Mg²⁺ ion concentrations are low, K⁺ concentrations are modest, and concentrations of macromolecular crowders and low-molecular-weight cosolutes are high. Herein, we apply a combination of biophysical and structure mapping techniques to tRNA to elucidate thermodynamic and functional principles that govern RNA folding under in vivo–like conditions. We show by thermal denaturation and SHAPE studies that tRNA folding cooperativity increases in physiologically low concentrations of Mg²⁺ (0.5–2 mM) and K⁺ (140 mM) if the solution is supplemented with physiological amounts (∼20%) of a water-soluble neutral macromolecular crowding agent such as PEG or dextran. Low-molecular-weight cosolutes show varying effects on tRNA folding cooperativity, increasing or decreasing it based on the identity of the cosolute. For those additives that increase folding cooperativity, the gain is manifested in sharpened two-state-like folding transitions for full-length tRNA over its secondary structural elements. Temperature-dependent SHAPE experiments in the absence and presence of crowders and cosolutes reveal extent of cooperative folding of tRNA on a nucleotide basis and are consistent with the melting studies. Mechanistically, crowding agents appear to promote cooperativity by stabilizing tertiary structure, while those low molecular cosolutes that promote cooperativity stabilize tertiary structure and/or destabilize secondary structure. Cooperative folding of functional RNA under physiological-like conditions parallels the behavior of many proteins and has implications for cellular RNA folding kinetics and evolution.