MultiMetEval: comparative and multi-objective analysis of genome-scale metabolic models

Comparative metabolic modelling is emerging as a novel field, supported by the development of reliable and standardized approaches for constructing genome-scale metabolic models in high throughput. New software solutions are needed to allow efficient comparative analysis of multiple models in the context of multiple cellular objectives. Here, we present the user-friendly software framework Multi-Metabolic Evaluator (MultiMetEval), built upon SurreyFBA, which allows the user to compose collections of metabolic models that together can be subjected to flux balance analysis. Additionally, MultiMetEval implements functionalities for multi-objective analysis by calculating the Pareto front between two cellular objectives. Using a previously generated dataset of 38 actinobacterial genome-scale metabolic models, we show how these approaches can lead to exciting novel insights. Firstly, after incorporating several pathways for the biosynthesis of natural products into each of these models, comparative flux balance analysis predicted that species like Streptomyces that harbour the highest diversity of secondary metabolite biosynthetic gene clusters in their genomes do not necessarily have the metabolic network topology most suitable for compound overproduction. Secondly, multi-objective analysis of biomass production and natural product biosynthesis in these actinobacteria shows that the well-studied occurrence of discrete metabolic switches during the change of cellular objectives is inherent to their metabolic network architecture. Comparative and multi-objective modelling can lead to insights that could not be obtained by normal flux balance analyses. MultiMetEval provides a powerful platform that makes these analyses straightforward for biologists. Sources and binaries of MultiMetEval are freely available from https://github.com/PiotrZakrzewski/MetEv​al/downloads

Characterization of growth and metabolism of the haloalkaliphile Natronomonas pharaonis

Natronomonas pharaonis is an archaeon adapted to two extreme conditions: high salt concentration and alkaline pH. It has become one of the model organisms for the study of extremophilic life. Here, we present a genome-scale, manually curated metabolic reconstruction for the microorganism. The reconstruction itself represents a knowledge base of the haloalkaliphile's metabolism and, as such, would greatly assist further investigations on archaeal pathways. In addition, we experimentally determined several parameters relevant to growth, including a characterization of the biomass composition and a quantification of carbon and oxygen consumption. Using the metabolic reconstruction and the experimental data, we formulated a constraints-based model which we used to analyze the behavior of the archaeon when grown on a single carbon source. Results of the analysis include the finding that Natronomonas pharaonis, when grown aerobically on acetate, uses a carbon to oxygen consumption ratio that is theoretically near-optimal with respect to growth and energy production. This supports the hypothesis that, under simple conditions, the microorganism optimizes its metabolism with respect to the two objectives. We also found that the archaeon has a very low carbon efficiency of only about 35%. This inefficiency is probably due to a very low P/O ratio as well as to the other difficulties posed by its extreme environment

Potentiating antibacterial activity by predictably enhancing endogenous microbial ROS production

The ever-increasing incidence of antibiotic-resistant infections combined with a weak pipeline of new antibiotics has created a global public health crisis1. Accordingly, novel strategies for enhancing our antibiotic arsenal are needed. As antibiotics kill bacteria in part by inducing reactive oxygen species (ROS)2–4, we reasoned that targeting microbial ROS production might potentiate antibiotic activity. Here we show that ROS production can be predictably enhanced in Escherichia coli, increasing the bacteria’s susceptibility to oxidative attack. We developed an ensemble, genome-scale metabolic modeling approach capable of predicting ROS production in E. coli. The metabolic network was systematically perturbed and its flux distribution analyzed to identify targets predicted to increase ROS production. In silico–predicted targets were experimentally validated and shown to confer increased susceptibility to oxidants. Validated targets also increased susceptibility to killing by antibiotics. This work establishes a systems-based method to tune ROS production in bacteria and demonstrates that increased microbial ROS production can potentiate killing by oxidants and antibiotics

Centralized Modularity of N-Linked Glycosylation Pathways in Mammalian Cells

Glycosylation is a highly complex process to produce a diverse repertoire of cellular glycans that are attached to proteins and lipids. Glycans are involved in fundamental biological processes, including protein folding and clearance, cell proliferation and apoptosis, development, immune responses, and pathogenesis. One of the major types of glycans, N-linked glycans, is formed by sequential attachments of monosaccharides to proteins by a limited number of enzymes. Many of these enzymes can accept multiple N-linked glycans as substrates, thereby generating a large number of glycan intermediates and their intermingled pathways. Motivated by the quantitative methods developed in complex network research, we investigated the large-scale organization of such N-linked glycosylation pathways in mammalian cells. The N-linked glycosylation pathways are extremely modular, and are composed of cohesive topological modules that directly branch from a common upstream pathway of glycan synthesis. This unique structural property allows the glycan production between modules to be controlled by the upstream region. Although the enzymes act on multiple glycan substrates, indicating cross-talk between modules, the impact of the cross-talk on the module-specific enhancement of glycan synthesis may be confined within a moderate range by transcription-level control. The findings of the present study provide experimentally-testable predictions for glycosylation processes, and may be applicable to therapeutic glycoprotein engineering

Optimal flux spaces of genome-scale stoichiometric models are determined by a few subnetworks

The metabolism of organisms can be studied with comprehensive stoichiometric models of their metabolic networks. Flux balance analysis (FBA) calculates optimal metabolic performance of stoichiometric models. However, detailed biological interpretation of FBA is limited because, in general, a huge number of flux patterns give rise to the same optimal performance. The complete description of the resulting optimal solution spaces was thus far a computationally intractable problem. Here we present CoPE-FBA: Comprehensive Polyhedra Enumeration Flux Balance Analysis, a computational method that solves this problem. CoPE-FBA indicates that the thousands to millions of optimal flux patterns result from a combinatorial explosion of flux patterns in just a few metabolic sub-networks. The entire optimal solution space can now be compactly described in terms of the topology of these sub-networks. CoPE-FBA simplifies the biological interpretation of stoichiometric models of metabolism, and provides a profound understanding of metabolic flexibility in optimal states

Development of a framework for metabolic pathway analysis-driven strain optimization methods

Genome-scale metabolic models (GSMMs) have become important assets for rational design of compound overproduction using microbial cell factories. Most computational strain optimization methods (CSOM) using GSMMs, while useful in metabolic engineering, rely on the definition of questionable cell objectives, leading to some bias. Metabolic pathway analysis approaches do not require an objective function. Though their use brings immediate advantages, it has mostly been restricted to small scale models due to computational demands. Additionally, their complex parameterization and lack of intuitive tools pose an important challenge towards making these widely available to the community. Recently, MCSEnumerator has extended the scale of these methods, namely regarding enumeration of minimal cut sets, now able to handle GSMMs. This work proposes a tool implementing this method as a Java library and a plugin within the OptFlux metabolic engineering platform providing a friendly user interface. A standard enumeration problem and pipeline applicable to GSMMs is proposed, making use by the community simpler. To highlight the potential of these approaches, we devised a case study for overproduction of succinate, providing a phenotype analysis of a selected strategy and comparing robustness with a selected solution from a bi-level CSOM.The authors thank the project “DeYeastLibrary—Designer yeast strain library optimized for metabolic engineering applications”, Ref. ERA-IB-2/0003/2013, funded by national funds through “Fundação para a Ciência e Tecnologia / Ministério da Ciência, Tecnologia e Ensino Superior”.info:eu-repo/semantics/publishedVersio

Future therapeutic targets in rheumatoid arthritis?

Rheumatoid arthritis (RA) is a chronic inflammatory disease characterized by persistent joint inflammation. Without adequate treatment, patients with RA will develop joint deformity and progressive functional impairment. With the implementation of treat-to-target strategies and availability of biologic therapies, the outcomes for patients with RA have significantly improved. However, the unmet need in the treatment of RA remains high as some patients do not respond sufficiently to the currently available agents, remission is not always achieved and refractory disease is not uncommon. With better understanding of the pathophysiology of RA, new therapeutic approaches are emerging. Apart from more selective Janus kinase inhibition, there is a great interest in the granulocyte macrophage-colony stimulating factor pathway, Bruton's tyrosine kinase pathway, phosphoinositide-3-kinase pathway, neural stimulation and dendritic cell-based therapeutics. In this review, we will discuss the therapeutic potential of these novel approaches

An integrated network visualization framework towards metabolic engineering applications

Background Over the last years, several methods for the phenotype simulation of microorganisms, under specified genetic and environmental conditions have been proposed, in the context of Metabolic Engineering (ME). These methods provided insight on the functioning of microbial metabolism and played a key role in the design of genetic modifications that can lead to strains of industrial interest. On the other hand, in the context of Systems Biology research, biological network visualization has reinforced its role as a core tool in understanding biological processes. However, it has been scarcely used to foster ME related methods, in spite of the acknowledged potential. Results In this work, an open-source software that aims to fill the gap between ME and metabolic network visualization is proposed, in the form of a plugin to the OptFlux ME platform. The framework is based on an abstract layer, where the network is represented as a bipartite graph containing minimal information about the underlying entities and their desired relative placement. The framework provides input/output support for networks specified in standard formats, such as XGMML, SBGN or SBML, providing a connection to genome-scale metabolic models. An user-interface makes it possible to edit, manipulate and query nodes in the network, providing tools to visualize diverse effects, including visual filters and aspect changing (e.g. colors, shapes and sizes). These tools are particularly interesting for ME, since they allow overlaying phenotype simulation results or elementary flux modes over the networks. Conclusions The framework and its source code are freely available, together with documentation and other resources, being illustrated with well documented case studies.This work is partially funded by ERDF - European Regional Development Fund through the COMPETE Programme (operational programme for competitiveness) and by National Funds through the FCT (Portuguese Foundation for Science and Technology) within project ref. COMPETE FCOMP-01-0124-FEDER-015079 and the FCT Strategic Project PEst-OE/EQB/LA0023/2013. The work of PV is funded by PhD grant ref. SFRH/BDE/51442/2011

Meneco, a Topology-Based Gap-Filling Tool Applicable to Degraded Genome-Wide Metabolic Networks

International audienceIncreasing amounts of sequence data are becoming available for a wide range of non-model organisms. Investigating and modelling the metabolic behaviour of those organisms is highly relevant to understand their biology and ecology. As sequences are often incomplete and poorly annotated, draft networks of their metabolism largely suffer from incompleteness. Appropriate gap-filling methods to identify and add missing reactions are therefore required to address this issue. However, current tools rely on phenotypic or taxonomic information, or are very sensitive to the stoichiometric balance of metabolic reactions, especially concerning the co-factors. This type of information is often not available or at least prone to errors for newly-explored organisms. Here we introduce Meneco, a tool dedicated to the topological gap-filling of genome-scale draft metabolic networks. Meneco reformulates gap-filling as a qualitative combinatorial optimization problem, omitting constraints raised by the stoichiometry of a metabolic network considered in other methods, and solves this problem using Answer Set Programming. Run on several artificial test sets gathering 10,800 degraded Escherichia coli networks Meneco was able to efficiently identify essential reactions missing in networks at high degradation rates, outperforming the stoichiometry-based tools in scalability. To demonstrate the utility of Meneco we applied it to two case studies. Its application to recent metabolic networks reconstructed for the brown algal model Ectocarpus siliculosus and an associated bacterium Candidatus Phaeomarinobacter ectocarpi revealed several candidate metabolic pathways for algal-bacterial interactions. Then Meneco was used to reconstruct, from transcriptomic and metabolomic data, the first metabolic network for the microalga Euglena mutabilis. These two case studies show that Meneco is a versatile tool to complete draft genome-scale metabolic networks produced from heterogeneous data, and to suggest relevant reactions that explain the metabolic capacity of a biological system

Multi-omic data integration elucidates Synechococcus adaptation mechanisms to fluctuations in light intensity and salinity

Synechococcus sp. PCC 7002 is a fast-growing cyanobacterium which flourishes in freshwater and marine environments, owing to its ability to tolerate high light intensity and a wide range of salinities. Harnessing the properties of cyanobacteria and understanding their metabolic efficiency has become an imperative goal in recent years owing to their potential to serve as biocatalysts for the production of renewable biofuels. To improve characterisation of metabolic networks, genome-scale models of metabolism can be integrated with multi-omic data to provide a more accurate representation of metabolic capability and refine phenotypic predictions. In this work, a heuristic pipeline is constructed for analysing a genome-scale metabolic model of Synechococcus sp. PCC 7002, which utilises flux balance analysis across multiple layers to observe flux response between conditions across four key pathways. Across various conditions, the detection of significant patterns and mechanisms to cope with fluctuations in light intensity and salinity provides insights into the maintenance of metabolic efficiency