Constraint-based Reconstruction Research Articles

Improved production of Taxol® precursors in S. cerevisiae using combinatorial in silico design and metabolic engineering

BackgroundIntegrated metabolic engineering approaches that combine system and synthetic biology tools enable the efficient design of microbial cell factories for synthesizing high-value products. In this study, we utilized in silico design algorithms on the yeast genome-scale model to predict genomic modifications that could enhance the production of early-step Taxol® in engineered Saccharomyces cerevisiae cells.ResultsUsing constraint-based reconstruction and analysis (COBRA) methods, we narrowed down the solution set of genomic modification candidates. We screened 17 genomic modifications, including nine gene deletions and eight gene overexpressions, through wet-lab studies to determine their impact on taxadiene production, the first metabolite in the Taxol® biosynthetic pathway. Under different cultivation conditions, most single genomic modifications resulted in increased taxadiene production. The strain named KM32, which contained four overexpressed genes (ILV2, TRR1, ADE13, and ECM31) involved in branched-chain amino acid biosynthesis, the thioredoxin system, de novo purine synthesis, and the pantothenate pathway, respectively, exhibited the best performance. KM32 achieved a 50% increase in taxadiene production, reaching 215 mg/L. Furthermore, KM32 produced the highest reported yields of taxa-4(20),11-dien-5α-ol (T5α-ol) at 43.65 mg/L and taxa-4(20),11-dien-5-α-yl acetate (T5αAc) at 26.2 mg/L among early-step Taxol® metabolites in S. cerevisiae.ConclusionsThis study highlights the effectiveness of computational and integrated approaches in identifying promising genomic modifications that can enhance the performance of yeast cell factories. By employing in silico design algorithms and wet-lab screening, we successfully improved taxadiene production in engineered S. cerevisiae strains. The best-performing strain, KM32, achieved substantial increases in taxadiene as well as production of T5α-ol and T5αAc. These findings emphasize the importance of using systematic and integrated strategies to develop efficient yeast cell factories, providing potential implications for the industrial production of high-value isoprenoids like Taxol®.

Mathematical reconstruction of the metabolic network in an in-vitro multiple myeloma model.

It is increasingly apparent that cancer cells, in addition to remodelling their metabolism to survive and proliferate, adapt and manipulate the metabolism of other cells. This property may be a telling sign that pre-clinical tumour metabolism studies exclusively utilising in-vitro mono-culture models could prove to be limited for uncovering novel metabolic targets able to translate into clinical therapies. Although this is increasingly recognised, and work towards addressing the issue is becoming routinary much remains poorly understood. For instance, knowledge regarding the biochemical mechanisms through which cancer cells manipulate non-cancerous cell metabolism, and the subsequent impact on their survival and proliferation remains limited. Additionally, the variations in these processes across different cancer types and progression stages, and their implications for therapy, also remain largely unexplored. This study employs an interdisciplinary approach that leverages the predictive power of mathematical modelling to enrich experimental findings. We develop a functional multicellular in-silico model that facilitates the qualitative and quantitative analysis of the metabolic network spawned by an in-vitro co-culture model of bone marrow mesenchymal stem- and myeloma cell lines. To procure this model, we devised a bespoke human genome constraint-based reconstruction workflow that combines aspects from the legacy mCADRE & Metabotools algorithms, the novel redHuman algorithm, along with 13C-metabolic flux analysis. Our workflow transforms the latest human metabolic network matrix (Recon3D) into two cell-specific models coupled with a metabolic network spanning a shared growth medium. When cross-validating our in-silico model against the in-vitro model, we found that the in-silico model successfully reproduces vital metabolic behaviours of its in-vitro counterpart; results include cell growth predictions, respiration rates, as well as support for observations which suggest cross-shuttling of redox-active metabolites between cells.

Cardinality optimization in constraint-based modelling: application to human metabolism.

Several applications in constraint-based modelling can be mathematically formulated as cardinality optimization problems involving the minimization or maximization of the number of nonzeros in a vector. These problems include testing for stoichiometric consistency, testing for flux consistency, testing for thermodynamic flux consistency, computing sparse solutions to flux balance analysis problems and computing the minimum number of constraints to relax to render an infeasible flux balance analysis problem feasible. Such cardinality optimization problems are computationally complex, with no known polynomial time algorithms capable of returning an exact and globally optimal solution. By approximating the zero-norm with nonconvex continuous functions, we reformulate a set of cardinality optimization problems in constraint-based modelling into a difference of convex functions. We implemented and numerically tested novel algorithms that approximately solve the reformulated problems using a sequence of convex programs. We applied these algorithms to various biochemical networks and demonstrate that our algorithms match or outperform existing related approaches. In particular, we illustrate the efficiency and practical utility of our algorithms for cardinality optimization problems that arise when extracting a model ready for thermodynamic flux balance analysis given a human metabolic reconstruction. Open source scripts to reproduce the results are here https://github.com/opencobra/COBRA.papers/2023_cardOpt with general purpose functions integrated within the COnstraint-Based Reconstruction and Analysis toolbox: https://github.com/opencobra/cobratoolbox.

A structured evaluation of genome-scale constraint-based modeling tools for microbial consortia.

Harnessing the power of microbial consortia is integral to a diverse range of sectors, from healthcare to biotechnology to environmental remediation. To fully realize this potential, it is critical to understand the mechanisms behind the interactions that structure microbial consortia and determine their functions. Constraint-based reconstruction and analysis (COBRA) approaches, employing genome-scale metabolic models (GEMs), have emerged as the state-of-the-art tool to simulate the behavior of microbial communities from their constituent genomes. In the last decade, many tools have been developed that use COBRA approaches to simulate multi-species consortia, under either steady-state, dynamic, or spatiotemporally varying scenarios. Yet, these tools have not been systematically evaluated regarding their software quality, most suitable application, and predictive power. Hence, it is uncertain which tools users should apply to their system and what are the most urgent directions that developers should take in the future to improve existing capacities. This study conducted a systematic evaluation of COBRA-based tools for microbial communities using datasets from two-member communities as test cases. First, we performed a qualitative assessment in which we evaluated 24 published tools based on a list of FAIR (Findability, Accessibility, Interoperability, and Reusability) features essential for software quality. Next, we quantitatively tested the predictions in a subset of 14 of these tools against experimental data from three different case studies: a) syngas fermentation by C. autoethanogenum and C. kluyveri for the static tools, b) glucose/xylose fermentation with engineered E. coli and S. cerevisiae for the dynamic tools, and c) a Petri dish of E. coli and S. enterica for tools incorporating spatiotemporal variation. Our results show varying performance levels of the best qualitatively assessed tools when examining the different categories of tools. The differences in the mathematical formulation of the approaches and their relation to the results were also discussed. Ultimately, we provide recommendations for refining future GEM microbial modeling tools.

Practical sampling of constraint-based models: Optimized thinning boosts CHRR performance.

Thinning is a sub-sampling technique to reduce the memory footprint of Markov chain Monte Carlo. Despite being commonly used, thinning is rarely considered efficient. For sampling constraint-based models, a highly relevant use-case in systems biology, we here demonstrate that thinning boosts computational and, thereby, sampling efficiencies of the widely used Coordinate Hit-and-Run with Rounding (CHRR) algorithm. By benchmarking CHRR with thinning with simplices and genome-scale metabolic networks of up to thousands of dimensions, we find a substantial increase in computational efficiency compared to unthinned CHRR, in our examples by orders of magnitude, as measured by the effective sample size per time (ESS/t), with performance gains growing with polytope (effective network) dimension. Using a set of benchmark models we derive a ready-to-apply guideline for tuning thinning to efficient and effective use of compute resources without requiring additional coding effort. Our guideline is validated using three (out-of-sample) large-scale networks and we show that it allows sampling convex polytopes uniformly to convergence in a fraction of time, thereby unlocking the rigorous investigation of hitherto intractable models. The derivation of our guideline is explained in detail, allowing future researchers to update it as needed as new model classes and more training data becomes available. CHRR with deliberate utilization of thinning thereby paves the way to keep pace with progressing model sizes derived with the constraint-based reconstruction and analysis (COBRA) tool set. Sampling and evaluation pipelines are available at https://jugit.fz-juelich.de/IBG-1/ModSim/fluxomics/chrrt.

A High-Quality Genome-Scale Model for Rhodococcus opacus Metabolism.

Rhodococcus opacus is a bacterium that has a high tolerance to aromatic compounds and can produce significant amounts of triacylglycerol (TAG). Here, we present iGR1773, the first genome-scale model (GSM) of R. opacus PD630 metabolism based on its genomic sequence and associated data. The model includes 1773 genes, 3025 reactions, and 1956 metabolites, was developed in a reproducible manner using CarveMe, and was evaluated through Metabolic Model tests (MEMOTE). We combine the model with two Constraint-Based Reconstruction and Analysis (COBRA) methods that use transcriptomics data to predict growth rates and fluxes: E-Flux2 and SPOT (Simplified Pearson Correlation with Transcriptomic data). Growth rates are best predicted by E-Flux2. Flux profiles are more accurately predicted by E-Flux2 than flux balance analysis (FBA) and parsimonious FBA (pFBA), when compared to 44 central carbon fluxes measured by 13C-Metabolic Flux Analysis (13C-MFA). Under glucose-fed conditions, E-Flux2 presents an R2 value of 0.54, while predictions based on pFBA had an inferior R2 of 0.28. We attribute this improved performance to the extra activity information provided by the transcriptomics data. For phenol-fed metabolism, in which the substrate first enters the TCA cycle, E-Flux2's flux predictions display a high R2 of 0.96 while pFBA showed an R2 of 0.93. We also show that glucose metabolism and phenol metabolism function with similar relative ATP maintenance costs. These findings demonstrate that iGR1773 can help the metabolic engineering community predict aromatic substrate utilization patterns and perform computational strain design.

Constraint-Based Reconstruction and Analyses of Metabolic Models: Open-Source Python Tools and Applications to Cancer.

The influence of metabolism on signaling, epigenetic markers, and transcription is highly complex yet important for understanding cancer physiology. Despite the development of high-resolution multi-omics technologies, it is difficult to infer metabolic activity from these indirect measurements. Fortunately, genome-scale metabolic models and constraint-based modeling provide a systems biology framework to investigate the metabolic states and define the genotype-phenotype associations by integrations of multi-omics data. Constraint-Based Reconstruction and Analysis (COBRA) methods are used to build and simulate metabolic networks using mathematical representations of biochemical reactions, gene-protein reaction associations, and physiological and biochemical constraints. These methods have led to advancements in metabolic reconstruction, network analysis, perturbation studies as well as prediction of metabolic state. Most computational tools for performing these analyses are written for MATLAB, a proprietary software. In order to increase accessibility and handle more complex datasets and models, community efforts have started to develop similar open-source tools in Python. To date there is a comprehensive set of tools in Python to perform various flux analyses and visualizations; however, there are still missing algorithms in some key areas. This review summarizes the availability of Python software for several components of COBRA methods and their applications in cancer metabolism. These tools are evolving rapidly and should offer a readily accessible, versatile way to model the intricacies of cancer metabolism for identifying cancer-specific metabolic features that constitute potential drug targets.

NCMW: A Python Package to Analyze Metabolic Interactions in the Nasal Microbiome.

The human upper respiratory tract is the reservoir of a diverse community of commensals and potential pathogens (pathobionts), including Streptococcus pneumoniae (pneumococcus), Haemophilus influenzae, Moraxella catarrhalis, and Staphylococcus aureus, which occasionally turn into pathogens causing infectious diseases, while the contribution of many nasal microorganisms to human health remains undiscovered. To better understand the composition of the nasal microbiome community, we create a workflow of the community model, which mimics the human nasal environment. To address this challenge, constraint-based reconstruction of biochemically accurate genome-scale metabolic models (GEMs) networks of microorganisms is mandatory. Our workflow applies constraint-based modeling (CBM), simulates the metabolism between species in a given microbiome, and facilitates generating novel hypotheses on microbial interactions. Utilizing this workflow, we hope to gain a better understanding of interactions from the metabolic modeling perspective. This article presents nasal community modeling workflow (NCMW)—a python package based on GEMs of species as a starting point for understanding the composition of the nasal microbiome community. The package is constructed as a step-by-step mathematical framework for metabolic modeling and analysis of the nasal microbial community. Using constraint-based models reduces the need for culturing species in vitro, a process that is not convenient in the environment of human noses. Availability: NCMW is freely available on the Python Package Index (PIP) via pip install NCMW. The source code, documentation, and usage examples (Jupyter Notebook and example files) are available at https://github.com/manuelgloeckler/ncmw.

Development of an integrating systems metabolic engineering and bioprocess modeling approach for rational strain improvement

Due to metabolic engineering advances, it is possible to produce a large number of bio-based chemical compounds. However, the downstream processing is challenging the development at an industrial scale. Selecting the proper recovery and separation method is critical in developing a cost-effective and competitive process. This study is focused on analyzing the production of succinic acid from glycerol through computational methods for strain design and determine the effects on the recovery and purification processes. Here, it was integrated the metabolic phenotype of different mutants obtained through constraint-based reconstruction and analysis modeling into two recoveries and purification processes by using a flash evaporator and a reactive extraction column. Two thousand forty-two mutants, predicted in silico, were phenotypically analyzed. Succinate was the main product; however, were also generated by-products such as acetate, formate, and lactate. Of those, twelve mutants were selected to create feed streams with different compositions and then downstream modeling. These models, after optimization, achieved a recovery rate > 93% and molar purity > 99%. Economic analysis suggests that when engineering a strain, the most important factors to consider are (i) the productivity (> 5 mmol/gDW h−1 of succinic acid) and (ii) the growth rate (> 0.10 gDW h−1).

COBREXA.jl: constraint-based reconstruction and exascale analysis.

SummaryCOBREXA.jl is a Julia package for scalable, high-performance constraint-based reconstruction and analysis of very large-scale biological models. Its primary purpose is to facilitate the integration of modern high performance computing environments with the processing and analysis of large-scale metabolic models of challenging complexity. We report the architecture of the package, and demonstrate how the design promotes analysis scalability on several use-cases with multi-organism community models.Availability and implementation https://doi.org/10.17881/ZKCR-BT30.Supplementary information Supplementary data are available at Bioinformatics online.

Curating COBRA Models of Microbial Metabolism.

Constraint-based reconstruction and analysis (COBRA) methods have been used for over 20years to generate genome-scale models of metabolism in biological systems. The COBRA models have been utilized to gain new insights into the biochemical conversions that occur within organisms and allow their survival and proliferation. Using these models, computational biologists can conduct a variety of different analyses such as examining network structures, predicting metabolic capabilities, resolving unexplained experimental observations, generating and testing new hypotheses, assessing the nutritional requirements of a biosystem and approximating its environmental niche, identifying missing enzymatic functions in the annotated genomes, and engineering desired metabolic capabilities in model organisms. This chapter details the protocol for developing curated system-level COBRA models of metabolism in microbes.

Topologically Consistent Reconstruction for Complex Indoor Structures from Point Clouds

Indoor structures are composed of ceilings, walls and floors that need to be modeled for a variety of applications. This paper proposes an approach to reconstructing models of indoor structures in complex environments. First, semantic pre-processing, including segmentation and occlusion construction, is applied to segment the input point clouds to generate semantic patches of structural primitives with uniform density. Then, a primitives extraction method with detected boundary is introduced to approximate both the mathematical surface and the boundary of the patches. Finally, a constraint-based model reconstruction is applied to achieve the final topologically consistent structural model. Under this framework, both the geometric and structural constraints are considered in a holistic manner to assure topologic regularity. Experiments were carried out with both synthetic and real-world datasets. The accuracy of the proposed method achieved an overall reconstruction quality of approximately 4.60 cm of root mean square error (RMSE) and 94.10% Intersection over Union (IoU) of the input point cloud. The development can be applied for structural reconstruction of various complex indoor environments.

DEMETER: efficient simultaneous curation of genome-scale reconstructions guided by experimental data and refined gene annotations.

MotivationManual curation of genome-scale reconstructions is laborious, yet existing automated curation tools do not typically take species-specific experimental and curated genomic data into account.ResultsWe developed Data-drivEn METabolic nEtwork Refinement (DEMETER), a Constraint-Based Reconstruction and Analysis (COBRA) Toolbox extension, which enables the efficient, simultaneous refinement of thousands of draft genome-scale reconstructions, while ensuring adherence to the quality standards in the field, agreement with available experimental data and refinement of pathways based on manually refined genome annotations.Availability and implementationDEMETER and tutorials are freely available at https://github.com/opencobra.Supplementary information Supplementary data are available at Bioinformatics online.

Genome-Scale Metabolic Modeling of the Human Microbiome in the Era of Personalized Medicine

The human microbiome plays an important role in human health and disease. Meta-omics analyses provide indispensable data for linking changes in microbiome composition and function to disease etiology. Yet, the lack of a mechanistic understanding of, e.g., microbiome-metabolome links hampers the translation of these findings into effective, novel therapeutics. Here, we propose metabolic modeling of microbial communities through constraint-based reconstruction and analysis (COBRA) as a complementary approach to meta-omics analyses. First, we highlight the importance of microbial metabolism in cardiometabolic diseases, inflammatory bowel disease, colorectal cancer, Alzheimer disease, and Parkinson disease. Next, we demonstrate that microbial community modeling can stratify patients and controls, mechanistically link microbes with fecal metabolites altered in disease, and identify host pathways affected by the microbiome. Finally, we outline our vision for COBRA modeling combined with meta-omics analyses and multivariate statistical analyses to inform and guide clinical trials, yield testable hypotheses, and ultimately propose novel dietary and therapeutic interventions.

MASSpy: Building, simulating, and visualizing dynamic biological models in Python using mass action kinetics.

Mathematical models of metabolic networks utilize simulation to study system-level mechanisms and functions. Various approaches have been used to model the steady state behavior of metabolic networks using genome-scale reconstructions, but formulating dynamic models from such reconstructions continues to be a key challenge. Here, we present the Mass Action Stoichiometric Simulation Python (MASSpy) package, an open-source computational framework for dynamic modeling of metabolism. MASSpy utilizes mass action kinetics and detailed chemical mechanisms to build dynamic models of complex biological processes. MASSpy adds dynamic modeling tools to the COnstraint-Based Reconstruction and Analysis Python (COBRApy) package to provide an unified framework for constraint-based and kinetic modeling of metabolic networks. MASSpy supports high-performance dynamic simulation through its implementation of libRoadRunner: the Systems Biology Markup Language (SBML) simulation engine. Three examples are provided to demonstrate how to use MASSpy: (1) a validation of the MASSpy modeling tool through dynamic simulation of detailed mechanisms of enzyme regulation; (2) a feature demonstration using a workflow for generating ensemble of kinetic models using Monte Carlo sampling to approximate missing numerical values of parameters and to quantify biological uncertainty, and (3) a case study in which MASSpy is utilized to overcome issues that arise when integrating experimental data with the computation of functional states of detailed biological mechanisms. MASSpy represents a powerful tool to address challenges that arise in dynamic modeling of metabolic networks, both at small and large scales.

Reconstruction of a Context-Specific Model Based on Genome-Scale Metabolic Simulation for Identification of Prochloraz Resistance Mechanisms in Penicillium digitatum.

Penicillium digitatum is the most destructive postharvest pathogen of citrus fruits, causing substantial economic losses. Prochloraz-resistant strains have emerged due to overuse of imidazole fungicides in agriculture. To study the prochloraz resistance mechanisms at the system level, a genome-scale metabolic model (GEM, iPD1512) of P. digitatum was reconstructed and constrained based on context-specific transcriptome data of the prochloraz-resistant strain, PdF6, from our previous work, a newly sequenced, context-specific transcriptome result of the major facilitator superfamily transporter-encoding gene mfs2 knockout mutant PdF6Δmfs2, and experimentally derived growth rate data. Through the model, iPD1512, the processes of prochloraz resistance in P. digitatum were well simulated. In detail, the growth rates of both wild-type and mutant P. digitatum under different prochloraz concentrations were simulated using constraint-based reconstruction and analysis. The growth rates of the mutant strains (sterol regulatory element-binding protein-encoding gene sreA knockout mutant PdF6ΔsreA and PdF6Δmfs2) were calculated and confirmed to be consistent with the simulation results. Furthermore, correlations between genes and prochloraz resistance were predicted and showed a great difference when compared with correlation results based on p-values from the hypothesis testing used by comparative transcriptomics. To sum up, in contrast to traditional transcriptome analysis, the GEM provides a systemic and dynamic drug resistance mechanism, which might help to detect some key upstream regulatory genes, but with small expression changes, and might provide more efficient targets to control prochloraz-resistant P. digitatum.

Medusa: Software to build and analyze ensembles of genome-scale metabolic network reconstructions.

Uncertainty in the structure and parameters of networks is ubiquitous across computational biology. In constraint-based reconstruction and analysis of metabolic networks, this uncertainty is present both during the reconstruction of networks and in simulations performed with them. Here, we present Medusa, a Python package for the generation and analysis of ensembles of genome-scale metabolic network reconstructions. Medusa builds on the COBRApy package for constraint-based reconstruction and analysis by compressing a set of models into a compact ensemble object, providing functions for the generation of ensembles using experimental data, and extending constraint-based analyses to ensemble scale. We demonstrate how Medusa can be used to generate ensembles and perform ensemble simulations, and how machine learning can be used in conjunction with Medusa to guide the curation of genome-scale metabolic network reconstructions. Medusa is available under the permissive MIT license from the Python Packaging Index (https://pypi.org) and from github (https://github.com/opencobra/Medusa), and comprehensive documentation is available at https://medusa.readthedocs.io/en/latest.

StructureDistiller: Structural relevance scoring identifies the most informative entries of a contact map

Protein folding and structure prediction are two sides of the same coin. Contact maps and the related techniques of constraint-based structure reconstruction can be considered as unifying aspects of both processes. We present the Structural Relevance (SR) score which quantifies the information content of individual contacts and residues in the context of the whole native structure. The physical process of protein folding is commonly characterized with spatial and temporal resolution: some residues are Early Folding while others are Highly Stable with respect to unfolding events. We employ the proposed SR score to demonstrate that folding initiation and structure stabilization are subprocesses realized by distinct sets of residues. The example of cytochrome c is used to demonstrate how StructureDistiller identifies the most important contacts needed for correct protein folding. This shows that entries of a contact map are not equally relevant for structural integrity. The proposed StructureDistiller algorithm identifies contacts with the highest information content; these entries convey unique constraints not captured by other contacts. Identification of the most informative contacts effectively doubles resilience toward contacts which are not observed in the native contact map. Furthermore, this knowledge increases reconstruction fidelity on sparse contact maps significantly by 0.4 Å.

Creation and analysis of biochemical constraint-based models using the COBRA Toolbox v.3.0.

Constraint-based reconstruction and analysis (COBRA) provides a molecular mechanistic framework for integrative analysis of experimental molecular systems biology data and quantitative prediction of physicochemically and biochemically feasible phenotypic states. The COBRA Toolbox is a comprehensive desktop software suite of interoperable COBRA methods. It has found widespread application in biology, biomedicine, and biotechnology because its functions can be flexibly combined to implement tailored COBRA protocols for any biochemical network. This protocol is an update to the COBRA Toolbox v.1.0 and v.2.0. Version 3.0 includes new methods for quality-controlled reconstruction, modeling, topological analysis, strain and experimental design, and network visualization, as well as network integration of chemoinformatic, metabolomic, transcriptomic, proteomic, and thermochemical data. New multi-lingual code integration also enables an expansion in COBRA application scope via high-precision, high-performance, and nonlinear numerical optimization solvers for multi-scale, multi-cellular, and reaction kinetic modeling, respectively. This protocol provides an overview of all these new features and can be adapted to generate and analyze constraint-based models in a wide variety of scenarios. The COBRA Toolbox v.3.0 provides an unparalleled depth of COBRA methods.

A simplified metabolic network reconstruction to promote understanding and development of flux balance analysis tools

GEnome-scale Network REconstructions (GENREs) mathematically describe metabolic reactions of an organism or a specific cell type. GENREs can be used with a number of constraint-based reconstruction and analysis (COBRA) methods to make computational predictions on how a system changes in different environments. We created a simplified GENRE (referred to as iSIM) that captures central energy metabolism with nine metabolic reactions to illustrate the use of and promote the understanding of GENREs and constraint-based methods. We demonstrate the simulation of single and double gene deletions, flux variability analysis (FVA), and test a number of metabolic tasks with the GENRE. Code to perform these analyses is provided in Python, R, and MATLAB. Finally, with iSIM as a guide, we demonstrate how inaccuracies in GENREs can limit their use in the interrogation of energy metabolism.