GeNeCK: a web server for gene network construction and visualization

Review on computer-assisted biosynthetic capacities elucidation to assess metabolic interactions and communication within microbial communities

Microbial communities thrive through interactions and communication, which are challenging to study as most microorganisms are not cultivable. To address this challenge, researchers focus on the extracellular space where communication events occur. Exometabolomics and interactome analysis provide insights into the molecules involved in communication and the dynamics of their interactions. Advances in sequencing technologies and computational methods enable the reconstruction of taxonomic and functional profiles of microbial communities using high-throughput multi-omics data. Network-based approaches, including community flux balance analysis, aim to model molecular interactions within and between communities. Despite these advances, challenges remain in computer-assisted biosynthetic capacities elucidation, requiring continued innovation and collaboration among diverse scientists. This review provides insights into the current state and future directions of computer-assisted biosynthetic capacities elucidation in studying microbial communities.

Bioinformatics tools and resources for cancer and application

ABSTRACT

Tumor bioinformatics plays an important role in cancer research and precision medicine. The primary focus of traditional cancer research has been molecular and clinical studies of a number of fundamental pathways and genes. In recent years, driven by breakthroughs in high-throughput technologies, large-scale cancer omics data have accumulated rapidly. How to effectively utilize and share these data is particularly important. To address this crucial task, many computational tools and databases have been developed over the past few years. To help researchers quickly learn and understand the functions of these tools, in this review, we summarize publicly available bioinformatics tools and resources for pan-cancer multi-omics analysis, regulatory analysis of tumorigenesis, tumor treatment and prognosis, immune infiltration analysis, immune repertoire analysis, cancer driver gene and driver mutation analysis, and cancer single-cell analysis, which may further help researchers find more suitable tools for their research.

Multioviz: an interactive platform for in silico perturbation and interrogation of gene regulatory networks

In this paper, we aim to build a platform that will help bridge the gap between high-dimensional computation and wet-lab experimentation by allowing users to interrogate genomic signatures at multiple molecular levels and identify best next actionable steps for downstream decision making. We introduce Multioviz: a publicly accessible R package and web application platform to easily perform in silico hypothesis testing of generated gene regulatory networks. We demonstrate the utility of Multioviz by conducting an end-to-end analysis in a statistical genetics application focused on measuring the effect of in silico perturbations of complex trait architecture. By using a real dataset from the Wellcome Trust Centre for Human Genetics, we both recapitulate previous findings and propose hypotheses about the genes involved in the percentage of immune CD8+ cells found in heterogeneous stocks of mice. Source code for the Multioviz R package is available at https://github.com/lcrawlab/multio-viz and an interactive version of the platform is available at https://multioviz.ccv.brown.edu/ .

Specific protein interactions between rice members of the GT43 and GT47 families form various central cores of putative xylan synthase complexes

Members of the glycosyltransferase (GT)43 and GT47 families have been associated with heteroxylan synthesis in both dicots and monocots and are thought to assemble into central cores of putative xylan synthase complexes (XSCs). Currently, it is unknown whether protein-protein interactions within these central cores are specific, how many such complexes exist, and whether these complexes are functionally redundant. Here, we used gene association network and co-expression approaches in rice to identify four OsGT43s and four OsGT47s that assemble into different GT43/GT47 complexes. Using two independent methods, we showed that (i) these GTs assemble into at least six unique complexes through specific protein-protein interactions and (ii) the proteins interact directly in vitro. Confocal microscopy showed that, when alone, all OsGT43s were retained in the endoplasmic reticulum (ER), while all OsGT47s were localized in the Golgi. co-expression of OsGT43s and OsGT47s displayed complexes that form in the ER but accumulate in Golgi. ER-to-Golgi trafficking appears to require interactions between OsGT43s and OsGT47s. Comparison of the central cores of the three putative rice OsXSCs to wheat, asparagus, and Arabidopsis XSCs, showed great variation in GT43/GT47 combinations, which makes the identification of orthologous central cores between grasses and dicots challenging. However, the emerging picture is that all central cores from these species seem to have at least one member of the IRX10/IRX10-L clade in the GT47 family in common, suggesting greater functional importance for this family in xylan synthesis. Our findings provide a new framework for future investigation of heteroxylan biosynthesis and function in monocots.

Mathematical modeling of regulatory networks of intracellular processes - Aims and selected methods

Regulatory networks structure and signaling pathways dynamics are uncovered in time- and resource consuming experimental work. However, it is increasingly supported by modeling, analytical and computational techniques as well as discrete mathematics and artificial intelligence applied to to extract knowledge from existing databases. This review is focused on mathematical modeling used to analyze dynamics and robustness of these networks. This paper presents a review of selected modeling methods that facilitate advances in molecular biology.

Identification and characterization of sex-dependent gene expression profile in glioblastoma

Glioblastoma (GBM) is the most lethal primary tumor in the human brain and lacks favorable treatment options. Sex differences in the outcome of GBM are broadly acknowledged, but the underlying molecular mechanisms remain largely unknown. To identify the sex-dependent critical genes in the progression of GBM, raw data from several microarray datasets with the same array platform were downloaded from the Gene Expression Omnibus (GEO) database. These datasets included tumorous and normal tissue from patients with GBM and crucial sex features. Then, the differentially expressed genes (DEGs) in female and male tumors were identified via bioinformatics analysis, respectively. Functional signatures of the identified DEGs were further annotated by Gene Ontology (GO) and pathway enrichment analyses. Venn diagram and functional protein-protein interaction (PPI) network analyses were performed to screen out the sex-specific DEGs. Survival analysis of patients with differences in the expression level of selected genes was then carried out using the data from The Cancer Genome Atlas (TCGA). Here, we showed that ECT2, AURKA, TYMS, CDK1, NCAPH, CENPU, OIP5, KIF14, ASPM, FBXO5, SGOL2, CASC5, SHCBP1, FN1, LOX, IGFBP3, CSPG4, and CD44 were enriched in female tumor samples, whereas TNFSF13B, CXCL10, CXCL8, CXCR4, TLR2, CCL2, and FCGR2A were enriched in male tumor samples. Among these key genes, interestingly, ECT2 was associated with increased an survival rate for female patients, whileTNFSF13B could be regarded as a potential marker of poor prognosis in male patients. These results suggested that sex differences in patients may be attributed to the heterogeneous gene activity, which might influence the oncogenesis and the outcomes of GBM.

Unraveling Plant-Pathogen Interactions in Cereals Using RNA-seq

Over the past two decades, there have been significant advancements in the realm of transcriptomics, or the study of genes and their expression. Modern RNA sequencing technologies and high-performance computing are creating a "big data" revolution that provides new opportunities to explore the interactions between cereals and pathogens that affect grain yield and food safety. These data are being used to annotate genes and gene variants, as well as identify differentially expressed genes and create global gene co-expression networks. Moreover, these data can unravel the complex interactions between pathogen and host and identify genes and pathways involved in these interactions. This information can then be used for disease mitigation and the development of crops with superior resistance.

Inferencing Bulk Tumor and Single-Cell Multi-Omics Regulatory Networks for Discovery of Biomarkers and Therapeutic Targets

There are insufficient accurate biomarkers and effective therapeutic targets in current cancer treatment. Multi-omics regulatory networks in patient bulk tumors and single cells can shed light on molecular disease mechanisms. Integration of multi-omics data with large-scale patient electronic medical records (EMRs) can lead to the discovery of biomarkers and therapeutic targets. In this review, multi-omics data harmonization methods were introduced, and common approaches to molecular network inference were summarized. Our Prediction Logic Boolean Implication Networks (PLBINs) have advantages over other methods in constructing genome-scale multi-omics networks in bulk tumors and single cells in terms of computational efficiency, scalability, and accuracy. Based on the constructed multi-modal regulatory networks, graph theory network centrality metrics can be used in the prioritization of candidates for discovering biomarkers and therapeutic targets. Our approach to integrating multi-omics profiles in a patient cohort with large-scale patient EMRs such as the SEER-Medicare cancer registry combined with extensive external validation can identify potential biomarkers applicable in large patient populations. These methodologies form a conceptually innovative framework to analyze various available information from research laboratories and healthcare systems, accelerating the discovery of biomarkers and therapeutic targets to ultimately improve cancer patient survival outcomes.

Angiogenesis goes computational - The future way forward to discover new angiogenic targets?

Multi-omics technologies are being increasingly utilized in angiogenesis research. Yet, computational methods have not been widely used for angiogenic target discovery and prioritization in this field, partly because (wet-lab) vascular biologists are insufficiently familiar with computational biology tools and the opportunities they may offer. With this review, written for vascular biologists who lack expertise in computational methods, we aspire to break boundaries between both fields and to illustrate the potential of these tools for future angiogenic target discovery. We provide a comprehensive survey of currently available computational approaches that may be useful in prioritizing candidate genes, predicting associated mechanisms, and identifying their specificity to endothelial cell subtypes. We specifically highlight tools that use flexible, machine learning frameworks for large-scale data integration and gene prioritization. For each purpose-oriented category of tools, we describe underlying conceptual principles, highlight interesting applications and discuss limitations. Finally, we will discuss challenges and recommend some guidelines which can help to optimize the process of accurate target discovery.

Gene regulatory network inference methodology for genomic and transcriptomic data acquired in genetically related heterozygote individuals

MOTIVATION

Inferring gene regulatory networks in non-independent genetically related panels is a methodological challenge. This hampers evolutionary and biological studies using heterozygote individuals such as in wild sunflower populations or cultivated hybrids.

RESULTS

First, we simulated 100 datasets of gene expressions and polymorphisms, displaying the same gene expression distributions, heterozygosities and heritabilities as in our dataset including 173 genes and 353 genotypes measured in sunflower hybrids. Secondly, we performed a meta-analysis based on six inference methods [least absolute shrinkage and selection operator (Lasso), Random Forests, Bayesian Networks, Markov Random Fields, Ordinary Least Square and fast inference of networks from directed regulation (Findr)] and selected the minimal density networks for better accuracy with 64 edges connecting 79 genes and 0.35 area under precision and recall (AUPR) score on average. We identified that triangles and mutual edges are prone to errors in the inferred networks. Applied on classical datasets without heterozygotes, our strategy produced a 0.65 AUPR score for one dataset of the DREAM5 Systems Genetics Challenge. Finally, we applied our method to an experimental dataset from sunflower hybrids. We successfully inferred a network composed of 105 genes connected by 106 putative regulations with a major connected component.

AVAILABILITY AND IMPLEMENTATION

Our inference methodology dedicated to genomic and transcriptomic data is available at https://forgemia.inra.fr/sunrise/inference_methods.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Review and assessment of Boolean approaches for inference of gene regulatory networks

Boolean descriptions of gene regulatory networks can provide an insight into interactions between genes. Boolean networks hold predictive power, are easy to understand, and can be used to simulate the observed networks in different scenarios. We review fundamental and state-of-the-art methods for inference of Boolean networks. We introduce a methodology for a straightforward evaluation of Boolean inference approaches based on the generation of evaluation datasets, application of selected inference methods, and evaluation of performance measures to guide the selection of the best method for a given inference problem. We demonstrate this procedure on inference methods REVEAL (REVerse Engineering ALgorithm), Best-Fit Extension, MIBNI (Mutual Information-based Boolean Network Inference), GABNI (Genetic Algorithm-based Boolean Network Inference) and ATEN (AND/OR Tree ENsemble algorithm), which infers Boolean descriptions of gene regulatory networks from discretised time series data. Boolean inference approaches tend to perform better in terms of dynamic accuracy, and slightly worse in terms of structural correctness. We believe that the proposed methodology and provided guidelines will help researchers to develop Boolean inference approaches with a good predictive capability while maintaining structural correctness and biological relevance.

NETISCE: a network-based tool for cell fate reprogramming

The search for effective therapeutic targets in fields like regenerative medicine and cancer research has generated interest in cell fate reprogramming. This cellular reprogramming paradigm can drive cells to a desired target state from any initial state. However, methods for identifying reprogramming targets remain limited for biological systems that lack large sets of experimental data or a dynamical characterization. We present NETISCE, a novel computational tool for identifying cell fate reprogramming targets in static networks. In combination with machine learning algorithms, NETISCE estimates the attractor landscape and predicts reprogramming targets using signal flow analysis and feedback vertex set control, respectively. Through validations in studies of cell fate reprogramming from developmental, stem cell, and cancer biology, we show that NETISCE can predict previously identified cell fate reprogramming targets and identify potentially novel combinations of targets. NETISCE extends cell fate reprogramming studies to larger-scale biological networks without the need for full model parameterization and can be implemented by experimental and computational biologists to identify parts of a biological system relevant to the desired reprogramming task.

Differential Expression of Genes Regulating Store-operated Calcium Entry in Conjunction With Mitochondrial Dynamics as Potential Biomarkers for Cancer: A Single-Cell RNA Analysis

Regulation of intracellular concentration of calcium levels is crucial for cell signaling, homeostasis, and in the pathology of diseases including cancer. Agonist-induced entry of calcium ions into the non-excitable cells is mediated by store-operated calcium channels (SOCs). This pathway is activated by the release of calcium ions from the endoplasmic reticulum and further regulated by the calcium uptake through mitochondria leading to calcium-dependent inactivation of calcium-release activated calcium channels (CARC). SOCs including stromal interaction molecules (STIM) and ORAI proteins have been implicated in tumor growth, progression, and metastasis. In the present study, we analyzed the mRNA and protein expression of genes mediating SOCs-STIM1, STIM2, ORAI1, ORAI2, ORAI3, TRPC1, TRPC3, TRPC4, TRPC5, TRPC6, TRPC7, TRPV1, TRPV2, TRPM1, and TRPM7 in head and neck squamous cell cancer (HNSC) patients using TCGA and CPTAC analysis. Further, our in silico analysis showed a significant correlation between the expression of SOCs and genes involved in the mitochondrial dynamics (MDGs) both at mRNA and protein levels. Protein-protein docking results showed lower binding energy for SOCs with MDGs. Subsequently, we validated these results using gene expression and single-cell RNA sequencing datasets retrieved from Gene Expression Omnibus (GEO). Single-cell gene expression analysis of HNSC tumor tissues revealed that SOCs expression is remarkably associated with the MDGs expression in both cancer and fibroblast cells.

Gene expression data of inflammatory mediators in apical periodontitis in 129 (wild type) and 5-lipoxygenase knockout mice

Apical periodontitis is an immune inflammatory response around periapical tissues as a result of pathogens invasion into the root canal. The host immunoinflammatory response could determine the progression of this disease, which involves the recruitment of immune cells, and the release of several cytokines in the lesion site. The 5-lipoxygenase pathway has been activated in some osteolytic diseases due to its capacity to interfere in the proliferation and differentiation of bone cells, including the osteoclasts. As mean to understand the inflammatory genes regulation in the apical periodontitis progression, we evaluated the network of 66 genes related to cytokines, chemokines and other inflammatory mediators and receptors in the wild-type (WT) and 5-lipoxygenase enzyme genetically deficient mice (KO). This article presents data not published but related to the research article “Effects of 5-lipoxygenase gene disruption on inflammation, osteoclastogenesis and bone resorption in polymicrobial apical periodontitis” .

MODalyseR-a novel software for inference of disease module hub regulators identified a putative multiple sclerosis regulator supported by independent eQTL data

Motivation

Network-based disease modules have proven to be a powerful concept for extracting knowledge about disease mechanisms, predicting for example disease risk factors and side effects of treatments. Plenty of tools exist for the purpose of module inference, but less effort has been put on simultaneously utilizing knowledge about regulatory mechanisms for predicting disease module hub regulators.

Results

We developed MODalyseR, a novel software for identifying disease module regulators and reducing modules to the most disease-associated genes. This pipeline integrates and extends previously published software packages MODifieR and ComHub and hereby provides a user-friendly network medicine framework combining the concepts of disease modules and hub regulators for precise disease gene identification from transcriptomics data. To demonstrate the usability of the tool, we designed a case study for multiple sclerosis that revealed IKZF1 as a promising hub regulator, which was supported by independent ChIP-seq data.

Availability and implementation

MODalyseR is available as a Docker image at https://hub.docker.com/r/ddeweerd/modalyser with user guide and installation instructions found at https://gustafsson-lab.gitlab.io/MODalyseR/.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

A Web Tool for Consensus Gene Regulatory Network Construction

Gene regulatory network (GRN) construction involves various steps of complex computational steps. This step-by-step procedure requires prior knowledge of programming languages such as R. Development of a web tool may reduce this complexity in the analysis steps which can be easy accessible for the user. In this study, a web tool for constructing consensus GRN by combining the outcomes obtained from four methods, namely, correlation, principal component regression, partial least square, and ridge regression, has been developed. We have designed the web tool with an interactive and user-friendly web page using the php programming language. We have used R script for the analysis steps which run in the background of the user interface. Users can upload gene expression data for constructing consensus GRN. The output obtained from analysis will be available in downloadable form in the result window of the web tool.

Tailored graphical lasso for data integration in gene network reconstruction

BACKGROUND

Identifying gene interactions is a topic of great importance in genomics, and approaches based on network models provide a powerful tool for studying these. Assuming a Gaussian graphical model, a gene association network may be estimated from multiomic data based on the non-zero entries of the inverse covariance matrix. Inferring such biological networks is challenging because of the high dimensionality of the problem, making traditional estimators unsuitable. The graphical lasso is constructed for the estimation of sparse inverse covariance matrices in such situations, using [Formula: see text]-penalization on the matrix entries. The weighted graphical lasso is an extension in which prior biological information from other sources is integrated into the model. There are however issues with this approach, as it naïvely forces the prior information into the network estimation, even if it is misleading or does not agree with the data at hand. Further, if an associated network based on other data is used as the prior, the method often fails to utilize the information effectively.

RESULTS

We propose a novel graphical lasso approach, the tailored graphical lasso, that aims to handle prior information of unknown accuracy more effectively. We provide an R package implementing the method, tailoredGlasso. Applying the method to both simulated and real multiomic data sets, we find that it outperforms the unweighted and weighted graphical lasso in terms of all performance measures we consider. In fact, the graphical lasso and weighted graphical lasso can be considered special cases of the tailored graphical lasso, and a parameter determined by the data measures the usefulness of the prior information. We also find that among a larger set of methods, the tailored graphical is the most suitable for network inference from high-dimensional data with prior information of unknown accuracy. With our method, mRNA data are demonstrated to provide highly useful prior information for protein-protein interaction networks.

CONCLUSIONS

The method we introduce utilizes useful prior information more effectively without involving any risk of loss of accuracy should the prior information be misleading.

In silico-driven analysis of the Glossina morsitans morsitans antennae transcriptome in response to repellent or attractant compounds

BACKGROUND

High-throughput sequencing generates large volumes of biological data that must be interpreted to make meaningful inference on the biological function. Problems arise due to the large number of characteristics p (dimensions) that describe each record [n] in the database. Feature selection using a subset of variables extracted from the large datasets is one of the approaches towards solving this problem.

METHODOLOGY

In this study we analyzed the transcriptome of Glossina morsitans morsitans (Tsetsefly) antennae after exposure to either a repellant (δ-nonalactone) or an attractant (ε-nonalactone). We identified 308 genes that were upregulated or downregulated due to exposure to a repellant (δ-nonalactone) or an attractant (ε-nonalactone) respectively. Weighted gene coexpression network analysis was used to cluster the genes into 12 modules and filter unconnected genes. Discretized and association rule mining was used to find association between genes thereby predicting the putative function of unannotated genes.

RESULTS AND DISCUSSION

Among the significantly expressed chemosensory genes (FDR < 0.05) in response to Ɛ-nonalactone were gustatory receptors (GrIA and Gr28b), ionotrophic receptors (Ir41a and Ir75a), odorant binding proteins (Obp99b, Obp99d, Obp59a and Obp28a) and the odorant receptor (Or67d). Several non-chemosensory genes with no assigned function in the NCBI database were co-expressed with the chemosensory genes. Exposure to a repellent (δ-nonalactone) did not show any significant change between the treatment and control samples. We generated a coexpression network with 276 edges and 130 nodes. Genes CAH3, Ahcy, Ir64a, Or67c, Ir8a and Or67a had node degree values above 11 and therefore could be regarded as the top hub genes in the network. Association rule mining showed a relation between various genes based on their appearance in the same itemsets as consequent and antecedent.

Inferring and analyzing gene regulatory networks from multi-factorial expression data: a complete and interactive suite

Background

High-throughput transcriptomic datasets are often examined to discover new actors and regulators of a biological response. To this end, graphical interfaces have been developed and allow a broad range of users to conduct standard analyses from RNA-seq data, even with little programming experience. Although existing solutions usually provide adequate procedures for normalization, exploration or differential expression, more advanced features, such as gene clustering or regulatory network inference, often miss or do not reflect current state of the art methodologies.

Results

We developed here a user interface called DIANE (Dashboard for the Inference and Analysis of Networks from Expression data) designed to harness the potential of multi-factorial expression datasets from any organisms through a precise set of methods. DIANE interactive workflow provides normalization, dimensionality reduction, differential expression and ontology enrichment. Gene clustering can be performed and explored via configurable Mixture Models, and Random Forests are used to infer gene regulatory networks. DIANE also includes a novel procedure to assess the statistical significance of regulator-target influence measures based on permutations for Random Forest importance metrics. All along the pipeline, session reports and results can be downloaded to ensure clear and reproducible analyses.

Conclusions

We demonstrate the value and the benefits of DIANE using a recently published data set describing the transcriptional response of Arabidopsis thaliana under the combination of temperature, drought and salinity perturbations. We show that DIANE can intuitively carry out informative exploration and statistical procedures with RNA-Seq data, perform model based gene expression profiles clustering and go further into gene network reconstruction, providing relevant candidate genes or signalling pathways to explore. DIANE is available as a web service (https://diane.bpmp.inrae.fr), or can be installed and locally launched as a complete R package.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-021-07659-2).

A deep learning-based model for screening and staging pneumoconiosis

This study aims to develop an artificial intelligence (AI)-based model to assist radiologists in pneumoconiosis screening and staging using chest radiographs. The model, based on chest radiographs, was developed using a training cohort and validated using an independent test cohort. Every image in the training and test datasets were labeled by experienced radiologists in a double-blinded fashion. The computational model started by segmenting the lung field into six subregions. Then, convolutional neural network classification model was used to predict the opacity level for each subregion respectively. Finally, the diagnosis for each subject (normal, stage I, II, or III pneumoconiosis) was determined by summarizing the subregion-based prediction results. For the independent test cohort, pneumoconiosis screening accuracy was 0.973, with both sensitivity and specificity greater than 0.97. The accuracy for pneumoconiosis staging was 0.927, better than that achieved by two groups of radiologists (0.87 and 0.84, respectively). This study develops a deep learning-based model for screening and staging of pneumoconiosis using man-annotated chest radiographs. The model outperformed two groups of radiologists in the accuracy of pneumoconiosis staging. This pioneer work demonstrates the feasibility and efficiency of AI-assisted radiography screening and diagnosis in occupational lung diseases.

Gaussian and Mixed Graphical Models as (multi-)omics data analysis tools

Gaussian Graphical Models (GGMs) are tools to infer dependencies between biological variables. Popular applications are the reconstruction of gene, protein, and metabolite association networks. GGMs are an exploratory research tool that can be useful to discover interesting relations between genes (functional clusters) or to identify therapeutically interesting genes, but do not necessarily infer a network in the mechanistic sense. Although GGMs are well investigated from a theoretical and applied perspective, important extensions are not well known within the biological community. GGMs assume, for instance, multivariate normal distributed data. If this assumption is violated Mixed Graphical Models (MGMs) can be the better choice. In this review, we provide the theoretical foundations of GGMs, present extensions such as MGMs or multi-class GGMs, and illustrate how those methods can provide insight in biological mechanisms. We summarize several applications and present user-friendly estimation software. This article is part of a Special Issue entitled: Transcriptional Profiles and Regulatory Gene Networks edited by Dr. Dr. Federico Manuel Giorgi and Dr. Shaun Mahony.

Network inference in systems biology: recent developments, challenges, and applications

One of the most interesting, difficult, and potentially useful topics in computational biology is the inference of gene regulatory networks (GRNs) from expression data. Although researchers have been working on this topic for more than a decade and much progress has been made, it remains an unsolved problem and even the most sophisticated inference algorithms are far from perfect. In this paper, we review the latest developments in network inference, including state-of-the-art algorithms like PIDC, Phixer, and more. We also discuss unsolved computational challenges, including the optimal combination of algorithms, integration of multiple data sources, and pseudo-temporal ordering of static expression data. Lastly, we discuss some exciting applications of network inference in cancer research, and provide a list of useful software tools for researchers hoping to conduct their own network inference analyses.

Repurposing Zileuton as a Depression Drug Using an AI and In Vitro Approach

Repurposing drugs to target M1 macrophages inflammatory response in depression constitutes a bright alternative for commonly used antidepressants. Depression is a significant type of mood disorder, where patients suffer from pathological disturbances associated with a proinflammatory M1 macrophage phenotype. Presently, the most commonly used antidepressants such as Zoloft and Citalopram can reduce inflammation, but suffer from dangerous side effects without offering specificity toward macrophages. We employed a new strategy for drug repurposing based on the integration of RNA-seq analysis and text mining using deep neural networks. Our system employs a Google semantic AI universal encoder to compute sentences embedding. Sentences similarity is calculated using a sorting function to identify drug compounds. Then sentence relevance is computed using a custom-built convolution differential network. Our system highlighted the NRF2 pathway as a critical drug target to reprogram M1 macrophage response toward an anti-inflammatory profile (M2). Using our approach, we were also able to predict that lipoxygenase inhibitor drug zileuton could modulate NRF2 pathway in vitro. Taken together, our results indicate that reorienting zileuton usage to modulate M1 macrophages could be a novel and safer therapeutic option for treating depression.

Key challenges facing data-driven multicellular systems biology

Increasingly sophisticated experiments, coupled with large-scale computational models, have the potential to systematically test biological hypotheses to drive our understanding of multicellular systems. In this short review, we explore key challenges that must be overcome to achieve robust, repeatable data-driven multicellular systems biology. If these challenges can be solved, we can grow beyond the current state of isolated tools and datasets to a community-driven ecosystem of interoperable data, software utilities, and computational modeling platforms. Progress is within our grasp, but it will take community (and financial) commitment.