Search PubMed⌕ Search
Home/Archive/PubMed

PubMed · 11734009

Towards a reliable objective function for multiple sequence alignments.

J D Thompson·F Plewniak·R Ripp·J C Thierry·O Poch
View original source ↗Export metadata ↓

Abstract

Multiple sequence alignment is a fundamental tool in a number of different domains in modern molecular biology, including functional and evolutionary studies of a protein family. Multiple alignments also play an essential role in the new integrated systems for genome annotation and analysis. Thus, the development of new multiple alignment scores and statistics is essential, in the spirit of the work dedicated to the evaluation of pairwise sequence alignments for database searching techniques. We present here norMD, a new objective scoring function for multiple sequence alignments. NorMD combines the advantages of the column-scoring techniques with the sensitivity of methods incorporating residue similarity scores. In addition, norMD incorporates ab initio sequence information, such as the number, length and similarity of the sequences to be aligned. The sensitivity and reliability of the norMD objective function is demonstrated using structural alignments in the SCOP and BAliBASE databases. The norMD scores are then applied to the multiple alignments of the complete sequences (MACS) detected by BlastP with E-value<10, for a set of 734 hypothetical proteins encoded by the Vibrio cholerae genome. Unrelated or badly aligned sequences were automatically removed from the MACS, leaving a high-quality multiple alignment which could be reliably exploited in a subsequent functional and/or structural annotation process. After removal of unreliable sequences, 176 (24 %) of the alignments contained at least one sequence with a functional annotation. 103 of these new matches were supported by significant hits to the Interpro domain and motif database.

Explore related subjects

Amino Acid MotifsAmino Acid SequenceArchaeal ProteinsBacterial ProteinsComputational BiologyDatabases, GeneticEukaryotic CellsGenome, BacterialGenomicsModels, MolecularMolecular Sequence DataProtein Structure, TertiaryReproducibility of ResultsResearch DesignSensitivity and SpecificitySequence AlignmentSoftwareVibrio choleraeJournal ArticleResearch Support, Non-U.S. Gov't

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

J D Thompson, F Plewniak, R Ripp, J C Thierry, O Poch. 2001-12-07. Towards a reliable objective function for multiple sequence alignments.. https://doi.org/10.1006/jmbi.2001.5187

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

KEEP EXPLORING

Related citations

PubMed2026

GAMMA: gap-aware motif mining under incomplete labeling with applications to MHC motifs.

MOTIVATION: Sequence motif identification is crucial for understanding molecular recognition, particularly in immune responses involving peptide binding to major histocompatibility complex (MHC) Class I molecules for antigen presentation to T cells. Traditionally, MHC Class I binding motifs are assumed to be contiguous and span nine amino acids. However, structural evidence suggests that binding may involve nonadjacent residues, challenging the assumptions of existing methods. RESULTS: In this study, we propose Gap-Aware Motif Mining Algorithm (GAMMA), a probabilistic framework designed to identify noncontiguous motifs under conditions of incomplete labeling. GAMMA employs Bayesian inference with Markov chain Monte Carlo sampling to jointly estimate motif parameters, binding locations, and the relative spacing between binding positions. Through extensive simulations and real-world applications to MHC Class I peptide datasets, GAMMA outperforms existing motif discovery tools such as GLAM2 in accurately localizing binding residues and identifying the underlying motifs. Notably, our results suggest that the true number of binding residues may be eight, fewer than the commonly assumed nine. In addition, for longer peptides, the model captures increased flexibility in the central region, consistent with structural observations that peptides may bulge in the middle. AVAILABILITY AND IMPLEMENTATION: The raw data and the source codes are available on GitHub (https://github.com/RanLIUaca/GAMMAmotif).

Amino Acid Motifs↗
PubMed2015

Signalling thresholds and negative B-cell selection in acute lymphoblastic leukaemia.

B cells are selected for an intermediate level of B-cell antigen receptor (BCR) signalling strength: attenuation below minimum (for example, non-functional BCR) or hyperactivation above maximum (for example, self-reactive BCR) thresholds of signalling strength causes negative selection. In &#x223c;25% of cases, acute lymphoblastic leukaemia (ALL) cells carry the oncogenic BCR-ABL1 tyrosine kinase (Philadelphia chromosome positive), which mimics constitutively active pre-BCR signalling. Current therapeutic approaches are largely focused on the development of more potent tyrosine kinase inhibitors to suppress oncogenic signalling below a minimum threshold for survival. We tested the hypothesis that targeted hyperactivation--above a maximum threshold--will engage a deletional checkpoint for removal of self-reactive B cells and selectively kill ALL cells. Here we find, by testing various components of proximal pre-BCR signalling in mouse BCR-ABL1 cells, that an incremental increase of Syk tyrosine kinase activity was required and sufficient to induce cell death. Hyperactive Syk was functionally equivalent to acute activation of a self-reactive BCR on ALL cells. Despite oncogenic transformation, this basic mechanism of negative selection was still functional in ALL cells. Unlike normal pre-B cells, patient-derived ALL cells express the inhibitory receptors PECAM1, CD300A and LAIR1 at high levels. Genetic studies revealed that Pecam1, Cd300a and Lair1 are critical to calibrate oncogenic signalling strength through recruitment of the inhibitory phosphatases Ptpn6 (ref. 7) and Inpp5d (ref. 8). Using a novel small-molecule inhibitor of INPP5D (also known as SHIP1), we demonstrated that pharmacological hyperactivation of SYK and engagement of negative B-cell selection represents a promising new strategy to overcome drug resistance in human ALL.

Amino Acid Motifs↗
PubMed2002

Selection of peptide entry motifs by bacterial surface display.

Surface display technologies have been established previously to select peptides and polypeptides that interact with purified immobilized ligands. In the present study, we designed and implemented a surface display-based technique to identify novel peptide motifs that mediate entry into eukaryotic cells. An Escherichia coli library expressing surface-displayed peptides was combined with eukaryotic cells and the gentamicin protection assay was performed to select recombinant E. coli, which were internalized into eukaryotic cells by virtue of the displayed peptides. To establish the proof of principle of this approach, the fibronectin-binding motifs of the fibronectin-binding protein A of Staphylococcus aureus were inserted into the E. coli FhuA protein. Surface expression of the fusion proteins was demonstrated by functional assays and by FACS analysis. The fibronectin-binding motifs were shown to mediate entry of the bacteria into non-phagocytic eukaryotic cells and brought about the preferential selection of these bacteria over E. coli expressing parental FhuA, with an enrichment of 100000-fold. Four entry sequences were selected and identified using an S. aureus library of peptides displayed in the FhuA protein on the surface of E. coli. These sequences included novel entry motifs as well as integrin-binding Arg-Gly-Asp (RGD) motifs and promoted a high degree of bacterial entry. Bacterial surface display is thus a powerful tool to effectively select and identify entry peptide motifs.

Amino Acid Motifs↗