The structure and function of omega loop A replacements in cytochrome c.

The structural and functional consequences of replacing omega-loop A (residues 18-32) in yeast iso-1-cytochrome c with the corresponding loop of Rhodospirillum rubrum cytochrome c2 have been examined. The three-dimensional structure of this loop replacement mutant RepA2 cytochrome c, and a second mutant RepA2(Val 20) cytochrome c in which residue 20 was back substituted to valine, were determined using X-ray diffraction techniques. A change in the molecular packing is evident in the RepA2 mutant protein, which has a phenylalanine at position 20, a residue considerably larger than the valine found in wild-type yeast iso-1-cytochrome c. The side chain of Phe 20 is redirected toward the molecular surface, altering the packing of this region of omega-loop A with the hydrophobic core of the protein. In the RepA2(Val 20) structure, omega-loop A contains a valine at position 20, which restores the original wild-type packing arrangement of the hydrophobic core. Also, as a result of omega-loop A replacement, residue 26 is changed from a histidine to asparagine, which results in displacements of the main-chain atoms near residue 44 to which residue 26 is hydrogen bonded. In vivo studies of the growth rate of the mutant strains on nonfermentable media indicate that the RepA2(Val 20) cytochrome c behaves much like the wild-type yeast iso-1 protein, whereas the stability and function of the RepA2 cytochrome c showed a temperature dependence. The midpoint reduction potential measured by cyclic voltammetry of the RepA2 mutant is 271 mV at 25 degrees C. This is 19 mV less than the wild-type and RepA2(Val 20) proteins (290 mV) and may result from disruption of the hydrophobic packing in the heme pocket and increased mobility of omega-loop A in RepA2 cytochrome c. The temperature dependence of the reduction potential is also greatly enhanced in the RepA2 protein.     

Minimization of cavity size ensures protein stability and folding: structures of Phe46-replaced bovine pancreatic RNase A

The Phe46 residue, located in the hydrophobic core of RNase A, was replaced with other hydrophobic residues, leucine, valine, or alanine, and their X-ray crystallographic structures were determined up to 1.50-1.80 A resolution in an attempt to examine the relationship between structural changes and conformational stability or folding kinetics. The backbone structure of F46L, F46V, and F46A was indistinguishable from that of the wild-type enzyme, retaining the correct active site structure. However, one water molecule was included in the hydrophobic core of F46A, forming two hydrogen bonds with the backbone peptide chain. The side chain of Met29 in F46V and F46A adopted two different conformations in an equal occupancy. A trapped water molecule and two conformations of Met29 represent changes that minimize the cavity volume. Nevertheless, the replacement of Phe46 with the above residues resulted in a marked decrease in both thermal stability and folding reaction. Thus, Phe46 ensures the thermal stability and the rapid and correct folding of RNase A by the role it plays in forming a highly packed, hydrophobic core.     

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.     

De novo design of the hydrophobic cores of proteins.

We have developed and experimentally tested a novel computational approach for the de novo design of hydrophobic cores. A pair of computer programs has been written, the first of which creates a "custom" rotamer library for potential hydrophobic residues, based on the backbone structure of the protein of interest. The second program uses a genetic algorithm to globally optimize for a low energy core sequence and structure, using the custom rotamer library as input. Success of the programs in predicting the sequences of native proteins indicates that they should be effective tools for protein design. Using these programs, we have designed and engineered several variants of the phage 434 cro protein, containing five, seven, or eight sequence changes in the hydrophobic core. As controls, we have produced a variant consisting of a randomly generated core with six sequence changes but equal volume relative to the native core and a variant with a "minimalist" core containing predominantly leucine residues. Two of the designs, including one with eight core sequence changes, have thermal stabilities comparable to the native protein, whereas the third design and the minimalist protein are significantly destabilized. The randomly designed control is completely unfolded under equivalent conditions. These results suggest that rational de novo design of hydrophobic cores is feasible, and stress the importance of specific packing interactions for the stability of proteins. A surprising aspect of the results is that all of the variants display highly cooperative thermal denaturation curves and reasonably dispersed NMR spectra. This suggests that the non-core residues of a protein play a significant role in determining the uniqueness of the folded structure.     

Native-specific stabilization of flavodoxin by the FMN cofactor: structural and thermodynamical explanation

Flavodoxins are useful models to investigate protein/cofactor interactions. The binding energy of the apoflavodoxin-FMN complex is high and therefore the holoflavodoxin is expected to be more stable than the apoprotein. This expectation has been challenged by reports on the stability of Desulfovibrio desulfuricans flavodoxin indicating that FMN binds to the unfolded polypeptide with similar affinity as to the native state, thus causing no net effect on protein stability. In previous work, we have analyzed in detail the stability of the apoflavodoxin from Anabaena PCC 7119 and the energetics of its functional complex with FMN. Here, we use the Anabaena holoprotein to directly investigate the contribution of the bound cofactor to protein stability through a detailed analysis of the chemical and thermal denaturation equilibria. Our data clearly shows that FMN binding largely stabilizes the protein towards both chemical and thermal denaturation, and that the stabilization observed at 25 degrees C in low ionic strength conditions is precisely the one expected if full release of the cofactor takes place upon flavodoxin unfolding. On the other hand, the binding of FMN to the native polypeptide is shown to simplify the thermal unfolding so that, while apoflavodoxin follows a three-state mechanism, the holoprotein unfolds in a two-state fashion. Comparison of the X-ray structure of native apoflavodoxin with the phi-structure of the thermal intermediate indicates that the increase in cooperativity driven by the cofactor originates in its preferential binding to the native state, which is a consequence of the disorganization in the intermediate of the FMN binding loops and of an adjacent longer loop.     

Conformational strictness required for maximum activity and stability of bovine pancreatic ribonuclease A as revealed by crystallographic study of three Phe120 mutants at 1.4 Å resolution

The replacement of Phe120 with other hydrophobic residues causes a decrease in the activity and thermal stability in ribonuclease A (RNase A). To explain this, the crystal structures of wild-type RNase A and three mutants--F120A, F120G, and F120W--were analyzed up to a 1.4 A resolution. Although the overall backbone structures of all mutant samples were nearly the same as that of wild-type RNase A, except for the C-terminal region of F120G with a high B-factor, two local conformational changes were observed at His119 in the mutants. First, His119 of the wild-type and F120W RNase A adopted an A position, whereas those of F120A and F120G adopted a B position, but the static crystallographic position did not reflect either the efficiency of transphosphorylation or the hydrolysis reaction. Second, His119 imidazole rings of all mutant enzymes were deviated from that of wild-type RNase A, and those of F120W and F120G appeared to be "inside out" compared with that of wild-type RNase A. Only approximately 1 A change in the distance between N(epsilon2) of His12 and N(delta1) of His119 causes a drastic decrease in k(cat), indicating that the active site requires the strict positioning of the catalytic residues. A good correlation between the change in total accessible surface area of the pockets on the surface of the mutant enzymes and enthalpy change in their thermal denaturation also indicates that the effects caused by the replacements are not localized but extend to remote regions of the protein molecule.     

Backbone-driven collapse in unfolded protein chains

Collapse of unfolded protein chains is an early event in folding. It affects structural properties of intrinsically disordered proteins, which take a considerable fraction of the human proteome. Collapse is generally believed to be driven by hydrophobic forces imposed by the presence of nonpolar amino acid side chains. Contributions from backbone hydrogen bonds to protein folding and stability, however, are controversial. To date, the experimental dissection of side-chain and backbone contributions has not yet been achieved because both types of interactions are integral parts of protein structure. Here, we realized this goal by applying mutagenesis and chemical modification on a set of disordered peptides and proteins. We measured the protein dimensions and kinetics of intra-chain diffusion of modified polypeptides at the level of individual molecules using fluorescence correlation spectroscopy, thereby avoiding artifacts commonly caused by aggregation of unfolded protein material in bulk. We found no contributions from side chains to collapse but, instead, identified backbone interactions as a source sufficient to form globules of native-like dimensions. The presence of backbone hydrogen bonds decreased polypeptide water solubility dramatically and accelerated the nanosecond kinetics of loop closure, in agreement with recent predictions from computer simulation. The presence of side chains, instead, slowed loop closure and modulated the dimensions of intrinsically disordered domains. It appeared that the transient formation of backbone interactions facilitates the diffusive search for productive conformations at the early stage of folding and within intrinsically disordered proteins.     

Hydrophobic Core Formation and Dehydration in Protein Folding Studied by Generalized-Ensemble Simulations

Despite its small size, chicken villin headpiece subdomain HP36 folds into the native structure with a stable hydrophobic core within several microseconds. How such a small protein keeps up its conformational stability and fast folding in solution is an important issue for understanding molecular mechanisms of protein folding. In this study, we performed multicanonical replica-exchange simulations of HP36 in explicit water, starting from a fully extended conformation. We observed at least five events of HP36 folding into nativelike conformations. The smallest backbone root mean-square deviation from the crystal structure was 1.1 Å. In the nativelike conformations, the stably formed hydrophobic core was fully dehydrated. Statistical analyses of the simulation trajectories show the following sequential events in folding of HP36: 1), Helix 3 is formed at the earliest stage; 2), the backbone and the side chains near the loop between Helices 2 and 3 take nativelike conformations; and 3), the side-chain packing at the hydrophobic core and the dehydration of the core side chains take place simultaneously at the later stage of folding. This sequence suggests that the initial folding nucleus is not necessarily the same as the hydrophobic core, consistent with a recent experimental ϕ-value analysis.     

Context-dependent effects of proline residues on the stability and folding pathway of ubiquitin.

Substitution of trans-proline at three positions in ubiquitin (residues 19, 37 and 38) produces significant context-dependent effects on protein stability (both stabilizing and destabilizing) that reflect changes to a combination of parameters including backbone flexibility, hydrophobic interactions, solvent accessibility to polar groups and intrinsic backbone conformational preferences. Kinetic analysis of the wild-type yeast protein reveals a predominant fast-folding phase which conforms to an apparent two-state folding model. Temperature-dependent studies of the refolding rate reveal thermodynamic details of the nature of the transition state for folding consistent with hydrophobic collapse providing the overall driving force. Br?nsted analysis of the refolding and unfolding rates of a family of mutants with a variety of side chain substitutions for P37 and P38 reveals that the two prolines, which are located in a surface loop adjacent to the C terminus of the main alpha-helix (residues 24-33), are not significantly structured in the transition state for folding and appear to be consolidated into the native structure only late in the folding process. We draw a similar conclusion regarding position 19 in the loop connecting the N-terminal beta-hairpin to the main alpha-helix. The proline residues of ubiquitin are passive spectators in the folding process, but influence protein stability in a variety of ways.     

Increasing the thermostability of staphylococcal nuclease: implications for the origin of protein thermostability

Seven hyper-stable multiple mutants have been constructed in staphylococcal nuclease by various combinations of eight different stabilizing single mutants. The stabilities of these multiple mutants determined by guanidine hydrochloride denaturation were 3.4 to 5.6 kcal/mol higher than that of the wild-type. Their thermal denaturation midpoint temperatures were 12.6 to 22.9 deg. C higher than that of the wild-type. These are among the greatest increases in protein stability and thermal denaturation midpoint temperature relative to the wild-type yet attained. There has been great interest in understanding how proteins found in thermophilic organisms are stabilized. One frequently cited theory is that the packing of hydrophobic side-chains is improved in the cores of proteins isolated from thermophiles when compared to proteins from mesophiles. The crystal structures of four single and five multiple stabilizing mutants of staphylococcal nuclease were solved to high resolution. No large overall structural change was found, with most changes localized around the sites of mutation. Rearrangements were observed in the packing of side-chains in the major hydrophobic core, although none of the mutations was in the core. It is surprising that detailed structural analysis showed that packing had improved, with the volume of the mutant protein's hydrophobic cores decreasing as protein stability increased. Further, the number of van der Waals interactions in the entire protein showed an experimentally significant increase correlated with increasing stability. These results indicate that optimization of packing follows as a natural consequence of increased protein thermostability and that good packing is not necessarily the proximate cause of high stability. Another popular theory is that thermostable proteins have more electrostatic and hydrogen bonding interactions and these are responsible for the high stabilities. The mutants here show that increased numbers of electrostatic and hydrogen bonding interactions are not obligatory for large increases in protein stability.     

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.     

Opposing structural changes in two symmetrical polypeptides bring about opposing changes to the thermal stability of a complex integral membrane protein

The relationship between membrane protein structure and thermal stability has been examined in the reaction centre from the bacterium Rhodobacter sphaeroides, a complex membrane protein comprising three polypeptide chains and 10 cofactors. The core of this protein exhibits an approximate twofold symmetry, the cofactors being held in two membrane-spanning branches by two polypeptides, termed L and M, that have very similar folds. In assays of the thermal stability of wild-type and mutant reaction centres embedded in the native bilayer membrane, replacement of a Phe at position 197 of the M polypeptide by His produced an increase in stability, whereas an opposing replacement of His by Phe at the symmetrical position 168 of the L-polypeptide produced a decrease in stability. In light of the known X-ray crystal structures of wild-type and mutant variants of this protein, and further mutagenesis, it is concluded that these stability changes result from the introduction or removal, respectively, of a hydrogen bond between the side-chain of the His and that of an Asn located two positions along the M or L polypeptide chain, in addition to a hydrogen bond between the His side-chain and an adjacent bacteriochlorophyll cofactor.     

Remodeling of the AB site of rat parvalbumin and oncomodulin into a canonical EF-hand.

Parvalbumin (PV) and the homologous protein oncomodulin (OM) contain three EF-hand motifs, but the first site (AB) cannot bind Ca2+. Here we aimed to recreate the putative ancestral proteins [D19-28E]PV and [D19-28E]OM by replacing the 10-residue-long nonfunctional loop in the AB site by a 12-residue canonical loop. To create an optical conformational probe we also expressed the homologs with a F102W replacement. Unexpectedly, in none of the proteins did the mutation reactivate the AB site. The AB-remodeled parvalbumins bind two Ca2+ ions with strong positive cooperativity (nH = 2) and moderate affinity ([Ca2+]0.5 = 2 microM), compared with [Ca2+]0.5 = 37 nM and nH = 1 for the wild-type protein. Increasing Mg2+ concentrations changed nH from 2 to 0.65, but without modification of the [Ca2+]0. 5-value. CD revealed that the Ca2+ and Mg2+ forms of the remodeled parvalbumins lost one-third of their alpha helix content compared with the Ca2+ form of wild-type parvalbumin. However, the microenvironment of single Trp residues in the hydrophobic cores, monitored using intrinsic fluorescence and difference optical density, is the same. The metal-free remodeled parvalbumins possess unfolded conformations. The AB-remodeled oncomodulins also bind two Ca2+ with [Ca2+]0.5 = 43 microM and nH = 1.45. Mg2+ does not affect Ca2+ binding. Again the Ca2+ forms display two-thirds of the alpha-helical content in the wild-type, while their core is still strongly hydrophobic as monitored by Trp and Tyr fluorescence. The metal-free oncomodulins are partially unfolded and seem not to possess a hydrophobic core. Our data indicate that AB-remodeled parvalbumin has the potential to regulate cell functions, whereas it is unlikely that [D19-28E]OM can play a regulatory role in vivo. The predicted evolution of the AB site from a canonical to an abortive EF-hand may have been dictated by the need for stronger interaction with Mg2+ and Ca2+, and a high conformational stability of the metal-free forms.     

Replacement of staphylococcal nuclease hydrophobic core residues with those from thermophilic homologues indicates packing is improved in some thermostable proteins

The importance of tight hydrophobic core packing in stabilizing proteins found in thermophilic organisms has been vigorously disputed. Here, portions of the cores found in three thermophilic homologues were transplanted into the core of staphylococcal nuclease, a protein of modest stability. Packing of the core was evaluated by comparing interaction energy of the three mutants to the comprehensive mutant library built up previously at these same sites in staphylococcal nuclease. It was found that the interaction energy of one thermophilic sequence is extraordinarily favorable and the interaction energies of other two transplanted thermophilic sequences are good, comparable to the interaction energies of mutant cores based on cores found in mesophilic homologues. As expected when transferring just a portion of the core sequence, the mutant proteins were destabilized overall relative to wild-type staphylococcal nuclease. The overall conclusion is that improvement of packing interactions is a mechanism to confer stability employed in some proteins from thermophiles, but not all.     

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.     

Akt/protein kinase B is regulated by autophosphorylation at the hypothetical PDK-2 site

The function of Akt (protein kinase B) is regulated by phosphorylation on two sites conserved within the AGC kinase family: the activation loop (Thr-308) in the kinase core and a hydrophobic phosphorylation site on the carboxyl terminus (Ser-473). Thr-308 is phosphorylated by the phosphoinositide-dependent kinase-1, (PDK-1), whereas the mechanism of phosphorylation of the hydrophobic site, tentatively referred to as the PDK-2 site, is unknown. Here we report that phosphorylation of the hydrophobic motif requires catalytically competent Akt. First we show that a kinase-inactive construct of Akt fails to incorporate phosphate at Ser-473 following IGF-1 stimulation in vivo but does incorporate phosphate at Thr-308 and a second carboxyl-terminal site, Thr-450; this ligand triggers the phosphorylation of both sites in wild-type enzyme. Neither does a catalytically inactive construct in which phosphorylation at the activation loop is blocked, T308A, become phosphorylated on the hydrophobic site in response to stimulation. Second, we show that Akt autophosphorylates on the hydrophobic site in vitro: phosphorylation of the activation loop by PDK-1 triggers the phosphorylation of the hydrophobic site in kinase-active, but not thermally inactivated, Akt alpha. Thus, Akt is regulated by autophosphorylation at the Ser-473 hydrophobic site.     

What makes a protein a protein? Hydrophobic core designs that specify stability and structural properties.

Here we describe how the systematic redesign of a protein's hydrophobic core alters its structure and stability. We have repacked the hydrophobic core of the four-helix-bundle protein, Rop, with altered packing patterns and various side chain shapes and sizes. Several designs reproduce the structure and native-like properties of the wild-type, while increasing the thermal stability. Other designs, either with similar sizes but different shapes, or with decreased sizes of the packing residues, destabilize the protein. Finally, overpacking the core with the larger side chains causes a loss of native-like structure. These results allow us to further define the roles of tight residue packing and the burial of hydrophobic surface area in the construction of native-like proteins.     

Increasing protein stability using a rational approach combining sequence homology and structural alignment: Stabilizing the WW domain

This study shows that a combination of sequence homology and structural information can be used to increase the stability of the WW domain by 2.5 kcal mol(-1) and increase the T(m) by 28 degrees C. Previous homology-based protein design efforts typically investigate positions with low sequence identity, whereas this study focuses on semi-conserved core residues and proximal residues, exploring their role(s) in mediating stabilizing interactions on the basis of structural considerations. The A20R and L30Y mutations allow increased hydrophobic interactions because of complimentary surfaces and an electrostatic interaction with a third residue adjacent to the ligand-binding hydrophobic cluster, increasing stability significantly beyond what additivity would predict for the single mutations. The D34T mutation situated in a pi-turn possibly disengages Asn31, allowing it to make up to three hydrogen bonds with the backbone in strand 1 and loop 2. The synergistic mutations A20R/L30Y in combination with the remotely located mutation D34T add together to create a hYap WW domain that is significantly more stable than any of the protein structures on which the design was based (Pin and FBP28 WW domains).     

Novel effects of FCCP [carbonyl cyanide p-(trifluoromethoxy)phenylhydrazone] on amyloid precursor protein processing

Amyloidogenic processing of the beta-amyloid precursor protein (APP) has been implicated in the pathology of Alzheimer's disease. Because it has been suggested that catabolic processing of the APP holoprotein occurs in acidic intracellular compartments, we studied the effects of the protonophore carbonyl cyanide p-(trifluoromethoxy)phenylhydrazone (FCCP) and the H+-ATPase inhibitor bafilomycin A1 on APP catabolism in human embryonic kidney 293 cells expressing either wild-type or "Swedish" mutant APP. Unlike bafilomycin A1, which inhibits beta-amyloid production in cells expressing mutant but not wild-type APP, FCCP inhibited beta-amyloid production in both cell types. Moreover, the effects of FCCP were independent of alterations in total cellular APP levels or APP maturation, and the concentrations used did not alter either cellular ATP levels or cell viability. Bafilomycin A1, which had no effect on beta-amyloid production in wild-type cells, inhibited endocytosis of fluorescent transferrin, whereas concentrations of FCCP that inhibited beta-amyloid production in these cells had no effect on endosomal function. Thus, in wild-type-expressing cells it appears that the beta-amyloid peptide is not produced in the classically defined endosome. Although bafilomycin A1 decreased beta-amyloid release from cells expressing mutant APP but not wild-type APP, it altered lysosomal function in both cell types, suggesting that in normal cells beta-amyloid is not produced in the lysosome. Although inhibition of beta-amyloid production by bafilomycin A1 in mutant cells may occur via changes in endosomal/lysosomal pH, our data suggest that FCCP inhibits wild-type beta-amyloid production by acting on a bafilomycin A1-insensitive acidic compartment that is distinct from either the endosome or the lysosome.