51
|
Yuan C, Goonetilleke EC, Unarta IC, Huang X. Incorporation efficiency and inhibition mechanism of 2'-substituted nucleotide analogs against SARS-CoV-2 RNA-dependent RNA polymerase. Phys Chem Chem Phys 2021; 23:20117-20128. [PMID: 34514487 DOI: 10.1039/d1cp03049c] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
The ongoing pandemic caused by SARS-CoV-2 emphasizes the need for effective therapeutics. Inhibition of SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) by nucleotide analogs provides a promising antiviral strategy. One common group of RdRp inhibitors, 2'-modified nucleotides, are reported to exhibit different behaviors in the SARS-CoV-2 RdRp transcription assay. Three of these analogs, 2'-O-methyl UTP, Sofosbuvir, and 2'-methyl CTP, act as effective inhibitors in previous biochemical experiments, while Gemcitabine and ara-UTP show no inhibitory activity. To understand the impact of the 2'-modification on their inhibitory effects, we conducted extensive molecular dynamics simulations and relative binding free energy calculations using the free energy perturbation method on SARS-CoV-2 replication-transcription complex (RTC) with these five nucleotide analogs. Our results reveal that the five nucleotide analogs display comparable binding affinities to SARS-CoV-2 RdRp and they can all be added to the nascent RNA chain. Moreover, we examine how the incorporation of these nucleotide triphosphate (NTP) analogs will impact the addition of the next nucleotide. Our results indicate that 2'-O-methyl UTP can weaken the binding of the subsequent NTP and consequently lead to partial chain termination. Additionally, Sofosbuvir and 2'-methyl CTP can cause immediate termination due to the strong steric hindrance introduced by their bulky 2'-methyl groups. In contrast, nucleotide analogs with smaller substitutions, such as the fluorine atoms and the ara-hydroxyl group in Gemcitabine and ara-UTP, have a marginal impact on the polymerization process. Our findings are consistent with experimental observations, and more importantly, shed light on the detailed molecular mechanism of SARS-CoV-2 RdRp inhibition by 2'-substituted nucleotide analogs, and may facilitate the rational design of antiviral agents to inhibit SARS-CoV-2 RdRp.
Collapse
Affiliation(s)
- Congmin Yuan
- Department of Chemistry, Centre of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong. .,Hong Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| | - Eshani C Goonetilleke
- Department of Chemistry, Centre of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong. .,Hong Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| | - Ilona Christy Unarta
- Department of Chemistry, Centre of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong. .,Hong Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| | - Xuhui Huang
- Department of Chemistry, Centre of Systems Biology and Human Health, State Key Laboratory of Molecular Neuroscience, The Hong Kong University of Science and Technology, Kowloon, Hong Kong. .,Hong Kong Center for Neurodegenerative Diseases, Hong Kong Science Park, Hong Kong
| |
Collapse
|
52
|
King E, Aitchison E, Li H, Luo R. Recent Developments in Free Energy Calculations for Drug Discovery. Front Mol Biosci 2021; 8:712085. [PMID: 34458321 PMCID: PMC8387144 DOI: 10.3389/fmolb.2021.712085] [Citation(s) in RCA: 49] [Impact Index Per Article: 16.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2021] [Accepted: 07/27/2021] [Indexed: 01/11/2023] Open
Abstract
The grand challenge in structure-based drug design is achieving accurate prediction of binding free energies. Molecular dynamics (MD) simulations enable modeling of conformational changes critical to the binding process, leading to calculation of thermodynamic quantities involved in estimation of binding affinities. With recent advancements in computing capability and predictive accuracy, MD based virtual screening has progressed from the domain of theoretical attempts to real application in drug development. Approaches including the Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA), Linear Interaction Energy (LIE), and alchemical methods have been broadly applied to model molecular recognition for drug discovery and lead optimization. Here we review the varied methodology of these approaches, developments enhancing simulation efficiency and reliability, remaining challenges hindering predictive performance, and applications to problems in the fields of medicine and biochemistry.
Collapse
Affiliation(s)
- Edward King
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA, United States
| | - Erick Aitchison
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA, United States
| | - Han Li
- Department of Chemical and Biomolecular Engineering, University of California, Irvine, CA, United States
| | - Ray Luo
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA, United States
- Department of Chemical and Biomolecular Engineering, University of California, Irvine, CA, United States
- Department of Materials Science and Engineering, University of California, Irvine, CA, United States
- Department of Biomedical Engineering, University of California, Irvine, CA, United States
| |
Collapse
|
53
|
Guterres H, Park SJ, Zhang H, Im W. CHARMM-GUI LBS Finder & Refiner for Ligand Binding Site Prediction and Refinement. J Chem Inf Model 2021; 61:3744-3751. [PMID: 34296608 DOI: 10.1021/acs.jcim.1c00561] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
A protein performs its task by binding a variety of ligands in its local region that is also known as the ligand-binding-site (LBS). Therefore, accurate prediction, characterization, and refinement of LBS can facilitate protein functional annotations and structure-based drug design. In this work, we present CHARMM-GUI LBS Finder & Refiner (https://www.charmm-gui.org/input/lbsfinder) that predicts potential LBS, offers interactive features for local LBS structure analysis, and prepares various molecular dynamics (MD) systems and inputs by setting up distance restraint potentials for LBS structure refinement. LBS Finder & Refiner supports 5 different commonly used simulation programs, such as NAMD, AMBER, GROMACS, GENESIS, and OpenMM, for LBS structure refinement together with hydrogen mass repartitioning. The capability of LBS Finder & Refiner is illustrated through LBS structure predictions and refinements of 48 modeled and 20 apo benchmark target proteins. Overall, successful LBS structure predictions and refinements are seen in our benchmark tests. We hope that LBS Finder & Refiner is useful to predict, characterize, and refine potential LBS on any given protein of interest.
Collapse
Affiliation(s)
- Hugo Guterres
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Sang-Jun Park
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Han Zhang
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Wonpil Im
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| |
Collapse
|
54
|
Salih AEM, Thissera B, Yaseen M, Hassane ASI, El-Seedi HR, Sayed AM, Rateb ME. Marine Sulfated Polysaccharides as Promising Antiviral Agents: A Comprehensive Report and Modeling Study Focusing on SARS CoV-2. Mar Drugs 2021; 19:406. [PMID: 34436245 PMCID: PMC8401819 DOI: 10.3390/md19080406] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2021] [Revised: 07/18/2021] [Accepted: 07/20/2021] [Indexed: 11/17/2022] Open
Abstract
SARS-CoV-2 (severe acute respiratory syndrome coronavirus-2) is a novel coronavirus strain that emerged at the end of 2019, causing millions of deaths so far. Despite enormous efforts being made through various drug discovery campaigns, there is still a desperate need for treatments with high efficacy and selectivity. Recently, marine sulfated polysaccharides (MSPs) have earned significant attention and are widely examined against many viral infections. This article attempted to produce a comprehensive report about MSPs from different marine sources alongside their antiviral effects against various viral species covering the last 25 years of research articles. Additionally, these reported MSPs were subjected to molecular docking and dynamic simulation experiments to ascertain potential interactions with both the receptor-binding domain (RBD) of SARS CoV-2's spike protein (S-protein) and human angiotensin-converting enzyme-2 (ACE2). The possible binding sites on both S-protein's RBD and ACE2 were determined based on how they bind to heparin, which has been reported to exhibit significant antiviral activity against SARS CoV-2 through binding to RBD, preventing the virus from affecting ACE2. Moreover, our modeling results illustrate that heparin can also bind to and block ACE2, acting as a competitor and protective agent against SARS CoV-2 infection. Nine of the investigated MSPs candidates exhibited promising results, taking into consideration the newly emerged SARS CoV-2 variants, of which five were not previously reported to exert antiviral activity against SARS CoV-2, including sulfated galactofucan (1), sulfated polymannuroguluronate (SPMG) (2), sulfated mannan (3), sulfated heterorhamnan (8), and chondroitin sulfate E (CS-E) (9). These results shed light on the importance of sulfated polysaccharides as potential SARS-CoV-2 inhibitors.
Collapse
Affiliation(s)
- Abdalla E. M. Salih
- School of Computing, Engineering & Physical Sciences, University of the West of Scotland, Paisley PA1 2BE, UK; (A.E.M.S.); (B.T.); (M.Y.); (A.S.I.H.)
| | - Bathini Thissera
- School of Computing, Engineering & Physical Sciences, University of the West of Scotland, Paisley PA1 2BE, UK; (A.E.M.S.); (B.T.); (M.Y.); (A.S.I.H.)
| | - Mohammed Yaseen
- School of Computing, Engineering & Physical Sciences, University of the West of Scotland, Paisley PA1 2BE, UK; (A.E.M.S.); (B.T.); (M.Y.); (A.S.I.H.)
| | - Ahmed S. I. Hassane
- School of Computing, Engineering & Physical Sciences, University of the West of Scotland, Paisley PA1 2BE, UK; (A.E.M.S.); (B.T.); (M.Y.); (A.S.I.H.)
- Aberdeen Royal Infirmary, Foresterhill Health Campus, Aberdeen AB25 2ZN, UK
| | - Hesham R. El-Seedi
- Pharmacognosy Group, Department of Pharmaceutical Biosciences, BMC, Uppsala University, Uppsala, Box 591, SE 751 24 Uppsala, Sweden;
- International Research Center for Food Nutrition and Safety, Jiangsu University, Zhenjiang 212013, China
- Department of Chemistry, Faculty of Science, Menoufia University, Shebin El-Kom 32512, Egypt
| | - Ahmed M. Sayed
- Department of Pharmacognosy, Faculty of Pharmacy, Nahda University, Beni-Suef 62513, Egypt
| | - Mostafa E. Rateb
- School of Computing, Engineering & Physical Sciences, University of the West of Scotland, Paisley PA1 2BE, UK; (A.E.M.S.); (B.T.); (M.Y.); (A.S.I.H.)
| |
Collapse
|
55
|
El-Hawary SS, Sayed AM, Issa MY, Ebrahim HS, Alaaeldin R, Elrehany MA, Abd El-Kadder EM, Abdelmohsen UR. Anti-Alzheimer chemical constituents of Morus macroura Miq.: chemical profiling, in silico and in vitro investigations. Food Funct 2021; 12:8078-8089. [PMID: 34286787 DOI: 10.1039/d1fo01177d] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Herein, we investigated both fruits and leaves of Morus macroura Miq. as a potential source of bioactive compounds against Alzheimer's disease (AD). LC-HRMS-assisted chemical profiling of its extracts showed that they are a rich source of diverse phytochemicals. Among the 29 identified compounds in both the fruit and leaf extracts, moracin D, chrysin, resveratrol, and ferulic acid were predicted to pass the human blood-brain barrier (BBB), and hence, reach their therapeutic targets in the brain. Subsequently, these compounds were subjected to a comprehensive pharmacophore-based screening for their protein targets relevant to AD using two independent software programs (i.e. Swiss Target Prediction and PharmMapper). The results of this initial virtual screening were further refined by a number of docking and molecular dynamic simulation experiments to suggest a number of crucial AD-related proteins (e.g. acetylcholine esterase, β-secretase, and monoamine oxidase) as potential targets for these compounds. Finally, in vitro testing was performed to validate the in silico investigation's results, where chrysin, resveratrol, and ferulic acid were found to inhibit the predicted AD-related enzymes with IC50 values comparable with those of the reference inhibitors. Additionally, they were able to inhibit the aggregation of amyloid-beta, one of the hallmarks in AD pathogenesis, and to exhibit considerable antioxidant capacity. Our results highlighted Morus macroura compounds as future anti-Alzheimer chemical leads.
Collapse
Affiliation(s)
- Seham S El-Hawary
- Pharmacognosy Department, Faculty of Pharmacy, Cairo University, Cairo, Egypt
| | | | | | | | | | | | | | | |
Collapse
|
56
|
Alhadrami HA, Sayed AM, Melebari SA, Khogeer AA, Abdulaal WH, Al-Fageeh MB, Algahtani M, Rateb ME. Targeting allosteric sites of human aromatase: a comprehensive in-silico and in-vitro workflow to find potential plant-based anti-breast cancer therapeutics. J Enzyme Inhib Med Chem 2021; 36:1334-1345. [PMID: 34139914 PMCID: PMC8759730 DOI: 10.1080/14756366.2021.1937145] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Recent findings suggested several allosteric pockets on human aromatase that could be utilised for the development of new modulators able to inhibit this enzyme in a new mechanism. Herein, we applied an integrated in-silico-based approach supported by in-vitro enzyme-based and cell-based validation assays to select the best leads able to target these allosteric binding sites from a small library of plant-derived natural products. Chrysin, apigenin, and resveratrol were found to be the best inhibitors targeting the enzyme’s substrate access channel and were able to produce a competitive inhibition with IC50 values ranged from 1.7 to 15.8 µM. Moreover, they showed a more potent antiproliferative effect against ER+ (MCF-7) than ER- one (MDA-MB-231) cell lines. On the other hand, both pomiferin and berberine were the best hits for the enzyme’s haem-proximal cavity producing a non-competitive inhibition (IC50 15.1 and 21.4 µM, respectively) and showed selective antiproliferative activity towards MCF-7 cell lines.
Collapse
Affiliation(s)
- Hani A Alhadrami
- Department of Medical Laboratory Technology, Faculty of Applied Medical Sciences, King Abdulaziz University, Jeddah, Saudi Arabia.,Molecular Diagnostic Lab, King Abdulaziz University Hospital, King Abdulaziz University, Jeddah, Saudi Arabia.,Molecular Diagnostic Unit, The Regional Laboratory in Makkah, Ministry of Health, Makkah, Kingdom of Saudi Arabia
| | - Ahmed M Sayed
- Department of Pharmacognosy, Faculty of Pharmacy, Nahda University, Beni Suef, Egypt
| | - Sami A Melebari
- Molecular Diagnostic Unit, The Regional Laboratory in Makkah, Ministry of Health, Makkah, Kingdom of Saudi Arabia
| | - Asem A Khogeer
- Plan and Research Department, General Directorate of Health Affairs, Makkah region, Ministry of Health, Makkah, Kingdom of Saudi Arabia
| | - Wesam H Abdulaal
- Department of Biochemistry, Faculty of Science, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Mohamed B Al-Fageeh
- General Directorate for Funds and Grants (GDFG), King Abdulaziz City for Science and Technology, Riyadh, Saudi Arabia
| | - Mohammad Algahtani
- Department of Laboratory and Blood Bank, Security Forces Hospital Program, Mecca, Saudi Arabia
| | - Mostafa E Rateb
- School of Computing, Engineering & Physical Sciences, University of the West of Scotland, Paisley, UK
| |
Collapse
|
57
|
Alhadrami HA, Sayed AM, Al-Khatabi H, Alhakamy NA, Rateb ME. Scaffold Hopping of α-Rubromycin Enables Direct Access to FDA-Approved Cromoglicic Acid as a SARS-CoV-2 M Pro Inhibitor. Pharmaceuticals (Basel) 2021; 14:541. [PMID: 34198933 PMCID: PMC8229550 DOI: 10.3390/ph14060541] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2021] [Revised: 06/01/2021] [Accepted: 06/03/2021] [Indexed: 12/23/2022] Open
Abstract
The COVID-19 pandemic is still active around the globe despite the newly introduced vaccines. Hence, finding effective medications or repurposing available ones could offer great help during this serious situation. During our anti-COVID-19 investigation of microbial natural products (MNPs), we came across α-rubromycin, an antibiotic derived from Streptomyces collinus ATCC19743, which was able to suppress the catalytic activity (IC50 = 5.4 µM and Ki = 3.22 µM) of one of the viral key enzymes (i.e., MPro). However, it showed high cytotoxicity toward normal human fibroblasts (CC50 = 16.7 µM). To reduce the cytotoxicity of this microbial metabolite, we utilized a number of in silico tools (ensemble docking, molecular dynamics simulation, binding free energy calculation) to propose a novel scaffold having the main pharmacophoric features to inhibit MPro with better drug-like properties and reduced/minimal toxicity. Nevertheless, reaching this novel scaffold synthetically is a time-consuming process, particularly at this critical time. Instead, this scaffold was used as a template to explore similar molecules among the FDA-approved medications that share its main pharmacophoric features with the aid of pharmacophore-based virtual screening software. As a result, cromoglicic acid (aka cromolyn) was found to be the best hit, which, upon in vitro MPro testing, was 4.5 times more potent (IC50 = 1.1 µM and Ki = 0.68 µM) than α-rubromycin, with minimal cytotoxicity toward normal human fibroblasts (CC50 > 100 µM). This report highlights the potential of MNPs in providing unprecedented scaffolds with a wide range of therapeutic efficacy. It also revealed the importance of cheminformatics tools in speeding up the drug discovery process, which is extremely important in such a critical situation.
Collapse
Affiliation(s)
- Hani A. Alhadrami
- Department of Medical Laboratory Technology, Faculty of Applied Medical Sciences, King Abdulaziz University, P.O. Box 80402, Jeddah 21589, Saudi Arabia; (H.A.A.); (H.A.-K.)
- Molecular Diagnostic Lab, King Abdulaziz University Hospital, King Abdulaziz University, P.O. Box 80402, Jeddah 21589, Saudi Arabia
| | - Ahmed M. Sayed
- Department of Pharmacognosy, Faculty of Pharmacy, Nahda University, Beni-Suef 62513, Egypt;
| | - Heba Al-Khatabi
- Department of Medical Laboratory Technology, Faculty of Applied Medical Sciences, King Abdulaziz University, P.O. Box 80402, Jeddah 21589, Saudi Arabia; (H.A.A.); (H.A.-K.)
- Center of Excellence in Genomic Medicine Research, King Abdulaziz University, Jeddah 21589, Saudi Arabia
| | - Nabil A. Alhakamy
- Department of Pharmaceutics, Faculty of Pharmacy, King Abdulaziz University, Jeddah 21589, Saudi Arabia;
| | - Mostafa E. Rateb
- School of Computing, Engineering & Physical Sciences, University of the West of Scotland, Paisley PA1 2BE, UK
| |
Collapse
|
58
|
Yagi K, Ito S, Sugita Y. Exploring the Minimum-Energy Pathways and Free-Energy Profiles of Enzymatic Reactions with QM/MM Calculations. J Phys Chem B 2021; 125:4701-4713. [PMID: 33914537 PMCID: PMC10986901 DOI: 10.1021/acs.jpcb.1c01862] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Understanding molecular mechanisms of enzymatic reactions is of vital importance in biochemistry and biophysics. Here, we introduce new functions of hybrid quantum mechanical/molecular mechanical (QM/MM) calculations in the GENESIS program to compute the minimum-energy pathways (MEPs) and free-energy profiles of enzymatic reactions. For this purpose, an interface in GENESIS is developed to utilize a highly parallel electronic structure program, QSimulate-QM (https://qsimulate.com), calling it as a shared library from GENESIS. Second, algorithms to search the MEP are implemented, combining the string method (E et al. J. Chem. Phys. 2007, 126, 164103) with the energy minimization of the buffer MM region. The method implemented in GENESIS is applied to an enzyme, triosephosphate isomerase, which converts dihyroxyacetone phosphate to glyceraldehyde 3-phosphate in four proton-transfer processes. QM/MM-molecular dynamics simulations show performances of greater than 1 ns/day with the density functional tight binding (DFTB), and 10-30 ps/day with the hybrid density functional theory, B3LYP-D3. These performances allow us to compute not only MEP but also the potential of mean force (PMF) of the enzymatic reactions using the QM/MM calculations. The barrier height obtained as 13 kcal mol-1 with B3LYP-D3 in the QM/MM calculation is in agreement with the experimental results. The impact of conformational sampling in PMF calculations and the level of electronic structure calculations (DFTB vs B3LYP-D3) suggests reliable computational protocols for enzymatic reactions without high computational costs.
Collapse
Affiliation(s)
- Kiyoshi Yagi
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Shingo Ito
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Yuji Sugita
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, 7-1-26 minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, 1-6-5 minatojima-Minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
59
|
Fu H, Chen H, Cai W, Shao X, Chipot C. BFEE2: Automated, Streamlined, and Accurate Absolute Binding Free-Energy Calculations. J Chem Inf Model 2021; 61:2116-2123. [PMID: 33906354 DOI: 10.1021/acs.jcim.1c00269] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Accurate absolute binding free-energy estimation in silico, following either an alchemical or a geometrical route, involves several subprocesses and requires the introduction of geometric restraints. Human intervention, for instance, to define the necessary collective variables, prepare the input files, monitor the simulation, and perform post-treatments is, however, tedious, cumbersome, and prone to errors. With the aim of automating and streamlining free-energy calculations, especially for nonexperts, version 2.0 of the binding free energy estimator (BFEE2) provides both standardized alchemical and geometrical workflows and obviates the need for extensive human intervention to guarantee complete reproducibility of the results. To achieve the largest gamut of protein-ligand and, more generally, of host-guest complexes, BFEE2 supports most academic force fields, such as CHARMM, Amber, OPLS, and GROMOS. Configurational files are generated in the NAMD and Gromacs formats, and all the post-treatments are performed in an automated fashion. Moreover, convergence of the free-energy calculation can be monitored from the intermediate files generated during the simulation. All in all, BFEE2 is a foolproof, versatile tool for accurate absolute binding free-energy calculations, assisting the end-user over a broad range of applications.
Collapse
Affiliation(s)
- Haohao Fu
- Research Center for Analytical Sciences, College of Chemistry, Nankai University, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, Tianjin 300071, China
| | - Haochuan Chen
- Research Center for Analytical Sciences, College of Chemistry, Nankai University, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, Tianjin 300071, China
| | - Wensheng Cai
- Research Center for Analytical Sciences, College of Chemistry, Nankai University, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, Tianjin 300071, China
| | - Xueguang Shao
- Research Center for Analytical Sciences, College of Chemistry, Nankai University, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, Tianjin 300071, China
| | - Christophe Chipot
- Laboratoire International Associé CNRS and University of Illinois at Urbana-Champaign, UMR n°7019, Université de Lorraine, BP 70239, F-54506 Vandœuvre-lès-Nancy, France.,Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois 61801, United States.,Theoretical and Computational Biophysics Group, Beckman Institute for Advanced Science and Technology, University of Illinois at Urban-Champaign, 405 North Mathews, Urbana, Illinois 61801, United States
| |
Collapse
|
60
|
Ren W, Dokainish HM, Shinobu A, Oshima H, Sugita Y. Unraveling the Coupling between Conformational Changes and Ligand Binding in Ribose Binding Protein Using Multiscale Molecular Dynamics and Free-Energy Calculations. J Phys Chem B 2021; 125:2898-2909. [PMID: 33728914 PMCID: PMC10954230 DOI: 10.1021/acs.jpcb.0c11600] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Conformational changes of proteins upon ligand binding are usually explained in terms of several mechanisms including the induced fit, conformational selection, or their mixtures. Due to the slow time scales, conventional molecular dynamics (cMD) simulations based on the atomistic models cannot easily simulate the open-to-closed conformational transition in proteins. In our previous study, we have developed an enhanced sampling scheme (generalized replica exchange with solute tempering selected surface charged residues: gREST_SSCR) for multidomain proteins and applied it to ligand-mediated conformational changes in the G134R mutant of ribose-binding protein (RBPG134R) in solution. The free-energy landscape (FEL) of RBPG134R in the presence of a ribose at the binding site included the open and closed states and two intermediates, open-like and closed-like forms. Only the open and open-like forms existed in the FEL without a ribose. In the current study, the coupling between the conformational changes and ligand binding is further investigated using coarse-grained MD, multiple atomistic cMD, and free-energy calculations. The ribose is easily dissociated from the binding site of wild-type RBP and RBPG134R in the cMD simulations starting from the open and open-like forms. In contrast, it is stable at the binding site in the simulations from the closed and closed-like forms. The free-energy calculations provide the binding affinities of different structures, supporting the results of cMD simulations. Importantly, cMD simulations from the closed-like structures reveal transitions toward the closed one in the presence of a bound ribose. On the basis of the computational results, we propose a molecular mechanism in which conformational selection and induced fit happen in the first and second halves of the open-to-closed transition in RBP, respectively.
Collapse
Affiliation(s)
- Weitong Ren
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Hisham M. Dokainish
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Ai Shinobu
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, Integrated Innovation Building 7F, 6-7-1 minatojima-minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Hiraku Oshima
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, Integrated Innovation Building 7F, 6-7-1 minatojima-minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Yuji Sugita
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Integrated Innovation Building 7F, 6-7-1 minatojima-minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, Integrated Innovation Building 7F, 6-7-1 minatojima-minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
61
|
Heinzelmann G, Gilson MK. Automation of absolute protein-ligand binding free energy calculations for docking refinement and compound evaluation. Sci Rep 2021; 11:1116. [PMID: 33441879 PMCID: PMC7806944 DOI: 10.1038/s41598-020-80769-1] [Citation(s) in RCA: 40] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2020] [Accepted: 12/24/2020] [Indexed: 02/06/2023] Open
Abstract
Absolute binding free energy calculations with explicit solvent molecular simulations can provide estimates of protein-ligand affinities, and thus reduce the time and costs needed to find new drug candidates. However, these calculations can be complex to implement and perform. Here, we introduce the software BAT.py, a Python tool that invokes the AMBER simulation package to automate the calculation of binding free energies for a protein with a series of ligands. The software supports the attach-pull-release (APR) and double decoupling (DD) binding free energy methods, as well as the simultaneous decoupling-recoupling (SDR) method, a variant of double decoupling that avoids numerical artifacts associated with charged ligands. We report encouraging initial test applications of this software both to re-rank docked poses and to estimate overall binding free energies. We also show that it is practical to carry out these calculations cheaply by using graphical processing units in common machines that can be built for this purpose. The combination of automation and low cost positions this procedure to be applied in a relatively high-throughput mode and thus stands to enable new applications in early-stage drug discovery.
Collapse
Affiliation(s)
- Germano Heinzelmann
- Departamento de Física, Universidade Federal de Santa Catarina, Florianópolis, Santa Catarina, Brazil.
| | - Michael K Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, USA
| |
Collapse
|
62
|
Guterres H, Park SJ, Jiang W, Im W. Ligand-Binding-Site Refinement to Generate Reliable Holo Protein Structure Conformations from Apo Structures. J Chem Inf Model 2020; 61:535-546. [PMID: 33337877 DOI: 10.1021/acs.jcim.0c01354] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The first important step in a structure-based virtual screening is the judicious selection of a receptor protein. In cases where the holo protein receptor structure is unavailable, significant reduction in virtual screening performance has been reported. In this work, we present a robust method to generate reliable holo protein structure conformations from apo structures using molecular dynamics (MD) simulation with restraints derived from holo structure binding-site templates. We perform benchmark tests on two different datasets: 40 structures from a directory of useful decoy-enhanced (DUD-E) and 84 structures from the Gunasekaran dataset. Our results show successful refinement of apo binding-site structures toward holo conformations in 82% of the test cases. In addition, virtual screening performance of 40 DUD-E structures is significantly improved using our MD-refined structures as receptors with an average enrichment factor (EF), an EF1% value of 6.2 compared to apo structures with 3.5. Docking of native ligands to the refined structures shows an average ligand root mean square deviation (RMSD) of 1.97 Å (DUD-E dataset and Gunasekaran dataset) relative to ligands in the holo crystal structures, which is comparable to the self-docking (i.e., docking of the native ligand back to its crystal structure receptor) average, 1.34 Å (DUD-E dataset) and 1.36 Å (Gunasekaran dataset). On the other hand, docking to the apo structures yields an average ligand RMSD of 3.65 Å (DUD-E) and 2.90 Å (Gunasekaran). These results indicate that our method is robust and can be useful to improve virtual screening performance of apo structures.
Collapse
Affiliation(s)
- Hugo Guterres
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Sang-Jun Park
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Wei Jiang
- Computational Science Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Wonpil Im
- Departments of Biological Sciences, Chemistry, Bioengineering, and Computer Science and Engineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| |
Collapse
|