EADock: docking of small molecules into protein active sites with a multiobjective evolutionary optimization.

In recent years, protein-ligand docking has become a powerful tool for drug development. Although several approaches suitable for high throughput screening are available, there is a need for methods able to identify binding modes with high accuracy. This accuracy is essential to reliably compute the binding free energy of the ligand. Such methods are needed when the binding mode of lead compounds is not determined experimentally but is needed for structure-based lead optimization. We present here a new docking software, called EADock, that aims at this goal. It uses an hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein-ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 A around the center of mass of the ligand position in the crystal structure, and on the contrary to other benchmarks, our algorithm was fed with optimized ligand positions up to 10 A root mean square deviation (RMSD) from the crystal structure, excluding the latter. This validation illustrates the efficiency of our sampling strategy, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 A, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures could be explained by the presence of crystal contacts in the experimental structure. Finally, the ability of EADock to accurately predict binding modes on a real application was illustrated by the successful docking of the RGD cyclic pentapeptide on the alphaVbeta3 integrin, starting far away from the binding pocket.

Soluble major histocompatibility complex-peptide octamers with impaired CD8 binding selectively induce Fas-dependent apoptosis.

Fluorescence-labeled soluble major histocompatibility complex class I-peptide "tetramers" constitute a powerful tool to detect and isolate antigen-specific CD8(+) T cells by flow cytometry. Conventional "tetramers" are prepared by refolding of heavy and light chains with a specific peptide, enzymatic biotinylation at an added C-terminal biotinylation sequence, and "tetramerization" by reaction with phycoerythrin- or allophycocyanin-labeled avidin derivatives. We show here that such preparations are heterogeneous and describe a new procedure that allows the preparation of homogeneous tetra- or octameric major histocompatibility complex-peptide complexes. These compounds were tested on T1 cytotoxic T lymphocytes (CTLs), which recognize the Plasmodium berghei circumsporzoite peptide 252-260 (SYIPSAEKI) containing photoreactive 4-azidobenzoic acid on Lys(259) in the context of H-2K(d). We report that mutation of the CD8 binding site of K(d) greatly impairs the binding of tetrameric but not octameric or multimeric K(d)-PbCS(ABA) complexes to CTLs. This mutation abolishes the ability of the octamer to elicit significant phosphorylation of CD3, intracellular calcium mobilization, and CTL degranulation. Remarkably, however, this octamer efficiently activates CTLs for Fas (CD95)-dependent apoptosis.

Relationship between estrogen receptor alpha location and gene induction reveals the importance of downstream sites and cofactors.

BACKGROUND: To understand cancer-related modifications to transcriptional programs requires detailed knowledge about the activation of signal-transduction pathways and gene expression programs. To investigate the mechanisms of target gene regulation by human estrogen receptor alpha (hERalpha), we combine extensive location and expression datasets with genomic sequence analysis. In particular, we study the influence of patterns of DNA occupancy by hERalpha on expression phenotypes. RESULTS: We find that strong ChIP-chip sites co-localize with strong hERalpha consensus sites and detect nucleotide bias near hERalpha sites. The localization of ChIP-chip sites relative to annotated genes shows that weak sites are enriched near transcription start sites, while stronger sites show no positional bias. Assessing the relationship between binding configurations and expression phenotypes, we find binding sites downstream of the transcription start site (TSS) to be equally good or better predictors of hERalpha-mediated expression as upstream sites. The study of FOX and SP1 cofactor sites near hERalpha ChIP sites shows that induced genes frequently have FOX or SP1 sites. Finally we integrate these multiple datasets to define a high confidence set of primary hERalpha target genes. CONCLUSION: Our results support the model of long-range interactions of hERalpha with the promoter-bound cofactor SP1 residing at the promoter of hERalpha target genes. FOX motifs co-occur with hERalpha motifs along responsive genes. Importantly we show that the spatial arrangement of sites near the start sites and within the full transcript is important in determining response to estrogen signaling.

Characterization of HbpR binding by site-directed mutagenesis of its DNA-binding site and by deletion of the effector domain.

In the presence of 2-hydroxybiphenyl, the enhancer binding protein, HbpR, activates the sigma54-dependent P(hbpC) promoter and controls the initial steps of 2-hydroxybiphenyl degradation in Pseudomonas azelaica. In the activation process, an oligomeric HbpR complex of unknown subunit composition binds to an operator region containing two imperfect palindromic sequences. Here, the HbpR-DNA binding interactions were investigated by site-directed mutagenesis of the operator region and by DNA-binding assays using purified HbpR. Mutations that disrupted the twofold symmetry in the palindromes did not affect the binding affinity of HbpR, but various mutations along a 60 bp region, and also outside the direct palindromic sequences, decreased the binding affinity. Footprints of HbpR on mutant operator fragments showed that a partial loss of binding contacts occurs, suggesting that the binding of one HbpR 'protomer' in the oligomeric complex is impaired whilst leaving the other contacts intact. An HbpR variant, devoid of its N-terminal sensing A-domain, was unable to activate transcription from the hbpC promoter while maintaining protection of the operator DNA in footprints. Wild-type HbpR was unable to activate transcription from the hbpC promoter when delta A-HbpR was expressed in the same cell, suggesting the formation of (repressing) hetero-oligomers. This model implies that HbpR can self-associate on its operator DNA without effector recognition or ATP binding. Furthermore, our findings suggest that the N-terminal sensing domain of HbpR is needed to activate the central ATPase domain rather than to repress a constitutively active C domain, as is the case for the related regulatory protein XylR.

Nuclear translocation and DNA recognition signals colocalized within the bZIP domain of cyclic adenosine 3',5'-monophosphate response element-binding protein CREB.

CREB is a cAMP-responsive nuclear DNA-binding protein that binds to cAMP response elements and stimulates gene transcription upon activation of the cAMP signalling pathway. The protein consists of an amino-terminal transcriptional transactivation domain and a carboxyl-terminal DNA-binding domain (bZIP domain) comprised of a basic region and a leucine zipper involved in DNA recognition and dimerization, respectively. Recently, we discovered a testis-specific transcript of CREB that contains an alternatively spliced exon encoding multiple stop codons. CREB encoded by this transcript is a truncated protein lacking the bZIP domain. We postulated that the antigen detected by CREB antiserum in the cytoplasm of germinal cells is the truncated CREB that must also lack its nuclear translocation signal (NTS). To test this hypothesis we prepared multiple expression plasmids encoding carboxyl-terminal deletions of CREB and transiently expressed them in COS-1 cells. By Western immunoblot analysis as well as immunocytochemistry of transfected cells, we show that CREB proteins truncated to amino acid 286 or shorter are sequestered in the cytoplasm, whereas a CREB of 295 amino acids is translocated into the nucleus. Chimeric CREBs containing a heterologous NTS fused to the first 248 or 261 amino acids of CREB are able to drive the translocation of the protein into the nucleus. Thus, the nine amino acids in the basic region involved in DNA recognition between positions 287 and 295 (RRKKKEYVK) of CREB contain the NTS. Further, mutation of the lysine at position 290 in CREB to an asparagine diminishes nuclear translocation of the protein.(ABSTRACT TRUNCATED AT 250 WORDS)

Ligand binding to heparan sulfate proteoglycans induces their aggregation and distribution along actin cytoskeleton

Cell surface heparan sulfate proteoglycans (HSPGs) participate in molecular events that regulate cell adhesion, migration, and proliferation. The present study demonstrates that soluble heparin-binding proteins or cross-linking antibodies induce the aggregation of cell surface HSPGs and their distribution along underlying actin filaments. Immunofluorescence and confocal microscopy and immunogold and electron microscopy indicate that, in the absence of ligands, HSPGs are irregularly distributed on the fibroblast cell surface, without any apparent codistribution with the actin cytoskeleton. In the presence of ligand (lipoprotein lipase) or antibodies against heparan sulfate, HSPGs aggregate and colocalize with the actin cytoskeleton. Triton X-100 extraction and immunoelectron microscopy have demonstrated that in this condition HSPGs were clustered and associated with the actin filaments. Crosslinking experiments that use biotinylated lipoprotein lipase have revealed three major proteoglycans as binding sites at the fibroblast cell surface. These cross-linked proteoglycans appeared in the Triton X-100 insoluble fraction. Platinum/carbon replicas of the fibroblast surface incubated either with lipoprotein lipase or antiheparan sulfate showed large aggregates of HSPGs regularly distributed along cytoplasmic fibers. Quantification of the spacing between HSPGs by confocal microscopy confirmed that the nonrandom distribution of HSPG aggregates along the actin cytoskeleton was induced by ligand binding. When cells were incubated either with lipoprotein lipase or antibodies against heparan sulfate, the distance between immunofluorescence spots was uniform. In contrast, the spacing between HSPGs on fixed cells not incubated with ligand was more variable. This highly organized spatial relationship between actin and proteoglycans suggests that cortical actin filaments could organize the molecular machinery involved in signal transduction and molecular movements on the cell surface that are triggered by heparin-binding proteins.

CNDOL: A fast and reliable method for the calculation of electronic properties of very large systems. Applications to retinal binding pocket in rhodopsin and gas phase porphine.

Very large molecular systems can be calculated with the so called CNDOL approximate Hamiltonians that have been developed by avoiding oversimplifications and only using a priori parameters and formulas from the simpler NDO methods. A new diagonal monoelectronic term named CNDOL/21 shows great consistency and easier SCF convergence when used together with an appropriate function for charge repulsion energies that is derived from traditional formulas. It is possible to obtain a priori molecular orbitals and electron excitation properties after the configuration interaction of single excited determinants with reliability, maintaining interpretative possibilities even being a simplified Hamiltonian. Tests with some unequivocal gas phase maxima of simple molecules (benzene, furfural, acetaldehyde, hexyl alcohol, methyl amine, 2,5 dimethyl 2,4 hexadiene, and ethyl sulfide) ratify the general quality of this approach in comparison with other methods. The calculation of large systems as porphine in gas phase and a model of the complete retinal binding pocket in rhodopsin with 622 basis functions on 280 atoms at the quantum mechanical level show reliability leading to a resulting first allowed transition in 483 nm, very similar to the known experimental value of 500 nm of "dark state." In this very important case, our model gives a central role in this excitation to a charge transfer from the neighboring Glu(-) counterion to the retinaldehyde polyene chain. Tests with gas phase maxima of some important molecules corroborate the reliability of CNDOL/2 Hamiltonians.

Building complexity in O2-binding copper complexes : site-selective metalation and intermolecular O2-binding at dicopper and heterometallic complexes derived from an unsymmetric ligand

A novel unsymmetric dinucleating ligand (LN3N4) combining a tridentate and a tetradentate binding sites linked through a m-xylyl spacer was synthesized as ligand scaffold for preparing homo- and dimetallic complexes, where the two metal ions are bound in two different coordination environments. Site-selective binding of different metal ions is demonstrated. LN3N4 is able to discriminate between CuI and a complementary metal (M′ = CuI, ZnII, FeII, CuII, or GaIII) so that pure heterodimetallic complexes with a general formula [CuIM′(LN3N4)]n+ are synthesized. Reaction of the dicopper(I) complex [CuI 2(LN3N4)]2+ with O2 leads to the formation of two different copper-dioxygen (Cu2O2) intermolecular species (O and TP) between two copper atoms located in the same site from different complex molecules. Taking advantage of this feature, reaction of the heterodimetallic complexes [CuM′(LN3N4)]n+ with O2 at low temperature is used as a tool to determine the final position of the CuI center in the system because only one of the two Cu2O2 species is formed

Selective destruction of nigrostriatal dopaminergic neurons does not alter [3H]-ryanodine binding in rat striatum

Dopamine nigrostriatal neurons are important for motor control and may contain a particularly dense population of ryanodine receptors involved in the control of dopamine release. To test this hypothesis, we used a classical model of unilateral selective lesion of these neurons in rats based on 6-hydroxydopamine (6-OHDA) injection into the substantia nigra. Binding of [3H]-GBR 12935, used as a presynaptic marker since it labels specifically the dopamine uptake complex, was dramatically decreased by 83-100% in striatum homogenates after 6-OHDA lesion. On the contrary, no reduction of [3H]-ryanodine binding was observed. The present data indicate that [3H]-ryanodine binding sites present in rat striatum are not preferentially localized in dopaminergic terminals.

Biological activity and binding of estradiol to SK-Mel 23 human melanoma cells

Patients expressing estradiol receptors in melanoma cells have been reported to have a better prognosis. We therefore decided to investigate the in vitro effects of ß-estradiol and tamoxifen on the growth and tyrosinase activity of SK-Mel 23 human melanoma cells. Twenty-four-hour treatment with 0.4 nM ß-estradiol inhibited cell proliferation in 30% (0.70 ± 0.03 x 10(5) cells) and increased tyrosinase activity in 50% (7130.5 ± 376.5 cpm/10(5) cells), as compared to untreated cells (1.0 ± 0.05 x 10(5) cells and 4769 ± 25.5 cpm/10(5) cells, respectively). Both responses were completely (100%) blocked by 1 µM tamoxifen. Higher concentrations (up to 1.6 nM) or longer treatments (up to 72 h) did not result in a larger effect of the hormone on proliferation or tyrosinase activity. Competition binding assays demonstrated the presence of binding sites to [2,4,6,7-³H]-ß-estradiol, and that the tritiated analogue was displaced by the unlabeled hormone (1 nM to 100 µM, Kd = 0.14 µM, maximal displacement of 93%) or by 10 µM tamoxifen (displacement of 60%). ß-estradiol also increased the phosphorylated state of two proteins of 16 and 46 kDa, after 4-h treatment, as determined by Western blot. The absorbance of each band was 1.9- and 4-fold the controls, respectively, as determined with Image-Pro Plus software. Shorter incubation periods with ß-estradiol did not enhance phosporylation; after 6-h treatment with the hormone, the two proteins returned to the control phosphorylation levels. The growth inhibition promoted by estradiol may explain the better prognosis of melanoma-bearing women as compared to men, and open new perspectives for drug therapy.

Thermodynamic evaluation and modeling of proton and water exchange associated with benzamidine and berenil binding to ß-trypsin

Serine-proteases are involved in vital processes in virtually all species. They are important targets for researchers studying the relationships between protein structure and activity, for the rational design of new pharmaceuticals. Trypsin was used as a model to assess a possible differential contribution of hydration water to the binding of two synthetic inhibitors. Thermodynamic parameters for the association of bovine ß-trypsin (homogeneous material, observed 23,294.4 ± 0.2 Da, theoretical 23,292.5 Da) with the inhibitors benzamidine and berenil at pH 8.0, 25ºC and with 25 mM CaCl2, were determined using isothermal titration calorimetry and the osmotic stress method. The association constant for berenil was about 12 times higher compared to the one for benzamidine (binding constants are K = 596,599 ± 25,057 and 49,513 ± 2,732 M-1, respectively; the number of binding sites is the same for both ligands, N = 0.99 ± 0.05). Apparently the driving force responsible for this large difference of affinity is not due to hydrophobic interactions because the variation in heat capacity (DCp), a characteristic signature of these interactions, was similar in both systems tested (-464.7 ± 23.9 and -477.1 ± 86.8 J K-1 mol-1 for berenil and benzamidine, respectively). The results also indicated that the enzyme has a net gain of about 21 water molecules regardless of the inhibitor tested. It was shown that the difference in affinity could be due to a larger number of interactions between berenil and the enzyme based on computational modeling. The data support the view that pharmaceuticals derived from benzamidine that enable hydrogen bond formation outside the catalytic binding pocket of ß-trypsin may result in more effective inhibitors.

Effects on prolactin secretion and binding to dopaminergic receptors in sleep-deprived lupus-prone mice

Sleep disturbances have far-reaching effects on the neuroendocrine and immune systems and may be linked to disease manifestation. Sleep deprivation can accelerate the onset of lupus in NZB/NZWF1 mice, an animal model of severe systemic lupus erythematosus. High prolactin (PRL) concentrations are involved in the pathogenesis of systemic lupus erythematosus in human beings, as well as in NZB/NZWF1 mice. We hypothesized that PRL could be involved in the earlier onset of the disease in sleep-deprived NZB/NZWF1 mice. We also investigated its binding to dopaminergic receptors, since PRL secretion is mainly controlled by dopamine. Female NZB/NZWF1 mice aged 9 weeks were deprived of sleep using the multiple platform method. Blood samples were taken for the determination of PRL concentrations and quantitative receptor autoradiography was used to map binding of the tritiated dopaminergic receptor ligands [³H]-SCH23390, [³H]-raclopride and [³H]-WIN35,428 to D1 and D2 dopaminergic receptors and dopamine transporter sites throughout the brain, respectively. Sleep deprivation induced a significant decrease in plasma PRL secretion (2.58 ± 0.95 ng/mL) compared with the control group (25.25 ± 9.18 ng/mL). The binding to D1 and D2 binding sites was not significantly affected by sleep deprivation; however, dopamine transporter binding was significantly increased in subdivisions of the caudate-putamen - posterior (16.52 ± 0.5 vs 14.44 ± 0.6), dorsolateral (18.84 ± 0.7 vs 15.97 ± 0.7) and ventrolateral (24.99 ± 0.5 vs 22.54 ± 0.7 µCi/g), in the sleep-deprived mice when compared to the control group. These results suggest that PRL is not the main mechanism involved in the earlier onset of the disease observed in sleep-deprived NZB/NZWF1 mice and the reduction of PRL concentrations after sleep deprivation may be mediated by modifications in the dopamine transporter sites of the caudate-putamen.

PPARγ induces growth inhibition and apoptosis through upregulation of insulin-like growth factor-binding protein-3 in gastric cancer cells

Peroxisome proliferator activator receptor-gamma (PPARγ) is a ligand-activated transcriptional factor involved in the carcinogenesis of various cancers. Insulin-like growth factor-binding protein-3 (IGFBP-3) is a tumor suppressor gene that has anti-apoptotic activity. The purpose of this study was to investigate the anticancer mechanism of PPARγ with respect to IGFBP-3. PPARγ was overexpressed in SNU-668 gastric cancer cells using an adenovirus gene transfer system. The cells in which PPARγ was overexpressed exhibited growth inhibition, induction of apoptosis, and a significant increase in IGFBP-3 expression. We investigated the underlying molecular mechanisms of PPARγ in SNU-668 cells using an IGFBP-3 promoter/luciferase reporter system. Luciferase activity was increased up to 15-fold in PPARγ transfected cells, suggesting that PPARγ may directly interact with IGFBP-3 promoter to induce its expression. Deletion analysis of the IGFBP-3 promoter showed that luciferase activity was markedly reduced in cells without putative p53-binding sites (-Δ1755, -Δ1795). This suggests that the critical PPARγ-response region is located within the p53-binding region of the IGFBP-3 promoter. We further demonstrated an increase in PPARγ-induced luciferase activity even in cells treated with siRNA to silence p53 expression. Taken together, these data suggest that PPARγ exhibits its anticancer effect by increasing IGFBP-3 expression, and that IGFBP-3 is a significant tumor suppressor.

Loss of noradrenaline transporter sites in frontal cortex of rats with acute (ischemic) liver failure.

There is increasing evidence that central noradrenaline (NA) transport mechanisms are implicated in the central nervous system complications of acute liver failure. In order to assess this possibility, binding sites for the high affinity NA transporter ligand [3H]-nisoxetine were measured by quantitative receptor autoradiography in the brains of rats with acute liver failure resulting from hepatic devascularization and in appropriate controls. In vivo microdialysis was used to measure extracellular brain concentrations of NA. Severe encephalopathy resulted in a significant loss of [3H]-nisoxetine sites in frontal cortex and a concomitant increase in extracellular brain concentrations of NA in rats with acute liver failure. A loss of transporter sites was also observed in thalamus of rats with acute liver failure. This loss of NA transporter sites could result from depletion of central NA stores due to a reserpine-like effect of ammonia which is known to accumulate to millimolar concentrations in brain in ischemic liver failure. Impaired NA transport and the consequent increase in synaptic concentrations and increased stimulation of neuronal and astrocytic noradrenergic receptors could be implicated in the pathogenesis of the encephalopathy and brain edema characteristic of acute liver failure.

Une nouvelle approche computationnelle pour la découverte des sites de fixation de facteurs de transcription à l’ADN, adaptée aux données de ChIP-chip et de ChIP-séquençage

Les facteurs de transcription sont des protéines spécialisées qui jouent un rôle important dans différents processus biologiques tel que la différenciation, le cycle cellulaire et la tumorigenèse. Ils régulent la transcription des gènes en se fixant sur des séquences d’ADN spécifiques (éléments cis-régulateurs). L’identification de ces éléments est une étape cruciale dans la compréhension des réseaux de régulation des gènes. Avec l’avènement des technologies de séquençage à haut débit, l’identification de tout les éléments fonctionnels dans les génomes, incluant gènes et éléments cis-régulateurs a connu une avancée considérable. Alors qu’on est arrivé à estimer le nombre de gènes chez différentes espèces, l’information sur les éléments qui contrôlent et orchestrent la régulation de ces gènes est encore mal définie. Grace aux techniques de ChIP-chip et de ChIP-séquençage il est possible d’identifier toutes les régions du génome qui sont liées par un facteur de transcription d’intérêt. Plusieurs approches computationnelles ont été développées pour prédire les sites fixés par les facteurs de transcription. Ces approches sont classées en deux catégories principales: les algorithmes énumératifs et probabilistes. Toutefois, plusieurs études ont montré que ces approches génèrent des taux élevés de faux négatifs et de faux positifs ce qui rend difficile l’interprétation des résultats et par conséquent leur validation expérimentale. Dans cette thèse, nous avons ciblé deux objectifs. Le premier objectif a été de développer une nouvelle approche pour la découverte des sites de fixation des facteurs de transcription à l’ADN (SAMD-ChIP) adaptée aux données de ChIP-chip et de ChIP-séquençage. Notre approche implémente un algorithme hybride qui combine les deux stratégies énumérative et probabiliste, afin d’exploiter les performances de chacune d’entre elles. Notre approche a montré ses performances, comparée aux outils de découvertes de motifs existants sur des jeux de données simulées et des jeux de données de ChIP-chip et de ChIP-séquençage. SAMD-ChIP présente aussi l’avantage d’exploiter les propriétés de distributions des sites liés par les facteurs de transcription autour du centre des régions liées afin de limiter la prédiction aux motifs qui sont enrichis dans une fenêtre de longueur fixe autour du centre de ces régions. Les facteurs de transcription agissent rarement seuls. Ils forment souvent des complexes pour interagir avec l’ADN pour réguler leurs gènes cibles. Ces interactions impliquent des facteurs de transcription dont les sites de fixation à l’ADN sont localisés proches les uns des autres ou bien médier par des boucles de chromatine. Notre deuxième objectif a été d’exploiter la proximité spatiale des sites liés par les facteurs de transcription dans les régions de ChIP-chip et de ChIP-séquençage pour développer une approche pour la prédiction des motifs composites (motifs composés par deux sites et séparés par un espacement de taille fixe). Nous avons testé ce module pour prédire la co-localisation entre les deux demi-sites ERE qui forment le site ERE, lié par le récepteur des œstrogènes ERα. Ce module a été incorporé à notre outil de découverte de motifs SAMD-ChIP.