1
|
Liggri PGV, Pérez-Garrido A, Tsitsanou KE, Dileep KV, Michaelakis A, Papachristos DP, Pérez-Sánchez H, Zographos SE. 2D finger-printing and molecular docking studies identified potent mosquito repellents targeting odorant binding protein 1. INSECT BIOCHEMISTRY AND MOLECULAR BIOLOGY 2023:103961. [PMID: 37217081 DOI: 10.1016/j.ibmb.2023.103961] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/08/2022] [Revised: 04/27/2023] [Accepted: 05/10/2023] [Indexed: 05/24/2023]
Abstract
Personal protection measures against the mosquitoes like the use of repellents constitute valuable tools in the effort to prevent the transmission of vector-borne diseases. Therefore, the discovery of novel repellent molecules which will be effective at lower concentrations and provide a longer duration of protection remains an urgent need. Mosquito Odorant-Binding Proteins (OBPs) involved in the initial steps of the olfactory signal transduction cascade have been recognized not only as passive carriers of odors and pheromones but also as the first molecular filter to discriminate semiochemicals, hence serving as molecular targets for the design of novel pest control agents. Among the three-dimensional structures of mosquito OBPs solved in the last decades, the OBP1 complexes with known repellents have been widely used as reference structures in docking analysis and molecular dynamics simulation studies for the structure-based discovery of new molecules with repellent activity. Herein, ten compounds known to be active against mosquitoes and/or displaying a binding affinity for Anopheles gambiae AgamOBP1 were used as queries in an in silico screening of over 96 million chemical samples in order to detect molecules with structural similarity. Further filtering of the acquired hits on the basis of toxicity, vapor pressure, and commercial availability resulted in 120 unique molecules that were subjected to molecular docking studies against OBP1. For seventeen potential OBP1-binders, the free energy of binding (FEB) and mode of interaction with the protein were further estimated by molecular docking simulations leading to the selection of eight molecules exhibiting the highest similarity with their parental compounds and favorable energy values. The in vitro determination of their binding affinity to AgamOBP1 and the evaluation of their repellent activity against female Aedes albopictus mosquitoes revealed that our combined ligand similarity screening and OBP1 structure-based molecular docking successfully detected three molecules with enhanced repellent properties. A novel DEET-like repellent with lower volatility (8.55 × 10-4 mmHg) but a higher binding affinity for OBP1 than DEET (1.35 × 10-3 mmHg). A highly active repellent molecule that is predicted to bind to the secondary Icaridin (sIC)-binding site of OBP1 with higher affinity than to the DEET-site and, therefore, represents a new scaffold to be exploited for the discovery of binders targeting multiple OBP sites. Finally, a third potent repellent exhibiting a high degree of volatility was found to be a strong DEET-site binder of OBP1 that could be used in slow-release formulations.
Collapse
Affiliation(s)
- Panagiota G V Liggri
- Institute of Chemical Biology, National Hellenic Research Foundation, 48 Vassileos Constantinou Avenue, 11635, Athens, Greece; Department of Biochemistry and Biotechnology, University of Thessaly, Biopolis, 41500, Larissa, Greece.
| | - Alfonso Pérez-Garrido
- Structural Bioinformatics and High Performance Computing Research Group (BIO-HPC), Universidad Católica de Murcia (UCAM), 30107, Spain
| | - Katerina E Tsitsanou
- Institute of Chemical Biology, National Hellenic Research Foundation, 48 Vassileos Constantinou Avenue, 11635, Athens, Greece
| | - Kalarickal V Dileep
- Laboratory for Computational and Structural Biology, Jubilee Centre for Medical Research, Jubilee Mission Medical College and Research Institute, Thrissur, Kerala, 680005, India
| | - Antonios Michaelakis
- Benaki Phytopathological Institute, Department of Entomology and Agricultural Zoology, 8 S Delta Str. 14561, Kifissia, Athens, Greece
| | - Dimitrios P Papachristos
- Benaki Phytopathological Institute, Department of Entomology and Agricultural Zoology, 8 S Delta Str. 14561, Kifissia, Athens, Greece
| | - Horacio Pérez-Sánchez
- Structural Bioinformatics and High Performance Computing Research Group (BIO-HPC), Universidad Católica de Murcia (UCAM), 30107, Spain.
| | - Spyros E Zographos
- Institute of Chemical Biology, National Hellenic Research Foundation, 48 Vassileos Constantinou Avenue, 11635, Athens, Greece.
| |
Collapse
|
2
|
Bairagya HR, Tasneem A, Rai GP, Reyaz S. New biochemical insights into the dynamics of water molecules at the GMP or IMP binding site of human GMPR enzyme: A molecular dynamics study. Proteins 2022; 90:200-217. [PMID: 34368983 DOI: 10.1002/prot.26207] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2020] [Revised: 06/24/2021] [Accepted: 07/22/2021] [Indexed: 12/31/2022]
Abstract
Human GMP reductase (hGMPR) enzyme is involved in a cellular metabolic pathway, converting GMP into IMP, and also it is an important target for anti-leukemic agents. Present computational investigations explain dynamical behavior of water molecules during the conformational transition process from GMP to IMP using molecular dynamics simulations. Residues at substrate-binding site of cancerous protein (PDB Id. 2C6Q) are mostly more dynamic in nature than the normal protein (PDB Id. 2BLE). Nineteen conserved water molecules are identified at the GMP/IMP binding site and are classified as (i) conserved stable dynamic and (ii) infrequent dynamic. Water molecules W11, W14, and W16 are classified as conserved stable dynamic due to their immobile character, whereas remaining water molecules (W1, W2, W3, W4, W5, W7, W8, W9, W10, W12, W13, W15, W17, W18, and W19) are infrequent with dynamic nature. Entrance or displacement of these infrequent water molecules at GMP/IMP sites may occur due to forward and backward movement of reference residues involving ligands. Four water molecules of hGMPR-I and nine water molecules of hGMPR-II are observed in repetitive transitions from GMP to IMP pathway, which indicates discrimination between two isoforms of hGMPRs. Water molecules in cancerous protein are more dynamic and unstable compared to normal protein. These water molecules execute rare dynamical events at GMP binding site and could assist in detailed understanding of conformational transitions that influence the hGMPR's biological functionality. The present study should be of interest to the experimental community engaged in leukemia research and drug discovery for CML cancer.
Collapse
Affiliation(s)
| | - Alvea Tasneem
- Department of Computer Science, Jamia Millia Islamia, New Delhi, India
| | - Gyan Prakash Rai
- Department of Computer Science, Jamia Millia Islamia, New Delhi, India
| | - Saima Reyaz
- Department of Computer Science, Jamia Millia Islamia, New Delhi, India
| |
Collapse
|
3
|
Lanrezac A, Férey N, Baaden M. Wielding the power of interactive molecular simulations. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1594] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- André Lanrezac
- CNRS, Laboratoire de Biochimie Théorique Université de Paris Paris France
| | - Nicolas Férey
- CNRS, Laboratoire interdisciplinaire des sciences du numérique Université Paris‐Saclay Orsay France
| | - Marc Baaden
- CNRS, Laboratoire de Biochimie Théorique Université de Paris Paris France
| |
Collapse
|
4
|
Reyaz S, Tasneem A, Rai GP, Bairagya HR. Investigation of structural analogs of hydroxychloroquine for SARS-CoV-2 main protease (Mpro): A computational drug discovery study. J Mol Graph Model 2021; 109:108021. [PMID: 34537554 PMCID: PMC8426616 DOI: 10.1016/j.jmgm.2021.108021] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2021] [Revised: 08/18/2021] [Accepted: 09/02/2021] [Indexed: 12/12/2022]
Abstract
The main protease (Mpro) is the key enzyme of nCOVID-19 and plays a decisive role that makes it an attractive drug target. Multiple analysis of crystal structures reveals the presence of W1, W2, and W3 water locations in the active site pocket of Mpro; W1 and W2 are unstable and are weakly bonded with protein in comparison to W3 of Mpro-native. So, we adopt the water displacement method to occupy W1 or W2 sites by triggering HCQ or its analogs to inactivate the enzyme. Virtual screening is employed to find out best analogs of HCQ, molecular docking is used for water displacement from catalytic region of Mpro, and finally, MD simulations are conducted for validation of these findings. The docking study reveals that W1 and W2 are occupied by respective atoms of ZINC28706440 whereas W2 by HCQ and indacaterol. Finally, MD results demonstrate (i) HCQ occupies W1 and W2 positions, but its analogs (indacaterol and ZINC28706440) are inadequate to retain either W1 or W2 (ii) His41 and Asp187 are stabilized by W3 in Mpro-native and His41, Cys145 and HCQ by W7 in ZINC28706440, and W4, W5, and W6 make water mediated bridge between indacaterol with His41. The structural, dynamical, and thermodynamic (WFP and J value) profiling parameters suggest that W3, W4, and W7 are prominent in their corresponding positions in comparison with W5 and W6. The final results conclude that ZINC28706440 may act as a best analog of HCQ with acceptable physico-chemical and toxicological scores and may further be synthesized for experimental validation.
Collapse
Affiliation(s)
- Saima Reyaz
- Department of Computer Science, Jamia Millia Islamia, New Delhi, 110025, India
| | - Alvea Tasneem
- Department of Computer Science, Jamia Millia Islamia, New Delhi, 110025, India
| | - Gyan Prakash Rai
- Department of Computer Science, Jamia Millia Islamia, New Delhi, 110025, India
| | - Hridoy R Bairagya
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, 110029, India.
| |
Collapse
|
5
|
Mukhopadhyay BP. Insights from molecular dynamics simulation of human ceruloplasmin (ferroxidase enzyme) binding with biogenic monoamines. Bioinformation 2019; 15:750-759. [PMID: 31831958 PMCID: PMC6900326 DOI: 10.6026/97320630015750] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2019] [Revised: 10/31/2019] [Accepted: 10/31/2019] [Indexed: 11/23/2022] Open
Abstract
Human ceruloplasmin (hCP) is a multi-copper oxidase with ferroxidase and amine oxidase activities. Molecular dynamics simulation (MDS) and docking analysis of biogenic monoamines with ceruloplasmin explain the role of Asp1025, Glu935, Glu272, Glu232 and Glu230 together with the binding site water molecules (referred as conserved water molecules) in the stabilization of neurotransmitter (Serotonin, Norepinephrine and Epinephrine) molecules within the binding cavity of hCP. Conserved water molecules are found at specific positions interacting with the protein structures that have sequence similarity. The ethylamine side chain nitrogen atom (N1) of neurotransmitter molecules interacts with water molecules in the binding cavity formed by Asp1025, Glu935 and Glu232 residues. These residues form an acidic triad mimicking a substrate binding cavity. The hydroxyl groups attached to the catechol ring of epinephrine and norepinephrine have been stabilized by Asp230 and Asp232 residues. Data suggests that the recognition of biogenic amines mediates through the N+(amine) ...Asp1025-His1026-CuCis-His path. The potential recognition path of biogenic monoamines to trinuclear copper cluster supported by active site water molecules (referred as conserved water molecules) is described in this report.
Collapse
|
6
|
Banerjee A, Bairagya HR, Mukhopadhyay BP, Nandi TK, Mishra DK. Conserved water mediated H-bonding dynamics of Ser117 and Thr119 residues in human transthyretin–thyroxin complexation: Inhibitor modeling study through docking and molecular dynamics simulation. J Mol Graph Model 2013; 44:70-80. [DOI: 10.1016/j.jmgm.2013.04.010] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2012] [Revised: 03/22/2013] [Accepted: 04/29/2013] [Indexed: 11/30/2022]
|
7
|
Chakrabarti B, Bairagya HR, Mishra DK, Chatterjee PK, Mukhopadhyay BP. Insight towards the conserved water mediated recognition of catalytic and structural Zn(+2) ions in human Matrix Metalloproteinase-8 enzyme: A study by MD-simulation methods. Bioinformation 2013; 9:126-33. [PMID: 23423544 PMCID: PMC3569599 DOI: 10.6026/97320630009126] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2013] [Accepted: 01/10/2013] [Indexed: 11/23/2022] Open
Abstract
Human matrix metalloproteinase-8 (hMMP-8) plays a important role in the progression of colorectal cancer, metastasis, multiple sclerosis and rheumetoid arthritis. Extensive MD-simulation of the PDB and solvated structures of hMMP-8 has revealed the presence of few conserved water molecules around the catalytic and structural zinc (ZnC and ZnS) ions. The coordination of two conserved water molecules (W and WS) to ZnS and the H-bonding interaction of WS to S151 have indicated the plausible involvement of that metal ion in the catalytic process. Beside this the coupling of ZnC and ZnS metal ions (ZnC - W(H) (W(1))…..W(2) ….H(162) - ZnS) through two conserved hydrophilic centers (occupied by water molecules) may also provide some rational on the recognition of two zinc ions which were separated by ~13 Å in their X-ray structures. This unique recognition of both the Zn(+2) ions in the enzyme through conserved water molecules may be implemented/ exploited for the design of antiproteolytic agent using water mimic drug design protocol.
Collapse
Affiliation(s)
- Bornali Chakrabarti
- Department of Chemistry, National Institute of Technology, Durgapur – 713 209, W.B. India
| | - Hridoy R Bairagya
- Department of Chemistry, National Institute of Technology, Durgapur – 713 209, W.B. India
| | - Deepak Kr Mishra
- Department of Chemistry, National Institute of Technology, Durgapur – 713 209, W.B. India
| | - Pradip Kumar Chatterjee
- Thermal Engineering Department, Central Mechanical Engineering Research Institute, Durgapur – 713209,W.B, India
| | - Bishnu P Mukhopadhyay
- Department of Chemistry, National Institute of Technology, Durgapur – 713 209, W.B. India
| |
Collapse
|
8
|
Nandi TK, Bairagya HR, Mukhopadhyay BP, Mallik P, Sukul D, Bera AK. Conserved water-mediated H-bonding dynamics of catalytic His159 and Asp158: insight into a possible acid–base coupled mechanism in plant thiol protease. J Mol Model 2011; 18:2633-44. [DOI: 10.1007/s00894-011-1277-z] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2010] [Accepted: 10/09/2011] [Indexed: 10/15/2022]
|
9
|
Bairagya HR, Mukhopadhyay BP, Bhattacharya S. Role of the conserved water molecules in the binding of inhibitor to IMPDH-II (human): A study on the water mimic inhibitor design. ACTA ACUST UNITED AC 2009. [DOI: 10.1016/j.theochem.2009.04.037] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|