Contribution of amino acid substitutions at two different interior positions to the conformational stability of human lysozyme.

To elucidate correlative relationships between structural change and thermodynamic stability in proteins, a series of mutant human lysozymes modified at two buried positions (Ile56 and Ile59) were examined. Their thermodynamic parameters of denaturation and crystal structures were studied by calorimetry and X-ray crystallography. The mutants at positions 56 and 59 exhibited different responses to a series of amino acid substitutions. The changes in stability due to substitutions showed a linear correlation with changes in hydrophobicity of substituted residues, having different slopes at each mutation site. However, the stability of each mutant was found to be represented by a unique equation involving physical properties calculated from mutant structures. By fitting present and previous stability data for mutant human lysozymes substituted at various positions to the equation, the magnitudes of the hydrophobicity of a carbon atom and the hydrophobicity of nitrogen and neutral oxygen atoms were found to be 0.178 and -0.013 kJ/mol.A(2), respectively. It was also found that the contribution of a hydrogen bond with a length of 3.0 A to protein stability was 5.1 kJ/mol and the entropy loss of newly introduction of a water molecules was 7.8 kJ/mol.     

Role of surface hydrophobic residues in the conformational stability of human lysozyme at three different positions

To evaluate the contribution of the amino acid residues on the surface of a protein to its stability, a series of hydrophobic mutant human lysozymes (Val to Gly, Ala, Leu, Ile, Met, and Phe) modified at three different positions on the surface, which are located in the alpha-helix (Val 110), the beta-sheet (Val 2), and the loop (Val 74), were constructed. Their thermodynamic parameters of denaturation and crystal structures were examined by calorimetry and by X-ray crystallography at 100 K, respectively. Differences in the denaturation Gibbs energy change between the wild-type and the hydrophobic mutant proteins ranged from 4.6 to -9.6 kJ/mol, 2.7 to -1.5 kJ/mol, and 3.6 to -0.2 kJ/mol at positions 2, 74, and 110, respectively. The identical substitution at different positions and different substitutions at the same position resulted in different degrees of stabilization. Changes in the stability of the mutant proteins could be evaluated by a unique equation considering the conformational changes due to the substitutions [Funahashi et al. (1999) Protein Eng. 12, 841-850]. For this calculation, secondary structural propensities were newly considered. However, some mutant proteins were not adapted to the equation. The hydration structures around the mutation sites of the exceptional mutant proteins were affected due to the substitutions. The stability changes in the exceptional mutant proteins could be explained by the formation or destruction of the hydration structures. These results suggest that the hydration structure mediated via hydrogen bonds covering the protein surface plays an important role in the conformational stability of the protein.     

Positive contribution of hydration structure on the surface of human lysozyme to the conformational stability

Water molecules make a hydration structure with the network of hydrogen bonds, covering on the surface of proteins. To quantitatively estimate the contribution of the hydration structure to protein stability, a series of hydrophilic mutant human lysozymes (Val to Ser, Tyr, Asp, Asn, and Arg) modified at three different positions on the surface, which are located in the alpha-helix (Val-110), the beta-sheet (Val-2), and the loop (Val-74), were constructed. Their thermodynamic parameters of denaturation and crystal structures were examined by calorimetry and by x-ray crystallography at 100 K, respectively. The introduced polar residues made hydrogen bonds with protein atoms and/or water molecules, sometimes changing the hydration structure around the mutation site. Changes in the stability of the mutant proteins can be evaluated by a unique equation that considers the conformational changes resulting from the substitutions. Using this analysis, the relationship between the changes in the stabilities and the hydration structures for mutant human lysozymes substituted on the surface could be quantitatively estimated. The analysis indicated that the hydration structure on protein surface plays an important role in determining the conformational stability of the protein.     

Understanding the Origins of Loss of Protein Function by Analyzing the Effects of Thousands of Variants on Activity and Abundance

Understanding and predicting how amino acid substitutions affect proteins are keys to our basic understanding of protein function and evolution. Amino acid changes may affect protein function in a number of ways including direct perturbations of activity or indirect effects on protein folding and stability. We have analyzed 6,749 experimentally determined variant effects from multiplexed assays on abundance and activity in two proteins (NUDT15 and PTEN) to quantify these effects and find that a third of the variants cause loss of function, and about half of loss-of-function variants also have low cellular abundance. We analyze the structural and mechanistic origins of loss of function and use the experimental data to find residues important for enzymatic activity. We performed computational analyses of protein stability and evolutionary conservation and show how we may predict positions where variants cause loss of activity or abundance. In this way, our results link thermodynamic stability and evolutionary conservation to experimental studies of different properties of protein fitness landscapes.     

A single disulfide bond restores thermodynamic and proteolytic stability to an extensively mutated protein

The potential for engineering stable proteins with multiple amino acid substitutions was explored. Eleven lysine, five methionine, two tryptophan, one glycine, and three threonine substitutions were simultaneously made in barley chymotrypsin inhibitor-2 (CI-2) to substantially improve the essential amino acid content of the protein. These substitutions were chosen based on the three-dimensional structure of CI-2 and an alignment of homologous sequences. The initial engineered protein folded into a wild-type-like structure, but had a free energy of unfolding of only 2.2 kcal/mol, considerably less than the wild-type value of 7.5 kcal/mol. Restoration of the lysine mutation at position 67 to the wild-type arginine increased the free energy of unfolding to 3.1 kcal/mol. Subsequent cysteine substitutions at positions 22 and 82 resulted in disulfide bond formation and a protein with nearly wild-type thermodynamic stability (7.0 kcal/mol). None of the engineered proteins retained inhibitory activity against chymotrypsin or elastase, and all had substantially reduced inhibitory activity against subtilisin. The proteolytic stabilities of the proteins correlated with their thermodynamic stabilities. Reduction of the disulfide bond resulted in substantial loss of both thermodynamic and proteolytic stabilities, confirming that the disulfide bond, and not merely the cysteine substitutions, was responsible for the increased stability. We conclude that it is possible to replace over a third of the residues in CI-2 with minimal disruption of stability and structural integrity.     

Experimental verification of the 'stability profile of mutant protein' (SPMP) data using mutant human lysozymes.

The stability profile of mutant protein (SPMP) (Ota,M., Kanaya,S. and Nishikawa,K., 1995, J. Mol. Biol., 248, 733-738) estimates the changes in conformational stability due to single amino acid substitutions using a pseudo-energy potential developed for evaluating structure-sequence compatibility in the structure prediction method, the 3D-1D compatibility evaluation. Nine mutant human lysozymes expected to significantly increase in stability from SPMP were constructed, in order to experimentally verify the reliability of SPMP. The thermodynamic parameters for denaturation and crystal structures of these mutant proteins were determined. One mutant protein was stabilized as expected, compared with the wild-type protein. However, the others were not stabilized even though the structural changes were subtle, indicating that SPMP overestimates the increase in stability or underestimates negative effects due to substitution. The stability changes in the other mutant human lysozymes previously reported were also analyzed by SPMP. The correlation of the stability changes between the experiment and prediction depended on the types of substitution: there were some correlations for proline mutants and cavity-creating mutants, but no correlation for mutants related to side-chain hydrogen bonds. The present results may indicate some additional factors that should be considered in the calculation of SPMP, suggesting that SPMP can be refined further.     

Contribution of the 30/36 hydrophobic contact at the C-terminus of the alpha-helix to the stability of the ubiquitin molecule

The contribution of the hydrophobic contact in the C-capping motif of the alpha-helix to the thermodynamic stability of the ubiquitin molecule has been analyzed. For this, 16 variants of ubiquitin containing the full combinatorial set of four nonpolar residues Val, Ile, Leu, and Phe at C4 (Ile30) and C' ' (Ile36) positions were generated. The secondary structure content as estimated using far-UV circular dichroism (CD) spectroscopy of all but Phe variants at position 30 did not show notable changes upon substitutions. The thermodynamic stability of these ubiquitin variants was measured using differential scanning calorimetry, and it was shown that all variants have lower stability as measured by decreases in the Gibbs energy. Since in some cases the decrease in stability was so dramatic that it rendered an unfolded protein, it was therefore concluded that, despite apparent preservation of the secondary structure, the 30/36 hydrophobic contact is essential for the stability of the ubiquitin molecule. The decrease in the Gibbs energy in many cases was found to be accompanied by a large (up to 25%) decrease in the enthalpy of unfolding, particularly significant in the variants containing Ile to Leu substitutions. This decrease in enthalpy of unfolding is proposed to be primarily the result of the perturbed packing interactions in the native state of the Ile --> Leu variants. The analysis of these data and comparison with effects of similar amino acid substitutions on the stability of other model systems suggest that Ile --> Leu substitutions cannot be isoenergetic at the buried site.     

Approaches to predicting effects of single amino acid substitutions on the function of a protein

The relative activities of 313 mutants of the gene V protein of bacteriophage f1, assayed in vivo, have been used to evaluate two approaches to predicting the effects of single amino acid substitutions on the function of a protein. First, we tested methods that only depend on the properties of the wild-type and substituting amino acids. None of the properties or measures of the functional equivalence of amino acids we tested, including the frequency of exchange of amino acids among homologous proteins as well as changes in side-chain size, hydrophobicity, and charge, were found to be more than weakly correlated with the activities of mutants. The principal reason for this poor correlation was found to be that the effect of a particular substitution varies considerably from site to site. We then tested an approach using the activities of several mutants with substitutions at a site to predict the activity of another mutant, and we find that this is a relatively good indicator of whether the other mutant at that site will be functional. A predictive scheme was developed that combines the weak information from the models depending on the properties of the wild-type and substituting amino acids with the stronger information from the tolerance of a site to substitution. Although this scheme requires no knowledge of the structure of a mutant protein, it is useful in predicting the activities of mutants.     

Globular protein stability: aspects of interest in protein turnover

The conformational stability of globular proteins is remarkably low. Under physiological conditions, the native globular conformation is only from 5 to 15 kcal/mole more stable than unfolded conformations. In addition, small changes in the structure of a protein such as removing one terminal residue or cleaving a single peptide bond frequently lead to a substantial decrease in the stability. Likewise, single substitutions in the amino acid sequence can increase or decrease the stability by several kilocalories per mole. The low conformational stability of globular proteins and the sensitivity to small changes in structure suggest a possible role for conformational stability in the intracellular degradation of proteins. Several lines of evidence from in vivo studies of protein degradation are consistent with this idea.     

Contributions of aromatic pairs to the folding and stability of long-lived human γD-crystallin

Human γD-crystallin (HγD-Crys) is a highly stable protein that remains folded in the eye lens for the majority of an individual's lifetime. HγD-Crys exhibits two homologous crystallin domains, each containing two Greek key motifs with eight β-strands. Six aromatic pairs (four Tyr/Tyr, one Tyr/Phe and one Phe/Phe) are present in the β-hairpin sequences of the Greek keys. Ultraviolet damage to the aromatic residues in lens crystallins may contribute to the genesis of cataract. Mutant proteins with these aromatic residues substituted with alanines were constructed and expressed in E. coli. All mutant proteins except F115A and F117A had lower thermal stability than the WT protein. In equilibrium experiments in guanidine hydrochloride (GuHCl), all mutant proteins had lower thermodynamic stability than the WT protein. N-terminal domain (N-td) substitutions shifted the N-td transition to lower GuHCl concentration, but the C-terminal domain (C-td) transition remained unaffected. C-td substitutions led to a more cooperative unfolding/refolding process, with both the N-td and C-td transitions shifted to lower GuHCl concentration. The aromatic pairs conserved for each Greek key motif (Greek key pairs) had larger contributions to both thermal stability and thermodynamic stability than the other pairs. Aromatic-aromatic interaction was estimated as 1.5-2.0 kcal/mol. In kinetic experiments, N-td substitutions accelerated the early phase of unfolding, while C-td substitutions accelerated the late phase, suggesting independent domain unfolding. Only substitutions of the second Greek key pair of each crystallin domain slowed refolding. The second Greek keys may provide nucleation sites during the folding of the double-Greek-key crystallin domains.     

Stabilization of myoglobin by multiple alanine substitutions in helical positions.

We have carried out a series of multiple Xaa-->Ala changes at nonadjacent surface positions in the sequence of sperm whale myoglobin. Although the corresponding single substitutions do not increase the thermal stability of the protein, multiple substitutions enhance the stability of the resulting myoglobins. The effect observed is an increase in the observed Tm (midpoint unfolding temperature) relative to that predicted from assuming additivity of the free energy changes corresponding to single mutations. The stabilization occurs in the presence of urea, as measured by the dependence of the unfolding temperature on urea concentration. The sites that have been altered occur in different helices and are not close in sequence or in the native structure of myoglobin. The observed effect is consistent with a role of multiple alanines in residual interactions in the unfolded state of the mutant proteins.     

Structural and thermodynamic analysis of the packing of two alpha-helices in bacteriophage T4 lysozyme

Packing interactions in bacteriophage T4 lysozyme were explored by determining the structural and thermodynamic effects of substitutions for Ala98 and neighboring residues. Ala98 is buried in the core of T4 lysozyme in the interface between two alpha-helices. The Ala98 to Val (A98V) replacement is a temperature-sensitive lesion that lowers the denaturation temperature of the protein by 15 degrees C (pH 3.0, delta delta G = -4.9 kcal/mol) and causes atoms within the two helices to move apart by up to 0.7 A. Additional structural shifts also occur throughout the C-terminal domain. In an attempt to compensate for the A98V replacement, substitutions were made for Val149 and Thr152, which make contact with residue 98. Site-directed mutagenesis was used to construct the multiple mutants A98V/T152S, A98V/V149C/T152S and the control mutants T152S, V149C and A98V/V149I/T152S. These proteins were crystallized, and their high-resolution X-ray crystal structures were determined. None of the second-site substitutions completely alleviates the destabilization or the structural changes caused by A98V. The changes in stability caused by the different mutations are not additive, reflecting both direct interactions between the sites and structural differences among the mutants. As an example, when Thr152 in wild-type lysozyme is replaced with serine, the protein is destabilized by 2.6 kcal/mol. Except for a small movement of Val94 toward the cavity created by removal of the methyl group, the structure of the T152S mutant is very similar to wild-type T4 lysozyme. In contrast, the same Thr152 to Ser replacement in the A98V background causes almost no change in stability. Although the structure of A98V/T152S remains similar to A98V, the combination of T152S with A98V allows relaxation of some of the strain introduced by the Ala98 to Val replacement. These studies show that removal of methyl groups by mutation can be stabilizing (Val98----Ala), neutral (Thr152----Ser in A98V) or destabilizing (Val149----Cys, Thr152----Ser). Such diverse thermodynamic effects are not accounted for by changes in buried surface area or free energies of transfer of wild-type and mutant side-chains. In general, the changes in protein stability caused by a mutation depend not only on changes in the free energy of transfer associated with the substitution, but also on the structural context within which the mutation occurs and on the ability of the surrounding structure to relax in response to the substitution.(ABSTRACT TRUNCATED AT 400 WORDS)     

Disease-associated mutations in the coil 2B domain of human lamin A/C affect structural properties that mediate dimerization and intermediate filament formation

The lamin proteins are essential components of the nuclear lamina of eukaryotic cells, that are involved in a complex association mechanism to attain a functional supermolecular structure. Mutations of the lamin A/C gene are associated with several different neuromuscular diseases, and the detailed effect of disease-associated amino acid substitutions on the structure and stability of human lamin dimers is yet unknown. Here we present a structural and thermodynamic characterization by means of molecular dynamics simulations of the effect of pathological mutations (S326T, R331P, R331Q, E347K, E358K, M371K, and R377H) on the association of the coil 2B domains that mediate lamin A/C oligomerization. The structures attained during the simulations, along with the quantification of the contribution of each residue to the dimerization energies, support a lamin association mechanism mediated by homophilic intermolecular interactions promoted by dissociative conformational changes at distinct positions in the coiled coil. The pathogenic mutations can both increase or decrease the stability of lamin A/C dimers, and a possible correlation between the effect of the amino acid substitutions and disease onset and severity is presented.     

iPTREE-STAB: interpretable decision tree based method for predicting protein stability changes upon mutations

We have developed a web server, iPTREE-STAB for discriminating the stability of proteins (stabilizing or destabilizing) and predicting their stability changes (delta deltaG) upon single amino acid substitutions from amino acid sequence. The discrimination and prediction are mainly based on decision tree coupled with adaptive boosting algorithm, and classification and regression tree, respectively, using three neighboring residues of the mutant site along N- and C-terminals. Our method showed an accuracy of 82% for discriminating the stabilizing and destabilizing mutants, and a correlation of 0.70 for predicting protein stability changes upon mutations. AVAILABILITY: http://bioinformatics.myweb.hinet.net/iptree.htm. SUPPLEMENTARY INFORMATION: Dataset and other details are given.     

Design and characterization of a hyperstable p16INK4a that restores Cdk4 binding activity when combined with oncogenic mutations

Cyclin-dependent kinase inhibitor p16(INK4a) is the founding member of the INK4 family of tumor suppressors capable of arresting mammalian cell division. Missense mutations in the p16(INK4a) gene (INK4a/CDKN2A/MTS1) are strongly linked to several types of human cancer. These mutations are evenly distributed throughout this small, ankyrin repeat protein and the majority of them disrupt the native secondary and/or tertiary structure, leading to protein unfolding, aggregation and loss of function. We report here the use of multiple stabilizing substitutions to increase the stability of p16(INK4a) and furthermore, to restore Cdk4 binding activity of several defective, cancer-related mutant proteins. Stabilizing substitutions were predicted using four different techniques. The three most effective substitutions were combined to create a hyperstable p16(INK4a) variant that is 1.4 kcal/mol more stable than wild-type. This engineered construct is monomeric in solution with wild-type-like secondary and tertiary structure and cyclin-dependent kinase 4 binding activity. Interestingly, these hyperstable substitutions, when combined with oncogenic mutations R24P, P81L or V126D, can significantly restore Cdk4 binding activity, despite the divergent features of each destabilizing mutation. Extensive biophysical studies indicate that the hyperstable substitutions enhance the binding activity of mutant p16 through several different mechanisms, including an increased amount of secondary structure and thermostability, reduction in exposed hydrophobic surface(s) and/or a reduced tendency to aggregate. This apparent global suppressor effect suggests that increasing the thermodynamic stability of p16 can be used as a general strategy to restore the biological activity to defective mutants of this important tumor suppressor protein.     

Hydrophobic core packing in the SH3 domain folding transition state.

How tightly packed is the hydrophobic core of a folding transition state structure? We have addressed this question by characterizing the effects on folding kinetics of > 40 substitutions of both large and small amino acids in the hydrophobic core of the Fyn SH3 domain. Our results show that residues at three positions, which we designate as the 'core folding nucleus', are tightly packed in the transition state, and substitutions at these positions cause the largest changes in the folding rate. The other six positions examined appear to be loosely packed; thus, substitutions at these positions with larger hydrophobic residues generally accelerate folding, presumably by increasing the rate of nonspecific hydrophobic collapse. Surprisingly, the folding rate can be greatly accelerated by residues that also significantly destabilize the native state structure. Furthermore, mutants with identical thermodynamic stability can differ by up to 55-fold in their folding rates. These results highlight the importance of hydrophobic core composition, as opposed to only topology, in determining the folding rate of a protein. They also provide a new explanation for the 'abnormal' phi-values observed in many protein folding kinetics studies.     

The relationship between conservation, thermodynamic stability, and function in the SH3 domain hydrophobic core

To investigate the relationships between sequence conservation, protein stability, and protein function, we have measured the thermodynamic stability, folding kinetics, and in vitro peptide-binding activity of a large number of single-site substitutions in the hydrophobic core of the Fyn SH3 domain. Comparison of these data to that derived from an analysis of a large alignment of SH3 domain sequences revealed a very good correlation between the distinct pattern of conservation observed at each core position and the thermodynamic stability of mutants. Conservation was also found to correlate well with the unfolding rates of mutants, but not to the folding rates, suggesting that evolution selects more strongly for optimal native state packing interactions than for maximal folding rates. Structural analysis suggests that residue-residue core packing interactions are very similar in all SH3 domains, which provides an explanation for the correlation between conservation and mutant stability effects studied in a single SH3 domain. We also demonstrate a correlation between stability and the in vivo activity of mutants, and between conservation and activity. However, the relationship between conservation and activity was very strong only for the three most conserved hydrophobic core positions. The weaker correlation between activity and conservation seen at the other seven core positions indicates that maintenance of protein stability is the dominant selective pressure at these positions. In general, the pattern of conservation at hydrophobic core positions appears to arise from conserved packing constraints, and can be effectively utilized to predict the destabilizing effects of amino acid substitutions.     

Equilibrium unfolding of mutant apomyoglobins with substitutions of conserved nonfunctional residues by alanine

The problems of protein aggregation and protein misfolding in the cell are connected with the appearance of many genetic diseases. Both processes can be a consequence of substitutions of certain amino acid residues in proteins. The substitutions can influence the protein stability and protein folding rates in both the intermediate and the native states. We have studied equilibrium urea unfolding of mutant forms of apomyoglobin with substitutions of conserved nonfunctional residues by Ala to estimate their influence on protein stability. These residues include Val10, Trp14, Ilel11, Leu115, Met131 and Leu135. Conformational transitions were monitored by intrinsic Trp fluorescence and by circular dichroism spectra in the far UV region. Free energy changes upon the transition from the native to intermediate state and from the intermediate to unfolded state were determined. It was shown that all substitutions used lead to an appreciable decrease of the apomyoglobin native state stability, whereas the stability of the intermediate state is affected substantially smaller.     

Detection of conserved physico-chemical characteristics of proteins by analyzing clusters of positions with co-ordinated substitutions.

MOTIVATION: It is known that the physico-chemical characteristics of proteins underlying specific folding of the polypeptide chain and the protein function are evolutionary conserved. Detection of such characteristics while analyzing homologous sequences would expand essentially the knowledge on protein function, structure, and evolution. These characteristics are maintained constant, in particular, by co-ordinated substitutions. In this process, the destabilizing effect of a substitution may be compensated by another substitution at a different position within the same protein, making the overall change in this protein characteristic insignificant. Consequently, the patterns of co-ordinated substitutions contain important information on conserved physico-chemical properties of proteins, requiring their investigation and development of the corresponding methods and software for correlation analysis of protein sequences available to a wide range of users. RESULTS: A software package for analyzing correlated amino acid substitutions at different positions within aligned protein sequences was developed. The approach implies searching for evolutionary conserved physico-chemical characteristics of proteins based on the information on the pairwise correlations of amino acid substitutions at different protein positions. The software was applied to analyze DNA-binding domains of the homeodomain class. As a result, two conservative physico-chemical characteristics preserved due to the co-ordinated substitutions at certain groups of positions in the protein sequence. Possible functional roles of these characteristics are discussed. AVAILABILITY: The program package is available at http://wwwmgs.bionet.nsc.ru/programs/CRASP/.     

Binding interactions of kistrin with platelet glycoprotein IIb-IIIa: Analysis by site-directed mutagenesis

The binding interactions between platelet fibrinogen receptor, glycoprotein (GP) IIb-IIIa, and kistrin, a snake venom disintegrin protein that contains the adhesion site recognition sequence Arg-Gly-Asp (RGD) and potently inhibits platelet aggregation, have been investigated by site-directed mutagenesis of a synthetic kistrin gene. Kistrin was expressed as a fusion protein in Escherichia coli under control of the alkaline phosphatase promoter. This construction included the stII signal sequence to direct secretion to the periplasmic space and one synthetic (Z) domain of Staphylococcal protein A to allow affinity purification using IgG Sepharose. Kistrin was cleaved from the Z-domain by site-specific proteolysis using a mutant subtilisin BPN' and purified by reverse-phase HPLC. This approach facilitated the rapid purification of a set of 43 alanine replacement mutants whose relative affinity for GP IIb-IIIa was measured by competition with immobilized kistrin and by inhibition of platelet aggregation in human platelet-rich plasma. Alanine replacements at R49, G50, and D51 led to weaker inhibitors of platelet aggregation by 90-fold, 2-fold, and >200-fold, respectively. The conservative D51E mutant was still >100-fold less potent whereas R49K had a minor effect (1.8-fold), implying the critical nature of the aspartate for high affinity binding. However, mutations outside of the RGD region led to proteins indistinguishable from kistrin, suggesting no substantial secondary binding interactions. Furthermore, reduced kistrin is not active. We therefore propose that a favorable conformation of the RGD region alone is responsible for the high affinity binding of kistrin to GP IIb-IIIa. © 1993 Wiley-Liss, Inc.