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.

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.

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.

Quantitative Analysis of Diepoxybutane Damage within the Chicken β-globin Domain

In industrial polymer and synthetic rubber production facilities, workers are exposed to 1,3-butadiene. This compound is converted in vivo to 1,2,3,4-diepoxybutane (DEB) and has been linked to increased incidences of cancer in these individuals. Carcinogenesis has been attributed to formation of DEB induced DNA interstrand cross-links. Previous studies have demonstrated that DEB cross-links deoxyguanosine residues within 5'-GNC sequences in synthetic DNA, in restriction fragments, and in defined sequence nucleosomes. The current study utilized the polymerase chain reaction (PCR) to examine DEB damage frequencies within nuclear genes, found within "open" regions of chromatin, as compared to regions of unexpressed sequence that reside in tightly packed, "closed" chromatin, to more closely model DEB reactivity in vivo. These initial studies have been performed in chicken liver homogenates. Preliminarily, we have found a dose-dependent DEB lesion-forming response within "open" chromatin. DEB appears to have little-to-no effect upon regions of "closed" chromatin.

Computer-aided design of a novel ligand for retinoic acid receptor in cancer chemotherapy

The isotypes of RAR and RXR are retinoic acid and retinoid X acid receptors, respectively, whose ligand-binding domain contains the ligand-dependent activation function, with distinct pharmacological targets for retinoids, involved in the treatment of various cancers and skin diseases. Due to the major challenge which cancer treatment and cure still imposes after many decades to the international scientific community, there is actually considerable interest in new ligands with increased bioactivity. We have focused on the retinoid acid receptor, which is considered an interesting target for drug design. In this work, we carried out density functional geometry optimizations, and different docking procedures. We performed screening in a large database (hundreds of thousands of molecules which we optimized at the AM1 level) yielding a set of potential bioactive ligands. A new ligand was selected and optimized at the B3LYP/6-31G* level. A flexible docking program was used to investigate the interactions between the receptor and the new ligand. The result of this work is compared with several crystallographic ligands of RAR. Our theoretically more bioactive new-ligand indicates stronger and more hydrogen bonds as well as hydrophobic interactions with the receptor. (c) 2005 Wiley Periodicals, Inc.

Sequence, evolution and ligand binding properties of mammalian Duffy antigen/receptor for chemokines

The Duffy antigen/receptor for chemokine, DARC, acts as a widely expressed promiscuous chemokine receptor and as the erythrocyte receptor for Plasmodium vivax. To gain insight into the evolution and structure/function relations of DARC, we analyzed the binding of anti-human Fy monoclonal antibodies (mAbs) and human chemokines to red blood cells (RBCs) from 11 nonhuman primates and two nonprimate mammals, and we elucidated the structures of the DARC genes from gorilla, gibbon, baboon, marmoset, tamarin, night monkey and cattle. CXCL-8 and CCL-5 chemokine binding analysis indicated that the promiscuous binding profile characteristic of DARC is conserved across species. Among three mAbs that detected the Fy6 epitope by flow cytometric analysis of human and chimpanzee RBCs, only one reacted with night monkey and squirrel monkey. Only chimpanzee RBCs bound a significant amount of the anti-Fy3 mAb. Fy3 was also poorly detected on RBCs from gorilla, baboon and rhesus monkey, but not from new world monkeys. Alignment of DARC homologous sequences allowed us to construct a phylogenetic tree in which all branchings were in accordance with current knowledge of primate phylogeny. Although DARC was expected to be under strong internal and external selection pressure, in order to maintain chemokine binding and avoid Plasmodium vivax binding, respectively, our present study did not provide arguments in favor of a selection pressure on the extracellular domains involved in ligand specificity. The amino acid variability of DARC-like polypeptides was found to be well correlated with the hydrophylicity indexes, with the highest divergence on the amino-terminal extracellular domain. Analysis of the deduced amino acid sequences highlighted the conservation of some amino acid residues, which should prove to be critical for the structural and functional properties of DARC.

Alteration of sFAS and sFAS ligand expression during canine visceral leishmaniosis

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Characterization of a novel five-transmembrane domain cholecystokinin-2 receptor splice variant identified in human tumors

The cholecystokinin-2 receptor (CCK2R), is expressed in cancers where it contributes to tumor progression. The CCK2R is over-expressed in a sub-set of tumors, allowing its use in tumor targeting with a radiolabel ligand. Since discrepancies between mRNA levels and CCK2R binding sites were noticed, we searched for abnormally spliced variants in tumors from various origins having been previously reported to frequently express cholecystokinin receptors, such as medullary thyroid carcinomas, gastrointestinal stromal tumors, leiomyomas and leiomyosarcomas, and gastroenteropancreatic tumors. A variant of the CCK2R coding for a putative five-transmembrane domains receptor has been cloned. This variant represented as much as 6% of CCK2R levels. Ectopic expression in COS-7 cells revealed that this variant lacks biological activity due to its sequestration in endoplasmic reticulum. When co-expressed with the CCK2R, this variant diminished membrane density of the CCK2R and CCK2R-mediated activity (phospholipase-C and ERK activation). In conclusion, a novel splice variant acting as a dominant negative on membrane density of the CCK2R may be of importance for the pathophysiology of certain tumors and for their in vivo CCK2R-targeting.

Posttranscriptional regulation of Fas (CD95) ligand killing activity by lipid rafts

Fas (CD95/Apo-1) ligand-mediated apoptosis induction of target cells is one of the major effector mechanisms by which cytotoxic lymphocytes (T cells and natural killer cells) kill their target cells. In T cells, Fas ligand expression is tightly regulated at a transcriptional level through the activation of a distinct set of transcription factors. Increasing evidence, however, supports an important role for posttranscriptional regulation of Fas ligand expression and activity. Lipid rafts are cholesterol- and sphingolipid-rich membrane microdomains, critically involved in the regulation of membrane receptor signaling complexes through the clustering and concentration of signaling molecules. Here, we now provide evidence that Fas ligand is constitutively localized in lipid rafts of FasL transfectants and primary T cells. Importantly, disruption of lipid rafts strongly reduces the apoptosis-inducing activity of Fas ligand. Localization to lipid rafts appears to be predominantly mediated by the characteristic cytoplasmic proline-rich domain of Fas ligand because mutations of this domain result in reduced recruitment to lipid rafts and attenuated Fas ligand killing activity. We conclude that Fas ligand clustering in lipid rafts represents an important control mechanism in the regulation of T cell-mediated cytotoxicity.

Crystal structure of the platelet glycoprotein Ib(alpha) N-terminal domain reveals an unmasking mechanism for receptor activation

Glycoprotein Ib (GPIb) is a platelet receptor with a critical role in mediating the arrest of platelets at sites of vascular damage. GPIb binds to the A1 domain of von Willebrand factor (vWF-A1) at high blood shear, initiating platelet adhesion and contributing to the formation of a thrombus. To investigate the molecular basis of GPIb regulation and ligand binding, we have determined the structure of the N-terminal domain of the GPIb(alpha) chain (residues 1-279). This structure is the first determined from the cell adhesion/signaling class of leucine-rich repeat (LRR) proteins and reveals the topology of the characteristic disulfide-bonded flanking regions. The fold consists of an N-terminal beta-hairpin, eight leucine-rich repeats, a disulfide-bonded loop, and a C-terminal anionic region. The structure also demonstrates a novel LRR motif in the form of an M-shaped arrangement of three tandem beta-turns. Negatively charged binding surfaces on the LRR concave face and anionic region indicate two-step binding kinetics to vWF-A1, which can be regulated by an unmasking mechanism involving conformational change of a key loop. Using molecular docking of the GPIb and vWF-A1 crystal structures, we were also able to model the GPIb.vWF-A1 complex.

Ligand binding induces Cbl-dependent EphB1 receptor degradation through the lysosomal pathway

Eph receptor tyrosine kinases play a critical role in embryonic patterning and angiogenesis. In the adult, they are involved in carcinogenesis and pathological neovascularization. However, the mechanisms underlying their role in tumor formation and metastasis remain to be defined. Here, we demonstrated that stimulation of EphB1 with ephrinB1/Fc led to a marked downregulation of EphB1 protein, a process blocked by the lysosomal inhibitor bafilomycin. Following ephrinB1 stimulation, the ubiquitin ligase Cbl was recruited by EphB1 and then phosphorylated. Both Cbl phosphorylation and EphB1 ubiquitination were blocked by the Src inhibitor PP2. Overexpression of wild-type Cbl, but not of 70Z mutant lacking ligase activity, enhanced EphB1 ubiquitination and degradation. This negative regulation required the tyrosine kinase activity of EphB1 as kinase-dead EphB1-K652R was resistant to Cbl. Glutathione S-transferase binding experiments showed that Cbl bound to EphB1 through its tyrosine kinase-binding domain. In aggregate, we demonstrated that Cbl induces the ubiquitination and lysosomal degradation of activated EphB1, a process requiring EphB1 and Src kinase activity. To our knowledge, this is the first study dissecting the molecular mechanisms leading to EphB1 downregulation, thus paving the way to new means of modulating their angiogenic and tumorigenic properties.

Nck recruitment to Eph receptor, EphB1/ELK, couples ligand activation to c-Jun kinase.

Eph family receptor tyrosine kinases signal axonal guidance, neuronal bundling, and angiogenesis; yet the signaling systems that couple these receptors to targeting and cell-cell assembly responses are incompletely defined. Functional links to regulators of cytoskeletal structure are anticipated based on receptor mediated cell-cell aggregation and migratory responses. We used two-hybrid interaction cloning to identify EphB1-interactive proteins. Six independent cDNAs encoding the SH2 domain of the adapter protein, Nck, were recovered in a screen of a murine embryonic library. We mapped the EphB1 subdomain that binds Nck and its Drosophila homologue, DOCK, to the juxtamembrane region. Within this subdomain, Tyr594 was required for Nck binding. In P19 embryonal carcinoma cells, activation of EphB1 (ELK) by its ligand, ephrin-B1/Fc, recruited Nck to native receptor complexes and activated c-Jun kinase (JNK/SAPK). Transient overexpression of mutant EphB1 receptors (Y594F) blocked Nck recruitment to EphB1, attenuated downstream JNK activation, and blocked cell attachment responses. These findings identify Nck as an important intermediary linking EphB1 signaling to JNK.

Lethal giant discs, a novel C2 domain protein, restricts Notch activation during endocytosis

The Notch signaling pathway plays a central role in metazoan growth and patterning, and its deregulation leads to many human diseases, including cancer. It is therefore important to understand the modes of Notch signaling regulation. Recent discoveries have demonstrated that mutations in conserved endosomal pathway components such as Erupted and Vps25 can ectopically activate Notch signaling in Drosophila. Mutations in the tumor suppressor lethal giant discs (lgd) display similar but even stronger and more specific Notch activation than in the erupted and vps25 mutant animals. This Notch activation in lgd mutant tissues causes hyperplastic overgrowth of the Drosophila imaginal discs, and the eventual lethality of the animal. However, the gene that encodes Lgd, and its function in the Notch pathway have not yet been identified. ^ I have found that Lgd is a novel, conserved C2 domain protein that regulates Notch trafficking. Lgd cell-autonomously restricts Notch signaling in the Drosophila wing disc to the target cells in the D/V boundary. The function of Lgd lies at or upstream of Notch S3 activation, but Lgd doesn't affect the binding affinities between Notch and Delta. Lgd is also not required for cis-inhibition of Notch signaling by ligands. Notch accumulates on the early endosome in lgd mutant cells and signals in a ligand-independent manner, a result that has previously been seen in endosomal pathway mutants. Interestingly, Notch activation in lgd mutant cells is dependent on the endosomal protein Hrs, and Lgd activity appears to be downstream of Hrs function in endocytosis. Taken together, my data identify Lgd as a novel tumor suppressor protein that regulates Notch signaling by targeting Notch for degradation or recycling. ^

Ligand-stimulated internalization of a Dictyostelium discoideum G protein-coupled cAMP receptor

The social amoeba, Dictyostelium discoideum, undergoes a remarkable starvation-induced program of development that transforms a population of unicellular amoebae into a fruiting body composed of resistant spores suspended on a stalk. During this development, secreted cAMP drives chemotaxis of the amoebae, leading to their aggregation, and subsequent differentiation and morphogenesis. Four sequentially expressed G protein-coupled receptors (GPCRs) for cAMP play critical roles in this process. The first of these, cAR1, is essential for aggregation as it mediates chemotaxis as well as the propagation of secreted cAMP waves throughout aggregating populations. Ligand-induced internalization has been shown to regulate a variety of GPCRs. However, little was known at the outset of this study about the role of internalization in the regulation of cAR1 function or, for that matter, in developmental systems in general. For this study, cAMP-induced cAR1 internalization was assessed by measuring (1) the reduction of cell surface binding sites for [ 3H]cAMP and (2) the redistribution of YFP-tagged receptors to the cell's interior, cAMP was found to induce little or no loss of ligand binding (LLB) in vegetative cells. However, the ability to induce LLB increased progressively over the initial 6 hrs of development, reaching ∼70% in cells undergoing aggregation. Despite these reductions in surface binding, detectable cAR1-YFP redistribution could be induced by cAMP only after the cells reached the mound stage (10 hrs) and was found to occur naturally by the ensuing slug stage (18 hrs). Site-directed substitution of a cluster of 5 serines in the receptor's cytoplasmic tail that was previously shown to be the principal site of cAMP-induced cAR1 phosphorylation impaired both LLB and receptor redistribution and furthermore resulted in mound-stage developmental arrest, suggesting that phosphorylation of cAR1 is a prerequisite for its internalization and that cAR1 internalization is required for post-aggregative development. To assess the involvement of clathrin mediated endocytosis, Dictyostelium cells lacking the clathrin light chain gene (clc-) or either of two dynamin genes were examined and found to be defective in LLB and, in the case of clc- cells, also cAR1 redistribution and turnover. Furthermore, cAR1 overexpression in clc- cells (like the serine mutant in wild-type cells) promoted developmental arrest in mounds. The mound-arrest phenotype was also recapitulated in a wild-type background by the specific expression of cAR1 in prestalk cells (but not prespore cells), suggesting that development depends critically on internalization and clearance of cAR1 from these cells. Persistent cAR1 expression following aggregation was found to be associated with aberrant expression of prestalk and prespore genes, which may adversely affect development in the prestalk cell lineage. The PI3 kinase-TORC2 signal transduction pathway, known to be important for Dictyostelium chemotaxis and internalization of yeast pheromone receptors, was examined using chemical inhibitors and null cells and found to be necessary for cAR1 internalization. In conclusion, cAR1 was shown to be similar to other GPCRs in that its internalization depends on phosphorylation of cytoplasmic domain serines, utilizes clathrin and dynamin, and involves the TORC2 complex. In addition, the findings presented here that cAR1 internalization is both developmentally regulated and required for normal development represent a novel regulatory paradigm that might pertain to other GPCRs known to play important roles in the development of humans and other metazoans. ^

Heregulin -beta 1 directly associates with RhoA via the immunoglobulin -like domain thereby modulating breast cancer metastasis-related properties

Heregulins constitute a family of growth factors belonging to the epidermal growth factor (EGF) family. Breast cancers that overexpress specific members of the EGF receptor family (EGFR, ErbB2, ErbB3, ErbB4) have increased metastatic potential, and Heregulin-β1 (HRGβ1), a ligand for ErbB3 and ErbB4, has also been shown to induce metastasis-related properties in breast cancer cells in vitro. The secreted form of the HRGβ1 is composed of five distinct structural domains, including the N-terminal domain, an immunoglobulin-like domain (IgG-like), a glycosylation domain, an EGF-like domain, and a β1-specific domain. Of these, the EGF-like domain is well characterized for its function in metastasis-related properties as well as its structure. However, the contributions of the other HRGβ1 domains in breast cancer metastasis remains unclear. ^ To investigate this, HRGβ1 proteins with targeted domain deletions were purified and subjected to assays for metastasis-related properties, including aggregation, invasion, activation of EGFR family members, and motility of breast cancer cells. These assays showed that retaining the EGF-like domain of HRGβ1 is important for activation of EGFRs. Interestingly, the HRGβ1 protein lacking the IgG-like domain (NGEB) led to a decrease in breast cancer cell motility, indicating the IgG-like domain modulates cell motility, an important step in cancer metastasis. ^ To understand the underlying mechanisms, I performed protein sequence and structural analysis of HRGβ1 and identified that the IgG-like domain of HRGβ1 shares sequence homology and three-dimensional structural similarity with the IgG-like domain of TRIO. TRIO is a cytoplasmic protein that directly associates with RhoA, a GTPase involved in cell reorganization and cell motility. Therefore, I hypothesized that HRGβ1 may translocate inside the breast cancer cells through receptor mediated endocytosis and bind to RhoA via its IgG-like domain. I show wild type HRGβ1 but not NGEB binds RhoA in vitro and in vivo, leading to RhoA activation. Inhibition of HRG-β1 internalization via endocytosis disrupted HRGβ1 binding to RhoA. Additionally, breast cancer cell motility induced by HRG-β1 is reduced after treatment with inhibitors to both endocytosis and RhoA function, similar to levels seen with NGEB treatment. ^ Thus, in addition to the well-known role of HRGβ1 as an extracellular stimulator of the EGFR family members, HRGβ1 also functions within the cell as a binding partner and activator of RhoA to modulate cancer cell motility. ^