A common feature pharmacophore for FDA-approved drugs inhibiting the Ebola virus

Sean Ekins; Joel S. Freundlich; Megan Coffee

doi:10.12688/f1000research.5741.1

Home Browse A common feature pharmacophore for FDA-approved drugs inhibiting the...

ALL Metrics

Views

Downloads

Get PDF

Get XML

Export

▬

✚

Research Note

A common feature pharmacophore for FDA-approved drugs inhibiting the Ebola virus

[version 1; peer review: 2 approved]

Sean Ekins ^1,2, Joel S. Freundlich³, Megan Coffee⁴

PUBLISHED 14 Nov 2014

Author details Author details

¹ Collaborative Drug Discovery, Burlingame, CA, 94010, USA
² Collaborations in Chemistry, Fuquay-Varina, NC, 27526, USA
³ Departments of Pharmacology & Physiology and Medicine, Center for Emerging and Reemerging Pathogens, UMDNJ - New Jersey Medical School, NJ, 07103, USA
⁴ Center for Infectious Diseases and Emerging Readiness, University of California, Berkeley, CA, 94720, USA

OPEN PEER REVIEW

REVIEWER STATUS

This article is included in the Emerging Diseases and Outbreaks gateway.

This article is included in the Ebola Virus collection.

Abstract

We are currently faced with a global infectious disease crisis which has been anticipated for decades. While many promising biotherapeutics are being tested, the search for a small molecule has yet to deliver an approved drug or therapeutic for the Ebola or similar filoviruses that cause haemorrhagic fever. Two recent high throughput screens published in 2013 did however identify several hits that progressed to animal studies that are FDA approved drugs used for other indications. The current computational analysis uses these molecules from two different structural classes to construct a common features pharmacophore. This ligand-based pharmacophore implicates a possible common target or mechanism that could be further explored. A recent structure based design project yielded nine co-crystal structures of pyrrolidinone inhibitors bound to the viral protein 35 (VP35). When receptor-ligand pharmacophores based on the analogs of these molecules and the protein structures were constructed, the molecular features partially overlapped with the common features of solely ligand-based pharmacophore models based on FDA approved drugs. These previously identified FDA approved drugs with activity against Ebola were therefore docked into this protein. The antimalarials chloroquine and amodiaquine docked favorably in VP35. We propose that these drugs identified to date as inhibitors of the Ebola virus may be targeting VP35. These computational models may provide preliminary insights into the molecular features that are responsible for their activity against Ebola virus in vitro and in vivo and we propose that this hypothesis could be readily tested.

Keywords

ebola virus, computational models, machine learning

Corresponding author: Sean Ekins

Competing interests: S.E. works for Collaborations in Chemistry, and consults for Collaborative Drug Discovery Inc.

Grant information: The author(s) declared that no grants were involved in supporting this work.

Copyright: © 2014 Ekins S et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Data associated with the article are available under the terms of the Creative Commons Zero "No rights reserved" data waiver (CC0 1.0 Public domain dedication).

How to cite: Ekins S, Freundlich JS and Coffee M. A common feature pharmacophore for FDA-approved drugs inhibiting the Ebola virus [version 1; peer review: 2 approved]. F1000Research 2014, 3:277 (https://doi.org/10.12688/f1000research.5741.1) First published: 14 Nov 2014, 3:277 (https://doi.org/10.12688/f1000research.5741.1) Latest published: 12 Dec 2014, 3:277 (https://doi.org/10.12688/f1000research.5741.2)

Introduction

The current Ebola virus (EBOV) crisis has demonstrated that globally we are not prepared to respond with therapeutics to treat existing infections or act as prophylactics as there is no Food and Drug Administration (FDA) or European Medicines Agency (EMEA) approved therapeutic. More importantly this suggests we should have been prepared for a pathogen which has been known about for nearly forty years. The current EBOV outbreak is already proving remarkably costly in terms of the mortality and financial ramifications^1,2. The best approaches to EBOV so far have relied on public health measures for containment³ which have been used in past outbreaks⁴. These lessons with EBOV will undoubtedly be important for the next virus outbreak⁵ but they also raise many questions⁶ which point to how little we know about these viruses in general, as well as how best to share knowledge openly⁷.

There have been a relatively small number of studies that have attempted to identify compounds active against EBOV. Two recent studies utilized high-throughput screens of a subset of FDA approved drugs against different EBOV strains (Zaire and Sudan) in vitro and in vivo. These independent reports suggested the promise of the antimalarials amodiaquine and chloroquine in one study⁸, while the selective estrogen receptor modulators (SERMs) clomiphene and toremifene were active in another⁹. Chloroquine to date has not progressed beyond the mouse EBOV model used in these studies. We hypothesized that we could use these four molecules to computationally define the features that are important for activity. The previous studies were not exhaustive screens of all FDA drugs and so we have taken this opportunity to suggest additional compounds. Looked at from another perspective “non-antiviral” drugs may be worth following up even though their molecular mechanism is unknown. These compounds may themselves have broad antiviral activity as reports describe modest inhibitory activity against other viruses^10–13.

Several studies have identified non-FDA approved drugs including an in silico docking approach to identify molecules targeting the viral Nedd4-PPxY interface¹⁴. These molecules were similar to the FDA benzimidazole and aminoquinoline^8,9 compounds that were active against EBOV. Another good example is the recent in silico docking of 5.4 million drug-like compounds docked in the viral protein VP35 protein¹⁵. This identified multiple pyrrolidinones which inhibit its polymerase cofactor activity¹⁵. The pyrrolidinones bind to an alpha helix which is proposed as important for viral function¹⁶. With the limited knowledge of small molecules and potential targets we have studied whether the FDA-approved drugs that are active in vitro and in vivo versus EBOV could be targeting VP35.

Methods

Common features pharmacophore for EBOV actives

Two papers from 2013 described compounds active as inhibitors of different EBOV strains in vitro and in vivo, namely amodiaquine and chloroquine in one study⁸, clomiphene and toremifene in another⁹. These active molecules were used as they have both in vitro and in vivo activity to build a common features pharmacophore with Discovery Studio 4.1 (Biovia, San Diego, CA) from 3D conformations of the molecules generated with the CAESAR algorithm. This identified key features. The pharmacophore was then used to search various databases (for which up to 100 molecule conformations with the FAST conformer generation method with the maximum energy threshold of 20 kcal/mol, were created). The pharmacophore was then used to search the Microsource Spectrum database (http://www.msdiscovery.com/spectrum.html) as well as the CDD FDA drugs dataset (https://www.collaborativedrug.com/pages/public_access). In both cases over 300 hits were retrieved initially. The van der Waals surface of amodiaquine (which was more potent than chloroquine⁸) was added to limit the number of hits retrieved^17–19.

Receptor-ligand pharmacophores for VP35

Receptor-ligand pharmacophores for the VP35 protein were generated from crystal structures (4IBB, 4IBC, 4IBD, 4IBE, 4IBF, 4IBG, 4IBI, 4IBJ, 4IBK) in the protein data bank PDB. Pharmacophores were constructed using the receptor-ligand pharmacophore generation protocol in Discovery Studio version 4.1 (Biovia, San Diego, CA) with a maximum number of pharmacophores (10), minimum features (4), and maximum number of features (6) as are described elsewhere²⁰.

in silico docking of molecules in VP35 structure

PDB 4IBI was used for docking using LibDock in Discovery Studio (Biovia, San Diego CA)²¹. The proposed binding site was centered on the ligand and a site sphere created (coordinates 2.14, 20.93, 1.71) with 9.45 Å diameter. The protocol included 10 hotspots and docking tolerance (0.25). The FAST conformation method was also used along with steepest descent minimization with CHARMm. Further parameters followed the default settings. The ligand VPL57 was removed from the binding site and re-docked. The four FDA approved drugs with activity against Ebola were docked in the structure from an sdf file. Molecules were visualized alongside the original ligand VPL57 and the 2D interaction plots generated.

Results

Dataset 1.Pharmacophores, receptor ligand models and docking data for FDA-approved drugs inhibiting the Ebola virus.

Data was downloaded sourced from Microsource Spectrum and CDD Drugs. Dataset includes sd files used to create the 3D database that was searched. Note that models only run on Discovery Studio.

Common features pharmacophore for EBOV actives

The pharmacophore was generated using the in vivo and in vitro active amodiaquine, chloroquine, clomiphene and toremifene (Supplemental Table 1) as these represent the most relevant FDA approved drugs to date. This pharmacophore consists of 4 hydrophobic features and a hydrogen bond acceptor feature (Figure 1). The pharmacophore with van der Waals surface was also used to search FDA drug various libraries (Supplemental Table 2 and Supplemental Table 3). The most interesting observations from this virtual screen are that various estradiol analogs score well (e.g. estradiol valerate Fit value 4.23). Previously estradiol was suggested to be active in the EBOV pseudotype assay in vitro⁸. In addition, dibucaine was also retrieved (Fit value 1.58) which was also active in the EBOV pseudotype assay⁸. Amodiaquine, chloroquine, clomiphene and toremifene can be used as positive controls for future screens. Because the original complete sets of FDA approved compounds screened are not publically accessible it is difficult to compare hit rates versus all compounds tested to date.

Figure 1. Pharmacophore based on 4 hits.

A. amodiaquine, B. chloroquine, C. clomiphene D. toremifene and E. Overlap showing all molecules in the van der Waals surface of amodiaquine.

Receptor-ligand pharmacophores for VP35

The nine receptor-ligand pharmacophores created all consisted of three to four hydrophobic features and one to two hydrogen bonding features (Table 1). Eight of these pharmacophores also had a negative ionizable feature. These suggest that the receptor-ligand based approach results in a general similarity across the nine structures, likely indicating the similar binding mode and importance of features for interfering with this generally hydrophobic pocket for protein-protein interactions.

Table 1. Pharmacophores for EBOV VP35 generated from crystal structures in the protein data bank PDB.

Pharmacophores were generated using the receptor-ligand pharmacophore generation protocol in Discovery Studio version 4.1 (Biovia, San Diego, CA) with minimum features (3) and maximum features (6). Pharmacophore features are Hydrophobic (H, cyan), Hydrogen bond acceptor (HBA, green), hydrogen bond donor (HBD, purple) and 1 negative ionizable (neg, blue). Excluded volumes (grey) were also automatically added. Further details on this approach are described elsewhere²⁰.

PDB	Pharmacophore features	Pharmacophore with ligand mapped
4IBB	4H, 1HBD, 1 neg ionizable
4IBC	3H, 2HBA, 1 neg ionizable
4IBD	4H, 1 HBA, 1 neg ionizable
4IBE	4H, 1HBA
4IBF	4H, 1 HBA, 1 neg ionizable
4IBG	3H, 2 HBA, 1 neg ionizable
4IBI	4H, 1HBA, 1 neg ionizable
4IBJ	4H, 1HBA, 1 neg ionizable
4IBK	4H, 1HBA, 1 neg ionizable

in silico docking of molecules in VP35 structure

Redocking the 4IBI ligand in the protein resulted in an RMSD of 3.02Å, which generally indicates the difficulty of predicting orientations for compounds binding in what is a relatively hydrophobic and shallow pocket (Figure S1). This molecule was ranked the 29^th pose and had a LibDock score of 86.62 (Figure S1). The four FDA approved drugs were docked into the VP35 structure 4IBI. All compounds docked similarly and overlapped with the co-crystal ligand (Figure 2). Amodiaquine and chloroquine had higher LibDock scores (> 90) than the 4IBI ligand, while clomiphene and toremifene had LibDock scores less than 70. All four FDA approved drugs bound similarly to the pyrrolidinone ligands in the pocket formed by residues from the α-helical and β-sheet subdomains¹⁵. We have highlighted proposed energetically favorable interactions of the antimalarial candidate binders with ILE295, LYS248 and GLN244, which scored favorably. Previously published studies suggested mutation of ILE295, LYS248 resulted in near-complete loss of binding activity¹⁵.

Figure 2. Docking FDA approved compounds in VP35 protein showing overlap with ligand (yellow) and 2D interaction diagram.

4IBI was used, 4IBI ligand VPL57 shown in yellow. A. Amodiaquine (grey) and 4IBI LibDock score 90.80, B. Chloroquine (grey) LibDock score 97.82, C. Clomiphene (grey) and 4IBI LibDock score 69.77, D. Toremifene (grey) and 4IBI LibDock score 68.11

Discussion

Our previous experience with common feature and quantitative pharmacophore models has demonstrated their value in predicting novel actives from collections of FDA approved drugs^22–27. Candidate predicted actives may be assessed by their Fit Value to the pharmacophore model. This score can be used to prioritize compounds for eventual testing. In the current study it was hypothesized that two different classes of compounds showing activity against EBOV in vitro and in vivo may share a common pharmacophore. Construction of this pharmacophore (Figure 1) indicated four hydrophobic features and a hydrogen bond acceptor feature. This pharmacophore (with an added van der Waals surface to limit the number of hits retrieved) was then used to screen and score other FDA drugs from a small database and identified 120 and 124 structures for future evaluation in vitro testing (Supplemental Table 2 and Supplemental Table 3). Out of these compounds estradiol and dibucaine had been previously described as active in in vitro EBOV assays. This suggested the pharmacophore could retrieve some structurally diverse classes of known hits⁸.

Recently identified co-crystal structures of the EBOV VP35 protein were used to derive receptor-ligand pharmacophores. These nine receptor-ligand pharmacophores suggested the importance of hydrophobic, hydrogen bonding and negative ionizable interactions to interfere with this protein-protein interaction (Table 1). Eight out of nine of the pharmacophores had one or more hydrogen bond acceptor feature. These pharmacophores are grossly similar to our ligand based pharmacophore (derived from four FDA approved drugs that inhibit EBOV), as both types of model had multiple hydrophobic features and at least one hydrogen bond acceptor. When we docked the antimalarials and SERMs into a representative VP35 structure these compounds were found to overlap with the X-ray ligand to differing extents. Amodiaquine and chloroquine had LibDock scores greater than 90 and higher than that for the redocked X-ray ligand. This indicated that VP35 may be a potential target for these two distinct classes of compounds. However, it is important to point out that we have not compared docking to other proteins in EBOV and it could also be possible that these molecules are active elsewhere as well as via other mechanisms than by specific binding to proteins^28,29. Further, VP35 may be a preferred target for the antimalarials while the SERMs are not predicted to bind as well as the X-ray ligand. The use of other docking and scoring methods may produce differences in the pose and predicted binding affinity, which could be of interest for further studies.

A combination of the promising efficacy of chloroquine (EC₅₀ 16 μM⁸) and amodiaquine (EC₅₀ 8.4 μM⁸) versus EBOV, their availability and likely low cost should prioritize their further laboratory exploration. Mechanistic studies against VP35 and possibly other proteins should also be pursued and may be enlightened by the observation that both of these compounds also have reported activity against other viruses. For example, chloroquine is active against human coronavirus OC43 (in vitro and in infected mice) as well as SARS (in vitro)^10,30,31, while amodiaquine also inhibits dengue virus 2 replication and infectivity in vitro¹¹.

Conclusions

In summary, this study has built on the previous publications that identified four FDA approved compounds active against different strains of EBOV^8,9. Our pharmacophore model for SERMs and aminoquinolines suggests that these compounds share multiple chemical features based on their overlap to the four hydrophobic features and a hydrogen bond acceptor (Figure 1E) and they may have a common mechanism or target. We suggest that VP35 may be the likely target based on the overlap of receptor-based pharmacophores and docking into the crystal structure. Amodiaquine and chloroquine score particularly well in terms of docking to VP35. If this is the case it could provide a means to follow up with other small molecule analogs and/or additional FDA approved drugs that could target this protein-protein interaction. As with our other tuberculosis-focused research^32,33, and computational approaches to repositioning compounds³⁴ we embrace the essentiality for computational predictions to be interrogated through rigorous experimental studies. For example at least two in silico docking studies screened commercially available compounds^14,15. We propose that docking FDA approved drugs could also be a viable first step to identifying potential compounds that could be used. We are actively seeking collaborators with experience with EBOV assays to enable further translational studies. We believe this computationally inspired approach may be applicable for other known infectious pathogens that do not have current treatments such as other viruses related to Ebola. Ultimately we need to be able to leverage such approaches to provide antivirals for future pathogens.

Data availability

F1000Research: Dataset 1. Pharmacophores, receptor ligand models and docking data for FDA-approved drugs inhibiting the Ebola virus, 10.5256/f1000research.5741.d38449³⁵.

The ligand-based pharmacophore was previously made available: http://figshare.com/articles/Ebola_active_cpds_pharmacophore/1190902.

The following PDB structures were used in this study (4IBB, 4IBC, 4IBD, 4IBE, 4IBF, 4IBG, 4IBI, 4IBJ, 4IBK).

For models and advice please contact Sean Ekins (ekinssean@yahoo.com).

Author contributions

S.E. and M.C. came up with the general idea for the study based on the published in vitro and in vivo data. All authors contributed to the collaborative writing of this project.

Competing interests

S.E. works for Collaborations in Chemistry, and consults for Collaborative Drug Discovery Inc.

Grant information

The author(s) declared that no grants were involved in supporting this work.

Acknowledgments

Dr. Christopher D. Southan, Dr. Peter Madrid and Dr. Nadia Litterman are acknowledged for discussions on Ebola. Biovia are kindly acknowledged for providing Discovery Studio. An earlier preliminary version of this pharmacophore was described previously: http://figshare.com/articles/A_pharmacophore_of_ebola_active_compounds/1190787.

Supplementary materials

Supplemental Table 1. Information for common features pharmacophore generation.

	Maximum Omittted features	Principal
chloroquine	0	1
amodiaquine	0	2
clomiphene	0	1
Toremifene	0	2

Supplemental Table 2. FDA drugs and common features pharmacophore. The dataset of 2643 molecules was downloaded from the CDD Public Access (https://www.collaborativedrug.com/pages/public_access) as an sdf and then a 3D database was created in Discovery Studio using FAST conformer generation with up to 255 conformations. The database was searched with the common feature pharmacophore developed from amodiaquine, chloroquine, clomiphene and toremifene. The search 3D database protocol was used with the Fast search method. In some cases the indication for the molecules is not described (ND).

Supplemental Table 3. Microsource Spectrum and common features pharmacophore. The dataset of 2311 molecules was provided by Microsource (http://www.msdiscovery.com/spectrum.html) as an sdf and then a 3D database was created in Discovery Studio using FAST conformer generation with up to 255 conformations. The database was searched with the common feature pharmacophore developed from amodiaquine, chloroquine, clomiphene and toremifene. The search 3D database protocol was used with the Fast search method.

Figure S1.Redocking VPL57 in 4IBI.

The 4IBI ligand was removed from the structure and redocked. The closest pose (grey) was ranked 29 with RMSD 3.02A and LibDock score 86.62 when compared to the actual ligand in 4IBI (yellow).

Faculty Opinions recommended

References

1. Butler D, Morello L: Ebola by the numbers: The size, spread and cost of an outbreak. Nature. 2014; 514(7522): 284–5. PubMed Abstract | Publisher Full Text
2. Piot P: Ebola’s perfect storm. Science. 2014; 345(6202): 1221. PubMed Abstract | Publisher Full Text
3. Trad MA, Fisher DA, Tambyah PA: Ebola in west Africa. Lancet Infect Dis. 2014; 14(11): 1045. PubMed Abstract | Publisher Full Text
4. Dhillon RS, Srikrishna D, Sachs J: Controlling Ebola: next steps. Lancet. 2014; 384(9952): 1409–1411. PubMed Abstract | Publisher Full Text
5. Anon: Call to action. Nature. 2014; 514(7524): 535–536. PubMed Abstract | Publisher Full Text
6. Check Hayden E: The Ebola questions. Nature. 2014; 514(7524): 554–7. PubMed Abstract | Publisher Full Text
7. Piot P: The F1000Research: Ebola article collection [v1; ref status: not peer reviewed, http://f1000r.es/4ot]. F1000Res. 2014; 3: 269. Publisher Full Text
8. Madrid PB, Chopra S, Manger ID, et al.: A systematic screen of FDA-approved drugs for inhibitors of biological threat agents. PLoS One. 2013; 8(4): e60579. PubMed Abstract | Publisher Full Text | Free Full Text
9. Johansen LM, Brannan JM, Delos SE, et al.: FDA-approved selective estrogen receptor modulators inhibit Ebola virus infection. Sci Transl Med. 2013; 5(190): 190ra79. PubMed Abstract | Publisher Full Text | Free Full Text
10. Keyaerts E, Li S, Vijgen L, et al.: Antiviral activity of chloroquine against human coronavirus OC43 infection in newborn mice. Antimicrob Agents Chemother. 2009; 53(8): 3416–21. PubMed Abstract | Publisher Full Text | Free Full Text
11. Boonyasuppayakorn S, Reichert ED, Manzano M, et al.: Amodiaquine, an antimalarial drug, inhibits dengue virus type 2 replication and infectivity. Antiviral Res. 2014; 106: 125–34. PubMed Abstract | Publisher Full Text
12. Acosta EG, Bruttomesso AC, Bisceglia JA, et al.: Dehydroepiandrosterone, epiandrosterone and synthetic derivatives inhibit Junin virus replication in vitro. Virus Res. 2008; 135(2): 203–12. PubMed Abstract | Publisher Full Text
13. Dyall J, Coleman CM, Hart BJ, et al.: Repurposing of clinically developed drugs for treatment of Middle East respiratory syndrome coronavirus infection. Antimicrob Agents Chemother. 2014; 58(8): 4885–93. PubMed Abstract | Publisher Full Text | Free Full Text
14. Han Z, Lu J, Liu Y, et al.: Small-molecule probes targeting the viral PPxY-host Nedd4 interface block egress of a broad range of RNA viruses. J Virol. 2014; 88(13): 7294–306. PubMed Abstract | Publisher Full Text | Free Full Text
15. Brown CS, Lee MS, Leung DW, et al.: In silico derived small molecules bind the filovirus VP35 protein and inhibit its polymerase cofactor activity. J Mol Biol. 2014; 426(10): 2045–58. PubMed Abstract | Publisher Full Text | Free Full Text
16. Chakraborty S, Rao B, Asgeirsson B, et al.: Characterizing alpha helical properties of Ebola viral proteins as potential targets for inhibition of alpha-helix mediated protein-protein interactions [v1; ref status: approved with reservations 1, http://f1000r.es/4lg]. F1000Res. 2014; 3: 251. Publisher Full Text
17. Lamichhane G, Freundlich JS, Ekins S, et al.: Essential metabolites of Mycobacterium tuberculosis and their mimics. MBio. 2011; 2(1): e00301–10. PubMed Abstract | Publisher Full Text | Free Full Text
18. Ekins S, Bradford J, Dole K, et al.: A collaborative database and computational models for tuberculosis drug discovery. Mol Biosyst. 2010; 6(5): 840–851. PubMed Abstract | Publisher Full Text
19. Zheng X, Ekins S, Raufman JP, et al.: Computational models for drug inhibition of the human apical sodium-dependent bile acid transporter. Mol Pharm. 2009; 6(5): 1591–1603. PubMed Abstract | Publisher Full Text | Free Full Text
20. Meslamani J, Li J, Sutter J, et al.: Protein-ligand-based pharmacophores: generation and utility assessment in computational ligand profiling. J Chem Inf Model. 2012; 52(4): 943–55. PubMed Abstract | Publisher Full Text
21. Rao SN, Head MS, Kulkarni A, et al.: Validation studies of the site-directed docking program LibDock. J Chem Inf Model. 2007; 47(6): 2159–71. PubMed Abstract | Publisher Full Text
22. Diao L, Ekins S, Polli JE: Novel inhibitors of human organic cation/carnitine transporter (hOCTN2) via computational modeling and in vitro testing. Pharm Res. 2009; 26(8): 1890–1900. PubMed Abstract | Publisher Full Text | Free Full Text
23. Zheng X, Ekins S, Raufman JP, et al.: Computational models for drug inhibition of the human apical sodium-dependent bile acid transporter. Mol Pharm. 2009; 6(5): 1591–603. PubMed Abstract | Publisher Full Text | Free Full Text
24. Diao L, Ekins S, Polli JE: Quantitative structure activity relationship for inhibition of human organic cation/carnitine transporter. Mol Pharm. 2010; 7(6): 2120–2131. PubMed Abstract | Publisher Full Text | Free Full Text
25. Ekins S, Diao L, Polli JE: A substrate pharmacophore for the human organic cation/carnitine transporter identifies compounds associated with rhabdomyolysis. Mol Pharm. 2012; 9(4): 905–913. PubMed Abstract | Publisher Full Text | Free Full Text
26. Dong Z, Ekins S, Polli JE: Structure-activity relationship for FDA approved drugs as inhibitors of the human sodium taurocholate cotransporting polypeptide (NTCP). Mol Pharm. 2013; 10(3): 1008–19. PubMed Abstract | Publisher Full Text | Free Full Text
27. Dong Z, Ekins S, Polli JE: Quantitative NTCP pharmacophore and lack of association between DILI and NTCP Inhibition. Eur J Pharm Sci. 2014; 66C: 1–9. PubMed Abstract | Publisher Full Text
28. Kazmi F, Hensley T, Pope C, et al.: Lysosomal sequestration (trapping) of lipophilic amine (cationic amphiphilic) drugs in immortalized human hepatocytes (Fa2N-4 cells). Drug Metab Dispos. 2013; 41(4): 897–905. PubMed Abstract | Publisher Full Text | Free Full Text
29. Nadanaciva S, Lu S, Gebhard DF, et al.: A high content screening assay for identifying lysosomotropic compounds. Toxicol In Vitro. 2011; 25(3): 715–23. PubMed Abstract | Publisher Full Text
30. de Wilde AH, Jochmans D, Posthuma CC, et al.: Screening of an FDA-approved compound library identifies four small-molecule inhibitors of Middle East respiratory syndrome coronavirus replication in cell culture. Antimicrob Agents Chemother. 2014; 58(8): 4875–84. PubMed Abstract | Publisher Full Text | Free Full Text
31. Vincent MJ, Bergeron E, Benjannet S, et al.: Chloroquine is a potent inhibitor of SARS coronavirus infection and spread. Virol J. 2005; 2: 69. PubMed Abstract | Publisher Full Text | Free Full Text
32. Ekins S, Freundlich JS: Computational models for tuberculosis drug discovery. Methods Mol Biol. 2013; 993: 245–62. PubMed Abstract | Publisher Full Text
33. Ekins S, Freundlich JS, Choi I, et al.: Computational databases, pathway and cheminformatics tools for tuberculosis drug discovery. Trends Microbiol. 2011; 19(2): 65–74. PubMed Abstract | Publisher Full Text | Free Full Text
34. Ekins S, Williams AJ, Krasowski MD, et al.: In silico repositioning of approved drugs for rare and neglected diseases. Drug Disc Today. 2011; 16(7–8): 298–310. PubMed Abstract | Publisher Full Text
35. Ekins S, Freundlich JS, Coffee M: Dataset 1 in: A common feature pharmacophore for FDA-approved drugs inhibiting the Ebola virus. F1000Research. 2014. Data Source

Comments on this article Comments (0)

Version 2

VERSION 2 PUBLISHED 14 Nov 2014

Author details Author details

Competing interests

S.E. works for Collaborations in Chemistry, and consults for Collaborative Drug Discovery Inc.

Grant information

The author(s) declared that no grants were involved in supporting this work.

Article Versions (2)

version 2

Revised

Published: 12 Dec 2014, 3:277

https://doi.org/10.12688/f1000research.5741.2

version 1

Published: 14 Nov 2014, 3:277

https://doi.org/10.12688/f1000research.5741.1

© 2014 Ekins S et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Data associated with the article are available under the terms of the Creative Commons Zero "No rights reserved" data waiver (CC0 1.0 Public domain dedication).

Download

Export To

metrics

	Views	Downloads
F1000Research	-	-
PubMed Central Data from PMC are received and updated monthly.	-	-

Citations

SEE MORE DETAILS

CITE

how to cite this article

Ekins S, Freundlich JS and Coffee M. A common feature pharmacophore for FDA-approved drugs inhibiting the Ebola virus [version 1; peer review: 2 approved] F1000Research 2014, 3:277 (https://doi.org/10.12688/f1000research.5741.1)

NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.

track

receive updates on this article

Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status: ?

Key to Reviewer Statuses VIEW HIDE

ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested

Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.

Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions

Version 1

VERSION 1

PUBLISHED 14 Nov 2014

Views

Reviewer Report 02 Dec 2014

Sandeep Chakraborty, Plant Sciences Department, University of California, Davis, CA, USA

Approved

https://doi.org/10.5256/f1000research.6137.r6708

Ekins et al. have presented a crisp and lucid manuscript on a very relevant topic. They have suggested a methodology to extract common features from four approved compounds that have recently been found to work against the Ebola virus (amodiaquine, chloroquine, clomiphene and toremifene), and define a pharma-
cophore, which has been used to search databases, and identify further compounds for in vitro and in vivo testing. The in silico methodology described here provides an excellent method to quickly screen known compounds for possible therapies against Ebola in particular, and other viruses in general.

Some minor comments.

The result that SERMs show lower scores for binding to VP35 is rationalized by the finding that ‘clomiphene and toremifene inhibit EBOV VLP entry with some specificity to GP’¹, and therefore does not probably inhibit VP35.
‘Chloroquine to date has not progressed beyond the mouse EBOV model used in these studies.’ This statement is not clear, does it mean that the others have progressed beyond the mouse EBOV model?
The first few compounds in Supplemental Table 2 should be part of the main manuscript as a table.
Color coding of pharmacophore features should be in Fig 1 too (it comes earlier than Table 1, where it is described).
The structures look better with a white background.
A 3 Å RMSD for redocking a given ligand is quite high². The authors should consider the use of other docking methods, as a comparison.
A table of Libdock scores would help easily analyze results (with a mention of whether a higher score is better, and the significance of a score).

References

1. Johansen LM, Brannan JM, Delos SE, Shoemaker CJ, et al.: FDA-approved selective estrogen receptor modulators inhibit Ebola virus infection.Sci Transl Med.2013; 5 (190): 190ra79 PubMed Abstract | Free Full Text | Publisher Full Text
2. Sándor M, Kiss R, Keseru GM: Virtual fragment docking by Glide: a validation study on 190 protein-fragment complexes.J Chem Inf Model.2010; 50 (6): 1165-1172 PubMed Abstract | Publisher Full Text | Reference Source

Competing Interests: No competing interests were disclosed.

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

CITE

Report a concern

Author Response 12 Dec 2014

Sean Ekins, Collaborative Drug Discovery, Burlingame, 94010, USA

12 Dec 2014

Author Response
Ekins et al. have presented a crisp and lucid manuscript on a very relevant topic. They have suggested a methodology to extract common features from four approved compounds that have ... Continue reading
Ekins et al. have presented a crisp and lucid manuscript on a very relevant topic. They have suggested a methodology to extract common features from four approved compounds that have recently been found to work against the Ebola virus (amodiaquine, chloroquine, clomiphene and toremifene), and define a pharma-
cophore, which has been used to search databases, and identify further compounds for in vitro and in vivo testing. The in silico methodology described here provides an excellent method to quickly screen known compounds for possible therapies against Ebola in particular, and other viruses in general.

Response: Thank you for your constructive comments.

Some minor comments.
The result that SERMs show lower scores for binding to VP35 is rationalized by the finding that ‘clomiphene and toremifene inhibit EBOV VLP entry with some specificity to GP’¹, and therefore does not probably inhibit VP35.

Response: While the lower docking scores are noted I do not think this necessarily excludes them from inhibiting, as we know docking scores may not be that accurate and in this case, docking was used to answer the question could they fit. It’s pretty clear that a wide variety of drugs could fit based on the binding site size and accessibility.

‘Chloroquine to date has not progressed beyond the mouse EBOV model used in these studies.’ This statement is not clear, does it mean that the others have progressed beyond the mouse EBOV model?

Response: From discussions with the author on the paper that described Chloroquine as active versus EBOV in vitro and in mouse, this work has not gone beyond the mouse model of EBOV.

The first few compounds in Supplemental Table 2 should be part of the main manuscript as a table.

Response: Because this data is available easily on the website, I do not see any benefits of taking these compounds out of this supplemental table and putting them into the body of the manuscript. It might also add more confusion cutting the table up.

Color coding of pharmacophore features should be in Fig 1 too (it comes earlier than Table 1, where it is described).

Response: Thank you – this has now been added.

The structures look better with a white background.

Response: I think this is a personal preference, the structures are clear in our opinion with a black background. I have not had this suggestion previously with other publications regarding the background color.

A 3 Å RMSD for redocking a given ligand is quite high². The authors should consider the use of other docking methods, as a comparison.

Response: The study was not intended as an exhaustive docking comparison, there are plenty of these in the literature as noted by the reviewer. I agreed the redocking RMSD was high, but I also provided some justification for the result (difficulty of predicting orientations for compounds binding in what is a relatively hydrophobic and shallow pocket). If others want to use different methods and perform a comparison for this target I would be supportive.

A table of Libdock scores would help easily analyze results (with a mention of whether a higher score is better, and the significance of a score).

Response: The Libdock scores for the best poses are in the ‘4IBI Libdock docking data best poses” file. A higher score is better and this has been added to the results section
The major concern with the manuscript is the use of proprietary software, and data formats, in the study, which makes it difficult for users to probe the resultant docked structures. Further, non-standard formats are subject to the existence of the company which uses it, and not a given in the future.

Response: All of the models were generated with the proprietary software Discovery studio, and all files have been provided. The comment about such software is true, while I support using open software, I have yet to find an open source pharmacophore tool as good as that in Discovery Studio to date. It is also more convenient to use this software generating pharmacophores, receptor-ligand pharmacophores and docking in the same place. The types of analysis I have described could be repeated with any software, open source or proprietary. My hope is that by making this work openly accessible others will be inspired to pursue computational approaches with EBOV. Perhaps the community could propose the use of standards for open pharmacophore files as well. By publishing in this journal we are making all our data open even though they are in proprietary formats, I do not think this should preclude publication.
Ekins et al. have presented a crisp and lucid manuscript on a very relevant topic. They have suggested a methodology to extract common features from four approved compounds that have recently been found to work against the Ebola virus (amodiaquine, chloroquine, clomiphene and toremifene), and define a pharma-
cophore, which has been used to search databases, and identify further compounds for in vitro and in vivo testing. The in silico methodology described here provides an excellent method to quickly screen known compounds for possible therapies against Ebola in particular, and other viruses in general.

Response: Thank you for your constructive comments.

Some minor comments.
The result that SERMs show lower scores for binding to VP35 is rationalized by the finding that ‘clomiphene and toremifene inhibit EBOV VLP entry with some specificity to GP’¹, and therefore does not probably inhibit VP35.

Response: While the lower docking scores are noted I do not think this necessarily excludes them from inhibiting, as we know docking scores may not be that accurate and in this case, docking was used to answer the question could they fit. It’s pretty clear that a wide variety of drugs could fit based on the binding site size and accessibility.

‘Chloroquine to date has not progressed beyond the mouse EBOV model used in these studies.’ This statement is not clear, does it mean that the others have progressed beyond the mouse EBOV model?

Response: From discussions with the author on the paper that described Chloroquine as active versus EBOV in vitro and in mouse, this work has not gone beyond the mouse model of EBOV.

The first few compounds in Supplemental Table 2 should be part of the main manuscript as a table.

Response: Because this data is available easily on the website, I do not see any benefits of taking these compounds out of this supplemental table and putting them into the body of the manuscript. It might also add more confusion cutting the table up.

Color coding of pharmacophore features should be in Fig 1 too (it comes earlier than Table 1, where it is described).

Response: Thank you – this has now been added.

The structures look better with a white background.

Response: I think this is a personal preference, the structures are clear in our opinion with a black background. I have not had this suggestion previously with other publications regarding the background color.

A 3 Å RMSD for redocking a given ligand is quite high². The authors should consider the use of other docking methods, as a comparison.

Response: The study was not intended as an exhaustive docking comparison, there are plenty of these in the literature as noted by the reviewer. I agreed the redocking RMSD was high, but I also provided some justification for the result (difficulty of predicting orientations for compounds binding in what is a relatively hydrophobic and shallow pocket). If others want to use different methods and perform a comparison for this target I would be supportive.

A table of Libdock scores would help easily analyze results (with a mention of whether a higher score is better, and the significance of a score).

Response: The Libdock scores for the best poses are in the ‘4IBI Libdock docking data best poses” file. A higher score is better and this has been added to the results section
The major concern with the manuscript is the use of proprietary software, and data formats, in the study, which makes it difficult for users to probe the resultant docked structures. Further, non-standard formats are subject to the existence of the company which uses it, and not a given in the future.

Response: All of the models were generated with the proprietary software Discovery studio, and all files have been provided. The comment about such software is true, while I support using open software, I have yet to find an open source pharmacophore tool as good as that in Discovery Studio to date. It is also more convenient to use this software generating pharmacophores, receptor-ligand pharmacophores and docking in the same place. The types of analysis I have described could be repeated with any software, open source or proprietary. My hope is that by making this work openly accessible others will be inspired to pursue computational approaches with EBOV. Perhaps the community could propose the use of standards for open pharmacophore files as well. By publishing in this journal we are making all our data open even though they are in proprietary formats, I do not think this should preclude publication.
Competing Interests: No competing interests were disclosed. Close
Report a concern
Respond or Comment

COMMENTS ON THIS REPORT

Author Response 12 Dec 2014

Sean Ekins, Collaborative Drug Discovery, Burlingame, 94010, USA

12 Dec 2014

Author Response
Ekins et al. have presented a crisp and lucid manuscript on a very relevant topic. They have suggested a methodology to extract common features from four approved compounds that have ... Continue reading
Ekins et al. have presented a crisp and lucid manuscript on a very relevant topic. They have suggested a methodology to extract common features from four approved compounds that have recently been found to work against the Ebola virus (amodiaquine, chloroquine, clomiphene and toremifene), and define a pharma-
cophore, which has been used to search databases, and identify further compounds for in vitro and in vivo testing. The in silico methodology described here provides an excellent method to quickly screen known compounds for possible therapies against Ebola in particular, and other viruses in general.

Response: Thank you for your constructive comments.

Some minor comments.
The result that SERMs show lower scores for binding to VP35 is rationalized by the finding that ‘clomiphene and toremifene inhibit EBOV VLP entry with some specificity to GP’¹, and therefore does not probably inhibit VP35.

Response: While the lower docking scores are noted I do not think this necessarily excludes them from inhibiting, as we know docking scores may not be that accurate and in this case, docking was used to answer the question could they fit. It’s pretty clear that a wide variety of drugs could fit based on the binding site size and accessibility.

‘Chloroquine to date has not progressed beyond the mouse EBOV model used in these studies.’ This statement is not clear, does it mean that the others have progressed beyond the mouse EBOV model?

Response: From discussions with the author on the paper that described Chloroquine as active versus EBOV in vitro and in mouse, this work has not gone beyond the mouse model of EBOV.

The first few compounds in Supplemental Table 2 should be part of the main manuscript as a table.

Response: Because this data is available easily on the website, I do not see any benefits of taking these compounds out of this supplemental table and putting them into the body of the manuscript. It might also add more confusion cutting the table up.

Color coding of pharmacophore features should be in Fig 1 too (it comes earlier than Table 1, where it is described).

Response: Thank you – this has now been added.

The structures look better with a white background.

Response: I think this is a personal preference, the structures are clear in our opinion with a black background. I have not had this suggestion previously with other publications regarding the background color.

A 3 Å RMSD for redocking a given ligand is quite high². The authors should consider the use of other docking methods, as a comparison.

Response: The study was not intended as an exhaustive docking comparison, there are plenty of these in the literature as noted by the reviewer. I agreed the redocking RMSD was high, but I also provided some justification for the result (difficulty of predicting orientations for compounds binding in what is a relatively hydrophobic and shallow pocket). If others want to use different methods and perform a comparison for this target I would be supportive.

A table of Libdock scores would help easily analyze results (with a mention of whether a higher score is better, and the significance of a score).

Response: The Libdock scores for the best poses are in the ‘4IBI Libdock docking data best poses” file. A higher score is better and this has been added to the results section
The major concern with the manuscript is the use of proprietary software, and data formats, in the study, which makes it difficult for users to probe the resultant docked structures. Further, non-standard formats are subject to the existence of the company which uses it, and not a given in the future.

Response: All of the models were generated with the proprietary software Discovery studio, and all files have been provided. The comment about such software is true, while I support using open software, I have yet to find an open source pharmacophore tool as good as that in Discovery Studio to date. It is also more convenient to use this software generating pharmacophores, receptor-ligand pharmacophores and docking in the same place. The types of analysis I have described could be repeated with any software, open source or proprietary. My hope is that by making this work openly accessible others will be inspired to pursue computational approaches with EBOV. Perhaps the community could propose the use of standards for open pharmacophore files as well. By publishing in this journal we are making all our data open even though they are in proprietary formats, I do not think this should preclude publication.
Ekins et al. have presented a crisp and lucid manuscript on a very relevant topic. They have suggested a methodology to extract common features from four approved compounds that have recently been found to work against the Ebola virus (amodiaquine, chloroquine, clomiphene and toremifene), and define a pharma-
cophore, which has been used to search databases, and identify further compounds for in vitro and in vivo testing. The in silico methodology described here provides an excellent method to quickly screen known compounds for possible therapies against Ebola in particular, and other viruses in general.

Response: Thank you for your constructive comments.

Some minor comments.
The result that SERMs show lower scores for binding to VP35 is rationalized by the finding that ‘clomiphene and toremifene inhibit EBOV VLP entry with some specificity to GP’¹, and therefore does not probably inhibit VP35.

Response: While the lower docking scores are noted I do not think this necessarily excludes them from inhibiting, as we know docking scores may not be that accurate and in this case, docking was used to answer the question could they fit. It’s pretty clear that a wide variety of drugs could fit based on the binding site size and accessibility.

‘Chloroquine to date has not progressed beyond the mouse EBOV model used in these studies.’ This statement is not clear, does it mean that the others have progressed beyond the mouse EBOV model?

Response: From discussions with the author on the paper that described Chloroquine as active versus EBOV in vitro and in mouse, this work has not gone beyond the mouse model of EBOV.

The first few compounds in Supplemental Table 2 should be part of the main manuscript as a table.

Response: Because this data is available easily on the website, I do not see any benefits of taking these compounds out of this supplemental table and putting them into the body of the manuscript. It might also add more confusion cutting the table up.

Color coding of pharmacophore features should be in Fig 1 too (it comes earlier than Table 1, where it is described).

Response: Thank you – this has now been added.

The structures look better with a white background.

Response: I think this is a personal preference, the structures are clear in our opinion with a black background. I have not had this suggestion previously with other publications regarding the background color.

A 3 Å RMSD for redocking a given ligand is quite high². The authors should consider the use of other docking methods, as a comparison.

Response: The study was not intended as an exhaustive docking comparison, there are plenty of these in the literature as noted by the reviewer. I agreed the redocking RMSD was high, but I also provided some justification for the result (difficulty of predicting orientations for compounds binding in what is a relatively hydrophobic and shallow pocket). If others want to use different methods and perform a comparison for this target I would be supportive.

A table of Libdock scores would help easily analyze results (with a mention of whether a higher score is better, and the significance of a score).

Response: The Libdock scores for the best poses are in the ‘4IBI Libdock docking data best poses” file. A higher score is better and this has been added to the results section
The major concern with the manuscript is the use of proprietary software, and data formats, in the study, which makes it difficult for users to probe the resultant docked structures. Further, non-standard formats are subject to the existence of the company which uses it, and not a given in the future.

Response: All of the models were generated with the proprietary software Discovery studio, and all files have been provided. The comment about such software is true, while I support using open software, I have yet to find an open source pharmacophore tool as good as that in Discovery Studio to date. It is also more convenient to use this software generating pharmacophores, receptor-ligand pharmacophores and docking in the same place. The types of analysis I have described could be repeated with any software, open source or proprietary. My hope is that by making this work openly accessible others will be inspired to pursue computational approaches with EBOV. Perhaps the community could propose the use of standards for open pharmacophore files as well. By publishing in this journal we are making all our data open even though they are in proprietary formats, I do not think this should preclude publication.
Competing Interests: No competing interests were disclosed. Close
Report a concern

Views

Reviewer Report 01 Dec 2014

Gloria Fuentes, Center for Life Sciences Technologies, RIKEN, Saitama, Japan

Approved

https://doi.org/10.5256/f1000research.6137.r6707

The current Ebola crisis in West Africa has shattered all expectations by continuing to grow months following the initial case. This has stimulated a massive and global emergency response, and it has challenged the health protocols and by extension, the efforts of scientific community.

The authors have carried out a computational analysis using several compounds detected in two previous high throughput screens to build a pharmacophore model. The key features of such a model are used to scan databases of small molecules. They have come up with a list of putative inhibitors. Their study will have a stronger scientific impact if the authors could elaborate more in suggestions on how the best ranked compounds will increase the binding affinity, based in the structural model VP30-inhibitors that they have built.

In parallel, they observed a highly overlapping between the motifs in the pharmacophore and those found in the crystal structures of several inhibitors of the viral protein 35. Based on this fact, and in their results in an in-silico docking, they propose that the most likely inhibitory mechanism for these compounds is the targeting of the protein-protein interaction involving this protein.

In this regard, the authors should extend their study to include different docking protocols, including different programs, in an attempt to verify their results. As they mention in the text (page 4), the redocking of the ligand in 4IBI to the protein does not show the crystal structure binding mode accurately. These different settings could help in a better prediction of the ligand orientations.

Concerning the docking and proposed mechanism, I wonder how different the other solved nucleocapsid proteins are from a structural and sequence point of view, in order to make the authors point that VP35 is indeed the target. Would it be possible to explore for surface patches with similar physico-chemical features? In the case of VP30, a potential binding pocket for small-molecule inhibitors has been suggested by Hartlieb et al. (2007). How good or bad the overlapping with the built pharmachophore is for this case?

To complete the structural understanding of the action of these compounds, a figure displaying their location on the protein surface as well as the binding site for RNA would clarify their role in the inhibition of protein-protein interactions.

In summary, the manuscript describes an interesting and fast approach to identify putative inhibitors for a currently serious target as Ebola virus. Although their results should be experimental validated to confirm their finding, this computational study and further extensions of it are of the great value.

Competing Interests: No competing interests were disclosed.

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

CITE

Report a concern

Author Response 12 Dec 2014

Sean Ekins, Collaborative Drug Discovery, Burlingame, 94010, USA

12 Dec 2014

Author Response
- The current Ebola crisis in West Africa has shattered all expectations by continuing to grow months following the initial case. This has stimulated a massive and global emergency response, and
... Continue reading
The current Ebola crisis in West Africa has shattered all expectations by continuing to grow months following the initial case. This has stimulated a massive and global emergency response, and it has challenged the health protocols and by extension, the efforts of scientific community.

The authors have carried out a computational analysis using several compounds detected in two previous high throughput screens to build a pharmacophore model. The key features of such a model are used to scan databases of small molecules. They have come up with a list of putative inhibitors. Their study will have a stronger scientific impact if the authors could elaborate more in suggestions on how the best ranked compounds will increase the binding affinity, based in the structural model VP30-inhibitors that they have built.

Response: To clarify, we have focused on VP35 not VP30. I am not aware of an X-ray structure with ligand bound for VP30. The same type of approach could certainly be pursued with other EBOV targets. We produced a common features pharmacophore for the 4 compounds, and after looking at the VP35 receptor-ligand pharmacophores proposed that there may be some overlap, and then this led to docking the 4 compounds in the X-ray structures. Our intent was not to design molecules but to use the available methods to perhaps infer a potential target/mechanism and then perhaps researchers would want to test the compounds. We do not have access to experimentally test these predictions, but this manuscript may lead to others doing this work perhaps. Whether one wants to use the Libdock score for (absolute) prediction of binding affinity interactions is debatable, rather this approach might help to limit or prioritize which compounds to test.

In parallel, they observed a highly overlapping between the motifs in the pharmacophore and those found in the crystal structures of several inhibitors of the viral protein 35. Based on this fact, and in their results in an in-silico docking, they propose that the most likely inhibitory mechanism for these compounds is the targeting of the protein-protein interaction involving this protein.

Response: VP35 may be a target for these compounds although we do not discount other targets or non-target related mechanisms.

In this regard, the authors should extend their study to include different docking protocols, including different programs, in an attempt to verify their results. As they mention in the text (page 4), the redocking of the ligand in 4IBI to the protein does not show the crystal structure binding mode accurately. These different settings could help in a better prediction of the ligand orientations.

Response: As explained earlier, our study is not intended to be an exhaustive evaluation of docking tools, we have used different computational approaches to suggest that the FDA drugs may have a common pharmacophore, which seems to be similar to that of the ligands co-crystallized with VP-35. Finally docking suggests they may fit into the pocket that the co-crystal ligands bind to. The work proposes that the compounds could fit in the binding site, but it is unclear what additional value more docking would add unless we were going to try to predict and then generate the X-ray structure of these FDA drugs. Certainly if experts in docking or crystallography want to pursue this target they can.

Concerning the docking and proposed mechanism, I wonder how different the other solved nucleocapsid proteins are from a structural and sequence point of view, in order to make the authors point that VP35 is indeed the target. Would it be possible to explore for surface patches with similar physico-chemical features? In the case of VP30, a potential binding pocket for small-molecule inhibitors has been suggested by Hartlieb et al. (2007). How good or bad the overlapping with the built pharmachophore is for this case?

Response: this is indeed a very good point. We are not experts on these proteins. I think the proposed work could be done, the difficulty may be that there is no crystal structure (that I can see) with a ligand bound that would be a useful guide to binding in this pocket and would be essential for a receptor-ligand pharmacophore to be built.

To complete the structural understanding of the action of these compounds, a figure displaying their location on the protein surface as well as the binding site for RNA would clarify their role in the inhibition of protein-protein interactions.

Response: we have now added Figure S2 which shows the molecules in the context of the full protein. They are in the site suggested in ref 15 as important for the nucleoprotein interaction.

In summary, the manuscript describes an interesting and fast approach to identify putative inhibitors for a currently serious target as Ebola virus. Although their results should be experimental validated to confirm their finding, this computational study and further extensions of it are of the great value.

Response: Thank you for your suggestions, we agree and would encourage other scientists to test whether these compounds are targeting VP35 or VP 30 as you propose, or having an alternative mechanism. I think we would also be happy to see any of these molecules progress into other animal models of EBOV.
The current Ebola crisis in West Africa has shattered all expectations by continuing to grow months following the initial case. This has stimulated a massive and global emergency response, and it has challenged the health protocols and by extension, the efforts of scientific community.

The authors have carried out a computational analysis using several compounds detected in two previous high throughput screens to build a pharmacophore model. The key features of such a model are used to scan databases of small molecules. They have come up with a list of putative inhibitors. Their study will have a stronger scientific impact if the authors could elaborate more in suggestions on how the best ranked compounds will increase the binding affinity, based in the structural model VP30-inhibitors that they have built.

Response: To clarify, we have focused on VP35 not VP30. I am not aware of an X-ray structure with ligand bound for VP30. The same type of approach could certainly be pursued with other EBOV targets. We produced a common features pharmacophore for the 4 compounds, and after looking at the VP35 receptor-ligand pharmacophores proposed that there may be some overlap, and then this led to docking the 4 compounds in the X-ray structures. Our intent was not to design molecules but to use the available methods to perhaps infer a potential target/mechanism and then perhaps researchers would want to test the compounds. We do not have access to experimentally test these predictions, but this manuscript may lead to others doing this work perhaps. Whether one wants to use the Libdock score for (absolute) prediction of binding affinity interactions is debatable, rather this approach might help to limit or prioritize which compounds to test.

In parallel, they observed a highly overlapping between the motifs in the pharmacophore and those found in the crystal structures of several inhibitors of the viral protein 35. Based on this fact, and in their results in an in-silico docking, they propose that the most likely inhibitory mechanism for these compounds is the targeting of the protein-protein interaction involving this protein.

Response: VP35 may be a target for these compounds although we do not discount other targets or non-target related mechanisms.

In this regard, the authors should extend their study to include different docking protocols, including different programs, in an attempt to verify their results. As they mention in the text (page 4), the redocking of the ligand in 4IBI to the protein does not show the crystal structure binding mode accurately. These different settings could help in a better prediction of the ligand orientations.

Response: As explained earlier, our study is not intended to be an exhaustive evaluation of docking tools, we have used different computational approaches to suggest that the FDA drugs may have a common pharmacophore, which seems to be similar to that of the ligands co-crystallized with VP-35. Finally docking suggests they may fit into the pocket that the co-crystal ligands bind to. The work proposes that the compounds could fit in the binding site, but it is unclear what additional value more docking would add unless we were going to try to predict and then generate the X-ray structure of these FDA drugs. Certainly if experts in docking or crystallography want to pursue this target they can.

Concerning the docking and proposed mechanism, I wonder how different the other solved nucleocapsid proteins are from a structural and sequence point of view, in order to make the authors point that VP35 is indeed the target. Would it be possible to explore for surface patches with similar physico-chemical features? In the case of VP30, a potential binding pocket for small-molecule inhibitors has been suggested by Hartlieb et al. (2007). How good or bad the overlapping with the built pharmachophore is for this case?

Response: this is indeed a very good point. We are not experts on these proteins. I think the proposed work could be done, the difficulty may be that there is no crystal structure (that I can see) with a ligand bound that would be a useful guide to binding in this pocket and would be essential for a receptor-ligand pharmacophore to be built.

To complete the structural understanding of the action of these compounds, a figure displaying their location on the protein surface as well as the binding site for RNA would clarify their role in the inhibition of protein-protein interactions.

Response: we have now added Figure S2 which shows the molecules in the context of the full protein. They are in the site suggested in ref 15 as important for the nucleoprotein interaction.

In summary, the manuscript describes an interesting and fast approach to identify putative inhibitors for a currently serious target as Ebola virus. Although their results should be experimental validated to confirm their finding, this computational study and further extensions of it are of the great value.

Response: Thank you for your suggestions, we agree and would encourage other scientists to test whether these compounds are targeting VP35 or VP 30 as you propose, or having an alternative mechanism. I think we would also be happy to see any of these molecules progress into other animal models of EBOV.
Competing Interests: No competing interests were disclosed. Close
Report a concern
Respond or Comment

COMMENTS ON THIS REPORT

Author Response 12 Dec 2014

Sean Ekins, Collaborative Drug Discovery, Burlingame, 94010, USA

12 Dec 2014

Author Response
- The current Ebola crisis in West Africa has shattered all expectations by continuing to grow months following the initial case. This has stimulated a massive and global emergency response, and
... Continue reading
The current Ebola crisis in West Africa has shattered all expectations by continuing to grow months following the initial case. This has stimulated a massive and global emergency response, and it has challenged the health protocols and by extension, the efforts of scientific community.

The authors have carried out a computational analysis using several compounds detected in two previous high throughput screens to build a pharmacophore model. The key features of such a model are used to scan databases of small molecules. They have come up with a list of putative inhibitors. Their study will have a stronger scientific impact if the authors could elaborate more in suggestions on how the best ranked compounds will increase the binding affinity, based in the structural model VP30-inhibitors that they have built.

Response: To clarify, we have focused on VP35 not VP30. I am not aware of an X-ray structure with ligand bound for VP30. The same type of approach could certainly be pursued with other EBOV targets. We produced a common features pharmacophore for the 4 compounds, and after looking at the VP35 receptor-ligand pharmacophores proposed that there may be some overlap, and then this led to docking the 4 compounds in the X-ray structures. Our intent was not to design molecules but to use the available methods to perhaps infer a potential target/mechanism and then perhaps researchers would want to test the compounds. We do not have access to experimentally test these predictions, but this manuscript may lead to others doing this work perhaps. Whether one wants to use the Libdock score for (absolute) prediction of binding affinity interactions is debatable, rather this approach might help to limit or prioritize which compounds to test.

In parallel, they observed a highly overlapping between the motifs in the pharmacophore and those found in the crystal structures of several inhibitors of the viral protein 35. Based on this fact, and in their results in an in-silico docking, they propose that the most likely inhibitory mechanism for these compounds is the targeting of the protein-protein interaction involving this protein.

Response: VP35 may be a target for these compounds although we do not discount other targets or non-target related mechanisms.

In this regard, the authors should extend their study to include different docking protocols, including different programs, in an attempt to verify their results. As they mention in the text (page 4), the redocking of the ligand in 4IBI to the protein does not show the crystal structure binding mode accurately. These different settings could help in a better prediction of the ligand orientations.

Response: As explained earlier, our study is not intended to be an exhaustive evaluation of docking tools, we have used different computational approaches to suggest that the FDA drugs may have a common pharmacophore, which seems to be similar to that of the ligands co-crystallized with VP-35. Finally docking suggests they may fit into the pocket that the co-crystal ligands bind to. The work proposes that the compounds could fit in the binding site, but it is unclear what additional value more docking would add unless we were going to try to predict and then generate the X-ray structure of these FDA drugs. Certainly if experts in docking or crystallography want to pursue this target they can.

Concerning the docking and proposed mechanism, I wonder how different the other solved nucleocapsid proteins are from a structural and sequence point of view, in order to make the authors point that VP35 is indeed the target. Would it be possible to explore for surface patches with similar physico-chemical features? In the case of VP30, a potential binding pocket for small-molecule inhibitors has been suggested by Hartlieb et al. (2007). How good or bad the overlapping with the built pharmachophore is for this case?

Response: this is indeed a very good point. We are not experts on these proteins. I think the proposed work could be done, the difficulty may be that there is no crystal structure (that I can see) with a ligand bound that would be a useful guide to binding in this pocket and would be essential for a receptor-ligand pharmacophore to be built.

To complete the structural understanding of the action of these compounds, a figure displaying their location on the protein surface as well as the binding site for RNA would clarify their role in the inhibition of protein-protein interactions.

Response: we have now added Figure S2 which shows the molecules in the context of the full protein. They are in the site suggested in ref 15 as important for the nucleoprotein interaction.

In summary, the manuscript describes an interesting and fast approach to identify putative inhibitors for a currently serious target as Ebola virus. Although their results should be experimental validated to confirm their finding, this computational study and further extensions of it are of the great value.

Response: Thank you for your suggestions, we agree and would encourage other scientists to test whether these compounds are targeting VP35 or VP 30 as you propose, or having an alternative mechanism. I think we would also be happy to see any of these molecules progress into other animal models of EBOV.
The current Ebola crisis in West Africa has shattered all expectations by continuing to grow months following the initial case. This has stimulated a massive and global emergency response, and it has challenged the health protocols and by extension, the efforts of scientific community.

The authors have carried out a computational analysis using several compounds detected in two previous high throughput screens to build a pharmacophore model. The key features of such a model are used to scan databases of small molecules. They have come up with a list of putative inhibitors. Their study will have a stronger scientific impact if the authors could elaborate more in suggestions on how the best ranked compounds will increase the binding affinity, based in the structural model VP30-inhibitors that they have built.

Response: To clarify, we have focused on VP35 not VP30. I am not aware of an X-ray structure with ligand bound for VP30. The same type of approach could certainly be pursued with other EBOV targets. We produced a common features pharmacophore for the 4 compounds, and after looking at the VP35 receptor-ligand pharmacophores proposed that there may be some overlap, and then this led to docking the 4 compounds in the X-ray structures. Our intent was not to design molecules but to use the available methods to perhaps infer a potential target/mechanism and then perhaps researchers would want to test the compounds. We do not have access to experimentally test these predictions, but this manuscript may lead to others doing this work perhaps. Whether one wants to use the Libdock score for (absolute) prediction of binding affinity interactions is debatable, rather this approach might help to limit or prioritize which compounds to test.

In parallel, they observed a highly overlapping between the motifs in the pharmacophore and those found in the crystal structures of several inhibitors of the viral protein 35. Based on this fact, and in their results in an in-silico docking, they propose that the most likely inhibitory mechanism for these compounds is the targeting of the protein-protein interaction involving this protein.

Response: VP35 may be a target for these compounds although we do not discount other targets or non-target related mechanisms.

In this regard, the authors should extend their study to include different docking protocols, including different programs, in an attempt to verify their results. As they mention in the text (page 4), the redocking of the ligand in 4IBI to the protein does not show the crystal structure binding mode accurately. These different settings could help in a better prediction of the ligand orientations.

Response: As explained earlier, our study is not intended to be an exhaustive evaluation of docking tools, we have used different computational approaches to suggest that the FDA drugs may have a common pharmacophore, which seems to be similar to that of the ligands co-crystallized with VP-35. Finally docking suggests they may fit into the pocket that the co-crystal ligands bind to. The work proposes that the compounds could fit in the binding site, but it is unclear what additional value more docking would add unless we were going to try to predict and then generate the X-ray structure of these FDA drugs. Certainly if experts in docking or crystallography want to pursue this target they can.

Concerning the docking and proposed mechanism, I wonder how different the other solved nucleocapsid proteins are from a structural and sequence point of view, in order to make the authors point that VP35 is indeed the target. Would it be possible to explore for surface patches with similar physico-chemical features? In the case of VP30, a potential binding pocket for small-molecule inhibitors has been suggested by Hartlieb et al. (2007). How good or bad the overlapping with the built pharmachophore is for this case?

Response: this is indeed a very good point. We are not experts on these proteins. I think the proposed work could be done, the difficulty may be that there is no crystal structure (that I can see) with a ligand bound that would be a useful guide to binding in this pocket and would be essential for a receptor-ligand pharmacophore to be built.

To complete the structural understanding of the action of these compounds, a figure displaying their location on the protein surface as well as the binding site for RNA would clarify their role in the inhibition of protein-protein interactions.

Response: we have now added Figure S2 which shows the molecules in the context of the full protein. They are in the site suggested in ref 15 as important for the nucleoprotein interaction.

In summary, the manuscript describes an interesting and fast approach to identify putative inhibitors for a currently serious target as Ebola virus. Although their results should be experimental validated to confirm their finding, this computational study and further extensions of it are of the great value.

Response: Thank you for your suggestions, we agree and would encourage other scientists to test whether these compounds are targeting VP35 or VP 30 as you propose, or having an alternative mechanism. I think we would also be happy to see any of these molecules progress into other animal models of EBOV.
Competing Interests: No competing interests were disclosed. Close
Report a concern

Comments on this article Comments (0)

Version 2

VERSION 2 PUBLISHED 14 Nov 2014

Open Peer Review

Reviewer Status

Reviewer Reports

	Invited Reviewers
	1	2
Version 2 (revision) 12 Dec 14
Version 1 14 Nov 14	read	read

Gloria Fuentes, RIKEN, Saitama, Japan
Sandeep Chakraborty, University of California, Davis, USA

Comments on this article

All Comments(0)

Add a comment

Browse by related subjects

Back to all reports

Reviewer Report

42 Views

02 Dec 2014 | for Version 1

Sandeep Chakraborty, Plant Sciences Department, University of California, Davis, CA, USA

42 Views Cite this report Responses(1)

Approved

The result that SERMs show lower scores for binding to VP35 is rationalized by the finding that ‘clomiphene and toremifene inhibit EBOV VLP entry with some specificity to GP’¹, and therefore does not probably inhibit VP35.
‘Chloroquine to date has not progressed beyond the mouse EBOV model used in these studies.’ This statement is not clear, does it mean that the others have progressed beyond the mouse EBOV model?
The first few compounds in Supplemental Table 2 should be part of the main manuscript as a table.
Color coding of pharmacophore features should be in Fig 1 too (it comes earlier than Table 1, where it is described).
The structures look better with a white background.
A 3 Å RMSD for redocking a given ligand is quite high². The authors should consider the use of other docking methods, as a comparison.
A table of Libdock scores would help easily analyze results (with a mention of whether a higher score is better, and the significance of a score).

References

Competing Interests

No competing interests were disclosed.

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

Respond to this report

Responses (1)

Author Response

12 Dec 2014

Sean Ekins, Collaborative Drug Discovery, Burlingame, 94010, USA

Ekins et al. have presented a crisp and lucid manuscript on a very relevant topic. They have suggested a methodology to extract common features from four approved compounds that have recently been found to work against the Ebola virus (amodiaquine, chloroquine, clomiphene and toremifene), and define a pharma-
cophore, which has been used to search databases, and identify further compounds for in vitro and in vivo testing. The in silico methodology described here provides an excellent method to quickly screen known compounds for possible therapies against Ebola in particular, and other viruses in general.

Response: Thank you for your constructive comments.

Some minor comments.

The result that SERMs show lower scores for binding to VP35 is rationalized by the finding that ‘clomiphene and toremifene inhibit EBOV VLP entry with some specificity to GP’¹, and therefore does not probably inhibit VP35.

Response: While the lower docking scores are noted I do not think this necessarily excludes them from inhibiting, as we know docking scores may not be that accurate and in this case, docking was used to answer the question could they fit. It’s pretty clear that a wide variety of drugs could fit based on the binding site size and accessibility.
‘Chloroquine to date has not progressed beyond the mouse EBOV model used in these studies.’ This statement is not clear, does it mean that the others have progressed beyond the mouse EBOV model?

Response: From discussions with the author on the paper that described Chloroquine as active versus EBOV in vitro and in mouse, this work has not gone beyond the mouse model of EBOV.
The first few compounds in Supplemental Table 2 should be part of the main manuscript as a table.

Response: Because this data is available easily on the website, I do not see any benefits of taking these compounds out of this supplemental table and putting them into the body of the manuscript. It might also add more confusion cutting the table up.
Color coding of pharmacophore features should be in Fig 1 too (it comes earlier than Table 1, where it is described).

Response: Thank you – this has now been added.
The structures look better with a white background.

Response: I think this is a personal preference, the structures are clear in our opinion with a black background. I have not had this suggestion previously with other publications regarding the background color.
A 3 Å RMSD for redocking a given ligand is quite high². The authors should consider the use of other docking methods, as a comparison.

Response: The study was not intended as an exhaustive docking comparison, there are plenty of these in the literature as noted by the reviewer. I agreed the redocking RMSD was high, but I also provided some justification for the result (difficulty of predicting orientations for compounds binding in what is a relatively hydrophobic and shallow pocket). If others want to use different methods and perform a comparison for this target I would be supportive.
A table of Libdock scores would help easily analyze results (with a mention of whether a higher score is better, and the significance of a score).

Response: The Libdock scores for the best poses are in the ‘4IBI Libdock docking data best poses” file. A higher score is better and this has been added to the results section

The major concern with the manuscript is the use of proprietary software, and data formats, in the study, which makes it difficult for users to probe the resultant docked structures. Further, non-standard formats are subject to the existence of the company which uses it, and not a given in the future.

Response: All of the models were generated with the proprietary software Discovery studio, and all files have been provided. The comment about such software is true, while I support using open software, I have yet to find an open source pharmacophore tool as good as that in Discovery Studio to date. It is also more convenient to use this software generating pharmacophores, receptor-ligand pharmacophores and docking in the same place. The types of analysis I have described could be repeated with any software, open source or proprietary. My hope is that by making this work openly accessible others will be inspired to pursue computational approaches with EBOV. Perhaps the community could propose the use of standards for open pharmacophore files as well. By publishing in this journal we are making all our data open even though they are in proprietary formats, I do not think this should preclude publication.

View more View less

Competing Interests

No competing interests were disclosed.

Back to all reports

Reviewer Report

39 Views

01 Dec 2014 | for Version 1

Gloria Fuentes, Center for Life Sciences Technologies, RIKEN, Saitama, Japan

39 Views Cite this report Responses(1)

Approved

Competing Interests

No competing interests were disclosed.

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

Respond to this report

Responses (1)

Author Response

12 Dec 2014

Sean Ekins, Collaborative Drug Discovery, Burlingame, 94010, USA

The current Ebola crisis in West Africa has shattered all expectations by continuing to grow months following the initial case. This has stimulated a massive and global emergency response, and it has challenged the health protocols and by extension, the efforts of scientific community.

The authors have carried out a computational analysis using several compounds detected in two previous high throughput screens to build a pharmacophore model. The key features of such a model are used to scan databases of small molecules. They have come up with a list of putative inhibitors. Their study will have a stronger scientific impact if the authors could elaborate more in suggestions on how the best ranked compounds will increase the binding affinity, based in the structural model VP30-inhibitors that they have built.

Response: To clarify, we have focused on VP35 not VP30. I am not aware of an X-ray structure with ligand bound for VP30. The same type of approach could certainly be pursued with other EBOV targets. We produced a common features pharmacophore for the 4 compounds, and after looking at the VP35 receptor-ligand pharmacophores proposed that there may be some overlap, and then this led to docking the 4 compounds in the X-ray structures. Our intent was not to design molecules but to use the available methods to perhaps infer a potential target/mechanism and then perhaps researchers would want to test the compounds. We do not have access to experimentally test these predictions, but this manuscript may lead to others doing this work perhaps. Whether one wants to use the Libdock score for (absolute) prediction of binding affinity interactions is debatable, rather this approach might help to limit or prioritize which compounds to test.
In parallel, they observed a highly overlapping between the motifs in the pharmacophore and those found in the crystal structures of several inhibitors of the viral protein 35. Based on this fact, and in their results in an in-silico docking, they propose that the most likely inhibitory mechanism for these compounds is the targeting of the protein-protein interaction involving this protein.

Response: VP35 may be a target for these compounds although we do not discount other targets or non-target related mechanisms.
In this regard, the authors should extend their study to include different docking protocols, including different programs, in an attempt to verify their results. As they mention in the text (page 4), the redocking of the ligand in 4IBI to the protein does not show the crystal structure binding mode accurately. These different settings could help in a better prediction of the ligand orientations.

Response: As explained earlier, our study is not intended to be an exhaustive evaluation of docking tools, we have used different computational approaches to suggest that the FDA drugs may have a common pharmacophore, which seems to be similar to that of the ligands co-crystallized with VP-35. Finally docking suggests they may fit into the pocket that the co-crystal ligands bind to. The work proposes that the compounds could fit in the binding site, but it is unclear what additional value more docking would add unless we were going to try to predict and then generate the X-ray structure of these FDA drugs. Certainly if experts in docking or crystallography want to pursue this target they can.
Concerning the docking and proposed mechanism, I wonder how different the other solved nucleocapsid proteins are from a structural and sequence point of view, in order to make the authors point that VP35 is indeed the target. Would it be possible to explore for surface patches with similar physico-chemical features? In the case of VP30, a potential binding pocket for small-molecule inhibitors has been suggested by Hartlieb et al. (2007). How good or bad the overlapping with the built pharmachophore is for this case?

Response: this is indeed a very good point. We are not experts on these proteins. I think the proposed work could be done, the difficulty may be that there is no crystal structure (that I can see) with a ligand bound that would be a useful guide to binding in this pocket and would be essential for a receptor-ligand pharmacophore to be built.
To complete the structural understanding of the action of these compounds, a figure displaying their location on the protein surface as well as the binding site for RNA would clarify their role in the inhibition of protein-protein interactions.

Response: we have now added Figure S2 which shows the molecules in the context of the full protein. They are in the site suggested in ref 15 as important for the nucleoprotein interaction.
In summary, the manuscript describes an interesting and fast approach to identify putative inhibitors for a currently serious target as Ebola virus. Although their results should be experimental validated to confirm their finding, this computational study and further extensions of it are of the great value.

Response: Thank you for your suggestions, we agree and would encourage other scientists to test whether these compounds are targeting VP35 or VP 30 as you propose, or having an alternative mechanism. I think we would also be happy to see any of these molecules progress into other animal models of EBOV.

View more View less

Competing Interests

No competing interests were disclosed.

Alongside their report, reviewers assign a status to the article:

Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested

Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.

Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions

Click here to access the data.

Downloaded data do not display as expected? Download the data

[1] 1. Butler D, Morello L: Ebola by the numbers: The size, spread and cost of an outbreak. Nature. 2014; 514(7522): 284–5. PubMed Abstract | Publisher Full Text

[2] 2. Piot P: Ebola’s perfect storm. Science. 2014; 345(6202): 1221. PubMed Abstract | Publisher Full Text

[3] 3. Trad MA, Fisher DA, Tambyah PA: Ebola in west Africa. Lancet Infect Dis. 2014; 14(11): 1045. PubMed Abstract | Publisher Full Text

[4] 4. Dhillon RS, Srikrishna D, Sachs J: Controlling Ebola: next steps. Lancet. 2014; 384(9952): 1409–1411. PubMed Abstract | Publisher Full Text

[5] 5. Anon: Call to action. Nature. 2014; 514(7524): 535–536. PubMed Abstract | Publisher Full Text

[6] 6. Check Hayden E: The Ebola questions. Nature. 2014; 514(7524): 554–7. PubMed Abstract | Publisher Full Text

[7] 7. Piot P: The F1000Research: Ebola article collection [v1; ref status: not peer reviewed, http://f1000r.es/4ot]. F1000Res. 2014; 3: 269. Publisher Full Text

[8] 8. Madrid PB, Chopra S, Manger ID, et al.: A systematic screen of FDA-approved drugs for inhibitors of biological threat agents. PLoS One. 2013; 8(4): e60579. PubMed Abstract | Publisher Full Text | Free Full Text

[9] 9. Johansen LM, Brannan JM, Delos SE, et al.: FDA-approved selective estrogen receptor modulators inhibit Ebola virus infection. Sci Transl Med. 2013; 5(190): 190ra79. PubMed Abstract | Publisher Full Text | Free Full Text

[10] 10. Keyaerts E, Li S, Vijgen L, et al.: Antiviral activity of chloroquine against human coronavirus OC43 infection in newborn mice. Antimicrob Agents Chemother. 2009; 53(8): 3416–21. PubMed Abstract | Publisher Full Text | Free Full Text

[11] 11. Boonyasuppayakorn S, Reichert ED, Manzano M, et al.: Amodiaquine, an antimalarial drug, inhibits dengue virus type 2 replication and infectivity. Antiviral Res. 2014; 106: 125–34. PubMed Abstract | Publisher Full Text

[12] 12. Acosta EG, Bruttomesso AC, Bisceglia JA, et al.: Dehydroepiandrosterone, epiandrosterone and synthetic derivatives inhibit Junin virus replication in vitro. Virus Res. 2008; 135(2): 203–12. PubMed Abstract | Publisher Full Text

[13] 13. Dyall J, Coleman CM, Hart BJ, et al.: Repurposing of clinically developed drugs for treatment of Middle East respiratory syndrome coronavirus infection. Antimicrob Agents Chemother. 2014; 58(8): 4885–93. PubMed Abstract | Publisher Full Text | Free Full Text

[14] 14. Han Z, Lu J, Liu Y, et al.: Small-molecule probes targeting the viral PPxY-host Nedd4 interface block egress of a broad range of RNA viruses. J Virol. 2014; 88(13): 7294–306. PubMed Abstract | Publisher Full Text | Free Full Text

[15] 15. Brown CS, Lee MS, Leung DW, et al.: In silico derived small molecules bind the filovirus VP35 protein and inhibit its polymerase cofactor activity. J Mol Biol. 2014; 426(10): 2045–58. PubMed Abstract | Publisher Full Text | Free Full Text

[16] 16. Chakraborty S, Rao B, Asgeirsson B, et al.: Characterizing alpha helical properties of Ebola viral proteins as potential targets for inhibition of alpha-helix mediated protein-protein interactions [v1; ref status: approved with reservations 1, http://f1000r.es/4lg]. F1000Res. 2014; 3: 251. Publisher Full Text

[17] 17. Lamichhane G, Freundlich JS, Ekins S, et al.: Essential metabolites of Mycobacterium tuberculosis and their mimics. MBio. 2011; 2(1): e00301–10. PubMed Abstract | Publisher Full Text | Free Full Text

[18] 18. Ekins S, Bradford J, Dole K, et al.: A collaborative database and computational models for tuberculosis drug discovery. Mol Biosyst. 2010; 6(5): 840–851. PubMed Abstract | Publisher Full Text

[19] 19. Zheng X, Ekins S, Raufman JP, et al.: Computational models for drug inhibition of the human apical sodium-dependent bile acid transporter. Mol Pharm. 2009; 6(5): 1591–1603. PubMed Abstract | Publisher Full Text | Free Full Text

[20] 20. Meslamani J, Li J, Sutter J, et al.: Protein-ligand-based pharmacophores: generation and utility assessment in computational ligand profiling. J Chem Inf Model. 2012; 52(4): 943–55. PubMed Abstract | Publisher Full Text

[21] 21. Rao SN, Head MS, Kulkarni A, et al.: Validation studies of the site-directed docking program LibDock. J Chem Inf Model. 2007; 47(6): 2159–71. PubMed Abstract | Publisher Full Text

[22] 22. Diao L, Ekins S, Polli JE: Novel inhibitors of human organic cation/carnitine transporter (hOCTN2) via computational modeling and in vitro testing. Pharm Res. 2009; 26(8): 1890–1900. PubMed Abstract | Publisher Full Text | Free Full Text

[23] 23. Zheng X, Ekins S, Raufman JP, et al.: Computational models for drug inhibition of the human apical sodium-dependent bile acid transporter. Mol Pharm. 2009; 6(5): 1591–603. PubMed Abstract | Publisher Full Text | Free Full Text

[24] 24. Diao L, Ekins S, Polli JE: Quantitative structure activity relationship for inhibition of human organic cation/carnitine transporter. Mol Pharm. 2010; 7(6): 2120–2131. PubMed Abstract | Publisher Full Text | Free Full Text

[25] 25. Ekins S, Diao L, Polli JE: A substrate pharmacophore for the human organic cation/carnitine transporter identifies compounds associated with rhabdomyolysis. Mol Pharm. 2012; 9(4): 905–913. PubMed Abstract | Publisher Full Text | Free Full Text

[26] 26. Dong Z, Ekins S, Polli JE: Structure-activity relationship for FDA approved drugs as inhibitors of the human sodium taurocholate cotransporting polypeptide (NTCP). Mol Pharm. 2013; 10(3): 1008–19. PubMed Abstract | Publisher Full Text | Free Full Text

[27] 27. Dong Z, Ekins S, Polli JE: Quantitative NTCP pharmacophore and lack of association between DILI and NTCP Inhibition. Eur J Pharm Sci. 2014; 66C: 1–9. PubMed Abstract | Publisher Full Text

[28] 28. Kazmi F, Hensley T, Pope C, et al.: Lysosomal sequestration (trapping) of lipophilic amine (cationic amphiphilic) drugs in immortalized human hepatocytes (Fa2N-4 cells). Drug Metab Dispos. 2013; 41(4): 897–905. PubMed Abstract | Publisher Full Text | Free Full Text

[29] 29. Nadanaciva S, Lu S, Gebhard DF, et al.: A high content screening assay for identifying lysosomotropic compounds. Toxicol In Vitro. 2011; 25(3): 715–23. PubMed Abstract | Publisher Full Text

[30] 30. de Wilde AH, Jochmans D, Posthuma CC, et al.: Screening of an FDA-approved compound library identifies four small-molecule inhibitors of Middle East respiratory syndrome coronavirus replication in cell culture. Antimicrob Agents Chemother. 2014; 58(8): 4875–84. PubMed Abstract | Publisher Full Text | Free Full Text

[31] 31. Vincent MJ, Bergeron E, Benjannet S, et al.: Chloroquine is a potent inhibitor of SARS coronavirus infection and spread. Virol J. 2005; 2: 69. PubMed Abstract | Publisher Full Text | Free Full Text

[32] 32. Ekins S, Freundlich JS: Computational models for tuberculosis drug discovery. Methods Mol Biol. 2013; 993: 245–62. PubMed Abstract | Publisher Full Text

[33] 33. Ekins S, Freundlich JS, Choi I, et al.: Computational databases, pathway and cheminformatics tools for tuberculosis drug discovery. Trends Microbiol. 2011; 19(2): 65–74. PubMed Abstract | Publisher Full Text | Free Full Text

[34] 34. Ekins S, Williams AJ, Krasowski MD, et al.: In silico repositioning of approved drugs for rare and neglected diseases. Drug Disc Today. 2011; 16(7–8): 298–310. PubMed Abstract | Publisher Full Text

[35] 35. Ekins S, Freundlich JS, Coffee M: Dataset 1 in: A common feature pharmacophore for FDA-approved drugs inhibiting the Ebola virus. F1000Research. 2014. Data Source

A common feature pharmacophore for FDA-approved drugs inhibiting the Ebola virus

Abstract

Keywords

Introduction

Methods

Common features pharmacophore for EBOV actives

Receptor-ligand pharmacophores for VP35

in silico docking of molecules in VP35 structure

Results

Common features pharmacophore for EBOV actives

Figure 1. Pharmacophore based on 4 hits.

Receptor-ligand pharmacophores for VP35

Table 1. Pharmacophores for EBOV VP35 generated from crystal structures in the protein data bank PDB.

in silico docking of molecules in VP35 structure

Figure 2. Docking FDA approved compounds in VP35 protein showing overlap with ligand (yellow) and 2D interaction diagram.

Discussion

Conclusions

Data availability

Author contributions

Competing interests

Grant information

Acknowledgments

Supplementary materials

Supplemental Table 1. Information for common features pharmacophore generation.

Figure S1.Redocking VPL57 in 4IBI.

References

Comments on this article Comments (0)

Open Peer Review

Comments on this article Comments (0)

Open Peer Review

Reviewer Status

Reviewer Reports

Comments on this article

Browse by related subjects

The problem

How to fix it

Competing Interests Policy

Stay Updated