Unfolding of the cold shock protein studied with biased molecular dynamics

The cold shock protein from Bacillus caldolyticus is a small beta-barrel protein that folds in a two-state mechanism. For the native protein and for several mutants, a wealth of experimental data are available on stability and folding, so that it is an optimal system to study this process. We compare data from unfolding simulations (trajectories of 5 and up to 12 ns) obtained with a bias potential at room temperature and from unbiased thermal unfolding simulations with experimental data. The unfolding patterns derived from the trajectories starting from different native-like conformations and subject to different unfolding conditions agree. The transition state found in the simulations of unfolding is close to the native structure in agreement with experiment. Moreover, a lower value of the free energy barrier of unfolding was found for the mutant R3E than for the mutant E46A and the native protein, as indicated by experimental data. The first unfolding event involves the three-stranded beta-sheet whose decomposition corresponds to the transition state. In contrast to conclusions drawn from experiments, we found that the two-stranded beta-strand forms the most stable substructure, which decomposes very late in the unfolding process. However, assuming that this structure forms very early in the folding process, our findings would not contradict the experiments but require a different interpretation of them.     

Protein folding kinetics provides a context-independent assessment of beta-strand propensity in the Fyn SH3 domain

Structural database-derived propensities for amino acids to adopt particular local protein structures, such as alpha-helix and beta-strand, have long been recognized and effectively exploited for the prediction of protein secondary structure. However, the experimental verification of database-derived propensities using mutagenesis studies has been problematic, especially for beta-strand propensities, because local structural preferences are often confounded by non-local interactions arising from formation of the native tertiary structure. Thus, the overall thermodynamic stability of a protein is not always altered in a predictable manner by changes in local structural propensity at a single position. In this study, we have undertaken an investigation of the relationship between beta-strand propensity and protein folding kinetics. By characterizing the effects of a wide variety of amino acid substitutions at two different beta-strand positions in an SH3 domain, we have found that the observed changes in protein folding rates are very well correlated to beta-strand propensities for almost all of the substitutions examined. In contrast, there is little correlation between propensities and unfolding rates. These data indicate that beta-strand conformation is well formed in the structured portion of the SH3 domain transition state, and that local structure propensity strongly influences the stability of the transition state. Since the transition state is known to be packed more loosely than the native state and likely lacks many of the non-local stabilizing interactions seen in the native state, we suggest that folding kinetics studies may generally provide an effective means for the experimental validation of database-derived local structural propensities.     

Localized nature of the transition-state structure in goat alpha-lactalbumin folding

To investigate whether the structure partially formed in the molten globule folding intermediate of goat alpha-lactalbumin is further organized in the transition state of folding, we constructed a number of mutant proteins and performed Phi-value analysis on them. For this purpose, we measured the equilibrium unfolding transitions and kinetic refolding and unfolding reactions of the mutants using equilibrium and stopped-flow kinetic circular dichroism techniques. The results show that the mutants with mutations located in the A-helix (V8A, L12A), the B-helix (V27A), the beta-domain (L52A, W60A), the C-helix (K93A, L96A), the C-D loop (Y103F), the D-helix (L105A, L110A), and the C-terminal 3(10)-helix (W118F), have low Phi-values, less than 0.2. On the other hand, D87N, which is located on the Ca(2+)-binding site, has a high Phi-value, 0.91, indicating that tight packing of the side-chain around Asp87 occurs in the transition state. One beta-domain mutant (I55V) and three C-helix mutants (I89V, V90A, and I95V) demonstrated intermediate Phi-values, between 0.4 and 0.7. These results indicate that the folding nucleus in the transition state of goat alpha-LA is not extensively distributed over the alpha-domain of the protein, but very localized in a region that contains the Ca(2+)-binding site and the interface between the C-helix and the beta-domain. This is apparently in contrast with the fact that the molten globule state of alpha-lactalbumin has a partially formed structure inside the alpha-domain. It is concluded that the specific docking of the alpha and beta-domains at a domain interface is necessary for this protein to organize its native structure from the molten globule intermediate.     

Unfolding transition state and intermediates of the tumor suppressor p16INK4a investigated by molecular dynamics simulations

The ankyrin repeat is one of the most common protein motifs and is involved in protein-protein interactions. It consists of 33 residues that assume a beta-hairpin helix-loop-helix fold. Mutagenesis and kinetic experiments (Phi-value analysis of the folding transition state) have shown that the tumor suppressor p16(INK4a), a four-repeat protein, unfolds sequentially starting from the two N-terminal repeats. Here, the flexibility of p16(INK4a) at room temperature and its unfolding mechanism at high temperature have been investigated by multiple molecular dynamics runs in explicit water for a total simulation time of 0.65 micros. The transition state ensemble (TSE) of p16(INK4a) was identified by monitoring both the deviation from the experimental Phi values and sudden conformational changes along the unfolding trajectories. Conformations in the TSE have a mainly unstructured second repeat whereas the other repeats are almost completely folded. A rigid-body displacement of the first repeat involving both a rotation and translation is observed in all molecular dynamics simulations at high temperature. The Trp(15), Pro(75), and Ala(76) side-chains are more buried in the TSE than the native state. The sequential unfolding starting at the second repeat is in agreement with the mutagenesis studies whereas the displacement of the first repeat and the presence of nonnative interactions at the TSE are simulation results which supplement the experimental data. Furthermore, the unfolding trajectories reveal the presence of two on-pathway intermediates with partial alpha-helical structure. Finally, on the basis of the available experimental and simulation results we suggest that in modular proteins the shift of the folding TSE toward the native structure upon reduction of the number of tandem repeats is consistent with the Hammond effect.     

On the roles of substrate binding and hinge unfolding in conformational changes of adenylate kinase

We characterized the conformational change of adenylate kinase (AK) between open and closed forms by conducting five all-atom molecular-dynamics simulations, each of 100 ns duration. Different initial structures and substrate binding configurations were used to probe the pathways of AK conformational change in explicit solvent, and no bias potential was applied. A complete closed-to-open and a partial open-to-closed transition were observed, demonstrating the direct impact of substrate-mediated interactions on shifting protein conformation. The sampled configurations suggest two possible pathways for connecting the open and closed structures of AK, affirming the prediction made based on available x-ray structures and earlier works of coarse-grained modeling. The trajectories of the all-atom molecular-dynamics simulations revealed the complexity of protein dynamics and the coupling between different domains during conformational change. Calculations of solvent density and density fluctuations surrounding AK did not show prominent variation during the transition between closed and open forms. Finally, we characterized the effects of local unfolding of an important hinge near Pro¹⁷⁷ on the closed-to-open transition of AK and identified a novel mechanism by which hinge unfolding modulates protein conformational change. The local unfolding of Pro¹⁷⁷ hinge induces alternative tertiary contacts that stabilize the closed structure and prevent the opening transition.     

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.     

Local and global cooperativity in the human alpha-lactalbumin molten globule

NMR spectroscopy has been used to follow the urea-induced unfolding of the low pH molten globule states of a single-disulfide variant of human alpha-lactalbumin ([28-111] alpha-LA) and of two mutants, each with a single proline substitution in a helix. [28-111] alpha-LA forms a molten globule very similar to that formed by the wild-type four-disulfide protein, and this variant has been used as a model for the alpha-lactalbumin (alpha-LA) molten globule in a number of studies. The urea-induced unfolding behavior of [28-111] alpha-LA is similar to that of the four-disulfide form of the protein, except that [28-111] alpha-LA is less stable and has greater cooperativity in the loss of different elements of structure. For one mutant, L11P, the helix containing the mutation is highly destabilized such that it is completely unfolded even in the absence of urea. By contrast, for the other mutant, Q117P, the helix containing the mutation retains its compact structure. Both mutations, however, show significant long-range destabilization of the overall fold showing that the molten globule state has a degree of global cooperativity. The results reveal that different permutations of three of the four major alpha-helices of the protein can form a stable, locally cooperative, compact structural core. Taken together, these findings demonstrate that the molten globule state of alpha-LA is an ensemble of conformations, with different subsets of structures linked by a range of long-range interactions.     

The structure of human apolipoprotein E2, E3 and E4 in solution. 2. Multidomain organization correlates with the stability of apoE structure

The stabilities toward thermal and chemical denaturation of three recombinant isoforms of human apolipoprotein E (r-apoE2, r-apoE3 and r-apoE4), human plasma apoE3, the recombinant amino-terminal (NT) and the carboxyl-terminal (CT) domains of plasma apoE3 at pH 7 were studied using near and far ultraviolet circular dichroism (UV CD), fluorescence and size-exclusion chromatography. By far UV CD, thermal unfolding was irreversible for the intact apoE isoforms and consisted of a single transition. The r-apoE3 was found to be less stable as compared to the plasma protein and the stability of recombinant isoforms was r-apoE4相似文献   

Role of solvation barriers in protein kinetic stability

The stability of several protein systems of interest has been shown to have a kinetic basis. Besides the obvious biotechnological implications, the general interest of understanding protein kinetic stability is emphasized by the fact that some emerging molecular approaches to the inhibition of amyloidogenesis focus on the increase of the kinetic stability of protein native states. Lipases are among the most important industrial enzymes. Here, we have studied the thermal denaturation of the wild-type form, four single-mutant variants and two highly stable, multiple-mutant variants of lipase from Thermomyces lanuginosa. In all cases, thermal denaturation was irreversible, kinetically controlled and conformed to the two-state irreversible model. This result supports that the novel molecular-dynamics-focused, directed-evolution approach involved in the preparation of the highly stable variants is successful likely because it addresses kinetic stability and, in particular, because heated molecular dynamics simulations possibly identify regions of disrupted native interactions in the transition state for irreversible denaturation. Furthermore, we find very large mutation effects on activation enthalpy and entropy, which were not accompanied by similarly large changes in kinetic urea m-value. From this we are led to conclude that these mutation effects are associated to some structural feature of the transition state for the irreversible denaturation process that is not linked to large changes in solvent accessibility. Recent computational studies have suggested the existence of solvation/desolvation barriers in at least some protein folding/unfolding processes. We thus propose that a solvation barrier (arising from the asynchrony between breaking of internal contacts and water penetration) may contribute to the kinetic stability of lipase from T. lanuginosa (and, possibly, to the kinetic stability of other proteins as well).     

Pressure denaturation of staphylococcal nuclease studied by neutron small-angle scattering and molecular simulation

We studied the pressure-induced folding/unfolding transition of staphylococcal nuclease (SN) over a pressure range of approximately 1-3 kilobars at 25 degrees C by small-angle neutron scattering and molecular dynamics simulations. We find that applying pressure leads to a twofold increase in the radius of gyration derived from the small-angle neutron scattering spectra, and P(r), the pair distance distribution function, broadens and shows a transition from a unimodal to a bimodal distribution as the protein unfolds. The results indicate that the globular structure of SN is retained across the folding/unfolding transition although this structure is less compact and elongated relative to the native structure. Pressure-induced unfolding is initiated in the molecular dynamics simulations by inserting water molecules into the protein interior and applying pressure. The P(r) calculated from these simulations likewise broadens and shows a similar unimodal-to-bimodal transition with increasing pressure. The simulations also reveal that the bimodal P(r) for the pressure-unfolded state arises as the protein expands and forms two subdomains that effectively diffuse apart during initial stages of unfolding. Hydrophobic contact maps derived from the simulations show that water insertions into the protein interior and the application of pressure together destabilize hydrophobic contacts between these two subdomains. The findings support a mechanism for the pressure-induced unfolding of SN in which water penetration into the hydrophobic core plays a central role.     

Unfolding kinetics of beta-lactoglobulin induced by surfactant and denaturant: a stopped-flow/fluorescence study

The beta-->alpha transition of beta-lactoglobulin, a globular protein abundant in the milk of several mammals, is investigated in this work. This transition, induced by the cationic surfactant dodecyltrimethylammonium chloride (DTAC), is accompanied by partial unfolding of the protein. In this work, unfolding of bovine beta-lactoglobulin in DTAC is compared with its unfolding induced by the chemical denaturant guanidine hydrochloride (GnHCl). The final protein states attained in the two media have quite different secondary structure: in DTAC the alpha-helical content increases, leading to the so-called alpha-state; in GnHCl the amount of ordered secondary-structure decreases, resulting in a random coil-rich final state (denatured, or D, state). To obtain information on both mechanistic routes, in DTAC and GnHCl, and to characterize intermediates, the kinetics of unfolding were investigated in the two media. Equilibrium and kinetic data show the partial accumulation of an on-pathway intermediate in each unfolding route: in DTAC, an intermediate (I(1)) with mostly native secondary structure but loose tertiary structure appears between the native (beta) and alpha-states; in GnHCl, another intermediate (I(2)) appears between states beta and D. Kinetic rate constants follow a linear Chevron-plot representation in GnHCl, but show a more complex mechanism in DTAC, which acts like a stronger binding species.     

A stopped-flow fluorescence study of the native and modified lysozyme

The protein folding kinetics of hen egg white lysozyme (HEWL) was studied using experimental and bioinformatics tools. The structure of the transition state in the unfolding pathway of lysozyme was determined with stopped-flow kinetics using intact HEWL and its chemically modified derivative, in which six lysine residues have been modified. The overall consistency of φ-value (φ ≈ 1) indicates that lysine side chains interactions are subject to breaking in the structure of the transition state. Following experimental evidences, multiple sequence alignment of lysozyme family in vertebrates and exact structural examination of lysozyme, showed that the α-helix in the structure of lysozyme has critical role in the unfolding kinetics.     

Structural and mechanistic exploration of acid resistance: kinetic stability facilitates evolution of extremophilic behavior

Kinetically stable proteins are unique in that their stability is determined solely by kinetic barriers rather than by thermodynamic equilibria. To better understand how kinetic stability promotes protein survival under extreme environmental conditions, we analyzed the unfolding behavior and determined the structure of Nocardiopsis alba Protease A (NAPase), an acid-resistant, kinetically stable protease, and compared these results with a neutrophilic homolog, α-lytic protease (αLP). Although NAPase and αLP have the same number of acid-titratable residues, kinetic studies revealed that the height of the unfolding free energy barrier for NAPase is less sensitive to acid than that of αLP, thereby accounting for NAPase's improved tolerance of low pH. A comparison of the αLP and NAPase structures identified multiple salt-bridges in the domain interface of αLP that were relocated to outer regions of NAPase, suggesting a novel mechanism of acid stability in which acid-sensitive electrostatic interactions are rearranged to similarly affect the energetics of both the native state and the unfolding transition state. An acid-stable variant of αLP in which a single interdomain salt-bridge is replaced with a corresponding intradomain NAPase salt-bridge shows a dramatic > 15-fold increase in acid resistance, providing further evidence for this hypothesis. These observations also led to a general model of the unfolding transition state structure for αLP protease family members in which the two domains separate from each other while remaining relatively intact themselves. These results illustrate the remarkable utility of kinetic stability as an evolutionary tool for developing longevity over a broad range of harsh conditions.     

Characterization of the unfolding pathway of the cell-cycle protein p13suc1 by molecular dynamics simulations: implications for domain swapping

BACKGROUND: The p13suc1 gene product is a member of the cks (cyclin-dependent protein kinase subunit) protein family and has been implicated in regulation of the cell cycle. Various crystal structures of suc1 are available, including a globular, monomeric form and a beta-strand exchanged dimer. It has been suggested that conversions between these forms, and perhaps others, may be important in the regulation of the cell cycle. RESULTS: We have undertaken molecular dynamics simulations of protein unfolding to investigate the conformational properties of suc1. Unfolding transition states were identified for each of four simulations. These states contain some native secondary structure, primarily helix alpha1 and the core of the beta sheet. The hydrophobic core is loosely packed. Further unfolding leads to an intermediate state that is slightly more expanded than the transition state, but with considerably fewer nonlocal, tertiary packing contacts and less secondary structure. The helices are fluctuating but partially formed in the denatured state and beta2 and beta4 remain associated. CONCLUSIONS: It appears that suc1 folds by a nucleation-condensation mechanism, similar to that observed for two-state folding proteins. However, suc1 forms an intermediate during unfolding and contains considerable residual structure in the denatured state. The stability of the beta2-beta4 residual structure is surprising, because beta4 is the strand involved in domain swapping. This stability suggests that the domain-swapping event, if physiologically relevant, may require the assistance of additional factors in vivo or occur early in the folding process.     

Protein unfolding transitions in an intrinsically unstable annexin domain: molecular dynamics simulation and comparison with nuclear magnetic resonance data

Unfolding transitions of an intrinsically unstable annexin domain and the unfolded state structure have been examined using multiple approximately 10-ns molecular dynamics simulations. Three main basins are observed in the configurational space: native-like state, compact partially unfolded or intermediate compact state, and the unfolded state. In the native-like state fluctuations are observed that are nonproductive for unfolding. During these fluctuations, after an initial loss of approximately 20% of the core residue native contacts, the core of the protein transiently completely refolds to the native state. The transition from the native-like basin to the partially unfolded compact state involves approximately 75% loss of native contacts but little change in the radius of gyration or core hydration properties. The intermediate state adopts for part of the time in one of the trajectories a novel highly compact salt-bridge stabilized structure that can be identified as a conformational trap. The intermediate-to-unfolded state transition is characterized by a large increase in the radius of gyration. After an initial relaxation the unfolded state recovers a native-like topology of the domain. The simulated unfolded state ensemble reproduces in detail experimental nuclear magnetic resonance data and leads to a convincing complete picture of the unfolded domain.     

Mesophile versus thermophile: insights into the structural mechanisms of kinetic stability

Obtaining detailed knowledge of folding intermediate and transition state (TS) structures is critical for understanding protein folding mechanisms. Comparisons between proteins adapted to survive extreme temperatures with their mesophilic homologs are likely to provide valuable information on the interactions relevant to the unfolding transition. For kinetically stable proteins such as alpha-lytic protease (alphaLP) and its family members, their large free energy barrier to unfolding is central to their biological function. To gain new insights into the mechanisms that underlie kinetic stability, we have determined the structure and high temperature unfolding kinetics of a thermophilic homolog, Thermobifida fusca protease A (TFPA). These studies led to the identification of a specific structural element bridging the N and C-terminal domains of the protease (the "domain bridge") proposed to be associated with the enhanced high temperature kinetic stability in TFPA. Mutagenesis experiments exchanging the TFPA domain bridge into alphaLP validate this hypothesis and illustrate key structural details that contribute to TFPA's increased kinetic thermostability. These results lead to an updated model for the unfolding transition state structure for this important class of proteases in which domain bridge undocking and unfolding occurs at or before the TS. The domain bridge appears to be a structural element that can modulate the degree of kinetic stability of the different members of this class of proteases.     

Cooperative unfolding of a metastable serpin to a molten globule suggests a link between functional and folding energy landscapes

Alpha-1 antitrypsin (alpha(1)-AT) is a member of the serpin class of protease inhibitors, and folds to a metastable state rather than its thermodynamically most stable native state. Upon cleavage by a target protease, alpha(1)-AT undergoes a dramatic conformational change to a stable form, translocating the bound protease more than 70 A to form an inhibitory protease-serpin complex. Numerous mutagenesis studies on serpins have demonstrated the trade-off between the stability of the metastable state on the one hand and the inhibitory efficiency on the other. Studies of the equilibrium unfolding of serpins provide insight into this connection between structural plasticity and metastability. We studied equilibrium unfolding of wild-type alpha(1)-AT using hydrogen-deuterium/exchange mass spectrometry to characterize the structure and the stability of an equilibrium intermediate that was observed in low concentrations of denaturant in earlier studies. Our results show that the intermediate observed at low concentrations of denaturant has no protection from hydrogen-deuterium exchange, indicating a lack of stable structure. Further, differential scanning calorimetry of alpha(1)-AT at low concentrations of denaturant shows no heat capacity peak during thermal denaturation, indicating that the transition from the intermediate to the unfolded state is not a cooperative first-order-like phase transition.. Our results show that the unfolding of alpha(1)-AT involves a cooperative transition to a molten globule form, followed by a non-cooperative transition to a random-coil form as more guanidine is added. Thus, the entire alpha(1)-AT molecule consists of one cooperative structural unit rather than multiple structural domains with different stabilities. Furthermore, our results together with previous mutagenesis studies suggest a possible link between an equilibrium molten globule and a functional intermediate that may be populated during the protease inhibition.