1
|
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
|
2
|
Wack JS, Brahm K, Babel P, Dalton JAR, Schmitz K. Effect of macrocyclization and tetramethylrhodamine labeling on chemokine binding peptides. J Pept Sci 2023:e3486. [PMID: 36843216 DOI: 10.1002/psc.3486] [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: 12/08/2022] [Revised: 02/05/2023] [Accepted: 02/21/2023] [Indexed: 02/28/2023]
Abstract
Receptor-derived peptides have played an important role in elucidating chemokine-receptor interactions. For the inflammatory chemokine CXC-class chemokine ligand 8 (CXCL8), a site II-mimetic peptide has been derived from parts of extracellular loops 2 and 3 and adjacent transmembrane helices of its receptor CXC-class chemokine receptor 1 (Helmer et al., RSC Adv., 2015, 5, 25657). The peptide sequence with a C-terminal glutamine did not bind to CXCL8, whereas one with a C-terminal glutamate did but with low micromolar affinity. We sought to improve the affinity and protease stability of the latter peptide through cyclization while also cyclizing the former for control purposes. To identify a cyclization strategy that permits a receptor-like interaction, we conducted a molecular dynamics simulation of CXCL8 in complex with full-length CXC-class chemokine receptor 1. We introduced a linker to provide an appropriate spacing between the termini and used an on-resin side-chain-to-tail cyclization strategy. Upon chemokine binding, the fluorescence intensity of the tetramethylrhodamine (TAMRA)-labeled cyclic peptides increased whereas the fluorescence anisotropy decreased. Additional molecular dynamics simulations indicated that the fluorophore interacts with the peptide macrocycle so that chemokine binding leads to its displacement and observed changes in fluorescence. Macrocyclization of both 18-amino acid-long peptides led to the same low micromolar affinity for CXCL8. Likewise, both TAMRA-labeled linear peptides interacted with CXCL8 with similar affinities. Interestingly, the linear TAMRA-labeled peptides were more resistant to tryptic digestion than the unlabeled counterparts, whereas the cyclized peptides were not degraded at all. We conclude that the TAMRA fluorophore tends to interact with peptides altering their protease stability and behavior in fluorescence-based assays.
Collapse
Affiliation(s)
- Julia S Wack
- Biological Chemistry, Clemens-Schöpf-Institute for Organic Chemistry and Biochemistry, Technical University of Darmstadt, Darmstadt, Germany
| | - Kevin Brahm
- Biological Chemistry, Clemens-Schöpf-Institute for Organic Chemistry and Biochemistry, Technical University of Darmstadt, Darmstadt, Germany
| | - Philipp Babel
- Computational Biology and Simulation, Technical University of Darmstadt, Darmstadt, Germany
| | - James A R Dalton
- Laboratory of Molecular Neuropharmacology and Bioinformatics, Unitat de Bioestadística, Institut de Neurociències, Universitat Autònoma de Barcelona, Bellaterra, Spain.,Ronin Institute, Montclair, New Jersey, USA
| | - Katja Schmitz
- Biological Chemistry, Clemens-Schöpf-Institute for Organic Chemistry and Biochemistry, Technical University of Darmstadt, Darmstadt, Germany
| |
Collapse
|
3
|
Thakur N, Ray AP, Sharp L, Jin B, Duong A, Pour NG, Obeng S, Wijesekara AV, Gao ZG, McCurdy CR, Jacobson KA, Lyman E, Eddy MT. Anionic phospholipids control mechanisms of GPCR-G protein recognition. Nat Commun 2023; 14:794. [PMID: 36781870 PMCID: PMC9925817 DOI: 10.1038/s41467-023-36425-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2022] [Accepted: 01/31/2023] [Indexed: 02/15/2023] Open
Abstract
G protein-coupled receptors (GPCRs) are embedded in phospholipids that strongly influence drug-stimulated signaling. Anionic lipids are particularly important for GPCR signaling complex formation, but a mechanism for this role is not understood. Using NMR spectroscopy, we explore the impact of anionic lipids on the function-related conformational equilibria of the human A2A adenosine receptor (A2AAR) in bilayers containing defined mixtures of zwitterionic and anionic phospholipids. Anionic lipids prime the receptor to form complexes with G proteins through a conformational selection process. Without anionic lipids, signaling complex formation proceeds through a less favorable induced fit mechanism. In computational models, anionic lipids mimic interactions between a G protein and positively charged residues in A2AAR at the receptor intracellular surface, stabilizing a pre-activated receptor conformation. Replacing these residues strikingly alters the receptor response to anionic lipids in experiments. High sequence conservation of the same residues among all GPCRs supports a general role for lipid-receptor charge complementarity in signaling.
Collapse
Affiliation(s)
- Naveen Thakur
- Department of Chemistry, College of Liberal Arts & Sciences, University of Florida, 126 Sisler Hall, Gainesville, FL, USA
| | - Arka P Ray
- Department of Chemistry, College of Liberal Arts & Sciences, University of Florida, 126 Sisler Hall, Gainesville, FL, USA
| | - Liam Sharp
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware, USA
| | - Beining Jin
- Department of Chemistry, College of Liberal Arts & Sciences, University of Florida, 126 Sisler Hall, Gainesville, FL, USA
| | - Alexander Duong
- Department of Chemistry, College of Liberal Arts & Sciences, University of Florida, 126 Sisler Hall, Gainesville, FL, USA
| | - Niloofar Gopal Pour
- Department of Chemistry, College of Liberal Arts & Sciences, University of Florida, 126 Sisler Hall, Gainesville, FL, USA
| | - Samuel Obeng
- Department of Pharmacodynamics, College of Pharmacy, University of Florida, Gainesville, FL, USA
- Department of Medicinal Chemistry, College of Pharmacy, University of Florida, Gainesville, FL, USA
| | - Anuradha V Wijesekara
- Department of Chemistry, College of Liberal Arts & Sciences, University of Florida, 126 Sisler Hall, Gainesville, FL, USA
| | - Zhan-Guo Gao
- Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland, 20892, USA
| | - Christopher R McCurdy
- Department of Medicinal Chemistry, College of Pharmacy, University of Florida, Gainesville, FL, USA
- Translational Drug Development Core, Clinical and Translational Sciences Institute, University of Florida, Gainesville, Florida, USA
| | - Kenneth A Jacobson
- Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland, 20892, USA
| | - Edward Lyman
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware, USA
- Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware, USA
| | - Matthew T Eddy
- Department of Chemistry, College of Liberal Arts & Sciences, University of Florida, 126 Sisler Hall, Gainesville, FL, USA.
| |
Collapse
|
4
|
Thakur N, Ray AP, Sharp L, Jin B, Duong A, Pour NG, Obeng S, Wijesekara AV, Gao ZG, McCurdy CR, Jacobson KA, Lyman E, Eddy MT. Anionic Phospholipids Control Mechanisms of GPCR-G Protein Recognition. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.01.11.523010. [PMID: 36711594 PMCID: PMC9882065 DOI: 10.1101/2023.01.11.523010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
G protein-coupled receptors (GPCRs) are embedded in phospholipids that strongly influence drug-stimulated signaling. Anionic lipids are particularly important for GPCR signaling complex formation, but a mechanism for this role is not understood. Using NMR spectroscopy, we visualized the impact of anionic lipids on the function-related conformational equilibria of the human A 2A adenosine receptor (A 2A AR) in bilayers containing defined mixtures of zwitterionic and anionic phospholipids. Anionic lipids primed the receptor to form complexes with G proteins through a conformational selection process. Without anionic lipids, signaling complex formation proceeded through a less favorable induced fit mechanism. In computational models, anionic lipids mimicked interactions between a G protein and positively charged residues in A 2A AR at the receptor intracellular surface, stabilizing a pre-activated receptor conformation. Replacing these residues strikingly altered the receptor response to anionic lipids in experiments. High sequence conservation of the same residues among all GPCRs supports a general role for lipid-receptor charge complementarity in signaling.
Collapse
Affiliation(s)
- Naveen Thakur
- Department of Chemistry, College of Liberal Arts & Sciences, University of Florida, 126 Sisler Hall, Gainesville, FL 32611, USA
| | - Arka P Ray
- Department of Chemistry, College of Liberal Arts & Sciences, University of Florida, 126 Sisler Hall, Gainesville, FL 32611, USA
| | - Liam Sharp
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware
| | - Beining Jin
- Department of Chemistry, College of Liberal Arts & Sciences, University of Florida, 126 Sisler Hall, Gainesville, FL 32611, USA
| | - Alexander Duong
- Department of Chemistry, College of Liberal Arts & Sciences, University of Florida, 126 Sisler Hall, Gainesville, FL 32611, USA
| | - Niloofar Gopal Pour
- Department of Chemistry, College of Liberal Arts & Sciences, University of Florida, 126 Sisler Hall, Gainesville, FL 32611, USA
| | - Samuel Obeng
- Department of Pharmacodynamics, College of Pharmacy, University of Florida, Gainesville, FL 32610, USA
- Department of Medicinal Chemistry, College of Pharmacy, University of Florida, Gainesville, FL 32610, USA
| | - Anuradha V Wijesekara
- Department of Chemistry, College of Liberal Arts & Sciences, University of Florida, 126 Sisler Hall, Gainesville, FL 32611, USA
| | - Zhan-Guo Gao
- Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Christopher R McCurdy
- Department of Medicinal Chemistry, College of Pharmacy, University of Florida, Gainesville, FL 32610, USA
- Translational Drug Development Core, Clinical and Translational Sciences Institute, University of Florida, Gainesville, Florida 32610, USA
| | - Kenneth A Jacobson
- Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland 20892, United States
| | - Edward Lyman
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware
- Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware
| | - Matthew T Eddy
- Department of Chemistry, College of Liberal Arts & Sciences, University of Florida, 126 Sisler Hall, Gainesville, FL 32611, USA
| |
Collapse
|
5
|
Kaczor AA, Wróbel TM, Bartuzi D. Allosteric Modulators of Dopamine D 2 Receptors for Fine-Tuning of Dopaminergic Neurotransmission in CNS Diseases: Overview, Pharmacology, Structural Aspects and Synthesis. Molecules 2022; 28:molecules28010178. [PMID: 36615372 PMCID: PMC9822192 DOI: 10.3390/molecules28010178] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2022] [Revised: 12/19/2022] [Accepted: 12/21/2022] [Indexed: 12/28/2022] Open
Abstract
Allosteric modulation of G protein-coupled receptors (GPCRs) is nowadays a hot topic in medicinal chemistry. Allosteric modulators, i.e., compounds which bind in a receptor site topologically distinct from orthosteric sites, exhibit a number of advantages. They are more selective, safer and display a ceiling effect which prevents overdosing. Allosteric modulators of dopamine D2 receptor are potential drugs against a number of psychiatric and neurological diseases, such as schizophrenia and Parkinson's disease. In this review, an insightful summary of current research on D2 receptor modulators is presented, ranging from their pharmacology and structural aspects of ligand-receptor interactions to their synthesis.
Collapse
Affiliation(s)
- Agnieszka A. Kaczor
- Department of Synthesis and Chemical Technology of Pharmaceutical Substances with Computer Modeling Laboratory, Faculty of Pharmacy, Medical University of Lublin, 4A Chodźki St., PL-20093 Lublin, Poland
- School of Pharmacy, University of Eastern Finland, Yliopistonranta 1, P.O. Box 1627, FI-70211 Kuopio, Finland
- Correspondence: ; Tel.: +48-81-448-72-73
| | - Tomasz M. Wróbel
- Department of Synthesis and Chemical Technology of Pharmaceutical Substances with Computer Modeling Laboratory, Faculty of Pharmacy, Medical University of Lublin, 4A Chodźki St., PL-20093 Lublin, Poland
| | - Damian Bartuzi
- Department of Synthesis and Chemical Technology of Pharmaceutical Substances with Computer Modeling Laboratory, Faculty of Pharmacy, Medical University of Lublin, 4A Chodźki St., PL-20093 Lublin, Poland
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, SE-75124 Uppsala, Sweden
| |
Collapse
|
6
|
Díaz O, Renault P, Giraldo J. Evaluating Allosteric Perturbations in Cannabinoid Receptor 1 by In Silico Single-Point Mutation. ACS OMEGA 2022; 7:37873-37884. [PMID: 36312415 PMCID: PMC9608382 DOI: 10.1021/acsomega.2c04980] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/04/2022] [Accepted: 09/28/2022] [Indexed: 06/16/2023]
Abstract
Cannabinoid receptor 1 (CB1) is a promising drug target involved in many physiological processes. Using atomistic molecular dynamics (MD) simulations, we examined the structural effect of F237L mutation on CB1, a mutation that has qualitatively similar effects to allosteric ligand ORG27569 binding. This mutation showed a global effect on CB1 conformations. Among the observed effects, TM6 outward movement and the conformational change of the NPxxY motif upon receptor activation by CB1 agonist CP55940 were hindered compared to wt CB1. Within the orthosteric binding site, CP55940 interactions with CB1 were altered. Our results revealed that allosteric perturbations introduced by the mutation had a global impact on receptor conformations, suggesting that the mutation site is a key region for allosteric modulation in CB1.
Collapse
Affiliation(s)
- Oscar Díaz
- Laboratory
of Molecular Neuropharmacology and Bioinformatics, Unitat de Bioestadística
and Institut de Neurociències, Universitat
Autònoma de Barcelona, Bellaterra 08193, Spain
- Instituto
de Salud Carlos III, Centro de Investigación
Biomédica en Red de Salud Mental (CIBERSAM), Madrid 28029, Spain
- Unitat
de Neurociència Traslacional, Parc Taulí Hospital Universitari,
Institut d’Investigació i Innovació Parc Taulí
(I3PT), Institut de Neurociències, Universitat Autònoma de Barcelona, Bellaterra 08193, Spain
| | - Pedro Renault
- Laboratory
of Molecular Neuropharmacology and Bioinformatics, Unitat de Bioestadística
and Institut de Neurociències, Universitat
Autònoma de Barcelona, Bellaterra 08193, Spain
- Instituto
de Salud Carlos III, Centro de Investigación
Biomédica en Red de Salud Mental (CIBERSAM), Madrid 28029, Spain
- Unitat
de Neurociència Traslacional, Parc Taulí Hospital Universitari,
Institut d’Investigació i Innovació Parc Taulí
(I3PT), Institut de Neurociències, Universitat Autònoma de Barcelona, Bellaterra 08193, Spain
| | - Jesús Giraldo
- Laboratory
of Molecular Neuropharmacology and Bioinformatics, Unitat de Bioestadística
and Institut de Neurociències, Universitat
Autònoma de Barcelona, Bellaterra 08193, Spain
- Instituto
de Salud Carlos III, Centro de Investigación
Biomédica en Red de Salud Mental (CIBERSAM), Madrid 28029, Spain
- Unitat
de Neurociència Traslacional, Parc Taulí Hospital Universitari,
Institut d’Investigació i Innovació Parc Taulí
(I3PT), Institut de Neurociències, Universitat Autònoma de Barcelona, Bellaterra 08193, Spain
| |
Collapse
|
7
|
Ding T, Karlov DS, Pino-Angeles A, Tikhonova IG. Intermolecular Interactions in G Protein-Coupled Receptor Allosteric Sites at the Membrane Interface from Molecular Dynamics Simulations and Quantum Chemical Calculations. J Chem Inf Model 2022; 62:4736-4747. [PMID: 36178787 PMCID: PMC9554917 DOI: 10.1021/acs.jcim.2c00788] [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] [Indexed: 11/30/2022]
Abstract
Allosteric modulators are called promising candidates in G protein-coupled receptor (GPCR) drug development by displaying subtype selectivity and more specific receptor modulation. Among the allosteric sites known to date, cavities at the receptor-lipid interface represent an uncharacteristic binding location that raises many questions about the ligand interactions and stability, the binding site structure, and how all of these are affected by lipid molecules. In this work, we analyze interactions in the allosteric sites of the PAR2, C5aR1, and GCGR receptors in three lipid compositions using molecular dynamics simulations. In addition, we performed quantum chemical calculations involving the symmetry-adapted perturbation theory (SAPT) and the natural population analysis to quantify the strength of intermolecular interactions. We show that besides classical hydrogen bonds, weak polar interactions such as O-HC, O-Br, and long-range electrostatics with the backbone amides contribute to the stability of allosteric modulators at the receptor-lipid interface. The allosteric cavities are detectable in various membrane compositions. The availability of polar atoms for interactions in such cavities can be assessed by water molecules from simulations. Although ligand-lipid interactions are weak, lipid tails play a role in ligand binding pose stability and the size of allosteric cavities. We discuss physicochemical aspects of ligand binding at the receptor-lipid interface and suggest a compound library enriched by weak donor groups for ligand search in such sites.
Collapse
Affiliation(s)
- Tianyi Ding
- School of Pharmacy, Medical Biology Centre, Queen's University Belfast, Belfast, Northern IrelandBT9 7BL, U.K
| | - Dmitry S Karlov
- School of Pharmacy, Medical Biology Centre, Queen's University Belfast, Belfast, Northern IrelandBT9 7BL, U.K
| | - Almudena Pino-Angeles
- School of Pharmacy, Medical Biology Centre, Queen's University Belfast, Belfast, Northern IrelandBT9 7BL, U.K
| | - Irina G Tikhonova
- School of Pharmacy, Medical Biology Centre, Queen's University Belfast, Belfast, Northern IrelandBT9 7BL, U.K
| |
Collapse
|
8
|
Bedart C, Renault N, Chavatte P, Porcherie A, Lachgar A, Capron M, Farce A. SINAPs: A Software Tool for Analysis and Visualization of Interaction Networks of Molecular Dynamics Simulations. J Chem Inf Model 2022; 62:1425-1436. [PMID: 35239339 PMCID: PMC8966674 DOI: 10.1021/acs.jcim.1c00854] [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] [Indexed: 11/29/2022]
Abstract
As long as the structural study of molecular mechanisms requires multiple molecular dynamics reflecting contrasted bioactive states, the subsequent analysis of molecular interaction networks remains a bottleneck to be fairly treated and requires a user-friendly 3D view of key interactions. Structural Interaction Network Analysis Protocols (SINAPs) is a proprietary python tool developed to (i) quickly solve key interactions able to distinguish two protein states, either from two sets of molecular dynamics simulations or from two crystallographic structures, and (ii) render a user-friendly 3D view of these key interactions through a plugin of UCSF Chimera, one of the most popular open-source viewing software for biomolecular systems. Through two case studies, glucose transporter-1 (GLUT-1) and A2A adenosine receptor (A2AR), SINAPs easily pinpointed key interactions observed experimentally and relevant for their bioactivities. This very effective tool was thus applied to identify the amino acids involved in the molecular enzymatic mechanisms ruling the activation of an immunomodulator drug candidate, P28 glutathione-S-transferase (P28GST). SINAPs is freely available at https://github.com/ParImmune/SINAPs.
Collapse
Affiliation(s)
- Corentin Bedart
- Univ.
Lille, Inserm, CHU Lille, U1286 - Infinite - Institute for Translational
Research in Inflammation, F-59000 Lille, France,Par’Immune,
Bio-incubateur Eurasanté, 70 rue du Dr. Yersin, 59120 Loos-Lez-Lille, France,
| | - Nicolas Renault
- Univ.
Lille, Inserm, CHU Lille, U1286 - Infinite - Institute for Translational
Research in Inflammation, F-59000 Lille, France
| | - Philippe Chavatte
- Univ.
Lille, Inserm, CHU Lille, U1286 - Infinite - Institute for Translational
Research in Inflammation, F-59000 Lille, France
| | - Adeline Porcherie
- Par’Immune,
Bio-incubateur Eurasanté, 70 rue du Dr. Yersin, 59120 Loos-Lez-Lille, France
| | - Abderrahim Lachgar
- Par’Immune,
Bio-incubateur Eurasanté, 70 rue du Dr. Yersin, 59120 Loos-Lez-Lille, France
| | - Monique Capron
- Univ.
Lille, Inserm, CHU Lille, U1286 - Infinite - Institute for Translational
Research in Inflammation, F-59000 Lille, France,Par’Immune,
Bio-incubateur Eurasanté, 70 rue du Dr. Yersin, 59120 Loos-Lez-Lille, France
| | - Amaury Farce
- Univ.
Lille, Inserm, CHU Lille, U1286 - Infinite - Institute for Translational
Research in Inflammation, F-59000 Lille, France,
| |
Collapse
|
9
|
Wang J, Bhattarai A, Do HN, Akhter S, Miao Y. Molecular Simulations and Drug Discovery of Adenosine Receptors. Molecules 2022; 27:2054. [PMID: 35408454 PMCID: PMC9000248 DOI: 10.3390/molecules27072054] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2022] [Revised: 03/18/2022] [Accepted: 03/20/2022] [Indexed: 02/02/2023] Open
Abstract
G protein-coupled receptors (GPCRs) represent the largest family of human membrane proteins. Four subtypes of adenosine receptors (ARs), the A1AR, A2AAR, A2BAR and A3AR, each with a unique pharmacological profile and distribution within the tissues in the human body, mediate many physiological functions and serve as critical drug targets for treating numerous human diseases including cancer, neuropathic pain, cardiac ischemia, stroke and diabetes. The A1AR and A3AR preferentially couple to the Gi/o proteins, while the A2AAR and A2BAR prefer coupling to the Gs proteins. Adenosine receptors were the first subclass of GPCRs that had experimental structures determined in complex with distinct G proteins. Here, we will review recent studies in molecular simulations and computer-aided drug discovery of the adenosine receptors and also highlight their future research opportunities.
Collapse
Affiliation(s)
| | | | | | | | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS 66047, USA; (J.W.); (A.B.); (H.N.D.); (S.A.)
| |
Collapse
|
10
|
Pluhackova K, Wilhelm FM, Müller DJ. Lipids and Phosphorylation Conjointly Modulate Complex Formation of β 2-Adrenergic Receptor and β-arrestin2. Front Cell Dev Biol 2022; 9:807913. [PMID: 35004696 PMCID: PMC8733679 DOI: 10.3389/fcell.2021.807913] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2021] [Accepted: 11/30/2021] [Indexed: 01/14/2023] Open
Abstract
G protein-coupled receptors (GPCRs) are the largest class of human membrane proteins that bind extracellular ligands at their orthosteric binding pocket to transmit signals to the cell interior. Ligand binding evokes conformational changes in GPCRs that trigger the binding of intracellular interaction partners (G proteins, G protein kinases, and arrestins), which initiate diverse cellular responses. It has become increasingly evident that the preference of a GPCR for a certain intracellular interaction partner is modulated by a diverse range of factors, e.g., ligands or lipids embedding the transmembrane receptor. Here, by means of molecular dynamics simulations of the β2-adrenergic receptor and β-arrestin2, we study how membrane lipids and receptor phosphorylation regulate GPCR-arrestin complex conformation and dynamics. We find that phosphorylation drives the receptor’s intracellular loop 3 (ICL3) away from a native negatively charged membrane surface to interact with arrestin. If the receptor is embedded in a neutral membrane, the phosphorylated ICL3 attaches to the membrane surface, which widely opens the receptor core. This opening, which is similar to the opening in the G protein-bound state, weakens the binding of arrestin. The loss of binding specificity is manifested by shallower arrestin insertion into the receptor core and higher dynamics of the receptor-arrestin complex. Our results show that receptor phosphorylation and the local membrane composition cooperatively fine-tune GPCR-mediated signal transduction. Moreover, the results suggest that deeper understanding of complex GPCR regulation mechanisms is necessary to discover novel pathways of pharmacological intervention.
Collapse
Affiliation(s)
- Kristyna Pluhackova
- Department of Biosystems Science and Engineering, Eidgenössische Technische Hochschule (ETH) Zurich, Basel, Switzerland
| | - Florian M Wilhelm
- Department of Biosystems Science and Engineering, Eidgenössische Technische Hochschule (ETH) Zurich, Basel, Switzerland
| | - Daniel J Müller
- Department of Biosystems Science and Engineering, Eidgenössische Technische Hochschule (ETH) Zurich, Basel, Switzerland
| |
Collapse
|
11
|
Róg T, Girych M, Bunker A. Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design. Pharmaceuticals (Basel) 2021; 14:1062. [PMID: 34681286 PMCID: PMC8537670 DOI: 10.3390/ph14101062] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2021] [Revised: 10/14/2021] [Accepted: 10/15/2021] [Indexed: 11/17/2022] Open
Abstract
We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.
Collapse
Affiliation(s)
- Tomasz Róg
- Department of Physics, University of Helsinki, 00014 Helsinki, Finland;
| | - Mykhailo Girych
- Department of Physics, University of Helsinki, 00014 Helsinki, Finland;
| | - Alex Bunker
- Drug Research Program, Division of Pharmaceutical Biosciences, Faculty of Pharmacy, University of Helsinki, 00014 Helsinki, Finland;
| |
Collapse
|
12
|
Insights into the Role of Membrane Lipids in the Structure, Function and Regulation of Integral Membrane Proteins. Int J Mol Sci 2021; 22:ijms22169026. [PMID: 34445730 PMCID: PMC8396450 DOI: 10.3390/ijms22169026] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2021] [Revised: 08/15/2021] [Accepted: 08/17/2021] [Indexed: 02/07/2023] Open
Abstract
Membrane proteins exist within the highly hydrophobic membranes surrounding cells and organelles, playing key roles in cellular function. It is becoming increasingly clear that the membrane does not just act as an appropriate environment for these proteins, but that the lipids that make up these membranes are essential for membrane protein structure and function. Recent technological advances in cryogenic electron microscopy and in advanced mass spectrometry methods, as well as the development of alternative membrane mimetic systems, have allowed experimental study of membrane protein–lipid complexes. These have been complemented by computational approaches, exploiting the ability of Molecular Dynamics simulations to allow exploration of membrane protein conformational changes in membranes with a defined lipid content. These studies have revealed the importance of lipids in stabilising the oligomeric forms of membrane proteins, mediating protein–protein interactions, maintaining a specific conformational state of a membrane protein and activity. Here we review some of the key recent advances in the field of membrane protein–lipid studies, with major emphasis on respiratory complexes, transporters, channels and G-protein coupled receptors.
Collapse
|
13
|
Ricarte A, Dalton JAR, Giraldo J. Structural Assessment of Agonist Efficacy in the μ-Opioid Receptor: Morphine and Fentanyl Elicit Different Activation Patterns. J Chem Inf Model 2021; 61:1251-1274. [PMID: 33448226 DOI: 10.1021/acs.jcim.0c00890] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Over the past two decades, the opioid epidemic in the United States and Canada has evidenced the need for a better understanding of the molecular mechanisms of medications used to fight pain. Morphine and fentanyl are widely used in opiate-mediated analgesia for the treatment of chronic pain. These compounds target the μ-opioid receptor (MOR), a class A G protein-coupled receptor (GPCR). In light of described higher efficacy of fentanyl with respect to morphine, we have performed independent μs-length unbiased molecular dynamics (MD) simulations of MOR complexes with each of these ligands, including the MOR antagonist naltrexone as a negative control. Consequently, MD simulations totaling 58 μs have been conducted to elucidate at the atomic level ligand-specific receptor activity and signal transmission in the MOR. In particular, we have identified stable binding poses of morphine and fentanyl, which interact differently with the MOR. Different ligand-receptor interaction landscapes directly induce sidechain conformational changes of orthosteric pocket residues: Asp1493.32, Tyr1503.33, Gln1262.60, and Lys2355.39. The induced conformations determine Asp1493.32-Tyr3287.43 sidechain-sidechain interactions and Trp2956.48-Ala2425.46 sidechain-backbone H-bond formations, as well as Met1533.36 conformational changes. In addition to differences in ligand binding, different intracellular receptor conformational changes are observed as morphine preferentially activates transmembrane (TM) helices: TM3 and TM5, while fentanyl preferentially activates TM6 and TM7. As conformational changes in TM6 and TM7 are widely described as being the most crucial aspect in GPCR activation, this may contribute to the greater efficacy of fentanyl over morphine. These computationally observed functional differences between fentanyl and morphine may provide new avenues for the design of safer but not weaker opioid drugs because it is desirable to increase the safety of medicines without sacrificing their efficacy.
Collapse
Affiliation(s)
- Adrián Ricarte
- Laboratory of Molecular Neuropharmacology and Bioinformatics, Unitat de Bioestadística and Institut de Neurociències, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain.,Unitat de Neurociència Traslacional, Parc Taulí Hospital Universitari, Institut d'Investigació i Innovació Parc Taulí (I3PT), Institut de Neurociències, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain.,Instituto de Salud Carlos III, Centro de Investigación Biomédica en Red de Salud Mental, CIBERSAM, 28029 Madrid, Spain
| | - James A R Dalton
- Laboratory of Molecular Neuropharmacology and Bioinformatics, Unitat de Bioestadística and Institut de Neurociències, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain.,Unitat de Neurociència Traslacional, Parc Taulí Hospital Universitari, Institut d'Investigació i Innovació Parc Taulí (I3PT), Institut de Neurociències, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain.,Instituto de Salud Carlos III, Centro de Investigación Biomédica en Red de Salud Mental, CIBERSAM, 28029 Madrid, Spain
| | - Jesús Giraldo
- Laboratory of Molecular Neuropharmacology and Bioinformatics, Unitat de Bioestadística and Institut de Neurociències, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain.,Unitat de Neurociència Traslacional, Parc Taulí Hospital Universitari, Institut d'Investigació i Innovació Parc Taulí (I3PT), Institut de Neurociències, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain.,Instituto de Salud Carlos III, Centro de Investigación Biomédica en Red de Salud Mental, CIBERSAM, 28029 Madrid, Spain
| |
Collapse
|
14
|
Renault P, Giraldo J. Dynamical Correlations Reveal Allosteric Sites in G Protein-Coupled Receptors. Int J Mol Sci 2020; 22:ijms22010187. [PMID: 33375427 PMCID: PMC7795036 DOI: 10.3390/ijms22010187] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2020] [Revised: 12/22/2020] [Accepted: 12/24/2020] [Indexed: 01/14/2023] Open
Abstract
G protein-coupled Receptors (GPCRs) play a central role in many physiological processes and, consequently, constitute important drug targets. In particular, the search for allosteric drugs has recently drawn attention, since they could be more selective and lead to fewer side effects. Accordingly, computational tools have been used to estimate the druggability of allosteric sites in these receptors. In spite of many successful results, the problem is still challenging, particularly the prediction of hydrophobic sites in the interface between the protein and the membrane. In this work, we propose a complementary approach, based on dynamical correlations. Our basic hypothesis was that allosteric sites are strongly coupled to regions of the receptor that undergo important conformational changes upon activation. Therefore, using ensembles of experimental structures, normal mode analysis and molecular dynamics simulations we calculated correlations between internal fluctuations of different sites and a collective variable describing the activation state of the receptor. Then, we ranked the sites based on the strength of their coupling to the collective dynamics. In the β2 adrenergic (β2AR), glucagon (GCGR) and M2 muscarinic receptors, this procedure allowed us to correctly identify known allosteric sites, suggesting it has predictive value. Our results indicate that this dynamics-based approach can be a complementary tool to the existing toolbox to characterize allosteric sites in GPCRs.
Collapse
Affiliation(s)
- Pedro Renault
- Laboratory of Molecular Neuropharmacology and Bioinformatics, Unitat de Bioestadística and Institut de Neurociències, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain;
- Unitat de Neurociència Traslacional, Parc Taulí Hospital Universitari, Institut d’Investigació i Innovació Parc Taulí (I3PT), Institut de Neurociències, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
- Instituto de Salud Carlos III, Centro de Investigación Biomédica en Red de Salud Mental, CIBERSAM, 08193 Bellaterra, Spain
| | - Jesús Giraldo
- Laboratory of Molecular Neuropharmacology and Bioinformatics, Unitat de Bioestadística and Institut de Neurociències, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain;
- Unitat de Neurociència Traslacional, Parc Taulí Hospital Universitari, Institut d’Investigació i Innovació Parc Taulí (I3PT), Institut de Neurociències, Universitat Autònoma de Barcelona, 08193 Bellaterra, Spain
- Instituto de Salud Carlos III, Centro de Investigación Biomédica en Red de Salud Mental, CIBERSAM, 08193 Bellaterra, Spain
- Correspondence:
| |
Collapse
|
15
|
Statistics for the analysis of molecular dynamics simulations: providing P values for agonist-dependent GPCR activation. Sci Rep 2020; 10:19942. [PMID: 33203907 PMCID: PMC7672096 DOI: 10.1038/s41598-020-77072-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2020] [Accepted: 11/05/2020] [Indexed: 02/07/2023] Open
Abstract
Molecular dynamics (MD) is the common computational technique for assessing efficacy of GPCR-bound ligands. Agonist efficacy measures the capability of the ligand-bound receptor of reaching the active state in comparison with the free receptor. In this respect, agonists, neutral antagonists and inverse agonists can be considered. A collection of MD simulations of both the ligand-bound and the free receptor are needed to provide reliable conclusions. Variability in the trajectories needs quantification and proper statistical tools for meaningful and non-subjective conclusions. Multiple-factor (time, ligand, lipid) ANOVA with repeated measurements on the time factor is proposed as a suitable statistical method for the analysis of agonist-dependent GPCR activation MD simulations. Inclusion of time factor in the ANOVA model is consistent with the time-dependent nature of MD. Ligand and lipid factors measure agonist and lipid influence on receptor activation. Previously reported MD simulations of adenosine A2a receptor (A2aR) are reanalyzed with this statistical method. TM6–TM3 and TM7–TM3 distances are selected as dependent variables in the ANOVA model. The ligand factor includes the presence or absence of adenosine whereas the lipid factor considers DOPC or DOPG lipids. Statistical analysis of MD simulations shows the efficacy of adenosine and the effect of the membrane lipid composition. Subsequent application of the statistical methodology to NECA A2aR agonist, with resulting P values in consistency with its pharmacological profile, suggests that the method is useful for ligand comparison and potentially for dynamic structure-based virtual screening.
Collapse
|
16
|
Reis RI, Moraes I. Probing Membrane Protein Assembly into Nanodiscs by In Situ Dynamic Light Scattering: A 2A Receptor as a Case Study. BIOLOGY 2020; 9:biology9110400. [PMID: 33202740 PMCID: PMC7696381 DOI: 10.3390/biology9110400] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/04/2020] [Revised: 11/02/2020] [Accepted: 11/12/2020] [Indexed: 12/20/2022]
Abstract
Simple Summary Located within the biological cell membranes, integral membrane proteins are responsible for a large variety of vital cellular processes. In humans, nearly a quarter of the genome codes integral membrane proteins, therefore malfunction of these proteins is associated with a variety of symptoms and diseases such as obesity, cancer and Parkinson’s disease. Clearly, knowledge of membrane proteins behaviour, in both structural and functional terms, is important not only in medicine but also in the design of better drugs with improved pharmaceutical properties. Nevertheless, much still remains unknown about these proteins, mainly because of the technical challenges associated with their production and stability in vitro once removed from their native lipidic environment. Recently, several membrane mimetic systems have been developed including nanodisc lipid particles. Nanodiscs are self-assembled lipidic structures that “trap” membrane proteins into a disc shaped phospholipid bilayer that is stabilised by a belt made of a protein know as membrane scaffold protein (MSP). Membrane proteins assembled into lipidic nanodiscs can maintain their structural and functional integrity and are compatible with most biophysical methods. Here we demonstrate the use of in situ dynamic light scattering as a high-throughput screening tool to assess the best conditions for nanodisc assembly and protein incorporation. Abstract Membrane proteins play a crucial role in cell physiology by participating in a variety of essential processes such as transport, signal transduction and cell communication. Hence, understanding their structure–function relationship is vital for the improvement of therapeutic treatments. Over the last decade, based on the development of detergents, amphipoles and styrene maleic-acid lipid particles (SMALPs), remarkable accomplishments have been made in the field of membrane protein structural biology. Nevertheless, there are still many drawbacks associated with protein–detergent complexes, depending on the protein in study or experimental application. Recently, newly developed membrane mimetic systems have become very popular for allowing a structural and functional characterisation of membrane proteins in vitro. The nanodisc technology is one such valuable tool, which provides a more native-like membrane environment than detergent micelles or liposomes. In addition, it is also compatible with many biophysical and biochemical methods. Here we describe the use of in situ dynamic light scattering to accurately and rapidly probe membrane proteins’ reconstitution into nanodiscs. The adenosine type 2A receptor (A2AR) was used as a case study.
Collapse
Affiliation(s)
- Rosana I. Reis
- National Physical Laboratory, Hampton Road, Teddington TW11 0LW, UK;
- Research Complex at Harwell Rutherford, Appleton Laboratory, Harwell, Oxford, Didcot OX11 0FA, UK
| | - Isabel Moraes
- National Physical Laboratory, Hampton Road, Teddington TW11 0LW, UK;
- Research Complex at Harwell Rutherford, Appleton Laboratory, Harwell, Oxford, Didcot OX11 0FA, UK
- Correspondence:
| |
Collapse
|
17
|
Torrens-Fontanals M, Stepniewski TM, Aranda-García D, Morales-Pastor A, Medel-Lacruz B, Selent J. How Do Molecular Dynamics Data Complement Static Structural Data of GPCRs. Int J Mol Sci 2020; 21:E5933. [PMID: 32824756 PMCID: PMC7460635 DOI: 10.3390/ijms21165933] [Citation(s) in RCA: 27] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2020] [Revised: 08/11/2020] [Accepted: 08/15/2020] [Indexed: 01/08/2023] Open
Abstract
G protein-coupled receptors (GPCRs) are implicated in nearly every physiological process in the human body and therefore represent an important drug targeting class. Advances in X-ray crystallography and cryo-electron microscopy (cryo-EM) have provided multiple static structures of GPCRs in complex with various signaling partners. However, GPCR functionality is largely determined by their flexibility and ability to transition between distinct structural conformations. Due to this dynamic nature, a static snapshot does not fully explain the complexity of GPCR signal transduction. Molecular dynamics (MD) simulations offer the opportunity to simulate the structural motions of biological processes at atomic resolution. Thus, this technique can incorporate the missing information on protein flexibility into experimentally solved structures. Here, we review the contribution of MD simulations to complement static structural data and to improve our understanding of GPCR physiology and pharmacology, as well as the challenges that still need to be overcome to reach the full potential of this technique.
Collapse
Affiliation(s)
- Mariona Torrens-Fontanals
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
| | - Tomasz Maciej Stepniewski
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
- InterAx Biotech AG, PARK innovAARE, 5234 Villigen, Switzerland
- Faculty of Chemistry, Biological and Chemical Research Centre, University of Warsaw, 02-093 Warsaw, Poland
| | - David Aranda-García
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
| | - Adrián Morales-Pastor
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
| | - Brian Medel-Lacruz
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
| | - Jana Selent
- Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Medical Research Institute (IMIM)—Department of Experimental and Health Sciences, Pompeu Fabra University (UPF), 08003 Barcelona, Spain; (M.T.-F.); (T.M.S.); (D.A.-G.); (A.M.-P.); (B.M.-L.)
| |
Collapse
|
18
|
Salmaso V, Jacobson KA. In Silico Drug Design for Purinergic GPCRs: Overview on Molecular Dynamics Applied to Adenosine and P2Y Receptors. Biomolecules 2020; 10:E812. [PMID: 32466404 PMCID: PMC7356333 DOI: 10.3390/biom10060812] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2020] [Revised: 05/21/2020] [Accepted: 05/22/2020] [Indexed: 12/14/2022] Open
Abstract
Molecular modeling has contributed to drug discovery for purinergic GPCRs, including adenosine receptors (ARs) and P2Y receptors (P2YRs). Experimental structures and homology modeling have proven to be useful in understanding and predicting structure activity relationships (SAR) of agonists and antagonists. This review provides an excursus on molecular dynamics (MD) simulations applied to ARs and P2YRs. The binding modes of newly synthesized A1AR- and A3AR-selective nucleoside derivatives, potentially of use against depression and inflammation, respectively, have been predicted to recapitulate their SAR and the species dependence of A3AR affinity. P2Y12R and P2Y1R crystallographic structures, respectively, have provided a detailed understanding of the recognition of anti-inflammatory P2Y14R antagonists and a large group of allosteric and orthosteric antagonists of P2Y1R, an antithrombotic and neuroprotective target. MD of A2AAR (an anticancer and neuroprotective target), A3AR, and P2Y1R has identified microswitches that are putatively involved in receptor activation. The approach pathways of different ligands toward A2AAR and P2Y1R binding sites have also been explored. A1AR, A2AAR, and A3AR were utilizes to study allosteric phenomena, but locating the binding site of structurally diverse allosteric modulators, such as an A3AR enhancer LUF6000, is challenging. Ligand residence time, a predictor of in vivo efficacy, and the structural role of water were investigated through A2AAR MD simulations. Thus, new MD and other modeling algorithms have contributed to purinergic GPCR drug discovery.
Collapse
Affiliation(s)
| | - Kenneth A. Jacobson
- Molecular Recognition Section, Laboratory of Bioorganic Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892, USA;
| |
Collapse
|