Odorant binding proteins of the red imported fire ant, Solenopsis invicta: an example of the problems facing the analysis of widely divergent proteins.

We describe the odorant binding proteins (OBPs) of the red imported fire ant, Solenopsis invicta, obtained from analyses of an EST library and separate 454 sequencing runs of two normalized cDNA libraries. We identified a total of 18 putative functional OBPs in this ant. A third of the fire ant OBPs are orthologs to honey bee OBPs. Another third of the OBPs belong to a lineage-specific expansion, which is a common feature of insect OBP evolution. Like other OBPs, the different fire ant OBPs share little sequence similarity (∼ 20%), rendering evolutionary analyses difficult. We discuss the resulting problems with sequence alignment, phylogenetic analysis, and tests of selection. As previously suggested, our results underscore the importance for careful exploration of the sensitivity to the effects of alignment methods for data comprising widely divergent sequences.

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)

Mutations in the TGFβ binding-protein-like domain 5 of FBN1 are responsible for acromicric and geleophysic dysplasias.

Geleophysic (GD) and acromicric dysplasia (AD) belong to the acromelic dysplasia group and are both characterized by severe short stature, short extremities, and stiff joints. Although AD has an unknown molecular basis, we have previously identified ADAMTSL2 mutations in a subset of GD patients. After exome sequencing in GD and AD cases, we selected fibrillin 1 (FBN1) as a candidate gene, even though mutations in this gene have been described in Marfan syndrome, which is characterized by tall stature and arachnodactyly. We identified 16 heterozygous FBN1 mutations that are all located in exons 41 and 42 and encode TGFβ-binding protein-like domain 5 (TB5) of FBN1 in 29 GD and AD cases. Microfibrillar network disorganization and enhanced TGFβ signaling were consistent features in GD and AD fibroblasts. Importantly, a direct interaction between ADAMTSL2 and FBN1 was demonstrated, suggesting a disruption of this interaction as the underlying mechanism of GD and AD phenotypes. Although enhanced TGFβ signaling caused by FBN1 mutations can trigger either Marfan syndrome or GD and AD, our findings support the fact that TB5 mutations in FBN1 are responsible for short stature phenotypes.

Disease-associated mutations in the actin-binding domain of filamin B cause cytoplasmic focal accumulations correlating with disease severity.

Dominant missense mutations in FLNB, encoding the actin-cross linking protein filamin B (FLNB), cause a broad range of skeletal dysplasias with varying severity by an unknown mechanism. Here these FLNB mutations are shown to cluster in exons encoding the actin-binding domain (ABD) and filamin repeats surrounding the flexible hinge 1 region of the FLNB rod domain. Despite being positioned in domains that bind actin, it is unknown if these mutations perturb cytoskeletal structure. Expression of several full-length FLNB constructs containing ABD mutations resulted in the appearance of actin-containing cytoplasmic focal accumulations of the substituted protein to a degree that was correlated with the severity of the associated phenotypes. In contrast, study of mutations leading to substitutions in the FLNB rod domain that result in the same phenotypes as ABD mutations demonstrated that with only one exception disease-associated substitutions, surrounding hinge 1 demonstrated no tendency to form actin-filamin foci. The exception, a substitution in filamin repeat 6, lies within a region previously implicated in filamin-actin binding. These data are consistent with mutations in the ABD conferring enhanced actin-binding activity but suggest that substitutions affecting repeats near the flexible hinge region of FLNB precipitate the same phenotypes through a different mechanism.

Inferring the mode of origin of polyploid species from next-generation sequence data.

Many eukaryote organisms are polyploid. However, despite their importance, evolutionary inference of polyploid origins and modes of inheritance has been limited by a need for analyses of allele segregation at multiple loci using crosses. The increasing availability of sequence data for nonmodel species now allows the application of established approaches for the analysis of genomic data in polyploids. Here, we ask whether approximate Bayesian computation (ABC), applied to realistic traditional and next-generation sequence data, allows correct inference of the evolutionary and demographic history of polyploids. Using simulations, we evaluate the robustness of evolutionary inference by ABC for tetraploid species as a function of the number of individuals and loci sampled, and the presence or absence of an outgroup. We find that ABC adequately retrieves the recent evolutionary history of polyploid species on the basis of both old and new sequencing technologies. The application of ABC to sequence data from diploid and polyploid species of the plant genus Capsella confirms its utility. Our analysis strongly supports an allopolyploid origin of C. bursa-pastoris about 80 000 years ago. This conclusion runs contrary to previous findings based on the same data set but using an alternative approach and is in agreement with recent findings based on whole-genome sequencing. Our results indicate that ABC is a promising and powerful method for revealing the evolution of polyploid species, without the need to attribute alleles to a homeologous chromosome pair. The approach can readily be extended to more complex scenarios involving higher ploidy levels.

T cell recognition of hapten. Anatomy of T cell receptor binding of a H-2kd-associated photoreactive peptide derivative.

To elucidate the structural basis of T cell recognition of hapten-modified antigenic peptides, we studied the interaction of the T1 T cell antigen receptor (TCR) with its ligand, the H-2Kd-bound Plasmodium berghei circumsporozoite peptide 252-260 (SYIPSAEKI) containing photoreactive 4-azidobenzoic acid (ABA) on P. berghei circumsporozoite Lys259. The photoaffinity-labeled TCR residue(s) were mapped as Tyr48 and/or Tyr50 of complementary determining region 2beta (CDR2beta). Other TCR-ligand contacts were identified by mutational analysis. Molecular modeling, based on crystallographic coordinates of closely related TCR and major histocompatibility complex I molecules, indicated that ABA binds strongly and specifically in a cavity between CDR3alpha and CDR2beta. We conclude that TCR expressing selective Vbeta and CDR3alpha sequences form a binding domain between CDR3alpha and CDR2beta that can accommodate nonpeptidic moieties conjugated at the C-terminal portion of peptides binding to major histocompatibility complex (MHC) encoded proteins.

Selective cooperation between fatty acid binding proteins and peroxisome proliferator-activated receptors in regulating transcription.

Lipophilic compounds such as retinoic acid and long-chain fatty acids regulate gene transcription by activating nuclear receptors such as retinoic acid receptors (RARs) and peroxisome proliferator-activated receptors (PPARs). These compounds also bind in cells to members of the family of intracellular lipid binding proteins, which includes cellular retinoic acid-binding proteins (CRABPs) and fatty acid binding proteins (FABPs). We previously reported that CRABP-II enhances the transcriptional activity of RAR by directly targeting retinoic acid to the receptor. Here, potential functional cooperation between FABPs and PPARs in regulating the transcriptional activities of their common ligands was investigated. We show that adipocyte FABP and keratinocyte FABP (A-FABP and K-FABP, respectively) selectively enhance the activities of PPARgamma and PPARbeta, respectively, and that these FABPs massively relocate to the nucleus in response to selective ligands for the PPAR isotype which they activate. We show further that A-FABP and K-FABP interact directly with PPARgamma and PPARbeta and that they do so in a receptor- and ligand-selective manner. Finally, the data demonstrate that the presence of high levels of K-FABP in keratinocytes is essential for PPARbeta-mediated induction of differentiation of these cells. Taken together, the data establish that A-FABP and K-FABP govern the transcriptional activities of their ligands by targeting them to cognate PPARs in the nucleus, thereby enabling PPARs to exert their biological functions.

The cAMP response element binding protein-2 (CREB-2) can interact with the C/EBP-homologous protein (CHOP).

cAMP response element binding protein-2 (CREB-2) is a basic leucine zipper (bZIP) factor that was originally described as a repressor of CRE-dependent transcription but that can also act as a transcriptional activator. Moreover, CREB-2 is able to function in association with the viral Tax protein as an activator of the human T-cell leukemia virus type I (HTLV-I) promoter. Here we show that CREB-2 is able to interact with C/EBP-homologous protein (CHOP), a bZIP transcription factor known to inhibit CAAT/enhancer-dependent transcription. Cotransfection of CHOP with CREB-2 results in decreased activation driven by the cellular CRE motif or the HTLV-I proximal Tax-responsive element, confirming that CREB-2 and CHOP can interact with each other in vivo.

Molecular characterization of a human matrix attachment region epigenetic regulator.

Matrix attachment regions (MAR) generally act as epigenetic regulatory sequences that increase gene expression, and they were proposed to partition chromosomes into loop-forming domains. However, their molecular mode of action remains poorly understood. Here, we assessed the possible contribution of the AT-rich core and adjacent transcription factor binding motifs to the transcription augmenting and anti-silencing effects of human MAR 1-68. Either flanking sequences together with the AT-rich core were required to obtain the full MAR effects. Shortened MAR derivatives retaining full MAR activity were constructed from combinations of the AT-rich sequence and multimerized transcription factor binding motifs, implying that both transcription factors and the AT-rich microsatellite sequence are required to mediate the MAR effect. Genomic analysis indicated that MAR AT-rich cores may be depleted of histones and enriched in RNA polymerase II, providing a molecular interpretation of their chromatin domain insulator and transcriptional augmentation activities.

Prospective study of insulin-like growth factor-I, insulin-like growth factor-binding protein 3, genetic variants in the IGF1 and IGFBP3 genes and risk of coronary artery disease.

Although experimental studies have suggested that insulin-like growth factor I (IGF-I) and its binding protein IGFBP-3 might have a role in the aetiology of coronary artery disease (CAD), the relevance of circulating IGFs and their binding proteins in the development of CAD in human populations is unclear. We conducted a nested case-control study, with a mean follow-up of six years, within the EPIC-Norfolk cohort to assess the association between circulating levels of IGF-I and IGFBP-3 and risk of CAD in up to 1,013 cases and 2,055 controls matched for age, sex and study enrolment date. After adjustment for cardiovascular risk factors, we found no association between circulating levels of IGF-I or IGFBP-3 and risk of CAD (odds ratio: 0.98 (95% Cl 0.90-1.06) per 1 SD increase in circulating IGF-I; odds ratio: 1.02 (95% Cl 0.94-1.12) for IGFBP-3). We examined associations between tagging single nucleotide polymorphisms (tSNPs) at the IGF1 and IGFBP3 loci and circulating IGF-I and IGFBP-3 levels in up to 1,133 cases and 2,223 controls and identified three tSNPs (rs1520220, rs3730204, rs2132571) that showed independent association with either circulating IGF-I or IGFBP-3 levels. In an assessment of 31 SNPs spanning the IGF1 or IGFBP3 loci, none were associated with risk of CAD in a meta-analysis that included EPIC-Norfolk and eight additional studies comprising up to 9,319 cases and 19,964 controls. Our results indicate that IGF-I and IGFBP-3 are unlikely to be importantly involved in the aetiology of CAD in human populations.

Quantitative modelling of transcription factor binding sites

Abstract One of the most important issues in molecular biology is to understand regulatory mechanisms that control gene expression. Gene expression is often regulated by proteins, called transcription factors which bind to short (5 to 20 base pairs),degenerate segments of DNA. Experimental efforts towards understanding the sequence specificity of transcription factors is laborious and expensive, but can be substantially accelerated with the use of computational predictions. This thesis describes the use of algorithms and resources for transcriptionfactor binding site analysis in addressing quantitative modelling, where probabilitic models are built to represent binding properties of a transcription factor and can be used to find new functional binding sites in genomes. Initially, an open-access database(HTPSELEX) was created, holding high quality binding sequences for two eukaryotic families of transcription factors namely CTF/NF1 and LEFT/TCF. The binding sequences were elucidated using a recently described experimental procedure called HTP-SELEX, that allows generation of large number (> 1000) of binding sites using mass sequencing technology. For each HTP-SELEX experiments we also provide accurate primary experimental information about the protein material used, details of the wet lab protocol, an archive of sequencing trace files, and assembled clone sequences of binding sequences. The database also offers reasonably large SELEX libraries obtained with conventional low-throughput protocols.The database is available at http://wwwisrec.isb-sib.ch/htpselex/ and and ftp://ftp.isrec.isb-sib.ch/pub/databases/htpselex. The Expectation-Maximisation(EM) algorithm is one the frequently used methods to estimate probabilistic models to represent the sequence specificity of transcription factors. We present computer simulations in order to estimate the precision of EM estimated models as a function of data set parameters(like length of initial sequences, number of initial sequences, percentage of nonbinding sequences). We observed a remarkable robustness of the EM algorithm with regard to length of training sequences and the degree of contamination. The HTPSELEX database and the benchmarked results of the EM algorithm formed part of the foundation for the subsequent project, where a statistical framework called hidden Markov model has been developed to represent sequence specificity of the transcription factors CTF/NF1 and LEF1/TCF using the HTP-SELEX experiment data. The hidden Markov model framework is capable of both predicting and classifying CTF/NF1 and LEF1/TCF binding sites. A covariance analysis of the binding sites revealed non-independent base preferences at different nucleotide positions, providing insight into the binding mechanism. We next tested the LEF1/TCF model by computing binding scores for a set of LEF1/TCF binding sequences for which relative affinities were determined experimentally using non-linear regression. The predicted and experimentally determined binding affinities were in good correlation.

Tumor sensitizing function and mechanism of action of a RasGAP derived peptide

Cancer is the second cause of death after cardio-vascular diseases in economically developed countries. Two of the most commonly used anti-cancer therapies are chemo and radiotherapy. Despite the remarkable advances made in term of delivery and specificity of these two anti-tumor regimens, their toxicity towards healthy tissue remains a limitation. A promising approach to overcome this obstacle would be the utilization of therapeutic peptides that specifically augment the sensitivity of tumoral cells to treatments. Lower therapeutical doses would then be required to kill malignant cells, limiting toxic effects on healthy tissues. It was previously shown in our laboratory that the caspase-3 generated fragment N2 of RasGAP is able to potentiate the genotoxin-induced apoptosis selectively in cancer cells. In this work we show that fragment N2 strictly requires a cytoplasmic localization to deliver its pro-apoptotic effect in genotoxin-treated cancer cells. The tumor sensitizing capacity of fragment N2 was found to reside within the 10 amino acid sequence 317-326. Our laboratory earlier demonstrated that a peptide corresponding to amino acids 317 to 326 of RasGAP fused to the TAT cell permeable moiety, called TAT-RasGAP317.326, is able to sensitize cancer cells, but not normal cells, to genotoxin-induced apoptosis. In the present study we describe the capacity of TAT-RasGAP 317.326 to sensitize tumors to both chemo and radiotherapy in an in vivo mouse model. The molecular mechanism underlying the TAT-RasGAP 317.326-mediated sensitization starts now to be elucidated. We demonstrate that G3BP1, an endoribonuclease binding to amino acids 317-326 of RasGAP, is not involved in the sensitization mechanism. We also provide evidence showing that TAT-RasGAP3 17-326 potentiates the genotoxin-mediated activation of Bax in a tBid-dependent manner. Altogether our results show that TAT-RasGAP 317.326 could be potentially used in cancer therapy as sensitizer, in order to improve the efficacy of chemo and radiotherapy and prolong the life expectancy of cancer patients. Moreover, the understanding of the TAT-RasGAP317.326 mode of action might help to unravel the mechanisms by which cancer cells resist to chemo and radiotherapy and therefore to design more targeted and efficient anti-tumoral strategies.

Glial hyaluronate-binding protein expression in aggregating brain cell cultures.

We analyzed the expression of glial hyaluronate-binding protein (GHAP), an integral component of the extracellular matrix, in aggregating brain cell cultures of fetal rat telencephalon using immunofluorescence. GHAP immunoreactivity appeared after 1 week in culture, simultaneous with the first deposits of myelin basic protein, and showed a development-dependent increase. Comparison of glia-enriched and neuron-enriched cultures showed that only glial cells express GHAP. Three peptide growth factors, epidermal growth factor, fibroblast growth factor and platelet-derived growth factor, which are known to stimulate the differentiation of glial cells, modulated the deposit of GHAP immunoreactivity. The 3-dimensional structure of aggregate cultures promoted GHAP deposition, suggesting that cell-cell interactions are required for extracellular matrix formation. Furthermore GHAP production seemed to depend on the developmental stage of the glial cells.

Inhibition of the shade avoidance response by formation of non-DNA binding bHLH heterodimers.

In shade-intolerant plants such as Arabidopsis, a reduction in the red/far-red (R/FR) ratio, indicative of competition from other plants, triggers a suite of responses known as the shade avoidance syndrome (SAS). The phytochrome photoreceptors measure the R/FR ratio and control the SAS. The phytochrome-interacting factors 4 and 5 (PIF4 and PIF5) are stabilized in the shade and are required for a full SAS, whereas the related bHLH factor HFR1 (long hypocotyl in FR light) is transcriptionally induced by shade and inhibits this response. Here we show that HFR1 interacts with PIF4 and PIF5 and limits their capacity to induce the expression of shade marker genes and to promote elongation growth. HFR1 directly inhibits these PIFs by forming non-DNA-binding heterodimers with PIF4 and PIF5. Our data indicate that PIF4 and PIF5 promote SAS by directly binding to G-boxes present in the promoter of shade marker genes, but their action is limited later in the shade when HFR1 accumulates and forms non-DNA-binding heterodimers. This negative feedback loop is important to limit the response of plants to shade.

Ligands for pheromone-sensing neurons are not conformationally activated odorant binding proteins.

Pheromones form an essential chemical language of intraspecific communication in many animals. How olfactory systems recognize pheromonal signals with both sensitivity and specificity is not well understood. An important in vivo paradigm for this process is the detection mechanism of the sex pheromone (Z)-11-octadecenyl acetate (cis-vaccenyl acetate [cVA]) in Drosophila melanogaster. cVA-evoked neuronal activation requires a secreted odorant binding protein, LUSH, the CD36-related transmembrane protein SNMP, and the odorant receptor OR67d. Crystallographic analysis has revealed that cVA-bound LUSH is conformationally distinct from apo (unliganded) LUSH. Recombinantly expressed mutant versions of LUSH predicted to enhance or diminish these structural changes produce corresponding alterations in spontaneous and/or cVA-evoked activity when infused into olfactory sensilla, leading to a model in which the ligand for pheromone receptors is not free cVA, but LUSH that is "conformationally activated" upon cVA binding. Here we present evidence that contradicts this model. First, we demonstrate that the same LUSH mutants expressed transgenically affect neither basal nor pheromone-evoked activity. Second, we compare the structures of apo LUSH, cVA/LUSH, and complexes of LUSH with non-pheromonal ligands and find no conformational property of cVA/LUSH that can explain its proposed unique activated state. Finally, we show that high concentrations of cVA can induce neuronal activity in the absence of LUSH, but not SNMP or OR67d. Our findings are not consistent with the model that the cVA/LUSH complex acts as the pheromone ligand, and suggest that pheromone molecules alone directly activate neuronal receptors.