many debilitating and often fatal disorders , including alzheimer s and parkinson s diseases , the transmissible spongiform encephalopathies , type ii diabetes , and a range of systemic amyloidoses , are associated with the deposition of normally soluble proteins as insoluble aggregates in various types of tissue.(1 ) although the sequences and native structures of the proteins involved in the different conditions vary greatly , the aggregates found in patients exhibit very similar biophysical properties .
these include their fibrillar nature and cross- structure , where -strands are oriented perpendicularly to the axis of the fibril.(1 ) in many cases the formation of amyloid fibrils requires the monomeric precursor protein to undergo a series of substantial conformational rearrangements .
a detailed characterization of the different species involved in these conformational transitions is crucial for developing an understanding of the mechanism behind misfolding diseases .
such species are also key therapeutic targets for developing drugs to treat or prevent amyloid - related diseases .
partially folded states , denatured states , disordered protein ensembles , and molten globules are possible precursors to amyloid fibrils or other aggregates and have therefore been the object of extensive research .
human lysozyme is a small , well - characterized globular protein with a native structure that can be divided into two domains : the  domain ( residues 1 to 38 , and 86 to 130 ) and the  domain ( residues 39 to 85 ) that primarily contain -helical and -sheet secondary structure , respectively(18 ) ( see figure 1a , b ) . in 1993 ,
pepys and co - workers reported that two variants of human lysozyme are responsible for a hereditary non - neuropathic systemic amyloidosis.(20 ) given the significant body of data on the structure , dynamics , and folding of both hen and human lysozymes , the link to human disease renders this protein an ideal system to probe the relationships between protein folding and aggregation .
the first two single - point mutations found to be linked to lysozyme amyloidosis were i56 t and d67h , and other amyloidogenic and non - amyloidogenic mutations have also been identified more recently ( for a review , see dumoulin et al.(21 ) ) .
the observation that only certain mutational variants form fibrils in vivo raises the question of which factors regulate the amyloidogenicity of this protein .
a comparative thermal denaturation study has shown that the two amyloidogenic variants i56 t and d67h are less stable than the wild - type protein ( wt)11abbreviations : wt , wild - type ; nmr , nuclear magnetic resonance ; hsqc , heteronuclear single quantum coherence ; uv , ultraviolet ; cd , circular dichroism ; dsc , differential scanning calorimetry . 
 , suggesting that a destabilization of the native state is , at least in part , responsible for the process of fibril formation .
indeed at ph 5.0 the midpoints of thermal unfolding of the i56 t and d67h variants were reported to be 12  2 c lower than that of the wt protein .
however , the naturally occurring and non - amyloidogenic t70n variant also shows a decrease in the unfolding temperature in vitro of approximately 4 c relative to the wt but is not observed to form detectable quantities of fibrils or to give rise to disease in individuals carrying this mutation .
in addition , under in vitro conditions , t70n does not form fibrils at rates comparable to the amyloidogenic variants .
it therefore appears that the disease - associated mutational variants need to be destabilized enough to allow detectable amyloid fibril formation to take place under physiological conditions but not so destabilized as to trigger clearance by the cell .
abbreviations : wt , wild - type ; nmr , nuclear magnetic resonance ; hsqc , heteronuclear single quantum coherence ; uv , ultraviolet ; cd , circular dichroism ; dsc , differential scanning calorimetry .
( a ) representation of the secondary structure of the native state of wt human lysozyme , indicating glycine residues ( using dssp(17 ) on pdb entry 1lz1(18 ) ) .
strands are represented by arrows , -helices ( a to d ) by large boxes , and 310 helices ( h1 to h3 ) by small boxes .
the disulfide bridges c6c128 , c30c116 , c65c81 , and c77c95 are shown in purple .
( b ) tertiary structure of wt human lysozyme with the same color coding as in panel a. all structures are drawn using vmd(19 ) with the secondary structure defined by dssp.(17 ) vmd only represents strands comprising three residues and more ; therefore , strand 3 ( residues 5960 ) is not shown .
the colored spheres represent residues 56 ( blue ) and 59 ( red ) , sites of the mutations i56 t and i59
( c ) normalized signal of the near - uv circular dichroism thermal unfolding traces obtained for the wt ( black ) , i59 t ( red ) and i56 t ( blue ) variants , at ph 1.2 .
hydrogen / deuterium ( h / d ) exchange experiments on the amyloidogenic variants i56 t and d67h have revealed the presence of a partially folded intermediate state that is transiently populated under native conditions and characterized by the unfolding of the -domain and one of the helices ( helix c of the -domain ) .
in addition , thermal unfolding experiments carried out in vitro under equilibrium conditions in the presence of a fluorescent dye sensitive to hydrophobic patches ( 1-anilinonaphthalene-8-sulfonic acid or 
ans ) indicate the presence of partially unfolded species near to the midpoint of thermal denaturation .
on the basis of these results , a reduction in the global cooperativity of the variants of lysozyme found in disease has emerged as a key determinant for their amyloidogenicity : the non - cooperative unfolding behavior leads to the formation of partially folded species that are able to associate to form highly organized fibrils , the deposition of which is characteristic of amyloid diseases . despite the importance of partially folded states in this mechanism ,
their transient nature and low populations at equilibrium under solution conditions commonly used in in vitro studies ( ph > 5 ) have precluded a direct characterization of these species at high resolution . in this paper , we use low ph conditions to provide such high - resolution structural evidence on the nature of the equilibrium thermal unfolding behavior of human lysozyme . using nmr spectroscopy in combination with a wide range of biophysical techniques , we studied three variants of human lysozyme that have variable stabilities and propensities to aggregate in vitro , namely the wt protein and the i56 t and i59 t mutants ( see figure 1).(29 ) as stated above , the i56 t variant is linked to familial lysozyme systemic amyloidosis , and the i59 t mutant has been designed to have properties that are intermediate between those of i56 t and the wt lysozymes .
our results reveal that under the acidic destabilizing conditions commonly used to form amyloid fibrils in vitro , human lysozyme is in a dynamic equilibrium between its native conformation and a broad ensemble of non - native conformations : the members of this non - native  denatured ensemble contain both molten globular and unfolded regions .
molten globular regions , in which secondary structure elements are present but tertiary interactions are disorganized to varying degrees , are predominant at lower temperatures and unfold gradually , and with local cooperativity , as they become increasingly disordered at higher temperatures .
the strategy used in this study exploits low ph conditions to destabilize the native state and to observe partially unfolded states of lysozyme under equilibrium conditions.(30 ) lowering the ph from 5.0 to 1.2 results in a substantial decrease in the stability of the native state of all three variants , thereby enabling a comparison of the progressive denaturation of the wt protein and the i59 t and i56 t variants at temperatures below 50 c ( see figure 1c ) .
hn backbone hsqc spectra recorded at 20 c are typical of a folded protein , confirming that the native structure is preserved at the lower temperatures in all three variants at ph 1.2 ( see figure 2 ) .
the temperature at which the tertiary structure is lost differs significantly between the three variants , showing that the i56 t and i59 t single - point mutations have a large and differential effect on the stability of the native state : 48.5  0.3 c for wt , 35.4  0.5 c for i59 t , and 22.6  0.3 c for i56 t lysozyme ( see figure 1c , and table s1 in the supporting information ) .
the proteins are stable for days at this ph at low temperature ( < 35 c ) , and for several hours at the higher temperatures , enabling detailed nmr analysis to be carried out on the chemically intact protein , and allowing the equilibrium unfolding of all three variants to be compared in detail without the addition of denaturant .
hsqc spectra of i59 t at ph 1.2 recorded at different temperatures . at 20 c ,
the protein is fully folded ( a ) , whereas at 50 c , the protein is fully unfolded ( c ) . near the midpoint of the unfolding transition , at 35 c , the protein shows resonances typical of both folded and unfolded states ( b ) .
the region boxed in panel ( c ) indicates the peaks corresponding to glycine residues in the unfolded state .
incubation of 1 mm solutions of wt lysozyme at temperatures near the midpoint of unfolding ( 47.5 c ) and with stirring leads to the essentially complete conversion of the soluble forms of the protein into amyloid fibrils within 45 days .
in order to accelerate the aggregation process and avoid prolonged incubation times at elevated temperatures and low ph , which can cause hydrolysis of the protein , aliquots ( 2% v / v ) of preformed fibril seeds were added to the solutions , resulting in more than 95% of the monomeric protein becoming incorporated in fibrils within one day of incubation .
the fibrils obtained contain non - hydrolyzed , full length monomeric protein as confirmed by sds - page electrophoresis and have all the morphologic characteristics of amyloid structures : transmission electron microscopy images show that they are long and unbranched , they bind thioflavin - t , an amyloid - specific dye , and they present a high degree of -sheet structure , as indicated by cd ( see figure s1 in the supporting information ) . the observation that amyloidogenic precursors are populated near the temperatures where tertiary structure is lost prompted us to study the nature and relative abundance of such species . monitoring
the unfolding transition of all three lysozyme variants by means of circular dichroism spectroscopy ( cd ) at different wavelengths reveals that the tertiary structure ( near - uv region , 270 nm ) is disrupted more readily , i.e. at lower temperatures , than the secondary structure ( far - uv region , 222 nm ) ( see figure 3a and table s1 , figure s3 , and figure s4a in the supporting information ) . the difference in melting temperatures detected by near- and far - uv cd increases as the stability of the native state decreases ; hence , the temperature difference between the two midpoints is largest for i56 t and smallest for the wt protein .
these data show that a simple two - state model fails to describe the complexity of the unfolding process .
in addition to the native and unfolded states , at least one more state exists with a secondary structure content similar to that of the native state but with a very significant loss of tertiary structure .
interestingly , the thermal unfolding of the secondary structure is essentially identical for the i56 t and i59 t variants , whereas the tertiary structure is more strongly destabilized for the i56 t than for the i59 t variant ( see figure s3 in the supporting information ) .
thermal unfolding of human lysozyme monitored using cd spectroscopy and dsc : a pseudo - two - state model .
( a ) equilibrium thermal unfolding of wt , i59 t , and i56 t lysozyme , followed by near - uv ( 270 nm , red ) and far - uv ( 222 nm , black ) cd at ph 1.2 .
thermal melts were fitted to a two - state model ( fits in solid lines ) . for more details , see experimental section and the supporting information .
( b ) evolution of excess heat capacity with temperature shown for wt ( black ) , i59 t ( red ) , and i56 t ( blue ) at ph 1.2 .
solid lines represent the experimental data , while dashed lines show curves simulated using the thermodynamic parameters obtained from the fitting the cd data to a three - state model ( see the supporting information ) .
the reversibility of the thermal denaturation monitored by dsc is shown in figure s2 of the supporting information .
( c ) proposed model : scheme of the native state n and the denatured ensemble 
( dashed green lines ) are derived by fitting the near - uv cd first - order transition to a pseudo - two - state model .
a three - state model defining native , intermediate , and unfolded states accounts essentially perfectly for the far- and near - uv cd data,(30 ) and allows the apparent population of each state to be determined for all three variants ( see figure s4a , b in the supporting information ) . in addition , knowledge of the populations of the three states enabled us to deconvolute the far- and near - uv cd data to extract separate spectra of the native , intermediate , and unfolded states ( see figure s4c in the supporting information ) . as expected , the intermediate lacks a well - defined tertiary structure , as indicated by the loss of the native near - uv cd signal , yet preserves an essentially native - like secondary structure signature in the far - uv cd region , indicative of molten globular character ( see figure s4c in the supporting information ) .
the population of the intermediate state and the temperature window where it is significantly populated both increase as the stability of the native state decreases . in order to investigate further the nature of the unfolding transition
dsc allows the identification of merged , multiple transitions that may not be well - resolved by spectroscopic techniques , provided that these transitions induce a detectable change in both heat capacity and enthalpy.(33 ) in the present case , the unique heat absorption peak observed by dsc is characteristic of a two - state transition , where the temperature midpoint of the unfolding transition is identical within error to that seen by near - uv cd ( see figure 3a , b , experimental section , and supporting information table s1 ) ; that is , the observed transition corresponds to the loss of tertiary structure . furthermore , simulated dsc thermograms , based on the thermodynamic parameters obtained from the three - state fit to the far- and near - uv cd data , show profiles that clearly differ from the experimental ones ( compare plain and dashed lines in figure 3b ) .
we therefore conclude that only one conformational transition is detected by dsc that , according to the three - state model , corresponds to the cooperative conversion of the native state into the intermediate state .
by contrast , the unfolding of the intermediate does not follow a cooperative first - order transition , but is instead best described as a continuous ( second - order ) 
therefore , the dsc data question the validity of a three - state model invoking two first - order transitions .
the three - state model defined by the cd data is still useful in an operational sense but needs to be improved with the intermediate and unfolded states replaced by an ensemble of denatured states separated by marginal energy barriers : in this improved model the 
denatured ensemble consists of a collection of non - native states possessing different degrees of unfolding .
the composition of this ensemble changes with temperature , being on average more molten globule - like at low temperatures , and more unfolded as the temperature is increased ( see figure 3c ) . to obtain populations for the native state and the  denatured ensemble ,
we have used the near - uv cd unfolding curve only and fitted it to this pseudo - two - state model ( dashed lines in figures 3 and 4 ) .
the model that we propose here is , therefore , a significant improvement on the three - state model used until now for human lysozyme .
this situation is reminiscent of the thermal unfolding of the molten globule of the structurally homologous protein -lactalbumin whose unfolding does not give rise to any noticeable heat absorption peak and for which a similar  denatured ensemble has been proposed .
black dots connected by lines represent the temperature dependence of individual hn hsqc cross - peak volumes of the native state ( left ) , normalized relative to the volume of each peak at 5 c , and the denatured ensemble ( right ) , normalized relative to the volume of each peak at 55 c ( see comment for wt in the experimental section ) .
all non - overlapping nmr cross - peaks are shown for the decay of the native state ( left ) , whereas the increases in cross - peak volumes of glycine residues only are shown for the denatured ensemble ( right ) .
glycine residues have been selected because they have all been assigned in the denatured state of all three variants ( see figure 5b ) .
increases in cross - peak volumes for 46 assigned residues of the denatured ensemble of i56 t are shown in figure s5 in the supporting information and show similar results , therefore justifying the fact that only glycines are considered here .
the decay of the integrated signal of selected native side chain methyl groups is shown in light blue circles in the left column ( see raw data in figure s6 in the supporting information ) .
full lines represent populations derived by simultaneously fitting far- and near - uv cd data to the three - state model : the populations of the native , intermediate , and unfolded states are indicated by the blue , red , and green lines , respectively ( the population of the native state of the i56 t variant is normalized to its value at 5 c ) .
dashed lines show the results obtained by fitting the near - uv cd data to the pseudo - two - state model : the populations of the native state and the 
denatured ensemble are shown in blue and green , respectively . to characterize the unfolding process at the level of individual residues , we recorded hn backbone hsqc spectra of all three variants over a wide range of temperatures , acquiring spectra at every 2.5 c between 5 and 55 c . at low temperature ,
the hn hsqc spectrum of each variant is characteristic of a folded protein with a large number of well resolved and dispersed cross - peaks ( see figure 2a ) . by studying the decrease in volume of non - overlapping backbone cross - peaks as a function of temperature , we observe that the resonances of residues distributed throughout the protein disappear in a concerted fashion , confirming that the unfolding of the native state is highly and globally cooperative ( see figure 4 , left column ) .
nmr resonances shifted upfield of 0.5 ppm in 1d nmr spectra are characteristic of buried methyl groups in folded proteins ; we observe that these resonances decrease in intensity with increasing temperature similarly to those of backbone amides , showing at the level of individual residues that the secondary and tertiary structures are lost simultaneously in the unfolding of the native state , in agreement with our previous work ( see light blue circles on the left column of figure 4 , experimental section , and raw data in figure s6 in the supporting information).(30 ) the relative decreases in the volume of the nmr signals corresponding to the native state agree very well with the population of this state derived from near - uv cd , showing that both techniques probe the disappearance of the highly cooperative native state ( see figure 4 , left column ) . at higher temperatures , hn hsqc cross - peaks collapse into the central region of the h spectrum with little chemical shift dispersion as expected for a denatured protein ( see figure 2c ) .
two sets of nmr signals , corresponding to folded and denatured species , are observed at temperatures near the denaturation midpoints identified by cd ( see figure 2b ) .
however , only a subset of the resonances from the denatured state is present in the nmr spectra recorded under these conditions , but additional resonances appear gradually as the temperature is increased further ( see figures 4 , right column , and figure 5 ) .
this behavior is exemplified by i56 t ( see figure 5a ) : at 30 c we observe only 26 backbone amide cross - peaks out of 130 residues . as discussed further below
, this temperature dependence is similar to that observed for the low ph molten globular state of -lactalbumin , where cross - peaks appear as the temperature is increased.(34 ) unfolding of the denatured ensemble of human lysozyme monitored by nmr spectroscopy .
( a ) hn hsqc spectra of the i56 t mutant recorded at different temperatures show the gradual appearance of the peaks from different residues .
the lowest contour level was adjusted so as not to display the weak signals of the native state at 30 c , which has a population of ca .
. the lowest contour level , the number of contours as well as the increments between them are constant across all temperatures . changing the lowest contour level displayed does not affect the results qualitatively and only slightly affects the temperature at which a peak first appears ; moreover , it affects all residues in a similar way , enabling robust comparisons to be made .
( b ) region of the hn hsqc spectra ( boxed in panel a ) containing all resonances from glycine residues in the denatured ensemble .
resonances from the locally folded and unfolded residues are shown in black and red , respectively .
the temperatures of disappearance of resonances from the native state , as well as appearance of resonances from the denatured ensemble , are summarized for all three variants in figure s7 in the supporting information .
the build - up of cross - peak intensity from the denatured state is markedly non - concerted , especially in the case of the i56 t and i59 t variants ( see figure 4 , right column ) .
moreover , the range of temperatures over which nmr cross - peaks from the denatured state appear is closely similar to the range defined by the population profiles of denatured conformations without native tertiary interactions ( dashed line ) and without secondary structure ( full line ) determined by cd spectroscopy .
this result indicates that nmr is able to distinguish between classes of residues with different unfolding behaviors , ranging from those that are fully unfolded at low temperatures where native tertiary structure is initially lost , to those that become fully unfolded only upon complete denaturation of the secondary structure elements .
residue specific assignments of the denatured state provide structural insights into the locally cooperative unfolding process ( see experimental section ) . to make direct comparisons between the unfolding behavior of the different variants , we focus on a well - resolved region of the spectrum that has been completely assigned in both the folded and denatured states for all three variants .
human lysozyme has eleven glycine residues and , as they are distributed throughout the native fold , these residues act as very valuable probes of the structural properties of different regions of sequence ( see figure 1 ) .
glycine resonances cluster in an isolated and well - dispersed region of the hn hsqc spectrum of non - native lysozyme ( see squared region in figure 5a ) thus enabling us to determine accurately the relative increases in cross - peak volumes of the denatured ensemble as the temperature is increased .
as shown in figure 4 , the relative volumes of the glycine cross - peaks of non - native lysozyme vary significantly over a wide range of temperatures , showing that the unfolding of non - native states follows a non - cooperative process .
calculation of the expected water exchange rates of amide protons of the unfolded form of human lysozyme ( http://hx2.med.upenn.edu/ ) reveals that the exchange process is very slow and that the magnitudes of the rates do not correlate with the temperature of the appearance of the corresponding peaks in the nmr spectrum , therefore ruling out the possibility that line broadening is due to exchange with water ( see figure s8 in the supporting information).(35 ) expanded plots of the glycine region of the hsqc spectra show the gradual appearance of resonances from the non - native states of the protein as the temperature is increased ( see figure 5b ) .
importantly , the temperature range over which glycine cross - peaks of non - native states appear is observed to differ between wt , i56 t , and i59 t ; the range increases as the stability of the native state decreases .
as can be seen from figure 5b , this range is around 5 c for wt , 10 c for i59 t , and 27.5 c for i56 t , directly reflecting the increasing degree of global non - cooperativity in the phase transition that accompanies the thermal unfolding of the three variants of lysozyme ( see figure 4 ) .
the partial assignment obtained for i56 t enables us to identify which parts of lysozyme unfold first .
figure 6 maps the differential and gradual appearance of resonances of locally unfolded states onto the native protein structure .
the results clearly show that the -hairpin is the first structural element to unfold , followed by the rest of the -domain , helix c , both 310 helices 1 and 2 , and some residues at the chain termini . the rest of the protein , which comprises helices a , b , d and 310 helix 3 , unfolds at higher temperatures , although residues located in loops between secondary structure elements tend to unfold at lower temperatures than the elements themselves .
for all three variants , the glycine residues that unfold at the lowest temperatures are g48 , g55 , g72 , g68 , and g105 , all of which are located in the -domain or in 310 helix 2 ( see figure 5b ) the results suggest that the unfolding process of the denatured ensemble is locally cooperative and structurally similar for all three variants , with the -domain , helix c , both 310 helices 1 and 2 , and some residues of the n- and c - termini unfolding first . however , the temperature window over which unfolding takes place , which indicates the degree of global non - cooperativity of the unfolding of the denatured ensemble , is modulated by the single - point mutations of the variants .
( a ) structural location of individual residues assigned in the hn hsqc spectrum of denatured lysozyme .
the color code indicates the temperature of the first appearance of the corresponding cross - peaks in the hn hsqc spectra of the denatured state of i56 t , according to figure 5 .
the temperature evolution of chemical shifts , as provided by the temperature series of hsqcs , provides additional information on the denatured ensemble ( see figure 7a ) .
the temperature coefficients of amide protons , defined as the slope of a linear fit of h amide chemical shift versus temperature , vary across glycine residues but do not correlate with the temperature at which they appear nor with the protein sequence ( see figure 7b ) .
values more positive than 4.5 ppb / k have been attributed to the existence of intramolecular hydrogen bonds .
this is the case for three glycines ( g16 , g37 , and g105 ) , located in the -domain of the protein , whereas the remaining glycines , located in both the - and -domains , have values that are more negative than 4.5 ppb / k .
these results suggest that part of the -domain could possess some residual structure at high temperature
. however , the interpretation of amide temperature coefficients is not straightforward , especially in the case of chemically exchanging systems like lysozyme at low ph.(39 ) temperature evolution of h and n chemical shifts of the denatured ensemble . ( a ) the glycine - rich region of the hsqc spectrum of the denatured ensembles of human lysozyme variants . in each case the top spectrum , in red , corresponds to a temperature of 50 c ; spectra were recorded every 2.5 c and are all shown with the same contour levels as in figure 5 .
( b ) temperature coefficients of h amide chemical shifts of glycine residues of the denatured ensemble .
the temperature coefficients are determined as the slope of a linear fit of chemical shifts versus temperature .
results are shown for the wt ( black ) , the i59 t ( red ) , and the i56 t ( green ) variants .
( c ) plot of residuals of h and n chemical shifts , after subtraction of the best first - order fit .
the p - value shown in the inset of each graph relates to the comparison of fits to first- and second - order polynomial functions .
only residues for which one chemical shift ( h or n ) has p < 0.05 are shown .
the temperature evolution of glycine h and n chemical shifts were also analyzed in terms of their deviation from linearity , which reports on possible conformational heterogeneity .
for all three variants , we used a p - value analysis to compare the fits to first- and second - order polynomial functions : we observe that nearly all glycine residues show a linear evolution of their h and n chemical shifts with temperature , i.e. , do not undergo any important structural changes in the denatured ensemble after becoming visible in the hsqc spectrum ( see residuals in figure 7c ) .
the exceptions are g48 in the i56 t and i59 t variants , and g72 in the i56 t variant , where both residues are located in the -domain of the protein : the p - values for g72 in i56 t , and g48 in i59 t , are significantly larger than those for g48 in i56 t , showing that the extent of conformational exchange is most significant for g48 in i56 t . moreover ,
deviation from linearity for g48 is larger in i56 t than in i59 t , and non - existent in the wt protein , suggesting that the extent of conformational heterogeneity experienced by g48 within the denatured state correlates with the decrease in native state stability ( a similar conclusion holds for g72 for which deviation from linearity is only observed for i56 t ) .
we should note here that the deviation from linearity is minimal and only becomes observable when performing a statistical test : therefore , although the conformational heterogeneity probed by this analysis is statistically significant , it is of rather small magnitude .
we attribute this conformational heterogeneity to slight differences in chemical shift , and therefore structure , of g48 and g72 between the denatured states at lower and higher temperature and is most pronounced for g48 in the i56 t variant , as this residue unfolds at lower temperatures .
overall , this analysis shows that , once resonances corresponding to the denatured ensemble appear in the nmr spectrum , their chemical shift temperature evolution is not indicative of major conformational rearrangements .
in this report , we have provided a detailed characterization of the energy landscape of monomeric human lysozyme at low ph with the aim of providing key structural insights into the nature of the amyloidogenic species populated under equilibrium conditions .
we present here a comparative study of the complete thermal unfolding transition of the three variants i56 t , i59 t , and wt under relatively modest temperatures and in the absence of denaturants .
using low ph conditions has allowed us to manipulate the folding energy landscape by greatly destabilizing the native state and , for the first time for human lysozyme , to observe directly non - native states that are populated in amyloidogenic conditions .
the present investigation is therefore distinct from previous studies based on h / d exchange techniques that indirectly detect partially folded species that are transiently populated as a result of fluctuations from the native state . using a combination of dsc , cd , and nmr spectroscopy
, we show that human lysozyme unfolds in a relatively complex process that involves a cooperative transition from the native structure into a denatured ensemble that unfolds progressively with increasing temperature ; the cooperative and first - order transition detected by near - uv cd and dsc corresponds essentially exclusively to the loss of tertiary structure . in the denatured ensemble ,
tertiary contacts are disrupted and secondary structure elements unfold independently , but in a locally cooperative manner . at lower temperatures , the properties of the denatured ensemble
are characteristic of a molten globule : it does not give rise to detectable nmr resonances for most of its residues , lacks tertiary contacts as shown by near - uv cd , and retains most of the secondary structure of the native state as seen by far - uv cd .
in addition , we observe that the nmr resonances characteristic of fully unfolded conformations do not appear simultaneously at a given temperature but rather become evident gradually ( see figures 4 and 5 ) .
these observations are consistent with a model in which unfolding proceeds through a continuum of states that are relatively close in energy and separated by marginal barriers , such that they interconvert relatively rapidly ( see figure 3c ) . in agreement with this model ,
the increase in cross - peak volumes of denatured resonances takes place at higher temperatures than the near - uv cd transition ( green dashed lines in figure 4 ) .
this behavior is consistent with the observation that no heat absorption peak is detected by dsc ; hence , unfolding of the molten globule is not associated with any 
our observations are in agreement with previous reports on the thermal unfolding of hen egg white lysozyme(43 ) as well as with numerous previous calorimetric and nmr studies of -lactalbumins , proteins that are homologous to lysozymes ; however , whereas -lactalbumin populates essentially 100% of a molten - globular species at 20 c and low ph , the energy landscape of human lysozyme at low ph is significantly more complex , as its native state is more stable and is therefore highly populated at low temperatures .
therefore , -lactalbumins seem to behave at low ph as extremely destabilized mutants of human lysozyme , a phenomenon that can be attributed to the loss of a stabilizing ca ion in the former group of proteins at low ph.(44 ) a major advance in our study has been to record nmr spectra describing the whole of the unfolding transition , yielding , in particular , a structural description of the states populated in the course of the locally cooperative unfolding of the denatured ensemble . residue - specific information from nmr spectroscopy has enabled us to map the unfolding process of the denatured ensemble onto the native structure , thereby revealing which parts of human lysozyme unfold first and are therefore of lower conformational stability . within the denatured state , each secondary structure element unfolds individually in a continuous transition , resulting in the formation of a wide ensemble of partially folded conformations .
this unfolding process is structurally similar for all three variants studied here , with the chain termini , the -domain ,  helix c , and 310 helix 2 from the -domain unfolding at temperatures lower than the rest of the protein ( comprising  helices a , b , and d and 310 helix 3 ) . within the -domain
, the i56 t variant shows that the -hairpin is the first structural element to unfold .
the regions of the protein that unfold first in the denatured ensemble are the ones that are transiently unfolded under native conditions at higher ph ( ph 5.0 , 37 c ) , as detected indirectly by h / d exchange techniques : both observations are therefore in agreement and suggest that these regions are of intrinsically lower stability than the remainder of the protein structure .
the temperature evolution of amide chemical shifts of the denatured ensemble suggests the presence of some intramolecular hydrogen bonding , and therefore a modest degree of residual structure in the -domain .
once structural elements of the denatured ensemble unfold , however , there is little conformational rearrangement within this locally unfolded conformation , as indicated by the fact that the chemical shift temperature evolution of most glycine residues is essentially linear .
interestingly , both the -domain and helix c have been found to form the core of human lysozyme fibrils prepared at ph 2.0.(31 ) as we demonstrate here that the -domain and helix c have marginal individual stability , being the first regions of the protein to unfold once the tertiary structure is disrupted , this unfolding process is likely to be at the origin of the formation of amyloid fibrils at low ph . the single - point mutations i59 t and i56 t not only affect the stability of the native state of the protein , as assessed by near - uv cd , therefore influencing the temperature at which the denatured ensemble starts to unfold , but also determine the extent of non - cooperativity in the unfolding of the denatured ensemble , as measured by the temperature window over which structural elements of the denatured ensemble progressively unfold .
the more destabilized the native state , the lower the temperature at which the denatured ensemble starts to be populated , and the more pronounced its non - cooperative unfolding .
interestingly , both the native state destabilization and the extent to which cooperativity is lost within the denatured ensemble correlate with the amyloidogenicity of the variant in other solution conditions ( ph 5.0 ) , as well as with the expression levels of lysozyme in pichia pastoris.(29 ) although the results of this study do not directly report on the structural properties of the amyloidogenic state of lysozyme under physiological conditions , they strongly indicate that the propensity of human lysozyme to aggregate is closely related to the shape of the energy landscape of the monomeric protein,(25 ) that has been probed in the present study .
the i56 t and i59 t variants and wt human lysozyme were expressed in aspergillus niger ( i56 t ) or pichia pastoris ( i59 t , wt ) and purified as described previously .
all experiments were carried out with the protein samples in 50 mm phosphate buffer adjusted to ph 1.2 with hcl .
circular dichroism experiments were carried out using a jasco j-810 spectropolarimeter equipped with a peltier holder .
spectra were recorded using cells of 0.1 and 1 cm path length for the far- and near - uv , respectively , with a protein concentration of ca .
heat - induced unfolding transitions were monitored from 5 to 90 c at a rate of 0.5 c / min or 1 c / min ( the data obtained were independent of the heating rate ) .
data points were acquired every 0.5 c with a response time of 2 s and a bandwidth of 1 nm . fitting to a two - state model was carried out after normalization of the unfolding transition ; the three mutants were fitted simultaneously , with common baselines for the native and the denatured states ( see equations in the supporting information ) .
fitting to a three - state model was carried out for each mutant separately , after normalization of both near- and far - uv transitions ( see equations in the supporting information).(33 ) the cp of the first transition was taken as 1.6 kcal / molk for i59 t and wt , and 1.3 kcal / mol.k for i56 t , according to previous studies.(46 ) because the theoretical cp calculated for the unfolding of the native state to the completely unfolded state is 1.6 kcal / molk , and therefore very close to the measured cp of the first transition , we assume that the cp for the second transition is much smaller than this value and arbitrarily set it at 0.2 kcal / molk for all three variants.(47 ) distributing the cp equally over the two transitions ( for example with cp = 0.8 kcal / mol.k for both transitions ) had very little effect on the fitting , validating our approach to the data analysis .
heat capacities were measured as a function of temperature using a vp - dsc microcalorimeter ( microcal llc , northampton , ma ) , with protein concentrations of ca .
40 mol / l and a scanning rate of 1 c / min , from 10 to 70 c .
the baseline obtained by heating a buffer solution was subtracted from the protein unfolding thermogram , and the excess heat capacity ( taking the native state as reference ) was obtained by subtracting a progressive baseline between the native and unfolded states .
the reversibility of the transition was assessed by overlaying the first and second thermal unfolding curves , the latter being obtained by reheating the sample after cooling inside the calorimetric cell .
the reversibility was greater than 90% for all variants ( see figure s1 in the supporting information ) .
vant hoff to calorimetric enthalpy ratios were close to unity for the i59 t and wt variants ( it was not possible to estimate it for i56 t because of the lack of native state baseline at low temperatures ) , and data could only be fitted to a two - state model ( see equations in the supporting information ) .
n ammonium sulfate and c methanol were used to label the protein produced in pichia pastoris with n and c isotopes , respectively .
fast hsqc schemes for hn 2d spectra were used for monitoring the unfolding process , and spectra were acquired from 5 to 55 c , every 2.5 c.(48 ) the backbone resonances of the folded state of the i59 t variant were assigned at 25 c using hnca , hncacb , and cbca(co)nh experiments.(49 ) conventional nmr characterization of the non - native ensemble of lysozyme , that requires several days of measurements at high temperature and low ph ( ph 1.2 ) , has previously been limited by the degradation of the sample . using a procedure based on a 3d version of the n zz - exchange experiments ,
however , enabled data collection to be achieved in 12 days : the denatured state of the i59 t variant was partially assigned ( ca .
82% ) , by observation of the transfer of magnetization at 35 c between cross - peaks arising from native and non - native states .
some resonances of the i56 t variant could be correlated with those of the i59 t variant by overlaying the spectra of the denatured states , resulting in the assignment of ca .
peak volumes of hsqc cross - peaks were determined assuming gaussian lineshapes , and overlapping peaks were discarded . in figure 4 , the traces of the increase in peak volumes for the wt are slightly to the left of the green dashed line , because normalization was carried out relative to a temperature of 55 c , where the far - uv for wt has achieved only 70% of its transition .
 the partial unfolding of human lysozyme underlies its conversion from the soluble state into amyloid fibrils observed in a fatal hereditary form of systemic amyloidosis . to understand the molecular origins of the disease , it is critical to characterize the structural and physicochemical properties of the amyloidogenic states of the protein . here 
 we provide a high - resolution view of the unfolding process at low ph for three different lysozyme variants , the wild - type protein and the mutants i56 t and i59 t , which show variable stabilities and propensities to aggregate in vitro . using a range of biophysical techniques that includes differential scanning calorimetry and nuclear magnetic resonance spectroscopy 
 , we demonstrate that thermal unfolding under amyloidogenic solution conditions involves a cooperative loss of native tertiary structure , followed by progressive unfolding of a compact , molten globule - like denatured state ensemble as the temperature is increased . the width of the temperature window over which the denatured ensemble progressively unfolds correlates with the relative amyloidogenicity and stability of these variants , and the region of lysozyme that unfolds first maps to that which forms the core of the amyloid fibrils formed under similar conditions . together , these results present a coherent picture at atomic resolution of the initial events underlying amyloid formation by a globular protein . 