Dynamic optimization of metabolic networks coupled with gene expression

The regulation of metabolic activity by tuning enzyme expression levels is crucial to sustain cellular growth in changing environments. Metabolic networks are often studied at steady state using constraint-based models and optimization techniques. However, metabolic adaptations driven by changes in gene expression cannot be analyzed by steady state models, as these do not account for temporal changes in biomass composition. Here we present a dynamic optimization framework that integrates the metabolic network with the dynamics of biomass production and composition, explicitly taking into account enzyme production costs and enzymatic capacity. In contrast to the established dynamic flux balance analysis, our approach allows predicting dynamic changes in both the metabolic fluxes and the biomass composition during metabolic adaptations. We applied our algorithm in two case studies: a minimal nutrient uptake network, and an abstraction of core metabolic processes in bacteria. In the minimal model, we show that the optimized uptake rates reproduce the empirical Monod growth for bacterial cultures. For the network of core metabolic processes, the dynamic optimization algorithm predicted commonly observed metabolic adaptations, such as a diauxic switch with a preference ranking for different nutrients, re-utilization of waste products after depletion of the original substrate, and metabolic adaptation to an impending nutrient depletion. These examples illustrate how dynamic adaptations of enzyme expression can be predicted solely from an optimization principle

Machine learning applied to enzyme turnover numbers reveals protein structural correlates and improves metabolic models.

Knowing the catalytic turnover numbers of enzymes is essential for understanding the growth rate, proteome composition, and physiology of organisms, but experimental data on enzyme turnover numbers is sparse and noisy. Here, we demonstrate that machine learning can successfully predict catalytic turnover numbers in Escherichia coli based on integrated data on enzyme biochemistry, protein structure, and network context. We identify a diverse set of features that are consistently predictive for both in vivo and in vitro enzyme turnover rates, revealing novel protein structural correlates of catalytic turnover. We use our predictions to parameterize two mechanistic genome-scale modelling frameworks for proteome-limited metabolism, leading to significantly higher accuracy in the prediction of quantitative proteome data than previous approaches. The presented machine learning models thus provide a valuable tool for understanding metabolism and the proteome at the genome scale, and elucidate structural, biochemical, and network properties that underlie enzyme kinetics

Decoding Complexity in Metabolic Networks using Integrated Mechanistic and Machine Learning Approaches

How can we get living cells to do what we want? What do they actually ‘want’? What ‘rules’ do they observe? How can we better understand and manipulate them? Answers to fundamental research questions like these are critical to overcoming bottlenecks in metabolic engineering and optimizing heterologous pathways for synthetic biology applications. Unfortunately, biological systems are too complex to be completely described by physicochemical modeling alone. In this research, I developed and applied integrated mechanistic and data-driven frameworks to help uncover the mysteries of cellular regulation and control. These tools provide a computational framework for seeking answers to pertinent biological questions. Four major tasks were accomplished. First, I developed innovative tools for key areas in the genome-to-phenome mapping pipeline. An efficient gap filling algorithm (called BoostGAPFILL) that integrates mechanistic and machine learning techniques was developed for the refinement of genome-scale metabolic network reconstructions. Genome-scale metabolic network reconstructions are finding ever increasing applications in metabolic engineering for industrial, medical and environmental purposes. Second, I designed a thermodynamics-based framework (called REMEP) for mutant phenotype prediction (integrating metabolomics, fluxomics and thermodynamics data). These tools will go a long way in improving the fidelity of model predictions of microbial cell factories. Third, I designed a data-driven framework for characterizing and predicting the effectiveness of metabolic engineering strategies. This involved building a knowledgebase of historical microbial cell factory performance from published literature. Advanced machine learning concepts, such as ensemble learning and data augmentation, were employed in combination with standard mechanistic models to develop a predictive platform for important industrial biotechnology metrics such as yield, titer, and productivity. Fourth, my modeling tools and skills have been used for case studies on fungal lipid metabolism analyses, E. coli resource allocation balances, reconstruction of the genome-scale metabolic network for a non-model species, R. opacus, as well as the rapid prediction of bacterial heterotrophic fluxomics. In the long run, this integrated modeling approach will significantly shorten the “design-build-test-learn” cycle of metabolic engineering, as well as provide a platform for biological discovery

Multi-omics integration accurately predicts cellular state in unexplored conditions for Escherichia coli.

A significant obstacle in training predictive cell models is the lack of integrated data sources. We develop semi-supervised normalization pipelines and perform experimental characterization (growth, transcriptional, proteome) to create Ecomics, a consistent, quality-controlled multi-omics compendium for Escherichia coli with cohesive meta-data information. We then use this resource to train a multi-scale model that integrates four omics layers to predict genome-wide concentrations and growth dynamics. The genetic and environmental ontology reconstructed from the omics data is substantially different and complementary to the genetic and chemical ontologies. The integration of different layers confers an incremental increase in the prediction performance, as does the information about the known gene regulatory and protein-protein interactions. The predictive performance of the model ranges from 0.54 to 0.87 for the various omics layers, which far exceeds various baselines. This work provides an integrative framework of omics-driven predictive modelling that is broadly applicable to guide biological discovery

Essential plasticity and redundancy of metabolism unveiled by synthetic lethality analysis

We unravel how functional plasticity and redundancy are essential mechanisms underlying the ability to survive of metabolic networks. We perform an exhaustive computational screening of synthetic lethal reaction pairs in Escherichia coli in a minimal medium and we find that synthetic lethal pairs divide in two different groups depending on whether the synthetic lethal interaction works as a backup or as a parallel use mechanism, the first corresponding to essential plasticity and the second to essential redundancy. In E. coli, the analysis of pathways entanglement through essential redundancy supports the view that synthetic lethality affects preferentially a single function or pathway. In contrast, essential plasticity, the dominant class, tends to be inter-pathway but strongly localized and unveils Cell Envelope Biosynthesis as an essential backup for Membrane Lipid Metabolism. When comparing E. coli and Mycoplasma pneumoniae, we find that the metabolic networks of the two organisms exhibit a large difference in the relative importance of plasticity and redundancy which is consistent with the conjecture that plasticity is a sophisticated mechanism that requires a complex organization. Finally, coessential reaction pairs are explored in different environmental conditions to uncover the interplay between the two mechanisms. We find that synthetic lethal interactions and their classification in plasticity and redundancy are basically insensitive to medium composition, and are highly conserved even when the environment is enriched with nonessential compounds or overconstrained to decrease maximum biomass formation.Comment: 22 pages, 4 figure

Coupling metabolic footprinting and flux balance analysis to predict how single gene knockouts perturb microbial metabolism

Tese de mestrado. Biologia (Bioinformática e Biologia Computacional). Universidade de Lisboa, Faculdade de Ciências, 2012The model organisms Caenorhabditis elegans and E. coli form one of the simplest gut microbe host interaction models. Interventions in the microbe that increase the host longevity including inhibition of folate synthesis have been reported previously. To find novel single gene knockouts with an effect on lifespan, a screen of the Keio collection of E. coli was undertaken, and some of the genes found are directly involved in metabolism. The next step in those specific cases is to understand how these mutations perturb metabolism systematically, so that hypotheses can be generated. For that, I employed dynamic Flux Balance Analysis (dFBA), a constraint-based modeling technique capable of simulating the dynamics of metabolism in a batch culture and making predictions about changes in intracellular flux distribution. Since the specificities of the C. elegans lifespan experiments demand us to culture microbes in conditions differing from most of the published literature on E. coli physiology, novel data must be acquired to characterize and make dFBA simulations as realistic as possible. To do this exchange fluxes were measured using quantitative H NMR Time-Resolved Metabolic Footprinting. Furthermore, I also investigate the combination of TReF and dFBA as a tool in microbial metabolism studies. These approaches were tested by comparing wild type E. coli with one of the knockout strains found, ΔmetL, a knockout of the metL gene which encodes a byfunctional enzyme involved in aspartate and threonine metabolism. I found that the strain exhibits a slower growth rate than the wild type. Model simulation results revealed that reduced homoserine and methionine synthesis, as well as impaired sulfur and folate metabolism are the main effects of this knockout and the reasons for the growth deficiency. These results indicate that there are common mechanisms of the lifespan extension between ΔmetL and inhibition of folate biosynthesis and that the flux balance analysis/metabolic footprinting approach can help us understand the nature of these mechanisms.Os organismos modelo Caenorhabditis elegans e E. coli formam um dos modelos mais simples de interacções entre micróbio do tracto digestivo e hospedeiro. Intervenções no micróbio capazes de aumentar a longevidade do hospedeiro, incluindo inibição de síntese de folatos, foram reportadas previamente. Para encontrar novas delecções génicas do micróbio capazes de aumentar a longevidade do hospedeiro, a colecção Keio de deleções génicas de E. coli foi rastreada. Alguns dos genes encontrados participam em processos metabólicos, e nesses casos, o próximpo passo é perceber como as deleções perturbam o metabolismo sistémicamente, para gerar hipóteses. Para isso, utilizo dynamic Flux Balance Analysis (dFBA), uma técnica de modelação metabólica capaz de fazer previsões sobre alterações na distribuição intracelular de fluxos. As especificidades das experiências de tempo de vida em C.elegans obrigam-nos a trabalhar em condições diferentes das usadas na maioria da literatura publicada em fisiologia de E. coli, e para dar o máximo realismo às simulações de dFBA novos dados foram adquiridos, utilizando H NMR Time-Resolved Metabolic Footprinting para medir fluxos de troca de metabolitos entre microorganismo e meio de cultura. A combinação de TReF e dFBA como ferramenta de estudo do metabolism microbiano é também investigada. Estas abordagens foram testadas ao comparar E. coli wild-type com uma das estirpes encontradas no rastreio, ΔmetL, knockout do gene metL, que codifica um enzima bifunctional participante no metabolismo de aspartato e treonina, e que exibe uma taxa de crescimento reduzida comparativamente ao wild-type. Os resultados das simulações revelaram que os principais efeitos da deleção deste gene, e as razões para a menor taxa de crescimento observada, são a produção reduzida de homoserina e metionina e os efeitos que provoca no metabolismo de folatos e enxofre. Estes resultados indicam que há mecanismos comuns na extensão da longevidade causada por esta deleção e inibição de síntese de folatos, e que a combinação metabolic footprinting/flux balance analysis pode ajudar-nos a compreender a natureza desses mecanismos

solveME: fast and reliable solution of nonlinear ME models.

BackgroundGenome-scale models of metabolism and macromolecular expression (ME) significantly expand the scope and predictive capabilities of constraint-based modeling. ME models present considerable computational challenges: they are much (>30 times) larger than corresponding metabolic reconstructions (M models), are multiscale, and growth maximization is a nonlinear programming (NLP) problem, mainly due to macromolecule dilution constraints.ResultsHere, we address these computational challenges. We develop a fast and numerically reliable solution method for growth maximization in ME models using a quad-precision NLP solver (Quad MINOS). Our method was up to 45 % faster than binary search for six significant digits in growth rate. We also develop a fast, quad-precision flux variability analysis that is accelerated (up to 60× speedup) via solver warm-starts. Finally, we employ the tools developed to investigate growth-coupled succinate overproduction, accounting for proteome constraints.ConclusionsJust as genome-scale metabolic reconstructions have become an invaluable tool for computational and systems biologists, we anticipate that these fast and numerically reliable ME solution methods will accelerate the wide-spread adoption of ME models for researchers in these fields

Network Analysis and Modeling in Systems Biology

This thesis is dedicated to the study and comprehension of biological networks at the molecular level. The objectives were to analyse their topology, integrate it in a genotype-phenotype analysis, develop richer mathematical descriptions for them, study their community structure and compare different methodologies for estimating their internal fluxes. The work presented in this document moves around three main axes. The first one is the biological. Which organisms were studied in this thesis? They range from the simplest biological agents, the viruses, in this case the Potyvirus genus to prokariotes such as Escherichia coli and complex eukariotes (Arabidopsis thaliana, Nicotiana benthamiana). The second axis refers to which biological networks were studied. Those are protein-protein interaction (PPIN) and metabolic networks (MN). The final axis relates to the mathematical and modelling tools used to generate knowledge from those networks. These tools can be classify in three main branches: graph theory, constraint-based modelling and multivariate statistics. The document is structured in six parts. The first part states the justification for the thesis, exposes a general thesis roadmap and enumerates its main contributions. In the second part important literature is reviewed, summarized and integrated. From the birth and development of Systems Biology to one of its most popular branches: biological network analysis. Particular focus is put on PPIN and MN and their structure, representations and features. Finally a general overview of the mathematical tools used is presented. The third, fourth and fifth parts represent the central work of this thesis. They deal respectively with genotypephenotype interaction and classical network analysis, constraint-based modelling methods comparison and modelling metabolic networks and community structure. Finally, in the sixth part the main conclusions of the thesis are summarized and enumerated. This thesis highlights the vital importance of studying biological entities as systems and how powerful and promising this integrated analysis is. Particularly, network analysis becomes a fundamental avenue of research to gain insight into those biological systems and to extract, integrate and display this new information. It generates knowledge from just data.Esta tesis está dedicada al estudio y comprensión de redes biológicas a nivel molecular. Los objetivos fueron analizar su topología, integrar esta en un análisis de genotipo-fenotipo, desarrollar descripciones matemáticas más completas para ellas, estudiar su estructura de comunidades y comparar diferentes metodologías para estimar sus flujos internos. El trabajo presentado en este documento gira entorno a tres ejes principales. El primero es el biológico. ¿Qué organismos han sido estudiados en esta tesis? Estos van desde los agentes biológicos mas simples, los virus, en este caso el género Potyvirus, hasta procariotas como Escherichia coli y eucariotas complejos (Arabidopsis thaliana, Nicotiana benthamiana). El segundo eje hace referencia a las redes biológicas estudiadas, que fueron redes de interacción de proteínas (PPIN) y redes metabólicas (MN). El eje final es el de las herramientas matemáticas y de modelización empleadas para interrogar esas redes. Estas herramientas pueden clasificarse en tres grandes grupos: teoría de grafos, modelización basada en restricciones y estadística multivariante. Este documento está estructurado en seis partes. La primera expone la justificación para la tesis, muestra un mapa visual de la misma y enumera sus contribuciones principales. En la segunda parte, la bibliografía relevante es revisada y resumida. Desde el nacimiento y desarrollo de la Biología de Sistemas hasta una de sus ramas más populares: el análisis de redes biomoleculares. Especial interés es puesto en PPIN y MN: su estructura, representación y características. Finalmente, un resumen general de las herramientas matemáticas usadas es presentado. Los capítulos tercero, cuarto y quinto representan el cuerpo central de esta tesis. Estos tratan respectivamente sobre la interacción de genotipo-fenotipo y análisis topolólogico clásico de redes, modelos basados en restricciones y modelización de redes metabólicas y su estructura de comunidades. Finalmente, en la sexta parte las principales conclusiones de la tesis son resumidas y expuestas. Esta tesis pone énfasis en la vital importancia de estudiar los fenómenos biológicos como sistemas y en la potencia y prometedor futuro de este análisis integrativo. En concreto el análisis de redes supone un camino de investigación fundamental para obtener conocimiento sobre estos sistemas biológicos y para extraer y mostrar información sobre los mismos. Este análisis genera conocimiento partiendo únicamente desde datos.Aquesta tesi està dedicada a l'estudi i comprensió de xarxes biològiques a nivell molecular. Els objectius van ser analitzar la seva topologia, integrar aquesta en una anàlisi de genotip-fenotip, desenvolupar descripcions matemàtiques més completes per a elles, estudiar la seva estructura de comunitats o modularitat i comparar diferents metodologies per estimar els fluxos interns. El treball presentat en aquest document gira entorn de tres eixos principals. El primer és el biològic. ¿Què organismes han estat estudiats en aquesta tesi? Aquests van des dels agents biològics mes simples, els virus, en aquest cas el gènere Potyvirus, fins procariotes com Escherichia coli i eucariotes complexos (Arabidopsis thaliana, Nicotiana benthamiana). El segon eix fa referència a les xarxes biològiques estudiades, que van ser les xarxes d'interacció de proteïnes (PPIN) i les xarxes metabòliques (MN). L'eix final és el de les eines matemàtiques i de modelització emprades per interrogar aquestes xarxes. Aquestes eines poden classificarse en tres grans grups: teoria de grafs, modelització basada en restriccions i estadística multivariant. Aquest document està estructurat en sis parts. La primera exposa la justificació per a la tesi, mostra un mapa visual de la mateixa i enumera les seves contribucions principals. A la segona part, la bibliografia rellevant és revisada i resumida. Des del naixement i desenvolupament de la Biologia de Sistemes fins a una de les seves branques més populars: l'anàlisi de xarxes moleculars. Especial interès és posat en PPIN i MN: la seva estructura, representació i característiques. Finalment, un resum general de les eines matemàtiques utilitzades és presentat. Els capítols tercer, quart i cinquè representen el cos central d'aquesta tesi. Aquests tracten respectivament sobre la interacció de genotip-fenotip i anàlisi topolólogico clàssic de xarxes, models basats en restriccions i modelització de xarxes metabòliques i la seva estructura de comunitats. Finalment, en la sisena part les principals conclusions de la tesi són resumides i exposades. Aquesta tesi posa èmfasi en la vital importància d'estudiar els fenòmens biològics com sistemes i en la potència i prometedor futur d'aquesta anàlisi integratiu. En concret l'anàlisi de xarxes suposa un camí d'investigació fonamental per obtenir coneixement sobre aquests sistemes biològics i per extreure i mostrar informació sobre els mateixos. Aquest anàlisi genera coneixement partint únicament des de dades.Bosque Chacón, G. (2017). Network Analysis and Modeling in Systems Biology [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/79082TESI