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.

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)

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.

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.

Aminoacetone, a putative endogenous source of methylglyoxal, causes oxidative stress and death to insulin-producing RINm5f cells

Aminoacetone (AA), triose phosphates, and acetone are putative endogenous sources of potentially cytotoxic and genotoxic methylglyoxal (MG), which has been reported to be augmented in the plasma of diabetic patients. In these patients, accumulation of MG derived from aminoacetone, a threonine and glycine catabolite, is inferred from the observed concomitant endothelial overexpression of circulating semicarbazide-sensitive amine oxidases. These copper-dependent enzymes catalyze the oxidation of primary amines, such as AA and methylamine, by molecular oxygen, to the corresponding aldehydes, NH4+ ion and H2O2. We recently reported that AA aerobic oxidation to MG also takes place immediately upon addition of catalytic amounts of copper and iron ions. Taking into account that (i) MG and H2O2 are reportedly cytotoxic to insulin-producing cell lineages such as RINm5f and that (ii) the metal-catalyzed oxidation of AA is propagated by O-2(center dot-) radical anion, we decided to investigate the possible pro-oxidant action of AA on these cells taken here as a reliable model system for pancreatic beta-cells. Indeed, we show that AA (0.10-5.0 mM) administration to RINm5f cultures induces cell death. Ferrous (50-300 mu M) and Fe3+ ion (100 mu M) addition to the cell cultures had no effect, whereas Cu2+ (5.0-100 mu M) significantly increased cell death. Supplementation of the AA- and Cu2+-containing culture medium with antioxidants, such as catalase (5.0 mu M), superoxide dismutase (SOD, 50 U/mL), and N-acetylcysteine (NAC, 5.0 mM) led to partial protection. mRNA expression of MnSOD, CuZnSOD, glutathione peroxidase, and glutathione reductase, but not of catalase, is higher in cells treated with AA (0.50-1.0 mM) plus Cu2+ ions (10-50 mu M) relative to control cultures. This may imply higher activity of antioxidant enzymes C, in RINm5f AA-treated cells. In addition, we have found that AA (0.50-1.0 mM) Plus Cu2+ (100 mu M) (i) increase RINm5f cytosolic calcium; (ii) promote DNA fragmentation; and (iii) increase the pro-apoptotic (Bax)/antiapoptotic (Bcl-2) ratio at the level of mRNA expression. In conclusion, although both normal and pathological concentrations of AA are probably much lower than those used here, it is tempting to propose that excess AA in diabetic patients may drive oxidative damage and eventually the death of pancreatic beta-cells.

Human endogenous RNAs: Implications for the immunomodulation of Toll-like receptor 3

Toll-like receptors (TLRs), a family of mammalian receptors, are able to recognize nucleic acids. TLR3 recognizes double-stranded (ds)RNA, a product of the replication of certain viruses. Polyinosinic-polycytidylic acid, referred to as poly(I:C), an analog of viral dsRNA, interacts with TLR3 thereby eliciting immunoinflammatory responses characteristic of viral infection or down-regulating the expression of chemokine receptor CXCR4. It is known that dsRNA also directly activates interferon (IFN)-induced enzymes, such as the RNA-dependent protein kinase (PKR). In the present study, the mRNA expression of TLR3, CXCR4, IFN gamma and PKR was investigated in a culture of peripheral blood mononuclear cells (PBMCs) stimulated with poly(I:C) and endogenous RNA from human PBMCs. No cytotoxic effect on the cells or on the proliferation of CD3(+), CD4(+) and CD8(+) cells was observed. TLR3 expression in the PBMCs in the presence of poly(I:C) was up-regulated 9.5-fold, and TLR3 expression in the PBMCs treated with endogenous RNA was down-regulated 1.8-fold (p=0.002). The same trend was observed for IFN gamma where in the presence of poly(I:C) an 8.7-fold increase was noted and in the presence of endogenous RNA a 3.1-fold decrease was observed. In the culture activated with poly(1:C), mRNA expression of CXCR4 increased 8.0-fold and expression of PKR increased 33.0-fold. Expression of these genes decreased in the culture treated with endogenous RNA when compared to the culture without stimulus. Thus, high expression of mRNA for TLR3, IFN gamma, CXCR4 and PKR was observed in the presence of poly(I:C) and low expression was observed in the cells cultured with endogenous RNA. In conclusion, TLR3 may play major physiological roles that are not in the context of viral infection. It is possible that RNA released from cells could contain enough double-stranded structures to regulate cell activation. The involvement of endogenous RNA in endogenous gene expression and its implications in the regulation thereof, are still being studied, and will have significant implications in the future.