Deciphering structure, function and mechanism of Plasmodium IspD homologs from their evolutionary imprints

Malaria is a life-threatening mosquito-borne blood disease caused by infection with Plasmodium parasites. Anti-malarial drug resistance is a global threat to control and eliminate malaria and therefore, it is very important to discover and evaluate new drug targets. The 2-C-methyl-d-erythritol 4-phosphate cytidylyltransferase (IspD) homolog is a second in vivo target for fosmidomycin within isoprenoid biosynthesis in malarial parasites. In the present study, we have deciphered the sequence-structure–function integrity of IspD homologs based on their evolutionary imprints. The function and catalytic mechanism of them were also intensively studied by using sequence-structure homology, molecular modeling, and docking approach. Results of our study indicated that substrate-binding and dimer interface motifs in their structures were extensively conserved and part of them closely related to eubacterial origins. Amino acid substitutions in their coiled-coil regions found to bring a radical change in secondary structural elements, which in turn may change the local structural environment. Arg or Asp was identified as a catalytic site in plasmodium IspD homologs, contributing a direct role in the cytidylyltransferase activity similar to bacterial IspD. Results of molecular docking studies demonstrated how anti-malarial drugs such as fosmidomycin and FR-900098 have competitively interacted with the substrate-binding site of these homologs. As shown by our analysis, species-specific evolutionary imprints in these homologs determine the sequence-structure–function-virulence integrity and binding site alterations in order to confer anti-malarial drug resistance.

[1]  R. O. Morgan,et al.  Deciphering function and mechanism of calcium-binding proteins from their evolutionary imprints. , 2006, Biochimica et biophysica acta.

[2]  Michael J. E. Sternberg,et al.  3DLigandSite: predicting ligand-binding sites using similar structures , 2010, Nucleic Acids Res..

[3]  W. V. Van Voorhis,et al.  Functional genetic analysis of the Plasmodium falciparum deoxyxylulose 5-phosphate reductoisomerase gene. , 2010, Molecular and biochemical parasitology.

[4]  W. Eisenreich,et al.  Biochemistry of the non-mevalonate isoprenoid pathway , 2011, Cellular and Molecular Life Sciences.

[5]  P. Rathod,et al.  Antimalarial Drug Resistance: Literature Review and Activities and Findings of the ICEMR Network , 2015, The American journal of tropical medicine and hygiene.

[6]  G. McFadden,et al.  Plasmodium falciparum apicoplast drugs: targets or off-targets? , 2012, Chemical reviews.

[7]  D. Cane,et al.  Structure of 4-diphosphocytidyl-2-C- methylerythritol synthetase involved in mevalonate- independent isoprenoid biosynthesis , 2001, Nature Structural Biology.

[8]  R. Altman,et al.  Computational analysis of Plasmodium falciparum metabolism: organizing genomic information to facilitate drug discovery. , 2004, Genome research.

[9]  Yang Zhang,et al.  Improving the physical realism and structural accuracy of protein models by a two-step atomic-level energy minimization. , 2011, Biophysical journal.

[10]  Yang Zhang,et al.  COFACTOR: an accurate comparative algorithm for structure-based protein function annotation , 2012, Nucleic Acids Res..

[11]  Christopher M. Armstrong,et al.  Plasmodium IspD (2-C-Methyl-D-erythritol 4-Phosphate Cytidyltransferase), an Essential and Druggable Antimalarial Target. , 2015, ACS infectious diseases.

[12]  G. Klebe,et al.  Experimental and Computational Active Site Mapping as a Starting Point to Fragment‐Based Lead Discovery , 2012, ChemMedChem.

[13]  H. Lichtenthaler,et al.  Inhibitors of the nonmevalonate pathway of isoprenoid biosynthesis as antimalarial drugs. , 1999, Science.

[14]  Johannes Söding,et al.  Protein homology detection by HMM?CHMM comparison , 2005, Bioinform..

[15]  Christopher J. Tonkin,et al.  Evolution of malaria parasite plastid targeting sequences , 2008, Proceedings of the National Academy of Sciences.

[16]  Marco Biasini,et al.  SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information , 2014, Nucleic Acids Res..

[17]  Xavier Soberón,et al.  Evolutionary walk between (β/α)(8) barrels: catalytic migration from triosephosphate isomerase to thiamin phosphate synthase. , 2012, Journal of molecular biology.

[18]  Amos Bairoch,et al.  ScanProsite: detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins , 2006, Nucleic Acids Res..

[19]  W. Eisenreich,et al.  Cytidine 5'-triphosphate-dependent biosynthesis of isoprenoids: YgbP protein of Escherichia coli catalyzes the formation of 4-diphosphocytidyl-2-C-methylerythritol. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[20]  Guido Tiana,et al.  Imprint of evolution on protein structures. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[21]  J. Thompson,et al.  CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. , 1994, Nucleic acids research.

[22]  M. Zhang,et al.  Biosynthesis of isoprenoids: characterization of a functionally active recombinant 2-C-methyl-D-erythritol 4-phosphate cytidyltransferase (IspD) from Mycobacterium tuberculosis H37Rv. , 2007, Journal of biochemistry and molecular biology.

[23]  W. Eisenreich,et al.  The crystal structure of a plant 2C‐methyl‐D‐erythritol 4‐phosphate cytidylyltransferase exhibits a distinct quaternary structure compared to bacterial homologues and a possible role in feedback regulation for cytidine monophosphate , 2006, The FEBS journal.

[24]  Christopher J. Tonkin,et al.  Dissecting Apicoplast Targeting in the Malaria Parasite Plasmodium falciparum , 2003, Science.

[25]  Yang Zhang,et al.  BSP‐SLIM: A blind low‐resolution ligand‐protein docking approach using predicted protein structures , 2012, Proteins.

[26]  J Zuegge,et al.  Deciphering apicoplast targeting signals--feature extraction from nuclear-encoded precursors of Plasmodium falciparum apicoplast proteins. , 2001, Gene.

[27]  Jason A. Papin,et al.  Novel Plasmodium falciparum metabolic network reconstruction identifies shifts associated with clinical antimalarial resistance , 2017, bioRxiv.

[28]  K. Sumathi,et al.  3dSS: 3D structural superposition , 2006, Nucleic Acids Res..

[29]  C. Hoogland,et al.  In The Proteomics Protocols Handbook , 2005 .

[30]  J. Thornton,et al.  Missing in action: enzyme functional annotations in biological databases. , 2009, Nature chemical biology.

[31]  Z. R. Li,et al.  PEARLS: Program for Energetic Analysis of Receptor-Ligand System , 2006, J. Chem. Inf. Model..

[32]  T. A. Hall,et al.  BIOEDIT: A USER-FRIENDLY BIOLOGICAL SEQUENCE ALIGNMENT EDITOR AND ANALYSIS PROGRAM FOR WINDOWS 95/98/ NT , 1999 .

[33]  J. Thompson,et al.  The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. , 1997, Nucleic acids research.

[34]  F. Diederich,et al.  Inhibitors of the herbicidal target IspD: allosteric site binding. , 2011, Angewandte Chemie.

[35]  S. Parija,et al.  Antimalarial drug resistance: An overview , 2016, Tropical parasitology.

[36]  K Henrick,et al.  Electronic Reprint Biological Crystallography Secondary-structure Matching (ssm), a New Tool for Fast Protein Structure Alignment in Three Dimensions Biological Crystallography Secondary-structure Matching (ssm), a New Tool for Fast Protein Structure Alignment in Three Dimensions , 2022 .

[37]  Narmada Thanki,et al.  CDD: a Conserved Domain Database for the functional annotation of proteins , 2010, Nucleic Acids Res..

[38]  Nidhi Singh,et al.  Targeting the methyl erythritol phosphate (MEP) pathway for novel antimalarial, antibacterial and herbicidal drug discovery: inhibition of 1-deoxy-D-xylulose-5-phosphate reductoisomerase (DXR) enzyme. , 2007, Current pharmaceutical design.

[39]  Wolfgang Eisenreich,et al.  Biosynthesis of isoprenoids: a bifunctional IspDF enzyme from Campylobacter jejuni. , 2004, European journal of biochemistry.

[40]  Janet M. Thornton,et al.  ProFunc: a server for predicting protein function from 3D structure , 2005, Nucleic Acids Res..

[41]  J M Thornton,et al.  LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions. , 1995, Protein engineering.

[42]  Yang Zhang,et al.  I-TASSER: a unified platform for automated protein structure and function prediction , 2010, Nature Protocols.

[43]  Christophe Geourjon,et al.  SOPMA: significant improvements in protein secondary structure prediction by consensus prediction from multiple alignments , 1995, Comput. Appl. Biosci..

[44]  D. Cane,et al.  Kinetic analysis of Escherichia coli 2-C-methyl-D-erythritol-4-phosphate cytidyltransferase, wild type and mutants, reveals roles of active site amino acids. , 2004, Biochemistry.

[45]  Neil Hall,et al.  Evolutionary pressures on apicoplast transit peptides. , 2004, Molecular biology and evolution.

[46]  K. M. Watts,et al.  A second target of the antimalarial and antibacterial agent fosmidomycin revealed by cellular metabolic profiling. , 2011, Biochemistry.

[47]  W. Hunter,et al.  Characterization of the Mycobacterium tuberculosis 4-Diphosphocytidyl-2-C-Methyl-d-Erythritol Synthase: Potential for Drug Development , 2007, Journal of bacteriology.

[48]  Tal Pupko,et al.  ConSurf 2010: calculating evolutionary conservation in sequence and structure of proteins and nucleic acids , 2010, Nucleic Acids Res..

[49]  M. Nei,et al.  MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. , 2011, Molecular biology and evolution.

[50]  Peer Bork,et al.  SMART 7: recent updates to the protein domain annotation resource , 2011, Nucleic Acids Res..

[51]  T. Jones,et al.  Structural and functional studies of mycobacterial IspD enzymes. , 2011, Acta crystallographica. Section D, Biological crystallography.

[52]  P. Babbitt,et al.  Divergent evolution of enzymatic function: mechanistically diverse superfamilies and functionally distinct suprafamilies. , 2001, Annual review of biochemistry.

[53]  W. Eisenreich,et al.  Hexameric Assembly of the Bifunctional Methylerythritol 2,4-Cyclodiphosphate Synthase and Protein-Protein Associations in the Deoxy-xylulose-dependent Pathway of Isoprenoid Precursor Biosynthesis* , 2004, Journal of Biological Chemistry.

[54]  J. Wiesner,et al.  FR-900098, an antimalarial development candidate that inhibits the non-mevalonate isoprenoid biosynthesis pathway, shows no evidence of acute toxicity and genotoxicity , 2016, Virulence.

[55]  J. Marles-Wright,et al.  Synthesis of Cdp-activated Ribitol For , 2008 .

[56]  W. Hunter,et al.  Structure of a tetragonal crystal form of Escherichia coli 2-C-methyl-D-erythritol 4-phosphate cytidylyltransferase. , 2003, Acta crystallographica. Section D, Biological crystallography.

[57]  H. Balaram,et al.  Triosephosphate isomerase from Plasmodium falciparum: the crystal structure provides insights into antimalarial drug design. , 1997, Structure.

[58]  G. McFadden,et al.  The evolution, metabolism and functions of the apicoplast , 2010, Philosophical Transactions of the Royal Society B: Biological Sciences.

[59]  J. Wiesner,et al.  Plasmodium falciparum: detection of the deoxyxylulose 5-phosphate reductoisomerase activity. , 2000, Experimental parasitology.

[60]  Ramaswamy Senthilkumar,et al.  PSAP: protein structure analysis package , 2007 .

[61]  B. Clantin,et al.  X-ray structure of papaya chitinase reveals the substrate binding mode of glycosyl hydrolase family 19 chitinases. , 2008, Biochemistry.

[62]  J. Rubio-Martínez,et al.  Homology modeling of Mycobacterium tuberculosis 2C-methyl-d-erythritol-4-phosphate cytidylyltransferase, the third enzyme in the MEP pathway for isoprenoid biosynthesis , 2010, Journal of molecular modeling.

[63]  D. Soldati-Favre,et al.  Metabolic pathways in the apicoplast of apicomplexa. , 2010, International review of cell and molecular biology.

[64]  A. Hemmerlin,et al.  Isoprenoid biosynthesis as a target for antibacterial and antiparasitic drugs: phosphonohydroxamic acids as inhibitors of deoxyxylulose phosphate reducto-isomerase. , 2005, The Biochemical journal.

[65]  Charles Elkan,et al.  Fitting a Mixture Model By Expectation Maximization To Discover Motifs In Biopolymer , 1994, ISMB.

[66]  Yu Wang,et al.  Identification of inhibitors for putative malaria drug targets among novel antimalarial compounds. , 2011, Molecular and biochemical parasitology.

[67]  R. Sabarinathan,et al.  Evolution, Homology Conservation, and Identification of Unique Sequence Signatures in GH19 Family Chitinases , 2010, Journal of Molecular Evolution.

[68]  R D Appel,et al.  Protein identification and analysis tools in the ExPASy server. , 1999, Methods in molecular biology.