This Article Abstract Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Gerland, U. Articles by Hwa, T. Search for Related Content PubMed PubMed Citation Articles by Gerland, U. Articles by Hwa, T.

Mechanically Probing the Folding Pathway of Single RNA Molecules

Ulrich Gerland ^*, Ralf Bundschuh  and Terence Hwa ^*

^* Department of Physics, University of California at San Diego, La Jolla, California 92093-0319; and Department of Physics, The Ohio State University, Columbus, Ohio 43210-1106

Correspondence: Address reprint requests to Ulrich Gerland, 9500 Gilman Dr., La Jolla, CA 92093 USA. Tel.: 858-534-7256; Fax: 858-534-7697; E-mail: gerland{at}physics.ucsd.edu.

	ABSTRACT

TOP ABSTRACT INTRODUCTION MODEL RESULTS SUMMARY AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
We study theoretically the denaturation of single RNA molecules by mechanical stretching, focusing on signatures of the (un)folding pathway in molecular fluctuations. Our model describes the interactions between nucleotides by incorporating the experimentally determined free energy rules for RNA secondary structure, whereas exterior single-stranded regions are modeled as freely jointed chains. For exemplary RNA sequences (hairpins and the Tetrahymena thermophila group I intron), we compute the quasiequilibrium fluctuations in the end-to-end distance as the molecule is unfolded by pulling on opposite ends. Unlike the average quasiequilibrium force-extension curves, these fluctuations reveal clear signatures from the unfolding of individual structural elements. We find that the resolution of these signatures depends on the spring constant of the force-measuring device, with an optimal value intermediate between very rigid and very soft. We compare and relate our results to recent experiments by Liphardt et al. (2001).

	INTRODUCTION

TOP ABSTRACT INTRODUCTION MODEL RESULTS SUMMARY AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
The advent and refinement of techniques to apply and measure forces on single molecules have allowed remarkable experiments probing the elastic properties and structural transitions of biomolecules such as DNA (see, e.g., Bockelmann et al., 1998; Maier et al., 2000; Smith et al., 1996) and proteins (see, e.g., Rief et al., 1997). Recently, single-molecule experiments using optical tweezers were also performed to study the unfolding of small RNA molecules under an applied force (Liphardt et al., 2001, 2002). RNA molecules, besides functioning as messengers of sequence information in the process of protein synthesis, also have many biological functions that intricately depend on a precisely folded RNA structure, e.g., as ribosomal RNA or as self-splicing introns (Cech, 1993). The formation of RNA structure (RNA folding) is therefore an important biophysical process, which, however, is currently not understood in sufficient detail (Tinoco and Bustamante, 1999). Pulling experiments (Liphardt et al., 2001, 2002), possibly in combination with single-molecule fluorescence methods (Zhuang et al., 2000, 2002), promise to reveal important aspects of RNA structure, its folding pathways and kinetics, and eventually its biological function.

The experiment of Liphardt et al. (2001) has shown that pulling on simple structural units of RNA, e.g., a single hairpin, yields characteristic features in the force-extension curve (FEC), which can be used to deduce the unfolding free energy and the size of the structural element. Moreover, the observation of end-to-end distance time traces at constant force revealed that the unfolding and refolding of these simple structures proceeds directly without intermediates, and the opening/closing rates could be extracted from the time traces. From the theoretical side, it is interesting to ask how far, in principle, the pulling approach could be pushed to study larger RNA molecules and how the resolution of such approaches depends on parameters of the experimental setup. Here, we expand our previous model for RNA pulling experiments (Gerland et al., 2001), and address these questions. The model incorporates the experimentally known free energy rules for RNA secondary structure (Walter et al., 1994) and a polymer model for the elastic properties of single-stranded RNA (ssRNA), but neglects pseudoknots and tertiary interactions, the energetics of which are currently poorly characterized.

Our model yields predictions for force-extension measurements, including fluctuations, and the mechanical (un)folding pathway for any given RNA sequence. Below, we begin with the P5ab hairpin used by Liphardt et al. (2001) and demonstrate that our model yields an FEC that is in semiquantitative agreement with the experimental curve. This agreement gives us confidence that our model is sufficiently realistic to permit its use to explore general questions regarding sequence-dependent signatures in mechanical single-molecule experiments on RNA.

We first address the question of intermediates in the unfolding pathway and show that, according to our model, a small modification in the sequence of the P5ab hairpin can change its two-state folding behavior and introduce a locally dominant intermediate. However, whether this intermediate state can be observed through quasiequilibrium fluctuation measurements critically depends on the experimental conditions: If the force-measuring device is a soft spring, the hairpin unfolds without any visible intermediates. Only when the force is measured with a stiff spring can an intermediate state be observed.

We then consider a larger RNA molecule, the group I intron of Tetrahymena thermophila with a sequence of 400 bases and a known secondary structure containing many individual elements (Cech, 1993). Previous theoretical work predicted that equilibrium FECs of large RNAs are smooth and display no secondary structure dependent features, due to a compensation effect between individual structural elements (Gerland et al., 2001). Here, we find that even in the absence of structure-based features in the FEC, measurements of the equilibrium fluctuations of the entire molecule can still be useful to obtain information on the (un)folding pathway. Also, for this longer molecule, a good choice of the stiffness of the force-measuring device is important for the determination of the pathway.

Finally, we compare the mechanical unfolding process studied here with the more conventional thermal unfolding. Again using the group I intron as an example, we find that the individual structural elements display significantly sharper opening transitions for force-induced than for thermal denaturation within our model. Indeed, in UV absorption experiments (Banerjee et al., 1993), thermal melting of the intron shows only one broad signature associated with the opening of the secondary structure.

	MODEL

TOP ABSTRACT INTRODUCTION MODEL RESULTS SUMMARY AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
We consider an experimental setup where the two ends of an RNA molecule are attached to a force- and extension-measuring device, e.g., an atomic force microscope or beads trapped by optical tweezers. For the sake of concreteness, we use the optical tweezer setup sketched in Fig. 1 a as an example here. In an idealized theoretical model, the RNA molecule (in essence, a highly nonlinear elastic element for the present context) is connected in series with a linear harmonic spring, and the total extension is externally controlled (see Fig. 1 b). The linear spring models the force-measuring device. E.g., for optical tweezers, its harmonic potential approximates the potential of the optical trap, and the spring constant is adjustable in the experiment through the laser intensity. The extension of the spring, R_s, corresponds to the position of the bead with respect to the minimum of the trapping potential, and is measured as a function of the total extension R_t. The average force acting on the RNA molecule and its average extension are

(1)

(2)

Here, ... denotes a thermal average over all accessible conformations of the RNA molecule and the spring at fixed total extension R_t. Throughout this article, we assume that the pulling experiment is performed slowly enough, such that the RNA molecule can fully sample its conformational space on timescales that are short compared to the external timescale of the imposed stretching process. This limit corresponds to the quasiequilibrium regime, where the work dissipated in the pulling process is negligible. The average ... then includes an average over all possible secondary structures of the molecule.

View larger version (28K): [in this window] [in a new window]   FIGURE 1  Sketch of the system considered. (a) In a typical experimental setup, the two ends of an RNA molecule are attached to beads, the position of which is controlled, e.g., with optical tweezers. (b) In a theoretical model, the potential of the optical trap can be approximated by a linear spring (with spring constant ). This spring is connected in series with a highly nonlinear elastic element (the RNA molecule, usually with an additional linker). The total extension, Rt, is externally controlled, whereas the individual extensions, Rs and R, undergo thermal fluctuations.

A model for an experimental setup such as the one sketched in Fig. 1 a requires two separate parts, one that contains only information on the free energies for the secondary structures of the RNA sequence, and another that describes the polymer properties of ssRNA (Gerland et al., 2001). The coupling between these two parts is through the total length of the exterior ssRNA segments of the molecule, i.e., the number m of bases that actually "feel" the applied force (see Fig. 1 a). For every fixed total extension R_t, the system has to find a compromise between lowering the RNA binding energy by decreasing m, and gaining entropy (plus lowering elastic energy) by increasing m. The secondary structure part is given by the partition function, Q(m), summed over all secondary structures of the RNA molecule with a fixed number m of exterior single-stranded bases. (We account for the width of stems within the exterior single strand by increasing the base count m by three for each stem.) The polymer properties of ssRNA then enter through a function W_tot(R_t;m), which denotes the total end-to-end distance distribution of an ssRNA molecule of m bases in series with the spring. (For simplicity, we restrict the spring to be collinear with the end-to-end distance of the RNA molecule.) The total partition function, Z(R_t), is then simply the convolution of Q(m) with W_tot(R_t;m),

(3)

The detailed calculation of Q(m) is described in Gerland et al. (2001). Briefly, the calculation takes the experimentally determined rules for the binding free energies of RNA secondary structures (Walter et al., 1994) into account, but neglects tertiary structure effects and pseudoknots. Our method is based on a recursive method to calculate the partition function (McCaskill, 1990) as implemented in the "Vienna RNA Package" (Hofacker et al., 1994), but introduces an additional recursion relation for the calculation of Q(m).

The total end-to-end distance distribution of an ssRNA molecule of m bases in series with the spring can be expressed as

(4)

where W_RNA(R;m) denotes the distribution of the ssRNA molecule alone. We model the ssRNA as an elastic freely jointed chain, i.e., as a chain of segments with unrestricted relative orientations, whereas the segment length r is constrained by a harmonic potential V(r) = (r - b)²/2. The average segment length b corresponds to the Kuhn length of a noninteracting ssRNA chain. The number of segments of the elastic freely jointed chain used to represent an RNA molecule with m bases is chosen as ml/b where l is the base-to-base distance of ssRNA. This yields where C is a normalization constant, and h is determined from (Gerland et al., 2001). Such a model has been shown to describe the elastic properties of ssDNA molecules (Maier et al., 2000; Montanari and Mézard, 2001). Since we are not aware of the corresponding data for RNA, we use the DNA values for the polymer parameters as obtained by Montanari and Mézard (2001) through fitting to the experiment of Maier et al. (2000), i.e., l = 0.7 nm, b = 1.9 nm, and (/k_BT)^-1/2 = 0.1 nm (we do not expect a large difference in the single-strand properties of DNA and RNA because of the high similarity between their chemical structures). We take the free energy parameters for RNA secondary structure as supplied with the Vienna package (version 1.3.1) at room temperature T = 25°C. The salt concentrations at which the free energy parameters were measured are [Na⁺] = 1 M and [Mg²⁺] = 0 M.

Observables
From the partition function (Eq. 3), we get the total free energy F(R_t) = -k_BT log Z(R_t), the average force f = F(R_t)/R_t, and through Eqs. 1 and 2 also the average extension R of the RNA molecule. The equilibrium fluctuations in the extension of the RNA molecule are given by (R)² = R² - R² = [ - ²F/R_t²]/ß². Using these relations, one can show (the derivation is straightforward after observing that that the variance is related to the derivative of the FEC as measured in the mixed ensemble through

(5)

This exact relation simply expresses the intuitive result that the extension of the molecule undergoes large thermal fluctuations in the flat regions of the FEC.

To study the unfolding pathway, we introduce the weight x_m(R_t) for the secondary structures with m exterior open bases at a given value of R_t,

(6)

which satisfy the normalization condition The binding probability of base i and j at given R_t may then be expressed as

(7)

where p_ij(m) denotes the basepairing probability in the ensemble of structures with fixed m, which can be calculated as described previously (Gerland et al., 2001).

	RESULTS

TOP ABSTRACT INTRODUCTION MODEL RESULTS SUMMARY AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
We first consider the P5ab hairpin, which is the simplest RNA molecule studied in the experiment of Liphardt et al. (2001), and is of a size that typically occurs as a structural subunit in larger RNAs. Fig. 2 a shows its native structure and the experimental FEC for P5ab with added linker, taken with permission from Liphardt et al. (2001). The superimposed solid line shows the FEC that results from our model as described above with a spring constant of = 0.2 pN/nm (the native structure of P5ab contains two nonstandard G-A basepairs for which no quantitative information on the binding energies is available. We therefore follow the suggestion of Liphardt et al. (2001) and approximate the binding energy of G-A basepairs by replacing them with G-U pairs). Here, we modeled the RNA-DNA hybrid linker, which in the experiment connects the RNA molecule to the beads, as a wormlike chain (WLC) with an experimentally known length of L = 320 nm. We used the WLC interpolation formula (Bustamante et al., 1994)

(8)

where the value for the persistence length, l_p = 3.57 nm, and the (experimentally unknown) origin of the distance scale were obtained by fitting to the experimental FEC at low forces. The shape of the experimental FEC is well captured by the theoretical curve, with a characteristic "hump" indicating the opening of the hairpin (the force drop in the hump occurs, because the released single-strand creates "slack" in the polymer). However, the force at which the hairpin opens is overestimated by the theoretical model. To assess the origin of this discrepancy, we applied a simple salt correction to the RNA-binding (free) energies (SantaLucia, 1998) to account for the different ionic conditions between the single-molecule experiment ([Na⁺] = 0.25 M) and the standard conditions ([Na⁺] = 1.0 M). The resulting FEC (dashed line in Fig. 2 a) yields near-quantitative agreement with the experiment, indicating that the difference in ionic conditions accounts for most of the discrepancy. The remaining discrepancy could be caused, for instance, by the freely jointed chain model for ssRNA or the unknown stacking energies for G-A basepairs. However, the good agreement we have achieved seems sufficient to justify the use of our model as a tool to study the general questions outlined in the Introduction (in the following, we keep the salt concentration at the standard conditions).

View larger version (23K): [in this window] [in a new window]   FIGURE 2  Force-extension curve for the P5ab hairpin with added linker. (a) Native structure of the hairpin, its FEC as experimentally determined by Liphardt et al. (2001), and the FEC resulting from our model for two different salt concentrations. (b) FEC for P5ab as results from our model for two different spring constants (solid and dashed lines; see main text for details). The dotted line represents the FEC for a nonbinding sequence of the same length as P5ab. (Inset) End-to-end distance fluctuations as a function of the average extension.

We now seek to characterize the unfolding process of the P5ab hairpin within our model. In the experiment of Liphardt et al. (2001), the time traces of the end-to-end extension at fixed forces around 14 pN show the characteristic pattern of a bistable system, i.e., the unfolding proceeds directly, without intermediates. To calculate theoretical time traces would require a full kinetic treatment of the RNA folding/unfolding process, which is beyond the scope of the present article. Instead, we determine the amplitude of the equilibrium fluctuations in the end-to-end extension, R, as a function of the mean force or extension. Experimentally, this amplitude would be determined as a time-average, where R(t) denotes the recorded time trace of the end-to-end extension and its average. Our theoretical prediction is for the infinite time limit, T , of this average. Note that the exact relation Eq. 5 directly links the equilibrium fluctuations to the derivative of the averaged FEC. In practice, where the averaging period T is finite and limited to a few minutes by instrumental drift (Liphardt et al., 2001), the amplitude R can be determined with much greater accuracy than the derivative of the time-averaged FEC, since taking the derivative strongly amplifies the statistical error. Therefore, the quasiequilibrium FEC and the equilibrium fluctuations, R, can effectively be regarded as two independent and complementary sources of information.

The inset of Fig. 2 b shows the fluctuations R as a function of the average extension, R, in the range of extensions over which the hairpin opens. We observe a single peak of the fluctuations at the force/extension where the transition takes place, which is consistent with the two-state behavior found experimentally. Physically, this peak is caused by continual kinetic fluctuations between the open and the closed state of the hairpin.

Hairpin with intermediate
Next, we address the question whether the simple two-state behavior observed for P5ab is a generic property for hairpins of this typical size. The discussion of this question also serves as a preparation for our ensuing study of a larger RNA with more complicated structure. To this end, we modify the sequence of the P5ab hairpin slightly by replacing the two G-A basepairs with U's on both strands, which leads to a small interior loop in the hairpin. This change does not significantly affect the FEC (see Fig. 3 a); however the fluctuations R take on a qualitatively different behavior, as shown in Fig. 3 b: For a soft spring, the R curve shows only a single peak as before, but a stiff spring yields two maxima with a pronounced minimum in between. At this minimum, the configurational distribution of the molecule is localized on an intermediate state, which strongly reduces the fluctuations in the extension. This is demonstrated explicitly in Fig. 4, which shows the probability distribution in m-space, i.e., x_m as given by Eq. 6 (see figure legend for details).

View larger version (14K): [in this window] [in a new window]   FIGURE 3  Calculated force-extension characteristics for the modified P5ab hairpin (G-A basepairs replaced by A-A) for a hard ( = 0.2 pN/nm) and a soft ( = 0.01 pN/nm) harmonic potential of the optical tweezer. (a) The average FEC, which resembles very much the FEC in Fig. 2 b. (b) The fluctuations of the extension as a function of the average extension. Although the curve is very similar to the inset of Fig. 2 b for a soft spring, it develops a pronounced minimum corresponding to an intermediate state for a hard spring. The arrows indicate the positions at which the xm distributions shown in Fig. 4 were calculated.

View larger version (36K): [in this window] [in a new window]   FIGURE 4  Probability distributions in m-space for the modified P5ab hairpin. The left column is calculated with a soft spring ( = 0.01 pN/nm), and the right column with a stiff spring ( = 0.2 pN/nm). The curves in the top, middle, and bottom rows are calculated for different average extensions R corresponding to the positions indicated by the arrows in Fig. 3 b (first peak, minimum, and second peak in the fluctuation curve, respectively). For the stiff spring, the distribution at the minimum is almost entirely localized on the intermediate state. In contrast, the distribution is spread out over all three states for the soft spring at the same average extension.

Intuitively, the fact that the appearance of the intermediate state depends on the spring constant is best understood by mapping the unzipping problem onto the problem of a particle with an attached spring in a random potential and coupled to a thermal bath (see Fig. 5). As the other end of the spring is moved in one direction, the particle performs a "stick-slip" motion (Bockelmann et al., 1998). With a stiff spring, the particle hops over short distances and tends to be localized, whereas the particle can make long jumps, sliding over many valleys, when the spring is soft.

View larger version (23K): [in this window] [in a new window]   FIGURE 5  Particle with an attached spring in the free energy landscape of the modified P5ab hairpin. The other end position of the spring is externally controlled (see Gerland et al., 2001, for the details of the mapping to the RNA unzipping problem). The spring in this picture represents the measurement device, the linker, and the ssRNA, and has an effective spring constant Its harmonic potential for = 0.01 pN/nm and = 0.2 pN/nm is shown as the dashed and dash-dotted lines, respectively. With the soft spring, the particle can jump from the closed to the open state immediately as soon as the force is strong enough to overcome the energy barrier. The hard spring has a steeper harmonic, potential and thus the particle will first jump into the intermediate state and later into the open state as the end of the spring is moved toward larger m.

View larger version (33K): [in this window] [in a new window]   FIGURE 6  Secondary structure of the Tetrahymena thermophila group I intron as obtained from the Vienna package. The overall structure, and almost all individual basepairs, are identical to the one of the inactive conformation as determined by Pan and Woodson (1998) (we follow the labeling introduced there).

How can one determine the (un)folding pathway in a single-molecule experiment? As we have seen, pulling on an isolated subunit produces a characteristic signature in the FEC (see Fig. 2). If every structural subunit in a larger RNA molecule would produce such a characteristic signature, one could directly read off the folding pathway from the FEC. However, our previous theoretical study predicted that the quasiequilibrium FEC of RNAs with many subunits will typically show no distinguishable signatures of individual subunits opening (Gerland et al., 2001). This prediction is consistent with the experimental observation of smooth FECs for long ssDNA molecules (Maier et al., 2000). Three physical effects are responsible for the smoothness: i), the "floppiness" of the external single strand and the linkers, ii), thermal fluctuations in the secondary structure, i.e., the contribution of suboptimal structures, and most importantly, iii), the fact that changes in the extension of individual structural subunits can compensate each other, since only the total end-to-end distance is measured (Gerland et al., 2001). As a result, quasiequilibrium FECs cannot be used to study the folding pathway (or the secondary structure) of larger RNA molecules.

The smoothing mechanisms (ii and iii) can be strongly suppressed by performing the experiments at large pulling speeds, which do not leave sufficient time for the molecule to sample the entire ensemble of different secondary structures with similar number of exterior unpaired bases, m, and free energy of folding. Hence, nonequilibrium FECs of larger RNA molecules should be very rugged and display signatures of individual structural elements opening (see also the discussion in Gerland et al., 2001). This was indeed observed recently in the Berkeley group (S. Dumont and I. Tinoco, Jr., private communication). From the theoretical side, the calculation of nonequilibrium FECs is challenging, since the simple particle-in-a-landscape picture of Fig. 5 breaks down in nonequilibrium: The free energy landscape is no longer well-defined. The secondary structures that are accessible for a molecule at a given extension and within a given time window depend on the present structure of the molecule and the detailed folding kinetics. Although important steps toward a kinetic theory of RNA folding have been taken (Isambert and Siggia, 2000), there is currently insufficient experimental information available to construct a full theoretical model.

Here, we consider instead the equilibrium fluctuations in the end-to-end distance of the entire intron, and explore, within our theoretical model, how much information on the folding pathway can be obtained from this. Although looking at the equilibrium fluctuations does not suppress any of the above-mentioned smoothing mechanisms, it does provide a greater resolution for the observation of structural transitions, as already discussed. We will see in the following that this enhanced resolution, which is dependent on the spring constant, is sufficient to reveal major steps of the unfolding process.

The quasiequilibrium FEC of the intron is shown in Fig. 7 a. Here, we have not added an additional linker in the calculation, since the intron is already fairly large and flexible, so that the effect of a comparatively short and stiff double-stranded linker is negligible. As expected, the FEC displays no signatures of the secondary structure, mostly due to "compensation" between different structural elements (Gerland et al., 2001). As for the P5ab hairpin above, we compute the equilibrium fluctuations R, which would correspond to the average standard deviation in an infinitely long experimental time trace R(t). Obviously, our averaged R contains less information than experimental R(t) traces. Nevertheless, as Fig. 7 b shows, already the equilibrium fluctuations R display interesting sequence-dependent features, i.e., local maxima and minima. The minima can again be interpreted as intermediate states along the (un)folding pathway, and the maxima as structural transitions. Note that the central peak close to R = 150 nm splits into two peaks as the stiffness of the spring is increased, indicating the appearance of a new intermediate on the (un)folding pathway. This illustrates again how the ability to resolve intermediate states depends on the stiffness of the force-measuring device. For the choice of the latter, a trade-off between resolution of intermediate structures and fluctuation amplitude has to be respected.

View larger version (25K): [in this window] [in a new window]   FIGURE 7  Calculated force-extension characteristics of the Tetrahymena thermophila group I intron. (a) FEC. (b) Fluctuations in the extension of the ribozyme for different stiffness of the force-measuring device. (c) Fractional opening for all of the structural elements at = 0.02 pN/nm. A black square indicates that the structural element is completely intact whereas a white square indicates that all basepairs of the element have been opened. It can be seen that the FEC lacks all structural features. In contrast, the fluctuations display pronounced peaks that coincide with the disappearance of specific structural elements.

By simultaneous inspection of Fig. 7, b and c, we can assign particular peaks in the fluctuation curve to the opening of particular elements. For instance, the central peak is associated with the parallel opening of P9.1, P9.1a, and P9.2, whereas the most stable elements P5b, P5c open at the final hump. Generally, elements that free a lot of single strand upon opening, such as a stem-loop structure with a large loop (e.g., P2.1), are associated with a distinct peak in the fluctuation curve. This peak reflects the large difference in extension between the open and closed state of the element.

To obtain the information displayed in Fig. 7 c in the experiment will be more laborious, but should be possible given the known secondary structure. Since at every mean extension (or force) only a small number (up to three) structural elements open in parallel, one can expect a small number of characteristic plateaus in the experimental R(t) traces. The identification of the plateaus with the corresponding elements can be obtained by cutting the molecule at several places between the known secondary structure elements and performing the same measurement on subsequences. A similar approach was taken for recent nonequilibrium experiments in the Berkeley group (S. Dumont and I. Tinoco, Jr., private communication). Obtaining the mechanical (un)folding pathway for a moderately sized RNA with known secondary structure therefore appears feasible by measuring the fluctuations.

It is interesting to compare the mechanical unfolding studied so far with thermal unfolding. The thermal unfolding process of the group I intron can be studied theoretically with the Vienna RNA Package (Hofacker et al., 1994). Fig. 8 b shows the fractional opening of all structural elements as in Fig. 7 c, but plotted against temperature. We observe that the individual stems show a less sharp transition from closed to open than in Fig. 7 c. As a result, the specific heat curve shown in Fig. 8 a displays only a single broad peak for the unfolding of the secondary structure (the specific heat is obtained from the total folding free energy, G, through -T ²(G)/T²). This is consistent with UV absorption experiments (Banerjee et al., 1993) that found only two broad peaks for the same RNA, where the low-temperature peak could be associated with melting of the tertiary structure and the high-temperature peak with the secondary structure.

View larger version (33K): [in this window] [in a new window]   FIGURE 8  Thermal unfolding of the intron. (a) Specific heat curve. (b) Fractional opening for all of the structural elements. Analogously to Fig. 7 c, a black square represents a fully formed structure whereas a white square represents a structure that has fully opened up.

	SUMMARY AND DISCUSSION

TOP ABSTRACT INTRODUCTION MODEL RESULTS SUMMARY AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Single-molecule experiments are particularly useful when the molecule under study can take on several conformations, e.g., folding intermediates, because in these cases bulk measurements usually yield only average properties, whereas single-molecule measurements can characterize each of the different conformations. Here, we predicted experimental signatures of such folding intermediates in the quasiequilibrium fluctuations of the end-to-end distance. We found that a local minimum in the fluctuations curve indicates a locally stable intermediate state in the average mechanical unfolding pathway (which is well-defined and identical to the (re)folding pathway). The number of local minima in the fluctuations curve can therefore serve as a lower estimate of the number of locally stable intermediate states in this pathway. Furthermore, these signatures could be used to reconstruct the (un)folding pathway of an RNA molecule with a known secondary structure by cutting the sequence at appropriate positions and individually probing different substructures of the molecule.

We also showed that the stiffness of the force-measuring device plays a crucial role in determining the resolution of the quasiequilibrium fluctuations curve. If the device is too soft, some local minima can be lost, whereas a too rigid device has a low fluctuation amplitude and therefore a low signal-to-noise ratio. The optimal choice for the stiffness (see Results) leads to the best estimate for the number of locally stable intermediate states along the pathway and increases the degree to which the pathway could be reconstructed with the "cutting approach". The dependence on the stiffness has a simple intuitive explanation within the statistical mechanics problem of a particle with an attached spring in a one-dimensional energy landscape (see Fig. 5): A stiff force-measuring device corresponds to a rigid spring, which tends to localize the particle in valleys of the landscape (i.e., locally stable intermediates) and forces it to make small jumps as it is pulled along the landscape. On the other hand, a soft spring allows the particle to spread out and make long jumps, effectively hiding the intermediates in the equilibrium curve.

We hope that this view, which also explains the smoothness of quasiequilibrium FECs for large RNAs (Gerland et al., 2001), will be useful for the design and interpretation of future experiments. Of course, experiments typically record not only the average fluctuations, but the detailed time traces of the end-to-end distance (Liphardt et al., 2001). These traces probe the kinetics of structural rearrangements in RNA, and could provide even more detailed resolution of intermediate states (if the kinetics is slow enough to be resolvable in the experiment). However, we expect that it is still helpful to make the choice for the stiffness of the force-measuring device as discussed above, to limit the number of states that contribute to each of the time traces.

Although single-molecule experiments on RNA have already produced remarkable results, there are many desirable future developments. For instance, it appears that the current approaches are not well suited to measure the secondary structure of an unknown molecule (e.g., in the cutting approach, one would not know where to cut and it would be too laborious to try all positions). Is there a direct way to mechanically measure the secondary structure? Also, could one detect signatures of pseudoknots or even tertiary interactions? We hope that theoretical models of the type presented here can be useful for the planning and design of new experimental approaches in the future.

	ACKNOWLEDGEMENTS

TOP ABSTRACT INTRODUCTION MODEL RESULTS SUMMARY AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Note: After completion of this work, we learned of a related theoretical study by S. Cocco, J.F. Marko, and R. Monasson (preprint, cond-mat/0207609). These authors present a kinetic model for the end-to-end distance fluctuations of stretched RNA molecules, which nicely complements our thermodynamic study. Very recently, A.F. Sauer-Budge, J.A. Nyamwanda, D.K. Lubensky, and D. Branton (preprint cond-mat/0209414) have experimentally studied the unzipping kinetics of double-stranded DNA with a short internal loop, and demonstrated the presence of an intermediate state comparable to the one discussed here for the modified P5ab hairpin.

We are grateful to S.M. Block, S. Dumont, H. Isambert, J. Liphardt, D.K. Lubensky, I. Tinoco, Jr., and S.A. Woodson for stimulating discussions.

U.G. was supported in part by a fellowship from the German Academic Exchange Service. R.B. and T.H. acknowledge support by the National Science Foundation through grant No. DMR-9971456, DBI-9970199, and the Beckmann Foundation.

Submitted on August 9, 2002; accepted for publication December 20, 2002.

	REFERENCES

TOP ABSTRACT INTRODUCTION MODEL RESULTS SUMMARY AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Banerjee, A. R., J. A. Jaeger, and D. H. Turner. 1993. Thermal unfolding of a group I ribozyme: the low-temperature transition is primarily disruption of tertiary structure. Biochemistry. 32:153–163.[Medline]

Bockelmann, U., B. Essevaz-Roulet, and F. Heslot. 1998. DNA strand separation studied by single molecule force measurements. Phys. Rev. E. 58:2386–2394.

Bustamante, C., J. F. Marko, E. D. Siggia, and S. B. Smith. 1994. Entropic elasticity of lambda-phage DNA. Science. 265:1599–1600.[Free Full Text]

Cech, T. R. 1993. Structure and mechanism of the large catalytic RNAs: group I and group II introns and ribonuclease P. In The RNA World. R. F. Gesteland and J. F. Atkins, editors. Cold Spring Harbor Laboratory Press, Plainview, NY. 239–269.

Fernandez, J. M., S. Chu, and A. F. Oberhauser. 2001. RNA structure. Pulling on hair(pins). Science. 292:653–654.[Free Full Text]

Gerland, U., R. Bundschuh, and T. Hwa. 2001. Force-induced denaturation of RNA. Biophys. J. 81:1324–1332.[Abstract/Free Full Text]

Heymann, B., and H. Grubmüller. 2000. Dynamic force spectroscopy of molecular adhesion bonds. Phys. Rev. Lett. 84:6126–6129.[Medline]

Hofacker, I. L., W. Fontana, P. F. Stadler, S. Bonhoeffer, M. Tacker, and P. Schuster. 1994. Fast folding and comparison of RNA secondary structures. Monatshefte f. Chemie. 125:167–188. Software available online at http://www.tbi.univie.ac.at/.

Isambert, H., and E. D. Siggia. 2000. Modeling RNA folding paths with pseudoknots: application to hepatitis delta virus ribozyme. Proc. Natl. Acad. Sci. USA. 97:6515–6520.[Abstract/Free Full Text]

Jaeger, A. A., M. Zuker, and D. H. Turner. 1990. Melting and chemical modification of a cyclized self-splicing group I intron: similarity of structures in 1 M Na⁺, in 10 mM Mg²⁺, and in the presence of substrate. Biochemistry. 29:10147–10158.[Medline]

Liphardt, J., B. Onoa, S. B. Smith, I. Tinoco, and C. Bustamante. 2001. Reversible unfolding of single RNA molecules by mechanical force. Science. 292:733–737.[Abstract/Free Full Text]

Liphardt, J., S. Dumont, S. B. Smith, I. Tinoco, and C. Bustamante. 2002. Equilibrium information from nonequilibrium measurements in an experimental test of Jarzynski's equality. Science. 296:1832–1835.[Abstract/Free Full Text]

Maier, B., D. Bensimon, and V. Croquette. 2000. Replication by a single DNA polymerase of a stretched single-stranded DNA. Proc. Natl. Acad. Sci. USA. 97:12002–12007.[Abstract/Free Full Text]

McCaskill, J. S. 1990. The equilibrium partition function and base pair binding probabilities for RNA secondary structure. Biopolymers. 29:1105–1119.[Medline]

Montanari, A., and M. Mézard. 2001. Hairpin formation and elongation of biomolecules. Phys. Rev. Lett. 86:2178–2181.[Medline]

Pan, J., and S. A. Woodson. 1998. Folding intermediates of a self-splicing RNA: mispairing of the catalytic core. J. Mol. Biol. 280:597–609.[Medline]

Rief, M., M. Gautel, F. Oesterhelt, J. M. Fernandez, and H. E. Gaub. 1997. Reversible unfolding of individual titin immunoglobulin domains by AFM. Science. 276:1109–1112.[Abstract/Free Full Text]

SantaLucia, J., Jr. 1998. A unified view of polymer, dumbbell, and oligonucleotide DNA nearest-neighbor thermodynamics. Proc. Natl. Acad. Sci. USA. 95:1460–1465.[Abstract/Free Full Text]

Smith, S. B., Y. Cui, and C. Bustamante. 1996. Overstretching B-DNA: the elastic response of individual double-stranded and single-stranded DNA molecules. Science. 271:795–799.[Abstract]

Tinoco, I., Jr., and C. Bustamante. 1999. How RNA folds. J. Mol. Biol. 293:271–281.[Medline]

Walter, A. E., D. H. Turner, J. Kim, M. H. Lyttle, P. Muller, D. H. Mathews, and M. Zuker. 1994. Coaxial stacking of helixes enhances binding of oligoribonucleotides and improves predictions of RNA folding. Proc. Natl. Acad. Sci. USA. 91:9218–9222.[Abstract/Free Full Text]

Zhuang, X., L. E. Bartley, H. P. Babcock, R. Russell, T. Ha, D. Herschlag, and S. Chu. 2000. A single-molecule study of RNA catalysis and folding. Science. 288:2048–2051.[Abstract/Free Full Text]

Zhuang, X., H. Kim, M. J. B. Pereira, H. P. Babcock, N. G. Walter, and S. Chu. 2002. Correlating structural dynamics and function in single ribozyme molecules. Science. 296:1473–1476.[Abstract/Free Full Text]

This article has been cited by other articles:

				C. Hyeon and D. Thirumalai Mechanical Unfolding of RNA: From Hairpins to Structures with Internal Multiloops Biophys. J., February 1, 2007; 92(3): 731 - 743. [Abstract] [Full Text] [PDF]

				N. W. Moore and T. L. Kuhl The Role of Flexible Tethers in Multiple Ligand-Receptor Bond Formation between Curved Surfaces Biophys. J., September 1, 2006; 91(5): 1675 - 1687. [Abstract] [Full Text] [PDF]

				C. Hyeon and D. Thirumalai Forced-Unfolding and Force-Quench Refolding of RNA Hairpins Biophys. J., May 15, 2006; 90(10): 3410 - 3427. [Abstract] [Full Text] [PDF]

				K. A. Walther, J. Brujic, H. Li, and J. M. Fernandez Sub-Angstrom Conformational Changes of a Single Molecule Captured by AFM Variance Analysis Biophys. J., May 15, 2006; 90(10): 3806 - 3812. [Abstract] [Full Text] [PDF]

				F. Liu, H. Tong, and Z.-c. Ou-Yang Force Unfolding Single RNAs Biophys. J., March 15, 2006; 90(6): 1895 - 1902. [Abstract] [Full Text] [PDF]

				M. Andersson, E. Fallman, B. E. Uhlin, and O. Axner A Sticky Chain Model of the Elongation and Unfolding of Escherichia coli P Pili under Stress Biophys. J., March 1, 2006; 90(5): 1521 - 1534. [Abstract] [Full Text] [PDF]

				C. Hyeon and D. Thirumalai Chemical Theory and Computation Special Feature: Mechanical unfolding of RNA hairpins PNAS, May 10, 2005; 102(19): 6789 - 6794. [Abstract] [Full Text] [PDF]

				M. Manosas and F. Ritort Thermodynamic and Kinetic Aspects of RNA Pulling Experiments Biophys. J., May 1, 2005; 88(5): 3224 - 3242. [Abstract] [Full Text] [PDF]

				F. Liu and Z.-c. Ou-Yang Monte Carlo Simulation for Single RNA Unfolding by Force Biophys. J., January 1, 2005; 88(1): 76 - 84. [Abstract] [Full Text] [PDF]

This Article Abstract Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Gerland, U. Articles by Hwa, T.