Grigory Kolesov

MIPS Arabidopsis thaliana Database (MAtDB): an integrated biological knowledge resource based on the first complete plant genome

Arabidopsis thaliana is the first plant for which the complete genome has been sequenced and published. Annotation of complex eukaryotic genomes requires more than the assignment of genetic elements to the sequence. Besides completing the list of genes, we need to discover their cellular roles, their regulation and their interactions in order to understand the workings of the whole plant. The MIPS Arabidopsis thaliana Database (MAtDB; http://mips.gsf.de/proj/thal/db) started out as a repository for genome sequence data in the European Scientists Sequencing Arabidopsis (ESSA) project and the Arabidopsis Genome Initiative. Our aim is to transform MAtDB into an integrated biological knowledge resource by integrating diverse data, tools, query and visualization capabilities and by creating a comprehensive resource for Arabidopsis as a reference model for other species, including crop plants.

The PEDANT genome database

The PEDANT genome database (http://pedant.gsf.de) provides exhaustive automatic analysis of genomic sequences by a large variety of established bioinformatics tools through a comprehensive Web-based user interface. One hundred and seventy seven completely sequenced and unfinished genomes have been processed so far, including large eukaryotic genomes (mouse, human) published recently. In this contribution, we describe the current status of the PEDANT database and novel analytical features added to the PEDANT server in 2002. Those include: (i) integration with the BioRS data retrieval system which allows fast text queries, (ii) pre-computed sequence clusters in each complete genome, (iii) a comprehensive set of tools for genome comparison, including genome comparison tables and protein function prediction based on genomic context, and (iv) computation and visualization of protein-protein interaction (PPI) networks based on experimental data. The availability of functional and structural predictions for 650 000 genomic proteins in well organized form makes PEDANT a useful resource for both functional and structural genomics.

Snipping polymorphisms from large EST collections in barley (Hordeum vulgare L.)

The public EST (expressed sequence tag) databases represent an enormous but heterogeneous repository of sequences, including many from a broad selection of plant species and a wide range of distinct varieties. The significant redundancy within large EST collections makes them an attractive resource for rapid pre-selection of candidate sequence polymorphisms. Here we present a strategy that allows rapid identification of candidate SNPs in barley (Hordeum vulgare L.) using publicly available EST databases. Analysis of 271,630 EST sequences from different cDNA libraries, representing 23 different barley varieties, resulted in the generation of 56,302 tentative consensus sequences. In all, 8171 of these unigene sequences are members of clusters with six or more ESTs. By applying a novel SNP detection algorithm (SNiPpER) to these sequences, we identified 3069 candidate inter-varietal SNPs. In order to verify these candidate SNPs, we selected a small subset of 63 present in 36 ESTs. Of the 63 SNPs selected, we were able to validate 54 (86%) using a direct sequencing approach. For further verification, 28 ESTs were mapped to distinct loci within the barley genome. The polymorphism information content (PIC) and nucleotide diversity (π) values of the SNPs identified by the SNiPpER algorithm are significantly higher than those that were obtained by random sequencing. This demonstrates the efficiency of our strategy for SNP identification and the cost-efficient development of EST-based SNP-markers.

How gene order is influenced by the biophysics of transcription regulation

What are the forces that shape the structure of prokaryotic genomes: the order of genes, their proximity, and their orientation? Coregulation and coordinated horizontal gene transfer are believed to promote the proximity of functionally related genes and the formation of operons. However, forces that influence the structure of the genome beyond the level of a single operon remain unknown. Here, we show that the biophysical mechanism by which regulatory proteins search for their sites on DNA can impose constraints on genome structure. Using simulations, we demonstrate that rapid and reliable gene regulation requires that the transcription factor (TF) gene be close to the site on DNA the TF has to bind, thus promoting the colocalization of TF genes and their targets on the genome. We use parameters that have been measured in recent experiments to estimate the relevant length and times scales of this process and demonstrate that the search for a cognate site may be prohibitively slow if a TF has a low copy number and is not colocalized. We also analyze TFs and their sites in a number of bacterial genomes, confirm that they are colocalized significantly more often than expected, and show that this observation cannot be attributed to the pressure for coregulation or formation of selfish gene clusters, thus supporting the role of the biophysical constraint in shaping the structure of prokaryotic genomes. Our results demonstrate how spatial organization can influence timing and noise in gene expression.

Protein knot server: detection of knots in protein structures.

KNOTS (http://knots.mit.edu) is a web server that detects knots in protein structures. Several protein structures have been reported to contain intricate knots. The physiological role of knots and their effect on folding and evolution is an area of active research. The user submits a PDB id or uploads a 3D protein structure in PDB or mmCIF format. The current implementation of the server uses the Alexander polynomial to detect knots. The results of the analysis that are presented to the user are the location of the knot in the structure, the type of the knot and an interactive visualization of the knot. The results can also be downloaded and viewed offline. The server also maintains a regularly updated list of known knots in protein structures.

GeneNet: a gene network database and its automated visualization.

MOTIVATION Gene networks that provide the regulation of physiological processes are the basic feature of organisms. Information regarding the regulation of gene expression and signal transduction pathways is increasing rapidly. However, the information is hard to formalize and systematize. Ways and means for automated visualization of the gene networks based on their formalized description are needed. RESULTS The object-oriented database GeneNet and the software for its automated visualization have been developed. The main principles of a formalized description of the gene network have been worked out. Antiviral response and erythropoiesis are provided as examples to show how this is achieved. The GeneNet graphical user interface written in Java provides automated generation of the gene network diagrams and allows visualization and exploration of the GeneNet database through the Internet. A system of filters allows the selection of particular components of the network for visualization. AVAILABILITY The GeneNet database and its graphical user interface are available at http://wwwmgs.bionet.nsc.ru/systems/MGL/GeneN et/ CONTACT fedor@bionet.nsc.ru

SNAPper: gene order predicts gene function

UNLABELLED SNAPper is a network service for predicting gene function based on the conservation of gene order. AVAILABILITY The SNAPper server is available at http://pedant.gsf.de/snapper. SNAPper-based functional predictions will soon be offered as part of the PEDANT genome analysis server http://pedant.gsf.de.

Using Evolutionary Information to Find Specificity-Determining and Co-evolving Residues

Intricate networks of protein interactions rely on the ability of a protein to recognize its targets: other proteins, ligands, and sites on DNA and RNA. To recognize other molecules, it was suggested that a protein uses a small set of specificity-determining residues (SDRs). How can one find these residues in proteins and distinguish them from other functionally important amino acids? A number of bioinformatics methods to predict SDRs have been developed in recent years. These methods use genomic information and multiple sequence alignments to identify positions exhibiting a specific pattern of conservation and variability. The challenge is to delineate the evolutionary pattern of SDRs from that of the active site residues and the residues responsible for formation of the proteins structure. The phylogenetic history of a protein family makes such analysis particularly hard. Here we present two methods for finding the SDRs and the co-evolving residues (CERs) in proteins. We use a Monte Carlo approach for statistical inference, allowing us to reveal specific evolutionary patterns of SDRs and CERs. We apply these methods to study specific recognition in the bacterial two-component system and in the class Ia aminoacyl-tRNA synthetases. Our results agree well with structural information and the experimental analyses of these systems. Our results point at the complex and distinct patterns characteristic of the evolution of specificity in these systems.

Preparation, Applications, and Digital Simulation of Carbon Interdigitated Array Electrodes

Carbon interdigitated array (IDA) electrodes with features sizes down to 1.2 μm were fabricated by controlled pyrolysis of patterned photoresist. Cyclic voltammetry of reversible redox species produced the expected steady-state currents. The collection efficiency depends on the IDA electrode spacing, which ranged from around 2.7 to 16.5 μm, with the smaller dimensions achieving higher collection efficiencies of up to 98%. The signal amplification because of redox cycling makes it possible to detect species at relatively low concentrations (10(-5) molar) and the small spacing allows detection of transient electrogenerated species with much shorter lifetimes (submillisecond). Digital simulation software that accounts for both the width and height of electrode elements as well as the electrode spacing was developed to model the IDA electrode response. The simulations are in quantitative agreement with experimental data for both a simple fast one electron redox reaction and an electron transfer with a following chemical reaction at the IDAs with larger gaps whereas currents measured for the smallest IDA electrodes, that were larger than the simulated currents, are attributed to convection from induced charge electrokinetic flow.

Electrogenerated Chemiluminescence of BODIPY, Ru(bpy)32+, and 9,10-Diphenylanthracene Using Interdigitated Array Electrodes

Interdigitated array electrodes (IDAs) were used to produce steady-state electrogenerated chemiluminescence (ECL) by annihilation of oxidized and reduced forms of a substituted boron dipyrromethene (BODIPY) dye, 9,10-diphenylanthracene (DPA), and ruthenium(II) tris(bypiridine) (Ru(bpy)3(2+)). Digital simulations were in good agreement with the experimentally obtained currents and light outputs. Coreactant experiments, using tri-n-propylamine and benzoyl peroxide as a sacrificial homogeneous reductant or oxidant, show currents corresponding to electrode reactions of the dyes and not the oxidation or reduction of the coreactants. The results show that interdigitated arrays can produce stable ECL where the light intensity is magnified due to the larger currents as a consequence of feedback between generator and collector electrodes in the IDA. The light output for ECL is around 100 times higher than that obtained with regular planar electrodes with similar area.