Search

RNA duplex stability depends strongly on ionic conditions, and inside cells RNAs are exposed to both monovalent and multivalent ions. Despite recent advances, we do not have general methods to quantitatively account for the effects of monovalent and multivalent ions on RNA stability, and the thermodynamic parameters for secondary structure prediction have only been derived at 1M [Na(+)]. Here, by mechanically unfolding and folding a 20 bp RNA hairpin using optical tweezers, we study the RNA thermodynamics and kinetics at different monovalent and mixed monovalent/Mg(2+) salt conditions. We measure the unfolding and folding rupture forces and apply Kramers theory to extract accurate information about the hairpin free energy landscape under tension at a wide range of ionic conditions. We obtain non-specific corrections for the free energy of formation of the RNA hairpin and measure how the distance of the transition state to the folded state changes with force and ionic strength. We experimentally validate the Tightly Bound Ion model and obtain values for the persistence length of ssRNA. Finally, we test the approximate rule by which the non-specific binding affinity of divalent cations at a given concentration is equivalent to that of monovalent cations taken at 100-fold concentration for small molecular constructs.

We derive analytical expressions for the bias of the Jarzynski free-energy estimator from N nonequilibrium work measurements, for a generic work distribution. To achieve this, we map the estimator onto the random energy model in a suitable scaling limit parametrized by (logN)/μ, where μ measures the width of the lower tail of the work distribution, and then compute the finite-N corrections to this limit with different approaches for different regimes of (logN)/μ. We show that these expressions describe accurately the bias for a wide class of work distributions and exploit them to build an improved free-energy estimator from unidirectional work measurements. We apply the method to optical tweezers unfolding and refolding experiments on DNA hairpins of varying loop size and dissipation, displaying both near-Gaussian and non-Gaussian work distributions.

Accurate knowledge of the thermodynamic properties of nucleic acids is crucial to predicting their structure and stability. To date most measurements of base-pair free energies in DNA are obtained in thermal denaturation experiments, which depend on several assumptions. Here we report measurements of the DNA base-pair free energies based on a simplified system, the mechanical unzipping of single DNA molecules. By combining experimental data with a physical model and an optimization algorithm for analysis, we measure the 10 unique nearest-neighbor base-pair free energies with 0.1 kcal mol(-1) precision over two orders of magnitude of monovalent salt concentration. We find an improved set of standard energy values compared with Unified Oligonucleotide energies and a unique set of 10 base-pair-specific salt-correction values. The latter are found to be strongest for AA/TT and weakest for CC/GG. Our unique energy values and salt corrections improve predictions of DNA unzipping forces and are fully compatible with melting temperatures for oligos. The method should make it possible to obtain free energies, enthalpies, and entropies in conditions not accessible by bulk methodologies.

A main goal of single-molecule experiments is to evaluate equilibrium free energy differences by applying fluctuation relations to repeated work measurements along irreversible processes. We quantify the error that is made in a free energy estimate by means of the Jarzynski equality when the accumulated work expended on the whole system (including the instrument) is erroneously replaced by the work transferred to the subsystem consisting of the sole molecular construct. We find that the error may be as large as 100%, depending on the number of experiments and on the bandwidth of the data acquisition apparatus. Our theoretical estimate is validated by numerical simulations and pulling experiments on DNA hairpins using optical tweezers.

We present a method for determining the free energy of coexisting states from irreversible work measurements. Our approach is based on a fluctuation relation that is valid for dissipative transformations in partially equilibrated systems. To illustrate the validity and usefulness of the approach, we use optical tweezers to determine the free energy branches of the native and unfolded states of a two-state molecule as a function of the pulling control parameter. We determine, within 0.6kBT accuracy, the transition point where the free energies of the native and the unfolded states are equal.

The Jarzynski equality and the fluctuation theorem relate equilibrium free energy differences to nonequilibrium measurements of the work. These relations extend to single-molecule experiments that have probed the finite-time thermodynamics of proteins and nucleic acids. The effects of experimental error and instrument noise have not been considered previously. Here, we present a Bayesian formalism for estimating free energy changes from nonequilibrium work measurements that compensates for instrument noise and combines data from multiple driving protocols. We reanalyze a recent set of experiments in which a single RNA hairpin is unfolded and refolded using optical tweezers at three different rates. Interestingly, the fastest and farthest-from-equilibrium measurements contain the least instrumental noise and, therefore, provide a more accurate estimate of the free energies than a few slow, more noisy, near-equilibrium measurements. The methods we propose here will extend the scope of single-molecule experiments; they can be used in the analysis of data from measurements with atomic force microscopy, optical, and magnetic tweezers.

Recent single-molecule force measurements on single-domain proteins have highlighted a three-state folding mechanism where a stabilized intermediate state (I) is observed on the folding trajectory between the stretched state and the native state. Here we investigate on-lattice protein-like heteropolymer models that lead to a three-state mechanism and show that force experiments can be useful to determine the structure of I. We have mostly found that I is composed of a core stabilized by a high number of native contacts, plus an unstructured extended chain. The lifetime of I is shown to be sensitive to modifications of the protein that spoil the core. We then propose three types of modifications--point mutations, cuts, and circular permutations--aiming at: (1) confirming the presence of the core and (2) determining its location, within one amino acid accuracy, along the polypeptide chain. We also propose force jump protocols aiming to probe the on/off-pathway nature of I., (2007 Wiley-Liss, Inc.)

Experimental variables of optical tweezers instrumentation that affect RNA folding/unfolding kinetics were investigated. A model RNA hairpin, P5ab, was attached to two micron-sized beads through hybrid RNA/DNA handles; one bead was trapped by dual-beam lasers and the other was held by a micropipette. Several experimental variables were changed while measuring the unfolding/refolding kinetics, including handle lengths, trap stiffness, and modes of force applied to the molecule. In constant-force mode where the tension applied to the RNA was maintained through feedback control, the measured rate coefficients varied within 40% when the handle lengths were changed by 10-fold (1.1-10.2 Kbp); they increased by two- to threefold when the trap stiffness was lowered to one-third (from 0.1 to 0.035 pN/nm). In the passive mode, without feedback control and where the force applied to the RNA varied in response to the end-to-end distance change of the tether, the RNA hopped between a high-force folded-state and a low-force unfolded-state. In this mode, the rates increased up to twofold with longer handles or softer traps. Overall, the measured rates remained with the same order-of-magnitude over the wide range of conditions studied. In the companion article on pages 3010-3021, we analyze how the measured kinetics parameters differ from the intrinsic molecular rates of the RNA, and thus how to obtain the molecular rates.

Most natural processes occur far from equilibrium and cannot be treated within the framework of classical thermodynamics. In 1998, Oono and Paniconi [Oono, Y. & Paniconi, M. (1998) Prog. Theor. Phys. Suppl. 130, 29-44] proposed a general phenomenological framework, steady-state thermodynamics, encompassing nonequilibrium steady states and transitions between such states. In 2001, Hatano and Sasa [Hatano, T. & Sasa, S. (2001) Phys. Rev. Lett. 86, 3463-3466] derived a testable prediction of this theory. Specifically, they were able to show that the exponential average of Y, a quantity similar to a dissipated work, should be equal to zero for arbitrary transitions between nonequilibrium steady states, -ln = 0. We have tested this strong prediction by measuring the dissipation and fluctuations of microspheres optically driven through water. We have found that -ln approximately 0 for three different nonequilibrium systems, supporting Hatano and Sasa's proposed extension of thermodynamics to arbitrary steady states and irreversible transitions.

In 1997, Jarzynski proved a remarkable equality that allows one to compute the equilibrium free-energy difference DeltaF between two states from the probability distribution of the nonequilibrium work W done on the system to switch between the states, e-DeltaF/kappaT =e-W/kappaT [Jarzynski, C. (1997) Phys. Rev. Lett. 87, 2690-2693]. The Jarzynski equality provides a powerful free-energy difference estimator from a set of N irreversible experiments and is closely related to free-energy perturbation, a common computational technique for estimating free-energy differences. Despite the many applications of the Jarzynski estimator, its behavior is only poorly understood. In this article we derive the large N limit for the Jarzynski estimator bias, variance, and mean square error that is correct for arbitrary perturbations. We then analyze the properties of the Jarzynski estimator for all N when the probability distribution of work values is Gaussian, as occurs, for example, in the near-equilibrium regime. This allows us to quantitatively compare it to two other free-energy difference estimators: the mean work estimator and the fluctuation-dissipation theorem estimator. We show that, for near-equilibrium switching, the Jarzynski estimator is always superior to the mean work estimator and is even superior to the fluctuation-dissipation estimator for small N. The Jarzynski-estimator bias is shown to be the dominant source of error in many cases. Our expression for the bias is used to develop a bias-corrected Jarzynski free-energy difference estimator in the near-equilibrium regime.