Enhanced dynamic flux variability analysis for improving growth and production rate in microbial strains

Metabolic engineering is highly demanded currently for the production of various useful compounds such as succinate and lactate that are very useful in food, pharmaceutical, fossil fuels, and energy industries. Gene or reaction deletion known as knockout is one of the strategies used in in silico metabolic engineering to change the metabolism of the chosen microbial cells to obtain the desired phenotypes. However, the size and complexity of the metabolic network are a challenge in determining the near-optimal set of genes to be knocked out in the metabolism due to the presence of competing pathway that interrupts the high production of desired metabolite, leading to low production rate and growth rate of the required microorganisms. In addition, the inefficiency of existing algorithms in reconstructing high growth rate and production rate becomes one of the issues to be solved. Therefore, this research proposes Dynamic Flux Variability Analysis (DFVA) algorithm to identify the best knockout reaction combination to improve the production of desired metabolites in microorganisms. Based on the experimental results, DFVA shows an improvement of growth rate of succinate and lactate by 12.06% and 47.16% respectively in E. coli and by 4.62% and 47.98% respectively in S. Cerevisae. Suggested reactions to be knocked out to improve the production of succinate and lactate have been identified and validated through the biological database

Efficient solution of ordinary differential equations with a parametric lexicographic linear program embedded

This work analyzes the initial value problem in ordinary differential equations with a parametric lexicographic linear program (LP) embedded. The LP is said to be embedded since the dynamics depend on the solution of the LP, which is in turn parameterized by the dynamic states. This problem formulation finds application in dynamic flux balance analysis, which serves as a modeling framework for industrial fermentation reactions. It is shown that the problem formulation can be intractable numerically, which arises from the fact that the LP induces an effective domain that may not be open. A numerical method is developed which reformulates the system so that it is defined on an open set. The result is a system of semi-explicit index-one differential algebraic equations, which can be solved with efficient and accurate methods. It is shown that this method addresses many of the issues stemming from the original problem’s intractability. The application of the method to examples of industrial fermentation processes demonstrates its effectiveness and efficiency

Genome-scale metabolic modeling of cyanbacteria: network structure, interactions, reconstruction and dynamics

2016 Fall.Includes bibliographical references.Metabolic network modeling, a field of systems biology and bioengineering, enhances the quantitative predictive understanding of cellular metabolism and thereby assists in the development of model-guided metabolic engineering strategies. Metabolic models use genome-scale network reconstructions, and combine it with mathematical methods for quantitative prediction. Metabolic system reconstructions, contain information on genes, enzymes, reactions, and metabolites, and are converted into two types of networks: (i) gene-enzyme-reaction, and (ii) reaction-metabolite. The former details the links between the genes that are known to code for metabolic enzymes, and the reaction pathways that the enzymes participate in. The latter details the chemical transformation of metabolites, step by step, into biomass and energy. The latter network is transformed into a system of equations and simulated using different methods. Prominent among these are constraint-based methods, especially Flux Balance Analysis, which utilizes linear programming tools to predict intracellular fluxes of single cells. Over the past 25 years, metabolic network modeling has had a range of applications in the fields of model-driven discovery, prediction of cellular phenotypes, analysis of biological network properties, multi-species interactions, engineering of microbes for product synthesis, and studying evolutionary processes. This thesis is concerned with the development and application of metabolic network modeling to cyanobacteria as well as E. coli. Chapter 1 is a brief survey of the past, present, and future of constraint-based modeling using flux balance analysis in systems biology. It includes discussion of (i) formulation, (ii) assumption, (iii) variety, (iv) availability, and (v) future directions in the field of constraint based modeling. Chapter 2, explores the enzyme-reaction networks of metabolic reconstructions belonging to various organisms; and finds that the distribution of the number of reactions an enzyme participates in, i.e. the enzyme-reaction distribution, is surprisingly similar. The role of this distribution in the robustness of the organism is also explored. Chapter 3, applies flux balance analysis on models of E. coli, Synechocystis sp. PCC6803, and C. reinhardtii to understand epistatic interactions between metabolic genes and pathways. We show that epistatic interactions are dependent on the environmental conditions, i.e. carbon source, carbon/oxygen ratio in E. coli, and light intensity in Synechocystis sp. PCC6803 and C. reinhardtii. Cyanobacteria are photosynthetic organisms and have great potential for metabolic engineering to produce commercially important chemicals such as biofuels, pharmaceuticals, and nutraceuticals. Chapter 4 presents our new genome scale reconstruction of the model cyanobacterium, Synechocystis sp. PCC6803, called iCJ816. This reconstruction was analyzed and compared to experimental studies, and used for predicting the capacity of the organism for (i) carbon dioxide remediation, and (ii) production of intracellular chemical species. Chapter 5 uses our new model iCJ816 for dynamic analysis under diurnal growth simulations. We discuss predictions of different optimization schemes, and present a scheme that qualitatively matches observations