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.

NMR structure and backbone dynamics of a concatemer of epidermal growth factor homology modules of the human low-density lipoprotein receptor

The ligand-binding region of the low-density lipoprotein (LDL) receptor is formed by seven N-terminal, imperfect, cysteine-rich (LB) modules. This segment is followed by an epidermal growth factor precursor homology domain with two N-terminal, tandem, EGF-like modules that are thought to participate in LDL binding and recycling of the endocytosed receptor to the cell surface. EGF-A and the concatemer, EGF-AB, of these modules were expressed in Escherichia coli. Correct protein folding of EGF-A and the concatemer EGF-AB was achieved in the presence or absence of calcium ions, in contrast to the LB modules, which require them for correct folding. Homonuclear and heteronuclear H-1-N-15 NMR spectroscopy at 17.6 T was used to determine the three-dimensional structure of the concatemer. Both modules are formed by two pairs of short, anti-parallel beta -strands. In the concatemer, these modules have a fixed relative orientation, stabilized by calcium ion-binding and hydrophobic interactions at the interface. N-15 longitudinal and transverse relaxation rates, and {H-1}-N-15 heteronuclear NOEs were used to derive a model-free description of the backbone dynamics of the molecule. The concatemer appears relatively rigid, particularly near the calcium ion-binding site at the module interface, with an average generalized order parameter of 0.85 +/- 0.11. Some mutations causing familial hypercholesterolemia may now be rationalized. Mutations of D41, D43 and E44 in the EGF-B calcium ion-binding region may affect the stability of the linker and thus the orientation of the tandem modules. The diminutive core also provides little structural stabilization, necessitating the presence of disulfide bonds. The structure and dynamics of EGF-AB contrast with the N-terminal LB modules, which require calcium ions both for folding to form the correct disulfide connectivities and for maintenance of the folded structure, and are connected by highly mobile linking peptides. (C) 2001 Academic Press.

Full activation of the T cell receptor requires both clustering and conformational changes at CD3.

T cell receptor (TCR-CD3) triggering involves both receptor clustering and conformational changes at the cytoplasmic tails of the CD3 subunits. The mechanism by which TCRalphabeta ligand binding confers conformational changes to CD3 is unknown. By using well-defined ligands, we showed that induction of the conformational change requires both multivalent engagement and the mobility restriction of the TCR-CD3 imposed by the plasma membrane. The conformational change is elicited by cooperative rearrangements of two TCR-CD3 complexes and does not require accompanying changes in the structure of the TCRalphabeta ectodomains. This conformational change at CD3 reverts upon ligand dissociation and is required for T cell activation. Thus, our permissive geometry model provides a molecular mechanism that rationalizes how the information of ligand binding to TCRalphabeta is transmitted to the CD3 subunits and to the intracellular signaling machinery.

Cooperative enhancement of specificity in a lattice of T cell receptors

Two of the most important models to account for the specificity and sensitivity of the T cell receptor (TCR) are the kinetic proofreading and serial ligation models. However, although kinetic proofreading provides a means for individual TCRs to measure accurately the length of time they are engaged and signal appropriately, the stochastic nature of ligand dissociation means the kinetic proofreading model implies that at high concentrations the response of the cell will be relatively nonspecific. Recent ligand experiments have revealed the phenomenon of both negative and positive crosstalk among neighboring TCRs. By using a Monte Carlo simulation of a lattice of TCRs, we integrate receptor crosstalk with the kinetic proofreading and serial ligation models and discover that receptor cooperativity can enhance T cell specificity significantly at a very modest cost to the sensitivity of the response.

Estimating the value of the metal-ligand bond dissociation enthalpy<D> (M-L) for adducts using empirical equations supported by TG data

In this work is presented and tested (for 106 adducts, mainly of the zinc group halides) two empirical equations supported in TG data to estimate the value of the metal-ligand bond dissociation enthalpy for adducts: <D> (M-O) = t i / g if t i < 420 K and <D> (M-O) = (t i / g ) - 7,75 . 10-2 . t i if t i > 420 K. In this empirical equations, t i is the thermodynamic temperature of the beginning of the thermal decomposition of the adduct, as determined by thermogravimetry, andg is a constant factor that is function of the metal halide considered and of the number of ligands, but is not dependant of the ligand itself. To half of the tested adducts the difference between experimental and calculated values was less than 5%. To about 80% of the tested adducts, the difference between the experimental (calorimetric) and the calculated (using the proposed equations) values are less than 15%.

Dissociation of thymic positive and negative selection in transgenic mice expressing major histocompatibility complex class I molecules exclusively on thymic cortical epithelial cells.

Thymic positive and negative selection of developing T lymphocytes confronts us with a paradox: How can a T-cell antigen receptor (TCR)-major histocompatibility complex (MHC)/peptide interaction in the former process lead to transduction of signals allowing for cell survival and in the latter induce programmed cell death or a hyporesponsive state known as anergy? One of the hypotheses put forward states that the outcome of a TCR-MHC/peptide interaction depends on the cell type presenting the selecting ligand to the developing thymocyte. Here we describe the development and lack of self-tolerance of CD8(+) T lymphocytes in transgenic mice expressing MHC class I molecules in the thymus exclusively on cortical epithelial cells. Despite the absence of MHC class I expression on professional antigen-presenting cells, normal numbers of CD8(+) cells were observed in the periphery. Upon specific activation, transgenic CD8(+) T cells efficiently lysed syngeneic MHC class I(+) targets in vitro and in vivo, indicating that thymic cortical epithelium (in contrast to medullary epithelium and antigen-presenting cells of hematopoietic origin) is incapable of tolerance induction. Thus, compartmentalization of the antigen-presenting cells involved in thymic positive selection and tolerance induction can (at least in part) explain the positive/negative selection paradox.

Mutations in the epithelial Na+ channel ENaC outer pore disrupt amiloride block by increasing its dissociation rate.

The epithelial Na+ channel ENaC mediates transepithelial Na+ transport in the distal kidney, the colon, and the lung and is a key element for the maintenance of Na+ balance and the regulation of blood pressure. Mutagenesis studies have identified residues alphaS583 and the homologous betaG525 and gammaG537 in the outer pore entrance that are critical for ENaC block by the K+-sparing diuretic amiloride. The aim of the present study was to determine first, whether these residues are part of the amiloride binding site, and second, whether they are general determinants of ENaC block by amiloride and its derivatives. Kinetic analysis of the association and dissociation rates of amiloride and benzamil to ENaC showed that mutation of residue alphaS583C and the homologous betaG525C increased the dissociation rate of the drugs from the binding site, with little changes in their association rate. Thus, these mutations destabilize the binding interaction between the blockers and the receptor on the channel, favoring the unbinding of the ligand. This strongly suggests that they are part of the binding site. Because mutations of alphaS583, betaG525, and gammaG537 have similar effects on amiloride, benzamil, and triamterene block, we conclude that these three ENaC blockers share a common receptor within the ion channel pore.

T cell receptor-ligand interactions: a conformational preequilibrium or an induced fit.

Kinetic parameters of T cell receptor (TCR) interactions with its ligand have been proposed to control T cell activation. Analysis of kinetic data obtained has so far produced conflicting insights; here, we offer a consideration of this problem. As a model system, association and dissociation of a soluble TCR (sT1) and its specific ligand, an azidobenzoic acid derivative of the peptide SYIPSAEK-(ABA)I (residues 252-260 from Plasmodium berghei circumsporozoite protein), bound to class I MHC H-2K(d)-encoded molecule (MHCp) were studied by surface plasmon resonance. The association time courses exhibited biphasic patterns. The fast and dominant phase was assigned to ligand association with the major fraction of TCR molecules, whereas the slow component was attributed to the presence of traces of TCR dimers. The association rate constant derived for the fast phase, assuming a reversible, single-step reaction mechanism, was relatively slow and markedly temperature-dependent, decreasing from 7.0 x 10(3) at 25 degrees C to 1.8 x 10(2) M(-1).s(-1) at 4 degrees C. Hence, it is suggested that these observed slow rate constants are the result of unresolved elementary steps of the process. Indeed, our analysis of the kinetic data shows that the time courses of TCR-MHCp interaction fit well to two different, yet closely related mechanisms, where an induced fit or a preequilibrium of two unbound TCR conformers are operational. These mechanisms may provide a rationale for the reported conformational flexibility of the TCR and its unusual ligand recognition properties, which combine high specificity with considerable crossreactivity.

Identification of Rare High-Avidity, Tumor-Reactive CD8+ T Cells by Monomeric TCR-Ligand Off-Rates Measurements on Living Cells.

The avidity of the T-cell receptor (TCR) for antigenic peptides presented by the peptide-MHC (pMHC) on cells is a key parameter for cell-mediated immunity. Yet a fundamental feature of most tumor antigen-specific CD8(+) T cells is that this avidity is low. In this study, we addressed the need to identify and select tumor-specific CD8(+) T cells of highest avidity, which are of the greatest interest for adoptive cell therapy in patients with cancer. To identify these rare cells, we developed a peptide-MHC multimer technology, which uses reversible Ni(2+)-nitrilotriacetic acid histidine tags (NTAmers). NTAmers are highly stable but upon imidazole addition, they decay rapidly to pMHC monomers, allowing flow-cytometric-based measurements of monomeric TCR-pMHC dissociation rates of living CD8(+) T cells on a wide avidity spectrum. We documented strong correlations between NTAmer kinetic results and those obtained by surface plasmon resonance. Using NTAmers that were deficient for CD8 binding to pMHC, we found that CD8 itself stabilized the TCR-pMHC complex, prolonging the dissociation half-life several fold. Notably, our NTAmer technology accurately predicted the function of large panels of tumor-specific T cells that were isolated prospectively from patients with cancer. Overall, our results demonstrated that NTAmers are effective tools to isolate rare high-avidity cytotoxic T cells from patients for use in adoptive therapies for cancer treatment.

The study of ligand binding specificities of the lipid binding proteins : recombinant human a-tocopherol transport protein (a-ttp), supernatant protein factor (spf) and S. cerevisiae Sec 14p for vitamin e (rrr-a-tocopherol) and other hydrophobic ligands.

One of the various functions of proteins in biological systems is the transport of small molecules, for this purpose proteins have naturally evolved special mechanisms to allow both ligand binding and its subsequent release to a target site; a process fundamental to many biological processes. Transport of Vitamin E (a-tocopherol), a lipid soluble antioxidant, to membranes helps in the protection of polyunsaturated fatty acids against peroxidative damage. In this research, the ligand binding characteristics of several members of the CRALTRIO family of lipid binding proteins was examined; the recombinant human a-Tocopherol Transfer Protein (a-TIP), Supernatant Protein Factor (SPF)ffocopherol Associated Protein (TAP), Cellular Retinaldehyde Binding Protein (CRALBP) and the phosphatidylinositol transfer protein from S. cerevisiae Sec 14p. Recombinant Sec 14p was expressed and purified from E. coli for comparison of tocopherol binding to the two other recombinant proteins postulated to traffic a-tocopherol. Competitive binding assays using [3H]-a-tocopherol and Lipidex-l000 resin allowed determination of the dissociation constants ~) of the CRAL-TRIO proteins for a-tocopherol and - 20 hydrophobic ligands for evaluation of the possible biological relevance of the binding interactions observed. The KIs (nM) for RRR-a-tocopherol are: a-TIP: 25.0, Sec 14p: 373, CRALBP: 528 and SPFffAP: 615. This indicates that all proteins recognize tocopherol but not with the same affinity. Sec 14p bound its native ligand PI with a KI of381 whereas SPFffAP bound PI (216) and y-tocopherol (268) similarly in contrast to the preferential binding ofRRR-a-tocopherol by a-TIP. Efforts to adequately represent biologically active SPFff AP involved investigation of tocopherol binding for several different recombinant proteins derived from different constructs and in the presence of different potential modulators (Ca+2, Mg+2, GTP and GDP); none of these conditions enhanced or inhibited a-tocopherol binding to SPF. This work suggests that only aTTP serves as the physiological mediator of a-tocopherol, yet structural homology between proteins allows common recognition of similar ligand features. In addition, several photo-affmity analogs of a-tocopherol were evaluated for their potential utility in further elucidation of a-TTP function or identification of novel tocopherol binding proteins.

Synthesis and characterization of BODIPY-α-tocopherol : a new ligand to explore the intracellular transfer of Vitamin E

Since its discovery nearly a century ago, a-tocopherol (vitamin E) research has been mainly focused on its ability to terminate the cycle of lipid peroxidation in membranes. Nitrobenzoxadiazole fluorescent analogues were made previously to study the intracellular transfer of vitamin E in cells. However, these molecules were reportedly susceptible to photobleaching while under illumination for transfer assays and microscopy. Here is reported the synthesis of a series of fluorescent analogues of vitamin E incorporating the more robust dipyrrometheneboron difluoride fluorophore (BDP-a-Tocs; Aex = 507 nm, Aem = 511 nm). C8-BDP-a-Toc 42c, having an eight-carbon chain between the chromanol and fluorophore, wa<; shown to bind specifically to a-tocopherol transfer protein with a dissociation constant of approximately 100 nM. Another fluorescent analogue of vitamin E with a thienyl derivative of BODIPY that is excited and fluoresces at longer wavelengths (Aex = 561 nm, Aem = 570 nm) is in development.

Oligomers of D-2 dopamine receptors - Evidence from ligand binding

There is increasing evidence that G protein-coupled receptors form oligomers and that this might be important for their function. We have studied this phenomenon for the D-2 dopamine receptor and have shown-using a variety of biochemical and biophysical techniques-that this receptor forms dimers or higher-order oligomers. Using ligand-binding studies, we have also found evidence that this oligomer formation has functional relevance. Thus, for the receptor expressed in either CHO cells or Sf 9 insect cells, the binding properties of several radioligands (in saturation, competition, and dissociation assays) do not conform to those expected for a monomeric receptor with a single binding site. We propose that the receptors exist in oligomers with homotropic and heterotropic negatively cooperative interactions between ligands.