DeLA-Drug: A Deep Learning Algorithm for Automated Design of Druglike Analogues.

In this paper, we present a deep learning algorithm for automated design of druglike analogues (DeLA-Drug), a recurrent neural network (RNN) model composed of two long short-term memory (LSTM) layers and conceived for data-driven generation of similar-to-bioactive compounds. DeLA-Drug captures the syntax of SMILES strings of more than 1 million compounds belonging to the ChEMBL28 database and, by employing a new strategy called sampling with substitutions (SWS), generates molecules starting from a single user-defined query compound. Remarkably, the algorithm preserves druglikeness and synthetic accessibility of the known bioactive compounds present in the ChEMBL28 repository. The absence of any time-demanding fine-tuning procedure enables DeLA-Drug to perform a fast generation of focused libraries for further high-throughput screening and makes it a suitable tool for performing de novo design even in low-data regimes. To provide a concrete idea of its applicability, DeLA-Drug was applied to the cannabinoid receptor subtype 2 (CB2R), a known target involved in different pathological conditions such as cancer and neurodegeneration. DeLA-Drug, available as a free web platform (http://www.ba.ic.cnr.it/softwareic/deladrugportal/), can help medicinal chemists interested in generating analogues of compounds already available in their laboratories and, for this reason, good candidates for an easy and low-cost synthesis.

Structure-Based Prediction of hERG-Related Cardiotoxicity: A Benchmark Study.

Drug-induced blockade of the human ether-à-go-go-related gene (hERG) channel is today considered the main cause of cardiotoxicity in postmarketing surveillance. Hence, several ligand-based approaches were developed in the last years and are currently employed in the early stages of a drug discovery process for in silico cardiac safety assessment of drug candidates. Herein, we present the first structure-based classifiers able to discern hERG binders from nonbinders. LASSO regularized support vector machines were applied to integrate docking scores and protein-ligand interaction fingerprints. A total of 396 models were trained and validated based on: (i) high-quality experimental bioactivity information returned by 8337 curated compounds extracted from ChEMBL (version 25) and (ii) structural predictor data. Molecular docking simulations were performed using GLIDE and GOLD software programs and four different hERG structural models, namely, the recently published structures obtained by cryoelectron microscopy (PDB codes: 5VA1 and 7CN1) and two published homology models selected for comparison. Interestingly, some classifiers return performances comparable to ligand-based models in terms of area under the ROC curve (AUCMAX = 0.86 ± 0.01) and negative predictive values (NPVMAX = 0.81 ± 0.01), thus putting forward the herein proposed computational workflow as a valuable tool for predicting hERG-related cardiotoxicity without the limitations of ligand-based models, typically affected by low interpretability and a limited applicability domain. From a methodological point of view, our study represents the first example of a successful integration of docking scores and protein-ligand interaction fingerprints (IFs) through a support vector machine (SVM) LASSO regularized strategy. Finally, the study highlights the importance of using hERG structural models accounting for ligand-induced fit effects and allowed us to select the best-performing protein conformation (made available in the Supporting Information, SI) to be employed for a reliable structure-based prediction of hERG-related cardiotoxicity.

Combined microRNA and mRNA expression analysis in pediatric multiple sclerosis: an integrated approach to uncover novel pathogenic mechanisms of the disease.

Multiple sclerosis (MS) is a complex disease of the CNS that usually affects young adults, although 3-5% of cases are diagnosed in childhood and adolescence (hence called pediatric MS, PedMS). Genetic predisposition, among other factors, seems to contribute to the risk of the onset, in pediatric as in adult ages, but few studies have investigated the genetic 'environmentally naïve' load of PedMS. The main goal of this study was to identify circulating markers (miRNAs), target genes (mRNAs) and functional pathways associated with PedMS; we also verified the impact of miRNAs on clinical features, i.e. disability and cognitive performances. The investigation was performed in 19 PedMS and 20 pediatric controls (PCs) using a High-Throughput Next-generation Sequencing (HT-NGS) approach followed by an integrated bioinformatics/biostatistics analysis. Twelve miRNAs were significantly upregulated (let-7a-5p, let-7b-5p, miR-25-3p, miR-125a-5p, miR-942-5p, miR-221-3p, miR-652-3p, miR-182-5p, miR-185-5p, miR-181a-5p, miR-320a, miR-99b-5p) and 1 miRNA was downregulated (miR-148b-3p) in PedMS compared with PCs. The interactions between the significant miRNAs and their targets uncovered predicted genes (i.e. TNFSF13B, TLR2, BACH2, KLF4) related to immunological functions, as well as genes involved in autophagy-related processes (i.e. ATG16L1, SORT1, LAMP2) and ATPase activity (i.e. ABCA1, GPX3). No significant molecular profiles were associated with any PedMS demographic/clinical features. Both miRNAs and mRNA expressions predicted the phenotypes (PedMS-PC) with an accuracy of 92% and 91%, respectively. In our view, this original strategy of contemporary miRNA/mRNA analysis may help to shed light in the genetic background of the disease, suggesting further molecular investigations in novel pathogenic mechanisms.

Functional Implications of MicroRNAs in Crohn's Disease Revealed by Integrating MicroRNA and Messenger RNA Expression Profiling.

Crohn's disease (CD) is a debilitating inflammatory bowel disease (IBD) that emerges due to the influence of genetic and environmental factors. microRNAs (miRNAs) have been identified in the tissue and sera of IBD patients and may play an important role in the induction of IBD. Our study aimed to identify differentially expressed miRNAs and miRNAs with the ability to alter transcriptome activity by comparing inflamed tissue samples with their non-inflamed counterparts. We studied changes in miRNA-mRNA interactions associated with CD by examining their differential co-expression relative to normal mucosa from the same patients. Correlation changes between the two conditions were incorporated into scores of predefined gene sets to identify biological processes with altered miRNA-mediated control. Our study identified 28 miRNAs differentially expressed (p-values < 0.01), of which 14 are up-regulated. Notably, our differential co-expression analysis highlights microRNAs (i.e., miR-4284, miR-3194 and miR-21) that have known functional interactions with key mechanisms implicated in IBD. Most of these miRNAs cannot be detected by differential expression analysis that do not take into account miRNA-mRNA interactions. The identification of differential miRNA-mRNA co-expression patterns will facilitate the investigation of the miRNA-mediated molecular mechanisms underlying CD pathogenesis and could suggest novel drug targets for validation.

Meta-Analysis of Differential Connectivity in Gene Co-Expression Networks in Multiple Sclerosis.

Differential gene expression analyses to investigate multiple sclerosis (MS) molecular pathogenesis cannot detect genes harboring genetic and/or epigenetic modifications that change the gene functions without affecting their expression. Differential co-expression network approaches may capture changes in functional interactions resulting from these alterations. We re-analyzed 595 mRNA arrays from publicly available datasets by studying changes in gene co-expression networks in MS and in response to interferon (IFN)-ß treatment. Interestingly, MS networks show a reduced connectivity relative to the healthy condition, and the treatment activates the transcription of genes and increases their connectivity in MS patients. Importantly, the analysis of changes in gene connectivity in MS patients provides new evidence of association for genes already implicated in MS by single-nucleotide polymorphism studies and that do not show differential expression. This is the case of amiloride-sensitive cation channel 1 neuronal (ACCN1) that shows a reduced number of interacting partners in MS networks, and it is known for its role in synaptic transmission and central nervous system (CNS) development. Furthermore, our study confirms a deregulation of the vitamin D system: among the transcription factors that potentially regulate the deregulated genes, we find TCF3 and SP1 that are both involved in vitamin D3-induced p27Kip1 expression. Unveiling differential network properties allows us to gain systems-level insights into disease mechanisms and may suggest putative targets for the treatment.

Transcriptome analysis reveals genes associated with stem cell activation by physical exercise in the dentate gyrus of aged p16Ink4a knockout mice.

Throughout adulthood neural stem cells divide in neurogenic niches-the dentate gyrus of the hippocampus and the subventricular zone-producing progenitor cells and new neurons. Stem cells self-renew, thus preserving their pool. Furthermore, the number of stem/progenitor cells in the neurogenic niches decreases with age. We have previously demonstrated that the cyclin-dependent kinase inhibitor p16Ink4a maintains, in aged mice, the pool of dentate gyrus stem cells by preventing their activation after a neurogenic stimulus such as exercise (running). We showed that, although p16Ink4a ablation by itself does not activate stem/progenitor cells, exercise strongly induced stem cell proliferation in p16Ink4a knockout dentate gyrus, but not in wild-type. As p16Ink4a regulates stem cell self-renewal during aging, we sought to profile the dentate gyrus transcriptome from p16Ink4a wild-type and knockout aged mice, either sedentary or running for 12 days. By pairwise comparisons of differentially expressed genes and by correlative analyses through the DESeq2 software, we identified genes regulated by p16Ink4a deletion, either without stimulus (running) added, or following running. The p16Ink4a knockout basic gene signature, i.e., in sedentary mice, involves upregulation of apoptotic, neuroinflammation- and synaptic activity-associated genes, suggesting a reactive cellular state. Conversely, another set of 106 genes we identified, whose differential expression specifically reflects the pattern of proliferative response of p16 knockout stem cells to running, are involved in processes that regulate stem cell activation, such as synaptic function, neurotransmitter metabolism, stem cell proliferation control, and reactive oxygen species level regulation. Moreover, we analyzed the regulation of these stem cell-specific genes after a second running stimulus. Surprisingly, the second running neither activated stem cell proliferation in the p16Ink4a knockout dentate gyrus nor changed the expression of these genes, confirming that they are correlated to the stem cell reactivity to stimulus, a process where they may play a role regulating stem cell activation.

Transcriptome Analysis in a Mouse Model of Premature Aging of Dentate Gyrus: Rescue of Alpha-Synuclein Deficit by Virus-Driven Expression or by Running Restores the Defective Neurogenesis.

The dentate gyrus of the hippocampus and the subventricular zone are neurogenic niches where neural stem and progenitor cells replicate throughout life to generate new neurons. The Btg1 gene maintains the stem cells of the neurogenic niches in quiescence. The deletion of Btg1 leads to an early transient increase of stem/progenitor cells division, followed, however, by a decrease during adulthood of their proliferative capability, accompanied by apoptosis. Since a physiological decrease of neurogenesis occurs during aging, the Btg1 knockout mouse may represent a model of neural aging. We have previously observed that the defective neurogenesis of the Btg1 knockout model is rescued by the powerful neurogenic stimulus of physical exercise (running). To identify genes responsible for stem and progenitor cells maintenance, we sought here to find genes underlying this premature neural aging, and whose deregulated expression could be rescued by running. Through RNA sequencing we analyzed the transcriptomic profiles of the dentate gyrus isolated from Btg1 wild-type or Btg1 knockout adult (2-month-old) mice submitted to physical exercise or sedentary. In Btg1 knockout mice, 545 genes were deregulated, relative to wild-type, while 2081 genes were deregulated by running. We identified 42 genes whose expression was not only down-regulated in the dentate gyrus of Btg1 knockout, but was also counter-regulated to control levels by running in Btg1 knockout mice, vs. sedentary. Among these 42 counter-regulated genes, alpha-synuclein (Snca), Fos, Arc and Npas4 showed significantly greater differential regulation. These genes control neural proliferation, apoptosis, plasticity and memory and are involved in aging. In particular, Snca expression decreases during aging. We tested, therefore, whether an Snca-expressing lentivirus, by rescuing the defective Snca levels in the dentate gyrus of Btg1 knockout mice, could also reverse the aging phenotype, in particular the defective neurogenesis. We found that the exogenous expression of Snca reversed the Btg1 knockout-dependent decrease of stem cell proliferation as well as the increase of progenitor cell apoptosis. This indicates that Snca has a functional role in the process of neural aging observed in this model, and also suggests that Snca acts as a positive regulator of stem cell maintenance.

Cannabinoid Receptor Subtype 2 (CB2R) in a Multitarget Approach: Perspective of an Innovative Strategy in Cancer and Neurodegeneration.

The cannabinoid receptor subtype 2 (CB2R) represents an interesting and new therapeutic target for its involvement in the first steps of neurodegeneration as well as in cancer onset and progression. Several studies, focused on different types of tumors, report a promising anticancer activity induced by CB2R agonists due to their ability to reduce inflammation and cell proliferation. Moreover, in neuroinflammation, the stimulation of CB2R, overexpressed in microglial cells, exerts beneficial effects in neurodegenerative disorders. With the aim to overcome current treatment limitations, new drugs can be developed by specifically modulating, together with CB2R, other targets involved in such multifactorial disorders. Building on successful case studies of already developed multitarget strategies involving CB2R, in this Perspective we aim at prompting the scientific community to consider new promising target associations involving HDACs (histone deacetylases) and σ receptors by employing modern approaches based on molecular hybridization, computational polypharmacology, and machine learning algorithms.

Comparative study of gene set enrichment methods.

BACKGROUND: The analysis of high-throughput gene expression data with respect to sets of genes rather than individual genes has many advantages. A variety of methods have been developed for assessing the enrichment of sets of genes with respect to differential expression. In this paper we provide a comparative study of four of these methods: Fisher's exact test, Gene Set Enrichment Analysis (GSEA), Random-Sets (RS), and Gene List Analysis with Prediction Accuracy (GLAPA). The first three methods use associative statistics, while the fourth uses predictive statistics. We first compare all four methods on simulated data sets to verify that Fisher's exact test is markedly worse than the other three approaches. We then validate the other three methods on seven real data sets with known genetic perturbations and then compare the methods on two cancer data sets where our a priori knowledge is limited. RESULTS: The simulation study highlights that none of the three method outperforms all others consistently. GSEA and RS are able to detect weak signals of deregulation and they perform differently when genes in a gene set are both differentially up and down regulated. GLAPA is more conservative and large differences between the two phenotypes are required to allow the method to detect differential deregulation in gene sets. This is due to the fact that the enrichment statistic in GLAPA is prediction error which is a stronger criteria than classical two sample statistic as used in RS and GSEA. This was reflected in the analysis on real data sets as GSEA and RS were seen to be significant for particular gene sets while GLAPA was not, suggesting a small effect size. We find that the rank of gene set enrichment induced by GLAPA is more similar to RS than GSEA. More importantly, the rankings of the three methods share significant overlap. CONCLUSION: The three methods considered in our study recover relevant gene sets known to be deregulated in the experimental conditions and pathologies analyzed. There are differences between the three methods and GSEA seems to be more consistent in finding enriched gene sets, although no method uniformly dominates over all data sets. Our analysis highlights the deep difference existing between associative and predictive methods for detecting enrichment and the use of both to better interpret results of pathway analysis. We close with suggestions for users of gene set methods.

Proteome Response of Staphylococcus xylosus DSM 20266T to Anaerobiosis and Nitrite Exposure.

The viability and competitiveness of Staphylococcus xylosus in meat mostly depend on the ability to adapt itself to rapid oxygen and nutrients depletion during meat fermentation. The utilization of nitrite instead of oxygen becomes a successful strategy for this strain to improve its performance in anaerobiosis; however, metabolic pathways of this strain underlying this adaptation, are partially known. The aim of this study was to provide an overview on proteomic changes of S. xylosus DSM 20266T cultured under anaerobiosis and nitrite exposure. Thus, two different cultures of this strain, supplemented or not with nitrite, were in vitro incubated in aerobiosis and anaerobiosis monitoring cell viability, pH, oxidation reduction potential and nitrite content. Protein extracts, obtained from cells, collected as nitrite content was depleted, were analyzed by 2DE/MALDI-TOF/TOF-MS. Results showed that DSM 20266T growth was significantly sustained by nitrite in anaerobiosis, whereas no differences were found in aerobiosis. Accordingly, nitrite content was depleted after 13 h only in anaerobiosis. At this time of sampling, a comparative proteomic analysis showed 45 differentially expressed proteins. Most differences were found between aerobic and anaerobic cultures without nitrite; the induction of glycolytic enzymes and glyoxylate cycle, the reduction of TCA enzymes, and acetate fermentation were found in anaerobiosis to produce ATP and maintain the cell redox balance. In anaerobic cultures the nitrite supplementation partially restored TCA cycle, and reduced the amount of glycolytic enzymes. These results were confirmed by phenotypic microarray that, for the first time, was carried out on cell previously adapted at the different growth conditions. Overall, metabolic changes were similar between aerobiosis and anaerobiosis NO2-adapted cells, whilst cells grown under anaerobiosis showed different assimilation profiles by confirming proteomic data; indeed, these latter extensively assimilated substrates addressed at both supplying glucose for glycolysis or fueling alternative pathways to TCA cycle. In conclusion, metabolic pathways underlying the ability of S. xylosus to adapt itself to oxygen starvation were revealed; the addition of nitrite allowed S. xylosus to take advantage of nitrite to this condition, restoring some metabolic pathway underlying aerobic behavior of the strain.

Biological and functional analysis of statistically significant pathways deregulated in colon cancer by using gene expression profiles.

Gene expression profiling offers a great opportunity for studying multi-factor diseases and for understanding the key role of genes in mechanisms which drive a normal cell to a cancer state. Single gene analysis is insufficient to describe the complex perturbations responsible for cancer onset, progression and invasion. A deeper understanding of the mechanisms of tumorigenesis can be reached focusing on deregulation of gene sets or pathways rather than on individual genes. We apply two known and statistically well founded methods for finding pathways and biological processes deregulated in pathological conditions by analyzing gene expression profiles. In particular, we measure the amount of deregulation and assess the statistical significance of predefined pathways belonging to a curated collection (Molecular Signature Database) in a colon cancer data set. We find that pathways strongly involved in different tumors are strictly connected with colon cancer. Moreover, our experimental results show that the study of complex diseases through pathway analysis is able to highlight genes weakly connected to the phenotype which may be difficult to detect by using classical univariate statistics. Our study shows the importance of using gene sets rather than single genes for understanding the main biological processes and pathways involved in colorectal cancer. Our analysis evidences that many of the genes involved in these pathways are strongly associated to colorectal tumorigenesis. In this new perspective, the focus shifts from finding differentially expressed genes to identifying biological processes, cellular functions and pathways perturbed in the phenotypic conditions by analyzing genes co-expressed in a given pathway as a whole, taking into account the possible interactions among them and, more importantly, the correlation of their expression with the phenotypical conditions.