Peroxisome proliferator-activated receptor structures: ligand specificity, molecular switch and interactions with regulators.

Peroxisome proliferator-activated receptors (PPARs) compose a family of nuclear receptors that mediate the effects of lipidic ligands at the transcriptional level. In this review, we highlight advances in the understanding of the PPAR ligand binding domain (LBD) structure at the atomic level. The overall structure of PPARs LBD is described, and important protein ligand interactions are presented. Structure-activity relationships between isotypes structures and ligand specificity are addressed. It is shown that the numerous experimental three-dimensional structures available, together with in silico simulations, help understanding the role played by the activating function-2 (AF-2) in PPARs activation and its underlying molecular mechanism. The relation between the PPARs constitutive activity and the intrinsic stability of the active conformation is discussed. Finally, the interactions of PPARs LBD with co-activators or co-repressors, as well as with the retinoid X receptor (RXR) are described and considered in relation to PPARs activation.

Structural analysis of TCR-ligand interactions studied on H-2Kd-restricted cloned CTL specific for a photoreactive peptide derivative.

To study the interaction of the TCR with its ligand, the complex of a MHC molecule and an antigenic peptide, we modified a TCR contact residue of a H-2Kd-restricted antigenic peptide with photoreactive 4-azidobenzoic acid. The photoreactive group was a critical component of the epitope recognized by CTL clones derived from mice immunized with such a peptide derivative. The majority of these clones expressed V beta 1-encoded beta chains that were paired with J alpha TA28-encoded alpha chains. For one of these TCR, the photoaffinity labeled sites were mapped on the alpha chain as a J alpha TA28-encoded tryptophan and on the beta chain as a residue of the C' strand of V beta 1. Molecular modeling of this TCR suggested the presence of a hydrophobic pocket that harbors this tryptophan as well as a tyrosine on the C' strand of V beta 1 between which the photoreactive side chain inserts. It is concluded that this avid binding principle may account for the preferential selection of V beta 1 and J alpha TA28-encoded TCR.

CD8 modulation of T-cell antigen receptor-ligand interactions on living cytotoxic T lymphocytes.

Thymocytes and class I major histocompatibility complex (MHC)-restricted cytotoxic T lymphocytes express predominantly heterodimeric alpha/beta CD8. By interacting with non-polymorphic regions of MHC class I molecules CD8 can mediate adhesion or by binding the same MHC molecules that interact with the T-cell antigen receptor (TCR) function as coreceptor in TCR-ligand binding and T-cell activation. Using TCR photoaffinity labelling with a soluble, monomeric photoreactive H-2Kd-peptide derivative complex, we report here that the avidity of TCR-ligand interactions on cloned cytotoxic T cells is very greatly strengthened by CD8. This is primarily explained by coordinate binding of ligand molecules by CD8 and TCR, because substitution of Asp 227 of Kd with Lys severely impaired the TCR-ligand binding on CD8+, but not CD8- cells. Kinetic studies on CD8+ and CD8- cells further showed that CD8 imposes distinct dynamics and a remarkable temperature dependence on TCR-ligand interactions. We propose that the ability of CD8 to act as coreceptor can be modulated by CD8-TCR interactions.

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.

Blind docking of 260 protein-ligand complexes with EADock 2.0.

Molecular docking softwares are one of the important tools of modern drug development pipelines. The promising achievements of the last 10 years emphasize the need for further improvement, as reflected by several recent publications (Leach et al., J Med Chem 2006, 49, 5851; Warren et al., J Med Chem 2006, 49, 5912). Our initial approach, EADock, showed a good performance in reproducing the experimental binding modes for a set of 37 different ligand-protein complexes (Grosdidier et al., Proteins 2007, 67, 1010). This article presents recent improvements regarding the scoring and sampling aspects over the initial implementation, as well as a new seeding procedure based on the detection of cavities, opening the door to blind docking with EADock. These enhancements were validated on 260 complexes taken from the high quality Ligand Protein Database [LPDB, (Roche et al., J Med Chem 2001, 44, 3592)]. Two issues were identified: first, the quality of the initial structures cannot be assumed and a manual inspection and/or a search in the literature are likely to be required to achieve the best performance. Second the description of interactions involving metal ions still has to be improved. Nonetheless, a remarkable success rate of 65% was achieved for a large scale blind docking assay, when considering only the top ranked binding mode and a success threshold of 2 A RMSD to the crystal structure. When looking at the five-top ranked binding modes, the success rate increases up to 76%. In a standard local docking assay, success rates of 75 and 83% were obtained, considering only the top ranked binding mode, or the five top binding modes, respectively.

Use of the FACTS solvation model for protein-ligand docking calculations. Application to EADock.

Protein-ligand docking has made important progress during the last decade and has become a powerful tool for drug development, opening the way to virtual high throughput screening and in silico structure-based ligand design. Despite the flattering picture that has been drawn, recent publications have shown that the docking problem is far from being solved, and that more developments are still needed to achieve high successful prediction rates and accuracy. Introducing an accurate description of the solvation effect upon binding is thought to be essential to achieve this goal. In particular, EADock uses the Generalized Born Molecular Volume 2 (GBMV2) solvent model, which has been shown to reproduce accurately the desolvation energies calculated by solving the Poisson equation. Here, the implementation of the Fast Analytical Continuum Treatment of Solvation (FACTS) as an implicit solvation model in small molecules docking calculations has been assessed using the EADock docking program. Our results strongly support the use of FACTS for docking. The success rates of EADock/FACTS and EADock/GBMV2 are similar, i.e. around 75% for local docking and 65% for blind docking. However, these results come at a much lower computational cost: FACTS is 10 times faster than GBMV2 in calculating the total electrostatic energy, and allows a speed up of EADock by a factor of 4. This study also supports the EADock development strategy relying on the CHARMM package for energy calculations, which enables straightforward implementation and testing of the latest developments in the field of Molecular Modeling.

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.

Electron microscopic visualization of protein-DNA interactions at the estrogen responsive element and in the first intron of the Xenopus laevis vitellogenin gene.

Stable protein-DNA complexes can be assembled in vitro at the 5' end of Xenopus laevis vitellogenin genes using extracts of nuclei from estrogen-induced frog liver and visualized by electron microscopy. Complexes at the three following sites can be identified on the gene B2: the transcription initiation site, the estrogen responsive element (ERE) and in the first intron. The complex at the transcription initiation site is stabilized by dinucleotides and thus represents a ternary transcription complex. The formation of the complexes at the two other sites is enhanced by estrogen and is reduced by tamoxifen, an antagonist of estrogen, while this latter effect is reversed by adding an excess of hormone. No sequence homology is apparent between the site containing the ERE and the binding site in intron I and functional tests in MCF-7 cells suggest that these two sites are not equivalent. Finally, we made use of previously characterized deletion mutants of the 5' flanking region of the gene B1, a close relative of the gene B2, to demonstrate that the 13-bp palindromic core element of the ERE is involved in the formation of the complexes observed upstream of the transcription initiation site.

Protein-protein interactions between adenovirus DNA polymerase and nuclear factor I mediate formation of the DNA replication preinitiation complex.

The in vitro adenovirus (Ad) DNA replication system provides an assay to study the interaction of viral and host replication proteins with the DNA template in the formation of the preinitiation complex. This initiation system requires in addition to the origin DNA sequences 1) Ad DNA polymerase (Pol), 2) Ad preterminal protein (pTP), the covalent acceptor for protein-primed DNA replication, and 3) nuclear factor I (NFI), a host cell protein identical to the CCAAT box-binding transcription factor. The interactions of these proteins were studied by coimmunoprecipitation and Ad origin DNA binding assays. The Ad Pol can bind to origin sequences only in the presence of another protein which can be either pTP or NFI. While NFI alone can bind to its origin recognition sequence, pTP does not specifically recognize DNA unless Ad Pol is present. Thus, protein-protein interactions are necessary for the targetting of either Ad Pol or pTP to the preinitiation complex. DNA footprinting demonstrated that the Ad DNA site recognized by the pTP.Pol complex was within the first 18 bases at the end of the template which constitutes the minimal origin of replication. Mutagenesis studies have defined the Ad Pol interaction site on NFI between amino acids 68-150, which overlaps the DNA binding and replication activation domain of this factor. A putative zinc finger on the Ad Pol has been mutated to a product that fails to bind the Ad origin sequences but still interacts with pTP. These results indicate that both protein-protein and protein-DNA interactions mediate specific recognition of the replication origin by Ad DNA polymerase.

Feedback on hypothalamic TRH transcription is dependent on thyroid hormone receptor N terminus.

The beta thyroid hormone receptor (TRbeta), but not TRalpha1, plays a specific role in mediating T(3)-dependent repression of hypothalamic TRH transcription. To investigate the structural basis of isoform specificity, we compared the transcriptional regulation and DNA binding obtained with chimeric and N-terminally deleted TRs. Using in vivo transfection assays to follow hypothalamic TRH transcription in the mouse brain, we found that TRbeta1 and chimeras with the TRbeta1 N terminus did not affect either transcriptional activation or repression from the rat TRH promoter, whereas N-terminally deleted TRbeta1 impaired T(3)-dependent repression. TRalpha1 or chimeras with the TRalpha1 N terminus reduced T(3)-independent transcriptional activation and blocked T(3)-dependent repression of transcription. Full deletion of the TRalpha1 N terminus restored ligand-independent activation of transcription. No TR isoform specificity was seen after transcription from a positive thyroid hormone response element. Gel mobility assays showed that all TRs tested bound specifically to the main negative thyroid hormone response element in the TRH promoter (site 4). Addition of neither steroid receptor coactivator 1 nor nuclear extracts from the hypothalamic paraventricular nuclei revealed any TR isoform specificity in binding to site 4. Thus N-terminal sequences specify TR T(3)-dependent repression of TRH transcription but not DNA recognition, emphasizing as yet unknown neuron-specific contributions to protein-promoter interactions in vivo.

NO rebinding to myoglobin: a reactive molecular dynamics study.

The rebinding of NO to myoglobin after photolysis is studied using the 'reactive molecular dynamics' method. In this approach the energy of the system is evaluated on two potential energy surfaces that include the heme-ligand interactions which change between liganded and unliganded myoglobin. This makes it possible to take into account in a simple way, the high dimensionality of the transition seam connecting the reactant and product states. The dynamics of the dissociated NO molecules are examined, and the geometrical and energetic properties of the transition seam are studied. Analysis of the frequency of recrossing shows that the height of the effective rebinding barrier is dependent on the time after photodissociation. This effect is due mainly to protein relaxation and may contribute to the experimentally observed non-exponential rebinding rate of NO, as has been suggested previously.

Genome-wide association study identifies two loci strongly affecting transferrin glycosylation.

Polysaccharide sidechains attached to proteins play important roles in cell-cell and receptor-ligand interactions. Variation in the carbohydrate component has been extensively studied for the iron transport protein transferrin, because serum levels of the transferrin isoforms asialotransferrin + disialotransferrin (carbohydrate-deficient transferrin, CDT) are used as biomarkers of excessive alcohol intake. We conducted a genome-wide association study to assess whether genetic factors affect CDT concentration in serum. CDT was measured in three population-based studies: one in Switzerland (CoLaus study, n = 5181) and two in Australia (n = 1509, n = 775). The first cohort was used as the discovery panel and the latter ones served as replication. Genome-wide single-nucleotide polymorphism (SNP) typing data were used to identify loci with significant associations with CDT as a percentage of total transferrin (CDT%). The top three SNPs in the discovery panel (rs2749097 near PGM1 on chromosome 1, and missense polymorphisms rs1049296, rs1799899 in TF on chromosome 3) were successfully replicated , yielding genome-wide significant combined association with CDT% (P = 1.9 × 10(-9), 4 × 10(-39), 5.5 × 10(-43), respectively) and explain 5.8% of the variation in CDT%. These allelic effects are postulated to be caused by variation in availability of glucose-1-phosphate as a precursor of the glycan (PGM1), and variation in transferrin (TF) structure.

Crystal packing modifies ligand binding affinity: the case of aldose reductase.

The relationship between the structures of protein-ligand complexes existing in the crystal and in solution, essential in the case of fragment-based screening by X-ray crystallography (FBS-X), has been often an object of controversy. To address this question, simultaneous co-crystallization and soaking of two inhibitors with different ratios, Fidarestat (FID; K(d) = 6.5 nM) and IDD594 (594; K(d) = 61 nM), which bind to h-aldose reductase (AR), have been performed. The subatomic resolution of the crystal structures allows the differentiation of both inhibitors, even when the structures are almost superposed. We have determined the occupation ratio in solution by mass spectrometry (MS) Occ(FID)/Occ(594) = 2.7 and by X-ray crystallography Occ(FID)/Occ(594) = 0.6. The occupancies in the crystal and in solution differ 4.6 times, implying that ligand binding potency is influenced by crystal contacts. A structural analysis shows that the Loop A (residues 122-130), which is exposed to the solvent, is flexible in solution, and is involved in packing contacts within the crystal. Furthermore, inhibitor 594 contacts the base of Loop A, stabilizing it, while inhibitor FID does not. This is shown by the difference in B-factors of the Loop A between the AR-594 and AR-FID complexes. A stable loop diminishes the entropic energy barrier to binding, favoring 594 versus FID. Therefore, the effect of the crystal environment should be taken into consideration in the X-ray diffraction analysis of ligand binding to proteins. This conclusion highlights the need for additional methodologies in the case of FBS-X to validate this powerful screening technique, which is widely used.

Fatty acids, eicosanoids, and hypolipidemic agents identified as ligands of peroxisome proliferator-activated receptors by coactivator-dependent receptor ligand assay.

Peroxisome proliferator-activated receptors (PPARs) are nuclear hormone receptors controlling the expression of genes involved in lipid homeostasis. PPARs activate gene transcription in response to a variety of compounds including hypolipidemic drugs as well as natural fatty acids. From the plethora of PPAR activators, Scatchard analysis of receptor-ligand interactions has thus far identified only four ligands. These are the chemotactic agent leukotriene B4 and the hypolipidemic drug Wy 14,643 for the alpha-subtype and a prostaglandin J2 metabolite and synthetic antidiabetic thiazolidinediones for the gamma-subtype. Based on the hypothesis that ligand binding to PPAR would induce interactions of the receptor with transcriptional coactivators, we have developed a novel ligand sensor assay, termed coactivator-dependent receptor ligand assay (CARLA). With CARLA we have screened several natural and synthetic candidate ligands and have identified naturally occurring fatty acids and metabolites as well as hypolipidemic drugs as bona fide ligands of the three PPAR subtypes from Xenopus laevis. Our results suggest that PPARs, by their ability to interact with a number of structurally diverse compounds, have acquired unique ligand-binding properties among the superfamily of nuclear receptors that are compatible with their biological activity.