Skip to main content

Support for phosphoinositol 3 kinase and mTOR inhibitors as treatment for lupus using in-silico drug-repurposing analysis



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 [1], the biological mechanisms behind disease development are poorly understood in general and it has been proposed that genetic and environmental factors are involved [2]. 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 [3]. Nevertheless, the multifactorial nature and the undefined etiology of this disease contribute to the absence of efficient treatments [4].

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 [5].

The Connectivity Map [6] 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 [7].

In this context, Lincscloud [8] 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. [9] explored the transcriptional connections between cAMP signaling and GPCR pathway-associated drug resistance candidates. Santagata et al. [10] revealed a strong connection between the HSF1 gene and compounds that inhibit protein translation, while Siavelis et al. [11] 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 [12] 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.

Drug-repurposing 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.

Fig. 1

Integrative drug-repurposing analysis. Fourteen signatures of SLE were obtained from 14 different datasets. Each signature was queried on the Lincscloud database and a set of drugs and knock-down and knock-in genes was obtained with similarity scores. The median similarity score and empirical p values were calculated to select significant results across all datasets. Bottom: summary interpretation of the positively and negatively correlated results. NCBI GEO National Center for Biotechnology Information Gene Expression Omnibus, SLE systemic lupus erythematosus

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 [13], ChEBI [14], and Therapeutic Target Database [15]. 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.

Fig. 2

Heatmap of significant drugs representing similarity scores for each drug in the results of each dataset. Rows: results of the different datasets used for the analysis. Datasets classified according to the blood cell type (see key). Columns: different drugs sorted decreasingly by the median of similarity scores, from left to right (Color figure online)

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.

Table 1 Drugs obtained and their biological targets

To further explore this result, we used information from the KEGG database [16] 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 [17], all important and impaired in SLE. We also obtained dual inhibitors of PI3K and mTOR such as NVP-BEZ235 [18]. 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 [19], or inhibitors of the mitogen-activated protein kinase (MAPK) signaling pathway [20].

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 [21], 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.

Table 2 Significant knock-down and knock-in genes obtained


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 [24], 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 [27]. In addition, Lohman et al. [28] 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 [29].

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 [30]. 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 [33] through impaired PKC activation. This pathway has also been proposed as a potential therapeutic target for rheumatoid arthritis [34]. 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 [3739]. 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 [40]. Indeed, PI3K and mTOR act in the same signaling cascade [38] promoting the interferon and cytokine signaling pathways [41].

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 [3739], 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.



Cyclin-dependent kinase


Histone deacetylase




Mitogen-activated protein kinase


Mammalian target of Rapamycin


National Centre for Biotechnology Information Gene Expression Omnibus database


Phosphoinositol 3 kinase


Protein kinase C


Systemic lupus erythematosus


Regulatory T cell


  1. 1.

    Crow MK. Type I, interferon in the pathogenesis of lupus. J Immunol. 2014;192:5459–68.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  2. 2.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Chambers SA, Rahman A, Isenberg DA. Treatment adherence and clinical outcome in systemic lupus erythematosus. Rheumatology. 2007;46:895–8.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    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.

    CAS  Article  Google Scholar 

  5. 5.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    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.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    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.

  8. 8.

    LINCS Consortium. Lincscloud.

  9. 9.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    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.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    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.

  12. 12.

    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.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    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.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    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.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    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.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    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.

    CAS  Article  PubMed Central  Google Scholar 

  18. 18.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    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.

    Article  PubMed  Google Scholar 

  20. 20.

    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.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. 21.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    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.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    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.

    Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Shao Y, Gao Z, Marks PA, Jiang X. Apoptotic and autophagic cell death induced by histone deacetylase inhibitors. PNAS. 2004;101:18030–5.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    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.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    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.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    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.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Brzostek-Racine S, Gordon C, Van Scoy S, Reich NC. The DNA damage response induces interferon. J Immunol. 2011;187:5336–45.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    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.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    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.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    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.

    CAS  Article  Google Scholar 

  33. 33.

    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.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    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.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    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.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    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.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    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.

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    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.

    Google Scholar 

  40. 40.

    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.

    Article  Google Scholar 

  41. 41.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references


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

Not applicable.

Authors’ contributions

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.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information



Corresponding authors

Correspondence to Pedro Carmona-Sáez or Marta E. Alarcón-Riquelme.

Additional files

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 (, 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 ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

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).

Download citation


  • Systemic lupus erythematosus
  • Drug repurposing
  • Lincscloud
  • Gene expression
  • Autoimmunity
  • Drug discovery