A structure-based model fails to probe the mechanical unfolding pathways of the titin I27 domain

Exploring the unfolding pathways of protein families using Elastic Network Model

We explore how a protein's native structure determines its unfolding process. We examine how the local structural features, like shear, and the global structural properties, like the number of soft modes, change during unfolding. Simulations are performed using a Gaussian Network Model (GNM) with bond breaking for both thermal and force-induced unfolding scenarios. We find that unfolding starts in areas of high shear in the native structure and progressively spreads to the low shear regions. Interestingly, analysis of single domain protein families (Chymotrypsin inhibitor and Barnase) reveal that proteins with distinct unfolding pathways exhibit divergent behavior of the number of soft modes during unfolding. This suggests that the number of soft modes might be a valuable tool for understanding thermal unfolding pathways. Additionally, we found a strong link between a protein's overall structural similarity (TM-score) and its unfolding pathways, highlighting the importance of the native structure in determining how a protein unfolds.

Interpretation of Single-Molecule Force Experiments on Proteins Using Normal Mode Analysis

Single-molecule force spectroscopy experiments allow protein folding and unfolding to be explored using mechanical force. Probably the most informative technique for interpreting the results of these experiments at the structural level makes use of steered molecular dynamics (MD) simulations, which can explicitly model the protein under load. Unfortunately, this technique is computationally expensive for many of the most interesting biological molecules. Here, we find that normal mode analysis (NMA), a significantly cheaper technique from a computational perspective, allows at least some of the insights provided by MD simulation to be gathered. We apply this technique to three non-homologous proteins that were previously studied by force spectroscopy: T4 lysozyme (T4L), Hsp70 and the glucocorticoid receptor domain (GCR). The NMA results for T4L and Hsp70 are compared with steered MD simulations conducted previously, and we find that we can recover the main results. For the GCR, which did not undergo MD simulation, our approach identifies substructures that correlate with experimentally identified unfolding intermediates. Overall, we find that NMA can make a valuable addition to the analysis toolkit for the structural analysis of single-molecule force experiments on proteins.

Modeling of Protein Structural Flexibility and Large-Scale Dynamics: Coarse-Grained Simulations and Elastic Network Models

Fluctuations of protein three-dimensional structures and large-scale conformational transitions are crucial for the biological function of proteins and their complexes. Experimental studies of such phenomena remain very challenging and therefore molecular modeling can be a good alternative or a valuable supporting tool for the investigation of large molecular systems and long-time events. In this minireview, we present two alternative approaches to the coarse-grained (CG) modeling of dynamic properties of protein systems. We discuss two CG representations of polypeptide chains used for Monte Carlo dynamics simulations of protein local dynamics and conformational transitions, and highly simplified structure-based elastic network models of protein flexibility. In contrast to classical all-atom molecular dynamics, the modeling strategies discussed here allow the quite accurate modeling of much larger systems and longer-time dynamic phenomena. We briefly describe the main features of these models and outline some of their applications, including modeling of near-native structure fluctuations, sampling of large regions of the protein conformational space, or possible support for the structure prediction of large proteins and their complexes.

Role of Resultant Dipole Moment in Mechanical Dissociation of Biological Complexes

Protein-peptide interactions play essential roles in many cellular processes and their structural characterization is the major focus of current experimental and theoretical research. Two decades ago, it was proposed to employ the steered molecular dynamics (SMD) to assess the strength of protein-peptide interactions. The idea behind using SMD simulations is that the mechanical stability can be used as a promising and an efficient alternative to computationally highly demanding estimation of binding affinity. However, mechanical stability defined as a peak in force-extension profile depends on the choice of the pulling direction. Here we propose an uncommon choice of the pulling direction along resultant dipole moment (RDM) vector, which has not been explored in SMD simulations so far. Using explicit solvent all-atom MD simulations, we apply SMD technique to probe mechanical resistance of ligand-receptor system pulled along two different vectors. A novel pulling direction—when ligand unbinds along the RDM vector—results in stronger forces compared to commonly used ligand unbinding along center of masses vector. Our observation that RDM is one of the factors influencing the mechanical stability of protein-peptide complex can be used to improve the ranking of binding affinities by using mechanical stability as an effective scoring function.

Kinetics and mechanical stability of the fibril state control fibril formation time of polypeptide chains: A computational study

Fibril formation resulting from protein misfolding and aggregation is a hallmark of several neurodegenerative diseases such as Alzheimer's and Parkinson's diseases. Despite much progress in the understanding of the protein aggregation process, the factors governing fibril formation rates and fibril stability have not been fully understood. Using lattice models, we have shown that the fibril formation time is controlled by the kinetic stability of the fibril state but not by its energy. Having performed all-atom explicit solvent molecular dynamics simulations with the GROMOS43a1 force field for full-length amyloid beta peptides Aβ₄₀ and Aβ₄₂ and truncated peptides, we demonstrated that kinetic stability can be accessed via mechanical stability in such a way that the higher the mechanical stability or the kinetic stability, the faster the fibril formation. This result opens up a new way for predicting fibril formation rates based on mechanical stability that may be easily estimated by steered molecular dynamics.

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.

As Simple As Possible, but Not Simpler: Exploring the Fidelity of Coarse-Grained Protein Models for Simulated Force Spectroscopy

Mechanical unfolding of a single domain of loop-truncated superoxide dismutase protein has been simulated via force spectroscopy techniques with both all-atom (AA) models and several coarse-grained models having different levels of resolution: A Gō model containing all heavy atoms in the protein (HA-Gō), the associative memory, water mediated, structure and energy model (AWSEM) which has 3 interaction sites per amino acid, and a Gō model containing only one interaction site per amino acid at the Cα position (Cα-Gō). To systematically compare results across models, the scales of time, energy, and force had to be suitably renormalized in each model. Surprisingly, the HA-Gō model gives the softest protein, exhibiting much smaller force peaks than all other models after the above renormalization. Clustering to render a structural taxonomy as the protein unfolds showed that the AA, HA-Gō, and Cα-Gō models exhibit a single pathway for early unfolding, which eventually bifurcates repeatedly to multiple branches only after the protein is about half-unfolded. The AWSEM model shows a single dominant unfolding pathway over the whole range of unfolding, in contrast to all other models. TM alignment, clustering analysis, and native contact maps show that the AWSEM pathway has however the most structural similarity to the AA model at high nativeness, but the least structural similarity to the AA model at low nativeness. In comparison to the AA model, the sequence of native contact breakage is best predicted by the HA-Gō model. All models consistently predict a similar unfolding mechanism for early force-induced unfolding events, but diverge in their predictions for late stage unfolding events when the protein is more significantly disordered.

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.