Average assignment method for predicting the stability of protein mutants

Prediction of protein stability upon amino acid substitutions is an important problem in molecular biology and it will be helpful for designing stable mutants. In this work, we have analyzed the stability of protein mutants using three different data sets of 1791, 1396, and 2204 mutants, respectively, for thermal stability (DeltaTm), free energy change due to thermal (DeltaDeltaG), and denaturant denaturations (DeltaDeltaGH2O), obtained from the ProTherm database. We have classified the mutants into 380 possible substitutions and assigned the stability of each mutant using the information obtained with similar type of mutations. We observed that this assignment could distinguish the stabilizing and destabilizing mutants to an accuracy of 70-80% at different measures of stability. Further, we have classified the mutants based on secondary structure and solvent accessibility (ASA) and observed that the classification significantly improved the accuracy of prediction. The classification of mutants based on helix, strand, and coil distinguished the stabilizing/destabilizing mutants at an average accuracy of 82% and the correlation is 0.56; information about the location of residues at the interior, partially buried, and surface regions of a protein correctly identified the stabilizing/destabilizing residues at an average accuracy of 81% and the correlation is 0.59. The nine subclassifications based on three secondary structures and solvent accessibilities improved the accuracy of assigning stabilizing/destabilizing mutants to an accuracy of 84-89% for the three data sets. Further, the present method is able to predict the free energy change (DeltaDeltaG) upon mutations within a deviation of 0.64 kcal/mol. We suggest that this method could be used for predicting the stability of protein mutants.     

Human-readable rule generator for integrating amino scid sequence information and stability of mutant proteins

Most of the bioinformatics tools developed for predicting mutant protein stability appear as a black box and the relationship between amino acid sequence/structure and stability is hidden to the users. We have addressed this problem and developed a human-readable rule generator for integrating the knowledge of amino acid sequence and experimental stability change upon single mutation. Using information about the original residue, substituted residue, and three neighboring residues, classification rules have been generated to discriminate the stabilizing and destabilizing mutants and explore the basis for experimental data. These rules are human readable, and hence, the method enhances the synergy between expert knowledge and computational system. Furthermore, the performance of the rules has been assessed on a nonredundant data set of 1,859 mutants and we obtained an accuracy of 80 percent using cross validation. The results showed that the method could be effectively used as a tool for both knowledge discovery and predicting mutant protein stability. We have developed a Web for classification rule generator and it is freely available at http://bioinformatics.myweb.hinet.net/irobot.htm.     

Prediction of protein mutant stability using classification and regression tool

Prediction of protein stability upon amino acid substitutions is an important problem in molecular biology and the solving of which would help for designing stable mutants. In this work, we have analyzed the stability of protein mutants using two different datasets of 1396 and 2204 mutants obtained from ProTherm database, respectively for free energy change due to thermal (DeltaDeltaG) and denaturant denaturations (DeltaDeltaG(H(2)O)). We have used a set of 48 physical, chemical energetic and conformational properties of amino acid residues and computed the difference of amino acid properties for each mutant in both sets of data. These differences in amino acid properties have been related to protein stability (DeltaDeltaG and DeltaDeltaG(H(2)O)) and are used to train with classification and regression tool for predicting the stability of protein mutants. Further, we have tested the method with 4 fold, 5 fold and 10 fold cross validation procedures. We found that the physical properties, shape and flexibility are important determinants of protein stability. The classification of mutants based on secondary structure (helix, strand, turn and coil) and solvent accessibility (buried, partially buried, partially exposed and exposed) distinguished the stabilizing/destabilizing mutants at an average accuracy of 81% and 80%, respectively for DeltaDeltaG and DeltaDeltaG(H(2)O). The correlation between the experimental and predicted stability change is 0.61 for DeltaDeltaG and 0.44 for DeltaDeltaG(H(2)O). Further, the free energy change due to the replacement of amino acid residue has been predicted within an average error of 1.08 kcal/mol and 1.37 kcal/mol for thermal and chemical denaturation, respectively. The relative importance of secondary structure and solvent accessibility, and the influence of the dataset on prediction of protein mutant stability have been discussed.     

The consensus concept for thermostability engineering of proteins: further proof of concept

Previously, we calculated a consensus amino acid sequence from 13 homologous fungal phytases. A synthetic gene was constructed and recombinantly expressed. Surprisingly, consensus phytase-1 was 15-26 degrees C more thermostable than all parent phytases used in its design [Lehmann et al. (2000)Protein Eng., 13, 49-57]. In the present study, inclusion of six further phytase sequences in the amino acid sequence alignment resulted in the replacement of 38 amino acid residues in either one or both of the new consensus phytases-10 and -11. Since consensus phytase-10, again, was 7.4 degrees C more thermostable than consensus phytase-1, the thermostability effects of most of the 38 amino acid substitutions were tested by site-directed mutagenesis. Both stabilizing and destabilizing mutations were identified, but all affected the stability of the enzyme by <3 degrees C. The combination of all stabilizing amino acid exchanges in a multiple mutant of consensus phytase-1 increased the unfolding temperature from 78.0 to 88.5 degrees C. Likewise, back-mutation of four destabilizing amino acids and introduction of an additional stabilizing amino acid in consensus phytase-10 further increased the unfolding temperature from 85.4 to 90.4 degrees C. The thermostabilization achieved is the result of a combination of slight improvements from multiple amino acid exchanges rather than being the effect of a single or of just a few dominating mutations that have been introduced by chance. The present findings support the general validity of the consensus concept for thermostability engineering of proteins.     

Thermodynamic consequences of burial of polar and non-polar amino acid residues in the protein interior

Effects of amino acid substitutions at four fully buried sites of the ubiquitin molecule on the thermodynamic parameters (enthalpy, Gibbs energy) of unfolding were evaluated experimentally using differential scanning calorimetry. The same set of substitutions has been incorporated at each of four sites. These substitutions have been designed to perturb packing (van der Waals) interactions, hydration, and/or hydrogen bonding. From the analysis of the thermodynamic parameters for these ubiquitin variants we conclude that: (i) packing of non-polar groups in the protein interior is favorable and is largely defined by a favorable enthalpy of van der Waals interactions. The removal of one methylene group from the protein interior will destabilize a protein by approximately 5 kJ/mol, and will decrease the enthalpy of a protein by 12 kJ/mol. (ii) Burial of polar groups in the non-polar interior of a protein is highly destabilizing, and the degree of destabilization depends on the relative polarity of this group. For example, burial of Thr side-chain in the non-polar interior will be less destabilizing than burial of Asn side-chain. This decrease in stability is defined by a large enthalpy of dehydration of polar groups upon burial. (iii) The destabilizing effect of dehydration of polar groups upon burial can be compensated if these buried polar groups form hydrogen bonding. The enthalpy of this hydrogen bonding will compensate for the unfavorable dehydration energy and as a result the effect will be energetically neutral or even slightly stabilizing.     

Amino acid composition and thermal stability of proteins

This study investigates the relationship between the thermal stability of a globular protein and its amino acid composition. The method deals with the relationship between the amino acid compositions and melting points in a set of proteins by computing single-residue and group correlations. Groups of residues are shown to stabilize or destabilize the molecule against temperature. The stabilizing group consists of polar-charged residues and nonpolar residues possessing high surrounding hydrophobicity. The polar-uncharged residues destabilize the molecule against temperature, serine being the most destabilizing residue. A very high cooperativity exists among the stabilizing nonpolar residues suggesting that their characteristic clustering inside the globule may enhance the thermostability of a protein. In small globular proteins which act as single cooperative units, the melting temperature remains mainly a function of amino acid composition, whereas in complex molecules it depends on other factors also.     

Destabilizing and stabilizing forces to assess equilibrium during everyday activities

Postural stability is essential to functional activities. This paper presents a new model of dynamic stability which takes into account both the equilibrium associated with the body position over the base of support (destabilizing force) and the effort the subject needs to produce to keep his/her centre of mass inside the base of support (stabilizing force). The ratio between these two forces (destabilizing over stabilizing) is calculated to provide an overall index of stability for an individual. Preliminary results from data collected during walking at preferred and maximal safe speed in four older adults (aged from 64 to 84 yr) showed that both forces are lower for subjects with reduced maximal gait speed. In addition, the stabilizing force increases by 2–3 times from preferred to maximal speed, while the destabilizing force barely changes with gait speed. Overall, the model through the index of stability attributes lower dynamic stability to subjects with lower maximal gait speed. These preliminary results call for larger-scale studies to pursue the development and validation of the model and its application to different functional tasks.     

Mitochondrial protein turnover: role of the precursor intermediate peptidase Oct1 in protein stabilization

Most mitochondrial proteins are encoded in the nucleus as precursor proteins and carry N-terminal presequences for import into the organelle. The vast majority of presequences are proteolytically removed by the mitochondrial processing peptidase (MPP) localized in the matrix. A subset of precursors with a characteristic amino acid motif is additionally processed by the mitochondrial intermediate peptidase (MIP) octapeptidyl aminopeptidase 1 (Oct1), which removes an octapeptide from the N-terminus of the precursor intermediate. However, the function of this second cleavage step is elusive. In this paper, we report the identification of a novel Oct1 substrate protein with an unusual cleavage motif. Inspection of the Oct1 substrates revealed that the N-termini of the intermediates typically carry a destabilizing amino acid residue according to the N-end rule of protein degradation, whereas mature proteins carry stabilizing N-terminal residues. We compared the stability of intermediate and mature forms of Oct1 substrate proteins in organello and in vivo and found that Oct1 cleavage increases the half-life of its substrate proteins, most likely by removing destabilizing amino acids at the intermediate's N-terminus. Thus Oct1 converts unstable precursor intermediates generated by MPP into stable mature proteins.     

Stabilizing and destabilizing clusters in the hydrophobic core of long two-stranded alpha-helical coiled-coils

Detailed sequence analyses of the hydrophobic core residues of two long two-stranded alpha-helical coiled-coils that differ dramatically in sequence, function, and length were performed (tropomyosin of 284 residues and the coiled-coil domain of the myosin rod of 1086 residues). Three types of regions were present in the hydrophobic core of both proteins: stabilizing clusters and destabilizing clusters, defined as three or more consecutive core residues of either stabilizing (Leu, Ile, Val, Met, Phe, and Tyr) or destabilizing (Gly, Ala, Cys, Ser, Thr, Asn, Gln, Asp, Glu, His, Arg, Lys, and Trp) residues, and intervening regions that consist of both stabilizing and destabilizing residues in the hydrophobic core but no clusters. Subsequently, we designed a series of two-stranded coiled-coils to determine what defines a destabilizing cluster and varied the length of the destabilizing cluster from 3 to 7 residues to determine the length effect of the destabilizing cluster on protein stability. The results showed a dramatic destabilization, caused by a single Leu to Ala substitution, on formation of a 3-residue destabilizing cluster (DeltaT(m) of 17-21 degrees C) regardless of the stability of the coiled-coil. Any further substitution of Leu to Ala that increased the size of the destabilizing cluster to 5 or 7 hydrophobic core residues in length had little effect on stability (DeltaT(m) of 1.4-2.8 degrees C). These results suggested that the contribution of Leu to protein stability is context-dependent on whether the hydrophobe is in a stabilizing cluster or its proximity to neighboring destabilizing and stabilizing clusters.     

Physicochemical consequences of amino acid variations that contribute to fibril formation by immunoglobulin light chains

The most common form of systemic amyloidosis originates from antibody light chains. The large number of amino acid variations that distinguish amyloidogenic from nonamyloidogenic light chain proteins has impeded our understanding of the structural basis of light-chain fibril formation. Moreover, even among the subset of human light chains that are amyloidogenic, many primary structure differences are found. We compared the thermodynamic stabilities of two recombinant kappa4 light-chain variable domains (V(L)s) derived from amyloidogenic light chains with a V(L) from a benign light chain. The amyloidogenic V(L)s were significantly less stable than the benign V(L). Furthermore, only the amyloidogenic V(L)s formed fibrils under native conditions in an in vitro fibril formation assay. We used site-directed mutagenesis to examine the consequences of individual amino acid substitutions found in the amyloidogenic V(L)s on stability and fibril formation capability. Both stabilizing and destabilizing mutations were found; however, only destabilizing mutations induced fibril formation in vitro. We found that fibril formation by the benign V(L) could be induced by low concentrations of a denaturant. This indicates that there are no structural or sequence-specific features of the benign V(L) that are incompatible with fibril formation, other than its greater stability. These studies demonstrate that the V(L) beta-domain structure is vulnerable to destabilizing mutations at a number of sites, including complementarity determining regions (CDRs), and that loss of variable domain stability is a major driving force in fibril formation.     

The bacterial N‐end rule pathway: expect the unexpected

The N‐end rule pathway is a highly conserved process that operates in many different organisms. It relates the metabolic stability of a protein to its N‐terminal amino acid. Consequently, amino acids are described as either ‘stabilizing’ or ‘destabilizing’. Destabilizing residues are organized into three hierarchical levels: primary, secondary, and in eukaryotes – tertiary. Secondary and tertiary destabilizing residues act as signals for the post‐translational modification of the target protein, ultimately resulting in the attachment of a primary destabilizing residue to the N‐terminus of the protein. Regardless of their origin, proteins containing N‐terminal primary destabilizing residues are recognized by a key component of the pathway. In prokaryotes, the recognition component is a specialized adaptor protein, known as ClpS, which delivers target proteins directly to the ClpAP protease for degradation. In contrast, eukaryotes use a family of E3 ligases, known as UBRs, to recognize and ubiquitylate their substrates resulting in their turnover by the 26S proteasome. While the physiological role of the N‐end rule pathway is largely understood in eukaryotes, progress on the bacterial pathway has been slow. However, new interest in this area of research has invigorated several recent advances, unlocking some of the secrets of this unique proteolytic pathway in prokaryotes.     

How protein stability and new functions trade off

Numerous studies have noted that the evolution of new enzymatic specificities is accompanied by loss of the protein's thermodynamic stability (DeltaDeltaG), thus suggesting a tradeoff between the acquisition of new enzymatic functions and stability. However, since most mutations are destabilizing (DeltaDeltaG>0), one should ask how destabilizing mutations that confer new or altered enzymatic functions relative to all other mutations are. We applied DeltaDeltaG computations by FoldX to analyze the effects of 548 mutations that arose from the directed evolution of 22 different enzymes. The stability effects, location, and type of function-altering mutations were compared to DeltaDeltaG changes arising from all possible point mutations in the same enzymes. We found that mutations that modulate enzymatic functions are mostly destabilizing (average DeltaDeltaG = +0.9 kcal/mol), and are almost as destabilizing as the "average" mutation in these enzymes (+1.3 kcal/mol). Although their stability effects are not as dramatic as in key catalytic residues, mutations that modify the substrate binding pockets, and thus mediate new enzymatic specificities, place a larger stability burden than surface mutations that underline neutral, non-adaptive evolutionary changes. How are the destabilizing effects of functional mutations balanced to enable adaptation? Our analysis also indicated that many mutations that appear in directed evolution variants with no obvious role in the new function exert stabilizing effects that may compensate for the destabilizing effects of the crucial function-altering mutations. Thus, the evolution of new enzymatic activities, both in nature and in the laboratory, is dependent on the compensatory, stabilizing effect of apparently "silent" mutations in regions of the protein that are irrelevant to its function.     

On the Stability of Parainfluenza Virus 5 F Proteins

The crystal structure of the F protein (prefusion form) of the paramyxovirus parainfluenza virus 5 (PIV5) WR isolate was determined. We investigated the basis by which point mutations affect fusion in PIV5 isolates W3A and WR, which differ by two residues in the F ectodomain. The P22 stabilizing site acts through a local conformational change and a hydrophobic pocket interaction, whereas the S443 destabilizing site appears sensitive to both conformational effects and amino acid charge/polarity changes.     

Protein stability and molecular adaptation to extreme conditions.

Proteins, due to the delicate balance of stabilizing and destabilizing interactions, are only marginally stable. Adaptation to extreme environments tends to shift the 'mesophilic' characteristics of proteins to the respective extremes of temperature, hydrostatic pressure, pH and salinity, such that, under the mutual physiological conditions, the molecular properties are similar regarding overall topology, flexibility and solvation. Enhanced intrinsic stability requires only minute local structural changes so that general strategies of stabilization cannot be established. Apart from mutative changes of amino-acid sequences, extrinsic factors (or cellular components) may be involved in 'extremophilic adaptation'. The molecular basis of acidophilic, alkalophilic and barophilic adaptation is still obscure. Mechanisms of enhanced thermal stability involve improved packing density, as well as specific local interactions. In halophiles, water and salt binding of the intrinsically stable protein inventory is accomplished by favoring acidic over basic amino acid residues and decreased hydrophobicity. General limits of viability are: (a) the susceptibility of the covalent structure of the polypeptide chain toward hydrolysis or hydrothermal degradation; (b) the competition of extreme solvent parameters with the weak electrostatic and hydrophobic interactions involved in protein stabilization; (c) perturbations of the folding and assembly of proteins; and (d) 'dislocation' of biochemical pathways due to effects of extreme conditions on the intricate network of metabolic reactions.     

Stability of 100 homo and heterotypic coiled-coil a-a' pairs for ten amino acids (A, L, I, V, N, K, S, T, E, and R)

We present the thermal stability monitored by circular dichroism (CD) spectroscopy at 222 nm of 100 heterodimers that contain all possible coiled-coil a-a' pairs for 10 amino acids (I, V, L, N, A, K S, T, E, and R). This includes the stability of 36 heterodimers for 6 amino acids (I, V, L, N, A, and K) previously described and 64 new heterodimers including the 4 amino acids (S, T, E, and R). We have calculated a double mutant alanine thermodynamic cycle to determine a-a' pair coupling energies to evaluate which a-a' pairs encourage specific dimerization partners. The four new homotypic a-a' pairs (T-T, S-S, R-R, E-E) are repulsive relative to A-A and have destabilizing coupling energies. Among the 90 heterotypic a-a' pairs, the stabilizing coupling energies contain lysine or arginine paired with either an aliphatic or a polar amino acid. The range in coupling energies for each amino acid reveals its potential to regulate dimerization specificity. The a-a' pairs containing isoleucine and asparagine have the greatest range in coupling energies and thus contribute dramatically to dimerization specificity, which is to encourage homodimerization. In contrast, the a-a' pairs containing charged amino acids (K, R, and E) show the least range in coupling energies and promiscuously encourage heterodimerization.     

Electrostatic interactions in leucine zippers: thermodynamic analysis of the contributions of Glu and His residues and the effect of mutating salt bridges

Electrostatic interactions play a complex role in stabilizing proteins. Here, we present a rigorous thermodynamic analysis of the contribution of individual Glu and His residues to the relative pH-dependent stability of the designed disulfide-linked leucine zipper AB(SS). The contribution of an ionized side-chain to the pH-dependent stability is related to the shift of the pK(a) induced by folding of the coiled coil structure. pK(a)(F) values of ten Glu and two His side-chains in folded AB(SS) and the corresponding pK(a)(U) values in unfolded peptides with partial sequences of AB(SS) were determined by 1H NMR spectroscopy: of four Glu residues not involved in ion pairing, two are destabilizing (-5.6 kJ mol(-1)) and two are interacting with the positive alpha-helix dipoles and are thus stabilizing (+3.8 kJ mol(-1)) in charged form. The two His residues positioned in the C-terminal moiety of AB(SS) interact with the negative alpha-helix dipoles resulting in net stabilization of the coiled coil conformation carrying charged His (-2.6 kJ mol(-1)). Of the six Glu residues involved in inter-helical salt bridges, three are destabilizing and three are stabilizing in charged form, the net contribution of salt-bridged Glu side-chains being destabilizing (-1.1 kJ mol(-1)). The sum of the individual contributions of protonated Glu and His to the higher stability of AB(SS) at acidic pH (-5.4 kJ mol(-1)) agrees with the difference in stability determined by thermal unfolding at pH 8 and pH 2 (-5.3 kJ mol(-1)). To confirm salt bridge formation, the positive charge of the basic partner residue of one stabilizing and one destabilizing Glu was removed by isosteric mutations (Lys-->norleucine, Arg-->norvaline). Both mutations destabilize the coiled coil conformation at neutral pH and increase the pK(a) of the formerly ion-paired Glu side-chain, verifying the formation of a salt bridge even in the case where a charged side-chain is destabilizing. Because removing charges by a double mutation cycle mainly discloses the immediate charge-charge effect, mutational analysis tends to overestimate the overall energetic contribution of salt bridges to protein stability.     

Evolution of Amino Acid Propensities under Stability-Mediated Epistasis

Site-specific amino acid preferences are influenced by the genetic background of the protein. The preferences for resident amino acids are expected to, on average, increase over time because of replacements at other sites—a nonadaptive phenomenon referred to as the “evolutionary Stokes shift.” Alternatively, decreases in resident amino acid propensity have recently been viewed as evidence of adaptations to external environmental changes. Using population genetics theory and thermodynamic stability constraints, we show that nonadaptive evolution can lead to both positive and negative shifts in propensities following the fixation of an amino acid, emphasizing that the detection of negative shifts is not conclusive evidence of adaptation. By examining propensity shifts from when an amino acid is first accepted at a site until it is subsequently replaced, we find that ≈50% of sites show a decrease in the propensity for the newly resident amino acid while the remaining sites show an increase. Furthermore, the distributions of the magnitudes of positive and negative shifts were comparable. Preferences were often conserved via a significant negative autocorrelation in propensity changes—increases in propensities often followed by decreases, and vice versa. Lastly, we explore the underlying mechanisms that lead propensities to fluctuate. We observe that stabilizing replacements increase the mutational tolerance at a site and in doing so decrease the propensity for the resident amino acid. In contrast, destabilizing substitutions result in more rugged fitness landscapes that tend to favor the resident amino acid. In summary, our results characterize propensity trajectories under nonadaptive stability-constrained evolution against which evidence of adaptations should be calibrated.     

Amino acid propensities for the collagen triple-helix

Determination of the tendencies of amino acids to form alpha-helical and beta-sheet structures has been important in clarifying stabilizing interactions, protein design, and the protein folding problem. In this study, we have determined for the first time a complete scale of amino acid propensities for another important protein motif: the collagen triple-helix conformation with its Gly-X-Y repeating sequence. Guest triplets of the form Gly-X-Hyp and Gly-Pro-Y are used to quantitate the conformational propensities of all 20 amino acids for the X and Y positions in the context of a (Gly-Pro-Hyp)(8) host peptide. The rankings for both the X and Y positions show the highly stabilizing nature of imino acids and the destabilizing effects of Gly and aromatic residues. Many residues show differing propensities in the X versus Y position, related to the nonequivalence of these positions in terms of interchain interactions and solvent exposure. The propensity of amino acids to adopt a polyproline II-like conformation plays a role in their triple-helix rankings, as shown by a moderate correlation of triple-helix propensity with frequency of occurrence in polyproline II-like regions. The high propensity of ionizable residues in the X position suggests the importance of interchain hydrogen bonding directly or through water to backbone carbonyls or hydroxyprolines. The low propensity of side chains with branching at the C(delta) in the Y position supports models suggesting these groups block solvent access to backbone C=O groups. These data provide a first step in defining sequence-dependent variations in local triple-helix stability and binding, and are important for a general understanding of side chain interactions in all proteins.     

Use of amber suppressors to investigate the thermostability of Bacillus licheniformis alpha-amylase. Amino acid replacements at 6 histidine residues reveal a critical position at His-133

A set of 12 Escherichia coli suppressor tRNAs, inserting different amino acids in response to an amber codon, has been used to create rapidly numerous protein variants of a thermostable amylase; by site-directed mutagenesis, amber mutations were first introduced into Bacillus licheniformis alpha-amylase gene at position His35, His133, His247, His293, His406, or His450; genes carrying one or two amber mutations were then expressed in the different suppressor strains, generating over 100 amylase variants with predicted amino acid changes that could be tested for thermostability. Within the detection limits of the assays, amino acid replacements at five histidine positions had no significant effect. In contrast, suppressed variants substituted at residue His133 clearly exhibited modified thermostability and could be either less stable or more stable than the wild-type amylase, depending on the amino acid inserted at this position; comparison of the variants indicates that the hydrophobicity of the substituting residue is an important but not a determinant factor of stabilization. The effect of the most stabilizing and destabilizing amino acid substitutions, His133 to Tyr and to Pro, respectively, were confirmed by introducing the corresponding missense mutations into the gene sequence. The advantages and limits of informational suppression in protein stability studies are discussed as well as structural features involved in the thermostability of B. licheniformis alpha-amylase.