Toxicology International (Formerly Indian Journal of Toxicology)

Open Access Open Access  Restricted Access Subscription Access
Open Access Open Access Open Access  Restricted Access Restricted Access Subscription Access

Toxicodynamics of Organophosphates with Human Acetylcholinesterase Interaction at Novel Site Trp-86 for Antidote Action


Affiliations
1 Amity Institute of Environmental Toxicology, Safety and Management, Institute of Engineering, Praca General Tiburcio 80, 22290-270, Rio de Janeiro, Brazil
2 Laboratory of Molecular Modeling Applied to Chemical and Biological Defense (LMACBD), Military Institute of Engineering, Praca General Tiburcio 80, 22290-270, Rio de Janeiro, Brazil
3 Amity Institute of Biotechnology, Amity University, Noida, India
     

   Subscribe/Renew Journal


Molecular Docking, site directed mutagenesis and molecular dynamic (MD) simulation approaches were used to explore mode of binding and inhibition for human acetylcholinesterase (hAChE) and organophosphates (OPs). More than 200 OPs molecules were investigated using Glide docking module of Schrodinger suite as co-crystal structure between two are not available in protein data bank. In initial screening Trp86 was found to be involved in maximum π-cation interaction on anionic subsite of hAChE other than Ser203 (catalytic site). With extra precision glide docking phoxim ethyl phosphonate (PEP) tops among 200 OPs based on glide docking score while interacted with Trp86, Gly121 and Ser203 whereas MMGBSA score shows less binding affinity than heptenophos and dichlorovos. Trp86 preferred π interaction with ring bearing OPs and hydrophobic interactions with smaller OPs without ring bearing structures. Site directed mutagenesis at Trp86 (Trp86 to Ala86) shown the deterioration of the binding site in terms of size reduction, loss of electrostatic and geometric stabilization in binding cavity and significant reduction in binding of OPs in preferred orientation. Dock score of both wild and mutated hAChE shows a perfect qualitative agreement (R2=64.1%) towards the study. Molecular dynamic simulation (GROMACS 4.5.5) of hAChE-PEP complex for 4×104 pico-second with SPC16 water system at 310K temperature explained the evident role of Trp86 in stabilizing the ligand at P-site of the enzyme. Asp74 and Tyr 124 were noticed in conveying H-bonds. Trp86 have shown consistent and stable distance between residues and ligand. Asp74 and Tyr124 appeared as important residues which contributed H bonds and Ser203 was expected to be closer for interaction to happen as it disappeared during simulation. Study suggests role of Trp86 on binding site is equally important to take consideration. As residue Trp86 plays significant role, it can be taken into consideration for the studies on development of more efficient antidotes to overcome the case of human poisoning.

Keywords

Docking, Molecular Dynamic Simulation, Human Acetylcholinesterase, Organophosphates.
User
Subscription Login to verify subscription
Notifications
Font Size

  • Boublik Y, Saint-Aguet P, Lougarre A, Arnaud M, Villatte F, Estrada-Mondaca S, Fournier D. Acetylcholinesterase engineering for detection of insecticide residues. Protein Engineering 2002; 15: 43-50.

  • Comroe JH, Todd J, Gammon GD, Leopold IH, Koelle GB, Bodansky O. The effect of di-isopropylfluorophosphate (DFP) upon patients with myasthenia gravis. American Journal of Medical Science 1946; 212: 641–651.

  • De-BleekerJL, DeReuckJL, Willems JL. Neurological aspects of organophosphate poisoning. ClinNeurolNeurosurg 1992; 94(2):93-103.

  • Friesner RA, Banks JL, Murphy RB, Halgren TA, Klicic JJ, Mainz DT, Repasky MP, Knoll EH, Shaw DE, Shelley M, Perry JK, Francis P, Shenkin PS. Glide: A New Approach for Rapid, Accurate Docking and Scoring. Method and Assessment of Docking Accuracy. J Med Chem. 2004; 47:1739– 1749

  • Gabriel Z, Soreq H. Termination and beyond: Acetylcholinesterase as a modulator of synaptic transmission. Cell and Tissue Research 2006; 326:655-669.

  • Getman DK, Eubanks JH, Camp S, Evans GA. Taylor P. The human gene encoding acetylcholinesterase is located on the long arm of chromosome 7. American Journal of Human Genetics 1992; 51:170–177.

  • Haddad L. Winchester J. Clinical management of poisoning and overdose. Philedelphia, WB Saunders.1983; 575-586.

  • 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 ChemSoc 1988; 110(6): 1657-1666.

  • Kryger G, Harel M, Giles M, Toker L, Velan B, Lazar A, Kronman C, Barak D,. Ariel N, Shafferman A, Silman I, Sussman JL. Structures of recombinant native and E202Q mutant human acetylcholinesterase complexed with the snake-venom toxin fasciculin-II. Acta Cryst 2000; D56:1385-1394

  • Kryger G, Silman I, Sussman JL. Structure of acetylcholinesterase complexed with E2020 (Aricept): Implications for the design of new anti-Alzheimer drugs. Structure 1999; 7:297–307.

  • Kua J, Zhang Y, Eslami AC, Butler JR, McCammon JA. Studying the roles of W86, E202, and Y337 in binding of acetylcholine to acetylcholinesterase using a combined molecular dynamics and multiple docking approach. Protein Science: A Publication of the Protein Society 2003;12(12):2675-2684.

  • Lo R, Chandar NB, Kesharwani MK, Jain A, Ganguly B. In Silico Studies in Probing the Role of Kinetic and Structural Effects of Different Drugs for the Reactivation of Tabun-Inhibited AChE. PLoS ONE 2013; 8(12): e79591

  • Manavalan P, Taylor P, Johnson Jr. WC. Circular dichroism studies of acetylcholine sterase conformation. Comparison of the 11 S and 5.6 S species and the differences induced by inhibitory ligands.BBA- Protein Structure and Molecular Enzymology 1985; 829(3):365-370.

  • Namba T, Nolte CT, Jackrel J. Grob D. Poisoning due to organophosphate insecticides.Acuteand chronic manifestations. Am J Med 1971; 50: 475–492.

  • Ordentlich A, Barak D, Kronman C, Flashner Y, Leitner M, Segall Y, Ariel N,Cohen S, Velan B, Shafferman A. Dissection of the human acetylcholinesterase active center determinants of substrate specificity. Identification of residues constituting the anionic site, the hydrophobic site, and the acyl pocket. J BiolChem 1993; 268(23):17083-95.

  • Quinn DM. Acetylcholinesterase: enzyme structure, reaction dynamics, and virtual transition states. Chem. Rev 1987;87:955-979

  • Ranjan A, Ghosh S, Chauhan A and Jindal T. Molecular docking and site directed mutagenic approach to investigate the role of Trp86 of human acetylcholinesterase with organophosphates. International Journal of Pharmaceutical Science and Research 2016; 7(9): 3802-09.

  • Ranjan A, Kumar A, Gulati K, Thakur S, Jindal T. Role of Aromatic Amino Acids in Stabilizing Organophosphate and Human Acetylcholinesterase Complex. Journal of current Pharma Research 2015; 6(4): 1632-1639.

  • Rider JA, Schulman S, Richtern RB, Moeller HD, DuBois KP. Treatment of myasthenia gravis with octamethyl pyrophosphoramide; preliminary report. Journal of American Medical Association 1951; 145(13): 967-972.

  • Sastry GM, Adzhigirey M, Day T, Annabhimoju R, Sherman W. Protein and ligand preparation: Parameters, protocols, and influence on virtual screening enrichments. J Comput Aided Mol Des 2013; 27(3): 221-234.

  • Schrödinger Release 2014-2: LigPrep, version 3.0, Schrödinger, LLC, New York, NY, 2014.

  • Shafferman A, Velan B, Ordentlich A, Kronman C, Grosfeld H, Leitner M, Flashner Y, Cohen S, Barak D, Ariel N. Substrate inhibition of acetylcholinesterase: residues affecting signal transduction from the surface to the catalytic center. The EMBO Journal 1992; 11 (10): 3561 -3568.

  • Small-Molecule Drug Discovery Suite 2014-2: QikProp, version 4.0, Schrödinger, LLC, New York, NY, 2014.

  • SussmanJL, Harel M, Frolow F, Oefner C, Goldman A, Toker L, Silman I. Atomic structure of acetylcholinesterase from Torpedo californica: a prototypic acetylcholine-binding protein. Science1991; 253: 872–879.

  • Worek F, Aurbek N. Thiermann H. Reactivation of organophosphateinhibited human AChE by combinations of obidoxime and HI-6 in vitro. J ApplToxicol2007; 27:582-588.

  • ZHENG QC, CHU HY, NIU RJ, SUN CC. Theoretical studies of interaction models of human acetylcholine esterase with different inhibitors. Science in China Chemistry. 2009; 52(11): 1911-1916


Abstract Views: 454

PDF Views: 0




  • Toxicodynamics of Organophosphates with Human Acetylcholinesterase Interaction at Novel Site Trp-86 for Antidote Action

Abstract Views: 454  |  PDF Views: 0

Authors

Anuj Ranjan
Amity Institute of Environmental Toxicology, Safety and Management, Institute of Engineering, Praca General Tiburcio 80, 22290-270, Rio de Janeiro, Brazil
Joyce S. F. D. De Almeida
Laboratory of Molecular Modeling Applied to Chemical and Biological Defense (LMACBD), Military Institute of Engineering, Praca General Tiburcio 80, 22290-270, Rio de Janeiro, Brazil
Tanos C. C. Franca
Laboratory of Molecular Modeling Applied to Chemical and Biological Defense (LMACBD), Military Institute of Engineering, Praca General Tiburcio 80, 22290-270, Rio de Janeiro, Brazil
Abhishek Chauhan
Amity Institute of Environmental Toxicology, Safety and Management, Institute of Engineering, Praca General Tiburcio 80, 22290-270, Rio de Janeiro, Brazil
Saradindu Ghosh
Amity Institute of Biotechnology, Amity University, Noida, India
Tanu Jindal
Amity Institute of Environmental Toxicology, Safety and Management, Institute of Engineering, Praca General Tiburcio 80, 22290-270, Rio de Janeiro, Brazil

Abstract


Molecular Docking, site directed mutagenesis and molecular dynamic (MD) simulation approaches were used to explore mode of binding and inhibition for human acetylcholinesterase (hAChE) and organophosphates (OPs). More than 200 OPs molecules were investigated using Glide docking module of Schrodinger suite as co-crystal structure between two are not available in protein data bank. In initial screening Trp86 was found to be involved in maximum π-cation interaction on anionic subsite of hAChE other than Ser203 (catalytic site). With extra precision glide docking phoxim ethyl phosphonate (PEP) tops among 200 OPs based on glide docking score while interacted with Trp86, Gly121 and Ser203 whereas MMGBSA score shows less binding affinity than heptenophos and dichlorovos. Trp86 preferred π interaction with ring bearing OPs and hydrophobic interactions with smaller OPs without ring bearing structures. Site directed mutagenesis at Trp86 (Trp86 to Ala86) shown the deterioration of the binding site in terms of size reduction, loss of electrostatic and geometric stabilization in binding cavity and significant reduction in binding of OPs in preferred orientation. Dock score of both wild and mutated hAChE shows a perfect qualitative agreement (R2=64.1%) towards the study. Molecular dynamic simulation (GROMACS 4.5.5) of hAChE-PEP complex for 4×104 pico-second with SPC16 water system at 310K temperature explained the evident role of Trp86 in stabilizing the ligand at P-site of the enzyme. Asp74 and Tyr 124 were noticed in conveying H-bonds. Trp86 have shown consistent and stable distance between residues and ligand. Asp74 and Tyr124 appeared as important residues which contributed H bonds and Ser203 was expected to be closer for interaction to happen as it disappeared during simulation. Study suggests role of Trp86 on binding site is equally important to take consideration. As residue Trp86 plays significant role, it can be taken into consideration for the studies on development of more efficient antidotes to overcome the case of human poisoning.

Keywords


Docking, Molecular Dynamic Simulation, Human Acetylcholinesterase, Organophosphates.

References





DOI: https://doi.org/10.22506/ti%2F2017%2Fv24%2Fi1%2F149042