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.

Computer-aided design and implementation of efficient biosynthetic pathways to produce high added-value products derived from tyrosine in Escherichia coli.

Developing efficient bioprocesses requires selecting the best biosynthetic pathways, which can be challenging and time-consuming due to the vast amount of data available in databases and literature. The extension of the shikimate pathway for the biosynthesis of commercially attractive molecules often involves promiscuous enzymes or lacks well-established routes. To address these challenges, we developed a computational workflow integrating enumeration/retrosynthesis algorithms, a toolbox for pathway analysis, enzyme selection tools, and a gene discovery pipeline, supported by manual curation and literature review. Our focus has been on implementing biosynthetic pathways for tyrosine-derived compounds, specifically L-3,4-dihydroxyphenylalanine (L-DOPA) and dopamine, with significant applications in health and nutrition. We selected one pathway to produce L-DOPA and two different pathways for dopamine-one already described in the literature and a novel pathway. Our goal was either to identify the most suitable gene candidates for expression in Escherichia coli for the known pathways or to discover innovative pathways. Although not all implemented pathways resulted in the accumulation of target compounds, in our shake-flask experiments we achieved a maximum L-DOPA titer of 0.71 g/L and dopamine titers of 0.29 and 0.21 g/L for known and novel pathways, respectively. In the case of L-DOPA, we utilized, for the first time, a mutant version of tyrosinase from Ralstonia solanacearum. Production of dopamine via the known biosynthesis route was accomplished by coupling the L-DOPA pathway with the expression of DOPA decarboxylase from Pseudomonas putida, resulting in a unique biosynthetic pathway never reported in literature before. In the context of the novel pathway, dopamine was produced using tyramine as the intermediate compound. To achieve this, tyrosine was initially converted into tyramine by expressing TDC from Levilactobacillus brevis, which, in turn, was converted into dopamine through the action of the enzyme encoded by ppoMP from Mucuna pruriens. This marks the first time that an alternative biosynthetic pathway for dopamine has been validated in microbes. These findings underscore the effectiveness of our computational workflow in facilitating pathway enumeration and selection, offering the potential to uncover novel biosynthetic routes, thus paving the way for other target compounds of biotechnological interest.

Anterior minimally invasive bridge-plate technique for treatment of humeral shaft nonunion.

BACKGROUND: The present study introduces a new surgical technique and the results of a case series of patients with humeral shaft nonunion. MATERIALS AND METHODS: Fifteen patients with diagnosis of diaphyseal nonunion of humerus were operated by a bridge-plate technique. A 4.5-mm plate is slid on the anterior surface of the humerus, submuscular to the brachial muscle. With the plate over the anterior surface of the humerus, screws are inserted from anterior to posterior on the ends of the plate. When there is a small bone gap, an iliac autologous graft is inserted. Minimum follow-up was 1 year. RESULTS: Bone healing was obtained in all patients: 1.5 months postoperatively in 11 patients, 2 months in 3 patients, and 3 months in 1 patient. There were no postoperative infections, there was one case with loosening of the screws and plate, and there were no nerve injuries. CONCLUSIONS: The present technique avoids wide dissection, radial nerve isolation, and periosteum stripping. The anterior minimally invasive bridge-plate technique for treatment of humeral shaft nonunion is a safe procedure and obtained bone healing in all patients in this series.

Turnover Dependent Phenotypic Simulation: A Quantitative Constraint-Based Simulation Method That Accommodates All Main Strain Design Strategies.

The uncertain relationship between genotype and phenotype can make strain engineering an arduous trial and error process. To identify promising gene targets faster, constraint-based modeling methodologies are often used, although they remain limited in their predictive power. Even though the search for gene knockouts is fairly established in constraint-based modeling, most strain design methods still model gene up/down-regulations by forcing the corresponding flux values to fixed levels without taking in consideration the availability of resources. Here, we present a constraint-based algorithm, the turnover dependent phenotypic simulation (TDPS) that quantitatively simulates phenotypes in a resource conscious manner. Unlike other available algorithms, TDPS does not force flux values and considers resource availability, using metabolite production turnovers as an indicator of metabolite abundance. TDPS can simulate up-regulation of metabolic reactions as well as the introduction of heterologous genes, alongside gene deletion and down-regulation scenarios. TDPS simulations were validated using engineered Saccharomyces cerevisiae strains available in the literature by comparing the simulated and experimental production yields of the target metabolite. For many of the strains evaluated, the experimental production yields were within the simulated intervals and the relative strain performance could be predicted with TDPS. However, the algorithm failed to predict some of the production changes observed experimentally, suggesting that further improvements are necessary. The results also showed that TDPS may be helpful in finding metabolic bottlenecks, but further experiments would be required to confirm these findings.

Discovery and implementation of a novel pathway for n-butanol production via 2-oxoglutarate.

BACKGROUND: One of the European Union directives indicates that 10% of all fuels must be bio-synthesized by 2020. In this regard, biobutanol-natively produced by clostridial strains-poses as a promising alternative biofuel. One possible approach to overcome the difficulties of the industrial exploration of the native producers is the expression of more suitable pathways in robust microorganisms such as Escherichia coli. The enumeration of novel pathways is a powerful tool, allowing to identify non-obvious combinations of enzymes to produce a target compound. RESULTS: This work describes the in silico driven design of E. coli strains able to produce butanol via 2-oxoglutarate by a novel pathway. This butanol pathway was generated by a hypergraph algorithm and selected from an initial set of 105,954 different routes by successively applying different filters, such as stoichiometric feasibility, size and novelty. The implementation of this pathway involved seven catalytic steps and required the insertion of nine heterologous genes from various sources in E. coli distributed in three plasmids. Expressing butanol genes in E. coli K12 and cultivation in High-Density Medium formulation seem to favor butanol accumulation via the 2-oxoglutarate pathway. The maximum butanol titer obtained was 85 ± 1 mg L-1 by cultivating the cells in bioreactors. CONCLUSIONS: In this work, we were able to successfully translate the computational analysis into in vivo applications, designing novel strains of E. coli able to produce n-butanol via an innovative pathway. Our results demonstrate that enumeration algorithms can broad the spectrum of butanol producing pathways. This validation encourages further research to other target compounds.

Analyzing and Designing Cell Factories with OptFlux.

OptFlux was launched in 2010 as the first open-source and user-friendly platform containing all the major methods for performing metabolic engineering tasks in silico. Main features included the possibility of performing microbial strain simulations with widely used methods such as Flux Balance Analysis and strain design using Evolutionary Algorithms. Since then, OptFlux suffered a major re-factoring to improve its efficiency and reliability, while many features were added in the form of novel plug-ins, such as the BioVisualizer and the over/under expression plug-ins. The current chapter described the main mathematical formulations of the major methods implemented within OptFlux, also providing a detailed guide on the usage of those functionalities.

Reconstruction of a genome-scale metabolic model for Actinobacillus succinogenes 130Z.

BACKGROUND: Actinobacillus succinogenes is a promising bacterial catalyst for the bioproduction of succinic acid from low-cost raw materials. In this work, a genome-scale metabolic model was reconstructed and used to assess the metabolic capabilities of this microorganism under producing conditions. RESULTS: The model, iBP722, was reconstructed based on the functional reannotation of the complete genome sequence of A. succinogenes 130Z and manual inspection of metabolic pathways, covering 1072 enzymatic reactions associated with 722 metabolic genes that involve 713 metabolites. The highly curated model was effective in capturing the growth of A. succinogenes on various carbon sources, as well as the SA production under various growth conditions with fair agreement between experimental and predicted data. Calculated flux distributions under different conditions show that a number of metabolic pathways are affected by the activity of some metabolic enzymes at key nodes in metabolism, including the transport mechanism of carbon sources and the ability to fix carbon dioxide. CONCLUSIONS: The established genome-scale metabolic model can be used for model-driven strain design and medium alteration to improve succinic acid yields.

Genome-Wide Semi-Automated Annotation of Transporter Systems.

Usually, transport reactions are added to genome-scale metabolic models (GSMMs) based on experimental data and literature. This approach does not allow associating specific genes with transport reactions, which impairs the ability of the model to predict effects of gene deletions. Novel methods for systematic genome-wide transporter functional annotation and their integration into GSMMs are therefore necessary. In this work, an automatic system to detect and classify all potential membrane transport proteins for a given genome and integrate the related reactions into GSMMs is proposed, based on the identification and classification of genes that encode transmembrane proteins. The Transport Reactions Annotation and Generation (TRIAGE) tool identifies the metabolites transported by each transmembrane protein and its transporter family. The localization of the carriers is also predicted and, consequently, their action is confined to a given membrane. The integration of the data provided by TRIAGE with highly curated models allowed the identification of new transport reactions. TRIAGE is included in the new release of merlin, a software tool previously developed by the authors, which expedites the GSMM reconstruction processes.

CBFA: phenotype prediction integrating metabolic models with constraints derived from experimental data.

BACKGROUND: Flux analysis methods lie at the core of Metabolic Engineering (ME), providing methods for phenotype simulation that allow the determination of flux distributions under different conditions. Although many constraint-based modeling software tools have been developed and published, none provides a free user-friendly application that makes available the full portfolio of flux analysis methods. RESULTS: This work presents Constraint-based Flux Analysis (CBFA), an open-source software application for flux analysis in metabolic models that implements several methods for phenotype prediction, allowing users to define constraints associated with measured fluxes and/or flux ratios, together with environmental conditions (e.g. media) and reaction/gene knockouts. CBFA identifies the set of applicable methods based on the constraints defined from user inputs, encompassing algebraic and constraint-based simulation methods. The integration of CBFA within the OptFlux framework for ME enables the utilization of different model formats and standards and the integration with complementary methods for phenotype simulation and visualization of results. CONCLUSIONS: A general-purpose and flexible application is proposed that is independent of the origin of the constraints defined for a given simulation. The aim is to provide a simple to use software tool focused on the application of several flux prediction methods.

An integrated computational environment for elementary modes analysis of biochemical networks.

Elementary flux modes (EFMs) have been claimed as one of the most promising approaches for pathway analysis. These are a set of vectors that emerge from the stoichiometric matrix of a biochemical network through the use of convex analysis. The computation of all EFMs of a given network is an NP-hard problem and existing algorithms do not scale well. Moreover, the analysis of results is difficult given the thousands or millions of possible modes generated. In this work, we propose a new plug-in, running on top of the OptFlux Metabolic Engineering workbench (Rocha et al., 2010), whose aims are to ease the analysis of these results and explore synergies among EFM analysis, phenotype simulation and strain optimisation. Two case studies are shown to illustrate the capabilities of the proposed tool.

A computational tool for the simulation and optimization of microbial strains accounting integrated metabolic/regulatory information.

BACKGROUND AND SCOPE: Recently, a number of methods and tools have been proposed to allow the use of genome-scale metabolic models for the phenotype simulation and optimization of microbial strains, within the field of Metabolic Engineering (ME). One of the limitations of most of these algorithms and tools is the fact that only metabolic information is taken into account, disregarding knowledge on regulatory events. IMPLEMENTATION AND PERFORMANCES: This work proposes a novel software tool that implements methods for the phenotype simulation and optimization of microbial strains using integrated models, encompassing both metabolic and regulatory information. This tool is developed as a plug-in that runs over OptFlux, a computational platform that aims to be a reference tool for the ME community. AVAILABILITY: The plug-in is made available in the OptFlux web site (www.optflux.org) together with examples and documentation.

OptFlux: an open-source software platform for in silico metabolic engineering.

BACKGROUND: Over the last few years a number of methods have been proposed for the phenotype simulation of microorganisms under different environmental and genetic conditions. These have been used as the basis to support the discovery of successful genetic modifications of the microbial metabolism to address industrial goals. However, the use of these methods has been restricted to bioinformaticians or other expert researchers. The main aim of this work is, therefore, to provide a user-friendly computational tool for Metabolic Engineering applications. RESULTS: OptFlux is an open-source and modular software aimed at being the reference computational application in the field. It is the first tool to incorporate strain optimization tasks, i.e., the identification of Metabolic Engineering targets, using Evolutionary Algorithms/Simulated Annealing metaheuristics or the previously proposed OptKnock algorithm. It also allows the use of stoichiometric metabolic models for (i) phenotype simulation of both wild-type and mutant organisms, using the methods of Flux Balance Analysis, Minimization of Metabolic Adjustment or Regulatory on/off Minimization of Metabolic flux changes, (ii) Metabolic Flux Analysis, computing the admissible flux space given a set of measured fluxes, and (iii) pathway analysis through the calculation of Elementary Flux Modes. OptFlux also contemplates several methods for model simplification and other pre-processing operations aimed at reducing the search space for optimization algorithms. The software supports importing/exporting to several flat file formats and it is compatible with the SBML standard. OptFlux has a visualization module that allows the analysis of the model structure that is compatible with the layout information of Cell Designer, allowing the superimposition of simulation results with the model graph. CONCLUSIONS: The OptFlux software is freely available, together with documentation and other resources, thus bridging the gap from research in strain optimization algorithms and the final users. It is a valuable platform for researchers in the field that have available a number of useful tools. Its open-source nature invites contributions by all those interested in making their methods available for the community. Given its plug-in based architecture it can be extended with new functionalities. Currently, several plug-ins are being developed, including network topology analysis tools and the integration with Boolean network based regulatory models.