Computational Reverse-Engineering of a Spider-Venom Derived Peptide Active Against Plasmodium falciparum SUB1

Background Psalmopeotoxin I (PcFK1), a protein of 33 aminoacids derived from the venom of the spider Psalmopoeus Cambridgei, is able to inhibit the growth of Plasmodium falciparum malaria parasites with an IC in the low micromolar range. PcFK1 was proposed to act as an ion channel inhibitor, although experimental validation of this mechanism is lacking. The surface loops of PcFK1 have some sequence similarity with the parasite protein sequences cleaved by PfSUB1, a subtilisin-like protease essential for egress of Plasmodium falciparum merozoites and invasion into erythrocytes. As PfSUB1 has emerged as an interesting drug target, we explored the hypothesis that PcFK1 targeted PfSUB1 enzymatic activity. Findings Molecular modeling and docking calculations showed that one loop could interact with the binding site of PfSUB1. The calculated free energy of binding averaged −5.01 kcal/mol, corresponding to a predicted low-medium micromolar constant of inhibition. PcFK1 inhibited the enzymatic activity of the recombinant PfSUB1 enzyme and the in vitro P.falciparum culture in a range compatible with our bioinformatics analysis. Using contact analysis and free energy decomposition we propose that residues A14 and Q15 are important in the interaction with PfSUB1. Conclusions Our computational reverse engineering supported the hypothesis that PcFK1 targeted PfSUB1, and this was confirmed by experimental evidence showing that PcFK1 inhibits PfSUB1 enzymatic activity. This outlines the usefulness of advanced bioinformatics tools to predict the function of a protein structure. The structural features of PcFK1 represent an interesting protein scaffold for future protein engineering.

[1]  C. Sander,et al.  Errors in protein structures , 1996, Nature.

[2]  M. Foley,et al.  Structures of phage-display peptides that bind to the malarial surface protein, apical membrane antigen 1, and block erythrocyte invasion. , 2003, Biochemistry.

[3]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[4]  J. Otlewski,et al.  Canonical protein inhibitors of serine proteases , 2003, Cellular and Molecular Life Sciences (CMLS).

[5]  J. Haynes,et al.  Quantitative assessment of antimalarial activity in vitro by a semiautomated microdilution technique , 1979, Antimicrobial Agents and Chemotherapy.

[6]  Morten Nielsen,et al.  CPH MODELS 2.0: X3M A COMPUTER PROGRAM TO EXTRACT 3D MODELS , 2002 .

[7]  J. Thornton,et al.  Stereochemical quality of protein structure coordinates , 1992, Proteins.

[8]  C. Lambert,et al.  ESyPred 3 D : Prediction of proteins 3 D structures , 2002 .

[9]  W. Trager,et al.  Human malaria parasites in continuous culture. , 1976, Science.

[10]  D. Case,et al.  Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes. , 2003, Journal of molecular biology.

[11]  Tanja Kortemme,et al.  Computer-aided design of functional protein interactions. , 2009, Nature chemical biology.

[12]  The UniProt Consortium,et al.  The Universal Protein Resource (UniProt) 2009 , 2008, Nucleic Acids Res..

[13]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[14]  Bruce Tidor,et al.  Progress in computational protein design. , 2007, Current opinion in biotechnology.

[15]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[16]  H. Ginsburg,et al.  In Vitro Antiplasmodium Effects of Dermaseptin S4 Derivatives , 2002, Antimicrobial Agents and Chemotherapy.

[17]  Sandor Vajda,et al.  Consensus alignment for reliable framework prediction in homology modeling , 2003, Bioinform..

[18]  W. Kabsch,et al.  Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical features , 1983, Biopolymers.

[19]  C. Withers-Martinez,et al.  Subcellular Discharge of a Serine Protease Mediates Release of Invasive Malaria Parasites from Host Erythrocytes , 2007, Cell.

[20]  Jing Zhang,et al.  Development and validation of flow cytometric measurement for parasitemia in cultures of P. falciparum vitally stained with YOYO‐1 , 2007, Cytometry. Part A : the journal of the International Society for Analytical Cytology.

[21]  D. Beveridge,et al.  Affinity and specificity of protein U1A-RNA complex formation based on an additive component free energy model. , 2007, Journal of molecular biology.

[22]  H. Carlson,et al.  Computational studies and peptidomimetic design for the human p53–MDM2 complex , 2004, Proteins.

[23]  T. Egan,et al.  Strategies to reverse drug resistance in malaria , 2007, Current opinion in infectious diseases.

[24]  Christophe G. Lambert,et al.  ESyPred3D: Prediction of proteins 3D structures , 2002, Bioinform..

[25]  M. Sternberg,et al.  Protein structure prediction on the Web: a case study using the Phyre server , 2009, Nature Protocols.

[26]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[27]  W. Delano.  The PyMOL Molecular Graphics System (2002) , 2002 .

[28]  D. Eisenberg,et al.  A method to identify protein sequences that fold into a known three-dimensional structure. , 1991, Science.

[29]  Cathy H. Wu,et al.  The Universal Protein Resource (UniProt) , 2004, Nucleic Acids Res..

[30]  W. Delano.  The PyMOL Molecular Graphics System , 2002 .

[31]  P. Baldacci,et al.  Drug inhibition of HDAC3 and epigenetic control of differentiation in Apicomplexa parasites , 2009, The Journal of experimental medicine.

[32]  D. Goldberg,et al.  Malaria parasite exit from the host erythrocyte: a two-step process requiring extraerythrocytic proteolysis. , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[33]  H. Bujard,et al.  A multifunctional serine protease primes the malaria parasite for red blood cell invasion , 2009, The EMBO journal.

[34]  F. Hackett,et al.  Functional Characterization of the Propeptide of Plasmodium falciparum Subtilisin-like Protease-1* , 2003, Journal of Biological Chemistry.

[35]  T. Yeates,et al.  Verification of protein structures: Patterns of nonbonded atomic interactions , 1993, Protein science : a publication of the Protein Society.

[36]  D. Craik,et al.  The cystine knot motif in toxins and implications for drug design. , 2001, Toxicon : official journal of the International Society on Toxinology.

[37]  D. Ojcius,et al.  Isolation and characterization of Psalmopeotoxin I and II: two novel antimalarial peptides from the venom of the tarantula Psalmopoeus cambridgei , 2004, FEBS letters.

[38]  G Vriend,et al.  WHAT IF: a molecular modeling and drug design program. , 1990, Journal of molecular graphics.

[39]  Dima Kozakov,et al.  Convergence and combination of methods in protein-protein docking. , 2009, Current opinion in structural biology.

[40]  T. Blundell,et al.  Comparative protein modelling by satisfaction of spatial restraints. , 1993, Journal of molecular biology.

[41]  Y. Shai,et al.  Selective Cytotoxicity of Dermaseptin S3 toward IntraerythrocyticPlasmodium falciparum and the Underlying Molecular Basis* , 1997, The Journal of Biological Chemistry.

[42]  G. Ellman,et al.  Tissue sulfhydryl groups. , 1959, Archives of biochemistry and biophysics.

[43]  S. Katiyar,et al.  Perspective in antimalarial chemotherapy. , 2003, Current medicinal chemistry.

[44]  M. Lepšík,et al.  Efficiency of a second‐generation HIV‐1 protease inhibitor studied by molecular dynamics and absolute binding free energy calculations , 2004, Proteins.

[45]  Manfred J. Sippl,et al.  Thirty years of environmental health research--and growing. , 1996, Nucleic Acids Res..

[46]  Fabrice Armougom,et al.  Expresso: automatic incorporation of structural information in multiple sequence alignments using 3D-Coffee , 2006, Nucleic Acids Res..

[47]  D. Case,et al.  Theory and applications of the generalized born solvation model in macromolecular simulations , 2000, Biopolymers.

[48]  B. Kuhlman,et al.  Computational design of affinity and specificity at protein-protein interfaces. , 2009, Current opinion in structural biology.

[49]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[50]  J. Otlewski,et al.  The many faces of protease–protein inhibitor interaction , 2005, The EMBO journal.

[51]  D. Carucci,et al.  Malaria research in the post-genomic era. , 2000, Parasitology today.

[52]  Wilfred F. van Gunsteren,et al.  Basic ingredients of free energy calculations: A review , 2009, J. Comput. Chem..

[53]  D. Kyle,et al.  Synthesis and in vitro studies of novel pyrimidinyl peptidomimetics as potential antimalarial therapeutic agents. , 2002, Journal of medicinal chemistry.

[54]  S. Wodak,et al.  Deviations from standard atomic volumes as a quality measure for protein crystal structures. , 1996, Journal of molecular biology.

[55]  Z. Weng,et al.  ZDOCK: An initial‐stage protein‐docking algorithm , 2003, Proteins.

[56]  M J Sternberg,et al.  Enhancement of protein modeling by human intervention in applying the automatic programs 3D‐JIGSAW and 3D‐PSSM , 2001, Proteins.

[57]  J. Camadro,et al.  Solution structure of PcFK1, a spider peptide active against Plasmodium falciparum , 2006, Protein science : a publication of the Protein Society.

[58]  Manuel C. Peitsch,et al.  SWISS-MODEL: an automated protein homology-modeling server , 2003, Nucleic Acids Res..

[59]  C. Withers-Martinez,et al.  Expression of recombinant Plasmodium falciparum subtilisin-like protease-1 in insect cells. Characterization, comparison with the parasite protease, and homology modeling. , 2002, The Journal of biological chemistry.

[60]  J. E. Hyde,et al.  Drug‐resistant malaria − an insight , 2007, The FEBS journal.