Abstract. Allosteric modulators of G protein-coupled receptors (GPCRs), which target at allosteric sites, have significant advantages against the corresponding orthosteric compounds including higher selectivity, improved chemical tractability or physicochemical properties, and reduced risk of receptor oversensitization. Bitopic ligands of GPCRs target both orthosteric and allosteric sites. Bitopic ligands can improve binding affinity, enhance subtype selectivity, stabilize receptors, and reduce side effects. Discovering allosteric modulators or bitopic ligands for GPCRs has become an emerging research area, in which the design of allosteric modulators is a key step in the detection of bitopic ligands. Radioligand binding and functional assays ([ 35 S]GTPγS and ERK1/2 phosphorylation) are used to test the effects for potential modulators or bitopic ligands. High-throughput screening (HTS) in combination with disulfide trapping and fragment-based screening are used to aid the discovery of the allosteric modulators or bitopic ligands of GPCRs. When used alone, these methods are costly and can often result in too many potential drug targets, including false positives. Alternatively, low-cost and efficient computational approaches are useful in drug discovery of novel allosteric modulators and bitopic ligands to help refine the number of targets and reduce the false-positive rates. This review summarizes the state-of-the-art computational methods for the discovery of modulators and bitopic ligands. The challenges and opportunities for future drug discovery are also discussed.

Abstract. Metabotropic glutamate receptors (mGluR) are mainly expressed in the central nervous system (CNS) and contain eight receptor subtypes, named mGluR1 to mGluR8. The crystal structures of mGluR1 and mGluR5 that are bound with the negative allosteric modulator (NAM) were reported recently. These structures provide a basic model for all class C of G-protein coupled receptors (GPCRs) and may aid in the design of new allosteric modulators for the treatment of CNS disorders. However, these structures are only combined with NAMs in the previous reports. The conformations that are bound with positive allosteric modulator (PAM) or agonist of mGluR1/5 remain unknown. Moreover, the structural information of the other six mGluRs and the comparisons of the mGluRs family have not been explored in terms of their binding pockets, the binding modes of different compounds, and important binding residues. With these crystal structures as the starting point, we built 3D structural models for six mGluRs by using homology modeling and molecular dynamics (MD) simulations. We systematically compared their allosteric binding sites/pockets, the important residues, and the selective residues by using a series of comparable dockings with both the NAM and the PAM. Our results show that several residues played important roles for the receptors' selectivity. The observations of detailed interactions between compounds and their correspondent receptors are congruent with the specificity and potency of derivatives or compounds bioassayed in vitro. We then carried out 100 ns MD simulations of mGluR5 (residue 26-832, formed by Venus Flytrap domain, a so-called cysteine-rich domain, and 7 transmembrane domains) bound with antagonist/NAM and with agonist/PAM. Our results show that both the NAM and the PAM seemed stable in class C GPCRs during the MD. However, the movements of Bionic lock,^of trans-membrane domains, and of some activation-related residues in 7 trans-membrane domains of mGluR5 were congruent with the findings in class A GPCRs. Finally, we selected nine representative bound structures to perform 30 ns MD simulations for validating the stabilities of interactions, respectively. All these bound structures kept stable during the MD simulations, indicating that the binding poses in this present work are reasonable. We provided new insight into better understanding of the structural and functional roles of the mGluRs family and facilitated the future structure-based design of novel ligands of mGluRs family with therapeutic potential.

Drug abuse is a serious problem worldwide. Recently, hallucinogens have been reported as a potential preventative and auxiliary therapy for substance abuse. However, the use of hallucinogens as a drug abuse treatment has potential risks, as the fundamental mechanisms of hallucinogens are not clear. So far, no scientific database is available for the mechanism research of hallucinogens. We constructed a hallucinogen-specific chemogenomics database by collecting chemicals, protein targets and pathways closely related to hallucinogens. This information, together with our established computational chemogenomics tools, such as TargetHunter and HTDocking, provided a one-step solution for the mechanism study of hallucinogens. We chose salvinorin A, a potent hallucinogen extracted from the plant Salvia divinorum, as an example to demonstrate the usability of our platform. With the help of HTDocking program, we predicted four novel targets for salvinorin A, including muscarinic acetylcholine receptor 2, cannabinoid receptor 1, cannabinoid receptor 2 and dopamine receptor 2. We looked into the interactions between salvinorin A and the predicted targets. The binding modes, pose and docking scores indicate that salvinorin A may interact with some of these predicted targets. Overall, our database enriched the information of systems pharmacological analysis, target identification and drug discovery for hallucinogens.