Evaluation of QSAR Equations for Virtual Screening.

Quantitative Structure Activity Relationship (QSAR) models can inform on the correlation between activities and structure-based molecular descriptors. This information is important for the understanding of the factors that govern molecular properties and for designing new compounds with favorable properties. Due to the large number of calculate-able descriptors and consequently, the much larger number of descriptors combinations, the derivation of QSAR models could be treated as an optimization problem. For continuous responses, metrics which are typically being optimized in this process are related to model performances on the training set, for example, R2 and QCV2. Similar metrics, calculated on an external set of data (e.g., QF1/F2/F32), are used to evaluate the performances of the final models. A common theme of these metrics is that they are context -" ignorant". In this work we propose that QSAR models should be evaluated based on their intended usage. More specifically, we argue that QSAR models developed for Virtual Screening (VS) should be derived and evaluated using a virtual screening-aware metric, e.g., an enrichment-based metric. To demonstrate this point, we have developed 21 Multiple Linear Regression (MLR) models for seven targets (three models per target), evaluated them first on validation sets and subsequently tested their performances on two additional test sets constructed to mimic small-scale virtual screening campaigns. As expected, we found no correlation between model performances evaluated by "classical" metrics, e.g., R2 and QF1/F2/F32 and the number of active compounds picked by the models from within a pool of random compounds. In particular, in some cases models with favorable R2 and/or QF1/F2/F32 values were unable to pick a single active compound from within the pool whereas in other cases, models with poor R2 and/or QF1/F2/F32 values performed well in the context of virtual screening. We also found no significant correlation between the number of active compounds correctly identified by the models in the training, validation and test sets. Next, we have developed a new algorithm for the derivation of MLR models by optimizing an enrichment-based metric and tested its performances on the same datasets. We found that the best models derived in this manner showed, in most cases, much more consistent results across the training, validation and test sets and outperformed the corresponding MLR models in most virtual screening tests. Finally, we demonstrated that when tested as binary classifiers, models derived for the same targets by the new algorithm outperformed Random Forest (RF) and Support Vector Machine (SVM)-based models across training/validation/test sets, in most cases. We attribute the better performances of the Enrichment Optimizer Algorithm (EOA) models in VS to better handling of inactive random compounds. Optimizing an enrichment-based metric is therefore a promising strategy for the derivation of QSAR models for classification and virtual screening.

Risk Prediction Method for Anticholinergic Action Using Auto-quantitative Structure-Activity Relationship and Docking Study with Molecular Operating Environment.

Lower urinary tract symptoms (LUTS) induced by anticholinergic drug action impair the QOL of patients and are associated with a poor prognosis. Therefore, it is expedient to develop methods of predicting the anticholinergic side effects of drugs, which we aimed to achieve in this study using a quantitative structure-activity relationship (QSAR) and docking study with molecular operations environment (MOE; Molecular Simulation Informatics Systems [MOLSIS], Inc.) In the QSAR simulation, the QSAR model built using the partial least squares regression (PLS) and genetic algorithm-multiple linear regression (GA-MLR) methods showed remarkable coefficient of determination (R2) and XR2 values. In the docking study, a specific relationship was identified between the adjusted docking score (-S) and bioactivity (pKi) values. In conclusion, the methods developed could be useful for in silico risk assessment of LUTS, and plans are potentially applicable to numerous drugs with anticholinergic activity that induce serious side effects, limiting their use.

M3 receptor modulates extracellular matrix synthesis via ERK1/2 signaling pathway in human bladder smooth muscle cells.

Extracellular matrix (ECM) accumulation plays a key role in the progression of bladder outlet obstruction (BOO). Muscarinic receptors have been widely reported to serve as pivotal regulators in lung tissue remodeling. However, the influence of them on human bladder smooth muscle cells (HBSMCs) and the underlying molecular mechanisms have not yet been evaluated. The purposes of the present study are to investigate the effect of muscarinic receptors on the synthesis of ECM in HBSMCs and the involvement of intracellular signal transducers. The results indicated that M1 -M5 muscarinic receptors were all encoded in HBSMCs. The expression rank order was M2 > M1 > M5 > M3 > M4 . The gene and protein expression of collagen I (COL1), TIMP-1, and TIMP-2 was carbachol (CCH) concentration-dependently enhanced. The synthesis of COL1 in the supernatant of cell culture medium was significantly elevated by exposure to CCH. The CCH-induced protein expression of COL1, TIMP-1, and TIMP-2, however, was obviously reduced by the pretreatment of muscarinic receptor antagonists, atropine, and M3 -preferring antagonist (1,1-dimethyl-4-diphenyl-acetoxypiperidinium iodide [4-DAMP]). Furthermore, ERK1/2 was activated by 100 µM CCH when compared with the control group and the pretreatment of ERK1/2 inhibitor significantly suppressed the synthesis of COL1 induced by 100 µM CCH. Besides, CCH-induced phosphorylation of ERK1/2 was remarkably restrained by the pretreatment of 4-DAMP. All in all, these findings demonstrated that M3 receptor can modulate extracellular matrix synthesis via the ERK1/2 signaling pathway, which may provide potential novel therapeutic targets for BOO.

Structure-Based Design and Discovery of New M2 Receptor Agonists.

Muscarinic receptor agonists are characterized by apparently strict restraints on their tertiary or quaternary amine and their distance to an ester or related center. On the basis of the active state crystal structure of the muscarinic M2 receptor in complex with iperoxo, we explored potential agonists that lacked the highly conserved functionalities of previously known ligands. Using structure-guided pharmacophore design followed by docking, we found two agonists (compounds 3 and 17), out of 19 docked and synthesized compounds, that fit the receptor well and were predicted to form a hydrogen-bond conserved among known agonists. Structural optimization led to compound 28, which was 4-fold more potent than its parent 3. Fortified by the discovery of this new scaffold, we sought a broader range of chemotypes by docking 2.2 million fragments, which revealed another three micromolar agonists unrelated either to 28 or known muscarinics. Even pockets as tightly defined and as deeply studied as that of the muscarinic reveal opportunities for the structure-based design and the discovery of new chemotypes.

M3 muscarinic acetylcholine receptor facilitates the endocytosis of mu opioid receptor mediated by morphine independently of the formation of heteromeric complexes.

Morphine inefficiency to induce the internalization of mu opioid (MOP) receptors observed in numerous experimental models constitutes a paradigm of G-protein coupled receptor (GPCR) functional selectivity. We recently described that activation of Gαq/11 proteins through 5-HT2A serotonin receptors co-expressed in the same cells facilitates MOP receptor endocytosis promoted by morphine. In order to explore whether a different Gαq/11 coupled GPCR would emulate this effect, a double stable Flp-In T-REx HEK293 cell line permanently expressing MOP-YFP receptors along with FLAG-M3-Cerulean receptors expressed in an inducible manner was generated. Fluorescence microscopy examination of these cells revealed a co-distribution of both receptors mainly compartmentalized in plasma membrane. Concurrent stimulation with carbachol and morphine promoted MOP receptor internalization, desensitization and down-regulation and this facilitation was not dependent on PKC activation. Co-immunoprecipitation experiments demonstrated that FLAG-M3-Cerulean/MOP-YFP receptors interact forming heteromeric complexes in a time depending manner, i.e. the strongest interaction was detected after 96h of FLAG-M3-Cerulean induced expression. Under these experimental conditions, treatment of cells with carbachol plus morphine resulted in the internalization of both receptors within separated endocytic vesicles as visualized by confocal microscopy. This trafficking segregation observed for FLAG-M3-Cerulean and MOP-YFP receptors upon agonist stimulation suggests that this protein-protein interaction presents temporal and dynamic properties. Moreover, MOP-YFP receptor internalization facilitated by FLAG-M3-Cerulean receptors is independent of the constitution of heteromeric complexes.

3-iodothyronamine (T1AM), a novel antagonist of muscarinic receptors.

3-iodothyronamine (T1AM) is a trace amine suspected to derive from thyroid hormone metabolism. T1AM was described as a ligand of G-protein coupled monoaminergic receptors, including trace amine associated receptors, suggesting the amine may exert a modulatory role on the monoaminergic transmission. Nothing is known on the possibility that T1AM could also modulate the cholinergic transmission interacting with muscarinic receptors. We evaluated whether T1AM (10nM-100µM) was able to i) displace [3H]-NMS (0.20nM) binding to membrane preparations from CHO cells stably transfected with human muscarinic receptor subtypes (M1-M5); ii) modify basal or acetylcholine induced pERK1/2 levels in CHO expressing the human muscarinic type 3 receptor subtype by Western blot iii) modify basal and carbachol-induced contraction of isolated rat urinary bladder. T1AM fitting within rat muscarinic type 3 receptor was simulated by Docking studies. T1AM recognized all muscarinic receptor subtypes (pKi values in the micromolar range). Interacting at type 3, T1AM reduced acetylcholine-increased pERK1/2 levels. T1AM reduced carbachol-induced contraction of the rat urinary bladder. The fenoxyl residue and the iodide ion were found essential for establishing contacts with the active site of the rat muscarinic type 3 receptor subtype. Our results indicate that T1AM binds at muscarinic receptors behaving as a weak, not selective, antagonist. This finding adds knowledge on the pharmacodynamics features of T1AM and it may prompt investigation on novel pharmacological effects of T1AM at conditions of hyper-activation of the muscarinic tone including the overactive urinary bladder.

Accelerated structure-based design of chemically diverse allosteric modulators of a muscarinic G protein-coupled receptor.

Design of ligands that provide receptor selectivity has emerged as a new paradigm for drug discovery of G protein-coupled receptors, and may, for certain families of receptors, only be achieved via identification of chemically diverse allosteric modulators. Here, the extracellular vestibule of the M2 muscarinic acetylcholine receptor (mAChR) is targeted for structure-based design of allosteric modulators. Accelerated molecular dynamics (aMD) simulations were performed to construct structural ensembles that account for the receptor flexibility. Compounds obtained from the National Cancer Institute (NCI) were docked to the receptor ensembles. Retrospective docking of known ligands showed that combining aMD simulations with Glide induced fit docking (IFD) provided much-improved enrichment factors, compared with the Glide virtual screening workflow. Glide IFD was thus applied in receptor ensemble docking, and 38 top-ranked NCI compounds were selected for experimental testing. In [(3)H]N-methylscopolamine radioligand dissociation assays, approximately half of the 38 lead compounds altered the radioligand dissociation rate, a hallmark of allosteric behavior. In further competition binding experiments, we identified 12 compounds with affinity of ≤30 µM. With final functional experiments on six selected compounds, we confirmed four of them as new negative allosteric modulators (NAMs) and one as positive allosteric modulator of agonist-mediated response at the M2 mAChR. Two of the NAMs showed subtype selectivity without significant effect at the M1 and M3 mAChRs. This study demonstrates an unprecedented successful structure-based approach to identify chemically diverse and selective GPCR allosteric modulators with outstanding potential for further structure-activity relationship studies.

Common molecular mechanism of the hepatic lesion and the cardiac parasympathetic regulation in chronic hepatitis C infection: a critical role for the muscarinic receptor type 3.

BACKGROUND: The pathophysiological overlapping between Sjorgen's Syndrome (SS) and HCV, presence of anti- muscarinic receptor type 3 (M3R) antibodies in SS, the role that M3R plays in the regulation of the heart rate, has led to the assumption that cardiovagal dysfunction in HCV patients is caused by anti-M3R antibodies elicited by HCV proteins or by their direct interaction with M3R. RESULTS: To identify HCV protein which possibly is crossreactive with M3R or which binds to this receptor, we performed the Informational Spectrum Method (ISM) analysis of the HCV proteome. This analysis revealed that NS5A protein represents the most probable interactor of M3R or that this viral protein could elicit antibodies which modulate function of this receptor. Further detailed structure/function analysis of NS5A and M3R performed by the ISM method extended with other Digital Signal processing (DSP) approaches revealed domains of these proteins which participate in their crossreactivity or in their direct interaction, representing promising diagnostic and therapeutic targets. CONCLUSIONS: Application of the ISM with other compatible bioinformatics methods offers new perspectives for identifying diagnostic and therapeutic targets for complicated forms of HCV and other viral infections. We show how the electron-ion interaction potential (EIIP) amino-acid scale used in the ISM combined with a robust, high performance hydrophobicity scale can provide new insights for understanding protein structure/function and protein-protein interactions.

Accelerated molecular dynamics simulations of ligand binding to a muscarinic G-protein-coupled receptor.

Elucidating the detailed process of ligand binding to a receptor is pharmaceutically important for identifying druggable binding sites. With the ability to provide atomistic detail, computational methods are well poised to study these processes. Here, accelerated molecular dynamics (aMD) is proposed to simulate processes of ligand binding to a G-protein-coupled receptor (GPCR), in this case the M3 muscarinic receptor, which is a target for treating many human diseases, including cancer, diabetes and obesity. Long-timescale aMD simulations were performed to observe the binding of three chemically diverse ligand molecules: antagonist tiotropium (TTP), partial agonist arecoline (ARc) and full agonist acetylcholine (ACh). In comparison with earlier microsecond-timescale conventional MD simulations, aMD greatly accelerated the binding of ACh to the receptor orthosteric ligand-binding site and the binding of TTP to an extracellular vestibule. Further aMD simulations also captured binding of ARc to the receptor orthosteric site. Additionally, all three ligands were observed to bind in the extracellular vestibule during their binding pathways, suggesting that it is a metastable binding site. This study demonstrates the applicability of aMD to protein-ligand binding, especially the drug recognition of GPCRs.

Allosteric effects of sodium ion binding on activation of the m3 muscarinic g-protein-coupled receptor.

G-protein-coupled receptors (GPCRs) are important membrane proteins that mediate cellular signaling and represent primary targets for about one-third of currently marketed drugs. Recent x-ray crystallographic studies identified distinct conformations of GPCRs in the active and inactive states. An allosteric sodium ion was found bound to a highly conserved D2.50 residue in inactive GPCRs, whereas the D2.50 allosteric pocket became collapsed in active GPCR structures. However, the dynamic mechanisms underlying these observations remain elusive. In this study, we aimed to understand the mechanistic effects of sodium ion binding on dynamic activation of the M3 muscarinic GPCR through long-timescale accelerated molecular dynamics (aMD) simulations. Results showed that with the D2.50 residue deprotonated, the M3 receptor is bound by an allosteric sodium ion and confined mostly in the inactive state with remarkably reduced flexibility. In contrast, the D2.50-protonated receptor does not exhibit sodium ion binding to the D2.50 allosteric site and samples a significantly larger conformational space. The receptor activation is captured and characterized by large-scale structural rearrangements of the transmembrane helices via dynamic hydrogen bond and salt bridge interactions. The residue motions are highly correlated during receptor activation. Further network analysis revealed that the allosteric signaling between residue D2.50 and key residues in the intracellular, extracellular, and orthosteric pockets is significantly weakened upon sodium ion binding.

Distinct Agonist Regulation of Muscarinic Acetylcholine M2-M3 Heteromers and Their Corresponding Homomers.

Each subtype of the muscarinic receptor family of G protein-coupled receptors is activated by similar concentrations of the neurotransmitter acetylcholine or closely related synthetic analogs such as carbachol. However, pharmacological selectivity can be generated by the introduction of a pair of mutations to produce Receptor Activated Solely by Synthetic Ligand (RASSL) forms of muscarinic receptors. These display loss of potency for acetylcholine/carbachol alongside a concurrent gain in potency for the ligand clozapine N-oxide. Co-expression of a form of wild type human M2 and a RASSL variant of the human M3 receptor resulted in concurrent detection of each of M2-M2 and M3-M3 homomers alongside M2-M3 heteromers at the surface of stably transfected Flp-In(TM) T-REx(TM) 293 cells. In this setting occupancy of the receptors with a muscarinic antagonist was without detectable effect on any of the muscarinic oligomers. However, selective agonist occupancy of the M2 receptor resulted in enhanced M2-M2 homomer interactions but decreased M2-M3 heteromer interactions. By contrast, selective activation of the M3 RASSL receptor did not significantly alter either M3-M3 homomer or M2-M3 heteromer interactions. Selectively targeting closely related receptor oligomers may provide novel therapeutic opportunities.

The molecular basis of oligomeric organization of the human M3 muscarinic acetylcholine receptor.

G protein-coupled receptors, including the M3 muscarinic acetylcholine receptor, can form homo-oligomers. However, the basis of these interactions and the overall organizational structure of such oligomers are poorly understood. Combinations of site-directed mutagenesis and homogenous time-resolved fluorescence resonance energy transfer studies that assessed interactions between receptor protomers at the surface of transfected cells indicated important contributions of regions of transmembrane domains I, IV, V, VI, and VII as well as intracellular helix VIII to the overall organization. Molecular modeling studies based on both these results and an X-ray structure of the inactive state of the M3 receptor bound by the antagonist/inverse agonist tiotropium were then employed. The results could be accommodated fully by models in which a proportion of the cell surface M3 receptor population is a tetramer with rhombic, but not linear, orientation. This is consistent with previous studies based on spectrally resolved, multiphoton fluorescence resonance energy transfer. Modeling studies furthermore suggest an important role for molecules of cholesterol at the dimer + dimer interface of the tetramer, which is consistent with the presence of cholesterol at key locations in many G protein-coupled receptor crystal structures. Mutants that displayed disrupted quaternary organization were often poorly expressed and showed immature N-glycosylation. Sustained treatment of cells expressing such mutants with the muscarinic receptor inverse agonist atropine increased cellular levels and restored both cell surface delivery and quaternary organization to many of the mutants. These observations suggest that organization as a tetramer may occur before plasma membrane delivery and may be a key step in cellular quality control assessment.

Helix 8 and the i3 loop of the muscarinic M3 receptor are crucial sites for its regulation by the Gß5-RGS7 complex.

The muscarinic M3 receptor (M3R) is a Gq-coupled receptor and is known to interact with many intracellular regulatory proteins. One of these molecules is Gß5-RGS7, the permanently associated heterodimer of G protein ß-subunit Gß5 and RGS7, a regulator of G protein signaling. Gß5-RGS7 can attenuate M3R-stimulated release of Ca(2+) from intracellular stores or enhance the influx of Ca(2+) across the plasma membrane. Here we show that deletion of amino acids 304-345 from the central portion of the i3 loop renders M3R insensitive to regulation by Gß5-RGS7. In addition to the i3 loop, interaction of M3R with Gß5-RGS7 requires helix 8. According to circular dichroism spectroscopy, the peptide corresponding to amino acids 548-567 in the C-terminus of M3R assumes an α-helical conformation. Substitution of Thr553 and Leu558 with Pro residues disrupts this α-helix and abolished binding to Gß5-RGS7. Introduction of the double Pro substitution into full-length M3R (M3R(TP/LP)) prevents trafficking of the receptor to the cell surface. Using atropine or other antagonists as pharmacologic chaperones, we were able to increase the level of surface expression of the TP/LP mutant to levels comparable to that of wild-type M3R. However, M3R-stimulated calcium signaling is still severely compromised. These results show that the interaction of M3R with Gß5-RGS7 requires helix 8 and the central portion of the i3 loop.

The first structure-activity relationship studies for designer receptors exclusively activated by designer drugs.

Over the past decade, two independent technologies have emerged and been widely adopted by the neuroscience community for remotely controlling neuronal activity: optogenetics which utilize engineered channelrhodopsin and other opsins, and chemogenetics which utilize engineered G protein-coupled receptors (Designer Receptors Exclusively Activated by Designer Drugs (DREADDs)) and other orthologous ligand-receptor pairs. Using directed molecular evolution, two types of DREADDs derived from human muscarinic acetylcholine receptors have been developed: hM3Dq which activates neuronal firing, and hM4Di which inhibits neuronal firing. Importantly, these DREADDs were not activated by the native ligand acetylcholine (ACh), but selectively activated by clozapine N-oxide (CNO), a pharmacologically inert ligand. CNO has been used extensively in rodent models to activate DREADDs, and although CNO is not subject to significant metabolic transformation in mice, a small fraction of CNO is apparently metabolized to clozapine in humans and guinea pigs, lessening the translational potential of DREADDs. To effectively translate the DREADD technology, the next generation of DREADD agonists are needed and a thorough understanding of structure-activity relationships (SARs) of DREADDs is required for developing such ligands. We therefore conducted the first SAR studies of hM3Dq. We explored multiple regions of the scaffold represented by CNO, identified interesting SAR trends, and discovered several compounds that are very potent hM3Dq agonists but do not activate the native human M3 receptor (hM3). We also discovered that the approved drug perlapine is a novel hM3Dq agonist with >10 000-fold selectivity for hM3Dq over hM3.

Enhancing the Reliability of GPCR Models by Accounting for Flexibility of Their Pro-Containing Helices: the Case of the Human mAChR1 Receptor.

To better investigate the GPCR structures, we have recently proposed to explore their flexibility by simulating the bending of their Pro-containing TM helices so generating a set of models (the so-called chimeras) which exhaustively combine the two conformations (bent and straight) of these helices. The primary objective of the study is to investigate whether such an approach can be exploited to enhance the reliability of the GPCR models generated by distant templates. The study was focused on the human mAChR1 receptor for which a presumably reliable model was generated using the congener mAChR3 as the template along with a second less reliable model based on the distant ß2-AR template. The second model was then utilized to produce the chimeras by combining the conformations of its Pro-containing helices (i.e., TM4, TM5, TM6 and TM7 with 16 modeled chimeras). The reliability of such chimeras was assessed by virtual screening campaigns as evaluated using a novel skewness metric where they surpassed the predictive power of the more reliable mAChR1 model. Finally, the virtual screening campaigns emphasize the opportunity of synergistically combining the scores of more chimeras using a specially developed tool which generates highly predictive consensus functions by maximizing the corresponding enrichment factors.

Muscarinic acetylcholine receptor X-ray structures: potential implications for drug development.

Muscarinic acetylcholine receptor antagonists are widely used as bronchodilating drugs in pulmonary medicine. The therapeutic efficacy of these agents depends on the blockade of M3 muscarinic receptors expressed on airway smooth muscle cells. All muscarinic antagonists currently used as bronchodilating agents show high affinity for all five muscarinic receptor subtypes, thus increasing the likelihood of unwanted side effects. Recent X-ray crystallographic studies have provided detailed structural information about the nature of the orthosteric muscarinic binding site (the conventional acetylcholine binding site) and an 'outer' receptor cavity that can bind allosteric (non-orthosteric) drugs. These new findings should guide the development of selective M3 receptor blockers that have little or no effect on other muscarinic receptor subtypes.

Novel insights into M3 muscarinic acetylcholine receptor physiology and structure.

Recent studies with M3 muscarinic acetylcholine receptor (M3R) mutant mice suggest that drugs selectively targeting this receptor subtype may prove useful for the treatment of various pathophysiological conditions. Moreover, the use of M3R-based designer G protein-coupled receptors (GPCRs) has provided novel insights into how Gq-coupled GPCRs can modulate whole-body glucose homeostasis by acting on specific peripheral cell types. More recently, we succeeded in using X-ray crystallography to determine the structure of the M3R bound to the bronchodilating drug tiotropium, a muscarinic antagonist (inverse agonist). This new structural information should facilitate the development of orthosteric or allosteric M3R-selective drugs that are predicted to have considerable therapeutic potential.

Novel structural and functional insights into M3 muscarinic receptor dimer/oligomer formation.

Class A G protein-coupled receptors (GPCRs) are able to form homodimers and/or oligomeric arrays. We recently proposed, based on bioluminescence resonance energy transfer studies with the M3 muscarinic receptor (M3R), a prototypic class A GPCR, that the M3R is able to form multiple, structurally distinct dimers that are probably transient in nature (McMillin, S. M., Heusel, M., Liu, T., Costanzi, S., and Wess, J. (2011) J. Biol. Chem. 286, 28584-28598). To provide more direct experimental support for this concept, we employed a disulfide cross-linking strategy to trap various M3R dimeric species present in a native lipid environment (transfected COS-7 cells). Disulfide cross-linking studies were carried out with many mutant M3Rs containing single cysteine (Cys) substitutions within two distinct cytoplasmic M3R regions, the C-terminal portion of the second intracellular loop (i2) and helix H8 (H8). The pattern of cross-links that we obtained, in combination with molecular modeling studies, was consistent with the existence of two structurally distinct M3R dimer interfaces, one involving i2/i2 contacts (TM4-TM5-i2 interface) and the other one characterized by H8-H8 interactions (TM1-TM2-H8 interface). Specific H8-H8 disulfide cross-links led to significant impairments in M3R-mediated G protein activation, suggesting that changes in the structural orientation or mobility of H8 are critical for efficient receptor-G protein coupling. Our findings provide novel structural and functional insights into the mechanisms involved in M3R dimerization (oligomerization). Because the M3R shows a high degree of sequence similarity with many other class A GPCRs, our findings should be of considerable general interest.

A novel experimental strategy to assess the metabolic effects of selective activation of a G(q)-coupled receptor in hepatocytes in vivo.

Increased hepatic glucose production is a key pathophysiological feature of type 2 diabetes. Like all other cell types, hepatocytes express many G protein-coupled receptors (GPCRs) that are linked to different functional classes of heterotrimeric G proteins. The important physiological functions mediated by G(s)-coupled hepatic glucagon receptors are well-documented. In contrast, little is known about the in vivo physiological roles of hepatocyte GPCRs that are linked to G proteins of the G(q) family. To address this issue, we established a transgenic mouse line (Hep-Rq mice) that expressed a G(q)-linked designer receptor (Rq) in a hepatocyte-selective fashion. Importantly, Rq could no longer bind endogenous ligands but could be selectively activated by a synthetic drug, clozapine-N-oxide. Clozapine-N-oxide treatment of Hep-Rq mice enabled us to determine the metabolic consequences caused by selective activation of a G(q)-coupled GPCR in hepatocytes in vivo. We found that acute Rq activation in vivo led to pronounced increases in blood glucose levels, resulting from increased rates of glycogen breakdown and gluconeogenesis. We also demonstrated that the expression of the V(1b) vasopressin receptor, a G(q)-coupled receptor expressed by hepatocytes, was drastically increased in livers of ob/ob mice, a mouse model of diabetes. Strikingly, treatment of ob/ob mice with a selective V(1b) receptor antagonist led to reduced glucose excursions in a pyruvate challenge test. Taken together, these findings underscore the importance of G(q)-coupled receptors in regulating hepatic glucose fluxes and suggest novel receptor targets for the treatment of type 2 diabetes.

The muscarinic M3 acetylcholine receptor exists as two differently sized complexes at the plasma membrane.

The literature on GPCR (G-protein-coupled receptor) homo-oligomerization encompasses conflicting views that range from interpretations that GPCRs must be monomeric, through comparatively newer proposals that they exist as dimers or higher-order oligomers, to suggestions that such quaternary structures are rather ephemeral or merely accidental and may serve no functional purpose. In the present study we use a novel method of FRET (Förster resonance energy transfer) spectrometry and controlled expression of energy donor-tagged species to show that M(3)Rs (muscarinic M(3) acetylcholine receptors) at the plasma membrane exist as stable dimeric complexes, a large fraction of which interact dynamically to form tetramers without the presence of trimers, pentamers, hexamers etc. That M(3)R dimeric units interact dynamically was also supported by co-immunoprecipitation of receptors synthesized at distinct times. On the basis of all these findings, we propose a conceptual framework that may reconcile the conflicting views on the quaternary structure of GPCRs.