DRMref Logo

Home

Download

Statistics

Landscape

Help

Contact

Center for Computational Systems Medicine
leaf

Dataset summary

leaf

Drug summary

leaf

Umap of single cell types and conditions (resistant and sensitive)

leaf

Comparison of cell composition between the resistant and sensitive groups (This dataset does not contain this module)

leaf

Comparison of intra-tumor heterogeneity (ITH) and epithelial-mesenchymal transition (EMT) scores of malignant cells between the resistant and sensitive groups

leaf

Difference of cell-cell interactions between the resistant and sensitive groups

Soverview

Differentially expressed genes (DEGs) between the resistant and sensitive groups for each cell type

Soverview

Mechanism analysis of drug resistance-related DEGs for each cell type in the pre-treatment samples (This dataset does not contain this module)

Soverview

Mechanism analysis of drug resistance-related DEGs for each cell type in the post-treatment samples

Soverview

Mechanism analysis of drug resistance-related DEGs for each cell type in one patient or cell lines (This dataset does not contain this module)

Soverview

MicroRNAs (miRNAs) regulating drug resistance-related DEGs for Malignant cells

Soverview

Motifs and transcription factors (TFs) regulating drug resistance-related DEGs for each cell type

     •

Differentially expressed genes act as transcription factors

Dataset: GSE207422_Sin

Dataset summary for GSE207422_Sin

Datast informationDatasetGSE207422_Sin
PMID36869384
Raw data ID/linkhttps://ngdc.cncb.ac.cn/gsa-human/browse/HRA001033
OrganismHomo sapiens
Sourcepatients
TissueTumor tissue
Cancer type level1Lung cancer
Cancer type level2EGFR/ALK mutation negative non-small cell lung cancer (NSCLC)
Regimensintilimab + carboplatin + (docetaxel or pemetrexed or emcitabine)
Drug typeImmunotherapy
Sample sizepost (resistant 5, sensitive 2)
Cell number22187
Extract protocolBD Rhapsody system
Data processing BD Rhapsody Whole Transcriptome Analysis (WTA) Pipeline
Public datePublic on Mar 21, 2023
DescriptionThis dataset has 7 patients with post-treatment samples.

Top

Drug summary for GSE207422_Sin

Drug nameDrugBank IDDrug typeTargets nameTargets uniprot ID
"Sintilimab"

DB15765

biotechPDCD1Q15116
"carboplatin" is not included in the drug list.
"Docetaxel"

DB01248

small moleculeTUBB1; MAP2; MAP4; MAPT; BCL2; NR1I2Q9H4B7; P11137; P27816; P10636; P10415; O75469
"Pemetrexed"

DB00642

small moleculeTYMS; ATIC; DHFR; GARTP04818; P31939; P00374; P22102
"emcitabine" is not included in the drug list.

Top

Umap of single cell types and conditions (resistant and sensitive)

check buttonUmap of cell types and conditions (resistant and sensitive). If this dataset has both pre-treatment and post-treatment samples, the cell types were annotated together. If this dataset is a cell line dataset, the UMAP of cell clusters and the cell ratio of drug-resistant and sensitive groups within each cluster will also be shown. (If the image exists, the user can click on it to enlarge it in a new window.)
Umap of cell typesboxplotUmap of conditionsboxplot

Top

Comparison of cell composition between the resistant and sensitive groups

check buttonIn the pre-treatment samples, the cell composition of the resistant and sensitive groups were compared.
* represents significant difference.
(If the image exists, the user can click on it to enlarge it in a new window.)

check buttonIn the post-treatment samples, the cell composition of the resistant and sensitive groups were compared.
* represents significant difference.
(If the image exists, the user can click on it to enlarge it in a new window.)

Top

Comparison of ITH and EMT scores of malignant cells between the resistant and sensitive groups

check buttonComparison of ITH (intra-tumoral heterogeneity) and EMT (epithelial-mesenchymal transition) scores of malignant cells between the resistant and sensitive groups. ITH represents intratumoral heterogeneity, EMT represents epithelial mesenchymal transition, pre represents pre-treatment, post represents post-treatment, post2 represents longer post-treatment. (If the image exists, the user can click on it to enlarge it in a new window.)
Comparison of ITH scoresboxplotComparison of EMT scoresboxplot

Top

Difference of cell-cell interactions between the resistant and sensitive groups

check buttonIn the pre-treatment samples, the difference in cell-cell interactions between the resistant and sensitive groups was compared. The numbers indicate the count of cell-cell interactions that were either higher or lower in the resistant group compared to the sensitive group. The dot plots of significant ligand-receptor pairs in the resistant group and sensitive group are subsequently displayed. (If the image exists, the user can click on it to enlarge it in a new window. Complete files can be downloaded from the download section.)

check buttonIn the post-treatment samples, the difference in cell-cell interactions between the resistant and sensitive groups was compared. The numbers indicate the count of cell-cell interactions that were either higher or lower in the resistant group compared to the sensitive group. The dot plots of significant ligand-receptor pairs in the resistant group and sensitive group are subsequently displayed. (If the image exists, the user can click on it to enlarge it in a new window. Complete files can be downloaded from the download section.)

resistant vs. sensitive

boxplot
Dot plot of significant ligand-receptor pairs in the resistant groupDot plot of significant ligand-receptor pairs in the sensitive group
boxplotboxplot

check buttonThis dataset only has one patient or is a cell line dataset, the resistant sample is the post-treatment sample and the sensitive sample is the pre-treatment sample. The difference in cell-cell interactions between the resistant and sensitive groups was compared. The numbers indicate the count of cell-cell interactions that were either higher or lower in the resistant group compared to the sensitive group. The dot plots of significant ligand-receptor pairs in the resistant group and sensitive group are subsequently displayed. (If the image exists, the user can click on it to enlarge it in a new window. Complete files can be downloaded from the download section.)

Top

Differentially expressed genes (DEGs) between the resistant and sensitive groups for each cell type

check buttonIn the pre-treatment samples, the differentially expressed genes (|avg_log2FC|>0.25 & p_val_adj<0.05) between the resistant and sensitive groups for each cell type.

check buttonIn the post-treatment samples, the differentially expressed genes (|avg_log2FC|>0.25 & p_val_adj<0.05) between the resistant and sensitive groups for each cell type.
Gene symbolGene idavg_log2FCp_valp_val_adjpct.1pct.2Cell typeTimepoint

FCGR3B

ENSG00000162747-0.792630.00e+000.00e+000.9090.985Neutrophilspost

JAML

ENSG00000160593-0.803680.00e+000.00e+000.2840.673Neutrophilspost

MYADM

ENSG00000179820-0.8176410.00e+000.00e+000.3340.667Neutrophilspost

LYZ

ENSG00000090382-0.9217510.00e+002.15e-420.3520.645Neutrophilspost

IER2

ENSG00000160888-0.9667070.00e+000.00e+000.5420.86Neutrophilspost

BTG2

ENSG00000159388-1.013990.00e+000.00e+000.7080.915Neutrophilspost

PROK2

ENSG00000163421-1.050641.51e-363.62e-320.1770.393Neutrophilspost

SCGB1A1

ENSG00000149021-1.14090.00e+000.00e+000.0590.537Neutrophilspost

NCOA4

ENSG00000266412-1.144350.00e+000.00e+000.2590.709Neutrophilspost

C5AR1

ENSG00000197405-1.205870.00e+000.00e+000.7160.936Neutrophilspost

RHOB

ENSG00000143878-1.251530.00e+000.00e+000.3110.675Neutrophilspost

MNDA

ENSG00000163563-1.324140.00e+000.00e+000.6530.943Neutrophilspost

KLF2

ENSG00000127528-1.328220.00e+000.00e+000.4520.866Neutrophilspost

RGS2

ENSG00000116741-1.341170.00e+000.00e+000.6920.958Neutrophilspost

PTGS2

ENSG00000073756-1.81880.00e+000.00e+000.2750.752Neutrophilspost

SPP1

ENSG000001187852.849575.04e-101.21e-050.3060.082Mono/Macropost

CTSL

ENSG000001350471.942573.86e-409.28e-360.8950.61Mono/Macropost

CTSB

ENSG000002851321.58150.00e+000.00e+000.9970.943Mono/Macropost

SDS

ENSG000001350941.559551.48e-293.54e-250.5310.057Mono/Macropost

CD163

ENSG000001775751.549460.00e+000.00e+000.8660.39Mono/Macropost

LGMN

ENSG000001006001.447912.02e-284.86e-240.6860.283Mono/Macropost

GLUL

ENSG000001358211.419770.00e+000.00e+000.9720.849Mono/Macropost

TMEM176B

ENSG000001065651.39992.19e-405.26e-360.8320.447Mono/Macropost

CCL2

ENSG000001086911.295321.95e-104.68e-060.3420.094Mono/Macropost

HIF1A

ENSG000001006441.23644.84e-321.16e-270.7080.277Mono/Macropost

SOD2

ENSG000002912371.217948.50e-242.04e-190.7640.472Mono/Macropost

SAT1

ENSG000001300661.158041.26e-442.87e-400.9840.805Mono/Macropost

TMEM176A

ENSG000000029331.143427.32e-351.76e-300.7550.314Mono/Macropost

PLIN2

ENSG000001478721.136423.89e-209.34e-160.810.541Mono/Macropost

RNASE1

ENSG000001295381.134471.43e-083.43e-040.2770.075Mono/Macropost

HAMP

ENSG000001056971.12243.51e-138.42e-090.3220.05Mono/Macropost

PLTP

ENSG000001009791.108672.52e-196.06e-150.5330.17Mono/Macropost

MT2A

ENSG000001251481.093273.00e-127.21e-080.5930.327Mono/Macropost

LILRB4

ENSG000002785551.040284.21e-381.01e-330.7530.264Mono/Macropost

IFI6

ENSG000001267091.038256.90e-211.66e-160.8710.642Mono/Macropost

PFKFB3

ENSG000001705251.031564.58e-301.10e-250.570.088Mono/Macropost

IFITM3

ENSG000001420890.9954162.16e-275.19e-230.7680.371Mono/Macropost

ADM

ENSG000001489260.9825512.08e-195.00e-150.4140.05Mono/Macropost

TIMP1

ENSG000001022650.9490884.33e-131.04e-080.5170.214Mono/Macropost

PSAP

ENSG000001977460.9315952.16e-405.20e-360.9970.956Mono/Macropost

PLAUR

ENSG000000114220.9040143.56e-188.55e-140.8410.629Mono/Macropost

VEGFA

ENSG000001127150.9028641.31e-223.15e-180.5110.088Mono/Macropost

CALR

ENSG000001792180.8977577.62e-291.83e-240.9480.755Mono/Macropost

IER3

ENSG000002350300.8920241.40e-133.37e-090.6350.327Mono/Macropost

ZFP36L1

ENSG000001856500.8666871.89e-174.54e-130.7930.503Mono/Macropost

G0S2

ENSG000001236890.8590132.04e-064.91e-020.3370.157Mono/Macropost

DDIT4

ENSG000001682090.857068.66e-182.08e-130.40.057Mono/Macropost

ARL4C

ENSG000001880420.8468792.23e-175.36e-130.5640.226Mono/Macropost

BASP1

ENSG000001767880.8367941.04e-212.51e-170.6690.314Mono/Macropost

CREG1

ENSG000001431620.8346381.05e-252.51e-210.8490.673Mono/Macropost
Page: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37

check buttonThis dataset only has one patient or is a cell line dataset, the resistant sample is the post-treatment sample and the sensitive sample is the pre-treatment sample. The following are differentially expressed genes (|avg_log2FC|>0.25 & p_val_adj<0.05) between the resistant and sensitive groups for each cell type.

Top

Mechanism analysis of drug resistance-related DEGs for each cell type in the pre-treatment samples

check buttonEnrichment results for 5 known drug-resistant mechanisms in malignant cells and 1 known mechanism in TME cells.The "Signature" column represents the differentially expressed genes (DEGs) used for enrichment. The "Geneset" column represents the genes included in this mechanism. The "Overlap" column represents the number of genes in the "Geneset" that are enriched in the "Signature". The "Enriched Genes" column represents the genes in the "Geneset" that are enriched in the "Signature". (If the image exists, the user can click on it to enlarge it in a new window.)

check buttonHallmark, KEGG and GOBP enrichment results for up-regulated genes in resistant group of pre-treatment samples. Only show the top 50 pathways with fdr<0.05. If certain cell types are not displayed, it means that there are no significant pathways based on DEGs in those cell types. (If the image exists, the user can click on it to enlarge it in a new window. Complete files can be downloaded from the download section.)

check buttonHallmark, KEGG and GOBP enrichment results for down-regulated genes in resistant group of pre-treatment samples. Only show the top 50 pathways with fdr<0.05. If certain cell types are not displayed, it means that there are no significant pathways based on DEGs in those cell types. (If the image exists, the user can click on it to enlarge it in a new window.Complete files can be downloaded from the download section.)

Top

Mechanism analysis of drug resistance-related DEGs for each cell type in the post-treatment samples

check buttonEnrichment results for 5 known drug resistance mechanisms in malignant cells and 1 known mechanism in TME cells.The "Signature" column represents the differentially expressed genes (DEGs) used for enrichment. The "Geneset" column represents the genes included in this mechanism. The "Overlap" column represents the number of genes in the "Geneset" that are enriched in the "Signature". The "Enriched Genes" column represents the genes in the "Geneset" that are enriched in the "Signature". (If the image exists, the user can click on it to enlarge it in a new window.)
Enrichment results for 5 known drug resistance mechanisms in malignant cells
boxplot
MechnismPvalFDRSignatureGenesetOverlapBackgroundCell typeEnriched Genes
Unusual Activation of Pro-survival Pathway1.10e-035.30e-031997811623467Malignant cellsCDKN1A,GSTP1,SLC2A1,HSPA1A,ID1,YWHAZ,FOXA1,GAPDH,LDHA,LRIG1,NOTCH3,NQO1,RAC1,UCP2,ENO1,ALDH1A1
Drug Inactivation by Structure Modification2.10e-014.50e-0119927123467Malignant cellsGSTP1
Irregularity in Drug Uptake and Drug Efflux2.70e-014.50e-0119937123467Malignant cellsSLC2A1
Epigenetic Alteration of DNA, RNA or Protein9.80e-011.00e+00199455123467Malignant cellsGAPDH
Aberration of the Drug's Therapeutic Target1.00e+001.00e+0019990023467Malignant cells

Enrichment results for 1 known drug resistance mechanisms in TME cells

boxplot
MechnismPvalFDRSignatureGenesetOverlapBackgroundCell typeEnriched Genes
Regulation by the Disease Microenvironment3.90e-013.90e-01168186223467B cellsTP53INP1,CXCR4
Regulation by the Disease Microenvironment2.00e-012.00e-01194186323467CD4+ T cellsCXCR4,FYN,SOCS3
Regulation by the Disease Microenvironment3.00e-013.00e-01349186423467CD8+ T cellsCXCR4,RHOB,SOCS3,MALAT1
Regulation by the Disease Microenvironment1.00e+001.00e+007186023467FibroblastsNA
Regulation by the Disease Microenvironment1.00e+001.00e+0041186023467Mast cellsNA
Regulation by the Disease Microenvironment3.00e-013.00e-01456186523467Mono_MacroVEGFA,CXCR4,SOCS3,VIM,FBP1
Regulation by the Disease Microenvironment3.50e-013.50e-01266186323467NeutrophilsPAK1,CXCR4,RHOB
Regulation by the Disease Microenvironment1.90e-021.90e-0227186223467NK cellsCX3CR1,CXCR4
Regulation by the Disease Microenvironment9.20e-029.20e-0264186223467pDCsCX3CR1,CXCR4
Regulation by the Disease Microenvironment1.00e+001.00e+0058186023467Plasma cellsNA

check buttonHallmark, KEGG and GOBP enrichment results for up-regulated genes in resistant group of post-treatment samples. Only show the top 50 pathways with fdr<0.05. If certain cell types are not displayed, it means that there are no significant pathways based on DEGs in those cell types. (If the image exists, the user can click on it to enlarge it in a new window. Complete files can be downloaded from the download section.)
Up-regulated HallmarkUp-regulated KEGGUp-regulated GO BP
boxplotboxplotboxplot

check buttonHallmark, KEGG and GOBP enrichment results for down-regulated genes in resistant group of post-treatment samples. Only show the top 50 pathways with fdr<0.05. If certain cell types are not displayed, it means that there are no significant pathways based on DEGs in those cell types. (If the image exists, the user can click on it to enlarge it in a new window. Complete files can be downloaded from the download section.)
Down-regulated HallmarkDown-regulated KEGGDown-regulated GO BP
There is no signaficant down-regulated pathways in resistant groupboxplotboxplot

Top

Mechanism analysis of drug resistance-related DEGs for each cell type in one patient or cell lines

check buttonThis dataset only has one patient or is a cell line dataset, the resistant sample is the post-treatment sample and the sensitive sample is the pre-treatment sample. Enrichment results for 5 known drug resistance mechanisms in malignant cells and 1 known mechanism in TME cells are shown here.The "Signature" column represents the differentially expressed genes (DEGs) used for enrichment. The "Geneset" column represents the genes included in this mechanism. The "Overlap" column represents the number of genes in the "Geneset" that are enriched in the "Signature". The "Enriched Genes" column represents the genes in the "Geneset" that are enriched in the "Signature". (If the image exists, the user can click on it to enlarge it in a new window.)

check buttonHallmark, KEGG and GOBP enrichment results for up-regulated genes in resistant group. Only show the top 50 pathways with fdr<0.05. If certain cell types are not displayed, it means that there are no significant pathways based on DEGs in those cell types. (If the image exists, the user can click on it to enlarge it in a new window. Complete files can be downloaded from the download section.)

check buttonHallmark, KEGG and GOBP enrichment results for down-regulated genes in resistant group. Only show the top 50 pathways with fdr<0.05. If certain cell types are not displayed, it means that there are no significant pathways based on DEGs in those cell types. (If the image exists, the user can click on it to enlarge it in a new window. Complete files can be downloaded from the download section.)

Top

MicroRNAs (miRNAs) regulating drug resistance-related DEGs for malignant cells

check buttonIn the pre-treatment samples, the miRNAs with a score greater than 80 that regulate the top10 abs(avg_log2FC) drug resistance-related DEGs for malignant cells are shown here. If the gene does not have any miRNAs with a score greater than 80, it will not be shown. Here only the names of miRNAs with a score greater than 90 are shown, but some names will not be shown because of too many overlaps. (Complete files containing all datasets and cell types can be downloaded from the download section.)

check buttonIn the post-treatment samples, the miRNAs with a score greater than 80 that regulate the top10 abs(avg_log2FC) drug resistance-related DEGs for Malignant cells are shown here. If the gene does not have any miRNAs with a score greater than 80, it will not be shown. Here only the names of miRNAs with a score greater than 90 are shown, but some names will not be shown because of too many overlaps. (Complete files containing all datasets and cell types can be downloaded from the download section.)
Regulating up-regulated DEGs in resistant groupRegulating down-regulated DEGs in resistant group
boxplotboxplot

check buttonThis dataset only has one patient or is a cell line dataset, the resistant sample is the post-treatment sample and the sensitive sample is the pre-treatment sample. The miRNAs with a score greater than 80 that regulate the top10 abs(avg_log2FC) drug resistance-related DEGs for Malignant cells are shown here. Here only the names of miRNAs with a score greater than 80 are shown, but some names will not be shown because of too many overlaps. If the gene does not have any miRNAs with a score greater than 90, it will not be shown. (Complete files containing all datasets and cell types can be downloaded from the download section.)

Top

Motifs and transcription factors (TFs) regulating drug resistance-related DEGs for each cell type

check buttonIn the pre-treatment samples, the motifs and TFs that regulate drug resistance-related DEGs for each cell type are shown here. (Complete files can be downloaded from the download section.)

check buttonIn the post-treatment samples, the motifs and TFs that regulate drug resistance-related DEGs for each cell type are shown here. (Complete files can be downloaded from the download section.)
Motifs and TFs regulating up-regulated DEGs in resistant group of the post-treatment samples
motifmotifNESAUCTF_highConfCell_type
motifhocomoco__TYY2_HUMAN.H11MO.0.D4.860.0798YY2 (directAnnotation). B cells
motiftransfac_pro__M020804.550.0762ISL2 (directAnnotation). B cells
motifmetacluster_141.24.330.0737E2F2; E2F3 (directAnnotation). B cells
motiftransfac_pro__M065004.190.0722ZNF614 (directAnnotation). B cells
motiftaipale_tf_pairs__ETV2_HES7_RSCGGAANNNNNNNCACGTGNN_CAP_repr4.180.072ETV2; HES7 (directAnnotation). B cells
motiftransfac_pro__M063014.060.0706ZNF404 (directAnnotation). B cells
motiftaipale__YY2_DBD_NNCCGCCATNW_repr3.850.0683YY2 (directAnnotation). B cells
motiftaipale_tf_pairs__ELK1_TEF_NSCGGAWNTTACGTAAN_CAP3.730.0669ELK1; TEF (directAnnotation). B cells
motiftfdimers__MD000323.660.0661GFI1; GFI1B; NKX3-2 (directAnnotation). B cells
motiftaipale_cyt_meth__JUND_NRTGACGCATN_eDBD_repr3.640.0658JUND (directAnnotation). B cells
Page: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43

Motifs and TFs regulating down-regulated DEGs in resistant group of the post-treatment samples
motifmotifNESAUCTF_highConfCell_type
motiftransfac_pro__M057315.270.232ZNF660 (directAnnotation). B cells
motiftransfac_pro__M060374.570.206ZNF197 (directAnnotation). B cells
motifmetacluster_103.24.480.203ZNF136; ZNF680 (directAnnotation). B cells
motiftransfac_pro__M055424.290.196SNAPC4 (directAnnotation). B cells
motiftransfac_pro__M065394.180.192ZNF483 (directAnnotation). B cells
motifmetacluster_52.133.940.183LEF1; LEF1; LEF1; LEF1; TCF7; TCF7; TCF7L1; TCF7L1 (directAnnotation). B cells
motiftfdimers__MD003753.820.179KLF4; PITX3 (directAnnotation). B cells
motiftransfac_pro__M066613.820.179ZNF426 (directAnnotation). B cells
motiftransfac_pro__M058213.490.167ZNF880 (directAnnotation). B cells
motiftaipale__ZNF524_full_ACCCTTGAACCC_repr3.470.166ZNF524 (directAnnotation). B cells
Page: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54

check buttonThis dataset only has one patient or is a cell line dataset, the resistant sample is the post-treatment sample and the sensitive sample is the pre-treatment sample. Motifs and TFs that regulate drug resistance-related DEGs for each cell type are shown here. (Complete files can be downloaded from the download section.)

Top

Differentially expressed genes act as transcription factors

check buttonThis table shows the differentially expressed genes that act as transcription factors. (Complete files can be downloaded from the download section.)
TFmotifExpressionCell typeTimepoint

KLF2

metacluster_163.1upB cellspost

KLF2

metacluster_170.2upB cellspost

CREM

transfac_pro__M00801upCD4+ T cellspost

IKZF2

tfdimers__MD00029upCD4+ T cellspost

TCF7

metacluster_52.13downCD4+ T cellspost

JUND

taipale_cyt_meth__JUND_NRTGACGCATN_eDBD_reprupCD8+ T cellspost

PRDM1

metacluster_2.9upCD8+ T cellspost

FOXP1

hocomoco__FOXP1_HUMAN.H11MO.0.AdownMast cellspost

CEBPD

taipale_tf_pairs__ETV2_CEBPD_RSCGGANNTTGCGYAAN_CAP_reprupMono/Macropost

FOSL2

metacluster_157.2upMono/Macropost
Page: 1




1"Sun X, Zhang Y, Li H, Zhou Y, Shi S, Chen Z, He X, Zhang H, Li F, Yin J, Mou M, Wang Y, Qiu Y, Zhu F. DRESIS: the first comprehensive landscape of drug resistance information. Nucleic Acids Res. 2023 Jan 6;51(D1):D1263-D1275. doi: 10.1093/nar/gkac812. PMID: 36243960; PMCID: PMC9825618."
2"Wishart DS, Feunang YD, Guo AC, Lo EJ, Marcu A, Grant JR, Sajed T, Johnson D, Li C, Sayeeda Z, Assempour N, Iynkkaran I, Liu Y, Maciejewski A, Gale N, Wilson A, Chin L, Cummings R, Le D, Pon A, Knox C, Wilson M. DrugBank 5.0: a major update to the DrugBank database for 2018. Nucleic Acids Res. 2018 Jan 4;46(D1):D1074-D1082. doi: 10.1093/nar/gkx1037. PMID: 29126136; PMCID: PMC5753335."