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

PubMed2002

The Celera Discovery System.

The Celera Discovery System (CDS) is a web-accessible research workbench for mining genomic and related biological information. Users have access to the human and mouse genome sequences with annotation presented in summary form in BioMolecule Reports for genes, transcripts and proteins. Over 40 additional databases are available, including sequence, mapping, mutation, genetic variation, mRNA expression, protein structure, motif and classification data. Data are accessible by browsing reports, through a variety of interactive graphical viewers, and by advanced query capability provided by the LION SRS search engine. A growing number of sequence analysis tools are available, including sequence similarity, pattern searching, multiple sequence alignment and Hidden Markov Model search. A user workspace keeps track of queries and analyses. CDS is widely used by the academic research community and requires a subscription for access. The system and academic pricing information are available at http://cds.celera.com.

Animals↗

PubMed2005

The ubiquitous protein domain EAL is a cyclic diguanylate-specific phosphodiesterase: enzymatically active and inactive EAL domains.

The EAL domain (also known as domain of unknown function 2 or DUF2) is a ubiquitous signal transduction protein domain in the Bacteria. Its involvement in hydrolysis of the novel second messenger cyclic dimeric GMP (c-di-GMP) was demonstrated in vivo but not in vitro. The EAL domain-containing protein Dos from Escherichia coli was reported to hydrolyze cyclic AMP (cAMP), implying that EAL domains have different substrate specificities. To investigate the biochemical activity of EAL, the E. coli EAL domain-containing protein YahA and its individual EAL domain were overexpressed, purified, and characterized in vitro. Both full-length YahA and the EAL domain hydrolyzed c-di-GMP into linear dimeric GMP, providing the first biochemical evidence that the EAL domain is sufficient for phosphodiesterase activity. This activity was c-di-GMP specific, optimal at alkaline pH, dependent on Mg(2+) or Mn(2+), strongly inhibited by Ca(2+), and independent of protein oligomerization. Linear dimeric GMP was shown to be 5'pGpG. The EAL domain from Dos was overexpressed, purified, and found to function as a c-di-GMP-specific phosphodiesterase, not as a cAMP-specific phosphodiesterase, in contrast to previous reports. The EAL domains can hydrolyze 5'pGpG into GMP, however, very slowly, thus implying that this activity is irrelevant in vivo. Therefore, c-di-GMP is the exclusive substrate of EAL. Multiple-sequence alignment revealed two groups of EAL domains hypothesized to correspond to enzymatically active and inactive domains. The domains in the latter group have mutations in residues conserved in the active domains. The enzymatic inactivity of EAL domains may explain their coexistence with GGDEF domains in proteins possessing c-di-GMP synthase (diguanulate cyclase) activity.

3',5'-Cyclic-GMP Phosphodiesterases↗

PubMed1995

Ribosomal protein L22 from Thermus thermophilus: sequencing, overexpression and crystallisation.

The gene for the ribosomal protein L22 from Thermus thermophilus has been sequenced and overexpressed in Escherichia coli. A multiple sequence alignment was carried out for all proteins of the L22 family reported so far. The recombinant protein was purified and crystallized. The crystals belong to the space group P2(1)2(1)2(1), with cell parameters of a = 32.6 A, b = 66.0 A, c = 67.8 A.

Amino Acid Sequence↗

PubMed2003

Molecular comparison of rat cyclic nucleotide phosphodiesterase 8 family: unique expression of PDE8B in rat brain.

Cyclic nucleotide phosphodiestease (PDE) type 8 is categorized into a family of 3-isobutyl-1-methylxanthine-insensitive PDE hydrolyzing cAMP with high affinity. We have isolated cDNAs encoding rat PDE8A and PDE8B from brain and testis, respectively. The sequence analysis demonstrated that rat PDE8A was a protein of 823 amino acid residues. Rat PDE8B protein was predicted as an N-terminal truncated form of 760 amino acid residues. Both of rat PDE8 proteins include REC, PAS and catalytic PDE domains. Tissue-specific expression patterns of rat PDE8A and PDE8B transcripts were demonstrated by Northern blot analysis. Rat PDE8A transcripts were rich in the liver and testis, and those of rat PDE8B were particularly abundant in the brain and were not expressed in the thyroid gland, while human thyroid gland contains PDE8B transcripts at a high level. Rat PDE8B transcripts were found in all brain regions other than cerebellum and shown to exist in the neuronal cells in in situ hybridization. Mouse PDE8B1 sequence was also identified by a database search and sequence alignment, revealing a protein of 885 amino acid residues, which is 99% and 96% identical to rat and human PDE8B1, respectively. As well as rat PDE8B, expression of mouse PDE8B transcripts was not confined to the thyroid gland. Species-dependent tissue expression pattern was quite unique features of PDE8B.

3',5'-Cyclic-AMP Phosphodiesterases↗

PubMed2003

Design and assessment of a fast algorithm for identifying specific probes for human and mouse genes.

MOTIVATION: Mammalian genomes are highly complex. To identify the unique sequences of each gene in a mammalian gene database containing tens of thousands of DNA sequences is a computation intensive task. With the advent of parallel genetic analysis methods such as microarrays and the availability of more and more whole genome sequences of organisms, an algorithm allowing speedy identification of the unique gene probes for functional studies of individual genes will be a very useful tool. RESULTS: We have developed a fast algorithm as well as a software program based on the algorithm for identifying gene specific probes of complex organisms. The algorithm was applied to the assemblies of gene sequences and was highly efficient for large databases such as the TIGR human THC and mouse TC databases. The results were assessed with the BLAST sequence alignment software. Two probe data sets have been compiled to contain specific probes for around 100 000 putative human gene transcripts and 70 000 putative mouse gene transcripts. AVAILABILITY: The gene specific probes for the putative human and mouse genes referenced in the TIGR gene indices are available at: ftp://genestamp.ibms.sinica.edu.tw/pub/SpecificP/. The software program and the source codes are available upon request.

Algorithms↗

PubMed2006

Truncated RIP3 (tRIP3) acts upstream of FADD to induce apoptosis in the human hepatocellular carcinoma cell line QGY-7703.

RIP3 (receptor-interacting protein 3) is a serine/threonine kinase that promotes apoptosis in various cell types. The C-terminal domain of RIP3 is critical for its apoptosis induction. In this study, we showed that a truncated form of RIP3 (tRIP3) containing only the unique C-terminal region (aa 224-518) induced significant apoptosis in human hepatocellular carcinoma cells QGY-7703. A FADD-dominant negative (FADD-DN) was shown to significantly block apoptosis induced by tRIP3. In contrast, the RIP3 dominant negative (RIP3-DN) was found unable to block FADD-induced apoptosis. Thus, we conclude that tRIP3 may function upstream of FADD to induce apoptosis in TNFR-1 signaling pathway. Additionally, sequence alignments of RIP3 with other death domain (DD)-containing proteins revealed that amino acids Leu 477 and Leu 488 in RIP3 are highly conserved, nonetheless, neither mutational change at Leu 477 nor at Leu 488 confers obvious effect on cell death, indicating that these two amino acids might not be critical for its pro-apoptotic activity as expected.

Adaptor Proteins, Signal Transducing↗

PubMed2006

Molecular characterization of a HMW glutenin subunit allele providing evidence for silencing of x-type gene on Glu-B1.

Understanding the molecular structure of high-molecular-weight glutenin subunit (HMW-GS) may provide useful evidence for the study on the improvement of quality of cultivated wheat and the evolution of Glu-1 alleles. Sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) shows that the subunits encoded by Glu-B1 were null, named 1Bxm, in a Triticum turgidum var. dicoccoides line PI94640. Primers based on the conserved regions in wheat HMW-GS gene promoter and coding sequences were used to amplify the genomic DNA of line PI94640. The PCR products were sequenced, and the total nucleotide sequence of 3,442 bp including upstream sequence of 1,070 bp was obtained. Compared with the reported gene sequences of Glu-1Bx alleles, the promoter region of the Glu-1Bxm showed close resemblance to 1Bx7. The Glu-1Bxm coding region differs from the other Glu-1Bx alleles for a deduced mature protein with only 212 residues, and a stop codon (TAA) at 637 bp downstream from the start codon was present, which was probably responsible for the silencing of x-type subunit genes at the Glu-B1 locus. Phylogenetic tree based on the nucleotide sequence alignment of HMW glutenin subunit genes showed that 1Bxm was the most ancient type of Glu-B1 alleles, suggesting that the evolution rates are different among Glu-1Bx genes. Further study on the contribution of the unique silenced Glu-B1 alleles to quality improvement was also discussed.

Alleles↗

PubMed2006

XRate: a fast prototyping, training and annotation tool for phylo-grammars.

BACKGROUND: Recent years have seen the emergence of genome annotation methods based on the phylo-grammar, a probabilistic model combining continuous-time Markov chains and stochastic grammars. Previously, phylo-grammars have required considerable effort to implement, limiting their adoption by computational biologists. RESULTS: We have developed an open source software tool, xrate, for working with reversible, irreversible or parametric substitution models combined with stochastic context-free grammars. xrate efficiently estimates maximum-likelihood parameters and phylogenetic trees using a novel "phylo-EM" algorithm that we describe. The grammar is specified in an external configuration file, allowing users to design new grammars, estimate rate parameters from training data and annotate multiple sequence alignments without the need to recompile code from source. We have used xrate to measure codon substitution rates and predict protein and RNA secondary structures. CONCLUSION: Our results demonstrate that xrate estimates biologically meaningful rates and makes predictions whose accuracy is comparable to that of more specialized tools.

Algorithms↗

PubMed2002

Structural aspects of aldehyde dehydrogenase that influence dimer-tetramer formation.

Aldehyde dehydrogenases are isolated as dimers or tetramers but have essentially identical structures. The homotetramer (ALDH1 or ALDH2) is a dimer of dimers (A-B + C-D). In the tetrameric enzyme, Ser500 from subunit "D" interacts with Arg84, a conserved residue, from subunit "A". In the dimeric ALDH3 form, the interaction cannot exist. It has been proposed that the formation of the tetramer is prevented by the presence of a C-terminal tail in ALDH3 that is not present in ALDH1 or 2. To understand the forces that maintain the tetramer, deletion of the tail in ALDH3, addition of different tails in ALDH1, and mutations of different residues located in the dimer-dimer interface were made. Gel filtration of the recombinantly expressed enzymes demonstrated that no change in their oligomerization occurred. Urea denaturation showed a diminution to the stability of the ALDH1 mutants. The K(m) for propionaldehyde was similar to that of the wild-type enzyme, but the K(m) for NAD was altered. A double mutant of D80G and S82A produced an enzyme with the ability to form dimers and tetramers in a protein-concentration-dependent manner. Though stable, this dimeric form was inactive. The tetramer exhibited 10% activity compared with the wild type. Sequence alignment demonstrated that the hydrophobic surface area is increased in the tetrameric enzymes. The hydrophobic surface seems to be the main force that drives the formation of tetramers. The results indicated that residues 80 and 82 are involved in maintaining the tetramer after its assembly but the C-terminal extension contributes to the overall stability of the assembled protein.

Aldehyde Dehydrogenase↗

PubMedDate not supplied

Identifying DNA and protein patterns with statistically significant alignments of multiple sequences.

MOTIVATION: Molecular biologists frequently can obtain interesting insight by aligning a set of related DNA, RNA or protein sequences. Such alignments can be used to determine either evolutionary or functional relationships. Our interest is in identifying functional relationships. Unless the sequences are very similar, it is necessary to have a specific strategy for measuring-or scoring-the relatedness of the aligned sequences. If the alignment is not known, one can be determined by finding an alignment that optimizes the scoring scheme. RESULTS: We describe four components to our approach for determining alignments of multiple sequences. First, we review a log-likelihood scoring scheme we call information content. Second, we describe two methods for estimating the P value of an individual information content score: (i) a method that combines a technique from large-deviation statistics with numerical calculations; (ii) a method that is exclusively numerical. Third, we describe how we count the number of possible alignments given the overall amount of sequence data. This count is multiplied by the P value to determine the expected frequency of an information content score and, thus, the statistical significance of the corresponding alignment. Statistical significance can be used to compare alignments having differing widths and containing differing numbers of sequences. Fourth, we describe a greedy algorithm for determining alignments of functionally related sequences. Finally, we test the accuracy of our P value calculations, and give an example of using our algorithm to identify binding sites for the Escherichia coli CRP protein. AVAILABILITY: Programs were developed under the UNIX operating system and are available by anonymous ftp from ftp://beagle.colorado.edu/pub/consensus.

Algorithms↗

PubMed1997

Identification of proteolipid from an extremely halophilic archaeon Halobacterium salinarum as an N,N'-dicyclohexyl-carbodiimide binding subunit of ATP synthase.

ATP synthesis in an extremely halophilic archaeon, Halobacterium salinarum, was inhibited by N-cyclohexyl-N'-[4-(dimethylamino)-alpha-naphthyl]carbodiimide (NCD-4), a fluorescent analog of N,N'-dicyclohexylcarbodiimide (DCCD). By tracing the fluorescent signal, a hydrophobic 8-kDa protein (proteolipid) was purified from the halobacterial membrane as one of the most DCCD-reactive proteins and its N-terminal amino acid sequence was determined. The gene encoding the proteolipid was found in the region upstream of the genes encoding the two major subunits of halobacterial A-type ATPase [K. Ihara and Y.Mukohata (1991) Arch. Biochem. Biophys. 286, 111-116]. Halobacterial proteolipid was more similar in size to the proteolipid of F-type ATPase than that of V-type ATPase. However, multiple amino acid sequence alignment of proteolipids showed a higher degree of relatedness between V-type and A-type ATPase proteolipids. Together with the recent finding of a triplicate proteolipid encoding gene from the methanogenic archaeon Methanococcus jannaschii [C. J. Bult et al. (1996) Science 273, 1058-1073], proteolipids from archaea seem to have diverse characteristics in comparison with those from eubacteria or from eukaryotes.

Adenosine Triphosphatases↗

PubMed1999

Identification of the RNA-binding sites of the triple gene block protein 1 of bamboo mosaic potexvirus.

The triple gene block protein 1 (TGBp1) encoded by open reading frame 2 of bamboo mosaic potexvirus (BaMV) was overexpressed in Escherichia coli and purified in order to test its RNA-binding activity. UV crosslinking assays revealed that the RNA-binding activity was present mainly in the soluble fraction of the refolded TGBp1. The binding activity was nonspecific and salt concentration-dependent: activity was present at 0-50 mM NaCl but was almost abolished at 200 mM. The RNA-binding domain was located by deletion mutagenesis to the N-terminal 3-24 amino acids of TGBp1. Sequence alignment analysis of the N-terminal 25 amino acids of the TGBp1 homologues of potexviruses identified three arginine residues. Arg-to-Ala substitution at any one of the three arginines eliminated most of the RNA-binding activity, indicating that they were all critical to the RNA-binding activity of the TGBp1 of BaMV.

Amino Acid Sequence↗

PubMed2004

SWAPSC: sliding window analysis procedure to detect selective constraints.

UNLABELLED: Sliding-window analysis procedure to detect selective constraints (SWAPSC) is a software system to dissect the constraints on the evolution of protein-coding genes. The program estimates rates of nucleotide substitutions at specific codon regions in each branch of a phylogenetic tree. The program uses several sets of simulated sequence alignments to estimate the probability of synonymous and non-synonymous nucleotide substitutions. Thereafter, a statistical analysis is conducted to determine the optimum window size to detect selective constraints. Finally, the optimum window size is slid along the real alignment and a test for significance of the estimated number of synonymous and non-synonymous nucleotide substitutions in each sliding step is conducted. A number of friendly useful output files is generated. AVAILABILITY: SWAPSC is available at http://www.may.ie/academic/biology/staff/mfmolecevolandbioinf.shtml distribution versions for both Linux and Windows operating systems are available, including manual and example files.

Algorithms↗

PubMed2006

Fast de novo peptide sequencing and spectral alignment via tree decomposition.

De novo sequencing and spectral alignment are computationally important for the prediction of new protein peptides via tandem mass spectrometry (MS/MS). Both approaches are established upon the problem of finding the longest antisymmetric path on formulated graphs. The problem is of high computational complexity and the prediction accuracy is compromised when given spectra involve noisy data, missing mass peaks, or post translational modifications (PTMs) and mutations. This paper introduces a graphical mechanism to describe relationships among mass peaks that, through graph tree decomposition, yields linear and quadratic time algorithms for optimal de novo sequencing and spectral alignment respectively. Our test results show that, in addition to high efficiency, the new algorithms can achieve desired prediction accuracy on spectra containing noisy peaks and PTMs while allowing the presence of both b-ions and y-ions.

Algorithms↗

PubMed1999

The evolution of starch-binding domain.

Amylolytic enzymes belonging to three distinct families of glycosidases (13, 14, 15) contain the starch-binding domain (SBD) positioned almost exclusively at the C-terminus. Detailed analysis of all available SBD sequences from 43 different amylases revealed its independent evolutionary behaviour with regard to the catalytic domains. In the evolutionary tree based on sequence alignment of the SBDs, taxonomy is respected so that fungi and actinomycetes form their own separate parts surrounded by bacteria that are also clustered according to taxonomy. The only known N-terminal SBD from Rhizopus oryzae glucoamylase is on the longest branch separated from all C-terminal SBDs. The 3-dimensional (3-D) structures of fungal glucoamylase and bacterial CGTase SBDs are compared and used to discuss the interesting SBD evolution.

Amino Acid Sequence↗

PubMed1997

cDNA cloning, tissue distribution, and identification of the catalytic triad of monoglyceride lipase. Evolutionary relationship to esterases, lysophospholipases, and haloperoxidases.

Monoglyceride lipase catalyzes the last step in the hydrolysis of stored triglycerides in the adipocyte and presumably also complements the action of lipoprotein lipase in degrading triglycerides from chylomicrons and very low density lipoproteins. Monoglyceride lipase was cloned from a mouse adipocyte cDNA library. The predicted amino acid sequence consisted of 302 amino acids, corresponding to a molecular weight of 33,218. The sequence showed no extensive homology to other known mammalian proteins, but a number of microbial proteins, including two bacterial lysophospholipases and a family of haloperoxidases, were found to be distantly related to this enzyme. By means of multiple sequence alignment and secondary structure prediction, the structural elements in monoglyceride lipase, as well as the putative catalytic triad, were identified. The residues of the proposed triad, Ser-122, in a GXSXG motif, Asp-239, and His-269, were confirmed by site-directed mutagenesis experiments. Northern blot analysis revealed that monoglyceride lipase is ubiquitously expressed among tissues, with a transcript size of about 4 kilobases.

Adipocytes↗

PubMed2006

The MIGenAS integrated bioinformatics toolkit for web-based sequence analysis.

We describe a versatile and extensible integrated bioinformatics toolkit for the analysis of biological sequences over the Internet. The web portal offers convenient interactive access to a growing pool of chainable bioinformatics software tools and databases that are centrally installed and maintained by the RZG. Currently, supported tasks comprise sequence similarity searches in public or user-supplied databases, computation and validation of multiple sequence alignments, phylogenetic analysis and protein-structure prediction. Individual tools can be seamlessly chained into pipelines allowing the user to conveniently process complex workflows without the necessity to take care of any format conversions or tedious parsing of intermediate results. The toolkit is part of the Max-Planck Integrated Gene Analysis System (MIGenAS) of the Max Planck Society available at www.migenas.org (click 'Start Toolkit').

Animals↗

PubMed1998

An enzymatically active chimeric HIV-1 reverse transcriptase (RT) with the RNase-H domain of murine leukemia virus RT exists as a monomer.

The existence of retroviral reverse transcriptases as monomers or dimers is rather intriguing. A classical example of the former is murine leukemia virus reverse transcriptase (MuLV RT), while human immunodeficiency virus type 1 (HIV-1) RT represents the latter. A careful scrutiny of the amino acid sequence alignment of the two enzymes pinpoints the region tentatively responsible for this phenomenon. We report here the construction of a chimeric enzyme containing the first 425 amino acid residues from the N-terminal domain of HIV-1 RT and 200 amino acid residues from the C-terminal domain of MuLV RT. The chimeric enzyme exists as a monomer with intact DNA polymerase and RNase-H functions.

Amino Acid Sequence↗