1
|
Pitaloka DAE, Arfan A, Ramadhan DSF, Chaidir L. Insights from the molecular mechanism of pyrazinamide to mutated pyrazinamidase linked to the pncA gene in clinical isolates of Mycobacterium tuberculosis. J Biomol Struct Dyn 2024; 42:759-765. [PMID: 37096659 DOI: 10.1080/07391102.2023.2195002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2022] [Accepted: 03/18/2023] [Indexed: 04/26/2023]
Abstract
This study aims to conduct a comprehensive molecular dynamics strategy to evaluate whether mutations found in pyrazinamide monoresistant (PZAMR) strains of Mycobacterium tuberculosis (MTB) can potentially reduce the effectiveness of pyrazinamide (PZA) for tuberculosis (TB) treatment. Five single point mutations of pyrazinamidase (PZAse), an enzyme which is responsible for the activation of prodrug PZA into pyrazinoic acid, found in MTB clinical isolates, namely His82Arg, Thr87Met, Ser66Pro, Ala171Val, and Pro62Leu, were analyzed by the dynamics simulations both in the apo state (unbound state) and in the PZA bound state. The results showed that the mutation of His82 to Arg, Thr87 to Met, and Ser66 to Pro in PZAse affects the coordination state of the Fe2+ ion, which is a cofactor required for enzyme activity. These mutations change the flexibility, stability, and fluctuation of His51, His57, and ASP49 amino acid residues around the Fe2+ ion, culminating in an unstable complex and dissociation of PZA from the PZAse binding site. However, mutations of Ala171 to Val and Pro62 to Leu were found to have no effect on the complex's stability. Based on the results, PZAse mutations of His82Arg, Thr87Met, and Ser66Pro culminated in weak binding affinity for PZA and caused significant structural deformations that led to PZA resistance. Future structural and functional studies, as well as investigations into other aspects of drug resistance in PZAse, will require experimental clarification.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Dian Ayu Eka Pitaloka
- Department of Pharmacology and Clinical Pharmacy, Faculty of Pharmacy, Universitas Padjadjaran, Sumedang, Indonesia
- Center for Translational Biomarker Research, Universitas Padjadjaran, Bandung, Indonesia
| | - Arfan Arfan
- Department of Medicinal Chemistry, Faculty of Pharmacy, Universitas Halu Oleo, Kendari, Indonesia
| | - Dwi Syah Fitra Ramadhan
- Department of Pharmaceutical Chemistry, Sekolah Tinggi Ilmu Kesehatan Mandala Waluya, Kendari, Indonesia
| | - Lidya Chaidir
- Center for Translational Biomarker Research, Universitas Padjadjaran, Bandung, Indonesia
- Department of Biomedical Sciences, Faculty of Medicine, Universitas Padjadjaran, Sumedang, Indonesia
| |
Collapse
|
2
|
The Hydrolysis Rate of Paraoxonase-1 Q and R Isoenzymes: An In Silico Study Based on In Vitro Data. Molecules 2022; 27:molecules27206780. [PMID: 36296373 PMCID: PMC9607273 DOI: 10.3390/molecules27206780] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2022] [Revised: 09/30/2022] [Accepted: 10/05/2022] [Indexed: 11/17/2022] Open
Abstract
Human serum paraoxonase-1 (PON1) is an important hydrolase-type enzyme found in numerous tissues. Notably, it can exist in two isozyme-forms, Q and R, that exhibit different activities. This study presents an in silico (QSAR, Docking, MD and QM/MM) study of a set of compounds on the activity towards the PON1 isoenzymes (QPON1 and RPON1). Different rates of reaction for the Q and R isoenzymes were analyzed by modelling the effect of Q192R mutation on active sites. It was concluded that the Q192R mutation is not even close to the active site, while it is still changing the geometry of it. Using the combined genetic algorithm with multiple linear regression (GA-MLR) technique, several QSAR models were developed and relative activity rates of the isozymes of PON1 explained. From these, two QSAR models were selected, one each for the QPON1 and RPON1. Best selected models are four-variable MLR models for both Q and R isozymes with squared correlation coefficient R2 values of 0.87 and 0.83, respectively. In addition, the applicability domain of the models was analyzed based on the Williams plot. The results were discussed in the light of the main factors that influence the hydrolysis activity of the PON1 isozymes.
Collapse
|
3
|
Fakih TM, Kurniawan F, Yusuf M, Mudasir M, Tjahjono DH. Molecular Dynamics of Cobalt Protoporphyrin Antagonism of the Cancer Suppressor REV-ERBβ. Molecules 2021; 26:3251. [PMID: 34071361 PMCID: PMC8198987 DOI: 10.3390/molecules26113251] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2021] [Revised: 05/15/2021] [Accepted: 05/26/2021] [Indexed: 11/17/2022] Open
Abstract
Nuclear receptor REV-ERBβ is an overexpressed oncoprotein that has been used as a target for cancer treatment. The metal-complex nature of its ligand, iron protoporphyrin IX (Heme), enables the REV-ERBβ to be used for multiple therapeutic modalities as a photonuclease, a photosensitizer, or a fluorescence imaging agent. The replacement of iron with cobalt as the metal center of protoporphyrin IX changes the ligand from an agonist to an antagonist of REV-ERBβ. The mechanism behind that phenomenon is still unclear, despite the availability of crystal structures of REV-ERBβ in complex with Heme and cobalt protoporphyrin IX (CoPP). This study used molecular dynamic simulations to compare the effects of REV-ERBβ binding to Heme and CoPP, respectively. The initial poses of Heme and CoPP in complex with agonist and antagonist forms of REV-ERBβ were predicted using molecular docking. The binding energies of each ligand were calculated using the MM/PBSA method. The computed binding affinity of Heme to REV-ERBβ was stronger than that of CoPP, in agreement with experimental results. CoPP altered the conformation of the ligand-binding site of REV-ERBβ, disrupting the binding site for nuclear receptor corepressor, which is required for REV-ERBβ to regulate the transcription of downstream target genes. Those results suggest that a subtle change in the metal center of porphyrin can change the behavior of porphyrin in cancer cell signaling. Therefore, modification of porphyrin-based agents for cancer therapy should be conducted carefully to avoid triggering unfavorable effects.
Collapse
Affiliation(s)
- Taufik Muhammad Fakih
- School of Pharmacy, Bandung Institute of Technology, Jalan Ganesha 10, Bandung 40135, Indonesia; (T.M.F.); (F.K.)
- Department of Pharmacy, Faculty of Mathematics and Natural Sciences, Universitas Islam Bandung, Jalan Rangga Gading 8, Bandung 40116, Indonesia
| | - Fransiska Kurniawan
- School of Pharmacy, Bandung Institute of Technology, Jalan Ganesha 10, Bandung 40135, Indonesia; (T.M.F.); (F.K.)
| | - Muhammad Yusuf
- Department of Chemistry, Faculty of Mathematics and Natural Sciences, Universitas Padjadjaran, Jalan Raya Bandung Sumedang Km 21, Sumedang 45363, Indonesia;
| | - Mudasir Mudasir
- Department of Chemistry, Faculty of Mathematics and Natural Sciences, Universitas Gadjah Mada, Sekip Utara BLS 21, Yogyakarta 55281, Indonesia;
| | - Daryono Hadi Tjahjono
- School of Pharmacy, Bandung Institute of Technology, Jalan Ganesha 10, Bandung 40135, Indonesia; (T.M.F.); (F.K.)
| |
Collapse
|
4
|
Docking-Based Virtual Screening and Molecular Dynamics Simulations of Quercetin Analogs as Enoyl-Acyl Carrier Protein Reductase (InhA) Inhibitors of Mycobacterium tuberculosis. Sci Pharm 2021. [DOI: 10.3390/scipharm89020020] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
The emergence of multidrug-resistant Mycobacterium tuberculosis (MTB) has become a major problem in treating tuberculosis (TB) and shows the need to develop new and efficient drugs for better TB control. This study aimed to use in silico techniques to discover potential inhibitors to the Enoyl-[acyl-carrier-protein] reductase (InhA), which controls mycobacterial cell wall construction. Initially, 391 quercetin analogs present in the KNApSAck_3D database were selected, filters were sequentially applied by docking-based virtual screening. After recategorizing the variables (bond energy prediction and molecular interaction, including hydrogen bond and hydrophobic bond), compounds C00013874, C00006532, and C00013887 were selected as hit ligands. These compounds showed great hydrophobic contributions, and for each hit ligand, 100 ns of molecular dynamic simulations were performed, and the binding free energy was calculated. C00013874 demonstrated the greatest capacity for the InhA enzyme inhibition with ΔGbind = −148.651 kcal/mol compare to NAD (native ligand) presented a ΔGbind = −87.570 kcal/mol. These data are preliminary studies and might be a suitable candidate for further experimental analysis.
Collapse
|
5
|
Corrigan RA, Qi G, Thiel AC, Lynn JR, Walker BD, Casavant TL, Lagardere L, Piquemal JP, Ponder JW, Ren P, Schnieders MJ. Implicit Solvents for the Polarizable Atomic Multipole AMOEBA Force Field. J Chem Theory Comput 2021; 17:2323-2341. [PMID: 33769814 DOI: 10.1021/acs.jctc.0c01286] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Computational protein design, ab initio protein/RNA folding, and protein-ligand screening can be too computationally demanding for explicit treatment of solvent. For these applications, implicit solvent offers a compelling alternative, which we describe here for the polarizable atomic multipole AMOEBA force field based on three treatments of continuum electrostatics: numerical solutions to the nonlinear and linearized versions of the Poisson-Boltzmann equation (PBE), the domain-decomposition conductor-like screening model (ddCOSMO) approximation to the PBE, and the analytic generalized Kirkwood (GK) approximation. The continuum electrostatics models are combined with a nonpolar estimator based on novel cavitation and dispersion terms. Electrostatic model parameters are numerically optimized using a least-squares style target function based on a library of 103 small-molecule solvation free energy differences. Mean signed errors for the adaptive Poisson-Boltzmann solver (APBS), ddCOSMO, and GK models are 0.05, 0.00, and 0.00 kcal/mol, respectively, while the mean unsigned errors are 0.70, 0.63, and 0.58 kcal/mol, respectively. Validation of the electrostatic response of the resulting implicit solvents, which are available in the Tinker (or Tinker-HP), OpenMM, and Force Field X software packages, is based on comparisons to explicit solvent simulations for a series of proteins and nucleic acids. Overall, the emergence of performative implicit solvent models for polarizable force fields opens the door to their use for folding and design applications.
Collapse
Affiliation(s)
- Rae A Corrigan
- Roy J Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa 52242, United States
| | - Guowei Qi
- Department of Biochemistry, University of Iowa, Iowa City, Iowa 52242, United States
| | - Andrew C Thiel
- Roy J Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa 52242, United States
| | - Jack R Lynn
- Roy J Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa 52242, United States
| | - Brandon D Walker
- Department of Biomedical Engineering, University of Texas in Austin, Austin, Texas 78712, United States
| | - Thomas L Casavant
- Roy J Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa 52242, United States
| | - Louis Lagardere
- Department of Chemistry, Sorbonne Université, F-75005 Paris, France
| | | | - Jay W Ponder
- Department of Chemistry, Washington University in St. Louis, St. Louis, Missouri 63130, United States
| | - Pengyu Ren
- Department of Biomedical Engineering, University of Texas in Austin, Austin, Texas 78712, United States
| | - Michael J Schnieders
- Roy J Carver Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa 52242, United States.,Department of Biochemistry, University of Iowa, Iowa City, Iowa 52242, United States
| |
Collapse
|
6
|
Forouzesh N, Mukhopadhyay A, Watson LT, Onufriev AV. Multidimensional Global Optimization and Robustness Analysis in the Context of Protein-Ligand Binding. J Chem Theory Comput 2020; 16:4669-4684. [PMID: 32450041 PMCID: PMC8594251 DOI: 10.1021/acs.jctc.0c00142] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Accuracy of protein-ligand binding free energy calculations utilizing implicit solvent models is critically affected by parameters of the underlying dielectric boundary, specifically, the atomic and water probe radii. Here, a global multidimensional optimization pipeline is developed to find optimal atomic radii specifically for protein-ligand binding calculations in implicit solvent. The computational pipeline has these three key components: (1) a massively parallel implementation of a deterministic global optimization algorithm (VTDIRECT95), (2) an accurate yet reasonably fast generalized Born implicit solvent model (GBNSR6), and (3) a novel robustness metric that helps distinguish between nearly degenerate local minima via a postprocessing step of the optimization. A graph-based "kT-connectivity" approach to explore and visualize the multidimensional energy landscape is proposed: local minima that can be reached from the global minimum without exceeding a given energy threshold (kT) are considered to be connected. As an illustration of the capabilities of the optimization pipeline, we apply it to find a global optimum in the space of just five radii: four atomic (O, H, N, and C) radii and water probe radius. The optimized radii, ρW = 1.37 Å, ρC = 1.40 Å, ρH = 1.55 Å, ρN = 2.35 Å, and ρO = 1.28 Å, lead to a closer agreement of electrostatic binding free energies with the explicit solvent reference than two commonly used sets of radii previously optimized for small molecules. At the same time, the ability of the optimizer to find the global optimum reveals fundamental limits of the common two-dielectric implicit solvation model: the computed electrostatic binding free energies are still almost 4 kcal/mol away from the explicit solvent reference. The proposed computational approach opens the possibility to further improve the accuracy of practical computational protocols for binding free energy calculations.
Collapse
Affiliation(s)
- Negin Forouzesh
- Department of Computer Science, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
| | - Abhishek Mukhopadhyay
- Department of Physics, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
| | - Layne T Watson
- Department of Computer Science, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
- Department of Mathematics, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
- Department of Aerospace and Ocean Engineering, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
- Center for Soft Matter and Biological Physics, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
| | - Alexey V Onufriev
- Department of Computer Science, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
- Department of Physics, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
- Center for Soft Matter and Biological Physics, Virginia Polytechnic Institute & State University, Blacksburg, Virginia 24061, United States
| |
Collapse
|
7
|
Perez CP, Elmore DE, Radhakrishnan ML. Computationally Modeling Electrostatic Binding Energetics in a Crowded, Dynamic Environment: Physical Insights from a Peptide–DNA System. J Phys Chem B 2019; 123:10718-10734. [DOI: 10.1021/acs.jpcb.9b09478] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
|
8
|
Nunes R, Vila-Viçosa D, Costa PJ. Tackling Halogenated Species with PBSA: Effect of Emulating the σ-Hole. J Chem Theory Comput 2019; 15:4241-4251. [PMID: 31142112 DOI: 10.1021/acs.jctc.9b00106] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
To model halogen-bond phenomena using classical force fields, an extra point (EP) of charge is frequently introduced at a given distance from the halogen (X) to emulate the σ-hole. The resulting molecular dynamics (MD) trajectories can be used in subsequent molecular mechanics (MM) combined with Poisson-Boltzmann and surface area calculations (PBSA) to estimate protein-ligand binding free energies (Δ Gbind). While EP addition improves the MM/MD description of halogen-containing systems, its effect on the calculation of solvation free energies (Δ Gsolv) using the PBSA approach is yet to be assessed. As the PBSA calculations depend, among other parameters, on the empirical assignment of radii (PB radii), a problematic issue arises, since standard halogen radii are smaller than the typical X···EP distances, thus placing the EP within the solvent dielectric. Herein, we took a common literature EP parametrization scheme, which uses X···EP = Rmin and RESP charges in the context of GAFF, and performed a comprehensive study on the performance of PBSA (using three different setups) in the calculation of Δ Gsolv values for 142 halogenated compounds (bearing Cl, Br, or I) for which the experimental values are known. By conducting an optimization (minimizing the error against experimental values), we provide a new optimized set of halogen PB radii, for each PBSA setup, that should be used in the context of the aforementioned scenario. A simultaneous optimization of PB radii and X···EP distances shows that a wide range of distance/radius pairs can be used without significant loss of accuracy, therefore laying the basis for expanding this halogen radii optimization strategy to other force fields and EP implementations. As ligand Δ Gsolv estimation is an important term in the determination of protein-ligand Δ Gbind, this work is particularly relevant in the framework of structure-based virtual screening and related computer-aided drug design routines.
Collapse
|
9
|
Aleksandrov A, Lin FY, Roux B, MacKerell AD. Combining the polarizable Drude force field with a continuum electrostatic Poisson-Boltzmann implicit solvation model. J Comput Chem 2018; 39:1707-1719. [PMID: 29737546 DOI: 10.1002/jcc.25345] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2017] [Revised: 02/26/2018] [Accepted: 04/08/2018] [Indexed: 12/13/2022]
Abstract
In this work, we have combined the polarizable force field based on the classical Drude oscillator with a continuum Poisson-Boltzmann/solvent-accessible surface area (PB/SASA) model. In practice, the positions of the Drude particles experiencing the solvent reaction field arising from the fixed charges and induced polarization of the solute must be optimized in a self-consistent manner. Here, we parameterized the model to reproduce experimental solvation free energies of a set of small molecules. The model reproduces well-experimental solvation free energies of 70 molecules, yielding a root mean square difference of 0.8 kcal/mol versus 2.5 kcal/mol for the CHARMM36 additive force field. The polarization work associated with the solute transfer from the gas-phase to the polar solvent, a term neglected in the framework of additive force fields, was found to make a large contribution to the total solvation free energy, comparable to the polar solute-solvent solvation contribution. The Drude PB/SASA also reproduces well the electronic polarization from the explicit solvent simulations of a small protein, BPTI. Model validation was based on comparisons with the experimental relative binding free energies of 371 single alanine mutations. With the Drude PB/SASA model the root mean square deviation between the predicted and experimental relative binding free energies is 3.35 kcal/mol, lower than 5.11 kcal/mol computed with the CHARMM36 additive force field. Overall, the results indicate that the main limitation of the Drude PB/SASA model is the inability of the SASA term to accurately capture non-polar solvation effects. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Alexey Aleksandrov
- Laboratoire d'Optique et Biosciences, CNRS, INSERM, Ecole Polytechnique, Palaiseau F-91128, France
| | - Fang-Yu Lin
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, Maryland 21201
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, Gordon Center for Integrative Science, 929 E57th Street, University of Chicago, Chicago, Illinois 60637
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, Maryland 21201
| |
Collapse
|