![]()
|
|
||||||||
J. Biol. Chem., Vol. 281, Issue 37, 27378-27388, September 15, 2006
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1
2






3
From the
Laboratory of Microbiology, Department of Agrotechnology and Food Sciences, Wageningen University, Hesselink van Suchtelenweg 4, 6703 CT Wageningen, the Netherlands, the
Biological and Environmental Systems Group, Department of Chemical and Process Engineering, University of Sheffield, Mappin Street, Sheffield S1 3JD, United Kingdom, the ¶Department of Biotechnology, Royal Institute of Technology (KTH), SE-106 91 Stockholm, Sweden, the ||Department of Molecular Evolution, Evolutionary Biology Center, Uppsala University, Norbyvägen 18C, SE-752 36 Uppsala, Sweden, the **Department of Biomolecular Mass Spectrometry, Bijvoet Center for Biomolecular Research and Utrecht Institute for Pharmaceutical Sciences, Utrecht University, Sorbonnelaan 16, 3584 CA Utrecht, the Netherlands, and the 
Laboratoire de Chimie Bioorganique et Bioinorganique, CNRS-UMR 8182, Institut de Chimie Moléculaire et des Matériaux d'Orsay (ICMMO), Université de Paris-Sud, Bâtiment 420, F-91405 Orsay, France
Received for publication, June 9, 2006 , and in revised form, July 17, 2006.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Sulfolobus spp. are obligatory aerobic Crenarchaea that are commonly found in acidic geothermal springs. Among the Archaea, this genus is well known for its broad saccharolytic capacity, which is reflected in their ability to utilize several pentoses and hexoses, as well as oligosaccharides and polysaccharides as a sole carbon and energy source (17). Although the catabolism of hexoses is well studied (reviewed in Ref. 18), the pathways for pentose degradation have neither been established in Sulfolobus solfataricus, nor in any other member of the Archaea (19).
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Growth of Sulfolobus Species
S. solfataricus P2 (DSM1617) was grown in media containing either 3 g/liter D-Ara or D-Glu as previously described (20).
Transcriptomics
Whole genome DNA microarrays containing gene-specific tags representing >90% of the S. solfataricus P2 genes (21) were used for global transcript profiling of cultures grown on D-Ara as compared with D-Glu. Total RNA extraction, cDNA synthesis and labeling, hybridization, and scanning were performed as previously described, as were data filtration, normalization, and statistical evaluation (22, 23).
Quantitative Proteomics
The proteome of S. solfataricus P2 was studied with a combination of two-dimensional gel electrophoresis, 15N metabolic labeling, and tandem mass spectrometry as previously described (24, 25). Two separate growth experiments were set up: 1) S. solfataricus with D-Ara as the carbon source and (14NH4)2SO4 as the nitrogen source; and 2) S. solfataricus with D-Glu as the carbon source and (15NH4)2SO4 as the nitrogen source. Next, the 14N and 15N cultures were mixed in equal amounts on the basis of optical density (A530) measurements, proteins were extracted and separated by two-dimensional gel electrophoresis. For the localization of proteins, a previously described two-dimensional gel electrophoresis reference map was used (23). Spots were excised from the gel, and peptides were quantified on the basis of their relative intensity in the time of flight mass spectrum, according to established methods (23).
Synthesis of Organic Compounds
D-Arabinonate was synthesized from D-Ara as previously described (26). The aldehyde 2,5-dioxopentanoate was synthesized from 1,4-dibromobutane according to reported procedures (27-29).
Gene Cloning and Protein Overexpression
The genes araDH (Sso1300), araD (Sso3124), kdaD (Sso3118), and dopDH (Sso3117) were amplified by PCR from genomic DNA using Pfu TURBO polymerase (Stratagene) and cloned in expression vector pET24d (Novagen) (supplemental materials Table 1). The resulting plasmids were harvested from Escherichia coli HB101 transformants by Miniprep (Qiagen), sequenced by Westburg genomics (Leusden, Netherlands), and transformed to E. coli expression strain BL21(DE3) containing the tRNA accessory plasmid pRIL (Stratagene).
All proteins were produced according to standard procedures in four 1-liter shaker flasks containing LB medium, but with some exceptions. When the culture A600 reached 0.5, the cultures were cold-shocked by placing them on ice for 30 min to induce host chaperones (20). After that, the expression was started by adding 0.5 mM isopropyl
-D-thiogalactopyranoside, and the cultures were incubated for 12-16 h at 37 °C after which they were spun down (10 min, 5000 x g, 4 °C). At the time of induction, the arabinose dehydrogenase (AraDH) and AraD overexpression cultures were supplemented with 0.25 mM ZnSO4 (30) and 20 mM MgCl2, respectively.
Protein Purification
Pelleted E. coli and S. solfataricus cells were resuspended in buffer and disrupted by sonication at 0 °C. Afterward, insoluble cell material was spun down (30 min, 26,500 x g, 4 °C) and the E. coli supernatants were subjected to heat treatment for 30 min at 75 °C. Denatured proteins were removed by centrifugation (30 min, 26,500 x g, 4 °C) yielding the heat-stable cell-free extract (HSCFE).
AraDHHSCFE in 20 mM Tris-HCl (pH 7.5) supplemented with 50 mM NaCl was applied to a 20-ml Matrex Red A affinity column (Amicon). After washing the bound protein with 2 column volumes of buffer, the recombinant protein was eluted by a linear gradient of 2 M NaCl.
AraDHSCFE in 50 mM HEPES-KOH (pH 8.0) supplemented with 50 mM NaCl was applied to a 70-ml Q-Sepharose Fast Flow (Amersham Biosciences) anion exchange column, and eluted in a 2 M NaCl gradient. Fractions containing the recombinant protein were pooled, concentrated with a 30-kDa cut-off filter (Vivaspin), and purified by size exclusion chromatography using a Superdex 200 HR 10/30 column (Amersham Biosciences) and 50 mM HEPES-KOH buffer (pH 8.0) supplemented with 100 mM NaCl as an eluent.
2-Keto-3-deoxy-D-arabinonate Dehydratase (KdaD)HSCFE in 25 mM NaPi buffer (pH 6.8) was applied to a 70-ml Q-Sepharose Fast Flow (Amersham Biosciences) anion exchange column. Flow-through fractions containing KdaD were collected, loaded onto a 46-ml Bio-Gel hydroxyapatite column (Bio-Rad), and eluted by a linear gradient of 0.5 M NaPi buffer (pH 6.8). Fractions containing the recombinant proteins were pooled, and dialyzed overnight in 50 mM HEPES-KOH (pH 8.0) supplemented with 0.5 mM dithiothreitol (DTT).
2,5-Dioxopentanoate Dehydrogenase (DopDH)HSCFE in 20 mM HEPES-KOH (pH 8.0) supplemented with 200 mM NaCl and 7.5 mM DTT was purified by affinity chromatography, as described for AraDH. Fractions containing the protein were pooled, concentrated using a 30-kDa cut-off membrane (Vivaspin), and purified by size exclusion chromatography as described for AraD.
Enzyme Assays
Unless stated otherwise, all enzymatic assays were performed in degassed 100 mM HEPES-KOH buffer (pH 7.5) at 70 °C. The optimal pH of catalysis was determined using a 30 mM citratephosphate-glycine buffer system that was adjusted in the range of pH 3-11 at 70 °C. Thermal inactivation assays were performed by incubating 50 µg/ml of enzyme at 70, 80, 85, and 90 °C and drawing aliquots at regular intervals during 2 h followed by a standard activity assay.
Dehydrogenase Assays
Sugar dehydrogenase activity was determined on a Hitachi U-1500 spectrophotometer in a continuous assay using 10 mM D- and L-arabinose, D- and L-xylose, D-ribose, D-lyxose, D- and L-fucose, D- and L-galactose, D-mannose, and D-glucose as a substrate, and 0.4 mM NAD+ or NADP+ as a cofactor. Aldehyde dehydrogenase reactions were performed using 5 mM 2,5-dioxopentanoate, glycolaldehyde, DL-glyceraldehyde, acetaldehyde, and propionaldehyde in the presence of 10 mM DTT. Initial enzymatic activity rates were obtained from the increase in absorption at 340 nm (A340), and calculated using a molar extinction coefficient of 6.22 mM-1 cm-1.
Dehydratase Assay
Standard reactions were performed using 10 mM potassium D-arabinonate in the presence of 1 mM MgCl2. The formation of 2-keto-3-deoxy-acid reaction products was determined with the thiobarbiturate assay at 549 nm using a molar extinction coefficient of 67.8 mM-1 cm-1 (31, 32). The effect of different divalent cations on enzymatic activity was investigated by a pre-treatment of the enzyme with 1 mM EDTA for 20 min at 70 °C, followed by a standard assay in the presence of 2 mM divalent metal ions.
Formation of 2-Oxoglutarate and Pyruvate
Enzyme reactions were performed with cell-free extract (CFE) from S. solfataricus cultures grown on either D-Ara or D-Glu, which were harvested at mid-exponential phase. The reaction was started by adding 25 µl of 3.5 mg/ml CFE to a mixture containing 10 mM potassium D-arabinonate, 1 mM MgCl2, and either 0.4 mM NAD+ or NADP+. After an incubation of 2 h at 75°C, the reactions were stopped by placing the tubes on ice. Identical reactions were set up in which the CFE was replaced by the purified enzymes AraD (4.2 µg), KdaD (13.4 µg), and DopDH (3.8 µg). The amount of 2-oxoglutarate in these mixtures was then determined by the reductive amination of 2-oxoglutarate to L-glutamate using purified recombinant Pyrococcus furiosus glutamate dehydrogenase at 60 °C (33). The detection reaction was started by the addition of 5 units of glutamate dehydrogenase to a sample that was supplemented with 10 mM NH4Cl and 0.12 mM NADPH. The formation of pyruvate was determined at 30 °C using 4 units of chicken heart lactate dehydrogenase and 0.1 mM NADH. The conversion of 2-oxoglutarate or pyruvate was continuously monitored on a Hitachi U-1500 spectrophotometer by following the decrease in A340 until substrate depletion occurred. Changes in concentrations of NAD(P)H were calculated as described above.
Determination of the Protein Oligomeric State
The oligomerization state of AraDH, AraD, KdaD, and DopDH was determined by nanoflow electrospray ionization mass spectrometry. For this, the protein was concentrated in the range of 5-15 µM and the buffer was exchanged to 50 or 200 mM ammonium acetate (pH 6.7 or 7.5) by using an Ultrafree 0.5-ml centrifugal filter device with a 5-kDa cut-off (Millipore). Protein samples were introduced into the nanoflow electrospray ionization source of a Micromass LCT mass spectrometer (Waters), which was modified for high mass operation and set in positive ion mode. Mass spectrometry experiments were performed under conditions of increased pressure in the interface region between the sample and extraction cone of 8 mbar by reducing the pumping capacity of the rotary pump (34, 35). Capillary and sample cone voltages were optimized for the different proteins and were in the range of 1.4-1.6 kV and 75-150 V, respectively.
Bioinformatic Analyses
Upstream sequences of the differentially expressed genes were extracted between -200 and +50 nucleotides relative to the open reading frame translation start site. These sequences were analyzed using the Gibbs Recursive Sampler algorithm (36). Possible sequence motifs were checked against all upstream sequences and the complete genome of S. solfataricus. A diagram of the sequence motif was created using the WebLogo server.
Protein sequences were retrieved from the National Center for Biotechnology Information (NCBI) and analyzed using PSI-BLAST on the non-redundant data base, and RPS-BLAST on the conserved domain data base. Multiple sequence alignments were built using either ClustalX or TCoffee software. Gene neighborhood analyses were performed using various webserver tools: STRING at the EMBL, Gene Ortholog Neighborhoods at the Integrated Microbial Genomes server of the Joint Genome Institute, and pinned regions at the ERGO bioinformatics suite.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
2 x 109 cells/ml (A600 2.5) with similar doubling times of around 6 h. Several enzyme activity assays were performed with CFEs from both cultures to establish a mode of D-Ara degradation (Fig. 1). A 12.3-fold higher NADP+-dependent D-Ara dehydrogenase activity (45.5 milliunits/mg) was detected in D-Ara CFE (Table 1), which indicated the presence of an inducible D-Ara dehydrogenase. D-Ara reductase, D-arabinitol dehydrogenase, and D-Ara isomerase activity were not detected. Activity assays using D-arabinonate indicated that D-Ara CFE contained a 13.9-fold higher D-arabinonate dehydratase activity (7.4 milliunits/mg) than D-Glu CFE (Table 1). Moreover, the multistep conversion of D-arabinonate to 2-OG could readily be demonstrated with D-Ara CFE in the presence of NADP+ (Fig. 2). The formation of pyruvate as one of the products from D-arabinonate was not observed, whereas control reactions with both CFEs and D-gluconate as a substrate did yield pyruvate (data not shown), indicating that the enzymes of the Entner-Doudoroff pathway were operative. In the final step of the pathway, D-Ara CFE contained a 3.6-fold higher activity (255 milliunits/mg) toward the aldehyde 2,5-dioxopentanoate (DOP) using NADP+ as a cofactor. The data suggest that S. solfataricus employs an inducible enzyme set that converts D-Ara into the tricarboxylic acid cycle intermediate 2-OG via the pentose oxidation pathway (Fig. 1C).
|
|
|
|
ProteomicsProtein expression in the soluble proteomes of D-Ara and D-Glu grown S. solfataricus cells was compared using a combination of two-dimensional gel electrophoresis, stable isotope labeling, and tandem mass spectrometry. By employing this strategy, five proteins were found with more than a 20-fold difference in expression level (supplemental materials Fig. 1, B-F), including the Ara-binding protein from the Ara ABC transporter (AraS, Sso3066) (39), Sso1300, Sso3124, Sso3118, and Sso3117 (Table 2). Interestingly, the difference in expression level of these genes at the protein level appeared to be more pronounced than at the mRNA level, which ranged from 3.4- to 16-fold. Three other proteins were also produced in higher amounts during growth on D-Ara, albeit only up to a 3-fold difference (Table 2). These were the isocitrate lyase (Sso1333) (45), the phosphoglycerate kinase (Sso0527) (46), and the malic enzyme (Sso2869) (47).
|
|
AraDHThe putative zinc-containing, medium-chain alcohol dehydrogenase encoded by Sso1300 was efficiently produced and purified using a single step of affinity chromatography (Fig. 4). The enzyme was most active on L-fucose (6-deoxy-L-galactose) (kcat 26.8 s-1), followed by D-Ara (kcat 23.8 s-1), using preferentially NADP+ (Km 0.04 ± 0.01 mM) over NAD+ (Km 1.25 ± 0.45 mM) as a cofactor. This enzyme was thus likely to account for the elevated D-Ara dehydrogenase activities in S. solfataricus CFE. AraDH could also oxidize L-galactose and the D-Ara C2-epimer D-ribose with similar rates (kcat 17.7 s-1) (Table 1). Enzyme activity toward other sugars remained below 7% of the highest activity. Similar substrate specificities and affinities have been found previously for mammalian and bacterial L-fucose or D-Ara dehydrogenases, although these enzymes prefer NAD+ as a cofactor (50, 51). AraDH was more than 50% active in a relatively narrow pH range from 7.3 to 9.3, with optimal catalysis proceeding at pH 8.2. The thermophilic nature of the enzyme is apparent from its optimal catalytic temperature of 91 °C. The enzyme maintained a half-life of 42 and 26 min at 85 and 90 °C, respectively, indicating that the enzyme is thermostable at physiological growth temperatures of S. solfataricus. Native mass spectrometry experiments showed that the intact recombinant AraDH has a molecular mass of 149,700 ± 24 Da. Comparing these data with the expected mass on the basis of the primary sequence (37,291 Da) clearly showed that the protein has a tetrameric structure and contains two zinc atoms per monomer. This is in good agreement with the tetrameric structure that has been reported for another alcohol dehydrogenase from S. solfataricus (Sso2536), which has a 33% identical protein sequence (30). This dehydrogenase, however, prefers aromatic or aliphatic alcohols as a substrate, and NAD+ over NADP+ as a cofactor. A structural study of AraDH is currently ongoing to explain the observed differences in substrate and cofactor selectivity.
Arabinonate Dehydratase (AraD)The protein encoded by gene Sso3124 was originally annotated as a member of the mandelate racemase and muconate lactonizing enzyme family. This superfamily, which additionally comprises of aldonic acid dehydratases, is mechanistically related by their common ability to abstract
-protons from carboxylic acids (52). Production of the enzyme in E. coli yielded
10% soluble recombinant protein, which was purified using anion exchange and size exclusion chromatography (Fig. 4). The enzyme was shown to catalyze the strictly Mg2+-dependent dehydration reaction of D-arabinonate to 2-keto-3-deoxy-D-arabinonate (KDA) (supplemental materials Fig. 2A). It is therefore conceivable that this enzyme is largely responsible for the increased levels of D-arabinonate dehydratase activity in S. solfataricus extracts. AraD displayed a maximum turnover rate of 1.8 s-1 at a substrate concentration of 8 mM, whereas higher substrate concentrations imposed severe inhibitory effects on the enzyme (supplemental materials Fig. 2B). No activity was measured with D-gluconate up to 20 mM. More than 50% enzyme activity was observed in a broad pH range of 5.2 to 10.2 with an optimum at pH 6.7 (Table 1). The enzyme was optimally active at 85 °C during which it maintained a half-life time of 18 min. Native mass spectrometry revealed that the protein had a molecular mass of 340,654 ± 63 Da, which corresponds well to an octameric protein assembly (expected monomeric mass is 42,437 Da). The native D-gluconate dehydratase from S. solfataricus (GnaD, Sso3198), which has a 23% identical protein sequence, was found to be an octamer as well (32). Interestingly, AraD was only produced as an octamer when the media was supplemented with 20 mM Mg2+ during protein overexpression. Without this divalent cation, the recombinant protein was inactive and appeared to be monomeric. Sequence alignment analysis as well as three-dimensional modeling based on a Agrobacterium tumefaciens protein with unknown function (Atu3453, Protein Data Bank code 1RVK) showed that Asp-199, Glu-225, and Glu-251 are likely to be involved in binding the divalent metal ion, which is required to stabilize the enolic reaction intermediate (52).
KdaDTo investigate the possible role of Sso3118, the protein was overproduced in E. coli, and subsequently purified (Fig. 4). Surprisingly, although the predicted pI of the enzyme is 5.9, the vast majority of protein did not bind to the anion exchange column at a pH of 8. Moreover, the protein had a tendency to precipitate, which could be reversed and effectively prevented by the addition of 0.5 mM DTT to all buffers. Native mass spectrometry under reducing conditions revealed that the protein had a molecular mass of 132,850 ± 47 Da, which corresponds with a tetrameric quaternary structure (expected monomeric mass of 33,143 Da). The catalytic activity of the protein was investigated by performing indirect enzyme assays using AraD with D-arabinonate as a substrate. A 50% decrease in the yield of KDA was observed when both enzymes were co-incubated in the presence of Mg2+, but this did not result in the formation of either 2-OG or pyruvate. Given the fact that D-arabinonate is converted to 2-OG in D-Ara CFE, this enzyme was anticipated to be responsible for the dehydration of D-KDA to the aldehyde DOP. However, due to the unavailability of D-KDA, it was not possible to show this in a direct enzyme assay. We therefore employed an indirect assay using AraD, the putative D-KDA dehydratase (KdaD) and the predicted aldehyde dehydrogenase. The results of this assay are described under "DopDH."
According to the Clusters of Orthologous Groups of proteins classification system, the putative KDA dehydratase belongs to COG3970. The catalytic domain of these proteins resembles that of the eukaryal fumarylacetoacetate hydrolase; an enzyme that catalyzes the Mg2+-orCa2+-dependent hydrolytic cleavage of fumarylacetoacetate to yield fumarate and acetoacetate as the final step of phenylalanine and tyrosine degradation (53). In humans, several mutations in the fumarylacetoacetate hydrolase gene will lead to hereditary tyrosinemia type I, which is mainly characterized by liver defects (54). Members of COG3970 are also homologous to the C-terminal decarboxylase domain of the bifunctional enzyme HpcE from E. coli, which in addition consists of an N-terminal isomerase domain (55). This enzyme is active in the homoprotocatechuate pathway of aromatic compounds and is responsible for the Mg2+-dependent decarboxylation of 2-oxo-5-carboxy-hept-3-ene-1,7-dioic acid to 2-hydroxy-hepta-2,4-diene-1,7-dioic acid and its subsequent isomerization to 2-oxo-hept-3-ene-1,7-dioic acid (55). Although the function of these enzyme classes is rather diverse, their structures have revealed similarities in terms of a fully conserved metal-ion binding site and a relatively conserved active site architecture. Multiple sequence alignment analysis of KdaD indicated the presence of a metal binding site consisting of Glu-143, Glu-145, and Asp-164, which may implicate a metal dependent activity as well. Further structural and kinetic studies of KdaD are currently ongoing.
DopDHThe putative aldehyde dehydrogenase encoded by Sso3117 was overproduced in E. coli, which resulted in the formation of
5% soluble protein. This protein fraction was purified using affinity and size exclusion chromatography (Fig. 4). From native mass spectrometry experiments we could determine a molecular mass of 210,110 Da, which is in reasonable agreement with the expected mass of the tetramer on the basis of the primary sequence (52,290 Da). The measured mass may be somewhat higher due to the binding of small molecules to the protein oligomer. The determined oligomerization state corresponds to that of the closely related aldehyde dehydrogenase ALDH-T from Geobacillus stearothermophilus (56). The aldehyde dehydrogenase was tested for the activity toward different aldehydes and cofactors (Table 1). This indicated that the enzyme preferred NADP+ over NAD+, and that it oxidized several hydrophilic aldehydes with the highest activity toward DOP followed by glycolaldehyde and DL-glyceraldehyde. More than 50% enzyme activity was observed in a pH range of 6.7-8.2, with an optimum at pH 7.8. The enzyme was also tested in conjunction with AraD and KdaD for the production of 2-OG or pyruvate. Similar to the activities in D-Ara CFE, these three enzymes were able to form 2-OG and not pyruvate, from D-arabinonate using preferably NADP+ as a cofactor (Fig. 2). Omission of either the cofactor, AraD, KdaD, or DopDH prevented the formation of 2-OG, indicating that all components were essential for the enzymatic conversions, and that KdaD was most likely responsible for the dehydration of D-KDA to DOP.
Extensive kinetic characterization of DopDH proved to be rather complicated, because the enzyme lost nearly all its activity within 1 day after its purification, even in the presence of high concentrations of reducing agents, such as DTT or
-mercaptoethanol. This could be due to the fact that this class of enzymes contains a catalytic cysteine residue (in DopDH Cys-293), which can become irreversibly oxidized, leading to a total loss of enzymatic activity. A rapid inactivation was also observed with ALDH-T from G. stearothermophilus (56).
Central Carbohydrate MetabolismSome central metabolic routes, such as the glycolysis, gluconeogenesis, and the tricarboxylic acid cycle have been studied extensively in S. solfataricus, Sulfolobus acidocaldarius, and other Archaea. The availability of their genome sequences (37, 57) as well as the genome sequence of Sulfolobus tokodaii (58), has recently allowed a reconstruction of the genes involved in these pathways (23). The effect of the introduction of excess 2-OG resulting from the D-Ara oxidative pathway led to the differential expression of only a few additional genes in these central carbon metabolic routes (Table 2; supplemental materials Fig. 3). The isocitrate lyase, the phosphoglycerate kinase, and the malic enzyme were up-regulated at the protein level under D-Ara conditions. The induction of the malic enzyme might indicate that the main proportion of 2-OG is converted to malate, which is then decarboxylated to pyruvate and acetyl-CoA, respectively, and is then fully oxidized to two molecules of CO2 in one round of the tricarboxylic acid cycle. Although this may seem energetically unfavorable, the net difference in yield between the full degradation of one molecule of D-Glu or D-Ara to CO2 is only one NADPH reduction equivalent in favor of D-Glu, because both degradation schemes lead to 6 reduced ferredoxins, 2 FADH2, 2 ATP, and 6 or 5 NADPH molecules, respectively. It is therefore not surprising that the growth rates under both conditions are similar. The phosphoglycerate kinase may be indicative of increased gluconeogenic activities that are required under D-Ara conditions. The isocitrate lyase is normally operative in the glyoxylate shunt, but high production levels of the enzyme have also been observed during growth on L-glutamate compared with D-Glu (25). Oxidative deamination of L-glutamate leads to the formation of 2-OG as well, which may inhibit the isocitrate dehydrogenase activity leading to an accumulation of isocitrate. This could trigger the production of the isocitrate lyase, which can bypass this step without the loss of CO2.
Pentose Oxidation Gene ClustersThe comprehensive analysis of conserved gene clustering in multiple genome sequences is becoming an increasingly important tool to predict functionally or physically associated proteins in prokaryotic cells (reviewed in Ref. 59). Genomic context analysis of the genes involved in the D-Ara oxidative pathway of S. solfataricus showed that kdaD and dopDH gene orthologs are often located adjacent in prokaryotic genomes. This finding supports the proposed enzymatic functions of an aldehyde producing and an aldehyde oxidative activity. In addition, the analysis uncovered the presence of putative pentose oxidative gene clusters in the genomes of several aerobic proteobacteria, such as members of the genera Burkholderia, Rhizobium, Bradyrhizobium, Agrobacterium, and Pseudomonas. In some cases, the presence of such a gene cluster correlates well with the ability of the organism to assimilate pentoses and with enzymatic activities present in cell extracts (7-11), whereas in other cases biochemical data is not available. Nonetheless, a few of these characteristic gene clusters have been demonstrated genetically to be linked to pentose degradation. Combined with the findings in S. solfataricus, this allows the identification of additional enzymatic components in the pentose oxidation pathway and prediction of their enzymatic functions (Fig. 5A).
A putative operon of five genes was found in the genome of the oligotrophic
-proteobacterium Caulobacter crescentus, which was 2.8-11.6-fold up-regulated during growth on D-Xyl as compared with D-Glu (60). Reporter fusion constructs of the CC0823 promoter to the
-galactosidase gene (lacZ) from E. coli confirmed that this promoter is highly induced during growth on D-Xyl, and repressed on D-Glu or proteinaceous media (60, 61). Moreover, the disruption of the CC0823 gene prevented the C. crescentus from growth on D-Xyl as a single carbon source (61).
A second pentose degradation gene cluster involved in L-Ara uptake and utilization was found on chromosome II of the pathogenic
-proteobacterium Burkholderia thailandensis. This cluster consisting of nine genes was proposed to be responsible for the L-Ara degradation to 2-OG (Fig. 5A) (62). Disruption of the araA, araC, araE, and araI genes led to an L-Ara negative phenotype. Reporter gene insertions showed that araC and araE gene expression was repressed during growth in D-Glu media, and induced in L-Ara media. The transfer of the gene cluster to the related bacterium B. pseudomallei enabled this organism to utilize L-Ara as a sole carbon source also (62). Interestingly, an L-Ara dehydrogenase with 80% sequence identity to AraE has recently been characterized from Azospirillum brasiliense (63); an organism that is known to degrade L-Ara to 2-OG (12). The flanking sequences of this gene revealed close homologs of the B. thailandensis araD and araE, which would indicate a similar gene cluster in A. brasiliense (63).
Apart from several bacteria, putative pentose oxidation clusters are also present in some Archaea. In the halophile Haloarcula marismortui, a gene cluster was found on chromosome I that seems to contain all of the necessary components for D-Xyl oxidation, including a gene that has been identified as a D-Xyl dehydrogenase (19) (Fig. 5A).
Components of the Pentose Oxidation PathwayCareful inspection of the different pentose oxidation gene clusters shows that the gene encoding the final enzymatic step, from DOP to 2-OG, is fully conserved between the different pentose oxidation gene clusters. The remaining analogous enzymatic steps that convert D-Ara, D-Xyl, or L-Ara into DOP are performed by enzymes from distinct COGs (Clusters of Orthologous Groups of proteins) (64) (Fig. 5, A and B, pentose panels). Whereas some of this variation in enzyme use may simply be explained by substrate differences, other variations may be due to the individual adaptation of existing enzymes with similar reaction chemistry, such as the pentose dehydrogenases.
|
-lactamases belonging to COG2220, which catalyze similar lactame-ring opening reactions (66). Other non-orthologous enzyme components of the pentose oxidation pathway include the pentonic acid dehydratases. Whereas the D-arabinonate dehydratase from S. solfataricus belongs to COG4948, the same function seems to be performed by members of COG0129 that are commonly annotated as dihydroxyacid dehydratases (IlvD) or 6-phosphogluconate dehydratases (Edd) (67). A member of this family has recently been characterized from S. solfataricus (DHAD, Sso3107), which revealed a broad substrate specificity for aldonic acids (68). However, this gene was not differentially expressed according to the transcriptome or proteome analysis.
The 2-keto-3-deoxy-D-arabinonate dehydratase (COG3970), or a member of the homologous COG0179, appears to be present in D-Ara and D-Xyl degradation gene clusters. Interestingly, in several Burkholderia species, and in A. brasiliense, this gene is replaced by a member of the dihydrodipicinolate synthase family (COG0329, B.th araD). Members of this family catalyze either aldolase or dehydratase reactions via a Schiff base-dependent reaction mechanism by a strictly conserved lysine residue. Interestingly, a detailed study of an L-KDA dehydratase involved in the L-Ara metabolism of P. saccharophila was reported a few decades ago, but unfortunately, neither the N-terminal sequence of the protein nor the gene sequence was determined (69, 70). The authors found that this enzyme was enantioselective for L-KDA (2-oxo-4(R), 5-dihydroxypentanoate), and that the reaction proceeds via a Schiff-base intermediate. The enzyme activity was not affected by the presence of 1 mM EDTA, which suggests a divalent metal-ion independent reaction. It seems likely that this enzyme is encoded by homologs of the B. thailandensis araD gene, and that the apparent enantioselectivity of this enzyme does not allow a function in the degradation of D-Ara or D-Xyl, which results in a 2-keto-3-deoxypentonic acid with the S-configuration (Fig. 5B).
The aldehyde dehydrogenase from COG1012 is fully conserved in the pentose oxidation gene clusters (Fig. 5A). Strikingly, close homologs of this gene can also be found in hexaric acid degradation gene clusters of Bacillus species (ycbC-ycbI) (71, 72) (Fig. 5A). The same holds for a gene cluster in Pseudomonas putida (PP1245-PP1249) that is likely to be involved in the breakdown of L-hydroxyproline, which is a major constituent of collagen and plant cell wall proteins (73, 74) (Fig. 5B). Apparently, because the degradation of both hexaric acids and L-hydroxyproline is also known to proceed through DOP (9), the genetic information for the conversion of DOP to 2-OG has been shared between multiple metabolic pathways during evolution (Fig. 5, A and B). Apart from the dopDH gene, orthologs of the D-glucarate dehydratase gene (ycbF, COG4948) are observed in the pentose degradation gene clusters of both S. solfataricus and H. marismortui, although remarkably, the keto-deoxy-acid dehydratase of COG0329 is found in all three pathways. In the hydroxyproline degradation pathway, this enzyme might function as a deaminase instead (75).
The apparent mosaic of orthologous and non-orthologous proteins involved in the pentose oxidation pathway suggests that some of these enzymatic steps may have evolved by recruitment events between enzymes from the hexaric acid or hydroxyproline degradation pathways, which also make use of DOP as an intermediate and produce 2-OG as the final product (76, 77). The low number of enzymes required, their common cofactor usage, and the large gain of obtaining the hub metabolite 2-OG as the end product of pentose oxidation, may have been the driving force in the creation of this pathway in aerobically respiring Bacteria and Archaea.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental Tables S1-S2, Figs. S1-S3, and additional references. ![]()
2 Supported by the University of Sheffield and the United Kingdom's Engineering and Physical Sciences Research Council. ![]()
3 Supported by Engineering and Physical Sciences Research Council, Advanced Research Fellowship GR/A11311/01 and Grant GR/S84347/01. ![]()
1 To whom correspondence should be addressed. Tel.: 31-317-483110; Fax: 31-317-483829; E-mail: stan.brouns{at}wur.nl.
4 The abbreviations used are: 2-OG, 2-oxoglutarate; HSCFE, heat-stable cell-free extract; CFE, cell-free extract; Ara, arabinose; KDA, 2-keto-3-deoxy-arabinonate; DOP, 2,5-dioxopentanoate; AraDH, arabinose dehydrogenase; AraD, arabinonate dehydratase; KdaD, 2-keto-3-deoxy-D-arabinonate dehydratase; DopDH, 2,5-dioxopentanoate dehydrogenase. ![]()
| REFERENCES |
|---|
|
|
|---|