A folding control element for tertiary collapse of a group II intron ribozyme

Ribozymes derived from the group II intron ai5gamma collapse to a compact intermediate, folding to the native state through a slow, direct pathway that is unperturbed by kinetic traps. Molecular collapse of ribozyme D135 requires high magnesium concentrations and is thought to involve a structural element in domain 1 (D1). We used nucleotide analog interference mapping, in combination with nondenaturing gel electrophoresis, to identify RNA substructures and functional groups that are essential for D135 tertiary collapse. This revealed that the most crucial atoms for compaction are located within a small section of D1 that includes the kappa and zeta elements. This small substructure controls specific collapse of the molecule and, in later steps of the folding pathway, it forms the docking site for catalytic D5. In this way, the stage is set for proper active site formation during the earliest steps of ribozyme folding.     

An obligate intermediate along the slow folding pathway of a group II intron ribozyme

Most RNA molecules collapse rapidly and reach the native state through a pathway that contains numerous traps and unproductive intermediates. The D135 group II intron ribozyme is unusual in that it can fold slowly and directly to the native state, despite its large size and structural complexity. Here we use hydroxyl radical footprinting and native gel analysis to monitor the timescale of tertiary structure collapse and to detect the presence of obligate intermediates along the folding pathway of D135. We find that structural collapse and native folding of Domain 1 precede assembly of the entire ribozyme, indicating that D1 contains an on-pathway intermediate to folding of the D135 ribozyme. Subsequent docking of Domains 3 and 5, for which D1 provides a preorganized scaffold, appears to be very fast and independent of one another. In contrast to other RNAs, the D135 ribozyme undergoes slow tertiary collapse to a compacted state, with a rate constant that is also limited by the formation D1. These findings provide a new paradigm for RNA folding and they underscore the diversity of RNA biophysical behaviors.     

Structural Rearrangements Linked to Global Folding Pathways of the Azoarcus Group I Ribozyme

Stable RNAs must fold into specific three-dimensional structures to be biologically active, yet many RNAs form metastable structures that compete with the native state. Our previous time-resolved footprinting experiments showed that Azoarcus group I ribozyme forms its tertiary structure rapidly (τ < 30 ms) without becoming significantly trapped in kinetic intermediates. Here, we use stopped-flow fluorescence spectroscopy to probe the global folding kinetics of a ribozyme containing 2-aminopurine in the loop of P9. The modified ribozyme was catalytically active and exhibited two equilibrium folding transitions centered at 0.3 and 1.6 mM Mg²⁺, consistent with previous results. Stopped-flow fluorescence revealed four kinetic folding transitions with observed rate constants of 100, 34, 1, and 0.1 s^− 1 at 37 °C. From comparison with time-resolved Fe(II)-ethylenediaminetetraacetic acid footprinting of the modified ribozyme under the same conditions, these folding transitions were assigned to formation of the I_C intermediate, tertiary folding and docking of the nicked P9 tetraloop, reorganization of the P3 pseudoknot, and refolding of nonnative conformers, respectively. The footprinting results show that 50-60% of the modified ribozyme folds in less than 30 ms, while the rest of the RNA population undergoes slow structural rearrangements that control the global folding rate. The results show how small perturbations to the structure of the RNA, such as a nick in P9, populate kinetic folding intermediates that are not observed in the natural ribozyme.     

Concerted folding of a Candida ribozyme into the catalytically active structure posterior to a rapid RNA compaction

Folding of the major population of Tetrahymena intron RNA into the catalytically active structure is trapped in a slow pathway. In this report, folding of Candida albicans intron was investigated using the trans-acting Ca.L-11 ribozyme as a model. We demonstrated that both the catalytic activity (k_obs) and compact folding equilibrium of Ca.L-11 are strongly dependent on Mg²⁺ at physiological concentrations, with both showing an Mg²⁺ Hill coefficient of 3. Formation of the compact structure of Ca.L-11 is shown to occur very rapidly, on a subsecond time scale similar to that of RNase T1 cleavage. Most of the ribozyme RNA population folds into the catalytically active structure with a rate constant of 2 min^–1 at 10 mM Mg²⁺; neither slower kinetics nor obvious Mg²⁺ inhibition is observed. These results suggest that folding of the Ca.L-11 ribozyme is initiated by a rapid magnesium-dependent RNA compaction, which is followed by a slower searching for the native contacts to form the catalytically active structure without interference from the long-lived trapped states. This model thus provides an ideal system to address a range of interesting aspects of RNA folding, such as conformational searching, ion binding and the role of productive intermediates.     

Group II intron folding under near-physiological conditions: collapsing to the near-native state

The folding of group II intron ribozymes has been studied extensively under optimal conditions for self-splicing in vitro (42 degrees C and high magnesium ion concentrations). In these cases, the ribozymes fold directly to the native state by an apparent two-state mechanism involving the formation of an obligate intermediate within intron domain 1. We have now characterized the folding pathway under near-physiological conditions. We observe that compaction of the RNA proceeds slowly to completion, even at low magnesium concentration (3 mM). Kinetic analysis shows that this compact species is a "near-native" intermediate state that is readily chased into the native state by the addition of high salt. Structural probing reveals that the near-native state represents a compact domain 1 scaffold that is not yet docked with the catalytic domains (D3 and D5). Interestingly, native ribozyme reverts to the near-native state upon reduction in magnesium concentration. Therefore, while the intron can sustain the intermediate state under physiological conditions, the native structure is not maintained and is likely to require stabilization by protein cofactors in vivo.     

Hinge stiffness is a barrier to RNA folding

Cation-mediated RNA folding from extended to compact, biologically active conformations relies on a temporal balance of forces. The Mg^2 +-mediated folding of the Tetrahymena thermophila ribozyme is characterized by rapid nonspecific collapse followed by tertiary-contact-induced compaction. This article focuses on an autonomously folding portion of the Tetrahymena ribozyme, its P4-P6 domain, in order to probe one facet of the rapid collapse: chain flexibility. The time evolution of P4-P6 folding was followed by global and local measures as a function of Mg^2 + concentration. While all concentrations of Mg^2 + studied are sufficient to screen the charge on the helices, the rates of compaction and tertiary contact formation diverge as the concentration of Mg^2 + increases; collapse is greatly accelerated by Mg^2 +, while tertiary contact formation is not. These studies highlight the importance of chain stiffness to RNA folding; at 10 mM Mg^2 +, a stiff hinge limits the rate of P4-P6 folding. At higher magnesium concentrations, the rate-limiting step shifts from hinge bending to tertiary contact formation.     

Monovalent and divalent promoted GAAA tetraloop-receptor tertiary interactions from freely diffusing single-molecule studies

Proper assembly of RNA into catalytically active three-dimensional structures requires multiple tertiary binding interactions, individual characterization of which is crucial to a detailed understanding of global RNA folding. This work focuses on single-molecule fluorescence studies of freely diffusing RNA constructs that isolate the GAAA tetraloop-receptor tertiary interaction. Freely diffusing conformational dynamics are explored as a function of Mg²⁺ and Na⁺ concentration, both of which promote facile docking, but with 500-fold different affinities. Systematic shifts in mean fluorescence resonance energy transfer efficiency values and line widths with increasing [Na⁺] are observed for the undocked species and can be interpreted with a Debye model in terms of electrostatic relaxation and increased flexibility in the RNA. Furthermore, we identify a 34 ± 2% fraction of freely diffusing RNA constructs remaining undocked even at saturating [Mg²⁺] levels, which agrees quantitatively with the 32 ± 1% fraction previously reported for immobilized constructs. This verifies that the kinetic heterogeneity observed in the docking rates is not the result of surface tethering. Finally, the K_D value and Hill coefficient for [Mg²⁺]-dependent docking decrease significantly for [Na⁺] = 25 mM vs. 125 mM, indicating Mg²⁺ and Na⁺ synergy in the RNA folding process.     

Mg2+-dependent folding of a Diels-Alderase ribozyme probed by single-molecule FRET analysis

Here, we report a single-molecule fluorescence resonance energy transfer (FRET) study of a Diels-Alderase (DAse) ribozyme, a 49-mer RNA with true catalytic properties. The DAse ribozyme was labeled with Cy3 and Cy5 as a FRET pair of dyes to observe intramolecular folding, which is a prerequisite for its recognition and turnover of two organic substrate molecules. FRET efficiency histograms and kinetic data were taken on a large number of surface-immobilized ribozyme molecules as a function of the Mg²⁺ concentration in the buffer solution. From these data, three separate states of the DAse ribozyme can be distinguished, the unfolded (U), intermediate (I) and folded (F) states. A thermodynamic model was developed to quantitatively analyze the dependence of these states on the Mg²⁺ concentration. The FRET data also provide information on structural properties. The I state shows a strongly cooperative compaction with increasing Mg²⁺ concentration that arises from association with several Mg²⁺ ions. This transition is followed by a second Mg²⁺-dependent cooperative transition to the F state. The observation of conformational heterogeneity and continuous fluctuations between the I and F states on the ~100ms timescale offers insight into the folding dynamics of this ribozyme.     

An alternative route for the folding of large RNAs: apparent two-state folding by a group II intron ribozyme

Despite a growing literature on the folding of RNA, our understanding of tertiary folding in large RNAs derives from studies on a small set of molecular examples, with primary focus on group I introns and RNase P RNA. To broaden the scope of RNA folding models and to better understand group II intron function, we have examined the tertiary folding of a ribozyme (D135) that is derived from the self-splicing ai5gamma intron from yeast mitochondria. The D135 ribozyme folds homogeneously and cooperatively into a compact, well-defined tertiary structure that includes all regions critical for active-site organization and substrate recognition. When D135 was treated with increasing concentrations of Mg(2+) and then subjected to hydroxyl radical footprinting, similar Mg(2+) dependencies were seen for internalization of all regions of the molecule, suggesting a highly cooperative folding behavior. In this work, we show that global folding and compaction of the molecule have the same magnesium dependence as the local folding previously observed. Furthermore, urea denaturation studies indicate highly cooperative unfolding of the ribozyme that is governed by thermodynamic parameters similar to those for forward folding. In fact, D135 folds homogeneously and cooperatively from the unfolded state to its native, active structure, thereby demonstrating functional reversibility in RNA folding. Taken together, the data are consistent with two-state folding of the D135 ribozyme, which is surprising given the size and multi-domain structure of the RNA. The findings establish that the accumulation of stable intermediates prior to formation of the native state is not a universal feature of RNA folding and that there is an alternative paradigm in which the folding landscape is relatively smooth, lacking rugged features that obstruct folding to the native state.     

Different Folding Pathways Taken by Highly Homologous Proteins, Goat α-Lactalbumin and Canine Milk Lysozyme

Is the folding pathway conserved in homologous proteins? To address this question, we compared the folding pathways of goat α-lactalbumin and canine milk lysozyme using equilibrium and kinetic circular dichroism spectroscopy. Both Ca²⁺-binding proteins have 41% sequence identity and essentially identical backbone structures. The Φ-value analysis, based on the effect of Ca²⁺ on the folding kinetics, showed that the Ca²⁺-binding site was well organized in the transition state in α-lactalbumin, although it was not yet organized in lysozyme. Equilibrium unfolding and hydrogen-exchange 2D NMR analysis of the molten globule intermediate also showed that different regions were stabilized in the two proteins. In α-lactalbumin, the Ca²⁺-binding site and the C-helix were weakly organized, whereas the A- and B-helices, both distant from the Ca²⁺-binding site, were well organized in lysozyme. The results thus provide an example of highly homologous proteins taking different folding pathways. To understand the molecular origin of this difference, we investigated the native three-dimensional structures of the proteins in terms of non-local contact clusters, a parameter based on the residue-residue contact map and known to be well correlated with the folding rate of non-two-state proteins. There were remarkable differences between the proteins in the distribution of the non-local contact clusters, and these differences provided a reasonable explanation of the observed difference in the folding initiation sites. In conclusion, the protein folding pathway is determined not only by the backbone topology but also by the specific side-chain interactions of contacting residues.     

Single-Molecule Studies of the Im7 Folding Landscape

Under appropriate conditions, the four-helical Im7 (immunity protein 7) folds from an ensemble of unfolded conformers to a highly compact native state via an on-pathway intermediate. Here, we investigate the unfolded, intermediate, and native states populated during folding using diffusion single-pair fluorescence resonance energy transfer by measuring the efficiency of energy transfer (or proximity or P ratio) between pairs of fluorophores introduced into the side chains of cysteine residues placed in the center of helices 1 and 4, 1 and 3, or 2 and 4. We show that while the native states of each variant give rise to a single narrow distribution with high P values, the distributions of the intermediates trapped at equilibrium (denoted I^eqm) are fitted by two Gaussian distributions. Modulation of the folding conditions from those that stabilize the intermediate to those that destabilize the intermediate enabled the distribution of lower P value to be assigned to the population of the unfolded ensemble in equilibrium with the intermediate state. The reduced stability of the I^eqm variants allowed analysis of the effect of denaturant concentration on the compaction and breadth of the unfolded state ensemble to be quantified from 0 to 6 M urea. Significant compaction is observed as the concentration of urea is decreased in both the presence and absence of sodium sulfate, as previously reported for a variety of proteins. In the presence of Na₂SO₄ in 0 M urea, the P value of the unfolded state ensemble approaches that of the native state. Concurrent with compaction, the ensemble displays increased peak width of P values, possibly reflecting a reduction in the rate of conformational exchange among iso-energetic unfolded, but compact conformations. The results provide new insights into the initial stages of folding of Im7 and suggest that the unfolded state is highly conformationally constrained at the outset of folding.     

Characterization of a folding intermediate from HIV-1 ribonuclease H.

The RNase H domain from HIV-1 (HIV RNase H) encodes an essential retroviral activity. Refolding of the isolated HIV RNase H domain shows a kinetic intermediate detectable by stopped-flow far UV circular dichroism and pulse-labeling H/D exchange. In this intermediate, strands 1, 4, and 5 as well as helices A and D appear to be structured. Compared to its homolog from Escherichia coli, the rate limiting step in refolding of HIV RNase H appears closer to the native state. We have modeled this kinetic intermediate using a C-terminal deletion fragment lacking helix E. Like the kinetic intermediate, this variant folds rapidly and shows a decrease in stability. We propose that inhibition of the docking of helix E to this folding intermediate may present a novel strategy for anti HIV-1 therapy.     

Thermal unfolding of apo- and holo-enolase from Saccharomyces cerevisiae: different mechanisms, similar activation enthalpies

Yeast enolase is stabilized by its natural cofactor Mg²⁺. This stabilization is ascribed to the reduced subunit dissociation of the holoprotein. Nevertheless, how Mg²⁺ alters the unfolding mechanism has yet to be fully characterized. Here, we investigate the role of Mg²⁺ in the denaturation mechanism and unfolding kinetics of yeast enolase. Apo-enolase unfolds through a three-state process (N₂ ↔ 2I → 2D). The intermediate species is described as a monomeric molten globule-like conformation that becomes noticeable in the presence of phosphate and is able to recover its native secondary structure when cooled down. Kinetic studies confirmed the presence of the intermediate species, even though it was not noticeable in the thermal scans. The cofactor increases the cooperativity of the unfolding transitions, while the intermediate species becomes less noticeable or nonexistent. Thus, holo-enolase follows a simple two-state mechanism (N₂ → 2D). Our results indicate smaller unfolding rate-constants in the presence of Mg²⁺, thus favoring the native state. The temperature dependence of the unfolding rates allowed us to calculate the activation enthalpies of denaturation. Interestingly, despite the different unfolding mechanisms of the apo and holo forms of enolase, they both have similar activation barriers of denaturation (185-190 kJ mol⁻¹).     

The structural stabilization of the κ three-way junction by Mg(II) represents the first step in the folding of a group II intron

Folding of group II introns is characterized by a first slow compaction of domain 1 (D1) followed by the rapid docking of other domains to this scaffold. D1 compaction initiates in a small subregion encompassing the κ and ζ elements. These two tertiary elements are also the major interaction sites with domain 5 to form the catalytic core. Here, we provide the first characterization of the structure adopted at an early folding step and show that the folding control element can be narrowed down to the three-way junction with the κ motif. In our nuclear magnetic resonance studies of this substructure derived from the yeast mitochondrial group II intron Sc.ai5γ, we show that a high affinity Mg(II) ion stabilizes the κ element and enables coaxial stacking between helices d′ and d′′, favoring a rigid duplex across the three-way junction. The κ-element folds into a stable GAAA-tetraloop motif and engages in A-minor interactions with helix d′. The addition of cobalt(III)hexammine reveals three distinct binding sites. The Mg(II)-promoted structural rearrangement and rigidification of the D1 core can be identified as the first micro-step of D1 folding.     

Extended Structures in RNA Folding Intermediates Are Due to Nonnative Interactions Rather than Electrostatic Repulsion

RNA folding occurs via a series of transitions between metastable intermediate states for Mg²⁺ concentrations below those needed to fold the native structure. In general, these folding intermediates are considerably less compact than their respective native states. Our previous work demonstrates that the major equilibrium intermediate of the 154-residue specificity domain (S-domain) of the Bacillus subtilis RNase P RNA is more extended than its native structure. We now investigate two models with falsifiable predictions regarding the origins of the extended intermediate structures in the S-domains of the B. subtilis and the Escherichia coli RNase P RNA that belong to different classes of P RNA and have distinct native structures. The first model explores the contribution of electrostatic repulsion, while the second model probes specific interactions in the core of the folding intermediate. Using small-angle X-ray scattering and Langevin dynamics simulations, we show that electrostatics plays only a minor role, whereas specific interactions largely account for the extended nature of the intermediate. Structural contacts in the core, including a nonnative base pair, help to stabilize the intermediate conformation. We conclude that RNA folding intermediates adopt extended conformations due to short-range, nonnative interactions rather than generic electrostatic repulsion of helical domains. These principles apply to other ribozymes and riboswitches that undergo functionally relevant conformational changes.     

Diffusive motions control the folding and unfolding kinetics of the apomyoglobin pH 4 molten globule intermediate

The sperm whale apomyoglobin pH 4 folding intermediate exists in two forms, Ia and Ib, that mimic transient kinetic intermediates in the folding of the native protein at pH 6. To characterize the nature of the kinetic barrier that controls the formation of the earliest intermediate Ia, we have investigated the effects of small viscogenic cosolvents on its folding and unfolding kinetics. The kinetics are measurable by stopped-flow fluorescence and follow a cooperative two-state model in the absence and presence of cosolvents. Small cosolvents stabilize Ia, but, by applying the isostability test to separate the viscogenic effect of the cosolvent from its stabilizing effect, we found that, in both folding and unfolding conditions, the apparent rate constant decreases when solvent viscosity increases. The unitary inverse dependence of the apparent rate constant on solvent viscosity indicates a diffusion-controlled reaction. This result is consistent with the hypothesis that folding of the apomyoglobin pH 4 intermediate obeys a diffusion-collision model. Additionally, the temperature dependence of the reaction rate at constant viscosity indicates that the formation of Ia is also controlled by an energy barrier. Linear free energy relationships show that the transition state of the U <==> Ia reaction is compact and buries 45% of the surface area that is buried in native apomyoglobin. We conclude that the transition state of the U <==> Ia reaction resembles that for the formation of native proteins; namely, it is dry and its compactness is closer to that of the folded (Ia) form than of the unfolded form.     

Conserved folding pathways of alpha-lactalbumin and lysozyme revealed by kinetic CD, fluorescence, NMR, and interrupted refolding experiments

In this report, it is shown by a combination of stopped-flow CD, fluorescence, and time-resolved NMR studies that the Ca^2 +-induced refolding of bovine α-lactalbumin (BLA) at constant denaturant concentration (4 M urea) exhibits triple-exponential kinetics. In order to distinguish between parallel folding pathways and a strictly sequential formation of the native state, interrupted refolding experiments were conducted. We show here that the Ca^2 +-induced refolding of BLA involves parallel pathways and the transient formation of a folding intermediate on the millisecond timescale. Our data furthermore suggest that the two structurally homologous proteins BLA and hen egg white lysozyme share a common folding mechanism. We provide evidence that the guiding role of long-range interactions in the unfolded state of lysozyme in mediating intersubdomain interactions during folding is replaced in the case of BLA by the Ca^2 +-binding site. Time-resolved NMR spectroscopy, in combination with fast ion release from caged compounds, enables the measurement of complex protein folding kinetics at protein concentrations as low as 100 μM and the concomitant detection of conformational transitions with rate constants of up to 8 s^− 1.     

Understanding the Role of Three-Dimensional Topology in Determining the?Folding Intermediates of Group I Introns

Many RNA molecules exert their biological function only after folding to unique three-dimensional structures. For long, noncoding RNA molecules, the complexity of finding the native topology can be a major impediment to correct folding to the biologically active structure. An RNA molecule may fold to a near-native structure but not be able to continue to the correct structure due to a topological barrier such as crossed strands or incorrectly stacked helices. Achieving the native conformation thus requires unfolding and refolding, resulting in a long-lived intermediate. We investigate the role of topology in the folding of two phylogenetically related catalytic group I introns, the Twort and Azoarcus group I ribozymes. The kinetic models describing the Mg²⁺-mediated folding of these ribozymes were previously determined by time-resolved hydroxyl (⋅OH) radical footprinting. Two intermediates formed by parallel intermediates were resolved for each RNA. These data and analytical ultracentrifugation compaction analyses are used herein to constrain coarse-grained models of these folding intermediates as we investigate the role of nonnative topology in dictating the lifetime of the intermediates. Starting from an ensemble of unfolded conformations, we folded the RNA molecules by progressively adding native constraints to subdomains of the RNA defined by the ⋅OH time-progress curves to simulate folding through the different kinetic pathways. We find that nonnative topologies (arrangement of helices) occur frequently in the folding simulations despite using only native constraints to drive the reaction, and that the initial conformation, rather than the folding pathway, is the major determinant of whether the RNA adopts nonnative topology during folding. From these analyses we conclude that biases in the initial conformation likely determine the relative flux through parallel RNA folding pathways.     

Structural and functional analysis of the phosphoryl transfer reaction mediated by the human small C‐terminal domain phosphatase,Scp1

Human small C‐terminal domain phosphatase 1 (Scp1) modulates the phosphorylation state of the C‐terminal domain (CTD) of eukaryotic RNA polymerase II (RNAP II), with preference for phosphorylated Ser5 in the tandem heptad repeats of the CTD. Additionally, Scp1 was identified as a conserved regulator of neuronal stem cell development. Scp1 is a member of haloacid dehalogenase (HAD) superfamily, whose catalysis depends on a Mg²⁺ ion and a DXDX(T/V) motif. The first Asp of the motif is identified as the nucleophile that is subject to phosphorylation leading to a phosphoryl‐aspartate intermediate. This high‐energy mixed anhydride intermediate is subsequently hydrolyzed to regenerate the enzyme. In the present study, we successfully captured the phosphoryl‐aspartate intermediate in the crystal structure of a Scp1D206A mutant soaked with para‐nitrophenyl phosphate (pNPP), providing strong evidence for the proposed mechanism. Furthermore, steady‐state kinetic analysis of a variety of Scp1 mutants revealed the importance of Asp206 in Mg²⁺ coordination mediated by a water molecule. Overall, we captured the snapshots of the phosphoryl transfer reaction at each stage of Scp1‐mediated catalysis. Through structural‐based sequence alignment, we show that the spatial position of the D206 side chain is strictly conserved throughout HAD family. Our results strongly suggest that Asp206 and its equivalent residues in other HAD family members play important structural and possible mechanistic roles.