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.

Differential effects of TR ligands on hormone dissociation rates: Evidence for multiple ligand entry/exit pathways

Some nuclear receptor (NR) ligands promote dissociation of radiolabeled bound hormone from the buried ligand binding cavity (LBC) more rapidly than excess unlabeled hormone itself This result was interpreted to mean that challenger ligands bind allosteric sites on the LBD to induce hormone dissociation, and recent findings indicate that ligands bind weakly to multiple sites on the LBD surface. Here we show, that a large fraction of thyroid hormone receptor (TR) ligands promote rapid dissociation (T(1/2) < 2 h) of , radiolabeled T(3) vs. T(3) (T(1/2), approximate to 5-7 h). We cannot discern relationships between this effect and ligand size, activity or affinity for TR beta. One ligand, GC-24, binds the TR LBC and (weakly) to the TR beta-LBD surface that mediates dimer/heterodimer interaction, but we cannot link this interaction to rapid T(3) dissociation. Instead, several lines of evidence suggest that the challenger ligand must interact with the buried LBC to promote rapid T(3) release. Since previous molecular dynamics simulations suggest that TR ligands leave the LBC by several routes, we propose that a subset of challenger ligands binds and stabilizes a partially unfolded intermediate state of TR that arises during T(3) release and that this effect enhances hormone dissociation. (C) 2009 Elsevier Ltd. 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.

Structural basis of GC-1 selectivity for thyroid hormone receptor isoforms

Background: Thyroid receptors, TRa and TR beta, are involved in important physiological functions such as metabolism, cholesterol level and heart activities. Whereas metabolism increase and cholesterol level lowering could be achieved by TR beta isoform activation, TRa activation affects heart rates. Therefore, beta-selective thyromimetics have been developed as promising drug-candidates for treatment of obesity and elevated cholesterol level. GC-1 [ 3,5-dimethyl-4-(4'-hydroxy- 3'-isopropylbenzyl)-phenoxy acetic acid] has ability to lower LDL cholesterol with 600-to 1400-fold more potency and approximately two-to threefold more efficacy than atorvastatin (Lipitor(C)) in studies in rats, mice and monkeys. Results: To investigate GC-1 specificity, we solved crystal structures and performed molecular dynamics simulations of both isoforms complexed with GC-1. Crystal structures reveal that, in TRa Arg228 is observed in multiple conformations, an effect triggered by the differences in the interactions between GC-1 and Ser277 or the corresponding asparagine (Asn331) of TR beta. The corresponding Arg282 of TR beta is observed in only one single stable conformation, interacting effectively with the ligand. Molecular dynamics support this model: our simulations show that the multiple conformations can be observed for the Arg228 in TR alpha, in which the ligand interacts either strongly with the ligand or with the Ser277 residue. In contrast, a single stable Arg282 conformation is observed for TR beta, in which it strongly interacts with both GC-1 and the Asn331. Conclusion: Our analysis suggests that the key factors for GC-1 selectivity are the presence of an oxyacetic acid ester oxygen and the absence of the amino group relative to T(3). These results shed light into the beta-selectivity of GC-1 and may assist the development of new compounds with potential as drug candidates to the treatment of hypercholesterolemia and obesity.

Conformational differences between the wild type and V30M mutant transthyretin modulate its binding to genistein: Implications to tetramer stability and ligand-binding

Transthyretin (TTR) is a tetrameric beta-sheet-rich transporter protein directly involved in human amyloid diseases. It was recently found that the isoflavone genistein (GEN) potently inhibits TTR amyloid fibril formation (Green et al., 2005) and is therefore a promising candidate for TTR amyloidosis treatment. Here we used structural and biophysical approaches to characterize genistein binding to the wild type (TTRwt) and to its most frequent amyloidogenic variant, the V30M mutant. In a dose-dependent manner, genistein elicited considerable increases in both mutant and TTRwt stability as demonstrated by high hydrostatic pressure (HHP) and acid-mediated dissociation/denaturation assays. TTR:GEN crystal complexes and isothermal titration calorimetry (ITC) experiments showed that the binding mechanisms of genistein to the TTRwt and to V30M are different and are dependent on apoTTR structure conformations. Furthermore, we could also identify potential allosteric movements caused by genistein binding to the wild type TTR that explains, at least in part, the frequently observed negatively cooperative process between the two sites of TTRwt when binding ligands. These findings show that TTR mutants may present different ligand recognition and therefore are of value in ligand design for inhibiting TTR amyloidosis. (C) 2010 Elsevier Inc. All rights reserved.

Spectroscopic characterization and investigation of the dynamic of charge compensation process of supramolecular films derived from tetra-2-pyridyl-1,4-pyrazine ligand

This work describes the infrared spectroscopy characterization and the charge compensation dynamics in supramolecular film FeTPPZFeCN derived from tetra-2-pyridyl-1,4-pyrazine (TPPZ) with hexacyanoferrate, as well as the hybrid film formed by FeTPPZFeCN and polypyrrole (PPy). For supramolecular film, it was found that anion flux is greater in a K+ containing solution than in Li+ solution, which seems to be due to the larger crystalline ionic radius of K+. The electroneutralization process is discussed in terms of electrostatic interactions between cations and metallic centers in the hosting matrix. The nature of the charge compensation process differs from others modified electrodes based on Prussian blue films, where only cations such as K+ participate in the electroneutralization process. In the case of FeTPPZFeCN/PPy hybrid film, the magnitude of the anions’s flux is also dependent on the identity of the anion of the supporting electrolyte.

Physical signatures of discontinuities of the time-dependent exchange-correlation potential

The exact exchange-correlation (XC) potential in time-dependent density-functional theory (TDDFT) is known to develop steps and discontinuities upon change of the particle number in spatially confined regions or isolated subsystems. We demonstrate that the self-interaction corrected adiabatic local-density approximation for the XC potential has this property, using the example of electron loss of a model quantum well system. We then study the influence of the XC potential discontinuity in a real-time simulation of a dissociation process of an asymmetric double quantum well system, and show that it dramatically affects the population of the resulting isolated single quantum wells. This indicates the importance of a proper account of the discontinuities in TDDFT descriptions of ionization, dissociation or charge transfer processes.

Human Airway Epithelial Cell Responses to Neisseria lactamica and Purified Porin via Toll-Like Receptor 2-Dependent Signaling

The human airway epithelium is constantly exposed to microbial products from colonizing organisms. Regulation of Toll-like receptor (TLR) expression and specific interactions with bacterial ligands is thought to mitigate exacerbation of inflammatory processes induced by the commensal flora in these cells. The genus Neisseria comprises pathogenic and commensal organisms that colonize the human nasopharynx. Neisseria lactamica is not associated with disease, but N. meningitidis occasionally invades the host, causing meningococcal disease and septicemia. Upon colonization of the airway epithelium, specific host cell receptors interact with numerous Neisseria components, including the PorB porin, at the immediate bacterial-host cell interface. This major outer membrane protein is expressed by all Neisseria strains, regardless of pathogenicity, but its amino acid sequence varies among strains, particularly in the surface-exposed regions. The interaction of Neisseria PorB with TLR2 is essential for driving TLR2/TLR1-dependent cellular responses and is thought to occur via the porin`s surface-exposed loop regions. Our studies show that N. lactamica PorB is a TLR2 ligand but its binding specificity for TLR2 is different from that of meningococcal PorB. Furthermore, N. lactamica PorB is a poor inducer of proinflammatory mediators and of TLR2 expression in human airway epithelial cells. These effects are reproduced by whole N. lactamica organisms. Since the responsiveness of human airway epithelial cells to colonizing bacteria is in part regulated via TLR2 expression and signaling, commensal organisms such as N. lactamica would benefit from expressing a product that induces low TLR2-dependent local inflammation, likely delaying or avoiding clearance by the host.

Peroxisome proliferator-activated receptor-gamma ligand, 15-deoxy-Delta A(12,14)-prostaglandin J(2), reduces neutrophil migration via a nitric oxide pathway

Ligands for peroxisome proliferator-activated receptor gamma (PPAR-gamma), such as 15-deoxy-Delta(12,14)-prostaglandin J(2) (15d-PGJ(2)) have been implicated as a new class of anti-inflammatory compounds with possible clinical applications. Based on this concept, this investigation was designed to determine the effect of 15d-PGJ(2)-mediated activation of PPAR-gamma ligand on neutrophil migration after an inflammatory stimulus and clarify the underlying molecular mechanisms using a mouse model of peritonitis. Our results demonstrated that 15d-PGJ(2) administration decreases leukocyte rolling and adhesion to the inflammated mesenteric tissues by a mechanism dependent on NO. Specifically, pharmacological inhibitors of NO synthase remarkably abrogated the 15d-PGJ(2)-mediated suppression of neutrophil migration to the inflammatory site. Moreover, inducible NOS(-/-) mice were not susceptible to 15d-PGJ(2)-mediated suppression of neutrophil migration to the inflammatory sites when compared with their wild type. In addition, 15d-PGJ(2)-mediated suppression of neutrophil migration appeared to be independent of the production of cytokines and chemokines, since their production were not significantly affected in the carrageenan-injected peritoneal cavities. Finally, up-regulation of carrageenan-triggered ICAM-I expression in the mesenteric microcirculation vessels was abrogated by pretreatment of wild-type mice with 15d-PGJ(2), whereas 15d-PGJ(2) inhibited F-actin rearrangement process in neutrophils. Taken together these findings demonstrated that 15d-PGJ(2) suppresses inflammation-initiated neutrophil migration in a mechanism dependent on NO production in mesenteric tissues.

Loops, Chains, Sheets, and Networks from Variable Coordination of Cu(hfac)(2) with a Flexibly Hinged Aminoxyl Radical Ligand

One pair of reactants, Cu(hfac)(2) = M and the hinge-flexible radical ligand 5-(3-N-tert-butyl-N-aminoxylphenyl)pyrimidine (3PPN = L), yields a diverse set of five coordination complexes: a cyclic loop M(2)L(1) dimer; a 1:1 cocrystal between an M(2)L(2) loop and an ML(2) fragment; a ID chain of M(2)L(2) loops linked by M; two 2D M(3)L(2) networks of (M-L)(n) chains crosslinked by M with different repeat length pitches; a 3D M(3)L(2) network of M(2)L(2) loops cross-linking (M-L)(n)-type chains with connectivity different from those in the 2D networks. Most of the higher dimensional complexes exhibit reversible, temperature-dependent spin-state conversion of high-temperature paramagnetic states to lower magnetic moment states having antiferromagnetic exchange within Cu-ON bonds upon cooling, with accompanying bond contraction. The 3D complex also exhibited antiferromagnetic exchange between Cu(II) ions linked in chains through pyrimidine rings.

The binding of synthetic triiodo L-thyronine analogs to human transthyretin: Molecular basis of cooperative and non-cooperative ligand recognition

Transthyretin (TTR) is a tetrameric beta-sheet-rich transporter protein directly involved in human amyloid diseases. Several classes of small molecules can bind to TTR delaying its amyloid fibril formation, thus being promising drug candidates to treat TTR amyloidoses. In the present study, we characterized the interactions of the synthetic triiodo L-thyronine analogs and thyroid hormone nuclear receptor TR beta-selecfive agonists GC-1 and GC-24 with the wild type and V30M variant of human transthyretin (TTR). To achieve this aim, we conducted in vitro TTR acid-mediated aggregation and isothermal titration calorimetry experiments and determined the TTR:GC-1 and TTR:GC-24 crystal structures. Our data indicate that both GC-1 and GC-24 bind to TTR in a non-cooperative manner and are good inhibitors of TTR aggregation, with dissociation constants for both hormone binding sites (HBS) in the low micromolar range. Analysis of the crystal structures of TTRwt:GC-1(24) complexes and their comparison with the TTRwt X-ray structure bound to its natural ligand thyroxine (T4) suggests, at the molecular level, the basis for the cooperative process displayed by T4 and the non-cooperative process provoked by both GC-1 and GC-24 during binding to TTR. (C) 2010 Elsevier Inc. All rights reserved.

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.

Gaining ligand selectivity in thyroid hormone receptors via entropy

Nuclear receptors are important targets for pharmaceuticals, but similarities between family members cause difficulties in obtaining highly selective compounds. Synthetic ligands that are selective for thyroid hormone (TH) receptor beta (TR beta) vs. TR alpha reduce cholesterol and fat without effects on heart rate; thus, it is important to understand TR beta-selective binding. Binding of 3 selective ligands (GC-1, KB141, and GC-24) is characterized at the atomic level; preferential binding depends on a nonconserved residue (Asn-331 beta) in the TR beta ligand-binding cavity (LBC), and GC-24 gains extra selectivity from insertion of a bulky side group into an extension of the LBC that only opens up with this ligand. Here we report that the natural TH 3,5,3`-triodothyroacetic acid (Triac) exhibits a previously unrecognized mechanism of TR beta selectivity. TR x-ray structures reveal better fit of ligand with the TR alpha LBC. The TR beta LBC, however, expands relative to TR alpha in the presence of Triac (549 angstrom(3) vs. 461 angstrom(3)), and molecular dynamics simulations reveal that water occupies the extra space. Increased solvation compensates for weaker interactions of ligand with TR beta and permits greater flexibility of the Triac carboxylate group in TR beta than in TR alpha. We propose that this effect results in lower entropic restraint and decreases free energy of interactions between Triac and TR beta, explaining subtype-selective binding. Similar effects could potentially be exploited in nuclear receptor drug design.

Can mass dissociation patterns of transition-metal complexes be predicted from electrochemical data?

The Cooks kinetic method has been very convenient to correlate the relative dissociation rates obtained by collision-induced fragmentation experiments with the energies of two related bonds in molecules and complexes in the gas phase. Reliable bond energy data are, however, not always available, particularly for polynuclear transition-metal complexes, such as the triruthenium acetate clusters of the general formula [Ru(3) (mu(3)-O)(mu-CH(3)COO)(6)(py)(2)(L)](+), where L = ring substituted N-heterocyclic ligands. Accordingly, their gas-phase collision-induced tandem mass spectrometry (CID MS/MS) dissociation patterns have been analyzed pursuing a relationship with the more easily accessible redox potentials (E(1/2)) and Lever`s E(L) parameters. In fact, excellent linear correlations of In(1/2A(L)/A(py)), where A(py) and A(L) are the abundance of the fragments retaining the pyridine (py) and L ligand, respectively, with E(1/2) and E(L) were found. This result shows that those electrochemical parameters are correlated with bond energies and can be used in the analysis of the dissociation data. Such modified Cooks method can be used, for example, to determine the electronic effects of substituents on the metal-ligand bonds for a series of transition-metal complexes. Copyright (C) 2008 John Wiley & Sons, Ltd.

Oncostatin M is a novel glucocorticoid-dependent neuroinflammatory factor that enhances oligodendrocyte precursor cell activity in demyelinated sites

The innate immune reaction to tissue injury is a natural process, which can have detrimental effects in the absence of negative feedbacks by glucocorticoids (GCs). Although acute lipopolysaccharide (LPS) challenge is relatively harmless to the brain parenchyma of adult animals, the endotoxin is highly neurotoxic in animals that are treated with the GC receptor antagonist RU486. This study investigated the role of cytokines of the gp130-related family in these effects, because they are essential components of the inflammatory process that provide survival signals to neurons. Intracerebral LPS injection stimulated expression of several members of this family of cytokines, but oncostatin M (Osm) was the unique ligand to be completely inhibited by the RU486 treatment. OSM receptor (Osmr) is expressed mainly in astrocytes and endothelial cells following LPS administration and GCs are directly responsible for its transcriptional activation in the presence of the endotoxin. In a mouse model of demyelination, exogenous OSM significantly modulated the expression of genes involved in the mobilization of oligodendrocyte precursor cells (OPCs), differentiation of oligodendrocyte, and production of myelin. In conclusion, the activation of OSM signaling is a mechanism activated by TLR4 in the presence of negative feedback by GCs on the innate immune system of the brain. OSM absence is associated with detrimental effects of LPS, whereas exogenous OSM favors repair response to demyelinated regions. (C) 2010 Elsevier Inc. All rights reserved.