SEARCH · PubMed Health

Results for “Sequence Alignment”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,045 records · Page 58Linked to original sources

PubMed2006

Genome BLAST distance phylogenies inferred from whole plastid and whole mitochondrion genome sequences.

BACKGROUND: Phylogenetic methods which do not rely on multiple sequence alignments are important tools in inferring trees directly from completely sequenced genomes. Here, we extend the recently described Genome BLAST Distance Phylogeny (GBDP) strategy to compute phylogenetic trees from all completely sequenced plastid genomes currently available and from a selection of mitochondrial genomes representing the major eukaryotic lineages. BLASTN, TBLASTX, or combinations of both are used to locate high-scoring segment pairs (HSPs) between two sequences from which pairwise similarities and distances are computed in different ways resulting in a total of 96 GBDP variants. The suitability of these distance formulae for phylogeny reconstruction is directly estimated by computing a recently described measure of "treelikeness", the so-called delta value, from the respective distance matrices. Additionally, we compare the trees inferred from these matrices using UPGMA, NJ, BIONJ, FastME, or STC, respectively, with the NCBI taxonomy tree of the taxa under study. RESULTS: Our results indicate that, at this taxonomic level, plastid genomes are much more valuable for inferring phylogenies than are mitochondrial genomes, and that distances based on breakpoints are of little use. Distances based on the proportion of "matched" HSP length to average genome length were best for tree estimation. Additionally we found that using TBLASTX instead of BLASTN and, particularly, combining TBLASTX and BLASTN leads to a small but significant increase in accuracy. Other factors do not significantly affect the phylogenetic outcome. The BIONJ algorithm results in phylogenies most in accordance with the current NCBI taxonomy, with NJ and FastME performing insignificantly worse, and STC performing as well if applied to high quality distance matrices. delta values are found to be a reliable predictor of phylogenetic accuracy. CONCLUSION: Using the most treelike distance matrices, as judged by their delta values, distance methods are able to recover all major plant lineages, and are more in accordance with Apicomplexa organelles being derived from "green" plastids than from plastids of the "red" type. GBDP-like methods can be used to reliably infer phylogenies from different kinds of genomic data. A framework is established to further develop and improve such methods. delta values are a topology-independent tool of general use for the development and assessment of distance methods for phylogenetic inference.

Algorithms↗

PubMed2006

Human-specific insertions and deletions inferred from mammalian genome sequences.

It has been suggested that insertions and deletions (indels) have contributed to the sequence divergence between the human and chimpanzee genomes more than do nucleotide changes (3% vs. 1.2%). However, although there have been studies of large indels between the two genomes, no systematic analysis of small indels (i.e., indels </= 100 bp) has been published. In this study, we first estimated that the false-positive rate of small indels inferred from human-chimpanzee pairwise sequence alignments is quite high, suggesting that the chimpanzee genome draft is not sufficiently accurate for our purpose. We have therefore inferred only human-specific indels using multiple sequence alignments of mammalian genomes. We identified >840,000 "small" indels, which affect >7000 UCSC-annotated human genes (>11,000 transcripts). These indels, however, amount to only approximately 0.21% sequence change in the human lineage for the regions compared, whereas in pseudogenes indels contribute to a sequence divergence of 1.40%, suggesting that most of the indels that occurred in genic regions have been eliminated. Functional analysis reveals that the genes whose coding exons have been affected by human-specific indels are enriched in transcription and translation regulatory activities but are underrepresented in catalytic and transporter activities, cellular and physiological processes, and extracellular region/matrix. This functional bias suggests that human-specific indels might have contributed to human unique traits by causing changes at the RNA and protein level.

Animals↗

PubMed1998

Sequence and structure conservation in a protein core.

In order to study structural aspects of sequence conservation in families of homologous proteins, we have analyzed structurally aligned sequences of 585 proteins grouped into 128 homologous families. The conservation of a residue in a family is defined as the average residue similarity in a given position of aligned sequences. The residue similarities were expressed in the form of log-odd substitution tables that take into account the environments of amino acids in three-dimensional structures. The protein core is defined as those residues that have less then 7% solvent accessibility. The density of a protein core is described in terms of atom packing, which is investigated as a criterion for residue substitution and conservation. Although there is no significant correlation between sequence conservation and average atom packing around nonpolar residues such as leucine, valine and isoleucine, a significant correlation is observed for polar residues in the protein core. This may be explained by the hydrogen bonds in which polar residues are involved; the better their protection from water access the more stable should be the structure in that position.

Amino Acid Sequence↗

PubMed1994

Multiple protein structure alignment.

A method was developed to compare protein structures and to combine them into a multiple structure consensus. Previous methods of multiple structure comparison have only concatenated pairwise alignments or produced a consensus structure by averaging coordinate sets. The current method is a fusion of the fast structure comparison program SSAP and the multiple sequence alignment program MULTAL. As in MULTAL, structures are progressively combined, producing intermediate consensus structures that are compared directly to each other and all remaining single structures. This leads to a hierarchic "condensation," continually evaluated in the light of the emerging conserved core regions. Following the SSAP approach, all interatomic vectors were retained with well-conserved regions distinguished by coherent vector bundles (the structural equivalent of a conserved sequence position). Each bundle of vectors is summarized by a resultant, whereas vector coherence is captured in an error term, which is the only distinction between conserved and variable positions. Resultant vectors are used directly in the comparison, which is weighted by their error values, giving greater importance to the matching of conserved positions. The resultant vectors and their errors can also be used directly in molecular modeling. Applications of the method were assessed by the quality of the resulting sequence alignments, phylogenetic tree construction, and databank scanning with the consensus. Visual assessment of the structural superpositions and consensus structure for various well-characterized families confirmed that the consensus had identified a reasonable core.

Amino Acid Sequence↗

PubMed2002

tGRAP, the G-protein coupled receptors mutant database.

The searchable mutant database tGRAP (previously called tinyGRAP) at the University of Tromsø contains data on mutated G-protein coupled receptors (GPCRs). All data have been extracted from scientific papers and entered manually into the database. The current version of the tGRAP mutant database (tGRAP.uit.no, release 10, April 2001) contains around 10 500 mutants extracted from almost 1400 research papers containing mutant data on five families of GPCRs, i.e. Family A, rhodopsin-like; Family B, secretin-like; Family C, metabotropic glutamate-like; Family D, pheromone; Family E, cAMP receptors. A query form provides rapid and simple access to relevant mutant information. In addition to this query form, a tool that enables the user to access mutation data via sequence alignments has been introduced. The ability to access mutant data from such alignments increases the usefulness of the mutant database and facilitates comparison of mutagenesis data between receptors. Moreover, this tool allows the construction of tailor-made sequence alignment views from any combination of receptors belonging to the same class. The database is available at http://tGRAP.uit.no/.

Amino Acid Sequence↗

PubMed2000

Molecular modeling, affinity labeling, and site-directed mutagenesis define the key points of interaction between the ligand-binding domain of the vitamin D nuclear receptor and 1 alpha,25-dihydroxyvitamin D3.

We have combined molecular modeling and classical structure-function techniques to define the interactions between the ligand-binding domain (LBD) of the vitamin D nuclear receptor (VDR) and its natural ligand, 1alpha,25-dihydroxyvitamin D(3) [1alpha,25-(OH)(2)D(3)]. The affinity analogue 1alpha,25-(OH)(2)D(3)-3-bromoacetate exclusively labeled Cys-288 in the VDR-LBD. Mutation of C288 to glycine abolished this affinity labeling, whereas the VDR-LBD mutants C337G and C369G (other conserved cysteines in the VDR-LBD) were labeled similarly to the wild-type protein. These results revealed that the A-ring 3-OH group docks next to C288 in the binding pocket. We further mutated M284 and W286 (separately creating M284A, M284S, W286A, and W286F) and caused severe loss of ligand binding, indicating the crucial role played by the contiguous segment between M284 and C288. Alignment of the VDR-LBD sequence with the sequences of nuclear receptor LBDs of known 3-D structure positioned M284 and W286 in the presumed beta-hairpin of the molecule, thereby identifying it as the region contacting the A-ring of 1alpha, 25-(OH)(2)D(3). From the multiple sequence alignment, we developed a homologous extension model of the VDR-LBD. The model has a canonical nuclear receptor fold with helices H1-H12 and a single beta hairpin but lacks the long insert (residues 161-221) between H2 and H3. We docked the alpha-conformation of the A-ring into the binding pocket first so as to incorporate the above-noted interacting residues. The model predicts hydrogen bonding contacts between ligand and protein at S237 and D299 as well as at the site of the natural mutation R274L. Mutation of S237 or D299 to alanine largely abolished ligand binding, whereas changing K302, a nonligand-contacting residue, to alanine left binding unaffected. In the "activation" helix 12, the model places V418 closest to the ligand, and, consistent with this prediction, the mutation V418S abolished ligand binding. The studies together have enabled us to identify 1alpha,25-(OH)(2)D(3)-binding motifs in the ligand-binding pocket of VDR.

Affinity Labels↗

PubMed2006

Pairwise alignment of protein interaction networks.

With an ever-increasing amount of available data on protein-protein interaction (PPI) networks and research revealing that these networks evolve at a modular level, discovery of conserved patterns in these networks becomes an important problem. Although available data on protein-protein interactions is currently limited, recently developed algorithms have been shown to convey novel biological insights through employment of elegant mathematical models. The main challenge in aligning PPI networks is to define a graph theoretical measure of similarity between graph structures that captures underlying biological phenomena accurately. In this respect, modeling of conservation and divergence of interactions, as well as the interpretation of resulting alignments, are important design parameters. In this paper, we develop a framework for comprehensive alignment of PPI networks, which is inspired by duplication/divergence models that focus on understanding the evolution of protein interactions. We propose a mathematical model that extends the concepts of match, mismatch, and gap in sequence alignment to that of match, mismatch, and duplication in network alignment and evaluates similarity between graph structures through a scoring function that accounts for evolutionary events. By relying on evolutionary models, the proposed framework facilitates interpretation of resulting alignments in terms of not only conservation but also divergence of modularity in PPI networks. Furthermore, as in the case of sequence alignment, our model allows flexibility in adjusting parameters to quantify underlying evolutionary relationships. Based on the proposed model, we formulate PPI network alignment as an optimization problem and present fast algorithms to solve this problem. Detailed experimental results from an implementation of the proposed framework show that our algorithm is able to discover conserved interaction patterns very effectively, in terms of both accuracies and computational cost.

Algorithms↗

PubMed2004

zPicture: dynamic alignment and visualization tool for analyzing conservation profiles.

Comparative sequence analysis has evolved as an essential technique for identifying functional coding and noncoding elements conserved throughout evolution. Here, we introduce zPicture, an interactive Web-based sequence alignment and visualization tool for dynamically generating conservation profiles and identifying evolutionarily conserved regions (ECRs). zPicture is highly flexible, because critical parameters can be modified interactively, allowing users to differentially predict ECRs in comparisons of sequences of different phylogenetic distances and evolutionary rates. We demonstrate the application of this module to identify a known regulatory element in the HOXD locus, in which functional ECRs are difficult to discern against the highly conserved genomic background. zPicture also facilitates transcription factor binding-site analysis via the rVista tool portal. We present an example of the HBB complex when zPicture/rVista combination specifically pinpoints to two ECRs containing GATA-1, NF-E2, and TAL1/E47 binding sites that were identified previously as transcriptional enhancers. In addition, zPicture is linked to the UCSC Genome Browser, allowing users to automatically extract sequences and gene annotations for any recorded locus. Finally, we describe how this tool can be efficiently applied to the analysis of nonvertebrate genomes, including those of microbial organisms.

Animals↗

PubMed2000

Functional importance of motif I of pseudouridine synthases: mutagenesis of aligned lysine and proline residues.

On the basis of sequence alignments, the pseudouridine synthases were grouped into four families that share no statistically significant global sequence similarity, though some common sequence motifs were discovered [Koonin, E. V. (1996) Nucleic Acids. Res. 24, 2411-2415; Gustafsson, C., Reid, R., Greene, P. J., and Santi, D. V. (1996) Nucleic Acids Res. 24, 3756-3762]. We have investigated the functional significance of these alignments by substituting the nearly invariant lysine and proline residues in Motif I of RluA and TruB, pseudouridine synthases belonging to different families. Contrary to our expectations, the altered enzymes display only very mild kinetic impairment. Substitution of the aligned lysine and proline residues does, however, reduce structural stability, consistent with a temperature sensitive phenotype that results from substitution of the cognate proline residue in Cbf5p, a yeast homologue of TruB [Zerbarjadian, Y., King, T., Fournier, M. J., Clarke, L., and Carbon, J. (1999) Mol. Cell. Biol. 19, 7461-7472]. Together, our data support a functional role for Motif I, as predicted by sequence alignments, though the effect of substituting the highly conserved residues was milder than we anticipated. By extrapolation, our findings also support the assignment of pseudouridine synthase function to certain physiologically important eukaryotic proteins that contain Motif I, including the human protein dyskerin, alteration of which leads to the disease dyskeratosis congenita.

Amino Acid Motifs↗

PubMed1999

Multiple alignment and sorting of peptides derived from phage-displayed random peptide libraries with polyclonal sera allows discrimination of relevant phagotopes.

Biopanning of phage-displayed random peptide libraries is a powerful technique for identifying peptides that mimic epitopes (mimotopes) for monoclonal antibodies (mAbs). However, peptides derived using polyclonal antisera may represent epitopes for a diverse range of antibodies. Hence following screening of phage libraries with polyclonal antisera, including autoimmune disease sera, a procedure is required to distinguish relevant from irrelevant phagotopes. We therefore applied the multiple sequence alignment algorithm PILEUP together with a matrix for scoring amino acid substitutions based on physicochemical properties to generate guide trees depicting relatedness of selected peptides. A random heptapeptide library was biopanned nine times using no selecting antibodies, immunoglobulin G (IgG) from sera of subjects with autoimmune diseases (primary biliary cirrhosis (PBC) and type 1 diabetes) and three murine ascites fluids that contained mAbs to overlapping epitope(s) on the Ross River Virus envelope protein 2. Peptides randomly sampled from the library were distributed throughout the guide tree of the total set of peptides whilst many of the peptides derived in the absence of selecting antibody aligned to a single cluster. Moreover peptides selected by different sources of IgG aligned to separate clusters, each with a different amino acid motif. These alignments were validated by testing all of the 53 phagotopes derived using IgG from PBC sera for reactivity by capture ELISA with antibodies affinity purified on the E2 subunit of the pyruvate dehydrogenase complex (PDC-E2), the major autoantigen in PBC: only those phagotopes that aligned to PBC-associated clusters were reactive. Hence the multiple sequence alignment procedure discriminates relevant from irrelevant phagotopes and thus a major difficulty with biopanning phage-displayed random peptide libraries with polyclonal antibodies is surmounted.

Algorithms↗

PubMed2006

PartTree: an algorithm to build an approximate tree from a large number of unaligned sequences.

MOTIVATION: To construct a multiple sequence alignment (MSA) of a large number (> approximately 10,000) of sequences, the calculation of a guide tree with a complexity of O(N2) to O(N3), where N is the number of sequences, is the most time-consuming process. RESULTS: To overcome this limitation, we have developed an approximate algorithm, PartTree, to construct a guide tree with an average time complexity of O(N log N). The new MSA method with the PartTree algorithm can align approximately 60,000 sequences in several minutes on a standard desktop computer. The loss of accuracy in MSA caused by this approximation was estimated to be several percent in benchmark tests using Pfam. AVAILABILITY: The present algorithm has been implemented in the MAFFT sequence alignment package (http://align.bmr.kyushu-u.ac.jp/mafft/software/). SUPPLEMENTARY INFORMATION: Supplementary information is available at Bioinformatics online.

Algorithms↗

PubMed2003

3MATRIX and 3MOTIF: a protein structure visualization system for conserved sequence motifs.

Computational methods such as sequence alignment and motif construction are useful in grouping related proteins into families, as well as helping to annotate new proteins of unknown function. These methods identify conserved amino acids in protein sequences, but cannot determine the specific functional or structural roles of conserved amino acids without additional study. In this work, we present 3MATRIX (http://3matrix.stanford.edu) and 3MOTIF (http://3motif.stanford.edu), a web-based sequence motif visualization system that displays sequence motif information in its appropriate three-dimensional (3D) context. This system is flexible in that users can enter sequences, keywords, structures or sequence motifs to generate visualizations. In 3MOTIF, users can search using discrete sequence motifs such as PROSITE patterns, eMOTIFs, or any other regular expression-like motif. Similarly, 3MATRIX accepts an eMATRIX position-specific scoring matrix, or will convert a multiple sequence alignment block into an eMATRIX for visualization. Each query motif is used to search the protein structure database for matches, in which the motif is then visually highlighted in three dimensions. Important properties of motifs such as sequence conservation and solvent accessible surface area are also displayed in the visualizations, using carefully chosen color shading schemes.

Amino Acid Motifs↗

PubMed1984

Fast optimal alignment.

We show how to speed up sequence alignment algorithms of the type introduced by Needleman and Wunsch (and generalized by Sellers and others). Faster alignment algorithms have been introduced, but always at the cost of possibly getting sub-optimal alignments. Our modification results in the optimal alignment still being found, often in 1/10 the usual time. What we do is reorder the computation of the usual alignment matrix so that the optimal alignment is ordinarily found when only a small fraction of the matrix is filled. The number of matrix elements which have to be computed is related to the distance between the sequences being aligned; the better the optimal alignment, the faster the algorithm runs.

Base Sequence↗

PubMed2001

Evidence for an independent radiation of endosymbiotic litostome ciliates within Australian marsupial herbivores.

The recent discovery of isotrichid-like ciliates occurring as endosymbionts in macropodid marsupials posed interesting questions in regard to both their phyletic origin (all previous records confined to eutherian mammals) and their morphological evolution (Australian forms possibly representing missing links between previously described genera). The SSU rRNA gene was sequenced for three species (Dasytricha dehorityi, D. dogieli, and Bitricha tasmaniensis) and aligned against representatives of all major ciliate classes. The Australian species did not group with the other isotrichid species but instead formed an independent radiation. Discrepancies between recent global phylogenies of the phylum Ciliophora were examined by manipulation of the aligned sequence data set. Sources of conflict between these studies did not stem from differences in outgroup choice or phylogenetic reconstruction methods. Differences in the application of confidence limits and primary sequence alignment have probably resulted in the reporting of spurious associations which are not supported by more conservative confidence or alignment methodology. At present, the ciliate subphylum Intramacronucleata is an unresolved polytomy which may be due to deficiencies in the SSU rRNA gene sequence dataset or indicate that the ciliates radiated into their extant classes by rapid burst-like evolution.

Animals↗

PubMed1992

Is the bacteriophage lambda lysozyme an evolutionary link or a hybrid between the C and V-type lysozymes? Homology analysis and detection of the catalytic amino acid residues.

The relationship between the bacteriophage lambda lysozyme (lambda L) and the C and V-type lysozymes has been investigated by sequence alignment, secondary structure prediction and pattern recognition methods. The alignment of the amino terminal part of lambda L with that of V-type lysozymes suggests that Glu19 is a residue essential for catalysis. Its mutation to Gln leads to a completely inactive enzyme. In the alignment of the sequence of lambda L with those of the C-type lysozymes a strongly homologous fragment of about 30 amino acid residues is detected. Taking into consideration this observation and the published structural alignments between C and V-type lysozymes, a repetition of the beta-sheet motif in lambda L is proposed. The multiple alignment draws the attention to a possible catalytic role for Asp34 that would be positioned in the middle of the second strand of the beta-sheet as in the C-type lysozymes. This role is confirmed by mutagenesis. The implications of these observations in terms of the evolutionary relationship between lambda L and the other lysozymes is discussed.

Amino Acid Sequence↗

PubMed1993

A protein alignment scoring system sensitive at all evolutionary distances.

Protein sequence alignments generally are constructed with the aid of a "substitution matrix" that specifies a score for aligning each pair of amino acids. Assuming a simple random protein model, it can be shown that any such matrix, when used for evaluating variable-length local alignments, is implicitly a "log-odds" matrix, with a specific probability distribution for amino acid pairs to which it is uniquely tailored. Given a model of protein evolution from which such distributions may be derived, a substitution matrix adapted to detecting relationships at any chosen evolutionary distance can be constructed. Because in a database search it generally is not known a priori what evolutionary distances will characterize the similarities found, it is necessary to employ an appropriate range of matrices in order not to overlook potential homologies. This paper formalizes this concept by defining a scoring system that is sensitive at all detectable evolutionary distances. The statistical behavior of this scoring system is analyzed, and it is shown that for a typical protein database search, estimating the originally unknown evolutionary distance appropriate to each alignment costs slightly over two bits of information, or somewhat less than a factor of five in statistical significance. A much greater cost may be incurred, however, if only a single substitution matrix, corresponding to the wrong evolutionary distance, is employed.

Algorithms↗

PubMed2007

Modeling the tertiary structure of the patatin domain of neuropathy target esterase.

Neuropathy target esterase (NTE) is a transmembrane protein of unknown function whose specific chemical modification by certain organophosphorus (OP) compounds leads to distal axonopathy. Therefore, solving the 3D structure of NTE would advance the understanding of its pathogenic and physiologic roles. In this study, the tertiary structures of the patatin (catalytic) domain and the N-terminal transmembrane domain of NTE were modeled using the crystal structures of patatin (PDB ID 1oxw) and moricin (PDB ID 1kv4) as templates. Sequence alignments and secondary structure predictions were obtained from the INUB server (Buffalo, NY). O and PyMol were used to build the PNTE and NTE TMD chains from these sequence alignments. The PNTE model was refined in the presence of water using the crystallography and NMR system, while the NTE TMD model was refined in vacuo using the GROMOS implementation in the Swiss PDB viewer. The modeled active site of NTE was found to consist of a Ser966-Asp1086 catalytic dyad, which is characteristic of phospholipase A2 enzymes. The Ser966 Ogamma was located 2.93 A from the Odelta2 of Asp1086. In addition, our NTE model was found to contain a single N-terminal transmembrane domain. This modeling effort provided structural and mechanistic predictions about the catalytic domain of NTE that are being verified via experimental techniques.

Amino Acid Sequence↗

PubMed2006

Characterization of OCA2 cDNA in different porcine breeds and analysis of its potential effect on skin pigmentation in a red Iberian strain.

Although the function of the OCA2 gene product has not been totally clarified, variation in OCA2 has been associated with skin and hair pigmentation in human and mouse. However, its contribution to skin colour in domestic species has not been reported. In this study, cDNA and intron 9 sequences of the porcine OCA2 gene have been characterized in several pig populations. The cDNA sequence alignment of 20 animals from eight porcine populations allowed the identification of 10 single nucleotide polymorphisms (SNPs); five of the 10 SNPs were non-synonymous. The intron 9 sequence alignment of 12 animals belonging to four pig populations revealed four additional SNPs. Skin colour variation was analysed in a red strain of Iberian pigs with segregation of three SNPs forming two OCA2 intragenic haplotypes. Results from this study provide evidence of a suggestive dominant effect of haplotypes on colour intensity and indicate an important contribution of additive polygenic effects (h2 = 0.56 +/- 0.21) to the variance of this trait.

Animals↗