Simon Penel

Tree pattern matching in phylogenetic trees: automatic search for orthologs or paralogs in homologous gene sequence databases

MOTIVATION Comparative sequence analysis is widely used to study genome function and evolution. This approach first requires the identification of homologous genes and then the interpretation of their homology relationships (orthology or paralogy). To provide help in this complex task, we developed three databases of homologous genes containing sequences, multiple alignments and phylogenetic trees: HOBACGEN, HOVERGEN and HOGENOM. In this paper, we present two new tools for automating the search for orthologs or paralogs in these databases. RESULTS First, we have developed and implemented an algorithm to infer speciation and duplication events by comparison of gene and species trees (tree reconciliation). Second, we have developed a general method to search in our databases the gene families for which the tree topology matches a peculiar tree pattern. This algorithm of unordered tree pattern matching has been implemented in the FamFetch graphical interface. With the help of a graphical editor, the user can specify the topology of the tree pattern, and set constraints on its nodes and leaves. Then, this pattern is compared with all the phylogenetic trees of the database, to retrieve the families in which one or several occurrences of this pattern are found. By specifying ad hoc patterns, it is therefore possible to identify orthologs in our databases.

Integr8 and Genome Reviews: integrated views of complete genomes and proteomes

Integr8 is a new web portal for exploring the biology of organisms with completely deciphered genomes. For over 190 species, Integr8 provides access to general information, recent publications, and a detailed statistical overview of the genome and proteome of the organism. The preparation of this analysis is supported through Genome Reviews, a new database of bacterial and archaeal DNA sequences in which annotation has been upgraded (compared to the original submission) through the integration of data from many sources, including the EMBL Nucleotide Sequence Database, the UniProt Knowledgebase, InterPro, CluSTr, GOA and HOGENOM. Integr8 also allows the users to customize their own interactive analysis, and to download both customized and prepared datasets for their own use. Integr8 is available at http://www.ebi.ac.uk/integr8.

Databases of homologous gene families for comparative genomics

BackgroundComparative genomics is a central step in many sequence analysis studies, from gene annotation and the identification of new functional regions in genomes, to the study of evolutionary processes at the molecular level (speciation, single gene or whole genome duplications, etc.) and phylogenetics. In that context, databases providing users high quality homologous families and sequence alignments as well as phylogenetic trees based on state of the art algorithms are becoming indispensable.MethodsWe developed an automated procedure allowing massive all-against-all similarity searches, gene clustering, multiple alignments computation, and phylogenetic trees construction and reconciliation. The application of this procedure to a very large set of sequences is possible through parallel computing on a large computer cluster.ResultsThree databases were developed using this procedure: HOVERGEN, HOGENOM and HOMOLENS. These databases share the same architecture but differ in their content. HOVERGEN contains sequences from vertebrates, HOGENOM is mainly devoted to completely sequenced microbial organisms, and HOMOLENS is devoted to metazoan genomes from Ensembl. Access to the databases is provided through Web query forms, a general retrieval system and a client-server graphical interface. The later can be used to perform tree-pattern based searches allowing, among other uses, to retrieve sets of orthologous genes. The three databases, as well as the software required to build and query them, can be used or downloaded from the PBIL (Pôle Bioinformatique Lyonnais) site at http://pbil.univ-lyon1.fr/.

Ultra-fast sequence clustering from similarity networks with SiLiX

BackgroundThe number of gene sequences that are available for comparative genomics approaches is increasing extremely quickly. A current challenge is to be able to handle this huge amount of sequences in order to build families of homologous sequences in a reasonable time.ResultsWe present the software package SiLiX that implements a novel method which reconsiders single linkage clustering with a graph theoretical approach. A parallel version of the algorithms is also presented. As a demonstration of the ability of our software, we clustered more than 3 millions sequences from about 2 billion BLAST hits in 7 minutes, with a high clustering quality, both in terms of sensitivity and specificity.ConclusionsComparing state-of-the-art software, SiLiX presents the best up-to-date capabilities to face the problem of clustering large collections of sequences. SiLiX is freely available at http://lbbe.univ-lyon1.fr/SiLiX.

Neutron crystallographic evidence of lipase-colipase complex activation by a micelle

The concept of lipase interfacial activation stems from the finding that the catalytic activity of most lipases depends on the aggregation state of their substrates. It is thought that activation involves the unmasking and structuring of the enzymes active site through conformational changes requiring the presence of oil‐in‐water droplets. Here, we present the neutron structure of the activated lipase–colipase–micelle complex as determined using the D2O/H2O contrast variation low resolution diffraction method. In the ternary complex, the disk‐shaped micelle interacts extensively with the concave face of colipase and the distal tip of the C‐terminal domain of lipase. Since the micelle‐ and substrate‐binding sites concern different regions of the protein complex, we conclude that lipase activation is not interfacial but occurs in the aqueous phase and is mediated by colipase and a micelle.

Pervasive positive selection on duplicated and nonduplicated vertebrate protein coding genes

A stringent branch-site codon model was used to detect positive selection in vertebrate evolution. We show that the test is robust to the large evolutionary distances involved. Positive selection was detected in 77% of 884 genes studied. Most positive selection concerns a few sites on a single branch of the phylogenetic tree: Between 0.9% and 4.7% of sites are affected by positive selection depending on the branches. No functional category was overrepresented among genes under positive selection. Surprisingly, whole genome duplication had no effect on the prevalence of positive selection, whether the fish-specific genome duplication or the two rounds at the origin of vertebrates. Thus positive selection has not been limited to a few gene classes, or to specific evolutionary events such as duplication, but has been pervasive during vertebrate evolution.

YodA from Escherichia coli is a Metal-binding, Lipocalin-like Protein

We have determined the crystal structure of YodA, an Escherichia coli protein of unknown function. YodA had been identified under conditions of cadmium stress, and we confirm that it binds metals such as cadmium and zinc. We have also found nickel bound in one of the crystal forms. YodA is composed of two domains: a main lipocalin/calycin-like domain and a helical domain. The principal metal-binding site lies on one side of the calycin domain, thus making YodA the first metal-binding lipocalin known. Our experiments suggest that YodA expression may be part of a more general stress response. From sequence analogy with the C-terminal domain of a metal-binding receptor of a member of bacterial ATP-binding cassette transporters, we propose a three-dimensional model for this receptor and suggest that YodA may have a receptor-type partner in E. coli.

Stabilizing nonpolar/polar side-chain interactions in the alpha-helix.

A simplistic, yet often used, view of protein stability is that amino acids attract other amino acids with similar polarity, whereas nonpolar and polar side chains repel. Here we show that nonpolar/polar interactions, namely Val or Ile bonding to Lys or Arg in α‐helices, can in fact be stabilizing. Residues spaced i, i + 4 in α‐helices are on the same face of the helix, with potential to favorably interact and stabilize the structure. We observe that the nonpolar/polar pairs Ile‐Lys, Ile‐Arg, and Val‐Lys occur in protein helices more often than expected when spaced i, i + 4. Partially helical peptides containing pairs of nonpolar/polar residues were synthesized. Controls with i, i + 5 spacing have the residues on opposite faces of the helix and are less helical than the test peptides with the i, i + 4 interactions. Experimental circular dichroism results were analyzed with helix‐coil theory to calculate the free energy for the interactions. All three stabilize the helix with ΔG between −0.14 and −0.32 kcal · mol−1. The interactions are hydrophobic with contacts between Val or Ile and the alkyl groups in Arg or Lys. Side chains such as Lys and Arg can thus interact favorably with both polar and nonpolar residues. Proteins 2001;45:449–455.

Integrated databanks access and sequence/structure analysis services at the PBIL

The World Wide Web server of the PBIL (Pôle Bioinformatique Lyonnais) provides on-line access to sequence databanks and to many tools of nucleic acid and protein sequence analyses. This server allows to query nucleotide sequence banks in the EMBL and GenBank formats and protein sequence banks in the SWISS-PROT and PIR formats. The query engine on which our data bank access is based is the ACNUC system. It allows the possibility to build complex queries to access functional zones of biological interest and to retrieve large sequence sets. Of special interest are the unique features provided by this system to query the data banks of gene families developed at the PBIL. The server also provides access to a wide range of sequence analysis methods: similarity search programs, multiple alignments, protein structure prediction and multivariate statistics. An originality of this server is the integration of these two aspects: sequence retrieval and sequence analysis. Indeed, thanks to the introduction of re-usable lists, it is possible to perform treatments on large sets of data. The PBIL server can be reached at: http://pbil.univ-lyon1.fr.

Detergent binding in trigonal crystals of OmpF porin from Escherichia coli

The structure of the detergent, ocytyl hydroxyethylsufoxide (C8(HE)SO), bound to the OmpF porin from E coli (in the trigonal crystal form) has been determined by neutron crystallography. Due to a dynamic exchange of detergent molecules with their environment they are not ordered on an atomic scale. The structure reported here is therefore at a resolution of approximately 16 A. The X-ray crystallographically determined structure of the protein provides a starting point for the neutron analysis in which the detergent is visualized primarily thanks to its high contrast against D2O. The structure shows the detergent to be located mainly in two areas. It forms toroidal annuli around each OmpF trimer, these annuli fusing to form a detergent belt surrounding a solvent filled column traversing the crystal. Those areas of the protein to which the detergent binds are formed almost exclusively of hydrophobic residues and form a band about 30 A high around the trimer. Its upper and lower bounds are defined by two bands of aromatic residues, tyrosines pointing away from the detergent belt and interacting with the polar headgroups while phenylalanines point inwards. This strongly suggests that the same areas define, in vivo, the location at which protein interacts with lipid. The hydrophobic moiety of detergent is also found mediating the hydrophobic protein-protein interactions at the interface between two trimers on the crystallographic two-fold axis.