Comparison of the structures and stabilities of coiled‐coil proteins containing hexafluoroleucine and t‐butylalanine provides insight into the stabilizing effects of highly fluorinated amino acid side‐chains

Highly fluorinated analogs of hydrophobic amino acids are well known to increase the stability of proteins toward thermal unfolding and chemical denaturation, but there is very little data on the structural consequences of fluorination. We have determined the structures and folding energies of three variants of a de novo designed 4‐helix bundle protein whose hydrophobic cores contain either hexafluoroleucine (hFLeu) or t‐butylalanine (tBAla). Although the buried hydrophobic surface area is the same for all three proteins, the incorporation of tBAla causes a rearrangement of the core packing, resulting in the formation of a destabilizing hydrophobic cavity at the center of the protein. In contrast, incorporation of hFLeu, causes no changes in core packing with respect to the structure of the nonfluorinated parent protein which contains only leucine in the core. These results support the idea that fluorinated residues are especially effective at stabilizing proteins because they closely mimic the shape of the natural residues they replace while increasing buried hydrophobic surface area.     

Hydrophobic effect on the stability and folding of a hyperthermophilic protein

Ribonuclease HII from hyperthermophile Thermococcus kodakaraensis (Tk-RNase HII) is a kinetically robust monomeric protein. The conformational stability and folding kinetics of Tk-RNase HII were measured for nine mutant proteins in which a buried larger hydrophobic side chain is replaced by a smaller one (Leu/Ile to Ala). The mutant proteins were destabilized by 8.9 to 22.0 kJ mol^− 1 as compared with the wild-type protein. The removal of each -CH₂- group burial decreased the stability by 5.1 kJ mol^− 1 on average in the mutant proteins of Tk-RNase HII examined. This is comparable with the value of 5.3 kJ mol^− 1 obtained from experiments for proteins from organisms growing at moderate temperature. We conclude that the hydrophobic residues buried inside protein molecules contribute to the stabilization of hyperthermophilic proteins to a similar extent as proteins at normal temperature. In the folding experiments, the mutant proteins of Tk-RNase HII examined exhibited faster unfolding compared with the wild-type protein. These results indicate that the buried hydrophobic residues strongly contribute to the kinetic robustness of Tk-RNase HII. This is the first report that provides a practical cause of slow unfolding of hyperthermostable proteins.     

Effects of the stabilization of the molten globule state on the folding mechanism of alpha-lactalbumin: a study of a chimera of bovine and human alpha-lactalbumin

The N-terminal half of the alpha-domain (residues 1 to 34) is more important for the stability of the acid-induced molten globule state of alpha-lactalbumin than the C-terminal half (residues 86 to 123). The refolding and unfolding kinetics of a chimera, in which the amino acid sequence of residues 1 to 34 was from human alpha-lactalbumin and the remainder of the sequence from bovine alpha-lactalbumin, were studied by stopped-flow tryptophan fluorescence spectroscopy. The chimeric protein refolded and unfolded substantially faster than bovine alpha-lactalbumin. The stability of the molten globule state formed by the chimera was greater than that of bovine alpha-lactalbumin, and the hydrophobic surface area buried inside of the molecule in the molten globule state was increased by the substitution of residues 1 to 34. Peptide fragments corresponding to the A- and B-helix of the chimera showed higher helix propensity than those of the bovine protein, indicating the contribution of local interactions to the high stability of the molten globule state of the chimera. Moreover, the substitution of residues 1-34 decreased the free energy level of the transition state and increased hydrophobic surface area buried inside of the molecule in the transition state. Our results indicate that local interactions as well as hydrophobic interactions formed in the molten globule state are important in guiding the subsequent structural formation of alpha-lactalbumin.     

Amino acid substitutions at protein–protein interfaces that modulate the oligomeric state

Despite similarities in their sequence and structure, there are a number of homologous proteins that adopt various oligomeric states. Comparisons of these homologous protein pairs, in terms of residue substitutions at the protein–protein interfaces, have provided fundamental characteristics that describe how proteins interact with each other. We have prepared a dataset composed of pairs of related proteins with different homo‐oligomeric states. Using the protein complexes, the interface residues were identified, and using structural alignments, the shadow‐interface residues have been defined as the surface residues that align with the interface residues. Subsequently, we investigated residue substitutions between the interfaces and the shadow interfaces. Based on the degree of the contributions to the interactions, the aligned sites of the interfaces and shadow interfaces were divided into primary and secondary sites; the primary sites are the focus of this work. The primary sites were further classified into two groups (i.e. exposed and buried) based on the degree to which the residue is buried within the shadow interfaces. Using these classifications, two simple mechanisms that mediate the oligomeric states were identified. In the primary‐exposed sites, the residues on the shadow interfaces are replaced by more hydrophobic or aromatic residues, which are physicochemically favored at protein–protein interfaces. In the primary‐buried sites, the residues on the shadow interfaces are replaced by larger residues that protrude into other proteins. These simple rules are satisfied in 23 out of 25 Structural Classification of Proteins (SCOP) families with a different‐oligomeric‐state pair, and thus represent a basic strategy for modulating protein associations and dissociations. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Sequence determinants of thermodynamic stability in a WW domain—An all‐β‐sheet protein

The stabilities of 66 sequence variants of the human Pin1 WW domain have been determined by equilibrium thermal denaturation experiments. All 34 residues composing the hPin1 WW three‐stranded β‐sheet structure could be replaced one at a time with at least one different natural or non‐natural amino acid residue without leading to an unfolded protein. Alanine substitutions at only four positions within the hPin1 WW domain lead to a partially or completely unfolded protein—in the absence of a physiological ligand. The side chains of these four residues form a conserved, partially solvent‐inaccessible, continuous hydrophobic minicore comprising the N‐ and C‐termini. Ala mutations at five other residues, three of which constitute the ligand binding patch on the concave side of the β‐sheet, significantly destabilize the hPin1 WW domain without leading to an unfolded protein. The remaining mutations affect protein stability only slightly, suggesting that only a small subset of side chain interactions within the hPin1 WW domain are mandatory for acquiring and maintaining a stable, cooperatively folded β‐sheet structure.     

Structural reorganization triggered by charging of Lys residues in the hydrophobic interior of a protein

Structural consequences of ionization of residues buried in the hydrophobic interior of proteins were examined systematically in 25 proteins with internal Lys residues. Crystal structures showed that the ionizable groups are buried. NMR spectroscopy showed that in 2 of 25 cases studied, the ionization of an internal Lys unfolded the protein globally. In five cases, the internal charge triggered localized changes in structure and dynamics, and in three cases, it promoted partial or local unfolding. Remarkably, in 15 proteins, the ionization of the internal Lys had no detectable structural consequences. Highly stable proteins appear to be inherently capable of withstanding the presence of charge in their hydrophobic interior, without the need for specialized structural adaptations. The extent of structural reorganization paralleled loosely with global thermodynamic stability, suggesting that structure-based pK(a) calculations for buried residues could be improved by calculation of thermodynamic stability and by enhanced conformational sampling.     

Study of the stability and unfolding mechanism of BBA1 by molecular dynamics simulations at different temperatures.

BBA1 is a designed protein that has only 23 residues. It is the smallest protein without disulfide bridges that has a well-defined tertiary structure in solution. We have performed unfolding molecular dynamics simulations on BBA1 and some of its mutants at 300, 330, 360, and 400 K to study their kinetic stability as well as the unfolding mechanism of BBA1. It was shown that the unfolding simulations can provide insights into the forces that stabilize the protein. Packing, hydrophobic interactions, and a salt bridge between Asp12 and Lys16 were found to be important to the protein's stability. The unfolding of BBA1 goes through two major steps: (1) disruption of the hydrophobic core and (2) unfolding of the helix. The beta-hairpin remains stable in the unfolding because of the high stability of the type II' turn connecting the two beta-strands.     

Boosting protein stability with the computational design of β‐sheet surfaces

β‐sheets often have one face packed against the core of the protein and the other facing solvent. Mutational studies have indicated that the solvent‐facing residues can contribute significantly to protein stability, and that the preferred amino acid at each sequence position is dependent on the precise structure of the protein backbone and the identity of the neighboring amino acids. This suggests that the most advantageous methods for designing β‐sheet surfaces will be approaches that take into account the multiple energetic factors at play including side chain rotamer preferences, van der Waals forces, electrostatics, and desolvation effects. Here, we show that the protein design software Rosetta, which models these energetic factors, can be used to dramatically increase protein stability by optimizing interactions on the surfaces of small β‐sheet proteins. Two design variants of the β‐sandwich protein from tenascin were made with 7 and 14 mutations respectively on its β‐sheet surfaces. These changes raised the thermal midpoint for unfolding from 45°C to 64°C and 74°C. Additionally, we tested an empirical approach based on increasing the number of potential salt bridges on the surfaces of the β‐sheets. This was not a robust strategy for increasing stability, as three of the four variants tested were unfolded.     

Role of conserved residues in structure and stability: tryptophans of human serum retinol-binding protein, a model for the lipocalin superfamily

Serum retinol binding protein (RBP) is a member of the lipocalin family, proteins with up-and-down beta-barrel folds, low levels of sequence identity, and diverse functions. Although tryptophan 24 of RBP is highly conserved among lipocalins, it does not play a direct role in activity. To determine if Trp24 and other conserved residues have roles in stability and/or folding, we investigated the effects of conservative substitutions for the four tryptophans and some adjacent residues on the structure, stability, and spectroscopic properties of apo-RBP. Crystal structures of recombinant human apo-RBP and of a mutant with substitutions for tryptophans 67 and 91 at 1.7 A and 2.0 A resolution, respectively, as well as stability measurements, indicate that these relatively exposed tryptophans have little influence on structure or stability. Although Trp105 is largely buried in the wall of the beta-barrel, it can be replaced with minor effects on stability to thermal and chemical unfolding. In contrast, substitutions of three different amino acids for Trp24 or replacement of Arg139, a conserved residue that interacts with Trp24, lead to similar large losses in stability and lower yields of native protein generated by in vitro folding. The results and the coordinated nature of natural substitutions at these sites support the idea that conserved residues in functionally divergent homologs have roles in stabilizing the native relative to misfolded structures. They also establish conditions for studies of the kinetics of folding and unfolding by identifying spectroscopic signals for monitoring the formation of different substructures.     

High polar content of long buried blocks of sequence in protein domains suggests selection against amyloidogenic non-polar sequences

Native protein structures achieve stability in part by burying hydrophobic side-chains. About 75% of all amino acid residues buried in protein interiors are non-polar. Buried residues are not uniformly distributed in protein sequences, but sometimes cluster as contiguous polypeptide stretches that run through the interior of protein domain structures. Such regions have an intrinsically high local sequence density of non-polar residues, creating a potential problem: local non-polar sequences also promote protein misfolding and aggregation into non-native structures such as the amyloid fibrils in Alzheimer's disease. Here we show that long buried blocks of sequence in protein domains of known structure have, on average, a lower content of non-polar amino acids (about 70%) than do isolated buried residues (about 80%). This trend is observed both in small and in large protein domains and is independent of secondary structure. Long, completely non-polar buried stretches containing many large side-chains are particularly avoided. Aspartate residues that are incorporated in long buried stretches were found to make fewer polar interactions than those in short stretches, hinting that they may be destabilizing to the native state. We suggest that evolutionary pressure is acting on non-native properties, causing buried polar residues to be placed at positions where they would break up aggregation-prone non-polar sequences, perhaps even at some cost to native state stability.     

Transient sampling of aggregation‐prone conformations causes pathogenic instability of a parkinsonian mutant of DJ‐1 at physiological temperature

Various missense mutations in the cytoprotective protein DJ‐1 cause rare forms of inherited parkinsonism. One mutation, M26I, diminishes DJ‐1 protein levels in the cell but does not result in large changes in the three‐dimensional structure or thermal stability of the protein. Therefore, the molecular defect that results in loss of M26I DJ‐1 protective function is unclear. Using NMR spectroscopy near physiological temperature, we found that the picosecond–nanosecond dynamics of wild‐type and M26I DJ‐1 are similar. In contrast, elevated amide hydrogen/deuterium exchange rates indicate that M26I DJ‐1 is more flexible than the wild‐type protein on longer timescales and that hydrophobic regions of M26I DJ‐1 are transiently exposed to solvent. Tryptophan fluorescence spectroscopy and thiol crosslinking analyzed by mass spectrometry also demonstrate that M26I DJ‐1 samples conformations that differ from the wild‐type protein at 37°C. These transiently sampled conformations are unstable and cause M26I DJ‐1 to aggregate in vitro at physiological temperature but not at lower temperatures. M26I DJ‐1 aggregation is correlated with pathogenicity, as the structurally similar but non‐pathogenic M26L mutation does not aggregate at 37°C. The onset of dynamically driven M26I DJ‐1 instability at physiological temperature resolves conflicting literature reports about the behavior of this disease‐associated mutant and illustrates the pitfalls of characterizing proteins exclusively at room temperature or below, as key aspects of their behavior may not be apparent.     

Filling small, empty protein cavities: structural and energetic consequences

Most proteins contain small cavities that can be filled by replacing cavity-lining residues by larger ones. Since shortening mutations in hydrophobic cores tend to destabilize proteins, it is expected that cavity-filling mutations may conversely increase protein stability. We have filled three small cavities in apoflavodoxin and determined by NMR and equilibrium unfolding analysis their impact in protein structure and stability. The smallest cavity (14 A3) has been filled, at two different positions, with a variety of residues and, in all cases, the mutant proteins are locally unfolded, their structure and energetics resembling those of an equilibrium intermediate of the thermal unfolding of the wild-type protein. In contrast, two slightly larger cavities of 20 A3 and 21 A3 have been filled with Val to Ile or Val to Leu mutations and the mutants preserve both the native fold and the equilibrium unfolding mechanism. From the known relationship, observed in shortening mutations, between stability changes and the differential hydrophobicity of the exchanged residues and the volume of the cavities, the filling of these apoflavodoxin cavities is expected to stabilize the protein by approximately 1.5 kcal mol(-1). However, both urea and thermal denaturation analysis reveal much more modest stabilizations, ranging from 0.0 kcal mol(-1) to 0.6 kcal mol(-1), which reflects that the accommodation of single extra methyl groups in small cavities requires some rearrangement, necessarily destabilizing, that lowers the expected theoretical stabilization. As the size of these cavities is representative of that of the typical small, empty cavities found in most proteins, it seems unlikely that filling this type of cavities will give rise to large stabilizations.     

Polar residues in the protein core of Escherichia coli thioredoxin are important for fold specificity

Most globular proteins contain a core of hydrophobic residues that are inaccessible to solvent in the folded state. In general, polar residues in the core are thermodynamically unfavorable except when they are able to form intramolecular hydrogen bonds. Compared to hydrophobic interactions, polar interactions are more directional in character and may aid in fold specificity. In a survey of 263 globular protein structures, we found a strong positive correlation between the number of polar residues at core positions and protein size. To probe the importance of buried polar residues, we experimentally tested the effects of hydrophobic mutations at the five polar core residues in Escherichia coli thioredoxin. Proteins with single hydrophobic mutations (D26I, C32A, C35A, T66L, and T77V) all have cooperative unfolding transitions like the wild type (wt), as determined by chemical denaturation. Relative to wt, D26I is more stable while the other point mutants are less stable. The combined 5-fold mutant protein (IAALV) is less stable than wt and has an unfolding transition that is substantially less cooperative than that of wt. NMR spectra as well as amide deuterium exchange indicate that IAALV is likely sampling a number of low-energy structures in the folded state, suggesting that polar residues in the core are important for specifying a well-folded native structure.     

Temperature-dependent downhill unfolding of ubiquitin. II. Modeling the free energy surface

To provide evidence for the interpretation of temperature‐dependent unfolding kinetics and the downhill unfolding scenario presented in the accompanying experimental article (Part I), the free energy surface of ubiquitin unfolding is calculated using statistical mechanical models of the Muñoz‐Eaton (ME) form. The models allow only two states for each amino acid residue, folded or unfolded, and permutations of these states generate an ensemble of microstates. One‐dimensional free energy curves are calculated using the number of folded residues as a reaction coordinate. The proposed sequential unfolding of ubiquitin's β‐sheet is tested by mapping the free energy onto two reaction coordinates inspired by the experiment as follows: the number of folded residues in ubiquitin's stable β‐strands I and II and those of the less stable strands III–V. Although the original ME model successfully captures folding features of zipper‐like one‐dimensional folders, it misses important tertiary interactions between residues that are far from each other in primary sequence. To take tertiary contacts into account, partially folded microstates based on a spherical growth model are included in the calculation and compared with the original model. By calculating the folding probability of each residue for a given point on the free energy surface, the unfolding pathway of ubiquitin is visualized. At low temperature, thermal unfolding occurs along a sequential unfolding pathway as follows: disruption of the β‐strands III–V followed by unfolding of the strands I and II. At high temperature, multiple unfolding routes are formed. The heterogeneity of the transition state explains the global nonexponential unfolding observed in the T‐jump experiment at high temperature. The calculation also reports a high stability for the α‐helix of ubiquitin. Proteins 2008. © 2008 Wiley‐Liss, Inc.     

Using physical features of protein core packing to distinguish real proteins from decoys

The ability to consistently distinguish real protein structures from computationally generated model decoys is not yet a solved problem. One route to distinguish real protein structures from decoys is to delineate the important physical features that specify a real protein. For example, it has long been appreciated that the hydrophobic cores of proteins contribute significantly to their stability. We used two sources to obtain datasets of decoys to compare with real protein structures: submissions to the biennial Critical Assessment of protein Structure Prediction competition, in which researchers attempt to predict the structure of a protein only knowing its amino acid sequence, and also decoys generated by 3DRobot, which have user‐specified global root‐mean‐squared deviations from experimentally determined structures. Our analysis revealed that both sets of decoys possess cores that do not recapitulate the key features that define real protein cores. In particular, the model structures appear more densely packed (because of energetically unfavorable atomic overlaps), contain too few residues in the core, and have improper distributions of hydrophobic residues throughout the structure. Based on these observations, we developed a feed‐forward neural network, which incorporates key physical features of protein cores, to predict how well a computational model recapitulates the real protein structure without knowledge of the structure of the target sequence. By identifying the important features of protein structure, our method is able to rank decoy structures with similar accuracy to that obtained by state‐of‐the‐art methods that incorporate many additional features. The small number of physical features makes our model interpretable, emphasizing the importance of protein packing and hydrophobicity in protein structure prediction.     

Understanding the folding and stability of a designed WW domain protein with replica exchange molecular dynamics simulations

WW domain proteins are usually regarded as simple models for understanding the folding mechanism of β-sheet. CC45 is an artificial protein that is capable of folding into the same structure as WW domain. In this article, the replica exchange molecular dynamics simulations are performed to investigate the folding mechanism of CC45. The analysis of thermal stability shows that β-hairpin 1 is more stable than β-hairpin 2 during the unfolding process. Free energy analysis shows that the unfolding of this protein substantially proceeds through solvating the smaller β-hairpin 2, followed by the unfolding of β-hairpin 1. We further propose the unfolding process of CC45 and the folding mechanism of two β-hairpins. These results are similar to the previous folding studies of formin binding protein 28 (FBP28). Compared with FBP28, it is found that CC45 has more aromatic residues in N-terminal loop, and these residues contact with C-terminal loop to form the outer hydrophobic core, which increases the stability of CC45. Knowledge about the stability and folding behaviour of CC45 may help in understanding the folding mechanisms of the β-sheet and in designing new WW domains.     

The formation of a native-like structure containing eight conserved hydrophobic residues is rate limiting in two-state protein folding of ACBP.

The acyl-coenzyme A-binding proteins (ACBPs) contain 26 highly conserved sequence positions. The majority of these have been mutated in the bovine protein, and their influence on the rate of two-state folding and unfolding has been measured. The results identify eight sequence positions, out of 24 probed, that are critical for fast productive folding. The residues are all hydrophobic and located in the interface between the N- and C-terminal helices. The results suggest that one specific site dominated by conserved hydrophobic residues forms the structure of the productive rate-determining folding step and that a sequential framework model can describe the protein folding reaction.     

Kinetic stability and crystal structure of the viral capsid protein SHP

SHP, the capsid-stabilizing protein of lambdoid phage 21, is highly resistant against denaturant-induced unfolding. We demonstrate that this high functional stability of SHP is due to a high kinetic stability with a half-life for unfolding of 25 days at zero denaturant, while the thermodynamic stability is not unusually high. Unfolding experiments demonstrated that the trimeric state (also observed in crystals and present on the phage capsid) of SHP is kinetically stable in solution, while the monomer intermediate unfolds very rapidly. We also determined the crystal structure of trimeric SHP at 1.5A resolution, which was compared to that of its functional homolog gpD. This explains how a tight network of H-bonds rigidifies crucial interpenetrating residues, leading to the observed extremely slow trimer dissociation or denaturation. Taken as a whole, our results provide molecular-level insights into natural strategies to achieve kinetic stability by taking advantage of protein oligomerization. Kinetic stability may be especially needed in phage capsids to allow survival in harsh environments.     

A new approach to the design of uniquely folded thermally stable proteins

A new computer program (CORE) is described that predicts core hydrophobic sequences of predetermined target protein structures. A novel scoring function is employed, which for the first time incorporates parameters directly correlated to free energies of unfolding (deltaGu), melting temperatures (Tm), and cooperativity. Metropolis-driven simulated annealing and low-temperature Monte Carlo sampling are used to optimize this score, generating sequences predicted to yield uniquely folded, stable proteins with cooperative unfolding transitions. The hydrophobic core residues of four natural proteins were predicted using CORE with the backbone structure and solvent exposed residues as input. In the two smaller proteins tested (Gbeta1, 11 core amino acids; 434 cro, 10 core amino acids), the native sequence was regenerated as well as the sequence of known thermally stable variants that exhibit cooperative denaturation transitions. Previously designed sequences of variants with lower thermal stability and weaker cooperativity were not predicted. In the two larger proteins tested (myoglobin, 32 core amino acids; methionine aminopeptidase, 63 core amino acids), sequences with corresponding side-chain conformations remarkably similar to that of native were predicted.     

Structural basis for unfolding pathway‐dependent stability of proteins: Vectorial unfolding versus global unfolding

Point mutations in proteins can have different effects on protein stability depending on the mechanism of unfolding. In the most interesting case of I27, the Ig‐like module of the muscle protein titin, one point mutation (Y9P) yields opposite effects on protein stability during denaturant‐induced “global unfolding” versus “vectorial unfolding” by mechanical pulling force or cellular unfolding systems. Here, we assessed the reason for the different effects of the Y9P mutation of I27 on the overall molecular stability and N‐terminal unraveling by NMR. We found that the Y9P mutation causes a conformational change that is transmitted through β‐sheet structures to reach the central hydrophobic core in the interior and alters its accessibility to bulk solvent, which leads to destabilization of the hydrophobic core. On the other hand, the Y9P mutation causes a bend in the backbone structure, which leads to the formation of a more stable N‐terminal structure probably through enhanced hydrophobic interactions.