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

PubMed1990

Homologous sugar transport proteins in Escherichia coli and their relatives in both prokaryotes and eukaryotes.

Separate proteins for proton-linked transport of D-xylose, L-arabinose, D-galactose, L-rhamnose and L-fucose into Escherichia coli are being studied. By cloning and sequencing the appropriate genes, the amino acid sequences of proteins for D-xylose/H+ symport (XylE), L-arabinose/H+ symport (AraE), and part of the protein for D-galactose/H+ symport (GalP) have been determined. These are homologous, with at least 28% identical amino acid residues conserved in the aligned sequences, although their primary sequences are not similar to those of other E. coli transport proteins for lactose, melibiose, or D-glucose. However, they are equally homologous to the passive D-glucose transport proteins from yeast, rat brain, rat adipocytes, human erythrocytes, human liver, and a human hepatoma cell line. The substrate specificity of GalP from E. coli is similar to that of the mammalian glucose transporters. Furthermore, the activities of GalP, AraE and the mammalian glucose transporters are all inhibited by cytochalasin B and N-ethylmaleimide. Conserved residues in the aligned sequences of the bacterial and mammalian transporters are identified, and the possible roles of some in sugar binding, cation binding, cytochalasin binding, and reaction with N-ethylmaleimide are discussed. Each protein is independently predicted to form 12 hydrophobic, membrane-spanning alpha-helices with a central hydrophilic segment, also comprised of alpha-helix. This unifying structural model of the sugar transporters shares features with other ion-linked transport proteins for citrate or tetracycline.

Amino Acid Sequence↗

PubMed1996

The C2 domain calcium-binding motif: structural and functional diversity.

The C2 domain is a Ca(2+)-binding motif of approximately 130 residues in length originally identified in the Ca(2+)-dependent isoforms of protein kinase C. Single and multiple copies of C2 domains have been identified in a growing number of eukaryotic signalling proteins that interact with cellular membranes and mediate a broad array of critical intracellular processes, including membrane trafficking, the generation of lipid-second messengers, activation of GTPases, and the control of protein phosphorylation. As a group, C2 domains display the remarkable property of binding a variety of different ligands and substrates, including Ca2+, phospholipids, inositol polyphosphates, and intracellular proteins. Expanding this functional diversity is the fact that not all proteins containing C2 domains are regulated by Ca2+, suggesting that some C2 domains may play a purely structural role or may have lost the ability to bind Ca2+. The present review summarizes the information currently available regarding the structure and function of the C2 domain and provides a novel sequence alignment of 65 C2 domain primary structures. This alignment predicts that C2 domains form two distinct topological folds, illustrated by the recent crystal structures of C2 domains from synaptotagmin 1 and phosphoinositide-specific phospholipase C-delta 1, respectively. The alignment highlights residues that may be critical to the C2 domain fold or required for Ca2+ binding and regulation.

Amino Acid Sequence↗

PubMed2005

Prediction of functional specificity determinants from protein sequences using log-likelihood ratios.

MOTIVATION: A number of methods have been developed to predict functional specificity determinants in protein families based on sequence information. Most of these methods rely on pre-defined functional subgroups. Manual subgroup definition is difficult because of the limited number of experimentally characterized subfamilies with differing specificity, while automatic subgroup partitioning using computational tools is a non-trivial task and does not always yield ideal results. RESULTS: We propose a new approach SPEL (specificity positions by evolutionary likelihood) to detect positions that are likely to be functional specificity determinants. SPEL, which does not require subgroup definition, takes a multiple sequence alignment of a protein family as the only input, and assigns a P-value to every position in the alignment. Positions with low P-values are likely to be important for functional specificity. An evolutionary tree is reconstructed during the calculation, and P-value estimation is based on a random model that involves evolutionary simulations. Evolutionary log-likelihood is chosen as a measure of amino acid distribution at a position. To illustrate the performance of the method, we carried out a detailed analysis of two protein families (LacI/PurR and G protein alpha subunit), and compared our method with two existing methods (evolutionary trace and mutual information based). All three methods were also compared on a set of protein families with known ligand-bound structures. AVAILABILITY: SPEL is freely available for non-commercial use. Its pre-compiled versions for several platforms and alignments used in this work are available at ftp://iole.swmed.edu/pub/SPEL/

Algorithms↗

PubMed2000

hNMP 200: a novel human common nuclear matrix protein combining structural and regulatory functions.

Previously we have reported about human nuclear matrix proteins (hNMPs) with increased reassembling and potential filament-forming capability [C. Gerner et al., 1999, J. Cell. Biochem. 74, 145-151]. Here, we cloned the cDNA of one of these proteins, hNMP 200, following partial amino acid sequencing of the novel 56-kDa nuclear protein. Sequence alignments show hNMP 200-related proteins in metazoans, plants, and yeast, the homologous Saccharomyces cerevisiae protein prp19 being an accessory, but essential, factor for pre-mRNA processing. Evidence for any enzymatic activity was not detected. However, the hNMP 200 primary sequence contained five consensus WD-repeat sequences, indicative of participation and regulatory function in larger protein assemblies. Northern blot analysis and 2D protein electrophoresis showed ubiquitous expression of hNMP 200 in a variety of cell types. (35)S labeling studies indicated a high metabolic stability of the protein. The hNMP 200 gene was assigned to chromosomal region 11q12.2. Confocal laser scanning microscopy revealed that the intracellular localization conformed with that reported for other structural nuclear proteins. In interphase cells, green fluorescent protein-tagged hNMP 200 was predominantly nucleoplasmic. Structures with speckled appearance extended through several sections of in situ-isolated nuclear matrices. During cell division hNMP 200 became irregularly distributed in prophase, sparing regions of condensing chromatin. In anaphase it was concentrated in the spindle midzone. The putative dual function of the novel NMP is discussed. Being a component of the nuclear framework, it may provide structural support for components of the RNA-processing machinery, thereby also modulating splicing activities.

Amino Acid Sequence↗

PubMed2000

Use of mitogenomic information in teleostean molecular phylogenetics: a tree-based exploration under the maximum-parsimony optimality criterion.

We explored the phylogenetic utility and limits of the individual and concatenated mitochondrial genes for reconstructing the higher-level relationships of teleosts, using the complete (or nearly complete) mitochondrial DNA sequences of eight teleosts (including three newly determined sequences), whose relative phylogenetic positions were noncontroversial. Maximum-parsimony analyses of the nucleotide and amino acid sequences of 13 protein-coding genes from the above eight teleosts, plus two outgroups (bichir and shark), indicated that all of the individual protein-coding genes, with the exception of ND5, failed to recover the expected phylogeny, although unambiguously aligned sequences from 22 concatenated transfer RNA (tRNA) genes (stem regions only) recovered the expected phylogeny successfully with moderate statistical support. The phylogenetic performance of the 13 protein-coding genes in recovering the expected phylogeny was roughly classified into five groups, viz. very good (ND5, ND4, COIII, COI), good (COII, cyt b), medium (ND3, ND2), poor (ND1, ATPase 6), and very poor (ND4L, ND6, ATPase 8). Although the universality of this observation was unclear, analysis of successive concatenation of the 13 protein-coding genes in the same ranking order revealed that the combined data sets comprising nucleotide sequences from the several top-ranked protein-coding genes (no 3rd codon positions) plus the 22 concatenated tRNA genes (stem regions only) best recovered the expected phylogeny, with all internal branches being supported by bootstrap values >90%. We conclude that judicious choice of mitochondrial genes and appropriate data weighting, in conjunction with purposeful taxonomic sampling, are prerequisites for resolving higher-level relationships in teleosts under the maximum-parsimony optimality criterion.

Animals↗

PubMedDate not supplied

Molecular cloning and characterization of a superoxide dismutase (sod) gene in Pneumocystis carinii.

This work reports the isolation and characterization of a gene encoding a superoxide dismutase (SOD, EC.1.15.1.1.) from Pneumocystis carinii derived from rat. Sense and antisense oligonucleotides, deduced from SOD amino acid sequences from a wide variety of organisms, allowed amplification of a 669 bp genomic DNA fragment specific to this P. carinii. RACE-PCR was used to obtain the major part of the complementary DNA; the 5'- and 3'-genomic regions were obtained respectively from a Mbol subgenomic library and from an amplified fragment using oligonucleotides designed from the cDNA sequence. Comparison of genomic and cDNA sequences showed an open reading frame of 660 bp interrupted by seven small introns. The deduced amino acid sequence contained 220 residues. Protein sequence alignment demonstrated the highest homology (50.5% identity; 70.3% similarity) with Saccharomyces cerevisiae manganese-SOD (MnSOD) suggesting that P. carinii SOD belongs to the mitochondrial MnSOD group. A putative targeting peptide found at the 5'-end of the P. carinii SOD sequence also suggested its mitochondrial localization.

Amino Acid Sequence↗

PubMed2005

Characterization of a mitochondrially targeted single-stranded DNA-binding protein in Arabidopsis thaliana.

A gene encoding a predicted mitochondrially targeted single-stranded DNA binding protein (mtSSB) was identified in the Arabidopsis thaliana genome sequence. This gene (At4g11060) codes for a protein of 201 amino acids, including a 28-residue putative mitochondrial targeting transit peptide. Protein sequence alignment shows high similarity between the mtSSB protein and single-stranded DNA binding proteins (SSB) from bacteria, including residues conserved for SSB function. Phylogenetic analysis indicates a close relationship between this protein and other mitochondrially targeted SSB proteins. The predicted targeting sequence was fused with the GFP coding region, and the organellar localization of the expressed fusion protein was determined. Specific targeting to mitochondria was observed in in-vitro import experiments and by transient expression of a GFP fusion construct in Arabidopsis leaves after microprojectile bombardment. The mature mtSSB coding region was overexpressed in Escherichia coli and the protein was purified for biochemical characterization. The purified protein binds single-stranded, but not double-stranded, DNA. MtSSB stimulates the homologous strand-exchange activity of E. coli RecA. These results indicate that mtSSB is a functional homologue of the E. coli SSB, and that it may play a role in mitochondrial DNA recombination.

Amino Acid Sequence↗

PubMed2004

Identification of two new HLA alleles: B3546 and B5611. How reliable are the published HLA-B intron 2 sequences?

Polymerase chain reaction-sequence-specific primer (PCR-SSP) typing for human leukocyte antigen (HLA)-B in a male 25-year-old Caucasian individual of Iranian origin and in a 42-year-old German Caucasian bone marrow donor revealed reaction patterns that did not agree with any known HLA specificity, thus suggesting in both cases the existence of a novel allele. Sequence-based typing (SBT) after allelic separation revealed the sequences of the new alleles HLA-B*5611 and B*3546. The sequence patterns of both new alleles might have been generated as the results of double crossing over, possibly over several generations. During the analysis of the HLA-B*3546 intron 2 sequence for possible crossing over points, a base insert, an additional G after position 700, was found. This insert was analyzed using SBT and PCR-SSP and was found to be present not only in all samples carrying B*35, but also in all HLA-B specificities tested. It appears that all known HLA-B alleles may contain a G insert at position 700 of intron 2, and that the published intron 2 sequence alignments of the HLA-B locus may contain errors at this position.

Adult↗

PubMed2004

Identification and analysis of a new family of bacterial serine proteinases.

A family of hypothetical proteins, identified predominantly from archaeal genomes, has been analyzed in order to understand its functional characteristics. Using extensive sequence similarity searches it is inferred that this family is remotely related (best sequence identity is 19%) to ClpP proteinases that belongs to serine proteinase class. This family of hypothetical proteins is referred to as SDH proteinase family based on conserved sequential order of Ser, Asp and His residues and predicted serine proteinase activity. Results of fold recognition of SDH family sequences confirmed the remote relationship between SDH proteinases and Clp proteinases and revealed similar tertiary location of putative catalytic triad residues critical for serine proteinase function. However, the best sequence alignment we could obtain suggests that while catalytic Ser is conserved across Clp and SDH proteinases the location of the other catalytic triad residues, namely, His and Asp are swapped in their amino acid alignment positions and hence in 3-D structure. The evidence of conserved catalytic triad suggests that SDH could be a new family of serine proteinases with the fold of Clp proteinase, however sharing the catalytic triad order of carboxypeptidase clan. Signal peptide sequence identified at the N-terminus of some of the homologues suggests that these might be secretory serine proteinases involved in cleavage of extracellular proteins while the remote homologues, ClpP proteinases, are known to work in intracellular environment.

Amino Acid Sequence↗

PubMed2004

GASP: Gapped Ancestral Sequence Prediction for proteins.

BACKGROUND: The prediction of ancestral protein sequences from multiple sequence alignments is useful for many bioinformatics analyses. Predicting ancestral sequences is not a simple procedure and relies on accurate alignments and phylogenies. Several algorithms exist based on Maximum Parsimony or Maximum Likelihood methods but many current implementations are unable to process residues with gaps, which may represent insertion/deletion (indel) events or sequence fragments. RESULTS: Here we present a new algorithm, GASP (Gapped Ancestral Sequence Prediction), for predicting ancestral sequences from phylogenetic trees and the corresponding multiple sequence alignments. Alignments may be of any size and contain gaps. GASP first assigns the positions of gaps in the phylogeny before using a likelihood-based approach centred on amino acid substitution matrices to assign ancestral amino acids. Important outgroup information is used by first working down from the tips of the tree to the root, using descendant data only to assign probabilities, and then working back up from the root to the tips using descendant and outgroup data to make predictions. GASP was tested on a number of simulated datasets based on real phylogenies. Prediction accuracy for ungapped data was similar to three alternative algorithms tested, with GASP performing better in some cases and worse in others. Adding simple insertions and deletions to the simulated data did not have a detrimental effect on GASP accuracy. CONCLUSIONS: GASP (Gapped Ancestral Sequence Prediction) will predict ancestral sequences from multiple protein alignments of any size. Although not as accurate in all cases as some of the more sophisticated maximum likelihood approaches, it can process a wide range of input phylogenies and will predict ancestral sequences for gapped and ungapped residues alike.

Amino Acid Sequence↗

PubMed2004

Vector NTI, a balanced all-in-one sequence analysis suite.

Vector NTI is a well-balanced desktop application integrated for molecular sequence analysis and biological data management. It has a centralised database and five application modules: Vector NTI, AlignX, BioAnnotator, ContigExpress and GenomBench. In this review, the features and functions available in this software are examined. These include database management, primer design, virtual cloning, alignments, sequence assembly, 3D molecular viewer and internet tools. Some problems encountered when using this software are also discussed. It is hoped that this review will introduce this software to more molecular biologists so they can make better-informed decisions when choosing computational tools to facilitate their everyday laboratory work. This tool can save time and enhance analysis but it requires some learning on the user's part and there are some issues that need to be addressed by the developer.

Algorithms↗

PubMed2003

Deciphering a novel thioredoxin-like fold family.

Sequence--and structure-based searching strategies have proven useful in the identification of remote homologs and have facilitated both structural and functional predictions of many uncharacterized protein families. We implement these strategies to predict the structure of and to classify a previously uncharacterized cluster of orthologs (COG3019) in the thioredoxin-like fold superfamily. The results of each searching method indicate that thioltransferases are the closest structural family to COG3019. We substantiate this conclusion using the ab initio structure prediction method rosetta, which generates a thioredoxin-like fold similar to that of the glutaredoxin-like thioltransferase (NrdH) for a COG3019 target sequence. This structural model contains the thiol-redox functional motif CYS-X-X-CYS in close proximity to other absolutely conserved COG3019 residues, defining a novel thioredoxin-like active site that potentially binds metal ions. Finally, the rosetta-derived model structure assists us in assembling a global multiple-sequence alignment of COG3019 with two other thioredoxin-like fold families, the thioltransferases and the bacterial arsenate reductases (ArsC).

Amino Acid Sequence↗

PubMed1997

Evolution of phosphagen kinase. VI. Isolation, characterization and cDNA-derived amino acid sequence of lombricine kinase from the earthworm Eisenia foetida, and identification of a possible candidate for the guanidine substrate recognition site.

Lombricine kinase (LK) from the body wall muscle of the earthworm Eisenia foetida was purified to homogeneity. The enzyme was shown to be a dimer consisting of 40 kDa subunits. The cDNA-derived amino acid sequence of 370 residues of Eisenia LK was determined. The validity of the sequence was supported by chemical sequencing of internal tryptic peptides. This is the first reported lombricine kinase amino acid sequence. Alignment of Eisenia LK with those of creatine kinases (CKs), arginine kinases (AKs) and glycocyamine kinase (GK) suggested a region displaying remarkable amino acid deletions (referred to GS region), as a possible candidate for guanidine substrate recognition site. A phylogenetic analysis using amino acid sequences of all four phosphagen kinases indicates that CK, GK and LK probably evolved from a common immediate ancestor protein.

Amino Acid Sequence↗

PubMed2005

The 1.1-angstrom structure of the spindle checkpoint protein Bub3p reveals functional regions.

Bub3p is a protein that mediates the spindle checkpoint, a signaling pathway that ensures correct chromosome segregation in organisms ranging from yeast to mammals. It is known to function by co-localizing at least two other proteins, Mad3p and the protein kinase Bub1p, to the kinetochore of chromosomes that are not properly attached to mitotic spindles, ultimately resulting in cell cycle arrest. Prior sequence analysis suggested that Bub3p was composed of three or four WD repeats (also known as WD40 and beta-transducin repeats), short sequence motifs appearing in clusters of 4-16 found in many hundreds of eukaryotic proteins that fold into four-stranded blade-like sheets. We have determined the crystal structure of Bub3p from Saccharomyces cerevisiae at 1.1 angstrom and a crystallographic R-factor of 15.3%, revealing seven authentic repeats. In light of this, it appears that many of these repeats therefore remain hidden in sequences of other proteins. Analysis of random and site-directed mutants identifies the surface of Bub3p involved in checkpoint function through binding of Bub1p and Mad3p. Sequence alignments indicate that these surfaces are mostly conserved across Bub3 proteins from diverse species. A structural comparison with other proteins containing WD repeats suggests that these folds may bind partner proteins using similar surface areas on the top and sides of the propeller. The sequences composing these regions are the most divergent within the repeat across all WD repeat proteins and could potentially be modulated to provide specificity in partner protein binding without perturbation of the core structure.

Amino Acid Sequence↗

PubMed2005

Prediction of functional class of novel plant proteins by a statistical learning method.

In plant genomes, the function of a substantial percentage of the putative protein-coding open reading frames (ORFs) is unknown. These ORFs have no significant sequence similarity to known proteins, which complicates the task of functional study of these proteins. Efforts are being made to explore methods that are complementary to, or may be used in combination with, sequence alignment and clustering methods. A web-based protein functional class prediction software, SVMProt, has shown some capability for predicting functional class of distantly related proteins. Here the usefulness of SVMProt for functional study of novel plant proteins is evaluated. To test SVMProt, 49 plant proteins (without a sequence homolog in the Swiss-Prot protein database, not in the SVMProt training set, and with functional indications provided in the literature) were selected from a comprehensive search of MEDLINE abstracts and Swiss-Prot databases in 1999-2004. These represent unique proteins the function of which, at present, cannot be confidently predicted by sequence alignment and clustering methods. The predicted functional class of 31 proteins was consistent, and that of four other proteins was weakly consistent, with published functions. Overall, the functional class of 71.4% of these proteins was consistent, or weakly consistent, with functional indications described in the literature. SVMProt shows a certain level of ability to provide useful hints about the functions of novel plant proteins with no similarity to known proteins.

Artificial Intelligence↗

PubMed2001

Structural, biochemical and electrostatic basis of serotype specificity in bovine enteroviruses.

We have performed immunostructural analyses of three closely related picornaviruses in order to gain understanding of the biochemical and structural basis of serotype specificity. We carried out sequence alignments of the capsid regions of three bovine enterovirus strains: VG-5-27 and M-4 from serotype 1 and PS-87 from serotype 2. Using our knowledge of the three dimensional and antigenic structure of strain VG-5-27 and the high levels of sequence identity between the strains, we have calculated the structures and solvent-accessible electrostatic potentials of the epitopes of all three viruses. We have demonstrated the viability of the molecular models of the epitopes of the M-4 and PS-87 strains. In each of the strains, we have explained the serotype specificities in terms of specific physical and chemical properties, and identified individual residues which are pivotal in determination of antibody recognition. These changes are in agreement with the known cross-reactivity of peptide and antiviral sera, showing that it is possible to derive structures for short variable sections of proteins of high sequence identity using molecular modelling which are significant in terms of biological function. We believe this study to be a novel approach in the analysis of virus serotype specificity.

Amino Acid Sequence↗

PubMed1996

Identification of the prolyl isomerase domain of Escherichia coli trigger factor.

E. coli trigger factor is a protein of 48 kDa which was recently identified as a ribosome-bound peptidyl-prolyl-cis/transisomerase (PPIase) capable of catalysing protein folding in vitro. We found trigger factor in association with nascent polypeptide chains, suggesting a function in the co-translational folding of proteins. Sequence comparisons revealed a homology of a segment of trigger factor with PPIases of the FK506 binding protein (FKBP) family. Here, we report on the purification of trigger factor and a domain assignment of its polypeptide chain by microsequencing and mass spectroscopy of proteolytic fragments. Two proteases generated fragments of 12-13 kDa molecular weight that encompass the predicted FKBP domain and possess PPIase activity in vitro. Sequence alignment of the known trigger factor proteins demonstrates a high degree of conservation within this central functional domain of the protein.

Amino Acid Isomerases↗

PubMed1998

Phylogenetic utility of the nuclear gene arginine decarboxylase: an example from Brassicaceae.

Arginine decarboxylase (ADC) is an important enzyme in the production of putrescine and polyamines in plants. It is encoded by a single or low-copy nuclear gene that lacks introns in sequences studied to date. The rate of Adc amino acid sequence evolution is similar to that of ndhF for the angiosperm family studied. Highly conserved regions provide several target sites for PCR priming and sequencing and aid in nucleotide and amino acid sequence alignment across a range of taxonomic levels, while a variable region provides an increased number of potentially informative characters relative to ndhF for the taxa surveyed. The utility of the Adc gene in plant molecular systematic studies is demonstrated by analysis of its partial nucleotide sequences obtained from 13 representatives of Brassicaceae and 3 outgroup taxa, 2 from the mustard oil clade (order Capparales) and 1 from the related order Malvales. Two copies of the Adc gene, Adc1 and Adc2, are found in all members of the Brassicaceae studied to data except the basal genus Aethionema. The resulting Adc gene tree provides robust phylogenetic data regarding relationships within the complex mustard family, as well as independent support for proposed tribal realignments based on other molecular data sets such as those from chloroplast DNA.

Arabidopsis↗