- Open Access
Metabolic engineering of Rhodopseudomonas palustris for the obligate reduction of n-butyrate to n-butanol
Biotechnology for Biofuelsvolume 10, Article number: 178 (2017)
Rhodopseudomonas palustris is a versatile microbe that encounters an innate redox imbalance while growing photoheterotrophically with reduced substrates. The resulting excess in reducing equivalents, together with ATP from photosynthesis, could be utilized to drive a wide range of bioconversions. The objective of this study was to genetically modify R. palustris to provide a pathway to reduce n-butyrate into n-butanol for maintaining redox balance.
Here, we constructed and expressed a plasmid-based pathway for n-butanol production from Clostridium acetobutylicum ATCC 824 in R. palustris. We maintained the environmental conditions in such a way that this pathway functioned as the obligate route to re-oxidize excess reducing equivalents, resulting in an innate selection pressure. The engineered strain of R. palustris grew under otherwise restrictive redox conditions and achieved concentrations of 1.5 mM n-butanol at a production rate of 0.03 g L−1 day−1 and a selectivity (i.e., products compared to the consumed substrate) of close to 40%. Since the theoretical maximum selectivity is 45%, the engineered strain converted close to its maximum selectivity.
The innate redox imbalance of R. palustris can be used to drive the reduction of n-butyrate into n-butanol after expression of a plasmid-based enzyme from a butanol-producing Clostridium strain.
Rhodopseudomonas palustris, which is a member of the purple non-sulfur bacteria group, is often noted for its metabolic versatility. One of these metabolic activities, which is H2 production by photofermentation, utilizes organic carbon under anaerobic conditions with light and has become an attractive application for bioenergy production . The ability to produce H2 at high concentrations is driven by an excess of reducing equivalents that accumulate during photoheterotrophic growth . Here, we aim to exploit the excess reducing equivalents generated during photoheterotrophic growth to drive an engineered reduction reaction for the obligate production of the liquid fuel n-butanol (hereafter butanol).
When growing on acetate, McKinlay and Harwood  reported that approximately half of all reducing equivalents that were generated were used for biosynthesis, while the other half required re-oxidation by some alternate pathway such as CO2 fixation, H2 evolution, or a combination of both. These two pathways were found to be complementary in this role, and when one route was made unavailable, the other facilitated the re-oxidation of all reducing equivalents with minimal changes occurring elsewhere in the central metabolism. However, when fed a more reduced substrate, such as n-butyrate (hereafter butyrate), it was observed that under non-H2 evolving conditions, R. palustris could only grow if exogenous CO2 was supplemented into the medium [2, 3]. This finding demonstrates that under these conditions, R. palustris generates more reducing equivalents than available electron sinks can accommodate, and ultimately suffers a lethal redox imbalance. Numerous studies have used these principles to improve H2 yields from photofermentation using R. palustris [4,5,6].
Recently, Fixen et al.  engineered a strain of R. palustris to utilize excess reducing equivalents to reduce CO2 into CH4 by remodeling the nitrogenase enzyme active site and expressing the modified gene from the chromosome. A NifA* strain of R. palustris was chosen as the host because it expresses the nitrogenase gene constitutively even in the presence of NH3. The result was a light-dependent R. palustris strain that produced CH4 with ATP coming from cyclic photophosphorylation and electrons coming from either an organic substrate or thiosulfate. This reduction is, however, not obligatory and a large excess of H2 was produced by the nitrogenase compared to CH4 . Another attractive route to redirect the excess reducing equivalents and to maintain redox balance when growing on butyrate is to produce the fuel butanol through an obligate reduction reaction. The advantage of the obligate reduction is the lack of side products to achieve a high selectivity (i.e., products compared to the consumed substrate).
Butanol is an attractive biofuel molecule due to its higher energy density, lower volatility, and complete intersolubility with traditional fuels compared to ethanol . Butanol is traditionally produced via fermentation of monosaccharides by certain strains of Clostridia. Clostridium acetobutylicum ATCC 824, which has been widely studied for this activity, is generally seen as the model organism for solvent production and previously has been used as a source of genes for engineered butanol production [9,10,11]. In C. acetobutylicum, reduction of butyrate proceeds through a two-phase pathway. First, butyrate is activated to butyryl-CoA, which can occur by two different routes: (1) a two-step butyrate kinase/phosphotransbutyrylase reaction (Fig. 1); or (2) a one-step reaction with butyrate/acetoacetyl-CoA transferase transferring a CoA group from acetoacetyl-CoA to butyrate. Second, butyryl-CoA is reduced to butanol through the butyraldehyde intermediate [12, 13]. The AdhE2 enzyme from C. acetobutylicum ATCC 824 (hereafter AdhE2824 and adhE2 824 for its gene) reduces butyryl-CoA to butyraldehyde and subsequently butanol (Fig. 1) . Analysis of AdhE2824 has revealed that the protein is a bifunctional NADH-dependent fusion protein with the N-terminal domain retaining highly conserved sequences of the aldehyde dehydrogenase family and the C-terminal domain retaining conserved sequences of the iron-containing alcohol dehydrogenase family . A second alcohol dehydrogenase enzyme BdhB has also been found to aid in the reduction of butyraldehyde to butanol in C. acetobutylicum, though its activity is not necessary for butanol production .
In planning to introduce the butanol production pathway together with the existing catabolism of butyrate by R. palustris, we established that before butyrate can be converted to two molecules of acetyl-CoA via anaerobic β-oxidation, butyryl-CoA is first formed through native butyrate kinase/phosphotransbutyrylase enzymes (Fig. 1). Under photoheterotrophic redox imbalance, the first oxidation reaction for the metabolism of butyrate (i.e., butyryl-CoA to crotonyl-CoA in Fig. 1) would likely be inhibited by the thermodynamic limitations of the reduced cofactor pool. Because of this, we anticipated the butyryl-CoA intermediate to accumulate and provide an ideal target for shunting excess reducing equivalents for the production of butanol by introducing only adhE2 824 (the bifunctional aldehyde/alcohol dehydrogenase gene) into R. palustris (Fig. 1). After necessary optimizations, the production of butanol indicated that this assumption was correct.
A direct introduction of plasmid-based adhE2 824 did not result in the production of AdhE2824 (and butanol production) due to the major incompatibilities in codon utilization preference between the C. acetobutylicum source and the R. palustris host. Therefore, we pursued two alternate attempts to introduce the activity of AdhE2 by using: (1) a gene homologous to adhE2 824 from R. palustris strain BisB18 (hereafter adhE BisB18 and AdhEBisB18 for its protein); and (2) a codon-optimized adhE2 824 (hereafter adhE2 opti and AdhE2opti for its protein). In all our growth experiments, we omitted HCO3 − (CO2) from, and added NH3 to, the growth medium to block the normal routes of shunting reducing equivalents, resulting in an obligate butyrate reduction. This strong selection pressure will only allow growth when butanol production occurs. In our hands, butanol production through the homologous AdhEBisB18 did not work, but AdhE2opti enabled recovered growth through butanol production. Subsequently, we were able to increase the butanol selectivity and production rate in further experiments, accomplishing our aim of generating an engineered reduction reaction for the obligate production of the liquid fuel butanol in R. palustris.
Media and culture conditions
Rhodopseudomonas palustris CGA009 was cultured in anaerobic butyrate Rhodospirillaceae medium (RM) , consisting of 1.64 g NaHCO3, 0.5 g KH2PO4, 0.5 g K2HPO4, 0.2 g MgSO4·7 H2O, 0.4 g NaCl, 0.05 g CaCl2·2 H2O, 0.01 g Fe-Citrate, 1.0 g (NH4)2SO4, 1.12 mL 10 M NaOH, and 1 mL trace metals per liter. The trace-metal solution consisted of 0.25 mL concentrated HCl, 70.0 mg ZnCl2, 100.0 mg MnCl2·4 H2O, 60.0 mg H3BO3, 200.0 mg CoCl2·6 H2O, 20.0 mg CuCl2·2 H2O, 20.0 mg NiCl2·6 H2O, and 40.0 mg Na2MoO4·2 H2O L−1. After autoclaving, 1.0 g sterile butyric acid was added as the substrate. No toxicity or inhibition of R. palustris growth was detected at this concentration (Additional file 1). 1 mL trace vitamins was added to the medium. The trace-vitamin solution consisted of 1.0 g p-aminobenzoic acid, 1.0 g thiamine, and 0.1 g biotin L−1. Sterile medium (50 mL) with a final pH of 7.0 was added to autoclaved serum bottles and sparged for 20 min with sterile-filtered 80:20 N2/CO2 gas. For the growth of engineered strains under restrictive butyrate-to-butanol conversion conditions, NaHCO3 was omitted from the medium, the addition of 10 M NaOH was increased to 2.3 mL total, and serum bottles were sparged for 20 min with N2 instead of N2/CO2. In addition, 200 μg mL−1 of kanamycin sulfate was added for retention of the plasmid in all growth conditions with engineered strains. All serum bottle cultures were grown in an environmental growth chamber (GC8-2VH, EGC, Chagrin Falls, OH) at 30 °C with 80 μmol of photosynthetically active photons s−1 m−2 (photons between 400 and 700 nm) illumination from both fluorescent and incandescent lamps.
DNA manipulation and cloning
We used the green fluorescence protein GFPmut2  as a reporter protein to screen for the activity of different expression vectors in R. palustris under aerobic conditions. We tested two different plasmids as candidate expression vectors: (1) the endogenously derived pMG105 ; and (2) the broad-host pBBR1MCS-2 . For detection by immunoreaction, we introduced a C-terminal 6X His tag prior to the termination site with PCR by the reverse primer (Table 1). GFPmut2 was inserted into the pMG105 and pBBR1MCS-2 vectors under the control of a range of promoters, including the endogenous phosphoenolpyruvate carboxykinase (ppckA) from R. palustris, arabinose (ara), and lactose (lac) promoters. The adhE BisB18 and adhE2 824 genes were amplified from R. palustris strain BisB18 and C. acetobutylicum ATCC 824, respectively, using Phusion HF polymerase (New England Biolabs, Ipswich, MA) and directionally ligated into the pBBR1MCS-2 expression vector using the XhoI/XbaI cloning sites. Restriction and ligation enzymes were purchased from Promega (Fitchburg, WI) and all constructs were verified with Sanger sequencing (Cornell University Genomics Facility, Ithaca, NY). Escherichia coli DH5α was used for routine cloning of constructs. Both E. coli and R. palustris were transformed by electroporation with a Bio-Rad Gene Pulser II Electroporator (Bio-Rad, Hercules, CA) with 1.8 kV, 200 Ω, 25 µF for E. coli, and 2.5 kV, 400 Ω, 25 µF for R. palustris. Electroporated cells were recovered in SOC for 1 and 3 h before plating on lysogeny broth (LB) agar plates with 50 and 200 μg mL−1 kanamycin sulfate for isolation of individual E. coli and R. palustris transformants, respectively.
To successfully produce AdhE2824 in R. palustris, codon optimization was pursued. Optimization design was performed in-house and largely utilized the one amino acid one codon approach, except in cases of multiple codon repeats and unfavorable energetics of mRNA secondary structure (Additional file 2). adhE2 opti with a C-terminal 6X His tag was synthesized by GenScript (Piscataway, NJ) with flanking XhoI/XbaI restriction sites that were used for directional cloning into pBBR1MCS-2.
SDS-PAGE and western blot detection
For both SDS-PAGE and subsequent western blot detection, 50 mL of late-log E. coli and R. palustris cells were harvested by centrifugation (10 min, 12,000×g, 4 °C) and resuspended in 1 mL 1X Laemmli sample buffer  with a protease inhibitor (Protease inhibitor cocktail, Promega, Fitchburg, WI). The suspension was lysed by ultrasonication on ice using three 30 s on 1 min off pulses of 20 W (Branson Sonifier 150, Emerson, Danbury, CT). After sonication, debris was removed by centrifugation and the total soluble protein concentration was determined by BCA assay (Pierce BCA protein assay kit, Thermo Scientific, Rockford, IL). Forty-five μL of 3 mg mL−1 soluble protein samples and 10 μL of Precision Plus Protein™ standard (Bio-Rad, Hercules, CA) were loaded into pockets of a Mini-PROTEAN® TGX gel (Bio-Rad, Hercules, CA) and run at 100 V for 1.5 h before blotting onto PVDF membranes at 100 V for 1 h. For the initial western blot, immunoreactive bands were developed using manufacturer-suggested protocols for the 1 h Western™ Standard Kit for rabbit (GenScript, Piscataway, NJ) and a 6X His primary antibody (Rabbit polyclonal, GeneTex, Irvine, CA). We used a similar protocol a second time to improve the sensitivity with these changes: we used 0.1 mL OD600 1.0 cells with 0.1 mL Laemmli sample buffer. The samples were boiled for 10 min before centrifugation. 15 µL of each sample was loaded into separate gel lanes. For development, an enhanced chemiluminescence protocol was used using the manufacturer-suggested protocol for promega ECL western blotting substrate (Fitchburg, WI). The film was allowed to incubate with the PVDF membrane for 30 min prior to development.
Bioreactor configuration and operation
Prior to operation, sealed glass bioreactors were autoclaved and filled with 350 mL sterile RM with 1 g L−1 butyrate and without HCO3 −. The bioreactors were equipped with a recirculating water jacket that was operated at 30 °C. Magnetic stir plates (IKA Works Inc., Wilmington, NC) provided stirring of the reactors at 450 rpm. The bioreactors were continuously sparged with sterile-filtered N2 gas at a flow rate of 0.1 L min−1 to maintain an anaerobic environment. After sparging for 1 h to ensure anaerobic reactor conditions, 1 mL of actively growing R. palustris cells at an OD600 of 0.1 was inoculated into the bioreactor. The preculture of R. palustris pBBR1MCS-2 lacp adhE2 opti was grown in anaerobic RM without HCO3 −, while the preculture of R. palustris pBBR1MCS-2 lacp GFPmut2 (control) was grown in aerobic RM without HCO3 −. Gas entered the liquid phase of each bioreactor through a sparging stone to increase the gas/liquid interface for butanol stripping. Off-gas was connected to 2 Friedrich’s condensers in parallel with recirculating water at 5 °C to condense the butanol that had been removed by gas stripping. Gas was released into the atmosphere after one passing through the condensers. The bioreactor volume was illuminated externally with a single 60-W incandescent lamp as used previously . The bioreactors were operated in batch-mode and liquid samples were periodically taken for OD600 and metabolite measurements.
Butyrate, butanol, and acetate were measured via HPLC (600 HPLC, Waters, Milford, MA) with a refractive index detector and an Aminex HPX-87H column (Bio-Rad, Hercules, CA). The column temperature was set to 60 °C, and a 5-mM sulfuric acid eluent at a flow rate of 0.6 mL min−1 was used as the mobile phase. Butanol and butyraldehyde were measured with a gas chromatograph (HP 5890, Hewlett Packard, Palo Alto, CA) with a 7673 autoinjector and flame ionization detector. The gas chromatograph contained a custom-made packed bed glass column of 1.8 m × 2 mm id (Supelco, Sigma-Aldrich, St. Louis, MO). The support matrix of the column was Chromosorb W/AW80 over 100 mesh; phases were preconditioned: phase A was 10% Carbowax-20M; phase B was 0.1% phosphoric acid. The inlet and detector temperatures were 220 and 240 °C, respectively. The column temperature profile was 100 °C for 2 min, a ramp of 40 °C min−1 to 180 °C with a 5-min hold. Butanol concentrations were verified by both HPLC and GC platforms.
Results and discussion
GFP screen of expression vectors in R. palustris CGA009
Different combinations of promoter and vector sequences were screened for highest levels of expression within R. palustris. The endogenous ppckA promoter of R. palustris is one of the few promoter sequences that had been previously validated to be active in R. palustris under gluconeogenic conditions . However, we only observed a minimal increase in fluorescence for both the pMG105 and pBBR1MCS-2 ppckAp GFPmut2 constructs compared to the empty vector controls (Fig. 2). While the pBBR1MCS-2 arap GFPmut2 construct yielded slightly higher levels compared to ppckAp, the pBBR1MCS-2 lacp GFPmut2 construct offered the best expression overall and was selected as the vector for all further expressions in R. palustris (Fig. 2). Expression from the lac promoter in R. palustris was not affected by the addition of the inducer IPTG and catabolite repression only showed a minor decrease in fluorescence (Additional file 3). This suggests that the lac regulatory response system utilized by E. coli is not functional in R. palustris. Although the pBBR1MCS-2 lac promoter system provided the best expression in R. palustris of all the combinations tested (Fig. 2), the fluorescent signal was still ~13 times lower than the signal from E. coli pBBR1MCS lacp GFPmut2 (data not shown).
Engineering adhE BisB18 in R. palustris
As determined by the lack of growth under restrictive conditions, the absence of butanol production (data not shown), and the lack of a signal in a western blot under permissive conditions (Fig. 3a, AdhEBisB18 Aerobic), direct introduction of adhE2 824 in the pBBR1MCS-2 lacp construct did not result in the production of AdhE2824 in R. palustris. This is likely due to the large GC discrepancy (30% vs. 65%) and codon utilization preference between C. acetobutylicum and R. palustris . To address this, we looked for homologous proteins within R. palustris. Indeed, endogenous alcohol/aldehyde dehydrogenases with homology to AdhE2824 have been found to function in the production of butanol from butyrate when overexpressed. This was previously demonstrated with the fucO gene in E. coli, where overexpression allowed production of butanol from an engineered reverse β-oxidation pathway where no butanol production was observed previously . In an attempt to enable butanol production in R. palustris, a BLAST search for an endogenous enzyme yielded the identification of RPC_4481 (AdhEBisB18) in R. palustris strain BisB18, which is an aldehyde/alcohol dehydrogenase with high homology to both domains of AdhE2824.
Before expressing adhE BisB18 as a candidate butanol-pathway gene, R. palustris strain BisB18 was screened for its ability to grow by producing butanol from butyrate under restrictive electron sink conditions (without HCO3 −). Although no growth or butanol production occurred (data not shown), it is possible that wild-type production of this protein under these conditions is either inactive or insufficient for an observable effect, similar to the FucO protein in E. coli . Because no experimental studies have previously characterized the activity of AdhEBisB18, its highly conserved primary sequence was further investigated for its structural homology to the verified butanol-producing AdhE2824.
Using Swiss-Pdbviewer/DeepView , three-dimensional structures of the individual protein domains for both AdhE2824 and AdhEBisB18 were generated. The predicted tertiary structure of both aldehyde/alcohol dehydrogenases appeared highly conserved with regard to overall shape, location, active site residues, channel volume, and surface charge (Additional file 4). This suggested that the two proteins likely recognized similar substrates. We, therefore, expressed adhE BisB18 in R. palustris using the pBBR1MCS-2 lacp construct. Transformants were inoculated into RM with butyrate and with and without HCO3 − to screen for viability and butanol production, respectively. Although transformants grew up with HCO3 −, no growth was observed without HCO3 − condition (data not shown), suggesting that either AdhEBisB18 was not successfully produced by R. palustris CGA009 or did not have butyryl-CoA-to-butanol conversion activity.
To investigate this further, a crude cell extract of the culture with HCO3 −, along with E. coli harboring the same construct, was analyzed via western blotting. Although the strong signal at the expected size of 94.5 kD for AdhEBisB18 for E. coli pBBR1MCS-2 lacp adhE BisB18 demonstrated that the protein was produced (Fig. 3a, AdhEBisB18), the lack of a signal from R. palustris pBBR1MCS-2 lacp adhE BisB18 suggested that the protein was either not produced or was produced at levels below detection (Fig. 3a, AdhEBisB18 Aerobic). Aerobic conditions were preferred since anaerobic cultures often had photosynthetic pigments overlapping the expected size for AdhE (Fig. 3a, AdhE2opti Anaerobic). The lack of growth and butanol production under restrictive conditions, however, demonstrated conclusively that R. palustris pBBR1MCS-2 lacp adhE BisB18 was not endowed with butyrate-to-butanol conversion functionality (data not shown).
Engineering adhE2 opti in R. palustris
Another route to circumvent the large GC discrepancy and codon utilization preference between C. acetobutylicum and R. palustris was to perform codon optimization on adhE2 824 to enable expression. The codon utilization table for R. palustris CGA009 was determined by analysis of all genomic coding sequences and was retrieved from the Kazusa database . The adhE2 opti gene resulted in an increase of the average codon utilization frequency in R. palustris from 23% for the original adhE2 824 to 64% for the optimized adhE2 opti (Additional file 5). Codon optimization resulted in a shift of the GC content from 32.5% for adhE2 824 to 62% for adhE2 opti. In addition, the presence of rare codons (<10% utilization frequency), which previously composed 47.4% of all codons in adhE2 824, was entirely eliminated. One disadvantage to this approach of codon optimization is that rare codons are known to induce ribosomal pausing, which has been demonstrated important for proper folding of proteins. Thus, with the elimination of all rare codons from adhE2 opti, the elimination of codon-dependent ribosomal pausing could lead to improper folding and aggregation of the target protein . For comparison with adhE2 opti, the codon utilization frequency of adhE BisB18 averaged 55% (Additional file 5), the gene possessed an overall 63% GC content, and the frequency of rare codons within the gene was 5%, demonstrating the utilization of rare codons for proper translation of proteins in endogenous alcohol/aldehyde genes (Additional file 6).
To test whether adhE2 opti would be expressed in R. palustris CGA009 and enable growth by conversion of butyrate to butanol, the gene was introduced using the pBBR1MCS-2 lacp construct. Transformants were then inoculated into RM with butyrate and with and without HCO3 − for growth screening and western blot analysis. This engineered strain was able to grow without HCO3 − (restrictive conditions) and it produced a measureable amount of butanol, demonstrating the successful production of AdhE2opti with a functional activity. While a detectable signal was observed on a western blot for the E. coli cloning culture, no signal was initially detected from R. palustris pBBR1MCS-2 lacp adhE2 opti (Fig. 3a, AdhE2opti Anaerobic and AdhE2opti Aerobic). We performed a second western blot using a more sensitive X-ray film imaging method. For aerobically grown R. palustris pBBR1MCS-2 lacp adhE2 opti, we observed a single signal at approximately 94 kD, demonstrating successful expression (Fig. 3b, AdhE2opti Aerobic). As stated above, GFPmut2 production from the pBBR1MCS lac promoter is much lower in R. palustris than in E. coli, so the lack of a signal in the majority of the samples is likely due to low protein production in R. palustris.
Butanol production rates and selectivity for R. palustris pBBR1MCS-2 lacp adhE2 opti
Because of the rescued growth and butanol production observed, R. palustris pBBR1MCS-2 lacp adhE2 opti was further investigated for its ability to grow under both permissive and restrictive conditions and these results were compared to R. palustris pBBR1MCS-2 lacp GFPmut2 as a control. R. palustris pBBR1MCS-2 lacp adhE2 opti was capable of growth without HCO3 − (restrictive conditions) (Fig. 4b), while this growth did not occur for the control strain. The growth for R. palustris pBBR1MCS-2 lacp adhE2 opti occurred on a long time scale and reached a relatively low maximum cell density of ~0.06 OD600 (Fig. 4b) at neutral pH levels (Fig. 4d). Following the onset of the stationary phase, consumption of butyrate and production of butanol still occurred (Fig. 4f, h), but the cultures only reached a maximum butanol concentration of <0.4 mM after extended incubation (Fig. 4h).
Although only low levels of butyrate were consumed (Fig. 4f), butanol production occurred with a 39.7 ± 7.6% selectivity based on moles of carbon. The theoretical maximum for butanol production is based on all electrons derived from complete oxidation of one molecule of butyrate reducing five molecules of butyrate to butanol (6 butyrate → 5 butanol + 4 CO2) for 83% conversion efficiency. However, based on previous empirical measurements of excess reducing equivalents of R. palustris growing with butyrate by McKinlay et al. , only an estimated 0.45 mol of butanol should be produced per mole of butyrate consumed for an effective 45% selectivity (6 butyrate → 2.7 butanol + biomass). Therefore, almost all reducing equivalents predicted to be in excess were being converted into butanol by R. palustris pBBR1MCS-2 lacp adhE2 opti (40% vs. 45%), which we had anticipated for an obligatory reaction.
To further investigate why R. palustris pBBR1MCS-2 lacp adhE2 opti showed low growth rates, R. palustris was screened for product (butanol) toxicity. However, no inhibition of growth was observed in R. palustris at butanol concentrations up to ~30 mM (Additional file 7), removing our concerns about a concentration of ~0.4 mM (Fig. 4h). In addition, we operated a bioreactor with product removal to investigate whether this would increase growth and butanol production rates by removing butanol from the substrate pool, thereby, improving the thermodynamics for butanol formation and eliminating any end-product inhibition. Sparging of N2 has been previously used to remove butanol from fermentation broths . Here, N2 gas was continuously sparged into 350-mL anaerobic bioreactors to remove butanol produced by R. palustris pBBR1MCS-2 lacp adhE2 opti in batch-mode for 6 consecutive weeks. The R. palustris pBBR1MCS-2 lacp adhE2 opti culture grew for that entire period up to an OD600 of 0.12 (Additional file 8: Fig. S6A), and achieved an average maximum butanol production rate of 0.0017 g L−1 day−1 (Additional file 8: Fig. S6D) with a maximum butanol concentration of 0.35 mM (Additional file 8: Fig. S6C) without depleting butyrate (Additional file 8: Fig. S6B). The lack of a considerable improvement compared to the earlier experiment showed that neither butanol toxicity nor thermodynamics nor end-product inhibition was a main reason for low growth rates. Finally, we found that the extraction system removed a majority of the butanol without removing any of the butyrate (Additional file 8: Fig. S6E, F).
The intermediate butyraldehyde within microbial cells is another metabolite that could have inhibited R. palustris, even though the maximum concentration in the supernatant was measured to be only 50 μM. However, due to the fusion nature of AdhE2opti, we anticipate that butyraldehyde would rapidly be converted into butanol with a high flux rate, resulting in a low effective concentration . Without butanol and butyraldehyde toxicity, we can only speculate about why growth and butanol production rates of our engineered strain are so low. Since AdhE2opti is provided as the sole route to maintaining redox balance, all excess electrons must be funneled through this pathway. When the capacity of the AdhE2opti pool is outweighed by the demand for regenerating oxidized reducing equivalents, this would introduce a bottleneck in the maximum rate of metabolism and growth. However, more work is required to be conclusive.
Some butanol production did also occur for R. palustris pBBR1MCS-2 lacp adhE2 opti with HCO3 − (permissive conditions) when exogenous CO2 was available for maintaining redox balance (Fig. 4a, g) at a neutral pH level (Fig. 4c). Because the preculture had originally been maintained without HCO3 −, a long lag phase and carryover of butanol producing activity was observed for R. palustris pBBR1MCS-2 lacp adhE2 opti (Fig. 4a), resulting in a maximum concentration of 0.37 mM with 11.5 ± 2.4% selectivity (Fig. 4e, g). Though the transient production of butanol was a surprise, this lower selectivity for the permissive condition was expected because more reducing equivalents could be used in central metabolism following the introduction of exogenous HCO3 −, resulting in much higher growth (~10× OD at 200 h), the complete depletion of butyrate (Fig. 4e), no production of acetate, and ultimate consumption of all butanol produced (Fig. 4g). This contrasts with the restrictive condition where all growth proceeds through the obligate engineered pathway, butyrate is not depleted, and all butanol produced remains in solution.
Optimizing butanol production by pregrowing R. palustris pBBR1MCS-2 lacp adhE2 opti
In an attempt to improve the maximum butanol concentrations produced by the engineered strain, 50 mL aerobic cultures was grown to mid-log phase, washed 3 times with sterile RM, resuspended in 0.5 mL RM, and used as an inoculum into 20-mL anaerobic cultures under restrictive conditions. This procedure ensured that the butyrate-to-butanol conversion was not catalyst limited as the initial OD600 of these cultures was >1. By using a concentrated inoculum, R. palustris pBBR1MCS-2 lacp adhE2 opti produced butanol concentrations of greater than 1.5 mM at volumetric butanol production rates of 0.034 g L−1 day−1 (Fig. 5d, f). All butanol was produced within 100 h of inoculation (Fig. 5d, f). Following 100 h, however, the R. palustris pBBR1MCS-2 lacp adhE2 opti culture continued to grow and consume butyrate (Fig. 5a, c) with the main carbon product switching from butanol to acetate (Fig. 5b). Acetate had not been detected before in any of the previous anaerobic experiments and in terms of maintaining redox balance, acetate is an unfavorable redox product compared to butanol. This metabolism reduced the butanol selectivity (Fig. 5e). We do not know exactly why this switch occurred, but it was likely a result from having a high biomass culture that underwent a rapid switch from aerobic to anaerobic conditions.
The introduction of adhE2 opti into R. palustris via a plasmid-based construct is an example of metabolically engineered fuel production that proceeds as an obligate strategy for growth. Although growth and production rates were low for this system, our selectivity of 39.7 ± 7.6% matched well with the empirical expectation of 45.0% under these conditions. Because the low growth rates and butanol production are likely constrained by redox flux through the engineered pathway, rational design and directed evolution could be used to further improve the performance of this activity. If improved, the system could be easily integrated into a simple bioreactor design with highly efficient selectivity and recovery. While this proof-of-concept for harnessing an obligate reducing metabolism to drive engineered reactions used butanol as a target molecule, this redox-based driving force could be used to achieve other reactions for the production of desirable chemicals that require reduction reactions.
McKinlay JB, Harwood CS. Carbon dioxide fixation as a central redox cofactor recycling mechanism in bacteria. Proc Natl Acad Sci USA. 2010;107(26):11669–75.
Hillmer P, Gest H. H2 metabolism in the photosynthetic bacterium Rhodopseudomonas capsulata: H2 production by growing cultures. J Bacteriol. 1977;129(2):724–31.
McKinlay JB, Harwood CS. Calvin cycle flux, pathway constraints, and substrate oxidation state together determine the H2 biofuel yield in photoheterotrophic bacteria. mBio. 2011;2(2):e00323.
Barbosa MJ, Rocha JMS, Tramper J, Wijffels RH. Acetate as a carbon source for hydrogen production by photosynthetic bacteria. J Biotechnol. 2001;85(1):25–33.
Gosse JL, Engel BJ, Rey FE, Harwood CS, Scriven LE, Flickinger MC. Hydrogen production by photoreactive nanoporous latex coatings of nongrowing Rhodopseudomonas palustris CGA009. Biotechnol Prog. 2007;23(1):124–30.
Chen CY, Saratale GD, Lee CM, Chen PC, Chang JS. Phototrophic hydrogen production in photobioreactors coupled with solar-energy-excited optical fibers. Int J Hydrog Energ. 2008;33(23):6886–95.
Fixen KR, Zheng Y, Harris DF, Shaw S, Yang Z-Y, Dean DR, Seefeldt LC, Harwood CS. Light-driven carbon dioxide reduction to methane by nitrogenase in a photosynthetic bacterium. Proc Natl Acad Sci USA. 2016;113(36):10163–7.
Jin C, Yao M, Liu H, Lee C-F, Ji J. Progress in the production and application of n-butanol as a biofuel. Renew Sustain Energ Rev. 2011;15(8):4080–106.
Nolling J, Breton G, Omelchenko MV, Makarova KS, Zeng QD, Gibson R, Lee HM, Dubois J, Qiu DY, Hitti J, et al. Genome sequence and comparative analysis of the solvent-producing bacterium Clostridium acetobutylicum. J Bacteriol. 2001;183(16):4823–38.
Inui M, Suda M, Kimura S, Yasuda K, Suzuki H, Toda H, Yamamoto S, Okino S, Suzuki N, Yukawa H. Expression of Clostridium acetobutylicum butanol synthetic genes in Escherichia coli. Appl Microbiol Biotechnol. 2008;77(6):1305–16.
Atsumi S, Cann AF, Connor MR, Shen CR, Smith KM, Brynildsen MP, Chou KJ, Hanai T, Liao JC. Metabolic engineering of Escherichia coli for 1-butanol production. Metab Eng. 2008;10(6):305–11.
Desai RP, Harris LM, Welker NE, Papoutsakis ET. Metabolic flux analysis elucidates the importance of the acid-formation pathways in regulating solvent production by Clostridium acetobutylicum. Metab Eng. 1999;1(3):206–13.
Sillers R, Al-Hinai MA, Papoutsakis ET. Aldehyde-alcohol dehydrogenase and/or thiolase overexpression coupled with CoA transferase downregulation lead to higher alcohol titers and selectivity in Clostridium acetobutylicum fermentations. Biotechnol Bioeng. 2009;102(1):38–49.
Dürre P, Fischer RJ, Kuhn A, Lorenz K, Schreiber W, Sturzenhofecker B, Ullmann S, Winzer K, Sauer U. Solventogenic enzymes of Clostridium acetobutylicum: catalytic properties, genetic organization, and transcriptional regulation. FEMS Microbiol Rev. 1995;17(3):251–62.
Nair RV, Bennett GN, Papoutsakis ET. Molecular characterization of an aldehyde/alcohol dehydrogenase gene from Clostridium acetobutylicum ATCC 824. J Bacteriol. 1994;176(3):871–85.
Fontaine L, Meynial-Salles I, Girbal L, Yang X, Croux C, Soucaille P. Molecular characterization and transcriptional analysis of adhE2, the gene encoding the NADH-dependent aldehyde/alcohol dehydrogenase responsible for butanol production in alcohologenic cultures of Clostridium acetobutylicum ATCC 824. J Bacteriol. 2002;184(3):821–30.
Lee CM, Chen PC, Wang CC, Tung YC. Photohydrogen production using purple nonsulfur bacteria with hydrogen fermentation reactor effluent. Int J Hydrog Energy. 2002;27(11–12):1309–13.
Cormack BP, Valdivia RH, Falkow S. FACS-optimized mutants of the green fluorescent protein (GFP). Gene. 1996;173(1):33–8.
Inui M, Roh JH, Zahn K, Yukawa H. Sequence analysis of the cryptic plasmid pMG101 from Rhodopseudomonas palustris and construction of stable cloning vectors. Appl Environ Microbiol. 2000;66(1):54–63.
Kovach ME, Elzer PH, Hill DS, Robertson GT, Farris MA, Roop RM, Peterson KM. Four new derivatives of the broad-host-range cloning vector pBBR1MCS, carrying different antibiotic-resistance cassettes. Gene. 1995;166(1):175–6.
Protocols CSH. Laemmli sample buffer (2X). Cold Spring Harbor Protocols. 2006;2006(1):pdb.rec10424.
Doud DFR, Angenent LT. Toward electrosynthesis with uncoupled extracellular electron uptake and metabolic growth: enhancing current uptake with Rhodopseudomonas palustris. Environ Sci Technol Lett. 2014;1(9):351–5.
Inui M, Nakata K, Roh JH, Zahn K, Yukawa H. Molecular and functional characterization of the Rhodopseudomonas palustris no. 7 phosphoenolpyruvate carboxykinase gene. J Bacteriol. 1999;181(9):2689–96.
Gustafsson C, Govindarajan S, Minshull J. Codon bias and heterologous protein expression. Trends Biotechnol. 2004;22(7):346–53.
Dellomonaco C, Clomburg JM, Miller EN, Gonzalez R. Engineered reversal of the beta-oxidation cycle for the synthesis of fuels and chemicals. Nature. 2011;476(7360):355–9.
Guex N, Peitsch MC. SWISS-MODEL and the Swiss-Pdb Viewer: an environment for comparative protein modeling. Electrophoresis. 1997;18(15):2714–23.
Nakamura Y, Gojobori T, Ikemura T. Codon usage tabulated from international DNA sequence databases: status for the year 2000. Nucleic Acids Res. 2000;28(1):292.
Rosano GL, Ceccarelli EA. Rare codon content affects the solubility of recombinant proteins in a codon bias-adjusted Escherichia coli strain. Microb Cell Fact. 2009;8:41.
Qureshi N, Blaschek HP. Recovery of butanol from fermentation broth by gas stripping. Renew Energy. 2001;22(4):557–64.
DeLisa MP, Conrado RJ. Synthetic metabolic pipelines. Nat Biotechnol. 2009;27(8):728–9.
DFRD designed and executed experiments, performed analysis, interpreted data, and drafted the manuscript. ECH executed experiments, performed analysis, and contributed to data interpretation. HR and BM advised on how to perform the experiments. GJ and LTA conceived the study and participated in its design. LTA coordinated the study and contributed to the drafting of the manuscript. All authors read and approved the final manuscript.
We would like to thank Dr. Meena Haribal, Dr. Lubna Richter, Prof. Dr. Matt DeLisa, Dr. Jason Boock, and Jen Schmidt for their expertise and advice (all at Cornell University). We would like to thank Prof. Dr. Caroline Harwood (University of Washington) for providing the R. palustris strain BisB18 and Wendy Hoose (Cornell University) for providing technical assistance.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.
This project and its researchers were funded in part through an Academic Venture Fund award from the Atkinson Center for a Sustainable Future to LTA and GJ, and the Advanced Research Projects Agency Energy Open 2012 program Project Number DE-AR0000312 to LTA. We acknowledge support by Deutsche Forschungsgemeinschaft and Open Access Publishing Fund of University of Tübingen.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.