PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2017; 7: 11401.
Published online 2017 September 12. doi:  10.1038/s41598-017-11924-4
PMCID: PMC5595859

From malaria to cancer: Computational drug repositioning of amodiaquine using PLIP interaction patterns

Abstract

Drug repositioning identifies new indications for known drugs. Here we report repositioning of the malaria drug amodiaquine as a potential anti-cancer agent. While most repositioning efforts emerge through serendipity, we have devised a computational approach, which exploits interaction patterns shared between compounds. As a test case, we took the anti-viral drug brivudine (BVDU), which also has anti-cancer activity, and defined ten interaction patterns using our tool PLIP. These patterns characterise BVDU’s interaction with its target s. Using PLIP we performed an in silico screen of all structural data currently available and identified the FDA approved malaria drug amodiaquine as a promising repositioning candidate. We validated our prediction by showing that amodiaquine suppresses chemoresistance in a multiple myeloma cancer cell line by inhibiting the chaperone function of the cancer target Hsp27. This work proves that PLIP interaction patterns are viable tools for computational repositioning and can provide search query information from a given drug and its target to identify structurally unrelated candidates, including drugs approved by the FDA, with a known safety and pharmacology profile. This approach has the potential to reduce costs and risks in drug development by predicting novel indications for known drugs and drug candidates.

Introduction

Pharmaceutical companies spend around $ 2.6 billion in developing a drug through to market approval1. To minimize risk and development time, drug repositioning moves experimental or approved drugs to new indications, so that data from previously conducted safety and pharmacology studies can be leveraged. There are several examples of repositioning success stories, such as sildenafil (Viagra), which was originally developed for heart disease and was repurposed for erectile dysfunction, the sedative thalidomide, which is now approved for treatment of multiple myeloma and leprosy2, or the cytotoxic anti-cancer agent gemcitabine, which was originally developed as an anti-viral. The link between anti-viral and anti-cancer effects is also demonstrated with the small molecule brivudine (BVDU), which is approved for treatment of herpes and which has been investigated for use in pancreatic cancer3.

It appears surprising that one drug should be a cure for two diseases. One explanation is that a drug can bind promiscuously, i.e. to multiple different targets. In previous work, we established that drug promiscuity correlates with shared binding sites across the drug’s multiple targets4. Thus, structural analyses of shared binding sites and drug-target interactions are promising approaches to drug repositioning.

Such analyses hinge on the availability of structural data. While structural data is not as abundant as sequence data, it is growing steadily. The Protein Data Bank (PDB) has more than doubled in size in the last seven years. Today, it contains 3D structures of over 1,200 different drug targets5 and more than 60% of all PDB structures contain proteins complexed with biologically relevant ligands4. The growing availability of data is complemented by an increasing number of tools and methods mostly focusing on binding pockets68 or ligands911. A third approach characterises the interaction of ligands and binding pockets. Here, we recently introduced the Protein-Ligand Interaction Profiler (PLIP)12, a tool for comprehensive detection of molecular contacts. In this paper, we show how an analysis of a known drug-target interaction with PLIP can define interaction patterns, which can then be run against the PDB (Fig. 1). Since the patterns are structure-invariant, the screen against PDB will reveal ligands with novel scaffolds and novel targets. To document the power of this novel in silico screening approach, we tested it on the herpes drug BVDU. After defining the patterns through which BVDU interacts with its target proteins, we screened for compounds matching these patterns and then tested the hits in vitro for their potency in inducing cell killing in cultured cancer cells and inhibiting the function of a drug resistance target (anti-cancer effect).

Figure 1
PLIP Screening. Screening for BVDU interaction patterns reduces an initial set with 13,704 compounds to 58 candidate compounds.

Results

In the 1980s, BVDU was first introduced as a treatment for Herpes zoster infection. BVDU is a thymidine analogue, composed of a nucleobase ring, a deoxyribose moiety, and a bromovinyl residue. In the infected cell, BVDU is phosphorylated by a viral thymidine kinase and then erroneously integrated into the viral genome13, terminating the viruses ability to replicate.

BVDU interaction patterns for PLIP screening

For our approach to structural drug repositioning it is important that crystal structure data for BVDU in complex with viral and non-viral kinases14 is available in PDB. Figure 2 shows in detail, how BVDU interacts with these kinases: BVDU’s nucleobase ring (A) engages in π-stacking (B–D) and features two prominent parallel hydrogen bonds (E–G). Its deoxyribose moiety has a hydrophobic contact (H) and a distal hydrogen bond (I) and the bromovinyl residue forms a halogen bond (J).

Figure 2
BVDU interaction patterns. Aromatic ring for π-stacking (A), double π-stacking to base ring (B), on opposite sites (C), in parallel (D), double hydrogen bonds to base ring (E), to the same residue (F), in parallel (G), distal hydrophobic ...

The ten interaction patterns (A–J) were derived from five different PDB structures originating from viruses, human, and fruitfly. Given the variety of species it can be expected that not all patterns are present in all complexes. Even within the same species, there might be different use of the patterns, since not all interactions may be essential. Figure 3 demonstrates this point by showing presence and absence of the patterns in five PDB structures. A single π-stacking to BVDU’s aromatic ring (A), double hydrogen bonds to the base ring (E) and the distal hydrophobic contact (H) are present in all structures, whereas the parallel hydrogen bonds (F, G) are present in nearly all, and the double π-stacking patterns (B–D) only in some structures. Interestingly, the halogen bond (J) is not present in all viral structures. Overall, analysis with PLIP reveals some distinct features and combinations among the ten patterns, and we screened the PDB for these ten patterns, ranking compounds by the number of matches.

Figure 3
BVDU interactions patterns from five target proteins. Ten patterns relating to π-stacking (AD), parallel hydrogen bonds (EG), hydrophobic contacts (H), hydrogen bond (I) and a halogen bond (J) and their presence in five targets ...

PLIP screening of PDB for BVDU interaction patterns

Our goal was to identify compounds which can form interactions in a similar manner to BVDU. We postulated that since these new compounds share interaction patterns with BVDU, they might also bind to the same target protein in cancer as BVDU. Screening the entire PDB for any ligands and targets, which interact similarly to BVDU, is not a trivial task as PDB comprises some 107,000 structures. After removing biologically irrelevant entries with BioLiP (see Methods), there are 170,219 complexes between 16,460 unique targets and 13,704 unique ligands. Each of these 170,219 complexes was analysed for the presence of the ten BVDU patterns and the results are shown in Fig. 4. There is no complex which contains all ten interaction patterns, but thymidine, the parent scaffold from which BVDU is derived, is among six ligands that contain nine patterns. Compounds were ranked by the number of patterns and annotated whether or not they were approved for therapy by the FDA. There are just 951 complexes containing six or more patterns and these complexes contain fewer than 250 compounds including 12 approved and 46 candidate drugs (i.e. experimental drugs or drugs in currently tested in clinical trials).

Figure 4
BVDU patterns effectively identify repositioning candidates. There are 247 compounds with six or more BVDU patterns including 12 approved drugs. One such drug is amodiaquine, an anti-malaria agent. As a control, thymidine, a BVDU analogue, satisfies nine ...

Prioritization of repurposing candidates

Three criteria were used to prioritize candidate compounds: Firstly, we sought compounds with a different scaffold to BVDU. Secondly, we used FDA-approval to prioritise the results, since these compounds will have been extensively researched and information on safety pharmacology, dosing and formulation may be available. Such information can significantly increase the value of potential repositioning candidates. Thirdly, we sought an improvement in potency over BVDU, meaning the drug should show a stronger inhibition of BVDU’s target in cancer than BVDU itself.

Regarding scaffold diversity, we expected many of the top-scoring hit compounds to be nucleobases, since BVDU is a thymidine analogue. As mentioned above, thymidine is ranked highest with nine matching patterns and thus serves as a positive control. Similarly, deoxy-adenosine, -guanosine, and -cytidine contain seven BVDU patterns and deoxy-uridine contains six patterns. Overall, the 58 drugs with six or more patterns contain 40 such nucleobases. In Fig. 5 the 58 drugs are grouped according to chemical similarity and this reveals two large clusters (red boxes), which comprise pyrimidine and purine scaffolds. The remaining groups are much smaller in size and comprise scaffolds different from BVDU. These 18 novel compounds are of particular interest.

Figure 5
Scaffold hopping from BVDU to amodiaquine. Pairwise chemical similarity of compounds (red = similar, yellow = dissimilar) arranged by chemical scaffold and clustered by chemical similarity within each group. The top group ...

The second priorization aspect is FDA-approval: Among the 58 drugs, 12 are FDA-approved drugs which contain between 6 and 8 BVDU patterns (see Fig. 6 and Table 1). Half of these 12 FDA-approved drugs have a different scaffold to BVDU, namely the aminopteridones (folic acid, tetrahydrobiopterin, triamterene), the flavone quercetin as well as the singletons pyridoxal phosphate and amodiaquine.

Figure 6
Approved drugs with six or more patterns. Twelve FDA-approved drugs across very different diseases and targets (panel A). Six are nucleobases like BVDU, but six are novel scaffolds. The patterns for π-stacking (AD) and double hydrogen ...
Table 1
Approved drugs with six or more patterns.

Finally, we evaluated the potential of hit compounds to bind to BVDU’s primary target Herpes simplex thymidine kinase, the starting point for our pattern-based screening. Heinrich et al.3 showed that BVDU also binds the human heat shock protein, Hsp273, and thereby blocks its anti-apoptotic activity allowing cytotoxic agents to reestablish their efficacy. To assess whether candidates could bind to these targets, we employed in silico docking with the widely used Autodock software (see Methods for details). We docked all 58 drugs to the relevant binding sites in the viral thymidine kinase and the a homology model of Hsp273 and found that chrysin, olaparib, amodiaquine, and acenaphthenoquinone ranked better than BVDU in both cases (see Supplementary Fig. 1 for details). Figure 7 shows how the software docked amodiaquine to the two targets and how the docking poses satisfy the π-stacking patterns A–D as well as the distal hydrophobic contact pattern. However, they do not satisfy the hydrogen bond patterns E–G. In this respect, amodiaquine binds the BVDU targets in the same manner as one of its targets, histamine methyltransferase.

Figure 7
Amodiaquine in complex with its target histamine methyltransferase (2AOU, left), with the in silico docked target thymidine kinase (1OSN, middle) and Hsp27 model (right). All contain the patterns A–D and H.

As a final step to test our structural repurposing approach with PLIP, we tested amodiaquine in vitro in biochemical and cell biological assays.

Amodiaquine as anti-cancer drug?

BVDU exerts its anti-cancer effect by suppressing the cancer cells’ ability to develop resistance to a chemotherapy. Heinrich et al.14 demonstrated this by treating a multiple myeloma cell line with bortezomib, a chemotherapeutic agent, at increasing doses, which creates a selective pressure for the cells to develop resistance. Co-treatment with bortezomib and BVDU14 inhibits cell growth effectively and allows bortezomib to reestablish its cytotoxic effect. They also show that BVDU achieves this effect only in combination with bortezomib and has no cytotoxic activity on its own. We tested amodiaquine in the same chemoresistance assay to document its anti-cancer potential (Fig. 8). The cancer cells were exposed to increasing doses of bortezomib over three cell culture passages (see Fig. 9 for details). At the beginning of the third passage on day 12, 100,000 cells are exposed to the treatment regimes (bortezomib on its own and bortezomib with amodiaquine). After one week, resistance to bortezomib is clearly shown as cells continue to grow and have multiplied eightfold to 800,000. In contrast, the co-treatment with amodiaquine leads to a significant reduction in cell number to 80,000 (10%).

Figure 8
Drug-resistant multiple myeloma cells continue to grow despite treatment with the cytostatic drug bortezomib. However, when combined with bortezomib, amodiaquine significantly reduces cell numbers. Chrysin achieves a similar effect.
Figure 9
Detailed results of the cell resistance assay. (a) Amodiaquine re-establishes chemosensitivity in resistant cancer cells. Cells are treated with increasing doses of the cytotoxic agent bortezomib and amodiaquine over three passages. There is a clear difference ...

Amodiaquine inhibits Hsp27’s chaperone function

In a biochemical assay, which elucidates how compounds inhibit the chaperone function of the heat shock protein Hsp27, Heinrich et al.3 showed that BVDU is a potent inhibitor. Straume et al.15 showed that breast cancer cells expressing Hsp27 become resistant to chemotherapy and that Hsp27 knockdown reestablished susceptibility. Furthermore, Bruey et al.16 showed that Hsp27 negatively regulates apoptosis via interaction with cytochrome C. Following this line of thought, Heinrich et al.14 devised an assay to show how BVDU impacts Hsp27’s function as a chaperone. Citrate synthase is used as a client protein of Hsp27 and misfolds at 43 °C, but in the presence of Hsp27 misfolding is inhibited. The chaperone activity of Hsp27 can be measured by the amount of misfolded client protein determined by capillary electrophoresis of precipitated protein. We tested amodiaquine in the Hsp27 chaperone assay and found that amodiaquine significantly inhibits Hsp27’s activity (see Methods for details) and that it was 43 more potent than BVDU. This activity is within the range of compounds tested in a previous study by Heinrich et al.14. Table 2 shows the results of the chaperone assay for the two tested compounds in comparison to lead compounds from the study by Heinrich et al.14.

Table 2
Results of the chaperone assay for amodiaquine, chrysin, and lead compounds from study by Heinrich et al.14 in comparison to the control BVDU.

Repositioning matrix

Amodiaquine shares key interaction patterns with BVDU and also binds Hsp27. But to what extend is the ability to form these interaction patterns deterministic for target protein binding? Does amodiaquine also bind to the Herpes thymidine kinase and does BVDU bind to amodiaquine’s target histamine methyltransferase? In fact, do all of the 58 drugs with six or more BVDU interaction patterns bind to all of the 43 PDB targets? It is unlikely that this holds in its entirety, but it is reasonable to assume that some of these potential “repositioning cases” are valid. A systematic in vitro validation of all 58 drugs against all 43 targets was not carried out, however, a survey of the literature and databases supports many of these predicted interactions, and we grouped these drugs into 19 drug classes and the 43 targets into 23 superfamilies (see Methods). Figure 10 reveals whether there is evidence for target binding from crystal structures (green) or from other sources such as text mining, biological assays, or pathways (orange).

Figure 10
Evidence for binding of 58 drugs to 43 targets (grouped by scaffold type and superfamily). Green indicates structural evidence, orange other evidence. For 28% of pairs there is binding evidence.

From all 437 possible binding combinations of proteins and compounds (Fig. 10), 41 (9%) could be confirmed by structural evidence (green). For 81 (18.5%) other drug-target pairings, there is evidence for binding from text mining, bioassays, screening or pathway data. Taken together, both sources provide binding evidence for 28% of all possible combinations. By arranging the data points by the number of shared targets (for compounds) and the number of shared ligands (for targets), clusters were produced in the matrix, which indicate that for some protein and ligand ensembles, evidence is available that they may bind to each other. One small distinct cluster (Fig. 10, lower right corner) is formed by dUTPase-like proteins, UDP-Glycosyltransferases, HAD-like protein, Nudix-, and (trans)glycosidases, together with nucleosides, nucleotides and singletons are supported by both structural and other evidence. A large cluster is visible in the top right of the matrix with 6 compound and 7 protein groups. Interestingly, this cluster also contains the chemical group of the hit compound amodiaquine and its target histamine methyltransferase (blue) as well as BVDU’s chemical group and the superfamily of BVDU’s target thymidine kinase (both marked in green font). The available evidence places the nucleosides as the compound group, which binds to the largest number of target superfamilies (17), while the P-loop containing nucleoside triphosphate hydrolases are known to bind to most (12) of the chemical classes from the set of 58 drugs.

Validation on a different scaffold: chrysin

Out of interest, we took a closer look at the highest-ranking compound not in the list of FDA-approved drugs to experimentally confirm another binding prediction with a different chemical scaffold (Fig. 10). Evidence supporting the potential for binding of the flavone chrysin to Hsp27 is particularly strong: Chrysin’s interaction with rabbit glycogen phosphorylase (PDB ID 3EBO) satisfies the BVDU consensus patterns A (base ring), E (double hydrogen bonds), and H (distal hydrophobic contact), as well as the patterns B-D for π-stacking. It also was ranked best in the in silico docking experiment. Thus, we tested chrysin in the chemoresistance assay (see Methods) and Fig. 8 shows that, similar to amodiaquine, chrysin in combination with bortezomib reestablishes cytostatic activity and reduces cell growth by 50% in comparison to cells only with bortezomib alone. Growth of the drug-resistant cells is not completely stopped (as with amodiaquine), but is substantially reduced. In the Hsp27 chaperone assay, chrysin performs 69 times better than the reference compound BVDU (see Methods and Table 2).

Discussion

The results we present in this work spark discussion on three points. Firstly, how does PLIP screening for drug repositioning compare to other computational approaches? Secondly, could the discovery of amodiaquine’s anti-cancer effect have been obtained differently? Thirdly, can the PLIP approach be applied to other indications?

PLIP utilizes interaction-based screening, in contrast to ligand- and target-based approaches, which are based on chemical or protein structure, respectively. Ligand-focused approaches811, 17, 18 usually consider molecular structure properties or fragment composition to find similar compounds. This is successful if the goal is the discovery of structurally related compounds. However, to identify novel scaffolds (ligand hopping), this method is not suitable. In our example, amodiaquine is structurally unrelated to BVDU and could therefore not have been found by these methods. Target-focused methods6, 7, 1923 rely on either geometrical analysis or chemical analysis of the target protein. Similar to ligand hopping, these methods do not lend themselves to target hopping. As an example, BVDU binds a thymidine kinase and amodiaquine a methyltransferase. These targets do not share any similarities, however, by focussing on the drug-protein interactions, PLIP can find a relationship between BVDU/thymidine kinase and amodiaquine/methyltransferase, since they interact in a very similar manner. This focus on interactions rather than on ligands or target is also implemented in other tools24, 25, however, these do not consider some interaction types (π-stacking and/or halogen bonds), which were necessary to define the BVDU patterns and finally to make the link to amodiaquine and other hit compounds.

Based on drug-target interaction patterns, PLIP identified the anti-cancer potential of the anti-malaria drug amodiaquine. We verified this anti-cancer effect by showing that multiple myeloma cancer cells, which have developed resistance to chemotherapy, become sensitized to chemotherapy again upon exposure to amodiaquine. One important question here is whether amodiaquine is also a cytotoxic agent or works via a specific mechanism on cell resistance. The resistance assay clearly shows that amodiaquine has a specific and synergistic effect with the cytotoxic agent bortezomib (see Fig. 9). Without the cytotoxic agent (control runs a + b), the cells can grow unhindered, so a cell toxicity of amodiaquine can be excluded. When administered in combination with bortezomib, however, the growth behaviour clearly changes due to a re-establishing of chemosensitivity by amodiaquine. The same is true for chrysin, although the effect on cell growth is not as strong.

We show that amodiaquine inhibits the chaperone function of the heat shock protein Hsp27, which plays a key role in the chemoresistance of some tumor types15. These findings are supported by Qiao et al.26, who treated cultured malignant melanoma cells with amodiaquine and found that it sensitized them to starvation- and chemotherapeutic-induced death. They also showed that heat shock proteins, Hsp70 and Hsp90, play an important role in this process, similar to Hsp27 in our studies.

Unlike most repositioned drugs, Qiao et al.26 used a rational approach to identify amodiaquine. The authors argued that autophagy and lysosomal degredation play an important role in cancer27, 28 and that, as many anti-malaria drugs are lysosomotropic, they set out to screen anti-malarial agents for their anti-cancer behaviour. Their approach to drug repositioning is nonetheless serendipitous, since the hypothesis was not generated systematically. However, it also shows how techniques such as text-mining, reasoning, mechanism of action insights, and database integration will play an important role in the systematization of drug repositioning.

In contrast to Qiao et al.26, PLIP’s identification of amodiaqine was based on an algorithmic processing of data on its target, methyltransferase, which shares its binding site interactions with BVDU’s targets Hsp27 and thymidine kinase. We obtained proof that amodiaquine inhibits Hsp27 function, but binding to thymidine kinase and any anti-viral activities were not tested. However, there is evidence that amodiaquine inhibits virus replication of flaviviridae including dengue virus29 and there is also evidence that it interferes with thymidine production30. Such circumstantial evidence supports future experiments to test amodiaquine’s anti-herpetic activity.

PLIP’s ranking of compounds by the number of matching interaction patterns was crucial and very efficient: 13,000 compounds in the PDB are reduced to just 247 compounds containing at least six BVDU patterns, including a mere 12 FDA-approved drugs, one being amodiaquine. This drug shows a double π-stacking sandwich, which is highly specific (0.5% of interactions in PDB), a hydrophobic contact and a halogen bond, but is missing the parallel hydrogen bonds. However, the latter is a consensus interaction pattern among the five structures (Fig. 3), which means that a search for this consensus would not have included amodiaquine. Thus, we believe, that a ranking by number of patterns is more advantageous than filtering by a defined set of patterns. Ranking by number of patterns becomes even more powerful if combined with a chemical similarity analysis to exclude structurally obvious hits and focus on compounds with novel scaffolds. In our study, 58 hit compounds (Fig. 4) contain 40 nucleobases, which are chemically similar and hence not novel, leaving only 18 compounds with a good match in interaction patterns, but highly dissimilar to BVDU and hence, are novel scaffolds.

In our previous work (Heinrich et al.14) we identified novel Hsp27 inhibitors by target hopping from thymidine kinase to Hsp27. We argued that if one thymidine kinase binder (BVDU) inhibits Hsp27, then so may others, and identified 249 thymidine kinase binders, that were reduced to 29 with better in silico binding affinity than BVDU. Six of these 29 were validated in vitro in the same chemoresistance and chaperone function assays used here for amodiaquine. Both approaches largely complement each other, since in our earlier work, we did not require structural data to identify the 249 compounds. Hence, there is only a small overlap of 10 compounds with our candidates. However, the success of our work14 means that all of the 43 targets we identified with PLIP can serve as a starting point to collect binders, which could then be ranked and experimentally validated.

Next to amodiaquine from our FDA-approved hit candidates, we selected chrysin, a natural compound, as a second candidate for validation, since the interaction pattern score from the PLIP screen was high. Result from the in vitro chemoresistance and chaperone assays were good, although not as potent as amodiaquine. Chrysin has shown beneficial effects in cancer cell lines3134 and additionally, a moderate antiviral activity was demonstrated against H. simplex, V. zoster, and human Enterovirus 713539. These findings are consistent with the hypothesis that chrysin binds Hsp27 and the herpes thymidine kinase.

However, chrysin also has anti-allergic and anti-inflammatory activity33 and contains a catechol group with promiscuous binding behaviour40. Many potential targets have been proposed34 and there is also structural evidence that its binding can be very different from the patterns discussed above. For example, chrysin binds to the transport protein transthyretin (PDB ID 4DES) with low affinity41 using almost exclusively hydrophobic contacts and just one hydrogen bond to a lysine at the entrance of the binding site. The in vitro effects of chrysin require very high micromolar concentrations4244, so while chrysin is an interesting compound, whose binding behaviour supports our structural approach to drug repositioning, in contrast to amodiaquine, chrysin’s pharmacological profile does not support further development as an anti-cancer lead substance.

To estimate how widely applicable PLIP screening is, consider the Repurposed Drug Database45 as a reference. This lists 233 repositioned drugs and for 56 of these (24%), there is structural data for the ligand and one therapeutic target. This means that for 24% of the repurposed drugs, there is similar coverage of structural data as there was for our case using BVDU patterns. Similarly, there is structural data for more than 25% of FDA-approved drugs (400 out of ca. 160046). Hence, while there is far less structural data than, e.g. sequence data, there is nonetheless a substantial amount, which is sufficient for structural approaches, such as PLIP, to be a viable tool beyond examples presented in this paper.

Concluding, we could show that our structure-based screening approach with PLIP interaction patterns identifies novel repositioning candidates for the cancer target Hsp27. Not only did we identify candidates structurally unrelated to the query drug BVDU (scaffold hopping), but also demonstrated that they show the desired inhibitory activity on the target protein Hsp27 and suppress chemoresistance. Especially the FDA-approved malaria drug amodiaquine, which emerged as the top hit from our screen and was subsequently validated, proves the potential of our approach for drug repositioning.

Methods

Interaction patterns

107,663 structures were downloaded from the Protein Data Bank (PDB)47 FTP Archive on Apr 2015. Each structure was analysed with PLIP v1.1.112 using standard settings, resulting in 408,877 complexes. These were filtered using BioLiP, a database for discrimination between biologically relevant ligand and artifacts48, to keep 170,219 complexes.

To construct BVDU interaction patterns, the PDB structures 2W0S, 2JAW, 1OSN, 2VQS, 1KI8 (which contain BVDU) were analysed with PLIP, resulting in fourteen complexes overall. The structure 4XSC was not available at that time (released date Dec 2016), but his absence does not imply any significant change in the analysis. Pattern definitions include one or two specific non-covalent interactions and their relative orientation (angle or distance ranges). The ten pattern are shown in Fig. 2 and their definitions listed in Table 3.

Table 3
Definitions of patterns A–J.

The 170,219 complexes were screened for the presence of the ten patterns and if a ligand was represented by several complexes, the complex with the most BVDU patterns was selected. Figure 4 shows the distribution of pattern numbers for these complexes. Of these complexes, 247 have six or more patterns and 58 of these are drugs, classified as experimental, clinical, or approved drugs based on the Therapeutic Target Database (TTD)49.

Chemical similarity

To prioritize the 58 compounds they were manually assigned to scaffold types and their pairwise chemical similarity was computed using the PubChem Score Matrix Service (pubchem.ncbi.nlm.nih.gov/score_matrix) with standard settings. Chemicals were downloaded from PubChem in SDF format. A heatmap of the pairwise scores was generated with the Heatplus package v2.16.0 in R 3.2.5. Compounds on the axes were grouped by scaffold and within each scaffold by hierarchical clustering on average distances.

Docking

Docking was performed with AutoDock 4.2 rigid body docking on the known BVDU binding site in 1OSN and in the model of Hsp27, as described by Heinrich et al.3. Structures were prepared with AutoDockTools v1.5.4 to assign atom types and partial charges. The docking area was defined in the respective binding pockets by a box of 60 × 60 × 60 Å for 1OSN and 50 × 60 × 50 Å for Hsp27 with 0.375 Å spacing. Lamarckian genetic algorithm with 150 randomly placed entities, 27,000 generations, 5,000,000 energy evaluations, a mutation rate of 0.02, elitism value 1, and a cross-over rate of 0.80 was used for the docking process, with a total of thirty runs per compound. The local search was performed using the Soils and Wets algorithm with 300 iterations per search. After the docking, pose clusters were generated using Root Mean Square Deviation (RMSD) values.

Chemoresistance assay

BVDU was a gift from Rudolf Fahrig, RESprotect GmbH, Dresden (Germany). Amodiaquine and chrysin were bought from Sigma-Aldrich. A multiple myeloma cell line (RPMI-8226) was obtained from the DSMZ (German Collection of Microorganisms and Cell Cultures). The cells were cultivated in RPMI 1640 medium, supplemented with 10% (v/v) fetal bovine serum in a humidified atmosphere at 37 °C and 5% CO 2. Cells in logarithmic growth phase were reseeded at a density of 100,000 cells per mL and incubated with the chemotherapeutic agent bortezomib (Velcade) together with the test compounds. Bortezomib was omitted for the control run. Cells were passaged regularly to prevent densities of more than 1,000,000 cells per mL. Bortezomib was added at an initial dose of 0.1 ng/mL, followed by increasing doses of 0.2 ng/mL (second passage) and 0.3 ng/mL (3rd passage). As a positive control, cells were incubated with bortezomib and 30 μM BVDU. For all tested compounds, non-toxic doses were previously determined.

Hsp27 chaperone function assay

Hsp27 is a chaperone for many client proteins including citrate synthase (CS). CS misfolds at temperatures above 43 °C, but in the presence of Hsp27, misfolding is reduced. Inhibitors binding to Hsp27 impair its chaperone function, leading to accumulation of denatured CS and the amount of misfolded protein can thus serve as a measure of Hsp27 inhibition. BVDU, controls, or test compounds were incubated in 40 mmol L-HEPES buffer with 1.44 μM CS and 481 nM Hsp27 at 43 °C and pH 7.4. Samples were taken after 30, 60 and 90 min and aggregated CS quantified using capillary electrophoresis. The relative inhibitory values were calculated by measuring the amount of misfolded CS and subtracting the value of the blank (DMSO buffer) and all values were normalized to the value for BVDU and the concentration of test compounds used.

Repositioning matrix

We constructed a “repositioning matrix” of the 58 compounds, satisfying six or more patterns, and their 43 targets. For each possible compound-target pair we reviewed the databases STITCH, Chembl, and PDB for binding evidence. STITCH v5.050 was accessed in Feb 2017, ChEMBL51 in Aug. 2016 and PDB47 in Apr. 2015. STITCH target IDs were mapped through the file full_uniprot_2_string.04_2015.tsv.gz. ChEMBL compounds were mapped by PubChem, ChEMBL targets by name.

Since the data points in the repositioning matrix are very scattered, we aggregated data by mapping target domains to their respective superfamilies using SCOPe52 (accessed in Oct 2016). This reduced the 43 targets to 23 superfamilies. The 58 compounds were mapped to MeSH chemicals and drugs categories, as available from PubChem (accessed in Oct 2016), which reduced 58 compounds to 19 classes.

Data availability

All primary data used for the analysis stems from public sources, which are described in the methods section of the paper. If required, identifiers are provided in the corresponding sections. PDB identifiers for the screening set as well as detailed information on hit compounds and their PubChem Compound Identifiers (CIDs) are provided as Supplementary Information. Source data for Fig. 5 is provided as Supplementary Information. All other relevant data are available from the corresponding author on request.

Code availability

The Protein-Ligand Interaction Profiler (PLIP) as the key software for the study is publicly available at github.com/ssalentin/plip under the open source Apache 2 license.

Electronic supplementary material

Acknowledgements

We thank Alex Mestiashvili for technical support and Petra Lennig for help with the chemical classification of hit compounds and fruitful discussions. Furthermore, thanks go to Gary Jennings for proofreading and extensive feedback on the manuscript.

Author Contributions

Author Contributions

M.S., S.S., M.A. designed the study, analysed the data, and wrote the manuscript. S.S., M.A. carried out the study. J.H., S.D., J.H., Y.Z. contributed data.

Notes

Competing Interests

The authors declare that they have no competing interests.

Footnotes

Sebastian Salentin and Melissa F. Adasme contributed equally to this work.

Electronic supplementary material

Supplementary information accompanies this paper at 10.1038/s41598-017-11924-4.

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

1. DiMasi JA, Grabowski HG, Hansen RW. Innovation in the pharmaceutical industry: new estimates of R&D costs. J Health Econ. 2016;47:20–33. doi: 10.1016/j.jhealeco.2016.01.012. [PubMed] [Cross Ref]
2. Kim JH, Scialli AR. Thalidomide: The tragedy of birth defects and the effective treatment of disease. Toxicol Sci. 2011;122(1):1–6. doi: 10.1093/toxsci/kfr088. [PubMed] [Cross Ref]
3. Heinrich JC, Tuukkanen A, Schroeder M, Fahrig T, Fahrig R. RP101 (BVDU) binds to heat shock protein HSP27 (HSPB1) and enhances survival in animals and pancreatic cancer patients. J Cancer Res Clin. 2011;137(9):1349–61. doi: 10.1007/s00432-011-1005-1. [PubMed] [Cross Ref]
4. Haupt VJ, et al. Computational drug repositioning by target hopping: a use case in chagas disease drug repositioning by target hopping. Curr Pharm Des. 2016;22(21):3124–34. doi: 10.2174/1381612822666160224143008. [PubMed] [Cross Ref]
5. Berman HM, et al. The Protein Data Bank. Nucleic Acids Res. 2000;28(1):235–42. doi: 10.1093/nar/28.1.235. [PMC free article] [PubMed] [Cross Ref]
6. Deng Z, Chuaqui C, Singh J. Structural interaction fingerprint (SIFt): a novel method for analyzing three-dimensional protein-ligand binding interactions. J Med Chem. 2004;47(2):337–44. doi: 10.1021/jm030331x. [PubMed] [Cross Ref]
7. Mpamhanga CP, Chen B, McLay IM, Willett P. Knowledge-based interaction fingerprint scoring: a simple method for improving the effectiveness of fast scoring functions. J Chem Inf Model. 2006;46(2):686–98. doi: 10.1021/ci050420d. [PubMed] [Cross Ref]
8. Sato T, Honma T, Yokoyama S. Combining machine learning and pharmacophore-based interaction fingerprint for in silico screening. J Chem Inf Model. 2010;50(1):170–85. doi: 10.1021/ci900382e. [PubMed] [Cross Ref]
9. Wolber G, Langer T. LigandScout: 3-D pharmacophores derived from protein-bound ligands and their use as virtual screening filters. J Chem Inf Model. 2005;45(1):160–9. doi: 10.1021/ci049885e. [PubMed] [Cross Ref]
10. Tan L, Lounkine E, Bajorath J. Similarity searching using fingerprints of molecular fragments involved in protein-ligand interactions. J Chem Inf Model. 2008;48(12):2308–12. doi: 10.1021/ci800322y. [PubMed] [Cross Ref]
11. Batista J, Tan L, Bajorath J. Atom-centered interacting fragments and similarity search Applications. J Chem Inf Model. 2010;50(1):79–86. doi: 10.1021/ci9004223. [PubMed] [Cross Ref]
12. Salentin S, Schreiber S, Haupt VJ, Adasme MF, Schroeder M. PLIP: fully automated protein – ligand interaction profiler. Nucleic Acids Res. 2015;43(W1):W443–W447. doi: 10.1093/nar/gkv315. [PMC free article] [PubMed] [Cross Ref]
13. DeClercq E. Discovery and development of BVDU (brivudin) as a therapeutic for the treatment of herpes zoster. Biochem Pharmacol. 2004;68:2301–2315. doi: 10.1016/j.bcp.2004.07.039. [PubMed] [Cross Ref]
14. Heinrich JC, Donakonda S, Haupt VJ, Lennig P. New HSP27 inhibitors efficiently down-regulate resistance development in cancer cells. Oncotarget. 2016;7(42):68156–68169. doi: 10.18632/oncotarget.11905. [PMC free article] [PubMed] [Cross Ref]
15. Straume O, et al. Suppression of heat shock protein 27 induces long-term dormancy in human breast cancer. PNAS. 2012;109(22):8699–8704. doi: 10.1073/pnas.1017909109. [PubMed] [Cross Ref]
16. Bruey JM, et al. Hsp27 negatively regulates cell death by interacting with cytochrome c. Nat Cell Biol. 2000;2(9):645–652. doi: 10.1038/35023595. [PubMed] [Cross Ref]
17. Crisman TJ, Sisay MT, Bajorath J. Ligand-target interaction-based weighting of substructures for virtual screening. J Chem Inf Model. 2008;48(10):1955–1964. doi: 10.1021/ci800229q. [PubMed] [Cross Ref]
18. Tan L, Bajorath J. Utilizing target–ligand interaction information in fingerprint searching for ligands of related targets. Chem Biol Drug Des. 2009;74(1):25–32. doi: 10.1111/j.1747-0285.2009.00829.x. [PubMed] [Cross Ref]
19. Kelly MD, Mancera RL. Expanded interaction fingerprint method for analyzing ligand binding modes in docking and structure-based drug design. J Chemn Inf Comp Sci. 2004;44(6):1942–51. doi: 10.1021/ci049870g. [PubMed] [Cross Ref]
20. Chuaqui C, Deng Z, Singh J. Interaction profiles of protein kinase-inhibitor complexes and their application to virtual screening. J Med Chem. 2005;48(1):121–33. doi: 10.1021/jm049312t. [PubMed] [Cross Ref]
21. Deng Z, Chuaqui C, Singh J. Knowledge-based design of target-focused libraries using protein-ligand interaction constraints. J Med Chem. 2006;49(2):490–500. doi: 10.1021/jm050381x. [PubMed] [Cross Ref]
22. Clark A, Labute P. 2D depiction of protein-ligand complexes. J Chem Inf Model. 2007;47:1933–1944. doi: 10.1021/ci7001473. [PubMed] [Cross Ref]
23. Marcou G, Rognan D. Optimizing fragment and scaffold docking by use of molecular interaction fingerprints. J Chem Inf Model. 2007;47(1):195–207. doi: 10.1021/ci600342e. [PubMed] [Cross Ref]
24. Pérez-Nueno VI, Rabal O, Borrell JI, Teixidó J. APIF: a new interaction fingerprint based on atom pairs and its application to virtual screening. J Chem Inf Model. 2009;49(5):1245–60. doi: 10.1021/ci900043r. [PubMed] [Cross Ref]
25. Desaphy J, Raimbaud E, Ducrot P, Rognan D. Encoding protein-ligand interaction patterns in fingerprints and graphs. J Chem Inf Model. 2013;53(3):623–37. doi: 10.1021/ci300566n. [PubMed] [Cross Ref]
26. Qiao S, et al. The antimalarial amodiaquine causes autophagic-lysosomal and proliferative blockade sensitizing human melanoma cells to starvation- and chemotherapy-induced cell death. Autophagy. 2013;9(12):2087–2102. doi: 10.4161/auto.26506. [PMC free article] [PubMed] [Cross Ref]
27. Kondo Y, Kanzawa T, Sawaya R, Kondo S. The role of autophagy in cancer development and response to therapy. Nat Rev Cancer. 2005;5(9):726–734. doi: 10.1038/nrc1692. [PubMed] [Cross Ref]
28. Rubinsztein DC, Codogno P, Levine B. Autophagy modulation as a potential therapeutic target for diverse diseases. Nat Rev Drug Discov. 2012;11(9):709–730. doi: 10.1038/nrd3802. [PMC free article] [PubMed] [Cross Ref]
29. Boonyasuppayakorn S, Reichert DE, Manzano M, Nagarajan K, Padmanabhan R. Amodiaquine, an antimalarial drug, inhibits dengue virus type 2 replication and infectivity. Antiviral Res. 2014;106:125–134. doi: 10.1016/j.antiviral.2014.03.014. [PMC free article] [PubMed] [Cross Ref]
30. Fischer AM, et al. Flaviviruses are sensitive to inhibition of thymidine synthesis pathways. J Virol. 2013;87(17):9411–9419. doi: 10.1128/JVI.00101-13. [PMC free article] [PubMed] [Cross Ref]
31. Zhang T, et al. Chrysin and its phosphate ester inhibit cell proliferation and induce apoptosis in hela cells. Bioorg Med Chem. 2004;12(23):6097–105. doi: 10.1016/j.bmc.2004.09.013. [PubMed] [Cross Ref]
32. Lirdprapamongkol K, et al. A flavonoid chrysin suppresses hypoxic survival and metastatic growth of mouse breast cancer cells. Oncol Rep. 2013;30(5):2357–2364. doi: 10.3892/or.2013.2667. [PubMed] [Cross Ref]
33. Kasala ER, et al. Chemopreventive and therapeutic potential of chrysin in cancer: mechanistic perspectives. Toxicol Lett. 2015;233(2):214–225. doi: 10.1016/j.toxlet.2015.01.008. [PubMed] [Cross Ref]
34. Ronnekleiv-Kelly S, et al. Aryl hydrocarbon receptor-dependent apoptotic cell death induced by the flavonoid chrysin in human colorectal cancer cells. Cancer Lett. 2016;370(1):91–9. doi: 10.1016/j.canlet.2015.10.014. [PMC free article] [PubMed] [Cross Ref]
35. Jabareen, A., Huleihil, M., Huleihel, M. Effect of extracts of passiflora edulis leaves on herpes viruses infection. J Virol Antivir Res. 2(2) (2013).
36. Lyu SY, Rhim JY, Park WB. Antiherpetic activities of flavonoids against herpes simplex virus type 1 (HSV-1) and type 2 (HSV-2) in vitro. Arch Pharm Res. 2005;28(11):1293–1301. doi: 10.1007/BF02978215. [PubMed] [Cross Ref]
37. Schnitzler P, et al. Antiviral Activity and Mode of Action of Propolis Extracts and Selected Compounds. Phytother Res. 2010;24(S1):S20–8. doi: 10.1002/ptr.2868. [PubMed] [Cross Ref]
38. Kurokawa M, et al. Efficacy of Brazilian propolis against herpes simplex virus type 1 infection in mice and their modes of antiherpetic efficacies. Evid Based Complement Alternat Med. 2011;2011:976196. [PMC free article] [PubMed]
39. Wang J, et al. Anti-enterovirus 71 effects of chrysin and its phosphate ester. PLoS One. 2014;9(3):e89668. doi: 10.1371/journal.pone.0089668. [PMC free article] [PubMed] [Cross Ref]
40. Baell JB. Feeling nature’s PAINS: natural products, natural product drugs, and pan assay interference compounds (PAINS) J Nat Prod. 2016;79(3):616–628. doi: 10.1021/acs.jnatprod.5b00947. [PubMed] [Cross Ref]
41. Trivella DBB, dos Reis CV, Lima LMTR, Foguel D, Polikarpov I. Flavonoid interactions with human transthyretin: combined structural and thermodynamic analysis. J Struct Biol. 2012;180(1):143–153. doi: 10.1016/j.jsb.2012.07.008. [PubMed] [Cross Ref]
42. Fu B, et al. Chrysin inhibits expression of hypoxia-inducible factor-1alpha through reducing hypoxia-inducible factor-1alpha stability and inhibiting its protein synthesis. Mol Cancer Ther. 2007;6(1):220–226. doi: 10.1158/1535-7163.MCT-06-0526. [PubMed] [Cross Ref]
43. Samarghandian S, Afshari JT, Davoodi S. Chrysin reduces proliferation and induces apoptosis in the human prostate cancer cell line pc-3. Clinics (Sao Paulo). 2011;66(6):1073–1079. doi: 10.1590/S1807-59322011000600026. [PMC free article] [PubMed] [Cross Ref]
44. Sun X, et al. The anticancer flavonoid chrysin induces the unfolded protein response in hepatoma cells. J Cell Mol Med. 2011;15(11):2389–2398. doi: 10.1111/j.1582-4934.2010.01244.x. [PMC free article] [PubMed] [Cross Ref]
45. Excelra Knowledge Services. Drug Repurposing Portal. drugrepurposingportal.com (2017).
46. Santos R, et al. A comprehensive map of molecular drug targets. Nat Rev Drug Discov. 2017;16:19–34. doi: 10.1038/nrd.2016.230. [PubMed] [Cross Ref]
47. Rose PW, et al. The RCSB Protein Data Bank: Views of structural biology for basic and applied research and education. Nucleic Acids Res. 2015;43(D1):D345–D356. doi: 10.1093/nar/gku1214. [PMC free article] [PubMed] [Cross Ref]
48. Yang J, Roy A, Zhang Y. BioLiP: a semi-manually curated database for biologically relevant ligand-protein interactions. Nucleic Acids Res. 2013;41(D1):D1096–D1103. doi: 10.1093/nar/gks966. [PMC free article] [PubMed] [Cross Ref]
49. Chen X, Ji ZL, Chen YZ. TTD: Therapeutic target database. Nucleic Acids Res. 2002;30(1):412–415. doi: 10.1093/nar/30.1.412. [PMC free article] [PubMed] [Cross Ref]
50. Kuhn M, et al. STITCH 4: Integration of protein-chemical interactions with user data. Nucleic Acids Res. 2014;42(D1):401–407. doi: 10.1093/nar/gkt1207. [PMC free article] [PubMed] [Cross Ref]
51. Bento AP, et al. The ChEMBL bioactivity database: an update. Nucleic Acids Res. 2014;42(D1):1083–1090. doi: 10.1093/nar/gkt1031. [PMC free article] [PubMed] [Cross Ref]
52. Lo Conte L, Brenner SE, Hubbard TJP, Chothia C, Murzin AG. SCOP database in 2002: refinements accommodate structural genomics. Nucleic Acids Res. 2002;30(1):264–267. doi: 10.1093/nar/30.1.264. [PMC free article] [PubMed] [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group