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

A mutant form of ERα associated with estrogen insensitivity affects the coupling between ligand binding and coactivator recruitment

A homozygous missense mutation in the gene encoding the estrogen receptor α (ERα) was previously identified in a female patient with estrogen insensitivity syndrome. We investigated the molecular features underlying the impaired transcriptional response of this mutant (ERα-Q375H) and four other missense mutations at this position designed to query alternative mechanisms. The identity of residue 375 greatly affected the sensitivity of the receptor to agonists without changing the ligand binding affinity. Instead, the mutations caused changes in the affinity of coactivator binding and alterations in the balance of coactivator and corepressor recruitment. Comparisons among the transcriptional regulatory responses of these six ERα genotypes to a set of ER agonists showed that both steric and electrostatic factors contributed to the functional deficits in gene regulatory activity of the mutant ERα proteins. ERα-coregulator peptide binding in vitro and RIME (rapid immunoprecipitation mass spectrometry of endogenous) analysis in cells showed that the degree of functional impairment paralleled changes in receptor-coregulator binding interactions. These findings uncover coupling between ligand binding and coregulator recruitment that affects the potency rather than the efficacy of the receptor response without substantially altering ligand binding affinity. This highlights a molecular mechanism for estrogen insensitivity syndrome involving mutations that perturb a bidirectional allosteric coupling between ligand binding and coregulator binding that determines receptor transcriptional output.

Ligand Pathways in Nuclear Receptors

Nuclear receptors (NRs) are ligand-inducible transcription factors that play an essential role in a multitude of physiological processes as well as diseases, rendering them attractive drug targets. Crystal structures revealed the binding site of NRs to be buried in the core of the protein, with no obvious route for ligands to access this cavity. The process of ligand binding is known to be an often-neglected contributor to the efficacy of drug candidates and is thought to influence the selectivity and specificity of NRs. While experimental methods generally fail to highlight the dynamic processes of ligand access or egress on the atomistic scale, computational methods have provided fundamental insight into the pathways connecting the buried binding pocket to the surrounding environment. Methods based on molecular dynamics (MD) and Monte Carlo simulations have been applied to identify pathways and quantify their capability to transport ligands. Here, we systematically review findings of more than 20 years of research in the field, including the applied methodology and controversies. Further, we establish a unified nomenclature to describe the pathways with respect to their location relative to protein secondary structure elements and summarize findings relevant to drug design. Lastly, we discuss the effect of NR interaction partners such as coactivators and corepressors, as well as mutations on the pathways.

The nuclear receptor superfamily: A structural perspective

Nuclear receptors (NRs) are a family of transcription factors that regulate numerous physiological processes such as metabolism, reproduction, inflammation, as well as the circadian rhythm. NRs sense changes in lipid metabolite levels to drive differential gene expression, producing distinct physiologic effects. This is an allosteric process whereby binding a cognate ligand and specific DNA sequences drives the recruitment of diverse transcriptional co-regulators at chromatin and ultimately transactivation or transrepression of target genes. Dysregulation of NR signaling leads to various malignances, metabolic disorders, and inflammatory disease. Given their important role in physiology and ability to respond to small lipophilic ligands, NRs have emerged as valuable therapeutic targets. Here, we summarize and discuss the recent progress on understanding the complex mechanism of action of NRs, primarily from a structural perspective. Finally, we suggest future studies to improve our understanding of NR signaling and better design drugs by integrating multiple structural and biophysical approaches.

Structural underpinnings of oestrogen receptor mutations in endocrine therapy resistance

Oestrogen receptor-α (ERα), a key driver of breast cancer, normally requires oestrogen for activation. Mutations that constitutively activate ERα without the need for hormone binding are frequently found in endocrine-therapy-resistant breast cancer metastases and are associated with poor patient outcomes. The location of these mutations in the ER ligand-binding domain and their impact on receptor conformation suggest that they subvert distinct mechanisms that normally maintain the low basal state of wild-type ERα in the absence of hormone. Such mutations provide opportunities to probe fundamental issues underlying ligand-mediated control of ERα activity. Instructive contrasts between these ERα mutations and those that arise in the androgen receptor (AR) during anti-androgen treatment of prostate cancer highlight differences in how activation functions in ERs and AR control receptor activity, how hormonal pressures (deprivation versus antagonism) drive the selection of phenotypically different mutants, how altered protein conformations can reduce antagonist potency and how altered ligand-receptor contacts can invert the response that a receptor has to an agonist ligand versus an antagonist ligand. A deeper understanding of how ligand regulation of receptor conformation is linked to receptor function offers a conceptual framework for developing new anti-oestrogens that might be more effective in preventing and treating breast cancer.

Structural Analysis of the Glucocorticoid Receptor Ligand-Binding Domain in Complex with Triamcinolone Acetonide and a Fragment of the Atypical Coregulator, Small Heterodimer Partner

The synthetic glucocorticoids (GCs) dexamethasone, mometasone furoate, and triamcinolone acetonide are pharmaceutical mainstays to treat chronic inflammatory diseases. These drugs bind to the glucocorticoid receptor (GR), a ligand-activated transcription factor and member of the nuclear receptor superfamily. The GR is widely recognized as a therapeutic target for its ability to counter proinflammatory signaling. Despite the popularity of GCs in the clinic, long-term use leads to numerous side effects, driving the need for new and improved drugs with less off-target pharmacology. X-ray crystal structures have played an important role in the drug-design process, permitting the characterization of robust structure-function relationships. However, steroid receptor ligand-binding domains (LBDs) are inherently unstable, and their crystallization requires extensive mutagenesis to enhance expression and crystallization. Here, we use an ancestral variant of GR as a tool to generate a high-resolution crystal structure of GR in complex with the potent glucocorticoid triamcinolone acetonide (TA) and a fragment of the small heterodimer partner (SHP). Using structural analysis, molecular dynamics, and biochemistry, we show that TA increases intramolecular contacts within the LBD to drive affinity and enhance stability of the receptor-ligand complex. These data support the emerging theme that ligand-induced receptor conformational dynamics at the mouth of the pocket play a major role in steroid receptor activation. This work also represents the first GR structure in complex with SHP, which has been suggested to play a role in modulating hepatic GR function.

Exploring the Structural Compliancy versus Specificity of the Estrogen Receptor Using Isomeric Three-Dimensional Ligands

The estrogen receptors (ERs) bind with high affinity to many structurally diverse ligands by significantly distorting the contours of their ligand-binding pockets. This raises a question: To what degree is ER able to distinguish between structurally related regioisomers and enantiomers? We have explored the structural compliance and specificity of ERα with a set of ligands having a 7-oxa-bicyclo[2.2.1]hept-5-ene sulfonate core and basic side chains typical of selective ER modulators (SERMs). These ligands have two regioisomers, each of which is a racemate of enantiomers. Using orthogonal protecting groups and chiral HPLC, we isolated all 4 isomers and assigned their absolute stereochemistry by X-ray analysis. The 1S,2R,4S isomer has a 80-170-fold higher affinity for ERα than the others, and it profiles as a partial agonist/antagonist in cellular reporter gene assays and in suppressing proliferation of MCF-7 breast cancer cells with subnanomolar potency, far exceeding that of the other isomers. It is the only isomer found bound to ERα by X-ray analysis after crystallization with four-isomer mixtures of closely related analogs. Thus, despite the general compliance of this receptor for binding a large variety of ligand structures, ER demonstrates marked structural specificity and stereospecificity by selecting a single component from a mixture of structurally related isomers to drive ER-regulated cellular activity. Our findings lay the necessary groundwork for seeking unique ER-mediated pharmacological profiles by rational structural perturbations of two different types of side chains in this unprecedented class of ER ligands, which may prove useful in developing more effective endocrine therapies for breast cancer.

The conformational dynamics of H2-H3n and S2-H6 in gating ligand entry into the buried binding cavity of vitamin D receptor

Crystal structures of holo vitamin D receptor (VDR) revealed a canonical conformation in which the ligand is entrapped in a hydrophobic cavity buried in the ligand-binding domain (LBD). The mousetrap model postulates that helix 12 is positioned away from the domain to expose the interior cavity. However, the extended form of helix 12 is likely due to artifacts during crystallization. In this study, we set out to investigate conformational dynamics of apo VDR using molecular dynamics simulation on microsecond timescale. Here we show the neighboring backbones of helix 2-helix 3n and beta strand 2-helix 6 of LBD, instead of the helix 12, undergo large-scale motion, possibly gating the entrance of ligand to the ligand binding domain. Docking analysis to the simulated open structure of VDR with the estimated free energy of -37.0 kJ/mol, would emphasise the role of H2-H3n and S2-H6 in facilitating the entrance of calcitriol to the LBD of VDR.

Unexpected Allosteric Network Contributes to LRH-1 Co-regulator Selectivity

Phospholipids (PLs) are unusual signaling hormones sensed by the nuclear receptor liver receptor homolog-1 (LRH-1), which has evolved a novel allosteric pathway to support appropriate interaction with co-regulators depending on ligand status. LRH-1 plays an important role in controlling lipid and cholesterol homeostasis and is a potential target for the treatment of metabolic and neoplastic diseases. Although the prospect of modulating LRH-1 via small molecules is exciting, the molecular mechanism linking PL structure to transcriptional co-regulator preference is unknown. Previous studies showed that binding to an activating PL ligand, such as dilauroylphosphatidylcholine, favors LRH-1's interaction with transcriptional co-activators to up-regulate gene expression. Both crystallographic and solution-based structural studies showed that dilauroylphosphatidylcholine binding drives unanticipated structural fluctuations outside of the canonical activation surface in an alternate activation function (AF) region, encompassing the β-sheet-H6 region of the protein. However, the mechanism by which dynamics in the alternate AF influences co-regulator selectivity remains elusive. Here, we pair x-ray crystallography with molecular modeling to identify an unexpected allosteric network that traverses the protein ligand binding pocket and links these two elements to dictate selectivity. We show that communication between the alternate AF region and classical AF2 is correlated with the strength of the co-regulator interaction. This work offers the first glimpse into the conformational dynamics that drive this unusual PL-mediated nuclear hormone receptor activation.

Neural-Network Scoring Functions Identify Structurally Novel Estrogen-Receptor Ligands

The magnitude of the investment required to bring a drug to the market hinders medical progress, requiring hundreds of millions of dollars and years of research and development. Any innovation that improves the efficiency of the drug-discovery process has the potential to accelerate the delivery of new treatments to countless patients in need. "Virtual screening," wherein molecules are first tested in silico in order to prioritize compounds for subsequent experimental testing, is one such innovation. Although the traditional scoring functions used in virtual screens have proven useful, improved accuracy requires novel approaches. In the current work, we use the estrogen receptor to demonstrate that neural networks are adept at identifying structurally novel small molecules that bind to a selected drug target, ultimately allowing experimentalists to test fewer compounds in the earliest stages of lead identification while obtaining higher hit rates. We describe 39 novel estrogen-receptor ligands identified in silico with experimentally determined Ki values ranging from 460 nM to 20 μM, presented here for the first time.

Dynamics of nuclear receptor Helix-12 switch of transcription activation by modeling time-resolved fluorescence anisotropy decays

Nuclear hormone receptors (NRs) are major targets for pharmaceutical development. Many experiments demonstrate that their C-terminal Helix (H12) is more flexible in the ligand-binding domains (LBDs) without ligand, this increased mobility being correlated with transcription repression and human diseases. Crystal structures have been obtained in which the H12 is extended, suggesting the possibility of large amplitude H12 motions in solution. However, these structures were interpreted as possible crystallographic artifacts, and thus the microscopic nature of H12 movements is not well known. To bridge the gap between experiments and molecular models and provide a definitive picture of H12 motions in solution, extensive molecular dynamics simulations of the peroxisome proliferator-activated receptor-γ LBD, in which the H12 was bound to a fluorescent probe, were performed. A direct comparison of the modeled anisotropy decays to time-resolved fluorescence anisotropy experiments was obtained. It is shown that the decay rates are dependent on the interactions of the probe with the surface of the protein, and display little correlation with the flexibility of the H12. Nevertheless, for the probe to interact with the surface of the LBD, the H12 must be folded over the body of the LBD. Therefore, the molecular mobility of the H12 should preserve the globularity of the LBD, so that ligand binding and dissociation occur by diffusion through the surface of a compact receptor. These results advance the comprehension of both ligand-bound and ligand-free receptor structures in solution, and also guide the interpretation of time-resolved anisotropy decays from a molecular perspective, particularly by the use of simulations.

Mutation of the palmitoylation site of estrogen receptor α in vivo reveals tissue-specific roles for membrane versus nuclear actions

Estrogen receptor alpha (ERα) activation functions AF-1 and AF-2 classically mediate gene transcription in response to estradiol (E2). A fraction of ERα is targeted to plasma membrane and elicits membrane-initiated steroid signaling (MISS), but the physiological roles of MISS in vivo are poorly understood. We therefore generated a mouse with a point mutation of the palmitoylation site of ERα (C451A-ERα) to obtain membrane-specific loss of function of ERα. The abrogation of membrane localization of ERα in vivo was confirmed in primary hepatocytes, and it resulted in female infertility with abnormal ovaries lacking corpora lutea and increase in luteinizing hormone levels. In contrast, E2 action in the uterus was preserved in C451A-ERα mice and endometrial epithelial proliferation was similar to wild type. However, E2 vascular actions such as rapid dilatation, acceleration of endothelial repair, and endothelial NO synthase phosphorylation were abrogated in C451A-ERα mice. A complementary mutant mouse lacking the transactivation function AF-2 of ERα (ERα-AF2(0)) provided selective loss of function of nuclear ERα actions. In ERα-AF2(0), the acceleration of endothelial repair in response to estrogen-dendrimer conjugate, which is a membrane-selective ER ligand, was unaltered, demonstrating integrity of MISS actions. In genome-wide analysis of uterine gene expression, the vast majority of E2-dependent gene regulation was abrogated in ERα-AF2(0), whereas in C451A-ERα it was nearly fully preserved, indicating that membrane-to-nuclear receptor cross-talk in vivo is modest in the uterus. Thus, this work genetically segregated membrane versus nuclear actions of a steroid hormone receptor and demonstrated their in vivo tissue-specific roles.

Phospholipid--driven gene regulation

Phospholipids (PLs), well known for their fundamental role in cellular structure, play critical signaling roles via their derivatives and cleavage products acting as second messengers in signaling cascades. Recent work has shown that intact PLs act as signaling molecules in their own right by modulating the activity of nuclear hormone transcription factors responsible for tuning genes involved in metabolism, lipid flux, steroid synthesis and inflammation. As such, PLs have been classified as novel hormones. This review highlights recent work in PL-driven gene regulation with a focus on the unique structural features of phospholipid-sensing transcription factors and what sets them apart from well known soluble phospholipid transporters.

Structural and mechanistic insights into bisphenols action provide guidelines for risk assessment and discovery of bisphenol A substitutes

Bisphenol A (BPA) is an industrial compound and a well known endocrine-disrupting chemical with estrogenic activity. The widespread exposure of individuals to BPA is suspected to affect a variety of physiological functions, including reproduction, development, and metabolism. Here we report that the mechanisms by which BPA and two congeners, bisphenol AF and bisphenol C (BPC), bind to and activate estrogen receptors (ER) α and β differ from that used by 17β-estradiol. We show that bisphenols act as partial agonists of ERs by activating the N-terminal activation function 1 regardless of their effect on the C-terminal activation function 2, which ranges from weak agonism (with BPA) to antagonism (with BPC). Crystallographic analysis of the interaction between bisphenols and ERs reveals two discrete binding modes, reflecting the different activities of compounds on ERs. BPA and 17β-estradiol bind to ERs in a similar fashion, whereas, with a phenol ring pointing toward the activation helix H12, the orientation of BPC accounts for the marked antagonist character of this compound. Based on structural data, we developed a protocol for in silico evaluation of the interaction between bisphenols and ERs or other members of the nuclear hormone receptor family, such as estrogen-related receptor γ and androgen receptor, which are two known main targets of bisphenols. Overall, this study provides a wealth of tools and information that could be used for the development of BPA substitutes devoid of nuclear hormone receptor-mediated activity and more generally for environmental risk assessment.

Tunable, mixed-resolution modeling using library-based Monte Carlo and graphics processing units

Building on our recently introduced library-based Monte Carlo (LBMC) approach, we describe a flexible protocol for mixed coarse-grained (CG)/all-atom (AA) simulation of proteins and ligands. In the present implementation of LBMC, protein side chain configurations are pre-calculated and stored in libraries, while bonded interactions along the backbone are treated explicitly. Because the AA side chain coordinates are maintained at minimal run-time cost, arbitrary sites and interaction terms can be turned on to create mixed-resolution models. For example, an AA region of interest such as a binding site can be coupled to a CG model for the rest of the protein. We have additionally developed a hybrid implementation of the generalized Born/surface area (GBSA) implicit solvent model suitable for mixed-resolution models, which in turn was ported to a graphics processing unit (GPU) for faster calculation. The new software was applied to study two systems: (i) the behavior of spin labels on the B1 domain of protein G (GB1) and (ii) docking of randomly initialized estradiol configurations to the ligand binding domain of the estrogen receptor (ERα). The performance of the GPU version of the code was also benchmarked in a number of additional systems.

Antidiabetic phospholipid-nuclear receptor complex reveals the mechanism for phospholipid-driven gene regulation

The human nuclear receptor liver receptor homolog-1 (LRH-1) has an important role in controlling lipid and cholesterol homeostasis and is a potential target for the treatment of diabetes and hepatic diseases. LRH-1 is known to bind phospholipids, but the role of phospholipids in controlling LRH-1 activation remains highly debated. Here we describe the structure of both apo LRH-1 and LRH-1 in complex with the antidiabetic phospholipid dilauroylphosphatidylcholine (DLPC). Together with hydrogen-deuterium exchange MS and functional data, our studies show that DLPC binding is a dynamic process that alters co-regulator selectivity. We show that the lipid-free receptor undergoes previously unrecognized structural fluctuations, allowing it to interact with widely expressed co-repressors. These observations enhance our understanding of LRH-1 regulation and highlight its importance as a new therapeutic target for controlling diabetes.

Deciphering modern glucocorticoid cross-pharmacology using ancestral corticosteroid receptors

Steroid receptors (SRs) are the largest family of metazoan transcription factors and control genes involved in development, endocrine signaling, reproduction, immunity, and cancer. The entire hormone receptor system is driven by a molecular switch triggered by the binding of small lipophilic ligands. This makes the SRs ideal pharmaceutical targets, yet even the best clinically approved synthetic steroidal agonists are prone to cross-reactivity and off-target pharmacology. The mechanism underlying this promiscuity is derived from the fact that SRs share common structural features derived from their evolutionary relationship. More often than not, rational attempts to probe SR drug selectivity via mutagenesis fail even when high quality structural and functional data are available due to the fact that important mutations often result in nonfunctional receptors. This highlights the fact that SRs suffer from instability, preventing in-depth mutational analysis and hampering crystallization of key receptor-ligand complexes. We have taken a unique approach to address this problem by using a resurrected ancestral protein to determine the structure of a previously intractable complex and identified the structural mechanisms that confer activation and selectivity for a widely used glucocorticoid, mometasone furoate. Moreover, we have identified a single residue located outside of the ligand-binding pocket that controls mometasone furoate antagonism versus agonism in the human mineralocorticoid receptor.

Binding of the ERα and ARNT1 AF2 domains to exon 21 of the SRC1 isoform SRC1e is essential for estrogen- and dioxin-related transcription

Steroid receptor co-activator 1 (SRC1) is a transcriptional co-activator of numerous transcription factors involving nuclear receptors. Aryl hydrocarbon receptor nuclear translocator 1 (ARNT1) is an obligatory transcriptional partner of the aryl hydrocarbon receptor (AhR) and hypoxia inducible factor-1α (HIF-1α), as well as a co-activator of estrogen receptors (ERs). To initiate transcription, the activation function 2 (AF2) domains of estrogen-activated ERs interact with LxxLL motifs in the nuclear receptor interaction domain (NID) of SRC1. Here we describe an estrogen and LxxLL domain-independent ERα AF2 binding to SRC1e exon 21. In addition, we found an AF2 domain in exon 16 of ARNT1 that also binds to SRC1e exon 21. Surprisingly, the interaction between SRC1e exon 21 and the AF2 domain of ERα functions as a crucial enhancer of estrogen-induced transcription. The binding of ARNT1 AF2 to SRC1e exon 21 enhances the transcriptional response to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), but the upregulation essentially depends on two cyclin destruction boxes (D-boxes), which are also located on exon 16 of ARNT1. Our findings reveal that a binding site for ERα and ARNT1 AF2 domains in the C-terminus of SRC1e upregulates estrogen- and TCDD-related responses in mammalian cells.

Structure of estradiol metal chelate and estrogen receptor complex: the basis for designing a new class of selective estrogen receptor modulators

Selective estrogen receptor modulators, such as 17β-estradiol derivatives bound to metal complexes, have been synthesized as targeted probes for the diagnosis and treatment of breast cancer. Here, we report the detailed 3D structure of estrogen receptor α ligand-binding domain (ERα-LBD) bound with a novel estradiol-derived metal complex, estradiol-pyridine tetra acetate europium(III), at 2.6 Å resolution. This structure provides important information pertinent to the design of novel functional ERα targeted probes for clinical applications.

Guanidinium chloride-induced spectral perturbations of 4,4'-dianilino-1,1'-binaphthyl-5,5'-disulfonic acid confound interpretation of data on molten globule states

We describe limitations in the use of 4,4'-dianilino-1,1'-binaphthyl-5,5'-disulfonic acid (bis-ANS) to examine unfolding intermediates associated with guanidinium chloride (GuHCl)-induced protein denaturation. Several studies have used alterations in fluorescence emission of bis-ANS to quantify the population of "molten globule" states. Our findings indicate that the observed changes in bis-ANS spectroscopic properties could originate from the interactions of bis-ANS and GuHCl and the aggregation of the dye at higher GuHCl concentrations. We posit that in the absence of additional complementary structural or spectroscopic measurements, the use of bis-ANS emission alone to monitor protein conformations can be misleading.

Synthesis and functional analysis of novel bivalent estrogens

The steroid hormone estrogen plays a critical role in female development and homeostasis. Estrogen mediates its effects through binding and activation of specific estrogen receptors alpha (ERalpha) and beta (ERbeta), members of the steroid/nuclear receptor family of ligand-induced transcription factors. Due to their intimate roles in genomic and nongenomic signaling pathways, these hormones and their receptors have been also implicated in the pathologies of a variety of cancers and metabolic disorders, and have been the target of large therapeutic development efforts. The binding of estrogen to its respective receptors initiates a cascade of events that include receptor dimerization, nuclear localization, DNA binding and recruitment of co-regulatory protein complexes. In this manuscript, we investigate the potential for manipulating steroid receptor gene expression activity through the development of bivalent steroid hormones that are predicted to facilitate hormone receptor dimerization events. Data are presented for the development and testing of novel estrogen dimers, linked through their C-17 moiety, that can activate estrogen receptor alpha (ERalpha)-mediated transcription events with efficacy and potency equal to or greater than that of ERalpha's cognate ligand, 17beta-estradiol. These bivalent estrogen structures open the door to the development of a variety of steroid therapeutics that could dramatically impact future drug development in this area.

Proposal for a role of the Hsp90/Hsp70-based chaperone machinery in making triage decisions when proteins undergo oxidative and toxic damage

The Hsp90/Hsp70-based chaperone machinery plays a well-established role in signaling protein function, trafficking and turnover. A number of recent observations also support the notion that Hsp90 and Hsp70 play key roles in the triage of damaged and aberrant proteins for degradation via the ubiquitin-proteasome pathway. In the mid-1990s, it was discovered that Hsp70 is required for ubiquitin-dependent degradation of short-lived and abnormal proteins, and it became clear that inhibition of Hsp90 uniformly leads to the proteasomal degradation of Hsp90 client proteins. Subsequently, CHIP and parkin were shown to be Hsp70-binding ubiquitin E3 ligases that direct ubiquitin-charged E2 enzymes to the Hsp70-bound client protein. Stabilization by Hsp90 reflects the interaction of the chaperone with the ligand binding cleft of the client protein. These hydrophobic clefts must be open to allow passage of ligands to binding sites in the protein interior, and they are inherent sites of conformational instability. Hsp90 stabilizes the open state of the cleft and prevents Hsp70-dependent ubiquitination. In the model we propose here, progressive oxidative events result in cleft opening as the initial step in protein unfolding, and as long as Hsp90 can interact to stabilize the cleft, it will buffer the effect of oxidative damage. When cleft opening is such that Hsp90 can no longer interact, Hsp70-dependent ubiquitination occurs. We summarize evidence that Hsp90 interacts very dynamically with a variety of proteins that are not classic Hsp90 clients, and we show that this dynamic cycling of Hsp90 with nitric oxide synthase protects against CHIP-mediated ubiquitination. Scientific interest to date has focused on stringent regulation of the classic client proteins, which have metastable clefts and are inherently short lived. But, the recognition that Hsp90 cycles dynamically with longer lived proteins with more stable clefts may permit extension of the triage model to the quality control of damaged proteins in general.

Ligand unbinding from the estrogen receptor: a computational study of pathways and ligand specificity

The estrogen receptor (ER) belongs to the nuclear receptor superfamily, whose members regulate important cellular events like development and metabolism. The ER functions as a transcription regulator and can be activated on ligand binding. Consequently, ligand binding and unbinding constitute fundamental processes in the regulation of genes. Even though both biochemical and structural data of ER are available, the actual mechanism of the ligand binding/unbinding remains elusive. We have performed computational studies on the unbinding mechanism of ERalpha and ERbeta, in the presence of cofactors and with ligands ranging from agonist to a full antagonist. Our results show that agonists or selective ER modulators can dissociate from the receptor through multiple pathways with minor effect on the receptor structure, whereas an antagonist requires larger conformational changes. Furthermore, a specific receptor/ligand combination can exhibit a pathway preference depending on character and conformation of both parts. Hence, it is possible that the binding/unbinding mechanism can explain ligand subtype specificity and thus have an impact in drug discovery.

Kinetic and thermodynamic characterization of dihydrotestosterone-induced conformational perturbations in androgen receptor ligand-binding domain

Ligand-induced conformational perturbations in androgen receptor (AR) are important in coactivator recruitment and transactivation. However, molecular rearrangements in AR ligand-binding domain (AR-LBD) associated with agonist binding and their kinetic and thermodynamic parameters are poorly understood. We used steady-state second-derivative absorption and emission spectroscopy, pressure and temperature perturbations, and 4,4'-bis-anilinonaphthalene 8-sulfonate (bis-ANS) partitioning to determine the kinetics and thermodynamics of the conformational changes in AR-LBD after dihydrotestosterone (DHT) binding. In presence of DHT, the second-derivative absorption spectrum showed a red shift and a change in peak-to-peak distance. Emission intensity increased upon DHT binding, and center of spectral mass was blue shifted, denoting conformational changes resulting in more hydrophobic environment for tyrosines and tryptophans within a more compact DHT-bound receptor. In pressure perturbation calorimetry, DHT-induced energetic stabilization increased the Gibbs free energy of unfolding to 8.4 +/- 1.3 kcal/mol from 3.5 +/- 1.6 kcal/mol. Bis-ANS partitioning studies revealed that upon DHT binding, AR-LBD underwent biphasic rearrangement with a high activation energy (13.4 kcal/mol). An initial, molten globule-like burst phase (k approximately 30 sec(-1)) with greater solvent accessibility was followed by rearrangement (k approximately 0.01 sec(-1)), leading to a more compact conformation than apo-AR-LBD. Molecular simulations demonstrated unique sensitivity of tyrosine and tryptophan residues during pressure unfolding with rearrangement of residues in the coactivator recruitment surfaces distant from the ligand-binding pocket. In conclusion, DHT binding leads to energetic stabilization of AR-LBD domain and substantial rearrangement of residues distant from the ligand-binding pocket. DHT binding to AR-LBD involves biphasic receptor rearrangement including population of a molten globule-like intermediate state.

Ligand dissociation from estrogen receptor is mediated by receptor dimerization: evidence from molecular dynamics simulations

Estrogen Receptor (ER) is an important target for pharmaceutical design. Like other ligand-dependent transcription factors, hormone binding regulates ER transcriptional activity. Nevertheless, the mechanisms by which ligands enter and leave ERs and other nuclear receptors remain poorly understood. Here, we report results of locally enhanced sampling molecular dynamics simulations to identify dissociation pathways of two ER ligands [the natural hormone 17beta-estradiol (E(2)) and the selective ER modulator raloxifene (RAL)] from the human ERalpha ligand-binding domain in monomeric and dimeric forms. E(2) dissociation occurs via three different pathways in ER monomers. One resembles the mousetrap mechanism (Path I), involving repositioning of helix 12 (H12), others involve the separation of H8 and H11 (Path II), and a variant of this pathway at the bottom of the ligand-binding domain (Path II'). RAL leaves the receptor through Path I and a Path I variant in which the ligand leaves the receptor through the loop region between H11 and H12 (Path I'). Remarkably, ER dimerization strongly suppresses Paths II and II' for E(2) dissociation and modifies RAL escape routes. We propose that differences in ligand release pathways detected in the simulations for ER monomers and dimers provide an explanation for previously observed effects of ER quaternary state on ligand dissociation rates and suggest that dimerization may play an important, and hitherto unexpected, role in regulation of ligand dissociation rates throughout the nuclear receptor family.

What is the relationship between the global structures of apo and holo proteins?

It is well known that ligand binding and release may induce a wide range of structural changes in a receptor protein, varying from small movements of loops or side chains in the binding pocket to large-scale domain hinge-bending and shear motions or even partial unfolding that facilitates the capture and release of a ligand. An interesting question is what in general are the conformational changes triggered by ligand binding? The aim of this work is analyze the magnitude of structural changes in a protein resulting from ligand binding to assess if the state of ligand binding needs to be included in template-based protein structure prediction algorithms. To address this issue, a nonredundant dataset of 521 paired protein structures in the ligand-free and ligand-bound form was created and used to estimate the degree of both local and global structure similarity between the apo and holo forms. In most cases, the proteins undergo relatively small conformational rearrangements of their tertiary structure upon ligand binding/release (most root-mean-square-deviations from native, RMSD, are <1 A). However, a clear difference was observed between single- and multiple-domain proteins. For the latter, RMSD changes greater than 1 A and sometimes larger were found for almost 1/3 of the cases; these are mainly associated with large-scale hinge-bending movements of entire domains. The changes in the mutual orientation of individual domains in multiple-domain proteins upon ligand binding were investigated using a mechanistic model based on mass-weighted principal axes as well as interface buried surface calculations. Some preferences toward the anticipated mechanism of protein domain movements are predictable based on the examination of just the ligand-free structural form. These results have applications to protein structure prediction, particularly in the context of protein domain assembly, if additional information concerning ligand binding is exploited.

Unbinding of retinoic acid from the retinoic acid receptor by random expulsion molecular dynamics

Unbinding pathways of retinoic acid (RA) bound to retinoic acid receptor (RAR) have been explored by the random expulsion molecular dynamics (REMD) method. Our results show that RA may exit the binding site of RAR through flexible regions close to the H1-H3 loop and beta-sheets, without displacing H12 from its agonist position. This result may explain kinetic differences between agonist and antagonist ligands observed for other nuclear receptors. The extended and flexible structure of RA initiated a methodological study in a simplified two-dimensional model system. The REMD force should in general be distributed to all atoms of the ligand to obtain the most unbiased results, but for a ligand which is tightly bound in the binding pocket through a strong electrostatic interaction, application of the REMD force on a single atom is preferred.

Defining the requirements for Hsp40 and Hsp70 in the Hsp90 chaperone pathway

The Hsp90 chaperoning pathway and its model client substrate, the progesterone receptor (PR), have been used extensively to study chaperone complex formation and maturation of a client substrate in a near native state. This chaperoning pathway can be reconstituted in vitro with the addition of five proteins plus ATP: Hsp40, Hsp70, Hop, Hsp90, and p23. The addition of these proteins is necessary to reconstitute hormone-binding capacity to the immuno-isolated PR. It was recently shown that the first step for the recognition of PR by this system is binding by Hsp40. We compared type I and type II Hsp40 proteins and created point mutations in Hsp40 and Hsp70 to understand the requirements for this first step. The type I proteins, Ydj1 and DjA1 (HDJ2), and a type II, DjB1 (HDJ1), act similarly in promoting hormone binding and Hsp70 association to PR, while having different binding characteristics to PR. Ydj1 and DjA1 bind tightly to PR whereas the binding of DjB1 apparently has rapid on and off rates and its binding cannot be observed by antibody pull-down methods using either purified proteins or cell lysates. Mutation studies indicate that client binding, interactions between Hsp40 and Hsp70, plus ATP hydrolysis by Hsp70 are all required to promote conformational maturation of PR via the Hsp90 pathway.

A biosensor for estrogenic substances using the quartz crystal microbalance

This article describes a biosensor that detects estrogenic substances using a quartz crystal microbalance with a genetically engineered construct of the hormone-binding domain of the alpha-estrogen receptor. The receptor was immobilized to a piezoelectric quartz crystal via a single exposed cysteine, forming a uniform orientation on the crystal surface. Our results illustrate that this sensor responds to a variety of ligands that are known to bind to the estrogen receptor. No response was observed for nonbinding substances such as testosterone and progesterone. The sensitive response of this biosensor to estrogenic substances results from changes in the structural rigidity of the immobilized receptor that occurs with ligand binding. Agonist and antagonist show different responses.

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.

Kinase-mediated transcription, activators of nongenotropic estrogen-like signaling (ANGELS), and osteoporosis: a different perspective on the HRT dilemma

Studies in bone, as well as other nonreproductive target tissues of sex steroid, like the cardiovascular and the central nervous system (CNS), have elucidated a previously unappreciated mechanism of sex steroid action involving the rapid activation of mitogen-activated protein kinases and/or phosphatidyl inositol 3 kinase, and consequent potent regulatory affects on the transcription of a set of genes that is distinct from that regulated through classic (genotropic) control of transcription. These actions stem from an unexpected function of the classic nuclear receptors outside the nucleus, most probably from receptor interactions within distinct signal transduction pathways in preassembled scaffolds. Importantly, these nongenotropic actions are mediated by the ligand-binding domain of the receptor and can be functionally dissociated from classic transcriptional activation with synthetic ligands, termed activators of nongenotropic estrogen-like signaling (ANGELS). We highlight this evidence and discuss its pharmacotherapeutic implications vis a vis the dilemmas posed by the recently appreciated shortfalls of postmenopausal hormone replacement therapy.

Coactivator proteins as determinants of estrogen receptor structure and function: spectroscopic evidence for a novel coactivator-stabilized receptor conformation

The direct regulation of gene transcription by nuclear receptors, such as the estrogen receptor (ER), involves not just ligand and DNA binding but the recruitment of coregulators. Typically, recruitment of p160 coactivator proteins to agonist-liganded ER is considered to be unidirectional, with ligand binding stabilizing an ER ligand binding domain (LBD) conformation that favors coactivator interaction. Using fluorophore-labeled ERalpha-LBDs, we present evidence for a pronounced stabilization of ER conformation that results from coactivator binding, manifest by decreased ER sensitivity to proteases and reduced conformational dynamics, as well as for the formation of a novel coactivator-stabilized (costabilized) receptor conformation, that can be conveniently monitored by the generation of an excimer emission from pyrene-labeled ERalpha-LBDs. This costabilized conformation may embody features required to support ER transcriptional activity. Different classes of coactivator proteins combine with estrogen agonists of different structure to elicit varying degrees of this receptor stabilization, and antagonists and coactivator binding inhibitors disfavor the costabilized conformation. Remarkably, high concentrations of coactivators engender this conformation even in apo- and antagonist-bound ERs (more so with selective ER modulators than with pure antagonists), providing an in vitro model for the development of resistance to hormone therapy in breast cancer.

The role of hsp90 in heme-dependent activation of apo-neuronal nitric-oxide synthase

Like other nitric-oxide synthase (NOS) enzymes, neuronal NOS (nNOS) turnover and activity are regulated by the ubiquitous protein chaperone hsp90. We have shown previously that nNOS expressed in Sf9 cells where endogenous heme levels are low is activated from the apo- to the holo-enzyme by addition of exogenous heme to the culture medium, and this activation is inhibited by radicicol, a specific inhibitor of hsp90 (Billecke, S. S., Bender, A. T., Kanelakis, K. C., Murphy, P. J. M., Lowe, E. R., Kamada, Y., Pratt, W. B., and Osawa, Y. (2002) J. Biol. Chem. 278, 15465-15468). In this work, we examine heme binding by apo-nNOS to form the active enzyme in a cell-free system. We show that cytosol from Sf9 cells facilitates heme-dependent apo-nNOS activation by promoting functional heme insertion into the enzyme. Sf9 cytosol also converts the glucocorticoid receptor (GR) to a state where the hydrophobic ligand binding cleft is open to access by steroid. Both cell-free heme activation of purified nNOS and activation of steroid binding activity of the immunopurified GR are inhibited by radicicol treatment of Sf9 cells prior to cytosol preparation, and addition of purified hsp90 to cytosol partially overcomes this inhibition. Although there is an hsp90-dependent machinery in Sf9 cytosol that facilitates heme binding by apo-nNOS, it is clearly different from the machinery that facilitates steroid binding by the GR. hsp90 regulation of apo-nNOS heme activation is very dynamic and requires higher concentrations of radicicol for its inhibition, whereas GR steroid binding is determined by assembly of stable GR.hsp90 heterocomplexes that are formed by a purified five-chaperone machinery that does not activate apo-nNOS.

The Hsp90 Cochaperone p23 Is the Limiting Component of the Multiprotein Hsp90/Hsp70-based Chaperone System in Vivo Where It Acts to Stabilize the Client Protein·Hsp90 Complex

A variety of signaling proteins form heterocomplexes with and are regulated by the heat shock protein chaperone hsp90. These complexes are formed by a multiprotein machinery, including hsp90 and hsp70 as essential and abundant components and Hop, hsp40, and p23 as non-essential cochaperones that are present in much lower abundance in cells. Overexpression of signaling proteins can overwhelm the capacity of this machinery to properly assemble heterocomplexes with hsp90. Here, we show that the limiting component of this assembly machinery in vitro in reticulocyte lysate and in vivo in Sf9 cells is p23. Only a fraction of glucocorticoid receptors (GR) overexpressed in Sf9 cells are in heterocomplex with hsp90 and have steroid binding activity, with the majority of the receptors present as both insoluble and cytosolic GR aggregates. Coexpression of p23 with the GR increases the proportion of cytosolic receptors that are in stable GR.hsp90 heterocomplexes with steroid binding activity, a strictly hsp90-dependent activity for the GR. Coexpression of p23 eliminates the insoluble GR aggregates and shifts the cytosolic receptor from very large aggregates without steroid binding activity to approximately 600-kDa heterocomplexes with steroid binding activity. These data lead us to conclude that p23 acts in vivo to stabilize hsp90 binding to client protein.

Visualization and mechanism of assembly of a glucocorticoid receptor.Hsp70 complex that is primed for subsequent Hsp90-dependent opening of the steroid binding cleft

A minimal system of five proteins, hsp90, hsp70, Hop, hsp40, and p23, assembles glucocorticoid receptor (GR).hsp90 heterocomplexes and causes the simultaneous opening of the steroid binding cleft to access by steroid. The first step in assembly is the ATP-dependent and hsp40 (YDJ-1)-dependent formation of a GR.hsp70 complex that primes the receptor for subsequent ATP-dependent activation by hsp90, Hop, and p23. This study focuses on three aspects of the GR priming reaction with hsp70. First, we have visualized the primed GR.hsp70 complexes by atomic force microscopy, and we find the most common stoichiometry to be 1:1, with some complexes of a size approximately 1:2 and a few complexes of larger size. Second, in a recent study of progesterone receptor priming, it was shown that hsp40 binds first, leading to the notion that it targets hsp70 to the receptor. We show here that hsp40 does not perform such a targeting function in priming the GR. Third, we focus on a short amino-terminal segment of the ligand binding domain that is required for GR.hsp90 heterocomplex assembly. By using two glutathione S-transferase (GST)/ligand binding domain fusions with (GST/520C) and without (GST/554C) hsp90 binding and steroid binding activity, we show that the priming step with hsp70 occurs with GST/554C, and it is the subsequent assembly step with hsp90 that is defective.

Regulation of signaling protein function and trafficking by the hsp90/hsp70-based chaperone machinery

Nearly 100 proteins are known to be regulated by hsp90. Most of these substrates or "client proteins" are involved in signal transduction, and they are brought into complex with hsp90 by a multiprotein hsp90/hsp70-based chaperone machinery. In addition to binding substrate proteins at the chaperone site(s), hsp90 binds cofactors at other sites that are part of the heterocomplex assembly machinery as well as immunophilins that connect assembled substrate*hsp90 complexes to protein-trafficking systems. In the 5 years since we last reviewed this subject, much has been learned about hsp90 structure, nucleotide-binding, and cochaperone interactions; the most important concept is that ATP hydrolysis by an intrinsic ATPase activity results in a conformational change in hsp90 that is required to induce conformational change in a substrate protein. The conformational change induced in steroid receptors is an opening of the steroid-binding cleft so that it can be accessed by steroid. We have now developed a minimal system of five purified proteins-hsp90, hsp70, Hop, hsp40, and p23- that assembles stable receptor*hsp90 heterocomplexes. An hsp90*Hop*hsp70*hsp40 complex opens the cleft in an ATP-dependent process to produce a receptor*hsp90 heterocomplex with hsp90 in its ATP-bound conformation, and p23 then interacts with the hsp90 to stabilize the complex. Stepwise assembly experiments have shown that hsp70 and hsp40 first interact with the receptor in an ATP-dependent reaction to produce a receptor*hsp70*hsp40 complex that is "primed" to be activated to the steroid-binding state in a second ATP-dependent step with hsp90, Hop, and p23. Successful use of the five-protein system with other substrates indicates that it can assemble signal protein*hsp90 heterocomplexes whether the substrate is a receptor, a protein kinase, or a transcription factor. This purified system should facilitate understanding of how eukaryotic hsp70 and hsp90 work together as essential components of a process that alters the conformations of substrate proteins to states that respond in signal transduction.

Mutations at positions 547-553 of rat glucocorticoid receptors reveal that hsp90 binding requires the presence, but not defined composition, of a seven-amino acid sequence at the amino terminus of the ligand binding domain

Glucocorticoid receptors (GRs) must heterocomplex with hsp90 to have an open steroid binding cleft that can be accessed by steroid. We reported that a seven-amino acid sequence (547-553 of rat GR) overlapping the amino-terminal end of the ligand binding domain is required for hsp90 binding to GR. We have now conducted saturation mutagenesis of this sequence, which appears to be part of the surface where the ligand binding cleft merges with the surface of the ligand binding domain. No single point mutation causes significant changes in any of a variety of biochemical and biological properties in addition to hsp90 binding. A triple mutation (P548A/T549A/V551A) increases by >100-fold the steroid concentration required for half-maximal induction without affecting the level of maximal induction or coactivator response. Interestingly, this triple mutant displays reduced binding of steroid and hsp90 in whole cells, but it possesses wild type affinity for steroid and normal hsp90 binding capacity under cell-free conditions. This phenotype of a dramatic shift in the dose response for transactivation would be expected from an increase in the rate of disassembly of the triple mutant GR.hsp90 heterocomplex in the cell. Mutation of the entire seven-amino acid region to CAAAAAC maintains the presence of a critical alpha-helical structure and heterocomplex formation with hsp90 but eliminates steroid binding and transcriptional activation, thus disconnecting hsp90 binding from opening of the ligand binding cleft and steroid binding.