Identifying functional metabolic shifts in heart failure with the integration of omics data and a heart-specific, genome-scale model

Deciphering and designing microbial communities by genome-scale metabolic modelling

Microbial communities are shaped by the complex interactions among organisms and the environment. Genome-scale metabolic models (GEMs) can provide deeper insights into the complexity and ecological properties of various microbial communities, revealing their intricate interactions. Many researchers have modified GEMs for the microbial communities based on specific needs. Thus, GEMs need to be comprehensively summarized to better understand the trends in their development. In this review, we summarized the key developments in deciphering and designing microbial communities using different GEMs. A timeline of selected highlights in GEMs indicated that this area is evolving from the single-strain level to the microbial community level. Then, we outlined a framework for constructing GEMs of microbial communities. We also summarized the models and resources of static and dynamic community-level GEMs. We focused on the role of external environmental and intracellular resources in shaping the assembly of microbial communities. Finally, we discussed the key challenges and future directions of GEMs, focusing on the integration of GEMs with quorum sensing mechanisms, microbial ecology interactions, machine learning algorithms, and automatic modeling, all of which contribute to consortia-based applications in different fields.

Modeling cardiomyocyte signaling and metabolism predicts genotype-to-phenotype mechanisms in hypertrophic cardiomyopathy

Familial hypertrophic cardiomyopathy (HCM) is a significant precursor of heart failure and sudden cardiac death, primarily caused by mutations in sarcomeric and structural proteins. Despite the extensive research on the HCM genotype, the complex and context-specific nature of many signaling and metabolic pathways linking the HCM genotype to phenotype has hindered therapeutic advancements for patients. Here, we have developed a computational model of HCM encompassing cardiomyocyte signaling and metabolic networks and their associated interactions. Utilizing a stochastic logic-based ODE approach, we linked cardiomyocyte signaling to the metabolic network through a gene regulatory network and post-translational modifications. We validated the model against published data on activities of signaling species in the HCM context and transcriptomes of two HCM mouse models (i.e., R403Q-αMyHC and R92W-TnT). Our model predicts that HCM mutation induces changes in metabolic functions such as ATP synthase deficiency and a transition from fatty acids to carbohydrate metabolism. The model indicated major shifts in glutamine-related metabolism and increased apoptosis after HCM-induced ATP synthase deficiency. We predicted that the transcription factors STAT, SRF, GATA4, TP53, and FoxO are the key regulators of cardiomyocyte hypertrophy and apoptosis in HCM in alignment with experiments. Moreover, we identified shared (e.g., activation of PGC1α by AMPK, and FHL1 by titin) and context-specific mechanisms (e.g., regulation of Ca2+ sensitivity by titin in HCM patients) that may control genotype-to-phenotype transition in HCM across different species or mutations. We also predicted potential combination drug targets for HCM (e.g., mavacamten plus ROS inhibitors) preventing or reversing HCM phenotype (i.e., hypertrophic growth, apoptosis, and metabolic remodeling) in cardiomyocytes. This study provides new insights into mechanisms linking genotype to phenotype in familial hypertrophic cardiomyopathy and offers a framework for assessing new treatments and exploring variations in HCM experimental models.

Effect of Low-Level Tragus Stimulation on Cardiac Metabolism in Heart Failure with Preserved Ejection Fraction: A Transcriptomics-Based Analysis

Abnormal cardiac metabolism precedes and contributes to structural changes in heart failure. Low-level tragus stimulation (LLTS) can attenuate structural remodeling in heart failure with preserved ejection fraction (HFpEF). The role of LLTS on cardiac metabolism is not known. Dahl salt-sensitive rats of 7 weeks of age were randomized into three groups: low salt (0.3% NaCl) diet (control group; n = 6), high salt diet (8% NaCl) with either LLTS (active group; n = 8), or sham stimulation (sham group; n = 5). Both active and sham groups received the high salt diet for 10 weeks with active LLTS or sham stimulation (20 Hz, 2 mA, 0.2 ms) for 30 min daily for the last 4 weeks. At the endpoint, left ventricular tissue was used for RNA sequencing and transcriptomic analysis. The Ingenuity Pathway Analysis tool (IPA) was used to identify canonical metabolic pathways and upstream regulators. Principal component analysis demonstrated overlapping expression of important metabolic genes between the LLTS, and control groups compared to the sham group. Canonical metabolic pathway analysis showed downregulation of the oxidative phosphorylation (Z-score: -4.707, control vs. sham) in HFpEF and LLTS improved the oxidative phosphorylation (Z-score = -2.309, active vs. sham). HFpEF was associated with the abnormalities of metabolic upstream regulators, including PPARGC1α, insulin receptor signaling, PPARα, PPARδ, PPARGC1β, the fatty acid transporter SLC27A2, and lysine-specific demethylase 5A (KDM5A). LLTS attenuated abnormal insulin receptor and KDM5A signaling. HFpEF is associated with abnormal cardiac metabolism. LLTS, by modulating the functioning of crucial upstream regulators, improves cardiac metabolism and mitochondrial oxidative phosphorylation.

Identifying metabolic adaptations characteristic of cardiotoxicity using paired transcriptomics and metabolomics data integrated with a computational model of heart metabolism

Improvements in the diagnosis and treatment of cancer have revealed long-term side effects of chemotherapeutics, particularly cardiotoxicity. Here, we present paired transcriptomics and metabolomics data characterizing in vitro cardiotoxicity to three compounds: 5-fluorouracil, acetaminophen, and doxorubicin. Standard gene enrichment and metabolomics approaches identify some commonly affected pathways and metabolites but are not able to readily identify metabolic adaptations in response to cardiotoxicity. The paired data was integrated with a genome-scale metabolic network reconstruction of the heart to identify shifted metabolic functions, unique metabolic reactions, and changes in flux in metabolic reactions in response to these compounds. Using this approach, we confirm previously seen changes in the p53 pathway by doxorubicin and RNA synthesis by 5-fluorouracil, we find evidence for an increase in phospholipid metabolism in response to acetaminophen, and we see a shift in central carbon metabolism suggesting an increase in metabolic demand after treatment with doxorubicin and 5-fluorouracil.

Physiological Indirect Response Model to Omics-Powered Quantitative Systems Pharmacology Model

Over the past several decades, mathematical modeling has been applied to increasingly wider scopes of questions in drug development. Accordingly, the range of modeling tools has also been evolving, as showcased by contributions of Jusko and colleagues: from basic pharmacokinetics/pharmacodynamics (PK/PD) modeling to today's platform-based approach of quantitative systems pharmacology (QSP) modeling. Aimed at understanding the mechanism of action of investigational drugs, QSP models characterize systemic effects by incorporating information about cellular signaling networks, which is often represented by omics data. In this perspective, we share a few examples illustrating approaches for the integration of omics into mechanistic QSP modeling. We briefly overview how the evolution of PK/PD modeling into QSP has been accompanied by an increase in available data and the complexity of mathematical methods that integrate it. We discuss current gaps and challenges of integrating omics data into QSP models and propose several potential areas where integrated QSP and omics modeling may benefit drug development.

How reliable are Chinese hamster ovary (CHO) cell genome-scale metabolic models?

Genome-scale metabolic models (GEMs) possess the power to revolutionize bioprocess and cell line engineering workflows thanks to their ability to predict and understand whole-cell metabolism in silico. Despite this potential, it is currently unclear how accurately GEMs can capture both intracellular metabolic states and extracellular phenotypes. Here, we investigate this knowledge gap to determine the reliability of current Chinese hamster ovary (CHO) cell metabolic models. We introduce a new GEM, iCHO2441, and create CHO-S and CHO-K1 specific GEMs. These are compared against iCHO1766, iCHO2048, and iCHO2291. Model predictions are assessed via comparison with experimentally measured growth rates, gene essentialities, amino acid auxotrophies, and 13 C intracellular reaction rates. Our results highlight that all CHO cell models are able to capture extracellular phenotypes and intracellular fluxes, with the updated GEM outperforming the original CHO cell GEM. Cell line-specific models were able to better capture extracellular phenotypes but failed to improve intracellular reaction rate predictions in this case. Ultimately, this work provides an updated CHO cell GEM to the community and lays a foundation for the development and assessment of next-generation flux analysis techniques, highlighting areas for model improvements.

Metabolic modeling of sex-specific liver tissue suggests mechanism of differences in toxicological responses

Male subjects in animal and human studies are disproportionately used for toxicological testing. This discrepancy is evidenced in clinical medicine where females are more likely than males to experience liver-related adverse events in response to xenobiotics. While previous work has shown gene expression differences between the sexes, there is a lack of systems-level approaches to understand the direct clinical impact of these differences. Here, we integrate gene expression data with metabolic network models to characterize the impact of transcriptional changes of metabolic genes in the context of sex differences and drug treatment. We used Tasks Inferred from Differential Expression (TIDEs), a reaction-centric approach to analyzing differences in gene expression, to discover that several metabolic pathways exhibit sex differences including glycolysis, fatty acid metabolism, nucleotide metabolism, and xenobiotics metabolism. When TIDEs is used to compare expression differences in treated and untreated hepatocytes, we find several subsystems with differential expression overlap with the sex-altered pathways such as fatty acid metabolism, purine and pyrimidine metabolism, and xenobiotics metabolism. Finally, using sex-specific transcriptomic data, we create individual and averaged male and female liver models and find differences in the pentose phosphate pathway and other metabolic pathways. These results suggest potential sex differences in the contribution of the pentose phosphate pathway to oxidative stress, and we recommend further research into how these reactions respond to hepatotoxic pharmaceuticals.

A contribution of metabolic engineering to addressing medical problems: Metabolic flux analysis

Metabolic engineering has served as a systematic discipline for industrial biotechnology as it has offered systematic tools and methods for strain development and bioprocess optimization. Because these metabolic engineering tools and methods are concerned with the biological network of a cell with emphasis on metabolic network, they have also been applied to a range of medical problems where better understanding of metabolism has also been perceived to be important. Metabolic flux analysis (MFA) is a unique systematic approach initially developed in the metabolic engineering community, and has proved its usefulness and potential when addressing a range of medical problems. In this regard, this review discusses the contribution of MFA to addressing medical problems. For this, we i) provide overview of the milestones of MFA, ii) define two main branches of MFA, namely constraint-based reconstruction and analysis (COBRA) and isotope-based MFA (iMFA), and iii) present successful examples of their medical applications, including characterizing the metabolism of diseased cells and pathogens, and identifying effective drug targets. Finally, synergistic interactions between metabolic engineering and biomedical sciences are discussed with respect to MFA.

Metabolic modeling of sex-specific tissue predicts mechanisms of differences in toxicological responses

Male subjects in animal and human studies are disproportionately used for toxicological testing. This discrepancy is evidenced in clinical medicine where females are more likely than males to experience liver-related adverse events in response to xenobiotics. While previous work has shown gene expression differences between the sexes, there is a lack of systems-level approaches to understand the direct clinical impact effect of these differences. Here, we integrate gene expression data with metabolic network models to characterize the impact of transcriptional changes of metabolic genes in the context of sex differences and drug treatment. We used Tasks Inferred from Differential Expression (TIDEs), a reaction-centric approach to analyzing differences in gene expression, to discover that androgen, ether lipid, glucocorticoid, tryptophan, and xenobiotic metabolism have more activity in the male liver, and serotonin, melatonin, pentose, glucuronate, and vitamin A metabolism have more activity in the female liver. When TIDEs is used to compare expression differences in treated and untreated hepatocytes, we see little response in those sex-altered subsystems, and the largest differences are in subsystems related to lipid metabolism. Finally, using sex-specific transcriptomic data, we create individual and averaged male and female liver models and find differences in the import of bile acids and salts. This result suggests that the sexually dimorphic behavior of the liver may be caused by differences in enterohepatic recirculation, and we suggest an investigation into sex-specific microbiome composition as an avenue of further research.

Author Summary

Male-bias in clinical testing of drugs has led to a disproportionate number of hepatotoxic events in women. Previous works use gene-by-gene differences in biological sex to explain this discrepancy, but there is little focus on the systematic interactions of these differences. To this end, we use a combination of gene expression data and metabolic modeling to compare metabolic activity between the male and female liver and treated and untreated hepatocytes. We find several subsystems with differential activity in each sex; however, when comparing these subsystems with those pathways altered by hepatotoxic agents, we find little overlap. To explore these differences on a reaction-by-reaction basis, we use the same sex-specific transcriptomic data to contextualize the previously published Human1 human cell metabolic model. In these models we find a difference in flux for the import of bile acids and salts, suggesting a potential difference in enterohepatic circulation. These findings can help guide future drug design, toxicological testing, and sex-specific research to better account for the entire human population.

Identifying metabolic shifts in Crohn's disease using' omics-driven contextualized computational metabolic network models

Crohn's disease (CD) is a chronic inflammatory disease of the gastrointestinal tract. A clear gap in our existing CD diagnostics and current disease management approaches is the lack of highly specific biomarkers that can be used to streamline or personalize disease management. Comprehensive profiling of metabolites holds promise; however, these high-dimensional profiles need to be reduced to have relevance in the context of CD. Machine learning approaches are optimally suited to bridge this gap in knowledge by contextualizing the metabolic alterations in CD using genome-scale metabolic network reconstructions. Our work presents a framework for studying altered metabolic reactions between patients with CD and controls using publicly available transcriptomic data and existing gene-driven metabolic network reconstructions. Additionally, we apply the same methods to patient-derived ileal enteroids to explore the utility of using this experimental in vitro platform for studying CD. Furthermore, we have piloted an untargeted metabolomics approach as a proof-of-concept validation strategy in human ileal mucosal tissue. These findings suggest that in silico metabolic modeling can potentially identify pathways of clinical relevance in CD, paving the way for the future discovery of novel diagnostic biomarkers and therapeutic targets.

Emerging computational paradigms to address the complex role of gut microbial metabolism in cardiovascular diseases

The human gut microbiota and its associated perturbations are implicated in a variety of cardiovascular diseases (CVDs). There is evidence that the structure and metabolic composition of the gut microbiome and some of its metabolites have mechanistic associations with several CVDs. Nevertheless, there is a need to unravel metabolic behavior and underlying mechanisms of microbiome-host interactions. This need is even more highlighted when considering that microbiome-secreted metabolites contributing to CVDs are the subject of intensive research to develop new prevention and therapeutic techniques. In addition to the application of high-throughput data used in microbiome-related studies, advanced computational tools enable us to integrate omics into different mathematical models, including constraint-based models, dynamic models, agent-based models, and machine learning tools, to build a holistic picture of metabolic pathological mechanisms. In this article, we aim to review and introduce state-of-the-art mathematical models and computational approaches addressing the link between the microbiome and CVDs.

Neuraminidases-Key Players in the Inflammatory Response after Pathophysiological Cardiac Stress and Potential New Therapeutic Targets in Cardiac Disease

Glycoproteins and glycolipids on the cell surfaces of vertebrates and higher invertebrates contain α-keto acid sugars called sialic acids, terminally attached to their glycan structures. The actual level of sialylation, regulated through enzymatic removal of the latter ones by NEU enzymes, highly affects protein-protein, cell-matrix and cell-cell interactions. Thus, their regulatory features affect a large number of different cell types, including those of the immune system. Research regarding NEUs within heart and vessels provides new insights of their involvement in the development of cardiovascular pathologies and identifies mechanisms on how inhibiting NEU enzymes can have a beneficial effect on cardiac remodelling and on a number of different cardiac diseases including CMs and atherosclerosis. In this regard, a multitude of clinical studies demonstrated the potential of N-acetylneuraminic acid (Neu5Ac) to serve as a biomarker following cardiac diseases. Anti-influenza drugs i.e., zanamivir and oseltamivir are viral NEU inhibitors, thus, they block the enzymatic activity of NEUs. When considering the improvement in cardiac function in several different cardiac disease animal models, which results from NEU reduction, the inhibition of NEU enzymes provides a new potential therapeutic treatment strategy to treat cardiac inflammatory pathologies, and thus, administrate cardioprotection.

Stable Isotopes for Tracing Cardiac Metabolism in Diseases

Although metabolic remodeling during cardiovascular diseases has been well-recognized for decades, the recent development of analytical platforms and mathematical tools has driven the emergence of assessing cardiac metabolism using tracers. Metabolism is a critical component of cellular functions and adaptation to stress. The pathogenesis of cardiovascular disease involves metabolic adaptation to maintain cardiac contractile function even in advanced disease stages. Stable-isotope tracer measurements are a powerful tool for measuring flux distributions at the whole organism level and assessing metabolic changes at a systems level in vivo. The goal of this review is to summarize techniques and concepts for in vivo or ex vivo stable isotope labeling in cardiovascular research, to highlight mathematical concepts and their limitations, to describe analytical methods at the tissue and single-cell level, and to discuss opportunities to leverage metabolic models to address important mechanistic questions relevant to all patients with cardiovascular disease.