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.

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.

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.

Neither bovine somatotropin nor growth hormone-releasing factor alters expression of thyroid hormone receptors in liver and mammary tissues

Physiological effects of thyroid hormones are mediated primarily by binding of triiodothyronine to specific nuclear receptors. Organ-specific changes in production of triiodothyronine from its prohormone, thyroxine, have been hypothesized to target the action of thyroid hormones on the mammary gland and play a role in mediating or augmenting a galactopoietic response to bovine somatotropin (bST). Additionally, tissue responsiveness to thyroid hormones may be altered by changes in the number or affinity of nuclear receptors for thyroid hormones. In the present study, effects of bST and bovine growth hormone-releasing factor (bGRF) on thyroid hormone receptors in liver and mammary gland were studied. Lactating Holstein cows received continuous infusions of bST or bGRF for 63 d or served as uninfused controls. Nuclei were isolated from harvested mammary and liver tissues and incubated with [(125)I]-triiodothyronine. Treatments did not alter the capacity or affinity of specific binding sites for triiodothyronine in liver or mammary nuclei. Evaluation of transcript abundance for thyroid hormone receptors showed that isoforms of thyroid hormone receptor or retinoid receptor (which may influence thyroid receptor action) expressed in the mammary gland were not altered by bST or bGRF treatment. Data do not support the hypothesis that administration of bST or bGRF alters sensitivity of mammary tissue by changing expression of thyroid hormone receptors.

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.

The Thyroid Hormone Receptor (TR) beta-Selective Agonist GC-1 Inhibits Proliferation But Induces Differentiation and TR beta mRNA Expression in Mouse and Rat Osteoblast-Like Cells

Previous studies showed anabolic effects of GC-1, a triiodothyronine (T3) analogue that is selective for both binding and activation functions of thyroid hormone receptor (TR) beta 1 over TR alpha 1, on bone tissue in vivo. The aim of this study was to investigate the responsiveness of rat (ROS17/2.8) and mouse (MC3T3-E1) osteoblast-like cells to GC-1. As expected, T3 inhibited cellular proliferation and stimulated mRNA expression of osteocalcin or alkaline phosphatase in both cell lineages. Whereas equimolar doses of T3 and GC-1 equally affected these parameters in ROS17/2.8 cells, the effects of GC-1 were more modest compared to those of T3 in MC3T3-E1 cells. Interestingly, we showed that there is higher expression of TR alpha 1 than TR beta 1 mRNA in rat (similar to 20-90%) and mouse (similar to 90-98%) cell lineages and that this difference is even higher in mouse cells, which highlights the importance of TR alpha 1 to bone physiology and may partially explain the modest effects of GC-1 in comparison with T3 in MC3T3-E1 cells. Nevertheless, we showed that TR beta 1 mRNA expression increases (similar to 2.8- to 4.3-fold) as osteoblastic cells undergo maturation, suggesting a key role of TR beta 1 in mediating T3 effects in the bone forming cells, especially in mature osteoblasts. It is noteworthy that T3 and GC-1 induced TR beta 1 mRNA expression to a similar extent in both cell lineages (similar to 2- to 4-fold), indicating that both ligands may modulate the responsiveness of osteoblasts to T3. Taken together, these data show that TR beta selective T3 analogues have the potential to directly induce the differentiation and activity of osteoblasts.

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)

Thyroid Hormone Increases TGF-beta 1 in Cardiomyocytes Cultures Independently of Angiotensin II Type 1 and Type 2 Receptors

TH-induced cardiac hypertrophy in vivo is accompanied by increased cardiac Transforming Growth Factor-beta 1 (TGF-beta 1) levels, which is mediated by Angiotensin II type 1 receptors (AT1R) and type 2 receptors (AT2R). However, the possible involvement of this factor in TH-induced cardiac hypertrophy is unknown. In this study we evaluated whether TH is able to modulate TGF-beta 1 in isolated cardiac, as well as the possible contribution of AT1R and AT2R in this response. The cardiac fibroblasts treated with T(3) did not show alteration on TGF-beta 1 expression. However, cardiomyocytes treated with T(3) presented an increase in TGF-beta 1 expression, as well as an increase in protein synthesis. The AT1R blockade prevented the T(3)-induced cardiomyocyte hypertrophy, while the AT2R blockage attenuated this response. The T(3)-induced increase on TGF-beta 1 expression in cardiomyocytes was not changed by the use of AT1R and AT2R blockers. These results indicate that Angiotensin II receptors are not implicated in T(3)-induced increase on TGF-beta expression and suggest that the trophic effects exerted by T(3) on cardiomyocytes are not dependent on the higher TGF-beta 1 levels, since the AT1R and AT2R blockers were able to attenuate the T(3)-induced cardiomyocyte hypertrophy but were not able to attenuate the increase on TGF-beta 1 levels promoted by T(3).

Dissecting the Relation between a Nuclear Receptor and GATA: Binding Affinity Studies of Thyroid Hormone Receptor and GATA2 on TSH beta Promoter

Background: Much is known about how genes regulated by nuclear receptors (NRs) are switched on in the presence of a ligand. However, the molecular mechanism for gene down-regulation by liganded NRs remains a conundrum. The interaction between two zinc-finger transcription factors, Nuclear Receptor and GATA, was described almost a decade ago as a strategy adopted by the cell to up-or down-regulate gene expression. More recently, cell-based assays have shown that the Zn-finger region of GATA2 (GATA2-Zf) has an important role in down-regulation of the thyrotropin gene (TSH beta) by liganded thyroid hormone receptor (TR). Methodology/Principal Findings: In an effort to better understand the mechanism that drives TSH beta down-regulation by a liganded TR and GATA2, we have carried out equilibrium binding assays using fluorescence anisotropy to study the interaction of recombinant TR and GATA2-Zf with regulatory elements present in the TSH beta promoter. Surprisingly, we observed that ligand (T3) weakens TR binding to a negative regulatory element (NRE) present in the TSH beta promoter. We also show that TR may interact with GATA2-Zf in the absence of ligand, but T3 is crucial for increasing the affinity of this complex for different GATA response elements (GATA-REs). Importantly, these results indicate that TR complex formation enhances DNA binding of the TR-GATA2 in a ligand-dependent manner. Conclusions: Our findings extend previous results obtained in vivo, further improving our understanding of how liganded nuclear receptors down-regulate gene transcription, with the cooperative binding of transcription factors to DNA forming the core of this process.

Structural basis of GC-1 selectivity for thyroid hormone receptor isoforms

Background: Thyroid receptors, TRa and TR beta, are involved in important physiological functions such as metabolism, cholesterol level and heart activities. Whereas metabolism increase and cholesterol level lowering could be achieved by TR beta isoform activation, TRa activation affects heart rates. Therefore, beta-selective thyromimetics have been developed as promising drug-candidates for treatment of obesity and elevated cholesterol level. GC-1 [ 3,5-dimethyl-4-(4'-hydroxy- 3'-isopropylbenzyl)-phenoxy acetic acid] has ability to lower LDL cholesterol with 600-to 1400-fold more potency and approximately two-to threefold more efficacy than atorvastatin (Lipitor(C)) in studies in rats, mice and monkeys. Results: To investigate GC-1 specificity, we solved crystal structures and performed molecular dynamics simulations of both isoforms complexed with GC-1. Crystal structures reveal that, in TRa Arg228 is observed in multiple conformations, an effect triggered by the differences in the interactions between GC-1 and Ser277 or the corresponding asparagine (Asn331) of TR beta. The corresponding Arg282 of TR beta is observed in only one single stable conformation, interacting effectively with the ligand. Molecular dynamics support this model: our simulations show that the multiple conformations can be observed for the Arg228 in TR alpha, in which the ligand interacts either strongly with the ligand or with the Ser277 residue. In contrast, a single stable Arg282 conformation is observed for TR beta, in which it strongly interacts with both GC-1 and the Asn331. Conclusion: Our analysis suggests that the key factors for GC-1 selectivity are the presence of an oxyacetic acid ester oxygen and the absence of the amino group relative to T(3). These results shed light into the beta-selectivity of GC-1 and may assist the development of new compounds with potential as drug candidates to the treatment of hypercholesterolemia and obesity.

The importance of hormone receptor analysis in osteosarcoma cells growth submitted to treatment with estrogen in association with thyroid hormone

Bone tumor incidence in women peaks at age 50-60, coinciding with the menopause. That estrogen (E2) and triiodothyronine (T3) interact in bone metabolism has been well established. However, few data on the action of these hormones are available. Our purpose was to determine the role of E2 and T3 in the expression of bone activity markers, namely alkaline phosphatase (AP) and receptor activator of nuclear factor kappa B ligand (RANKL). Two osteosarcoma cell lines: MG-63 (which has both estrogen (ER) and thyroid hormone (TR) receptors) and SaOs-29 (ER receptors only) were treated with infraphysiological E2 associated with T3 at infraphysiological, physiological, and supraphysiological concentrations. Real-time RT-PCR was used for expression analysis. Our results show that, in MG-63 cells, infraphysiological E2 associated with supraphysiological T3 increases AP expression and decreases RANKL expression, while infraphysiological E2 associated with either physiological or supraphysiological T3 decreases both AP and RANKL expression. On the other hand, in SaOs-2 cells, the same hormone combinations had no significant effect on the markers` expression. Thus, the analysis of hormone receptors was shown to be crucial for the assessment of tumor potential growth in the face of hormonal changes. Special care should be provided to patients with T3 and E2 hormone receptors that may increase tumor growth. Copyright (C) 2007 John Wiley & Sons, Ltd.

Angiotensin type 1 receptor mediates thyroid hormone-induced cardiomyocyte hypertrophy through the Akt/GSK-3 beta/mTOR signaling pathway

Several studies have implicated the renin angiotensin system in the cardiac hypertrophy induced by thyroid hormone. However, whether Angiotensin type 1 receptor (AT(1)R) is critically required to the development of T(3)-induced cardiomyocyte hypertrophy as well as whether the intracellular mechanisms that are triggered by AT(1)R are able to contribute to this hypertrophy model is unknown. To address these questions, we employed a selective small interfering RNA (siRNA, 50 nM) or an AT(1)R blocker (Losartan, 1 mu M) to evaluate the specific role of this receptor in primary cultures of neonatal cardiomyocytes submitted to T(3) (10 nM) treatment. The cardiomyocytes transfected with the AT(1)R siRNA presented reduced mRNA (90%, P < 0.001) and protein (70%, P < 0.001) expression of AT(1)R. The AT(1)R silencing and the AT(1)R blockade totally prevented the T(3)-induced cardiomyocyte hypertrophy, as evidenced by lower mRNA expression of atrial natriuretic factor (66%, P < 0.01) and skeletal alpha-actin (170%, P < 0.01) as well as by reduction in protein synthesis (85%, P < 0.001). The cardiomyocytes treated with T(3) demonstrated a rapid activation of Akt/GSK-3 beta/mTOR signaling pathway, which was completely inhibited by the use of PI3K inhibitors (LY294002, 10 mu M and Wortmannin, 200 nM). In addition, we demonstrated that the AT(1)R mediated the T(3)-induced activation of Akt/GSK-3 beta/mTOR signaling pathway, since the AT(1)R silencing and the AT(1)R blockade attenuated or totally prevented the activation of this signaling pathway. We also reported that local Angiotensin I/II (Ang I/II) levels (120%, P < 0.05) and the AT(1)R expression (180%, P < 0.05) were rapidly increased by T(3) treatment. These data demonstrate for the first time that the AT(1)R is a critical mediator to the T(3)-induced cardiomyocyte hypertrophy as well as to the activation of Akt/GSK-3 beta/mTOR signaling pathway. These results represent a new insight into the mechanism of T(3)-induced cardiomyocyte hypertrophy, indicating that the Ang I/II-AT(1)R-Akt/GSK-3 beta/mTOR pathway corresponds to a potential mediator of the trophic effect exerted by T(3) in cardiomyocytes.

Recognition by the Thyroid Hormone Receptor of Canonical DNA Response Elements

To shed more light on the molecular requirements for recognition of thyroid response elements (TRES) by thyroid receptors (TRs), we compared the specific aspects of DNA TRE recognition by different TR constructs. Using fluorescence anisotropy, we performed a detailed and hierarchical study of TR-TRE binding. This wits done by comparing the binding affinities of three different TR constructs for four different TRE DNA elements, including palindromic sequences and direct repeats (F2, PAL, DR-1, and DR-4) as well as their interactions with nonspecific DNA sequences. The effect of MgCl(2) on suppressing of nonselective DNA binding to TR was also investigated. Furthermore, we determined the dissociation constants of the hTR beta DBD (DNA binding domain) and hTR beta DBD-LBD (DNA binding and ligand binding domains) for specific TRES. We found that a minimum DNA recognition peptide derived from DBD (H1TR) is sufficient for recognition and interaction with TREs, whereas scrambled DNA sequences were unrecognized. Additionally, we determined that the TR DBD binds to F2, PAL, and DR-4 with high affinity and similar K(d) values. The TR DBD-LBD recognizes all the tested TRES but binds preferentially to F2, with even higher affinity. Finally, our results demonstrate the important role played by LBDs in modulating TR-DNA binding.

Deiodinase-mediated thyroid hormone inactivation minimizes thyroid hormone signaling in the early development of fetal skeleton

Thyroid hormone (TH) plays a key role on post-natal bone development and metabolism, while its relevance during fetal bone development is uncertain. To Study this, pregnant once were made hypothyroid and fetuses harvested at embryonic days (E) 12.5, 14.5, 16.5 and 18.5. Despite a marked reduction in fetal tissue concentration of both T4 and T3, bone development, as assessed at the distal epiphyseal growth plate of the femur and vertebra, was largely preserved Lip to E16.5. Only at E18.5, the hypothyroid fetuses exhibited a reduction in femoral type I and type X collagen and osteocalcin mRNA levels, in the length and area of the proliferative and hypertrophic zones, in the number of chondrocytes per proliferative column, and in the number of hypertrophic chondrocyres, in addition to a slight delay in endochondral and intramembranous ossification. This Suggests that LIP to E 16.5, thyroid hormone signaling in bone is kept to a minimum. In fact, measuring the expression level of the activating and inactivating iodothyronine deiodinases (D2 and D3) helped understand how this is achieved. D3 mRNA was readily detected as early as E14.5 and its expression decreased markedly (similar to 10-fold) at E18.5, and even more at 14 days after birth (P14). In contrast. D2 mRNA expression increased significantly by E18.5 and markedly (similar to 2.5-fold) by P14. The reciprocal expression levels of D2 and D3 genes during early bone development along with the absence of a hypothyroidism-induced bone phenotype at this time Suggest that coordinated reciprocal deiodinase expression keeps thyroid hormone signaling in bone to very low levels at this early stage of bone development. (c) 2008 Elsevier Inc. All rights reserved.

Role of Halogen Bonds in Thyroid Hormone Receptor Selectivity: Pharmacophore-Based 3D-QSSR Studies

Most physiological effects of thyroid hormones are mediated by the two thyroid hormone receptor subtypes, TR alpha and TR beta. Several pharmacological effects mediated by TR beta might be beneficial in important medical conditions such as obesity, hypercholesterolemia and diabetes, and selective TR beta activation may elicit these effects while maintaining an acceptable safety profile, To understand the molecular determinants of affinity and subtype selectivity of TR ligands, we have successfully employed a ligand- and structure-guided pharmacophore-based approach to obtain the molecular alignment of a large series of thyromimetics. Statistically reliable three-dimensional quantitative structure-activity relationship (3D-QSAR) and three-dimensional quantitative structure-selectivity relationship (3D-QSSR) models were obtained using the comparative molecular field analysis (CoMFA) method, and the visual analyses of the contour maps drew attention to a number of possible opportunities for the development of analogs with improved affinity and selectivity. Furthermore, the 3D-QSSR analysis allowed the identification of a novel and previously unmentioned halogen bond, bringing new insights to the mechanism of activity and selectivity of thyromimetics.