Previous Article | Next Article ![]()
Applied and Environmental Microbiology, February 2006, p. 1588-1594, Vol. 72, No. 2
0099-2240/06/$08.00+0 doi:10.1128/AEM.72.2.1588-1594.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Seung-Pyo Hong,2,
Jae Jun Song,1
Su-Jin Kim,1
Mi-Sun Kwak,2 and
Moon-Hee Sung2,3*
Laboratory of Microbial Function, KRIBB, Daejeon,1 Bioleaders Corp., Daejeon,2 Kookmin University, Seoul, Korea3
Received 5 July 2005/ Accepted 11 November 2005
|
|
|---|
-ketoglutarate decreased by half. The mutant was very close to the wild-type in thermal stability, indicating that the mutations did not disturb the overall structure of the enzyme. Homology modeling also suggested that the two tyrosine residues in the EYcY sequence from the Geobacillus D-AATs had a
/
interaction that was replaceable with the salt bridge interaction between the arginine and aspartate residues in the LRcD sequence. |
|
|---|
-keto acid substrates into their respective D-amino acids, some of which are indispensable for bacteria as peptidoglycan components of the cell wall (25). As such, the enzyme has been applied as a catalyst to produce optically pure D-amino acids (1, 3, 30) that act as intermediates of semisynthetic antibiotics, bioactive peptides, and other physiologically active compounds (20, 21). D-AAT activity is found in various gram-positive bacteria, including Bacillus (8, 28, 29), Staphylococcus (22), and Listeria (31), and yet recent biotechnology studies have mainly focused on thermophilic or mesophilic Bacillus enzymes due to their high catalytic activity and broad substrate specificity (4, 9, 10, 23). For example, the D-AAT from thermophilic Bacillus sp. strain YM1 was remarkable in its activity and thermal stability and showed a high identity of 67% with a mesophilic Bacillus sphaericus enzyme, while representing a limited identity of less than 50% with other Bacillus enzymes.
16S rRNA analyses are useful for comparing phenotypically close and yet genetically different microorganisms (6). For instance, in the phylogenetic analysis of bacilli, the thermophilic YM1 has been assigned to genetic group II, together with B. sphaericus, whereas most other thermophilic bacillus species, such as B. stearothermophilus, B. thermodenitrificans, and B. kaustophilus, have been positioned in genetic group V (2), recently renamed as the genus Geobacillus (17, 19).
Consequently, Geobacillus species would seem to be a remarkable resource for new D-AATs with unique sequences and enzymatic properties. Accordingly, the present study presents new thermostable D-AATs from the genus Geobacillus, including a characterization of their genetic and catalytic properties, and mutational studies of the conservative ß5-
5 loop (bearing the cofactor-binding lysine) (26). In addition, homology modeling is conducted to analyze the functional relevance of the distinctive loop structure of the Geobacillus D-AATs.
|
|
|---|
DNA manipulation and mutagenesis.
The genomic DNA of the Geobacillus strains was partially digested with Sau3AI. DNA fragments of 3 to 10 kb were then isolated by centrifugation on a sucrose gradient (5 to 40% [wt/vol]) for 20 h at 25,000 rpm in a Beckman SW40 rotor and ligated into the BamHI site of pUC118 at 16°C for 12 h with a T4 DNA ligase. Thereafter, the E. coli WM335 was transformed with the ligation mixture with electroporation.
The site specific mutagenesis to introduce the desired mutations into the target DNA sequences was performed by using the megaprimer PCR method. The mutagenic internal primer 5-ACAAGAGATGTCCGCTGGCTACGTTGCGATATTAAGAGTTTAAATCTTCTA-3 is designed to bear an LRcD residue (underlined) instead of the wild-type sequence, EYcY, whereas the C-terminal primer 5-GCCGGATCCTTATTTTGCGTTTTTGACAGC-3 is designed to have a BamHI site (underlined). The product of the first PCR with the mutagenic primer and C-terminal primer was purified and then used as a megaprimer for a second PCR, along with the N-terminal primer 5-GCATTAAAGCTGTACGTACTA-3. The final PCR product contained the desired mutation and 3-terminal BamHI site in the DNA sequence. Plasmid pHCE19T(II) was used for the expression of the LRcD mutant: the plasmid vector was digested with NcoI, blunt ended by Klenow treatment, and sequentially digested with BamHI. The resulting vector was ligated with the second PCR product by a blunt-cohesive ligation at 16°C with a T4 DNA ligase.
Expression and purification.
E. coli XL1-Blue cells bearing the plasmid pDSK2 or pHKLS23 were cultivated at 37°C for 16 h in 1 liter of LB medium containing 100 µg of ampicillin/ml. After being harvested by centrifugation, the cells were disrupted by sonification in a standard buffer, including 30 mM Tris-HCl (pH 8.0), 0.01% ß-mercaptoethanol, and 0.05 mM pyridoxal-5'-phosphate. The active D-AAT was recovered from the supernatant of the cell lysate and incubated at 60°C for 20 min to remove heat-labile E. coli proteins. The resulting enzyme solution was then loaded onto a Resource Q ion exchange (Pharmacia, Sweden), washed with the standard buffer, and eluted by using a potassium chloride gradient from 0 to 0.5 M. Next, the active fractions were collected, adjusted to include 1.7 M ammonium sulfate, and loaded onto a Phenyl Superose (Pharmacia, Sweden). The elution was carried out by using a reverse gradient of ammonium sulfate from 1.7 to 0 M, and then the active fractions were pooled and dialyzed against the standard buffer and stored in a deep freezer. All of the column procedures were carried out by using an AKTA System (Amersham Bioscience, Sweden) at room temperature.
Sequence analysis.
The nucleotide sequences of the D-AAT and 16S rRNA genes were determined at Solgent Co. (Daejon, South Korea) by primer walking on both DNA strands and then analyzed by using the Vector NTI suite (Informax, Maryland).
The PCR amplification of the 16S rRNA gene was performed with two universal primers, as previously described (33). The 16S rRNA gene sequences were then aligned with reference sequences from related taxa by using CLUSTAL W software (32). Gaps at the 5' and 3' ends of the alignment were omitted from further analyses. A phylogenetic tree was constructed by using the neighbor-joining method (as implemented with the NEIGHBOR program of the PHYLIP package) (24).
Homology modeling and structural analysis.
The D-AAT structures of the PDB entries 1DAA, 2DAA, and 3DAA were aligned with the 3-D Align module from DS modeling 1.1 (Accelrys, San Diego, CA). The resulting consensus structure was then used as the template for the 2-D Align with the D-AAT sequences. Three-dimensional homology models of the SK1 D-AAT were generated by using the Modeler module (7), and the model with the best loop conformations determined by using the Profiles-3-D verification method of the Protein Health module for further modeling (14). A coenzyme analog, D-cycloserine pyridoxalphosphate, was extracted from the 2DAA and fitted to the PLP-binding site of each monomer. Energy minimization and molecular dynamics (MD) to optimize the structure were carried out by using the DS CHARMm module in vacuo, using the steepest descent method first, followed by the conjugate gradient method (15). During the minimization process, the protein backbone was restrained by a harmonic constraint. The lowest energy conformer resulting from the MD run was then finally relaxed without restraints.
Enzyme and protein assay.
The D-AAT activity was calculated from the pyruvate formation rate using a coupling assay with lactate dehydrogenase (Roche Diagnostics, Switzerland) or the salicylaldehyde method (5). The assay mixture contained 0.05 mM pyridoxal-5'-phosphate, 0.5 mM NADH, 10 µg of lactate dehydrogenase/ml, and various concentrations of D-alanine and
-ketoglutarate in a 100 mM Tris-HCl buffer (pH 8.5). One unit of enzyme was defined as the activity catalyzing the formation of 1 µmol of pyruvate per min at 50°C. The protein concentration was determined by the Bradford method with bovine serum albumin as the standard.
|
|
|---|
-ketoglutarate as the substrates. The screening of D-AAT activity at 55°C identified two strains with the highest activity, subsequently designated LK1 and LK2. However, these highly active strains were unable to grow at 65°C, the temperature at which most Geobacillus species grow, and in microscopic and physiological examinations were found to be very close to the well studied Bacillus sp. strain YM1 (12). Thus, three D-AAT-producers (SK1, KLS1, and SD1) growing at 65°C were selected for further experiments. The average activity of the selected strains, 0.02 U/mg of protein, was approximately one-fifth the level of the highly active Bacillus enzymes growing at 55°C. A Western blotting analysis with a polyclonal YM1 D-AAT antibody exhibited strong immuno-signals toward the LK1 and LK2 proteins (Fig. 1, lanes 1 to 3), coinciding with their microbial similarity to the YM1 strain. However, no immuno-signals were detected against SK1, KLS1, and SD1 (lanes 4 to 7 in Fig. 1). Since the depth of the immunoblotting reflects the relevancy between the primary sequences, the three D-AATs isolated at 65°C were determined as sequentially distant from the YM1 D-AAT.
![]() View larger version (33K): [in a new window] |
FIG. 1. Comparison of D-amino acid aminotransferases from thermophilic Bacilli strains using Western blotting analyses. Lanes 1 to 7 are the immunostains for the cell extracts from different Bacillus strains: 1, YM1; 2, LK1; 3, LK2; 4, SK1; 5, KLS1; 6, KL1; and 7, SD1. The polyclonal antibody was prepared from an antiserum boosted by the YM1 D-AAT. Lanes 8 to 11 are the purified D-AATs: 8 and 10, Bacillus sp. strain YM1; 9 and 11, G. toebii SK1. Lanes 10 and 11 are the sodium dodecyl sulfate-polyacrylamide gel electrophoresis analysis of the purified YM1 and SK1 proteins.
|
Cloning and expression of D-AATs.
For a genetic complementation assay, Escherichia coli WM335 cells containing the genomic DNA of the Geobacillus strains were spread on LB-ampicillin plates without the supplementation of D-glutamate (13). After being incubated overnight, one or two colonies were picked from each plate and transferred to fresh culture media to isolate the plasmids. The resulting pUC118 plasmids were retransformed into E. coli XL1-Blue cells to examine the D-AAT activity. A clone pDSK2 (insert size 3.6 kb) from the SK1 library showed the highest activity of 35 U per mg of protein (1,750-fold increase from the original activity of the SK1 strain), the activity attributed to the 1.7-kb EcoRI-XhoI fragment based on a restriction analysis. Thus, the 1.7-kb fragment was subcloned into pBluescript II SK to construct pDSK231. Meanwhile, a clone pKLS1 was selected from the KLS1 library, and restriction analyses attributed the enzyme activity to the 1.86-kb HindIII-SacI fragment. Thus, the 1.86-kb fragment was subcloned into pBluescript II SK to construct pHKLS23.
Nucleotide sequencing of the 1.7-kb insert in pDSK231 showed an 864-bp open reading frame (gb: DQ100074) that encoded a protein of 288 amino acid residues with a molecular weight of 33,198. These results were in accordance with the results from sodium dodecyl sulfate-polyacrylamide gel electrophoresis. The 1.86-kb insert in pHKLS23 was also sequenced and the open-reading frame (gb: DQ100075) homologous to the SK1 D-AAT was designated KLS1 D-AAT.
Sequence analysis and structural alignment.
The alignment of the D-AAT sequences showed that almost all of the residues with key roles in the YM1 enzyme (26) were conserved in the Geobacillus enzymes, including the catalytic K145 (K149 in the SK1 D-AAT), carboxylate traps that define the position and orientation of the substrate-binding Y31 (Y35), R98 (R103), and H100 (H105), and cofactor-binding residues R50 (R54), I204 (I208), T241 (T245), E177 (E181), and R138 (R142), all shown in Fig. 2. The loop from Ser240 to Ser243, thought to be important in substrate specificity, was also conserved.
![]() View larger version (127K): [in a new window] |
FIG. 2. Sequence alignment of Bacillus D-AATs. Abbreviations: G.SK1, G. toebii SK1; G.KLS1, Geobacillus sp. strain KLS1; G.kau, G. kaustophilus; B.YM1, Bacillus sp. strain YM1, B.sph, B. sphaericus; B.sub, B. subtilis; B.lic, B. licheniformis; B.cer, B. cereus; B.ant, B. anthracis; B.hal, B. halodurans; L.mon, Listeria monocytogenes; S.hae, Staphylococcus haemolyticus. The conserved residues are shaded dark, while similar residues are shaded light. The asterisk indicates the lysine residue that forms a sciff base with pyridoxal 5'-phosphate.
|
5 covering the catalytic Lys residue (26). Although the sequence of the loop, 138-RWLRCDIKSLNLL-150, has been well conserved in the D-AAT sequences identified to date, it was distinctively changed in the Geobacillus D-AATs, where L140, R141, and D143 in the YM1 D-AAT were substituted with residues 144E, 145Y, and 147Y, respectively (Fig. 2). The overall sequence of the SK1 D-AAT exhibited a 90.2% identity with the KLS1 enzyme and 71.5% identity with a hypothetical correspondent from the G. kaustophilus genome. However, the identity level with other D-AATs, including Bacillus sp. strain YM1, was lower than 45%.
When the alignments of the Bacilli D-AATs were imported into the PHYLIP program, the Geobacillus enzymes from the present study were located in the same branch of the tree phenogram (Fig. S2 in the supplemental material). The enzymes from YM1 and B. sphaericus in group II were located in the same group, revealing a 67% identity with each other, whereas the enzymes from B. cereus and B. anthracis, with an 89.7%-identity, formed a separate branch to the enzymes from B. subtilis and B. licheniformis that belong to the same genetic group I.
Purification and characterization of kinetic properties.
The E. coli XL1-Blue cells bearing the plasmid pDSK2 exhibited a thick protein band with a molecular mass of 34 kDa. The expressed protein was estimated to be more than 30% of the total E. coli proteins when quantified with an image analyzer (Bio-Rad, California). After heat treatment at 60°C for 20 min and two sequential chromatographies (ion exchange and hydrophobic), the SK1 D-AAT was purified to homogeneity (lane 11 in Fig. 1) at a recovery yield of 59%. The purification protocols were also successfully applied to the site-directed mutant of the SK1 D-AAT, the KLS1 enzyme, and YM1 enzyme expressed in E. coli XL1-Blue (23).
The purified D-AATs from the SK1 and KLS1 strains showed absorption maxima at 330 and 415 nm, corresponding to pyridoxal-5-phosphate bound to the active site lysine residue (23). The molecular masses of the enzymes were also estimated to be around 70 kDa by gel permeation chromatography on a Superose 12 column (Pharmacia), indicating that the enzymes consisted of two identical subunits.
The kinetic characterization of the purified enzymes was performed based on triplicate experiments at 50°C, while changing the concentrations of D-alanine and
-ketoglutarate from 1 to 40 mM. Lineweaver-Burk plots of the bi-reactant reaction showed characteristic parallel lines for a ping-pong-bi-bi reaction. The catalytic rate of the D-AATs from SK1 and KLS1 was determined to be 47 and 55 s1, respectively, which was close to one-third of the highly active YM1 D-AAT (Table 1).
|
View this table: [in a new window] |
TABLE 1. Kinetic parameters of D-AATs
|
-KG values were remarkably smaller for the Geobacillus enzymes. In particular, the Km,
-KG for the KLS1 enzyme was 12 times smaller than that for the YM1 enzyme. Meanwhile, for the Geobacillus enzymes, the Km,D-Ala values were much larger than the Km,
-KG values, in distinct contrast to the same constants for the YM1 enzyme.
Investigations of the substrate preference of the SK1 D-AAT revealed that D-alanine, D-aminobutyrate, and D-aspartate were the best amino donors for
-ketoglutarate (Table 2), whereas D-serine, D-asparagine, and D-norvaline also served as amino donors, but to a lesser extent. Moreover, investigation of the amino acceptors showed that
-ketoglutarate,
-ketobutyrate, and
-ketoisovalrate were the preferred amino acceptors.
|
View this table: [in a new window] |
TABLE 2. Substrate specificity of D-AAT from Geobacillus toebii SK1
|
![]() View larger version (24K): [in a new window] |
FIG. 3. Effect of temperature and pH on SK1 D-AAT. The open symbols indicate the native enzyme, while the closed symbols with dotted lines indicate the mutant enzyme with an LRcD sequence instead of 144-EYcY-147. (A) The native and mutant enzymes were incubated at different temperatures for 20 min in a 0.1 M Tris-HCl buffer (pH 8.5), and the remaining activities determined at 50°C to evaluate the respective thermal stabilities. (B) The native and mutant enzymes were assayed at different reaction temperatures for 20 min. The concentrations of pyruvate from D-alanine were measured by using the salicylaldehyde method. (C) The native enzyme was incubated in different buffers for 1 h at 50°C, and the remaining activities were assayed to evaluate the pH stability. (D) The native enzyme was assayed for 20 min at different pHs with the following 0.1 M buffers: sodium acetate buffer (pH 4.0 to 6.0), potassium phosphate buffer (pH 6.0 to 8.0), Tris-HCl buffer (pH 7.5 to 9.0), and N-cyclohexyl-3-aminopropanesulfonic acid buffer (pH 9.0 to 12.0).
|
5 loop sequence.
-ketoglutarate increased twofold (Table 1). Although the thermal stability of the mutant did not change very much, its maximum activity was identified at 60°C, and the change in activity with temperature was less steep than that for the native enzyme, as represented in Fig. 3A and B (solid symbols for the mutant). The increase in enzyme activity with temperature was interpreted by the Arrhenius equation [ln(k) = ln(A) (Ea/RT)] to compare the activation energies (Ea) for the native and mutant enzymes (Fig. S3 in the supplemental material). From the slope of ln(k) versus 1/T, the activation energies for the native and mutant enzymes were calculated as 29 and 16 kJ/mol, respectively.
Homology modeling of SK1 D-AAT.
Comparative modeling with DS Modeling 1.1 generated hypothetical structures for the wild and mutant type SK1-DAATs (Fig. 4A), which only differed in the nonaligned loop region. The models displaying the highest stereochemical quality were then determined by using the Protein Health module. The resulting model contained no cis-proline, and the final RMS values for the C
trace of the wild and mutant enzymes after 200 ps of MD were 1.2 and 1.0 Å, respectively. The estimated structures were very stable during the MD simulations, and at least 95% of the backbone
and
dihedral angles for the residues were located within the allowed regions of the Ramachadran plot (70% core).
![]() View larger version (46K): [in a new window] |
FIG. 4. Homology model for G. toebii SK1 D-AAT. (A) 3-D Align analysis of the monomeric units of the wild and mutant enzymes; (B) magnified structures of the 143-wEYcYik-149 loop. The loop structure bears the cofactor binding lysine (underlined) and is located at the crevice between the N and C domains. The 143W is related to the intersubunit packing to form a catalytic dimer. The capital letters in EYcY indicate the residues substituted for LRcD in the mutant enzyme. The structures of the native and mutant enzymes are represented in green and blue, respectively. The gray molecule shows the structure of D-cycloserine pyridoxalphosphate, adopted from the PDB-entry 2DAA.
|
/
interaction between the aromatic rings (18). Plus, the salt-bridge interaction between the arginine and aspartate was reconstituted at a 2.7-Å distance, using the molecular modeling and dynamics for the LRcD mutant of the SK1 enzyme. |
|
|---|
The sequence and structural alignments revealed that most of the residues with key roles in the YM1 enzyme were conserved in the SK1 D-AAT, as described in Results. In previous structural studies of the YM1 D-AAT, the 240-STTS-243 loop has been interpreted as the substrate-binding pocket, and was conserved as 244-STTS-247 in the SK1 D-AAT. An mutational study performed to broaden the substrate specificity by replacing the loop with a SVSS sequence, as found in the B. sphaericus enzyme with a broader substrate range, resulted in a 30% lower activity compared to the original and no significant expansion in the substrate specificity was observed (unpublished data).
In previous crystallographic studies of the YM1 D-AAT, the long loop between ß5 and
5 was located in a crevice between the N and C domains and included the cofactor binding lysine residue (26). Plus, mutational and biochemical analyses have indicated that the 139W in the loop plays an important role in subunit packing to form a catalytic dimer, and any mutations replacing the 139W with a nonaromatic residue were detrimental to the enzyme activity (11, 16).
In the primary sequence of the SK1 D-AAT, the EYcY sequence was located between the 143W (139W in YM1 D-AAT) and the cofactor binding lysine residue, 149K (145K). Therefore, substituting LRcD for EYcY was expected to lower the activity and stability of the enzyme by disturbing the active site or subunit interface. However, the mutant enzyme in the present study maintained its stability and even increased its activity up to 68%, while the Michaelis constant for
-ketoglutarate increased twofold (Table 1). Therefore, these results suggest that the LRcD sequence of highly active homologue is recruited effectively and fitted into the structure of the SK1 enzyme. When the profiles in Fig. 3B were replotted based on Arrhenius' equation (Fig. S3 in the supplemental material), the mutant enzyme showed a decreased activation energy (Ea = 16 kJ/mol) compared to that of the wild type (Ea = 29 kJ/mol). The lower activation energy and decreased affinity to substrates as exhibited by the mutant may imply an increased flexibility in the local region adjacent to the active site.
The arginine and aspartate in the LRcD sequence of the YM1 D-AAT were analyzed to form a stable salt bridge within a 3-Å distance (26), whereas homology modeling of the SK1 D-AAT suggested that the two tyrosine residues in the EYcY sequence of the SK1 enzyme were located in a proximate distance to form a
/
interaction (18). Accordingly, the catalytic results and molecular simulation imply that the salt bridge in the YM1 D-AAT and presumed aromatic interaction in the Geobacillus D-AATs are interchangeable with each other, despite the significant difference in molecular size and chemical properties.
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
S.-G.L. and S.-P.H. contributed equally to this study. ![]()
|
|
|---|
-keto acids with Escherichia coli cells expressing heterologous genes. Appl. Environ. Microbiol. 61:4651-4656.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»