Model.it: building three dimensional DNA models from sequence data

SUMMARY: A WWW server is described for creating 3D models of canonical or bent DNA starting from sequence data. Predicted DNA trajectory is first computed based on a choice of di- and tri-nucleotide models (M.G. Munteanu et al., Trends Biochem. Sci. 23, 341-347, 1998); an atomic model is then constructed and optionally energy-minimized with constrained molecular dynamics. The data are presented as a standard PDB file, directly viewable on the user's PC using any molecule manipulation program. AVAILABILITY: The model.it server is freely available at http://www.icgeb.trieste.it/dna/ CONTACT: kristian@icgeb.trieste.it; pongor@icgeb.trieste.it SUPPLEMENTARY INFORMATION: a series of help files is available at the above address.     

The domain-server: direct prediction of protein domain-homologies from BLAST search.

RESULTS: A WWW server for protein domain homology prediction, based on BLAST search and a simple data-mining algorithm (Hegyi,H. and Pongor,S. (1993) Comput. Appl. Biosci., 9, 371-372), was constructed providing a tabulated list and a graphic plot of similarities. AVAILABILITY: http://www.icgeb.trieste.it/domain. Mirror site is available at http://sbase.abc.hu/domain. A standalone programme will be available on request. SUPPLEMENTARY INFORMATION: A series of help files is available at the above addresses.     

Geno viewer, a SAM/BAM viewer tool

The ever evolving Next Generation Sequencing technology is calling for new and innovative ways of data processing and visualization. Following a detailed survey of the current needs of researchers and service providers, the authors have developed GenoViewer: a highly user-friendly, easy-to-operate SAM/BAM viewer and aligner tool. GenoViewer enables fast and efficient NGS assembly browsing, analysis and read mapping. It is highly customized, making it suitable for a wide range of NGS related tasks. Due to its relatively simple architecture, it is easy to add specialised visualization functionalities, facilitating further customised data analysis. The software's source code is freely available; it is open for project and task-specific modifications. AVAILABILITY: The database is available for free at http://www.genoviewer.com/     

Matrix2png: a utility for visualizing matrix data

We describe a simple software tool, 'matrix2png', for creating color images of matrix data. Originally designed with the display of microarray data sets in mind, it is a general tool that can be used to make simple visualizations of matrices for use in figures, web pages, slide presentations and the like. It can also be used to generate images 'on the fly' in web applications. Both continuous-valued and discrete-valued (categorical) data sets can be displayed. Many options are available to the user, including the colors used, the display of row and column labels, and scale bars. In this note we describe some of matrix2png's features and describe some places it has been useful in the authors' work. AVAILABILITY: A simple web interface is available, and Unix binaries are available from http://microarray.cpmc.columbia.edu/matrix2png. Source code is available on request.     

Pedstrip: extracting a maximal subset of available,unrelated individuals from a pedigree

SUMMARY: Certain types of genetic analysis are simplified by assembling a collection of unrelated individuals, e.g. case-control experiments. If a family study is being performed then it will be necessary to extract subsets of unrelated, available individuals from pedigress. Our program provides an optimal method for performing this task. AVAILABILITY: The software is available, free of charge, on request from Oxagen Ltd. SUPPLEMENTARY INFORMATION: http://www.oxagen.co.uk     

Dotlet: diagonal plots in a web browser

SUMMARY: Dotlet is a program for comparing sequences by the diagonal plot method. It is designed to be platform-independent and to run in a Web browser, thus enabling the majority of researchers to use it. AVAILABILITY: The applet can be tested at http://www.isrec.isb-sib.ch/java/dotlet/ Dotlet.html, and the source code is available upon request. CONTACT: Thomas.Junier Marco.Pagni @isrec.unil.ch SUPPLEMENTARY: The full documentation about d o t l e t is available from the above URL.     

MMDB: Entrez's 3D-structure database

Three-dimensional structures are now known within most protein families and it is likely, when searching a sequence database, that one will identify a homolog of known structure. The goal of Entrez's 3D-structure database is to make structure information and the functional annotation it can provide easily accessible to molecular biologists. To this end, Entrez's search engine provides several powerful features: (i) links between databases, for example between a protein's sequence and structure; (ii) pre-computed sequence and structure neighbors; and (iii) structure and sequence/structure alignment visualization. Here, we focus on a new feature of Entrez's Molecular Modeling Database (MMDB): Graphical summaries of the biological annotation available for each 3D structure, based on the results of automated comparative analysis. MMDB is available at: http://www.ncbi.nlm.nih.gov/Entrez/structure.html.     

Virgil database for rich links (1999 update).

With so many databases available for research in the Human Genome Project, it is crucial to efficiently relate information from different resources. For that purpose, we maintain Virgil, a database of rich links for data browsing, data analysis and database interconnection. Virgil current version contains more than 40 000 rich links from five major databases: SWISS-PROT, GenBank, PDB, GDB and OMIM. Materials described in this paper are available from http://www.infobiogen.fr/services/virgil/     

Eigen-R2 for dissecting variation in high-dimensional studies

We provide a new statistical algorithm and software package called 'eigen-R(2)' for dissecting the variation of a high-dimensional biological dataset with respect to other measured variables of interest. We apply eigen-R(2) to two real-life examples and compare it with simply averaging R(2) over many features. AVAILABILITY: An R-package eigenR2 is available at http://www.genomine.org/eigenr2/ and will be made publicly available via Bioconductor.     

Codelink: an R package for analysis of GE healthcare gene expression bioarrays

MOTIVATION: Microarray-based expression profiles have become a standard methodology in any high-throughput analysis. Several commercial platforms are available, each with its strengths and weaknesses. The R platform for statistical analysis and graphics is a powerful environment for the analysis of microarray data, because it has many integrated statistical methods available as well as the specialized microarray analysis project Bioconductor. Many packages have been added in the last few years increasing the range of possible analysis. Here, we report the availability of a package for reading and analyzing data from GE Healthcare Gene Expression Bioarrays within the R environment. AVAILABILITY: The software is implemented in the R language, is open source and available for download free of charge through the Bioconductor (http://www.bioconductor.org) project.     

CGHcall: calling aberrations for array CGH tumor profiles

CGHcall achieves high calling accuracy for array CGH data by effective use of breakpoint information from segmentation and by inclusion of several biological concepts that are ignored by existing algorithms. The algorithm is validated for simulated and verified real array CGH data. By incorporating more than three classes, CGHcall improves detection of single copy gains and amplifications. Moreover, it allows effective inclusion of chromosome arm information. AVAILABILITY: An R-package (GUI), a manual and an example data set are available at http://www.few.vu.nl/~mavdwiel/CGHcall.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

TransTerm, the translational signal database, extended to include full coding sequences and untranslated regions.

TransTerm is a database of mRNA sequences and parameters useful for detecting translational control signals in general. TransTerm-98 has been expanded beyond previous years to include full coding sequences and UTRs, while retaining the original small contexts about the coding sequence start- and stop-codons. The database contains more than 130 000 non-redundant coding sequences with associated untranslated regions (UTRs) from over 450 species. This includes the complete genomes of 12 prokaryotic and one eukaryotic organism. Several coding sequence parameters are available: coding sequence length, Nc, GC3 and, when it is computable, Codon Adaptation Index (CAI). Codon usage tables and summaries of start- and stop-codon contexts are also included. TransTerm-98 has both a relational database form with a WWW interface and a flatfile format, also available by Internet browser. TransTerm is available at: http://biochem.otago.ac.nz:800/Transterm/homepage.h tml     

Align-m--a new algorithm for multiple alignment of highly divergent sequences

MOTIVATION: Multiple alignment of highly divergent sequences is a challenging problem for which available programs tend to show poor performance. Generally, this is due to a scoring function that does not describe biological reality accurately enough or a heuristic that cannot explore solution space efficiently enough. In this respect, we present a new program, Align-m, that uses a non-progressive local approach to guide a global alignment. RESULTS: Two large test sets were used that represent the entire SCOP classification and cover sequence similarities between 0 and 50% identity. Performance was compared with the publicly available algorithms ClustalW, T-Coffee and DiAlign. In general, Align-m has comparable or slightly higher accuracy in terms of correctly aligned residues, especially for distantly related sequences. Importantly, it aligns much fewer residues incorrectly, with average differences of over 15% compared with some of the other algorithms. AVAILABILITY: Align-m and the test sets are available at http://bioinformatics.vub.ac.be     

Robust estimation of the false discovery rate

MOTIVATION: Presently available methods that use p-values to estimate or control the false discovery rate (FDR) implicitly assume that p-values are continuously distributed and based on two-sided tests. Therefore, it is difficult to reliably estimate the FDR when p-values are discrete or based on one-sided tests. RESULTS: A simple and robust method to estimate the FDR is proposed. The proposed method does not rely on implicit assumptions that tests are two-sided or yield continuously distributed p-values. The proposed method is proven to be conservative and have desirable large-sample properties. In addition, the proposed method was among the best performers across a series of 'real data simulations' comparing the performance of five currently available methods. AVAILABILITY: Libraries of S-plus and R routines to implement the method are freely available from www.stjuderesearch.org/depts/biostats.     

RALEE--RNA ALignment editor in Emacs

SUMMARY: Production of high quality multiple sequence alignments of structured RNAs relies on an iterative combination of manual editing and structure prediction. An essential feature of an RNA alignment editor is the facility to mark-up the alignment based on how it matches a given secondary structure prediction, but few available alignment editors offer such a feature. The RALEE (RNA ALignment Editor in Emacs) tool provides a simple environment for RNA multiple sequence alignment editing, including structure-specific colour schemes, utilizing helper applications for structure prediction and many more conventional editing functions. This is accomplished by extending the commonly used text editor, Emacs, which is available for Linux, most UNIX systems, Windows and Mac OS. AVAILABILITY: The ELISP source code for RALEE is freely available from http://www.sanger.ac.uk/Users/sgj/ralee/ along with documentation and examples. CONTACT: sgj@sanger.ac.uk     

ToolShop: prerelease inspections for protein structure prediction servers.

The ToolShop server offers a possibility to compare a protein tertiary structure prediction server with other popular servers before releasing it to the public. The comparison is conducted on a set of 203 proteins and the collected models are compared with over 20 other programs using various assessment procedures. The evaluation lasts circa one week. AVAILABILITY: The ToolShop server is available at http://BioInfo.PL/ToolShop/. The administrator should be contacted to couple the tested server to the evaluation suite. CONTACT: leszek@bioinfo.pl SUPPLEMENTARY INFORMATION: The evaluation procedures are similar to those implemented in the continuous online server evaluation program, LiveBench. Additional information is available from its homepage (http://BioInfo.PL/LiveBench/).     

Enhancing scatterplots with smoothed densities

MOTIVATION: Scatterplots of microarray data generally contain a very large number of dots, making it difficult to get a good impression of their distribution in dense areas. RESULTS: We present a fast and simple algorithm for two-dimensional histogram smoothing, to visually enhance scatterplots. AVAILABILITY: Functions for Matlab and R are available from the corresponding author.     

Statistical adjustment of signal censoring in gene expression experiments

MOTIVATION: Numerical output of spotted microarrays displays censoring of pixel intensities at some software dependent threshold. This reduces the quality of gene expression data, because it seriously violates the linearity of expression with respect to signal intensity. Statistical methods based on typically available spot summaries together with some parametric assumptions can suggest ways to correct for this defect. RESULTS: A maximum likelihood approach is suggested together with a sensible approximation to the joint density of the mean, median and variance-which are typically available to the biological end-user. The method 'corrects' the gene expression values for pixel censoring. A by-product of our approach is a comparison between several two-parameter models for pixel intensity values. It suggests that pixels separated by one or two other pixels can be considered independent draws from a Lognormal or a Gamma distribution. AVAILABILITY: The R/S-Plus code is available at http://www.stats.gla.ac.uk/~microarray/software.