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.

Crystallographic comparison of the estrogen and progesterone receptor’s ligand binding domains

The 2.8-Å crystal structure of the complex formed by estradiol and the human estrogen receptor-α ligand binding domain (hERαLBD) is described and compared with the recently reported structure of the progesterone complex of the human progesterone receptor ligand binding domain, as well as with similar structures of steroid/nuclear receptor LBDs solved elsewhere. The hormone-bound hERαLBD forms a distinctly different and probably more physiologically important dimer interface than its progesterone counterpart. A comparison of the specificity determinants of hormone binding reveals a common structural theme of mutually supported van der Waals and hydrogen-bonded interactions involving highly conserved residues. The previously suggested mechanism by which the estrogen receptor distinguishes estradiol’s unique 3-hydroxy group from the 3-keto function of most other steroids is now described in atomic detail. Mapping of mutagenesis results points to a coactivator-binding surface that includes the region around the “signature sequence” as well as helix 12, where the ligand-dependent conformation of the activation function 2 core is similar in all previously solved steroid/nuclear receptor LBDs. A peculiar crystal packing event displaces helix 12 in the hERαLBD reported here, suggesting a higher degree of dynamic variability than expected for this critical substructure.

NMR and mutagenesis evidence for an I domain allosteric site that regulates lymphocyte function-associated antigen 1 ligand binding

The leukocyte integrin, lymphocyte function-associated antigen 1 (LFA-1) (CD11a/CD18), mediates cell adhesion and signaling in inflammatory and immune responses. To support these functions, LFA-1 must convert from a resting to an activated state that avidly binds its ligands such as intercellular adhesion molecule 1 (ICAM-1). Biochemical and x-ray studies of the Mac-1 (CD11b/CD18) I domain suggest that integrin activation could involve a conformational change of the C-terminal α-helix. We report the use of NMR spectroscopy to identify CD11a I domain residues whose resonances are affected by binding to ICAM-1. We observed two distinct sites in the CD11a I domain that were affected. As expected from previous mutagenesis studies, a cluster of residues localized around the metal ion-dependent adhesion site (MIDAS) was severely perturbed on ICAM-1 binding. A second cluster of residues distal to the MIDAS that included the C-terminal α-helix of the CD11a I domain was also affected. Substitution of residues in the core of this second I domain site resulted in constitutively active LFA-1 binding to ICAM-1. Binding data indicates that none of the 20 substitution mutants we tested at this second site form an essential ICAM-1 binding interface. We also demonstrate that residues in the I domain linker sequences can regulate LFA-1 binding. These results indicate that LFA-1 binding to ICAM-1 is regulated by an I domain allosteric site (IDAS) and that this site is structurally linked to the MIDAS.

Calcium is essential for the structural integrity of the cysteine-rich, ligand-binding repeat of the low-density lipoprotein receptor

Seven cysteine-rich repeats form the ligand-binding region of the low-density lipoprotein (LDL) receptor. Each of these repeats is assumed to bind a calcium ion, which is needed for association of the receptor with its ligands, LDL and beta-VLDL. The effects of metal ions on the folding of the reduced N-terminal cysteine-rich repeat have been examined by using reverse-phase high-performance liquid chromatography to follow the formation of fully oxidized isomers with different disulfide connectivities. in the absence of calcium many of the 15 possible isomers formed on oxidation, whereas in its presence the predominant product at equilibrium had the native disulfide bond connectivities. Other metals were far less effective at directing disulfide bond formation: Mn2+ partly mimicked the action of Ca2+, but Ba2+, Sr2+, and Mg2+ had little effect. This metal-ion specificity was also observed in two-dimensional H-1 NMR spectral studies: only Ca2+ induced the native three-dimensional fold. The two paramagnetic ions, Gd3+ and Mn2+, and Cd2+ did not promote adoption of a well-defined structure, and the two paramagnetic ions did not displace calcium ions. The location of calcium ion binding sites in the repeat was also explored by NMR spectroscopy. The absence of chemical shift changes for the side chain proton resonances of Asp26, Asp36, and Glu37 from pH 3.9 to 6.8 in the presence of calcium ions and their proximal location in the NMR structures implicated these side chains as calcium ligands. Deuterium exchange NMR experiments also revealed a network of hydrogen bonds that stabilizes the putative calcium-binding loop.

Differential roles of T cell receptor alpha and beta chains in ligand binding among H-2Kd-restricted cytolytic T lymphocyte clones specific for a photoreactive Plasmodium berghei circumsporozoite peptide derivative.

To study the interaction of T cell receptor with its ligand, a complex of a major histocompatibility complex molecule and a peptide, we derived H-2Kd-restricted cytolytic T lymphocyte clones from mice immunized with a Plasmodium berghei circumsporozoite peptide (PbCS) 252-260 (SYIPSAEKI) derivative containing photoreactive Nepsilon-[4-azidobenzoyl] lysine in place of Pro-255. This residue and Lys-259 were essential parts of the epitope recognized by these clones. Most of the clones expressed BV1S1A1 encoded beta chains along with specific complementary determining region (CDR) 3beta regions but diverse alpha chain sequences. Surprisingly, all T cell receptors were preferentially photoaffinity labeled on the alpha chain. For a representative T cell receptor, the photoaffinity labeled site was located in the Valpha C-strand. Computer modeling suggested the presence of a hydrophobic pocket, which is formed by parts of the Valpha/Jalpha C-, F-, and G-strands and adjacent CDR3alpha residues and structured to be able to avidly bind the photoreactive ligand side chain. We previously found that a T cell receptor specific for a PbCS peptide derivative containing this photoreactive side chain in position 259 similarly used a hydrophobic pocket located between the junctional CDR3 loops. We propose that this nonpolar domain in these locations allow T cell receptors to avidly and specifically bind epitopes containing non-peptidic side chains.

The gene for the ligand binding chain of the human interferon gamma receptor.

In order to characterize the gene encoding the ligand binding (1(st); alpha) chain of the human IFN-gamma receptor, two overlapping cosmid clones were analyzed. The gene spans over 25 kilobases (kb) of the genomic DNA and has seven exons. The extracellular domain is encoded by exons 1 to 5 and by part of exon 6. The transmembrane region is also encoded by exon 6. Exon 7 encodes the intracellular domain and the 3' untranslated portion. The gene was located on chromosome 6q23.1, as determined by in situ hybridization. The 4 kb region upstream (5') of the gene was sequenced and analyzed for promoter activity. No consensus-matching TATA or CAAT boxes in the 5' region were found. Potential binding sites for Sp1, AP-1, AP-2, and CREB nuclear factors were identified. Compatible with the presence of the Sp1/AP-2 sites and the lack of TATA box, S1-nuclease mapping experiments showed multiple transcription initiation sites. Promoter activity of the 5' flanking region was analyzed with two different reporter genes: the Escherichia coli chloramphenicol acetyltransferase and human growth hormone. The smallest 5' region of the gene that still had full promoter activity was 692 base pairs in length. In addition, we found sequences belonging to the oldest family of Alu repeats, 2 - 3 kb upstream of the gene, which could be useful for genetic studies.

Alpha2-adrenoceptors: Structure and ligand binding properties at the molecular level

Alpha₂-Adrenoceptors: structure and ligand binding properties at the molecular level The mouse is the most frequently used animal model in biomedical research, but the use of zebrafish as a model organism to mimic human diseases is on the increase. Therefore it is considered important to understand their pharmacological differences from humans also at the molecular level. The zebrafish Alpha_2-adrenoceptors were expressed in mammalian cells and the binding affinities of 20 diverse ligands were determined and compared to the corresponding human receptors. The pharmacological properties of the human and zebrafish Alpha_2--adrenoceptors were found to be quite well conserved. Receptor models based on the crystal structures of bovine rhodopsin and the human Beta₂-adrenoceptor revealed that most structural differences between the paralogous and orthologous Alpha_2--adrenoceptors were located within the second extracellular loop (XL2). Reciprocal mutations were generated in the mouse and human Alpha_2--adrenoceptors. Ligand binding experiments revealed that substitutions in XL2 reversed the binding profiles of the human and mouse Alpha_2--adrenoceptors for yohimbine, rauwolscine and RS-79948-197, evidence for a role for XL2 in the determination of species-specific ligand binding. Previous mutagenesis studies had not been able to explain the subtype preference of several large Alpha_2--adrenoceptor antagonists. We prepared chimaeric Alpha_2--adrenoceptors where the first transmembrane (TM1) domain was exchanged between the three human Alpha_2--adrenoceptor subtypes. The binding affinities of spiperone, spiroxatrine and chlorpromazine were observed to be significantly improved by TM1 substitutions of the Alpha_2a--adrenoceptor. Docking simulations indicated that indirect effects, such as allosteric modulation, are more likely to be involved in this phenomenon rather than specific side-chain interactions between ligands and receptors.

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.

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.

Multiple loop structures critical for ligand binding of the integrin α4 subunit in the upper face of the β-propeller mode 1

A non-I-domain integrin, α4β1, recognizes vascular cell adhesion molecule 1 (VCAM-1) and the IIICS portion of fibronectin. To localize regions of α4 critical for ligand binding, we swapped several predicted loops within or near the putative ligand-binding site of α4 (which spans repeats 2–5 of the seven N-terminal repeats) with the corresponding regions of α5. Swapping residues 112–131 in repeat 2, or residues 237–247 in repeat 4, completely blocked adhesion to immobilized VCAM-1 and connecting segment 1 (CS-1) peptide. However, swapping residues 40–52 in repeat 1, residues 151–164 in repeat 3, or residues 282–288 (which contain a putative cation binding motif) in repeat 5 did not affect or only slightly reduced adhesion to these ligands. The binding of several function-blocking antibodies is blocked by swapping residues 112–131, 151–164, and 186–191 (which contain previously identified residues critical for ligand binding, Tyr-187 and Gly-190). These results are consistent with the recently published β-propeller folding model of the integrin α4 subunit [Springer, T. A. (1997) Proc. Natl. Acad. Sci. USA 94, 65–72], in which seven four-stranded β-sheets are arranged in a torus around a pseudosymmetric axis. The regions of α4 critical for ligand binding are adjacent to each other and are located in the upper face, the predicted ligand-binding site, of the β-propeller model, although they are not adjacent in the primary structure.

High-affinity binding of the Drosophila Numb phosphotyrosine-binding domain to peptides containing a Gly-Pro-(p)Tyr motif

The phosphotyrosine-binding (PTB) domain is a recently identified protein module that has been characterized as binding to phosphopeptides containing an NPXpY motif (X = any amino acid). We describe here a novel peptide sequence recognized by the PTB domain from Drosophila Numb (dNumb), a protein involved in cell fate determination and asymmetric cell division during the development of the Drosophila nervous system. Using a Tyr-oriented peptide library to screen for ligands, the dNumb PTB domain was found to bind selectively to peptides containing a YIGPYφ motif (φ represents a hydrophobic residue). A synthetic peptide containing this sequence bound specifically to the isolated dNumb PTB domain in solution with a dissociation constant (Kd) of 5.78 ± 0.74 μM. Interestingly, the affinity of this peptide for the dNumb PTB domain was increased (Kd = 1.41 ± 0.10 μM) when the second tyrosine in the sequence was phosphorylated. Amino acid substitution studies of the phosphopeptide demonstrated that a core motif of sequence GP(p)Y is required for high-affinity binding to the dNumb PTB domain. Nuclear magnetic resonance experiments performed on isotopically labeled protein complexed with either Tyr- or pTyr-containing peptides suggest that the same set of amino acids in the dNumb PTB domain is involved in binding both phosphorylated and nonphosphorylated forms of the peptide. The in vitro selectivity of the dNumb PTB domain is therefore markedly different from those of the Shc and IRS-1 PTB domains, in that it interacts preferentially with a GP(p)Y motif, rather than NPXpY, and does not absolutely require ligand phosphorylation for binding. Our results suggest that the PTB domain is a versatile protein module, capable of exhibiting varied binding specificities.

Mapping rat megalin: the second cluster of ligand binding repeats contains a 46-amino acid pathogenic epitope involved in the formation of immune deposits in Heymann nephritis.

Megalin (gp330), an epithelial endocytic receptor, is a major target antigen of Heymann nephritis (HN), an autoimmune disease in rats. To elucidate the mechanisms of HN, we have mapped a pathogenic epitope in megalin that binds anti-megalin antibodies. We focused our attention on four clusters of cysteine-rich, low density lipoprotein receptor (LDLR) ligand binding repeats in the extracellular domain of megalin because they represent putative ligand binding regions and therefore would be expected to be exposed in vivo and to be able to bind circulating antibodies. Rat megalin cDNA fragments I through IV encoding the first through fourth clusters of ligand-binding repeats, respectively, were expressed in a baculovirus system. All four expression products were detected by immunoblotting with two antisera capable of inducing passive HN (pHN). When antibodies eluted from glomeruli of rats with pHN were used for immunoblotting, only the expression product encoded by fragment II was detected. This indicates that the second cluster of LDLR ligand binding repeats is directly involved in binding anti-megalin antibodies and in the induction of pHN. To narrow the major epitope in this domain, fragment II was used to prepare proteins sequentially truncated from the C- and N-terminal ends by in vitro translation. Analysis of the truncated translation products by immunoprecipitation with anti-megalin IgG revealed that the fifth ligand-binding repeat (amino acids 1160-1205) contains the major epitope recognized. This suggests that a 46-amino acid sequence in the second cluster of LDLR ligand binding repeats contains a major pathogenic epitope that plays a key role in pHN. Identification of this epitope will facilitate studies on the pathogenesis of HN.

An inhalational anesthetic binding domain in the nicotinic acetylcholine receptor.

To determine inhalational anesthetic binding domains on a ligand-gated ion channel, I used halothane direct photoaffinity labeling of the nicotinic acetylcholine receptor (nAChR) in native Torpedo membranes. [14C]Halothane photoaffinity labeling of both the native Torpedo membranes and the isolated nAChR was saturable, with Kd values within the clinically relevant range. All phospholipids were labeled, with greater than 95% of the label in the acyl chain region. Electrophoresis of labeled nAChR demonstrated no significant subunit selectivity for halothane incorporation. Within the alpha-subunit, greater than 90% of label was found in the endoprotease Glu-C digestion fragments which contain the four transmembrane regions, and the pattern was different from that reported for photoactivatable phospholipid binding to the nAChR. Unlabeled halothane reduced labeling more than did isoflurane, suggesting differences in the binding domains for inhalational anesthetics in the nAChR. These data suggest multiple similar binding domains for halothane in the transmembrane region of the nAChR.

Identification of residues that control specific binding of the Shc phosphotyrosine-binding domain to phosphotyrosine sites.

The Shc adaptor protein contains two phosphotyrosine [Tyr(P)]binding modules--an N-terminal Tyr(P) binding (PTB) domain and a C-terminal Src homology 2 (SH2) domain. We have compared the ability of the Shc PTB domain to bind the receptors for nerve growth factor and insulin, both of which contain juxtamembrane Asn-Pro-Xaa-Tyr(P) motifs implicated in PTB binding. The Shc PTB domain binds with high affinity to a phosphopeptide corresponding to the nerve growth factor receptor Tyr-490 autophosphorylation site. Analysis of individual residues within this motif indicates that the Asn at position -3 [with respect to Tyr(P)], in addition to Tyr(P), is critical for PTB binding, while the Pro at position -2 plays a less significant role. A hydrophobic amino acid 5 residues N-terminal to the Tyr(P) is also essential for high-affinity binding. In contrast, the Shc PTB domain does not bind stably to the Asn-Pro-Xaa-Tyr(P) site at Tyr-960 in the activated insulin receptor, which has a polar residue (Ser) at position -5. Substitution of this Ser at position -5 with Ile markedly increased binding of the insulin receptor Tyr-960 phosphopeptide to the PTB domain. These results suggest that while the Shc PTB domain recognizes a core sequence of Asn-Pro-Xaa-Tyr(P), its binding affinity is modulated by more N-terminal residues in the ligand, which therefore contribute to the specificity of PTB-receptor interactions. An analysis of residues in the Shc PTB domain required for binding to Tyr(P) sites identified a specific and evolutionarily conserved Arg (Arg-175) that is uniquely important for ligand binding and is potentially involved in Tyr(P) recognition.