Previous Article | Next Article ![]()
Journal of Bacteriology, July 2007, p. 5265-5275, Vol. 189, No. 14
0021-9193/07/$08.00+0 doi:10.1128/JB.00352-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Department of Microbiology and Institute for Genomic Biology,1 Department of Biochemistry, University of Illinois, Urbana, Illinois,2 Department of Microbiology and Immunology, The University of Oklahoma Health Sciences Center, Oklahoma City, Oklahoma3
Received 8 March 2007/ Accepted 1 May 2007
|
|
|---|
|
|
|---|
Inhibition of cell wall synthesis remains an effective approach for preventing bacterial growth (53, 54). For most eubacteria, an important constituent of the peptidoglycan cell wall is D-glutamate (38, 39, 41, 44, 45). D-Glutamate is not typically available in the environment but instead is generated by the enzyme glutamate racemase (E.C. 5.1.1.3), which catalyzes the reversible stereoinversion of L-glutamate (5, 14, 31). Insights into the cofactor-independent amino acid racemases have begun to emerge from biochemical studies of enzymes isolated from several organisms, including Bacillus subtilis, Bacillus pumilus, Bacillus sphaericus, Escherichia coli, Lactobacillus fermentum, Lactobacillus brevis, Aquifex pyrophilus, Staphylococcus haemolyticus, Brevibacterium lactofermentum, and Mycobacterium tuberculosis (1, 2, 5, 10, 14, 18, 19, 21, 28, 29, 33, 35, 40, 47, 58, 61), as well as the recently described crystal structure of RacE-D-glutamate from B. subtilis (43). Several studies have identified glutamate racemase as an essential gene in B. subtilis and E. coli, which has led to the prediction that glutamate racemase activity is important for peptidoglycan biosynthesis in these organisms (14, 31). Because glutamate racemases are not found in mammals, these enzymes have emerged as excellent targets for the design of a new class of antibacterial agents (3, 12, 22, 28, 43, 59).
Analogous to the peptidoglycan of other eubacteria, D-glutamate is predicted to be an important constituent of B. anthracis (48). However, in B. anthracis D-glutamate is also the sole component of the poly-
-D-glutamic acid (PDGA) capsule (24), an important virulence factor that is required for dissemination in a murine model of inhalational anthrax (16) and is presumed to be required for disease in humans as well (16, 36, 46). Therefore, in addition to its proposed role in cell wall biosynthesis, glutamate racemase is also predicted to be the major source of D-glutamate for PDGA capsule synthesis in B. anthracis (9). In contrast to most bacteria that possess only one glutamate racemase gene (5, 10, 14, 18, 19, 26, 29, 33, 35, 40, 61), the B. anthracis genome contains two genes (BAS0806 and BAS4379) predicted to encode glutamate racemases, which are designated racE1 and racE2. Recently, inactivation of racE2 was reported to cause a severe growth defect in B. anthracis, while inactivation of racE1 only moderately inhibited growth (52). However, the underlying reasons for these growth phenotypes, especially for the racE2 mutant, were not identified and thus cannot be attributed at this time to defects in peptidoglycan synthesis resulting from insufficient D-glutamate availability (52). Further complicating the interpretation of the phenotypes reported for the racE1 and racE2 mutants is uncertainty concerning whether either of these two genes encodes functional enzymes with glutamate racemase activity. Thus, validation of RacE1 or RacE2 as a potential drug target awaits elucidation of the enzymatic and biochemical properties of these proteins.
Here, we characterized and compared recombinant wild-type and mutant forms of RacE1 and RacE2 from B. anthracis. The racE1 and racE2 genes were cloned from B. anthracis Sterne 7702 and expressed as recombinant proteins in E. coli. These studies revealed that in cell-free assays, B. anthracis RacE1 and RacE2 both catalyze stereoinversion of L-glutamate to D-glutamate. However, the two enzymes differ in ways that could influence future inhibitor development.
|
|
|---|
Multiple-sequence alignment. To compare the primary amino acid sequences of B. subtilis RacE and B. anthracis RacE1 and RacE2, a multiple-sequence alignment was generated. Amino acid sequences for B. anthracis Sterne RacE1 (BAS0806), B. anthracis Sterne RacE2 (BAS4379), and B. subtilis RacE (BSU2835) were aligned using ClustalW (http://align.genome.jp/). The multiple-sequence alignment file was then entered into the ESPript V2.2 alignment program (http://mail.bic.nus.edu.sg/ESPript/cgi-bin/ESPript.cgi). Secondary structural elements and solvent accessibility indices were imported from the crystal structure of B. subtilis 168 RacE (Protein Data Bank accession no. 1ZUW) (43).
Cloning of racE1 and racE2 and preparation of expression strains. To obtain genomic DNA, B. anthracis Sterne 7702 was cultivated at 37°C with aeration on a rotary shaker in brain heart infusion medium (3.7% Bacto brain heart infusion, Millipore deionized water, 0.5% glycerol shaker). DNA was isolated from mid-log-phase cultures and was purified using the DNeasy tissue kit. B. anthracis Sterne racE1 (BAS0806) and B. anthracis Sterne racE2 (BAS4379) were PCR amplified using primers corresponding to the 5' and 3' ends of each gene (Table 1). These primers were engineered such that 5' XhoI (racE1) or SalI (racE2) and 3' BamHI (racE1 and racE2) restriction sites were incorporated. Each PCR product was purified using the PCR purification kit. The purified amplicons were incubated with XhoI (racE1) or SalI (racE2) and BamHI to generate directional annealing sites. The amplicons were then ligated with pET-15b to replace the XhoI-BamHI fragment within the polylinker region. The ligation mixtures were introduced into E. coli XL1-Blue by electroporation. The integrity of each gene from individual clones was confirmed by DNA sequencing. pET-15b-racE1 and pET-15b-racE2 were isolated using the plasmid miniprep kit and introduced by electroporation into E. coli T7 lysogen BL21(DE3).
|
View this table: [in a new window] |
TABLE 1. Cloning of B. anthracis glutamate racemase genes, using B. anthracis Sterne 7702(pXO1+/pXO2)
|
Protein concentrations were quantified by absorbance spectroscopy using the method described by Gill and von Hippel (20). Briefly, the extinction coefficient at 280 nm was calculated based upon the amino acid composition of the polyhistidine-tagged recombinant protein. The resulting extinction coefficients (17,420 M1 cm1 for RacE1 and 21,430 M1 cm1 for RacE2) were then used to determine the concentration of the appropriate protein solution utilizing a Multiskan Spectrum UV spectrophotometer from Thermo Electron Company (Waltham, MA). The protein stocks were stored at 20°C in storage buffer (50 mM Tris, 100 mM NaCl, 0.2 mM DTT; pH 8.0) with 20% glycerol at a final concentration of 10 mg/ml for up to 2 months.
Preliminary experiments revealed no detectable differences in the kinetic properties of RacE1 or RacE2 in the presence or in the absence of the amino-terminal polyhistidine fusion peptides (data not shown). Based on these data, purified proteins with amino-terminal polyhistidine fusion peptides were utilized for all of the experiments described in this study.
CD spectroscopy of purified RacE1 and RacE2. Circular dichroism (CD) spectra were collected for RacE1 and RacE2 in the far-UV range utilizing a J-720 CD spectropolarimeter from JASCO (Easton, MD). A cylindrical cuvette with a total volume of 350 µl and a path length of 0.1 cm was used for each assay. The CD spectra of RacE1 (2.7 µM) and RacE2 (2.4 µM) in optically clear borate buffer (50 mM potassium borate, pH 8.0) were recorded from 190 to 260 nm at a scan rate of 50 nm/s with a 1-nm wavelength step and with five accumulations.
Data acquisition was coordinated using the JASCO Spectra Manager v1.54A software. Raw data files were uploaded onto the DICHROWEB online server (http://www.cryst.bbk.ac.uk/cdweb/html/home.html) and analyzed using the CDSSTR algorithm with reference set 4, which is optimized for the analysis of data recorded in the range from 190 to 240 nm (34).
Size exclusion chromatography. Size exclusion chromatography was conducted using an AKTA Purifier 900 fast protein liquid chromatography (FPLC) system equipped with a Superdex 200 10/300 GL size exclusion column and a UV detector, all obtained from Amersham Pharmacia Biotech (Little Chalfont, United Kingdom). RacE1 (5 mg/ml; 100 µl), RacE2 (5 mg/ml; 100 µl), or a gel filtration standard mixture was injected onto the column preequilbrated with a potassium borate buffer (50 mM boric acid, 100 mM KCl, 0.2 mM DTT; pH 8.0) liquid phase at a flow rate of 0.5 ml/min. Standard curves were generated by plotting the log of the molecular weights (provided by the supplier) of the gel filtration standards versus retention times. Experimental retention times were used to calculate the apparent molecular weights of RacE1 and RacE2 from the standard curve.
Racemization assays. Enzyme-catalyzed stereoisomerization of D- or L-glutamate was assayed using a J-720 CD spectropolarimeter from JASCO (Easton, MD). A thermostat-equipped cylindrical cuvette with a capacity of 700 µl and a path length of 1 cm was used for each assay. D- or L-glutamate (5 mM) in optically clear potassium borate buffer (50 mM boric acid, 100 mM KCl, 0.2 mM DTT; pH 8.0) was incubated at 25°C in the absence or presence of RacE1 or RacE2 at a concentration of 0.08, 0.33, or 1.3 µM. D- or L-glutamate stereoisomerization was monitored by recording the CD signal at 217 nm. Data acquisition was performed using the JASCO Spectra Manager v1.54A software, and a nonlinear curve fit was applied using GraphPad Prism V4.03 from GraphPad Software (San Diego, CA).
Determination of steady-state kinetic parameters. Racemase assays were carried out as described above, with the exception that the D-glutamate concentration was varied from 0.2 to 10 mM while the concentration of L-glutamate was varied from 5 to 200 mM. In addition, a 700-µl cuvette with a 1-cm path length was used for reactions with substrate concentrations less than 5 mM, and a 350-µl cuvette with a 0.1-cm path length was used for reactions with higher substrate concentrations. Reactions were initiated by addition of RacE1 (0.78 µM) or RacE2 (0.78 µM), and the levels of D- or L-glutamate were monitored by recording the CD signal at 215 nm for D- or L-glutamate concentrations from 0.2 to 10 mM and at 225 nm for L-glutamate concentrations from 30 to 200 mM at 25°C. Data acquisition was performed using the JASCO Spectra Manager v1.54A software, and a nonlinear curve fit was applied using GraphPad Prism V4.03 from GraphPad Software (San Diego, CA).
pH rate profile.
The stereoisomerization of glutamate in the L
D direction by RacE1 and RacE2 in buffers having various pH values was assayed using a J-720 CD spectropolarimeter from JASCO (Easton, MD). A cylindrical cuvette with a total volume of 350 µl and a path length of 0.1 cm was used for each assay. Seven different buffers spanning a pH range from 6.5 to 9.5 with increments of 0.5 pH unit were prepared. To maintain a well-buffered system, the following phosphate and borate buffer formulations were utilized: for pH 6.5, 7.0, and 7.5, 50 mM potassium phosphate, 100 mM KCl, 200 mM L-glutamate, 0.2 mM DTT; and for pH 8.0, 8.5, 9.0, and 9.5, 50 mM boric acid, 100 mM KCl, 200 mM L-glutamate, 0.2 mM DTT. Each buffer was prepared with 200 mM L-glutamate, which yielded fully saturating conditions for RacE1 and nearly saturating conditions (83% saturation) for RacE2, so that the initial rate data would report true kcat values. Reactions were initiated by addition of RacE1 (0.78 µM) or RacE2 (0.78 µM), and the levels of L-glutamate were monitored by recording the CD signal at 225 nm at 25°C. Data acquisition was performed using the JASCO Spectra Manager v1.54A software, and a user-defined curve fit was applied using GraphPad Prism V4.03 from GraphPad Software (San Diego, CA).
Homology models. The homology models for RacE1 and RacE2 were constructed using The Chemical Computing Group's Molecular Operating Environment (MOE) 2006.08. The template for both models was the B. subtilis RacE-D-glutamate structure (Protein Data Bank accession no. 1ZUW), which was aligned with the sequences for RacE1 and RacE2 using the Blosum62 substitution matrix. Ten intermediate homology models resulting from permutational selection of different loop candidates and side chain rotamers were built for RacE1 and RacE2. The intermediate model which scored best according to a packing evaluation function was chosen as the final model. Each of the intermediate models was subjected to a degree of energy minimization using the force field MMFF94x, with a distance-dependent dielectric (i.e., it simulated the polar environment of water).
Docking of compound 69 to RacE1 and RacE2. A conformational database was generated for (2R,4R)-2-amino-4-(2-benzo[b]thienyl)methyl pentanedioic acid (compound 69) (12) by using a stochastic conformational search, as implemented in MOE 2006.08 (Chemical Computing Group, Inc.). This program employs a variation of the method of Ferguson and Raber (17), in which bonds are randomly rotated, rather than using perturbation of Cartesian coordinates. The force field was MMFF94x. Minimization was performed for each conformation up to a root mean square gradient of 0.001. Any two conformers were considered identical if their optimal heavy atom root mean square superposition distance was less than a tolerance value of 0.1 Å. All conformations with an energy greater than 7 kcal/mol were excluded from the database. This yielded a conformational database of 17 unique conformations of compound 69, which were used in the docking procedure. The docking of compound 69 into RacE1 and RacE2 was performed with the Dock function in MOE 2006.08, using the Alpha Triangle placement method, and the London dG Scoring method for free energy estimation.
Site-directed mutagenesis. Mutagenesis was performed using the QuikChange mutagenesis kit from Stratagene (La Jolla, CA). First, complementary mutagenic primers (Table 2) were engineered with the desired mutation in the center of the primer and 10 to 15 bases of correct sequence on either side. Reaction mixtures were prepared as described in the QuikChange protocol with pET-15b-racE1 or pET-15b-racE2 as the plasmid template for generation of the RacE1 and RacE2 mutants. After cycling of the reaction mixture 18 times in a thermal cycler, the resulting mixture was digested with DpnI, and the resulting DNA was transformed into supercompetent E. coli XL1-Blue cells. The resulting mutant plasmids were isolated, and the entire gene was sequenced to ensure that the appropriate mutations were introduced, while the rest of the gene sequences remained unchanged.
|
View this table: [in a new window] |
TABLE 2. Primer sequences used for mutagenesis
|
|
|
|---|
![]() View larger version (55K): [in a new window] |
FIG. 1. B. anthracis RacE1 and RacE2 possess significant sequence homology to B. subtilis RacE. Primary amino acid sequences for B. anthracis RacE1 and RacE2 and B. subtilis RacE were aligned utilizing ESPript V.2.2 (http://espript.ibcp.fr/ESPript/ESPript/). Secondary structural elements and solvent accessibility indices were imported from the crystal structure of B. subtilis 168 RacE (Protein Data Bank accession no. 1ZUW) (43). For solvent accessibility indices, the darker shaded regions indicate sections in the folded proteins that are exposed to solvent.
|
![]() View larger version (29K): [in a new window] |
FIG. 2. RacE1 and RacE2 are both functional enzymes in a cell-free system. (A) Purification of recombinant RacE1 and RacE2. The RacE1 and RacE2 clarified lysates and fractions from nickel-chelate affinity chromatography were analyzed by 12% SDS-PAGE, followed by Coomassie brilliant blue G-250 staining. (B) RacE1 and RacE2 secondary structure. CD spectra in the far-UV region (190 to 260 nm) were recorded for RacE1 and RacE2 (2.7 and 2.4 µM, respectively, both in 50 mM potassium borate buffer; pH 8.0). (C) Gel filtration chromatography. The sizes of purified RacE1 and RacE2 were determined by size exclusion FPLC. The molecular weights (MW) of RacE1 and RacE2 were calculated from the retention times of the peak absorbance by comparison with calibration standards having known molecular weights. (D) Racemization of glutamate in the L D direction. RacE1 and RacE2 were assessed for the capacity to convert L-glutamate to the corresponding D enantiomer by using CD to directly observe the loss of L-glutamate as it was converted to D-glutamate. L-Glutamate (5 mM) was incubated in the absence or presence of RacE1 or RacE2 (1.3, 0.33, or 0.08 µM), and the CD signal was recorded for 2.25 h. (E) Racemization of glutamate in the D L direction. RacE1 and RacE2 were assessed for the capacity to catalyze the forward (L-glutamate D-glutamate) and reverse (D-glutamate L-glutamate) reactions by CD spectroscopy. L-Glutamate (5 mM) or D-glutamate (5 mM) was incubated in the presence of RacE1 (0.31 µM) or RacE2 (0.31 µM), and the CD signal was recorded for 2.25 h. For panels A to E, at least three separate experiments were performed. For each independent experiment, we used RacE1 or RacE2 from one of three independent enzyme preparations, as well as assay reagents from one of three independent preparations. For panels A to E, representative data from a single experiment are shown. In panel C, the molecular weights are reported as the means ± standard deviations from three independent experiments. mAU, milliabsorbance units.
|
-helix, ß-sheet, and ß-turn secondary structural elements. Moreover, comparison of the CD spectra revealed that the relative percentages of
-helix, ß-sheet, and ß-turn secondary structural elements were nearly identical for RacE1 and RacE2 (Fig. 2B and Table 3). These data suggest that when expressed as recombinant proteins, RacE1 and RacE2 had similar overall secondary structural properties. |
View this table: [in a new window] |
TABLE 3. Analysis of CD spectra using DICHROWEBa
|
Purified RacE1 and RacE2 are both functional in cell-free assays. To evaluate whether racE1 and racE2 encode functional glutamate racemase enzymes, RacE1 and RacE2 were assessed to determine their capacities to convert L-glutamate to the corresponding D enantiomer by using CD to directly observe the loss of L-glutamate as it was converted to D-glutamate. These experiments revealed that in the absence of RacE1 or RacE2, stereoisomerization of L-glutamate was not detectable (Fig. 2D). In contrast, consumption of L-glutamate was readily observed in the presence of either RacE1 or RacE2. Moreover, the rate of L-glutamate consumption increased as a function of RacE1 or RacE2 concentration. These data indicated that under cell-free and highly defined conditions, RacE1 and RacE2 both catalyzed the conversion of L-glutamate to D-glutamate in a concentration-dependent fashion. These results also established, for the first time, that despite reported differences in phenotypes of the null mutants (52), racE1 and racE2 both encode functional enzymes that have glutamate racemase activity.
RacE1 and RacE2 both catalyze the reverse reaction: conversion of D-glutamate to L-glutamate. In addition to catalyzing the conversion of L-glutamate to the corresponding D enantiomer, glutamate racemases also catalyze the reverse reaction, conversion of D-glutamate to the corresponding L enantiomer. To evaluate whether RacE1 and RacE2 share this canonical property of the glutamate racemase family, we used CD to directly measure the stereoisomerization of D-glutamate to the corresponding L enantiomer. The experiments revealed a time-dependent increase in the CD signal corresponding to the loss of D-glutamate in the presence of either RacE1 or RacE2 (Fig. 2E). Under the conditions of the reaction, the initial rate of stereoisomerization of D-glutamate was higher for RacE2 than for RacE1. In comparison, the rates of stereoisomerization of L-glutamate were similar for RacE2 and RacE1. Finally, these experiments revealed that for both RacE1 and RacE2, the rate of D-glutamate conversion is lower than the rate of L-glutamate conversion. These results indicated that RacE1 and RacE2 share one of the canonical properties of glutamate racemases, which is the capacity to catalyze the stereoisomerization of either glutamate enantiomer.
Steady-state kinetic analysis of RacE1 and RacE2.
To compare the catalytic properties of RacE1 and RacE2 in more detail, we analyzed the steady-state kinetic parameters of the two enzymes in the presence of D- or L-glutamate. RacE1 or RacE2 was incubated in a potassium borate buffer in the presence of various concentrations of D-glutamate or L-glutamate, and the change in magnitude of the CD signal was monitored. Initial rate data were obtained for RacE1 and RacE2 for each of the substrate concentrations, and plots of the rate of glutamate racemization (nmol/s) versus substrate concentration (mM) were generated for RacE1 and RacE2 in the presence of L- and D-glutamate (Fig. 3). Steady-state kinetic parameters for RacE1 and RacE2 in the presence of L- or D-glutamate were obtained by applying a nonlinear curve fit to the data (Table 4). In the L
D direction, differences in the individual kinetic parameters for RacE1 and RacE2 were statistically significant (for RacE1 kcat = 12 ± 0.6 s1 and for RacE2 kcat = 18 ± 0.6 s1 [P = 0.0003]; for RacE1 Km = 19 ± 4 mM and for RacE2 Km = 38 ± 6 mM [P = 0.01]). In the D
L direction, differences in the individual kinetic parameters for RacE1 and RacE2 were not statistically significant (for RacE1 kcat = 1.8 ± 0.1 s1 and for RacE2 kcat = 2.0 ± 0.1 s1 [P = 0.07]; for RacE1 Km = 1.0 ± 0.2 mM and for RacE2 Km = 0.77 ± 0.1 mM [P = 0.1]). These findings are in contrast to those for B. subtilis, in which there are approximately 100-fold differences in the catalytic efficiencies of RacE and YrpC (1, 2). Overall, these data indicated that RacE1 and RacE2 have similar, but not identical, steady-state kinetic properties under cell-free conditions.
![]() View larger version (16K): [in a new window] |
FIG. 3. Steady-state kinetic analysis of RacE1 and RacE2. RacE1 (0.78 µM) or RacE2 (0.78 µM) was incubated in a potassium borate buffer (50 mM boric acid, 100 mM KCl, 0.2 mM DTT; pH 8.0) in the presence of various concentrations of D-glutamate (0.1 to 10 mM) (A) or L-glutamate (5 to 200 mM) (B), and the change in magnitude of the CD signal was monitored. (A) Steady-state kinetic parameters for the racemization of glutamate in the L D direction by RacE1 and RacE2. The data are expressed as initial rate of racemization for RacE1 or RacE2 as a function of the L-glutamate concentration. Steady-state kinetic parameters for RacE1 and RacE2 in the presence of L-glutamate were obtained by applying a nonlinear curve fit to the data. (B) Steady-state kinetic parameters for the racemization of glutamate in the D L direction by RacE1 and RacE2. The data are expressed as the initial rate of racemization for RacE1 and RacE2 as a function of the D-glutamate concentration. Steady-state kinetic parameters for RacE1 or RacE2 in the presence of D-glutamate were obtained by applying a nonlinear least-squares regression utilizing GraphPad Prisim V4.03. For all studies, at least three independent experiments were performed. For each independent experiment, we used RacE1 or RacE2 from one of three independent enzyme preparations, as well as assay reagents from one of three independent preparations. The symbols indicate the means of the data from three independent experiments, and the error bars indicate the standard deviations of the means.
|
|
View this table: [in a new window] |
TABLE 4. Kinetic parameters for RacE1 and RacE2
|
D direction catalyzed by RacE1 or RacE2 are both pH dependent, with pH optima of 8.2 and 8.1, respectively. The kcat-versus-pH data fit a bell-shaped curve (equation 1) representative of an enzyme utilizing two ionizable side chains for catalysis (Fig. 4):
![]() | (1) |
![]() View larger version (9K): [in a new window] |
FIG. 4. pH dependence of RacE1 and RacE2 catalysis. RacE1 and RacE2 were assayed in buffers having pH values ranging from 6.5 to 9.5 with increments of 0.5 pH unit. Each buffer was prepared with a concentration of glutamate (200 mM) that was fully saturating for RacE1 and nearly saturating (83% saturating) for RacE2, so that the initial rate data would report true kcat values. Initial rate data were then obtained for RacE1 (0.78 µM) and RacE2 (0.78 µM) in each of the different buffer formulations by measuring the change in magnitude of the CD signal over time. The data are expressed as the turnover number (kcat) as a function of the reaction pH. The symbols indicate the means of the data from three independent experiments. For each independent experiment, we used RacE1 or RacE2 from one of three independent enzyme preparations, as well as assay reagents from one of three independent preparations. The error bars indicate standard deviations of the means.
|
root mean square deviation (RMSD) values of 7.5 Å and 2.6, respectively, indicating that RacE2 may share significant structural similarities with B. subtilis RacE, while RacE1 may possess a more divergent structural arrangement. The homology models for both RacE1 and RacE2 revealed that two cysteine residues (C77 and C188 for RacE1 and C74 and C185 for RacE2) are predicted to be in close proximity to the predicted location of the D-glutamate substrate (Fig. 5A), suggesting that these two cysteine residues may be important for the enzymatic activities of RacE1 or RacE2.
![]() View larger version (52K): [in a new window] |
FIG. 5. Site-directed mutagenesis reveals active site residues important for catalysis in both RacE1 and RacE2. (A) Three-dimensional homology models. Homology models for RacE1 and RacE2 were constructed using the Chemical Computing Group's MOE 2006.08. The template for both models was the B. subtilis glutamate racemase-D-glutamate structure (Protein Data Bank accession no. 1ZUW), which was aligned with the sequences for B. anthracis RacE1 (BAS0806) and RacE2 (BAS4379) using the Blosum62 substitution matrix. (B) Demonstration of residues important for catalysis. The two putative catalytic cysteine residues in RacE1 (Cys77 and Cys188) and RacE2 (Cys74 and Cys185) were independently changed to alanine by site-directed mutagenesis and assessed for the capacity to racemize glutamate in the L D direction. The four mutant enzymes (0.31 µM) or two wild-type enzymes (0.31 µM) were independently incubated in the presence of L-glutamate (5 mM). The differential absorption of CD light by glutamate was constantly monitored for 2.25 h utilizing a J-720 CD spectropolarimeter. (C) Chymotrypsin sensitivity patterns. Chymotrypsin protease sensitivity patterns were generated for the wild-type and mutant forms of RacE1 and RacE2. Wild-type RacE1 (0.13 mM), RacE1 C77A (0.13 mM), and RacE1 C188A (0.13 mM) were incubated in a Tris buffer (50 mM Tris-HCl, 100 mM NaCl, 2 mM DTT; pH 8.0) with various concentrations of chymotrypsin (lane A, 0 µg/ml; lane B, 53 µg/ml; lane C, 213 µg/ml; lane D, 640 µg/ml) at 4°C. Wild-type RacE2 (0.13 mM), RacE2 C74A (0.13 mM), and RacE2 C185A (0.13 mM) were incubated in a Tris buffer (50 mM Tris-HCl, 100 mM NaCl, 2 mM DTT; pH 8.0) with various concentrations of chymotrypsin (lane A, 0 µg/ml; lane B, 10 µg/ml; lane C, 40 µg/ml; lane D, 160 µg/ml) at 4°C. After incubation for 1 h, the reactions were stopped by addition of SDS sample buffer, and the samples were electrophoresed on a 16% SDS-polyacrylamide gel and stained with Coomassie brilliant blue G-250. The experiments in panels B and C were performed three separate times. For each independent experiment, we used wild-type or mutant forms of RacE1 or RacE2 from one of three independent enzyme preparations, as well as assay reagents from one of three independent preparations. Representative data from a single experiment are shown. WT, wild type.
|
|
|
|---|
A recent study reported a significantly more severe growth defect for a B. anthracis racE2 deletion mutant than for a racE1 deletion mutant (52), suggesting that there are fundamental differences between these two genes. However, it was not demonstrated that the growth defects were due specifically to reduced stereoisomerization of L-glutamate to D-glutamate by the mutant bacilli. Moreover, at the time of this study, neither racE1 nor racE2 had been demonstrated to encode functional glutamate racemases. The overall objective of the current study was therefore twofold: to establish if either racE1 or racE2 encodes a functional glutamate racemase and, if so, to explore whether there are fundamental differences in the biochemical properties between RacE1 and RacE2, which could provide insights into the growth phenotype differences between the racE1 and racE2 mutant strains (52). Notably, the two glutamate racemase isoenzymes in B. subtilis (RacE and YrpC) were demonstrated to have markedly different catalytic properties, with RacE exhibiting 100-fold-higher catalytic efficiency than YrpC (1, 2). In addition, racE is an essential gene, while yrpC is dispensable for growth in rich medium, supporting the idea that in B. subtilis only one of two glutamate racemases is necessary for converting L-glutamate to the D-glutamate required for rapid proliferation (31).
Our data indicated that racE1 and racE2 both encode functional glutamate racemases. In addition, we demonstrated that in a highly defined, cell-free system, RacE1 and RacE2 both catalyze the stereoisomerization of glutamate. In contrast to the disparate properties demonstrated for B. subtilis RacE and YrpC (1, 2), steady-state kinetic analysis identified B. anthracis as the first organism harboring genes encoding two glutamate racemases that in cell-free assays have similar, although not identical, catalytic properties. Analysis of the pH dependence of L-glutamate racemization suggested that RacE1 and RacE2 each possess two active site residues with titratable side chains having nearly identical pKa values. Furthermore, based on homology models that we generated for RacE1 and RacE2, we used directed mutagenesis to demonstrate the importance of two cysteine residues predicted to be in the active sites of both enzymes. Taken together, these results suggested that RacE1 and RacE2 may share similar active site features. However, a more detailed comparison of the predicted active site geometries of the RacE1 and RacE2 homology models (based on the B. subtilis RacE-D-glutamate crystal structure) suggested that several conserved active site residues may be positioned differently relative to bound D-glutamate (Fig. 6A). For example, the distances from the nitrogen atom of D-glutamate to the analogous oxygen atoms of two conserved active site side chains varied for RacE1 (3.1 Å for T189 and 3.0 Å for S15) and RacE2 (3.4 Å for T186 and 2.6 Å for S12). Furthermore, the proximities of the
-carbon of D-glutamate to the analogous sulfur atoms of RacE1 C77 (3.1 Å) and RacE2 C74 (3.4 Å) were different.
![]() View larger version (55K): [in a new window] |
FIG. 6. Three-dimensional homology models reveal differences in RacE1 and RacE2 fine active site features. Three-dimensional homology models for RacE1 and RacE2 were constructed using the Chemical Computing Group's MOE 2006.08. The template for both models was the B. subtilis RacE-D-glutamate structure (Protein Data Bank accession no. 1ZUW), which was aligned with the sequences for B. anthracis RacE1 (BAS0806) and RacE2 (BAS4379) using the Blosum62 substitution matrix. (A) RacE1 and RacE2 exhibit differences in the spatial arrangement of active site residues. Residues predicted to be important for catalysis in B. anthracis RacE1 (green) and RacE2 (purple) and B. subtilis RacE (gray) were aligned. The RMSD of the superposition of C atoms between B. anthracis RacE1 and B. subtilis RacE is 7.5 Å, while the RMSD between B. anthracis RacE2 and B. subtilis RacE is 2.6 Å. (B and C) Differences in inhibitor docking to the active sites of RacE1 and RacE2. The glutamate racemase inhibitor compound 69 (12) is shown docked within the active site of B. anthracis RacE1 (B) and overlaid with the active site of B. anthracis RacE2 (C). Residues that are present at the entrance to the hydrophobic binding pocket of B. anthracis RacE1 (green) and RacE2 (purple) are shown. Compound 69 was unable to dock within the active site of RacE2, presumably due to the larger side chain of RacE2 V149, which aligns with A152 in RacE1.
|
The results of these studies suggest that differences in the racE1 and racE2 phenotypes reported earlier are unlikely to be due solely to differences between the intrinsic catalytic efficiencies of the two enzymes (52). However, we cannot rule out the possibility that the catalytic properties of RacE1 and RacE2 may be substantially different within the bacterium. Several other possible levels of regulation of racE1 and racE2 inside the organism may contribute to the observed phenotypic differences in the respective deletion mutants. Transcript levels of virulence factors including the PDGA capsule are regulated by the global transcriptional regulator atxA (8, 11, 32, 60). The increased production of PDGA capsule that occurs in response to CO2 induction may lead to a larger demand for cytoplasmic levels of D-glutamate; thus, it is plausible to speculate that racE1 and racE2 could be differentially regulated at the transcriptional level. In addition, the relative importance or roles of RacE1 and RacE2 may be associated with differential localization of the two proteins. Differential intracellular localization of proteins has been shown to be important and tightly coupled with the life cycle of B. subtilis (49-51). Putative localization signals were not computationally identified for RacE1 or RacE2 (data not shown). An alternative possibility is that the enzymatic activities of RacE1 and RacE2 are regulated at the posttranslational level. Indeed, the enzymatic activity of glutamate racemase from E. coli is regulated by a peptidoglycan precursor (15, 27), a function that raises the possibility that a metabolite or peptidoglycan precursor could regulate the enzymatic activity of RacE1 or RacE2. Finally, our experiments revealed that RacE1 and RacE2 had distinct properties in solution when they were analyzed by gel filtration. The source of the unexpected apparent molecular weight for RacE2 in solution (51.5 x 103 ± 5.2 x 103) (Fig. 2C) is currently unknown but could be one or more heretofore unrecognized factors, including polydispersity, the binding of cofactors, posttranslational modifications, etc. Moreover, the significance of the apparent quaternary structural differences between RacE1 and RacE2 in solution is unclear, as is whether these differences are associated with regulation of enzymatic activities within the cell. Notably, within the family of glutamate racemases, quaternary structure differences apparently exist, as enzymes from B. subtilis (YrpC) (1), B. pumilus (33), L. fermentum (13, 19), and Pediococcus pentosaceus (37) have been reported to be monomers, while a dimeric form has been reported for A. pyrophilus (29) and there are conflicting reports about the quaternary structure of the E. coli enzyme (15, 62). Furthermore, RacE from B. subtilis has been reported to exist in equilibrium between a monomer and a dimer (57). Studies are currently under way in our laboratory to further explore the cellular roles of RacE1 and RacE2 and to elucidate the basis of regulation of their cellular activities.
In summary, we have demonstrated that racE1 and racE2 encode functional glutamate racemases, suggesting that the enzymatic activities of both enzymes should be targeted for inhibition. Although we have demonstrated experimentally that RacE1 and RacE2 possess similar, but not identical, enzymatic properties and catalytic residues, several differences in active site features are predicted for RacE1 and RacE2, suggesting that both active sites need to be considered when drugs for effective inhibition of glutamate racemase activity in B. anthracis are designed.
Published ahead of print on 11 May 2007. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»