Interaction network among functional drug groups
© Lee et al.; licensee BioMed Central Ltd. 2013
Published: 16 October 2013
More attention has been being paid to combinatorial effects of drugs to treat complex diseases or to avoid adverse combinations of drug cocktail. Although drug interaction information has been increasingly accumulated, a novel approach like network-based method is needed to analyse that information systematically and intuitively
Beyond focussing on drug-drug interactions, we examined interactions between functional drug groups. In this work, functional drug groups were defined based on the Anatomical Therapeutic Chemical (ATC) Classification System. We defined criteria whether two functional drug groups are related. Then we constructed the interaction network of drug groups. The resulting network provides intuitive interpretations. We further constructed another network based on interaction sharing ratio of the first network. Subsequent analysis of the networks showed that some features of drugs can be well described by this kind of interaction even for the case of structurally dissimilar drugs.
Our networks in this work provide intuitive insights into interactions among drug groups rather than those among single drugs. In addition, information on these interactions can be used as a useful source to describe mechanisms and features of drugs.
Recently, various drug-related properties including therapeutic target, off-target, activity, toxicity, pharmacokinetic properties and side effect have been successfully described by some elegant methods such as analysis of chemical structures , docking simulation [2–4], chemical-genetic profile [5–7], connectivity map , target-sharing information, and side-effect similarity . To go further and to treat complex disease like cancer, these recent advances in drug discovery make people pay more attention to a novel approach named polyphamacology  than to the design of a particular drug which could target disease-causing genes, because most of complex diseases result from vast range of biological abnormalities. This attention to polypharmacology also involves increasing attention to detect rational combination of older drugs to treat novel diseases. In this paradigm, significant attention to drug-drug interaction (DDI) is necessarily critical to detect ideal drug combinations to have synergistic clinical effects .
Another field, where DDI information would be significant, is combinatorial drug toxicology . In some cases, DDIs exhibit adverse drug reactions (ADRs)  or critical threat to patients with multiple medications [14–16]. For example, about 4% of causes of death of cancer patients were originated from DDIs . The risk of toxicity caused by unwanted polypharmacology is emerging threat to public health . Thus, it is quite important to detect and predict DDIs in both drug discovery and drug toxicology.
It is essential to use system-based approach, or network-based in order to detect beneficial and harmful drug combination, because DDI necessarily involves the complexity. Many recent approaches regarding DDI have utilized information from networks among human diseases, biology, and chemistry [19, 20] such methods as prediction of DDIs through protein-protein interaction (PPI) network , prediction of ADRs by integrating PPI network and drug structures  and a method using drug-target interaction network . There have been also some researches on DDI using network analysis. Hu et al. constructed DDI network of 966 drugs and simply reported some properties of their DDI networks . More advanced work by Xu et al. (drug cocktail network) developed drug combination predictor based on DDI network . However, analysis of systematic DDI relationship in terms of functional context has not been made so far.
In this study, we perform systematic analysis of DDIs based on functional group of drugs and construct the network which consists of interactions among the drug groups (DGs). The classification we used are based on the anatomical therapeutic chemical (ATC) classification system . The resulting DG-DG interaction network provides more contextual and intuitive view on DDI. In addition, the secondary network, whose edge represents that two DGs share a number of DG-DG interactions, is also constructed to detect a set of DGs showing similar DG-DG interaction patterns. We show that the DGs that have similar interactions share many similar drug features, suggesting that DDIs contain the information about drug mechanisms. Moreover, we question whether DDI information can be used to infer drug mechanisms by investigating common drug features for interacting drug pairs. The results show that some drug features such as metabolizing enzyme, drug function, and target pathway are well-described by DDI even for the case that chemical structure similarity is low.
Results and discussion
Statistically significant DG-DG interactions
Like this example, our analysis on drug-group interactions can provide the information about not only the meaningful drug-group interactions but also missing drug-drug interactions.
Systematic interaction map among functional DGs
Highly-linked drug groups
Number of links
Non-selective monoamine reuptake inhibitors
Vitamin K antagonists
Barbiturates and derivatives
Monoamine oxidase inhibitors, non-selective
Beta blocking agents, non-selective
Beta blocking agents, selective
BETA BLOCKING AGENTS
Tetracycline and derivatives
Selective serotonin reuptake inhibitors
On the other hand, the drugs in N03AB group which had the second largest degree were mostly related to many different cytochrome P450 mechanism (substrate, inhibitor, or inducer of CYP2C8, CYP2C19, CYP3A4, CYP2B6, CYP1A2, CYP2D6, and CYP2C9), which seems to be a potential common DDI mechanisms of DG N03AB. In addition, the drugs in J01MA such as fluoroquinolone (quinolone antibacterials) had DDIs with calcium, magnesium, zinc, and aluminium by formation of non-absorbable complexes.
Characteristics of DG-DG networks
Characteristic path length
Average number of neighbours
Secondary network of DGs sharing similar DG-DG interaction patterns
Network modules that consist of DGs sharing drug-interactions
Functional modules in drug-interaction sharing network
(Density × Number of Nodes)
Number of nodes
Number of edges
D10AA, R01AD, D07XA, D07XB, A01AC, L02AA, D07AB, D07AC, D07AA, S03BA
R03AA, N02CX, G02CA, S01EA, R03CB, R03AC, R03AB
S02BA, C05AA, S01CB, S01BA, H02AB
G04BA, A12AA, A12CC, B05XA
J01DC, J01DD, M03AC
D06AX, J01GB, S01AA
N03AA, N07AA, N05CA
The first module in Figure 4 (Cluster 1 in Table 3) contained the largest number of DGs, and, expectedly, most of them had the same therapeutic class which is corticosteroids. Although the majority of the module is composed of D07xx (xx denotes ambiguous ATC subclass), the other ATC classes L02AA (estrogens), A01AC (corticosteroids for local oral treatment), and S03BA (corticosteroids and anti-infectives in combination) also belonged to the similar therapeutic classes. The second module contained seven DGs. Most of them are sympathomimetics drugs that mimic the effects of transmitter substances of the sympathetic nervous system. Similar to the first module, the major DGs are R03xx (drugs for obstructive airway diseases) and the other DGs such as G02CA (sympathomimetics, labour repressants), S01EA (Sympathomimetics in glaucoma therapy), and N02CX (antimigraine preparations) also had the similar drug effects. In addition, their drug-interactions are mainly described by antagonism which can be an evidence of the similar molecular mechanisms of the DGs in this module. The module which consists of A12AA (calcium), A12CC (magnesium), B05XA (electrolyte solution), and G04BA (acidifier) showed the common drug-interaction forming non-absorbable complex, which might be due to the common ionization property.
Among the triangularly-linked modules, N05CA (barbiturates, antiepileptics), N07AA (anticholinesterases), and N03AA (barbiturates and derivatives, psycholeptics) commonly had an effect on nervous system but had different therapeutic subclasses. The other triangular groups, J01GB (other aminoglycoside antibacterials), D06AX (Other antibiotics for topical use), and S01AA (antibiotics, ophthalmologicals) were also commonly the antibiotics in spite of the different target organisms. Interestingly, two betalactam antibacterial drug-groups (J01DD and J01DC) were connected to muscle relaxants (M03AC). By literature search, we found that these two drug groups had the unexpected cross-reactivity.
In summary, the results suggested that information on DG-DG interaction partner sharing can be one useful feature to infer the mechanisms of DDI and therapeutic effects, and also, to reposition the existing drugs for another use.
DDI contains information on various drug features
From the analysis of the DG-DG network, we hypothesized that the drugs sharing common Drug-DG interaction partners might have similar drug mechanisms. To test this hypothesis, the Drug-DG partner sharing ratio was calculated for every drug pair to investigate common drug features between them. This sharing ratio was calculated by the same way as constructing Figure 4.
Instead of using Drug-Drug interaction, considering Drug-DG interactions in the calculation of the sharing ratio between drugs could have a beneficial effect for the reason that holistic method can compensate the noise caused by almost identical drugs.
In addition, features of each drug (e.g. prediction of drug target) were predicted by those of interaction-sharing drug (ISD) based on our hypothesis that an enriched common feature of ISDs is also likely to be a feature of the query drug. After collecting all ISDs, the features of the collected drugs were scored by frequency.
Simply, when the most frequent common feature for each drug was assigned to the feature of query drug, those of 80% were correctly predicted.
Drug interaction-sharing but dissimilar drugs
The relationship between chemical structure similarity and the feature-matching ratio for the interaction-sharing drugs was analysed with nine different drug features in Figure 5. Figures for all drug features were shown in Additional file 1 and Additional file 2. In Figure 6, feature-matching ratios (vertical axis) were plotted along the threshold of structural similarity (horizontal axis). For example, 0.5 in X-axis means that only the drug pairs whose structural similarity is less than 0.5 are only considered to calculate feature-matching ratios. 0.5 in Y-axis means that 50% of pairs share the same value of corresponding feature. It has to be noted that only the drugs which had the same type of annotations were considered for calculating the matching-ratio. The results showed that feature-matching ratio for most drug features were decreased as structurally dissimilar drugs were considered. Obviously 'phage1 metabolizing enzymes' seemed to be an important common factor for sharing drug-interactions regardless of structure similarity. However, interestingly, some drug features (e.g. 'general drug function', 'target reaction', and 'target pathway') maintained high matching ratio higher than 60% even when the structure similarity was around 0.4, which suggested that sharing drug-interactions contain not only information of pharmacodynamics such as 'target pathway' and 'target domain function', but also pharmacokinetic information such as metabolizing enzyme and transporter. In particular, 'target pathway' information showed greater than 80% matching ratio when the structure similarity was even ambiguous (0.4~0.5) while 'drug target' identity was very low.
On the other hand, the 'target essentiality' showed almost 100% matching ratio independent on structural similarity. Taken all together, we can conclude that the information on drug-drug interaction could be another unique source to describe various drug mechanisms with other chemical similarity measures.
Are DG-DG networks more useful than drug-drug interaction network?
In this work, we have showed that a DG-DG interaction network, its supplemental network and neighbour sharing information are intuitive and useful. This may misleads people thinking DG-DG interaction network gives more accurate information than DDI network does. Grouping drugs is a kind of abstraction, and the main purpose of our work is to show the concept of DG and to provide intuitive information. Thus, we expect complementary use of both DG-DG and drug-drug interactions will show more accurate and sensitive inference.
We used one ATC classification to bundle drugs, but this classification is not a unique method to categorize drugs. Our main concept can be applied with any kinds of drug classification (e.g. chemical type). More sophisticated way to define Drug Group will be able to improve DG-DG network and to provide novel kinds of information. We suggest that constructing DG-DG networks with different definitions of DG is new way to view novel dimension of drug interaction. Further studies are welcomed to validate.
In this study, we carried out systematic analysis on functional DG-DG interactions which provide more contextual and intuitive view about DDI. From the interaction map, we also constructed the secondary network which consists of DGs sharing drug interactions. The detected modules of the network represent the similar functional DGs in spite of different annotation of therapeutic class, suggesting that drug interactions contain the information about mechanisms of drugs. In addition the usefulness of our work in drug repositioning was shown with the example of betalactam antibacterial drugs and muscle relaxants. Moreover, we questioned whether DG-DG interaction information can be used to infer drug's mechanisms. The results show that some drug features such as metabolizing enzyme, drug function, and target pathway are well described by interactions even for the case that two drugs have low structural similarity. Thus, we expect that information on DG-DG interaction can be utilized as a novel useful source to describe drug's mechanisms.
Collection of drug-drug interaction and drug features
Descriptions on drug-drug interaction were retrieved from 'Drug_Interactions' records in DrugBank database . To define DDI pairs from the information, various drug names were unified by mapping synonyms, brand names generic name of each drug to a unique DrugBank ID. The descriptions that could not be mapped by DrugBank ID were discarded. As a result, 10,759 DDIs of 1,075 drugs were detected. In addition, various features of drugs and drug targets such as phage1 metabolizing enzyme, transporter, specific/general drug functions, drug category, ATC codes, drug reaction, affected organism; target information: drug target, pathway, cellular location, locus, domain function, and essentiality were also collected with the drugs.
Functional grouping of drugs by ATC code
Functional grouping of drugs was performed by mapping ATC codes into drugs. Anatomical Therapeutic Chemical (ATC) Classification System divides drugs into 5 different levels according to the organ or system on which they act and/or their therapeutic and chemical characteristics. One drug can have more than one code. Based on the drug classification system and its hierarchical structure, functional drug groups were defined, thus, each DG corresponds to each ATC code.
Statistical significance of DG-DG interactions
To calculate the statistical significance of the number of drug interactions between two ATC groups, ten thousands of random pairs of drug groups were generated for each unique combination of DG pairs and we counted how many drug interactions exists in the each random pair. Then p-value was calculated based on these distributions.
Interaction network between functional DGs of various levels
Significant DG-DG interactions (p-value < 0.0001) were considered to construct DG-DG interaction network. Note that ATC codes in different levels could be defined as a significant interaction because we tested all ATC codes regardless of the hierarchical structure. In addition, due to this hierarchical structure, some group interactions seemed to be redundant. For example, even though C01AA group was linked to both of C03BA and C03B, the interaction between C01AA and C03B just stemmed from more specific interaction between C01AA and C03BA. All of the redundant interactions were removed in constructing DG interaction map.
Interaction partner sharing ratio between DGs and between drugs
The interaction partner sharing ratio, how much interaction profile is similar, was measured between drug-groups or drugs. Specifically, drug groups (or drugs) were represented by the set of their interactive ATC groups (i.e. bit string). Based on the representation, all possible DG pairs (or drug pairs) were compared for measuring how much they are sharing the interactive ATC groups. The similarity measures used in this procedure were hyper geometric p-value measuring overlapping significance (p-value < 0.01) and Tanimoto coefficient (ration of common bits to union bits > 0.75). The DG pairs or drug pairs satisfying the two criteria were used to construct the secondary DG network based on the tendency sharing drug-interaction profiles.
Network module detection
Molecular Complex Detection (MCODE) was applied to identify densely connected modules in the secondary DG-DG network based on partner-sharing ratio. The method uses graph-theoretic clustering based on vertex weighting by local neighbourhood density and outward traversal from a locally dense seed node. Module score of each cluster was defined by the ratio of the number of edges to the number of nodes in each cluster.
Calculation of chemical structure similarity
Chemical structure similarity between drugs was based on 881-bit PubChem  fingerprints calculated by PaDEL software . Specifically, the fingerprint was calculated for each drug, and Tanimoto coefficient (ratio of intersection-bits to union-bits) was used as a 'chemical structural similarity' in this study.
We thank all members of Bioinformatics and Computational Biology Laboratory at KAIST for helpful discussions. This work was supported by a grant of the Korean Health Technology R&D Project, Ministry of Health & Welfare, Republic of Korea (A120017) and the National Research Foundation of Korea (NRF) grant funded by the Korea Government, the Ministry of Science, ICT & Future Planning (MSIP) [2008-2000220 & 2009-0086964].
The publication costs for this article were funded by a grant of the Korean Health Technology R&D Project, Ministry of Health & Welfare, Republic of Korea (A120017).
This article has been published as part of BMC Systems Biology Volume 7 Supplement 3, 2013: Twelfth International Conference on Bioinformatics (InCoB2013): Systems Biology. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcsystbiol/supplements/7/S3.
- Bender A, Scheiber J, Glick M, Davies JW, Azzaoui K, Hamon J, Urban L, Whitebread S, Jenkins JL: Analysis of pharmacology data and the prediction of adverse drug reactions and off-target effects from chemical structure. ChemMedChem. 2007, 2 (6): 861-873. 10.1002/cmdc.200700026.View ArticlePubMedGoogle Scholar
- Lee M, Kim D: Large-scale reverse docking profiles and their applications. BMC bioinformatics. 2012, 13 (Suppl 17): S6-PubMed CentralPubMedGoogle Scholar
- Yang L, Wang K, Chen J, Jegga AG, Luo H, Shi L, Wan C, Guo X, Qin S, He G: Exploring off-targets and off-systems for adverse drug reactions via chemical-protein interactome--clozapine-induced agranulocytosis as a case study. PLoS computational biology. 2011, 7 (3): e1002016-10.1371/journal.pcbi.1002016.PubMed CentralView ArticlePubMedGoogle Scholar
- Lee J, Xia Y, Son MY, Jin G, Seol B, Kim MJ, Son MJ, Do M, Lee M, Kim D: A novel small molecule facilitates the reprogramming of human somatic cells into a pluripotent state and supports the maintenance of an undifferentiated state of human pluripotent stem cells. Angewandte Chemie. 2012, 51 (50): 12509-12513. 10.1002/anie.201206691.View ArticlePubMedGoogle Scholar
- Lee M, Han S, Chang H, Kwak YS, Weller DM, Kim D: FitSearch: a robust way to interpret a yeast fitness profile in terms of drug's mode-of-action. BMC genomics. 2013, 14 (Suppl 1): S6-PubMed CentralPubMedGoogle Scholar
- Parsons AB, Lopez A, Givoni IE, Williams DE, Gray CA, Porter J, Chua G, Sopko R, Brost RL, Ho CH: Exploring the mode-of-action of bioactive compounds by chemical-genetic profiling in yeast. Cell. 2006, 126 (3): 611-625. 10.1016/j.cell.2006.06.040.View ArticlePubMedGoogle Scholar
- Han S, Lee M, Chang H, Nam M, Park HO, Kwak YS, Ha HJ, Kim D, Hwang SO, Hoe KL: Construction of the first compendium of chemical-genetic profiles in the fission yeast Schizosaccharomyces pombe and comparative compendium approach. Biochemical and biophysical research communications. 2013Google Scholar
- Lamb J, Crawford ED, Peck D, Modell JW, Blat IC, Wrobel MJ, Lerner J, Brunet JP, Subramanian A, Ross KN: The Connectivity Map: using gene-expression signatures to connect small molecules, genes, and disease. Science. 2006, 313 (5795): 1929-1935. 10.1126/science.1132939.View ArticlePubMedGoogle Scholar
- Campillos M, Kuhn M, Gavin AC, Jensen LJ, Bork P: Drug target identification using side-effect similarity. Science. 2008, 321 (5886): 263-266. 10.1126/science.1158140.View ArticlePubMedGoogle Scholar
- Morphy R: Selectively nonselective kinase inhibition: striking the right balance. Journal of medicinal chemistry. 2010, 53 (4): 1413-1437. 10.1021/jm901132v.View ArticlePubMedGoogle Scholar
- Chou TC: Theoretical basis, experimental design, and computerized simulation of synergism and antagonism in drug combination studies. Pharmacological reviews. 2006, 58 (3): 621-681. 10.1124/pr.58.3.10.View ArticlePubMedGoogle Scholar
- Bai JP, Abernethy DR: Systems pharmacology to predict drug toxicity: integration across levels of biological organization. Annual review of pharmacology and toxicology. 2013, 53: 451-473. 10.1146/annurev-pharmtox-011112-140248.View ArticlePubMedGoogle Scholar
- Kohler GI, Bode-Boger SM, Busse R, Hoopmann M, Welte T, Boger RH: Drug-drug interactions in medical patients: effects of in-hospital treatment and relation to multiple drug use. International journal of clinical pharmacology and therapeutics. 2000, 38 (11): 504-513. 10.5414/CPP38504.View ArticlePubMedGoogle Scholar
- Kuhlman JJ, McCaulley R, Valouch TJ, Behonick GS: Fentanyl use, misuse, and abuse: a summary of 23 postmortem cases. Journal of analytical toxicology. 2003, 27 (7): 499-504. 10.1093/jat/27.7.499.View ArticlePubMedGoogle Scholar
- Piecoro LT, Browning SR, Prince TS, Ranz TT, Scutchfield FD: A database analysis of potentially inappropriate drug use in an elderly medicaid population. Pharmacotherapy. 2000, 20 (2): 221-228. 10.1592/phco.184.108.40.206779.View ArticlePubMedGoogle Scholar
- Forrester MB: Ingestions of hydrocodone, carisoprodol, and alprazolam in combination reported to Texas poison centers. Journal of addictive diseases. 2011, 30 (2): 110-115. 10.1080/10550887.2011.554778.View ArticlePubMedGoogle Scholar
- Buajordet I, Ebbesen J, Erikssen J, Brors O, Hilberg T: Fatal adverse drug events: the paradox of drug treatment. Journal of internal medicine. 2001, 250 (4): 327-341. 10.1046/j.1365-2796.2001.00892.x.View ArticlePubMedGoogle Scholar
- Percha B, Altman RB: Informatics confronts drug-drug interactions. Trends in pharmacological sciences. 2013, 34 (3): 178-184. 10.1016/j.tips.2013.01.006.View ArticlePubMedGoogle Scholar
- Schadt EE, Friend SH, Shaywitz DA: A network view of disease and compound screening. Nature reviews Drug discovery. 2009, 8 (4): 286-295. 10.1038/nrd2826.View ArticlePubMedGoogle Scholar
- Dudley JT, Schadt E, Sirota M, Butte AJ, Ashley E: Drug discovery in a multidimensional world: systems, patterns, and networks. Journal of cardiovascular translational research. 2010, 3 (5): 438-447. 10.1007/s12265-010-9214-6.View ArticlePubMedGoogle Scholar
- Huang J, Niu C, Green CD, Yang L, Mei H, Han JD: Systematic prediction of pharmacodynamic drug-drug interactions through protein-protein-interaction network. PLoS Comput Biol. 2013, 9 (3): e1002998-10.1371/journal.pcbi.1002998.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang LC, Wu X, Chen JY: Predicting adverse drug reaction profiles by integrating protein interaction networks with drug structures. Proteomics. 2013, 13 (2): 313-324. 10.1002/pmic.201200337.View ArticlePubMedGoogle Scholar
- Azuaje F: Drug interaction networks: an introduction to translational and clinical applications. Cardiovascular research. 2013, 97 (4): 631-641. 10.1093/cvr/cvs289.View ArticlePubMedGoogle Scholar
- Hu TM, Hayton WL: Architecture of the drug-drug interaction network. J Clin Pharm Ther. 2011, 36 (2): 135-143. 10.1111/j.1365-2710.2009.01103.x.View ArticlePubMedGoogle Scholar
- Xu KJ, Song J, Zhao XM: The drug cocktail network. BMC systems biology. 2012, 6 (Suppl 1): S5-10.1186/1752-0509-6-S1-S5.PubMed CentralView ArticlePubMedGoogle Scholar
- WHO Collaborating Center for Drug Statistics Methodology: The anatomical therapeutic chemical (ATC) classification system. Guidelines for ATC Classification and DDD Assignment. 1998, Oslo, 12-29.Google Scholar
- Bader GD, Hogue CW: An automated method for finding molecular complexes in large protein interaction networks. BMC bioinformatics. 2003, 4: 2-10.1186/1471-2105-4-2.PubMed CentralView ArticlePubMedGoogle Scholar
- Knox C, Law V, Jewison T, Liu P, Ly S, Frolkis A, Pon A, Banco K, Mak C, Neveu V: DrugBank 3.0: a comprehensive resource for 'omics' research on drugs. Nucleic Acids Res. 2011, 39 (Database): D1035-1041. 10.1093/nar/gkq1126.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Y, Xiao J, Suzek TO, Zhang J, Wang J, Bryant SH: PubChem: a public information system for analyzing bioactivities of small molecules. Nucleic Acids Res. 2009, 37 (Web Server): W623-633. 10.1093/nar/gkp456.PubMed CentralView ArticlePubMedGoogle Scholar
- Yap CW: PaDEL-descriptor: an open source software to calculate molecular descriptors and fingerprints. J Comput Chem. 2011, 32 (7): 1466-1474. 10.1002/jcc.21707.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.