1
|
Bairagya HR. Dynamics of nucleoplasm in human leukemia cells: A thrust towards designing anti-leukemic agents. J Mol Graph Model 2024; 131:108807. [PMID: 38908255 DOI: 10.1016/j.jmgm.2024.108807] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2024] [Revised: 04/20/2024] [Accepted: 06/02/2024] [Indexed: 06/24/2024]
Abstract
The human inosine monophosphate dehydrogenase (hIMPDH) is a metabolic enzyme that possesses a unique ability to self-assemble into higher-order structures, forming cytoophidia. The hIMPDH II isoform is more active in chronic myeloid leukemia (CML) cancer cells, making it a promising target for anti-leukemic therapy. However, the structural details and molecular mechanisms of the dynamics of hIMPDHcytoophidia assembly in vitro need to be better understood, and it is crucial to reconstitute the computational nucleoplasm model with cytophilic-like polymers in vitro to characterize their structure and function. Finally, a computational model and its dynamics of the nucleoplasm for CML cells have been proposed in this short review. This research on nucleoplasm aims to aid the scientific community's understanding of how metabolic enzymes like hIMPDH function in cancer and normal cells. However, validating and justifying the computational results from modeling and simulation with experimental data is essential. The new insights gained from this research could explain the structure/topology, geometrical, and electronic consequences of hIMPDH inhibitors on leukemic and normal cells. They could lead to further advancements in the knowledge of nucleoplasmic chemical reaction dynamics.
Collapse
Affiliation(s)
- Hridoy R Bairagya
- Computational Drug Design and Bio-molecular Simulation Lab, Department of Bioinformatics, Maulana Abul Kalam Azad University of Technology, West Bengal, 741249, India.
| |
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.5] [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
|
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: 7] [Impact Index Per Article: 2.3] [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
|
4
|
Bairagya HR, Tasneem A, Rai GP, Reyaz S. Structural and Dynamical Impact of Water Molecules at Substrate- or Product-Binding Sites in Human GMPR Enzyme: A Study by Molecular Dynamics Simulations. J Phys Chem B 2021; 125:1351-1362. [PMID: 33369428 DOI: 10.1021/acs.jpcb.0c08818] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Human guanosine monophosphate reductase (hGMPR) enzyme maintains the intracellular balance between adenine and guanine nucleotide pools, and it is an excellent target for the design of isoform-specific antileukemic agents. In the present study, we have investigated solvation properties of substrate GMP or product inosine-5'-monophosphate (IMP)-binding pocket of hGMPR by employing molecular dynamics simulations on conformations A (substrate GMP), B [substrate GMP with cofactor nicotinamide adenine dinucleotide phosphate (NDP)], C (product IMP with cofactor NDP), and D (product IMP). Nineteen water sites are identified precisely; they are responsible for the catalytic activity of this site, control structural and dynamical integrity, and electronic consequences of GMP or IMP in the binding site of hGMPR. The water sites of category-1 (W1, W4, W5, W6, W13, and W15) in normal protein and category-2 (W2, W3, W7, W8, W10, W17, and W18) in cancerous protein are unique and stabilize the guanosine or inosine group of GMP or IMP for participation in the enzymatic reaction, whereas the remaining water centers either stabilize pentose sugar ribose or the phosphate group of GMP or IMP. Furthermore, water sites of category-4 (W11, W14, and W16) appear to be conserved in all conformations during the entire simulation. The GMP-binding site in cancerous protein 2C6Q is significantly expanded, and its dynamics are very different from normal protein 2BLE. Furthermore, unique interactions of GMP(N1)···W2···Asp129/Asn158, IMP(N1)···W3···Glu289, and IMP(O6)···W10···Ser270 might be used in a water mimic drug design for hGMPR-II. In this context, water finding probability, relative interaction energy (J) associated with water site W, entropy, and topologies of these three water sites are thermodynamically acceptable for the water displacement method by the modified ligand. Hence, their positions in the catalytic pocket may also facilitate future drug discovery for chronic myelogenous leukemia by the design of appropriately oriented chemical groups that may displace these water molecules to mimic their structural, electronic, and thermodynamic properties.
Collapse
Affiliation(s)
- Hridoy R Bairagya
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K
| | - Alvea Tasneem
- Department of Computer Science, Jamia Millia Islamia, Jamia Nagar, Okhla, New Delhi 110025, India
| | - Gyan Prakash Rai
- Department of Computer Science, Jamia Millia Islamia, Jamia Nagar, Okhla, New Delhi 110025, India
| | - Saima Reyaz
- Department of Computer Science, Jamia Millia Islamia, Jamia Nagar, Okhla, New Delhi 110025, India
| |
Collapse
|
5
|
Mukhopadhyay BP. Putative role of conserved water molecules in the hydration and inter-domain recognition of mono nuclear copper centers in O2-bound human ceruloplasmin: A comparative study between X-ray and MD simulated structures. Bioinformation 2019; 15:402-411. [PMID: 31312077 PMCID: PMC6614124 DOI: 10.6026/97320630015402] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2019] [Accepted: 06/02/2019] [Indexed: 11/23/2022] Open
Abstract
Human Ceruloplasmin (hCP) is an unique multicopper oxidase which involves in different biological functions e.g., iron metabolism, copper transportation, biogenic amine oxidation ,and its malfunction causes Wilson's and Menkes diseases. MD- simulation studies of O2- bound solvated structure have revealed the role of several conserved/ semi-conserved water molecules in the hydration of type-I copper centers and their involvement to recognition dynamics of these metal centers. In O2- bound structure, hydration potentiality of CuRS (Cu1106) type-I copper center is observed to be unique, where two water molecules (W1-W2) are interacting with the metal sites, which was not found in X-ray structures of hCP. Generally, in the interdomain recognition of CuCys-His to CuRS, CuRS to CuPR and CuPR to CuCys-His centers, the copper bound His-residue of one domain interacts with the Glu-residue of other complementary domain through conserved/ semi-conserved (W3 to W5) water- mediated hydrogen bonds (Cu-His...W...Glu), however direct salt-bridge (Cu-His...Glu) interaction were observed in the X- ray structures. The MD- simulated and X- ray structures have indicated some possibilities on the Cu-His...W...Glu ↔ Cu-His...Glu transition during the interdomain recognition of type-I copper centers, which may have some importance in biology and chemistry of ceruloplasmin.
Collapse
Affiliation(s)
- Bishnu P. Mukhopadhyay
- Department of Chemistry National Institute of Technology-Durgapur, West Bengal, Durgapur - 713209, India
| |
Collapse
|
6
|
Putative role of invariant water molecules in the X-ray structures of family G fungal endoxylanases. J Biosci 2018. [DOI: 10.1007/s12038-018-9752-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
7
|
Dasgupta S, Mukherjee S, Mukhopadhyay BP, Banerjee A, Mishra DK. Recognition dynamics of dopamine to human Monoamine oxidase B: role of Leu171/Gln206 and conserved water molecules in the active site cavity. J Biomol Struct Dyn 2017; 36:1439-1462. [PMID: 28460566 DOI: 10.1080/07391102.2017.1325405] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Subrata Dasgupta
- Department of Chemistry, National Institute of Technology-Durgapur, Durgapur 713209, West Bengal, India
| | - Soumita Mukherjee
- Department of Chemistry, National Institute of Technology-Durgapur, Durgapur 713209, West Bengal, India
| | - Bishnu P Mukhopadhyay
- Department of Chemistry, National Institute of Technology-Durgapur, Durgapur 713209, West Bengal, India
| | - Avik Banerjee
- Department of Chemistry, National Institute of Technology-Durgapur, Durgapur 713209, West Bengal, India
| | - Deepak K Mishra
- Department of Chemistry, National Institute of Technology-Durgapur, Durgapur 713209, West Bengal, India
| |
Collapse
|
8
|
Chakrabarti B, Bairagya HR, Mukhopadhyay BP, Sekar K. New biochemical insight of conserved water molecules at catalytic and structural Zn 2+ ions in human matrix metalloproteinase-I: a study by MD-simulation. J Mol Model 2017; 23:57. [PMID: 28161785 DOI: 10.1007/s00894-017-3210-6] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2015] [Accepted: 01/09/2017] [Indexed: 12/01/2022]
Abstract
Human matrix metalloproteinase (MMP)-1 or collagenase-1 plays a significant role in embryonic development, tissue remodeling, and is also involved in several diseases like arthritis, metastasis, etc. Molecular dynamics simulation studies on hMMP-1 X-ray structures (PDB Id. 1CGE, 1CGF, 1CGL, 1HFC, and 2TCL) suggest that the three conserved water molecules (WH/1, WI, WS) are coordinated with catalytic zinc (ZnC), and one water molecule (W) is associated at structural zinc ion (ZnS). Transition of the coordination geometry around ZnC from tetrahedral to octahedral and tetrahedral to trigonal bipyramidal at ZnS are also observed during the dynamics. Recognition of two zinc ions through water mediated bridges (ZnC - WH (W1)…W2….H183 - ZnS) and stabilization of secondary coordination zone around the metal ions indicates the possibility of ZnC…ZnS coupled catalytic mechanism in hMMP-I. This study not only reveals a functionally important role of conserved water molecules in hMMP-I but also highlights the involvement of other non catalytic residues, such as S172 and D170 in the catalytic mechanism. The results obtained in this study could be relevant for importance of conserved water mediated recognition site of the sequence residue id. 202(RWTNNFREY)210, interaction of W(tryptophan)203 to zinc bound histidine, their influence on the water molecules that are involved in bridging between ZnC and ZnS, and structure-based design of specific hMMP inhibitors. Graphical abstract Water mediated recognition of structural and catalytic zinc ions of hMMP-1 structure (MD simulatated conformation).
Collapse
Affiliation(s)
- Bornali Chakrabarti
- Department of Chemistry, National Institute of Technology-Durgapur, Durgapur, 713209, India
| | - Hridoy R Bairagya
- Department of Chemistry, National Institute of Technology-Durgapur, Durgapur, 713209, India.,Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge, CB21EW, UK
| | - Bishnu P Mukhopadhyay
- Department of Chemistry, National Institute of Technology-Durgapur, Durgapur, 713209, India.
| | - K Sekar
- Department of Computational and Data Sciences Indian Institute of Science, Laboratory for Structural Biology and Bio-computing, Bangalore, 560 012, India
| |
Collapse
|
9
|
Bairagya HR, Bansal M. New insight into the architecture of oxy-anion pocket in unliganded conformation of GAT domains: A MD-simulation study. Proteins 2016; 84:360-73. [PMID: 26756917 DOI: 10.1002/prot.24983] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2015] [Revised: 12/06/2015] [Accepted: 12/19/2015] [Indexed: 11/11/2022]
Abstract
Human Guanine Monophosphate Synthetase (hGMPS) converts XMP to GMP, and acts as a bifunctional enzyme with N-terminal "glutaminase" (GAT) and C-terminal "synthetase" domain. The enzyme is identified as a potential target for anti-cancer and immunosuppressive therapies. GAT domain of enzyme plays central role in metabolism, and contains conserved catalytic residues Cys104, His190, and Glu192. MD simulation studies on GAT domain suggest that position of oxyanion in unliganded conformation is occupied by one conserved water molecule (W1), which also stabilizes that pocket. This position is occupied by a negatively charged atom of the substrate or ligand in ligand bound crystal structures. In fact, MD simulation study of Ser75 to Val indicates that W1 conserved water molecule is stabilized by Ser75, while Thr152, and His190 also act as anchor residues to maintain appropriate architecture of oxyanion pocket through water mediated H-bond interactions. Possibly, four conserved water molecules stabilize oxyanion hole in unliganded state, but they vacate these positions when the enzyme (hGMPS)-substrate complex is formed. Thus this study not only reveals functionally important role of conserved water molecules in GAT domain, but also highlights essential role of other non-catalytic residues such as Ser75 and Thr152 in this enzymatic domain. The results from this computational study could be of interest to experimental community and provide a testable hypothesis for experimental validation. Conserved sites of water molecules near and at oxyanion hole highlight structural importance of water molecules and suggest a rethink of the conventional definition of chemical geometry of inhibitor binding site.
Collapse
Affiliation(s)
- Hridoy R Bairagya
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore, Karnataka, 560012, India
| | - Manju Bansal
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore, Karnataka, 560012, India
| |
Collapse
|