Protein Science
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Protein Science (2004), 13:2628-2638. Published by Cold Spring Harbor Laboratory Press. Copyright © 2004 The Protein Society
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Sandberg, A.
Right arrow Articles by Karlsson, B. G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Sandberg, A.
Right arrow Articles by Karlsson, B. G.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Reddit   Add to Technorati  
What's this?

Apo-azurin folds via an intermediate that resembles the molten-globule

Anders Sandberg1, Johan Leckner3 and B. Göran Karlsson2,3

1 Department of Chemistry and
2 Swedish NMR Centre, Göteborg University, SE 405 30 Göteborg, Sweden
3 Department of Chemistry and Bioscience, Chalmers University of Technology, SE 405 30 Göteborg, Sweden

Reprint requests to: B. Göran Karlsson, Swedish NMR Centre, Göteborg University, Box 465, SE 405 30 Göteborg, Sweden; e-mail: goran{at}nmr.se; fax: +46-31-773-3880.

(RECEIVED May 7, 2004; FINAL REVISION July 1, 2004; ACCEPTED July 2, 2004)


    Abstract
 TOP
 Abstract
 Introduction
 Results
 Discussion
 Materials and methods
 References
 
The folding of Pseudomonas aeruginosa apo-azurin was investigated with the intent of identifying putative intermediates. Two apo-mutants were constructed by replacing the main metal-binding ligand C112 with a serine (C112S) and an alanine (C112A). The guanidinium-induced unfolding free energies ({Delta}GU–NH2O) of the C112S and C112A mutants were measured to 36.8 ± 1 kJ mole–1 and 26.1 ± 1 kJ mole–1, respectively, and the m-value of the transition to 23.5 ± 0.7 kJ mole–1 M–1. The difference in folding free energy ({Delta}{Delta}GU–NH2O) is largely attributed to the intramolecular hydrogen bonding properties of the serine O{gamma} in the C112S mutant, which is lacking in the C112A structure. Furthermore, only the unfolding rates differ between the two mutants, thus pointing to the energy of the native state as the source of the observed {Delta} {Delta}GU–NH2O. This also indicates that the formation of the hydrogen bonds present in C112S but absent in C112A is a late event in the folding of the apo-protein, thus suggesting that formation of the metal-binding site occurs after the rate-limiting formation of the transition state. In both mutants we also noted a burst-phase intermediate. Because this intermediate was capable of binding 1-anilinonaphtalene-8-sulfonate (ANS), as were an acid-induced species at pH 2.6, we ascribe it molten globule-like status. However, despite the presence of an intermediate, the folding of apo-azurin C112S is well approximated by a two-state kinetic mechanism.

Keywords: protein folding; intermediate; {beta}-sandwich; Greek key; azurin; cupredoxin

Abbreviations: ANS, 1-anilinonaphtalene-8-sulfonate • {beta}TS, Tanford {beta}-value for the transition state • {beta}I, Tanford {beta}-value for the intermediate state • CD, circular dichroism • [D], denaturant concentration • [D]50%, mid-point of denaturation • GdmCl, guanidinium chloride • I, intermediate state • N, native state • TS, transition state • U, unfolded state.

Article and publication are at http://www.proteinscience.org/cgi/doi/10.1110/ps.04848204.


    Introduction
 TOP
 Abstract
 Introduction
 Results
 Discussion
 Materials and methods
 References
 
Anfinsen demonstrated over 40 years ago that the spontaneous acquisition of the three-dimensional structure of proteins is an inherent property of the polypeptide chain (Haber and Anfinsen 1962; Anfinsen 1973). This implies that protein folding is strictly under thermodynamic control, resulting in an end state representing the global energetic minimum. However, as Levinthal (1968) succinctly pointed out, because of the vast number of conformations available to a protein, the search for the native state must be a nonrandom process involving kinetically accessible pathways that directs the folding by restricting conformational space. Indeed, it would typically take 1042 years for a 100-residue protein to randomly search all possible conformations. A crucial feature of the observation made by Levinthal (1968) is that the native state need not be the most thermodynamically stable structure, but only the structure reached through folding via the fastest pathway. Since then, much effort has been spent on the identification and characterization of such pathways, as well as on the factors that govern the thermodynamic stability of proteins.

Recently, there has been considerable focus on proteins dominated by long-range interactions, such as structures with primarily {beta}-secondary structural elements. The all-{beta} group of proteins is dominated by {beta}-barrels and {beta}-sandwich structures which usually have Greek-key topology (Chothia et al. 1997). Its prevalence is an indication that the folding of the polypeptide into a coherent all-{beta} structural entity is limited by geometrical and/or kinetic-mechanistic factors. Clarke and coworkers (Clarke et al. 1999; Fowler and Clarke 2001) undertook the first extensive experimental folding analysis of the Greek-key immunoglobulin-like domains, where they propose a common nucleation-condensation mechanism for this group of proteins. They concluded that the folding nucleus is centered on a common structural core made up of nonlocal residues that are also important in defining the structure of these proteins. Intriguingly, nearly all (94%) {beta}-sandwich-like proteins do contain an invariant supersecondary substructure defining the core which is even more conserved than the Greek-key topology (Kister et al. 2002). Hence, analysis of the folding behavior of {beta}-sandwich proteins other than the Ig-like domains may reveal whether there also is a "super pathway" in the folding of these proteins. Herein lies our main interest in the folding of the cupredoxins, which is a structurally homogeneous but sequentially inhomogeneous family of proteins.

Our model protein for analysis of the folding of Greek-key proteins is the cupredoxin azurin from Pseudomonas aeruginosa. This is a fairly small protein of 14 kDa with 128 residues wrapped into an eight-stranded antiparallel {beta}-sandwich structure (Fig. 1Go). The native protein binds a copper ion (but can also accommodate other metals) with a 5[N2S(OS)] trigonal bipyramidal coordination at the top of the tightly wound structure (Gray et al. 2000). The high hydrophobicity of the azurin core is indicated by the fact that the only tryptophan (W48) has a structured emission in the native protein resembling that of indole in cyclohexane. Its emission maximum at 308 nm is probably the least red-shifted tryptophan fluorescence found in proteins (Eftink 1991). Hence, the integrity of the native tertiary structure can uniquely be monitored by the tryptophan fluorescence at this wavelength with great sensitivity. In addition, there is a disulfide bond that stabilizes one of the two sheets, but there are no cis prolines in this protein.



View larger version (36K):
[in this window]
[in a new window]
 
Figure 1. The structure of azurin with a schematic presentation of its Greek-key topological fold. The metal ion is coordinated in a plane by C112, H117, and H46 (also shown in the structure). In addition, there are two axial ligands of weaker character, G45 and M121, that complete the 5[N2S(OS)] trigonal bipyramidal coordination geometry. * indicates the position of the W48 fluorescent probe in the topology. This residue is also shown in the structure. A pattern found in ancient Greek pottery depicted at the bottom gave this topological motif its name (Richardson 1977). The structure was drawn using Molscript (Kraulis 1991).

 
Previous studies on proteins from the cupredoxin family have pointed out the existence of kinetic intermediates in apo-pseudoazurin (Capaldi et al. 1999; Reader et al. 2001) and apo-plastocyanin (Koide et al. 1993). Similarly, other {beta}-sandwich proteins, for example, the fibronectin type III module and other Ig-like domains (Clarke et al. 1999), also display at least three-state kinetics. In contrast, there are several reports stating that apo-azurin folds via a two-state mechanism (see Pozdnyakova and Wittung-Stafshede 2003 and references therein). We find it rather unlikely for a protein of such complex topology dominated by long-range interactions to fold by a two-state mechanism, and thus decided to investigate this further. In the present study we analyzed the folding properties of two apo forms of this protein. In order to eliminate distortion of the data by metal contamination and/or cysteine oxidation of the apo form, we constructed mutants in which the metal ligating C112 residue was replaced with serine and alanine. These azurin mutants were then subjected to equilibrium and kinetic folding studies using guanidinium chloride (GdmCl) as denaturant.

As expected for a protein of such size and complex topology, we found the existence of multiphasic kinetics and molten globule-like burst phases in these two apo-mutants. Thus, we obtained evidence suggesting that apo-azurin does not fold by a two-state mechanism with observed curvature being attributed to movement of the transition state (TS) alone (Pozdnyakova and Wittung-Stafshede 2003), although a two-state approximation is a fairly reasonable approximation under the physicochemical conditions studies here. We also conclude, based on {phi}-value analysis, that the formation of the metal site is a late event in the folding pathway.


    Results
 TOP
 Abstract
 Introduction
 Results
 Discussion
 Materials and methods
 References
 
Equilibrium analysis of C112S and C112A
Equilibrium folding curves for the two apo mutants is shown in Figure 2Go, and the parameters obtained by least-squares fitting of a two-state model to the data are presented in Table 1Go. Table 1Go also contains parameters from an apo-wild-type (WT) control experiment (data not shown). The transitions are fully reversible and conform to two-state behavior irrespective of the probe monitored, indicating an apparent global conformational change in GdmCl. The mU–N values are also similar for both mutants, as well as for the apo-WT protein, with an average of 23.5 ± 0.7 kJ mole–1 M–1 for the seven independent measurements. The C112S mutant is the more stable of the two mutants, with a midpoint of denaturation ([D]50%) centered at 1.55 M and an average free energy of folding of 36.8 ± 1 kJ mole–1 for the three different measurements. This indicates a slight destabilization of C112S with respect to the apo-WT protein by 3.0 ± 2 kJ mole–1. Replacing the serine with an alanine, that is, the C112A mutant, apparently destabilizes the C112S protein by 9.9 ± 1 kJ mole–1 by shifting the midpoint of denaturation to 1.13 M.



View larger version (20K):
[in this window]
[in a new window]
 
Figure 2. Equilibrium denaturation curves of C112S (A) and C112A (B). Fluorescence emission was collected at 308 nm using an excitation wavelength of 295 nm (triangles). Far-UV CD was collected at 220 nm (circles) and near-UV CD at 281 nm (squares). The protein concentration was 5 µM for the fluorescence, 2.5 µM for the far-UV CD, and 0.1 mM in near-UV CD measurements. The buffer was 0.5 mM K phosphate (pH 7.0).

 

View this table:
[in this window]
[in a new window]
 
Table 1. Equilibrium folding data
 
Biphasic kinetics of C112S and C112A folding
The refolding traces of C112S are biphasic at all concentrations of GdmCl ([D]) in both circular dichroism (CD) and fluorescence experiments. This contrasts previous reports of monophasic refolding kinetics for apo-azurin (Pozdnyakova et al. 2001; Pozdnyakova and Wittung-Stafshede 2003). Figure 3Go shows a typical trace of C112S refolded in 1.40 M GdmCl as monitored by fluorescence. At [D]<<[D]50% the amplitude of the minor phase (k2) is typically less than 10% of the total amplitude (Fig. 4Go). The minor phase reaches a maximum in relative total amplitude of 20% at ~1.40 M GdmCl for both fluorescence and CD data. This minor phase is also evident in the unfolding limbs within the transition region up to ~2.0 M GdmCl. At [D]<<[D]50%, however, the unfolding traces are all monophasic.



View larger version (41K):
[in this window]
[in a new window]
 
Figure 3. Refolding of C112S in 1.4 M GdmCl from a 3.0 M GdmCl solution (A) with resultant residuals for mono-exponential (B) and a bi-exponential (C) fits. The unfolded protein (20 µM) was diluted 1:10 with refolding buffer. The buffer was 0.5 mM K phosphate (pH 7.0).

 


View larger version (14K):
[in this window]
[in a new window]
 
Figure 4. Relative amplitudes of the minor phase of the fluorescence (triangles) and far-UV CD (squares) data for C112S folding kinetics.

 
Observed nonlinearity in relative amplitudes of the minor and major phases (k1) was not due to concentration-dependent aggregation effects, as judged from refolding from 3 M to 1.4 M GdmCl with end concentrations of protein ranging from 1.5 to 30 µM (data not shown). Interestingly, the population folding via the slower pathway reaches a maximum at a GdmCl concentration that equals the kinetic [D]50% of this minor phase (cf. Figs. 4Go and 5A,BGo), and this [D]50% does not equal the [D]50% of the major phase. Similarly, in the C112A mutant, a minor phase with nonlinearity in its dependence of k2 on [D] is also present (Fig. 5CGo). The relative amplitudes of the minor and major phases are similar to that for C112S, although the maximum of the minor phase (also 20% of the total phase) is shifted to lower GdmCl concentrations (1.0 M; data not shown). Furthermore, there were no changes in the relative amplitudes of the two phases when refolding C112S from different start concentrations of denaturant (the initial conditions test; see Wallace and Matthews 2002), so the slow phase is probably not a sequential process. The latter experiment also rules out the possibility that significant residual structure is a cause of the observed biphasic kinetics.



View larger version (22K):
[in this window]
[in a new window]
 
Figure 5. Chevron plots and amplitude data for C112S and C112A folding kinetics. (A) Fluorescence data for C112S with refolding (filled triangles) and unfolding (open triangles) kinetics of the major phase (k1) complemented with ANS data (open squares). The solid line is the least-squares fit to equation 4 (see Materials and Methods), and the broken line represents the two-state approximation using equation 2. The refolding (filled diamonds) and unfolding (open diamonds) of the minor phase (k2) is also shown. (B) The far-UV CD data for C112S with symbols as in panel A. The unfolding data of the minor phase were of too poor quality to be of any value (not shown). (C) The fluorescence data for C112A with symbols as in A. (D) The normalized extrapolated start (open triangles) and end (filled triangles) values for C112S folding as monitored by fluorescence, and the start (open squares) and end (filled squares) values for the unfolding. (E) The normalized extrapolated start and end values for C112S far-UV CD data denoted as in D. (F) The extrapolated start and end values for the C112A fluorescence data with symbols as in D.

 
Multistate folding of C112S
Chevron plots (ln kobs versus [GdmCl]) of the major and minor phases are shown in Figure 5Go, A for the fluorescence data and B for the CD data. The dependence of the natural logarithm of the rate constant associated with the major phase (k1) on GdmCl concentration is clearly nonlinear for both mutants. This curvature was not caused by transient aggregation (Silow and Oliveberg 1997b) as determined by additional refolding experiments at 1, 2, 5, and 10 µM protein concentrations (data not shown). The electrolyte nature of the denaturant as causing the curvature (Monera et al. 1994; Makhatadze 1999) was ruled out by preliminary data obtained by urea denaturation (A. Sandberg, unpubl.). Nonlinearity in the activation free energies of folding has been observed in a number of proteins and is usually interpreted in terms of movement of the TS across a broad energy barrier as the physicochemical properties of the solvent change upon addition of denaturant (Matouschek and Fersht 1993; Silow and Oliveberg 1997a). Alternatively, downward curvature in the Chevron plots may also be due to a denaturant-dependent switch between two or a few distinct TSs on a sequential pathway (Sanchez and Kiefhaber 2003).

Inspection of the extrapolated start values of the fluorescence refolding traces as shown in Figure 5DGo reveals the presence of a burst-phase intermediate becoming accumulated at denaturant concentrations below 0.70 M GdmCl. This burst phase is not evident in the CD data (Fig. 5EGo), possibly because this data is generally noisier. The existence of refolding intermediates is also implicated from the additional curvature in the refolding limb (Cf) of the major phase compared to the unfolding limb (Cu) in the Chevron plots. For the fluorescence data, the refolding limb had a squared parameter of –(1.45 ± 0.2)[D]2, whereas for the unfolding limb the same value was – (0.99 ± 0.2)[D]2. This additional curvature of –0.46 ± 0.2 is not significant enough to allow for equation 4 to converge properly, however. For the CD data the additional curvature in the refolding limb was well within experimental error.

Using the rate of unfolding kinetics and the equilibrium data, it is possible to calculate the expected refolding kinetics by using the requirements of a two-state transition (Jackson and Fersht 1991a). The expected folding rate in water is determined by:


(1)

and the full expression for the expected two-state refolding rates is (Jackson and Fersht 1991a; Silow and Oliveberg 1997a):


(2)

Hence, using this expression to calculate a theoretical Chevron profile (as kobs = kf, calc + ku) one expects the kinetic refolding m-value and curvature to be equal to those obtained from the kinetic unfolding data (here mu and Cu, respectively) because unfolding usually approximates two-state. This is shown in Figure 5, A and BGo, for the major phases of the C112S mutant as dashed lines. In our opinion, a two-state mechanism with curvature-ascribed movement of the TS is a good approximation for the folding of this mutant, although it appears that intermediates may be present.

Kinetic analysis of C112A folding
In order to probe the degree of involvement of the copper site in the folding of apo-azurin, fluorescence folding kinetics were collected for the C112A mutant. The results are nearly identical to those obtained with C112S, although the C112A mutant is somewhat more unstable, as reflected in the increased unfolding rate (see Table 2Go). There is also a more pronounced curvature in the dependence of the natural logarithm of the refolding rate constant k1 on denaturant concentration when fitting equation 3 to the data: The re-folding limb had a parameter of – (2.35 ± 0.2)[D]2 defining the curvature, which is significantly more than the – (0.88 ± 0.1)[D]2 for the unfolding limb. This additional curvature of –1.46 ± 0.3 is also evident from the two-state approximation of refolding data (dashed line in Fig. 5CGo). A function describing the folding behavior with intermediates (equation 4), which ascribes all of the additional curvature to the accumulation of an intermediate, did converge properly in this case without affecting the quality of the fit (data not shown; this fit is virtually indistinguishable from the solid line in Fig. 5CGo), presumably because of the more pronounced curvature (parameters in Table 2Go). A burst phase is also evident from the extrapolated start values presented in Figure 5FGo. This intermediate starts to accumulate at approximately the same GdmCl concentration as does the C112S intermediate (~0.70 M).


View this table:
[in this window]
[in a new window]
 
Table 2. Kinetic data for the major folding phase
 
Kinetic and equilibrium analysis of C112S and C112A molten globules
1-anilinonaphtalene-8-sulfonate (ANS) is an extrinsic probe for analyzing hydrophobic sites on proteins (Stryer 1965). It has been used extensively to monitor conformational changes and partially folded states in proteins because of its tendency to bind to hydrophobic pockets of proteins, but usually not to the native or denatured states (Semisotnov et al. 1991). Protein species capable of binding ANS are usually indicative of the molten globule state (Semisotnov et al. 1991; Ptitsyn 1995). These are partly folded compact structures that preserve a protein core and native-like secondary structure, but lack the tight packing of side chains and, hence, are more flexible than the native state (Arai and Kuwajima 2000). As the burst-phase intermediate observed during azurin refolding has some features consistent with this definition, and because an ANS-binding species has been demonstrated during apo-pseudoazurin refolding (Reader et al. 2001), we decided to investigate whether this intermediate is capable of binding ANS. Indeed, transient ANS binding follows the appearance of the burst-phase intermediate closely for both mutants (Fig. 5Go), and its amplitude increases exponentially (or sigmoidally) at denaturant concentrations under approximately 0.70 M GdmCl (data not shown). Furthermore, these refolding traces could only be described satisfactorily by a triple-exponential decay (Fig. 6Go). Notably, both mutants have very similar ANS-binding kinetics. In addition, there was no detectable ANS fluorescence at [D] above {approx}0.7 M GdmCl or in any of the unfolding traces measured.



View larger version (18K):
[in this window]
[in a new window]
 
Figure 6. Kinetics collected for refolding of C112S (filled symbols) and C112A (open symbols) monitored by ANS fluorescence. (A) The kinetics of the major phase (circles) and the two minor phases (squares and diamonds). (B) The relative amplitudes associated with the major phase (circles) and minor phases (squares and diamonds).

 
Acid-induced structural changes in many proteins are consistent with the definition of the molten globule, and they also seem to resemble their kinetic counterparts (Arai and Kuwajima 2000). Azurin displays such a molten-globule at pH 2.6, depicted in Figure 7Go, with a strong increase in fluorescence quantum yield of the ANS probe at this pH (Fig. 7A,BGo insets). The fluorescence spectrum of the C112S mutant at pH 2.6 resembles a mixture of native and denatured states in having a significantly red-shifted yet structured emission. A red shift is expected upon contact of the indole side chain with hydrogen-bonding partners as the hydrophobic core relaxes around the fluorophore, but the blue-shifted emission suggests a persistent hydrophobic environment. Molten globules may indeed have significant hydrophobic interactions resistant to solvent penetration (Denisov et al. 1999). The C112A mutant, on the other hand, has an unstructured fluorescence spectrum at pH 2.6 that is almost identical to solvent-exposed tryptophan with a broad peak at 355 nm, and is thus indistinguishable from denatured protein. That these mutants differ in this respect may be due to the different hydrogen-bonding properties of the serine and alanine residues of the copper site, which is located directly above the tryptophan (Fig. 1Go). The original spectra of both mutants are regained upon increasing the pH to 7.0 again (data not shown).



View larger version (26K):
[in this window]
[in a new window]
 
Figure 7. Fluorescence spectrum of native (pH 7.0) and acid-induced (pH 2.6) molten globule-like states of azurin C112S (A) and C112A (B). The insets show fluorescence spectra of the ANS probe bound to both states.

 

    Discussion
 TOP
 Abstract
 Introduction
 Results
 Discussion
 Materials and methods
 References
 
Structure and stability
We investigated the folding properties of two mutants of P. aeruginosa azurin, C112S and C112A, that have been deprived the main metal-coordinating ligand C112. The choice of constructing these apo-mutants instead of using the wild-type apo-protein was based on the desire to dismiss the possibility of metal contamination from laboratory equipment and chemicals, as this protein has a high affinity for various metal ions. Of special interest in this respect is the C112-Cu2+ interaction, which is of near-covalent nature in this protein (Solomon et al. 1992, 1998). In addition, the C112 residue is prone to transition-metal catalyzed oxidation in the denatured state (Sandberg et al. 2002). Both of these factors are difficult to control, but are eliminated by removal of the cysteine. The C112S mutant was characterized previously and was shown not to coordinate metal ions (Hoitink 1993). High-resolution crystal structures confirm the conformational integrity of both mutants and their inability to coordinate metal ions (M. Ökvist, A. Sandberg, B.G. Karlsson, and L. Sjölin, unpubl.).

Compared to wild-type apo protein, the effect of the cysteine to serine substitution is very small, resulting in a {Delta} {Delta}GU–N of only 3.0 ± 2 kJ mole–1. A larger effect was seen when removing the hydrogen-bonding properties of the 112 residue: 9.9 ± 1 kJ mole–1 was lost upon the serine to alanine substitution. We note a significant discrepancy of 12 ± 2 kJ mole between our data for apo-WT and data previously published for the same species (Pozdnyakova et al. 2001). We believe our data to be a more correct estimate, as only 10 data points were used in defining the transition curve in the earlier publication, of which only three fell within the actual transition region. Our data are also in accord with that obtained in this study for the C112S mutant with several different optical probes, and this species should realistically be nearly isoenergetic with apo-WT.

Folding pathway
The importance of the hydrogen bonds to the Ser112 O{gamma} in the folding process can be elucidated with {phi}-value analysis (Fersht 1999). However, it is important to ascertain that the removal of these hydrogen bonds does not affect the TS. An estimate of the relative exposure of the hydrophobic surfaces of the two TSs relative to the folded state can be obtained from the Tanford {beta}-value ({beta}TS). For both mutants the {beta}TS obtained for the fluorescence data is close to 0.4, indicating that significant compactness is lacking, and suggesting a very similar nature of the TS in both mutants.

A comparison of the folding limbs in Figure 5Go clearly demonstrates that nearly all of the difference in folding stability between the mutants lies on the kinetic unfolding limb. Indeed, the {phi}f-value for the fluorescence data at 0.7 M, where the contribution of folding intermediate is negligible, is virtually zero ({phi}f = 0.04). This means that the hydrogen bonding to the O{gamma} of Ser112 is a late event which occurs after the formation of the TS in the folding pathway. Based on this, it follows that the metal site is not a stabilizing factor in early folding events, but rather forms very late in the folding process.

Folding by the minor phase
There is a significant heterogeneity in the refolding reaction of both mutants. That both fluorescence and CD report similarly on the slow-folding population indicates that it is not an artifact of fluorescence, nor is it a local folding event. Evidence of this refolding heterogeneity comes mainly from three observations. First, there is biexponential kinetics at all [D]. Second, this results in two vertically and horizontally displaced Chevron curves. Third, there is a nonlinear dependence of the relative amplitudes of the two phases. The initial conditions test (Wallace and Matthews 2002) indicates no dependence on start conditions of folding, which implies that there is probably not a significant equilibrium between different macroscopic compact states in the denaturing solutions used for the refolding studies, despite 2–3 M GdmCl being a bad solvent for polypeptide chains (Dill and Shortle 1991). Although most of these observations may be attributed to the existence of a parallel folding pathway (Wallace and Matthews 2002), the horizontally displaced Chevron profile of the slow phase relative to the fast phase cannot. Instead, it suggests the presence of an alternative species folding in parallel.

The occurrence of a slow kinetic phase (k2) with a maximum relative amplitude of ~20% somewhat displaced from the midpoint of the U {leftrightarrows} N folding transition of the faster kinetic phase (k1) is consistent with isomerization of cis proline residues being rate-limiting for the folding of a fraction of the population (Kiefhaber et al. 1992). However, that the rate of the isomerization reaction is dependent on [D], and the observation that its amplitude decreases in magnitude under increasingly stronger native conditions, indicate that the rate of conversion of cis into trans configuration is increased by the conformational strain exerted by the native or intermediate structure. Such folding-assisted proline isomerization has been observed previously, for example in the refolding kinetics of chymotrypsin inhibitor 2 (Jackson and Fersht 1991b). As the appearance of the burst-phase intermediate occurs only after the relative amplitude of the slow phase has lost most of its magnitude, the molten globule-like intermediate appears not to be a result of incorrect isomers impeding the formation of the correct tertiary interactions.

Three-state folding of the major phase
Interestingly, the burst-phase intermediate in the fluorescence implies a compact core with high hydrophobicity, because fluorescence at 308 nm is unusual in proteins as it indicates a very small Stokes shift (Eftink 1991). Even though a burst phase is not resolvable in the noisier CD data, such significant core packing also implies a satisfied hydrogen-bonding pattern in this region and, hence, significant secondary structure. This is corroborated by its ability to bind ANS, which selectively monitors the I {leftrightarrows} N transition. Burst-phase intermediates are, however, controversial in that they may be a result of experimental artifacts (Krantz et al. 2002), but the concomitant formation of a species capable of binding ANS suggests that it is a proper intermediate. It is also possible that the acid-induced molten globules identified in this study (Fig. 7Go) resemble this kinetic intermediate, as has been demonstrated for several other proteins (Arai and Kuwajima 2000), but the lack of fluorescence of this state for the C112A mutant suggests that it is more compact in the high-salt conditions of guanidinium refolding experiments. Because of the identical burst appearance of a fluorescent intermediate and the identical disappearance of the molten-globule intermediate in both of these mutants, there is no reason to suspect that the intermediate differs between the two mutants. This reasoning is corroborated by the {phi}-value analysis, which points to the copper site as being detached from the folding of the rest of the {beta}-sandwich structure. From our data it is not discernible whether the burst-phase species and the molten globule-like species are the same, but judging from their co-appearance they probably are.

The kinetic parameters from fitting indicate that there is curvature of comparable magnitude in the unfolding limbs of both mutants (Table 2Go). Whether this movement of the TS is due to continuous Hammond behavior across broad energy barriers or to a discrete change in the rate-limiting step of the reaction remains to be investigated. Nevertheless, for the fluorescence data of both mutants there is additional curvature in the refolding limbs that, in light of the presence of intermediates, may not be attributed to the movement of the TS along the reaction coordinate a priori. This warrants analysis of these Chevron profiles with rate equations that take intermediates into account (Parker et al. 1995). However, the "kink" resulting from the shift in the rate-limiting step as the intermediate accumulates is not well resolved in the C112S mutant, and a three-state analysis of the kinetic parameters is hampered by this fact. The C112A mutant is, on the other hand, less stable and kinetic data can be collected at even stronger native-like conditions (0.2 M GdmCl as opposed to 0.3 M GdmCl for C112S), which may be one of the factors contributing to the kink being better resolved in this mutant. It would thus appear from the data presented above that only the C112A mutant reports qualitatively on the three-state folding mechanism of azurin. As the Tanford {beta}-value can be regarded as a reaction coordinate by reflecting the progression of compaction during folding, an energy diagram for the U {leftrightarrows} I {leftrightarrows} N scheme can be constructed for the folding of this mutant (Fig. 8Go).



View larger version (17K):
[in this window]
[in a new window]
 
Figure 8. An energy diagram for the major species of the plausible three-state U {leftrightarrows} I {leftrightarrows} N folding mechanism of C112A plotted against the {beta}-Tanford reaction coordinate. The Tanford {beta}-values for the transition state (TS) and intermediate (I) species were taken according to equations 6 and 7, respectively, and the free energy was calculated using equation 9. When the denaturant concentration approaches 0.7 M GdmCl, the intermediate becomes destabilized and is no longer rate-limiting for the folding reaction (illustrated by the dashed line) which thus approximates a two-state reaction.

 
The data presented in this study thus strongly suggest that both apo-azurin mutants, and hence very likely also apo-WT, fold by a three-state mechanism. As azurin starts populating the intermediate at ~0.7 M GdmCl, the kinetics are slowed down accordingly, giving rise to the additional curvature in the refolding limb, although this step is not as well resolved in the C112S mutant as it is in the C112A derivative. The similarity of the kinetics for both mutants monitored by ANS-binding, however, suggests to us that extrapolation of these data into H2O corresponds to the speed limit of azurin folding (76 sec–1 and 90 sec–1 for the C112S and C112A mutants, respectively). In other words, we believe that the rate-limiting step in the folding of azurin is the conversion of the molten globule into the native state, although for the set of proteins analyzed here the speed limit of approximately 83 sec–1 (the average of both mutants) is only approached, not reached. Admittedly, treating the data with an off-pathway model (kf = KU–I x kin) does not give folding rates that are altogether unreasonable, but they are nevertheless highly unlikely for the complex topology and size of azurin (Plaxco et al. 1998). In addition, it significantly exceeds the expected two-state folding rate obtained by extrapolation. Taken together, the results presented herein are in agreement with what is to be expected for the folding of proteins exceeding ~100 residues (Jackson 1998), and considering that the unusually high meq value for azurin reflects a large core of high hydrophobicity, a rapid collapse into a molten globule-like intermediate is to be expected. Because this compact state has native-like fluorescence at 308 nm, it is probably not a pre-molten globule or a denatured state at physiological conditions, but rather a proper intermediate with native-like features and ANS-binding properties.


    Materials and methods
 TOP
 Abstract
 Introduction
 Results
 Discussion
 Materials and methods
 References
 
Chemicals
The GdmCl was of highest grade from Sigma or Apollo Scientific. The concentration of GdmCl was determined by refractive index (Nozaki 1972). All other chemicals were of reagent grade from Fluka.

Expression and purification
The two azurin mutants in this study (C112S and C112A) have been described (Sandberg et al. 2002). They were expressed and purified by a method analogous to wild-type azurin, which is described by Karlsson et al. (1989), with the exception that the cultivation medium was not supplemented with 100 µM CuSO4 and all additions of CuSO4 during the purification procedure were omitted. Expression yields of the apo-mutants were comparable to wild-type protein using the cited protocol, that is, approaching 100 mg per L culture.

Equilibrium data collection and analysis
Fluorescence equilibrium denaturation was monitored by excitation at 295 nm and emission at 308 nm on a SPEX Fluorolog {tau}2 spectrofluorometer at 20°C. The concentration of protein was 5 µM in 0.5 mM potassium phosphate buffer (pH 7.0). Different concentrations of GdmCl were obtained by successive addition of azurin denatured in 7 M GdmCl to a 0 M GdmCl solution of native azurin. Each point was incubated at 20°C for several minutes prior to collection of the emission.

An apo-WT equilibrium curve (data not shown) was measured by first incubating 30 µM Zn2+-coordinated WT azurin (the A280:A628 absorption ratio for this sample (which reflects the amount of bound copper) was zero, thus indicating 100% Zn2+) in 5.0 M GdmCl and 1.0 mM EDTA for ~20 min. The formed apo-azurin was then refolded by a 1:10 dilution to native conditions (0.5 M GdmCl) prior to titration by successive addition of GdmCl, with 1 mM EDTA present throughout the experiment. All holo-species were converted into the apo form by this treatment as judged from the obtained equilibrium curve. Dilution of the fluorescence intensity by the GdmCl additions was corrected for prior to curve fitting. The buffer was 0.5 mM potassium phosphate (pH 7.0), and the temperature was 20°C. Excitation was set to 295 nm, and the emission was collected at 308 nm.

Equilibrium CD data were obtained on a Jasco J-810 spectropolarimeter thermostated at 20°C. All points were prepared separately and incubated at room temperature overnight before data collection. Far-UV CD was collected at 220 nm using a 3-mm path-length cuvette, and near-UV CD was monitored at 281 nm using a 1-cm cuvette. The protein concentration was 2.5 µM for the far-UV measurements and 0.10 mM for the near-UV measurements. The buffer was 0.5 mM K phosphate (pH 7.0).

Parameters describing the equilibrium unfolding transitions were obtained as described (Santoro and Bolen 1988; Clarke and Fersht 1993) using the Igor Pro program (Wavemetrics).

The acid-induced molten globule-like states and their ability to bind ANS were analyzed in a SPEX Fluorolog {tau}2 spectrofluorometer thermostated at 20°C. The protein concentration was 5 µM, the ANS concentration was 0.5 mM, and the buffer was 0.5 mM potassium phosphate. Spectra of azurin were recorded by excitation at 295 nm and emission collected from 300 to 460 nm, and spectra of ANS were recorded by excitation at 388 nm and emission collected from 400 to 550 nm. The pH was set to 7.0 for the native azurin references, and to 2.6 for the acid-induced intermediate states.

Kinetic data collection and analysis
Kinetic unfolding and refolding traces were obtained using a BioLogic µSFM-20/MPS-52 stopped-flow apparatus. The instrumental dead time was 10.7 msec in all experiments. For the fluorescence, the observation head was connected to a Bio-Logic MOS-250 spectrometer, and the emission was recorded at 308 ± 2 nm with excitation set to 295 ± 2 nm. However, for the initial conditions test and the concentration-dependency test, the fluorescence was collected through an interference filter at 310 ± 2 nm (Melles Griot) using a Bio-Logic MOS-450 spectrometer coupled to the same stopped-flow instrument. CD data were collected at 220 nm also using the Bio-Logic MOS-450 spectrometer.

Refolding and unfolding of azurin monitored by fluorescence were analyzed by a 1:10 dilution of a 20 µM azurin solution to different end concentrations of denaturant. For the concentration-dependency test, the start concentration of protein was instead made 15, 20, 50, 80, 110, 140, 170, 200, 250, and 300 µM. Fluorescence refolding was initiated from a 3.0 M and a 2.0 M GdmCl solution for the C112S and C112A mutants, respectively, except for the initial conditions test of C112S, which was refolded in 0.5 M from different start conditions (4.0, 3.7, 3.4, 3.1, 2.9, 2.6, 2.3, 2.0, and 1.7 M GdmCl). The CD measurements were measured similarly, albeit with a starting concentration of 200 µM azurin. Prior to refolding, the denatured protein was incubated at room temperature for at least 2 h.

For the refolding experiments in the presence of ANS, both buffers to be mixed were made 0.5 mM ANS. The excitation was set to 388 nm, and the emission at 480 ± 5 nm was collected through an interference filter (Melles Griot) using the Bio-Logic MOS-450 configuration of the instrument. The end concentration of protein was 5 µM in these experiments, and the buffer was 0.5 mM potassium phosphate (pH 7.0).

The unfolding and refolding traces were averaged (typically five traces for the fluorescence and 15 traces for the CD data) and fitted to an exponential decay with one to three exponents using the Igor Pro program (Wavemetrics). The start values were obtained by extrapolation to time zero. Chevron plots (ln kobs versus [GdmCl]) were constructed, and a complete description of the observed rate constants for the major phases was obtained by least squares fitting of the following equation to the experimental data:


(3)

The curvature parameter Cf,u reflects the nonlinear dependence of the activation free energy of folding and unfolding on denaturant concentration. The Chevron profiles were also analyzed in terms of a U {leftrightarrows} I {leftrightarrows} N three-state model assuming a rapid U {leftrightarrows} I pre-equilibrium (Parker et al. 1995):


(4)

Here, kni denotes unfolding from the native state to the intermediate state, and kin denotes folding from the intermediate state to the native state. The m values associated with these transitions are similarly denoted. KU–I is equivalent to the rapid pre-equilibrium constant of the I {leftrightarrows} U transition. Because of this rapid pre-equilibrium, an on-pathway U {leftrightarrows} I {leftrightarrows} N model is kinetically indistinguishable from an off-pathway I {leftrightarrows} U {leftrightarrows} N model (Fersht 1999). However, a distinction between on- and off-pathway models can sometimes be made on the basis of how reasonable the ln kobs H2O is for an off-pathway model for which kobsH20 = KU–I x kinH2O (Capaldi et al. 1999). The curvature in the dependence of ln kobs on [D] is reflected in the C[D]2 parameter. This curvature was set equal in both unfolding and refolding limbs during the curve fitting, thereby ascribing additional curvature to refolding intermediates becoming accumulated.

The free energies of folding from the kinetic data were taken as –RTxln(ku/kf) or –RTxln(kni/kin) and the Tanford {beta}-values, reflecting buried surface area of the TS relative N, were calculated as:


(5)

For the three-state model, the same parameter was instead calculated as:


(6)

The corresponding parameter relating the intermediate (I) to N was obtained by:


(7)

All m parameters were taken at the equilibrium [D]50% as their respective first derivatives in folding free energies on denaturant concentration ({delta}{Delta}G/{delta} [D]) on account of the curvature in the Chevron plot (Fersht 1999).

The {phi}f value is a measure of the relative energetic sensitivity of the transition state and the unfolded state as resulting from the mutation (Fersht 1999). It is calculated as:


(8)

Here, {Delta}{Delta}GU–TS = RT ln kf/kf'(the prime denotes mutant parameter), where kf is the folding rate. The height of the TS barrier in Figure 8Go was estimated using the folding rate and the following equation (Fersht 1999):


(9)

where kB is the Boltzmann constant and h is the Planck constant.


    Acknowledgments
 
This work was funded by the Swedish Research Council and the Carl Trygger Foundation. Mats Ökvist at the Department of Chemistry, Göteborg University, Sweden, is acknowledged for the C112S and C112A PDB coordinates.

The publication costs of this article were defrayed in part by payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 USC section 1734 solely to indicate this fact.


    References
 TOP
 Abstract
 Introduction
 Results
 Discussion
 Materials and methods
 References
 
Anfinsen, C.B. 1973. Principles that govern the folding of protein chains. Science 181: 223–230.[Free Full Text]

Arai, M. and Kuwajima, K. 2000. Role of the molten globule state in protein folding. Adv. Protein Chem. 53: 209–282.[Medline]

Capaldi, A.P., Ferguson, S.J., and Radford, S.E. 1999. The Greek key protein apo-pseudoazurin folds through an obligate on-pathway intermediate. J. Mol. Biol. 286: 1621–1632.[CrossRef][Medline]

Chothia, C., Hubbard, T., Brenner, S., Barns, H., and Murzin, A. 1997. Protein folds in the all-{beta} and all-{alpha} classes. Annu. Rev. Biophys. Biomol. Struct. 26: 597–627.[CrossRef][Medline]

Clarke, J. and Fersht, A.R. 1993. Engineered disulfide bonds as probes of the folding pathway of barnase, increasing the stability of proteins against the rate of denaturation. Biochemistry 32: 4322–4329.[CrossRef][Medline]

Clarke, J., Cota, E., Fowler, S.B., and Hamill, S.J. 1999. Folding studies of immunoglobulin-like {beta}-sandwich proteins suggest that they share a common folding pathway. Structure Fold. Des. 7: 1145–1153.[Medline]

Denisov, V.P., Jonsson, B.H., and Halle, B. 1999. Hydration of denatured and molten globule proteins. Nat. Struct. Biol. 6: 253–260.[CrossRef][Medline]

Dill, K.A. and Shortle, D. 1991. Denatured states of proteins. Annu. Rev. Biochem. 60: 795–825.[CrossRef][Medline]

Eftink, M.R. 1991. Fluorescence techniques for studying protein structure. Methods Biochem. Anal. 35: 127–205.[CrossRef][Medline]

Fersht, A. 1999. Structure and mechanism in protein science. A guide to enzyme catalysis and protein folding. W.H. Freeman and Company, New York.

Fowler, S.B. and Clarke, J. 2001. Mapping the folding pathway of an immunoglobulin domain, structural detail from Phi value analysis and movement of the transition state. Structure (Camb.) 9: 355–366.

Gray, H.B., Malmstrom, B.G., and Williams, R.J. 2000. Copper coordination in blue proteins. J. Biol. Inorg. Chem. 5: 551–559.[CrossRef][Medline]

Haber, E. and Anfinsen, C.B. 1962. Side-chain interactions governing the pairing of half-cystine residues in ribonuclease. J. Biol. Chem. 237: 1839–1844.[Free Full Text]

Hoitink, C.W. 1993. "Engineering of the coppersite of azurin from Alcaligenes denitrificans." Ph.D. thesis, Leiden University, The Netherlands.

Jackson, S.E. 1998. How do small single-domain proteins fold? Fold. Des. 3: R81–R91.[CrossRef][Medline]

Jackson, S.E. and Fersht, A.R. 1991a. Folding of chymotrypsin inhibitor 2. 1. Evidence for a two-state transition. Biochemistry 30: 10428–10435.[CrossRef][Medline]

———. 1991b. Folding of chymotrypsin inhibitor 2. 2. Influence of proline isomerization on the folding kinetics and thermodynamic characterization of the transition state of folding. Biochemistry 30: 10436–10443.[CrossRef][Medline]

Karlsson, B.G., Pascher, T., Nordling, M., Arvidsson, R.H., and Lundberg, L.G. 1989. Expression of the blue copper protein azurin from Pseudomonas aeruginosa in Escherichia coli. FEBS Lett. 246: 211–217.[CrossRef][Medline]

Kiefhaber, T., Kohler, H.H., and Schmid, F.X. 1992. Kinetic coupling between protein folding and prolyl isomerization. I. Theoretical models. J. Mol. Biol. 224: 217–229.[CrossRef][Medline]

Kister, A.E., Finkelstein, A.V., and Gelfand, I.M. 2002. Common features in structures and sequences of sandwich-like proteins. Proc. Natl. Acad. Sci. 99: 14137–14141.[Abstract/Free Full Text]

Koide, S., Dyson, H.J., and Wright, P.E. 1993. Characterization of a folding intermediate of apoplastocyanin trapped by proline isomerization. Biochemistry 32: 12299–12310.[CrossRef][Medline]

Krantz, B.A., Mayne, L., Rumbley, J., Englander, S.W., and Sosnick, T.R. 2002. Fast and slow intermediate accumulation and the initial barrier mechanism in protein folding. J. Mol. Biol. 324: 359–371.[CrossRef][Medline]

Kraulis, P.J. 1991. Molscript—A program to produce both detailed and schematic plots of protein structures. J. Appl. Crystallogr. 24: 946–950.[CrossRef]

Levinthal, C. 1968. Are there pathways for protein folding? J. Chim. Phys. 65: 44–45.[CrossRef]

Makhatadze, G.I. 1999. Thermodynamics of protein interactions with urea and guanidinium hydrochloride. J. Phys. Chem. B 103: 4781–4785.[CrossRef]

Matouschek, A. and Fersht, A.R. 1993. Application of physical organic chemistry to engineered mutants of proteins, Hammond postulate behavior in the transition state of protein folding. Proc. Natl. Acad. Sci. 90: 7814–7818.[Abstract/Free Full Text]

Monera, O.D., Kay, C.M., and Hodges, R.S. 1994. Protein denaturation with guanidine hydrochloride or urea provides a different estimate of stability depending on the contributions of electrostatic interactions. Protein Sci. 3: 1984–1991.[Abstract]

Nozaki, Y. 1972. The preparation of guanidine hydrochloride. Methods Enzymol. 26: 43–50.

Parker, M.J., Spencer, J., and Clarke, A.R. 1995. An integrated kinetic analysis of intermediates and transition states in protein folding reactions. J. Mol. Biol. 253: 771–786.[CrossRef][Medline]

Plaxco, K.W., Simons, K.T., and Baker, D. 1998. Contact order, transition state placement and the refolding rates of single domain proteins. J. Mol. Biol. 277: 985–994.[CrossRef][Medline]

Pozdnyakova, I. and Wittung-Stafshede, P. 2003. Approaching the speed limit for Greek Key {beta}-barrel formation, transition-state movement tunes folding rate of zinc-substituted azurin. Biochim. Biophys. Acta 1651: 1–4.[Medline]

Pozdnyakova, I., Guidry, J., and Wittung-Stafshede, P. 2001. Copper stabilizes azurin by decreasing the unfolding rate. Arch. Biochem. Biophys. 390: 146–148.[CrossRef][Medline]

Ptitsyn, O.B. 1995. Molten globule and protein folding. Adv. Protein Chem. 47: 83–229.[Medline]

Reader, J.S., Van Nuland, N.A., Thompson, G.S., Ferguson, S.J., Dobson, C.M., and Radford, S.E. 2001. A partially folded intermediate species of the {beta}-sheet protein apo-pseudoazurin is trapped during proline-limited folding. Protein Sci. 10: 1216–1224.[Abstract/Free Full Text]

Richardson, J.S. 1977. {beta}-Sheet topology and the relatedness of proteins. Nature 268: 495–500.[CrossRef][Medline]

Sanchez, I.E. and Kiefhaber, T. 2003. Evidence for sequential barriers and obligatory intermediates in apparent two-state protein folding. J. Mol. Biol. 325: 367–376.[CrossRef][Medline]

Sandberg, A., Leckner, J., Shi, Y., Schwarz, F.P., and Karlsson, B.G. 2002. Effects of metal ligation and oxygen on the reversibility of the thermal denaturation of Pseudomonas aeruginosa azurin. Biochemistry 41: 1060–1069.[CrossRef][Medline]

Santoro, M.M. and Bolen, D.W. 1988. Unfolding free energy changes determined by the linear extrapolation method. 1. Unfolding of phenylmethanesulfonyl {alpha}-chymotrypsin using different denaturants. Biochemistry 27: 8063–8068.[CrossRef][Medline]

Semisotnov, G.V., Rodionova, N.A., Razgulyaev, O.I., Uversky, V.N., Gripas, A.F., and Gilmanshin, R.I. 1991. Study of the molten globule intermediate state in protein folding by a hydrophobic fluorescent probe. Biopolymers 31: 119–128.[CrossRef][Medline]

Silow, M. and Oliveberg, M. 1997a. High-energy channeling in protein folding. Biochemistry 36: 7633–7637.[CrossRef][Medline]

Silow, M. and Oliveberg, M. 1997b. Transient aggregates in protein folding are easily mistaken for folding intermediates. Proc. Natl. Acad. Sci. 94: 6084–6086.[Abstract/Free Full Text]

Solomon, E.I., Baldwin, M.J., and Lowery, M.D. 1992. Electronic-structures of active-sites in copper proteins—Contributions to reactivity. Chem. Rev. 92: 521–542.[CrossRef]

Solomon, E.I., LaCroix, L.B., and Randall, D.W. 1998. Electronic structure contributions to function in bioinorganic chemistry, the blue copper active site. Pure Appl. Chem. 70: 799–808.

Stryer, L. 1965. The interaction of a naphthalene dye with apomyoglobin and apohemoglobin. A fluorescent probe of nonpolar binding sites. J. Mol. Biol. 13: 482–495.[Medline]

Wallace, L.A. and Matthews, C.R. 2002. Sequential vs. parallel protein-folding mechanisms, experimental tests for complex folding reactions. Biophys. Chem. 101–102: 113–131.


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us   Add to Digg Digg   Add to Reddit Reddit   Add to Technorati Technorati    What's this?


This article has been cited by other articles:


Home page
Biophys. JHome page
K. Tveen Jensen, G. Strambini, M. Gonnelli, J. Broos, and J. B. Jackson
Mutations in Transhydrogenase Change the Fluorescence Emission State of TRP72 from 1La to 1Lb
Biophys. J., October 1, 2008; 95(7): 3419 - 3428.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
G. B. Strambini and M. Gonnelli
Protein Stability in Ice
Biophys. J., March 15, 2007; 92(6): 2131 - 2138.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow