Mechanical properties of the domains of titin in a Go-like model

Connecting thermal and mechanical protein (un)folding landscapes

Molecular dynamics simulations supplement single-molecule pulling experiments by providing the possibility of examining the full free energy landscape using many coordinates. Here, we use an all-atom structure-based model to study the force and temperature dependence of the unfolding of the protein filamin by applying force at both termini. The unfolding time-force relation τ(F) indicates that the force-induced unfolding behavior of filamin can be characterized into three regimes: barrier-limited low- and intermediate-force regimes, and a barrierless high-force regime. Slope changes of τ(F) separate the three regimes. We show that the behavior of τ(F) can be understood from a two-dimensional free energy landscape projected onto the extension X and the fraction of native contacts Q. In the low-force regime, the unfolding rate is roughly force-independent due to the small (even negative) separation in X between the native ensemble and transition state ensemble (TSE). In the intermediate-force regime, force sufficiently separates the TSE from the native ensemble such that τ(F) roughly follows an exponential relation. This regime is typically explored by pulling experiments. While X may fail to resolve the TSE due to overlap with the unfolded ensemble just below the folding temperature, the overlap is minimal at lower temperatures where experiments are likely to be conducted. The TSE becomes increasingly structured with force, whereas the average order of structural events during unfolding remains roughly unchanged. The high-force regime is characterized by barrierless unfolding, and the unfolding time approaches a limit of ∼10 μs for the highest forces we studied. Finally, a combination of X and Q is shown to be a good reaction coordinate for almost the entire force range.

Cooperativity in thermal and force-induced protein unfolding: integration of crack propagation and network elasticity models

We investigate force-induced and temperature-induced unfolding of proteins using the combination of a gaussian network model and a crack propagation model based on "bond"-breaking independent events. We assume the existence of threshold values for the mean strain and strain fluctuations that dictate bond rupture. Surprisingly, we find that this stepwise process usually leads to a few cooperative, first-order-like, transitions in which several bonds break simultaneously, reminiscent of the "avalanches" seen in disordered networks.

Mechanical stability of multidomain proteins and novel mechanical clamps

We estimate the size of mechanostability for 318 multidomain proteins which are single-chain and contain up to 1021 amino acids. We predict existence of novel types of mechanical clamps in which interdomain contacts play an essential role. Mechanical clamps are structural regions which are the primary source of a protein's resistance to pulling. Among these clamps there is one that opposes tensile stress due to two domains swinging apart. This movement strains and then ruptures the contacts that hold the two domains together. Another clamp also involves tensile stress but it originates from an immobilization of a structural region by a surrounding knot-loop (without involving any disulfide bonds). Still another mechanism involves shear between helical regions belonging to two domains. We also consider the amyloid-prone cystatin C which provides an example of a two-chain 3D domain-swapped protein. We predict that this protein should withstand remarkably large stress, perhaps of order 800 pN, when inducing a shearing strain. The survey is generated through molecular dynamics simulations performed within a structure-based coarse grained model.

Energy landscape, structure and rate effects on strength properties of alpha-helical proteins

The strength of protein domains is crucial to identify the mechanical role of protein domains in biological processes such as mechanotransduction, tissue mechanics and tissue remodeling. Whereas the concept of strength has been widely investigated for engineered materials, the strength of fundamental protein material building blocks and how it depends on structural parameters such as the chemical bonding, the protein filament length and the timescale of observation or deformation velocity remains poorly understood. Here we report a systematic analysis of the influence of key parameters that define the energy landscape of the strength properties of alpha-helical protein domains, including energy barriers, unfolding and refolding distances, the locations of folded and unfolded states, as well as variations of the length and pulling velocity of alpha-helical protein filaments. The analysis is facilitated by the development of a double-well mesoscale potential formulation, utilized here to carry out a systematic numerical analysis of the behavior of alpha-helices. We compare the results against widely used protein strength models based on the Bell model, one of the simplest models used to characterize the strength of protein filaments. We find that, whereas Bell-type models are a reasonable approximation to describe the rupture of alpha-helical protein domains for a certain range of pulling speeds and values of energy barriers, the model ceases to hold for very large energy barriers and for very small pulling speeds, in agreement with earlier findings. We conclude with an application of our mesoscale model to investigate the effect of the length of alpha-helices on their mechanical strength. We find a weakening effect as the length of alpha-helical proteins increases, followed by an asymptotic regime in which the strength remains constant. We compare strand lengths found in biological proteins with the scaling law of strength versus alpha-helix filament length. The mesoscale model reported here is generally applicable to other protein filaments that feature a serial array of domains that unfold under applied strain, where a similar length-dependent strength could be observed.

Model for stretching and unfolding the giant multidomain muscle protein using single-molecule force spectroscopy

Single-molecule manipulation has allowed the forced unfolding of multidomain proteins. Here we outline a theory that not only explains these experiments but also points out a number of difficulties in their interpretation and makes suggestions for further experiments. For titin we reproduce force-extension curves, the dependence of break force on pulling speed, and break-force distributions and also validate two common experimental views: Unfolding titin Ig domains can be explained as stepwise increases in contour length, and increasing force peaks in native Ig sequences represent a hierarchy of bond strengths. Our theory is valid for essentially any molecule that can be unfolded in atomic force microscopy; as a further example, we present force-extension curves for the unfolding of RNA hairpins.

Hydrogen bond networks determine emergent mechanical and thermodynamic properties across a protein family

Background

Gram-negative bacteria use periplasmic-binding proteins (bPBP) to transport nutrients through the periplasm. Despite immense diversity within the recognized substrates, all members of the family share a common fold that includes two domains that are separated by a conserved hinge. The hinge allows the protein to cycle between open (apo) and closed (ligated) conformations. Conformational changes within the proteins depend on a complex interplay of mechanical and thermodynamic response, which is manifested as an increase in thermal stability and decrease of flexibility upon ligand binding.

Results

We use a distance constraint model (DCM) to quantify the give and take between thermodynamic stability and mechanical flexibility across the bPBP family. Quantitative stability/flexibility relationships (QSFR) are readily evaluated because the DCM links mechanical and thermodynamic properties. We have previously demonstrated that QSFR is moderately conserved across a mesophilic/thermophilic RNase H pair, whereas the observed variance indicated that different enthalpy-entropy mechanisms allow similar mechanical response at their respective melting temperatures. Our predictions of heat capacity and free energy show marked diversity across the bPBP family. While backbone flexibility metrics are mostly conserved, cooperativity correlation (long-range couplings) also demonstrate considerable amount of variation. Upon ligand removal, heat capacity, melting point, and mechanical rigidity are, as expected, lowered. Nevertheless, significant differences are found in molecular cooperativity correlations that can be explained by the detailed nature of the hydrogen bond network.

Conclusion

Non-trivial mechanical and thermodynamic variation across the family is explained by differences within the underlying H-bond networks. The mechanism is simple; variation within the H-bond networks result in altered mechanical linkage properties that directly affect intrinsic flexibility. Moreover, varying numbers of H-bonds and their strengths control the likelihood for energetic fluctuations as H-bonds break and reform, thus directly affecting thermodynamic properties. Consequently, these results demonstrate how unexpected large differences, especially within cooperativity correlation, emerge from subtle differences within the underlying H-bond network. This inference is consistent with well-known results that show allosteric response within a family generally varies significantly. Identifying the hydrogen bond network as a critical determining factor for these large variances may lead to new methods that can predict such effects.

Selection of optimal variants of Gō-like models of proteins through studies of stretching

The Gō-like models of proteins are constructed based on the knowledge of the native conformation. However, there are many possible choices of a Hamiltonian for which the ground state coincides with the native state. Here, we propose to use experimental data on protein stretching to determine what choices are most adequate physically. This criterion is motivated by the fact that stretching processes usually start with the native structure, in the vicinity of which the Gō-like models should work the best. Our selection procedure is applied to 62 different versions of the Gō model and is based on 28 proteins. We consider different potentials, contact maps, local stiffness energies, and energy scales--uniform and nonuniform. In the latter case, the strength of the nonuniformity was governed either by specificity or by properties related to positioning of the side groups. Among them is the simplest variant: uniform couplings with no i, i + 2 contacts. This choice also leads to good folding properties in most cases. We elucidate relationship between the local stiffness described by a potential which involves local chirality and the one which involves dihedral and bond angles. The latter stiffness improves folding but there is little difference between them when it comes to stretching.

An entropic perspective of protein stability on surfaces

The interaction of proteins with surfaces regulates numerous processes in nature, science, and technology. In many applications, it is desirable to place proteins on surfaces in an active state, and tethering represents one manner in which to accomplish this. However, a clear understanding of how tether placement and design affects protein activity is lacking. Available theoretical models predict that proteins will be stabilized when tethered to substrates. Such models suggest that the surface reduces the number of states accessible to the unfolded state of the protein, thereby reducing the entropic cost of folding on the surface compared to the bulk case. Recent studies, however, have shown that this stabilization is not always seen. The purpose of this article is to determine the validity of the theory with a thorough thermodynamic analysis of the folding of peptides attached to surfaces. Configuration-temperature-density-of-states Monte Carlo simulations are used to examine the behavior of four different peptides of different secondary and tertiary structure. It is found that the surface does reduce the entropic cost of folding for tethered peptides, as the theory suggests. This effect, however, does not always translate into improved stability because the surface may also have a destabilizing enthalpic effect. The theory neglects this effect and assumes that the enthalpy of folding is the same on and off the surface. Both the enthalpic and entropic contributions to the stability are found to be topology- and tether-placement-specific; we show that stability cannot be predicted a priori. A detailed analysis of the folding of protein A shows how the same protein can be both stabilized and destabilized on a surface depending upon how the tethering enhances or hinders the ability of the peptide to form correct tertiary structures.

Predicting the order in which contacts are broken during single molecule protein stretching experiments

We combine two methods to enable the prediction of the order in which contacts are broken under external stretching forces in single molecule experiments. These two methods are Gō-like models and elastic network models. The Gō-like models have shown remarkable success in representing many aspects of protein behavior, including the reproduction of experimental data obtained from atomic force microscopy. The simple elastic network models are often used successfully to predict the fluctuations of residues around their mean positions, comparing favorably with the experimentally measured crystallographic B-factors. The behavior of biomolecules under external forces has been demonstrated to depend principally on their elastic properties and the overall shape of their structure. We have studied in detail the muscle protein titin and green fluorescent protein and tested for ten other proteins. First, we stretch the proteins computationally by performing stochastic dynamics simulations with the Gō-like model. We obtain the force-displacement curves and unfolding scenarios of possible mechanical unfolding. We then use the elastic network model to calculate temperature factors (B-factors) and compare the slowest modes of motion for the stretched proteins and compare them with the predicted order of breaking contacts between residues in the Gō-like model. Our results show that a simple Gaussian network model is able to predict contacts that break in the next time stage of stretching. Additionally, we have found that the contact disruption is strictly correlated with the highest force exerted by the backbone on these residues. Our prediction of bond-breaking agrees well with the unfolding scenario obtained with the Gō-like model. We anticipate that this method will be a useful new tool for interpreting stretching experiments.

Toward a molecular understanding of the anisotropic response of proteins to external forces: insights from elastic network models

With recent advances in single-molecule manipulation techniques, it is now possible to measure the mechanical resistance of proteins to external pulling forces applied at specific positions. Remarkably, such recent studies demonstrated that the pulling/stretching forces required to initiate unfolding vary considerably depending on the location of the application of the forces, unraveling residue/position-specific response of proteins to uniaxial tension. Here we show that coarse-grained elastic network models based on the topology of interresidue contacts in the native state can satisfactory explain the relative sizes of such stretching forces exerted on different residue pairs. Despite their simplicity, such models presumably capture a fundamental property that dominates the observed behavior: deformations that can be accommodated by the relatively lower frequency modes of motions intrinsically favored by the structure require weaker forces and vice versa. The mechanical response of proteins to external stress is therefore shown to correlate with the anisotropic fluctuation dynamics intrinsically accessible in the folded state. The dependence on the overall fold implies that evolutionarily related proteins sharing common structural features tend to possess similar mechanical properties. However, the theory cannot explain the differences observed in a number of structurally similar but sequentially distant domains, such as the fibronectin domains.

Stretching to understand proteins - a survey of the protein data bank

We make a survey of resistance of 7510 proteins to mechanical stretching at constant speed as studied within a coarse-grained molecular dynamics model. We correlate the maximum force of resistance with the native structure, predict proteins which should be especially strong, and identify the nature of their force clamps.

Stretching of proteins in a uniform flow

Stretching of a protein by a fluid flow is compared to that in a force-clamp apparatus. The comparison is made within a simple topology-based dynamical model of a protein in which the effects of the flow are implemented using Langevin dynamics. We demonstrate that unfolding induced by a uniform flow shows a richer behavior than that in the force clamp. The dynamics of unfolding is found to depend strongly on the selection of the amino acid, usually one of the termini, which is anchored. These features offer potentially wider diagnostic tools to investigate structure of proteins compared to experiments based on the atomic force microscopy.

Mechanical unfolding of ubiquitin molecules

Mechanical stretching of ubiquitin and of its several repeats are studied through molecular-dynamics simulations. A Go-type model [H. Abe and N. Go, Biopolymers 20, 1013 (1981)] with a realistic contact map and with Lennard-Jones contact interactions is used. The model qualitatively reproduces the experimentally observed differences between force-extension patterns obtained on polyubiquitins stretched by various linkages. The terminal-to-terminal stretching of polyubiquitin results in peak forces similar to those measured for titin-based polyproteins and of a magnitude that matches measurements. Consistent with the experimental measurements, the simulated peak forces depend on the pulling speed logarithmically when thermal fluctuations are explicitly introduced. These results validate the application of topology-based models in the study of the mechanical stretching of proteins.

Protein folding in a force clamp

Kinetics of folding of a protein held in a force clamp are compared to an unconstrained folding. The comparison is made within a simple topology-based dynamical model of ubiquitin. We demonstrate that the experimentally observed variations in the end-to-end distance reflect microscopic events during folding. However, the folding scenarios in and out of the force clamp are distinct.

Mechanical unfolding revisited through a simple but realistic model

Single-molecule experiments and their application to probe the mechanical resistance and related properties of proteins provide a new dimension in our knowledge of these important and complex biological molecules. Single-molecule techniques may not have yet overridden solution experiments as a method of choice to characterize biophysical and biological properties of proteins, but have stimulated a debate and contributed considerably to bridge theory and experiment. Here we demonstrate this latter contribution by illustrating the reach of some theoretical findings using a solvable but nontrivial molecular model whose properties are analogous to those of the corresponding experimental systems. In particular, we show the relationship between the thermodynamic and the mechanical properties of a protein. The simulations presented here also illustrate how forced and spontaneous unfolding occur through different pathways and that folding and unfolding rates at equilibrium cannot in general be obtained from forced unfolding experiments or simulations. We also study the relationship between the energy surface and the mechanical resistance of a protein and show how a simple analysis of the native state can predict much of the mechanical properties of a protein.

Ising-like model for protein mechanical unfolding

The mechanical unfolding of proteins is studied by extending the Wako-Saitô-Muñoz-Eaton model. This model is generalized by including an external force, and its thermodynamics turns out to be exactly solvable. We consider two molecules, the 27th immunoglobulin domain of titin and protein PIN1. We determine equilibrium force-extension curves for the titin and study the mechanical unfolding of this molecule, finding good agreement with experiments. By using an extended form of the Jarzynski equality, we compute the free energy landscape of the PIN1 as a function of the molecule length.

Pulling single bacteriorhodopsin out of a membrane: Comparison of simulation and experiment

Mechanical unfolding of single bacteriorhodopsins from a membrane bilayer is studied using molecular dynamics simulations. The initial conformation of the lipid membrane is determined through all-atom simulations and then its coarse-grained representation is used in the studies of stretching. A Go-like model with a realistic contact map and with Lennard-Jones contact interactions is applied to model the protein-membrane system. The model qualitatively reproduces the experimentally observed differences between force-extension patterns obtained on bacteriorhodopsin at different temperatures and predicts a lack of symmetry in the choice of the terminus to pull by. It also illustrates the decisive role of the interactions of the protein with the membrane in determining the force pattern and thus the stability of transmembrane proteins.

Thermal unfolding of proteins

Thermal unfolding of proteins is compared to folding and mechanical stretching in a simple topology-based dynamical model. We define the unfolding time and demonstrate its low-temperature divergence. Below a characteristic temperature, contacts break at separate time scales and unfolding proceeds approximately in a way reverse to folding. Features in these scenarios agree with experiments and atomic simulations on titin.

Site-Specific Cross-Linking of Proteins through Tyrosine Hexahistidine Tags

The genetic addition of hexahistidine (H(6)) tags is widely used to isolate recombinant proteins by immobilized metal-affinity chromatography (IMAC). Addition of a tyrosine residue to H(6) tags enabled proteins to be covalently cross-linked under mild conditions in a manner similar to the natural, site-specific cross-linking of tyrosines into dityrosine. A series of seven hexahistidine tags with tyrosines placed in various positions (H(6)Y tags) were added to the amino terminus of the I28 immunoglobulin domain of the human cardiac titin. The H(6)Y-tagged I28 dimerized in the presence of excess Ni(2+) with a K(D) of 200 microM. Treatment of Ni(2+)-dimerized H(6)Y-I28 with an oxidant, monoperoxyphthalic acid (MMPP) or sodium sulfite, resulted in covalent protein multimerization through chelated Ni(2+)-catalyzed cross-linking of the Y residues engineered into the H(6) tag. The protein oligomerization was observed by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS PAGE). The presence of dityrosine in the cross-linked proteins was confirmed by fluorescence emission at 410 nm. Proteins lacking the Y residue in the H(6) tag treated with the same oxidative conditions did not cross-link or exhibit dityrosine fluorescence, despite the presence of an endogenous Y residue. The method may have potential uses in other protein conjugation applications such as protein labeling and interfacial immobilization of proteins on artificial surfaces.

Mechanical resistance of proteins explained using simple molecular models

Recent experiments have demonstrated that proteins unfold when two atoms are mechanically pulled apart, and that this process is different to when heated or when a chemical denaturant is added to the solution. Experiments have also shown that the response of proteins to external forces is very diverse, some of them being "hard," and others "soft." Mechanical resistance originates from the presence of barriers on the energy landscape; together, experiment and simulation have demonstrated that unfolding occurs through alternative pathways when different pairs of atoms undergo mechanical extension. Here we use simulation to probe the mechanical resistance of six structurally diverse proteins when pulled in different directions. For this, we use two very different models: a detailed, transferable one, and a coarse-grained, structure-based one. The coarse-grained model gives results that are surprisingly similar to the detailed one and qualitatively agree with experiment; i.e., the mechanical resistance of different proteins or of a single protein pulled in different directions can be predicted by simulation. The results demonstrate the importance of pulling direction relative to the local topology in determining mechanical stability, and rationalize the effect of the location of importation/degradation tags on the rates of mitochondrial import or protein degradation in vivo.