Journal of Healthcare Engineering

Journal of Healthcare Engineering / 2021 / Article

Research Article | Open Access

Volume 2021 |Article ID 6674744 | https://doi.org/10.1155/2021/6674744

Rong-Guo Yu, Jia-Yu Zhang, Zhen-Tao Liu, You-Guang Zhuo, Hai-Yang Wang, Jie Ye, Nannan Liu, Yi-Yuan Zhang, "Text Mining-Based Drug Discovery in Osteoarthritis", Journal of Healthcare Engineering, vol. 2021, Article ID 6674744, 14 pages, 2021. https://doi.org/10.1155/2021/6674744

Text Mining-Based Drug Discovery in Osteoarthritis

Academic Editor: Tiago H. Silva
Received29 Dec 2020
Revised08 Mar 2021
Accepted01 Apr 2021
Published14 Apr 2021

Abstract

Background. Osteoarthritis (OA) is a chronic and degenerative joint disease, which causes stiffness, pain, and decreased function. At the early stage of OA, nonsteroidal anti-inflammatory drugs (NSAIDs) are considered the first-line treatment. However, the efficacy and utility of available drug therapies are limited. We aim to use bioinformatics to identify potential genes and drugs associated with OA. Methods. The genes related to OA and NSAIDs therapy were determined by text mining. Then, the common genes were performed for GO, KEGG pathway analysis, and protein-protein interaction (PPI) network analysis. Using the MCODE plugin-obtained hub genes, the expression levels of hub genes were verified using quantitative real-time polymerase chain reaction (qRT-PCR). The confirmed genes were queried in the Drug Gene Interaction Database to determine potential genes and drugs. Results. The qRT-PCR result showed that the expression level of 15 genes was significantly increased in OA samples. Finally, eight potential genes were targetable to a total of 53 drugs, twenty-one of which have been employed to treat OA and 32 drugs have not yet been used in OA. Conclusions. The 15 genes (including PTGS2, NLRP3, MMP9, IL1RN, CCL2, TNF, IL10, CD40, IL6, NGF, TP53, RELA, BCL2L1, VEGFA, and NOTCH1) and 32 drugs, which have not been used in OA but approved by the FDA for other diseases, could be potential genes and drugs, respectively, to improve OA treatment. Additionally, those methods provided tremendous opportunities to facilitate drug repositioning efforts and study novel target pharmacology in the pharmaceutical industry.

1. Introduction

OA is a chronic and degenerative joint disease that occurs commonly in the elderly, and the main concomitant symptoms were cartilage loss, subchondral bone sclerosis, synovitis, and pain [1, 2]. The incidence of OA is increasing around the world due to the aging population and the growing number of obese individuals [3]. According to the statistics, the incidence of OA can reach nearly half of people over age 65 and 80% of those over age 80 [4].

There are different management methods between the end stage and early stage for OA. Traditionally, the treatments of end-stage OA for improving function and relieving pain have been joint operation [5], which may have long-term problems and potential complications, such as pain, infection, hemorrhage, and arthrofibrosis [6]. Therefore, early treatment is particularly important. Comprehensive clinical practice guidelines by the American Academy of Orthopaedic Surgeons held that NSAIDs were appropriate for all nonsurgical treatment OA patients [7]. Using the drug of NSAIDs in the treatment of OA has also been approved by the American College of Rheumatology and other academic institutions [8].

NSAIDs include traditional NSAIDs and new COX-2 inhibitors; the former blocks COX-1 and COX-2 but the new one only blocks the COX-2 enzyme. The formation of prostaglandinH2 (PGH2) by the metabolic process of the COX enzyme is metabolized by prostaglandin E (PGE) synthase into PGE2, which is one of the significant inflammation mediator [9]. COX-1, found in most tissues, is encoded by prostaglandin-endoperoxide synthase 1 (PTGS1) and COX-2, induced by various cytokines and growth factors, is encoded by prostaglandin-endoperoxide synthase 2 (PTGS2) [10]. COX-1 NSAIDs are associated with more gastrointestinal side effects than COX-2 [11], but likely have a lower rate of cardiovascular events [12]. In clinical practice, diclofenac, one kind of NSAIDs (inhibiting COX-1 and COX-2), is used as the most effective for OA treatment [4] and increases the rating of the risk of a cardiovascular event by fourfold, stroke by threefold, and all-cause death by twofold [13]. A long-term use of Ibuprofen (also inhibiting COX-1 and COX-2) shows more than three times the incidence of stroke complications compared with placebo [13]. Although there are significant limitations and health risks, approximately 65% of patients are provided oral NSAIDs to control OA [14]. After consuming NSAIDs chronically, most patients still report persistent pain and disability, and more than half will receive a prescription to stop using within one year [15]. Given all that, the securities of NSAIDs are problems that have to be paid attention to. It is still limited to research into drug therapy, and we need further studies to fill the vacancies in the development of novel therapeutic drugs.

Discovering new drug therapies is typically a 7–12 years process. With the development of computer technology and biotechnology, the article type of text mining for new drugs is becoming more common worldwide. Besides, abundant clinical trials and investment of billions of dollars are recommended but what we got are unpredictable returns on investment [16]. However, it is a less costly and possibly quicker alternative that the existing drugs are beyond their original intent to treat a novel situation [17]. Text mining of biomedical literature is a well-established approach in revealing new associations between pathologies and genes [18]. In recent years, approximately 30% of the newly Food and Drug Administration (FDA) approved vaccines and drugs are repurposing discoveries in the US [19]. As a classic example, sildenafil (Viagra), originally targeted towards angina, is used in the treatment of erectile dysfunction in men in 1998 [20].

In this study, we use bioinformatics methods to investigate novel drugs for OA. Firstly, we made a preliminary list of candidate genes for further analysis by mining text. Next, we analyzed the characteristics of genes selected by way of GO and KEGG enrichment analysis. Then, we acquired the association between those genes by PPI analysis and generated focused target genes set by the MCODE plugin. Finally, candidate key drugs were discovered from the final genes set through the drug-gene interaction analysis.

2. Materials and Methods

2.1. Text Mining

The web-based service pubmed2ensembl (http://pubmed2ensembl.ls.manchester.ac.uk/) was used to perform text mining. The link is a publicly available resource database, including biological articles from MEDLINE life science journals and online books, and linking over 2,000,000 literature to nearly 150,000 genes from 50 species [21]. All gene names, related to the search concepts found in the available literature, are extracted from pubmed2ensembl when a query is performed [21]. In this study, we performed queries using the search terms (Osteoarthritis and NSAIDs) for producing two gene lists. Common unique genes between the two lists were then used to proceed to the next steps.

2.2. Biological Process and Pathway Analysis

We used the FUNRICH (http://www.funrich.org) for a GO enrichment analysis, including biological process (BP), cellular component (CC), molecular function (MF). FUNRICH is a stand-alone software tool whose one kind of function is for the functional enrichment of genes and proteins. The KEGG pathway analysis was performed with annotations (chemokine activity, catalytic activity, transcription factor activity, receptor activity, growth factor, receptor binding, and cytokine activity) from the DAVID (Database for Annotation, Visualization, and Integrated Discovery), which contained all the gene terms we mined. DAVID (https://david.ncifcrf.gov/summary.jsp) is a web-based tool, providing functional annotation tools for the researcher to understand biological meaning among a large list of genes.

2.3. Protein Interaction Network

We used the STRING database (http://string-db.org) for acquiring a key estimation of protein-protein functional associations. The study reported that this database could comprehensively integrate the PPI of selected genes. Then, we used the MCODE (molecular complex detection) plugin of Cytoscape (version 3.7.1) to seek a gene list of the tightest interaction network among all common genes. The parameters of MCODE were used: node score cutoff ≥0.05, degree cutoff ≥2, maximum depth = 100, and k-score ≥2 [22].

2.4. Quantitative Real-Time- (qRT-) PCR Array Validation

To validate the findings of bioinformatics analysis, synovial tissues from patients were harvested for qRT-PCR validation between the OA group (n = 8) and the control group (n = 8). The use of verbal consent was approved by Fuzhou Second Hospital affiliated to Xiamen University, and oral consent was obtained by each subject. Total RNA from the lesion of synovial tissue was extracted using TRIzol reagent (TAKARA, Dalian, China). RNA samples from total RNA were reverse-transcribed to cDNA using the Revert Aid First Strand cDNA Synthesis Kit (TaKaRa, Japan) according to the manufacturer’s instructions. The conditions were set as follows: 95°C for 30 sec for denaturation, cycling stage 95 °C/10 s, and 60 °C/0.5 minutes (40 cycles). GAPDH was then handled as an internal reference. Relative mRNA expression was calculated using the 2-ΔΔCt method. values < 0.05 indicated a significant difference. The sequences of primers are presented in Table 1.


GenesSequences
ForwardReverse

PTGS2CCGAGGTGTATGTATGAGTGGGGAGGGAGAAATTAATGGG
NLRP3CACAGTCTAGTTGGGAAGACGAGACCATGTCTTCCATCAC
MMP13CTTCCTCTTCTTGAGCTGGGCAAGATACTCTACCTCTGC
MMP9ATGTCTAAGGAGGGGAGATCCTCAAACTCCTGAGCTTCAG
IL1RNCTGGTATGTTCCTGTGGAAGGGATAGGAGGCAGGACTAAT
CCL2GTAATCAAAGAGAGGGTGGGCTCATTACCTCTGCAGAGTC
CXCL8CAAGCTTCTAGGACAAGAGCCTTACCTTCACACAGAGCTG
TNFRSF1BGGGCTGTCAGACTAGAATTGCAGAGTATATGGCATGACCC
TNFGAAGACACTCAGGGAAAGAGTGATTAGAGAGAGGTCCCTG
CASP3CAGGCAGAGGGATCTTTATGCACTAGGTATGAAGGTCAGC
IL10GCTTCTTAATCCCTGGGATCCTCAGTCCTGGTCTTCTTTC
CD40CCAAAGGTCTCATAGCTAGCGTTCTCACTCCCCATAGTTC
MAPK8GAGGTCAGATGCTCTAAGTGCTATATTCCAGGCTAGGTGG
STAT3CTGAAAGACTTCATGGGAGGCCTCACTTCCCAAGTTAGAC
IL6CGTCCGTAGTTTCCTTCTAGAGGGGGAGAATACTACTCAC
NGFCTCTTGTTGAAGGGGGATAGCTCTGAGGTGCTCCTATACT
TP53GATGAGTCCTCTCTGAGTCAGTCCTAACATCCCCATCATC
RELACTATGAAGTAGCCGCTACAGGAAGTGATTAGGGACCCATC
BCL2L1TAGACCCAGACCTTCGTAAGGAAGGGAGAGAAAGAGCTTC
VEGFACAGTGCTAGGAGGAATTTCCCAGTGATCATCCTTTCCCTG
NOTCH1GGCGAGAATTATCTGTCCTCGGATCAATTACTCAAGCCCC
HMGB1CAGGAGATCAGCAGAGATTGCTGGGACCAGGATATGTAGA

2.5. Drug-Gene Interactions

The web-based Drug-Gene Interaction Database (DGIdb) (version 3.0, http://www.dgidb.org) was used to search drug-gene interactions of the confirmed genes, which might cause associations with small organic compounds or drugs. The DGIdb aggregates drug-gene interaction data from 27 kinds of sources, including PharmGKB, DrugBank, NCBI Entrez, ChEMBL, PubChem, Ensembl, and other authoritative databases. It organized the results of the search so that drug-gene interaction was available to match with relevant gene entry by its source link [23].

3. Results

3.1. Text Mining and Biological Process

Based on the strategy of data investigation described in Figure 1, from text-mining searches by pubmed2ensembl, we determined 550 unique genes (Figure 2), related to both OA and NSAIDs (we do not rule it out that some of them may have more than one name). All gene names are shown in Appendix S1. In the process of GO biological annotations, we are based on the characteristics of the results from data analysis in FUNRICH. The cellular-components annotations revealed that most of the genes are expressed in the Cytoplasm (51.1%), plasma membrane (46.6%), and extracellular (32.5%), and other details are seen in Figure 3. The GO enrichment analysis of molecular functions (transmembrane receptor protein tyrosine kinase activity, G-protein coupled receptor activity, receptor activity, growth factor activity, receptor binding, chemokine activity, protein serine/threonine kinase activity, and cytokine activity) demonstrated that 7.5% of the identified genes had G-protein coupled receptor activity (Figure 4). Biological process annotations showed cell communication (35.1%) and signal transduction (39%) were the most highly enriched terms (Figure 5).

3.2. Pathway Analysis

For the procedure of KEGG pathways enrichment analysis, more enriched genes and relatively low values showed that the pathways are particularly relevant to OA and NSAIDs. The four most enriched biological KEGG pathway were (1) pathways in cancer (5.49E-36); (2) TNF signaling pathway (2.45E-26); and (3) proteoglycans in cancer (1.86E-21), PI3K-Akt signaling pathway (6.99E-21), containing 101, 46, 56, and 74 genes related to pathways enrichment analysis, respectively, and other highly enriched pathways including FoxO signaling pathway, HIF-1 signaling pathway, MAPK signaling pathway, prolactin signaling pathway, osteoclast differentiation, apoptosis, chemokine signaling pathway, and Toll-like receptor signaling pathway. To more intuitively present the details, we conducted mapping with R statistical software (version 3.6.1) (Figure 6).

3.3. Protein-Protein Interaction

The protein-protein interaction analysis was performed using STRING and we set the following parameters: (1) minimum required interaction score: medium confidence (0.400) and (2) display simplifications: hide disconnected nodes in the network. Network Stats showed the results: the number of nodes: 452; the number of edges: 11095; average node degree: 49.1; local clustering coefficient: 0.529; expected number of edges: 4532; and PPI enrichment value: <1.0e-16 (Figure 7). Then, the data format “.tsv” was imported from the STRING EXPORTs channel to the MCODE plugin of Cytoscape and the following parameters were set: node score cutoff: 0. 05; K- -core: 2; max. depth: 100. Based on these criteria, we selected 22 central genes which formed the tightest module network, including PTGS2, NLRP3, MMP13, MMP9, IL1RN, CCL2, CXCL8, TNFRSF1B, TNF, CASP3, IL10, CD40, MAPK8, STAT3, IL6, NGF, TP53, RELA, BCL2L1, VEGFA, NOTCH1, and HMGB1 (Figure 8). We used DAVID to perform KEGG pathways enrichment analysis of central genes and picked the top 10 most enriched biological KEGG pathways to list out in Table 2 (Table 2). From Table 2, it is not difficult to find that all of those genes are involved in the enriched pathways.


PathwaysGenes in query setTotal genes in genomePercentage (%) valueGenes

TNF signaling pathway91080.2678571.73E-10CASP3, TNFRSF1B, IL6, TNF, CCL2, PTGS2, RELA, MMP9, MAPK8
NOD-like receptor signaling pathway71660.2083335.53E-09IL6, TNF, CCL2, RELA, CXCL8, MAPK8, NLRP3
Pathways in cancer115150.3273811.92E-08CASP3, IL6, PTGS2, RELA, MMP9, VEGFA, TP53, CXCL8, MAPK8, BCL2L1, STAT3
Apoptosis61350.1785715.33E-07CASP3, TNF, RELA, TP53, BCL2L1, NGF
NF-kappa B signaling pathway6930.1785712.91E-06TNF, PTGS2, RELA, CXCL8, BCL2L1, CD40
Toll-like receptor signaling pathway61020.1785717.74E-06IL6, TNF, RELA, CXCL8, MAPK8, CD40
Amyotrophic lateral sclerosis (ALS)5500.148818.84E-06CASP3, TNFRSF1B, TNF, TP53, BCL2L1
Inflammatory bowel disease (IBD)5620.148812.38E-05IL6, TNF, RELA, IL10, STAT3
Cytokine-cytokine receptor interaction72630.2083333.37E-05TNFRSF1B, IL6, TNF, CCL2, CXCL8, CD40, IL10
Adipocytokine signaling pathway5690.148813.40E-05TNFRSF1B, TNF, RELA, MAPK8, STAT3

To make sure that only the most enriched annotations were selected, a value cutoff () was set. Among the most significantly enriched pathway annotations above the cutoff, those most relevant to OA pathology based on the available literature and research were selected. A total of 10 pathways were picked out from the analysis of enriched pathway annotations.
3.4. qRT-PCR Validation of the Hub Genes

A qRT-PCR approach was used to detect the expression levels of 22 potential genes. The verification result showed that the expression levels of PTGS2, NLRP3, MMP9, IL1RN, CCL2, TNF, IL10, CD40, IL6, NGF, TP53, RELA, BCL2L1, VEGFA, and NOTCH1 were significantly increased in osteoarthritis samples () (Figure 9), which confirmed the analytical results of bioinformatics were reliable in this study.

3.5. Drug-Gene Interactions

The final confirmed genes were used to conduct drug-gene interaction analysis, and a list of 53 drugs was regarded as potentially meeting possible drug therapy for OA (Table 3). Potential gene targets of the drugs in this list are PTGS2 (35 drugs), TNF (7 drugs), VEGFA (5 drugs), MMP9 (2 drugs), and CCL2, IL1RN, IL6, and NGF (1 drug each). 21 of 53 drugs have been employed to treat OA and 32 drugs (including 16 drugs target COX-2) have not yet been used in OA. The interrelation of the 32 drugs with genes and pathways is shown in Figure 10.


NumberDrugDescriptionGeneDrug-gene interactionScoreApproved by FDAApproved use in OAClinical trial phase in O AReference (PubMed ID)

1CelecoxibCOX-2PTGS2Inhibitor12YesYes412093311|12086292
2EtodolacMainly COX-2PTGS2Inhibitor11YesYes212824918|11009046
3OxaprozinMainly COX-1PTGS2Inhibitor10YesYes219338579|9831331
4MeloxicamCOX-2PTGS2Inhibitor11YesYes410567199|10220944
5IcosapentEicosanoidsPTGS2Inhibitor7YesNoNot available12573452|16190133
6Aminosalicylic acidAntimycobacterial agentPTGS2Inhibitor6YesNoNot available16855178|12463455
7MesalamineAnti-inflammatory agentPTGS2Inhibitor9YesNoNot available16855178|12463455
8IndomethacinMainly COX-1PTGS2Inhibitor8YesYes415668944|15730717
9NabumetoneCOX-2PTGS2Inhibitor11YesYes411304699|12047490
10TenoxicamPTG synthetase inhibitorPTGS2Inhibitor10YesYes211563332|16245223
11LenalidomideImmunomodulatory drugPTGS2Modulator3YesNoNot available15598423|12720152
12RofecoxibCOX-2PTGS2Inhibitor15YesYes310580458|11398914
13PiroxicamCOX-1PTGS2Inhibitor10YesYes211952155|11785774
14SulindacPTG synthetase inhibitorPTGS2Inhibitor9YesYes211118042|10372826
15Mefenamic acidPTG synthetase inhibitorPTGS2Inhibitor11YesYes210393680|15792781
16NaproxenPTG synthetase inhibitorPTGS2Inhibitor11YesYes417604186|17612049
17SulfasalazinePTG synthetase inhibitorPTGS2Inhibitor7YesNoNot available16855178|12463455
18PhenylbutazonePTG synthetase inhibitorPTGS2Inhibitor9YesNoNot available15939622|15489888
19DiflunisalPTG synthetase inhibitorPTGS2Inhibitor11YesYesNot available11315375|8737748
20SuprofenPTG synthetase inhibitorPTGS2Inhibitor7YesNoNot available11885959|17139284
21Salicylic acidPTG synthetase inhibitorPTGS2Inhibitor10YesYesNot available15035793|15348270
22AspirinMainly COX-1PTGS2Inhibitor12YesYesNot available17522398|17181859
23BromfenacPTG synthetase inhibitorPTGS2Inhibitor6YesNoNot available16982289|16846546
24KetoprofenCOX-2PTGS2Inhibitor11YesYes314513718|11729362
25Balsalazide disodiumPTG synthetase inhibitorPTGS2Inhibitor7YesNoNot available17981262|12950415
26LumiracoxibCOX-2PTGS2Inhibitor10YesYes314965322|15017615
27SalsalatePTG synthetase inhibitorPTGS2Inhibitor7YesYesNot available9711054|10452868
28Ginseng, AsianAdaptogenPTGS2Inhibitor4YesNoNot available17436372|15753395
29AntrafenineCyclooxygenase inhibitorPTGS2Inhibitor7YesNoNot available17604186|17612049
30AntipyrineCyclooxygenase inhibitorPTGS2Inhibitor2YesNoNot available11695253
31Tiaprofenic acidPTG synthetase inhibitorPTGS2Inhibitor8YesYesNot available2277128|11752352
32EtoricoxibCOX-2PTGS2Inhibitor12YesYes417573128|17164136
33ResveratrolPolyphenolic phytoalexinPTGS2Inhibitor2YesNo311752352
34NepafenacCOX-1 and COX-2PTGS2Inhibitor5YesNoNot available10850857
35ParecoxibCOX-2PTGS2Inhibitor3YesNoNot available10794682
36CaptoprilInhibitor of ACEMMP9Inhibitor5YesNoNot available12381651|19057128
37GlucosaminePrecursor of lipidsMMP9Antagonist6YesYes412405690|16616490
38RilonaceptInterleukin 1 inhibitorIL1RNBinder3YesNoNot available23319019|12815153
39DanazolGonadotropin inhibitorCCL2Inhibitor3YesNoNot available11056242|11216879
40EtanerceptFusion proteinTNFAntibody12YesNoNot available10375846|10405518
41AdalimumabMonoclonal antibodyTNFAntibody12YesNoNot available12044041|15022409
42InfliximabTNF-α blockerTNFInhibitor17YesNoNot available16456024|16052578
43InamrinonePDE3 inhibitorTNFInhibitor6YesNoNot available11805217|1611705
44GolimumabMonoclonal antibodyTNFAntibody6YesNoNot available21079302
45Certolizumab pegolTNF-α inhibitorTNFNeutralizer6YesNoNot available23620660
46PomalidomideCOX-2TNFInhibitor2YesNoNot available22917017
47SiltuximabMonoclonal antibodyIL6Antagonist4YesNoNot available8823310
48ClenbuterolBeta(2) agonistNGFStimulator6YesNoNot available10525174|10403429
49AfliberceptRecombinant fusion proteinVEGFAInhibitor7YesNoNot available22813448
50CarvedilolBeta adrenoceptor blockerVEGFAOther5YesNoNot available15071347|15942707
51Dalteparin sodiumAnticoagulantVEGFAInhibitor4YesNoNot available17692905|21091776
52RanibizumabVEGFA antagonistVEGFAInhibitor14YesNoNot available18046235|18054637
53BevacizumabVEG inhibitorVEGFAInhibitor10YesNoNot available15705858|11752352

Fifty-three drugs which met the criteria of targeting one of the candidate genes by an appropriate interaction were collected in the final list. Whether the drug has been approved by the FDA and whether it has been approved for OA are documented in the table. ACE: angiotensin-converting enzyme; TNF: tumor necrosis factor; PDE: phosphodiesterase; PTG: prostaglandin; VEGF: vascular endothelial growth factor.

4. Discussion

As the obesity rate and ages of the world population increase, the morbidity of OA is increasing especially in older adults. OA is one of the most common conditions causing decreased function, joint pain, and stiffness. To relieve symptoms by the improvement of joint function and reduction of pain, the classical treatment of OA is still using NSAIDs [24]. However, the long-term use of NSAIDs sometimes may cause headaches, gastrointestinal disturbances, and cardiovascular side effects. To address this question, we utilized bioinformatics tools to identify existing potential NSAIDs for OA treatment. Consequently, we identified eight targeting potential genes and 53 drugs associated with OA (Table 3).

The exact mechanism of OA pathogenesis remains to be elucidated, but the inflammation and cartilage degradation is the main pathological feature of OA [25]. The mechanism at a molecular level of inflammation and cartilage destruction involves multiple factors, such as matrix metalloproteinases (MMPs), COX-2, and ADAMTS (a disintegrin and metalloproteinases with thrombospondin motifs), in articular chondrocytes [26]. Mechanical stress, proinflammatory cytokines (including tumor necrosis factor-α (TNF-α), and interleukin-6 (IL-6), IL-17, IL-1β) can promote COX-2 and ADAMTS expression and accelerate the destruction of cartilage and OA development. PGE2 that is associated with the activity of TNF-α and IL-1β can lead to a large reduction in proteoglycan content in cartilage tissue [27]. The MMPs (especially 1 and 13) play major roles in cartilage matrix degradation and result in focal lesions in the articular surface [28]. Several studies reported increased expression of MMP-13 in OA. This difference may be caused as follows: PCR stochasticity, primer biases, multiple species in individual samples, and pooling of samples exerts too many uncertainties that could bias the results. Additionally, the sample size was insufficient. Proinflammatory cytokines can increase in nitric oxide (NO) and PGE2 synthesis and release [29]. Previous studies have shown that the anti-inflammatory effects of NSAIDs mainly depend on the function of inhibiting COX, reducing the generation of prostaglandins, which are the major mediator of pain and inflammatory response [30, 31].

In our results, the PTGS2 gene responsible for encoding COX-2 is corresponding to the most diverse class of drugs. Among them, 20 drugs are not applied for the treatment of OA but have been approved by the FDA to treat other diseases and obtained excellent outcomes. The mesalamine (5-ASA) compounds are commonly well tolerated and rare complications [32]. In the European [33] and the United States [34] guidelines, 5-ASA was recommended as a frontline agent for treating patients with mildly to moderately active ulcerative colitis. From the molecular structures’ point of view, the intensity of halogenation is positively correlated with drug therapy strength. It is considered that such a difference may also allow bromfenac (containing a bromine atom) wider ocular distribution [35]. Bromfenac is better than both dexamethasone and fluorometholone in preventing cystoid macular edema [36]. Although the molecular mechanism of sulfasalazine remained ambiguous, the drug is known to play a role that is immunomodulatory and anti-inflammatory [37]. It has also been the widely used treatment of other diseases, such as and ulcerative colitis, Crohn’s disease, and ankylosing spondylitis [34, 38, 39].

The platelet-derived growth factor superfamily secreted vascular endothelial growth factor A (VEGFA), a glycoprotein that mediates glomerular endothelial cell biological functions [40]. Glomerular visceral epithelial cells, known as podocytes, secrete VEGFA, a glycoprotein that mediates glomerular endothelial cell biological functions [40]. In the early stages of diabetes, podocytes increasing the expression of Vegfa can lead to both insulin deficiency and resistance in rodents [41]. Furthermore, Vegfa is also irreplaceable to induce the formation of cartilage. A significant positive correlation between the overexpression of Vegfa in osteochondrogenic cells and bone mass is also observed [42]. At the early stage of bone development, Vegfa was regarded as a regulator of both osteoblast differentiation and perichondrial vascularity [43].

There were also some limitations to the study. Firstly, the databases we used may be limited to get accurate information on genes concerning function or role. With the databases becoming more comprehensive, the results of the analysis will be more accurate. Although we have validated results in multiple biological databases, further molecular biological experiments of the large sample are indispensable to confirm the reliability of the results. We cannot ensure that a given gene is known to all relations of existing drugs and drug-gene interactions. Therefore, it is possible that we missed or ignored the drugs which could potentially be useful.

5. Conclusion

In conclusion, we have presented a way to investigate the potential candidate drugs which target the genes/pathways relevant to OA. With the databases and analytical tools evolved and improved, such a method may be used routinely to explore more diseases and to facilitate drug repositioning efforts in the pharmaceutical industry. From the present analysis, we identified a total of 53 potential drugs, of which 32 drugs targeted other OA-relevant pathways but have not yet been used in OA, which could give us a clue to novel therapeutic targets as potential treatments for OA. However, further direct evidence by molecular biological experiments is necessary to make clear their association with OA.

Abbreviations

OA:Osteoarthritis
NSAIDs:Nonsteroidal anti-inflammatory drugs
GO:Gene ontology
KEGG:Kyoto Encyclopedia of Genes and Genomes
PPI:Protein-protein interaction
COX:Cyclooxygenase
PGH2:Prostaglandin H2
PGE:Prostaglandin E
PTGS1:Prostaglandin-endoperoxide synthase
PTGS2:Prostaglandin-endoperoxide synthase
FDA:Food and Drug Administration
BP:Biological process
CC:Cellular component
MF:Molecular function
DGIdb:Drug-Gene Interaction Database
MMPs:Matrix metalloproteinases
TNF-α:Tumor necrosis factor-α
IL:Interleukin
NO:Nitric oxide
VEGFA:Vascular endothelial growth factor A.

Data Availability

The data used to support the findings of this study are included within the supplementary information file.

Conflicts of Interest

The authors declare no conflicts of interest with respect to the research, authorship, and/or publication of this article.

Acknowledgments

This study was supported by the Project of Innovation Platform for Fuzhou Health and Family Planning Commission (no. 2017-S-wp1).

Supplementary Materials

All gene names are shown in Appendix S1. (Supplementary Materials)

References

  1. C. Y. Wenham and P. G. Conaghan, “The role of synovitis in osteoarthritis,” Therapeutic Advances in Musculoskeletal Disease, vol. 2, no. 6, pp. 349–359, 2010. View at: Publisher Site | Google Scholar
  2. D. Lajeunesse, F. Massicotte, J.-P. Pelletier, and J. Martel-Pelletier, “Subchondral bone sclerosis in osteoarthritis: not just an innocent bystander,” Modern Rheumatology, vol. 13, no. 1, pp. 7–14, 2003. View at: Publisher Site | Google Scholar
  3. J. Bedson and P. R. Croft, “The discordance between clinical and radiographic knee osteoarthritis: a systematic search and summary of the literature,” BMC Musculoskeletal Disorders, vol. 9, no. 1, p. 116, 2008. View at: Publisher Site | Google Scholar
  4. B. R. da Costa, S. Reichenbach, N. Keller et al., “Effectiveness of non-steroidal anti-inflammatory drugs for the treatment of pain in knee and hip osteoarthritis: a network meta-analysis,” The Lancet, vol. 390, no. 10090, pp. e21–e33, 2017. View at: Publisher Site | Google Scholar
  5. S. Glyn-Jones, A. J. R. Palmer, R. Agricola et al., “Osteoarthritis,” The Lancet, vol. 386, no. 9991, pp. 376–387, 2015. View at: Publisher Site | Google Scholar
  6. M. Pitta, C. I. Esposito, Z. Li, Y.-Y. Lee, T. M. Wright, and D. E. Padgett, “Failure after modern total knee arthroplasty: a prospective study of 18,065 knees,” The Journal of Arthroplasty, vol. 33, no. 2, pp. 407–414, 2018. View at: Publisher Site | Google Scholar
  7. D. S. Jevsevar, “Treatment of osteoarthritis of the knee: evidence-based guideline, 2nd edition,” Journal of the American Academy of Orthopaedic Surgeons, vol. 21, no. 9, pp. 571–576, 2013. View at: Publisher Site | Google Scholar
  8. M. C. Hochberg, R. D. Altman, K. T. April et al., “American College of Rheumatology 2012 recommendations for the use of nonpharmacologic and pharmacologic therapies in osteoarthritis of the hand, hip, and knee,” Arthritis Care & Research, vol. 64, no. 4, pp. 465–474, 2012. View at: Publisher Site | Google Scholar
  9. J. R. Vane, “Inhibition of prostaglandin synthesis as a mechanism of action for aspirin-like drugs,” Nature New Biology, vol. 231, no. 25, pp. 232–235, 1971. View at: Publisher Site | Google Scholar
  10. M. C. Zweers, T. N. de Boer, J. van Roon, J. W. Bijlsma, F. P. Lafeber, and S. C. Mastbergen, “Celecoxib: considerations regarding its potential disease-modifying properties in osteoarthritis,” Arthritis Research & Therapy, vol. 13, no. 5, p. 239, 2011. View at: Publisher Site | Google Scholar
  11. C. Sostres, C. J. Gargallo, M. T. Arroyo, and A. Lanas, “Adverse effects of non-steroidal anti-inflammatory drugs (NSAIDs, aspirin and coxibs) on upper gastrointestinal tract,” Best Practice and Research Clinical Gastroenterology, vol. 24, no. 2, pp. 121–132, 2010. View at: Publisher Site | Google Scholar
  12. P. M. Kearney, C. Baigent, J. Godwin, H. Halls, J. R. Emberson, and C. Patrono, “Do selective cyclo-oxygenase-2 inhibitors and traditional non-steroidal anti-inflammatory drugs increase the risk of atherothrombosis? Meta-analysis of randomised trials,” BMJ, vol. 332, no. 7553, pp. 1302–1308, 2006. View at: Publisher Site | Google Scholar
  13. S. Trelle, S. Reichenbach, S. Wandel et al., “Cardiovascular safety of non-steroidal anti-inflammatory drugs: network meta-analysis,” BMJ, vol. 342, Article ID c7086, 2011. View at: Publisher Site | Google Scholar
  14. M. Gore, K.-S. Tai, A. Sadosky, D. Leslie, and B. R. Stacey, “Use and costs of prescription medications and alternative treatments in patients with osteoarthritis and chronic low back pain in community-based settings,” Pain Practice, vol. 12, no. 7, pp. 550–560, 2012. View at: Publisher Site | Google Scholar
  15. D. L. Scott, H. Berry, H. Capell et al., “The long‐term effects of non‐steroidal anti‐inflammatory drugs in osteoarthritis of the knee: a randomized placebo‐controlled trial,” Rheumatology, vol. 39, no. 10, pp. 1095–1101, 2000. View at: Publisher Site | Google Scholar
  16. C. P. Adams and V. V. Brantner, “Estimating the cost of new drug development: is it really $802 million?” Health Affairs, vol. 25, no. 2, pp. 420–428, 2006. View at: Publisher Site | Google Scholar
  17. S. Moosavinasab, J. Patterson, R. Strouse et al., “RE:fine drugs: an interactive dashboard to access drug repurposing opportunities,” Database, vol. 2016, Article ID baw083, 2016. View at: Publisher Site | Google Scholar
  18. E. Mosca, G. Bertoli, E. Piscitelli et al., “Identification of functionally related genes using data mining and data integration: a breast cancer case study,” BMC Bioinformatics, vol. 10, no. S12, p. S8, 2009. View at: Publisher Site | Google Scholar
  19. G. Jin and S. T. C. Wong, “Toward better drug repositioning: prioritizing and integrating existing methods into efficient pipelines,” Drug Discovery Today, vol. 19, no. 5, pp. 637–644, 2014. View at: Publisher Site | Google Scholar
  20. N. Novac, “Challenges and opportunities of drug repositioning,” Trends in Pharmacological Sciences, vol. 34, no. 5, pp. 267–272, 2013. View at: Publisher Site | Google Scholar
  21. J. Baran, M. Gerner, M. Haeussler, G. Nenadic, and C. M. Bergman, “pubmed2ensembl: a resource for mining the biological literature on genes,” PLoS One, vol. 6, no. 9, p. e24716, 2011. View at: Publisher Site | Google Scholar
  22. R. Saito, M. E. Smoot, K. Ono et al., “A travel guide to Cytoscape plugins,” Nature Methods, vol. 9, no. 11, p. 1069, 2012. View at: Publisher Site | Google Scholar
  23. A. H. Wagner, A. C. Coffman, B. J. Ainscough et al., “DGIdb 2.0: mining clinically relevant drug-gene interactions,” Nucleic Acids Research, vol. 44, no. D1, pp. D1036–D1044, 2016. View at: Publisher Site | Google Scholar
  24. R. R. Bannuru, M. C. Osani, E. E. Vaysbrot et al., “OARSI guidelines for the non-surgical management of knee, hip, and polyarticular osteoarthritis,” Osteoarthritis and Cartilage, vol. 27, no. 11, pp. 1578–1589, 2019. View at: Publisher Site | Google Scholar
  25. W. B. Van den Berg, “Osteoarthritis year 2010 in review: pathomechanisms,” Osteoarthritis and Cartilage, vol. 19, no. 4, pp. 338–341, 2011. View at: Publisher Site | Google Scholar
  26. M. B. Goldring and M. Otero, “Inflammation in osteoarthritis,” Current Opinion in Rheumatology, vol. 23, no. 5, p. 471, 2011. View at: Publisher Site | Google Scholar
  27. S. C. Mastbergen, J. W. J. Bijlsma, and F. P. J. G. Lafeber, “Synthesis and release of human cartilage matrix proteoglycans are differently regulated by nitric oxide and prostaglandin-E2,” Annals of the Rheumatic Diseases, vol. 67, no. 1, pp. 52–58, 2008. View at: Publisher Site | Google Scholar
  28. J. Martel-Pelletier, C. Boileau, J.-P. Pelletier, and P. J. Roughley, “Cartilage in normal and osteoarthritis conditions,” Best Practice & Research Clinical Rheumatology, vol. 22, no. 2, pp. 351–384, 2008. View at: Publisher Site | Google Scholar
  29. M. B. Goldring, “Osteoarthritis and cartilage: the role of cytokines,” Current Rheumatology Reports, vol. 2, no. 6, pp. 459–465, 2000. View at: Publisher Site | Google Scholar
  30. E. V. Tchetina, J. A. Di Battista, D. J. Zukor, J. Antoniou, and A. R. Poole, “Prostaglandin PGE2 at very low concentrations suppresses collagen cleavage in cultured human osteoarthritic articular cartilage: this involves a decrease in expression of proinflammatory genes, collagenases and COL10A1, a gene linked to chondrocyte hypertrophy,” Arthritis Research & Therapy, vol. 9, no. 4, p. R75, 2007. View at: Publisher Site | Google Scholar
  31. M. Attur, H. E. Al-Mussawir, J. Patel et al., “Prostaglandin E2 exerts catabolic effects in osteoarthritis cartilage: evidence for signaling via the EP4 receptor,” The Journal of Immunology, vol. 181, no. 7, pp. 5082–5088, 2008. View at: Publisher Site | Google Scholar
  32. J. H. Chidlow Jr., D. Shukla, M. B. Grisham, and C. G. Kevil, “Pathogenic angiogenesis in IBD and experimental colitis: new ideas and therapeutic avenues,” American Journal of Physiology-Gastrointestinal and Liver Physiology, vol. 293, no. 1, pp. G5–G18, 2007. View at: Publisher Site | Google Scholar
  33. A. Dignass, J. O. Lindsay, A. Sturm et al., “Second European evidence-based consensus on the diagnosis and management of ulcerative colitis part 2: current management,” Journal of Crohn’s and Colitis, vol. 6, no. 10, pp. 991–1030, 2012. View at: Publisher Site | Google Scholar
  34. A. Kornbluth and D. B. Sachar, “Ulcerative colitis practice guidelines in adults: American college of gastroenterology, practice parameters committee,” American Journal of Gastroenterology, vol. 105, no. 3, pp. 501–523, 2010. View at: Publisher Site | Google Scholar
  35. T. Kida, S. Kozai, H. Takahashi, M. Isaka, H. Tokushige, and T. Sakamoto, “Pharmacokinetics and efficacy of topically applied nonsteroidal anti-inflammatory drugs in retinochoroidal tissues in rabbits,” PLoS One, vol. 9, no. 5, Article ID e96481, 2014. View at: Publisher Site | Google Scholar
  36. Q. W. Wang, K. Yao, W. Xu et al., “Bromfenac sodium 0.1%, fluorometholone 0.1% and dexamethasone 0.1% for control of ocular inflammation and prevention of cystoid macular edema after phacoemulsification,” Ophthalmologica. Journal International D’ophtalmologie. International Journal of Ophthalmology. Zeitschrift Fur Augenheilkunde, vol. 229, no. 4, pp. 187–194, 2013. View at: Publisher Site | Google Scholar
  37. G. Plosker and K. Croom, “Sulfasalazine: a review of its use in the management of rheumatoid arthritis,” Drugs, vol. 65, no. 18, p. 2591, 2005. View at: Publisher Site | Google Scholar
  38. N. Akkoc, S. van der Linden, and M. A. Khan, “Ankylosing spondylitis and symptom-modifying vs disease-modifying therapy,” Best Practice & Research Clinical Rheumatology, vol. 20, no. 3, pp. 539–557, 2006. View at: Publisher Site | Google Scholar
  39. G. R. Lichtenstein, E. V. Loftus, K. L. Isaacs, M. D. Regueiro, L. B. Gerson, and B. E. Sands, “ACG clinical guideline: management of Crohn’s disease in adults,” American Journal of Gastroenterology, vol. 113, no. 4, pp. 481–517, 2018. View at: Publisher Site | Google Scholar
  40. V. Eremina, S. Cui, H. Gerber et al., “Vascular endothelial growth factor a signaling in the podocyte-endothelial compartment is required for mesangial cell migration and survival,” Journal of the American Society of Nephrology, vol. 17, no. 3, pp. 724–735, 2006. View at: Publisher Site | Google Scholar
  41. A. Flyvbjerg, F. Dagnaes-Hansen, A. S. De Vriese, B. F. Schrijvers, R. G. Tilton, and R. Rasch, “Amelioration of long-term renal changes in obese type 2 diabetic mice by a neutralizing vascular endothelial growth factor antibody,” Diabetes, vol. 51, no. 10, pp. 3090–3094, 2002. View at: Publisher Site | Google Scholar
  42. C. Maes, S. Goossens, S. Bartunkova et al., “Increased skeletal VEGF enhances β-catenin activity and results in excessively ossified bones,” The EMBO Journal, vol. 29, no. 2, pp. 424–441, 2010. View at: Publisher Site | Google Scholar
  43. X. Duan, Y. Murata, Y. Liu, C. Nicolae, B. R. Olsen, and A. D. Berendsen, “Vegfa regulates perichondrial vascularity and osteoblast differentiation in bone development,” Development, vol. 142, no. 11, pp. 1984–1991, 2015. View at: Publisher Site | Google Scholar

Copyright © 2021 Rong-Guo Yu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views328
Downloads658
Citations

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.