PubMed Health⌕ Search

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

FoldMiner: structural motif discovery using an improved superposition algorithm.

We report an unsupervised structural motif discovery algorithm, FoldMiner, which is able to detect global and local motifs in a database of proteins without the need for multiple structure or sequence alignments and without relying on prior classification of proteins into families. Motifs, which are discovered from pairwise superpositions of a query structure to a database of targets, are described probabilistically in terms of the conservation of each secondary structure element's position and are used to improve detection of distant structural relationships. During each iteration of the algorithm, the motif is defined from the current set of homologs and is used both to recruit additional homologous structures and to discard false positives. FoldMiner thus achieves high specificity and sensitivity by distinguishing between homologous and nonhomologous structures by the regions of the query to which they align. We find that when two proteins of the same fold are aligned, highly conserved secondary structure elements in one protein tend to align to highly conserved elements in the second protein, suggesting that FoldMiner consistently identifies the same motif in members of a fold. Structural alignments are performed by an improved superposition algorithm, LOCK 2, which detects distant structural relationships by placing increased emphasis on the alignment of secondary structure elements. LOCK 2 obeys several properties essential in automated analysis of protein structure: It is symmetric, its alignments of secondary structure elements are transitive, its alignments of residues display a high degree of transitivity, and its scoring system is empirically found to behave as a metric.

Algorithms↗

Enhanced genome annotation using structural profiles in the program 3D-PSSM.

A method (three-dimensional position-specific scoring matrix, 3D-PSSM) to recognise remote protein sequence homologues is described. The method combines the power of multiple sequence profiles with knowledge of protein structure to provide enhanced recognition and thus functional assignment of newly sequenced genomes. The method uses structural alignments of homologous proteins of similar three-dimensional structure in the structural classification of proteins (SCOP) database to obtain a structural equivalence of residues. These equivalences are used to extend multiply aligned sequences obtained by standard sequence searches. The resulting large superfamily-based multiple alignment is converted into a PSSM. Combined with secondary structure matching and solvation potentials, 3D-PSSM can recognise structural and functional relationships beyond state-of-the-art sequence methods. In a cross-validated benchmark on 136 homologous relationships unambiguously undetectable by position-specific iterated basic local alignment search tool (PSI-Blast), 3D-PSSM can confidently assign 18 %. The method was applied to the remaining unassigned regions of the Mycoplasma genitalium genome and an additional 13 regions were assigned with 95 % confidence. 3D-PSSM is available to the community as a web server: http://www.bmm.icnet.uk/servers/3dpssm

Algorithms↗

Physical mapping of the bovine immunoglobulin heavy chain constant region gene locus.

Bovine antibodies have recently attracted increasing attention, as they have been shown to exhibit prophylactic and therapeutic properties in selected infectious diseases in humans. In the present study, we have isolated bacterial artificial chromosomes and cosmid clones containing the bovine JH, mu, delta, gamma 1, gamma 2, gamma 3, epsilon, and alpha genes, which allowed us to make a contig of the genes within the bovine IGHC locus. The genes are arranged in a 5'-JH-7 kb-mu-5 kb-delta-33 kb-gamma 3-20 kb-gamma 1-34 kb-gamma 2-20 kb-epsilon- 13 kb-alpha-3' order, spanning approximately 150 kb DNA. Examination of the bovine germline JH locus revealed six JH segments, two of which, JH1 and JH2, were shown to be functional although there was a strong preference for expression of the former. Sequence alignment of the bovine 5' E mu enhancer core region with those of other mammals, demonstrated an absence of the mu E3 motif and a shortened spacer between the mu A and mu B sites within the bovine E mu enhancer core region. Furthermore, the essential sequence element for class switching, switch mu, spanning approximately 3-kb repetitive sequence and abundant in the switch region motifs CTGGG (187 repeats) and CTGAG (127 repeats), was identified immediately upstream of the mu gene. A further sequence comparison revealed that the bovine IGHC genes display an extensive polymorphism leading to expression of multiple antibody allotypes.

Amino Acid Sequence↗

Analysis in non-human primates reveals that the ancestral Band 3 gene encodes Dib and the Band 3-Memphis phenotype.

BACKGROUND: The anion exchanger, Band 3, carries antigens in the Diego blood group system, and can carry the Band 3-Memphis phenotype. Although Di(b) is of high prevalence and Band 3-Memphis is of low prevalence in humans, it has been suggested that both are on the ancestral gene. We determined the orthologue nucleotide sequences corresponding to these two polymorphic sites, Di(a)/Di(b) (2561T > C; Leu854Pro) and Band 3-Memphis(166A > G; Lys56Glu) in several nonhuman primates. METHODS: Genomic DNA was extracted from blood samples of great apes, lesser apes, old world monkeys, new world monkeys and prosimians. PCR amplifications were done with primer pairs that were located in the flanking intronic regions of Exon 4 and Exon 19; and the amplified products were sequenced. RESULTS: Amino acid sequence alignment of nonhuman primates band 3 with that of human showed extensive homologies. In exon 4, Glu56Lys polymorphic site showed Glu similar to Band 3-Memphis type and in exon 19, Leu854Pro polymorphic site showed Pro indicating Di(b) phenotype. CONCLUSIONS: The nonhuman primates have nucleotide sequences of Di(b)(2561C) in cis to Band 3-Memphis (166G), which is consistent with the assertion that the Di(b) and Band 3-Memphis phenotype represents the ancestral Band 3 gene.

Amino Acid Sequence↗

Characterization of an atypical superoxide dismutase from Sinorhizobium meliloti.

Sinorhizobium meliloti Rm5000 is an aerobic bacterium that can live free in the soil or in symbiosis with the roots of leguminous plants. A single detectable superoxide dismutase (SOD) was found in free-living growth conditions. The corresponding gene was isolated from a genomic library by using a sod fragment amplified by PCR from degenerate primers as a probe. The sodA gene was located in the chromosome. It is transcribed monocistronically and encodes a 200-amino-acid protein with a theoretical M(r) of 22,430 and pI of 5. 8. S. meliloti SOD complemented a deficient E. coli mutant, restoring aerobic growth of a sodA sodB recA strain, when the gene was expressed from the synthetic tac promoter but not from its own promoter. Amino acid sequence alignment showed great similarity with Fe-containing SODs (FeSODs), but the enzyme was not inactivated by H(2)O(2). The native enzyme was purified and found to be a dimeric protein, with a specific activity of 4,000 U/mg. Despite its Fe-type sequence, atomic absorption spectroscopy showed manganese to be the cofactor (0.75 mol of manganese and 0.24 mol of iron per mol of monomer). The apoenzyme was prepared from crude extracts of S. meliloti. Activity was restored by dialysis against either MnCl(2) or Fe(NH(4))(2)(SO(4))(2), demonstrating the cambialistic nature of the S. meliloti SOD. The recovered activity with manganese was sevenfold higher than with iron. Both reconstituted enzymes were resistant to H(2)O(2). Sequence comparison with 70 FeSODs and MnSODs indicates that S. meliloti SOD contains several atypical residues at specific sites that might account for the activation by manganese and resistance to H(2)O(2) of this unusual Fe-type SOD.

Amino Acid Sequence↗

Comparative Studies on Crystal Structures of Four Members of Cytochrome P450 Superfamily.

Secondary and steric structures and hydropathy plots of the 4 crystals of P450cam, P450terp, P450eryF and P450BM3 were compared to illustrate the structural conservation of cytochrome P450 superfamily proteins. Although sequence identities of four P450s are generally low (19%-26%), their topology is quite similar. All four structures have 13 alpha-helices and beta1-beta4 sheets in common. Four crystal structures were superimpossed by root-mean-square (RMS) fit of the prophyrin ring carbon atoms of prosthetic group heme to obtain the structure-based sequence alignment of four proteins. The RMS deviations of Calpha distance of each motifs were analyzed by hierarchical cluster analysis. The structural subsets were divided into four categories of structural conservation: the most conserved region, the less conserved region, the less variable region and the variable region. The first two groups (56.9 percentage of the aligned positions that have no gaps) include all the interior structures and active site residues. All four P450 proteins have the common hydrophobic and hydrophilic segments by hydropathy plots analyses. All the comparison results of P450 protein crystal structures provided the basis for structure-based sequence alignment of cytochrome P450 proteins.

Journal Article↗

P-cats: prediction of catalytic residues in proteins from their tertiary structures.

UNLABELLED: P-cats is a web server that predicts the catalytic residues in proteins from the atomic coordinates. P-cats receives a coordinate file of the tertiary structure and sends out analytical results via e-mail. The reply contains a summary and two URLs to allow the user to examine the conserved residues: one for interactive images of the prediction results and the other for a graphical view of the multiple sequence alignment. AVAILABILITY: P-cats is freely available at http://p-cats.hgc.jp/p-cats CONTACT: kino@ims.u-tokyo.ac.jp

Algorithms↗

Biochemical characterization and mechanism of action of a thermostable beta-glucosidase purified from Thermoascus aurantiacus.

An extracellular beta-glucosidase from Thermoascus aurantiacus was purified to homogeneity by DEAE-Sepharose, Ultrogel AcA 44 and Mono-P column chromatography. The enzyme was a homotrimer, with a monomer molecular mass of 120 kDa; only the trimer was optimally active at 80 degrees C and at pH 4.5. At 90 degrees C, the enzyme showed 70% of its optimal activity. It was stable at pH 5.2 and at temperatures up to 70 degrees C for 48 h, but stability decreased above 70 degrees C and at pH values above and below 5.0. The enzyme hydrolysed aryl and alkyl beta-d-glucosides and cello-oligosaccharides, and was specific for substrates with a beta-glycosidic linkage. The hydroxy groups at positions 2, 4 and 6 of a glucose residue at the non-reducing end of a disaccharide appeared to be essential for catalysis. The enzyme had the lowest K(m) towards p-nitrophenyl beta-d-glucoside (0.1137 mM) and the highest k(cat) towards cellobiose and beta,beta-trehalose (17052 min(-1)). It released one glucose unit at a time from the non-reducing end of cello-oligosaccharides, and the rate of hydrolysis decreased with an increase in chain length. Glucose and d-delta-gluconolactone inhibited the beta-glucosidase competitively, with K(i) values of 0.29 mM and 8.3 nM respectively, while methanol, ethanol and propan-2-ol activated the enzyme. The enzyme catalysed the synthesis of methyl, ethyl and propyl beta-d-glucosides in the presence of methanol, ethanol and propan-2-ol respectively with either glucose or cellobiose, although cellobiose was preferred. An acidic pH favoured hydrolysis and transglycosylation, but high concentrations of alcohols favoured the latter reaction. The stereochemistry of cellobiose hydrolysis revealed that beta-glucosidase from T. aurantiacus is a retaining glycosidase, while N-terminal amino acid sequence alignment indicated that it is a member of glycoside hydrolase family 3.

Alcohols↗

Complete structure determination of N-acetyl-D-galactosamine-binding mistletoe lectin-3 from Viscum album L. album.

The primary structure of the B chain of the N-acetyl-D-galactosamine-recognizing mistletoe lectin-3 (ML-3B) has been deduced from proteolytic digest peptides of the purified glycoprotein, their HPLC-separation and Edman degradation and confirmation of the peptide sequences by MALDI-MS. ML-3B consists of 262 amino acid residues including 10 cysteine moieties. The structure and linkage of the carbohydrate side chains, connected to two N-glycosylation sites at positions Asn(95) and Asn(135) of the lectin, were determined by a combination of glycosidase treatment and MALDI-MS of corresponding glycopeptide fragments. The sequence alignment reveals a high homology with other B chains of type-II RIPs, although there are remarkable differences in the D-galactose-specific mistletoe lectin-1B chain. The recently published primary structure of the mistletoe lectin-3A chain1 and the now available primary sequence of the 3B chain allowed the construction of a preliminary homology model of ML-3. The model demonstrates, unequivocally, that ML-3 is a member of the type-II RIP family with rigid conservation of the enzymatic active site of the A chain and an identical overall protein fold. Specific amino acid residue exchanges and the different glycosylation pattern in comparison with ML-1 are discussed and related to the properties of the two glycoproteins. The knowledge of the complete primary structure of mistletoe lectin-3 is a major contribution towards more insight into the mechanism of the biological activity of commercial mistletoe preparations.

Acetylgalactosamine↗

Sequence and hydropathy profile analysis of two classes of secondary transporters.

A structural class in the MemGen classification of membrane proteins is a set of evolutionary related proteins sharing a similar global fold. A structural class contains both closely related pairs of proteins for which homology is clear from sequence comparison and very distantly related pairs, for which it is not possible to establish homology based on sequence similarity alone. In the latter case the evolutionary link is based on hydropathy profile analysis. Here, we use these evolutionary related sets of proteins to analyze the relationship between E-values in BLAST searches, sequence similarities in multiple sequence alignments and structural similarities in hydropathy profile analyses. Two structural classes of secondary transporters termed ST[3], which includes the Ion Transporter (IT) superfamily and ST[4], which includes the DAACS family (TC# 2.A.23) were extracted from the NCBI protein database. ST[3] contains 2051 unique sequences distributed over 32 families and 59 subfamilies. ST[4] is a smaller class containing 399 unique sequences distributed over 2 families and 7 subfamilies. One subfamily in ST[4] contains a new class of binding protein dependent secondary transporters. Comparison of the averaged hydropathy profiles of the subfamilies in ST[3] and ST[4] revealed that the two classes represent different folds. Divergence of the sequences in ST[4] is much smaller than observed in ST[3], suggesting different constraints on the proteins during evolution. Analysis of the correlation between the evolutionary relationship of pairs of proteins in a class and the BLAST E-value revealed that: (i) the BLAST algorithm is unable to pick up the majority of the links between proteins in structural class ST[3], (ii) "low complexity filtering" and "composition based statistics" improve the specificity, but strongly reduce the sensitivity of BLAST searches for distantly related proteins, indicating that these filters are too stringent for the proteins analyzed, and (iii) the E-value cut-off, which may be used to evaluate evolutionary significance of a hit in a BLAST search is very different for the two structural classes of membrane proteins.

Algorithms↗

Two histidines are essential for the activity of the beta-galactosidase from Lactobacillus delbrückii subsp. bulgaricus.

Twelve tyrosines, ten glutamates, and five histidines were individually substituted for phenylalanine, glutamine, and asparagine, respectively, in the beta-galactosidase from Lactobacillus delbrückii subsp. bulgaricus. Only Y509, E464, H351, and H534 appear to be essential for the activity of the enzyme. The residues Y509 and E464 are homologous to Y503 and E461, respectively, of the Escherichia coli lacZ beta-galactosidase which have been shown to be necessary for its activity. Surprisingly, a number of amino acids that are highly conserved in the sequence alignments of nine beta-galactosidases and four beta-glucuronidases are not essential for enzyme activity.

Amino Acid Sequence↗

Model structure for the human blood coagulation agent beta-factor XIIa.

An improvement to the human blood coagulation agent beta-factor XIIa three-dimensional model is proposed. The sequence alignment as well as the modeling procedures are presented and the minimized energy of the new model is reported before and after solvation of the active center.

Amino Acid Sequence↗

TreeDet: a web server to explore sequence space.

The TreeDet (Tree Determinant) Server is the first release of a system designed to integrate results from methods that predict functional sites in protein families. These methods take into account the relation between sequence conservation and evolutionary importance. TreeDet fully analyses the space of protein sequences in either user-uploaded or automatically generated multiple sequence alignments. The methods implemented in the server represent three main classes of methods for the detection of family-dependent conserved positions, a tree-based method, a correlation based method and a method that employs a principal component analyses coupled to a cluster algorithm. An additional method is provided to highlight the reliability of the position in the alignments. The server is available at http://www.pdg.cnb.uam.es/servers/treedet.

Amino Acid Sequence↗

Quantitative estimates of sequence divergence for comparative analyses of mammalian genomes.

Comparative sequence analyses on a collection of carefully chosen mammalian genomes could facilitate identification of functional elements within the human genome and allow quantification of evolutionary constraint at the single nucleotide level. High-resolution quantification would be informative for determining the distribution of important positions within functional elements and for evaluating the relative importance of nucleotide sites that carry single nucleotide polymorphisms (SNPs). Because the level of resolution in comparative sequence analyses is a direct function of sequence diversity, we propose that the information content of a candidate mammalian genome be defined as the sequence divergence it would add relative to already-sequenced genomes. We show that reliable estimates of genomic sequence divergence can be obtained from small genomic regions. On the basis of a multiple sequence alignment of approximately 1.4 megabases each from eight mammals, we generate such estimates for five unsequenced mammals. Estimates of the neutral divergence in these data suggest that a small number of diverse mammalian genomes in addition to human, mouse, and rat would allow single nucleotide resolution in comparative sequence analyses.

Animals↗

A molecular method for a qualitative analysis of potentially coding sequences of DNA.

Total sequence phylogenies have low information content. Ordinary misconceptions are that character quality can be ignored and that relying on computer algorithms is enough. Despite widespread preference for a posteriori methods of character evaluation, a priori methods are necessary to produce transformation series that are independent of tree topologies. We propose a stepwise qualitative method for analyzing protein sequences. Informative codons are selected, alternative amino acid transformation series are analyzed, and most parsimonious transformations are hypothesized. We conduct four phylogenetic analyses of philodryanine snakes. The tree based on all nucleotides produces least resolution. Trees based on the exclusion of third positions, on an asymmetric step matrix, and on our protocol, produce similar results. Our method eliminates noise by hypothesizing explicit transformation series for each informative protein-coding amino acid. This approaches qualitative methods for morphological data, in which only characters successfully interpreted in a phylogenetic context are used in cladistic analyses. The method allows utilizing character information contained in the original sequence alignment and, therefore, has higher resolution in inferring a phylogenetic tree than some traditional methods (such as distance methods).

Algorithms↗

Inferring property selection pressure from positional residue conservation.

In this study, we attempt to understand and explain positional selection pressure in terms of underlying physical and chemical properties. We propose a set of constraining assumptions about how these pressures behave, then describe a procedure for analysing and explaining the distribution of residues at a particular position in a multiple sequence alignment. In contrast to previous approaches, our model takes into account both amino acid frequencies and a large number of physical-chemical properties. By analysing each property separately, it is possible to identify positions where distinct conservation patterns are present. In addition, the model can easily incorporate sequence weights that adjust for bias in the sample sequences. Finally, a test of statistical significance is provided for our conservation measure. The applicability of this method is demonstrated on two HIV-1 proteins: Nef and Env. The tools, data and results presented in this article are available at http://flan.blm.cs.cmu.edu.

Algorithms↗

Evolutionary conservation in protein folding kinetics.

The sequence and structural conservation of folding transition states have been predicted on theoretical grounds. Using homologous sequence alignments of proteins previously characterized via coupled mutagenesis/kinetics studies, we tested these predictions experimentally. Only one of the six appropriately characterized proteins exhibits a statistically significant correlation between residues' roles in transition state structure and their evolutionary conservation. However, a significant correlation is observed between the contributions of individual sequence positions to the transition state structure across a set of homologous proteins. Thus the structure of the folding transition state ensemble appears to be more highly conserved than the specific interactions that stabilize it.

Animals↗

Molecular modeling of family GH16 glycoside hydrolases: potential roles for xyloglucan transglucosylases/hydrolases in cell wall modification in the poaceae.

Family GH16 glycoside hydrolases can be assigned to five subgroups according to their substrate specificities, including xyloglucan transglucosylases/hydrolases (XTHs), (1,3)-beta-galactanases, (1,4)-beta-galactanases/kappa-carrageenases, "nonspecific" (1,3/1,3;1,4)-beta-D-glucan endohydrolases, and (1,3;1,4)-beta-D-glucan endohydrolases. A structured family GH16 glycoside hydrolase database has been constructed (http://www.ghdb.uni-stuttgart.de) and provides multiple sequence alignments with functionally annotated amino acid residues and phylogenetic trees. The database has been used for homology modeling of seven glycoside hydrolases from the GH16 family with various substrate specificities, based on structural coordinates for (1,3;1,4)-beta-D-glucan endohydrolases and a kappa-carrageenase. In combination with multiple sequence alignments, the models predict the three-dimensional (3D) dispositions of amino acid residues in the substrate-binding and catalytic sites of XTHs and (1,3/1,3;1,4)-beta-d-glucan endohydrolases; there is no structural information available in the databases for the latter group of enzymes. Models of the XTHs, compared with the recently determined structure of a Populus tremulos x tremuloides XTH, reveal similarities with the active sites of family GH11 (1,4)-beta-D-xylan endohydrolases. From a biological viewpoint, the classification, molecular modeling and a new 3D structure of the P. tremulos x tremuloides XTH establish structural and evolutionary connections between XTHs, (1,3;1,4)-beta-D-glucan endohydrolases and xylan endohydrolases. These findings raise the possibility that XTHs from higher plants could be active not only on cell wall xyloglucans, but also on (1,3;1,4)-beta-D-glucans and arabinoxylans, which are major components of walls in grasses. A role for XTHs in (1,3;1,4)-beta-D-glucan and arabinoxylan modification would be consistent with the apparent overrepresentation of XTH sequences in cereal expressed sequence tags databases.

Amino Acid Sequence↗