A Method for Finding Metabolic Pathways Using Atomic Group Tracking

Application of coincidence index in the discovery of co-expressed metabolic pathways

Analyzing transcription data requires intensive statistical analysis to obtain useful biological information and knowledge. A significant portion of this data is affected by random noise or even noise intrinsic to the modeling of the experiment. Without robust treatment, the data might not be explored thoroughly, and incorrect conclusions could be drawn. Examining the correlation between gene expression profiles is one way bioinformaticians extract information from transcriptomic experiments. However, the correlation measurements traditionally used have worrisome shortcomings that need to be addressed. This paper compares five already published and experimented-with correlation measurements to the newly developed coincidence index, a similarity measurement that combines Jaccard and interiority indexes and generalizes them to be applied to vectors containing real values. We used microarray and RNA-Seq data from the archaeonHalobacterium salinarumand the bacteriumEscherichia coli, respectively, to evaluate the capacity of each correlation/similarity measurement. The utilized method explores the co-expressed metabolic pathways by measuring the correlations between the expression levels of enzymes that share metabolites, represented in the form of a weighted graph. It then searches for local maxima in this graph using a simulated annealing algorithm. We demonstrate that the coincidence index extracts larger, more comprehensive, and more statistically significant pathways for microarray experiments. In RNA-Seq experiments, the results are more limited, but the coincidence index managed the largest percentage of significant components in the graph.

PyMiner: A method for metabolic pathway design based on the uniform similarity of substrate-product pairs and conditional search

Metabolic pathway design is an essential step in the course of constructing an efficient microbial cell factory to produce high value-added chemicals. Meanwhile, the computational design of biologically meaningful metabolic pathways has been attracting much attention to produce natural and non-natural products. However, there has been a lack of effective methods to perform metabolic network reduction automatically. In addition, comprehensive evaluation indexes for metabolic pathway are still relatively scarce. Here, we define a novel uniform similarity to calculate the main substrate-product pairs of known biochemical reactions, and develop further an efficient metabolic pathway design tool named PyMiner. As a result, the redundant information of general metabolic network (GMN) is eliminated, and the number of substrate-product pairs is shown to decrease by 81.62% on average. Considering that the nodes in the extracted metabolic network (EMN) constructed in this work is large in scale but imbalanced in distribution, we establish a conditional search strategy (CSS) that cuts search time in 90.6% cases. Compared with state-of-the-art methods, PyMiner shows obvious advantages and demonstrates equivalent or better performance on 95% cases of experimentally verified pathways. Consequently, PyMiner is a practical and effective tool for metabolic pathway design.

Revisiting steroidogenesis and its role in immune regulation with the advanced tools and technologies

Historically tools and technologies facilitated scientific discoveries. Steroid hormone research is not an exception. Unfortunately, the dramatic advancement of the field faded this research area and flagged it as a solved topic. However, it should have been the opposite. The area should glitter with its strong foundation and attract next-generation scientists. Over the past century, a myriad of new facts on biochemistry, molecular biology, cell biology, physiology and pathology of the steroid hormones was discovered. Several innovations were made and translated into life-saving treatment strategies such as synthetic steroids, and inhibitors of steroidogenesis and steroid signaling. Steroid molecules exhibit their diverse effects on cell metabolism, salt and water balance, development and function of the reproductive system, pregnancy, and immune-cell function. Despite vigorous research, the molecular basis of the immunomodulatory effect of steroids is still mysterious. The recent excitement on local extra-glandular steroidogenesis in regulating inflammation and immunity is revitalizing the topic with a new perspective. Therefore, here we review the role of steroidogenesis in regulating inflammation and immunity, discuss the unresolved questions, and how this area can bring another golden age of steroid hormone research with the development of new tools and technologies and advancement of the scientific methods.

Finding metabolic pathways in large networks through atom-conserving substrate-product pairs

Motivation

Finding biosynthetic pathways is essential for metabolic engineering of organisms to produce chemicals, biodegradation prediction of pollutants and drugs, and for the elucidation of bioproduction pathways of secondary metabolites. A key step in biosynthetic pathway design is the extraction of novel metabolic pathways from big networks that integrate known biological, as well as novel, predicted biotransformations. However, the efficient analysis and the navigation of big biochemical networks remain a challenge.

Results

Here, we propose the construction of searchable graph representations of metabolic networks. Each reaction is decomposed into pairs of reactants and products, and each pair is assigned a weight, which is calculated from the number of conserved atoms between the reactant and the product molecule. We test our method on a biochemical network that spans 6546 known enzymatic reactions to show how our approach elegantly extracts biologically relevant metabolic pathways from biochemical networks, and how the proposed network structure enables the application of efficient graph search algorithms that improve navigation and pathway identification in big metabolic networks. The weighted reactant–product pairs of an example network and the corresponding graph search algorithm are available online. The proposed method extracts metabolic pathways fast and reliably from big biochemical networks, which is inherently important for all applications involving the engineering of metabolic networks.

Availability and implementation

https://github.com/EPFL-LCSB/nicepath.

Supplementary information

Supplementary data are available at Bioinformatics online.

Finding branched pathways in metabolic network via atom group tracking

Finding non-standard or new metabolic pathways has important applications in metabolic engineering, synthetic biology and the analysis and reconstruction of metabolic networks. Branched metabolic pathways dominate in metabolic networks and depict a more comprehensive picture of metabolism compared to linear pathways. Although progress has been developed to find branched metabolic pathways, few efforts have been made in identifying branched metabolic pathways via atom group tracking. In this paper, we present a pathfinding method called BPFinder for finding branched metabolic pathways by atom group tracking, which aims to guide the synthetic design of metabolic pathways. BPFinder enumerates linear metabolic pathways by tracking the movements of atom groups in metabolic network and merges the linear atom group conserving pathways into branched pathways. Two merging rules based on the structure of conserved atom groups are proposed to accurately merge the branched compounds of linear pathways to identify branched pathways. Furthermore, the integrated information of compound similarity, thermodynamic feasibility and conserved atom groups is also used to rank the pathfinding results for feasible branched pathways. Experimental results show that BPFinder is more capable of recovering known branched metabolic pathways as compared to other existing methods, and is able to return biologically relevant branched pathways and discover alternative branched pathways of biochemical interest. The online server of BPFinder is available at http://114.215.129.245:8080/atomic/. The program, source code and data can be downloaded from https://github.com/hyr0771/BPFinder.

Computational search of branched metabolic pathways is a fundamental problem in metabolic engineering and metabolic network analysis, which provides a systematic way of understanding the metabolism and discovering alternative pathways for synthesis of useful biomolecules. We propose BPFinder, a novel computational approach to identify branched metabolic pathways via atom group tracking. Different from other pathfinding methods using atom tracking, BPFinder tracks the movement of atom groups in metabolic network to find linear atom group conserving pathways, and merge the found linear pathways by the selected branched compounds to generate branched pathways. Based on the structure of conserved atom groups in branched compounds, we design two merging rules for branched compounds: overlapping rule and non-overlapping rule. The user can flexibly adopt these rules to accurately find the branched pathways that contain overlapping/non-overlapping conserved atom groups. BPFinder also enables the user to combine the information of compound similarity, Gibbs free energy of reactions, and conserved atom groups to sort resulting pathways. Compared with other existing methods, BPFinder can more accurately recover the known branched pathways. The alternative branched pathways returned by BPFinder reveal that the user can flexibly utilize our proposed merging rules to discover biochemically meaningful pathways of interest.

NetFlow: A tool for isolating carbon flows in genome-scale metabolic networks

Genome-scale stoichiometric models (GSMs) have been widely utilized to predict and understand cellular metabolism. GSMs and the flux predictions resulting from them have proven indispensable to fields ranging from metabolic engineering to human disease. Nonetheless, it is challenging to parse these flux predictions due to the inherent size and complexity of the GSMs. Several previous approaches have reduced this complexity by identifying key pathways contained within the genome-scale flux predictions. However, a reduction method that overlays carbon atom transitions on stoichiometry and flux predictions is lacking. To fill this gap, we developed NetFlow, an algorithm that leverages genome-scale carbon mapping to extract and quantitatively distinguish biologically relevant metabolic pathways from a given genome-scale flux prediction. NetFlow extends prior approaches by utilizing both full carbon mapping and context-specific flux predictions. Thus, NetFlow is uniquely able to quantitatively distinguish between biologically relevant pathways of carbon flow within the given flux map. NetFlow simulates ¹³C isotope labeling experiments to calculate the extent of carbon exchange, or carbon yield, between every metabolite in the given GSM. Based on the carbon yield, the carbon flow to or from any metabolite or between any pair of metabolites of interest can be isolated and readily visualized. The resulting pathways are much easier to interpret, which enables an in-depth mechanistic understanding of the metabolic phenotype of interest. Here, we first demonstrate NetFlow with a simple network. We then depict the utility of NetFlow on a model of central carbon metabolism in E. coli. Specifically, we isolated the production pathway for succinate synthesis in this model and the metabolic mechanism driving the predicted increase in succinate yield in a double knockout of E. coli. Finally, we describe the application of NetFlow to a GSM of lycopene-producing E. coli, which enabled the rapid identification of the mechanisms behind the measured increases in lycopene production following single, double, and triple knockouts.

PhDSeeker: Pheromone-Directed Seeker for metabolic pathways

Manually finding relationship networks among compounds can be a hard and time-consuming task. However, this process is fundamental when looking for a metabolic pathway that explains how multiple compounds are related, to identify relevant pathways in organisms, filling gaps on metabolic networks, or when new mechanisms for the synthesis of important compounds are sought. Here, we present PhDSeeker, a new tool for the automatic search of metabolic pathways. This tool is able to relate simultaneously several compounds. Furthermore, its flexibility allows it to be easily configured for addressing a wide range of situations. Solutions found are provided not only in plain text but also as interactive representations that can be analyzed in a web browser. Source code is available at https://github.com/sinc-lab/phdseeker. A web service is also available at https://sinc.unl.edu.ar/web-demo/phds/. Several fully documented study cases, including their settings and solutions files, are also provided as Supplementary Material.

A system architecture for parallel analysis of flux-balanced metabolic pathways

Elementary flux mode (EFM) analysis is a well-studied method in constraint-based modeling of metabolic networks. In EFM analysis, a network is decomposed into minimal functional pathways based on the assumption of balanced metabolic fluxes. In this paper, a system architecture is proposed that approximately models the functionality of metabolic networks. The AND/OR graph model is used to represent the metabolic network and each processing element in the system emulates the functionality of a metabolite. The system is implemented on a graphics processing unit (GPU) as the hardware platform using CUDA environment. The proposed architecture takes advantage of the inherent parallelism in the network structure in terms of both pathway and metabolite traversal. The function of each element is defined such that it can find flux-balanced pathways. Pathways in both small and large metabolic networks are applied to the proposed architecture and the results are discussed.

MAPPS: A Web-Based Tool for Metabolic Pathway Prediction and Network Analysis in the Postgenomic Era

Comparative and evolutionary analyses of metabolic networks have a wide range of applications, ranging from research into metabolic evolution through to practical applications in drug development, synthetic biology, and biodegradation. We present MAPPS: Metabolic network Analysis and Pathway Prediction Server (https://mapps.lums.edu.pk), a web-based tool to study functions and evolution of metabolic networks using traditional and 'omics data sets. MAPPS provides diverse functionalities including an interactive interface, graphical visualization of results, pathway prediction and network comparison, identification of potential drug targets, in silico metabolic engineering, host-microbe interactions, and ancestral network building. Importantly, MAPPS also allows users to upload custom data, thus enabling metabolic analyses on draft and custom genomes, and has an 'omics pipeline to filter pathway results, making it relevant in today's postgenomic era.

Nicotinamide attenuates the decrease in dendritic spine density in hippocampal primary neurons from 5xFAD mice, an Alzheimer's disease animal model

Alzheimer’s disease (AD) is the most common neurodegenerative disease characterized by memory loss and the presence of amyloid plaques and neurofibrillary tangles in the patients’ brains. In this study, we investigated the alterations in metabolite profiles of the hippocampal tissues from 6, 8, and 12 month-old wild-type (WT) and 5xfamiliar AD (5xFAD) mice, an AD mouse model harboring 5 early-onset familiar AD mutations, which shows memory loss from approximately 5 months of age, by exploiting the untargeted metabolomics profiling. We found that nicotinamide and adenosine monophosphate levels have been significantly decreased while lysophosphatidylcholine (LysoPC) (16:0), LysoPC (18:0), and lysophosphatidylethanolamine (LysoPE) (16:0) levels have been significantly increased in the hippocampi from 5xFAD mice at 8 months or 12 months of age, compared to those from age-matched wild-type mice. In the present study, we focused on the role of nicotinamide and examined if replenishment of nicotinamide exerts attenuating effects on the reduction in dendritic spine density in hippocampal primary neurons from 5xFAD mice. Treatment with nicotinamide attenuated the deficits in spine density in the hippocampal primary neurons derived from 5xFAD mice, indicating a potential role of nicotinamide in the pathogenesis of AD. Taken together, these findings suggest that the decreased hippocampal nicotinamide level could be linked with AD pathogenesis and be a useful therapeutic target for AD.

Improving the organization and interactivity of metabolic pathfinding with precomputed pathways

Background

The rapid growth of available knowledge on metabolic processes across thousands of species continues to expand the possibilities of producing chemicals by combining pathways found in different species. Several computational search algorithms have been developed for automating the identification of possible heterologous pathways; however, these searches may return thousands of pathway results. Although the large number of results are in part due to the large number of possible compounds and reactions, a subset of core reaction modules is repeatedly observed in pathway results across multiple searches, suggesting that some subpaths between common compounds were more consistently explored than others.To reduce the resources spent on searching the same metabolic space, a new meta-algorithm for metabolic pathfinding, Hub Pathway search with Atom Tracking (HPAT), was developed to take advantage of a precomputed network of subpath modules. To investigate the efficacy of this method, we created a table describing a network of common hub metabolites and how they are biochemically connected and only offloaded searches to and from this hub network onto an interactive webserver capable of visualizing the resulting pathways.

Results

A test set of nineteen known pathways taken from literature and metabolic databases were used to evaluate if HPAT was capable of identifying known pathways. HPAT found the exact pathway for eleven of the nineteen test cases using a diverse set of precomputed subpaths, whereas a comparable pathfinding search algorithm that does not use precomputed subpaths found only seven of the nineteen test cases. The capability of HPAT to find novel pathways was demonstrated by its ability to identify novel 3-hydroxypropanoate (3-HP) synthesis pathways. As for pathway visualization, the new interactive pathway filters enable a reduction of the number of displayed pathways from hundreds down to less than ten pathways in several test cases, illustrating their utility in reducing the amount of presented information while retaining pathways of interest.

Conclusions

This work presents the first step in incorporating a precomputed subpath network into metabolic pathfinding and demonstrates how this leads to a concise, interactive visualization of pathway results. The modular nature of metabolic pathways is exploited to facilitate efficient discovery of alternate pathways.

Detecting list-colored graph motifs in biological networks using branch-and-bound strategy

In this work, we study the list-colored graph motif problem, which was introduced to detect functional motifs in biological networks. Given a multi-set M of colors as the query motif and a list-colored graph G where each vertex in G is associated with a set of colors, the aim of this problem is to find a sub-graph of G whose vertex set is colored exactly as motif M. To solve this problem, we present a heuristic method to efficiently and accurately detect list-colored graph motifs in biological networks using branch-and-bound strategy. We transform the detection of list-colored graph motif to the search of connected induced sub-graphs in list-colored graph, where the vertices in the sub-graph are assigned to distinctive colors of query motif. This transformation enables our method to accurately discover the occurrences of query motif without enumerating and verifying all sub-graphs. Furthermore, a new initial vertex selection strategy based on the colors of vertices is proposed to accurately determine the search scope of motifs. Experiments conducted on metabolic networks and protein-interaction networks demonstrate that our method can achieve better performance in accuracy and efficiency in comparison to other existing methods.

Metabolic pathways synthesis based on ant colony optimization

One of the current challenges in bioinformatics is to discover new ways to transform a set of compounds into specific products. The usual approach is finding the reactions to synthesize a particular product, from a given substrate, by means of classical searching algorithms. However, they have three main limitations: difficulty in handling large amounts of reactions and compounds; absence of a step that verifies the availability of substrates; and inability to find branched pathways. We present here a novel bio-inspired algorithm for synthesizing linear and branched metabolic pathways. It allows relating several compounds simultaneously, ensuring the availability of substrates for every reaction in the solution. Comparisons with classical searching algorithms and other recent metaheuristic approaches show clear advantages of this proposal, fully recovering well-known pathways. Furthermore, solutions found can be analyzed in a simple way through graphical representations on the web.

A Prospective Study on the Fermentation Landscape of Gaseous Substrates to Biorenewables Using Methanosarcina acetivorans Metabolic Model

The abundance of methane in shale gas and of other gases such as carbon monoxide, hydrogen, and carbon dioxide as chemical process byproducts has motivated the use of gas fermentation for bioproduction. Recent advances in metabolic engineering and synthetic biology allow for engineering of microbes metabolizing a variety of chemicals including gaseous feeds into a number of biorenewables and transportation liquid fuels. New computational tools enable the systematic exploration of all feasible conversion alternatives. Here we computationally assessed all thermodynamically feasible ways of co-utilizing CH₄, CO, and CO₂ using ferric as terminal electron acceptor for the production of all key precursor metabolites. We identified the thermodynamically feasible co-utilization ratio ranges of CH₄, CO, and CO₂ toward production of the target metabolite(s) as a function of ferric uptake. A revised version of the iMAC868 genome-scale metabolic model of Methanosarcina acetivorans was chosen to assess co-utilization of CH₄, CO, and CO₂ and their conversion into selected target products using the optStoic pathway design tool. This revised version contains the latest information on electron flow mechanisms by the methanogen while supplied with methane as the sole carbon source. The interplay between different gas co-utilization ratios and the energetics of reverse methanogenesis were also analyzed using the same metabolic model.

A pathway for every product? Tools to discover and design plant metabolism

The vast diversity of plant natural products is a powerful indication of the biosynthetic capacity of plant metabolism. Synthetic biology seeks to capitalize on this ability by understanding and reconfiguring the biosynthetic pathways that generate this diversity to produce novel products with improved efficiency. Here we review the algorithms and databases that presently support the design and manipulation of metabolic pathways in plants, starting from metabolic models of native biosynthetic pathways, progressing to novel combinations of known reactions, and finally proposing new reactions that may be carried out by existing enzymes. We show how these tools are useful for proposing new pathways as well as identifying side reactions that may affect engineering goals.

Reconstructing Phylogeny by Aligning Multiple Metabolic Pathways Using Functional Module Mapping

Comparison of metabolic pathways provides a systematic way for understanding the evolutionary and phylogenetic relationships in systems biology. Although a number of phylogenetic methods have been developed, few efforts have been made to provide a unified phylogenetic framework that sufficiently reflects the metabolic features of organisms. In this paper, we propose a phylogenetic framework that characterizes the metabolic features of organisms by aligning multiple metabolic pathways using functional module mapping. Our method transforms the alignment of multiple metabolic pathways into constructing the union graph of pathways, builds mappings between functional modules of pathways in the union graph, and infers phylogenetic relationships among organisms based on module mappings. Experimental results show that the use of functional module mapping enables us to correctly categorize organisms into main categories with specific metabolic characteristics. Traditional genome-based phylogenetic methods can reconstruct phylogenetic relationships, whereas our method can offer in-depth metabolic analysis for phylogenetic reconstruction, which can add insights into traditional phyletic reconstruction. The results also demonstrate that our phylogenetic trees are closer to the classic classifications in comparison to existing classification methods using metabolic pathway data.

A review of computational tools for design and reconstruction of metabolic pathways

Metabolic pathways reflect an organism's chemical repertoire and hence their elucidation and design have been a primary goal in metabolic engineering. Various computational methods have been developed to design novel metabolic pathways while taking into account several prerequisites such as pathway stoichiometry, thermodynamics, host compatibility, and enzyme availability. The choice of the method is often determined by the nature of the metabolites of interest and preferred host organism, along with computational complexity and availability of software tools. In this paper, we review different computational approaches used to design metabolic pathways based on the reaction network representation of the database (i.e., graph or stoichiometric matrix) and the search algorithm (i.e., graph search, flux balance analysis, or retrosynthetic search). We also put forth a systematic workflow that can be implemented in projects requiring pathway design and highlight current limitations and obstacles in computational pathway design.

A review of parameters and heuristics for guiding metabolic pathfinding

Recent developments in metabolic engineering have led to the successful biosynthesis of valuable products, such as the precursor of the antimalarial compound, artemisinin, and opioid precursor, thebaine. Synthesizing these traditionally plant-derived compounds in genetically modified yeast cells introduces the possibility of significantly reducing the total time and resources required for their production, and in turn, allows these valuable compounds to become cheaper and more readily available. Most biosynthesis pathways used in metabolic engineering applications have been discovered manually, requiring a tedious search of existing literature and metabolic databases. However, the recent rapid development of available metabolic information has enabled the development of automated approaches for identifying novel pathways. Computer-assisted pathfinding has the potential to save biochemists time in the initial discovery steps of metabolic engineering. In this paper, we review the parameters and heuristics used to guide the search in recent pathfinding algorithms. These parameters and heuristics capture information on the metabolic network structure, compound structures, reaction features, and organism-specificity of pathways. No one metabolic pathfinding algorithm or search parameter stands out as the best to use broadly for solving the pathfinding problem, as each method and parameter has its own strengths and shortcomings. As assisted pathfinding approaches continue to become more sophisticated, the development of better methods for visualizing pathway results and integrating these results into existing metabolic engineering practices is also important for encouraging wider use of these pathfinding methods.

Network design and analysis for multi-enzyme biocatalysis

BACKGROUND

As more and more biological reaction data become available, the full exploration of the enzymatic potential for the synthesis of valuable products opens up exciting new opportunities but is becoming increasingly complex. The manual design of multi-step biosynthesis routes involving enzymes from different organisms is very challenging. To harness the full enzymatic potential, we developed a computational tool for the directed design of biosynthetic production pathways for multi-step catalysis with in vitro enzyme cascades, cell hydrolysates and permeabilized cells.

RESULTS

We present a method which encompasses the reconstruction of a genome-scale pan-organism metabolic network, path-finding and the ranking of the resulting pathway candidates for proposing suitable synthesis pathways. The network is based on reaction and reaction pair data from the Kyoto Encyclopedia of Genes and Genomes (KEGG) and the thermodynamics calculator eQuilibrator. The pan-organism network is especially useful for finding the most suitable pathway to a target metabolite from a thermodynamic or economic standpoint. However, our method can be used with any network reconstruction, e.g. for a specific organism. We implemented a path-finding algorithm based on a mixed-integer linear program (MILP) which takes into account both topology and stoichiometry of the underlying network. Unlike other methods we do not specify a single starting metabolite, but our algorithm searches for pathways starting from arbitrary start metabolites to a target product of interest. Using a set of biochemical ranking criteria including pathway length, thermodynamics and other biological characteristics such as number of heterologous enzymes or cofactor requirement, it is possible to obtain well-designed meaningful pathway alternatives. In addition, a thermodynamic profile, the overall reactant balance and potential side reactions as well as an SBML file for visualization are generated for each pathway alternative.

CONCLUSION

We present an in silico tool for the design of multi-enzyme biosynthetic production pathways starting from a pan-organism network. The method is highly customizable and each module can be adapted to the focus of the project at hand. This method is directly applicable for (i) in vitro enzyme cascades, (ii) cell hydrolysates and (iii) permeabilized cells.