Role of disulfide bonds in goose-type lysozyme

The role of the two disulfide bonds (Cys4-Cys60 and Cys18-Cys29) in the activity and stability of goose-type (G-type) lysozyme was investigated using ostrich egg-white lysozyme as a model. Each of the two disulfide bonds was deleted separately or simultaneously by substituting both Cys residues with either Ser or Ala. No remarkable differences in secondary structure or catalytic activity were observed between the wild-type and mutant proteins. However, thermal and guanidine hydrochloride unfolding experiments revealed that the stabilities of mutants lacking one or both of the disulfide bonds were significantly decreased relative to those of the wild-type. The destabilization energies of mutant proteins agreed well with those predicted from entropic effects in the denatured state. The effects of deleting each disulfide bond on protein stability were found to be approximately additive, indicating that the individual disulfide bonds contribute to the stability of G-type lysozyme in an independent manner. Under reducing conditions, the thermal stability of the wild-type was decreased to a level nearly equivalent to that of a Cys-free mutant (C4S/C18S/C29S/C60S) in which all Cys residues were replaced by Ser. Moreover, the optimum temperature of the catalytic activity for the Cys-free mutant was downshifted by about 20 degrees C as compared with that of the wild-type. These results indicate that the formation of the two disulfide bonds is not essential for the correct folding into the catalytically active conformation, but is crucial for the structural stability of G-type lysozyme.     

Role of disulfide bonds in the stability of recombinant manganese peroxidase.

Phanerochaete chrysosporium manganese peroxidase (MnP) [isoenzyme H4] was engineered with additional disulfide bonds to provide structural reinforcement to the proximal and distal calcium-binding sites. This rational protein engineering investigated the effects of multiple disulfide bonds on the stabilization of the enzyme heme environment and oxidase activity. Stabilization of the heme environment was monitored by UV-visible spectroscopy based on the electronic state of the alkaline transition species of ferric and ferrous enzyme. The optical spectral data confirm an alkaline transition to hexacoordinate, low-spin heme species for native and wild-type MnP and show that the location of the engineered disulfide bonds in the protein can have significant effects on the electronic state of the enzyme. The addition of a single disulfide bond in the distal region of MnP resulted in an enzyme that maintained a pentacoordinate, high-spin heme at pH 9.0, whereas MnP with multiple engineered disulfide bonds did not exhibit an increase in stability of the pentacoordinate, high-spin state of the enzyme at alkaline pH. The mutant enzymes were assessed for increased stability by incubation at high pH. In comparison to wild-type MnP, enzymes containing engineered disulfide bonds in the distal and proximal regions of the protein retained greater levels of activity when restored to physiological pH. Additionally, when assayed for oxidase activity at pH 9.0, proteins containing engineered disulfide bonds exhibited slower rates of inactivation than wild-type MnP.     

Disulfide bonds and the stability of globular proteins.

An understanding of the forces that contribute to stability is pivotal in solving the protein-folding problem. Classical theory suggests that disulfide bonds stabilize proteins by reducing the entropy of the denatured state. More recent theories have attempted to expand this idea, suggesting that in addition to configurational entropic effects, enthalpic and native-state effects occur and cannot be neglected. Experimental thermodynamic evidence is examined from two sources: (1) the disruption of naturally occurring disulfides, and (2) the insertion of novel disulfides. The data confirm that enthalpic and native-state effects are often significant. The experimental changes in free energy are compared to those predicted by different theories. The differences between theory and experiment are large near 300 K and do not lend support to any of the current theories regarding the stabilization of proteins by disulfide bonds. This observation is a result of not only deficiencies in the theoretical models but also from difficulties in determining the effects of disulfide bonds on protein stability against the backdrop of numerous subtle stabilizing factors (in both the native and denatured states), which they may also affect.     

Stabilization of an immunoglobulin fold domain by an engineered disulfide bond at the buried hydrophobic region

We report for the first time the stabilization of an immunoglobulin fold domain by an engineered disulfide bond. In the llama single-domain antibody, which has human chorionic gonadotropin as its specific antigen, Ala49 and Ile70 are buried in the structure. A mutant with an artificial disulfide bond at this position showed a 10 degrees C higher midpoint temperature of thermal unfolding than that without the extra disulfide bond. The modified domains exhibited an antigen binding affinity comparable with that of the wild-type domain. Ala49 and Ile70 are conserved in camel and llama single-domain antibody frameworks. Therefore, domains against different antigens are expected to be stabilized by the engineered disulfide bond examined here. In addition to the effect of the loop constraints in the unfolded state, thermodynamic analysis indicated that internal interaction and hydration also control the stability of domains with disulfide bonds. The change in physical properties resulting from mutation often causes unpredictable and destabilizing effects on these interactions. The introduction of a hydrophobic cystine into the hydrophobic region maintains the hydrophobicity of the protein and is expected to minimize the unfavorable mutational effects.     

Disulfide bond effects on protein stability: designed variants of Cucurbita maxima trypsin inhibitor-V

Attempts to increase protein stability by insertion of novel disulfide bonds have not always been successful. According to the two current models, cross-links enhance stability mainly through denatured state effects. We have investigated the effects of removal and addition of disulfide cross-links, protein flexibility in the vicinity of a cross-link, and disulfide loop size on the stability of Cucurbita maxima trypsin inhibitor-V (CMTI-V; 7 kD) by differential scanning calorimetry. CMTI-V offers the advantage of a large, flexible, and solvent-exposed loop not involved in extensive intra-molecular interactions. We have uncovered a negative correlation between retention time in hydrophobic column chromatography, a measure of protein hydrophobicity, and melting temperature (T(m)), an indicator of native state stabilization, for CMTI-V and its variants. In conjunction with the complete set of thermodynamic parameters of denaturation, this has led to the following deductions: (1) In the less stable, disulfide-removed C3S/C48S (Delta Delta G(d)(50 degrees C) = -4 kcal/mole; Delta T(m) = -22 degrees C), the native state is destabilized more than the denatured state; this also applies to the less-stable CMTI-V* (Delta Delta G(d)(50 degrees C) = -3 kcal/mole; Delta T(m) = -11 degrees C), in which the disulfide-containing loop is opened by specific hydrolysis of the Lys(44)-Asp(45) peptide bond; (2) In the less stable, disulfide-inserted E38C/W54C (Delta Delta G(d)(50 degrees C) = -1 kcal/mole; Delta T(m) = +2 degrees C), the denatured state is more stabilized than the native state; and (3) In the more stable, disulfide-engineered V42C/R52C (Delta Delta G(d)(50 degrees C) = +1 kcal/mole; Delta T(m) = +17 degrees C), the native state is more stabilized than the denatured state. These results show that a cross-link stabilizes both native and denatured states, and differential stabilization of the two states causes either loss or gain in protein stability. Removal of hydrogen bonds in the same flexible region of CMTI-V resulted in less destabilization despite larger changes in the enthalpy and entropy of denaturation. The effect of a cross-link on the denatured state of CMTI-V was estimated directly by means of a four-state thermodynamic cycle consisting of native and denatured states of CMTI-V and CMTI-V*. Overall, the results show that an enthalpy-entropy compensation accompanies disulfide bond effects and protein stabilization is profoundly modulated by altered hydrophobicity of both native and denatured states, altered flexibility near the cross-link, and residual structure in the denatured state.     

Multi-template approach to modeling engineered disulfide bonds

The key issue for disulfide bond engineering is to select the most appropriate location in the protein. By surveying the structure of experimentally engineered disulfide bonds, we found about half of them that have geometry incompatible with any native disulfide bond geometry. To improve the current prediction methods that tend to apply either ideal geometrical or energetical criteria to single three-dimensional structures, we have combined a novel computational protocol with the usage of multiple protein structures to take into account protein backbone flexibility. The multiple structures can be selected from either independently determined crystal structures for identical proteins, models of nuclear magnetic resonance experiments, or crystal structures of homology-related proteins. We have validated our approach by comparing the predictions with known disulfide bonds. The accuracy of prediction for native disulfide bonds reaches 99.6%. In a more stringent test on the reported engineered disulfide bonds, we have obtained a success rate of 93%. Our protocol also determines the oxido-reduction state of a predicted disulfide bond and the corresponding mutational cost. From the energy ranking, the user can easily choose top predicted sites for mutagenesis experiments. Our method provides information about local stability of the engineered disulfide bond surroundings.     

Denaturation and unfolding of human anaphylatoxin C3a: an unusually low covalent stability of its native disulfide bonds

The complement C3a anaphylatoxin is a major molecular mediator of innate immunity. It is a potent activator of mast cells, basophils and eosinophils and causes smooth muscle contraction. Structurally, C3a is a relatively small protein (77 amino acids) comprising a N-terminal domain connected by 3 native disulfide bonds and a helical C-terminal segment. The structural stability of C3a has been investigated here using three different methods: Disulfide scrambling; Differential CD spectroscopy; and Reductive unfolding. Two uncommon features regarding the stability of C3a and the structure of denatured C3a have been observed in this study. (a) There is an unusual disconnection between the conformational stability of C3a and the covalent stability of its three native disulfide bonds that is not seen with other disulfide proteins. As measured by both methods of disulfide scrambling and differential CD spectroscopy, the native C3a exhibits a global conformational stability that is comparable to numerous proteins with similar size and disulfide content, all with mid-point denaturation of [GdmCl]_1/2 at 3.4-5 M. These proteins include hirudin, tick anticoagulant protein and leech carboxypeptidase inhibitor. However, the native disulfide bonds of C3a is 150-1000 fold less stable than those proteins as evaluated by the method of reductive unfolding. The 3 native disulfide bonds of C3a can be collectively and quantitatively reduced with as low as 1 mM of dithiothreitol within 5 min. The fragility of the native disulfide bonds of C3a has not yet been observed with other native disulfide proteins. (b) Using the method of disulfide scrambling, denatured C3a was shown to consist of diverse isomers adopting varied extent of unfolding. Among them, the most extensively unfolded isomer of denatured C3a is found to assume beads-form disulfide pattern, comprising Cys³⁶-Cys⁴⁹ and two disulfide bonds formed by two pair of consecutive cysteines, Cys²²-Cys²³ and Cys⁵⁶-Cys⁵⁷, a unique disulfide structure of polypeptide that has not been documented previously.     

Protein engineering of disulfide bonds in subtilisin BPN'

Five single-disulfide mutants were studied in subtilisin BPN', a cysteine-free, secreted serine protease from Bacillus amyloliquefaciens. The disulfides were engineered between residues 26-232, 29-119, 36-210, 41-80, and 148-243. These bonds connected a variety of secondary structural elements, located in buried or exposed positions at least 10 A from the catalytic Ser-221, and linked residues that were separated by 39 up to 206 amino acids. All disulfide bonds formed in the enzyme when the expressed protein was secreted from Bacillus subtilis, and the disulfides had only minor effects on the enzyme kinetics. Although these disulfide bonds varied by over 50-fold in their equilibrium constants for reduction with dithiothreitol, there was no correlation between the strength of the disulfide bond and the stability it imparted to the enzyme to irreversible inactivation. In some cases, the disulfide-bonded protein was stabilized greatly relative to its reduced counterpart. However, no disulfide mutant was substantially more stable than wild-type subtilisin BPN'. Some of these results can be rationalized by destabilizing effects of the cysteine mutations that disrupt interactions present in the folded enzyme structure. It is also possible that the rate of irreversible inactivation depends upon the kinetics and not the thermodynamics of unfolding and so the entropically stabilizing effect expected from a disulfide bond may not apply.     

Enhancing the thermal stability of avidin. Introduction of disulfide bridges between subunit interfaces

In this study we showed that tetrameric chicken avidin can be stabilized by introducing intermonomeric disulfide bridges between its subunits. These covalent bonds had no major effects on the biotin binding properties of the respective mutants. Moreover, one of the mutants (Avd-ccci) maintained its tetrameric integrity even in denaturing conditions. The new avidin forms Avd-ci and Avd-ccci, which have native --> denatured transition midpoints (T(m)) of 98.6 and 94.7 degrees C, respectively, in the absence of biotin, will find use in applications where extreme stability or minimal leakage of subunits is required. Furthermore, we showed that the intramonomeric disulfide bridges found in the wild-type avidin affect its stability. The mutant Avd-nc, in which this bridge was removed, had a lower T(m) in the absence of biotin than the wild-type avidin but showed comparable stability in the presence of biotin.     

Significance of the four carboxyl terminal amino acid residues of bovine pancreatic ribonuclease A for structural folding

The C-terminal amino acid residues of bovine pancreatic ribonuclease A (RNase A) form a core structure in the initial stage of the folding process that leads to the formation of the tertiary structure. In this paper, roles of the C-terminal four amino acids in the structure, function, and refolding were studied by use of recombinant mutant enzymes in which these residues were deleted or replaced. Purified mutant enzymes were analyzed for their secondary structure, thermal stability, and ability to regenerate from the denatured and reduced state. The C-terminal deleted mutant enzymes showed lower hydrolytic activity for C>p and nearly identical CD spectra compared with the wild-type enzyme. The rate of recovery of activity was significantly different among the C-terminal deleted mutant enzymes when air oxidation was employed in the absence of GSH and GSSG: the rates decreased in the order of des-124-, des-(123-124)-, and des-(122-124)-RNase A. It is noteworthy that the regeneration rates of mutant RNase A in the presence of GSH and GSSG were nearly the same. Des-(121-124)-RNase A failed to recover activity both in the presence and absence of glutathione, due to the mismatched formation of disulfide bonds. The mutant enzyme in which all of the C-terminal four amino acid residues were replaced by alanine residues showed lower hydrolytic activity and an indistinguishable CD spectrum compared with the wild-type enzyme, and also recovered its activity from the denatured and reduced state by air oxidation. The D121 mutant enzymes showed decreased hydrolytic activity and identical CD spectra compared with the wild type. The recovery rates of activity of D121A and D121K were determined to be lower than that of the wild-type enzyme, while the rate of recovery of D121E was comparable to that of the wild type. The C-terminal amino acids play a significant role in the formation of the correct disulfide bonds during the refolding process, and the interaction of amino acid residues and the existence of the main chain around the C-terminal region are both important for achieving the efficient packing of the RNase A molecule.     

Kinetic analysis of molecular dynamics simulations reveals changes in the denatured state and switch of folding pathways upon single-point mutation of a beta-sheet miniprotein

The effects of a single-point mutation on folding thermodynamics and kinetics are usually interpreted by focusing on the native structure and the transition state. Here, the entire conformational spaces of a 20-residue three-stranded antiparallel beta-sheet peptide (double hairpin) and of its single-point mutant W10V are sampled close to the melting temperature by equilibrium folding-unfolding molecular dynamics simulations for a total of 40 micros. The folded state as well as the most populated free energy basins in the denatured state are isolated by grouping conformations according to fast relaxation at equilibrium. Such kinetic analysis provides more detailed and useful information than a simple projection of the free energy. The W10V mutant has the same native structure as the wild type peptide, and similar folding rate and stability. In the denatured state, the N-terminal hairpin is about 20% more structured in W10V than the wild type mainly because of van der Waals interactions. Notably, the W10V mutation influences also the van der Waals energy at the transition state ensemble causing a shift in the ratio of fluxes between two different transition state regions on parallel folding pathways corresponding to nucleation at either of the two beta-hairpins. Previous experimental studies have focused on the effects of denaturant-dependent or temperature-dependent changes in the structure of the denatured state. The atomistic simulations show that a single-point mutation in the central strand of a beta-sheet peptide results in remarkable changes in the topography of the denatured state ensemble. These changes modulate the relative accessibility of parallel folding pathways because of kinetic partitioning of the denatured state. Therefore, the observed dependence of the folding process on the starting ensemble raises questions on the biological significance of in vitro folding studies under strongly denaturing conditions.     

Dynamics and thermodynamic properties of CXCL7 chemokine

Chemokines form a family of signaling proteins mainly responsible for directing the traffic of leukocytes, where their biological activity can be modulated by their oligomerization state. We characterize the dynamics and thermodynamic stability of monomer and homodimer structures of CXCL7, one of the most abundant platelet chemokines, using experimental methods that include circular dichroism (CD) and nuclear magnetic resonance (NMR) spectroscopy, and computational methods that include the anisotropic network model (ANM), molecular dynamics (MD) simulations and the distance constraint model (DCM). A consistent picture emerges for the effects of dimerization and Cys5‐Cys31 and Cys7‐Cys47 disulfide bonds formation. The presence of disulfide bonds is not critical for maintaining structural stability in the monomer or dimer, but the monomer is destabilized more than the dimer upon removal of disulfide bonds. Disulfide bonds play a key role in shaping the characteristics of native state dynamics. The combined analysis shows that upon dimerization flexibly correlated motions are induced between the 30s and 50s loop within each monomer and across the dimer interface. Interestingly, the greatest gain in flexibility upon dimerization occurs when both disulfide bonds are present, and the homodimer is least stable relative to its two monomers. These results suggest that the highly conserved disulfide bonds in chemokines facilitate a structural mechanism that is tuned to optimally distinguish functional characteristics between monomer and dimer. Proteins 2015; 83:1987–2007. © 2015 Wiley Periodicals, Inc.     

Effect of an alternative disulfide bond on the structure, stability, and folding of human lysozyme

Human lysozyme has four disulfide bonds, one of which, Cys65-Cys81, is included in a long loop of the beta-domain. A cysteine-scanning mutagenesis in which the position of Cys65 was shifted within a continuous segment from positions 61 to 67, with fixed Cys81, has previously shown that only the mutant W64CC65A, which has a nonnative Cys64-Cys81 disulfide, can be correctly folded and secreted by yeast. Here, using the W64CC65A mutant, we investigated the effects of an alternative disulfide bond on the structure, stability, and folding of human lysozyme using circular dichroism (CD) and fluorescence spectroscopy combined with a stopped-flow technique. Although the mutant is expected to have a different main-chain structure from that of the wild-type protein around the loop region, far- and near-UV CD spectra show that the native state of the mutant has tightly packed side chains and secondary structure similar to that of the wild-type. Guanidine hydrochloride-induced equilibrium unfolding transition of the mutant is reversible, showing high stability and cooperativity of folding. In the kinetic folding reaction, both proteins accumulate a similar burst-phase intermediate having pronounced secondary structure within the dead time of the measurement and fold into the native structure by means of a similar folding mechanism. Both the kinetic refolding and unfolding reactions of the mutant protein are faster than those of the wild-type, but the increase in the unfolding rate is larger than that of the refolding rate. The Gibbs' free-energy diagrams obtained from the kinetic analysis suggest that the structure around the loop region in the beta-domain of human lysozyme is formed after the transition state of folding, and thus, the effect of the alternative disulfide bond on the structure, stability, and folding of human lysozyme appears mainly in the native state.     

The disulfide structure of denatured epidermal growth factor: preparation of scrambled disulfide isomers

The conformational stability of human epidermal growth factor (EGF) and the structure of denatured EGF were investigated using the technique of disulfide scrambling. Under denaturing conditions and in the presence of a thiol catalyst, the native EGF denatures by shuffling its three native disulfide bonds and converts to a mixture of scrambled isomers. Analysis by HPLC reveals that the denatured EGF is composed of about 10 fractions of scrambled isomers. The heterogeneity varies under different denaturing conditions, with the heat-denatured samples exhibiting the highest degree of heterogeneity. The disulfide structures of eight major scrambled isomers of EGF were determined. The most predominant isomer adopts the bead-form structure with disulfide bonds bridged by three pairs of neighboring cysteines: Cys⁶-Cys¹⁴, Cys²⁰-Cys³¹, and Cys³³-Cys⁴². The denaturation curve of EGF is determined by the relative yield of the scrambled and native species of EGF. EGF is a highly stable molecule and can be effectively denatured only by guanidine chloride at a concentration of greater than 4–5 M. At 8 M urea, less than 16% of the native EGF was denatured. The unusual conformational stability of EGF was compared with that of eight different disulfide proteins that were similarly characterized by the method of disulfide scrambling.     

Disulfide bonds and protein folding

The applications of disulfide-bond chemistry to studies of protein folding, structure, and stability are reviewed and illustrated with bovine pancreatic ribonuclease A (RNase A). After surveying the general properties and advantages of disulfide-bond studies, we illustrate the mechanism of reductive unfolding with RNase A, and discuss its application to probing structural fluctuations in folded proteins. The oxidative folding of RNase A is then described, focusing on the role of structure formation in the regeneration of the native disulfide bonds. The development of structure and conformational order in the disulfide intermediates during oxidative folding is characterized. Partially folded disulfide species are not observed, indicating that disulfide-coupled folding is highly cooperative. Contrary to the predictions of "rugged funnel" models of protein folding, misfolded disulfide species are also not observed despite the potentially stabilizing effect of many nonnative disulfide bonds. The mechanism of regenerating the native disulfide bonds suggests an analogous scenario for conformational folding. Finally, engineered covalent cross-links may be used to assay for the association of protein segments in the folding transition state, as illustrated with RNase A.     

Unfolding and pH studies on manganese peroxidase: role of heme and calcium on secondary structure stability

The present study characterizes the unfolding and folding processes of recombinant manganese peroxidase. This enzyme contains five disulfide bonds, two calcium ions, and one heme prosthetic group. Circular dichroism in the far UV was used to monitor global changes of the protein secondary structure, whereas UV-visible spectroscopy of the Soret band provided information about local changes in the heme cavity. The effects of reducing agents, oxidizing agents, and denaturants on this process were investigated. In addition to affecting the secondary structure content, these factors also affect the binding of the heme and the calcium ions, both of which have a significant effect on the folding process. Our results also show that denaturants induce irreversible changes, which are most likely due to the inability of the denatured protein to rebind either calcium or the heme. Breaking of disulfide bonds by 30 mM dithiothreitol causes complete unfolding of recombinant manganese peroxidase. The unfolding process was also studied at low and high pH, where the protein reaches the final unfolded state through two different intermediate states. The data also indicate that only the acidic folding-unfolding process is reversible. Our results indicate a complex synergistic relationship between the secondary structure content, the tertiary structure arrangement, and the binding of the heme and the calcium ions and disulfide bridge formation.     

Structural and thermodynamic consequences of removal of a conserved disulfide bond from equine beta-lactoglobulin

A disulfide bond between cysteine 66 and cysteine 160 of equine beta-lactoglobulin was removed by substituting cysteine residues with alanine. This disulfide bond is conserved across the lipocalin family. The conformation and stability of the disulfide-deleted mutant protein was investigated by circular dichroism. The mutant protein assumes a native-like structure under physiological conditions and assumes a helix-rich molten globule structure at acid pH or at moderate concentrations of urea as the wild-type protein does. The urea-induced unfolding experiment shows that the stability of the native conformation was reduced but that of the molten globule intermediate is not significantly changed at pH 4 by removal of the disulfide bond. On the other hand, the molten globule at acid pH was destabilized by removal of the disulfide bond. This difference in the stabilizing effect of the disulfide bond was interpreted by the effect of the disulfide in keeping the molecule compact against the electrostatic repulsion at acid pH. In contrast to the wild-type protein, the circular dichroism spectrum in the molten globule state at acid pH depends on anion concentration, suggesting that the expansion of the molecule through electrostatic repulsion induces alpha-helices as observed in the cold denatured state of the wild-type protein.     

Effects of disulfide bonds on folding behavior and mechanism of the beta-sheet protein tendamistat

Tendamistat, a small disulfide-bonded beta-sheet protein, and its three single/double-disulfide mutants are investigated by using a modified Gō-like model, aiming to understand the folding mechanism of disulfide-bonded protein as well as the effects of removal of disulfide bond on the folding process. Our simulations show that tendamistat and its two single-disulfide mutants are all two-state folders, consistent with the experimental observations. It is found that the disulfide bonds as well as three hydrogen bonds between the N-terminal loop-0 and strand-6 are of significant importance for the folding of tendamistat. Without these interactions, their two-state behaviors become unstable and the predictions of the model are inconsistent with experiments. In addition, the effect of disulfide bonds on the folding process are studied by comparing the wild-type tendamistat and its two mutants; it is found that the removal of either of the C11-C27 or C45-C73 disulfide bond leads to a large decrease in the thermodynamical stability and loss of structure in the unfolded state, and the effect of the former is stronger than that of the later. These simulation results are in good agreement with experiments and, thus, validate our model. Based on the same model, the detailed folding pathways of the wild-type tendamistat and two mutants are studied, and the effect of disulfide bonds on the folding kinetics are discussed. The obtained results provide a detailed folding picture of these proteins and complement experimental findings. Finally, the folding nuclei predicted to be existent in this protein tendamistat as well as its mutants are firstly identified in this work. The positions of the nucleus are consistent with those argued in experimental studies. Therefore, a nucleation/growth folding mechanism that can explain the two-state folding manner is clearly characterized. Moreover, the effect by the removal of each disulfide bond on the folding thermodynamics and dynamics can also be well interpreted from their influence on the folding nucleus. The implementation of this work indicates that the modified Gō-like model really describes the folding behavior of protein tendamistat and could be used to study the folding of other disulfide-bonded proteins.     

The structure of denatured alpha-lactalbumin elucidated by the technique of disulfide scrambling: fractionation of conformational isomers of alpha-lactalbumin

The structure of denatured alpha-lactalbumin (alpha-LA) has been characterized using the method of disulfide scrambling. Under denaturing conditions (urea, guanidine hydrochloride, guanidine thiocyanate, organic solvent or elevated temperature) and in the presence of thiol initiator, alpha-LA denatures by shuffling its four native disulfide bonds and converts to a mixture of fully oxidized scrambled structures. Analysis by reversed-phase HPLC reveals that the denatured alpha-LA comprises a minimum of 45 fractions of scrambled isomers. Among them, six well populated isomers have been isolated and structurally characterized. Their relative concentrations, which represent the fingerprinting of the denatured alpha-LA, vary substantially under different denaturing conditions. These results permit independent plotting of the denaturation and unfolding curves of alpha-LA. Most importantly, unique isomers of partially unfolded alpha-LA were shown to populate at mild and selected denaturing conditions. Organic solvent disrupts preferentially the hydrophobic alpha-helical domain, generating a predominant isomer containing two native disulfide bonds at the beta-sheet domain and two scrambled disulfide bonds at the alpha-helical region. Thermal denaturation selectively unfolds the beta-sheet domain of alpha-LA, producing a prevalent isomer that exhibits structural characteristics of the molten globule state of alpha-LA.     

Uncovering Specific Electrostatic Interactions in the Denatured States of Proteins

The stability and folding of proteins are modulated by energetically significant interactions in the denatured state that is in equilibrium with the native state. These interactions remain largely invisible to current experimental techniques, however, due to the sparse population and conformational heterogeneity of the denatured-state ensemble under folding conditions. Molecular dynamics simulations using physics-based force fields can in principle offer atomistic details of the denatured state. However, practical applications are plagued with the lack of rigorous means to validate microscopic information and deficiencies in force fields and solvent models. This study presents a method based on coupled titration and molecular dynamics sampling of the denatured state starting from the extended sequence under native conditions. The resulting denatured-state pK_as allow for the prediction of experimental observables such as pH- and mutation-induced stability changes. I show the capability and use of the method by investigating the electrostatic interactions in the denatured states of wild-type and K12M mutant of NTL9 protein. This study shows that the major errors in electrostatics can be identified by validating the titration properties of the fragment peptides derived from the sequence of the intact protein. Consistent with experimental evidence, our simulations show a significantly depressed pK_a for Asp⁸ in the denatured state of wild-type, which is due to a nonnative interaction between Asp⁸ and Lys¹². Interestingly, the simulation also shows a nonnative interaction between Asp⁸ and Glu⁴⁸ in the denatured state of the mutant. I believe the presented method is general and can be applied to extract and validate microscopic electrostatics of the entire folding energy landscape.