Anupam Chowdhury

An integrated computational and experimental study for overproducing fatty acids in Escherichia coli

Increasing demands for petroleum have stimulated sustainable ways to produce chemicals and biofuels. Specifically, fatty acids of varying chain lengths (C₆-C₁₆) naturally synthesized in many organisms are promising starting points for the catalytic production of industrial chemicals and diesel-like biofuels. However, bio-production of fatty acids from plants and other microbial production hosts relies heavily on manipulating tightly regulated fatty acid biosynthetic pathways. In addition, precursors for fatty acids are used along other central metabolic pathways for the production of amino acids and biomass, which further complicates the engineering of microbial hosts for higher yields. Here, we demonstrate an iterative metabolic engineering effort that integrates computationally driven predictions and metabolic flux analysis techniques to meet this challenge. The OptForce procedure was used for suggesting and prioritizing genetic manipulations that overproduce fatty acids of different chain lengths from C₆ to C₁₆ starting with wild-type E. coli. We identified some common but mostly chain-specific genetic interventions alluding to the possibility of fine-tuning overproduction for specific fatty acid chain lengths. In accordance with the OptForce prioritization of interventions, fabZ and acyl-ACP thioesterase were upregulated and fadD was deleted to arrive at a strain that produces 1.70 g/L and 0.14 g fatty acid/g glucose (∼39% maximum theoretical yield) of C₁₄₋₁₆ fatty acids in minimal M9 medium. These results highlight the benefit of using computational strain design and flux analysis tools in the design of recombinant strains of E. coli to produce free fatty acids.

k-OptForce: Integrating Kinetics with Flux Balance Analysis for Strain Design

Computational strain design protocols aim at the system-wide identification of intervention strategies for the enhanced production of biochemicals in microorganisms. Existing approaches relying solely on stoichiometry and rudimentary constraint-based regulation overlook the effects of metabolite concentrations and substrate-level enzyme regulation while identifying metabolic interventions. In this paper, we introduce k-OptForce, which integrates the available kinetic descriptions of metabolic steps with stoichiometric models to sharpen the prediction of intervention strategies for improving the bio-production of a chemical of interest. It enables identification of a minimal set of interventions comprised of both enzymatic parameter changes (for reactions with available kinetics) and reaction flux changes (for reactions with only stoichiometric information). Application of k-OptForce to the overproduction of L-serine in E. coli and triacetic acid lactone (TAL) in S. cerevisiae revealed that the identified interventions tend to cause less dramatic rearrangements of the flux distribution so as not to violate concentration bounds. In some cases the incorporation of kinetic information leads to the need for additional interventions as kinetic expressions render stoichiometry-only derived interventions infeasible by violating concentration bounds, whereas in other cases the kinetic expressions impart flux changes that favor the overproduction of the target product thereby requiring fewer direct interventions. A sensitivity analysis on metabolite concentrations shows that the required number of interventions can be significantly affected by changing the imposed bounds on metabolite concentrations. Furthermore, k-OptForce was capable of finding non-intuitive interventions aiming at alleviating the substrate-level inhibition of key enzymes in order to enhance the flux towards the product of interest, which cannot be captured by stoichiometry-alone analysis. This study paves the way for the integrated analysis of kinetic and stoichiometric models and enables elucidating system-wide metabolic interventions while capturing regulatory and kinetic effects.

Recent advances in the reconstruction of metabolic models and integration of omics data

With the ever-accelerating pace of genome sequencing and annotation information generation, the development of computational pipelines for the rapid reconstruction of high-quality metabolic networks has received significant attention. Herein, we review the available biological databases and automated/semi-automated reconstruction tools. In addition, we describe available methodologies for the integration of high-throughput omics data to increase metabolic phenotype prediction accuracy. Data heterogeneity and lack of better integration of metabolic reconstruction pipelines with omics data generation protocols have hampered rapid progress thus far.

Systems Metabolic Engineering Design: Fatty Acid Production as an Emerging Case Study

Increasing demand for petroleum has stimulated industry to develop sustainable production of chemicals and biofuels using microbial cell factories. Fatty acids of chain lengths from C6 to C16 are propitious intermediates for the catalytic synthesis of industrial chemicals and diesel-like biofuels. The abundance of genetic information available for Escherichia coli and specifically, fatty acid metabolism in E. coli, supports this bacterium as a promising host for engineering a biocatalyst for the microbial production of fatty acids. Recent successes rooted in different features of systems metabolic engineering in the strain design of high-yielding medium chain fatty acid producing E. coli strains provide an emerging case study of design methods for effective strain design. Classical metabolic engineering and synthetic biology approaches enabled different and distinct design paths towards a high-yielding strain. Here we highlight a rational strain design process in systems biology, an integrated computational and experimental approach for carboxylic acid production, as an alternative method. Additional challenges inherent in achieving an optimal strain for commercialization of medium chain-length fatty acids will likely require a collection of strategies from systems metabolic engineering. Not only will the continued advancement in systems metabolic engineering result in these highly productive strains more quickly, this knowledge will extend more rapidly the carboxylic acid platform to the microbial production of carboxylic acids with alternate chain-lengths and functionalities.

Succinate Overproduction: A Case Study of Computational Strain Design Using a Comprehensive Escherichia coli Kinetic Model

Computational strain-design prediction accuracy has been the focus for many recent efforts through the selective integration of kinetic information into metabolic models. In general, kinetic model prediction quality is determined by the range and scope of genetic and/or environmental perturbations used during parameterization. In this effort, we apply the k-OptForce procedure on a kinetic model of E. coli core metabolism constructed using the Ensemble Modeling (EM) method and parameterized using multiple mutant strains data under aerobic respiration with glucose as the carbon source. Minimal interventions are identified that improve succinate yield under both aerobic and anaerobic conditions to test the fidelity of model predictions under both genetic and environmental perturbations. Under aerobic condition, k-OptForce identifies interventions that match existing experimental strategies while pointing at a number of unexplored flux re-directions such as routing glyoxylate flux through the glycerate metabolism to improve succinate yield. Many of the identified interventions rely on the kinetic descriptions that would not be discoverable by a purely stoichiometric description. In contrast, under fermentative (anaerobic) condition, k-OptForce fails to identify key interventions including up-regulation of anaplerotic reactions and elimination of competitive fermentative products. This is due to the fact that the pathways activated under anaerobic condition were not properly parameterized as only aerobic flux data were used in the model construction. This study shed light on the importance of condition-specific model parameterization and provides insight on how to augment kinetic models so as to correctly respond to multiple environmental perturbations.

Improving prediction fidelity of cellular metabolism with kinetic descriptions

Several modeling frameworks for describing and redirecting cellular metabolism have been developed keeping pace with the rapid development in high-throughput data generation and advances in metabolic engineering techniques. The incorporation of kinetic information within stoichiometry-only modeling techniques offers potential advantages for improved phenotype prediction and consequently more precise computational strain design. In addition to substrate-level kinetic regulatory information, the integration of a number of additional layers of regulation at the transcription, translation, and post-translation levels is sought after by many research groups. However, the practical integration of these complex biological processes into a unified framework amenable to design remains an ongoing challenge.

Designing overall stoichiometric conversions and intervening metabolic reactions.

Existing computational tools for de novo metabolic pathway assembly, either based on mixed integer linear programming techniques or graph-search applications, generally only find linear pathways connecting the source to the target metabolite. The overall stoichiometry of conversion along with alternate co-reactant (or co-product) combinations is not part of the pathway design. Therefore, global carbon and energy efficiency is in essence fixed with no opportunities to identify more efficient routes for recycling carbon flux closer to the thermodynamic limit. Here, we introduce a two-stage computational procedure that both identifies the optimum overall stoichiometry (i.e., optStoic) and selects for (non-)native reactions (i.e., minRxn/minFlux) that maximize carbon, energy or price efficiency while satisfying thermodynamic feasibility requirements. Implementation for recent pathway design studies identified non-intuitive designs with improved efficiencies. Specifically, multiple alternatives for non-oxidative glycolysis are generated and non-intuitive ways of co-utilizing carbon dioxide with methanol are revealed for the production of C2+ metabolites with higher carbon efficiency.

Using Gene Essentiality and Synthetic Lethality Information to Correct Yeast and CHO Cell Genome-Scale Models.

Essentiality (ES) and Synthetic Lethality (SL) information identify combination of genes whose deletion inhibits cell growth. This information is important for both identifying drug targets for tumor and pathogenic bacteria suppression and for flagging and avoiding gene deletions that are non-viable in biotechnology. In this study, we performed a comprehensive ES and SL analysis of two important eukaryotic models (S. cerevisiae and CHO cells) using a bilevel optimization approach introduced earlier. Information gleaned from this study is used to propose specific model changes to remedy inconsistent with data model predictions. Even for the highly curated Yeast 7.11 model we identified 50 changes (metabolic and GPR) leading to the correct prediction of an additional 28% of essential genes and 36% of synthetic lethals along with a 53% reduction in the erroneous identification of essential genes. Due to the paucity of mutant growth phenotype data only 12 changes were made for the CHO 1.2 model leading to an additional correctly predicted 11 essential and eight non-essential genes. Overall, we find that CHO 1.2 was 76% less accurate than the Yeast 7.11 metabolic model in predicting essential genes. Based on this analysis, 14 (single and double deletion) maximally informative experiments are suggested to improve the CHO cell model by using information from a mouse metabolic model. This analysis demonstrates the importance of single and multiple knockout phenotypes in assessing and improving model reconstructions. The advent of techniques such as CRISPR opens the door for the global assessment of eukaryotic models.

Advances in de novo strain design using integrated systems and synthetic biology tools.

Recent efforts in expanding the range of biofuel and biorenewable molecules using microbial production hosts have focused on the introduction of non-native pathways in model organisms and the bio-prospecting of non-model organisms with desirable features. Current challenges lie in the assembly and coordinated expression of the (non-)native pathways and the elimination of competing pathways and undesirable regulation. Several systems and synthetic biology approaches providing contrasting top-down and bottom-up strategies, respectively, have been developed. In this review, we discuss recent advances in both in silico and experimental approaches for metabolic pathway design and engineering, with a critical assessment of their merits and remaining challenges.

Multilevel engineering of the upstream module of aromatic amino acid biosynthesis in Saccharomyces cerevisiae for high production of polymer and drug precursors

A multilevel approach was implemented in Saccharomyces cerevisiae to optimize the precursor module of the aromatic amino acid biosynthesis pathway, which is a rich resource for synthesizing a great variety of chemicals ranging from polymer precursor, to nutraceuticals and pain-relief drugs. To facilitate the discovery of novel targets to enhance the pathway flux, we incorporated the computational tool YEASTRACT for predicting novel transcriptional repressors and OptForce strain-design for identifying non-intuitive pathway interventions. The multilevel approach consisted of (i) relieving the pathway from strong transcriptional repression, (ii) removing competing pathways to ensure high carbon capture, and (iii) rewiring precursor pathways to increase the carbon funneling to the desired target. The combination of these interventions led to the establishment of a S. cerevisiae strain with shikimic acid (SA) titer reaching as high as 2.5gL-1, 7-fold higher than the base strain. Further expansion of the platform led to the titer of 2.7gL-1 of muconic acid (MA) and its intermediate protocatechuic acid (PCA) together. Both the SA and MA production platforms demonstrated increases in titer and yield nearly 300% from the previously reported, highest-producing S. cerevisiae strains. Further examination elucidated the diverged impacts of disrupting the oxidative branch (ZWF1) of the pentose phosphate pathway on the titers of desired products belonging to different portions of the pathway. The investigation of other non-intuitive interventions like the deletion of the Pho13 enzyme also revealed the important role of the transaldolase in determining the fate of the carbon flux in the pathways of study. This integrative approach identified novel determinants at both transcriptional and metabolic levels that constrain the flux entering the aromatic amino acid pathway. In the future, this platform can be readily used for engineering the downstream modules toward the production of important plant-sourced aromatic secondary metabolites.