Reconstruction of a catalogue of genome-scale metabolic models with enzymatic constraints using GECKO 2.0

Genome-scale metabolic models (GEMs) have been widely used for quantitative exploration of the relation between genotype and phenotype. Streamlined integration of enzyme constraints and proteomics data into such models was first enabled by the GECKO toolbox, allowing the study of phenotypes constrained by protein limitations. Here, we upgrade the toolbox in order to enhance models with enzyme and proteomics constraints for any organism with a compatible GEM reconstruction. With this, enzyme-constrained models for the budding yeasts Saccharomyces cerevisiae, Yarrowia lipolytica and Kluyveromyces marxianus are generated to study their long-term adaptation to several stress factors by incorporation of proteomics data. Predictions reveal that upregulation and high saturation of enzymes in amino acid metabolism are common across organisms and conditions, suggesting the relevance of metabolic robustness in contrast to optimal protein utilization as a cellular objective for microbial growth under stress and nutrient-limited conditions. The functionality of GECKO is expanded with an automated framework for continuous and version-controlled update of enzyme-constrained GEMs, also producing such models for Escherichia coli and Homo sapiens. In this work, we facilitate the utilization of enzyme-constrained GEMs in basic science, metabolic engineering and synthetic biology purposes

METABOLIC MODELING AND OMICS-INTEGRATIVE ANALYSIS OF SINGLE AND MULTI-ORGANISM SYSTEMS: DISCOVERY AND REDESIGN

Computations and modeling have emerged as indispensable tools that drive the process of understanding, discovery, and redesign of biological systems. With the accelerating pace of genome sequencing and annotation information generation, the development of computational pipelines for the rapid reconstruction of high-quality genome-scale metabolic networks has received significant attention. These models provide a rich tapestry for computational tools to quantitatively assess the metabolic phenotypes for various systems-level studies and to develop engineering interventions at the DNA, RNA, or enzymatic level by careful tuning in the biophysical modeling frameworks. in silico genome-scale metabolic modeling algorithms based on the concept of optimization, along with the incorporation of multi-level omics information, provides a diverse array of toolboxes for new discovery in the metabolism of living organisms (which includes single-cell microbes, plants, animals, and microbial ecosystems) and allows for the reprogramming of metabolism for desired output(s). Throughout my doctoral research, I used genome-scale metabolic models and omics-integrative analysis tools to study how microbes, plants, animal, and microbial ecosystems respond or adapt to diverse environmental cues, and how to leverage the knowledge gleaned from that to answer important biological questions. Each chapter in this dissertation will provide a detailed description of the methodology, results, and conclusions from one specific research project. The research works presented in this dissertation represent important foundational advance in Systems Biology and are crucial for sustainable development in food, pharmaceuticals and bioproduction of the future. Advisor: Rajib Sah

Creation and Application of Various Tools for the Reconstruction, Curation, and Analysis of Genome-Scale Models of Metabolism

Systems biology uses mathematics tools, modeling, and analysis for holistic understanding and design of biological systems, allowing the investigation of metabolism and the generation of actionable hypotheses based on model analyses. Detailed here are several systems biology tools for model reconstruction, curation, analysis, and application through synthetic biology. The first, OptFill, is a holistic (whole model) and conservative (minimizing change) tool to aid in genome-scale model (GSM) reconstructions by filling metabolic gaps caused by lack of system knowledge. This is accomplished through Mixed Integer Linear Programming (MILP), one step of which may also be independently used as an additional curation tool. OptFill is applied to a GSM reconstruction of the melanized fungus Exophiala dermatitidis, which underwent various analyses investigating pigmentogenesis and similarity to human melanogenesis. Analysis suggest that carotenoids serve a currently unknown function in E. dermatitidis and that E. dermatitidis could serve as a model of human melanocytes for biomedical applications. Next, a new approach to dynamic Flux Balance Analysis (dFBA) is detailed, the Optimization- and Runge-Kutta- based Approach (ORKA). The ORKA is applied to the model plant Arabidopsis thaliana to show its ability to recreate in vivo observations. The analyzed model is more detailed than previous models, encompassing a larger time scale, modeling more tissues, and with higher accuracy. Finally, a pair of tools, the Eukaryotic Genetic Circuit Design (EuGeneCiD) and Modeling (EuGeneCiM) tools, is introduced which can aid in the design and modeling of synthetic biology applications hypothesized using systems biology. These tools bring a computational approach to synthetic biology, and are applied to Arabidopsis thaliana to design thousands of potential two-input genetic circuits which satisfy 27 different input and logic gate combinations. EuGeneCiM is further used to model a repressilator circuit. Efforts are ongoing to disseminate these tools to maximize their impact on the field of systems biology. Future research will include further investigation of E. dermatitidis through modeling and expanding my expertise to kinetic models of metabolism. Advisor: Rajib Sah

Yeast Systems Biotechnology for Production of Value-Added Biochemicals

Ph.DDOCTOR OF PHILOSOPH

Understanding and engineering cofactor metabolism: a quest for improved biochemical production

In the field of metabolic engineering, where cells are treated as “factories” that synthesise industrial compounds, it is essential for cell metabolism to accommodate the energy and re- dox cofactor demands of synthetic pathways. A balanced supply and consumption of ATP and NAD(P)H directly influences biotechnological performance. This study develops computational and experimental frameworks to explore how ATP and NAD(P)H limit the yield of synthetic pathways during bioproduction. Constraint-based modelling was used to develop a novel computational protocol, CBA (Co- factor Balance Assessment), which tracks how ATP and NAD(P)H contribute to cell target production, as opposed to cell maintenance, biomass and waste release, in the presence of a synthetic pathway. Using butanol pathways (a non-native product in E.coli) with varying cofactor demands, CBA discerned the network-wide effects of cofactor variations on butanol yield. Results indicate that yields could be boosted by up to 13% if the introduced pathway is balanced both in terms of energy and redox. CBA simplified cofactor balance assessments and provided insights into how to improve the efficiency of recombinant strains. Physiological and metabolic responses to cofactor perturbations were also experimentally assessed. The predominant phenotypes of strains harbouring the ATP synthase and PCK knockouts included high glycolytic flux, lower biomass and ATP. These strains were used to improve ethanol production, resulting in yields 10% and 29% higher than the WT overproducing ethanol and reaching over 70% of the theoretical maximum. The low but positive ATP yields boosted ethanol production and minimised unrestricted growth. This research posits that early-stage in silico cofactor usage profiling serves as an instrument to select better performing pathways. Significant yield improvements can be achieved experimentally with a small number of cofactor-driven modifications that reduce the waste of cofactors, illustrating the potential of these strains as platforms to improve bioproduction of cofactor-neutral or cofactor-surplus synthetic pathways.Open Acces

A systems biology understanding of protein constraints in the metabolism of budding yeasts

Fermentation technologies, such as bread making and production of alcoholic beverages, have been crucial for development of humanity throughout history. Saccharomyces cerevisiae provides a natural platform for this, due to its capability to transform sugars into ethanol. This, and other yeasts, are now used for production of pharmaceuticals, including insulin and artemisinic acid, flavors, fragrances, nutraceuticals, and fuel precursors. In this thesis, different systems biology methods were developed to study interactions between metabolism, enzymatic capabilities, and regulation of gene expression in budding yeasts. In paper I, a study of three different yeast species (S. cerevisiae, Yarrowia lipolytica and Kluyveromyces marxianus), exposed to multiple conditions, was carried out to understand their adaptation to environmental stress. Paper II revises the use of genome-scale metabolic models (GEMs) for the study and directed engineering of diverse yeast species. Additionally, 45 GEMs for different yeasts were collected, analyzed, and tested. In paper III, GECKO 2.0, a toolbox for integration of enzymatic constraints and proteomics data into GEMs, was developed and used for reconstruction of enzyme-constrained models (ecGEMs) for three yeast species and model organisms. Proteomics data and ecGEMs were used to further characterize the impact of environmental stress over metabolism of budding yeasts. On paper IV, gene engineering targets for increased accumulation of heme in S. cerevisiae cells were predicted with an ecGEM. Predictions were experimentally validated, yielding a 70-fold increase in intracellular heme. The prediction method was systematized and applied to the production of 102 chemicals in S. cerevisiae (Paper V). Results highlighted general principles for systems metabolic engineering and enabled understanding of the role of protein limitations in bio-based chemical production. Paper VI presents a hybrid model integrating an enzyme-constrained metabolic network, coupled to a gene regulatory model of nutrient-sensing mechanisms in S. cerevisiae. This model improves prediction of protein expression patterns while providing a rational connection between metabolism and the use of nutrients from the environment.This thesis demonstrates that integration of multiple systems biology approaches is valuable for understanding the connection of cell physiology at different levels, and provides tools for directed engineering of cells for the benefit of society

Path2Models: large-scale generation of computational models from biochemical pathway maps

Background: Systems biology projects and omics technologies have led to a growing number of biochemical pathway models and reconstructions. However, the majority of these models are still created de novo, based on literature mining and the manual processing of pathway data. Results: To increase the efficiency of model creation, the Path2Models project has automatically generated mathematical models from pathway representations using a suite of freely available software. Data sources include KEGG, BioCarta, MetaCyc and SABIO-RK. Depending on the source data, three types of models are provided: kinetic, logical and constraint-based. Models from over 2 600 organisms are encoded consistently in SBML, and are made freely available through BioModels Database at http://www.ebi.ac.uk/biomodels-main/path2models. Each model contains the list of participants, their interactions, the relevant mathematical constructs, and initial parameter values. Most models are also available as easy-to-understand graphical SBGN maps. Conclusions: To date, the project has resulted in more than 140 000 freely available models. Such a resource can tremendously accelerate the development of mathematical models by providing initial starting models for simulation and analysis, which can be subsequently curated and further parameterized

Modelling and multiobjective optimization for simulation of cyanobacterial metabolism

The present thesis is devoted to the development of models and algorithms to improve metabolic simulations of cyanobacterial metabolism. Cyanobacteria are photosynthetic bacteria of great biotechnological interest to the development of sustainable bio-based manufacturing processes. For this purpose, it is fundamental to understand metabolic behaviour of these organisms, and constraint-based metabolic modelling techniques offer a platform for analysis and assessment of cell's metabolic functionality. Reliable simulations are needed to enhance the applicability of the results, and this is the main goal of this thesis. This dissertation has been structured in three parts. The first part is devoted to introduce needed fundamentals of the disciplines that are combined in this work: metabolic modelling, cyanobacterial metabolism and multi-objective optimisation. In the second part the reconstruction and update of metabolic models of two cyanobacterial strains is addressed. These models are then used to perform metabolic simulations with the application of the classic Flux Balance Analysis (FBA) methodology. The studies conducted in this part are useful to illustrate the uses and applications of metabolic simulations for the analysis of living organisms. And at the same time they serve to identify important limitations of classic simulation techniques based on mono-objective linear optimisation that motivate the search of new strategies. Finally, in the third part a novel approach is defined based on the application of multi-objective optimisation procedures to metabolic modelling. Main steps in the definition of multi-objective problem and the description of an optimisation algorithm that ensure the applicability of the obtained results, as well as the multi-criteria analysis of the solutions are covered. The resulting tool allows the definition of non-linear objective functions and constraints, as well as the analysis of multiple Pareto-optimal solutions. It avoids some of the main drawbacks of classic methodologies, leading to more flexible simulations and more realistic results. Overall this thesis contributes to the advance in the study of cyanobacterial metabolism by means of definition of models and strategies that improve plasticity and predictive capacities of metabolic simulations.La presente tesis está dedicada al desarrollo de modelos y algoritmos para mejorar las simulaciones metabólicas de cianobacterias. Las cianobacterias son bacterias fotosintéticas de gran interés biotecnológico para el desarrollo de bioprocesos productivos sostenibles. Para este propósito, es fundamental entender el comportamiento metabólico de estos organismos, y el modelado metabólico basado en restricciones ofrece una plataforma para el análisis y la evaluación de las funcionalidades metabólicas de las células. Se necesitan simulaciones fidedignas para aumentar la aplicabilidad de los resultados, y este es el objetivo principal de esta tesis. Esta disertación se ha estructurado en tres partes. La primera parte está dedicada a introducir los fundamentos necesarios de las disciplinas que se combinan en este trabajo: el modelado metabólico, el metabolismo de cianobacterias, y la optimización multiobjetivo. En la segunda parte, se encara la reconstrucción y la actualización de los modelos metabólicos de dos cepas de cianobacterias. Estos modelos se usan después para llevar a cabo simulaciones metabólicas con la aplicación de la metodología clásica Flux Balance Analysis (FBA). Los estudios realizados en esta parte son útiles para ilustrar los usos y aplicaciones de las simulaciones metabólicas para el análisis de los organismos vivos. Y al mismo tiempo sirven para identificar importantes limitaciones de las técnicas clásicas de simulación basadas en optimización lineal mono-objetivo que motivan la búsqueda de nuevas estrategias. Finalmente, en la tercera parte, se define una nueva aproximación basada en la aplicación al modelado metabólico de procedimientos de optimización multiobjetivo. Se cubren los principales pasos en la definición de un problema multiobjetivo y la descripción de un algoritmo de optimización que aseguren la aplicabilidad de los resultados obtenidos, así como el análisis multi-criterio de las soluciones. La herramienta resultante permite la definición de funciones objetivo y restricciones no lineales, así como el análisis de múltiples soluciones en el sentido de Pareto. Esta herramienta evita algunos de los principales inconvenientes de las metodologías clásicas, lo que lleva a obtener simulaciones más flexibles y resultados más realistas. En conjunto, esta tesis contribuye al avance en el estudio del metabolismo de cianobacterias por medio de la definición de modelos y estrategias que mejoran la plasticidad y las capacidades predictivas de las simulaciones metabólicas.La present tesi està dedicada al desenvolupament de models i algorismes per a millorar les simulacions metabòliques de cianobacteris. Els cianobacteris són bacteris fotosintètics de gran interés biotecnològic per al desenvolupament de bioprocessos productius sostenibles. Per a aquest propòsit, és fonamental entendre el comportament metabòlic d'aquests organismes, i el modelatge metabòlic basat en restriccions ofereix una plataforma per a l'anàlisi i l'avaluació de les funcionalitats metabòliques de les cèl·lules. Es necessiten simulacions fidedignes per a augmentar l'aplicabilitat dels resultats, i aquest és l'objectiu principal d'aquesta tesi. Aquesta dissertació s'ha estructurat en tres parts. La primera part està dedicada a introduir els fonaments necessaris de les disciplines que es combinen en aquest treball: el modelatge metabòlic, el metabolisme de cianobacteris i l'optimització multiobjectiu. En la segona part, s'adreça la reconstrucció i l'actualització dels models metabòlics de dos soques de cianobacteris. Aquests models s'empren després per a portar a terme simulacions metabòliques amb l'aplicació de la metodologia clàssica Flux Balance Analysis (FBA). Els estudis realitzats en aquesta part són útils per a il·lustrar els usos i aplicacions de les simulacions metabòliques per a l'anàlisi dels organismes vius. I al mateix temps serveixen per a identificar importants limitacions de les tècniques clàssiques de simulació basades en optimització lineal mono-objectiu que motiven la cerca de noves estratègies. Finalment, en la tercera part, es defineix una nova aproximació basada en l'aplicació al modelatge metabòlic de procediments d'optimització multiobjectiu. Es cobreixen els principals passos en la definició d'un problema multiobjectiu i la descripció d'un algorisme d'optimització que asseguren l'aplicabilitat dels resultats obtinguts, així com l'anàlisi multi-criteri de les solucions. La ferramenta resultant permet la definició de funcions objectiu i restriccions no lineals, així com l'anàlisi de múltiples solucions òptimes en el sentit de Pareto. Aquesta ferramenta evita alguns dels principals inconvenients de les metodologies clàssiques, el que porta a obtenir simulacions més flexibles i resultats més realistes. En conjunt, aquesta tesi contribueix a l'avanç en l'estudi del metabolisme de cianobacteris per mitjà de la definició de models i estratègies que milloren la plasticitat i les capacitats predictives de les simulacions metabòliques.Siurana Paula, M. (2017). Modelling and multiobjective optimization for simulation of cyanobacterial metabolism [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/9057

Computational strategies for a system-level understanding of metabolism

Cell metabolism is the biochemical machinery that provides energy and building blocks to sustain life. Understanding its fine regulation is of pivotal relevance in several fields, from metabolic engineering applications to the treatment of metabolic disorders and cancer. Sophisticated computational approaches are needed to unravel the complexity of metabolism. To this aim, a plethora of methods have been developed, yet it is generally hard to identify which computational strategy is most suited for the investigation of a specific aspect of metabolism. This review provides an up-to-date description of the computational methods available for the analysis of metabolic pathways, discussing their main advantages and drawbacks. In particular, attention is devoted to the identification of the appropriate scale and level of accuracy in the reconstruction of metabolic networks, and to the inference of model structure and parameters, especially when dealing with a shortage of experimental measurements. The choice of the proper computational methods to derive in silico data is then addressed, including topological analyses, constraint-based modeling and simulation of the system dynamics. A description of some computational approaches to gain new biological knowledge or to formulate hypotheses is finally provided