GPR6 Structural Insights: Homology Model Construction and Docking Studies

GPR6 is an orphan G protein-coupled receptor that has been associated with the cannabinoid family because of its recognition of a sub-set of cannabinoid ligands. The high abundance of GPR6 in the central nervous system, along with high constitutive activity and a link to several neurodegenerative diseases make GPR6 a promising biological target. In fact, diverse research groups have demonstrated that GPR6 represents a possible target for the treatment of neurodegenerative disorders such as Parkinson’s disease, Alzheimer’s disease, and Huntington’s disease. Several patents have claimed the use of a wide range of pyrazine derivatives as GPR6 inverse agonists for the treatment of Parkinson’s disease symptoms and other dyskinesia syndromes. However, the full pharmacological importance of GPR6 has not yet been fully explored due to the lack of high potency, readily available ligands targeting GPR6. The long-term goal of the present study is to develop such ligands. In this paper, we describe our initial steps towards this goal. A human GPR6 homology model was constructed using a suite of computational techniques. This model permitted the identification of unique GPR6 structural features and the exploration of the GPR6 binding crevice. A subset of patented pyrazine analogs were docked in the resultant GPR6 inactive state model to validate the model, rationalize the structure-activity relationships from the reported patents and identify the key residues in the binding crevice for ligand recognition. We will take this structural knowledge into the next phase of GPR6 project, in which scaffold hopping will be used to design new GPR6 ligands.

[1]  Alexander D. MacKerell,et al.  CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field , 2015, Journal of chemical theory and computation.

[2]  Patrick Oeckl,et al.  G-protein coupled receptor 6 deficiency alters striatal dopamine and cAMP concentrations and reduces dyskinesia in a mouse model of Parkinson's disease , 2014, Experimental Neurology.

[3]  G Vassart,et al.  Molecular cloning of an orphan G-protein-coupled receptor that constitutively activates adenylate cyclase. , 1995, The Biochemical journal.

[4]  A. Tenner,et al.  C1q-induced LRP1B and GPR6 Proteins Expressed Early in Alzheimer Disease Mouse Models, Are Essential for the C1q-mediated Protection against Amyloid-β Neurotoxicity* , 2012, The Journal of Biological Chemistry.

[5]  Leonardo Pardo,et al.  Serine and Threonine Residues Bend α-Helices in the χ1 = g− Conformation , 2000 .

[6]  Hyeon Joo,et al.  OPM database and PPM web server: resources for positioning of proteins in membranes , 2011, Nucleic Acids Res..

[7]  Tudor I. Oprea,et al.  A comprehensive map of molecular drug targets , 2016, Nature Reviews Drug Discovery.

[8]  Z. Song,et al.  Novel inverse agonists for the orphan G protein-coupled receptor 6 , 2018, Heliyon.

[9]  Mary Kay Lobo,et al.  Genetic control of instrumental conditioning by striatopallidal neuron–specific S1P receptor Gpr6 , 2007, Nature Neuroscience.

[10]  David S. Goodsell,et al.  The RCSB Protein Data Bank: redesigned web site and web services , 2010, Nucleic Acids Res..

[11]  W. Hagmann,et al.  Conformational analysis and receptor docking of N-[(1S,2S)-3-(4-chlorophenyl)-2-(3-cyanophenyl)-1-methylpropyl]-2-methyl-2-{[5-(trifluoromethyl)pyridin-2-yl]oxy}propanamide (taranabant, MK-0364), a novel, acyclic cannabinoid-1 receptor inverse agonist. , 2008, Journal of medicinal chemistry.

[12]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[13]  Tod D Romo,et al.  A Lipid Pathway for Ligand Binding Is Necessary for a Cannabinoid G Protein-coupled Receptor* , 2010, The Journal of Biological Chemistry.

[14]  Hugh Rosen,et al.  Crystal Structure of a Lipid G Protein–Coupled Receptor , 2012, Science.

[15]  A. Sali,et al.  Modeling of loops in protein structures , 2000, Protein science : a publication of the Protein Society.

[16]  Leonardo Pardo,et al.  Ser and Thr Residues Modulate the Conformation of Pro-Kinked Transmembrane α-Helices , 2004 .

[17]  R. Aronstam,et al.  Constitutive Activity among Orphan Class-A G Protein Coupled Receptors , 2015, PloS one.

[18]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 1. Generalized Born , 2012, Journal of chemical theory and computation.

[19]  Evangelia Kotsikorou,et al.  The importance of hydrogen bonding and aromatic stacking to the affinity and efficacy of cannabinoid receptor CB2 antagonist, 5-(4-chloro-3-methylphenyl)-1-[(4-methylphenyl)methyl]-N-[(1S,2S,4R)-1,3,3-trimethylbicyclo[2.2.1]hept-2-yl]-1H-pyrazole-3-carboxamide (SR144528). , 2013, Journal of medicinal chemistry.

[20]  H. Schiöth,et al.  The Repertoire of G-Protein–Coupled Receptors in Fully Sequenced Genomes , 2005, Molecular Pharmacology.

[21]  Ben M. Webb,et al.  Comparative Protein Structure Modeling Using Modeller , 2006, Current protocols in bioinformatics.

[22]  Patricia H. Reggio,et al.  Computational methods in drug design: Modeling G protein-coupled receptor monomers, dimers, and oligomers , 2006, The AAPS Journal.

[23]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[24]  H. Schiöth,et al.  The G-protein-coupled receptors in the human genome form five main families. Phylogenetic analysis, paralogon groups, and fingerprints. , 2003, Molecular pharmacology.

[25]  J. Olson,et al.  Regional and cellular gene expression changes in human Huntington's disease brain. , 2006, Human molecular genetics.

[26]  B. Honig,et al.  A hierarchical approach to all‐atom protein loop prediction , 2004, Proteins.

[27]  J. Ballesteros,et al.  Analysis and refinement of criteria for predicting the structure and relative orientations of transmembranal helical domains. , 1992, Biophysical journal.

[28]  J. Ballesteros,et al.  [19] Integrated methods for the construction of three-dimensional models and computational probing of structure-function relations in G protein-coupled receptors , 1995 .

[29]  Patrick Oeckl,et al.  Increased susceptibility of G-protein coupled receptor 6 deficient mice to MPTP neurotoxicity , 2016, Neuroscience.

[30]  T. Blundell,et al.  Comparative protein modelling by satisfaction of spatial restraints. , 1993, Journal of molecular biology.

[31]  J. Ballesteros,et al.  Activation of the β2-Adrenergic Receptor Involves Disruption of an Ionic Lock between the Cytoplasmic Ends of Transmembrane Segments 3 and 6* , 2001, The Journal of Biological Chemistry.

[32]  Robert Abel,et al.  OPLS3e: Extending Force Field Coverage for Drug-Like Small Molecules. , 2019, Journal of chemical theory and computation.

[33]  Alexander D. MacKerell,et al.  Importance of the CMAP correction to the CHARMM22 protein force field: dynamics of hen lysozyme. , 2006, Biophysical journal.

[34]  Jie Yin,et al.  High-resolution crystal structure of the human CB1 cannabinoid receptor , 2016, Nature.

[35]  G. Pasternak,et al.  Mechanistic Insights into the Allosteric Modulation of Opioid Receptors by Sodium Ions , 2014, Biochemistry.

[36]  Robert J. Lefkowitz,et al.  Transduction of Receptor Signals by ß-Arrestins , 2005, Science.

[37]  Matthew P. Repasky,et al.  Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. , 2004, Journal of medicinal chemistry.

[38]  Z. Song,et al.  GPR3 and GPR6, novel molecular targets for cannabidiol. , 2017, Biochemical and biophysical research communications.

[39]  Alexander D. MacKerell,et al.  Update of the CHARMM all-atom additive force field for lipids: validation on six lipid types. , 2010, The journal of physical chemistry. B.

[40]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 2. Explicit Solvent Particle Mesh Ewald. , 2013, Journal of chemical theory and computation.

[41]  T. Bonner,et al.  Molecular cloning of a novel candidate G protein‐coupled receptor from rat brain , 1994, FEBS letters.

[42]  Evi Kostenis,et al.  Sphingosine 1-phosphate is a ligand of the human gpr3, gpr6 and gpr12 family of constitutively active G protein-coupled receptors. , 2002, Cellular signalling.

[43]  P. Morales,et al.  Towards a better understanding of the cannabinoid-related orphan receptors GPR3, GPR6, and GPR12 , 2018, Drug metabolism reviews.

[44]  P. Morales,et al.  Methods for the Development of In Silico GPCR Models. , 2017, Methods in enzymology.

[45]  Frank Guarnieri,et al.  Conformational memories and a simulated annealing program that learns: Application to LTB4 , 1995, J. Comput. Chem..

[46]  M. Wheatley,et al.  Glycosylation of G-protein-coupled receptors for hormones central to normal reproductive functioning: its occurrence and role. , 1999, Human reproduction update.

[47]  B. L. de Groot,et al.  CHARMM36m: an improved force field for folded and intrinsically disordered proteins , 2016, Nature Methods.

[48]  D. Engelman,et al.  Sequence motifs, polar interactions and conformational changes in helical membrane proteins. , 2003, Current opinion in structural biology.

[49]  J Andrew McCammon,et al.  Allosteric effects of sodium ion binding on activation of the m3 muscarinic g-protein-coupled receptor. , 2015, Biophysical journal.

[50]  Eric Trinquet,et al.  Structural insights into biased G protein-coupled receptor signaling revealed by fluorescence spectroscopy , 2012, Proceedings of the National Academy of Sciences.

[51]  K. Merienne,et al.  Altered enhancer transcription underlies Huntington’s disease striatal transcriptional signature , 2017, Scientific Reports.

[52]  Jennifer L. Knight,et al.  OPLS3: A Force Field Providing Broad Coverage of Drug-like Small Molecules and Proteins. , 2016, Journal of chemical theory and computation.

[53]  R. Stevens,et al.  Crystal Structure of the Human Cannabinoid Receptor CB1 , 2016, Cell.

[54]  Z. Song,et al.  Discovery of endogenous inverse agonists for G protein-coupled receptor 6. , 2019, Biochemical and biophysical research communications.

[55]  Leonardo Pardo,et al.  Structural models of class a G protein-coupled receptors as a tool for drug design: insights on transmembrane bundle plasticity. , 2007, Current topics in medicinal chemistry.

[56]  H. Schaller,et al.  Sphingosine-1-phosphate is a high-affinity ligand for the G protein-coupled receptor GPR6 from mouse and induces intracellular Ca2+ release by activating the sphingosine-kinase pathway. , 2003, Biochemical and biophysical research communications.

[57]  Vadim Cherezov,et al.  Allosteric sodium in class A GPCR signaling. , 2014, Trends in biochemical sciences.

[58]  Yasuyuki Kihara,et al.  Crystal Structure of Antagonist Bound Human Lysophosphatidic Acid Receptor 1 , 2015, Cell.

[59]  A. Thathiah,et al.  The role of GPCRs in neurodegenerative diseases: avenues for therapeutic intervention , 2017, Current opinion in pharmacology.

[60]  R. Friesner,et al.  The VSGB 2.0 model: A next generation energy model for high resolution protein structure modeling , 2011, Proteins.

[61]  Nagarajan Vaidehi,et al.  Position of helical kinks in membrane protein crystal structures and the accuracy of computational prediction. , 2009, Journal of molecular graphics & modelling.

[62]  A. Makriyannis,et al.  Novel Insights into CB1 Cannabinoid Receptor Signaling: A Key Interaction Identified between the Extracellular-3 Loop and Transmembrane Helix 2 , 2013, The Journal of Pharmacology and Experimental Therapeutics.

[63]  David Rodríguez,et al.  The role of a sodium ion binding site in the allosteric modulation of the A(2A) adenosine G protein-coupled receptor. , 2013, Structure.

[64]  H. Horiuchi.  Seven-transmembrane receptors , 2015 .

[65]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[66]  H. Weinstein,et al.  Conformational Memories and the Exploration of Biologically Relevant Peptide Conformations: An Illustration for the Gonadotropin-Releasing Hormone , 1996 .

[67]  W. Sherman,et al.  Prediction of Absolute Solvation Free Energies using Molecular Dynamics Free Energy Perturbation and the OPLS Force Field. , 2010, Journal of chemical theory and computation.

[68]  H Weinstein,et al.  Prokink: a protocol for numerical evaluation of helix distortions by proline. , 2000, Protein engineering.

[69]  W. Im,et al.  Automated Builder and Database of Protein/Membrane Complexes for Molecular Dynamics Simulations , 2007, PloS one.

[70]  C. Breivogel,et al.  Investigation of non-CB1, non-CB2 WIN55212-2-sensitive G-protein-coupled receptors in the brains of mammals, birds, and amphibians , 2018, Journal of receptor and signal transduction research.

[71]  D. Hurst,et al.  Identification of Crucial Amino Acid Residues Involved in Agonist Signaling at the GPR55 Receptor. , 2017, Biochemistry.

[72]  F. Cordes,et al.  Proline-induced distortions of transmembrane helices. , 2002, Journal of molecular biology.

[73]  Adam J Pawson,et al.  Class A Orphans (version 2019.5) in the IUPHAR/BPS Guide to Pharmacology Database , 2019, IUPHAR/BPS Guide to Pharmacology CITE.

[74]  Mohit Kumar,et al.  Novel Structural Insights into GPCR-β-Arrestin Interaction and Signaling. , 2017, Trends in cell biology.

[75]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[76]  Kurt Wüthrich,et al.  Biased Signaling Pathways in β2-Adrenergic Receptor Characterized by 19F-NMR , 2012, Science.

[77]  Mark S.P. Sansom,et al.  Hinges, swivels and switches: the role of prolines in signalling via transmembrane α-helices , 2000 .