Molecular dynamics simulations reveal multiple pathways of ligand dissociation from thyroid hormone receptors

Nuclear receptor (NR) ligands occupy a pocket that lies within the core of the NR ligand-binding domain (LBD), and most NR LBDs lack obvious entry/exit routes upon the protein surface. Thus, significant NR conformational rearrangements must accompany ligand binding and release. The precise nature of these processes, however, remains poorly understood. Here, we utilize locally enhanced sampling (LES) molecular dynamics computer simulations to predict molecular motions of x-ray structures of thyroid hormone receptor (TR) LBDs and determine events that permit ligand escape. We find that the natural ligand 3,5,3'-triiodo-L-thyronine (T(3)) dissociates from the TRalpha1 LBD along three competing pathways generated through i), opening of helix (H) 12; ii), separation of H8 and H11 and the Omega-loop between H2 and H3; and iii), opening of H2 and H3, and the intervening beta-strand. Similar pathways are involved in dissociation of T(3) and the TRbeta-selective ligand GC24 from TRbeta; the TR agonist IH5 from the alpha- and beta-TR forms; and Triac from two natural human TRbeta mutants, A317T and A234T, but are detected with different frequencies in simulations performed with the different structures. Path I was previously suggested to represent a major pathway for NR ligand dissociation. We propose here that Paths II and III are also likely ligand escape routes for TRs and other NRs. We also propose that different escape paths are preferred in different situations, implying that it will be possible to design NR ligands that only associate stably with their cognate receptors in specific cellular contexts.     

Conformational adaptation of nuclear receptor ligand binding domains to agonists: potential for novel approaches to ligand design

Ligands occupy the core of nuclear receptor (NR) ligand binding domains (LBDs) and modulate NR function. X-ray structures of NR LBDs reveal most NR agonists fill the enclosed pocket and promote packing of C-terminal helix 12 (H12), whereas the pockets of unliganded NR LBDs differ. Here, we review evidence that NR pockets rearrange to accommodate different agonists. Some thyroid hormone receptor (TR) ligands with 5′ extensions designed to perturb H12 act as antagonists, but many are agonists. One mode of adaptation is seen in a TR/thyroxine complex; the pocket expands to accommodate a 5′ iodine extension. Crystals of other NR LBDs reveal that the pocket can expand or contract and some agonists do not fill the pocket. A TRβ structure in complex with an isoform selective drug (GC-24) reveals another mode of adaptation; the LBD hydrophobic interior opens to accommodate a bulky 3′ benzyl extension. We suggest that placement of extensions on NR agonists will highlight unexpected areas of flexibility within LBDs that could accommodate extensions; thereby enhancing the selectivity of agonist binding to particular NRs. Finally, agonists that induce similar LBD structures differ in their activities and we discuss strategies to reveal subtle structural differences responsible for these effects.     

Molecular basis for dimer formation of TRbeta variant D355R

Protein quality and stability are critical during protein purification for X-ray crystallography. A target protein that is easy to manipulate and crystallize becomes a valuable product useful for high-throughput crystallography for drug design and discovery. In this work, a single surface mutation, D355R, was shown to be crucial for converting the modestly stable monomeric ligand binding domain of the human thyroid hormone receptor (TR LBD) into a stable dimer. The structure of D335R TR LBD mutant was solved using X-ray crystallography and refined to 2.2 A resolution with R(free)/R values of 24.5/21.7. The crystal asymmetric unit reveals the TR dimer with two molecules of the hormone-bound LBD related by twofold symmetry. The ionic interface between the two LBDs comprises residues within loop H10-H11 and loop H6-H7 as well as the C-terminal halves of helices 8 of both protomers. Direct intermolecular contacts formed between the introduced residue Arg 355 of one TR molecule and Glu 324 of the second molecule become a part of the extended dimerization interface of 1330 A(2) characteristic for a strong complex assembly that is additionally strengthened by buffer solutes.     

Hormone selectivity in thyroid hormone receptors

Separate genes encode thyroid hormone receptor subtypes TRalpha (NR1A1) and TRbeta (NR1A2). Products from each of these contribute to hormone action, but the subtypes differ in tissue distribution and physiological response. Compounds that discriminate between these subtypes in vivo may be useful in treating important medical problems such as obesity and hypercholesterolemia. We previously determined the crystal structure of the rat (r) TRalpha ligand-binding domain (LBD). In the present study, we determined the crystal structure of the rTRalpha LBD in a complex with an additional ligand, Triac (3,5, 3'-triiodothyroacetic acid), and two crystal structures of the human (h) TRbeta receptor LBD in a complex with either Triac or a TRbeta-selective compound, GC-1 [3,5-dimethyl-4-(4'-hydroy-3'-isopropylbenzyl)-phenoxy acetic acid]. The rTRalpha and hTRbeta LBDs show close structural similarity. However, the hTRbeta structures extend into the DNA-binding domain and allow definition of a structural "hinge" region of only three amino acids. The two TR subtypes differ in the loop between helices 1 and 3, which could affect both ligand recognition and the effects of ligand in binding coactivators and corepressors. The two subtypes also differ in a single amino acid residue in the hormone-binding pocket, Asn (TRbeta) for Ser (TRalpha). Studies here with TRs in which the subtype-specific residue is exchanged suggest that most of the selectivity in binding derives from this amino acid difference. The flexibility of the polar region in the TRbeta receptor, combined with differential recognition of the chemical group at the 1-carbon position, seems to stabilize the complex with GC-1 and contribute to its beta-selectivity. These results suggest a strategy for development of subtype-specific compounds involving modifications of the ligand at the 1-position.     

Molecular modeling, affinity labeling, and site-directed mutagenesis define the key points of interaction between the ligand-binding domain of the vitamin D nuclear receptor and 1 alpha,25-dihydroxyvitamin D3

We have combined molecular modeling and classical structure-function techniques to define the interactions between the ligand-binding domain (LBD) of the vitamin D nuclear receptor (VDR) and its natural ligand, 1alpha,25-dihydroxyvitamin D(3) [1alpha,25-(OH)(2)D(3)]. The affinity analogue 1alpha,25-(OH)(2)D(3)-3-bromoacetate exclusively labeled Cys-288 in the VDR-LBD. Mutation of C288 to glycine abolished this affinity labeling, whereas the VDR-LBD mutants C337G and C369G (other conserved cysteines in the VDR-LBD) were labeled similarly to the wild-type protein. These results revealed that the A-ring 3-OH group docks next to C288 in the binding pocket. We further mutated M284 and W286 (separately creating M284A, M284S, W286A, and W286F) and caused severe loss of ligand binding, indicating the crucial role played by the contiguous segment between M284 and C288. Alignment of the VDR-LBD sequence with the sequences of nuclear receptor LBDs of known 3-D structure positioned M284 and W286 in the presumed beta-hairpin of the molecule, thereby identifying it as the region contacting the A-ring of 1alpha, 25-(OH)(2)D(3). From the multiple sequence alignment, we developed a homologous extension model of the VDR-LBD. The model has a canonical nuclear receptor fold with helices H1-H12 and a single beta hairpin but lacks the long insert (residues 161-221) between H2 and H3. We docked the alpha-conformation of the A-ring into the binding pocket first so as to incorporate the above-noted interacting residues. The model predicts hydrogen bonding contacts between ligand and protein at S237 and D299 as well as at the site of the natural mutation R274L. Mutation of S237 or D299 to alanine largely abolished ligand binding, whereas changing K302, a nonligand-contacting residue, to alanine left binding unaffected. In the "activation" helix 12, the model places V418 closest to the ligand, and, consistent with this prediction, the mutation V418S abolished ligand binding. The studies together have enabled us to identify 1alpha,25-(OH)(2)D(3)-binding motifs in the ligand-binding pocket of VDR.     

Design of thyroid hormone receptor antagonists from first principles

It is desirable to obtain TR antagonists for treatment of hyperthyroidism and other conditions. We have designed TR antagonists from first principles based on TR crystal structures. Since agonist ligands are buried in the fold of the TR ligand binding domain (LBD), we reasoned that ligands that resemble agonists with large extensions should bind the LBD, but would prevent its folding into an active conformation. In particular, we predicted that extensions at the 5′ aryl position of ligand should reposition helix (H) 12, which forms part of the co-activator binding surface, and thereby inhibit TR activity. We have found that some synthetic ligands with 5′ aryl ring extensions behave as antagonists (DIBRT, NH-3), or partial antagonists (GC-14, NH-4). Moreover, one compound (NH-3) represents the first potent TR antagonist with nanomolar affinity that also inhibits TR action in an animal model. However, the properties of the ligands also reveal unexpected aspects of TR behavior. While nuclear receptor antagonists generally promote binding of co-repressors, NH-3 blocks co-activator binding and also prevents co-repressor binding. More surprisingly, many compounds with extensions behave as full or partial agonists. We present hypotheses to explain both behaviors in terms of dynamic equilibrium of H12 position.     

Low-resolution structures of thyroid hormone receptor dimers and tetramers in solution

High-resolution X-ray structures of thyroid hormone (TH) receptor (TR) DNA and ligand binding domains (DBD and LBD) have yielded significant insights into TR action. Nevertheless, the TR DBD and LBD act in concert to mediate TH effects upon gene expression, and TRs form multiple oligomers; however, structures of full-length TRs or DBD-LBD constructs that would clarify these influences are not available. Here, we report low-resolution X-ray structures of the TRbeta DBD-LBD construct in solution which define the shape of dimers and tetramers and likely positions of the DBDs and LBDs. The holo TRbeta DBD-LBD construct forms a homodimer with LBD-DBD pairs in close contact and DBDs protruding from the base in the same direction. The DBDs are connected to the LBDs by crossed extended D domains. The apo hTRbeta DBD-LBD construct forms tetramers that resemble bulged cylinders with pairs of LBD dimers in a head-to-head arrangement with DBD pairs packed tightly against the LBD core. Overall, there are similarities with our previous low-resolution structures of retinoid X receptors, but TRs exhibit two unique features. First, TR DBDs are closely juxtaposed in the dimer and tetramer forms. Second, TR DBDs are closely packed against LBDs in the tetramer, but not the dimer. These findings suggest that TRs may be able to engage in hitherto unknown interdomain interactions and that the D domain must rearrange in different oligomeric forms. Finally, the data corroborate our suggestion that apo TRs form tetramers in solution which dissociate into dimers upon hormone binding.     

Medium chain fatty acids are selective peroxisome proliferator activated receptor (PPAR) γ activators and pan-PPAR partial agonists

Thiazolidinediones (TZDs) act through peroxisome proliferator activated receptor (PPAR) γ to increase insulin sensitivity in type 2 diabetes (T2DM), but deleterious effects of these ligands mean that selective modulators with improved clinical profiles are needed. We obtained a crystal structure of PPARγ ligand binding domain (LBD) and found that the ligand binding pocket (LBP) is occupied by bacterial medium chain fatty acids (MCFAs). We verified that MCFAs (C8-C10) bind the PPARγ LBD in vitro and showed that they are low-potency partial agonists that display assay-specific actions relative to TZDs; they act as very weak partial agonists in transfections with PPARγ LBD, stronger partial agonists with full length PPARγ and exhibit full blockade of PPARγ phosphorylation by cyclin-dependent kinase 5 (cdk5), linked to reversal of adipose tissue insulin resistance. MCFAs that bind PPARγ also antagonize TZD-dependent adipogenesis in vitro. X-ray structure B-factor analysis and molecular dynamics (MD) simulations suggest that MCFAs weakly stabilize C-terminal activation helix (H) 12 relative to TZDs and this effect is highly dependent on chain length. By contrast, MCFAs preferentially stabilize the H2-H3/β-sheet region and the helix (H) 11-H12 loop relative to TZDs and we propose that MCFA assay-specific actions are linked to their unique binding mode and suggest that it may be possible to identify selective PPARγ modulators with useful clinical profiles among natural products.     

Probing conformational changes in the estrogen receptor: evidence for a partially unfolded intermediate facilitating ligand binding and release

Because the ligand bound to the ligand-binding domain (LBD) of nuclear hormone receptors is completely enveloped by protein, it is thought that the process of ligand binding or unbinding must involve a significant conformational change of this domain. We have used the intrinsic tryptophan fluorescence of the estrogen receptor-alpha (ERalpha) or estrogen receptor-beta (ERbeta) LBD, as well as bis-anilinonaphthalenesulfonate (bis-ANS), a probe for accessible interior regions of protein, to follow the guanidine-hydrochloride (Gua-HCl)-induced unfolding of this domain. In both cases, we find that the ER-LBD unfolding follows a two-phase process. At low Gua-HCl, the ER-LBD undergoes partial unfolding, whereas at high Gua-HCl, this domain undergoes a global unfolding, with bis-ANS binding preferentially to the partially unfolded state. The partially unfolded state of the ERalpha-LBD induced by denaturant does not bind ligand stably, but it may resemble an intermediate that this domain accesses transiently under native conditions that allow ligands to enter or exit the ligand-binding pocket.     

Population Shift Mechanism for Partial Agonism of AMPA Receptor

α-amino-3-hydroxy-5-methyl-4-isoaxazolepropionic acid (AMPA) ionotropic glutamate receptors mediate fast excitatory neurotransmission in the central nervous system, and their dysfunction is associated with neurological diseases. Glutamate binding to ligand-binding domains (LBDs) of AMPA receptors induces channel opening in the transmembrane domains of the receptors. The T686A mutation reduces glutamate efficacy so that the glutamate behaves as a partial agonist. The crystal structures of wild-type and mutant LBDs are very similar and cannot account for the observed behavior. To elucidate the molecular mechanism inducing partial agonism of the T686A mutant, we computed the free-energy landscapes governing GluA2 LBD closure using replica-exchange umbrella sampling simulations. A semiclosed state, not observed in crystal structures, appears in the mutant during simulation. In this state, the LBD cleft opens slightly because of breaking of interlobe hydrogen bonds, reducing the efficiency of channel opening. The energy difference between the LBD closed and semiclosed states is small, and transitions between the two states would occur by thermal fluctuations. Evidently, glutamate binding to the T686A mutant induces a population shift from a closed to a semiclosed state, explaining the partial agonism in the AMPA receptor.