Extreme stabilization of a thermolysin-like protease by an engineered disulfide bond

Rational engineering of enzyme stability

During the past 15 years there has been a continuous flow of reports describing proteins stabilized by the introduction of mutations. These reports span a period from pioneering rational design work on small enzymes such as T4 lysozyme and barnase to protein design, and directed evolution. Concomitantly, the purification and characterization of naturally occurring hyperstable proteins has added to our understanding of protein stability. Along the way, many strategies for rational protein stabilization have been proposed, some of which (e.g. entropic stabilization by introduction of prolines or disulfide bridges) have reasonable success rates. On the other hand, comparative studies and efforts in directed evolution have revealed that there are many mutational strategies that lead to high stability, some of which are not easy to define and rationalize. Recent developments in the field include increasing awareness of the importance of the protein surface for stability, as well as the notion that normally a very limited number of mutations can yield a large increase in stability. Another development concerns the notion that there is a fundamental difference between the "laboratory stability" of small pure proteins that unfold reversibly and completely at high temperatures and "industrial stability", which is usually governed by partial unfolding processes followed by some kind of irreversible inactivation process (e.g. aggregation). Provided that one has sufficient knowledge of the mechanism of thermal inactivation, successful and efficient rational stabilization of enzymes can be achieved.

Do Proteins Always Benefit from a Stability Increase? Relevant and Residual Stabilisation in a Three-state Protein by Charge Optimisation

The vast majority of our knowledge on protein stability arises from the study of simple two-state models. However, proteins displaying equilibrium intermediates under certain conditions abound and it is unclear whether the energetics of native/intermediate equilibria is well represented in current knowledge. We consider here that the overall conformational stability of three-state proteins is made of a "relevant" term and a "residual" one, corresponding to the free energy differences of the native to intermediate (N-to-I) and intermediate to denatured (I-to-D) equilibria, respectively. The N-to-I free energy difference is considered to be the relevant stability because protein-unfolding intermediates are likely devoid of biological activity. We use surface charge optimisation to first increase the overall (N-to-D) stability of a model three-state protein (apoflavodoxin) and then investigate whether the stabilisation obtained is realised into relevant or into residual stability. Most of the mutations designed from electrostatic calculations or from simple sequence conservation analysis produce large increases in the overall stability of the protein. However, in most cases, this simply leads to similarly large increases of the residual stability. Two mutations, nevertheless, show a different trend and increase the relevant stability of the protein substantially. When all the mutations are mapped onto the structure of the apoflavodoxin thermal-unfolding intermediate (obtained independently by equilibrium phi-analysis and NMR) they cluster perfectly so that the mutations increasing the relevant stability appear in the small unstructured region of the intermediate and the others in the native-like region. This illustrates the need for specific investigation of N-to-I equilibria and the structure of protein intermediates, and indicates that it is possible to rationally stabilise a protein against partial unfolding once the structure of the intermediate conformation is known, even if at low resolution.

Structural Basis of Protein Kinetic Stability: Resistance to Sodium Dodecyl Sulfate Suggests a Central Role for Rigidity and a Bias Toward β-Sheet Structure

The term kinetic stability is used to describe proteins that are trapped in a specific conformation because of an unusually high-unfolding barrier that results in very slow unfolding rates. Motivated by the observation that some proteins are resistant to sodium dodecyl sulfate (SDS)-induced denaturation, an attempt was made to determine whether this property is a result of kinetic stability. We studied many proteins, including a few kinetically stable proteins known to be resistant to SDS. The resistance to SDS-induced denaturation was investigated by comparing the migration on polyacrylamide gels of identical boiled and unboiled protein samples containing SDS. On the basis of the different migration of these samples, eight proteins emerged as being resistant to SDS. The kinetic stability of these proteins was confirmed by their slow unfolding rate upon incubation in guanidine hydrochloride. Further studies showed that these proteins were also extremely resistant to proteolysis by proteinase K, suggesting that a common mechanism may account for their resistance to SDS and proteolytic cleavage. Together, these observations suggest that a rigid protein structure may be the physical basis for kinetic stability and that resistance to SDS may serve as a simple assay for identifying proteins whose native conformations are kinetically trapped. Remarkably, most of the kinetically stable SDS-resistant proteins in this study are oligomeric beta-sheet proteins, suggesting a bias of these types of structures toward kinetic stability.

Prediction of the bonding states of cysteines Using the support vector machines based on multiple feature vectors and cysteine state sequences

The support vector machine (SVM) method is used to predict the bonding states of cysteines. Besides using local descriptors such as the local sequences, we include global information, such as amino acid compositions and the patterns of the states of cysteines (bonded or nonbonded), or cysteine state sequences, of the proteins. We found that SVM based on local sequences or global amino acid compositions yielded similar prediction accuracies for the data set comprising 4136 cysteine-containing segments extracted from 969 nonhomologous proteins. However, the SVM method based on multiple feature vectors (combining local sequences and global amino acid compositions) significantly improves the prediction accuracy, from 80% to 86%. If coupled with cysteine state sequences, SVM based on multiple feature vectors yields 90% in overall prediction accuracy and a 0.77 Matthews correlation coefficient, around 10% and 22% higher than the corresponding values obtained by SVM based on local sequence information.

Crystal structures and structural stabilities of the disulfide bond-deficient soybean proglycinin mutants C12G and C88S

The constituent subunits of seed storage protein 11S globulin have two disulfide bonds that are common among 11S globulins from legume and nonlegume seeds. In the case of the A1aB1b subunit of soybean 11S globulin, glycinin, Cys12-Cys45 and Cys88-Cys298 are observed by X-ray crystallography. The significance of these two disulfide bonds for structural stability was investigated by mutagenesis of Cys12 to Gly and of Cys88 to Ser. The disulfide bond-deficient mutants C12G and C88S could form the correct conformations identical to that of the wild-type proglycinin except in the vicinities of the mutation sites C12 and C88 as shown by their crystal structures. Thermal stability monitored by differential scanning calorimetry of the mutants indicated that the contribution of these disulfide bonds to the thermal stability of proglycinin A1aB1b is low, although there is a small difference in the extent of the contribution between the two disulfide bonds (Cys12-Cys45 > Cys88-Cys298). The contribution of Cys88-Cys298 to the resistance of proglycinin A1aB1b to proteinase digestion is higher than that of Cys12-Cys45. Possible effects of structure on the different properties of C12G and C88S are discussed.

Kinetic analysis of enhanced thermal stability of an alkaline protease with engineered twin disulfide bridges and calcium-dependent stability

The thermal stability of a cysteine-free alkaline protease (Alp) secreted by the eukaryote Aspergillus oryzae was improved both by the introduction of engineered twin disulfide bridges (Cys-69/Cys-101 and Cys-169/Cys-200), newly constructed as part of this study, and by the addition of calcium ions. We performed an extensive kinetic analysis of the increased thermal stability of the mutants as well as the role of calcium dependence. The thermodynamic activation parameters for irreversible thermal inactivation, the activation free energy (deltaG), the activation enthalpy (deltaH), and the activation entropy (deltaS) were determined from absolute reaction rate theory. The values of deltaH and deltaS were significantly and concomitantly increased as a result of introducing the twin disulfide bridges, for which the increase in the value of deltaH outweighed that of deltaS, resulting in significant increases in the value of deltaG. The enhancement of the thermal stability obtained by introducing the twin disulfide bridges is an example of the so-called low-temperature stabilization of enzymes. The stabilizing effect of calcium ions on wild-type Alp is similar to the results we obtained by introducing the engineered twin disulfide bridges.

Stabilizing the integrin alpha M inserted domain in alternative conformations with a range of engineered disulfide bonds

Conformational movement of the C-terminal alpha7 helix in the integrin inserted (I) domain, a major ligand-binding domain that adopts an alpha/beta Rossmann fold, has been proposed to allosterically regulate ligand-binding activity. Disulfide bonds were engineered here to reversibly lock the position of the alpha7 helix in one of two alternative conformations seen in crystal structures, termed open and closed. Our results show that pairs of residues with Cbeta atoms farther apart than optimal for disulfide bond stereochemistry can be successfully replaced by cysteine, suggesting that backbone movement accommodates disulfide formation. We also find more success with substituting partially exposed than buried residues. Disulfides stabilizing the open conformation resulted in constitutively active alphaMbeta2 heterodimers and isolated alphaM inserted domains, which were reverted to an inactive form by dithiothreitol reduction. By contrast, a disulfide stabilizing the closed conformation resulted in inactive alphaMbeta2 that was resistant to activation but became activatable after dithiothreitol treatment.

The influence of intramolecular bridges on the dynamics of a protein folding reaction

Thirteen versions of a beta-sheet protein have been constructed, each with a single, surface-exposed disulfide bridge. A comparison of folding kinetics, in oxidizing and reducing conditions, is used to elucidate the order in which beta-strands become associated during the folding process and, hence, the relationship between topology and folding dynamics. In common with the wild-type molecule, all the proteins fold through a two-step (three state) mechanism with a rapidly formed intermediate which slowly converts to the native state. In a majority of cases, the bridge is seen to stabilize the folded state, and for five of the modified proteins, the additional stability is greater than 3 kcal/mol. Surprisingly, cross-links which connect beta-strands which are distant in sequence predominantly stabilize the rapidly formed intermediate state, suggesting that these strand-strand interactions occur in the initial stages of folding. Cross-links which stabilize local hairpins have their major influence on the second, rate-determining step leading to significant enhancements in the folding rate. We find that enhancement of the folding rate in the second, rate-limiting step is correlated with a reduction in contact order in the same way as in naturally occurring proteins of different folds. The large increases in native-state stability resulting from the insertion of disulfide bridges on the surface of beta-sheet structures have implications for enhancing the robustness of proteins by molecular engineering.

Selection of mutations for increased protein stability

There are many ways to select mutations that increase the stability of proteins, including rational design, functional screening of randomly generated mutant libraries, and comparison of naturally occurring homologous proteins. The protein engineer's toolbox is expanding and the number of successful examples of engineered protein stability is increasing. Still, the selection of thermostable mutations is not a standard process. Selection is complicated by lack of knowledge of the process that leads to thermal inactivation and by the fact that proteins employ a large variety of structural tricks to achieve stability.

Structural basis of thermostability. Analysis of stabilizing mutations in subtilisin BPN'

The crystal structures of two thermally stabilized subtilisin BPN' variants, S63 and S88, are reported here at 1.8 and 1.9 A resolution, respectively. The micromolar affinity calcium binding site (site A) has been deleted (Delta75-83) in these variants, enabling the activity and thermostability measurements in chelating conditions. Each of the variants includes mutations known previously to increase the thermostability of calcium-independent subtilisin in addition to new stabilizing mutations. S63 has eight amino acid replacements: D41A, M50F, A73L, Q206W, Y217K, N218S, S221C, and Q271E. S63 has 75-fold greater stability than wild type subtilisin in chelating conditions (10 mm EDTA). The other variant, S88, has ten site-specific changes: Q2K, S3C, P5S, K43N, M50F, A73L, Q206C, Y217K, N218S, and Q271E. The two new cysteines form a disulfide bond, and S88 has 1000 times greater stability than wild type subtilisin in chelating conditions. Comparisons of the two new crystal structures (S63 in space group P2(1) with A cell constants 41.2, 78.1, 36.7, and beta = 114.6 degrees and S88 in space group P2(1)2(1)2(1) with cell constants 54.2, 60.4, and 82.7) with previous structures of subtilisin BPN' reveal that the principal changes are in the N-terminal region. The structural bases of the stabilization effects of the new mutations Q2K, S3C, P5S, D41A, Q206C, and Q206W are generally apparent. The effects are attributed to the new disulfide cross-link and to improved hydrophobic packing, new hydrogen bonds, and other rearrangements in the N-terminal region.

The effects of modifying the surface charge on the catalytic activity of a thermolysin-like protease

The impact of long range electrostatic interactions on catalysis in the thermolysin-like protease from Bacillus stearothermophilus was studied by analyzing the effects of inserting or removing charges on the protein surface. Various mutations were introduced at six different positions, and double-mutant cycle analysis was used to study the extent to which mutational effects were interdependent. The effects of single point mutations on the k(cat)/K(m) were non-additive, even in cases where the point mutations were located 10 A or more from the active site Zn(2+) and separated from each other by up to 25 A. This shows that catalysis is affected by large electrostatic networks that involve major parts of the enzyme. The interdependence of mutations at positions as much as 25 A apart in space also indicates that other effects, such as active site dynamics, play an important role in determining active site electrostatics. Several mutations yielded a significant increase in the activity, the most active (quadruple) mutant being almost four times as active as the wild type. In some cases the shape of the pH-activity profile was changed significantly. Remarkably, large changes in the pH-optimum were not observed.

Molecular dynamics simulations as a tool for improving protein stability

Haloalkane dehalogenase (DhlA) was used as a model protein to explore the possibility to use molecular dynamics (MD) simulations as a tool to identify flexible regions in proteins that can serve as a target for stability enhancement by introduction of a disulfide bond. DhlA consists of two domains: an alpha/beta-hydrolase fold main domain and a cap domain composed of five alpha-helices. MD simulations of DhlA showed high mobility in a helix-loop-helix region in the cap domain, involving residues 184-211. A disulfide cross-link was engineered between residue 201 of this flexible region and residue 16 of the main domain. The mutant enzyme showed substantial changes in both thermal and urea denaturation. The oxidized form of the mutant enzyme showed an increase of the apparent transition temperature from 47.5 to 52.5 degrees C, whereas the T(m,app) of the reduced mutant decreased by more than 8 degrees C compared to the wild-type enzyme. Urea denaturation results showed a similar trend. Measurement of the kinetic stability showed that the introduction of the disulfide bond caused a decrease in activation free energy of unfolding of 0.43 kcal mol(-1) compared to the wild-type enzyme and also indicated that the helix-loop-helix region was involved early in the unfolding process. The results show that MD simulations are capable of identifying mobile protein domains that can successfully be used as a target for stability enhancement by the introduction of a disulfide cross-link.

Stabilization of a (betaalpha)8-barrel protein by an engineered disulfide bridge

The aim of this study was to increase the stability of the thermolabile (betaalpha)8-barrel enzyme indoleglycerol phosphate synthase from Escherichia coli by the introduction of disulfide bridges. For the design of such variants, we selected two out of 12 candidates, in which newly introduced cysteines potentially form optimal disulfide bonds. These variants avoid short-range connections, substitutions near catalytic residues, and crosslinks between the new and the three parental cysteines. The variant linking residues 3 and 189 fastens the N-terminus to the (betaalpha)8-barrel. The rate of thermal inactivation at 50 degrees C of this variant with a closed disulfide bridge is 65-fold slower than that of the reference dithiol form, but only 13-fold slower than that of the parental protein. The near-ultraviolet CD spectrum, the reactivity of parental buried cysteines with Ellman's reagent as well as the decreased turnover number indicate that the protein structure is rigidified. To confirm these data, we have solved the X-ray structure to 2.1-A resolution. The second variant was designed to crosslink the terminal modules betaalpha1 and betaalpha8. However, not even the dithiol form acquired the native fold, possibly because one of the targeted residues is solvent-inaccessible in the parental protein.

The effect of changing the hydrophobic S1' subsite of thermolysin-like proteases on substrate specificity

The hydrophobic S1' subsite is one of the major determinants of the substrate specificity of thermolysin and related M4 family proteases. In the thermolysin-like protease (TLP) produced by Bacillus stearothermophilus (TLP-ste), the hydrophobic S1' subsite is mainly formed by Phe130, Phe133, Val139 and Leu202. In the present study, we have examined the effects of replacing Leu202 by smaller (Gly, Ala, Val) and larger (Phe, Tyr) hydrophobic residues. The mutational effects showed that the wild-type S1' pocket is optimal for binding leucine side chains. Reduction of the size of residue 202 resulted in a higher efficiency towards substrates with Phe in the P1' position. Rather unexpectedly, the Leu202-->Phe and Leu202-->Tyr mutations, which were expected to decrease the size of the S1' subsite, resulted in a large increase in activity towards dipeptide substrates with Phe in the P1' position. This is probably due to the fact that 202Phe and 202Tyr adopt a second possible rotamer that opens up the subsite compared to Leu202, and also favours interactions with the substrate. To validate these results, we constructed variants of thermolysin with changes in the S1' subsite. Thermolysin and TLP-ste variants with identical S1' subsites were highly similar in terms of their preference for Phe vs. Leu in the P1' position.

Enhancement of the thermal stability of pyroglutamyl peptidase I by introduction of an intersubunit disulfide bond

From the comparison of the three-dimensional structure of mesophilic pyroglutamyl peptidase from Bacillus amyloliquefaciens and the thermophilic enzyme from Thermococcus litoralis, the intersubunit disulfide bond was estimated to be one of the factors for thermal stability. Since Ser185 was corresponded to Cys190 of the thermophilic enzyme by sequence alignment, the Ser185 residue was replaced with cysteine by site-directed mutagenesis. The S185C mutant enzyme appeared to form a disulfide bond, which was confirmed by SDS-PAGE with and without 2-mercaptoethanol. The mutant enzyme showed a catalytic efficiency equivalent to that of the wild-type enzyme for hydrolysis of a synthetic peptide substrate. However, the thermal stability of the S185C mutant was found to be 30 degrees C higher than that of wild-type. Thus the introduction of a disulfide bond enhanced thermal stability without changing the catalytic efficiency of the enzyme.

Differentiation between conformational and autoproteolytic stability of the neutral protease from Bacillus stearothermophilus containing an engineered disulfide bond

The introduction of a disulfide bond into the neutral protease from Bacillus stearothermophilus by the double mutation G8C/N60C had resulted in an extremely thermostable enzyme with a half-life of 35.9 min at 92.5 degrees C [Mansfeld, J., Vriend, G., Dijkstra, B.W., Veltman, O.R., van den Burg, B., Venema, G., Ulbrich-Hofmann, R. & Eijsink, V.G. (1997) J. Biol. Chem. 272, 11152-11156]. The study in guanidine hydrochloride of this enzyme and the respective wild-type enzyme allowed us to distinguish between the stability toward global unfolding and autoproteolysis. At low protease concentrations (20 microg.mL-1) and short periods of incubation with guanidine hydrochloride (5 min), transition curves without the interference by autoproteolysis could be derived from fluorescence emission measurements. The effect of the disulfide bond on the global unfolding of the protein proved to be smaller than expected. In contrast, the measurement of autoproteolysis at higher protein concentrations (100 microg.mL-1) by quantitative evaluation of the bands of intact protein on SDS/PAGE revealed a strong stabilization toward autoproteolytic degradation by the disulfide bond. The rate of autoproteolysis in guanidine hydrochloride was found to be much lower than that of thermal denaturation, which can be attributed to the inhibition of the proteases by this denaturant. The results suggest that the disulfide bond stabilizes the protease against autoproteolysis more than against global unfolding. Autoproteolysis starts as soon as the cleavage sites in flexible external structural regions become accessible. It is suggested that the stabilizing effect of the disulfide bond is caused by the fixation of the crucial loop region 56-69 or by hindrance of the primary cleavage in this region by the amino acid exchanges.

Reversibly locking a protein fold in an active conformation with a disulfide bond: integrin alphaL I domains with high affinity and antagonist activity in vivo

The integrin alphaLbeta2 has three different domains in its headpiece that have been suggested to either bind ligand or to regulate ligand binding. One of these, the inserted or I domain, has a fold similar to that of small G proteins. The I domain of the alphaM and alpha2 subunits has been crystallized in both open and closed conformations; however, the alphaL I domain has been crystallized in only the closed conformation. We hypothesized that the alphaL domain also would have an open conformation, and that this would be the ligand binding conformation. Therefore, we introduced pairs of cysteine residues to form disulfides that would lock the alphaL I domain in either the open or closed conformation. Locking the I domain open resulted in a 9,000-fold increase in affinity to intercellular adhesion molecule-1 (ICAM-1), which was reversed by disulfide reduction. By contrast, the affinity of the locked closed conformer was similar to wild type. Binding completely depended on Mg(2+). Orders of affinity were ICAM-1 > ICAM-2 > ICAM-3. The k(on), k(off), and K(D) values for the locked open I domain were within 1.5-fold of values previously determined for the alphaLbeta2 complex, showing that the I domain is sufficient for full affinity binding to ICAM-1. The locked open I domain antagonized alphaLbeta2-dependent adhesion in vitro, lymphocyte homing in vivo, and firm adhesion but not rolling on high endothelial venules. The ability to reversibly lock a protein fold in an active conformation with dramatically increased affinity opens vistas in therapeutics and proteomics.

Hyperthermophilic enzymes: sources, uses, and molecular mechanisms for thermostability

Enzymes synthesized by hyperthermophiles (bacteria and archaea with optimal growth temperatures of > 80 degrees C), also called hyperthermophilic enzymes, are typically thermostable (i.e., resistant to irreversible inactivation at high temperatures) and are optimally active at high temperatures. These enzymes share the same catalytic mechanisms with their mesophilic counterparts. When cloned and expressed in mesophilic hosts, hyperthermophilic enzymes usually retain their thermal properties, indicating that these properties are genetically encoded. Sequence alignments, amino acid content comparisons, crystal structure comparisons, and mutagenesis experiments indicate that hyperthermophilic enzymes are, indeed, very similar to their mesophilic homologues. No single mechanism is responsible for the remarkable stability of hyperthermophilic enzymes. Increased thermostability must be found, instead, in a small number of highly specific alterations that often do not obey any obvious traffic rules. After briefly discussing the diversity of hyperthermophilic organisms, this review concentrates on the remarkable thermostability of their enzymes. The biochemical and molecular properties of hyperthermophilic enzymes are described. Mechanisms responsible for protein inactivation are reviewed. The molecular mechanisms involved in protein thermostabilization are discussed, including ion pairs, hydrogen bonds, hydrophobic interactions, disulfide bridges, packing, decrease of the entropy of unfolding, and intersubunit interactions. Finally, current uses and potential applications of thermophilic and hyperthermophilic enzymes as research reagents and as catalysts for industrial processes are described.

Substrate specificity in the highly heterogeneous M4 peptidase family is determined by a small subset of amino acids

The members of the M4 peptidase family are involved in processes as diverse as pathogenicity and industrial applications. For the first time a number of M4 family members, also known as thermolysin-like proteases, has been characterized with an identical substrate set and a uniform set of assay conditions. Characterization with peptide substrates as well as high performance liquid chromatography analysis of beta-casein digests shows that the M4 family is a homogeneous family in terms of catalysis, even though there is a significant degree of amino acid sequence variation. The results of this study show that differences in substrate specificity within the M4 family do not correlate with overall sequence differences but depend on a small number of identifiable amino acids. Indeed, molecular modeling followed by site-directed mutagenesis of one of the substrate binding pocket residues of the thermolysin-like proteases of Bacillus stearothermophilus converted the catalytic characteristics of this variant into that of thermolysin.

Signal peptide-dependent protein transport in Bacillus subtilis: a genome-based survey of the secretome

One of the most salient features of Bacillus subtilis and related bacilli is their natural capacity to secrete a variety of proteins into their environment, frequently to high concentrations. This has led to the commercial exploitation of bacilli as major "cell factories" for secreted enzymes. The recent sequencing of the genome of B. subtilis has provided major new impulse for analysis of the molecular mechanisms underlying protein secretion by this organism. Most importantly, the genome sequence has allowed predictions about the composition of the secretome, which includes both the pathways for protein transport and the secreted proteins. The present survey of the secretome describes four distinct pathways for protein export from the cytoplasm and approximately 300 proteins with the potential to be exported. By far the largest number of exported proteins are predicted to follow the major "Sec" pathway for protein secretion. In contrast, the twin-arginine translocation "Tat" pathway, a type IV prepilin-like export pathway for competence development, and ATP-binding cassette transporters can be regarded as "special-purpose" pathways, through which only a few proteins are transported. The properties of distinct classes of amino-terminal signal peptides, directing proteins into the various protein transport pathways, as well as the major components of each pathway are discussed. The predictions and comparisons in this review pinpoint important differences as well as similarities between protein transport systems in B. subtilis and other well-studied organisms, such as Escherichia coli and the yeast Saccharomyces cerevisiae. Thus, they may serve as a lead for future research and applications.

What can disulfide bonds tell us about protein energetics, function and folding: simulations and bioninformatics analysis

We study the impact of disulfide bonds on protein stability and folding. Using lattice model simulations, we show that formation of a disulfide bond stabilizes a protein to an extent that depends on the distance along the chain between linked cysteine residues. However, the impact of disulfide bonds on folding kinetics varies broadly, from acceleration when disulfides are introduced in or close to the folding nucleus, to slowing when disulfides are introduced outside the nucleus. Having established the effect of disulfide bonds on stability, we study the correlation between the number of disulfide bonds and the composition of certain amino acid classes with the goal to use it as a statistical probe into factors that contribute to stability of proteins. We find that the number of disulfides is negatively correlated with aliphatic hydrophobic but not aromatic content. It is surprising that we observe a strong correlation of disulfide content with polar (Q,S,T,N) amino acid content and a strong negative correlation with charged (E,D,K,R) content. These findings provide insights into factors that determine protein stability and principles of protein design as well as possible relations of disulfide bonds and protein function.

Novel disulfide engineering in human carbonic anhydrase II using the PAIRWISE side-chain geometry database

An analysis of the pairwise side-chain packing geometries of cysteine residues observed in high-resolution protein crystal structures indicates that cysteine pairs have pronounced orientational preferences due to the geometric constraints of disulfide bond formation. A potential function was generated from these observations and used to evaluate models for novel disulfide bonds in human carbonic anhydrase II (HCAII). Three double-cysteine variants of HCAII were purified and the effective concentrations of their thiol groups were determined by titrations with glutathione and dithiothreitol. The effects of the cysteine mutations on the native state structure and stability were characterized by circular dichroism, enzymatic activity, sulfonamide binding, and guanidine hydrochloride titration. These analyses indicate that the PAIRWISE potential is a good predictor of the strength of the disulfide bond itself, but the overall structural and thermodynamic effects on the protein are complicated by additional factors. In particular, the effects of cysteine substitutions on the native state and the stabilization of compact nonnative states by the disulfide can override any stabilizing effect of the cross-link.

Single-point amino acid substitutions at the 119th residue of thermolysin and their pressure-induced activation

The effect of amino acid substitution at the 119th site of thermolysin (TLN) on the pressure activation behavior of this enzyme was studied for four mutants at pressures < 300 MPa. For Q119Q, Q119N and Q119R, the highest activation was observed to be over 30 times that at atmospheric pressure and the activation volumes (deltaV++) were about -75 ml/mol. However, we obtained only 10 times higher activation for Q119E and Q119D (deltaV++ approximately -60 ml/mol). The intrinsic fluorescence of TLN changed at pressures > 300 MPa, and the latter two mutants showed a smaller deltaGapp and deltaVapp of transition than the wild type. These results are discussed with respect to the hydration change in the enzyme protein around the substituted region.

DNA shuffling of subgenomic sequences of subtilisin

DNA family shuffling of 26 protease genes was used to create a library of chimeric proteases that was screened for four distinct enzymatic properties. Multiple clones were identified that were significantly improved over any of the parental enzymes for each individual property. Family shuffling, also known as molecular breeding, efficiently created all of the combinations of parental properties, producing a great diversity of property combinations in the progeny enzymes. Thus, molecular breeding, like classical breeding, is a powerful tool for recombining existing diversity to tailor biological systems for multiple functional parameters.

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.

Probing the unfolding region in a thermolysin-like protease by site-specific immobilization

Protein stabilization by immobilization has been proposed to be most effective if the protein is attached to the carrier at that region where unfolding is initiated. To probe this hypothesis, we have studied the effects of site-specific immobilization on the thermal stability of mutants of the thermolysin-like protease from Bacillus stearothermophilus (TLP-ste). This enzyme was chosen because previous studies had revealed which parts of the molecule are likely to be involved in the early steps of thermal unfolding. Cysteine residues were introduced by site-directed mutagenesis into various positions of a cysteine-free variant of TLP-ste. The mutant enzymes were immobilized in a site-specific manner onto Activated Thiol-Sepharose. Two mutants (T56C, S65C) having their cysteine in the proposed unfolding region of TLP-ste showed a 9- and 12-fold increase in half-lives at 75 degrees C due to immobilization. The stabilization by immobilization was even larger (33-fold) for the T56C/S65C double mutant enzyme. In contrast, mutants containing cysteines in other parts of the TLP-ste molecule (N181C, S218C, T299C) showed only small increases in half-lives due to immobilization (maximum 2.5-fold). Thus, the stabilization obtained by immobilization was strongly dependent on the site of attachment. It was largest when TLP-ste was fixed to the carrier through its postulated unfolding region. The concept of the unfolding region may be of general use for the design of strategies to stabilize proteins.

Heat stability of the potato tuber ADP-glucose pyrophosphorylase: role of Cys residue 12 in the small subunit

Most of the ADP-glucose pyrophosphorylases from different sources are stable to a heat treatment. We found that in the potato (Solanum tuberosum L.) tuber enzyme, the intermolecular disulfide bridge located between Cys12 of the small subunits is responsible for the stability at 60 degrees C. When this unique disulfide bond is cleaved the enzyme is stable up to 40 degrees C. Mutation of Cys12 in the small subunit into either Ala or Ser yielded enzymes with stability similar to the reduced form of the wild type. Concurrently, the enzyme with a truncated small subunit on the N-terminal was stable only up to 40 degrees C. Thus, the N-terminal is important for the stability of the enzyme because of the presence of a disulfide bond.

Engineering out motion: introduction of a de novo disulfide bond and a salt bridge designed to close a dynamic cleft on the surface of cytochrome b5

A previous molecular dynamics (MD) simulation of cytochrome b5 (cyt b5) at 25 degrees C displayed localized dynamics on the surface of the protein giving rise to the periodic formation of a cleft that provides access to the heme through a protected hydrophobic channel [Storch and Daggett (1995) Biochemistry 34, 9682]. Here we describe the production and testing of mutants designed to prevent the cleft from opening using a combination of experimental and theoretical techniques. Two mutants have been designed to close the surface cleft: S18D to introduce a salt bridge and S18C:R47C to incorporate a disulfide bond. The putative cleft forms between two separate cores of the protein: one is structural in nature and can be monitored through the fluorescence of Trp 22, and the other binds the heme prosthetic group and can be tracked via heme absorbance. An increase in motion localized to the cleft region was observed for each protein, except for the disulfide-containing variant, in MD simulations at 50 degrees C compared to simulations at 25 degrees C. For the disulfide-containing variant, the cleft remained closed. Both urea and temperature denaturation curves were nearly identical for wild-type and mutant proteins when heme absorbance was monitored. In contrast, fluorescence studies revealed oxidized S18C:R47C to be considerably more stable based on the midpoints of the denaturation transitions, Tm and U1/2. Moreover, the fluorescence changes for each protein were complete at approximately 50 degrees C and a urea concentration of approximately 3.9 M, significantly below the temperature and urea concentration (62 degrees C, 5 M urea) required to observe heme release. In addition, solvent accessibility based on acrylamide quenching of Trp 22 was lower in the S18C:R47C mutant, particularly at 50 degrees C, before heme release [presented in the accompanying paper (58)]. The results suggest that a constraining disulfide bond can be designed to inhibit dynamic cleft formation on the surface of cyt b5. Located near the heme, the native dynamics of the cleft may be functionally important for protein-protein recognition and/or complex stabilization.

Structural and mutagenesis studies of leishmania triosephosphate isomerase: a point mutation can convert a mesophilic enzyme into a superstable enzyme without losing catalytic power

The dimeric enzyme triosephosphate isomerase (TIM) has a very tight and rigid dimer interface. At this interface a critical hydrogen bond is formed between the main chain oxygen atom of the catalytic residue Lys13 and the completely buried side chain of Gln65 (of the same subunit). The sequence of Leishmania mexicana TIM, closely related to Trypanosoma brucei TIM (68% sequence identity), shows that this highly conserved glutamine has been replaced by a glutamate. Therefore, the 1.8 A crystal structure of leishmania TIM (at pH 5.9) was determined. The comparison with the structure of trypanosomal TIM shows no rearrangements in the vicinity of Glu65, suggesting that its side chain is protonated and is hydrogen bonded to the main chain oxygen of Lys13. Ionization of this glutamic acid side chain causes a pH-dependent decrease in the thermal stability of leishmania TIM. The presence of this glutamate, also in its protonated state, disrupts to some extent the conserved hydrogen bond network, as seen in all other TIMs. Restoration of the hydrogen bonding network by its mutation to glutamine in the E65Q variant of leishmania TIM results in much higher stability; for example, at pH 7, the apparent melting temperature increases by 26 degrees C (57 degrees C for leishmania TIM to 83 degrees C for the E65Q variant). This mutation does not affect the kinetic properties, showing that even point mutations can convert a mesophilic enzyme into a superstable enzyme without losing catalytic power at the mesophilic temperature.

Early steps in the unfolding of thermolysin-like proteases

Several series of site-directed mutations in thermolysin-like proteases are presented that show remarkable nonadditivity in their effect on thermal stability. A simple model is proposed that relates this nonadditivity to the occurrence of independent partial unfolding processes that occur in parallel at elevated temperatures. To prove this model, a thermolysin-like protease was designed in which two mutations located approximately 35 A apart in the structure individually exert small stabilizing effects of 2.3 and 4. 1 degreesC, respectively, but when combined stabilize the protease by 14.6 degreesC. This overadditivity, which follows directly from the model, confirms that unfolding of this engineered protease starts in parallel at two different regions of the protein.

Probing catalytic hinge bending motions in thermolysin-like proteases by glycine --> alanine mutations

The active site of thermolysin-like proteases (TLPs) is located at the bottom of a cleft between the N- and C-terminal domains. Crystallographic studies have shown that the active-site cleft is more closed in ligand-binding TLPs than in ligand-free TLPs. Accordingly, it has been proposed that TLPs undergo a hinge-bending motion during catalysis resulting in "closure" and "opening" of the active-site cleft. Two hinge regions have been proposed. One is located around a conserved glycine 78; the second involves residues 135 and 136. The importance of conserved glycine residues in these hinge regions was studied experimentally by analyzing the effects of Gly --> Ala mutations on catalytic activity. Eight such mutations were made in the TLP of Bacillus stearothermophilus (TLP-ste) and their effects on activity toward casein and various peptide substrates were determined. Only the Gly78Ala, Gly136Ala, and Gly135Ala + Gly136Ala mutants decreased catalytic activity significantly. These mutants displayed a reduction in kcat/Km for 3-(2-furylacryloyl)-L-glycyl-L-leucine amide of 73%, 62%, and 96%, respectively. Comparisons of effects on kcat/Km for various substrates with effects on the Ki for phosphoramidon suggested that the mutation at position 78 primarily had an effect on substrate binding, whereas the mutations at positions 135 and 136 primarily influence kcat. The apparent importance of conserved glycine residues in proposed hinge-bending regions for TLP activity supports the idea that hinge-bending is an essential part of catalysis.

A single calcium binding site is crucial for the calcium-dependent thermal stability of thermolysin-like proteases

Thermostable thermolysin-like proteases (TLPs), such as the TLP of Bacillus stearothermophilus CU-21 (TLP-ste), bind calcium in one double (Ca1,2) and two single (Ca3, Ca4) calcium binding sites. The single sites are absent in thermolabile TLPs, suggesting that they are determinants of (variation in) TLP stability. Mutations in the Ca3 and Ca4 sites of TLP-ste indeed reduced thermal stability, but only mutations in the Ca3 site affected the calcium-dependence of stability. The predominant effect of the Ca3 site results from the fact that the Ca3 site is part of a region of TLP-ste, which unfolding is crucial for thermal inactivation. Thermal inactivation is not caused by the absence of calcium from the Ca3 site per se, but rather by unfolding of a region of TLP-ste for which stability depends on the occupancy of the Ca3 site. In accordance with this concept is the observation that the effects of mutations in the Ca3 site could be compensated by stabilizing mutations near this site. In addition, it was observed that the contribution of calcium binding to the Ca3 was substantially reduced in extremely stable TLP-ste variants containing multiple stabilizing mutations in the Ca3 region. Apparently, in these latter variants, unfolding of the Ca3 region contributes little to the overall process of thermal inactivation.

Engineering an enzyme to resist boiling

In recent years, many efforts have been made to isolate enzymes from extremophilic organisms in the hope to unravel the structural basis for hyperstability and to obtain hyperstable biocatalysts. Here we show how a moderately stable enzyme (a thermolysin-like protease from Bacillus stearothermophilus, TLP-ste) can be made hyperstable by a limited number of mutations. The mutational strategy included replacing residues in TLP-ste by residues found at equivalent positions in naturally occurring, more thermostable variants, as well as rationally designed mutations. Thus, an extremely stable 8-fold mutant enzyme was obtained that was able to function at 100 degrees C and in the presence of denaturing agents. This 8-fold mutant contained a relatively large number of mutations whose stabilizing effect is generally considered to result from a reduction of the entropy of the unfolded state ("rigidifying" mutations such as Gly --> Ala, Ala --> Pro, and the introduction of a disulfide bridge). Remarkably, whereas hyperstable enzymes isolated from natural sources often have reduced activity at low temperatures, the 8-fold mutant displayed wild-type-like activity at 37 degrees C.

Mutational analysis of a surface area that is critical for the thermal stability of thermolysin-like proteases

Site-directed mutagenesis was used to assess the contribution of individual residues and a bound calcium in the 55-69 region of the thermolysin-like protease of Bacillus stearothermophilus (TLP-ste) to thermal stability. The importance of the 55-69 region was reflected by finding that almost all mutations had drastic effects on stability. These effects (both stabilizing and destabilizing) were obtained by mutations affecting main chain flexibility, as well as by mutations affecting the interaction between the 55-69 region and the rest of the protease molecule. The calcium-dependency of stability could be largely abolished by mutating one of its ligands (Asp57 or Asp59). In the case of the Asp57-->Ser mutation, the accompanying loss in stability was modest compared with the effects of other destabilizing mutations or the effects of (combinations of) stabilizing mutations. The detailed knowledge of the stability-determining region of TLP-ste permits effective rational design of stabilizing mutations, which, presumably, are also useful for related TLP such as thermolysin. This is demonstrated by the successful design of a stabilizing salt bridge involving residues 65 and 11.