Temperature-Dependent Estimation of Gibbs Energies Using an Updated Group-Contribution Method

Multi-scale reactor designs extend the physical limits of CO 2 fixation

CO 2 valorization is a promising strategy for climate adaptation and transitioning towards a circular carbon economy. Here, we present a multi-scale, integrated systems approach for designing biomanufacturing systems that can utilizeCO 2 as a feedstock, focusing on the Wood-Ljungdahl and reductive glycine pathways. This approach relies on first principles, coupling the optimization of pathway and process variables. We examine theCO 2 -fixation capacity of both pathways in single- and multi-compartment reactor systems, demonstrating that the reductive glycine pathway has the potential to fixCO 2 at significantly higher rates than photosynthetic organisms. We show that small differences in the energy-dissipative and stoichiometric structures of carbon-fixation pathways could significantly impact optimal designs and feasible design spaces. Our first-principle, systems-level approach quantifies these differences and uncovers strategies to expand the design space and extend the physical limits of carbon fixation, offering insights into pathway selection and process configurations for efficient biomanufacturing.

Bottom-up parameterization of enzyme rate constants: Reconciling inconsistent data

Kinetic models of metabolism are promising platforms for studying complex metabolic systems and designing production strains. Given the availability of enzyme kinetic data from historical experiments and machine learning estimation tools, a straightforward modeling approach is to assemble kinetic data enzyme by enzyme until a desired scale is reached. However, this type of 'bottom up' parameterization of kinetic models has been difficult due to a number of issues including gaps in kinetic parameters, the complexity of enzyme mechanisms, inconsistencies between parameters obtained from different sources, and in vitro-in vivo differences. Here, we present a computational workflow for the robust estimation of kinetic parameters for detailed mass action enzyme models while taking into account parameter uncertainty. The resulting software package, termed MASSef (the Mass Action Stoichiometry Simulation Enzyme Fitting package), can handle standard 'macroscopic' kinetic parameters, including Km, kcat, Ki, Keq, and nh, as well as diverse reaction mechanisms defined in terms of mass action reactions and 'microscopic' rate constants. We provide three enzyme case studies demonstrating that this approach can identify and reconcile inconsistent data either within in vitro experiments or between in vitro and in vivo enzyme function. We further demonstrate how parameterized enzyme modules can be used to assemble pathway-scale kinetic models consistent with in vivo behavior. This work builds on the legacy of knowledge on kinetic behavior of enzymes by enabling robust parameterization of enzyme kinetic models at scale utilizing the abundance of historical literature data and machine learning parameter estimates.

Improving pathway prediction accuracy of constraints-based metabolic network models by treating enzymes as microcompartments

Metabolic network models have become increasingly precise and accurate as the most widespread and practical digital representations of living cells. The prediction functions were significantly expanded by integrating cellular resources and abiotic constraints in recent years. However, if unreasonable modeling methods were adopted due to a lack of consideration of biological knowledge, the conflicts between stoichiometric and other constraints, such as thermodynamic feasibility and enzyme resource availability, would lead to distorted predictions. In this work, we investigated a prediction anomaly of EcoETM, a constraints-based metabolic network model, and introduced the idea of enzyme compartmentalization into the analysis process. Through rational combination of reactions, we avoid the false prediction of pathway feasibility caused by the unrealistic assumption of free intermediate metabolites. This allowed us to correct the pathway structures of l-serine and l-tryptophan. A specific analysis explains the application method of the EcoETM-like model and demonstrates its potential and value in correcting the prediction results in pathway structure by resolving the conflict between different constraints and incorporating the evolved roles of enzymes as reaction compartments. Notably, this work also reveals the trade-off between product yield and thermodynamic feasibility. Our work is of great value for the structural improvement of constraints-based models.

Warm Cells, Hot Mitochondria: Achievements and Problems of Ultralocal Thermometry

Temperature is a crucial regulator of the rate and direction of biochemical reactions and cell processes. The recent data indicating the presence of local thermal gradients associated with the sites of high-rate thermogenesis, on the one hand, demonstrate the possibility for the existence of "thermal signaling" in a cell and, on the other, are criticized on the basis of thermodynamic calculations and models. Here, we review the main thermometric techniques and sensors developed for the determination of temperature inside living cells and diverse intracellular compartments. A comparative analysis is conducted of the results obtained using these methods for the cytosol, nucleus, endo-/sarcoplasmic reticulum, and mitochondria, as well as their biological consistency. Special attention is given to the limitations, possible sources of errors and ambiguities of the sensor's responses. The issue of biological temperature limits in cells and organelles is considered. It is concluded that the elaboration of experimental protocols for ultralocal temperature measurements that take into account both the characteristics of biological systems, as well as the properties and limitations of each type of sensor is of critical importance for the generation of reliable results and further progress in this field.

Model-driven insights into the effects of temperature on metabolism

Temperature affects cellular processes at different spatiotemporal scales, and identifying the genetic and molecular mechanisms underlying temperature responses paves the way to develop approaches for mitigating the effects of future climate scenarios. A systems view of the effects of temperature on cellular physiology can be obtained by focusing on metabolism since: (i) its functions depend on transcription and translation and (ii) its outcomes support organisms' development, growth, and reproduction. Here we provide a systematic review of modelling efforts directed at investigating temperature effects on properties of single biochemical reactions, system-level traits, metabolic subsystems, and whole-cell metabolism across different prokaryotes and eukaryotes. We compare and contrast computational approaches and theories that facilitate modelling of temperature effects on key properties of enzymes and their consideration in constraint-based as well as kinetic models of metabolism. In addition, we provide a summary of insights from computational approaches, facilitating integration of omics data from temperature-modulated experiments with models of metabolic networks, and review the resulting biotechnological applications. Lastly, we provide a perspective on how different types of metabolic modelling can profit from developments in machine learning and models of different cellular layers to improve model-driven insights into the effects of temperature relevant for biotechnological applications.

Real-Time Temperature Sensing Using a Ratiometric Dual Fluorescent Protein Biosensor

Accurate temperature control within biological and chemical reaction samples and instrument calibration are essential to the diagnostic, pharmaceutical and chemical industries. This is particularly challenging for microlitre-scale reactions typically used in real-time PCR applications and differential scanning fluorometry. Here, we describe the development of a simple, inexpensive ratiometric dual fluorescent protein temperature biosensor (DFPTB). A combination of cycle three green fluorescent protein and a monomeric red fluorescent protein enabled the quantification of relative temperature changes and the identification of temperature discrepancies across a wide temperature range of 4-70 °C. The maximal sensitivity of 6.7% °C^-1 and precision of 0.1 °C were achieved in a biologically relevant temperature range of 25-42 °C in standard phosphate-buffered saline conditions at a pH of 7.2. Good temperature sensitivity was achieved in a variety of biological buffers and pH ranging from 4.8 to 9.1. The DFPTB can be used in either purified or mixed bacteria-encapsulated formats, paving the way for in vitro and in vivo applications for topologically precise temperature measurements.

Reconstruction and Analysis of Thermodynamically Constrained Models Reveal Metabolic Responses of a Deep-Sea Bacterium to Temperature Perturbations

Microbial acclimation to different temperature conditions can involve broad changes in cell composition and metabolic efficiency. A systems-level view of these metabolic responses in nonmesophilic organisms, however, is currently missing. In this study, thermodynamically constrained genome-scale models were applied to simulate the metabolic responses of a deep-sea psychrophilic bacterium, Shewanella psychrophila WP2, under suboptimal (4°C), optimal (15°C), and supraoptimal (20°C) growth temperatures. The models were calibrated with experimentally determined growth rates of WP2. Gibbs free energy change of reactions (Δ_rG'), metabolic fluxes, and metabolite concentrations were predicted using random simulations to characterize temperature-dependent changes in the metabolism. The modeling revealed the highest metabolic efficiency at the optimal temperature, and it suggested distinct patterns of ATP production and consumption that could lead to lower metabolic efficiency under suboptimal or supraoptimal temperatures. The modeling also predicted rearrangement of fluxes through multiple metabolic pathways, including the glycolysis pathway, Entner-Doudoroff pathway, tricarboxylic acid (TCA) cycle, and electron transport system, and these predictions were corroborated through comparisons to WP2 transcriptomes. Furthermore, predictions of metabolite concentrations revealed the potential conservation of reducing equivalents and ATP in the suboptimal temperature, consistent with experimental observations from other psychrophiles. Taken together, the WP2 models provided mechanistic insights into the metabolism of a psychrophile in response to different temperatures. IMPORTANCE Metabolic flexibility is a central component of any organism's ability to survive and adapt to changes in environmental conditions. This study represents the first application of thermodynamically constrained genome-scale models in simulating the metabolic responses of a deep-sea psychrophilic bacterium to various temperatures. The models predicted differences in metabolic efficiency that were attributed to changes in metabolic pathway utilization and metabolite concentration during growth under optimal and nonoptimal temperatures. Experimental growth measurements were used for model calibration, and temperature-dependent transcriptomic changes corroborated the model-predicted rearrangement of metabolic fluxes. Overall, this study highlights the utility of modeling approaches in studying the temperature-driven metabolic responses of an extremophilic organism.

eQuilibrator 3.0: a database solution for thermodynamic constant estimation

eQuilibrator (equilibrator.weizmann.ac.il) is a database of biochemical equilibrium constants and Gibbs free energies, originally designed as a web-based interface. While the website now counts around 1,000 distinct monthly users, its design could not accommodate larger compound databases and it lacked a scalable Application Programming Interface (API) for integration into other tools developed by the systems biology community. Here, we report on the recent updates to the database as well as the addition of a new Python-based interface to eQuilibrator that adds many new features such as a 100-fold larger compound database, the ability to add novel compounds, improvements in speed and memory use, and correction for Mg²⁺ ion concentrations. Moreover, the new interface can compute the covariance matrix of the uncertainty between estimates, for which we show the advantages and describe the application in metabolic modelling. We foresee that these improvements will make thermodynamic modelling more accessible and facilitate the integration of eQuilibrator into other software platforms.

dGPredictor: Automated fragmentation method for metabolic reaction free energy prediction and de novo pathway design

Group contribution (GC) methods are conventionally used in thermodynamics analysis of metabolic pathways to estimate the standard Gibbs energy change (Δ_rG′^o) of enzymatic reactions from limited experimental measurements. However, these methods are limited by their dependence on manually curated groups and inability to capture stereochemical information, leading to low reaction coverage. Herein, we introduce an automated molecular fingerprint-based thermodynamic analysis tool called dGPredictor that enables the consideration of stereochemistry within metabolite structures and thus increases reaction coverage. dGPredictor has comparable prediction accuracy compared to existing GC methods and can capture Gibbs energy changes for isomerase and transferase reactions, which exhibit no overall group changes. We also demonstrate dGPredictor’s ability to predict the Gibbs energy change for novel reactions and seamless integration within de novo metabolic pathway design tools such as novoStoic for safeguarding against the inclusion of reaction steps with infeasible directionalities. To facilitate easy access to dGPredictor, we developed a graphical user interface to predict the standard Gibbs energy change for reactions at various pH and ionic strengths. The tool allows customized user input of known metabolites as KEGG IDs and novel metabolites as InChI strings (https://github.com/maranasgroup/dGPredictor).

The standard Gibbs energy change is commonly used to check for the feasibility of enzyme-catalyzed reactions as thermodynamics plays a crucial role in pathway design for biochemical synthesis. The group contribution methods using expert-defined functional groups have been extensively used for estimating standard Gibbs energy change. Here, we introduce a molecular fingerprint-based thermodynamic tool, dGPredictor, that enables distinguishing between (stereo)isomers in metabolic reactions leading to improved reaction coverage and comparable prediction accuracy as GC methods. dGPredictor can also be used alongside de novo pathway design tools to ensure the correct directionality of chosen reaction steps. We applied and tested dGPredictor on reactions from the KEGG database and applied it to screen an isobutanol synthesis pathway design. An open-source, user-friendly web interface is provided to facilitate easy access for standard Gibbs energy change of reactions at different pH values. (https://github.com/maranasgroup/dGPredictor).

Integrating thermodynamic and enzymatic constraints into genome-scale metabolic models

Stoichiometric genome-scale metabolic network models (GEMs) have been widely used to predict metabolic phenotypes. In addition to stoichiometric ratios, other constraints such as enzyme availability and thermodynamic feasibility can also limit the phenotype solution space. Extended GEM models considering either enzymatic or thermodynamic constraints have been shown to improve prediction accuracy. In this paper, we propose a novel method that integrates both enzymatic and thermodynamic constraints in a single Pyomo modeling framework (ETGEMs). We applied this method to construct the EcoETM (E. coli metabolic model with enzymatic and thermodynamic constraints). Using this model, we calculated the optimal pathways for cellular growth and the production of 22 metabolites. When comparing the results with those of iML1515 and models with one of the two constraints, we observed that many thermodynamically unfavorable and/or high enzyme cost pathways were excluded from EcoETM. For example, the synthesis pathway of carbamoyl-phosphate (Cbp) from iML1515 is both thermodynamically unfavorable and enzymatically costly. After introducing the new constraints, the production pathways and yields of several Cbp-derived products (e.g. L-arginine, orotate) calculated using EcoETM were more realistic. The results of this study demonstrate the great application potential of metabolic models with multiple constraints for pathway analysis and phenotype prediction.

Quantum Mechanical Methods Predict Accurate Thermodynamics of Biochemical Reactions

Thermodynamics plays a crucial role in regulating the metabolic processes in all living organisms. Accurate determination of biochemical and biophysical properties is important to understand, analyze, and synthetically design such metabolic processes for engineered systems. In this work, we extensively performed first-principles quantum mechanical calculations to assess its accuracy in estimating free energy of biochemical reactions and developed automated quantum-chemistry (QC) pipeline (https://appdev.kbase.us/narrative/45710) for the prediction of thermodynamics parameters of biochemical reactions. We benchmark the QC methods based on density functional theory (DFT) against different basis sets, solvation models, pH, and exchange-correlation functionals using the known thermodynamic properties from the NIST database. Our results show that QC calculations when combined with simple calibration yield a mean absolute error in the range of 1.60-2.27 kcal/mol for different exchange-correlation functionals, which is comparable to the error in the experimental measurements. This accuracy over a diverse set of metabolic reactions is unprecedented and near the benchmark chemical accuracy of 1 kcal/mol that is usually desired from DFT calculations.

Probabilistic thermodynamic analysis of metabolic networks

MOTIVATION

Random sampling of metabolic fluxes can provide a comprehensive description of the capabilities of a metabolic network. However, current sampling approaches do not model thermodynamics explicitly, leading to inaccurate predictions of an organism's potential or actual metabolic operations.

RESULTS

We present a probabilistic framework combining thermodynamic quantities with steady-state flux constraints to analyze the properties of a metabolic network. It includes methods for probabilistic metabolic optimization and for joint sampling of thermodynamic and flux spaces. Applied to a model of Escherichia coli, we use the methods to reveal known and novel mechanisms of substrate channeling, and to accurately predict reaction directions and metabolite concentrations. Interestingly, predicted flux distributions are multimodal, leading to discrete hypotheses on E.coli's metabolic capabilities.

AVAILABILITY AND IMPLEMENTATION

Python and MATLAB packages available at https://gitlab.com/csb.ethz/pta.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Modelling Cell Metabolism: A Review on Constraint-Based Steady-State and Kinetic Approaches

Studying cell metabolism serves a plethora of objectives such as the enhancement of bioprocess performance, and advancement in the understanding of cell biology, of drug target discovery, and in metabolic therapy. Remarkable successes in these fields emerged from heuristics approaches, for instance, with the introduction of effective strategies for genetic modifications, drug developments and optimization of bioprocess management. However, heuristics approaches have showed significant shortcomings, such as to describe regulation of metabolic pathways and to extrapolate experimental conditions. In the specific case of bioprocess management, such shortcomings limit their capacity to increase product quality, while maintaining desirable productivity and reproducibility levels. For instance, since heuristics approaches are not capable of prediction of the cellular functions under varying experimental conditions, they may lead to sub-optimal processes. Also, such approaches used for bioprocess control often fail in regulating a process under unexpected variations of external conditions. Therefore, methodologies inspired by the systematic mathematical formulation of cell metabolism have been used to address such drawbacks and achieve robust reproducible results. Mathematical modelling approaches are effective for both the characterization of the cell physiology, and the estimation of metabolic pathways utilization, thus allowing to characterize a cell population metabolic behavior. In this article, we present a review on methodology used and promising mathematical modelling approaches, focusing primarily to investigate metabolic events and regulation. Proceeding from a topological representation of the metabolic networks, we first present the metabolic modelling approaches that investigate cell metabolism at steady state, complying to the constraints imposed by mass conservation law and thermodynamics of reactions reversibility. Constraint-based models (CBMs) are reviewed highlighting the set of assumed optimality functions for reaction pathways. We explore models simulating cell growth dynamics, by expanding flux balance models developed at steady state. Then, discussing a change of metabolic modelling paradigm, we describe dynamic kinetic models that are based on the mathematical representation of the mechanistic description of nonlinear enzyme activities. In such approaches metabolic pathway regulations are considered explicitly as a function of the activity of other components of metabolic networks and possibly far from the metabolic steady state. We have also assessed the significance of metabolic model parameterization in kinetic models, summarizing a standard parameter estimation procedure frequently employed in kinetic metabolic modelling literature. Finally, some optimization practices used for the parameter estimation are reviewed.

Physicochemical and metabolic constraints for thermodynamics-based stoichiometric modelling under mesophilic growth conditions

Metabolic engineering in the post-genomic era is characterised by the development of new methods for metabolomics and fluxomics, supported by the integration of genetic engineering tools and mathematical modelling. Particularly, constraint-based stoichiometric models have been widely studied: (i) flux balance analysis (FBA) (in silico), and (ii) metabolic flux analysis (MFA) (in vivo). Recent studies have enabled the incorporation of thermodynamics and metabolomics data to improve the predictive capabilities of these approaches. However, an in-depth comparison and evaluation of these methods is lacking. This study presents a thorough analysis of two different in silico methods tested against experimental data (metabolomics and ¹³C-MFA) for the mesophile Escherichia coli. In particular, a modified version of the recently published matTFA toolbox was created, providing a broader range of physicochemical parameters. Validating against experimental data allowed the determination of the best physicochemical parameters to perform the TFA (Thermodynamics-based Flux Analysis). An analysis of flux pattern changes in the central carbon metabolism between ¹³C-MFA and TFA highlighted the limited capabilities of both approaches for elucidating the anaplerotic fluxes. In addition, a method based on centrality measures was suggested to identify important metabolites that (if quantified) would allow to further constrain the TFA. Finally, this study emphasised the need for standardisation in the fluxomics community: novel approaches are frequently released but a thorough comparison with currently accepted methods is not always performed.

Biotechnology has benefitted from the development of high throughput methods characterising living systems at different levels (e.g. concerning genes or proteins), allowing the industrial production of chemical commodities. Recently, focus has been placed on determining reaction rates (or metabolic fluxes) in the metabolic network of certain microorganisms, in order to identify bottlenecks hindering their exploitation. Two main approaches are commonly used, termed metabolic flux analysis (MFA) and flux balance analysis (FBA), based on measuring and estimating fluxes, respectively. While the influence of thermodynamics in living systems was accepted several decades ago, its application to study biochemical networks has only recently been enabled. In this sense, a multitude of different approaches constraining well-established modelling methods with thermodynamics has been suggested. However, physicochemical parameters are generally not properly adjusted to the experimental conditions, which might affect their predictive capabilities. In this study, we have explored the reliability of currently available tools by investigating the impact of varying said parameters in the simulation of metabolic fluxes and metabolite concentration values. Additionally, our in-depth analysis allowed us to highlight limitations and potential solutions that should be considered in future studies.

The Expanding Computational Toolbox for Engineering Microbial Phenotypes at the Genome Scale

Microbial strains are being engineered for an increasingly diverse array of applications, from chemical production to human health. While traditional engineering disciplines are driven by predictive design tools, these tools have been difficult to build for biological design due to the complexity of biological systems and many unknowns of their quantitative behavior. However, due to many recent advances, the gap between design in biology and other engineering fields is closing. In this work, we discuss promising areas of development of computational tools for engineering microbial strains. We define five frontiers of active research: (1) Constraint-based modeling and metabolic network reconstruction, (2) Kinetics and thermodynamic modeling, (3) Protein structure analysis, (4) Genome sequence analysis, and (5) Regulatory network analysis. Experimental and machine learning drivers have enabled these methods to improve by leaps and bounds in both scope and accuracy. Modern strain design projects will require these tools to be comprehensively applied to the entire cell and efficiently integrated within a single workflow. We expect that these frontiers, enabled by the ongoing revolution of big data science, will drive forward more advanced and powerful strain engineering strategies.

Recent advances in constraint and machine learning-based metabolic modeling by leveraging stoichiometric balances, thermodynamic feasibility and kinetic law formalisms

Understanding the governing principles behind organisms' metabolism and growth underpins their effective deployment as bioproduction chassis. A central objective of metabolic modeling is predicting how metabolism and growth are affected by both external environmental factors and internal genotypic perturbations. The fundamental concepts of reaction stoichiometry, thermodynamics, and mass action kinetics have emerged as the foundational principles of many modeling frameworks designed to describe how and why organisms allocate resources towards both growth and bioproduction. This review focuses on the latest algorithmic advancements that have integrated these foundational principles into increasingly sophisticated quantitative frameworks.

Thermodynamic Favorability of End Products of Anaerobic Glucose Metabolism

The eQuilibrator component contribution method allows calculation of the overall Gibbs energy changes for conversion of glucose to a wide range of final products in the absence of other oxidants. Values are presented for all possible combinations of products with up to three carbons and selected others. The most negative Gibbs energy change is for the formation of graphite and water (-499 kJ mol-1) followed by CH4 and CO2 (-430 kJ mol-1), the observed final products of anaerobic digestion. Other favored products (with various combinations having Gibbs energy changes between -300 and -367 kJ mol-1) are short-chain alkanes, fatty acids, dicarboxylic acids, and even hexane and benzene. The most familiar products, lactate and ethanol + CO2, are less favored (Gibbs energy changes of -206 and -265 kJ mol-1 respectively). The values presented offer an interesting perspective on observed metabolism and its evolutionary origins as well as on cells engineered for biotechnological purposes.

Thermodynamically consistent estimation of Gibbs free energy from data: data reconciliation approach

MOTIVATION

Thermodynamic analysis of biological reaction networks requires the availability of accurate and consistent values of Gibbs free energies of reaction and formation. These Gibbs energies can be measured directly via the careful design of experiments or can be computed from the curated Gibbs free energy databases. However, the computed Gibbs free energies of reactions and formations do not satisfy the thermodynamic constraints due to the compounding effect of measurement errors in the experimental data. The propagation of these errors can lead to a false prediction of pathway feasibility and uncertainty in the estimation of thermodynamic parameters.

RESULTS

This work proposes a data reconciliation framework for thermodynamically consistent estimation of Gibbs free energies of reaction, formation and group contributions from experimental data. In this framework, we formulate constrained optimization problems that reduce measurement errors and their effects on the estimation of Gibbs energies such that the thermodynamic constraints are satisfied. When a subset of Gibbs free energies of formations is unavailable, it is shown that the accuracy of their resulting estimates is better than that of existing empirical prediction methods. Moreover, we also show that the estimation of group contributions can be improved using this approach. Further, we provide guidelines based on this approach for performing systematic experiments to estimate unknown Gibbs formation energies.

AVAILABILITY AND IMPLEMENTATION

The MATLAB code for the executing the proposed algorithm is available for free on the GitHub repository: https://github.com/samansalike/DR-thermo.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Near-equilibrium glycolysis supports metabolic homeostasis and energy yield

Glycolysis plays a central role in producing ATP and biomass. Its control principles, however, remain incompletely understood. Here, we develop a method that combines ²H and ¹³C tracers to determine glycolytic thermodynamics. Using this method, we show that, in conditions and organisms with relatively slow fluxes, multiple steps in glycolysis are near to equilibrium, reflecting spare enzyme capacity. In Escherichia coli, nitrogen or phosphorus upshift rapidly increases the thermodynamic driving force, deploying the spare enzyme capacity to increase flux. Similarly, respiration inhibition in mammalian cells rapidly increases both glycolytic flux and the thermodynamic driving force. The thermodynamic shift allows flux to increase with only small metabolite concentration changes. Finally, we find that the cellulose-degrading anaerobe Clostridium cellulolyticum exhibits slow, near-equilibrium glycolysis due to the use of pyrophosphate rather than ATP for fructose-bisphosphate production, resulting in enhanced per-glucose ATP yield. Thus, near-equilibrium steps of glycolysis promote both rapid flux adaptation and energy efficiency.

Interrogating metabolism as an electron flow system

Metabolism is generally considered as a neatly organised system of modular pathways, shaped by evolution under selection for optimal cellular growth. This view falls short of explaining and predicting a number of key observations about the structure and dynamics of metabolism. We highlight these limitations of a pathway-centric view on metabolism and summarise studies suggesting how these could be overcome by viewing metabolism as a thermodynamically and kinetically constrained, dynamical flow system. Such a systems-level, first-principles based view of metabolism can open up new avenues of metabolic engineering and cures for metabolic diseases and allow better insights to a myriad of physiological processes that are ultimately linked to metabolism. Towards further developing this view, we call for a closer interaction among physical and biological disciplines and an increased use of electrochemical and biophysical approaches to interrogate cellular metabolism together with the microenvironment in which it exists.

Machine learning applied to enzyme turnover numbers reveals protein structural correlates and improves metabolic models

Knowing the catalytic turnover numbers of enzymes is essential for understanding the growth rate, proteome composition, and physiology of organisms, but experimental data on enzyme turnover numbers is sparse and noisy. Here, we demonstrate that machine learning can successfully predict catalytic turnover numbers in Escherichia coli based on integrated data on enzyme biochemistry, protein structure, and network context. We identify a diverse set of features that are consistently predictive for both in vivo and in vitro enzyme turnover rates, revealing novel protein structural correlates of catalytic turnover. We use our predictions to parameterize two mechanistic genome-scale modelling frameworks for proteome-limited metabolism, leading to significantly higher accuracy in the prediction of quantitative proteome data than previous approaches. The presented machine learning models thus provide a valuable tool for understanding metabolism and the proteome at the genome scale, and elucidate structural, biochemical, and network properties that underlie enzyme kinetics.

Thermodynamic favorability and pathway yield as evolutionary tradeoffs in biosynthetic pathway choice

The structure of the metabolic network contains myriad organism-specific variations across the tree of life, but the selection basis for pathway choices in different organisms is not well understood. Here, we examined the metabolic capabilities with respect to cofactor use and pathway thermodynamics of all sequenced organisms in the Kyoto Encyclopedia of Genes and Genomes Database. We found that (i) many biomass precursors have alternate synthesis routes that vary substantially in thermodynamic favorability and energy cost, creating tradeoffs that may be subject to selection pressure; (ii) alternative pathways in amino acid synthesis are characteristically distinguished by the use of biosynthetically unnecessary acyl-CoA cleavage; (iii) distinct choices preferring thermodynamic-favorable or cofactor-use-efficient pathways exist widely among organisms; (iv) cofactor-use-efficient pathways tend to have a greater yield advantage under anaerobic conditions specifically; and (v) lysine biosynthesis in particular exhibits temperature-dependent thermodynamics and corresponding differential pathway choice by thermophiles. These findings present a view on the evolution of metabolic network structure that highlights a key role of pathway thermodynamics and cofactor use in determining organism pathway choices.