## New Research In

### Physical Sciences

### Social Sciences

#### Featured Portals

#### Articles by Topic

### Biological Sciences

#### Featured Portals

#### Articles by Topic

- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology

# Free-energy inference from partial work measurements in small systems

Edited* by Attila Szabo, National Institutes of Health, Bethesda, MD, and approved June 23, 2014 (received for review October 25, 2013)

## Significance

Fluctuation relations (FRs) provide general results about work (or total entropy production) distributions in nonequilibrium systems. However, in many cases the full work is not measurable and only partial work measurements are possible. The latter do not fulfill a FR and cannot be used to extract free-energy differences from irreversible work measurements. We show how FRs can be used to infer the full work distribution from partial work measurements. We illustrate the inference process, using dual-trap optical tweezers where two forces (one per trap) are measured, and we can access both the full and the partial work distributions. We derive a set of results of direct interest to single-molecule scientists and, more generally, to physicists and biophysicists.

## Abstract

Fluctuation relations (FRs) are among the few existing general results in nonequilibrium systems. Their verification requires the measurement of the total work performed on a system. Nevertheless in many cases only a partial measurement of the work is possible. Here we consider FRs in dual-trap optical tweezers where two different forces (one per trap) are measured. With this setup we perform pulling experiments on single molecules by moving one trap relative to the other. We demonstrate that work should be measured using the force exerted by the trap that is moved. The force that is measured in the trap at rest fails to provide the full dissipation in the system, leading to a (incorrect) work definition that does not satisfy the FR. The implications to single-molecule experiments and free-energy measurements are discussed. In the case of symmetric setups a second work definition, based on differential force measurements, is introduced. This definition is best suited to measure free energies as it shows faster convergence of estimators. We discuss measurements using the (incorrect) work definition as an example of partial work measurement. We show how to infer the full work distribution from the partial one via the FR. The inference process does also yield quantitative information, e.g., the hydrodynamic drag on the dumbbell. Results are also obtained for asymmetric dual-trap setups. We suggest that this kind of inference could represent a previously unidentified and general application of FRs to extract information about irreversible processes in small systems.

Fluctuation relations (FRs) are mathematical equations connecting non equilibrium work measurements to equilibrium free-energy differences. FRs, such as the Jarzynski equality (JE) or the Crooks fluctuation relation (CFR), have become a valuable tool in single-molecule biophysics where they are used to measure folding free energies from irreversible pulling experiments (1, 2). Such measurements have been carried out with laser optical tweezers on different nucleic acid structures such as hairpins (3⇓⇓–6), G quadruplexes (7, 8), and proteins (9⇓⇓–12) and with atomic force microscopes on proteins (13) and bimolecular complexes (14). An important issue regarding FRs is the correct definition of work, which rests on the correct identification of configurational variables and control parameters. In the single-trap optical tweezers configuration this issue has been thoroughly discussed (17⇓–19).

The situation of how to correctly measure work in small systems becomes subtle when there are different forces applied to the system. In this case theory gives the prescription to correctly define the work (*W*_{Γ}) for a given trajectory (Γ): Integrate the generalized force (*f*_{λ}) (conjugated to the control parameter, *λ*) over *λ* along Γ, *A*). In this case a dumbbell formed by a molecule tethered between two optically trapped beads is manipulated by moving one trap relative to the other. In this setup two different forces (one per trap) can be measured and at least two different work definitions are possible. In equilibrium conditions, i.e., when the traps are not moved, both forces are equivalent: The forces acting on each bead have equal magnitude and opposite sign. On the contrary, in pulling experiments, where one trap is at rest (with respect to water) while the other is moved, the two forces become inequivalent. This is so because the center of mass of the dumbbell drifts and the beads are affected by different viscous drags (purple arrows in Fig. 1*B*). In such conditions theory prescribes that the full thermodynamic work must be defined on the force measured at the moving trap whereas the force measured in the trap at rest (with respect to water) leads to a partial work measurement that, as we show below, entails a systematic error in free-energy estimates. The difference between both works equals the dissipation by the center of mass of the dumbbell, which is correctly accounted for only in the correct work definition.

In this paper we combine theory and experiments in a dual-trap setup to demonstrate several results. First, we show that if the wrong work definition is used, free-energy estimates will be flawed. The error is especially severe in the case of unidirectional work estimates (e.g., JE) whereas it influences bidirectional estimates to a lesser extent. This fact is not purely academic: Measuring the force in the trap at rest is experimentally easier in dual-trap setups and, in fact, many groups choose to do so (5, 20, 21). For example, if different lasers are used for trapping and detection, measuring the force in the moving trap poses the additional challenge of keeping the trapping and detection lasers aligned while moving them. Second, we show how it is possible, by using the CFR, to infer the full work distribution from partial work measurements. We elucidate the inference process in our dual-trap setup by showing how to reconstruct the correct work distribution (i.e., the one we would have measured in the moving optical trap) from partial work measurements in the wrong optical trap (i.e., the one at rest with respect to water) by using the CFR. In particular, for symmetric setups the correct work distribution can be directly inferred by simply shifting the partial work distribution. In asymmetric setups inference is still possible in the framework of a Gaussian approximation, but the knowledge of some equilibrium properties of the system is still required. This type of inference should be seen as an example of a more general application of FRs that aims at extracting information about the total entropy production of a nonequilibrium system from partial entropy production measurements. This allows us to determine the average power dissipated by the center of mass of the dumbbell, from which we can extract the corresponding hydrodynamic coefficient, thereby avoiding direct hydrodynamic measurements. Moreover, we argue this type of inference might find applicability to future biophysical experiments where the sources of entropy production are not directly measurable, e.g., ATP-dependent motor translocation where the hydrolysis reaction cycle cannot be followed one ATP at a time. Finally we show how, in symmetric setups, the work definition satisfying the CFR is not unique. In particular the differential work, based on differential force measurements (22), still satisfies the CFR and leads to the least-biased free-energy estimates. The distinguishing feature of this work definition is that it completely filters out the dissipation due to the motion of the center of mass of the dumbbell. Throughout this paper, and for simplicity and pedagogical reasons, most of the derivations, calculations, and experiments are shown for symmetric dual-trap setups whereas the asymmetric case is discussed toward the end of the paper.

## The Model

In dual-trap setups a molecule is stretched by two optical traps, the control parameter *λ* being the trap-to-trap distance. Let *A* and *B* denote the two optical traps. When one trap (say trap *A*) is moved with respect to the bath while trap *B* is at rest, suitable configurational variables are the positions of the beads, both measured from the center of the trap at rest (trap *B*). We denote these variables by *y*_{A} and *y*_{B}, Fig. 1*B*. The total energy of the system is composed of three terms,*U*_{m} describes the properties of the tether. However, one could measure the positions of the beads and the trap-to-trap distance in the moving frame of trap *A* (*x*_{A}, *x*_{B}, and −*λ* in Fig. 1*B*), and the potential energy in Eq. **1** would be written as*y*_{A} − *x*_{A} = *y*_{B} − *x*_{B} = *λ* (*x*_{A} and *x*_{B} are negative). Central to our analysis is the equation connecting the potential *U* and the work performed on the system by changing the control parameter (1):**1** we get*f*_{A} = −*k*_{A}(*y*_{A} − *λ*). Inserting *U*′ instead of *U* in [3] gives*W* and *W*′ are remarkably different. In fact, from the reference frame of trap *A*, the bath is seen to flow with velocity *B*). Because of this flow an experiment in which trap *A* is moved is not Galilean equivalent to one in which trap *B* is moved. In the presence of a flow the connection between potential and work, Eq. **3**, is not valid anymore. In fact, thermodynamic work measurements must be based on the force measured in the trap being moved. This fact has been discussed in ref. 23, yet the implications to single-molecule experiments have never been pointed out. Summarizing, if trap *A* is moved and *B* is at rest with respect to water, using the work *W* (Eq. **4**) in the JE leads to correct free-energy estimates whereas using *W*′ (Eq. **5**) in the JE leads to a systematic error. Below we quantify such error in detail. The difference between *W* and *W*′ can be readily discussed in symmetric setups (*k*_{A} = *k*_{B} = *k*) where calculations are much simpler. To do this we switch to a new coordinate system: *x*_{+} is the position of the geometric center of the dumbbell, whereas *x*_{−} is the differential coordinate (22). In this coordinate system the potential Eq. **1** reads *U*(*x*_{−}, *x*_{+}; *λ*) = *U*^{−}(*x*_{−}; *λ*) + *U*^{+}(*x*_{+}; *λ*):*U*_{+} associated to *x*_{+} is that of a moving trap, a problem that has been addressed both with experiments and with theory (24, 25), whereas the potential energy term *U*_{−} associated to *x*_{−} corresponds to pulling experiments performed using a single trap and a fixed point. The dumbbell is in contact with an isothermal bath where the equilibrium state is described by the Boltzmann distribution and the corresponding partition function (we assume a weak system–environment coupling, a situation satisfied in our experimental conditions, *SI Text*, sections S1 and S2). Consequently, the 2 df are uncoupled and the total partition function for the system factorizes,*β* = (*K*_{B}*T*)^{−1}, *T* being the temperature and *K*_{B} being Boltzmann constant. As a consequence free-energy changes in the system can be decomposed into two contributions:*W*^{−} contains the work done in stretching the molecule whereas *W*^{+} is pure dissipation due to the movement of the center of mass of the dumbbell. Note that*W* and *W*′ is entirely due to *W*^{+}. The JE holds for *W*, the standard work definition, so that**8** and **9** into Eq. **10**, we get Δ*G*^{+} + Δ*G*^{−} = −*β*^{−1}〈exp(−*β*(*W*^{+} + *W*^{−}))〉. In symmetric setups *W*^{+} and *W*^{−} are independently distributed random variables (*SI Text*, sections S1 and S3) and we can conclude that*W* and *W*′, we get two different free-energy estimates: Δ*G* (Eq. **10**) and Δ*G*′ = −*β*^{−1}log〈exp(−*βW*′)〉. The error *W*′ instead of *W* can be quantified as**. 9** and **10** and again using the fact *W*^{+} and *W*^{−} are independently distributed random variables, we get 〈exp(−*βW*)〉 = 〈exp(−*β*(*W*^{+} + *W*^{−}))〉 = 〈exp(−*βW*^{+})〉〈exp(−*βW*^{−})〉 and similarly for *W*′. As a consequence, **12** is reduced to*x*_{+} is subject to a quadratic potential (Eq. **6**), we expect *W*^{+} to be a Gaussian random variable. This is not true in general for *W*^{−} given that *x*_{−} is affected by the nonlinear term *U*_{m}. For Gaussian random variables we have*W*^{+}. Moreover, dragging a trapped bead in a fluid causes no free-energy change, so that**14** and **15** into Eq. **13**, we get**16** gives *G*′ is lower than Δ*G* (Eq. **12**). Interestingly enough, using *W*′ instead of *W* in the JE leads to free-energy estimates in apparent violation of the second law. The error on free-energy estimates obtained using *W*′ instead of *W* is proportional to the mean work performed on the center of the dumbbell. This mean work 〈*W*^{+}〉 is just the mean friction force times the total trap displacement Δ*λ*,*γ*_{+} is the friction coefficient of the drag force opposing the movement of the geometric center of the dumbbell. The value of *γ*_{+} can be independently obtained from equilibrium measurements (26) (*SI Text*, sections S4 and S5).

## Differential Work Measurements

Eqs. **9** and **15** show that free-energy estimates based on the standard work *W* and the differential work *W*^{−} are equivalent:*W*^{+} and *W*^{−} are independent random variables. The case of asymmetric setups is discussed further below. Therefore, *W*^{−} can be used for free-energy determination, as has been done in ref. 11, although without discussion. Eq. **18** does hold only when the number of work measurements, *N*, tends to infinity. In all practical cases we deal with finite *N* and the Jarzynski estimator is biased (27, 28). The bias is strongly linked to the typical dissipation *D*_{typ} and a reliable estimate of free-energy differences requires a number of work measurements, which scales as *N* ≃ exp(*D*_{typ}) (29, 30), so that even a small reduction in *D*_{typ} entails a considerable improvement in the convergence of free-energy estimators. Moreover, the bias is superadditive. Let us consider for simplicity Gaussian work distributions. In this case the bias, *N* limit is a function of the variance of the distribution *σ*^{2} and of *N* (27):*σ* and is superadditive; i.e., **19** was derived under strong assumptions, superadditivity also holds for other theoretical expressions for the bias and has been checked in our experimental data (see below). Let us introduce the following Jarzynski estimators for finite *N*,*W* = *W*^{+} + *W*^{−}, superadditivity guarantees**24** differential work measurements always improve the convergence of free-energy estimates in dual-trap setups. This is especially important in all those cases in which bidirectional methods (e.g., the CFR) cannot be used and one has to use unidirectional methods.

## Pulling on dsDNA

The theory discussed so far has been put to test in a series of pulling experiments performed in a recently developed dual-trap optical tweezers setup that directly measures force in each trap (31, 32). The setup can move the two optical traps independently and measure their relative position with subnanometer accuracy, giving direct access to both *W* and *W*′. In these experiments 3-kb dsDNA tethers (≃1 μm in contour length) were stretched between 1 pN and 3 pN (Fig. 1*C*) in a symmetric dual-trap setup (*k*_{A} = *k*_{B} = 0.02 pN/nm), using 4-μm silica beads as force probes. The experiments were performed moving one of the two traps (trap *A*) with respect to the laboratory frame and leaving trap *B* at rest. All experiments were performed in PBS buffer (pH 7.4, 1 M NaCl). We chose cyclical protocols (*T*_{c} is the total duration of the cyclic protocol). The excursion of the control parameter, *P*_{F}) and time-reversed (*P*_{R}) protocols:*P*_{F} = *P*_{R} = *P* and Δ*G* = 0 so that the CFR takes the form*D* can be directly tested in cases where negative dissipation events are observed. In Fig. 2*A* we show measured work histograms (solid points, left-hand side of Eq. **27**) and reconstructed histograms (open points, right-hand side of Eq. **27**). If Eq. **27** is fulfilled, then the measured and reconstructed histograms match each other. A quantitative measure of the deviation from Eq. **27** can be obtained from the ratio *P*(*D*)/*P*(−*D*), as shown in Fig. 2*B*. Experimental data show that *D* fulfills the FR whereas *D*′ does not. In Fig. 2*C* we show that the probability density functions of *D*^{+} and *D*^{−}, with**27**. *D*^{−} is just the differential work, *W*^{−}, Eq. **9** evaluated on a cyclic protocol, whereas *D*^{+} is the dissipation due to the movement of the center of mass of the dumbbell. Summarizing, although in general *W* is the only observable we expect to fulfill a FR, in symmetric setups two new FRs emerge, for *W*^{−} and *W*^{+}. In Fig. 3 *A* and *B* we compare the predictions of Eqs. **12** and **17** with experimental results for different pulling speeds, *λ*. Eq. **17** must be used to correct free-energy estimates obtained in all those dual-trap setups that do not measure the force applied by the trap that is being moved (as in refs. 5, 20, and 21). The advantages of using *W*^{−} in free-energy estimates are shown in Fig. 3*C*. There we show the convergence of the Jarzynski estimator with sample size for the cycles in Fig. 1*C*. Being evaluated over cycles, the expected free-energy change is zero. The convergence of the estimator is faster for *W*^{−} than for *W* in the three cases (Eq. **24**). The effect is enhanced in our experiments by the high pulling speed (in the range 1–7 μm) and by the large bead radius (2 μm). Let us note that due to the finite lifetime of molecular tethers and unavoidable drift effects, raising the pulling speed is a convenient strategy to improve the quality of free-energy estimates. Similar results have been found also at low pulling speeds where, again, *W*′ does not satisfy the CFR (*SI Text*, section S7).

## Experiments on DNA Hairpins

Fluctuation theorems are used to extract folding free energies for nucleic acid secondary structures or proteins. We further tested the different work definitions by performing pulling experiments on a 20-bp DNA hairpin (Fig. 4*A*) at a 0.96 ± 0.02-μm/s pulling speed in the same dual-trap setup as in the previous dsDNA experiments. In this case the work performed during the unfolding and that performed during the refolding of the molecule were considered separately, as is customary for free-energy determination. In Fig. 4*B* we present forward and reverse work histograms for *W*, *W*′, and *W*^{+}. Again *W* and *W*^{−} both fulfill the CFR but *W* shows higher dissipation than *W*^{−}, resulting in slower convergence of unidirectional free-energy estimators (Fig. 4*C*). The difference between unidirectional free-energy estimates based on *W* and *W*^{−} is in this case ≃1*K*_{B}*T*. As previously discussed for double-stranded DNA (Fig. 2*A*), *W*′ does not fulfill the CFR and, as a consequence, unidirectional free-energy estimates based on *W*′ are flawed. In our experimental conditions the error committed by using the wrong work definition is again positive and equal to *K*_{B}*T*. As previously discussed this leads to a negative average dissipated work, apparently violating the second law. It must be noted that the difference in free-energy estimates based on *W*^{−} and *W* is a finite-size effect, whose magnitude decreases when an increasing number of work measurements are considered. On the contrary the error committed by using *W*′ does not vanish by increasing the number of work measurements. It can be noted from Fig. 4*B* that, although *W*′ does not fulfill the CFR and gives wrong unidirectional estimates, its forward and reverse distributions apparently cross at *W*′ = Δ*G* within the experimental error. Although this could be used for free-energy determination, the result should be taken with caution as we have no general proof that this should happen in all cases.

## Free-Energy Inference from Partial Work Measurements

The JE and CFR are statements on the statistics of the total dissipation in irreversible thermodynamic transformations. The need to measure the total dissipation limits the range of applicability of these and other FRs. For example, testing FRs concerning the dynamics of molecular motors would need the simultaneous measurement of both the work performed by the motor and the number of hydrolyzed ATPs. Here we demonstrate that, at least in some cases, a different approach is possible. Let us start by considering the simple case of symmetric setups as developed in the previous sections. In the experiments discussed so far there are two sources of dissipation that we were able to measure and characterize separately: the motion of the dumbbell and the dissipation of the differential coordinate. We already learned that *W* satisfies a FR whereas *W*′ does not. Moreover, we know that *W* = (*W*_{−} + *W*_{+})/2 and *W*′ = (*W*_{−} − *W*_{+})/2 and that, being *W*^{+} and *W*^{−} uncorrelated random variables, *W* and *W*′ have the same variance. Imagine now to have only partial information on the system. For example, one could be able to measure force only in the trap at rest, as many experimental setups do. With this information, and in absence of any guiding principle, no statement about the total entropy production is possible. We will find such a guiding principle if we assume the FR to hold for *W*. Knowing that *W* and *W*′ have the same variance, we just have to shift the work distribution *P*′(*W*′) by *W* = *W*′ + Δ to get a new distribution that satisfies the CFR. In practice this is done starting from the set of *W*′ values and tuning the value of Δ (Fig. 5*A*) until *P*(*W*) = *P*′(*W* − Δ) fulfills the CFR (Fig. 5*B*). In the case of the hairpin, this same shifting procedure is operated for both forward and reverse work distributions. Again the value of Δ is tuned (Fig. 5*C*) until the CFR symmetry is recovered (Fig. 5*D*). The unique value of Δ that restores the validity of the CFR equals the average work dissipated by the motion of center of mass of the dumbbell, giving the hydrodynamic coefficient *γ*_{+} via Eq. **17**. Let us note that, once the work distribution in the correct trap (i.e., the moving trap) has been recovered, then we could also extract the correct free-energy difference (the value of Δ*G*, Fig. 5*E*) and infer the distribution for the differential work *W*_{−} by deconvolution.

The extension of this analysis to the asymmetric case is more complex but equally interesting (Fig. 6*A*). The decomposition of *W*, *W*′ in *W*^{+} and *W*^{−} (Eq. **9**) is still possible although *W*^{+}, *W*^{−} are not uncorrelated variables anymore and neither *W*^{+} nor *W*^{−} satisfies a FR. In this general case only *W* satisfies a FR (but not *W*′, *W*^{+}, or *W*^{−}). Remarkably enough, in the framework of a Gaussian approximation, it is still possible to infer the correct work distribution *P*(*W*) out of partial work *W*′ measurements. The analysis is presented in *SI Text*, section S6. In this case it is enough to know the trap and molecular stiffnesses *k*_{A}, *k*_{B}, *k*_{m}, i.e., some equilibrium properties of the system, for a successful inference. To reconstruct *P*(*W*) both the mean and the variance of *P*′(*W*′) must be changed, which can be achieved by doing a convolution between the *P*′(*W*′) and a normal distribution: *P*′(*W*′) there are infinitely many choices of Δ and Σ that yield a *P*_{Δ,Σ}(*W*) satisfying the CFR. Indeed, let us suppose that the pair Δ*, Σ* is such that *ϕ* (Fig. 6*B*). In this situation the inference cannot rest on the CFR alone. Explicit calculations in the Gaussian case (*SI Text*, section S6) show that variances (*σ*^{2}) and means (〈…〉) of *P*(*W*) and *P*′(*W*′) are related by an asymmetry factor (AF),*SI Text*, section S6). Knowing the AF allows us to select the unique pair Δ, Σ such that AF = Σ^{2}/Δ with *P*_{Δ,Σ}(*W*) satisfying the CFR. The inference procedure can be described with a very simple formula (*SI Text*, section S6). The key idea is to proceed as previously done in the case of symmetric setups by just shifting the mean of *P*′(*W*′) by a parameter *δ* until the CFR is satisfied; i.e., Δ = *δ*, Σ = 0 (Fig. 6*B*, *Center*). From the values of AF and *δ* we can reconstruct *P*(*W*) by using the formulas*C* and *D* for cyclic dsDNA pulling experiments. For noncyclic pulls (Δ*G* ≠ 0) the procedure can be easily generalized in the line of what has been shown for the case of the hairpin in the symmetric setup (Fig. 5*C*).

## Discussion

FRs are among the few general exact results in nonequilibrium statistical mechanics. Their validity has been already extensively tested in different systems, ranging from single molecules to single-electron transistors, and in different conditions (steady-state dynamics, irreversible transformations between steady states, and transient nonequilibrium states). At the present stage, the main widespread application of FR is free-energy recovery from nonequilibrium pulling experiments in the single-molecule field. What we are presenting here is an application of FR for inference. All FRs are statements about the statistics of the total entropy production in a system plus the environment. If some part of the entropy production is missed or inadequately considered, FRs will in general not hold. This is why, for irreversible transformations between equilibrium states, we have a FR for the dissipated work (which is the total entropy production) but not for the dissipated heat (which is just the entropy production in the environment). The main tenet is now that the violation of FRs in a given setting provides useful information: It is evidence that some contribution to the total entropy production is being missed. We have given rigorous examples in which the violation of FRs can be used to characterize the missing entropy production. Remarkably, in our model system, one could even replace the moving trap by a moving micropipette, an object lacking any measurement capability, and still infer the work distribution exerted by that object on the molecular system [this extremely asymmetric setup would still be described by Eq. **30**, with *k*_{A} → ∞ and AF = 4*K*_{B}*Tk*_{m}/(*k*_{B} + 2*k*_{m})]. These results open the exciting prospect of extending and applying these ideas to steady-state systems, such as molecular motors, to extract useful information about their mechanochemical cycle.

## Conclusions

To give a clear and definite meaning to free-energy inference we have discussed irreversible transformations between equilibrium states performed with dual-trap optical tweezers. In these experiments a molecular tether is attached between two beads that are manipulated with two optical traps. The irreversible transformation is performed by increasing the trap-to-trap distance at a finite speed. In this kind of transformation the dissipated work equals the total entropy production, leading to our first result: In pulling experiments work *W* must be defined on the force measured in the trap that is moved with respect to the thermal bath. The force measured in the trap at rest gives rise to a work definition, *W*′, which does not satisfy the FR and is unsuitable to extract free-energy differences. We have called *W*′ a partial work measurement because it misses part of the total dissipation. This result is of direct interest to experimentalists: Many optical tweezers setups are designed so that they can measure only *W*′. We have thus imagined a situation in which *W*′ is measurable and *W* is not and asked the question, Can we infer the distribution of *W* from that of *W*′? If the question is asked in full generality, without any system-specific information, the answer is probably negative. Knowing only the extent of violation of the FR will be of little use; in general, some additional system-specific information will be needed for a successful inference. Here we discussed free-energy inference in the framework of a Gaussian approximation; the extent to which such inference is generally possible should be the subject of future studies. Let us summarize our main results:

• A symmetry of the system can be crucial for the inference. For left–right symmetric systems (as exemplified in our symmetric dual-trap setup) the

*P*(*W*) can be inferred from*P*′(*W*′) just by imposing that the former satisfies the CFR. When symmetry considerations cannot be used, the knowledge of some equilibrium properties of the system may suffice to successfully guide the inference (such as the stiffnesses for the asymmetric setup).• The inference process can be used to recover the full dissipation spectrum plus additional information about the hidden entropy source. In our specific dual-trap example

*W*′ does not account for the dissipation due to the movement of the center of mass of the dumbbell and the inference procedure can be seen as a method to measure the associated hydrodynamic drag.• We stress the benefits of using symmetric dumbbells in single-molecule manipulation. In this case an alternative work definition, the differential work

*W*^{−}, fulfills the CFR and is thus suitable for free-energy measurements. Given that*W*^{−}is less influenced by dissipation than*W*, switching from*W*to*W*^{−}ensures faster convergence of unidirectional free-energy estimates. For asymmetric setups*W*^{−}does not satisfy a FR anymore (only*W*does) and cannot be used to extract free-energy differences.

A deep understanding of how to correctly define and measure thermodynamic work in small systems (a long debated question in the past 20 y) is not just a fine detail for experimentalists and theorists working in the single-molecule field, but an essential question pertaining to all areas of modern science concerned with energy transfer processes at the nanoscale. The added feature of free-energy inference discussed in this paper paves the way to apply FRs to new problems and contexts. This remains among the most interesting open problems in this exciting field.

## Methods

### Buffers and DNA Substrates.

All experiments were performed in PBS buffer, 1 M NaCl at 25 °C; 1 mg/mL BSA was added to passivate the surfaces and avoid nonspecific interactions. The dsDNA tether was obtained by ligating a 1-kb segment to a biotin-labeled oligo at one end and a digoxigenin labelled oligo at the other end. The DNA hairpin used in the experiments has short (20 bp) molecular handles and was synthesized by hybridization and ligating three different oligos. One oligo is biotin labeled and a second is dig labeled. Details of the synthesis procedure are given in ref. 33.

### Optical Tweezers Assay.

Measurements were performed with a highly stable miniaturized laser tweezers in the dual-trap mode (32). This instrument directly measures forces by linear momentum conservation. In all experiments we used silica beads with 4-μm diameter, which give a maximum trapping force around 20 pN. Data are acquired at 1 kHz.

## Acknowledgments

We thank A. Alemany and M. Palassini for a critical reading of the manuscript. F.R. is supported by an Institució Catalana de Recerca i Estudis Avançats Academia 2008 grant. The research leading to these results has received funding from the European Union Seventh Framework Programme (FP7/2007-2013) under Grant 308850 INFERNOS (information, fluctuations, and energy control in small systems).

## Footnotes

- ↵
^{1}To whom correspondence should be addressed. Email: fritort{at}gmail.com.

Author contributions: M.R.-C. and F.R. designed research; M.R.-C. performed research; M.R.-C. analyzed data; and M.R.-C. and F.R. wrote the paper.

The authors declare no conflict of interest.

↵*This Direct Submission article had a prearranged editor.

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1320006111/-/DCSupplemental.

## References

- ↵.
- Hummer G,
- Szabo A

- ↵.
- Hummer G,
- Szabo A

- ↵.
- Liphardt J,
- Dumont S,
- Smith SB,
- Tinoco I Jr,
- Bustamante C

- ↵
- ↵
- ↵
- ↵
- ↵.
- Dhakal S,
- et al.

- ↵.
- Cecconi C,
- Shank EA,
- Bustamante C,
- Marqusee S

- ↵
- ↵.
- Gebhardt JC,
- Bornschlögl T,
- Rief M

- ↵.
- Yu H,
- et al.

- ↵
- ↵
- ↵
- ↵.
- Douarche F,
- Ciliberto S,
- Petrosyan A

- ↵
- ↵.
- Klages R,
- Just W,
- Jarzynski C

- Alemany A,
- Ribezzi-Crivellari M,
- Ritort F

- ↵
- ↵.
- Peterman EJG,
- Wuite GJL

- Cisse I,
- Mangeol P,
- Bockelmann U

- ↵.
- Moffitt JR,
- Chemla YR,
- Izhaky D,
- Bustamante C

- ↵
- ↵.
- Mazonka O,
- Jarzynski C

- ↵
- ↵
- ↵.
- Gore J,
- Ritort F,
- Bustamante C

- ↵
- ↵.
- Ritort F,
- Bustamante C,
- Tinoco I Jr

- ↵
- ↵
- ↵
- ↵

## Citation Manager Formats

## Sign up for Article Alerts

## Article Classifications

- Biological Sciences
- Biophysics and Computational Biology

- Physical Sciences
- Physics