*Corresponding Author:
Pratima Tatke
Department of Pharmaceutical Chemistry, CU Shah College of Pharmacy, Shreemati Nathibai Damodar Thackersey (SNDT) Women’s University, Mumbai, Maharashtra 400049, India
E-mail: drpratimatatke@gmail.com
Date of Received 12 March 2020
Date of Revision 01 December 2021
Date of Acceptance 18 November 2022
Indian J Pharm Sci 2022;84(6):1398-1406  

This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 3.0 License, which allows others to remix, tweak, and build upon the work non-commercially, as long as the author is credited and the new creations are licensed under the identical terms

Abstract

Glycyrrhiza glabra (licorice), is used for therapeutic purposes since ages. Neuropathic pain is a disease with complex pathophysiology, involving dysfunction in multiple ion channel receptors such as sodium, transient receptor potential channels, neurotransmitters, inflammatory mediators and reactive oxygen species. Network pharmacology is newly emerging approach dealing with understanding mechanism of actions of druggable molecules with the multiple targets. This approach is suitable to deal with AYUSH based medicines such as extracts containing many bioactives, acting on multiple targets. The aim of the present study is to construct a network of licorice to predict its actions against neuropathic pain. Various knowledge databases were used for collecting information regarding bioactives, their targets and genes associated with these targets or pathophysiology of neuropathic pain. These databases include Dr. Duke's Phytochemical and Ethnobotanical Database, PubMed, Google scholar, PubChem, BindingDB, Kyoto encyclopedia of genes and genomes, DisGeNET and Therapeutic Target Database. The network was constructed using Cytoscape3.7.2. 19 bioactive were shortlisted by extensive search having targets on humans and rats using BindingDB. These targets were searched on Kyoto encyclopedia of genes and genomes, DisGeNET for their association with neuropathic pain. Targets were shown to be part of pathways of neuropathic pain, helping to predict mechanism of action of bioactives. Hence, this approach can prove helpful in modern pharmacological research.

Keywords

Network pharmacology, licorice, neuropathic pain, transient receptor potential, herbal extracts

Network pharmacology is newly emerging approach which deals with understanding mechanism of actions of drug like molecules with its multiple Targets (TARs). This method uses various databases to retrieve information of Bioactives (BAs), their TARs and involvement of the TARs in the pathophysiology of disease. It also includes using software Cytoscape, to construct a network highlighting interrelation between the medicinal plant, its BAs and their TARs. The modern drug discovery process follows one drug-one target approach. This approach well suited for synthetic medicines having a single active ingredient. But there are limitations of this approach when it comes to deal with AYUSH medicines. Ayurvedic drugs as well as extracts of medicinal plants consist of many BAs, which would be acting on multiple TARs at the same time, producing synergistic actions[1]. Many disorders like diabetes, cardiovascular diseases, obesity and even Neuropathic Pain (NP) have complex pathophysiology, which could be treated well with multi targeted approach. Network pharmacology can be an important tool for understanding complex mechanism of action of various BAs present in herbal extracts/polyherbal medicines[2].

Glycyrrhiza glabra (G. glabra) also known as licorice and sweet wood, is native to the Mediterranean and certain areas of Asia. Historically, the dried rhizome and root of this plant were employed medicinally by the Egyptian, Chinese, Greek and Indian civilizations as an expectorant and carminative[3]. Licorice is known to have antibacterial, antiviral, antioxidant, antiinflammatory, anti-tussive, anti-ulcer, anticancer and hepatoprotective activities[4].

NP has been described as “the worst of tortures which a nerve wound may experience”[5]. NP is prominently experienced as a complication of diabetes, herpes-zoster infection, Human Immunodeficiency Virus (HIV) infection and cancer[6]. The pathophysiology of NP is complex. Peripheral NP is divided into two types depending on its pathophysiology, spontaneous pain (stimulus-independent pain) and evoked pain (stimulus dependent pain)[7]. Stimulus dependent pain results due to over activity of nociceptors. These receptors get stimulated by various stimuli such as mechanical stimuli, thermal stimuli and chemical stimuli. Mechanosensitive nociceptors are activated by pressure and mechanical stress. Thermosensitive receptors are sensitive to heat above 43°-45°[8]. Inflammation within the nerve tends to release variety of endogenous and exogenous substances like prostaglandins, bradykinin, histamine, serotonin, Tumor Necrosis Factor-alpha (TNF-α) and nerve growth factor, which also stimulates nociceptors. Nerve injury also induces up regulation of receptors of TRP family, such as TRPV1[8]. Formation of Advanced Glycosylated End products (AGEs) is related to diabetic neuropathy. It is also been found that in diabetic neuropathy, AGEs and Receptors for AGE (RAGE) interactions induce oxidative stress, result in up regulation of Nuclear Factor-Kappa B (NFkB) and various NF-kB-mediated proinflammatory genes, inducing altered pain sensation[9]. Apart from sodium and alpha adrenoceptors, calcium channels are also shown to be involved in pathophysiology of NP. The aim of this paper is to evaluate phytoconstituents of licorice by network pharmacological techniques to predict their mechanism of action against NP.

Materials and Methods

The information regarding BAs of licorice, its TARs involved in NP, genes associated with NP were collected from various databases as follows.

BAs:

The list of BAs from licorice was obtained from Dr. Duke's Phytochemical and Ethno botanical databases[10]. (Version is 1.9.12.6-Beta_2019-02-06). This website provides information about different chemicals present in plant. The website also gives list of chemicals in different parts of plants and it also provides a comprehensive list of chemicals with their [11]. For the current study, list of BAs from roots and rhizomes of licorice was downloaded.

The search of BAs was further narrowed down using PubMed and Google Scholar. Articles dealing with Glycyrrhiza uralensis, Glycyrrhiza inflata and Glycyrrhiza echinata were excluded from search as scope of this study were limited to G. glabra. Articles describing chemical analysis of aerial parts of licorice were excluded from search as the aim of the study is restricted to phytoconstituent of licorice root. Further molecular structures of these BAs were searched from databases like PubChem and Chem Spider. The names of BAs, their molecular formula, PubChem IDs were noted down and its Three Dimensional (3D) structures were downloaded in Standard Delay Format (SDF).

In order to search TARs for the BAs, it is necessary to look into the pharmacokinetics of the BAs. Hence these BAs were first searched in a database known as Swiss-Absorption, Distribution, Metabolism and Excretion (ADME). Swiss-ADME is a free web tool developed by Swiss Institute of Bioinformatics. This database helps in predicting about physicochemical properties, lipophilicity and water solubility. It can also predict pharmacokinetic parameters like GI absorption, blood brain barrier permeation, effect of hepatic metabolic enzymes and skin permeation rate. The database can also predict “drug likeness” and “lead likeness” of the compound[12]. The BAs which are predicted to be bioavailable as per Swiss-ADME were considered for TAR search.

TARs:

The TARs for the BAs were searched using Binding Database (DB). BindingDB is a freely available database which focuses mainly on interactions of proteins which can be considered as TARs with small drug like molecules[13]. The compounds can be searched in this database using SD Files or SMILES format. The TARs for the BAs were searched with the help of the tab, “Find my compound’s target”, by submitting the structures in SD files format. BindingDB searches the TARs for drug like molecules based upon structural similarities between the ligand and the target[14]. To get the optimum list of TARs the similarity search was adopted for 85 % of similarity i.e., the ligand and TARs are ≥85 % structurally similar to each other. For each of the obtained target, the information was collected as; target name, number of hits, percentage similarity, UniProt ID, UniProt name, name of the organism and gene associated with target. BindingDB is further connected to UniProt/SwissProt, Kyoto Encyclopedia of Genes and Genomes (KEGG), DrugBank and Google Scholar. UniProtKB was used for collecting UniProt ID, UniProt name, associated gene and name of organism. The BAs were also searched in database ChEMBL. ChEMBL is another database which gives information about TARs of drug like compounds. The BA was searched in this database using name and list of TARs was collected. The common TARs were excluded using filter in Excel file. TARs present in human beings and rats were included in final TARs list.

KEGG:

KEGG pathway database is a database giving detailed information about various pathways involved in pathogenesis of various diseases[15,16]. KEGG has not enlisted any separate pathway for NP. Hence the TARs were manually searched in pathways related to pathophysiology of NP i.e., regulation of TRP channels, arachidonic acid pathway and inflammation etc.

DisGeNET:

DisGeNET is a discovery platform containing one of the largest publicly available collections of genes and variants associated to human diseases. This database integrates data from repositories, GWAS catalogues, animal models and the scientific literature[16]. The version 6.0 was used for this study. This version contains 628 685 Gene-Disease Associations (GDAs), between 17 549 genes and 24 166 diseases, disorders, traits and clinical or abnormal human phenotypes[17]. Gene IDs which were collected from UniProt were submitted to this database, to check if any of these TARs are associated genetically with NP. TARs associated with peripheral neuropathy, sensory neuropathy, diabetic neuropathy, pain, neuralgia, sciatic neuropathy, inflammation were collected and a combined list was prepared.

Therapeutic Target Database (TTD):

TTD is a database which provides information about various therapeutic proteins and nucleic acid TARs of the targeted disease[18]. It also provides pathway information and the corresponding drugs directed at each of these TARs. This database also has links to relevant databases containing information about target function, sequence, 3D structures, ligand binding properties, enzyme nomenclature and drug structure, therapeutic class, clinical development status. All information provided in this database is fully referenced[19]. The list of therapeutic TARs for NP was collected from this database. Out of different collected TARs, two TARs were found to be successful therapeutic TARs for NP.

Network construction:

Network pharmacology uses various databases to retrieve information of BAs, their TARs and involvement of the TARs in the pathophysiology of disease. It also includes using software Cytoscape, to construct a network highlighting interrelation between the medicinal plant, its BAs and their TARs. Cytoscape, a Java-based open-source software platform, is a useful tool for visualizing molecular interaction networks and integrating them with any type of attribute data[20]. In addition to the basic set of features for data integration, analysis and visualization, additional features are available in the form of apps, including network and molecular profiling analysis and links with other databases. Other software, such as Gephi, an exploration platform for networks and complex systems and cell illustrator, a Java-based tool specialized in biological processes and systems can also be used for building networks[21].

A network consists of nodes and edges. The entities in which interaction is studied are called as nodes and lines connecting nodes are called as edges. The nodes of the networks in the current study include licorice, the BAs of licorice, the TARs of BAs (genes) and various neuropathy related pathways associated with those TARs. The networks were constructed using Cytoscape 3.7.2, Cytoscape is useful open-source software platform for visualizing complex networks and integrating them with any type of attribute data[22]. Excel format was used to input data in Cytoscape and to develop a network. Cytoscape network analyzer tool was used to analyze the networks.

Results and Discussion

The list of BAs from licorice was obtained from Dr. Duke's phytochemical and ethno botanical databases. The plant name G. glabra was chosen. The database gave list of 538 compounds present in licorice, out which 217 have different activities. The search of BAs from rhizomes and roots retrieved 13 and 130 BAs respectively. Further extensive search was carried out using PubMed and Google scholar and articles detailing about component analysis of licorice extracts. The articles enlisting constituents of ethanol, methanol and ethyl acetate extracts of licorice using analytical techniques like high-performance liquid chromatography, liquid chromatography–mass spectrometry, highperformance thin-layer chromatography, electrospray ionization mass spectrometry were searched and list of BAs was prepared. A comprehensive list of 36 BAs was made. These BAs were first searched in Swiss-ADME database in order to predict their bioavailability and drug likeliness as given in Table 1.

Sl. no Name of phytoconstituent Molecular formula Molecular weight (g/mol) GI absorption Drug-likeness (whether following Lipinski Rule)
1 18βGA C30H46O4 470.68 High Yes; 1 violation: MLOGP>4.15
2 Isoliquiritegenin C15H12O4 256.25 High Yes; 0 violation
3 Isoliquiritin C21H22O9 418.39 Low Yes; 1 violation: NH or OH>5
4 Kumatakenin C17H14O6 314.29 High Yes; 0 violation
5 Licoflavone B C25H26O4 390.47 High Yes; 0 violation
6 Licoflavonol C20H18O6 354.35 High Yes; 0 violation
7 Licoflavone A C20H18O4 322.35 High Yes; 0 violation
8 Neoisoliquiritine C21H22O9 418.4 Low Yes; 1 violation: NH or OH>5
9 Neoliquiritin C21H22O9 418.4 Low Yes; 0 violation
10 Liquiritin C21H22O9 418.4 Low Yes; 0 violation
11 Glycyrrhizin C42H62O16 822.93 Low No; 3 violations: MW>500, N or O>10, NH or OH>5
12 Glabrene C20H18O4 322.4 High Yes; 0 violation
13 Licochalcone C21H22O4 338.4 High Yes; 0 violation
14 Glabridin C20H20O4 324.37 High Yes; 0 violation
15 Glabrol C25H28O4 392.49 High Yes; 0 violation
16 Glycycoumarin C21H20O6 368.4 High Yes; 0 violation
17 Glisoflavone C21H20O6 368.4 High Yes; 0 violation
18 Licoisoflavanone C20H18O6 354.4 High Yes; 0 violation
19 Isowighteone C20H18O5 338.4 High Yes; 0 violation
20 Isoangustone A C25H26O6 422.5 Low Yes; 0 violation
21 Glyasperin D C22H26O5 370.4 High Yes; 0 violation
22 Glicophenone C20H22O6 358.4 High Yes; 0 violation
23 Liquiritin C21H22O9 418.4 Low Yes; 0 violation
24 Liquiritigenin C15H12O4 256.25 High Yes; 0 violation
25 Licoisoflavone A C20H18O6 354.4 High Yes; 0 violation
26 Glabridin C20H20O4 324.4 High Yes; 0 violation
27 Licoricone C22H22O6 382.4 High Yes; 0 violation
28 Neoliquiritin 2''-apioside C26H30O13 550.5 Low No; 3 violations: MW>500, N or O>10, NH or OH>5
29 Licorice glycoside A C36H38O16 726.7 Low No; 3 violations: MW>500, N or O>10, NH or OH>5
30 Isopropyl apiosylglucoside C14H26O10 354.35 Low Yes; 1 violation: NH or OH>5
31 Glycyrrhizaisoflavone C C21H20O6 368.4 High Yes; 0 violation
32 Glycyrrhizaisoflavone C C21H20O6 368.4 High Yes; 0 violation
33 Liquiritic acid C30H46O4 470.7 High Yes; 1 violation: MLOGP>4.15
34 Licoagroside B C18H24O12 432.4 Low Yes; 1 violation: N or O>10
35 Licorice glycoside B C35H36O15 696.6 Low No; 3 violations: MW>500, N or O>10, NH or OH>5
36 Licorice saponin A3 C48H72O21 985.1 Low No; 3 violations: MW>500, N or O>10, NH or OH>5
37 Licorice saponin G2 C42H62O17 838.9 Low No; 3 violations: MW>500, N or O>10, NH or OH>5

Table 1: Predicted Bioavailability of Licorice Phytoconstituents as per Swiss-Adme Database

The BAs listed in Table 1 were searched in BindingDB for the TARs, using the parameters mentioned above. Out of 36 BAs, 19 were found to be having TARs in pathways related to pathophysiology of NP and diabetic neuropathy. The network of licorice and its BAs were prepared using Cytoscape, as shown in fig. 1. The list of BAs with their TARs is given in Table 2.

pharmaceutical-sciences-network

Fig. 1: Network of licorice and BAs with their PubChem IDs

Sl. no Name of bioactive TARs involved in NP
1 18βGA HSD11B1, HSD11B2, Ar, CYP19A1, CACNA1C, PTGES, PRKCH
2 Isoliquiritigenin AKR1B1, Alox5, ABCG2, BACE1, KCNA3
3 Kumatakenin Adora1,ADORA3,Adora2a,Akr1b1,ABCG2,CA2, HSD17B2, NOX4
4 Licoflavone B CYP19A1, PPARG, Htr2a
5 Licoflavonol Adora1, ADORA3, Adora2a, ABCB1
6 Licoflavone A Htr2a, Adora1, Adora2a
7 Neoisoliquiritine AKR1B1
8 Glycyrrhizin HSD11B1, HSD11B2
9 Isowighteone PTPN1
10 Isoangustone A PTPN2
11 Liquiritigenin Ar, CA2, HSP90AB1, RIPK1, PSMB5
12 Licoisoflavone A PTPN1
13 Licoricone PTPN1
14 Licorice glycoside A AKR1B1
15 Isopropyl apiosylglucoside FGF1, ABL1, CCNA2, FGF2, HDAC10 HDC3, HDC5, HDC8, HDC9
16 Liquiritic acid HSD11B1, HSD11B2, HSD11B3, HSD11B4, CACNA1C, PTGES, PRKCH
17 Licorice glycoside B AKR1B1
18 Licorice saponin A3 HSD11B1, HSD11B2
19 Licorice saponin G2 HSD11B1, HSD11B2

Table 2: List of Bas and their Tars Involved in Np (Gene Names)

Each of the TARs was further searched in the KEGG database for their involvement in neuropathy related pathophysiology such as inflammation (TNF-α, bradykinin, serotonin, prostaglandin, histamine) arachidonic acid pathway, Na channel blocking, TRP channel regulation, mitogen-activated protein kinase pathway etc. The TARs were searched on DisGeNET for their association with NP. The list of TARs for NP was also retrieved from TTD. As per KEGG database, TARs Protein Kinase C Eta (PRKCH) and 5-Hydroxytryptamine Receptor 2A (HTR2A) were found to be associated with regulation of TRP channels. TARs Prostaglandin E Synthase (PTGES) and Arachidonate 5-Lipoxygenase (ALOX5) were found to be associated with arachidonic pathway. TARs Protein Tyrosine Phosphatase Non-Receptor Type (PTPN)-1, PTPN2 and Nicotinamide Adenine Dinucleotide Phosphate Oxidase 4 (NOX4) were found to be associated with insulin resistance and AGE-RAGE pathway involved in diabetic complications such as diabetic neuropathy. Two TARs with UniProt name, adenosine receptor A1 and adenosine receptor A2a are included in list of TARs on TTD as successful TARs. Genes Androgen Receptor (AR), PTGES, Aldo-Keto Reductase Family 1 Member B (AKR1B1), Adenosine Triphosphate (ATP) binding cassette subfamily G member 2, potassium voltage-gated channel subfamily A member 3, adenosine A1 receptor, carbonic anhydrase 2, NOX4, peroxisome proliferator activated receptor gamma, HTR2A, ATP binding cassette subfamily B member 1 were found to be associated with NP as per DisGeNet database. TARs Hydroxysteroid 11-Beta Dehydrogenase 1 (HSD11B1), HSD11B1, HSD11B2, acetylcholinesterase, heat shock protein 90 alpha family class B member 1, receptor interacting serine/threonine kinase 1, proteasome 20S subunit beta 5, Fibroblast Growth Factor (FGF) 1, ABL proto-oncogene 1, non-receptor tyrosine kinase, cyclin A2, FGF2, histone deacetylase 10, Histidine Decarboxylase (HDC) 3, HDC5, HDC8, HDC9, PTGES, Calcium Voltage-Gated Channel Subunit Alpha1 C (CACNA1C) were found be novel TARs for treatment of NP as per literature search[23-26]. The network showing the plant licorice, its BAs and TARs associated with are shown in fig. 2.

pharmaceutical-sciences-licorice

Fig. 2: Network showing licorice, its BAs and TARs associated with it

Finally, a network of licorice BAs and their TARs involved in pathophysiology of NP was constructed fig. 3. These networks were constructed using Excel data sheets. TARs having involvement with pathophysiology of NP only were included in the network. Other TARs were excluded.

pharmaceutical-sciences-their

Fig. 3: Network of licorice with its BAs, TARs and their association with NP

Further the network was analyzed with the help of tools available with network analyzer and taking help of various layout styles. The network was viewed in degree sorted circle layout fig. 4. Node degree of a node is defined as number of edges linked to a node. Degree sorted circle layout arranges a network depending on node degree distribution. The node with highest degree appears at 6’o clock position, with biggest circle. The nodes with decreasing degrees appear in anticlockwise direction, with decreasing size of circle. The network analysis is shown in Table 3.

pharmaceutical-sciences-involve

Fig. 4: Degree sorted network of licorice, showing its BAs and TARs involved in NP

Licorice Network
Number of nodes 69
Network density 0.042
Network heterogeneity 0.862
Self-loops 0
Network diameter 8
Network radius 4
Isolated nodes 0
Average number of neighbors 2.841

Table 3: Analysis of Licorice Network Using Cytoscape Network Analyzer

The current study generated network provides a molecular overview of antinociceptive effect of licorice by bringing the BAs, their TARs, interacting proteins and pain pathways together. Licorice is known to have various active constituents with proven preclinical anti-inflammatory and antioxidant activities. Chen et al., has reported mechanism of action licorice by network pharmacological techniques. As per the network reported by them, licorice constituents can work on TARs involved in anti-inflammatory process, neuroactive and vasoactive amine metabolisms, hormone transcription regulation and glyceraldehyde oxidoreductase activity[27].

From the degree sorted circle layout, it can be seen that 18β Glycyrrhetinic Acid (GA), glycyrrhizin, liquiritin and liquiritigenin are few of the important constituents of the licorice extract. 18β GA has shown sodium channel blocking activity and antiinflammatory activity[28,29]. Liquiritin and liquiritigenin have shown to reduce lipopolysaccharide induced release pro-inflammatory mediators such as TNF-α, Interleukin (IL)-1β and IL-6, in various pre-clinical models[30]. The network pharmacology technique, helped to predict various TARs through which BAs are working to attenuate NP. Phytoconstituent 18βGA was found to work through HSD11B1, HSD11B2, CACNA1C, PTGES and PRKCH receptors. HSD11B1, HSD11B2 are associated with enzyme 11β-Hydroxysteroid Dehydrogenase (11-βOHSD). It is reported that 18βGA is a potent inhibitor of 11-βOHSD. It is postulated that its antiinflammatory activity is the expression of reduced conversion of cortisol to its biologically inactive product cortisone[31]. A research report by Kiso et al., also claims involvement of HSD11B1, HSD11B2 in pathophysiology of NP, further straightening the claim[23]. CACNA1C is associated with voltage dependent L-Type Calcium Channels (LTCs). LTCs are mostly post-synaptic channels regulating neuronal firing and gene expression. Their involvement in pain processes was also recognized recently. When applied locally in the dorsal horn of the spinal cord, the three families of LTC pharmacological blockers–dihydropyridines, phenylalkylamines and benzothiazepines proved effective in altering short-term sensitization to pain, inflammationinduced hyper excitability and neuropathy-induced allodynia[32]. Target PRKCH is associated with TRP channels. TRP channels are activated in damaging pressure, extreme temperature and toxic chemicals. TRP channels including TRPV1, TRPA1 and TRPM8 play a central role in the sensitization of nociceptive transduction as mentioned earlier. TRPM8 is also listed as a successful target in TTD. Target PTGES is associated with Prostaglandin E synthase. Prostaglandin E2 (PGE2) is considered to be the principal proinflammatory prostanoid and to play an important role in nociceptive processing and sensitization in the spinal cord as well as in the periphery by acting through four PGE receptor. Blockade of PGE2 have shown to attenuate symptoms of NP[24].

Apart from these TARs, other phytoconstituent were shown to interact with two more TARs, namely PTPN1 and AKRB1, associated with insulin resistance and AR pathway respectively. During hyperglycemia, these TARs are associated with insulin resistance and aldose reductase pathway. This pathway is responsible for the reduction of glucose to sorbitol via the enzyme AR and subsequent accumulation of sorbitol in cells. Excessive accumulation of sorbitol in peripheral leads to decrease in nerve conduction, which is a prime symptom of neuropathy. AR inhibitors have shown to produce effect in case diabetic neuropathy. Hence association BAs with target AKR1B1 could be beneficial for treating diabetic neuropathy [33].

Network pharmacology approach helps to predict mechanism of action of drug like molecules. This technique is helpful for studying medicinal plants or AYUSH medicines having multiple ingredients and acting on various TARs. This study reports potential of licorice in treatment of NP.

Network pharmacology is a modern approach to understand the multi-target working of Indian traditional medicinal system. One of the major hindrances in acceptance of herbal therapies in global market is proof of efficacy with respect of mechanism of action. This proof can be obtained by in vitro and in vivo studies of medicinal plants, but these methods are time consuming and costly. This in silico mechanistic study offers a rapid and reliable way of predicting BA-TAR interaction and a lead molecule for further potential pharmacological activity or drug discovery.

Evaluation of plants using this approach will help in understanding various therapeutic uses, with providing insight into its mechanism of action. Network pharmacology approach can also be looked as a preliminary screening technique to predict the TARs/pathways. The extract would act on and can help in reducing number of animals in pre-clinical studies. Further, inclusion of databases like Search Tool for the Retrieval of Interacting Genes (STRING), into this approach can also help to predict the protein-protein interactions of TARs, giving deeper insights into molecular mechanisms. This approach is cost-effective, time saving and as well as alternative to animal experimentation, hence can be a boon to modern pharmacological research.

This article attempted to evaluate effectiveness of licorice constituents for NP by applying network pharmacology. Licorice contains various phytoconstituent working on multiple TARs involved in NP. In our laboratory, we have tested ethanol, methanol and ethyl acetate extracts of licorice in preclinical models of NP and they have shown beneficial effects[34]. Hence this plant should be explored further for treatment of NP.

Acknowledgements

The authors are grateful to Dr. Uma Chandran and Dr. Bhushan Patwardhan, Interdisciplinary School of Sciences, SPPU, Pune, Maharashtra, for their support and guidance.

Conflict of Interest

The authors declared no conflict of interests.

References