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

PubMed2006

Evolution of EV71 genogroup in Taiwan from 1998 to 2005: an emerging of subgenogroup C4 of EV71.

In Taiwan, enterovirus 71 (EV71) has played an important role in severe enterovirus-related cases every year since the devastating outbreak in 1998. Three genogroups A, B, C occur worldwide; with the B and C genogroups being subdivided into B1-B4 and C1-C4 subgenogroups respectively. To understand the mutation of the EV71 genogroup in Taiwan before and after 1998, a total of 54 worldwide strains were studied including 41 Taiwanese strains obtained in 1986 and 1998-2004. A fragment of 207 bp of the VP4 region was amplified and sequenced. Genetic analysis was performed using MEGA software (version 3.0) for the nucleotide sequence alignment and phylogenetic analysis. In Taiwan, the subgenogroup B1 was predominant before 1998 while subgenogroup C2 was the major etiologic group in 1998 outbreak. A minor etiologic group outbreak in 1998, subgenogroup B4, became predominant during the period from 1999 to 2003. In this study, subgenogroup C4 emerged and became predominant in 2004 in Taiwan. The nucleotide differences between B1 and C2, C2 and B4, B4 and C4 were 20%-26%, 19%-27%, 18%-22%, respectively. Nucleotide sequence alignment revealed 67 substitutions. Most of the substitutions (62/67) were silent mutations. This is the first report about the emergence of EV71 subgenogroup C4 in Taiwan.

Base Sequence↗

PubMed2002

Protein sequence threading: Averaging over structures.

Multiple sequence alignments are a routine tool in protein fold recognition, but multiple structure alignments are computationally less cooperative. This work describes a method for protein sequence threading and sequence-to-structure alignments that uses multiple aligned structures, the aim being to improve models from protein threading calculations. Sequences are aligned into a field due to corresponding sites in homologous proteins. On the basis of a test set of more than 570 protein pairs, the procedure does improve alignment quality, although no more than averaging over sequences. For the force field tested, the benefit of structure averaging is smaller than that of adding sequence similarity terms or a contribution from secondary structure predictions. Although there is a significant improvement in the quality of sequence-to-structure alignments, this does not directly translate to an immediate improvement in fold recognition capability.

Animals↗

PubMed1993

Phyletic relationships of protein structures based on spatial preference of residues.

A structure-based scoring matrix MDPRE was derived from amino acid spatial preferences in protein structures. Sequence alignment and evolutionary studies by using MDPRE matrix gave similar results as those from ordinary sequence and structure alignments. It is interesting that a matrix derived from structure data solely could give comparable alignment results, strongly indicating the intimate connection between protein sequences and structures. The branch order and length from this approach were close to those obtained by a structure comparison method. Thus, by applying this structure-based matrix, the trees obtained should reflect evolutionary characteristics of protein structure. This approach takes advantage over a direct structure comparison in that (1) only a sequence and MDPRE matrix are needed, making it simple and widely applicable (especially in the absence of 3-dimensional protein structure data); (2) an established algorithm for sequence alignment and tree building could be employed, providing opportunities for direct comparison between matrices from different methodologies. One of the most striking features of this method is its capability to detect protein structure homologies when the sequence identities are low. This was well reflected in the given examples of the alignment of dinucleotide-binding domains.

Amino Acid Sequence↗

PubMed2005

REVCOM: a robust Bayesian method for evolutionary rate estimation.

MOTIVATION: Evolutionary conservation estimated from a multiple sequence alignment is a powerful indicator of the functional significance of a residue and helps to predict active sites, ligand binding sites, and protein interaction interfaces. Many algorithms that calculate conservation work well, provided an accurate and balanced alignment is used. However, such a strong dependence on the alignment makes the results highly variable. We attempted to improve the conservation prediction algorithm by making it more robust and less sensitive to (1) local alignment errors, (2) overrepresentation of sequences in some branches and (3) occasional presence of unrelated sequences. RESULTS: A novel method is presented for robust constrained Bayesian estimation of evolutionary rates that avoids overfitting independent rates and satisfies the above requirements. The method is evaluated and compared with an entropy-based conservation measure on a set of 1494 protein interfaces. We demonstrated that approximately 62% of the analyzed protein interfaces are more conserved than the remaining surface at the 5% significance level. A consistent method to incorporate alignment reliability is proposed and demonstrated to reduce arbitrary variation of calculated rates upon inclusion of distantly related or unrelated sequences into the alignment.

Algorithms↗

PubMed2002

The Pfam protein families database.

Pfam is a large collection of protein multiple sequence alignments and profile hidden Markov models. Pfam is available on the World Wide Web in the UK at http://www.sanger.ac.uk/Software/Pfam/, in Sweden at http://www.cgb.ki.se/Pfam/, in France at http://pfam.jouy.inra.fr/ and in the US at http://pfam.wustl.edu/. The latest version (6.6) of Pfam contains 3071 families, which match 69% of proteins in SWISS-PROT 39 and TrEMBL 14. Structural data, where available, have been utilised to ensure that Pfam families correspond with structural domains, and to improve domain-based annotation. Predictions of non-domain regions are now also included. In addition to secondary structure, Pfam multiple sequence alignments now contain active site residue mark-up. New search tools, including taxonomy search and domain query, greatly add to the functionality and usability of the Pfam resource.

Animals↗

PubMed1997

Inferring the conformation of RNA base pairs and triples from patterns of sequence variation.

The success of comparative analysis in resolving RNA secondary structure and numerous tertiary interactions relies on the presence of base covariations. Although the majority of base covariations in aligned sequences is associated to Watson-Crick base pairs, many involve non-canonical or restricted base pair exchanges (e.g. only G:C/A:U), reflecting more specific structural constraints. We have developed a computer program that determines potential base pairing conformations for a given set of paired nucleotides in a sequence alignment. This program (ISOPAIR) assumes that the base pair conformation is maintained through sequence variation without significantly affecting the path of the sugar-phosphate backbone. ISOPAIR identifies such 'isomorphic' structures for any set of input base pair or base triple sequences. The program was applied to base pairs and triples with known structures and sequence exchanges. In several instances, isomorphic structures were correctly identified with ISOPAIR. Thus, ISOPAIR is useful when assessing non-canonical base pair conformations in comparative analysis. ISOPAIR applications are limited to those cases where unusual base pair exchanges indeed reflect a non-canonical conformation.

Base Composition↗

PubMed2005

A new method for identification of protein (sub)families in a set of proteins based on hydropathy distribution in proteins.

Structural similarity among proteins is reflected in the distribution of hydropathicity along the amino acids in the protein sequence. Similarities in the hydropathy distributions are obvious for homologous proteins within a protein family. They also were observed for proteins with related structures, even when sequence similarities were undetectable. Here we present a novel method that employs the hydropathy distribution in proteins for identification of (sub)families in a set of (homologous) proteins. We represent proteins as points in a generalized hydropathy space, represented by vectors of specifically defined features. The features are derived from hydropathy of the individual amino acids. Projection of this space onto principal axes reveals groups of proteins with related hydropathy distributions. The groups identified correspond well to families of structurally and functionally related proteins. We found that this method accurately identifies protein families in a set of proteins, or subfamilies in a set of homologous proteins. Our results show that protein families can be identified by the analysis of hydropathy distribution, without the need for sequence alignment.

Algorithms↗

PubMed1996

Compared chemical properties of dermonecrotic and lethal toxins from spiders of the genus Loxosceles (Araneae).

Loxosceles spider venom usually causes a typical dermonecrotic lesion in bitten patients, but it may also cause systemic effects that may be lethal. Gel filtration on Sephadex G-100 of Loxosceles gaucho, L. laeta, or L. intermedia spider venoms resulted in three fractions (A, containing higher molecular mass components. B containing intermediate molecular mass components, and C with lower molecular mass components). The dermonecrotic and lethal activities were detected exclusively in fraction A of all three species. Analysis by SDS-PAGE showed that the major protein contained in fraction A has molecular weight approximately 35 kDa in L. gaucho and L. intermedia, but 32 kDa in L. laeta venom. These toxins were isolated from venoms of L. gaucho, L. laeta, and L. intermedia by SDS-PAGE followed by blotting to PVDF membrane and sequencing. A database search showed a high level of identity between each toxin and a fragment of the L. reclusa (North American spider) toxin. A multiple sequence alignment of the Loxosceles toxins showed many common identical residues in their N-terminal sequences. Identities ranged from 50.0% (L. gaucho and L. reclusa) to 61.1% (L. intermedia and L. reclusa). The purified toxins were also submitted to capillary electrophoresis peptide mapping after in situ partial hydrolysis of the blotted samples. The results obtained suggest that L. intermedia protein is more similar to L. laeta toxin than L. gaucho toxin and revealed a smaller homology between L. intermedia and L. gaucho. Altogether these findings suggest that the toxins responsible for most important activities of venoms of Loxosceles species have a molecular mass of 32-35 kDa and are probably homologous proteins.

Amino Acid Sequence↗

PubMed2005

Confidence intervals of similarity values determined for cloned SSU rRNA genes from environmental samples.

The goal of this research was to investigate the influence of the error rate of sequence determination on the differentiation of cloned SSU rRNA gene sequences for assessment of community structure. SSU rRNA cloned sequences from groundwater samples that represent different bacterial divisions were sequenced multiple times with the same sequencing primer. From comparison of sequence alignments with unedited data, confidence intervals were obtained from both a 'double binomial' model of sequence comparison and by non-parametric methods. The results indicated that similarity values below 0.9946 are likely derived from dissimilar sequences at a confidence level of 0.95, and not sequencing errors. The results confirmed that screening by direct sequence determination could be reliably used to differentiate at the species level. However, given sequencing errors comparable to those seen in this study, sequences with similarities above 0.9946 should be treated as the same sequence if a 95% confidence is desired.

Base Sequence↗

PubMed2004

The Selenophosphate synthetase gene from Leishmania major.

Selenophosphate synthetase coding sequence was cloned from Leishmania major by RT-PCR amplification. The DNA sequence was found to have an open reading frame encoding protein with 398 amino acids and does not have in-frame UGA codon. The deduced amino acid sequence indicates that it has cysteine residue instead of selenocysteine at the active site of enzyme. Amino acid sequence alignment of Selenophosphate synthetase from parasite with the human enzyme showed approximately 45% homology. The sequences also indicated presence of conserved amino acid residues and motifs that are present in mammalian Selenophosphate synthetase. Southern analysis done with restriction enzyme digested genomic DNA and pulse filed separated chromosome suggests that L. major genome contain a single copy of Selenophosphate synthetase sequence. Expression analysis by Northern analysis and RT-PCR indicated Selenophosphate synthetase mRNA is present in promastigote and amastigote stages of parasite.

Amino Acid Sequence↗

PubMed2006

LGICdb: a manually curated sequence database after the genomes.

Ligand-gated ion channels form transmembrane ionic pores controlled by the binding of chemicals. The LGICdb aims to be a non-redundant, manually curated resource offering access to the large number of subunits composing extracellularly activated ligand-gated ion channels, such as nicotinic, ATP, GABA and glutamate ionotropic receptors. Composed of more than 500 human curated entries, the XML native database has been relocated in 2004 to the EBI. Its facilities have been enhanced with a new search system, customized multiple sequence alignments and manipulation of protein structures (http://www.ebi.ac.uk/compneur-srv/LGICdb/). Despite the vast improvement of general sequence resources, the LGICdb still provide sequences unavailable elsewhere.

Databases, Protein↗

PubMedDate not supplied

Isolation and characterization of somatolactin genes from two cold water marine teleosts, lumpfish (Cyclopterus lumpus) and halibut (Hippoglossus hippoglossus).

Somatolactin (SL) is a newly discovered pituitary hormone structurally related to both growth hormone and prolactin. We isolated the SL gene from two cold water marine teleosts, lumpfish (Cyclopterus lumpus) and halibut (Hippoglossus hippoglossus). The SL genes for lumpfish were isolated from a complementary DNA (cDNA) library by random gene tagging and sequencing. Two full-size clones of 1.2 (lfSL-6) and 1.6 (lfSL-9) kilobases (kb) were characterized. The SL genes for halibut were identified from its cDNA library using lumpfish SL gene as a probe. Two full-size SL clones of 1.1 and 1.6 kb were identified in halibut. All these clones encoded full-size SL (229 amino acids for lumpfish SL, 230 amino acids for halibut SL), including 7 cysteines and two and three potential glycosylation sites in halibut and lumpfish, respectively. The differences in length of these clones are attributed to their differences in 3' untranslated sequences. Although the halibut SL genes contain the consesus polyadenylation signal AATAAA, the lumpfish SL gene has the ATTAAA sequence in the larger clone (lfSL-9). None of the consensus or the alternative polyadenylation signals were found in the shorter clone (lfSL-6). Protein homology and DNA sequence alignments of SL genes from all five known marine fishes suggest that the SL gene sequence is highly conserved.

Amino Acid Sequence↗

PubMed2006

Prediction of cis/trans isomerization in proteins using PSI-BLAST profiles and secondary structure information.

BACKGROUND: The majority of peptide bonds in proteins are found to occur in the trans conformation. However, for proline residues, a considerable fraction of Prolyl peptide bonds adopt the cis form. Proline cis/trans isomerization is known to play a critical role in protein folding, splicing, cell signaling and transmembrane active transport. Accurate prediction of proline cis/trans isomerization in proteins would have many important applications towards the understanding of protein structure and function. RESULTS: In this paper, we propose a new approach to predict the proline cis/trans isomerization in proteins using support vector machine (SVM). The preliminary results indicated that using Radial Basis Function (RBF) kernels could lead to better prediction performance than that of polynomial and linear kernel functions. We used single sequence information of different local window sizes, amino acid compositions of different local sequences, multiple sequence alignment obtained from PSI-BLAST and the secondary structure information predicted by PSIPRED. We explored these different sequence encoding schemes in order to investigate their effects on the prediction performance. The training and testing of this approach was performed on a newly enlarged dataset of 2424 non-homologous proteins determined by X-Ray diffraction method using 5-fold cross-validation. Selecting the window size 11 provided the best performance for determining the proline cis/trans isomerization based on the single amino acid sequence. It was found that using multiple sequence alignments in the form of PSI-BLAST profiles could significantly improve the prediction performance, the prediction accuracy increased from 62.8% with single sequence to 69.8% and Matthews Correlation Coefficient (MCC) improved from 0.26 with single local sequence to 0.40. Furthermore, if coupled with the predicted secondary structure information by PSIPRED, our method yielded a prediction accuracy of 71.5% and MCC of 0.43, 9% and 0.17 higher than the accuracy achieved based on the singe sequence information, respectively. CONCLUSION: A new method has been developed to predict the proline cis/trans isomerization in proteins based on support vector machine, which used the single amino acid sequence with different local window sizes, the amino acid compositions of local sequence flanking centered proline residues, the position-specific scoring matrices (PSSMs) extracted by PSI-BLAST and the predicted secondary structures generated by PSIPRED. The successful application of SVM approach in this study reinforced that SVM is a powerful tool in predicting proline cis/trans isomerization in proteins and biological sequence analysis.

Databases, Protein↗

PubMed1994

On near-optimal alignments of biological sequences.

A near-optimal alignment between a pair of sequences is an alignment whose score lies within the neighborhood of the optimal score. We present an efficient method for representing all alignments whose score is within any given delta from the optimal score. The representation is a compact graph that makes it easy to impose additional biological constraints and select one desirable alignment from the large set of alignments. We study the combinatorial nature of near-optimal alignments, and define a set of "canonical" near-optimal alignments. We then show how to enumerate near-optimal alignments efficiently in order of their score, and count their number. When applied to comparisons of two distantly related proteins, near-optimal alignments reveal that the most conserved regions among the near-optimal alignments are the highly structured regions in the proteins. We also show that by counting the number of near optimal alignments as a function of the distance from the optimal score, we can select a good set of parameters that best constraints the biologically relevant alignments.

Algorithms↗

PubMed1998

Proximal and distal histidines in thyroid peroxidase: relation to the alternatively spliced form, TPO-2.

The distal and proximal histidines in thyroid peroxidase (TPO), located by amino acid sequence alignment with their known counterparts in myeloperoxidase, are His 239 and His 494, respectively. These histidines lie outside the 57 amino acid peptide (residues 533-589) that is absent in the alternatively spliced form, TPO-2. However, asparagine 579, which very likely forms a stabilizing hydrogen bond with the proximal histidine in TPO, lies within the missing peptide region. The absence of Asn 579 from TPO-2 may be at least partially responsible for the reported lack of activity of this form of the enzyme. Formation of TPO compound I may also depend on Arg 396, based on analogy with the catalytic mechanism previously proposed for the more widely studied plant and fungal peroxidases. A multiple sequence alignment prepared with five mammalian and five invertebrate peroxidases shows complete conservation of Arg 396, as well as residues corresponding to His 239, His 494, and Asn 579 in TPO. The animal peroxidases comprise a family of homologous proteins that differ markedly from the plant/fungal/bacterial peroxidases in primary, secondary, and tertiary structure, yet share with them a common function. Animal peroxidases probably arose independently of the plant/fungal/bacterial peroxidase superfamily and most likely belong to a different gene family. The relation between animal and nonanimal peroxidases may represent an example of convergent evolution to a common enzymatic mechanism.

Alternative Splicing↗

PubMed1977

Secondary-structure predictions of calcium-binding proteins.

The known tertiary structure of carp muscle parvalbumin is consistent with an "EF-hand" architecture (helix-loop-helix) for each calcium-ion binding site. Primary-sequence alignments have indicated four EF hands in rabbit skeletal muscle troponin C and in rabbit myosin alkali light chains. Five secondary-structure prediction methods, based on amino acid sequence only, have been fully computerized and used to calculate joint prediction histograms for several calcium-binding proteins. The joint histogram can suggest directly the extent and sequence of the helical- and loop-structural elements, as well as any secondary structural distortions or evolutionary developments. Since the histogram predicted well the length and sequence of secondary structural elements in carp muscle parvalbumin, it seemed reasonable to calculate the joint distribution for other proteins that might bind calcium through the EF-hand configuration. The histograms indicated the four EF-hand regions speculated fro rabbit skeletal muscle troponin C but suggested only three such hands in bovine cardiac muscle troponin C and with a distorted fourth hand. Considerable secondary structural distortion is postulated for the alkali light chains. Possible EF configurations consistent with the histogram results are speculated for Escherichia coli acyl-carrier protein and bovine prothrombin fragment 1, which have been shown to bind calcium. The secondary-structure-prediction algorithms appear to be a useful adjunct to sequence-alignment techniques, especially in cases where the primary sequence homology is weak or the evolutionary distance is large.

Amino Acid Sequence↗

PubMed2005

Mutational effects at the tetramerization site of nonerythroid alpha spectrin.

Spectrin, a prominent cytoskeletal protein, exerts its fundamental role in cellular function by forming a sub-membrane filamentous network. An essential aspect of spectrin network formation is the tetramerization of spectrin alphabeta heterodimers. We used laboratory methods, the yeast two-hybrid system and random mutagenesis, to investigate, for the first time, effects of amino acid mutations on tetramerization of nonerythroid (brain) spectrin (fodrin). Based on high sequence homology with erythroid spectrin, we assume the putative tetramerization region of nonerythroid alpha-spectrin at the N-terminal region. We introduced mutations in the region consisting of residues 1-45 and studied mutational effects on spectrin alphabeta association to form tetramers. We detected single, double, and triple mutations involving 24 residues in this region. These amino acid mutations of nonerythroid alpha-spectrin exhibit full, partial, or no effect on the association with nonerythroid beta-spectrin. Single amino acid mutations in the region of residues 1-9 (D2Y, G5V, V6D, and V8M) did not affect the association. However, seven single mutations (I15F, I15N, R18G, V22D, R25P, Y26N, and R28P) affected the alphabeta association. These mutations were clustered in the region predicted by sequence alignment to be crucial in nonerythroid alpha-spectrin for tetramerization, a region that spanned residues 12-36, corresponding to the partial domain Helix C' (residues 21-45) in erythroid alpha-spectrin. In addition, two other mutations, one upstream and one downstream of this region at positions 10 (E10D) and 37 (R37P), also affected the alphabeta association. Our results implied nonerythroid alpha-spectrin partial domain helix may be longer than Helix C' (residues 21-45 and a total of 25 residues) in erythroid alpha-spectrin and spanned at least residues 10-37. It is interesting to note that seven out of these nine single mutations (I15F, I15N, R18G, V22D, R25P, Y26N, R37P) were at the a, d, e or g heptad positions based on sequence alignment with erythroid alpha-spectrin. Four of the mutated residues (I15, R18, V22, R25) are conserved in both erythroid and nonerythroid spectrin. These positions were previously identified as hot spots in erythroid alpha-spectrin that lead to severe hematological symptoms. This study clearly demonstrated that single mutation in a region predicted to be critical functionally in nonerythroid alpha-spectrin indeed leads to functional abnormalities and may lead to neurological disorders.

Amino Acid Sequence↗

PubMed1975

The primary structure of human plasma high density apolipoprotein glutamine I (ApoA-I). II. The amino acid sequence and alignment of cyanogen bromide fragments IV, III, and I.

Apolipoprotein glutamine I (apoLP-Gln-I or apoA-I) is the major protein constituent of the human plasma high density lipoproteins. Cleavage of this protein with cyanogen bromide yields four fragments, designated in the order of elution from Bio-Gel P-30 as CNBr I, II, III, and IV. In the first paper in this series, the amino acid sequence of the NH2-terminal fragment, CNBr II, is reported. In the present study, CNBr IV, III, and I, containing, respectively, 25, 36, and 94 amino acids were sequenced by conventional means. To establish the alignment of the cyanogen bromide fragments, apoLP-Gln-I was digested with trypsin and two of the three methionine-containing tryptic peptides were isolated. The amino acid sequence of apoLP-Gln-I is as follows: (see article). With the complete amino acid sequence available, a CPK space-filling model of apoLP-Gln-I was constructed. The protein was placed into an alpha helical conformation wherever the primary structure permitted. Thirteen helical regions were identified. These regions account for 70% of the amino acid residues of the protein. Each helix is characterized as having two faces (amphipathic). One is a polar face that occupies approximately 180 degrees of the surface of the helix and the other is a hydrophobic face that occupies the other 180 degrees of the helical surface. Similar amphipathic helices have been identified previously in the other lipoprotein-proteins that have known sequences. It is suggested that the amphipathic helical regions of apoLP-Gln-I are important in the binding of phospholipids in high density lipoproteins.

Amino Acid Sequence↗