Transition path times reveal memory effects and anomalous diffusion in the dynamics of protein folding

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.

Spectral Map for Slow Collective Variables, Markovian Dynamics, and Transition State Ensembles

Understanding the behavior of complex molecular systems is a fundamental problem in physical chemistry. To describe the long-time dynamics of such systems, which is responsible for their most informative characteristics, we can identify a few slow collective variables (CVs) while treating the remaining fast variables as thermal noise. This enables us to simplify the dynamics and treat it as diffusion in a free-energy landscape spanned by slow CVs, effectively rendering the dynamics Markovian. Our recent statistical learning technique, spectral map [Rydzewski, J. J. Phys. Chem. Lett. 2023, 14(22), 5216-5220], explores this strategy to learn slow CVs by maximizing a spectral gap of a transition matrix. In this work, we introduce several advancements into our framework, using a high-dimensional reversible folding process of a protein as an example. We implement an algorithm for coarse-graining Markov transition matrices to partition the reduced space of slow CVs kinetically and use it to define a transition state ensemble. We show that slow CVs learned by spectral map closely approach the Markovian limit for an overdamped diffusion. We demonstrate that coordinate-dependent diffusion coefficients only slightly affect the constructed free-energy landscapes. Finally, we present how spectral maps can be used to quantify the importance of features and compare slow CVs with structural descriptors commonly used in protein folding. Overall, we demonstrate that a single slow CV learned by spectral map can be used as a physical reaction coordinate to capture essential characteristics of protein folding.

Data-Driven Learning of the Generalized Langevin Equation with State-Dependent Memory

We present a data-driven method to learn stochastic reduced models of complex systems that retain a state-dependent memory beyond the standard generalized Langevin equation with a homogeneous kernel. The constructed model naturally encodes the heterogeneous energy dissipation by jointly learning a set of state features and the non-Markovian coupling among the features. Numerical results demonstrate the limitation of the standard generalized Langevin equation and the essential role of the broadly overlooked state-dependency nature in predicting molecule kinetics related to conformation relaxation and transition.

Construction of Coarse-Grained Molecular Dynamics with Many-Body Non-Markovian Memory

We introduce a machine-learning-based coarse-grained molecular dynamics model that faithfully retains the many-body nature of the intermolecular dissipative interactions. Unlike the common empirical coarse-grained models, the present model is constructed based on the Mori-Zwanzig formalism and naturally inherits the heterogeneous state-dependent memory term rather than matching the mean-field metrics such as the velocity autocorrelation function. Numerical results show that preserving the many-body nature of the memory term is crucial for predicting the collective transport and diffusion processes, where empirical forms generally show limitations.

Coordinate-Dependent Drift-Diffusion Reveals the Kinetic Intermediate Traps of Top7-Based Proteins

The computer-designed Top7 served as a scaffold to produce immunoreactive proteins by grafting of the 2F5 HIV-1 antibody epitope (Top7-2F5) followed by biotinylation (Top7-2F5-biotin). The resulting nonimmunoglobulin affinity proteins were effective in inducing and detecting the HIV-1 antibody. However, the grafted Top7-2F5 design led to protein aggregation, as opposed to the soluble biotinylated Top7-2F5-biotin. The structure-based model predicted that the thermodynamic cooperativity of Top7 increases after grafting and biotin-labeling, reducing their intermediate state populations. In this work, the folding kinetic traps that might contribute to the aggregation propensity are investigated by the diffusion theory. Since the engineered proteins have similar sequence and structural homology, they served as protein models to study the kinetic intermediate traps that were uncovered by characterizing the position-dependent drift-velocity (v(Q)) and the diffusion (D(Q)) coefficients. These coordinate-dependent coefficients were taken into account to obtain the folding and transition path times over the free energy transition states containing the intermediate kinetic traps. This analysis may be useful to predict the aggregated kinetic traps of scaffold-epitope proteins that might compose novel diagnostic and therapeutic platforms.

Pair-Reaction Dynamics in Water: Competition of Memory, Potential Shape, and Inertial Effects

When described by a one-dimensional reaction coordinate, pair-reaction rates in a solvent depend, in addition to the potential barrier height and the friction coefficient, on the potential shape, the effective mass, and the friction relaxation spectrum, but a rate theory that accurately accounts for all of these effects does not exist. After a review of classical reaction-rate theories, we show how to extract all parameters of the generalized Langevin equation (GLE) and, in particular, the friction memory function from molecular dynamics (MD) simulations of two prototypical pair reactions in water, the dissociation of NaCl and of two methane molecules. The memory exhibits multiple time scales and, for NaCl, pronounced oscillatory components. Simulations of the GLE by Markovian embedding techniques accurately reproduce the pair-reaction kinetics from MD simulations without any fitting parameters, which confirms the accuracy of the approximative form of the GLE and of the parameter extraction techniques. By modification of the GLE parameters, we investigate the relative importance of memory, mass, and potential shape effects. Neglect of memory slows down NaCl and methane dissociation by roughly a factor of 2; neglect of mass accelerates reactions by a similar factor, and the harmonic approximation of the potential shape gives rise to slight acceleration. This partial error cancellation explains why Kramers' theory, which neglects memory effects and treats the potential shape in harmonic approximation, describes reaction rates better than more sophisticated theories. In essence, all three effects, friction memory, inertia, and the potential shape nonharmonicity, are important to quantitatively describe pair-reaction kinetics in water.

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.

Conformational dynamics modulating electron transfer

Diffusional dynamics of the donor-acceptor distance are responsible for the appearance of a new time scale of diffusion over the distance of electronic tunneling in electron-transfer reactions. The distance dynamics compete with the medium polarization dynamics in the dynamics-controlled electron-transfer kinetics. The pre-exponential factor of the electron-transfer rate constant switches, at the crossover distance, between a distance-independent, dynamics-controlled plateau and exponential distance decay. The crossover between two regimes is controlled by an effective relaxation time slowed down by a factor exponentially depending on the variance of the donor-acceptor displacement. Flexible donor-acceptor complexes must show a greater tendency for dynamics-controlled electron transfer. Energy chains based on electron transport are best designed by placing the redox cofactors near the crossover distance.

Role of ergodicity, aging, and Gaussianity in resolving the origins of biomolecule subdiffusion

The internal motions of biomolecules are essential to their function. Although biological macromolecules conventionally show subdiffusive dynamics, only recently has subdiffusion been associated with non-ergodicity. These findings have stimulated new questions in biophysics and statistical mechanics. Is non-ergodic subdiffusion a general strategy shared by biomolecules? What underlying mechanisms are responsible for it? Here, we performed extensive molecular dynamics (MD) simulations to characterize the internal dynamics of six different biomolecules, ranging from single or double-stranded DNA, a single domain protein (KRAS), two globular proteins (PGK and SHP2), to an intrinsically disordered protein (SNAP-25). We found that the subdiffusive behavior of these biomolecules falls into two classes. The internal motion of the first three cases is ergodic subdiffusion and can be interpreted by fractional Brownian motion (FBM), while the latter three cases involve non-ergodic subdiffusion and can be modeled by mixed origins of continuous-time random walk (CTRW) and FBM.

Analysis of the Passage Times for Unfolding/Folding of the Adenine Riboswitch Aptamer

The conformational transitions of the adenosine deaminase A-riboswitch aptamer both with and without ligand binding are investigated within the tenets of the generalized Langevin equation in a complex viscoelastic cellular environment. Steered molecular dynamics (SMD) simulations are performed to evaluate and compare the results of the first passage times (FPTs) with those obtained from the theory for the unfold and fold transitions of the aptamer. The results of the distribution of Kramers's FPT reveal that the unfold-fold transitions are faster and hence more probable as compared to the fold-unfold transitions of the riboswitch aptamer for both ligand-bound and -unbound states. The transition path time is lower than Kramers's FPT for the riboswitch aptamer as the transition path times for the unfold-fold transition of both without and with ligand binding are insensitive to the details of the exact mechanism of the transition events. However, Kramers's FPTs show varied distributions which correspond to different transition pathways, unlike the transition path times. The mean FPT increases with an increase in the complexity of the cellular environment. The results of Kramers's FPT, transition path time distribution, and mean FPT obtained from our calculations qualitatively match with those obtained from the SMD simulations. Analytically derived values of the mean transition path time show good quantitative agreement with those estimated from the single-molecule force spectroscopy experiments for higher barrier heights.

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.

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.

Effect of ligand binding on riboswitch folding: Theory and simulations

The effect of ligand binding on the conformational transitions of the add A-riboswitch in cellular environments is investigated theoretically within the framework of the generalized Langevin equation combined with steered molecular dynamics simulations. Results for the transition path time distribution provide an estimate of the transit times, which are difficult to determine experimentally. The time for the conformational transitions of the riboswitch aptamer is longer for the ligand bound state as compared to that of the unbound one. The transition path time of the riboswitch follows a counterintuitive trend as it decreases with an increase in the barrier height. The mean transition path time of either transitions of the riboswitch in the ligand bound/unbound state increases with an increase in the complexity of the surrounding environment due to the caging effect. The results of the probability density function, transition path time distribution, and mean transition path time obtained from the theory qualitatively agree with those obtained from the simulations and with earlier experimental and theoretical studies.

Non-Markovian modeling of protein folding

Protein-folding kinetics is often described as Markovian (i.e., memoryless) diffusion in a one-dimensional free energy landscape, governed by an instantaneous friction coefficient that is fitted to reproduce experimental or simulated folding times. For the α-helix forming polypeptide alanine9 and a specific reaction coordinate that consists of the summed native hydrogen-bond lengths, we demonstrate that the friction extracted from molecular dynamics simulations exhibits significant memory with a decay time that is in the nanosecond range and thus, of the same order as the folding and unfolding times. Our non-Markovian modeling not only reproduces the molecular dynamics simulations accurately but also demonstrates that memory friction effects lead to anomalous and drastically accelerated protein kinetics.

We extract the folding free energy landscape and the time-dependent friction function, the two ingredients of the generalized Langevin equation (GLE), from explicit-water molecular dynamics (MD) simulations of the α-helix forming polypeptide alanine9 for a one-dimensional reaction coordinate based on the sum of the native H-bond distances. Folding and unfolding times from numerical integration of the GLE agree accurately with MD results, which demonstrate the robustness of our GLE-based non-Markovian model. In contrast, Markovian models do not accurately describe the peptide kinetics and in particular, cannot reproduce the folding and unfolding kinetics simultaneously, even if a spatially dependent friction profile is used. Analysis of the GLE demonstrates that memory effects in the friction significantly speed up peptide folding and unfolding kinetics, as predicted by the Grote–Hynes theory, and are the cause of anomalous diffusion in configuration space. Our methods are applicable to any reaction coordinate and in principle, also to experimental trajectories from single-molecule experiments. Our results demonstrate that a consistent description of protein-folding dynamics must account for memory friction effects.

Effect of Memory and Inertial Contribution on Transition-Time Distributions: Theory and Simulations

Transition paths refer to the time taken by molecules to cross a barrier separating two molecular conformations. In this work, we study how memory, as well as inertial contribution in the dynamics along a reaction coordinate, can affect the distribution of the transition-path time. We use a simple model of dynamics governed by a generalized Langevin equation with a power-law memory along with the inertial term, which was neglected in previous studies, where memory effects were explored only in the overdamped limit. We derive an approximate expression for the transit-time distribution and discuss our results for the short- and long-time limits and also compare it with known results in the high friction (overdamped) limit as well as in the Markovian limit. We have developed a numerical algorithm to test our theoretical results against extensive numerical simulations.

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.

Piecewise All-Atom SMD Simulations Reveal Key Secondary Structures in Luciferase Unfolding Pathway

Although the folding of single-domain proteins is well characterized theoretically and experimentally, the folding of large multidomain proteins is less well known. Firefly luciferase, a 550 residue three-domain protein, has been commonly used as a substrate to study chaperone reactions and as a model system for the study of folding of long polypeptide chains, including related phenomena such as cotranslational folding. Despite being characterized by various experimental techniques, the atomic-level contributions of various secondary structures of luciferase to its fold's mechanical stability remain unknown. Here, we developed a piecewise approach for all-atom steered molecular dynamics simulations to examine specific secondary structures that resist mechanical unfolding while minimizing the amount of computational resources required by the large water box of standard all-atom steered molecular dynamics simulations. We validated the robustness of this approach with a small NI3C protein and used our approach to elucidate the specific secondary structures that provide the largest contributions to luciferase mechanostability. In doing so, we show that piecewise all-atom steered molecular dynamics simulations can provide novel atomic resolution details regarding mechanostability and can serve as a platform for novel mutagenesis studies as well as a point for comparison with high-resolution force spectroscopy experiments.

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.

Disordered proteins follow diverse transition paths as they fold and bind to a partner

Transition paths of macromolecular conformational changes such as protein folding are predicted to be heterogeneous. However, experimental characterization of the diversity of transition paths is extremely challenging because it requires measuring more than one distance during individual transitions. In this work, we used fast three-color single-molecule Förster resonance energy transfer spectroscopy to obtain the distribution of binding transition paths of a disordered protein. About half of the transitions follow a path involving strong non-native electrostatic interactions, resulting in a transition time of 300 to 800 microseconds. The remaining half follow more diverse paths characterized by weaker electrostatic interactions and more than 10 times shorter transition path times. The chain flexibility and non-native interactions make diverse binding pathways possible, allowing disordered proteins to bind faster than folded proteins.

Kinetics of Loop Closure in Disordered Proteins: Theory vs Simulations vs Experiments

We study intrachain dynamics of intrinsically disordered proteins, as manifested by the time scales of loop formation, using atomistic simulations, experiment-parametrized coarse-grained models, and one-dimensional theories assuming Markov or non-Markov dynamics along the reaction coordinate. Despite the generally non-Markov character of monomer dynamics in polymers, we find that the simplest model of one-dimensional diffusion along the reaction coordinate (equated to the distance between the loop-forming monomers) well captures the mean first passage times to loop closure measured in coarse-grained and atomistic simulations, which, in turn, agree with the experimental values. This justifies use of the one-dimensional diffusion model in interpretation of experimental data. At the same time, the transition path times for loop closure in longer polypeptide chains show significant non-Markov effects; at intermediate times, these effects are better captured by the generalized Langevin equation model. At long times, however, atomistic simulations predict long tails in the distributions of transition path times, which are at odds with both the one-dimensional diffusion model and the generalized Langevin equation model.

Inferring effective forces for Langevin dynamics using Gaussian processes

Effective forces derived from experimental or in silico molecular dynamics time traces are critical in developing reduced and computationally efficient descriptions of otherwise complex dynamical problems. This helps motivate why it is important to develop methods to efficiently learn effective forces from time series data. A number of methods already exist to do this when data are plentiful but otherwise fail for sparse datasets or datasets where some regions of phase space are undersampled. In addition, any method developed to learn effective forces from time series data should be minimally a priori committal as to the shape of the effective force profile, exploit every data point without reducing data quality through any form of binning or pre-processing, and provide full credible intervals (error bars) about the prediction for the entirety of the effective force curve. Here, we propose a generalization of the Gaussian process, a key tool in Bayesian nonparametric inference and machine learning, which meets all of the above criteria in learning effective forces for the first time.

Transition path times in asymmetric barriers

Biomolecular conformational transitions are usually modeled as barrier crossings in a free energy landscape. The transition paths connect two local free energy minima and transition path times (TPT) are the actual durations of the crossing events. The simplest model employed to analyze TPT and to fit empirical data is that of a stochastic particle crossing a parabolic barrier. Motivated by some disagreement between the value of the barrier height obtained from the TPT distributions as compared to the value obtained from kinetic and thermodynamic analyses, we investigate here TPT for barriers which deviate from the symmetric parabolic shape. We introduce a continuous set of potentials, that starting from a parabolic shape, can be made increasingly asymmetric by tuning a single parameter. The TPT distributions obtained in the asymmetric case are very well-fitted by distributions generated by parabolic barriers. The fits, however, provide values for the barrier heights and diffusion coefficients which deviate from the original input values. We show how these findings can be understood from the analysis of the eigenvalues spectrum of the Fokker-Planck equation and highlight connections with experimental results.

Extension-Dependent Drift Velocity and Diffusion (DrDiff) Directly Reconstructs the Folding Free Energy Landscape of Atomic Force Microscopy Experiments

Two equilibrium force microscopy trajectories [q(t)] of high-precision single-molecule spectroscopy assays were analyzed: the pulling of an HIV RNA hairpin and of a 3-aa sequence of the bacteriorhodopsin membrane protein. Both present hundreds of two-state folding transitions, and their free-energy [F(q)] landscapes were previously obtained by deconvolving time signals with the inverse Boltzmann and p_fold methods. In this letter, the two F profiles were reconstructed directly from the measured time-series by the drift-diffusion (DrDiff) framework that characterized the effective conformational drift-velocity [v(q)] and diffusion [D(q)] coefficients. The two thermodynamic F profiles reconstructed with DrDiff directly from q(t) were in good agreement with those previously obtained from the deconvolved time signals. q(t) trajectories simulated with a two-dimensional framework in which the diffusion coefficient of the pulling setup (q coordinate) differed from the molecule (x coordinate) were also analyzed by DrDiff. The performance in reconstructing F was investigated in different conditions of diffusion anisotropy in the simulated time-series using Brownian dynamics. In addition, recently developed theories were used in order to evaluate the quality of the analysis performed in the experimental time series: the memory effects and the intrinsic biomolecular dynamic properties after connecting the probe to the molecule. With the 2-dimensional diffusive models and the additional analyses, it is proposed that the different physical regimes imposed by the stiffer probes of the two biomolecules will have an impact in the measured extension-dependent D and, thus, in the reconstruction of F by DrDiff. Stiffer AFM probes may reflect the molecular behavior more faithfully and reconstruction of F might be more successful. The reported quantities extracted directly from q(t) highlights the current state of the biomolecule characterization by force spectroscopy experiments: it is still challenging despite the recent advances, yet it is very promising.

Transition path dynamics in the binding of intrinsically disordered proteins: A simulation study

Association of proteins and other biopolymers is a ubiquitous process in living systems. Recent single-molecule measurements probe the dynamics of association in unprecedented detail by measuring the properties of association transition paths, i.e., short segments of molecular trajectories between the time the proteins are close enough to interact and the formation of the final complex. Interpretation of such measurements requires adequate models for describing the dynamics of experimental observables. In an effort to develop such models, here we report a simulation study of the association dynamics of two oppositely charged, disordered polymers. We mimic experimental measurements by monitoring intermonomer distances, which we treat as "experimental reaction coordinates." While the dynamics of the distance between the centers of mass of the molecules is found to be memoryless and diffusive, the dynamics of the experimental reaction coordinates displays significant memory and can be described by a generalized Langevin equation with a memory kernel. We compute the most commonly measured property of transition paths, the distribution of the transition path time, and show that, despite the non-Markovianity of the underlying dynamics, it is well approximated as one-dimensional diffusion in the potential of mean force provided that an apparent value of the diffusion coefficient is used. This apparent value is intermediate between the slow (low frequency) and fast (high frequency) limits of the memory kernel. We have further studied how the mean transition path time depends on the ionic strength and found only weak dependence despite strong electrostatic attraction between the polymers.

Understanding biochemical processes in the presence of sub-diffusive behavior of biomolecules in solution and living cells

In order to maintain cellular function, biomolecules like protein, DNA, and RNAs have to diffuse to the target spaces within the cell. Changes in the cytosolic microenvironment or in the nucleus during the fulfillment of these cellular processes affect their mobility, folding, and stability thereby impacting the transient or stable interactions with their adjacent neighbors in the organized and dynamic cellular interior. Using classical Brownian motion to elucidate the diffusion behavior of these biomolecules is hard considering their complex nature. The understanding of biomolecular diffusion inside cells still remains elusive due to the lack of a proper model that can be extrapolated to these cases. In this review, we have comprehensively addressed the progresses in this field, laying emphasis on the different aspects of anomalous diffusion in the different biochemical reactions in cell interior. These experiment-based models help to explain the diffusion behavior of biomolecules in the cytosolic and nuclear microenvironment. Moreover, since understanding of biochemical reactions within living cellular system is our main focus, we coupled the experimental observations with the concept of sub-diffusion from in vitro to in vivo condition. We believe that the pairing between the understanding of complex behavior and structure-function paradigm of biological molecules would take us forward by one step in order to solve the puzzle around diseases caused by cellular dysfunction.

Memory effects in single-molecule force spectroscopy measurements of biomolecular folding

Folding is generally assumed to be a Markov process, without memory. When the molecular motion is coupled to that of a probe as in single-molecule force spectroscopy (SMFS) experiments, however, theory predicts that the coupling to a second Markov process should induce memory when monitoring a projection of the full multi-dimensional motion onto a reduced coordinate. We developed a method to evaluate the time constant of the induced memory from its effects on the autocorrelation function, which can be readily determined from experimental data. Applying this method to both simulated SMFS measurements and experimental trajectories of DNA hairpin folding measured by optical tweezers as a model system, we validated the prediction that the linker induces memory. For these measurements, the timescale of the induced memory was found to be similar to the time required for the force probe to respond to changes in the molecule, and in the regime where the experimentally observed dynamics were not significantly perturbed by probe-molecule coupling artifacts. Memory effects are thus a general feature of SMFS measurements induced by the mechanical connection between the molecule and force probe that should be considered when interpreting experimental data.

Probing microscopic conformational dynamics in folding reactions by measuring transition paths

Transition paths comprise those parts of a folding trajectory where the molecule passes through the high-energy transition states separating folded and unfolded conformations. The transition states determine the folding kinetics and mechanism but are difficult to observe because of their brief duration. Single-molecule experiments have in recent years begun to characterize transition paths in folding reactions, allowing the microscopic conformational dynamics that occur as a molecule traverses the energy barriers to be probed directly. Here we review single-molecule fluorescence and force spectroscopy measurements of transition-path properties, including the time taken to traverse the paths, the local velocity along them, the path shapes, and the variability within these measurements reflecting differences between individual barrier crossings. We discuss how these measurements have been related to theories of folding as diffusion over an energy landscape to deduce properties such as the diffusion coefficient, and how they are being combined with simulations to obtain enhanced atomistic understanding of folding. The richly detailed information available from transition path measurements holds great promise for improved understanding of microscopic mechanisms in folding.

Exact Solutions for Distributions of First-Passage, Direct-Transit, and Looping Times in Symmetric Cusp Potential Barriers and Wells

For a particle diffusing in one dimension, the distribution of its first-passage time from point a to point b is determined by the durations of the particle trajectories that start from point a and are terminated as soon as they touch point b for the first time. Any such trajectory consists of looping and direct-transit segments. The latter is the final part of the trajectory that leaves point a and goes to point b without returning to point a. The rest of the trajectory is the looping segment that makes numerous loops which begin and end at the same point a without touching point b. In this article we discuss general relations between the first-passage time distribution and those for the durations of the two segments. These general relations allow us to find exact solutions for the Laplace transforms of the distributions of the first-passage, direct-transit, and looping times for transitions between two points separated by a symmetric cusp potential barrier or well of arbitrary height and depth, respectively. The obtained Laplace transforms are inverted numerically, leading to nontrivial dependences of the resulting distributions on the barrier height and the well depth.

Generalized Langevin Equation as a Model for Barrier Crossing Dynamics in Biomolecular Folding

Conformational memory in single-molecule dynamics has attracted recent attention and, in particular, has been invoked as a possible explanation of some of the intriguing properties of transition paths observed in single-molecule force spectroscopy (SMFS) studies. Here we study one candidate for a non-Markovian model that can account for conformational memory, the generalized Langevin equation with a friction force that depends not only on the instantaneous velocity but also on the velocities in the past. The memory in this model is determined by a time-dependent friction memory kernel. We propose a method for extracting this kernel directly from an experimental signal and illustrate its feasibility by applying it to a generalized Rouse model of a SMFS experiment, where the memory kernel is known exactly. Using the same model, we further study how memory affects various statistical properties of transition paths observed in SMFS experiments and evaluate the performance of recent approximate analytical theories of non-Markovian dynamics of barrier crossing. We argue that the same type of analysis can be applied to recent single-molecule observations of transition paths in protein and DNA folding.

Dynamics of Disordered Proteins under Confinement: Memory Effects and Internal Friction

Many proteins are disordered under physiological conditions. How efficiently they can search for their cellular targets and how fast they can fold upon target binding is determined by their intrinsic dynamics, which have thus attracted much recent attention. Experiments and molecular simulations show that the inherent reconfiguration timescale for unfolded proteins has a solvent friction component and an internal friction component, and the microscopic origin of the latter, along with its proper mathematical description, has been a topic of considerable debate. Internal friction varies across different proteins of comparable length and increases with decreasing denaturant concentration, showing that it depends on how compact the protein is. Here we report on a systematic atomistic simulation study of how confinement, which induces a more compact unfolded state, affects dynamics and friction in disordered peptides. We find that the average reconfiguration timescales increase exponentially as the peptide's spatial dimensions are reduced; at the same time, confinement broadens the spectrum of relaxation timescales exhibited by the peptide. There are two important implications of this broadening: First, it limits applicability of the common Rouse and Zimm models with internal friction, as those models attempt to capture internal friction effects by introducing a single internal friction timescale. Second, the long-tailed distribution of relaxation times leads to anomalous diffusion effects in the dynamics of intramolecular distances. Analysis and interpretation of experimental signals from various measurements that probe intramolecular protein dynamics (such as single-molecule fluorescence correlation spectroscopy and single-molecule force spectroscopy) rely on the assumption of diffusive dynamics along the distances being probed; hence, our results suggest the need for more general models allowing for anomalous diffusion effects.

Accurate Protein-Folding Transition-Path Statistics from a Simple Free-Energy Landscape

A central goal of protein-folding theory is to predict the stochastic dynamics of transition paths-the rare trajectories that transit between the folded and unfolded ensembles-using only thermodynamic information, such as a low-dimensional equilibrium free-energy landscape. However, commonly used one-dimensional landscapes typically fall short of this aim, because an empirical coordinate-dependent diffusion coefficient has to be fit to transition-path trajectory data in order to reproduce the transition-path dynamics. We show that an alternative, first-principles free-energy landscape predicts transition-path statistics that agree well with simulations and single-molecule experiments without requiring dynamical data as an input. This "topological configuration" model assumes that distinct, native-like substructures assemble on a time scale that is slower than native-contact formation but faster than the folding of the entire protein. Using only equilibrium simulation data to determine the free energies of these coarse-grained intermediate states, we predict a broad distribution of transition-path transit times that agrees well with the transition-path durations observed in simulations. We further show that both the distribution of finite-time displacements on a one-dimensional order parameter and the ensemble of transition-path trajectories generated by the model are consistent with the simulated transition paths. These results indicate that a landscape based on transient folding intermediates, which are often hidden by one-dimensional projections, can form the basis of a predictive model of protein-folding transition-path dynamics.

Communication: Coordinate-dependent diffusivity from single molecule trajectories

Single-molecule observations of biomolecular folding are commonly interpreted using the model of one-dimensional diffusion along a reaction coordinate, with a coordinate-independent diffusion coefficient. Recent analysis, however, suggests that more general models are required to account for single-molecule measurements performed with high temporal resolution. Here, we consider one such generalization: a model where the diffusion coefficient can be an arbitrary function of the reaction coordinate. Assuming Brownian dynamics along this coordinate, we derive an exact expression for the coordinate-dependent diffusivity in terms of the splitting probability within an arbitrarily chosen interval and the mean transition path time between the interval boundaries. This formula can be used to estimate the effective diffusion coefficient along a reaction coordinate directly from single-molecule trajectories.

Single-Molecule Test for Markovianity of the Dynamics along a Reaction Coordinate

In an effort to answer the much-debated question of whether the time evolution of common experimental observables can be described as one-dimensional diffusion in the potential of mean force, we propose a simple criterion that allows one to test whether the Markov assumption is applicable to a single-molecule trajectory x( t). This test does not involve fitting of the data to any presupposed model and can be applied to experimental data with relatively low temporal resolution.

Slow domain reconfiguration causes power-law kinetics in a two-state enzyme

Protein dynamics are typically captured well by rate equations that predict exponential decays for two-state reactions. Here, we describe a remarkable exception. The electron-transfer enzyme quiescin sulfhydryl oxidase (QSOX), a natural fusion of two functionally distinct domains, switches between open- and closed-domain arrangements with apparent power-law kinetics. Using single-molecule FRET experiments on time scales from nanoseconds to milliseconds, we show that the unusual open-close kinetics results from slow sampling of an ensemble of disordered domain orientations. While substrate accelerates the kinetics, thus suggesting a substrate-induced switch to an alternative free energy landscape of the enzyme, the power-law behavior is also preserved upon electron load. Our results show that the slow sampling of open conformers is caused by a variety of interdomain interactions that imply a rugged free energy landscape, thus providing a generic mechanism for dynamic disorder in multidomain enzymes.

The influence of absorbing boundary conditions on the transition path time statistics

A novel analytical expression, explicitly taking into account absorbing boundaries, exactly describes TPT distributions for particles crossing a parabolic potential.

Protein folding transition path times from single molecule FRET

The transition path is the tiny segment of a single molecule trajectory when the free energy barrier between states is crossed and for protein folding contains all of the information about the self-assembly mechanism. As a first step toward obtaining structural information during the transition path from experiments, single molecule FRET spectroscopy has been used to determine average transition path times from a photon-by-photon analysis of fluorescence trajectories. These results, obtained for several different proteins, have already provided new and demanding tests that support both the accuracy of all-atom molecular dynamics simulations and the basic postulates of energy landscape theory of protein folding.