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

Structure and polymorphism of the Chironomus thummi gene encoding special lobe-specific silk protein, ssp160.

cDNA encoding Chironomus thummi ssp160 was used to isolate a genomic clone that hybridized in situ to band A2b on polytene chromosome IV, the site of the ssp160 gene. DNA sequencing, primer extension and gene/cDNA nucleotide sequence alignment revealed the gene contains six exons and five introns; 70% of ssp160 is encoded in exon 3. Variations between cDNA and gene sequences led to the design of a polymerase chain reaction, restriction fragment length polymorphism assay that was subsequently used to demonstrate the existence of polymorphic alleles whose distribution varied between geographically separated populations of larvae. The polymorphism is associated with codon deletions in a six-amino-acid repeat containing an N-linked glycosylation motif. These deletions may have resulted from slipped-strand mispairing during DNA replication.

Amino Acid Sequence↗

Automatic prediction of protein domains from sequence information using a hybrid learning system.

MOTIVATION: We describe a novel method for detecting the domain structure of a protein from sequence information alone. The method is based on analyzing multiple sequence alignments that are derived from a database search. Multiple measures are defined to quantify the domain information content of each position along the sequence and are combined into a single predictor using a neural network. The output is further smoothed and post-processed using a probabilistic model to predict the most likely transition positions between domains. RESULTS: The method was assessed using the domain definitions in SCOP and CATH for proteins of known structure and was compared with several other existing methods. Our method performs well both in terms of accuracy and sensitivity. It improves significantly over the best methods available, even some of the semi-manual ones, while being fully automatic. Our method can also be used to suggest and verify domain partitions based on structural data. A few examples of predicted domain definitions and alternative partitions, as suggested by our method, are also discussed. AVAILABILITY: An online domain-prediction server is available at http://biozon.org/tools/domains/

Algorithms↗

Primary structure of locust opsins: a speculative model which may account for ultraviolet wavelength light detection.

The sequences of two locust opsins have been determined by dideoxy nucleotide sequencing of PCR products from cDNA derived from eyecup tissue. The opsins (Lo1 and Lo2) are encoded by 381 and 380 amino acid residues, respectively, with hydropathy profiles and placement of key amino acid residues suggestive of a typical seven-transmembrane rhodopsin structure. The sequence alignment of Lo1 reveals significant homology to mantid opsin. These opsins contain retinal as their visual chromophore and have similarity to the Rh1 type sequences from Drosophila and Calliphora which use 3-hydroxy retinal. Lo2 is most closely related to the Rh3/4 type of visual pigments from Drosophila. The retinal-based opsins show reduced numbers of charged amino acids in the loop region connecting transmembrane segments V and VI compared to the 3-hydroxy retinal opsins. Sequence alignment of all the known insect visual pigments has shown that only those with maximal sensitivity in the blue/UV spectral range, Lo2 and the Rh3/4 opsins of Drosophila, have three charged amino acids in transmembrane segments II, IV and VII. The charged residue in transmembrane VII is two helical turns away from the positively charged Schiff base and could act directly as a counterion to it. From the secondary structure analysis of opsin, the two charged residues in transmembrane II and IV would be in close proximity to form a dipole. These polar motifs in Lo2 and Rh3/Rh4 could act in wavelength modulation of short wavelength sensitive pigments and substantiate the proposed external two-point charge model which accounts for the spectral sensitivity of visual pigments [Honig, B., Dinur, U., Nakanishi, K., Balogh-Nair, V., Gawinowicz, M.A. and Motto, M. (1979). Journal of the American Chemical Society, 101, 7084-7086].

Amino Acid Sequence↗

A symmetric-iterated multiple alignment of protein sequences.

A new symmetric-iterative method for multiple alignment of protein sequences is presented. The method can be described as a combination of motif finding and dynamic programming procedures. It uses each sequence as a standard to which all sequences are aligned based on the significant segment pair alignment (SSPA) protocol. Sequences are further matched using a reduced scoring threshold to provide fillers and extensions between highly significant segment pair matches. The method produces alignment blocks that accommodate indels and are separated by variable-length unaligned segments. Construction of consensus sequences is iterative, assigning greater weights to more distantly related sequences. A consensus sequence and various measures of conservation at each aligned position can be used for comparisons between protein families, for data base searches, and for analysis of functional and evolutionary features. The method is illustrated on the extended family of prokaryotic and eukaryotic RecA-like sequences. The RecA-like sequences reveal extended alignments among eubacterial RecA and separately among eukaryotic/archaebacterial Rad51/RadA. Eleven conserved blocks are common to both groups, two of them encompassing the ATP-binding A and B-sites. Among the most conserved positions are glycine residues. For example, they occur twice as doublets putatively serving as hinge connections that provide opportunity for alternative structural conformations. Also several charged/polar residues are highly conserved, probably consequent upon the extensive intermonomer interactions in RecA/Rad51 filament formation and possibly relevant protein-protein and protein-nucleic acid interactions.

Algorithms↗

Yet another numbering scheme for immunoglobulin variable domains: an automatic modeling and analysis tool.

A common residue numbering scheme for all immunoglobulin variable domains (immunoglobulin light chain lambda (V(lambda)) and kappa (V(kappa)) variable domains, heavy chain variable domains (V(H)) and T-cell receptor alpha (V(alpha)), beta (V(beta)), gamma (V(gamma)) and delta (V(delta)) variable domains) has been devised. Based on the spatial alignment of known three-dimensional structures of immunoglobulin domains, it places the alignment gaps in a way that minimizes the average deviation from the averaged structure of the aligned domains. This residue numbering scheme was applied to the immunoglobulin variable domain structures in the PDB database to automate the extraction of information on structural variations in homologous positions of the different molecules. A number of methods are presented that allow the automated projection of information derived from individual structures or from the comparison of multi-structure alignments onto a graphical representation of the sequence alignment.

Amino Acid Sequence↗

Crystal structure of Rv2118c: an AdoMet-dependent methyltransferase from Mycobacterium tuberculosis H37Rv.

Rv2118c belongs to the class of conserved hypothetical proteins from Mycobacterium tuberculosis H37Rv. The crystal structure of Rv2118c in complex with S-adenosyl-l-methionine (AdoMet) has been determined at 1.98 A resolution. The crystallographic asymmetric unit consists of a monomer, but symmetry-related subunits interact extensively, leading to a tetrameric structure. The structure of the monomer can be divided functionally into two domains: the larger catalytic C-terminal domain that binds the cofactor AdoMet and is involved in the transfer of methyl group from AdoMet to the substrate and a smaller N-terminal domain. The structure of the catalytic domain is very similar to that of other AdoMet-dependent methyltransferases. The N-terminal domain is primarily a beta-structure with a fold not found in other methyltransferases of known structure. Database searches reveal a conserved family of Rv2118c-like proteins from various organisms. Multiple sequence alignments show several regions of high sequence similarity (motifs) in this family of proteins. Structure analysis and homology to yeast Gcd14p suggest that Rv2118c could be an RNA methyltransferase, but further studies are required to establish its functional role conclusively.

Amino Acid Motifs↗

Response regulators of bacterial signal transduction systems: selective domain shuffling during evolution.

Response regulators of bacterial sensory transduction systems generally consist of receiver module domains covalently linked to effector domains. The effector domains include DNA binding and/or catalytic units that are regulated by sensor kinase-catalyzed aspartyl phosphorylation within their receiver modules. Most receiver modules are associated with three distinct families of DNA binding domains, but some are associated with other types of DNA binding domains, with methylated chemotaxis protein (MCP) demethylases, or with sensor kinases. A few exist as independent entities which regulate their target systems by noncovalent interactions. In this study the molecular phylogenies of the receiver modules and effector domains of 49 fully sequenced response regulators and their homologues were determined. The three major, evolutionarily distinct, DNA binding domains found in response regulators were evaluated for their phylogenetic relatedness, and the phylogenetic trees obtained for these domains were compared with those for the receiver modules. Members of one family (family 1) of DNA binding domains are linked to large ATPase domains which usually function cooperatively in the activation of E. coli sigma 54-dependent promoters or their equivalents in other bacteria. Members of a second family (family 2) always function in conjunction with the E. coli sigma 70 or its equivalent in other bacteria. A third family of DNA binding domains (family 3) functions by an uncharacterized mechanism involving more than one sigma factor. These three domain families utilize distinct helix-turn-helix motifs for DNA binding. The phylogenetic tree of the receiver modules revealed three major and several minor clusters of these domains. The three major receiver module clusters (clusters 1, 2, and 3) generally function with the three major families of DNA binding domains (families 1, 2, and 3, respectively) to comprise three classes of response regulators (classes 1, 2, and 3), although several exceptions exist. The minor clusters of receiver modules were usually, but not always, associated with other types of effector domains. Finally, several receiver modules did not fit into a cluster. It was concluded that receiver modules usually diverged from common ancestral protein domains together with the corresponding effector domains, although domain shuffling, due to intragenic splicing and fusion, must have occurred during the evolution of some of these proteins. Multiple sequence alignments of the 49 receiver modules and their various types of effector domains, together with other homologous domains, allowed definition of regions of striking sequence similarity and degrees of conservation of specific residues. Sequence data were correlated with structure/function when such information was available.(ABSTRACT TRUNCATED AT 250 WORDS)

Adenosine Triphosphatases↗

Thrombin-like effect of an important green pit viper toxin, albolabrin: a bioinformatic study.

The green pit viper venom has a major effect on the hematological system. Clinical features of venomous snakebites vary from asymptomatic to fatal bleeding. The venom is found to have a thrombin-like effect in vitro. Here, the author performs a bioinformatic analysis on the green pit viper venom focusing on its thrombin-like effect. Sequence comparison between green pit viper venom, albolabrin and thrombin was performed. In addition, the author performed a search for other human proteins closely relating to the thrombin and created a multiple sequence alignment phylogenetic tree to present the family tree of the thrombin, albolabrin and those proteins recorded in the genomic database. In conclusion, the comparative sequence analysis between green pit viper venom and thrombin gives several identities. The reported relationship on the phylogenetic tree can match with the reported in vivo function of green pit viper. Explanations on the effect of green pit viper toxin on the hemostasis can be derived from this study. Furthermore, future researches based on the reported identities can be expected.

Amino Acid Sequence↗

Identification and characterization of a cathepsin L-like cysteine protease from Taenia solium metacestode.

Taenia solium metacestode, a larval pork tapeworm, is a causative agent of neurocysticercosis, one of the most common parasitic diseases in the human central nervous system. In this study, we identified a cDNA encoding for a cathepsin L-like cysteine protease from the T. solium metacestode (TsCL-1) and characterized the biochemical properties of the recombinant enzyme. The cloned cDNA of 1216 bp encoded 339 amino acids with an approximate molecular weight of 37.6 kDa which containing a typical signal peptide sequence (17 amino acids), a pro-domain (106 amino acids), and a mature domain (216 amino acids). Sequence alignments of TsCL-1 showed low sequence similarity of 27.3-44.6 to cathepsin L-like cysteine proteases from other helminth parasites, but the similarity was increased to 35.9-55.0 when compared to mature domains. The bacterially expressed recombinant protein (rTsCL-1) did not show enzyme activity; however, the rTsCL-1 expressed in Pichia pastoris showed typical biochemical characteristics of cysteine proteases. It degraded human immunoglobulin G (IgG) and bovine serum albumin (BSA), but not collagen. Western blot analysis of the rTsCL-1 showed antigenicity against the sera from patients with cysticercosis, sparganosis or fascioliasis, but weak or no antigenicity against the sera from patients with paragonimiasis or clonorchiasis.

Amino Acid Sequence↗

Glutamine-binding protein from Escherichia coli specifically binds a wheat gliadin peptide. 2. Resonance energy transfer studies suggest a new sensing approach for an easy detection of wheat gliadin.

In this work is presented the first attempt to develop a fluorescence assay for detection of traces of gluten in food by utilizing the recombinant glutamine-binding protein (GlnBP) from E. coli. We found that GlnBP specifically binds the sequence of amino acids present both in gliadin and other prolamines classified as toxic for celiac patients. Affinity chromatography experiments together with mass spectrometry experiments demonstrated that GlnBP can bind the following amino acid sequence XXQPQPQQQQQQQQQQQQL. Sequence alignment experiments pointed out that this sequence is exclusively representative of the gliadin and the other prolamines considered toxic for celiac patients. These findings suggest the development of a competitive resonance energy transfer (RET) assay for an easy and rapid detection of this sequence in raw and cooked food.

Amino Acid Sequence↗

Follow-up study of intrahost HIV type 2 variability reveals discontinuous evolution of C2V3 sequences.

Sequence change dynamics in the highly polymorphic C2V3 region of the envelope glycoprotein gp105 coding sequence of HIV-2 was studied for two HIV-2 seropositive Guinean individuals over a 5-year period. Proviral DNA was amplified by PCR from uncultured peripheral blood mononuclear cells (PBMC), cloned, and sequenced. Amino acid sequence alignment and phylogenetic analysis showed distinct viral populations for the three collected samples (1992, 1995, and 1997) for both individuals. A discontinuous replacement of sequences suggests activation of latent viruses present in reservoirs, a process that has been documented for the HIV-1 infection but unreported so far for HIV-2. The nucleotide sequences presented in this work have been assigned accession numbers AJ400276 to AJ400319.

Amino Acid Sequence↗

QuasiMotiFinder: protein annotation by searching for evolutionarily conserved motif-like patterns.

Sequence signature databases such as PROSITE, which include amino acid segments that are indicative of a protein's function, are useful for protein annotation. Lamentably, the annotation is not always accurate. A signature may be falsely detected in a protein that does not carry out the associated function (false positive prediction, FP) or may be overlooked in a protein that does carry out the function (false negative prediction, FN). A new approach has emerged in which a signature is replaced with a sequence profile, calculated based on multiple sequence alignment (MSA) of homologous proteins that share the same function. This approach, which is superior to the simple pattern search, essentially searches with the sequence of the query protein against an MSA library. We suggest here an alternative approach, implemented in the QuasiMotiFinder web server (http://quasimotifinder.tau.ac.il/), which is based on a search with an MSA of homologous query proteins against the original PROSITE signatures. The explicit use of the average evolutionary conservation of the signature in the query proteins significantly reduces the rate of FP prediction compared with the simple pattern search. QuasiMotiFinder also has a reduced rate of FN prediction compared with simple pattern searches, since the traditional search for precise signatures has been replaced by a permissive search for signature-like patterns that are physicochemically similar to known signatures. Overall, QuasiMotiFinder and the profile search are comparable to each other in terms of performance. They are also complementary to each other in that signatures that are falsely detected in (or overlooked by) one may be correctly detected by the other.

Amino Acid Motifs↗

Predicting the transmembrane secondary structure of ligand-gated ion channels.

Recent mutational analyses of ligand-gated ion channels (LGICs) have demonstrated a plausible site of anesthetic action within their transmembrane domains. Although there is a consensus that the transmembrane domain is formed from four membrane-spanning segments, the secondary structure of these segments is not known. We utilized 10 state-of-the-art bioinformatics techniques to predict the transmembrane topology of the tetrameric regions within six members of the LGIC family that are relevant to anesthetic action. They are the human forms of the GABA alpha 1 receptor, the glycine alpha 1 receptor, the 5HT3 serotonin receptor, the nicotinic AChR alpha 4 and alpha 7 receptors and the Torpedo nAChR alpha 1 receptor. The algorithms utilized were HMMTOP, TMHMM, TMPred, PHDhtm, DAS, TMFinder, SOSUI, TMAP, MEMSAT and TOPPred2. The resulting predictions were superimposed on to a multiple sequence alignment of the six amino acid sequences created using the CLUSTAL W algorithm. There was a clear statistical consensus for the presence of four alpha helices in those regions experimentally thought to span the membrane. The consensus of 10 topology prediction techniques supports the hypothesis that the transmembrane subunits of the LGICs are tetrameric bundles of alpha helices.

Amino Acid Sequence↗

Identification of a gag protein epitope conserved among all four groups of primate immunodeficiency viruses by using monoclonal antibodies.

Five monoclonal antibodies (MAbs) were raised against the gag proteins of simian immunodeficiency virus (SIV) from African green monkey (SIVagmTYO-7). Two MAbs reacted with the matrix protein p17 and the other three with the core protein p24. Studies on the cross-reactivity of the MAbs revealed that the anti-p24 MAbs detected an epitope shared by the viruses belonging to the human immunodeficiency virus type 2 (HIV-2)/SIVmac group and SIVagmTYO-7 and SIVagmTYO-5. The anti-p17 MAbs recognized an epitope present on all these viruses and on SIVagmTYO-1, HIV-1 and SIVmnd. This finding demonstrates for the first time that the matrix protein, p17 or p18, respectively, of all nine HIV and SIV isolates tested in this study expresses at least one conserved immunogenic epitope recognized serologically. By using synthetic peptides, this epitope was identified at the N terminus of p17. Furthermore, this epitope was analysed by multiple sequence alignments of the peptide with homologous sequences of HIV and SIV p17.

Amino Acid Sequence↗

Recognition of another member of the malignant catarrhal fever virus group: an endemic gammaherpesvirus in domestic goats.

A novel gammaherpesvirus in goats that is herein tentatively designated as caprine herpesvirus-2 was identified based on the sequence of a fragment from the herpesvirus DNA polymerase gene. Sequence alignment analysis revealed that the virus sequence isolated from goats was 67% identical to the homologous sequence from alcelaphine herpesvirus-1, 71% identical to ovine herpesvirus-2 and 73% identical to a recently recognized herpesvirus causing malignant catarrhal fever in white-tailed deer. Combined serological and PCR-survey data demonstrated that this virus is endemic in goats and its transmission pattern may be similar to that of ovine herpesvirus-2 in sheep.

Amino Acid Sequence↗

The Adaptive Evolution Database (TAED).

BACKGROUND: Developing an understanding of the molecular basis for the divergence of species lies at the heart of biology. The Adaptive Evolution Database (TAED) serves as a starting point to link events that occur at the same time in the evolutionary history (tree of life) of species, based upon coding sequence evolution analyzed with the Master Catalog. The Master Catalog is a collection of evolutionary models, including multiple sequence alignments, phylogenetic trees, and reconstructed ancestral sequences, for all independently evolving protein sequence modules encoded by genes in GenBank [1]. RESULTS: We have estimated from these models the ratio of nonsynonymous to synonymous nucleotide substitution (Ka/Ks), for each branch in their respective evolutionary trees of every subtree containing only chordata or only embryophyta proteins. Branches with high Ka/Ks values represent candidate episodes in the history of the family where the protein may have undergone positive selection, a phenomenon in molecular evolution where the mutant form of a gene must have conferred more fitness than the ancestral form. Such episodes are frequently associated with change in function. We have found that an unexpectedly large number of families (between 10 and 20% of those families examined) have at least one branch with a notably high Ka/Ks value (putative adaptive evolution). As a resource for biologists wishing to understand the interaction between protein sequences and the Darwinian processes that shape these sequences, we have collected these into The Adaptive Evolution Database (TAED). CONCLUSIONS: Placed in a phylogenetic perspective, candidate genes that are undergoing evolution at the same time in the same lineage can be viewed together. This framework based upon coding sequence evolution can be readily expanded to include other types of evolution. In its present form, TAED provides a resource for bioinformaticists interested in data mining and for experimental evolutionists seeking candidate examples of adaptive evolution for further experimental study.

Animals↗

The adaptive evolution database (TAED).

BACKGROUND: The Master Catalog is a collection of evolutionary families, including multiple sequence alignments, phylogenetic trees and reconstructed ancestral sequences, for all protein-sequence modules encoded by genes in GenBank. It can therefore support large-scale genomic surveys, of which we present here The Adaptive Evolution Database (TAED). In TAED, potential examples of positive adaptation are identified by high values for the normalized ratio of nonsynonymous to synonymous nucleotide substitution rates (KA/KS values) on branches of an evolutionary tree between nodes representing reconstructed ancestral sequences. RESULTS: Evolutionary trees and reconstructed ancestral sequences were extracted from the Master Catalog for every subtree containing proteins from the Chordata only or the Embryophyta only. Branches with high KA/KS values were identified. These represent candidate episodes in the history of the protein family when the protein may have undergone positive selection, where the mutant form conferred more fitness than the ancestral form. Such episodes are frequently associated with change in function. An unexpectedly large number of families (between 10% and 20% of those families examined) were found to have at least one branch with high KA/KS values above arbitrarily chosen cut-offs (1 and 0.6). Most of these survived a robustness test and were collected into TAED. CONCLUSIONS: TAED is a raw resource for bioinformaticists interested in data mining and for experimental evolutionists seeking candidate examples of adaptive evolution for further experimental study. It can be expanded to include other evolutionary information (for example changes in gene regulation or splicing) placed in a phylogenetic perspective.

Adaptation, Physiological↗

A strategy for the identification of T-cell epitopes on Leishmania cysteine proteinases.

In this study computational analysis was used to compile sequence alignments, construct a dendrogram and calculate physical data in order to predict potential T-cell epitopes of the Leishmania cysteine proteinase. Using multiple alignment of human and Leishmania proteinase sequences deposited on data bank sequences, it was possible to predict that the extreme C-terminus of cysteine proteinase (Cyspep, 355-444) contained three peptides (pI 361-370, pII 415-422 and pIII 431-444) with charge score, hydrophobicity and isoelectric points compatible for human leucocyte-associated antigen (HLA) class II binding. The prediction was confirmed in vitro through the ability of synthetic peptides corresponding to the predicted regions to stimulate peripheral blood mononuclear cells of patients with leishmaniasis.

Algorithms↗