PubMed HealthSearch

SEARCH · PubMed Health

Results for “Multiple 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 235 records · Page 13Linked to original sources

Evolutionary relationships among putative RNA-dependent RNA polymerases encoded by a mitochondrial virus-like RNA in the Dutch elm disease fungus, Ophiostoma novo-ulmi, by other viruses and virus-like RNAs and by the Arabidopsis mitochondrial genome.

The nucleotide sequence (2617 nucleotides) of virus-like double-stranded (ds) RNA 3a in a diseased isolate, Log1/3-8d2 (Ld), of the ascomycete fungus Ophiostoma novo-ulmi has been determined. One strand of the dsRNA contains an open reading frame (ORF) with the potential to encode a protein of 718 amino acids, and the complementary strand contains two smaller ORFs with the potential to encode proteins of 178 and 182 amino acids, respectively. The large ORF contains 12 UGA codons which code for tryptophan in ascomycete mitochondria and has a codon bias typical of mitochondrial genes, consistent with the localization of Ld dsRNAs within the mitochondria. The amino acid sequence contains motifs characteristic of RNA-dependent RNA polymerases (RdRps). This putative RdRp was shown to be related to putative RdRps of mitochondrial dsRNAs of another ascomycete and a basidiomycete fungus and also to a putative RdRp encoded by the mitochondrial genome of Arabidopsis thaliana. In multiple sequence alignments, the fungal mitochondrial dsRNA-encoded RdRp-like proteins formed a cluster, ancestrally related to the RdRps of the yeast 20S and 23S RNA replicons and of the positive-stranded RNA bacteriophages of the Leviviridae family, but distinct from RdRps of other families and genera of fungal RNA viruses and related plant and animal RNA viruses. Northern blot analysis with RNA 3a strand-specific probes indicated that nucleic acid extracts of Ld contain more single-stranded (positive-stranded) RNA than dsRNA, consistent with an evolutionary relationship between RNA 3a and positive-stranded RNA phages.

Amino Acid Sequence

Protein engineering in the alpha-amylase family: catalytic mechanism, substrate specificity, and stability.

Most starch hydrolases and related enzymes belong to the alpha-amylase family which contains a characteristic catalytic (beta/alpha)8-barrel domain. Currently known primary structures that have sequence similarities represent 18 different specificities, including starch branching enzyme. Crystal structures have been reported in three of these enzyme classes: the alpha-amylases, the cyclodextrin glucanotransferases, and the oligo-1,6-glucosidases. Throughout the alpha-amylase family, only eight amino acid residues are invariant, seven at the active site and a glycine in a short turn. However, comparison of three-dimensional models with a multiple sequence alignment suggests that the diversity in specificity arises by variation in substrate binding at the beta-->alpha loops. Designed mutations thus have enhanced transferase activity and altered the oligosaccharide product patterns of alpha-amylases, changed the distribution of alpha-, beta- and gamma-cyclodextrin production by cyclodextrin glucanotransferases, and shifted the relative alpha-1,4:alpha-1,6 dual-bond specificity of neopullulanase. Barley alpha-amylase isozyme hybrids and Bacillus alpha-amylases demonstrate the impact of a small domain B protruding from the (beta/alpha)8-scaffold on the function and stability. Prospects for rational engineering in this family include important members of plant origin, such as alpha-amylase, starch branching and debranching enzymes, and amylomaltase.

Amino Acid Sequence

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

The evolutionary divergence of neurotransmitter receptors and second-messenger pathways.

Members of the superfamily of G-protein-coupled neurotransmitter receptors have a conserved secondary structure, a moderate and reasonably steady rate of sequence change, and usually lack introns within the coding sequence. These properties are advantageous for evolutionary studies. The duplication and divergence of the genes in this gene family led to the formation of distinct neurotransmitter pathways and may have facilitated the evolution of complex nervous systems. I have analyzed this evolutionary divergence by quantitative multiple sequence alignment, bootstrap resampling, and statistical analysis of 49 adrenergic, muscarinic cholinergic, dopamine, and octopamine receptor sequences from 12 animal species. The results indicate that the first event to occur within this gene family was the divergence of the catecholamine receptors from the muscarinic acetylcholine receptors, which occurred prior to the divergence of the arthropod and vertebrate lineages. Subsequently, the ability to activate specific second-messenger pathways diverged independently in both the muscarinic and the catecholamine receptors. This appears to have occurred after the divergence of the arthropod and vertebrate lineages but before the divergence of the avian and mammalian lineages. However, the second-messenger pathways activated by adrenergic and dopamine receptors did not diverge independently. Rather, the ability of the catecholamine receptors to bind to specific ligands, such as epinephrine, norepinephrine, dopamine, or octopamine, was repeatedly modified in evolutionary history, and in some cases was modified after the divergence of the second-messenger pathways.

Adenylyl Cyclases

Comparison of lantibiotic gene clusters and encoded proteins.

Lantibiotics form a group of modified peptides with unique structures, containing post-translationally modified amino acids such as dehydrated and lanthionine residues. In the gram-positive bacteria that secrete these lantibiotics, the gene clusters flanking the structural genes for various linear (type A) lantibiotics have recently been characterized. The best studied representatives are those of nisin (nis), subtilin (spa), epidermin (epi), Pep5 (pep), cytolysin (cyl), lactocin S (las) and lacticin 481 (lct). Comparison of the lantibiotic gene clusters shows that they contain conserved genes that probably encode similar functions. The nis, spa, epi and pep clusters contain lanB and lanC genes that are presumed to code for two types of enzymes that have been implicated in the modification reactions characteristic of all lantibiotics, i.e. dehydration and thio-ether ring formation. The cyl, las and lct gene clusters have no homologue of the lanB gene, but they do contain a much larger lanM gene that is the lanC gene homologue. Most lantibiotic gene clusters contain a lanP gene encoding a serine protease that is presumably involved in the proteolytic processing of the prelantibiotics. All clusters contain a lanT gene encoding an ABC transporter likely to be involved in the export of (precursors of) the lantibiotics. The lanE, lanF and lanG genes in the nis, spa and epi clusters encode another transport system that is possibly involved in self-protection. In the nisin and subtilin gene clusters two tandem genes, lanR and lanK, have been located that code for a two-component regulatory system. Finally, non-homologous genes are found in some lantibiotic gene clusters. The nisI and spaI genes encode lipoproteins that are involved in immunity, the pepI gene encodes a membrane-located immunity protein, and epiD encodes an enzyme involved in a post-translational modification found only in the C-terminus of epidermin. Several genes of unknown function are also found in the las gene cluster. A database has been assembled for all putative gene products of type A lantibiotic gene clusters. Database searches, multiple sequence alignment and secondary structure prediction have been used to identify conserved sequence segments in the LanB, LanC, LanE, LanF, LanG, LanK, LanM, LanP, LanR and LanT gene products that may be essential for structure and function. This database allows for a rapid screening of newly determined sequences in lantibiotic gene clusters.

ATP-Binding Cassette Transporters

Identification of novel homologues of three low molecular weight subunits of the mitochondrial bc1 complex.

Large-scale random cDNA sequencing projects have been started for several organisms and are a valuable tool for the analysis of quantitative and qualitative aspects of gene expression. However, the reliability of the obtained data is limited as most of the clones are only partially analysed on one strand. As a consequence the sequence entries derived from random cDNA sequencing projects usually comprise incomplete open reading frames. They nevertheless define complete and reliable coding sequences, if two prerequisites are fulfilled: (i) the clones encode very small proteins, and (ii) the clones have a high frequency in the cDNA-banks. The present study describes the use of cDNA databases for the identification of homologues of three low-molecular-weight subunits of the mitochondrial bc1 complex, termed the QCR6, QCR9 and QCR10 proteins. These polypeptides are only characterized for a small number of organisms, have a scarcely defined function and exhibit a low degree of structural conservation if compared between different species. Several clones were identified for each polypeptide by searches with TBLASTN using the known sequences as probes. Most of the database entries contain complete open reading frames and sequencing queries could be excluded due to the abundancy of the clones. Multiple sequence alignments are presented for all three polypeptides and consensus sequences are given which may provide a basis for the investigation of the proteins by site-directed mutagenesis.

Animals

Taxonomic relationships between distinct potato virus Y isolates based on detailed comparisons of the viral coat proteins and 3'-nontranslated regions.

Detailed comparisons were made of the sequences of the coat protein (CP) cistrons and 3'-nontranslated regions (3'-NTR) of 21 (geographically) distinct isolates of potato virus Y (PVY) and a virus isolate initially described as pepper mottle virus (PepMoV). Multiple sequence alignments and phylogenetic relationships based on these alignments resulted into a subgrouping of virus isolates which largely corresponded with the historical strain differentiation based on biological criteria as host range, symptomatology and serology. Virus isolates belonging to the same subgroup shared a number of characteristic CP amino acid and 3'-NTR nucleotide residues indicating that, by using sequences from the 3'-terminal region of the potyvirus genome, a distinction could be made between different isolates of one virus species as well as between different virus species. RNA secondary structure analysis of the 3'-NTR of twelve PVY isolates revealed four major stem-loop structures of which, surprisingly, the loop sequences gave a similar clustering of isolates as resulting from the overall comparisons of CP and 3'-NTR sequences. This implies a biological significance of these structural elements.

Amino Acid Sequence

Predicted structure of the adenovirus DNA binding protein.

The DNA sequence of a portion of the MAV1 SmaI-D fragment coding for the C-terminal 147 amino acids of the adenoviral DNA-binding protein (DBP) has been determined. A multiple sequence alignment was constructed of the MAV1 fragment and the DBPs of Ad.2, 4, 5, 7, 12, 40, and 41 to examine the degree of conservation of features that have been mapped on the Ad.2 DBP and to identify further conserved features. The less conserved N-terminal segment of the protein contains two nuclear localization signals and two acidic regions, the host range region, and all of the 11 phosphorylation sites. The highly conserved C-terminal segment contains a potential leucine zipper and zinc finger motifs. These sequence features were mapped onto a predicted secondary structure of the Ad.2 DBP.

Adenoviridae

Sequence of the avian adenovirus FAV 1 (CELO) DNA encoding the hexon-associated protein pVI and hexon.

The genomic region of the avian adenovirus FAV1 (CELO) encoding the precursor to virion structural protein VI (p VI) and the major capsid protein hexon has been sequenced. The 223-unit sequence of the CELO pVI protein has two potential Ad endoproteinase cleavage sites and a conserved C-terminal sequence including the Cys residue supposedly involved in endoproteinase activation. The CELO hexon gene sequence predicts a 942-residue protein (106.7 kDa). Multiple sequence alignment with other six known hexon protein sequences (human, bovine, murine, and avian) reveals high overall homology. The identity is highest in the regions corresponding to the pedestals which from the base of the hexon, and lowest in the regions corresponding to the loops which are exposed on the outer surface of the virion.

Adenoviridae

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

The evolution of rhodopsins and neurotransmitter receptors.

Rhodopsins share a limited number of amino acid identities with a variety of other integral membrane proteins. Most of these proteins have seven putative transmembrane segments and are likely to play a role in transmembrane signaling. We have undertaken a systematic series of comparisons of primary and secondary structure in order to clarify the functional and evolutionary significance of these sequence similarities. On the basis of consistently high similarity scores, we find that the most internally consistent definition of the rhodopsin gene family would include vertebrate rhodopsins, alpha- and beta-adrenergic receptors, M1 and M2 muscarinic acetylcholine receptors, substance K receptors, and insect rhodopsins, while excluding bacteriorhodopsin, the mass human oncogene, vertebrate and insect nicotinic acetylcholine receptors, and the yeast STE2 and STE3 peptide receptors. The rhodopsin gene family is highly diverged at the primary sequence level but has maintained a conserved secondary structure, including a previously unidentified hierarchy of transmembrane segment hydrophobicity. We have developed new computer algorithms for progressive multiple sequence alignment and the analysis of local conservation of protein domains, and we have used these algorithms to examine the phylogeny of the rhodopsin gene family and the changing domains of sequence conservation. The results show striking differences and similarities in the conserved domains in each of the three main branches of the rhodopsin gene family, and indicate that color vision arose independently in the lines of descent leading to modern humans and fruit flies.

Algorithms

Common origin of arthropod tyrosinase, arthropod hemocyanin, insect hexamerin, and dipteran arylphorin receptor.

Dipteran arylphorin receptors, insect hexamerins, cheliceratan and crustacean hemocyanins, and crustacean and insect tyrosinases display significant sequence similarities. We have undertaken a systematic comparison of primary and secondary structures of these proteins. On the basis of multiple sequence alignments the phylogeny of these proteins was investigated. Hexamerin subunits, hemocyanin subunits, and tyrosinases share extensive similarities throughout the entire amino acid sequence. Our studies suggest the origin of arthropod hemocyanins from ancient tyrosinase-like proteins. Insect hexamerins likely evolved from hemocyanins of ancient crustaceans, supporting the proposed sister-group position of these subphyla. Arylphorin receptors, responsible for incorporation of hexamerins into the larval fat body of diptera, are related to hexamerins, hemocyanins, and tyrosinase. The receptor sequences display extensive similarities to the first and third domains of hemocyanins and hexamerins. In the middle region only limited amino acid conservation was observed. Elements important for hexamer formation are deleted in the receptors. Phylogenetic analysis indicated that dipteran arylphorin receptors diverged from ancient hexamerins, probably early in insect evolution.

Amino Acid Sequence

Structural and functional relationships of human DNA polymerases.

A continuing theme of our laboratory has been the understanding of human DNA polymerases at the structural level. We have purified DNA polymerases delta, epsilon and alpha from human placenta. Monoclonal antibodies to these polymerases were isolated and used as tools to study their immunochemical relationships. These studies have shown that while DNA polymerases delta, epsilon and alpha are discrete proteins, they must share common structural features by virtue of the ability of several of our monoclonal antibodies to exhibit cross-reactivity. A second approach we have taken is the molecular cloning of human DNA polymerase delta and epsilon. We have cloned the DNA polymerase delta cDNA, and this has allowed us to compare its primary structure to those of human polymerase alpha and other members of this polymerase family. Multiple sequence alignments have revealed that human DNA polymerase delta is also closely related to the herpes virus family of DNA polymerases. In situ hybridization has shown that the human DNA polymerase delta gene is localized to chromosome 19 q13.3-q13.4. In order to further determine the functional regions of the DNA polymerase delta structure we are currently expressing human pol delta in E. coli and baculovirus systems. Other work in our laboratory is directed toward examining the expression of DNA polymerase delta during the cell cycle.

Amino Acid Sequence

The evolution of hexamerins and the phylogeny of insects.

The evolutionary relationships among arthropod hemocyanins and insect hexamerins were investigated. A multiple sequence alignment of 12 hemocyanin and 31 hexamerin subunits was constructed and used for studying sequence conservation and protein phylogeny. Although hexamerins and hemocyanins belong to a highly divergent protein superfamily and only 18 amino acid positions are identical in all the sequences, the core structures of the three protein domains are well conserved. Under the assumption of maximum parsimony, a phylogenetic tree was obtained that matches perfectly the assumed phylogeny of the insect orders. An interesting common clade of the hymenopteran and coleopteran hexamerins was observed. In most insect orders, several paralogous hexamerin subclasses were identified that diversified after the splitting of the major insect orders. The dipteran arylphorin/LSP-1-like hexamerins were subject to closer examination, demonstrating hexamerin gene amplification and gene loss in the brachyceran Diptera. The hexamerin receptors, which belong to the hexamerin/hemocyanin superfamily, diverged early in insect evolution, before the radiation of the winged insects. After the elimination of some rapidly or slowly evolving sequences, a linearized phylogenetic tree of the hexamerins was constructed under the assumption of a molecular clock. The inferred time scale of hexamerin evolution, which dates back to the Carboniferous, agrees with the available paleontological data and reveals some previously unknown divergence times among and within the insect orders.

Amino Acid Sequence

Cloning of the homogentisate 1,2-dioxygenase gene, the key enzyme of alkaptonuria in mouse.

We determined 48 amino acid residues from five peptides from the homogeneous monomer of homogentisate 1,2-dioxygenase (HGO; E.C. 1.13. 11.15) of mouse liver. After digestion with trypsin, peptides were separated by reversed phase chromatography and amino acid sequenced. The deduced codon sequence of three peptides was used to derive degenerated oligomeres. By combining these oligos, we were able to amplify fragments from 100 to 300 bases (b) from mouse liver cDNA by polymerase chain reaction after reverse transcription (RT-PCR). A fragment of 200 b was cloned and used as a probe to screen a mouse liver cDNA library. One clone from this library contained the complete cDNA-insert for HGO as determined by sequencing. The cDNA encodes for a protein of 50 kDa, as predicted. The cDNA of mouse HGO has an overall identity of 41% to the corresponding gene hmgA from Aspergillus. Sequence similarities to human expressed sequence tags (EST) clones ranged from 70% to 20%. The positions of 122 conserved amino acids could be determined by multiple sequence alignment. We identified one first intron of 928 b in the mouse gene. The gene for HGO seems to be expressed in various tissues, as shown by RT-PCR on different cDNAs. FISH experiments with the whole murine cDNA as probe clearly revealed signals at the human chromosomal band 3q13. 3-q21. This corresponds well to the previous assignment of the locus for the human alkaptonuria gene (AKU) to the same chromosomal region by multipoint linkage analysis. We therefore conclude that the HGO cDNA encodes the gene responsible for alkaptonuria.

Alkaptonuria

High sequence similarity within ras exons 1 and 2 in different mammalian species and phylogenetic divergence of the ras gene family.

We have determined the canine and feline N-, K-, and H-ras gene sequences from position +23 to +270 covering exons I and II which contain the mutational hot spot codons 12, 13, and 61. The results were used to assess the degree of similarity between ras gene DNA regions containing the critical domains affected in neoplastic disorders in different mammalian species. The comparative analyses performed included human, canine, feline, murine, rattine, and, whenever possible, bovine, leporine (rabbit), porcelline (guinea pig), and mesocricetine (hamster) ras gene sequences within the region of interest. Comparison of feline and canine nucleotide sequences with the corresponding regions in human DNA revealed a sequence similarity greater than 85% to the human sequence. Contemporaneous analysis of previously published ras DNA sequences from other mammalian species showed a similar degree of homology to human DNA. Most nucleotide differences observed represented synonymous changes without effect on the amino acid sequence of the respective proteins. For assessment of the phylogenetic evolution of ras gene family, a maximum parsimony dendrogram based on multiple sequence alignment of the common region of exons I and II in the N-, K-, and H-ras genes was constructed. Interestingly, a higher substitution rate among the H-ras genes became apparent, indicating accelerated sequence evolution within this particular clade. The most parsimonious tree clearly shows that the duplications giving rise to the three ras genes must have occurred before the mammalian radiation.

Animals

A helix-turn-helix DNA-binding motif predicted for transposases of DNA transposons.

A helix-turn-helix (HTH) DNA-binding motif is identified in transposase sequences in Tc1, mariner and pogo DNA transposum. The findings are supported by results of various sequence analysis methods. Tc1 transposases are also predicted to contain another DNA-binding region. These findings are in accord with experimental evidence obtained from Tc1A, Tc3A and pogo transposases. The pogo family transposases, but not the pogo-type transcription factors, contain the HTH motif, suggesting that HTH structures are essential for Tc1/mariner/pogo transposition. Analysis of multiple sequence alignments enabled the identification of the HTH motif in distantly related protein sequences.

Amino Acid Sequence

Molecular characterization of morphologically typical human calicivirus Sapporo.

Human calicivirus Sapporo (SV) has typical calicivirus morphology and causes acute gastroenteritis in children. The nucleotide sequence of 3.2 kb of the 3' end of SV was determined from a cloned cDNA. The 3' end of the SV genome is predicted to encode the RNA-dependent RNA polymerase region, the capsid protein and two small open reading frames. The nonstructural and capsid protein coding sequences in the SV genome are fused in a single open reading frame. The organization of these proteins in the SV sequence is similar to that of rabbit hemorrhagic disease virus and the recently described Manchester virus, and distinct from the genome organization of the prototype human calicivirus, Norwalk virus, that lacks typical calicivirus morphology and has been described as a small round structured virus (SRSV). Sequence analysis of the predicted capsid region showed that the SV capsid is longer by approximately 30 amino acids than the capsid of any of the SRSVs, and multiple sequence alignments showed that these additional amino acids are located in the variable region of the capsid protein. Expression of the capsid protein of SV in insect cells resulted in the self-assembly of virus-like particles that have a morphology similar to that of the native virus. This result shows that calicivirus morphology is determined by the primary sequence of the capsid protein.

Amino Acid Sequence