Comparison of solution structures and stabilities of native, partially unfolded and partially refolded pepsin

A zymogen-derived protein, pepsin, appears to be incapable of folding to the native state without the presence of the prosegment. To better understand the nature of the irreversible denaturation of pepsin, the present study reports on the characterization of the stability and low-resolution tertiary and secondary structures of native, alkaline unfolded and acid refolded porcine pepsin. Through a combination of small-angle neutron scattering (SANS), CD, and DSC, acid refolded pepsin (Rp) was shown to have secondary and tertiary structures intermediate between the alkaline denatured and native forms but was found to be thermodynamically stable relative to the native state. It was also observed that the acid refolded state of pepsin was dependent on the protein concentration during refolding because CD and SANS data revealed that both the secondary and tertiary structures of concentrated-refolded pepsin (>10 mg/mL) (CRp) were native-like, in contrast to the intermediate nature of Rp, refolded under dilute concentration (<10 mg/mL). Despite a native-like conformation, CRp was more stable and had substantially reduced activity compared to that of the native state, suggesting that the protein was misfolded. It is proposed that the stable but misfolded, acid-refolded states are evidence that pepsin in its native conformation was metastable. Furthermore, the disruption of the active site cleft in the denatured states could be discerned by modeling of the SANS data.     

Equilibrium unfolding of DLC8 monomer by urea and guanidine hydrochloride: Distinctive global and residue level features

We present circular dichroism (CD), steady state fluorescence and multidimensional NMR investigations on the equilibrium unfolding of monomeric dynein light chain protein (DLC8) by urea and guanidine hydrochloride (GdnHCl). Quantitative analysis of the CD and fluorescence denaturation curves reveals that urea unfolding is a two-state process, whereas guanidine unfolding is more complex. NMR investigations in the native state and in the near native states created by low denaturant concentrations enabled residue level characterization of the early structural and dynamic perturbations by the two denaturants. Firstly, (15)N transverse relaxation rates in the native state indicate that the regions around N10, Q27, the loop between beta2 and beta4 strands, and K87 at the C-terminal are potential unfolding initiation sites in the protein. Amide and (15)N chemical shift perturbations indicate different accessibilities of the residues along the chain and help identify locations of the early perturbations by the two denaturants. Guanidine and urea are seen to interact at several sites some of which are different in the two cases. Notable among the common interaction site is that around K87 which is in close proximity to W54 on the protein structure, but the interaction modes of the two denaturants are different. The secondary chemical shifts indicate that the structural perturbation by 1M urea is small, compared to that by guanidine which is more encompassing over the length of the chain. The probable (phi, psi) changes at the individual residues have been calculated using the TALOS algorithm. It appears that the helices in the protein are significantly perturbed by guanidine. Further, comparison of the spectral density functions of the native and the two near native states in the two denaturants implicate greater loosening of the structure by guanidine as compared to that by urea, even though the structures are still in the native state ensemble. These differences in the early perturbations of the native state structure and dynamics by the two denaturants might direct the protein along different pathways, as the unfolding progresses on further increasing the denaturant concentration.     

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.     

Exploring Early Stages of the Chemical Unfolding of Proteins at the Proteome Scale

After decades of using urea as denaturant, the kinetic role of this molecule in the unfolding process is still undefined: does urea actively induce protein unfolding or passively stabilize the unfolded state? By analyzing a set of 30 proteins (representative of all native folds) through extensive molecular dynamics simulations in denaturant (using a range of force-fields), we derived robust rules for urea unfolding that are valid at the proteome level. Irrespective of the protein fold, presence or absence of disulphide bridges, and secondary structure composition, urea concentrates in the first solvation shell of quasi-native proteins, but with a density lower than that of the fully unfolded state. The presence of urea does not alter the spontaneous vibration pattern of proteins. In fact, it reduces the magnitude of such vibrations, leading to a counterintuitive slow down of the atomic-motions that opposes unfolding. Urea stickiness and slow diffusion is, however, crucial for unfolding. Long residence urea molecules placed around the hydrophobic core are crucial to stabilize partially open structures generated by thermal fluctuations. Our simulations indicate that although urea does not favor the formation of partially open microstates, it is not a mere spectator of unfolding that simply displaces to the right of the folded←→unfolded equilibrium. On the contrary, urea actively favors unfolding: it selects and stabilizes partially unfolded microstates, slowly driving the protein conformational ensemble far from the native one and also from the conformations sampled during thermal unfolding.     

Is protein unfolding the reverse of protein folding? A lattice simulation analysis.

Simulations and experiments that monitor protein unfolding under denaturing conditions are commonly employed to study the mechanism by which a protein folds to its native state in a physiological environment. Due to the differences in conditions and the complexity of the reaction, unfolding is not necessarily the reverse of folding. To assess the relevance of temperature initiated unfolding studies to the folding problem, we compare the folding and unfolding of a 125-residue protein model by Monte Carlo dynamics at two temperatures; the lower one corresponds to the range used in T -jump experiments and the higher one to the range used in unfolding simulations of all-atom models. The trajectories that lead from the native state to the denatured state at these elevated temperatures are less diverse than those observed in the folding simulations. At the lower temperature, the system unfolds through a mandatory intermediate that corresponds to a local free energy minimum. At the higher temperature, no such intermediate is observed, but a similar pathway is followed. The structures contributing to the unfolding pathways resemble most closely those that make up the "fast track" of folding. The transition state for unfolding at the lower temperature (above Tm) is determined and is found to be more structured than the transition state for folding below the melting temperature. This shift towards the native state is consistent with the Hammond postulate. The implications for unfolding simulations of higher resolution models and for unfolding experiments of proteins are discussed.     

Characterization of the denatured states distribution of neocarzinostatin by small-angle neutron scattering and differential scanning calorimetry

The denatured states of a small globular protein, apo-neocarzinostatin (NCS), have been characterized using several techniques. Structural properties were investigated by optical spectroscopy techniques and small-angle neutron scattering (SANS), as a function of guanidinium chloride (GdmCl) concentration. SANS experiments show that in heavy water, the protein keeps its native size at GdmCl concentrations below 2.5 M. A sharp transition occurs at about 3.6 M GdmCl, and NCS behaves like an excluded volume chain above 5 M. The same behavior is observed in deuterated buffer by fluorescence and circular dichroism measurements. For the H(2)O buffer, the transition occurs with lower concentration of denaturant, the shift being about 0.6 M. 8-Anilino-1-naphthalenesulfonate (ANS) was used as a hydrophobic fluorescent probe for studying the early stages of protein unfolding. Protein denaturation modifies the fluorescence intensity of ANS, a maximum of intensity being detected close to 2 M GdmCl in hydrogenated buffer, which shows the existence of at least one intermediate state populated at the beginning of the unfolding pathway. Differential scanning calorimetry (DSC) was used to obtain thermodynamic values for NCS denaturation. The melting curves recorded between 20 and 90 degrees C in the presence of various GdmCl concentrations (0-3 M) cannot be explained by a simple two-state model. Altogether, the data presented in this paper suggest that before unfolding the protein explores a distribution of states which is centered around compact states at denaturant concentrations below 2 M in H(2)O, and then shifts to less structured states by increasing denaturant concentrations.     

Comparative Fourier transform infrared spectroscopy study of cold-, pressure-, and heat-induced unfolding and aggregation of myoglobin

We studied the cold unfolding of myoglobin with Fourier transform infrared spectroscopy and compared it with pressure and heat unfolding. Because protein aggregation is a phenomenon with medical as well as biotechnological implications, we were interested in both the structural changes as well as the aggregation behavior of the respective unfolded states. The cold- and pressure-induced unfolding both yield a partially unfolded state characterized by a persistent amount of secondary structure, in which a stable core of G and H helices is preserved. In this respect the cold- and pressure-unfolded states show a resemblance with an early folding intermediate of myoglobin. In contrast, the heat unfolding results in the formation of the infrared bands typical of intermolecular antiparallel beta-sheet aggregation. This implies a transformation of alpha-helix into intermolecular beta-sheet. H/2H-exchange data suggest that the helices are first unfolded and then form intermolecular beta-sheets. The pressure and cold unfolded states do not give rise to the intermolecular aggregation bands that are typical for the infrared spectra of many heat-unfolded proteins. This suggests that the pathways of the cold and pressure unfolding are substantially different from that of the heat unfolding. After return to ambient conditions the cold- or pressure-treated proteins adopt a partially refolded conformation. This aggregates at a lower temperature (32 degrees C) than the native state (74 degrees C).     

High-pressure SANS and fluorescence unfolding study of calmodulin

Apo-calmodulin, a small soluble mainly α protein, is a calcium-dependent protein activator. Calcium binding affects the calmodulin conformation but also its stability. Calcium free form unfolds between 40 and 80 °C, whereas the calcium-saturated form is stable up to temperatures as high as 100 °C, forbidding comparison of the thermal unfolding pathways of the two forms. Thus, this paper focuses especially on the conformation of pressure-induced unfolding states of both forms of calmodulin, by combining small-angle neutron scattering (SANS) with biophysical techniques such as tyrosines and ANS fluorescence. In contrast to heat denaturation (Gibrat et al., BBA, 2012), the pressure denaturation of calmodulin is reversible up to pressures of 3000 bar (300 MPa). A pressure-induced compact intermediate state has been found for the two calmodulin forms, but their unfolding pathways are different. A domain compaction and an increase of the ANS fluorescence of holo form have been evidenced. On the contrary, a domain dilatation and an ANS fluorescence decrease have been found for the apo form. The pressure induced an increase of the interdomain distance for both calmodulin forms, suggesting that the central linker of calmodulin is flexible in solution.     

Conformational stability of human erythrocyte transglutaminase. Patterns of thermal unfolding at acid and alkaline pH.

Tissue-type transglutaminase is irreversibly inactivated during heat treatment. The rate of inactivation is low at pH 7.5; it increases slightly at acid pH (6.1) but much more at alkaline pH (9.0-9.5), suggesting that specific effects take place in the alkaline range, possibly in relation to decreased stability of the transition-state intermediate as pH is raised above 9.0. Differential scanning calorimetry experiments indicate that thermal unfolding of the protein occurs with two separate transitions, involving independent regions of the enzyme. They are assigned to domains 1 and 2 and domains 3 and 4, respectively, by a combination of calorimetric and spectroscopic techniques. When considering the effects of pH, we noted that transglutaminase was unfolded via different pathways at the different pH values considered. At acid pH, the whole structure of the protein was lost irreversibly, with massive aggregation. At neutral and, even more so, at alkaline pH, aggregation was absent (or very limited at high protein concentration) and the loss of secondary structure was dependent on the ionization state of crucial lysine residues. Unfolding at pH 9.5 apparently chiefly involved the N-terminal region, as testified by changes in protein intrinsic fluorescence. In addition, the C-terminal region was destabilized at each pH value tested during thermal unfolding, as shown by digestion with V8 proteinase, which is inactive on the native protein. Evidence was obtained that the N-terminal and C-terminal regions interact with each other in determining the structure of the native protein.     

The heat capacity of proteins

The heat capacity plays a major role in the determination of the energetics of protein folding and molecular recognition. As such, a better understanding of this thermodynamic parameter and its structural origin will provide new insights for the development of better molecular design strategies. In this paper we have analyzed the absolute heat capacity of proteins in different conformations. The results of these studies indicate that three major terms account for the absolute heat capacity of a protein: (1) one term that depends only on the primary or covalent structure of a protein and contains contributions from vibrational frequencies arising from the stretching and bending modes of each valence bond and internal rotations; (2) a term that contains the contributions of noncovalent interactions arising from secondary and tertiary structure; and (3) a term that contains the contributions of hydration. For a typical globular protein in solution the bulk of the heat capacity at 25°C is given by the covalent structure term (close to 85% of the total). The hydration term contributes about 15 and 40% to the total heat capacity of the native and unfolded states, respectively. The contribution of non-covalent structure to the total heat capacity of the native state is positive but very small and does not amount to more than 3% at 25°C. The change in heat capacity upon unfolding is primarily given by the increase in the hydration term (about 95%) and to a much lesser extent by the loss of noncovalent interactions (up to ～5%). It is demonstrated that a single universal mathematical function can be used to represent the partial molar heat capacity of the native and unfolded states of proteins in solution. This function can be experimentally written in terms of the molecular weight, the polar and apolar solvent accessible surface areas, and the total area buried from the solvent. This unique function accurately predicts the different magnitude and temperature dependences of the heat capacity of both the native and unfolded states, and therefore of the heat capacity changes associated with folding/unfolding transitions. © 1995 Wiley-Liss, Inc.     

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.     

Structure of an Unfolding Intermediate of an RRM Domain of ETR-3 Reveals Its Native-like Fold

Prevalence of one or more partially folded intermediates during protein unfolding with different secondary and ternary conformations has been identified as an integral character of protein unfolding. These transition-state species need to be characterized structurally for elucidation of their folding pathways. We have determined the three-dimensional structure of an intermediate state with increased conformational space sampling under urea-denaturing condition. The protein unfolds completely at 10 M urea but retains residual secondary structural propensities with restricted motion. Here, we describe the native state, observable intermediate state, and unfolded state for ETR-3 RRM-3, which has canonical RRM fold. These observations can shed more light on unfolding events for RRM-containing proteins.     

Alternatively folded states of an immunoglobulin

Well-defined, non-native protein structures of low stability have been increasingly observed as intermediates in protein folding or as equilibrium structures populated under specific solvent conditions. These intermediate structures, frequently referred to as molten globule states, are characterized by the presence of secondary structure, a lack of significant tertiary contacts, increased hydrophobicity and partial specific volume as compared to native structures, and low cooperativity in thermal unfolding. The present study demonstrates that under acidic conditions (pH less than 3) the antibody MAK33 can assume a folded stable conformation. This A-state is characterized by a high degree of secondary structure, increased hydrophobicity, a native-like maximum wavelength of fluorescence emission, and a tendency toward slow aggregation. A prominent feature of this low-pH conformation is the stability against denaturant and thermal unfolding that is manifested in highly cooperative reversible phase transitions indicative of the existence of well-defined tertiary contacts. These thermodynamic results are corroborated by the kinetics of folding from the completely unfolded chain to the alternatively folded state at pH 2. The given data suggest that MAK33 at pH 2 adopts a cooperative structure that differs from the native immunoglobulin fold at pH 7. This alternatively folded state exhibits certain characteristics of the molten globule but differs distinctly from it by its extraordinary structural stability that is characteristic for native protein structures.     

Biophysical study of thermal denaturation of apo-calmodulin: dynamics of native and unfolded states

Apo-calmodulin, a small, mainly α, soluble protein is a calcium-dependent protein activator. This article presents a study of internal dynamics of native and thermal unfolded apo-calmodulin, using quasi-elastic neutron scattering. This technique can probe protein internal dynamics in the picosecond timescale and in the nanometer length-scale. It appears that a dynamical transition is associated with thermal denaturation of apo-calmodulin. This dynamical transition goes together with a decrease of the confinement of hydrogen atoms, a decrease of immobile protons proportion and an increase of dynamical heterogeneity. The comparison of native and unfolded states dynamics suggests that the dynamics of protein atoms is more influenced by their distance to the backbone than by their solvent exposure.     

Titration Properties and Thermodynamics of the Transition State for Folding: Comparison of Two-state and Multi-state Folding Pathways

CI2 folds and unfolds as a single cooperative unit by simple two-state kinetics, which enables the properties of the transition state to be measured from both the forward and backward rate constants. We have examined how the free energy of the transition state for the folding of chymotrypsin inhibitor 2 (CI2) changes with pH and temperature. In addition to the standard thermodynamic quantities, we have measured the overall acid-titration properties of the transition state and its heat capacity relative to both the denatured and native states. We were able to determine the latter by a method analogous to a well-established procedure for measuring the change in heat capacity for equilibrium unfolding: the enthalpy of activation of unfolding at different values of acid pH were plotted against the average temperature of each determination. Our results show that the transition state of CI2 has lost most of the electrostatic and van der Waals' interactions that are found in the native state, but it remains compact and this prevents water molecules from entering some parts of the hydrophobic core. The properties of the transition state of CI2 are then compared with the major folding transition state of the larger protein barnase, which folds by a multi-state mechanism, with the accumulation of a partly structured intermediate (D^physor I). CI2 folds from a largely unstructured denatured state under physiological conditionsviaa transition state which is compact but relatively uniformly unstructured, with tertiary and secondary structure being formed in parallel. We term this an expanded pathway. Conversely, barnase folds from a largely structured denatured state in which elements of structure are well formed through a transition state that has islands of folded elements of structure. We term this a compact pathway. These two pathways may correspond to the two extreme ends of a continuous spectrum of protein folding mechanisms. Although the properties of the two transition states are very different, the activation barrier for folding (D^phys→3 ) is very similar for both proteins.     

Thermal effects in stretching of Go-like models of titin and secondary structures

The effect of temperature on mechanical unfolding of proteins is studied using a Go-like model with a realistic contact map and Lennard-Jones contact interactions. The behavior of the I27 domain of titin and its serial repeats is contrasted to that of simple secondary structures. In all cases, thermal fluctuations accelerate the unraveling process, decreasing the unfolding force nearly linearly at low temperatures. However, differences in bonding geometry lead to different sensitivity to temperature and different changes in the unfolding pattern. Due to its special native-state geometry, titin is much more thermally and elastically stable than the secondary structures. At low temperatures, serial repeats of titin show a parallel unfolding of all domains to an intermediate state, followed by serial unfolding of the domains. At high temperatures, all domains unfold simultaneously, and the unfolding distance decreases monotonically with the contact order, that is, the sequence distance between the amino acids that form the native contact.