Helix stability and hydrophobicity in the folding mechanism of the bacterial immunity protein Im9

Recent models suggest that the mechanism of protein folding is determined by the balance between the stability of secondary structural elements and the hydrophobicity of the sequence. Here we determine the role of these factors in the folding kinetics of Im9* by altering the secondary structure propensity or hydrophobicity of helices I, II or IV by the substitution of residues at solvent exposed sites. The folding kinetics of each variant were measured at pH 7.0 and 10 degrees C, under which conditions wild-type Im9* folds with two-state kinetics. We show that increasing the helicity of these sequences in regions known to be structured in the folding intermediate of Im7*, switches the folding of Im9* from a two- to three-state mechanism. By contrast, increasing the hydrophobicity of helices I or IV has no effect on the kinetic folding mechanism. Interestingly, however, increasing the hydrophobicity of solvent-exposed residues in helix II stabilizes the folding intermediate and the rate-limiting transition state, consistent with the view that this helix makes significant non-native interactions during folding. The results highlight the generic importance of intermediates in folding and show that such species can be populated by increasing helical propensity or by stabilizing inter-helix contacts through non-native interactions.     

Fluorescence resonance energy transfer analysis of the folding pathway of Engrailed Homeodomain

The Engrailed Homeodomain folds on the microsecond time scale via an intermediate that is experimentally well characterised using structural Engrailed-Homeodomain mimics. Here, we analysed directly the changes in distance between key residues during the kinetics of unfolding and at equilibrium using fluorescence resonance energy transfer (FRET). Trp was the donor and 5-(((acetylamino)ethyl)amino) naphthalene-1-sulphate, the acceptor, substituted in positions that caused little change in stability. Distances calculated for the native state were in good agreement with those derived from the NMR structure. The distances between the N- and C-termini of Helix I and of Helix III increased, then decreased and finally increased again with increasing GdmCl concentration on equilibrium denaturation. This behaviour implied that there was a folding intermediate on the folding pathway and that this intermediate was populated at low concentrations of GdmCl concentration ( approximately 1 M). We analysed the changes in distance during temperature-jump relaxation kinetics, using a qualitative and very conservative procedure that drew conclusions only when changes in fluorescence of mutants containing either the donor or the acceptor alone would not obscure the change in the FRET signal when both donor and acceptor were present. The distance changes obtained under equilibrium and kinetic measurements were self-consistent and also consistent with the known high-resolution structures of the mimics of the folding intermediates. We showed that for analysing distances in disordered ensembles, it is important to use FRET probes with a critical distance close to the average separation in the ensemble. Otherwise, average distances could be over or underestimated.     

Kinetic folding mechanism of PDZ2 from PTP-BL

PDZ domains represent a large family of protein-interaction modules associated with a variety of unrelated proteins with different functions. We report a complete characterization of the kinetic folding mechanism of a fluorescent variant of PDZ2 from PTP-BL, investigated under a variety of different experimental conditions. For this purpose, we engineered a fluorescent variant of this protein Y43W (called pseudo-wild-type, pWT43). The results suggest the presence of a high-energy intermediate in the folding of PDZ2, as revealed by a pronounced non-linear dependence of the unfolding rate constant on denaturant concentration. Such an intermediate may or may not be detectable depending on the experimental conditions, giving rise to apparent two-state folding under stabilizing conditions (e.g. in the presence of sodium sulfate). Interestingly, even under these conditions, three-state folding can be restored by selectively destabilizing the native-like rate-limiting barrier by one specific mutation (V44A). Finally, we show that data taken on pWT43 under different experimental conditions (e.g. different pH values from 2.1 to 8.0 or in the presence of a stabilizing salt) and also data on a site-directed conservative mutant can be rationalized in terms of a simple reaction scheme involving a single set of intermediates and transition states.     

Folding and conformational studies on SCR1-3 domains of human complement receptor 1

Short consensus repeats SCR3 and SCR1-3 are soluble recombinantproteins, consisting of the third and first three N-terminaldomains of complement receptor 1, respectively, which retainsome anti-complement activity. The conformational stabilitiesand folding/unfolding of SCR3 and SCR1-3 have been studied usingcircular dichroism and equilibrium and pre-equilibrium fluorescencespectroscopy. Denaturation by guanidinium hydrochloride (GdnHCl)is rapid and completely reversible. Reduction of disulphidebridges in the folded proteins by ß-mercaptoethanolleads to an increase in fluorescence intensity. The fluorescenceintensity of the folded proteins is {small tilde}7.5% of thatof the respective unfolded proteins. The data can be approximatedto a two-state transition between native and denatured formsof the proteins. SCR3 has a conformational stability in waterof 12–13 kJ/mol whereas that of SCR1-3 is 19.5–19.9kJ/mol depending upon the technique utilized. The heat capacitychange associated with the unfolding of SCR1-3 was obtainedby a series of GdnHCl unfolding experiments over a range oftemperatures and was found to be 6.6 kJ/K.mol or 33.8 J/K.mol_residue.The refolding process of SCR3 was found to be simple, describedby a single exponential equation, whereas that of SCR1-3 wasfound to be complex and could be fitted to a double exponentialequation indicating the presence of folding intermediates.     

Mapping the folding pathway of the transmembrane protein DsbB by protein engineering

The four-helical transmembrane protein DsbB (disulfide bond reducing protein B) folds and unfolds reversibly in mixed anionic/non-ionic micelles, consisting of an unfolding intermediate I and a rate-limiting transition state (TS) between I and the denatured state D. Here, I describe the analysis of the folding behavior of 12 different alanine-scanning mutants of DsbB. For all mutants, TS is as compact as D and there is an accelerating increase in compaction as the protein proceeds to I and the native state. This unusual pattern of consolidation may reflect significant amounts of secondary structure in D, analogous to a classical folding intermediate. Unexpectedly, an increase in apolar surface area upon mutation is stabilizing whereas an increase in polar surface area is destabilizing. This effect is probably dominated by the effect of the mutations on the structure of the denatured state. I observe clear Hammond postulate behavior, in which a destabilization of I moves it closer to D. Φ-Value analysis indicates that in TS, a folding nucleus consisting of two to three residues with Φ-values of > 0.5 forms at one end of the transmembrane helices, which expands to include residues closer to the middle of the protein in I. Thus, folding proceeds from a highly polarized starting point.     

Conservation of mechanism, variation of rate: folding kinetics of three homologous four-helix bundle proteins

The amino acid sequence of a protein determines both its final folded structure and the folding mechanism by which this structure is attained. The differences in folding behaviour between homologous proteins provide direct insights into the factors that influence both thermodynamic and kinetic properties. Here, we present a comprehensive thermodynamic and kinetic analysis of three homologous homodimeric four-helix bundle proteins. Previous studies with one member of this family, Rop, revealed that both its folding and unfolding behaviour were interesting and unusual: Rop folds (k(0)(f) = 29 s(-1)) and unfolds (k(0)(u) = 6 x 10(-7) s(-1)) extremely slowly for a protein of its size that contains neither prolines nor disulphides in its folded structure. The homologues we discuss have significantly different stabilities and rates of folding and unfolding. However, the rate of protein folding directly correlates with stability for these homologous proteins: proteins with higher stability fold faster. Moreover, in spite of possessing differing thermodynamic and kinetic properties, the proteins all share a similar folding and unfolding mechanism. We discuss the properties of these naturally occurring Rop homologues in relation to previously characterized designed variants of Rop.     

Interactions between immunoglobulin-like and catalytic modules in Clostridium thermocellum cellulosomal cellobiohydrolase CbhA

Cellobiohydrolase CbhA from Clostridium thermocellum cellulosome is a multi-modular protein composed starting from the N-terminus of a carbohydrate-binding module (CBM) of family 4, an immunoglobulin(Ig)-like module, a catalytic module of family 9 glycoside hydrolases (GH9), X1(1) and X1(2) modules, a CBM of family 3 and a dockerin module. Deletion of the Ig-like module from the Ig-GH9 construct results in complete inactivation of the GH9 module. The crystal structure of the Ig-GH9 module pair reveals the existence of an extensive module interface composed of over 40 amino acid residues of both modules and maintained through a large number of hydrophilic and hydrophobic interactions. To investigate the importance of these interactions between the two modules, we compared the secondary and tertiary structures and thermostabilities of the individual Ig-like and GH9 modules and the Ig-GH9 module pair using both circular dichroism (CD) spectroscopy and differential scanning calorimetry (DSC). Thr230, Asp262 and Asp264 of the Ig-like module are located in the module interface of the Ig-GH9 module pair and are suggested to be important in 'communication' between the modules. These residues were mutated to alanyl residues. The structure, stability and catalytic properties of the native Ig-GH9 and its D264A and T230A/D262A mutants were compared. The results indicate that despite being able to fold relatively independently, the Ig-like and GH9 modules interact and these interactions affect the final fold and stability of each module. Mutations of one or two amino acid residues lead to destabilization and change of the mechanism of thermal unfolding of the polypeptides. The enzymatic properties of native Ig-GH9, D264A and T230A/D262A mutants are similar. The results indicate that inactivation of the GH9 module occurs as a result of multiple structural disturbances finally affecting the topology of the catalytic center.     

The importance of loop length in the folding of an immunoglobulin domain

Immunoglobulin (Ig)-like proteins have been shown to fold following formation of a nucleus comprising interactions between residues that are distant in the primary sequence. What role do the loops connecting these nucleus residues play? Here, the importance of loops connecting beta-strands in different sheets of the Ig fold is investigated, by insertion of five glycine residues into the B-C loop of an Ig domain from human titin, TI I27. The folding pathway of this elongated 'pseudo wild-type' TI I27 is probed using protein engineering and Phi-value analysis. The Phi-values calculated for mutants within the pseudo wild-type protein indicate that the folding nucleus in wild-type TI I27 is conserved, supporting the hypothesis that the inter-sheet loop is not critical to the formation of a long-range folding nucleus.     

Library screening studies to investigate substrate specificity in the reaction catalyzed by cholesterol oxidase

We tested whether it is possible to alter the substrate specificity of cholesterol oxidase for similarly sized sterols, i.e. cholesterol, beta-sitosterol and stigmasterol. Using existing X-ray crystal structures, we made a model of the predicted Michaelis complex of cholesterol and cholesterol oxidase. Based on this model, we identified five residues that are in direct contact with the steroid tail, Met58, Leu82, Val85, Met365 and Phe433. We prepared seven mutant libraries that contained the codon NYS (N = A, C, G, T; Y = C, Y; S = C, G) at one, two or three of the targeted positions by cassette mutagenesis. The libraries were screened for catalytic activity against three different sterols under k(cat)(*)/K(m)(*) conditions with 25 mol% sterol/DOPC unilamellar vesicles. The results of our screens suggest that specific packing interactions are not realized in the transition state of binding and that loss of active site water may be the predominant source of binding energy.     

The stability and unfolding of an IgG binding protein based upon the B domain of protein A from Staphylococcus aureus probed by tryptophan substitution and fluorescence spectroscopy

The stability and unfolding of an immunoglobulin (Ig) G bindingprotein based upon the B domain of protein A (SpAB) from Staphylococcusaureus were studied by substituting tryptophan residues at strategiclocations within each of the three a-helical regions (al-a3)of the domain. The role of the C-terminal helix, a3, was investigatedby generating two protein constructs, one corresponding to thecomplete SpAB, the other lacking a part of ct3; the Trp substitutionswere made in both one-and two-domain versions of each of theseconstructs. The fluorescence properties of each of the single-tryptophanmutants were studied in the native state and as a function ofguanidine-HCl-mediated unfolding, and their IgG binding activitieswere determined by a competitive enzyme-linked immunosorbentassay. The free energies of folding and of binding to IgG foreach mutant were compared with those for the native domains.The effect of each substitution upon the overall structure andupon the IgG binding interface was modelled by molecular graphicsand energy minimization. These studies indicate that (i) 3 contributesto the overall stability of the domain and to the formationof the IgG binding site in l and 2, and (ii) al unfolds first,followed by 2 and 3 together.     

Enhancing the thermal stability of mitochondrial cytochrome b5 by introducing a structural motif characteristic of the less stable microsomal isoform

Outer mitochondrial membrane cytochrome b5 (OM b5) is the most thermostable cytochrome b5 isoform presently known. Herein, we show that OM b5 thermal stability is substantially enhanced by swapping an apparently invariant motif in its heme-independent folding core with the corresponding motif characteristic of its less stable evolutionary relative, microsomal cytochrome b5 (Mc b5). The motif swap involved replacing two residues, Arg15 with His and Glu20 with Ser, thereby introducing a Glu11-His15-Ser20 H-bonding triad on the protein surface along with a His15/Trp22 pi-stacking interaction. The ferric and ferrous forms of the OM b5 R15H/E20S double mutant have thermal denaturation midpoints (Tm values) of approximately 93 degrees C and approximately 104 degrees C, respectively. A 15 degrees C increase in apoprotein Tm plays a key role in the holoprotein thermal stability enhancement, and is achieved by one of the most common natural mechanisms for stabilization of thermophilic versus mesophilic proteins: raising the unfolding free energy along the entire stability curve.     

Introduction of internal cysteines as conformational probes in yeast phosphoglycerate kinase

Several mutants of yeast phosphoglycerate kinase, each containingonly one internal cysteine residue, were constructed from asingle mutant devoid of cysteine. These cysteines were introducedas local conformational probes in selected buried positions.The enzyme activity, conformational characteristics and stabilityindicated that the mutations introduced only small perturbationsin the molecule. The folding–unfolding process mediatedby guanidine hydrochloride under equilibrium conditions wasstudied by following the variations in ellipticity and the reactivityof the cysteine residue towards 5,5'-dithiobis(nitrobenzoate).The process was found to be reversible except for mutant C97A,V49C,suggesting that this region located in helix I might be crucialin determining an intermediate on the folding pathway. The transitionsobtained by the two signals did not coincide, indicating thatthe local structures, in several parts inside the molecule,are more sensitive to the denaturant than the overall conformation.     

The Trp-cage: optimizing the stability of a globular miniprotein

The Trp-cage, as the smallest miniprotein, remains the subject of numerous computational and experimental studies of protein folding dynamics and pathways. The original Trp-cage (NLYIQWLKDGGPSSGRPPPS, Tm = 42 degrees C) can be significantly stabilized by mutations; melting points as high as 64 degrees C are reported. In helical portions of the structure, each allowed replacement of Leu, Ile, Lys or Ser residues by Ala results in a 1.5 (+/-0.35) kJ/mol fold stabilization. No changes in structure or fluxionality of the core results upon stabilization. Contrary to the initial hypothesis, specific Pro/Trp interactions are not essential for core formation. The entropic advantage of Pro versus Ala (DeltaDeltaS(U) = 11 +/- 2 J/mol K) was measured at the solvent-exposed P17 site. Pro-Ala mutations at two of the three prolines (P12 and P18) that encage the indole ring result in less fold destabilization (2.3-3.4 kJ/mol). However, a P19A mutation reduces fold stability by 16 kJ/mol reflecting a favorable Y3/P19 interaction as well as Trp burial. The Y3/P19 hydrophobic staple interaction defines the folding motif as an 18-residue unit. Other stabilizing features that have been identified include a solvent-exposed Arg/Asp salt bridge (3.4-6 kJ/mol) and a buried H-bonded Ser side chain ( approximately 10 kJ/mol).     

Protein stability for single substitution mutants and the extent of local compactness in the denatured state

The stability changes caused by single amino acid substitutionsare studied by a simple, empirical method which takes accountof the free energy change in the compact denatured state aswell as in the native state. The conformational free energyis estimated from effective inter-residue contact energies,as evaluated in our previous study. When this method is applied,with a simple assumption about the compactness of the denaturedstate, for single amino acid replacements at Glu49 of the tryptophansynthase subunit and at Ile3 of bacteriophage T4 lysozyme,the estimates of the unfolding Gibbs free energy changes correlatewell with observed values, especially for hydrophobic aminoacids, and it also yields the same magnitudes of energy as theobserved values for both proteins. When it is also applied foramino acid replacements at various positions to estimate theaverage number of contacts at each position in the denaturedstate from the observed value of unfolding free energy change,those values for replacements with Gly and Ala at the same residueposition in staphylococcal nuclease correlate well with eachother. The estimated numbers of contacts indicate that the proteinis not fully expanded in the denatured state and also that thecompact denatured state may have a substantially native-liketopology, like the molten globule state, in that there is aweak correlation between the estimated average number of contactsat each residue position in the denatured state and the numberof contacts in the native structure. These results provide somefurther evidence that the inter-residue contact energies asapplied here (i) properly reflect actual inter-residue interactionsand (ii) can be considered to be a pairwise hydrophobicity scale.Also, the results indicate that characterization of the denaturedstate is critical to understanding the folding process.     

Glycine 85 of the trp-repressor of E.coli is important in forming the hydrophobic tryptophan binding pocket: experimental and computational approaches

Experimental and computational analyses were performed on thecorepressor (L-tryptophan) binding site of the trp-repressorof Escherichia coli to investigate the ligandprotein interactions.Gly 85, one of the residues forming the hydrophobic pocket ofthe binding site, was systematically replaced with Ala, Val,Leu and Trp by cassette mutagenesis. Biochemical characterizationshowed that all these mutations caused significant decreasesin tryptophan binding activity. Free energy perturbation calculationswere performed for the mutants and were consistent with theexperimental results. The lack of a side chain at position 85was concluded to be essential for binding the corepressor; thestructure of the binding pocket was suggested to be tight inthe vicinity of Gly85.     

Cumulative stabilizing effects of glycine to alanine substitutions in Bacillus subtilis neutral protease

Oligonucleotide-directed mutagenesis has been used to replaceglycine residues by alanine in neutral protease from Bacillussubtilis. One Gly to Ala substitution (G147A) was located ina helical region of the protein, while the other (G189A) wasin a loop. The effects of mutational substitutions on the functional,conformational and stability properties of the enzyme have beeninvestigated using enzymatic assays and spectroscopic measurements.Single substitutions of both G1y147 and Gly189 with Ala residuesaffect the enzyme kinetic properties using synthetic peptidesas substrates. When Gly replacements were concurrently introducedat both positions, the kinetic characteristics of the doublemutant were roughly intermediate between those of the two singlemutants, and similar to those of the wild-type protease. Bothmutants G147A and G189A were found to be more stable towardsirreversible thermal inactivation/unfolding than the wild-typespecies. Moreover, the stabilizing effect of the Gly to Alasubstitution was roughly additive in the double mutant G147A/G189A,which shows a 3.2°C increase in T_m with respect to the wild-typeprotein. These findings indicate that the Gly to Ala substitutioncan be used as a strategy to stabilize globular proteins. Thepossible mechanisms of protein stabilization are also discussed.     

Use of a quantitative structure-property relationship to design larger model proteins that fold rapidly

A quantitative structure–property relationship (QSPR)was used to design model protein sequences that fold repeatedlyand relatively rapidly to stable target structures. The specificmodel was a 125-residue heteropolymer chain subject to MonteCarlo dynamics on a simple cubic lattice. The QSPR was derivedfrom an analysis of a database of 200 sequences by a statisticalmethod that uses a genetic algorithm to select the sequenceattributes that are most important for folding and a neuralnetwork to determine the corresponding functional dependenceof folding ability on the chosen attributes. The QSPR dependson the number of anti-parallel sheet contacts, the energy gapbetween the native state and quasi-continuous part of the spectrumand the total energy of the contacts between surface residues.Two Monte Carlo procedures were used in series to optimize boththe target structures and the sequences. We generated 20 fullyoptimized sequences and 60 partially optimized control sequencesand tested each for its ability to fold in dynamic MC simulations.Although sequences in which either the number of anti-parallelsheet contacts or the energy of the surface residues is non-optimalare capable of folding almost as well as fully optimized ones,sequences in which only the energy gap is optimized fold markedlymore slowly. Implications of the results for the design of proteinsare discussed.     

Folding core predictions from network models of proteins

Two different computational methods are employed to predict protein folding nuclei from native state structures, one based on an elastic network (EN) model and the other on a constraint network model of freely rotating rods. Three sets of folding cores are predicted with these models, and their correlation against the slow exchange folding cores identified by native state hydrogen-deuterium exchange (HX) experiments is used to test each method. These three folding core predictions rely on differences in the underlying models and relative importance of global or local motions for protein unfolding/folding reactions. For non-specific residue interactions, we use the Gaussian Network Model (GNM) to identify folding cores in the limits of two classes of motions, shortly referred to as global and local. The global mode minima from GNM represent the residues with the greatest potential for coordinating collective motions and are explored as potential folding nuclei. Additionally, the fast mode peaks that have previously been labeled as the kinetically hot residues are identified as a second folding core set dependent on local interactions. Finally, a third folding core set is defined by the most stable residues in a simulated thermal denaturation procedure of the FIRST software. This method uses an all-atomic analysis of the rigidity and flexibility of protein structures, which includes specific hydrophobic, polar and charged interactions. Comparison of the three folding core sets to HX data indicate that the fast mode peak residues determined by the GNM and the rigid folding cores of FIRST provide statistically significant enhancements over random correlation. The role of specific interactions in protein folding is also investigated by contrasting the differences between these two network-based computational methods.     

Stabilization of TRAIL, an all-beta-sheet multimeric protein, using computational redesign

Protein thermal stability is important for therapeutic proteins, both influencing the pharmacokinetic and pharmacodynamic properties and for stability during production and shelf-life of the final product. In this paper we show the redesign of a therapeutically interesting trimeric all-beta-sheet protein, the cytokine TRAIL (tumor necrosis factor-related apoptosis-inducing ligand), yielding variants with improved thermal stability. A combination of tumor necrosis factor (TNF) ligand family alignment information and the computational design algorithm PERLA was used to propose several mutants with improved thermal stability. The design was focused on non-conserved residues only, thus reducing the use of computational resources. Several of the proposed mutants showed a significant increase in thermal stability as experimentally monitored by far-UV circular dichroism thermal denaturation. Stabilization of the biologically active trimer was achieved by monomer subunit or monomer-monomer interface modifications. A double mutant showed an increase in apparent T(m) of 8 degrees C in comparison with wild-type TRAIL and remained biologically active after incubation at 73 degrees C for 1 h. To our knowledge, this is the first study that improves the stability of a large multimeric beta-sheet protein structure by computational redesign. A similar approach can be used to alter the characteristics of other multimeric proteins, including other TNF ligand family members.     

{beta}-sheet folding of fragment (16-36) of bovine pancreatic trypsin inhibitor as predicted by Monte Carlo simulated annealing

A tertiary structure prediction is described using Monte Carlosimulated annealing for the peptide fragment corresponding toresidues 16–36 of bovine pancreatic trypsin inhibitor(BPTI). The simulation starts with randomly chosen initial conformationsand is performed without imposing experimental constraints usingenergy functions given for generic interatomic interactions.Out of 20 simulation trials, seven conformations show a sheet-likestructure—two strands connected by a turn—althoughthis sheet-like structure is not as rigid as that observed innative BPTI. It is also shown that these conformations are mostlylooped and exhibit a native- like right-handed twist. Unlikethe case with the C-peptide of RNase A, no conspicuous -helicalstructure is found in any of the final conformations obtainedin the simulation. However, the lowest-energy conformation doesnot resemble exactly the native structure. This indicates thatthe rigid ß-sheet conformation of native BPTI merelycorresponds to a local minimum of the energy function if thefragment with residues 16–36 is isolated from the nativeprotein. A statistical analysis of all 20 final conformationssuggests that the tendency for the peptide segments to formextended ß-strands is strong for those with residues18–24, and moderate for those with residues 30–35.The segment of residues 25–29 does not tend to form anydefinite structure. In native BPTI, the former segments areinvolved in the ß-sheet and the latter in the turn.A folding scenario is also speculated from this analysis.