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,531 records · Page 85Linked to original sources

PubMed2026

Exploring differences across pangenome-graph representations using Escherichia coli O157:H7 as a model.

Pangenome graphs are increasingly used to represent population-scale bacterial diversity, yet construction methods span fundamentally different representation paradigms whose outputs and sensitivities to assembly quality remain poorly quantified. We systematically reviewed microbial pangenome graph tools and benchmarked seven representative methods spanning gene-cluster, compacted coloured de Bruijn graph, one hybrid approach and one multiple sequence alignment method. Using a repeat-rich Escherichia coli O157:H7 dataset with complete genomes and matched short-read data, we constructed graphs from identical inputs and observed orders-of-magnitude differences in graph size and fragmentation, indicating that global topology is driven by representation strategy. Varying completeness composition revealed that assembly fragmentation is a first-order determinant of graph structure: gene-cluster graphs contracted as draft assemblies replaced complete genomes, whereas compacted coloured de Bruijn graphs expanded, with distinct degree-prevalence fingerprints across tools. In contrast, the multiple sequence alignment method could not be evaluated across fragmented inputs because it did not run reliably on draft-assembly datasets. Computational cost mirrored these shifts and depended strongly on completeness composition, including a pronounced runtime penalty for one compacted coloured de Bruijn graph implementation on all-draft inputs. Finally, analysis of Shiga toxin loci showed that pangenome-level reconciliation by gene-cluster-based tools does not reliably correct assembly artefacts at challenging multi-copy genes and that performance varies by locus. Together, these findings show that pangenome graphs are representation-dependent models of bacterial diversity, and that, in this repeat-rich O157:H7 benchmark dataset, assembly completeness is a primary determinant of their topology, scalability, and locus-level accuracy.

Escherichia coli O157↗

PubMed1999

Cellobiose-6-phosphate hydrolase (CelF) of Escherichia coli: characterization and assignment to the unusual family 4 of glycosylhydrolases.

The gene celF of the cryptic cel operon of Escherichia coli has been cloned, and the encoded 6-phospho-beta-glucosidase (cellobiose-6-phosphate [6P] hydrolase; CelF [EC 3.2.1.86]) has been expressed and purified in a catalytically active state. Among phospho-beta-glycosidases, CelF exhibits unique requirements for a divalent metal ion and NAD(+) for activity and, by sequence alignment, is assigned to family 4 of the glycosylhydrolase superfamily. CelF hydrolyzed a variety of P-beta-glucosides, including cellobiose-6P, salicin-6P, arbutin-6P, gentiobiose-6P, methyl-beta-glucoside-6P, and the chromogenic analog, p-nitrophenyl-beta-D-glucopyranoside-6P. In the absence of a metal ion and NAD(+), purified CelF was rapidly and irreversibly inactivated. The functional roles of the cofactors have not been established, but NAD(+) appears not to be a reactant and there is no evidence for reduction of the nucleotide during substrate cleavage. In solution, native CelF exists as a homotetramer (M(w), approximately 200,000) composed of noncovalently linked subunits, and this oligomeric structure is maintained independently of the presence or absence of a metal ion. The molecular weight of the CelF monomer (M(r), approximately 50,000), estimated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis, is in agreement with that calculated from the amino acid sequence of the polypeptide (450 residues; M(r) = 50,512). Comparative sequence alignments provide tentative identification of the NAD(+)-binding domain (residues 7 to 40) and catalytically important glutamyl residues (Glu(112) and Glu(356)) of CelF.

Amino Acid Sequence↗

PubMed2001

Molecular modeling of the extracellular domain of the RET receptor tyrosine kinase reveals multiple cadherin-like domains and a calcium-binding site.

Using bioinformatic tools, mutagenesis, and binding studies, we have investigated the structural organization of the extracellular region of the RET receptor tyrosine kinase, a functional receptor for glial cell line-derived neurotrophic factor (GDNF). Multiple sequence alignments of seven vertebrate sequences and one invertebrate RET sequence delineated four distinct N-terminal domains, each of about 110 residues, containing many of the consensus motifs of the cadherin fold. Based on these alignments and the crystal structures of epithelial and neural cadherins, we have generated molecular models of each of the four cadherin-like domains in the extracellular region of human RET. The modeled structures represent realistic models from both energetic and geometrical points of view and are consistent with previous observations gathered from biochemical analyses of the effects of Hirschsprung's disease mutations affecting the folding and stability of the RET molecule, as well as our own site-directed mutagenesis studies of RET cadherin-like domain 1. We have also investigated the role of Ca(2+) in ligand binding by RET and found that Ca(2+) ions are required for RET binding to GDNF but not for GDNF binding to the GFRalpha1 co-receptor. In agreement with these results, RET, but not GFRalpha1, was found to bind Ca(2+) directly. Our results indicate that the overall architecture of the extracellular region of RET is more closely related to cadherins than previously thought. The models of the cadherin-like domains of human RET represent valuable tools with which to guide future site-directed mutagenesis studies aimed at identifying residues involved in ligand binding and receptor activation.

Amino Acid Sequence↗

PubMed2003

Calpain silencing by a reversible intrinsic mechanism.

Uncontrolled activation of calpain can lead to necrotic cell death and irreversible tissue damage. We have discovered an intrinsic mechanism whereby the autolysis-generated protease core fragment of calpain is inactivated through the inherent instability of a key alpha-helix. This auto-inactivation state was captured by the 1.9 A Ca(2+)-bound structure of the protease core from m-calpain, and sequence alignments suggest that it applies to about half of the calpain isoforms. Intact calpain large subunits are also subject to this inhibition, which can be prevented through assembly of the heterodimers. Other isoforms or their released cores are not silenced by this mechanism and might contribute to calpain patho-physiologies.

Amino Acid Sequence↗

PubMed1992

Phylogenetic classification of human papillomaviruses: correlation with clinical manifestations.

Human papillomaviruses (HPVs) are a heterogeneous group of small dsDNA viruses which cause a variety of proliferative epithelial lesions at specific anatomical sites. Although more than 65 different virus types have been cloned and characterized, no uniform classification system exists. In order to classify HPV DNA types, phylogenetic trees were constructed based on nucleotide sequence alignments using parsimony and distance matrix algorithms. The resulting phylogenetic trees provide a classification of the HPVs into specific groups encompassing the known tissue tropism and oncogenic potential of each HPV type. The implications of a phylogenetic taxonomy on the diagnostic detection of HPVs and the concept of different HPV species are discussed.

Algorithms↗

PubMed2004

Primaclade--a flexible tool to find conserved PCR primers across multiple species.

Primaclade is a web-based application that accepts a multiple species nucleotide alignment file as input and identifies a set of polymerase chain reaction (PCR) primers that will bind across the alignment. Primaclade iteratively runs the Primer3 application for each alignment sequence and collates the results. Primaclade creates an HTML results page that recaps the original alignment, provides a consensus sequence and lists primers for each alignment area, with primers color-coded to reflect the level of degeneracy in the primer.

Algorithms↗

PubMed2025

SegMantX: A Novel Tool for Detecting DNA Duplications Uncovers Prevalent Duplications in Plasmids.

Segmental duplications play an important role in genome evolution via their contribution to copy-number variation, gene-family diversification, and the emergence of novel functions. The detection of segmental duplications is challenging due to heterogeneous amelioration of sequence similarity among duplicates, which hinders the reconstruction of continuous sequence alignment. Here we introduce SegMantX, a novel approach for the identification of diverged segmental duplications in prokaryote genomes using local alignment chaining. In this approach, local alignments resulting from a preliminary sequence similarity search (e.g. BLASTn) are chained into continuous segments. Evaluating the performance of SegMantX using simulated sequences shows that the tool can detect diverged duplications beyond the sensitivity limits of standard alignment-based methods. Applying SegMantX to 6,784 enterobacterial plasmids, we find that 65% plasmids contain duplicated regions and gene duplications, most of which correspond either to dispersed, noncoding regions or duplicated mobile genetic elements (MGEs; e.g. transposons and insertion sequences). Furthermore, we demonstrate the applicability of SegMantX for the identification of diverged gene transfers between replicons and plasmid hybridization events. Our findings highlight MGEs as drivers of segmental duplications in plasmid evolution, leading to the amplification of their cargo genes, including antibiotic resistance genes. SegMantX provides a powerful framework for reconstructing diverged segmental duplications and other alignment problems.

Plasmids↗

PubMed1998

Identification and characterization of the KlCMD1 gene encoding Kluyveromyces lactis calmodulin.

The KlCMD1 gene was isolated from a Kluyveromyces lactis genomic library as a suppressor of the Saccharomyces cerevisiae temperature-sensitive mutant spc110-124, an allele previously shown to be suppressed by elevated copy number of the S. cerevisiae calmodulin gene CMD1. The KlCMD1 gene encodes a polypeptide which is 95% identical to S. cerevisiae calmodulin and 55% identical to calmodulin from Schizosaccharomyces pombe. Complementation of a S. cerevisiae cdm1 deletion mutant by KlCMD1 demonstrates that this gene encodes a functional calmodulin homologue. Multiple sequence alignment of calmodulins from yeast and multicellular eukaryotes shows that the K. lactis and S. cerevisiae calmodulins are considerably more closely related to each other than to other calmodulins, most of which have four functional Ca2+-binding EF hand domains. Thus like its S. cerevisiae counterpart Cmd1p, the KlCMD1 product is predicted to form only three Ca2+-binding motifs.

Amino Acid Sequence↗

PubMed2000

Molecular typing of Toxoplasma gondii strains by GRA6 gene sequence analysis.

The utility of sequence polymorphisms in the dense granule antigen GRA6 gene as typing markers for Toxoplasma gondii was investigated. The coding region of GRA6 was amplified, sequenced and compared for 30 Toxoplasma strains from eight different zymodemes (Z1-Z8). Sequence alignment identified nucleotide polymorphisms at 24 positions out of 690 bp, which correlated with murine-virulence. Types I, II, and III could be distinguished from each other on the basis of three, 10, and six variable positions, respectively. Two deletions of 15 bp and 3 bp existed in the avirulent (type II) strains. With one exception, all polymorphic positions resulted in amino acid substitutions, and the two gaps of 15 bp and 3 bp caused the deletion of six amino acids in type II strains. Intra-specific polymorphisms were also found in the virulent group. A high degree of sequence polymorphism correlating with the phenotypes of T. gondii strains points to the GRA6 gene being a good marker for strain characterisation and typing of the isolates of this apicomplexan. The large variety of amino acid changes supports the view that the GRA6 protein plays an important role in the antigenicity and pathogenicity of T. gondii. The existence of polymorphic restriction sites for endonuclease MseI was used to develop a PCR-RFLP method which could simply differentiate the three different groups (types I, II, III) of T. gondii.

Animals↗

PubMed2005

Transactivation of Abl by the Crk II adapter protein requires a PNAY sequence in the Crk C-terminal SH3 domain.

To gain a better understanding of how Crk II regulates the function of the Abl tyrosine kinase, we explored the function of the C-terminal linker and SH3 domain, a region of Crk II that is still poorly understood. Molecular modeling, tryptophan fluorescence, and covariation sequence alignment indicate that the Crk-SH3-C has a unique binding groove and RT loop not observed in typical SH3 domains. Based on these models, we made a series of mutations in the linker and in residues predicted to destabilize the putative binding pocket and RT loop. In Abl transactivation assays, Y222F and P225A mutations in the linker resulted in strong transactivation of Abl by Crk II. However, mutations predicted to be at the surface of the Crk SH3-C were not activators of Abl. Interestingly, combinations of activating mutations of Crk II with mutations in the highly conserved PNAY sequence in the SH3-C inactivated the activating mutations, suggesting that the SH3-C is necessary for activation. Our data provide insight into the role of highly conserved residues in the Crk-SH3-C, suggesting a mechanism for how the linker and the Crk-SH3-C function in the transactivation of the Abl tyrosine kinase.

Amino Acid Sequence↗

PubMed2005

Aflatoxin biosynthesis gene clusters and flanking regions.

AIMS: To compare the biosynthetic gene cluster sequences of the main aflatoxin (AF)-producing Aspergillus species. METHODS AND RESULTS: Sequencing was on fosmid clones selected by homology to Aspergillus parasiticus sequence. Alignments revealed that gene order is conserved among AF gene clusters of Aspergillus nomius, A. parasiticus, two sclerotial morphotypes of Aspergillus flavus, and an unnamed Aspergillus sp. Phylogenetic relationships were established using the maximum likelihood method implemented in PAUP. Based on the Eurotiomycete/Sordariomycete divergence time, the A. flavus-type cluster has been maintained for at least 25 million years. Such conservation of the genes and gene order reflects strong selective constraints on rearrangement. Phylogenetic comparison of individual genes in the cluster indicated that ver-1, which has homology to a melanin biosynthesis gene, experienced selective forces distinct from the other pathway genes. Sequences upstream of the polyketide synthase-encoding gene vary among the species, but a four-gene sugar utilization cluster at the distal end is conserved, indicating a functional relationship between the two adjacent clusters. CONCLUSIONS: The high conservation of cluster components needed for AF production suggests there is an adaptive value for AFs in character-shaping niches important to those taxa. SIGNIFICANCE AND IMPACT OF THE STUDY: This is the first comparison of the complete nucleotide sequences of gene clusters harbouring the AF biosynthesis genes of the main AF-producing species. Such a comparison will aid in understanding how AF biosynthesis is regulated in experimental and natural environments.

Aflatoxins↗

PubMed2000

Identification and characterization of a shrimp white spot syndrome virus (WSSV) gene that encodes a novel chimeric polypeptide of cellular-type thymidine kinase and thymidylate kinase.

From previously constructed genomic libraries of a Taiwan WSSV isolate, a putative WSSV tk-tmk gene was identified. Uniquely, the open reading frame (ORF) of this gene was predicted to encode a novel chimeric protein of 388 amino acids with significant homology to two proteins: thymidine kinase (TK) and thymidylate kinase (TMK). Northern blot analysis with a WSSV tk-tmk-specific riboprobe detected a major transcript of 1.6 kb. When healthy adult Penaeus monodon shrimp were inoculated with WSSV, the tk-tmk gene transcript was first detected by RT-PCR analysis at 4 h postinfection and transcription levels continued to increase over the first 18 h. The gene's major in vitro transcription and translation product, equivalent to the predicted size (43 kDa), is a single chimeric protein that includes both the TK and TMK functional motifs. Evidence for phylogenetic analysis and sequence alignment suggested that the gene may have resulted from the fusion of a cellular-type TK gene and a cellular-type TMK gene. Its unique arrangement may also provide a valuable gene marker for WSSV.

Amino Acid Sequence↗

PubMed2004

Genomic and phylogenic comparisons of the alpha-globin and beta-globin intergenic sequences between zebra fish (Danio rerio) and six closely related Cyprinindae species.

The alpha-globin and beta-globin genes of the zebrafish are tightly linked on the same chromosome in a 3'-5' and 5'-3' configuration, respectively. Although the location of the controlling sequences has been mapped to the intergenic region, analysis determining the uniqueness of this unusual arrangement to zebrafish has not been undertaken. To explore this, we isolated, sequenced, and phylogenetically analyzed the intergenic region between globin gene families of seven Cyprinindae species including zebrafish. These species were grouped into an in group (immediate relatives, not so distant relatives), and an out group (distant relative). Cellulose acetate electrophoresis of hemoglobin (Hb) detected multiple variants in each species, but a band with electrophoretic mobility (EM) of 6.7 x 10(-5) cm(2).volt(-1).sec(-1) was shared between species. Polymerase chain reaction (PCR) amplification of the intergenic globin gene region also detected a 1.0-kb fragment that was repeated in the in group and a 1.2-kb fragment in the out group. Sequence comparison confirmed that the genetic orientation and controlling sequences location were conserved throughout this region in all seven species. This phylogenic footprinting indicated that the configuration was not exclusive to zebrafish. To confirm sequence alignment, maximum parsimony phylogenic analysis, was performed. Only one member of that group the giant danio, was not closely clustered, being located almost equidistance between the immediate relative and the species of the other clusters. This may represent an ancestral configuration prior to transposition of the alpha globin and beta-globin genes families to nonsynteny.

Algorithms↗

PubMed1998

A computer program for aligning a cDNA sequence with a genomic DNA sequence.

We address the problem of efficiently aligning a transcribed and spliced DNA sequence with a genomic sequence containing that gene, allowing for introns in the genomic sequence and a relatively small number of sequencing errors. A freely available computer program, described herein, solves the problem for a 100-kb genomic sequence in a few seconds on a workstation.

Algorithms↗

PubMed1999

Ab initio folding of proteins using restraints derived from evolutionary information.

We present our predictions in the ab initio structure prediction category of CASP3. Eleven targets were folded, using a method based on a Monte Carlo search driven by secondary and tertiary restraints derived from multiple sequence alignments. Our results can be qualitatively summarized as follows: The global fold can be considered "correct" for targets 65 and 74, "almost correct" for targets 64, 75, and 77, "half-correct" for target 79, and "wrong" for targets 52, 56, 59, and 63. Target 72 has not yet been solved experimentally. On average, for small helical and alpha/beta proteins (on the order of 110 residues or smaller), the method predicted low resolution structures with a reasonably good prediction of the global topology. Most encouraging is that in some situations, such as with target 75 and, particularly, target 77, the method can predict a substantial portion of a rare or even a novel fold. However, the current method still fails on some beta proteins, proteins over the 110-residue threshold, and sequences in which only a poor multiple sequence alignment can be built. On the other hand, for small proteins, the method gives results of quality at least similar to that of threading, with the advantage of not being restricted to known folds in the protein database. Overall, these results indicate that some progress has been made on the ab initio protein folding problem. Detailed information about our results can be obtained by connecting to http:/(/)www.bioinformatics.danforthcenter.org/+ ++CASP3.

Algorithms↗

PubMed2005

Transmembrane domain prediction and consensus sequence identification of the oligopeptide transport family.

Few polytopic membrane proteins have had their topology determined experimentally. Often, researchers turn to an algorithm to predict where the transmembrane domains might lie. Here we use a consensus method, using six different transmembrane domain prediction algorithms on six members of the oligopeptide transport family, all of which have been experimentally characterized. PSI-BLAST results indicate that the six chosen oligopeptide transport family members are distributed throughout most branches of the phylogram, suggesting that these members represent a broad view of the oligopeptide transport family. We combined the prediction algorithms with a multiple sequence alignment, and consensus transmembrane domains were assigned not only based on algorithmic output, but also based on conserved familial motifs found by analysis of the PSI-BLAST results. The consensus method combined with the "charge-difference rule" yields a model topology for the family containing 12 transmembrane domains with the N- and C-termini facing extracellular.

Algorithms↗

PubMed1998

Identification of enhancer-like elements in the ribosomal intergenic spacer of Odontophrynus americanus 2n and 4n (Amphibia, Anura).

Repetitive elements are found in the ribosomal intergenic spacer (IGS) of most organisms. A particularly complex pattern of internal repetition occurs in the IGSs of O. americanus 2n and 4n, which are composed of several types of BamHI subrepeats (B-SRs). The most repetitive one is approximately 87 bp long, and is highly represented in the IGS variants of these amphibians. Sequence analyses of six diploid and two tetraploid B-SRs show 87% and 86% homology, respectively, and related secondary structure predictions. The comparison of the 2n and 4n B-SR sequences aligned with the 81 bp enhancer of Xenopus laevis reveals 36% homology. Furthermore, other B-SR features like size, number, and secondary structures resemble those of Xenopus enhancers, suggesting that B-SRs may function as regulators of O. americanus rDNA transcription. The present data also corroborate the close evolutionary relationship between 2n and 4n O. americanus species.

Animals↗

PubMed2006

The WWWH of remote homolog detection: the state of the art.

The detection of remote homolog pairs of proteins using computational methods is a pivotal problem in structural bioinformatics, aiming to compute protein folds on the basis of information in the database of known structures. In the last 25 years, several methods have been developed to tackle this problem, based on different approaches including sequence-sequence alignments and/or structure comparison. In this article, we will briefly discuss When, Why, Where and How (WWWH) to perform remote homology search, reviewing some of the most widely adopted computational approaches. The specific aim is highlighting the basic criteria implemented by different research groups and commenting on the status of the art as well as on still-open questions.

Algorithms↗