Activated rate processes in the multidimensional case. Consideration of recrossings in the multidimensional Kramers problem with anisotropic friction

Mechanobiology: protein refolding under force

The application of direct force to a protein enables to probe wide regions of its energy surface through conformational transitions as unfolding, extending, recoiling, collapsing, and refolding. While unfolding under force typically displayed a two-state behavior, refolding under force, from highly extended unfolded states, displayed a more complex behavior. The first recording of protein refolding at a force quench step displayed an initial rapid elastic recoil, followed by a plateau phase at some extension, concluding with a collapse to a final state, at which refolding occurred. These findings stirred a lively discussion, which led to further experimental and theoretical investigation of this behavior. It was demonstrated that the polymeric chain of the unfolded protein is required to fully collapse to a globular conformation for the maturation of native structure. This behavior was modeled using one-dimensional free energy landscape over the end-to-end length reaction coordinate, the collective measured variable. However, at low forces, conformational space is not well captured by such models, and using two-dimensional energy surfaces provides further insight into the dynamics of this process. This work reviews the main concepts of protein refolding under constant force, which is essential for understanding how mechanotransducing proteins operate in vivo.

Transition paths in single-molecule force spectroscopy

In a typical single-molecule force spectroscopy experiment, the ends of the molecule of interest are connected by long polymer linkers to a pair of mesoscopic beads trapped in the focus of two laser beams. At constant force load, the total extension, i.e., the end-to-end distance of the molecule plus linkers, is measured as a function of time. In the simplest systems, the measured extension fluctuates about two values characteristic of folded and unfolded states, with occasional transitions between them. We have recently shown that molecular (un)folding rates can be recovered from such trajectories, with a small linker correction, as long as the characteristic time of the bead fluctuations is shorter than the residence time in the unfolded (folded) state. Here, we show that accurate measurements of the molecular transition path times require an even faster apparatus response. Transition paths, the trajectory segments in which the molecule (un)folds, are properly resolved only if the beads fluctuate more rapidly than the end-to-end distance of the molecule. Therefore, over a wide regime, the measured rates may be meaningful but not the transition path times. Analytic expressions for the measured mean transition path times are obtained for systems diffusing anisotropically on a two-dimensional free energy surface. The transition path times depend on the properties both of the molecule and of the pulling device.

Communication: Does force spectroscopy of biomolecules probe their intrinsic dynamic properties?

In single-molecule pulling experiments, the molecule of interest is attached to a much larger object such as an atomic force microscope tip or a micrometer sized bead. The measured dynamics of molecular transitions is therefore affected by the hydrodynamic drag on the pulling instrument itself. By considering the transitions within the combined system (the molecule and the instrument), it is shown here that two distinct physical regimes exist: when the intrinsic stiffness of the molecule is greater than that of the linker connecting the molecule to the pulling setup then the pulling experiment probes the intrinsic dynamics of the molecule with only relatively small (and quantifiable) corrections resulting from the pulling setup. In contrast, when the stiffness of the linker exceeds that of the molecule, the molecular transition in question involves concerted motion of the molecule and the pulling setup and the hydrodynamic drag on the pulling instrument becomes the dominant source of friction along the molecular reaction coordinate. An analytical formula interpolating between these two cases is further derived. These results explain recent conflicting observations where some single-molecule pulling measurements report anomalously low diffusion coefficients along molecular reaction coordinates while others do not.

Multidimensional reaction rate theory with anisotropic diffusion

An analytical expression is derived for the rate constant that describes diffusive transitions between two deep wells of a multidimensional potential. The expression, in contrast to the Kramers-Langer formula for the rate constant, is valid even when the diffusion is highly anisotropic. Our approach is based on a variational principle for the reactive flux and uses a trial function for the splitting probability or commitor. The theoretical result is validated by Brownian dynamics simulations.

On artifacts in single-molecule force spectroscopy

In typical force spectroscopy experiments, a small biomolecule is attached to a soft polymer linker that is pulled with a relatively large bead or cantilever. At constant force, the total extension stochastically changes between two (or more) values, indicating that the biomolecule undergoes transitions between two (or several) conformational states. In this paper, we consider the influence of the dynamics of the linker and mesoscopic pulling device on the force-dependent rate of the conformational transition extracted from the time dependence of the total extension, and the distribution of rupture forces in force-clamp and force-ramp experiments, respectively. For these different experiments, we derive analytic expressions for the observables that account for the mechanical response and dynamics of the pulling device and linker. Possible artifacts arise when the characteristic times of the pulling device and linker become comparable to, or slower than, the lifetimes of the metastable conformational states, and when the highly anharmonic regime of stretched linkers is probed at high forces. We also revisit the problem of relating force-clamp and force-ramp experiments, and identify a linker and loading rate-dependent correction to the rates extracted from the latter. The theory provides a framework for both the design and the quantitative analysis of force spectroscopy experiments by highlighting, and correcting for, factors that complicate their interpretation.

Extracting intrinsic dynamic parameters of biomolecular folding from single-molecule force spectroscopy experiments

Single-molecule studies in which a mechanical force is transmitted to the molecule of interest and the molecular extension or position is monitored as a function of time are versatile tools for probing the dynamics of protein folding, stepping of molecular motors, and other biomolecular processes involving activated barrier crossing. One complication in interpreting such studies, however, is the fact that the typical size of a force probe (e.g., a dielectric bead in optical tweezers or the atomic force microscope tip/cantilever assembly) is much larger than the molecule itself, and so the observed molecular motion is affected by the hydrodynamic drag on the probe. This presents the experimenter with a nontrivial task of deconvolving the intrinsic molecular parameters, such as the intrinsic free energy barrier and the effective diffusion coefficient exhibited while crossing the barrier from the experimental signal. Here we focus on the dynamical aspect of this task and show how the intrinsic diffusion coefficient along the molecular reaction coordinate can be inferred from single-molecule measurements of the rates of biomolecular folding and unfolding. We show that the feasibility of accomplishing this task is strongly dependent on the relationship between the intrinsic molecular elasticity and that of the linker connecting the molecule to the force probe and identify the optimal range of instrumental parameters allowing determination of instrument-free molecular dynamics.

Free energy profiles from single-molecule pulling experiments

Nonequilibrium pulling experiments provide detailed information about the thermodynamic and kinetic properties of molecules. We show that unperturbed free energy profiles as a function of molecular extension can be obtained rigorously from such experiments without using work-weighted position histograms. An inverse Weierstrass transform is used to relate the system free energy obtained from the Jarzynski equality directly to the underlying molecular free energy surface. An accurate approximation for the free energy surface is obtained by using the method of steepest descent to evaluate the inverse transform. The formalism is applied to simulated data obtained from a kinetic model of RNA folding, in which the dynamics consists of jumping between linker-dominated folded and unfolded free energy surfaces.

Optimal flexibility for conformational transitions in macromolecules

Conformational transitions in macromolecular complexes often involve the reorientation of leverlike structures. Using a simple theoretical model, we show that the rate of such transitions is drastically enhanced if the lever is bendable, e.g., at a localized hinge. Surprisingly, the transition is fastest with an intermediate flexibility of the hinge. In this intermediate regime, the transition rate is also least sensitive to the amount of "cargo" attached to the lever arm, which could be exploited by molecular motors. To explain this effect, we generalize the Kramers-Langer theory for multidimensional barrier crossing to configuration-dependent mobility matrices.

Reaction rate theory: what it was, where is it today, and where is it going?

A brief history is presented, outlining the development of rate theory during the past century. Starting from Arrhenius [Z. Phys. Chem. 4, 226 (1889)], we follow especially the formulation of transition state theory by Wigner [Z. Phys. Chem. Abt. B 19, 203 (1932)] and Eyring [J. Chem. Phys. 3, 107 (1935)]. Transition state theory (TST) made it possible to obtain quick estimates for reaction rates for a broad variety of processes even during the days when sophisticated computers were not available. Arrhenius' suggestion that a transition state exists which is intermediate between reactants and products was central to the development of rate theory. Although Wigner gave an abstract definition of the transition state as a surface of minimal unidirectional flux, it took almost half of a century until the transition state was precisely defined by Pechukas [Dynamics of Molecular Collisions B, edited by W. H. Miller (Plenum, New York, 1976)], but even this only in the realm of classical mechanics. Eyring, considered by many to be the father of TST, never resolved the question as to the definition of the activation energy for which Arrhenius became famous. In 1978, Chandler [J. Chem. Phys. 68, 2959 (1978)] finally showed that especially when considering condensed phases, the activation energy is a free energy, it is the barrier height in the potential of mean force felt by the reacting system. Parallel to the development of rate theory in the chemistry community, Kramers published in 1940 [Physica (Amsterdam) 7, 284 (1940)] a seminal paper on the relation between Einstein's theory of Brownian motion [Einstein, Ann. Phys. 17, 549 (1905)] and rate theory. Kramers' paper provided a solution for the effect of friction on reaction rates but left us also with some challenges. He could not derive a uniform expression for the rate, valid for all values of the friction coefficient, known as the Kramers turnover problem. He also did not establish the connection between his approach and the TST developed by the chemistry community. For many years, Kramers' theory was considered as providing a dynamic correction to the thermodynamic TST. Both of these questions were resolved in the 1980s when Pollak [J. Chem. Phys. 85, 865 (1986)] showed that Kramers' expression in the moderate to strong friction regime could be derived from TST, provided that the bath, which is the source of the friction, is handled at the same level as the system which is observed. This then led to the Mel'nikov-Pollak-Grabert-Hanggi [Mel'nikov and Meshkov, J. Chem. Phys. 85, 1018 (1986); Pollak, Grabert, and Hanggi, ibid. 91, 4073 (1989)] solution of the turnover problem posed by Kramers. Although classical rate theory reached a high level of maturity, its quantum analog leaves the theorist with serious challenges to this very day. As noted by Wigner [Trans. Faraday Soc. 34, 29 (1938)], TST is an inherently classical theory. A definite quantum TST has not been formulated to date although some very useful approximate quantum rate theories have been invented. The successes and challenges facing quantum rate theory are outlined. An open problem which is being investigated intensively is rate theory away from equilibrium. TST is no longer valid and cannot even serve as a conceptual guide for understanding the critical factors which determine rates away from equilibrium. The nonequilibrium quantum theory is even less well developed than the classical, and suffers from the fact that even today, we do not know how to solve the real time quantum dynamics for systems with "many" degrees of freedom.

Multidimensional and memory effects on diffusion of a particle

The diffusion of an overdamped Brownian particle in the two-dimensional (2D) channel bounded periodically by a parabola is studied, where the particle is subject to an additive white or colored noise. The diffusion rate constant D* of the particle is evaluated by the quasi-2D approximation and the effective potential approach, and the theoretical result is compared with the Langevin simulation. The properties of the diffusion rate constant are stressed for weak and strong noise cases. It is shown that, in an entropy channel, the value of D* in units of Q decreases with increasing intensity of the colored noise. In the presence of energetic barriers, a nonmonotonic behavior of the reduced diffusion rate constant D*Q-1 as a function of the noise intensity is shown.