Model-guided dynamic control of essential metabolic nodes boosts acetyl-coenzyme A–dependent bioproduction in rewired Pseudomonas putida

Pseudomonas putida is evolutionarily endowed with features relevant for bioproduction, especially under harsh operating conditions. The rich metabolic versatility of this species, however, comes at the price of limited formation of acetyl-coenzyme A (CoA) from sugar substrates. Since acetyl-CoA is a key metabolic precursor for a number of added-value products, in this work we deployed an in silico-guided rewiring program of central carbon metabolism for upgrading P. putida as a host for acetyl-CoA–dependent bioproduction. An updated kinetic model, integrating fluxomics and metabolomics datasets in addition to manually-curated information of enzyme mechanisms, identified targets that would lead to increased acetyl-CoA levels. Based on these predictions, a set of plasmids based on clustered regularly interspaced short palindromic repeats (CRISPR) and dead CRISPR-associated protein 9 (dCas9) was constructed to silence genes by CRISPR interference (CRISPRi). Dynamic reduction of gene expression of two key targets (gltA, encoding citrate synthase, and the essential accA gene, encoding subunit A of the acetyl-CoA carboxylase complex) mediated an 8-fold increase in the acetyl-CoA content of rewired P. putida. Poly(3-hydroxybutyrate) (PHB) was adopted as a proxy of acetyl-CoA availability, and two synthetic pathways were engineered for biopolymer accumulation. By including cell morphology as an extra target for the CRISPRi approach, fully rewired P. putida strains programmed for PHB accumulation had a 5-fold increase in PHB titers in bioreactor cultures using glucose. Thus, the strategy described herein allowed for rationally redirecting metabolic fluxes in P. putida from central metabolism towards product biosynthesis—especially relevant when deletion of essential pathways is not an option.     

metaGEM: reconstruction of genome scale metabolic models directly from metagenomes

Metagenomic analyses of microbial communities have revealed a large degree of interspecies and intraspecies genetic diversity through the reconstruction of metagenome assembled genomes (MAGs). Yet, metabolic modeling efforts mainly rely on reference genomes as the starting point for reconstruction and simulation of genome scale metabolic models (GEMs), neglecting the immense intra- and inter-species diversity present in microbial communities. Here, we present metaGEM (https://github.com/franciscozorrilla/metaGEM), an end-to-end pipeline enabling metabolic modeling of multi-species communities directly from metagenomes. The pipeline automates all steps from the extraction of context-specific prokaryotic GEMs from MAGs to community level flux balance analysis (FBA) simulations. To demonstrate the capabilities of metaGEM, we analyzed 483 samples spanning lab culture, human gut, plant-associated, soil, and ocean metagenomes, reconstructing over 14,000 GEMs. We show that GEMs reconstructed from metagenomes have fully represented metabolism comparable to isolated genomes. We demonstrate that metagenomic GEMs capture intraspecies metabolic diversity and identify potential differences in the progression of type 2 diabetes at the level of gut bacterial metabolic exchanges. Overall, metaGEM enables FBA-ready metabolic model reconstruction directly from metagenomes, provides a resource of metabolic models, and showcases community-level modeling of microbiomes associated with disease conditions allowing generation of mechanistic hypotheses.     

In Silico Models for Dynamic Connected Cell Cultures Mimicking Hepatocyte-Endothelial Cell-Adipocyte Interaction Circle

The biochemistry of a system made up of three kinds of cell is virtually impossible to work out without the use of in silico models. Here, we deal with homeostatic balance phenomena from a metabolic point of view and we present a new computational model merging three single-cell models, already available from our research group: the first model reproduced the metabolic behaviour of a hepatocyte, the second one represented an endothelial cell, and the third one described an adipocyte. Multiple interconnections were created among these three models in order to mimic the main physiological interactions that are known for the examined cell phenotypes. The ultimate aim was to recreate the accomplishment of the homeostatic balance as it was observed for an in vitro connected three-culture system concerning glucose and lipid metabolism in the presence of the medium flow. The whole model was based on a modular approach and on a set of nonlinear differential equations implemented in Simulink, applying Michaelis-Menten kinetic laws and some energy balance considerations to the studied metabolic pathways. Our in silico model was then validated against experimental datasets coming from literature about the cited in vitro model. The agreement between simulated and experimental results was good and the behaviour of the connected culture system was reproduced through an adequate parameter evaluation. The developed model may help other researchers to investigate further about integrated metabolism and the regulation mechanisms underlying the physiological homeostasis.     

多组学前沿技术专刊序言

后基因组时代,基因组学、转录组学、蛋白质组学及代谢组学等技术应用日趋广泛,功能注释成为生命科学研究的中心任务,多组学整合分析成为全面解析生物学机理的主要手段。本专刊邀请了国内多组学领域的相关专家学者介绍了基因组学、转录组学、蛋白质组学及代谢组学等领域最新进展和应用成果,收录了相关文章28篇,以供从事多组学研究的科研工作者参考。     

iSCHRUNK – In Silico Approach to Characterization and Reduction of Uncertainty in the Kinetic Models of Genome-scale Metabolic Networks

Accurate determination of physiological states of cellular metabolism requires detailed information about metabolic fluxes, metabolite concentrations and distribution of enzyme states. Integration of fluxomics and metabolomics data, and thermodynamics-based metabolic flux analysis contribute to improved understanding of steady-state properties of metabolism. However, knowledge about kinetics and enzyme activities though essential for quantitative understanding of metabolic dynamics remains scarce and involves uncertainty. Here, we present a computational methodology that allow us to determine and quantify the kinetic parameters that correspond to a certain physiology as it is described by a given metabolic flux profile and a given metabolite concentration vector. Though we initially determine kinetic parameters that involve a high degree of uncertainty, through the use of kinetic modeling and machine learning principles we are able to obtain more accurate ranges of kinetic parameters, and hence we are able to reduce the uncertainty in the model analysis. We computed the distribution of kinetic parameters for glucose-fed E. coli producing 1,4-butanediol and we discovered that the observed physiological state corresponds to a narrow range of kinetic parameters of only a few enzymes, whereas the kinetic parameters of other enzymes can vary widely. Furthermore, this analysis suggests which are the enzymes that should be manipulated in order to engineer the reference state of the cell in a desired way. The proposed approach also sets up the foundations of a novel type of approaches for efficient, non-asymptotic, uniform sampling of solution spaces.     

Whole-genome metabolic model of <Emphasis Type="Italic">Trichoderma reesei</Emphasis> built by comparative reconstruction

Background

Trichoderma reesei is one of the main sources of biomass-hydrolyzing enzymes for the biotechnology industry. There is a need for improving its enzyme production efficiency. The use of metabolic modeling for the simulation and prediction of this organism’s metabolism is potentially a valuable tool for improving its capabilities. An accurate metabolic model is needed to perform metabolic modeling analysis.

Results

A whole-genome metabolic model of T. reesei has been reconstructed together with metabolic models of 55 related species using the metabolic model reconstruction algorithm CoReCo. The previously published CoReCo method has been improved to obtain better quality models. The main improvements are the creation of a unified database of reactions and compounds and the use of reaction directions as constraints in the gap-filling step of the algorithm. In addition, the biomass composition of T. reesei has been measured experimentally to build and include a specific biomass equation in the model.

Conclusions

The improvements presented in this work on the CoReCo pipeline for metabolic model reconstruction resulted in higher-quality metabolic models compared with previous versions. A metabolic model of T. reesei has been created and is publicly available in the BIOMODELS database. The model contains a biomass equation, reaction boundaries and uptake/export reactions which make it ready for simulation. To validate the model, we dem1onstrate that the model is able to predict biomass production accurately and no stoichiometrically infeasible yields are detected. The new T. reesei model is ready to be used for simulations of protein production processes.

     

Flux balance analysis of the ammonia-oxidizing bacterium Nitrosomonas europaea ATCC19718 unravels specific metabolic activities while degrading toxic compounds

The ammonia-oxidizing bacterium Nitrosomonas europaea has been widely recognized as an important player in the nitrogen cycle as well as one of the most abundant members in microbial communities for the treatment of industrial or sewage wastewater. Its natural metabolic versatility and extraordinary ability to degrade environmental pollutants (e.g., aromatic hydrocarbons such as benzene and toluene) enable it to thrive under various harsh environmental conditions. Constraint-based metabolic models constructed from genome sequences enable quantitative insight into the central and specialized metabolism within a target organism. These genome-scale models have been utilized to understand, optimize, and design new strategies for improved bioprocesses. Reduced modeling approaches have been used to elucidate Nitrosomonas europaea metabolism at a pathway level. However, genome-scale knowledge about the simultaneous oxidation of ammonia and pollutant metabolism of N. europaea remains limited. Here, we describe the reconstruction, manual curation, and validation of the genome-scale metabolic model for N. europaea, iGC535. This reconstruction is the most accurate metabolic model for a nitrifying organism to date, reaching an average prediction accuracy of over 90% under several growth conditions. The manually curated model can predict phenotypes under chemolithotrophic and chemolithoorganotrophic conditions while oxidating methane and wastewater pollutants. Calculated flux distributions under different trophic conditions show that several key pathways are affected by the type of carbon source available, including central carbon metabolism and energy production.