In silico insights into the dimer structure and deiodinase activity of type III iodothyronine deiodinase from bioinformatics, molecular dynamics simulations, and QM/MM calculations.

The homodimeric family of iodothyronine deiodinases (Dios) regioselectively remove iodine from thyroid hormones. Currently, structural data has only been reported for the monomer of the mus type III thioredoxin (Trx) fold catalytic domain (Dio3Trx), but the mode of dimerization has not yet been determined. Various groups have proposed dimer structures that are similar to the A-type and B-type dimerization modes of peroxiredoxins. Computational methods are used to compare the sequence of Dio3Trx to related proteins known to form A-type and B-type dimers. Sequence analysis and in silico protein-protein docking methods suggest that Dio3Trx is more consistent with proteins that adopt B-type dimerization. Molecular dynamics (MD) simulations of the refined Dio3Trx dimer constructed using the SymmDock and GalaxyRefineComplex databases indicate stable dimer formation along the ß4α3 interface consistent with other Trx fold B-type dimers. Free energy calculations show that the dimer is stabilized by interdimer interactions between the ß-sheets and α-helices. A comparison of MD simulations of the apo and thyroxine-bound dimers suggests that the active site binding pocket is not affected by dimerization. Determination of the transition state for deiodination of thyroxine from the monomer structure using QM/MM methods provides an activation barrier consistent with previous small model DFT studies.Communicated by Ramaswamy H. Sarma.

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 (T4) through an I⋯Se halogen bond to Sec170 and the amino acid group with a polar cleft. The Dio3-T4 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.

A Halogen Bonding Perspective on Iodothyronine Deiodinase Activity.

Iodothyronine deiodinases (Dios) are involved in the regioselective removal of iodine from thyroid hormones (THs). Deiodination is essential to maintain TH homeostasis, and disruption can have detrimental effects. Halogen bonding (XB) to the selenium of the selenocysteine (Sec) residue in the Dio active site has been proposed to contribute to the mechanism for iodine removal. Polybrominated diphenyl ethers (PBDEs) and polychlorinated biphenyls (PCBs) are known disruptors of various pathways of the endocrine system. Experimental evidence shows PBDEs and their hydroxylated metabolites (OH-BDEs) can inhibit Dio, while data regarding PCB inhibition are limited. These xenobiotics could inhibit Dio activity by competitively binding to the active site Sec through XB to prevent deiodination. XB interactions calculated using density functional theory (DFT) of THs, PBDEs, and PCBs to a methyl selenolate (MeSe-) arrange XB strengths in the order THs > PBDEs > PCBs in agreement with known XB trends. THs have the lowest energy C-X*-type unoccupied orbitals and overlap with the Se lp donor leads to high donor-acceptor energies and the greatest activation of the C-X bond. The higher energy C-Br* and C-Cl* orbitals similarly result in weaker donor-acceptor complexes and less activation of the C-X bond. Comparison of the I···Se interactions for the TH group suggest that a threshold XB strength may be required for dehalogenation. Only highly brominated PBDEs have binding energies in the same range as THs, suggesting that these compounds may inhibit Dio and undergo debromination. While these small models provide insight on the I···Se XB interaction itself, interactions with other active site residues are governed by regioselective preferences observed in Dios.

Halogen Bonding Interactions of Polychlorinated Biphenyls and the Potential for Thyroid Disruption.

Polychlorinated biphenyl (PCB) flame retardants are persistent pollutants and inhibit neurodevelopment, particularly in the early stages of life. Halogen bonding (XB) to the iodothyronine deiodinases (Dio) that modulate thyroid hormones (THs) is a potential mechanism for endocrine disruption. Cl⋅⋅⋅Se XB interactions of PCBs with SeMe- , a small model of the Dio active site selenocysteine, are compared with previous results on polybrominated diphenylethers (PBDEs) and THs using density functional theory. PCBs generally display weaker XB interactions compared to PBDEs and THs, consistent with the dependence of XB strength on the size of the halogen (I>Br>Cl). PCBs also do not meet a proposed energy threshold for substrates to undergo dehalogenation, suggesting they may behave as competitive inhibitors of Dio in addition to other mechanisms of endocrine disruption. XB interactions in PCBs are position-dependent, with ortho interactions slightly more favorable than meta and para interactions, suggesting that PCBs may have a greater effect on certain classes of Dio. Flexibility of PCBs around the biphenyl C-C bond is limited by ortho substitutions relative to the biphenyl linkage, which may contribute to the ability to inhibit Dio and other TH-related proteins.

Halogen-Bonding Interactions of Polybrominated Diphenyl Ethers and Thyroid Hormone Derivatives: A Potential Mechanism for the Inhibition of Iodothyronine Deiodinase.

Polybrominated diphenyl ethers (PBDEs) increase flame resistance in consumer goods, but these compounds and their hydroxylated derivatives (OH-BDEs) impair normal thyroid function. Halogen bonding (XB) of PBDEs to an active site selenocysteine may prevent iodothyronine deiodinase(Dio)-catalyzed activation/deactivation of thyroid hormone (TH) derivatives. In this study, we compare the strength of the XB interactions of TH derivatives, iodine-based contrast agents and PBDEs/OH-BDEs with a methylselenolate model of the Dio active site using density functional theory calculations. The strength of the XB interaction depends upon the acceptor halide, the position of the halide, the number of ring substituents, and the proximity of hydroxyl groups to the XB site. The weaker Se⋅⋅⋅Br interactions relative to Se⋅⋅⋅I interactions are consistent with a model of competitive inhibition that blocks binding of THs at elevated PBDE/OH-BDE concentrations. XB interactions were generally more favorable at ortho and meta positions and in substrates with more electron-withdrawing substituents. PBDEs/OH-BDEs that mimic the binding behavior of THs, that is, containing ortho and meta bromides and adjacent hydroxyl groups, may be the most effective inhibitors. Highly-brominated PBDEs/OH-BDEs have comparable interaction energies to THs and may undergo debromination. These results may also suggest that XB strength must exceed a threshold value in order for PBDEs/OH-BDEs to undergo nucleophilic attack by Dio.

Structure of 'linkerless' hydroxamic acid inhibitor-HDAC8 complex confirms the formation of an isoform-specific subpocket.

Histone deacetylases (HDACs) catalyze the hydrolysis of acetylated lysine side chains in histone and non-histone proteins, and play a critical role in the regulation of many biological processes, including cell differentiation, proliferation, senescence, and apoptosis. Aberrant HDAC activity is associated with cancer, making these enzymes important targets for drug design. In general, HDAC inhibitors (HDACi) block the proliferation of tumor cells by inducing cell differentiation, cell cycle arrest, and/or apoptosis, and comprise some of the leading therapies in cancer treatments. To date, four HDACi have been FDA approved for the treatment of cancers: suberoylanilide hydroxamic acid (SAHA, Vorinostat, Zolinza®), romidepsin (FK228, Istodax®), belinostat (Beleodaq®), and panobinostat (Farydak®). Most current inhibitors are pan-HDACi, and non-selectively target a number of HDAC isoforms. Six previously reported HDACi were rationally designed, however, to target a unique sub-pocket found only in HDAC8. While these inhibitors were indeed potent against HDAC8, and even demonstrated specificity for HDAC8 over HDACs 1 and 6, there were no structural data to confirm the mode of binding. Here we report the X-ray crystal structure of Compound 6 complexed with HDAC8 to 1.98Å resolution. We also describe the use of molecular docking studies to explore the binding interactions of the other 5 related HDACi. Our studies confirm that the HDACi induce the formation of and bind in the HDAC8-specific subpocket, offering insights into isoform-specific inhibition.