Surface point mutations that significantly alter the structure and stability of a protein's denatured state

Cosolute and Crowding Effects on a Side-By-Side Protein Dimer

The effects of small (∼102 Da) and larger (>103 Da) cosolutes on the equilibrium stability of monomeric globular proteins are broadly understood, excluding volume stabilizes proteins and chemical interactions are stabilizing when repulsive, but destabilizing when attractive. Proteins, however, rarely work alone. Here, we investigate the effects of small and large cosolutes on the equilibrium stability of the simplest defined protein-protein interactions, the side-by-side homodimer formed by the A34F variant of the 56-residue B1 domain of protein G. We used 19F nuclear magnetic resonance spectroscopy to quantify the effects of urea, trimethylamine oxide, Ficoll, and more physiologically relevant cosolutes on the dimer dissociation constant. The data reveal the same stabilizing and destabilizing influences from chemical interactions as observed in studies of protein stability. Results with more physiologically relevant molecules such as bovine serum albumin, lysozyme, and reconstituted Escherichia coli cytosol reflect the importance of chemical interactions between these cosolutes and the test protein. Our study serves as a stepping-stone to a more complete understanding of crowding effects on protein-protein interactions.

Consistent View of Polypeptide Chain Expansion in Chemical Denaturants from Multiple Experimental Methods

There has been a long-standing controversy regarding the effect of chemical denaturants on the dimensions of unfolded and intrinsically disordered proteins: A wide range of experimental techniques suggest that polypeptide chains expand with increasing denaturant concentration, but several studies using small-angle X-ray scattering (SAXS) have reported no such increase of the radius of gyration (Rg). This inconsistency challenges our current understanding of the mechanism of chemical denaturants, which are widely employed to investigate protein folding and stability. Here, we use a combination of single-molecule Förster resonance energy transfer (FRET), SAXS, dynamic light scattering (DLS), and two-focus fluorescence correlation spectroscopy (2f-FCS) to characterize the denaturant dependence of the unfolded state of the spectrin domain R17 and the intrinsically disordered protein ACTR in two different denaturants. Standard analysis of the primary data clearly indicates an expansion of the unfolded state with increasing denaturant concentration irrespective of the protein, denaturant, or experimental method used. This is the first case in which SAXS and FRET have yielded even qualitatively consistent results regarding expansion in denaturant when applied to the same proteins. To more directly illustrate this self-consistency, we used both SAXS and FRET data in a Bayesian procedure to refine structural ensembles representative of the observed unfolded state. This analysis demonstrates that both of these experimental probes are compatible with a common ensemble of protein configurations for each denaturant concentration. Furthermore, the resulting ensembles reproduce the trend of increasing hydrodynamic radius with denaturant concentration obtained by 2f-FCS and DLS. We were thus able to reconcile the results from all four experimental techniques quantitatively, to obtain a comprehensive structural picture of denaturant-induced unfolded state expansion, and to identify the most likely sources of earlier discrepancies.

Small-angle X-ray scattering and single-molecule FRET spectroscopy produce highly divergent views of the low-denaturant unfolded state

The results of more than a dozen single-molecule Förster resonance energy transfer (smFRET) experiments suggest that chemically unfolded polypeptides invariably collapse from an expanded random coil to more compact dimensions as the denaturant concentration is reduced. In sharp contrast, small-angle X-ray scattering (SAXS) studies suggest that, at least for single-domain proteins at non-zero denaturant concentrations, such compaction may be rare. Here, we explore this discrepancy by studying protein L, a protein previously studied by SAXS (at 5 °C), which suggested fixed unfolded-state dimensions from 1.4 to 5 M guanidine hydrochloride (GuHCl), and by smFRET (at 25 °C), which suggested that, in contrast, the chain contracts by 15-30% over this same denaturant range. Repeating the earlier SAXS study under the same conditions employed in the smFRET studies, we observe little, if any, evidence that the unfolded state of protein L contracts as the concentration of GuHCl is reduced. For example, scattering profiles (and thus the shape and dimensions) collected within ∼4 ms after dilution to as low as 0.67 M GuHCl are effectively indistinguishable from those observed at equilibrium at higher denaturant. Our results thus argue that the disagreement between SAXS and smFRET is statistically significant and that the experimental evidence in favor of obligate polypeptide collapse at low denaturant cannot be considered conclusive yet.

How, when and why proteins collapse: the relation to folding

Unfolded proteins under strongly denaturing conditions are highly expanded. However, when the conditions are more close to native, an unfolded protein may collapse to a compact globular structure distinct from the folded state. This transition is akin to the coil-globule transition of homopolymers. Single-molecule FRET experiments have been particularly conducive in revealing the collapsed state under conditions of coexistence with the folded state. The collapse can be even more readily observed in natively unfolded proteins. Time-resolved studies, using FRET and small-angle scattering, have shown that the collapse transition is a very fast event, probably occurring on the submicrosecond time scale. The forces driving collapse are likely to involve both hydrophobic and backbone interactions. The loss of configurational entropy during collapse makes the unfolded state less stable compared to the folded state, thus facilitating folding.

Analysis of electrostatic interactions in the denatured state ensemble of the N-terminal domain of L9 under native conditions

The pH dependence of protein stability is defined by the difference in the number of protons bound to the folded state and to the denatured state ensemble (DSE) as a function of pH. In many cases, the protonation behavior can be described as the sum of a set of independently titrating residues; in this case, the pH dependence of stability reflects differences in folded and DSE pK(a)'s. pH dependent stability studies have shown that there are energetically important interactions involving charged residues in the DSE of the N-terminal domain of L9 (NTL9), which affect significantly the stability of the protein. The DSE of wild type NTL9 cannot be directly characterized under native conditions because of its high stability. A destabilized double mutant of NTL9, V3AI4A, significantly populates the folded state and the DSE in the absence of denaturant. The two states are in slow exchange on the nuclear magnetic resonance time scale, and diffusion measurements indicate that the DSE is compact. The DSE pK(a)'s of all of the acidic residues were directly determined. The DSE pK(a) of Asp8 and Asp23 are depressed relative to model compounds values. Use of the mutant DSE pK(a)'s together with known native state pK(a)'s leads to a significantly improved agreement between the measured pH dependent stability and that predicted by the Tanford-Wyman linkage relationship. An analysis of the literature suggests that DSE interactions involving charged residues are relatively common and should be considered in discussions of protein stability.

Role of solvation effects in protein denaturation: from thermodynamics to single molecules and back

Protein stability often is studied in vitro through the use of urea and guanidinium chloride, chemical cosolvents that disrupt protein native structure. Much controversy still surrounds the underlying mechanism by which these molecules denature proteins. Here we review current thinking on various aspects of chemical denaturation. We begin by discussing classic models of protein folding and how the effects of denaturants may fit into this picture through their modulation of the collapse, or coil-globule transition, which typically precedes folding. Subsequently, we examine recent molecular dynamics simulations that have shed new light on the possible microscopic origins of the solvation effects brought on by denaturants. It seems likely that both denaturants operate by facilitating solvation of hydrophobic regions of proteins. Finally, we present recent single-molecule fluorescence studies of denatured proteins, the analysis of which corroborates the role of denaturants in shifting the equilibrium of the coil-globule transition.

Molecular origin of constant m-values, denatured state collapse, and residue-dependent transition midpoints in globular proteins

Experiments show that for many two-state folders the free energy of the native state, DeltaG(ND)([C]), changes linearly as the denaturant concentration, [C], is varied. The slope {m = [dDeltaG(ND)([C])]/(d[C])}, is nearly constant. According to the transfer model, the m-value is associated with the difference in the surface area between the native (N) and denatured (D) state, which should be a function of DeltaR(g)(2), the difference in the square of the radius of gyration between the D and N states. Single-molecule experiments show that the R(g) of the structurally heterogeneous denatured state undergoes an equilibrium collapse transition as [C] decreases, which implies m also should be [C]-dependent. We resolve the conundrum between constant m-values and [C]-dependent changes in R(g) using molecular simulations of a coarse-grained representation of protein L, and the molecular transfer model, for which the equilibrium folding can be accurately calculated as a function of denaturant (urea) concentration. In agreement with experiment, we find that over a large range of denaturant concentration (>3 M) the m-value is a constant, whereas under strongly renaturing conditions (<3 M), it depends on [C]. The m-value is a constant above [C] > 3 M because the [C]-dependent changes in the surface area of the backbone groups, which make the largest contribution to m, are relatively narrow in the denatured state. The burial of the backbone and hydrophobic side chains gives rise to substantial surface area changes below [C] < 3 M, leading to collapse in the denatured state of protein L. Dissection of the contribution of various amino acids to the total surface area change with [C] shows that both the sequence context and residual structure are important. There are [C]-dependent variations in the surface area for chemically identical groups such as the backbone or Ala. Consequently, the midpoints of transition of individual residues vary significantly (which we call the Holtzer effect) even though global folding can be described as an all-or-none transition. The collapse is specific in nature, resulting in the formation of compact structures with appreciable populations of nativelike secondary structural elements. The collapse transition is driven by the loss of favorable residue-solvent interactions and a concomitant increase in the strength of intrapeptide interactions with a decreasing [C]. The strength of these interactions is nonuniformly distributed throughout the structure of protein L. Certain secondary structure elements have stronger [C]-dependent interactions than others in the denatured state.

Collapse transition in proteins

The coil-globule transition, a tenet of the physics of polymers, has been identified in recent years as an important unresolved aspect of the initial stages of the folding of proteins. We describe the basics of the collapse transition, starting with homopolymers and continuing with proteins. Studies of denatured-state collapse under equilibrium are then presented. An emphasis is placed on single-molecule fluorescence experiments, which are particularly useful for measuring properties of the denatured state even under conditions of coexistence with the folded state. Attempts to understand the dynamics of collapse, both theoretically and experimentally, are then described. Only an upper limit for the rate of collapse has been obtained so far. Improvements in experimental and theoretical methodology are likely to continue to push our understanding of the importance of the denatured-state thermodynamics and dynamics for protein folding in the coming years.

Electrostatic interactions in the denatured state ensemble: their effect upon protein folding and protein stability

It is now recognized that the denatured state ensemble (DSE) of proteins can contain significant amounts of structure, particularly under native conditions. Well-studied examples include small units of hydrogen bonded secondary structure, particularly helices or turns as well as hydrophobic clusters. Other types of interactions are less well characterized and it has often been assumed that electrostatic interactions play at most a minor role in the DSE. However, recent studies have shown that both favorable and unfavorable electrostatic interactions can be formed in the DSE. These can include surprisingly specific non-native interactions that can even persist in the transition state for protein folding. DSE electrostatic interactions can be energetically significant and their modulation either by mutation or by varying solution conditions can have a major impact upon protein stability. pH dependent stability studies have shown that electrostatic interactions can contribute up to 4 kcal mol(-1) to the stability of the DSE.

Folding pathway of the b1 domain of protein G explored by multiscale modeling

The understanding of the folding mechanisms of single-domain proteins is an essential step in the understanding of protein folding in general. Recently, we developed a mesoscopic CA-CB side-chain protein model, which was successfully applied in protein structure prediction, studies of protein thermodynamics, and modeling of protein complexes. In this research, this model is employed in a detailed characterization of the folding process of a simple globular protein, the B1 domain of IgG-binding protein G (GB1). There is a vast body of experimental facts and theoretical findings for this protein. Performing unbiased, ab initio simulations, we demonstrated that the GB1 folding proceeds via the formation of an extended folding nucleus, followed by slow structure fine-tuning. Remarkably, a subset of native interactions drives the folding from the very beginning. The emerging comprehensive picture of GB1 folding perfectly matches and extends the previous experimental and theoretical studies.

Protein folding kinetics provides a context-independent assessment of beta-strand propensity in the Fyn SH3 domain

Structural database-derived propensities for amino acids to adopt particular local protein structures, such as alpha-helix and beta-strand, have long been recognized and effectively exploited for the prediction of protein secondary structure. However, the experimental verification of database-derived propensities using mutagenesis studies has been problematic, especially for beta-strand propensities, because local structural preferences are often confounded by non-local interactions arising from formation of the native tertiary structure. Thus, the overall thermodynamic stability of a protein is not always altered in a predictable manner by changes in local structural propensity at a single position. In this study, we have undertaken an investigation of the relationship between beta-strand propensity and protein folding kinetics. By characterizing the effects of a wide variety of amino acid substitutions at two different beta-strand positions in an SH3 domain, we have found that the observed changes in protein folding rates are very well correlated to beta-strand propensities for almost all of the substitutions examined. In contrast, there is little correlation between propensities and unfolding rates. These data indicate that beta-strand conformation is well formed in the structured portion of the SH3 domain transition state, and that local structure propensity strongly influences the stability of the transition state. Since the transition state is known to be packed more loosely than the native state and likely lacks many of the non-local stabilizing interactions seen in the native state, we suggest that folding kinetics studies may generally provide an effective means for the experimental validation of database-derived local structural propensities.

Coil-globule transition in the denatured state of a small protein

Upon transfer from strongly denaturing to native conditions, proteins undergo a collapse that either precedes folding or occurs simultaneously with it. This collapse is similar to the well known coil-globule transition of polymers. Here we employ single-molecule fluorescence methods to fully characterize the equilibrium coil-globule transition in the denatured state of the IgG-binding domain of protein L. By using FRET measurements on freely diffusing individual molecules, we determine the radius of gyration of the protein, which shows a gradual expansion as the concentration of the denaturant, guanidinium hydrochloride, is increased all the way up to 7 M. This expansion is observed also in fluorescence correlation spectroscopy measurements of the hydrodynamic radius of the protein. We analyze the radius of gyration measurements using the theory of the coil-globule transition of Sanchez [Sanchez, I. C. (1979) Macromolecules 12, 980-988], which balances the excluded volume entropy of the chain with the average interresidue interaction energy. In particular, we calculate the solvation energy of the denatured protein, a property that is not readily accessible in other experiments. The dependence of this energy on denaturant concentration is nonlinear, contrasting with the common linear extrapolation method used to describe denaturation energy. Interestingly, a fit to the binding model of chemical denaturation suggests a single denaturant binding site per protein residue. The size of the denatured protein under native conditions can be extrapolated from the data as well, showing that the fully collapsed state of protein is only approximately 10% larger than the folded state.

Characterization of Intra-molecular Distances and Site-specific Dynamics in Chemically Unfolded Barstar: Evidence for Denaturant-dependent Non-random Structure

The structure and dynamics of the unfolded form of a protein are expected to play critical roles in determining folding pathways. In this study, the urea and guanidine hydrochloride (GdnHCl)-unfolded forms of the small protein barstar were explored by time-resolved fluorescence techniques. Barstar was labeled specifically with thionitrobenzoate (TNB), by coupling it to the thiol side-chain of a cysteine residue at one of the following positions on the sequence: 14, 25, 40, 42, 62, 82 and 89, in single cysteine-containing mutant proteins. Seven intra-molecular distances (R(DA)) under unfolding conditions were estimated from measurements of time-resolved fluorescence resonance energy transfer between the donor Trp53 and the non-fluorescent acceptor TNB coupled to one of the seven cysteine side-chains. The unfolded protein chain expands with an increase in the concentration of the denaturants. The extent of expansion was found to be non-uniform, with different intra-molecular distances expanding to different extents. In general, shorter distances were found to expand less when compared to longer spans. The extent of expansion was higher in the case of GdnHCl when compared to urea. A comparison of the measured values of R(DA) with those derived from a model based on excluded volume, revealed that while shorter spans showed good agreement, the experimental values of R(DA) of longer spans were smaller when compared to the theoretical values. Sequence-specific flexibility of the polypeptide was determined by time-resolved fluorescence anisotropy decay measurements on acrylodan or 1,5-IAEDANS labeled single cysteine-containing proteins under unfolding conditions. Rotational dynamics derived from these measurements indicated that the level of flexibility increased with increase in the concentration of denaturants and showed a graded increase towards the C-terminal end. Taken together, these results appear to indicate the presence of specific non-random coil structures and show that the deviation from random coil structure is different for the two denaturants.

Random-coil behavior and the dimensions of chemically unfolded proteins

Spectroscopic studies have identified a number of proteins that appear to retain significant residual structure under even strongly denaturing conditions. Intrinsic viscosity, hydrodynamic radii, and small-angle x-ray scattering studies, in contrast, indicate that the dimensions of most chemically denatured proteins scale with polypeptide length by means of the power-law relationship expected for random-coil behavior. Here we further explore this discrepancy by expanding the length range of characterized denatured-state radii of gyration (R(G)) and by reexamining proteins that reportedly do not fit the expected dimensional scaling. We find that only 2 of 28 crosslink-free, prosthetic-group-free, chemically denatured polypeptides deviate significantly from a power-law relationship with polymer length. The R(G) of the remaining 26 polypeptides, which range from 16 to 549 residues, are well fitted (r(2) = 0.988) by a power-law relationship with a best-fit exponent, 0.598 +/- 0.028, coinciding closely with the 0.588 predicted for an excluded volume random coil. Therefore, it appears that the mean dimensions of the large majority of chemically denatured proteins are effectively indistinguishable from the mean dimensions of a random-coil ensemble.

Dramatic acceleration of protein folding by stabilization of a nonnative backbone conformation

Through a mutagenic investigation of Gly-48, a highly conserved position in the Src homology 3 domain, we have discovered a series of amino acid substitutions that are highly destabilizing, yet dramatically accelerate protein folding, some up to 10-fold compared with the wild-type rate. The unique folding properties of these mutants allowed for accurate measurement of their folding and unfolding rates in water with no denaturant by using an NMR spin relaxation dispersion technique. A strong correlation was found between beta-sheet propensity and the folding rates of the Gly-48 mutants, even though Gly-48 lies in an unusual non-beta-strand backbone conformation in the native state. This finding indicates that the accelerated folding rates of the Gly-48 mutants are the result of stabilization of a nonnative beta-strand conformation in the transition-state structure at this position, thus providing the first, to our knowledge, experimentally elucidated example of a mechanism by which folding can occur fastest through a nonnative conformation. We also demonstrate that residues that are most stabilizing in the transition-state structure are most destabilizing in the native state, and also cause the greatest reductions in in vitro functional activity. These data indicate that the unusual native conformation of the Gly-48 position is important for function, and that evolutionary selection for function can result in a domain that folds at a rate far below the maximum possible.

Non-linear rate-equilibrium free energy relationships and Hammond behavior in protein folding

Non-linear rate-equilibrium relationships upon mutation or changes in solvent conditions are frequently observed in protein folding reactions and are usually interpreted in terms of Hammond behavior. Here we first give a general overview over the concept of transition state movements in chemical reactions and discuss its application to protein folding. We then show examples for genuine Hammond behavior and for apparent transition state movements caused by other effects like changes in the rate-limiting step of the folding reaction or ground state effects, i.e. structural changes in either the native state or the unfolded state. These examples show that apparent transition state movements can easily be mistaken for Hammond behavior. We describe experimental tests using self- and cross-interaction parameters to distinguish between structural changes in a single transition state following Hammond behavior and apparent transition state movements caused by other effects.

Extended application of gel-permeation chromatography by spin column

Separation of small volumes of proteins from unbound ligands or reequilibration with buffer by passing through a 1-ml Sephadex G-50 column under mild centrifugal force is a popular technique. Here it has been demonstrated that other Sephadex matrix could similarly be used for complete or partial separation of protein molecules. Proteins to be eluted at void volume are recovered near quantitatively, while others are partly or almost completely retained depending on molecular size. Calibration curves using standard proteins of Mw 12.5 to 440 kDa with Sephadex G-50-G-200 representing recovery versus molecular weight show profiles as expected from the fractionation ranges of the column matrix. The procedure may be applied to follow protein association-dissociation reactions if the molecular weights of the species concerned are known and a proper matrix exists for separating them. Equilibrium unfolding transitions constructed with model proteins in presence of 0-8M urea using recovery as an index correspond to profiles obtained from other physical measurements. This may be a convenient approach to follow change of protein hydrodynamic volume quickly when a parallel methodology is not readily available.

Prediction of folding rates and transition-state placement from native-state geometry

A variety of experimental and theoretical studies have established that the folding process of monomeric proteins is strongly influenced by the topology of the native state. In particular, folding times have been shown to correlate well with the contact order, a measure of contact locality. Our investigation focuses on identifying additional topologic properties that correlate with experimentally measurable quantities, such as folding rates and transition-state placement, for both two- and three-state folders. The validation against data from 40 experiments shows that a particular topological property that measures the interdependence of contacts, termed cliquishness or clustering coefficient, can account with statistically significant accuracy both for the transition state placement and especially for folding rates. The observed correlations can be further improved by optimally combining the distinct topological information captured by cliquishness and contact order.

Computational simulation of the statistical properties of unfolded proteins

A simple Monte Carlo method was used to generate ensembles of simulated polypeptide conformations that are restricted only by steric repulsion. The models used for these simulations were based on the sequences of four real proteins, ranging in size from 26 to 268 amino acid residues, and included all non-hydrogen atoms. Two sets of calculations were performed, one that included only intra-residue steric repulsion terms and those between adjacent residues, and one that included repulsion terms between all possible atom pairs, so as to explicitly account for the excluded volume effect. Excluded volume was found to increase the average radius of gyration of the chains by 20-40%, with the expansion factor increasing with chain length. Contrary to recent suggestions, however, the excluded volume effect did not greatly restrict the distribution of dihedral angles or favor native-like topologies. The average dimensions of the ensembles calculated with excluded volume were consistent with those measured experimentally for unfolded proteins of similar sizes under denaturing conditions, without introducing any adjustable scaling factor. The simulations also reproduced experimentally determined effective concentrations for the formation of disulfide bonds in reduced and unfolded proteins. The statistically generated ensembles included significant numbers of conformations that were nearly as compact as the corresponding native proteins, as well as many that were as accessible to solvent as a fully extended chain. On the other hand, conformations with as much buried surface area as the native proteins were very rare, as were highly extended conformations. These results suggest that the overall properties of unfolded proteins can be usefully described by a random coil model and that an unfolded polypeptide can undergo significant collapse while losing only a relatively small fraction of its conformational entropy.

Exploring the potential of the monobody scaffold: effects of loop elongation on the stability of a fibronectin type III domain

The tenth fibronectin type III domain of human fibronectin (FNfn10) is a small, monomeric beta-sandwich protein, similar to the immunoglobulins. We have developed small antibody mimics, 'monobodies', using FNfn10 as a scaffold. We initially altered two loops of FNfn10 that are structurally equivalent to two of the hypervariable loops of the immunoglobulin domain. In order to assess the possibility of utilizing other loops in FNfn10 for target binding, we determined the effects of the elongation of each loop on the conformational stability of FNfn10. We found that all six loops of FNfn10 allowed the introduction of four glycine residues while retaining the global fold. Insertions in the AB and FG loops exhibited very small degrees of destabilization, comparable to or less than predicted entropic penalties due to the elongation, suggesting the absence of stabilizing interactions in these loops in wild-type FNfn10. Insertions in the BC, CD and DE loops, respectively, resulted in modest destabilization. In contrast, the EF loop elongation was highly destabilizing, consistent with previous studies showing the presence of stabilizing interactions in this loop. These results suggest that all loops, except for the EF loop, can be used for engineering a binding site, thus demonstrating excellent properties of the monobody scaffold.

Glutamine 53 is a gatekeeper residue in the FK506 binding protein

The effect of non-random conformational averaging in the urea-unfolded state on the folding pathway has been investigated in a variant of the FK506 binding protein with three additional residues at the amino terminus (FKBP(*)). Three mutations (asparagine, aspartate, and threonine) were introduced into position Q53 to enhance formation of non-native helix observed in this part of the protein in the urea-unfolded state. NMR analysis showed minor structural changes in the native state of each mutant, but additional medium-range alphaN(i,i+2) of each mutant nuclear Overhauser enhancements were observed in the urea-unfolded state that were not in FKBP(*), indicating that the mutations had a more substantial effect on the unfolded state ensemble than on the native state ensemble. Isothermal equilibrium denaturation measurements showed that the Q53T and Q53D mutants were destabilized, whereas the Q53N mutant was stabilized relative to FKBP(*) with little change in the equilibrium m values. The unfolding rates of Q53N and Q53T were similar to that of FKBP(*), but Q53D unfolded twice as fast as FKBP(*). In contrast, the mutations had a more pronounced effect on the refolding kinetics. Q53N refolded slightly faster and exhibited a kinetic folding intermediate similar to that of FKBP(*). The Q53D and Q53T mutants also refolded faster than FKBP(*) but lacked the folding intermediate, indicating that these mutants experienced a different folding trajectory and transition state than FKBP(*) and Q53N. The refolding kinetic Phi values were 0.74, 1.4 and 7.9 for Q53N, Q53T, and Q53D, respectively. The data point to Q53 functioning as a gatekeeper residue in the folding of FKBP(*). This study shows that perturbing the unfolded state ensemble via mutagenesis can provide insights into residues that play important roles in the folding pathway, and represents an attractive strategy for mapping the high-energy portions of the folding energy landscape.

Model for calculation of electrostatic interactions in unfolded proteins

An approach for the calculation of electrostatic interactions and titration properties of unfolded polypeptide chains (denatured proteins) is proposed. It is based on a simple representation of the denatured proteins as a state with titratable sites distributed on the surface of a sphere, radius of which is assumed to be equal to the radius of gyration, R(g), of an unfolded molecule. Distances between the charges, d, obey constraints arising from the protein sequence. Criteria for evaluation of the parameters R(g) and d were obtained from computer simulations on a polypeptide consisting of 20 identical amino acids (polylysine). The model was applied for calculation of titration curves of denatured barnase and staphylococcal nuclease. It was demonstrated that the approach proposed gives considerably better agreement with the experimental data than the commonly used null approximation. It was also found that titration properties of denatured proteins are slightly, but distinguishably influenced by the amino-acid sequence of the protein.

NMR and SAXS characterization of the denatured state of the chemotactic protein CheY: implications for protein folding initiation

The denatured state of a double mutant of the chemotactic protein CheY (F14N/V83T) has been analyzed in the presence of 5 M urea, using small angle X-ray scattering (SAXS) and heteronuclear magnetic resonance. SAXS studies show that the denatured protein follows a wormlike chain model. Its backbone can be described as a chain composed of rigid elements connected by flexible links. A comparison of the contour length obtained for the chain at 5 M urea with the one expected for a fully expanded chain suggests that approximately 25% of the residues are involved in residual structures. Conformational shifts of the alpha-protons, heteronuclear (15)N-[(1)H] NOEs and (15)N relaxation properties have been used to identify some regions in the protein that deviate from a random coil behavior. According to these NMR data, the protein can be divided into two subdomains, which largely coincide with the two folding subunits identified in a previous kinetic study of the folding of the protein. The first of these subdomains, spanning residues 1-70, is shown here to exhibit a restricted mobility as compared to the rest of the protein. Two regions, one in each subdomain, were identified as deviating from the random coil chemical shifts. Peptides corresponding to these sequences were characterized by NMR and their backbone (1)H chemical shifts were compared to those in the intact protein under identical denaturing conditions. For the region located in the first subdomain, this comparison shows that the observed deviation from random coil parameters is caused by interactions with the rest of the molecule. The restricted flexibility of the first subdomain and the transient collapse detected in that subunit are consistent with the conclusions obtained by applying the protein engineering method to the characterization of the folding reaction transition state.

Amino-acid substitutions at the fully exposed P1 site of bovine pancreatic trypsin inhibitor affect its stability

It is widely accepted that solvent-exposed sites in proteins play only a negligible role in determining protein energetics. In this paper we show that amino acid substitutions at the fully exposed Lys15 in bovine pancreatic trypsin inhibitor (BPTI) influenced the CD- and DSC-monitored stability: The T(den) difference between the least (P1 Trp) and the most stable (P1 His) mutant is 11.2 degrees C at pH 2.0. The DeltaH(den) versus T(den) plot for all the variants at three pH values (2.0, 2.5, 3.0) is linear (DeltaC(p,den) = 0.41 kcal* mole(-1) * K(-1); 1 cal = 4.18 J) leading to a DeltaG(den) difference of 2.1 kcal*mole(-1). Thermal denaturation of the variants monitored by CD signal at pH 2.0 in the presence of 6 M GdmCl again showed differences in their stability, albeit somewhat smaller (DeltaT(den) =7.1 degrees C). Selective reduction of the Cys14-Cys 38 disulfide bond, which is located in the vicinity of the P1 position did not eliminate the stability differences. A correlation analysis of the P1 stability with different properties of amino acids suggests that two mechanisms may be responsible for the observed stability differences: the reverse hydrophobic effect and amino acid propensities to occur in nonoptimal dihedral angles adopted by the P1 position. The former effect operates at the denatured state level and causes a drop in protein stability for hydrophobic side chains, due to their decreased exposure upon denaturation. The latter factor influences the native state energetics and results from intrinsic properties of amino acids in a way similar to those observed for secondary structure propensities. In conclusion, our results suggest that the protein-stability-derived secondary structure propensity scales should be taken with more caution.

Mechanisms of protein folding

The strong correlation between protein folding rates and the contact order suggests that folding rates are largely determined by the topology of the native structure. However, for a given topology, there may be several possible low free energy paths to the native state and the path that is chosen (the lowest free energy path) may depend on differences in interaction energies and local free energies of ordering in different parts of the structure. For larger proteins whose folding is assisted by chaperones, such as the Escherichia coli chaperonin GroEL, advances have been made in understanding both the aspects of an unfolded protein that GroEL recognizes and the mode of binding to the chaperonin. The possibility that GroEL can remove non-native proteins from kinetic traps by unfolding them either during polypeptide binding to the chaperonin or during the subsequent ATP-dependent formation of folding-active complexes with the co-chaperonin GroES has also been explored.

Modulating protein folding rates in vivo and in vitro by side-chain interactions between the parallel beta strands of green fluorescent protein

We have identified pairs of residues across the two parallel beta strands of green fluorescent protein that facilitate native strand register of the surface-exposed beta barrel. After constructing a suitable host environment around two guest residues, minimizing interactions of the guest residues with surrounding side-chains yet maintaining the wild-type protein structure and the chromophore environment, we introduced a library of cross-strand pairings by cassette mutagenesis. Colonies of Escherichia coli transformed with the library differ in intracellular fluorescence. Most of the fluorescent pairs have predominantly charged and polar guest site residues. The magnitude and the rate of fluorescence acquisition in vivo from transformed E. coli cells varies among the mutants despite comparable levels of protein expression. Spectroscopic measurements of purified mutants show that the native protein structure is maintained. Kinetic studies using purified protein with fully matured chromophores demonstrate that the mutants span a 10-fold range in folding rates with undetectable differences in unfolding rates. Thus, green fluorescent protein provides an ideal system for monitoring determinants of in vivo protein folding. Cross-strand pairings affect both protein stability and folding kinetics by favoring the formation of native strand register preferentially to non-native strand alignments.

Effects of side-chain characteristics on stability and oligomerization state of a de novo-designed model coiled-coil: 20 amino acid substitutions in position "d"

We describe the de novo design and biophysical characterization of a model coiled-coil protein in which we have systematically substituted 20 different amino acid residues in the central "d" position. The model protein consists of two identical 38 residue polypeptide chains covalently linked at their N termini via a disulfide bridge. The hydrophobic core contained Val and Ile residues at positions "a" and Leu residues at positions "d". This core allowed for the formation of both two-stranded and three-stranded coiled-coils in benign buffer, depending on the substitution at position "d". The structure of each analog was analyzed by CD spectroscopy and their relative stability determined by chemical denaturation using GdnHCI (all analogs denatured from the two-stranded state). The oligomeric state(s) was determined by high-performance size-exclusion chromatography and sedimentation equilibrium analysis in benign medium. Our results showed a thermodynamic stability order (in order of decreasing stability) of: Leu, Met, Ile, Tyr, Phe, Val, Gln, Ala, Trp, Asn, His, Thr, Lys, Ser, Asp, Glu, Arg, Orn, and Gly. The Pro analog prevented coiled-coil formation. The overall stability range was 7.4 kcal/mol from the lowest to the highest analog, indicating the importance of the hydrophobic core and the dramatic effect a single substitution in the core can have upon the stability of the protein fold. In general, the side-chain contribution to the level of stability correlated with side-chain hydrophobicity. Molecular modelling studies, however, showed that packing effects could explain deviations from a direct correlation. In regards to oligomerization state, eight analogs demonstrated the ability to populate exclusively one oligomerization state in benign buffer (0.1 M KCl, 0.05 M K(2)PO(4)(pH 7)). Ile and Val (the beta-branched residues) induced the three-stranded oligomerization state, whereas Tyr, Lys, Arg, Orn, Glu and Asp induced the two-stranded state. Asn, Gln, Ser, Ala, Gly, Phe, Leu, Met and Trp analogs were indiscriminate and populated two-stranded and three-stranded states. Comparison of these results with similar substitutions in position "a" highlights the positional effects of individual residues in defining the stability and numbers of polypeptide chains occurring in a coiled-coil structure. Overall, these results in conjunction with other work now generate a relative thermodynamic stability scale for 19 naturally occurring amino acid residues in either an "a" or "d" position of a two-stranded coiled-coil. Thus, these results will aid in the de novo design of new coiled-coil structures, a better understanding of their structure/function relationships and the design of algorithms to predict the presence of coiled-coils within native protein sequences.

First principles prediction of protein folding rates

Experimental studies have demonstrated that many small, single-domain proteins fold via simple two-state kinetics. We present a first principles approach for predicting these experimentally determined folding rates. Our approach is based on a nucleation-condensation folding mechanism, where the rate-limiting step is a random, diffusive search for the native tertiary topology. To estimate the rates of folding for various proteins via this mechanism, we first determine the probability of randomly sampling a conformation with the native fold topology. Next, we convert these probabilities into folding rates by estimating the rate that a protein samples different topologies during diffusive folding. This topology-sampling rate is calculated using the Einstein diffusion equation in conjunction with an experimentally determined intra-protein diffusion constant. We have applied our prediction method to the 21 topologically distinct small proteins for which two-state rate data is available. For the 18 beta-sheet and mixed alpha-beta native proteins, we predict folding rates within an average factor of 4, even though the experimental rates vary by a factor of approximately 4 x 10(4). Interestingly, the experimental folding rates for the three four-helix bundle proteins are significantly underestimated by this approach, suggesting that proteins with significant helical content may fold by a faster, alternative mechanism. This method can be applied to any protein for which the structure is known and hence can be used to predict the folding rates of many proteins prior to experiment.

A thermodynamic analysis of a family of small globular proteins: SH3 domains

The stability and folding thermodynamics of two SH3-domains, belonging to Fyn and Abl proteins, have been studied by scanning calorimetry and urea-induced unfolding. They undergo an essentially two-state unfolding with parameters similar to those of the previously studied alpha-spectrin SH3 domain. The correlations between the thermodynamic parameters (heat capacity increment, delta Cp,U, the proportionality factor, m, and the Gibbs energy, delta Gw298) of unfolding and some integral structural parameters, such as polar and non-polar areas exposed upon domain denaturation, have been analyzed. The experimental data on delta Cp,U and the m-factor of the linear extrapolation model (LEM) obey the simple empirical correlations deduced elsewhere. The Gibbs energies calculated from the DSC data were compared with those found by fitting urea-unfolding curves to the LEM and the denaturant-binding model (DBM). The delta Gw298 values found with DBM correlate better with the DSC data, while those obtained with LEM are systematically smaller. The systematic difference between the parameters calculated with LEM and DBM are explained by an inherent imperfection of the LEM.

Conformations of peptide fragments from the FK506 binding protein: comparison with the native and urea-unfolded states

The helix-forming tendency of seven peptide fragments corresponding with the entire sequence of the FK506 binding protein (FKBP) has been investigated in aqueous buffer and in 2,2,2-trifluoroethanol (TFE) using CD and NMR spectroscopy. All fragments exhibited random coil conformations in aqueous buffer, whereas the amount of helix induced in the peptide fragments by TFE varied. The fragment with the highest degree of helicity in TFE corresponded with the single (alpha-helix in native FKBP. Fragments corresponding with beta-strands 2 and 3 also exhibited strong propensity towards helix formation. In contrast, the fragment corresponding with beta-strand 1 did not form helix in TFE. The inherent helix-forming tendencies are interpreted in light of the native structure to suggest possible folding nucleation sites. Conformational sampling in each peptide fragment was also compared with that observed in urea-denatured FKBP. With the exception of the fragment corresponding with beta-strand 2, the formation of helical structures in the peptide fragments in TFE was correlated with the observation of turn and/or helix conformers in urea-unfolded FKBP. Surprisingly, peptide fragments in aqueous solution were less structured than the corresponding regions in urea-denatured FKBP. The conformational differences between the peptide fragments and unfolded FKBP were not due to the urea buffer or to differences in their rotational correlation times. We conclude that local amino acid interactions are not generally sufficient to account for the formation of non-random conformations in unfolded FKBP. Formation of non-random structures in unfolded FKBP may require stabilization of incipient turn or helical conformations through transient contact with non-local non-polar residues.

Structure of the transition state in the folding process of human procarboxypeptidase A2 activation domain

The transition state for the folding pathway of the activation domain of human procarboxypeptidase A2 (ADA2h) has been analyzed by the protein engineering approach. Recombinant ADA2h is an 81-residue globular domain with no disulfide bridges or cis-prolyl bonds, which follows a two-state folding transition. Its native fold is arranged in two alpha-helices packing against a four-stranded beta-sheet. Application of the protein engineering analysis for 20 single-point mutants spread throughout the whole sequence indicates that the transition state for this molecule is quite compact, possessing some secondary structure and a hydrophobic core in the process of being consolidated. The core (folding nucleus) is made by the packing of alpha-helix 2 and the two central beta-strands. The other two strands, at the edges of the beta-sheet, and alpha-helix 1 seem to be completely unfolded. These results, together with previous analysis of ADA2h with either of its two alpha-helices stabilized through improved local interactions, suggest that alpha-helix 1 does not contribute to the folding nucleus, even though it is partially folded in the denatured state under native conditions. On the other hand, alpha-helix 2 folds partly in the transition state and is part of the folding nucleus. It is suggested that a good strategy to improve folding speed in proteins would be to stabilize the helices that are not folded in the denatured state but are partly present in the transition state. Comparison with other proteins shows that there is no clear relationship between fold and/or size with folding speed and level of structure in the transition state of proteins.

Synergy between simulation and experiment in describing the energy landscape of protein folding

Experimental data from protein engineering studies and NMR spectroscopy have been used by theoreticians to develop algorithms for helix propensity and to benchmark computer simulations of folding pathways and energy landscapes. Molecular dynamic simulations of the unfolding of chymotrypsin inhibitor 2 (CI2) have provided detailed structural models of the transition state ensemble for unfolding/folding of the protein. We now have used the simulated transition state structures to design faster folding mutants of CI2. The models pinpoint a number of unfavorable local interactions at the carboxyl terminus of the single alpha-helix and in the protease-binding loop region of CI2. By removing these interactions or replacing them with stabilizing ones, we have increased the rate of folding of the protein up to 40-fold (tau = 0.4 ms). This correspondence, and other examples of agreement between experiment and theory in general, Phi-values and molecular dynamics simulations, in particular, suggest that significant progress has been made toward describing complete folding pathways at atomic resolution by combining experiment and simulation.

Formation of a denatured dimer limits the thermal stability of Arc repressor

The thermal stability of the Arc repressor dimer normally increases with concentration because protein folding and subunit association are thermodynamically coupled. At Arc concentrations above 100 microM, however, thermal denaturation remains reversible and cooperative but tm does not continue to increase. In this concentration regime, thermally denatured Arc shows significantly reduced secondary structure and no evidence of a tightly packed core, but light scattering and fluorescence polarization studies indicate that the protein is dimeric. Higher order denatured oligomers are not observed and the stability of the non-native dimer is reduced by Arc mutations, indicating that non-native dimerization involves specific interactions between Arc subunits.

Insights into biomolecular function from small-angle scattering

Recent advances in neutron and X-ray sources and instrumentation, new and improved scattering techniques, and molecular biology techniques, which have permitted facile preparation of samples, have each led to new opportunities in using small-angle scattering to study the conformations and interactions of biological macromolecules in solution as a function of their properties. For example, new instrumentation on synchrotron sources has facilitated time-resolved studies that yield insights into protein folding. More powerful neutron sources, combined with molecular biology tools that isotopically label samples, have facilitated studies of biomolecular interactions, including those involving active enzymes.

The domains in gammaB-crystallin: identical fold-different stabilities

gammaB-crystallin from vertebrate eye lens is an all beta-sheet two-domain protein with a high degree of intrachain symmetry. Its N and C-terminal domains show high levels of sequence similarity and structural identity. In natural gammaB-crystallin, the domains fold independently. The recombinantly expressed isolated domains are stable monomeric proteins, which do not associate spontaneously to form a gammaB-like dimer. In contrast to their identical folding topology, the two domains obviously follow different folding mechanisms. While the two-state model is valid for the C-terminal domain, the folding behaviour of the N-terminal domain is more complex. The stability of the C-terminal domain is strongly dependent on pH. At pH 2, the C-terminal domain in its isolated form is significantly less stable than within the gammaB-molecule. In contrast, the isolated N-terminal domain does not differ in its stability from the N-terminal domain in wild-type gammaB-crystallin. The strongly decreased stability of the C-terminal domain at acid pH allowed a dissection of the intrinsic stabilities of the domains and their interactions in gammaB-crystallin. At pH 2, domain interactions contribute -16 kJ/mol to the overall stability of gammaB-crystallin.