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.

Mechanisms of Action of Cassiae Semen for Weight Management: A Computational Molecular Docking Study of Serotonin Receptor 5-HT2C.

Overweight and obesity is a growing global health concern. Current management of obesity includes lifestyle intervention, bariatric surgery and medication. The serotonin receptor, 5-HT2C, is known to mediate satiety, appetite and consumption behaviour. Lorcaserin, an appetite control drug, has demonstrated efficacy in appetite control by targeting 5-HT2C but causes undesirable side effects. This study aimed to explore the potential usage of Cassiae semen (CS), a well-known traditional Chinese medicine used to treat obesity. A computational molecular docking study was performed to determine the binding mechanism of CS compounds to the 5-HT2C receptors in both active, agonist-bound and inactive, antagonist-bound conformations. By comparing binding poses and predicted relative binding affinities towards the active or inactive forms of the receptor, we hypothesise that two of the CS compounds studied may be potent agonists which may mimic the appetite suppression effects of lorcaserin: obtusifoliol and cassiaside B2. Furthermore, two ligands, beta-sitosterol and juglanin, were predicted to bind favourably to 5-HT2C outside of the known agonist binding pocket in the active receptor, suggesting that such ligands may serve as positive allosteric modulators of 5-HT2C receptor function. Overall, this study proposed several CS compounds which may be responsible for exerting anti-obesity effects via appetite suppression by 5-HT2C receptor activation.

Structural basis of ligand recognition at the human MT1 melatonin receptor.

Melatonin (N-acetyl-5-methoxytryptamine) is a neurohormone that maintains circadian rhythms1 by synchronization to environmental cues and is involved in diverse physiological processes2 such as the regulation of blood pressure and core body temperature, oncogenesis, and immune function3. Melatonin is formed in the pineal gland in a light-regulated manner4 by enzymatic conversion from 5-hydroxytryptamine (5-HT or serotonin), and modulates sleep and wakefulness5 by activating two high-affinity G-protein-coupled receptors, type 1A (MT1) and type 1B (MT2)3,6. Shift work, travel, and ubiquitous artificial lighting can disrupt natural circadian rhythms; as a result, sleep disorders affect a substantial population in modern society and pose a considerable economic burden7. Over-the-counter melatonin is widely used to alleviate jet lag and as a safer alternative to benzodiazepines and other sleeping aids8,9, and is one of the most popular supplements in the United States10. Here, we present high-resolution room-temperature X-ray free electron laser (XFEL) structures of MT1 in complex with four agonists: the insomnia drug ramelteon11, two melatonin analogues, and the mixed melatonin-serotonin antidepressant agomelatine12,13. The structure of MT2 is described in an accompanying paper14. Although the MT1 and 5-HT receptors have similar endogenous ligands, and agomelatine acts on both receptors, the receptors differ markedly in the structure and composition of their ligand pockets; in MT1, access to the ligand pocket is tightly sealed from solvent by extracellular loop 2, leaving only a narrow channel between transmembrane helices IV and V that connects it to the lipid bilayer. The binding site is extremely compact, and ligands interact with MT1 mainly by strong aromatic stacking with Phe179 and auxiliary hydrogen bonds with Asn162 and Gln181. Our structures provide an unexpected example of atypical ligand entry for a non-lipid receptor, lay the molecular foundation of ligand recognition by melatonin receptors, and will facilitate the design of future tool compounds and therapeutic agents, while their comparison to 5-HT receptors yields insights into the evolution and polypharmacology of G-protein-coupled receptors.

The role of 5-HT2c receptor on corticosterone-mediated food intake.

Corticosterone plays an important role in feeding behavior. However, its mechanism remains unclear. Therefore, the present study aimed to investigate the effect of corticosterone on feeding behavior. In this study, cumulative food intake was increased by acute corticosterone administration in a dose-dependent manner. Administration of the 5-HT2c receptor agonist m-chlorophenylpiperazin (mCPP) reversed the effect of corticosterone on food intake. The anorectic effects of mCPP were also blocked by the 5-HT2c receptor antagonist RS102221 in corticosterone-treated mice. Both corticosterone and mCPP increased c-Fos expression in hypothalamic nuclei, but not the nucleus of the solitary tract. RS102221 inhibited c-Fos expression induced by mCPP, but not corticosterone. In addition, mCPP had little effect on TH and POMC levels in the hypothalamus. Furthermore, mCPP antagonized decreasing effect of the leptin produced by corticosterone. Taken together, our findings suggest that 5-HT2c receptors and leptin may be involved in the effects of corticosterone-induced hyperphagia.

Predictive in silico studies of human 5-hydroxytryptamine receptor subtype 2B (5-HT2B) and valvular heart disease.

Serotonin (5-hydroxytryptamine, 5-HT) receptors are neuromodulator neurotransmitter receptors which when activated trigger a signal transduction cascade within cells resulting in cell-cell communication. 5-hydroxytryptamine receptor 2B (5-HT2B) is a subtype of the seven members of 5-hydroxytrytamine receptors family which is the largest member of the super family of 7-transmembrane G-protein coupled receptors (GPCRs). Not only do 5-HT receptors play physiological roles in the cardiovascular system, gastrointestinal and endocrine function as well as the central nervous system, but they also play a role in behavioral functions. In particular 5-HT2B receptor is widely spread with regards to its distribution throughout bodily tissues and is expressed at high levels in the lungs, peripheral tissues, liver, kidneys and prostate, just to name a few. Hence 5-HT2B participates in multiple biological functions including CNS regulation, regulation of gastrointestinal motality, cardiovascular regulation and 5-HT transport system regulation. While 5-HT2B is a viable drug target and has therapeutic indications for treating obesity, psychosis, Parkinson's disease etc. there is a growing concern regarding adverse drug reactions, specifically valvulopathy associated with 5-HT2B agonists. Due to the sequence homology experienced by 5-HT2 subtypes there is also a concern regarding the off-target effects of 5-HT2A and 5-HT2C agonists. The concepts of sensitivity and subtype selectivity are of paramount importance and now can be tackled with the aid of in silico studies, especially cheminformatics, to develop models to predict valvulopathy associated toxicity of drug candidates prior to clinical trials. This review has highlighted three in silico approaches thus far that have been successful in either predicting 5-HT2B toxicity of molecules or identifying important interactions between 5-HT2B and drug molecules that bring about valvulopathy related toxicities.

PTEN-5-HT2C coupling: a new target for treating drug addiction.

It is well known that the ventral tegmental area (VTA) is a brain region in which virtually all abused drugs exert rewarding effects by activating its dopamine neurons. We recently found that the tumour suppressor enzyme phosphatase and tensin homologue deleted on chromosome 10 (PTEN) directly interacts to a region in the third intracellular loop (3L4F) of serotonin 5-HT2C receptors (5-HT2cR) in the rat VTA. PTEN limits agonist-induced 5-HT2cR phosphorylation via its protein phosphatase activity. Systemic or intra-amygdaloid application of the interfering peptide Tat-3L4F is able to disrupt PTEN coupling with 5-HT2cR in the rat VTA, resulting both in a suppression of the increased firing rate of VTA dopaminergic neurons induced by Delta 9-tetrahydrocannabinol (THC), the psychoactive ingredient of marijuana, and in a blockade of the conditioned place preference induced by THC and nicotine [Ji, S.P. et al. (2006). Nat. Med., 12: 324-329]. Because the blockade effects of Tat-3L4F peptide on the conditioned preference could be achieved by the suppression of Tat-3L4F peptide on the rewarding and/or learning/memory mechanisms associated with conditioned place preference, we recently explored whether Tat-3L4F can affect learning and memory. We observed that Tat-3L4F did not produce significant effects on spatial learning and memory in a Morris water maze test, thus indicating that Tat-3L4F can effectively suppress the rewarding effects induced by drugs of abuse.

Inhibition of serotonin 5-hydroxytryptamine2c receptor function through heterodimerization: receptor dimers bind two molecules of ligand and one G-protein.

Although dimerization appears to be a common property of G-protein-coupled receptors (GPCRs), it remains unclear whether a GPCR dimer binds one or two molecules of ligand and whether ligand binding results in activation of one or two G-proteins when measured using functional assays in intact living cells. Previously, we demonstrated that serotonin 5-hydroxytryptamine2C (5-HT2C) receptors form homodimers (Herrick-Davis, K., Grinde, E., and Mazurkiewicz, J. (2004) Biochemistry 43, 13963-13971). In the present study, an inactive 5-HT(2C) receptor was created and coexpressed with wild-type 5-HT2C receptors to determine whether dimerization regulates receptor function and to determine the ligand/dimer/G-protein stoichiometry in living cells. Mutagenesis of Ser138 to Arg (S138R) produced a 5-HT2C receptor incapable of binding ligand or stimulating inositol phosphate (IP) signaling. Confocal fluorescence imaging revealed plasma membrane expression of yellow fluorescent protein-tagged S138R receptors. Expression of wild-type 5-HT2C receptors in an S138R-expressing stable cell line had no effect on ligand binding to wild-type 5-HT2C receptors, but inhibited basal and 5-HT-stimulated IP signaling as well as constitutive and 5-HT-stimulated endocytosis of wild-type 5-HT2C receptors. M1 muscarinic receptor activation of IP production was normal in the S138R-expressing cells. Heterodimerization of S138R with wild-type 5-HT2C receptors was visualized in living cells using confocal fluorescence resonance energy transfer (FRET). FRET was dependent on the donor/acceptor ratio and independent of the receptor expression level. Therefore, inactive 5-HT2C receptors inhibit wild-type 5-HT2C receptor function by forming nonfunctional heterodimers expressed on the plasma membrane. These results are consistent with a model in which one GPCR dimer binds two molecules of ligand and one G-protein and indicate that dimerization is essential for 5-HT receptor function.

Novel, highly potent, selective 5-HT2A/D2 receptor antagonists as potential atypical antipsychotics.

The discovery of N-substituted-pyridoindolines and their binding affinities at the 5-HT(2A), 5-HT(2C) and D(2) receptors, and in vivo efficacy as 5-HT(2A) antagonists is described. The structure-activity relationship of a series of core tetracyclic derivatives with varying butyrophenone sidechains is also discussed. This study has led to the identification of potent, orally bioavailable 5-HT(2A)/D(2) receptor dual antagonists as potential atypical antipsychotics.