Integrative system biology and mathematical modeling of genetic networks identifies shared biomarkers for obesity and diabetes

<abstract> <p>Obesity and type 2 and diabetes mellitus (T2D) are two dual epidemics whose shared genetic pathological mechanisms are still far from being fully understood. Therefore, this study is aimed at discovering key genes, molecular mechanisms, and new drug targets for obesity and T2D by analyzing the genome wide gene expression data with different computational biology approaches. In this study, the RNA-sequencing data of isolated primary human adipocytes from individuals who are lean, obese, and T2D was analyzed by an integrated framework consisting of gene expression, protein interaction network (PIN), tissue specificity, and druggability approaches. Our findings show a total of 1932 unique differentially expressed genes (DEGs) across the diabetes versus obese group comparison (p≤0.05). The PIN analysis of these 1932 DEGs identified 190 high centrality network (HCN) genes, which were annotated against 3367 GO terms and functional pathways, like response to insulin signaling, phosphorylation, lipid metabolism, glucose metabolism, etc. (p≤0.05). By applying additional PIN and topological parameters to 190 HCN genes, we further mapped 25 high confidence genes, functionally connected with diabetes and obesity traits. Interestingly, <italic>ERBB2, FN1, FYN, HSPA1A, HBA1</italic>, and <italic>ITGB1</italic> genes were found to be tractable by small chemicals, antibodies, and/or enzyme molecules. In conclusion, our study highlights the potential of computational biology methods in correlating expression data to topological parameters, functional relationships, and druggability characteristics of the candidate genes involved in complex metabolic disorders with a common etiological basis.</p> </abstract>

The Reasonable Effectiveness of Randomness in Scalable and Integrative Gene Regulatory Network Inference and Beyond

Gene regulation is orchestrated by a vast number of molecules, including transcription factors and co-factors, chromatin regulators, as well as epigenetic mechanisms, and it has been shown that transcriptional misregulation, e.g., caused by mutations in regulatory sequences, is responsible for a plethora of diseases, including cancer, developmental or neurological disorders. As a consequence, decoding the architecture of gene regulatory networks has become one of the most important tasks in modern (computational) biology. However, to advance our understanding of the mechanisms involved in the transcriptional apparatus, we need scalable approaches that can deal with the increasing number of large-scale, high-resolution, biological datasets. In particular, such approaches need to be capable of efficiently integrating and exploiting the biological and technological heterogeneity of such datasets in order to best infer the underlying, highly dynamic regulatory networks, often in the absence of sufficient ground truth data for model training or testing. With respect to scalability, randomized approaches have proven to be a promising alternative to deterministic methods in computational biology. As an example, one of the top performing algorithms in a community challenge on gene regulatory network inference from transcriptomic data is based on a random forest regression model. In this concise survey, we aim to highlight how randomized methods may serve as a highly valuable tool, in particular, with increasing amounts of large-scale, biological experiments and datasets being collected. Given the complexity and interdisciplinary nature of the gene regulatory network inference problem, we hope our survey maybe helpful to both computational and biological scientists. It is our aim to provide a starting point for a dialogue about the concepts, benefits, and caveats of the toolbox of randomized methods, since unravelling the intricate web of highly dynamic, regulatory events will be one fundamental step in understanding the mechanisms of life and eventually developing efficient therapies to treat and cure diseases.

Current and emerging tools of computational biology to improve the detoxification of mycotoxins

Biological organisms carry a rich potential for removing toxins from our environment, but identifying suitable candidates and improving them remain challenging. We explore the use of computational tools to discover strains and enzymes that detoxify harmful compounds. In particular, we will focus on mycotoxins—fungi-produced toxins that contaminate food and feed—and biological enzymes that are capable of rendering them less harmful. We discuss the use of established and novel computational tools to complement existing empirical data in three directions: discovering the prospect of detoxification among underexplored organisms, finding important cellular processes that contribute to detoxification, and improving the performance of detoxifying enzymes. We hope to create a synergistic conversation between researchers in computational biology and those in the bioremediation field. We showcase open bioremediation questions where computational researchers can contribute and highlight relevant existing and emerging computational tools that could benefit bioremediation researchers.

CMash: fast, multi-resolution estimation of k-mer-based Jaccard and containment indices

K-mer based methods are used ubiquitously in the field of computational biology. However, determining the optimal value of k for a specific application often remains heuristic. Simply reconstructing a new k-mer set with another k-mer size is computationally expensive, especially in metagenomic analysis where data sets are large. Here, we introduce a hashing-based technique that leverages a kind of bottom-m sketch as well as a k-mer ternary search tree (KTST) to obtain k-mer based similarity estimates for a range of k values. By truncating k-mers stored in a pre-built KTST with a large k = kmax value, we can simultaneously obtain k-mer based estimates for all k values up to kmax. This truncation approach circumvents the reconstruction of new k-mer sets when changing k values, making analysis more time and space-efficient. For example, we show that when using a KTST to estimate the containment index between a RefSeq-based microbial reference database and simulated metagenome data for 10 values of k, the running time is close to 10x faster compared to a classic MinHash approach while using less than one-fifth the space to store the data structure. A python implementation of this method, CMash, is available at https://github.com/dkoslicki/CMash. The reproduction of all experiments presented herein can be accessed via https://github.com/KoslickiLab/CMASH-reproducibles.

HLAncPred: A method for predicting promiscuous non-classical HLA binding sites

In the last two decades, ample of methods have been developed to predict the classical HLA binders in an antigen. In contrast, limited attempts have been made to develop methods for predicting binders for non-classical HLA; due to the scarcity of sufficient experimental data and lack of community interest. Of Note, non-classical HLA plays a crucial immunomodulatory role and regulates various immune responses. Recent studies revealed that non-classical HLA (HLA-E & HLA-G) based immunotherapies have many advantages over classical HLA based-immunotherapy, particularly against COVID-19. In order to facilitate the scientific community, we have developed an artificial intelligence-based method for predicting binders of non-classical HLA alleles (HLA-G and HLA-E). All the models were trained and tested on experimentally validated data obtained from the recent release of IEDB. The machine learning based-models achieved more than 0.98 AUC for HLA-G alleles on validation or independent dataset. Similarly, our models achieved the highest AUC of 0.96 and 0.88 on the validation dataset for HLA-E*01:01, HLA-E*01:03, respectively. We have summarized the models developed in the past for non-classical HLA binders and compared with the models developed in this study. Moreover, we have also predicted the non-classical HLA binders in the spike protein of different variants of virus causing COVID-19 including omicron (B.1.1.529) to facilitate the community. One of the major challenges in the field of immunotherapy is to identify the promiscuous binders or antigenic regions that can bind to a large number of HLA alleles. In order to predict the promiscuous binders for the non-classical HLA alleles, we developed a web server HLAncPred (https://webs.iiitd.edu.in/raghava/hlancpred), and a standalone package.Key PointsNon-classical HLAs play immunomodulatory roles in the immune system.HLA-E restricted T-cell therapy may reduce COVID-19 associated cytokine storm.In silico models developed for predicting binders for HLA-G and HLA-E.Identification of non-classical HLA binders in strains of coronavirusA webserver for predicting promiscuous binders for non-classical HLA allelesAuthor’s BiographyAnjali Dhall is currently working as Ph.D. in Bioinformatics from Department of Computational Biology, Indraprastha Institute of Information Technology, New Delhi, India.Sumeet Patiyal is currently working as Ph.D. in Bioinformatics from Department of Computational Biology, Indraprastha Institute of Information Technology, New Delhi, India.Gajendra P. S. Raghava is currently working as Professor and Head of Department of Computational Biology, Indraprastha Institute of Information Technology, New Delhi, India.

Targeting CK2 mediated signaling to impair/tackle SARS-CoV-2 infection: a computational biology approach

Background Similarities in the hijacking mechanisms used by SARS-CoV-2 and several types of cancer, suggest the repurposing of cancer drugs to treat Covid-19. CK2 kinase antagonists have been proposed for cancer treatment. A recent study in cells infected with SARS-CoV-2 found a significant CK2 kinase activity, and the use of a CK2 inhibitor showed antiviral responses. CIGB-300, originally designed as an anticancer peptide, is an antagonist of CK2 kinase activity that binds to the CK2 phospho-acceptor sites. Recent preliminary results show the antiviral activity of CIGB-300 using a surrogate model of coronavirus. Here we present a computational biology study that provides evidence, at the molecular level, of how CIGB-300 may interfere with the SARS-CoV-2 life cycle within infected human cells. Methods Sequence analyses and data from phosphorylation studies were combined to predict infection-induced molecular mechanisms that can be interfered by CIGB-300. Next, we integrated data from multi-omics studies and data focusing on the antagonistic effect on the CK2 kinase activity of CIGB-300. A combination of network and functional enrichment analyses was used. Results Firstly, from the SARS-CoV studies, we inferred the potential incidence of CIGB-300 in SARS-CoV-2 interference on the immune response. Afterwards, from the analysis of multiple omics data, we proposed the action of CIGB-300 from the early stages of viral infections perturbing the virus hijacking of RNA splicing machinery. We also predicted the interference of CIGB-300 in virus-host interactions that are responsible for the high infectivity and the particular immune response to SARS-CoV-2 infection. Furthermore, we provided evidence of how CIGB-300 may participate in the attenuation of phenotypes related to muscle, bleeding, coagulation and respiratory disorders. Conclusions Our computational analysis proposes putative molecular mechanisms that support the antiviral activity of CIGB-300.