Protein mechanical unfolding: A model with binary variables

The Wako-Saitô-Muñoz-Eaton Model for Predicting Protein Folding and Dynamics

Despite the recent advances in the prediction of protein structures by deep neutral networks, the elucidation of protein-folding mechanisms remains challenging. A promising theory for describing protein folding is a coarse-grained statistical mechanical model called the Wako-Saitô-Muñoz-Eaton (WSME) model. The model can calculate the free-energy landscapes of proteins based on a three-dimensional structure with low computational complexity, thereby providing a comprehensive understanding of the folding pathways and the structure and stability of the intermediates and transition states involved in the folding reaction. In this review, we summarize previous and recent studies on protein folding and dynamics performed using the WSME model and discuss future challenges and prospects. The WSME model successfully predicted the folding mechanisms of small single-domain proteins and the effects of amino-acid substitutions on protein stability and folding in a manner that was consistent with experimental results. Furthermore, extended versions of the WSME model were applied to predict the folding mechanisms of multi-domain proteins and the conformational changes associated with protein function. Thus, the WSME model may contribute significantly to solving the protein-folding problem and is expected to be useful for predicting protein folding, stability, and dynamics in basic research and in industrial and medical applications.

Emergence of effective temperatures in an out-of-equilibrium model of biopolymer folding

We investigate the possibility of extending the notion of temperature in a stochastic model for the RNA or protein folding driven out of equilibrium. We simulate the dynamics of a small RNA hairpin subject to an external pulling force, which is time-dependent. First, we consider a fluctuation-dissipation relation (FDR) whereby we verify that various effective temperatures can be obtained for different observables, only when the slowest intrinsic relaxation timescale of the system regulates the dynamics of the system. Then, we introduce a different nonequilibrium temperature, which is defined from the rate of heat exchanged with a weakly interacting thermal bath. Notably, this "kinetic" temperature can be defined for any frequency of the external switching force. We also discuss and compare the behavior of these two emerging parameters, by discriminating the time-delayed nature of the FDR temperature from the instantaneous character of the kinetic temperature. The validity of our numerics are corroborated by a simple four-state Markov model which describes the long-time behavior of the RNA molecule.

Mapping the Topography of a Protein Energy Landscape

Protein energy landscapes are highly complex, yet the vast majority of states within them tend to be invisible to experimentalists. Here, using site-directed mutagenesis and exploiting the simplicity of tandem-repeat protein structures, we delineate a network of these states and the routes between them. We show that our target, gankyrin, a 226-residue 7-ankyrin-repeat protein, can access two alternative (un)folding pathways. We resolve intermediates as well as transition states, constituting a comprehensive series of snapshots that map early and late stages of the two pathways and show both to be polarized such that the repeat array progressively unravels from one end of the molecule or the other. Strikingly, we find that the protein folds via one pathway but unfolds via a different one. The origins of this behavior can be rationalized using the numerical results of a simple statistical mechanics model that allows us to visualize the equilibrium behavior as well as single-molecule folding/unfolding trajectories, thereby filling in the gaps that are not accessible to direct experimental observation. Our study highlights the complexity of repeat-protein folding arising from their symmetrical structures; at the same time, however, this structural simplicity enables us to dissect the complexity and thereby map the precise topography of the energy landscape in full breadth and remarkable detail. That we can recapitulate the key features of the folding mechanism by computational analysis of the native structure alone will help toward the ultimate goal of designed amino-acid sequences with made-to-measure folding mechanisms-the Holy Grail of protein folding.

Worm-like Ising model for protein mechanical unfolding under the effect of osmolytes

We show via single-molecule mechanical unfolding experiments that the osmolyte glycerol stabilizes the native state of the human cardiac I27 titin module against unfolding without shifting its unfolding transition state on the mechanical reaction coordinate. Taken together with similar findings on the immunoglobulin-binding domain of streptococcal protein G (GB1), these experimental results suggest that osmolytes act on proteins through a common mechanism that does not entail a shift of their unfolding transition state. We investigate the above common mechanism via an Ising-like model for protein mechanical unfolding that adds worm-like-chain behavior to a recent generalization of the Wako-Saitô-Muñoz-Eaton model with support for group-transfer free energies. The thermodynamics of the model are exactly solvable, while protein kinetics under mechanical tension can be simulated via Monte Carlo algorithms. Notably, our force-clamp and velocity-clamp simulations exhibit no shift in the position of the unfolding transition state of GB1 and I27 under the effect of various osmolytes. The excellent agreement between experiment and simulation strongly suggests that osmolytes do not assume a structural role at the mechanical unfolding transition state of proteins, acting instead by adjusting the solvent quality for the protein chain analyte.

Effects of confinement on thermal stability and folding kinetics in a simple Ising-like model

In a cellular environment, confinement and macromolecular crowding play an important role in thermal stability and folding kinetics of a protein. We have resorted to a generalized version of the Wako-Saitô-Muñoz-Eaton model for protein folding to study the behavior of six different protein structures confined between two walls. Changing the distance 2R between the walls, we found, in accordance with previous studies, two confinement regimes: starting from large R and decreasing R, confinement first enhances the stability of the folded state as long as this is compact and until a given value of R; then a further decrease of R leads to a decrease of folding temperature and folding rate. We found that in the low confinement regime, both unfolding temperatures and logarithm of folding rates scale as R(-γ) where γ values lie in between 1.42 and 2.35.

Direction-dependent mechanical unfolding and green fluorescent protein as a force sensor

An Ising-like model of proteins is used to investigate the mechanical unfolding of the green fluorescent protein along different directions. When the protein is pulled from its ends, we recover the major and minor unfolding pathways observed in experiments. Upon varying the pulling direction, we find the correct order of magnitude and ranking of the unfolding forces. Exploiting the direction dependence of the unfolding force at equilibrium, we propose a force sensor whose luminescence depends on the applied force.

Analysis of the equilibrium and kinetics of the ankyrin repeat protein myotrophin

We apply the Wako-Saito-Muñoz-Eaton model to the study of myotrophin, a small ankyrin repeat protein, whose folding equilibrium and kinetics have been recently characterized experimentally. The model, which is a native-centric with binary variables, provides a finer microscopic detail than the Ising model that has been recently applied to some different repeat proteins, while being still amenable for an exact solution. In partial agreement with the experiments, our results reveal a weakly three-state equilibrium and a two-state-like kinetics of the wild-type protein despite the presence of a nontrivial free-energy profile. These features appear to be related to a careful "design" of the free-energy landscape, so that mutations can alter this picture, stabilizing some intermediates and changing the position of the rate-limiting step. Also, the experimental findings of two alternative pathways, an N-terminal and a C-terminal one, are qualitatively confirmed, even if the variations in the rates upon the experimental mutations cannot be quantitatively reproduced. Interestingly, the folding and unfolding pathways appear to be different, even if closely related: a property that is not generally considered in the phenomenological interpretation of the experimental data.

Exact nonequilibrium work generating function for a small classical system

We obtain the exact nonequilibrium work generating function (NEWGF) for a small system consisting of a massive Brownian particle connected to internal and external springs. The external work is provided to the system for a finite-time interval. The Jarzynski equality, obtained in this case directly from the NEWGF, is shown to be valid for the present model, in an exact way regardless of the rate of external work.

Unfolding times for proteins in a force clamp

The escape process from the native valley for proteins subjected to a constant stretching force is examined using a model for a beta barrel. For a wide range of forces, the unfolding dynamics can be treated as one-dimensional diffusion, parametrized in terms of the end-to-end distance. In particular, the escape times can be evaluated as first passage times for a Brownian particle moving on the protein free-energy landscape, using the Smoluchowski equation. At strong forces, the unfolding process can be viewed as a diffusive drift away from the native state, while at weak forces thermal activation is the relevant mechanism. An escape-time analysis within this approach reveals a crossover from an exponential to an inverse Gaussian escape-time distribution upon passing from weak to strong forces. Moreover, a single expression valid at weak and strong forces can be devised both for the average unfolding time as well as for the corresponding variance. The analysis offers a possible explanation of recent experimental findings for the proteins ddFLN4 and ubiquitin.

Nearly symmetrical proteins: folding pathways and transition states

The folding pathways of the B domain of protein A have been the subject of many experimental and computational studies. Based on a statistical mechanical model, it has been suggested that the native state symmetry leads to multiple pathways, highly dependent on temperature and denaturant concentration. Experiments, however, have not confirmed this scenario. By considering four nearly symmetrical proteins, one of them being the above molecule, here we show that, if contact energies are properly taken into account, a different picture emerges from kinetic simulations of the above-mentioned model. This is characterized by a dominant folding pathway, which is consistent with the most recent experimental results. Given the simplicity of the model, we also report on a direct sampling of the transition state.

Equilibrium properties and force-driven unfolding pathways of RNA molecules

The mechanical unfolding of a simple RNA hairpin and of a 236-base portion of the Tetrahymena thermophila ribozyme is studied by means of an Ising-like model. Phase diagrams and free energy landscapes are computed exactly and suggest a simple two-state behavior for the hairpin and the presence of intermediate states for the ribozyme. Nonequilibrium simulations give the possible unfolding pathways for the ribozyme, and the dominant pathway corresponds to the experimentally observed one.

Changing the mechanical unfolding pathway of FnIII10 by tuning the pulling strength

We investigate the mechanical unfolding of the tenth type III domain from fibronectin (FnIII(10)) both at constant force and at constant pulling velocity, by all-atom Monte Carlo simulations. We observe both apparent two-state unfolding and several unfolding pathways involving one of three major, mutually exclusive intermediate states. All three major intermediates lack two of seven native beta-strands, and share a quite similar extension. The unfolding behavior is found to depend strongly on the pulling conditions. In particular, we observe large variations in the relative frequencies of occurrence for the intermediates. At low constant force or low constant velocity, all three major intermediates occur with a significant frequency. At high constant force or high constant velocity, one of them, with the N- and C-terminal beta-strands detached, dominates over the other two. Using the extended Jarzynski equality, we also estimate the equilibrium free-energy landscape, calculated as a function of chain extension. The application of a constant pulling force leads to a free-energy profile with three major local minima. Two of these correspond to the native and fully unfolded states, respectively, whereas the third one can be associated with the major unfolding intermediates.

Mechanical unfolding and refolding pathways of ubiquitin

Mechanical unfolding and refolding of ubiquitin are studied by Monte Carlo simulations of a Gō model with binary variables. The exponential dependence of the time constants on the force is verified, and folding and unfolding lengths are computed, with good agreement with experimental results. Furthermore, the model exhibits intermediate kinetic states, as observed in experiments. Unfolding and refolding pathways and intermediate states, obtained by tracing single secondary structure elements, are consistent with simulations of previous all-atom models and with the experimentally observed step sizes.