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.

Gaining ligand selectivity in thyroid hormone receptors via entropy

Nuclear receptors are important targets for pharmaceuticals, but similarities between family members cause difficulties in obtaining highly selective compounds. Synthetic ligands that are selective for thyroid hormone (TH) receptor beta (TR beta) vs. TR alpha reduce cholesterol and fat without effects on heart rate; thus, it is important to understand TR beta-selective binding. Binding of 3 selective ligands (GC-1, KB141, and GC-24) is characterized at the atomic level; preferential binding depends on a nonconserved residue (Asn-331 beta) in the TR beta ligand-binding cavity (LBC), and GC-24 gains extra selectivity from insertion of a bulky side group into an extension of the LBC that only opens up with this ligand. Here we report that the natural TH 3,5,3`-triodothyroacetic acid (Triac) exhibits a previously unrecognized mechanism of TR beta selectivity. TR x-ray structures reveal better fit of ligand with the TR alpha LBC. The TR beta LBC, however, expands relative to TR alpha in the presence of Triac (549 angstrom(3) vs. 461 angstrom(3)), and molecular dynamics simulations reveal that water occupies the extra space. Increased solvation compensates for weaker interactions of ligand with TR beta and permits greater flexibility of the Triac carboxylate group in TR beta than in TR alpha. We propose that this effect results in lower entropic restraint and decreases free energy of interactions between Triac and TR beta, explaining subtype-selective binding. Similar effects could potentially be exploited in nuclear receptor drug design.

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.

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.

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.

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)

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.

TLR4 and CD14 receptors expressed in rat pineal gland trigger NFKB pathway

Nuclear factor-kappa B (NFKB), a pivotal player in inflammatory responses, is constitutively expressed in the pineal gland. Corticosterone inhibits pineal NFKB leading to an enhancement of melatonin production, while tumor necrosis factor (TNF) leads to inhibition of Aa-nat transcription and the production of N-acetylserotonin in cultured glands. The reduction in nocturnal melatonin surge favors the mounting of the inflammatory response. Despite these data, there is no clear evidence of the ability of the pineal gland to recognize molecules that signal infection. This study investigated whether the rat pineal gland expresses receptors for lipopolysaccharide (LPS), the endotoxin from the membranes of Gram-negative bacteria, and to establish the mechanism of action of LPS. Here, we show that pineal glands possess both CD14 and toll-like receptor 4 (TLR4), membrane proteins that bind LPS and trigger the NFKB pathway. LPS induced the nuclear translocation of p50/p50 and p50/RELA dimers and the synthesis of TNF. The maximal expression of TNF in cultured glands coincides with an increase in the expression of TNF receptor 1 (TNFR1) in isolated pinealocytes. In addition, LPS inhibited the synthesis of N-acetylserotonin and melatonin. Therefore, the pineal gland transduces Gram-negative endotoxin stimulation by producing TNF and inhibiting melatonin synthesis. Here, we provide evidence to reinforce the idea of an immune-pineal axis, showing that the pineal gland is a constitutive player in the innate immune response.

Expression of Circadian Clock and Melatonin Receptors within Cultured Rat Cardiomyocytes

Melatonin, the pineal gland hormone, provides entrainment of many circadian rhythms to the ambient light/dark cycle. Recently, cardiovascular studies have demostrated melatonin interactions with many physiological processes and diseases, such as hypertension and cardiopathologies. Although membrane melatonin receptors (MT1, MT2) and the transcriptional factor ROR alpha have been reported to be expressed in the heart, there is no evidence of the cell-type expressing receptors as well as the possible role of melatonin on the expression of the circadian clock of cardiomyocytes, which play an important role in cardiac metabolism and function. Therefore, the aim of this study was to evaluate the mRNA and protein expressions of MT1, MT2, and ROR alpha and to determine whether melatonin directly influences expression of circadian clocks within cultured rat cardiomyocytes. Adult rat cardiomyocyte cultures were created, and the cells were stimulated with 1 nM melatonin or vehicle. Gene expressions were assayed by real-time polymerase chain reaction (PCR). The mRNA and protein expressions of membrane melatonin receptors and RORa were established within adult rat cardiomyocytes. Two hours of melatonin stimulation did not alter the expression pattern of the analyzed genes. However, given at the proper time, melatonin kept Rev-erb alpha expression chronically high, specifically 12 h after melatonin treatment, avoiding the rhythmic decline of Rev-erb alpha mRNA. The blockage of MT1 and MT2 by luzindole did not alter the observed melatonin-induced expression of Rev-erb alpha mRNA, suggesting the nonparticipation of MT1 and MT2 on the melatonin effect within cardiomyocytes. It is possible to speculate that melatonin, in adult rat cardiomyocytes, may play an important role in the light signal transduction to peripheral organs, such as the heart, modulating its intrinsic rhythmicity. (Author correspondence: cipolla@icb.usp.br)

Amyloid beta-peptide activates nuclear factor-kappa B through an N-methyl-D-aspartate signaling pathway in cultured cerebellar cells

Amyloid P-peptide (A beta) likely causes functional alterations in neurons well prior to their death. Nuclear factor-kappa B (NF-kappa B), a transcription factor that is known to play important roles in cell survival and apoptosis, has been shown to be modulated by A beta in neurons and glia, but the mechanism is unknown. Because A beta has also been shown to enhance activation of N-methyl-D-aspartate (NMDA) receptors, we investigated the role of NMDA receptor-mediated intracellular signaling pathways in A beta-induced NF-kappa B activation in primary cultured rat cerebellar cells. Cells were treated with different concentrations of A beta 1-40 (1 or 2 mu M) for different periods (6, 12, or 24 hr). MK-801 (NMDA antagonist), manumycin A and FTase inhibitor 1 (farnesyltransferase inhibitors), PP1 (Src-family tyrosine kinase inhibitor), PD98059 [mitogen-activated protein kinase (MAPK) inhibitor], and LY294002 [phosphatidylinositol 3-kinase (PI3-k) inhibitor] were added 20 min before A beta treatment of the cells. A beta induced a time- and concentration-dependent activation of NF-kappa B (1 mu M, 12 hr); both p50/p65 and p50/p50 NF-kappa B dimers were involved. This activation was abolished by MK-801 and attenuated by manumycin A, FTase inhibitor 1, PP1, PD98059, and LY294002. AP at 1 mu M increased the expression of inhibitory protein I kappa B, brain-derived neurotrophic factor, inducible nitric oxide synthase, tumor necrosis factor-alpha, and interleukin-1 beta as shown by RTPCR assays. Collectively, these findings suggest that AP activates NF-kappa B by an NMDA-Src-Ras-like protein through MAPK and PI3-k pathways in cultured cerebellar cells. This pathway may mediate an adaptive, neuroprotective response to A beta. (c) 2007 Wiley-Liss, Inc.

TRPV1 receptors are involved in protein nitration and Muller cell reaction in the acutely axotomized rat retina

We report here the protein expression of TRPV1 receptor in axotomized rat retinas and its possible participation in mechanisms involved in retinal ganglion cell (RGC) death. Adult rats were subjected to unilateral, intraorbital axotomy of the optic nerve, and the retinal tissue was removed for further processing. TRPV1 total protein expression decreased progressively after optic nerve transection, reaching 66.2% of control values 21 days after axotomy. The number of cells labeled for TRPV1 in the remnant GCL decreased after 21 days post-lesion (to 63%). Fluoro-jade B staining demonstrated that the activation of TRPV1 in acutely-lesioned eyes elicited more intense neuronal degeneration in the GCL and in the inner nuclear layer than in sham-operated retinas. A single intraocular injection of capsazepine (100 mu M), a TRPV1 antagonist, 5 days after optic nerve lesion, decreased the number of GFAP-expressing Muller cells (72.5% of control values) and also decreased protein nitration in the retinal vitreal margin (75.7% of control values), but did not affect lipid peroxidation. Furthermore, retinal explants were treated with capsaicin (100 mu M), and remarkable protein nitration was then present, which was reduced by blockers of the constitutive and inducible nitric oxide synthases (7-NI and aminoguanidine, respectively). TRPV1 activation also increased GFAP expression, which was reverted by both TRPV1 antagonism with capsazepine and by 7-NI and aminoguanidine. Given that Muller cells do not express TRPV1, we suppose that the increased GFAP expression in these cells might be elicited by TRPV1 activation and by its indirect effect upon nitric oxide overproduction and peroxynitrite formation. We incubated Fluorogold pre-labeled retinal explants in the presence of capsazepine (1 mu M) during 48 h. The numbers of surviving RGCs stained with fluorogold and the numbers of apoptotic cells in the GCL detected with TUNEL were similar in lesioned and control retinas. We conclude that TRPV1 receptor expression decreased after optic nerve injury due to death of TRPV1-containing cells. Furthermore, these data indicate that TRPV1 might be involved in intrinsic protein nitration and Muller cell reaction observed after optic nerve injury. (C) 2010 Elsevier Ltd. All rights reserved.

Ontogenetic expression of the vanilloid receptors TRPV1 and TRPV2 in the rat retina

The present study aimed to analyze the gene and protein expression and the pattern of distribution of the vanilloid receptors TRPV1 and TRPV2 in the developing rat retina. During the early phases of development, TRPV1 was found mainly in the neuroblastic layer of the retina and in the pigmented epithelium. In the adult, TRPV1 was found in microglial cells, blood vessels, astrocytes and in neuronal structures, namely synaptic boutons of both retina] plexiform layers, as well as in cell bodies of the inner nuclear layer and the ganglion cell layer. The pattern of distribution of TRPV1 was mainly punctate, and there was higher TRPV1 labeling in the peripheral retina than in central regions. TRPV2 expression was quite distinct. its expression was virtually undetectable by immunoblotting before P1, and that receptor was found by immunohistochemistry only by postnatal day 15 (PI 5). RNA and protein analysis showed that the adult levels are only reached by P60, which includes small processes in the retinal plexiform layers, and labeled cellular bodies in the inner nuclear layer and the ganglion cell layer. There was no overlapping between the signal observed for both receptors. in conclusion, our results showed that the patterns of distribution of TRPV1 and TRPV2 are different during the development of the rat retina, suggesting that they have specific roles in both visual processing and in providing specific cues to neural development. (C) 2009 ISDN. Published by Elsevier Ltd. All rights reserved.

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.