Finding undetected protein associations in cell signaling by belief propagation

External information propagates in the cell mainly through signaling cascades and transcriptional activation, allowing it to react to a wide spectrum of environmental changes. High throughput experiments identify numerous molecular components of such cascades that may, however, interact through unknown partners. Some of them may be detected using data coming from the integration of a protein-protein interaction network and mRNA expression profiles. This inference problem can be mapped onto the problem of finding appropriate optimal connected subgraphs of a network defined by these datasets. The optimization procedure turns out to be computationally intractable in general. Here we present a new distributed algorithm for this task, inspired from statistical physics, and apply this scheme to alpha factor and drug perturbations data in yeast. We identify the role of the COS8 protein, a member of a gene family of previously unknown function, and validate the results by genetic experiments. The algorithm we present is specially suited for very large datasets, can run in parallel, and can be adapted to other problems in systems biology. On renowned benchmarks it outperforms other algorithms in the field.Comment: 6 pages, 3 figures, 1 table, Supporting Informatio

Mechanistic interplay between ceramide and insulin resistance

Recent research adds to a growing body of literature on the essential role of ceramides in glucose homeostasis and insulin signaling, while the mechanistic interplay between various components of ceramide metabolism remains to be quantified. We present an extended model of C16:0 ceramide production through both the de novo synthesis and the salvage pathways. We verify our model with a combination of published models and independent experimental data. In silico experiments of the behavior of ceramide and related bioactive lipids in accordance with the observed transcriptomic changes in obese/diabetic murine macrophages at 5 and 16 weeks support the observation of insulin resistance only at the later phase. Our analysis suggests the pivotal role of ceramide synthase, serine palmitoyltransferase and dihydroceramide desaturase involved in the de novo synthesis and the salvage pathways in influencing insulin resistance versus its regulation

Discovery of Sexual Dimorphisms in Metabolic and Genetic Biomarkers

Metabolomic profiling and the integration of whole-genome genetic association data has proven to be a powerful tool to comprehensively explore gene regulatory networks and to investigate the effects of genetic variation at the molecular level. Serum metabolite concentrations allow a direct readout of biological processes, and association of specific metabolomic signatures with complex diseases such as Alzheimer's disease and cardiovascular and metabolic disorders has been shown. There are well-known correlations between sex and the incidence, prevalence, age of onset, symptoms, and severity of a disease, as well as the reaction to drugs. However, most of the studies published so far did not consider the role of sexual dimorphism and did not analyse their data stratified by gender. This study investigated sex-specific differences of serum metabolite concentrations and their underlying genetic determination. For discovery and replication we used more than 3,300 independent individuals from KORA F3 and F4 with metabolite measurements of 131 metabolites, including amino acids, phosphatidylcholines, sphingomyelins, acylcarnitines, and C6-sugars. A linear regression approach revealed significant concentration differences between males and females for 102 out of 131 metabolites (p-values<3.8 x 10(-4); Bonferroni-corrected threshold). Sex-specific genome-wide association studies (GWAS) showed genome-wide significant differences in beta-estimates for SNPs in the CPS1 locus (carbamoyl-phosphate synthase 1, significance level: p<3.8 x 10(-10); Bonferroni-corrected threshold) for glycine. We showed that the metabolite profiles of males and females are significantly different and, furthermore, that specific genetic variants in metabolism-related genes depict sexual dimorphism. Our study provides new important insights into sex-specific differences of cell regulatory processes and underscores that studies should consider sex-specific effects in design and interpretation

Discovery of Sexual Dimorphisms in Metabolic and Genetic Biomarkers

Metabolomic profiling and the integration of whole-genome genetic association data has proven to be a powerful tool to comprehensively explore gene regulatory networks and to investigate the effects of genetic variation at the molecular level. Serum metabolite concentrations allow a direct readout of biological processes, and association of specific metabolomic signatures with complex diseases such as Alzheimer's disease and cardiovascular and metabolic disorders has been shown. There are well-known correlations between sex and the incidence, prevalence, age of onset, symptoms, and severity of a disease, as well as the reaction to drugs. However, most of the studies published so far did not consider the role of sexual dimorphism and did not analyse their data stratified by gender. This study investigated sex-specific differences of serum metabolite concentrations and their underlying genetic determination. For discovery and replication we used more than 3,300 independent individuals from KORA F3 and F4 with metabolite measurements of 131 metabolites, including amino acids, phosphatidylcholines, sphingomyelins, acylcarnitines, and C6-sugars. A linear regression approach revealed significant concentration differences between males and females for 102 out of 131 metabolites (p-values<3.8 x 10(-4); Bonferroni-corrected threshold). Sex-specific genome-wide association studies (GWAS) showed genome-wide significant differences in beta-estimates for SNPs in the CPS1 locus (carbamoyl-phosphate synthase 1, significance level: p<3.8 x 10(-10); Bonferroni-corrected threshold) for glycine. We showed that the metabolite profiles of males and females are significantly different and, furthermore, that specific genetic variants in metabolism-related genes depict sexual dimorphism. Our study provides new important insights into sex-specific differences of cell regulatory processes and underscores that studies should consider sex-specific effects in design and interpretation

DCcov: Repositioning of Drugs and Drug Combinations for SARS-CoV-2 Infected Lung through Constraint-Based Modelling

The 2019 coronavirus disease (COVID-19) became a worldwide pandemic with currently no effective antiviral drug except treatments for symptomatic therapy. Flux balance analysis is an efficient method to analyze metabolic networks. It allows optimizing for a metabolic function and thus e.g., predicting the growth rate of a specific cell or the production rate of a metabolite of interest. Here flux balance analysis was applied on human lung cells infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to reposition metabolic drugs and drug combinations against the replication of the SARS-CoV-2 virus within the host tissue. Making use of expression data sets of infected lung tissue, genome-scale COVID-19-specific metabolic models were reconstructed. Then host-specific essential genes and gene-pairs were determined through in-silico knockouts that permit reducing the viral biomass production without affecting the host biomass. Key pathways that are associated with COVID-19 severity in lung tissue are related to oxidative stress, as well as ferroptosis, sphingolipid metabolism, cysteine metabolism, and fat digestion. By in-silico screening of FDA approved drugs on the putative disease-specific essential genes and gene-pairs, 45 drugs and 99 drug combinations were predicted as promising candidates for COVID-19 focused drug repositioning (https://github.com/sysbiolux/DCcov). Among the 45 drug candidates, six antiviral drugs were found and seven drugs that are being tested in clinical trials against COVID-19. Other drugs like gemcitabine, rosuvastatin and acetylcysteine, and drug combinations like azathioprine-pemetrexed might offer new chances for treating COVID-19

Application of bioinformatics in studies of sphingolipid biosynthesis

The studies in this dissertation demonstrate that the gene expression pathway maps are useful tools to notice alteration in different branches of sphingolipid biosynthesis pathway based on microarray and other transcriptomic analysis. To facilitate the integrative analysis of gene expression and sphingolipid amounts, updated pathway maps were prepared using an open access visualization tool, Pathvisio v1.1. The datasets were formatted using Perl scripts and visualized with the aid of color coded pathway diagrams. Comparative analysis of transcriptomics and sphingolipid alterations from experimental studies and published literature revealed 72.8 % correlation between mRNA and sphingolipid differences (p-value < 0.0001 by the Fisher's exact test).The high correlation between gene expression differences and sphingolipid alterations highlights the application of this tool to evaluate molecular changes associate with sphingolipid alterations as well as predict differences in specific metabolites that can be experimentally verified using sensitive approaches such as mass spectrometry. In addition, bioinformatics sequence analysis was used to identify transcripts for sphingolipid biosynthesis enzyme 3-ketosphinganine reductase, and homology modeling studies helped in the evaluation of a cell line defective in sphingolipid metabolism due to mutation in the enzyme serine palmitoyltransferase, the first enzyme of de novo biosynthesis pathway. Hence, the combination of different bioinformatics approaches, including protein and DNA sequence analysis, structure modeling and pathway diagrams can provide valuable inputs for biochemical and molecular studies of sphingolipid metabolism.Ph.D.Committee Chair: Alfred H Merrill Jr; Committee Member: Cameron Sullards; Committee Member: I King Jordan; Committee Member: Marion B. Sewer; Committee Member: Stephen C Harve

Lipidomics in Health and Diseases - Beyond the Analysis of Lipids

published_or_final_versio

Integrated human/SARS-CoV-2 metabolic models present novel treatment strategies against COVID-19

The coronavirus disease 2019 (COVID-19) pandemic caused by the new coronavirus (SARS-CoV-2) is currently responsible for more than 3 million deaths in 219 countries across the world and with more than 140 million cases. The absence of FDA-approved drugs against SARS-CoV-2 has highlighted an urgent need to design new drugs. We developed an integrated model of the human cell and SARS-CoV-2 to provide insight into the virus'' pathogenic mechanism and support current therapeutic strategies. We show the biochemical reactions required for the growth and general maintenance of the human cell, first, in its healthy state. We then demonstrate how the entry of SARS-CoV-2 into the human cell causes biochemical and structural changes, leading to a change of cell functions or cell death. A new computational method that predicts 20 unique reactions as drug targets from our models and provides a platform for future studies on viral entry inhibition, immune regulation, and drug optimisation strategies. The model is available in BioModels (https://www.ebi.ac.uk/biomodels/MODEL2007210001) and the software tool, findCPcli, that implements the computational method is available at https://github.com/findCP/findCPcli. © 2021 Bannerman et al