Previous Article |
Table of Contents
| Next Article
BIOLOGICAL SCIENCES / BIOPHYSICS
Real-time control of the energy landscape by force directs the folding of RNA molecules
,

Departments of *Chemistry and ¶Physics and Molecular and Cell Biology, ||Howard Hughes Medical Institute, University of California, Berkeley, CA 94720
Contributed by Ignacio Tinoco, Jr., March 9, 2007 (received for review January 30, 2007)
| Abstract |
|---|
|
|
|---|
mechanical force | misfolding | optical tweezers | RNA folding | single molecule
The promiscuous nature of secondary interactions in RNA begs answers for three questions: How do RNAs find their native structure? How do they avoid the kinetic traps on their rugged folding-energy landscape? And, how can RNAs switch between alternate structures in response to changes in the cellular environment, as happens, for example, during transcription attenuation in bacteria (7)? To address these questions, we need to follow the folding trajectories of RNA molecules as they diffuse over their energy landscape in search of stable conformations. Because of the structural polymorphism and manifold pathways of RNA, single-molecule approaches (8–12) are particularly useful for answering these questions by identifying folding intermediates and distinguishing native pathways from those that misfold.
Here we use optical tweezers to unfold and control the refolding of single RNA molecules by force (11, 13). Two micrometer-sized beads are attached to the ends of an RNA molecule and are manipulated by optical tweezers to exert force on the RNA (Fig. 1a). The unfolding and refolding of the RNA molecule, induced by controlled application of force, is monitored by nanometer changes in the end-to-end distance of the molecule. In a typical "pulling" experiment, force is increased or decreased at a fixed rate (the loading and unloading rate, respectively). Significantly, the applied tension slows the hairpin folding kinetics from an order of 106 s–1 at zero force (14) to a more measurable range from 102 to 10–2 s–1 at its transition force such that the process of the folding can be monitored in real time. The force (F)-dependent un/refolding rate constant (kF) can be described by an Arrhenius-type equation (15):
|
|
where k0 is the apparent rate constant of the mechanical unfolding process at zero force, and kBT is Boltzmann's constant times temperature in Kelvin. X
, the distance to the transition state, is positive for unfolding and negative for refolding. Hence, we can control the un/refolding rate constants by manipulating the tension on a molecule.
|
To explore the possibility of alternative folding pathways, individual TAR molecules were stretched and relaxed repeatedly at approximately fixed loading and unloading rates (pN/s). In a typical cycle, the force was first ramped down from 20 to <2 pN (Fig. 1c, red) before it was increased again to 20 pN (blue). These experiments revealed three different types of trajectories, each associated with distinctive refolding characteristics: single-step refolding, multiple-step refolding, and misfolding. The single-step refolding mostly occurs at force relaxation rates <1 pN/s. In this case the force-extension curve displays a single transition (zip) at 10–13 pN with a decrease in extension of
17 nm. This behavior is consistent with a single-step refolding process from a single strand to the native structure. When pulled, the natively folded TAR always displayed a sudden length increase of
18 nm that appeared as a single rip at
13–16 pN. In contrast, when the tension on the molecule was allowed to relax at a rate of 1.5 pN/s, >95% of the refolding curves show multiple-step refolding, displaying a decreasing extension with characteristic fluctuations (Fig. 1c, black arrows) followed by a small zip. This zip corresponds to the folding of an average of 26 nt, which is half the value for the refolding of the entire molecule. Subsequent pulling of such molecules displayed unfolding characteristics similar to those observed in the slow relaxation curves, indicating that these RNAs had folded into the native structure. The back-and-forth oscillations in the force-extension curve before the zip are likely successive refolding and unfolding events as the molecule moves into and out of alternative, competing folding pathways. Once an RNA molecule, following this search process, folds into intermediate structures that contain only native contacts, the remaining helix can form cooperatively, as manifested by the observed zip.
At unloading rates >2 pN/s, we observed another type of relaxation curve that showed no zipping at all; at low force the extension of the RNA never returned to that of the native hairpin (Fig. 1c, misfolding). When the force on the molecule was subsequently increased, the extension of the RNA increased monotonically, retracing the relaxation curve until a zip transition suddenly shortened the end-to-end distance of the molecule, against the increasing force. This transition, not observed previously, is a refolding event in that the extension of the molecule is shortened, yet it occurs as the force is increasing, which normally causes unfolding. The magnitude of this zip shows a broad distribution with a mean of
30 single-stranded nucleotides becoming paired (Fig. 2a), indicating that a range of partially folded structures existed before the transition. As the force was further raised, the force-extension curve always showed a rip similar to the unfolding of the native hairpin. We conclude that the structures at low force are misfolded, and that the application of force allowed the RNA molecules to form the native hairpin in the zip transition. Accordingly, we term this zip a "rescue" transition.
|
12 pN (19). To further rule out this possibility, we kept the misfolded structures at zero tension for 30 min. Subsequent pulls showed the characteristic zip–rip signature (in that order), clearly indicating that these nonnative structures are stable at low forces and that forces >6 pN are required for their conversion to the native hairpin. Such forces were needed to break the nonnative interactions that trapped the RNA structure, allowing the molecule to search for its correct folding pathway and effectively rescuing it from the misfolded state.
Assuming that the rate-liming step in the rescue is disruption of nonnative base pairs, we can determine the kinetics for rescue. The distribution of rescue forces was fitted to the following equation (27),
|
|
where [M(F, r)] is the fraction of misfolded molecule at each force (F), and the loading rate is (r); A is a factor reflecting both the rescue-rate constant at zero force and instrumental effects (11); and X
represents the average distance to the transition state of the native fold from the misfolded structures. We obtained X
of 6 ± 1 nm (Fig. 2d) from fitting the data shown in Fig. 2c to Eq. 2.
Which steps in the folding process of the native hairpin are perturbed by the fast force relaxation? High force favors the longer single-strand forms over base-paired structures. As the force is relaxed, new intramolecular contacts and partly folded structures become accessible to the molecule. The effect of relaxing the force, therefore, is to expand and change continuously the topography of the energy surface over which the molecules can diffuse. The folding of a hairpin involves nucleation and helix formation, both of which are influenced by force. Formation of the first base pair closing the native loop has to compete with many other possible loop-closing base pairs. When the molecule is allowed to equilibrate, either at constant force near its F1/2 or when the force is relaxed very slowly relative to the rate of forming the correct base pairs, the single strand folds into the native hairpin by sequentially base-pairing from the native loop; the refolding appears to be two-state. In contrast, when the force is relaxed fast, alternative loop closures and nonnative base pairings become accessible to the molecule before the native contacts form. As the force continues to decrease fast, some of these nonnative loop closures will not have sufficient time to open but instead lead to propagation of nonnative helices at lower forces. Such a structural evolution is reflected by the fluctuations displayed in the multiple-step refolding trajectories. Probabilities of both forming and unfolding of these misfolded structures depend on the instantaneous value of the force. When the force is decreased faster than the rate at which misfolded structures can unfold, the molecule becomes kinetically trapped in metastable states.
To reduce or avoid misfolding, the topography of the energy landscape needs to be modified in such a way that the channels leading to the native fold become more favorable than alternate pathways. This situation can be achieved by engineering the molecular structure such that the correct folding is speeded up and/or formation of the misfolded structures is slowed. To test these hypotheses, we measured misfolding of two mutant TAR hairpins. First, we tried to facilitate the correct folding. In the native structure (Fig. 1b), formation of the bulge-closing base pair is energetically unfavorable (28) and is a slow step in the helix formation. To eliminate this slow step, we constructed a mutant hairpin, TAR
b, in which the 3-nt bulge is deleted. As predicted, the bulgeless mutant is more stable than TAR, with an F1/2 of
16 pN. Moreover, TAR
b folds faster than TAR by four orders of magnitude (Fig. 3d) and shows no misfolding (Fig. 4).
|
|
|
How can a misfolded RNA molecule escape from the kinetic trap and find its minimum energy structure? To gain insight into this question, we monitored the structural evolution of misfolded RNAs to the native fold. First, force was decreased on an unfolded molecule from 20 pN to between 7 and 9.5 pN at a rate of
11 pN/s to induce misfolding; in nearly half of the trajectories, the molecule was not completely folded when the set force was reached. Then, the force was held constant until the native fold was reached, which was indicated by a decrease in extension (Fig. 3a). Typically, the extension of the misfolded structures remained relatively constant for a few seconds until a zipping event occurred (Fig. 3b). The decrease in extension of this zip,
Xrescue2, shows a broad distribution (Fig. 3c), indicating that many different partially folded structures, with an average of 10 ± 3 bp, had formed after the fast force relaxation. At the rescue forces, the single strand would fold into the native hairpin in milliseconds or less, yet lifetimes of these partially folded structures are at least tens of seconds (Fig. 3d), indicating that these structures must be off-pathway, misfolded species. After the zipping, the extension varied little up to 5 min, indicating that the native hairpin is stable at this force and that the mechanical rescue of the native structure from the metastable structures is irreversible. Subsequent pulling verified the native folding (Fig. 3a, gray curve).
The lifetimes of the misfolded structures at each force fit a single exponential distribution, from which we calculated the rescue-rate constant at force F, krescue2(F). The force dependence of the rescue-rate constant is fitted to Eq. 1, in which X
represents the average distance to the transition state of the native fold from the misfolded structures. The rescue rate increases with the force, yielding a positive value of X
(4.9 ± 0.8 nm; Fig. 3d), characteristic of a helix-unfolding transition (30). This value also agrees with the X
obtained from pulling experiments (Fig. 2d). Rescue below 7 pN took at least a few minutes and was rarely observed, consistent with the rescue-force distribution observed in the force-ramp experiments (Fig. 2b). However, between 9.5 and 11 pN, both the native folding rate and the rescue rate are so fast that misfolded structures are practically not observable.
In a classic study of RNA misfolding, the activity of a misfolded tRNA was restored by heating the sample in the presence of Mg2+ (6). Presumably, raising the temperature, just as increasing the force, disrupts the nonnative conformation and thereby allows correct folding. Here, the rescue rate was raised by increasing the force applied to the molecule, confirming that the rate-limiting step in the rescue is the disruption of incorrect, trapped structures. An important question, then, is: Must the trapped molecule unfold completely to gain access to the native fold? If, during rescue, the misfolded RNAs were first unfolded into a single strand, the extension would increase by
10 nm on average. However, such an unfolding step was not observed before the rescue zip (Fig. 3). To further investigate the nature of the molecular rearrangements that take place during rescue, we characterized the fluctuations of the molecular extension in the constant-force rescue experiments. The root-mean-square deviation of the molecular extension before the appearance of the zip is 0.7 ± 0.2 nm, just slightly higher than the value observed after the native hairpin formed (0.5 ± 0.1 nm). The lack of a distinctive unfolding step suggests that a series of partial unfolding and refolding events in small steps are responsible for gradually converting the RNA into an on-pathway intermediate, at which time the zipping transition occurred. It is notable that in the multiple-step refolding observed at high relaxation rates, many small-step unfolding/refolding transitions occur before the final zip (Fig. 1c). Taken together these observations suggest that the structural rearrangements that lead the molecule from misfolded structures to the native state occur through pathways made up of multiple small unfolding/refolding steps and not through a path that requires the complete unfolding of the RNA. A similar mechanism has been proposed for the conformational switch of a ribozyme (31).
We have shown how direct mechanical manipulation can be used to control the folding pathway of a single molecule. In particular, we showed that decreasing the force applied to the molecule can be used to change, in real time, the number of configurations accessible to the molecule; lowering the force increases the number of structures and interactions accessible to the molecule. Moreover, by controlling the rate at which the force is relaxed, we effectively modify the folding-energy landscape (15, 21, 32) and the alternative pathways over which the molecule diffuses in its search for the native state. The relaxation rate, then, is the control parameter that determines how many and how fast these pathways become accessible to the molecule, the degree of kinetic competition between these pathways, and, just as importantly, how long the molecule is allowed to follow a particular pathway before other competing paths become available.
As shown here in the case of TAR, by direct manipulation and modification of the force applied to the molecule, it is possible to suddenly "tilt" the energy landscape of the molecule, increasing its chance to attain off-pathway states and even allow the molecule to switch back to the on-pathway amid the folding. We know of no other experimental method that permits such control over the detailed kinetics of refolding of a molecule in real time. Combined with mutagenesis, we have characterized here what the crucial elements are in TAR that favor misfolding (i.e., the three-base bulge and the sequence at the 5' end). More elaborate schemes of force relaxation interspersed with constant force events could be used to characterize the molecular processes leading to various intermediates, to trapping, and to structural rearrangement of the RNA.
To maintain cellular homeostasis, living cells must ensure the proper folding of millions of copies of RNAs. Deletion of RNA chaperones has been implicated in cell death after UV irradiation and in autoimmune disease (33). Mounting evidence, accumulated in recent years, suggests that RNA folding in vivo is an assisted process through the participation of specialized proteins (34–40). The mechanism by which RNA chaperone proteins facilitate correct RNA folding remains unclear (3, 41). The results presented above suggest two possible scenarios. Proteins may simply facilitate the attainment of the folded structure by binding to specific regions of the unfolded structure and sterically interfering with the formation of alternative competing pathways on the folding-energy surface. Alternatively, motor proteins, such as helicases (17), generate and apply force to actively induce the mechanical disruption of RNA molecules trapped in nonnative folds in the cell. The methods described here make it possible to investigate the mechanisms used by proteins or other molecular components to steer and facilitate the attainment of the native fold through the rugged RNA landscapes.
| Materials and Methods |
|---|
|
|
|---|
500-bp dsDNA handles corresponding to the regions flanking the hairpin were amplified by PCRs and modified by either biotin or digoxigenin at each of the termini. The DNA handles were annealed to the RNA. Properly annealed molecules can be tethered between two microspheres coated with streptavidin and anti-digoxigenin antibody, respectively (11) (Fig. 1a). Optical Tweezers. Dual-beam optical tweezers (13) were used to study the RNA folding. The streptavidin-coated bead was placed on a micropipette mounted on a piezoelectric flexture stage. The anti-digoxigenin antibody-coated bead was held by the optical trap. Force is determined by the displacement of the trapped bead from the center of the trap. The change in the extension of the molecule results from the movement of both the trapped bead and the bead on the micropipette.
Manipulation of Single Molecules. In the force-ramp experiments, the piezoelectric stage moved the micropipette at a constant rate (nm/s), which changed the tension on the RNA with an approximately constant (un)loading rate (pN/s) between 3 and 20 pN. When constant force was implemented, a feedback mechanism using a proportional, integrative, and differential algorithm was used to maintain the tension. In all experiments, force and extension were recorded at 100 Hz.
Folding Condition. The folding was studied in 10 mM Hepes (pH 8.0), 100 mM KCl, and 0.1 mM EDTA; 0.05% NaN3 was added to prevent the growth of bacteria. All folding reactions were studied at 22°C.
Conversion of
X.
To compare the increase in end-to-end extension,
X, at different forces, each
X is converted to the number of single-stranded ribonucleotides at the corresponding force. The conversion is calculated by using the worm-like-chain interpolation formula (42) with a single-stranded RNA contour length of 0.59 nm/nt and persistent length of 1 nm.
| Acknowledgements |
|---|
|
|
|---|
| Footnotes |
|---|
Abbreviations: TAR, transactivation response region.
To whom correspondence may be addressed. E-mail: intinoco{at}lbl.gov or panli{at}albany.edu
Author contributions: P.T.X.L., C.B., and I.T. designed research; P.T.X.L. performed research; P.T.X.L. and I.T. analyzed data; and P.T.X.L., C.B., and I.T. wrote the paper.
Present address: Department of Biological Sciences, University at Albany, State University of New York, Albany, NY 12222. ![]()
The authors declare no conflict of interest.
© 2007 by The National Academy of Sciences of the USA
| References |
|---|
|
|
|---|
This article has been cited by other articles in HighWire Press-hosted journals:
![]() |
Q. Peng and H. Li Atomic force microscopy reveals parallel mechanical unfolding pathways of T4 lysozyme: Evidence for a kinetic partitioning mechanism PNAS, February 12, 2008; 105(6): 1885 - 1890. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Mickler, R. I. Dima, H. Dietz, C. Hyeon, D. Thirumalai, and M. Rief Revealing the bifurcation in the unfolding pathways of GFP by using single-molecule experiments and simulations PNAS, December 18, 2007; 104(51): 20268 - 20273. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Chen, J.-D. Wen, and I. Tinoco Jr Single-molecule mechanical unfolding and folding of a pseudoknot in human telomerase RNA RNA, December 1, 2007; 13(12): 2175 - 2188. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||