Desk-top analysis of the structural stability of various point mutations introduced into ribonuclease H

Strategy for Stabilization of CutA1 Proteins Due to Ion-Ion Interactions at Temperatures of over 100 °C

In order to elucidate the contribution of charged residues to protein stabilization at temperatures of over 100 °C, we constructed many mutants of the CutA1 protein ( EcCutA1) from Escherichia coli. The goal was to see if one can achieve the same stability as for a CutA1 from hyperthermophile Pyrococcus horikoshii that has the denaturation temperature near 150 °C. The hydrophobic mutant of EcCutA1 ( Ec0VV) with denaturation temperature ( T_d) of 113.2 °C was used as a template for mutations. The highest T_d of Ec0VV mutants substituted by a single charged residue was 118.4 °C. Multiple ion mutants were also constructed by combination of single mutants and found to have an increased thermostability. The highest stability of multiple mutants was a mutant substituted by nine charged residues that had a T_d of 142.2 °C. To evaluate the energy of ion-ion interactions of mutant proteins, we used the structural ensemble obtained by a molecular dynamics simulation at 300 K. The T_d of ionic mutants linearly increases with the increments of the computed energy of ion-ion interactions for ionic mutant proteins even up to the temperatures near 140 °C, suggesting that ion-ion interactions cumulatively contribute to the stabilization of a protein at high temperatures.

Stabilization of E. coli Ribonuclease HI by the ‘stability profile of mutant protein’ (SPMP)-inspired random and non-random mutagenesis

The change in the structural stability of Escherichia coli ribonuclease HI (RNase HI) due to single amino acid substitutions has been estimated computationally by the stability profile of mutant protein (SPMP) [Ota, M., Kanaya, S. Nishikawa, K., 1995. Desk-top analysis of the structural stability of various point mutations introduced into ribonuclease H. J. Mol. Biol. 248, 733-738]. As well, an effective strategy using random mutagenesis and genetic selection has been developed to obtain E. coli RNase HI mutants with enhanced thermostability [Haruki, M., Noguchi, E., Akasako, A., Oobatake, M., Itaya, M., Kanaya, S., 1994. A novel strategy for stabilization of Escherichia coli ribonuclease HI involving a screen for an intragenic suppressor of carboxyl-terminal deletions. J. Biol. Chem. 269, 26904-26911]. In this study, both methods were combined: random mutations were individually introduced to Lys99-Val101 on the N-terminus of the alpha-helix IV and the preceding beta-turn, where substitutions of other amino acid residues were expected to significantly increase the stability from SPMP, and then followed by genetic selection. Val101 to Ala, Gln, and Arg mutations were selected by genetic selection. The Val101-->Ala mutation increased the thermal stability of E. coli RNase HI by 2.0 degrees C in Tm at pH 5.5, whereas the Val101-->Gln and Val101-->Arg mutations decreased the thermostability. Separately, the Lys99-->Pro and Asn100-->Gly mutations were also introduced directly. The Lys99-->Pro mutation increased the thermostability of E. coli RNase HI by 1.8 degrees C in Tm at pH 5.5, whereas the Asn100-->Gly mutation decreased the thermostability by 17 degrees C. In addition, the Lys99-->Pro mutation altered the dependence of the enzymatic activity on divalent metal ions.

Stabilization mechanism of the tryptophan synthase alpha-subunit from Thermus thermophilus HB8: X-ray crystallographic analysis and calorimetry

In order to elucidate the thermo-stabilization mechanism of the tryptophan synthase alpha-subunit from the extreme thermophile Thermus thermophilus HB8 (Tt-alpha-subunit), its crystal structure was determined and its stability was examined using DSC. The results were compared to those of other orthologs from mesophilic and hyperthermophilic organisms. The denaturation temperature of the Tt-alpha-subunit was higher than that of the alpha-subunit from S. typhimurium (St-alpha-subunit) but lower than that of the alpha-subunit from P. furiosus (Pf-alpha-subunit). Specific denaturation enthalpy and specific denaturation heat capacity values of the Tt-alpha-subunit were the lowest among the three proteins, suggesting that entropy effects are responsible for the stabilization of the Tt-alpha-subunit. Based on a structural comparison with the St-alpha-subunit, two deletions in loop regions, an increase in the number of ion pairs and a decrease in cavity volume seem to be responsible for the stabilization of the Tt-alpha-subunit. The results of structural comparison suggest that the native structure of the Tt-alpha-subunit is better adapted to an ideally stable structure than that of the St-alpha-subunit, but worse than that of the Pf-alpha-subunit. The results of calorimetry suggest that the residual structure of the Tt-alpha-subunit in the denatured state contributes to the stabilization.

Design of λ Cro Fold: Solution Structure of a Monomeric Variant of the De Novo Protein

One of the classical DNA-binding proteins, bacteriophage lambda Cro, forms a homodimer with a unique fold of alpha-helices and beta-sheets. We have computationally designed an artificial sequence of 60 amino acid residues to stabilize the backbone tertiary structure of the lambda Cro dimer by simulated annealing using knowledge-based structure-sequence compatibility functions. The designed amino acid sequence has 25% identity with that of natural lambda Cro and preserves Phe58, which is important for formation of the stably folded structure of lambda Cro. The designed dimer protein and its monomeric variant, which was redesigned by the insertion of a beta-hairpin sequence at the C-terminal region to prevent dimerization, were synthesized and biochemically characterized to be well folded. The designed protein was monomeric under a wide range of protein concentrations and its solution structure was determined by NMR spectroscopy. The solved structure is similar to that of a monomeric variant of natural lambda Cro with a root-mean-square deviation of the polypeptide backbones at 2.1A and has a well-packed protein core. Thus, our knowledge-based functions provide approximate but essential relationships between amino acid sequences and protein structures, and are useful for finding novel sequences that are foldable into a given target structure.

Prediction of protein thermostability with a direction- and distance-dependent knowledge-based potential

The increasing use of enzymes in industrial processes and the importance of understanding protein folding and stability have led to several attempts to predict and quantify the effect of every possible amino acid exchange (mutation) on the thermostability of proteins. In this article we describe a knowledge-based discrimination function that acts as a fast and reliable guide in protein engineering and optimization. The function used consists of two parts, a pairwise energy function based on a distance- and direction-dependent atomic description of the amino acid environment, and a torsion angle energy function. In a first step a training set of 11 proteins including 646 mutant proteins with experimentally determined thermostability was used to optimize the knowledge-based energy functions. The resulting potential function was then tested using a test mutant database consisting of 918 various point mutations introduced in 27 proteins. The best correlation coefficient obtained for the experimental data and the predicted thermostability for the training set is r = 0.81 (561 data points). A total of 76% of the mutations could be predicted correctly as being either stabilizing or destabilizing. The results for the test set are r = 0.74 (747 data points) and 72%, respectively. The global correlation over the combined data (1308 mutants) obtained is 0.78.

The crystal structure of the tryptophan synthase beta subunit from the hyperthermophile Pyrococcus furiosus. Investigation of stabilization factors

The structure of the tryptophan synthase beta2 subunit (Pfbeta2) from the hyperthermophile, Pyrococcus furiosus, was determined by X-ray crystallographic analysis at 2.2 A resolution, and its stability was examined by DSC. This is the first report of the X-ray structure of the tryptophan synthase beta2 subunit alone, although the structure of the tryptophan synthase alpha2beta2 complex from Salmonella typhimurium has already been reported. The structure of Pfbeta2 was essentially similar to that of the beta2 subunit (Stbeta2) in the alpha2beta2 complex from S. typhimurium. The sequence alignment with secondary structures of Pfbeta and Stbeta in monomeric form showed that six residues in the N-terminal region and three residues in the C-terminal region were deleted in Pfbeta, and one residue at Pro366 of Stbeta and at Ile63 of Pfbeta was inserted. The denaturation temperature of Pfbeta2 was higher by 35 degrees C than the reported values from mesophiles at approximately pH 8. On the basis of structural information on both proteins, the analyses of the contributions of each stabilization factor indicate that: (a) the higher stability of Pfbeta2 is not caused by either a hydrophobic interaction or an increase in ion pairs; (b) the number of hydrogen bonds involved in the main chains of Pfbeta is greater by about 10% than that of Stbeta, indicating that the secondary structures of Pfbeta are more stabilized than those of Stbeta and (c) the sequence of Pfbeta seems to be better fitted to an ideally stable structure than that of Stbeta, as assessed from X-ray structure data.

Large-scale prediction of protein geometry and stability changes for arbitrary single point mutations

We have developed a method to both predict the geometry and the relative stability of point mutants that may be used for arbitrary mutations. The geometry optimization procedure was first tested on a new benchmark of 2141 ordered pairs of X-ray crystal structures of proteins that differ by a single point mutation, the largest data set to date. An empirical energy function, which includes terms representing the energy contributions of the folded and denatured proteins and uses the predicted mutant side chain conformation, was fit to a training set consisting of half of a diverse set of 1816 experimental stability values for single point mutations in 81 different proteins. The data included a substantial number of small to large residue mutations not considered by previous prediction studies. After removing 22 (approximately 2%) outliers, the stability calculation gave a standard deviation of 1.08 kcal/mol with a correlation coefficient of 0.82. The prediction method was then tested on the remaining half of the experimental data, giving a standard deviation of 1.10 kcal/mol and covariance of 0.66 for 97% of the test set. A regression fit of the energy function to a subset of 137 mutants, for which both native and mutant structures were available, gave a prediction error comparable to that for the complete training set with predicted side chain conformations. We found that about half of the variation is due to conformation-independent residue contributions. Finally, a fit to the experimental stability data using these residue parameters exclusively suggests guidelines for improving protein stability in the absence of detailed structure information.

Effect of mutations at Glu160 and Val198 on the thermostability of lactate oxidase

We have obtained two types of thermostable mutant lactate oxidase - one that exhibited an E-to-G point mutation at position 160 (E160G) through error-prone PCR-based random mutagenesis, and another that exhibited an E-to-G mutation at position 160 and a V-to-I mutation at position 198 (E160G/V198I) through DNA shuffling-based random mutagenesis - both of which we have previously reported. Our molecular modeling of lactate oxidase suggests that the substitution of G for E at position 160 reduces the electrostatic repulsion between the negative charges of E160 and E130 in the (beta/alpha)8 barrel structure, but a thermal-inactivation experiment on the five kinds of single-mutant lactate oxidase at position 160 (E160A, E160Q, E160H, E160R, and E160K) showed that the side-chain volume of the amino acid at position 160 mainly contributes to the thermostability of lactate oxidase. We also produced V198I single-mutant lactate oxidase through site-directed mutagenesis, and analysed the thermostability of wild-type, V198I, E160G, and E160G/V198I lactate oxidase enzymes. The half-life of E160G/V198I lactate oxidase at 70 degrees C was about three times longer than that of E160G lactate oxidase, and was about 20 times longer than that of wild-type lactate oxidase. In contrast, the thermostability of the V198I lactate oxidase was almost identical to that of wild-type lactate oxidase. This indicates that the V198I mutation alone does not affect lactate oxidase thermostability, but does affect it when combined with the E160G mutation.

Prediction of catalytic residues in enzymes based on known tertiary structure, stability profile, and sequence conservation

The catalytic or functionally important residues of a protein are known to exist in evolutionarily constrained regions. However, the patterns of residue conservation alone are sometimes not very informative, depending on the homologous sequences available for a given query protein. Here, we present an integrated method to locate the catalytic residues in an enzyme from its sequence and structure. Mutations of functional residues usually decrease the activity, but concurrently often increase stability. Also, catalytic residues tend to occupy partially buried sites in holes or clefts on the molecular surface. After confirming these general tendencies by carrying out statistical analyses on 49 representative enzymes, these data together with amino acid conservation were evaluated. This novel method exhibited better sensitivity in the prediction accuracy than traditional methods that consider only the residue conservation. We applied it to some so-called "hypothetical" proteins, with known structures but undefined functions. The relationships among the catalytic, conserved, and destabilizing residues in enzymatic proteins are discussed.

Knowledge-based potential defined for a rotamer library to design protein sequences

A knowledge-based potential for a rotamer library was developed to design protein sequences. Protein side-chain conformations are represented by 56 templates. Each of their fitness to a given structural site-environment is evaluated by a combined function of the three knowledge-based terms, i.e. two-body side-chain packing, one-body hydration and local conformation. The number of matches between the native sequence and the structural site-environment in the database and that of the virtually settled mismatches, counted in advance, were transformed into the energy scores. In the best-14 test (assessment for the reproduction ability of the native rotamer on its structural site within a quarter of 56 fitness rank positions), the structural stability analysis on mutants of human and T4 lysozymes and the inverse-folding search by a structure profile against the sequence database, this function performs better than the function deduced with the conventional normalization and our previously developed function. Targeting various structural motifs, de novo sequence design was conducted with the function. The sequences thus obtained exhibit reasonable molecular masses and hydrophobic/hydrophilic patterns similar to the native sequences of the target and act as if they were the homologs to the target proteins in BLASTP search. This significant improvement is discussed in terms of the reference state for normalization and the crucial role of short-range repulsion to prohibit residue bumps.

PoPMuSiC, an algorithm for predicting protein mutant stability changes: application to prion proteins

A novel tool for computer-aided design of single-site mutations in proteins and peptides is presented. It proceeds by performing in silico all possible point mutations in a given protein or protein region and estimating the stability changes with linear combinations of database-derived potentials, whose coefficients depend on the solvent accessibility of the mutated residues. Upon completion, it yields a list of the most stabilizing, destabilizing or neutral mutations. This tool is applied to mouse, hamster and human prion proteins to identify the point mutations that are the most likely to stabilize their cellular form. The selected mutations are essentially located in the second helix, which presents an intrinsic preference to form beta-structures, with the best mutations being T183-->F, T192-->A and Q186-->A. The T183 mutation is predicted to be by far the most stabilizing one, but should be considered with care as it blocks the glycosylation of N181 and this blockade is known to favor the cellular to scrapie conversion. Furthermore, following the hypothesis that the first helix might induce the formation of hydrophilic beta-aggregates, several mutations that are neutral with respect to the structure's stability but improve the helix hydrophobicity are selected, among which is E146-->L. These mutations are intended as good candidates to undergo experimental tests.

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.

Feasibility in the inverse protein folding protocol

Methods for protein structure (3D)-sequence (1D) compatibility evaluation (threading) have been developed during the past decade. The protocol in which a sequence can recognize its compatible structure in the structural library (i.e., the fold recognition or the forward-folding search) is available for the structure prediction of new proteins. However, the reverse protocol, in which a structure recognizes its homologous sequences among a sequence database, named the inverse-folding search, is a more difficult application. In this study, we have investigated the feasibility of the latter approach. A structural library, composed of about 400 well-resolved structures with mutually dissimilar sequences, was prepared, and 163 of them had remote homologs in the library. We examined whether they could correctly seek their homologs by both forward- and inverse-folding searches. The results showed that the inverse-folding protocol is more effective than the forward-folding protocol, once the reference states of the compatibility functions are appropriately adjusted. This adjustment only slightly affects the ability of the forward-folding search. We noticed that the scoring, in which a given sequence is re-mounted onto a structure according to the 3D-1D alignment determined by the dynamic programming method, is only effective in the forward-folding protocol and not in the inverse-folding protocol. Namely, the inverse-folding search works significantly better with the score given by the 3D-1D alignment per se, rather than that obtained by the re-mounting. The implications of these results are discussed.

Structural requirement of highly-conserved residues in globins

Globins have remarkable sequence diversity, and yet maintain a common fold. In spite of the diversity, there are highly-conserved residues at several sites. The conserved residues were examined in terms of the structural stability, by employing the pseudo-energy functions of the structure/sequence compatibility method. The fitness of each residue type to the structural environment was evaluated at seven highly-conserved sites: the Leu (at the B10 site), Phe (CD1), and Leu (F4) residues were found to fit their respective sites due to hydrophobic interactions; Pro (C2) stabilizes the N-terminal edge of an alpha-helical structure; and Phe (CD4) is stabilized by backbone hydrogen-bonding to Phe (CD1). On the other hand, the other two residues, His (E7) and His (F8), are poorly suited to the sites from a structural viewpoint, suggesting that their conservation clearly results from a heme-related functional requirement. The invariant Phe residue (CD1) has been suggested to be important for supporting the heme. The present analysis revealed that this residue is also well suited to the site in terms of energy.

Comparison of protein structures using 3D profile alignment

A novel method for protein structure comparison using 3D profile alignment is presented. The 3D profile is a position-dependent scoring matrix derived from three-dimensional structures and is basically used to estimate sequence-structure compatibility for prediction of protein structure. Our idea is to compare two 3D profiles using a dynamic programming algorithm to obtain optimal alignment and a similarity score between them. When the 3D profile of hemoglobin was compared with each of the profiles in the library, which contained 325 profiles of representative structures, all the profiles of other globins were detected with relatively high scores, and proteins in the same structural class followed the globins. Exhaustive comparison of 3D profiles in the library was also performed to depict protein relatedness in the structure space. Using multidimensional scaling, a planar projection of points in the protein structure space revealed an overall grouping in terms of structural classes, i.e., all-alpha, all-beta, alpha/beta, and alpha+beta. These results differ in implication from those obtained by the conventional structure-structure comparison method. Differences are discussed with respect to the structural divergence of proteins in the course of molecular evolution.

Thermal stability of Escherichia coli ribonuclease HI and its active site mutants in the presence and absence of the Mg2+ ion. Proposal of a novel catalytic role for Glu48

Escherichia coli ribonuclease HI, which requires divalent cations (Mg2+ or Mn2+) for activity, was thermostabilized by 2.6-3.0 kcal/mol in the presence of the Mg2+, Mn2+, or Ca2+ ion, probably because the negative charge repulsion around the active site was canceled upon the binding of these metal ions. The dissociation constants were determined to be 0.71 mM for Mg2+, 0.035 mM for Mn2+, and 0.16 mM for Ca2+. Likewise, various active site mutants at Asp10, Glu48, Asp70, or Asp134 were thermostabilized by 0.4-3.0 kcal/mol in the presence of the Mg2+ ion, suggesting that this ion binds to these mutant proteins as well. The dissociation constants of Mg2+ were determined to be 9.8 mM for D10N, 1.1 mM for E48Q, 18.8 mM for D70N, and 1.8 mM for D134N. Thus, the mutation of Asp10 or Asp70 to Asn considerably impairs the Mg2+ binding, whereas the mutation of Glu48 to Gln or Asp134 to Asn does not. Comparison of the thermal stability of the mutant proteins with that of the wild-type protein in the absence of the Mg2+ ion suggests that the negative charge repulsion between Asp10 and Asp70 is responsible for the binding of the metal cofactor. Glu48 may be required to anchor a water molecule, which functions as a general acid.