PubMed Health⌕ Search

PubMed · 15201059

3DCoffee: combining protein sequences and structures within multiple sequence alignments.

Abstract

Most bioinformatics analyses require the assembly of a multiple sequence alignment. It has long been suspected that structural information can help to improve the quality of these alignments, yet the effect of combining sequences and structures has not been evaluated systematically. We developed 3DCoffee, a novel method for combining protein sequences and structures in order to generate high-quality multiple sequence alignments. 3DCoffee is based on TCoffee version 2.00, and uses a mixture of pairwise sequence alignments and pairwise structure comparison methods to generate multiple sequence alignments. We benchmarked 3DCoffee using a subset of HOMSTRAD, the collection of reference structural alignments. We found that combining TCoffee with the threading program Fugue makes it possible to improve the accuracy of our HOMSTRAD dataset by four percentage points when using one structure only per dataset. Using two structures yields an improvement of ten percentage points. The measures carried out on HOM39, a HOMSTRAD subset composed of distantly related sequences, show a linear correlation between multiple sequence alignment accuracy and the ratio of number of provided structure to total number of sequences. Our results suggest that in the case of distantly related sequences, a single structure may not be enough for computing an accurate multiple sequence alignment.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Orla O'Sullivan, Karsten Suhre, Chantal Abergel, Desmond G Higgins, Cédric Notredame. 2004-07-02. 3DCoffee: combining protein sequences and structures within multiple sequence alignments.. https://doi.org/10.1016/j.jmb.2004.04.058

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

Crystal structures of glutamine:phenylpyruvate aminotransferase from Thermus thermophilus HB8: induced fit and substrate recognition.

The following three-dimensional structures of three forms of glutamine:phenylpyruvate aminotransferase from Thermus thermophilus HB8 have been determined and represent the first x-ray analysis of the enzyme: the unliganded pyridoxal 5'-phosphate form at 1.9 A resolution and two complexes with 3-phenylpropionate and alpha-keto-gamma-methylthiobutyrate at 2.35 and 2.6 A resolution, respectively. The enzyme shows high activity toward phenylalanine, tyrosine, tryptophan, kynurenine, methionine, and glutamine. The enzyme is a homodimer, and each subunit is divided into an N-terminal arm and small and large domains. Based on its folding, the enzyme belongs to fold type I, aminotransferase subclass Ib. The subclass I aminotransferases whose structures have so far been determined exhibit a large movement of the small domain region upon binding of a substrate. Similarly, the glutamine:phenylpyruvate aminotransferase undergoes a large movement in part of the small domain to close the active site. The active-site pocket has a shape and size suitable to enclose the side chain of an aromatic amino acid or that of methionine. The inner side of the pocket is mostly hydrophobic, but also has polar sites. The kynurenine complex generated by computer modeling fits the pocket of the enzyme and its hydrophilic groups interact with the polar sites of the pocket.

Protein Conformation↗

Identifying residue-residue clashes in protein hybrids by using a second-order mean-field approach.

In this article, a second-order mean-field-based approach is introduced for characterizing the complete set of residue-residue couplings consistent with a given protein structure. This information is subsequently used to classify protein hybrids with respect to their potential to be functional based on the presenceabsence and severity of clashing residue-residue interactions. First, atomistic representations of both the native and denatured states are used to calculate rotamer-backbone, rotamer-intrinsic, and rotamer-rotamer conformational energies. Next, this complete conformational energy table is coupled with a second-order mean-field description to elucidate the probabilities of all possible rotamer-rotamer combinations in a minimum Helmholtz free-energy ensemble. Computational results for the dihydrofolate reductase family reveal correlation in substitution patterns between not only contacting but also distal second-order structural elements. Residue-residue clashes in hybrid proteins are quantified by contrasting the ensemble probabilities of protein hybrids against the ones of the original parental sequences. Good agreement with experimental data is demonstrated by superimposing these clashes against the functional crossover profiles of bidirectional incremental truncation libraries for Escherichia coli and human glycinamide ribonucleotide transformylases.

Protein Conformation↗