- Research article
- Open Access
Support for phosphoinositol 3 kinase and mTOR inhibitors as treatment for lupus using in-silico drug-repurposing analysis
Arthritis Research & Therapy volume 19, Article number: 54 (2017)
Systemic lupus erythematosus (SLE) is an autoimmune disease with few treatment options. Current therapies are not fully effective and show highly variable responses. In this regard, large efforts have focused on developing more effective therapeutic strategies. Drug repurposing based on the comparison of gene expression signatures is an effective technique for the identification of new therapeutic approaches. Here we present a drug-repurposing exploratory analysis using gene expression signatures from SLE patients to discover potential new drug candidates and target genes.
We collected a compendium of gene expression signatures comprising peripheral blood cells and different separate blood cell types from SLE patients. The Lincscloud database was mined to link SLE signatures with drugs, gene knock-down, and knock-in expression signatures. The derived dataset was analyzed in order to identify compounds, genes, and pathways that were significantly correlated with SLE gene expression signatures.
We obtained a list of drugs that showed an inverse correlation with SLE gene expression signatures as well as a set of potential target genes and their associated biological pathways. The list includes drugs never or little studied in the context of SLE treatment, as well as recently studied compounds.
Our exploratory analysis provides evidence that phosphoinositol 3 kinase and mammalian target of rapamycin (mTOR) inhibitors could be potential therapeutic options in SLE worth further future testing.
Systemic lupus erythematosus (SLE) is an autoimmune disorder in which the immune system produces autoantibodies against its own cells and tissues leading to chronic inflammation and organ damage. Although some biological pathways are well known to be altered in lupus, such as the type I interferon (IFN) pathway , the biological mechanisms behind disease development are poorly understood in general and it has been proposed that genetic and environmental factors are involved . There are many classes of drugs commonly used for SLE treatment, such as corticosteroids, immunosuppressants, nonsteroidal anti-inflammatory drugs, or specific monoclonal antibodies directed against cell surface receptors or cytokines . Nevertheless, the multifactorial nature and the undefined etiology of this disease contribute to the absence of efficient treatments .
In the last decade, the widespread use of high-throughput technologies such as gene expression microarrays has enabled access to large collections of gene expression databases that can be exploited for a wide range of applications. In this context, in-silico drug-repurposing analysis based on gene expression data allows us to identify new therapeutic applications for drugs used in other contexts. This technique compares the disease gene expression signature against a large collection of profiles derived from different compounds, measuring the degree of similarity among them. A positive similarity score means that the compound produces a similar gene expression pattern to that of the disease. In the same way, a negative similarity score represents the opposite; that is, the overexpressed genes in the disease appear underexpressed in the drug signature and vice versa. This evidences that the effect of the drug on transcription is opposite to the effect of the disease, and it is reasonable to hypothesize that the drug might be able to reverse the disease gene expression program and the phenotype itself .
The Connectivity Map  was a pioneer tool that implemented this approach. Since its publication, many studies have proven the potential of this type of analysis to discover new treatments for different diseases such as several types of cancer, muscle atrophy, or inflammatory bowel disease, among others .
In this context, Lincscloud  has been deployed recently as the successor to the Connectivity Map. This database contains genetic profiles derived from a larger number of drugs and also includes knock-down and knock-in gene experiments, where whole gene expression profiles are measured after inhibiting or overexpressing a single gene. During the last few years there has been an increasing interest in the application of this approach for drug repurposing or target predictions. For example, Johannessen et al.  explored the transcriptional connections between cAMP signaling and GPCR pathway-associated drug resistance candidates. Santagata et al.  revealed a strong connection between the HSF1 gene and compounds that inhibit protein translation, while Siavelis et al.  proposed new treatments for Alzheimer’s disease.
In this work we performed a drug-repurposing analysis using a collection of gene expression signatures derived from previously published studies of SLE patients and gene expression signatures derived from Lincscloud. This analysis allowed us to establish a set of drug candidates that reverse the SLE signatures and a set of genetic targets, as well as new pharmacological paths in SLE.
Processing gene expression data
We mined the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database  to retrieve gene expression datasets from SLE patients. We selected experiments performed in any blood tissue, with case and healthy samples, without any treatment applied in the case of in-vitro samples, and each experiment with more than four replicates. To purposely obtain a heterogeneous dataset we searched for gene expression data from adult and juvenile SLE performed in different microarray platforms. By doing this we considered the patterns conserved across all SLE cases removing differences between SLE clinical types or microarray platform-dependent biases.
Each gene expression dataset was downloaded and processed independently using the R statistical environment. Genes with a high percentage of missing values (more than 15% across samples) were filtered out and remaining missing values were imputed using the average expression values within each group (case or control) of each dataset. We annotated probes to gene symbol identifiers, data were transformed to a logarithm scale, and the median expression value was computed for probes corresponding to the same gene. Differential expression analysis was performed between controls and cases for each dataset using the limma R package. Next we discarded genes presenting p > 0.05, and the top 500 most overexpressed and underexpressed genes were selected as the SLE genetic signature from each dataset to be used for further analysis.
For each independent SLE signature we performed a query on the Lincscloud database and retrieved the list of drugs and knock-in and knock-down genes with high similarity scores. We used as a similarity score the “best score 4” value, which is the proposed threshold in Lincscloud and is calculated as the mean connectivity score across the four cell lines in which the drug or perturbagen connected most strongly to the query.
To integrate the results from each independent SLE signature, a unique dataset was created where rows represent drugs and columns represent SLE signatures, and each entry of the matrix is the similarity score (best score 4 values) between drugs and SLE signatures. For each drug (row) we calculated the median similarity score across all SLE signatures. To evaluate whether equal or better scores could be obtained by chance, an empirical p value was calculated generating 10,000 random datasets permuting rows and columns in the original set of data. We then computed the p value as the fraction of permutations having a similarity score equal to or higher than (in absolute value) the observed score. Significant drugs were then selected if they presented p < 0.05 and showed a median similarity score > 80. The same procedure was applied to knock-in and knock-down gene expression signatures (see Fig. 1). The results obtained are therefore independent of the cell lines’ inherent gene expression patterns but are consistent with the patterns that are common to all of the SLE signatures.
Drug-target enrichment analysis
To evaluate whether some drug targets were significantly enriched in the list of obtained drugs we downloaded drug-target information from DrugBank , ChEBI , and Therapeutic Target Database . Data files from these three databases were parsed and an annotation file was created with information for 131,162 drugs (including synonymous names) and their biological targets. With this information, we associated target genes to the list of drugs in Lincscloud and our list of significant drugs. For drugs without target information in these resources we carefully revised the information available from compound manufacturer catalogs and the associated literature. Drugs without any information in the literature or in databases were discarded from the drug-target analysis.
Fisher’s exact test was applied to evaluate what target genes were statistically overrepresented in the list of significant drugs with respect to the total set of annotated drugs.
Analysis of gene expression signatures
After careful exploration we found 10 datasets of SLE in the NCBI GEO, two of which contained samples from juvenile SLE patients. Some of the datasets contained samples from different tissues, which we treated as independent datasets in our analysis. Thus, we identified 14 different tissue-specific datasets that passed the initial filters (see Additional file 1: Sheets 1 and 3). These datasets comprised a total of 327 SLE samples and 173 healthy controls. Each dataset was subjected to quality control and processed as described in Methods, generating 14 individual signatures including different blood tissues (see Additional file 1: Sheet 2).
Connections between SLE and drug gene expression signatures
Our analysis yielded 61 drugs that were significantly associated with the SLE signatures, 40 with similar gene expression patterns and 21 with opposite patterns (see Fig. 2 and Additional file 1: Sheet 4). Some of these compounds have been associated previously with SLE but some others have not been described in this context and hence could be new potential drug candidates (see Discussion). We used the information from DrugBank, ChEBI, and Therapeutic Target Database to annotate target genes for each drug and classify these compounds into groups with the same target.
The analysis of targets common across the list of drugs yielded three sets with similar gene expression signatures that showed significant p values, including topoisomerase II inhibitors, histone deacetylase (HDAC) inhibitors, and PKC activators, as well as three groups with negative scores, where we found phosphoinositol 3 kinase (PI3K) inhibitors, cyclin-dependent kinase (CDK) inhibitors, and mammalian target of rapamycin (mTOR) inhibitors (see Table 1). Five different compounds were PI3K inhibitors, providing the most significant p value in the enrichment analysis.
To further explore this result, we used information from the KEGG database  to construct a network of the PI3K signaling pathway (see Additional file 2). Interestingly, we found that most of the other drug targets, such as IGF, Rho, mTOR, or CDK, were also playing important roles in the PI3K signaling pathway. PI3K regulates important processes such as cell survival, immune proliferation, anti-apoptotic pathways of immune cells, and immune response linked to interferon signaling and cytokine signaling pathways , all important and impaired in SLE. We also obtained dual inhibitors of PI3K and mTOR such as NVP-BEZ235 . Other drug targets of the PI3K signaling pathway have been related with SLE or other SLE-like disorders, such as CDK inhibitors, recently proposed to be used for treatment of some autoimmune disorders , or inhibitors of the mitogen-activated protein kinase (MAPK) signaling pathway .
Study of gene effect-caused profiles
We obtained seven knock-in and 90 knock-down genes with a positive similarity score that produce an SLE-like profile, and 50 knock-down genes with a negative similarity score (see Table 2 and Additional file 1: Sheet 4) that reverse the SLE profile (genes up-regulated in the disease signature are down-regulated in the drug signature, and vice versa). Many genes have been already described in SLE, such as CD40 , interferon-related genes, and translation initiation factors, such as EIEF4 [22, 23]. Additional functional analyses with these genes are described in Additional file 3. Interestingly and in agreement with our previous analysis, we found that the gene expression signature associated with knock-down genes such as PI3K or IGF1R show a negative similarity score. That is, the inhibition of these genes could reverse the gene expression profile induced by SLE. This is consistent with the fact that gene expression profiles of drugs which inhibit these genes showed a negative score with respect to the SLE signatures.
In this study we performed a systematic screening for drugs or genes that induced similar or opposite gene expression programs to signatures from SLE patients. We integrated signatures from different blood cell populations and SLE subtypes in order to identify consistent and conserved profiles, reducing considerably the false positive ratio. In this analysis, we found 40 drugs (see Additional file 1: Sheet 4) with a positive similarity score, which induces changes similar to the SLE phenotype. In this set of compounds, HDAC, topoisomerase II, and PKC were the more significant targets. Many of these drug targets are key factors in biological processes that are altered in SLE. For example, HDAC inhibitors have been related to impairment of immune processes described in lupus, such as autophagy , although there is contradictory information about the role of some isoforms of HDAC in the immune system [25, 26]. A recent study shows that HDAC inhibitors may be suitable for treatment of autoimmunity, while primary responses to the same inhibitors were greatly impaired, probably explaining the contradiction between the positive similarity score we obtained and the potential use of HDAC inhibitors in SLE . In addition, Lohman et al.  showed that HDAC inhibitors have anti-inflammatory activity which is inversely correlated with dose, amplifying the production of inflammatory mediators at concentration > 3 μM. In another context, treatment of human cells with topoisomerase II inhibitors such as etoposide has been shown to induce interferon-stimulated genes .
Other positively correlated compounds are phorbol-12-myristate-13-acetate and ingenol, the former of which has been used to stimulate the immune response and the interferon signaling pathway . These drugs are protein kinase C (PKC) activators, a protein with some isoforms associated with SLE. In this context, the use of PKC inhibitors has been proposed as treatment for autoimmune disorders [31, 32] due to their induced increase in proliferation of regulatory T cells (Tregs). In addition, deficient MEK/ERK signaling pathway is related to SLE and cytokine generation  through impaired PKC activation. This pathway has also been proposed as a potential therapeutic target for rheumatoid arthritis . Another compound with a high positive similarity score was LE-135, which is a retinoic acid receptor inhibitor. The use of retinoic acid has been also related to an improvement in SLE recovering the Treg balance [35, 36].
Attending to drugs with negative similarity scores, we identified 21 compounds that induce opposite gene expression programs with respect to SLE signatures (see Additional file 1: Sheet 4). Almost all of them act in the same processes, down-regulating the immune response and the proliferation of immune cells. PI3K was the most significant in the target enrichment analysis, due to a set of PI3K inhibitors. PI3K inhibitors have been reported to ameliorate the effects of SLE and other autoimmune disorders in animal models [37–39]. In addition, mTOR was also found as a significantly enriched target associated with mTOR inhibitors such as NVP-BEZ235, AZD8055, TGX115, or Ku0063794.
Recent experimental evidence suggests that mTOR inhibitors may provide a new therapeutic strategy for the treatment of SLE patients . Indeed, PI3K and mTOR act in the same signaling cascade  promoting the interferon and cytokine signaling pathways .
Complementarily, the analysis of gene-caused profiles defined a set of genes – both described and not previously described in SLE – that could play an important role in the development of the disease. Some of these were interferon-related genes, transcriptional and translational factors, and a set of biological pathways related to these genes including the PI3K and the insulin signaling pathways, immune response, or transcriptional and translational processes (Additional file 2). These results are highly consistent with the analyzed list of drugs and also support that the inhibition of PI3K signaling could improve the SLE phenotype. The evidence presented here should lead not only to testing of PI3K inhibitors as potential SLE treatment, but also to actively testing any other compound obtained, such as the insulin growth factor receptor inhibitors that crosstalk with the PI3K and mTOR pathways or the Rho kinase inhibitors.
Although the Lincscloud database contains mostly experiments carried out in cancer cell lines, the integration of different SLE signatures and the inclusion of summarized drug signatures from different cell populations enable one to establish global associations based on ubiquitous expression across different cell lines. In-silico analyses are often exploratory studies and should be confirmed by in-vitro or in-vivo experiments. In this sense, previous experiments already provide evidence that PI3K inhibitors ameliorate the SLE phenotype in animal models [37–39], and that of other autoimmune disorders, although these drugs are not used clinically. Our results would therefore provide further support for the inhibition of the PI3K signaling pathway to treat SLE.
We performed an integrative in-silico drug-repurposing exploratory analysis based on comparing gene expression data of SLE against gene expression profiles produced by perturbagens from the Lincscloud database. Our analysis is designed to reduce the biases of using different microarray platforms and the heterogeneity of SLE, leading to discovery of conserved genetic patterns across different disease states or cell types. We identified a set of pathways related to biological processes impaired in SLE, compounds, and drug targets with potential therapeutic interest for SLE treatment. Based on the results, we highlighted PI3K and mTOR as good candidates and PI3K signaling pathway inhibitors as potential treatment options that are interesting enough to be further explored, although we described other targets that could also be further evaluated to test their effect in improving the phenotype of SLE, such as PKC, MAPK, or other specific kinases. This type of analysis has seldom been performed for autoimmune diseases and can provide novel therapeutic approaches for heterogeneous and multifactorial disorders, such as SLE.
Mitogen-activated protein kinase
Mammalian target of Rapamycin
- NCBI GEO:
National Centre for Biotechnology Information Gene Expression Omnibus database
Phosphoinositol 3 kinase
Protein kinase C
Systemic lupus erythematosus
Regulatory T cell
Crow MK. Type I, interferon in the pathogenesis of lupus. J Immunol. 2014;192:5459–68.
Bentham J, Morris DL, Cunninghame Graham DS, Pinder CL, Tombleson P, Behrens TW, et al. Genetic association analyses implicate aberrant regulation of innate and adaptive immunity genes in the pathogenesis of systemic lupus erythematosus. Nat Genet. 2015;47:1457–64.
Chambers SA, Rahman A, Isenberg DA. Treatment adherence and clinical outcome in systemic lupus erythematosus. Rheumatology. 2007;46:895–8.
Muangchan C, van Vollenhoven RF, Bernatsky SR, Smith CD, Hudson M, Inanç M, et al. Treatment algorithms in systemic lupus erythematosus. Arthritis Care Res (Hoboken). 2015;67:1237–45.
Iorio F, Rittman T, Ge H, Menden M, Saez-Rodriguez J. Transcriptional data: a new gateway to drug repositioning? Drug Discov Today. 2013;18:350–7.
Lamb J, Crawford ED, Peck D, Modell JW, Blat IC, Wrobel MJ, et al. The Connectivity Map: using gene-expression signatures to connect small molecules, genes, and disease. Science. 2006;313:1929–35.
Vilar S, Hripcsak G. The role of drug profiles as similarity metrics: applications to repurposing, adverse effects detection and drug-drug interactions. Brief Bioinform. 2016; bbw048.
LINCS Consortium. Lincscloud. http://lincsproject.org/.
Johannessen CM, Johnson LA, Piccioni F, Townes A, Frederick DT, Donahue MK, et al. A melanocyte lineage program confers resistance to MAP kinase pathway inhibition. Nature. 2013;504:138–42.
Santagata S, Mendillo ML, Tang Y, Subramanian A, Perley CC, Roche SP, et al. Tight coordination of protein translation and HSF1 activation supports the anabolic malignant state. Science. 2013;341:1238303.
Siavelis JC, Bourdakou MM, Athanasiadis EI, Spyrou GM, Nikita KS. Bioinformatics methods in drug repurposing for Alzheimer’s disease. Brief Bioinform. 2015;17:322–35.
Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: archive for functional genomics data sets—update. Nucleic Acids Res. 2013;41(Database issue):D991–5.
Wishart DS, Knox C, Guo AC, Shrivastava S, Hassanali M, Stothard P, et al. DrugBank: a comprehensive resource for in silico drug discovery and exploration. Nucleic Acids Res. 2006;34(Database issue):D668–72.
Hastings J, de Matos P, Dekker A, Ennis M, Harsha B, Kale N, et al. The ChEBI reference database and ontology for biologically relevant chemistry: enhancements for 2013. Nucl Acids Res. 2013;41:D456–63.
Zhu F, Shi Z, Qin C, Tao L, Liu X, Xu F, et al. Therapeutic target database update 2012: a resource for facilitating target-oriented drug discovery. Nucleic Acids Res. 2012;40(Database issue):D1128–36.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Guiducci C, Ghirelli C, Marloie-Provost M-A, Matray T, Coffman RL, Liu Y-J, et al. PI3K is critical for the nuclear translocation of IRF-7 and type I IFN production by human plasmacytoid predendritic cells in response to TLR activation. J Exp Med. 2008;205:315–22.
Mukherjee B, Tomimatsu N, Amancherla K, Camacho CV, Pichamoorthy N, Burma S. The dual PI3K/mTOR inhibitor NVP-BEZ235 is a potent inhibitor of ATM- and DNA-PKCs-mediated DNA damage responses. Neoplasia. 2012;14:34–43.
Xia Y, Lin L-Y, Liu M-L, Wang Z, Hong H-H, Guo X-G, et al. Selective inhibition of CDK7 ameliorates experimental arthritis in mice. Clin Exp Med. 2014;15:269–75.
Qu H, Bian W, Xu Y. A novel NF-κB inhibitor, DHMEQ, ameliorates pristane-induced lupus in mice. Exp Ther Med. 2014;8:100–4.
Desai-Mehta A, Lu L, Ramsey-Goldman R, Datta SK. Hyperexpression of CD40 ligand by B and T cells in human lupus and its role in pathogenic autoantibody production. J Clin Invest. 1996;97:2063–73.
Wu Y-Y, Kumar R, Haque MS, Castillejo-López C, Alarcón-Riquelme ME. BANK1 controls CpG-induced IL-6 secretion via a p38 and MNK1/2/eIF4E translation initiation pathway. J Immunol. 2013;191:6110–6.
Toro-Domínguez D, Carmona-Sáez P, Alarcón-Riquelme ME. Shared signatures between rheumatoid arthritis, systemic lupus erythematosus and Sjögren’s syndrome uncovered through gene expression meta-analysis. Arthritis Res Ther. 2014;16:489.
Shao Y, Gao Z, Marks PA, Jiang X. Apoptotic and autophagic cell death induced by histone deacetylase inhibitors. PNAS. 2004;101:18030–5.
Choi SW, Gatza E, Hou G, Sun Y, Whitfield J, Song Y, et al. Histone deacetylase inhibition regulates inflammation and enhances Tregs after allogeneic hematopoietic cell transplantation in humans. Blood. 2015;125:815–9.
Wang L, Liu Y, Han R, Beier UH, Bhatti TR, Akimova T, et al. FOXP3+ regulatory T cell development and function require histone/protein deacetylase 3. J Clin Invest. 2015;125:1111–23.
Waibel M, Christiansen AJ, Hibbs ML, Shortt J, Jones SA, Simpson I, et al. Manipulation of B-cell responses with histone deacetylase inhibitors. Nat Commun. 2015;6:6838.
Lohman R-J, Iyer A, Fairlie TJ, Cotterell A, Gupta P, Reid RC, et al. Differential anti-inflammatory activity of HDAC inhibitors in human macrophages and rat arthritis. J Pharmacol Exp Ther. 2016;356:387–96.
Brzostek-Racine S, Gordon C, Van Scoy S, Reich NC. The DNA damage response induces interferon. J Immunol. 2011;187:5336–45.
Lund ME, To J, O’Brien BA, Donnelly S. The choice of phorbol 12-myristate 13-acetate differentiation protocol influences the response of THP-1 macrophages to a pro-inflammatory stimulus. J Immunol Methods. 2016;430:64–70.
Oleksyn D, Pulvino M, Zhao J, Misra R, Vosoughi A, Jenks S, et al. Protein kinase Cβ is required for lupus development in Sle mice. Arthritis Rheum. 2013;65:1022–31.
Gray RD, Lucas CD, MacKellar A, Li F, Hiersemenzel K, Haslett C, et al. Activation of conventional protein kinase C (PKC) is critical in the generation of human neutrophil extracellular traps. J Inflamm (Lond). 2013;10:12.
Gorelik G, Fang JY, Wu A, Sawalha AH, Richardson B. Impaired T cell protein kinase Cδ activation decreases ERK pathway signaling in idiopathic and hydralazine-induced lupus. J Immunol. 2007;179:5553–63.
Thiel MJ, Schaefer CJ, Lesch ME, Mobley JL, Dudley DT, Tecle H, et al. Central role of the MEK/ERK MAP kinase pathway in a mouse model of rheumatoid arthritis: potential proinflammatory mechanisms. Arthritis Rheum. 2007;56:3347–57.
Kinoshita K, Yoo B-S, Nozaki Y, Sugiyama M, Ikoma S, Ohno M, et al. Retinoic acid reduces autoimmune renal injury and increases survival in NZB/W F1 mice. J Immunol. 2003;170:5793–8.
Xiao S, Jin H, Korn T, Liu SM, Oukka M, Lim B, et al. Retinoic acid increases Foxp3+ regulatory T cells and inhibits development of Th17 cells by enhancing TGF-β-driven Smad3 signaling and inhibiting IL-6 and IL-23 receptor expression. J Immunol. 2008;181:2277–84.
Winkler DG, Faia KL, DiNitto JP, Ali JA, White KF, Brophy EE, et al. PI3K-δ and PI3K-γ Inhibition by IPI-145 Abrogates immune responses and suppresses activity in autoimmune and inflammatory disease models. Chem Biol. 2013;20:1364–74.
Weichhart T, Säemann MD. The PI3K/Akt/mTOR pathway in innate immune cells: emerging therapeutic applications. Ann Rheum Dis. 2008;67 Suppl 3:iii70–4.
Yin CC, Trigunaite A, Dinh D, Yiu Y, Lannutti BJ, Stein PL. Targeting the phosphatidylinositol 3-kinase signaling pathway to inhibit pathogenic B cells in systemic lupus erythematosus (SLE). J Immunol. 2016;196(1 Supplement):210. 8.
Gu Z, Tan W, Ji J, Feng G, Meng Y, Da Z, et al. Rapamycin reverses the senescent phenotype and improves immuno-regulation of mesenchymal stem cells from MRL/lpr mice and systemic lupus erythematosus patients through inhibition of the mTOR signaling pathway. Aging (Albany NY). 2016;8:1102–14.
Perl A, Fernandez DR, Telarico T, Doherty E, Francis L, Phillips PE. T-cell and B-cell signaling biomarkers and treatment targets in lupus. Curr Opin Rheumatol. 2009;21:454–64.
The authors thank the Swedish Association Against Rheumatism (Reumatikerförbundet) and the Gustaf den Ve: 80th-year Foundation of Sweden for support.
Funding for the work presented received support from the Innovative Medicines Initiative Joint Undertaking (IMI-JU) under grant agreement n°115565, resources for which are composed of financial contribution from the European Union’s Seventh Framework Programme (FP7/2007-2013) and EFPIA companies’ in-kind contribution.
Availability of data and materials
PC-S and MEA-R conceived and supervised the project. PC-S and DT-D designed the study and the data analysis pipeline. DT-D implemented the drug-repurposing analysis workflow and performed all of the analysis. PC-S, MEA-R, and DT-D contributed to the interpretation of the data. All authors wrote and revised the manuscript critically. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Additional file 1:
presents information about quality control of the data, SLE signatures, and significant results obtained: Sheet 1 shows information about datasets including their GEO identifiers, the SLE state, the cell type, and the microarray platform for each, number of control and case samples, PubMed identifier, and date of publication; Sheet 2 shows results of the quality control based on the percentage of missing values and the number of significant genes in each dataset of SLE; Sheet 3 shows the signatures of each dataset used to query on Lincscloud, significant genes sorted by fold change in each signature; and Sheet 4 shows the lists of drugs and knock-in and knockdown genes with similarity scores, median of similarity scores across datasets, and significance values (XLSX 182 kb)
Additional file 2:
is a figure showing the PI3K molecular signaling pathway. Plot constructed based on the information of different PI3K interaction graphs from the KEGG database. Red, drug targets with positive similarity scores; blue, drug targets with negative similarity scores (PDF 373 kb)
Additional file 3:
shows a description of the functional analysis of gene targets obtained and their significance, and is divided into three sections: methods, results, and references. (PDF 116 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Toro-Domínguez, D., Carmona-Sáez, P. & Alarcón-Riquelme, M.E. Support for phosphoinositol 3 kinase and mTOR inhibitors as treatment for lupus using in-silico drug-repurposing analysis. Arthritis Res Ther 19, 54 (2017). https://doi.org/10.1186/s13075-017-1263-7
- Systemic lupus erythematosus
- Drug repurposing
- Gene expression
- Drug discovery