Jole Costanza

Robust design of microbial strains

MOTIVATION Metabolic engineering algorithms provide means to optimize a biological process leading to the improvement of a biotechnological interesting molecule. Therefore, it is important to understand how to act in a metabolic pathway in order to have the best results in terms of productions. In this work, we present a computational framework that searches for optimal and robust microbial strains that are able to produce target molecules. Our framework performs three tasks: it evaluates the parameter sensitivity of the microbial model, searches for the optimal genetic or fluxes design and finally calculates the robustness of the microbial strains. We are capable to combine the exploration of species, reactions, pathways and knockout parameter spaces with the Pareto-optimality principle. RESULTS Our framework provides also theoretical and practical guidelines for design automation. The statistical cross comparison of our new optimization procedures, performed with respect to currently widely used algorithms for bacteria (e.g. Escherichia coli) over different multiple functions, reveals good performances over a variety of biotechnological products. AVAILABILITY http://www.dmi.unict.it/nicosia/pathDesign.html. CONTACT nicosia@dmi.unict.it or pl219@cam.ac.uk SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Pareto Optimality in Organelle Energy Metabolism Analysis

In low and high eukaryotes, energy is collected or transformed in compartments, the organelles. The rich variety of size, characteristics, and density of the organelles makes it difficult to build a general picture. In this paper, we make use of the Pareto-front analysis to investigate the optimization of energy metabolism in mitochondria and chloroplasts. Using the Pareto optimality principle, we compare models of organelle metabolism on the basis of single- and multiobjective optimization, approximation techniques (the Bayesian Automatic Relevance Determination), robustness, and pathway sensitivity analysis. Finally, we report the first analysis of the metabolic model for the hydrogenosome of Trichomonas vaginalis, which is found in several protozoan parasites. Our analysis has shown the importance of the Pareto optimality for such comparison and for insights into the evolution of the metabolism from cytoplasmic to organelle bound, involving a model order reduction. We report that Pareto fronts represent an asymptotic analysis useful to describe the metabolism of an organism aimed at maximizing concurrently two or more metabolite concentrations.

Pareto Optimal Design for Synthetic Biology

Recent advances in synthetic biology call for robust, flexible and efficient in silico optimization methodologies. We present a Pareto design approach for the bi-level optimization problem associated to the overproduction of specific metabolites in Escherichia coli. Our method efficiently explores the high dimensional genetic manipulation space, finding a number of trade-offs between synthetic and biological objectives, hence furnishing a deeper biological insight to the addressed problem and important results for industrial purposes. We demonstrate the computational capabilities of our Pareto-oriented approach comparing it with state-of-the-art heuristics in the overproduction problems of i) 1,4-butanediol, ii) myristoyl-CoA, i ii) malonyl-CoA , iv) acetate and v) succinate. We show that our algorithms are able to gracefully adapt and scale to more complex models and more biologically-relevant simulations of the genetic manipulations allowed. The Results obtained for 1,4-butanediol overproduction significantly outperform results previously obtained, in terms of 1,4-butanediol to biomass formation ratio and knock-out costs. In particular overproduction percentage is of +662.7%, from 1.425 mmolh-1gDW-1 (wild type) to 10.869 mmolh-1gDW-1, with a knockout cost of 6. Whereas, Pareto-optimal designs we have found in fatty acid optimizations strictly dominate the ones obtained by the other methodologies, e.g., biomass and myristoyl-CoA exportation improvement of +21.43% (0.17 h-1) and +5.19% (1.62 mmolh-1gDW-1), respectively. Furthermore CPU time required by our heuristic approach is more than halved. Finally we implement pathway oriented sensitivity analysis, epsilon-dominance analysis and robustness analysis to enhance our biological understanding of the problem and to improve the optimization algorithm capabilities.

Efficient Behavior of Photosynthetic Organelles via Pareto Optimality, Identifiability, and Sensitivity Analysis

In this work, we develop methodologies for analyzing and cross comparing metabolic models. We investigate three important metabolic networks to discuss the complexity of biological organization of organisms, modeling, and system properties. In particular, we analyze these metabolic networks because of their biotechnological and basic science importance: the photosynthetic carbon metabolism in a general leaf, the Rhodobacter spheroides bacterium, and the Chlamydomonas reinhardtii alga. We adopt single- and multi-objective optimization algorithms to maximize the CO 2 uptake rate and the production of metabolites of industrial interest or for ecological purposes. We focus both on the level of genes (e.g., finding genetic manipulations to increase the production of one or more metabolites) and on finding concentration enzymes for improving the CO 2 consumption. We find that R. spheroides is able to absorb an amount of CO 2 until 57.452 mmol h (-1) gDW (-1) , while C. reinhardtii obtains a maximum of 6.7331. We report that the Pareto front analysis proves extremely useful to compare different organisms, as well as providing the possibility to investigate them with the same framework. By using the sensitivity and robustness analysis, our framework identifies the most sensitive and fragile components of the biological systems we take into account, allowing us to compare their models. We adopt the identifiability analysis to detect functional relations among enzymes; we observe that RuBisCO, GAPDH, and FBPase belong to the same functional group, as suggested also by the sensitivity analysis.

Multi-Target Analysis and Design of Mitochondrial Metabolism.

Analyzing and optimizing biological models is often identified as a research priority in biomedical engineering. An important feature of a model should be the ability to find the best condition in which an organism has to be grown in order to reach specific optimal output values chosen by the researcher. In this work, we take into account a mitochondrial model analyzed with flux-balance analysis. The optimal design and assessment of these models is achieved through single- and/or multi-objective optimization techniques driven by epsilon-dominance and identifiability analysis. Our optimization algorithm searches for the values of the flux rates that optimize multiple cellular functions simultaneously. The optimization of the fluxes of the metabolic network includes not only input fluxes, but also internal fluxes. A faster convergence process with robust candidate solutions is permitted by a relaxed Pareto dominance, regulating the granularity of the approximation of the desired Pareto front. We find that the maximum ATP production is linked to a total consumption of NADH, and reaching the maximum amount of NADH leads to an increasing request of NADH from the external environment. Furthermore, the identifiability analysis characterizes the type and the stage of three monogenic diseases. Finally, we propose a new methodology to extend any constraint-based model using protein abundances.

Identification of Sensitive Enzymes in the Photosynthetic Carbon Metabolism

Understanding and optimizing the CO(2) fixation process would allow human beings to address better current energy and biotechnology issues. We focused on modeling the C(3) photosynthetic Carbon metabolism pathway with the aim of identifying the minimal set of enzymes whose biotechnological alteration could allow a functional re-engineering of the pathway. To achieve this result we merged in a single powerful pipe-line Sensitivity Analysis (SA), Single- (SO) and Multi-Objective Optimization (MO), and Robustness Analysis (RA). By using our recently developed multipurpose optimization algorithms (PAO and PMO2) here we extend our work exploring a large combinatorial solution space and most importantly, here we present an important reduction of the problem search space. From the initial number of 23 enzymes we have identified 11 enzymes whose targeting in the C(3) photosynthetic Carbon metabolism would provide about 90% of the overall functional optimization. Both in terms of maximal CO(2) Uptake and minimal Nitrogen consumption, these 11 sensitive enzymes are confirmed to play a key role. Finally we present a RA to confirm our findings.

LowMACA: exploiting protein family analysis for the identification of rare driver mutations in cancer

BackgroundThe increasing availability of resequencing data has led to a better understanding of the most important genes in cancer development. Nevertheless, the mutational landscape of many tumor types is heterogeneous and encompasses a long tail of potential driver genes that are systematically excluded by currently available methods due to the low frequency of their mutations. We developed LowMACA (Low frequency Mutations Analysis via Consensus Alignment), a method that combines the mutations of various proteins sharing the same functional domains to identify conserved residues that harbor clustered mutations in multiple sequence alignments. LowMACA is designed to visualize and statistically assess potential driver genes through the identification of their mutational hotspots.ResultsWe analyzed the Ras superfamily exploiting the known driver mutations of the trio K-N-HRAS, identifying new putative driver mutations and genes belonging to less known members of the Rho, Rab and Rheb subfamilies. Furthermore, we applied the same concept to a list of known and candidate driver genes, and observed that low confidence genes show similar patterns of mutation compared to high confidence genes of the same protein family.ConclusionsLowMACA is a software for the identification of gain-of-function mutations in putative oncogenic families, increasing the amount of information on functional domains and their possible role in cancer. In this context LowMACA emphasizes the role of genes mutated at low frequency otherwise undetectable by classical single gene analysis.LowMACA is an R package available at http://www.bioconductor.org/packages/release/bioc/html/LowMACA.html. It is also available as a GUI standalone downloadable at: https://cgsb.genomics.iit.it/wiki/projects/LowMACA

Analysis and design of molecular machines

Abstract Biologically inspired computation has been recently used with mathematical models towards the design of new synthetic organisms. In this work, we use Pareto optimality to optimize these organisms in a multi-objective fashion. We infer the best knockout strategies to perform specific tasks in bacteria, which involve concurrent maximization/minimization of multiple functions (codomain) and optimization of several decision variables (domain). Furthermore, we propose and exploit a mapping between the metabolism and a register machine. We show that optimized bacteria have computational capability and act as molecular Turing machines programmed using a Pareto optimal solution. Finally, we investigate communication between bacteria as a means to evaluate their computational capability. We report that the density and gradient of the Pareto curve are useful tools to compare models and understand their structure, while modelling organisms as computers proves useful to carry out computation using biological machines with specific input–output conditions, as well as to estimate the bacterial computational effort for specific tasks.

Pareto epsilon-dominance and identifiable solutions for BioCAD modeling

We propose a framework to design metabolic pathways in which many objectives are optimized simultaneously. This allows to characterize the energy signature in models of algal and mitochondrial metabolism. The optimal design and assessment of the model is achieved through a multi-objective optimization technique driven by epsilon-dominance and identifiability analysis. A faster convergence process with robust candidate solutions is permitted by a relaxed Pareto dominance, regulating the granularity of the approximation of the Pareto front. Our framework is also suitable for black-box analysis, enabling to investigate and optimize any biological pathway modeled with ODEs, DAEs, FBA and GPR.

Multi-objective optimisation, sensitivity and robustness analysis in FBA modelling

In this work, we propose a computational framework to design in silico robust bacteria able to overproduce multiple metabolites. To this end, we search the optimal genetic manipulations, in terms of knockout, which also guarantee the growth of the organism. We introduce a multi-objective optimisation algorithm, called Genetic Design through Multi-Objective (GDMO), and test it in several organisms to maximise the production of key intermediate metabolites such as succinate and acetate. We obtain a vast set of Pareto optimal solutions; each of them represents an organism strain. For each solution, we evaluate the fragility by calculating three robustness indexes and by exploring reactions and metabolite interactions. Finally, we perform the Sensitivity Analysis of the metabolic model, which finds the inputs with the highest influence on the outputs of the model. We show that our methodology provides effective vision of the achievable synthetic strain landscape and a powerful design pipeline.