Previous Article | Next Article ![]()
Applied and Environmental Microbiology, May 2008, p. 3229-3241, Vol. 74, No. 10
0099-2240/08/$08.00+0 doi:10.1128/AEM.02750-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Department of Chemical Engineering, University of California, Berkeley, California 94720-1462,1 Department of Bioengineering, University of California, Berkeley, California 94720-1762,2 Synthetic Biology Department, Physical Biosciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 947203
Received 6 December 2007/ Accepted 4 March 2008
|
|
|---|
9-18:1) in the growth medium. This work demonstrates the utility of using transcriptomic and metabolomic methods to optimize synthetic biological systems. |
|
|---|
Previously, Escherichia coli was engineered to produce amorpha-4,11-diene, the sesquiterpene precursor to the antimalarial compound artemisinin, using a multigene, heterologous pathway (44). Artemisinin, a sesquiterpene lactone endoperoxide extracted from Artemisia annua L (family Asteraceae; commonly known as sweet wormwood), is highly effective against multidrug-resistant Plasmodium spp. and is the key component of the artemisinin-based combination therapies now endorsed as the frontline therapy to treat the disease (2). However, the drug is in short supply and unaffordable to most malaria sufferers. Although total synthesis of artemisinin is difficult and costly, the semisynthesis of artemisinin or any derivative from microbially sourced artemisinic acid, its immediate precursor, represents a cost-effective, environmentally friendly, high-quality, and reliable source of artemisinin for use in artemisinin-based combination therapies. During the design optimization of the amorpha-4,11-diene pathway in E. coli, it was discovered that high expression of the first three genes in this synthetic pathway (atoB, ERG13, and tHMG1) resulted in a severe growth defect. The enzymes encoded by these three genes convert acetyl-coenzyme A (acetyl-CoA) to mevalonate via two intermediate compounds, acetoacetyl-CoA and 3-hydroxy-3-methylglutaryl-CoA (HMG-CoA) (Fig. 1A). While acetoacetyl-CoA is a compound endogenous to E. coli, HMG-CoA is not, and the growth inhibition was shown to correlate with the intracellular accumulation of this exogenous small molecule (51). However, the reason for its toxicity, and thus a solution to the problem, was not known.
![]() View larger version (27K): [in a new window] |
FIG. 1. (A) The heterologous pathway produces mevalonate from acetyl-CoA in three biochemical steps. The data presented in this article show that the growth inhibition associated with the pathway intermediate HMG-CoA is due to inhibition of malonyl-CoA:ACP transacylase (FabD) in E. coli's type II FAB pathway. The inhibition of the FAB pathway invokes a generalized stress response in the mevalonate-producing strain. As can be seen, the heterologous intermediate is similar in structure to malonyl-CoA, the native substrate for FabD. Abbreviations: A-CoA, acetyl-CoA; AA-CoA, acetoacetyl-CoA; HMG-CoA, hydroxymethylglutaryl-CoA; G6P, glucose-6-phosphate; FDP, fructose-1,6-bisphosphate; G3P, glyceraldehyde-3-phosphate; DHAP, dihydroxyacetone phosphate; PEP, phosphoenolpyruvate. (B) Genes of the heterologous mevalonate-producing pathway in E. coli as well as the inactive-pathway control, which due to a point mutation in ERG13, cannot produce HMG-CoA or mevalonate. Abbreviations: atoB, acetoacetyl-CoA thiolase gene; Erg13, HMG-CoA synthase gene; tHMGR, truncated HMG-CoA reductase gene; Erg13(C159A), inactive HMG-CoA synthase gene.
|
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Strains used in this study
|
The low-copy-number plasmid pBAD33MevT contains the previously constructed MevT operon under regulated control of the arabinose-inducible PBAD (44). The operon consists of the genes atoB from E. coli, encoding acetoacetyl-CoA thiolase, and ERG13 and tHMG1 from Saccharomyces cerevisiae, encoding HMG-CoA synthase and truncated HMG-CoA reductase 1, respectively, whose expression allows the in vivo conversion of acetyl-CoA to mevalonate.
In order to investigate the source of toxicity in cells expressing the MevT operon, an inactive-pathway operon was constructed by replacing the catalytic cysteine of S. cerevisiae HMGS (at amino acid position 159), encoded by pBAD33MevT, with an alanine (51). Expression of the modified operon produced active acetoacetyl-CoA thiolase, active HMG-CoA reductase, and full-length, catalytically inactive HMG-CoA synthase. The resulting plasmid was pMevT(C159A) (51).
Medium components and chemicals were purchased from Sigma-Aldrich (St. Louis, MO) and Fisher Scientific (Pittsburgh, PA). For propagation of engineered E. coli strains, Luria broth with Miller's modification (Sigma) was used with appropriate antibiotics for plasmid selection, and 0.06% glucose was used for the repression of the PBAD promoter system. Studies to characterize the cell physiology, metabolite levels, and mRNA transcript levels of the engineered strains were performed using defined C medium (30) supplemented with 3.4% glycerol, all individual amino acids per the work of Neidhardt et al. (46), 4.5 µg/ml thiamine-HCl, micronutrients, and 50 µg/ml (each) chloramphenicol and carbenicillin.
Fatty acids for medium supplementation were purchased from Sigma and employed in defined medium at a concentration of 100 µg/ml in the presence of 400 µg/ml Brij detergent (Sigma).
Growth, metabolite, and mRNA transcript analysis of engineered cells.
Cell growth, metabolite concentrations, and mRNA transcript levels of E. coli DP10 expressing genes encoding the active and inactive MevT operons were assayed. Starter cultures of E. coli DP10 harboring pBAD33MevT and pBAD18 (empty vector) or pMevT(C159A) and pBAD18 were inoculated from single colonies and incubated overnight at 37°C in defined C medium supplemented with 0.06% glucose (to repress PBAD) and antibiotics. Overnight starter cultures were diluted to an optical density at 600 nm (OD600) of 0.05 in fresh medium, incubated at 37°C with continuous shaking, and induced with the addition of 1.33 mM (0.02%) arabinose at an OD600 of
0.2 to 0.3. During experiments involving metabolite or transcript level analysis, the engineered strains were incubated in baffled shake flasks. Samples were taken at multiple time points during the course of the experiments to assay for cell growth, mevalonate production, transcripts, and intracellular metabolites. Samples for microarray analysis were taken preinduction and during logarithmic growth to capture the early response dynamics of the host to the HMG-CoA toxicity. The OD600 of shake flask samples was measured using a UV spectrophotometer (Beckman), and values for dry cell weight (DCW) were calculated using an equivalence factor of 1 g DCW/liter to an OD600 of 2.5 (4). Experiments were repeated in duplicate to confirm trends in growth and metabolite concentrations.
To screen for fatty acids that alleviate HMG-CoA-associated growth inhibition, cultures were grown in 96-well microtiter plates at 37°C, using defined C medium as described above. The optical density was measured using a microtiter plate reader (SpectraMax; Molecular Devices).
Transcript analysis sample preparation.
Biomass designated for transcript analysis was snap-frozen in liquid nitrogen immediately and stored at –80°C until analysis. RNAs were extracted from the cell samples by use of an RNeasy Midi kit (Qiagen). Using 40-µg aliquots of extracted RNA from each sample point, prelabeled cDNA was synthesized using random-primed reverse transcription in a 40-µl volume containing 12.5 µg primers (Invitrogen, Carlsbad, CA), 1x reverse transcriptase buffer (Invitrogen, Carlsbad, CA), 0.01 mM dithiothreitol (Invitrogen, Carlsbad, CA), 1 unit/µl Superase-In (Ambion, Austin, TX), a 0.5 mM concentration (each) of dATP, dCTP, and dGTP (Invitrogen, Carlsbad, CA), 0.1 mM dTTP (Invitrogen, Carlsbad, CA), 0.4 mM amino-allyl-dUTP (Ambion, Austin, TX), and 10 units/µl Superscript II (Invitrogen, Carlsbad, CA), following the enzyme manufacturer's instructions. The cDNA was base hydrolyzed in 100 mM NaOH-10 mM EDTA at 65°C for 10 minutes and then neutralized in HEPES, pH 7.0, at a final concentration of 500 mM. The Tris remaining in the cDNA suspension was removed by three buffer exchange spins using Micron YM-30 columns (Millipore) and eluted in a final volume of 15 µl of water. The cDNA was then labeled using either Alexa 555 or Alexa 647 (Invitrogen, Carlsbad, CA) following the manufacturer's protocol.
Microarray hybridization.
Glass microarrays printed with full-length double-stranded DNA prepared by PCR or 70-mer oligonucleotides (Operon) designed to probe every open reading frame of E. coli MG1655 were hybridized in a Tecan hybridization station with
6 to 10 µg of labeled cDNA per channel of detection. The hybridization program included a prehybridization step (5x SSC [1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate]-0.2% sodium dodecyl sulfate [SDS]-1% bovine serum albumin, 42°C, 60 min), a 15-hour hybridization step (Ambion Hyb solution 3, 40°C, medium agitation), two low-stringency washes (1x SSC-0.2% SDS, 42°C, 2 min [each]), two high-stringency washes (0.1x SSC-0.2% SDS, 25°C, 2 min [each]), and two final washes (0.1x SSC, 25°C, 2 min [each]). Following hybridization, the slides were scanned with an Axon 4500 scanner.
Transcriptional profile data analysis.
The raw scans were globally normalized using Genepix software and then exported to SNOMAD (9) for loess normalization to correct for any hybridization artifacts and to generate local Z scores (see reference 52 for a review of local Z score use in microarray analysis). The local Z scores generated by SNOMAD as well as by serial analysis for microarray (SAM) (60) software were used as guides to determine biologically significant gene expression changes among the replicate hybridization data sets. This list of significant genes was then mined using hierarchal clustering (13) (Cluster 3.0) to determine a base set of clusters in each data set. Once a base set of clusters was chosen, k-means clustering was also used to search the data set for temporal patterns in gene expression (Cluster 3.0).
GC-MS quantification of mevalonate.
The mevalonate (mevalonic acid) concentration in cultures of engineered E. coli was determined by gas chromatography-mass spectrometry (GC-MS) analysis as described previously (51). Briefly, culture aliquots were acidified with HCl in glass GC vials to convert mevalonate to mevalonic acid lactone. The acidified cultures were extracted with an equal volume of ethyl acetate containing (–)-trans-caryophyllene (internal standard), and the organic layer was diluted into fresh ethyl acetate prior to GC-MS analysis. Ethyl acetate extracts were analyzed using GC-MS by scanning for only ions m/z 71 and 58, corresponding to mevalonic acid lactone, and ions m/z 189 and 204, corresponding to (–)-trans-caryophyllene. The retention time, mass spectrum, and concentration of extracted mevalonic acid lactone were confirmed using commercial DL-mevalonic acid lactone (Sigma).
Intracellular metabolite extraction and analysis.
The concentrations of intracellular acyl-CoAs and adenylate pool were determined by liquid chromatography-MS (LC-MS) analysis of trichloroacetic acid (TCA) culture extracts as described previously (51). To simultaneously and rapidly quench cellular metabolism, isolate E. coli cells from growth medium, and extract metabolites, cells were centrifuged through a layer of silicone oil into a denser solution of TCA. TCA extracts, neutralized with tri-n-octylamine, were analyzed using an Agilent 1100 series LC-MS employing electrospray ionization and operating in positive mode. The following selected ions corresponding to the protonated molecular ion of each metabolite were monitored: for ATP, m/z 508; for ADP, m/z 428; for AMP, m/z 348; for CoA, m/z 768; for acetyl-CoA, m/z 810; for propionyl-CoA, m/z 824; for crotonyl-CoA, m/z 836; for acetoacetyl-CoA, m/z 852; for malonyl-CoA, m/z 854; for succinyl-CoA, m/z 868; for methylmalonyl-CoA, m/z 868; and for HMG-CoA, m/z 912. Retention times, mass spectra, and concentrations of extracted metabolites were confirmed using commercial standards (Sigma). The adenylate energy charge of each strain was calculated from the adenylate pool measurement as defined by Atkinson (3), as follows:
![]() |
Analysis of cellular fatty acid composition.
The fatty acid composition and the fatty acid fraction of engineered cells were determined by fatty acid methyl ester (FAME) analysis of lyophilized cell pellets (55). Sample aliquots of the bacterial cultures were centrifuged to pellet the cells, and the biomass was snap-frozen in liquid nitrogen. The frozen pellets were then lyophilized for 24 h at –80°C under vacuum, using a freeze drier. FAME analysis of lyophilized cell pellets was performed by Microbial ID (Newark, DE).
|
|
|---|
Malonyl-CoA accumulates during HMG-CoA-mediated growth inhibition.
The growth and acyl-CoA profiles of the active-pathway (mevalonate-producing) strain and the inactive-pathway control strain were analyzed by LC-MS to monitor metabolite changes that correlate with HMG-CoA growth inhibition (Fig. 2). In comparison to the inactive-pathway control, growth of E. coli expressing MevT from pBAD33MevT was inhibited (Fig. 2A), and the growth inhibition correlated directly with the accumulation of HMG-CoA (Fig. 2D), as previously reported (51). This pathway intermediate accumulated to
125 nmol/g DCW in the active-pathway strain during log-phase growth and did not accumulate in the inactive-pathway control strain (Fig. 2D).
![]() View larger version (26K): [in a new window] |
FIG. 2. Growth and metabolite time course profiles for inactive-pathway control strain [E. coli DP10 harboring pMevT(C159A) and pBAD18] ( ) and the growth-inhibited, mevalonate-producing strain (E. coli DP10 harboring pBAD33MevT and pBAD18) ( ). Arabinose induction was done at 0.1 g DCW/liter, and the first sample for metabolite extraction was taken at 0.5 h postinduction. Shown are cell growth, in g DCW/liter (A), mevalonate concentration (B), intracellular malonyl-CoA (C), intracellular HMG-CoA (D), intracellular free CoA (E), intracellular acetyl-CoA (F), intracellular acetoacetyl-CoA (G), and energy charge (H).
|
34% of its total acyl-CoA pool to malonyl-CoA, which represented a >4-fold increase in the intracellular concentration of this fatty acid precursor compared to that in the inactive-pathway control. Acetyl-CoA, an important central metabolite and precursor to mevalonate, was observed to be twofold lower in the active-pathway strain (Fig. 2F). This decrease in acetyl-CoA concentration following arabinose induction of the MevT genes can be explained partially by diversion of carbon into the engineered pathway. The anomalous accumulation of malonyl-CoA most likely played a role in the lower concentration of acetyl-CoA observed in the active-pathway strain as well, suggesting that the normal regulation of fatty acid metabolite levels was impaired by high concentrations of intracellular HMG-CoA. The accumulation of malonyl-CoA in E. coli DP10 containing pBAD33MevT may also account for the slightly lower level of acetoacetyl-CoA (Fig. 2G) than that in E. coli DP10 harboring pMevT(C159A). While it has been reported that low free CoA levels can inhibit protein production and growth (34), the free CoA concentration in the mevalonate-producing strain was similar to that in the inactive-pathway control strain (Fig. 2E).
The remaining acyl-CoAs that were tracked (propionyl-CoA and succinyl-CoA) did not vary significantly among strains. The adenylate energy charges of all cultures were >0.80 during early logarithmic-phase growth, indicating rapid quenching of cellular metabolism and negligible degradation of adenylates or acyl-CoAs during extraction from the cell (8).
FAB is altered by HMG-CoA accumulation.
DNA microarray analysis was performed on the growth-inhibited, mevalonate-producing strain and the inactive-pathway control strain to study the early transcriptional response to HMG-CoA toxicity. RNA was isolated from biomass sampled from each culture just prior to induction with arabinose and at 1 and 3 hours postinduction; the RNA was used to synthesize labeled cDNA. In order to fully elucidate the transcriptional response to HMG-CoA accumulation, the following three separate expression profiles (with technical replicates) were generated by two-color hybridization of the cDNA to DNA microarrays: profile A, the 1- and 3-hour-postinduction transcript profile for the inactive-pathway control [E. coli DP10 containing pMevT(C159A) and pBAD18] relative to a preinduction profile; profile B, the mevalonate-producing strain (E. coli DP10 containing pBAD33MevT and pBAD18) 1- and 3-hour-postinduction transcript profile relative to a preinduction profile; and profile C, the mevalonate-producing strain's 1- and 3-hour-postinduction transcript profile relative to the inactive-pathway strain's 1- and 3-hour-postinduction transcript profile.
There was a consistent up-regulation observed in the expression of the β-ketoacyl-ACP synthase I (encoded by fabB), both over the growth time course of the active-pathway strain (postinduction compared to preinduction) (data not shown) and in the cross-strain comparison (active-pathway strain versus inactive-pathway control strain) (Fig. 3). There was also an up-regulation in transcription of malonyl-CoA-ACP transacylase (encoded by fabD) observed in the
3-h postinduction sample, with a smaller up-regulation of fabH. Both fabD and fabH share a common promoter, though there can be extensive posttranscriptional processing of these transcripts (10). FabD is the only enzyme in E. coli that interacts with malonyl-CoA directly, and this occurs during the transfer of the malonyl moiety from CoA to the acyl carrier protein (ACP) (10). Additionally, there were several genes involved in the initial steps of FAB whose expression was up-regulated in mevalonate-producing E. coli, including those encoding biotin biosynthetic enzymes and acetyl-CoA carboxylase (Fig. 3).
![]() View larger version (31K): [in a new window] |
FIG. 3. Transcript profiles of the initial steps of type II FAB in Escherichia coli. Malonyl-CoA is synthesized from acetyl-CoA by the action of acetyl-CoA carboxylase, a heterotetramer composed of subunits encoded by accABCD. The malonate moiety is transferred from CoA to ACP by the action of malonyl-CoA:ACP transacylase (FabD). See the work of Magnuson et al. for a full review of the E. coli FAB pathway (43). Also shown (inset) are the expression values and Z scores (in parentheses) for FAB genes that exhibited a biologically significant up-regulation in the mevalonate-producing strain (E. coli DP10 containing pBAD33MevT and pBAD18) relative to the inactive-pathway control strain [E. coli DP10 containing pMevT(C159A) and pBAD18] in the microarray analysis. Also shown are the expression profiles for the CFA synthase gene (cfa) as well as two genes in the biotin synthesis pathway, the gene encoding 7,8-diaminopelargonic acid synthase (bioA) and one encoding biotin synthase (bioB).
|
![]() View larger version (26K): [in a new window] |
FIG. 4. Cellular fatty acid compositions of the mevalonate-producing strain (E. coli DP10 containing pBAD33MevT and pBAD18) and the inactive-pathway strain [E. coli DP10 containing pMevT(C159A) and pBAD18] obtained by FAME analysis. Open symbols in each graph represent the inactive-pathway strain, and solid symbols represent the mevalonate-producing strain. (A) Total fatty acids as a percentage of DCW. (B) UFA as a percentage of total fatty acids in each strain. and , 16-carbon UFA (includes 16-carbon CFA); and , 18-carbon UFA (includes 18-carbon CFA); and and , CFA (both 16-carbon and 18-carbon CFA). (C) SFA as a percentage of total fatty acids. and , 16-carbon SFA; and , 14-carbon SFA; and and , 12-carbon SFA.
|
The accumulation of HMG-CoA induces a cascade of stress responses.
E. coli has evolved several stress response regulons that allow the organism to adapt rapidly to environmental changes. The transcriptional modulation of these regulons (box plots in Fig. 5) provided insight into the toxicity associated with HMG-CoA accumulation in the mevalonate-producing strain. These box plots show the regulation of groups of related genes as a whole based upon EcoCyc (40) and gene ontology classification.
![]() View larger version (30K): [in a new window] |
FIG. 5. (A) Box plots of transcript expression Z scores for the genes in the 32 heat shock regulon, oxidative stress-regulated genes (OxyR and SoxS oxidative stress regulons), and osmotically induced operons at 1 hour and 3 hours postinduction. The box plots include the Z scores for all genes detected in the respective functional category based upon EcoCyc (40) annotation and gene ontology classification. The local Z score generated from the SNOMAD analysis represents a change value weighted to account for the confidence in the data, and any gene with a Z score with an absolute value of >1.96 was considered differentially expressed (9). The box plots show the values of the greatest down-regulation, lower quartile (Q1), median, upper quartile (Q3), and largest up-regulation observed for each profile for the regulon or gene group indicated. There were three expression profiles generated in this experiment, as follows: profile A, the 1- and 3-hour-postinduction transcript profile for the inactive-pathway control [E. coli DP10 containing pMevT(C159A) and pBAD18] relative to a preinduction profile; profile B, the mevalonate-producing strain (E. coli DP10 containing pBAD33MevT and pBAD18) 1- and 3-hour-postinduction transcript profile relative to a preinduction profile; and profile C, the mevalonate-producing strain's 1- and 3-hour-postinduction transcript profile relative to the inactive-pathway strain's 1- and 3-hour-postinduction transcript profile. A strong activation of the oxidative and osmotic stress regulons was observed in the active-pathway strain but not in the inactive-pathway control strain. The heat shock regulon was activated early in the active-pathway strain, but overall expression of the regulon was lowered at the 3-hour time point. The heat shock regulon remained highly activated in the inactive-pathway control strain. (B) Time course of expression of genes encoding the 30S and 50S ribosomal proteins as well as the ribosome modulation factor (rmf) observed in the cross-strain microarray analysis (mevalonate-producing strain versus the inactive-pathway control strain). The log2(expression ratio) is defined such that positive values (green) represent up-regulation and negative values (red) represent down-regulation in the mevalonate-producing strain compared to the inactive-pathway control. There was a significant down-regulation of ribosomal protein genes in E. coli DP10 harboring pBAD33MevT (active-pathway strain) compared to their expression in E. coli DP10 harboring pMevT(C159A) (the inactive-pathway control strain).
|
|
View this table: [in a new window] |
TABLE 2. Transcript expression ratios of genes with significant differential expression in the heat shock, osmotic stress, and oxidative stress responsesa
|
The expression of the mevalonate pathway also elicited a moderate heat shock response consistent with heterologous protein production. As shown in the time course expression profile box plots of
32-regulated genes in cells harboring the active pathway (Fig. 5A), expression of many members of the heat shock regulon was induced at 1 h postinduction, but then expression of most of these genes decreased at 3 h postinduction. The heat shock regulon genes most highly expressed in the active-pathway strain were those encoding the chaperones ClpB, DnaK, GroEL, and GroES as well as the inclusion body proteins IbpA and IbpB. In contrast, the inactive-pathway control continued to express genes in the heat shock regulon at high levels 3 hours after induction, and the same trend was observed when the two strain's transcriptional profiles were compared directly.
Coincident with expression of the heat shock regulon, there was a strong, two- to fivefold down-regulation of the genes encoding the ribosomal proteins observed in the active-pathway strain at 1 and 3 hours postinduction (time course clustering) (Fig. 5B). A strong heat shock response in E. coli is known to result in a down-regulation of the ribosomal protein genes (12, 53), yet the inactive-pathway strain, which had a much stronger up-regulation of the heat shock genes (Fig. 5A, profile A), did not exhibit a similar down-regulation of the ribosomal genes. Interestingly, the decreased transcription of the ribosomal protein genes in the active-pathway strain correlated with the fourfold up-regulation of the gene encoding the ribosome modulation factor, rmf. This gene, whose product dimerizes the 70S ribosome and reduces translational capacity, is generally expressed during the transition to stationary phase or under conditions of stress (33).
Thus, the two strains both exhibited an activation of the heat shock regulon, but each had a quite different response dynamic. There was a strong expression of this regulon in the active-pathway strain 1 hour after induction, but that response decayed by the third hour. The inactive-pathway control exhibited a steadily increasing induction of the heat shock genes. Heterologous protein production often increases the transcription of heat shock genes encoding molecular chaperones and proteases (12). Thus, the strong down-regulation of the ribosomal proteins, likely part of a coordinated response in the HMG-CoA-stressed cells, limited heterologous protein synthesis and, in turn, led to the general down-regulation of the heat shock regulon observed for the time course of cells harboring pBAD33MevT. Conversely, since the inactive-pathway strain was not stressed by HMG-CoA accumulation, it was able to maintain a higher rate of protein synthesis, which induced increased transcription of heat shock regulon genes, as observed at the later time point. This could explain the observation previously made that HMG-CoA reductase activity is higher in the inactive-pathway strain (51).
Medium supplementation tests.
Based on the differential fatty acid auxotrophies reported in the literature for E. coli strains with various mutations of FAB genes, we designed a panel of medium supplementation tests to screen for fatty acids that could relieve the HMG-CoA-associated growth defect. E. coli strain DP10 was transformed with either pBAD33MevT (active pathway), pMevT(C159A) (inactive-pathway control), or pBAD33 (empty vector control) and grown in a defined medium (C medium) in a 96-well plate growth assay. The growth medium was supplemented with 100 µg/ml of oleic acid (cis-
9-18:1), palmitoleic acid (cis-
9-16:1), palmitic acid (16:0), a combination (16:0 plus cis-
9-16:1), or no supplement at all. Both control strains grew well in all medium formulations tested (Fig. 6).
![]() View larger version (22K): [in a new window] |
FIG. 6. Effect of fatty acid supplementation on growth of mevalonate-producing and control strains of E. coli DP10. (A) Growth of the active-pathway strain [E. coli DP10(pBAD33MevT)] in a defined medium supplemented with 100 µg/ml of oleic acid (cis- 9-18:1) ( ), palmitoleic acid (cis- 9-16:1) ( ), palmitic acid (16:0) ( ), a combination of fatty acids (16:0 plus cis- 9-16:1) (+), or no fatty acid supplement ( ) in a 96-well incubator/plate reader monitored continuously at OD600. The growth of the inactive-pathway strain [E. coli DP10 carrying pMevT(C159A)] in nonsupplemented medium ( ) is shown as well. Growth profiles of the inactive-pathway control strain and an empty vector control [E. coli DP10(pBAD33)] were equivalent in medium supplemented with the various fatty acids. (B) Shake flask growth of E. coli DP10(pBAD33MevT) with () and without ( ) a palmitic acid (16:0) supplement compared to that of the inactive-pathway control [E. coli DP10 carrying pMevT(C159A)] ( ).
|
|
|
|---|
It has previously been reported that there is significant conversion of the acyl-CoA pool to malonyl-CoA in E. coli when the initial steps of FAB are inhibited by chemicals, such as the antibiotics thiolactomycin (28, 29) and cerrulenin (17, 29), or genetic methods, such as disruption of protein function (58) and mutation (36). E. coli coordinates FAB with lipid requirements (and hence growth), and the accumulation of long-chain acyl-ACPs is a key signal of slowing growth, which results in the coordinated down-regulation of fabB, fabA, and accAB (27, 28, 35, 36) as well as the inhibition of enzymatic activity (27). During this down-regulation of FAB, excess malonyl-CoA is recycled back to the acetyl-CoA pool through the actions of KAS1 (fabB) and malonyl-CoA-ACP transacylase (fabD) (28). Thus, the impressive accumulation of malonyl-CoA observed in the mevalonate-producing cells is an indication that this complex regulation has been disrupted.
Supplementation of the growth medium with certain fatty acids abolished the HMG-CoA-induced growth inhibition in the active-pathway strain by removing the need for fatty acid anabolism and confirmed the hypothesis that FAB was inhibited, either directly or indirectly, by HMG-CoA accumulation. The results from the entire panel of fatty acid supplements offer insight as to where HMG-CoA may inhibit the FAB pathway. E. coli requires both SFA and UFA to maintain proper membrane fluidity. Harder et al. observed that temperature-sensitive fabB mutants were auxotrophic for trans-UFA or a combination of SFA and cis-UFA at the nonpermissive temperature (25), while fabD mutants grew at the nonpermissive temperature only when the medium was supplemented with 16:0 SFA or 18:1 UFA but not 16:1 UFA (26). Additionally, when this fabD mutant was grown in nonsupplemented medium at a temperature that only partially inhibited growth, the membrane lipids were enriched for 16:1 SFAs. Since the supplementation of UFA or a combination of SFA and UFA failed to relieve the HMG-CoA-associated toxicity, FabB does not appear to be the target of HMG-CoA inhibition in the FAB pathway. Instead, our results suggest that FabD was the target, since inhibited growth of the active-pathway strain in nonsupplemented medium enriched UFA and the addition of palmitic acid (16:0) (and, to a lesser extent, oleic acid [cis-
9-18:1]) markedly improved the growth of this strain.
FabD is the only enzyme in E. coli that interacts with malonyl-CoA directly, and this occurs during the transfer of the malonyl moiety from CoA to ACP (10). This reaction is critical to FAB since it provides the two-carbon units needed to extend fatty acids. The reaction mechanism is mediated by an active-site serine that attacks the thioester carbonyl and releases free CoA, leaving a malonyl-FabD complex. Next, the malonyl-FabD ester carbonyl is attacked by the phosphopantetheinyl thiol of ACP, which ultimately yields malonyl-ACP (47). HMG-CoA is structurally very similar to malonyl-CoA and may interfere with substrate binding to the active site. Supporting this hypothesis is the evidence that another acyl-CoA, acetyl-CoA, has been shown to be a weak inhibitor of malonyl-CoA-ACP transacylase activity (38, 47).
The up-regulation of accBC, fabB, fabD, and fabH observed in the growth-inhibited, mevalonate-producing strain is consistent with the hypothesis that HMG-CoA inhibits FAB. This inhibition reduced the availability of long-chain acyl-ACPs, which was sensed by the heterologous host as a decoupling of the rates of FAB and growth. The resulting up-regulation of FAB genes was the host cell's attempt to maintain growth by up-regulating the flux of carbon into fatty acid anabolism. A similar transcriptional response to chemical inhibition of FAB has been documented for Mycobacterium tuberculosis, where exposure to thiolactomycin elicited an up-regulation of the genes encoding β-ketoacyl-ACP synthase I (fabB), malonyl-CoA-ACP transacylase (fabD), and alkyl hydroperoxide reductase C (ahpC) (7), just as observed in the mevalonate-producing strain.
The up-regulation of FAB genes allowed the heterologous host to overcome the most deleterious effect of HMG-CoA, a complete cessation of growth, which is similar in principle to the resistance to thiolactomycin that E. coli gains when fabB is expressed on a multicopy plasmid (11). Yet while the up-regulation of the FAB system was sufficient to allow slow growth in the presence of high levels of cytosolic HMG-CoA, there were consequences to the cell. Jackowski and Rock observed that inhibition of the initial steps of FAB increased the UFA/SFA ratio in the membrane (34) because the inhibition reduced the availability of fatty acids, which, in turn, induced expression of the β-ketoacyl-ACP synthase I (fabB). This enzyme has multiple roles in E. coli's FAB pathway, including both SFA synthesis and the diversion of saturated long-chain acyl-ACPs into the desaturation pathway that produces UFA. Indeed, the increase in FAB gene expression observed in the HMG-CoA-inhibited cells resulted in a coordinate enrichment of UFA in the cell membrane and an overall decrease in fatty acids as a percentage of DCW.
The altered UFA/SFA ratio in the membrane that resulted from HMG-CoA inhibition of fatty acid anabolism in the mevalonate-producing strain could explain the cascade of osmotic, oxidative, and heat shock stress responses observed in that strain's transcriptional profile. Recent studies of the osmosensor ProP in E. coli have discovered that membrane composition is a key signal governing the long-term response to alterations in medium osmolality (59). The oxidative stress response in the active-pathway strain was also an indication of membrane-associated stress, since H2O2 is mostly a by-product of the respiratory electron transport chain in E. coli (20), and both degradation of membrane integrity and osmotic stress have been reported to induce transcription of H2O2 defense genes (14, 57). Finally, the coordinated heat shock response to H2O2 is well established (15, 16), and the sum of all these stress responses appeared ultimately to result in a down-regulation of the translational machinery in the mevalonate-producing strain. Thus, the HMG-CoA-induced changes in membrane composition resulted in a coordinated modulation of many stress response regulons, which may be the ultimate cause for the slow growth observed in HMG-CoA-stressed cells.
In summary, the model for the growth inhibition observed in E. coli expressing the MevT operon is one where the cell's type II fatty acid anabolism was impeded, most likely by inhibition of FabD activity, which in turn limited the availability of long-chain acyl-ACPs. This reduced availability increased the transcription of genes in the initial steps of fatty acid synthesis, which overcame the HMG-CoA inhibition enough to allow slow growth. As a consequence of low levels of acyl-ACP and increased expression of FabB activity, however, the membrane lipids were enriched in UFA. This in turn altered the membrane structural properties and induced the transcription of genes associated with osmotic, oxidative, and heat shock stress. The sum of these responses ultimately limited the growth rate in the mevalonate-producing culture. It should be noted that a similar deleterious interaction does not appear to occur during accumulation of HMG-CoA in S. cerevisiae (56), which has a nondisassociated, type I FAB pathway, and interestingly, the type II pathway antibiotic thiolactomycin is not active against the S. cerevisiae type I FAB pathway. This report highlights the unexpected interactions that can occur when a novel, heterologous biochemical pathway is engineered into a host organism as well as the utility of a systems biology approach to the design problems that inevitably arise in metabolic engineering.
We acknowledge Jack Newman for his assistance.
Published ahead of print on 14 March 2008. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»