Analysing the essential proteins set of Plasmodium falciparum PF3D7 for novel drug targets identification against malaria
Malaria Journal volume 20, Article number: 335 (2021)
Plasmodium falciparum is an obligate intracellular parasite of humans that causes malaria. Falciparum malaria is a major public health threat to human life responsible for high mortality. Currently, the risk of multi-drug resistance of P. falciparum is rapidly increasing. There is a need to address new anti-malarial therapeutics strategies to combat the drug-resistance threat.
The P. falciparum essential proteins were retrieved from the recently published studies. These proteins were initially scanned against human host and its gut microbiome proteome sets by comparative proteomics analyses. The human host non-homologs essential proteins of P. falciparum were additionally analysed for druggability potential via in silico methods to possibly identify novel therapeutic targets. Finally, the PfAp4AH target was prioritized for pharmacophore modelling based virtual screening and molecular docking analyses to identify potent inhibitors from drug-like compounds databases.
The analyses identified six P. falciparum essential and human host non-homolog proteins that follow the key druggability features. These druggable targets have not been catalogued so far in the Drugbank repository. These prioritized proteins seem novel and promising drug targets against P. falciparum due to their key protein–protein interactions features in pathogen-specific biological pathways and to hold appropriate drug-like molecule binding pockets. The pharmacophore features based virtual screening of Pharmit resource predicted a lead compound i.e. MolPort-045–917-542 as a promising inhibitor of PfAp4AH among prioritized targets.
The prioritized protein targets may worthy to test in malarial drug discovery programme to overcome the anti-malarial resistance issues. The in-vitro and in-vivo studies might be promising for additional validation of these prioritized lists of drug targets against malaria.
Malaria is a life-threatening infectious disease caused by parasitic protozoan plasmodium. It is a vector borne disease, transmitted to humans through the bite of infected carrier female Anopheles mosquitoes. Among five parasite species that cause malaria in humans, two species, Plasmodium falciparum and Plasmodium vivax, have the greatest threat to human life . The P. falciparum, a unicellular protozoan, belongs to the family Plasmodiidae and lies in the phylum Apicomplexa . Plasmodium falciparum alone is responsible for almost all malaria-inflicted deaths in sub-Saharan Africa, with the continent bearing over 90% of the global P. falciparum burden [3,4,5]. Asia is second to Africa in terms of malaria prevalence. In 2019, the World Health Organization (WHO) estimated 229 million malaria cases and about 409,000 deaths due to malaria worldwide every year . More than 85% of confirmed recorded cases and deaths in Asia occurred in India, Indonesia, Myanmar, and Pakistan .
The P. falciparum resistance is reported to many approved anti-malarial drugs, including the chloroquine and artemisinin . Resistance to chloroquine was first observed in Thailand in 1957 and the Colombian-Venezuelan border in 1959 . Drug resistance had spread across sub-Saharan Africa by 1988, and today chloroquine is no longer effective in almost all parts of the world . Specific polymorphisms in the P. falciparum chloroquine resistance transporter (PfCRT) are reported in association with chloroquine resistance . Likewise, the artemisinin resistance was first documented in 2008 in the Thailand-Cambodia border regions [12, 13]. The artemisinin resistance associated with delayed parasitic clearance after three days of artemisinin monotherapy. Several studies have reported that artemisinin resistance emerged due to a polymorphism in the pfk13 gene [14, 15].
The indispensable proteins of P. falciparum recently explored from well experimental approaches in some studies. These repositories are promising to identify suitable targets to overcome the drug-resistant P. falciparum infection [16, 17]. Zhang et al.  experimentally analysed 5399 genes and identified 2680 as essential for optimal growth of P. falciparum during asexual blood stages. These essential genes coding the P. falciparum vital proteins, including drug targets and potential vaccine candidates. Besides, there are over 1,000 of Plasmodium-conserved essential genes with unknown biological functions so far. In the current study, these essential proteins were retrieved from the two recently published studies and assessed for druggable potential based on comparative proteomics, protein–protein interactions and drug-like molecules binding potential.
The methodological layout of the current study is depicted in Fig. 1.
Retrieval of non-paralogous essential proteins
The essential proteins of P. falciparum strain 3D7 were retrieved from two recently published studies [16, 17]. The paralogous proteins were removed by CD-HIT clustering analysis with ≥ 80% sequence similarity cutoff .
Human host non-homologous and virulent proteins identification
The P. falciparum essential proteins non-homolog to human host proteome were identified by comparative sequence analyses via BLASTp tool . The threshold values of ≤ 35% query coverage and sequence identity were set during this analysis . The proteins having significant similarity with human proteome were discarded and the remaining non-homologs were shortlisted for further analysis. The non-homology search against human gut microbiota proteins sequences was also carried out with a threshold cutoff, i.e. E value ≥ 0.001 [21, 22]. The Victors database was screened for P. falciparum 3D7 virulent proteins annotation. The Victors database contains 5304 virulent proteins data from various parasites including P. falciparum .
Drugbank database scanning
The shortlisted essential proteins of P. falciparum from above analyses were scanned against the Drugbank database to identify novel drug targets with ≤ 60% query coverage and percent identity threshold of BLASTp .
Structure homologs search
The proteins data bank (PDB) was screened to identify the homologous 3D structures of pathogenic proteins . The pathogen proteins were BLAST against the entire PDB database entries [22, 26]. The pathogen sequences having ≥ 60% homology were modeled with Swiss Model  and verified by ERRAT  and RAMPAGE .
The prioritized list of essential proteins shortlisted from above analyses were tested for druggability potential. The drug-like molecules binding pockets of the targets were identified by PockDrug-server . The subcellular localization was performed with CELLO v.2.5 . The molecular weight of the proteins were also accounted during their druggability assessment . The PPI analysis was performed by the STRING database and the Hub proteins were identified based on node degree (K ≥ 5) that represents the number of interactions .
Pharmacophore based virtual screening
Pharmit server was employed to design a pharmacophore model using the 3D structure of P. falciparum diadenosine tetraphosphate hydrolase (PfAp4AH) enzyme (PDB ID 5CFJ) . Pharmit provides interactive screening of millions of chemical compounds from built-in databases, i.e. Molprot, ChEMBL, ZINC, and PubChem. The Pharmit server is based on a pharmacophore model using AutoDock Vina scoring function . The pharmacophore model was built using seven features, i.e. two hydrogen donors, two hydrogen acceptors, two hydrophobic features, and one aromatic feature. The screening results were then minimized to a significant level based on Pharmit scoring and RMSD values to obtain the best possible inhibitors out of millions of drug-like compounds. The top ten hits identified based on Pharmit screening were then docked using the CB-Dock tool to verify the ligand-binding sites . Discovery studio was used for protein–ligand interaction visualization .
Subtractive proteomic analyses
Essential genes perform key cellular functions for the survival of pathogens . The P. falciparum strain 3D7 essential genes information were obtained from the recently published articles [16, 17] and total 3380 essential proteins were identified. These proteins sequences were retrieved from Uniprot. Seven paralogous protein sequences were excluded by CD-HIT analysis and the remaining 3373 sequences were considered for downstream analysis (Additional file 1: Table S1). The non-paralogous protein sequences were subjected to BLASTp against human proteome as well as human gut microbiome proteome data with threshold parameters. The gut flora is helpful to the host in many ways like vitamins biosynthesis and absorption of short chain-fatty acids . The unintended inhibition of gut microbe leads to a decrease in gut flora and colonization of pathogenic bacteria in the host gut . The comparative sequence analysis based screening (see, methodology section) against human proteome identified total 2852 P. falciparum proteins non-homolog to human host. Further screening of these P. falciparum essential and human host non-homolog proteins against human gut microbiome proteins database identified 2669 homologs and 183 non-homologs proteins (Fig. 1; Additional file 2: Table S2; Additional file 3: Table S3). The five among these 183 were annotated as P. falciparum 3D7 virulent proteins during screening of Victors database (Table 1). The 183 P. falciparum essential proteins identified as human as well as human gut proteome non-homologs were prioritized for downstream analyses. The homology screening of these 183 proteins based on set threshold against Drugbank repository inferred no homology with already reported drug targets deposited in drugbank database. However, 38 established anti-malarial targets from Alexander et al.  (Additional file 4: Table S4) were additionally catalogued to check the recovery of previously reported targets in current study. Among approved anti-malarial targets, the P. falciparum 2-C-methyl-D-erythritol 4-phosphate cytidylyltransferase and P. falciparum subtilisin-like protease-1 were listed within the 183 prioritized targets (Additional file 2: Table S2). However, when the cutoff criteria was relaxed to 35% sequence identity and 35% query coverage to screen the human gut microbiome database, then additional 16 already approved anti-malarial targets are recovered (Additional file 4: Table S4). The recovery of these established anti-malarial targets somehow validates the strategy been acquired in current study.
The Drugbank non-homologous proteins were prioritized for downstream druggability analyses. The subcellular localization is one of the key aspect of druggability and the cytoplasmic proteins are considered as suitable drug targets [40, 41]. The 15 P. falciparum proteins among shortlisted prioritized targets were annotated as cytoplasmic proteins. The proteins 3D structures identified by Swiss model were validated with the ERRAT tool with quality factor score of > 50, which is accepted as high quality model . Ramachandran plot identified 80–90% of modelled proteins residues in the allowed region assuring good quality structure modelling of the target proteins (Table 2). Finally, six (06) proteins were prioritized on the basis of (i) pockdrug probability score ˃0.5 , (ii) ERRAT quality factor ≥ 90 , and (iii) protein–protein interaction node degree i.e. K ≥ 5  (Table 2; Fig. 2). These six prioritized targets are speculating to hold promising druggable pockets to anchor small drug-like molecules and act as indispensable hub proteins in P. falciparum metabolic network.
Pharmacophore based virtual screening
Among finally shortlisted targets, the PfAp4AH was prioritized for drug-like compounds virtual screening. The pharmacophore model designed from the three-dimensional structure of ligand binding pocket of PfAp4AH showed seven features, i.e. two hydrogen donors, two hydrogen acceptors, two hydrophobic features, and one aromatic feature shown in Fig. 3A. Top 10 hits acquired from Pharmit resource screening were prioritized based on score and RMSD values (Table 3). These top hits compounds were then docked within the ligand binding pocket of PfAp4AH enzyme to predict their binding conformation. The 3D structure of the PfAp4AH enzyme reported with bound substrate (i.e. ATP) within active-site. The site comprises of seven key residues, i.e. Tyr87, Lys94 Ser135, Pro133, His43, Lys48, and Glu115, responsible for ATP hydrolysis . This active-site was utilized in the current study for pharmacophore-based virtual screening. The three key residues of this active-site, i.e. Tyr87, Glu115, and Pro133 are found to primarily engage with the binding of top hit compound (i.e. MolPort-045–917-542) (Fig. 3B). The conserved Glu115 residue is reported as especially important for hydrolysis reaction of PfAp4AH . Besides, the docking analyses predicted that all the top screened compounds explicitly anchor to this active site.
In the current study, the 3380 essentially reported proteins of P. falciparum strain 3D7 were analysed to address potent novel druggable targets. These proteins were analysed based on their non-homology with the human host as well as human gut microbiome proteome. The targets were additionally shortlisted based on strict threshold criteria of basic druggability features. Among the shortlisted targets (Table 2), the protein, Bis(5'-nucleosyl)-tetraphosphatase (asymmetrical) (C0H4F3_PLAF7) also known as asymmetrical diadenosine 5′,5″-P1,P4-tetraphosphate hydrolase (PfAp4A) (EC 184.108.40.206) enzyme participates in pyrimidine and purine metabolism [44, 45]. The PfAp4A hydrolase exhibited high-temperature stability even at 60 °C . Previously in few studies, the PfAp4A is also tested as a potential drug target against P. falciparum [47, 48].
The protein serine/threonine-protein phosphatase (Q8I2N2_PLAF7) was also found among finally shortlisted target that involve in regulation of many cellular signaling pathways by catalysing the removal of phosphate group from target enzymes. This enzyme plays a central role in the functional regulation and control of different genes related to the cell cycle . The phosphorylation regulates several primary steps in P. falciparum's diverse life cycles. Many of the kinases and phosphatases as well as their substrates are specific to parasites, making eventually the phosphorylation event as a viable target for anti-parasitic action . The protein phosphatase-1, a type of PfPP, involve in the mitotic division of P. falciparum and plays an important role in the liberation of merozoites. Prior studies on P. falciparum revealed that the activity of PfPP1 is more important as compared to protein phosphatase 2A (PP2A) . This also verified by transcriptomic analysis, where the PfPP1 transcript levels reported higher than PP2A after 24 h of RBC infection . The okadaic acid (OA), a toxin initially isolated from a marine sponge, i.e. Halichondria okadai has been identified as a selective inhibitor of serine/threonine protein phosphatases (PPPs) and reported to strongly inhibits the PP1, 2A, and 2B in-vitro . Out of 30 examined protein phosphatase, the 16 protein phosphatases along with PP1 and putative phosphatases seem to be important for blood-stage parasites . Moreover, some studies also showed that PfPP1 is indispensable for blood-stage parasite survival . Many phosphatases play key roles in the pathological pathways, and their inactivation may help to prevent or postpone the emergence of human diseases. Therefore, the potent inhibitors for such phosphatases might be of great therapeutic benefit.
The enzyme cytochrome b5 Reductase (cb5r) (Q8I599_PLAF7) plays a role in fatty acid elongation, cholesterol biosynthesis, and cytochrome P450-mediated detoxification of xenobiotics . This protein has been thoroughly studied in mammals, but still needs to be characterized in microorganisms, such as fungi and parasites, including P. falciparum. There is a close phylogenetic relationship between the plant and P. falciparum cb5r proteins. The plant cb5r has been identified as a novel herbicidal target . This protein reported essential for P. falciparum survival and was found human host non-homolog and possibly be a potent therapeutic target, thereby might be a worthy candidate for drug development against malaria.
The vacuolar protein sorting-associated protein 29 (VPS29) (Q8IM27_PLAF7) is involved in the essential metabolic process of proteins translocation to the subcellular organelles. The P. falciparum sort and traffic newly synthesized proteins to target intracellular organelles as well as beyond the plasma membrane into the host cell in some cases . The P. falciparum VPS29 (i.e. PfVPS29) is the functional component in the assembly of the retromer complex . During the PPI analysis, the pfVPS29 showed direct interactions with other retromer complex components i.e. PfVPS26, VPS9, VPS10 as shown in Fig. 2. The PfVPS29 is located in the cytosol and highly expressed in early trophozoite and schizont stages . Inhibiting the activity of PfVPS29 may lead to the disassembling of the retromer complex and possibly halt the protein sorting function of the P. falciparum.
The multifunctional methyltransferase subunit (Q8IM19_PLAF7) have methyltransferase activity during post-translational modifications, chromatin remodeling and protein heterodimerization activity . The protein methyltransferases (PMTs) have been linked to the pathogenesis of a variety of diseases, including human cancers, inflammatory diseases, metabolic diseases, and neurodegenerative diseases. The PMTs are highly attractive among the histone-modifying enzymes and act as drug targets [59, 60]. However, to date no study been conducted about the inhibition of P. falciparum methyltransferase.
The RuvB-like helicase (Q8ID85_PLAF7) also shortlisted as therapeutic target in the current study. The RuvB-like helicase function like ATP- dependent helicases. It has a vital role in the cell cycle and transcription [61,62,63]. The RUVBL proteins (RUVBL1 & 2) are known to regulate various essential cellular processes in different organisms like Saccharomyces cerevisiae, Drosophila melanogaster and Caenorhabditis elegans [64,65,66]. Three types of RuvB, i.e., PfRuvB1, PfRuvB2, and PfRuvB3 are present in the P. falciparum. The PfRuvB1 possesses ssDNA-stimulated ATPase activity and function as a helicase that unwind the DNA in a 5' to 3' direction . The PfRuvB2 function similar to PfRuvB1, however, its helicase activity is comparatively weak. The PfRuvB3 function only as ATPase with no helicase activity during schizont/merozoits or interaerythocytic mitosis . During the developmental stages of the parasite, the PfRuvB1 and PfRuvB2 are expressed in the asexual phase, while the PfRuvB3 expresses only during the schizont stage, where intraerythrocytic mitosis of P. falciparum occurs . The PfRvuBL3 protein is a true homolog of yeast RuvBL2. Since in yeast, the RuvBL proteins are found extremely indispensable for survival and known to regulate the transcription of almost 5% of yeast genes . The RuvB- like helicases are suitable drug targets to control malaria due to their essentiality for pathogen and non-homology with human host proteome. It is reported that helicases are required for the proliferation of bacteria, viruses and Plasmodium, and inhibiting the DNA unwinding activity reduces the replication of these pathogens in cell cultures and animal models [69,70,71]. The PfRuvB1 ATPase activity is formerly reported to be inhibited by actinomycin, novobiocin, and ethidium bromide .
Among the shortlisted targets, the refined 3D structure of PfAp4AH is available in PDB. It was, therefore, prioritized for drug-like compounds screening based on in silico drug discovery approaches to address potent inhibitors. The PfAp4AH enzyme regulates the levels of signaling molecules, i.e. diadenosine tetraphosphate (Ap4A) by hydrolyzing it to ATP and AMP. This enzyme is localized at the infected RBC membrane in the subpopulation of infected cells . The Ap4A and Ap5A molecules are the chief substrates of PfAp4AH enzyme and key mediators of cellular communication and function through purinergic receptors . Hence, signaling mediated by these molecules within RBCs is of special interest in malaria . The purinergic signaling event is reported to play an important role in parasite invasion . Limited or no comprehensive studies are available about the PfAp4AH inhibition and drug-like compound screening . However, this target been considered worthy for drug discovery and inhibition in other species, such as Mycobacterium tuberculosis .
The recently published essential proteins of P. falciparum were utilized and the comparative proteomics analyses along with in silico druggability approaches were employed to identify novel and suitable drug targets against P. falciparum. The study based on comparative sequence analysis, updated biological databases scanning and multi-direction druggability analyses. This ultimately prioritized and addressed several novel druggable targets against P. falciparum infection not highlighted before. The drug targets and the drug-like compounds prioritized in the current study would be worthy to devise new strategies to combat the P. falciparum drug resistance issues.
Availability of data and materials
The supplementary data relevant to this study is provided.
World Health Organization
Plasmodium falciparum Chloroquine resistant transporter
- Pfk13 gene:
Plasmodium falciparum Kelch-13 gene
Basic Local Alignment Search Tool
Diadenosine tetraphosphate Hydrolase
Protein Data Bank
Plasmodium falciparum Protein phosphatase
Bakhraibah AO. Malaria in modern day: a review article. Austr J Basic Appl Sci. 2018;12:73–5.
Rich S, Leendertz FH, Xu G. LeBreton M, Djoko CF, Aminake MN, et al The origin of malignant malaria. Proc Natl Acad Sci USA. 2009;106:14902–7.
Perkins DJ, Were T, Davenport GC, Kempaiah P, Hittner JB, Ong’echa JM. Severe malarial anemia: innate immunity and pathogenesis. Int J Biol Sci. 2011;7:1427.
Perlmann P, Troye-Blomberg M. Malaria blood-stage infection and its control by the immune system. Folia Biol (Praha). 2000;46:210–8.
Baingana FK, Bos ER. Changing patterns of disease and mortality in Sub-Saharan Africa: an overview. In: Jamison DT, Feachem RG, et al., editors. Disease and mortality in sub-Saharan Africa. 2nd ed. Washington: International Bank for Reconstruction and Development, The World Bank; 2006.
WHO. World malaria report. 20 years of global progress and challenges. Geneva: World Health Organization; 2020. p. 2020.
Bhatia R, Rastogi RM, Ortega L. Malaria successes and challenges in Asia. J Vector Borne Dis. 2013;50:239–47.
Zuber JA, Takala-Harrison S. Multidrug-resistant malaria and the impact of mass drug administration. Infect Drug Resist. 2018;11:299–306.
Lim P, Chy S, Ariey F, Incardona S, Chim P, Sem R, et al. Pfcrt polymorphism and chloroquine resistance in Plasmodium falciparum strains isolated in Cambodia. Antimicrob Agents Chemother. 2003;47:87–94.
Trape JF. The public health impact of chloroquine resistance in Africa. Am J Trop Med Hyg. 2001;64((1_Suppl)):12–7.
Hyde JE. Drug-resistant malaria−an insight. FEBS J. 2007;274:4688–98.
Dondorp AM, Nosten F, Yi P, Das D, Phyo AP, Tarning J, et al. Artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2009;361:455–67.
Noedl H, Se Y, Schaecher K, Smith BL, Socheat D, Fukuda MM. Evidence of artemisinin-resistant malaria in western Cambodia. N Engl J Med. 2008;359:2619–20.
Ouji M, Augereau J-M, Paloque L, Benoit-Vical F. Plasmodium falciparum resistance to artemisinin-based combination therapies: a sword of Damocles in the path toward malaria elimination. Parasite. 2018;25:24.
Dong Y, Wang J, Sun A, Deng Y, Chen M, Xu Y, et al. Genetic association between the Pfk13 gene mutation and artemisinin resistance phenotype in Plasmodium falciparum isolates from Yunnan Province. China Malar J. 2018;17:478.
Zhang M, Wang C, Otto TD, Oberstaller J, Liao X, Adapa SR, et al. Uncovering the essential genes of the human malaria parasite Plasmodium falciparum by saturation mutagenesis. Science. 2018;360:eaap7847.
Bushell E, Gomes AR, Sanderson T, Anar B, Girling G, Herd C, et al. Functional profiling of a Plasmodium genome reveals an abundance of essential genes. Cell. 2017;170:260–72.
Huang Y, Niu B, Gao Y, Fu L, Li W. CD-HIT Suite: a web server for clustering and comparing biological sequences. Bioinformatics. 2010;26:680–2.
Johnson M, Zaretskaya I, Raytselis Y, Merezhuk Y, McGinnis S, Madden TL. NCBI BLAST: a better web interface. Nucleic Acids Res. 2008;36((suppl_2)):W5–9.
Mukherjee S, Gangopadhyay K, Mukherjee SB. Identification of potential new vaccine candidates in Salmonella typhi using reverse vaccinology and subtractive genomics-based approach. bioRxiv. 2020. https://doi.org/10.1101/521518.
Hasan MA, Khan MA, Sharmin T, Mazumder MHH, Chowdhury AS. Identification of putative drug targets in Vancomycin-resistant Staphylococcus aureus (VRSA) using computer aided protein data analysis. Gene. 2016;575:132–43.
Aslam M, Shehroz M, Shah M, Khan MA, Afridi SG, Khan A. Potential druggable proteins and chimeric vaccine construct prioritization against Brucella melitensis from species core genome data. Genomics. 2020;112:1734–45.
Sayers S, Li L, Ong E, Deng S, Fu G, Lin Y, et al. Victors: a web-based knowledge base of virulence factors in human and animal pathogens. Nucleic Acids Res. 2019;47:D693–700.
Wishart DS, Feunang YD, Guo AC, Lo EJ, Marcu A, Grant JR, et al. DrugBank 5.0: a major update to the DrugBank database for 2018. Nucleic Acids Res. 2018;46:D1074–82.
Berman HM, Battistuz T, Bhat TN, Bluhm WF, Bourne PE, Burkhardt K, et al. The protein data bank. Acta Crystallogr D Biol Crystallogr. 2002;58:899–907.
Uddin R, Jamil F. Prioritization of potential drug targets against P. aeruginosa by core proteomic analysis using computational subtractive genomics and Protein-Protein interaction network. Comput Biol Chem. 2018;74:115–22.
Waterhouse A, Bertoni M, Bienert S, Studer G, Tauriello G, Gumienny R, et al. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 2018;46:W296–303.
Dym O, Eisenberg D, Yeates T. ERRAT: Detection of errors in protein models. Int Tables Crystallogr. 2012. https://doi.org/10.1107/97809553602060000881.
Kleywegt GJ. Jones TA. Phi/psi-chology: Ramachandran revisited. Structure. 1996;4:1395–400.
Hussein HA, Borrel A, Geneix C, Petitjean M, Regad L, Camproux A-C. PockDrug-Server: a new web server for predicting pocket druggability on holo and apo proteins. Nucleic Acids Res. 2015;43:W436–42.
Yu CS, Chen YC, Lu CH, Hwang JK. Prediction of protein subcellular localization. Proteins. 2006;64:643–51.
Artimo P, Jonnalagedda M, Arnold K, Baratin D, Csardi G, De Castro E, et al. ExPASy: SIB bioinformatics resource portal. Nucleic Acids Res. 2012;40:W597–603.
Mering CV, Huynen M, Jaeggi D, Schmidt S, Bork P, Snel B. STRING: a database of predicted functional associations between proteins. Nucleic Acids Res. 2003;31:258–61.
Sunseri J, Koes DR. Pharmit: interactive exploration of chemical space. Nucleic Acids Res. 2016;44:W442–8.
Liu Y, Grimm M, Dai W, Hou M, Xiao ZX, Cao Y. CB-Dock: a web server for cavity detection-guided protein–ligand blind docking. Acta Pharmacol Sinica. 2020;41:138–44.
Discovery Studio Visualizer Software Version 4.0. 2012. http://www.accelrys.com.
Oberstaller J, Otto TD, Rayner JC, Adams JH. Essential genes of the parasitic Apicomplexa. Trends Parasitol. 2021;37:304–16.
Rowland I, Gibson G, Heinken A, Scott K, Swann J, Thiele I, et al. Gut microbiota functions: metabolism of nutrients and other food components. Eur J Nutr. 2018;57:1–24.
Yamashiro Y. Gut microbiota in health and disease. Ann Nutr Metabol. 2017;71:242–6.
Solanki V, Tiwari V. Subtractive proteomics to identify novel drug targets and reverse vaccinology for the development of chimeric vaccine against Acinetobacter baumannii. Sci Rep. 2018;8:9044.
Hadizadeh M, Tabatabaiepour SN, Tabatabaiepour SZ, Hosseini Nave H, Mohammadi M, Sohrabi SM. Genome-wide identification of potential drug target in enterobacteriaceae family: a homology-based method. Microb Drug Resist. 2018;24:8–17.
Messaoudi A, Belguith H, Hamida JB. Homology modeling and virtual screening approaches to identify potent inhibitors of VEB-1 β-lactamase. Theor Biol Med Model. 2013;10:22.
Sharma A, Yogavel M, Sharma A. Structural and functional attributes of malaria parasite diadenosine tetraphosphate hydrolase. Sci Rep. 2016;6:19981.
Gaudet P, Livstone MS, Lewis SE, Thomas PD. Phylogenetic-based propagation of functional annotations within the Gene Ontology consortium. Brief Bioinforma. 2011;12:449–62.
Guranowski A, Jakubowski H, Holler E. Catabolism of diadenosine 5’, 5 “‘-P1, P4-tetraphosphate in procaryotes Purification and properties of diadenosine 5’, 5”’-P1, P4-tetraphosphate (symmetrical) pyrophosphohydrolase from Escherichia coli K12. J Biol Chem. 1983;258:14784–9.
Osman W, Endo S, Oh-Hashi K, Kitamura Y, Kitade Y. Molecular characterization and mutational analysis of recombinant diadenosine 5′, 5 ″-P1, P4-tetraphosphate hydrolase from Plasmodium falciparum. Biol Pharm Bull. 2012;35:1191–6.
Guranowski A, Starzyńska E, Pietrowska-Borek M, Rejman D, Blackburn GM. Novel diadenosine polyphosphate analogs with oxymethylene bridges replacing oxygen in the polyphosphate chain: potential substrates and/or inhibitors of Ap4A hydrolases. FEBS J. 2009;276:1546–53.
Guranowski A, Starzyska E, McLennan AG, Baraniak J, Stec WJ. Adenosine-5’-O-phosphorylated and adenosine-5’-O-phosphorothioylated polyols as strong inhibitors of (symmetrical) and (asymmetrical) dinucleoside tetraphosphatases. Biochem J. 2003;373:635–40.
Shi Y. Serine/threonine phosphatases: mechanism through structure. Cell. 2009;139:468–84.
Yang C, Arrizabalaga G. The serine/threonine phosphatases of apicomplexan parasites. Mol Microbiol. 2017;106:1–21.
Yokoyama D, Saito-Ito A, Asao N, Tanabe K, Yamamoto M, Matsumura T. Modulation of the Growth of Plasmodium falciparum in Vitroby protein serine/threonine phosphatase inhibitors. Biochem Biophys Res Commun. 1998;247:18–23.
Khalife J, Pierrot C. Phosphatases are emerging as novel druggable targets in Plasmodium. Future Microbiol. 2016;11:603–6.
Dounay A, Forsyth C. Okadaic acid: the archetypal serine/threonine protein phosphatase inhibitor. Curr Med Chem. 2002;9:1939–80.
Zhang M, Mishra S, Sakthivel R, Fontoura BM, Nussenzweig V. UIS2: a unique phosphatase required for the development of Plasmodium liver stages. PLoS Pathog. 2016;12:e1005370.
Malvisi L. Functional characterization of cytochrome b5 reductase and its electron acceptor cytochrome b5 in Plasmodium falciparum. Dissertation, University of South Florida; 2009.
Maier AG, Cooke BM, Cowman AF, Tilley L. Malaria parasite proteins that remodel the host erythrocyte. Nat Rev Microbiol. 2009;7:341–54.
Iqbal MS, Siddiqui AA, Banerjee C, Nag S, Mazumder S, De R, et al. Detection of retromer assembly in Plasmodium falciparum by immunosensing coupled to Surface Plasmon Resonance. Biochim Biophys Acta Proteins Proteom. 2018;1866:722–30.
Iqbal MS, Siddiqui AA, Alam A, Goyal M, Banerjee C, Sarkar S, et al. Expression, purification and characterization of Plasmodium falciparum vacuolar protein sorting 29. Protein Expr Purif. 2016;120:7–15.
Copeland RA, Solomon ME, Richon VM. Protein methyltransferases as a target class for drug discovery. Nat Rev Drug Discov. 2009;8:724–32.
Brecher M, Chen H, Liu B, Banavali NK, Jones SA, Zhang J, et al. Novel broad spectrum inhibitors targeting the flavivirus methyltransferase. PLoS ONE. 2015;10:0130062.
Ahmad M, Tuteja R. Plasmodium falciparum RuvB proteins: emerging importance and expectations beyond cell cycle progression. Commun Integr Biol. 2012;5:350–61.
Ahmad M, Tuteja R. Plasmodium falciparum RuvB1 is an active DNA helicase and translocates in the 5′–3′ direction. Gene. 2013;515:99–109.
Jha S, Dutta A. RVB1/RVB2: running rings around molecular biology. Mol Cell. 2009;34:521–33.
Qiu X-B, Lin Y-L, Thome KC, Pian P, Schlegel BP, Weremowicz S, et al. An eukaryotic RuvB-like protein (RUVBL1) essential for growth. J Biol Chem. 1998;273:27786–93.
Bauer A, Chauvet S, Huber O, Usseglio F, Rothbächer U, Aragnol D, et al. Pontin52 and Reptin52 function as antagonistic regulators of β-catenin signalling activity. EMBO J. 2000;19:6121–30.
Matias PM, Gorynia S, Donner P, Carrondo MA. Crystal structure of the human AAA+ protein RuvBL1. J Biol Chem. 2006;281:38918–29.
Ahmad M, Tuteja R. Plasmodium falciparum RuvB2 translocates in 5′–3′ direction, relocalizes during schizont stage and its enzymatic activities are up regulated by RuvB3 of the same complex. Biochim Biophys Acta Proteins Proteom. 2013;1834:2795–811.
Ahmad M, Singh S, Afrin F, Tuteja R. Novel RuvB nuclear ATPase is specific to intraerythrocytic mitosis during schizogony of Plasmodium falciparum. Mol Biochem Parasitol. 2012;185:58–65.
Borowski P, Lang M, Haag A, Schmitz H, Choe J, Chen H-M, et al. Characterization of imidazo[4, 5-d]pyridazine nucleosides as modulators of unwinding reaction mediated by West Nile virus nucleoside triphosphatase/helicase: evidence for activity on the level of substrate and/or enzyme. Antimicrob Agents Chemother. 2002;46:1231–9.
Frick DN. Helicases as antiviral drug targets. Drug News Perspect. 2003;16:355–62.
Tuteja R. Helicases− feasible antimalarial drug target for Plasmodium falciparum. FEBS J. 2007;274:4699–704.
Ahmad M, Tarique M, Afrin F, Tuteja N, Tuteja R. Identification of inhibitors of Plasmodium falciparum RuvB1 helicase using biochemical assays. Protoplasma. 2015;252:117–25.
McLennan A. The Nudix hydrolase superfamily. Cell Mol Life Sci. 2006;63:123–43.
Mildvan A, Xia Z, Azurmendi H, Saraswat V, Legler P, Massiah M, et al. Structures and mechanisms of Nudix hydrolases. Arch Biochem Biophys. 2005;433:129–43.
Levano-Garcia J, Dluzewski AR, Markus RP, Garcia CRS. Purinergic signalling is involved in the malaria parasite Plasmodium falciparum invasion to red blood cells. Purinergic Signal. 2010;6:365–72.
Branson KM, Mertens HD, Swarbrick JD, Fletcher JI, Kedzierski L, Gayler KR, et al. Discovery of inhibitors of lupin diadenosine 5 tetraphosphate hydrolase by virtual screening ′, 5′′′-P 1, P 4. Biochemistry. 2009;48:7614–20.
Alexander SP, Kelly E, Mathie A, Peters JA, Veale EL, Armstrong JF, et al. The concise guide to pharmacology 2019/20: introduction and other protein targets. Br J Pharmacol. 2019;176:S1-20.
Ethics approval and consent to participate
Consent for publication
All authors read and agreed to publish the study.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Plasmodium falciparum strain 3D7 non paralogous essential proteins.
Plasmodium falciparum strain 3D7 proteins, non-homologous to human as well as human’s gut flora proteins.
Plasmodium falciparum 3D7 essential proteins homologs to human gut microbiome proteome.
The approved antimalarial targets obtained from Alexander, S. P., et al. (2019) . The accession IDs in Bold show approved antimalarial targets recovered in the list of 183 prioritized targets (Table S2), while accession IDs in italic and asterisk* are approved targets recovered by relaxing the human microbiome database screening criteria to 35% sequence identity and query coverage.
About this article
Cite this article
Ali, F., Wali, H., Jan, S. et al. Analysing the essential proteins set of Plasmodium falciparum PF3D7 for novel drug targets identification against malaria. Malar J 20, 335 (2021). https://doi.org/10.1186/s12936-021-03865-1