Anthony P. Burgard

Metabolic engineering of Escherichia coli for direct production of 1,4-butanediol

1,4-Butanediol (BDO) is an important commodity chemical used to manufacture over 2.5 million tons annually of valuable polymers, and it is currently produced exclusively through feedstocks derived from oil and natural gas. Herein we report what are to our knowledge the first direct biocatalytic routes to BDO from renewable carbohydrate feedstocks, leading to a strain of Escherichia coli capable of producing 18 g l(-1) of this highly reduced, non-natural chemical. A pathway-identification algorithm elucidated multiple pathways for the biosynthesis of BDO from common metabolic intermediates. Guided by a genome-scale metabolic model, we engineered the E. coli host to enhance anaerobic operation of the oxidative tricarboxylic acid cycle, thereby generating reducing power to drive the BDO pathway. The organism produced BDO from glucose, xylose, sucrose and biomass-derived mixed sugar streams. This work demonstrates a systems-based metabolic engineering approach to strain design and development that can enable new bioprocesses for commodity chemicals that are not naturally produced by living cells.

Minimal Reaction Sets for Escherichia coli Metabolism under Different Growth Requirements and Uptake Environments

A computational procedure for identifying the minimal set of metabolic reactions capable of supporting various growth rates on different substrates is introduced and applied to a flux balance model of the Escherichia colimetabolic network. This task is posed mathematically as a generalized network optimization problem. The minimal reaction sets capable of supporting specified growth rates are determined for two different uptake conditions: (i) limiting the uptake of organic material to a single organic component (e.g., glucose or acetate) and (ii) allowing the importation of any metabolite with available cellular transport reactions. We find that minimal reaction network sets are highly dependent on the uptake environment and the growth requirements imposed on the network. Specifically, we predict that the E. coli network, as described by the flux balance model, requires 224 metabolic reactions to support growth on a glucose‐only medium and 229 for an acetate‐only medium, while only 122 reactions enable growth on a specially engineered growth medium.

Geobacter sulfurreducens strain engineered for increased rates of respiration

Geobacter species are among the most effective microorganisms known for the bioremediation of radioactive and toxic metals in contaminated subsurface environments and for converting organic compounds to electricity in microbial fuel cells. However, faster rates of electron transfer could aid in optimizing these processes. Therefore, the Optknock strain design methodology was applied in an iterative manner to the constraint-based, in silico model of Geobacter sulfurreducens to identify gene deletions predicted to increase respiration rates. The common factor in the Optknock predictions was that each resulted in a predicted increase in the cellular ATP demand, either by creating ATP-consuming futile cycles or decreasing the availability of reducing equivalents and inorganic phosphate for ATP biosynthesis. The in silico model predicted that increasing the ATP demand would result in higher fluxes of acetate through the TCA cycle and higher rates of NADPH oxidation coupled with decreases in flux in reactions that funnel acetate toward biosynthetic pathways. A strain of G. sulfurreducens was constructed in which the hydrolytic, F(1) portion of the membrane-bound F(0)F(1) (H(+))-ATP synthase complex was expressed when IPTG was added to the medium. Induction of the ATP drain decreased the ATP content of the cell by more than half. The cells with the ATP drain had higher rates of respiration, slower growth rates, and a lower cell yield. Genome-wide analysis of gene transcript levels indicated that when the higher rate of respiration was induced transcript levels were higher for genes involved in energy metabolism, especially in those encoding TCA cycle enzymes, subunits of the NADH dehydrogenase, and proteins involved in electron acceptor reduction. This was accompanied by lower transcript levels for genes encoding proteins involved in amino acid biosynthesis, cell growth, and motility. Several changes in gene expression that involve processes not included in the in silico model were also detected, including increased expression of a number of redox-active proteins, such as c-type cytochromes and a putative multicopper outer-surface protein. The results demonstrate that it is possible to genetically engineer increased respiration rates in G. sulfurreducens in accordance with predictions from in silico metabolic modeling. To our knowledge, this is the first report of metabolic engineering to increase the respiratory rate of a microorganism.

Applications of metabolic modeling to drive bioprocess development for the production of value-added chemicals

Increasing numbers of value added chemicals are being produced using microbial fermentation strategies. Computational modeling and simulation of microbial metabolism is rapidly becoming an enabling technology that is driving a new paradigm to accelerate the bioprocess development cycle. In particular, constraint-based modeling and the development of genome-scale models of industrial microbes are finding increasing utility across many phases of the bioprocess development workflow. Herein, we review and discuss the requirements and trends in the industrial application of this technology as we build toward integrated computational/experimental platforms for bioprocess engineering. Specifically we cover the following topics: (1) genome-scale models as genetically and biochemically consistent representations of metabolic networks; (2) the ability of these models to predict, assess, and interpret metabolic physiology and flux states of metabolism; (3) the model-guided integrative analysis of high throughput ‘omics’ data; (4) the reconciliation and analysis of on- and off-line fermentation data as well as flux tracing data; (5) model-aided strain design strategies and the integration of calculated biotransformation routes; and (6) control and optimization of the fermentation processes. Collectively, constraint-based modeling strategies are impacting the iterative characterization of metabolic flux states throughout the bioprocess development cycle, while also driving metabolic engineering strategies and fermentation optimization.

Development of a commercial scale process for production of 1,4-butanediol from sugar.

A sustainable bioprocess for the production of 1,4-butanediol (BDO) from carbohydrate feedstocks was developed. BDO is a chemical intermediate that goes into a variety of products including automotive parts, electronics, and apparel, and is currently manufactured commercially through energy-intensive petrochemical processes using fossil raw materials. This review highlights the development of an Escherichia coli strain and an overall process that successfully performed at commercial scale for direct production of bio-BDO from dextrose. Achieving such high level performance required an integrated technology platform enabling detailed engineering of enzyme, pathway, metabolic network, and organism, as well as development of effective fermentation and downstream recovery processes.

An integrated biotechnology platform for developing sustainable chemical processes.

Genomatica has established an integrated computational/experimental metabolic engineering platform to design, create, and optimize novel high performance organisms and bioprocesses. Here we present our platform and its use to develop E. coli strains for production of the industrial chemical 1,4-butanediol (BDO) from sugars. A series of examples are given to demonstrate how a rational approach to strain engineering, including carefully designed diagnostic experiments, provided critical insights about pathway bottlenecks, byproducts, expression balancing, and commercial robustness, leading to a superior BDO production strain and process.

Do genome-scale models need exact solvers or clearer standards?

Constraint‐based analysis of genome‐scale models (GEMs) arose shortly after the first genome sequences became available. As numerous reviews of the field show, this approach and methodology has proven to be successful in studying a wide range of biological phenomena (McCloskey et al, 2013; Bordbar et al, 2014). However, efforts to expand the user base are impeded by hurdles in correctly formulating these problems to obtain numerical solutions. In particular, in a study entitled “An exact arithmetic toolbox for a consistent and reproducible structural analysis of metabolic network models” (Chindelevitch et al, 2014), the authors apply an exact solver to 88 genome‐scale constraint‐based models of metabolism. The authors claim that COBRA calculations (Orth et al, 2010) are inconsistent with their results and that many published and actively used (Lee et al, 2007; McCloskey et al, 2013) genome‐scale models do support cellular growth in existing studies only because of numerical errors. They base these broad claims on two observations: (i) three reconstructions (iAF1260, iIT341, and iNJ661) compute feasibly in COBRA, but are infeasible when exact numerical algorithms are used by their software (entitled MONGOOSE); (ii) linear programs generated by MONGOOSE for iIT341 were submitted to the NEOS Server (a Web site that runs linear programs through various solvers) and gave inconsistent results. They further claim that a large percentage of these COBRA models are actually unable to produce biomass flux. Here, we demonstrate that the claims made by Chindelevitch et al (2014) stem from an incorrect parsing of models from files rather than actual problems with numerical error or COBRA computations.

Identification of metabolic engineering targets for the enhancement of 1,4-butanediol production in recombinant E. coli using large-scale kinetic models

Rational metabolic engineering methods are increasingly employed in designing the commercially viable processes for the production of chemicals relevant to pharmaceutical, biotechnology, and food and beverage industries. With the growing availability of omics data and of methodologies capable to integrate the available data into models, mathematical modeling and computational analysis are becoming important in designing recombinant cellular organisms and optimizing cell performance with respect to desired criteria. In this contribution, we used the computational framework ORACLE (Optimization and Risk Analysis of Complex Living Entities) to analyze the physiology of recombinant Escherichia coli producing 1,4-butanediol (BDO) and to identify potential strategies for improved production of BDO. The framework allowed us to integrate data across multiple levels and to construct a population of large-scale kinetic models despite the lack of available information about kinetic properties of every enzyme in the metabolic pathways. We analyzed these models and we found that the enzymes that primarily control the fluxes leading to BDO production are part of central glycolysis, the lower branch of tricarboxylic acid (TCA) cycle and the novel BDO production route. Interestingly, among the enzymes between the glucose uptake and the BDO pathway, the enzymes belonging to the lower branch of TCA cycle have been identified as the most important for improving BDO production and yield. We also quantified the effects of changes of the target enzymes on other intracellular states like energy charge, cofactor levels, redox state, cellular growth, and byproduct formation. Independent earlier experiments on this strain confirmed that the computationally obtained conclusions are consistent with the experimentally tested designs, and the findings of the present studies can provide guidance for future work on strain improvement. Overall, these studies demonstrate the potential and effectiveness of ORACLE for the accelerated design of microbial cell factories.