Tensorial elastic network model for protein dynamics: Integration of the anisotropic network model with bond‐bending and twist elasticities

We present a tensorial elastic network model (TNM) to describe the equilibrium fluctuations of proteins near their native fold structure. The model combines the anisotropic network model (ANM), bond bending elasticity, and backbone twist elasticity, and can predict both the isotropic fluctuations, similar to the Gaussian network model (GNM), and anisotropic fluctuations, similar to the ANM. TNM performs equally well for B‐factor predictions as GNM and predicts the anisotropy of B‐factors better than ANM. The model also outperforms the ANM in its predictability of the complete anisotropic displacement parameters. Proteins 2012; © 2012 Wiley Periodicals, Inc.     

All‐atom and coarse‐grained simulations of the forced unfolding pathways of the SNARE complex

The SNARE complex, consisting of three proteins (VAMP2, syntaxin, and SNAP‐25), is thought to drive membrane fusion by assembling into a four‐helix bundle through a zippering process. In support of the above zippering model, a recent single‐molecule optical tweezers experiment by Gao et al. revealed a sequential unzipping of SNARE along VAMP2 in the order of the linker domain → the C‐terminal domain → the N‐terminal domain. To offer detailed structural insights to this unzipping process, we have performed all‐atom and coarse‐grained steered molecular dynamics (sMD) simulations of the forced unfolding pathways of SNARE using different models and force fields. Our findings are summarized as follows: First, the sMD simulations based on either an all‐atom force field (with an implicit solvent model) or a coarse‐grained Go model were unable to capture the forced unfolding pathway of SNARE as observed by Gao et al., which may be attributed to insufficient simulation time and inaccurate force fields. Second, the sMD simulations based on a reparameterized coarse‐grained model (i.e., modified elastic network model) were able to predict a sequential unzipping of SNARE in good agreement with the findings by Gao et al. The key to this success is to reparameterize the intrahelix and interhelix nonbonded force constants against the pair‐wise residue–residue distance fluctuations collected from all‐atom MD simulations of SNARE. Therefore, our finding supports the importance of accurately describing the inherent dynamics/flexibility of SNARE (in the absence of force), in order to correctly simulate its unfolding behaviors under force. This study has established a useful computational framework for future studies of the zippering function of SNARE and its perturbations by point mutations with amino‐acid level of details, and more generally the forced unfolding pathways of other helix bundle proteins. Proteins 2014; 82:1376–1386. © 2014 Wiley Periodicals, Inc.     

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.     

Predicting Atomic Details of the Unfolding Pathway for YibK,a Knotted Protein from the SPOUT Superfamily

Abstract

Several protein structures have been reported to contain intricate knots of the polypeptide backbone but the mechanism of the (un)folding process of knotted proteins remains unknown. The members of the SPOUT superfamily of RNA methyltransferases are some of the most intensely studied systems for investigation of the knot formation and function. YibK (whose biochemical function remains unknown) is the representative protein of the SPOUT superfamily. This protein exhibits a deep trefoil knot at the C-terminus.

We conducted an extensive computational analysis of the unfolding process for the monomeric form of YibK. In order to predict the (un)folding pathway of YibK, we have calculated the order of secondary structure disassembly using UNFOLD, and performed thermal unfolding simulations using classical Molecular Dynamics (MD), as well as simulations employing reduced representation of the peptide chain using either MD with the UNRES method or the Monte Carlo (MC) unfolding with the REFINER method.

Results obtained from all methods used in this work are in qualitative agreement. We found that YibK unfolds through four intermediate states. The trefoil knot in YibK disappears at the end of the unfolding process, long after the protein loses its native topology. We observed that the C-terminus leaves the knotting loop folded into a hairpin-like structure, in agreement with the results of coarse-grained simulation reported earlier. We propose that the folding pathway of YibK corresponds to the reversed sequence of events observed in the unfolding pathway elucidated in this study. Thus, we predict that the knot formation is the slowest part of the YibK folding process.     

Soft Vibrational Modes Predict Breaking Events during Force-Induced Protein Unfolding

We investigate the correlation between soft vibrational modes and unfolding events in simulated force spectroscopy of proteins. Unfolding trajectories are obtained from molecular dynamics simulations of a Gō model of a monomer of a mutant of superoxide dismutase 1 protein containing all heavy atoms in the protein, and a normal mode analysis is performed based on the anisotropic network model. We show that a softness map constructed from the superposition of the amplitudes of localized soft modes correlates with unfolding events at different stages of the unfolding process. Soft residues are up to eight times more likely to undergo disruption of native structure than the average amino acid. The memory of the softness map is retained for extensions of up to several nanometers, but decorrelates more rapidly during force drops.     

Essential dynamics of the cold denaturation: pressure and temperature effects in yeast frataxin

The cold denaturation of globular proteins is a process that can be caused by increasing pressure or decreasing the temperature. Currently, the action mechanism of this process has not been clearly understood, raising an interesting debate on the matter. We have studied the process of cold denaturation using molecular dynamics simulations of the frataxin system Yfh1, which has a dynamic experimental characterization of unfolding at low and high temperatures. The frataxin model here studied allows a comparative analysis using experimental data. Furthermore, we monitored the cold denaturation process of frataxin and also investigated the effect under the high‐pressure regime. For a better understanding of the dynamics and structural properties of the cold denaturation, we also analyzed the MD trajectories using essentials dynamic. The results indicate that changes in the structure of water by the effect of pressure and low temperatures destabilize the hydrophobic interaction modifying the solvation and the system volume leading to protein denaturation. Proteins 2016; 85:125–136. © 2016 Wiley Periodicals, Inc.     

Inferring a weighted elastic network from partial unfolding with coarse‐grained simulations

A number of studies have demonstrated that simple elastic network models can reproduce experimental B‐factors, providing insights into the structure–function properties of proteins. Here, we report a study on how to improve an elastic network model and explore its performance by predicting the experimental B‐factors. Elastic network models are built on the experimental coordinates, and they only take the pairs of atoms within a given cutoff distance r_c into account. These models describe the interactions by elastic springs with the same force constant. We have developed a method based on numerical simulations with a simple coarse‐grained force field, to attribute weights to these spring constants. This method considers the time that two atoms remain connected in the network during partial unfolding, establishing a means of measuring the strength of each link. We examined two different coarse‐grained force fields and explored the computation of these weights by unfolding the native structures. Proteins 2014; 82:119–129. © 2013 Wiley Periodicals, Inc.     

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.     

High temperature unfolding of a truncated hemoglobin by molecular dynamics simulation

Heme containing proteins are associated with peroxidase activity. The proteins like hemoglobin, myoglobins, cytochrome c and micro-peroxidase other than peroxidases have been shown to exhibit weak peroxidase-like activity. This weak peroxidase–like activity in hemoglobin-like molecules is due to heme moiety. We conducted molecular dynamics (MD) studies to decipher the unfolding path of Ba-Glb (a truncated hemoglobin from Bacillus anthracis) and the role of heme moiety to its unfolding path. The similar unfolding path is also observed in vitro by UV/VIS spectroscopy. The data confirmed that the unfolding of Ba-Glb follows a three state process with a meta-stable (intermediate) state between the native and unfolded conformations. The present study is supported by several unfolding parameters like root-mean-square-deviation (RMSD), dictionary of protein secondary structure (DSSP), and free energy landscape. Understanding the structure of hemoglobin like proteins in unicellular dreaded pathogens like B. anthracis will pave way for newer drug discovery targets and in the disease management of anthrax.     

Effect of circular permutations on transient partial unfolding in proteins

Under native conditions, proteins can undergo transient partial unfolding, which may cause proteins to misfold or aggregate. A change in sequence connectivity by circular permutation may affect the energetics of transient partial unfolding in proteins without altering the three‐dimensional structures. Using Escherichia coli dihydrofolate reductase (DHFR) as a model system, we investigated how circular permutation affects transient partial unfolding in proteins. We constructed three circular permutants, CP18, CP37, and CP87, with the new N‐termini at residue 18, 37, and 87, respectively, and probed transient partial unfolding by native‐state proteolysis. The new termini in CP18, CP37, and CP87 are within, near, and distal to the Met20 loop, which is known to be dynamic and also part of the region that undergoes transient unfolding in wild‐type DHFR. The stabilities of both native and partially unfolded forms of CP18 are similar to those of wild‐type DHFR, suggesting that the influence of introducing new termini in a dynamic region to the protein is minimal. CP37 has a significantly more accessible partially unfolded form than wild‐type DHFR, demonstrating that introducing new termini near a dynamic region may promote transient partial unfolding. CP87 has significantly destabilized native and partially unfolded forms, confirming that modification of the folded region in a partially unfolded form destabilizes the partially unfolded form similar to the native form. Our findings provide valuable guidelines to control transient partial unfolding in designing circular permutants in proteins.     

The Mechanical Properties of PCNA: Implications for the Loading and Function of a DNA Sliding Clamp

Sliding clamps are toroidal proteins that encircle DNA and act as mobile platforms for DNA replication and repair machinery. To be loaded onto DNA, the eukaryotic sliding clamp Proliferating Cell Nuclear Antigen (PCNA) must be splayed open at one of the subunit-subunit interfaces by the ATP-dependent clamp loader, Replication Factor C, whose clamp-interacting sites form a right-handed spiral. Earlier molecular dynamics (MD) studies suggested that when PCNA opens, it preferentially adopts a right-handed spiral to match the spiral of the clamp loader. Here, analysis of considerably longer MD simulations shows that although the opened form of PCNA can achieve conformations matching the helical pitch of Replication Factor C, it is not biased toward a right-handed spiral structure. A coarse-grained elastic model was also built; its strong correspondence to the all-atom MD simulations of PCNA suggests that the behavior of the open clamp is primarily due to elastic deformation governed by the topology of the clamp domains. The elastic model was further used to construct the energy landscape of the opened PCNA clamp, including conformations that would allow binding to the clamp loader and loading onto double-stranded DNA. A picture of PCNA emerges of a rather flexible protein that, once opened, is mechanically compliant in the clamp opening process.     

Probing possible downhill folding: native contact topology likely places a significant constraint on the folding cooperativity of proteins with approximately 40 residues

Experiments point to appreciable variations in folding cooperativity among natural proteins with approximately 40 residues, indicating that the behaviors of these proteins are valuable for delineating the contributing factors to cooperative folding. To explore the role of native topology in a protein's propensity to fold cooperatively and how native topology might constrain the degree of cooperativity achievable by a given set of physical interactions, we compared folding/unfolding kinetics simulated using three classes of native-centric C^α chain models with different interaction schemes. The approach was applied to two homologous 45-residue fragments from the peripheral subunit-binding domain family and a 39-residue fragment of the N-terminal domain of ribosomal protein L9. Free-energy profiles as functions of native contact number were computed to assess the heights of thermodynamic barriers to folding. In addition, chevron plots of folding/unfolding rates were constructed as functions of native stability to facilitate comparison with available experimental data. Although common Gō-like models with pairwise Lennard-Jones-type interactions generally fold less cooperatively than real proteins, the rank ordering of cooperativity predicted by these models is consistent with experiment for the proteins investigated, showing increasing folding cooperativity with increasing nonlocality of a protein's native contacts. Models that account for water-expulsion (desolvation) barriers and models with many-body (nonadditive) interactions generally entail higher degrees of folding cooperativity indicated by more linear model chevron plots, but the rank ordering of cooperativity remains unchanged. A robust, experimentally valid rank ordering of model folding cooperativity independent of the multiple native-centric interaction schemes tested here argues that native topology places significant constraints on how cooperatively a protein can fold.     

Unfolding of globular proteins: monte carlo dynamics of a realistic reduced model

Reduced lattice models of proteins and Monte Carlo dynamics were used to simulate the initial stages of the unfolding of several proteins of various structural types, and the results were compared to experiment. The models semiquantitatively reproduce the approximate order of events of unfolding as well as subtle mutation effects and effects resulting from differences in sequences of similar folds. The short-time mobility of particular residues, observed in simulations, correlates with the crystallographic temperature factor. The main factor controlling unfolding is the native state topology, with sequence playing a less important role. The correlation with various experiments, especially for sequence-specific effects, strongly suggests that properly designed reduced models of proteins can be used for qualitative studies (or prediction) of protein unfolding pathways.     

Folding and stability of sweet protein single-chain monellin. An insight to protein engineering

Engineered single-chain monellin (SCM) proteins were constructed by recombinant technology without disrupting the topology and sweet activity of native protein. Data from 8-anilinonaphthalene-1-sulfonic acid fluorescence, size-exclusion chromatography, and heteronuclear NMR strongly suggest the presence of a folding intermediate at 1.5 m GdnHCl for SCM protein. The structural feature of the folding intermediate from NMR data reveals that the secondary structures became mostly unstable, and protein experiences a dynamic equilibrium between native and unfolded state. All backbone amide protons exchange within 10 min, which imply that no stable hydrogen bonds exist in the secondary structural regions in the folding intermediate. From equilibrium unfolding and mutagenesis studies, the unfolding transition midpoints of mutant proteins gradually shifted toward lower denaturant concentration, indicating stability reductions of mutant proteins. Our results suggest that stability and folding pathways of SCM proteins could be regulated by a combined study of spectroscopy and mutagenesis, and these studies will provide useful information for understanding the folding kinetics of novel engineered proteins.     

Molecular force modulation spectroscopy revealing the dynamic response of single bacteriorhodopsins

Recent advances in atomic force microscopy allowed globular and membrane proteins to be mechanically unfolded on a single-molecule level. Presented is an extension to the existing force spectroscopy experiments. While unfolding single bacteriorhodopsins from native purple membranes, small oscillation amplitudes (6-9 nm) were supplied to the vertical displacement of the cantilever at a frequency of 3 kHz. The phase and amplitude response of the cantilever-protein system was converted to reveal the elastic (conservative) and viscous (dissipative) contributions to the unfolding process. The elastic response (stiffness) of the extended parts of the protein were in the range of a few tens pN/nm and could be well described by the derivative of the wormlike chain model. Discrete events in the viscous response coincided with the unfolding of single secondary structure elements and were in the range of 1 microNs/m. In addition, these force modulation spectroscopy experiments revealed novel mechanical unfolding intermediates of bacteriorhodopsin. We found that kinks result in a loss of unfolding cooperativity in transmembrane helices. Reconstructing force-distance spectra by the integration of amplitude-distance spectra verified their position, offering a novel approach to detect intermediates during the forced unfolding of single proteins.     

An analysis of the influence of protein intrinsic dynamical properties on its thermal unfolding behavior

The influence of the protein topology-encoded dynamical properties on its thermal unfolding motions was studied in the present work. The intrinsic dynamics of protein topology was obtained by the anisotropic network model (ANM). The ANM has been largely used to investigate protein collective functional motions, but it is not well elucidated if this model can also reveal the preferred large-scale motions during protein unfolding. A small protein barnase is used as a typical case study to explore the relationship between protein topology-encoded dynamics and its unfolding motions. Three thermal unfolding simulations at 500 K were performed for barnase and the entire unfolding trajectories were sampled and partitioned into several windows. For each window, the preferred unfolding motions were investigated by essential dynamics analysis, and then associated with the intrinsic dynamical properties of the starting conformation in this window, which is detected by ANM. The results show that only a few slow normal modes imposed by protein structure are sufficient to give a significant overlap with the preferred unfolding motions. Especially, the large amplitude unfolding movements, which imply that the protein jumps out of a local energy basin, can be well described by a single or several ANM slow modes. Besides the global motions, it is also found that the local residual fluctuations encoded in protein structure are highly correlated with those in the protein unfolding process. Furthermore, we also investigated the relationship between protein intrinsic flexibility and its unfolding events. The results show that the intrinsic flexible regions tend to unfold early. Several early unfolding events can be predicted by analysis of protein structural flexibility. These results imply that protein structure-encoded dynamical properties have significant influences on protein unfolding motions.     

Secondary and tertiary structure elasticity of titin Z1Z2 and a titin chain model

The giant protein titin, which is responsible for passive elasticity in muscle fibers, is built from approximately 300 regular immunoglobulin-like (Ig) domains and FN-III repeats. While the soft elasticity derived from its entropic regions, as well as the stiff mechanical resistance derived from the unfolding of the secondary structure elements of Ig- and FN-III domains have been studied extensively, less is known about the mechanical elasticity stemming from the orientation of neighboring domains relative to each other. Here we address the dynamics and energetics of interdomain arrangement of two adjacent Ig-domains of titin, Z1, and Z2, using molecular dynamics (MD) simulations. The simulations reveal conformational flexibility, due to the domain-domain geometry, that lends an intermediate force elasticity to titin. We employ adaptive biasing force MD simulations to calculate the energy required to bend the Z1Z2 tandem open to identify energetically feasible interdomain arrangements of the Z1 and Z2 domains. The finding is cast into a stochastic model for Z1Z2 interdomain elasticity that is generalized to a multiple domain chain replicating many Z1Z2-like units and representing a long titin segment. The elastic properties of this chain suggest that titin derives so-called tertiary structure elasticity from bending and twisting of its domains. Finally, we employ steered molecular dynamics simulations to stretch individual Z1 and Z2 domains and characterize the so-called secondary structure elasticity of the two domains. Our study suggests that titin's overall elastic response at weak force stems from a soft entropic spring behavior (not described here), from tertiary structure elasticity with an elastic spring constant of approximately 0.001-1 pN/A and, at strong forces, from secondary structure elasticity.     

Changes in the Fabric and Compliance Tensors of Cancellous Bone due to Trabecular Surface Remodeling,Predicted by a Digital Image-based Model

Fabric and compliance tensors of a cube of cancellous bone with a complicated three-dimensional trabecular structure were obtained for trabecular surface remodeling by using a digital image-based model combined with a large-scale finite element method. Using mean intercept length and a homogenization method, the fabric and compliance tensors were determined for the trabecular structure obtained in the computer remodeling simulation. The tensorial quantities obtained indicated that anisotropic structural changes occur in cancellous bone adapting to the compressive loading condition. There were good correlations between the fabric tensor, bone volume fraction, and compliance tensor in the remodeling process. The result demonstrates that changes in the structural and mechanical properties of cancellous bone are essentially anisotropic and should be expressed by tensorial quantities.