Mechanically probing the folding pathway of single RNA molecules

Self-Attractive Semiflexible Polymers under an External Force Field

The dynamical response of a tethered semiflexible polymer with self-attractive interactions and subjected to an external force field is numerically investigated by varying stiffness and self-interaction strength. The chain is confined in two spatial dimensions and placed in contact with a heat bath described by the Brownian multi-particle collision method. For strong self-attraction the equilibrium conformations range from compact structures to double-stranded chains, and to rods when increasing the stiffness. Under the external field at small rigidities, the initial close-packed chain is continuously unwound by the force before being completely elongated. For double-stranded conformations the transition from the folded state to the open one is sharp being steeper for larger stiffnesses. The discontinuity in the transition appears in the force-extension relation, as well as in the probability distribution function of the gyration radius. The relative deformation with respect to the equilibrium case along the direction normal to the force is found to decay as the inverse of the applied force.

Predicting Cotranscriptional Folding Kinetics For Riboswitch

On the basis of a helix-based transition rate model, we developed a new method for sampling cotranscriptional RNA conformational ensemble and the prediction of cotranscriptional folding kinetics. Applications to E. coli. SRP RNA and pbuE riboswitch indicate that the model may provide reliable predictions for the cotranscriptional folding pathways and population kinetics. For E. coli. SRP RNA, the predicted population kinetics and the folding pathway are consistent with the SHAPE profiles in the recent cotranscriptional SHAPE-seq experiments. For the pbuE riboswitch, the model predicts the transcriptional termination efficiency as a function of the force. The theoretical results show (a) a force-induced transition from the aptamer (antiterminator) to the terminator structure and (b) the different folding pathways for the riboswitch with and without the ligand (adenine). More specifically, without adenine, the aptamer structure emerges as a short-lived kinetic transient state instead of a thermodynamically stable intermediate state. Furthermore, from the predicted extension-time curves, the model identifies a series of conformational switches in the pulling process, where the predicted relative residence times for the different structures are in accordance with the experimental data. The model may provide a new tool for quantitative predictions of cotranscriptional folding kinetics, and results can offer useful insights into cotranscriptional folding-related RNA functions such as regulation of gene expression with riboswitches.

Statistical mechanics of the Huxley-Simmons model

The chemomechanical model of Huxley and Simmons (HS) [A. F. Huxley and R. M. Simmons, Nature 233, 533 (1971)NATUAS0028-083610.1038/233533a0] provides a paradigmatic description of mechanically induced collective conformational changes relevant in a variety of biological contexts, from muscles power stroke and hair cell gating to integrin binding and hairpin unzipping. We develop a statistical mechanical perspective on the HS model by exploiting a formal analogy with a paramagnetic Ising model. We first study the equilibrium HS model with a finite number of elements and compute explicitly its mechanical and thermal properties. To model kinetics, we derive a master equation and solve it for several loading protocols. The developed formalism is applicable to a broad range of allosteric systems with mean-field interactions.

Elasticity of a semiflexible filament with a discontinuous tension due to a cross-link or a molecular motor

We analyze the stretching elasticity of a wormlike chain with a tension discontinuity resulting from a Hookean spring connecting its backbone to a fixed point. The elasticity of isolated semiflexible filaments has been the subject in a significant body of literature, primarily because of its relevance to the mechanics of biological matter. In real systems, however, these filaments are usually part of supramolecular structures involving cross-linkers or molecular motors, which cause tension discontinuities. Our model is intended as a minimal structural element incorporating such a discontinuity. We obtain analytical results in the weakly bending limit of the filament, concerning its force-extension relation and the response of the two parts in which the filament is divided by the spring. For a small tension discontinuity, the linear response of the filament extension to this discontinuity strongly depends on the external tension. For large external tension f, the spring force contributes a subdominant correction ∼1/f^{3/2} to the well-known ∼1/sqrt[f]-dependence of the end-to-end extension.

Calculation of a fluctuating entropic force by phase space sampling

A polymer chain pinned in space exerts a fluctuating force on the pin point in thermal equilibrium. The average of such fluctuating force is well understood from statistical mechanics as an entropic force, but little is known about the underlying force distribution. Here, we introduce two phase space sampling methods that can produce the equilibrium distribution of instantaneous forces exerted by a terminally pinned polymer. In these methods, both the positions and momenta of mass points representing a freely jointed chain are perturbed in accordance with the spatial constraints and the Boltzmann distribution of total energy. The constraint force for each conformation and momentum is calculated using Lagrangian dynamics. Using terminally pinned chains in space and on a surface, we show that the force distribution is highly asymmetric with both tensile and compressive forces. Most importantly, the mean of the distribution, which is equal to the entropic force, is not the most probable force even for long chains. Our work provides insights into the mechanistic origin of entropic forces, and an efficient computational tool for unbiased sampling of the phase space of a constrained system.

Elastic properties and secondary structure formation of single-stranded DNA at monovalent and divalent salt conditions

Single-stranded DNA (ssDNA) plays a major role in several biological processes. It is therefore of fundamental interest to understand how the elastic response and the formation of secondary structures are modulated by the interplay between base pairing and electrostatic interactions. Here we measure force-extension curves (FECs) of ssDNA molecules in optical tweezers set up over two orders of magnitude of monovalent and divalent salt conditions, and obtain its elastic parameters by fitting the FECs to semiflexible models of polymers. For both monovalent and divalent salts, we find that the electrostatic contribution to the persistence length is proportional to the Debye screening length, varying as the inverse of the square root of cation concentration. The intrinsic persistence length is equal to 0.7 nm for both types of salts, and the effectivity of divalent cations in screening electrostatic interactions appears to be 100-fold as compared with monovalent salt, in line with what has been recently reported for single-stranded RNA. Finally, we propose an analysis of the FECs using a model that accounts for the effective thickness of the filament at low salt condition and a simple phenomenological description that quantifies the formation of non-specific secondary structure at low forces.

Quantifying transient 3D dynamical phenomena of single mRNA particles in live yeast cell measurements

Single-particle tracking (SPT) has been extensively used to obtain information about diffusion and directed motion in a wide range of biological applications. Recently, new methods have appeared for obtaining precise (10s of nm) spatial information in three dimensions (3D) with high temporal resolution (measurements obtained every 4 ms), which promise to more accurately sense the true dynamical behavior in the natural 3D cellular environment. Despite the quantitative 3D tracking information, the range of mathematical methods for extracting information about the underlying system has been limited mostly to mean-squared displacement analysis and other techniques not accounting for complex 3D kinetic interactions. There is a great need for new analysis tools aiming to more fully extract the biological information content from in vivo SPT measurements. High-resolution SPT experimental data has enormous potential to objectively scrutinize various proposed mechanistic schemes arising from theoretical biophysics and cell biology. At the same time, methods for rigorously checking the statistical consistency of both model assumptions and estimated parameters against observed experimental data (i.e., goodness-of-fit tests) have not received great attention. We demonstrate methods enabling (1) estimation of the parameters of 3D stochastic differential equation (SDE) models of the underlying dynamics given only one trajectory; and (2) construction of hypothesis tests checking the consistency of the fitted model with the observed trajectory so that extracted parameters are not overinterpreted (the tools are applicable to linear or nonlinear SDEs calibrated from nonstationary time series data). The approach is demonstrated on high-resolution 3D trajectories of single ARG3 mRNA particles in yeast cells in order to show the power of the methods in detecting signatures of transient directed transport. The methods presented are generally relevant to a wide variety of 2D and 3D SPT tracking applications.

Single molecules in an extension clamp: extracting rates and activation barriers

When a macromolecule, held at a fixed end-to-end separation, undergoes conformational rearrangements, the fluctuating mechanical force generated by the molecule can be used as a reporter of the molecular conformational dynamics. We present an analytical framework for extracting the intrinsic rates of conformational transitions and the locations and heights of the rate-limiting barriers from such extension clamp measurements. The unique nature of the bias imposed by the extension clamp on the activation barriers allows access to biomolecular transitions currently not accessible in pulling experiments. A mapping rule is established between the outputs of different types of experiments, providing information about poorly accessible regions of the molecular landscape.

Theory for RNA folding, stretching, and melting including loops and salt

Secondary structure formation of nucleic acids strongly depends on salt concentration and temperature. We develop a theory for RNA folding that correctly accounts for sequence effects, the entropic contributions associated with loop formation, and salt effects. Using an iterative expression for the partition function that neglects pseudoknots, we calculate folding free energies and minimum free energy configurations based on the experimentally derived basepairing free energies. The configurational entropy of loop formation is modeled by the asymptotic expression -clnm, where m is the length of the loop and c the loop exponent, which is an adjustable constant. Salt effects enter in two ways: first, we derive salt-induced modifications of the free energy parameters for describing basepairing, and second, we include the electrostatic free energy for loop formation. Both effects are modeled on the Debye-Hückel level including counterion condensation. We validate our theory for two different RNA sequences. For tRNA-phe, the resultant heat capacity curves for thermal denaturation at various salt concentrations accurately reproduce experimental results. For the P5ab RNA hairpin, we derive the global phase diagram in the three-dimensional space spanned by temperature, stretching force, and salt concentration and obtain good agreement with the experimentally determined critical unfolding force. We show that for a proper description of RNA melting and stretching, both salt and loop entropy effects are needed.

Temperature and chemical denaturant dependence of forced unfolding of titin I27

Single-molecule force measurement opens a new door for investigating detailed biomolecular interactions and their thermodynamic properties by pulling molecules apart while monitoring the force exerted on them. Recent advances in the nonequilibrium work theorem allows one to determine the free-energy landscapes of these events. Such information is valuable for understanding processes such as protein and RNA folding and receptor-ligand binding. Here, we used force as a physical parameter under the traditional chemical and temperature denaturing environment to alter the protein folding energy landscape and compared the change in the unfolding free-energy barrier of the I27 domain of human cardiac titin. We found that the trends in protein unfolding free-energy barriers are consistent for single-molecule force measurements and bulk chemical and temperature studies. The results suggest that the information from single-molecule pulling experiments are meaningful and useful for understanding the mechanism of folding of titin I27.

Equilibrium properties and force-driven unfolding pathways of RNA molecules

The mechanical unfolding of a simple RNA hairpin and of a 236-base portion of the Tetrahymena thermophila ribozyme is studied by means of an Ising-like model. Phase diagrams and free energy landscapes are computed exactly and suggest a simple two-state behavior for the hairpin and the presence of intermediate states for the ribozyme. Nonequilibrium simulations give the possible unfolding pathways for the ribozyme, and the dominant pathway corresponds to the experimentally observed one.

A new computational approach for mechanical folding kinetics of RNA hairpins

Based on an ensemble of kinetically accessible conformations, we propose a new analytical model for RNA folding kinetics. The model gives populational kinetics, kinetic rates, transition states, and pathways from the rate matrix. Applications of the new kinetic model to mechanical folding of RNA hairpins such as trans-activation-responsive RNA reveal distinct kinetic behaviors in different force regimes, from zero force to forces much stronger than the critical force for the folding-unfolding transition. In the absence of force or a low force, folding can be initiated (nucleated) at any position by forming the first base stack and there exist many pathways for the folding process. In contrast, for a higher force, the folding/unfolding would predominantly proceed along a single zipping/unzipping pathway. Studies for different hairpin-forming sequences indicate that depending on the nucleotide sequence, a kinetic intermediate can emerge in the low force regime but disappear in high force regime, and a new kinetic intermediate, which is absent in the low and high force regimes, can emerge in the medium force range. Variations of the force lead to changes in folding cooperativity and rate-limiting steps. The predicted network of pathways for trans-activation-responsive RNA suggests two parallel dominant pathways. The rate-limiting folding steps (at f = 8 pN) are the formation of specific basepairs that are 2-4 basepairs away from the loop. At a higher force (f = 11 pN), the folding rate is controlled by the formation of the bulge loop. The predicted rates and transition states are in good agreement with the experimental data for a broad force regime.

Quantifying DNA melting transitions using single-molecule force spectroscopy

We stretched a DNA molecule using atomic force microscope and quantified the mechanical properties associated with B and S forms of double-stranded DNA (dsDNA), molten DNA, and single-stranded DNA (ssDNA). We also fit overdamped diffusion models to the AFM time series and used these models to extract additional kinetic information about the system. Our analysis provides additional evidence supporting the view that S-DNA is a stable intermediate encountered during dsDNA melting by mechanical force. In addition, we demonstrated that the estimated diffusion models can detect dynamical signatures of conformational degrees of freedom not directly observed in experiments.

RNA folding: conformational statistics, folding kinetics, and ion electrostatics

RNA folding is a remarkably complex problem that involves ion-mediated electrostatic interaction, conformational entropy, base pairing and stacking, and noncanonical interactions. During the past decade, results from a variety of experimental and theoretical studies pointed to (a) the potential ion correlation effect in Mg2+-RNA interactions, (b) the rugged energy landscapes and multistate RNA folding kinetics even for small RNA systems such as hairpins and pseudoknots, (c) the intraloop interactions and sequence-dependent loop free energy, and (d) the strong nonadditivity of chain entropy in RNA pseudoknot and other tertiary folds. Several related issues, which have not been thoroughly resolved, require combined approaches with thermodynamic and kinetic experiments, statistical mechanical modeling, and all-atom computer simulations.

Predicting ribosomal frameshifting efficiency

Many retroviruses use -1 ribosomal frameshifting as part of the mechanism in translational control of viral protein synthesis. Quantitative prediction of the efficiency of -1 frameshifting is crucial for understanding the viral gene expression. Here we investigate the free energy landscape for a minimal -1 programmed ribosomal frameshifting machinery, including the codon-anticodon base pairs at the slippery site, the downstream messenger RNA structure and the spacer between the slippery site and the downstream structure. The free energy landscape analysis leads to a quantitative relationship between the frameshifting efficiency and the tension force generated during the movement of codon-anticodon complexes, which may occur in the A/T to A/A accommodation process or the translocation process. The analysis shows no consistent correlation between frameshifting efficiency and global stability of the downstream mRNA structure.

Compressing a rigid filament: buckling and cyclization

We study elastic properties of rigid filaments modeled as stiff chains shorter than their persistence length. By rigid filaments we mean that fluctuations around the optimal filament shape are weak and that low-order expansions (quadratic or quartic) in the deviation from the optimal shape are sufficient to describe them. Our main interest lies in the profiles of force vs. projected filament length, closure probability and weakly buckled states. Results may be relevant to experiments on self-assembled biological (microtubules, actin filaments) and synthetic (organo-gelators) filaments, carbon nanotubes and polymers grafted with strongly repelling side chains, some of which are discussed here.

Folding thermodynamics of pseudoknotted chain conformations

We develop a statistical mechanical framework for the folding thermodynamics of pseudoknotted structures. As applications of the theory, we investigate the folding stability and the free energy landscapes for both the thermal and the mechanical unfolding of pseudoknotted chains. For the mechanical unfolding process, we predict the force-extension curves, from which we can obtain the information about structural transitions in the unfolding process. In general, a pseudoknotted structure unfolds through multiple structural transitions. The interplay between the helix stems and the loops plays an important role in the folding stability of pseudoknots. For instance, variations in loop sizes can lead to the destabilization of some intermediate states and change the (equilibrium) folding pathways (e.g., two helix stems unfold either cooperatively or sequentially). In both thermal and mechanical unfolding, depending on the nucleotide sequence, misfolded intermediate states can emerge in the folding process. In addition, thermal and mechanical unfoldings often have different (equilibrium) pathways. For example, for certain sequences, the misfolded intermediates, which generally have longer tails, can fold, unfold, and refold again in the pulling process, which means that these intermediates can switch between two different average end-end extensions.

Pathways and kinetic barriers in mechanical unfolding and refolding of RNA and proteins

Using self-organized polymer models, we predict mechanical unfolding and refolding pathways of ribozymes, and the green fluorescent protein. In agreement with experiments, there are between six and eight unfolding transitions in the Tetrahymena ribozyme. Depending on the loading rate, the number of rips in the force-ramp unfolding of the Azoarcus ribozymes is between two and four. Force-quench refolding of the P4-P6 subdomain of the Tetrahymena ribozyme occurs through a compact intermediate. Subsequent formation of tertiary contacts between helices P5b-P6a and P5a/P5c-P4 leads to the native state. The force-quench refolding pathways agree with ensemble experiments. In the dominant unfolding route, the N-terminal alpha helix of GFP unravels first, followed by disruption of the N terminus beta strand. There is a third intermediate that involves disruption of three other strands. In accord with experiments, the force-quench refolding pathway of GFP is hierarchic, with the rate-limiting step being the closure of the barrel.

Mechanical unfolding of RNA: from hairpins to structures with internal multiloops

Mechanical unfolding of RNA structures, ranging from hairpins to ribozymes, using laser optical tweezer experiments have begun to reveal the features of the energy landscape that cannot be easily explored using conventional experiments. Upon application of constant force (f), RNA hairpins undergo cooperative transitions from folded to unfolded states whereas subdomains of ribozymes unravel one at a time. Here, we use a self-organized polymer model and Brownian dynamics simulations to probe mechanical unfolding at constant force and constant-loading rate of four RNA structures of varying complexity. For simple hairpins, such as P5GA, application of constant force or constant loading rate results in bistable cooperative transitions between folded and unfolded states without populating any intermediates. The transition state location (DeltaxFTS) changes dramatically as the loading rate is varied. At loading rates comparable to those used in laser optical tweezer experiments, the hairpin is plastic, with DeltaxFTS being midway between folded and unfolded states; whereas at high loading rates, DeltaxFTS moves close to the folded state, i.e., RNA is brittle. For the 29-nucleotide TAR RNA with the three-nucleotide bulge, unfolding occurs in a nearly two-state manner with an occasional pause in a high free energy metastable state. Forced unfolding of the 55 nucleotides of the Hepatitis IRES domain IIa, which has a distorted L-shaped structure, results in well-populated stable intermediates. The most stable force-stabilized intermediate represents straightening of the L-shaped structure. For these structures, the unfolding pathways can be predicted using the contact map of the native structures. Unfolding of a RNA motif with internal multiloop, namely, the 109-nucleotide prohead RNA that is part of the 29 DNA packaging motor, at constant value of rf occurs with three distinct rips that represent unraveling of the paired helices. The rips represent kinetic barriers to unfolding. Our work shows 1), the response of RNA to force is largely determined by the native structure; and 2), only by probing mechanical unfolding over a wide range of forces can the underlying energy landscape be fully explored.

The role of flexible tethers in multiple ligand-receptor bond formation between curved surfaces

Ligands mounted to surfaces via extensible tethers are present in nature and represent a growing class of molecules used to engineer adhesion in drug targeting, biosensing, self-assembling nanostructures, and in other biophysical research. Using a continuum approach with geometric and thermodynamic arguments, we derive a number of analytical expressions that relate key properties of single-tethered ligand-receptor interactions to multiple bond formation between curved surfaces. The theoretical predictions are in good agreement with measurements made with the surface forces apparatus. We establish that, when ligated, many tethers commonly used in biophysical research exhibit a discrete binding range that can be accurately measured with force spectroscopy. The distribution of bound ligated tethers is independent of the surfaces' interaction radius, R. The bridging force scales linearly with R, the tether's effective spring constant and grafting density, and with the ligand-receptor bond energy when the surfaces are in direct contact. These results are contrasted to bridging forces that evolve between plane-parallel geometries. Last, we show how our simple analytical reductions can be used to predict adhesive forces for STEALTH liposomes and other targeted and self-assembled nanoparticles.

Single-molecule experiments in biological physics: methods and applications

I review single-molecule experiments (SMEs) in biological physics. Recent technological developments have provided the tools to design and build scientific instruments of high enough sensitivity and precision to manipulate and visualize individual molecules and measure microscopic forces. Using SMEs it is possible to manipulate molecules one at a time and measure distributions describing molecular properties, characterize the kinetics of biomolecular reactions and detect molecular intermediates. SMEs provide additional information about thermodynamics and kinetics of biomolecular processes. This complements information obtained in traditional bulk assays. In SMEs it is also possible to measure small energies and detect large Brownian deviations in biomolecular reactions, thereby offering new methods and systems to scrutinize the basic foundations of statistical mechanics. This review is written at a very introductory level, emphasizing the importance of SMEs to scientists interested in knowing the common playground of ideas and the interdisciplinary topics accessible by these techniques. The review discusses SMEs from an experimental perspective, first exposing the most common experimental methodologies and later presenting various molecular systems where such techniques have been applied. I briefly discuss experimental techniques such as atomic-force microscopy (AFM), laser optical tweezers (LOTs), magnetic tweezers (MTs), biomembrane force probes (BFPs) and single-molecule fluorescence (SMF). I then present several applications of SME to the study of nucleic acids (DNA, RNA and DNA condensation) and proteins (protein-protein interactions, protein folding and molecular motors). Finally, I discuss applications of SMEs to the study of the nonequilibrium thermodynamics of small systems and the experimental verification of fluctuation theorems. I conclude with a discussion of open questions and future perspectives.

A puzzle in DNA biophysics

In this issue, Lee et al. report the experimental temperature-dependence of the unzipping force for two natural DNA sequences. For both sequences, the curves show an anomaly at temperatures around 40 degrees C. In this brief contribution, we stress that the anomaly is not easily explained within the established theoretical models for the biophysics of DNA. As this puzzle questions our basic understanding of DNA, it must be resolved, most likely by a combination of additional experiments and new theoretical work.

Sub-angstrom conformational changes of a single molecule captured by AFM variance analysis

A system's equilibrium variance can be analyzed to probe its underlying dynamics at higher resolution. Here, using single-molecule atomic-force microscope techniques, we show how the variance in the length of a single dextran molecule can be used to establish thermodynamic equilibrium and to detect conformational changes not directly observable with other methods. Dextran is comprised of a chain of pyranose rings that each undergoes an Angstrom-scale transition from a chair to boat conformation under a stretching force. Our analysis of the variance of the molecule's fluctuations verifies equilibrium throughout the force-extension curve, consistent with the expected thermodynamic ensemble. This validates further analysis of the variance in the transition region, which reveals an intermediate conformation between the chair and the boat on the sub-Angstrom scale. Our test of thermal equilibrium as well as our variance analysis can be readily extended to a wide variety of molecules, including proteins.

Forced-unfolding and force-quench refolding of RNA hairpins

Nanomanipulation of individual RNA molecules, using laser optical tweezers, has made it possible to infer the major features of their energy landscape. Time-dependent mechanical unfolding trajectories, measured at a constant stretching force (f(S)) of simple RNA structures (hairpins and three-helix junctions) sandwiched between RNA/DNA hybrid handles show that they unfold in a reversible all-or-none manner. To provide a molecular interpretation of the experiments we use a general coarse-grained off-lattice Gō-like model, in which each nucleotide is represented using three interaction sites. Using the coarse-grained model we have explored forced-unfolding of RNA hairpin as a function of f(S) and the loading rate (r(f)). The simulations and theoretical analysis have been done both with and without the handles that are explicitly modeled by semiflexible polymer chains. The mechanisms and timescales for denaturation by temperature jump and mechanical unfolding are vastly different. The directed perturbation of the native state by f(S) results in a sequential unfolding of the hairpin starting from their ends, whereas thermal denaturation occurs stochastically. From the dependence of the unfolding rates on r(f) and f(S) we show that the position of the unfolding transition state is not a constant but moves dramatically as either r(f) or f(S) is changed. The transition-state movements are interpreted by adopting the Hammond postulate for forced-unfolding. Forced-unfolding simulations of RNA, with handles attached to the two ends, show that the value of the unfolding force increases (especially at high pulling speeds) as the length of the handles increases. The pathways for refolding of RNA from stretched initial conformation, upon quenching f(S) to the quench force f(Q), are highly heterogeneous. The refolding times, upon force-quench, are at least an order-of-magnitude greater than those obtained by temperature-quench. The long f(Q)-dependent refolding times starting from fully stretched states are analyzed using a model that accounts for the microscopic steps in the rate-limiting step, which involves the trans to gauche transitions of the dihedral angles in the GAAA tetraloop. The simulations with explicit molecular model for the handles show that the dynamics of force-quench refolding is strongly dependent on the interplay of their contour length and persistence length and the RNA persistence length. Using the generality of our results, we also make a number of precise experimentally testable predictions.

Force unfolding single RNAs

We develop a continue time Monte Carlo algorithm to simulate single RNAs unfolded by a time-dependent external force on the secondary structure level. Two recent unfolding RNA experiments carried out by Bustamante group are mainly investigated. We find that, in contrast to popular two-state assumption about the RNAs free energy landscape along the molecular extension, the molecules used in the experiments do not present apparent energy barriers. The strong cooperative folding and unfolding transitions of the RNAs observed in the experiments and in our simulations arise from the interaction of the molecules and the light trap. In addition, we also investigate the properties of Jarzynski's remarkable equality, whose experimental test has received considerable attention.

A sticky chain model of the elongation and unfolding of Escherichia coli P pili under stress

A model of the elongation of P pili expressed by uropathogenic Escherichia coli exposed to stress is presented. The model is based upon the sticky chain concept, which is based upon Hooke's law for elongation of the layer-to-layer and head-to-tail bonds between neighboring units in the PapA rod and a kinetic description of the opening and closing of bonds, described by rate equations and an energy landscape model. It provides an accurate description of the elongation behavior of P pili under stress and supports a hypothesis that the PapA rod shows all three basic stereotypes of elongation/unfolding: elongation of bonds in parallel, the zipper mode of unfolding, and elongation and unfolding of bonds in series. The two first elongation regions are dominated by a cooperative bond opening, in which each bond is influenced by its neighbor, whereas the third region can be described by individual bond opening, in which the bonds open and close randomly. A methodology for a swift extraction of model parameters from force-versus-elongation measurements performed under equilibrium conditions is derived. Entities such as the free energy, the stiffness, the elastic elongation, the opening length of the various bonds, and the number of PapA units in the rod are determined.

Translocation of structured polynucleotides through nanopores

We investigate theoretically the translocation of structured RNA/DNA molecules through narrow pores which allow single but not double strands to pass. The unzipping of basepaired regions within the molecules presents significant kinetic barriers for the translocation process. We show that this circumstance may be exploited to determine the full basepairing pattern of polynucleotides, including RNA pseudoknots. The crucial requirement is that the translocation dynamics (i.e. the length of the translocated molecular segment) needs to be recorded as a function of time with a spatial resolution of a few nucleotides. This could be achieved, for instance, by applying a mechanical driving force for translocation and recording force-extension curves (FECs) with a device such as an atomic force microscope or optical tweezers. Our analysis suggests that, with this added spatial resolution, nanopores could be transformed into a powerful experimental tool to study the folding of nucleic acids.

Thermodynamic and kinetic aspects of RNA pulling experiments

Recent single-molecule pulling experiments have shown how it is possible to manipulate RNA molecules using laser tweezers. In this article we investigate a minimal model for the experimental setup which includes an RNA molecule connected to two polymers (handles) and a bead trapped in the optical potential and attached to one of the handles. We start by considering the case of small single-domain RNA molecules, which unfold in a cooperative way. The model qualitatively reproduces the experimental results and allows us to investigate the influence of the bead and handles on the unfolding reaction. A main ingredient of the model is to consider the appropriate statistical ensemble and the corresponding thermodynamic potential describing thermal fluctuations in the system. We then investigate several questions relevant to extract thermodynamic information from experimental data. The kinetics of unfolding is also studied by introducing a dynamical model. Finally, we apply the model to the more general problem of a multidomain RNA molecule with Mg(2+) tertiary contacts that unfolds in a sequential way.

Mechanical unfolding of RNA hairpins

Mechanical unfolding trajectories, generated by applying constant force in optical-tweezer experiments, show that RNA hairpins and the P5abc subdomain of the group I intron unfold reversibly. We use coarse-grained Go-like models for RNA hairpins to explore forced unfolding over a broad range of temperatures. A number of predictions that are amenable to experimental tests are made. At the critical force, the hairpin jumps between folded and unfolded conformations without populating any discernible intermediates. The phase diagram in the force-temperature (f, T) plane shows that the hairpin unfolds by an all-or-none process. The cooperativity of the unfolding transition increases dramatically at low temperatures. Free energy of stability, obtained from time averages of mechanical unfolding trajectories, coincides with ensemble averages, which establishes ergodicity. The hopping time between the native basin of attraction (NBA) and the unfolded basin increases dramatically along the phase boundary. Thermal unfolding is stochastic, whereas mechanical unfolding occurs in "quantized steps" with great variations in the step lengths. Refolding times, upon force quench, from stretched states to the NBA are at least an order of magnitude greater than folding times by temperature quench. Upon force quench from stretched states, the NBA is reached in at least three stages. In the initial stages, the mean end-to-end distance decreases nearly continuously, and there is a sudden transition to the NBA only in the last stage. Because of the generality of the results, we propose that similar behavior should be observed in force quench refolding of proteins.

Theoretical studies of the kinetics of mechanical unfolding of cross-linked polymer chains and their implications for single-molecule pulling experiments

We have used kinetic Monte Carlo simulations to study the kinetics of unfolding of cross-linked polymer chains under mechanical loading. As the ends of a chain are pulled apart, the force transmitted by each cross-link increases until it ruptures. The stochastic cross-link rupture process is assumed to be governed by first order kinetics with a rate that depends exponentially on the transmitted force. We have performed random searches to identify optimal cross-link configurations whose unfolding requires a large applied force (measure of strength) and/or large dissipated energy (measure of toughness). We found that such optimal chains always involve cross-links arranged to form parallel strands. The location of those optimal strands generally depends on the loading rate. Optimal chains with a small number of cross-links were found to be almost as strong and tough as optimal chains with a large number of cross-links. Furthermore, optimality of chains with a small number of cross-links can be easily destroyed by adding cross-links at random. The present findings are relevant for the interpretation of single molecule force probe spectroscopy studies of the mechanical unfolding of "load-bearing" proteins, whose native topology often involves parallel strand arrangements similar to the optimal configurations identified in the study.

Unfolding and unzipping of single-stranded DNA by stretching

We present a theoretical study of single-stranded DNA under stretching. Within the proposed framework, the effects of base pairing on the mechanical response of the molecule can be studied in combination with an arbitrary underlying model of chain elasticity. In a generic case, we show that the stretching curve of single-stranded DNA exhibits two distinct features: the second-order "unfolding" phase transition, and a sharp crossover, reminiscent of the first-order "unzipping" transition in double-stranded DNA. We apply the theory to the particular cases of wormlike chain and freely jointed chain models, and discuss the universal and model-dependent features of the mechanical response of single-stranded DNA. In particular, we show that variation of the width of the unzipping crossover with interaction strength is very sensitive to the energetics of hairpin loops. This opens another way of testing the elastic properties of ssDNA.

Monte Carlo simulation for single RNA unfolding by force

Using polymer elastic theory and known RNA free energies, we construct a Monte Carlo algorithm to simulate the single RNA folding and unfolding by mechanical force on the secondary structure level. For the constant force ensemble, we simulate the force-extension curves of the P5ab, P5abc deltaA, and P5abc molecules in equilibrium. For the constant extension ensemble, we focus on the mechanical behaviors of the RNA P5ab molecule, which include the unfolding force dependence on the pulling speed, the force-hysteresis phenomenon, and the coincidence of stretching-relaxing force-curves in thermal equilibrium. We particularly simulate the time traces of the end-to-end distance of the P5ab under the constant force in equilibrium, which also have been recorded in the recent experiment. The reaction rate constants for the folding and unfolding are calculated. Our results show that the agreement between the simulation and the experimental measurements is satisfactory.

Force as a useful variable in reactions: unfolding RNA

The effect of force on the thermodynamics and kinetics of reactions is described. The key parameters are the difference in end-to-end distance between reactant and product for thermodynamics, and the distance to the transition state for kinetics. I focus the review on experimental results on force unfolding of RNA. Methods to measure Gibbs free energies and kinetics for reversible and irreversible reactions are described. The use of the worm-like-chain model to calculate the effects of force on thermodynamics and kinetics is illustrated with simple models. The main purpose of the review is to describe the simple experiments that have been done so far, and to encourage more people to enter a field that is new and full of opportunities.

Stretching of homopolymers and contact order

Mechanical stretching of self-interacting homopolymers is studied through molecular dynamics simulations in which the polymers are pulled with constant speed. At temperatures below the compactification temperature, the force-extension curves show a plateau that corresponds to the situation in which the polymer adopts "ball-string" configurations. The dependence of rupture distances on contact order and the effects of temperature are similar to those found in the case of model proteins. The dependence of behavior on the pulling speed is logarithmic. In the entropic limit, above the compactification temperature, the rupturing of contacts shows a monotonic decrease of extension with the contact order. The attainment of this limit depends on the system size and the pulling speed.

RNA folding and unfolding

Single-molecule studies of RNA folding and unfolding are providing impressive details of the intermediates that occur and their rates of interconversion. The folding and unfolding of RNA are controlled by varying the concentration of magnesium ions and measuring fluorescence energy transfer, or by applying force to the RNA and measuring the end-to-end distance. The hierarchical nature of RNA folding - first secondary structure, then tertiary structure - makes the process susceptible to analysis and prediction.