Structure-based virtual screening of MT2 melatonin receptor: influence of template choice and structural refinement

Melatonin receptor structure and signaling

Melatonin (5-methoxy-N-acetyltryptamine) binds with high affinity and specificity to membrane receptors. Several receptor subtypes exist in different species, of which the mammalian MT1 and MT2 receptors are the best-characterized. They are members of the G protein-coupled receptor superfamily, preferentially coupling to Gi/o proteins but also to other G proteins in a cell-context-depending manner. In this review, experts on melatonin receptors will summarize the current state of the field. We briefly report on the discovery and classification of melatonin receptors, then focus on the molecular structure of human MT1 and MT2 receptors and highlight the importance of molecular simulations to identify new ligands and to understand the structural dynamics of these receptors. We then describe the state-of-the-art of the intracellular signaling pathways activated by melatonin receptors and their complexes. Brief statements on the molecular toolbox available for melatonin receptor studies and future perspectives will round-up this review.

How good are AlphaFold models for docking-based virtual screening?

A crucial component in structure-based drug discovery is the availability of high-quality three-dimensional structures of the protein target. Whenever experimental structures were not available, homology modeling has been, so far, the method of choice. Recently, AlphaFold (AF), an artificial-intelligence-based protein structure prediction method, has shown impressive results in terms of model accuracy. This outstanding success prompted us to evaluate how accurate AF models are from the perspective of docking-based drug discovery. We compared the high-throughput docking (HTD) performance of AF models with their corresponding experimental PDB structures using a benchmark set of 22 targets. The AF models showed consistently worse performance using four docking programs and two consensus techniques. Although AlphaFold shows a remarkable ability to predict protein architecture, this might not be enough to guarantee that AF models can be reliably used for HTD, and post-modeling refinement strategies might be key to increase the chances of success.

In silico drug discovery of melatonin receptor ligands with therapeutic potential

INTRODUCTION

The neurohormone melatonin (N-acetyl-5-methoxytryptamine) regulates circadian rhythms exerting a variety of effects in the central nervous system and in periphery. These activities are mainly mediated by activation of MT₁ and MT₂ GPCRs. MT₁/MT₂ agonist compounds are used clinically for insomnia, depression, and circadian rhythm disturbances.

AREA COVERED

The following review describes the design strategies that have led to the identification of melatonin receptor ligands, guided by in silico approaches and molecular modeling. Initial ligand-based design, mainly relying on pharmacophore modeling and 3D-QSAR, has been flanked by structure-based virtual screening, given the recent availability of MT₁ and MT₂ crystal structures. Receptor ligands with different activity profiles, agonist/antagonist and subtype-selective compounds, are available.

EXPERT OPINION

An insight on the pharmacological characterization and therapeutic perspectives for relevant ligands is provided. In silico drug discovery has been instrumental in the design of novel ligands targeting melatonin receptors. Ligand-based approaches has led to the construction of a solid framework defining structure-activity relationships to obtain compounds with a tailored pharmacological profile. Structure-based techniques could integrate previous knowledge and provide compounds with novel chemotypes and pharmacological activity as drug candidates for disease conditions in which melatonin receptor ligands are currently being investigated, including cancer and pain.

Structural Refinement of 2,4-Thiazolidinedione Derivatives as New Anticancer Agents Able to Modulate the BAG3 Protein

The multidomain BAG3 protein is a member of the BAG (Bcl-2-associated athanogene) family of co-chaperones, involved in a wide range of protein-protein interactions crucial for many key cellular pathways, including autophagy, cytoskeletal dynamics, and apoptosis. Basal expression of BAG3 is elevated in several tumor cell lines, where it promotes cell survival signaling and apoptosis resistance through the interaction with many protein partners. In addition, its role as a key player of several hallmarks of cancer, such as metastasis, angiogenesis, autophagy activation, and apoptosis inhibition, has been established. Due to its involvement in malignant transformation, BAG3 has emerged as a potential and effective biological target to control multiple cancer-related signaling pathways. Recently, by using a multidisciplinary approach we reported the first synthetic BAG3 modulator interfering with its BAG domain (BD), based on a 2,4-thiazolidinedione scaffold and endowed with significant anti-proliferative activity. Here, a further in silico-driven selection of a 2,4-thiazolidinedione-based compound was performed. Thanks to a straightforward synthesis, relevant binding affinity for the BAG3BD domain, and attractive biological activities, this novel generation of compounds is of great interest for the development of further BAG3 binders, as well as for the elucidation of the biological roles of this protein in tumors. Specifically, we found compound 6 as a new BAG3 modulator with a relevant antiproliferative effect on two different cancer cell lines (IC₅₀: A375 = 19.36 μM; HeLa = 18.67 μM).

Efficiency of Homology Modeling Assisted Molecular Docking in G-protein Coupled Receptors

BACKGROUND

Molecular docking is in regular practice to assess ligand affinity on a target protein crystal structure. In the absence of protein crystal structure, the homology modeling or comparative modeling is the best alternative to elucidate the relationship details between a ligand and protein at the molecular level. The development of accurate homology modeling (HM) and its integration with molecular docking (MD) is essential for successful, rational drug discovery.

OBJECTIVE

The G-protein coupled receptors (GPCRs) are attractive therapeutic targets due to their immense role in human pharmacology. The GPCRs are membrane-bound proteins with the complex constitution, and the understanding of their activation and inactivation mechanisms is quite challenging. Over the past decade, there has been a rapid expansion in the number of solved G-protein-coupled receptor (GPCR) crystal structures; however, the majority of the GPCR structures remain unsolved. In this context, HM guided MD has been widely used for structure-based drug design (SBDD) of GPCRs.

METHODS

The focus of this review is on the recent (i) developments on HM supported GPCR drug discovery in the absence of GPCR crystal structures and (ii) application of HM in understanding the ligand interactions at the binding site, virtual screening, determining receptor subtype selectivity and receptor behaviour in comparison with GPCR crystal structures.

RESULTS

The HM in GPCRs has been extremely challenging due to the scarcity in template structures. In such a scenario, it is difficult to get accurate HM that can facilitate understanding of the ligand-receptor interactions. This problem has been alleviated to some extent by developing refined HM based on incorporating active /inactive ligand information and inducing protein flexibility. In some cases, HM proteins were found to outscore crystal structures.

CONCLUSION

The developments in HM have been highly operative to gain insights about the ligand interaction at the binding site and receptor functioning at the molecular level. Thus, HM guided molecular docking may be useful for rational drug discovery for the GPCRs mediated diseases.

Can molecular dynamics simulations improve the structural accuracy and virtual screening performance of GPCR models?

The determination of G protein-coupled receptor (GPCR) structures at atomic resolution has improved understanding of cellular signaling and will accelerate the development of new drug candidates. However, experimental structures still remain unavailable for a majority of the GPCR family. GPCR structures and their interactions with ligands can also be modelled computationally, but such predictions have limited accuracy. In this work, we explored if molecular dynamics (MD) simulations could be used to refine the accuracy of in silico models of receptor-ligand complexes that were submitted to a community-wide assessment of GPCR structure prediction (GPCR Dock). Two simulation protocols were used to refine 30 models of the D3 dopamine receptor (D3R) in complex with an antagonist. Close to 60 μs of simulation time was generated and the resulting MD refined models were compared to a D3R crystal structure. In the MD simulations, the receptor models generally drifted further away from the crystal structure conformation. However, MD refinement was able to improve the accuracy of the ligand binding mode. The best refinement protocol improved agreement with the experimentally observed ligand binding mode for a majority of the models. Receptor structures with improved virtual screening performance, which was assessed by molecular docking of ligands and decoys, could also be identified among the MD refined models. Application of weak restraints to the transmembrane helixes in the MD simulations further improved predictions of the ligand binding mode and second extracellular loop. These results provide guidelines for application of MD refinement in prediction of GPCR-ligand complexes and directions for further method development.

Pharmacological Actions of Carbamate Insecticides at Mammalian Melatonin Receptors

Integrated in silico chemical clustering and melatonin receptor molecular modeling combined with in vitro 2-[125I]-iodomelatonin competition binding were used to identify carbamate insecticides with affinity for human melatonin receptor 1 (hMT1) and human melatonin receptor 2 (hMT2). Saturation and kinetic binding studies with 2-[125I]-iodomelatonin revealed lead carbamates (carbaryl, fenobucarb, bendiocarb, carbofuran) to be orthosteric ligands with antagonist apparent efficacy at hMT1 and agonist apparent efficacy at hMT2 Furthermore, using quantitative receptor autoradiography in coronal brain slices from C3H/HeN mice, carbaryl, fenobucarb, and bendiocarb competed for 2-[125I]-iodomelatonin binding in the suprachiasmatic nucleus (SCN), paraventricular nucleus of the thalamus (PVT), and pars tuberalis (PT) with affinities similar to those determined for the hMT1 receptor. Carbaryl (10 mg/kg i.p.) administered in vivo also competed ex vivo for 2-[125I]-iodomelatonin binding to the SCN, PVT, and PT, demonstrating the ability to reach brain melatonin receptors in C3H/HeN mice. Furthermore, the same dose of carbaryl given to C3H/HeN mice in constant dark for three consecutive days at subjective dusk (circadian time 10) phase-advanced circadian activity rhythms (mean = 0.91 hours) similar to melatonin (mean = 1.12 hours) when compared with vehicle (mean = 0.04 hours). Carbaryl-mediated phase shift of overt circadian activity rhythm onset is likely mediated via interactions with SCN melatonin receptors. Based on the pharmacological actions of carbaryl and other carbamate insecticides at melatonin receptors, exposure may modulate time-of-day information conveyed to the master biologic clock relevant to adverse health outcomes. SIGNIFICANCE STATEMENT: In silico chemical clustering and molecular modeling in conjunction with in vitro bioassays identified several carbamate insecticides (i.e., carbaryl, carbofuran, fenobucarb, bendiocarb) as pharmacologically active orthosteric melatonin receptor 1 and 2 ligands. This work further demonstrated that carbaryl competes for melatonin receptor binding in the master biological clock (suprachiasmatic nucleus) and phase-advances overt circadian activity rhythms in C3H/HeN mice, supporting the relevance of circadian effects when interpreting toxicological findings related to carbamate insecticide exposure.

A Virtual Screening Platform Identifies Chloroethylagelastatin A as a Potential Ribosomal Inhibitor

Chloroethylagelastatin A (CEAA) is an analogue of agelastatin A (AA), a natural alkaloid derived from a marine sponge. It is under development for therapeutic use against brain tumors as it has excellent central nervous system (CNS) penetration and pre-clinical therapeutic activity against brain tumors. Recently, AA was shown to inhibit protein synthesis by binding to the ribosomal A-site. In this study, we developed a novel virtual screening platform to perform a comprehensive screening of various AA analogues showing that AA analogues with proven therapeutic activity including CEAA have significant ribosomal binding capacity whereas therapeutically inactive analogues show poor ribosomal binding and revealing structural fingerprint features essential for drug-ribosome interactions. In particular, CEAA was found to have greater ribosomal binding capacity than AA. Biological tests showed that CEAA binds the ribosome and contributes to protein synthesis inhibition. Our findings suggest that CEAA may possess ribosomal inhibitor activity and that our virtual screening platform may be a useful tool in discovery and development of novel ribosomal inhibitors.

Melatonin receptors in Atlantic salmon stimulate cAMP levels in heterologous cell lines and show season-dependent daily variations in pituitary expression levels

The hormone melatonin connects environmental cues, such as photoperiod and temperature, with a number of physiological and behavioural processes, including seasonal reproduction, through binding to their cognate receptors. This study reports the structural, functional and physiological characterization of five high-affinity melatonin receptors (Mtnr1aaα, Mtnr1aaβ, Mtnr1ab, Mtnr1al, Mtnr1b) in Atlantic salmon. Phylogenetic analysis clustered salmon melatonin receptors into three monophyletic groups, Mtnr1A, Mtnr1Al and Mtnr1B, but no functional representative of the Mtnr1C group. Contrary to previous studies in vertebrates, pharmacological characterization of four receptors in COS-7, CHO and SH-SY5Y cell lines (Mtnr1Aaα, Mtnr1Aaβ, Mtnr1Ab, Mtnr1B) showed induction of intracellular cAMP levels following 2-iodomelatonin or melatonin exposure. No consistent response was measured after N-acetyl-serotonin or serotonin exposure. Melatonin receptor genes were expressed at all levels of the hypothalamo-pituitary-gonad axis, with three genes (mtnr1aaβ, mtnr1ab and mtnr1b) detected in the pituitary. Pituitary receptors displayed daily fluctuations in mRNA levels during spring, prior to the onset of gonadal maturation, but not in autumn, strongly implying a direct involvement of melatonin in seasonal processes regulated by the pituitary. To the best of our knowledge, this is the first report of cAMP induction mediated via melatonin receptors in a teleost species.

Discovery of small molecule agonists targeting neuropeptide Y4 receptor using homology modeling and virtual screening

Neuropeptide Y4 receptor has the most significant effect on body weight and fat mass in its physiological functions, and the activation of Y4 receptor has explicit role on losing weight. The Y4 receptor has been successfully applied in the development of anti-obesity agent, thus representing a potential therapeutic target for obesity treatment. Here, we reported the first discovery of small molecule agonists targeting Y4 receptor: three Y4 receptor models with active and inactive conformations were built, each model was submitted following structure-based virtual screening, and finally six hits were identified as Y4 receptor agonists. These results confirm the reliability of the constructed Y4 receptor models and the proposed computational strategy for investigating novel Y4 receptor agonists. These new small molecule Y4 receptor agonists will contribute to the further development of Y4 agonists as potential therapeutics and functional probes.

Polybenzyls from Gastrodia elata, their agonistic effects on melatonin receptors and structure-activity relationships

Gastrodia elata is a famous traditional Chinese herb with medicinal and edible application. In this study, nine polybenzyls (1-9), including six new ones (2-5, 7 and 9), were isolated from the EtOAc extract of G. elata. Five compounds 1, 3, 4, 6 and 8 were found to activate melatonin receptors. Especially, compound 1 showed agonistic effects on MT₁ and MT₂ receptors with EC₅₀ values of 237 and 244 μM. For better understanding their structure-activity relationships (SARs), ten polybenzyl analogs were further synthesized and assayed for their activities on melatonin receptors. Preliminary SARs study suggested that two para-hydroxy groups were the key pharmacophore for maintaining activity. Molecular docking simulations verified that compound 1 could strongly interact with MT₂ receptor by bonding to Phe 118, Gly 121, His 208, Try 294 and Ala 297 residues.

Discovery of 3-hydroxy-3-pyrrolin-2-one-based mPGES-1 inhibitors using a multi-step virtual screening protocol

Targeting microsomal prostaglandin E₂ synthase-1 (mPGES-1) represents an efficient strategy for the development of novel drugs against inflammation and cancer with potentially reduced side effects. With this aim, a virtual screening was performed on a large library of commercially available molecules using the X-ray structure of mPGES-1 co-complexed with a potent inhibitor. Combining fast ligand-based shape alignment, molecular docking experiments, and qualitative analysis of the binding poses, a small set of molecules was selected for the subsequent steps of validation of the biological activity. Compounds 2 and 3, bearing the 3-hydroxy-3-pyrrolin-2-one nucleus, showed mPGES-1-inhibitory activity in the low micromolar range. These data highlighted the applicability of the reported virtual screening protocol for the selection of new mPGES-1 inhibitors as promising anti-inflammatory/anti-cancer drugs.

Docking studies for melatonin receptors

INTRODUCTION

Melatonin is a neurohormone that controls many relevant physiological processes beyond the control of circadian rhythms. Melatonin's actions are carried out by two main types of melatonin receptors; MT₁ and MT₂. These receptors are important, and not just because of the biological actions of its natural agonist; but also, because melatonin analogues can improve or antagonize their biological effect. Area covered: The following article describes the importance of melatonin as a biologically relevant molecule. It also defines the receptors for this substance, as well as the second messengers coupled to these receptors. Lastly, the article describes the amino acid residues involved in the docking process in both MT₁ and MT₂ melatonin receptors. Expert opinion: The biological actions of melatonin and their interpretations are becoming more relevant and therefore require the development of new pharmacological tools. Understanding the second messenger mechanisms involved in melatonin actions, as well as the characteristics of the docking of this molecule to MT₁ and MT₂ melatonin receptors, will permit the development of more selective agonists and antagonists which will help us to better understand this molecule as well to develop new therapeutic compounds.

Carbamate Insecticides Target Human Melatonin Receptors

Carbaryl (1-naphthyl methylcarbamate) and carbofuran (2,3-dihydro-2,2-dimethyl-7-benzofuranyl methylcarbamate) are among the most toxic insecticides, implicated in a variety of diseases including diabetes and cancer among others. Using an integrated pharmacoinformatics based screening approach, we have identified these insecticides to be structural mimics of the neurohormone melatonin and were able to bind to the putative melatonin binding sites in MT₁ and MT₂ melatonin receptors in silico. Carbaryl and carbofuran then were tested for competition with 2-[¹²⁵I]-iodomelatonin (300 pM) binding to hMT₁ or hMT₂ receptors stably expressed in CHO cells. Carbaryl and carbofuran showed higher affinity for competition with 2-[¹²⁵I]-iodomelatonin binding to the hMT₂ compared to the hMT₁ melatonin receptor (33 and 35-fold difference, respectively) as predicted by the molecular modeling. In the presence of GTP (100 μM), which decouples the G-protein linked receptors to modulate signaling, the apparent efficacy of carbaryl and carbofuran for 2-[¹²⁵I]-iodomelatonin binding for the hMT₁ melatonin receptor was not affected but significantly decreased for the hMT₂ melatonin receptor compatible with receptor antagonist/inverse agonist and agonist efficacy, respectively. Altogether, our data points to a potentially new mechanism through which carbamate insecticides carbaryl and carbofuran could impact human health by altering the homeostatic balance of key regulatory processes by directly binding to melatonin receptors.

Update on melatonin receptors: IUPHAR Review 20

Melatonin receptors are seven transmembrane-spanning proteins belonging to the GPCR superfamily. In mammals, two melatonin receptor subtypes exist - MT1 and MT2 - encoded by the MTNR1A and MTNR1B genes respectively. The current review provides an update on melatonin receptors by the corresponding subcommittee of the International Union of Basic and Clinical Pharmacology. We will highlight recent developments of melatonin receptor ligands, including radioligands, and give an update on the latest phenotyping results of melatonin receptor knockout mice. The current status and perspectives of the structure of melatonin receptor will be summarized. The physiological importance of melatonin receptor dimers and biologically important and type 2 diabetes-associated genetic variants of melatonin receptors will be discussed. The role of melatonin receptors in physiology and disease will be further exemplified by their functions in the immune system and the CNS. Finally, antioxidant and free radical scavenger properties of melatonin and its relation to melatonin receptors will be critically addressed.

Comparative Analysis of Virtual Screening Approaches in the Search for Novel EphA2 Receptor Antagonists

The EphA2 receptor and its ephrin-A1 ligand form a key cell communication system, which has been found overexpressed in many cancer types and involved in tumor growth. Recent medicinal chemistry efforts have identified bile acid derivatives as low micromolar binders of the EphA2 receptor. However, these compounds suffer from poor physicochemical properties, hampering their use in vivo. The identification of compounds able to disrupt the EphA2-ephrin-A1 complex lacking the bile acid scaffold may lead to new pharmacological tools suitable for in vivo studies. To identify the most promising virtual screening (VS) protocol aimed at finding novel EphA2 antagonists, we investigated the ability of both ligand-based and structure-based approaches to retrieve known EphA2 antagonists from libraries of decoys with similar molecular properties. While ligand-based VSs were conducted using UniPR129 and ephrin-A1 ligand as reference structures, structure-based VSs were performed with Glide, using the X-ray structure of the EphA2 receptor/ephrin-A1 complex. A comparison of enrichment factors showed that ligand-based approaches outperformed the structure-based ones, suggesting ligand-based methods using the G-H loop of ephrin-A1 ligand as template as the most promising protocols to search for novel EphA2 antagonists.

Novel N-Acetyl Bioisosteres of Melatonin: Melatonergic Receptor Pharmacology, Physicochemical Studies, and Phenotypic Assessment of Their Neurogenic Potential

Herein we present a new family of melatonin-based compounds, in which the acetamido group of melatonin has been bioisosterically replaced by a series of reversed amides and azoles, such as oxazole, 1,2,4-oxadiazole, and 1,3,4-oxadiazole, as well as other related five-membered heterocycles, namely, 1,3,4-oxadiazol(thio)ones, 1,3,4-triazol(thio)ones, and an 1,3,4-thiadiazole. New compounds were fully characterized at melatonin receptors (MT1R and MT2R), and results were rationalized by superimposition studies of their structures to the bioactive conformation of melatonin. We also found that several of these melatonin-based compounds promoted differentiation of rat neural stem cells to a neuronal phenotype in vitro, in some cases to a higher extent than melatonin. This unique profile constitutes the starting point for further pharmacological studies to assess the mechanistic pathways and the relevance of neurogenesis induced by melatonin-related structures.

Docking and Virtual Screening Strategies for GPCR Drug Discovery

Progress in structure determination of G protein-coupled receptors (GPCRs) has made it possible to apply structure-based drug design (SBDD) methods to this pharmaceutically important target class. The quality of GPCR structures available for SBDD projects fall on a spectrum ranging from high resolution crystal structures (<2 Å), where all water molecules in the binding pocket are resolved, to lower resolution (>3 Å) where some protein residues are not resolved, and finally to homology models that are built using distantly related templates. Each GPCR project involves a distinct set of opportunities and challenges, and requires different approaches to model the interaction between the receptor and the ligands. In this review we will discuss docking and virtual screening to GPCRs, and highlight several refinement and post-processing steps that can be used to improve the accuracy of these calculations. Several examples are discussed that illustrate specific steps that can be taken to improve upon the docking and virtual screening accuracy. While GPCRs are a unique target class, many of the methods and strategies outlined in this review are general and therefore applicable to other protein families.

Expanding the horizons of G protein-coupled receptor structure-based ligand discovery and optimization using homology models

We show the key role of structural homology models in GPCR structure-based lead discovery and optimization, highlighting methodological aspects, recent progress and future directions.

Structure and putative signaling mechanism of Protease activated receptor 2 (PAR2) - a promising target for breast cancer

Experimental evidences have observed enhanced expression of protease activated receptor 2 (PAR2) in breast cancer consistently. However, it is not yet recognized as an important therapeutic target for breast cancer as the primary molecular mechanisms of its activation are not yet well-defined. Nevertheless, recent reports on the mechanism of GPCR activation and signaling have given new insights to GPCR functioning. In the light of these details, we attempted to understand PAR2 structure & function using molecular modeling techniques. In this work, we generated averaged representative stable models of PAR2, using protease activated receptor 1 (PAR1) as a template and selected conformation based on their binding affinity with PAR2 specific agonist, GB110. Further, the selected model was used for studying the binding affinity of putative ligands. The selected ligands were based on a recent publication on phylogenetic analysis of Class A rhodopsin family of GPCRs. This study reports putative ligands, their interacting residues, binding affinity and molecular dynamics simulation studies on PAR2-ligand complexes. The results reported from this study would be useful for researchers and academicians to investigate PAR2 function as its physiological role is still hypothetical. Further, this information may provide a novel therapeutic scheme to manage breast cancer.

An exhaustive yet simple virtual screening campaign against Sortase A from multiple drug resistant Staphylococcus aureus

Methicillin resistant Staphylococcus aureus (MRSA) is one of the challenging bacterial pathogen due to its acquired resistance to the β lactam antibiotics. The Sortase A is an enzyme of Gram-positive bacteria including S. aureus to anchor surface proteins to the cell wall. Sortase A is well studied enzyme and considered as the drug target against MRSA. Sortase A plays active role in anchoring the virulence proteins on the cell wall of the Gram-positive bacteria. The inhibition of Sortase A activity results in the separation of S. aureus from the host cells and ultimately alleviation of the infection. Here, we adapted a structure-based virtual screening protocol which helped in identification of novel potential inhibitors of Sortase A. The protocol involved the docking of a chemical library of druglike compounds with the Sortase A binding site represented by multiple crystal structures. The compounds were ranked by multiple scoring functions and shortlisted for future experimental screening. The method resulted in shortlisting of three compounds as potential novel inhibitors of Sortase A out of a large chemical library. The high rankings of shortlisted compounds estimated by multiple scoring functions showed their binding potential with Sortase A. The results are proved to be a simple yet efficient choice of structure-based virtual screening. The identified compounds are druglike and show high rankings among all set protocols of the virtual screening. We hope that the study would eventually help to expedite the discovery of novel drug candidates against MRSA.

Homology modeling of human muscarinic acetylcholine receptors

We have developed homology models of the acetylcholine muscarinic receptors M₁R-M₅R, based on the β₂-adrenergic receptor crystal as the template. This is the first report of homology modeling of all five subtypes of acetylcholine muscarinic receptors with binding sites optimized for ligand binding. The models were evaluated for their ability to discriminate between muscarinic antagonists and decoy compounds using virtual screening using enrichment factors, area under the ROC curve (AUC), and an early enrichment measure, LogAUC. The models produce rational binding modes of docked ligands as well as good enrichment capacity when tested against property-matched decoy libraries, which demonstrates their unbiased predictive ability. To test the relative effects of homology model template selection and the binding site optimization procedure, we generated and evaluated a naïve M₂R model, using the M₃R crystal structure as a template. Our results confirm previous findings that binding site optimization using ligand(s) active at a particular receptor, i.e. including functional knowledge into the model building process, has a more pronounced effect on model quality than target-template sequence similarity. The optimized M₁R-M₅R homology models are made available as part of the Supporting Information to allow researchers to use these structures, compare them to their own results, and thus advance the development of better modeling approaches.

A molecular and chemical perspective in defining melatonin receptor subtype selectivity

Melatonin is primarily synthesized and secreted by the pineal gland during darkness in a normal diurnal cycle. In addition to its intrinsic antioxidant property, the neurohormone has renowned regulatory roles in the control of circadian rhythm and exerts its physiological actions primarily by interacting with the G protein-coupled MT₁ and MT₂ transmembrane receptors. The two melatonin receptor subtypes display identical ligand binding characteristics and mediate a myriad of signaling pathways, including adenylyl cyclase inhibition, phospholipase C stimulation and the regulation of other effector molecules. Both MT₁ and MT₂ receptors are widely expressed in the central nervous system as well as many peripheral tissues, but each receptor subtype can be linked to specific functional responses at the target tissue. Given the broad therapeutic implications of melatonin receptors in chronobiology, immunomodulation, endocrine regulation, reproductive functions and cancer development, drug discovery and development programs have been directed at identifying chemical molecules that bind to the two melatonin receptor subtypes. However, all of the melatoninergics in the market act on both subtypes of melatonin receptors without significant selectivity. To facilitate the design and development of novel therapeutic agents, it is necessary to understand the intrinsic differences between MT₁ and MT₂ that determine ligand binding, functional efficacy, and signaling specificity. This review summarizes our current knowledge in differentiating MT₁ and MT₂ receptors and their signaling capacities. The use of homology modeling in the mapping of the ligand-binding pocket will be described. Identification of conserved and distinct residues will be tremendously useful in the design of highly selective ligands.