Raphaël Leplae

Assessment of CAPRI predictions in rounds 3–5 shows progress in docking procedures

The current status of docking procedures for predicting protein–protein interactions starting from their three‐dimensional (3D) structure is reassessed by evaluating blind predictions, performed during 2003–2004 as part of Rounds 3–5 of the community‐wide experiment on Critical Assessment of PRedicted Interactions (CAPRI). Ten newly determined structures of protein–protein complexes were used as targets for these rounds. They comprised 2 enzyme–inhibitor complexes, 2 antigen–antibody complexes, 2 complexes involved in cellular signaling, 2 homo‐oligomers, and a complex between 2 components of the bacterial cellulosome. For most targets, the predictors were given the experimental structures of 1 unbound and 1 bound component, with the latter in a random orientation. For some, the structure of the free component was derived from that of a related protein, requiring the use of homology modeling. In some of the targets, significant differences in conformation were displayed between the bound and unbound components, representing a major challenge for the docking procedures. For 1 target, predictions could not go to completion. In total, 1866 predictions submitted by 30 groups were evaluated. Over one‐third of these groups applied completely novel docking algorithms and scoring functions, with several of them specifically addressing the challenge of dealing with side‐chain and backbone flexibility. The quality of the predicted interactions was evaluated by comparison to the experimental structures of the targets, made available for the evaluation, using the well‐agreed‐upon criteria used previously. Twenty‐four groups, which for the first time included an automatic Web server, produced predictions ranking from acceptable to highly accurate for all targets, including those where the structures of the bound and unbound forms differed substantially. These results and a brief survey of the methods used by participants of CAPRI Rounds 3–5 suggest that genuine progress in the performance of docking methods is being achieved, with CAPRI acting as the catalyst. Proteins 2005;60:150–169.

Assessment of blind predictions of protein–protein interactions: Current status of docking methods

The current status of docking procedures for predicting protein–protein interactions starting from their three‐dimensional structure is assessed from a first major evaluation of blind predictions. This evaluation was performed as part of a communitywide experiment on Critical Assessment of PRedicted Interactions (CAPRI). Seven newly determined structures of protein–protein complexes were available as targets for this experiment. These were the complexes between a kinase and its protein substrate, between a T‐cell receptor β‐chain and a superantigen, and five antigen‐antibody complexes. For each target, the predictors were given the experimental structures of the free components, or of one free and one bound component in a random orientation. The structure of the complex was revealed only at the time of the evaluation. A total of 465 predictions submitted by 19 groups were evaluated. These groups used a wide range of algorithms and scoring functions, some of which were completely novel. The quality of the predicted interactions was evaluated by comparing residue–residue contacts and interface residues to those in the X‐ray structures and by analyzing the fit of the ligand molecules (the smaller of the two proteins in the complex) or of interface residues only, in the predicted versus target complexes. A total of 14 groups produced predictions, ranking from acceptable to highly accurate for five of the seven targets. The use of available biochemical and biological information, and in one instance structural information, played a key role in achieving this result. It was essential for identifying the native binding modes for the five correctly predicted targets, including the kinase‐substrate complex where the enzyme changes conformation on association. But it was also the cause for missing the correct solution for the two remaining unpredicted targets, which involve unexpected antigen‐antibody binding modes. Overall, this analysis reveals genuine progress in docking procedures but also illustrates the remaining serious limitations and points out the need for better scoring functions and more effective ways for handling conformational flexibility. Proteins 2003;52:51–67.

Diversity of bacterial type II toxin–antitoxin systems: a comprehensive search and functional analysis of novel families

Type II toxin–antitoxin (TA) systems are generally composed of two genes organized in an operon, encoding a labile antitoxin and a stable toxin. They were first discovered on plasmids where they contribute to plasmid stability by a phenomenon denoted as ‘addiction’, and subsequently in bacterial chromosomes. To discover novel families of antitoxins and toxins, we developed a bioinformatics approach based on the ‘guilt by association’ principle. Extensive experimental validation in Escherichia coli of predicted antitoxins and toxins increased significantly the number of validated systems and defined novel toxin and antitoxin families. Our data suggest that toxin families as well as antitoxin families originate from distinct ancestors that were assembled multiple times during evolution. Toxin and antitoxin families found on plasmids tend to be promiscuous and widespread, indicating that TA systems move through horizontal gene transfer. We propose that due to their addictive properties, TA systems are likely to be maintained in chromosomes even though they do not necessarily confer an advantage to their bacterial hosts. Therefore, addiction might play a major role in the evolutionary success of TA systems both on mobile genetic elements and in bacterial chromosomes.

ACLAME: A CLAssification of Mobile genetic Elements

The ACLAME database (http://aclame.ulb.ac.be) is a collection and classification of prokaryotic mobile genetic elements (MGEs) from various sources, comprising all known phage genomes, plasmids and transposons. In addition to providing information on the full genomes and genetic entities, it aims to build a comprehensive classification of the functional modules of MGEs at the protein, gene and higher levels. This first version contains a comprehensive classification of 5069 proteins from 119 DNA bacteriophages into over 400 functional families. This classification was produced automatically using TRIBE-MCL, a graph-theory-based Markov clustering algorithm that uses sequence measures as input, and then manually curated. Manual curation was aided by consulting annotations available in public databases retrieved through additional sequence similarity searches using Psi-Blast and Hidden Markov Models. The database is publicly accessible and open to expert volunteers willing to participate in its curation. Its web interface allows browsing as well as querying the classification. The main objectives are to collect and organize in a rational way the complexity inherent to MGEs, to extend and improve the inadequate annotation currently associated with MGEs and to screen known genomes for the validation and discovery of new MGEs.

ACLAME: A CLAssification of Mobile genetic Elements, update 2010

The ACLAME database is dedicated to the collection, analysis and classification of sequenced mobile genetic elements (MGEs, in particular phages and plasmids). In addition to providing information on the MGEs content, classifications are available at various levels of organization. At the gene/protein level, families group similar sequences that are expected to share the same function. Families of four or more proteins are manually assigned with a functional annotation using the GeneOntology and the locally developed ontology MeGO dedicated to MGEs. At the genome level, evolutionary cohesive modules group sets of protein families shared among MGEs. At the population level, networks display the reticulate evolutionary relationships among MGEs. To increase the coverage of the phage sequence space, ACLAME version 0.4 incorporates 760 high-quality predicted prophages selected from the Prophinder database. Most of the data can be downloaded from the freely accessible ACLAME web site (http://aclame.ulb.ac.be). The BLAST interface for querying the database has been extended and numerous tools for in-depth analysis of the results have been added.

Prophinder: a computational tool for prophage prediction in prokaryotic genomes

UNLABELLED Prophinder is a prophage prediction tool coupled with a prediction database, a web server and web service. Predicted prophages will help to fill the gaps in the current sparse phage sequence space, which should cover an estimated 100 million species. Systematic and reliable predictions will enable further studies of prophages contribution to the bacteriophage gene pool and to better understand gene shuffling between prophages and phages infecting the same host. AVAILABILITY Softare is available at http://aclame.ulb.ac.be/prophinder

Evidence for copper homeostasis function of metallothionein (MT3) in the hyperaccumulator Thlaspi caerulescens

Metallothioneins chelate metals and consequently may be a control point of metal homeostasis. Homologous to type 3 metallothioneins, TcMT3 cDNA was identified in the Cd/Zn hyperaccumulator, Thlaspi caerulescens. TcMT3 amino acid sequence showed modifications in the Cys positions when compared with its Arabidopsis orthologue. A structural model established that the MT3 carboxyterminal domain is similar to the β domain of animal metallothioneins and predicts a smaller cavity to chelate metals for A. thaliana than for T. caerulescens. Functional testing in yeast and Northern blot analysis added further evidence for adaptative variations of MT3 for the maintenance of Cu homeostasis in a metal hyperaccumulator.

Variations in plant metallothioneins: the heavy metal hyperaccumulator Thlaspi caerulescens as a study case

Plant metallothioneins (MTs) are extremely diverse and are thought to be involved in metal homeostasis or detoxification. Thlaspi caerulescens is a model Zn/Cd hyperaccumulator and thus constitutes an ideal system to study the variability of these MTs. Two T. caerulescens cDNAs (accession: 665511; accession: 665515), that are highly homologous to type 1 and type 2 Arabidopsis thaliana MTs, have been isolated using a functional screen for plant cDNAs that confer Cd tolerance to yeast. However, TcMT1 has a much shorter N-terminal domain than that of A. thaliana and so lacks Cys motifs conserved through all the plant MTs classified as type 1. A systematic search in plant databases allowed the detection of MT-related sequences. Sixty-four percent fulfil the criteria for MT classification described in Cobbett and Goldsbrough (2002) and further extend our knowledge about other conserved residues that might play an important role in plant MT structure. In addition, 34% of the total MT-related sequences cannot be classified strictly as they display modifications in the conserved residues according to the current plant MTs’ classification. The significance of this variability in plant MT sequences is discussed. Functional complementation in yeast was used to assess whether these variations may alter the MTs’ function in T. caerulescens. Regulation of the expression of MTs in T. caerulescens was also investigated. TcMT1 and TcMT2 display higher expression in T. caerulescens than in A. thaliana. Moreover, their differential expression patterns in organs and in response to metal exposure, suggest that the two types of MTs may have diverse roles and functions in T. caerulescens.

The Biphenyl- and 4-Chlorobiphenyl-Catabolic Transposon Tn4371, a Member of a New Family of Genomic Islands Related to IncP and Ti Plasmids

ABSTRACT The nucleotide sequence of the biphenyl catabolic transposon Tn4371 has been completed and analyzed. It confirmed that the element has a mosaic structure made of several building blocks. In addition to previously identified genes coding for a tyrosine recombinase related to phage integrases and for biphenyl degradation enzymes very similar to those of Achromobacter georgiopolitanum KKS102, Tn4371 carries many plasmid-related genes involved in replication, partition, and other, as-yet-unknown, plasmid functions. One gene cluster contains most of the genes required to express a type IV secretion-mating pair formation apparatus coupled with a TraG ATPase, all of which are related to those found on IncP and Ti plasmids. Orthologues of all Tn4371 plasmid-related genes and of the tyrosine recombinase gene were found, with a very similar organization, in the chromosome of Ralstonia solanacearum and on the yet-to-be-determined genomic sequences of Erwinia chrysanthemi and Azotobacter vinelandii. In each of these chromosomal segments, conserved segments were separated by different groups of genes, which also differed from the Tn4371 bph genes. The conserved blocks of genes were also identified, in at least two copies, in the chromosome of Ralstonia metallidurans CH34. Tn4371 thus appears to represent a new family of potentially mobile genomic islands with a broad host range since they reside in a wide range of soil proteobacteria, including plant pathogens.

Genetic diversity of Group A Streptococcus M protein: implications for typing and vaccine development

Group A Streptococci (GAS) are classified into 180 emm-types based on the sequence of the amino-terminal hyper-variable region of the M surface protein. The genetic relatedness of the whole surface-exposed part of M was investigated in well-characterized Belgian and Brazilian GAS isolates which belong to different epidemiological and clinical landscapes. Despite a small number of different emm-types and an apparent low diversity in the Belgian isolates (as revealed by the emm-typing method), our data showed that the overall genetic diversity of the M proteins was higher than expected. On the contrary, the M proteins from the Brazilian isolates were genetically highly related. Since M is a multi-functional protein, an analysis of the whole surface-exposed sequence that takes into account the different functional domains may provide tools for typing as well as for analyzing the molecular mechanisms of M virulence or defining vaccine strategies.