Structure-Based Approach for the Study of Estrogen Receptor Binding Affinity and Subtype Selectivity

Estrogens exert important physiological effects through the modulation of two human estrogen receptor (hER) subtypes, alpa (hER alpha) and beta (hER beta). Because the levels and relative proportion of hER alpha and hER beta differ significantly in different target cells, selective hER ligands could target specific tissues or pathways regulated by one receptor subtype without affecting the other. To understand the structural and chemical basis by which small molecule modulators are able to discriminate between the two subtypes, we have applied three-dimensional target-based approaches employing a series of potent hER-ligands. Comparative molecular field analysis (CoMFA) studies were applied to a data set of 81 hER modulators, for which binding affinity values were collected for both hER alpha and hER beta. Significant statistical coefficients were obtained (hER alpha, q(2) = 0.76; hER beta, q(2) = 0.70), indicating the internal consistency of the models. The generated models were validated using external test sets, and the predicted values were in good agreement with the experimental results. Five hER crystal structures were used in GRID/PCA investigations to generate molecular interaction fields (MIF) maps. hER alpha and hER beta were separated using one factor. The resulting 3D information was integrated with the aim of revealing the most relevant structural features involved in hER subtype selectivity. The final QSAR and GRID/PCA models and the information gathered from 3D contour maps should be useful for the design or novel hER modulators with improved selectivity.

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 17 beta-estradiol (E-2) and the selective ER modulator raloxifene (RAL)] from the human ER alpha 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.

Estrogen receptor: Structural differences and potential implications on selectivity examined by the GRID/CPCA approach

Selective Estrogen Receptor Modulators ( SERMs) have been developed, but the selectivity towards the subtypes ( ER or ER is not well understood. Based on three-dimensional structural properties of ligand binding domains, a model that takes into account this aspect was developed via molecular interaction fields and consensus principal component analysis (GRID/CPCA).

Role of alpha 7 Nicotinic Acetylcholine Receptor in Calcium Signaling Induced by Prion Protein Interaction with Stress-inducible Protein 1

The prion protein (PrP(C)) is a conserved glycosylphosphatidyl-inositol-anchored cell surface protein expressed by neurons and other cells. Stress-inducible protein 1 (STI1) binds PrP(C) extracellularly, and this activated signaling complex promotes neuronal differentiation and neuroprotection via the extracellular signal-regulated kinase 1 and 2 (ERK1/2) and cAMP-dependent protein kinase 1 (PKA) pathways. However, the mechanism by which the PrPC-STI1 interaction transduces extracellular signals to the intracellular environment is unknown. We found that in hippocampal neurons, STI1-PrP(C) engagement induces an increase in intracellular Ca(2+) levels. This effect was not detected in PrP(C)-null neurons or wild-type neurons treated with an STI1 mutant unable to bind PrP(C). Using a best candidate approach to test for potential channels involved in Ca(2+) influx evoked by STI1-PrP(C), we found that alpha-bungarotoxin, a specific inhibitor for alpha 7 nicotinic acetylcholine receptor (alpha 7nAChR), was able to block PrP(C)-STI1-mediated signaling, neuroprotection, and neuritogenesis. Importantly, when alpha 7nAChR was transfected into HEK 293 cells, it formed a functional complex with PrP(C) and allowed reconstitution of signaling by PrP(C)-STI1 interaction. These results indicate that STI1 can interact with the PrP(C).alpha 7nAChR complex to promote signaling and provide a novel potential target for modulation of the effects of prion protein in neurodegenerative diseases.

Analysis of Agonist and Antagonist Effects on Thyroid Hormone Receptor Conformation by Hydrogen/Deuterium Exchange

Thyroid hormone receptors (TRs) are ligand-gated transcription factors with critical roles in development and metabolism. Although x-ray structures of TR ligand-binding domains (LBDs) with agonists are available, comparable structures without ligand (apo-TR) or with antagonists are not. It remains important to understand apo-LBD conformation and the way that it rearranges with ligands to develop better TR pharmaceuticals. In this study, we conducted hydrogen/deuterium exchange on TR LBDs with or without agonist (T(3)) or antagonist (NH(3)). Both ligands reduce deuterium incorporation into LBD amide hydrogens, implying tighter overall folding of the domain. As predicted, mass spectroscopic analysis of individual proteolytic peptides after hydrogen/deuterium exchange reveals that ligand increases the degree of solvent protection of regions close to the buried ligand-binding pocket. However, there is also extensive ligand protection of other regions, including the dimer surface at H10-H11, providing evidence for allosteric communication between the ligand-binding pocket and distant interaction surfaces. Surprisingly, C-terminal activation helix H12, which is known to alter position with ligand, remains relatively protected from solvent in all conditions suggesting that it is packed against the LBD irrespective of the presence or type of ligand. T(3), but not NH(3), increases accessibility of the upper part of H3-H5 to solvent, and we propose that TR H12 interacts with this region in apo-TR and that this interaction is blocked by T(3) but not NH(3.) We present data from site-directed mutagenesis experiments and molecular dynamics simulations that lend support to this structural model of apo-TR and its ligand-dependent conformational changes. (Molecular Endocrinology 25: 15-31, 2011)

Structural modeling of high-affinity thyroid receptor-ligand complexes

Understanding the molecular basis of the binding modes of natural and synthetic ligands to nuclear receptors is fundamental to our comprehension of the activation mechanism of this important class of hormone regulated transcription factors and to the development of new ligands. Thyroid hormone receptors (TR) are particularly important targets for pharmaceuticals development because TRs are associated with the regulation of metabolic rates, body weight, and circulating levels of cholesterol and triglycerides in humans. While several high-affinity ligands are known, structural information is only partially available. In this work we obtain structural models of several TR-ligand complexes with unknown structure by docking high affinity ligands to the receptors` ligand binding domain with subsequent relaxation by molecular dynamics simulations. The binding modes of these ligands are discussed providing novel insights into the development of TR ligands. The experimental binding free energies are reasonably well-reproduced from the proposed models using a simple linear interaction energy free-energy calculation scheme.

Recognition by the Thyroid Hormone Receptor of Canonical DNA Response Elements

To shed more light on the molecular requirements for recognition of thyroid response elements (TRES) by thyroid receptors (TRs), we compared the specific aspects of DNA TRE recognition by different TR constructs. Using fluorescence anisotropy, we performed a detailed and hierarchical study of TR-TRE binding. This wits done by comparing the binding affinities of three different TR constructs for four different TRE DNA elements, including palindromic sequences and direct repeats (F2, PAL, DR-1, and DR-4) as well as their interactions with nonspecific DNA sequences. The effect of MgCl(2) on suppressing of nonselective DNA binding to TR was also investigated. Furthermore, we determined the dissociation constants of the hTR beta DBD (DNA binding domain) and hTR beta DBD-LBD (DNA binding and ligand binding domains) for specific TRES. We found that a minimum DNA recognition peptide derived from DBD (H1TR) is sufficient for recognition and interaction with TREs, whereas scrambled DNA sequences were unrecognized. Additionally, we determined that the TR DBD binds to F2, PAL, and DR-4 with high affinity and similar K(d) values. The TR DBD-LBD recognizes all the tested TRES but binds preferentially to F2, with even higher affinity. Finally, our results demonstrate the important role played by LBDs in modulating TR-DNA binding.

Chemical identity and connections of medial preoptic area neurons expressing melanin-concentrating hormone during lactation

Lactation is an energy-demanding process characterized by massive food and water consumption, cessation of the reproductive cycle and induction of maternal behavior. During lactation, melanin-concentrating hormone (MCH) mRNA and peptide expression are increased in the medial preoptic area (MPO) and in the anterior paraventricular nucleus of the hypothalamus. Here we show that MCH neurons in the MPO coexpress the GABA synthesizing enzyme GAD-67 mRNA. We also show that MCH neurons in the MPO of female rats are innervated by neuropeptides that control energy homeostasis including agouti-related protein (AgRP), alpha-melanocyte stimulating hormone (alpha MSH) and cocaine- and amphetamine-regulated transcript (CART). Most of these inputs originate from the arcuate nucleus neurons. Additionally, using injections of retrograde tracers we found that CART neurons in the ventral premammillary nucleus also innervate the MPO. We then assessed the projections of the female MPO using injections of anterograde tracers. The MPO densely innervates hypothalamic nuclei related to reproductive control including the anteroventral periventricular nucleus, the ventrolateral subdivision of the ventromedial nucleus (VMHvl) and the ventral premammillary nucleus (PMV). We found that the density of MCH-ir fibers is increased in the VMHvl and PMV during lactation. Our findings suggest that the expression of MCH in the MPO may be induced by changing levels of neuropeptides involved in metabolic control. These MCH/GABA neurons may, in turn, participate in the suppression of cyclic reproductive function and/or sexual behavior during lactation through projections to reproductive control sites. (C) 2009 Elsevier B.V. All rights reserved.

Synergistic acceleration of thyroid hormone degradation by phenobarbital and the PPAR alpha agonist WY14643 in rat hepatocytes

Energy balance is maintained by controlling both energy intake and energy expenditure. Thyroid hormones play a crucial role in regulating energy expenditure. Their levels are adjusted by a tight feed back-control led regulation of thyroid hormone production/incretion and by their hepatic metabolism. Thyroid hormone degradation has previously been shown to be enhanced by treatment with phenobarbital or other antiepileptic drugs due to a CAR-dependent induction of phase 11 enzymes of xenobiotic metabolism. We have recently shown, that PPAR alpha agonists synergize with phenobarbital to induce another prototypical CAR target gene, CYP2B1. Therefore, it was tested whether a PPAR alpha agonist could enhance the phenobarbital-dependent acceleration of thyroid hormone elimination. In primary cultures of rat hepatocytes the apparent half-life of T3 was reduced after induction with a combination of phenobarbital and the PPARa agonist WY14643 to a larger extent than after induction with either Compound alone. The synergistic reduction of the half-life could be attributed to a synergistic induction of CAR and the CAR target genes that code for enzymes and transporters involved in the hepatic elimination of T3, such as OATP1A1, OATP1A3, UGT1A3 and UCT1A10. The PPAR alpha-dependent CAR induction and the subsequent induction of T3-eliminating enzymes might be of physiological significance for the fasting-incluced reduction in energy expenditure by fatty acids as natural PPARa ligands. The synergism of the PPAR alpha agonist WY14643 and phenobarbital in inducing thyroid hormone breakdown might serve as a paradigm for the synergistic disruption of endocrine control by other combinations of xenobiotics. (C) 2009 Elsevier Inc. All rights reserved.

Only subtle protein conformational adaptations are required for ligand binding to thyroid hormone receptors: Simulations using a novel multipoint steered molecular dynamics approach

Thyroid hormone receptors (TR) are hormone-dependent transcription regulators that play a major role in human health, development, and metabolic functions. The thyroid hormone resistance syndrome, diabetes, obesity, and some types of cancer are just a few examples of important diseases that are related to TR malfunctioning, particularly impaired hormone binding. Ligand binding to and dissociation from the receptor ultimately control gene transcription and, thus, detailed knowledge of binding and release mechanisms are fundamental for the comprehension of the receptor`s biological function and development of pharmaceuticals. In this work, we present the first computational study of ligand entry into the ligand binding domain (LBD) of a nuclear receptor. We report molecular dynamics simulations of ligand binding to TRs using a generalization of the steered molecular dynamics technique designed to perform single-molecule pulling simulations along arbitrarily nonlinear driving pathways. We show that only gentle protein movements and conformational adaptations are required for ligand entry into the LBDs and that the magnitude of the forces applied to assist ligand binding are of the order of the forces involved in ligand dissociation. Our simulations suggest an alternative view for the mechanisms ligand binding and dissociation of ligands from nuclear receptors in which ligands can simply diffuse through the protein surface to reach proper positioning within the binding pocket. The proposed picture indicates that the large-amplitude protein motions suggested by the apo- and holo-RXR alpha crystallographic structures are not required, reconciling conformational changes of LBDs required for ligand entry with other nuclear receptors apo-structures that resemble the ligand-bound LBDs.

On the Denaturation Mechanisms of the Ligand Binding Domain of Thyroid Hormone Receptors

The ligand binding domain (LBD) of nuclear hormone receptors adopts a very compact, mostly alpha-helical structure that binds specific ligands with very high affinity. We use circular dichroism spectroscopy and high-temperature molecular dynamics Simulations to investigate unfolding of the LBDs of thyroid hormone receptors (TRs). A molecular description of the denaturation mechanisms is obtained by molecular dynamics Simulations of the TR alpha and TR beta LBDs in the absence and in the presence of the natural ligand Triac. The Simulations Show that the thermal unfolding of the LBD starts with the loss of native contacts and secondary Structure elements, while the Structure remains essentially compact, resembling a molten globule state. This differs From most protein denaturation simulations reported to date and suggests that the folding mechanism may start with the hydrophobic collapse of the TR LBDs. Our results reveal that the stabilities of the LBDs of the TR alpha and TR beta Subtypes are affected to different degrees by the binding of the isoform selective ligand Triac and that ligand binding confers protection against thermal denaturation and unfolding in a subtype specific manner. Our Simulations indicate two mechanisms by which the ligand stabilizes the LBD: (1) by enhancing the interactions between H8 and H 11, and the interaction of the region between H I and the Omega-loop with the core of the LBD, and (2) by shielding the hydrophobic H6 from hydration.

Probing conformational changes in orphan nuclear receptor: The NGFI-B intermediate is a partially unfolded dimer

Human nerve growth factor-induced B (NGFI-B) is a member of the NR4A subfamily of orphan nuclear receptors (NRs). Lacking identified ligands, orphan NRs show particular co-regulator proteins binding properties, different from other NRs, and they might have a non-classical quaternary organization. A body of evidence suggests that NRs recognition of and binding to ligands, DNA, homo- and heterodimerization partners and co-regulator proteins involve significant conformational changes of the NR ligand-binding domains (LBDs). To shed light on largely unknown biophysical properties of NGFI-B, here we studied structural organization and unfolding properties of NGFI-B ligand (like)-binding domain induced by chemical perturbation. Our results show that NGFI-B LBD undergoes a two-state guanidine hydrochloride (GndHCl) induced denaturation, as judged by changes in the a-helical content of the protein monitored by circular dichroism spectroscopy (CD). In contrast, changes in the tertiary structure of NGFI-B LBD, reported by intrinsic fluorescence, reveal a clear intermediate state. Additionally, SAXS results demonstrate that the intermediate observed by intrinsic fluorescence is a partially folded homodimeric structure, which further unfolds without dissociation at higher GndHCl concentrations. This partially unfolded dimeric assembly of NGFI-B LBD might resemble an intermediate that this domain access momentarily in the native state upon interactions with functional partners. (C) 2008 Elsevier B.V. All rights reserved.

Adenosine Modulatesa alpha(2)-Adrenergic Receptors within Specific Subnuclei of the Nucleus Tractus Solitarius in Normotensive and Spontaneously Hypertensive Rats

Adenosine Is known to modulate neuronal activity within the nucleus tractus solitarius (NTS). The modulatory effect of adenosine A, receptors (A(1R)) on alpha(2)-adrenoceptors (Adr(2R)) was evaluated using quantitative radioautography within NTS subnuclei and using neuronal culture of normotensive (WKY) and spontaneously hypertensive rats (SHR). Radioautography was used in a saturation experiment to measure Adr2R binding parameters (B(max), K(d)) In the presence of 3 different concentrations of N(6)-cyclopentyladenosine (CPA), an A(1R) agonist. Neuronal culture confirmed our radioautographic results. [(3)H]RX821002, an Adr(2R) antagonist, was used as a ligand for both approaches. The dorsomedial/dorsolateral subnucleus of WKY showed an increase in B(max) values (21%) Induced by 10 nmol/L of CPA. However, the subpostremal subnucleus showed a decrease in Kd values (24%) induced by 10 nmol/L of CPA. SHR showed the same pattern of changes as WKY within the same subnuclei; however, the modulatory effect of CPA was induced by I nmol/L (increased B(max), 17%; decreased K(d), 26%). Cell culture confirmed these results, because 10(-5) and 10(-7) mol/L of CPA promoted an Increase in [3H]RX821002 binding of WKY (53%) and SHR cells (48%), respectively. DPCPX, an AIR antagonist, was used to block the modulatory effect promoted by CPA with respect to Adr2R binding. In conclusion, our study shows for the first time an interaction between A(1R) that increases the binding of Adr2R within specific subnuclei of the NTS. This may be important In understanding the complex autonomic response induced by adenosine within the NTS. In addition, changes in interactions between receptors might be relevant to understanding the development of hypertension. (Hypertens Res 2008; 31: 2177-2186)

Adenosine modulates alpha2-adrenergic receptors through a phospholipase C pathway in brainstem cell culture of rats

Adenosine acts in the nucleus tractus solitarii (NTS), one of the main brain sites related to cardiovascular control. In the present study we show that A(1) adenosine receptor (A(1R)) activation promotes an increase on alpha(2)-adrenoceptor (Alpha(2R)) binding in brainstem cell culture from newborn rats. We investigated the intracellular cascade involved in such modulatory process using different intracellular signaling molecule inhibitors as well as calcium chelators. Phospholipase C, protein kinase Ca(2+)-dependent, IP(3) receptor and intracellular calcium were shown to participate in A(1R)/Alpha(2R) interaction. In conclusion, this result might be important to understand the role of adenosine within the NTS regarding autonomic cardiovascular control. (C) 2009 Elsevier B.V. All rights reserved.

The geodiamolide H, derived from Brazilian sponge Geodia corticostylifera, regulates actin cytoskeleton, migration and invasion of breast cancer cells cultured in three-dimensional environment

We are investigating effects of the depsipeptide geodiamolide H, isolated from the Brazilian sponge Geodia corticostylifera, on cancer cell lines grown in 3D environment. As shown previously geodiamolide H disrupts actin cytoskeleton in both sea urchin eggs and breast cancer cell monolayers. We used a normal mammary epithelial cell line MCF 10A that in 3D assay results formation of polarized spheroids. We also used cell lines derived from breast tumors with different degrees of differentiation: MCF7 positive for estrogen receptor and the Hs578T, negative for hormone receptors. Cells were placed on top of Matrigel. Spheroids obtained from these cultures were treated with geodiamolide H. Control and treated samples were analyzed by light and confocal microscopy. Geodiamolide H dramatically affected the poorly differentiated and aggressive Hs578T cell line. The peptide reverted HsS78T malignant phenotype to polarized spheroid-like structures. MCF7 cells treated by geodiamolide H exhibited polarization compared to controls. Geodiamolide H induced striking phenotypic modifications in Hs578T cell line and disruption of actin cytoskeleton. We investigated effects of geodiamolide H on migration and invasion of Hs578T cells. Time-lapse microscopy showed that the peptide inhibited migration of these cells in a dose-dependent manner. Furthermore invasion assays revealed that geodiamolide H induced a 30% decrease on invasive behavior of Hs578T cells. Our results suggest that geodiamolide H inhibits migration and invasion of Hs578T cells probably through modifications in actin cytoskeleton. The fact that normal cell lines were not affected by treatment with geodiamolide H stimulates new studies towards therapeutic use for this peptide.