Related Articles

Single-molecule force spectroscopy is providing unique, and sometimes unexpected, insights into the free-energy landscapes of proteins. Despite the complexity of the free-energy landscapes revealed by mechanical probes, forced unfolding experiments are often analyzed using one-dimensional models that predict a logarithmic dependence of the unfolding force on the pulling velocity. We previously found that the unfolding force of the protein filamin at low pulling speed did not decrease logarithmically with the pulling speed. Here we present results from a large number of unfolding simulations of a coarse-grain model of the protein filamin under a broad range of constant forces. These show that a two-path model is physically plausible and produces a deviation from the behavior predicted by one-dimensional models analogous to that observed experimentally. We also show that the analysis of the distributions of unfolding forces (p[F]) contains crucial and exploitable information, and that a proper description of the unfolding of single-domain proteins needs to account for the intrinsic multidimensionality of the underlying free-energy landscape, especially when the applied perturbation is small.

Single-molecule force spectroscopy is providing unique, and sometimes unexpected, insights into the free-energy landscapes of proteins. Despite the complexity of the free-energy landscapes revealed by mechanical probes, forced unfolding experiments are often analyzed using one-dimensional models that predict a logarithmic dependence of the unfolding force on the pulling velocity. We previously found that the unfolding force of the protein filamin at low pulling speed did not decrease logarithmically with the pulling speed. Here we present results from a large number of unfolding simulations of a coarse-grain model of the protein filamin under a broad range of constant forces. These show that a two-path model is physically plausible and produces a deviation from the behavior predicted by one-dimensional models analogous to that observed experimentally. We also show that the analysis of the distributions of unfolding forces (p[F]) contains crucial and exploitable information, and that a proper description of the unfolding of single-domain proteins needs to account for the intrinsic multidimensionality of the underlying free-energy landscape, especially when the applied perturbation is small.

Summary

Experimental observation has led to the commonly held view that native state protein topology is the principle determinant of mechanical strength. However, the PKD domains of polycystin-1 challenge this assumption: they are stronger than predicted from their native structure. Molecular dynamics simulations suggest that force induces rearrangement to an intermediate structure, with nonnative hydrogen bonds, that resists unfolding. Here we test this hypothesis directly by introducing mutations designed to prevent formation of these nonnative interactions. We find that these mutations, which only moderately destabilize the native state, reduce the mechanical stability dramatically. The results demonstrate that nonnative interactions impart significant mechanical stability, necessary for the mechanosensor function of polycystin-1. Remarkably, such nonnative interactions result from force-induced conformational change: the PKD domain is strengthened by the application of force.

All-atom explicit-solvent molecular dynamics simulations are used to pull with extremely large constant force (750–3000 pN) on three small proteins. The introduction of a nondimensional timescale permits direct comparison of unfolding across all forces. A crossover force of approximately 1100 pN divides unfolding dynamics into two regimes. At higher forces, residues sequentially unfold from the pulling end while maintaining the remainder of the protein force-free. Measurements of hydrodynamic viscous stresses are made easy by the high speeds of unfolding. Using an exact low-Reynolds-number scaling, these measurements can be extrapolated to provide, for the first time, an estimate of the hydrodynamic force on low-force unfolding. Below 1100 pN, but surprisingly still at extremely large applied force, intermediate states and cooperative unfoldings as seen at much lower forces are observed. The force-insensitive persistence of these structures indicates that decomposition into unfolded fragments requires a large fluctuation. This finding suggests how proteins are constructed to resist transient high force. The progression of helix and sheet unfolding is also found to be insensitive to force. The force-insensitivity of key aspects of unfolding opens the possibility that numerical simulations can be accelerated by high applied force while still maintaining critical features of unfolding.

Summary

Discrete molecular dynamics (DMD) is a rapid sampling method used in protein folding and aggregation studies. Until now, DMD was used to perform simulations of simplified protein models in conjunction with structure-based force fields. Here, we develop an all-atom protein model and a transferable force field featuring packing, solvation, and environment-dependent hydrogen bond interactions. Using the replica exchange method, we perform folding simulations of six small proteins (20–60 residues) with distinct native structures. In all cases, native or near-native states are reached in simulations. For three small proteins, multiple folding transitions are observed and the computationally-characterized thermodynamics are in quantitative agreement with experiments. The predictive power of all-atom DMD highlights the importance of environment-dependent hydrogen bond interactions in modeling protein folding. The developed approach can be used for accurate and rapid sampling of conformational spaces of proteins and protein-protein complexes, and applied to protein engineering and design of protein-protein interactions.

Reports on the ultrastructure of cells as well as biochemical data have, for several years, been indicating a connection between caveolae and the actin cytoskeleton. Here, using a yeast two-hybrid approach, we have identified the F-actin cross-linking protein filamin as a ligand for the caveolae-associated protein caveolin-1. Binding of caveolin-1 to filamin involved the N-terminal region of caveolin-1 and the C terminus of filamin close to the filamin-dimerization domain. In in vitro binding assays, recombinant caveolin-1 bound to both nonmuscle and muscle filamin, indicating that the interaction might not be cell type specific. With the use of confocal microscopy, colocalization of caveolin-1 and filamin was observed in elongated patches at the plasma membrane. Remarkably, when stress fiber formation was induced with Rho-stimulating Escherichia coli cytotoxic necrotizing factor 1, the caveolin-1–positive structures became coaligned with stress fibers, indicating that there was a physical link connecting them. Immunogold double-labeling electron microscopy confirmed that caveolin-1–labeled racemose caveolae clusters were positive for filamin. The actin network, therefore, seems to be directly involved in the spatial organization of caveolin-1–associated membrane domains.

The role of mechanical force in cellular processes is increasingly revealed by single molecule experiments and simulations of force-induced transitions in proteins. How the applied force propagates within proteins determines their mechanical behavior yet remains largely unknown. We present a new method based on molecular dynamics simulations to disclose the distribution of strain in protein structures, here for the newly determined high-resolution crystal structure of I27, a titin immunoglobulin (IG) domain. We obtain a sparse, spatially connected, and highly anisotropic mechanical network. This allows us to detect load-bearing motifs composed of interstrand hydrogen bonds and hydrophobic core interactions, including parts distal to the site to which force was applied. The role of the force distribution pattern for mechanical stability is tested by in silico unfolding of I27 mutants. We then compare the observed force pattern to the sparse network of coevolved residues found in this family. We find a remarkable overlap, suggesting the force distribution to reflect constraints for the evolutionary design of mechanical resistance in the IG family. The force distribution analysis provides a molecular interpretation of coevolution and opens the road to the study of the mechanism of signal propagation in proteins in general.

Author Summary

Many biological processes such as cell proliferation and signaling are guided by mechanical stress. Proteins as the molecular machinery behind these processes are reacting to or withstanding mechanical forces in specific ways. How mechanical stress propagates through proteins to induce a certain mechanical response is currently unknown. We here present a new method that detects force distribution in proteins, reminiscent of computational approaches used to engineer macroscopic structures. The method is based on molecular dynamics simulations during which we calculate changes in interatomic forces, here caused by pulling the protein. We apply this method to the extremely robust immunoglobulin domain of the muscle protein titin and obtain the mechanical network of stress-bearing elements spanning the protein scaffold. Mutations in this region were shown to result in a significant loss of mechanical stability. We then ask how the remarkable mechanical stability has been designed during evolution. To this end, we compare the force distribution network with a network of coevolved residues found in the immunoglobulin family. Both networks show a remarkable overlap, thereby suggesting that the observed stress propagation pattern reflects constraints used during evolution to render this protein mechanically robust.

Granular matter is rigid when jammed, and flows under external loads. Here temperature- and force-unfolding molecular dynamics stimulations are used to demonstrate that proteins display features of jamming, characterized by a force distribution peak on folding and a slowdown of stress relaxation.

Grains and glasses, widely different materials, arrest their motions upon decreasing temperature and external load, respectively, in common ways, leading to a universal jamming phase diagram conjecture. However, unified theories are lacking, mainly because of the disparate nature of the particle interactions. Here we demonstrate that folded proteins exhibit signatures common to both glassiness and jamming by using temperature- and force-unfolding molecular dynamics simulations. Upon folding, proteins develop a peak in the interatomic force distributions that falls on a universal curve with experimentally measured forces on jammed grains and droplets. Dynamical signatures are found as a dramatic slowdown of stress relaxation upon folding. Together with granular similarities, folding is tied not just to the jamming transition, but a more nuanced picture of anisotropy, preparation protocol and internal interactions emerges. Results have implications for designing stable polymers and can open avenues to link protein folding to jamming theory.

The early stages of the thermal unfolding of apoflavodoxin have been determined by using atomistic multi microsecond-scale molecular dynamics (MD) simulations complemented with a variety of experimental techniques. Results strongly suggest that the intermediate is reached very early in the thermal unfolding process and that it has the properties of an “activated” form of the native state, where thermal fluctuations in the loops break loop-loop contacts. The unrestrained loops gain then kinetic energy corrupting short secondary structure elements without corrupting the core of the protein. The MD-derived ensembles agree with experimental observables and draw a picture of the intermediate state inconsistent with a well-defined structure and characteristic of a typical partially disordered protein. Our results allow us to speculate that proteins with a well packed core connected by long loops might behave as partially disordered proteins under native conditions, or alternatively behave as three state folders. Small details in the sequence, easily tunable by evolution, can yield to one or the other type of proteins.

Author Summary

A simplistic view of protein structure tends to emphasize the opposition between the native state and the denatured ensemble of unfolded conformations. In addition to these extreme conformations, proteins subjected to a variety of perturbations often populate alternative partly unfolded conformations, some of which are close in energy to the native state and, accordingly, can be populated under native or quasi-native conditions. There is increasing evidence that these “perturbed” conformations participate in protein function or, in some cases, are related to the outcome of folding diseases. We have used the “state of the art” molecular dynamics combined with a variety of experimental techniques to characterize for the first time, to our knowledge, the thermal intermediate of a three-state folding protein (apoflavodoxin). Based on our results we have been able to suggest a general mechanism of thermal unfolding in complex proteins and to determine interesting links between thermal intermediates and partially unfolded proteins.

von Willebrand factor (VWF) multimers mediate primary adhesion and aggregation of platelets. VWF potency critically depends on multimer size, which is regulated by a feedback mechanism involving shear-induced unfolding of the VWF-A2 domain and cleavage by the metalloprotease ADAMTS-13. Here we report crystallographic and single-molecule optical tweezers data on VWF-A2 providing mechanistic insight into calcium-mediated stabilization of the native conformation that protects A2 from cleavage by ADAMTS-13. Unfolding of A2 requires higher forces when calcium is present and primarily proceeds through a mechanically stable intermediate with non-native calcium coordination. Calcium further accelerates refolding markedly, in particular, under applied load. We propose that calcium improves force sensing by allowing reversible force switching under physiologically relevant hydrodynamic conditions. Our data show for the first time the relevance of metal coordination for mechanical properties of a protein involved in mechanosensing.

von Willebrand factor (VWF) multimers mediate primary adhesion and aggregation of platelets. Jakobi et al. reveal a calcium-binding site in the VWF-A2 domain, and show that calcium binding encourages folding of the protein and has a role in mechanosensing.

Both folded and unfolded conformations should be observed for a protein at its melting temperature, or Tm, where the ΔG between these states is zero. In an all-atom molecular dynamics simulation of chymotrypsin inhibitor 2 (CI2) at its experimental Tm, the protein rapidly loses its low-temperature native structure, it then unfolds before refolding to a stable, native-like conformation. The initial unfolding follows the unfolding pathway described previously for higher temperature simulations: the hydrophobic core is disrupted, the β-sheet pulls apart and the α-helix unravels. The unfolded state reached under these conditions maintains a kernel of structure in the form of a nonnative hydrophobic cluster. Refolding simply reverses this path, the side chain interactions shift, the helix refolds, and the native packing and hydrogen bonds are recovered. The end result of this refolding is not the initial crystal structure; it contains the proper topology and the majority of the native contacts, but the structure is expanded and the contacts are long. We believe this state to be the native one at elevated temperature and the change in volume and contact lengths is consistent with experimental studies of other native proteins at elevated temperature and the chemical denaturant equivalent of Tm.

We have studied the mechanical properties of the immunoglobulin-binding domain of protein G at the atomic level under stretching at constant velocity using molecular dynamics simulations. We have found that the unfolding process can occur either in a single step or through intermediate states. Analysis of the trajectories from the molecular dynamic simulations showed that the mechanical unfolding of the immunoglobulin-binding domain of protein G is triggered by the separation of the terminal β-strands and the order in which the secondary-structure elements break is practically the same in two- and multi-state events and at the different extension velocities studied. It is seen from our analysis of 24 trajectories that the theoretical pathway of mechanical unfolding for the immunoglobulin-binding domain of protein G does not coincide with that proposed in denaturant studies in the absence of force.

Ankyrin repeat proteins are elastic materials that unfold and refold sequentially, repeat by repeat, under force. Herein we use atomistic molecular dynamics to compare the mechanical properties of the 7-ankyrin-repeat oncoprotein Gankyrin in isolation and in complex with its binding partner S6-C. We show that the bound S6-C greatly increases the resistance of Gankyrin to mechanical stress. The effect is specific to those repeats of Gankyrin directly in contact with S6-C, and the mechanical ‘hot spots’ of the interaction map to the same repeats as the thermodynamic hot spots. A consequence of stepwise nature of unfolding and the localized nature of ligand binding is that it impacts on all aspects of the protein's mechanical behavior, including the order of repeat unfolding, the diversity of unfolding pathways accessed, the nature of partially unfolded intermediates, the forces required and the work transferred to the system to unfold the whole protein and its parts. Stepwise unfolding thus provides the means to buffer repeat proteins and their binding partners from mechanical stress in the cell. Our results illustrate how ligand binding can control the mechanical response of proteins. The data also point to a cellular mechano-switching mechanism whereby binding between two partner macromolecules is regulated by mechanical stress.

Author Summary

Here we use molecular dynamics simulation to compare the mechanical properties of the 7-ankyrin-repeat oncoprotein Gankyrin in isolation and in complex with binding partner S6-C. Tandem repeat proteins like Gankyrin comprise tandem arrays of small structural motifs that pack linearly to produce elongated architectures. They are elastic, mechanically weak molecules and they unfold and refold repeat by repeat under force. We show that S6-C binding greatly increases the resistance of Gankyrin to mechanical stress. The enhanced mechanical stability is specific to those ankyrin repeats in contact with S6-C, and the localized nature of the effect results in fundamental changes in the way the protein responds to force. Thus, the forced unfolding of isolated Gankryin involves a diverse set of pathways with a preference for a C- to N-terminus unfolding mechanism whereas this diversity is reduced upon complex formation with the central repeats, which are those most tightly bound to the ligand, tending to unfold last. Our study shows how stepwise unfolding can buffer repeat proteins and their binding partners from mechanical stress in the cell. It also points to a mechano-switching mechanism whereby binding between two partner macromolecules is regulated by mechanical stress.

One of the applications of Molecular Dynamics (MD) simulations is to explore the energetic barriers to mechanical unfolding of proteins such as occurs in response to the mechanical pulling of single molecules in Atomic Force Microscopy (AFM) experiments. Although Steered Molecular Dynamics simulations have provided microscopic details of the unfolding process during the pulling, the simulated forces required for unfolding are typically far in excess of the measured values. To rectify this, we have developed the Pulsed Unconstrained Fluctuating Forces (PUFF) method, which induces constant-momentum motions by applying forces directly to the instantaneous velocity of selected atoms in a protein system. The driving forces are applied in pulses, which allows the system to relax between pulses, resulting in more accurate unfolding force estimations than in previous methods. In the cases of titin, ubiquitin and e2lip3, the PUFF trajectories produce force fluctuations that agree quantitatively with AFM experiments. Another useful property of PUFF is that simulations get trapped if the target momentum is too low, simplifying the discovery and analysis of unfolding intermediates.

We report high temperature molecular dynamics simulations of the unfolding of the TRPZ1 peptide using an explicit model for the solvent. The system has been simulated for a total of 6 μs with 100-ns minimal continuous stretches of trajectory. The populated states along the simulations are identified by monitoring multiple observables, probing both the structure and the flexibility of the conformations. Several unfolding and refolding transition pathways are sampled and analyzed. The unfolding process of the peptide occurs in two steps because of the accumulation of a metastable on-pathway intermediate state stabilized by two native backbone hydrogen bonds assisted by nonnative hydrophobic interactions between the tryptophan side chains. Analysis of the un/folding kinetics and classical commitment probability calculations on the conformations extracted from the transition pathways show that the rate-limiting step for unfolding is the disruption of the ordered native hydrophobic packing (Trp-zip motif) leading from the native to the intermediate state. But, the speed of the folding process is mainly determined by the transition from the completely unfolded state to the intermediate and specifically by the closure of the hairpin loop driven by formation of two native backbone hydrogen bonds and hydrophobic contacts between tryptophan residues. The temperature dependence of the unfolding time provides an estimate of the unfolding activation enthalpy that is in agreement with experiments. The unfolding time extrapolated to room temperature is in agreement with the experimental data as well, thus providing a further validation to the analysis reported here.

Amyloid -protein (A) is central to the pathology of Alzheimer's disease. A 5% difference in the primary structure of the two predominant alloforms, A and A, results in distinct assembly pathways and toxicity properties. Discrete molecular dynamics (DMD) studies of A and A assembly resulted in alloform-specific oligomer size distributions consistent with experimental findings. Here, a large ensemble of DMD–derived A and A monomers and dimers was subjected to fully atomistic molecular dynamics (MD) simulations using the OPLS-AA force field combined with two water models, SPCE and TIP3P. The resulting all-atom conformations were slightly larger, less compact, had similar turn and lower -strand propensities than those predicted by DMD. Fully atomistic A and A monomers populated qualitatively similar free energy landscapes. In contrast, the free energy landscape of A dimers indicated a larger conformational variability in comparison to that of A dimers. A dimers were characterized by an increased flexibility in the N-terminal region D1-R5 and a larger solvent exposure of charged amino acids relative to A dimers. Of the three positively charged amino acids, R5 was the most and K16 the least involved in salt bridge formation. This result was independent of the water model, alloform, and assembly state. Overall, salt bridge propensities increased upon dimer formation. An exception was the salt bridge propensity of K28, which decreased upon formation of A dimers and was significantly lower than in A dimers. The potential relevance of the three positively charged amino acids in mediating the A oligomer toxicity is discussed in the light of available experimental data.

When analyzing single-molecule data, a low-dimensional set of system observables typically serve as the observational data. We calibrate stochastic dynamical models from time series that record such observables (our focus throughout is on a molecule’s end-to-end distance). Numerical techniques for quantifying noise from multiple time scales in a single trajectory, including experimental instrument and inherent thermal noise, are demonstrated. The techniques are applied to study time series coming from both simulations and experiments associated with the nonequilibrium mechanical unfolding of titin’s I27 domain. The estimated models can be used for several purposes: (1) detect dynamical signatures of “rare events” by analyzing the effective diffusion and force as a function of the monitored observable; (2) quantify the influence that experimentally unobservable conformational degrees of freedom have on the dynamics of the monitored observable; (3) quantitatively compare the inherent thermal noise to other noise sources, e.g. instrument noise, variation induced by conformational heterogeneity, etc.; (4) simulate random quantities associated with repeated experiments; (5) apply pathwise (i.e. trajectory-wise) hypothesis tests to assess the goodness-of-fit of models and even detect conformational transitions in noisy signals. These items are all illustrated with several examples.

We use Langevin dynamics to investigate the role played by the recently discovered force-induced entropic energy barrier on the two-state hopping phenomena that has been observed in single RNA, DNA and protein molecules placed under a stretching force. Simple considerations about the free energy of a molecule readily show that the application of force introduces an entropic barrier separating the collapsed state of the molecule, from a force-driven extended conformation. A notable characteristic of the force induced barrier is its long distances to transition state, up to tens of nanometers, which renders the kinetics of crossing this barrier highly sensitive to an applied force. Langevin dynamics across such force induced barriers readily demonstrates the hopping behavior observed for a variety of single molecules placed under force. Such hopping is frequently interpreted as a manifestation of two-state folding/unfolding reactions observed in bulk experiments. However, given that such barriers do not exist at zero force these reactions do not take place at all in bulk.

The folding of the B-domain of staphylococcal protein A has been studied by coarse-grained canonical and multiplexed replica-exchange molecular dynamics simulations with the UNRES force field in a broad range of temperatures (270K ≤ T ≤ 350K). In canonical simulations, the folding was found to occur either directly to the native state or through kinetic traps, mainly the topological mirror image of the native three-helix bundle. The latter folding scenario was observed more frequently at low temperatures. With increase of temperature, the frequency of the transitions between the folded and misfolded/unfolded states increased and the folded state became more diffuse with conformations exhibiting increased root-mean-square deviations from the experimental structure (from about 4 Å at T = 300K to 8.7 Å at T = 325K). An analysis of the equilibrium conformational ensemble determined from multiplexed replica exchange simulations at the folding-transition temperature (Tf = 325K) showed that the conformational ensemble at this temperature is a collection of conformations with residual secondary structures, which possess native or near-native clusters of nonpolar residues in place, and not a 50%-50% mixture of fully-folded and fully-unfolded conformations. These findings contradict the quasi-chemical picture of two- or multi-state protein folding, which assumes an equilibrium between the folded, unfolded, and intermediate states, with equilibrium shifting with temperature but with the native conformations remaining essentially unchanged. Our results also suggest that long-range hydrophobic contacts are the essential factor to keep the structure of a protein thermally stable.

While several experimental techniques now exist for characterizing protein unfolded states, all-atom simulation of unfolded states has been challenging due to the long time scales and conformational sampling required. We address this problem by using a combination of accelerated calculations on graphics processor units and distributed computing to simulate tens of thousands of molecular dynamics trajectories each up to ~10 μs (for a total aggregate simulation time of 127 milliseconds). We used this approach in conjunction with Trp-Cys contact quenching experiments to characterize the unfolded structure and dynamics of protein L. We employed a polymer theory method to make quantitative comparisons between high temperature simulated and chemically denatured experimental ensembles, and find that reaction-limited quenching rates calculated from simulation agree remarkably well with experiment. In both experiment and simulation, we find that unfolded state intramolecular diffusion rates are very slow compared to highly denatured chains, and that a single-residue mutation can significantly alter unfolded state dynamics and structure. This work suggests a view of the unfolded state in which surprisingly low diffusion rates could limit folding, and opens the door for all-atom molecular simulation to be a useful predictive tool for characterizing protein unfolded states along with experiments that directly measure intramolecular diffusion.

The hemostatic function of von Willebrand factor is downregulated by the metalloprotease ADAMTS13, which cleaves at a unique site normally buried in the A2 domain. Exposure of the proteolytic site is induced in the wild-type by shear stress as von Willebrand factor circulates in blood. Mutations in the A2 domain, which increase its susceptibility to cleavage, cause type 2A von Willebrand disease. In this study, molecular dynamics simulations suggest that the A2 domain unfolds under tensile force progressively through a series of steps. The simulation results also indicated that three type 2A mutations in the C-terminal half of the A2 domain, L1657I, I1628T and E1638K, destabilize the native state fold of the protein. Furthermore, all three type 2A mutations lowered in silico the tensile force necessary to undock the C-terminal helix 6 from the rest of the A2 domain, the first event in the unfolding pathway. The mutations F1520A, I1651A and A1661G were also predicted by simulations to destabilize the A2 domain and facilitate exposure of the cleavage site. Recombinant A2 domain proteins were expressed and cleavage assays were performed with the wild-type and single-point mutants. All three type 2A and two of the three predicted mutations exhibited increased rate of cleavage by ADAMTS13. These results confirm that destabilization of the helix 6 in the A2 domain facilitates exposure of the cleavage site and increases the rate of cleavage by ADAMTS13.

The tumor suppressor p53 plays a crucial role in the cell cycle checkpoints, DNA repair, and apoptosis. p53 consists of a natively unfolded N-terminal region (NTR), central DNA binding domain (DBD), C-terminal tetramerization domain, and regulatory region. In this paper, the interactions between the DBD and the NTR, and between the DBD and DNA were investigated by measuring changes in the mechanical unfolding trajectory of the DBD using atomic force microscopy (AFM)-based single molecule force spectroscopy. In the absence of DNA, the DBD (94–293, 200 amino acids (AA)) showed two different mechanical unfolding patterns. One indicated the existence of an unfolding intermediate consisting of approximately 60 AA, and the other showed a 100 AA intermediate. The DBD with the NTR did not show such unfolding patterns, but heterogeneous unfolding force peaks were observed. Of the heterogeneous patterns, we observed a high frequency of force peaks indicating the unfolding of a domain consisting of 220 AA, which is apparently larger than that of a sole DBD. This observation implies that a part of NTR binds to the DBD, and the mechanical unfolding happens not solely on the DBD but also accompanying a part of NTR. When DNA is bound, the mechanical unfolding trajectory of p53NTR+DBD showed a different pattern from that without DNA. The pattern was similar to that of the DBD alone, but two consecutive unfolding force peaks corresponding to 60 and 100 AA sub-domains were observed. These results indicate that interactions with the NTR or DNA alter the mechanical stability of DBD and result in drastic changes in the mechanical unfolding trajectory of the DBD.

Cells severely depleted of filamins were observed to have numerous motility-related defects, including a defect in endoplasmic spreading; smaller, more dynamic focal adhesions; and an inability to sustain high levels of traction force. The endoplasm as a separate mechanical unit spread by pulling forces is also discussed.

Cell motility is an essential process that depends on a coherent, cross-linked actin cytoskeleton that physically coordinates the actions of numerous structural and signaling molecules. The actin cross-linking protein, filamin (Fln), has been implicated in the support of three-dimensional cortical actin networks capable of both maintaining cellular integrity and withstanding large forces. Although numerous studies have examined cells lacking one of the multiple Fln isoforms, compensatory mechanisms can mask novel phenotypes only observable by further Fln depletion. Indeed, shRNA-mediated knockdown of FlnA in FlnB–/– mouse embryonic fibroblasts (MEFs) causes a novel endoplasmic spreading deficiency as detected by endoplasmic reticulum markers. Microtubule (MT) extension rates are also decreased but not by peripheral actin flow, because this is also decreased in the Fln-depleted system. Additionally, Fln-depleted MEFs exhibit decreased adhesion stability that appears in increased ruffling of the cell edge, reduced adhesion size, transient traction forces, and decreased stress fibers. FlnA–/– MEFs, but not FlnB–/– MEFs, also show a moderate defect in endoplasm spreading, characterized by initial extension followed by abrupt retractions and stress fiber fracture. FlnA localizes to actin linkages surrounding the endoplasm, adhesions, and stress fibers. Thus we suggest that Flns have a major role in the maintenance of actin-based mechanical linkages that enable endoplasmic spreading and MT extension as well as sustained traction forces and mature focal adhesions.

The principle of microscopic reversibility states that at equilibrium the number of molecules entering a state by a given path must equal those exiting the state via the same path under identical conditions, or in structural terms, that the conformations along the two pathways are the same. There has been some indirect evidence indicating that protein folding is such a process, but there have been few conclusive findings. In this study, we performed molecular dynamics simulations of an ultra-fast unfolding and folding protein at its melting temperature in order to observe, on an atom-by-atom basis, the pathways the protein followed as it unfolded and folded within a continuous trajectory. In a total of 0.67 μs of simulation in water, we found 6 transient denaturing events near the melting temperature (323K and 330K) and an additional refolding event following a previously identified unfolding event at high-temperature (373K). In each case unfolding and refolding transition state ensembles were identified, and they agreed well with experiment based on comparison of S- and Φ-values. Based on several structural properties, these 13 transition state ensembles agreed very well with each other and with 4 previously identified transition states from high-temperature denaturing simulations. Thus, not only were the unfolding and refolding transition states part of the same ensemble, but in five of the seven cases, the pathway the protein took as it unfolded was nearly identical to the subsequent refolding pathway. These events provide compelling evidence that protein folding is a microscopically reversible process. In the other two cases, the folding and unfolding transition states were remarkably similar to each other, but the paths deviated.

Fibronectin polymerization is essential for the development and repair of the extracellular matrix. Consequently, deciphering the mechanism of fibronectin fibril formation is of immense interest. Fibronectin fibrillogenesis is driven by cell-traction forces that mechanically unfold particular modules within fibronectin. Previously, mechanical unfolding of fibronectin has been modeled by applying tensile forces at the N- and C-termini of fibronectin domains; however, physiological loading is likely focused on the solvent-exposed RGD loop in the 10th type-III repeat of fibronectin (10FNIII), which mediates binding to cell-surface integrin receptors. In this work we used steered molecular dynamics to study the mechanical unfolding of 10FNIII under tensile force applied at this RGD site. We demonstrate that mechanically unfolding 10FNIII by pulling at the RGD site requires less work than unfolding by pulling at the N- and C- termini. Moreover, pulling at the N- and C-termini leads to 10FNIII unfolding along several pathways while pulling on the RGD site leads to a single exclusive unfolding pathway that includes a partially unfolded intermediate with exposed hydrophobic N-terminal β-strands – residues that may facilitate fibronectin self-association. Additional mechanical unfolding triggers an essential arginine residue, which is required for high affinity binding to integrins, to move to a position far from the integrin binding site. This cell traction-induced conformational change may promote cell detachment after important partially unfolded kinetic intermediates are formed. These data suggest a novel mechanism that explains how cell-mediated forces promote fibronectin fibrillogenesis and how cell surface integrins detach from newly forming fibrils. This process enables cells to bind and unfold additional fibronectin modules – a method that propagates matrix assembly.