Search PubMed⌕ 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

The effect of box shape on the dynamic properties of proteins simulated under periodic boundary conditions.

The effect of the box shape on the dynamic behavior of proteins simulated under periodic boundary conditions is evaluated. In particular, the influence of simulation boxes defined by the near-densest lattice packing (NDLP) in conjunction with rotational constraints is compared to that of standard box types without these constraints. Three different proteins of varying size, shape, and secondary structure content were examined in the study. The statistical significance of differences in RMSD, radius of gyration, solvent-accessible surface, number of hydrogen bonds, and secondary structure content between proteins, box types, and the application or not of rotational constraints has been assessed. Furthermore, the differences in the collective modes for each protein between different boxes and the application or not of rotational constraints have been examined. In total 105 simulations were performed, and the results compared using a three-way multivariate analysis of variance (MANOVA) for properties derived from the trajectories and a three-way univariate analysis of variance (ANOVA) for collective modes. It is shown that application of roto-translational constraints does not have a statistically significant effect on the results obtained from the different simulations. However, the choice of simulation box was found to have a small (5-10%), but statistically significant effect on the behavior of two of the three proteins included in the study.

Protein Conformation↗

Defining long-range order and local disorder in native alpha-synuclein using residual dipolar couplings.

Intrinsically unstructured proteins play key biochemical roles in a vast range of normal and pathological processes. To study these systems, it is necessary to invoke an ensemble of rapidly interconverting conformations. Residual dipolar couplings (RDCs) are particularly powerful probes of the behavior of unfolded proteins, reporting on time and ensemble-averaged conformations up to and beyond the millisecond time scale. In this study, we present a novel interpretation of RDCs in unfolded systems that simultaneously defines long-range structural order and local conformational sampling. This approach is used to describe the structure and dynamics of alpha-Synuclein (alphaS), a protein that is strongly implicated in the development of Parkinson's disease (PD), allowing unambiguous detection of strongly populated conformers containing long-range contacts between the N- and C-terminal domains. The structural model combines two features required for the description of alphaS in solution: local conformational fluctuation based on random sampling of residue-specific phi/psi distributions, and long-range contacts induced by the presence of nonbonding interactions between domains that are distant in primary sequence. Both aspects are found to be necessary for the reproduction of the nonaveraged RDCs from alphaS. Although RDCs have previously been shown to report on local conformational preferences in unstructured proteins, this study shows the additional sensitivity of these measurements to the presence of long-range order in highly flexible systems.

Protein Conformation↗

Molecular origin of time-dependent fluorescence shifts in proteins.

Time-resolved fluorescence spectroscopy is used increasingly to probe molecular motions at the aqueous interfaces of biological macromolecules and membranes. By recording the time variation of the fluorescence frequency, thermal atomic fluctuations in the vicinity of the chromophore can be probed. From such fluorescence Stokes shift (FSS) experiments, it has been inferred that water motions in the hydration layer are slowed down by 1-3 orders of magnitude. To provide a more secure foundation for the interpretation of FSS data, we use molecular dynamics simulations to examine the molecular origin of the FSS from a tryptophan residue in a protein. By using linear response theory to decompose the FSS into its water and protein components, we find that the water component dominates the static FSS but decays rapidly. Thus, after a few picoseconds, the FSS essentially reflects protein dynamics, including the self-motion of the chromophore. Because of its collective nature, the FSS response is insensitive to the motion of individual water molecules. Collective water displacement by slowly fluctuating protein groups introduces a long-time tail in the water autocorrelation function, but this dynamic coupling is hardly manifested in the observed FSS. Our analysis reconciles FSS data with the picture of a highly dynamic hydration layer, derived mainly from magnetic relaxation dispersion and simulation studies, and calls for a revision of previous interpretations of FSS decays in terms of slow hydration dynamics at biomolecular and other interfaces.

Protein Conformation↗