H-2Kd-restricted antigenic peptides share a simple binding motif.

We have defined structural features that are apparently important for the binding of four different, unrelated antigenic epitopes to the same major histocompatibility complex (MHC) class I molecule, H-2Kd. The four epitopes are recognized in the form of synthetic peptides by cytotoxic T lymphocytes of the appropriate specificity. By analysis of the relative potency of truncated peptides, we demonstrated that for each of the four epitopes, optimal antigenic activity was present in a peptide of 9 or 10 amino acid residues. A comparison of the relative competitor activity of the different-length peptides in a functional competition assay, as well as in a direct binding assay based on photoaffinity labeling of the Kd molecule, indicated that the enhanced potency of the peptides upon reduction in length was most likely due to a higher affinity of the shorter peptides for the Kd molecule. A remarkably simple motif that appears to be important for the specific binding of Kd-restricted peptides was identified by the analysis of peptides containing amino acid substitutions or deletions. The motif consists of two elements, a Tyr in the second position relative to the NH2 terminus and a hydrophobic residue with a large aliphatic side chain (Leu, Ile, or Val) at the COOH-terminal end of the optimal 9- or 10-mer peptides. We demonstrated that a simple peptide analogue (AYP6L) that incorporates the motif can effectively and specifically interact with the Kd molecule. Moreover, all of the additional Kd-restricted epitopes defined thus far in the literature contain the motif, and it may thus be useful for the prediction of new epitopes recognized by T cells in the context of this MHC class I molecule.

Simultaneous co-expression of memory- and effector-related genes by individual human CD8 T cells depends on antigen specificity and differentiation

Purpose/Objective: Phenotypic and functional T cell properties are usually analyzed at the level of defined cell populations. However, large differences between individual T cells may have important functional consequences. To answer this issue, we performed highly sensitive single-cell gene expression profiling, which allows the direct ex vivo characterization of individual virus- and tumor-specific T cells from healthy donors and melanoma patients. Materials and methods: HLA-A*0201-positive patients with stage III/ IV metastatic melanoma were included in a phase I clinical trial (LUD- 00-018). Patients received monthly low-dose of the Melan-AMART- 1 26_35 unmodified natural (EAAGIGILTV) or the analog A27L (ELAGIGILTV) peptides, mixed CPG and IFA. Individual effector memory CD28+ (EM28+) and EM28- tetramer-specific CD8pos T cells were sorted by flow cytometer. Following direct cell lysis and reverse transcription, the resulting cDNA was precipitated and globally amplified. Semi-quantitative PCR was used for gene expression and TCR BV repertoire analyses. Results: We have previously shown that vaccination with the natural Melan-A peptide induced T cells with superior effector functions as compared to the analog peptide optimized for enhanced HLA binding. Here we found that natural peptide vaccination induced EM28+ T cells with frequent co-expression of both memory/homing-associated genes (CD27, IL7R, EOMES, CXCR3 and CCR5) and effector-related genes (IFNG, KLRD1, PRF1 and GZMB), comparable to protective EBV- and CMV-specific T cells. In contrast, memory/homing- and effectorassociated genes were less frequently co-expressed after vaccination with the analog peptide. Conclusions: These findings reveal a previously unknown level of gene expression diversity among vaccine- and virus-specific T cells with the simultaneous co-expression of multiple memory/homing- and effector- related genes by the same cell. Such broad functional gene expression signatures within antigen-specific T cells may be critical for mounting efficient responses to pathogens or tumors. In summary, direct ex vivo high-resolution molecular characterization of individual T cells provides key insights into the processes shaping the functional properties of tumor- and virus-specific T cells.

Molecular insights for optimizing T cell receptor specificity against cancer.

Cytotoxic CD8 T cells mediate immunity to pathogens and they are able to eliminate malignant cells. Immunity to viruses and bacteria primarily involves CD8 T cells bearing high affinity T cell receptors (TCRs), which are specific to pathogen-derived (non-self) antigens. Given the thorough elimination of high affinity self/tumor-antigen reactive T cells by central and peripheral tolerance mechanisms, anti-cancer immunity mostly depends on TCRs with intermediate-to-low affinity for self-antigens. Because of this, a promising novel therapeutic approach to increase the efficacy of tumor-reactive T cells is to engineer their TCRs, with the aim to enhance their binding kinetics to pMHC complexes, or to directly manipulate the TCR-signaling cascades. Such manipulations require a detailed knowledge on how pMHC-TCR and co-receptors binding kinetics impact the T cell response. In this review, we present the current knowledge in this field. We discuss future challenges in identifying and targeting the molecular mechanisms to enhance the function of natural or TCR-affinity optimized T cells, and we provide perspectives for the development of protective anti-tumor T cell responses.

Specific binding of antigenic peptides to cell-associated MHC class I molecules.

T lymphocytes recognize antigen in the form of peptides that associate with specific alleles of class I or class II major histocompatibility (MHC) molecules. By contrast with the clear MHC allele-specific binding of peptides to purified class II molecules purified solubilized class I molecules either bind relatively poorly or show degenerate specificity. Using photo-affinity labelling, we demonstrate here the specific interaction of peptides with cell-associated MHC class I molecules and show that this involves metabolically active processes.

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.

Neutralization of Macrophage Migration Inhibitory Factor (MIF) by Fully Human Antibodies Correlates with Their Specificity for the β-Sheet Structure of MIF.

The macrophage migration inhibitory factor (MIF) is a proinflammatory cytokine that recently emerged as an attractive therapeutic target for a variety of diseases. A diverse panel of fully human anti-MIF antibodies was generated by selection from a phage display library and extensively analyzed in vitro. Epitope mapping studies identified antibodies specific for linear as well as structural epitopes. Experimental animal studies revealed that only those antibodies binding epitopes within amino acids 50-68 or 86-102 of the MIF molecule exerted protective effects in models of sepsis or contact hypersensitivity. Within the MIF protein, these two binding regions form a β-sheet structure that includes the MIF oxidoreductase motif. We therefore conclude that this β-sheet structure is a crucial region for MIF activity and a promising target for anti-MIF antibody therapy.

Molecular modeling study of the binding and signaling properties of the TCRpMHC complex

The activation of the specific immune response against tumor cells is based on the recognition by the CD8+ Cytotoxic Τ Lymphocytes (CTL), of antigenic peptides (p) presented at the surface of the cell by the class I major histocompatibility complex (MHC). The ability of the so-called T-Cell Receptors (TCR) to discriminate between self and non-self peptides constitutes the most important specific control mechanism against infected cells. The TCR/pMHC interaction has been the subject of much attention in cancer therapy since the design of the adoptive transfer approach, in which Τ lymphocytes presenting an interesting response against tumor cells are extracted from the patient, expanded in vitro, and reinfused after immunodepletion, possibly leading to cancer regression. In the last decade, major progress has been achieved by the introduction of engineered lypmhocytes. In the meantime, the understanding of the molecular aspects of the TCRpMHC interaction has become essential to guide in vitro and in vivo studies. In 1996, the determination of the first structure of a TCRpMHC complex by X-ray crystallography revealed the molecular basis of the interaction. Since then, molecular modeling techniques have taken advantage of crystal structures to study the conformational space of the complex, and understand the specificity of the recognition of the pMHC by the TCR. In the meantime, experimental techniques used to determine the sequences of TCR that bind to a pMHC complex have been used intensively, leading to the collection of large repertoires of TCR sequences that are specific for a given pMHC. There is a growing need for computational approaches capable of predicting the molecular interactions that occur upon TCR/pMHC binding without relying on the time consuming resolution of a crystal structure. This work presents new approaches to analyze the molecular principles that govern the recognition of the pMHC by the TCR and the subsequent activation of the T-cell. We first introduce TCRep 3D, a new method to model and study the structural properties of TCR repertoires, based on homology and ab initio modeling. We discuss the methodology in details, and demonstrate that it outperforms state of the art modeling methods in predicting relevant TCR conformations. Two successful applications of TCRep 3D that supported experimental studies on TCR repertoires are presented. Second, we present a rigid body study of TCRpMHC complexes that gives a fair insight on the TCR approach towards pMHC. We show that the binding mode of the TCR is correctly described by long-distance interactions. Finally, the last section is dedicated to a detailed analysis of an experimental hydrogen exchange study, which suggests that some regions of the constant domain of the TCR are subject to conformational changes upon binding to the pMHC. We propose a hypothesis of the structural signaling of TCR molecules leading to the activation of the T-cell. It is based on the analysis of correlated motions in the TCRpMHC structure. - L'activation de la réponse immunitaire spécifique dirigée contre les cellules tumorales est basée sur la reconnaissance par les Lymphocytes Τ Cytotoxiques (CTL), d'un peptide antigénique (p) présenté à la suface de la cellule par le complexe majeur d'histocompatibilité de classe I (MHC). La capacité des récepteurs des lymphocytes (TCR) à distinguer les peptides endogènes des peptides étrangers constitue le mécanisme de contrôle le plus important dirigé contre les cellules infectées. L'interaction entre le TCR et le pMHC est le sujet de beaucoup d'attention dans la thérapie du cancer, depuis la conception de la méthode de transfer adoptif: les lymphocytes capables d'une réponse importante contre les cellules tumorales sont extraits du patient, amplifiés in vitro, et réintroduits après immunosuppression. Il peut en résulter une régression du cancer. Ces dix dernières années, d'importants progrès ont été réalisés grâce à l'introduction de lymphocytes modifiés par génie génétique. En parallèle, la compréhension du TCRpMHC au niveau moléculaire est donc devenue essentielle pour soutenir les études in vitro et in vivo. En 1996, l'obtention de la première structure du complexe TCRpMHC à l'aide de la cristallographie par rayons X a révélé les bases moléculaires de l'interaction. Depuis lors, les techniques de modélisation moléculaire ont exploité les structures expérimentales pour comprendre la spécificité de la reconnaissance du pMHC par le TCR. Dans le même temps, de nouvelles techniques expérimentales permettant de déterminer la séquence de TCR spécifiques envers un pMHC donné, ont été largement exploitées. Ainsi, d'importants répertoires de TCR sont devenus disponibles, et il est plus que jamais nécessaire de développer des approches informatiques capables de prédire les interactions moléculaires qui ont lieu lors de la liaison du TCR au pMHC, et ce sans dépendre systématiquement de la résolution d'une structure cristalline. Ce mémoire présente une nouvelle approche pour analyser les principes moléculaires régissant la reconnaissance du pMHC par le TCR, et l'activation du lymphocyte qui en résulte. Dans un premier temps, nous présentons TCRep 3D, une nouvelle méthode basée sur les modélisations par homologie et ab initio, pour l'étude de propriétés structurales des répertoires de TCR. Le procédé est discuté en détails et comparé à des approches standard. Nous démontrons ainsi que TCRep 3D est le plus performant pour prédire des conformations pertinentes du TCR. Deux applications à des études expérimentales des répertoires TCR sont ensuite présentées. Dans la seconde partie de ce travail nous présentons une étude de complexes TCRpMHC qui donne un aperçu intéressant du mécanisme d'approche du pMHC par le TCR. Finalement, la dernière section se concentre sur l'analyse détaillée d'une étude expérimentale basée sur les échanges deuterium/hydrogène, dont les résultats révèlent que certaines régions clés du domaine constant du TCR sont sujettes à un changement conformationnel lors de la liaison au pMHC. Nous proposons une hypothèse pour la signalisation structurelle des TCR, menant à l'activation du lymphocyte. Celle-ci est basée sur l'analyse des mouvements corrélés observés dans la structure du TCRpMHC.

DNA binding properties of peroxisome proliferator-activated receptor subtypes on various natural peroxisome proliferator response elements. Importance of the 5'-flanking region.

The three subtypes of the peroxisome proliferator-activated receptors (PPARalpha, beta/delta, and gamma) form heterodimers with the 9-cis-retinoic acid receptor (RXR) and bind to a common consensus response element, which consists of a direct repeat of two hexanucleotides spaced by one nucleotide (DR1). As a first step toward understanding the molecular mechanisms determining PPAR subtype specificity, we evaluated by electrophoretic mobility shift assays the binding properties of the three PPAR subtypes, in association with either RXRalpha or RXRgamma, on 16 natural PPAR response elements (PPREs). The main results are as follows. (i) PPARgamma in combination with either RXRalpha or RXRgamma binds more strongly than PPARalpha or PPARbeta to all natural PPREs tested. (ii) The binding of PPAR to strong elements is reinforced if the heterodimerization partner is RXRgamma. In contrast, weak elements favor RXRalpha as heterodimerization partner. (iii) The ordering of the 16 natural PPREs from strong to weak elements does not depend on the core DR1 sequence, which has a relatively uniform degree of conservation, but correlates with the number of identities of the 5'-flanking nucleotides with respect to a consensus element. This 5'-flanking sequence is essential for PPARalpha binding and thus contributes to subtype specificity. As a demonstration of this, the PPARgamma-specific element ARE6 PPRE is able to bind PPARalpha only if its 5'-flanking region is exchanged with that of the more promiscuous HMG PPRE.

Monoclonal antibodies against idiotypic determinant(s) of the T cell receptor from HPB-ALL cells induce IL2 production in Jurkat cells without apparent evidence of binding.

Two monoclonal antibodies (mAb) directed against idiotypic determinants of the T cell receptor (anti-Ti) from HPB-ALL cells induce interleukin 2 (IL2) production in Jurkat T cells without evidence of binding to these cells as judged by fluorescence-activated cell sorter (FACS) analysis, indirect antibody-binding radioimmunoassay and direct binding studies with 125I-labeled mAb. The IL2 response induced by these mAb observed both in the presence and absence of phorbol myristate acetate was in the range of that obtained when Jurkat cells were stimulated with phytohemagglutinin or anti-T3 mAb (Leu 4). The idiotypic specificity of the two anti-HPB-ALL Ti mAb was demonstrated by several criteria. Both mAb bound specifically to HPB-ALL cells as determined by radioimmunoassay or FACS analysis but not with 8 other T cell lines. The anti-HPB-ALL Ti mAb precipitated a disulfide-linked heterodimer of 85 kDa only from 125I-labeled HPB-ALL cells and not from other cell lines tested. Incubation of HPB-ALL cells with anti-T3 abrogated the expression of T3 and induced co-modulation of the idiotypic structures detected by the two anti-HPB-ALL Ti mAb. Conversely, incubation of HPB-ALL cells with either one of the anti-Ti mAb abrogated the expression of T3 and of the idiotypic structures. Our results suggest that mAb with an apparent unique specificity for the receptor of the immunizing T cell line HPB-ALL can activate Jurkat cells by a very weak cross-reaction with these cells, which is not detectable by conventional binding tests.

Integrin-mediated adhesion and soluble ligand binding stabilize COX-2 protein levels in endothelial cells by inducing expression and preventing degradation.

Cyclooxygenase-2 (COX-2), a key enzyme in prostaglandin synthesis, is highly expressed during inflammation and cellular transformation and promotes tumor progression and angiogenesis. We have previously demonstrated that endothelial cell COX-2 is required for integrin alphaVbeta3-dependent activation of Rac-1 and Cdc-42 and for endothelial cell spreading, migration, and angiogenesis (Dormond, O., Foletti, A., Paroz, C., and Ruegg, C. (2001) Nat. Med. 7, 1041-1047; Dormond, O., Bezzi, M., Mariotti, A., and Ruegg, C. (2002) J. Biol. Chem. 277, 45838-45846). In this study, we addressed the question of whether integrin-mediated cell adhesion may regulate COX-2 expression in endothelial cells. We report that cell detachment from the substrate caused rapid degradation of COX-2 protein in human umbilical vein endothelial cells (HUVEC) independent of serum stimulation. This effect was prevented by broad inhibition of cellular proteinases and by neutralizing lysosomal activity but not by inhibiting the proteasome. HUVEC adhesion to laminin, collagen I, fibronectin, or vitronectin induced rapid COX-2 protein expression with peak levels reached within 2 h and increased COX-2-dependent prostaglandin E2 production. In contrast, nonspecific adhesion to poly-L-lysine was ineffective in inducing COX-2 expression. Furthermore, the addition of matrix proteins in solution promoted COX-2 protein expression in suspended or poly-L-lysine-attached HUVEC. Adhesion-induced COX-2 expression was strongly suppressed by pharmacological inhibition of c-Src, phosphatidylinositol 3-kinase, p38, extracellular-regulated kinase 1/2, and, to a lesser extent, protein kinase C and by the inhibition of mRNA or protein synthesis. In conclusion, this work demonstrates that integrin-mediated cell adhesion and soluble integrin ligands contribute to maintaining COX-2 steady-state levels in endothelial cells by the combined prevention of lysosomal-dependent degradation and the stimulation of mRNA synthesis involving multiple signaling pathways.

Crystal structure of yeast peroxisomal multifunctional enzyme: structural basis for substrate specificity of (3R)-hydroxyacyl-CoA dehydrogenase units.

(3R)-hydroxyacyl-CoA dehydrogenase is part of multifunctional enzyme type 2 (MFE-2) of peroxisomal fatty acid beta-oxidation. The MFE-2 protein from yeasts contains in the same polypeptide chain two dehydrogenases (A and B), which possess difference in substrate specificity. The crystal structure of Candida tropicalis (3R)-hydroxyacyl-CoA dehydrogenase AB heterodimer, consisting of dehydrogenase A and B, determined at the resolution of 2.2A, shows overall similarity with the prototypic counterpart from rat, but also important differences that explain the substrate specificity differences observed. Docking studies suggest that dehydrogenase A binds the hydrophobic fatty acyl chain of a medium-chain-length ((3R)-OH-C10) substrate as bent into the binding pocket, whereas the short-chain substrates are dislocated by two mechanisms: (i) a short-chain-length 3-hydroxyacyl group ((3R)-OH-C4) does not reach the hydrophobic contacts needed for anchoring the substrate into the active site; and (ii) Leu44 in the loop above the NAD(+) cofactor attracts short-chain-length substrates away from the active site. Dehydrogenase B, which can use a (3R)-OH-C4 substrate, has a more shallow binding pocket and the substrate is correctly placed for catalysis. Based on the current structure, and together with the structure of the 2-enoyl-CoA hydratase 2 unit of yeast MFE-2 it becomes obvious that in yeast and mammalian MFE-2s, despite basically identical functional domains, the assembly of these domains into a mature, dimeric multifunctional enzyme is very different.

Simultaneous coexpression of memory-related and effector-related genes by individual human CD8 T cells depends on antigen specificity and differentiation.

Phenotypic and functional cell properties are usually analyzed at the level of defined cell populations but not single cells. Yet, large differences between individual cells may have important functional consequences. It is likely that T-cell-mediated immunity depends on the polyfunctionality of individual T cells, rather than the sum of functions of responding T-cell subpopulations. We performed highly sensitive single-cell gene expression profiling, allowing the direct ex vivo characterization of individual virus-specific and tumor-specific T cells from healthy donors and melanoma patients. We have previously shown that vaccination with the natural tumor peptide Melan-A-induced T cells with superior effector functions as compared with vaccination with the analog peptide optimized for enhanced HLA-A*0201 binding. Here we found that natural peptide vaccination induced tumor-reactive CD8 T cells with frequent coexpression of both memory/homing-associated genes (CD27, IL7R, EOMES, CXCR3, and CCR5) and effector-related genes (IFNG, KLRD1, PRF1, and GZMB), comparable with protective Epstein-Barr virus-specific and cytomegalovirus-specific T cells. In contrast, memory/homing-associated and effector-associated genes were less frequently coexpressed after vaccination with the analog peptide. Remarkably, these findings reveal a previously unknown level of gene expression diversity among vaccine-specific and virus-specific T cells with the simultaneous coexpression of multiple memory/homing-related and effector-related genes by the same cell. Such broad functional gene expression signatures within antigen-specific T cells may be critical for mounting efficient responses to pathogens or tumors. In summary, direct ex vivo high-resolution molecular characterization of individual T cells provides key insights into the processes shaping the functional properties of tumor-specific and virus-specific T cells.

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.

Swine monoclonal antibodies of high affinity and specificity to carcinoembryonic antigen.

To avoid the exclusive use of rodent monoclonal antibodies (MAbs) in patients for the detection of tumors by immunoscintigraphy and for radioimmunotherapy, swine MAbs were produced that are directed against carcinoembryonic antigen (CEA). Spleen cells from 2 pigs immunized with purified colon carcinoma CEA were fused with a nonsecreting mouse myeloma cell line by conventional methods, except that a particularly long immunization protocol and large amounts of spleen and myeloma cells were used. Of 1,200 growing hybrids tested, 20 were found initially to produce antibodies binding to radiolabeled CEA. Seven stable clones producing anti-CEA MAbs for more than 6 months were derived from these hybrids by repeated subcloning. The pig origin of the seven MAbs was demonstrated in a solid-phase CEA enzyme immunoassay where anti-pig immunoglobin (Ig) antibodies coupled to peroxidase gave a positive reaction while anti-mouse Ig antibodies were entirely negative. All swine MAbs were of the IgG isotype. Three anti-CEA MAbs showed no cross-reactivity with granulocytes, while four others gave various degrees of reactivity with different granulocyte glycoproteins. Competitive binding to CEA performed for two purified swine MAbs showed that they recognized two different epitopes. The affinity constants measured for these two MAbs by Scatchard plot on purified CEA were high (1.2 X 10(9) and 1.2 X 10(10) liter/mol). One of the MAbs was tested in vivo for tumor localization by injection, after radiolabeling, in nude mice bearing human colon carcinoma xenograft. High ratios of tumor to normal tissue were obtained with mean values of 10.5 for intact MAbs and of 26.8 for F(ab')2 fragments of the porcine MAb. The results showed that heterofusion with this particular protocol can be used to produce swine MAbs of high affinity and specificity for a well-defined tumor marker. These reagents may have an important clinical utility, particularly in patients who became sensitized to mouse immunoglobulins.

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.