![[Feedback]](/icons/banner/feedbackACT.gif){kind=icon}
![[For Subscribers]](/icons/banner/subscriberACT.gif){kind=icon}
![[Archive]](/icons/banner/archiveACT.gif){kind=icon}
![[Search]](/icons/banner/searchACT.gif){kind=icon}
![[Contents]](/icons/banner/tocACT.gif){kind=icon}
Nicola Zamboni,2
Uwe Sauer,2 and
Stéphane Aymerich1*
Microbiologie et Génétique Moléculaire, INRA (UMR1238) CNRS (UMR2585) and INAP-G, F-78850 Thiverval-Grignon, France,1 Institute for Molecular Systems Biology, ETH Zurich, 8093 Zurich, Switzerland2
Received 31 January 2006/ Accepted 16 April 2006
 |
ABSTRACT |
|---|
 Top Abstract IntroductionMaterials and MethodsResultsDiscussionReferences |
|---|
|
INTRODUCTION |
|---|
|
TopAbstractIntroductionMaterials and MethodsResultsDiscussionReferences |
|---|
|
Bioinformatic analysis of the Bacillus subtilis genome has revealed the presence of four paralogous genes encoding putative malic enzymes: maeA (formerly ywkA), malS, ytsJ, and mleA. On the basis of protein sequence comparisons, these putative malic enzymes can be distributed into two groups, MaeA and MalS on the one hand and YtsJ and MleA on the other (8). In a previous study, we have demonstrated that MaeA possesses malic enzyme activity and exhibits a clear preference for NAD+ as a coenzyme. We have also shown that expression of maeA is induced by the presence of malate. This transcriptional regulation is mediated by the two-component system MalK-MalR (formerly YufL-YufM) (8). These results confirmed the original observation of the stimulation of the synthesis of an NADP/NAD-dependent malic enzyme activity in the presence of malate (7). Furthermore, the two-component system MalK-MalR has been shown to control the malate-dependent expression of two other genes, maeN and yflS, encoding malate transporters (32). No information is available concerning the other three putative malic enzyme genes, except that mleA has been predicted to encode a malolactic rather than a malic enzyme, and thus to catalyze the conversion of malate into lactate, on the basis of its putative cotranscription with mleN, which encodes a malate-lactate antiporter (34). The two malolactic enzymes from lactic acid bacteria (Lactococcus lactis and Oenococcus oeni) characterized to date clearly belong to the malic enzyme family (4).
The current study was designed to characterize the respective function of each putative malic enzyme gene of B. subtilis. The expression pattern of these genes was analyzed; the corresponding gene products were purified, and their enzymatic activities were characterized. To investigate the physiological role of these genes, single- and multiple-mutant strains were constructed and their phenotypes were studied under different growth conditions. These experiments have revealed a major and complex physiological role for ytsJ.
|
MATERIALS AND METHODS |
|---|
|
TopAbstractIntroductionMaterials and MethodsResultsDiscussionReferences |
|---|
(lac-pro) supE thi hsd5/F' traD36 proA+B+ lacIq lacZM15] was used as a general cloning host, and E. coli strain M15(pREP4) (QIAGEN) was used for overproduction of His6-MalS, His6-YtsJ, His6-MleA, and His6-MaeA. E. coli strains were grown in Luria-Bertani broth (LB) supplemented with antibiotics when necessary (ampicillin,100 mg/liter; kanamycin, 25 mg/liter). A conventional calcium shock procedure was used for transformation (28). The B. subtilis strains used in this work are listed in Table 1. Standard procedures were used to transform B. subtilis (1). B. subtilis strains were grown in LB except when stated otherwise. Antibiotics for selection were added at 5 mg/liter (chloramphenicol), 0.4 mg/liter (erythromycin), 100 mg/liter (spectinomycin), 5 mg/liter (kanamycin), or 0.2 mg/liter (phleomycin). Cultures were performed at 37°C. Growth tests were performed in C minimal medium (2) [K2HPO4 at 70 mM, KH2PO4 at 30 mM, (NH4)2SO4 at 25 mM, MgSO4 at 0.5 mM, MnSO4 at 0.01 mM, ferric ammonium citrate at 22 mg/liter] supplemented with tryptophan at 50 mg/liter and carbon sources (10 g/liter). The cultures were inoculated from precultures in the same medium (excepted for GM1622, which was precultured in C glucose). Cultures for ß-galactosidase assays were performed in CQTHC minimal medium (C minimal medium supplemented with tryptophan 50 at mg/liter, glutamine at 0.5 g/liter, and casein hydrolysate at 0.5 g/liter) supplemented with carbon sources (10 g/liter).
|
Strain construction. Fusions of ytsJ, malS, ytnP, and mleN promoter fragments to a promoterless lacZ gene were constructed with the vector pDG1661 (18). The different promoter fragments to be tested were obtained by PCR with primers GLO03 and GLO04 for ytsJ, GLO06 and GLO07 for malS, GLO05 and GLO07 for ytnP, and GLO08 and GLO09 for mleN (Table S1 in the supplemental material). These fragments were inserted between the EcoRI and BamHI sites of pDG1661 to generate plasmids pGL00, pGL01, pGL02, and pGL03, respectively. These plasmids were used to insert the reporter constructs into the B. subtilis amyE locus to generate strains GM1600, GM1601, GM1602, and GM1603, respectively (Table 1). Comparisons with lacZ fusions at the homologous locus (for ytsJ and ytnP-malS) indicated that the promoter fragments designed for construction of strains GM1600 and GM1602 contain all of the cis-acting elements for regulated expression of ytsJ and ytnP, respectively; this was confirmed by disruption tests and reverse transcription-PCR analyses for the ytsJ promoter-containing fragment (data not shown). The mleN-mleA promoter fragment used contains the 260 bp preceding the mleN start codon; it thus includes the last 126 bp of the upstream asnB gene and the transcription terminator of the asnA-asnB operon.
Strain GM1608 contains a disruption of the ytsJ open reading frame (ORF) obtained by the insertion of a pMUTIN derivative plasmid containing an internal fragment of the ytsJ ORF (pHV30253). The in vivo replacement by homologous recombination of a kanamycin resistance cassette for the lacZ and erythromycin resistance genes of pMUTIN was then obtained with pEC23, a plasmid specifically designed for this purpose and constructed by P. Stragier. Strain GM1632 was constructed by interruption of the mleA ORF by simple crossover integration of pGL12, a derivative of vector pJH101 (14) containing an internal fragment of the mleA ORF.
ß-Galactosidase assays. ß-Galactosidase activities were measured by the method of Miller (25) on cell extracts prepared by lysozyme treatment and centrifugation. Protein concentration was determined by the Bradford method (3) with the Bio-Rad protein assay solution (Bio-Rad). One unit of ß-galactosidase activity is defined as the amount of enzyme that produces 1 nmol of o-nitrophenol/min at 20°C. No endogenous ß-galactosidase activity in parental strain 168CA, which is known to express ß-galactosidase at low levels only during sporulation (5), was detected under the conditions used here.
Expression and purification of His6-tagged proteins. His6-tagged proteins were overexpressed and purified with the QIAexpress kit (QIAGEN). ytsJ, malS, and mleA coding sequences were integrated into E. coli expression vector pQE-30 with corresponding primers (Table S1 in the supplemental material) in order to produce fusion proteins with six histidines at the N terminus. The resulting plasmids, pEC05, pEC06, and pEC07, respectively, were transformed into E. coli strain M15(pREP4) to generate the strains used for production. E. coli strain GTDp4, previously constructed in our laboratory (8), allowed overexpression of His6-MaeA. All of the following procedures were performed as described in the QIAexpress manual. His6-tagged proteins were stored in a buffer containing Tris at 50 mM (pH 8), NaCl at 150 mM, dithiothreitol at 1 mM, and EDTA at 0.5 mM.
Determination of kinetic parameters of YtsJ, MalS, and MleA. YtsJ, MalS, and MleA malic enzyme activities were tested at 37°C by spectrophotometrically monitoring, at 340 nm, NAD(P)H formation during the reductive decarboxylation of malate [malate plus NAD(P) into pyruvate plus NAD(P)H]. The reaction mixture was composed of 25 nmol of protein in 1 ml of buffer containing divalent cations (Tris-HCl at 50 mM [pH 8], MgCl2 at 10 mM, MnCl2 at 10 mM, KCl at 50 mM, ß-mercaptoethanol at 10 mM). For determination of the Km for malate, both cofactors, NAD+ and NADP+, were added to the reaction mixture at 5 mM. For determination of the Km for each cofactor, the malate concentration was raised to 40 mM. For both cofactors, a molar extinction coefficient of 6.22 x 106 cm2 mol–1 was used for calculations. The catalytic constant (kcat) is expressed per subunit. The initial rate data of the reaction were fitted to the Michaelis-Menten relationship by least-squares analysis to determine Km and kcat. Demonstration of malic enzyme activity was assessed in a coupled reaction: at the end of the first reaction (end of NAD reduction), 40 U of strictly NAD-dependent lactate dehydrogenase from Bacillus stearothermophilus (Sigma) and 10 mM CaCl2 were added to the reaction mixture. The oxidation of NADH, if produced by the first conversion of malate to pyruvate, was then monitored. Demonstration of the ability of each protein to decarboxylate OAA was realized by monitoring spectrophotometrically, at 430 nm, the disappearance of OAA after addition of YtsJ, MalS, MleA, or MaeA (1 U is defined as the activity leading to a decrease of 1 U of optical density at 430 nm [OD430]/min).
Overexpression of homologous and heterologous malic enzymes and transhydrogenase systems in B. subtilis.
The ytsJ, maeA, malS, mleA, sfcA, maeB, udhA, and pntAB genes were overexpressed in vivo in the GM1608 strain. The ORF of each gene has been amplified by PCR with corresponding primers and inserted into the pND23 replicative expression vector (6) downstream of an isopropyl-ß-D-thiogalactopyranoside (IPTG)-inducible promoter and an optimized ribosome binding site, between the BamHI and SphI sites for ytsJ and malS, between the BamHI and XhoI sites for udhA, between the BamHI and EcoRI sites for maeB, and between the HindIII and BamHI sites for maeA, mleA, sfcA, and pntAB. Plasmids pGL04, pGL05, pGL06, pGL07, pGL08, pGL09, pGL10, and pGL11, containing the ytsJ, malS, mleA, maeA, sfcA, maeB, udhA, and pntAB ORFs, respectively, were transformed into B. subtilis strain GM1608 to generate strains GM1660, GM1661, GM1662, GM1663, GM1664, GM1665, GM1666, and GM1667, respectively. The N-terminal domain of maeB was amplified with primers lgO75 and lgO89 or with primers lgO75 and lgO74, providing fragments encoding the first 427 or 489 amino acids of MaeB (maeB1 or maeB2, respectively). The longest fragment ends at the internal SphI restriction site of maeB (nucleotides 1369 to 1374). These fragments were inserted between the BamHI and SphI sites of vector pND23, and the resulting plasmids were transformed into strain GM1608 to give strains GM1669 and GM1670 (Table 1). Furthermore, native vector pND23 was transformed into strain GM1608 to give control strain GM1658. Complementation tests were performed with C minimal medium containing malate (1%), and overexpression was induced by addition of IPTG (final concentration, 0.5 mM).
Assay of NADP-dependent malic enzyme activity in total cell soluble extracts. Cells pellets of strains GM1668, GM1669, GM1670, GM1659, and 168CA were obtained from 2-ml aliquots of B. subtilis cultures growing exponentially in C minimal medium containing malate (10 g/liter). Cells were resuspended in 80 µl of Tris HCl, pH 7.4, containing 0.025 mM MgCl2. Lysozyme (0.8 mg) was added, and cells were incubated at 37°C for 30 min and then centrifuged for 10 min at 4°C and 12,000 rpm. Supernatants were assayed for malic enzyme activity in the presence of NADP as described for purified proteins.
Metabolic flux ratio analysis. Cultures were grown in Erlenmeyer flasks containing 30 ml of M9 minimal medium (21) [Na2HPO4 · 7H2O at 12.8 g/liter, KH2PO4 at 3 g/liter, NaCl at 0.5 g/liter, NH4Cl at 1 g/liter, MgSO4 at 1 mM, CaCl2 at 0.1 mM] with either 1 g/liter [U-13C]glucose and 4 g/liter naturally labeled glucose or 5 g/liter [1-13C]glucose. Flasks were inoculated at 1:100 with exponentially growing precultures. Culture aliquots were harvested at an OD600 of 1 to 1.2 and then hydrolyzed and derivatized as described elsewhere (35). Derivatized amino acids were then analyzed on a series 8000 gas chromatograph combined with an MD 800 mass spectrometer (Fisons Instruments, Beverly, MA). The gas chromatography-mass spectrometry-derived mass isotope distributions of proteinogenic amino acids were corrected for naturally occurring isotopes (15). The corrected mass distributions were related to the in vivo metabolic activities by previously described algebraic equations and statistical data treatment, which quantified several ratios of fluxes through converging reactions and pathways to the synthesis of five intracellular metabolites (15, 16) with the software FiatFlux (36).
|
RESULTS |
|---|
|
TopAbstractIntroductionMaterials and MethodsResultsDiscussionReferences |
|---|
|
Once transported into the cells, malate can be converted into pyruvate by a malic enzyme(s) or can be utilized in the Krebs cycle by malate dehydrogenase conversion to OAA. To test whether the residual growth on malate of strain GM1621, with the four putative malic enzyme-encoding genes deleted, relied on malate dehydrogenase activity, an mdh mutation was introduced into GM1621, thus generating quintuple-mutant stain GM1622. The growth tests reported in Fig. 2B indicate that this strain is totally unable to grow in malate minimal medium. This result confirmed that only the two pathways involving malic enzyme activity on the one hand and malate dehydrogenase activity on the other hand for malate utilization exist in B. subtilis. We also noted that ytsJ inactivation reduced the growth rate of an mdh mutant strain on glucose (doubling time = 196 ± 8 min versus 124 ± 4 min; Fig. 2C), even when aspartate (0.05%) was added to the medium, a supplementation that almost completely suppresses the growth defect of an mdh mutant strain (22). Further introductions of maeA, mleA, and malS mutations very slightly increased the growth defect in this medium (doubling time = 217 ± 10 min). These results confirmed that YtsJ has a major role that is probably distinct from that of the other malic enzymes MaeA, MalS, and MleA.
Expression pattern of the maeA, malS, ytsJ, and mleA genes. In order to gain information on the physiological role of the malic enzyme-like genes, their expression pattern was investigated. The expression of maeA, which had been studied in our previous work, is specifically induced in the presence of malate via the two-component system MalK-MalR (formerly YufL-YufM) (8). To obtain a kinetic description of the expression of these genes, lacZ transcriptional reporter gene fusions were constructed and inserted into the chromosome at the amyE ectopic locus. Measurements of ß-galactosidase activity synthesized by strain GM1600 (amyE'::PytsJ::lacZ) during growth on glucose or malate minimal medium revealed fairly constant expression from the beginning of the exponential phase until the beginning of the stationary phase of growth, with a twofold higher level on malate than on glucose (Fig. 3A). By contrast, in complex medium, expression dropped at the end of the exponential phase of growth. The malS ORF starts 85 bp downstream of the stop codon of the ytnP ORF, and no putative terminator sequence can be identified within this intergenic region, suggesting cotranscription of the two ORFs. Indeed, reverse transcription-PCR experiments clearly demonstrated that mRNAs overlapping ytnP and malS are synthesized (data not shown). Furthermore, a fusion between the ca. 500 bp preceding the malS ORF with the lacZ reporter gene resulted in barely detectable ß-galactosidase activity (Fig. 3B). In contrast, the fusion between the ytnP promoter region and lacZ in strain GM1602 (Fig. 3C) resulted in moderate and rather constant ß-galactosidase activity at a similar level on glucose or malate minimal medium or in complex medium.
|
YtsJ, MalS, and MleA exhibit malic enzyme activity, like MaeA. Our previous study had demonstrated that MaeA exhibits malic enzyme activity and preferentially utilizes NAD as a cofactor. MleA and MalS have been predicted to exhibit malolactic and malic enzyme activities, respectively (24, 34). To determine whether the three paralogous genes ytsJ, malS, and mleA do encode malic enzymes, the corresponding proteins were purified and their specific enzymatic activities were characterized. Vectors designed to overexpress the His-tagged YtsJ, MalS, and MleA proteins were introduced into E. coli strain M15(pREP4). Overexpression was induced by addition of IPTG, and the His6-YtsJ, His6-MalS, and His6-MleA proteins were purified to apparent homogeneity as described in Materials and Methods (MaeA was also purified and tested again here for more accurate comparisons; the quantitative results were similar to those reported in reference 8). The capacity of these proteins to reduce NAD+ or NADP+ in the presence of malate was then assayed. These tests demonstrated that YtsJ, MalS, and MleA exhibit a strong ability to reduce either NAD+ or NADP+ that strictly depends on the presence of malate. To demonstrate that this reducing activity is associated with the conversion of malate into pyruvate, i.e., that the proteins exhibit malic enzyme activity and not malate dehydrogenase or malolactic enzyme activity, a secondary assay with lactate dehydrogenase, an enzyme specifically active on pyruvate, was coupled to the first assay. Indeed, addition of lactate dehydrogenase after the end of the first reaction (no further increase in OD340) with either YtsJ, MalS, or MleA led to rapid reoxidation of NADH (data not shown). This could result only from oxidative conversion of pyruvate into lactate and thus demonstrated that pyruvate is the product of the first reaction. Kinetic parameters of YtsJ, MalS, and MleA malic enzyme activities (Table 2) were then determined as described in Materials and Methods. The Km values of YtsJ and MleA for malate were very similar (1.55 and 1.56 mM, respectively) and about twofold lower than those of MalS and MaeA (3.52 and 3.95 mM, respectively). Km, kcat, and catalytic efficiency (kcat/Km) for NAD+ and NADP+ were calculated for the four proteins. No activity in the presence of NADP+ could be detected for MleA. This malic enzyme is thus strictly NAD dependent. The kcat/Km values of MalS, 90 s–1 mM–1 for NAD and 9 s–1 mM–1 for NADP+, were very close to those of MaeA (30 s–1 mM–1 and 8 s–1 mM–1, respectively). Thus, MalS is, like MaeA, a malic enzyme with a dual specificity for NAD and NADP but with a strong preference for NAD. By contrast, the catalytic efficiency of YtsJ was much higher with NADP+ than with NAD+ (Table 2). Thus, YtsJ is also a malic enzyme with a dual specificity for NAD and NADP but with a very strong preference for NADP.
|
Occurrence of malic enzyme genes in bacteria.
Having established that B. subtilis possesses four paralogous genes encoding malic enzyme activity, we wondered whether such a situation is common in other bacterial species. We therefore determined how many putative malic enzyme genes are present in each sequenced bacterial genome. According to pairwise sequence comparisons, two groups of genes could clearly be distinguished: a ytsJ-like group and a maeA-like group (Fig. S1 in the supplemental material). Our analysis revealed that most of the sequenced species contain at least one malic enzyme-like gene (Table S2 in the supplemental material). Of these, most have at least one ytsJ-like gene. Only a few species have three or more malic enzyme paralogs (essentially bacilli and ß- or 
-proteobacteria).
A more detailed analysis of the sequenced Firmicutes species was then performed. Only three other Bacillus species or strains were found also to possess four malic enzyme genes. Four classes could be distinguished on the basis of sequence similarities and chromosomal context. ytsJ-like genes linked to dnaE, as in B. subtilis, form class 1. ytsJ-like genes linked to the malK-malR two-component genes (which control maeA induction by malate in B. subtilis) form class 2. ytsJ-like genes linked neither to dnaE nor to malK-malR form class 3. maeA-like genes linked to an aspartate lyase gene (as is a ytsJ-like gene, mleA, in B. subtilis) form class 4. The malic enzyme genes present in the Firmicutes species were assigned to one of these classes (Table S3 in the supplemental material).
Three main conclusions come from these analyses. (i) The bacterial species that possess at least one malic enzyme gene nearly always possess a ytsJ-like gene (175 of 196 species). (ii) Exchanges between the two groups of malic enzyme genes are frequently detected in members of the division Firmicutes (whereas a maeA-like gene is linked to malK-malR in B. subtilis and B. licheniformis, a ytsJ-like gene is associated with these two-component genes in other bacilli and in Enterococcus faecalis; an maeA-like gene, and not a ytsJ-like gene as in B. subtilis, is linked to an aspartate lyase gene in B. anthracis, B. thuringiensis, and Streptococcus mutans). (iii) Very few bacterial species (7 of 196) possess four putative malic enzyme genes, and only B. licheniformis, a species very closely related to B. subtilis, possesses four putative malic enzyme genes with the same genetic organization as in B. subtilis.
In vivo metabolic fluxes in malic enzyme mutants. Previous 13C-based flux analyses revealed small but significant in vivo malic enzyme fluxes even during growth on glucose, when they are not required (38, 37). To identify the responsible isoform, we grew malic enzyme mutants on 13C-labeled glucose. Intracellular ratios of metabolic fluxes were then quantified from the mass isotope distribution in protein-bound amino acids (15). In vivo malic enzyme activity is quantified by the molar fraction of pyruvate originating from malate, where the remaining fraction originates from pyruvate kinase (Table 3). Since the labeling pattern of malate is not directly accessible but inferred from glutamate (synthesized from 2-oxoglutarate) and aspartate (synthesized from OAA), the exchange between malate and OAA cannot be quantified. Hence, lower and upper boundaries for the fraction of pyruvate originating from malate are given which correspond to the extreme cases of absent or complete equilibration of the malate and OAA pools, respectively (15).
|
YtsJ cannot be substituted for by the E. coli NADP-dependent malic enzyme for efficient growth on malate. Our phenotypic analysis showed that a ytsJ deletion cannot be compensated for by the presence of the three other malic enzyme-encoding genes for malate utilization. Three hypotheses, which are not mutually exclusive, can be proposed to explain this phenomenon: (i) none of the other malic enzyme genes is sufficiently expressed under the conditions of the test, (ii) the unique NADP preference of YtsJ is necessary for fully efficient malate utilization, and (iii) YtsJ has, besides the catalysis of malate decarboxylation into pyruvate, another role necessary for malate utilization that cannot be fulfilled by any of the other three malic enzymes.
To test the first hypothesis, we constructed replicative vectors designed to overexpress malS, mleA, or maeA under the control of the IPTG-inducible Pspac promoter and introduced them into strain GM1608, generating GM1661, GM1662, or GM1663, respectively. As a control, the ytsJ gene was introduced into the same expression vector and the resulting plasmid was used to transform GM1608, creating GM1660. The growth rate of the transformed strains was measured in malate minimal medium in the presence of IPTG and compared to that of wild-type strain 168CA containing the empty replicative expression vector (GM1659). Strain GM1660 showed a roughly wild-type growth rate (Fig. 4), while GM1661 to GM1663 showed the same growth defect as the parental strain (data not shown). This indicated that overexpression of maeA, malS, or mleA cannot compensate for the ytsJ deletion and that the phenotype of the ytsJ mutant is indeed caused by the absence of YtsJ synthesis and not by any polar effect of the deletion.
|
Overexpression of the E. coli UdhA transhydrogenase partially suppresses the growth defect of a B. subtilis ytsJ mutant strain. The above results did not indicate an important role for malic enzyme-produced NADPH per se but left open a potential contribution of YtsJ to the maintenance of the NADPH-NADH balance. In the B. subtilis genome, no putative transhydrogenase-encoding gene has been identified nor is any other enzymatic system(s) known that could fulfill a similar function, i.e., interconversion of NADH and NADPH. To test this hypothesis, the two E. coli transhydrogenases, PntAB and UdhA, that catalyze the production of NADPH from NADH and the opposite reaction, respectively (29), were overexpressed from pND23 in the ytsJ mutant GM1608. While pntAB overexpression did not affect the growth rate of the ytsJ mutant strain, overexpression of udhA (GM1666) partially suppressed the growth defect on malate (Fig. 5). On malate, the doubling times of the wild type, the ytsJ mutant, and its UdhA-overexpressing derivative were 65, 130, and 105 (± 6) min, respectively. Overexpression of udhA partially relieved the metabolic defect caused by the ytsJ knockout but had no effect on the growth rate of a wild-type strain (data not shown).
|
|
DISCUSSION |
|---|
|
TopAbstractIntroductionMaterials and MethodsResultsDiscussionReferences |
|---|
All of the malic enzymes from eubacteria and eukaryotes characterized to date belong to the same protein family, together with the bacterial malolactic enzymes. Numerous biochemical and structural studies of eukaryotic malic enzymes (reviewed in reference 4) have allowed not only description of the catalytic mechanism but also identification of the residues involved in substrate and metal ion binding, as well as in catalysis. These residues are highly conserved in the four malic enzymes of B. subtilis (Fig. S1 in the supplemental material). Such extensive studies have not yet allowed unequivocal identification of the residues involved in preferential binding of NAD or NADP. However a careful examination of the different biochemical and structural data available led Chang and Tong (4) to identify two residues (Ser346 and Lys362 in the pigeon NADP-dependent malic enzyme) that would be directly responsible for NADP specificity. It is interesting that these residues are conserved in YtsJ but not in the other three malic enzymes of B. subtilis (Fig. S1 in the supplemental material), in agreement with the cofactor preferences of these enzymes characterized in this study.
What is the physiological function of malic enzymes? While expressed during growth on malate, they are not strictly necessary because malate can also be utilized via OAA conversion to PEP by using malate dehydrogenase and PEP-carboxykinase (PckA). Indeed, an E. coli double-mutant strain with both malic enzymes deleted is clearly but moderately affected in growth on malate (33). Instead, malic enzyme activity is assumed to increase the efficiency of malate utilization by contributing to a correct balance among the three key metabolites OAA, PEP, and pyruvate and, in the case of NADP-dependent malic enzyme activity, by producing the main source of reducing power for anabolism, NADPH (30).
Our phenotypic analysis clearly established that YtsJ plays an important role in efficient growth on malate that cannot be fulfilled by any of the other three malic enzymes, even when these are artificially overexpressed. This indicates distinct and synergistic roles for YtsJ compared to MalS, YwkA, and MleA (which can substitute for each other; Fig. 2). The four Km values for malate were very similar (Table 2) and thus can hardly explain these distinct roles. Similarly, the catalytic efficiencies (with NADP for YtsJ and with NAD for MaeA and MalS) were in the same range. Since YtsJ was the only NADP-dependent enzyme, we reasoned that its singular physiological role could be linked to this property. Overexpression of NADP-dependent MaeB from E. coli, however, did not complement the ytsJ B. subtilis mutant phenotype, despite wild-type-like NADP-dependent malic enzyme activity synthesized in the transformed strain. This result suggests that the role of YtsJ during growth on malate relies not only on its NADP-dependent malic enzyme activity. Nevertheless, we cannot exclude the possibility that the absence of complementation is due to kinetic parameters and/or to allosteric effectors of MaeB that are significantly different from those of YtsJ.
Much to our surprise, overexpression of the soluble E. coli transhydrogenase UdhA partially suppressed the growth defect of the ytsJ knockout on malate. The sole function of this enzyme is fully reversible conversion of NADPH and NADH. During growth of E. coli on glucose, the physiological role of UdhA is conversion of NADPH to NADH, while the membrane-associated transhydrogenase PntAB catalyzes the opposite reaction (29). During growth of the ytsJ mutant on malate, however, we expect massive production of NADH because all of the malate must be catabolized through either NADH-producing malate dehydrogenase or the NADH-producing malic enzyme and pyruvate dehydrogenase. Thus, it is not unreasonable to hypothesize that the cofactor concentrations are different from those during growth on glucose, such that UdhA could indeed catalyze NADPH production from NADH. The fact that pntAB overexpression does not complement the ytsJ phenotype—although it should, on the basis of its function in E. coli—might be related to its incorrect folding or the inability to form an active membrane-associated complex upon heterologous expression.
Since the partial complementation of UdhA overexpression can only be explained by its known involvement in cofactor balancing, we assume that at least one important physiological role of YtsJ is related to NADPH metabolism. This hypothesis, however, remains to be verified by cofactor concentration and detailed flux analyses. The current data suggest that YtsJ might have another function that is independent from the NADP-dependent conversion of malate to pyruvate, and this would explain why UdhA complements only partially. This additional function is only necessary under conditions with high malate-to-pyruvate fluxes, i.e., growth on malate, fumarate, or succinate, or in mdh mutants but not during growth on glucose.
The large majority of sequenced bacterial species have at least one malic enzyme-like gene, and most of these species have at least one ytsJ-like gene (Table S2 in the supplemental material). This higher conservation of ytsJ, which is particularly clear among bacilli, staphylococci, and clostridia, might indicate that a YtsJ-like enzyme also has a major physiological role in species other than B. subtilis. A second ytsJ homolog is present in bacilli other than B. subtilis and B. licheniformis which is linked to malK-malR two-component genes, whereas a maeA gene is associated with malK-malR in B. subtilis and B. licheniformis (Table S3 in the supplemental material). This exchange between the two different groups of malic enzyme genes suggests that a ytsJ-like gene can fulfill the specific function associated with MalK-MalR (while a maeA-like gene is unable to back up the major function of ytsJ in B. subtilis). Finally, B. anthracis, B. subtilis, and B. licheniformis appear to be singular in having four malic enzyme paralogs (together with only Streptomyces avermitilis, Colwellia psychrerythraea, and Photobacterium profundum of more than 300 species considered) but precise selective advantages conferred by this property and ensuring the conservation of this multigenic family remain to be clarified.
Whereas maeA transcription is specifically induced in the presence of malate in the growth medium (8), the ytsJ gene was transcribed at a high level under all of the conditions tested. This transcriptional pattern of ytsJ is in accordance with the major physiological role of the corresponding enzyme for efficient utilization of different C4 dicarboxylates, and not only malate, for growth but also more generally in central carbon metabolism (Table 3). On the contrary, the very weak transcription of malS or mleA in glucose or malate minimal medium suggests that they could have a physiological significance under specific conditions not explored in this study. The spectacular increase in mleA transcription in complex medium only during the exponential phase of growth could illustrate such specific conditions. The identification of the parameter responsible for high mleA transcription will give insights into the physiological conditions under which mleA is important. Our enzymatic tests showed that, in contrast to what has been proposed (34), MleA is not a malolactic enzyme; actually, a malolactic enzyme converts malate into lactate without release of pyruvate. It should also be noted that the two bacterial malolactic enzymes demonstrated enzymatically have a sequence much more similar to that of MaeA or MalS than to that of MleA or YtsJ. However, MleA could have a particularly significant role in acting together with lactate dehydrogenase to catalyze a coupled reaction equivalent to a malolactic enzyme reaction. This could explain its cotranscription with mleN, which encodes a malate-lactate antiporter (34).
In conclusion, B. subtilis possesses four paralogous genes encoding malic enzyme activity that are not physiologically equivalent. On the contrary, YtsJ exhibits a singular specificity for NADP as a cofactor and plays a major physiological role for which the three other malic enzyme cannot compensate. Our results indicate that this major role would be related to NADPH metabolism, but complementary experiments are necessary to determine its precise contribution to this metabolism.
|
ACKNOWLEDGMENTS |
|---|
Guillaume Lerondel is on a fellowship from the French Ministère de la Recherche et de la Technologie. Nicola Zamboni was supported by DSM Nutritional Products. Thierry Doan was and Stéphane Aymerich is on the staff at INRA. This work was supported by EU Quality of Life and Management of Living Resources program grant QL GZ-CT-1999-O1455.
|
FOOTNOTES |
|---|
Supplemental material for this article may be found at http://jb.asm.org/.
Present address: Department of Microbiology, Harvard Medical School, Boston, MA 02115.
|
REFERENCES |
|---|
|
TopAbstractIntroductionMaterials and MethodsResultsDiscussionReferences |
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
; GM1608 (ytsJ), 
; GM1621 (ytsJ maeA malS mleA), 
; GM1655 (maeA malS mleA), •; GM1614 (mdh), 
; GM1615 (ytsJ mdh), x; GM1622 (ytsJ maeA malS mleA mdh), 
; GM1617 (ytsJ maeA malS mdh), 
. Experiments were repeated at least three times, and the results of one representative experiment are shown.
