On the Denaturation Mechanisms of the Ligand Binding Domain of Thyroid Hormone Receptors

The ligand binding domain (LBD) of nuclear hormone receptors adopts a very compact, mostly alpha-helical structure that binds specific ligands with very high affinity. We use circular dichroism spectroscopy and high-temperature molecular dynamics Simulations to investigate unfolding of the LBDs of thyroid hormone receptors (TRs). A molecular description of the denaturation mechanisms is obtained by molecular dynamics Simulations of the TR alpha and TR beta LBDs in the absence and in the presence of the natural ligand Triac. The Simulations Show that the thermal unfolding of the LBD starts with the loss of native contacts and secondary Structure elements, while the Structure remains essentially compact, resembling a molten globule state. This differs From most protein denaturation simulations reported to date and suggests that the folding mechanism may start with the hydrophobic collapse of the TR LBDs. Our results reveal that the stabilities of the LBDs of the TR alpha and TR beta Subtypes are affected to different degrees by the binding of the isoform selective ligand Triac and that ligand binding confers protection against thermal denaturation and unfolding in a subtype specific manner. Our Simulations indicate two mechanisms by which the ligand stabilizes the LBD: (1) by enhancing the interactions between H8 and H 11, and the interaction of the region between H I and the Omega-loop with the core of the LBD, and (2) by shielding the hydrophobic H6 from hydration.

Only subtle protein conformational adaptations are required for ligand binding to thyroid hormone receptors: Simulations using a novel multipoint steered molecular dynamics approach

Thyroid hormone receptors (TR) are hormone-dependent transcription regulators that play a major role in human health, development, and metabolic functions. The thyroid hormone resistance syndrome, diabetes, obesity, and some types of cancer are just a few examples of important diseases that are related to TR malfunctioning, particularly impaired hormone binding. Ligand binding to and dissociation from the receptor ultimately control gene transcription and, thus, detailed knowledge of binding and release mechanisms are fundamental for the comprehension of the receptor`s biological function and development of pharmaceuticals. In this work, we present the first computational study of ligand entry into the ligand binding domain (LBD) of a nuclear receptor. We report molecular dynamics simulations of ligand binding to TRs using a generalization of the steered molecular dynamics technique designed to perform single-molecule pulling simulations along arbitrarily nonlinear driving pathways. We show that only gentle protein movements and conformational adaptations are required for ligand entry into the LBDs and that the magnitude of the forces applied to assist ligand binding are of the order of the forces involved in ligand dissociation. Our simulations suggest an alternative view for the mechanisms ligand binding and dissociation of ligands from nuclear receptors in which ligands can simply diffuse through the protein surface to reach proper positioning within the binding pocket. The proposed picture indicates that the large-amplitude protein motions suggested by the apo- and holo-RXR alpha crystallographic structures are not required, reconciling conformational changes of LBDs required for ligand entry with other nuclear receptors apo-structures that resemble the ligand-bound LBDs.

Identification of a novel ligand binding motif in the transthyretin channel

The design of therapeutic compounds targeting transthyretin (TTR) is challenging due to the low specificity of interaction in the hormone binding site. Such feature is highlighted by the interactions of TTR with diclofenac, a compound with high affinity for TTR, in two dissimilar modes, as evidenced by crystal structure of the complex. We report here structural analysis of the interactions of TTR with two small molecules, 1-amino-5-naphthalene sulfonate (1,5-AmNS) and 1-anilino-8-naphthalene sulfonate (1,8-ANS). Crystal structure of TTR: 1,8-ANS complex reveals a peculiar interaction, through the stacking of the naphthalene ring between the side-chain of Lys15 and Leu17. The sulfonate moiety provides additional interaction with Lys15` and a water-mediated hydrogen bond with Thr119`. The uniqueness of this mode of ligand recognition is corroborated by the crystal structure of TTR in complex with the weak analogue 1,5-AmNS, the binding of which is driven mainly by hydrophobic partition and one electrostatic interaction between the sulfonate group and the Lys15. The ligand binding motif unraveled by 1,8-ANS may open new possibilities to treat TTR amyloid diseases by the elucidation of novel candidates for a more specific pharmacophoric pattern. (C) 2009 Published by Elsevier Ltd.

Ligand induced interaction of thyroid hormone receptor beta with its coregulators

Thyroid hormones exert most of their physiological effects through two thyroid hormone receptor (TR) subtypes, TR alpha and TR beta, which associate with many transcriptional coregulators to mediate activation or repression of target genes. The search for selective TR beta ligands has been stimulated by the finding that several pharmacological actions mediated by TR beta might be beneficial in medical conditions such as obesity, hypercholesterolemia and diabetes. Here, we present a new methodology which employs surface plasmon resonance to investigate the interactions between TR beta ligand binding domain (LBD) complexes and peptides derived from the nuclear receptor interaction motifs of two of its coregulators, SRC2 and DAX1. The effect of several TR beta ligands, including the TR beta selective agonist GC-I and the TR beta selective antagonist NH-3, were investigated. We also determined the kinetic rate constants for the interaction of TR beta-T3 with both coregulators, and accessed the thermodynamic parameters for the interaction with DAX1. Our findings Suggest that flexibility plays an important role in the interaction between the receptor and its coregulators. and point out important aspects of experimental design that should be addressed when using TR beta LBD and its agonists. Furthermore, the methodology described here may be useful for the identification of new TR beta ligands. (C) 2008 Elsevier Ltd. All rights reserved.

Ligand dissociation from estrogen receptor is mediated by receptor dimerization: Evidence from molecular dynamics simulations

Estrogen Receptor (ER) is an important target for pharmaceutical design. Like other ligand-dependent transcription factors, hormone binding regulates ER transcriptional activity. Nevertheless, the mechanisms by which ligands enter and leave ERs and other nuclear receptors remain poorly understood. Here, we report results of locally enhanced sampling molecular dynamics simulations to identify dissociation pathways of two ER ligands [the natural hormone 17 beta-estradiol (E-2) and the selective ER modulator raloxifene (RAL)] from the human ER alpha ligand-binding domain in monomeric and dimeric forms. E-2 dissociation occurs via three different pathways in ER monomers. One resembles the mousetrap mechanism (Path I), involving repositioning of helix 12 (H12), others involve the separation of H8 and H11 (Path II), and a variant of this pathway at the bottom of the ligand-binding domain (Path II`). RAL leaves the receptor through Path I and a Path I variant in which the ligand leaves the receptor through the loop region between H11 and H12 (Path I`). Remarkably, ER dimerization strongly suppresses Paths II and II` for E-2 dissociation and modifies RAL escape routes. We propose that differences in ligand release pathways detected in the simulations for ER monomers and dimers provide an explanation for previously observed effects of ER quaternary state on ligand dissociation rates and suggest that dimerization may play an important, and hitherto unexpected, role in regulation of ligand dissociation rates throughout the nuclear receptor family.

Structural modeling of high-affinity thyroid receptor-ligand complexes

Understanding the molecular basis of the binding modes of natural and synthetic ligands to nuclear receptors is fundamental to our comprehension of the activation mechanism of this important class of hormone regulated transcription factors and to the development of new ligands. Thyroid hormone receptors (TR) are particularly important targets for pharmaceuticals development because TRs are associated with the regulation of metabolic rates, body weight, and circulating levels of cholesterol and triglycerides in humans. While several high-affinity ligands are known, structural information is only partially available. In this work we obtain structural models of several TR-ligand complexes with unknown structure by docking high affinity ligands to the receptors` ligand binding domain with subsequent relaxation by molecular dynamics simulations. The binding modes of these ligands are discussed providing novel insights into the development of TR ligands. The experimental binding free energies are reasonably well-reproduced from the proposed models using a simple linear interaction energy free-energy calculation scheme.

Conformational differences between the wild type and V30M mutant transthyretin modulate its binding to genistein: Implications to tetramer stability and ligand-binding

Transthyretin (TTR) is a tetrameric beta-sheet-rich transporter protein directly involved in human amyloid diseases. It was recently found that the isoflavone genistein (GEN) potently inhibits TTR amyloid fibril formation (Green et al., 2005) and is therefore a promising candidate for TTR amyloidosis treatment. Here we used structural and biophysical approaches to characterize genistein binding to the wild type (TTRwt) and to its most frequent amyloidogenic variant, the V30M mutant. In a dose-dependent manner, genistein elicited considerable increases in both mutant and TTRwt stability as demonstrated by high hydrostatic pressure (HHP) and acid-mediated dissociation/denaturation assays. TTR:GEN crystal complexes and isothermal titration calorimetry (ITC) experiments showed that the binding mechanisms of genistein to the TTRwt and to V30M are different and are dependent on apoTTR structure conformations. Furthermore, we could also identify potential allosteric movements caused by genistein binding to the wild type TTR that explains, at least in part, the frequently observed negatively cooperative process between the two sites of TTRwt when binding ligands. These findings show that TTR mutants may present different ligand recognition and therefore are of value in ligand design for inhibiting TTR amyloidosis. (C) 2010 Elsevier Inc. All rights reserved.

Mineralocorticoid Receptor Mutations Differentially Affect Individual Gene Expression Profiles in Pseudohypoaldosteronism Type 1

Context: Type 1 pseudohypoaldosteronism (PHA1), a primary form of mineralocorticoid resistance, isdueto inactivating mutations of the NR3C2 gene, coding for the mineralocorticoid receptor (MR). Objective: The objective of the study was to assess whether different NR3C2 mutations have distinct effects on the pattern of MR-dependent transcriptional regulation of aldosterone-regulated genes. Design and Methods: Four MR mutations affecting residues in the ligand binding domain, identified in families with PHA1, were tested. MR proteins generated by site-directed mutagenesis were analyzed for their binding to aldosterone and were transiently transfected into renal cells to explore the functional effects on the transcriptional activity of the receptors by cis-trans-cotrans-activation assays and by measuring the induction of endogenous gene transcription. Results: Binding assays showed very low or absent aldosterone binding for mutants MR(877Pro), MR(848Pro), and MR(947stop) and decreased affinity for aldosterone of MR(843Pro). Compared with wildtype MR, the mutations p.Leu843Pro and p.Leu877Pro displayed half-maximal aldosterone-dependent transactivation of reporter genes driven by mouse mammary tumor virus or glucocorticoid response element-2 dependent promoters, whereas MR(848Pro) and MR(947stop) nearly or completely lost transcriptional activity. Although MR(848Pro) and MR(947stop) were also incapable of inducing aldosterone-dependent gene expression ofendogenoussgk1, GILZ, NDRG2, and SCNN1A, MR(843Pro) retained complete transcriptional activity on sgk1 and GILZ gene expression, and MR(877Pro) negatively affected the expression of sgk1, NDRG2, and SCNN1A. Conclusions: Our data demonstrate that MR mutations differentially affect individual gene expression in a promoter-dependent manner. Investigation of differential gene expression profiles in PHA1 may allow a better understanding of the molecular substrate of phenotypic variability and to elucidate pathogenic mechanisms underlying the disease. (J Clin Endocrinol Metab 96: E519-E527, 2011)

Analysis of Agonist and Antagonist Effects on Thyroid Hormone Receptor Conformation by Hydrogen/Deuterium Exchange

Thyroid hormone receptors (TRs) are ligand-gated transcription factors with critical roles in development and metabolism. Although x-ray structures of TR ligand-binding domains (LBDs) with agonists are available, comparable structures without ligand (apo-TR) or with antagonists are not. It remains important to understand apo-LBD conformation and the way that it rearranges with ligands to develop better TR pharmaceuticals. In this study, we conducted hydrogen/deuterium exchange on TR LBDs with or without agonist (T(3)) or antagonist (NH(3)). Both ligands reduce deuterium incorporation into LBD amide hydrogens, implying tighter overall folding of the domain. As predicted, mass spectroscopic analysis of individual proteolytic peptides after hydrogen/deuterium exchange reveals that ligand increases the degree of solvent protection of regions close to the buried ligand-binding pocket. However, there is also extensive ligand protection of other regions, including the dimer surface at H10-H11, providing evidence for allosteric communication between the ligand-binding pocket and distant interaction surfaces. Surprisingly, C-terminal activation helix H12, which is known to alter position with ligand, remains relatively protected from solvent in all conditions suggesting that it is packed against the LBD irrespective of the presence or type of ligand. T(3), but not NH(3), increases accessibility of the upper part of H3-H5 to solvent, and we propose that TR H12 interacts with this region in apo-TR and that this interaction is blocked by T(3) but not NH(3.) We present data from site-directed mutagenesis experiments and molecular dynamics simulations that lend support to this structural model of apo-TR and its ligand-dependent conformational changes. (Molecular Endocrinology 25: 15-31, 2011)

Probing conformational changes in orphan nuclear receptor: The NGFI-B intermediate is a partially unfolded dimer

Human nerve growth factor-induced B (NGFI-B) is a member of the NR4A subfamily of orphan nuclear receptors (NRs). Lacking identified ligands, orphan NRs show particular co-regulator proteins binding properties, different from other NRs, and they might have a non-classical quaternary organization. A body of evidence suggests that NRs recognition of and binding to ligands, DNA, homo- and heterodimerization partners and co-regulator proteins involve significant conformational changes of the NR ligand-binding domains (LBDs). To shed light on largely unknown biophysical properties of NGFI-B, here we studied structural organization and unfolding properties of NGFI-B ligand (like)-binding domain induced by chemical perturbation. Our results show that NGFI-B LBD undergoes a two-state guanidine hydrochloride (GndHCl) induced denaturation, as judged by changes in the a-helical content of the protein monitored by circular dichroism spectroscopy (CD). In contrast, changes in the tertiary structure of NGFI-B LBD, reported by intrinsic fluorescence, reveal a clear intermediate state. Additionally, SAXS results demonstrate that the intermediate observed by intrinsic fluorescence is a partially folded homodimeric structure, which further unfolds without dissociation at higher GndHCl concentrations. This partially unfolded dimeric assembly of NGFI-B LBD might resemble an intermediate that this domain access momentarily in the native state upon interactions with functional partners. (C) 2008 Elsevier B.V. All rights reserved.

Actions of octocoral and tobacco cembranoids on nicotinic receptors

Nicotinic acetylcholine receptors (AChRs) are pentameric proteins that form agonist-gated cation channels through the plasma membrane. AChR agonists and antagonists are potential candidates for the treatment of neurodegenerative diseases. Cembranoids are naturally occurring diterpenoids that contain a 14-carbon ring. These diterpenoids interact with AChRs in complex ways: as irreversible inhibitors at the agonist sites, as noncompetitive inhibitors, or as positive modulators, but no cembranoid was ever shown to have agonistic activity on AChRs. The cembranoid eupalmerin acetate displays positive modulation of agonist-induced currents in the muscle-type AChR and in the related gamma-aminobutyric acid (GABA) type A receptor. Moreover, cembranoids display important biological effects, many of them mediated by nicotinic receptors. Cembranoids from tobacco are neuroprotective through a nicotinic anti-apoptotic mechanism preventing excitotoxic neuronal death which in part could result from anti-inflammatory properties of cembranoids. Moreover, tobacco cembranoids also have anti-inflammatory properties which could enhance their neuroprotective properties. Cembranoids from tobacco affect nicotine-related behavior: they increase the transient initial ataxia caused by first nicotine injection into naive rats and inhibit the expression of locomotor sensitization to repeated injections of nicotine. In addition, cembranoids are known to act as anti-tumor compounds. In conclusion, cembranoids provide a promising source of lead drugs for many clinical areas, including neuroprotection, smoking-cessation, and anti-cancer therapies. (C) 2009 Elsevier Ltd. All rights reserved.

Estrogen receptor: Structural differences and potential implications on selectivity examined by the GRID/CPCA approach

Selective Estrogen Receptor Modulators ( SERMs) have been developed, but the selectivity towards the subtypes ( ER or ER is not well understood. Based on three-dimensional structural properties of ligand binding domains, a model that takes into account this aspect was developed via molecular interaction fields and consensus principal component analysis (GRID/CPCA).

Genetic variability and natural selection at the ligand domain of the Duffy binding protein in brazilian Plasmodium vivax populations

Background: Plasmodium vivax malaria is a major public health challenge in Latin America, Asia and Oceania, with 130-435 million clinical cases per year worldwide. Invasion of host blood cells by P. vivax mainly depends on a type I membrane protein called Duffy binding protein (PvDBP). The erythrocyte-binding motif of PvDBP is a 170 amino-acid stretch located in its cysteine-rich region II (PvDBP(II)), which is the most variable segment of the protein. Methods: To test whether diversifying natural selection has shaped the nucleotide diversity of PvDBP(II) in Brazilian populations, this region was sequenced in 122 isolates from six different geographic areas. A Bayesian method was applied to test for the action of natural selection under a population genetic model that incorporates recombination. The analysis was integrated with a structural model of PvDBP(II), and T-and B-cell epitopes were localized on the 3-D structure. Results: The results suggest that: (i) recombination plays an important role in determining the haplotype structure of PvDBP(II), and (ii) PvDBP(II) appears to contain neutrally evolving codons as well as codons evolving under natural selection. Diversifying selection preferentially acts on sites identified as epitopes, particularly on amino acid residues 417, 419, and 424, which show strong linkage disequilibrium. Conclusions: This study shows that some polymorphisms of PvDBP(II) are present near the erythrocyte-binding domain and might serve to elude antibodies that inhibit cell invasion. Therefore, these polymorphisms should be taken into account when designing vaccines aimed at eliciting antibodies to inhibit erythrocyte invasion.

Essential role of TM V and VI for binding the C-terminal sequences of Des-Arg-kinins

In the kallikrein-kinin and renin-angiotensin systems the main receptors, B-1 and B-2 (kinin receptors) and AT(1) and AT(2) (angiotensin receptors) respectively, are seven-transmembrane domain G-protein-coupled receptors. Considering that the B, agonists Des-Arg(9)-BK (Arg-Pro-Pro-Gly-Phe-Ser-Pro-Phe), Lys-desArg(9)-BK or Des-Arg(10)-KD (Lys-Arg-Pro-Pro-Gly-Phe-Ser-Pro-Phe) and the AT, agonist (Asp-Arg-Val-Tyr-lle-His-Pro-Phe) have the same two residues at the C-terminal region (i.e. Pro-Phe), we hypothesized that TM V and TM VI of the B-1 receptor could play an essential role in agonist binding and activity, being these regions receptor sites for binding the C-terminal sequences of Des-Arg-kinins similarly to that observed to AT, receptor. To investigate this hypothesis, we replaced Arg(212) for Ala at the top of the TM V and the sequence 274-282 (CPYHFFAFL) in TM VI of the rat kinin B, receptor by the 32 receptor homologous sequence, 289-297 (FPFQISTFL) and subsequently analyzed the consequences of these mutations by competition binding and functional assays. Despite correct expression, observed at the mRNA and protein level by RT-PCR and confocal microscopy, respectively, no agonist binding and function was verified for the mutated receptors. Therefore, our results suggest an important role for Arg(212) in the TM V and a region of TM VI of rat B, receptor in the interaction with the C-terminal residues of Des-Arg-kinins, similar to that observed with AngII. (c) 2007 Elsevier B.V. All rights reserved.

Differential effects of TR ligands on hormone dissociation rates: Evidence for multiple ligand entry/exit pathways

Some nuclear receptor (NR) ligands promote dissociation of radiolabeled bound hormone from the buried ligand binding cavity (LBC) more rapidly than excess unlabeled hormone itself This result was interpreted to mean that challenger ligands bind allosteric sites on the LBD to induce hormone dissociation, and recent findings indicate that ligands bind weakly to multiple sites on the LBD surface. Here we show, that a large fraction of thyroid hormone receptor (TR) ligands promote rapid dissociation (T(1/2) < 2 h) of , radiolabeled T(3) vs. T(3) (T(1/2), approximate to 5-7 h). We cannot discern relationships between this effect and ligand size, activity or affinity for TR beta. One ligand, GC-24, binds the TR LBC and (weakly) to the TR beta-LBD surface that mediates dimer/heterodimer interaction, but we cannot link this interaction to rapid T(3) dissociation. Instead, several lines of evidence suggest that the challenger ligand must interact with the buried LBC to promote rapid T(3) release. Since previous molecular dynamics simulations suggest that TR ligands leave the LBC by several routes, we propose that a subset of challenger ligands binds and stabilizes a partially unfolded intermediate state of TR that arises during T(3) release and that this effect enhances hormone dissociation. (C) 2009 Elsevier Ltd. All rights reserved.