Consistent Estimation of Gibbs Energy Using Component Contributions

Standard Gibbs energies of reactions are increasingly being used in metabolic modeling for applying thermodynamic constraints on reaction rates, metabolite concentrations and kinetic parameters. The increasing scope and diversity of metabolic models has led scientists to look for genome-scale solutions that can estimate the standard Gibbs energy of all the reactions in metabolism. Group contribution methods greatly increase coverage, albeit at the price of decreased precision. We present here a way to combine the estimations of group contribution with the more accurate reactant contributions by decomposing each reaction into two parts and applying one of the methods on each of them. This method gives priority to the reactant contributions over group contributions while guaranteeing that all estimations will be consistent, i.e. will not violate the first law of thermodynamics. We show that there is a significant increase in the accuracy of our estimations compared to standard group contribution. Specifically, our cross-validation results show an 80% reduction in the median absolute residual for reactions that can be derived by reactant contributions only. We provide the full framework and source code for deriving estimates of standard reaction Gibbs energy, as well as confidence intervals, and believe this will facilitate the wide use of thermodynamic data for a better understanding of metabolism

Comparative evaluation of atom mapping algorithms for balanced metabolic reactions: application to Recon 3D

The mechanism of each chemical reaction in a metabolic network can be represented as a set of atom mappings, each of which relates an atom in a substrate metabolite to an atom of the same element in a product metabolite. Genome-scale metabolic network reconstructions typically represent biochemistry at the level of reaction stoichiometry. However, a more detailed representation at the underlying level of atom mappings opens the possibility for a broader range of biological, biomedical and biotechnological applications than with stoichiometry alone. Complete manual acquisition of atom mapping data for a genome-scale metabolic network is a laborious process. However, many algorithms exist to predict atom mappings. How do their predictions compare to each other and to manually curated atom mappings? For more than four thousand metabolic reactions in the latest human metabolic reconstruction, Recon 3D, we compared the atom mappings predicted by six atom mapping algorithms. We also compared these predictions to those obtained by manual curation of atom mappings for over five hundred reactions distributed among all top level Enzyme Commission number classes. Five of the evaluated algorithms had similarly high prediction accuracy of over 91% when compared to manually curated atom mapped reactions. On average, the accuracy of the prediction was highest for reactions catalysed by oxidoreductases and lowest for reactions catalysed by ligases. In addition to prediction accuracy, the algorithms were evaluated on their accessibility, their advanced features, such as the ability to identify equivalent atoms, and their ability to map hydrogen atoms. In addition to prediction accuracy, we found that software accessibility and advanced features were fundamental to the selection of an atom mapping algorithm in practice

The Virtual Metabolic Human database: integrating human and gut microbiome metabolism with nutrition and disease

A multitude of factors contribute to complex diseases and can be measured with ‘omics’ methods. Databases facilitate data interpretation for underlying mechanisms. Here, we describe the Virtual Metabolic Human (VMH, www.vmh.life) database encapsulating current knowledge of human metabolism within five interlinked resources ‘Human metabolism’, ‘Gut microbiome’, ‘Disease’, ‘Nutrition’, and ‘ReconMaps’. The VMH captures 5180 unique metabolites, 17 730 unique reactions, 3695 human genes, 255 Mendelian diseases, 818 microbes, 632 685 microbial genes and 8790 food items. The VMH’s unique features are (i) the hosting of the metabolic reconstructions of human and gut microbes amenable for metabolic modeling; (ii) seven human metabolic maps for data visualization; (iii) a nutrition designer; (iv) a user-friendly webpage and application-programming interface to access its content; (v) user feedback option for community engagement and (vi) the connection of its entities to 57 other web resources. The VMH represents a novel, interdisciplinary database for data interpretation and hypothesis generation to the biomedical community

Quantitative Assignment of Reaction Directionality in a Multicompartmental Human Metabolic Reconstruction

Reaction directionality is a key constraint in the modeling of genome-scale metabolic networks. We thermodynamically constrained reaction directionality in a multicompartmental genome-scale model of human metabolism, Recon 1, by calculating, in vivo, standard transformed reaction Gibbs energy as a function of compartment-specific pH, electrical potential, and ionic strength. We show that compartmental pH is an important determinant of thermodynamically determined reaction directionality. The effects of pH on transport reaction thermodynamics are only seen to their full extent when metabolites are represented as pseudoisomer groups of multiple protonated species. We accurately predict the irreversibility of 387 reactions, with detailed propagation of uncertainty in input data, and manually curate the literature to resolve conflicting directionality assignments. In at least half of all cases, a prediction of a reversible reaction directionality is due to the paucity of compartment-specific quantitative metabolomic data, with remaining cases due to uncertainty in estimation of standard reaction Gibbs energy. This study points to the pressing need for 1), quantitative metabolomic data, and 2), experimental measurement of thermochemical properties for human metabolites

Comparative evaluation of open source software for mapping between metabolite identifiers in metabolic network reconstructions: application to Recon 2

BACKGROUND: An important step in the reconstruction of a metabolic network is annotation of metabolites. Metabolites are generally annotated with various database or structure based identifiers. Metabolite annotations in metabolic reconstructions may be incorrect or incomplete and thus need to be updated prior to their use. Genome-scale metabolic reconstructions generally include hundreds of metabolites. Manually updating annotations is therefore highly laborious. This prompted us to look for open-source software applications that could facilitate automatic updating of annotations by mapping between available metabolite identifiers. We identified three applications developed for the metabolomics and chemical informatics communities as potential solutions. The applications were MetMask, the Chemical Translation System, and UniChem. The first implements a “metabolite masking” strategy for mapping between identifiers whereas the latter two implement different versions of an InChI based strategy. Here we evaluated the suitability of these applications for the task of mapping between metabolite identifiers in genome-scale metabolic reconstructions. We applied the best suited application to updating identifiers in Recon 2, the latest reconstruction of human metabolism. RESULTS: All three applications enabled partially automatic updating of metabolite identifiers, but significant manual effort was still required to fully update identifiers. We were able to reduce this manual effort by searching for new identifiers using multiple types of information about metabolites. When multiple types of information were combined, the Chemical Translation System enabled us to update over 3,500 metabolite identifiers in Recon 2. All but approximately 200 identifiers were updated automatically. CONCLUSIONS: We found that an InChI based application such as the Chemical Translation System was better suited to the task of mapping between metabolite identifiers in genome-scale metabolic reconstructions. We identified several features, however, that could be added to such an application in order to tailor it to this task

Modeling the effects of commonly used drugs on human metabolism

Metabolism contributes significantly to the pharmacokinetics and pharmacodynamics of a drug. In addition, diet and genetics have a profound effect on cellular metabolism under health and disease. Herein, we assembled a comprehensive, literature-based drug metabolic reconstruction of the 18 most highly prescribed drug groups including statins, antihypertensives, immunosuppressants, and analgesics. This reconstruction captures in detail our current understanding of their absorption, intra-cellular distribution, metabolism, and elimination. We combined this drug module with the most comprehensive reconstruction of human metabolism, Recon 2, yielding Recon2_DM1796, which accounts for 2803 metabolites and 8161 reactions. By defining 50 specific drug objectives that captured the overall drug metabolism of these compounds, we investigated effects of dietary composition and inherited metabolic disorders on drug metabolism and drug-drug interactions. Our main findings include (i) shift in dietary patterns significantly affect statins and acetaminophen metabolism, (ii) disturbed statin metabolism contributes to the clinical phenotype of mitochondrial energy disorders, and (iii) the interaction between statins and cyclosporine can be explained by several common metabolic and transport pathways other than the previously established CYP3A4 connection. This work holds the potential for studying adverse drug reactions and designing patient-specific therapies

Comparing prenatal screening experiences of Icelandic women who received false-positive and true-negative first-trimester combined screening results in Iceland in 2012-2016.

To access publisher's full text version of this article click on the hyperlink belowFirst-trimester combined screening (FTS) has been offered to all pregnant women in Iceland since 2003. Individuals with high-risk FTS results are offered an invasive test option with a ≤1% risk of fetal loss. This study gives insight into the prenatal screening and diagnosis experiences and preferences of 101 women who underwent FTS in Iceland in the years 2012-2016, comparing the experience of those who received false-positive FTS results to those who received true-negative results. Retrospective patient-reported anxiety levels at the time of receiving FTS results were significantly higher in those who received false-positive results compared to those who received true-negative results. For a subset of these participants, the anxiety lasted through pregnancy, and for a smaller subset, it lasted even longer. Non-invasive prenatal testing (NIPT) is currently not offered in Iceland, aside from the rare exceptional case. Given the extremely low false-positive rates of NIPT, we believe NIPT is worth considering as Iceland's standard first-tier screening method for trisomy 13, 18, and 21. We believe the findings of this study are beneficial not only for Iceland but also for other countries where FTS is the first-tier prenatal screening method or the only offered test. Additionally, only 21% of participants in our study reported that they had heard of NIPT, which emphasizes the need for comprehensive NIPT pretest information to be available prior to its uptake to ensure informed and autonomous decision-making.Department of Genetics and Molecular Medicine at Landspitali University Hospita

CHRR: Coordi- nate hit-and-run with rounding for uniform sampling of constraint-based models

In constraint-based metabolic modelling, physical and biochemical constraints define a polyhedral convex set of feasible flux vectors. Uniform sampling of this set provides an unbiased characterization of the metabolic capabilities of a biochemical network. However, reliable uniform sampling of genome-scale biochemical networks is challenging due to their high dimensionality and inherent anisotropy. Here, we present an implementation of a new sampling algorithm, coordinate hit-and-run with rounding (CHRR). This algorithm is based on the provably efficient hit-and-run random walk and crucially uses a preprocessing step to round the anisotropic flux set. CHRR provably converges to a uniform stationary sampling distribution. We apply it to metabolic networks of increasing dimensionality. We show that it converges several times faster than a popular artificial centering hit-and-run algorithm, enabling reliable and tractable sampling of genome-scale biochemical networks.FWN – Publicaties zonder aanstelling Universiteit Leide

Krabbamein í ristli og endaþarmi - yfirlitsgrein

Colorectal cancer is the third most common cancer in the Western hemisphere and the incidence increases with increasing age. Most colorectal cancers are localized with or without lymph node metastases. Up to 20% of patients present with metastatic disease, most commonly to the liver. Surgery is the only curative therapy for localized colorectal cancer and adjuvant chemotherapy is usually recommended for patients with lymph node metastases. Surgery, radiation therapy and chemotherapy are the key components of rectal cancer therapy. Selected patients with recurrent and metastatic disease can be salvaged with surgery but chemotherapy remains the mainstay of therapy for advanced colorectal cancer. Substantial progress has been observed in the treatment of metastatic colorectal cancer in recent years

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

Motivation: 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.Results: 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.Availability and implementation: 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