Search NASA⌕ Search

SEARCH · Search NASA

Results for “metabolic modeling”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Metabolic Model of the Nitrogen-Fixing Obligate Aerobe Azotobacter vinelandii Predicts Its Adaptation to Oxygen Concentration and Metal Availability

There is considerable interest in promoting biological nitrogen fixation (BNF) as a mechanism to reduce the inputs of nitrogenous fertilizers in agriculture, but considerable fundamental knowledge gaps still need to be addressed. BNF is catalyzed by nitrogenase, which requires a large input of energy in the form of ATP and low potential electrons. Diazotrophs that respire aerobically have an advantage in meeting the ATP demands of BNF but face challenges in protecting nitrogenase from inactivation by oxygen. Here, we constructed a genome-scale metabolic model of the nitrogen-fixing bacterium Azotobacter vinelandii, which uses a complex respiratory protection mechanism to consume oxygen at a high rate to keep intracellular conditions microaerobic. Our model accurately predicts growth rate under high oxygen and substrate concentrations, consistent with a large electron flux directed to the respiratory protection mechanism. While a partially decoupled electron transport chain compensates for some of the energy imbalance under high-oxygen conditions, it does not account for all substrate intake, leading to increased maintenance rates. Interestingly, the respiratory protection mechanism is required for accurate predictions even when ammonia is supplemented during growth, suggesting that the respiratory protection mechanism might be a core principle of metabolism and not just used for nitrogenase protection. We have also shown that rearrangement of flux through the electron transport system allows A. vinelandii to adapt to different oxygen concentrations, metal availability, and genetic disruption, which cause an ammonia excretion phenotype. Accurately determining the energy balance in an aerobic nitrogen-fixing metabolic model is required for future engineering approaches.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Recent advances in constraint and machine learning-based metabolic modeling by leveraging stoichiometric balances, thermodynamic feasibility and kinetic law formalisms

Understanding the governing principles behind organisms’ metabolism and growth underpins their effective deployment as bioproduction chassis. A central objective of metabolic modeling is predicting how metabolism and growth are affected by both external environmental factors and internal genotypic perturbations. The fundamental concepts of reaction stoichiometry, thermodynamics, and mass action kinetics have emerged as the foundational principles of many modeling frameworks designed to describe how and why organisms allocate resources towards both growth and bioproduction. Furthermore, this review focuses on the latest algorithmic advancements that have integrated these foundational principles into increasingly sophisticated quantitative frameworks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Quantifying microbial roles in environmental iron oxidation via an integrated kinetics, `omics and metabolic modeling study (Final Report)

Iron oxyhydroxides are extremely reactive components of environmental systems, and therefore exert a strong influence on biogeochemical cycles. These oxyhydroxides strongly adsorb many biologically-relevant elements, including organic carbon and phosphate, as well as a wide range of metals including uranium and actinide species. Thus, the formation mechanism of iron oxyhydroxides are key to understanding both nutrient and contaminant cycling. Microorganisms can catalyze iron oxidation and promote the formation of Fe biominerals and thus are increasingly recognized as important players in biogeochemical cycling. However, it is completely unknown how much of environmental iron oxidation is biologically mediated versus abiotic, and various challenges in studying microbial iron oxidation have hindered accurate incorporation into hydrobiogeochemical models. The overarching goal of our work was to quantify and constrain microbial iron oxidation rates and use ‘omics to gain insight into the controls on this process, while developing tools to enable integration of biotic iron oxidation into hydrobiogeochemical models. Our work focused on the Savannah River Site (SRS) in South Carolina, where extensive microbial iron oxidation has been observed. At Tims Branch, part of the Argonne National Laboratory Wetland Hydrobiogeochemistry Science Focus Area (Argonne SFA), where groundwater discharges into a stream, iron-oxidizing microbial mats form and appear to be a major sink of uranium. In the wetlands that surround Tims Branch, there are wide swaths of iron microbial mats and flocs (mobilized mat). We measured biotic and abiotic iron oxidation rates using mats and water sampled from these sites and found that iron oxidation is primarily carried out by chemolithotrophic microorganisms. The resulting rate constants can be incorporated into models. These mats were characterized by metagenomics and metatranscriptomics, which showed that aerobic chemolithotrophs were the dominant iron-oxidizing bacteria (FeOB), and these included Gallionellaceae and Leptothrix, and possibly Rhodoferax, which is known as an Fe-reducer but may also oxidize Fe(II). This demonstrated that diverse FeOB can coexist and suggests that there are a range of niches and therefore drivers of chemolithotrophic iron oxidation. Analysis of reconstructed genomes strongly suggests that a major factor in diversity is carbon source, as genomes contained varied pathways for autotrophy and heterotrophy. We performed an in-depth analysis of Leptothrix ochracea genomes, since this sheath-former is one of the primary mat builders, yet its physiology remained unresolved. A combination of genomics, transcriptomics, and metabolic modeling suggest that L. ochracea grows mixotrophically using a combination of Fe(II) and organics for energy and both inorganic and organic carbon to create biomass. This contrasts with the largely autotrophic Gallionellaceae (Gallionella, Sideroxydans, and Ferriphaselus) also present in the mats and flocs. Remarkably, multiple FeOB, both Leptothrix and Gallionellaceae, showed activity in response to Fe(II) in live mat incubations. We tracked the gene expression of individual MAGs to Fe(II) and found that various autotrophic and heterotrophic FeOB responded to Fe(II), increasing expression of both carbon fixation and organic utilization genes. The results of the integrated field, kinetics, and omics studies give detailed insight into 1) the taxa that oxidize Fe, and 2) how they connect Fe, C, and N cycles. Towards the goal of connecting omics data to hydrobiogeochemical models, we worked with the KBase team to create a template metabolic model for chemolithotrophic iron oxidation. We initially modeled the well-characterized isolate Gallionellaceae Sideroxydans lithotrophicus, and also applied the model to the mixotroph L. ochracea. In all, we have characterized diverse FeOB in a representative wetland system and solved key problems that enable better incorporation of iron-oxidizing microbes into hydrobiogeochemical models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Addressing uncertainty in genome-scale metabolic model reconstruction and analysis

The reconstruction and analysis of genome-scale metabolic models constitutes a powerful systems biology approach, with applications ranging from basic understanding of genotype-phenotype mapping to solving biomedical and environmental problems. However, the biological insight obtained from these models is limited by multiple heterogeneous sources of uncertainty, which are often difficult to quantify. Here we review the major sources of uncertainty and survey existing approaches developed for representing and addressing them. A unified formal characterization of these uncertainties through probabilistic approaches and ensemble modeling will facilitate convergence towards consistent reconstruction pipelines, improved data integration algorithms, and more accurate assessment of predictive capacity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Integrative teaching of metabolic modeling and flux analysis with interactive python modules

Abstract The modeling of rates of biochemical reactions—fluxes—in metabolic networks is widely used for both basic biological research and biotechnological applications. A number of different modeling methods have been developed to estimate and predict fluxes, including kinetic and constraint‐based (Metabolic Flux Analysis and flux balance analysis) approaches. Although different resources exist for teaching these methods individually, to‐date no resources have been developed to teach these approaches in an integrative way that equips learners with an understanding of each modeling paradigm, how they relate to one another, and the information that can be gleaned from each. We have developed a series of modeling simulations in Python to teach kinetic modeling, metabolic control analysis, 13C‐metabolic flux analysis, and flux balance analysis. These simulations are presented in a series of interactive notebooks with guided lesson plans and associated lecture notes. Learners assimilate key principles using models of simple metabolic networks by running simulations, generating and using data, and making and validating predictions about the effects of modifying model parameters. We used these simulations as the hands‐on computer laboratory component of a four‐day metabolic modeling workshop and participant survey results showed improvements in learners' self‐assessed competence and confidence in understanding and applying metabolic modeling techniques after having attended the workshop. The resources provided can be incorporated in their entirety or individually into courses and workshops on bioengineering and metabolic modeling at the undergraduate, graduate, or postgraduate level.

Kaste, Joshua A. M.↗
DOE OSTI2022

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↗
DOE OSTI2022

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↗
DOE OSTI2020

The ModelSEED Biochemistry Database for the integration of metabolic annotations and the reconstruction, comparison and analysis of metabolic models for plants, fungi and microbes

Abstract For over 10 years, ModelSEED has been a primary resource for the construction of draft genome-scale metabolic models based on annotated microbial or plant genomes. Now being released, the biochemistry database serves as the foundation of biochemical data underlying ModelSEED and KBase. The biochemistry database embodies several properties that, taken together, distinguish it from other published biochemistry resources by: (i) including compartmentalization, transport reactions, charged molecules and proton balancing on reactions; (ii) being extensible by the user community, with all data stored in GitHub; and (iii) design as a biochemical ‘Rosetta Stone’ to facilitate comparison and integration of annotations from many different tools and databases. The database was constructed by combining chemical data from many resources, applying standard transformations, identifying redundancies and computing thermodynamic properties. The ModelSEED biochemistry is continually tested using flux balance analysis to ensure the biochemical network is modeling-ready and capable of simulating diverse phenotypes. Ontologies can be designed to aid in comparing and reconciling metabolic reconstructions that differ in how they represent various metabolic pathways. ModelSEED now includes 33,978 compounds and 36,645 reactions, available as a set of extensible files on GitHub, and available to search at https://modelseed.org and KBase.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Advances in genome-scale metabolic models of industrially important fungi

Many fungal species have been used in industrial production for biofuels and bioproducts. Developing strains with better performance in biomanufacturing requires systematic understanding of cellular metabolism. Genome-scale metabolic models (GEMs) offer a comprehensive view of interconnected pathways and a mathematical framework for downstream analysis. Recently, GEMs have been developed or updated in several industrially important fungi. Some of them incorporate enzyme constraints, enabling improved predictions of cell states and proteome allocation. In this report we provide an overview of these newly developed GEMs and computational methods that facilitate construction of enzyme-constrained GEMs and utilize flux predictions from GEMs. Furthermore, we highlight the pivotal roles of these GEMs in iterative design-build-test- learn cycles, ultimately advancing the field of fungal biomanufacturing.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Sorghum bicolor Core Metabolism Model

Upregulation of triacylglycerols (TAGs) in vegetative plant tissues such as leaves has the potential to drastically increase the energy density and biomass yield of bioenergy crops. In this context, constraint-based analysis has the promise to improve metabolic engineering strategies. Here we present a core metabolism model for the C4 biomass crop Sorghum bicolor (iTJC1414) along with a minimal model for photosynthetic CO2 assimilation, sucrose and TAG biosynthesis in C3 plants. Extending iTJC1414 to a four-cell diel model we simulate C4 photosynthesis in mature leaves with the principal photo-assimilatory product being replaced by TAG produced at different levels. Independent of specific pathways and per unit carbon assimilated, energy content and biosynthetic demands in reducing equivalents are about 1.3 to 1.4 times higher for TAG than for sucrose. For plant generic pathways, ATP- and NADPH-demands per CO2 assimilated are higher by 1.3- and 1.5-fold, respectively. If the photosynthetic supply in ATP and NADPH in iTJC1414 is adjusted to be balanced for sucrose as the sole photo-assimilatory product, overproduction of TAG is predicted to cause a substantial surplus in photosynthetic ATP. This means that if TAG synthesis was the sole photo-assimilatory process, there could be an energy imbalance that might impede the process. Adjusting iTJC1414 to a photo-assimilatory rate that approximates field conditions, we predict possible daily rates of TAG accumulation, dependent on varying ratios of carbon partitioning between exported assimilates and accumulated oil droplets (TAG, oleosin) and in dependence of activation of futile cycles of TAG synthesis and degradation. We find that, based on the capacity of leaves for photosynthetic synthesis of exported assimilates, mature leaves should be able to reach a 20% level of TAG per dry weight within one month if only 5% of the photosynthetic net assimilation can be allocated into oil droplets. From this we conclude that high TAG levels should be achievable if TAG synthesis is induced only during a final phase of the plant life cycle.

Feedstock Production↗
DOE OSTI2021

Elucidation of trophic interactions in an unusual single-cell nitrogen-fixing symbiosis using metabolic modeling

Marine nitrogen-fixing microorganisms are an important source of fixed nitrogen in oceanic ecosystems. The colonial cyanobacterium Trichodesmium and diatom symbionts were thought to be the primary contributors to oceanic N 2 fixation until the discovery of the unusual uncultivated symbiotic cyanobacterium UCYN-A ( Candidatus Atelocyanobacterium thalassa ). UCYN-A has atypical metabolic characteristics lacking the oxygen-evolving photosystem II, the tricarboxylic acid cycle, the carbon-fixation enzyme RuBisCo and de novo biosynthetic pathways for a number of amino acids and nucleotides. Therefore, it is obligately symbiotic with its single-celled haptophyte algal host. UCYN-A receives fixed carbon from its host and returns fixed nitrogen, but further insights into this symbiosis are precluded by both UCYN-A and its host being uncultured. In order to investigate how this syntrophy is coordinated, we reconstructed bottom-up genome-scale metabolic models of UCYN-A and its algal partner to explore possible trophic scenarios, focusing on nitrogen fixation and biomass synthesis. Since both partners are uncultivated and only the genome sequence of UCYN-A is available, we used the phylogenetically related Chrysochromulina tobin as a proxy for the host. Through the use of flux balance analysis (FBA), we determined the minimal set of metabolites and biochemical functions that must be shared between the two organisms to ensure viability and growth. We quantitatively investigated the metabolic characteristics that facilitate daytime N 2 fixation in UCYN-A and possible oxygen-scavenging mechanisms needed to create an anaerobic environment to allow nitrogenase to function. This is the first application of an FBA framework to examine the tight metabolic coupling between uncultivated microbes in marine symbiotic communities and provides a roadmap for future efforts focusing on such specialized systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Integration of genome-scale metabolic model with biorefinery process model reveals market-competitive carbon-negative sustainable aviation fuel utilizing microbial cell mass lipids and biogenic CO 2

Producing scalable, economically viable, low-carbon biofuels or biochemicals hinges on more efficient bioconversion processes. While microbial conversion can offer robust solutions, the native microbial growth process often redirects a large fraction of carbon to CO 2 and cell mass. By integrating genome-scale metabolic models with techno-economic and life cycle assessment models, this study analyzes the effects of converting cell mass lipids to hydrocarbon fuels, and CO 2 to methanol on the facility’s costs and life-cycle carbon footprint. Results show that upgrading microbial lipids or both microbial lipids and CO 2 using renewable hydrogen produces carbon-negative bisabolene. Additionally, on-site electrolytic hydrogen production offers a supply of pure oxygen to use in place of air for bioconversion and fuel combustion in the boiler. To reach cost parity with conventional jet fuel, renewable hydrogen needs to be produced at less than $\$2.2$ to $\$3.1$/kg, with a bisabolene yield of 80% of the theoretical yield, along with cell mass and CO 2 yields of 22 wt% and 54 wt%, respectively. The economic combination of cell mass, CO 2 , and bisabolene yields demonstrated in this study provides practical insights for prioritizing research, selecting suitable hosts, and determining necessary engineered production levels.

09 BIOMASS FUELS↗
DOE OSTI2022

YEASTRACT+: a portal for the exploitation of global transcription regulation and metabolic model data in yeast biotechnology and pathogenesis

Abstract YEASTRACT+ (http://yeastract-plus.org/) is a tool for the analysis, prediction and modelling of transcription regulatory data at the gene and genomic levels in yeasts. It incorporates three integrated databases: YEASTRACT (http://yeastract-plus.org/yeastract/), PathoYeastract (http://yeastract-plus.org/pathoyeastract/) and NCYeastract (http://yeastract-plus.org/ncyeastract/), focused on Saccharomyces cerevisiae, pathogenic yeasts of the Candida genus, and non-conventional yeasts of biotechnological relevance. In this release, YEASTRACT+ offers upgraded information on transcription regulation for the ten previously incorporated yeast species, while extending the database to another pathogenic yeast, Candida auris. Since the last release of YEASTRACT+ (January 2020), a fourth database has been integrated. CommunityYeastract (http://yeastract-plus.org/community/) offers a platform for the creation, use, and future update of YEASTRACT-like databases for any yeast of the users’ choice. CommunityYeastract currently provides information for two Saccharomyces boulardii strains, Rhodotorula toruloides NP11 oleaginous yeast, and Schizosaccharomyces pombe 972h-. In addition, YEASTRACT+ portal currently gathers 304 547 documented regulatory associations between transcription factors (TF) and target genes and 480 DNA binding sites, considering 2771 TFs from 11 yeast species. A new set of tools, currently implemented for S. cerevisiae and C. albicans, is further offered, combining regulatory information with genome-scale metabolic models to provide predictions on the most promising transcription factors to be exploited in cell factory optimisation or to be used as novel drug targets. The expansion of these new tools to the remaining YEASTRACT+ species is ongoing.

Teixeira, Miguel Cacho (ORCID:0000000256766174)↗
DOE OSTI2023

Accurate flux predictions using tissue-specific gene expression in plant metabolic modeling

The accurate prediction of complex phenotypes such as metabolic fluxes in living systems is a grand challenge for systems biology and central to efficiently identifying biotechnological interventions that can address pressing industrial needs. The application of gene expression data to improve the accuracy of metabolic flux predictions using mechanistic modeling methods such as flux balance analysis (FBA) has not been previously demonstrated in multi-tissue systems, despite their biotechnological importance. We hypothesized that a method for generating metabolic flux predictions informed by relative expression levels between tissues would improve prediction accuracy. Relative gene expression levels derived from multiple transcriptomic and proteomic datasets were integrated into FBA predictions of a multi-tissue, diel model of Arabidopsis thaliana’s central metabolism. This integration dramatically improved the agreement of flux predictions with experimentally based flux maps from 13 C metabolic flux analysis compared with a standard parsimonious FBA approach. Disagreement between FBA predictions and MFA flux maps was measured using weighted averaged percent error values, and for parsimonious FBA this was 169%–180% for high light conditions and 94%–103% for low light conditions, depending on the gene expression dataset used. This fell to 10%-13% and 9%-11% upon incorporating expression data into the modeling process, which also substantially altered the predicted carbon and energy economy of the plant.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A Genome-Scale Metabolic Model of Anabaena 33047 to Guide Genetic Modifications to Overproduce Nylon Monomers

Nitrogen fixing-cyanobacteria can significantly improve the economic feasibility of cyanobacterial production processes by eliminating the requirement for reduced nitrogen. Anabaena sp. ATCC 33047 is a marine, heterocyst forming, nitrogen fixing cyanobacteria with a very short doubling time of 3.8 h. We developed a comprehensive genome-scale metabolic (GSM) model, iAnC892, for this organism using annotations and content obtained from multiple databases. iAnC892 describes both the vegetative and heterocyst cell types found in the filaments of Anabaena sp. ATCC 33047. iAnC892 includes 953 unique reactions and accounts for the annotation of 892 genes. Comparison of iAnC892 reaction content with the GSM of Anabaena sp. PCC 7120 revealed that there are 109 reactions including uptake hydrogenase, pyruvate decarboxylase, and pyruvate-formate lyase unique to iAnC892. iAnC892 enabled the analysis of energy production pathways in the heterocyst by allowing the cell specific deactivation of light dependent electron transport chain and glucose-6-phosphate metabolizing pathways. The analysis revealed the importance of light dependent electron transport in generating ATP and NADPH at the required ratio for optimal N2 fixation. When used alongside the strain design algorithm, OptForce, iAnC892 recapitulated several of the experimentally successful genetic intervention strategies that over produced valerolactam and caprolactam precursors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Investigating the Chemolithoautotrophic and Formate Metabolism of Nitrospira moscoviensis by Constraint-Based Metabolic Modeling and 13 C-Tracer Analysis

Nitrospira spp. are globally abundant nitrifying bacteria in soil and aquatic ecosystems and in wastewater treatment plants, where they control the oxidation of nitrite to nitrate. Despite their critical contribution to nitrogen cycling across diverse environments, detailed understanding of their metabolic network and prediction of their function under different environmental conditions remains a major challenge.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Modification and analysis of context-specific genome-scale metabolic models: methane-utilizing microbial chassis as a case study

ABSTRACT Context-specific genome-scale model (CS-GSM) reconstruction is becoming an efficient strategy for integrating and cross-comparing experimental multi-scale data to explore the relationship between cellular genotypes, facilitating fundamental or applied research discoveries. However, the application of CS modeling for non-conventional microbes is still challenging. Here, we present a graphical user interface that integrates COBRApy, EscherPy, and RIPTiDe, Python-based tools within the BioUML platform, and streamlines the reconstruction and interrogation of the CS genome-scale metabolic frameworks via Jupyter Notebook. The approach was tested using -omics data collected for Methylotuvimicrobium alcaliphilum 20Z R , a prominent microbial chassis for methane capturing and valorization. We optimized the previously reconstructed whole genome-scale metabolic network by adjusting the flux distribution using gene expression data. The outputs of the automatically reconstructed CS metabolic network were comparable to manually optimized i IA409 models for Ca-growth conditions. However, the CS model questions the reversibility of the phosphoketolase pathway and suggests higher flux via primary oxidation pathways. The model also highlighted unresolved carbon partitioning between assimilatory and catabolic pathways at the formaldehyde-formate node. Only a very few genes and only one enzyme with a predicted function in C1 metabolism, a homolog of the formaldehyde oxidation enzyme ( fae1-2 ), showed a significant change in expression in La-growth conditions. The CS-GSM predictions agreed with the experimental measurements under the assumption that the Fae1-2 is a part of the tetrahydrofolate-linked pathway. The cellular roles of the tungsten (W)-dependent formate dehydrogenase ( fdhAB ) and fae homologs ( fae1-2 and fae3 ) were investigated via mutagenesis. The phenotype of the f dhAB mutant followed the model prediction. Furthermore, a more significant reduction of the biomass yield was observed during growth in La-supplemented media, confirming a higher flux through formate. M. alcaliphilum 20Z R mutants lacking fae1-2 did not display any significant defects in methane or methanol-dependent growth. However, contrary to fae1, the fae1-2 homolog failed to restore the formaldehyde-activating enzyme function in complementation tests. Overall, the presented data suggest that the developed computational workflow supports the reconstruction and validation of CS-GSM networks of non-model microbes. IMPORTANCE The interrogation of various types of data is a routine strategy to explore the relationship between genotype and phenotype. An efficient approach for integrating and cross-comparing experimental multi-scale data in the context of whole-genome-based metabolic network reconstruction becomes a powerful tool that facilitates fundamental and applied research discoveries. The present study describes the reconstruction of a context-specific (CS) model for the methane-utilizing bacterium, Methylotuvimicrobium alcaliphilum 20Z R . M. alcaliphilum 20Z R is becoming an attractive microbial platform for the production of biofuels, chemicals, pharmaceuticals, and bio-sorbents for capturing atmospheric methane. We demonstrate that this pipeline can help reconstruct metabolic models that are similar to manually curated networks. Furthermore, the model is able to highlight previously overlooked pathways, thus advancing fundamental knowledge of non-model microbial systems or promoting their development toward biotechnological or environmental implementations.

Kulyashov, M. A.↗