Characterizing Strain Variation in Engineered E. coli Using a Multi-Omics-Based Workflow

Understanding the complex interactions that occur between heterologous and native biochemical pathways represents a major challenge in metabolic engineering and synthetic biology. We present a workflow that integrates metabolomics, proteomics, and genome-scale models of Escherichia coli metabolism to study the effects of introducing a heterologous pathway into a microbial host. This workflow incorporates complementary approaches from computational systems biology, metabolic engineering, and synthetic biology; provides molecular insight into how the host organism microenvironment changes due to pathway engineering; and demonstrates how biological mechanisms underlying strain variation can be exploited as an engineering strategy to increase product yield. As a proof of concept, we present the analysis of eight engineered strains producing three biofuels: isopentenol, limonene, and bisabolene. Application of this workflow identified the roles of candidate genes, pathways, and biochemical reactions in observed experimental phenomena and facilitated the construction of a mutant strain with improved productivity. The contributed workflow is available as an open-source tool in the form of iPython notebooks.

Impact of elevated nitrate on sulfate-reducing bacteria: a comparative study of Desulfovibrio vulgaris

Sulfate-reducing bacteria have been extensively studied for their potential in heavy-metal bioremediation. However, the occurrence of elevated nitrate in contaminated environments has been shown to inhibit sulfate reduction activity. Although the inhibition has been suggested to result from the competition with nitrate-reducing bacteria, the possibility of direct inhibition of sulfate reducers by elevated nitrate needs to be explored. Using Desulfovibrio vulgaris as a model sulfate-reducing bacterium, functional genomics analysis reveals that osmotic stress contributed to growth inhibition by nitrate as shown by the upregulation of the glycine/betaine transporter genes and the relief of nitrate inhibition by osmoprotectants. The observation that significant growth inhibition was effected by 70 mM NaNO(3) but not by 70 mM NaCl suggests the presence of inhibitory mechanisms in addition to osmotic stress. The differential expression of genes characteristic of nitrite stress responses, such as the hybrid cluster protein gene, under nitrate stress condition further indicates that nitrate stress response by D. vulgaris was linked to components of both osmotic and nitrite stress responses. The involvement of the oxidative stress response pathway, however, might be the result of a more general stress response. Given the low similarities between the response profiles to nitrate and other stresses, less-defined stress response pathways could also be important in nitrate stress, which might involve the shift in energy metabolism. The involvement of nitrite stress response upon exposure to nitrate may provide detoxification mechanisms for nitrite, which is inhibitory to sulfate-reducing bacteria, produced by microbial nitrate reduction as a metabolic intermediate and may enhance the survival of sulfate-reducing bacteria in environments with elevated nitrate level.

Redesigning enzymes based on adaptive evolution for optimal function in synthetic metabolic pathways

Nature has balanced most metabolic pathways such that no one enzyme in the pathway controls the flux through that pathway. However, unnatural or nonnative, constructed metabolic pathways may have limited product flux due to unfavorable in vivo properties of one or more enzymes in the pathway. One such example is the mevalonate-based isoprenoid biosynthetic pathway that we previously reconstructed in Escherichia coli. We have used a probable mechanism of adaptive evolution to engineer the in vivo properties of two enzymes (3-hydroxy-3-methylglutaryl-CoA reductase [tHMGR] and many terpene synthases) in this pathway and thereby eliminate or minimize the bottleneck created by these inefficient or nonfunctional enzymes. Here, we demonstrate how we significantly improved the productivity (by approximately 1000 fold) of this reconstructed biosynthetic pathway using this strategy. We anticipate that this strategy will find broad applicability in the functional construction (or reconstruction) of biological pathways in heterologous hosts.

ClusterCAD: a computational platform for type I modular polyketide synthase design

ClusterCAD is a web-based toolkit designed to leverage the collinear structure and deterministic logic of type I modular polyketide synthases (PKSs) for synthetic biology applications. The unique organization of these megasynthases, combined with the diversity of their catalytic domain building blocks, has fueled an interest in harnessing the biosynthetic potential of PKSs for the microbial production of both novel natural product analogs and industrially relevant small molecules. However, a limited theoretical understanding of the determinants of PKS fold and function poses a substantial barrier to the design of active variants, and identifying strategies to reliably construct functional PKS chimeras remains an active area of research. In this work, we formalize a paradigm for the design of PKS chimeras and introduce ClusterCAD as a computational platform to streamline and simplify the process of designing experiments to test strategies for engineering PKS variants. ClusterCAD provides chemical structures with stereochemistry for the intermediates generated by each PKS module, as well as sequence- and structure-based search tools that allow users to identify modules based either on amino acid sequence or on the chemical structure of the cognate polyketide intermediate. ClusterCAD can be accessed at https://clustercad.jbei.org and at http://clustercad.igb.uci.edu.

Amplification of HMG-CoA reductase production enhances carotenoid accumulation in Neurospora crassa

Neurospora crassa, a filamentous fungus, naturally produces the carotenoids lycopene and neurosporaxanthin. To increase the carbon flux through the carotenoid biosynthetic pathway, the 1658-bp region of the HMG1 gene encoding the catalytic domain (cHMG1) of 3-hydroxy-3-methylglutaryl coenzyme A (HMG-CoA) reductase of Saccharomyces cerevisiae was expressed in N. crassa under control of the strong, constitutive glyceraldehyde-3-phosphate dehydrogenase (GPD) promoter and the inducible alcohol dehydrogenase (alcA) promoter. Overexpressing cHMG1 under control of the GPD promoter increased lycopene and neurosporaxanthin production 6- and 1.5-fold, respectively, relative to the wild-type strain. Over-expression of cHMG1 under control of the alcA promoter increased production of lycopene and neurosporaxanthin 3- and 2-fold, respectively.

Production of jet fuel precursor monoterpenoids from engineered Escherichia coli

Monoterpenes (C₁₀ isoprenoids) are the main components of essential oils and are possible precursors for many commodity chemicals and high energy density fuels. Monoterpenes are synthesized from geranyl diphosphate (GPP), which is also the precursor for the biosynthesis of farnesyl diphosphate (FPP). FPP biosynthesis diverts the carbon flux from monoterpene production to C₁₅ products and quinone biosynthesis. In this study, we tested a chromosomal mutation of Escherichia coli's native FPP synthase (IspA) to improve GPP availability for the production of monoterpenes using a heterologous mevalonate pathway. Monoterpene production at high levels required not only optimization of GPP production but also a basal level of FPP to maintain growth. The optimized strains produced two jet fuel precursor monoterpenoids 1,8-cineole and linalool at the titer of 653 mg/L and 505 mg/L, respectively, in batch cultures with 1% glucose. The engineered strains developed in this work provide useful resources for the production of high-value monoterpenes. Biotechnol. Bioeng. 2017;114: 1703-1712. © 2017 Wiley Periodicals, Inc.

A three-dimensional, stochastic simulation of biofilm growth and transport-related factors that affect structure

Biofilm structural heterogeneity affects a broad range of microbially catalysed processes. Solute transport limitation and autoinhibitor production, two factors that contribute to heterogeneous biofilm development, were investigated using BacMIST, a computer simulation model. BacMIST combines a cellular automaton algorithm for biofilm growth with Brownian diffusion for solute transport. The simulation represented the growth of microbial unit cells in a three-dimensional domain modelled after a repeating section of a constant depth film fermenter. The simulation was implemented to analyse the effects of various levels of transport limitation on a growing single-species biofilm. In a system with rapid solute diffusion, cells throughout the biofilm grew at their maximum rate, and no solute gradient was formed over the biofilm thickness. In increasingly transport-limited systems, the rapidly growing fraction of the biofilm population decreased, and was found exclusively at the biofilm–liquid interface. Trans-biofilm growth substrate gradients also deepened with increasing transport limitation. Autoinhibitory biofilm growth was simulated for various rates of microbially produced inhibitor transport. Inhibitor transport rates affected both the biofilm population dynamics and the resulting biofilm structures. The formation of networks of void spaces in slow-growing regions of the biofilm and the development of columns in the fast-growing regions suggested a possible mechanism for the microscopically observed evolution of channels in biofilms.

Whole-cell biocatalytic and de novo production of alkanes from free fatty acids in Saccharomyces cerevisiae

Rapid global industrialization in the past decades has led to extensive utilization of fossil fuels, which resulted in pressing environmental problems due to excessive carbon emission. This prompted increasing interest in developing advanced biofuels with higher energy density to substitute fossil fuels and bio‐alkane has gained attention as an ideal drop‐in fuel candidate. Production of alkanes in bacteria has been widely studied but studies on the utilization of the robust yeast host, Saccharomyces cerevisiae, for alkane biosynthesis have been lacking. In this proof‐of‐principle study, we present the unprecedented engineering of S. cerevisiae for conversion of free fatty acids to alkanes. A fatty acid α‐dioxygenase from Oryza sativa (rice) was expressed in S. cerevisiae to transform C_12–18 free fatty acids to C_11–17 aldehydes. Co‐expression of a cyanobacterial aldehyde deformylating oxygenase converted the aldehydes to the desired alkanes. We demonstrated the versatility of the pathway by performing whole‐cell biocatalytic conversion of exogenous free fatty acid feedstocks into alkanes as well as introducing the pathway into a free fatty acid overproducer for de novo production of alkanes from simple sugar. The results from this work are anticipated to advance the development of yeast hosts for alkane production. Biotechnol. Bioeng. 2017;114: 232–237. © 2016 The Authors. Biotechnology and Bioengineering Published by Wiley Periodicals, Inc.

Recent applications of synthetic biology tools for yeast metabolic engineering

The last 20 years of metabolic engineering has enabled bio-based production of fuels and chemicals from renewable carbon sources using cost-effective bioprocesses. Much of this work has been accomplished using engineered microorganisms that act as chemical factories. Although the time required to engineer microbial chemical factories has steadily decreased, improvement is still needed. Through the development of synthetic biology tools for key microbial hosts, it should be possible to further decrease the development times and improve the reliability of the resulting microorganism. Together with continuous decreases in price and improvements in DNA synthesis, assembly and sequencing, synthetic biology tools will rationalize time-consuming strain engineering, improve control of metabolic fluxes, and diversify screening assays for cellular metabolism. This review outlines some recently developed synthetic biology tools and their application to improve production of chemicals and fuels in yeast. Finally, we provide a perspective for the challenges that lie ahead.

Flux analysis of central metabolic pathways in Geobacter metallireducens during reduction of soluble Fe(III)-nitrilotriacetic acid

We analyzed the carbon fluxes in the central metabolism of Geobacter metallireducens strain GS-15 using 13C isotopomer modeling. Acetate labeled in the first or second position was the sole carbon source, and Fe-nitrilotriacetic acid was the sole terminal electron acceptor. The measured labeled acetate uptake rate was 21 mmol/g (dry weight)/h in the exponential growth phase. The resulting isotope labeling pattern of amino acids allowed an accurate determination of the in vivo global metabolic reaction rates (fluxes) through the central metabolic pathways using a computational isotopomer model. The tracer experiments showed that G. metallireducens contained complete biosynthesis pathways for essential metabolism, and this strain might also have an unusual isoleucine biosynthesis route (using acetyl coenzyme A and pyruvate as the precursors). The model indicated that over 90% of the acetate was completely oxidized to CO2 via a complete tricarboxylic acid cycle while reducing iron. Pyruvate carboxylase and phosphoenolpyruvate (PEP) carboxykinase were present under these conditions, but enzymes in the glyoxylate shunt and malic enzyme were absent. Gluconeogenesis and the pentose phosphate pathway were mainly employed for biosynthesis and accounted for less than 3% of total carbon consumption. The model also indicated surprisingly high reversibility in the reaction between oxoglutarate and succinate. This step operates close to the thermodynamic equilibrium, possibly because succinate is synthesized via a transferase reaction, and the conversion of oxoglutarate to succinate is a rate-limiting step for carbon metabolism. These findings enable a better understanding of the relationship between genome annotation and extant metabolic pathways in G. metallireducens.

Aerobic sulfide production and cadmium precipitation by Escherichia coli expressing the Treponema denticola cysteine desulfhydrase gene

The cysteine desulfhydrase gene of Treponema denticola was over-expressed in Escherichia coli to produce sulfide under aerobic conditions and to precipitate metal sulfide complexes on the cell wall. When grown in a defined salts medium supplemented with cadmium and cysteine, E. coli producing cysteine desulfhydrase secreted sulfide and removed nearly all of the cadmium from solution after 48 h. A control strain produced significantly less sulfide and removed significantly less cadmium. Measurement of acid-labile sulfide and energy dispersive X-ray spectroscopy indicated that cadmium was precipitated as cadmium sulfide. Without supplemental cysteine, both the E. coli producing cysteine desulfhydrase and the control E. coli demonstrated minimal cadmium removal.

Intracellular cellobiose metabolism and its applications in lignocellulose-based biorefineries

Complete hydrolysis of cellulose has been a key characteristic of biomass technology because of the limitation of industrial production hosts to use cellodextrin, the partial hydrolysis product of cellulose. Cellobiose, a β-1,4-linked glucose dimer, is a major cellodextrin of the enzymatic hydrolysis (via endoglucanase and exoglucanase) of cellulose. Conversion of cellobiose to glucose is executed by β-glucosidase. The complete extracellular hydrolysis of celluloses has several critical barriers in biomass technology. An alternative bioengineering strategy to make the bioprocessing less challenging is to engineer microbes with the abilities to hydrolyze and assimilate the cellulosic-hydrolysate cellodextrin. Microorganisms engineered to metabolize cellobiose rather than the monomeric glucose can provide several advantages for lignocellulose-based biorefineries. This review describes the recent advances and challenges in engineering efficient intracellular cellobiose metabolism in industrial hosts. This review also describes the limitations of and future prospectives in engineering intracellular cellobiose metabolism.

Study of nitrate stress in Desulfovibrio vulgaris Hildenborough using iTRAQ proteomics

The response of Desulfovibrio vulgaris Hildenborough (DvH), a sulphate-reducing bacterium, to nitrate stress was examined using quantitative proteomic analysis. DvH was stressed with 105 mM sodium nitrate (NaNO(3)), a level that caused a 50% inhibition in growth. The protein profile of stressed cells was compared with that of cells grown in the absence of nitrate using the iTRAQ peptide labelling strategy and tandem liquid chromatography separation coupled with mass spectrometry (quadrupole time-of-flight) detection. A total of 737 unique proteins were identified by two or more peptides, representing 22% of the total DvH proteome and spanning every functional category. The results indicate that this was a mild stress, as proteins involved in central metabolism and the sulphate reduction pathway were unperturbed. Proteins involved in the nitrate reduction pathway increased. Increases seen in transport systems for proline, glycine-betaine and glutamate indicate that the NaNO(3) exposure led to both salt stress and nitrate stress. Up-regulation observed in oxidative stress response proteins (Rbr, RbO, etc.) and a large number of ABC transport systems as well as in iron-sulphur-cluster-containing proteins, however, appear to be specific to nitrate exposure. Finally, a number of hypothetical proteins were among the most significant changers, indicating that there may be unknown mechanisms initiated upon nitrate stress in DvH.

Effect of gene location, mRNA secondary structures, and RNase sites on expression of two genes in an engineered operon

The effects of endoribonuclease sites, secondary structures in mRNA, and gene placement on protein production and mRNA stability and steady-state levels were tested in a dual-gene operon containing the genes encoding beta-galactosidase (lacZ) from Escherichia coli and green fluorescent protein (gfp) from Aequorea victoria. Two previously identified RNase E sites were placed separately between the coding regions to direct cleavage in this area and produce two secondary transcripts, each containing a single-gene coding region. Novel secondary structures were engineered into the 3' and 5' ends of each of the coding regions to protect the transcript from inactivation by endoribonucleases (5' hairpins) and degradation by exoribonucleases (3' hairpins). In addition, the effects of relative gene placement were examined by switching the locations of the two coding regions. Depending on the particular secondary structures and RNase E sites placed between the genes the relative steady-state transcript and protein levels encoded by the two reporter genes could be changed up to 2.5-fold and 4-fold, respectively. By changing gene location and incorporating secondary structures and RNase E sites the relative steady-state transcript and protein levels encoded by the two reporter genes could be changed up to 100-fold and 750-fold, respectively.

Leveraging microbial biosynthetic pathways for the generation of 'drop-in' biofuels

Advances in retooling microorganisms have enabled bioproduction of 'drop-in' biofuels, fuels that are compatible with existing spark-ignition, compression-ignition, and gas-turbine engines. As the majority of petroleum consumption in the United States consists of gasoline (47%), diesel fuel and heating oil (21%), and jet fuel (8%), 'drop-in' biofuels that replace these petrochemical sources are particularly attractive. In this review, we discuss the application of aldehyde decarbonylases to produce gasoline substitutes from fatty acid products, a recently crystallized reductase that could hydrogenate jet fuel precursors from terpene synthases, and the exquisite control of polyketide synthases to produce biofuels with desired physical properties (e.g., lower freezing points). With our increased understanding of biosynthetic logic of metabolic pathways, we discuss the unique advantages of fatty acid, terpene, and polyketide synthases for the production of bio-based gasoline, diesel and jet fuel.

An Orthogonal and pH-Tunable Sensor-Selector for Muconic Acid Biosynthesis in Yeast

Microbes offer enormous potential for production of industrially relevant chemicals and therapeutics, yet the rapid identification of high-producing microbes from large genetic libraries is a major bottleneck in modern cell factory development. Here, we develop and apply a synthetic selection system in Saccharomyces cerevisiae that couples the concentration of muconic acid, a plastic precursor, to cell fitness by using the prokaryotic transcriptional regulator BenM driving an antibiotic resistance gene. We show that the sensor-selector does not affect production nor fitness, and find that tuning pH of the cultivation medium limits the rise of nonproducing cheaters. We apply the sensor-selector to selectively enrich for best-producing variants out of a large library of muconic acid production strains, and identify an isolate that produces more than 2 g/L muconic acid in a bioreactor. We expect that this sensor-selector can aid the development of other synthetic selection systems based on allosteric transcription factors.

Genome-scale metabolic rewiring improves titers rates and yields of the non-native product indigoidine at scale

High titer, rate, yield (TRY), and scalability are challenging metrics to achieve due to trade-offs between carbon use for growth and production. To achieve these metrics, we take the minimal cut set (MCS) approach that predicts metabolic reactions for elimination to couple metabolite production strongly with growth. We compute MCS solution-sets for a non-native product indigoidine, a sustainable pigment, in Pseudomonas putida KT2440, an emerging industrial microbe. From the 63 solution-sets, our omics guided process identifies one experimentally feasible solution requiring 14 simultaneous reaction interventions. We implement a total of 14 genes knockdowns using multiplex-CRISPRi. MCS-based solution shifts production from stationary to exponential phase. We achieve 25.6 g/L, 0.22 g/l/h, and ~50% maximum theoretical yield (0.33 g indigoidine/g glucose). These phenotypes are maintained from batch to fed-batch mode, and across scales (100-ml shake flasks, 250-ml ambr®, and 2-L bioreactors).

Development of a Transcription Factor-Based Lactam Biosensor

Lactams are an important class of commodity chemicals used in the manufacture of nylons, with millions of tons produced every year. Biological production of lactams could be greatly improved by high-throughput sensors for lactam biosynthesis. To identify biosensors of lactams, we applied a chemoinformatic approach inspired by small molecule drug discovery. We define this approach as analogue generation toward catabolizable chemicals or AGTC. We discovered a lactam biosensor based on the ChnR/Pb transcription factor-promoter pair. The microbial biosensor is capable of sensing ε-caprolactam, δ-valerolactam, and butyrolactam in a dose-dependent manner. The biosensor has sufficient specificity to discriminate against lactam biosynthetic intermediates and therefore could potentially be applied for high-throughput metabolic engineering for industrially important high titer lactam biosynthesis.

Catabolite repression of the propionate catabolic genes in Escherichia coli and Salmonella enterica: evidence for involvement of the cyclic AMP receptor protein

Previous studies with Salmonella enterica serovar Typhimurium LT2 demonstrated that transcriptional activation of the prpBCDE operon requires the function of transcription factor PrpR, sigma-54, and IHF. In this study, we found that transcription from the prpBCDE and prpR promoters was down-regulated by the addition of glucose or glycerol, indicating that these genes may be regulated by the cyclic AMP (cAMP)-cAMP receptor protein (CRP) complex. Targeted mutagenesis of a putative CRP-binding site in the promoter region between prpR and prpBCDE suggested that these genes are under the control of CRP. Furthermore, cells with defects in cya or crp exhibited reduced transcriptional activation of prpR and prpBCDE in Escherichia coli. These results demonstrate that propionate metabolism is subject to catabolite repression by the global transcriptional regulator CRP and that this regulation is effected through control of both the regulator gene prpR and the prpBCDE operon itself. The unique properties of the regulation of these two divergent promoters may have important implications for mechanisms of CRP-dependent catabolite repression acting in conjunction with a member of the sigma-54 family of transcriptional activators.