pK(a) values for the unfolded state under native conditions explain the pH-dependent stability of PGB1

Understanding the role of electrostatics in protein stability requires knowledge of these interactions in both the folded and unfolded states. Electrostatic interactions can be probed experimentally by characterizing ionization equilibria of titrating groups, parameterized as pK_a values. However, pK_a values of the unfolded state are rarely accessible under native conditions, where the unfolded state has a very low population. Here, we report pK_a values under nondenaturing conditions for two unfolded fragments of the protein G B1 domain that mimic the unfolded state of the intact protein. pK_a values were determined for carboxyl groups by monitoring their pH-dependent ¹³C chemical shifts. Monte Carlo simulations using a Gaussian chain model provide corrections for changes in electrostatic interactions that arise from fragmentation of the protein. Most pK_a values for the unfolded state agree well with model values, but some residues show significant perturbations that can be rationalized by local electrostatic interactions. The pH-dependent stability was calculated from the experimental pK_a values of the folded and unfolded states and compared to experimental stability data. The use of experimental pK_a values for the unfolded state results in significantly improved agreement with experimental data, as compared to calculations based on model data alone.     

The effects of pK(a) tuning on the thermodynamics and kinetics of folding: design of a solvent-shielded carboxylate pair at the a-position of a coiled-coil

The tuning of the pK_a of ionizable residues plays a critical role in various protein functions, such as ligand-binding, catalysis, and allostery. Proteins harness the free energy of folding to position ionizable groups in highly specific environments that strongly affect their pK_a values. To investigate the interplay among protein folding kinetics, thermodynamics, and pK_a modulation, we introduced a pair of Asp residues at neighboring interior positions of a coiled-coil. A single Asp residue was replaced for an Asn side chain at the a-position of the coiled-coil from GCN4, which was also crosslinked at the C-terminus via a flexible disulfide bond. The thermodynamic and kinetic stability of the system was measured by circular dichroism and stopped-flow fluorescence as a function of pH and concentration of guanidine HCl. Both sets of data are consistent with a two-state equilibrium between fully folded and unfolded forms. Distinct pK_a values of 6.3 and 5.35 are assigned to the first and second protonation of the Asp pair; together they represent an energetic difference of 5 kcal/mol relative to the protonation of two Asp residues with unperturbed pK_a values. Analysis of the rate data as a function of pH and denaturant concentration allowed calculation of the kinetic constants for the conformational transitions of the peptide with the Asp residues in the doubly protonated, singly protonated, and unprotonated forms. The doubly and singly protonated forms fold rapidly, and a ϕ-value analysis shows that their contribution to folding occurs subsequent to the transition state ensemble for folding. By contrast, the doubly charged state shows a reduced rate of folding and a ϕ-value near 0.5 indicative of a repulsive interaction, and possibly also heterogeneity in the transition state ensemble.     

pH-Dependent Conformational Changes in Proteins and Their Effect on Experimental pKas: The Case of Nitrophorin 4

The acid-base behavior of amino acids is an important subject of study due to their prominent role in enzyme catalysis, substrate binding and protein structure. Due to interactions with the protein environment, their pK_as can be shifted from their solution values and, if a protein has two stable conformations, it is possible for a residue to have different “microscopic”, conformation-dependent pK_a values. In those cases, interpretation of experimental measurements of the pK_a is complicated by the coupling between pH, protonation state and protein conformation. We explored these issues using Nitrophorin 4 (NP4), a protein that releases NO in a pH sensitive manner. At pH 5.5 NP4 is in a closed conformation where NO is tightly bound, while at pH 7.5 Asp30 becomes deprotonated, causing the conformation to change to an open state from which NO can easily escape. Using constant pH molecular dynamics we found two distinct microscopic Asp30 pK_as: 8.5 in the closed structure and 4.3 in the open structure. Using a four-state model, we then related the obtained microscopic values to the experimentally observed “apparent” pK_a, obtaining a value of 6.5, in excellent agreement with experimental data. This value must be interpreted as the pH at which the closed to open population transition takes place. More generally, our results show that it is possible to relate microscopic structure dependent pKa values to experimentally observed ensemble dependent apparent pK_as and that the insight gained in the relatively simple case of NP4 can be useful in several more complex cases involving a pH dependent transition, of great biochemical interest.     

Electrostatics of folded and unfolded bovine β-lactoglobulin

We report on electrophoretic, spectroscopic, and computational studies aimed at clarifying, at atomic resolution, the electrostatics of folded and unfolded bovine β-lactoglobulin (BLG) with a detailed characterization of the specific aminoacids involved. The procedures we used involved denaturant gradient gel electrophoresis, isoelectric focusing, electrophoretic titration curves, circular dichroism and fluorescence spectra in the presence of increasing concentrations of urea (up to 8 M), electrostatics computations and low-mode molecular dynamics. Discrepancy between electrophoretic and spectroscopic evidence suggests that changes in mobility induced by urea are not just the result of changes in gyration radius upon unfolding. Electrophoretic titration curves run across a pH range of 3.5–9 in the presence of urea suggest that more than one aminoacid residue may have anomalous pK _a value in native BLG. Detailed computational studies indicate a shift in pKa of Glu44, Glu89, and Glu114, mainly due to changes in global and local desolvation. For His161, the formation of hydrogen bond(s) could add up to desolvation contributions. However, since His161 is at the C terminus, the end-effect associated to the solvated form strongly influences its pK _a value with extreme variation between crystal structures on one side and NMR or low-mode molecular dynamics structures on the other. The urea concentration effective in BLG unfolding depends on pH, with higher stability of the protein at lower pH.     

Improving the accuracy of protein pKa calculations: Conformational averaging versus the average structure

Several methods for including the conformational flexibility of proteins in the calculation of titration curves are compared. The methods use the linearized Poisson-Boltzmann equation to calculate the electrostatic free energies of solvation and are applied to bovine pancreatic trypsin inhibitor (BPTI) and hen egg-white lysozyme (HEWL). An ensemble of conformations is generated by a molecular dynamics simulation of the proteins with explicit solvent. The average titration curve of the ensemble is calculated in three different ways: an average structure is used for the pK_a calculation; the electrostatic interaction free energies are averaged and used for the pK_a calculation; and the titration curve for each structure is calculated and the curves are averaged. The three averaging methods give very similar results and improve the pK_a values to approximately the same degree. This suggests, in contrast to implications from other work, that the observed improvement of pK_a values in the present studies is due not to averaging over an ensemble of structures, but rather to the generation of a single properly averaged structure for the pK_a calculation. Proteins 33:145–158, 1998. © 1998 Wiley-Liss, Inc.     

Mechanism of Long-Chain Free Fatty Acid Protonation at the Membrane-Water Interface

Long-chain free fatty acids (FFAs) play an important role in several physiological and pathological processes such as lipid fusion, adjustments of membrane permeability and fluidity, and the regulation of enzyme and protein activities. FFA-facilitated membrane proton transport (flip-flop) and FFA-dependent proton transport by membrane proteins (e.g., mitochondrial uncoupling proteins) are governed by the difference between FFA’s intrinsic pK_a value and the pH in the immediate membrane vicinity. Thus far, a quantitative understanding of the process has been hampered, because the pK_a value shifts upon moving the FFA from the aqueous solution into the membrane. For the same FFA, pK_a values between 5 and 10.5 were reported. Here, we systematically evaluated the dependence of pK_a values on chain length and number of double bonds by measuring the ζ-potential of liposomes reconstituted with FFA at different pH values. The experimentally obtained intrinsic pK_a values (6.25, 6.93, and 7.28 for DOPC membranes) increased with FFA chain length (C16, C18, and C20), indicating that the hydrophobic energy of transfer into the bilayer is an important pK_a determinant. The observed pK_a decrease in DOPC with increasing number of FFA double bonds (7.28, 6.49, 6.16, and 6.13 for C20:0, C20:1, C20:2, and C20:4, respectively) is in line with a decrease in transfer energy. Molecular dynamic simulations revealed that the ionized carboxylic group of the FFAs occupied a fixed position in the bilayer independent of chain length, underlining the importance of Born energy. We conclude that pK_a is determined by the interplay between the energetic costs for 1) burying the charged moiety into the lipid bilayer and 2) transferring the hydrophobic protonated FFA into the bilayer.     

The transient accumulation of the signaling state of photoactive yellow protein is controlled by the external pH

The signaling state of the photoreceptor photoactive yellow protein is the long-lived intermediate I₂′. The pH dependence of the equilibrium between the transient photocycle intermediates I₂ and I₂′ was investigated. The formation of I₂′ from I₂ is accompanied by a major conformational change. The kinetics and intermediates of the photocycle and of the photoreversal were measured by transient absorption spectroscopy from pH 4.6 to 8.4. Singular value decomposition (SVD) analysis of the data at pH 7 showed the presence of three spectrally distinguishable species: I₁, I₂, and I₂′. Their spectra were determined using the extrapolated difference method. I₂ and I₂′ have electronic absorption spectra, with maxima at 370 ± 5 and 350 ± 5 nm, respectively. Formation of the signaling state is thus associated with a change in the environment of the protonated chromophore. The time courses of the I₁, I₂, and I₂′ intermediates were determined from the wavelength-dependent transient absorbance changes at each pH, assuming that their spectra are pH-independent. After the formation of I₂′ (~2 ms), these three intermediates are in equilibrium and decay together to the initial dark state. The equilibrium between I₂ and I₂′ is pH dependent with a pK_a of 6.4 and with I₂′ the main species above this pK_a. Measurements of the pH dependence of the photoreversal kinetics with a second flash of 355 nm at a delay of 20 ms confirm this pK_a value. I₂ and I₂′ are photoreversed with reversal times of ~55 μs and several hundred microseconds, respectively. The corresponding signal amplitudes are pH dependent with a pK_a of ~6.1. Photoreversal from I₂′ dominates above the pK_a. The transient accumulation of I₂′, the active state of photoactive yellow protein, is thus controlled by the proton concentration. The rate constant k₃ for the recovery to the initial dark state also has a pK_a of ~6.3. This equality of the equilibrium and kinetic pK_a values is not accidental and suggests that k₃ is proportional to [I₂′].     

Different molten globule-like folding intermediates of hen egg white lysozyme induced by high pH and tertiary butanol

We have provided evidence that hen egg white lysozyme (HEWL) existed in alpha helical and beta structure dominated molten globule (MG) states at high pH and in the presence of tertiary butanol, respectively. Circular dichroism (CD), intrinsic fluorescence, ANS binding and acrylamide-induced fluorescence quenching techniques have been used to investigate alkali-induced unfolding of HEWL and the effect of tertiary butanol on the alkaline-induced state. At pH 12.75, HEWL existed as molten globule like intermediate. The observed MG-like intermediate was characterized by (i) retention of 77% of the native secondary structure, (ii) enhanced binding of ANS (approximately 5 times) compared to native and completely unfolded state, (iii) loss of the tertiary structure as indicated by the tertiary structural probes (near-UV, CD and Intrinsic fluorescence) and (iv) acrylamide quenching studies showed that MG state has compactness intermediate between native and completely unfolded states. Moreover, structural properties of the protein at isoelectric point (pI) and denatured states have also been described. We have also shown that in the presence of 45% tertiary butanol (t-butanol), HEWL at pH 7.0 and 11.0 (pI 11.0) existed in helical structure without much affecting tertiary structure. Interestingly, MG state of HEWL at pH 12.7 transformed into another MG state (MG2) at 20% t-butanol (v/v), in which secondary structure is mainly beta sheets. On further increasing the t-butanol concentration alpha helix was found to reform. We have proposed that formation of both alpha helical and beta sheet dominated intermediate may be possible in the folding pathway of alpha + beta protein.     

Mechanism of acridine uptake by submitochondrial particles.

Acridines were compared regarding their ability to be taken up by submitochondrial particles under energized conditions. pH dependence of uptake was explored, and it was found that acridines fell into three classes independently of their pK_a value: acridines which are not taken up, acridines taken up at all pH values, and acridines taken up only at alkaline pH. Partition measurements between H₂O and chloroform phase showed a similar pattern, and affinity for the organic phase seemed to parallel uptake. Acridines which are taken up by submitochondrial particles at acidic pH under energized conditions despite a high pK_a value could also be extracted into chloroform at acidic pH, thus implying that the dye's uncharged form has a high affinity for the organic phase. By supplementing the aqueous medium with lipophilic anions, the dye may also be extracted in its charged form. The data support a mechanism for acridine uptake in which diffusion of the uncharged form across the membranes is an obligatory step. Some previously reported inhibitory anion effects on uptake may be explained by ion pair formation, which allows release of the accumulated charged form.     

Backbone relaxation coupled to the ionization of internal groups in proteins: a self-guided Langevin dynamics study

Pathways of structural relaxation triggered by ionization of internal groups in staphylococcal nuclease (SNase) were studied through multiple self-guided Langevin dynamics (SGLD) simulations. Circular dichroism, steady-state Trp fluorescence, and nuclear magnetic resonance spectroscopy have shown previously that variants of SNase with internal Glu, Asp, and Lys at positions 66 or 92, and Arg at position 66, exhibit local reorganization or global unfolding when the internal ionizable group is charged. Except for Arg-66, these internal ionizable groups have unusual pK_a values and are neutral at physiological pH. The structural trends observed in the simulations are in general agreement with experimental observations. The I92D variant, which unfolds globally upon ionization of Asp-92, in simulations often exhibits extensive hydration of the protein core, and sometimes also significant perturbations of the β-barrel. In the crystal structure of the V66R variant, the β1 strand from the β-barrel is domain-swapped; in the simulations, the β1 strand is sometimes partially released. The V66K variant, which in solutions shows reorganization of six residues at the C-terminus of helix α1 and perturbations in the β-barrel structure, exhibits fraying of three residues of helix α1 in one simulation, and perturbations and partial unfolding of three β-strands in a few other simulations. In sharp contrast, very small structural changes were observed in simulations of the wild-type protein. The simulations indicate that charging of internal groups frequently triggers penetration of water into the protein interior. The pK_a values of Asp-92 and Arg-66 calculated with continuum methods on SGLD-relaxed structures reached the normal values in most simulations. Detailed analysis of accuracy and performance of SGLD demonstrates that SGLD outperforms LD in sampling of alternative protein conformations without loss of the accuracy and level of detail characteristic of regular LD.     

Prediction of Secondary Ionization of the Phosphate Group in Phosphotyrosine Peptides

A computational approach, based on a continuum molecular electrostatics model, for the calculation of the pK_a values of secondary ionization of the phosphate group in phenyl phosphate derivatives is described. The method uses the ESP atomic charges of the mono-anionic and di-anionic forms of the ionizable phosphate group, computed with the use of the density functional method, and applies a new concept of the model group, being the reference state for the pK_a calculations. Both conformational flexibility and tautomeric degrees of freedom are taken into account in the calculations. The method was parameterized using experimentally available pK_a values of four derivatives of phenyl phosphates, and phosphotyrosine. Subsequently this parameterization was used to predict pK_a of the phosphate group in a short peptide Gly-Gly-Tyr(P)-Ala, and in a longer peptide consisting of 12 residues, the latter in water, and in a complex with a protein—phospholipase. The agreement between the computed and the experimental pK_a values is better than ±0.3 pH units for the optimized solute dielectric constant of 11–13. This approach is promising and its extension to other phospho-amino acids is in progress.     

Rates of structural fluctuations of lysozyme in the range of thermal unfolding transition

The hydrogen–deuterium exchange reaction for the tryptophan residues in lysozyme have been followed in 4.5M LiBr at pH 7.2 in the temperature range of the unfolding transition by measuring the transmittance change at 293 nm. The exchange reaction proceeded in three phases at low temperature for native protein. The first and the second phases were ascribed to the H-D exchange reactions of three relatively exposed tryptophan residues on the molecular surface. The third phase corresponded to the H-D exchange reaction of the three tryptophan residues buried in the interior of the molecule. The H-D exchange reaction proceeded in two phases near the melting temperature and in a single phase at high temperature, where almost all molecules are unfolded. The H-D exchange of three tryptophan residues buried in folded molecules was caused by fluctuation between the folded and unfolded structure of the protein molecule. The rates of such a fluctuation were determined from the rates of the exchange reaction at various temperatures. These rates agreed very well with those determined from the temperature-jump method. This means that a protein molecule in solution fluctuates between the N- and D-states at every temperature within the transition region, where the N-form is the tightly folded native structure and the D-form the randomly coiled chain. From measurements of thermal unfolding of ester-108-lysozyme and the binding constant of (NAG)₃ to ester-108-lysozyme, it was found that almost all cross-linked molecules are in the folded state near 50°C and pH 7.2 in 4.5M LiBr, where intact molecules are unfolded. We also studied the H-D exchange reaction of ester-108-lysozyme. In the temperature region of 43–50°C, about 70% of the exchangeable tryptophan residues of ester-108-lysozyme were exchanged within 1 s immediately after the mixing of D₂O, in spite of the fact that almost all molecules are in the folded state. This was considered the premelting of the surface of a corss-linked molecule.     

Realistic modeling of the denatured states of proteins allows accurate calculations of the pH dependence of protein stability

Computational techniques based on continuum electrostatics treatments have been successful in predicting and interpreting the pKa values of ionizable amino acids in folded proteins. Despite this progress, efforts to reproduce the pH-dependence of protein stability have met with only limited success: agreement with experimental results has been only qualitative. It has been argued previously that the most likely reason for discrepancies is the presence of residual electrostatic interactions in the unfolded state, which cause pKa values to be shifted from their model compound values. Here we show that by constructing atomistic models of the unfolded state with a simple molecular mechanics protocol that uses the native state as a starting point, much improved reproduction of pH effects on protein stability can be obtained. In contrast, when a fully extended model of the unfolded state is used, no such improvement is obtained, a result that suggests that local interactions with residues nearby in the sequence are not sufficient to properly account for the pKa shifts in the unfolded state. In comparison to model compound values, the pKa values of acidic residues in "native-like" unfolded states are typically found to be shifted downwards by approximately 0.3 pH unit, in good agreement with the average downward shift deduced from experimental measurements. Given its success in the present situation, the protocol employed here for developing simple models of the unfolded state may prove useful in other computer simulation applications.     

Rheostat Re-Wired: Alternative Hypotheses for the Control of Thioredoxin Reduction Potentials

Thioredoxins are small soluble proteins that contain a redox-active disulfide (CXXC). These disulfides are tuned to oxidizing or reducing potentials depending on the function of the thioredoxin within the cell. The mechanism by which the potential is tuned has been controversial, with two main hypotheses: first, that redox potential (E_m) is specifically governed by a molecular ‘rheostat’—the XX amino acids, which influence the Cys pK_a values, and thereby, E_m; and second, the overall thermodynamics of protein folding stability regulates the potential. Here, we use protein film voltammetry (PFV) to measure the pH dependence of the redox potentials of a series of wild-type and mutant archaeal Trxs, PFV and glutathionine-equilibrium to corroborate the measured potentials, the fluorescence probe BADAN to measure pK_a values, guanidinium-based denaturation to measure protein unfolding, and X-ray crystallography to provide a structural basis for our functional analyses. We find that when these archaeal thioredoxins are probed directly using PFV, both the high and low potential thioredoxins display consistent 2H⁺:2e^- coupling over a physiological pH range, in conflict with the conventional ‘rheostat’ model. Instead, folding measurements reveals an excellent correlation to reduction potentials, supporting the second hypothesis and revealing the molecular mechanism of reduction potential control in the ubiquitous Trx family.     

Direct observation of multiple protonation states in recombinant human purple acid phosphatase

To date, most spectroscopic studies on mammalian purple acid phosphatases (PAPs) have been performed at a single pH, typically pH 5. The catalytic activity of these enzymes is, however, pH dependent, with optimal pH values of 5.5–6.2 (depending on the form). For example, the pH optimum of PAPs isolated as single polypeptides is around pH 5.5, which is substantially lower that of proteolytically cleaved PAPs (ca. pH 6.2). In addition, the catalytic activity of single polypeptide PAPs at their optimal pH values is four to fivefold lower than that of the proteolytically cleaved enzymes. In order to elucidate the chemical basis for the pH dependence of these enzymes, the spectroscopic properties of both the single polypeptide and proteolytically cleaved forms of recombinant human PAP (recHPAP) and their complexes with inhibitory anions have been examined over the pH range 4 to 8. The EPR spectra of both forms of recHPAP are pH dependent and show the presence of three species: an inactive low pH form (pH<pK _a,1), an active form (pK _a,1<pH<pK _a,2), and an inactive high pH form (pH>pK _a,2). The pK _a,1 values observed by EPR for the single polypeptide and proteolytically cleaved forms are similar to those previously observed in kinetics studies. The spectroscopic properties of the enzyme–phosphate complex (which should mimic the enzyme–substrate complex), the enzyme–fluoride complex, and the enzyme–fluoride–phosphate complex (which should mimic the ternary enzyme–substrate–hydroxide complex) were also examined. EPR spectra show that phosphate binds to the diiron center of the proteolytically cleaved form of the enzyme, but not to that of the single polypeptide form. EPR spectra also show that fluoride binds only to the low pH form of the enzymes, in which it presumably replaces a coordinated water molecule. The binding of fluoride and phosphate to form a ternary complex appears to be cooperative.Electronic Supplementary Material Supplementary material is available for this article at     

Electrostatic contributions to the stability of the GCN4 leucine zipper structure

Ion pairs are ubiquitous in X-ray structures of coiled coils, and mutagenesis of charged residues can result in large stability losses. By contrast, pK_a values determined by NMR in solution often predict only small contributions to stability from charge interactions. To help reconcile these results we used triple-resonance NMR to determine pK_a values for all groups that ionize between pH 1 and 13 in the 33 residue leucine zipper fragment, GCN4p. In addition to the native state we also determined comprehensive pK_a values for two models of the GCN4p denatured state: the protein in 6 M urea, and unfolded peptide fragments of the protein in water. Only residues that form ion pairs in multiple X-ray structures of GCN4p gave large pK_a differences between the native and denatured states. Moreover, electrostatic contributions to stability were not equivalent for oppositely charged partners in ion pairs, suggesting that the interactions between a charge and its environment are as important as those within the ion pair. The pH dependence of protein stability calculated from NMR-derived pK_a values agreed with the stability profile measured from equilibrium urea-unfolding experiments as a function of pH. The stability profile was also reproduced with structure-based continuum electrostatic calculations, although contributions to stability were overestimated at the extremes of pH. We consider potential sources of errors in the calculations, and how pK_a predictions could be improved. Our results show that although hydrophobic packing and hydrogen bonding have dominant roles, electrostatic interactions also make significant contributions to the stability of the coiled coil.     

Electrostatic contribution of surface charge residues to the stability of a thermophilic protein: benchmarking experimental and predicted pKa values

Optimization of the surface charges is a promising strategy for increasing thermostability of proteins. Electrostatic contribution of ionizable groups to the protein stability can be estimated from the differences between the pKa values in the folded and unfolded states of a protein. Using this pKa-shift approach, we experimentally measured the electrostatic contribution of all aspartate and glutamate residues to the stability of a thermophilic ribosomal protein L30e from Thermococcus celer. The pKa values in the unfolded state were found to be similar to model compound pKas. The pKa values in both the folded and unfolded states obtained at 298 and 333 K were similar, suggesting that electrostatic contribution of ionizable groups to the protein stability were insensitive to temperature changes. The experimental pKa values for the L30e protein in the folded state were used as a benchmark to test the robustness of pKa prediction by various computational methods such as H++, MCCE, MEAD, pKD, PropKa, and UHBD. Although the predicted pKa values were affected by crystal contacts that may alter the side-chain conformation of surface charged residues, most computational methods performed well, with correlation coefficients between experimental and calculated pKa values ranging from 0.49 to 0.91 (p<0.01). The changes in protein stability derived from the experimental pKa-shift approach correlate well (r = 0.81) with those obtained from stability measurements of charge-to-alanine substituted variants of the L30e protein. Our results demonstrate that the knowledge of the pKa values in the folded state provides sufficient rationale for the redesign of protein surface charges leading to improved protein stability.     

Thermodynamics and Kinetics of Single-Chain Monellin Folding with Structural Insights into Specific Collapse in the Denatured State Ensemble

Proteins, which behave as random coils in high denaturant concentrations undergo collapse transition similar to polymers on denaturant dilution. We study collapse in the denatured ensemble of single-chain monellin (MNEI) using a coarse-grained protein model and molecular dynamics simulations. The model is validated by quantitatively comparing the computed guanidinium chloride and pH-dependent thermodynamic properties of MNEI folding with the experiments. The computed properties such as the fraction of the protein in the folded state and radius of gyration (R_g) as function of [GuHCl] are in good agreement with the experiments. The folded state of MNEI is destabilized with an increase in pH due to the deprotonation of the residues Glu24 and Cys42. On decreasing [GuHCl], the protein in the unfolded ensemble showed specific compaction. The R_g of the protein decreased steadily with [GuHCl] dilution due to increase in the number of native contacts in all the secondary structural elements present in the protein. MNEI folding kinetics is complex with multiple folding pathways and transiently stable intermediates are populated in these pathways. In strong stabilizing conditions, the protein in the unfolded ensemble showed transition to a more compact unfolded state where R_g decreased by ≈ 17% due to the formation of specific native contacts in the protein. The intermediate populated in the dominant MNEI folding pathway satisfies the structural features of the dry molten globule inferred from experiments.     

Probing residue-level unfolding during lysozyme precipitation

We have employed nuclear magnetic resonance (NMR) measurements of hydrogen exchange to identify residue-level conformational changes in hen egg white lysozyme (HEWL) as induced by salt precipitation. Deuterated HEWL was dissolved into a phosphate (H2O) buffer and precipitated at pH 2.1 upon addition of solid KSCN or (ND4)2SO4, allowing isotope labeling of unfolded regions. After 1 h, each precipitate was then dissolved at pH 3.8 to initiate refolding and preserve labeling and subsequently purified for NMR analysis. HEWL precipitated by 1.0 M KSCN exhibited increased hydrogen exchange at 14 residues out of 42 normally well-protected in the native state. Of the affected residues, 9 were situated in the beta-sheet/loop domain. A similar, though less extensive, effect was observed at 0.2 M KSCN. Precipitation by 1.2 M (ND4)2SO4 resulted in none of the changes detected with KSCN. The popularity of ammonium sulfate as a precipitant is thus supported by this observed preservation of structural integrity. KSCN, in comparison, produced partial unfolding of specific regions in HEWL due most likely to known preferential interactions between -SCN and proteins. The severity of unfolding increased with KSCN concentration such that, at 1.0 M KSCN, almost the entire beta-sheet/loop domain of HEWL was disrupted. Even so, a portion of the HEWL core encompassed by three alpha-helices remained intact, possibly facilitating precipitate dissolution.