Identification of potential inhibitors of Fasciola gigantica thioredoxin1: computational screening, molecular dynamics simulation, and binding free energy studies

Fasciola gigantica is the causative organism of fascioliasis and is responsible for major economic losses in livestock production globally. F. gigantica thioredoxin1 (FgTrx1) is an important redox-active enzyme involved in maintaining the redox homeostasis in the cell. To identify a potential anti-fasciolid compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,740) against the FgTrx1 structure. The ligands were docked against FgTrx1 and 309 ligands were found to have better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 30 compounds were found to fit well for re-docking studies. After refinement by molecular docking and drug-likeness analysis, three potential inhibitors (ZINC15970091, ZINC9312362, and ZINC9312661) were identified. These three ligands were further subjected to molecular dynamics simulation (MDS) to compare the dynamics and stability of the protein structure after binding of the ligands. The binding free energy analyses were calculated to determine the intermolecular interactions. The results suggested that the two compounds had a binding free energy of –82.237, and –109.52 kJ.mol−1 for compounds with IDs ZINC9312362 and ZINC9312661, respectively. These predicted compounds displayed considerable pharmacological and structural properties to be drug candidates. We concluded that these two compounds could be potential drug candidates to fight against F. gigantica parasites.

[1]  W. Brown,et al.  Bovine T cell responses to recombinant thioredoxin of Fasciola hepatica. , 1999, Veterinary parasitology.

[2]  Rohit Shukla,et al.  Distant Phe345 mutation compromises the stability and activity of Mycobacterium tuberculosis isocitrate lyase by modulating its structural flexibility , 2017, Scientific Reports.

[3]  R. Zadoks,et al.  Confirmation of triclabendazole resistance in liver fluke in the UK , 2012, Veterinary Record.

[4]  G. Coles Anthelmintic Activity of Triclabendazole , 1986, Journal of Helminthology.

[5]  T. Tripathi,et al.  UDP-N-Acetylglucosamine enolpyruvyl transferase (MurA) of Acinetobacter baumannii (AbMurA): Structural and functional properties. , 2017, International journal of biological macromolecules.

[6]  A. Holmgren,et al.  The thioredoxin antioxidant system. , 2014, Free radical biology & medicine.

[7]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[8]  A. Holmgren,et al.  Thiol redox control via thioredoxin and glutaredoxin systems. , 2005, Biochemical Society transactions.

[9]  Ryan G. Coleman,et al.  ZINC: A Free Tool to Discover Chemistry for Biology , 2012, J. Chem. Inf. Model..

[10]  F. Giannerini,et al.  Role of protein -SH groups in redox homeostasis--the erythrocyte as a model system. , 1998, Archives of biochemistry and biophysics.

[11]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997, J. Comput. Chem..

[12]  A. Holmgren,et al.  Antioxidant function of thioredoxin and glutaredoxin systems. , 2000, Antioxidants & redox signaling.

[13]  H. Toet,et al.  Liver fluke vaccines in ruminants: strategies, progress and future opportunities. , 2014, International journal for parasitology.

[14]  H. Raadsma,et al.  Improving animal and human health through understanding liver fluke immunology , 2010, Parasite immunology.

[15]  P. Brindley,et al.  Apoptosis of cholangiocytes modulated by thioredoxin of carcinogenic liver fluke. , 2015, The international journal of biochemistry & cell biology.

[16]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[17]  A. Holmgren,et al.  Glutaredoxins: glutathione-dependent redox enzymes with functions far beyond a simple thioredoxin backup system. , 2004, Antioxidants & redox signaling.

[18]  Rajendra Kumar,et al.  g_mmpbsa - A GROMACS Tool for High-Throughput MM-PBSA Calculations , 2014, J. Chem. Inf. Model..

[19]  T. Tripathi,et al.  A combined biochemical and computational studies of the rho-class glutathione s-transferase sll1545 of Synechocystis PCC 6803. , 2017, International journal of biological macromolecules.

[20]  A. Loukas,et al.  The secreted and surface proteomes of the adult stage of the carcinogenic human liver fluke Opisthorchis viverrini , 2010, Proteomics.

[21]  Chris Oostenbrink,et al.  A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force‐field parameter sets 53A5 and 53A6 , 2004, J. Comput. Chem..

[22]  Xavier Robert,et al.  Deciphering key features in protein structures with the new ENDscript server , 2014, Nucleic Acids Res..

[23]  D S Goodsell,et al.  Automated docking of flexible ligands: Applications of autodock , 1996, Journal of molecular recognition : JMR.

[24]  The Fasciola hepatica thioredoxin: High resolution structure reveals two oxidation states. , 2008, Molecular and biochemical parasitology.

[25]  Janusz M. Bujnicki,et al.  MetaMQAP: A meta-server for the quality assessment of protein models , 2008, BMC Bioinformatics.

[26]  Identification and expression of Fasciola gigantica thioredoxin , 2014, Parasitology Research.

[27]  M. L. Jones,et al.  PDBsum: a Web-based database of summaries and analyses of all PDB structures. , 1997, Trends in biochemical sciences.

[28]  Ankita Gupta,et al.  Purification and characterization of two-domain glutaredoxin in the parasitic helminth Fasciola gigantica. , 2017, Parasitology international.

[29]  D. Eisenberg,et al.  Assessment of protein models with three-dimensional profiles , 1992, Nature.

[30]  J. Esteban,et al.  HUMAN FASCIOLASIS IN BOLIVIA: A GENERAL ANALYSIS AND A CRITICAL REVIEW OF EXISTING DATA , 1995 .

[31]  T. Beddoe,et al.  Current Threat of Triclabendazole Resistance in Fasciola hepatica. , 2016, Trends in parasitology.

[32]  T. Tripathi,et al.  Preferential regeneration of thioredoxin from parasitic flatworm Fasciola gigantica using glutathione system. , 2015, International journal of biological macromolecules.

[33]  Heather M. Alger,et al.  Molecular and enzymatic characterisation of Schistosoma mansoni thioredoxin. , 2002, International journal for parasitology.

[34]  A. Holmgren,et al.  Thioredoxins and glutaredoxins as facilitators of protein folding. , 2008, Biochimica et biophysica acta.

[35]  Meng Bian,et al.  Identification and immunological characterization of thioredoxin transmembrane-related protein from Clonorchis sinensis , 2013, Parasitology Research.

[36]  A. Loukas,et al.  Characterization of the antioxidant enzyme, thioredoxin peroxidase, from the carcinogenic human liver fluke, Opisthorchis viverrini. , 2008, Molecular and biochemical parasitology.

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

[38]  Peng Huang,et al.  Redox regulation of cell survival. , 2008, Antioxidants & redox signaling.

[39]  T. Tripathi,et al.  Unusual thiol-based redox metabolism of parasitic flukes. , 2017, Parasitology international.

[40]  F. Parra,et al.  Heterologous expression and functional characterization of thioredoxin from Fasciola hepatica , 2001, Parasitology Research.

[41]  A. W. Schüttelkopf,et al.  PRODRG: a tool for high-throughput crystallography of protein-ligand complexes. , 2004, Acta crystallographica. Section D, Biological crystallography.

[42]  Andrej Sali,et al.  Comparative Protein Structure Modeling Using MODELLER , 2014, Current protocols in bioinformatics.

[43]  D. Shi,et al.  Expression, characterization and crystal structure of thioredoxin from Schistosoma japonicum , 2015, Parasitology.

[44]  Nathan A. Baker,et al.  Electrostatics of nanosystems: Application to microtubules and the ribosome , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[45]  Chiu-Hsieh Hsu,et al.  The antitumor thioredoxin-1 inhibitor PX-12 (1-methylpropyl 2-imidazolyl disulfide) decreases thioredoxin-1 and VEGF levels in cancer patient plasma. , 2006, The Journal of laboratory and clinical medicine.

[46]  A. Loukas,et al.  Molecular expression and enzymatic characterization of thioredoxin from the carcinogenic human liver fluke Opisthorchis viverrini. , 2012, Parasitology international.

[47]  E. Novellino,et al.  Exploring the biological consequences of conformational changes in aspartame models containing constrained analogues of phenylalanine , 2016, Journal of enzyme inhibition and medicinal chemistry.

[48]  A. Holmgren,et al.  Thioredoxin and glutaredoxin systems. , 2019, The Journal of biological chemistry.

[49]  Claus Jacob,et al.  Metal and redox modulation of cysteine protein function. , 2003, Chemistry & biology.

[50]  Conrad C. Huang,et al.  UCSF Chimera—A visualization system for exploratory research and analysis , 2004, J. Comput. Chem..

[51]  B. Mannervik The enzymes of glutathione metabolism: an overview. , 1987, Biochemical Society transactions.

[52]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[53]  K. Becker,et al.  Thioredoxin and glutathione systems in Plasmodium falciparum. , 2012, International journal of medical microbiology : IJMM.

[54]  M. Ricaurte,et al.  Fasciola hepatica and Schistosoma mansoni: Identification of Common Proteins by Comparative Proteomic Analysis , 2011, The Journal of parasitology.

[55]  M. Wang,et al.  Fasciola hepatica , 2013, Endoscopy.

[56]  Banchob Sripa,et al.  Liver Fluke Induces Cholangiocarcinoma , 2007, PLoS medicine.

[57]  A. Holmgren,et al.  Glutaredoxin systems. , 2008, Biochimica et biophysica acta.

[58]  René Thomsen,et al.  MolDock: a new technique for high-accuracy molecular docking. , 2006, Journal of medicinal chemistry.

[59]  Peter M. Kasson,et al.  GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit , 2013, Bioinform..

[60]  Arthur J. Olson,et al.  AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading , 2009, J. Comput. Chem..

[61]  Elias S. J. Arnér,et al.  Reactive oxygen species, antioxidants, and the mammalian thioredoxin system. , 2001, Free radical biology & medicine.

[62]  Jie Shen,et al.  admetSAR: A Comprehensive Source and Free Tool for Assessment of Chemical ADMET Properties , 2012, J. Chem. Inf. Model..

[63]  C. Coustau,et al.  Excretory-secretory products of larval Fasciola hepatica investigated using a two-dimensional proteomic approach. , 2008, Molecular and biochemical parasitology.

[64]  M. Gilson,et al.  Calculation of protein-ligand binding affinities. , 2007, Annual review of biophysics and biomolecular structure.

[65]  A. Miele,et al.  Structural and functional characterization of Schistosoma mansoni Thioredoxin , 2011, Protein science : a publication of the Protein Society.

[66]  A. Holmgren,et al.  Thioredoxin structure and mechanism: conformational changes on oxidation of the active-site sulfhydryls to a disulfide. , 1995, Structure.

[67]  H. Eklund,et al.  Three-dimensional structure of Escherichia coli thioredoxin-S2 to 2.8 A resolution. , 1975, Proceedings of the National Academy of Sciences of the United States of America.

[68]  E. Herrero,et al.  Monothiol glutaredoxins: a common domain for multiple functions , 2007, Cellular and Molecular Life Sciences.

[69]  J. Gutteridge,et al.  Thiols in cellular redox signalling and control. , 2001, Current medicinal chemistry.

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

[71]  J. Williams Fascioliasis , 1961 .

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

[73]  J. Barrett,et al.  Proteomic analysis of Fasciola hepatica excretory‐secretory products , 2001, Proteomics.

[74]  Ben M. Webb,et al.  Comparative Protein Structure Modeling Using MODELLER , 2016, Current protocols in bioinformatics.

[75]  R. Gillies,et al.  The thioredoxin-1 inhibitor 1-methylpropyl 2-imidazolyl disulfide (PX-12) decreases vascular permeability in tumor xenografts monitored by dynamic contrast enhanced magnetic resonance imaging. , 2005, Clinical cancer research : an official journal of the American Association for Cancer Research.

[76]  S. Mas‐Coma,et al.  Fascioliasis and other plant-borne trematode zoonoses. , 2005, International journal for parasitology.

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