Proteins Fold by Subdiffusion of the Order Parameter

Emergence of Memory in Equilibrium versus Nonequilibrium Systems

Experiments often probe observables that correspond to low-dimensional projections of high-dimensional dynamics. In such situations distinct microscopic configurations become lumped into the same observable state. It is well known that correlations between the observable and the hidden degrees of freedom give rise to memory effects. However, how and under which conditions these correlations emerge remain poorly understood. Here we shed light on two fundamentally different scenarios of the emergence of memory in minimal stationary systems, where observed and hidden degrees of freedom either evolve cooperatively or are coupled by a hidden nonequilibrium current. In the reversible setting the strongest memory manifests when the timescales of hidden and observed dynamics overlap, whereas, strikingly, in the driven setting maximal memory emerges under a clear timescale separation. Our results hint at the possibility of fundamental differences in the way memory emerges in equilibrium versus driven systems that may be utilized as a "diagnostic" of the underlying hidden transport mechanism.

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.

Origin of subdiffusions in proteins: Insight from peptide systems

Subdiffusive kinetics are popular in proteins and peptides as observed in experiments and simulations. For protein systems with diverse interactions, are there multiple mechanisms to produce the common subdiffusion behavior? To approach this problem, long trajectories of two model peptides are simulated to study the mechanism of subdiffusion and the relations with their interactions. The free-energy profiles and the subdiffusive kinetics are observed for these two peptides. A hierarchical plateau analysis is employed to extract the features of the landscape from the mean square of displacement. The mechanism of subdiffusions can be postulated by comparing the exponents by simulations with those based on various models. The results indicate that the mechanisms of these two peptides are different and are related to the characteristics of their energy landscapes. The subdiffusion of the flexible peptide is mainly caused by depth distribution of traps on the energy landscape, while the subdiffusion of the helical peptide is attributed to the fractal topology of local minima on the landscape. The emergence of these different mechanisms reflects different kinetic scenarios in peptide systems though the peptides behave in a similar way of diffusion. To confirm these ideas, the transition networks between various conformations of these peptides are generated. Based on the network description, the controlled kinetics based only on the topology of the networks are calculated and compared with the results based on simulations. For the flexible peptide, the feature of controlled diffusion is distinct from that of simulation, and for the helical peptide, two kinds of kinetics have a similar exponent of subdiffusion. These results further exemplify the importance of the landscape topology in the kinetics of structural proteins and the effect of depth distribution of traps for the subdiffusion of disordered peptides.

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.

Stochastic diffusion framework determines the free-energy landscape and rate from single-molecule trajectory

A theoretical stochastic diffusion framework is developed that characterizes the position-dependent diffusion coefficient [D(Q)] and drift velocity [ v (Q)] by analysing single-molecule time traces [Q(t)]. The free-energy landscape [F(Q)] that governs the dynamics is reconstructed with the calculated D and v . There are many computational tools that perform this task in which some are computationaly demanding, difficult to run, and, most of the time, not directly available to the community. This is a first attempt to implement the simplified stochastic diffusion framework that is fast, easy to run in a Python environment, and available to be extended as needed. It does not require adjustable parameters, inference methods, or sampling bias such as Monte Carlo Bayesian estimators or umbrella samplings. The stochastic framework was applied in the protein-like lattice model with Monte Carlo simulations, which accurately predicted the folding rates with the coordinate-dependent D and F plugged into Kramers' theory. The results were compared with two other independently developed methodologies (the Bayesian analysis and fep1D algorithm) presenting a good match, which confirms its validity. This theoretical framework might be useful in determining the free-energy and rates by providing time series only from biological or condensed-phase systems. The code is freely available at https://github.com/ronaldolab/stochastic_diffusion.

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.

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

Recent single-molecule experiments probed transition paths of biomolecular folding and, in particular, measured the time biomolecules spend while crossing their free energy barriers. A surprising finding from these studies is that the transition barriers crossed by transition paths, as inferred from experimentally observed transition path times, are often lower than the independently determined free energy barriers. Here we explore memory effects leading to anomalous diffusion as a possible origin of this discrepancy. Our analysis of several molecular dynamics trajectories shows that the dynamics of common reaction coordinates used to describe protein folding is subdiffusive, at least at sufficiently short times. We capture this effect using a one-dimensional fractional Brownian motion (FBM) model, in which the system undergoes a subdiffusive process in the presence of a potential of mean force, and show that this model yields much broader distributions of transition path times with stretched exponential long-time tails. Without any adjustable parameters, these distributions agree well with the transition path times computed directly from protein trajectories. We further discuss how the FBM model can be tested experimentally.

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.

Quantifying the Sources of Kinetic Frustration in Folding Simulations of Small Proteins

Experiments and atomistic simulations of polypeptides have revealed structural intermediates that promote or inhibit conformational transitions to the native state during folding. We invoke a concept of “kinetic frustration” to quantify the prevalence and impact of these behaviors on folding rates within a large set of atomistic simulation data for 10 fast-folding proteins, where each protein’s conformational space is represented as a Markov state model of conformational transitions. Our graph theoretic approach addresses what conformational features correlate with folding inhibition and therefore permits comparison among features within a single protein network and also more generally between proteins. Nonnative contacts and nonnative secondary structure formation can thus be quantitatively implicated in inhibiting folding for several of the tested peptides.

Exploring multi-dimensional coordinate-dependent diffusion dynamics on the energy landscape of protein conformation change

We explore the multi-dimensional diffusion dynamics of protein conformational change. We found in general that the diffusion is anisotropic and inhomogeneous. The directional and positional dependence of diffusion have significant impacts on the protein conformational kinetics: the dominant kinetic path of conformational change is shifted from the naively expected steepest decent gradient paths. The kinetic transition state is shifted away from the transition state. The effective kinetic free energy barrier height, determining the kinetic rate of the conformational change, is shifted away from the one estimated from the thermodynamic free energy barrier. The shift of the transition state in position and value will modify the phi value analysis for identification of hot residues and interactions responsible for conformational dynamics. Ongoing and future experiments can test the predictions of the model.

Configuration-dependent diffusion dynamics of downhill and two-state protein folding

Configuration-dependent diffusion (CDD) is important for protein folding kinetics with small thermodynamic barriers. CDD can be even more crucial in downhill folding without thermodynamic barriers. We explored the CDD of a downhill protein (BBL), and a two-state protein (CI2). The hidden kinetic barriers due to CDD were revealed. The increased ~1 k(B)T kinetic barrier is in line with experimental value based on other fast folding proteins. Compared to that of CI2, the effective free-energy profile of BBL is found to be significantly influenced by CDD, and the kinetics are totally determined by diffusion. These findings are consistent with both earlier bulk and single-molecule fluorescence measurements. In addition, we found the temperature dependence of CDD. We also found that the ratio of folding transition temperature against optimal kinetic folding temperature can provide both a quantitative measure for the underlying landscape topography and an indicator for the possible appearance of downhill folding. Our study can help for a better understanding of the role of diffusion in protein folding dynamics.

Diffusion models of protein folding

In theory and in the analysis of experiments, protein folding is often described as diffusion along a single coordinate. We explore here the application of a one-dimensional diffusion model to interpret simulations of protein folding, where the parameters of a model that "best" describes the simulation trajectories are determined using a Bayesian analysis. We discuss the requirements for such a model to be a good approximation to the global dynamics, and several methods for testing its accuracy. For example, one test considers the effect of an added bias potential on the fitted free energies and diffusion coefficients. Such a bias may also be used to extend our approach to determining parameters for the model to systems that would not normally explore the full coordinate range on accessible time scales. Alternatively, the propagators predicted from the model at different "lag" times may be compared with observations from simulation. We then present some applications of the model to protein folding, including Kramers-like turnover in folding rates of coarse-grained models, the effect of non-native interactions on folding, and the effect of the chosen coordinate on the observed position-dependence of the diffusion coefficients. Lastly, we consider how our results are useful for the interpretation of experiments, and how this type of Bayesian analysis may eventually be applied directly to analyse experimental data.

Fractal symmetry of protein interior: what have we learned?

The application of fractal dimension-based constructs to probe the protein interior dates back to the development of the concept of fractal dimension itself. Numerous approaches have been tried and tested over a course of (almost) 30 years with the aim of elucidating the various facets of symmetry of self-similarity prevalent in the protein interior. In the last 5 years especially, there has been a startling upsurge of research that innovatively stretches the limits of fractal-based studies to present an array of unexpected results on the biophysical properties of protein interior. In this article, we introduce readers to the fundamentals of fractals, reviewing the commonality (and the lack of it) between these approaches before exploring the patterns in the results that they produced. Clustering the approaches in major schools of protein self-similarity studies, we describe the evolution of fractal dimension-based methodologies. The genealogy of approaches (and results) presented here portrays a clear picture of the contemporary state of fractal-based studies in the context of the protein interior. To underline the utility of fractal dimension-based measures further, we have performed a correlation dimension analysis on all of the available non-redundant protein structures, both at the level of an individual protein and at the level of structural domains. In this investigation, we were able to separately quantify the self-similar symmetries in spatial correlation patterns amongst peptide-dipole units, charged amino acids, residues with the π-electron cloud and hydrophobic amino acids. The results revealed that electrostatic environments in the interiors of proteins belonging to 'α/α toroid' (all-α class) and 'PLP-dependent transferase-like' domains (α/β class) are highly conducive. In contrast, the interiors of 'zinc finger design' ('designed proteins') and 'knottins' ('small proteins') were identified as folds with the least conducive electrostatic environments. The fold 'conotoxins' (peptides) could be unambiguously identified as one type with the least stability. The same analyses revealed that peptide-dipoles in the α/β class of proteins, in general, are more correlated to each other than are the peptide-dipoles in proteins belonging to the all-α class. Highly favorable electrostatic milieu in the interiors of TIM-barrel, α/β-hydrolase structures could explain their remarkably conserved (evolutionary) stability from a new light. Finally, we point out certain inherent limitations of fractal constructs before attempting to identify the areas and problems where the implementation of fractal dimension-based constructs can be of paramount help to unearth latent information on protein structural properties.

Turbulent phenomena in protein folding

Protein folding and hydrodynamic turbulence are two long-standing challenges, in molecular biophysics and fluid dynamics, respectively. The theories of these phenomena have been developed independently and used different formalisms. Here we show that the protein folding flows can be surprisingly similar to turbulent fluid flows. Studying a benchmark model protein (an SH3 domain), we have found that the flows for the slow folding trajectories of the protein, in which a partly formed N- and C-terminal β sheet hinders the RT loop from attaching to the protein core, have many properties of turbulent flows of a fluid. The flows are analyzed in a three-dimensional (3D) space of collective variables, which are the numbers of native contacts between the terminal β strands, between the RT loop and the protein core, and the rest of the native contacts. We have found that the flows have fractal nature and are filled with 3D eddies; the latter contain strange attractors, at which the tracer flow paths behave as saddle trajectories. Two regions of the space increment have been observed, in which the flux variations are self-similar with the scaling exponent h=1/3, in surprising agreement with the Kolmogorov inertial range theory of turbulence. In one region, the cascade of protein rearrangements is directed from larger to smaller scales (net folding), and in the other, it is oppositely directed (net unfolding). Folding flows for the fast trajectories are essentially "laminar" and do not have the property of self-similarity. Based on the results of our study, we infer, and support this inference by simulations, that the origin of the similarity between the protein folding and turbulent motion of a fluid is in a cascade mechanism of structural transformations in the systems that underlies these phenomena.

Is protein folding sub-diffusive?

Protein folding dynamics is often described as diffusion on a free energy surface considered as a function of one or few reaction coordinates. However, a growing number of experiments and models show that, when projected onto a reaction coordinate, protein dynamics is sub-diffusive. This raises the question as to whether the conventionally used diffusive description of the dynamics is adequate. Here, we numerically construct the optimum reaction coordinate for a long equilibrium folding trajectory of a Go model of a -repressor protein. The trajectory projected onto this coordinate exhibits diffusive dynamics, while the dynamics of the same trajectory projected onto a sub-optimal reaction coordinate is sub-diffusive. We show that the higher the (cut-based) free energy profile for the putative reaction coordinate, the more diffusive the dynamics become when projected on this coordinate. The results suggest that whether the projected dynamics is diffusive or sub-diffusive depends on the chosen reaction coordinate. Protein folding can be described as diffusion on the free energy surface as function of the optimum reaction coordinate. And conversely, the conventional reaction coordinates, even though they might be based on physical intuition, are often sub-optimal and, hence, show sub-diffusive dynamics.

To understand dynamics of complex systems with many degrees of freedom, one often projects it onto one or several collective variables. Protein folding, the complex, concerted motion of a protein chain towards a unique three-dimensional structure, is one example of where such reduction of complexity is useful. It is usually assumed that the projected dynamics is diffusive. However, many experiments and simulations have shown that the projected dynamics is sub-diffusive, i.e., the mean square displacement grows slower than linear with time. It means that the dynamics has a memory; that the free energy surface together with diffusion coefficient do not properly define the dynamics; and that such projections cannot be used to accurately describe dynamics. Here, we show that if one carefully constructs the reaction coordinate by optimizing (maximizing) its free energy profile, one can use a simple (memory-less) diffusive description. Loosely speaking, when the complex dynamics is projected onto a simple coordinate, all the complexity of the original dynamics goes into the memory of the projected dynamics. If the dynamics is projected onto the (complex) optimum reaction coordinate, all the complexity of the original dynamics is in the reaction coordinate, and the projected dynamics is simple.

Coordinate and time-dependent diffusion dynamics in protein folding

We developed both analytical and simulation methods to explore the diffusion dynamics in protein folding. We found the diffusion as a quantitative measure of escape from local traps along the protein folding funnel with chosen reaction coordinates has two remarkable effects on kinetics. At a fixed coordinate, local escape time depends on the distribution of barriers around it, therefore the diffusion is often time distributed. On the other hand, the environments (local escape barriers) change along the coordinates, therefore diffusion is coordinate dependent. The effects of time-dependent diffusion on folding can lead to non-exponential kinetics and non-Poisson statistics of folding time distribution. The effects of coordinate dependent diffusion on folding can lead to the change of the kinetic barrier height as well as the position of the corresponding transition state and therefore modify the folding kinetic rates as well as the kinetic routes. Our analytical models for folding are based on a generalized Fokker-Planck diffusion equation with diffusion coefficient both dependent on coordinate and time. Our simulation for folding are based on structure-based folding models with a specific fast folding protein CspTm studied experimentally on diffusion and folding with single molecules. The coordinate and time-dependent diffusion are especially important to be considered in fast folding and single molecule studies, when there is a small or no free energy barrier and kinetics is controlled by diffusion while underlying statistics of kinetics become important. Including the coordinate dependence of diffusion will challenge the transition state theory of protein folding. The classical transition state theory will have to be modified to be consistent. The more detailed folding mechanistic studies involving phi value analysis based on the classical transition state theory will also have to be quantitatively modified. Complex kinetics with multiple time scales may allow us not only to explore the folding kinetics but also probe the local landscape and barrier height distribution with single-molecule experiments.