Mechanical unfolding of a titin Ig domain: structure of transition state revealed by combining atomic force microscopy,protein engineering and molecular dynamics simulations

Titin I27 shows a high resistance to unfolding when subject to external force. To investigate the molecular basis of this mechanical stability, protein engineering Phi-value analysis has been combined with atomic force microscopy to investigate the structure of the barrier to forced unfolding. The results indicate that the transition state for forced unfolding is significantly structured, since highly destabilising mutations in the core do not affect the force required to unfold the protein. As has been shown before, mechanical strength lies in the region of the A' and G-strands but, contrary to previous suggestions, the results indicate clearly that side-chain interactions play a significant role in maintaining mechanical stability. Since Phi-values calculated from molecular dynamics simulations are the same as those determined experimentally, we can, with confidence, use the molecular dynamics simulations to analyse the structure of the transition state in detail, and are able to show loss of interactions between the A' and G-strands with associated A-B and E-F loops in the transition state. The key event is not a simple case of loss of hydrogen bonding interactions between the A' and G-strands alone. Comparison with Phi-values from traditional folding studies shows differences between the force and "no-force" transition states but, nevertheless, the region important for kinetic stability is the same in both cases. This explains the correspondence between hierarchy of kinetic stability (measured in stopped-flow denaturant studies) and mechanical strength in these titin domains.     

Mechanical unfolding of a titin Ig domain: structure of unfolding intermediate revealed by combining AFM,molecular dynamics simulations,NMR and protein engineering

The mechanical unfolding of an immunoglobulin domain from the human muscle protein titin (TI I27) has been shown to proceed via a metastable intermediate in which the A-strand is detached. The structure and properties of this intermediate are characterised in this study. A conservative destabilising mutation in the A-strand has no effect on the unfolding force, nor the dependence of the unfolding force on the pulling speed, indicating that the unfolding forces measured in an AFM experiment are those required for the unfolding of the intermediate and not the native state. A mutant of TI I27 with the A-strand deleted (TI I27-A) is studied by NMR and standard biophysical techniques, combined with protein engineering. Molecular dynamics simulations show TI I27-A to be a good model for the intermediate. It has a structure very similar to the native state, and is surprisingly stable. Comparison with a Phi-value analysis of the unfolding pathway clearly shows that the protein unfolds by a different pathway under an applied force than on addition of denaturant.     

Folding/Unfolding Kinetics of Apomyoglobin

The equilibrium and kinetic folding/unfolding of apomyoglobin (ApoMb) were studied at pH 6.2, 11 °C by recording tryptophan fluorescence. The equilibrium unfolding of ApoMb in the presence of urea was shown to involve accumulation of an intermediate state, which had a higher fluorescence intensity as compared with the native and unfolded states. The folding proceeded through two kinetic phases, a rapid transition from the unfolded to the intermediate state and a slow transition from the intermediate to the native state. The accumulation of the kinetic intermediate state was observed in a wide range of urea concentrations. The intermediate was detected even in the region corresponding to the unfolding limb of the chevron plot. Urea concentration dependence was obtained for the observed folding/unfolding rate. The shape of the dependence was compared with that of two-state proteins characterized by a direct transition from the unfolded to the native state.     

Complex stability of single proteins explored by forced unfolding experiments

In the last decade atomic force microscopy has been used to measure the mechanical stability of single proteins. These force spectroscopy experiments have shown that many water-soluble and membrane proteins unfold via one or more intermediates. Recently, Li and co-workers found a linear correlation between the unfolding force of the native state and the intermediate in fibronectin, which they suggested indicated the presence of a molecular memory or multiple unfolding pathways (1). Here, we apply two independent methods in combination with Monte Carlo simulations to analyze the unfolding of alpha-helices E and D of bacteriorhodopsin (BR). We show that correlation analysis of unfolding forces is very sensitive to errors in force calibration of the instrument. In contrast, a comparison of relative forces provides a robust measure for the stability of unfolding intermediates. The proposed approach detects three energetically different states of alpha-helices E and D in trimeric BR. These states are not observed for monomeric BR and indicate that substantial information is hidden in forced unfolding experiments of single proteins.     

Structural details,pathways, and energetics of unfolding apomyoglobin

Protein folding is often difficult to characterize experimentally because of the transience of intermediate states, and the complexity of the protein-solvent system. Atomistic simulations, which could provide more detailed information, have had to employ highly simplified models or high temperatures, to cope with the long time scales of unfolding; direct simulation of folding is even more problematic. We report a fully atomistic simulation of the acid-induced unfolding of apomyoglobin in which the protonation of acidic side-chains to simulate low pH is sufficient to induce unfolding at room temperature with no added biasing forces or other unusual conditions; and the trajectory is validated by comparison to experimental characterization of intermediate states. Novel insights provided by their analysis include: characterization of a dry swollen globule state forming a barrier to initial unfolding or final folding; observation of cooperativity in secondary and tertiary structure formation and its explanation in terms of dielectric environments; and structural details of the intermediate and the completely unfolded states. These insights involve time scales and levels of structural detail that are presently beyond the range of experiment, but come within reach through the simulation methods described here. An implicit solvation model is used to analyze the energetics of protein folding at various pH and ionic strength values, and a reasonable estimate of folding free energy is obtained. Electrostatic interactions are found to disfavor folding.     

Complex Stability of Single Proteins Explored by Forced Unfolding Experiments

In the last decade atomic force microscopy has been used to measure the mechanical stability of single proteins. These force spectroscopy experiments have shown that many water-soluble and membrane proteins unfold via one or more intermediates. Recently, Li and co-workers found a linear correlation between the unfolding force of the native state and the intermediate in fibronectin, which they suggested indicated the presence of a molecular memory or multiple unfolding pathways (1). Here, we apply two independent methods in combination with Monte Carlo simulations to analyze the unfolding of α-helices E and D of bacteriorhodopsin (BR). We show that correlation analysis of unfolding forces is very sensitive to errors in force calibration of the instrument. In contrast, a comparison of relative forces provides a robust measure for the stability of unfolding intermediates. The proposed approach detects three energetically different states of α-helices E and D in trimeric BR. These states are not observed for monomeric BR and indicate that substantial information is hidden in forced unfolding experiments of single proteins.     

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.     

Mechanical unfolding of TNfn3: the unfolding pathway of a fnIII domain probed by protein engineering, AFM and MD simulation

Protein engineering Phi-value analysis combined with single molecule atomic force microscopy (AFM) was used to probe the molecular basis for the mechanical stability of TNfn3, the third fibronectin type III domain from human tenascin. This approach has been adopted previously to solve the forced unfolding pathway of a titin immunoglobulin domain, TI I27. TNfn3 and TI I27 are members of different protein superfamilies and have no sequence identity but they have the same beta-sandwich structure consisting of two antiparallel beta-sheets. TNfn3, however, unfolds at significantly lower forces than TI I27. We compare the response of these proteins to mechanical force. Mutational analysis shows that, as is the case with TI I27, TNfn3 unfolds via a force-stabilised intermediate. The key event in forced unfolding in TI I27 is largely the breaking of hydrogen bonds and hydrophobic interactions between the A' and G-strands. The mechanical Phi-value analysis and molecular dynamics simulations reported here reveal that significantly more of the TNfn3 molecule contributes to its resistance to force. Both AFM experimental data and molecular dynamics simulations suggest that the rate-limiting step of TNfn3 forced unfolding reflects a transition from the extended early intermediate to an aligned intermediate state. As well as losses of interactions of the A and G-strands and associated loops there are rearrangements throughout the core. As was the case for TI I27, the forced unfolding pathway of TNfn3 is different from that observed in denaturant studies in the absence of force.     

The remarkable mechanical strength of polycystin-1 supports a direct role in mechanotransduction

Polycystin-1 is a large membrane-associated protein that interacts with polycystin-2 in the primary cilia of renal epithelial cells to form a mechanosensitive ion channel. Bending of the cilia induces calcium flow into the cells, mediated by the polycystin complex. Antibodies to polycystin-1 and polycystin-2 abolish this activation. Based on this, it has been suggested that the extracellular region of polycystin-1, which has a number of putative binding domains, may act as a mechanosensor. A large proportion of the extracellular region of polycystin-1 consists of beta-sandwich PKD domains in tandem array. We use atomic force microscopy to investigate the mechanical properties of the PKD domains of polycystin-1. We show that these domains, despite having a low thermodynamic stability, exhibit a remarkable mechanical strength, similar to that of immunoglobulin domains in the giant muscle protein titin. In agreement with the experimental results molecular dynamics simulations performed at low constant force show that the first PKD domain of polycystin (PKDd1) has a similar unfolding time as titin I27, under the same conditions. The simulations suggest that the basis for this mechanical stability is the formation of a force-stabilised intermediate. Our results suggest that these domains will remain folded under external force supporting the hypothesis that polycystin-1 could act as a mechanosensor, detecting changes in fluid flow in the kidney tubule.     

Diffusing and colliding: the atomic level folding/unfolding pathway of a small helical protein

Proteins with ultra-fast folding/unfolding kinetics are excellent candidates for study by molecular dynamics. Here, we describe such simulations of a three helix bundle protein, the engrailed homeodomain (En-HD), which folds via the diffusion-collision model. The unfolding pathway of En-HD was characterized by seven simulations of the protein and 12 simulations of its helical fragments yielding over 1.1 micros of simulation time in water. Various conformational states along the unfolding pathway were identified. There is the compact native-like transition state, a U-shaped helical intermediate and an unfolded state with dynamic helical segments. Each of these states is in good agreement with experimental data. Examining these states as well as the transitions between them, we find the role of long-range tertiary contacts, specifically salt-bridges, important in the folding/unfolding pathway. In the folding direction, charged residues form long-range tertiary contacts before the hydrophobic core is formed. The formation of HII is assisted by a specific salt-bridge and by non-specific (fluctuating) tertiary contacts, which we call contact-assisted helix formation. Salt-bridges persist as the protein approaches the transition state, stabilizing HII until the hydrophobic core is formed. To complement this information, simulations of fragments of En-HD illustrate the helical propensities of the individual segments. By thermal denaturation, HII proved to be the least stable helix, unfolding in less than 450 ps at high temperature. We observed the low helical propensity of C-terminal residues from HIII in fragment simulations which, when compared to En-HD unfolding simulations, link the unraveling of HIII to the initial event that drives the unfolding of En-HD.     

嗜温冷休克蛋白力致去折叠研究

嗜温冷休克蛋白拥有一个由五个β股形成的反平行β桶结构,目前已被用于蛋白质去折叠的研究。当使用机械力对嗜温冷休克蛋白进行拉伸研究时,发现嗜温冷休克蛋白的去折叠过程具有明显的中间态。在常速和常力两种情况下对嗜温冷休克蛋白进行拉伸分子动力学模拟,发现其在两种情况下具有相同的去折叠次序,即C端β片层首先去折叠,随后N端β片层去折叠;同时这两种模拟都表现出明确的中间态。研究结果表明,嗜温冷休克蛋白抵抗外力作用除了依赖链间氢键外,分子内的静电相互作用也发挥着重要的作用。     

Direct observation of microscopic reversibility in single-molecule protein folding

Both folded and unfolded conformations should be observed for a protein at its melting temperature (T(m)), where DeltaG between these states is zero. In an all-atom molecular dynamics simulation of chymotrypsin inhibitor 2 (CI2) at its experimental T(m), 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 beta-sheet pulls apart and the alpha-helix unravels. The unfolded state reached under these conditions maintains a kernel of structure in the form of a non-native 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 to be the native state 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 T(m).     

Molecular basis of cooperativity in protein folding IV. Core: A general cooperative folding model

The cooperative nature of the protein folding process is independent of the characteristic fold and the specific secondary structure attributes of a globular protein. A general folding/unfolding model should, therefore, be based upon structural features that transcend the peculiarities of α-helices, β-sheets, and other structural motifs found in proteins. The studies presented in this paper suggest that a single structural characteristic common to all globular proteins is essential for cooperative folding. The formation of a partly folded state from the native state results in the exposure to solvent of two distinct regions: (1) the portions of the protein that are unfolded; and (2) the “complementary surfaces,” located in the regions of the protein that remain folded. The cooperative character of the folding/unfolding transition is determined largely by the energetics of exposing complementary surface regions to the solvent. By definition, complementary regions are present only in partly folded states; they are absent from the native and unfolded states. An unfavorable free energy lowers the probability of partly folded states and increases the cooperativity of the transition. In this paper we present a mathematical formulation of this behavior and develop a general cooperative folding/unfolding model, termed the “complementary region” (CORE) model. This model successfully reproduces the main properties of folding/unfolding transitions without limiting the number of partly folded states accessible to the protein, thereby permitting a systematic examination of the structural and solvent conditions under which intermediates become populated. It is shown that the CORE model predicts two-state folding/unfolding behavior, even though the two-state character is not assumed in the model. © 1993 Wiley-Liss, Inc.     

Observing the osmophobic effect in action at the single molecule level

Protecting osmolytes are widespread small organic molecules able to stabilize the folded state of most proteins against various denaturing stresses in vivo. The osmophobic model explains thermodynamically their action through a preferential exclusion of the osmolyte molecules from the protein surface, thus favoring the formation of intrapeptide hydrogen bonds. Few works addressed the influence of protecting osmolytes on the protein unfolding transition state and kinetics. Among those, previous single molecule force spectroscopy experiments evidenced a complexation of the protecting osmolyte molecules at the unfolding transition state of the protein, in apparent contradiction with the osmophobic nature of the protein backbone. We present single-molecule evidence that glycerol, which is a ubiquitous protecting osmolyte, stabilizes a globular protein against mechanical unfolding without binding into its unfolding transition state structure. We show experimentally that glycerol does not change the position of the unfolding transition state as projected onto the mechanical reaction coordinate. Moreover, we compute theoretically the projection of the unfolding transition state onto two other common reaction coordinates, that is, the number of native peptide bonds and the weighted number of native contacts. To that end, we augment an analytic Ising-like protein model with support for group-transfer free energies. Using this model, we find again that the position of the unfolding transition state does not change in the presence of glycerol, giving further support to the conclusions based on the single-molecule experiments.     

Structure of stable protein folding intermediates by equilibrium phi-analysis: the apoflavodoxin thermal intermediate

Protein intermediates in equilibrium with native states may play important roles in protein dynamics but, in cases, can initiate harmful aggregation events. Investigating equilibrium protein intermediates is thus important for understanding protein behaviour (useful or pernicious) but it is hampered by difficulties in gathering structural information. We show here that the phi-analysis techniques developed to investigate transition states of protein folding can be extended to determine low-resolution three-dimensional structures of protein equilibrium intermediates. The analysis proposed is based solely on equilibrium data and is illustrated by determination of the structure of the apoflavodoxin thermal unfolding intermediate. In this conformation, a large part of the protein remains close to natively folded, but a 40 residue region is clearly unfolded. This structure is fully consistent with the NMR data gathered on an apoflavodoxin mutant designed specifically to stabilise the intermediate. The structure shows that the folded region of the intermediate is much larger than the proton slow-exchange core at 25 degrees C. It also reveals that the unfolded region is made of elements whose packing surface is more polar than average. In addition, it constitutes a useful guide to rationally stabilise the native state relative to the intermediate state, a far from trivial task.     

Kinetic studies of folding of the B-domain of staphylococcal protein A with molecular dynamics and a united-residue (UNRES) model of polypeptide chains

Langevin dynamics is used with our physics-based united-residue (UNRES) force field to study the folding pathways of the B-domain of staphylococcal protein A (1BDD (alpha; 46 residues)). With 400 trajectories of protein A started from the extended state (to gather meaningful statistics), and simulated for more than 35 ns each, 380 of them folded to the native structure. The simulations were carried out at the optimal folding temperature of protein A with this force field. To the best of our knowledge, this is the first simulation study of protein-folding kinetics with a physics-based force field in which reliable statistics can be gathered. In all the simulations, the C-terminal alpha-helix forms first. The ensemble of the native basin has an average RMSD value of 4 A from the native structure. There is a stable intermediate along the folding pathway, in which the N-terminal alpha-helix is unfolded; this intermediate appears on the way to the native structure in less than one-fourth of the folding pathways, while the remaining ones proceed directly to the native state. Non-native structures persist until the end of the simulations, but the native-like structures dominate. To express the kinetics of protein A folding quantitatively, two observables were used: (i) the average alpha-helix content (averaged over all trajectories within a given time window); and (ii) the fraction of conformations (averaged over all trajectories within a given time window) with Calpha RMSD values from the native structure less than 5 A (fraction of completely folded structures). The alpha-helix content grows quickly with time, and its variation fits well to a single-exponential term, suggesting fast two-state kinetics. On the other hand, the fraction of folded structures changes more slowly with time and fits to a sum of two exponentials, in agreement with the appearance of the intermediate, found when analyzing the folding pathways. This observation demonstrates that different qualitative and quantitative conclusions about folding kinetics can be drawn depending on which observable is monitored.     

Three-state protein folding: experimental determination of free-energy profile

When considering protein folding with a transient intermediate, a difficulty arises as to determination of the rates of separate transitions. Here we overcome this problem, using the kinetic studies of the unfolding/refolding reactions of the three-state protein apomyoglobin as a model. Amplitudes of the protein refolding kinetic burst phase corresponding to the transition from the unfolded (U) to intermediate (I) state, that occurs prior to the native state (N) formation, allow us to estimate relative populations of the rapidly converting states at various final urea concentrations. On the basis of these proportions, a complicated experimental chevron plot has been deconvolved into the urea-dependent rates of the I<-->N and U<-->N transitions to give the dependence of free energies of the main transition state and of all three (N, I, and U) stable states on urea concentration.     

The folding pathway of spectrin R17 from experiment and simulation: using experimentally validated MD simulations to characterize States hinted at by experiment

We present an experimental and computational analysis of the folding pathway of the 17th domain of chicken brain alpha-spectrin, R17. Wild-type R17 folds in a two-state manner and the chevron plot (plot of the logarithm of the observed rate constant against concentration of urea) shows essentially linear folding and unfolding arms. A number of mutant proteins, however, show a change in slope of the unfolding arm at high concentration of denaturant, hinting at complexity in the folding landscape. Through a combination of mutational studies and high temperature molecular dynamics simulations we show that the folding of R17 can be described by a model with two sequential transition states separated by an intermediate species. The rate limiting transition state for folding in water has been characterized both through experimental Phi-value analysis and by simulation. In contrast, a detailed analysis of the transition state predicted to dominate under highly denaturing conditions is only possible by simulation.     

Detection and structure determination of an equilibrium unfolding intermediate of Rd-apocytochrome b562: native fold with non-native hydrophobic interactions

The absence of detectable kinetic and equilibrium folding intermediates by optical probes is commonly taken to indicate that protein folding is a two-state process. However, for some small proteins with apparent two-state behavior, unfolding intermediates have been identified in native-state hydrogen exchange or kinetic unfolding experiments monitored by nuclear magnetic resonance. Rd-apocytochrome b(562), a four-helix bundle, is one such protein. Here, we found another unfolding intermediate for Rd-apocytochrome b(562). It is based on a cooperative transition of (15)N chemical shifts of amide protons as a function of urea concentrations before the global unfolding. We have solved the high-resolution structure of the protein at 2.8 M urea, which is after this cooperative transition but before the global unfolding. All four helices remained intact, but a number of hydrophobic core residues repacked. This intermediate provides a possible structural interpretation for the kinetic unfolding intermediates observed using nuclear magnetic resonance methods for several proteins and has important implications for theoretical studies of protein folding.     

Thermal unfolding simulations of a multimeric protein--transition state and unfolding pathways

The folding of an oligomeric protein poses an extra challenge to the folding problem because the protein not only has to fold correctly; it has to avoid nonproductive aggregation. We have carried out over 100 molecular dynamics simulations using an implicit solvation model at different temperatures to study the unfolding of one of the smallest known tetramers, p53 tetramerization domain (p53tet). We found that unfolding started with disruption of the native tetrameric hydrophobic core. The transition state for the tetramer to dimer transition was characterized as a diverse ensemble of different structures using Phi value analysis in quantitative agreement with experimental data. Despite the diversity, the ensemble was still native-like with common features such as partially exposed tetramer hydrophobic core and shifts in the dimer-dimer arrangements. After passing the transition state, the secondary and tertiary structures continued to unfold until the primary dimers broke free. The free dimer had little secondary structure left and the final free monomers were random-coil like. Both the transition states and the unfolding pathways from these trajectories were very diverse, in agreement with the new view of protein folding. The multiple simulations showed that the folding of p53tet is a mixture of the framework and nucleation-condensation mechanisms and the folding is coupled to the complex formation. We have also calculated the entropy and effective energy for the different states along the unfolding pathway and found that the tetramerization is stabilized by hydrophobic interactions.