1
|
Mishra V, Pathak AK, Bandyopadhyay T. Binding of human serum albumin with uranyl ion at various pH: an all atom molecular dynamics study. J Biomol Struct Dyn 2023; 41:7318-7328. [PMID: 36099177 DOI: 10.1080/07391102.2022.2120080] [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: 07/28/2022] [Accepted: 08/26/2022] [Indexed: 10/14/2022]
Abstract
Uranium is routinely handled in various stages of nuclear fuel cycle and its association with human serum albumin (HSA) has been reported in literature, however, their binding characteristics still remains obscure. The present study aims to understand interaction of uranium with HSA by employing all atom molecular dynamics simulation of the HSA-metal ion complex. His67, His247 and Asp249 residues constitute the major binding site of HSA, which capture the uranyl ion (UO22+). A total of six sets of initial coordinates are used for Zn2+-HSA and UO22+-HSA system at pH = 4, 7.4 and 9, respectively. Enhance sampling method, namely, well-tempered meta-dynamics (WT-MtD) is employed to study the binding and un-binding processes of UO22+ and Zn2+ ions. Potential of mean force (PMF) profiles are generated for all the six sets of complexes from the converged WT-MtD run. Various basins and barriers are observed along the (un)binding pathways. Hydrogen bond dynamics and short-range Coulomb interactions are evaluated from the equilibrium run at each basins and barriers for both the ions at all pH values. The binding of UO22+ ion with HSA is the result of the dynamical balance between UO22+-HSA and UO22+-water short range Coulomb interactions. Zn2+ ion interact more strongly than UO22+ at all pH through short range Coulomb interactions. PMF values further concludes that UO22+ cannot associate to the Zn2+ bound HSA protein but can be captured by free HSA at all pH values i.e. endosomal, alkaline and physiological pH.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Vijayakriti Mishra
- Radiation Safety Systems Division, Bhabha Atomic Research Centre, Mumbai, India
- Homi Bhabha National Institute, Mumbai, India
| | - Arup Kumar Pathak
- Homi Bhabha National Institute, Mumbai, India
- Chemistry Division, Bhabha Atomic Research Centre, Mumbai, India
| | - Tusar Bandyopadhyay
- Homi Bhabha National Institute, Mumbai, India
- Chemistry Division, Bhabha Atomic Research Centre, Mumbai, India
| |
Collapse
|
2
|
In silico drug design and molecular docking of novel amidophosphonates and sulfamidophosphonates as inhibitors of urokinase-type plasminogen activator. J INDIAN CHEM SOC 2022. [DOI: 10.1016/j.jics.2022.100650] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
3
|
Zorbaz T, Malinak D, Hofmanova T, Maraković N, Žunec S, Hrvat NM, Andrys R, Psotka M, Zandona A, Svobodova J, Prchal L, Fingler S, Katalinić M, Kovarik Z, Musilek K. Halogen substituents enhance oxime nucleophilicity for reactivation of cholinesterases inhibited by nerve agents. Eur J Med Chem 2022; 238:114377. [DOI: 10.1016/j.ejmech.2022.114377] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2022] [Revised: 03/28/2022] [Accepted: 04/09/2022] [Indexed: 11/03/2022]
|
4
|
Ngo V, Li H, MacKerell AD, Allen TW, Roux B, Noskov S. Polarization Effects in Water-Mediated Selective Cation Transport across a Narrow Transmembrane Channel. J Chem Theory Comput 2021; 17:1726-1741. [PMID: 33539082 DOI: 10.1021/acs.jctc.0c00968] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Despite the progress in modeling complex molecular systems of ever-increasing complexity, a quantitatively accurate computational treatment of ion permeation through narrow membrane channels remains challenging. An important factor to reach this goal is induced electronic polarization, which is likely to impact the permeation rate of small ions through narrow molecular pores. In this work, we extended the recently developed polarizable force field based on the classical Drude oscillators to assess the role of induced polarization effects on the energetics of sodium and potassium ion transport across the gramicidin A (gA) ion channel. The inclusion of induced polarization lowers barriers present in 1D potential of mean force (PMF) for cation permeation by ∼50% compared to those obtained with the additive force field. Conductance properties calculated with 1D PMFs from Drude simulations are in better agreement with experimental results. Polarization of single-file water molecules and protein atoms forming the narrow pore has a direct impact on the free-energy barriers and cation-specific solid-state NMR chemical shifts. Sensitivity analysis indicates that small changes to water-channel interactions can alter the free energy barrier for ion permeation. These results, illustrating polarization effects present in the complex electrostatic environment of the gA channel, have broad implications for revising proposed mechanisms of ion permeation and selectivity in a variety of ion channels.
Collapse
Affiliation(s)
- Van Ngo
- Centre for Molecular Simulation, Department of Biological Sciences, University of Calgary, Calgary, Alberta T2N1N4, Canada.,Center for Nonlinear Studies, Los Alamos National Lab, Los Alamos, New Mexico 87544, United States
| | - Hui Li
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois 60637, United States
| | - Alexander D MacKerell
- Computer-Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland 21201, United States
| | - Toby W Allen
- School of Science, RMIT University, Melbourne, VIC 3001, Australia
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois 60637, United States
| | - Sergei Noskov
- Centre for Molecular Simulation, Department of Biological Sciences, University of Calgary, Calgary, Alberta T2N1N4, Canada
| |
Collapse
|
5
|
Mishra L, Sundararajan M, Bandyopadhyay T. MD simulation reveals differential binding of Cm(III) and Th(IV) with serum transferrin at acidic pH. Proteins 2020; 89:193-206. [PMID: 32892408 DOI: 10.1002/prot.26006] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2019] [Revised: 07/13/2020] [Accepted: 08/27/2020] [Indexed: 01/21/2023]
Abstract
The iron carrier human serum transferrin (sTf) is known to transport other metals, including some actinides (An). Radiotoxic An are routinely involved in the nuclear fuel cycle and the possibility of their accidental exposure cannot be ruled out. Understanding An interaction with sTf assumes a greater significance for the development of safe and efficacious chelators for their removal from the blood stream. Here we report several 100 ns equilibrium MD simulations of Cm(III)- and Th(IV)-loaded sTf at various protonation states of the protein to explore the possibility of the two An ions release and speciation. The results demonstrate variation in protonation state of dilysine pair (K206 and K296) and the tyrosine (Y188) residue is necessary for the opening of Cm(III)-bound protein and the release of the ion. For the tetravalent thorium, protonation of dilysine pair suffices to cause conformational changes of protein. However, in none of the protonation states, Th(IV) releases from sTf because of its strong electrostatic interaction with D63 in the first shell of the sTf binding cleft. Analysis of hydrogen bond, water bridge, and the evaluation of potential of mean forces of the An ions' release from sTf, substantiate the differential behavior of Cm(III) and Th(IV) at endosomal pH. The results provide insight in the regulation of Cm(III) and Th(IV) bioavailability that may prove useful for effective design of their decorporating agents and as well may help the future design of radiotherapy based on tetravalent ions.
Collapse
Affiliation(s)
- Lokpati Mishra
- Radiation Safety Systems Division, Bhabha Atomic Research Centre, Mumbai, India.,Department of Chemical Sciences, Homi Bhabha National Institute, Mumbai, India
| | | | - Tusar Bandyopadhyay
- Department of Chemical Sciences, Homi Bhabha National Institute, Mumbai, India.,Theoretical Chemistry Section, Bhabha Atomic Research Centre, Mumbai, India
| |
Collapse
|
6
|
Mishra L, Sundararajan M, Bandyopadhyay T. Molecular dynamics simulations of plutonium binding and its decorporation from the binding-cleft of human serum transferrin. J Biol Inorg Chem 2020; 25:213-231. [PMID: 31980924 DOI: 10.1007/s00775-020-01753-8] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2019] [Accepted: 01/03/2020] [Indexed: 11/29/2022]
Abstract
The possibility of plutonium (Pu) intake by radiation workers can not be ruled out. Transportation of Pu(IV) to various organs/cells is mainly carried through iron-carrying protein, serum transferrin (sTf), by receptor-mediated endocytosis. Understanding the Pu-sTf interaction is a primary step toward future design of its decorporating agents. We report MD simulations of Pu(IV) binding with sTf and look out for its decorporation at extracellular pH using suitable ligands. MD simulations were carried out in polarizable water environment at different protonation states of the protein. Results unravel the binding motif of Pu(IV): (1) sTf binds the ion in closed conformation at extracellular serum pH with carbonate as synergistic anions, (2) change in protonation state of dilysine (K206 and K296)-trigger and that of the carbonate ion at acidic endosomal pH is found to cause conformational changes of protein, conducive for the heavy ion to be released, although; (3) strong electrostatic interaction between D63 in the binding-cleft and Pu(IV) is found not to ever set free the ion. In an endeavour to decorporate Pu(IV), fragmented molecular form of hydroxypyridinone (HOPO) and catechol (CAM)-based ligands are docked at the binding site (BS) of the protein and metadynamics simulations are conducted. Pu(IV) binding at BS is found to be so strong that it was not detached from BS with the docked HOPO. However, for the identical set of simulation parameters, CAM is found to facilitate dislodging the heavy ion from the protein's binding influence. Differential behaviour of the two chelators is further explored. Fragmented molecular form of hydroxy-pyridinone (HOPO) and catecholamide (CAM) ligands were docked at the binding-site (BS) of human serum transferrin (sTf) to explore their feasibility as plausible Pu(IV) decorporating agents by employing metadynamics method. CAM was found to dislodge Pu from the sTf BS, while HOPO could not.
Collapse
Affiliation(s)
- Lokpati Mishra
- Radiation Safety Systems Division, Bhabha Atomic Research Centre, Mumbai, 400 085, India.,Homi Bhabha National Institute, Mumbai, 400094, India
| | - Mahesh Sundararajan
- Theoretical Chemistry Section, Bhabha Atomic Research Centre, Mumbai, 400 085, India
| | - Tusar Bandyopadhyay
- Homi Bhabha National Institute, Mumbai, 400094, India. .,Theoretical Chemistry Section, Bhabha Atomic Research Centre, Mumbai, 400 085, India.
| |
Collapse
|
7
|
Influence of gauche effect on uncharged oxime reactivators for the reactivation of tabun-inhibited AChE: quantum chemical and steered molecular dynamics studies. J Comput Aided Mol Des 2018; 32:793-807. [DOI: 10.1007/s10822-018-0130-1] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2018] [Accepted: 07/04/2018] [Indexed: 02/06/2023]
|
8
|
Pathak AK. Effect of pH on the hinge region of influenza viral protein: a combined constant pH and well-tempered molecular dynamics study. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:195101. [PMID: 29578453 DOI: 10.1088/1361-648x/aab98c] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Despite the knowledge that the influenza protein, hemagglutinin, undergoes a large conformational change at low pH during the process of fusion with the host cell, its molecular mechanism remains elusive. The present constant pH molecular dynamics (CpHMD) study identifies the residues responsible for large conformational change in acidic condition. Based on the pKa calculations, it is predicted that His-106 is much more responsible for the large conformational change than any other residues in the hinge region of hemagglutinin protein. Potential of mean force profile from well-tempered meta-dynamics (WT-MtD) simulation is also generated along the folding pathway by considering radius of gyration (R gyr) as a collective variable (CV). It is very clear from the present WT-MtD study, that the initial bending starts at that hinge region, which may trigger other conformational changes. Both the protein-protein and protein-water HB time correlation functions are monitored along the folding pathway. The protein-protein (full or hinge region) HB time correlation functions are always found to be stronger than those of the protein-water time correlation functions. The dynamical balance between protein-protein and protein-water HB interactions favors the stabilization of the folded state.
Collapse
Affiliation(s)
- Arup Kumar Pathak
- Theoretical Chemistry Section, Bhabha Atomic Research Centre, Mumbai 400085, India
| |
Collapse
|
9
|
Pathak AK, Bandyopadhyay T. Dynamic Mechanism of a Fluorinated Oxime Reactivator Unbinding from AChE Gorge in Polarizable Water. J Phys Chem B 2018; 122:3876-3888. [DOI: 10.1021/acs.jpcb.8b01171] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Arup K. Pathak
- Theoretical Chemistry Section, Bhabha Atomic Research Centre, Mumbai 400 085, India
| | - Tusar Bandyopadhyay
- Theoretical Chemistry Section, Bhabha Atomic Research Centre, Mumbai 400 085, India
| |
Collapse
|
10
|
Duan J, Hu C, Guo J, Guo L, Sun J, Zhao Z. A molecular dynamics study of the complete binding process of meropenem to New Delhi metallo-β-lactamase 1. Phys Chem Chem Phys 2018; 20:6409-6420. [PMID: 29442101 DOI: 10.1039/c7cp07459j] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The mechanism of substrate hydrolysis of New Delhi metallo-β-lactamase 1 (NDM-1) has been reported, but the process in which NDM-1 captures and transports the substrate into its active center remains unknown. In this study, we investigated the process of the substrate entry into the NDM-1 activity center through long unguided molecular dynamics simulations using meropenem as the substrate. A total of 550 individual simulations were performed, each of which for 200 ns, and 110 of them showed enzyme-substrate binding events. The results reveal three categories of relatively persistent and noteworthy enzyme-substrate binding configurations, which we call configurations A, B, and C. We performed binding free energy calculations of the enzyme-substrate complexes of different configurations using the molecular mechanics Poisson-Boltzmann surface area method. The role of each residue of the active site in binding the substrate was investigated using energy decomposition analysis. The simulated trajectories provide a continuous atomic-level view of the entire binding process, revealing potentially valuable regions where the enzyme and the substrate interact persistently and five possible pathways of the substrate entering into the active center, which were validated using well-tempered metadynamics. These findings provide important insights into the binding mechanism of meropenem to NDM-1, which may provide new prospects for the design of novel metallo-β-lactamase inhibitors and enzyme-resistant antibiotics.
Collapse
Affiliation(s)
- Juan Duan
- Department of Microbiology and Immunology of Guangdong Medical University, No. 2 West Wenming Road, Zhanjiang City, Guangdong Province 524023, China.
| | | | | | | | | | | |
Collapse
|
11
|
Xu Y, Cheng S, Sussman JL, Silman I, Jiang H. Computational Studies on Acetylcholinesterases. Molecules 2017; 22:molecules22081324. [PMID: 28796192 PMCID: PMC6152020 DOI: 10.3390/molecules22081324] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2017] [Revised: 08/07/2017] [Accepted: 08/07/2017] [Indexed: 01/18/2023] Open
Abstract
Functions of biomolecules, in particular enzymes, are usually modulated by structural fluctuations. This is especially the case in a gated diffusion-controlled reaction catalyzed by an enzyme such as acetylcholinesterase. The catalytic triad of acetylcholinesterase is located at the bottom of a long and narrow gorge, but it catalyzes the extremely rapid hydrolysis of the neurotransmitter, acetylcholine, with a reaction rate close to the diffusion-controlled limit. Computational modeling and simulation have produced considerable advances in exploring the dynamical and conformational properties of biomolecules, not only aiding in interpreting the experimental data, but also providing insights into the internal motions of the biomolecule at the atomic level. Given the remarkably high catalytic efficiency and the importance of acetylcholinesterase in drug development, great efforts have been made to understand the dynamics associated with its functions by use of various computational methods. Here, we present a comprehensive overview of recent computational studies on acetylcholinesterase, expanding our views of the enzyme from a microstate of a single structure to conformational ensembles, strengthening our understanding of the integration of structure, dynamics and function associated with the enzyme, and promoting the structure-based and/or mechanism-based design of new inhibitors for it.
Collapse
Affiliation(s)
- Yechun Xu
- CAS Key Laboratory of Receptor Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences (CAS), Shanghai 201203, China.
| | - Shanmei Cheng
- CAS Key Laboratory of Receptor Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences (CAS), Shanghai 201203, China.
| | - Joel L Sussman
- Israel Structural Proteomics Center, Weizmann Institute of Science, Rehovot 76100, Israel.
- Department of Structural Biology, Weizmann Institute of Science, Rehovot 76100, Israel.
| | - Israel Silman
- Department of Neurobiology, Weizmann Institute of Science, Rehovot 76100, Israel.
| | - Hualiang Jiang
- CAS Key Laboratory of Receptor Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences (CAS), Shanghai 201203, China.
| |
Collapse
|