Comparing methods for metabolic network analysis and an application to metabolic engineering

Metabolic Connectome and Its Role in the Prediction, Diagnosis, and Treatment of Complex Diseases

The interconnectivity of advanced biological systems is essential for their proper functioning. In modern connectomics, biological entities such as proteins, genes, RNA, DNA, and metabolites are often represented as nodes, while the physical, biochemical, or functional interactions between them are represented as edges. Among these entities, metabolites are particularly significant as they exhibit a closer relationship to an organism's phenotype compared to genes or proteins. Moreover, the metabolome has the ability to amplify small proteomic and transcriptomic changes, even those from minor genomic changes. Metabolic networks, which consist of complex systems comprising hundreds of metabolites and their interactions, play a critical role in biological research by mediating energy conversion and chemical reactions within cells. This review provides an introduction to common metabolic network models and their construction methods. It also explores the diverse applications of metabolic networks in elucidating disease mechanisms, predicting and diagnosing diseases, and facilitating drug development. Additionally, it discusses potential future directions for research in metabolic networks. Ultimately, this review serves as a valuable reference for researchers interested in metabolic network modeling, analysis, and their applications.

Design and thermodynamic analysis of a pathway enabling anaerobic production of poly-3-hydroxybutyrate in Escherichia coli

Utilizing anaerobic metabolisms for the production of biotechnologically relevant products presents potential advantages, such as increased yields and reduced energy dissipation. However, lower energy dissipation may indicate that certain reactions are operating closer to their thermodynamic equilibrium. While stoichiometric analyses and genetic modifications are frequently employed in metabolic engineering, the use of thermodynamic tools to evaluate the feasibility of planned interventions is less documented. In this study, we propose a novel metabolic engineering strategy to achieve an efficient anaerobic production of poly-(R)-3-hydroxybutyrate (PHB) in the model organism Escherichia coli. Our approach involves re-routing of two-thirds of the glycolytic flux through non-oxidative glycolysis and coupling PHB synthesis with NADH re-oxidation. We complemented our stoichiometric analysis with various thermodynamic approaches to assess the feasibility and the bottlenecks in the proposed engineered pathway. According to our calculations, the main thermodynamic bottleneck are the reactions catalyzed by the acetoacetyl-CoA β-ketothiolase (EC 2.3.1.9) and the acetoacetyl-CoA reductase (EC 1.1.1.36). Furthermore, we calculated thermodynamically consistent sets of kinetic parameters to determine the enzyme amounts required for sustaining the conversion fluxes. In the case of the engineered conversion route, the protein pool necessary to sustain the desired fluxes could account for 20% of the whole cell dry weight.

Metabolomics and modelling approaches for systems metabolic engineering

Metabolic engineering involves the manipulation of microbes to produce desirable compounds through genetic engineering or synthetic biology approaches. Metabolomics involves the quantitation of intracellular and extracellular metabolites, where mass spectrometry and nuclear magnetic resonance based analytical instrumentation are often used. Here, the experimental designs, sample preparations, metabolite quenching and extraction are essential to the quantitative metabolomics workflow. The resultant metabolomics data can then be used with computational modelling approaches, such as kinetic and constraint-based modelling, to better understand underlying mechanisms and bottlenecks in the synthesis of desired compounds, thereby accelerating research through systems metabolic engineering. Constraint-based models, such as genome scale models, have been used successfully to enhance the yield of desired compounds from engineered microbes, however, unlike kinetic or dynamic models, constraint-based models do not incorporate regulatory effects. Nevertheless, the lack of time-series metabolomic data generation has hindered the usefulness of dynamic models till today. In this review, we show that improvements in automation, dynamic real-time analysis and high throughput workflows can drive the generation of more quality data for dynamic models through time-series metabolomics data generation. Spatial metabolomics also has the potential to be used as a complementary approach to conventional metabolomics, as it provides information on the localization of metabolites. However, more effort must be undertaken to identify metabolites from spatial metabolomics data derived through imaging mass spectrometry, where machine learning approaches could prove useful. On the other hand, single-cell metabolomics has also seen rapid growth, where understanding cell-cell heterogeneity can provide more insights into efficient metabolic engineering of microbes. Moving forward, with potential improvements in automation, dynamic real-time analysis, high throughput workflows, and spatial metabolomics, more data can be produced and studied using machine learning algorithms, in conjunction with dynamic models, to generate qualitative and quantitative predictions to advance metabolic engineering efforts.

In Silico Design Strategies for the Production of Target Chemical Compounds Using Iterative Single-Level Linear Programming Problems

The optimization of metabolic reaction modifications for the production of target compounds is a complex computational problem whose execution time increases exponentially with the number of metabolic reactions. Therefore, practical technologies are needed to identify reaction deletion combinations to minimize computing times and promote the production of target compounds by modifying intracellular metabolism. In this paper, a practical metabolic design technology named AERITH is proposed for high-throughput target compound production. This method can optimize the production of compounds of interest while maximizing cell growth. With this approach, an appropriate combination of metabolic reaction deletions can be identified by solving a simple linear programming problem. Using a standard CPU, the computation time could be as low as 1 min per compound, and the system can even handle large metabolic models. AERITH was implemented in MATLAB and is freely available for non-profit use.

Automated Exploration of Prebiotic Chemical Reaction Space: Progress and Perspectives

Prebiotic chemistry often involves the study of complex systems of chemical reactions that form large networks with a large number of diverse species. Such complex systems may have given rise to emergent phenomena that ultimately led to the origin of life on Earth. The environmental conditions and processes involved in this emergence may not be fully recapitulable, making it difficult for experimentalists to study prebiotic systems in laboratory simulations. Computational chemistry offers efficient ways to study such chemical systems and identify the ones most likely to display complex properties associated with life. Here, we review tools and techniques for modelling prebiotic chemical reaction networks and outline possible ways to identify self-replicating features that are central to many origin-of-life models.

Physicochemical and metabolic constraints for thermodynamics-based stoichiometric modelling under mesophilic growth conditions

Metabolic engineering in the post-genomic era is characterised by the development of new methods for metabolomics and fluxomics, supported by the integration of genetic engineering tools and mathematical modelling. Particularly, constraint-based stoichiometric models have been widely studied: (i) flux balance analysis (FBA) (in silico), and (ii) metabolic flux analysis (MFA) (in vivo). Recent studies have enabled the incorporation of thermodynamics and metabolomics data to improve the predictive capabilities of these approaches. However, an in-depth comparison and evaluation of these methods is lacking. This study presents a thorough analysis of two different in silico methods tested against experimental data (metabolomics and ¹³C-MFA) for the mesophile Escherichia coli. In particular, a modified version of the recently published matTFA toolbox was created, providing a broader range of physicochemical parameters. Validating against experimental data allowed the determination of the best physicochemical parameters to perform the TFA (Thermodynamics-based Flux Analysis). An analysis of flux pattern changes in the central carbon metabolism between ¹³C-MFA and TFA highlighted the limited capabilities of both approaches for elucidating the anaplerotic fluxes. In addition, a method based on centrality measures was suggested to identify important metabolites that (if quantified) would allow to further constrain the TFA. Finally, this study emphasised the need for standardisation in the fluxomics community: novel approaches are frequently released but a thorough comparison with currently accepted methods is not always performed.

Biotechnology has benefitted from the development of high throughput methods characterising living systems at different levels (e.g. concerning genes or proteins), allowing the industrial production of chemical commodities. Recently, focus has been placed on determining reaction rates (or metabolic fluxes) in the metabolic network of certain microorganisms, in order to identify bottlenecks hindering their exploitation. Two main approaches are commonly used, termed metabolic flux analysis (MFA) and flux balance analysis (FBA), based on measuring and estimating fluxes, respectively. While the influence of thermodynamics in living systems was accepted several decades ago, its application to study biochemical networks has only recently been enabled. In this sense, a multitude of different approaches constraining well-established modelling methods with thermodynamics has been suggested. However, physicochemical parameters are generally not properly adjusted to the experimental conditions, which might affect their predictive capabilities. In this study, we have explored the reliability of currently available tools by investigating the impact of varying said parameters in the simulation of metabolic fluxes and metabolite concentration values. Additionally, our in-depth analysis allowed us to highlight limitations and potential solutions that should be considered in future studies.

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.

Comparison of Optimization-Modelling Methods for Metabolites Production in Escherichia coli

The metabolic network is the reconstruction of the metabolic pathway of an organism that is used to represent the interaction between enzymes and metabolites in genome level. Meanwhile, metabolic engineering is a process that modifies the metabolic network of a cell to increase the production of metabolites. However, the metabolic networks are too complex that cause problem in identifying near-optimal knockout genes/reactions for maximizing the metabolite’s production. Therefore, through constraint-based modelling, various metaheuristic algorithms have been improvised to optimize the desired phenotypes. In this paper, PSOMOMA was compared with CSMOMA and ABCMOMA for maximizing the production of succinic acid in E. coli. Furthermore, the results obtained from PSOMOMA were validated with results from the wet lab experiment.

Construction, comparison and evolution of networks in life sciences and other disciplines

Network approaches have become pervasive in many research fields. They allow for a more comprehensive understanding of complex relationships between entities as well as their group-level properties and dynamics. Many networks change over time, be it within seconds or millions of years, depending on the nature of the network. Our focus will be on comparative network analyses in life sciences, where deciphering temporal network changes is a core interest of molecular, ecological, neuropsychological and evolutionary biologists. Further, we will take a journey through different disciplines, such as social sciences, finance and computational gastronomy, to present commonalities and differences in how networks change and can be analysed. Finally, we envision how borrowing ideas from these disciplines could enrich the future of life science research.

A genome-scale metabolic model of potato late blight suggests a photosynthesis suppression mechanism

BACKGROUND

Phytophthora infestans is a plant pathogen that causes an important plant disease known as late blight in potato plants (Solanum tuberosum) and several other solanaceous hosts. This disease is the main factor affecting potato crop production worldwide. In spite of the importance of the disease, the molecular mechanisms underlying the compatibility between the pathogen and its hosts are still unknown.

RESULTS

To explain the metabolic response of late blight, specifically photosynthesis inhibition in infected plants, we reconstructed a genome-scale metabolic network of the S. tuberosum leaf, PstM1. This metabolic network simulates the effect of this disease in the leaf metabolism. PstM1 accounts for 2751 genes, 1113 metabolic functions, 1773 gene-protein-reaction associations and 1938 metabolites involved in 2072 reactions. The optimization of the model for biomass synthesis maximization in three infection time points suggested a suppression of the photosynthetic capacity related to the decrease of metabolic flux in light reactions and carbon fixation reactions. In addition, a variation pattern in the flux of carboxylation to oxygenation reactions catalyzed by RuBisCO was also identified, likely to be associated to a defense response in the compatible interaction between P. infestans and S. tuberosum.

CONCLUSIONS

In this work, we introduced simultaneously the first metabolic network of S. tuberosum and the first genome-scale metabolic model of the compatible interaction of a plant with P. infestans.

Unlocking conserved and diverged metabolic characteristics in cassava carbon assimilation via comparative genomics approach

Globally, cassava is an important source of starch, which is synthesized through carbon assimilation in cellular metabolism whereby harvested atmospheric carbon is assimilated into macromolecules. Although the carbon assimilation pathway is highly conserved across species, metabolic phenotypes could differ in composition, type, and quantity. To unravel the metabolic complexity and advantage of cassava over other starch crops, in terms of starch production, we investigated the carbon assimilation mechanisms in cassava through genome-based pathway reconstruction and comparative network analysis. First, MeRecon - the carbon assimilation pathway of cassava was reconstructed based upon six plant templates: Arabidopsis, rice, maize, castor bean, potato, and turnip. MeRecon, available at http://bml.sbi.kmutt.ac.th/MeRecon, comprises 259 reactions (199 EC numbers), 1,052 proteins (870 genes) and 259 metabolites in eight sub-metabolisms. Analysis of MeRecon and the carbon assimilation pathways of the plant templates revealed the overall topology is highly conserved, but variations at sub metabolism level were found in relation to complexity underlying each biochemical reaction, such as numbers of responsible enzymatic proteins and their evolved functions, which likely explain the distinct metabolic phenotype. Thus, this study provides insights into the network characteristics and mechanisms that regulate the synthesis of metabolic phenotypes of cassava.

Finite-time parametric identification for the model representing the metabolic and genetic regulatory effects of sequential aerobic respiration and anaerobic fermentation processes in Escherichia coli

Mathematical modelling applied to biological systems allows for the inferring of changes in the dynamic behaviour of organisms associated with variations in the environment. Models based on ordinary differential equations are most commonly used because of their ability to describe the mechanisms of biological systems such as transcription. The disadvantage of using this approach is that there is a large number of parameters involved and that it is difficult to obtain them experimentally. This study presents an algorithm to obtain a finite-time parameter characterization of the model used to describe changes in the metabolic behaviour of Escherichia coli associated with environmental changes. In this scheme, super-twisting algorithm was proposed to recover the derivative of all the proteins and mRNA of E. coli associated to changes in the concentration of oxygen available in the growth media. The 75 identified parameters in this study maintain the biological coherence of the system and they were estimated with no more than 20% error with respect to the real ones included in the proposed model.

Metabolomics analysis: Finding out metabolic building blocks

In this paper we propose a new methodology for the analysis of metabolic networks. We use the notion of strongly connected components of a graph, called in this context metabolic building blocks. Every strongly connected component is contracted to a single node in such a way that the resulting graph is a directed acyclic graph, called a metabolic DAG, with a considerably reduced number of nodes. The property of being a directed acyclic graph brings out a background graph topology that reveals the connectivity of the metabolic network, as well as bridges, isolated nodes and cut nodes. Altogether, it becomes a key information for the discovery of functional metabolic relations. Our methodology has been applied to the glycolysis and the purine metabolic pathways for all organisms in the KEGG database, although it is general enough to work on any database. As expected, using the metabolic DAGs formalism, a considerable reduction on the size of the metabolic networks has been obtained, specially in the case of the purine pathway due to its relative larger size. As a proof of concept, from the information captured by a metabolic DAG and its corresponding metabolic building blocks, we obtain the core of the glycolysis pathway and the core of the purine metabolism pathway and detect some essential metabolic building blocks that reveal the key reactions in both pathways. Finally, the application of our methodology to the glycolysis pathway and the purine metabolism pathway reproduce the tree of life for the whole set of the organisms represented in the KEGG database which supports the utility of this research.

Integrated omics analysis of specialized metabolism in medicinal plants

Medicinal plants are a rich source of highly diverse specialized metabolites with important pharmacological properties. Until recently, plant biologists were limited in their ability to explore the biosynthetic pathways of these metabolites, mainly due to the scarcity of plant genomics resources. However, recent advances in high-throughput large-scale analytical methods have enabled plant biologists to discover biosynthetic pathways for important plant-based medicinal metabolites. The reduced cost of generating omics datasets and the development of computational tools for their analysis and integration have led to the elucidation of biosynthetic pathways of several bioactive metabolites of plant origin. These discoveries have inspired synthetic biology approaches to develop microbial systems to produce bioactive metabolites originating from plants, an alternative sustainable source of medicinally important chemicals. Since the demand for medicinal compounds are increasing with the world's population, understanding the complete biosynthesis of specialized metabolites becomes important to identify or develop reliable sources in the future. Here, we review the contributions of major omics approaches and their integration to our understanding of the biosynthetic pathways of bioactive metabolites. We briefly discuss different approaches for integrating omics datasets to extract biologically relevant knowledge and the application of omics datasets in the construction and reconstruction of metabolic models.

Genomics of lactic acid bacteria: Current status and potential applications

Lactic acid bacteria (LAB) are widely used for the production of a variety of foods and feed raw materials where they contribute to flavor and texture of the fermented products. In addition, specific LAB strains are considered as probiotic due to their health-promoting effects in consumers. Recently, the genome sequencing of LAB is booming and the increased amount of published genomics data brings unprecedented opportunity for us to reveal the important traits of LAB. This review describes the recent progress on LAB genomics and special emphasis is placed on understanding the industry-related physiological features based on genomics analysis. Moreover, strategies to engineer metabolic capacity and stress tolerance of LAB with improved industrial performance are also discussed.

A comparative review of computational methods for pathway perturbation analysis: dynamical and topological perspectives

Stem cells offer great promise within the field of regenerative medicine but despite encouraging results, the large scale use of stem cells for therapeutic applications still faces challenges when it comes to controlling signaling pathway responses with respect to environmental perturbations.

CompNet: a GUI based tool for comparison of multiple biological interaction networks

Background

Network visualization and analysis tools aid in better understanding of complex biological systems. Furthermore, to understand the differences in behaviour of system(s) under various environmental conditions (e.g. stress, infection), comparing multiple networks becomes necessary. Such comparisons between multiple networks may help in asserting causation and in identifying key components of the studied biological system(s). Although many available network comparison methods exist, which employ techniques like network alignment and querying to compute pair-wise similarity between selected networks, most of them have limited features with respect to interactive visual comparison of multiple networks.

Results

In this paper, we present CompNet - a graphical user interface based network comparison tool, which allows visual comparison of multiple networks based on various network metrics. CompNet allows interactive visualization of the union, intersection and/or complement regions of a selected set of networks. Different visualization features (e.g. pie-nodes, edge-pie matrix, etc.) aid in easy identification of the key nodes/interactions and their significance across the compared networks. The tool also allows one to perform network comparisons on the basis of neighbourhood architecture of constituent nodes and community compositions, a feature particularly useful while analyzing biological networks. To demonstrate the utility of CompNet, we have compared a (time-series) human gene-expression dataset, post-infection by two strains of Mycobacterium tuberculosis, overlaid on the human protein-protein interaction network. Using various functionalities of CompNet not only allowed us to comprehend changes in interaction patterns over the course of infection, but also helped in inferring the probable fates of the host cells upon infection by the two strains.

Conclusions

CompNet is expected to be a valuable visual data mining tool and is freely available for academic use from http://metagenomics.atc.tcs.com/compnet/ or http://121.241.184.233/compnet/

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1013-x) contains supplementary material, which is available to authorized users.

Agent-based spatiotemporal simulation of biomolecular systems within the open source MASON framework

Agent-based modelling is being used to represent biological systems with increasing frequency and success. This paper presents the implementation of a new tool for biomolecular reaction modelling in the open source Multiagent Simulator of Neighborhoods framework. The rationale behind this new tool is the necessity to describe interactions at the molecular level to be able to grasp emergent and meaningful biological behaviour. We are particularly interested in characterising and quantifying the various effects that facilitate biocatalysis. Enzymes may display high specificity for their substrates and this information is crucial to the engineering and optimisation of bioprocesses. Simulation results demonstrate that molecule distributions, reaction rate parameters, and structural parameters can be adjusted separately in the simulation allowing a comprehensive study of individual effects in the context of realistic cell environments. While higher percentage of collisions with occurrence of reaction increases the affinity of the enzyme to the substrate, a faster reaction (i.e., turnover number) leads to a smaller number of time steps. Slower diffusion rates and molecular crowding (physical hurdles) decrease the collision rate of reactants, hence reducing the reaction rate, as expected. Also, the random distribution of molecules affects the results significantly.

Double and multiple knockout simulations for genome-scale metabolic network reconstructions

Background

Constraint-based modeling of genome-scale metabolic network reconstructions has become a widely used approach in computational biology. Flux coupling analysis is a constraint-based method that analyses the impact of single reaction knockouts on other reactions in the network.

Results

We present an extension of flux coupling analysis for double and multiple gene or reaction knockouts, and develop corresponding algorithms for an in silico simulation. To evaluate our method, we perform a full single and double knockout analysis on a selection of genome-scale metabolic network reconstructions and compare the results.

Software

A prototype implementation of double knockout simulation is available at http://hoverboard.io/L4FC.

Multi -omics and metabolic modelling pipelines: challenges and tools for systems microbiology

Integrated -omics approaches are quickly spreading across microbiology research labs, leading to (i) the possibility of detecting previously hidden features of microbial cells like multi-scale spatial organization and (ii) tracing molecular components across multiple cellular functional states. This promises to reduce the knowledge gap between genotype and phenotype and poses new challenges for computational microbiologists. We underline how the capability to unravel the complexity of microbial life will strongly depend on the integration of the huge and diverse amount of information that can be derived today from -omics experiments. In this work, we present opportunities and challenges of multi -omics data integration in current systems biology pipelines. We here discuss which layers of biological information are important for biotechnological and clinical purposes, with a special focus on bacterial metabolism and modelling procedures. A general review of the most recent computational tools for performing large-scale datasets integration is also presented, together with a possible framework to guide the design of systems biology experiments by microbiologists.

Predictive sulfur metabolism - a field in flux

The key role of sulfur metabolites in response to biotic and abiotic stress in plants, as well as their importance in diet and health has led to a significant interest and effort in trying to understand and manipulate the production of relevant compounds. Metabolic engineering utilizes a set of theoretical tools to help rationally design modifications that enhance the production of a desired metabolite. Such approaches have proven their value in bacterial systems, however, the paucity of success stories to date in plants, suggests that challenges remain. Here, we review the most commonly used methods for understanding metabolic flux, focusing on the sulfur assimilatory pathway. We highlight known issues with both experimental and theoretical approaches, as well as presenting recent methods for integrating different modeling strategies, and progress toward an understanding of flux at the whole plant level.

A genome-scale metabolic flux model of Escherichia coli K-12 derived from the EcoCyc database

BACKGROUND

Constraint-based models of Escherichia coli metabolic flux have played a key role in computational studies of cellular metabolism at the genome scale. We sought to develop a next-generation constraint-based E. coli model that achieved improved phenotypic prediction accuracy while being frequently updated and easy to use. We also sought to compare model predictions with experimental data to highlight open questions in E. coli biology.

RESULTS

We present EcoCyc-18.0-GEM, a genome-scale model of the E. coli K-12 MG1655 metabolic network. The model is automatically generated from the current state of EcoCyc using the MetaFlux software, enabling the release of multiple model updates per year. EcoCyc-18.0-GEM encompasses 1445 genes, 2286 unique metabolic reactions, and 1453 unique metabolites. We demonstrate a three-part validation of the model that breaks new ground in breadth and accuracy: (i) Comparison of simulated growth in aerobic and anaerobic glucose culture with experimental results from chemostat culture and simulation results from the E. coli modeling literature. (ii) Essentiality prediction for the 1445 genes represented in the model, in which EcoCyc-18.0-GEM achieves an improved accuracy of 95.2% in predicting the growth phenotype of experimental gene knockouts. (iii) Nutrient utilization predictions under 431 different media conditions, for which the model achieves an overall accuracy of 80.7%. The model's derivation from EcoCyc enables query and visualization via the EcoCyc website, facilitating model reuse and validation by inspection. We present an extensive investigation of disagreements between EcoCyc-18.0-GEM predictions and experimental data to highlight areas of interest to E. coli modelers and experimentalists, including 70 incorrect predictions of gene essentiality on glucose, 80 incorrect predictions of gene essentiality on glycerol, and 83 incorrect predictions of nutrient utilization.

CONCLUSION

Significant advantages can be derived from the combination of model organism databases and flux balance modeling represented by MetaFlux. Interpretation of the EcoCyc database as a flux balance model results in a highly accurate metabolic model and provides a rigorous consistency check for information stored in the database.

Computer-aided design for metabolic engineering

The development and application of biotechnology-based strategies has had a great socio-economical impact and is likely to play a crucial role in the foundation of more sustainable and efficient industrial processes. Within biotechnology, metabolic engineering aims at the directed improvement of cellular properties, often with the goal of synthesizing a target chemical compound. The use of computer-aided design (CAD) tools, along with the continuously emerging advanced genetic engineering techniques have allowed metabolic engineering to broaden and streamline the process of heterologous compound-production. In this work, we review the CAD tools available for metabolic engineering with an emphasis, on retrosynthesis methodologies. Recent advances in genetic engineering strategies for pathway implementation and optimization are also reviewed as well as a range of bionalytical tools to validate in silico predictions. A case study applying retrosynthesis is presented as an experimental verification of the output from Retropath, the first complete automated computational pipeline applicable to metabolic engineering. Applying this CAD pipeline, together with genetic reassembly and optimization of culture conditions led to improved production of the plant flavonoid pinocembrin. Coupling CAD tools with advanced genetic engineering strategies and bioprocess optimization is crucial for enhanced product yields and will be of great value for the development of non-natural products through sustainable biotechnological processes.

Network analysis: a new approach to study endocrine disorders

Systems biology is the study of the interactions that occur between the components of individual cells - including genes, proteins, transcription factors, small molecules, and metabolites, and their relationships to complex physiological and pathological processes. The application of systems biology to medicine promises rapid advances in both our understanding of disease and the development of novel treatment options. Network biology has emerged as the primary tool for studying systems biology as it utilises the mathematical analysis of the relationships between connected objects in a biological system and allows the integration of varied 'omic' datasets (including genomics, metabolomics, proteomics, etc.). Analysis of network biology generates interactome models to infer and assess function; to understand mechanisms, and to prioritise candidates for further investigation. This review provides an overview of network methods used to support this research and an insight into current applications of network analysis applied to endocrinology. A wide spectrum of endocrine disorders are included ranging from congenital hyperinsulinism in infancy, through childhood developmental and growth disorders, to the development of metabolic diseases in early and late adulthood, such as obesity and obesity-related pathologies. In addition to providing a deeper understanding of diseases processes, network biology is also central to the development of personalised treatment strategies which will integrate pharmacogenomics with systems biology of the individual.

FastPros: screening of reaction knockout strategies for metabolic engineering

Motivation: Although constraint-based flux analysis of knockout strains has facilitated the production of desirable metabolites in microbes, current screening methods have placed a limitation on the number knockouts that can be simultaneously analyzed.

Results: Here, we propose a novel screening method named FastPros. In this method, the potential of a given reaction knockout for production of a specific metabolite is evaluated by shadow pricing of the constraint in the flux balance analysis, which generates a screening score to obtain candidate knockout sets. To evaluate the performance of FastPros, we screened knockout sets to produce each metabolite in the entire Escherichia coli metabolic network. We found that 75% of these metabolites could be produced under biomass maximization conditions by adding up to 25 reaction knockouts. Furthermore, we demonstrated that using FastPros in tandem with another screening method, OptKnock, could further improve target metabolite productivity.

Availability and implementation: Source code is freely available at http://www-shimizu.ist.osaka-u.ac.jp/shimizu_lab/FastPros/, implemented in MATLAB and COBRA toolbox.

Contact:chikara.furusawa@riken.jp or shimizu@ist.osaka-u.ac.jp

Supplementary information:Supplementary data are available at Bioinformatics online.