Bridging substrate intake kinetics and bacterial growth phenotypes with flux balance analysis incorporating proteome allocation

DLTKcat: deep learning-based prediction of temperature-dependent enzyme turnover rates

The enzyme turnover rate, ${k}_{cat}$, quantifies enzyme kinetics by indicating the maximum efficiency of enzyme catalysis. Despite its importance, ${k}_{cat}$ values remain scarce in databases for most organisms, primarily because of the cost of experimental measurements. To predict ${k}_{cat}$ and account for its strong temperature dependence, DLTKcat was developed in this study and demonstrated superior performance (log10-scale root mean squared error = 0.88, R-squared = 0.66) than previously published models. Through two case studies, DLTKcat showed its ability to predict the effects of protein sequence mutations and temperature changes on ${k}_{cat}$ values. Although its quantitative accuracy is not high enough yet to model the responses of cellular metabolism to temperature changes, DLTKcat has the potential to eventually become a computational tool to describe the temperature dependence of biological systems.

Simple Growth–Metabolism Relations Are Revealed by Conserved Patterns of Heat Flow from Cultured Microorganisms

Quantitative analyses of cell replication address the connection between metabolism and growth. Various growth models approximate time-dependent cell numbers in culture media, but physiological implications of the parametrizations are vague. In contrast, isothermal microcalorimetry (IMC) measures with unprecedented sensitivity the heat (enthalpy) release via chemical turnover in metabolizing cells. Hence, the metabolic activity can be studied independently of modeling the time-dependence of cell numbers. Unexpectedly, IMC traces of various origins exhibit conserved patterns when expressed in the enthalpy domain rather than the time domain, as exemplified by cultures of Lactococcus lactis (prokaryote), Trypanosoma congolese (protozoan) and non-growing Brassica napus (plant) cells. The data comply extraordinarily well with a dynamic Langmuir adsorption reaction model of nutrient uptake and catalytic turnover generalized here to the non-constancy of catalytic capacity. Formal relations to Michaelis–Menten kinetics and common analytical growth models are briefly discussed. The proposed formalism reproduces the “life span” of cultured microorganisms from exponential growth to metabolic decline by a succession of distinct metabolic phases following remarkably simple nutrient–metabolism relations. The analysis enables the development of advanced enzyme network models of unbalanced growth and has fundamental consequences for the derivation of toxicity measures and the transferability of metabolic activity data between laboratories.

Resource-allocation constraint governs structure and function of microbial communities in metabolic modeling

Predictive modeling tools for assessing microbial communities are important for realizing transformative capabilities of microbiomes in agriculture, ecology, and medicine. Constraint-based community-scale metabolic modeling is unique in its potential for making mechanistic predictions regarding both the structure and function of microbial communities. However, accessing this potential requires an understanding of key physicochemical constraints, which are typically considered on a per-species basis. What is needed is a means of incorporating global constraints relevant to microbial ecology into community models. Resource-allocation constraint, which describes how limited resources should be distributed to different cellular processes, sets limits on the efficiency of metabolic and ecological processes. In this study, we investigate the implications of resource-allocation constraints in community-scale metabolic modeling through a simple mechanism-agnostic implementation of resource-allocation constraints directly at the flux level. By systematically performing single-, two-, and multi-species growth simulations, we show that resource-allocation constraints are indispensable for predicting the structure and function of microbial communities. Our findings call for a scalable workflow for implementing a mechanistic version of resource-allocation constraints to ultimately harness the full potential of community-scale metabolic modeling tools.

A metabolic modeling platform for the computation of microbial ecosystems in time and space (COMETS)

Genome-scale stoichiometric modeling of metabolism has become a standard systems biology tool for modeling cellular physiology and growth. Extensions of this approach are emerging as a valuable avenue for predicting, understanding and designing microbial communities. Computation of microbial ecosystems in time and space (COMETS) extends dynamic flux balance analysis to generate simulations of multiple microbial species in molecularly complex and spatially structured environments. Here we describe how to best use and apply the most recent version of COMETS, which incorporates a more accurate biophysical model of microbial biomass expansion upon growth, evolutionary dynamics and extracellular enzyme activity modules. In addition to a command-line option, COMETS includes user-friendly Python and MATLAB interfaces compatible with the well-established COBRA models and methods, as well as comprehensive documentation and tutorials. This protocol provides a detailed guideline for installing, testing and applying COMETS to different scenarios, generating simulations that take from a few minutes to several days to run, with broad applicability to microbial communities across biomes and scales.

Integrating Systems and Synthetic Biology to Understand and Engineer Microbiomes

Microbiomes are complex and ubiquitous networks of microorganisms whose seemingly limitless chemical transformations could be harnessed to benefit agriculture, medicine, and biotechnology. The spatial and temporal changes in microbiome composition and function are influenced by a multitude of molecular and ecological factors. This complexity yields both versatility and challenges in designing synthetic microbiomes and perturbing natural microbiomes in controlled, predictable ways. In this review, we describe factors that give rise to emergent spatial and temporal microbiome properties and the meta-omics and computational modeling tools that can be used to understand microbiomes at the cellular and system levels. We also describe strategies for designing and engineering microbiomes to enhance or build novel functions. Throughout the review, we discuss key knowledge and technology gaps for elucidating the networks and deciphering key control points for microbiome engineering, and highlight examples where multiple omics and modeling approaches can be integrated to address these gaps.

Integration of enzyme constraints in a genome-scale metabolic model of Aspergillus niger improves phenotype predictions

Background

Genome-scale metabolic model (GSMM) is a powerful tool for the study of cellular metabolic characteristics. With the development of multi-omics measurement techniques in recent years, new methods that integrating multi-omics data into the GSMM show promising effects on the predicted results. It does not only improve the accuracy of phenotype prediction but also enhances the reliability of the model for simulating complex biochemical phenomena, which can promote theoretical breakthroughs for specific gene target identification or better understanding the cell metabolism on the system level.

Results

Based on the basic GSMM model iHL1210 of Aspergillus niger, we integrated large-scale enzyme kinetics and proteomics data to establish a GSMM based on enzyme constraints, termed a GEM with Enzymatic Constraints using Kinetic and Omics data (GECKO). The results show that enzyme constraints effectively improve the model’s phenotype prediction ability, and extended the model’s potential to guide target gene identification through predicting metabolic phenotype changes of A. niger by simulating gene knockout. In addition, enzyme constraints significantly reduced the solution space of the model, i.e., flux variability over 40.10% metabolic reactions were significantly reduced. The new model showed also versatility in other aspects, like estimating large-scale \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_{{cat}}$$\end{document}kcat values, predicting the differential expression of enzymes under different growth conditions.

Conclusions

This study shows that incorporating enzymes’ abundance information into GSMM is very effective for improving model performance with A. niger. Enzyme-constrained model can be used as a powerful tool for predicting the metabolic phenotype of A. niger by incorporating proteome data. In the foreseeable future, with the fast development of measurement techniques, and more precise and rich proteomics quantitative data being obtained for A. niger, the enzyme-constrained GSMM model will show greater application space on the system level.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-021-01614-2.

A steady-state model of microbial acclimation to substrate limitation

Microbes acclimate to changes in substrate availability by altering the number of transporters on the cell surface, however there is some disagreement on just how. We revisit the physics of substrate uptake and consider the steady-state scenario whereby cells have acclimated to maximize fitness. Flux balance analysis of a stoichiometric model of Escherichia coli was used in conjunction with quantitative proteomics data and molecular modeling of membrane transporters to reconcile these opposing views. An emergent feature of the proposed model is a critical substrate concentration S*, which delineates two rate limits. At concentrations above S*, transporter abundance can be regulated to maintain uptake rates as demanded by maximal growth rates, whereas below S*, uptake rates are strictly diffusion limited. In certain scenarios, the proposed model can take on a qualitatively different shape from the familiar hyperbolic kinetics curves, instead resembling the long-forgotten Blackman kinetics.