MRS glucose mapping and PET joining forces: re-evaluation of the lumped constant in the rat brain under isoflurane anaesthesia.

Although numerous positron emission tomography (PET) studies with (18) F-fluoro-deoxyglucose (FDG) have reported quantitative results on cerebral glucose kinetics and consumption, there is a large variation between the absolute values found in the literature. One of the underlying causes is the inconsistent use of the lumped constants (LCs), the derivation of which is often based on multiple assumptions that render absolute numbers imprecise and errors hard to quantify. We combined a kinetic FDG-PET study with magnetic resonance spectroscopic imaging (MRSI) of glucose dynamics in Sprague-Dawley rats to obtain a more comprehensive view of brain glucose kinetics and determine a reliable value for the LC under isoflurane anaesthesia. Maps of Tmax /CMRglc derived from MRSI data and Tmax determined from PET kinetic modelling allowed to obtain an LC-independent CMRglc . The LC was estimated to range from 0.33 ± 0.07 in retrosplenial cortex to 0.44 ± 0.05 in hippocampus, yielding CMRglc between 62 ± 14 and 54 ± 11 μmol/min/100 g, respectively. These newly determined LCs for four distinct areas in the rat brain under isoflurane anaesthesia provide means of comparing the growing amount of FDG-PET data available from translational studies.

Binding free energy differences in a TCR-peptide-MHC complex induced by a peptide mutation: a simulation analysis.

Recognition by the T-cell receptor (TCR) of immunogenic peptides presented by class I major histocompatibility complexes (MHCs) is the determining event in the specific cellular immune response against virus-infected cells or tumor cells. It is of great interest, therefore, to elucidate the molecular principles upon which the selectivity of a TCR is based. These principles can in turn be used to design therapeutic approaches, such as peptide-based immunotherapies of cancer. In this study, free energy simulation methods are used to analyze the binding free energy difference of a particular TCR (A6) for a wild-type peptide (Tax) and a mutant peptide (Tax P6A), both presented in HLA A2. The computed free energy difference is 2.9 kcal/mol, in good agreement with the experimental value. This makes possible the use of the simulation results for obtaining an understanding of the origin of the free energy difference which was not available from the experimental results. A free energy component analysis makes possible the decomposition of the free energy difference between the binding of the wild-type and mutant peptide into its components. Of particular interest is the fact that better solvation of the mutant peptide when bound to the MHC molecule is an important contribution to the greater affinity of the TCR for the latter. The results make possible identification of the residues of the TCR which are important for the selectivity. This provides an understanding of the molecular principles that govern the recognition. The possibility of using free energy simulations in designing peptide derivatives for cancer immunotherapy is briefly discussed.

Heterologously expressed Staphylococcus aureus fibronectin-binding proteins are sufficient for invasion of host cells.

Staphylococcus aureus invasion of mammalian cells, including epithelial, endothelial, and fibroblastic cells, critically depends on fibronectin bridging between S. aureus fibronectin-binding proteins (FnBPs) and the host fibronectin receptor integrin alpha(5)beta(1) (B. Sinha et al., Cell. Microbiol. 1:101-117, 1999). However, it is unknown whether this mechanism is sufficient for S. aureus invasion. To address this question, various S. aureus adhesins (FnBPA, FnBPB, and clumping factor [ClfA]) were expressed in Staphylococcus carnosus and Lactococcus lactis subsp. cremoris. Both noninvasive gram-positive microorganisms are genetically distinct from S. aureus, lack any known S. aureus surface protein, and do not bind fibronectin. Transformants of S. carnosus and L. lactis harboring plasmids coding for various S. aureus surface proteins (FnBPA, FnBPB, and ClfA) functionally expressed adhesins (as determined by bacterial clumping in plasma, specific latex agglutination, Western ligand blotting, and binding to immobilized and soluble fibronectin). FnBPA or FnBPB but not of ClfA conferred invasiveness to S. carnosus and L. lactis. Invasion of 293 cells by transformants was comparable to that of strongly invasive S. aureus strain Cowan 1. Binding of soluble and immobilized fibronectin paralleled invasiveness, demonstrating that the amount of accessible surface FnBPs is rate limiting. Thus, S. aureus FnBPs confer invasiveness to noninvasive, apathogenic gram-positive cocci. Furthermore, FnBP-coated polystyrene beads were internalized by 293 cells, demonstrating that FnBPs are sufficient for invasion of host cells without the need for (S. aureus-specific) coreceptors.

A role for the transcriptional repressor ICER in the molecular pathogenesis of type 2 diabetes

AbstractType 2 diabetes (T2D) is a metabolic disease which affects more than 200 millions people worldwide. The progression of this affection reaches nowadays epidemic proportions, owing to the constant augmentation in the frequency of overweight, obesity and sedentary. The pathogenesis of T2D is characterized by reduction in the action of insulin on its target tissues - an alteration referred as insulin resistance - and pancreatic β-cell dysfunction. This latter deterioration is defined by impairment in insulin biosynthesis and secretion, and a loss of β-cell mass by apoptosis. Environmental factors related to T2D, such as chronic elevation in glucose and free fatty acids levels, inflammatory cytokines and pro-atherogenic oxidized low- density lipoproteins (LDL), contribute to the loss of pancreatic β-cell function.In this study, we have demonstrated that the transcription factor Inducible Cyclic AMP Early Repressor (ICER) participates to the progression of both β-cell dysfunction and insulin resistance. The expression of this factor is driven by an alternative promoter and ICER protein represents therefore a truncated product of the Cyclic AMP Response Element Modulator (CREM) family which lacks transactivation domain. Consequently, the transcription factor ICER acts as a passive repressor which reduces expression of genes controlled by the cyclic AMP and Cyclic AMP Response Element Binding protein (CREB) pathway.In insulin-secreting cells, the accumulation of reactive oxygen species caused by environmental factors and notably oxidized LDL - a process known as oxidative stress - induces the transcription factor ICER. This transcriptional repressor hampers the secretory capacity of β-cells by silencing key genes of the exocytotic machinery. In addition, the factor ICER reduces the expression of the scaffold protein Islet Brain 1 (IB 1 ), thereby favouring the activation of the c-Jun N-terminal Kinase (JNK) pathway. This triggering alters in turn insulin biosynthesis and survival capacities of pancreatic β-cells.In the adipose tissue of mice and human subjects suffering from obesity, the transcription factor ICER contributes to the alteration in insulin action. The loss in ICER protein in these tissues induces a constant activation of the CREB pathway and the subsequent expression of the Activating Transcription Factor 3 (ATF3). In turn, this repressor reduces the transcript levels of the glucose transporter GLUT4 and the insulin-sensitizer peptide adiponectin, thereby contributing to the diminution in insulin action.In conclusion, these data shed light on the important role of the transcriptional repressor ICER in the pathogenesis of T2D, which contributes to both alteration in β-cell function and aggravation of insulin resistance. Consequently, a better understanding of the molecular mechanisms responsible for the alterations in ICER levels is required and could lead to develop new therapeutic strategies for the treatment of T2D.RésuméLe diabète de type 2 (DT2) est une maladie métabolique qui affecte plus de 200 millions de personnes dans le monde. La progression de cette affection atteint aujourd'hui des proportions épidémiques imputables à l'augmentation rapide dans les fréquences du surpoids, de l'obésité et de la sédentarité. La pathogenèse du DT2 se caractérise par une diminution de l'action de l'insuline sur ses tissus cibles - un processus nommé insulino-résistance - ainsi qu'une dysfonction des cellules β pancréatiques sécrétrices d'insuline. Cette dernière détérioration se définit par une réduction de la capacité de synthèse et de sécrétion de l'insuline et mène finalement à une perte de la masse de cellules β par apoptose. Des facteurs environnementaux fréquemment associés au DT2, tels l'élévation chronique des taux plasmatiques de glucose et d'acides gras libres, les cytokines pro-inflammatoires et les lipoprotéines de faible densité (LDL) oxydées, contribuent à la perte de fonction des cellules β pancréatiques.Dans cette étude, nous avons démontré que le facteur de transcription « Inducible Cyclic AMP Early Repressor » (ICER) participe à la progression de la dysfonction des cellules β pancréatiques et au développement de Pinsulino-résistance. Son expression étant gouvernée par un promoteur alternatif, la protéine d'ICER représente un produit tronqué de la famille des «Cyclic AMP Response Element Modulator » (CREM), sans domaine de transactivation. Par conséquent, le facteur ICER agit comme un répresseur passif qui réduit l'expression des gènes contrôlés par la voie de l'AMP cyclique et des « Cyclic AMP Response Element Binding protein » (CREB).Dans les cellules sécrétrices d'insuline, l'accumulation de radicaux d'oxygène libres, soutenue par les facteurs environnementaux et notamment les LDL oxydées - un processus appelé stress oxydatif- induit de manière ininterrompue le facteur de transcription ICER. Ainsi activé, ce répresseur transcriptionnel altère la capacité sécrétoire des cellules β en bloquant l'expression de gènes clés de la machinerie d'exocytose. En outre, le facteur ICER favorise l'activation de la cascade de signalisation « c-Jun N- terminal Kinase » (JNK) en réduisant l'expression de la protéine « Islet Brain 1 » (IB1), altérant ainsi les fonctions de biosynthèse de l'insuline et de survie des cellules β pancréatiques.Dans le tissu adipeux des souris et des sujets humains souffrant d'obésité, le facteur de transcription ICER contribue à l'altération de la réponse à l'insuline. La disparition de la protéine ICER dans ces tissus entraîne une activation persistante de la voie de signalisation des CREB et une induction du facteur de transcription « Activating Transcription Factor 3 » (ATF3). A son tour, le répresseur ATF3 inhibe l'expression du transporteur de glucose GLUT4 et du peptide adipocytaire insulino-sensibilisateur adiponectine, contribuant ainsi à la diminution de l'action de l'insuline en conditions d'obésité.En conclusion, à la lumière de ces résultats, le répresseur transcriptionnel ICER apparaît comme un facteur important dans la pathogenèse du DT2, en participant à la perte de fonction des cellules β pancréatiques et à l'aggravation de l'insulino-résistance. Par conséquent, l'étude des mécanismes moléculaires responsables de l'altération des niveaux du facteur ICER pourrait permettre le développement de nouvelles stratégies de traitement du DT2.Résumé didactiqueL'énergie nécessaire au bon fonctionnement de l'organisme est fournie par l'alimentation, notamment sous forme de sucres (glucides). Ceux-ci sont dégradés en glucose, lequel sera distribué aux différents organes par la circulation sanguine. Après un repas, le niveau de glucose sanguin, nommé glycémie, s'élève et favorise la sécrétion d'une hormone appelée insuline par les cellules β du pancréas. L'insuline permet, à son tour, aux organes, tels le foie, les muscles et le tissu adipeux de capter et d'utiliser le glucose ; la glycémie retrouve ainsi son niveau basai.Le diabète de type 2 (DT2) est une maladie métabolique qui affecte plus de 200 millions de personnes dans le monde. Le développement de cette affection est causée par deux processus pathologiques. D'une part, les quantités d'insuline secrétée par les cellules β pancréatiques, ainsi que la survie de ces cellules sont réduites, un phénomène connu sous le nom de dysfonction des cellules β. D'autre part, la sensibilité des tissus à l'insuline se trouve diminuée. Cette dernière altération, l'insulino-résistance, empêche le transport et l'utilisation du glucose par les tissus et mène à une accumulation de ce sucre dans le sang. Cette stagnation de glucose dans le compartiment sanguin est appelée hyperglycémie et favorise l'apparition des complications secondaires du diabète, telles que les maladies cardiovasculaires, l'insuffisance rénale, la cécité et la perte de sensibilité des extrémités.Dans cette étude, nous avons démontré que le facteur ICER qui contrôle spécifiquement l'expression de certains gènes, contribue non seulement à la dysfonction des cellules β, mais aussi au développement de l'insulino-résistance. En effet, dans les cellules β pancréatiques en conditions diabétiques, l'activation du facteur ICER altère la capacité de synthèse et de sécrétion d'insuline et réduit la survie ces cellules.Dans le tissu adipeux des souris et des sujets humains souffrant d'obésité, le facteur ICER contribue à la perte de sensibilité à l'insuline. La disparition d'ICER altère l'expression de la protéine qui capte le glucose, le transoprteur GLUT4, et l'hormone adipocytaire favorisant la sensibilité à l'insuline, nommée adiponectine. Ainsi, la perte d'ICER participe à la réduction de la captation de glucose par le tissue adipeux et au développement de l'insulino-résistance au cours de l'obésité.En conclusion, à la lumière de ces résultats, le facteur ICER apparaît comme un contributeur important à la progression du DT2, en soutenant la dysfonction des cellules β pancréatiques et l'aggravation de l'insulino-résistance. Par conséquent, l'étude des mécanismes responsables de la dérégulation du facteur ICER pourrait permettre le développement de nouvelles stratégies de traitement du DT2.

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.

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.

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.