Sulfone-based human liver pyruvate kinase inhibitors - Design, synthesis and in vitro bioactivity

Non-alcoholic fatty liver disease (NAFLD) is a prevalent pathological condition characterised by the accumulation of fat in the liver. Almost one-third of the global population is affected by NAFLD, making it a significant health concern. However, despite its prevalence, there is currently no approved drug specifically designed for the treatment of NAFLD. To address this critical gap, researchers have been investigating potential targets for NAFLD drug development. One promising candidate is the liver isoform of pyruvate kinase (PKL). In recent studies, Urolithin C, an allosteric inhibitor of PKL, has emerged as a potential lead compound for therapeutic intervention. Building upon this knowledge, our team has conducted a comprehensive structure-activity relationship of Urolithin C. In this work, we have employed a scaffold-hopping approach, modifying the urolithin structure by replacing the urolithin carbonyl with a sulfone moiety. Our structure-activity relationship analysis has identified the sulfone group as particularly favourable for potent PKL inhibition. Additionally, we have found that the presence of catechol moieties on the two aromatic rings further improves the inhibitory activity. The most promising inhibitor from this new series displayed nanomolar inhibition, boasting an IC50 value of 0.07 μM. This level of potency rivals that of urolithin D and significantly surpasses the effectiveness of urolithin C by an order of magnitude. To better understand the molecular interactions underlying this inhibition, we obtained the crystal structure of one of the inhibitors complexed with PKL. This structural insight served as a valuable reference point, aiding us in the design of inhibitors.

Sequential sequencing by synthesis and the next-generation sequencing revolution

The impact of next-generation sequencing (NGS) cannot be overestimated. The technology has transformed the field of life science, contributing to a dramatic expansion in our understanding of human health and disease and our understanding of biology and ecology. The vast majority of the major NGS systems today are based on the concept of 'sequencing by synthesis' (SBS) with sequential detection of nucleotide incorporation using an engineered DNA polymerase. Based on this strategy, various alternative platforms have been developed, including the use of either native nucleotides or reversible terminators and different strategies for the attachment of DNA to a solid support. In this review, some of the key concepts leading to this remarkable development are discussed.

Single-cell spatial transcriptome reveals cell-type organization in the macaque cortex

Elucidating the cellular organization of the cerebral cortex is critical for understanding brain structure and function. Using large-scale single-nucleus RNA sequencing and spatial transcriptomic analysis of 143 macaque cortical regions, we obtained a comprehensive atlas of 264 transcriptome-defined cortical cell types and mapped their spatial distribution across the entire cortex. We characterized the cortical layer and region preferences of glutamatergic, GABAergic, and non-neuronal cell types, as well as regional differences in cell-type composition and neighborhood complexity. Notably, we discovered a relationship between the regional distribution of various cell types and the region's hierarchical level in the visual and somatosensory systems. Cross-species comparison of transcriptomic data from human, macaque, and mouse cortices further revealed primate-specific cell types that are enriched in layer 4, with their marker genes expressed in a region-dependent manner. Our data provide a cellular and molecular basis for understanding the evolution, development, aging, and pathogenesis of the primate brain.

Exploration of Novel Urolithin C Derivatives as Non-Competitive Inhibitors of Liver Pyruvate Kinase

The inhibition of liver pyruvate kinase could be beneficial to halt or reverse non-alcoholic fatty liver disease (NAFLD), a progressive accumulation of fat in the liver that can lead eventually to cirrhosis. Recently, urolithin C has been reported as a new scaffold for the development of allosteric inhibitors of liver pyruvate kinase (PKL). In this work, a comprehensive structure-activity analysis of urolithin C was carried out. More than 50 analogues were synthesized and tested regarding the chemical features responsible for the desired activity. These data could pave the way to the development of more potent and selective PKL allosteric inhibitors.

Discovery of drug-omics associations in type 2 diabetes with generative deep-learning models

The application of multiple omics technologies in biomedical cohorts has the potential to reveal patient-level disease characteristics and individualized response to treatment. However, the scale and heterogeneous nature of multi-modal data makes integration and inference a non-trivial task. We developed a deep-learning-based framework, multi-omics variational autoencoders (MOVE), to integrate such data and applied it to a cohort of 789 people with newly diagnosed type 2 diabetes with deep multi-omics phenotyping from the DIRECT consortium. Using in silico perturbations, we identified drug-omics associations across the multi-modal datasets for the 20 most prevalent drugs given to people with type 2 diabetes with substantially higher sensitivity than univariate statistical tests. From these, we among others, identified novel associations between metformin and the gut microbiota as well as opposite molecular responses for the two statins, simvastatin and atorvastatin. We used the associations to quantify drug-drug similarities, assess the degree of polypharmacy and conclude that drug effects are distributed across the multi-omics modalities.

Metabolic modelling of the human gut microbiome in type 2 diabetes patients in response to metformin treatment

The human gut microbiome has been associated with several metabolic disorders including type 2 diabetes mellitus. Understanding metabolic changes in the gut microbiome is important to elucidate the role of gut bacteria in regulating host metabolism. Here, we used available metagenomics data from a metformin study, together with genome-scale metabolic modelling of the key bacteria in individual and community-level to investigate the mechanistic role of the gut microbiome in response to metformin. Individual modelling predicted that species that are increased after metformin treatment have higher growth rates in comparison to species that are decreased after metformin treatment. Gut microbial enrichment analysis showed prior to metformin treatment pathways related to the hypoglycemic effect were enriched. Our observations highlight how the key bacterial species after metformin treatment have commensal and competing behavior, and how their cellular metabolism changes due to different nutritional environment. Integrating different diets showed there were specific microbial alterations between different diets. These results show the importance of the nutritional environment and how dietary guidelines may improve drug efficiency through the gut microbiota.

Serendipitous Identification of a Covalent Activator of Liver Pyruvate Kinase

Enzymes are effective biological catalysts that accelerate almost all metabolic reactions in living organisms. Synthetic modulators of enzymes are useful tools for the study of enzymatic reactions and can provide starting points for the design of new drugs. Here, we report on the discovery of a class of biologically active compounds that covalently modifies lysine residues in human liver pyruvate kinase (PKL), leading to allosteric activation of the enzyme (EC₅₀ =0.29 μM). Surprisingly, the allosteric activation control point resides on the lysine residue K282 present in the catalytic site of PKL. These findings were confirmed by structural data, MS/MS experiments, and molecular modelling studies. Altogether, our study provides a molecular basis for the activation mechanism and establishes a framework for further development of human liver pyruvate kinase covalent activators.

Integrative functional analysis uncovers metabolic differences between Candida species

Candida species are a dominant constituent of the human mycobiome and associated with the development of several diseases. Understanding the Candida species metabolism could provide key insights into their ability to cause pathogenesis. Here, we have developed the BioFung database, providing an efficient annotation of protein-encoding genes. Along, with BioFung, using carbohydrate-active enzyme (CAZymes) analysis, we have uncovered core and accessory features across Candida species demonstrating plasticity, adaption to the environment and acquired features. We show a greater importance of amino acid metabolism, as functional analysis revealed that all Candida species can employ amino acid metabolism. However, metabolomics revealed that only a specific cluster of species (AGAu species—C. albicans, C. glabrata and C. auris) utilised amino acid metabolism including arginine, cysteine, and methionine metabolism potentially improving their competitive fitness in pathogenesis. We further identified critical metabolic pathways in the AGAu cluster with biomarkers and anti-fungal target potential in the CAZyme profile, polyamine, choline and fatty acid biosynthesis pathways. This study, combining genomic analysis, and validation with gene expression and metabolomics, highlights the metabolic diversity with AGAu species that underlies their remarkable ability to dominate they mycobiome and cause disease.

Metabolic differences between Candida species are uncovered using the BioFung database alongside genomic and metabolic analysis.

Genome-scale metabolic modelling of the human gut microbiome reveals changes of the glyoxylate and dicarboxylate metabolism in metabolic disorders

The human gut microbiome has been associated with metabolic disorders including obesity, type 2 diabetes, and atherosclerosis. Understanding the contribution of microbiome metabolic changes is important for elucidating the role of gut bacteria in regulating metabolism. We used available metagenomics data from these metabolic disorders, together with genome-scale metabolic modeling of key bacteria in the individual and community-level to investigate the mechanistic role of the gut microbiome in metabolic diseases. Modeling predicted increased levels of glutamate consumption along with the production of ammonia, arginine, and proline in gut bacteria common across the disorders. Abundance profiles and network-dependent analysis identified the enrichment of tartrate dehydrogenase in the disorders. Moreover, independent plasma metabolite levels showed associations between metabolites including proline and tyrosine and an increased tartrate metabolism in healthy obese individuals. We, therefore, propose that an increased tartrate metabolism could be a significant mediator of the microbiome metabolic changes in metabolic disorders.

•

Metagenomic analysis highlights key common bacterial species across metabolic diseases

•

Metabolic models showed higher levels of acetate produced by disease enriched bacteria

•

Reaction analysis revealed increases in the glyoxylate and dicarboxylate pathway

•

Metabolomics and modeling analysis showed the potential role of tartrate metabolism

Cell transcriptomic atlas of the non-human primate Macaca fascicularis

Studying tissue composition and function in non-human primates (NHPs) is crucial to understand the nature of our own species. Here we present a large-scale cell transcriptomic atlas that encompasses over 1 million cells from 45 tissues of the adult NHP Macaca fascicularis. This dataset provides a vast annotated resource to study a species phylogenetically close to humans. To demonstrate the utility of the atlas, we have reconstructed the cell-cell interaction networks that drive Wnt signalling across the body, mapped the distribution of receptors and co-receptors for viruses causing human infectious diseases, and intersected our data with human genetic disease orthologues to establish potential clinical associations. Our M. fascicularis cell atlas constitutes an essential reference for future studies in humans and NHPs.

Prediction of drug candidates for clear cell renal cell carcinoma using a systems biology-based drug repositioning approach

BACKGROUND

The response rates of the clinical chemotherapies are still low in clear cell renal cell carcinoma (ccRCC). Computational drug repositioning is a promising strategy to discover new uses for existing drugs to treat patients who cannot get benefits from clinical drugs.

METHODS

We proposed a systematic approach which included the target prediction based on the co-expression network analysis of transcriptomics profiles of ccRCC patients and drug repositioning for cancer treatment based on the analysis of shRNA- and drug-perturbed signature profiles of human kidney cell line.

FINDINGS

First, based on the gene co-expression network analysis, we identified two types of gene modules in ccRCC, which significantly enriched with unfavorable and favorable signatures indicating poor and good survival outcomes of patients, respectively. Then, we selected four genes, BUB1B, RRM2, ASF1B and CCNB2, as the potential drug targets based on the topology analysis of modules. Further, we repurposed three most effective drugs for each target by applying the proposed drug repositioning approach. Finally, we evaluated the effects of repurposed drugs using an in vitro model and observed that these drugs inhibited the protein levels of their corresponding target genes and cell viability.

INTERPRETATION

These findings proved the usefulness and efficiency of our approach to improve the drug repositioning researches for cancer treatment and precision medicine.

FUNDING

This study was funded by Knut and Alice Wallenberg Foundation and Bash Biotech Inc., San Diego, CA, USA.

A Gene Co-Expression Network-Based Drug Repositioning Approach Identifies Candidates for Treatment of Hepatocellular Carcinoma

Simple Summary

Hepatocellular carcinoma (HCC) is the most common malignancy of liver cancer. However, treatment of HCC is still severely limited due to limitation of drug therapy. We aimed to screen more possible target genes and candidate drugs for HCC, exploring the possibility of drug treatments from systems biological perspective. We identified ten candidate target genes, which are hub genes in HCC co-expression networks, which also possess significant prognostic value in two independent HCC cohorts. The rationality of these target genes was well demonstrated through variety analyses of patient expression profiles. We then screened candidate drugs for target genes and finally identified withaferin-a and mitoxantrone as the candidate drug for HCC treatment. The drug effectiveness was validated in in vitro model and computational analysis, providing more evidence for our drug repositioning method and results.

Abstract

Hepatocellular carcinoma (HCC) is a malignant liver cancer that continues to increase deaths worldwide owing to limited therapies and treatments. Computational drug repurposing is a promising strategy to discover potential indications of existing drugs. In this study, we present a systematic drug repositioning method based on comprehensive integration of molecular signatures in liver cancer tissue and cell lines. First, we identify robust prognostic genes and two gene co-expression modules enriched in unfavorable prognostic genes based on two independent HCC cohorts, which showed great consistency in functional and network topology. Then, we screen 10 genes as potential target genes for HCC on the bias of network topology analysis in these two modules. Further, we perform a drug repositioning method by integrating the shRNA and drug perturbation of liver cancer cell lines and identifying potential drugs for every target gene. Finally, we evaluate the effects of the candidate drugs through an in vitro model and observe that two identified drugs inhibited the protein levels of their corresponding target genes and cell migration, also showing great binding affinity in protein docking analysis. Our study demonstrates the usefulness and efficiency of network-based drug repositioning approach to discover potential drugs for cancer treatment and precision medicine approach.

Genome-wide annotation of protein-coding genes in pig

BACKGROUND

There is a need for functional genome-wide annotation of the protein-coding genes to get a deeper understanding of mammalian biology. Here, a new annotation strategy is introduced based on dimensionality reduction and density-based clustering of whole-body co-expression patterns. This strategy has been used to explore the gene expression landscape in pig, and we present a whole-body map of all protein-coding genes in all major pig tissues and organs.

RESULTS

An open-access pig expression map ( www.rnaatlas.org ) is presented based on the expression of 350 samples across 98 well-defined pig tissues divided into 44 tissue groups. A new UMAP-based classification scheme is introduced, in which all protein-coding genes are stratified into tissue expression clusters based on body-wide expression profiles. The distribution and tissue specificity of all 22,342 protein-coding pig genes are presented.

CONCLUSIONS

Here, we present a new genome-wide annotation strategy based on dimensionality reduction and density-based clustering. A genome-wide resource of the transcriptome map across all major tissues and organs in pig is presented, and the data is available as an open-access resource ( www.rnaatlas.org ), including a comparison to the expression of human orthologs.

Host-mycobiome metabolic interactions in health and disease

Fungal communities (mycobiome) have an important role in sustaining the resilience of complex microbial communities and maintenance of homeostasis. The mycobiome remains relatively unexplored compared to the bacteriome despite increasing evidence highlighting their contribution to host-microbiome interactions in health and disease. Despite being a small proportion of the total species, fungi constitute a large proportion of the biomass within the human microbiome and thus serve as a potential target for metabolic reprogramming in pathogenesis and disease mechanism. Metabolites produced by fungi shape host niches, induce immune tolerance and changes in their levels prelude changes associated with metabolic diseases and cancer. Given the complexity of microbial interactions, studying the metabolic interplay of the mycobiome with both host and microbiome is a demanding but crucial task. However, genome-scale modelling and synthetic biology can provide an integrative platform that allows elucidation of the multifaceted interactions between mycobiome, microbiome and host. The inferences gained from understanding mycobiome interplay with other organisms can delineate the key role of the mycobiome in pathophysiology and reveal its role in human disease.

Spatiotemporal characterization of cellular tau pathology in the human locus coeruleus-pericoerulear complex by three-dimensional imaging

Tau pathology of the noradrenergic locus coeruleus (LC) is a hallmark of several age-related neurodegenerative disorders, including Alzheimer's disease. However, a comprehensive neuropathological examination of the LC is difficult due to its small size and rod-like shape. To investigate the LC cytoarchitecture and tau cytoskeletal pathology in relation to possible propagation patterns of disease-associated tau in an unprecedented large-scale three-dimensional view, we utilized volume immunostaining and optical clearing technology combined with light sheet fluorescence microscopy. We examined AT8⁺ pathological tau in the LC/pericoerulear region of 20 brains from Braak neurofibrillary tangle (NFT) stage 0-6. We demonstrate an intriguing morphological complexity and heterogeneity of AT8⁺ cellular structures in the LC, representing various intracellular stages of NFT maturation and their diverse transition forms. We describe novel morphologies of neuronal tau pathology such as AT8⁺ cells with fine filamentous somatic protrusions or with disintegrating soma. We show that gradual dendritic atrophy is the first morphological sign of the degeneration of tangle-bearing neurons, even preceding axonal lesions. Interestingly, irrespective of the Braak NFT stage, tau pathology is more advanced in the dorsal LC that preferentially projects to vulnerable forebrain regions in Alzheimer's disease, like the hippocampus or neocortical areas, compared to the ventral LC projecting to the cerebellum and medulla. Moreover, already in the precortical Braak 0 stage, 3D analysis reveals clustering tendency and dendro-dendritic close appositions of AT8⁺ LC neurons, AT8⁺ long axons of NFT-bearing cells that join the ascending dorsal noradrenergic bundle after leaving the LC, as well as AT8⁺ processes of NFT-bearing LC neurons that target the 4th ventricle wall. Our study suggests that the unique cytoarchitecture, comprised of a densely packed and dendritically extensively interconnected neuronal network with long projections, makes the human LC to be an ideal anatomical template for early accumulation and trans-neuronal spreading of hyperphosphorylated tau.

Plasma proteins elevated in severe asthma despite oral steroid use and unrelated to Type-2 inflammation

Rationale

Asthma phenotyping requires novel biomarker discovery.

Objectives

To identify plasma biomarkers associated with asthma phenotypes by application of a new proteomic panel to samples from two well-characterised cohorts of severe (SA) and mild-to-moderate (MMA) asthmatics, COPD subjects and healthy controls (HCs).

Methods

An antibody-based array targeting 177 proteins predominantly involved in pathways relevant to inflammation, lipid metabolism, signal transduction and extracellular matrix was applied to plasma from 525 asthmatics and HCs in the U-BIOPRED cohort, and 142 subjects with asthma and COPD from the validation cohort BIOAIR. Effects of oral corticosteroids (OCS) were determined by a 2-week, placebo-controlled OCS trial in BIOAIR, and confirmed by relation to objective OCS measures in U-BIOPRED.

Results

In U-BIOPRED, 110 proteins were significantly different, mostly elevated, in SA compared to MMA and HCs. 10 proteins were elevated in SA versus MMA in both U-BIOPRED and BIOAIR (alpha-1-antichymotrypsin, apolipoprotein-E, complement component 9, complement factor I, macrophage inflammatory protein-3, interleukin-6, sphingomyelin phosphodiesterase 3, TNF receptor superfamily member 11a, transforming growth factor-β and glutathione S-transferase). OCS treatment decreased most proteins, yet differences between SA and MMA remained following correction for OCS use. Consensus clustering of U-BIOPRED protein data yielded six clusters associated with asthma control, quality of life, blood neutrophils, high-sensitivity C-reactive protein and body mass index, but not Type-2 inflammatory biomarkers. The mast cell specific enzyme carboxypeptidase A3 was one major contributor to cluster differentiation.

Conclusions

The plasma proteomic panel revealed previously unexplored yet potentially useful Type-2-independent biomarkers and validated several proteins with established involvement in the pathophysiology of SA.

Application of new proteomic panel in two established European asthma cohorts identifies plasma proteins associated with disease severity independently of Type-2 inflammation, suggesting potentially useful novel biomarkers and therapeutic targets.https://bit.ly/3jtTq5m

Multi-omics approaches for revealing the complexity of cardiovascular disease

The development and progression of cardiovascular disease (CVD) can mainly be attributed to the narrowing of blood vessels caused by atherosclerosis and thrombosis, which induces organ damage that will result in end-organ dysfunction characterized by events such as myocardial infarction or stroke. It is also essential to consider other contributory factors to CVD, including cardiac remodelling caused by cardiomyopathies and co-morbidities with other diseases such as chronic kidney disease. Besides, there is a growing amount of evidence linking the gut microbiota to CVD through several metabolic pathways. Hence, it is of utmost importance to decipher the underlying molecular mechanisms associated with these disease states to elucidate the development and progression of CVD. A wide array of systems biology approaches incorporating multi-omics data have emerged as an invaluable tool in establishing alterations in specific cell types and identifying modifications in signalling events that promote disease development. Here, we review recent studies that apply multi-omics approaches to further understand the underlying causes of CVD and provide possible treatment strategies by identifying novel drug targets and biomarkers. We also discuss very recent advances in gut microbiota research with an emphasis on how diet and microbial composition can impact the development of CVD. Finally, we present various biological network analyses and other independent studies that have been employed for providing mechanistic explanation and developing treatment strategies for end-stage CVD, namely myocardial infarction and stroke.

iNetModels 2.0: an interactive visualization and database of multi-omics data

It is essential to reveal the associations between various omics data for a comprehensive understanding of the altered biological process in human wellness and disease. To date, very few studies have focused on collecting and exhibiting multi-omics associations in a single database. Here, we present iNetModels, an interactive database and visualization platform of Multi-Omics Biological Networks (MOBNs). This platform describes the associations between the clinical chemistry, anthropometric parameters, plasma proteomics, plasma metabolomics, as well as metagenomics for oral and gut microbiome obtained from the same individuals. Moreover, iNetModels includes tissue- and cancer-specific Gene Co-expression Networks (GCNs) for exploring the connections between the specific genes. This platform allows the user to interactively explore a single feature's association with other omics data and customize its particular context (e.g. male/female specific). The users can also register their data for sharing and visualization of the MOBNs and GCNs. Moreover, iNetModels allows users who do not have a bioinformatics background to facilitate human wellness and disease research. iNetModels can be accessed freely at https://inetmodels.com without any limitation.

iNetModels 2.0: an interactive visualization and database of multi-omics data

It is essential to reveal the associations between various omics data for a comprehensive understanding of the altered biological process in human wellness and disease. To date, very few studies have focused on collecting and exhibiting multi-omics associations in a single database. Here, we present iNetModels, an interactive database and visualization platform of Multi-Omics Biological Networks (MOBNs). This platform describes the associations between the clinical chemistry, anthropometric parameters, plasma proteomics, plasma metabolomics, as well as metagenomics for oral and gut microbiome obtained from the same individuals. Moreover, iNetModels includes tissue- and cancer-specific Gene Co-expression Networks (GCNs) for exploring the connections between the specific genes. This platform allows the user to interactively explore a single feature's association with other omics data and customize its particular context (e.g. male/female specific). The users can also register their data for sharing and visualization of the MOBNs and GCNs. Moreover, iNetModels allows users who do not have a bioinformatics background to facilitate human wellness and disease research. iNetModels can be accessed freely at https://inetmodels.com without any limitation.

Stratification of patients with clear cell renal cell carcinoma to facilitate drug repositioning

Clear cell renal cell carcinoma (ccRCC) is the most common histological type of kidney cancer and has high heterogeneity. Stratification of ccRCC is important since distinct subtypes differ in prognosis and treatment. Here, we applied a systems biology approach to stratify ccRCC into three molecular subtypes with different mRNA expression patterns and prognosis of patients. Further, we developed a set of biomarkers that could robustly classify the patients into each of the three subtypes and predict the prognosis of patients. Then, we reconstructed subtype-specific metabolic models and performed essential gene analysis to identify the potential drug targets. We identified four drug targets, including SOAT1, CRLS1, and ACACB, essential in all the three subtypes and GPD2, exclusively essential to subtype 1. Finally, we repositioned mitotane, an FDA-approved SOAT1 inhibitor, to treat ccRCC and showed that it decreased tumor cell viability and inhibited tumor cell growth based on in vitro experiments.

•

Three consistent molecular ccRCC subtypes were found to guide patients' prognoses

•

REOs-based biomarker was developed to robustly classify patients at individual level

•

SOAT1 is identified as a common drug target for all ccRCC subtypes

•

Mitotane was repositioned treatment of ccRCC via inhibiting SOAT1

Informing Pharmacokinetic Models With Physiological Data: Oral Population Modeling of L-Serine in Humans

To determine how to set optimal oral L-serine (serine) dose levels for a clinical trial, existing literature was surveyed. Data sufficient to set the dose was inadequate, and so an (n = 10) phase I-A calibration trial was performed, administering serine with and without other oral agents. We analyzed the trial and the literature data using pharmacokinetic (PK) modeling and statistical analysis. The therapeutic goal is to modulate specific serine-related metabolic pathways in the liver using the lowest possible dose which gives the desired effect since the upper bound was expected to be limited by toxicity. A standard PK approach, in which a common model structure was selected using a fit to data, yielded a model with a single central compartment corresponding to plasma, clearance from that compartment, and an endogenous source of serine. To improve conditioning, a parametric structure was changed to estimate ratios (bioavailability over volume, for example). Model fit quality was improved and the uncertainty in estimated parameters was reduced. Because of the particular interest in the fate of serine, the model was used to estimate whether serine is consumed in the gut, absorbed by the liver, or entered the blood in either a free state, or in a protein- or tissue-bound state that is not measured by our assay. The PK model structure was set up to represent relevant physiology, and this quantitative systems biology approach allowed a broader set of physiological data to be used to narrow parameter and prediction confidence intervals, and to better understand the biological meaning of the data. The model results allowed us to determine the optimal human dose for future trials, including a trial design component including IV and tracer studies. A key contribution is that we were able to use human physiological data from the literature to inform the PK model and to set reasonable bounds on parameters, and to improve model conditioning. Leveraging literature data produced a more predictive, useful model.

Human Immune System Variation during 1 Year

The human immune system varies extensively between individuals, but variation within individuals over time has not been well characterized. Systems-level analyses allow for simultaneous quantification of many interacting immune system components and the inference of global regulatory principles. Here, we present a longitudinal, systems-level analysis in 99 healthy adults 50 to 65 years of age and sampled every third month for 1 year. We describe the structure of interindividual variation and characterize extreme phenotypes along a principal curve. From coordinated measurement fluctuations, we infer relationships between 115 immune cell populations and 750 plasma proteins constituting the blood immune system. While most individuals have stable immune systems, the degree of longitudinal variability is an individual feature. The most variable individuals, in the absence of overt infections, exhibited differences in markers of metabolic health suggestive of a possible link between metabolic and immunologic homeostatic regulation.

Systematic analysis of gut microbiome reveals the role of bacterial folate and homocysteine metabolism in Parkinson's disease

Parkinson's disease (PD) is the most common progressive neurological disorder compromising motor functions. However, nonmotor symptoms, such as gastrointestinal (GI) dysfunction, precede those affecting movement. Evidence of an early involvement of the GI tract and enteric nervous system highlights the need for better understanding of the role of gut microbiota in GI complications in PD. Here, we investigate the gut microbiome of patients with PD using metagenomics and serum metabolomics. We integrate these data using metabolic modeling and construct an integrative correlation network giving insight into key microbial species linked with disease severity, GI dysfunction, and age of patients with PD. Functional analysis reveals an increased microbial capability to degrade mucin and host glycans in PD. Personalized community-level metabolic modeling reveals the microbial contribution to folate deficiency and hyperhomocysteinemia observed in patients with PD. The metabolic modeling approach could be applied to uncover gut microbial metabolic contributions to PD pathophysiology.

Applications of Genome-Wide Screening and Systems Biology Approaches in Drug Repositioning

Modern drug discovery through de novo drug discovery entails high financial costs, low success rates, and lengthy trial periods. Drug repositioning presents a suitable approach for overcoming these issues by re-evaluating biological targets and modes of action of approved drugs. Coupling high-throughput technologies with genome-wide essentiality screens, network analysis, genome-scale metabolic modeling, and machine learning techniques enables the proposal of new drug-target signatures and uncovers unanticipated modes of action for available drugs. Here, we discuss the current issues associated with drug repositioning in light of curated high-throughput multi-omic databases, genome-wide screening technologies, and their application in systems biology/medicine approaches.

Integration of molecular profiles in a longitudinal wellness profiling cohort

An important aspect of precision medicine is to probe the stability in molecular profiles among healthy individuals over time. Here, we sample a longitudinal wellness cohort with 100 healthy individuals and analyze blood molecular profiles including proteomics, transcriptomics, lipidomics, metabolomics, autoantibodies and immune cell profiling, complemented with gut microbiota composition and routine clinical chemistry. Overall, our results show high variation between individuals across different molecular readouts, while the intra-individual baseline variation is low. The analyses show that each individual has a unique and stable plasma protein profile throughout the study period and that many individuals also show distinct profiles with regards to the other omics datasets, with strong underlying connections between the blood proteome and the clinical chemistry parameters. In conclusion, the results support an individual-based definition of health and show that comprehensive omics profiling in a longitudinal manner is a path forward for precision medicine.

An important aspect of precision medicine is to probe the stability in molecular profiles among healthy individuals over time. Here, the authors sample a longitudinal wellness cohort and analyse blood molecular profiles as well as gut microbiota composition.

Current Status of COVID-19 Therapies and Drug Repositioning Applications

The rapid and global spread of a new human coronavirus (SARS-CoV-2) has produced an immediate urgency to discover promising targets for the treatment of COVID-19. Drug repositioning is an attractive approach that can facilitate the drug discovery process by repurposing existing pharmaceuticals to treat illnesses other than their primary indications. Here, we review current information concerning the global health issue of COVID-19 including promising approved drugs and ongoing clinical trials for prospective treatment options. In addition, we describe computational approaches to be used in drug repurposing and highlight examples of in silico studies of drug development efforts against SARS-CoV-2.

Systems Biology Approaches to Understand the Host-Microbiome Interactions in Neurodegenerative Diseases

Neurodegenerative diseases (NDDs) comprise a broad range of progressive neurological disorders with multifactorial etiology contributing to disease pathophysiology. Evidence of the microbiome involvement in the gut-brain axis urges the interest in understanding metabolic interactions between the microbiota and host physiology in NDDs. Systems Biology offers a holistic integrative approach to study the interplay between the different biologic systems as part of a whole, and may elucidate the host–microbiome interactions in NDDs. We reviewed direct and indirect pathways through which the microbiota can modulate the bidirectional communication of the gut-brain axis, and explored the evidence of microbial dysbiosis in Alzheimer’s and Parkinson’s diseases. As the gut microbiota being strongly affected by diet, the potential approaches to targeting the human microbiota through diet for the stimulation of neuroprotective microbial-metabolites secretion were described. We explored the potential of Genome-scale metabolic models (GEMs) to infer microbe-microbe and host-microbe interactions and to identify the microbiome contribution to disease development or prevention. Finally, a systemic approach based on GEMs and ‘omics integration, that would allow the design of sustainable personalized anti-inflammatory diets in NDDs prevention, through the modulation of gut microbiota was described.

A systems biology approach for studying neurodegenerative diseases

Neurodegenerative diseases (NDDs), such as Alzheimer's (AD) and Parkinson's (PD), are among the leading causes of lost years of healthy life and exert a great strain on public healthcare systems. Despite being first described more than a century ago, no effective cure exists for AD or PD. Although extensively characterised at the molecular level, traditional neurodegeneration research remains marred by narrow-sense approaches surrounding amyloid β (Aβ), tau, and α-synuclein (α-syn). A systems biology approach enables the integration of multi-omics data and informs discovery of biomarkers, drug targets, and treatment strategies. Here, we present a comprehensive timeline of high-throughput data collection, and associated biotechnological advancements and computational analysis related to AD and PD. We hereby propose that a philosophical change in the definitions of AD and PD is now needed.

Improvement in the Current Therapies for Hepatocellular Carcinoma Using a Systems Medicine Approach

Hepatocellular carcinoma (HCC) is one of the leading causes of cancer-related death primarily due to the lack of effective targeted therapies. Despite the distinct morphological and phenotypic patterns of HCC, treatment strategies are restricted to relatively homogeneous therapies, including multitargeted tyrosine kinase inhibitors and immune checkpoint inhibitors. Therefore, more effective therapy options are needed to target dysregulated metabolic and molecular pathways in HCC. Integrative genomic profiling of HCC patients provides insight into the most frequently mutated genes and molecular targets, including telomerase reverse transcriptase, the TP53 gene, and the Wnt/β-catenin signaling pathway oncogene (CTNNB1). Moreover, emerging techniques, such as genome-scale metabolic models may elucidate the underlying cancer-specific metabolism, which allows for the discovery of potential drug targets and identification of biomarkers. De novo lipogenesis has been revealed as consistently upregulated since it is required for cell proliferation in all HCC patients. The metabolic network-driven stratification of HCC patients in terms of redox responses, utilization of metabolites, and subtype-specific pathways may have clinical implications to drive the development of personalized medicine. In this review, the current and emerging therapeutic targets in light of molecular approaches and metabolic network-based strategies are summarized, prompting effective treatment of HCC patients.

Addressing the heterogeneity in liver diseases using biological networks

The abnormalities in human metabolism have been implicated in the progression of several complex human diseases, including certain cancers. Hence, deciphering the underlying molecular mechanisms associated with metabolic reprogramming in a disease state can greatly assist in elucidating the disease aetiology. An invaluable tool for establishing connections between global metabolic reprogramming and disease development is the genome-scale metabolic model (GEM). Here, we review recent work on the reconstruction of cell/tissue-type and cancer-specific GEMs and their use in identifying metabolic changes occurring in response to liver disease development, stratification of the heterogeneous disease population and discovery of novel drug targets and biomarkers. We also discuss how GEMs can be integrated with other biological networks for generating more comprehensive cell/tissue models. In addition, we review the various biological network analyses that have been employed for the development of efficient treatment strategies. Finally, we present three case studies in which independent studies converged on conclusions underlying liver disease.

Classification of clear cell renal cell carcinoma based on PKM alternative splicing

Clear cell renal cell carcinoma (ccRCC) accounts for 70-80% of kidney cancer diagnoses and displays high molecular and histologic heterogeneity. Hence, it is necessary to reveal the underlying molecular mechanisms involved in progression of ccRCC to better stratify the patients and design effective treatment strategies. Here, we analyzed the survival outcome of ccRCC patients as a consequence of the differential expression of four transcript isoforms of the pyruvate kinase muscle type (PKM). We first extracted a classification biomarker consisting of eight gene pairs whose within-sample relative expression orderings (REOs) could be used to robustly classify the patients into two groups with distinct molecular characteristics and survival outcomes. Next, we validated our findings in a validation cohort and an independent Japanese ccRCC cohort. We finally performed drug repositioning analysis based on transcriptomic expression profiles of drug-perturbed cancer cell lines and proposed that paracetamol, nizatidine, dimethadione and conessine can be repurposed to treat the patients in one of the subtype of ccRCC whereas chenodeoxycholic acid, fenoterol and hexylcaine can be repurposed to treat the patients in the other subtype.

A genome-wide transcriptomic analysis of protein-coding genes in human blood cells

Blood is the predominant source for molecular analyses in humans, both in clinical and research settings. It is the target for many therapeutic strategies, emphasizing the need for comprehensive molecular maps of the cells constituting human blood. In this study, we performed a genome-wide transcriptomic analysis of protein-coding genes in sorted blood immune cell populations to characterize the expression levels of each individual gene across the blood cell types. All data are presented in an interactive, open-access Blood Atlas as part of the Human Protein Atlas and are integrated with expression profiles across all major tissues to provide spatial classification of all protein-coding genes. This allows for a genome-wide exploration of the expression profiles across human immune cell populations and all major human tissues and organs.

Systems biology perspective for studying the gut microbiota in human physiology and liver diseases

The advancement in high-throughput sequencing technologies and systems biology approaches have revolutionized our understanding of biological systems and opened a new path to investigate unacknowledged biological phenomena. In parallel, the field of human microbiome research has greatly evolved and the relative contribution of the gut microbiome to health and disease have been systematically explored. This review provides an overview of the network-based and translational systems biology-based studies focusing on the function and composition of gut microbiota. We also discussed the association between the gut microbiome and the overall human physiology, as well as hepatic diseases and other metabolic disorders.

Systems biology based drug repositioning for development of cancer therapy

Drug repositioning is a powerful method that can assists the conventional drug discovery process by using existing drugs for treatment of a disease rather than its original indication. The first examples of repurposed drugs were discovered serendipitously, however data accumulated by high-throughput screenings and advancements in computational biology methods have paved the way for rational drug repositioning methods. As chemotherapeutic agents have notorious side effects that significantly reduce quality of life, drug repositioning promises repurposed noncancer drugs with little or tolerable adverse effects for cancer patients. Here, we review current drug-related data types and databases including some examples of web-based drug repositioning tools. Next, we describe systems biology approaches to be used in drug repositioning for effective cancer therapy. Finally, we highlight examples of mostly repurposed drugs for cancer treatment and provide an overview of future expectations in the field for development of effective treatment strategies.

TCSBN: a database of tissue and cancer specific biological networks

Biological networks provide new opportunities for understanding the cellular biology in both health and disease states. We generated tissue specific integrated networks (INs) for liver, muscle and adipose tissues by integrating metabolic, regulatory and protein-protein interaction networks. We also generated human co-expression networks (CNs) for 46 normal tissues and 17 cancers to explore the functional relationships between genes as well as their relationships with biological functions, and investigate the overlap between functional and physical interactions provided by CNs and INs, respectively. These networks can be employed in the analysis of omics data, provide detailed insight into disease mechanisms by identifying the key biological components and eventually can be used in the development of efficient treatment strategies. Moreover, comparative analysis of the networks may allow for the identification of tissue-specific targets that can be used in the development of drugs with the minimum toxic effect to other human tissues. These context-specific INs and CNs are presented in an interactive website http://inetmodels.com without any limitation.

Broad Views of Non-alcoholic Fatty Liver Disease

Multi-omics multi-tissue data are used to interpret genome-wide association study results from mice to identify key driver genes of non-alcoholic fatty liver disease. Non-alcoholic fatty liver disease (NAFLD) is the accumulation of fat (steatosis) in the liver due to causes other than excessive alcohol consumption. The disease may progress to more severe forms of liver diseases, including non-alcoholic steatohepatitis, cirrhosis, and hepatocellular carcinoma. In this issue of Cell Systems, Krishnan et al. (2018) reveal mechanisms underlying NAFLD by generating multi-omics data using liver and adipose tissues obtained from the Hybrid Mouse Diversity Panel, consisting of 113 mouse strains with various degrees of NAFLD. The study identified key driver genes of NAFLD that can be used in the development of efficient treatment strategies and illustrates the potential utility of systematic analysis of multi-layer biological networks.

Screening a Resource of Recombinant Protein Fragments for Targeted Proteomics

The availability of proteomics resources hosting protein and peptide standards, as well as the data describing their analytical performances, will continue to enhance our current capabilities to develop targeted proteomics methods for quantitative biology. This study describes the analysis of a resource of 26,840 individually purified recombinant protein fragments corresponding to more than 16,000 human protein-coding genes. The resource was screened to identify proteotypic peptides suitable for targeted proteomics efforts, and we report LC-MS/MS assay coordinates for more than 25,000 proteotypic peptides, corresponding to more than 10,000 unique proteins. Additionally, peptide formation and digestion kinetics were, for a subset of the standards, monitored using a time-course protocol involving parallel digestion of isotope-labeled recombinant protein standards and endogenous human plasma proteins. We show that the strategy by adding isotope-labeled recombinant proteins before trypsin digestion enables short digestion protocols (≤60 min) with robust quantitative precision. In a proof-of-concept study, we quantified 23 proteins in human plasma using assay parameters defined in our study and used the standards to describe distinct clusters of individuals linked to different levels of LPA, APOE, SERPINA5, and TFRC. In summary, we describe the use and utility of a resource of recombinant proteins to identify proteotypic peptides useful for targeted proteomics assay development.

Multi-Omic Data Interpretation to Repurpose Subtype Specific Drug Candidates for Breast Cancer

Triple-negative breast cancer (TNBC), which is largely synonymous with the basal-like molecular subtype, is the 5th leading cause of cancer deaths for women in the United States. The overall prognosis for TNBC patients remains poor given that few treatment options exist; including targeted therapies (not FDA approved), and multi-agent chemotherapy as standard-of-care treatment. TNBC like other complex diseases is governed by the perturbations of the complex interaction networks thereby elucidating the underlying molecular mechanisms of this disease in the context of network principles, which have the potential to identify targets for drug development. Here, we present an integrated "omics" approach based on the use of transcriptome and interactome data to identify dynamic/active protein-protein interaction networks (PPINs) in TNBC patients. We have identified three highly connected modules, EED, DHX9, and AURKA, which are extremely activated in TNBC tumors compared to both normal tissues and other breast cancer subtypes. Based on the functional analyses, we propose that these modules are potential drivers of proliferation and, as such, should be considered candidate molecular targets for drug development or drug repositioning in TNBC. Consistent with this argument, we repurposed steroids, anti-inflammatory agents, anti-infective agents, cardiovascular agents for patients with basal-like breast cancer. Finally, we have performed essential metabolite analysis on personalized genome-scale metabolic models and found that metabolites such as sphingosine-1-phosphate and cholesterol-sulfate have utmost importance in TNBC tumor growth.

Discovery of therapeutic agents for prostate cancer using genome-scale metabolic modeling and drug repositioning

BACKGROUND

Genome-scale metabolic models (GEMs) offer insights into cancer metabolism and have been used to identify potential biomarkers and drug targets. Drug repositioning is a time- and cost-effective method of drug discovery that can be applied together with GEMs for effective cancer treatment.

METHODS

In this study, we reconstruct a prostate cancer (PRAD)-specific GEM for exploring prostate cancer metabolism and also repurposing new therapeutic agents that can be used in development of effective cancer treatment. We integrate global gene expression profiling of cell lines with >1000 different drugs through the use of prostate cancer GEM and predict possible drug-gene interactions.

FINDINGS

We identify the key reactions with altered fluxes based on the gene expression changes and predict the potential drug effect in prostate cancer treatment. We find that sulfamethoxypyridazine, azlocillin, hydroflumethiazide, and ifenprodil can be repurposed for the treatment of prostate cancer based on an in silico cell viability assay. Finally, we validate the effect of ifenprodil using an in vitro cell assay and show its inhibitory effect on a prostate cancer cell line.

INTERPRETATION

Our approach demonstate how GEMs can be used to predict therapeutic agents for cancer treatment based on drug repositioning. Besides, it paved a way and shed a light on the applicability of computational models to real-world biomedical or pharmaceutical problems.

A deep proteome and transcriptome abundance atlas of 29 healthy human tissues

Genome-, transcriptome- and proteome-wide measurements provide insights into how biological systems are regulated. However, fundamental aspects relating to which human proteins exist, where they are expressed and in which quantities are not fully understood. Therefore, we generated a quantitative proteome and transcriptome abundance atlas of 29 paired healthy human tissues from the Human Protein Atlas project representing human genes by 18,072 transcripts and 13,640 proteins including 37 without prior protein-level evidence. The analysis revealed that hundreds of proteins, particularly in testis, could not be detected even for highly expressed mRNAs, that few proteins show tissue-specific expression, that strong differences between mRNA and protein quantities within and across tissues exist and that protein expression is often more stable across tissues than that of transcripts. Only 238 of 9,848 amino acid variants found by exome sequencing could be confidently detected at the protein level showing that proteogenomics remains challenging, needs better computational methods and requires rigorous validation. Many uses of this resource can be envisaged including the study of gene/protein expression regulation and biomarker specificity evaluation.

Discovery of KIRREL as a biomarker for prognostic stratification of patients with thin melanoma

There is a great unmet clinical need to identify patients with thin primary cutaneous melanomas (T1, Breslow thickness ≤ 1 mm) who have a high risk for tumour recurrence and death from melanoma. Kin of IRRE-like protein 1 (KIRREL/NEPH1) is expressed in podocytes and involved in glomerular filtration. Screening in the Human Protein Atlas portal revealed a particularly high expression of KIRREL in melanoma, both at the mRNA and protein levels. In this study, we followed up on these findings and examined the prognostic value of KIRREL in a population-based cohort. Immunohistochemical expression of KIRREL was examined in tissue microarrays with a subset of primary tumours and paired lymph node metastases from an original cohort of 268 incident cases of melanoma in the Malmö Diet and Cancer study. KIRREL mRNA expression was examined in 103 melanoma cases in The Cancer Genome Atlas (TCGA). Membranous/cytoplasmic expression of KIRREL was detected in 158/185 (85.4%) primary tumours and 18/19 (94.7%) metastases. High expression of KIRREL was significantly associated with several unfavourable clinicopathological factors. High KIRREL protein expression was an independent factor of reduced recurrence free and melanoma specific survival, particularly in thin melanomas, even outperforming absolute thickness and ulceration (HR = 30.85; 95% CI 1.54-616.36 and HR = 6.32 95% CI 1.19-33.65). High mRNA levels of KIRREL were not significantly associated with survival in TCGA. In conclusion, KIRREL is not only a novel potential diagnostic marker for melanoma, but may also be a useful prognostic biomarker for improved stratification of patients with thin melanoma. These findings may be of high clinical relevance and therefore merit further validation.

Pyruvate kinase L/R is a regulator of lipid metabolism and mitochondrial function

The pathogenesis of non-alcoholic fatty liver disease (NAFLD) and hepatocellular carcinoma (HCC) has been associated with altered expression of liver-specific genes including pyruvate kinase liver and red blood cell (PKLR), patatin-like phospholipase domain containing 3 (PNPLA3) and proprotein convertase subtilisin/kexin type 9 (PCSK9). Here, we inhibited and overexpressed the expression of these three genes in HepG2 cells, generated RNA-seq data before and after perturbation and revealed the altered global biological functions with the modulation of these genes using integrated network (IN) analysis. We found that modulation of these genes effects the total triglycerides levels within the cells and viability of the cells. Next, we generated IN for HepG2 cells, identified reporter transcription factors based on IN and found that the modulation of these genes affects key metabolic pathways associated with lipid metabolism (steroid biosynthesis, PPAR signalling pathway, fatty acid synthesis and oxidation) and cancer development (DNA replication, cell cycle and p53 signalling) involved in the progression of NAFLD and HCC. Finally, we observed that inhibition of PKLR lead to decreased glucose uptake and decreased mitochondrial activity in HepG2 cells. Hence, our systems level analysis indicated that PKLR can be targeted for development efficient treatment strategy for NAFLD and HCC.

Characterization of heterogeneous redox responses in hepatocellular carcinoma patients using network analysis

BACKGROUND

Redox metabolism is often considered a potential target for cancer treatment, but a systematic examination of redox responses in hepatocellular carcinoma (HCC) is missing.

METHODS

Here, we employed systems biology and biological network analyses to reveal key roles of genes associated with redox metabolism in HCC by integrating multi-omics data.

FINDINGS

We found that several redox genes, including 25 novel potential prognostic genes, are significantly co-expressed with liver-specific genes and genes associated with immunity and inflammation. Based on an integrative analysis, we found that HCC tumors display antagonistic behaviors in redox responses. The two HCC groups are associated with altered fatty acid, amino acid, drug and hormone metabolism, differentiation, proliferation, and NADPH-independent vs -dependent antioxidant defenses. Redox behavior varies with known tumor subtypes and progression, affecting patient survival. These antagonistic responses are also displayed at the protein and metabolite level and were validated in several independent cohorts. We finally showed the differential redox behavior using mice transcriptomics in HCC and noncancerous tissues and associated with hypoxic features of the two redox gene groups.

INTERPRETATION

Our integrative approaches highlighted mechanistic differences among tumors and allowed the identification of a survival signature and several potential therapeutic targets for the treatment of HCC.

Metabolic Network-Based Identification and Prioritization of Anticancer Targets Based on Expression Data in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is a deadly form of liver cancer with high mortality worldwide. Unfortunately, the large heterogeneity of this disease makes it difficult to develop effective treatment strategies. Cellular network analyses have been employed to study heterogeneity in cancer, and to identify potential therapeutic targets. However, the existing approaches do not consider metabolic growth requirements, i.e., biological network functionality, to rank candidate targets while preventing toxicity to non-cancerous tissues. Here, we developed an algorithm to overcome these issues based on integration of gene expression data, genome-scale metabolic models, network controllability, and dispensability, as well as toxicity analysis. This method thus predicts and ranks potential anticancer non-toxic controlling metabolite and gene targets. Our algorithm encompasses both objective-driven and-independent tasks, and uses network topology to finally rank the predicted therapeutic targets. We employed this algorithm to the analysis of transcriptomic data for 50 HCC patients with both cancerous and non-cancerous samples. We identified several potential targets that would prevent cell growth, including 74 anticancer metabolites, and 3 gene targets (PRKACA, PGS1, and CRLS1). The predicted anticancer metabolites showed good agreement with existing FDA-approved cancer drugs, and the 3 genes were experimentally validated by performing experiments in HepG2 and Hep3B liver cancer cell lines. Our observations indicate that our novel approach successfully identifies therapeutic targets for effective treatment of cancer. This approach may also be applied to any cancer type that has tumor and non-tumor gene or protein expression data.

Understanding the Representative Gut Microbiota Dysbiosis in Metformin-Treated Type 2 Diabetes Patients Using Genome-Scale Metabolic Modeling

Dysbiosis in the gut microbiome composition may be promoted by therapeutic drugs such as metformin, the world’s most prescribed antidiabetic drug. Under metformin treatment, disturbances of the intestinal microbes lead to increased abundance of Escherichia spp., Akkermansia muciniphila, Subdoligranulum variabile and decreased abundance of Intestinibacter bartlettii. This alteration may potentially lead to adverse effects on the host metabolism, with the depletion of butyrate producer genus. However, an increased production of butyrate and propionate was verified in metformin-treated Type 2 diabetes (T2D) patients. The mechanisms underlying these nutritional alterations and their relation with gut microbiota dysbiosis remain unclear. Here, we used Genome-scale Metabolic Models of the representative gut bacteria Escherichia spp., I. bartlettii, A. muciniphila, and S. variabile to elucidate their bacterial metabolism and its effect on intestinal nutrient pool, including macronutrients (e.g., amino acids and short chain fatty acids), minerals and chemical elements (e.g., iron and oxygen). We applied flux balance analysis (FBA) coupled with synthetic lethality analysis interactions to identify combinations of reactions and extracellular nutrients whose absence prevents growth. Our analyses suggest that Escherichia sp. is the bacteria least vulnerable to nutrient availability. We have also examined bacterial contribution to extracellular nutrients including short chain fatty acids, amino acids, and gasses. For instance, Escherichia sp. and S. variabile may contribute to the production of important short chain fatty acids (e.g., acetate and butyrate, respectively) involved in the host physiology under aerobic and anaerobic conditions. We have also identified pathway susceptibility to nutrient availability and reaction changes among the four bacteria using both FBA and flux variability analysis. For instance, lipopolysaccharide synthesis, nucleotide sugar metabolism, and amino acid metabolism are pathways susceptible to changes in Escherichia sp. and A. muciniphila. Our observations highlight important commensal and competing behavior, and their association with cellular metabolism for prevalent gut microbes. The results of our analysis have potential important implications for development of new therapeutic approaches in T2D patients through the development of prebiotics, probiotics, or postbiotics.

Systems biology in hepatology: approaches and applications

Detailed insights into the biological functions of the liver and an understanding of its crosstalk with other human tissues and the gut microbiota can be used to develop novel strategies for the prevention and treatment of liver-associated diseases, including fatty liver disease, cirrhosis, hepatocellular carcinoma and type 2 diabetes mellitus. Biological network models, including metabolic, transcriptional regulatory, protein-protein interaction, signalling and co-expression networks, can provide a scaffold for studying the biological pathways operating in the liver in connection with disease development in a systematic manner. Here, we review studies in which biological network models were used to integrate multiomics data to advance our understanding of the pathophysiological responses of complex liver diseases. We also discuss how this mechanistic approach can contribute to the discovery of potential biomarkers and novel drug targets, which might lead to the design of targeted and improved treatment strategies. Finally, we present a roadmap for the successful integration of models of the liver and other human tissues with the gut microbiota to simulate whole-body metabolic functions in health and disease.