Search

Trade-offs between traits are present across different levels of biological systems and ultimately reflect constraints imposed by physicochemical laws and the structure of underlying biochemical networks. Yet, mechanistic explanation of how trade-offs between molecular traits arise and how they relate to optimization of fitness-related traits remains elusive. Here, we introduce the concept of relative flux trade-offs and propose a constraint-based approach, termed FluTOr, to identify metabolic reactions whose fluxes are in relative trade-off with respect to an optimized fitness-related cellular task, like growth. We then employed FluTOr to identify relative flux trade-offs in the genome-scale metabolic networks of Escherichia coli, Saccharomyces cerevisiae, and Arabidopsis thaliana. For the metabolic models of E. coli and S. cerevisiae we showed that: (i) the identified relative flux trade-offs depend on the carbon source used and that (ii) reactions that participated in relative trade-offs in both species were implicated in cofactor biosynthesis. In contrast to the two microorganisms, the relative flux trade-offs for the metabolic model of A. thaliana did not depend on the available nitrogen sources, reflecting the differences in the underlying metabolic network as well as the considered environments. Lastly, the established connection between relative flux trade-offs allowed us to identify overexpression targets that can be used to optimize fitness-related traits. Altogether, our computational approach and findings demonstrate how relative flux trade-offs can shape optimization of metabolic tasks, important in biotechnological applications.

The chemistry underpinning microbial interactions provides an integrative framework for linking the activities of individual microbes, microbial communities, plants, and their environments. Currently, we know very little about the functions of genes and metabolites within these communities because genome annotations and functions are derived from the minority of microbes that have been propagated in the laboratory. Yet the diversity, complexity, inaccessibility, and irreproducibility of native microbial consortia limit our ability to interpret chemical signaling and map metabolic networks. In this perspective, we contend that standardized laboratory ecosystems are needed to dissect the chemistry of soil microbiomes. We argue that dissemination and application of standardized laboratory ecosystems will be transformative for the field, much like how model organisms have played critical roles in advancing biochemistry and molecular and cellular biology. Community consensus on fabricated ecosystems ("EcoFABs") along with protocols and data standards will integrate efforts and enable rapid improvements in our understanding of the biochemical ecology of microbial communities.

Lichens, algae and fungi-based symbiotic associations, are sources of many important secondary metabolites, such as antibiotics, anti-inflammatory, antioxidants, and anticancer agents. Wide range of experiments based on in vivo and in vitro studies revealed that lichens are a rich treasure of anti-cancer compounds. Lichen extracts and isolated lichen compounds can interact with all biological entities currently identified to be responsible for tumor development. The critical ways to control the cancer development include induction of cell cycle arrests, blocking communication of growth factors, activation of anti-tumor immunity, inhibition of tumor-friendly inflammation, inhibition of tumor metastasis, and suppressing chromosome dysfunction. Also, lichen-based compounds induce the killing of cells by the process of apoptosis, autophagy, and necrosis, that inturn positively modulates metabolic networks of cells against uncontrolled cell division. Many lichen-based compounds have proven to possess potential anti-cancer activity against a wide range of cancer cells, either alone or in conjunction with other anti-cancer compounds. This review primarily emphasizes on an updated account of the repository of secondary metabolites reported in lichens. Besides, we discuss the anti-cancer potential and possible mechanism of the most frequently reported secondary metabolites derived from lichens. [ABSTRACT FROM AUTHOR]

The metabolic networks underlying skotomorphogenesis in seedlings remain relatively unknown. On the basis of our previous study on the folate metabolism in seedlings grown in darkness, the plastidial folylpolyglutamate synthetase gene (AtDFB) T‐DNA insertion Arabidopsis thaliana mutant (atdfb‐3) was examined. Under the nitrate‐sufficient condition, the mutant exhibited deficient folate metabolism and hypocotyl elongation, which affected skotomorphogenesis. Further analyses revealed changes to multiple intermediate metabolites related to carbon and nitrogen metabolism in the etiolated atdfb‐3 seedlings. Specifically, the sugar, polyol, and fatty acid contents decreased in the atdfb‐3 mutant under the nitrate‐sufficient condition, whereas the abundance of various organic acids and amino acids increased. In response to nitrate‐limited stress, multiple metabolites, including sugars, polyols, fatty acids, organic acids, and amino acids, accumulated more in the mutant than in the wild‐type control. The differences in the contents of multiple metabolites between the atdfb‐3 and wild‐type seedlings decreased following the addition of exogenous 5‐F‐THF under both nitrogen conditions. Additionally, the mutant accumulated high levels of one‐carbon metabolites, such as Cys, S‐adenosylmethionine, and S‐adenosylhomocysteine, under both nitrogen conditions. Thus, our data demonstrated that the perturbed folate metabolism in the atdfb‐3 seedlings, which was caused by the loss‐of‐function mutation to AtDFB, probably altered carbon and nitrogen metabolism, thereby modulating skotomorphogenesis. Furthermore, the study findings provide new evidence of the links among folate metabolism, metabolic networks, and skotomorphogenesis. [ABSTRACT FROM AUTHOR]

Background: Sophisticated methods to properly pre-process and analyze the increasing collection of single-cell RNA sequencing (scRNA-seq) data are increasingly being developed. On the contrary, the best practices to integrate these data into metabolic networks, aiming at describing metabolic phenotypes within a heterogeneous cell population, have been poorly investigated. In this regard, a critical factor is the presence of false zero values in reactions essential for a fundamental metabolic function, such as biomass or energy production. Here, we investigate the role of denoising strategies in mitigating this problem. Methods: We applied state-of-the-art denoising strategies - namely MAGIC, ENHANCE, and SAVER - on three public scRNA-seq datasets. We then associated a metabolic flux distribution with every single cell by embedding its noise-free transcriptomics profile in the constraints of the optimization of a core metabolic model. Finally, we used the obtained single-cell optimal metabolic fluxes as features for cluster analysis. We compared the results obtained with different techniques, and with or without the use of denoising. We also investigated the possibility of applying denoising directly on the Reaction Activity Scores, which are metabolic features extracted from the read counts, rather than on the read counts. Results: We show that denoising of transcriptomics data improves the clustering of single cells. We also illustrate that denoising restores important metabolic properties, such as the correlation between cell cycle phase and biomass accumulation, and between the RAS scores of reactions belonging to the same metabolic pathway. We show that MAGIC performs better than ENHANCE and SAVER, and that, denoising applied directly on the RAS matrix could be an effective alternative in removing false zero values from essential metabolic reactions. Conclusions: Our results indicate that including denoising as a pre-processing operation represents a milestone to integrate scRNA-seq data into Flux Balance Analysis simulations and to perform single-cell cluster analysis with a focus on metabolic phenotypes. [ABSTRACT FROM AUTHOR]

Plant growth and development are characterized by systematic and continuous processes, each involving intricate metabolic coordination mechanisms. Mathematical models are essential tools for investigating plant growth and development, metabolic regulation networks, and growth patterns across different stages. These models offer insights into secondary metabolism patterns in plants and the roles of metabolites. The proliferation of data related to plant genomics, transcriptomics, proteomics, and metabolomics in the last decade has underscored the growing importance of mathematical modeling in this field. This review aims to elucidate the principles and types of metabolic models employed in studying plant secondary metabolism, their strengths, and limitations. Furthermore, the application of mathematical models in various plant systems biology subfields will be discussed. Lastly, the review will outline how mathematical models can be harnessed to address research questions in this context., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Gao and Zhao.)

Perhaps the most fundamental model in synthetic and systems biology for inferring pathways in metabolic reaction networks is a metabolic factory : a system of reactions that starts from a set of source compounds and produces a set of target molecules, while conserving or not depleting intermediate metabolites. Finding a shortest factory-that minimizes a sum of real-valued weights on its reactions to infer the most likely pathway-is NP-complete. The current state-of-the-art for shortest factories solves a mixed-integer linear program with a major drawback: it requires the user to set a critical parameter, where too large a value can make optimal solutions infeasible, while too small a value can yield degenerate solutions due to numerical error. We present the first robust algorithm for optimal factories that is both parameter-free (relieving the user from determining a parameter setting) and degeneracy-free (guaranteeing it finds an optimal nondegenerate solution). We also give for the first time a complete characterization of the graph-theoretic structure of shortest factories, that reveals an important class of degenerate solutions which was overlooked and potentially output by the prior state-of-the-art.We show degeneracy is precisely due to invalid stoichiometries in reactions, and provide an efficient algorithm for identifying all such misannotations in a metabolic network. In addition we settle the relationship between the two established pathway models of hyperpaths and factories by proving hyperpaths actually comprise a subclass of factories. Comprehensive experiments over all instances from the standard metabolic reaction databases in the literature demonstrate our parameter-free exact algorithm is fast in practice , quickly finding optimal factories in large real-world networks containing thousands of reactions. A preliminary implementation of our robust algorithm for shortest factories in a new tool called Freeia is available free for research use at http://freeia.cs.arizona.edu.

We present a tool for multi-omics data analysis that enables simultaneous visualization of up to four types of omics data on organism-scale metabolic network diagrams. The tool's interactive web-based metabolic charts depict the metabolic reactions, pathways, and metabolites of a single organism as described in a metabolic pathway database for that organism; the charts are constructed using automated graphical layout algorithms. The multi-omics visualization facility paints each individual omics dataset onto a different "visual channel" of the metabolic-network diagram. For example, a transcriptomics dataset might be displayed by coloring the reaction arrows within the metabolic chart, while a companion proteomics dataset is displayed as reaction arrow thicknesses, and a complementary metabolomics dataset is displayed as metabolite node colors. Once the network diagrams are painted with omics data, semantic zooming provides more details within the diagram as the user zooms in. Datasets containing multiple time points can be displayed in an animated fashion. The tool will also graph data values for individual reactions or metabolites designated by the user. The user can interactively adjust the mapping from data value ranges to the displayed colors and thicknesses to provide more informative diagrams., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Swart, Caspi, Paley and Karp.)

Abstract The metabolic networks underlying skotomorphogenesis in seedlings remain relatively unknown. On the basis of our previous study on the folate metabolism in seedlings grown in darkness, the plastidial folylpolyglutamate synthetase gene (AtDFB) T‐DNA insertion Arabidopsis thaliana mutant (atdfb‐3) was examined. Under the nitrate‐sufficient condition, the mutant exhibited deficient folate metabolism and hypocotyl elongation, which affected skotomorphogenesis. Further analyses revealed changes to multiple intermediate metabolites related to carbon and nitrogen metabolism in the etiolated atdfb‐3 seedlings. Specifically, the sugar, polyol, and fatty acid contents decreased in the atdfb‐3 mutant under the nitrate‐sufficient condition, whereas the abundance of various organic acids and amino acids increased. In response to nitrate‐limited stress, multiple metabolites, including sugars, polyols, fatty acids, organic acids, and amino acids, accumulated more in the mutant than in the wild‐type control. The differences in the contents of multiple metabolites between the atdfb‐3 and wild‐type seedlings decreased following the addition of exogenous 5‐F‐THF under both nitrogen conditions. Additionally, the mutant accumulated high levels of one‐carbon metabolites, such as Cys, S‐adenosylmethionine, and S‐adenosylhomocysteine, under both nitrogen conditions. Thus, our data demonstrated that the perturbed folate metabolism in the atdfb‐3 seedlings, which was caused by the loss‐of‐function mutation to AtDFB, probably altered carbon and nitrogen metabolism, thereby modulating skotomorphogenesis. Furthermore, the study findings provide new evidence of the links among folate metabolism, metabolic networks, and skotomorphogenesis.

Background: Recent studies have provided insights into the persistence and succession of microbes aboard the International Space Station (ISS), notably the dominance of Klebsiella pneumoniae. However, the interactions between the various microbes aboard the ISS and how they shape the microbiome remain to be clearly understood. In this study, we apply a computational approach to predict possible metabolic interactions in the ISS microbiome and shed further light on its organization. Results: Through a combination of a systems-based graph-theoretical approach, and a constraint-based community metabolic modeling approach, we demonstrated several key interactions in the ISS microbiome. These complementary approaches provided insights into the metabolic interactions and dependencies present amongst various microbes in a community, highlighting key interactions and keystone species. Our results showed that the presence of K. pneumoniae is beneficial to many other microorganisms it coexists with, notably those from the Pantoea genus. Species belonging to the Enterobacteriaceae family were often found to be the most beneficial for the survival of other microorganisms in the ISS microbiome. However, K. pneumoniae was found to exhibit parasitic and amensalistic interactions with Aspergillus and Penicillium species, respectively. To prove this metabolic prediction, K. pneumoniae and Aspergillus fumigatus were co-cultured under normal and simulated microgravity, where K. pneumoniae cells showed parasitic characteristics to the fungus. The electron micrography revealed that the presence of K. pneumoniae compromised the morphology of fungal conidia and degenerated its biofilm-forming structures. Conclusion: Our study underscores the importance of K. pneumoniae in the ISS, and its potential positive and negative interactions with other microbes, including potential pathogens. This integrated modeling approach, combined with experiments, demonstrates the potential for understanding the organization of other such microbiomes, unravelling key organisms and their interdependencies. -HMRmrRFAwPDF1qD73q3qY Video Abstract [ABSTRACT FROM AUTHOR]

Two novel natural products, the polyketide cuniculene and the peptide antibiotic aquimarin, were recently discovered from the marine bacterial genus Aquimarina. However, the diversity of the secondary metabolite biosynthetic gene clusters (SM-BGCs) in Aquimarina genomes indicates a far greater biosynthetic potential. In this study, nine representative Aquimarina strains were tested for antimicrobial activity against diverse human-pathogenic and marine microorganisms and subjected to metabolomic and genomic profiling. We found an inhibitory activity of most Aquimarina strains against Candida glabrata and marine Vibrio and Alphaproteobacteria species. Aquimarina sp. Aq135 and Aquimarina muelleri crude extracts showed particularly promising antimicrobial activities, amongst others against methicillin-resistant Staphylococcus aureus. The metabolomic and functional genomic profiles of Aquimarina spp. followed similar patterns and were shaped by phylogeny. SM-BGC and metabolomics networks suggest the presence of novel polyketides and peptides, including cyclic depsipeptide-related compounds. Moreover, exploration of the 'Sponge Microbiome Project' dataset revealed that Aquimarina spp. possess low-abundance distributions worldwide across multiple marine biotopes. Our study emphasizes the relevance of this member of the microbial rare biosphere as a promising source of novel natural products. We predict that future metabologenomics studies of Aquimarina species will expand the spectrum of known secondary metabolites and bioactivities from marine ecosystems. [ABSTRACT FROM AUTHOR]

Genome-scale metabolism can best be described as a highly interconnected network of biochemical reactions and metabolites. The flow of metabolites, i.e., flux, throughout these networks can be predicted and analyzed using approaches such as flux balance analysis (FBA). By knowing the network topology and employing only a few simple assumptions, FBA can efficiently predict metabolic functions at the genome scale as well as microbial phenotypes. The network topology is represented in the form of genome-scale metabolic models (GEMs), which provide a direct mapping between network structure and function via the enzyme-coding genes and corresponding metabolic capacity. Recently, the role of protein limitations in shaping metabolic phenotypes have been extensively studied following the reconstruction of enzyme-constrained GEMs. This framework has been shown to significantly improve the accuracy of predicting microbial phenotypes, and it has demonstrated that a global limitation in protein availability can prompt the ubiquitous metabolic strategy of overflow metabolism. Being one of the most abundant and differentially expressed proteome sectors, metabolic proteins constitute a major cellular demand on proteinogenic amino acids. However, little is known about the impact and sensitivity of amino acid availability with regards to genome-scale metabolism. Here, we explore these aspects by extending on the enzyme-constrained GEM framework by also accounting for the usage of amino acids in expressing the metabolic proteome. Including amino acids in an enzyme-constrained GEM of Saccharomyces cerevisiae, we demonstrate that the expanded model is capable of accurately reproducing experimental amino acid levels. We further show that the metabolic proteome exerts variable demands on amino acid supplies in a condition-dependent manner, suggesting that S. cerevisiae must have evolved to efficiently fine-tune the synthesis of amino acids for expressing its metabolic proteins in response to changes in the external environment. Finally, our results demonstrate how the metabolic network of S. cerevisiae is robust towards perturbations of individual amino acids, while simultaneously being highly sensitive when the relative amino acid availability is set to mimic a priori distributions of both yeast and non-yeast origins.

Abstract Background Computational mathematical models of biological and biomedical systems have been successfully applied to advance our understanding of various regulatory processes, metabolic fluxes, effects of drug therapies, and disease evolution and transmission. Unfortunately, despite community efforts leading to the development of SBML and the BioModels database, many published models have not been fully exploited, largely due to a lack of proper documentation or the dependence on proprietary software. To facilitate the reuse and further development of systems biology and systems medicine models, an open-source toolbox that makes the overall process of model construction more consistent, understandable, transparent, and reproducible is desired. Results and discussion We provide an update on the development of modelbase, a free, expandable Python package for constructing and analysing ordinary differential equation-based mathematical models of dynamic systems. It provides intuitive and unified methods to construct and solve these systems. Significantly expanded visualisation methods allow for convenient analysis of the structural and dynamic properties of models. After specifying reaction stoichiometries and rate equations modelbase can automatically assemble the associated system of differential equations. A newly provided library of common kinetic rate laws reduces the repetitiveness of the computer programming code. modelbase is also fully compatible with SBML. Previous versions provided functions for the automatic construction of networks for isotope labelling studies. Now, using user-provided label maps, modelbase v1.2.3 streamlines the expansion of classic models to their isotope-specific versions. Finally, the library of previously published models implemented in modelbase is growing continuously. Ranging from photosynthesis to tumour cell growth to viral infection evolution, all these models are now available in a transparent, reusable and unified format through modelbase. Conclusion With this new Python software package, which is written in currently one of the most popular programming languages, the user can develop new models and actively profit from the work of others. modelbase enables reproducing and replicating models in a consistent, tractable and expandable manner. Moreover, the expansion of models to their isotopic label-specific versions enables simulating label propagation, thus providing quantitative information regarding network topology and metabolic fluxes.

Abstract Genome-scale metabolic models of microorganisms are powerful frameworks to predict phenotypes from an organism’s genotype. While manual reconstructions are laborious, automated reconstructions often fail to recapitulate known metabolic processes. Here we present gapseq ( https://github.com/jotech/gapseq ), a new tool to predict metabolic pathways and automatically reconstruct microbial metabolic models using a curated reaction database and a novel gap-filling algorithm. On the basis of scientific literature and experimental data for 14,931 bacterial phenotypes, we demonstrate that gapseq outperforms state-of-the-art tools in predicting enzyme activity, carbon source utilisation, fermentation products, and metabolic interactions within microbial communities.

Abstract Background Significant efforts have been made in building large-scale kinetic models of cellular metabolism in the past two decades. However, most kinetic models published to date, remain focused around central carbon pathways or are built around ad hoc reduced models without clear justification on their derivation and usage. Systematic algorithms exist for reducing genome-scale metabolic reconstructions to build thermodynamically feasible and consistently reduced stoichiometric models. However, it is important to study how network complexity affects conclusions derived from large-scale kinetic models built around consistently reduced models before we can apply them to study biological systems. Results We reduced the iJO1366 Escherichia Coli genome-scale metabolic reconstruction systematically to build three stoichiometric models of different size. Since the reduced models are expansions around the core subsystems for which the reduction was performed, the models are nested. We present a method for scaling up the flux profile and the concentration vector reference steady-states from the smallest model to the larger ones, whilst preserving maximum equivalency. Populations of kinetic models, preserving similarity in kinetic parameters, were built around the reference steady-states and their metabolic sensitivity coefficients (MSCs) were computed. The MSCs were sensitive to the model complexity. We proposed a metric for measuring the sensitivity of MSCs to these structural changes. Conclusions We proposed for the first time a workflow for scaling up the size of kinetic models while preserving equivalency between the kinetic models. Using this workflow, we demonstrate that model complexity in terms of networks size has significant impact on sensitivity characteristics of kinetic models. Therefore, it is essential to account for the effects of network complexity when constructing kinetic models. The presented metric for measuring MSC sensitivity to structural changes can guide modelers and experimentalists in improving model quality and guide synthetic biology and metabolic engineering. Our proposed workflow enables the testing of the suitability of a kinetic model for answering certain study-specific questions. We argue that the model-based metabolic design targets that are common across models of different size are of higher confidence, while those that are different could be the objective of investigations for model improvement.

Strain optimization aims to overproduce valuable metabolites by leveraging an understanding of biological systems, including metabolic networks and gene regulatory networks (GRNs). Accordingly, researchers proposed integrating metabolic networks and GRNs to be analyzed simultaneously. The proposed algorithms from 2002 to 2021 were rFBA, SR-FBA, iFBA, PROM, PROM2.0, TIGER, BeReTa, CoRegFlux, IDREAM, TRFBA, OptRAM, TRIMER, and PRIME. Each algorithm has different characteristics. Thus, using the appropriate algorithm for designing strain optimization is essential. Therefore, a critical review was conducted by synthesizing and analyzing the existing algorithms. Five aspects are discussed in this review: the strategic approaches, model of GRNs, source of GRNs, optimization, supplementary methods, and the programming language used. Based on the review, several algorithms were better at modeling integrated regulatory-metabolic networks with high confidence, i.e., PROM, PROM2.0, and TRFBA. A simulation was applied to six strains. The results show that PROM2.0 best predicted the production rate and time complexity. However, the model is heavily influenced by the quality and quantity of the gene expression data. In addition, there are inconsistencies between GRNs and the gene expression data. Thus, this review also discussed future work based on GRNs and gene expression data. [ABSTRACT FROM AUTHOR]

Normal aging causes profound changes of structural degeneration and glucose hypometabolism in the human brain, even in the absence of disease. In recent years, with the extensive exploration of the topological characteristics of the human brain, related studies in rats have begun to investigate. However, age-related alterations of topological properties in individual brain metabolic network of rats remain unknown. In this study, a total of 48 healthy female Sprague–Dawley (SD) rats were used, including 24 young rats and 24 aged rats. We used Jensen-Shannon Divergence Similarity Estimation (JSSE) method for constructing individual metabolic networks to explore age-related topological properties and rich-club organization changes. Compared with the young rats, the aged rats showed significantly decreased clustering coefficient (Cp) and local efficiency (E loc ) across the whole-brain metabolic network. In terms of changes in local network measures, degree (D) and nodal efficiency (E nod ) of left posterior dorsal hippocampus, and E nod of left olfactory tubercle were higher in the aged rats than in the young rats. About the rich-club analysis, the existence of rich-club organization in individual brain metabolic networks of rats was demonstrated. In addition, our findings further confirmed that rich-club connections were susceptible to aging. Relative to the young rats, the overall strength of rich-club connections was significantly reduced in the aged rats, while the overall strength of feeder and local connections was significantly increased. These findings demonstrated the age-related reorganization principle of the brain structure and improved our understanding of brain alternations during aging. [ABSTRACT FROM AUTHOR]

Mathematical modeling of metabolic networks is a powerful approach to investigate the underlying principles of metabolism and growth. Such approaches include, among others, differential-equation-based modeling of metabolic systems, constraint-based modeling and metabolic network expansion of metabolic networks. Most of these methods are well established and are implemented in numerous software packages, but these are scattered between different programming languages, packages and syntaxes. This complicates establishing straight forward pipelines integrating model construction and simulation. We present a Python package moped that serves as an integrative hub for reproducible construction, modification, curation and analysis of metabolic models. moped supports draft reconstruction of models directly from genome/proteome sequences and pathway/genome databases utilizing GPR annotations, providing a completely reproducible model construction and curation process within executable Python scripts. Alternatively, existing models published in SBML format can be easily imported. Models are represented as Python objects, for which a wide spectrum of easy-to-use modification and analysis methods exist. The model structure can be manually altered by adding, removing or modifying reactions, and gap-filling reactions can be found and inspected. This greatly supports the development of draft models, as well as the curation and testing of models. Moreover, moped provides several analysis methods, in particular including the calculation of biosynthetic capacities using metabolic network expansion. The integration with other Python-based tools is facilitated through various model export options. For example, a model can be directly converted into a CobraPy object for constraint-based analyses. moped is a fully documented and expandable Python package. We demonstrate the capability to serve as a hub for integrating reproducible model construction and curation, database import, metabolic network expansion and export for constraint-based analyses. [ABSTRACT FROM AUTHOR]

Mycobacterium tuberculosis, the causative agent of tuberculosis, is composed of several lineages characterized by a genome identity higher than 99%. Although the majority of the lineages are associated with humans, at least four lineages are adapted to other mammals, including different M. tuberculosis ecotypes. Host specificity is associated with higher virulence in its preferred host in ecotypes such as M. bovis. Deciphering what determines the preference of the host can reveal host-specific virulence patterns. However, it is not clear which genomic determinants might be influencing host specificity. In this study, we apply a combination of unsupervised and supervised classification methods on genomic data of ~27,000 M. tuberculosis clinical isolates to decipher host-specific genomic determinants. Host-specific genomic signatures are scarce beyond known lineage-specific mutations. Therefore, we integrated lineage-specific mutations into the iEK1011 2.0 genome-scale metabolic model to obtain lineage-specific versions of it. Flux distributions sampled from the solution spaces of these models can be accurately separated according to host association. This separation correlated with differences in cell wall processes, lipid, amino acid and carbon metabolic subsystems. These differences were observable when more than 95% of the samples had a specific growth rate significantly lower than the maximum achievable by the models. This suggests that these differences might manifest at low growth rate settings, such as the restrictive conditions M. tuberculosis suffers during macrophage infection. [ABSTRACT FROM AUTHOR]

The ever-increasing importance of structured data in different applications, especially in the biomedical field, has driven the need for reducing its complexity through projections into a more manageable space. The latest methods for learning features on graphs focus mainly on the neighborhood of nodes and edges. Methods capable of providing a representation that looks beyond the single node neighborhood are kernel graphs. However, they produce handcrafted features unaccustomed with a generalized model. To reduce this gap, in this work we propose a neural embedding framework, based on probability distribution representations of graphs, named Netpro2vec. The goal is to look at basic node descriptions other than the degree, such as those induced by the Transition Matrix and Node Distance Distribution. Netpro2vec provides embeddings completely independent from the task and nature of the data. The framework is evaluated on synthetic and various real biomedical network datasets through a comprehensive experimental classification phase and is compared to well-known competitors. [ABSTRACT FROM AUTHOR]

Human serine hydroxymethyltransferase (SHMT) regulates the serine-glycine one carbon metabolism and plays a role in cancer metabolic reprogramming. Two SHMT isozymes are acting in the cell: SHMT1 encoding the cytoplasmic isozyme, and SHMT2 encoding the mitochondrial one. Here we present a molecular model built on experimental data reporting the interaction between SHMT1 protein and SHMT2 mRNA, recently discovered in lung cancer cells. Using a stochastic dynamic model, we show that RNA moieties dynamically regulate serine and glycine concentration, shaping the system behaviour. For the first time we observe an active functional role of the RNA in the regulation of the serine-glycine metabolism and availability, which unravels a complex layer of regulation that cancer cells exploit to fine tune amino acids availability according to their metabolic needs. The quantitative model, complemented by an experimental validation in the lung adenocarcinoma cell line H1299, exploits RNA molecules as metabolic switches of the SHMT1 activity. Our results pave the way for the development of RNA-based molecules able to unbalance serine metabolism in cancer cells.

Subject of Research. The joint clustering method in network and correlation context is designed to identify active modules in metabolic graphs based on transcriptomic data represented by a large number of samples. The active modules obtained by this method describe the dynamic metabolic regulation in all samples of the analyzed dataset. The paper presents modifications of the proposed method for application on real data. Method. For results stability study the modified method was repeatedly run on real data with small variations of the initial parameters. For result analysis, several metrics were formulated that display modules similarity and representation under various start-up conditions. Main Results. The analysis results are sufficiently robust. For the most modules, their profiles are detected well in the noisy data, and the most genes are also preserved. Practical Relevance. The results of the presented study have shown that the modified method analyzes successfully real data by producing active modules that are stable and easy in interpretation.

Normal aging causes profound changes of structural degeneration and glucose hypometabolism in the human brain, even in the absence of disease. In recent years, with the extensive exploration of the topological characteristics of the human brain, related studies in rats have begun to investigate. However, age-related alterations of topological properties in individual brain metabolic network of rats remain unknown. In this study, a total of 48 healthy female Sprague–Dawley (SD) rats were used, including 24 young rats and 24 aged rats. We used Jensen-Shannon Divergence Similarity Estimation (JSSE) method for constructing individual metabolic networks to explore age-related topological properties and rich-club organization changes. Compared with the young rats, the aged rats showed significantly decreased clustering coefficient (Cp) and local efficiency (Eloc) across the whole-brain metabolic network. In terms of changes in local network measures, degree (D) and nodal efficiency (Enod) of left posterior dorsal hippocampus, and Enod of left olfactory tubercle were higher in the aged rats than in the young rats. About the rich-club analysis, the existence of rich-club organization in individual brain metabolic networks of rats was demonstrated. In addition, our findings further confirmed that rich-club connections were susceptible to aging. Relative to the young rats, the overall strength of rich-club connections was significantly reduced in the aged rats, while the overall strength of feeder and local connections was significantly increased. These findings demonstrated the age-related reorganization principle of the brain structure and improved our understanding of brain alternations during aging.

Motivation: The increasing availability of metabolomic data and their analysis are improving the understanding of cellular mechanisms and how biological systems respond to different perturbations. Currently, there is a need for novel computational methods that facilitate the analysis and integration of increasing volume of available data. Results: In this paper, we present Totoro a new constraint-based approach that integrates quantitative non-targeted metabolomic data of two different metabolic states into genome-wide metabolic models and predicts reactions that were most likely active during the transient state. We applied Totoro to real data of three different growth experiments (pulses of glucose, pyruvate, succinate) from Escherichia coli and we were able to predict known active pathways and gather new insights on the different metabolisms related to each substrate. We used both the E. coli core and the iJO1366 models to demonstrate that our approach is applicable to both smaller and larger networks. Availability: Totoro is an open source method (available at https://gitlab.inria.fr/erable/totoro) suitable for any organism with an available metabolic model. It is implemented in C++ and depends on IBM CPLEX which is freely available for academic purposes. [ABSTRACT FROM AUTHOR]

Hevea brasiliensis (rubber tree) is a large tree species of the Euphorbiaceae family with inestimable economic importance. Rubber tree breeding programs currently aim to improve growth and production, and the use of early genotype selection technologies can accelerate such processes, mainly with the incorporation of genomic tools, such as marker-assisted selection (MAS). However, few quantitative trait loci (QTLs) have been used successfully in MAS for complex characteristics. Recent research shows the efficiency of genome-wide association studies (GWAS) for locating QTL regions in different populations. In this way, the integration of GWAS, RNA-sequencing (RNA-Seq) methodologies, coexpression networks and enzyme networks can provide a better understanding of the molecular relationships involved in the definition of the phenotypes of interest, supplying research support for the development of appropriate genomic based strategies for breeding. In this context, this work presents the potential of using combined multiomics to decipher the mechanisms of genotype and phenotype associations involved in the growth of rubber trees. Using GWAS from a genotyping-by-sequencing (GBS) Hevea population, we were able to identify molecular markers in QTL regions with a main effect on rubber tree plant growth under constant water stress. The underlying genes were evaluated and incorporated into a gene coexpression network modelled with an assembled RNA-Seq-based transcriptome of the species, where novel gene relationships were estimated and evaluated through in silico methodologies, including an estimated enzymatic network. From all these analyses, we were able to estimate not only the main genes involved in defining the phenotype but also the interactions between a core of genes related to rubber tree growth at the transcriptional and translational levels. This work was the first to integrate multiomics analysis into the in-depth investigation of rubber tree plant growth, producing useful data for future genetic studies in the species and enhancing the efficiency of the species improvement programs. [ABSTRACT FROM AUTHOR]

Equilibrium bifurcations arise from sign changes of Jacobian determinants, as parameters are varied. Therefore we address the Jacobian determinant for metabolic networks with general reaction kinetics. Our approach is based on the concept of Child Selections: each (mother) metabolite is mapped, injectively, to one of those (child) reactions that it drives as an input. Our analysis distinguishes reaction network Jacobians with constant sign from the bifurcation case, where that sign depends on specific reaction rates. In particular, we distinguish "good" Child Selections, which do not affect the sign, from more interesting and mischievous "bad" children, which gang up towards sign changes, instability, and bifurcation.

Abstract Background The concept of minimal cut sets (MCS) has become an important mathematical framework for analyzing and (re)designing metabolic networks. However, the calculation of MCS in genome-scale metabolic models is a complex computational problem. The development of duality-based algorithms in the last years allowed the enumeration of thousands of MCS in genome-scale networks by solving mixed-integer linear problems (MILP). A recent advancement in this field was the introduction of the MCS2 approach. In contrast to the Farkas-lemma-based dual system used in earlier studies, the MCS2 approach employs a more condensed representation of the dual system based on the nullspace of the stoichiometric matrix, which, due to its reduced dimension, holds promise to further enhance MCS computations. Results In this work, we introduce several new variants and modifications of duality-based MCS algorithms and benchmark their effects on the overall performance. As one major result, we generalize the original MCS2 approach (which was limited to blocking the operation of certain target reactions) to the most general case of MCS computations with arbitrary target and desired regions. Building upon these developments, we introduce a new MILP variant which allows maximal flexibility in the formulation of MCS problems and fully leverages the reduced size of the nullspace-based dual system. With a comprehensive set of benchmarks, we show that the MILP with the nullspace-based dual system outperforms the MILP with the Farkas-lemma-based dual system speeding up MCS computation with an averaged factor of approximately 2.5. We furthermore present several simplifications in the formulation of constraints, mainly related to binary variables, which further enhance the performance of MCS-related MILP. However, the benchmarks also reveal that some highly condensed formulations of constraints, especially on reversible reactions, may lead to worse behavior when compared to variants with a larger number of (more explicit) constraints and involved variables. Conclusions Our results further enhance the algorithmic toolbox for MCS calculations and are of general importance for theoretical developments as well as for practical applications of the MCS framework.

The study of microbiomes by sequencing has revealed a plethora of correlations between microbial community composition and various life-history characteristics of the corresponding host species. However, inferring causation from correlation is often hampered by the sheer compositional complexity of microbiomes, even in simple organisms. Synthetic communities offer an effective approach to infer cause-effect relationships in host-microbiome systems. Yet the available communities suffer from several drawbacks, such as artificial (thus non-natural) choice of microbes, microbe-host mismatch (e.g., human microbes in gnotobiotic mice), or hosts lacking genetic tractability. Here we introduce CeMbio, a simplified natural Caenorhabditis elegans microbiota derived from our previous meta-analysis of the natural microbiome of this nematode. The CeMbio resource is amenable to all strengths of the C. elegans model system, strains included are readily culturable, they all colonize the worm gut individually, and comprise a robust community that distinctly affects nematode life-history. Several tools have additionally been developed for the CeMbio strains, including diagnostic PCR primers, completely sequenced genomes, and metabolic network models. With CeMbio, we provide a versatile resource and toolbox for the in-depth dissection of naturally relevant host-microbiome interactions in C. elegans.

Abstract Background Biological networks are representative of the diverse molecular interactions that occur within cells. Some of the commonly studied biological networks are modeled through protein-protein interactions, gene regulatory, and metabolic pathways. Among these, metabolic networks are probably the most studied, as they directly influence all physiological processes. Exploration of biochemical pathways using multigraph representation is important in understanding complex regulatory mechanisms. Feature extraction and clustering of these networks enable grouping of samples obtained from different biological specimens. Clustering techniques separate networks depending on their mutual similarity. Results We present a clustering analysis on tissue-specific metabolic networks for single samples from three primary tumor sites: breast, lung, and kidney cancer. The metabolic networks were obtained by integrating genome scale metabolic models with gene expression data. We performed network simplification to reduce the computational time needed for the computation of network distances. We empirically proved that networks clustering can characterize groups of patients in multiple conditions. Conclusions We provide a computational methodology to explore and characterize the metabolic landscape of tumors, thus providing a general methodology to integrate analytic metabolic models with gene expression data. This method represents a first attempt in clustering large scale metabolic networks. Moreover, this approach gives the possibility to get valuable information on what are the effects of different conditions on the overall metabolism.

Increased technological methods have enabled the investigation of biology at nanoscale levels. Such systems require the use of computational methods to comprehend the complex interactions that occur. The dynamics of metabolic systems have been traditionally described utilizing differential equations without fully capturing the heterogeneity of biological systems. Stochastic modeling approaches have recently emerged with the capacity to incorporate the statistical properties of such systems. However, the processing of stochastic algorithms is a computationally intensive task with intrinsic limitations. Alternatively, the queueing theory approach, historically used in the evaluation of telecommunication networks, can significantly reduce the computational power required to generate simulated results while simultaneously reducing the expansion of errors. We present here the application of queueing theory to simulate stochastic metabolic networks with high efficiency. With the use of glycolysis as a well understood biological model, we demonstrate the power of the proposed modeling methods discussed herein. Furthermore, we describe the simulation and pharmacological inhibition of glycolysis to provide an example of modeling capabilities.

Motivation: The increasing availability of metabolomic data and their analysis are improving the understanding of cellular mechanisms and how biological systems respond to different perturbations. Currently, there is a need for novel computational methods that facilitate the analysis and integration of increasing volume of available data.Results: In this paper, we present Totoro a new constraint-based approach that integrates quantitative non-targeted metabolomic data of two different metabolic states into genome-wide metabolic models and predicts reactions that were most likely active during the transient state. We applied Totoro to real data of three different growth experiments (pulses of glucose, pyruvate, succinate) from Escherichia coli and we were able to predict known active pathways and gather new insights on the different metabolisms related to each substrate. We used both the E. coli core and the iJO1366 models to demonstrate that our approach is applicable to both smaller and larger networks.Availability:Totoro is an open source method (available at https://gitlab.inria.fr/erable/totoro) suitable for any organism with an available metabolic model. It is implemented in C++ and depends on IBM CPLEX which is freely available for academic purposes.

Studying the correlation between microbiome metabolism and flavor of fermented foods has garnered significant attention recently. Understanding the contribution of metabolic regulation and environmental stress to microecosystems is essential for exploring the mechanisms of action of traditional fermented foods. Here, the interaction between microbial communities was investigated using a Xiaoqu fermentation system, processed as “simulative microecosystems,” in which starters were composed of Rhizopus-specific species/strains, Meyerozyma guilliermondii, and Bacillus licheniformis. The differences between community succession and metabolites were also explored. The results indicated that Rhizopus species/strain specificity affected starch hydrolyzation, resulting in a remarkable difference in the type and content of organic acids. This further suggested that the differences in nutrient abundance and organic acids influenced the colonization of microorganisms in the fermentation system, thereby influencing the succession of their communities. The fungi in the community predominantly originated from starters, whereas the bacteria were derived from both the environment and starter. Environmentally colonized microbes were the major contributors to the co-occurrence network and were strongly correlated with network. Regional characteristics of fermented foods were closely related to environmental microbes. These results contribute to the understanding of microbial assembly and flavor metabolism in fermented foods and provide strategies for quality regulation.

Hevea brasiliensis (rubber tree) is a large tree species of the Euphorbiaceae family with inestimable economic importance. Rubber tree breeding programs currently aim to improve growth and production, and the use of early genotype selection technologies can accelerate such processes, mainly with the incorporation of genomic tools, such as marker-assisted selection (MAS). However, few quantitative trait loci (QTLs) have been used successfully in MAS for complex characteristics. Recent research shows the efficiency of genome-wide association studies (GWAS) for locating QTL regions in different populations. In this way, the integration of GWAS, RNA-sequencing (RNA-Seq) methodologies, coexpression networks and enzyme networks can provide a better understanding of the molecular relationships involved in the definition of the phenotypes of interest, supplying research support for the development of appropriate genomic based strategies for breeding. In this context, this work presents the potential of using combined multiomics to decipher the mechanisms of genotype and phenotype associations involved in the growth of rubber trees. Using GWAS from a genotyping-by-sequencing (GBS) Hevea population, we were able to identify molecular markers in QTL regions with a main effect on rubber tree plant growth under constant water stress. The underlying genes were evaluated and incorporated into a gene coexpression network modelled with an assembled RNA-Seq-based transcriptome of the species, where novel gene relationships were estimated and evaluated through in silico methodologies, including an estimated enzymatic network. From all these analyses, we were able to estimate not only the main genes involved in defining the phenotype but also the interactions between a core of genes related to rubber tree growth at the transcriptional and translational levels. This work was the first to integrate multiomics analysis into the in-depth investigation of rubber tree plant growth, producing useful data for future genetic studies in the species and enhancing the efficiency of the species improvement programs.

• Abstract interpretation allows gene knock out predictions based on elementary modes. • Qualitative analysis of metabolic networks realized through abstract interpretation. • Elementary modes provide Boolean analysis of reaction networks. Abstract interpretation was proposed for predicting changes of reaction networks with partial kinetic information in systems biology. This requires to compute the set of difference abstractions of a system of linear equations under nonlinear constraints. We present the first practical algorithm that can compute the difference abstractions of linear equation systems exactly. We also present a new heuristics based on minimal support consequences for overapproximating the set of difference abstractions. Our algorithms rely on elementary modes, first-order definitions, and finite domain constraint programming. We implemented our algorithms and applied them to change prediction in systems biology. It turns out experimentally that the new heuristics is often exact in practice, while outperforming the exact algorithm. This journal article extends on a paper published at the 17th International Conference on Computational Methods in Systems Biology (CMSB'2019) [1]. [ABSTRACT FROM AUTHOR]

Metabolism is the set of biochemical reactions of an organism that enables it to assimilate nutrients from its environment and to generate building blocks for growth and proliferation. It forms a complex network that is intertwined with the many molecular and cellular processes that take place within cells. Systems biology aims to capture the complexity of cells, organisms, or communities by reconstructing models based on information gathered by high-throughput analyses (omics data) and prior knowledge. One type of model is a genome-scale metabolic model (GEM) that allows studying the distributions of metabolic fluxes, i.e., the "mass-flow" through the network of biochemical reactions. GEMs are nowadays widely applied and have been reconstructed for various microbial pathogens, either in a free-living state or in interaction with their hosts, with the aim to gain insight into mechanisms of pathogenicity. In this review, we first introduce the principles of systems biology and GEMs. We then describe how metabolic modeling can contribute to unraveling microbial pathogenesis and host–pathogen interactions, with a specific focus on oomycete plant pathogens and in particular Phytophthora infestans. Subsequently, we review achievements obtained so far and identify and discuss potential pitfalls of current models. Finally, we propose a workflow for reconstructing high-quality GEMs and elaborate on the resources needed to advance a system biology approach aimed at untangling the intimate interactions between plants and pathogens. [ABSTRACT FROM AUTHOR]

Metabolic conversions allow organisms to produce essential metabolites from the available nutrients in an environment, frequently requiring metabolic exchanges among co-inhabiting organisms. Here, we applied genomic-based simulations for exploring tri-trophic interactions among the sap-feeding insect whitefly (Bemisia tabaci), its host-plants, and symbiotic bacteria. The simplicity of this ecosystem allows capturing the interacting organisms (based on genomic data) and the environmental content (based on metabolomics data). Simulations explored the metabolic capacities of insect-symbiont combinations under environments representing natural phloem. Predictions were correlated with experimental data on the dynamics of symbionts under different diets. Simulation outcomes depict a puzzle of three-layer origins (plant-insect-symbionts) for the source of essential metabolites across habitats and stratify interactions enabling the whitefly to feed on diverse hosts. In parallel to simulations, natural and artificial feeding experiments provide supporting evidence for an environment-based effect on symbiont dynamics. Based on simulations, a decrease in the relative abundance of a symbiont can be associated with a loss of fitness advantage due to an environmental excess in amino-acids whose production in a deprived environment used to depend on the symbiont. The study demonstrates that genomic-based predictions can bridge environment and community dynamics and guide the design of symbiont manipulation strategies. [ABSTRACT FROM AUTHOR]

Pregnancy generally alters the balance of maternal metabolism, but the molecular profiles in early pregnancy and associated factors of folate supplementation in pregnant women remains incompletely understood. Untargeted metabonomics based on high-performance liquid chromatography-high-resolution mass spectrometry integrated with multivariate metabolic pathway analysis were applied to characterize metabolite profiles and associated factors of folate supplements in early pregnancy. The metabolic baseline of early pregnancy was determined by metabolic analysis of 510 serum samples from 131 non-pregnant and 379 pregnant healthy Chinese women. The pathophysiology of adaptive reactions and metabolic challenges induced by folate supplementation in early pregnancy was further compared between pregnant women with (n = 168) and without (n = 184) folate supplements. Compared with non-pregnant participants, 106 metabolites, majority of which are related to amino acids and lysophosphatidylcholine/phosphatidylcholine, and 13 metabolic pathways were significantly changed in early pregnancy. The supplementation of folate in early pregnancy induced marked changes in N-acyl ethanolamine 22:0, N-acyl taurine 18:2, glycerophosphoserine 44:1 and 8,11,14-eicosatrienoate, proline, and aminoimidazole ribotide levels. During early pregnancy, the metabolism of amino acids significantly changes to meet the physiological requirements of pregnant women. Folate intake may change glucose and lipid metabolism. These findings provide a comprehensive landscape for understanding the basic characteristics and gestational metabolic networks of early pregnancy and folate supplementation. This study provides a basis for further research into the relationship between metabolic markers and pregnancy diseases. This study protocol was registered on www.ClinicalTrials.gov , NCT03651934, on August 29, 2018 (prior to recruitment). [ABSTRACT FROM AUTHOR]

URAT1 (slc22a12) was identified as the transporter responsible for renal reabsorption of the medically important compound, uric acid. However, subsequent studies have indicated that other transporters make contributions to this process, and that URAT1 transports other organic anions besides urate (including several in common with the closely related multi-specific renal organic anion transporters, OAT1 (slc22a6) and OAT3 (slc22a8)). These findings raise the possibility that urate transport is not the sole physiological function of URAT1. We previously characterized mice null for the murine ortholog of URAT1 (mURAT1; previously cloned as RST), finding a relatively modest decrement in urate reabsorptive capacity. Nevertheless, there were shifts in the plasma and urinary concentrations of multiple small molecules, suggesting significant metabolic changes in the knockouts. Although these molecules remain unidentified, here we have computationally delineated the biochemical networks consistent with transcriptomic data from the null mice. These analyses suggest alterations in the handling of not only urate but also other putative URAT1 substrates comprising intermediates in nucleotide, carbohydrate, and steroid metabolism. Moreover, the analyses indicate changes in multiple other pathways, including those relating to the metabolism of glycosaminoglycans, methionine, and coenzyme A, possibly reflecting downstream effects of URAT1 loss. Taken together with the available substrate and metabolomic data for the other OATs, our findings suggest that the transport and biochemical functions of URAT1 overlap those of OAT1 and OAT3, and could contribute to our understanding of the relationship between uric acid and the various metabolic disorders to which it has been linked.

High fructose consumption has been implicated with deleterious effects on human health, including hyperlipidemia elicited through de novo lipogenesis. However, more global effects of fructose on cellular metabolism have not been elucidated. In order to explore the metabolic impact of fructose-containing nutrients, we applied both GC-TOF and HILIC-QTOF mass spectrometry metabolomic strategies using extracts from cultured HepG2 cells exposed to fructose, glucose, or fructose + glucose. Cellular responses were analyzed in a time-dependent manner, incubated in media containing 5.5 mM glucose + 5.0 mM fructose in comparison to controls incubated in media containing either 5.5 mM glucose or 10.5 mM glucose. Mass spectrometry identified 156 unique known metabolites and a large number of unknown compounds, which revealed metabolite changes due to both utilization of fructose and high-carbohydrate loads independent of hexose structure. Fructose was shown to be partially converted to sorbitol, and generated higher levels of fructose-1-phosphate as a precursor for glycolytic intermediates. Differentially regulated ratios of 3-phosphoglycerate to serine pathway intermediates in high fructose media indicated a diversion of carbon backbones away from energy metabolism. Additionally, high fructose conditions changed levels of complex lipids toward phosphatidylethanolamines. Patterns of acylcarnitines in response to high hexose exposure (10.5 mM glucose or glucose/fructose combination) suggested a reduction in mitochondrial beta-oxidation.