1
|
Pettini F, Spiga O, Furini S, Fusi F. Electrophysiology, molecular modelling, and functional analysis of the effects of dietary quercetin and flavonoid analogues on K ir6.1 channels in rat stomach fundus smooth muscle. Biochem Pharmacol 2024; 220:115969. [PMID: 38086489 DOI: 10.1016/j.bcp.2023.115969] [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: 11/01/2023] [Revised: 12/01/2023] [Accepted: 12/04/2023] [Indexed: 12/17/2023]
Abstract
Flavonoids, ubiquitously distributed in the plant world, are regularly ingested with diets rich in fruit, vegetables, wine, and tea. During digestion, they are partially absorbed in the stomach. The present work aimed to assess the in vitro effects of quercetin and ten structurally related flavonoids on the rat gastric fundus smooth muscle, focussing on ATP-dependent K+ (Kir6.1) channels, which play a central role in the regulation of resting membrane potential, membrane excitability and, consequently, of gastric motility. Whole-cell currents through Kir6.1 channels (IKir6.1) were recorded with the patch-clamp technique and the mechanical activity of gastric fundus smooth muscle strips was studied under isometric conditions. Galangin ≈ tamarixetin > quercetin > kaempferol > isorhamnetin ≈ luteolin ≈ fisetin > (±)-taxifolin inhibited pinacidil-evoked, glibenclamide-sensitive IKir6.1 in a concentration-dependent manner. Morin, rutin, and myricetin were ineffective. The steric hindrance of the molecule and the number and position of hydroxyl groups on the B ring played an important role in the activity of the molecule. Molecular docking simulations revealed a possible binding site for flavonoids in the C-terminal domain of the Kir6.1 channel subunit SUR2B, in a flexible loop formed by residues 251 to 254 of chains C and D. Galangin and tamarixetin, but not rutin relaxed both high K+- and carbachol-induced contraction of fundus strips in a concentration-dependent manner. Furthermore, both flavonoids shifted to the right the concentration-relaxation curves to either pinacidil or L-cysteine constructed in strips pre-contracted by high K+, rutin being ineffective. In conclusion, IKir6.1 inhibition exerted by dietary flavonoids might counterbalance their myorelaxant activity, affect gastric accommodation or, at least, some stages of digestion.
Collapse
Affiliation(s)
- Francesco Pettini
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, Via A. Moro 2, 53100 Siena, Italy.
| | - Ottavia Spiga
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, Via A. Moro 2, 53100 Siena, Italy.
| | - Simone Furini
- Department of Electrical, Electronic and Information Engineering ″Guglielmo Marconi", University of Bologna, via dell'Università 50, 47521, Cesena, (FC), Italy.
| | - Fabio Fusi
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, Via A. Moro 2, 53100 Siena, Italy.
| |
Collapse
|
2
|
Ugurlu SY, McDonald D, Lei H, Jones AM, Li S, Tong HY, Butler MS, He S. Cobdock: an accurate and practical machine learning-based consensus blind docking method. J Cheminform 2024; 16:5. [PMID: 38212855 PMCID: PMC10785400 DOI: 10.1186/s13321-023-00793-x] [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/27/2023] [Accepted: 12/10/2023] [Indexed: 01/13/2024] Open
Abstract
Probing the surface of proteins to predict the binding site and binding affinity for a given small molecule is a critical but challenging task in drug discovery. Blind docking addresses this issue by performing docking on binding regions randomly sampled from the entire protein surface. However, compared with local docking, blind docking is less accurate and reliable because the docking space is too largetly sampled. Cavity detection-guided blind docking methods improved the accuracy by using cavity detection (also known as binding site detection) tools to guide the docking procedure. However, it is worth noting that the performance of these methods heavily relies on the quality of the cavity detection tool. This constraint, namely the dependence on a single cavity detection tool, significantly impacts the overall performance of cavity detection-guided methods. To overcome this limitation, we proposed Consensus Blind Dock (CoBDock), a novel blind, parallel docking method that uses machine learning algorithms to integrate docking and cavity detection results to improve not only binding site identification but also pose prediction accuracy. Our experiments on several datasets, including PDBBind 2020, ADS, MTi, DUD-E, and CASF-2016, showed that CoBDock has better binding site and binding mode performance than other state-of-the-art cavity detector tools and blind docking methods.
Collapse
Affiliation(s)
- Sadettin Y Ugurlu
- School of Computer Science, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK
| | | | - Huangshu Lei
- YaoPharma Co. Ltd., 100 Xingguang Avenue, Renhe Town, Yubei District, Chongqing, 401121, People's Republic of China
| | - Alan M Jones
- School of Pharmacy, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK
| | - Shu Li
- Centre for Artificial Intelligence Driven Drug Discovery, Macao Polytechnic University, R. de Luís Gonzaga Gomes, Macao, 5HV2+CP8, China
| | - Henry Y Tong
- Centre for Artificial Intelligence Driven Drug Discovery, Macao Polytechnic University, R. de Luís Gonzaga Gomes, Macao, 5HV2+CP8, China
| | | | - Shan He
- School of Computer Science, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK.
- AIA Insights Ltd, Birmingham, UK.
| |
Collapse
|
3
|
Mueller NPF, Carloni P, Alfonso-Prieto M. Molecular determinants of acrylamide neurotoxicity through covalent docking. Front Pharmacol 2023; 14:1125871. [PMID: 36937867 PMCID: PMC10018202 DOI: 10.3389/fphar.2023.1125871] [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: 12/16/2022] [Accepted: 02/17/2023] [Indexed: 03/06/2023] Open
Abstract
Acrylamide (ACR) is formed during food processing by Maillard reaction between sugars and proteins at high temperatures. It is also used in many industries, from water waste treatment to manufacture of paper, fabrics, dyes and cosmetics. Unfortunately, cumulative exposure to acrylamide, either from diet or at the workplace, may result in neurotoxicity. Such adverse effects arise from covalent adducts formed between acrylamide and cysteine residues of several neuronal proteins via a Michael addition reaction. The molecular determinants of acrylamide reactivity and its impact on protein function are not completely understood. Here we have compiled a list of acrylamide protein targets reported so far in the literature in connection with neurotoxicity and performed a systematic covalent docking study. Our results indicate that acrylamide binding to cysteine is favored in the presence of nearby positively charged amino acids, such as lysines and arginines. For proteins with more than one reactive Cys, docking scores were able to discriminate between the primary ACR modification site and secondary sites modified only at high ACR concentrations. Therefore, docking scores emerge as a potential filter to predict Cys reactivity against acrylamide. Inspection of the ACR-protein complex structures provides insights into the putative functional consequences of ACR modification, especially for non-enzyme proteins. Based on our study, covalent docking is a promising computational tool to predict other potential protein targets mediating acrylamide neurotoxicity.
Collapse
Affiliation(s)
- Nicolas Pierre Friedrich Mueller
- Institute for Advanced Simulations IAS-5, Institute of Neuroscience and Medicine INM-9, Computational Biomedicine, Forschungszentrum Jülich, Jülich, Germany
- Faculty of Mathematics and Natural Sciences, Heinrich Heine University Düsseldorf, Düsseldorf, Germany
| | - Paolo Carloni
- Institute for Advanced Simulations IAS-5, Institute of Neuroscience and Medicine INM-9, Computational Biomedicine, Forschungszentrum Jülich, Jülich, Germany
- Department of Physics, RWTH Aachen University, Aachen, Germany
| | - Mercedes Alfonso-Prieto
- Institute for Advanced Simulations IAS-5, Institute of Neuroscience and Medicine INM-9, Computational Biomedicine, Forschungszentrum Jülich, Jülich, Germany
- Cécile and Oskar Vogt Institute for Brain Research, Medical Faculty, University Hospital Düsseldorf, Heinrich Heine University Düsseldorf, Düsseldorf, Germany
- *Correspondence: Mercedes Alfonso-Prieto,
| |
Collapse
|
4
|
Buchanan CJ, Gaunt B, Harrison PJ, Yang Y, Liu J, Khan A, Giltrap AM, Le Bas A, Ward PN, Gupta K, Dumoux M, Tan TK, Schimaski L, Daga S, Picchiotti N, Baldassarri M, Benetti E, Fallerini C, Fava F, Giliberti A, Koukos PI, Davy MJ, Lakshminarayanan A, Xue X, Papadakis G, Deimel LP, Casablancas-Antràs V, Claridge TDW, Bonvin AMJJ, Sattentau QJ, Furini S, Gori M, Huo J, Owens RJ, Schaffitzel C, Berger I, Renieri A, Naismith JH, Baldwin AJ, Davis BG. Pathogen-sugar interactions revealed by universal saturation transfer analysis. Science 2022; 377:eabm3125. [PMID: 35737812 DOI: 10.1126/science.abm3125] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Many pathogens exploit host cell-surface glycans. However, precise analyses of glycan ligands binding with heavily modified pathogen proteins can be confounded by overlapping sugar signals and/or compounded with known experimental constraints. Universal saturation transfer analysis (uSTA) builds on existing nuclear magnetic resonance spectroscopy to provide an automated workflow for quantitating protein-ligand interactions. uSTA reveals that early-pandemic, B-origin-lineage severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike trimer binds sialoside sugars in an "end-on" manner. uSTA-guided modeling and a high-resolution cryo-electron microscopy structure implicate the spike N-terminal domain (NTD) and confirm end-on binding. This finding rationalizes the effect of NTD mutations that abolish sugar binding in SARS-CoV-2 variants of concern. Together with genetic variance analyses in early pandemic patient cohorts, this binding implicates a sialylated polylactosamine motif found on tetraantennary N-linked glycoproteins deep in the human lung as potentially relevant to virulence and/or zoonosis.
Collapse
Affiliation(s)
- Charles J Buchanan
- Rosalind Franklin Institute, Harwell Science and Innovation Campus, Oxford OX11 0FA, UK.,Department of Chemistry, University of Oxford, Oxford OX1 3TA, UK.,Kavli Institute of Nanoscience Discovery, University of Oxford, Oxford OX1 3QU, UK
| | - Ben Gaunt
- Rosalind Franklin Institute, Harwell Science and Innovation Campus, Oxford OX11 0FA, UK
| | - Peter J Harrison
- Division of Structural Biology, Wellcome Centre for Human Genetics, University of Oxford, Headington, Oxford OX3 7BN, UK.,Diamond Light Source, Harwell Science and Innovation Campus, Oxfordshire, UK
| | - Yun Yang
- Rosalind Franklin Institute, Harwell Science and Innovation Campus, Oxford OX11 0FA, UK.,Division of Structural Biology, Wellcome Centre for Human Genetics, University of Oxford, Headington, Oxford OX3 7BN, UK
| | - Jiwei Liu
- Rosalind Franklin Institute, Harwell Science and Innovation Campus, Oxford OX11 0FA, UK
| | - Aziz Khan
- Rosalind Franklin Institute, Harwell Science and Innovation Campus, Oxford OX11 0FA, UK.,Department of Chemistry, University of Oxford, Oxford OX1 3TA, UK
| | - Andrew M Giltrap
- Rosalind Franklin Institute, Harwell Science and Innovation Campus, Oxford OX11 0FA, UK.,Department of Chemistry, University of Oxford, Oxford OX1 3TA, UK
| | - Audrey Le Bas
- Rosalind Franklin Institute, Harwell Science and Innovation Campus, Oxford OX11 0FA, UK.,Division of Structural Biology, Wellcome Centre for Human Genetics, University of Oxford, Headington, Oxford OX3 7BN, UK
| | - Philip N Ward
- Division of Structural Biology, Wellcome Centre for Human Genetics, University of Oxford, Headington, Oxford OX3 7BN, UK
| | - Kapil Gupta
- Max Planck Bristol Centre for Minimal Biology, University of Bristol, Bristol, UK
| | - Maud Dumoux
- Rosalind Franklin Institute, Harwell Science and Innovation Campus, Oxford OX11 0FA, UK
| | - Tiong Kit Tan
- MRC Human Immunology Unit, Weatherall Institute of Molecular Medicine, University of Oxford, John Radcliffe Hospital, Oxford, UK
| | - Lisa Schimaski
- MRC Human Immunology Unit, Weatherall Institute of Molecular Medicine, University of Oxford, John Radcliffe Hospital, Oxford, UK
| | - Sergio Daga
- Medical Genetics, University of Siena, Siena, Italy.,Med Biotech Hub and Competence Center, Department of Medical Biotechnologies, University of Siena, Siena, Italy
| | - Nicola Picchiotti
- Department of Information Engineering and Mathematics, University of Siena, Siena, Italy.,Department of Mathematics, University of Pavia, Pavia, Italy
| | - Margherita Baldassarri
- Medical Genetics, University of Siena, Siena, Italy.,Med Biotech Hub and Competence Center, Department of Medical Biotechnologies, University of Siena, Siena, Italy
| | - Elisa Benetti
- Med Biotech Hub and Competence Center, Department of Medical Biotechnologies, University of Siena, Siena, Italy
| | - Chiara Fallerini
- Medical Genetics, University of Siena, Siena, Italy.,Med Biotech Hub and Competence Center, Department of Medical Biotechnologies, University of Siena, Siena, Italy
| | - Francesca Fava
- Medical Genetics, University of Siena, Siena, Italy.,Med Biotech Hub and Competence Center, Department of Medical Biotechnologies, University of Siena, Siena, Italy.,Genetica Medica, Azienda Ospedaliera Universitaria Senese, Siena, Italy
| | - Annarita Giliberti
- Medical Genetics, University of Siena, Siena, Italy.,Med Biotech Hub and Competence Center, Department of Medical Biotechnologies, University of Siena, Siena, Italy
| | - Panagiotis I Koukos
- Bijvoet Centre for Biomolecular Research, Faculty of Science, Utrecht University, Utrecht, Netherlands
| | - Matthew J Davy
- Rosalind Franklin Institute, Harwell Science and Innovation Campus, Oxford OX11 0FA, UK
| | - Abirami Lakshminarayanan
- Rosalind Franklin Institute, Harwell Science and Innovation Campus, Oxford OX11 0FA, UK.,Department of Chemistry, University of Oxford, Oxford OX1 3TA, UK
| | - Xiaochao Xue
- Department of Chemistry, University of Oxford, Oxford OX1 3TA, UK.,Sir William Dunn School of Pathology, Oxford, UK
| | | | | | - Virgínia Casablancas-Antràs
- Department of Chemistry, University of Oxford, Oxford OX1 3TA, UK.,Kavli Institute of Nanoscience Discovery, University of Oxford, Oxford OX1 3QU, UK
| | | | - Alexandre M J J Bonvin
- Bijvoet Centre for Biomolecular Research, Faculty of Science, Utrecht University, Utrecht, Netherlands
| | | | - Simone Furini
- Med Biotech Hub and Competence Center, Department of Medical Biotechnologies, University of Siena, Siena, Italy
| | - Marco Gori
- Department of Information Engineering and Mathematics, University of Siena, Siena, Italy.,Maasai, I3S CNRS, Université Côte d'Azur, Nice, France
| | - Jiandong Huo
- Rosalind Franklin Institute, Harwell Science and Innovation Campus, Oxford OX11 0FA, UK.,Division of Structural Biology, Wellcome Centre for Human Genetics, University of Oxford, Headington, Oxford OX3 7BN, UK
| | - Raymond J Owens
- Rosalind Franklin Institute, Harwell Science and Innovation Campus, Oxford OX11 0FA, UK.,Division of Structural Biology, Wellcome Centre for Human Genetics, University of Oxford, Headington, Oxford OX3 7BN, UK
| | - Christiane Schaffitzel
- Bijvoet Centre for Biomolecular Research, Faculty of Science, Utrecht University, Utrecht, Netherlands
| | - Imre Berger
- Bijvoet Centre for Biomolecular Research, Faculty of Science, Utrecht University, Utrecht, Netherlands
| | - Alessandra Renieri
- Medical Genetics, University of Siena, Siena, Italy.,Med Biotech Hub and Competence Center, Department of Medical Biotechnologies, University of Siena, Siena, Italy.,Genetica Medica, Azienda Ospedaliera Universitaria Senese, Siena, Italy
| | | | - James H Naismith
- Rosalind Franklin Institute, Harwell Science and Innovation Campus, Oxford OX11 0FA, UK.,Division of Structural Biology, Wellcome Centre for Human Genetics, University of Oxford, Headington, Oxford OX3 7BN, UK
| | - Andrew J Baldwin
- Rosalind Franklin Institute, Harwell Science and Innovation Campus, Oxford OX11 0FA, UK.,Department of Chemistry, University of Oxford, Oxford OX1 3TA, UK.,Kavli Institute of Nanoscience Discovery, University of Oxford, Oxford OX1 3QU, UK
| | - Benjamin G Davis
- Rosalind Franklin Institute, Harwell Science and Innovation Campus, Oxford OX11 0FA, UK.,Department of Chemistry, University of Oxford, Oxford OX1 3TA, UK.,Department of Pharmacology, University of Oxford, Oxford OX1 3QT, UK
| |
Collapse
|
5
|
Sergeeva OA, Mazur K, Reiner-Link D, Lutsenko K, Haas HL, Alfonso-Prieto M, Stark H. OLHA (N α-oleoylhistamine) modulates activity of mouse brain histaminergic neurons. Neuropharmacology 2022; 215:109167. [PMID: 35750238 DOI: 10.1016/j.neuropharm.2022.109167] [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: 01/28/2022] [Revised: 06/01/2022] [Accepted: 06/13/2022] [Indexed: 11/26/2022]
Abstract
Histaminergic (HA) neurons are located in the tuberomamillary nucleus (TMN) of the posterior hypothalamus, from where they project throughout the whole brain to control wakefulness. We examined the effects of Nα-oleoylhistamine (OLHA), a non-enzymatic condensation product of oleic acid (OLA) and histamine, on activity of mouse HA neurons in brain slices. OLHA bidirectionally modulated the firing of HA neurons. At 10 nM OLHA inhibited or had no action, whereas at 1 μM it evoked excitatory and inhibitory responses. Inhibition was not seen in presence of the histamine receptor H3 (H3R) antagonist clobenpropit and in calcium-free medium. Pre-incubation with a histamine-reuptake blocker prevented the decrease in firing by OLHA. OLHA-evoked increase in firing (EC50 ∼44 nM) was insensitive to blockers of cannabinoid 1 and 2 receptors and of the capsaicin receptor, but was significantly impaired by the peroxisome proliferator-activated receptor-alpha (PPAR-alpha) antagonist MK886, which suppressed also the rise in intracellular calcium level caused by OLHA. The OLHA-evoked excitation was mimicked by synthetic PPAR-alpha agonists (gemfibrozil and GW7647) and was abolished by the PKA inhibitor H-89. The H3R affinity (Ki) for histamine, measured in HEK293 cells with stable expression of human H3R, was higher than for OLHA (Ki: 42 vs 310 nM, respectively). Expression of PPAR-alpha was not different between TMN regions of males and females, responses to OLHA did not differ. Molecular modelling of PPAR-alpha bound to either OLHA or OEA showed similar binding energies. These findings shed light on a novel biotransformation product of histamine which may play a role in health and disease.
Collapse
Affiliation(s)
- Olga A Sergeeva
- Institute of Neural and Sensory Physiology, Medical Faculty, Heinrich Heine University Düsseldorf, 40225, Düsseldorf, Germany; Institute of Clinical Neurosciences and Medical Psychology, Medical Faculty, Heinrich Heine University Düsseldorf, 40225, Düsseldorf, Germany.
| | - Karolina Mazur
- Institute of Clinical Neurosciences and Medical Psychology, Medical Faculty, Heinrich Heine University Düsseldorf, 40225, Düsseldorf, Germany
| | - David Reiner-Link
- Institute of Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225, Düsseldorf, Germany
| | - Kiril Lutsenko
- Institute of Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225, Düsseldorf, Germany
| | - Helmut L Haas
- Institute of Neural and Sensory Physiology, Medical Faculty, Heinrich Heine University Düsseldorf, 40225, Düsseldorf, Germany
| | - Mercedes Alfonso-Prieto
- Cécile and Oskar Vogt Institute for Brain Research, Medical Faculty, Heinrich Heine University Düsseldorf, 40225, Düsseldorf, Germany; Computational Biomedicine, Institute for Advanced Simulation IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, 52425, Jülich, Germany
| | - Holger Stark
- Institute of Pharmaceutical and Medicinal Chemistry, Heinrich Heine University Düsseldorf, 40225, Düsseldorf, Germany
| |
Collapse
|
6
|
Koukos PI, Réau M, Bonvin AMJJ. Shape-Restrained Modeling of Protein-Small-Molecule Complexes with High Ambiguity Driven DOCKing. J Chem Inf Model 2021; 61:4807-4818. [PMID: 34436890 PMCID: PMC8479858 DOI: 10.1021/acs.jcim.1c00796] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Small-molecule docking remains one of the most valuable computational techniques for the structure prediction of protein-small-molecule complexes. It allows us to study the interactions between compounds and the protein receptors they target at atomic detail in a timely and efficient manner. Here, we present a new protocol in HADDOCK (High Ambiguity Driven DOCKing), our integrative modeling platform, which incorporates homology information for both receptor and compounds. It makes use of HADDOCK's unique ability to integrate information in the simulation to drive it toward conformations, which agree with the provided data. The focal point is the use of shape restraints derived from homologous compounds bound to the target receptors. We have developed two protocols: in the first, the shape is composed of dummy atom beads based on the position of the heavy atoms of the homologous template compound, whereas in the second, the shape is additionally annotated with pharmacophore data for some or all beads. For both protocols, ambiguous distance restraints are subsequently defined between those beads and the heavy atoms of the ligand to be docked. We have benchmarked the performance of these protocols with a fully unbound version of the widely used DUD-E (Database of Useful Decoys-Enhanced) dataset. In this unbound docking scenario, our template/shape-based docking protocol reaches an overall success rate of 81% when a reliable template can be identified (which was the case for 99 out of 102 complexes in the DUD-E dataset), which is close to the best results reported for bound docking on the DUD-E dataset.
Collapse
Affiliation(s)
- Panagiotis I Koukos
- Computational Structural Biology Group, Department of Chemistry, Faculty of Science, Utrecht University, Utrecht 3584CH, The Netherlands
| | - Manon Réau
- Computational Structural Biology Group, Department of Chemistry, Faculty of Science, Utrecht University, Utrecht 3584CH, The Netherlands
| | - Alexandre M J J Bonvin
- Computational Structural Biology Group, Department of Chemistry, Faculty of Science, Utrecht University, Utrecht 3584CH, The Netherlands
| |
Collapse
|
7
|
Erguven M, Karakulak T, Diril MK, Karaca E. How Far Are We from the Rapid Prediction of Drug Resistance Arising Due to Kinase Mutations? ACS OMEGA 2021; 6:1254-1265. [PMID: 33490784 PMCID: PMC7818309 DOI: 10.1021/acsomega.0c04672] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/23/2020] [Accepted: 12/11/2020] [Indexed: 06/12/2023]
Abstract
In all living organisms, protein kinases regulate various cell signaling events through phosphorylation. The phosphorylation occurs upon transferring an ATP's terminal phosphate to a target residue. Because of the central role of protein kinases in several proliferative pathways, point mutations occurring within the kinase's ATP-binding site can lead to a constitutively active enzyme, and ultimately, to cancer. A select set of these point mutations can also make the enzyme drug resistant toward the available kinase inhibitors. Because of technical and economical limitations, rapid experimental exploration of the impact of these mutations remains to be a challenge. This underscores the importance of kinase-ligand binding affinity prediction tools that are poised to measure the efficacy of inhibitors in the presence of kinase mutations. To this end, here, we compare the performances of six web-based scoring tools (DSX-ONLINE, KDEEP, HADDOCK2.2, PDBePISA, Pose&Rank, and PRODIGY-LIG) in assessing the impact of kinase mutations on their interactions with their inhibitors. This assessment is carried out on a new structure-based BINDKIN benchmark we compiled. BINDKIN contains wild-type and mutant structure pairs of kinase-inhibitor complexes, together with their corresponding experimental binding affinities (in the form of IC50, K d, and K i). The performance of various web servers over BINDKIN shows that they cannot predict the binding affinities (ΔGs) of wild-type and mutant cases directly. Still, they could catch whether a mutation improves or worsens the ligand binding (ΔΔGs) where the highest Pearson's R correlation coefficient is reached by DSX-ONLINE over the K i dataset. When homology models are used instead of K i-associated crystal structures, DSX-ONLINE loses its predictive capacity. These results highlight that there is room to improve the available scoring functions to estimate the impact of protein kinase point mutations on inhibitor binding. The BINDKIN benchmark with all related results is freely accessible online (https://github.com/CSB-KaracaLab/BINDKIN).
Collapse
Affiliation(s)
- Mehmet Erguven
- Izmir
Biomedicine and Genome Center, 35330 Izmir, Turkey
- Izmir
International Biomedicine and Genome Institute, Dokuz Eylul University, 35340 Izmir, Turkey
| | - Tülay Karakulak
- Izmir
Biomedicine and Genome Center, 35330 Izmir, Turkey
- Izmir
International Biomedicine and Genome Institute, Dokuz Eylul University, 35340 Izmir, Turkey
| | - M. Kasim Diril
- Izmir
Biomedicine and Genome Center, 35330 Izmir, Turkey
- Izmir
International Biomedicine and Genome Institute, Dokuz Eylul University, 35340 Izmir, Turkey
| | - Ezgi Karaca
- Izmir
Biomedicine and Genome Center, 35330 Izmir, Turkey
- Izmir
International Biomedicine and Genome Institute, Dokuz Eylul University, 35340 Izmir, Turkey
| |
Collapse
|
8
|
Gill SC, Mobley DL. Reversibly Sampling Conformations and Binding Modes Using Molecular Darting. J Chem Theory Comput 2021; 17:302-314. [PMID: 33289558 PMCID: PMC8121195 DOI: 10.1021/acs.jctc.0c00752] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Sampling multiple binding modes of a ligand in a single molecular dynamics simulation is difficult. A given ligand may have many internal degrees of freedom, along with many different ways it might orient itself in a binding site or across several binding sites, all of which might be separated by large energy barriers. We have developed a novel Monte Carlo move called molecular darting (MolDarting) to reversibly sample between predefined binding modes of a ligand. Here, we couple this with nonequilibrium candidate Monte Carlo (NCMC) to improve acceptance of moves. We apply this technique to a simple dipeptide system, a ligand binding to T4 lysozyme L99A, and ligand binding to HIV integrase to test this new method. We observe significant increases in acceptance compared to uniformly sampling the internal and rotational/translational degrees of freedom in these systems.
Collapse
Affiliation(s)
- Samuel C Gill
- Department of Chemistry, University of California, Irvine, California 92617, United States
| | - David L Mobley
- Department of Chemistry, University of California, Irvine, California 92617, United States
- Department of Pharmaceutical Sciences, University of California, Irvine, California 92617, United States
| |
Collapse
|
9
|
Schneider J, Korshunova K, Si Chaib Z, Giorgetti A, Alfonso-Prieto M, Carloni P. Ligand Pose Predictions for Human G Protein-Coupled Receptors: Insights from the Amber-Based Hybrid Molecular Mechanics/Coarse-Grained Approach. J Chem Inf Model 2020; 60:5103-5116. [PMID: 32786708 DOI: 10.1021/acs.jcim.0c00661] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Human G protein-coupled receptors (hGPCRs) are the most frequent targets of Food and Drug Administration (FDA)-approved drugs. Structural bioinformatics, along with molecular simulation, can support structure-based drug design targeting hGPCRs. In this context, several years ago, we developed a hybrid molecular mechanics (MM)/coarse-grained (CG) approach to predict ligand poses in low-resolution hGPCR models. The approach was based on the GROMOS96 43A1 and PRODRG united-atom force fields for the MM part. Here, we present a new MM/CG implementation using, instead, the Amber 14SB and GAFF all-atom potentials for proteins and ligands, respectively. The new implementation outperforms the previous one, as shown by a variety of applications on models of hGPCR/ligand complexes at different resolutions, and it is also more user-friendly. Thus, it emerges as a useful tool to predict poses in low-resolution models and provides insights into ligand binding similarly to all-atom molecular dynamics, albeit at a lower computational cost.
Collapse
Affiliation(s)
- Jakob Schneider
- Computational Biomedicine, Institute for Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany.,Department of Physics, RWTH Aachen University, 52074 Aachen, Germany.,JARA-Institute: Molecular Neuroscience and Neuroimaging, Institute for Neuroscience and Medicine INM-11/JARA-BRAIN Institute JBI-2, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany
| | - Ksenia Korshunova
- Computational Biomedicine, Institute for Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany.,Department of Physics, RWTH Aachen University, 52074 Aachen, Germany
| | - Zeineb Si Chaib
- Computational Biomedicine, Institute for Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany.,RWTH Aachen University, 52062 Aachen, Germany
| | - Alejandro Giorgetti
- Computational Biomedicine, Institute for Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany.,Department of Biotechnology, University of Verona, 37314 Verona, Italy.,JARA-HPC, IAS-5/INM-9 Computational Biomedicine, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany
| | - Mercedes Alfonso-Prieto
- Computational Biomedicine, Institute for Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany.,JARA-HPC, IAS-5/INM-9 Computational Biomedicine, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany.,Cecile and Oskar Vogt Institute for Brain Research, Medical Faculty, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| | - Paolo Carloni
- Computational Biomedicine, Institute for Advanced Simulations IAS-5/Institute for Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany.,Department of Physics, RWTH Aachen University, 52074 Aachen, Germany.,JARA-Institute: Molecular Neuroscience and Neuroimaging, Institute for Neuroscience and Medicine INM-11/JARA-BRAIN Institute JBI-2, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany.,JARA-HPC, IAS-5/INM-9 Computational Biomedicine, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany
| |
Collapse
|
10
|
Protein-ligand binding with the coarse-grained Martini model. Nat Commun 2020; 11:3714. [PMID: 32709852 PMCID: PMC7382508 DOI: 10.1038/s41467-020-17437-5] [Citation(s) in RCA: 96] [Impact Index Per Article: 24.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2019] [Accepted: 06/29/2020] [Indexed: 02/06/2023] Open
Abstract
The detailed understanding of the binding of small molecules to proteins is the key for the development of novel drugs or to increase the acceptance of substrates by enzymes. Nowadays, computer-aided design of protein–ligand binding is an important tool to accomplish this task. Current approaches typically rely on high-throughput docking essays or computationally expensive atomistic molecular dynamics simulations. Here, we present an approach to use the recently re-parametrized coarse-grained Martini model to perform unbiased millisecond sampling of protein–ligand interactions of small drug-like molecules. Remarkably, we achieve high accuracy without the need of any a priori knowledge of binding pockets or pathways. Our approach is applied to a range of systems from the well-characterized T4 lysozyme over members of the GPCR family and nuclear receptors to a variety of enzymes. The presented results open the way to high-throughput screening of ligand libraries or protein mutations using the coarse-grained Martini model. Computer-aided design of protein-ligand binding is important for the development of novel drugs. Here authors present an approach to use the recently re-parametrized coarse-grained Martini model to perform unbiased millisecond sampling of protein-ligand binding interactions of small drug-like molecules.
Collapse
|
11
|
Balhorn R, Balhorn MC, Balakrishnan K, Rebhun RB. The small molecule antibody mimic SH7139 targets a family of HLA-DRs expressed by B-cell lymphomas and other solid cancers. J Drug Target 2020; 28:1124-1136. [PMID: 32588667 DOI: 10.1080/1061186x.2020.1787418] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Selective high-affinity ligands (SHALs) belong to a novel class of small-molecule cancer therapeutics that function as targeted prodrugs. SH7139, the most advanced of the SHAL drugs designed to bind to a unique β-subunit structural epitope located on HLA-DR10, has exhibited exceptional preclinical efficacy and safety profiles. A comparison of SH7139 and SH7129, a biotin derivative of the drug developed for use as a diagnostic, showed the incorporation of a biotin tag did not alter the SHALs ability to target or kill HLA-DR10 expressing Raji cells. The use of SH7129 in an immuno-histochemical type assay to stain peripheral blood mononuclear cells (PBMCs) obtained from individuals expressing specific HLA-DRB1 alleles has also revealed that in addition to HLA-DR10, seven other more commonly expressed HLA-DRs are targeted by the drug. Computational dockings of the SHAL's recognition ligands to a number of HLA-DR structures explain, in part, why the targeting domains of SH7129 and SH7139 bind to some HLA-DRs but not others. The results also substantiate the selectivity of SH7129 and suggest it may prove useful as a companion diagnostic for pre-screening biopsy samples to identify those patients whose tumours should respond to SH7139 therapy.
Collapse
Affiliation(s)
| | | | - Karuppiah Balakrishnan
- Department of Immunology, School of Biological Sciences, Madurai Kamaraj University, Madurai, India
| | - Robert B Rebhun
- The Comparative Cancer Center, University of California, Davis, Davis, CA, USA
| |
Collapse
|
12
|
Don CG, Smieško M. In Silico Pharmacogenetics CYP2D6 Study Focused on the Pharmacovigilance of Herbal Antidepressants. Front Pharmacol 2020; 11:683. [PMID: 32477141 PMCID: PMC7237870 DOI: 10.3389/fphar.2020.00683] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2020] [Accepted: 04/27/2020] [Indexed: 01/27/2023] Open
Abstract
The annual increase in depression worldwide together with an upward trend in the use of alternative medicine as treatment asks for developing reliable safety profiles of herbal based medicine. A considerable risk on adverse reactions exists when herbal remedies are combined with prescription medication. Around 25% of the drugs, including many antidepressants, depend on the activity of CYP2D6 for their metabolism and corresponding efficacy. Therefore, probing CYP2D6 inhibition by the active substances in herbal based medicine within the wild-type enzyme and clinically relevant allelic variants is crucial to avoid toxicity issues. In this in silico study several compounds with herbal origin suggested to have antidepressant activity were analyzed on their CYP2D6 wild-type and CYP2D6*53 inhibition potential using molecular docking. In addition, several pharmacokinetic properties were evaluated to assess their probability to cross the blood brain barrier and subsequently reach sufficient brain bioavailability for the modulation of central nervous system targets as well as characteristics which may hint toward potential safety issues.
Collapse
|
13
|
Parks CD, Gaieb Z, Chiu M, Yang H, Shao C, Walters WP, Jansen JM, McGaughey G, Lewis RA, Bembenek SD, Ameriks MK, Mirzadegan T, Burley SK, Amaro RE, Gilson MK. D3R grand challenge 4: blind prediction of protein-ligand poses, affinity rankings, and relative binding free energies. J Comput Aided Mol Des 2020; 34:99-119. [PMID: 31974851 PMCID: PMC7261493 DOI: 10.1007/s10822-020-00289-y] [Citation(s) in RCA: 60] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2019] [Accepted: 01/13/2020] [Indexed: 12/11/2022]
Abstract
The Drug Design Data Resource (D3R) aims to identify best practice methods for computer aided drug design through blinded ligand pose prediction and affinity challenges. Herein, we report on the results of Grand Challenge 4 (GC4). GC4 focused on proteins beta secretase 1 and Cathepsin S, and was run in an analogous manner to prior challenges. In Stage 1, participant ability to predict the pose and affinity of BACE1 ligands were assessed. Following the completion of Stage 1, all BACE1 co-crystal structures were released, and Stage 2 tested affinity rankings with co-crystal structures. We provide an analysis of the results and discuss insights into determined best practice methods.
Collapse
Affiliation(s)
- Conor D Parks
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Zied Gaieb
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Michael Chiu
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Huanwang Yang
- RCSB Protein Data Bank, Institute for Quantitative Biomedicine, Rutgers, The State University of New Jersey, New Brunswick, NJ, 08903, USA
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Chenghua Shao
- RCSB Protein Data Bank, Institute for Quantitative Biomedicine, Rutgers, The State University of New Jersey, New Brunswick, NJ, 08903, USA
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, CA, 92093, USA
| | | | - Johanna M Jansen
- Novartis Institutes for BioMedical Research, Emeryville, CA, 94608, USA
| | | | - Richard A Lewis
- Novartis Institutes for BioMedical Research, Novartis Pharma AG, 4002, Basel, Switzerland
| | | | | | | | - Stephen K Burley
- RCSB Protein Data Bank, Institute for Quantitative Biomedicine, Rutgers, The State University of New Jersey, New Brunswick, NJ, 08903, USA
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Rommie E Amaro
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA, 92093, USA.
- Department of Chemistry and Biochemistry, UC San Diego, La Jolla, CA, 92093-0340, USA.
| | - Michael K Gilson
- Drug Design Data Resource, University of California, San Diego, La Jolla, CA, 92093, USA.
- Skaggs School of Pharmacy and Pharmaceutical Sciences, UC San Diego, 9500 Gilman Drive, MC0751, La Jolla, CA, 92093, USA.
| |
Collapse
|
14
|
Kadukova M, Chupin V, Grudinin S. Docking rigid macrocycles using Convex-PL, AutoDock Vina, and RDKit in the D3R Grand Challenge 4. J Comput Aided Mol Des 2019; 34:191-200. [PMID: 31784861 DOI: 10.1007/s10822-019-00263-3] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2019] [Accepted: 11/22/2019] [Indexed: 12/15/2022]
Abstract
The D3R Grand Challenge 4 provided a brilliant opportunity to test macrocyclic docking protocols on a diverse high-quality experimental data. We participated in both pose and affinity prediction exercises. Overall, we aimed to use an automated structure-based docking pipeline built around a set of tools developed in our team. This exercise again demonstrated a crucial importance of the correct local ligand geometry for the overall success of docking. Starting from the second part of the pose prediction stage, we developed a stable pipeline for sampling macrocycle conformers. This resulted in the subangstrom average precision of our pose predictions. In the affinity prediction exercise we obtained average results. However, we could improve these when using docking poses submitted by the best predictors. Our docking tools including the Convex-PL scoring function are available at https://team.inria.fr/nano-d/software/.
Collapse
Affiliation(s)
- Maria Kadukova
- Univ. Grenoble Alpes, CNRS, Inria, Grenoble INP, LJK, 38000, Grenoble, France
- Moscow Institute of Physics and Technology, Dolgoprudny, Russia, 141700
| | - Vladimir Chupin
- Moscow Institute of Physics and Technology, Dolgoprudny, Russia, 141700
| | - Sergei Grudinin
- Univ. Grenoble Alpes, CNRS, Inria, Grenoble INP, LJK, 38000, Grenoble, France.
| |
Collapse
|
15
|
Exploring fragment-based target-specific ranking protocol with machine learning on cathepsin S. J Comput Aided Mol Des 2019; 33:1095-1105. [PMID: 31729618 DOI: 10.1007/s10822-019-00247-3] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2019] [Accepted: 11/02/2019] [Indexed: 12/12/2022]
Abstract
Cathepsin S (CatS), a member of cysteine cathepsin proteases, has been well studied due to its significant role in many pathological processes, including arthritis, cancer and cardiovascular diseases. CatS inhibitors have been included in D3R-GC3 for both docking pose prediction and affinity ranking, and in D3R-GC4 for binding affinity ranking. The difficulties posed by CatS inhibitors in D3R mainly come from three aspects: large size, high flexibility and similar chemical structures. We have participated in GC4; our best submitted model, which employs a similarity-based alignment docking and Vina scoring protocol, yielded Kendall's τ of 0.23 for 459 binders in GC4. In our further explorations with machine learning, by curating a CatS specific training set, adopting a similarity-based constrained docking method as well as an arm-based fragmentation strategy which can describe large inhibitors in a locality-sensitive fashion, our best structure-based ranking protocol can achieve Kendall's τ of 0.52 for all binders in GC4. In this exploration process, we have demonstrated the importance of training data, docking approaches and fragmentation strategies in inhibitor-ranking protocol development with machine learning.
Collapse
|
16
|
Coupling enhanced sampling of the apo-receptor with template-based ligand conformers selection: performance in pose prediction in the D3R Grand Challenge 4. J Comput Aided Mol Des 2019; 34:149-162. [DOI: 10.1007/s10822-019-00244-6] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2019] [Accepted: 10/30/2019] [Indexed: 12/20/2022]
|
17
|
Basciu A, Malloci G, Pietrucci F, Bonvin AMJJ, Vargiu AV. Holo-like and Druggable Protein Conformations from Enhanced Sampling of Binding Pocket Volume and Shape. J Chem Inf Model 2019; 59:1515-1528. [PMID: 30883122 DOI: 10.1021/acs.jcim.8b00730] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Understanding molecular recognition of small molecules by proteins in atomistic detail is key for drug design. Molecular docking is a widely used computational method to mimic ligand-protein association in silico. However, predicting conformational changes occurring in proteins upon ligand binding is still a major challenge. Ensemble docking approaches address this issue by considering a set of different conformations of the protein obtained either experimentally or from computer simulations, e.g., molecular dynamics. However, holo structures prone to host (the correct) ligands are generally poorly sampled by standard molecular dynamics simulations of the apo protein. In order to address this limitation, we introduce a computational approach based on metadynamics simulations called ensemble docking with enhanced sampling of pocket shape (EDES) that allows holo-like conformations of proteins to be generated by exploiting only their apo structures. This is achieved by defining a set of collective variables that effectively sample different shapes of the binding site, ultimately mimicking the steric effect due to the ligand. We assessed the method on three challenging proteins undergoing different extents of conformational changes upon ligand binding. In all cases our protocol generates a significant fraction of structures featuring a low RMSD from the experimental holo geometry. Moreover, ensemble docking calculations using those conformations yielded in all cases native-like poses among the top-ranked ones.
Collapse
Affiliation(s)
- Andrea Basciu
- Dipartimento di Fisica , Università di Cagliari, Cittadella Universitaria , I- 09042 Monserrato (CA) , Italy
| | - Giuliano Malloci
- Dipartimento di Fisica , Università di Cagliari, Cittadella Universitaria , I- 09042 Monserrato (CA) , Italy
| | - Fabio Pietrucci
- Sorbonne Université , Muséum National d'Histoire Naturelle, UMR CNRS 7590, IRD, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, IMPMC , F-75005 Paris , France
| | - Alexandre M J J Bonvin
- Bijvoet Center for Biomolecular Research, Faculty of Science - Chemistry , Utrecht University , Padualaan 8 , 3584 CH Utrecht , The Netherlands
| | - Attilio V Vargiu
- Dipartimento di Fisica , Università di Cagliari, Cittadella Universitaria , I- 09042 Monserrato (CA) , Italy.,Bijvoet Center for Biomolecular Research, Faculty of Science - Chemistry , Utrecht University , Padualaan 8 , 3584 CH Utrecht , The Netherlands
| |
Collapse
|