Beta(1)-selective agonist (-)-1-(3,4-dimethoxyphenetylamino)-3-(3,4-dihydroxy)-2-propanol [(-)-RO363] differentially interacts with key amino acids responsible for beta(1)-selective binding in resting and active states.

(-)-1-(3,4-Dimethoxyphenetylamino)-3-(3,4-dihydroxy)-2-propanol [(-)-RO363] is a highly selective beta(1)-adrenergic receptor (beta(1)AR) agonist. To study the binding site of beta(1)-selective agonist, chimeric beta(1)/beta(2)ARs and Ala-substituted beta(1)ARs were constructed. Several key residues of beta(1)AR [Leu(110) and Thr(117) in transmembrane domain (TMD) 2], and Phe(359) in TMD 7] were found to be responsible for beta(1)-selective binding of (-)-RO363, as determined by competitive binding. Based on these results, we built a three-dimensional model of the binding domain for (-)-RO363. The model indicated that TMD 2 and TMD 7 of beta(1)AR form a binding pocket; the methoxyphenyl group of N-substituent of (-)-RO363 seems to locate within the cavity surrounded by Leu(110), Thr(117), and Phe(359). The amino acids Leu(110) and Phe(359) interact with the phenyl ring of (-)-RO363, whereas Thr(117) forms hydrogen bond with the methoxy group of (-)-RO363. To examine the interaction of these residues with beta(1)AR in an active state, each of the amino acids was changed to Ala in a constitutively active (CA)-beta(1)AR mutant. The degree of decrease in the affinity of CA-beta(1)AR for (-)-RO363 was essentially the same as that of wild-type beta(1)AR when mutated at Leu(110) and Thr(117). However, the affinity was decreased in Ala-substituted mutant of Phe(359) compared with that of wild-type beta(1)AR. These results indicated that Leu(110) and Thr(117) are necessary for the initial binding of (-)-RO363 with beta(1)-selectivity, and interaction of Phe(359) with the N-substituent of (-)-RO363 in an active state is stronger than in the resting state.

Decreased binding of peptides-MHC class I (pMHC) multimeric complexes to CD8 affects their binding avidity for the TCR but does not significantly impact on pMHC/TCR dissociation rate.

The CD8 coreceptor plays a crucial role in both T cell development in the thymus and in the activation of mature T cells in response to Ag-specific stimulation. In this study we used soluble peptides-MHC class I (pMHC) multimeric complexes bearing mutations in the CD8 binding site that impair their binding to the MHC, together with altered peptide ligands, to assess the impact of CD8 on pMHC binding to the TCR. Our data support a model in which CD8 promotes the binding of TCR to pMHC. However, once the pMHC/TCR complex is formed, the TCR dominates the pMHC/TCR dissociation rates. As a consequence of these molecular interactions, under physiologic conditions CD8 plays a key role in complex formation, resulting in the enhancement of CD8 T cell functions whose specificity, however, is determined by the TCR.

Ezrin directly interacts with the alpha1b-adrenergic receptor and plays a role in receptor recycling.

Using the yeast two-hybrid system, we identified ezrin as a protein interacting with the C-tail of the alpha1b-adrenergic receptor (AR). The interaction was shown to occur in vitro between the receptor C-tail and the N-terminal portion of ezrin, or Four-point-one ERM (FERM) domain. The alpha1b-AR/ezrin interaction occurred inside the cells as shown by the finding that the transfected alpha1b-AR and FERM domain or ezrin could be coimmunoprecipitated from human embryonic kidney 293 cell extracts. Mutational analysis of the alpha1b-AR revealed that the binding site for ezrin involves a stretch of at least four arginines on the receptor C-tail. The results from both receptor biotinylation and immunofluorescence experiments indicated that the FERM domain impaired alpha1b-AR recycling to the plasma membrane without affecting receptor internalization. The dominant negative effect of the FERM domain, which relies on its ability to mask the ezrin binding site for actin, was mimicked by treatment of cells with cytochalasin D, an actin depolymerizing agent. A receptor mutant (DeltaR8) lacking its binding site in the C-tail for ezrin displayed delayed receptor recycling. These findings identify ezrin as a new protein directly interacting with a G protein-coupled receptor and demonstrate the direct implication of ezrin in GPCR trafficking via an actin-dependent mechanism.

ATP release due to Thy-1-integrin binding induces P2X7-mediated calcium entry required for focal adhesion formation.

Thy-1, an abundant mammalian glycoprotein, interacts with αvβ3 integrin and syndecan-4 in astrocytes and thus triggers signaling events that involve RhoA and its effector p160ROCK, thereby increasing astrocyte adhesion to the extracellular matrix. The signaling cascade includes calcium-dependent activation of protein kinase Cα upstream of Rho; however, what causes the intracellular calcium transients required to promote adhesion remains unclear. Purinergic P2X7 receptors are important for astrocyte function and form large non-selective cation pores upon binding to their ligand, ATP. Thus, we evaluated whether the intracellular calcium required for Thy-1-induced cell adhesion stems from influx mediated by ATP-activated P2X7 receptors. Results show that adhesion induced by the fusion protein Thy-1-Fc was preceded by both ATP release and sustained intracellular calcium elevation. Elimination of extracellular ATP with Apyrase, chelation of extracellular calcium with EGTA, or inhibition of P2X7 with oxidized ATP, all individually blocked intracellular calcium increase and Thy-1-stimulated adhesion. Moreover, Thy-1 mutated in the integrin-binding site did not trigger ATP release, and silencing of P2X7 with specific siRNA blocked Thy-1-induced adhesion. This study is the first to demonstrate a functional link between αvβ3 integrin and P2X7 receptors, and to reveal an important, hitherto unanticipated, role for P2X7 in calcium-dependent signaling required for Thy-1-stimulated astrocyte adhesion.

Manganese-induced integrin affinity maturation promotes recruitment of alpha V beta 3 integrin to focal adhesions in endothelial cells: evidence for a role of phosphatidylinositol 3-kinase and Src.

Integrin activity is controlled by changes in affinity (i.e. ligand binding) and avidity (i.e. receptor clustering). Little is known, however, about the effect of affinity maturation on integrin avidity and on the associated signaling pathways. To study the effect of affinity maturation on integrin avidity, we stimulated human umbilical vein endothelial cells (HUVEC) with MnCl(2) to increase integrin affinity and monitored clustering of beta 1 and beta 3 integrins. In unstimulated HUVEC, beta 1 integrins were present in fibrillar adhesions, while alpha V beta 3 was detected in peripheral focal adhesions. Clustered beta 1 and beta 3 integrins expressed high affinity/ligand-induced binding site (LIBS) epitopes. MnCl(2)-stimulation promoted focal adhesion and actin stress fiber formation at the basal surface of the cells, and strongly enhanced mAb LM609 staining and expression of beta 3 high affinity/LIBS epitopes at focal adhesions. MnCl(2)-induced alpha V beta 3 clustering was blocked by a soluble RGD peptide, by wortmannin and LY294002, two pharmacological inhibitors of phosphatidylinositol 3-kinase (PI 3-K), and by over-expressing a dominant negative PI 3-K mutant protein. Conversely, over-expression of active PI 3-K and pharmacological inhibiton of Src with PP2 and CGP77675, enhanced basal and manganese-induced alpha V beta 3 clustering. Transient increased phosphorylation of protein kinase B/Akt, a direct target of PI 3K, occurred upon manganese stimulation. MnCl(2) did not alter beta 1 integrin distribution or beta1 high-affinity/LIBS epitope expression. Based on these results, we conclude that MnCl(2)-induced alpha V beta 3 integrin affinity maturation stimulates focal adhesion and actin stress fiber formation, and promotes recruitment of high affinity alpha V beta 3 to focal adhesions. Affinity-modulated alpha V beta 3 clustering requires PI3-K signaling and is negatively regulate by Src.

Posttranscriptional repression of GacS/GacA-controlled genes by the RNA-binding protein RsmE acting together with RsmA in the biocontrol strain Pseudomonas fluorescens CHA0.

In the plant-beneficial soil bacterium Pseudomonas fluorescens CHA0, the production of biocontrol factors (antifungal secondary metabolites and exoenzymes) is controlled at a posttranscriptional level by the GacS/GacA signal transduction pathway involving RNA-binding protein RsmA as a key regulatory element. This protein is assumed to bind to the ribosome-binding site of target mRNAs and to block their translation. RsmA-mediated repression is relieved at the end of exponential growth by two GacS/GacA-controlled regulatory RNAs RsmY and RsmZ, which bind and sequester the RsmA protein. A gene (rsmE) encoding a 64-amino-acid RsmA homolog was identified and characterized in strain CHA0. Overexpression of rsmE strongly reduced the expression of target genes (hcnA, for a hydrogen cyanide synthase subunit; aprA, for the main exoprotease; and phlA, for a component of 2,4-diacetylphloroglucinol biosynthesis). Single null mutations in either rsmA or rsmE resulted in a slight increase in the expression of hcnA, aprA, and phlA. By contrast, an rsmA rsmE double mutation led to strongly increased and advanced expression of these target genes and completely suppressed a gacS mutation. Both the RsmE and RsmA levels increased with increasing cell population densities in strain CHA0; however, the amount of RsmA showed less variability during growth. Expression of rsmE was controlled positively by GacA and negatively by RsmA and RsmE. Mobility shift assays demonstrated specific binding of RsmE to RsmY and RsmZ RNAs. The transcription and stability of both regulatory RNAs were strongly reduced in the rsmA rsmE double mutant. In conclusion, RsmA and RsmE together account for maximal repression in the GacS/GacA cascade of strain CHA0.

Artificial zinc finger fusions targeting Sp1-binding sites and the trans-activator-responsive element potently repress transcription and replication of HIV-1.

Tat activates transcription by interacting with Sp1, NF-kappaB, positive transcription elongation factor b, and trans-activator-responsive element (TAR). Tat and Sp1 play major roles in transcription by protein-protein interactions at human immunodeficiency virus, type 1 (HIV-1) long terminal repeat. Sp1 activates transcription by interacting with cyclin T1 in the absence of Tat. To disrupt the transcription activation by Tat and Sp1, we fused Sp1-inhibiting polypeptides, zinc finger polypeptide, and the TAR-binding mutant Tat (TatdMt) together. A designed or natural zinc finger and Tat mutant fusion was used to target the fusion to the key regulatory sites (GC box and TAR) on the long terminal repeat and nascent short transcripts to disrupt the molecular interaction that normally result in robust transcription. The designed zinc finger and TatdMt fusions were targeted to the TAR, and they potently repressed both transcription and replication of HIV-1. The Sp1-inhibiting POZ domain, TatdMt, and zinc fingers are key functional domains important in repression of transcription and replication. The designed artificial zinc fingers were targeted to the high affinity Sp1-binding site, and by being fused with TatdMt and POZ domain, they strongly block both Sp1-cyclin T1-dependent transcription and Tat-dependent transcription, even in the presence of excess expressed Tat.

Homeodomain binding motifs modulate the probability of odorant receptor gene choice in transgenic mice.

Odorant receptor (OR) genes constitute with 1200 members the largest gene family in the mouse genome. A mature olfactory sensory neuron (OSN) is thought to express just one OR gene, and from one allele. The cell bodies of OSNs that express a given OR gene display a mosaic pattern within a particular region of the main olfactory epithelium. The mechanisms and cis-acting DNA elements that regulate the expression of one OR gene per OSN - OR gene choice - remain poorly understood. Here, we describe a reporter assay to identify minimal promoters for OR genes in transgenic mice, which are produced by the conventional method of pronuclear injection of DNA. The promoter transgenes are devoid of an OR coding sequence, and instead drive expression of the axonal marker tau-β-galactosidase. For four mouse OR genes (M71, M72, MOR23, and P3) and one human OR gene (hM72), a mosaic, OSN-specific pattern of reporter expression can be obtained in transgenic mice with contiguous DNA segments of only ~300 bp that are centered around the transcription start site (TSS). The ~150bp region upstream of the TSS contains three conserved sequence motifs, including homeodomain (HD) binding sites. Such HD binding sites are also present in the H and P elements, DNA sequences that are known to strongly influence OR gene expression. When a 19mer encompassing a HD binding site from the P element is multimerized nine times and added upstream of a MOR23 minigene that contains the MOR23 coding region, we observe a dramatic increase in the number of transgene-expressing founders and lines and in the number of labeled OSNs. By contrast, a nine times multimerized 19mer with a mutant HD binding site does not have these effects. We hypothesize that HD binding sites in the H and P elements and in OR promoters modulate the probability of OR gene choice.

High-throughput SELEX SAGE method for quantitative modeling of transcription-factor binding sites.

The ability to determine the location and relative strength of all transcription-factor binding sites in a genome is important both for a comprehensive understanding of gene regulation and for effective promoter engineering in biotechnological applications. Here we present a bioinformatically driven experimental method to accurately define the DNA-binding sequence specificity of transcription factors. A generalized profile was used as a predictive quantitative model for binding sites, and its parameters were estimated from in vitro-selected ligands using standard hidden Markov model training algorithms. Computer simulations showed that several thousand low- to medium-affinity sequences are required to generate a profile of desired accuracy. To produce data on this scale, we applied high-throughput genomics methods to the biochemical problem addressed here. A method combining systematic evolution of ligands by exponential enrichment (SELEX) and serial analysis of gene expression (SAGE) protocols was coupled to an automated quality-controlled sequence extraction procedure based on Phred quality scores. This allowed the sequencing of a database of more than 10,000 potential DNA ligands for the CTF/NFI transcription factor. The resulting binding-site model defines the sequence specificity of this protein with a high degree of accuracy not achieved earlier and thereby makes it possible to identify previously unknown regulatory sequences in genomic DNA. A covariance analysis of the selected sites revealed non-independent base preferences at different nucleotide positions, providing insight into the binding mechanism.

Synergistic action of GA-binding protein and glucocorticoid receptor in transcription from the mouse mammary tumor virus promoter.

B lymphocytes are among the first cells to be infected by mouse mammary tumor virus (MMTV), and they play a crucial role in its life cycle. To study transcriptional regulation of MMTV in B cells, we have analyzed two areas of the long terminal repeat (LTR) next to the glucocorticoid receptor binding site, fp1 (at position -139 to -146 from the cap site) and fp2 (at -157 to -164). Both showed B-cell-specific protection in DNase I in vitro footprinting assays and contain binding sites for Ets transcription factors, a large family of proteins involved in cell proliferation and differentiation and oncogenic transformation. In gel retardation assays, fp1 and fp2 bound the heterodimeric Ets factor GA-binding protein (GABP) present in B-cell nuclear extracts, which was identified by various criteria: formation of dimers and tetramers, sensitivity to pro-oxidant conditions, inhibition of binding by specific antisera, and comigration of complexes with those formed by recombinant GABP. Mutations which prevented complex formation in vitro abolished glucocorticoid-stimulated transcription from an MMTV LTR linked to a reporter gene in transiently transfected B-cell lines, whereas they did not affect the basal level. Exogenously expressed GABP resulted in an increased level of hormone response of the LTR reporter plasmid and produced a synergistic effect with the coexpressed glucocorticoid receptor, indicating cooperation between the two. This is the first example of GABP cooperation with a steroid receptor, providing the opportunity for studying the integration of their intracellular signaling pathways.

Organisation and role of actin microfilaments during cellular growth and morphogenesis in the moss Physcomitrella patens

ABSTRACT The network of actin cytoskeleton is composed of actin filaments (F-actin) that are made by polymerisation of actin monomers and actin binding proteins. It is required for growth and morphogenesis of eukaryotic cells. The labelling of F-actin with constitutively expressed GFP-Talin (Kost et al., 1998) reveals the organisation of cellular actin networks in plants. Due to the lack of information on actin cytoskeleton through gametophytic development of the model moss plant Physcornitrella patens, stable transgenic lines overexpressing GFP-Talin were generated to detect F-actin structures. It is shown that the 35S promoter driven expression is not suitable for F-actin labelling in all cells. When it is replaced by the inducible heat-shock promoter Gmhsp17.3 from soybean, one hour mild heat stress at 37°C followed by recovery at 25°C is enough to induce efficient and transient labelling in all tissues without altering cellular morphology. The optimal observations of F-actin structures at different stages of moss development can be done between 12-18 hours after the induction. By using confocal microscopy, we demonstrate that stellated actin arrays were densely accumulated at the growing tip in regenerating protoplasts, apical protonemal cells and rhizoids and connected with a fine dispersed F-actin mesh. Following three-dimensional growth, the cortical star-like structures are widespread in the meristematic cells of developing bud and young gametophores. On the contrary, undulating networks of actin cables are found at the final stage of cell differentiation. During redifferentiation of mature leaf cells into protonemal filaments the rather stagnant web of actin cables is replaced by diffuse actin meshwork. In eukaryotes, nucleation of the actin monomers prior to their polymerization is driven by the seven-subunit ARP2/3 complex and formins. We cloned the gene encoding the ARP3 subunit of P. patens and generated arp3 mutants of the moss through gene disruption. The knockout of ARP3 affects the elongation of chloronemal cells and blocks further differentiation of caulonemal cells and rhizoids, and the gametophores are slightly stunted compared to wild-type. The arp mutants were created in the heat-shock inducible GFP-Talin strains allowing us to visualise a disorganised actin network and a lack of star-like actin cytoskeleton arrays. We conclude that ARP2/3 dependent nucleation of actin filaments is critical for the growth of filamentous cells, which in turn influences moss colonization. In complementation assays, the overexpression of Physcomitrella and Arab idopsis ARP3 genes in the moss arp3 mutant results in full recovery of wild type phenotype. In contrast the ARP3 subunit of fission yeast is not able to complement the moss arp3 mutant of moss indicating that regulation of the ARP2/3 dependent actin nucleation diverged in different kingdoms. RESUME Le réseau d'actine est composé de filaments de F-actine et d'un ensemble de protéines s'y attachant (Actin binding proteins). Le réseau d'actine est nécessaire à la croissance et à la morphogenèse de toutes les cellules eucaryotes. Chez les plantes, le marquage ainsi que l'étude de l'organisation du réseau d'actine ont été réalisés en utilisant une fusion GFP-Talin (Kost et al., 1998) exprimée sous le control d'un promoteur constitutif. Afin d'étudier les structures F-actine dans les cellules de Physcomitrella Patens et pour combler le manque d'information sur le développement des gamétophores, des lignées transgéniques stables surexprimant GFP-Talin ont été crées. Nous avons démontré que l'utilisation du promoteur 35S est inadéquate pour le marquage complet et homogène des filaments d'actine dans toutes les cellules de P. patens. Par contre, l'utilisation du promoteur inductible Gmhsp17.3 nous a permis de réaliser un marquage transitoire et général dans tous les tissus de la mousse. Une heure de choc thermique à 37°C suivis d'un temps de récupération de 12-18h à 25°C sont les conditions optimales (sans dommages cellulaires) pour l'observation des structures F-actine à différentes étapes de développement de la mousse. En utilisant la microscopie confocale, nous avons observé l'existence de structures F-actine accumulées en forme d'étoiles. Ces structures, qui sont liées au réseau de microfilaments d'actine, ont été observées dans les protoplastes en régénération, les cellules des protonema apicales ainsi que dans les rhizoïdes. En suivant la croissance tridimensionnelle, ces structures en étoiles ont été observées dans les cellules meristématiques des bourgeons et des jeunes gamétophores. Par contre, dans les cellules différentiées ces structures laissent place à des réseaux de câbles épais. Nous avons également remarqué que durant la redifferentiation des cellules foliaires le réseau de câbles de F-actine est remplacé par un réseau de F-actine diffus. Dans les cellules eucaryotes, la nucléation des filaments d'actirie précédant leur polymérisation est contrôlé par sept sous unités du complexe ARP2/3 et par des formines. Nous avons isolé le gène codant pour la sous unité ARP3 de P. patens et nous avons crée des mutants arp3 par intégration ciblée (Knockout). L'élongation des cellules chloronema est clairement affectée dans les mutants arp3. La différentiation des caulonemata et des rhizoïdes est bloquée et les gametophores sont légèrement plus courts comparé au type sauvage. A fin d'étudier l'organisation des filaments d'actines dans les mutants arp3, nous avons aussi réalisé un arp3-knockout dans la lignée Hsp-GFP-Talin. La nouvelle lignée générée nous a permis de visualiser une désorganisation du réseau d'actine et une absence complète de structures de F-actine accumulée en forme d'étoiles. Les résultats obtenus nous amènent à conclure que la nucléation (ARP2/3 dépendante) des filaments d'actine est indispensable à la croissance des cellules filamenteuses. Par conséquent, les filaments d'actine semblent avoir un rôle dans la colonisation des milieux par les mousses. Nous avons également procédé à des essais de complémentation du mutant arp3. La surexpression des gènes ARP3 de Physcomitrella et d'Arabidopsis dans les cellules du mutant arp3 rétabli complètement le phénotype WT. Par contre, le gène ARP3 des levures n'est pas suffisant pour complémenter la même mutation dans les cellules de mousses. Ce résultat démontre que les mécanismes de régulation de la nucléation des filaments d'actine (ARP2/3 dépendante) sont différents entre les différents groupes d'eucaryotes.

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.

Experimental analysis and computer prediction of CTF/NFI transcription factor DNA binding sites.

Accurate prediction of transcription factor binding sites is needed to unravel the function and regulation of genes discovered in genome sequencing projects. To evaluate current computer prediction tools, we have begun a systematic study of the sequence-specific DNA-binding of a transcription factor belonging to the CTF/NFI family. Using a systematic collection of rationally designed oligonucleotides combined with an in vitro DNA binding assay, we found that the sequence specificity of this protein cannot be represented by a simple consensus sequence or weight matrix. For instance, CTF/NFI uses a flexible DNA binding mode that allows for variations of the binding site length. From the experimental data, we derived a novel prediction method using a generalised profile as a binding site predictor. Experimental evaluation of the generalised profile indicated that it accurately predicts the binding affinity of the transcription factor to natural or synthetic DNA sequences. Furthermore, the in vitro measured binding affinities of a subset of oligonucleotides were found to correlate with their transcriptional activities in transfected cells. The combined computational-experimental approach exemplified in this work thus resulted in an accurate prediction method for CTF/NFI binding sites potentially functioning as regulatory regions in vivo.

Transcriptional repression of peroxisome proliferator-activated receptor beta/delta in murine keratinocytes by CCAAT/enhancer-binding proteins.

The roles of peroxisome proliferator-activated receptors (PPARs) and CCAAT/enhancer-binding proteins (C/EBPs) in keratinocyte and sebocyte differentiation suggest that both families of transcription factors closely interact in the skin. Initial characterization of the mouse PPARbeta promoter revealed an AP-1 site that is crucial for the regulation of PPARbeta expression in response to inflammatory cytokines in the skin. We now present evidence for a novel regulatory mechanism of the expression of the PPARbeta gene by which two members of the C/EBP family of transcription factors inhibit its basal promoter activity in mouse keratinocytes. We first demonstrate that C/EBPalpha and C/EBPbeta, but not C/EBPdelta, inhibit the expression of PPARbeta through the recruitment of a transcriptional repressor complex containing HDAC-1 to a specific C/EBP binding site on the PPARbeta promoter. Consistent with this repression, the expression patterns of PPARbeta and C/EBPs are mutually exclusive in keratinocytes of the interfollicular epidermis and hair follicles in mouse developing skin. This work reveals the importance of the regulatory interplay between PPARbeta and C/EBP transcription factors in the control of proliferation and differentiation in this organ. Such insights are crucial for the understanding of the molecular control regulating the balance between proliferation and differentiation in many cell types including keratinocytes.

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.