BACKGROUND: Ticks counteract host inflammatory responses by secreting proteins from their saliva that compete for histamine binding. Among these tick salivary proteins are lipocalins, antiparallel beta-barrel proteins that sequester small molecules. A tick salivary lipocalin has been structurally resolved and experimentally shown to efficiently compete for histamine with its native receptor (e.g., H1 histamine receptor). To date, molecular dynamics simulations focus on protein-protein and protein-ligand interactions, but there are currently no studies for simultaneous ligand exploration between two competing proteins. METHODS: Aided by state-of-the-art, high-throughput computational methods, the current study simulated and analyzed the dynamics of competitive histamine binding at the tick-host interface using the available crystal structures of both the tick salivary lipocalin histamine-binding protein from Rhipicephalus appendiculatus and the human histamine receptor 1. RESULTS: The attraction towards the tick salivary lipocalin seems to depend on the protonated (adding a hydrogen ion) state of histamine since the current study shows that as histamine becomes more protonated it increases its exploration for the tick salivary lipocalin. This implies that during tick feeding, histamine may need to be protonated for the tick salivary lipocalin to efficiently sequester it in order to counteract inflammation. Additionally, the beta-hairpin loops (at both ends of the tick salivary lipocalin barrel) were reported to have a functional role in sequestering histamine and the results in the current study concur and provide evidence for this hypothesis. These beta-hairpin loops of the tick salivary lipocalin possess more acidic residues than a structurally similar but functionally unrelated lipocalin from the butterfly, Pieris brassicae; comparative results indicate these acidic residues may be responsible for the ability of the tick lipocalin to out-compete the native (H1) receptor for histamine. CONCLUSIONS: Three explanatory types of data can be obtained from the current study: (i) the dynamics of multiple binding sites, (ii) competition between two proteins for a ligand, and (iii) the intrinsic molecular components involved in the competition. These data can provide further insight at the atomic level of the host-tick interface that cannot be experimentally determined. Additionally, the methods used in this study can be applied in rationally designing drugs.

Institute of Parasitology Biology Centre of the Academy of Sciences of the Czech Republic 37005 České Budějovice Czech Republic

Zobrazit více v PubMed

Falus A. Histamine and Inflammation (Medical Intelligence Unit) Austin, TX: R.G Landes Company; 1994.

Mans BJ, Andersen JF, Francischetti IMB, Valenzuela JG, Schwan TG, Pham VM, Garfield MK, Hammer CH, Ribeiro JMC. Comparative sialomics between hard and soft ticks: Implications for the evolution of blood-feeding behavior. Insect Biochem Mol Biol. 2008;38(1):42–58. doi: 10.1016/j.ibmb.2007.09.003. PubMed DOI PMC

Konnai S, Nishikado H, Yamada S, Imamura S, Ito T, Onuma M, Murata S, Ohashi K. Molecular identification and expression analysis of lipocalins from blood feeding taiga tick, Ixodes persulcatus Schulze. Exp Parasitol. 2011;127(2):467–474. doi: 10.1016/j.exppara.2010.10.002. PubMed DOI

Heekin A, Guerrero F, Bendele K, Saldivar L, Scoles G, Gondro C, Nene V, Djikeng A, Brayton K. Analysis of Babesia bovis infection-induced gene expression changes in larvae from the cattle tick, Rhipicephalus (Boophilus) microplus. Parasit Vectors. 2012;5(1):162. doi: 10.1186/1756-3305-5-162. PubMed DOI PMC

Sangamnatdej S, Paesen GC, Slovak M, Nuttall PA. A high affinity serotonin- and histamine-binding lipocalin from tick saliva. Insect Mol Biol. 2002;11(1):79–86. doi: 10.1046/j.0962-1075.2001.00311.x. PubMed DOI

Paesen GC, Adams PL, Harlos K, Nuttall PA, Stuart DI. Tick histamine-binding proteins: isolation, cloning, and three-dimensional structure. Mol Cell. 1999;3(5):661–671. doi: 10.1016/S1097-2765(00)80359-7. PubMed DOI

Garcia G, Gardinassi L, Ribeiro J, Anatriello E, Ferreira B, Moreira HN, Mafra C, Martins M, Szabo MP, de Miranda-Santos IK, Maruyama S. The sialotranscriptome of Amblyomma triste, Amblyomma parvum and Amblyomma cajennense ticks, uncovered by 454-based RNA-seq. Parasit Vectors. 2014;7(1):430. doi: 10.1186/1756-3305-7-430. PubMed DOI PMC

Borhani D, Shaw D. The future of molecular dynamics simulations in drug discovery. J Comput Aided Mol Des. 2012;26(1):15–26. doi: 10.1007/s10822-011-9517-y. PubMed DOI PMC

Arkhipov A, Shan Y, Kim ET, Shaw DE. Membrane interaction of bound ligands contributes to the negative binding cooperativity of the EGF receptor. PLoS Comput Biol. 2014;10(7):e1003742. doi: 10.1371/journal.pcbi.1003742. PubMed DOI PMC

Dror RO, Arlow DH, Maragakis P, Mildorf TJ, Pan AC, Xu H, Borhani DW, Shaw DE. Activation mechanism of the β-adrenergic receptor. Proc Natl Acad Sci U S A. 2011;108(46):18684–18689. doi: 10.1073/pnas.1110499108. PubMed DOI PMC

Oliveira A, Allegri A, Bidon-Chanal A, Knipp M, Roitberg AE, Abbruzzetti S, Viappiani C, Luque FJ. Kinetics and computational studies of ligand migration in nitrophorin 7 and its Δ1–3 mutant. Biochim Biophys Acta Proteins Proteom. 2013;1834(9):1711–1721. doi: 10.1016/j.bbapap.2013.04.009. PubMed DOI

Kondrashov DA, Roberts SA, Weichsel A, Montfort WR. Protein functional cycle viewed at atomic resolution: conformational change and mobility in nitrophorin 4 as a function of pH and NO binding. Biochemistry (Mosc) 2004;43:13637–13647. doi: 10.1021/bi0483155. PubMed DOI

Borrelli KW, Vitalis A, Alcantara R, Guallar V. PELE: protein energy landscape exploration. A novel Monte Carlo based technique. J Chem Theory Comput. 2005;1(6):1304–1311. doi: 10.1021/ct0501811. PubMed DOI

Madadkar-Sobhani A, Guallar V. PELE web server: atomistic study of biomolecular systems at your fingertips. Nucleic Acids Res. 2013;41:W322–W328. doi: 10.1093/nar/gkt454. PubMed DOI PMC

Bowers KJ, Chow E, Huageng X, Dror RO, Eastwood MP, Gregersen BA, Klepeis JL, Kolossvary I, Moraes MA, Sacerdoti FD, Salmon JK, Yibing S, Shaw DE. Proceedings of the ACM/IEEE Conference on Supercomputing (SC06): 11-17 November 2006. Tampa, Florida. New York, NY: ACM/IEEE; 2006. Scalable Algorithms for Molecular Dynamics Simulations on Commodity Clusters; pp. 43–56.

Li X, Jacobson MP, Zhu K, Zhao S, Friesner RA. Assignment of polar states for protein amino acid residues using an interaction cluster decomposition algorithm and its application to high resolution protein structure modeling. Proteins Struct Funct Bioinf. 2007;66(4):824–837. doi: 10.1002/prot.21125. PubMed DOI

Jacobson MP, Pincus DL, Rapp CS, Day TJ, Honig B, Shaw DE, Friesner RA. A hierarchical approach to all-atom protein loop prediction. Proteins Struct Funct Bioinf. 2004;55(2):351–367. doi: 10.1002/prot.10613. PubMed DOI

Atilgan AR, Durell SR, Jernigan RL, Demirel MC, Keskin O, Bahar I. Anisotropy of fluctuation dynamics of proteins with an elastic network model. Biophys J. 2001;80(1):505–515. doi: 10.1016/S0006-3495(01)76033-X. PubMed DOI PMC

Jacobson MP, Friesner RA, Xiang Z, Honig B. On the role of the crystal environment in determining protein side-chain conformations. J Mol Biol. 2002;320(3):597–608. doi: 10.1016/S0022-2836(02)00470-9. PubMed DOI

Still WC, Tempczyk A, Hawley RC, Hendrickson T. Semianalytical treatment of solvation for molecular mechanics and dynamics. J Am Chem Soc. 1990;112(16):6127–6129. doi: 10.1021/ja00172a038. DOI

Jorgensen WL, Tirado-Rives J. The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. J Am Chem Soc. 1988;110(6):1657–1666. doi: 10.1021/ja00214a001. PubMed DOI

Koshland DE. Application of a theory of enzyme specificity to protein synthesis. Proc Natl Acad Sci U S A. 1958;44(2):98–104. doi: 10.1073/pnas.44.2.98. PubMed DOI PMC

Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE. UCSF Chimera—A visualization system for exploratory research and analysis. J Comput Chem. 2004;25(13):1605–1612. doi: 10.1002/jcc.20084. PubMed DOI

Humphrey W, Dalke A, Schulten K. VMD - Visual Molecular Dynamics. J Mol Graph. 1996;14:33–38. doi: 10.1016/0263-7855(96)00018-5. PubMed DOI

Shimamura T, Shiroishi M, Weyand S, Tsujimoto H, Winter G, Katritch V, Abagyan R, Cherezov V, Liu W, Han GW, Kobayashi T, Stevens RC, Iwata S. Structure of the human histamine H1 receptor complex with doxepin. Nature. 2011;475(7354):65–70. doi: 10.1038/nature10236. PubMed DOI PMC

Driver AG, Kukoly CA, Bennett TE. Expression of histamine H1 receptors on cultured histiocytic lymphoma cells. Biochem Pharmacol. 1989;38(18):3083–3091. doi: 10.1016/0006-2952(89)90019-1. PubMed DOI

Wikel S. Ticks and tick-borne pathogens at the cutaneous interface: host defenses, tick countermeasures, and a suitable environment for pathogen establishment. Front Microbiol. 2013;4:337. doi: 10.3389/fmicb.2013.00337. PubMed DOI PMC

Paesen GC, Siebold C, Harlos K, Peacey MF, Nuttall PA, Stuart DI. A tick protein with a modified kunitz fold inhibits human Tryptase. J Mol Biol. 2007;368(4):1172–1186. doi: 10.1016/j.jmb.2007.03.011. PubMed DOI

Sonenshine DE, Roe RM. Biology of Ticks. New York, NY (USA): Oxford University Press; 2013.

Farr M, Wainwright A, Salmon M, Hollywell CA, Bacon PA. Platelets in the synovial fluid of patients with rheumatoid arthritis. Rheumatol Int. 1984;4(1):13–17. doi: 10.1007/BF00683878. PubMed DOI

Valenzuela JG. High-throughput approaches to study salivary proteins and genes from vectors of disease. Insect Biochem Mol Biol. 2002;32(10):1199–1209. doi: 10.1016/S0965-1748(02)00083-8. PubMed DOI

Francischetti IMB, Sá-Nunes A, Mans BJ, Santos IM, Ribeiro JMC. The role of saliva in tick feeding. Front Biosci. 2010;14:2051–2088. PubMed PMC

Substrate prediction of Ixodes ricinus salivary lipocalins differentially expressed during Borrelia afzelii infection