1
|
Del Campo CMZM, Nicolson GL, Sfera A. Neurolipidomics in schizophrenia: A not so well-oiled machine. Neuropharmacology 2024; 260:110117. [PMID: 39153730 DOI: 10.1016/j.neuropharm.2024.110117] [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: 04/17/2024] [Revised: 08/03/2024] [Accepted: 08/14/2024] [Indexed: 08/19/2024]
Abstract
Most patients with schizophrenia (SCZ) do not exhibit violent behaviors and are more likely to be victims rather than perpetrators of violent acts. However, a subgroup of forensic detainees with SCZ exhibit tendencies to engage in criminal violations. Although numerous models have been proposed, ranging from substance use, serotonin transporter gene, and cognitive dysfunction, the molecular underpinnings of violence in SCZ patients remains elusive. Lithium and clozapine have established anti-aggression properties and recent studies have linked low cholesterol levels and ultraviolet (UV) radiation with human aggression, while vitamin D3 reduces violent behaviors. A recent study found that vitamin D3, omega-3 fatty acids, magnesium, and zinc lower aggression in forensic population. In this review article, we take a closer look at aryl hydrocarbon receptor (AhR) and the dysfunctional lipidome in neuronal membranes, with emphasis on cholesterol and vitamin D3 depletion, as sources of aggressive behavior. We also discuss modalities to increase the fluidity of neuronal double layer via membrane lipid replacement (MLR) and natural or synthetic compounds. This article is part of the Special Issue on "Personality Disorders".
Collapse
Affiliation(s)
| | - Garth L Nicolson
- Department of Molecular Pathology, The Institute for Molecular Medicine, Huntington Beach, CA, 92647, USA
| | - Adonis Sfera
- Patton State Hospital, Loma Linda University, Department of Psychiatry, University of California, Riverside, USA.
| |
Collapse
|
2
|
Chan GKL. Spiers Memorial Lecture: Quantum chemistry, classical heuristics, and quantum advantage. Faraday Discuss 2024. [PMID: 39258407 DOI: 10.1039/d4fd00141a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/12/2024]
Abstract
We describe the problems of quantum chemistry, the intuition behind classical heuristic methods used to solve them, a conjectured form of the classical complexity of quantum chemistry problems, and the subsequent opportunities for quantum advantage. This article is written for both quantum chemists and quantum information theorists. In particular, we attempt to summarize the domain of quantum chemistry problems as well as the chemical intuition that is applied to solve them within concrete statements (such as a classical heuristic cost conjecture) in the hope that this may stimulate future analysis.
Collapse
Affiliation(s)
- Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California, USA
| |
Collapse
|
3
|
Ramdhani D, Watabe H, Hardianto A, Janitra RS. Complexation of 3p- C-NETA with radiometal ions: A density functional theory study for targeted radioimmunotherapy. Heliyon 2024; 10:e34875. [PMID: 39144950 PMCID: PMC11320446 DOI: 10.1016/j.heliyon.2024.e34875] [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: 12/18/2023] [Revised: 07/13/2024] [Accepted: 07/17/2024] [Indexed: 08/16/2024] Open
Abstract
Bifunctional chelators (BFCs) are vital in the design of effective radiopharmaceuticals, as they are able to bind to both a radiometal ion and a targeting vector. The 3p-C-NETA or 4-[2-(bis-carboxy-methylamino)-5-(4-nitrophenyl)-entyl])-7-carboxymethyl-[1,4,7]tri-azonan-1-yl acetic acid is a novel and promising BFC, developed for diagnostic and therapeutic purposes. The binding affinity between the BFC and radiometal ion significantly impacts their effectiveness. Predicting the equilibrium constants for the formation of 1:1 radiometals/chelator complexes (log K1 values) is crucial for designing BFCs with improved affinity and selectivity for radiometals. The purpose of this study is to evaluate the complexation of Ga3+, Tb3+, Bi3+, and Ac3+ radiometal ions with 3p-C-NETA using density functional theory (B3LYP and M06-HF functional) and 6-311G(d)/SDD basis sets, where the 1,4,7,10-tetrazacyclodecane-1,4,7,10-tetracetic acid (DOTA) was employed as a benchmark. Formation of the [Ac3+(3p-C-NETA)(H2O)]- complexes is predicted to be markedly less stable compared to the other complexes, exhibiting the lowest chemical hardness and the highest chemical softness. Additionally, the chelation stability of the complexes is mainly determined by ligand-ion and ion-water interactions, which depend on the atomic charge and atomic radius of the metal ion.
Collapse
Affiliation(s)
- Danni Ramdhani
- Department of Pharmaceutical Analysis and Medicinal Chemistry, Faculty of Pharmacy, Universitas Padjadjaran, Sumedang, Indonesia
- Division of Radiation Protection and Safety Control, Cyclotron and Radioisotope Center (CYRIC), Tohoku University, Sendai, Japan
| | - Hiroshi Watabe
- Division of Radiation Protection and Safety Control, Cyclotron and Radioisotope Center (CYRIC), Tohoku University, Sendai, Japan
| | - Ari Hardianto
- Department of Chemistry, Faculty of Mathematics and Natural Sciences, Universitas Padjadjaran, Sumedang, Indonesia
| | - Regaputra S. Janitra
- Department of Chemistry, Faculty of Mathematics and Natural Sciences, Universitas Padjadjaran, Sumedang, Indonesia
| |
Collapse
|
4
|
Prat A, Abdel Aty H, Bastas O, Kamuntavičius G, Paquet T, Norvaišas P, Gasparotto P, Tal R. HydraScreen: A Generalizable Structure-Based Deep Learning Approach to Drug Discovery. J Chem Inf Model 2024; 64:5817-5831. [PMID: 39037942 DOI: 10.1021/acs.jcim.4c00481] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/24/2024]
Abstract
We propose HydraScreen, a deep-learning framework for safe and robust accelerated drug discovery. HydraScreen utilizes a state-of-the-art 3D convolutional neural network designed for the effective representation of molecular structures and interactions in protein-ligand binding. We designed an end-to-end pipeline for high-throughput screening and lead optimization, targeting applications in structure-based drug design. We assessed our approach using established public benchmarks based on the CASF-2016 core set, achieving top-tier results in affinity and pose prediction (Pearson's r = 0.86, RMSE = 1.15, Top-1 = 0.95). We introduced a novel approach for interaction profiling, aimed at detecting potential biases within both the model and data sets. This approach not only enhanced interpretability but also reinforced the impartiality of our methodology. Finally, we demonstrated HydraScreen's ability to generalize effectively across novel proteins and ligands through a temporal split. We also provide insights into potential avenues for future development aimed at enhancing the robustness of machine learning scoring functions. HydraScreen (accessible at http://hydrascreen.ro5.ai/paper) provides a user-friendly GUI and a public API, facilitating the easy-access assessment of protein-ligand complexes.
Collapse
Affiliation(s)
- Alvaro Prat
- AI Chemistry, Ro5 2801 Gateway Drive, Irving, 75063 Texas, United States
| | - Hisham Abdel Aty
- AI Chemistry, Ro5 2801 Gateway Drive, Irving, 75063 Texas, United States
| | - Orestis Bastas
- AI Chemistry, Ro5 2801 Gateway Drive, Irving, 75063 Texas, United States
| | | | - Tanya Paquet
- AI Chemistry, Ro5 2801 Gateway Drive, Irving, 75063 Texas, United States
| | - Povilas Norvaišas
- AI Chemistry, Ro5 2801 Gateway Drive, Irving, 75063 Texas, United States
| | - Piero Gasparotto
- AI Chemistry, Ro5 2801 Gateway Drive, Irving, 75063 Texas, United States
| | - Roy Tal
- AI Chemistry, Ro5 2801 Gateway Drive, Irving, 75063 Texas, United States
| |
Collapse
|
5
|
Ginex T, Vázquez J, Estarellas C, Luque FJ. Quantum mechanical-based strategies in drug discovery: Finding the pace to new challenges in drug design. Curr Opin Struct Biol 2024; 87:102870. [PMID: 38914031 DOI: 10.1016/j.sbi.2024.102870] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2024] [Revised: 06/02/2024] [Accepted: 06/04/2024] [Indexed: 06/26/2024]
Abstract
The expansion of the chemical space to tangible libraries containing billions of synthesizable molecules opens exciting opportunities for drug discovery, but also challenges the power of computer-aided drug design to prioritize the best candidates. This directly hits quantum mechanics (QM) methods, which provide chemically accurate properties, but subject to small-sized systems. Preserving accuracy while optimizing the computational cost is at the heart of many efforts to develop high-quality, efficient QM-based strategies, reflected in refined algorithms and computational approaches. The design of QM-tailored physics-based force fields and the coupling of QM with machine learning, in conjunction with the computing performance of supercomputing resources, will enhance the ability to use these methods in drug discovery. The challenge is formidable, but we will undoubtedly see impressive advances that will define a new era.
Collapse
Affiliation(s)
- Tiziana Ginex
- Pharmacelera, Parc Científic de Barcelona (PCB), Baldiri Reixac 4-8, 08028 Barcelona, Spain
| | - Javier Vázquez
- Pharmacelera, Parc Científic de Barcelona (PCB), Baldiri Reixac 4-8, 08028 Barcelona, Spain; Departament de Nutrició, Ciències de l'Alimentació i Gastronomia, Universitat de Barcelona, Institut de Biomedicina (IBUB), 08921 Santa Coloma de Gramenet, Spain; Institut de Biomedicina (IBUB), 08921 Santa Coloma de Gramenet, Spain
| | - Carolina Estarellas
- Departament de Nutrició, Ciències de l'Alimentació i Gastronomia, Universitat de Barcelona, Institut de Biomedicina (IBUB), 08921 Santa Coloma de Gramenet, Spain; Institut de Química Teòrica i Computacional (IQTCUB), 08921 Santa Coloma de Gramenet, Spain
| | - F Javier Luque
- Departament de Nutrició, Ciències de l'Alimentació i Gastronomia, Universitat de Barcelona, Institut de Biomedicina (IBUB), 08921 Santa Coloma de Gramenet, Spain; Institut de Biomedicina (IBUB), 08921 Santa Coloma de Gramenet, Spain; Institut de Química Teòrica i Computacional (IQTCUB), 08921 Santa Coloma de Gramenet, Spain.
| |
Collapse
|
6
|
Thai QM, Nguyen TH, Phung HTT, Pham MQ, Pham NKT, Horng JT, Ngo ST. MedChemExpress compounds prevent neuraminidase N1 via physics- and knowledge-based methods. RSC Adv 2024; 14:18950-18956. [PMID: 38873542 PMCID: PMC11167619 DOI: 10.1039/d4ra02661f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2024] [Accepted: 06/07/2024] [Indexed: 06/15/2024] Open
Abstract
Influenza A viruses spread out worldwide, causing several global concerns. Hence, discovering neuraminidase inhibitors to prevent the influenza A virus is of great interest. In this work, a machine learning model was employed to evaluate the ligand-binding affinity of ca. 10 000 compounds from the MedChemExpress (MCE) database for inhibiting neuraminidase. Atomistic simulations, including molecular docking and molecular dynamics simulations, then confirmed the ligand-binding affinity. Furthermore, we clarified the physical insights into the binding process of ligands to neuraminidase. It was found that five compounds, including micronomicin, didesmethyl cariprazine, argatroban, Kgp-IN-1, and AY 9944, are able to inhibit neuraminidase N1 of the influenza A virus. Ten residues, including Glu119, Asp151, Arg152, Trp179, Gln228, Glu277, Glu278, Arg293, Asn295, and Tyr402, may be very important in controlling the ligand-binding process to N1.
Collapse
Affiliation(s)
- Quynh Mai Thai
- Laboratory of Biophysics, Institute for Advanced Study in Technology, Ton Duc Thang University Ho Chi Minh City Vietnam
- Faculty of Pharmacy, Ton Duc Thang University Ho Chi Minh City Vietnam
| | - Trung Hai Nguyen
- Laboratory of Biophysics, Institute for Advanced Study in Technology, Ton Duc Thang University Ho Chi Minh City Vietnam
- Faculty of Pharmacy, Ton Duc Thang University Ho Chi Minh City Vietnam
| | | | - Minh Quan Pham
- Institute of Natural Products Chemistry, Vietnam Academy of Science and Technology Hanoi Vietnam
- Graduate University of Science and Technology, Vietnam Academy of Science and Technology Hanoi Vietnam
| | - Nguyen Kim Tuyen Pham
- Faculty of Environment, Sai Gon University 273 An Duong Vuong, Ward 3, District 5 Ho Chi Minh City Vietnam
| | - Jim-Tong Horng
- Department of Biochemistry and Molecular Biology, College of Medicine, Chang Gung University Kweishan Taoyuan Taiwan
| | - Son Tung Ngo
- Laboratory of Biophysics, Institute for Advanced Study in Technology, Ton Duc Thang University Ho Chi Minh City Vietnam
- Faculty of Pharmacy, Ton Duc Thang University Ho Chi Minh City Vietnam
| |
Collapse
|
7
|
Wang M, Mei Y, Ryde U. Convergence criteria for single-step free-energy calculations: the relation between the Π bias measure and the sample variance. Chem Sci 2024; 15:8786-8799. [PMID: 38873060 PMCID: PMC11168088 DOI: 10.1039/d4sc00140k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2024] [Accepted: 05/08/2024] [Indexed: 06/15/2024] Open
Abstract
Free energy calculations play a crucial role in simulating chemical processes, enzymatic reactions, and drug design. However, assessing the reliability and convergence of these calculations remains a challenge. This study focuses on single-step free-energy calculations using thermodynamic perturbation. It explores how the sample distributions influence the estimated results and evaluates the reliability of various convergence criteria, including Kofke's bias measure Π and the standard deviation of the energy difference ΔU, σ ΔU . The findings reveal that for Gaussian distributions, there is a straightforward relationship between Π and σ ΔU , free energies can be accurately approximated using a second-order cumulant expansion, and reliable results are attainable for σ ΔU up to 25 kcal mol-1. However, interpreting non-Gaussian distributions is more complex. If the distribution is skewed towards more positive values than a Gaussian, converging the free energy becomes easier, rendering standard convergence criteria overly stringent. Conversely, distributions that are skewed towards more negative values than a Gaussian present greater challenges in achieving convergence, making standard criteria unreliable. We propose a practical approach to assess the convergence of estimated free energies.
Collapse
Affiliation(s)
- Meiting Wang
- School of Medical Engineering & Henan International Joint Laboratory of Neural Information Analysis and Drug Intelligent Design, Xinxiang Medical University Xinxiang 453003 China
- Department of Computational Chemistry, Lund University, Chemical Centre P.O. Box 124 SE-221 00 Lund Sweden
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University Shanghai 200241 China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai Shanghai 200062 China
- Collaborative Innovation Center of Extreme Optics, Shanxi University Taiyuan Shanxi 030006 China
| | - Ulf Ryde
- Department of Computational Chemistry, Lund University, Chemical Centre P.O. Box 124 SE-221 00 Lund Sweden
| |
Collapse
|
8
|
Kairys V, Baranauskiene L, Kazlauskiene M, Zubrienė A, Petrauskas V, Matulis D, Kazlauskas E. Recent advances in computational and experimental protein-ligand affinity determination techniques. Expert Opin Drug Discov 2024; 19:649-670. [PMID: 38715415 DOI: 10.1080/17460441.2024.2349169] [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: 03/18/2024] [Accepted: 04/25/2024] [Indexed: 05/22/2024]
Abstract
INTRODUCTION Modern drug discovery revolves around designing ligands that target the chosen biomolecule, typically proteins. For this, the evaluation of affinities of putative ligands is crucial. This has given rise to a multitude of dedicated computational and experimental methods that are constantly being developed and improved. AREAS COVERED In this review, the authors reassess both the industry mainstays and the newest trends among the methods for protein - small-molecule affinity determination. They discuss both computational affinity predictions and experimental techniques, describing their basic principles, main limitations, and advantages. Together, this serves as initial guide to the currently most popular and cutting-edge ligand-binding assays employed in rational drug design. EXPERT OPINION The affinity determination methods continue to develop toward miniaturization, high-throughput, and in-cell application. Moreover, the availability of data analysis tools has been constantly increasing. Nevertheless, cross-verification of data using at least two different techniques and careful result interpretation remain of utmost importance.
Collapse
Affiliation(s)
- Visvaldas Kairys
- Department of Bioinformatics, Institute of Biotechnology, Life Sciences Center, Vilnius University, Vilnius, Lithuania
| | - Lina Baranauskiene
- Department of Biothermodynamics and Drug Design, Institute of Biotechnology, Life Sciences Center, Vilnius University, Vilnius, Lithuania
| | | | - Asta Zubrienė
- Department of Biothermodynamics and Drug Design, Institute of Biotechnology, Life Sciences Center, Vilnius University, Vilnius, Lithuania
| | - Vytautas Petrauskas
- Department of Biothermodynamics and Drug Design, Institute of Biotechnology, Life Sciences Center, Vilnius University, Vilnius, Lithuania
| | - Daumantas Matulis
- Department of Biothermodynamics and Drug Design, Institute of Biotechnology, Life Sciences Center, Vilnius University, Vilnius, Lithuania
| | - Egidijus Kazlauskas
- Department of Biothermodynamics and Drug Design, Institute of Biotechnology, Life Sciences Center, Vilnius University, Vilnius, Lithuania
| |
Collapse
|
9
|
Chan B, Dawson W, Nakajima T. Sorting drug conformers in enzyme active sites: the XTB way. Phys Chem Chem Phys 2024; 26:12610-12618. [PMID: 38597505 DOI: 10.1039/d4cp00930d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/11/2024]
Abstract
In the present study, we have used the MEI196 set of interaction energies to investigate low-cost computational chemistry approaches for the calculation of binding between a molecule and its environment. Density functional theory (DFT) methods, when used with the vDZP basis set, yield good agreement with the reference energies. On the other hand, semi-empirical methods are less accurate as expected. By examining different groups of systems within MEI196 that contain species of a similar nature, we find that chemical similarity leads to cancellation of errors in the calculation of relative binding energies. Importantly, the semi-empirical method GFN1-xTB (XTB1) yields reasonable results for this purpose. We have thus further assessed the performance of XTB1 for calculating relative energies of docking poses of substrates in enzyme active sites represented by cluster models or within the ONIOM protocol. The results support the observations on error cancellation. This paves the way for the use of XTB1 in parts of large-scale virtual screening workflows to accelerate the drug discovery process.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki 852-8521, Japan.
- RIKEN Center for Computational Science, 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe, 650-0047, Japan
| | - William Dawson
- RIKEN Center for Computational Science, 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe, 650-0047, Japan
| | - Takahito Nakajima
- RIKEN Center for Computational Science, 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe, 650-0047, Japan
| |
Collapse
|
10
|
Pecina A, Fanfrlík J, Lepšík M, Řezáč J. SQM2.20: Semiempirical quantum-mechanical scoring function yields DFT-quality protein-ligand binding affinity predictions in minutes. Nat Commun 2024; 15:1127. [PMID: 38321025 PMCID: PMC10847445 DOI: 10.1038/s41467-024-45431-8] [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: 07/20/2023] [Accepted: 01/24/2024] [Indexed: 02/08/2024] Open
Abstract
Accurate estimation of protein-ligand binding affinity is the cornerstone of computer-aided drug design. We present a universal physics-based scoring function, named SQM2.20, addressing key terms of binding free energy using semiempirical quantum-mechanical computational methods. SQM2.20 incorporates the latest methodological advances while remaining computationally efficient even for systems with thousands of atoms. To validate it rigorously, we have compiled and made available the PL-REX benchmark dataset consisting of high-resolution crystal structures and reliable experimental affinities for ten diverse protein targets. Comparative assessments demonstrate that SQM2.20 outperforms other scoring methods and reaches a level of accuracy similar to much more expensive DFT calculations. In the PL-REX dataset, it achieves excellent correlation with experimental data (average R2 = 0.69) and exhibits consistent performance across all targets. In contrast to DFT, SQM2.20 provides affinity predictions in minutes, making it suitable for practical applications in hit identification or lead optimization.
Collapse
Affiliation(s)
- Adam Pecina
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Prague, Czech Republic
| | - Jindřich Fanfrlík
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Prague, Czech Republic
| | - Martin Lepšík
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Prague, Czech Republic
| | - Jan Řezáč
- Institute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Prague, Czech Republic.
| |
Collapse
|
11
|
Amaya-Rodriguez CA, Carvajal-Zamorano K, Bustos D, Alegría-Arcos M, Castillo K. A journey from molecule to physiology and in silico tools for drug discovery targeting the transient receptor potential vanilloid type 1 (TRPV1) channel. Front Pharmacol 2024; 14:1251061. [PMID: 38328578 PMCID: PMC10847257 DOI: 10.3389/fphar.2023.1251061] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Accepted: 12/14/2023] [Indexed: 02/09/2024] Open
Abstract
The heat and capsaicin receptor TRPV1 channel is widely expressed in nerve terminals of dorsal root ganglia (DRGs) and trigeminal ganglia innervating the body and face, respectively, as well as in other tissues and organs including central nervous system. The TRPV1 channel is a versatile receptor that detects harmful heat, pain, and various internal and external ligands. Hence, it operates as a polymodal sensory channel. Many pathological conditions including neuroinflammation, cancer, psychiatric disorders, and pathological pain, are linked to the abnormal functioning of the TRPV1 in peripheral tissues. Intense biomedical research is underway to discover compounds that can modulate the channel and provide pain relief. The molecular mechanisms underlying temperature sensing remain largely unknown, although they are closely linked to pain transduction. Prolonged exposure to capsaicin generates analgesia, hence numerous capsaicin analogs have been developed to discover efficient analgesics for pain relief. The emergence of in silico tools offered significant techniques for molecular modeling and machine learning algorithms to indentify druggable sites in the channel and for repositioning of current drugs aimed at TRPV1. Here we recapitulate the physiological and pathophysiological functions of the TRPV1 channel, including structural models obtained through cryo-EM, pharmacological compounds tested on TRPV1, and the in silico tools for drug discovery and repositioning.
Collapse
Affiliation(s)
- Cesar A. Amaya-Rodriguez
- Centro Interdisciplinario de Neurociencia de Valparaíso, Facultad de Ciencias, Universidad de Valparaíso, Valparaíso, Chile
- Departamento de Fisiología y Comportamiento Animal, Facultad de Ciencias Naturales, Exactas y Tecnología, Universidad de Panamá, Ciudad de Panamá, Panamá
| | - Karina Carvajal-Zamorano
- Centro Interdisciplinario de Neurociencia de Valparaíso, Facultad de Ciencias, Universidad de Valparaíso, Valparaíso, Chile
| | - Daniel Bustos
- Centro de Investigación de Estudios Avanzados del Maule (CIEAM), Vicerrectoría de Investigación y Postgrado Universidad Católica del Maule, Talca, Chile
- Laboratorio de Bioinformática y Química Computacional, Departamento de Medicina Traslacional, Facultad de Medicina, Universidad Católica del Maule, Talca, Chile
| | - Melissa Alegría-Arcos
- Núcleo de Investigación en Data Science, Facultad de Ingeniería y Negocios, Universidad de las Américas, Santiago, Chile
| | - Karen Castillo
- Centro Interdisciplinario de Neurociencia de Valparaíso, Facultad de Ciencias, Universidad de Valparaíso, Valparaíso, Chile
- Centro de Investigación de Estudios Avanzados del Maule (CIEAM), Vicerrectoría de Investigación y Postgrado Universidad Católica del Maule, Talca, Chile
| |
Collapse
|
12
|
da Rocha JM, Campos DMDO, Esmaile SC, Menezes GDL, Bezerra KS, da Silva RA, Junior EDDS, Tayyeb JZ, Akash S, Fulco UL, Alqahtani T, Oliveira JIN. Quantum biochemical analysis of the binding interactions between a potential inhibitory drug and the Ebola viral glycoprotein. J Biomol Struct Dyn 2024:1-17. [PMID: 38258414 DOI: 10.1080/07391102.2024.2305314] [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: 09/19/2023] [Accepted: 01/08/2024] [Indexed: 01/24/2024]
Abstract
Ebola virus disease (EVD) causes outbreaks and epidemics in West Africa that persist until today. The envelope glycoprotein of Ebola virus (GP) consists of two subunits, GP1 and GP2, and plays a key role in anchoring or fusing the virus to the host cell in its active form on the virion surface. Toremifene (TOR) is a ligand that mainly acts as an estrogen receptor antagonist; however, a recent study showed a strong and efficient interaction with GP. In this context, we aimed to evaluate the energetic affinity features involved in the interaction between GP and toremifene by computer simulation techniques using the Molecular Fractionation Method with Conjugate Caps (MFCC) scheme and quantum-mechanical (QM) calculations, as well as missense mutations to assess protein stability. We identified ASP522, GLU100, TYR517, THR519, LEU186, LEU515 as the most attractive residues in the EBOV glycoprotein structure that form the binding pocket. We divided toremifene into three regions and evaluated that region i was more important than region iii and region ii for the formation of the TOR-GP1/GP2 complex, which might control the molecular remodeling process of TOR. The mutations that caused more destabilization were ARG134, LEU515, TYR517 and ARG559, while those that caused stabilization were GLU523 and ASP522. TYR517 is a critical residue for the binding of TOR, and is highly conserved among EBOV species. Our results may help to elucidate the mechanism of drug action on the GP protein of the Ebola virus and subsequently develop new pharmacological approaches against EVD.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Jaerdyson M da Rocha
- Department of Biophysics and Pharmacology, Bioscience Center, Federal University of Rio Grande do Norte, Natal, RN, Brazil
| | - Daniel M de O Campos
- Department of Biophysics and Pharmacology, Bioscience Center, Federal University of Rio Grande do Norte, Natal, RN, Brazil
| | - Stephany C Esmaile
- Department of Biophysics and Pharmacology, Bioscience Center, Federal University of Rio Grande do Norte, Natal, RN, Brazil
| | - Gabriela de L Menezes
- Department of Biophysics and Pharmacology, Bioscience Center, Federal University of Rio Grande do Norte, Natal, RN, Brazil
| | - Katyanna S Bezerra
- Department of Biophysics and Pharmacology, Bioscience Center, Federal University of Rio Grande do Norte, Natal, RN, Brazil
| | - Roosevelt A da Silva
- Core Collaboratives of BioSistemas, Special Unit of Exact Sciences, Federal University of Jataí, Jataí, GO, Brazil
| | - Edilson D da S Junior
- Department of Biophysics and Pharmacology, Bioscience Center, Federal University of Rio Grande do Norte, Natal, RN, Brazil
| | - Jehad Zuhair Tayyeb
- Department of Clinical Biochemistry, College of Medicine, University of Jeddah, Jeddah, Saudi Arabia
| | - Shopnil Akash
- Department of Pharmacy, Faculty of Allied Health Sciences, Daffodil International University, Birulia, Ashulia, Dhaka, Bangladesh
| | - Umberto L Fulco
- Department of Biophysics and Pharmacology, Bioscience Center, Federal University of Rio Grande do Norte, Natal, RN, Brazil
| | - Taha Alqahtani
- Department of Pharmacology, College of Pharmacy, King Khalid University, Abha, Saudi Arabia
| | - Jonas I N Oliveira
- Department of Biophysics and Pharmacology, Bioscience Center, Federal University of Rio Grande do Norte, Natal, RN, Brazil
| |
Collapse
|
13
|
Lima Neto JX, Bezerra KS, Barbosa ED, Araujo RL, Galvão DS, Lyra ML, Oliveira JIN, Akash S, Jardan YAB, Nafidi HA, Bourhia M, Fulco UL. Investigation of protein-protein interactions and hotspot region on the NSP7-NSP8 binding site in NSP12 of SARS-CoV-2. Front Mol Biosci 2024; 10:1325588. [PMID: 38304231 PMCID: PMC10830813 DOI: 10.3389/fmolb.2023.1325588] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2023] [Accepted: 12/22/2023] [Indexed: 02/03/2024] Open
Abstract
Background: The RNA-dependent RNA polymerase (RdRp) complex, essential in viral transcription and replication, is a key target for antiviral therapeutics. The core unit of RdRp comprises the nonstructural protein NSP12, with NSP7 and two copies of NSP8 (NSP81 and NSP82) binding to NSP12 to enhance its affinity for viral RNA and polymerase activity. Notably, the interfaces between these subunits are highly conserved, simplifying the design of molecules that can disrupt their interaction. Methods: We conducted a detailed quantum biochemical analysis to characterize the interactions within the NSP12-NSP7, NSP12-NSP81, and NSP12-NSP82 dimers. Our objective was to ascertain the contribution of individual amino acids to these protein-protein interactions, pinpointing hotspot regions crucial for complex stability. Results: The analysis revealed that the NSP12-NSP81 complex possessed the highest total interaction energy (TIE), with 14 pairs of residues demonstrating significant energetic contributions. In contrast, the NSP12-NSP7 complex exhibited substantial interactions in 8 residue pairs, while the NSP12-NSP82 complex had only one pair showing notable interaction. The study highlighted the importance of hydrogen bonds and π-alkyl interactions in maintaining these complexes. Intriguingly, introducing the RNA sequence with Remdesivir into the complex resulted in negligible alterations in both interaction energy and geometric configuration. Conclusion: Our comprehensive analysis of the RdRp complex at the protein-protein interface provides invaluable insights into interaction dynamics and energetics. These findings can guide the design of small molecules or peptide/peptidomimetic ligands to disrupt these critical interactions, offering a strategic pathway for developing effective antiviral drugs.
Collapse
Affiliation(s)
- José Xavier Lima Neto
- Department of Biophysics and Pharmacology, Bioscience Center, Federal University of Rio Grande do Norte, Natal, Brazil
| | - Katyanna Sales Bezerra
- Department of Biophysics and Pharmacology, Bioscience Center, Federal University of Rio Grande do Norte, Natal, Brazil
| | - Emmanuel Duarte Barbosa
- Department of Biophysics and Pharmacology, Bioscience Center, Federal University of Rio Grande do Norte, Natal, Brazil
| | - Roniel Lima Araujo
- Department of Biophysics and Pharmacology, Bioscience Center, Federal University of Rio Grande do Norte, Natal, Brazil
| | | | | | - Jonas Ivan Nobre Oliveira
- Department of Biophysics and Pharmacology, Bioscience Center, Federal University of Rio Grande do Norte, Natal, Brazil
| | - Shopnil Akash
- Department of Pharmacy, Daffodil International University, Dhaka, Bangladesh
| | - Yousef A. Bin Jardan
- Department of Pharmaceutics, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
| | - Hiba-Allah Nafidi
- Department of Food Science, Faculty of Agricultural and Food Sciences, Laval University, Quebec City, QC, Canada
| | - Mohammed Bourhia
- Department of Chemistry and Biochemistry, Faculty of Medicine and Pharmacy, Ibn Zohr University, Laayoune, Morocco
| | - Umberto Laino Fulco
- Department of Biophysics and Pharmacology, Bioscience Center, Federal University of Rio Grande do Norte, Natal, Brazil
| |
Collapse
|
14
|
Mohebbinia Z, Firouzi R, Karimi-Jafari MH. Improving protein-ligand docking results using the Semiempirical quantum mechanics: testing on the PDBbind 2016 core set. J Biomol Struct Dyn 2024:1-11. [PMID: 38165642 DOI: 10.1080/07391102.2023.2299742] [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: 10/30/2023] [Accepted: 12/20/2023] [Indexed: 01/04/2024]
Abstract
Molecular docking techniques are routinely employed for predicting ligand binding conformations and affinities in the in silico phase of the drug design and development process. In this study, a reliable semiempirical quantum mechanics (SQM) method, PM7, was employed for geometry optimization of top-ranked poses obtained from two widely used docking programs, AutoDock4 and AutoDock Vina. The PDBbind core set (version 2016), which contains high-quality crystal protein - ligand complexes with their corresponding experimental binding affinities, was used as an initial dataset in this research. It was shown that docking pose optimization improves the accuracy of pose predictions and is very useful for the refinement of docked complexes via removing clashes between ligands and proteins. It was also demonstrated that AutoDock Vina achieves a higher sampling power than AutoDock4 in generating accurate ligand poses (RMSD ≤ 2.0 Å), while AutoDock4 exhibits a better ranking power than AutoDock Vina. Finally, a new protocol based on a combination of the results obtained from the two docking programs was proposed for structure-based virtual screening studies, which benefits from the robust sampling abilities of AutoDock Vina and the reliable ranking performance of AutoDock4.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Zainab Mohebbinia
- Department of Physical Chemistry, Chemistry and Chemical Engineering Research Center of Iran, Tehran, Iran
| | - Rohoullah Firouzi
- Department of Physical Chemistry, Chemistry and Chemical Engineering Research Center of Iran, Tehran, Iran
| | | |
Collapse
|
15
|
Shakour N, Hoseinpoor S, Rajabian F, Azimi SG, Iranshahi M, Sadeghi-Aliabadi H, Hadizadeh F. Discovery of non-peptide GLP-1r natural agonists for enhancing coronary safety in type 2 diabetes patients. J Biomol Struct Dyn 2024:1-18. [PMID: 38165453 DOI: 10.1080/07391102.2023.2298734] [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: 09/21/2023] [Accepted: 12/17/2023] [Indexed: 01/03/2024]
Abstract
This study explores the computational discovery of non-peptide agonists targeting the Glucagon-Like Peptide-1 Receptor (GLP-1R) to enhance the safety of major coronary outcomes in individuals affected by Type 2 Diabetes. The objective is to identify novel compounds that can activate the GLP-1R pathway without the limitations associated with peptide agonists. Type 2 diabetes mellitus (T2DM) is associated with an increased risk of cardiovascular disease (CVD) and mortality, which is attributed to the accumulation of fat in organs, including the heart. Glucagon-like peptide-1 receptor agonists (GLP-1RAs) are frequently used to manage T2DM and could potentially offer cardiovascular benefits. Therefore, this study examines non-peptide agonists of GLP-1R to improve coronary safety in type 2 diabetes patients. After rigorous assessments, two standout candidates were identified, with natural compound 12 emerging as the most promising. This study represents a notable advancement in enhancing the management of coronary outcomes among individuals with type 2 diabetes. The computational methodology employed successfully pinpointed potential GLP-1R natural agonists, providing optimism for the development of safer and more effective therapeutic interventions. Although computational methodologies have provided crucial insights, realizing the full potential of these compounds requires extensive experimental investigations, crucial in advancing therapeutic strategies for this critical patient population.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Neda Shakour
- Department of Medicinal Chemistry, School of Pharmacy, Mashhad University of Medical Sciences, Mashhad, Iran
- Student Research Committee, Faculty of Medicine, Mashhad University of Medical Sciences, Mashhad, Iran
| | - Saeideh Hoseinpoor
- Department of Biochemistry and Biophysics, Faculty of Sciences, Mashhad Branch, Islamic Azad University, Mashhad, Iran
- International UNESCO Center for Health-Related Basic Sciences and Human Nutrition, Mashhad University of Medical Sciences, Mashhad, Iran
| | - Fatemeh Rajabian
- Department of Pharmacodynamics and Toxicology, School of Pharmacy, Mashhad University of Medical Sciences, Mashhad, Iran
| | - Sabikeh G Azimi
- Department of Chemistry, Faculty of Sciences, University of Birjand, Birjand, Iran
| | - Mehrdad Iranshahi
- Biotechnology Research Center, Pharmaceutical Technology Institute, Mashhad University of Medical Sciences, Mashhad, Iran
| | - Hojjat Sadeghi-Aliabadi
- Department of Pharmaceutical Chemistry, Isfahan University of Medical Sciences, Isfahan, Iran
| | - Farzin Hadizadeh
- Department of Medicinal Chemistry, School of Pharmacy, Mashhad University of Medical Sciences, Mashhad, Iran
- Biotechnology Research Center, Pharmaceutical Technology Institute, Mashhad University of Medical Sciences, Mashhad, Iran
| |
Collapse
|
16
|
Rubina, Moin ST. Attempting Well-Tempered Funnel Metadynamics Simulations for the Evaluation of the Binding Kinetics of Methionine Aminopeptidase-II Inhibitors. J Chem Inf Model 2023; 63:7729-7743. [PMID: 38059911 DOI: 10.1021/acs.jcim.3c01130] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/08/2023]
Abstract
Understanding the unbinding kinetics of protein-ligand complexes is considered a significant approach for the design of ligands with desired specificity and safety. In recent years, enhanced sampling methods have emerged as effective tools for studying the unbinding kinetics of protein-ligand complexes at the atomistic level. MetAP-II is a target for the treatment of cancer for which not a single effective drug is available yet. The identification of the dissociation rate of ligands from the complexes often serves as a better predictor for in vivo efficacy than the ligands' binding affinity. Here, funnel-based restraint well-tempered metadynamics simulations were applied to predict the residence time of two ligands bound to MetAP-II, along with the ligand association and dissociation mechanism involving the identification of the binding hotspot during ligand egress. The ligand-egressing route revealed by metadynamics simulations also correlated with the identified pathways from the CAVER analysis and by the enhanced sampling simulation using PLUMED. Ligand 1 formed a strong H-bond interaction with GLU364 estimating a higher residence time of 28.22 ± 5.29 ns in contrast to ligand 2 with a residence time of 19.05 ± 3.58 ns, which easily dissociated from the binding pocket of MetAP-II. The results obtained from the simulations were consistent to reveal ligand 1 being superior to ligand 2; however, the experimental data related to residence time were close for both ligands, and no kinetic data were available for ligand 2. The current study could be considered the first attempt to apply an enhanced sampling method for the evaluation of the binding kinetics and thermodynamics of two different classes of ligands to a binuclear metalloprotein.
Collapse
Affiliation(s)
- Rubina
- Third World Center for Science and Technology H.E.J. Research Institute of Chemistry, International Center for Chemical and Biological Science University of Karachi, Karachi 75270, Pakistan
| | - Syed Tarique Moin
- Third World Center for Science and Technology H.E.J. Research Institute of Chemistry, International Center for Chemical and Biological Science University of Karachi, Karachi 75270, Pakistan
| |
Collapse
|
17
|
Yuan Y, Cui Q. Accurate and Efficient Multilevel Free Energy Simulations with Neural Network-Assisted Enhanced Sampling. J Chem Theory Comput 2023; 19:5394-5406. [PMID: 37527495 PMCID: PMC10810721 DOI: 10.1021/acs.jctc.3c00591] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/03/2023]
Abstract
Free energy differences (ΔF) are essential to quantitative characterization and understanding of chemical and biological processes. Their direct estimation with an accurate quantum mechanical potential is of great interest and yet impractical due to high computational cost and incompatibility with typical alchemical free energy protocols. One promising solution is the multilevel free energy simulation in which the estimate of ΔF at an inexpensive low level of theory is combined with the correction toward a higher level of theory. The poor configurational overlap generally expected between the two levels of theory, however, presents a major challenge. We overcome this challenge by using a deep neural network model and enhanced sampling simulations. An adversarial autoencoder is used to identify a low-dimensional (latent) space that compactly represents the degrees of freedom that encode the distinct distributions at the two levels of theory. Enhanced sampling in this latent space is then used to drive the sampling of configurations that predominantly contribute to the free energy correction. Results for both gas phase and condensed phase systems demonstrate that this data-driven approach offers high accuracy and efficiency with great potential for scalability to complex systems.
Collapse
Affiliation(s)
- Yuchen Yuan
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| | - Qiang Cui
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
- Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
- Department of Biomedical Engineering, Boston University, 44 Cummington Mall, Boston, Massachusetts 02215, United States
| |
Collapse
|
18
|
Liu J, Wan J, Ren Y, Shao X, Xu X, Rao L. DOX_BDW: Incorporating Solvation and Desolvation Effects of Cavity Water into Nonfitting Protein-Ligand Binding Affinity Prediction. J Chem Inf Model 2023; 63:4850-4863. [PMID: 37539963 DOI: 10.1021/acs.jcim.3c00776] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/05/2023]
Abstract
Accurate prediction of the protein-ligand binding affinity (PLBA) with an affordable cost is one of the ultimate goals in the field of structure-based drug design (SBDD), as well as a great challenge in the computational and theoretical chemistry. Herein, we have systematically addressed the complicated solvation and desolvation effects on the PLBA brought by the difference of the explicit water in the protein cavity before and after ligands bind to the protein-binding site. Based on the new solvation model, a nonfitting method at the first-principles level for the PLBA prediction was developed by taking the bridging and displaced water (BDW) molecules into account simultaneously. The newly developed method, DOX_BDW, was validated against a total of 765 noncovalent and covalent protein-ligand binding pairs, including the CASF2016 core set, Cov_2022 covalent binding testing set, and six testing sets for the hit and lead compound optimization (HLO) simulation. In all of the testing sets, the DOX_BDW method was able to produce PLBA predictions that were strongly correlated with the corresponding experimental data (R = 0.66-0.85). The overall performance of DOX_BDW is better than the current empirical scoring functions that are heavily parameterized. DOX_BDW is particularly outstanding for the covalent binding situation, implying the need for considering an electronic structure in covalent drug design. Furthermore, the method is especially recommended to be used in the HLO scenario of SBDD, where hundreds of similar derivatives need to be screened and refined. The computational cost of DOX_BDW is affordable, and its accuracy is remarkable.
Collapse
Affiliation(s)
- Jiaqi Liu
- National Key Laboratory of Green Pesticide, Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, Hubei International Scientific and Technological Cooperation Base of Pesticide and Green Synthesis, College of Chemistry, Central China Normal University, Wuhan 43009, People's Republic of China
| | - Jian Wan
- National Key Laboratory of Green Pesticide, Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, Hubei International Scientific and Technological Cooperation Base of Pesticide and Green Synthesis, College of Chemistry, Central China Normal University, Wuhan 43009, People's Republic of China
| | - Yanliang Ren
- National Key Laboratory of Green Pesticide, Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, Hubei International Scientific and Technological Cooperation Base of Pesticide and Green Synthesis, College of Chemistry, Central China Normal University, Wuhan 43009, People's Republic of China
| | - Xubo Shao
- National Key Laboratory of Green Pesticide, Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, Hubei International Scientific and Technological Cooperation Base of Pesticide and Green Synthesis, College of Chemistry, Central China Normal University, Wuhan 43009, People's Republic of China
| | - Xin Xu
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Ministry of Education (MOE) Laboratory for Computational Physical Science, Department of Chemistry, Fudan University, Shanghai 200433, People's Republic of China
| | - Li Rao
- National Key Laboratory of Green Pesticide, Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, Hubei International Scientific and Technological Cooperation Base of Pesticide and Green Synthesis, College of Chemistry, Central China Normal University, Wuhan 43009, People's Republic of China
| |
Collapse
|
19
|
Vasile S, Roos K. Understanding the Structure-Activity Relationship through Density Functional Theory: A Simple Method Predicts Relative Binding Free Energies of Metalloenzyme Fragment-like Inhibitors. ACS OMEGA 2023; 8:21438-21449. [PMID: 37360476 PMCID: PMC10285960 DOI: 10.1021/acsomega.2c08156] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/23/2022] [Accepted: 03/29/2023] [Indexed: 06/28/2023]
Abstract
Despite being involved in several human diseases, metalloenzymes are targeted by a small percentage of FDA-approved drugs. Development of novel and efficient inhibitors is required, as the chemical space of metal binding groups (MBGs) is currently limited to four main classes. The use of computational chemistry methods in drug discovery has gained momentum thanks to accurate estimates of binding modes and binding free energies of ligands to receptors. However, exact predictions of binding free energies in metalloenzymes are challenging due to the occurrence of nonclassical phenomena and interactions that common force field-based methods are unable to correctly describe. In this regard, we applied density functional theory (DFT) to predict the binding free energies and to understand the structure-activity relationship of metalloenzyme fragment-like inhibitors. We tested this method on a set of small-molecule inhibitors with different electronic properties and coordinating two Mn2+ ions in the binding site of the influenza RNA polymerase PAN endonuclease. We modeled the binding site using only atoms from the first coordination shell, hence reducing the computational cost. Thanks to the explicit treatment of electrons by DFT, we highlighted the main contributions to the binding free energies and the electronic features differentiating strong and weak inhibitors, achieving good qualitative correlation with the experimentally determined affinities. By introducing automated docking, we explored alternative ways to coordinate the metal centers and we identified 70% of the highest affinity inhibitors. This methodology provides a fast and predictive tool for the identification of key features of metalloenzyme MBGs, which can be useful for the design of new and efficient drugs targeting these ubiquitous proteins.
Collapse
|
20
|
Bao X, Liu X, Dou R, Xu S, Liu D, Luo J, Gong X, Wong CF, Zhou B. How are N-methylcarbamates encapsulated by β-cyclodextrin: insight into the binding mechanism. Phys Chem Chem Phys 2023; 25:13923-13932. [PMID: 37184134 DOI: 10.1039/d3cp01252b] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
Guest molecules containing chromophore groups encapsulated by β-cyclodextrin (β-CD) generate circular dichroism (CD) signals, which enables a preliminary prediction of their binding modes. However, the accurate determination of the representative binding conformation (RC) remains a challenging task due to the complex conformational space of these host-guest systems. Here, we combine a molecular dynamics/quantum mechanics/continuum solvent model (MD/QM/CSM) with induced circular dichroism (ICD) data (N. L. Pacioni, A. B. Pierini and A. V. Veglia, Spectrochim. Acta A Mol. Biomol. Spectrosc., 2013, 103, 319-324.) to explore the binding mechanism of β-CD with four N-methylcarbamate molecules: promecarb (PC), bendiocarb (BC), carbaryl (CY) and carbofuran (CF). In aqueous solution, their stability decreases as: PC > BC > CY > CF. Comparing the ECD spectra computed from TD-DFT with the ICD data can help eliminate many common binding configurations and identify the RC. The host-guest binding affinities (BAs) estimated using a ONIOM2(B971:PM6)/SMD model reproduce the measured binding trend, reveal the competition between the non-covalent interaction and solvent effect and explain the large difference in their binding modes. We also examine the fluctuations in the computed BA using similar structures.
Collapse
Affiliation(s)
- Xiaofang Bao
- Computational Institute for Molecules and Materials, School of Chemistry and Chemical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China.
| | - Xiao Liu
- Computational Institute for Molecules and Materials, School of Chemistry and Chemical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China.
| | - Ran Dou
- Computational Institute for Molecules and Materials, School of Chemistry and Chemical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China.
| | - Sen Xu
- Computational Institute for Molecules and Materials, School of Chemistry and Chemical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China.
| | - Dabin Liu
- Computational Institute for Molecules and Materials, School of Chemistry and Chemical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China.
| | - Jun Luo
- Computational Institute for Molecules and Materials, School of Chemistry and Chemical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China.
| | - Xuedong Gong
- Computational Institute for Molecules and Materials, School of Chemistry and Chemical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China.
| | - Chung F Wong
- Department of Chemistry and Biochemistry and Center for Nanoscience, University of Missouri-Saint Louis, One University Boulevard, Saint Louis, MO 63121, USA
| | - Baojing Zhou
- Computational Institute for Molecules and Materials, School of Chemistry and Chemical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China.
| |
Collapse
|
21
|
Gurushankar K, Jeyaseelan SC, Grishina M, Siswanto I, Tiwari R, Puspaningsih NNT. Density Functional Theory, Molecular Dynamics and AlteQ Studies Approaches of Baimantuoluoamide A and Baimantuoluoamide B to Identify Potential Inhibitors of M pro Proteins: a Novel Target for the Treatment of SARS COVID-19. JETP LETTERS 2023; 117:1-10. [PMID: 37360903 PMCID: PMC10184967 DOI: 10.1134/s0021364023600039] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/02/2023] [Revised: 04/12/2023] [Accepted: 04/13/2023] [Indexed: 06/28/2023]
Abstract
COVID-19 has resulted in epidemi conditions over the world. Despite efforts by scientists from all over the world to develop an effective va ine against this virus, there is presently no recognized cure for COVID-19. The most succeed treatments for various ailments come from natural components found in medicinal plants, which are also rucial for the development of new medications. This study intends to understand the role of the baimantuoluoamide A and baimantuoluoamide B molecules in the treatment of Covid19. Initially, density functional theory (DFT) used to explore their electronic potentials along with the Becke3-Lee-Yang-Parr (B3LYP) 6-311 + G(d, p) basis set. A number of characteristics, including the energy gap, hardness, local softness, electronegativity, and electrophilicity, have also been calculated to discuss the reactivity of mole ules. Using natural bond orbital, the title compound's bioactive nature and stability were investigated. Further, both compounds potential inhibitors with main protease (Mpro) proteins, molecular dynamics simulations and AlteQ investigations also studied. Supplementary Information The online version contains supplementary material available at 10.1134/S0021364023600039.
Collapse
Affiliation(s)
- K. Gurushankar
- Laboratory of Computational Modeling of Drugs, Higher Medical and Biological School, South Ural State University, 454080 Chelyabinsk, Russia
- Department of Physics, Kalasalingam Academy of Research and Education, 626126 Krishnankoil, Tamilnadu India
| | - S. Ch. Jeyaseelan
- Post Graduate & Research Department of Physics, N.M.S.S.V.N. College, 625019 Madurai, Tamilnadu India
- Post Graduate Department of Physics, Mannar Thirumalai Naciker College, 625004 Madurai, Tamilnadu India
| | - M. Grishina
- Laboratory of Computational Modeling of Drugs, Higher Medical and Biological School, South Ural State University, 454080 Chelyabinsk, Russia
| | - I. Siswanto
- Bioinformati Laboratory, UCoE Research Center for Bio-Molecule Engineering Universitas Airlangga, 60115 Surabaya, Indonesia
| | - R. Tiwari
- Department of Physics, Coordinator Research and Development Cell, Dr CV Raman University, 495113 Kargi Kota, Bilaspur CG India
| | - N. N. T. Puspaningsih
- Department of Chemistry, Faculty of Science and Technology, Universitas Airlangga, 60115 Surabaya, Indonesia
| |
Collapse
|
22
|
Vornweg JR, Wolter M, Jacob CR. A simple and consistent quantum-chemical fragmentation scheme for proteins that includes two-body contributions. J Comput Chem 2023; 44:1634-1644. [PMID: 37171574 DOI: 10.1002/jcc.27114] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Revised: 03/28/2023] [Accepted: 03/30/2023] [Indexed: 05/13/2023]
Abstract
The Molecular Fractionation with Conjugate Caps (MFCC) method is a popular fragmentation method for the quantum-chemical treatment of proteins. However, it does not account for interactions between the amino acid fragments, such as intramolecular hydrogen bonding. Here, we present a combination of the MFCC fragmentation scheme with a second-order many-body expansion (MBE) that consistently accounts for all fragment-fragment, fragment-cap, and cap-cap interactions, while retaining the overall simplicity of the MFCC scheme with its chemically meaningful fragments. We show that with the resulting MFCC-MBE(2) scheme, the errors in the total energies of selected polypeptides and proteins can be reduced by up to one order of magnitude and relative energies of different protein conformers can be predicted accurately.
Collapse
Affiliation(s)
- Johannes R Vornweg
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Braunschweig, Germany
| | - Mario Wolter
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Braunschweig, Germany
| | - Christoph R Jacob
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Braunschweig, Germany
| |
Collapse
|
23
|
Ariai J, Gellrich U. The entropic penalty for associative reactions and their physical treatment during routine computations. Phys Chem Chem Phys 2023; 25:14005-14015. [PMID: 37161492 DOI: 10.1039/d3cp00970j] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
A systematic study of the entropic penalty for associative reactions is presented. It is shown that computed solution-phase Gibbs free energies typically overestimate entropic contributions. This entropic penalty for associative reactions in solution, i.e., if the number of particles decreases along the reaction coordinate (sum of stoichiometric numbers ), originates from the insufficient treatment of entropic effects by implicit solvent models. We propose an additive correction scheme to Gibbs free energies that is suitable for routine applications by non-expert users. This correction is based on Garza's formalism for the solution-phase entropy [A. J. Garza, J. Chem. Theory Comput., 2019, 15, 3204.] that is physically sound and embedded into an efficient black-box type algorithm. To critically evaluate the entropic penalty and its proposed treatment, we compiled an experimental benchmark set of 31 ΔrG and 22 in 15 different solvents. Using a representative best-practice computational protocol (at wave function theory (WFT) based DLPNO-CCSD(T) and density functional theory (DFT) based revDSD-PBEP86-D4 level with an implicit solvent model), we determined a sizeable entropic penalty ranging from 2-11 kcal mol-1. Using the correction scheme presented herein, the entropic penalty is corrected to the chemical accuracy of ≤1 kcal mol-1 (WFT and DFT). The same applies to at the WFT level. Barriers at the DFT level are overestimated by 2 kcal mol-1 (classic) and underestimated by 2 kcal mol-1 (corrected). This effect is attributed to the finding that barriers computed at the DFT level are systematically 2-3 kcal mol-1 lower than barriers obtained with WFT.
Collapse
Affiliation(s)
- Jama Ariai
- Institute of Organic Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany.
| | - Urs Gellrich
- Institute of Organic Chemistry, Justus Liebig University, Heinrich-Buff-Ring 17, 35392 Giessen, Germany.
| |
Collapse
|
24
|
Hok L, Vianello R. Selective Deuteration Improves the Affinity of Adenosine A 2A Receptor Ligands: A Computational Case Study with Istradefylline and Caffeine. J Chem Inf Model 2023; 63:3138-3149. [PMID: 37155356 DOI: 10.1021/acs.jcim.3c00424] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
We used a range of computational techniques to assess the effect of selective C-H deuteration on the antagonist istradefylline affinity for the adenosine A2A receptor, which was discussed relative to its structural analogue caffeine, a well-known and likely the most widely used stimulant. The obtained results revealed that smaller caffeine shows high receptor flexibility and exchanges between two distinct poses, which agrees with crystallographic data. In contrast, the additional C8-trans-styryl fragment in istradefylline locks the ligand within a uniform binding pose, while contributing to the affinity through the C-H···π and π···π contacts with surface residues, which, together with its much lower hydration prior to binding, enhances the affinity over caffeine. In addition, the aromatic C8-unit shows a higher deuteration sensitivity over the xanthine part, so when both of its methoxy groups are d6-deuterated, the affinity improvement is -0.4 kcal mol-1, which surpasses the overall affinity gain of -0.3 kcal mol-1 in the perdeuterated d9-caffeine. Yet, the latter predicts around 1.7-fold potency increase, being relevant for its pharmaceutical implementations, and also those within the coffee and energy drink production industries. Still, the full potential of our strategy is achieved in polydeuterated d19-istradefylline, whose A2A affinity improves by -0.6 kcal mol-1, signifying a 2.8-fold potency increase that strongly promotes it as a potential synthetic target. This knowledge supports deuterium application in drug design, and while the literature already reports about over 20 deuterated drugs currently in the clinical development, it is easily foreseen that more examples will hit the market in the years to come. With this in mind, we propose that the devised computational methodology, involving the ONIOM division of the QM region for the ligand and the MM region for its environment, with an implicit quantization of nuclear motions relevant for the H/D exchange, allows fast and efficient estimates of the binding isotope effects in any biological system.
Collapse
Affiliation(s)
- Lucija Hok
- Laboratory for the Computational Design and Synthesis of Functional Materials, Division of Organic Chemistry and Biochemistry, Ruđer Bošković Institute, 10000 Zagreb, Croatia
| | - Robert Vianello
- Laboratory for the Computational Design and Synthesis of Functional Materials, Division of Organic Chemistry and Biochemistry, Ruđer Bošković Institute, 10000 Zagreb, Croatia
| |
Collapse
|
25
|
Meelua W, Wanjai T, Thinkumrob N, Oláh J, Cairns JRK, Hannongbua S, Ryde U, Jitonnom J. A computational study of the reaction mechanism and stereospecificity of dihydropyrimidinase. Phys Chem Chem Phys 2023; 25:8767-8778. [PMID: 36912034 DOI: 10.1039/d2cp05262h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/06/2023]
Abstract
Dihydropyrimidinase (DHPase) is a key enzyme in the pyrimidine pathway, the catabolic route for synthesis of β-amino acids. It catalyses the reversible conversion of 5,6-dihydrouracil (DHU) or 5,6-dihydrothymine (DHT) to the corresponding N-carbamoyl-β-amino acids. This enzyme has the potential to be used as a tool in the production of β-amino acids. Here, the reaction mechanism and origin of stereospecificity of DHPases from Saccharomyces kluyveri and Sinorhizobium meliloti CECT4114 were investigated and compared using a quantum mechanical cluster approach based on density functional theory. Two models of the enzyme active site were designed from the X-ray crystal structure of the native enzyme: a small cluster to characterize the mechanism and the stationary points and a large model to probe the stereospecificity and the role of stereo-gate-loop (SGL) residues. It is shown that a hydroxide ion first performs a nucleophilic attack on the substrate, followed by the abstraction of a proton by Asp358, which occurs concertedly with protonation of the ring nitrogen by the same residue. For the DHT substrate, the enzyme displays a preference for the L-configuration, in good agreement with experimental observation. Comparison of the reaction energetics of the two models reveals the importance of SGL residues in the stereospecificity of catalysis. The role of the conserved Tyr172 residue in transition-state stabilization is confirmed as the Tyr172Phe mutation increases the activation barrier of the reaction by ∼8 kcal mol-1. A detailed understanding of the catalytic mechanism of the enzyme could offer insight for engineering in order to enhance its activity and substrate scope.
Collapse
Affiliation(s)
- Wijitra Meelua
- Demonstration School, University of Phayao, Phayao 56000, Thailand
- Unit of Excellence in Computational Molecular Science and Catalysis, and Division of Chemistry, School of Science, University of Phayao, Phayao 56000, Thailand.
| | - Tanchanok Wanjai
- Unit of Excellence in Computational Molecular Science and Catalysis, and Division of Chemistry, School of Science, University of Phayao, Phayao 56000, Thailand.
| | - Natechanok Thinkumrob
- Unit of Excellence in Computational Molecular Science and Catalysis, and Division of Chemistry, School of Science, University of Phayao, Phayao 56000, Thailand.
| | - Julianna Oláh
- Department of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics, Műegyetem rakpart 3, Budapest H-1111, Hungary
| | - James R Ketudat Cairns
- Center for Biomolecular Structure, Function and Application and School of Chemistry, Institute of Science, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand
| | - Supa Hannongbua
- Department of Chemistry, Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, Lund SE-221 00, Sweden
| | - Jitrayut Jitonnom
- Unit of Excellence in Computational Molecular Science and Catalysis, and Division of Chemistry, School of Science, University of Phayao, Phayao 56000, Thailand.
| |
Collapse
|
26
|
Nakamura S, Akaki T, Nishiwaki K, Nakatani M, Kawase Y, Takahashi Y, Nakanishi I. System truncation accelerates binding affinity calculations with the fragment molecular orbital method: A benchmark study. J Comput Chem 2023; 44:824-831. [PMID: 36444861 DOI: 10.1002/jcc.27044] [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: 08/07/2022] [Revised: 11/01/2022] [Accepted: 11/03/2022] [Indexed: 11/30/2022]
Abstract
The fragment molecular orbital (FMO) method is a fast quantum-mechanical method that divides systems into pieces of fragments and performs ab initio calculations. The system truncation enables further speed improvement. In this article, we systematically study the effects of system truncations on binding affinity calculations obtained with FMO in combination with either the polarizable continuum model (FMO/PCM) or in combination with the Møller-Plesset method (FMO-MP2). We have used five protein complexes with ligands of several charged states. The calculated binding energies of the size variants of the truncated system, including only a restricted number of atoms around the ligand, are compared to the energy obtained from a full system. The result shows that the systems could be truncated to a radius of 8 Å from neutral ligands within an error of 0.7 kcal/mol, and 12 Å from charged ligands within an error of 1.1 kcal/mol for calculating the binding energy in solution.
Collapse
Affiliation(s)
- Shinya Nakamura
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Tatsuo Akaki
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan.,Chemical Research Laboratories, Central Pharmaceutical Research Institute, Japan Tobacco Inc., Osaka, Japan
| | - Keiji Nishiwaki
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Midori Nakatani
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Yuji Kawase
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Yuki Takahashi
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Isao Nakanishi
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| |
Collapse
|
27
|
Losev TV, Gerasimov IS, Panova MV, Lisov AA, Abdyusheva YR, Rusina PV, Zaletskaya E, Stroganov OV, Medvedev MG, Novikov FN. Quantum Mechanical-Cluster Approach to Solve the Bioisosteric Replacement Problem in Drug Design. J Chem Inf Model 2023; 63:1239-1248. [PMID: 36763797 DOI: 10.1021/acs.jcim.2c01212] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/12/2023]
Abstract
Bioisosteres are molecules that differ in substituents but still have very similar shapes. Bioisosteric replacements are ubiquitous in modern drug design, where they are used to alter metabolism, change bioavailability, or modify activity of the lead compound. Prediction of relative affinities of bioisosteres with computational methods is a long-standing task; however, the very shape closeness makes bioisosteric substitutions almost intractable for computational methods, which use standard force fields. Here, we design a quantum mechanical (QM)-cluster approach based on the GFN2-xTB semi-empirical quantum-chemical method and apply it to a set of H → F bioisosteric replacements. The proposed methodology enables advanced prediction of biological activity change upon bioisosteric substitution of -H with -F, with the standard deviation of 0.60 kcal/mol, surpassing the ChemPLP scoring function (0.83 kcal/mol), and making QM-based ΔΔG estimation comparable to ∼0.42 kcal/mol standard deviation of in vitro experiment. The speed of the method and lack of tunable parameters makes it affordable in current drug research.
Collapse
Affiliation(s)
- Timofey V Losev
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation.,Department of Chemistry, Lomonosov Moscow State University, Leninskie Gory 1/3, 119991 Moscow, Russian Federation.,A.N. Nesmeyanov Institute of Organoelement Compounds of Russian Academy of Sciences, Vavilov Str. 28, 119991 Moscow, Russian Federation
| | - Igor S Gerasimov
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation.,Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Maria V Panova
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation
| | - Alexey A Lisov
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation
| | - Yana R Abdyusheva
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation.,National Research University Higher School of Economics, Myasnitskaya Street 20, 101000 Moscow, Russian Federation
| | - Polina V Rusina
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation
| | - Eugenia Zaletskaya
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation.,National Research University Higher School of Economics, Myasnitskaya Street 20, 101000 Moscow, Russian Federation
| | - Oleg V Stroganov
- BioMolTech Corp., 226 York Mills Rd, Toronto, Ontario M2L 1L1, Canada
| | - Michael G Medvedev
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation
| | - Fedor N Novikov
- N.D. Zelinsky Institute of Organic Chemistry of Russian Academy of Sciences, Leninsky prospect 47, 119991 Moscow, Russian Federation.,National Research University Higher School of Economics, Myasnitskaya Street 20, 101000 Moscow, Russian Federation
| |
Collapse
|
28
|
Wang D, Li W, Dong X, Li H, Hu L. TFRegNCI: Interpretable Noncovalent Interaction Correction Multimodal Based on Transformer Encoder Fusion. J Chem Inf Model 2023; 63:782-793. [PMID: 36652718 DOI: 10.1021/acs.jcim.2c01283] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
The interpretability is an important issue for end-to-end learning models. Motivated by computer vision algorithms, an interpretable noncovalent interaction (NCI) correction multimodal (TFRegNCI) is proposed for NCI prediction. TFRegNCI is based on RegNet feature extraction and a transformer encoder fusion strategy. RegNet is a network design paradigm that mainly focuses on local features. Meanwhile, the Vision Transformer is also leveraged for feature extraction, because it can capture global features better than RegNet while lowering the computational cost. Using a transformer encoder as the fusion strategy rather than multilayer perceptron can enhance model performance, due to its emphasis on important features with less parameters. Therefore, the proposed TFRegNCI achieved high accurate prediction (mean absolute error of ∼0.1 kcal/mol) comparing with the coupled cluster single double (triple) (CCSD(T)) benchmark. To further improve the model efficiency, TFRegNCI applies two-dimensional (2D) inputs transformed from three-dimensional (3D) electron density cubes, which saves time (30%), while the model accuracy remains. To improve model interpretability, a visualization module, Gradient-weighted Regression Activation Mapping (Grad-RAM) has been embedded. Grad-RAM is promoted from the classification algorithm, Gradient-weighted Class Activation Mapping, to perform feature visualization for the regression task. With Grad-RAM, the visual location map for features in deep learning models can be displayed. The feature map visualizations suggest that the 2D model has the similar performance as the 3D model, because of equally effective feature extractions from electron density. Moreover, the valid feature region on the location map by the 3D model is consistent with the NCIPLOT NCI isosurface. It is confirmed that the model does extract significant features related to the NCI interaction. The interpretable analyses are carried out through molecular orbital contribution on effective features. Thereby, the proposed model is likely to be a promising tool to reveal some essential information on NCIs, with regard to the level of electronic theory.
Collapse
Affiliation(s)
- Donghan Wang
- School of Information Science and Technology, Northeast Normal University, Changchun130117, China
| | - Wenze Li
- College of Computer and Information Engineering, Henan Normal University, Henan, Xinxiang453007, China
| | - Xu Dong
- School of Information Science and Technology, Northeast Normal University, Changchun130117, China
| | - Hongzhi Li
- School of Information Science and Technology, Northeast Normal University, Changchun130117, China
| | - LiHong Hu
- School of Information Science and Technology, Northeast Normal University, Changchun130117, China
| |
Collapse
|
29
|
Han Y, Wang Z, Chen A, Ali I, Cai J, Ye S, Wei Z, Li J. A deep transfer learning-based protocol accelerates full quantum mechanics calculation of protein. Brief Bioinform 2023; 24:6901901. [PMID: 36516300 DOI: 10.1093/bib/bbac532] [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: 08/12/2022] [Revised: 10/07/2022] [Accepted: 11/07/2022] [Indexed: 12/15/2022] Open
Abstract
Effective full quantum mechanics (FQM) calculation of protein remains a grand challenge and of great interest in computational biology with substantial applications in drug discovery, protein dynamic simulation and protein folding. However, the huge computational complexity of the existing QM methods impends their applications in large systems. Here, we design a transfer-learning-based deep learning (TDL) protocol for effective FQM calculations (TDL-FQM) on proteins. By incorporating a transfer-learning algorithm into deep neural network (DNN), the TDL-FQM protocol is capable of performing calculations at any given accuracy using models trained from small datasets with high-precision and knowledge learned from large amount of low-level calculations. The high-level double-hybrid DFT functional and high-level quality of basis set is used in this work as a case study to evaluate the performance of TDL-FQM, where the selected 15 proteins are predicted to have a mean absolute error of 0.01 kcal/mol/atom for potential energy and an average root mean square error of 1.47 kcal/mol/$ {\rm A^{^{ \!\!\!o}}} $ for atomic forces. The proposed TDL-FQM approach accelerates the FQM calculation more than thirty thousand times faster in average and presents more significant benefits in efficiency as the size of protein increases. The ability to learn knowledge from one task to solve related problems demonstrates that the proposed TDL-FQM overcomes the limitation of standard DNN and has a strong power to predict proteins with high precision, which solves the challenge of high precision prediction in large chemical and biological systems.
Collapse
Affiliation(s)
- Yanqiang Han
- National Key Laboratory of Science and Technology on Micro/Nano Fabrication, Shanghai Jiao Tong University, Shanghai 200240, China.,Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Zhilong Wang
- National Key Laboratory of Science and Technology on Micro/Nano Fabrication, Shanghai Jiao Tong University, Shanghai 200240, China.,Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China
| | - An Chen
- National Key Laboratory of Science and Technology on Micro/Nano Fabrication, Shanghai Jiao Tong University, Shanghai 200240, China.,Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Imran Ali
- National Key Laboratory of Science and Technology on Micro/Nano Fabrication, Shanghai Jiao Tong University, Shanghai 200240, China.,Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Junfei Cai
- National Key Laboratory of Science and Technology on Micro/Nano Fabrication, Shanghai Jiao Tong University, Shanghai 200240, China.,Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Simin Ye
- National Key Laboratory of Science and Technology on Micro/Nano Fabrication, Shanghai Jiao Tong University, Shanghai 200240, China.,Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China
| | - Zhiyun Wei
- Shanghai Key Laboratory of Maternal Fetal Medicine, Shanghai First Maternity and Infant Hospital, School of Medicine, Tongji University, Shanghai 200092, China
| | - Jinjin Li
- National Key Laboratory of Science and Technology on Micro/Nano Fabrication, Shanghai Jiao Tong University, Shanghai 200240, China.,Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China
| |
Collapse
|
30
|
Nguyen TH, Tam NM, Tuan MV, Zhan P, Vu VV, Quang DT, Ngo ST. Searching for potential inhibitors of SARS-COV-2 main protease using supervised learning and perturbation calculations. Chem Phys 2023; 564:111709. [PMID: 36188488 PMCID: PMC9511900 DOI: 10.1016/j.chemphys.2022.111709] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2022] [Revised: 08/11/2022] [Accepted: 09/21/2022] [Indexed: 11/28/2022]
Abstract
Inhibiting the biological activity of SARS-CoV-2 Mpro can prevent viral replication. In this context, a hybrid approach using knowledge- and physics-based methods was proposed to characterize potential inhibitors for SARS-CoV-2 Mpro. Initially, supervised machine learning (ML) models were trained to predict a ligand-binding affinity of ca. 2 million compounds with the correlation on a test set of R = 0.748 ± 0.044 . Atomistic simulations were then used to refine the outcome of the ML model. Using LIE/FEP calculations, nine compounds from the top 100 ML inhibitors were suggested to bind well to the protease with the domination of van der Waals interactions. Furthermore, the binding affinity of these compounds is also higher than that of nirmatrelvir, which was recently approved by the US FDA to treat COVID-19. In addition, the ligands altered the catalytic triad Cys145 - His41 - Asp187, possibly disturbing the biological activity of SARS-CoV-2.
Collapse
Affiliation(s)
- Trung Hai Nguyen
- Laboratory of Theoretical and Computational Biophysics, Advanced Institute of Materials Science, Ton Duc Thang University, Ho Chi Minh City, Viet Nam
- Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Viet Nam
| | - Nguyen Minh Tam
- Laboratory of Theoretical and Computational Biophysics, Advanced Institute of Materials Science, Ton Duc Thang University, Ho Chi Minh City, Viet Nam
- Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Viet Nam
| | - Mai Van Tuan
- Department of Microbiology, Hue Central Hospital, Hue City, Viet Nam
| | - Peng Zhan
- Department of Medicinal Chemistry, Key Laboratory of Chemical Biology (Ministry of Education), School of Pharmaceutical Sciences, Cheeloo College of Medicine, Shandong University, 44 West Culture Road, 250012 Jinan, Shandong, PR China
| | - Van V Vu
- NTT Hi-Tech Institute, Nguyen Tat Thanh University, Ho Chi Minh City, Viet Nam
| | - Duong Tuan Quang
- Department of Chemistry, Hue University, Thua Thien Hue Province, Hue City, Viet Nam
| | - Son Tung Ngo
- Laboratory of Theoretical and Computational Biophysics, Advanced Institute of Materials Science, Ton Duc Thang University, Ho Chi Minh City, Viet Nam
- Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Viet Nam
| |
Collapse
|
31
|
Molecular dynamics-based insight of VEGFR-2 kinase domain: a combined study of pharmacophore modeling and molecular docking and dynamics. J Mol Model 2022; 29:17. [PMID: 36550239 DOI: 10.1007/s00894-022-05427-x] [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: 06/19/2022] [Accepted: 12/17/2022] [Indexed: 12/24/2022]
Abstract
BACKGROUND Inhibition of vascular endothelial growth factor receptor 2 (VEGFR-2) tyrosine kinase by small molecules has become a promising target in the treatment of cancer. OBJECTIVE In this study, we approached pharmacophore modeling coupled with a structure-based virtual screening workflow to identify the potent inhibitors. METHODS The top selected hit compounds have been rescored using the MM/GBSA approach. To understand the molecular reactivity, electronic properties, and stability of those inhibitors, we have employed density functional theory and molecular dynamics. Following that, the best 21 hit compounds have been further post-processed with a Quantum ligand partial charge-based rescoring process and further validated by implementing molecular dynamics simulation. RESULTS The ten hit compounds have been hypothesized and considered as potent inhibitors of VEGFR-2 tyrosine kinase. This study also signifies the contribution of QM-based ligand partial charge, which is more accurate in predicting reliable free binding energy and filtering large ligand libraries to hit optimization, rather than assigning those of the force field-based method. From the binding pattern analysis of all the complexes, amino acids, such as Glu885, Cys919, Cys1045, Thr916, Thr919, and Asp1046, were found to have comprehensive interaction with the hit compounds. CONCLUSION Hence, this could prove to be useful as a potential inhibition site of the VEGFR-2 tyrosine kinase domain for future researchers. Moreover, this study also emphasizes the conformational changes upon ATP binding, based on either the receptor's rigidity or flexibility.
Collapse
|
32
|
Ngo ST, Nguyen TH, Tung NT, Vu VV, Pham MQ, Mai BK. Characterizing the ligand-binding affinity toward SARS-CoV-2 Mpro via physics- and knowledge-based approaches. Phys Chem Chem Phys 2022; 24:29266-29278. [PMID: 36449268 DOI: 10.1039/d2cp04476e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Computational approaches, including physics- and knowledge-based methods, have commonly been used to determine the ligand-binding affinity toward SARS-CoV-2 main protease (Mpro or 3CLpro). Strong binding ligands can thus be suggested as potential inhibitors for blocking the biological activity of the protease. In this context, this paper aims to provide a short review of computational approaches that have recently been applied in the search for inhibitor candidates of Mpro. In particular, molecular docking and molecular dynamics (MD) simulations are usually combined to predict the binding affinity of thousands of compounds. Quantitative structure-activity relationship (QSAR) is the least computationally demanding and therefore can be used for large chemical collections of ligands. However, its accuracy may not be high. Moreover, the quantum mechanics/molecular mechanics (QM/MM) method is most commonly used for covalently binding inhibitors, which also play an important role in inhibiting the activity of SARS-CoV-2. Furthermore, machine learning (ML) models can significantly increase the searching space of ligands with high accuracy for binding affinity prediction. Physical insights into the binding process can then be confirmed via physics-based calculations. Integration of ML models into computational chemistry provides many more benefits and can lead to new therapies sooner.
Collapse
Affiliation(s)
- Son Tung Ngo
- Laboratory of Theoretical and Computational Biophysics, Advanced Institute of Materials Science, Ton Duc Thang University, Ho Chi Minh City, Vietnam. .,Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Trung Hai Nguyen
- Laboratory of Theoretical and Computational Biophysics, Advanced Institute of Materials Science, Ton Duc Thang University, Ho Chi Minh City, Vietnam. .,Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Nguyen Thanh Tung
- Institute of Materials Science, Vietnam Academy of Science and Technology, Hanoi, Vietnam. .,Graduate University of Science and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Van V Vu
- NTT Hi-Tech Institute, Nguyen Tat Thanh University, Ho Chi Minh City, Vietnam
| | - Minh Quan Pham
- Graduate University of Science and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam.,Institute of Natural Products Chemistry, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Binh Khanh Mai
- Department of Chemistry, University of Pittsburgh, Pittsburgh, PA, USA
| |
Collapse
|
33
|
Dutkiewicz Z. Computational methods for calculation of protein-ligand binding affinities in structure-based drug design. PHYSICAL SCIENCES REVIEWS 2022. [DOI: 10.1515/psr-2020-0034] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Abstract
Drug design is an expensive and time-consuming process. Any method that allows reducing the time the costs of the drug development project can have great practical value for the pharmaceutical industry. In structure-based drug design, affinity prediction methods are of great importance. The majority of methods used to predict binding free energy in protein-ligand complexes use molecular mechanics methods. However, many limitations of these methods in describing interactions exist. An attempt to go beyond these limits is the application of quantum-mechanical description for all or only part of the analyzed system. However, the extensive use of quantum mechanical (QM) approaches in drug discovery is still a demanding challenge. This chapter briefly reviews selected methods used to calculate protein-ligand binding affinity applied in virtual screening (VS), rescoring of docked poses, and lead optimization stage, including QM methods based on molecular simulations.
Collapse
Affiliation(s)
- Zbigniew Dutkiewicz
- Department of Chemical Technology of Drugs , Poznan University of Medical Sciences , ul. Grunwaldzka 6 , 60-780 Poznań , Poznan , 60-780, Poland
| |
Collapse
|
34
|
Al-Ansi AY, Lin Z. MDO: A Computational Protocol for Prediction of Flexible Enzyme-Ligand Binding Mode. Curr Comput Aided Drug Des 2022; 18:CAD-EPUB-125919. [PMID: 36043706 DOI: 10.2174/1573409918666220827151546] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2022] [Revised: 06/15/2022] [Accepted: 07/05/2022] [Indexed: 11/22/2022]
Abstract
AIM Developing a method for use in computer aided drug design Background: Predicting the structure of enzyme-ligand binding mode is essential for understanding the properties, functions, and mechanisms of the bio-complex, but is rather difficult due to the enormous sampling space involved. OBJECTIVE Accurate prediction of enzyme-ligand binding mode conformation. METHOD A new computational protocol, MDO, is proposed for finding the structure of ligand binding pose. MDO consists of sampling enzyme sidechain conformations via molecular dynamics simulation of enzyme-ligand system and clustering of the enzyme configurations, sampling ligand binding poses via molecular docking and clustering of the ligand conformations, and the optimal ligand binding pose prediction via geometry optimization and ranking by the ONIOM method. MDO is tested on 15 enzyme-ligand complexes with known accurate structures. RESULTS The success rate of MDO predictions, with RMSD < 2 Å, is 67%, substantially higher than the 40% success rate of conventional methods. The MDO success rate can be increased to 83% if the ONIOM calculations are applied only for the starting poses with ligands inside the binding cavities. CONCLUSION The MDO protocol provides high quality enzyme-ligand binding mode prediction with reasonable computational cost. The MDO protocol is recommended for use in the structure-based drug design.
Collapse
Affiliation(s)
- Amar Y Al-Ansi
- Hefei National Laboratory for Physical Sciences at Microscale & CAS Key Laboratory of Strongly-Coupled Quantum Matter Physics, Department of Physics, University of Science and Technology of China, Hefei 230026, China
- Department of Physics, Sana'a University, Sana'a, Yemen
| | - Zijing Lin
- Hefei National Laboratory for Physical Sciences at Microscale & CAS Key Laboratory of Strongly-Coupled Quantum Matter Physics, Department of Physics, University of Science and Technology of China, Hefei 230026, China
| |
Collapse
|
35
|
Multistep orthophosphate release tunes actomyosin energy transduction. Nat Commun 2022; 13:4575. [PMID: 35931685 PMCID: PMC9356070 DOI: 10.1038/s41467-022-32110-9] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2021] [Accepted: 07/13/2022] [Indexed: 11/29/2022] Open
Abstract
Muscle contraction and a range of critical cellular functions rely on force-producing interactions between myosin motors and actin filaments, powered by turnover of adenosine triphosphate (ATP). The relationship between release of the ATP hydrolysis product ortophosphate (Pi) from the myosin active site and the force-generating structural change, the power-stroke, remains enigmatic despite its central role in energy transduction. Here, we present a model with multistep Pi-release that unifies current conflicting views while also revealing additional complexities of potential functional importance. The model is based on our evidence from kinetics, molecular modelling and single molecule fluorescence studies of Pi binding outside the active site. It is also consistent with high-speed atomic force microscopy movies of single myosin II molecules without Pi at the active site, showing consecutive snapshots of pre- and post-power stroke conformations. In addition to revealing critical features of energy transduction by actomyosin, the results suggest enzymatic mechanisms of potentially general relevance. Release of the ATP hydrolysis product orthophosphate (Pi) from the myosin active site is central in force generation but is poorly understood. Here, Moretto et al. present evidence for multistep Pi-release reconciling apparently contradictory results.
Collapse
|
36
|
Nguyen TH, Tran PT, Pham NQA, Hoang VH, Hiep DM, Ngo ST. Identifying Possible AChE Inhibitors from Drug-like Molecules via Machine Learning and Experimental Studies. ACS OMEGA 2022; 7:20673-20682. [PMID: 35755364 PMCID: PMC9219098 DOI: 10.1021/acsomega.2c00908] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/14/2022] [Accepted: 05/27/2022] [Indexed: 05/30/2023]
Abstract
Acetylcholinesterase (AChE) is one of the most important drug targets for Alzheimer's disease (AD) treatment. In this work, a machine learning model was trained to rapidly and accurately screen large chemical databases for the potential inhibitors of AChE. The obtained results were then validated via in vitro enzyme assay. Moreover, atomistic simulations including molecular docking and molecular dynamics simulations were then used to understand molecular insights into the binding process of ligands to AChE. In particular, two compounds including benzyl trifluoromethyl ketone and trifluoromethylstyryl ketone were indicated as highly potent inhibitors of AChE because they established IC50 values of 0.51 and 0.33 μM, respectively. The obtained IC50 of two compounds is significantly lower than that of galantamine (2.10 μM). The predicted log(BB) suggests that the compounds may be able to traverse the blood-brain barrier. A good agreement between computational and experimental studies was observed, indicating that the hybrid approach can enhance AD therapy.
Collapse
Affiliation(s)
- Trung Hai Nguyen
- Laboratory
of Theoretical and Computational Biophysics, Advanced Institute of
Materials Science, Ton Duc Thang
University, Ho Chi Minh City, Vietnam
- Faculty
of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Phuong-Thao Tran
- Hanoi
University of Pharmacy, 13-15 Le Thanh Tong, Hanoi 008404, Vietnam
| | - Ngoc Quynh Anh Pham
- Faculty
of Chemical Engineering, Ho Chi Minh City
University of Technology (HCMUT), Ho Chi Minh City 700000, Vietnam
| | - Van-Hai Hoang
- Faculty
of Pharmacy, Phenikka University, Hanoi 008404, Vietnam
- Phenikka
Institute for Advanced Study, Phenikka University, Hanoi 008404, Vietnam
| | - Dinh Minh Hiep
- Department
of Agriculture and Rural Development, Ho Chi Minh City 700000, Vietnam
| | - Son Tung Ngo
- Laboratory
of Theoretical and Computational Biophysics, Advanced Institute of
Materials Science, Ton Duc Thang
University, Ho Chi Minh City, Vietnam
- Faculty
of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| |
Collapse
|
37
|
Combining classical molecular docking with self-consistent charge density-functional tight-binding computations for the efficient and quality prediction of ligand binding structure. JOURNAL OF CHEMICAL RESEARCH 2022. [DOI: 10.1177/17475198221101999] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
To improve the successful prediction rate of the existing molecular docking methods, a new docking approach is proposed that consists of three steps: generating an ensemble of docked poses with a conventional docking method, performing clustering analysis of the ensemble to select the representative poses, and optimizing the representative structures with a low-cost quantum mechanics method. Three quantum mechanics methods, self-consistent charge density-functional tight-binding, ONIOM(DFT:PM6), and ONIOM(SCC-DFTB:PM6), are tested on 18 ligand-receptor bio-complexes. The rate of successful binding pose predictions by the proposed self-consistent charge density-functional tight-binding docking method is the highest, at 67%. The self-consistent charge density-functional tight-binding docking method should be useful for the structure-based drug design.
Collapse
|
38
|
Zhou Y, Jiang Y, Chen SJ. RNA-ligand molecular docking: advances and challenges. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2022; 12:e1571. [PMID: 37293430 PMCID: PMC10250017 DOI: 10.1002/wcms.1571] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/26/2021] [Accepted: 07/20/2021] [Indexed: 12/16/2022]
Abstract
With rapid advances in computer algorithms and hardware, fast and accurate virtual screening has led to a drastic acceleration in selecting potent small molecules as drug candidates. Computational modeling of RNA-small molecule interactions has become an indispensable tool for RNA-targeted drug discovery. The current models for RNA-ligand binding have mainly focused on the docking-and-scoring method. Accurate docking and scoring should tackle four crucial problems: (1) conformational flexibility of ligand, (2) conformational flexibility of RNA, (3) efficient sampling of binding sites and binding poses, and (4) accurate scoring of different binding modes. Moreover, compared with the problem of protein-ligand docking, predicting ligand binding to RNA, a negatively charged polymer, is further complicated by additional effects such as metal ion effects. Thermodynamic models based on physics-based and knowledge-based scoring functions have shown highly encouraging success in predicting ligand binding poses and binding affinities. Recently, kinetic models for ligand binding have further suggested that including dissociation kinetics (residence time) in ligand docking would result in improved performance in estimating in vivo drug efficacy. More recently, the rise of deep-learning approaches has led to new tools for predicting RNA-small molecule binding. In this review, we present an overview of the recently developed computational methods for RNA-ligand docking and their advantages and disadvantages.
Collapse
Affiliation(s)
- Yuanzhe Zhou
- Department of Physics and Astronomy, Department of Biochemistry, Institute of Data Sciences and Informatics, University of Missouri, Columbia, MO 65211-7010, USA
| | - Yangwei Jiang
- Department of Physics and Astronomy, Department of Biochemistry, Institute of Data Sciences and Informatics, University of Missouri, Columbia, MO 65211-7010, USA
| | - Shi-Jie Chen
- Department of Physics and Astronomy, Department of Biochemistry, Institute of Data Sciences and Informatics, University of Missouri, Columbia, MO 65211-7010, USA
| |
Collapse
|
39
|
Thai QM, Pham TNH, Hiep DM, Pham MQ, Tran PT, Nguyen TH, Ngo ST. Searching for AChE inhibitors from natural compounds by using machine learning and atomistic simulations. J Mol Graph Model 2022; 115:108230. [DOI: 10.1016/j.jmgm.2022.108230] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2022] [Revised: 05/19/2022] [Accepted: 05/20/2022] [Indexed: 12/14/2022]
|
40
|
Inverse Mixed-Solvent Molecular Dynamics for Visualization of the Residue Interaction Profile of Molecular Probes. Int J Mol Sci 2022; 23:ijms23094749. [PMID: 35563139 PMCID: PMC9103889 DOI: 10.3390/ijms23094749] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2022] [Revised: 04/18/2022] [Accepted: 04/23/2022] [Indexed: 02/01/2023] Open
Abstract
To ensure efficiency in discovery and development, the application of computational technology is essential. Although virtual screening techniques are widely applied in the early stages of drug discovery research, the computational methods used in lead optimization to improve activity and reduce the toxicity of compounds are still evolving. In this study, we propose a method to construct the residue interaction profile of the chemical structure used in the lead optimization by performing “inverse” mixed-solvent molecular dynamics (MSMD) simulation. Contrary to constructing a protein-based, atom interaction profile, we constructed a probe-based, protein residue interaction profile using MSMD trajectories. It provides us the profile of the preferred protein environments of probes without co-crystallized structures. We assessed the method using three probes: benzamidine, catechol, and benzene. As a result, the residue interaction profile of each probe obtained by MSMD was a reasonable physicochemical description of the general non-covalent interaction. Moreover, comparison with the X-ray structure containing each probe as a ligand shows that the map of the interaction profile matches the arrangement of amino acid residues in the X-ray structure.
Collapse
|
41
|
Chan B, Dawson W, Nakajima T. Searching for a Reliable Density Functional for Molecule-Environment Interactions, Found B97M-V/def2-mTZVP. J Phys Chem A 2022; 126:2397-2406. [PMID: 35390254 DOI: 10.1021/acs.jpca.2c02032] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In the present study, we have examined density functional theory methods for the calculation of the interaction energy between a small molecule and its environment. For simple systems such as a neutral solute in a neutral solvent, good accuracy can be attained using low-cost "3c" methods, in particular r2SCAN-3c. When part(s) of the system is charged, the accurate computation of the interactions is more challenging. In these cases, we find the B97M-V/def2-mTZVP method to agree well with reference values; it also shows good accuracy for the more straightforward neutral systems. Thus, B97M-V/def2-mTZVP provides a means for accurate and low-cost computation of interaction energies, notably the binding between a substrate or a drug molecule and an enzyme, which may facilitate rational drug design.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki 852-8521, Japan
| | - William Dawson
- RIKEN Center for Computational Science, 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe 650-0047, Japan
| | - Takahito Nakajima
- RIKEN Center for Computational Science, 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe 650-0047, Japan
| |
Collapse
|
42
|
Wei L, Chen Y, Liu J, Rao L, Ren Y, Xu X, Wan J. Cov_DOX: A Method for Structure Prediction of Covalent Protein-Ligand Bindings. J Med Chem 2022; 65:5528-5538. [PMID: 35353519 DOI: 10.1021/acs.jmedchem.1c02007] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A handful of molecular docking tools have been extended to enable a covalent docking. However, all of them face the challenge brought by the covalent bond between proteins and ligands. Many covalent drug design scenarios still heavily rely on demanding crystallographic experiments for accurate binding structures. Aiming at filling the gap between covalent dockings and crystallographic experiments, we develop and validate a hybrid method, dubbed as Cov_DOX, in this work. Cov_DOX achieves an overall success rate of 81% with RMSD < 2 Å for the Top 1 pose prediction in the validation against a test set including 405 crystal structures for covalent protein-ligand complexes, covering various types of the warhead chemistry and receptors. Such accuracy is not far from the much more demanding crystallographic experiments, in sharp contrast to the performance of the covalent docking front runners (success rate: 40-60%).
Collapse
Affiliation(s)
- Lin Wei
- Hubei International Scientific and Technological Cooperation Base of Pesticide and Green Synthesis, Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 43009, China
| | - Yaru Chen
- Hubei International Scientific and Technological Cooperation Base of Pesticide and Green Synthesis, Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 43009, China
| | - Jiaqi Liu
- Hubei International Scientific and Technological Cooperation Base of Pesticide and Green Synthesis, Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 43009, China
| | - Li Rao
- Hubei International Scientific and Technological Cooperation Base of Pesticide and Green Synthesis, Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 43009, China
| | - Yanliang Ren
- Hubei International Scientific and Technological Cooperation Base of Pesticide and Green Synthesis, Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 43009, China
| | - Xin Xu
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Ministry of Education (MOE) Laboratory for Computational Physical Science, Department of Chemistry, Fudan University, Shanghai 200433, People's Republic of China
| | - Jian Wan
- Hubei International Scientific and Technological Cooperation Base of Pesticide and Green Synthesis, Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 43009, China
| |
Collapse
|
43
|
López R, Díaz N, Francisco E, Martín-Pendás A, Suárez D. QM/MM Energy Decomposition Using the Interacting Quantum Atoms Approach. J Chem Inf Model 2022; 62:1510-1524. [PMID: 35212531 PMCID: PMC8965874 DOI: 10.1021/acs.jcim.1c01372] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The interacting quantum atoms (IQA) method decomposes the quantum mechanical (QM) energy of a molecular system in terms of one- and two-center (atomic) contributions within the context of the quantum theory of atoms in molecules. Here, we demonstrate that IQA, enhanced with molecular mechanics (MM) and Poisson-Boltzmann surface-area (PBSA) solvation methods, is naturally extended to the realm of hybrid QM/MM methodologies, yielding intra- and inter-residue energy terms that characterize all kinds of covalent and noncovalent bonding interactions. To test the robustness of this approach, both metal-water interactions and QM/MM boundary artifacts are characterized in terms of the IQA descriptors derived from QM regions of varying size in Zn(II)- and Mg(II)-water clusters. In addition, we analyze a homologous series of inhibitors in complex with a matrix metalloproteinase (MMP-12) by carrying out QM/MM-PBSA calculations on their crystallographic structures followed by IQA energy decomposition. Overall, these applications not only show the advantages of the IQA QM/MM approach but also address some of the challenges lying ahead for expanding the QM/MM methodology.
Collapse
Affiliation(s)
- Roberto López
- Departamento de Química y Física Aplicadas, Universidad de León, Facultad de Biología, Campus de Vegazana s/n, 24071 León (Castilla y León), Spain
| | - Natalia Díaz
- Departamento de Química Física y Analítica, Universidad de Oviedo, Facultad de Química, Julián Clavería 8, 33006 Oviedo (Asturias), Spain
| | - Evelio Francisco
- Departamento de Química Física y Analítica, Universidad de Oviedo, Facultad de Química, Julián Clavería 8, 33006 Oviedo (Asturias), Spain
| | - Angel Martín-Pendás
- Departamento de Química Física y Analítica, Universidad de Oviedo, Facultad de Química, Julián Clavería 8, 33006 Oviedo (Asturias), Spain
| | - Dimas Suárez
- Departamento de Química Física y Analítica, Universidad de Oviedo, Facultad de Química, Julián Clavería 8, 33006 Oviedo (Asturias), Spain
| |
Collapse
|
44
|
Trapl D, Krupička M, Višňovský V, Hozzová J, Ol'ha J, Křenek A, Spiwok V. Property Map Collective Variable as a Useful Tool for a Force Field Correction. J Chem Inf Model 2022; 62:567-576. [PMID: 35112877 DOI: 10.1021/acs.jcim.1c00651] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The accuracy of biomolecular simulations depends on the accuracy of an empirical molecular mechanics potential known as a force field: a set of parameters and expressions to estimate the potential from atomic coordinates. Accurate parametrization of force fields for small organic molecules is a challenge due to their high diversity. One of the possible approaches is to apply a correction to the existing force fields. Here, we propose an approach to estimate the density functional theory (DFT)-derived force field correction which is calculated during the run of molecular dynamics without significantly affecting its speed. Using the formula known as a property map collective variable, we approximate the force field correction by a weighted average of this force field correction calculated only for a small series of reference structures. To validate this method, we used seven AMBER force fields, and we show how it is possible to convert one force field to behave like the other one. We also present the force field correction for the important anticancer drug Imatinib as a use case example. Our method appears to be suitable for adjusting the force field for general drug-like molecules. We provide a pipeline that generates the correction; this pipeline is available at https://pmcvff-correction.cerit-sc.cz/.
Collapse
Affiliation(s)
- Dalibor Trapl
- Department of Biochemistry and Microbiology, University of Chemistry and Technology, Technická 5, Prague 6 166 28, Czech Republic
| | - Martin Krupička
- Department of Organic Chemistry, University of Chemistry and Technology, Technická 5, Prague 6 166 28, Czech Republic
| | - Vladimír Višňovský
- Institute of Computer Science, Masaryk University, Botanická 554/68a, Brno 602 00, Czech Republic
| | - Jana Hozzová
- Institute of Computer Science, Masaryk University, Botanická 554/68a, Brno 602 00, Czech Republic
| | - Jaroslav Ol'ha
- Institute of Computer Science, Masaryk University, Botanická 554/68a, Brno 602 00, Czech Republic
| | - Aleš Křenek
- Institute of Computer Science, Masaryk University, Botanická 554/68a, Brno 602 00, Czech Republic
| | - Vojtěch Spiwok
- Department of Biochemistry and Microbiology, University of Chemistry and Technology, Technická 5, Prague 6 166 28, Czech Republic
| |
Collapse
|
45
|
Qu X, Dong L, Si Y, Zhao Y, Wang Q, Su P, Wang B. Reliable Prediction of the Protein-Ligand Binding Affinity Using a Charge Penetration Corrected AMOEBA Force Field: A Case Study of Drug Resistance Mutations in Abl Kinase. J Chem Theory Comput 2022; 18:1692-1700. [PMID: 35107298 DOI: 10.1021/acs.jctc.1c01005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Protein mutations that directly impair drug binding are related to therapeutic resistance, and accurate prediction of their impact on drug binding would benefit drug design and clinical practice. Here, we have developed a scoring strategy that predicts the effect of the mutations on the protein-ligand binding affinity. In view of the critical importance of electrostatics in protein-ligand interactions, the charge penetration corrected AMOEBA force field (AMOEBA_CP model) was employed to improve the accuracy of the calculated electrostatic energy. We calculated the electrostatic energy using an energy decomposition analysis scheme based on the generalized Kohn-Sham (GKS-EDA). The AMOEBA_CP model was validated by a protein-fragment-ligand complex data set (Abl236) constructed from the co-crystal structures of the cancer target Abl kinase with six inhibitors. To predict ligand binding affinity changes upon protein mutation of Abl kinase, we used sampling protocol with multistep simulated annealing to search conformations of mutant proteins. The scoring strategy based on AMOEBA_CP model has achieved considerable performance in predicting resistance for 8 kinase inhibitors across 144 clinically identified point mutations. Overall, this study illustrates that the AMOEBA_CP model, which accurately treats electrostatics through penetration correction, enables the accurate prediction of the mutation-induced variation of protein-ligand binding affinity.
Collapse
Affiliation(s)
- Xiaoyang Qu
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, P. R. China
| | - Lina Dong
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, P. R. China
| | - Yubing Si
- College of Chemistry, Zhengzhou University, Zhengzhou 450001, P. R. China
| | - Yuan Zhao
- The Key Laboratory of Natural Medicine and Immuno-Engineering, Henan University, Kaifeng 475004, P. R. China
| | - Qiantao Wang
- Key Laboratory of Drug-Targeting and Drug Delivery System of the Education Ministry and Sichuan Province, Sichuan Engineering Laboratory for Plant-Sourced Drug and Sichuan Research Center for Drug Precision Industrial Technology, West China School of Pharmacy, Sichuan University, Chengdu 610041, P. R. China
| | - Peifeng Su
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, P. R. China
| | - Binju Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, P. R. China
| |
Collapse
|
46
|
Decomposition of the interaction energy of several flavonoids with Escherichia coli DNA Gyr using the SAPT (DFT) method: The relation between the interaction energy components, ligand structure, and biological activity. Biochim Biophys Acta Gen Subj 2022; 1866:130111. [DOI: 10.1016/j.bbagen.2022.130111] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2021] [Revised: 01/19/2022] [Accepted: 02/07/2022] [Indexed: 12/28/2022]
|
47
|
Umbrella Sampling-Based Method to Compute Ligand-Binding Affinity. METHODS IN MOLECULAR BIOLOGY (CLIFTON, N.J.) 2022; 2385:313-323. [PMID: 34888726 DOI: 10.1007/978-1-0716-1767-0_14] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Many proteins have a solvent-exposed binding cleft, which permits their inhibitors to bind and unbind without significant protein conformation transforms. The binding/unbinding pathways of these protein-inhibitor complexes can be rather straightforwardly sampled by using umbrella sampling (US) simulation methods. During a US simulation, the Cα atoms of the protein are restrained via a harmonic force. The potential of mean force (PMF) along the binding pathway can be estimated by using the weighted histogram analysis method (WHAM). The binding affinity is then computed as the difference in PMF between the binding and unbinding states.
Collapse
|
48
|
Penaloza-Amion M, C Rêgo CR, Wenzel W. Local Electronic Charge Transfer in the Helical Induction of Cis-Transoid Poly(4-carboxyphenyl)acetylene by Chiral Amines. J Chem Inf Model 2022; 62:544-552. [PMID: 35080886 DOI: 10.1021/acs.jcim.1c01347] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Understanding the phenomena that lead to the formation of a specific helicity in helical polymers remains a challenge even today. Various polymers have been shown to assume different helical screw-senses depending on different stimuli. Acid-base chiral amines, for example, can induce helical conformations on cis-transoid poly(4-carboxyphenyl)acetylene yielding high-intensity circular dichroism signals. There have been many experimental attempts to elucidate the driving forces involved, but the induction process remains unclear. Here, we investigate the mechanism of helical polymer formation by both Molecular Dynamics (MD) and Density Functional Theory (DFT) approaches. We find that DFT calculations and the dissociation energies between 4 monomer polymers and amines show a clear trend in the affinity of R and S conformers with clockwise and counterclockwise polymer screw-senses, respectively. The charge analysis revealed that the local charge transfer effect plays a crucial role that leads to the helical polymer-amine induction.
Collapse
Affiliation(s)
- Montserrat Penaloza-Amion
- Institute of Nanotechnology Hermann-von-Helmholtz-Platz, Karlsruhe Institute of Technology, 76021 Karlsruhe, Germany
| | - Celso R C Rêgo
- Institute of Nanotechnology Hermann-von-Helmholtz-Platz, Karlsruhe Institute of Technology, 76021 Karlsruhe, Germany
| | - Wolfgang Wenzel
- Institute of Nanotechnology Hermann-von-Helmholtz-Platz, Karlsruhe Institute of Technology, 76021 Karlsruhe, Germany
| |
Collapse
|
49
|
Gundelach L, Fox T, Tautermann CS, Skylaris CK. BRD4: quantum mechanical protein–ligand binding free energies using the full-protein DFT-based QM-PBSA method. Phys Chem Chem Phys 2022; 24:25240-25249. [PMID: 36222107 DOI: 10.1039/d2cp03705j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Fully quantum mechanical approaches to calculating protein–ligand free energies of binding have the potential to reduce empiricism and explicitly account for all physical interactions responsible for protein–ligand binding.
Collapse
Affiliation(s)
- Lennart Gundelach
- University of Southampton, Faculty of Engineering Science and Mathematics, Chemistry, University Road, Southampton, SO17 1BJ, UK
| | - Thomas Fox
- Boehringer Ingelheim Pharma GmbH & Co KG, Medicinal Chemistry, Birkendorfer Str 65, 88397, Biberach, Germany
| | - Christofer S. Tautermann
- Boehringer Ingelheim Pharma GmbH & Co KG, Medicinal Chemistry, Birkendorfer Str 65, 88397, Biberach, Germany
| | - Chris-Kriton Skylaris
- University of Southampton, Faculty of Engineering Science and Mathematics, Chemistry, University Road, Southampton, SO17 1BJ, UK
| |
Collapse
|
50
|
Pham TNH, Nguyen TH, Tam NM, Y Vu T, Pham NT, Huy NT, Mai BK, Tung NT, Pham MQ, V Vu V, Ngo ST. Improving ligand-ranking of AutoDock Vina by changing the empirical parameters. J Comput Chem 2021; 43:160-169. [PMID: 34716930 DOI: 10.1002/jcc.26779] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2021] [Revised: 10/10/2021] [Accepted: 10/14/2021] [Indexed: 01/09/2023]
Abstract
AutoDock Vina (Vina) achieved a very high docking-success rate, p ^ , but give a rather low correlation coefficient, R , for binding affinity with respect to experiments. This low correlation can be an obstacle for ranking of ligand-binding affinity, which is the main objective of docking simulations. In this context, we evaluated the dependence of Vina R coefficient upon its empirical parameters. R is affected more by changing the gauss2 and rotation than other terms. The docking-success rate p ^ is sensitive to the alterations of the gauss1, gauss2, repulsion, and hydrogen bond parameters. Based on our benchmarks, the parameter set1 has been suggested to be the most optimal. The testing study over 800 complexes indicated that the modified Vina provided higher correlation with experiment R set 1 = 0.556 ± 0.025 compared with R Default = 0.493 ± 0.028 obtained by the original Vina and R Vina 1.2 = 0.503 ± 0.029 by Vina version 1.2. Besides, the modified Vina can be also applied more widely, giving R ≥ 0.500 for 32/48 targets, compared with the default package, giving R ≥ 0.500 for 31/48 targets. In addition, validation calculations for 1036 complexes obtained from version 2019 of PDBbind refined structures showed that the set1 of parameters gave higher correlation coefficient ( R set 1 = 0.617 ± 0.017 ) than the default package ( R Default = 0.543 ± 0.020 ) and Vina version 1.2 ( R Vina 1.2 = 0.540 ± 0.020 ). The version of Vina with set1 of parameters can be downloaded at https://github.com/sontungngo/mvina. The outcomes would enhance the ranking of ligand-binding affinity using Autodock Vina.
Collapse
Affiliation(s)
- T Ngoc Han Pham
- Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Trung Hai Nguyen
- Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City, Vietnam.,Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Nguyen Minh Tam
- Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam.,Computational Chemistry Research Group, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Thien Y Vu
- Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Nhat Truong Pham
- Faculty of Electrical and Electronics Engineering, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Nguyen Truong Huy
- Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Binh Khanh Mai
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania, USA
| | - Nguyen Thanh Tung
- Institute of Materials Science, Vietnam Academy of Science and Technology, Hanoi, Vietnam.,Graduate University of Science and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Minh Quan Pham
- Graduate University of Science and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam.,Institute of Natural Products Chemistry, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Van V Vu
- NTT Hi-Tech Institute, Nguyen Tat Thanh University, Ho Chi Minh City, Vietnam
| | - Son Tung Ngo
- Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City, Vietnam.,Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| |
Collapse
|