A simple electrostatic criterion for predicting the thermal stability of proteins

Angel Mozo-Villarías1,2, Juan Cedano3 and Enrique Querol3

1 Departament de Ciències Mèdiques Bàsiques, Facultat de Medicina, Universitat de Lleida, Avda. Rovira Roure 44, 25198 Lleida and 3 Institut de Biotecnologia i Biomedicina, Departament de Bioquímica i Biologia Molecular, Universitat Autònoma de Barcelona, Bellaterra, Barcelona, Spain

2 To whom correspondence should be addressed. E-mail: angel.mozo{at}cmb.udl.es


    Abstract
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
The enhancement of protein thermostability is an important issue for both basic science and biotechnology purposes. We have developed a thermostability criterion for a protein in terms of a quasi-electric dipole moment (contributed by its charged residues) defined for an electric charge distribution whose total charge is not zero. It was found that minimization of the modulus of this dipole moment increased its thermal stability, as demonstrated by surveying these values in pairs of mesostable–thermostable homologous proteins and in mutations described in the literature. The analysis of these observations provides criteria for thermostabilization of a protein, by computing its dipole profile. This profile is obtained by direct substitution of each amino acid of the sequence by either a positive, negative or neutral amino acid, followed by a recalculation of the dipole moment. As an experimental example, these criteria were applied to a ß-glucanase to enhance its thermal stability.

Keywords: electric dipole moment/electrostatic interactions/hydrophobicity/thermostability


    Introduction
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
In recent years, there has been increasing interest in unraveling the specific mechanisms used by Nature to render proteins stable in the face of thermal inactivation. The interest in this issue involves questions that range from basic science to the design and manufacture of more efficient enzymes for industrial processes. At present, when confronted with the need to determine which amino acid is to be changed in order to enhance the thermal stability of a protein, most researchers still tend to opt for mainly casuistic criteria, since it is still difficult to follow definitive rules (Querol et al., 1996Go).

Of all the interactions that take place within proteins, electrostatic and hydrophobic are the ones that have probably captured the attention of most researchers (Querol et al., 1996Go; Mozo-Villarías and Querol, 2000Go; Sterner and Liebl, 2001Go; Vieille and Zeikus, 2001Go). As far as electrostatic forces are concerned, an increasing number of papers have established the enormous importance (both in number and constitution) of ion pair networks in maintaining thermal stability (Bashford, 1991Go; Honig and Nicholls, 1995Go; Vogt et al., 1997Go; Chin et al., 1999Go; Xiao and Honig, 1999Go). The idea that emerges from these observations is that of a protein which works as an electric network, mostly at the surface and around a hydrophobic core. The stronger and more complex the net is, the most resistant the protein is to destabilization. A large number of computer software tools based on different numerical approximations analyzing electrical interactions have demonstrated their utility for computing and establishing electrical characteristics, such as the distributions of electric field and electric potential in proteins (Honig and Nicholls, 1995Go). The present work suggests the use of magnitudes that characterize electrical properties of a protein in order to help determine which points of the protein are most sensitive in maintaining its 3D structure.

In this paper, we describe common traits shown by studies of protein thermostabilization involving electrical interactions. It must be noted that when extracting protein thermal information from the literature, it is difficult to choose interesting cases in a coherent manner, owing to the different criteria used by different researchers when reporting their results in terms of the thermal characteristics of the protein under study. The transition temperature, Tm, is considered to be the standard parameter for describing the thermal stability of a protein, yet many researchers prefer to present their results in terms of thermoresistance, which is the half-life of residual activity at a given temperature and therefore not necessarily related to Tm. This is a frequent case in studies involving enzymes for biotechnological purposes. Other authors simply provide information about optimal temperatures at which original organisms live, while there are those who directly provide values for {Delta}{Delta}G, that is, the change in free energy associated with thermostabilization.

The analysis presented in this paper provides an approximate but very simple way to monitor the outcome of any prediction of thermostabilization when a given mutation involving charged residues is proposed.


    Materials and methods
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
Theory

Hydrophobic density calculations. The center-of-mass of each protein was determined using the coordinates provided by the Protein Data Bank (PDB). The coordinates for each atom were re-calculated and used for subsequent calculations. In this way, the center-of-mass of the protein was the origin of coordinates for all subsequent computations.

Hydrophobic tensors equivalent to the inertia tensor were calculated for positive hydrophobicities (h+i) and negative hydrophobicities (h-i), according to the Eisenberg hydrophobicity scale (Eisenberg et al., 1982Go). On this scale, positive values indicate hydrophobicity and negative values indicate hydrophylicity. These tensors were subsequently diagonalized: . The radii of the equivalent ellipsoids were calculated from their respective diagonal elements. Thus, the radii of the ellipsoid corresponding to the positive hydrophobicity amino acids were calculated as


where are the diagonal elements of the negative electric charge distribution and is the sum of all the positive hydrophobicity amino acids in the protein. Subscript i denotes the three dimensions of space x, y, z. Similarly, for negative hydrophobicity were calculated.

The positive hydrophobicity ellipsoid is of particular interest in this study, since it provides a measure equivalent to the density of the hydrophobic packing of a protein, . This density is computed as

Dipole moment vectors. The classical definition of the electric dipole moment vector or first-order electric moment is given

by , where is the position vector of residue j with charge qj. Unless the total charge of the protein is zero, is not independent of the origin of coordinates. Since in most proteins the total charge is not zero, it is useful to define two quasi-dipole moments, and , as and , where and are the total negative and positive charges ( is the total charge) and


are the positive and negative centroids of the charge distributions, respectively.

Note that and differ, as , and when

.

Both and are independent of the origin of coordinates and therefore they are intrinsic parameters of the protein. As a consequence, these quasi-dipole moments provide a useful characterization of the dipolar nature of the charge distribution of a protein. Since the information provided by is essentially

the same as that of , in what follows is used to characterize the charge distribution in a protein and its modulus is designed as P and expressed in debye. We prefer to use these quasi-moments rather than the origin-dependent true moment since the purpose of this study was to compare moments among proteins with different thermal properties as described below and the use of D could give rise to meaningless conclusions.

All these magnitudes were easily calculated for the proteins described below using a desktop computer and employing a very simple function written in Mathematica.

Mesostable–thermostable protein pairs. Coordinates for mesostable–thermostable pairs of proteins, whose thermal properties are described in the literature, were obtained from the PDB. These proteins are listed in Table IGo. Here we report only those proteins whose sequences are complete. In some cases, proteins in their native state are complexes of several monomers, but only the coordinates of the monomers are reported in the PDB. In such cases, the oligomers were reconstructed by means of the Swiss-PDBViewer software, using their oligomeric homologues as templates. Proteins which lacked important fragments are not reported. In a few cases, mesostable–thermostable pairs also included a hyperthermostable member whose parameters were also calculated for comparison.


View this table:
[in this window]
[in a new window]
 
Table I. Electrical and hydrophobic characteristics in mesostable–thermostable pairs of proteins
 
Sets of mutants. Table IIGo shows significant cases of mutations that were collected according to the same criteria used for thermostable-mesostable pairs, in which {Delta}Tm >= 1°C. This collection of mutations was obtained from a previously reported mutant database (Querol et al., 1996Go), together with later reports on mutations that have appeared in the literature in recent years. When reporting, we tried to give preference to cases in which mutations resulted in crystallized proteins whose 3D structure could be resolved and thus have a PDB entry code. If not, the position of the mutant amino acid was taken to be the same as that of the mutated and an energy minimization was performed using the Swiss-PDBViewer software.


View this table:
[in this window]
[in a new window]
 
Table II. Electrical and hydrophobic characteristics in mutated proteins
 
Most of the mutations reported in the literature and referred to in this paper provide a wealth of mutants for each protein. Owing to space restrictions we report only the most significant mutants here. A complete list of all the mutants identified is available from the authors upon request.

Experimental

ß-Glucanase and its mutant N207D were obtained and purified by using experimental procedures described elsewhere (Pons et al., 1997Go). Transition temperatures (Tm) upon heating were measured spectrofluorimetrically on ß-glucanase wild-type and N207D mutant. Experiments were carried out with a Shimadzu RF2000 spectrofluorimeter.


    Results
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
Mesostable–thermostable pairs

Sets of mesostable–thermostable pairs of proteins were obtained by pooling all cases recently reported in the literature (Vogt et al., 1997Go; Karshikoff and Ladenstein, 1998Go; Xiao and Honig, 1999Go; Knochel et al., 2000Go; Kumar et al., 2000Go; Sterner and Liebl, 2001Go; Vieille and Zeikus, 2001Go; Klein et al., 2002; Nicholson et al., 2002Go; Robic et al., 2002Go). In the cases reported here, we did not take into account particular atomic or structural considerations. The only condition imposed was that the proteins reported should not lack any fragments or elements that could be of importance in the interpretation of electrical results. The electrical parameters described in the Materials and methods section were calculated for 37 groups of mesostable–thermostable (and/or hyperthermostable) proteins and are shown in Table IGo, together with their thermal and mechanical characteristics. For the computation of the electrical parameters, the charged amino acids were used [+1 for arginines and lysines, -1 for aspartates and glutamates and +0.5 for histidines (Spassov et al., 1997Go)]. Regardless of the thermal properties of each protein, it was found that the sizes of the equivalent electric ellipsoids (both positive and negative) were, as expected, somewhat larger than those of the inertia ellipsoids, since most charged amino acids tend to be located in the surface of the protein. Of the 37 pairs of proteins reported, 27 (73%) had their quasi-dipole moment value (P) lower than that of their mesostable counterpart. It was also observed that many cases also show an increase in {rho}+ concomitant with higher thermostability. Only two families (CheY and methionine aminopeptidase) showed both opposite effects simultaneously: increase in P and decrease in {rho}+ for higher thermostability.

We estimated the error in these parameters from observed differences in values obtained in different subunits of polymeric crystals. For example, we found that in glutamate dehydrogenase, the variability was of the order of 2% for these calculated parameters.

Point mutations

The considerations explained in the preceding section also apply for the proteins reported in Table IIGo and their respective point mutations (Dao-Pin et al., 1991Go; Meiering et al., 1992Go; Zhang et al., 1992Go; Kanaya et al., 1996Go; Predki et al., 1996Go; White et al., 1996Go; Malakauskas and Mayo, 1998Go; Grimsley et al., 1999Go; Lebbnik et al., 1999Go; Mikami et al., 1999Go; Peterson et al., 1999Go; Wray et al., 1999Go; Zhu et al., 1999Go; Arrizubieta and Polaina, 2000Go; Chen et al., 2000Go; González-Blasco et al., 2000Go; Ha et al., 2000Go; Hasegawa et al., 2000Go; Tchan et al., 2000Go; Martin et al., 2001Go; Ohmura et al., 2001Go; Pedone et al., 2001Go; Perl and Schmid, 2001Go; Shaw et al., 2001Go; Stewart et al., 2001Go; Sung et al., 2001Go; Trejo et al., 2001Go; Almog et al., 2002Go; Pokkuluri et al., 2002Go). This table reports 55 mutations carried out in 30 proteins involving electric changes. In 43 of the cases reported (78%) mutations showed a decrease in P when thermostability increased or, conversely, an increase in P when thermostability decreased, which is common among reported mutations.

Relative electric dipole moment versus transition temperature

In order to illustrate the sensitivity of the thermal stability of a protein with respect to its quasi-dipole moment, the relative change of P was plotted against the increase in Tm, {Delta}Tm. This relative change, {Delta}P/Pmeso, is calculated by taking the value of P of the mesostable protein as reference, in cases dealing with meso–thermostable pairs or the value of the wild-type when dealing with point mutations. Figure 1Go shows this relative change in P vs {Delta}Tm for those cases whose variations of thermal properties are reported in terms of Tm, for both meso–thermostable pairs of proteins and point mutations. Both plots show linear correlations between both increments, {Delta}P/Pmeso and {Delta}Tm. The linear fit characteristics of both populations are displayed in Table IIIGo. It should be noted that both fits show similar slope values around -0.011. The level of significance for both sets of proteins is p < 0.05 for meso–thermo pairs (n = 18) and p < 0.01 for mutated proteins (n = 39).



View larger version (13K):
[in this window]
[in a new window]
 
Fig. 1. The relative change of the quasi-electric dipole moment ({Delta}P/Pmeso) is plotted against the increment in transition temperature {Delta}Tm. Closed circles represent {Delta}P/Pmeso values from Table IGo; open circles are values taken from Table IIGo. Straight lines represent linear fits to data in Table IGo (continuous line) and Table IIGo (broken line). Dotted line represents the fitting of {Delta}P/Pmeso for mesostable–thermostable pairs of proteins calculated according to the method of Felder et al. (Felder et al., 1998Go).

 

View this table:
[in this window]
[in a new window]
 
Table III. Characteristics of linear plots of relative changes in the quasi-electric dipole moment with changes in Tm (Figure 1Go)
 
Similar calculations were also carried out using P->n instead

of P->p and very similar results were obtained in terms of {Delta}P/Pmeso.

Detailed analysis and profiles of some examples

Cold shock protein. In order to illustrate the use of these electric dipole moments in predicting thermostabilizing mutations, profiles of P were calculated throughout the cold shock protein sequence. Figure 2Go shows profiles of computed P for the mesostable–thermostable pair, when each amino acid is replaced by either a negative charge, a neutral amino acid or a positive charge throughout its sequences. Those substitutions in which wild-type P values were reduced should, in principle, be considered as candidates for stabilizing mutations. For example, mutations E3R and E66L produced a decrease in P. Perl and Schmidt performed these mutations on the mesostable cold shock protein 1csp with a resulting increase in Tm (Perl and Schmidt, 2001). A corresponding decrease in Tm was observed when reversed mutations were carried out on the thermostable cold shock protein 1c90. An increase in P was observed for R3E and L66E mutations.



View larger version (43K):
[in this window]
[in a new window]
 
Fig. 2. Quasi-electric dipole moment profiles of cold shock protein. Dotted lines, mesostable 1csp; continuous lines, thermostable 1c90. Profiles were obtained by replacing each amino acid by either a negative charge (red squares), a neutral amino acid (green circles) or a positive charge (blue triangles) and by recalculating the value of the dipole moment P after each change. At the beginning and end of the 1csp sequence, mutations to positive charges render values of P similar to those in native 1c90. Perl and Schmidt found a concomitant increase in Tm to values similar to 1c90 in their E3R and E66L mutants (Perl and Schmidt, 2001). Similarly, mutations to negative charges in the N- and C-domains of 1c90 (R3E and L66E) gave values of P similar to those in native 1csp.

 
Glucanase. It is clear that not all of the substitutions suggested by the P profile should be regarded as mutation candidates, as many of these amino acids may be important for protein functions. In the case of ß-glucanase, the resulting profile shows several minima corresponding to residues/stretches in positions 9–10, 50, 55, 80, 99, 104–109, 128–130, 145, 178 and 207 (Figure 3Go). Careful analysis of the protein structure by Rasmol and previous work on ß-glucanase mutants carried out in our laboratories (Juncosa et al., 1994Go; Pons et al., 1997Go) enabled us to exclude some potential substitutions as potential targets for replacement by a charged residue. These deductions were reached on the basis of structural/functional conservation: (a) positive/negative charged residues electrostaticaly interacting with a neighboring negative/positive charged residues/ions: Lys80 with Glu46; Lys155 with Asp150; Glu9 with Ca2+; Asp50 with Ca2+; Asp99 with Lys178; (b) residues belonging to the active site (general acid and base catalysts, nucleophiles or residues from the active site cleft), which we prefer not to replace or which lie on a key secondary structure: 104–109, 128–130 and 178–180 (Juncosa et al., 1994Go; Pons et al., 1997Go); (c) the remaining amino acids contain two putative targets: His148 and Asn207. Asp was chosen as a substitute for Asn207, because they are isosteric. For this mutation N207D, the value of P was found to change from 12.3 (wild-type) to 10.9 D, whereas Tm increased from 71.4 to 75.3°C.



View larger version (44K):
[in this window]
[in a new window]
 
Fig. 3. Quasi-electric dipole moment profile of ß-glucanase. Profiles were obtained as described in Figure 2Go. Substitutions to different charges are represented as follows: red circles, changes to negative charges; green squares, changes to neutral amino acids; blue circles, changes to positive charges. The straight line represents the original value of the pseudo-dipole moment: 12.26 D.

 
It should also be noted that, upon profiling and as shown empirically in our previous work, most of the positions initially considered as good candidates for mutations lay in external regions of the protein structure (Querol et al., 1996Go).


    Discussion
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
It is important to stress that the electrical characteristics of a protein, especially in terms of its dipolar nature, are determined not only by its charged amino acids, but also by its peptide bonds and secondary structures. However, we found that for any given protein, the total dipole moment contributed by these elements (excluding the charged amino acids) has a much lower value than that contributed by the charged amino acids and, what is most important in this study, its value does not change substantially among the mesostable and thermostable varieties or between the wild-type protein and the mutant.

We found an inverse linear correlation between the relative change in dipole moment and transition temperature increase. Such a correlation should be expected since the attractive Coulomb force, which is a cohesive force, exerted by both positive and negative centroids depends on the inverse square distance between them. This finding is in good agreement with current ideas about the role played by the electrostatic networks of ion pairs and their optimization in thermostable proteins, as pointed out by several investigators (Spassov et al., 1994Go; Vogt et al., 1997Go; Xiao and Honig, 1999Go). According to these authors, such a clustering of opposite charges brings the centers of positive and negative charges close in space and thus reduces the magnitude of the electrostatic dipole moment of the protein. A similar correlation is found when the mesostable–thermostable pairs of Table IGo are used to compute dipole moments according to the method described by Felder et al.(1998)Go by means of the true concept of dipole moment taking the precaution of being systematic with the choice of the origin of coordinates. However, the correlation coefficient found with their dipole moments was lower than those reported in Table IIIGo.

In the calculations of the electrostatic dipole moment, it was assumed that all ionizable groups are fully charged except for histidines that were considered to have +0.5 charge. Given the large set of proteins being considered, we propose that this assumption is acceptable in the light of recent results on the possibility that a point mutation might not change the net charge of the protein. Alexov et al.(2000)Go, working with replacements of a neutral group with ionizable in the reaction centers of several species of Rhodobacter, found that the replacements might not change the charge of the protein simply because the ionizable group could be uncharged owing to the specific environment in the protein or it might alter the charge of its neighbor.

It is worth mentioning that the attractive electrostatic force also depends on the inverse of the dielectric constant of the protein, which is also likely to decrease upon increasing the hydrophobic density. The difficulty in obtaining values for the dielectric constant of a protein has already been pointed out (Pitera et al. 2001Go; Schutz and Warshel, 2001Go), being a function of the hydrophobic amino acid composition of the protein. The link between the electrical and hydrodynamic forces operating in a protein is certainly an unresolved complex issue and one in which determination of the distribution of the dielectric constant in the protein is crucial. It may not even be possible to propose general rules to describe their mutual influence, since they may depend on the context of each protein. It is necessary to rely on the opportunistic nature of this equilibrium in which other forces, like hydrogen bonding, solvation and others that have been excluded from this study, are also present. However, in a first approximation, an increase in positive hydrophobic density, {rho}+, should imply a decrease in the average dielectric constant of the core of the protein and this would contribute to an increase in the strength of the attractive electric force, even in cases in which dipole moment increases. There are some cases reported in this paper in which we see thermal stabilizations in spite of a relative increase in P. At this point we can only suggest that when the core is more densely packed (as expressed by an increase in the positive hydrophobicity density, {Delta}{rho}+), the effect of increasing the separation between electric centroids (increase in P) may be compensated by an increase in the dielectric constant of the medium.

Another important aspect to be taken into account is the fact that the stabilization factors suggested here (minimization of P and maximization of {rho}+) are defined in a specific structural context involving relatively small structural variations within a protein. When major changes are to be considered, such as those simulated by Petrey and Honig (Petrey and Honig, 2000Go), where the sequence of a given protein was superimposed on a totally different fold, the effects described in this paper should not be expected necessarily to hold.

The interest in this analysis lies mostly in the way it effectively provides both electric and hydrophobic equivalent spheroids for the study protein, without the need for a detailed knowledge of atomic peculiarities such as steric constraints, the existence of water-filled or empty cavities, etc. This circumstance imposes certain limitations on the applicability of this analysis in the case of very asymmetric proteins, as exemplified in the case of aminotransferase (not reported here). In this protein, the association of an ellipsoid with such an irregular structure (‘S’-shaped) causes a considerable distortion of the model, which considers a protein as a spheroid possessing a compact hydrophobic nucleus, enclosed in an electrical layer.

Although the present state of knowledge does not allow a quantitative prediction of the degree of protein stability from its sequence and/or its 3D structure, or the amount of acquired stabilization (or destabilization) when a given mutation takes place, it is possible to check if putative specific replacements are likely to enhance the thermal stability of a protein.


    Acknowledgments
 
This research was supported by grants BIO2000-647 and BIO2001-2046 from the CICYT (Ministerio de Educación y Ciencia, Spain), by the Center de Refèrencia R+D de Biotecnologia de la Generalitrat de Catalunya and by La Paeria (Lleida City Hall). We thank Drs S.Marqusee and C.Klein for providing unpublished coordinates used in this study, Mr Malcolm Hayes for correcting the English and Mr Carmel Bonet for useful suggestions and comments.


    References
 Top
 Abstract
 Introduction
 Materials and methods
 Results
 Discussion
 References
 
Alexov,E., Miksovska,J., Baciou,L., Schiffer,M., Hanson,D., Sebba,P. and Gunner,M. (2000) Biochemistry, 39, 5940–5952.[CrossRef][ISI][Medline]

Almog,O., Gallagher, D.T., Ladner,J., Strausberg, S., Alexander,P., Bryan,P. and Gililand,G. (2002) J. Biol. Chem., 277, 27553–27558.[Abstract/Free Full Text]

Arrizubieta,M.J. and Polaina,J. (2000) J. Mol. Biol., 275, 28843–28848.

Bashford,D. (1991) Curr. Opin. Struct. Biol., 1, 175–184.[CrossRef]

Chen,J., Lu,Z., Sakon,J. and Stites,W. (2000) J. Mol. Biol., 303, 125–130.[CrossRef][ISI][Medline]

Chin,K., Sharp,K., Honig,B. and Pyle,A.-M. (1999) Nature Struct. Biol., 6, 1055–1061.[CrossRef][ISI][Medline]

Dao-Pin,S., Soderlind,E., Baase,W., Wozniak,J., Sauer,U. and Matthews,B. (1991) J. Mol. Biol., 221, 873–887.[CrossRef][ISI][Medline]

Eisenberg, D., Weiss, R., Terwilliger, T. and Wilcox, W. (1982) Faraday Symp. Chem. Soc., 17, 109–120.

Felder,C., Botti,S., Lifson,S., Silman,I. and Sussman,J. (1998) J. Mol. Graphics Modelling, 15, 318–327.[CrossRef][ISI]

González-Blasco,G., Sanz-Aparicio,J., González,B., Hermoso,J. and Polaina,J. (2000) J. Chem. Biol., 275, 13780–13712.[Abstract/Free Full Text]

Grimsley,G., Shaw,K., Fee,L., Alston,R., Huyghes-Despointes,B., Thurlkill,R., Scholz,R. and Pace,N. (1999) Protein Sci., 8, 1843–1849.[Abstract]

Ha,N.C., Oh,B., Shin,S., Kim,H.-J., Oh,T.-K., Kim,Y.-O., Choi,K. and Oh,B.-H. (2000) Nature Struct. Biol., 7, 147–153.[CrossRef][ISI][Medline]

Hasegawa,J., Uchiyama,S., Tanimoto,Y., Mizutani,M., Kobayashi,Y., Sambongi,Y. and Igarashi, Y. (2000) J. Biol. Chem., 275, 37824–37828.[Abstract/Free Full Text]

Honig,B. and Nicholls,A. (1995) Science, 268, 1144–1149.[ISI][Medline]

Juncosa, M., Pons, J., Planas, A. and Querol, E. (1994) BioTechniques, 16, 829–824.

Kanaya,S., Oobatake,M. and Liu,Y. (1996) J. Biol. Chem., 271, 32729–32736.[Abstract/Free Full Text]

Karshikoff, A. and Ladenstein, R. (1998) Protein Eng., 11, 867–872.[Abstract]

Klein,C., Georges,G., Kunkele,K.-P., Huber,R., Engh,R. and Hansen,S. (2001) J. Biol. Chem., 276, 37390–37401.[Abstract/Free Full Text]

Knochel,T., Pappenberger,A., Jansonius,J. and Kirschner,K. (2000). J. Biol. Chem., 277, 8626–8634.[Abstract/Free Full Text]

Kumar,S., Tsai,C.-J. and Nussinov,R. (2000) Protein Eng., 13, 179–191.[Abstract/Free Full Text]

Lebbnik,J., Knapp,S., van der Oost,J., Rice,D. and de Vos,W. (1999) J. Mol. Biol., 289, 357–369.[CrossRef][ISI][Medline]

Malakauskas,S. and Mayo,S. (1998) Nature Struct. Biol. 5, 470–475.[ISI][Medline]

Marshall,S., Morgan,C. and Mayo,S. (2000) J. Mol. Biol., 316, 189–199.[CrossRef]

Martin,A., Sieber,V. and Schmid,F. (2001) J. Mol. Biol., 309, 717–726.[CrossRef][ISI][Medline]

Meiering,E., Serrano,L. and Fersht,A. (1992) J. Mol. Biol., 225, 585–589.[ISI][Medline]

Mikami,B., Yoo,H.-J. and Yoshigi,N. (1999) J. Mol. Biol., 285, 1235–1243.[CrossRef][ISI][Medline]

Mozo-Villarias,A. and Querol,E. (2000) Curr. Top. Biochem. Res., 2, 19–31.

Nicholson,E., Mo,H., Prusiner,S., Cohen,F. and Merqusee,S. (2002) J. Mol. Biol., 316, 807–15.[CrossRef][ISI][Medline]

Ohmura,T., Ueda,T., Ootsuka,K., Saito,M. and Imoto,T. (2001) Protein Sci., 10, 313–320.[Abstract/Free Full Text]

Pedone,E., Saviano,M., Rossi,M. and Bartolucci,S. (2001) Protein Eng., 14, 255–260.[Abstract/Free Full Text]

Perl,D. and Schmid,F. (2001) J. Mol. Biol., 313, 343–357.[CrossRef][ISI][Medline]

Peterson,R., Nicholson,E., Tharpar,R., Klevit,R. and Scholtz,M. (1999) J. Mol. Biol., 286, 1609–1619.[CrossRef][ISI][Medline]

Petrey,D. and Honig,B. (2000) Protein Sci., 9, 2181–2191.[Abstract]

Pitera, J, Falta,M. and van Gunsteren,W. (2001) Biophys. J., 80, 2546–2555.[Abstract/Free Full Text]

Pokkuluri,P.R., Gu,M., Cai,X., Raffen,R., Stevens,F.J. and Schiffer,M. (2002) Protein Sci., 11, 1687–1694.[Abstract/Free Full Text]

Pons,J., Querol,E. and Planas,A. (1997) J. Biol. Chem., 272, 13006–13012.[Abstract/Free Full Text]

Predki,P., Agrawal,V., Brunger,A. and Regan,L. (1996) Nature Struct. Biol., 3, 54–58.[ISI][Medline]

Querol,E., Perez-Pons,J.A. and Mozo-Villarias,A. (1996) Protein Eng., 9, 265–271.[Abstract]

Robic,S., Berger,J. and Marquesee,S. (2002). Protein Sci., 11, 381–389.[Abstract/Free Full Text]

Schutz,C. and Warshel,A. (2001) Proteins: Struct. Funct. Genet., 44, 400–417.[CrossRef][ISI][Medline]

Shaw,K., Grimsley,G., Yakovlev,G., Makarov,A. and Pace,N. (2001) Protein Sci., 10, 1206–1215.[Abstract/Free Full Text]

Spassov,V., Karshikoff,A. and Ladenstein,R. (1994). Protein Sci., 3, 1556–1569.[Abstract/Free Full Text]

Spassov,V.Z., Ladenstein,R. and Karshikoff,A. (1997) Protein Sci., 6, 1190–1196.[Abstract/Free Full Text]

Sterner,R. and Liebl,W. (2001) Crit. Rev. Biochem. Mol. Biol., 36, 39–106.[Abstract/Free Full Text]

Stewart,R., Varghese,J., Garrett,T., Hog,P. and Fincher,G. (2001) Protein Eng., 14, 245–253.[Abstract/Free Full Text]

Sung,Y.-H., Hong,H.-D., Cheong,C., Kim,J., Cho,J., Kim,Y.R. and Lee,W. (2001) J. Biol. Chem., 276, 44229–44238.[Abstract/Free Full Text]

Tchan,M., Choy,K., Mackay,J., Lyons,A., Bains,N. and Weiss,A. (2000) J. Biol. Chem., 275, 37454–37461.[Abstract/Free Full Text]

Trejo,F., Gelpí,J., Ferrer,A., Boronat,A., Busquets,M. and Cortés,A. (2001) Protein Eng., 14, 911–917.[Abstract/Free Full Text]

Vieille,C. and Zeikus,G. (2001) Microbiol. Mol. Biol. Rev., 65, 1–43.[Abstract/Free Full Text]

Vogt,G., Woell,S. and Argos,P. (1997) J. Mol. Biol., 269, 631–643.[CrossRef][ISI][Medline]

White,P., Squirrell,D., Arnaud,P., Lowe,C. and Murray,J. (1996) Biochem. J., 319, 343–350.[ISI][Medline]

Wray,J., Baase,W., Lindstrom,J., Weaver,L., Poteete,A. and Matthews,B. (1999) J. Mol. Biol., 292, 1111–1120.[CrossRef][ISI][Medline]

Xiao,L. and Honig,B. (1999) J. Mol. Biol., 289, 1435–1444.[CrossRef][ISI][Medline]

Zhang,J., Liu,Z.-P., Jones,A., Gierasch,L. and Sambrook,J. (1992) Proteins: Struct. Funct. Genet., 13, 87–99.[ISI][Medline]

Zhu,G.-P., Xu,C.-X., Teng,M.-K., Tao,L.-M., Zhu,X.-Y., Wu,C.-J., Hang,J., Niu,L.-W. and Wang,Y.-Zh. (1999) Protein Eng., 12, 635–638.[Abstract/Free Full Text]

Received July 23, 2002; revised January 20, 2003; accepted February 5, 2003.