International Journal of Medicinal Chemistry
Volume 2016 (2016), Article ID 2741038, 13 pages
http://dx.doi.org/10.1155/2016/2741038

Research Article

In Silico Designing and Analysis of Inhibitors against Target Protein Identified through Host-Pathogen Protein Interactions in Malaria

Monika Samant,¹ Nidhi Chadha,² Anjani K. Tiwari,² and Yasha Hasija¹

¹Department of Biotechnology, Delhi Technological University, Main Bawana Road, Shahbad Daulatpur, Delhi 110042, India
²Division of Cyclotron and Radiopharmaceutical Sciences, Institute of Nuclear Medicine and Allied Sciences, Brig. S. K. Mazumdar Road, Delhi 110054, India

Received 9 October 2015; Revised 13 November 2015; Accepted 17 November 2015

Academic Editor: Armando Rossello

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

Abstract.

Malaria, a life-threatening blood disease, has been a major concern in the field of healthcare. One of the severe forms of malaria is caused by the parasite Plasmodium falciparum which is initiated through protein interactions of pathogen with the host proteins. It is essential to analyse the protein-protein interactions among the host and pathogen for better understanding of the process and characterizing specific molecular mechanisms involved in pathogen persistence and survival. In this study, a complete protein-protein interaction network of human host and Plasmodium falciparum has been generated by integration of the experimental data and computationally predicting interactions using the interolog method. The interacting proteins were filtered according to their biological significance and functional roles. α-tubulin was identified as a potential protein target and inhibitors were designed against it by modification of amiprophos methyl. Docking and binding affinity analysis showed two modified inhibitors exhibiting better docking scores of −10.5 kcal/mol and −10.43 kcal/mol and an improved binding affinity of −83.80 kJ/mol and −98.16 kJ/mol with the target. These inhibitors can further be tested and validated in vivo for their properties as an antimalarial drug.

1. Introduction

Malaria, one of the most distressing diseases, is caused by the parasitic protozoan Plasmodium falciparum. It takes away millions of lives with the rate increasing each growing year. According to WHO’s Factsheet on the World Malaria Report 2013, 1.2 billion people out of a total of an estimated 3.4 billion are at a high risk of malaria. Malaria is highly prevalent in sub-Saharan Africa where 90% of all malaria deaths occur (WHO 2013). A lot of research has been going on in the field of malarial therapeutics. Nowadays, there are a wide variety of antimalarial drugs, such as chloroquine and artemisinin, and strategies available for the control and treatment of malaria [1–3].

Despite clinical researches in the field of infectious diseases, it remains to be a major problem in the worldwide health issue [4–6]. Exploring the infection process in detail can help in deciphering the mechanisms that govern and control it. In the process of evolution, pathogens have evolved an infection mechanism and humans have evolved immune responses as defense mechanism. A majority of host-pathogen interactions are governed by specific protein-protein interactions [7–9]. To obtain a deep understanding of the infection process, specific interactions between the host and pathogen need to be studied [10–15]. Host-pathogen protein interactions are typically studied using conventional small-scale methods which focuses on single protein at a time. Few methods for large-scale discovery have also been developed such as yeast two-hybrid experiments which allow more comprehensive identification but are expensive and time consuming. Computational methods are less time consuming and cost-effective and hence are a better alternative for the prediction of protein interactions [10, 16–20].

Several studies in the field of computational prediction of host-pathogen interactions and drug discovery have shown significant results. Dyer et al. developed a computational method to identify protein interactions which uses known PPIs to identify functional domains in interacting pairs and Bayesian statistics for assessing the probability of interaction of two proteins. In this study they observed that interacting pairs were coexpressed [10]. Krishnadev and Srinivasan adopted a homology based approach for identification of protein interactions between Plasmodium falciparum and its human host [21]. Lee et al. used an ortholog approach to predict protein interactions in human-Plasmodium system. They compared these interactions with Bayesian and structure based approaches [22]. Durmuş Tekir et al. created Pathogen-Host Interaction Search Tool (PHISTO), a web accessible tool that provides up-to-date information of experimentally verified data on protein interactions. The tool offers integrated visualization of pathogen-host interaction networks, BLAST search, text mining for detecting missing experimental methods, and graph-theoretical analysis of targeted human proteins [23]. Rapanoel et al. adopted interolog method for prediction of protein interactions between Mycobacterium tuberculosis and its human host. This method predicted interactions using experimentally known intraspecies and interspecies interactions and filtered proteins on several parameters, such as cellular location and cellular function, to confirm the practicality of the predictions. Function analysis of the predicted interactions is carried out to analyse the role of proteins in infection process [24]. These studies have contributed to the knowledge of protein interactions through different methods, but a lot is still unknown in the field of malarial therapeutics.

In this study, a complete protein interaction network between human host and Plasmodium falciparum has been developed by integration of experimental and computational methods. Experimental interactions are obtained from PHISTO and STRING. Interolog method is adopted, which hypothesizes that a set of two proteins, each from different species, can be predicted as possible interactions if their respective homologs are found to be interacting in any single species. Interaction network is used as a platform for identification of potential drug target.

Microtubule is a heterodimer consisting of two subunits, that is, α-tubulin and β-tubulin. These subunits bind to each other and make a small subunit which polymerizes to make complete microtubule. It is very important for structural integrity of cells [25, 26]. After the analysis of predicted protein interactions, α-tubulin, a validated target in the malarial infection [27], was found to be one of the highest interacting proteins in malarial infection. The eventual aim of this study is to design an efficient drug molecule for the target. Several tubulin mitotic inhibitors such as benzimidazole, dinitroanilines are already present in the literature which interact with these proteins and hinder the infection process. We have included the study in Figures 10(a) and 10(b). Amiprophos methyl prevents erythrocytic schizogony and blocks mitosis in Plasmodium falciparum infection and results in abnormal microtubule accumulation. This suggests that amiprophos methyl is worthy of investigation for its antimalarial potential. Amiprophos methyl is a validated tubulin inhibitor in reference studies and is found to have least mammalian toxicity [27, 28]. Therefore, derivatives of amiprophos methyl were designed by carrying out modifications at preferred locations with several functional groups. Finally, a molecule is identified which has better binding affinity than the reference molecule which can be considered as a novel drug molecule.

2. Methodology

Whole methodology is divided into two major parts as shown in Figure 1. First part includes identification of target molecule through protein interaction network while the second part involves docking analysis of the target with designed inhibitors.

(a)

(b)

Figure 1: Flowchart of methodology. (a) Identification of the target molecule by filtering and analysis of the protein interactions predicted by BIPS. (b) Docking and binding affinity analysis of the target protein and modified inhibitors.

2.1. Host-Pathogen Interactions

The complete protein interaction network was developed by integration of experimentally found interactions, interaction data fetched from protein interaction database STRING, and computationally predicted interactions.

The computational prediction of host-pathogen protein interactions was carried out through BIPS, that is, BIANA interolog prediction server (http://sbi.imim.es/web/index.php/research/servers/bips) [29]. The whole proteome of Plasmodium falciparum, downloaded from UniProt (http://www.uniprot.org/), was taken as an input by the server and the predicting protein partner was restricted to Homo sapiens. The predicted interactions were annotated for their cellular location, gene ontologies, and functional processes, using the Database for Annotation, Visualization, and Integrated Discovery (DAVID), Functional Annotation Chart tool [30, 31], and Panther [32]. These protein interactions were then functionally annotated for the functional processes and pathways involved.

The experimental protein interactions were obtained from PHISTO, that is, Pathogen-Host Interaction Search Tool (http://www.phisto.org/). The PHISTO server has a drop-down menu where one can select the type of pathogen for which experimentally known interactions with human are desired [23].

STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) is a database for known and predicted protein interactions for a large number of organisms. The complete database was downloaded and data mining was performed in the following two steps:

(1)The interactions with a sum score more than 900 were filtered.(2)Interacting proteins were then filtered on the basis of participating organism (Homo sapiens-Plasmodium falciparum interactions).

Protein interactions from the three sources were integrated to develop a complete human-Plasmodium interaction network. Protein interactions were analysed to find most significant proteins involved in the infection process. Pathogen proteins were then shortlisted on the basis of the number of their interacting partners. Highest interacting pathogen proteins were shortlisted for further analysis. The shortlisted proteins were then studied in literature for their functionality, role in the pathogen, and the infection mechanism. A potential drug target was found on the basis of above-mentioned analysis.

2.2. Ligand Dataset Preparation

Finalized drug target is a known target for cancer; hence it was searched for its present inhibitors in literature. Out of all already known inhibitors, the one with least known toxicity and highest absorption was finalized and its derivatives were designed by addition of specific chemical groups at the preferred positions according to the present literature. Reference inhibitor molecule chosen was amiprophos methyl. A library of the designed inhibitors was prepared by MarvinSketch. Designing of inhibitors was performed by addition of specific functional groups at single and multiple locations in many combinations [33].

2.3. Docking Studies

Rigid docking was carried out between protein target and designed inhibitors using GLIDE docking program of Schrödinger (Schrödinger Release 2014-2). The results of docking were analysed followed by prioritization of functional groups and positions on the chemical structure. More molecules were designed by combining the prioritized position and functional groups. Prime docking molecules were then taken for flexible docking by induced fit docking in Schrödinger (Schrödinger Release 2014-2) (Figure 9). Most suitable docking inhibitor is considered as a potential drug molecule [34–37].

The final molecule is checked for its binding affinity and ADME properties [38] (Table 2).

3. Result and Discussion

The whole proteome of Plasmodium falciparum was downloaded from UniProt which contained 5353 proteins; 157 reviewed and 5196 unreviewed proteins. From BIPS, a total of 2381 interspecies interactions were obtained. These interactions were annotated for their cellular location, gene ontology, and functional role and then filtered upon these parameters. Functional annotation of the predicted pathogen proteins was carried out using DAVID and Panther. Panther identified 31 pathogen proteins and classified them on the basis of biological process, molecular function, pathways, and protein class. It was observed that most of the interacting pathogen proteins are involved in metabolic processes and have binding as their molecular function. These proteins fall into nucleic acid binding class of proteins and thus it was concluded that they have an important function in DNA replication and cell survival. There was no majority observed in the pathway classification of these proteins but there was a similarity: they all were involved in signalling and disease pathways (Figure 2).

(a)

(b)

(c)

(d)

Figure 2: Pie chart representation of pathogen protein classified on the basis of (a) biological process, (b) molecular function, (c) protein class, and (d) pathway.

The predicted protein interactions were used to generate host-pathogen interaction network (Figure 3). The predicted interactions were then analysed for the number of interacting partners of pathogen. Out of all the protein interactions, highly interacting proteins were shortlisted. Human proteins were not considered as a drug target because those proteins might have a role as an essential component in the biological processes. So, targeting human proteins have a high risk of drug toxicity. Human proteins which were interacting with more than three pathogen proteins were usually structural and assembly proteins such as actin, tubulin, and histone. Most of these human proteins were histone proteins. This shows that parasite infection affects mostly nuclear and cell assembly proteins in humans. Therefore, these proteins were left out because they are important proteins in human cells. If human proteins or proteins similar to human proteins are targeted, the problem of drug toxicity after drug delivery might persist. Hence, pathogen proteins are considered for drug targeting to remove the possibility of negative effects in the host [39]. It was observed that most of the highly interacting proteins were involved in structural assembly of the pathogen such as actin, tubulin, and histone. Already present literatures of these proteins were studied to understand their functional role in the pathogen.

Figure 3: Protein-protein interaction network of human host and Plasmodium falciparum. Different colours and shapes are representing proteins of different species.

As a result of this analysis, α-tubulin was finalized as an important protein involved in the infection process. α-tubulin is a monomer unit which polymerizes to carry out several critically important roles throughout parasite life cycle [40, 41]. In parasite, they form mitotic spindle during cell division and even slight disruption of microtubule causes a severe impact on the viability of parasite. Plasmodium falciparum infects the host and an initial contact occurs between merozoite and erythrocyte. α-tubulin is also present in humans but tubulins of humans and plasmodium have significant differences in their mode of function. Present studies have shown that antitubulins show very low or no cross-reactivity to mammalian tubulins [28]. In experimental studies it was confirmed that invasion was decreased and eventually stopped when merozoites were exposed to tubulin inhibitors [27]. Experimental studies have demonstrated that microtubules were disrupted on exposure to antitubulin agent indicating the role of intact microtubule in merozoite invasion. Microtubule is found in many stages of malaria parasite, validating it as a potential drug target. Detailed examination of merozoites in erythrocyte invasion identified targeting α-tubulin as a potential approach for malaria therapy [42–44].

Amiprophos methyl (Figure 4) is an antimitotic herbicide and a known inhibitor for α-tubulin. It is found to be a promising molecule because of its low mammalian toxicity. It was reported in reference studies that amiprophos methyl has better specificity for pathogen proteins and has no binding site in human tubulin protein [29, 30, 44–46].

Figure 4: Molecular structure of amiprophos methyl.

Glide docking analysis of the target protein with the designed inhibitors and reference molecule was carried out. Glide docking performs significantly better than other docking programs with accuracy of more than 90% [34–37]. Amiprophos methyl was considered as a reference molecule with a docking score of 4.43. Among the designed inhibitors, molecules with best docking score were analysed for the functional molecules and their positions in that molecule. Molecules were then designed with the functional groups in combinations at the preferred positions. It was found that electronegative groups at positions 2, 3, and 9 showed significant improvements in docking score. A molecule with CF₃ at position 4, piperidine at position 9, and OH at position 3 showed the best docking score of −8.14 at site 2. This molecule had 83% better docking score than the reference molecule. The molecules with the best docking score were shortlisted and flexible docking was carried out (Figure 5). Highest docking score of −10.5 kcal/mol and −10.43 kcal/mol was obtained in flexible docking (Table 3).

(a)

(b)

Figure 5: Molecular structure of final shortlisted molecules.

These final molecules are analysed for their interaction with residues at site II using LigPlot [46]. 3-(Ethoxy[(piperidin-1-yl)amino]phosphoryl} oxy)-2-(hydroxymethyl)-6-(trifluoromethyl)phenol forms two hydrogen bonds with Glu(22) and Tyr(83). 5-(Ethoxy[(piperidin-1-yl)amino]phosphoryl} oxy)-4-(hydroxymethyl)-2-(trifluoromethyl) benzene-1,3-diol forms four hydrogen bonds with surrounding residue, that is, two with Arg(229), two with Thr(82), and one with Glu(77). The reference molecule was forming only a single hydrogen bond with Trp(21). The hydrogen bonds of these molecules with surrounding residues govern their stability and hence new molecules have better stability than reference molecule (Figures 6, 7, and 8).

(a)

(b)

(c)

Figure 6: Interaction of amiprophos methyl with target protein at site II (a) and (b) docked reference molecule in the target protein at the site II with blue colour showing hydrogen bind between ligand and Tyr(21) of protein. (c) LigPlot analysis of protein and ligand interactions.

(a)

(b)

(c)

Figure 7: Interaction of 3-(ethoxy[(piperidin-1-yl)amino]phosphoryl} oxy)-2-(hydroxymethyl)-6-(trifluoromethyl) phenol with target protein at site II (a) and (b) docked molecule in the target protein at the site II with different colours showing different polarity and charges. Pink, blue, green, and cyan represent negative charge and positive charge and hydrophobic and polar residues, respectively. (c) LigPlot analysis of protein and ligand interactions.

(a)

(b)

(c)

Figure 8: Interaction of 5-(ethoxy[(piperidin-1-yl)amino]phosphoryl} oxy)-4-(hydroxymethyl)-2-(trifluoromethyl) benzene-1,3-diol with target protein at site II (a) and (b) docked molecule in the target protein at the site II with different colours showing different polarity and charges. Pink, blue, green, and cyan represent negative charge, positive charge and hydrophobic and polar residues, respectively. (c) LigPlot analysis of protein and ligand interactions.

Figure 9: Molecular structure of the molecules shortlisted for flexible docking and binding affinity analysis.

(a)

(b)

Figure 10: (a) Chemical structure of benzimidazole. (b) Chemical structure of dinitroanilines.

The two final listed molecules have the highest docking score and low binding energy than the reference molecule which shows that it has better binding properties (Table 1). From ADME property analysis of the two molecules, it was observed that these have high values of human oral absorption, zero violations in Lipinski’s rule of five, and zero violation in Jorgensen’s rule of three and all the drug-ability determining properties lie under range (Table 2).

Table 1: Comparison of binding energy scores of reference and potential drug molecules.

Mol. number
H-bonding energy
Coulomb energy
Covalent binding energy
Pi-pi packing energy
Lipophilic energy
Binding energy
Electrosolvation energy
Van der Waals energy

Ref.
−0.33
3.52
4.77
−0.17
−42.46
−80.96
5.05
−51.35

A
−2.35
−19.23
2.09
−3.63
−37.46
−83.79
21.08
−44.34

B
−3.78
−66.20
5.94
0
−32.50
−98.16
37.69
−39.30

Table 2: Comparison between QikProp scores of final molecules and reference molecule for analysis of ADME properties.

Mol. number
# stars
MW
Human oral absorption
Percent human oral absorption
Rule Of five
Rule Of three

Ref. mol.
0
304.3
3
100
0
0

A
0
398.32
3
80.59
0
0

B
0
414.32
3
69.87
0
0

Table 3: Glide scores and induced fit scores of the short listed molecules for flexible docking and binding affinity analysis.

Mol. number
Glide score
IFD score

1
−10.5
−762.72

2
−10.43
−762.29

3
−10.24
−766.93

4
−9.65
−768.15

5
−9.25
−761.87

6
−8.99
−760.97

7
−8.46
−764.92

8
−7.61
−767.01

9
−6.68
−757.47

10
−6.65
−763.28

Better binding affinity and ADME property analysis of the final molecule confirm its potential to be used as a drug molecule for further analysis and trials.

4. Conclusion

From the predicted host-pathogen PPIs, the present study concludes that most of the host proteins with which pathogen protein interactions are structural proteins such as actin, tubulin, and histone. Most of the pathogen proteins involved in the infection process are structural and assembly proteins and most of the host proteins are either structural proteins or nuclear assembly proteins. Hence, the pathogen caused infection by targeting nuclear assembly proteins and thereby inhibiting the host cell to function properly. α-tubulin of pathogen is targeted for development of antimalarial agent for malarial treatment. Derivatives of herbicide having antimalarial property were developed and molecule with better binding affinity and ADME property was obtained. It was observed that molecules with electronegative groups have better binding properties than original molecule.

3-(Ethoxy[(piperidin-1-yl)amino]phosphoryl} oxy)-2-(hydroxymethyl)-6-(trifluoromethyl)phenol and 5-(ethoxy[(piperidin-1-yl)amino]phosphoryl} oxy)-4-(hydroxymethyl)-2-(trifluoromethyl)benzene-1,3-diol were the two best molecules which can be considered as drug molecules for in vivo analysis and validation.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The authors would like to thank Kirti Bhadhadhara for scientific discussions.

References

M. H. Gelb, “Drug discovery for malaria: a very challenging and timely endeavor,” Current Opinion in Chemical Biology, vol. 11, no. 4, pp. 440–445, 2007.

T. N. C. Wells, “Discovering and developing new medicines for malaria control and elimination,” Infectious Disorders—Drug Targets, vol. 13, no. 4, pp. 292–302, 2013.

K. G. Le Roch, D.-W. D. Chung, and N. Ponts, “Genomics and integrated systems biology in Plasmodium falciparum: a path to malaria control and eradication,” Parasite Immunology, vol. 34, no. 2-3, pp. 50–60, 2012.

R. W. Snow, C. A. Guerra, A. M. Noor, H. Y. Myint, and S. I. Hay, “The global distribution of clinical episodes of Plasmodium falciparum malaria,” Nature, vol. 434, no. 7030, pp. 214–217, 2005.

W. P. O'Meara, J. N. Mangeni, R. Steketee, and B. Greenwood, “Changes in the burden of malaria in sub-Saharan Africa,” The Lancet Infectious Diseases, vol. 10, no. 8, pp. 545–555, 2010.

S. Foster and M. Phillips, “Economics and its contribution to the fight against malaria,” Annals of Tropical Medicine and Parasitology, vol. 92, no. 4, pp. 391–398, 1998.

K. I. Barnes, D. N. Durrheim, F. Little et al., “Effect of artemether-lumefantrine policy and improved vector control on malaria burden in KwaZulu-Natal, South Africa,” PLoS Medicine, vol. 2, no. 11, article e330, 2005.

A. Casadevall and L.-A. Pirofski, “Host-pathogen interactions: redefining the basic concepts of virulence and pathogenicity,” Infection and Immunity, vol. 67, no. 8, pp. 3703–3713, 1999.

A. Casadevall and L.-A. Pirofski, “Host-pathogen interactions: basic concepts of microbial commensalism, colonization, infection, and disease,” Infection and Immunity, vol. 68, no. 12, pp. 6511–6518, 2000.

M. D. Dyer, T. M. Murali, and B. W. Sobral, “Computational prediction of host-pathogen protein-protein interactions,” Bioinformatics, vol. 23, no. 13, pp. i159–i166, 2007.

J. Mufunda, P. Nyarango, A. Usman et al., “Roll back malaria—an African success story in Eritrea,” South African Medical Journal, vol. 97, no. 1, pp. 46–50, 2007.

E. A. Franzosa, S. Garamszegi, and Y. Xia, “Toward a three-dimensional view of protein networks between species,” Frontiers in Microbiology, vol. 3, article 428, 2012.

J. Zahiri, J. H. Bozorgmehr, and A. Masoudi-Nejad, “Computational prediction of protein-protein interaction networks: algo-rithms and resources,” Current Genomics, vol. 14, no. 6, pp. 397–414, 2013.

N. Yosef, M. Kupiec, E. Ruppin, and R. Sharan, “A complex-centric view of protein network evolution,” Nucleic Acids Research, vol. 37, no. 12, article e88, 2009.

P. Uetz, L. Glot, G. Cagney et al., “A comprehensive analysis of protein-protein interactions in Saccharomyces cerevisiae,” Nature, vol. 403, no. 6770, pp. 623–627, 2000.

J. Mendez-Rios and P. Uetz, “Global approaches to study protein-protein interactions among viruses and hosts,” Future Microbiology, vol. 5, no. 2, pp. 289–301, 2010.

P. V. Desai, A. Patny, J. Gut et al., “Identification of novel parasitic cysteine protease inhibitors by use of virtual screening. 2. The available chemical directory,” Journal of Medicinal Chemistry, vol. 49, no. 5, pp. 1576–1584, 2006.

L. R. Matthews, P. Vaglio, J. Reboul et al., “Identification of potential interaction networks using sequence-based searches for conserved protein-protein interactions or ‘interologs’,” Genome Research, vol. 11, no. 12, pp. 2120–2126, 2001.

Y. Qi, O. Tastan, J. G. Carbonell, J. Klein-Seetharaman, and J. Weston, “Semi-supervised multi-task learning for predicting interactions between HIV-1 and human proteins,” Bioinformatics, vol. 26, no. 18, Article ID btq394, pp. i645–i652, 2010.

S. Wuchty, “Computational prediction of host-parasite protein interactions between P. falciparum and H. sapiens,” PLoS ONE, vol. 6, no. 11, Article ID e26960, 2011.

O. Krishnadev and N. Srinivasan, “A data integration approach to predict host-pathogen protein-protein interactions: application to recognize protein interactions between human and a malarial parasite,” In Silico Biology, vol. 8, no. 3-4, pp. 235–250, 2008.

S.-A. Lee, C.-H. Chan, C.-H. Tsai et al., “Ortholog-based protein-protein interaction prediction and its application to inter-species interactions,” BMC Bioinformatics, vol. 9, supplement 12, article S11, 2008.

S. Durmuş Tekir, T. Çakir, E. Ardiç et al., “PHISTO: pathogen-host interaction search tool,” Bioinformatics, vol. 29, no. 10, pp. 1357–1358, 2013.

H. A. Rapanoel, G. K. Mazandu, and N. J. Mulder, “Predicting and analyzing interactions between Mycobacterium tuberculosis and its human host,” PLoS ONE, vol. 8, no. 7, Article ID e67472, 2013.

J.-C. Barale, T. Blisnick, H. Fujioka et al., “Plasmodium falciparum subtilisin-like protease 2, a merozoite candidate for the merozoite surface protein 1–42 maturase,” Proceedings of the National Academy of Sciences of the United States of America, vol. 96, no. 11, pp. 6445–6450, 1999.

G. S. Dow, A. Armson, M. R. Boddy et al., “Plasmodium: assessment of the antimalarial potential of trifluralin and related compounds using a rat model of malaria, Rattus norvegicus,” Experimental Parasitology, vol. 100, no. 3, pp. 155–160, 2002.

B. J. Fennell, J. A. Naughton, E. Dempsey, and A. Bell, “Cellular and molecular actions of dinitroaniline and phosphorothioamidate herbicides on Plasmodium falciparum: tubulin as a specific antimalarial target,” Molecular and Biochemical Parasitology, vol. 145, no. 2, pp. 226–238, 2006.

C. Mara, E. Dempsey, A. Bell, and J. W. Barlow, “Synthesis and evaluation of phosphoramidate and phosphorothioamidate analogues of amiprophos methyl as potential antimalarial agents,” Bioorganic and Medicinal Chemistry Letters, vol. 21, no. 20, pp. 6180–6183, 2011.

J. Garcia-Garcia, S. Schleker, J. Klein-Seetharaman, and B. Oliva, “BIPS: BIANA Interolog Prediction Server. A tool for protein-protein interaction inference,” Nucleic Acids Research, vol. 40, no. 1, pp. W147–W151, 2012.

D. W. Huang, B. T. Sherman, and R. A. Lempicki, “Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists,” Nucleic Acids Research, vol. 37, no. 1, pp. 1–13, 2009.

D. W. Huang, B. T. Sherman, and R. A. Lempicki, “Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources,” Nature Protocols, vol. 4, no. 1, pp. 44–57, 2009.

H. Mi, A. Muruganujan, and P. D. Thomas, “PANTHER in 2013: modeling the evolution of gene function, and other gene attributes, in the context of phylogenetic trees,” Nucleic Acids Research, vol. 41, no. 1, pp. D377–D386, 2013.

G. Madhavi Sastry, M. Adzhigirey, T. Day, R. Annabhimoju, and W. Sherman, “Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments,” Journal of Computer-Aided Molecular Design, vol. 27, no. 3, pp. 221–234, 2013.

R. A. Friesner, J. L. Banks, R. B. Murphy et al., “Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy,” Journal of Medicinal Chemistry, vol. 47, no. 7, pp. 1739–1749, 2004.

R. A. Friesner, R. B. Murphy, M. P. Repasky et al., “Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes,” Journal of Medicinal Chemistry, vol. 49, no. 21, pp. 6177–6196, 2006.

D. A. Gschwend, A. C. Good, and I. D. Kuntz, “Molecular docking towards drug discovery,” Journal of Molecular Recognition, vol. 9, no. 2, pp. 175–186, 1996.

T. A. Halgren, R. B. Murphy, R. A. Friesner et al., “Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening,” The Journal of Medicinal Chemistry, vol. 47, no. 7, pp. 1750–1759, 2004.

J. C. Shelley, A. Cholleti, L. L. Frye, J. R. Greenwood, M. R. Timlin, and M. Uchimaya, “Epik: a software program for pKa prediction and protonation state generation for drug-like molecules,” Journal of Computer-Aided Molecular Design, vol. 21, no. 12, pp. 681–691, 2007.

H. N. Chua, W.-K. Sung, and L. Wong, “Exploiting indirect neighbours and topological weight to predict protein function from protein-protein interactions,” Bioinformatics, vol. 22, no. 13, pp. 1623–1630, 2006.

D. J. Rawlings, H. Fujioka, M. Fried, D. B. Keister, M. Aikawa, and D. C. Kaslow, “α-Tubulin II is a male-specific protein in Plasmodium falciparum,” Molecular and Biochemical Parasitology, vol. 56, no. 2, pp. 239–250, 1992.

N. Gerald, B. Mahajan, and S. Kumar, “Mitosis in the human malaria parasite Plasmodium falciparum,” Eukaryotic Cell, vol. 10, no. 4, pp. 474–482, 2011.

T. W. A. Kooij, C. J. Janse, and A. P. Waters, “Plasmodium post-genomics: better the bug you know?” Nature Reviews Microbiology, vol. 4, no. 5, pp. 344–357, 2006.

C. Mara, E. Dempsey, A. Bell, and J. W. Barlow, “Synthesis and evaluation of phenoxyoxazaphospholidine, phenoxyoxazaphosphinane, and benzodioxaphosphininamine sulfides and related compounds as potential anti-malarial agents,” Bioorganic and Medicinal Chemistry Letters, vol. 23, no. 12, pp. 3580–3583, 2013.

E. Lacey, “Mode of action of benzimidazoles,” Parasitology Today, vol. 6, no. 4, pp. 112–115, 1990.

L. M. MacDonald, A. Armson, R. C. A. Thompson, and J. A. Reynoldson, “Characterisation of benzimidazole binding with recombinant tubulin from Giardia duodenalis, Encephalitozoon intestinalis, and Cryptosporidium parvum,” Molecular and Biochemical Parasitology, vol. 138, no. 1, pp. 89–96, 2004.

A. C. Wallace, R. A. Laskowski, and J. M. Thornton, “LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions,” Protein Engineering, vol. 8, no. 2, pp. 127–134, 1995.