Structural Rearrangements in the Thyroid Hormone Receptor Hinge Domain and Their Putative Role in the Receptor Function

Molecular Dynamic Simulation To Reveal the Mechanism Underlying MGL-3196 Resistance to Thyroxine Receptor Beta

Thyroxine receptor beta (TRβ) is a ligand-dependent nuclear receptor that participates in regulating multiple biological processes, particularly playing an important role in lipid metabolism regulation. TRβ is currently a popular therapeutic target for nonalcoholic steatohepatitis (NASH), while no drugs have been approved to treat this disease. MGL-3196 (Resmetirom) is the first TRβ agonist that has succeeded in phase III clinical trials for the treatment of NASH; therefore, studying its molecular mechanism of action is of great significance. In this study, we employed molecular dynamic simulation to investigate the interaction mode between MGL-3196 and TRβ at the all-atom level. More importantly, by comparing the binding patterns of MGL-3196 in several prevalent TRβ mutants, it was identified that the mutations R243Q and H435R located, respectively, around and within the ligand-binding pocket of TRβ cause TRβ to be insensitive to MGL-3196. This indicates that patients with NASH carrying these two mutations may exhibit resistance to the medication of MGL-3196, thereby highlighting the potential impact of TRβ mutations on TRβ-targeted treatment of NASH and beyond.

Ligand binding free energy evaluation by Monte Carlo Recursion

The correct evaluation of ligand binding free energies by computational methods is still a very challenging active area of research. The most employed methods for these calculations can be roughly classified into four groups: (i) the fastest and less accurate methods, such as molecular docking, designed to sample a large number of molecules and rapidly rank them according to the potential binding energy; (ii) the second class of methods use a thermodynamic ensemble, typically generated by molecular dynamics, to analyze the endpoints of the thermodynamic cycle for binding and extract differences, in the so-called 'end-point' methods; (iii) the third class of methods is based on the Zwanzig relationship and computes the free energy difference after a chemical change of the system (alchemical methods); and (iv) methods based on biased simulations, such as metadynamics, for example. These methods require increased computational power and as expected, result in increased accuracy for the determination of the strength of binding. Here, we describe an intermediate approach, based on the Monte Carlo Recursion (MCR) method first developed by Harold Scheraga. In this method, the system is sampled at increasing effective temperatures, and the free energy of the system is assessed from a series of terms W(b,T), computed from Monte Carlo (MC) averages at each iteration. We show the application of the MCR for ligand binding with datasets of guest-hosts systems (N = 75) and we observed that a good correlation is obtained between experimental data and the binding energies computed with MCR. We also compared the experimental data with an end-point calculation from equilibrium Monte Carlo calculations that allowed us to conclude that the lower-energy (lower-temperature) terms in the calculation are the most relevant to the estimation of the binding energies, resulting in similar correlations between MCR and MC data and the experimental values. On the other hand, the MCR method provides a reasonable view of the binding energy funnel, with possible connections with the ligand binding kinetics, as well. The codes developed for this analysis are publicly available on GitHub as a part of the LiBELa/MCLiBELa project (https://github.com/alessandronascimento/LiBELa).

Discovery of Highly Potent and Selective Thyroid Hormone Receptor β Agonists for the Treatment of Nonalcoholic Steatohepatitis

Nonalcoholic steatohepatitis (NASH) is a progressive stage of nonalcoholic fatty liver disease (NAFLD) and is characterized by steatosis, inflammation, hepatocyte ballooning, and fibrosis. While there are currently no approved therapies for NASH, the thyroid hormone receptor β (THR-β), primarily expressed in the liver, is emerging as an effective molecular target for the treatment of NASH. However, the adverse cardiac and bone effects mediated by thyroid hormone receptor α (THR-α) need to be minimized. Herein, we reported the discovery of a series of novel THR-β agonists featuring pyrrolo[3,2-b]pyridin-5-one skeletons based on structure-based drug design. Further optimization led to compound 15, which exhibited higher potency and selectivity for THR-β over THR-α compared to clinical drug MGL-3196. More significantly, an excellent liver-to-serum ratio of 93:1 was observed for compound 15. We believe that the high hepatic concentration of compound 15 may result in no cardiotoxicity.

A detailed theoretical exploration on the THR-β binding affinities and antioxidant activity of some halogenated bisphenols

Natural halogenated phenolic compounds are unique bioactive structures which share features and physicochemical properties with thyroid hormones, who are essential regulators of neurological development and metabolism processes. Also, these structures can be used as natural antioxidants to minimize the diseases related to oxidative stress. In this work, the binding affinity of 32 natural and synthetic halogenated bisphenols were investigated on thyroid hormone receptor-β (THR-β) using the molecular docking, MM/GBSA, molecular dynamics, and a rigorous three-layer ONIOM ((M06-2X/6-31G*:PM6:AMBER) calculation. The desirable potency is observed for binding of selected compounds to THR-β. The Met313, Asn331, and His435 are the main interacting residues in the binding cavity which involved in the hydrogen and halogen bond interactions with the ligands. The most potent candidate on binding to the active site of THR-β is presented with respect to the results of mentioned calculations. Moreover, the antioxidant activity of compounds has been investigated using the quantum mechanical calculations. The electrostatic potential surfaces illustrate well the antioxidant capacity of compounds. The halogen substituents increase the antioxidant activity of the most stable conformers. The position and number of OH groups are crucial factors which affect the activity, whereas two adjacent hydroxyl groups enhance the antioxidant activity of selected compounds.Communicated by Ramaswamy H. Sarma.

Thyroid hormone resistance: Mechanisms and therapeutic development

As an essential primary hormone, thyroid hormone (TH) is indispensable for human growth, development and metabolism. Impairment of TH function in several aspects, including TH synthesis, activation, transportation and receptor-dependent transactivation, can eventually lead to thyroid hormone resistance syndrome (RTH). RTH is a rare syndrome that manifests as a reduced target cell response to TH signaling. The majority of RTH cases are related to thyroid hormone receptor β (TRβ) mutations, and only a few RTH cases are associated with thyroid hormone receptor α (TRα) mutations or other causes. Patients with RTH suffer from goiter, mental retardation, short stature and bradycardia or tachycardia. To date, approximately 170 mutated TRβ variants and more than 20 mutated TRα variants at the amino acid level have been reported in RTH patients. In addition to these mutated proteins, some TR isoforms can also reduce TH function by competing with primary TRs for TRE and RXR binding. Fortunately, different treatments for RTH have been explored with structure-activity relationship (SAR) studies and drug design, and among these treatments. With thyromimetic potency but biochemical properties that differ from those of primary TH (T3 and T4), these TH analogs can bypass specific defective transporters or reactive mutant TRs. However, these compounds must be carefully applied to avoid over activating TRα, which is associated with more severe heart impairment. The structural mechanisms of mutation-induced RTH in the TR ligand-binding domain are summarized in this review. Furthermore, strategies to overcome this resistance for therapeutic development are also discussed.

Main Factors Involved in Thyroid Hormone Action

The thyroid hormone receptors are the mediators of a multitude of actions by the thyroid hormones in cells. Most thyroid hormone activities require interaction with nuclear receptors to bind DNA and regulate the expression of target genes. In addition to genomic regulation, thyroid hormones function via activation of specific cytosolic pathways, bypassing interaction with nuclear DNA. In the present work, we reviewed the most recent literature on the characteristics and roles of different factors involved in thyroid hormone function in particular, we discuss the genomic activity of thyroid hormone receptors in the nucleus and the functions of different thyroid hormone receptor isoforms in the cytosol. Furthermore, we describe the integrin αvβ3-mediated thyroid hormone signaling pathway and its rapid nongenomic action in the cell. We furthermore reviewed the thyroid hormone transporters enabling the uptake of thyroid hormones in the cell, and we also include a paragraph on the proteins that mediate thyroid receptors’ shuttling from the nucleus to the cytosol.

The clinical significance of the glucocorticoid receptors: Genetics and epigenetics

The impacts of glucocorticoids (GCs) are mainly mediated by a nuclear receptor (GR) existing in almost every tissue. The GR regulates a wide range of physiological functions, including inflammation, cell metabolism, and differentiation playing a major role in cellular responses to GCs and stress. Therefore, the dysregulation or disruption of GR can cause deficiencies in the adaptation to stress and the preservation of homeostasis. The number of GR polymorphisms associated with different diseases has been mounting per year. Tackling these clinical complications obliges a comprehensive understanding of the molecular network action of GCs at the level of the GR structure and its signaling pathways. Beyond genetic variation in the GR gene, epigenetic changes can enhance our understanding of causal factors involved in the development of diseases and identifying biomarkers. In this review, we highlight the relationships of GC receptor gene polymorphisms and epigenetics with different diseases.

Molecular modeling and molecular dynamics simulation studies on thyroid hormone receptor from Rattus norvegicus: role of conserved water molecules

Thyroid hormone receptor (THR) belongs to the nuclear receptor (NR) superfamily that is activated by binding of appropriate ligand molecules (thyroid hormones). These receptors directly bind to specific DNA sequences for gene expression, which is essential for metabolism, homeostasis, and the development of organisms, making it an important drug target. Extensive MD-simulation studies of triiodothyronine (T3) docked modeled rnTHRβ1 structures have indicated the presence of twelve conserved water molecules at the DNA-DBD (DNA binding domain) interface. The W1-W5 water centers have been involved in the recognition between the A-chain of DBD to C-chain of DNA, W6 and W7 mediated the interaction between A-chain of DBD and D-chain of DNA, W8 and W9 recognized the B-chain of DBD and C-chain of DNA, and W9-W12 centers conjugated the residues of B-chain of DBD to D-chain of DNA through hydrogen bonds. The conformation flexibility of Phe272 and Met313 residues in the absence of T3 at the LBD (ligand-binding domain) region have been observed and reported.

Nuclear hormone and peptide hormone therapeutics for NAFLD and NASH

Background

Non-alcoholic steatohepatitis (NASH) is a spectrum of histological liver pathologies ranging from hepatocyte fat accumulation, hepatocellular ballooning, lobular inflammation, and pericellular fibrosis. Based on early investigations, it was discovered that visceral fat accumulation, hepatic insulin resistance, and atherogenic dyslipidemia are pathological triggers for NASH progression. As these pathogenic features are common with obesity, type 2 diabetes (T2D), and atherosclerosis, therapies that target dysregulated core metabolic pathways may hold promise for treating NASH, particularly as first-line treatments.

Scope of Review

In this review, the latest clinical data on nuclear hormone- and peptide hormone-based drug candidates for NASH are reviewed and contextualized, culminating with a discovery research perspective on emerging combinatorial therapeutic approaches that merge nuclear and peptide strategies.

Major Conclusion

Several drug candidates targeting the metabolic complications of NASH have shown promise in early clinical trials, albeit with unique benefits and challenges, but questions remain regarding their translation to larger and longer clinical trials, as well as their utility in a more diseased patient population. Promising polypharmacological approaches can potentially overcome some of these perceived challenges, as has been suggested in preclinical models, but deeper characterizations are required to fully evaluate these opportunities.

•

Despite no approved treatments for NASH, several drug candidates have shown promise in early clinical trials.

•

Therapies targeting metabolic pathologies of NASH have shown efficacy to reduce hepatic fat content and improve fibrosis.

•

Many of these therapies have been rationally designed to mimic nuclear hormone or peptide hormone action.

•

Despite provocative preclinical findings of nuclear and peptide hormone combination, clinical translation remains unproven.

Thyroxine binding to type III iodothyronine deiodinase

Iodothyronine deiodinases (Dios) are important selenoproteins that control the concentration of the active thyroid hormone (TH) triiodothyronine through regioselective deiodination. The X-ray structure of a truncated monomer of Type III Dio (Dio3), which deiodinates TH inner rings through a selenocysteine (Sec) residue, revealed a thioredoxin-fold catalytic domain supplemented with an unstructured Ω-loop. Loop dynamics are driven by interactions of the conserved Trp207 with solvent in multi-microsecond molecular dynamics simulations of the Dio3 thioredoxin(Trx)-fold domain. Hydrogen bonding interactions of Glu200 with residues conserved across the Dio family anchor the loop’s N-terminus to the active site Ser-Cys-Thr-Sec sequence. A key long-lived loop conformation coincides with the opening of a cryptic pocket that accommodates thyroxine (T₄) through an I⋯Se halogen bond to Sec170 and the amino acid group with a polar cleft. The Dio3-T₄ complex is stabilized by an I⋯O halogen bond between an outer ring iodine and Asp211, consistent with Dio3 selectivity for inner ring deiodination. Non-conservation of residues, such as Asp211, in other Dio types in the flexible portion of the loop sequence suggests a mechanism for regioselectivity through Dio type-specific loop conformations. Cys168 is proposed to attack the selenenyl iodide intermediate to regenerate Dio3 based upon structural comparison with related Trx-fold proteins.

Secretory expression of thyroid hormone receptor using transgenic silkworms and its DNA binding activity

Silkworms are economically important insects that have the ability to produce large amounts of silk. They have mass breeding methods and silk glands, which are specialized tissues that secrete silk fibroin and sericin. Thus, the production of recombinant proteins in a transgenic silkworm system is a promising approach. We developed a silkworm, Bombyx mori, as a host expression insect for recombinant proteins and successfully produced different proteins including antibodies, glycoproteins, and membrane receptors. The thyroid hormone receptor (TR) is a regulatory factor for many physiological phenomena. It is a lipophilic protein that has DNA-binding and ligand-binding domains. Based on our previous experiences, it was inferred that the recombinant TR easily formed aggregates and precipitates which is potentially due to an unstructured hinge domain. We applied the silkworm expression system to produce mice TRβ1 that was fused with glutathione S-transferase. Using 160 larvae, the yield of the recombinant GST-TRβ was approximately 4 mg, and the purified GST-TRβ completely retained its physiological activity. Our results indicated that the recombinant TRβ was secreted extracellularly using the silk fibroin signal peptide sequence. Moreover, we found that the expression system of silkworms was applicable to nuclear proteins.

Selective thyroid hormone receptor β agonists with oxadiazolone acid isosteres

Use of the oxadiazolone acid isostere in triiodothyronine analogs yielded potent and selective agonists for the thyroid hormone receptor β. Selected examples showed good in-vivo efficacy in a rat hypercholesterolemic model. One compound was further profiled in a diet-induced mouse model of nonalcoholic steatohepatitis (NASH) and showed robust target engagement and significant histological improvements in both liver steatosis and fibrosis.

Two Novel Cases of Resistance to Thyroid Hormone Due to THRA Mutation

Resistance to thyroid hormone alpha (RTHα) is a rare and under-recognized genetic disease caused by mutations of THRA, the gene encoding thyroid hormone receptor α1 (TRα1). We report here two novel THRA missense mutations (M259T, T273A) in patients with RTHα. We combined biochemical and cellular assays with in silico modeling to assess the capacity of mutant TRα1 to bind triiodothyronine (T3), to heterodimerize with RXR, to interact with transcriptional coregulators, and to transduce a T3 transcriptional response. M259T, and to a lower extent T273A, reduces the affinity of TRα1 for T3. Their negative influence is only reverted by large excess of T3. The severity of the two novel RTHα cases originates from a reduction in the binding affinity of TRα1 mutants to T3 and thus correlates with the incapacity of corepressors to dissociate from TRα1 mutants in the presence of T3.

Genomic and Non-Genomic Mechanisms of Action of Thyroid Hormones and Their Catabolite 3,5-Diiodo-L-Thyronine in Mammals

Since the realization that the cellular homologs of a gene found in the retrovirus that contributes to erythroblastosis in birds (v-erbA), i.e. the proto-oncogene c-erbA encodes the nuclear receptors for thyroid hormones (THs), most of the interest for THs focalized on their ability to control gene transcription. It was found, indeed, that, by regulating gene expression in many tissues, these hormones could mediate critical events both in development and in adult organisms. Among their effects, much attention was given to their ability to increase energy expenditure, and they were early proposed as anti-obesity drugs. However, their clinical use has been strongly challenged by the concomitant onset of toxic effects, especially on the heart. Notably, it has been clearly demonstrated that, besides their direct action on transcription (genomic effects), THs also have non-genomic effects, mediated by cell membrane and/or mitochondrial binding sites, and sometimes triggered by their endogenous catabolites. Among these latter molecules, 3,5-diiodo-L-thyronine (3,5-T2) has been attracting increasing interest because some of its metabolic effects are similar to those induced by T3, but it seems to be safer. The main target of 3,5-T2 appears to be the mitochondria, and it has been hypothesized that, by acting mainly on mitochondrial function and oxidative stress, 3,5-T2 might prevent and revert tissue damages and hepatic steatosis induced by a hyper-lipid diet, while concomitantly reducing the circulating levels of low density lipoproteins (LDL) and triglycerides. Besides a summary concerning general metabolism of THs, as well as their genomic and non-genomic effects, herein we will discuss resistance to THs and the possible mechanisms of action of 3,5-T2, also in relation to its possible clinical use as a drug.

Thyroid Hormone Transporters

Thyroid hormone transporters at the plasma membrane govern intracellular bioavailability of thyroid hormone. Monocarboxylate transporter (MCT) 8 and MCT10, organic anion transporting polypeptide (OATP) 1C1, and SLC17A4 are currently known as transporters displaying the highest specificity toward thyroid hormones. Structure-function studies using homology modeling and mutational screens have led to better understanding of the molecular basis of thyroid hormone transport. Mutations in MCT8 and in OATP1C1 have been associated with clinical disorders. Different animal models have provided insight into the functional role of thyroid hormone transporters, in particular MCT8. Different treatment strategies for MCT8 deficiency have been explored, of which thyroid hormone analogue therapy is currently applied in patients. Future studies may reveal the identity of as-yet-undiscovered thyroid hormone transporters. Complementary studies employing animal and human models will provide further insight into the role of transporters in health and disease. (Endocrine Reviews 41: 1 - 55, 2020).

The Race to Bash NASH: Emerging Targets and Drug Development in a Complex Liver Disease

Nonalcoholic steatohepatitis (NASH) is a severe form of nonalcoholic fatty liver disease (NAFLD) characterized by liver steatosis, inflammation, and hepatocellular damage. NASH is a serious condition that can progress to cirrhosis, liver failure, and hepatocellular carcinoma. The association of NASH with obesity, type 2 diabetes mellitus, and dyslipidemia has led to an emerging picture of NASH as the liver manifestation of metabolic syndrome. Although diet and exercise can dramatically improve NASH outcomes, significant lifestyle changes can be challenging to sustain. Pharmaceutical therapies could be an important addition to care, but currently none are approved for NASH. Here, we review the most promising targets for NASH treatment, along with the most advanced therapeutics in development. These include targets involved in metabolism (e.g., sugar, lipid, and cholesterol metabolism), inflammation, and fibrosis. Ultimately, combination therapies addressing multiple aspects of NASH pathogenesis are expected to provide benefit for patients.

The In Vitro Functional Impairment of Thyroid Hormone Receptor Alpha 1 Isoform Mutants Is Mainly Dictated by Reduced Ligand Sensitivity

Background: Thyroid hormone (TH) acts on TH receptors (TRs) and regulates gene transcription by binding of TRs to TH response elements (TREs) in target gene promoters. The transcriptional activity of TRs is modulated by interactions with TR-coregulatory proteins. Mutations in TRα cause resistance to thyroid hormone alpha (RTHα). In this study, we analyzed if, beyond reduced triiodothyronine (T3) affinity, altered interactions with cofactors or different TREs could account for the differential impaired transcriptional activity of different mutants. Methods: We evaluated four mutants derived from patients (D211G, M256T, A263S, and R384H) and three artificial mutants at equivalent positions in patients with RTHβ (T223A, L287V, and P398H). The in vitro transcriptional activity was evaluated on TRE-luciferase reporters (DR4, IR0, and ER6). The affinity for T3 and interaction with coregulatory proteins (nuclear receptor corepressor 1 [NCoR1] and steroid receptor coactivator 1 [SRC1]) were also determined. Results: We found that the affinity for T3 was significantly reduced for all mutants, except for TRα1-T223A. The reduction in the T3 sensitivity of the transcriptional activity on three TREs, the dissociation of the corepressor NCoR1, and the association of the coactivator SRC1 recruitment for each mutant correlated with the reduced affinity for T3. We did not observe mutation-specific alterations in interactions with cofactors or TREs. Conclusions: In summary, the degree of impaired transcriptional activity of mutants is mainly determined by their reduced affinity for T3.

Revealing a Mutant-Induced Receptor Allosteric Mechanism for the Thyroid Hormone Resistance

Resistance to thyroid hormone (RTH) is a clinical disorder without specific and effective therapeutic strategy, partly due to the lack of structural mechanisms for the defective ligand binding by mutated thyroid hormone receptors (THRs). We herein uncovered the prescription drug roxadustat as a novel THRβ-selective ligand with therapeutic potentials in treating RTH, thereby providing a small molecule tool enabling the first probe into the structural mechanisms of RTH. Despite a wide distribution of the receptor mutation sites, different THRβ mutants induce allosteric conformational modulation on the same His435 residue, which disrupts a critical hydrogen bond required for the binding of thyroid hormones. Interestingly, roxadustat retains hydrophobic interactions with THRβ via its unique phenyl extension, enabling the rescue of the activity of the THRβ mutants. Our study thus reveals a critical receptor allosterism mechanism for RTH by mutant THRβ, providing a new and viable therapeutic strategy for the treatment of RTH.

•

We identified a novel THR ligand that effectively binds to THRβ mutants

•

Structures revealed mechanisms for the RTH controlled by a key residue switch

•

Roxadustat retains unique hydrophobic interactions with THRβ mutants

•

We provide a promising approach to design THR ligands in treating RTH

Acetylation modulates thyroid hormone receptor intracellular localization and intranuclear mobility

The thyroid hormone receptor (TR) undergoes nucleocytoplasmic shuttling, but is primarily nuclear-localized and mediates expression of genes involved in development and homeostasis. Given the proximity of TR acetylation and sumoylation sites to nuclear localization (NLS) and nuclear export signals, we investigated their role in regulating intracellular localization. The nuclear/cytosolic fluorescence ratio (N/C) of fluorescent protein-tagged acetylation mimic, nonacetylation mimic, and sumoylation-deficient TR was quantified in transfected mammalian cells. While nonacetylation mimic and sumoylation-deficient TRs displayed wild-type N/C, the acetylation mimic's N/C was significantly lower. Importins that interact with wild-type TR also interact with acetylation and nonacetylation mimics, suggesting factors other than reduced importin binding alter nuclear localization. FRAP analysis showed wild-type intranuclear dynamics of acetylation mimic and sumoylation-deficient TRs, whereas the nonacetylation mimic had significantly reduced mobility and transcriptional activity. Acetyltransferase CBP/p300 inhibition enhanced TR's nuclear localization, further suggesting that nonacetylation correlates with nuclear retention, while acetylation promotes cytosolic localization.

Insight Into Molecular Determinants of T3 vs T4 Recognition From Mutations in Thyroid Hormone Receptor α and β

CONTEXT

The two major forms of circulating thyroid hormones (THs) are T3 and T4. T3 is regarded as the biologically active hormone because it binds to TH receptors (TRs) with greater affinity than T4. However, it is currently unclear what structural mechanisms underlie this difference in affinity.

OBJECTIVE

Prompted by the identification of a novel M256T mutation in a resistance to TH (RTH)α patient, we investigated Met256 in TRα1 and the corresponding residue (Met310) in TRβ1, residues previously predicted by crystallographic studies in discrimination of T3 vs T4.

METHODS

Clinical characterization of the RTHα patient and molecular studies (in silico protein modeling, radioligand binding, transactivation, and receptor-cofactor studies) were performed.

RESULTS

Structural modeling of the TRα1-M256T mutant showed that distortion of the hydrophobic niche to accommodate the outer ring of ligand was more pronounced for T3 than T4, suggesting that this substitution has little impact on the affinity for T4. In agreement with the model, TRα1-M256T selectively reduced the affinity for T3. Also, unlike other naturally occurring TRα mutations, TRα1-M256T had a differential impact on T3- vs T4-dependent transcriptional activation. TRα1-M256A and TRβ1-M310T mutants exhibited similar discordance for T3 vs T4.

CONCLUSIONS

Met256-TRα1/Met310-TRβ1 strongly potentiates the affinity of TRs for T3, thereby largely determining that T3 is the bioactive hormone rather than T4. These observations provide insight into the molecular basis for underlying the different affinity of TRs for T3 vs T4, delineating a fundamental principle of TH signaling.

Replacement of Protein Binding-Site Waters Contributes to Favorable Halogen Bond Interactions

Halogen bond interaction between a protein electronegative atom and a ligand halogen atom is increasingly attracting attention in the field of structure-based drug design. Nevertheless, gaps in understanding make it desirable to better examine the role of forces governing the formation of favorable halogen bond interactions, and the development of effective and efficient computational approaches to "design in" favorable halogen bond interactions in lead optimization process are warranted. Here, we analyzed the binding-site water properties of crystal structures with characterized halogen bond interactions between ligand halogen atoms and protein backbone carbonyl groups and, thus, found that halogen atoms involved in halogen bond interactions frequently replace calculated binding-site waters upon ligand binding. Moreover, we observed that the preferential directionality of halogen bond interactions aligns well with the orientations of these replaced waters, and these replaced waters exhibited differential energetic characteristics as compared to waters that are displaced by halogen atoms that do not form halogen bond interactions. Our discovery that replacement of calculated binding-site waters contributes to the formation of favorable halogen bond interactions suggests a practical approach for rational drug design utilizing halogen bond interactions with protein backbone carbonyl groups.

Echo dephasing and heat capacity from constrained and unconstrained dynamics of triiodothyronine nuclear receptor protein

The objective of this study is to observe the echo feature curves, vibrational dephasing, and heat capacity of a protein-hormone system taking thyroid hormone receptor-beta (THR-β) as an example. Constrained and unconstrained molecular dynamics simulations are performed by implementing the theory of velocity reassignments to probe the phase coherent state in terms of echo pulses. The constrained vibrations are incorporated by adjusting rigid bonds to all hydrogen atoms with an integrator parameter of 2 fs/step in order to reduce the degrees of freedom whereas 1 fs/step is used in the free vibrations of the atomic cluster. The nature of temperature auto-correlation functions changes so that echo feature curves also show a distinct nature in the cases of constrained and unconstrained vibrations. There is a large variation in kinetic temperature and internal potential energy in the echo time zone. The temperature rate of change of internal potential energy is the main contributor to the heat capacity of the native state protein-hormone system. The heat capacity of proteins estimated from this technique is in good agreement with the values from experiments. This study shows that triiodothyronine (T3) hormone makes some differences in heat capacity upon binding to the THR-β ligand binding domain (LBD). The physical properties of unliganded THR-β and T3-bound THR-β LBD in the cases of constrained and unconstrained dynamics are observed distinctly under the effect of anharmonicity on the phase coherent state of normal modes and the dephasing time lies in a range of 0.6-0.8 ps when the systems are perturbed suddenly.

Data on biosynthesis of BPAF glucuronide, enzyme kinetics of BPAF glucuronidation, and molecular modeling

Bisphenol AF (BPAF) is in the body mainly metabolized to the corresponding bisphenol AF glucuronide (BPAF-G). While BPAF-G is not commercially available, enzyme-assisted synthesis of BPAF-G using the human recombinant enzyme UGT2A1, purification of BPAF-G by solid phase extraction and semi-preparative HPLC and chemical characterization of BPAF-G by NMR and LC-MS/MS were performed and are described here. Furthermore, BPAF glucuronidation kinetics with the UGT enzymes that showed the highest glucuronidation activity in previous studies (i.e hepatic UGTs 1A3, 2B7, and 2B17, intestinal UGT 1A10 and UGT 2A1 that is present in airways) was performed and data is presented. Hepatic enzymes exhibited high affinities toward BPAF, while extrahepatic UGTs 2A1 and 1A10 showed the high v_max values (3.3 and 3.0 nmol/min/mg, respectively). To understand molecular interactions of BPA, BPAF and BPAF-G with ligand biding sites of several nuclear receptors, molecular modeling was performed and data on the binding modes of BPAF, BPA, and BPAF-G in the ligand-binding sites of nuclear receptors are presented.

This article is related to “Endocrine activities and adipogenic effects of bisphenol AF and its main metabolite” (Skledar et al., 2019).

Pathological Interactions Between Mutant Thyroid Hormone Receptors and Corepressors and Their Modulation by a Thyroid Hormone Analogue with Therapeutic Potential

BACKGROUND

Thyroid hormone receptors (TRs) are tightly regulated by the corepressors nuclear receptor corepressor (NCoR) and silencing mediator of retinoic acid and thyroid hormone receptors. Three conserved corepressor/NR signature box motifs (CoRNR1-3) forming the nuclear receptor interaction domain have been identified in these corepressors. Whereas TRs regulate multiple normal physiological and developmental pathways, mutations in TRs can result in endocrine diseases and be associated with cancers due to impairment of corepressor release. Three mutants that are located in helix H11 of TRs are of special interest: TRα-M388I, a mutant associated with the development of renal clear cell carcinomas (RCCCs), and TRβ-Δ430 and TRβ-Δ432, two deletion mutants causing resistance to thyroid hormone syndrome.

METHODS

Several cell-based and biophysical methods were used to measure the affinity between wild-type and mutant TRα and TRβ and all the CoRNR motifs from corepressors to quantify the effects of different thyroid hormone analogues on these interactions. This study was coupled with the measurement of interactions between wild-type and mutant TRs in the context of a heterodimer with RXR to a NCoR fragment in the presence of the same ligands. Structural insights into the binding mode of corepressors to TRs were assessed in parallel by nuclear magnetic resonance spectroscopy.

RESULTS

The study shows that TRs interact more avidly with the silencing mediator of retinoic acid and thyroid hormone receptors than with NCoR peptides, and that TRα binds most avidly to S-CoRNR3, whereas TRβ binds preferentially to S-CoRNR2. In the studied TR mutants, a transfer of the CoRNR-specificity toward CoRNR1 was observed, coupled with a significant increase in the binding strength. In contrast to 3,5,3'-triiodothyronine (T3), the agonist TRIAC and the antagonist NH-3 were very efficient at dissociating the abnormally strong interactions between mutant TRβs and corepressors. A strong impairment of T3-binding for TRβ mutants was shown compared to TRIAC and NH-3 and could explain the different efficiencies of the different ligands in releasing corepressors from the studied TRβ mutants. Consequently, TRIAC was found to be more effective than T3 in facilitating coactivator recruitment and decreasing the dominant activity of TRβ-Δ430.

CONCLUSION

This study helps to clarify the specific interaction surfaces involved in the pathologic phenotype of TR mutants and demonstrates that TRIAC is a potential therapeutic agent for patients suffering from resistance to thyroid hormone syndromes.

Role of Leucine 341 in Thyroid Hormone Receptor Beta Revealed by a Novel Mutation Causing Thyroid Hormone Resistance

Leucine 341 has been predicted from crystal structure as an important residue for thyroid hormone receptor beta (TRβ) function, but this has never been confirmed in functional studies. Here, a novel p.L341V mutation as a cause of resistance to TRβ is described, suggesting an important role for L341 in TRβ function. In silico and in vitro studies confirmed that substituting L341 with valine and other non-polar amino acids impairs sensitivity of TRβ for triiodothyronine to various degrees, depending on their side-chain size and orientation.

Single Nucleotide Variations of the Human GR Gene Manifested as Pathologic Mutations or Polymorphisms

The human genome contains numerous single nucleotide variations, and the human glucocorticoid receptor (GR) gene harbors ∼450 of these genetic changes. Among them, extremely rare, nonsynonymous variants, known as pathologic GR gene mutations, develop a characteristic pathologic condition, familial/sporadic generalized glucocorticoid resistance syndrome, by replacing the amino acids critical for GR protein structure and functions, whereas others, known as pathologic polymorphisms, develop mild manifestations recognized mainly at population bases by changing the GR activities slightly. Recent progress on the structural analysis to the GR protein and subsequent computer-based structural simulation revealed details of the molecular defects caused by such pathologic GR gene mutations, including their impact on the receptor interaction to ligands, nuclear receptor coactivators (NCoAs) or DNA glucocorticoid response elements (GREs). Indeed, those found in the GR ligand-binding domain significantly damage protein structure of the ligand-binding pocket and/or the activation function-2 transactivation domain and change their molecular interaction to glucocorticoids or the LxxLL signature motif of NCoAs. Two mutations found in GR DNA-binding domain also affect interaction of the mutant receptors to GRE DNA by affecting the critical amino acid for the interaction or changing local hydrophobic circumstance. In this review, I discuss recent findings on the structural simulation of the pathologic GR mutants in connection to their functional and clinical impacts, along with a brief explanation to recent research achievement on the GR polymorphisms.

The Affinity of Brominated Phenolic Compounds for Human and Zebrafish Thyroid Receptor β: Influence of Chemical Structure

Brominated phenolic compounds (BPCs) are found in the environment, and in human and wildlife tissues, and some are considered to have endocrine disrupting activities. The goal of this study was to determine how structural differences of 3 BPC classes impact binding affinities for the thyroid receptor beta (TRβ) in humans and zebrafish. BPC classes included halogenated bisphenol A derivatives, halogenated oxidative transformation products of 2,2',4,4'-tetrabromodiphenyl ether (BDE-47), and brominated phenols. Affinities were assessed using recombinant TRβ protein in competitive binding assays with 125I-triiodothyronine (125I-T3) as the radioligand. Zebrafish and human TRβ displayed similar binding affinities for T3 (Ki = 0.40 and 0.49 nM) and thyroxine (T4, Ki = 6.7 and 6.8 nM). TRβ affinity increased with increasing halogen mass and atomic radius for both species, with the iodinated compounds having the highest affinity within their compound classes. Increasing halogen mass and radius increases the molecular weight, volume, and hydrophobicity of a compound, which are all highly correlated with increasing affinity. TRβ affinity also increased with the degree of halogenation for both species. Human TRβ displayed higher binding affinities for the halogenate bisphenol A compounds, whereas zebrafish TRβ displayed higher affinities for 2,4,6-trichlorophenol and 2,4,6-trifluorophenol. Observed species differences may be related to amino acid differences within the ligand binding domains. Overall, structural variations impact TRβ affinities in a similar manner, supporting the use of zebrafish as a model for TRβ disruption. Further studies are necessary to investigate how the identified structural modifications impact downstream receptor activities and potential in vivo effects.

Thyroid hormone receptor localization in target tissues

The thyroid hormone receptors, TRα1, TRβ1 and other subtypes, are members of the nuclear receptor superfamily that mediate the action of thyroid hormone signaling in numerous tissues to regulate important physiological and developmental processes. Their most well-characterized role is as ligand-dependent transcription factors; TRs bind thyroid hormone response elements in the presence or absence of thyroid hormone to facilitate the expression of target genes. Although primarily residing in the nucleus, TRα1 and TRβ1 shuttle rapidly between the nucleus and cytoplasm. We have identified multiple nuclear localization signals and nuclear export signals within TRα1 and TRβ1 that interact with importins and exportins, respectively, to mediate translocation across the nuclear envelope. More recently, enigmatic cytoplasmic functions have been ascribed to other TR subtypes, expanding the diversity of the cellular response to thyroid hormone. By integrating data on localization signal motifs, this review provides an overview of the complex interplay between TR's dynamic transport pathways and thyroid hormone signaling activities. We examine the variation in TR subtype response to thyroid hormone signaling, and what is currently known about regulation of the variety of tissue-specific localization patterns, including targeting to the nucleus, the mitochondria and the inner surface of the plasma membrane.

Protein Disulfide Isomerase Modulates the Activation of Thyroid Hormone Receptors

Thyroid hormone receptors (TRs) are responsible for mediating thyroid hormone (T3 and T4) actions at a cellular level. They belong to the nuclear receptor (NR) superfamily and execute their main functions inside the cell nuclei as hormone-regulated transcription factors. These receptors also exhibit so-called "non-classic" actions, for which other cellular proteins, apart from coregulators inside nuclei, regulate their activity. Aiming to find alternative pathways of TR modulation, we searched for interacting proteins and found that PDIA1 interacts with TRβ in a yeast two-hybrid screening assay. The functional implications of PDIA1-TR interactions are still unclear; however, our co-immunoprecipitation (co-IP) and fluorescence assay results showed that PDI was able to bind both TR isoforms in vitro. Moreover, T3 appears to have no important role in these interactions in cellular assays, where PDIA1 was able to regulate transcription of TRα and TRβ-mediated genes in different ways depending on the promoter region and on the TR isoform involved. Although PDIA1 appears to act as a coregulator, it binds to a TR surface that does not interfere with coactivator binding. However, the TR:PDIA1 complex affinity and activation are different depending on the TR isoform. Such differences may reflect the structural organization of the PDIA1:TR complex, as shown by models depicting an interaction interface with exposed cysteines from both proteins, suggesting that PDIA1 might modulate TR by its thiol reductase/isomerase activity.

Novel thyroid hormone analogues, enzyme inhibitors and mimetics, and their action

Thyroid hormones (THs) play key roles in modulating the overall metabolism of the body, protein synthesis, fat metabolism, neuronal and bone growth, and cardiovascular as well as renal functions. In this review, we discuss on the thyroid hormone synthesis and activation, thyroid hormone receptors (TRs) and mechanism of action, applications of thyroid hormone analogues, particularly the compounds that are selective ligands for TRβ receptors, or enzyme inhibitors for the treatment of thyroidal disorders with a specific focus on thyroid peroxidase and iodothyronine deiodinases. We also discuss on the development of small-molecule deiodinase mimetics and their mechanism of deiodination, as these compounds have the potential to regulate the thyroid hormone levels.

Structural aspects of thyroid hormone binding to proteins and competitive interactions with natural and synthetic compounds

Thyroid hormones and their metabolites constitute a vast class of related iodothyronine compounds that contribute to the regulation of metabolic activity and cell differentiation. They are in turn transported, transformed and recognized as signaling molecules through binding to a variety of proteins from a wide range of evolutionary unrelated protein families, which renders these proteins and their iodothyronine binding sites an example for extensive convergent evolution. In this review, we will briefly summarize what is known about iodothyronine binding sites in proteins, the modes of protein/iodothyronine interaction, and the ligand conformations. We will then discuss physiological and synthetic compounds, including popular drugs and food components, that can interfere with iodothyronine binding and recognition by these proteins. The discussion also includes compounds persisting in the environment and acting as endocrine disrupting chemicals.

Naturally Occurring Amino Acids in Helix 10 of the Thyroid Hormone Receptor Mediate Isoform-Specific TH Gene Regulation

Thyroid hormone (TH) action is mediated by the products of two genes, TH receptor (THR)α (THRA) and THRβ (THRB) that encode several closely related receptor isoforms with differing tissue distributions. The vast majority of THR isoform-specific effects are thought to be due to tissue-specific differences in THR isoform expression levels. We investigated the alternative hypothesis that intrinsic functional differences among THR isoforms mediate these tissue-specific effects. To achieve the same level of expression of each isoform, we created tagged THR isoforms and tested their DNA and functional properties in vitro. We found significant homodimerization and functional differences among the THR isoforms. THRA1 was unable to form homodimers on direct repeat separated by 4 bp DNA elements and was also defective in TH-dependent repression of Tshb and Rxrg in a thyrotroph cell line, TαT1.1. In contrast, THRB2 was both homodimer sufficient and fully functional on these negatively regulated genes. Using domain exchanges and individual amino acid switches between THRA1 and THRB2, we identified three amino acids in helix 10 of the THRB2 ligand-binding domain that are required for negative regulation and are absent in THRA1.

Nuclear Import and Export of the Thyroid Hormone Receptor

The thyroid hormone receptors, TRα1 and TRβ1, are members of the nuclear receptor superfamily that forms one of the most abundant classes of transcription factors in multicellular organisms. Although primarily localized to the nucleus, TRα1 and TRβ1 shuttle rapidly between the nucleus and cytoplasm. The fine balance between nuclear import and export of TRs has emerged as a critical control point for modulating thyroid hormone-responsive gene expression. Mutagenesis studies have defined two nuclear localization signal (NLS) motifs that direct nuclear import of TRα1: NLS-1 in the hinge domain and NLS-2 in the N-terminal A/B domain. Three nuclear export signal (NES) motifs reside in the ligand-binding domain. A combined approach of shRNA-mediated knockdown and coimmunoprecipitation assays revealed that nuclear entry of TRα1 is facilitated by importin 7, likely through interactions with NLS-2, and importin β1 and the adapter importin α1 interacting with both NLS-1 and NLS-2. Interestingly, TRβ1 lacks NLS-2 and nuclear import depends solely on the importin α1/β1 heterodimer. Heterokaryon and fluorescence recovery after photobleaching shuttling assays identified multiple exportins that play a role in nuclear export of TRα1, including CRM1 (exportin 1), and exportins 4, 5, and 7. Even single amino acid changes in TRs dramatically alter their intracellular distribution patterns. We conclude that mutations within NLS and NES motifs affect nuclear shuttling activity, and propose that TR mislocalization contributes to the development of some types of cancer and Resistance to Thyroid Hormone syndrome.

Membrane-traversing mechanism of thyroid hormone transport by monocarboxylate transporter 8

Monocarboxylate transporter 8 (MCT8) mediates thyroid hormone (TH) transport across the plasma membrane in many cell types. In order to better understand its mechanism, we have generated three new MCT8 homology models based on sugar transporters XylE in the intracellular opened (PDB ID: 4aj4) and the extracellular partly occluded (PDB ID: 4gby) conformations as well as FucP (PDB ID: 3o7q) and GLUT3 (PDB ID: 4zwc) in the fully extracellular opened conformation. T3-docking studies from both sides revealed interactions with His192, His415, Arg445 and Asp498 as previously identified. Selected mutations revealed further transport-sensitive positions mainly at the discontinuous transmembrane helices TMH7 and 10. Lys418 is potentially involved in neutralising the charge of the TH substrate because it can be replaced by charged, but not by uncharged, amino acids. The side chain of Thr503 was hypothesised to stabilise a helix break at TMH10 that undergoes a prominent local shift during the transport cycle. A T503V mutation accordingly affected transport. The aromatic Tyr419, the polar Ser313 and Ser314 as well as the charged Glu422 and Glu423 lining the transport channel have been studied. Based on related sugar transporters, we suggest an alternating access mechanism for MCT8 involving a series of amino acid positions previously and newly identified as critical for transport.

Human EGF-derived direct and reverse short linear motifs: conformational dynamics insight into the receptor-binding residues

Short linear motifs (SLiMs) have been recognized to perform diverse functions in a variety of regulatory proteins through the involvement in protein-protein interactions, signal transduction, cell cycle regulation, protein secretion, etc. However, detailed molecular mechanisms underlying their functions including roles of definite amino acid residues remain obscure. In our previous studies, we demonstrated that conformational dynamics of amino acid residues in oligopeptides derived from regulatory proteins such as alpha-fetoprotein (AFP), carcino-embryonic antigen (CEA), and pregnancy specific β1-glycoproteins (PSGs) contributes greatly to their biological activities. In the present work, we revealed the 22-member linear modules composed of direct and reverse AFP_14-20-like heptapeptide motifs linked by CxxGY/FxGx consensus motif within epidermal growth factor (EGF), growth factors of EGF family and numerous regulatory proteins containing EGF-like modules. We showed, first, the existence of similarity in amino acid signatures of both direct and reverse motifs in terms of their physicochemical properties. Second, molecular dynamics (MD) simulation study demonstrated that key receptor-binding residues in human EGF in the aligned positions of the direct and reverse motifs may have similar distribution of conformational probability densities and dynamic behavior despite their distinct physicochemical properties. Third, we found that the length of a polypeptide chain (from 7 to 53 residues) has no effect, while disulfide bridging and backbone direction significantly influence the conformational distribution and dynamics of the residues. Our data may contribute to the atomic level structure-function analysis and protein structure decoding; additionally, they may provide a basis for novel protein/peptide engineering and peptide-mimetic drug design.

In vitro assessment of thyroid hormone receptor activity of four organophosphate esters

Previous animal experiments have implied that organophosphate esters (OPEs) have a disruption effect on the thyroid endocrine system. However, knowledge of the toxicological mechanism remains limited. In this study, the activities of four OPEs have been characterized against the thyroid hormone (TH) nuclear receptor (TR) using two in vitro models, with the aim of evaluating their toxicity mechanisms towards the TR. The results of a TH-dependent cell proliferation assay showed that tris(2-chloro-1-(chloromethyl)ethyl)phosphate (TDCPP) could induce cell growth, while the other three OPEs had no effect. The results of a luciferase reporter gene assay revealed that all four of the OPEs tested in the current study showed agonistic activity towards TRβ, with TDCPP being the most potent one. Moreover, molecular docking revealed that all the tested OPEs could fit into the ligand binding pocket of TRβ, with TDCPP binding more effectively than the other three OPEs. Taken together, these data suggest that OPEs might disrupt the thyroid endocrine system via a mechanism involving the activation of TR.

Experimental Binding Energies in Supramolecular Complexes

On the basis of many literature measurements, a critical overview is given on essential noncovalent interactions in synthetic supramolecular complexes, accompanied by analyses with selected proteins. The methods, which can be applied to derive binding increments for single noncovalent interactions, start with the evaluation of consistency and additivity with a sufficiently large number of different host-guest complexes by applying linear free energy relations. Other strategies involve the use of double mutant cycles, of molecular balances, of dynamic combinatorial libraries, and of crystal structures. Promises and limitations of these strategies are discussed. Most of the analyses stem from solution studies, but a few also from gas phase. The empirically derived interactions are then presented on the basis of selected complexes with respect to ion pairing, hydrogen bonding, electrostatic contributions, halogen bonding, π-π-stacking, dispersive forces, cation-π and anion-π interactions, and contributions from the hydrophobic effect. Cooperativity in host-guest complexes as well as in self-assembly, and entropy factors are briefly highlighted. Tables with typical values for single noncovalent free energies and polarity parameters are in the Supporting Information.

New insights into the structure and mechanism of iodothyronine deiodinases

Iodothyronine deiodinases are a family of enzymes that remove specific iodine atoms from one of the two aromatic rings in thyroid hormones (THs). They thereby fine-tune local TH concentrations and cellular TH signaling. Deiodinases catalyze a remarkable biochemical reaction, i.e., the reductive elimination of a halogenide from an aromatic ring. In metazoans, deiodinases depend on the rare amino acid selenocysteine. The recent solution of the first experimental structure of a deiodinase catalytic domain allowed for a reappraisal of the many mechanistic and mutagenesis data that had been accumulated over more than 30 years. Hence, the structure generates new impetus for research directed at understanding catalytic mechanism, substrate specificity, and regulation of deiodinases. This review will focus on structural and mechanistic aspects of iodothyronine deiodinases and briefly compare these enzymes with dehalogenases, which catalyze related reactions. A general mechanism for the selenium-dependent deiodinase reaction will be described, which integrates the mouse deiodinase 3 crystal structure and biochemical studies. We will summarize further, sometimes isoform-specific molecular features of deiodinase catalysis and regulation, and we will then discuss available compounds for modulating deiodinase activity for therapeutic purposes.

Androgen receptor: structure, role in prostate cancer and drug discovery

Androgens and androgen receptors (AR) play a pivotal role in expression of the male phenotype. Several diseases, such as androgen insensitivity syndrome (AIS) and prostate cancer, are associated with alterations in AR functions. Indeed, androgen blockade by drugs that prevent the production of androgens and/or block the action of the AR inhibits prostate cancer growth. However, resistance to these drugs often occurs after 2–3 years as the patients develop castration-resistant prostate cancer (CRPC). In CRPC, a functional AR remains a key regulator. Early studies focused on the functional domains of the AR and its crucial role in the pathology. The elucidation of the structures of the AR DNA binding domain (DBD) and ligand binding domain (LBD) provides a new framework for understanding the functions of this receptor and leads to the development of rational drug design for the treatment of prostate cancer. An overview of androgen receptor structure and activity, its actions in prostate cancer, and how structural information and high-throughput screening have been or can be used for drug discovery are provided herein.