Search

Metabolic reprogramming is a general feature of cancer cells. Regrettably, the comprehensive quantification of metabolites in biological specimens does not promptly translate into knowledge on the utilization of metabolic pathways. By estimating fluxes across metabolic pathways, computational models hold the promise to bridge this gap between data and biological functionality. These models currently portray the average behavior of cell populations however, masking the inherent heterogeneity that is part and parcel of tumorigenesis as much as drug resistance. To remove this limitation, we propose single-cell Flux Balance Analysis (scFBA) as a computational framework to translate single-cell transcriptomes into single-cell fluxomes. We show that the integration of single-cell RNA-seq profiles of cells derived from lung adenocarcinoma and breast cancer patients into a multi-scale stoichiometric model of a cancer cell population: significantly 1) reduces the space of feasible single-cell fluxomes; 2) allows to identify clusters of cells with different growth rates within the population; 3) points out the possible metabolic interactions among cells via exchange of metabolites. The scFBA suite of MATLAB functions is available at , as well as the case study datasets. [ABSTRACT FROM AUTHOR]

Continuous cultures of mammalian cells are complex systems displaying hallmark phenomena of nonlinear dynamics, such as multi-stability, hysteresis, as well as sharp transitions between different metabolic states. In this context mathematical models may suggest control strategies to steer the system towards desired states. Although even clonal populations are known to exhibit cell-to-cell variability, most of the currently studied models assume that the population is homogeneous. To overcome this limitation, we use the maximum entropy principle to model the phenotypic distribution of cells in a chemostat as a function of the dilution rate. We consider the coupling between cell metabolism and extracellular variables describing the state of the bioreactor and take into account the impact of toxic byproduct accumulation on cell viability. We present a formal solution for the stationary state of the chemostat and show how to apply it in two examples. First, a simplified model of cell metabolism where the exact solution is tractable, and then a genome-scale metabolic network of the Chinese hamster ovary (CHO) cell line. Along the way we discuss several consequences of heterogeneity, such as: qualitative changes in the dynamical landscape of the system, increasing concentrations of byproducts that vanish in the homogeneous case, and larger population sizes. [ABSTRACT FROM AUTHOR]

Today we live in a world of digital objects and digital technology; industry and humanities as well as technologies are truly in the midst of a digital environment driven by ICT and cyber informatics. A digital ecosystem can be defined as a digital environment populated by interacting and competing digital species. Digital species have autonomous, proactive and adaptive behaviors, regulated by peer-to-peer interactions without central control point. An interconnecting architecture with few highly connected nodes (hubs) and many low connected nodes has a scale- free architecture. A new bio-inspired analysis methodology (BIAM) environment, an investigation strategy for information flow, fault and error tolerance detection in digital ecosystems based on a scale-free architecture is presented in this paper. In order to extract the information about modules and digital species role, the analysis methodology, inspired by metabolic network working, implements a set of three interacting techniques, i.e., topological analysis, flux balance analysis and extreme pathway analysis. Highly connected nodes, intermodule connectors and ultra-peripheral nodes can be identified by evaluating their impact on digital ecosystems behavior and addressing their strengthen, fault tolerance and protection countermeasures. Two real case studies of ecosystems have been analyzed in order to test the functionalities of the proposed (BIAM) environment and the goodness of this approach. [ABSTRACT FROM AUTHOR]

Cellular systems shift metabolic states by adjusting gene expression and enzyme activities to adapt to physiological and environmental changes. Biochemical and genetic studies are identifying how metabolic regulation affects the selection of metabolic phenotypes. However, how metabolism influences its regulatory architecture still remains unexplored. We present a new method of extreme pathway analysis (the minimal set of conically independent metabolic pathways) to deduce regulatory structures from pure pathway information. Applying our method to metabolic networks of human red blood cells and Escherichia coli, we shed light on how metabolic regulation are organized by showing which reactions within metabolic networks are more prone to transcriptional or allosteric regulation. Applied to a human genome-scale metabolic system, our method detects disease-associated reactions. Thus, our study deepens the understanding of the organizing principle of cellular metabolic regulation and may contribute to metabolic engineering, synthetic biology, and disease treatment. [ABSTRACT FROM AUTHOR]

Constraint-based modeling (CBM) is increasingly used to analyze the metabolism of complex microbial communities involved in ecology, biomedicine, and various biotechnological processes. While CBM is an established framework for studying the metabolism of single species with linear stoichiometric models, CBM of communities with balanced growth is more complicated, not only due to the larger size of the multi-species metabolic network but also because of the bilinear nature of the resulting community models. Moreover, the solution space of these community models often contains biologically unrealistic solutions, which, even with model linearization and under application of certain objective functions, cannot easily be excluded. Here we present RedCom, a new approach to build reduced community models in which the metabolisms of the participating organisms are represented by net conversions computed from the respective single-species networks. By discarding (single-species) net conversions that violate a minimality criterion in the exchange fluxes, it is ensured that unrealistic solutions in the community model are excluded where a species altruistically synthesizes large amounts of byproducts (instead of biomass) to fulfill the requirements of other species. We employed the RedCom approach for modeling communities of up to nine organisms involved in typical degradation steps of anaerobic digestion in biogas plants. Compared to full (bilinear and linearized) community models, we found that the reduced community models obtained with RedCom are not only much smaller but allow, also in the largest model with nine species, extensive calculations required to fully characterize the solutions space and to reveal key properties of communities with maximum methane yield and production rates. Furthermore, the predictive power of the reduced community models is significantly larger because they predict much smaller ranges of feasible community compositions and exchange fluxes still being consistent with measurements obtained from enrichment cultures. For an enrichment culture for growth on ethanol, we also used metaproteomic data to further constrain the solution space of the community models. Both model and proteomic data indicated a dominance of acetoclastic methanogens (Methanosarcinales) and Desulfovibrionales being the least abundant group in this microbial community. [ABSTRACT FROM AUTHOR]

Gut microbiota and human relationships are strictly connected to each other. What we eat reflects our body-mind connection and synchronizes with people around us. However, how this impacts on gut microbiota and, conversely, how gut bacteria influence our dietary behaviors has not been explored yet. To quantify the complex dynamics of this interplay between gut and human behaviors we explore the “gut-human behavior axis” and its evolutionary dynamics in a real-world scenario represented by the social multiplex network. We consider a dual type of similarity, homophily and gut similarity, other than psychological and unconscious biases. We analyze the dynamics of social and gut microbial communities, quantifying the impact of human behaviors on diets and gut microbial composition and, backwards, through a control mechanism. Meal timing mechanisms and “chrono-nutrition” play a crucial role in feeding behaviors, along with the quality and quantity of food intake. Considering a population of shift workers, we explore the dynamic interplay between their eating behaviors and gut microbiota, modeling the social dynamics of chrono-nutrition in a multiplex network. Our findings allow us to quantify the relation between human behaviors and gut microbiota through the methodological introduction of gut metabolic modeling and statistical estimators, able to capture their dynamic interplay. Moreover, we find that the timing of gut microbial communities is slower than social interactions and shift-working, and the impact of shift-working on the dynamics of chrono-nutrition is a fluctuation of strategies with a major propensity for defection (e.g. high-fat meals). A deeper understanding of the relation between gut microbiota and the dietary behavioral patterns, by embedding also the related social aspects, allows improving the overall knowledge about metabolic models and their implications for human health, opening the possibility to design promising social therapeutic dietary interventions. [ABSTRACT FROM AUTHOR]

Cellular functions are shaped by reaction networks whose dynamics are determined by the concentrations of underlying components. However, cellular mechanisms ensuring that a component’s concentration resides in a given range remain elusive. We present network properties which suffice to identify components whose concentration ranges can be efficiently computed in mass-action metabolic networks. We show that the derived ranges are in excellent agreement with simulations from a detailed kinetic metabolic model of Escherichia coli. We demonstrate that the approach can be used with genome-scale metabolic models to arrive at predictions concordant with measurements from Escherichia coli under different growth scenarios. By application to 14 genome-scale metabolic models from diverse species, our approach specifies the cellular determinants of concentration ranges that can be effectively employed to make predictions for a variety of biotechnological and medical applications. [ABSTRACT FROM AUTHOR]

Growth rate is one of the most important and most complex phenotypic characteristics of unicellular microorganisms, which determines the genetic mutations that dominate at the population level, and ultimately whether the population will survive. Translating changes at the genetic level to their growth-rate consequences remains a subject of intense interest, since such a mapping could rationally direct experiments to optimize antibiotic efficacy or bioreactor productivity. In this work, we directly map transcriptional profiles to growth rates by gathering published gene-expression data from Escherichia coli and Saccharomyces cerevisiae with corresponding growth-rate measurements. Using a machine-learning technique called k-nearest-neighbors regression, we build a model which predicts growth rate from gene expression. By exploiting the correlated nature of gene expression and sparsifying the model, we capture 81% of the variance in growth rate of the E. coli dataset, while reducing the number of features from >4,000 to 9. In S. cerevisiae, we account for 89% of the variance in growth rate, while reducing from >5,500 dimensions to 18. Such a model provides a basis for selecting successful strategies from among the combinatorial number of experimental possibilities when attempting to optimize complex phenotypic traits like growth rate. [ABSTRACT FROM AUTHOR]

Metabolic networks comprise thousands of enzymatic reactions functioning in a controlled manner and have been shaped by natural selection. Thanks to the genome data, the footprints of adaptive (positive) selection are detectable, and the strength of purifying selection can be measured. This has made possible to know where, in the metabolic network, adaptive selection has acted and where purifying selection is more or less strong and efficient. We have carried out a comprehensive molecular evolutionary study of all the genes involved in the human metabolism. We investigated the type and strength of the selective pressures that acted on the enzyme-coding genes belonging to metabolic pathways during the divergence of primates and rodents. Then, we related those selective pressures to the functional and topological characteristics of the pathways. We have used DNA sequences of all enzymes (956) of the metabolic pathways comprised in the HumanCyc database, using genome data for humans and five other mammalian species. We have found that the evolution of metabolic genes is primarily constrained by the layer of the metabolism in which the genes participate: while genes encoding enzymes of the inner core of metabolism are much conserved, those encoding enzymes participating in the outer layer, mediating the interaction with the environment, are evolutionarily less constrained and more plastic, having experienced faster functional evolution. Genes that have been targeted by adaptive selection are endowed by higher out-degree centralities than non-adaptive genes, while genes with high in-degree centralities are under stronger purifying selection. When the position along the pathway is considered, a funnel-like distribution of the strength of the purifying selection is found. Genes at bottom positions are highly preserved by purifying selection, whereas genes at top positions, catalyzing the first steps, are open to evolutionary changes. These results show how functional and topological characteristics of metabolic pathways contribute to shape the patterns of evolutionary pressures driven by natural selection and how pathway network structure matters in the evolutionary process that shapes the evolution of the system. [ABSTRACT FROM AUTHOR]

Metabolic biomarkers may play an important role in the diagnosis, prognostication and assessment of response to pharmacological therapy in complex diseases. The process of discovering new metabolic biomarkers is a non-trivial task which involves a number of bioanalytical processing steps coupled with a computational approach for the search, prioritization and verification of new biomarker candidates. Kinetic analysis provides an additional dimension of complexity in time-series data, allowing for a more precise interpretation of biomarker dynamics in terms of molecular interaction and pathway modulation. A novel network-based computational strategy for the discovery of putative dynamic biomarker candidates is presented, enabling the identification and verification of unexpected metabolic signatures in complex diseases such as myocardial infarction. The novelty of the proposed method lies in combining metabolic time-series data into a superimposed graph representation, highlighting the strength of the underlying kinetic interaction of preselected analytes. Using this approach, we were able to confirm known metabolic signatures and also identify new candidates such as carnosine and glycocholic acid, and pathways that have been previously associated with cardiovascular or related diseases. This computational strategy may serve as a complementary tool for the discovery of dynamic metabolic or proteomic biomarkers in the field of clinical medicine. [ABSTRACT FROM AUTHOR]

Defining the essential gene components for a system to be considered alive is a crucial step toward the synthesis of artificial life. Fifteen years ago, Gil and coworkers proposed the core of a putative minimal bacterial genome, which would provide the capability to achieve metabolic homeostasis, reproduce, and evolve to a bacterium in an ideally controlled environment. They also proposed a simplified metabolic chart capable of providing energy and basic components for a minimal living cell. For this work, we have identified the components of the minimal metabolic network based on the aforementioned studies, associated them to the KEGG database and, by applying the MetaDAG methodology, determined its Metabolic Building Blocks (MBB) and reconstructed its metabolic Directed Acyclic Graph (m-DAG). The reaction graph of this metabolic network consists of 80 compounds and 98 reactions, while its m-DAG has 36 MBBs. Additionally, we identified 12 essential reactions in the m-DAG that are critical for maintaining the connectivity of this network. In a similar manner, we reconstructed the m-DAG of JCVI-syn3.0, which is an artificially designed and manufactured viable cell whose genome arose by minimizing the one from Mycoplasma mycoides JCVI-syn1.0, and of “Candidatus Nasuia deltocephalinicola”, the bacteria with the smallest natural genome known to date. The comparison of the m-DAGs derived from a theoretical, an artificial, and a natural genome denote slightly different lifestyles, with a consistent core metabolism. The MetaDAG methodology we employ uses homogeneous descriptors and identifiers from the KEGG database, so that comparisons between bacterial strains are not only easy but also suitable for many research fields. The modeling of m-DAGs based on minimal metabolisms can be the first step for the synthesis and manipulation of minimal cells.

Structural analysis of constraint-based metabolic network models attempts to find the network’s properties by searching for subsets of suitable modes or Elementary Flux Modes (EFMs). One useful approach is based on Linear Program (LP) techniques, which introduce an objective function to convert the stoichiometric and thermodynamic constraints into a linear program (LP), using additional constraints to generate different nontrivial modes. This work introduces FLFS-FC (Fixed Length Function Sampling with Flux Coupling), a new approach to increase the efficiency of generation of large sets of different EFMs for the network. FLFS-FC is based on the importance of the length of the objective functions used in the associated LP problem and the imposition of additional negative constraints. Our proposal overrides some of the known drawbacks associated with the EFM extraction, such as the appearance of unfeasible problems or multiple repeated solutions arising from different LP problems.

Metabolic conversions allow organisms to produce a set of essential metabolites from the available nutrients in an environment, frequently requiring metabolic exchanges among co-inhabiting organisms. Genomic-based metabolic simulations are being increasingly applied for exploring metabolic capacities, considering different environments and different combinations of microorganisms. NetMet is a web-based tool and a software package for predicting the metabolic performances of microorganisms and their corresponding combinations in user-defined environments. The algorithm takes, as input, lists of (i) species-specific enzymatic reactions (EC numbers), and (ii) relevant metabolic environments. The algorithm generates, as output, lists of (i) compounds that individual species can produce in each given environment, and (ii) compounds that are predicted to be produced through complementary interactions. The tool is demonstrated in two case studies. First, we compared the metabolic capacities of different haplotypes of the obligatory fruit and vegetable pathogen Candidatus Liberibacter solanacearum to those of their culturable taxonomic relative Liberibacter crescens. Second, we demonstrated the potential production of complementary metabolites by pairwise combinations of co-occurring endosymbionts of the plant phloem-feeding whitefly Bemisia tabaci.

The integration of high-throughput data to build predictive computational models of cellular metabolism is a major challenge of systems biology. These models are needed to predict cellular responses to genetic and environmental perturbations. Typically, this response involves both metabolic regulations related to the kinetic properties of enzymes and a genetic regulation affecting their concentrations. Thus, the integration of the transcriptional regulatory information is required to improve the accuracy and predictive ability of metabolic models. Integrative modeling is of primary importance to guide the search for various applications such as discovering novel potential drug targets to develop efficient therapeutic strategies for various diseases. In this paper, we propose an integrative predictive model based on techniques combining semantic web, probabilistic modeling, and constraint-based modeling methods. We applied our approach to human cancer metabolism to predict in silico the growth response of specific cancer cells under approved drug effects. Our method has proven successful in predicting the biomass rates of human liver cancer cells under drug-induced transcriptional perturbations.

The modelbase package is a free expandable Python package for building and analysing dynamic mathematical models of biological systems. Originally it was designed for the simulation of metabolic systems, but it can be used for virtually any deterministic chemical processes. modelbase provides easy construction methods to define reactions and their rates. Based on the rates and stoichiometries, the system of differential equations is assembled automatically. modelbase minimises the constraints imposed on the user, allowing for easy and dynamic access to all variables, including derived ones, in a convenient manner. A simple incorporation of algebraic equations is, for example, convenient to study systems with rapid equilibrium or quasi steady-state approximations. Moreover, modelbase provides construction methods that automatically build all isotope-specific versions of a particular reaction, making it a convenient tool to analyse non-steady state isotope-labelling experiments. Funding statement: This work was financially supported by the Deutsche Forschungsgemeinschaft “Cluster of Excellence on Plant Sciences” CEPLAS (EXC 1028).

Determining how microbial communities organize and function at the ecosystem level is essential to understanding and predicting how they will respond to environmental change. Mathematical models can be used to describe these communities, but properly representing all the biological interactions in extremely diverse natural microbial ecosystems in a mathematical model is challenging. We examine a complementary approach based on the maximum entropy production (MEP) principle, which proposes that systems with many degrees of freedom will likely organize to maximize the rate of free energy dissipation. In this study, we develop an MEP model to describe biogeochemistry observed in Siders Pond, a phosphate limited meromictic system located in Falmouth, MA that exhibits steep chemical gradients due to density-driven stratification that supports anaerobic photosynthesis as well as microbial communities that catalyze redox cycles involving O, N, S, Fe, and Mn. The MEP model uses a metabolic network to represent microbial redox reactions, where biomass allocation and reaction rates are determined by solving an optimization problem that maximizes entropy production over time, and a 1D vertical profile constrained by an advection-dispersion-reaction model. We introduce a new approach for modeling phototrophy and explicitly represent oxygenic photoautotrophs, photoheterotrophs and anoxygenic photoautotrophs. The metabolic network also includes reactions for aerobic organoheterotrophic bacteria, sulfate reducing bacteria, sulfide oxidizing bacteria and aerobic and anaerobic grazers. Model results were compared to observations of biogeochemical constituents collected over a 24 h period at 8 depths at a single 15 m deep station in Siders Pond. Maximizing entropy production over long (3 day) intervals produced results more similar to field observations than short (0.25 day) interval optimizations, which support the importance of temporal strategies for maximizing entropy production over time. Furthermore, we found that entropy production must be maximized locally instead of globally where energy potentials are degraded quickly by abiotic processes, such as light absorption by water. This combination of field observations and modeling results indicate that natural microbial systems can be modeled by using the maximum entropy production principle applied over time and space using many fewer parameters than conventional models.

ABSTRACT 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.

In Salmonella enterica, 2-aminoacrylate (2AA) is a reactive enamine intermediate generated during a number of biochemical reactions. When the 2-iminobutanoate/2-iminopropanoate deaminase (RidA; EC: 3.5.99.10) is eliminated, 2AA accumulates and inhibits the activity of multiple pyridoxal 5’-phosphate(PLP)-dependent enzymes. In this study, untargeted proton nuclear magnetic resonance (1H NMR) metabolomics and transcriptomics data were used to uncover the global metabolic response of S. enterica to the accumulation of 2AA. The data showed that elimination of RidA perturbed folate and branched chain amino acid metabolism. Many of the resulting perturbations were consistent with the known effect of 2AA stress, while other results suggested additional potential enzyme targets of 2AA-dependent damage. The majority of transcriptional and metabolic changes appeared to be the consequence of downstream effects on the metabolic network, since they were not directly attributable to a PLP-dependent enzyme. In total, the results highlighted the complexity of changes stemming from multiple perturbations of the metabolic network, and suggested hypotheses that will be valuable in future studies of the RidA paradigm of endogenous 2AA stress.

The identification of therapeutic strategies exploiting the metabolic alterations of malignant cells is a relevant area in cancer research. Here, we discuss a novel computational method, based on the COBRA (COnstraint-Based Reconstruction and Analysis) framework for metabolic networks, to perform this task. Current and future steps are presented.

Host cells provide protection from the environment for intracellular pathogens, but acquisition of nutrients and sensing of changes in the host cell bring new challenges. Krishnan et al. have recently developed a comprehensive genome-scale metabolic model (GEM) of the Toxoplasma gondii metabolic network that incorporates genetic, transcriptomic, and metabolomic data. [ABSTRACT FROM AUTHOR]

Subnetwork mining is an essential issue in the analysis of biological, social and communication networks. Recent applications require the simultaneous mining of several networks on the same or a similar vertex set. That is, one searches for subnetworks fulfilling different properties in each input network. We study the case that the input consists of a directed graph D and an undirected graph G on the same vertex set, and the sought pattern is a path P in D whose vertex set induces a connected subgraph of G. In this context, three concrete problems arise, depending on whether the existence of P is questioned or whether the length of P is to be optimized: in that case, one can search for a longest path or (maybe less intuitively) a shortest one. These problems have immediate applications in biological networks and predictable applications in social, information and communication networks. We study the classic and parameterized complexity of the problem, thus identifying polynomial and NP-complete cases, as well as fixed-parameter tractable and W[1]-hard cases. We also propose two enumeration algorithms that we evaluate on synthetic and biological data.

Females and males are known to have different abilities to cope with stress and disease. This study was designed to investigate the effect of sex on properties of a complex interlinked network constructed of central biochemical metabolites. The study involved the blood collection and analysis of a large set of blood metabolic markers from a total of 236 healthy participants, which included 140 females and 96 males. Metabolic profiling yielded concentrations of 168 metabolites for each subject. A differential correlation network analysis approach was developed for this study that allowed detection and characterization of interconnection differences in metabolites in males and females. Through topological analysis of the differential network that depicted metabolite differences in the sexes, we identified metabolites with high centralities in this network. These key metabolites were identified as 10 phosphatidylcholines (PCaaC34:4, PCaaC36:6, PCaaC34:3, PCaaC42:2, PCaeC38:1, PCaeC38:2, PCaaC40:1, PCaeC34:1, PC aa C32:1 and PC aa C40:6) and 4 acylcarnitines (C3-OH, C7-DC, C3 and C0). Identification of these metabolites may help further studies of sex-specific differences in the metabolome that may underlie different responses to stress and disease in males and females. [ABSTRACT FROM AUTHOR]

Many link prediction methods have been developed to infer unobserved links or predict missing links based on the observed network structure that is always incomplete and subject to interfering noise. Thus, the performance of existing methods is usually limited in that their computation depends only on input graph structures, and they do not consider external information. The effects of social influence and homophily suggest that both network structure and node attribute information should help to resolve the task of link prediction. This work proposes SASNMF, a link prediction unified framework based on non-negative matrix factorization that considers not only graph structure but also the internal and external auxiliary information, which refers to both the node attributes and the structural latent feature information extracted from the network. Furthermore, three different combinations of internal and external information are proposed and input into the framework to solve the link prediction problem. Extensive experimental results on thirteen real networks, five node attribute networks and eight non-attribute networks show that the proposed framework has competitive performance compared with benchmark methods and state-of-the-art methods, indicating the superiority of the presented algorithm. [ABSTRACT FROM AUTHOR]

Essential metabolic reactions are shaping constituents of metabolic networks, enabling viable and distinct phenotypes across diverse life forms. Here we analyse and compare modelling predictions of essential metabolic functions with experimental data and thereby identify core metabolic pathways in prokaryotes. Simulations of 15 manually curated genome-scale metabolic models were integrated with 36 large-scale gene essentiality datasets encompassing a wide variety of species of bacteria and archaea. Conservation of metabolic genes was estimated by analysing 79 representative genomes from all the branches of the prokaryotic tree of life. We find that essentiality patterns reflect phylogenetic relations both for modelling and experimental data, which correlate highly at the pathway level. Genes that are essential for several species tend to be highly conserved as opposed to non-essential genes which may be conserved or not. The tRNA-charging module is highlighted as ancestral and with high centrality in the networks, followed closely by cofactor metabolism, pointing to an early information processing system supplied by organic cofactors. The results, which point to model improvements and also indicate faults in the experimental data, should be relevant to the study of centrality in metabolic networks and ancient metabolism but also to metabolic engineering with prokaryotes. [ABSTRACT FROM AUTHOR]

Carbon dioxide enters metabolism via six known CO2 fixation pathways, of which only one is linear, exergonic in the direction of CO2‐assimilation, and present in both bacterial and archaeal anaerobes – the Wood‐Ljungdahl (WL) or reductive acetyl‐CoA pathway. Carbon monoxide (CO) plays a central role in the WL pathway as an energy rich intermediate. Here, we scan the major biochemical reaction databases for reactions involving CO and CO2. We identified 415 reactions corresponding to enzyme commission (EC) numbers involving CO2, which are non‐randomly distributed across different biochemical pathways. Their taxonomic distribution, reversibility under physiological conditions, cofactors and prosthetic groups are summarized. In contrast to CO2, only 15 reaction classes involving CO were detected. Closer inspection reveals that CO interfaces with metabolism and the carbon cycle at only two enzymes: anaerobic carbon monoxide dehydrogenase (CODH), a Ni‐ and Fe‐containing enzyme that generates CO for CO2 fixation in the WL pathway, and aerobic CODH, a Mo‐ and Cu‐containing enzyme that oxidizes environmental CO as an electron source. The CO‐dependent reaction of the WL pathway involves carbonyl insertion into a methyl carbon‐nickel at the Ni‐Fe‐S A‐cluster of acetyl‐CoA synthase (ACS). It appears that no alternative mechanisms to the CO‐dependent reaction of ACS have evolved in nearly 4 billion years, indicating an ancient and mechanistically essential role for CO at the onset of metabolism. Carbon monoxide (CO) is essential for CO2 fixation in many microbes that inhabit strictly anaerobic environments. Here, we scan major biochemical databases for the involvement of CO and CO2 in metabolism. While CO2 participates in hundreds of reactions, CO enters metabolism at only two. CO2 assimilation via CO is one of biochemistry's most ancient and conserved reactions. CO is apparently special. [ABSTRACT FROM AUTHOR]

Evidence from the literature keeps highlighting the impact of mutualistic bacterial communities of the gut microbiota on human health. The gut microbita is a complex ecosystem of symbiotic bacteria which contributes to mammalian host biology by processing, otherwise, indigestible nutrients, supplying essential metabolites, and contributing to modulate its immune system. Advances in sequencing technologies have enabled structural analysis of the human gut microbiota and allowed detection of changes in gut bacterial composition in several common diseases, including cardiometabolic disorders. Biological signals sent by the gut microbiota to the host, including microbial metabolites and pro-inflammatory molecules, mediate microbiome-host genome cross-talk. This rapidly expanding line of research can identify disease-causing and disease-predictive microbial metabolite biomarkers, which can be translated into novel biodiagnostic tests, dietary supplements, and nutritional interventions for personalized therapeutic developments in common diseases. Here, we review results from the most significant studies dealing with the association of products from the gut microbial metabolism with cardiometabolic disorders. We underline the importance of these postbiotic biomarkers in the diagnosis and treatment of human disorders. [ABSTRACT FROM AUTHOR]

Effective stratification of cancer patients on the basis of their molecular make-up is a key open challenge. Given the altered and heterogenous nature of cancer metabolism, we here propose to use the overall expression of central carbon metabolism as biomarker to characterize groups of patients with important characteristics, such as response to ad-hoc therapeutic strategies and survival expectancy. To this end, we here introduce the data integration framework named Metabolic Reaction Enrichment Analysis (MaREA), which strives to characterize the metabolic deregulations that distinguish cancer phenotypes, by projecting RNA-seq data onto metabolic networks, without requiring metabolic measurements. MaREA computes a score for each network reaction, based on the expression of the set of genes encoding for the associated enzyme(s). The scores are first used as features for cluster analysis and then to rank and visualize in an organized fashion the metabolic deregulations that distinguish cancer sub-types. We applied our method to recent lung and breast cancer RNA-seq datasets from The Cancer Genome Atlas and we were able to identify subgroups of patients with significant differences in survival expectancy. We show how the prognostic power of MaREA improves when an extracted and further curated core model focusing on central carbon metabolism is used rather than the genome-wide reference network. The visualization of the metabolic differences between the groups with best and worst prognosis allowed to identify and analyze key metabolic properties related to cancer aggressiveness. Some of these properties are shared across different cancer (sub) types, e.g., the up-regulation of nucleic acid and amino acid synthesis, whereas some other appear to be tumor-specific, such as the up- or down-regulation of the phosphoenolpyruvate carboxykinase reaction, which display different patterns in distinct tumor (sub)types. These results might be soon employed to deliver highly automated diagnostic and prognostic strategies for cancer patients. [ABSTRACT FROM AUTHOR]

Over the last years, several genome‐scale metabolic models (GEMs) and kinetic models of Escherichia coli were published. Their predictive performance varies according to the evaluation metric considered, the computational simulation methods used, and the type/quality of experimental data available. However, the GEM approach is often not compared with the kinetic modeling framework. Also, the different genome‐scale reconstruction versions and simulation methods of mutant phenotypes are usually not validated to predict intracellular fluxes using large experimental datasets. Here, we intended to (i) systematically evaluate the prediction performance of three E. coli GEMs (iJR904, iAF1260, and iJO1366) available in the literature according to predictive growth metrics (intracellular flux distribution); (ii) assess the reliability of a E. coli GEM in the prediction of gene knockout phenotypes when different simulation methods (parsimonious flux balance analysis, Minimization of Metabolic Adjustment, linear version of MoMA, Regulatory on/off minimization, and Minimization of Metabolites Balance) are used; and finally (iii) investigate the flux distribution predictive power of the constrained‐based modeling approach (selected stoichiometric GEM) and compare it with the kinetic modeling approach (two published kinetic models) for E. coli central metabolism, in order to assess their accuracy. Results show that the phenotype predictions were not significantly sensitive to the metabolic models, although the GEM iAF1260 was more accurate in the prediction of central carbon fluxes at low dilution rates. Furthermore, we observed that the choice of the appropriate simulation method of mutant phenotypes depends on the biological question to be addressed. In terms of the two modeling approaches, none outperformed the other for all the tested scenarios. © 2018 American Institute of Chemical Engineers Biotechnol. Prog., 34:1344–1354, 2018 [ABSTRACT FROM AUTHOR]

Even across genomes of the same species, prokaryotes exhibit remarkable flexibility in gene content. We do not know whether this flexible or “accessory” content is mostly neutral or adaptive, largely due to the lack of explicit analyses of accessory gene function. Here, across 96 diverse prokaryotic species, I show that a considerable fraction (~40%) of accessory genomes harbours beneficial metabolic functions. These functions take two forms: (1) they significantly expand the biosynthetic potential of individual strains, and (2) they help reduce strain-specific metabolic auxotrophies via intra-species metabolic exchanges. I find that the potential of both these functions increases with increasing genome flexibility. Together, these results are consistent with a significant adaptive role for prokaryotic pangenomes. [ABSTRACT FROM AUTHOR]

RAVEN is a commonly used MATLAB toolbox for genome-scale metabolic model (GEM) reconstruction, curation and constraint-based modelling and simulation. Here we present RAVEN Toolbox 2.0 with major enhancements, including: (i) de novo reconstruction of GEMs based on the MetaCyc pathway database; (ii) a redesigned KEGG-based reconstruction pipeline; (iii) convergence of reconstructions from various sources; (iv) improved performance, usability, and compatibility with the COBRA Toolbox. Capabilities of RAVEN 2.0 are here illustrated through de novo reconstruction of GEMs for the antibiotic-producing bacterium Streptomyces coelicolor. Comparison of the automated de novo reconstructions with the iMK1208 model, a previously published high-quality S. coelicolor GEM, exemplifies that RAVEN 2.0 can capture most of the manually curated model. The generated de novo reconstruction is subsequently used to curate iMK1208 resulting in Sco4, the most comprehensive GEM of S. coelicolor, with increased coverage of both primary and secondary metabolism. This increased coverage allows the use of Sco4 to predict novel genome editing targets for optimized secondary metabolites production. As such, we demonstrate that RAVEN 2.0 can be used not only for de novo GEM reconstruction, but also for curating existing models based on up-to-date databases. Both RAVEN 2.0 and Sco4 are distributed through GitHub to facilitate usage and further development by the community ( and ). [ABSTRACT FROM AUTHOR]

The genetic regulatory network (GRN) plays a key role in controlling the response of the cell to changes in the environment. Although the structure of GRNs has been the subject of many studies, their large scale structure in the light of feedbacks from the metabolic network (MN) has received relatively little attention. Here we study the causal structure of the GRNs, namely the chain of influence of one component on the other, taking into account feedback from the MN. First we consider the GRNs of E. coli and B. subtilis without feedback from MN and illustrate their causal structure. Next we augment the GRNs with feedback from their respective MNs by including (a) links from genes coding for enzymes to metabolites produced or consumed in reactions catalyzed by those enzymes and (b) links from metabolites to genes coding for transcription factors whose transcriptional activity the metabolites alter by binding to them. We find that the inclusion of feedback from MN into GRN significantly affects its causal structure, in particular the number of levels and relative positions of nodes in the hierarchy, and the number and size of the strongly connected components (SCCs). We then study the functional significance of the SCCs. For this we identify condition specific feedbacks from the MN into the GRN by retaining only those enzymes that are essential for growth in specific environmental conditions simulated via the technique of flux balance analysis (FBA). We find that the SCCs of the GRN augmented by these feedbacks can be ascribed specific functional roles in the organism. Our algorithmic approach thus reveals relatively autonomous subsystems with specific functionality, or regulatory modules in the organism. This automated approach could be useful in identifying biologically relevant modules in other organisms for which network data is available, but whose biology is less well studied. [ABSTRACT FROM AUTHOR]

The aim of this study was to investigate the physiological and proteomic changes in the embryo and endosperm during seed germination in the elite Chinese bread wheat cultivar Zhengmai 366. Upon imbibition, seed size and water content increased rapidly, followed by a series of metabolic changes including increases in soluble sugar content and α-amylase activity, a decrease in starch content, and a rapid increase in plant hormones. In total, 57 and 45 differentially accumulated proteins (DAPs) from the embryo and endosperm, respectively, were identified at five germination stages (0, 6, 12, 18, and 24 h). Principal component analysis revealed a significant proteome difference between embryo and endosperm as well as the different germination stages. The largest proteome changes occurred 24 h after seed imbibition. Embryo DAP spots were mainly involved in energy metabolism, amino acid metabolism, stress/defense, and protein metabolism; those from the endosperm were primarily related to storage protein and carbohydrate metabolism. Protein-protein interaction analysis revealed a complicated interaction network between energy-related proteins and other proteins. Metabolic pathway analysis highlighted complex regulatory networks in the embryo and endosperm that regulate wheat seed germination. These results provide new insights into the molecular mechanisms of seed germination. [ABSTRACT FROM AUTHOR]

Escherichia coli BL21(DE3) is an industrial model microbe for the mass-production of bioproducts such as biofuels, biorefineries, and recombinant proteins. However, despite its important role in scientific research and biotechnological applications, a high-quality metabolic network model for metabolic engineering is yet to be developed. Here, we present the comprehensive metabolic network model of E. coli BL21(DE3), named iHK1487, based on the latest genome reannotation and phenome analysis. The metabolic model consists of 1,164 unique metabolites, 2,701 metabolic reactions, and 1,487 genes. The model was validated and improved by comparing the simulation results with phenome data from phenotype microarray tests. Previous transcriptome profile data was incorporated during model reconstruction, and flux prediction was simulated using the model. iHK1487 was simulated to explore the metabolic features of BL21(DE3) such as broad spectrum amino acid utilization and enhanced flux through the upper glycolytic pathway and TCA cycle. iHK1487 will contribute to systematic understanding of cellular physiology and metabolism of E. coli BL21(DE3) and highlight its biotechnological applications. [ABSTRACT FROM AUTHOR]

One of the marvels of biology is the phenotypic plasticity of microorganisms. It allows them to maintain high growth rates across conditions. Studies suggest that cells can express metabolic enzymes at tuned concentrations through adjustment of gene expression. The associated transcription factors are often regulated by intracellular metabolites. Here we study metabolite-mediated regulation of metabolic-gene expression that maximises metabolic fluxes across conditions. We developed an adaptive control theory, qORAC (for ‘Specific Flux (q) Optimization by Robust Adaptive Control’), and illustrate it with several examples of metabolic pathways. The key feature of the theory is that it does not require knowledge of the regulatory network, only of the metabolic part. We derive that maximal metabolic flux can be maintained in the face of varying N environmental parameters only if the number of transcription-factor binding metabolites is at least equal to N. The controlling circuits appear to require simple biochemical kinetics. We conclude that microorganisms likely can achieve maximal rates in metabolic pathways, in the face of environmental changes. [ABSTRACT FROM AUTHOR]