Thermostabilization by replacement of specific residues with lysine in a Bacillus alkaline cellulase: building a structural model and implications of newly formed double intrahelical salt bridges

Tadahiro Ozawa1, Yoshihiro Hakamada1, Yuji Hatada1, Tohru Kobayashi1, Tsuyoshi Shirai2 and Susumu Ito1,3

1 Tochigi Research Laboratories, Kao Corporation, 2606 Akabane, Ichikai, Haga, Tochigi 321-3497 and 2 Department of Biotechnology and Biomaterial Chemistry, Graduate School of Engineering, Nagoya University, Chikusa-ku, Nagoya 464-8603, Japan


    Abstract
 Top
 Abstract
 Introduction
 Materials and methods
 Results and discussion
 References
 
An alkaline, mesophilic endo-1,4-ß-glucanase from alkaliphilic Bacillus sp. strain KSM-64 was significantly thermostabilized by replacement of both Asn179 and Asp194 with lysine by site-directed mutagenesis. Structural remodeling of the mutant enzyme newly generated by the double mutation suggested that Glu175->Lys179 and Glu190->Lys194 were the most plausible ion pairs, both of which involved side chains at the i and i + 4 positions on the {alpha}4-helix from Glu175 to Ser195. By molecular dynamics simulations, the N{zeta} hydrogens of Lys179 and Lys194 were found to coordinate with the carbonyl O{varepsilon}1 and O{varepsilon}2 of Glu175 and the carbonyl O{varepsilon}1 of Glu190, respectively, with distances of around 2 Å for all. These results confirm that the formation of these double intrahelical ion pairs (salt bridges) is responsible for the thermostabilization by the double mutation.

Keywords: Bacillus/cellulase/intrahelix/ion pair/salt bridge/thermostability


    Introduction
 Top
 Abstract
 Introduction
 Materials and methods
 Results and discussion
 References
 
Alkaline endo-1,4-ß-glucanases (Egls) (EC 3.2.1.4) from alkaliphilic Bacilli have been studied extensively for use as effective detergent additives (Hoshino and Ito, 1997Go; Ito, 1997Go; Horikoshi, 1999Go), including Egl-K for detergent powders (Yoshimatsu et al., 1990Go) and Egl-64 for liquid detergents (Shikata et al., 1990Go). A mutant Egl with higher thermostability is desirable because the most serious problem with detergent enzymes is thermal inactivation during isolation, granulation and storage. As is the case for detergent enzymes (Ito et al., 1998Go), one of the most pursued goals in protein engineering and biotechnology is the design of mutant proteins with increased thermostability. Many attempts have been made to relate the properties of enhanced thermostability to features inherent in the amino acid sequences and 3D structures of mesophilic and thermophilic enzymes (Vogt et al., 1997Go; Wallon et al., 1997Go; Závodszky et al., 1998Go; Declerck et al., 2000Go). The enhancement of thermostability of certain proteins is brought about by increased internal hydrophobicity, increased polar surface area, shortening or deletion of surface loops, improved helix capping, additional ion pairs, increased occurrence of Pro residues, decreased helix-destabilizing residues and so on (e.g. Argos et al., 1979Go; Imanaka et al., 1986Go; Serrano and Fersht, 1989Go; Davies et al., 1993Go; Takano et al., 1995Go; Haney et al., 1997Go; Russell et al., 1997Go; Shirai et al., 1997bGo; Watanabe et al., 1997Go). It is very informative that a net increase in ionic interactions is often observed in thermophilic enzymes (Vogt et al., 1997Go; Wallon et al., 1997Go; Hashimoto et al., 1999Go; Kumar et al., 2000Go).

Recently, we found a thermostable, alkaline Egl (Egl-237) from the alkaliphilic Bacillus isolate KSM-S237 (Hakamada et al., 1997Go) and sequenced the cloned gene of the enzyme (Hakamada et al., 2000Go). The deduced amino acid sequence of Egl-237 showed high homology with those of Egl-K (Ozaki et al., 1990Go) and Egl-64 (Sumitomo et al., 1992Go). Based on the difference of a few amino acids, we constructed several chimeric genes from Egl-64 and Egl-237, and Lys179 and Lys194 residues in Egl-237, not conserved in Egl-64, were found to contribute to thermostability after site-directed mutagenesis. As a result, the double mutation Asn179->Lys (NK)/Asp194->Lys (DK) (NK/DK) significantly improved thermostability, the level of which was as high as that of the wild-type Egl-237 (Hakamada et al., 2001Go). One of the most striking features of thermostable enzymes, when compared with mesophilic enzymes, is the decrease in the number of Lys residues, which are mainly replaced with Arg and Glu residues (Argos et al., 1979Go; Davies et al., 1993Go; Shirai et al., 1997aGo). However, the thermal stabilization of mesophilic Egl-64 by the additional Lys residues instead of Arg residues is clearly contrary to these results with thermophilic enzymes. In spite of the gained basic Lys residue(s), the pI value of the mutant Egl-64 enzymes estimated by electrofocusing was slightly higher than or similar to that of the wild-type level. This suggests that the thermal stabilization achieved with the mutations involves the formation of new ion pairs, such as salt bridges, with some carboxylate residues localized on the enzyme surfaces. Although the 3D structures of glycosyl hydrolase family 5 cellulases have been determined (Ducros et al., 1995Go; Sakon et al., 1996Go; Davies et al., 1998Go), their amino acid sequences and catalytic properties were too different from those of Egl-64 and Egl-237 in the same family to analyze the roles of the thermostabilizing Lys residues. Recently, we succeeded in crystallizing (Shirai et al., 1997bGo) and determining the 3D structure of the truncated Egl-K, designated Egl-Kt (Ozaki et al., 1995Go), similar to the two alkaline Bacillus Egls. Here we show the precise involvement of the specific Lys residues in and the mechanisms for the thermal stabilization of Egl-64 by homology modeling and implicating by molecular dynamics (MD) simulation.


    Materials and methods
 Top
 Abstract
 Introduction
 Materials and methods
 Results and discussion
 References
 
Bacterial strains and enzyme production

The alkaliphilic Bacillus sp. strain KSM-64 of our stock cultures was used as the source of the gene encoding Egl-64 (Sumitomo et al., 1992Go). Constructed plasmids were used individually for transformation of B.subtilis ISW1214 to exoproduce wild-type, NK, DK or NK/DK mutant enzyme in the medium described previously (Sumitomo et al., 1995Go; Hakamada et al., 2001Go).

Purification and assay of enzymes

The recombinant enzymes produced by B.subtilis were purified by chromatography on a DEAE-Toyopearl 650C column to homogeneity as judged by SDS–PAGE. Egl activity was measured at 40°C in 0.1 M glycine–NaOH buffer (pH 9.0) with 1.0% CMC as substrate by the dinitrosalicylic acid procedure as described (Hakamada et al., 2000Go).

Homology modeling of tertiary structures of Egl

The structural models of Egl-64 and its engineered protein (NK/DK) were constructed by the homology modeling method, based on the 3D structure of Egl-Kt (PDB code No. 1go1) and the deduced amino acid sequence of Egl-64 (Sumitomo et al., 1992Go). All data sets were processed on a Silicon Graphics Indigo2 workstation using the InsightII/Discover software package (Molecular Simulation). MD simulations were calculated in an AMBER forcefield, where a dielectric constant {varepsilon} of 1 and the partial atomic charges computed at pH 7 were used. The MD simulations were performed in 4000 fs after linear heating from 0 to 330 K within 1000 fs and the MD trajectory was saved at 20 fs intervals.

Gene accession numbers

The original nucleotide sequence data of Egl-64 have been deposited in the DDBJ, EMBL and GenBank databases under the accession No. M84963.


    Results and discussion
 Top
 Abstract
 Introduction
 Materials and methods
 Results and discussion
 References
 
Cumulative thermostabilization by mutations NK and DK

The thermostability of the NK, DK and NK/DK mutant enzymes was examined by heating each at various temperatures for 30 min in 0.1 M glycine–NaOH buffer (pH 9.0). As shown in Figure 1Go, the replacement of Asp194 with lysine (DK) thermostabilized Egl-64 significantly at high temperatures, as described (Hakamada et al., 2001Go). The stabilizing effect of the double mutation NK/DK was reproducibly cumulative, although the single NK mutation improved the thermostability only marginally. At 70°C, the half-lives (t1/2) due to thermal inactivation increased in the order wild-type Egl-64 (9 min), NK mutant (15 min), DK mutant (38 min) and NK/DK mutant (41 min).



View larger version (23K):
[in this window]
[in a new window]
 
Fig. 1. Effect of temperature on the stability of wild-type Egl-64 and its mutant enzymes. Each enzyme had previously been incubated for 30 min at the indicated temperatures in 0.1 M glycine–NaOH buffer (pH 9.0). The residual activities were measured under the standard conditions of enzyme assay. The original activity of each enzyme was taken as 100%. {circ}, Wild-type Egl-64; {blacktriangleup}, NK mutant; •, DK mutant; {blacksquare}, NK/DK mutant.

 
Structure-based alignment of amino acid sequences of Egl-64 and Egl-Kt

To allow accurate homology modeling and to understand better the mutagenesis results, the amino acid sequence of Egl-64 was aligned with that of Egl-Kt, as shown in Figure 2Go. High homology was observed within the suitably aligned sequences of both enzymes (73.6% identity and 94.4% similarity). The secondary structural elements, {alpha}-helices and ß-strands, including catalytic residues (Hakamada et al. 2000Go), in the 3D structure of Egl-Kt are included in Figure 2Go.



View larger version (64K):
[in this window]
[in a new window]
 
Fig. 2. Structure-based amino acid sequence alignment of Egl-Kt and Egl-64. Sequences shown and residue numbers correspond to the mature proteins of Egl-K and Egl-64. The N- and C-termini of Egl-Kt correspond to Ala228 and Leu584, respectively, of mature Egl-K (Ozaki et al., 1995Go). Amino acids are shown as the single-letter codes and asterisks above the sequences indicate the conserved residues. Open circles indicate the amino acid residues of Egl-64, replaced with lysine. Solid circles show the catalytic Glu residues (Hakamada et al. 2000Go). Underlines denote {alpha}-helices and double-underlines represent ß-strands.

 
Homology model building

To clarify the mechanism of thermostabilization of Egl-64 by the NK/DK mutation, we built a model of the enzyme, based on the 3D structure of Egl-Kt (PDB code No. 1go1), as shown in Figure 3Go. Before building the model, the amino acid sequences of both enzymes were appropriately aligned (Figure 2Go). The validity of the modeled structure was confirmed using the 3D-1D, XPLORE and PROCHECK programs. Owing to the high sequence homology, the modeled structure of Egl-64 (Figure 3BGo) closely resembles the 3D structure of Egl-Kt (Figure 3AGo). Essentially, the model of Egl-64 contains an {alpha}/ß-barrel core, as in the cases of the 3D structures of neutral Egls belonging to the glycosyl hydrolase family 5 (Ducros et al., 1995Go; Sakon et al., 1996Go; Davies et al., 1998Go).



View larger version (44K):
[in this window]
[in a new window]
 
Fig. 3. Schematic top views of (A) Egl-Kt and (B) Egl-64. Coordinates of Egl-Kt were obtained from the PDB using the accession number 1go1. The solid ribbon diagrams were drawn using the program WebLab Viewer. Helices and strands in the core barrel structure are colored red and cyan, respectively. The catalytic Glu residues are shown in the CPK representation.

 
Generation of intrahelical salt bridges by double mutation

In the modeled structure of the NK/DK mutant enzyme, the replacements of Lys residues at positions 179 and 194 form ion pairs with the side chains of Glu175 and Glu190, respectively, yielding double intrahelical ion pairs, as shown in Figure 4Go. Lys194 appears to interact also with Asp224 occurring on a loop between Asn220 and His226.



View larger version (19K):
[in this window]
[in a new window]
 
Fig. 4. Coordination of the Lys residues at positions (A) 179 and (B) 194 with the surrounding carboxylic acid residues in (A) NK and (B) DKmutant enzymes, with oxygen atoms in red, hydrogen atoms in gray, nitrogen atoms in purple and main chains in green. Broken lines represent hydrogen bonds.

 
To examine the difference in the ionic interactions, the ion pair-forming Lys residues were MD simulated. Glu175, Glu176, Glu183 and Asp219 surrounding Lys179 and Glu137, Glu190, Asp199 and Asp224 surrounding Lys194 were included in the simulations. In the results, the N{zeta} hydrogens on Lys179 interacted with the carbonyl O{varepsilon}1 and O{varepsilon}2 on Glu175, as shown in Figure 5AGo. In the MD trajectory, a minimum potential energy was obtained after 5000 fs (5000th step). In the case of the N{zeta} hydrogen on Lys194, the distance to the carbonyl O{varepsilon}1 on Glu190 was around 2 Å, whereas the distance to the carbonyl O{delta}2 on Asp224 was small but frequently >3 Å (Figure 5BGo). Although the computed distances are shorter than expected, these results indicate that the most plausible ion pairs between Lys179 and Glu175 and between Lys194 and Glu190 are generated by the NK/DK mutation. Each pair involves side chains at the i and i + 4 positions of the single {alpha}4-helix from Glu175 to Ser195 in the mutant Egl-64. Marqusee and Baldwin (1987) showed that peptides containing (i + 4) Glu–Lys forms more stable helices than peptides containing (i + 3) Glu–Lys at extremes of pH (2 and 12) and also at pH 7. Hence we propose that the improved thermostability of the double mutation is due to the double intrahelical ion pairs (salt bridges). Moreover, the ion pair, the N{zeta} Lys194 hydrogen and the carbonyl Asp224 O{delta}2 may also be involved in the thermostabilization.



View larger version (15K):
[in this window]
[in a new window]
 
Fig. 5. Deviation of distances (A) between Lys179 (N{zeta}) and Glu175 (O{varepsilon}1 and O{varepsilon}2) in NK mutant and (B) between Lys194 (N{zeta}) and Glu190 (O{varepsilon}1 and O{varepsilon}2) or Asp224 (O{delta}1 and O{delta}2) in DK mutant as a function of time. The simulations were conducted in 4000 fs at 330 K (4000 steps) after linear heating from 0 to 330 K within 1000 fs (1000 steps). Each simulation was repeated five times and typical results are shown. Identifications of the lines used for the course of the deviation of distance between each ion pair are indicated.

 
The Lys residues gained at positions 179 and 194 are exposed to the hydrophilic enzyme surfaces (Hakamada et al., 2001Go). Generally, solvent-exposed salt bridges appear to play little role in the stabilization of proteins as observed in the cases of T4 lysozyme (Dao-pin et al., 1991Go) and barnase (Serrano et al., 1990Go; Sali et al., 1991Go). In fact, individual salt bridges often do not contribute significantly to protein stability (Dao-pin et al., 1991Go; Hendsch and Tidor, 1994Go; Wimley et al., 1996Go). However, the cumulative or cooperative effects of a number of salt bridges or ion pairs may contribute substantially to stability enhancement in thermophilic proteins (Chan et al., 1995Go; Hennig et al., 1995Go; Yip et al., 1995Go; Vetriani et al., 1998Go; Hashimoto et al., 1999Go) and synthetic {alpha}-helical peptides (Kohn et al., 1997Go). The salt bridge Lys179–Glu175 in the less stabilized NK mutant is located on the N-terminus of the {alpha}4-helix immediately after a long random coil from Ala158 to Asn174 (see Figure 2Go). Therefore, the contribution of this ionic interaction may be limited and small owing to the mobility around this region. In contrast, the salt bridge Lsy194–Glu190 in the highly stabilized DK mutant is located on the C-terminal end of the helix connected immediately to the relatively rigid ß-sheet architecture Ile203–Gly206 (ß5) with the short interruption Gly196–Ile202. The locations at both N- and C-terminal positions of the gained two intrahelical ion pairs may contribute, as if capping both terminals, to reinforce the stability of the {alpha}4-helix conformation, as illustrated in Figure 6Go.



View larger version (54K):
[in this window]
[in a new window]
 
Fig. 6. Model structure of Egl-64 showing the catalytic residues and the salt bridges generated by Lys substitutions. (A) Overall model structure of the NK/DK mutant enzyme; (B) close-up view showing the region around the intrahelical double salt bridges. The solid ribbon diagrams were drawn using the program WebLab Viewer and {alpha}-helices and ß-strands around the {alpha}4-helix region are colored red and cyan, respectively. The catalytic Glu residues are shown in the CPK representation.

 


    Notes
 
3 To whom correspondence should be addressed. E-mail: itos{at}jamstec.go.jp Back


    References
 Top
 Abstract
 Introduction
 Materials and methods
 Results and discussion
 References
 
Argos,P., Rossmann,M.G., Grau,U.M., Suborn,H., Frank,G. and Tratschin,J.D. (1979) Biochemistry, 18, 5698–5703.[ISI][Medline]

Chan,M.K., Mukund,S., Kletzin,A., Adams,M.W. and Rees,D.C. (1995) Science, 267, 1463–1469.[ISI][Medline]

Dao-pin,S., Anderson,D.E., Baase,W.A., Dahlquist,F.W. and Matthews,B.W. (1991) Biochemistry, 30, 11521–11529.[ISI][Medline]

Davies,G.J., Gambli,S.J., Littlechild,J.A. and Watson,H.C. (1993) Proteins: Struct. Funct. Genet., 15, 283–289.[ISI][Medline]

Davies,G.J., Dauter,M., Brzozowski,A.M., Bjørnvad,M.S., Andersen,K.V. and Schülein,M. (1998) Biochemistry, 37, 1926–1932.[ISI][Medline]

Declerck,N., Machius,M., Wiegand,G., Huber,R. and Gaillardin,C. (2000) J. Mol. Biol., 301, 1041–1057.[ISI][Medline]

Ducros,V., Czjzek,M., Beleïch,A., Gaudin,C., Fierobe,H.-P., Beleïch,J.-P., Davies,G.J. and Haser,R. (1995) Structure, 3, 939–949.[ISI][Medline]

Hakamada,Y., Koike,K., Yoshimatsu,T., Mori,H., Kobayashi,T. and Ito,S. (1997) Extremophiles, 1, 151–156.[ISI][Medline]

Hakamada,Y., Hatada,Y., Koike,K., Yoshimatsu,T., Kawai,S., Kobayashi,T. and Ito,S. (2000) Biosci. Biotechnol. Biochem., 64, 2281–2289.[ISI][Medline]

Hakamada,Y., Hatada,Y., Ozawa,T., Ozaki,K., Kobayashi,T. and Ito,S. (2001) FEMS Microbiol. Lett., 195, 67–72.[ISI][Medline]

Haney,P., Konisky,J., Koretke,K.K., Luthey-Schluten,Z. and Wolynes,P.G. (1997) Proteins, 28, 117–130.[ISI][Medline]

Hashimoto,H., Inoue,T., Nishioka,M., Fujisawa,S., Takagi,M., Imanaka,T. and Kai,Y. (1999) J. Mol. Biol., 292, 707–716.[ISI][Medline]

Hendsch,Z.S. and Tidor,B. (1994) Protein Sci., 3, 211–226.[Abstract/Free Full Text]

Hennig,M., Darimont,B., Sterner,R., Kirschner,K. and Jansonius,J.N. (1995) Structure, 3, 1295–1306.[ISI][Medline]

Horikoshi,K. (1999) Microbiol. Mol. Biol. Rev., 63, 735–750.[Abstract/Free Full Text]

Hoshino,E. and Ito,S. (1997) In Van Ee,J.H., Misset,O. and Baas,E. (eds), Enzymes in Detergency. Marcel Dekker, New York, pp. 149–174.

Imanaka,T., Shibazaki,M. and Takagi,M. (1986) Nature, 324, 695–697.[ISI][Medline]

Ito,S. (1997) Extremophiles, 1, 61–66.[ISI][Medline]

Ito,S., Kobayashi,T., Ara,K., Ozaki,K., Kawai,S. and Hatada,Y. (1998) Extremophiles, 2, 185–190.[ISI][Medline]

Kohn,W.D., Kay,C.M. and Hodges,R.S. (1997) J. Mol. Biol., 267, 1039–1052.[ISI][Medline]

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

Marqusee,S. and Baldwin,R.L. (1987) Proc. Natl Acad. Sci. USA, 84, 8898–8902.[Abstract]

Ozaki,K., Shikata,S., Kawai,S., Ito,S. and Okamoto,K. (1990) J. Gen. Microbiol., 136, 1327–1334.[ISI][Medline]

Ozaki,K., Hayashi,Y., Sumitomo,N., Kawai,S. and Ito,S. (1995) Biosci. Biotechnol. Biochem., 59, 1613–1618.[ISI][Medline]

Russell,R.J.M., Ferguson,J.M.C., Hough,D.W. and Taylor,G.L. (1997) Biochemistry, 36, 9983–9994.[ISI][Medline]

Sakon,J., Adney,W.S., Himmel, M.E., Thomas,S.R. and Karplus,P.A. (1996) Biochemistry, 35, 10648–10660.[ISI][Medline]

Sali,D., Bycroft,M. and Fersht,A.R. (1991) J. Mol. Biol., 220, 779–788.[ISI][Medline]

Serrano,L. and Fersht,A.R. (1989) Nature, 342, 296–299.[ISI][Medline]

Serrano,L., Horovitz,A., Avron,B., Bycroft,M. and Fersht,A.R. (1990) Biochemistry, 29, 9343–9352.[ISI][Medline]

Shikata,S., Saeki,K., Okoshi,H., Yoshimatsu,T., Ozaki,K., Kawai,S. and Ito,S. (1990) Agric. Biol. Chem., 54, 91–96.[ISI]

Shirai,T., Suzuki,A., Yamane,T., Ashida,T., Kobayashi,T., Hitomi,J. and Ito,S. (1997a) Protein Eng., 10, 627–634.[Abstract]

Shirai,T., Yamane,T., Hidaka,T., Kuyama,K., Suzuki,A., Ashida,T., Ozaki,K. and Ito,S. (1997b) J. Biochem., 122, 683–685.[Abstract]

Sumitomo,N., Ozaki,K., Kawai,S. and Ito,S. (1992) Biosci. Biotechnol. Biochem., 56, 872–877.[ISI][Medline]

Sumitomo,N., Ozaki,K., Hitomi,J., Kawaminami S., Kobayashi,T., Kawai,S. and Ito,S. (1995) Biosci. Biotechnol. Biochem., 59, 2172–2176.[ISI][Medline]

Takano,K., Ogasahara,K., Kaneda,H., Yamagata,Y., Fujii,S., Kanaya,E., Kikuchi,M., Oobatake,M. and Yutani,K. (1995) J. Mol. Biol., 254, 62–76.[ISI][Medline]

Vetriani,C., Maeder,D.L., Tollida,N., Yip,K.S., Stillman,T.J., Britton,K.L., Rice,D.W., Klump,H.H. and Robb,F.T. (1998) Proc. Natl Acad. Sci. USA, 95, 12300–12305.[Abstract/Free Full Text]

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

Wallon,G., Kryger,G., Lovett,S.T., Oshima,T., Ringe,D. and Petsko,G.A. (1997) J. Mol. Biol., 266, 1016–1031.[ISI][Medline]

Watanabe,K., Hata,Y., Kizaki,H., Katsube,Y. and Suzuki,Y. (1997) J. Mol. Biol., 269, 142–153.[ISI][Medline]

Wimley,W.C., Gawrisch,K., Creamer,T.P. and White.S.H. (1996) Proc. Natl Acad. Sci. USA, 93, 2985–2990.[Abstract/Free Full Text]

Yip,K.S. et al. (1995) Structure, 3, 1147–1158.[ISI][Medline]

Yoshimatsu,T., Ozaki,K., Shikata,S., Ohta,Y., Koike,K., Kawai,S. and Ito,S.(1990) J. Gen. Microbiol., 136, 1973–1979.[ISI]

Závodszky,P., Kardos,J., Svingor,Á. and Petsko,G.A. (1998) Proc. Natl Acad. Sci. USA, 95, 7406–7411.[Abstract/Free Full Text]

Received January 29, 2001; revised March 23, 2001; accepted April 10, 2001.