ExPlain: finding upstream drug targets in disease gene regulatory networks.

Different signal transduction pathways leading to the activation of transcription factors (TFs) converge at key molecules that master the regulation of many cellular processes. Such crossroads of signalling networks often appear as "Achilles Heels" causing a disease when not functioning properly. Novel computational tools are needed for analysis of the gene expression data in the context of signal transduction and gene regulatory pathways and for identification of the key nodes in the networks. An integrated computational system, ExPlain (www.biobase.de) was developed for causal interpretation of gene expression data and identification of key signalling molecules. The system utilizes data from two databases (TRANSFAC and TRANSPATH) and integrates two programs: (1) Composite Module Analyst (CMA) analyses 5'-upstream regions of co-expressed genes and applies a genetic algorithm to reveal composite modules (CMs) consisting of co-occurring single TF binding sites and composite elements; (2) ArrayAnalyzer is a fast network search engine that analyses signal transduction networks controlling the activities of the corresponding TFs and seeks key molecules responsible for the observed concerted gene activation. ExPlain system was applied to microarray data on inflammatory bowel diseases (IBD). The results obtained suggest a number of highly interesting biological hypotheses about molecular mechanisms of pathological genetic disregulation.

Composite Module Analyst: identification of transcription factor binding site combinations using genetic algorithm.

Composite Module Analyst (CMA) is a novel software tool aiming to identify promoter-enhancer models based on the composition of transcription factor (TF) binding sites and their pairs. CMA is closely interconnected with the TRANSFAC database. In particular, CMA uses the positional weight matrix (PWM) library collected in TRANSFAC and therefore provides the possibility to search for a large variety of different TF binding sites. We model the structure of the long gene regulatory regions by a Boolean function that joins several local modules, each consisting of co-localized TF binding sites. Having as an input a set of co-regulated genes, CMA builds the promoter model and optimizes the parameters of the model automatically by applying a genetic-regression algorithm. We use a multicomponent fitness function of the algorithm which includes several statistical criteria in a weighted linear function. We show examples of successful application of CMA to a microarray data on transcription profiling of TNF-alpha stimulated primary human endothelial cells. The CMA web server is freely accessible at http://www.gene-regulation.com/pub/programs/cma/CMA.html. An advanced version of CMA is also a part of the commercial system ExPlaintrade mark (www.biobase.de) designed for causal analysis of gene expression data.

Composite Module Analyst: a fitness-based tool for identification of transcription factor binding site combinations.

MOTIVATION: Functionally related genes involved in the same molecular-genetic, biochemical or physiological process are often regulated coordinately. Such regulation is provided by precisely organized binding of a multiplicity of special proteins [transcription factors (TFs)] to their target sites (cis-elements) in regulatory regions of genes. Cis-element combinations provide a structural basis for the generation of unique patterns of gene expression. RESULTS: Here we present a new approach for defining promoter models based on the composition of TF binding sites and their pairs. We utilize a multicomponent fitness function for selection of the promoter model that fits best to the observed gene expression profile. We demonstrate examples of successful application of the fitness function with the help of a genetic algorithm for the analysis of functionally related or co-expressed genes as well as testing on simulated and permutated data. AVAILABILITY: The CMA program is freely available for non-commercial users. URL http://www.gene-regulation.com/pub/programs.html#CMAnalyst. It is also a part of the commercial system ExPlain (www.biobase.de) designed for causal analysis of gene expression data..

TRANSFAC and its module TRANSCompel: transcriptional gene regulation in eukaryotes.

The TRANSFAC database on transcription factors, their binding sites, nucleotide distribution matrices and regulated genes as well as the complementing database TRANSCompel on composite elements have been further enhanced on various levels. A new web interface with different search options and integrated versions of Match and Patch provides increased functionality for TRANSFAC. The list of databases which are linked to the common GENE table of TRANSFAC and TRANSCompel has been extended by: Ensembl, UniGene, EntrezGene, HumanPSD and TRANSPRO. Standard gene names from HGNC, MGI and RGD, are included for human, mouse and rat genes, respectively. With the help of InterProScan, Pfam, SMART and PROSITE domains are assigned automatically to the protein sequences of the transcription factors. TRANSCompel contains now, in addition to the COMPEL table, a separate table for detailed information on the experimental EVIDENCE on which the composite elements are based. Finally, for TRANSFAC, in respect of data growth, in particular the gain of Drosophila transcription factor binding sites (by courtesy of the Drosophila DNase I footprint database) and of Arabidopsis factors (by courtesy of DATF, Database of Arabidopsis Transcription Factors) has to be stressed. The here described public releases, TRANSFAC 7.0 and TRANSCompel 7.0, are accessible under http://www.gene-regulation.com/pub/databases.html.

Recognition of multiple patterns in unaligned sets of sequences: comparison of kernel clustering method with other methods.

MOTIVATION: Transcription factor binding sites often differ significantly in their primary sequence and can hardly be aligned. Often one set of sites can contain several subsets of sequences that follow not just one but several different patterns. There is a need for sensitive methods to reveal multiple patterns in unaligned sets of sequences. RESULTS: We developed a novel method for analysis of unaligned sets of sequences based on kernel estimation. The method is able to reveal 'multiple local patterns'-a set of weight matrices. Every weight matrix characterizes a pattern that can be found in a significant subset of sequences under analysis. The method developed has been compared with several other methods of pattern discovery such as Gibbs sampling, MEME, CONSENSUS, MULTIPROFILER and PROJECTION. The kernel method showed the best performance in terms of how close the revealed weight matrices are to the original ones. We applied the kernel method to analyze three samples of promoters (cell-cycle, T-cells and muscle-specific). We compared the multiple patterns revealed with the TRANSFAC library of weight matrices and found a strong similarity to several weight matrices for transcription factors known to be involved in the mentioned specific gene regulation. AVAILABILITY: The program is available for on-line use at: http://www.biobase.de/cgi-bin/biobase/cbs2/bin/template.cgi?template=cbscall.html

The Helmholtz Network for Bioinformatics: an integrative web portal for bioinformatics resources.

SUMMARY: The Helmholtz Network for Bioinformatics (HNB) is a joint venture of eleven German bioinformatics research groups that offers convenient access to numerous bioinformatics resources through a single web portal. The 'Guided Solution Finder' which is available through the HNB portal helps users to locate the appropriate resources to answer their queries by employing a detailed, tree-like questionnaire. Furthermore, automated complex tool cascades ('tasks'), involving resources located on different servers, have been implemented, allowing users to perform comprehensive data analyses without the requirement of further manual intervention for data transfer and re-formatting. Currently, automated cascades for the analysis of regulatory DNA segments as well as for the prediction of protein functional properties are provided. AVAILABILITY: The HNB portal is available at http://www.hnbioinfo.de

[Modeling dynamics of gene net, regulating the cell cycle in mammalian cells]. / Modelirovanie dinamiki gennykh setei, reguliruiushchikh kletochnyi tsikl v kletakh mlekopitaiushchikh.

The study of the molecular mechanisms determining cellular programs of proliferation, differentiation, and apoptosis is currently attracting much attention. Recent studies have demonstrated that the system of cell-cycle control based on the transcriptional regulation of the expression of specific genes is responsible for the transition between programs. These groups of functionally connected genes from so-called gene networks characterized by numerous feedbacks and a complex behavioral dynamics. Computer simulation methods have been applied to studying the dynamics of gene networks regulating the cell cycle of vertebrates. The data on the regulation of the key genes obtained from the CYCLE-TRRD database have been used as a basis to construct gene networks of different degrees of complexity controlling the G1/S transition, one of the most important stages of the cell cycle. The behavior dynamics of the model constructed has been analyzed. Two qualitatively different functional modes of the system has been obtained. It has also been shown that the transition between these modes depends on the duration of the proliferation signal. It has also been demonstrated that the additional feedback from factor E2F to genes c-fos and c-jun, which was predicted earlier based on the computer analysis of promoters, plays an important role in the transition of the cell to the S phase.

MATCH: A tool for searching transcription factor binding sites in DNA sequences.

Match is a weight matrix-based tool for searching putative transcription factor binding sites in DNA sequences. Match is closely interconnected and distributed together with the TRANSFAC database. In particular, Match uses the matrix library collected in TRANSFAC and therefore provides the possibility to search for a great variety of different transcription factor binding sites. Several sets of optimised matrix cut-off values are built in the system to provide a variety of search modes of different stringency. The user may construct and save his/her specific user profiles which are selected subsets of matrices including default or user-defined cut-off values. Furthermore a number of tissue-specific profiles are provided that were compiled by the TRANSFAC team. A public version of the Match tool is available at: http://www.gene-regulation.com/pub/programs.html#match. The same program with a different web interface can be found at http://compel.bionet.nsc.ru/Match/Match.html. An advanced version of the tool called Match Professional is available at http://www.biobase.de.

TRANSFAC: transcriptional regulation, from patterns to profiles.

The TRANSFAC database on eukaryotic transcriptional regulation, comprising data on transcription factors, their target genes and regulatory binding sites, has been extended and further developed, both in number of entries and in the scope and structure of the collected data. Structured fields for expression patterns have been introduced for transcription factors from human and mouse, using the CYTOMER database on anatomical structures and developmental stages. The functionality of Match, a tool for matrix-based search of transcription factor binding sites, has been enhanced. For instance, the program now comes along with a number of tissue-(or state-)specific profiles and new profiles can be created and modified with Match Profiler. The GENE table was extended and gained in importance, containing amongst others links to LocusLink, RefSeq and OMIM now. Further, (direct) links between factor and target gene on one hand and between gene and encoded factor on the other hand were introduced. The TRANSFAC public release is available at http://www.gene-regulation.com. For yeast an additional release including the latest data was made available separately as TRANSFAC Saccharomyces Module (TSM) at http://transfac.gbf.de. For CYTOMER free download versions are available at http://www.biobase.de:8080/index.html.

Evaluation of sequence motifs found in scaffold/matrix-attached regions (S/MARs).

Based on the contents of the database S/MARt DB, the most comprehensive data collection of scaffold/matrix-attached regions (S/MARs) publicly available thus far, we initiated a systematic evaluation of the stored data. By analyzing the 245 S/MAR sequences presently described in this database, we found that the S/MARs contained in this collection are generally AT-rich, with certain significant exceptions. Comparative analyses showed that most of the AT-rich motifs which were found to be enriched in S/MARs are also enriched in randomized S/MAR sequences of the same AT content. Some sequence patterns previously suggested to be characteristic for S/MARs were also investigated, among them potential binding sites for homeodomain transcription factors. Even though hexanucleotides containing the core motif of homeodomain factors were frequently observed in S/MARs, only a few potential binding sites for these factors were found enriched when compared with regulatory regions or exon sequences. All our analyses indicated that, on average, the observed frequency of motifs in S/MAR elements is largely influenced by the AT content. Our results can serve as a guideline for further improvements in the definition of S/MARs, which are now believed to constitute the functional coordinate system for genomic regulatory regions.

Automatic annotation of genomic regulatory sequences by searching for composite clusters.

A new method was developed for revealing of composite clusters of cis-elements in promoters of eukaryotic genes that are functionally related or coexpressed. A software system "ClusterScan" have been created that enables: (i) to train system on representative samples of promoters to reveal cis-elements that tend to cluster, (ii) to train system on a number of samples of functionally related promoters to identify functionally coupled transcription factors; (iii) to provide tools for searching of this clusters in genomic sequences to identify and functionally characterize regulatory regions in genome. A number of training samples of different functional and structural groups of promoters were analysed. Search for composite clusters in human chromosomes 21 and 22 reveals a number of interesting examples. Finally, a decision tree system was constructed to classify promoters of several functionally related gene groups. The decision tree system enables to identify new promoters and computationally predict their possible function.

The TRANSPATH signal transduction database: a knowledge base on signal transduction networks.

UNLABELLED: TRANSPATH is an information system on gene-regulatory pathways, and an extension module to the TRANSFAC database system (Wingender et al., Nucleic Acids Res., 28, 316-319, 2000). It focuses on pathways involved in the regulation of transcription factors in different species, mainly human, mouse and rat. Elements of the relevant signal transduction pathways like complexes, signaling molecules, and their states are stored together with information about their interaction in an object-oriented database. The database interface provides clickable maps and automatically generated pathway cascades as additional ways to explore the data. All information is validated with references to the original publications. Also, references to other databases are provided (TRANSFAC, SWISS-PROT, EMBL, PubMed and others). AVAILABILITY: The database is available over (http://transpath.gbf.de) for interactive perusal. As an exchange format for the data, eXtensible Markup Language (XML) flatfiles and a Document Type Definition (DTD) are provided.

Computer-assisted identification of cell cycle-related genes: new targets for E2F transcription factors.

The processes that take place during development and differentiation are directed through coordinated regulation of expression of a large number of genes. One such gene regulatory network provides cell cycle control in eukaryotic organisms. In this work, we have studied the structural features of the 5' regulatory regions of cell cycle-related genes. We developed a new method for identifying composite substructures (modules) in regulatory regions of genes consisting of a binding site for a key transcription factor and additional contextual motifs: potential targets for other transcription factors that may synergistically regulate gene transcription. Applying this method to cell cycle-related promoters, we created a program for context-specific identification of binding sites for transcription factors of the E2F family which are key regulators of the cell cycle. We found that E2F composite modules are found at a high frequency and in close proximity to the start of transcription in cell cycle-related promoters in comparison with other promoters. Using this information, we then searched for E2F sites in genomic sequences with the goal of identifying new genes which play important roles in controlling cell proliferation, differentiation and apoptosis. Using a chromatin immunoprecipitation assay, we then experimentally verified the binding of E2F in vivo to the promoters predicted by the computer-assisted methods. Our identification of new E2F target genes provides new insight into gene regulatory networks and provides a framework for continued analysis of the role of contextual promoter features in transcriptional regulation. The tools described are available at http://compel.bionet.nsc.ru/FunSite/SiteScan.html.

Database-assisted promoter analysis.

The analysis of regulatory sequences is greatly facilitated by database-assisted bioinformatic approaches. The TRANSFAC database contains information on transcription factors and their origins, functional properties and sequence-specific binding activities. Software tools enable us to screen the database with a given DNA sequence for interacting transcription factors. If a regulatory function is already attributed to this sequence then the database-assisted identification of binding sites for proteins or protein classes and subsequent experimental verification might establish functionally relevant sites within this sequence. The binding transcription factors and interacting factors might already be present in the database.

The TRANSFAC system on gene expression regulation.

The TRANSFAC database on transcription factors and their DNA-binding sites and profiles (http://www.gene-regulation.de/) has been quantitatively extended and supplemented by a number of modules. These modules give information about pathologically relevant mutations in regulatory regions and transcription factor genes (PathoDB), scaffold/matrix attached regions (S/MARt DB), signal transduction (TRANSPATH) and gene expression sources (CYTOMER). Altogether, these distinct database modules constitute the TRANSFAC system. They are accompanied by a number of program routines for identifying potential transcription factor binding sites or for localizing individual components in the regulatory network of a cell.

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.

COMPEL: a database on composite regulatory elements providing combinatorial transcriptional regulation.

COMPEL is a database on composite regulatory elements, the basic structures of combinatorial regulation. Composite regulatory elements contain two closely situated binding sites for distinct transcription factors and represent minimal functional units providing combinatorial transcriptional regulation. Both specific factor-DNA and factor-factor interactions contribute to the function of composite elements (CEs). Information about the structure of known CEs and specific gene regulation achieved through such CEs appears to be extremely useful for promoter prediction, for gene function prediction and for applied gene engineering as well. The structure of the relational model of COMPEL is determined by the concept of molecular structure and regulatory role of CEs. Based on the set of a particular CE, a program has been developed for searching potential CEs in gene regulatory regions. WWW search and browse routines were developed for COMPEL release 3.0. The COMPEL database equipped with the search and browse tools is available at http://compel.bionet.nsc.ru/. The program for prediction of potential CEs of NFAT type is available at http://compel.bionet.nsc. ru/FunSite.html and http://transfac.gbf.de/dbsearch/funsitep/ s_comp.html

TRANSFAC: an integrated system for gene expression regulation.

TRANSFAC is a database on transcription factors, their genomic binding sites and DNA-binding profiles (http://transfac.gbf.de/TRANSFAC/). Its content has been enhanced, in particular by information about training sequences used for the construction of nucleotide matrices as well as by data on plant sites and factors. Moreover, TRANSFAC has been extended by two new modules: PathoDB provides data on pathologically relevant mutations in regulatory regions and transcription factor genes, whereas S/MARt DB compiles features of scaffold/matrix attached regions (S/MARs) and the proteins binding to them. Additionally, the databases TRANSPATH, about signal transduction, and CYTOMER, about organs and cell types, have been extended and are increasingly integrated with the TRANSFAC data sources.

Integrated databases and computer systems for studying eukaryotic gene expression.

MOTIVATION: The goal of the work was to develop a WWW-oriented computer system providing a maximal integration of informational and software resources on the regulation of gene expression and navigation through them. Rapid growth of the variety and volume of information accumulated in the databases on regulation of gene expression necessarily requires the development of computer systems for automated discovery of the knowledge that can be further used for analysis of regulatory genomic sequences. RESULTS: The GeneExpress system developed includes the following major informational and software modules: (1) Transcription Regulation (TRRD) module, which contains the databases on transcription regulatory regions of eukaryotic genes and TRRD Viewer for data visualization; (2) Site Activity Prediction (ACTIVITY), the module for analysis of functional site activity and its prediction; (3) Site Recognition module, which comprises (a) B-DNA-VIDEO system for detecting the conformational and physicochemical properties of DNA sites significant for their recognition, (b) Consensus and Weight Matrices (ConsFrec) and (c) Transcription Factor Binding Sites Recognition (TFBSR) systems for detecting conservative contextual regions of functional sites and their recognition; (4) Gene Networks (GeneNet), which contains an object-oriented database accumulating the data on gene networks and signal transduction pathways, and the Java-based Viewer for exploration and visualization of the GeneNet information; (5) mRNA Translation (Leader mRNA), designed to analyze structural and contextual properties of mRNA 5'-untranslated regions (5'-UTRs) and predict their translation efficiency; (6) other program modules designed to study the structure-function organization of regulatory genomic sequences and regulatory proteins. AVAILABILITY: GeneExpress is available at http://wwwmgs.bionet.nsc. ru/systems/GeneExpress/ and the links to the mirror site(s) can be found at http://wwwmgs.bionet.nsc.ru/mgs/links/mirrors.html+ ++.

Recognition of NFATp/AP-1 composite elements within genes induced upon the activation of immune cells.

Composite elements are regulatory modules of promoters or enhancers that consist of binding sites of two different but synergizing transcription factors. A well-studied example is nuclear factors of activated T-cell (NFAT) sites which are composite elements of a NFATp/c and an activating protein 1 (AP-1) binding site. We have developed a computational approach to identify potential NFAT target genes which (a) comprises an improved method to scan for individual NFAT composite elements; (b) considers positional effects relative to transcription start sites; and (c) involves cluster analysis of potential NFAT composite elements. All three steps progressively helpX?ed to discriminate T-cell-specific promoter sequences against other functional regions (coding and intronic sequences) of the same genes, against promoters of muscle-specific genes or against random sequences. Using this approach, we identified potential NFAT composite elements in promoters of cytokine genes and their receptors as well as in promoters of genes for AP-1 family members, Ca2+-binding proteins and some other components of the regulatory network operating in activated T-cells and other immune cells. The method developed can be adapted to characterize and identify other composite elements as well. The program for recognition NFAT composite elements is available through the World Wide Web (http://compel.bionet.nsc.ru/FunSite/CompelScan. html and http://transfac.gbf.de/dbsearch/funsitep/s _comp.html).