Molecular modeling of the GABAC receptor ligand-binding domain

We have constructed a molecular model of the ligand-binding domain of the GABA(C) receptor, which is a member of the Cys-loop ligand-gated ion channel family. The extracellular domains of these receptors share similar sequence homology (20%) with Limnaea acetylcholine-binding protein for which an X-ray crystal structure is available. We used this structure as a template for homology modeling of the GABA(C) receptor extracellular domain using FUGUE and MODELLER software. FlexX was then used to dock GABA into the receptor ligand-binding site, resulting in three alternative energetically favorable orientations. Residues located no more than 5 A from the docked GABA were identified for each model; of these, three were found to be common to all models with 14 others present only in certain models. Using data from experimental studies, we propose that the most likely orientation of GABA is with its amine close to Y198, and its carboxylate close to R104. These studies have therefore provided a model of the ligand-binding domain, which will be useful for both GABA(C) and GABA(A) receptor studies, and have also yielded an experimentally testable hypothesis of the location of GABA in the binding pocket. [Figure: see text].     

Prediction of 5-HT3 receptor agonist-binding residues using homology modeling

5-HT(3) receptors demonstrate significant structural and functional homology to other members of the Cys-loop ligand-gated ion channel superfamily. The extracellular domains of these receptors share similar sequence homology (approximately 20%) with Limnaea acetylcholine binding protein, for which an x-ray crystal structure is available. We used this structure as a template for computer-based homology modeling of the 5-HT(3) receptor extracellular domain. AutoDock software was used to dock 5-HT into the putative 5-HT(3) receptor ligand-binding site, resulting in seven alternative energetically favorable models. Residues located no more than 5 A from the docked 5-HT were identified for each model; of these, 12 were found to be common to all seven models with five others present in only certain models. Some docking models reflected the cation-pi interaction previously demonstrated for W183, and data from these and other studies were used to define our preferred models.     

Binding mode of ecdysone agonists to the receptor: comparative modeling and docking studies

Three-dimensional structure models of the ligand-binding domain of the ecdysone receptor of Heliothis virescens were built by the homology modeling technique from the crystal structures of nuclear receptors. Two models were created based both on known ligand-binding domain structures of the receptors with the highest sequence identity to the ecdysone receptor, and on those of steroid hormone receptors. The latter model, which was found to have better stereochemical quality and be in good agreement with the binding of the steroidal framework of the endogenous agonist 20-hydroxyecdysone, was used for docking studies. The docking of 20-hydroxyecdysone to the receptor model revealed that the ligand molecule can interact with the receptor in a similar manner to other steroid hormone-receptor complexes. The docking of a dibenzoylhydrazine agonist, chromafenozide, was performed based on the correspondences between the molecule and 20-dydroxyecdysone expected by molecular comparison. The interactions of the ligands with the receptor in the complexes modeled were investigated and found to be consistent with known structure-activity relationships.     

Towards the understanding of tetrahydroquinolines action in Aedes aegypti: larvicide or adulticide?

Aedes aegypti is an important vector of arboviruses such as dengue, yellow fever, chikungunya and Zika. Among the various types of insecticides used to combat this vector, the insect growth regulators have been developed and recommended for control of their larvae. In this work compounds with proven regulatory action, tetrahydroquinolines will be studied. These regulators act on the hormones responsible for the insect development. Ecdysone, one of the main hormones involved in this process has a specific receptor (EcR), where tetrahydroquinolines derivatives can bind, disrupting the normal action of this hormone, because they have structure similar to hormone 20-hydroxyecdysone (20E). In addition, studies show that this class of compounds interacts strongly in the potassium channel activated by calcium (BK channel). Thus, the goal is to study the action of compounds (tetrahydroquinolines) as insecticides and evaluate their larvicidal action (action on the ecdysone receptor) or adulticide (action on the BK channel) through homology modelling techniques, molecular docking and molecular dynamics simulations and aiming to propose a compound that presents both actions (larvicide / adulticide).     

Homology-modeled ligand-binding domains of medaka estrogen receptors and androgen receptors: A model system for the study of reproduction

Estrogen and androgen and their receptors play critical roles in physiological processes such as sexual differentiation and development. Using the available structural models for the human estrogen receptors alpha and beta and androgen receptor as templates, we designed in silico agonist and antagonist models of medaka estrogen receptor (meER) alpha, beta-1, and beta-2, and androgen receptor (meAR) alpha and beta. Using these models, we studied (1) the structural relationship between the ligand-binding domains (LBDs) of ERs and ARs of human and medaka, and (2) whether medaka ER and AR can be potential models for studying the ligand-binding activities of various agonists and antagonists of these receptors by docking analysis. A high level of conservation was observed between the sequences of the ligand-binding domains of meERα and huERα, meERβ1 and huERβ, meERβ2, and huERβ with 62.8%, 66.4%, and 65.1% identity, respectively. The sequence conservation between meARα and huAR, meARβ, and huAR was found with 70.1% and 61.0% of identity, respectively. Thirty-three selected endocrine disrupting chemicals (EDCs), including both agonists and antagonists, were docked into the LBD of ER and AR, and the corresponding docking score for medaka models and human templates were calculated. In order to confirm the conservation of the overall geometry and the binding pocket, the backbone root mean square deviation (RMSD) for Cα atoms was derived from the structure superposition of all 10 medaka homology models to the six human templates. Our results suggested conformational conservation between the ERs and ARs of medaka and human, Thus, medaka could be highly useful as a model system for studies involving estrogen and androgen interaction with their receptors.     

Protein‐based virtual screening of chemical databases. II. Are homology models of g‐protein coupled receptors suitable targets?

The aim of the current study is to investigate whether homology models of G-Protein-Coupled Receptors (GPCRs) that are based on bovine rhodopsin are reliable enough to be used for virtual screening of chemical databases. Starting from the recently described 2.8 A-resolution X-ray structure of bovine rhodopsin, homology models of an "antagonist-bound" form of three human GPCRs (dopamine D3 receptor, muscarinic M1 receptor, vasopressin V1a receptor) were constructed. The homology models were used to screen three-dimensional databases using three different docking programs (Dock, FlexX, Gold) in combination with seven scoring functions (ChemScore, Dock, FlexX, Fresno, Gold, Pmf, Score). Rhodopsin-based homology models turned out to be suitable, indeed, for virtual screening since known antagonists seeded in the test databases could be distinguished from randomly chosen molecules. However, such models are not accurate enough for retrieving known agonists. To generate receptor models better suited for agonist screening, we developed a new knowledge- and pharmacophore-based modeling procedure that might partly simulate the conformational changes occurring in the active site during receptor activation. Receptor coordinates generated by this new procedure are now suitable for agonist screening. We thus propose two alternative strategies for the virtual screening of GPCR ligands, relying on a different set of receptor coordinates (antagonist-bound and agonist-bound states).     

Sequence annotation of nuclear receptor ligand-binding domains by automated homology modeling

The quality of three-dimensional homology models derived from protein sequences provides an independent measure of the suitability of a protein sequence for a certain fold. We have used automated homology modeling and model assessment tools to identify putative nuclear hormone receptor ligand-binding domains in the genome of Caenorhabditis elegans. Our results indicate that the availability of multiple crystal structures is crucial to obtaining useful models in this receptor family. The majority of annotated mammalian nuclear hormone receptors could be assigned to a ligand-binding domain fold by using the best model derived from any of four template structures. This strategy also assigned the ligand-binding domain fold to a number of C.elegans. sequences without prior annotation. Interestingly, the retinoic acid receptor crystal structure contributed most to the number of sequences that could be assigned to a ligand-binding domain fold. Several causes for this can be suggested, including the high quality of this protein structure in terms of our assessment tools, similarity between the biological function or ligand of this receptor and the modeled genes and gene duplication in C.elegans.     

Multivalency - a way to enhance binding avidities and bioactivity - preliminary applications to EPO.

Multivalency has advantages over monovalency for binding interactions and even for activity. In particular, avidity is higher since the off-rate of a multivalent species is much slower than that of a monomer. This is particularly profitable for ligand-binding receptors that require dimerization for activity, such as the receptor of erythropoietin (EPOR). Peptides that mimic the action of erythropoietin (EPO) have been described with no sequence similarity with the human hormone: erythropoietin mimetic peptide (EMP) and EPO receptor peptide (ERP). These two peptides have similar activity but interact through different sites on the EPOR. Here, we describe the construction of several new synthetic homo- and hetero-dimers based on EMP-ERP sequences. To link the monomeric molecules together, several monodisperse polyamide linkers of different lengths were synthesized with dialdehyde functionalities. Chemoselective oxime chemistry was used to obtain homogeneous constructs. Certain chemical incompatibilities were dealt with via a protection approach. The oximes are stable under normal conditions and so lend themselves to biological testing.     

Using nondenaturing mass spectrometry to detect fortuitous ligands in orphan nuclear receptors

Nondenaturing electrospray mass spectrometry (ESI-MS) has been used to reveal the presence of potential ligands in the ligand-binding domain (LBD) of orphan nuclear receptors. This new approach, based on supramolecular mass spectrometry, allowed the detection and identification of fortuitous ligands for the retinoic acid-related orphan receptor beta (RORbeta) and the ultraspiracle protein (USP). These fortuitous ligands were specifically captured from the host cell with the proper stoichiometry. After organic extraction, these molecules have been characterized by classic analytical methods and identified as stearic acid for RORbeta and a phosphatidylethanolamine (PE) for USP, as confirmed by crystallography. These molecules act as "fillers" and may not be the physiological ligands, but they prove to be essential to stabilize the active conformation of the LBD, enabling its crystallization. The resulting crystal structures provide a detailed picture of the ligand-binding pocket, allowing the design of highly specific synthetic ligands that can be used to characterize the function of orphan nuclear receptors. An additional advantage of this new method is that it is not based on a functional test and that it can detect low-affinity ligands.     

T2384, a novel antidiabetic agent with unique peroxisome proliferator-activated receptor gamma binding properties

The nuclear hormone receptor peroxisome proliferator-activated receptor gamma (PPARgamma) plays central roles in adipogenesis and glucose homeostasis and is the molecular target for the thiazolidinedione (TZD) class of antidiabetic drugs. Activation of PPARgamma by TZDs improves insulin sensitivity; however, this is accompanied by the induction of several undesirable side effects. We have identified a novel synthetic PPARgamma ligand, T2384, to explore the biological activities associated with occupying different regions of the receptor ligand-binding pocket. X-ray crystallography studies revealed that T2384 can adopt two distinct binding modes, which we have termed "U" and "S", interacting with the ligand-binding pocket of PPARgamma primarily via hydrophobic contacts that are distinct from full agonists. The different binding modes occupied by T2384 induced distinct patterns of coregulatory protein interaction with PPARgamma in vitro and displayed unique receptor function in cell-based activity assays. We speculate that these unique biochemical and cellular activities may be responsible for the novel in vivo profile observed in animals treated systemically with T2384. When administered to diabetic KKAy mice, T2384 rapidly improved insulin sensitivity in the absence of weight gain, hemodilution, and anemia characteristics of treatment with rosiglitazone (a TZD). Moreover, upon coadministration with rosiglitazone, T2384 was able to antagonize the side effects induced by rosiglitazone treatment alone while retaining robust effects on glucose disposal. These results are consistent with the hypothesis that interactions between ligands and specific regions of the receptor ligand-binding pocket might selectively trigger a subset of receptor-mediated biological responses leading to the improvement of insulin sensitivity, without eliciting less desirable responses associated with full activation of the receptor. We suggest that T2384 may represent a prototype for a novel class of PPARgamma ligand and, furthermore, that molecules sharing some of these properties would be useful for treatment of type 2 diabetes.     

Three-dimensional structure and function study on the active region in the extracellular ligand-binding domain of human IL-6 receptor

In this study the three-dimensional (3-D) model of the ligand-binding domain (V106-P322) of human interleukin-6 receptor (hIL-6 R) was constructed by computer-guided homology modeling technique using the crystal structure of the ligand-binding domain (K52-L251) of human growth hormone receptor (hGHR) as templet. Furthermore, the active binding region of the 3-D model of hIL-6R with the ligand (hIL-6) was predicted. In light of the structural characteristics of the active region, a hydrophobic pocket shielded by two hydrophilic residues (E115 and E505) of the region was identified by a combination of molecular modelling and the site-directed or double-site mutation of the twelve crucial residues in the ligand-binding domain of hIL-6R (V106-P322). We observed and analyzed the effects of these mutants on the spatial conformation of the pocket-like region of hIL-6 R. The results indicated that any site-directed mutation of the five Cys residues (four conservative Cys residues: Cys121, Cys132, Cys165, Cys176; near membrane Cys residue: Cys193) or each double-site mutation of the five residues in WSEWS motif of hIL-6R (V106-P322) makes the corresponding spatial conformation of the pocket region block the linkage between hIL-6 R and hIL-6. However, the influence of the site-directed mutation of Cys211 and Cys277 individually on the conformation of the pocket region benefits the interaction between hIL-6R and hIL-6. Our study suggests that the predicted hydrophobic pocket in the 3-D model of hIL-6R (V106-P322) is the critical molecular basis for the binding of hIL-6R with its ligand, and the active pocket may be used as a target for designing small hIL-6R-inhibiting molecules in our further study.     

Crystal structure of the ligand-binding domain of the ultraspiracle protein USP, the ortholog of retinoid X receptors in insects

The major postembryonic developmental events happening in insect life, including molting and metamorphosis, are regulated and coordinated temporally by pulses of ecdysone. The biological activity of this steroid hormone is mediated by two nuclear receptors: the ecdysone receptor (EcR) and the Ultraspiracle protein (USP). The crystal structure of the ligand-binding domain from the lepidopteran Heliothis virescens USP reported here shows that the loop connecting helices H1 and H3 precludes the canonical agonist conformation. The key residues that stabilize this unique loop conformation are strictly conserved within the lepidopteran USP family. The presence of an unexpected bound ligand that drives an unusual antagonist conformation confirms the induced-fit mechanism accompanying the ligand binding. The ligand-binding pocket exhibits a retinoid X receptor-like anchoring part near a conserved arginine, which could interact with a USP ligand functional group. The structure of this receptor provides the template for designing inhibitors, which could be utilized as a novel type of environmentally safe insecticides.     

Structural basis for ligand recognition and activation of RAGE

The receptor for advanced glycation end products (RAGE) is a pattern recognition receptor involved in?inflammatory processes and is associated with diabetic complications, tumor outgrowth, and neurodegenerative disorders. RAGE induces cellular signaling events upon binding of a variety of ligands, such as glycated proteins, amyloid-β, HMGB1, and S100 proteins. The X-ray crystal structure of the VC1 ligand-binding region of the human RAGE ectodomain was determined at 1.85?? resolution. The VC1 ligand-binding surface was mapped onto the structure from titrations with S100B monitored by heteronuclear NMR spectroscopy. These NMR chemical shift perturbations were used as input for restrained docking calculations to generate a model for the VC1-S100B complex. Together, the arrangement of VC1 molecules in the crystal and complementary biochemical studies suggest a role for self-association in RAGE function. Our results enhance understanding of the functional outcomes of S100 protein binding to RAGE and provide insight into mechanistic models for how the receptor is activated.     

Structural determinants of imidacloprid-based nicotinic acetylcholine receptor inhibitors identified using 3D-QSAR, docking and molecular dynamics

Nicotinic acetylcholine receptor (nAChR) is a target for insect-selective neonicotinoid insecticides (NNs), exemplified by imidacloprid (IMI). In the present study, 78 IMI derivatives reported as inhibitors of Drosophila melanogaster nAChR (Dm-nAChR) and Musca domestica nAChR (Md-nAChR) were used for three-dimensional quantitative structure-activity relationship (3D-QSAR) studies. Two optimal models with good predictive power were obtained: Q(2) = 0.64, R(2)(pred) = 0.72 for Dm-nAChR, and Q(2) = 0.63, R(2)(pred) = 0.62 for Md-nAChR. In addition, homology modeling, molecular dynamic (MD) simulation, and molecular docking also showed that amino acids located within loops A, C, D and E play key roles in the interaction of Dm-/Md-nAChR with NNs. This is highly consistent with the results of graphical analysis of 3D-QSAR contour plots. Mutation analysis also implicates the Y/S mutation within loop B as being associated closely with NN resistance in Drosophila and Musca. The results obtained lead to a better understanding not only of interactions between these antagonists and Dm-/Md-nAChR, but also of the essential features that should be considered when designing novel inhibitors with desired activities.     

Small molecule inhibitors as probes for estrogen and androgen receptor action

Because activated estrogen (ER) and androgen (AR) receptors stimulate cell proliferation in breast and prostate cancer, inhibiting their actions represents a major therapeutic goal. Most efforts to modulate ER and AR activity have focused on inhibiting the synthesis of estrogens or androgens or on the identification of small molecules that act by competing with agonist hormones for binding in the ligand-binding pocket of the receptor. An alternative approach is to implement screens for small molecule inhibitors that target other sites in the pathway of steroid receptor action. Many of these second-site inhibitors directly target ER or AR; others have still unknown sites of action. Small molecule inhibitors that target second sites represent new leads with clinical potential; they serve as novel modulators of receptor action; and they can reveal new and as yet unidentified interactions and pathways that modulate ER and AR action.