Variance sum rule for entropy production.

Entropy production is the hallmark of nonequilibrium physics, quantifying irreversibility, dissipation, and the efficiency of energy transduction processes. Despite many efforts, its measurement at the nanoscale remains challenging. We introduce a variance sum rule (VSR) for displacement and force variances that permits us to measure the entropy production rate σ in nonequilibrium steady states. We first illustrate it for directly measurable forces, such as an active Brownian particle in an optical trap. We then apply the VSR to flickering experiments in human red blood cells. We find that σ is spatially heterogeneous with a finite correlation length, and its average value agrees with calorimetry measurements. The VSR paves the way to derive σ using force spectroscopy and time-resolved imaging in living and active matter.

Detection of single DNA mismatches by force spectroscopy in short DNA hairpins.

Identification of defective DNA structures is a difficult task, since small differences in base-pair bonding are hidden in the local structural variability of a generally random base-pair sequence. Defects, such as base mismatches, missing bases, crosslinks, and so on, occur in DNA with high frequency and must be efficiently identified and repaired to avoid dire consequences such as genetic mutations. Here, we focus on the detection of base mismatches, which is local deviations from the ideal Watson-Crick pairing rule, which may typically originate from DNA replication process, foreign chemical attack, or ionizing radiation. Experimental detection of a mismatch defect demands the ability to measure slight deviations in the free energy and molecular structure. We introduce different mismatches in short DNA hairpins (10 or 20 base pairs plus a 4-base loop) sandwiched between dsDNA handles to be used in single-molecule force spectroscopy with optical tweezers. We perform both hopping and force-pulling experiments to measure the excess free energies and deduce the characteristic kinetic signatures of the mismatch from the force-distance curves. All-atom molecular dynamics simulations lend support to the detailed interpretation of the experimental data. Such measurements, at the lowest sensitivity limits of this experimental technique, demonstrate the capability of identifying the presence of mismatches in a random complementary dsDNA sequence and provide lower bounds for the ability to distinguish different structural defects.

Experimental evidence of symmetry breaking of transition-path times.

While thermal rates of state transitions in classical systems have been studied for almost a century, associated transition-path times have only recently received attention. Uphill and downhill transition paths between states at different free energies should be statistically indistinguishable. Here, we systematically investigate transition-path-time symmetry and report evidence of its breakdown on the molecular- and meso-scale out of equilibrium. In automated Brownian dynamics experiments, we establish first-passage-time symmetries of colloids driven by femtoNewton forces in holographically-created optical landscapes confined within microchannels. Conversely, we show that transitions which couple in a path-dependent manner to fluctuating forces exhibit asymmetry. We reproduce this asymmetry in folding transitions of DNA-hairpins driven out of equilibrium and suggest a topological mechanism of symmetry breakdown. Our results are relevant to measurements that capture a single coordinate in a multidimensional free energy landscape, as encountered in electrophysiology and single-molecule fluorescence experiments.

Derivation of the spin-glass order parameter from stochastic thermodynamics.

A fluctuation relation is derived to extract the order parameter function q(x) in weakly ergodic systems. The relation is based on measuring and classifying entropy production fluctuations according to the value of the overlap q between configurations. For a fixed value of q, entropy production fluctuations are Gaussian distributed allowing us to derive the quasi-FDT so characteristic of aging systems. The theory is validated by extracting the q(x) in various types of glassy models. It might be generally applicable to other nonequilibrium systems and experimental small systems.

Force feedback effects on single molecule hopping and pulling experiments.

Single-molecule experiments with optical tweezers have become an important tool to study the properties and mechanisms of biological systems, such as cells and nucleic acids. In particular, force unzipping experiments have been used to extract the thermodynamics and kinetics of folding and unfolding reactions. In hopping experiments, a molecule executes transitions between the unfolded and folded states at a preset value of the force [constant force mode (CFM) under force feedback] or trap position [passive mode (PM) without feedback] and the force-dependent kinetic rates extracted from the lifetime of each state (CFM) and the rupture force distributions (PM) using the Bell-Evans model. However, hopping experiments in the CFM are known to overestimate molecular distances and folding free energies for fast transitions compared to the response time of the feedback. In contrast, kinetic rate measurements from pulling experiments have been mostly done in the PM while the CFM is seldom implemented in pulling protocols. Here, we carry out hopping and pulling experiments in a short DNA hairpin in the PM and CFM at three different temperatures (6 °C, 25 °C, and 45 °C) exhibiting largely varying kinetic rates. As expected, we find that equilibrium hopping experiments in the CFM and PM perform well at 6 °C (where kinetics are slow), whereas the CFM overestimates molecular parameters at 45 °C (where kinetics are fast). In contrast, nonequilibrium pulling experiments perform well in both modes at all temperatures. This demonstrates that the same kind of feedback algorithm in the CFM leads to more reliable determination of the folding reaction parameters in irreversible pulling experiments.

Control of force through feedback in small driven systems.

Controlling a time-dependent force applied to single molecules or colloidal particles is crucial for many types of experiments. Since in optical tweezers the primary controlled variable is the position of the trap, imposing a target force requires an active feedback process. We analyze this feedback process for the paradigmatic case of a nonequilibrium steady state generated by a dichotomous force protocol, first theoretically for a colloidal particle in a harmonic trap and then with both simulations and experiments for a long DNA hairpin. For the first setup, we find there is an optimal feedback gain separating monotonic from oscillatory response, whereas a too strong feedback leads to an instability. For the DNA molecule, reaching the target force requires substantial feedback gain since weak feedback cannot overcome the tendency to relax towards the equilibrium force.

Counter-propagating dual-trap optical tweezers based on linear momentum conservation.

We present a dual-trap optical tweezers setup which directly measures forces using linear momentum conservation. The setup uses a counter-propagating geometry, which allows momentum measurement on each beam separately. The experimental advantages of this setup include low drift due to all-optical manipulation, and a robust calibration (independent of the features of the trapped object or buffer medium) due to the force measurement method. Although this design does not attain the high-resolution of some co-propagating setups, we show that it can be used to perform different single molecule measurements: fluctuation-based molecular stiffness characterization at different forces and hopping experiments on molecular hairpins. Remarkably, in our setup it is possible to manipulate very short tethers (such as molecular hairpins with short handles) down to the limit where beads are almost in contact. The setup is used to illustrate a novel method for measuring the stiffness of optical traps and tethers on the basis of equilibrium force fluctuations, i.e., without the need of measuring the force vs molecular extension curve. This method is of general interest for dual trap optical tweezers setups and can be extended to setups which do not directly measure forces.

Fluctuation relation for weakly ergodic aging systems.

A fluctuation relation for aging systems is introduced and verified by extensive numerical simulations. It is based on the hypothesis of partial equilibration over phase-space regions in a scenario of entropy-driven relaxation. The relation provides a simple alternative method, amenable of experimental implementation, to measure replica symmetry breaking parameters in aging systems. The connection with the effective temperatures obtained from the fluctuation-dissipation theorem is discussed.

Force spectroscopy with dual-trap optical tweezers: molecular stiffness measurements and coupled fluctuations analysis.

Dual-trap optical tweezers are often used in high-resolution measurements in single-molecule biophysics. Such measurements can be hindered by the presence of extraneous noise sources, the most prominent of which is the coupling of fluctuations along different spatial directions, which may affect any optical tweezers setup. In this article, we analyze, both from the theoretical and the experimental points of view, the most common source for these couplings in dual-trap optical-tweezers setups: the misalignment of traps and tether. We give criteria to distinguish different kinds of misalignment, to estimate their quantitative relevance and to include them in the data analysis. The experimental data is obtained in a, to our knowledge, novel dual-trap optical-tweezers setup that directly measures forces. In the case in which misalignment is negligible, we provide a method to measure the stiffness of traps and tether based on variance analysis. This method can be seen as a calibration technique valid beyond the linear trap region. Our analysis is then employed to measure the persistence length of dsDNA tethers of three different lengths spanning two orders of magnitude. The effective persistence length of such tethers is shown to decrease with the contour length, in accordance with previous studies.

Non-specific binding of Na+ and Mg2+ to RNA determined by force spectroscopy methods.

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.

Improving signal/noise resolution in single-molecule experiments using molecular constructs with short handles.

We investigate unfolding/folding force kinetics in DNA hairpins exhibiting two and three states with newly designed short dsDNA handles (29 bp) using optical tweezers. We show how the higher stiffness of the molecular setup moderately enhances the signal/noise ratio (SNR) in hopping experiments as compared to conventional long-handled constructs (≅700 bp). The shorter construct results in a signal of higher SNR and slower folding/unfolding kinetics, thereby facilitating the detection of otherwise fast structural transitions. A novel analysis, as far as we are aware, of the elastic properties of the molecular setup, based on high-bandwidth measurements of force fluctuations along the folded branch, reveals that the highest SNR that can be achieved with short handles is potentially limited by the marked reduction of the effective persistence length and stretch modulus of the short linker complex.

Negative fluctuation-dissipation ratios in the backgammon model.

We analyze fluctuation-dissipation relations in the backgammon model: a system that displays glassy behavior at zero temperature due to the existence of entropy barriers. We study local and global fluctuation relations for the different observables in the model. For the case of a global perturbation we find a unique negative fluctuation-dissipation ratio that is independent of the observable and which diverges linearly with the waiting time. This result suggests that a negative effective temperature can be observed in glassy systems even in the absence of thermally activated processes.

Thermal activation and ATP dependence of the cytoskeleton remodeling dynamics.

The cytoskeleton (CSK) is a nonequilibrium polymer network that uses hydrolyzable sources of free energy such as adenosine triphosphate (ATP) to remodel its internal structure. As in inert nonequilibrium soft materials, CSK remodeling has been associated with structural rearrangements driven by energy-activated processes. We carry out particle tracking and traction microscopy measurements of alveolar epithelial cells at various temperatures and ATP concentrations. We provide the first experimental evidence that the remodeling dynamics of the CSK is driven by structural rearrangements over free-energy barriers induced by thermally activated forces mediated by ATP. The measured activation energy of these forces is approximately 40k_{B}T_{r} ( k_{B} being the Boltzmann constant and T_{r} being the room temperature). Our experiments provide clues to understand the analogy between the dynamics of the living CSK and that of inert nonequilibrium soft materials.

Statistical properties of metastable intermediates in DNA unzipping.

We unzip DNA molecules using optical tweezers and determine the sizes of the cooperatively unzipping and zipping regions separating consecutive metastable intermediates along the unzipping pathway. Sizes are found to be distributed following a power law, ranging from one base pair up to more than a hundred base pairs. We find that a large fraction of unzipping regions smaller than 10 bp are seldom detected because of the high compliance of the released single stranded DNA. We show how the compliance of a single nucleotide sets a limit value around 0.1 N/m for the stiffness of any local force probe aiming to discriminate one base pair at a time in DNA unzipping experiments.

Force-induced misfolding in RNA.

RNA folding is a kinetic process governed by the competition of a large number of structures stabilized by the transient formation of base pairs that may induce complex folding pathways and the formation of misfolded structures. Despite its importance in modern biophysics, the current understanding of RNA folding kinetics is limited by the complex interplay between the weak base pair interactions that stabilize the native structure and the disordering effect of thermal forces. The possibility of mechanically pulling individual molecules offers a new perspective to understand the folding of nucleic acids. Here we investigate the folding and misfolding mechanism in RNA secondary structures pulled by mechanical forces. We introduce a model based on the identification of the minimal set of structures that reproduce the patterns of force-extension curves obtained in single molecule experiments. The model requires only two fitting parameters: the attempt frequency at the level of individual base pairs and a parameter associated to a free-energy correction that accounts for the configurational entropy of an exponentially large number of neglected secondary structures. We apply the model to interpret results recently obtained in pulling experiments in the three-helix junction S15 RNA molecule (RNAS15). We show that RNAS15 undergoes force-induced misfolding where force favors the formation of a stable non-native hairpin. The model reproduces the pattern of unfolding and refolding force-extension curves, the distribution of breakage forces, and the misfolding probability obtained in the experiments.

Force unfolding kinetics of RNA using optical tweezers. II. Modeling experiments.

By exerting mechanical force, it is possible to unfold/refold RNA molecules one at a time. In a small range of forces, an RNA molecule can hop between the folded and the unfolded state with force-dependent kinetic rates. Here, we introduce a mesoscopic model to analyze the hopping kinetics of RNA hairpins in an optical tweezers setup. The model includes different elements of the experimental setup (beads, handles, and RNA sequence) and limitations of the instrument (time lag of the force-feedback mechanism and finite bandwidth of data acquisition). We investigated the influence of the instrument on the measured hopping rates. Results from the model are in good agreement with the experiments reported in the companion article. The comparison between theory and experiments allowed us to infer the values of the intrinsic molecular rates of the RNA hairpin alone and to search for the optimal experimental conditions to do the measurements. We conclude that the longest handles and softest traps that allow detection of the folding/unfolding signal (handles approximately 5-10 Kbp and traps approximately 0.03 pN/nm) represent the best conditions to obtain the intrinsic molecular rates. The methodology and rationale presented here can be applied to other experimental setups and other molecules.

Force-dependent fragility in RNA hairpins.

We apply Kramers theory to investigate the dissociation of multiple bonds under mechanical force and interpret experimental results for the unfolding and refolding force distributions of an RNA hairpin pulled at different loading rates using laser tweezers. We identify two different kinetic regimes depending on the range of forces explored during the unfolding and refolding process. The present approach extends the range of validity of the two-states approximation by providing a theoretical framework to reconstruct free-energy landscapes and identify force-induced structural changes in molecular transition states using single molecule pulling experiments. The method should be applicable to RNA hairpins with multiple kinetic barriers.

Condensation transition in DNA-polyaminoamide dendrimer fibers studied using optical tweezers.

When mixed together, DNA and polyaminoamide dendrimers form fibers that condense into a compact structure. We use optical tweezers to pull condensed fibers and investigate the decondensation transition by measuring force-extension curves (FECs). A characteristic force plateau (around 10 pN) and hysteresis between the pulling and relaxation cycles are observed for different dendrimer sizes, indicating the existence of a first-order transition between two phases (condensed and extended) of the fiber. Upon salt variation FECs change noticeably confirming that electrostatic forces drive the condensation transition. We propose a simple model for the decondensing transition that qualitatively reproduces the FECs and which is confirmed by atomic force microscopy images.

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.

Universality of fluctuation-dissipation ratios: the ferromagnetic model.

We calculate analytically the fluctuation-dissipation ratio (FDR) for Ising ferromagnets quenched to criticality, both for the long-range model and its short-range analog in the limit of large dimension. Our exact solution shows that, for both models, if the system is unmagnetized while if the initial magnetization is nonzero. This indicates that two different classes of critical coarsening dynamics need to be distinguished depending on the initial conditions, each with its own nontrivial FDR. We also analyze the dependence of the FDR on whether local and global observables are used. These results clarify how a proper local FDR (and the corresponding effective temperature) should be defined in long-range models in order to avoid spurious inconsistencies and maintain the expected correspondence between local and global results; global observables turn out to be far more robust tools for detecting nonequilibrium FDRs.