GeneReg: a constraint-based approach for design of feasible metabolic engineering strategies at the gene level.
Ontology highlight
ABSTRACT: Motivation
Large-scale metabolic models are widely used to design metabolic engineering strategies for diverse biotechnological applications. However, the existing computational approaches focus on alteration of reaction fluxes and often neglect the manipulations of gene expression to implement these strategies.Results
Here, we find that the association of genes with multiple reactions leads to infeasibility of engineering strategies at the flux level, since they require contradicting manipulations of gene expression. Moreover, we identify that all of the existing approaches to design gene knockout strategies do not ensure that the resulting design may also require other gene alterations, such as up- or downregulations, to match the desired flux distribution. To address these issues, we propose a constraint-based approach, termed GeneReg, that facilitates the design of feasible metabolic engineering strategies at the gene level and that is readily applicable to large-scale metabolic networks. We show that GeneReg can identify feasible strategies to overproduce ethanol in Escherichia coli and lactate in Saccharomyces cerevisiae, but overproduction of the TCA cycle intermediates is not feasible in five organisms used as cell factories under default growth conditions. Therefore, GeneReg points at the need to couple gene regulation and metabolism to design rational metabolic engineering strategies.Availability and implementation
https://github.com/MonaRazaghi/GeneReg.Supplementary information
Supplementary data are available at Bioinformatics online.
SUBMITTER: Razaghi-Moghadam Z
PROVIDER: S-EPMC8289378 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:While numerous computational methods have been developed that use genome-scale models to propose mutants for the purpose of metabolic engineering, they generally compare mutants based on a single criteria (e.g., production rate at a mutant's maximum growth rate). As such, these approaches remain limited in their ability to include multiple complex engineering constraints. To address this shortcoming, we have developed feasible space and shadow price constraint (FaceCon and ShadowCon) modules that can be added to existing mixed integer linear adaptive evolution metabolic engineering algorithms, such as OptKnock and OptORF. These modules allow strain designs to be identified amongst a set of multiple metabolic engineering algorithm solutions that are capable of high chemical production while also satisfying additional design criteria. We describe the various module implementations and their potential applications to the field of metabolic engineering. We then incorporated these modules into the OptORF metabolic engineering algorithm. Using an Escherichia coli genome-scale model (iJO1366), we generated different strain designs for the anaerobic production of ethanol from glucose, thus demonstrating the tractability and potential utility of these modules in metabolic engineering algorithms.
| S-EPMC4249821 | biostudies-literature
Project description:BackgroundConstraint-based computational approaches, such as flux balance analysis (FBA), have proven successful in modeling genome-level metabolic behavior for conditions where a set of simple cellular objectives can be clearly articulated. Recently, the necessity to expand the current range of constraint-based methods to incorporate high-throughput experimental data has been acknowledged by the proposal of several methods. However, these methods have rarely been used to address cellular metabolic responses to some relevant perturbations such as antimicrobial or temperature-induced stress. Here, we present a new method for combining gene-expression data with FBA (GX-FBA) that allows modeling of genome-level metabolic response to a broad range of environmental perturbations within a constraint-based framework. The method uses mRNA expression data to guide hierarchical regulation of cellular metabolism subject to the interconnectivity of the metabolic network.ResultsWe applied GX-FBA to a genome-scale model of metabolism in the gram negative bacterium Yersinia pestis and analyzed its metabolic response to (i) variations in temperature known to induce virulence, and (ii) antibiotic stress. Without imposition of any a priori behavioral constraints, our results show strong agreement with reported phenotypes. Our analyses also lead to novel insights into how Y. pestis uses metabolic adjustments to counter different forms of stress.ConclusionsComparisons of GX-FBA predicted metabolic states with fluxomic measurements and different reported post-stress phenotypes suggest that mass conservation constraints and network connectivity can be an effective representative of metabolic flux regulation in constraint-based models. We believe that our approach will be of aid in the in silico evaluation of cellular goals under different conditions and can be used for a variety of analyses such as identification of potential drug targets and their action.
| S-EPMC3572438 | biostudies-literature
Project description:MOTIVATION:Flux balance analysis (FBA) based bilevel optimization has been a great success in redesigning metabolic networks for biochemical overproduction. To date, many computational approaches have been developed to solve the resulting bilevel optimization problems. However, most of them are of limited use due to biased optimality principle, poor scalability with the size of metabolic networks, potential numeric issues or low quantity of design solutions in a single run. RESULTS:Here, we have employed a network interdiction model free of growth optimality assumptions, a special case of bilevel optimization, for computational strain design and have developed a hybrid Benders algorithm (HBA) that deals with complicating binary variables in the model, thereby achieving high efficiency without numeric issues in search of best design strategies. More importantly, HBA can list solutions that meet users' production requirements during the search, making it possible to obtain numerous design strategies at a small runtime overhead (typically ?1?h, e.g. studied in this article). AVAILABILITY AND IMPLEMENTATION:Source code implemented in the MATALAB Cobratoolbox is freely available at https://github.com/chang88ye/NIHBA. CONTACT:math4neu@gmail.com or natalio.krasnogor@ncl.ac.uk. SUPPLEMENTARY INFORMATION:Supplementary data are available at Bioinformatics online.
| S-EPMC7267835 | biostudies-literature
Project description:Background:Pseudomonas putida is a promising candidate for the industrial production of biofuels and biochemicals because of its high tolerance to toxic compounds and its ability to grow on a wide variety of substrates. Engineering this organism for improved performances and predicting metabolic responses upon genetic perturbations requires reliable descriptions of its metabolism in the form of stoichiometric and kinetic models. Results:In this work, we developed kinetic models of P. putida to predict the metabolic phenotypes and design metabolic engineering interventions for the production of biochemicals. The developed kinetic models contain 775 reactions and 245 metabolites. Furthermore, we introduce here a novel set of constraints within thermodynamics-based flux analysis that allow for considering concentrations of metabolites that exist in several compartments as separate entities. We started by a gap-filling and thermodynamic curation of iJN1411, the genome-scale model of P. putida KT2440. We then systematically reduced the curated iJN1411 model, and we created three core stoichiometric models of different complexity that describe the central carbon metabolism of P. putida. Using the medium complexity core model as a scaffold, we generated populations of large-scale kinetic models for two studies. In the first study, the developed kinetic models successfully captured the experimentally observed metabolic responses to several single-gene knockouts of a wild-type strain of P. putida KT2440 growing on glucose. In the second study, we used the developed models to propose metabolic engineering interventions for improved robustness of this organism to the stress condition of increased ATP demand. Conclusions:The study demonstrates the potential and predictive capabilities of the kinetic models that allow for rational design and optimization of recombinant P. putida strains for improved production of biofuels and biochemicals. The curated genome-scale model of P. putida together with the developed large-scale stoichiometric and kinetic models represents a significant resource for researchers in industry and academia.
| S-EPMC7048048 | biostudies-literature
Project description:BackgroundHydrogenobyrinic acid is a key intermediate of the de-novo aerobic biosynthesis pathway of vitamin B12. The introduction of a heterologous de novo vitamin B12 biosynthesis pathway in Escherichia coli offers an alternative approach for its production. Although E. coli avoids major limitations that currently faced by industrial producers of vitamin B12, such as long growth cycles, the insufficient supply of hydrogenobyrinic acid restricts industrial vitamin B12 production.ResultsBy designing combinatorial ribosomal binding site libraries of the hemABCD genes in vivo, we found that their optimal relative translational initiation rates are 10:1:1:5. The transcriptional coordination of the uroporphyrinogen III biosynthetic module was realized by promoter engineering of the hemABCD operon. Knockdown of competitive heme and siroheme biosynthesis pathways by RBS engineering enhanced the hydrogenobyrinic acid titer to 20.54 and 15.85 mg L-1, respectively. Combined fine-tuning of the heme and siroheme biosynthetic pathways enhanced the hydrogenobyrinic acid titer to 22.57 mg L-1, representing a remarkable increase of 1356.13% compared with the original strain FH215-HBA.ConclusionsThrough multi-level metabolic engineering strategies, we achieved the metabolic balance of the uroporphyrinogen III biosynthesis pathway, eliminated toxicity due to by-product accumulation, and finally achieved a high HBA titer of 22.57 mg L-1 in E. coli. This lays the foundation for high-yield production of vitamin B12 in E. coli and will hopefully accelerate its industrial production.
| S-EPMC7268678 | biostudies-literature
Project description:In the growing field of metabolic engineering, where cells are treated as 'factories' that synthesize industrial compounds, it is essential to consider the ability of the cells' native metabolism to accommodate the demands of synthetic pathways, as these pathways will alter the homeostasis of cellular energy and electron metabolism. From the breakdown of substrate, microorganisms activate and reduce key co-factors such as ATP and NAD(P)H, which subsequently need to be hydrolysed and oxidized, respectively, in order to restore cellular balance. A balanced supply and consumption of such co-factors, here termed co-factor balance, will influence biotechnological performance. To aid the strain selection and design process, we used stoichiometric modelling (FBA, pFBA, FVA and MOMA) and the Escherichia coli (E.coli) core stoichiometric model to investigate the network-wide effect of butanol and butanol precursor production pathways differing in energy and electron demand on product yield. An FBA-based co-factor balance assessment (CBA) algorithm was developed to track and categorise how ATP and NAD(P)H pools are affected in the presence of a new pathway. CBA was compared to the balance calculations proposed by Dugar et al. (Nature Biotechnol. 29 (12), 1074-1078). Predicted solutions were compromised by excessively underdetermined systems, displaying greater flexibility in the range of reaction fluxes than experimentally measured by 13C-metabolic flux analysis (MFA) and the appearance of unrealistic futile co-factor cycles. With the assumption that futile cycles are tightly regulated in reality, the FBA models were manually constrained in a step-wise manner. Solutions with minimal futile cycling diverted surplus energy and electrons towards biomass formation. As an alternative, the use of loopless FBA or constraining the models with measured flux ranges were tried but did not prevent futile co-factor cycles. The results highlight the need to account for co-factor imbalance and confirm that better-balanced pathways with minimal diversion of surplus towards biomass formation present the highest theoretical yield. The analysis also suggests that ATP and NAD(P)H balancing cannot be assessed in isolation from each other, or even from the balance of additional co-factors such as AMP and ADP. We conclude that, through revealing the source of co-factor imbalance CBA can facilitate pathway and host selection when designing new biocatalysts for implementation by metabolic engineering.
| S-EPMC7440669 | biostudies-literature
Project description:BackgroundThe optimization of metabolic rates (as linear objective functions) represents the methodical core of flux-balance analysis techniques which have become a standard tool for the study of genome-scale metabolic models. Besides (growth and synthesis) rates, metabolic yields are key parameters for the characterization of biochemical transformation processes, especially in the context of biotechnological applications. However, yields are ratios of rates, and hence the optimization of yields (as nonlinear objective functions) under arbitrary linear constraints is not possible with current flux-balance analysis techniques. Despite the fundamental importance of yields in constraint-based modeling, a comprehensive mathematical framework for yield optimization is still missing.ResultsWe present a mathematical theory that allows one to systematically compute and analyze yield-optimal solutions of metabolic models under arbitrary linear constraints. In particular, we formulate yield optimization as a linear-fractional program. For practical computations, we transform the linear-fractional yield optimization problem to a (higher-dimensional) linear problem. Its solutions determine the solutions of the original problem and can be used to predict yield-optimal flux distributions in genome-scale metabolic models. For the theoretical analysis, we consider the linear-fractional problem directly. Most importantly, we show that the yield-optimal solution set (like the rate-optimal solution set) is determined by (yield-optimal) elementary flux vectors of the underlying metabolic model. However, yield- and rate-optimal solutions may differ from each other, and hence optimal (biomass or product) yields are not necessarily obtained at solutions with optimal (growth or synthesis) rates. Moreover, we discuss phase planes/production envelopes and yield spaces, in particular, we prove that yield spaces are convex and provide algorithms for their computation. We illustrate our findings by a small example and demonstrate their relevance for metabolic engineering with realistic models of E. coli.ConclusionsWe develop a comprehensive mathematical framework for yield optimization in metabolic models. Our theory is particularly useful for the study and rational modification of cell factories designed under given yield and/or rate requirements.
| S-EPMC5992331 | biostudies-literature
Project description:Antibody immunodominance refers to the preferential and asymmetric elicitation of antibodies against specific epitopes on a complex protein antigen. Traditional vaccination approaches for rapidly evolving pathogens have had limited success in part because of this phenomenon, as elicited antibodies preferentially target highly variable regions of antigens, and thus do not confer long lasting protection. While antibodies targeting functionally conserved epitopes have the potential to be broadly protective, they often make up a minority of the overall repertoire. Here, we discuss recent protein engineering strategies used to favorably alter patterns of immunodominance, and selectively focus antibody responses toward broadly protective epitopes in the pursuit of next-generation vaccines for rapidly evolving pathogens.
| S-EPMC8683408 | biostudies-literature
Project description:Although bacteria with multipartite genomes are prevalent, our knowledge of the mechanisms maintaining their genome is very limited, and much remains to be learned about the structural and functional interrelationships of multiple chromosomes. Owing to its bi-chromosomal genome architecture and its importance in public health, Vibrio cholerae, the causative agent of cholera, has become a preferred model to study bacteria with multipartite genomes. However, most in vivo studies in V. cholerae have been hampered by its genome architecture, as it is difficult to give phenotypes to a specific chromosome. This difficulty was surmounted using a unique and powerful strategy based on massive rearrangement of prokaryotic genomes. We developed a site-specific recombination-based engineering tool, which allows targeted, oriented, and reciprocal DNA exchanges. Using this genetic tool, we obtained a panel of V. cholerae mutants with various genome configurations: one with a single chromosome, one with two chromosomes of equal size, and one with both chromosomes controlled by identical origins. We used these synthetic strains to address several biological questions--the specific case of the essentiality of Dam methylation in V. cholerae and the general question concerning bacteria carrying circular chromosomes--by looking at the effect of chromosome size on topological issues. In this article, we show that Dam, RctB, and ParA2/ParB2 are strictly essential for chrII origin maintenance, and we formally demonstrate that the formation of chromosome dimers increases exponentially with chromosome size.
| S-EPMC3257285 | biostudies-literature
Project description:MotivationThe computer-aided design of metabolic intervention strategies has become a key component of an integrated metabolic engineering approach and a broad range of methods and algorithms has been developed for this task. Many of these algorithms enforce coupling of growth with product synthesis and may return thousands of possible intervention strategies from which the most suitable strategy must then be selected.ResultsThis work focuses on how to evaluate and rank, in a meaningful way, a given pool of computed metabolic engineering strategies for growth-coupled product synthesis. Apart from straightforward criteria, such as a preferably small number of necessary interventions, a reasonable growth rate and a high product yield, we present several new criteria useful to pick the most suitable intervention strategy. Among others, we investigate the robustness of the intervention strategies by searching for metabolites that may disrupt growth coupling when accumulated or secreted and by checking whether the interventions interrupt pathways at their origin (preferable) or at downstream steps. We also assess thermodynamic properties of the pathway(s) favored by the intervention strategy. Furthermore, strategies that have a significant overlap with alternative solutions are ranked higher because they provide flexibility in implementation. We also introduce the notion of equivalence classes for grouping intervention strategies with identical solution spaces. Our ranking procedure involves in total ten criteria and we demonstrate its applicability by assessing knockout-based intervention strategies computed in a genome-scale model of E.coli for the growth-coupled synthesis of l-methionine and of the heterologous product 1,4-butanediol.Availability and implementationThe MATLAB scripts that were used to characterize and rank the example intervention strategies are available at http://www2.mpi-magdeburg.mpg.de/projects/cna/etcdownloads.html.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC6735923 | biostudies-literature