Broad distributions of transition-path times are fingerprints of multidimensionality of the underlying free energy landscapes

Accelerating Molecular Dynamics through Informed Resetting

We present a procedure for enhanced sampling of molecular dynamics simulations through informed stochastic resetting. Many phenomena, such as protein folding and crystal nucleation, occur over time scales inaccessible in standard simulations. We recently showed that stochastic resetting can accelerate molecular simulations that exhibit broad transition time distributions. However, standard stochastic resetting does not exploit any information about the reaction progress. For a model system and chignolin in explicit water, we demonstrate that an informed resetting protocol leads to greater accelerations than standard stochastic resetting in molecular dynamics and Metadynamics simulations. This is achieved by resetting only when a certain condition is met, e.g., when the distance from the target along the reaction coordinate is larger than some threshold. We use these accelerated simulations to infer important kinetic observables such as the unbiased mean first-passage time and direct transit time. For the latter, Metadynamics with informed resetting leads to speedups of 2-3 orders of magnitude over unbiased simulations with relative errors of only ∼35-70%. Our work significantly extends the applicability of stochastic resetting for enhanced sampling of molecular simulations.

Transition Path Dynamics of Non-Markovian Systems across a Rough Potential Barrier

Transition paths refer to rare events in physics, chemistry, and biology where the molecules cross barriers separating stable molecular conformations. The conventional analysis of the transition path times employs a diffusive and memoryless transition over a smooth potential barrier. However, it is widely acknowledged that the free energy profile between two minima in biomolecular processes is inherently not smooth. In this article, we discuss a theoretical model with a parabolic rough potential barrier and obtain analytical results of the transition path distribution and mean transition path times by incorporating absorbing boundary conditions across the boundaries under the driving of Gaussian white noise. Further, the influence of anomalous dynamics in rough potential driven by a power-law memory kernel is analyzed by deriving a time-dependent scaled diffusion coefficient that coarse-grains the effects of roughness, and the system's dynamics is reduced to a scaled diffusion on a smooth potential. Our theoretical results are tested and validated against numerical simulations. The findings of our study show the influence of the boundary conditions, barrier height, barrier roughness, and memory effect on the transition path time distributions in a rough potential, and the validity of the scaling diffusion coefficient has been discussed.

The significance of fuzzy boundaries of the barrier regions in single-molecule measurements of failed barrier crossing attempts

A recent ground-breaking experimental study [Lyons et al., Phys. Rev. X 14(1), 011017 (2024)] reports on measuring the temporal duration and the spatial extent of failed attempts to cross an activation barrier (i.e., "loops") for a folding transition in a single molecule and for a Brownian particle trapped within a bistable potential. Within the model of diffusive dynamics, however, both of these quantities are, on average, exactly zero because of the recrossings of the barrier region boundary. That is, an observer endowed with infinite spatial and temporal resolution would find that finite loops do not exist (or, more precisely, form a set of measure zero). Here we develop a description of the experiment that takes the "fuzziness" of the boundaries caused by finite experimental resolution into account and show how the experimental uncertainty of localizing the point, in time and space, where the barrier is crossed leads to observable distributions of loop times and sizes. Although these distributions generally depend on the experimental resolution, this dependence, in certain cases, may amount to a simple resolution-dependent factor and, therefore, the experiments do probe inherent properties of barrier crossing dynamics.

General theory for localizing the where and when of entropy production meets single-molecule experiments

The laws of thermodynamics apply to biophysical systems on the nanoscale as described by the framework of stochastic thermodynamics. This theory provides universal, exact relations for quantities like work, which have been verified in experiments where a fully resolved description allows direct access to such quantities. Complementary studies consider partially hidden, coarse-grained descriptions, in which the mean entropy production typically is not directly accessible but can be bounded in terms of observable quantities. Going beyond the mean, we introduce a fluctuating entropy production that applies to individual trajectories in a coarse-grained description under time-dependent driving. Thus, this concept is applicable to the broad and experimentally significant class of driven systems in which not all relevant states can be resolved. We provide a paradigmatic example by studying an experimentally verified protein unfolding process. As a consequence, the entire distribution of the coarse-grained entropy production rather than merely its mean retains spatial and temporal information about the microscopic process. In particular, we obtain a bound on the distribution of the physical entropy production of individual unfolding events.

The significance of fuzzy boundaries of the barrier regions in single-molecule measurements of failed barrier crossing attempts

A recent experimental study reports on measuring the temporal duration and the spatial extent of failed attempts to cross an activation barrier (i.e., "loops") for a folding transition in a single molecule and for a Brownian particle trapped within a bistable potential. Within the model of diffusive dynamics, however, both of these quantities are, on the average, exactly zero because of the recrossings of the barrier region boundary. That is, an observer endowed with infinite spatial and temporal resolution would find that finite loops do not exist (or, more precisely, form a set of measure zero). Here we develop a description of the experiment that takes finite experimental resolution into account and show how the experimental uncertainty of localizing the point, in time and space, where the barrier is crossed leads to observable distributions of loop times and sizes. Although these distributions generally depend on the experimental resolution, this dependence, in certain cases, may amount to a simple resolution-dependent factor and thus the experiments do probe inherent properties of barrier crossing dynamics.

Assessing eco-geographic influences on COVID-19 transmission: a global analysis

COVID-19 has been massively transmitted for almost 3 years, and its multiple variants have caused serious health problems and an economic crisis. Our goal was to identify the influencing factors that reduce the threshold of disease transmission and to analyze the epidemiological patterns of COVID-19. This study served as an early assessment of the epidemiological characteristics of COVID-19 using the MaxEnt species distribution algorithm using the maximum entropy model. The transmission of COVID-19 was evaluated based on human factors and environmental variables, including climate, terrain and vegetation, along with COVID-19 daily confirmed case location data. The results of the SDM model indicate that population density was the major factor influencing the spread of COVID-19. Altitude, land cover and climatic factor showed low impact. We identified a set of practical, high-resolution, multi-factor-based maximum entropy ecological niche risk prediction systems to assess the transmission risk of the COVID-19 epidemic globally. This study provided a comprehensive analysis of various factors influencing the transmission of COVID-19, incorporating both human and environmental variables. These findings emphasize the role of different types of influencing variables in disease transmission, which could have implications for global health regulations and preparedness strategies for future outbreaks.

Time-Resolved Statistics of Snippets as General Framework for Model-Free Entropy Estimators

Irreversibility is commonly quantified by entropy production. An external observer can estimate it through measuring an observable that is antisymmetric under time reversal like a current. We introduce a general framework that allows us to infer a lower bound on entropy production through measuring the time-resolved statistics of events with any symmetry under time reversal, in particular, time-symmetric instantaneous events. We emphasize Markovianity as a property of certain events rather than of the full system and introduce an operationally accessible criterion for this weakened Markov property. Conceptually, the approach is based on snippets as particular sections of trajectories between two Markovian events, for which a generalized detailed balance relation is discussed.

The effect of time resolution on the observed first passage times in diffusive dynamics

Single-molecule and single-particle tracking experiments are typically unable to resolve fine details of thermal motion at short timescales where trajectories are continuous. We show that, when a diffusive trajectory xt is sampled at finite time intervals δt, the resulting error in measuring the first passage time to a given domain can exceed the time resolution of the measurement by more than an order of magnitude. Such surprisingly large errors originate from the fact that the trajectory may enter and exit the domain while being unobserved, thereby lengthening the apparent first passage time by an amount that is larger than δt. Such systematic errors are particularly important in single-molecule studies of barrier crossing dynamics. We show that the correct first passage times, as well as other properties of the trajectories such as splitting probabilities, can be recovered via a stochastic algorithm that reintroduces unobserved first passage events probabilistically.

Fick-Jacobs description and first passage dynamics for diffusion in a channel under stochastic resetting

The transport of particles through channels is of paramount importance in physics, chemistry, and surface science due to its broad real world applications. Much insight can be gained by observing the transition paths of a particle through a channel and collecting statistics on the lifetimes in the channel or the escape probabilities from the channel. In this paper, we consider the diffusive transport through a narrow conical channel of a Brownian particle subject to intermittent dynamics, namely, stochastic resetting. As such, resetting brings the particle back to a desired location from where it resumes its diffusive phase. To this end, we extend the Fick-Jacobs theory of channel-facilitated diffusive transport to resetting-induced transport. Exact expressions for the conditional mean first passage times, escape probabilities, and the total average lifetime in the channel are obtained, and their behavior as a function of the resetting rate is highlighted. It is shown that resetting can expedite the transport through the channel-rigorous constraints for such conditions are then illustrated. Furthermore, we observe that a carefully chosen resetting rate can render the average lifetime of the particle inside the channel minimal. Interestingly, the optimal rate undergoes continuous and discontinuous transitions as some relevant system parameters are varied. The validity of our one-dimensional analysis and the corresponding theoretical predictions is supported by three-dimensional Brownian dynamics simulations. We thus believe that resetting can be useful to facilitate particle transport across biological membranes-a phenomenon that can spearhead further theoretical and experimental studies.

Challenges in Inferring the Directionality of Active Molecular Processes from Single-Molecule Fluorescence Resonance Energy Transfer Trajectories

We discuss some of the practical challenges that one faces in using stochastic thermodynamics to infer directionality of molecular machines from experimental single-molecule trajectories. Because of the limited spatiotemporal resolution of single-molecule experiments and because both forward and backward transitions between the same pairs of states cannot always be detected, differentiating between the forward and backward directions of, e.g., an ATP-consuming molecular machine that operates periodically, turns out to be a nontrivial task. Using a simple extension of a Markov-state model that is commonly employed to analyze single-molecule transition-path measurements, we illustrate how irreversibility can be hidden from such measurements but in some cases can be uncovered when non-Markov effects in low-dimensional single-molecule trajectories are considered.

Computing committors via Mahalanobis diffusion maps with enhanced sampling data

The study of phenomena such as protein folding and conformational changes in molecules is a central theme in chemical physics. Molecular dynamics (MD) simulation is the primary tool for the study of transition processes in biomolecules, but it is hampered by a huge timescale gap between the processes of interest and atomic vibrations that dictate the time step size. Therefore, it is imperative to combine MD simulations with other techniques in order to quantify the transition processes taking place on large timescales. In this work, the diffusion map with Mahalanobis kernel, a meshless approach for approximating the Backward Kolmogorov Operator (BKO) in collective variables, is upgraded to incorporate standard enhanced sampling techniques, such as metadynamics. The resulting algorithm, which we call the target measure Mahalanobis diffusion map (tm-mmap), is suitable for a moderate number of collective variables in which one can approximate the diffusion tensor and free energy. Imposing appropriate boundary conditions allows use of the approximated BKO to solve for the committor function and utilization of transition path theory to find the reactive current delineating the transition channels and the transition rate. The proposed algorithm, tm-mmap, is tested on the two-dimensional Moro-Cardin two-well system with position-dependent diffusion coefficient and on alanine dipeptide in two collective variables where the committor, the reactive current, and the transition rate are compared to those computed by the finite element method (FEM). Finally, tm-mmap is applied to alanine dipeptide in four collective variables where the use of finite elements is infeasible.

Single-molecule mechanical studies of chaperones and their clients

Single-molecule force spectroscopy provides access to the mechanics of biomolecules. Recently, magnetic and laser optical tweezers were applied in the studies of chaperones and their interaction with protein clients. Various aspects of the chaperone-client interactions can be revealed based on the mechanical probing strategies. First, when a chaperone is probed under load, one can examine the inner workings of the chaperone while it interacts with and works on the client protein. Second, when protein clients are probed under load, the action of chaperones on folding clients can be studied in great detail. Such client folding studies have given direct access to observing actions of chaperones in real-time, like foldase, unfoldase, and holdase activity. In this review, we introduce the various single molecule mechanical techniques and summarize recent single molecule mechanical studies on heat shock proteins, chaperone-mediated folding on the ribosome, SNARE folding, and studies of chaperones involved in the folding of membrane proteins. An outlook on significant future developments is given.

Microscopic origin of diffusive dynamics in the context of transition path time distributions for protein folding and unfolding

Experimentally measured transition path time distributions are usually analyzed theoretically in terms of a diffusion equation over a free energy barrier. It is though well understood that the free energy profile separating the folded and unfolded states of a protein is characterized as a transition through many stable micro-states which exist between the folded and unfolded states. Why is it then justified to model the transition path dynamics in terms of a diffusion equation, namely the Smoluchowski equation (SE)? In principle, van Kampen has shown that a nearest neighbor Markov chain of thermal jumps between neighboring microstates will lead in a continuum limit to the SE, such that the friction coefficient is proportional to the mean residence time in each micro-state. However, the practical question of how many microstates are needed to justify modeling the transition path dynamics in terms of an SE has not been addressed. This is a central topic of this paper where we compare numerical results for transition paths based on the diffusion equation on the one hand and the nearest neighbor Markov jump model on the other. Comparison of the transition path time distributions shows that one needs at least a few dozen microstates to obtain reasonable agreement between the two approaches. Using the Markov nearest neighbor model one also obtains good agreement with the experimentally measured transition path time distributions for a DNA hairpin and PrP protein. As found previously when using the diffusion equation, the Markov chain model used here also reproduces the experimentally measured long time tail and confirms that the transition path barrier height is ∼3k_BT. This study indicates that in the future, when attempting to model experimentally measured transition path time distributions, one should perhaps prefer a nearest neighbor Markov model which is well defined also for rough energy landscapes. Such studies can also shed light on the minimal number of microstates needed to unravel the experimental data.

The Effect of Time Resolution on Apparent Transition Path Times Observed in Single-Molecule Studies of Biomolecules

Single-molecule experiments have now achieved a time resolution allowing observation of transition paths, the brief trajectory segments where the molecule undergoing an unfolding or folding transition enters the energetically or entropically unfavorable barrier region from the folded/unfolded side and exits to the unfolded/folded side, thereby completing the transition. This resolution, however, is yet insufficient to identify the precise entrance/exit events that mark the beginning and the end of a transition path: the nature of the diffusive dynamics is such that a molecular trajectory will recross the boundary between the barrier region and the folded/unfolded state, multiple times, at a time scale much shorter than that of the typical experimental resolution. Here we use theory and Brownian dynamics simulations to show that, as a result of such recrossings, the apparent transition path times are generally longer than the true ones. We quantify this effect using a simple model where the observed dynamics is a moving average of the true dynamics and discuss experimental implications of our results.

Direct Observation of Single-Protein Transition State Passage by Nanopore Ionic Current Jumps

Conformational transitions of proteins are governed by chemical kinetics, often toggled by passage through an activated state separating two conformational ensembles. The passage time of a protein through the activated state can be too fast to be detected by single-molecule experiments without the aid of viscogenic agents. Here, we use high-bandwidth nanopore measurements to resolve microsecond-duration transitions that occur between conformational states of individual protein molecules partly blocking pore current. We measure the transition state passage time between folded and unfolded states of a two-state λ_6-85 mutant and between metastable intermediates and the unfolded state of the multistate folder cytochrome c. Consistent with the principle of microscopic reversibility, the transition state passage time is the same for the forward and backward reactions. A passage time distribution whose tail is broader than a single exponential observed in cytochrome c suggests a multidimensional energy landscape for this protein.

Heterogeneity in Protein Folding and Unfolding Reactions

Proteins have dynamic structures that undergo chain motions on time scales spanning from picoseconds to seconds. Resolving the resultant conformational heterogeneity is essential for gaining accurate insight into fundamental mechanistic aspects of the protein folding reaction. The use of high-resolution structural probes, sensitive to population distributions, has begun to enable the resolution of site-specific conformational heterogeneity at different stages of the folding reaction. Different states populated during protein folding, including the unfolded state, collapsed intermediate states, and even the native state, are found to possess significant conformational heterogeneity. Heterogeneity in protein folding and unfolding reactions originates from the reduced cooperativity of various kinds of physicochemical interactions between various structural elements of a protein, and between a protein and solvent. Heterogeneity may arise because of functional or evolutionary constraints. Conformational substates within the unfolded state and the collapsed intermediates that exchange at rates slower than the subsequent folding steps give rise to heterogeneity on the protein folding pathways. Multiple folding pathways are likely to represent distinct sequences of structure formation. Insight into the nature of the energy barriers separating different conformational states populated during (un)folding can also be obtained by resolving heterogeneity.

On distributions of barrier crossing times as observed in single-molecule studies of biomolecules

Single-molecule experiments that monitor time evolution of molecular observables in real time have expanded beyond measuring transition rates toward measuring distributions of times of various molecular events. Of particular interest is the first-passage time for making a transition from one molecular configuration ( a ) to another ( b ) and conditional first-passage times such as the transition path time, which is the first-passage time from a to b conditional upon not leaving the transition region intervening between a and b . Another experimentally accessible (but not yet studied experimentally) observable is the conditional exit time, i.e., the time to leave the transition region through a specified boundary. The distributions of such times contain a wealth of mechanistic information about the transitions in question. Here, we use the first and the second (and, if desired, higher) moments of these distributions to characterize their relative width for the model in which the experimental observable undergoes Brownian motion in a potential of mean force. We show that although the distributions of transition path times are always narrower than exponential (in that the ratio of the standard deviation to the distribution's mean is always less than 1), distributions of first-passage times and of conditional exit times can be either narrow or broad, in some cases displaying long power-law tails. The conditional exit time studied here provides a generalization of the transition path time that also allows one to characterize the temporal scales of failed barrier crossing attempts.

Protein Folding Dynamics as Diffusion on a Free Energy Surface: Rate Equation Terms, Transition Paths, and Analysis of Single-Molecule Photon Trajectories

The rates of protein (un)folding are often described as diffusion on the projection of a hyperdimensional energy landscape onto a few (ideally one) order parameters. Testing such an approximation by experiment requires resolving the reactive transition paths of individual molecules, which is now becoming feasible with advanced single-molecule spectroscopic techniques. This has also sparked the interest of theorists in better understanding reactive transition paths. Here we focus on these issues aiming to establish (i) practical guidelines for the mechanistic interpretation of transition path times (TPT) and (ii) methods to extract the free energy surface and protein dynamics from the maximum likelihood analysis of photon trajectories (MLA-PT). We represent the (un)folding rates as diffusion on a 1D free energy surface with the FRET efficiency as a reaction coordinate proxy. We then perform diffusive kinetic simulations on surfaces with two minima and a barrier, but with different shapes (curvatures, barrier height, and symmetry), coupled to stochastic simulations of photon emissions that reproduce current SM-FRET experiments. From the analysis of transition paths, we find that the TPT is inversely proportional to the barrier height (difference in free energy between minimum and barrier top) for any given surface shape, and that dividing the TPT into climb and descent segments provides key information about the barrier's symmetry. We also find that the original MLA-PT procedure used to determine the TPT from experiments underestimates its value, particularly for the cases with smaller barriers (e.g., fast folders), and we suggest a simple strategy to correct for this bias. Importantly, we also demonstrate that photon trajectories contain enough information to extract the 1D free energy surface's shape and dynamics (if TPT is >4-5-fold longer than the interphoton time) using the MLA-PT directly implemented with a diffusive free energy surface model. When dealing with real (unknown) experimental data, the comparison between the likelihoods of the free energy surface and discrete kinetic three-state models can be used to evaluate the statistical significance of the estimated free energy surface.

Exposure of Von Willebrand Factor Cleavage Site in A1A2A3-Fragment under Extreme Hydrodynamic Shear

Von Willebrand Factor (vWf) is a giant multimeric extracellular blood plasma involved in hemostasis. In this work we present multi-scale simulations of its three-domains fragment A1A2A3. These three domains are essential for the functional regulation of vWf. Namely the A2 domain hosts the site where the protease ADAMTS13 cleavages the multimeric vWf allowing for its length control that prevents thrombotic conditions. The exposure of the cleavage site follows the elongation/unfolding of the domain that is caused by an increased shear stress in blood. By deploying Lattice Boltzmann molecular dynamics simulations based on the OPEP coarse-grained model for proteins, we investigated at molecular level the unfolding of the A2 domain under the action of a perturbing shear flow. We described the structural steps of this unfolding that mainly concerns the β-strand structures of the domain, and we compared the process occurring under shear with that produced by the action of a directional pulling force, a typical condition of single molecule experiments. We observe, that under the action of shear flow, the competition among the elongational and rotational components of the fluid field leads to a complex behaviour of the domain, where elongated structures can be followed by partially collapsed melted globule structures with a very different degree of exposure of the cleavage site. Our simulations pose the base for the development of a multi-scale in-silico description of vWf dynamics and functionality in physiological conditions, including high resolution details for molecular relevant events, e.g., the binding to platelets and collagen during coagulation or thrombosis.

What can we learn from transition path time distributions for protein folding and unfolding?

Recent advances in experimental measurements of transition path time distributions have raised intriguing theoretical questions. The present interpretation of the experimental data indicates a small value of the fitted transition path barrier height as compared to the barrier height of the unfolded to folded transition. Secondly, as shown in this paper, it is essential to analyse the experimental data using absorbing boundary conditions at the end points used to determine the transition paths. Such an analysis reveals long time tails that have thus far eluded quantitative theoretical interpretation. Is this due to uncertainty in the experimental data or does it call for a rethinking of the theoretical interpretation? A detailed study of the transition path time distribution using a diffusive model leads to the following conclusions. a. The present experimental data is not accurate enough to discern between functional forms of bell shaped free energy barriers. b. Long time tails indicate the possible existence of a "trap" in the transition path region. c. The "trap" may be considered as a well in the free energy surface. d. The long time tail is quite sensitive to the form of the trap so that future measurements of the long time tail as a function of the location of the end points of the transition path may make it possible to not only determine the well depth but also to distinguish between different functional forms for the free energy surface. e. Introduction of a well along the transition path leads to good fits with the experimental data provided that the transition path barrier height is ∼3k_BT, substantially higher than the estimates of ∼1k_BT based on bell shaped functions. The results presented here negate the need of introducing multi-dimensional effects, free energy barrier asymmetry, sub-diffusive memory kernels or systematic ruggedness to explain the experimentally measured data.

Nonparametric Analysis of Nonequilibrium Simulations

We extend the nonparametric framework of reaction coordinate optimization to nonequilibrium ensembles of (short) trajectories. For example, we show how, starting from such an ensemble, one can obtain an equilibrium free-energy profile along the committor, which can be used to determine important properties of the dynamics exactly. A new adaptive sampling approach, the transition-state ensemble enrichment, is suggested, which samples the configuration space by "growing" committor segments toward each other starting from the boundary states. This framework is suggested as a general tool, alternative to the Markov state models, for a rigorous and accurate analysis of simulations of large biomolecular systems, as it has the following attractive properties. It is immune to the curse of dimensionality, does not require system-specific information, can approximate arbitrary reaction coordinates with high accuracy, and has sensitive and rigorous criteria to test optimality and convergence. The approaches are illustrated on a 50-dimensional model system and a realistic protein folding trajectory.

Energy Landscape of Ubiquitin Is Weakly Multidimensional

Single molecule pulling experiments report time-dependent changes in the extension (X) of a biomolecule as a function of the applied force (f). By fitting the data to one-dimensional analytical models of the energy landscape, we can extract the hopping rates between the folded and unfolded states in two-state folders as well as the height and the location of the transition state (TS). Although this approach is remarkably insightful, there are cases for which the energy landscape is multidimensional (catch bonds being the most prominent). To assess if the unfolding energy landscape in small single domain proteins could be one-dimensional, we simulated force-induced unfolding of ubiquitin (Ub) using the coarse-grained self-organized polymer-side chain (SOP-SC) model. Brownian dynamics simulations using the SOP-SC model reveal that the Ub energy landscape is weakly multidimensional (WMD), governed predominantly by a single barrier. The unfolding pathway is confined to a narrow reaction pathway that could be described as diffusion in a quasi-1D X-dependent free energy profile. However, a granular analysis using the P_fold analysis, which does not assume any form for the reaction coordinate, shows that X alone does not account for the height and, more importantly, the location of the TS. The f-dependent TS location moves toward the folded state as f increases, in accord with the Hammond postulate. Our study shows that, in addition to analyzing the f-dependent hopping rates, the transition state ensemble must also be determined without resorting to X as a reaction coordinate to describe the unfolding energy landscapes of single domain proteins, especially if they are only WMD.

Reconstruction of mechanical unfolding and refolding pathways of proteins with atomic force spectroscopy and computer simulations

Most proteins in proteomes are large, typically consist of more than one domain and are structurally complex. This often makes studying their mechanical unfolding pathways challenging. Proteins composed of tandem repeat domains are a subgroup of multi-domain proteins that, when stretched, display a saw-tooth pattern in their mechanical unfolding force extension profiles due to their repetitive structure. However, the assignment of force peaks to specific repeats undergoing mechanical unraveling is complicated because all repeats are similar and they interact with their neighbors and form a contiguous tertiary structure. Here, we describe in detail a combination of experimental and computational single-molecule force spectroscopy methods that proved useful for examining the mechanical unfolding and refolding pathways of ankyrin repeat proteins. Specifically, we explain and delineate the use of atomic force microscope-based single molecule force spectroscopy (SMFS) to record the mechanical unfolding behavior of ankyrin repeat proteins and capture their unusually strong refolding propensity that is responsible for generating impressive refolding force peaks. We also describe Coarse Grain Steered Molecular Dynamic (CG-SMD) simulations which complement the experimental observations and provide insights in understanding the unfolding and refolding of these proteins. In addition, we advocate the use of novel coiled-coils-based mechanical polypeptide probes which we developed to demonstrate the vectorial character of folding and refolding of these repeat proteins. The combination of AFM-based SMFS on native and CC-equipped proteins with CG-SMD simulations is powerful not only for ankyrin repeat polypeptides, but also for other repeat proteins and more generally to various multidomain, non-repetitive proteins with complex topologies.

From folding to function: complex macromolecular reactions unraveled one-by-one with optical tweezers

Single-molecule manipulation with optical tweezers has uncovered macromolecular behaviour hidden to other experimental techniques. Recent instrumental improvements have made it possible to expand the range of systems accessible to optical tweezers. Beyond focusing on the folding and structural changes of isolated single molecules, optical tweezers studies have evolved into unraveling the basic principles of complex molecular processes such as co-translational folding on the ribosome, kinase activation dynamics, ligand-receptor binding, chaperone-assisted protein folding, and even dynamics of intrinsically disordered proteins (IDPs). In this mini-review, we illustrate the methodological principles of optical tweezers before highlighting recent advances in studying complex protein conformational dynamics - from protein synthesis to physiological function - as well as emerging future issues that are beginning to be addressed with novel approaches.

Localized potential well vs binding site: Mapping solute dynamics in a membrane channel onto one-dimensional description

In the one-dimensional description, the interaction of a solute molecule with the channel wall is characterized by the potential of mean force U(x), where the x-coordinate is measured along the channel axis. When the molecule can reversibly bind to certain amino acid(s) of the protein forming the channel, this results in a localized well in the potential U(x). Alternatively, this binding can be modeled by introducing a discrete localized site, in addition to the continuum of states along x. Although both models may predict identical equilibrium distributions of the coordinate x, there is a fundamental difference between the two: in the first model, the molecule passing through the channel unavoidably visits the potential well, while in the latter, it may traverse the channel without being trapped at the discrete site. Here, we show that when the two models are parameterized to have the same thermodynamic properties, they automatically yield identical translocation probabilities and mean translocation times, yet they predict qualitatively different shapes of the translocation time distribution. Specifically, the potential well model yields a narrower distribution than the model with a discrete site, a difference that can be quantified by the distribution's coefficient of variation. This coefficient turns out to be always smaller than unity in the potential well model, whereas it may exceed unity when a discrete trapping site is present. Analysis of the translocation time distribution beyond its mean thus offers a way to differentiate between distinct translocation mechanisms.

Estimating transition path times and shapes from single-molecule photon trajectories: A simulation analysis

In a two-state molecular system, transition paths comprise the portions of trajectories during which the system transits from one stable state to the other. Because of their low population, it is essentially impossible to obtain information on transition paths from experiments on a large sample of molecules. However, single-molecule experiments such as laser optical tweezers or Förster resonance energy transfer (FRET) spectroscopy have allowed transition-path durations to be estimated. Here, we use molecular simulations to test the methodology for obtaining information on transition paths in single-molecule FRET by generating photon trajectories from the distance trajectories obtained in the simulation. Encouragingly, we find that this maximum likelihood analysis yields transition-path times within a factor of 2-4 of the values estimated using a good coordinate for folding, but tends to systematically underestimate them. The underestimation can be attributed partly to the fact that the large changes in the end-end distance occur mostly early in a folding trajectory. However, even if the transfer efficiency is a good reaction coordinate for folding, the assumption that the transition-path shape is a step function still leads to an underestimation of the transition-path time as defined here. We find that allowing more flexibility in the form of the transition path model allows more accurate transition-path times to be extracted and points the way toward further improvements in methods for estimating transition-path time and transition-path shape.

Barrier Crossing Dynamics from Single-Molecule Measurements

Chemists visualize chemical reactions as motion along one-dimensional "reaction coordinates" over free energy barriers. Various rate theories, such as transition state theory and the Kramers theory of diffusive barrier crossing, differ in their assumptions regarding the mathematical specifics of this motion. Direct experimental observation of the motion along reaction coordinates requires single-molecule experiments performed with unprecedented time resolution. Toward this goal, recent single-molecule studies achieved time resolution sufficient to catch biomolecules in the act of crossing free energy barriers as they fold, bind to their targets, or undergo other large structural changes, offering a window into the elusive reaction "mechanisms". This Perspective describes what we can learn (and what we have already learned) about barrier crossing dynamics through synergy of single-molecule experiments, theory, and molecular simulations. In particular, I will discuss how emerging experimental data can be used to answer several questions of principle. For example, is motion along the reaction coordinate diffusive, is there conformational memory, and is reduction to just one degree of freedom to represent the reaction mechanism justified? It turns out that these questions can be formulated as experimentally testable mathematical inequalities, and their application to experimental and simulated data has already led to a number of insights. I will also discuss open issues and current challenges in this fast evolving field of research.

Broad distributions of transition-path times are fingerprints of multidimensionality of the underlying free energy landscapes

Recent single-molecule experiments have observed transition paths, i.e., brief events where molecules (particularly biomolecules) are caught in the act of surmounting activation barriers. Such measurements offer unprecedented mechanistic insights into the dynamics of biomolecular folding and binding, molecular machines, and biological membrane channels. A key challenge to these studies is to infer the complex details of the multidimensional energy landscape traversed by the transition paths from inherently low-dimensional experimental signals. A common minimalist model attempting to do so is that of one-dimensional diffusion along a reaction coordinate, yet its validity has been called into question. Here, we show that the distribution of the transition path time, which is a common experimental observable, can be used to differentiate between the dynamics described by models of one-dimensional diffusion from the dynamics in which multidimensionality is essential. Specifically, we prove that the coefficient of variation obtained from this distribution cannot possibly exceed 1 for any one-dimensional diffusive model, no matter how rugged its underlying free energy landscape is: In other words, this distribution cannot be broader than the single-exponential one. Thus, a coefficient of variation exceeding 1 is a fingerprint of multidimensional dynamics. Analysis of transition paths in atomistic simulations of proteins shows that this coefficient often exceeds 1, signifying essential multidimensionality of those systems.