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.     

核雌激素受体配体结合域的晶体结构特性与功能

雌激素或类雌激素活性物质通过细胞核雌激素受体（nuclear estrogen receptor, nER）通路发挥相应的生理性作用。当这些配体被nER的配体结合域（ligand binding domain, LBD）识别后进入疏水性配体结合空腔内并引起受体构象发生改变,使得原先处于高度活动性的helix 12(H12)被固定从而进一步稳定空腔结构|同时nER也能通过招募一系列辅助调节因子及其他共调节蛋白质,最终调控基因转录。但是,由于不同的配体和受体结合形成的晶体结构并不完全相同,导致这些复合体具有不同的性质,从而影响基因的转录活性。本文综述了nER配体结合域及结合配体后形成的相应晶体结构与活性以及不同配体对受体结构和基因转录的影响。     

Tyrosine and peptide reorientational mobility in polymer solutions: Time-dependent fluorescence anisotropy measurements

The ability of peptides to form biologically active conformations that bind to receptors is governed by their dynamics and their propensity to form stable structures. Such factors are consequently important in the design of peptide drugs. Moreover, the stability of such peptides depends on interactions of the peptide with the surrounding matrix. In this article, we study the effect of the polymer poly(vinyl pyrrolidone) (PVP) on the mobility and orientational dynamics of tyrosine and a model peptide, Val-Tyr-Pro-Asn-Gly-Ala (VYPNGA) in glycerol-water solutions. Orientational dynamics are studied experimentally by time-resolved fluorescence anisotropy decays of tyrosine. The presence of PVP leads to the possibility of a distribution of environments for the peptide. The orientational dynamics of tyrosine show that the probe molecule experiences two very different environments. In one, tyrosine rotational motion is weakly coupled to PVP, while in the other, tyrosine interacts strongly with PVP leading to much slower rotational times. The dynamics of VYPNGA are more complex. Fast intramolecular, localized reorientations of the tyrosine are detected. The temperature dependence of the reorientational dynamics of the tyrosine side chain reveal that these motions are shielded from solvent friction. In contrast, global motions of the peptide are severely restricted by PVP, suggesting the ability of the polymer to restrict peptide mobility.     

Exploring the binding site structure of the PPAR gamma ligand-binding domain by computational solvent mapping

Solvent mapping moves molecular probes, small organic molecules containing various functional groups, around the protein surface, finds favorable positions, clusters the conformations, and ranks the clusters based on the average free energy. Using at least six different solvents as probes, the probes cluster in major pockets of the functional site, providing detailed and reliable information on the amino acid residues that are important for ligand binding. Solvent mapping was applied to 12 structures of the peroxisome proliferator activated receptor gamma (PPARgamma) ligand-binding domain (LBD), including 2 structures without a ligand, 2 structures with a partial agonist, and 8 structures with a PPAR agonist bound. The analysis revealed 10 binding "hot spots", 4 in the ligand-binding pocket, 2 in the coactivator-binding region, 1 in the dimerization domain, 2 around the ligand entrance site, and 1 minor site without a known function. Mapping is a major source of information on the role and cooperativity of these sites. It shows that large portions of the ligand-binding site are already formed in the PPARgamma apostructure, but an important pocket near the AF-2 transactivation domain becomes accessible only in structures that are cocrystallized with strong agonists. Conformational changes were seen in several other sites, including one involved in the stabilization of the LBD and two others at the region of the coactivator binding. The number of probe clusters retained by these sites depends on the properties of the bound agonist, providing information on the origin of correlations between ligand and coactivator binding.     

Protein rotational dynamics investigated with a dual EPR/optical molecular probe. Spin-labeled eosin.

An acyl spin-label derivative of 5-aminoeosin (5-SLE) was chemically synthesized and employed in studies of rotational dynamics of the free probe and of the probe when bound noncovalently to bovine serum albumin using the spectroscopic techniques of fluorescence anisotropy decay and electron paramagnetic resonance (EPR) and their long-lifetime counterparts phosphorescence anisotropy decay and saturation transfer EPR. Previous work (Beth, A. H., Cobb, C. E., and J. M. Beechem, 1992. Synthesis and characterization of a combined fluorescence, phosphorescence, and electron paramagnetic resonance probe. Society of Photo-Optical Instrumentation Engineers. Time-Resolved Laser Spectroscopy III. 504-512) has shown that the spin-label moiety only slightly altered the fluorescence and phosphorescence lifetimes and quantum yields of 5-SLE when compared with 5-SLE whose nitroxide had been reduced with ascorbate and with the diamagnetic homolog 5-acetyleosin. In the present work, we have utilized time-resolved fluorescence anisotropy decay and linear EPR spectroscopies to observe and quantitate the psec motions of 5-SLE in solution and the nsec motions of the 5-SLE-bovine serum albumin complex. Time-resolved phosphorescence anisotropy decay and saturation transfer EPR studies have been carried out to observe and quantitate the microseconds motions of the 5-SLE-albumin complex in glycerol/buffer solutions of varying viscosity. These latter studies have enabled a rigorous comparison of rotational correlation times obtained from these complementary techniques to be made with a single probe. The studies described demonstrate that it is possible to employ a single molecular probe to carry out the full range of fluorescence, phosphorescence, EPR, and saturation transfer EPR studies. It is anticipated that "dual" molecular probes of this general type will significantly enhance capabilities for extracting dynamics and structural information from macromolecules and their functional assemblies.     

C-terminus of a long alpha-neurotoxin is highly mobile when bound to the nicotinic acetylcholine receptor: a time-resolved fluorescence anisotropy approach

To better understand how alpha-neurotoxins interact with the acetylcholine receptor, four fluorescein isothiocyanate derivatives of the siamemsis alpha-cobratoxin were prepared (conjugated to the epsilon-amino group in Lys(23), Lys(35), Lys(49), or Lys(69)) and the time-resolved fluorescence anisotropy of each conjugate was measured free in solution and bound to the Torpedo acetylcholine receptor. All the conjugated reporter groups displayed a high and comparable level of mobility free in solution. When receptor bound, on the other hand, significant differences in the conformational dynamics of the reporter groups were observed with the C-terminal Lys(69) derivative displaying by far the greatest mobility strongly suggesting that the C-terminal domain of the bound neurotoxin is highly mobile and does not participate in the toxin-nAChR binding surface. Additionally, this study demonstrates the utility of time-resolved fluorescence anisotropy to characterize the interaction of heteroproteins.     

Microsecond rotational motions of eosin-labeled myosin measured by time-resolved anisotropy of absorption and phosphorescence

We have studied submicrosecond and microsecond rotational motions within the contractile protein myosin by observing the time-resolved anisotropy of both absorption and emission from the long-lived triplet state of eosin-5-iodoacetamide covalently bound to a specific site on the myosin head. These results, reporting anisotropy data up to 50 microseconds after excitation, extend by two orders of magnitude the time range of data on time-resolved site-specific probe motion in myosin. Optical and enzymatic analyses of the labeled myosin and its chymotryptic digests show that more than 95% of the probe is specifically attached to sulfhydryl-1 (SH1) on the myosin head. In a solution of labeled subfragment-1 (S-1) at 4 degrees C, absorption anisotropy at 0.1 microseconds after a laser pulse is about 0.27. This anisotropy decays exponentially with a rotational correlation time of 210 ns, in good agreement with the theoretical prediction for end-over-end tumbling of S-1, and with times determined previously by fluorescence and electron paramagnetic resonance. In aqueous glycerol solutions, this correlation time is proportional to viscosity/temperature in the microsecond time range. Furthermore, binding to actin greatly restricts probe motion. Thus the bound eosin is a reliable probe of myosin-head rotational motion in the submicrosecond and microsecond time ranges. Our submicrosecond data for myosin monomers (correlation time 400 ns) also agree with previous results using other techniques, but we also detect a previously unresolvable slower decay component (correlation time 2.6 microseconds), indicating that the faster motions are restricted in amplitude. This restriction is not consistent with the commonly accepted free-swivel model of S-1 attachment in myosin. In synthetic thick filaments of myosin, both fast (700 ns) and slow (5 microseconds) components of anisotropy decay are observed. In contrast to the data for monomers, the anisotropy of filaments has a substantial residual component (26% of the initial anisotropy) that does not decay to zero even at times as long as 50 microseconds, implying significant restriction in overall rotational amplitude. This result is consistent with motion restricted to a cone half-angle of about 50 degrees. The combined results are consistent with a model in which myosin has two principal sites of segmental flexibility, one giving rise to submicrosecond motions (possibly corresponding to the junction between S-1 and S-2) and the other giving rise to microsecond motions (possibly corresponding to the junction between S-2 and light meromyosin).(ABSTRACT TRUNCATED AT 400 WORDS)     

Enhanced resolution of fluorescence anisotropy decays by simultaneous analysis of progressively quenched samples. Applications to anisotropic rotations and to protein dynamics.

Enhanced resolution of rapid and complex anisotropy decays was obtained by measurement and analysis of data from progressively quenched samples. Collisional quenching by acrylamide was used to vary the mean decay time of indole or of the tryptophan fluorescence from melittin. Anisotropy decays were obtained from the frequency-response of the polarized emission at frequencies from 4 to 2,000 MHz. Quenching increases the fraction of the total emission, which occurs on the subnanosecond timescale, and thereby provides increased information on picosecond rotational motions or local motions in proteins. For monoexponential subnanosecond anisotropy decays, enhanced resolution is obtained by measurement of the most highly quenched samples. For complex anisotropy decays, such as those due to both local motions and overall protein rotational diffusion, superior resolution is obtained by simultaneous analysis of data from quenched and unquenched samples. We demonstrate that measurement of quenched samples greatly reduces the uncertainty of the 50-ps correlation time of indole in water at 20 degrees C, and allows resolution of the anisotropic rotation of indole with correlation times of 140 and 720 ps. The method was applied to melittin in the monomeric and tetrameric forms. With increased quenching, the anisotropy data showed decreasing contributions from overall protein rotation and increased contribution from picosecond tryptophan motions. The tryptophan residues in both the monomeric and the tetrameric forms of melittin displayed substantial local motions with correlation times near 0.16 and 0.06 ns, respectively. The amplitude of the local motion is twofold less in the tetramer. These highly resolved anisotropy decays should be valuable for comparison with molecular dynamics simulations of melittin.