Work distribution in manipulated single biomolecules

We consider the relation between the microscopic and effective descriptions of the unfolding experiment on a model polypeptide. We evaluate the probability distribution function of the performed work by Monte Carlo simulations and compare it with that obtained by evaluating the work distribution generating function on an effective Brownian motion model tailored to reproduce exactly the equilibrium properties. The agreement is satisfactory for fast protocols, but deteriorates for slower ones, hinting at the existence of processes on several time scales even in such a simple system.

💡 Research Summary

The paper investigates how microscopic and effective descriptions compare in the context of single‑molecule pulling experiments on a model polypeptide. Using Monte Carlo simulations, the authors generate a large ensemble of non‑equilibrium unfolding trajectories in which an external force protocol stretches the chain. For each trajectory the work performed on the molecule is recorded, and from thousands of repetitions a probability distribution function (PDF) of the work is constructed. This “microscopic” PDF exhibits pronounced non‑Gaussian features, especially in the tails, reflecting the underlying complex energy landscape of the peptide.

In parallel, the authors construct a one‑dimensional effective Brownian‑motion model that is calibrated to reproduce exactly the equilibrium free‑energy profile of the same peptide. Within this reduced description the work generating function can be evaluated analytically (or numerically) and, by inverse Laplace transformation, the corresponding work PDF is obtained. Because the effective model is built to match the static thermodynamic properties, any discrepancy between the two PDFs must arise from dynamical differences.

The central comparison is performed for a set of pulling protocols ranging from very fast (tens of microseconds) to relatively slow (milliseconds). For fast protocols the system is driven far from equilibrium and the unfolding proceeds essentially in a single dominant step. Under these conditions the microscopic Monte Carlo PDF and the effective‑model PDF overlap almost perfectly; the reduced description captures the essential stochastic dynamics. As the pulling speed is reduced, however, the peptide has time to explore intermediate conformations, to undergo multiple barrier crossings, and to relax partially between successive events. These processes introduce several characteristic time scales that are not represented in the simple one‑dimensional Langevin equation. Consequently the work PDFs begin to diverge: the microscopic distribution broadens, its tails become more asymmetric, and the effective model underestimates the probability of large positive or negative work values.

The authors interpret this divergence as evidence that even in a minimal polypeptide system the unfolding dynamics cannot be fully captured by a single effective coordinate. Multi‑scale dynamics, memory effects, and non‑Markovian behavior become relevant for slower driving. This finding has practical implications for the application of non‑equilibrium work relations such as the Jarzynski equality. Since the equality relies on accurate sampling of the work distribution’s exponential average, insufficient representation of the tails (as occurs with the effective model at slow rates) can lead to systematic errors in the estimated free‑energy difference.

Overall, the study demonstrates that the validity of coarse‑grained stochastic models depends critically on the time scale of the external perturbation. Fast pulling experiments can be interpreted reliably with simple effective Brownian dynamics, whereas slow protocols demand a more detailed microscopic treatment that retains the multiple relaxation pathways of the molecule. The work therefore provides a clear benchmark for assessing when reduced models are adequate and when they must be supplemented by richer dynamical descriptions in single‑molecule biophysics.