Difference and Influence of Inactive and Active States of Cannabinoid Receptor Subtype CB2: From Conformation to Drug Discovery

Cannabinoid receptor 2 (CB2), a G protein-coupled receptor (GPCR), is a promising target for the treatment of neuropathic pain, osteoporosis, immune system, cancer, and drug abuse. The lack of an experimental three-dimensional CB2 structure has hindered not only the development of studies of conformational differences between the inactive and active CB2 but also the rational discovery of novel functional compounds targeting CB2. In this work, we constructed models of both inactive and active CB2 by homology modeling. Then we conducted two comparative 100 ns molecular dynamics (MD) simulations on the two systems-the active CB2 bound with both the agonist and G protein and the inactive CB2 bound with inverse agonist-to analyze the conformational difference of CB2 proteins and the key residues involved in molecular recognition. Our results showed that the inactive CB2 and the inverse agonist remained stable during the MD simulation. However, during the MD simulations, we observed dynamical details about the breakdown of the "ionic lock" between R131(3.50) and D240(6.30) as well as the outward/inward movements of transmembrane domains of the active CB2 that bind with G proteins and agonist (TM5, TM6, and TM7). All of these results are congruent with the experimental data and recent reports. Moreover, our results indicate that W258(6.48) in TM6 and residues in TM4 (V164(4.56)-L169(4.61)) contribute greatly to the binding of the agonist on the basis of the binding energy decomposition, while residues S180-F183 in extracellular loop 2 (ECL2) may be of importance in recognition of the inverse agonist. Furthermore, pharmacophore modeling and virtual screening were carried out for the inactive and active CB2 models in parallel. Among all 10 hits, two compounds exhibited novel scaffolds and can be used as novel chemical probes for future studies of CB2. Importantly, our studies show that the hits obtained from the inactive CB2 model mainly act as inverse agonist(s) or neutral antagonist(s) at low concentration. Moreover, the hit from the active CB2 model also behaves as a neutral antagonist at low concentration. Our studies provide new insight leading to a better understanding of the structural and conformational differences between two states of CB2 and illuminate the effects of structure on virtual screening and drug design.

[1]  K. Palczewski,et al.  Crystal Structure of Rhodopsin: A G‐Protein‐Coupled Receptor , 2002, Chembiochem : a European journal of chemical biology.

[2]  Zhiwei Feng,et al.  Design and activity of AP endonuclease-1 inhibitors , 2015, Journal of chemical biology.

[3]  Daniel Fleischer,et al.  Mutation Studies of Ser7.39 and Ser2.60 in the Human CB1 Cannabinoid Receptor: Evidence for a Serine-Induced Bend in CB1 Transmembrane Helix 7 , 2007, Molecular Pharmacology.

[4]  Peng Yang,et al.  Structural Insight into Tetrameric hTRPV1 from Homology Modeling, Molecular Docking, Molecular Dynamics Simulation, Virtual Screening, and Bioassay Validations , 2015, J. Chem. Inf. Model..

[5]  Tiziano Tuccinardi,et al.  Cannabinoid CB2/CB1 selectivity. Receptor modeling and automated docking analysis. , 2006, Journal of medicinal chemistry.

[6]  Alessandro Pedretti,et al.  VEGA: a versatile program to convert, handle and visualize molecular structure on Windows-based PCs. , 2002, Journal of molecular graphics & modelling.

[7]  J. Thornton,et al.  PROCHECK: a program to check the stereochemical quality of protein structures , 1993 .

[8]  W. C. Still,et al.  Approximate atomic surfaces from linear combinations of pairwise overlaps (LCPO) , 1999 .

[9]  B. Kobilka,et al.  Energy landscapes as a tool to integrate GPCR structure, dynamics, and function. , 2010, Physiology.

[10]  Ajay N. Jain Scoring noncovalent protein-ligand interactions: A continuous differentiable function tuned to compute binding affinities , 1996, J. Comput. Aided Mol. Des..

[11]  B. Kobilka,et al.  The role of protein dynamics in GPCR function: insights from the β2AR and rhodopsin. , 2014, Current opinion in cell biology.

[12]  Krzysztof Palczewski,et al.  The Significance of G Protein-Coupled Receptor Crystallography for Drug Discovery , 2011, Pharmacological Reviews.

[13]  Alexandros Makriyannis,et al.  Pharmacotherapeutic targeting of the endocannabinoid signaling system: Drugs for obesity and the metabolic syndrome , 2008, Physiology & Behavior.

[14]  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.

[15]  John W Huffman,et al.  Recent advances in the development of selective ligands for the cannabinoid CB(2) receptor. , 2008, Current topics in medicinal chemistry.

[16]  Jan H. Jensen,et al.  Improved Treatment of Ligands and Coupling Effects in Empirical Calculation and Rationalization of pKa Values. , 2011, Journal of chemical theory and computation.

[17]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[18]  Tingjun Hou,et al.  Studies on the Interactions between β2 Adrenergic Receptor and Gs Protein by Molecular Dynamics Simulations , 2012, J. Chem. Inf. Model..

[19]  S. Sprang,et al.  The structure of the G protein heterotrimer Giα1 β 1 γ 2 , 1995, Cell.

[20]  D. Case,et al.  Theory and applications of the generalized born solvation model in macromolecular simulations , 2000, Biopolymers.

[21]  S. Rasmussen,et al.  Crystal Structure of the β2Adrenergic Receptor-Gs protein complex , 2011, Nature.

[22]  Cody J. Wenthur,et al.  Drugs for allosteric sites on receptors. , 2014, Annual review of pharmacology and toxicology.

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

[24]  A. Sali,et al.  Comparative protein structure modeling of genes and genomes. , 2000, Annual review of biophysics and biomolecular structure.

[25]  Zheng-Xiong Xi,et al.  Cannabinoid CB2 receptors modulate midbrain dopamine neuronal activity and dopamine-related behavior in mice , 2014, Proceedings of the National Academy of Sciences.

[26]  D. Case,et al.  Generalized born models of macromolecular solvation effects. , 2000, Annual review of physical chemistry.

[27]  Laxmikant V. Kale,et al.  NAMD2: Greater Scalability for Parallel Molecular Dynamics , 1999 .

[28]  Patricia H. Reggio,et al.  The Difference between the CB1 and CB2Cannabinoid Receptors at Position 5.46 Is Crucial for the Selectivity of WIN55212-2 for CB2 , 1999 .

[29]  Youyong Li,et al.  Discovery of novel inhibitors targeting the macrophage migration inhibitory factor via structure-based virtual screening and bioassays. , 2014, Journal of medicinal chemistry.

[30]  Xiang-Qun Xie,et al.  Beta-caryophyllene is a dietary cannabinoid , 2008, Proceedings of the National Academy of Sciences.

[31]  Junmei Wang,et al.  GPCR Structure-Based Virtual Screening Approach for CB2 Antagonist Search , 2007, J. Chem. Inf. Model..

[32]  Xiang-Qun Xie,et al.  Alkylamides from Echinacea Are a New Class of Cannabinomimetics , 2006, Journal of Biological Chemistry.

[33]  Zhiwei Feng,et al.  Selectivity and activation of dopamine D3R from molecular dynamics , 2012, Journal of Molecular Modeling.

[34]  P. Casellas,et al.  Signaling pathway associated with stimulation of CB2 peripheral cannabinoid receptor. Involvement of both mitogen-activated protein kinase and induction of Krox-24 expression. , 1996, European journal of biochemistry.

[35]  P Ferrara,et al.  Mutational analysis and molecular modelling of the antagonist SR 144528 binding site on the human cannabinoid CB(2) receptor. , 2000, European journal of pharmacology.

[36]  R. Stevens,et al.  High-Resolution Crystal Structure of an Engineered Human β2-Adrenergic G Protein–Coupled Receptor , 2007, Science.

[37]  Ichio Shimada,et al.  Efficacy of the β2-adrenergic receptor is determined by conformational equilibrium in the transmembrane region , 2012, Nature Communications.

[38]  Dow P Hurst,et al.  Residues accessible in the binding-site crevice of transmembrane helix 6 of the CB2 cannabinoid receptor. , 2008, Biochemistry.

[39]  Tingjun Hou,et al.  Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics Simulations , 2011, J. Chem. Inf. Model..

[40]  M. Babu,et al.  Molecular signatures of G-protein-coupled receptors , 2013, Nature.

[41]  Danfeng Shi,et al.  Ligand induced change of β2 adrenergic receptor from active to inactive conformation and its implication for the closed/open state of the water channel: insight from molecular dynamics simulation, free energy calculation and Markov state model analysis. , 2014, Physical chemistry chemical physics : PCCP.

[42]  F Guarnieri,et al.  Activation of the cannabinoid CB1 receptor may involve a W6 48/F3 36 rotamer toggle switch. , 2002, The journal of peptide research : official journal of the American Peptide Society.

[43]  John E. Stone,et al.  Using VMD: An Introductory Tutorial , 2008, Current protocols in bioinformatics.

[44]  Xiang-Qun Xie,et al.  3D structural model of the G‐protein‐coupled cannabinoid CB2 receptor , 2003, Proteins.

[45]  A. Poso,et al.  Targeting the cannabinoid CB2 receptor: modelling and structural determinants of CB2 selective ligands , 2008, British journal of pharmacology.

[46]  Alan Grossfield,et al.  Structural Basis of G Protein-coupled Receptor-Gi Protein Interaction , 2014, The Journal of Biological Chemistry.

[47]  Manfred J. Sippl,et al.  Thirty years of environmental health research--and growing. , 1996, Nucleic Acids Res..

[48]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[49]  Xiang-Qun Xie,et al.  Mutagenesis and computer modeling studies of a GPCR conserved residue W5.43(194) in ligand recognition and signal transduction for CB2 receptor. , 2011, International immunopharmacology.

[50]  Nagarajan Vaidehi,et al.  Agonist-induced conformational changes in bovine rhodopsin: insight into activation of G-protein-coupled receptors. , 2008, Journal of molecular biology.

[51]  Dan Li,et al.  Exploring the binding mechanisms of MIF to CXCR2 using theoretical approaches. , 2015, Physical chemistry chemical physics : PCCP.

[52]  Dan Li,et al.  Binding mechanisms of 1,4-dihydropyridine derivatives to L-type calcium channel Cav1.2: a molecular modeling study. , 2016, Molecular bioSystems.

[53]  Maxime Louet,et al.  GDP Release Preferentially Occurs on the Phosphate Side in Heterotrimeric G-proteins , 2012, PLoS Comput. Biol..

[54]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[55]  Pankaj R Daga,et al.  Homology modeling and molecular dynamics simulations of the active state of the nociceptin receptor reveal new insights intoagonist binding and activation , 2012, Proteins.

[56]  Qin Tong,et al.  Novel triaryl sulfonamide derivatives as selective cannabinoid receptor 2 inverse agonists and osteoclast inhibitors: discovery, optimization, and biological evaluation. , 2013, Journal of medicinal chemistry.

[57]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[58]  Qin Tong,et al.  Lead discovery, chemistry optimization, and biological evaluation studies of novel biamide derivatives as CB2 receptor inverse agonists and osteoclast inhibitors. , 2012, Journal of medicinal chemistry.

[59]  Peng Yang,et al.  Modeling, Molecular Dynamics Simulation, and Mutation Validation for Structure of Cannabinoid Receptor 2 Based on Known Crystal Structures of GPCRs , 2014, J. Chem. Inf. Model..

[60]  Philippe Chavatte,et al.  Novel 4-oxo-1,4-dihydroquinoline-3-carboxamide derivatives as new CB2 cannabinoid receptors agonists: synthesis, pharmacological properties and molecular modeling. , 2006, Journal of medicinal chemistry.

[61]  B. Wallace,et al.  HOLE: a program for the analysis of the pore dimensions of ion channel structural models. , 1996, Journal of molecular graphics.

[62]  K. Sharp,et al.  Accurate Calculation of Hydration Free Energies Using Macroscopic Solvent Models , 1994 .

[63]  Tingjun Hou,et al.  Molecular dynamics and free energy studies on the wild-type and double mutant HIV-1 protease complexed with amprenavir and two amprenavir-related inhibitors: mechanism for binding and drug resistance. , 2007, Journal of medicinal chemistry.