Non-native local interactions in protein folding and stability: introducing a helical tendency in the all beta-sheet alpha-spectrin SH3 domain

Roles of beta-turns in protein folding: from peptide models to protein engineering

Reverse turns are a major class of protein secondary structure; they represent sites of chain reversal and thus sites where the globular character of a protein is created. It has been speculated for many years that turns may nucleate the formation of structure in protein folding, as their propensity to occur will favor the approximation of their flanking regions and their general tendency to be hydrophilic will favor their disposition at the solvent-accessible surface. Reverse turns are local features, and it is therefore not surprising that their structural properties have been extensively studied using peptide models. In this article, we review research on peptide models of turns to test the hypothesis that the propensities of turns to form in short peptides will relate to the roles of corresponding sequences in protein folding. Turns with significant stability as isolated entities should actively promote the folding of a protein, and by contrast, turn sequences that merely allow the chain to adopt conformations required for chain reversal are predicted to be passive in the folding mechanism. We discuss results of protein engineering studies of the roles of turn residues in folding mechanisms. Factors that correlate with the importance of turns in folding indeed include their intrinsic stability, as well as their topological context and their participation in hydrophobic networks within the protein's structure.

Tricalcium Phosphate Delayed Release Formulation for Oral Delivery of Insulin: A Proof-of-Concept Study

Several attempts have been made for delivering insulin orally utilizing several polymers with varying degrees of effectiveness. A major obstacle associated with polymeric delivery system for protein or polypeptide drugs is the poor retention of the structure and its biological activity of encapsulated proteins particularly for the unstable insulin. Calcium phosphate ceramic is considered highly compatible to protein or peptide drugs, particularly insulin. Therefore, an attempt has been made to load insulin in tricalcium phosphate (TCP) microspheres and coat with pH sensitive polymer of methacrylate derivative, and to study the stability and conformational variations of loaded insulin, and finally the biological activity of the formulation in diabetic rats. TCP microspheres were prepared by a standard procedure. Human insulin was loaded in to these porous microspheres by diffusion filling and coated with Eudragit S100. This was subjected to in vitro release studies in simulated fluids and the stability and conformational variations of the released insulin were studied using photon correlation spectroscopy and circular dichroism (CD). Biological activity of the formulation was studied on induced diabetic rats. Insulin released in the intestinal fluid (SIF) maintained the native conformation without showing any aggregation. A dose dependent reduction of blood glucose level (BGL) was achieved in streptozotocin induced diabetic Wistar rats, demonstrating its biological activity. It has been established from this preliminary study that insulin loaded in to TCP microspheres is highly compatible with no degradation or loss of biological activity of loaded insulin. The TCP microsphere based delayed release formulation of insulin has effected a decrease in elevated glucose level in induced diabetic rats, establishing its feasibility towards the development of a noninvasive delivery device.

Engineering Diverse Changes in β-Turn Propensities in the N-Terminal β-Hairpin of Ubiquitin Reveals Significant Effects on Stability and Kinetics but a Robust Folding Transition State

Using the N-terminal 17-residue beta-hairpin of ubiquitin as a "host" for mutational studies, we have investigated the influence of the beta-turn sequence on protein stability and folding kinetics by replacing the native G-bulged turn (TLTGK) with more flexible analogues (TG3K and TG5K) and a series of four-residue type I' beta-turn sequences, commonly found in beta-hairpins. Although a statistical analysis of type I' turns demonstrates residue preferences at specific sites, the frequency of occurrence appears to only broadly correlate with experimentally determined protein stabilities. The subsequent engineering of context-dependent non-native tertiary contacts involving turn residues is shown to produce large changes in stability. Relatively few point mutations have been described that probe secondary structure formation in ubiquitin in a manner that is independent of tertiary contacts. To this end, we have used the more rigorous rate-equilibrium free energy relationship (Leffler analysis), rather than the two-point phi value analysis, to show for a family of engineered beta-turn mutants that stability (range of approximately 20 kJ/mol) and folding kinetics (190-fold variation in refolding rate) are linearly correlated (alpha(f) = 0.74 +/- 0.08). The data are consistent with a transition state that is robust with regard to a wide range of statistically favored and disfavored beta-turn mutations and implicate a loosely assembled beta-hairpin as a key template in transition state stabilization with the beta-turn playing a central role.

Extending the folding nucleus of ubiquitin with an independently folding beta-hairpin finger: hurdles to rapid folding arising from the stabilisation of local interactions

The N-terminal beta-hairpin sequence of ubiquitin has been implicated as a folding nucleation site. To extend and stabilise the ubiquitin folding nucleus, we have inserted an autonomously folding 14-residue peptide sequence beta4 which in isolation forms a highly populated beta-hairpin (>70%) stabilised by local interactions. NMR structural analysis of the ubiquitin mutant (Ubeta4) shows that the hairpin finger is fully structured and stabilises ubiquitin by approximately 8kJmol(-1). Protein engineering and kinetic (phi(F)-value) analysis of a series of Ubeta4 mutants shows that the hairpin extension of Ubeta4 is also significantly populated in the transition state (phi(F)-values >0.7) and has the effect of templating the formation of native contacts in the folding nucleus of ubiquitin. However, at low denaturant concentrations the chevron plot of Ubeta4 shows a small deviation from linearity (roll-over effect), indicative of the population of a compact collapsed state, which appears to arise from over-stabilisation of local interactions. Destabilising mutations within the native hairpin sequence and within the engineered hairpin extension, but not elsewhere, eliminate this non-linearity and restore apparent two-state behaviour. The pitfall to stabilising local interactions is to present hurdles to the rapid and efficient folding of small proteins down a smooth folding funnel by trapping partially folded or misfolded states that must unfold or rearrange before refolding.

Non-linear effects of temperature and urea on the thermodynamics and kinetics of folding and unfolding of hisactophilin

Extensive measurements and analysis of thermodynamic stability and kinetics of urea-induced unfolding and folding of hisactophilin are reported for 5-50 degrees C, at pH 6.7. Under these conditions hisactophilin has moderate thermodynamic stability, and equilibrium and kinetic data are well fit by a two-state transition between the native and the denatured states. Equilibrium and kinetic m values decrease with increasing temperature, and decrease with increasing denaturant concentration. The betaF values at different temperatures and urea concentrations are quite constant, however, at about 0.7. This suggests that the transition state for hisactophilin unfolding is native-like and changes little with changing solution conditions, consistent with a narrow free energy profile for the transition state. The activation enthalpy and entropy of unfolding are unusually low for hisactophilin, as is also the case for the corresponding equilibrium parameters. Conventional Arrhenius and Eyring plots for both folding and unfolding are markedly non-linear, but these plots become linear for constant DeltaG/T contours. The Gibbs free energy changes for structural changes in hisactophilin have a non-linear denaturant dependence that is comparable to non-linearities observed for many other proteins. These non-linearities can be fit for many proteins using a variation of the Tanford model, incorporating empirical quadratic denaturant dependencies for Gibbs free energies of transfer of amino acid constituents from water to urea, and changes in fractional solvent accessible surface area of protein constituents based on the known protein structures. Noteworthy exceptions that are not well fit include amyloidogenic proteins and large proteins, which may form intermediates. The model is easily implemented and should be widely applicable to analysis of urea-induced structural transitions in proteins.

Domain-specific incorporation of noninvasive optical probes into recombinant proteins

An integrated approach is described that allows the domain-specific incorporation of optical probes into large recombinant proteins. The strategy is the combination of two existing techniques, expressed protein ligation (EPL) and in vivo amino acid replacement of tryptophans with tryptophan (Trp) analogues. The Src homology 3 (SH3) domain from the c-Crk-I adaptor protein has been labeled with a Trp analogue, 7-azatryptophan (7AW), using Escherichia coli Trp auxotrophs. Structural, biochemical, and thermodynamic studies show that incorporation of 7AW does not significantly perturb the structure or function of the isolated domain. Ligation of the 7AW-labeled SH3 domain to the c-Crk-I Src homology 2 (SH2) domain, via EPL, generated the multidomain protein, c-Crk-I, with a domain-specific label. Studies of this labeled protein show that the biochemical and thermodynamic properties of the SH3 domain do not change within the context of a larger multidomain protein. The technology described here is likely to be a useful tool in enhancing our understanding of the behavior of modular domains in their natural context, within multidomain proteins.

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.

Short amino acid stretches can mediate amyloid formation in globular proteins: the Src homology 3 (SH3) case

Protein misfolding and deposition underlie an increasing number of debilitating human disorders. We have shown that model proteins unrelated to disease, such as the Src homology 3 (SH3) domain of the p58alpha subunit of bovine phosphatidyl-inositol-3'-kinase (PI3-SH3), can be converted in vitro into assemblies with structural and cytotoxic properties similar to those of pathological aggregates. By contrast, homologous proteins, such as alpha-spectrin-SH3, lack the capability of forming amyloid fibrils at a measurable rate under any of the conditions we have so far examined. However, transplanting a small sequence stretch (6 aa) from PI3-SH3 to alpha-spectrin-SH3, comprising residues of the diverging turn and adjacent RT loop, creates an amyloidogenic protein closely similar in its behavior to the original PI3-SH3. Analysis of specific PI3-SH3 mutants further confirms the involvement of this region in conferring amyloidogenic properties to this domain. Moreover, the inclusion in this stretch of two consensus residues favored in SH3 sequences substantially inhibits aggregation. These findings show that short specific amino acid stretches can act as mediators or facilitators in the incorporation of globular proteins into amyloid structures, and they support the suggestion that natural protein sequences have evolved in part to code for structural characteristics other than those included in the native fold, such as avoidance of aggregation.

Stability and Folding Kinetics of a Ubiquitin Mutant with a Strong Propensity for Nonnative β-Hairpin Conformation in the Unfolded State

A F45W mutant of yeast ubiquitin has been used as a model system to examine the effects of nonnative local interactions on protein folding and stability. Mutating the native TLTGK G-bulged type I turn in the N-terminal beta-hairpin to NPDG stabilizes a nonnative beta-strand alignment in the isolated peptide fragment. However, NMR structural analysis of the native and mutant proteins shows that the NPDG mutant is forced to adopt the native beta-strand alignment and an unfavorable type I NPDG turn. The mutant is significantly less stable (approximately 9 kJ mol(-1)) and folds 30 times slower than the native sequence, demonstrating that local interactions can modulate protein stability and that attainment of a nativelike beta-hairpin conformation in the transition state ensemble is frustrated by the turn mutations. Surprising, alcoholic cosolvents [5-10% (v/v) TFE] are shown to accelerate the folding rate of the NPDG mutant. We conclude, backed-up by NMR data on the peptide fragments, that even though nonnative states in the denatured ensemble are highly populated and their stability further enhanced in the presence of cosolvents, the simultaneous increase in the proportion of nativelike secondary structure (hairpin or helix), in rapid equilibrium with nonnative states, is sufficient to accelerate the folding process. It is evident that modulating local interactions and increasing nonnative secondary structure propensities can change protein stability and folding kinetics. However, nonlocal contacts formed in the global cooperative folding event appear to determine structural specificity.

A thermodynamic and kinetic analysis of the folding pathway of an SH3 domain entropically stabilised by a redesigned hydrophobic core

The folding thermodynamics and kinetics of the alpha-spectrin SH3 domain with a redesigned hydrophobic core have been studied. The introduction of five replacements, A11V, V23L, M25V, V44I and V58L, resulted in an increase of 16% in the overall volume of the side-chains forming the hydrophobic core but caused no remarkable changes to the positions of the backbone atoms. Judging by the scanning calorimetry data, the increased stability of the folded structure of the new SH3-variant is caused by entropic factors, since the changes in heat capacity and enthalpy upon the unfolding of the wild-type and mutant proteins were identical at 298 K. It appears that the design process resulted in an increase in burying both the hydrophobic and hydrophilic surfaces, which resulted in a compensatory effect upon the changes in heat capacity and enthalpy. Kinetic analysis shows that both the folding and unfolding rate constants are higher for the new variant, suggesting that its transition state becomes more stable compared to the folded and unfolded states. The phi(double dagger-U) values found for a number of side-chains are slightly lower than those of the wild-type protein, indicating that although the transition state ensemble (TSE) did not change overall, it has moved towards a more denatured conformation, in accordance with Hammond's postulate. Thus, the acceleration of the folding-unfolding reactions is caused mainly by an improvement in the specific and/or non-specific hydrophobic interactions within the TSE rather than by changes in the contact order. Experimental evidence showing that the TSE changes globally according to its hydrophobic content suggests that hydrophobicity may modulate the kinetic behaviour and also the folding pathway of a protein.

Evidence for sequential barriers and obligatory intermediates in apparent two-state protein folding

Many small proteins fold fast and without detectable intermediates. This is frequently taken as evidence against the importance of partially folded states, which often transiently accumulate during folding of larger proteins. To get insight into the properties of free energy barriers in protein folding we analyzed experimental data from 23 proteins that were reported to show non-linear activation free-energy relationships. These non-linearities are generally interpreted in terms of broad transition barrier regions with a large number of energetically similar states. Our results argue against the presence of a single broad barrier region. They rather indicate that the non-linearities are caused by sequential folding pathways with consecutive distinct barriers and a few obligatory high-energy intermediates. In contrast to a broad barrier model the sequential model gives a consistent picture of the folding barriers for different variants of the same protein and when folding of a single protein is analyzed under different solvent conditions. The sequential model is also able to explain changes from linear to non-linear free energy relationships and from apparent two-state folding to folding through populated intermediates upon single point mutations or changes in the experimental conditions. These results suggest that the apparent discrepancy between two-state and multi-state folding originates in the relative stability of the intermediates, which argues for the importance of partially folded states in protein folding.

Mechanism of fast protein folding

An explosion of in vitro experimental data on the folding of proteins has revealed many examples of folding in the millisecond or faster timescale, often occurring in the absence of stable intermediate states. We review experimental methods for measuring fast protein folding kinetics, and then discuss various analytical models used to interpret these data. Finally, we classify general mechanisms that have been proposed to explain fast protein folding into two catagories, heterogeneous and homogeneous, reflecting the nature of the transition state. One heterogeneous mechanism, the diffusion-collision mechanism, can be used to interpret experimental data for a number of proteins.

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.

Unspecific hydrophobic stabilization of folding transition states

Here we present a method for determining the inference of non-native conformations in the folding of a small domain, alpha-spectrin Src homology 3 domain. This method relies on the preservation of all native interactions after Tyr/Phe exchanges in solvent-exposed, contact-free positions. Minor changes in solvent exposure and free energy of the denatured ensemble are in agreement with the reverse hydrophobic effect, as the Tyr/Phe mutations slightly change the polypeptide hydrophilic/hydrophobic balance. Interestingly, more important Gibbs energy variations are observed in the transition state ensemble (TSE). Considering the small changes induced by the H/OH replacements, the observed energy variations in the TSE are rather notable, but of a magnitude that would remain undetected under regular mutations that alter the folded structure free energy. Hydrophobic residues outside of the folding nucleus contribute to the stability of the TSE in an unspecific nonlinear manner, producing a significant acceleration of both unfolding and refolding rates, with little effect on stability. These results suggest that sectors of the protein transiently reside in non-native areas of the landscape during folding, with implications in the reading of phi values from protein engineering experiments. Contrary to previous proposals, the principle that emerges is that non-native contacts, or conformations, could be beneficial in evolution and design of some fast folding proteins.

Speeding protein folding beyond the G(o) model: how a little frustration sometimes helps

By perturbing a G(o) model toward a realistic protein Hamiltonian by adding non-native interactions, we find that the folding rate is in general enhanced as ruggedness is initially increased, as long as the protein is sufficiently large and flexible. Eventually, the rate drops rapidly toward zero when ruggedness significantly slows conformational transitions. Energy landscape arguments for thermodynamics and kinetics are coupled with a treatment of non-native collapse to elucidate this effect.

The pseudomolecule method and the structure of globular proteins. II. The example of ribonuclease F1 and T1

The pseudomolecule approach to the structure of globular proteins in which a small number of water molecules are incorporated into the "molecule" is tested again by comparing the ribbon of hydrogen bonds in two proteins, ribonuclease F1 and T1. These two molecules are 59% homologous and have the same backbone conformation both globally and locally. The two ribbons of hydrogen bonds that cover the whole of the backbone are conserved with an accuracy of some 95% providing that allowance is made for the intrusion into one of the pair of such extra factors as the presence of adducts or metal ions, the insertions and the absence of a few water molecules from one of the x-ray data sets. Without these corrections, the conservation of the ribbon is some 85%. There are 35 conserved hydrogen-bonding residues, nearly all of which show many unions to the backbone or interactions with the active site. There are 36 point mutations that involve one or two hydrogen-bonding side chains and nearly all of these have either none or one hydrogen bond to the backbone. These are minor contributors to the ribbon of hydrogen bonds. Of the 71 residues involved in these two categories, all but six fit into the pseudomolecular picture of the structure of globular proteins. The remaining 30 residues almost all contain conserved hydrocarbon side chains that may have a second order effect on the structure through their space filling effects.

Computer simulations aimed at structure prediction of supersecondary motifs in proteins

It is well established that protein structures are more conserved than protein sequences. One-third of all known protein structures can be classified into ten protein folds, which themselves are composed mainly of alpha-helical hairpin, beta hairpin, and betaalphabeta supersecondary structural elements. In this study, we explore the ability of a recent Monte Carlo-based procedure to generate the 3D structures of eight polypeptides that correspond to units of supersecondary structure and three-stranded antiparallel beta sheet. Starting from extended or misfolded compact conformations, all Monte Carlo simulations show significant success in predicting the native topology using a simplified chain representation and an energy model optimized on other structures. Preliminary results on model peptides from nucleotide binding proteins suggest that this simple protein folding model can help clarify the relation between sequence and topology.

Computer-assisted re-design of spectrin SH3 residue clusters

We have developed a protein design computer program, called Perla, which performs searches in sequence space to uncover optimal amino acid sequences for desired protein three-dimensional structures. Optimal sequences are localised at the minima of a sequence-structure energy landscape defined using a complex scoring function (an all-atom molecular mechanics force field plus statistical terms including entropy and solvation) measured with respect to a reference state simulating a denatured protein. Sequence choices eventually optimise side chain packing, secondary structure propensities, and hydrogen bonding and electrostatics interactions. Perla was used to re-design clusters of residues of the SH3 domain of alpha-spectrin. Several mutant proteins were produced and characterised. Some of our designed proteins have significantly higher stabilities (stability enhancements about 0.25, 0.70 and 1.0 kcal mol(-1)) than the wild-type protein. These successful protein re-designs, and similar examples found in the literature, establish the quality of the structure-based computational approach to protein design.

Bergerac-SH3: "frustation" induced by stabilizing the folding nucleus

The influence of an inserted exogenous independent folding element on the thermodynamics and folding properties of SH3 domain from alpha-spectrin has been investigated by creating a fused form between this small all-beta domain and a stable beta-hairpin (BH19). NMR analysis of synthetic peptides shows that insertion of BH19 nucleates formation of the original natural beta-hairpin (distal loop) that is part of the SH3 folding nucleus. The resulting protein (Bergerac-SHH) is more stable, folds faster and contains an elongated hairpin protruding from the globular domain as determined by 2D-NMR. "Protein engineering" analysis of the inserted region shows that it is folded in the transition state. Interestingly, stabilisation by insertion of the distal loop region results in the appearance of a compact intermediate revealed by a curved chevron plot at low denaturant concentration. This effect is eliminated at low salt concentrations by a single mutation of a hydrophobic residue within BH19 sequence, which is most probably involved in non-native interactions. Local stabilisation by enlargement and reinforcement of the folding nucleus, global compaction by the addition of salt and non-native interactions are shown to contribute to the observed deviation from the two-state behaviour.

Dramatic stabilization of an SH3 domain by a single substitution: roles of the folded and unfolded states

The N-terminal SH3 domain of the Drosophila drk protein (drkN SH3) exists in equilibrium between folded and unfolded states under non-denaturing buffer conditions. In order to examine the origins of this instability, we have made mutations in the domain and characterized the thermodynamics and kinetics of folding. Results of substitutions of negatively charged residues to neutral amino acid residues suggest that the large electrostatic potential of the domain does not play a dominant role in the instability of the domain. Sequence alignment of a large number of SH3 domains reveals that the drkN SH3 domain has a threonine (T22) at a position corresponding to an otherwise highly conserved glycine residue in the diverging beta-turn connecting the beta3 and beta4 strands. Mutation of T22 to glycine results in significant stabilization of the drkN SH3 domain by 2.5 kcal/mole. To further characterize the basis for the stabilization of the T22 mutant relative to wild-type, we made additional mutant proteins with substitutions of residue T22. A strong correlation is seen between protein stability or folding rate and propensity for native beta-turn structure at this position. Correlation of folding rates with AGADIR predictions of non-native helical structure in the diverging turn region, along with our previous NMR evidence for non-native structure in this region of the unfolded state of the drkN SH3 domain, suggests that the free energy of the unfolded state also plays a role in stability. This result highlights the importance of both folded and unfolded states for understanding protein stability.

Thermodynamic and structural characterization of Asn and Ala residues in the disallowed II' region of the Ramachandran plot

Residue Asn47 at position L1 of a type II' beta-turn of the alpha-spectrin SH3 domain is located in a disallowed region of the Ramachandran plot (phi = 56 +/- 12, psi = -118 +/- 17). Therefore, it is expected that replacement of Asn47 by Gly should result in a considerable stabilization of the protein. Thermodynamic analysis of the N47G and N47A mutants shows that the change in free energy is small (approximately 0.7 kcal/mol; approximately 3 kJ/mol) and comparable to that found when mutating a Gly to Ala in a alpha-helix or beta-sheet. X-ray structural analysis of these mutants shows that the conformation of the beta-turn does not change upon mutation and, therefore, that there is no relaxation of the structure, nor is there any gain or loss of interactions that could explain the small energy change. Our results indicate that the energetic definition of II' region of the Ramachandran plot (phi = 60 +/- 30, psi = -115 +/- 15) should be revised for at least Ala and Asn in structure validation and protein design.

The duplication of an eight-residue helical stretch in Staphylococcal nuclease is not helical: a model for evolutionary change

A common method of evolutionary change is gene duplication, followed by other events that lead to new function, decoration of folds, oligomerization, or other changes. As part of a study on the potential for evolutionary change created by duplicated sequences, we have carried out a crystallographic study on a mutant of Staphylococcal nuclease in which residues 55-62 have been duplicated in a wild-type variant termed PHS. In the parental protein (PHS) these residues form the first two turns of a helix running from residue 54 to 68 (hereafter designated as helix I). The crystal structure of the mutant is very similar to that of the parental, with helix I being unaltered. The duplicated residues are accommodated by expanding an existing loop N-terminal to helix I. In addition, circular dichroism (CD) studies have been carried out on a parental peptide containing helix I with six flanking residues at each terminus (residues 48-74) and on the same peptide expanded by the duplication, as a function of 2,2,2-trifluoroethanol (TFE) concentration. Each peptide possesses only modest helical propensity in solution. Our data, which is different from what was observed in T4 lysozyme, show that the conformation of the duplicated sequence is determined by a balance of sequential and longer-range effects. Thus duplicating sequence need not mean duplicating structure. Proteins 2000;40:465-472.

Circular dichroic investigation of the native and non-native conformational states of the growth factor receptor-binding protein 2 N-terminal src homology domain 3: effect of binding to a proline-rich peptide from guanine nucleotide exchange factor

SH3 (src homology domain 3) domains are small protein modules that interact with proline-rich peptides. The structure of the N-terminal SH3 domain from growth factor receptor-binding protein 2 (Grb2), an adapter protein in the intracellular signaling pathway to Ras, was investigated by circular dichroic (CD) spectroscopy. The compact native beta-barrel conformation, previously elucidated by NMR spectroscopy, was largely predominant at pH = 4.8, in the absence of salt. From the structural changes induced by varying pH, ionic strength, temperature, or hydrophobicity of the environment, evidence for the existence of distinct nonnative conformations was obtained in the far- and near-UV domains. Along the free energy scale, these appear to distribute into two conformational ensembles, depending on the extent of structural and thermodynamic differences compared to the native conformation. The first ensemble consists of non-native conformations with a nativelike secondary structure, and the second is composed of partially unfolded conformations having short alpha-helical fragments or turnlike motifs in their nonnative secondary structure. Most of the observed nonnative conformations exist in mild or nondenaturing conditions. They probably have distinct compactness of their inner structure, depending on the strength of nonlocal interactions, but only the native all-beta conformation possesses a condensed protein exterior, appropriate for the binding to the VPPPVPPRRR decapeptide from Sos. Upon binding, the native conformation undergoes a local tertiary structure change in a hydrophobic pocket at the binding site. This is accompanied by the PP-II helix folding of the proline-rich peptide. Interestingly, in the near-UV domain, a significant change in the spectral contribution of an aromatic exciton was observed, thus allowing quantitative tracking of the binding process.

Similarities between the spectrin SH3 domain denatured state and its folding transition state

We have expanded our description of the energy landscape for folding of the SH3 domain of chicken alpha-spectrin by a detailed structural characterization of its denatured state ensemble (DSE). This DSE is significantly populated under mildly acidic conditions in equilibrium with the folded state. Evidence from heteronuclear nuclear magnetic resonance (NMR) experiments on (2)H, (15)N-labeled protein suggests the presence of conformers whose residual structure bears some resemblence to the structure of the folding transition state of this protein. NMR analysis in a mutant with an engineered, non-native alpha-helical tendency shows a significant amount of local non-native structure in the mutant, while the overall characteristics of the DSE are unchanged. Comparison with recent theoretical predictions of SH3 domain folding reactions reveals an interesting correlation with the predicted early events. Based on these results and recent data from other systems, we propose that the DSE of a protein will resemble the intermediate or transition state of its nearest rate-limiting step, as a consequence of simple energetic and kinetic principles.

Is folding of beta-lactoglobulin non-hierarchic? Intermediate with native-like beta-sheet and non-native alpha-helix

The refolding of beta-lactoglobulin, a beta-barrel protein consisting of beta strands betaA-betaI and one major helix, is unusual because non-native alpha-helices are formed at the beginning of the process. We studied the refolding kinetics of bovine beta-lactoglobulin A at pH 3 using the stopped-flow circular dichroism and manual H/(2)H exchange pulse labeling coupled with heteronuclear NMR. The protection pattern from the H/(2)H exchange of the native state indicated the presence of a stable hydrophobic core consisting of betaF, betaG and betaH strands. The protection pattern of the kinetic intermediate obtained about one second after initiating the reaction was compared with that of the native state. In this relatively late kinetic intermediate, which still contains some non-native helical structure, the disulfide-bonded beta-hairpin made up of betaG and betaH strands was formed, but the rest of the molecule was fluctuating, where the non-native alpha-helices may reside. Subsequently, the core beta-sheet extends, accompanied by a further alpha-helix to beta-sheet transition. Thus, the refolding of beta-lactoglobulin exhibits two elements: the critical role of the core beta-sheet is consistent with the hierarchic mechanism, whereas the alpha-helix to beta-sheet transition suggests the non-hierarchic mechanism.

Exploring local and non-local interactions for protein stability by structural motif engineering

In order to probe the relative contribution of local and non-local interactions to the thermodynamic stability of proteins, we have devised an experimental approach based on a combination of motif engineering and sequence shuffling. Candidate chain segments in an immunoglobulin V(L) domain were identified whose conformation is proposed to be dominated by non-local interactions. Locally interacting structural motifs of a different conformation were then constructed as replacements, by introducing motif consensus sequences. We find that all nine replacements we constructed systematically reduce the folding cooperativity. By comparing this destabilising effect with the folding transitions of shuffled sequences for three of these motifs, we estimate the contribution of local, native interactions to the free energy of folding. Our results suggest that local and non-local interactions contribute to stability by an approximately equal amount, but that local interactions stabilise by increasing the resistance to denaturation while non-local interactions increase folding cooperativity. The systematic loss of stability by sequence shuffling in these host-guest experiments suggests that the designed interactions indeed are present in the native state, thus consensus sequence engineering may be a useful tool in structure design, but non-local interactions must be taken into account for global stability engineering. Statistical approaches are powerful tools for engineering protein structure and stability, but an analysis based on local sequence propensities alone does not adequately represent the balance of sequence and context in protein structures.

Synthesis, folding, and structure of the beta-turn mimic modified B1 domain of streptococcal protein G

The mechanism of beta-sheet formation remains a fundamental issue in our understanding of the protein folding process, but is hampered by the often encountered kinetic competition between folding and aggregation. The role of local versus nonlocal interactions has been probed traditionally by mutagenesis of both turn and strand residues. Recently, rigid organic molecules that impose a correct chain reversal have been introduced in several small peptides to isolate the importance of the long-range interactions. Here, we present the incorporation of a well-studied beta-turn mimic, designated as the dibenzofuran-based (DBF) amino acid, in the B1 domain of streptococcal protein G (B1G), and compare our results with those obtained upon insertion of the same mimic into the N-terminal beta-hairpin of B1G (O Melnyk et al., 1998, Lett Pept Sci 5:147-150). The DBF-B1G domain conserves the structure and the functional and thermodynamical properties of the native protein, whereas the modified peptide does not adopt a native-like conformation. The nature of the DBF flanking residues in the modified B1G domain prevents the beta-turn mimic from acting as a strong beta-sheet nucleator, which reinforces the idea that the native beta-hairpin formation is not driven by the beta-turn formation, but by tertiary interactions.

A tale of two secondary structure elements: when a beta-hairpin becomes an alpha-helix

In this work, we have analyzed the relative importance of secondary versus tertiary interactions in stabilizing and guiding protein folding. For this purpose, we have designed four different mutants to replace the alpha-helix of the GB1 domain by a sequence with strong beta-hairpin propensity in isolation. In particular, we have chosen the sequence of the second beta-hairpin of the GB1 domain, which populates the native conformation in aqueous solution to a significant extent. The resulting protein has roughly 30 % of its sequence duplicated and maintains the 3D-structure of the wild-type protein, but with lower stability (up to -5 kcal/mol). The loss of intrinsic helix stability accounts for about 80 % of the decrease in free energy, illustrating the importance of local interactions in protein stability. Interestingly enough, all the mutant proteins, included the one with the duplicated beta-hairpin sequence, fold with similar rates as the GB1 domain. Essentially, it is the nature of the rate-limiting step in the folding reaction that determines whether a particular interaction will speed up, or not, the folding rates. While local contacts are important in determining protein stability, residues involved in tertiary contacts in combination with the topology of the native fold, seem to be responsible for the specificity of protein structures. Proteins with non-native secondary structure tendencies can adopt stable folds and be as efficient in folding as those proteins with native-like propensities.

Instability, stabilization, and formulation of liquid protein pharmaceuticals

One of the most challenging tasks in the development of protein pharmaceuticals is to deal with physical and chemical instabilities of proteins. Protein instability is one of the major reasons why protein pharmaceuticals are administered traditionally through injection rather than taken orally like most small chemical drugs. Protein pharmaceuticals usually have to be stored under cold conditions or freeze-dried to achieve an acceptable shelf life. To understand and maximize the stability of protein pharmaceuticals or any other usable proteins such as catalytic enzymes, many studies have been conducted, especially in the past two decades. These studies have covered many areas such as protein folding and unfolding/denaturation, mechanisms of chemical and physical instabilities of proteins, and various means of stabilizing proteins in aqueous or solid state and under various processing conditions such as freeze-thawing and drying. This article reviews these investigations and achievements in recent years and discusses the basic behavior of proteins, their instabilities, and stabilization in aqueous state in relation to the development of liquid protein pharmaceuticals.

Structural characterization of an engineered tandem repeat contrasts the importance of context and sequence in protein folding

To test a different approach to understanding the relationship between the sequence of part of a protein and its conformation in the overall folded structure, the amino acid sequence corresponding to an alpha-helix of T4 lysozyme was duplicated in tandem. The presence of such a sequence repeat provides the protein with "choices" during folding. The mutant protein folds with almost wild-type stability, is active, and crystallizes in two different space groups, one isomorphous with wild type and the other with two molecules in the asymmetric unit. The fold of the mutant is essentially the same in all cases, showing that the inserted segment has a well-defined structure. More than half of the inserted residues are themselves helical and extend the helix present in the wild-type protein. Participation of additional duplicated residues in this helix would have required major disruption of the parent structure. The results clearly show that the residues within the duplicated sequence tend to maintain a helical conformation even though the packing interactions with the remainder of the protein are different from those of the original helix. It supports the hypothesis that the structures of individual alpha-helices are determined predominantly by the nature of the amino acids within the helix, rather than the structural environment provided by the rest of the protein.

Exploring the temperature-pressure phase diagram of staphylococcal nuclease

The temperature dependence of the pressure-induced equilibrium unfolding of staphylococcal nuclease (Snase) was determined by fluorescence of the single tryptophan residue, FTIR absorption for the amide I' and tyrosine O-H bands, and small-angle X-ray scattering (SAXS). The results from these three techniques were similar, although the stability as measured by fluorescence was slightly lower than that measured by FTIR and SAXS. The resulting phase diagram exhibits the well-known curvature for heat and cold denaturation of proteins, due to the large decrease in heat capacity upon folding. The volume change for unfolding became less negative with increasing temperatures, consistent with a larger thermal expansivity for the unfolded state than for the folded state. Fluorescence-detected pressure-jump kinetics measurements revealed that the curvature in the phase diagram is due primarily to the rate constant for folding, indicating a loss in heat capacity for the transition state relative to the unfolded state. The similar temperature dependence of the equilibrium and activation volume changes for folding indicates that the thermal expansivities of the folded and transition states are similar. This, along with the fact that the activation volume for folding is positive over the temperature range examined, the nonlinear dependence of the folding rate constant upon temperature implicates significant dehydration in the rate-limiting step for folding of Snase.

Thermodynamic analysis of alpha-spectrin SH3 and two of its circular permutants with different loop lengths: discerning the reasons for rapid folding in proteins

The temperature dependences of the unfolding-refolding reaction of a shorter version of the alpha-spectrin SH3 domain (PWT) used as a reference and of two circular permutants (with different poly-Gly loop lengths at the newly created fused loop) have been measured by differential scanning microcalorimetry and stopped-flow kinetics, to characterize the thermodynamic nature of the transition and native states. Differential scanning calorimetry results show that all these species do not belong to the same temperature dependency of heat effect. The family of the N47-D48s circular permutant (with 0-6 Gly inserted at the fused-loop) shows a higher enthalpy as happens with the PWT domain. The wild type (WT) and the S19-P20s permutant family have a more similar behavior although the second is far less stable. The crystallographic structure of the PWT shows a hairpin formation in the region corresponding to the unstructured N-terminus tail of the WT, explaining the enthalpic difference. There is a very good correlation between the calorimetric changes and the structural differences between the WT, PWT, and two circular permutants that suggests that their unfolded state cannot be too different. Elongation of the fused loop in the two permutants, taking as a reference the protein with one inserted Gly, results in a small Gibbs energy change of entropic origin as theoretically expected. Eyring plots of the unfolding and refolding semireactions show different behaviors for PWT, S19-P20s, and N47-D48s in agreement with previous studies indicating that they have different transition states. The SH3 transition state is relatively close to the native state with regard to changes in heat capacity and entropy, indicating a high degree of compactness and order. Regarding the differences in thermodynamic parameters, it seems that rapid folding could be achieved in proteins by decreasing the entropic barrier.

High populations of non-native structures in the denatured state are compatible with the formation of the native folded state

The structures of the denatured states of the spectrin SH3 domain and a mutant designed to have a non-native helical tendency at the N terminus have been analyzed under mild acidic denaturing conditions by nuclear magnetic resonance methods with improved resolution. The wild-type denatured state has little residual structure. However, the denatured state of the mutant has an approximately 50% populated helical structure from residues 2 to 14, a region that forms part of the beta-sheet structure in the folded state. Comparison with a peptide corresponding to the same sequence shows that the helix is stabilized in the whole domain, likely by non-local interactions with other parts of the protein as suggested by changes in a region far from the mutated sequence. These results demonstrate that high populations of non-native secondary structure elements in the denatured state are compatible with the formation of the native folded structure.

A test of the relationship between sequence and structure in proteins: excision of the heme binding site in apocytochrome b5

The water-soluble domain of rat hepatic holocytochrome b5 is an alphabeta protein containing elements of secondary structure in the sequence beta1-alpha1-beta4-beta3-alpha2-alpha3-beta5- alpha4-alpha5-beta2-alpha6. The heme group is enclosed by four helices, a2, a3, a4, and a5. To test the hypothesis that a small b hemoprotein can be constructed in two parts, one forming the heme site, the other an organizing scaffold, a protein fragment corresponding to beta1-alpha1-beta4-beta3-lambda-beta2-alpha6 was prepared, where lambda is a seven-residue linker bypassing the heme binding site. The fragment ("abridged b5") was found to contain alpha and beta secondary structure by circular dichroism spectroscopy and tertiary structure by Trp fluorescence emission spectroscopy. NMR data revealed a species with spectral properties similar to those of the full-length apoprotein. This folded form is in slow equilibrium on the chemical shift time scale with other less folded species. Thermal denaturation, as monitored by circular dichroism, absorption, and fluorescence spectroscopy, as well as size-exclusion chromatography-fast protein liquid chromatography (SEC-FPLC), confirmed the coexistence of at least two distinct conformational ensembles. It was concluded that the protein fragment is capable of adopting a specific fold likely related to that of cytochrome b5, but does not achieve high thermodynamic stability and cooperativity. Abridged b5 demonstrates that the spliced sequence contains the information necessary to fold the protein. It suggests that the dominating influence to restrict the conformational space searched by the chain is structural propensities at a local level rather than internal packing. The sequence also holds the properties necessary to generate a barrier to unfolding.

Obligatory steps in protein folding and the conformational diversity of the transition state

We have analyzed the existence of obligatory steps in the folding reaction of the alpha-spectrin SH3 domain by mutating Asp 48 (D48G), which is at position i+3 of an isolated two-residue type II' beta-turn. Calorimetry and X-ray analysis show an entropic stabilizing effect resulting from local changes at the dihedral angles of the beta-turn. Kinetic analysis of D48G shows that this beta-turn is fully formed in the transition state, while there is no evidence of its formation in an isolated fragment. Introduction of several mutations in the D48G protein reveals that the local stabilization has not significantly altered the transition state ensemble. All these results, together with previous analysis of other alpha-spectrin and src SH3 mutants, indicate that: (i) in the folding reaction there could be obligatory steps which are not necessarily part of the folding nucleus; (ii) transition state ensembles in beta-sheet proteins could be quite defined and conformationally restricted ('mechanic folding nucleus'); and (iii) transition state ensembles in some proteins could be evolutionarily conserved.

Important role of hydrogen bonds in the structurally polarized transition state for folding of the src SH3 domain

Experimental and theoretical studies on the folding of small proteins such as the chymotrypsin inhibitor 2 (CI-2) and the P22 Arc repressor suggest that the folding transition state is an expanded version of the native state with most interactions partially formed. Here we report that this picture does not hold generally: a hydrogen bond network involving two beta-turns and an adjacent hydrophobic cluster appear to be formed in the folding transition state of the src SH3 domain, while the remainder of the polypeptide chain is largely unstructured. Comparison with data on other small proteins suggests that this structural polarization is a consequence of the topology of the SH3 domain fold. The non-uniform distribution of structure in the folding transition state provides a challenging test for computational models of the folding process.

Design of a 20-amino acid, three-stranded beta-sheet protein

A 20-residue protein (named Betanova) forming a monomeric, three-stranded, antiparallel beta sheet was designed using a structural backbone template and an iterative hierarchical approach. Structural and physicochemical characterization show that the beta-sheet conformation is stabilized by specific tertiary interactions and that the protein exhibits a cooperative two-state folding-unfolding transition, which is a hallmark of natural proteins. The Betanova molecule constitutes a tractable model system to aid in the understanding of beta-sheet formation, including beta-sheet aggregation and amyloid fibril formation.

A conformational equilibrium in a protein fragment caused by two consecutive capping boxes: 1H-, 13C-NMR, and mutational analysis

The conformational properties of an 18 residues peptide spanning the entire sequence, L1KTPA5QFDAD10ELRAA15MKG, of the first helix (A-helix) of domain 2 of annexin I, were thoroughly investigated. This fragment exhibits several singular features, and in particular, two successive potential capping boxes, T3xxQ6 and D8xxE11. The former corresponds to the native hydrogen bond network stabilizing the alpha helix N-terminus in the protein; the latter is a non-native capping box able to break the helix at residue D8, and is observed in the domain 2 partially folded state. Using 2D-NMR techniques, we showed that two main populations of conformers coexist in aqueous solution. The first corresponds to a single helix extending from T3 to K17. The second corresponds to a broken helix at residue Ds. Four mutants, T3A, F7A, D8A, and E11A, were designed to further analyze the role of key amino acids in the equilibrium between the two ensembles of conformers. The sensitivity of NMR parameters to account for the variations in the populations of conformers was evaluated for each peptide. Our data show the delta13Calpha chemical shift to be the most relevant parameter. We used it to estimate the population ratio in the various peptides between the two main ensembles of conformers, the full helix and the broken helix. For the WT, E11A, and F7A peptides, these ratios are respectively 35/65, 60/40, 60/40. Our results were compared to the data obtained from helix/coil transition algorithms.

Characterisation of urea-denatured states of an immunoglobulin superfamily domain by heteronuclear NMR

The structural and dynamic properties of an immunoglobulin superfamily domain (IgSF), Ig 18', have been characterised by NMR at 285 K, in the presence of 4.2 M and 6.0 M urea, respectively. Analysis of chemical shift deviations, 3JHNHalpha coupling constants, sequential NOE pattern, and 15N relaxation data reveals that although the two urea-denatured states are highly disordered, some local turn-like residual structures do exist. Moreover, some distinct differences between the properties of the two denatured states are observed. In 4.2 M urea-denatured Ig 18', regions 80-83 and 86-92 adopt turn-like conformations, furthermore, region 84-93 is involved in slow exchange processes that occur on a micro- to millisecond time-scale. In the 6.0 M urea-denatured state, these turn-like conformations are less occupied, and chemical exchange processes in region 84-93 are largely reduced. In contrast, region 32-36 has persistent turn-like structures in both urea-denatured states. Some correlation between the spectral density function at 0 frequency, Jeff(0), for the urea-denatured states and the secondary structure elements of the folded state have been observed. Except for the terminal regions, residues corresponding to beta-strands have higher Jeff(0) values compared to residues corresponding to loops. The characterisation and comparison of the two urea-denatured states highlight residues that possess properties that may be crucial for the initiation of folding of this domain.

Conservation of rapid two-state folding in mesophilic, thermophilic and hyperthermophilic cold shock proteins

The cold shock protein CspB from Bacillus subtilis is only marginally stable, but it folds extremely fast in a simple N reversible U two-state reaction. The corresponding cold shock proteins from the thermophile Bacillus caldolyticus and the hyperthermophile Thermotoga maritima show strongly increased conformational stabilities, but unchanged very fast two-state refolding kinetics. The absence of intermediates in the folding of B. subtilis CspB is thus not a corollary of its low stability. Rather, two-state folding and an unusually native-like activated state of folding seem to be inherent properties of these small all-beta proteins. There is no link between stability and folding rate, and numerous sequence positions exist which can be varied to modulate the stability without affecting the rate and mechanism of folding.

Kinetic studies of beta-sheet protein folding

New studies have shown that folding of beta-sheet proteins can occur with and without intermediates, with fast to slow refolding rates and late to very late transition states. These experiments demonstrate that, despite early speculation to the contrary, beta-sheet protein folding does not appear to be fundamentally different from that of helical and mixed alpha, beta proteins.