Search NASA⌕ Search

SEARCH · Search NASA

Results for “genome-scale”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

iNovo479: Metabolic Modeling Provides a Roadmap to Optimize Bioproduct Yield from Deconstructed Lignin Aromatics by Novosphingobium aromaticivorans

Lignin is an abundant renewable source of aromatics and precursors for the production of other organic chemicals. However, lignin is a heterogeneous polymer, so the mixture of aromatics released during its depolymerization can make its conversion to chemicals challenging. Microbes are a potential solution to this challenge, as some can catabolize multiple aromatic substrates into one product. Novosphingobium aromaticivorans has this ability, and its use as a bacterial chassis for lignin valorization could be improved by the ability to predict product yields based on thermodynamic and metabolic inputs. In this work, we built a genome-scale metabolic model of N. aromaticivorans, iNovo479, to guide the engineering of strains for aromatic conversion into products. iNovo479 predicted product yields from single or multiple aromatics, and the impact of combinations of aromatic and non-aromatic substrates on product yields. We show that enzyme reactions from other organisms can be added to iNovo479 to predict the feasibility and profitability of producing additional products by engineered strains. Thus, we conclude that iNovo479 can help guide the design of bacteria to convert lignin aromatics into valuable chemicals.

59 BASIC BIOLOGICAL SCIENCES↗

Uses of Multi-Objective Flux Analysis for Optimization of Microbial Production of Secondary Metabolites

Secondary metabolites are not essential for the growth of microorganisms, but they play a critical role in how microbes interact with their surroundings. In addition to this important ecological role, secondary metabolites also have a variety of agricultural, medicinal, and industrial uses, and thus the examination of secondary metabolism of plants and microbes is a growing scientific field. While the chemical production of certain secondary metabolites is possible, industrial-scale microbial production is a green and economically attractive alternative. This is even more true, given the advances in bioengineering that allow us to alter the workings of microbes in order to increase their production of compounds of interest. This type of engineering requires detailed knowledge of the “chassis” organism’s metabolism. Since the resources and the catalytic capacity of enzymes in microbes is finite, it is important to examine the tradeoffs between various bioprocesses in an engineered system and alter its working in a manner that minimally perturbs the robustness of the system while allowing for the maximum production of a product of interest. The in silico multi-objective analysis of metabolism using genome-scale models is an ideal method for such examinations.

59 BASIC BIOLOGICAL SCIENCES↗

CompLaB v1.0: a scalable pore-scale model for flow, biogeochemistry, microbial metabolism, and biofilm dynamics

Abstract. Microbial activity and chemical reactions in porous media depend on the local conditions at the pore scale and can involve complex feedback with fluid flow and mass transport. We present a modeling framework that quantitatively accounts for the interactions between the bio(geo)chemical and physical processes and that can integrate genome-scale microbial metabolic information into a dynamically changing, spatially explicit representation of environmental conditions. The model couples a lattice Boltzmann implementation of Navier–Stokes (flow) and advection–diffusion-reaction (mass conservation) equations. Reaction formulations can include both kinetic rate expressions and flux balance analysis, thereby integrating reactive transport modeling and systems biology. We also show that the use of surrogate models such as neural network representations of in silico cell models can speed up computations significantly, facilitating applications to complex environmental systems. Parallelization enables simulations that resolve heterogeneity at multiple scales, and a cellular automaton module provides additional capabilities to simulate biofilm dynamics. The code thus constitutes a platform suitable for a range of environmental, engineering and – potentially – medical applications, in particular ones that involve the simulation of microbial dynamics.

58 GEOSCIENCES↗

Conjugation-based genome engineering enables rapid prototyping and bioproduction in non-model bacteria

Abstract Non-model bacteria offer unique metabolic capabilities for sustainable bioproduction, yet their limited genetic accessibility hinders systematic strain development. Here we present conjugation-based serine recombinase-assisted genome engineering (cSAGE), a broad-host-range platform that enables predictable, iterative genomic integration in transformation-resistant bacteria. cSAGE combines conjugative DNA delivery, standardized low-copy vectors, orthogonal recombinases, and modular genetic parts to support rapid pathway assembly and cross-host benchmarking. Using purple nonsulfur bacteria as a testbed, we integrate promoter engineering, multi-payload genome modification, and genome-scale metabolic modeling to empirically evaluate host-dependent pathway performance. Applying this workflow, we identify strain-specific differences in photosynthetic conversion of lignin-derived p -coumarate to the thermoplastic precursor p -vinylphenol. By enabling genome engineering and functional comparison across diverse bacteria using a single plasmid system, cSAGE provides a general framework for non-model strain prototyping and biotransformation discovery.

Guzman, Michael S. [Department of Chemical Enginee↗

Metabolic model-based ecological modeling for probiotic design

The microbial community composition in the human gut has a profound effect on human health. This observation has lead to extensive use of microbiome therapies, including over-the-counter ‘probiotic’ treatments intended to alter the composition of the microbiome. Despite so much promise and commercial interest, the factors that contribute to the success or failure of microbiome-targeted treatments remain unclear. We investigate the biotic interactions that lead to successful engraftment of a novel bacterial strain introduced to the microbiome as in probiotic treatments. We use pairwise genome-scale metabolic modeling with a generalized resource allocation constraint to build a network of interactions between taxa that appear in an experimental engraftment study. We create induced sub-graphs using the taxa present in individual samples and assess the likelihood of invader engraftment based on network structure. To do so, we use a generalized Lotka-Volterra model, which we show has strong ability to predict if a particular invader or probiotic will successfully engraft into an individual’s microbiome. Furthermore, we show that the mechanistic nature of the model is useful for revealing which microbe-microbe interactions potentially drive engraftment.

59 BASIC BIOLOGICAL SCIENCES↗

The Genetics of Aerotolerant Growth in Zymomonas mobilis

Reduced genome bacteria are genetically simplified systems that facilitate biological study and industrial use. The free-living Alphaproteobacterium, Zymomonas mobilis, has a naturally reduced genome containing fewer than 2000 protein coding genes. Despite its small genome, Z. mobilis thrives in diverse conditions including the presence or absence of atmospheric oxygen. However, insufficient characterization of essential and conditionally essential genes has limited broader adoption of Z. mobilis as a model Alphaproteobacterium. Here, we use genome-scale CRISPRi-seq to systematically identify and characterize Z. mobilis genes that are conditionally essential for aerotolerant or anaerobic growth, or are generally essential across both conditions. Comparative genomics revealed that the essentiality of most "generally essential" genes was shared between Z. mobilis and other Alphaproteobacteria, validating Z. mobilis as reduced genome model. Among conditionally essential genes, we found that the DNA repair gene, recJ, was critical only for aerobic growth but reduced the mutation rate under both conditions. Further, we show that genes encoding the F1FO ATP synthase and Rnf respiratory complex are required for anaerobic growth of Z. mobilis. Combining CRISPRi partial knockdowns with metabolomics and membrane potential measurements, we determined that the ATP synthase generates membrane potential that is consumed by Rnf to power downstream processes. Rnf knockdown strains accumulated isoprenoid biosynthesis intermediates, suggesting a key role for Rnf in powering essential biosynthetic reactions. Our work establishes Z. mobilis as a streamlined model for alphaproteobacterial genetics, has broad implications in bacterial energy coupling, and informs Z. mobilis genome manipulation for optimized production of valuable isoprenoid-based bioproducts

Alphaproteobacteria↗

Zymomonas mobilis oxidative stress transcriptomics

Zymomonas mobilis is an important bioenergy organism that has potential to produce biofuels, including ethanol, in high volumes. Here we examined the response of Zymomonas mobilis to various oxidative stresses using genome-scale transcriptomics data. We first examined the transcrpit abundance in WT aerobic growth compared to aerobic grown in paraquat, which forms superoxide. Under anaerobic growth conditions we compared WT Zymomonas mobilis with strains grown in media lacking iron as well as strains lacking iron that were treated with the iron chelator DIP before collection. Finally we examined transcript abundance in cells lacking ZMO_0422 (Rrf2 family transcription factor homolog) and ZMO_1411 (Fur homolog) grown under anaerobic conditions. Overall design: Transcriptomic analysis of WT, a deletion of ZMO_0422, and a deletion of ZMO_1411 in Zymomonas mobilis ZM4 under aerobic and anaerobic growth conditions along with various oxidative stresses: Paraquate addition, No Iron, and hydrogen peroxide addition.

aerobic↗

Modeling the Activity of Single Genes

The central dogma of molecular biology states that information is stored in DNA, transcribed to messenger RNA (mRNA) and then translated into proteins. This picture is significantly augmentated when we consider the action of certain proteins in regulating transcription. These transcription factors provide a feedback pathway by which genes can regulate one another's expression as mRNA and then as protein. To review: DNA, RNA and proteins have different functions. DNA is the molecular storehouse of genetic information. When cells divide, the DNA is replicated, so that each daughter cell maintains the same genetic information as the mother cell. RNA acts as a go-between from DNA to proteins. Only a single copy of DNA is present, but multiple copies of the same piece of RNA may be present, allowing cells to make huge amounts of protein. In eukaryotes (organisms with a nucleus), DNA is found in the nucleus only. RNA is copied in the nucleus then translocates(moves) outside the nucleus, where it is transcribed into proteins. Along the way, the RNA may be spliced, i.e., may have pieces cut out. RNA then attaches to ribosomes and is translated to proteins. Proteins are the machinery of the cell other than DNA and RNA, all the complex molecules of the cell are proteins. Proteins are specialized machines, each of which fulfills its own task, which may be transporting oxygen, catalyzing reactions, or responding to extracellular signals, just to name a few. One of the more interesting functions a protein may have is binding directly or indirectly to DNA to perform transcriptional regulation, thus forming a closed feedback loop of gene regulation. The structure of DNA and the central dogma were understood in the 50s; in the early 80s it became possible to make arbitrary modifications to DNA and use cellular machinery to transcribe and translate the resulting genes; more recently, genomes (i.e., the complete DNA sequence) of many organisms have been sequenced. This large-scale sequencing began with simple organisms, viruses and bacteria, progressed to eukaryotes such as yeast, and more recently (1998) progressed to a multi-cellular animal, the nematode Caenorhabditis elegans. Sequencers have now moved on to the fruit fly Drosophila melanogaster, whose sequence is slated for completion by the end of 1999. The human genome project is expected to determine the complete sequence of all 3 billion bases of human DNA within the next five years. In the wake of genome-scale sequencing, further instrumentation is being developed to assay gene expression and function on a comparably large scale. Much of the work in computational biology focuses on computational tools used in sequencing, finding genes that are related to a particular gene, finding which parts of the DNA code for proteins and which do not, understanding what proteins will be formed from a given length of DNA, predicting how the proteins will fold from a one-dimensional structure into a three dimensional structure, and so on. Much less computational work has been done regarding the function of proteins. One reason for this is that different proteins function very differently, and so work on protein function is very specific to certain classes of proteins. There are, for example, proteins such enzymes that catalyze various intracellular reactions, receptors that respond to extracellular signals and ion channels that regulate the flow of charged particles into and out of the cell. In this chapter, we will consider a particular class of proteins called transcription factors(TFs), which are responsible for regulating when a certain gene is expressed in a certain cell, which cells it is express in, and how much is expressed. Understanding these processes will involve developing a deeper understanding of transcription, translation, and the cellular processes that control those processes. All of these elements fall under the aegis of gene regulation or more narrowly transcriptional regulation. Some of the key questions in gene regulation are: What genes are expressed in a certain cell at a certain time? How does gene expression differ from cell to cell in a multicellular organism? Which proteins act as transcription factors, i.e., are important in regulating gene expression? From questions like these, we hope to understand which genes are important for various macroscopic processes. Nearly all of the cells of a multicellular organism contain the same DNA. Yet this same genetic information yields a large number of different cell types. The fundamental difference between a neuron and a liver cell, for example, is which genes are expressed. Thus understanding gene regulation is an important step in understanding development. Furthermore, understanding the usual genes that are expressed in cells may give important clues about various diseases. Some diseases, such as sickle cell anemia and cystic fibrosis, are caused by defects in single, non-regulatory genes; others, such as certain cancers, are caused when the cellular control circuitry malfunctions - an understanding of these diseases will involve pathways of multiple interacting gene products. There are numerous challenges in the area of understanding and modeling gene regulation. First and foremost, biologists would like to develop a deeper understanding of the processes involved, including which genes and families of genes are important, how they interact, etc. From a computation point of view, there has been embarrassingly little work done. In this chapter there are many areas in which we can phrase meaningful, non-trivial computational questions, but questions that have not been addressed. Some of these are purely computational (what is a good algorithm for dealing with a model of type X) and others are more mathematical (given a system with certain characteristics, what sort of model can one use? How does one find biochemical parameters from system-level behavior using as few experiments as possible?). In addition to biological and algorithmic problems, there is also the ever-present issue of theoretical biology - what general principles can be derived from these systems, what can one do with models other than just simulate time-courses, what can be deduced about a class of systems without knowing all the details? The fundamental challenge to computationalists and theorists is to add value to the biology - to use models, modeling techniques and algorithms to understand the biology in new ways.

Mjolsness, Eric↗

The Transcriptional Response of Diverse Saccharomyces Cerevisiae Strains to Simulated Microgravity

Spaceflight imposes multiple stresses on biological systems resulting in genome-scale adaptations. Understanding these adaptations and their underlying molecular mechanisms is important to clarifying and reducing the risks associated with spaceflight. One such risk is infection by microbes present in spacecraft and their associated systems and inhabitants. This risk is compounded by results suggesting that some microbes may exhibit increased virulence after exposure to spaceflight conditions. The yeast, S. cerevisiae, is a powerful microbial model system, and its response to spaceflight has been studied for decades. However, to date, these studies have utilized common lab strains. Yet studies on trait variation in S. cerevisiae demonstrate that these lab strains are not representative of wild yeast and instead respond to environmental stimuli in an atypical manner. Thus, it is not clear how transferable these results are to the wild S. cerevisiae strains likely to be encountered during spaceflight. To determine if diverse S. cerevisiae strains exhibit a conserved response to simulated microgravity, we will utilize a collection of 100 S. cerevisiae strains isolated from clinical, environmental and industrial settings. We will place selected S. cerevisiae strains in simulated microgravity using a high-aspect rotating vessel (HARV) and document their transcriptional response by RNA-sequencing and quantify similarities and differences between strains. Our research will have a strong impact on the understanding of how genetic diversity of microorganisms effects their response to spaceflight, and will serve as a platform for further studies.

yeast↗

The Transcriptional Response of Diverse Saccharomyces Cerevisiae Strains to Simulated Microgravity

Spaceflight imposes multiple stresses on biological systems resulting in genome-scale adaptations. Understanding these adaptations and their underlying molecular mechanisms is important to clarifying and reducing the risks associated with spaceflight. One such risk is infection by microbes present in spacecraft and their associated systems and inhabitants. This risk is compounded by results suggesting that some microbes may exhibit increased virulence after exposure to spaceflight conditions. The yeast, S. cerevisiae, is a powerful microbial model system, and it's response to spaceflight has been studied for decades. However, to date, these studies have utilized common lab strains. Yet studies on trait variation in S. cerevisiae demonstrate that these lab strains are not representative of wild yeast and instead respond to environmental stimuli in an atypical manner. Thus, it is not clear how transferable these results are to the wild S. cerevisiae strains likely to be encountered during spaceflight. To determine if diverse S. cerevisiae strains exhibit a conserved response to simulated microgravity, we will utilize a collection of 100 S. cerevisiae strains isolated from clinical, environmental and industrial settings. We will place selected S. cerevisiae strains in simulated microgravity using a high-aspect rotating vessel (HARV) and document their transcriptional response by RNA-sequencing and quantify similarities and differences between strains. Our research will have a strong impact on the understanding of how genetic diversity of microorganisms effects their response to spaceflight, and will serve as a platform for further studies.

genetic diversity↗

The Transcriptional Response of Diverse Saccharomyces cerevisiae Strains to Simulated Microgravity

Spaceflight imposes multiple stresses on biological systems resulting in genome-scale adaptations. Understanding these adaptations and their underlying molecular mechanisms is important to clarifying and reducing the risks associated with spaceflight. One such risk is infection by microbes present in spacecraft and their associated systems and inhabitants. This risk is compounded by results suggesting that some microbes may exhibit increased virulence after exposure to spaceflight conditions. The yeast, S. cerevisiae, is a powerful microbial model system, and it's response to spaceflight has been studied for decades. However, to date, these studies have utilized common lab strains. Yet studies on trait variation in S. cerevisiae demonstrate that these lab strains are not representative of wild yeast and instead respond to environmental stimuli in an a typical manner. Thus, it is not clear how transferable these results are to the wild S. cerevisiae strains likely to be encountered during spaceflight. To determine if diverse S. cerevisiae strains exhibit a conserved response to simulated microgravity, we will utilize a collection of 100 S. cerevisiae strains isolated from clinical, environmental and industrial settings. We will place selected S. cerevisiae strains in simulated microgravity using a high-aspect rotating vessel (HARV) and document their transcriptional response by RNA-sequencing and quantify similarities and differences between strains. Our research will have a strong impact on the understanding of how genetic diversity of microorganisms effects their response to spaceflight, and will serve as a platform for further studies.

spaceflight↗

Application of the Metabolic Modeling Pipeline in KBase to Categorize Reactions, Predict Essential Genes, and Predict Pathways in an Isolate Genome

The DOE Systems Biology Knowledgebase (KBase) platform offers a range of powerful tools for the reconstruction, refinement, and analysis of genome-scale metabolic models built from microbial isolate genomes. In this chapter, we describe and demonstrate these tools in action with an analysis of isoprene production in the Bacillus subtilis DSM genome. Two different methods are applied to build initial metabolic models for the DSM genome, then the models are gapfilled in three different growth conditions. Next, flux balance analysis (FBA) and flux variability analysis (FVA) techniques are applied to both study the growth of these models in minimal media and classify reactions within each model based on essentiality and functionality. The models are applied with the FBA method to predict essential genes, which are then compared to an updated list of essential genes obtained for B. subtilis 168, a very similar strain to the DSM isolate. The models are also applied to simulate Biolog growth conditions, and these results are compared with Biolog data collected for B. subtilis 168. Finally, the DSM metabolic models are applied to explore the pathways and genes responsible for producing isoprene in this strain. These studies demonstrate the accuracy and utility of models generated from the KBase pipelines, as well as exploring the tools available for analyzing these models.

Allen, Benjamin↗

IMAGINE BioSecurity: Mesocosm-Based Methods to Evaluate Biocontainment Strategies and Impact of Industrial Microbes Upon Native Ecosystems

Project Goals: The Integrative Modeling and Genome-scale Engineering for Biosystems Security (IMAGINE BioSecurity) SFA project seeks to establish an understanding of the behavior of engineered microbes in controlled versus environmental conditions to predictively devise new strategies for responding to biological escape. To this end, the IMAGINE Team has established a plant-soil mesocosm platform to track and quantify the fate of industrial microbes in environmental systems and assess the efficacy of biocontainment constraints upon genetically engineered microbe escape frequency and the impact of industrial microbes upon native ecological microbiomes. Abstract Text: Genetically modified industrial production microbes and their associated bioproducts have emerged as an integral component of a sustainable bioeconomy. However, the rapid development of these innovative technologies raises biosecurity concerns, namely, the risk of environmental escape. Thus, the realization of a bioeconomy hinges not only on the development and deployment of microbial production hosts, but also on the development of secure biosystems and biocontainment designs. Current laboratory-based biocontainment testing systems do not accurately reflect complexities found in natural environments, necessitating an environmentally relevant analysis pipeline that allows for the detection of rare escapees, the effect of associated bio-products, and the impact on native ecologies. To this end, we have developed an approach that utilizes soil mesocosms and integrated systems analyses to evaluate the efficacy of novel biocontainment strategies and to assess the impact of production systems upon terrestrial microbiome dynamics. We demonstrate the utility of this approach by modeling a contamination with industrial microbial chasses versus their biocontained counterparts. Here we demonstrate the broad utility of this system by highlighting findings from both strains of Saccharomyces cerevisiae that are contained with an inducible toxin anti-toxin system, and stains of Escherichia coli that are contained via genomic recoding. The resultant data demonstrate that this system has broad utility across diverse microbial chassis and biocontainment strategies, enables us to track the fate of our contaminating microbe with high sensitivity in the soil, as well as monitor broader impacts of the perturbation on the underlying soil system. The findings presented here support the use of this mesocosm-based approach to assess the environmental impact of industrial microbes and to validate biocontainment strategies.

BASIC BIOLOGICAL SCIENCES,INORGANIC, ORGANIC, PHYS↗

An automated workflow that generates atom mappings for large‐scale metabolic models and its application to Arabidopsis thaliana

SUMMARY Quantification of reaction fluxes of metabolic networks can help us understand how the integration of different metabolic pathways determines cellular functions. Yet, intracellular fluxes cannot be measured directly but are estimated with metabolic flux analysis (MFA), which relies on the patterns of isotope labeling of metabolites in the network. The application of MFA also requires a stoichiometric model with atom mappings that are currently not available for the majority of large‐scale metabolic network models, particularly of plants. While automated approaches such as the Reaction Decoder Toolkit (RDT) can produce atom mappings for individual reactions, tracing the flow of individual atoms of the entire reactions across a metabolic model remains challenging. Here we establish an automated workflow to obtain reliable atom mappings for large‐scale metabolic models by refining the outcome of RDT, and apply the workflow to metabolic models of Arabidopsis thaliana . We demonstrate the accuracy of RDT through a comparative analysis with atom mappings from a large database of biochemical reactions, MetaCyc. We further show the utility of our automated workflow by simulating 15 N isotope enrichment and identifying nitrogen (N)‐containing metabolites which show enrichment patterns that are informative for flux estimation in future 15 N‐MFA studies of A. thaliana . The automated workflow established in this study can be readily expanded to other species for which metabolic models have been established and the resulting atom mappings will facilitate MFA and graph‐theoretic structural analyses with large‐scale metabolic networks.

59 BASIC BIOLOGICAL SCIENCES↗

Integration of physiologically relevant photosynthetic energy flows into whole genome models of light‐driven metabolism

SUMMARY Characterizing photosynthetic productivity is necessary to understand the ecological contributions and biotechnology potential of plants, algae, and cyanobacteria. Light capture efficiency and photophysiology have long been characterized by measurements of chlorophyll fluorescence dynamics. However, these investigations typically do not consider the metabolic network downstream of light harvesting. By contrast, genome‐scale metabolic models capture species‐specific metabolic capabilities but have yet to incorporate the rapid regulation of the light harvesting apparatus. Here, we combine chlorophyll fluorescence parameters defining photosynthetic and non‐photosynthetic yield of absorbed light energy with a metabolic model of the pennate diatom Phaeodactylum tricornutum. This integration increases the model predictive accuracy regarding growth rate, intracellular oxygen production and consumption, and metabolic pathway usage. Through the quantification of excess electron transport, we uncover the sequential activation of non‐radiative energy dissipation processes, cross‐compartment electron shuttling, and non‐photochemical quenching as the rapid photoacclimation strategy in P. tricornutum. Interestingly, the photon absorption thresholds that trigger the transition between these mechanisms were consistent at low and high incident photon fluxes. We use this understanding to explore engineering strategies for rerouting cellular resources and excess light energy towards bioproducts in silico . Overall, we present a methodology for incorporating a common, informative data type into computational models of light‐driven metabolism and show its utilization within the design–build–test–learn cycle for engineering of photosynthetic organisms.

59 BASIC BIOLOGICAL SCIENCES↗

Metabolic model predictions enable targeted microbiome manipulation through precision prebiotics

While numerous health-beneficial interactions between host and microbiota have been identified, there is still a lack of targeted approaches for modulating these interactions. Thus, we here identify precision prebiotics that specifically modulate the abundance of a microbiome member species of interest. In the first step, we show that defining precision prebiotics by compounds that are only taken up by the target species but no other species in a community is usually not possible due to overlapping metabolic niches. Subsequently, we use metabolic modeling to identify precision prebiotics for a two-member Caenorhabditis elegans microbiome community comprising the immune-protective target species Pseudomonas lurida MYb11 and the persistent colonizer Ochrobactrum vermis MYb71. We experimentally confirm four of the predicted precision prebiotics, L-serine, L-threonine, D-mannitol, and γ-aminobutyric acid, to specifically increase the abundance of MYb11. L-serine was further assessed in vivo, leading to an increase in MYb11 abundance also in the worm host. Overall, our findings demonstrate that metabolic modeling is an effective tool for the design of precision prebiotics as an important cornerstone for future microbiome-targeted therapies.

59 BASIC BIOLOGICAL SCIENCES↗

Evaluating E. coli genome‐scale metabolic model accuracy with high‐throughput mutant fitness data

Abstract The Escherichia coli genome‐scale metabolic model (GEM) is an exemplar systems biology model for the simulation of cellular metabolism. Experimental validation of model predictions is essential to pinpoint uncertainty and ensure continued development of accurate models. Here, we quantified the accuracy of four subsequent E. coli GEMs using published mutant fitness data across thousands of genes and 25 different carbon sources. This evaluation demonstrated the utility of the area under a precision–recall curve relative to alternative accuracy metrics. An analysis of errors in the latest (iML1515) model identified several vitamins/cofactors that are likely available to mutants despite being absent from the experimental growth medium and highlighted isoenzyme gene‐protein‐reaction mapping as a key source of inaccurate predictions. A machine learning approach further identified metabolic fluxes through hydrogen ion exchange and specific central metabolism branch points as important determinants of model accuracy. This work outlines improved practices for the assessment of GEM accuracy with high‐throughput mutant fitness data and highlights promising areas for future model refinement in E. coli and beyond.

59 BASIC BIOLOGICAL SCIENCES↗