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

PubMed1991

A molecular phylogeny of dinoflagellate protists (pyrrhophyta) inferred from the sequence of 24S rRNA divergent domains D1 and D8.

The sequence of two divergent domains (D1 and D8) from dinoflagellate 24S large subunit rRNA was determined by primer extension using total RNA as template. Nucleotide sequence alignments over 401 bases have been analyzed in order to investigate phylogenetic relationships within this highly divergent and taxonomically controversial group of protists of the division Pyrrhophyta. Data are provided confirming that dinoflagellates represent a monophyletic group. For 11 out of the 13 investigated laboratory grown species, an additional domain (D2) could not be completely sequenced by reverse transcription because of a hidden break located near its 3'-terminus. Two sets of sequence alignments were used to infer dinoflagellate phylogeny. The first [199 nucleotides (nt)] included conservative sequences flanking the D1 and D8 divergent domains. It was used to reconstruct a broad evolutionary tree for the dinoflagellates, which was rooted using Tetrahymena thermophila as the outgroup. To confirm the tree topology, and mainly the branchings leading to closely related species, a second alignment (401 nt) was considered, which included the D1 and D8 variable sequences in addition to the more conserved flanking regions. Species that showed sequence similarities with other species lower than 60% on average (Knuc values higher than 0.550) were removed from this analysis. A coherent and convincing evolutionary pattern was obtained for the dinoflagellates, also confirmed by the position of the hidden break within the D2 domain, which appears to be group specific. The reconstructed phylogeny indicates that the early emergence of Oxyrrhis marina preceded that of most Peridiniales, a large order of thecate species, whereas the unarmored Gymnodiniales appeared more recently, along with members of the Prorocentrales characterized by two thecal plates. In addition, the emergence of heterotrophic species preceded that of photosynthetic species. These results provide new perspectives on proposed evolutionary trees for the dinoflagellates based on morphology, biology, and fossil records.

Animals↗

PubMed2003

Conserved structural elements in glutathione transferase homologues encoded in the genome of Escherichia coli.

Multiple sequence alignments of the eight glutathione (GSH) transferase homologues encoded in the genome of Escherichia coli were used to define a consensus sequence for the proteins. The consensus sequence was analyzed in the context of the three-dimensional structure of the gst gene product (EGST) obtained from two different crystal forms of the enzyme. The enzyme consists of two domains. The N-terminal region (domain I) has a thioredoxin-like alpha/beta-fold, while the C-terminal domain (domain II) is all alpha-helical. The majority of the consensus residues (12/17) reside in the N-terminal domain. Fifteen of the 17 residues are involved in hydrophobic core interactions, turns, or electrostatic interactions between the two domains. The results suggest that all of the homologues retain a well-defined group of structural elements both in and between the N-terminal alpha/beta domain and the C-terminal domain. The conservation of two key residues for the recognition motif for the gamma-glutamyl-portion of GSH indicates that the homologues may interact with GSH or GSH analogues such as glutathionylspermidine or alpha-amino acids. The genome context of two of the homologues forms the basis for a hypothesis that the b2989 and yibF gene products are involved in glutathionylspermidine and selenium biochemistry, respectively.

Amino Acid Sequence↗

PubMed2002

The evolution of the vertebrate beta-globin gene promoter.

Complexity analysis is capable of highlighting those gross evolutionary changes in gene promoter regions (loosely termed "promoter shuffling") that are undetectable by conventional DNA sequence alignment. Complexity analysis was therefore used here to identify the modular components (blocks) of the orthologous beta-globin gene promoter sequences of 22 vertebrate species, from zebrafish to humans. Considerable variation between the beta-globin gene promoters was apparent in terms of block presence/absence, copy number, and relative location. Some sequence blocks appear to be ubiquitous, whereas others are restricted to a specific taxon. Block similarities were also evident between the promoters of the paralogous human beta-like globin genes. It may be inferred that a wide variety of different mutational mechanisms have operated upon the beta-globin gene promoter over evolutionary time. Because these include gross changes such as deletion, duplication, amplification, elongation, contraction, and fusion, as well as the steady accumulation of single base-pair substitutions, it is clear that some redefinition of the term "promoter shuffling" is required. This notwithstanding, and as previously described for the vertebrate growth hormone gene promoter, the modular structure of the beta-globin promoter region and those of its paralogous counterparts have continually been rearranged into new combinations through the alteration, or shuffling, of preexisting blocks. Some of these changes may have had no influence on promoter function, but others could have altered either the level of gene expression or the responsiveness of the promoter to external stimuli. The comparative study of vertebrate beta-globin gene promoter regions described here confirms the generality of the phenomenon of sequence block shuffling and thus supports the view that it could have played an important role in the evolution of differential gene expression.

Amino Acid Sequence↗

PubMed2004

ConSSeq: a web-based application for analysis of amino acid conservation based on HSSP database and within context of structure.

SUMMARY: A web-based application to analyze protein amino acids conservation-Consensus Sequence (ConSSeq) is presented. ConSSeq graphically represents information about amino acid conservation based on sequence alignments reported in homology-derived structures of proteins. Beyond the relative entropy for each position in the alignment, ConSSeq also presents the consensus sequence and information about the amino acids, which are predominant at each position of the alignment. ConSSeq is part of the STING Millennium Suite and is implemented as a Java Applet. AVAILABILITY: http://sms.cbi.cnptia.embrapa.br/SMS/STINGm/consseq/, http://trantor.bioc.columbia.edu/SMS/STINGm/consseq/, http://mirrors.rcsb.org//SMS/STINGm/consseq/, http://www.es.embnet.org/SMS/STINGm/consseq/ and http://www.ar.embnet.org/SMS/STINGm/consseq/

Amino Acids↗

PubMed2005

transAlign: using amino acids to facilitate the multiple alignment of protein-coding DNA sequences.

BACKGROUND: Alignments of homologous DNA sequences are crucial for comparative genomics and phylogenetic analysis. However, multiple alignment represents a computationally difficult problem. For protein-coding DNA sequences, it is more advantageous in terms of both speed and accuracy to align the amino-acid sequences specified by the DNA sequences rather than the DNA sequences themselves. Many implementations making use of this concept of "translated alignments" are incomplete in the sense that they require the user to manually translate the DNA sequences and to perform the amino-acid alignment. As such, they are not well suited to large-scale automated alignments of large and/or numerous DNA data sets. RESULTS: transAlign is an open-source Perl script that aligns protein-coding DNA sequences via their amino-acid translations to take advantage of the superior multiple-alignment capabilities and speed of an amino-acid alignment. It operates by translating each DNA sequence into its corresponding amino-acid sequence, passing the entire matrix to ClustalW for alignment, and then back-translating the resulting amino-acid alignment to derive the aligned DNA sequences. In the translation step, transAlign determines the optimal orientation and reading frame for each DNA sequence according to the desired genetic code. It also checks for apparent frame shifts in the DNA sequences and can handle frame-shifted sequences in one of three ways (delete, align as amino acids regardless, or profile align as DNA). As a set of comparative benchmarks derived from six protein-coding genes for mammals shows, the strategy implemented in transAlign always improves the speed and usually the apparent accuracy of the alignment of protein-coding DNA sequences. CONCLUSION: transAlign represents one of few full and cross-platform implementations of the concept of translated alignments. Both the advantages accruing from performing a translated alignment and the suite of user-definable options available in the program mean that transAlign is ideally suited for large-scale automated alignments of very large and/or very numerous protein-coding DNA data sets. However, the good performance offered by the program also translates to the alignment of any set of protein-coding sequences. transAlign, including the source code, is freely available at http://www.tierzucht.tum.de/Bininda-Emonds/ (under "Programs").

Algorithms↗

PubMed2006

[Molecular cloning of an amidase gene from Nocardia sp. and its expressionin Escherichia coli].

The amidase of Nocardia sp. is one of important industrial enzymes. Based on DNA and protein sequence alignment from different strains, a new gene of amidase was successfully cloned from Nocardia YS-2002, which is widely used for industrial production of acrylamide in China. DNA sequence analyses showed that the 1466bp cloned-fragment contains promoter, open reading frame and terminating-palindrome. Protein sequence alignment and phylogenetic tree analyses showed that the amidase coming from Nocardia sp. YS-2002 is a kind of specialamidase, without the typical conserved sequence of the amidases. Enzymatic characteristics predictions indicated that the molecular weight and pI of the new amidase is approximately 38.05 kD and 4.88, respectively, and it would be stable when heterogeneously expressed in E. coli. By inserting the ORF of the amidase into plasmid pET-28a(+), a recombinant strain, pEAB, was selected using E. coli BL21(DE3) as the host. SDS-PAGE analyses of both the whole cells and ultrasonic-treated cells confirmed the feasibility of the heterogeneous expression of amidase in the recombinant E. coli. But the activity of amidase in E. coli BL21(DE3) not more than 0.5 u/mg, because most of the enzymes expressed were formed as inclusion bodies.

Amidohydrolases↗

PubMed2006

FeatureMap3D--a tool to map protein features and sequence conservation onto homologous structures in the PDB.

FeatureMap3D is a web-based tool that maps protein features onto 3D structures. The user provides sequences annotated with any feature of interest, such as post-translational modifications, protease cleavage sites or exonic structure and FeatureMap3D will then search the Protein Data Bank (PDB) for structures of homologous proteins. The results are displayed both as an annotated sequence alignment, where the user-provided annotations as well as the sequence conservation between the query and the target sequence are displayed, and also as a publication-quality image of the 3D protein structure with the selected features and sequence conservation enhanced. The results are also returned in a readily parsable text format as well as a PyMol (http://pymol.sourceforge.net/) script file, which allows the user to easily modify the protein structure image to suit a specific purpose. FeatureMap3D can also be used without sequence annotation, to evaluate the quality of the alignment of the input sequences to the most homologous structures in the PDB, through the sequence conservation colored 3D structure visualization tool. FeatureMap3D is available at: http://www.cbs.dtu.dk/services/FeatureMap3D/.

Amino Acid Sequence↗

PubMed2001

Comparing protein structures: a Gaussian-based approach to the three-dimensional structural similarity of proteins.

This study describes a new method for comparing three-dimensional protein structures based on an optimal alignment of their steric fields. The method is based upon the use of spherical Gaussian functions located on individual atoms. This representation generates a flexible description of the underlying fold geometry of proteins that can be adjusted by changing the 'width' of the Gaussians. Reducing the width sharpens the representation and leads to a more 'atomlike' description; increasing the width creates a fuzzier representation that preserves the general shape features of the chain fold but with a consequent loss in atomic resolution. The width used in this study is based upon the features of individual atoms and provides a representation that is quite robust with respect to the variety of geometric features typically encountered in the alignment process. In addition, a post-alignment analysis is performed that generates sequence alignments from the corresponding structure alignments. An example, based on five mammalian and fungal matrix metalloproteinase crystal structures (human fibroblast collagenase, neutrophil collagenase, stromelysin, astacin, and adamalysin), illustrates a number of features of the Gaussian-based approach.

Amino Acid Sequence↗

PubMed1994

Molecular phylogenetics of cupped oysters based on partial 28S rRNA gene sequences.

Partial sequences of 28S-like rDNA were amplified using PCR and sequenced for eight species of oyster and one species of mussel. Phylogenetic relationships among seven species of Crassostreinid oyster were inferred from aligned sequences by parsimony and maximum-likelihood methods. Of the 315 sites that varied, 90 were phylogenetically informative in parsimony analysis. Inference by maximum parsimony (MP) is consistent with maximum-likelihood (ML) analysis for the major lineages, yielding a tree with the topology (Mytilus edulis (Ostrea edulis ((Crassostrea rivularis (C. belcheri, C. gigas))(C. virginica, C. rhizophorae, Saccostrea cuccullata, S. commercialis)))). MP and ML analyses resolved the systematic relationships of the Saccostrea and Atlantic Crassostrea differently such that a polytomy linking these four taxa is preferred with the data available. Molecular data support a later divergence of the tropical Pacific Saccostrea from a common ancestor of the Atlantic Crassostrea species. Molecular data from domains D1, D2, and partial D3 of the 28S rDNA supply sufficient phylogenetic information to determine systematic relationships among the extant oyster taxa, from the major species groups to the family level, thus providing valuable characters that are able to supplement the paucity of morphological characters so far recognized.

Animals↗

PubMed1996

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↗

PubMed2005

A secondary structural model of the 28S rRNA expansion segments D2 and D3 for Chalcidoid wasps (Hymenoptera: Chalcidoidea).

We analyze the secondary structure of two expansion segments (D2, D3) of the 28S ribosomal (rRNA)-encoding gene region from 527 chalcidoid wasp taxa (Hymenoptera: Chalcidoidea) representing 18 of the 19 extant families. The sequences are compared in a multiple sequence alignment, with secondary structure inferred primarily from the evidence of compensatory base changes in conserved helices of the rRNA molecules. This covariation analysis yielded 36 helices that are composed of base pairs exhibiting positional covariation. Several additional regions are also involved in hydrogen bonding, and they form highly variable base-pairing patterns across the alignment. These are identified as regions of expansion and contraction or regions of slipped-strand compensation. Additionally, 31 single-stranded locales are characterized as regions of ambiguous alignment based on the difficulty in assigning positional homology in the presence of multiple adjacent indels. Based on comparative analysis of these sequences, the largest genetic study on any hymenopteran group to date, we report an annotated secondary structural model for the D2, D3 expansion segments that will prove useful in assigning positional nucleotide homology for phylogeny reconstruction in these and closely related apocritan taxa.

Animals↗

PubMed1998

Incremental threading optimization (TITO) to help alignment and modelling of remote homologues.

MOTIVATION: Protein sequence comparison has become a major tool for biologists. Various methods have been designed in order to reveal even remote homologies among proteins. Sequence analyses are powerful, but the lower the identity score, the more skill and time are required to perform them accurately. While necessary to confirm the sequence comparison through the potential of mean force, sequence alignment and molecular modelling are still tedious and time-consuming tasks. RESULTS: In order to help sequence threading onto known three-dimensional (3D) structures, we have developed a new program, named TITO (Tool for Incremental Threading Optimization), which uses a multiple alignment to validate and help the refinement of a sequence/structure comparison. It measures the compatibility of a family of related sequences with a known 3D structure. The coordinates of the common core are extracted, allowing further molecular modelling and 3D survey. AVAILABILITY: TITO is available by anonymous ftp at: ftp://lmcp.jussieu. fr/pub/tito CONTACT: labesse@lmcp.jussieu.fr SUPPLEMENTARY INFORMATION: http://www.lmcp.jussieu.fr/ approximately labesse/TITO/TITO.html

Algorithms↗

PubMed2004

A new criterion and method for amino acid classification.

It is accepted that many evolutionary changes of amino acid sequence in proteins are conservative: the replacement of one amino acid by another residue has a far greater chance of being accepted if the two residues have similar properties. It is difficult, however, to identify relevant physicochemical properties that capture this similarity. In this paper we introduce a criterion that determines similarity from an evolutionary point of view. Our criterion is based on the description of protein evolution by a Markov process and the corresponding matrix of instantaneous replacement rates. It is inspired by the conductance, a quantity that reflects the strength of mixing in a Markov process. Furthermore we introduce a method to divide the 20 amino acid residues into subsets that achieve good scores with our criterion. The criterion has the time-invariance property that different time distances of the same amino acid replacement rate matrix lead to the same grouping; but different rate matrices lead to different groupings. Therefore it can be used as an automated method to compare matrices derived from consideration of different types of proteins, or from parts of proteins sharing different structural or functional features. We present the groupings resulting from two standard matrices used in sequence alignment and phylogenetic tree estimation.

Algorithms↗

PubMed1994

The structure of E.coli soluble inorganic pyrophosphatase at 2.7 A resolution.

The structure of E.coli soluble inorganic pyrophosphatase has been refined at 2.7 A resolution to an R-factor of 20.9%. The overall fold of the molecule is essentially the same as yeast pyrophosphatase, except that yeast pyrophosphatase is longer at both the N- and C-termini. Escherichia coli pyrophosphatase is a mixed alpha + beta protein with a complicated topology. The active site cavity, which is also very similar to the yeast enzyme, is formed by seven beta-strands and an alpha-helix and has a rather asymmetric distribution of charged residues. Our structure-based alignment extends and improves upon earlier sequence alignment studies; it shows that probably no more than 14, not 15-17 charged and polar residues are part of the conserved enzyme mechanism of pyrophosphatases. Six of these conserved residues, at the bottom of the active site cavity, form a tight group centred on Asp70 and probably bind the two essential Mg2+ ions. The others, more spreadout and more positively charged, presumably bind substrate. Escherichia coli pyrophosphatase has an extra aspartate residue in the active site cavity, which may explain why the two enzymes bind divalent cation differently. Based on the structure, we have identified a sequence motif that seems to occur only in soluble inorganic pyrophosphatases.

Amino Acid Sequence↗

PubMed2003

Protein family annotation in a multiple alignment viewer.

SUMMARY: The Pfaat protein family alignment annotation tool is a Java-based multiple sequence alignment editor and viewer designed for protein family analysis. The application merges display features such as dendrograms, secondary and tertiary protein structure with SRS retrieval, subgroup comparison, and extensive user-annotation capabilities. AVAILABILITY: The program and source code are freely available from the authors under the GNU General Public License at http://www.pfizerdtc.com

Amino Acid Sequence↗

PubMed1993

The sequence of alignment of microtubules, focal contacts and actin filaments in fibroblasts spreading on smooth and grooved titanium substrata.

Contact guidance refers to the reactions of cells with the topography of their substratum. Current hypotheses on the mechanism of contact guidance focus on the dynamic behaviour of the cytoskeletal components, but most observations have been made on cells that have already become oriented with topographic features of the substratum. The purpose of this study was to examine the sequence in which microtubules, focal contacts and microfilament bundles become aligned to the substratum topography as fibroblasts spread on grooved substrata. Human gingival fibroblasts were trypsinized and seeded onto grooved titanium surfaces produced by micromachining, as well as onto control smooth surfaces. After observation and photography of the spreading cells at times up to 6 hours, the cells were fixed and exposed to one or more of the following antibodies or fluorescent stains: phallacidin to stain actin filaments, monoclonal anti-tubulin, monoclonal anti-vinculin, anti-mouse IgG labelled with Texas-Red or FITC, and/or an aldehyde-reactive stain to identify the cell outline. The cells were photographed and cell area, shape and orientation were calculated. Cells were also examined with confocal microscopy to obtain optical sections so that cell height as well as the precise locations of the cytoskeletal components with respect to the vertical dimension of the grooved substrata could be determined. Microtubules were the first element to become oriented parallel to the direction of the grooves and were first aligned at the bottom of the grooves. This alignment of microtubules was evident as early as 20 minutes after plating and preceded the orientation of the cell as a whole. Aligned actin microfilament bundles were not observed until 40-60 minutes and were observed first at the wall-ridge edges. At early times, focal contacts were distributed radially, but only after 3 hours did the majority of cells demonstrate aligned focal contacts. If the first cytoskeletal component to become aligned is the prime determinant of cell orientation, then these data suggest that microtubules in human gingival fibroblasts may determine cell orientation on grooved titanium surfaces. By analogy with microtubule behaviour in other systems, we suggest that microtubule orientation on grooved substrata may occur as a result of the substratum establishing shear-free planes.

Actin Cytoskeleton↗

PubMed2004

Molecular identification and chromosomal localization of genes encoding Triticum aestivum xylanase inhibitor I-like proteins in cereals.

TAXI ( Triticum aestivum xylanase inhibitor) proteins are present in wheat flour and are known to inhibit glycosyl hydrolase family 11 endoxylanases, enzymes which are commonly applied in grain processing. Here, we describe the PCR-based molecular identification of genes encoding endoxylanase inhibitors HVXI and SCXI, the TAXI-like proteins from barley ( Hordeum vulgare) and rye ( Secale cereale) respectively. The HVXI coding sequence encodes a mature protein of 384 amino acids preceded by a 19 amino acid long signal sequence. SCXI-II/III has an open reading frame encoding a signal peptide of 21 amino acids and a mature protein of 375 amino acids. As for TAXI-I, no introns were detected in the untranslated regions and coding sequences identified. These newly identified sequences allowed us to perform a multiple sequence alignment with TAXI-I and similar proteins. Rice TAXI-type proteins clustered together with the cereal endoxylanase inhibitors. Dicotyledonous proteins with sequence similarity to TAXI-I, including the tomato xyloglucan-specific endoglucanase inhibiting protein, formed a different clade. The TAXI-type proteins may hence be part of a superfamily of proteins all involved in plant responses to biotic or abiotic stress and for which a function as glycosyl hydrolase inhibitors can be suggested. The chromosomal localization of the TAXI-I gene identified on wheat chromosome 3B, of the SCXI-II/III gene identified on rye chromosome 6R, and the presence of a cluster of TAXI-like genes on rice chromosome 1, allowed us to assign the location of TAXI-like genes to the wheat-rye translocation area 3BL/6RL characterized by RFLP markers XGlb33 and Xpsr454 and isozyme Est-5. In rice, RFLP marker C1310S corresponds to a TAXI-like protein encoding sequence.

Amino Acid Sequence↗

PubMed1996

Molecular analysis of the Pseudomonas aeruginosa genes, ruvA, ruvB and ruvC, involved in processing of homologous recombination intermediates.

In Escherichia coli, the products of the ruvA, ruvB and ruvC genes are all involved in the processing of recombination intermediates (Holliday structures) into recombinant molecules. We cloned a 9.4-kb DNA fragment from Pscudomonas aeruginosa PAO1 in a plasmid by functional complementation of the UV sensitivity of an E. coli strain with ruvABC deleted. In P. aeruginosa, the ruv region seemed to form a non-SOS regulated single operon consisting of orf26-ruvC-ruvA-ruvB, while in this region of E. coli, ruvA and ruvB form an SOS-regulated operon, orf26 and ruvC form a non-SOS operon, and these two operons are split by orf23. The deduced amino acid sequences of P. aeruginosa RuvA, RuvB and RuvC proteins were 55, 72 and 55% identical to those of the corresponding E. coli Ruv proteins. The individual ruv genes of P. aeruginosa complemented the corresponding single ruv mutations of E. coli, suggesting that the P. aeruginosa Ruv proteins can interact functionally with their E. coli Ruv partners in forming heterologous complexes. The sequence alignments of the Ruv proteins were extended by incorporation of data about the putative ruv genes obtained from data banks, and the RuvB sequences were conspicuously more conserved than the RuvA and RuvC sequences.

Amino Acid Sequence↗