Similarity of force-induced unfolding of apomyoglobin to its chemical-induced unfolding: an atomistic molecular dynamics simulation approach

We have compared force-induced unfolding with traditional unfolding methods using apomyoglobin as a model protein. Using molecular dynamics simulation, we have investigated the structural stability as a function of the degree of mechanical perturbation. Both anisotropic perturbation by stretching two terminal atoms and isotropic perturbation by increasing the radius of gyration of the protein show the same key event of force-induced unfolding. Our primary results show that the native structure of apomyoglobin becomes destabilized against the mechanical perturbation as soon as the interhelical packing between the G and H helices is broken, suggesting that our simulation results share a common feature with the experimental observation that the interhelical contact is more important for the folding of apomyoglobin than the stability of individual helices. This finding is further confirmed by simulating both helix destabilizing and interhelical packing destabilizing mutants.     

Controlled unfolding and refolding of a single sodium-proton antiporter using atomic force microscopy

Single-molecule force-spectroscopy was employed to unfold and refold single sodium-proton antiporters (NhaA) of Escherichia coli from membrane patches. Although transmembrane alpha-helices and extracellular polypeptide loops exhibited sufficient stability to individually establish potential barriers against unfolding, two helices predominantly unfolded pairwise, thereby acting as one structural unit. Many of the potential barriers were detected unfolding NhaA either from the C-terminal or the N-terminal end. It was found that some molecular interactions stabilizing secondary structural elements were directional, while others were not. Additionally, some interactions appeared to occur between the secondary structural elements. After unfolding ten of the 12 helices, the extracted polypeptide was allowed to refold back into the membrane. After five seconds, the refolded polypeptide established all secondary structure elements of the native protein. One helical pair showed a characteristic spring like "snap in" into its folded conformation, while the refolding process of other helices was not detected in particular. Additionally, individual helices required characteristic periods of time to fold. Correlating these results with the primary structure of NhaA allowed us to obtain the first insights into how potential barriers establish and determine the folding kinetics of the secondary structure elements.     

Role of water on unfolding kinetics of helical peptides studied by molecular dynamics simulations

Molecular dynamics simulations have been carried out with four polypeptides, Ala13, Val(13), Ser13, and Ala4Gly5Ala4, in vacuo and with explicit hydration. The unfolding of the polypeptides, which are initially fully alpha-helix in conformation, has been monitored during trajectories of 0.3 ns at 350 K. A rank of Ala < Val < Ser < Gly is found in the order of increasing rate of unwinding. The unfolding of Ala13 and Val(13) is completed in hundreds of picoseconds, while that of Ser13 is about one order of magnitude faster. The helix content of the peptide containing glycine residues falls to zero within a few picoseconds. Ramachandran plots indicate quite distinct equilibrium distributions and time evolution of dihedral angles in water and in vacuum for each residue type. The unfolding of polyalanine and polyvaline helices is accelerated due to solvation. In contrast, polyserine is more stable in water compared to vacuum, because its side chains can form intramolecular hydrogen bonds with the backbone more readily in vacuum, which disrupts the helix. Distribution functions of the spatial and angular position of water molecules in the proximity of the polypeptide backbone polar groups reveal the stabilization of the coiled structures by hydration. The transition from helix to coil is characterized by the appearance of a new peak in the probability distribution at a specific location characteristic of hydrogen bond formation between water and backbone polar groups. No significant insertion of water molecules is observed at the precise onset of unwinding, while (i, i+3) hydrogen bond formation is frequently detected at the initiation of alpha-helix unwinding.     

Understanding beta-hairpin formation by molecular dynamics simulations of unfolding

We have studied the mechanism of formation of a 16-residue beta-hairpin from the protein GB1 using molecular dynamics simulations in an aqueous environment. The analysis of unfolding trajectories at high temperatures suggests a refolding pathway consisting of several transient intermediates. The changes in the interaction energies of residues are related with the structural changes during the unfolding of the hairpin. The electrostatic energies of the residues in the turn region are found to be responsible for the transition between the folded state and the hydrophobic core state. The van der Waals interaction energies of the residues in the hydrophobic core reflect the behavior of the radius of gyration of the core region. We have examined the opposing influences of the protein-protein (PP) energy, which favors the native state, and the protein-solvent (PS) energy, which favors unfolding, in the formation of the beta-hairpin structure. It is found that the behavior of the electrostatic components of PP and PS energies reflects the structural changes associated with the loss of backbone hydrogen bonding. Relative changes in the PP and PS van der Waals interactions are related with the disruption of the hydrophobic core of a protein. The results of the simulations support the hydrophobic collapse mechanism of beta-hairpin folding.     

Folding-unfolding of goat alpha-lactalbumin studied by stopped-flow circular dichroism and molecular dynamics simulations

Folding reaction of goat alpha-lactalbumin has been studied by stopped-flow circular dichroism and molecular dynamics simulations. The effects of four single mutations and a double mutation on the stability of the protein under a native condition were studied. The mutations were introduced into residues located at a hydrophobic core in the alpha-domain of the molecule. Here we show that an amino acid substitution (T29I) increases the native-state stability of goat alpha-lactalbumin against the guanidine hydrochloride-induced unfolding by 3.5 kcal/mol. Kinetic refolding and unfolding of wild-type and mutant goat alpha-lactalbumin measured by stopped-flow circular dichroism showed that the local structure around the Thr29 side chain was not constructed in the transition state of the folding reaction. To characterize the local structural change around the Thr29 side chain to an atomic level of resolution, we performed high-temperature (at 400 K and 600 K) molecular dynamics simulations and studied the structural change at an initial stage of unfolding observed in the simulation trajectories. The Thr29 portion of the molecule experienced structural disruption accompanied with the loss of inter-residue contacts and with the water molecule penetration in the 400-K simulation as well as in four of the six 600-K simulations. Disruption of the N-terminal portion was also observed and was consistent with the results of kinetic refolding/unfolding experiments shown in our previous report.     

Probing the folding and unfolding dynamics of secondary and tertiary structures in a three-helix bundle protein

Fast relaxation kinetics studies of the B-domain of staphylococcal protein A were performed to characterize the folding and unfolding of this small three-helix bundle protein. The relaxation kinetics were initiated using a laser-induced temperature jump and probed using time-resolved infrared spectroscopy. The kinetics monitored within the amide I' absorbance of the polypeptide backbone exhibit two distinct kinetics phases with nanosecond and microsecond relaxation times. The fast kinetics relaxation time is close to the diffusion limits placed on protein folding reactions. The fast kinetics phase is dominated by the relaxation of the solvated helix (nu = 1632 cm(-1)), which reports on the fast relaxation of the individual helices. The slow kinetics phase follows the cooperative relaxation of the native helical bundle core that is monitored by both solvated (nu = 1632 cm(-1)) and buried helical IR bands (nu = 1652 cm(-1)). The folding rates of the slow kinetics phase calculated over an extended temperature range indicate that the core formation of this protein follows a pathway that is energetically downhill. The unfolding rates are much more strongly temperature-dependent indicating an activated process with a large energy barrier. These results provide significant insight into the primary process of protein folding and suggest that fast formation of helices can drive the folding of helical proteins.     

Differential stability of beta-sheets and alpha-helices in beta-lactamase: a high temperature molecular dynamics study of unfolding intermediates.

beta-Lactamase, which catalyzes beta-lactam antibiotics, is prototypical of large alpha/beta proteins with a scaffolding formed by strong noncovalent interactions. Experimentally, the enzyme is well characterized, and intermediates that are slightly less compact and having nearly the same content of secondary structure have been identified in the folding pathway. In the present study, high temperature molecular dynamics simulations have been carried out on the native enzyme in solution. Analysis of these results in terms of root mean square fluctuations in cartesian and [phi, psi] space, backbone dihedral angles and secondary structural hydrogen bonds forms the basis for an investigation of the topology of partially unfolded states of beta-lactamase. A differential stability has been observed for alpha-helices and beta-sheets upon thermal denaturation to putative unfolding intermediates. These observations contribute to an understanding of the folding/unfolding processes of beta-lactamases in particular, and other alpha/beta proteins in general.     

Computer simulations of protein folding by targeted molecular dynamics

We have performed 128 folding and 45 unfolding molecular dynamics runs of chymotrypsin inhibitor 2 (CI2) with an implicit solvation model for a total simulation time of 0.4 microseconds. Folding requires that the three-dimensional structure of the native state is known. It was simulated at 300 K by supplementing the force field with a harmonic restraint which acts on the root-mean-square deviation and allows to decrease the distance to the target conformation. High temperature and/or the harmonic restraint were used to induce unfolding. Of the 62 folding simulations started from random conformations, 31 reached the native structure, while the success rate was 83% for the 66 trajectories which began from conformations unfolded by high-temperature dynamics. A funnel-like energy landscape is observed for unfolding at 475 K, while the unfolding runs at 300 K and 375 K as well as most of the folding trajectories have an almost flat energy landscape for conformations with less than about 50% of native contacts formed. The sequence of events, i.e., secondary and tertiary structure formation, is similar in all folding and unfolding simulations, despite the diversity of the pathways. Previous unfolding simulations of CI2 performed with different force fields showed a similar sequence of events. These results suggest that the topology of the native state plays an important role in the folding process.     

Thermal unfolding simulations of apo-calmodulin using leap-dynamics

The simulation method leap-dynamics (LD) has been applied to protein thermal unfolding simulations to investigate domain-specific unfolding behavior. Thermal unfolding simulations of the 148-residue protein apo-calmodulin with implicit solvent were performed at temperatures 290 K, 325 K, and 360 K and compared with the corresponding molecular dynamics trajectories in terms of a number of calculated conformational parameters. The main experimental results of unfolding are reproduced in showing the lower stability of the C-domain: at 290 K, both the N- and C-domains are essentially stable; at 325 K, the C-domain unfolds, whereas the N-domain remains folded; and at 360 K, both domains unfold extensively. This behavior could not be reproduced by molecular dynamics simulations alone under the same conditions. These results show an encouraging degree of convergence between experiment and LD simulation. The simulations are able to describe the overall plasticity of the apo-calmodulin structure and to reveal details such as reversible folding/unfolding events within single helices. The results show that by using the combined application of a fast and efficient sampling routine with a detailed molecular dynamics force field, unfolding simulations of proteins at atomic resolution are within the scope of current computational power.     

Unfolding a linker between helical repeats

In many multi-repeat proteins, linkers between repeats have little secondary structure and place few constraints on folding or unfolding. However, the large family of spectrin-like proteins, including alpha-actinin, spectrin, and dystrophin, share three-helix bundle, spectrin repeats that appear in crystal structures to be linked by long helices. All of these proteins are regularly subjected to mechanical stress. Recent single molecule atomic force microscopy (AFM) experiments demonstrate not only forced unfolding but also simultaneous unfolding of tandem repeats at finite frequency, which suggests that the contiguous helix between spectrin repeats can propagate a cooperative helix-to-coil transition. Here, we address what happens atomistically to the linker under stress by steered molecular dynamics simulations of tandem spectrin repeats in explicit water. The results for alpha-actinin repeats reveal rate-dependent pathways, with one pathway showing that the linker between repeats unfolds, which may explain the single-repeat unfolding pathway observed in AFM experiments. A second pathway preserves the structural integrity of the linker, which explains the tandem-repeat unfolding event. Unfolding of the linker begins with a splay distortion of proximal loops away from hydrophobic contacts with the linker. This is followed by linker destabilization and unwinding with increased hydration of the backbone. The end result is an unfolded helix that mechanically decouples tandem repeats. Molecularly detailed insights obtained here aid in understanding the mechanical coupling of domain stability in spectrin family proteins.     

Helix nucleation kinetics from molecular simulations in explicit solvent

We study the reversible folding/unfolding of short Ala and Gly-based peptides by molecular dynamics simulations of all-atom models in explicit water solvent. A kinetic analysis shows that the formation of a first alpha-helical turn occurs within 0.1-1 ns, in agreement with the analyses of laser temperature jump experiments. The unfolding times exhibit Arrhenius temperature dependence. For a rapidly nucleating all-Ala peptide, the helix nucleation time depends only weakly on temperature. For a peptide with enthalpically competing turn-like structures, helix nucleation exhibits an Arrhenius temperature dependence, corresponding to the unfolding of enthalpic traps in the coil ensemble. An analysis of structures in a "transition-state ensemble" shows that helix-to-coil transitions occur predominantly through breaking of hydrogen bonds at the helix ends, particularly at the C-terminus. The temperature dependence of the transition-state ensemble and the corresponding folding/unfolding pathways illustrate that folding mechanisms can change with temperature, possibly complicating the interpretation of high-temperature unfolding simulations. The timescale of helix formation is an essential factor in molecular models of protein folding. The rapid helix nucleation observed here suggests that transient helices form early in the folding event.     

Probing the energy landscape of the membrane protein bacteriorhodopsin

The folding and stability of transmembrane proteins is a fundamental and unsolved biological problem. Here, single bacteriorhodopsin molecules were mechanically unfolded from native purple membranes using atomic force microscopy and force spectroscopy. The energy landscape of individual transmembrane alpha helices and polypeptide loops was mapped by monitoring the pulling speed dependence of the unfolding forces and applying Monte Carlo simulations. Single helices formed independently stable units stabilized by a single potential barrier. Mechanical unfolding of the helices was triggered by 3.9-7.7 A extension, while natural unfolding rates were of the order of 10(-3) s(-1). Besides acting as individually stable units, helices associated pairwise, establishing a collective potential barrier. The unfolding pathways of individual proteins reflect distinct pulling speed-dependent unfolding routes in their energy landscapes. These observations support the two-stage model of membrane protein folding in which alpha helices insert into the membrane as stable units and then assemble into the functional protein.     

Thermal‐ and urea‐induced unfolding processes of glutathione S‐transferase by molecular dynamics simulation

The Schistosoma juponicum 26 kDa glutathione S‐transferase (sj26GST) consists of the N‐terminal domain (N‐domain), containing three alpha‐helices (named H1‐H3) and four anti‐parallel beta‐strands (S1‐S4), and the C‐terminal domain (C‐domain), comprising five alpha‐helices (named H4‐H8). In present work, molecular dynamics simulations and fluorescence spectroscopic were used to gain insights into the unfolding process of sj26GST. The molecular dynamics simulations on sj26GST subunit both in water and in 8 M urea were carried out at 300 K, 400 K and 500 K, respectively. Spectroscopic measurements were employed to monitor structural changes. Molecular dynamics simulations of sj26GST subunit induced by urea and temperature showed that the initial unfolding step of sj26GST both in water and urea occurred on N‐domain, involving the disruption of helices H2, H3 and strands S3 and S4, whereas H6 was the last region exposed to solution and was the last helix to unfold. Moreover, simulations analyses combining with fluorescence and circular dichroism spectra indicated that N‐domain could not fold independent, suggesting that correct folding of N‐domain depended on its interactions with C‐domain. We further proposed that the folding of GSTs could begin with the hydrophobic collapse of C‐domain whose H4, H5, H6 and H7 could move close to each other and form a hydrophobic core, especially H6 wrapped in the hydrophobic center and beginning spontaneous formation of the helix. S3, S4, H3, and H2 could form in the wake of the interaction between C‐domain and N‐domain. The paper can offer insights into the molecular mechanism of GSTs unfolding. © 2014 Wiley Periodicals, Inc. Biopolymers 103: 247–259, 2015.     

Rates of unfolding, rather than refolding, determine thermal stabilities of thermophilic, mesophilic, and psychrotrophic 3-isopropylmalate dehydrogenases

The relationship between the thermal stability of proteins and rates of unfolding and refolding is still an open issue. The data are very scarce, especially for proteins with complex structure. Here, time-dependent denaturation-renaturation experiments on Thermus thermophilus, Escherichia coli, and Vibrio sp. I5 3-isopropylmalate dehydrogenases (IPMDHs) of different heat stabilities are presented. Unfolding, as monitored by several methods, occurs in a single first-order step with half-times of approximately 1 h, several minutes, and few seconds for the thermophilic, mesophilic, and psychrotrophic enzymes, respectively. The binding of Mn*IPM (the manganese complex of 3-isopropylmalate) markedly reduces the rates of unfolding; this effect is more prominent for the less stable enzyme variants. Refolding is a two-step or multistep first-order process involving an inactive intermediate(s). The restoration of the native structure and reactivation take place with a half-time of a few minutes for all three IPMDHs. Thus, the comparative experimental unfolding-refolding studies of the three IPMDHs with different thermostabilities have revealed a close relationship between thermostability and unfolding rate. Structural analysis has shown that the differences in the molecular contacts between selected nonconserved residues are responsible for the different rates of unfolding. On the other hand, the folding rates might be correlated with the absolute contact order, which does not significantly vary between IPMDHs with different thermostabilities. On the basis of our observations, folding rates appear to be dictated by global structural characteristics (such as native topology, i.e., contact order) rather than by thermodynamic stability.     

Mechanical unfolding of a beta-hairpin using molecular dynamics

Single-molecule mechanical unfolding experiments have the potential to provide insights into the details of protein folding pathways. To investigate the relationship between force-extension unfolding curves and microscopic events, we performed molecular dynamics simulations of the mechanical unfolding of the C-terminal hairpin of protein G. We have studied the dependence of the unfolding pathway on pulling speed, cantilever stiffness, and attachment points. Under conditions that generate low forces, the unfolding trajectory mimics the untethered, thermally accessible pathway previously proposed based on high-temperature studies. In this stepwise pathway, complete breakdown of backbone hydrogen bonds precedes dissociation of the hydrophobic cluster. Under more extreme conditions, the cluster and hydrogen bonds break simultaneously. Transitions between folding intermediates can be identified in our simulations as features of the calculated force-extension curves.     

Unfolding of hen egg lysozyme by molecular dynamics simulations at 300K: insight into the role of the interdomain interface

We present the results of two 1.2 ns molecular dynamics (MD) unfolding simulations on hen egg lysozyme in water at 300K, performed using a new procedure called PEDC (Path Exploration With Distance Constraints). This procedure allows exploration of low energy structures as a function of increasing RMSD from the native structure, and offers especially the possibility of extensive exploration of the conformational space during the initial unfolding stages. The two independent MD simulations gave similar chronology of unfolding events: disruption of the active site, kinking of helix C, partial unfolding of the three-stranded beta-sheet to a two-stranded sheet (during which the helices A, B, and D remain to a great extent native), and finally unfolding of the beta-domain and partial unfolding of the alpha-domain in which hydrophobic clusters persist. We show particularly that the loss of hydrophobic contacts between the beta-sheet turn residues Leu55 and Ile56 and the hydrobic patch of the alpha-domain destabilizes the beta-domain and leads to its unfolding, suggesting that the correct embedding of these residues in the alpha-beta interface may constitute the rate limiting step in folding. These results are in accord with experimental observations on the folding/unfolding behavior of hen egg lysozyme at room temperature. They would also explain the loss of stability and the tendency to aggregation observed for the mutant Leu55Thr, and the slow refolding kinetics observed in the analogous amyloidogenic variant of human lysozyme.     

Structural stability of disulfide mutants of basic pancreatic trypsin inhibitor: A molecular dynamics study

The structure and folding of basic pancreatic trypsin inhibitor (BPTI) has been studied extensively by experimental means. We report a computer simulation study of the structural stability of various disulfide mutants of BPTI, involving eight 250-psec molecular dynamics simulations of the proteins in water, with and without a phosphate counterion. The presence of the latter alters the relative stability of the single disulfide species [5–55] and [30–51]. This conclusion can explain results of mutational studies and the conservation of residues in homologues of BPTI, and suggests a possible role of ions in stabilizing one intermediate over another in unfolding or folding processes. © 1996 Wiley-Liss, Inc.     

Mass-weighted molecular dynamics simulation and conformational analysis of polypeptide.

Atomic motions in protein molecules have been studied by molecular dynamics (MD) simulations; dynamics simulation methods have also been employed in conformational studies of polypeptide molecules. It was found that when atomic masses are weighted, the molecular dynamics method can significantly increase the sampling of dihedral conformation space in such studies, compared to a conventional MD simulation of the same total simulation time length. Herein the theoretical study of molecular conformation sampling by the molecular dynamics-based simulation method in which atomic masses are weighted is reported in detail; moreover, a numerical scheme for analyzing the extensive conformational sampling in the simulation of a tetrapeptide amide molecule is presented. From numerical analyses of the mass-weighted molecular dynamics trajectories of backbone dihedral angles, low-resolution structures covering the entire backbone dihedral conformation space of the molecule were determined, and the distribution of rotationally stable conformations in this space were analyzed quantitatively. The theoretical analyses based on the computer simulation and numerical analytical methods suggest that distinctive regimes in the conformational space of the peptide molecule can be identified.