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

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

Multi-objective optimization of genome-scale metabolic models: the case of ethanol production

Ethanol is among the largest fermentation product used worldwide, accounting for more than 90% of all biofuel produced in the last decade. However current production methods of ethanol are unable to meet the requirements of increasing global demand, because of low yields on glucose sources. In this work, we present an in silico multi-objective optimization and analyses of eight genome-scale metabolic networks for the overproduction of ethanol within the engineered cell. We introduce MOME (multi-objective metabolic engineering) algorithm, that models both gene knockouts and enzymes up and down regulation using the Redirector framework. In a multi-step approach, MOME tackles the multi-objective optimization of biomass and ethanol production in the engineered strain; and performs genetic design and clustering analyses on the optimization results. We find in silico E. coli Pareto optimal strains with a knockout cost of 14 characterized by an ethanol production up to 19.74mmolgDW−1h−1 (+832.88% with respect to wild-type) and biomass production of 0.02h−1 (−98.06% ). The analyses on E. coli highlighted a single knockout strategy producing 16.49mmolgDW−1h−1 (+679.29% ) ethanol, with biomass equals to 0.23h−1 (−77.45% ). We also discuss results obtained by applying MOME to metabolic models of: (i) S. aureus; (ii) S. enterica; (iii) Y. pestis; (iv) S. cerevisiae; (v) C. reinhardtii; (vi) Y. lipolytica. We finally present a set of simulations in which constrains over essential genes and minimum allowable biomass were included. A bound over the maximum allowable biomass was also added, along with other settings representing rich media compositions. In the same conditions the maximum improvement in ethanol production is +195.24%

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

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

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

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

APOLO-Bari, an internet-based program for longitudinal support of bariatric surgery patients: study protocol for a randomized controlled trial

Background: Despite evidence of successful weight loss for bariatric surgery patients, some patients experience considerable weight regain over the long term. Given the strong association between post-surgery health behaviors and outcomes, aftercare intervention to address key behaviors appears to be a reasonable relapse-prevention strategy. As the burden of obesity rates increases in healthcare centers, an internet-based program appears to be a reasonable strategy for supporting bariatric surgery patients in the long term. The primary purpose of the current project is to develop and test the efficacy and perceived utility of APOLO-Bari.Methods/design: This study is a randomized control trial, which will be conducted in two hospital centers in the North of Portugal; it includes a control group receiving treatment as usual and an intervention group receiving the APOLO-Bari program for one year in addition to treatment as usual. A total of 180 male and female participants who underwent bariatric surgery (gastric sleeve or gastric bypass surgery) for 12 to 20 months will be recruited. Both groups will complete a similar set of questionnaires at baseline, every 4 months until the end of the intervention, and at 6 and 12 months follow-up. Assessment includes anthropometric variables and psychological self-report measures. The primary outcome measure will be weight regain measured at the end of treatment, and at 6 and 12 months follow-up. The secondary aims are to test the cost-effectiveness of the intervention and to investigate psychological predictors and trajectories of weight regain. APOLO-Bari was developed to address the weight regain problem in the bariatric population by offering additional guidance to bariatric patients during the postoperative period. The program includes: (a) a psychoeducational cognitive-behavioral-based self-help manual, (b) a weekly feedback messaging system that sends a feedback statement related to information reported by the participant, and (c) interactive chat sessions scheduled witThis research was partially supported by the Fundacao para a Ciencia e a Tecnologia through a European Union COMPETE program grant to Eva Conceicao (IF/01219/2014 and PTDC/MHC-PCL/4974/2012), a doctoral scholarship to Ana Pinto-Bastos (SFRH/BD/104159/2014), a doctoral scholarship to Sofia Ramalho (SFRH/BD/104182/2014), and a postdoctoral scholarship to Ana Rita Vaz (SFRH/BPD/94490/2013), co-financed by FEDER under the PT2020 Partnership Agreement (UID/PSI/01662/2013).info:eu-repo/semantics/publishedVersio

Preconcentration-enhanced electrochemical detection of paraoxon in food and environmental samples using reduced graphene oxide-modified disposable sensors

Organophosphates, such as paraoxon, are widely used as insecticides in agriculture, making their detection in environmental and food samples crucial due to their high toxicity. This study presents the development of an electrochemical sensor for the detection of paraoxon, using a screen-printed carbon electrode (SPCE) modified with electrochemically reduced graphene oxide (rGO). The modification enhanced the sensor's electrical conductivity and electrochemical performance. A novel preconcentration approach, involving potential pulses at -1.0 and 0.0 V, was employed to improve the adsorption of paraoxon on the electrode surface. Detection was performed by square wave voltammetry, and under optimized conditions, the rGO-SPCE sensor exhibited a linear range from 1.0 to 30 μmol L-1, with detection and quantification limits of 0.26 and 0.86 μmol L-1, respectively. The sensor demonstrated excellent repeatability (RSD = 4.22%), reproducibility (RSD = 7.14%), and selectivity (RSD -1) with minimal matrix effects. This approach offers a simple, low-cost, and rapid method for paraoxon detection in water and food samples

An Intermittent Live Cell Imaging Screen for siRNA Enhancers and Suppressors of a Kinesin-5 Inhibitor

Kinesin-5 (also known as Eg5, KSP and Kif11) is required for assembly of a bipolar mitotic spindle. Small molecule inhibitors of Kinesin-5, developed as potential anti-cancer drugs, arrest cell in mitosis and promote apoptosis of cancer cells. We performed a genome-wide siRNA screen for enhancers and suppressors of a Kinesin-5 inhibitor in human cells to elucidate cellular responses, and thus identify factors that might predict drug sensitivity in cancers. Because the drug's actions play out over several days, we developed an intermittent imaging screen. Live HeLa cells expressing GFP-tagged histone H2B were imaged at 0, 24 and 48 hours after drug addition, and images were analyzed using open-source software that incorporates machine learning. This screen effectively identified siRNAs that caused increased mitotic arrest at low drug concentrations (enhancers), and vice versa (suppressors), and we report siRNAs that caused both effects. We then classified the effect of siRNAs for 15 genes where 3 or 4 out of 4 siRNA oligos tested were suppressors as assessed by time lapse imaging, and by testing for suppression of mitotic arrest in taxol and nocodazole. This identified 4 phenotypic classes of drug suppressors, which included known and novel genes. Our methodology should be applicable to other screens, and the suppressor and enhancer genes we identified may open new lines of research into mitosis and checkpoint biology

Integrating Flux Balance Analysis into Kinetic Models to Decipher the Dynamic Metabolism of Shewanella oneidensis MR-1

Shewanella oneidensis MR-1 sequentially utilizes lactate and its waste products (pyruvate and acetate) during batch culture. To decipher MR-1 metabolism, we integrated genome-scale flux balance analysis (FBA) into a multiple-substrate Monod model to perform the dynamic flux balance analysis (dFBA). The dFBA employed a static optimization approach (SOA) by dividing the batch time into small intervals (i.e., ∼400 mini-FBAs), then the Monod model provided time-dependent inflow/outflow fluxes to constrain the mini-FBAs to profile the pseudo-steady-state fluxes in each time interval. The mini-FBAs used a dual-objective function (a weighted combination of “maximizing growth rate” and “minimizing overall flux”) to capture trade-offs between optimal growth and minimal enzyme usage. By fitting the experimental data, a bi-level optimization of dFBA revealed that the optimal weight in the dual-objective function was time-dependent: the objective function was constant in the early growth stage, while the functional weight of minimal enzyme usage increased significantly when lactate became scarce. The dFBA profiled biologically meaningful dynamic MR-1 metabolisms: 1. the oxidative TCA cycle fluxes increased initially and then decreased in the late growth stage; 2. fluxes in the pentose phosphate pathway and gluconeogenesis were stable in the exponential growth period; and 3. the glyoxylate shunt was up-regulated when acetate became the main carbon source for MR-1 growth