KvSNP: accurately predicting the effect of genetic variants in voltage-gated potassium channels.

MOTIVATION: Non-synonymous single nucleotide polymorphisms (nsSNPs) in voltage-gated potassium (Kv) channels cause diseases with potentially fatal consequences in seemingly healthy individuals. Identifying disease-causing genetic variation will aid presymptomatic diagnosis and treatment of such disorders. NsSNP-effect predictors are hypothesized to perform best when developed for specific gene families. We, thus, created KvSNP: a method that assigns a disease-causing probability to Kv-channel nsSNPs. RESULTS: KvSNP outperforms popular non gene-family-specific methods (SNPs&GO, SIFT and Polyphen) in predicting the disease potential of Kv-channel variants, according to all tested metrics (accuracy, Matthews correlation coefficient and area under receiver operator characteristic curve). Most significantly, it increases the separation of the median predicted disease probabilities between benign and disease-causing SNPs by 26% on the next-best competitor. KvSNP has ranked 172 uncharacterized Kv-channel nsSNPs by disease-causing probability. AVAILABILITY AND IMPLEMENTATION: KvSNP, a WEKA implementation is available at www.bioinformatics.leeds.ac.uk/KvDB/KvSNP.html. CONTACT: d.r.westhead@leeds.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Combining the interactome and deleterious SNP predictions to improve disease gene identification.

A method has been developed for the prediction of proteins involved in genetic disorders. This involved combining deleterious SNP prediction with a system based on protein interactions and phenotype distances; this is the first time that deleterious SNP prediction has been used to make predictions across linkage-intervals. At each step we tested and selected the best procedure, revealing that the computationally expensive method of assigning medical meta-terms to create a phenotype distance matrix was outperformed by a simple word counting technique. We carried out in-depth benchmarking with increasingly stringent data sets, reaching precision values of up to 75% (19% recall) for 10-Mb linkage-intervals (averaging 100 genes). For the most stringent (worst-case) data we attained an overall recall of 6%, yet still achieved precision values of up to 90% (4% recall). At all levels of stringency and precision the addition of predicted deleterious SNPs was shown to increase recall.

Protein structure prediction.

Genome sequencing projects continue to provide a flood of new protein sequences, and prediction methods remain an important means of adding structural information. Recently, there have been advances in secondary structure prediction, which feed, in turn, into improved fold recognition algorithms. Finally, there have been technical improvements in comparative modelling, and studies of the expected accuracy of three-dimensional structural models built by this method.

Automated derivation and refinement of sequence length patterns for protein sequences using evolutionary computation.

Several stratagems are used in protein bioinformatics for the classification of proteins based on sequence, structure or function. We explore the concept of a minimal signature embedded in a sequence that defines the likely position of a protein in a classification. Specifically, we address the derivation of sparse profiles for the G-protein coupled receptor (GPCR) clan of integral membrane proteins. We present an evolutionary algorithm (EA) for the derivation of sparse profiles (signatures) without the need to supply a multiple alignment. We also apply an evolution strategy (ES) to the problem of pattern and profile refinement. Patterns were derived for the GPCR 'superfamily' and GPCR families 1-3 individually from starting populations of randomly generated signatures, using a database of integral membrane protein sequences and an objective function using a modified receiver operator characteristic (ROC) statistic. The signature derived for the family 1 GPCR sequences was shown to perform very well in a stringent cross-validation test, detecting 76% of unseen GPCR sequences at 5% error. Application of the ES refinement method to a signature developed by a previously described method [Sadowski, M.I., Parish, J.H., 2003. Automated generation and refinement of protein signatures: case study with G-protein coupled receptors. Bioinformatics 19, 727-734] resulted in a 6% increase of coverage for 5% error as measured in the validation test. We note that there might be a limit to this or any classification of proteins based on patterns or schemata.

Protein structural topology: Automated analysis and diagrammatic representation.

The topology of a protein structure is a highly simplified description of its fold including only the sequence of secondary structure elements, and their relative spatial positions and approximate orientations. This information can be embodied in a two-dimensional diagram of protein topology, called a TOPS cartoon. These cartoons are useful for the understanding of particular folds and making comparisons between folds. Here we describe a new algorithm for the production of TOPS cartoons, which is more robust than those previously available, and has a much higher success rate. This algorithm has been used to produce a database of protein topology cartoons that covers most of the data bank of known protein structures.

PRO_LIGAND: an approach to de novo molecular design. 2. Design of novel molecules from molecular field analysis (MFA) models and pharmacophores.

A computational approach for molecular design, PRO_LIGAND, has been developed within the PROMETHEUS molecular design and simulation system in order to provide a unified framework for the de novo generation of diverse molecules which are either similar or complementary to a specified target. In this instance, the target is a pharmacophore derived from a series of active structures either by a novel interpretation of molecular field analysis data or by a pharmacophore-mapping procedure based on clique detection. After a brief introduction to PRO_LIGAND, a detailed description is given of the two pharmacophore generation procedures and their abilities are demonstrated by the elucidation of pharmacophores for steroid binding and ACE inhibition, respectively. As a further indication of its efficacy in aiding the rational drug design process, PRO_LIGAND is then employed to build novel organic molecules to satisfy the physicochemical constraints implied by the pharmacophores.

A physiologically based mathematical model of dermal absorption in man.

A sound understanding of the mechanisms determining percutaneous absorption is necessary for toxicological risk assessment of chemicals contacting the skin. As part of a programme investigating these mechanisms we have developed a physiologically based mathematical model. The structure of the model parallels the multi-layer structure of the skin, with separate surface, stratum corneum and viable tissue layers. It simulates the effects of partitioning and diffusive transport between the sub-layers, and metabolism in the viable epidermis. In addition the model describes removal processes on the surface of the skin, including the effects of washing and desquamation, and rubbing off onto clothing. This model is applied to data on the penetration of the herbicide fluazifop-butyl through human skin in vivo and in vitro. Part of this dataset is used to estimate unknown model parameter values and the remainder is used to provide a partial validation of the model. Only a small fraction of the applied dose was absorbed through the skin; most of it was removed by washing or onto clothing. The model provides a quantitative description of these loss processes on the skin surface.

Depletion of RUNX1/ETO in t(8;21) AML cells leads to genome-wide changes in chromatin structure and transcription factor binding.

The t(8;21) translocation fuses the DNA-binding domain of the hematopoietic master regulator RUNX1 to the ETO protein. The resultant RUNX1/ETO fusion protein is a leukemia-initiating transcription factor that interferes with RUNX1 function. The result of this interference is a block in differentiation and, finally, the development of acute myeloid leukemia (AML). To obtain insights into RUNX1/ETO-dependant alterations of the epigenetic landscape, we measured genome-wide RUNX1- and RUNX1/ETO-bound regions in t(8;21) cells and assessed to what extent the effects of RUNX1/ETO on the epigenome depend on its continued expression in established leukemic cells. To this end, we determined dynamic alterations of histone acetylation, RNA Polymerase II binding and RUNX1 occupancy in the presence or absence of RUNX1/ETO using a knockdown approach. Combined global assessments of chromatin accessibility and kinetic gene expression data show that RUNX1/ETO controls the expression of important regulators of hematopoietic differentiation and self-renewal. We show that selective removal of RUNX1/ETO leads to a widespread reversal of epigenetic reprogramming and a genome-wide redistribution of RUNX1 binding, resulting in the inhibition of leukemic proliferation and self-renewal, and the induction of differentiation. This demonstrates that RUNX1/ETO represents a pivotal therapeutic target in AML.

Sequence relationships in the legume lectin fold and other jelly rolls.

Distant sequence relationships in proteins containing the beta jelly-roll fold were investigated using sensitive sequence comparison methods, including PSI-BLAST and Hidden Markov Models. A relationship was identified between the rmlC-like and phosphomannose isomerase SCOP (version 1.53) superfamilies, which were merged in the most recent SCOP release. No other distant sequence relationships linking jelly roll superfamilies were found.

A comparative study of machine-learning methods to predict the effects of single nucleotide polymorphisms on protein function.

MOTIVATION: The large volume of single nucleotide polymorphism data now available motivates the development of methods for distinguishing neutral changes from those which have real biological effects. Here, two different machine-learning methods, decision trees and support vector machines (SVMs), are applied for the first time to this problem. In common with most other methods, only non-synonymous changes in protein coding regions of the genome are considered. RESULTS: In detailed cross-validation analysis, both learning methods are shown to compete well with existing methods, and to out-perform them in some key tests. SVMs show better generalization performance, but decision trees have the advantage of generating interpretable rules with robust estimates of prediction confidence. It is shown that the inclusion of protein structure information produces more accurate methods, in agreement with other recent studies, and the effect of using predicted rather than actual structure is evaluated. AVAILABILITY: Software is available on request from the authors.

Evolutionary algorithms in computer-aided molecular design.

In recent years, search and optimisation algorithms inspired by evolutionary processes have been applied with marked success to a wide variety of problems in diverse fields of study. In this review, we survey the growing application of these 'evolutionary algorithms' in one such area: computer-aided molecular design. In the course of the review, we seek to summarise the work to date and to indicate where evolutionary algorithms have met with success and where they have not fared so well. In addition to this, we also attempt to discern some future trends in both the basic research concerning these algorithms and their application to the elucidation, design and modelling of chemical and biochemical structures.

Multiple structural alignment for distantly related all beta structures using TOPS pattern discovery and simulated annealing.

Topsalign is a method that will structurally align diverse protein structures, for example, structural alignment of protein superfolds. All proteins within a superfold share the same fold but often have very low sequence identity and different biological and biochemical functions. There is often significant structural diversity around the common scaffold of secondary structure elements of the fold. Topsalign uses topological descriptions of proteins. A pattern discovery algorithm identifies equivalent secondary structure elements between a set of proteins and these are used to produce an initial multiple structure alignment. Simulated annealing is used to optimize the alignment. The output of Topsalign is a multiple structure-based sequence alignment and a 3D superposition of the structures. This method has been tested on three superfolds: the beta jelly roll, TIM (alpha/beta) barrel and the OB fold. Topsalign outperforms established methods on very diverse structures. Despite the pattern discovery working only on beta strand secondary structure elements, Topsalign is shown to align TIM (alpha/beta) barrel superfamilies, which contain both alpha helices and beta strands.

Petri Net representations in systems biology.

The mathematical structures known as Petri Nets have recently become the focus of much research effort in both the structural and quantitative analysis of all kinds of biological networks. This review provides a very brief summary of these interesting new research directions.

A comparison of heuristic search algorithms for molecular docking.

This paper describes the implementation and comparison of four heuristic search algorithms (genetic algorithm, evolutionary programming, simulated annealing and tabu search) and a random search procedure for flexible molecular docking. To our knowledge, this is the first application of the tabu search algorithm in this area. The algorithms are compared using a recently described fast molecular recognition potential function and a diverse set of five protein-ligand systems. Statistical analysis of the results indicates that overall the genetic algorithm performs best in terms of the median energy of the solutions located. However, tabu search shows a better performance in terms of locating solutions close to the crystallographic ligand conformation. These results suggest that a hybrid search algorithm may give superior results to any of the algorithms alone.

AI-based algorithms for protein surface comparisons.

Many current methods for protein analysis depend on the detection of similarity in either the primary sequence, or the overall tertiary structure (the Calpha atoms of the protein backbone). These common sequences or structures may imply similar functional characteristics or active properties. Active sites and ligand binding sites usually occur on or near the surface of the protein; so similarly shaped surface regions could imply similar functions. We investigate various methods for describing the shape properties of protein surfaces and for comparing them. Our current work uses algorithms from computer vision to describe the protein surfaces, and methods from graph theory to compare the surface regions. Early results indicate that we can successfully match a family of related ligand binding sites, and find their similarly shaped surface regions. This method of surface analysis could be extended to help identify unknown surface regions for possible ligand binding or active sites.

Active-site-directed 3D database searching: pharmacophore extraction and validation of hits.

Two new computational tools, PRO_PHARMEX and PRO_SCOPE, for use in active-site-directed searching of 3D databases are described. PRO_PHARMEX is a flexible, graphics-based program facilitating the extraction of pharmacophores from the active site of a target macromolecule. These pharmacophores can then be used to search a variety of databases for novel lead compounds. Such searches can often generate many 'hits' of varying quality. To aid the user in setting priorities for purchase, synthesis or testing, PRO_SCOPE can be used to dock molecules rapidly back into the active site and to assign them a score using an empirical scoring function correlated to the free energy of binding. To illustrate how these tools can add value to existing 3D database software, their use in the design of novel thrombin inhibitors is described.

PRO-LIGAND: an approach to de novo molecular design. 3. A genetic algorithm for structure refinement.

Recently, the development of computer programs which permit the de novo design of molecular structures satisfying a set of steric and chemical constraints has become a burgeoning area of research and many operational systems have been reported in the literature. Experience with PRO-LIGAND-the de novo design methodology embodied in our in-house molecular design and simulation system PRO-METHEUS-has suggested that the addition of a genetic algorithm (GA) structure refinement procedure can 'add value' to an already useful tool. Starting with the set of designed molecules as an initial population, the GA can combine features from both high- and low-scoring structures and, over a number of generations, produce individuals of better score than any of the starting structures. This paper describes how we have implemented such a procedure and demonstrates its efficacy in improving two sets of molecules generated by different de novo design projects.

Flexible docking using Tabu search and an empirical estimate of binding affinity.

This article describes the implementation of a new docking approach. The method uses a Tabu search methodology to dock flexibly ligand molecules into rigid receptor structures. It uses an empirical objective function with a small number of physically based terms derived from fitting experimental binding affinities for crystallographic complexes. This means that docking energies produced by the searching algorithm provide direct estimates of the binding affinities of the ligands. The method has been tested on 50 ligand-receptor complexes for which the experimental binding affinity and binding geometry are known. All water molecules are removed from the structures and ligand molecules are minimized in vacuo before docking. The lowest energy geometry produced by the docking protocol is within 1.5 A root-mean square of the experimental binding mode for 86% of the complexes. The lowest energies produced by the docking are in fair agreement with the known free energies of binding for the ligands.

New tools and resources for analysing protein structures and their interactions.

The determination of protein structures has furthered our understanding of how various proteins perform their functions. With the large number of structures currently available in the PDB, it is necessary to be able to easily study these proteins in detail. Here new software tools are presented which aim to facilitate this analysis; these include the PDBsum WWW site which provides a summary description of all PDB entries, the programs TOPS and NUCPLOT to plot schematic diagrams representing protein topology and DNA-binding interactions, SAS a WWW-based sequence-analysis tool incorporating structural data, and WWW servers for the analysis of protein-protein interfaces and analyses of over 300 haem-binding proteins.