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

PubMed2002

CINEMA-MX: a modular multiple alignment editor.

UNLABELLED: Analyzing and visualizing multiple sequence alignments is a common task in many areas of molecular biology and bioinformatics. Many tools exist for this purpose, but are not easily customizable for specific in-house uses. Here we report the development of an editor, CINEMA-MX, that addresses these issues. CINEMA-MX is highly modular and configurable, and we present examples to illustrate its extensibility. AVAILABILITY: The program and full source code, which are available from http://www.bioinf.man.ac.uk/dbbrowser/cinema-mx, are being released under a combination of the LGPL and GPL, for Unix or Windows platforms.

Computer Graphics↗

PubMed1994

Malate dehydrogenase: a model for structure, evolution, and catalysis.

Malate dehydrogenases are widely distributed and alignment of the amino acid sequences show that the enzyme has diverged into 2 main phylogenetic groups. Multiple amino acid sequence alignments of malate dehydrogenases also show that there is a low degree of primary structural similarity, apart from in several positions crucial for nucleotide binding, catalysis, and the subunit interface. The 3-dimensional structures of several malate dehydrogenases are similar, despite their low amino acid sequence identity. The coenzyme specificity of malate dehydrogenase may be modulated by substitution of a single residue, as can the substrate specificity. The mechanism of catalysis of malate dehydrogenase is similar to that of lactate dehydrogenase, an enzyme with which it shares a similar 3-dimensional structure. Substitution of a single amino acid residue of a lactate dehydrogenase changes the enzyme specificity to that of a malate dehydrogenase, but a similar substitution in a malate dehydrogenase resulted in relaxation of the high degree of specificity for oxaloacetate. Knowledge of the 3-dimensional structures of malate and lactate dehydrogenases allows the redesign of enzymes by rational rather than random mutation and may have important commercial implications.

Amino Acid Sequence↗

PubMed1996

Molecular cloning of the mouse gene coding for carbonic anhydrase IV.

Carbonic anhydrase IV (CA IV) is expressed on apical surfaces of renal tubular epithelium and endothelium of specialized capillary beds. It plays a key role in bicarbonate reabsorption in kidney and in CO2 transport in other tissues. The human cDNA and genomic sequences have been cloned and characterized. Here we report the cloning and characterization of the entire mouse CA IV gene (contained in two overlapping lambda clones), which should enable generation of targeting constructs for disrupting the mouse CA IV gene to produce mouse models for in vivo analysis of CA IV gene function. The gene is approximately 8.2 kb long and contains eight exons ranging from 54 to 434 bp in length. The first exon (exon 1a) encodes the signal sequence. Exons 1b through 7 encode the remaining coding sequences. Exon 7 encodes the C terminus of the membrane-associated protein, as well as the 242-bp 3' untranslated sequence. The nucleotide sequence alignment between mouse and rat CA IV cDNAs reveals 84% identity. The nucleotide sequence alignment between mouse and human CA IV shows 69% identity in the coding region and all of the exon-intron boundaries are conserved, as are the sizes of the introns. The corresponding mouse and human exons are similar, except for the length of the untranslated regions in exons 1a and 7 and two small insertion/deletion events in exons 1a and 4. The 5' flanking region of the mouse gene (-300 to -1) is GC rich and contains 16 CpG dinucleotides. A TATA box sequence and several transcription factor binding sequences are identified upstream of exon 1a. Comparison of the nucleotide sequences surrounding the TATA box (-300 to -1) between mouse and human CA IV genes revealed 70% identity, indicating that regulatory sequences are as highly conserved as coding sequences between mouse and human CA IV genes.

Animals↗

PubMed1991

A predicted three-dimensional structure of human cytochrome P450: implications for substrate specificity.

A three-dimensional structure for human cytochrome P450IA1 was predicted based on the crystal coordinates of cytochrome P450cam from Pseudomonas putida. As there was only 15% residue identity between the two enzymes, additional information was used to establish an accurate sequence alignment that is a prerequisite for model building. Twelve representative eukaryotic sequences were aligned and a net prediction of secondary structure was matched against the known alpha-helices and beta-sheets of P450cam. The cam secondary structure provided a fixed main-chain framework onto which loops of appropriate length from the human P450IA1 structure were added. The model-built structure of the human cytochrome conformed to the requirements for the segregation of polar and nonpolar residues between the core and the surface. The first 44 residues of human cytochrome P450 could not be built into the model and sequence analysis suggested that residues 1-26 formed a single membrane-spanning segment. Examination of the sequences of cytochrome P450s from distinct gene families suggested specific residues that could account for the differences in substrate specificity. A major substrate for P450IA1, 3-methyl-cholanthrene, was fitted into the proposed active site and this planar aromatic molecule could be accommodated into the available cavity. Residues that are likely to interact with the haem were identified. The sequence similarity between 59 eukaryotic enzymes was represented as a dendrogram that in general clustered according to gene family. Until a crystallographic structure is available, this model-building study identifies potential residues in cytochrome P450s important in the function of these enzymes and these residues are candidates for site-directed mutagenesis.

Amino Acid Sequence↗

PubMed1992

A new approach for displaying identities and differences among aligned amino acid sequences.

An algorithm is presented for computing degrees of sequence conservation found among aligned amino acid sequences. Sequence identities are calculated for each position of an alignment and average identity values of neighboring positions are figured. The average identity value of the whole alignment is chosen as a limit to discriminate between well and less conserved sequence sections. A second algorithm is given to calculate the degree of divergence of individual sequences compared to the other sequences of the alignment. The approach is easy to use on microcomputers and gives an exact picture of sequence identities and differences in order to determine, first, protein regions of high functional or structural importance among homologous proteins, and, second, significant differences of single sequences that may contribute to individual properties of the analysed protein. The method is illustrated by an example analysing a sequence alignment of higher plant nitrate reductases.

Algorithms↗

PubMed1995

VISTAS: a package for VIsualizing STructures and sequences of proteins.

VISTAS is a suite of programs for protein sequence and structure analysis. The system allows the simultaneous display, in separate windows, of multiple sequence alignments, of known or model 3D structures, and of 2D graphic representations of sequence and/or alignment properties. The displays are fully integrated, and therefore manipulations in one window can be reflected in each of the others. Beyond its display facilities, VISTAS brings together a number of existing tools under a single, user-friendly umbrella: these include a fully functional interactive color alignment procedure, conserved motif selection, a range of database-scanning routines, and interactive access to the OWL composite sequence database and to the PRINTS protein fingerprint database. Exploration of the sequence database is thus straightforward, and predefined structural motifs from the fingerprint database may be readily visualized. Of particular note is the ability to calculate conservation criteria from sequence alignments and to display the information in a 3D context: this renders VISTAS a powerful tool for aiding mutagenesis studies and for facilitating refinement of molecular models.

Amino Acid Sequence↗

PubMed1991

A sequence assembly and editing program for efficient management of large projects.

We describe a sequence assembly and editing program for managing large and small projects. It is being used to sequence complete cosmids and has substantially reduced the time taken to process the data. In addition to handling conventionally derived sequences it can use data obtained from Applied Biosystems,Inc. 373A and Pharmacia A.L.F. fluorescent sequencing machines. Readings are assembled automatically. All editing is performed using a mouse operated contig editor that displays aligned sequences and their traces together on the screen. The editor, which can be used on single contigs or for joining contigs, permits rapid movement along the aligned sequences. Insertions, deletions and replacements can be made in individual aligned readings and global changes can be made by editing the consensus. All changes are recorded. A click on a mouse button will display the traces covering the current cursor position, hence allowing quick resolution of problems. Another function automatically moves the cursor to the next unresolved character. The editor also provides facilities for annotating the sequences. Typical annotations include flagging the positions of primers used for walking, or for marking sites, such as compressions, that have caused problems during sequencing. Graphical displays aid the assessment of progress.

Base Sequence↗

PubMed2004

Bioinformatics of the glycoside hydrolase family 57 and identification of catalytic residues in amylopullulanase from Thermococcus hydrothermalis.

Fifty-nine amino acid sequences belonging to family 57 (GH-57) of the glycoside hydrolases were collected using the CAZy server, Pfam database and blast tools. Owing to the sequence heterogeneity of the GH-57 members, sequence alignments were performed using mainly manual methods. Likewise, five conserved regions were identified, which are postulated to be GH-57 consensus motifs. In the 659 amino acid-long 4-alpha-glucanotransferase from Thermococcus litoralis, these motifs correspond to 13_HQP (region I), 76_GQLEIV (region II), 120_WLTERV (region III), 212_HDDGEKFGVW (region IV), and 350_AQCNDAYWH (region V). The third and fourth conserved regions contain the catalytic nucleophile and the proton donor, respectively. Based on our sequence alignment, residues Glu291 and Asp394 were proposed as the nucleophile and proton donor, respectively, in a GH-57 amylopullulanase from Thermococcus hydrothermalis. To validate this prediction, site-directed mutagenesis was performed. The results of this work reveal that both residues are critical for the pullulanolytic and amylolytic activities of the amylopullulanase. Therefore, these data support the prediction and strongly suggest that the bifunctionality of the amylopullulanase is determined by a single catalytic centre. Despite this positive validation, our alignment also reveals that certain GH-57 members do not possess the Glu and Asp corresponding to the predicted GH-57 catalytic residues. However, the sequences concerned by this anomaly encode putative proteins for which no biochemical or enzymatic data are yet available. Finally, the evolutionary trees generated for GH-57 reveal that the entire family can be divided into several subfamilies that may reflect the different enzyme specificities.

Amino Acid Sequence↗

PubMed2004

Computational analysis of evolution and conservation in a protein superfamily.

Many gene superfamilies have hundreds or thousands of members and hence pose a significant challenge when performing a large-scale phylogenetic analysis. Derivation of the most accurate alignment possible and inference of evolutionary relationships (with an appropriate measure of confidence) are significant "bottlenecks" in the process. A generally applicable strategy is outlined for identifying and aligning sequences, performing simple analysis of the resulting alignment, and inferring evolutionary relationships. Reference is made to the serpin superfamily. The 'partition cluster' method, a relatively rapid technique for extracting underlying associations from phylogenetic bootstrap trees, is also presented.

Computational Biology↗

PubMed2005

A simple tool to explore the distance distribution of correlated mutations in proteins.

The analysis of correlated mutations in protein sequence alignments is of considerable interest, since it may provide useful energetic and even structural information (ideally, residue contacts). However, a number of recent experimental studies support the existence of long-distance communication in proteins, a fact that may lead to correlation between distant residues. We introduce in this work a simple statistical procedure to describe the relation structure--alignments on the basis of the residue--residue distance dependence of the number of residue couples over given thresholds of a correlation measure (such as a covariance value). This procedure may lead to clear pictures of the distance distribution of correlated mutations and may provide a simple but efficient tool to explore the different structural features that are reflected in the sequence alignments.

Amino Acids↗

PubMed2006

Sequence representation and prediction of protein secondary structure for structural motifs in twilight zone proteins.

Characterizing and classifying regularities in protein structure is an important element in uncovering the mechanisms that regulate protein structure, function and evolution. Recent research concentrates on analysis of structural motifs that can be used to describe larger, fold-sized structures based on homologous primary sequences. At the same time, accuracy of secondary protein structure prediction based on multiple sequence alignment drops significantly when low homology (twilight zone) sequences are considered. To this end, this paper addresses a problem of providing an alternative sequences representation that would improve ability to distinguish secondary structures for the twilight zone sequences without using alignment. We consider a novel classification problem, in which, structural motifs, referred to as structural fragments (SFs) are defined as uniform strand, helix and coil fragments. Classification of SFs allows to design novel sequence representations, and to investigate which other factors and prediction algorithms may result in the improved discrimination. Comprehensive experimental results show that statistically significant improvement in classification accuracy can be achieved by: (1) improving sequence representations, and (2) removing possible noise on the terminal residues in the SFs. Combining these two approaches reduces the error rate on average by 15% when compared to classification using standard representation and noisy information on the terminal residues, bringing the classification accuracy to over 70%. Finally, we show that certain prediction algorithms, such as neural networks and boosted decision trees, are superior to other algorithms.

Algorithms↗

PubMed2000

Characterization of the structure and expression of a highly conserved ribosomal protein gene, L9, from pea.

The eukaryotic ribosomal protein (RP) L9 is highly conserved in nature, and its gene is expressed to high levels in the actively growing tissues of pea. The transcriptional activity of the gene is highest in root, cambial and shoot meristems and immature tissues of the plant. Promoter deletion analysis using constructs employing the reporter gene gus were stably transferred into tobacco and revealed that the fully functional promoter is found in the first 316bp upstream from the start codon. Transgenic pea plants carrying one of these constructs show that translational efficiency mirrors gene transcription; gene expression appears to be developmentally regulated at the level of transcription. The coding region of the gene shares 80% amino acid homology with Arabidopsis and 76% homology with rice. Comparisons of the gene structure to that of the human, fruit fly, yeast, and Arabidopsis homologues reveal a close relationship in both promoter structure and intron insertion sites with the Arabidopsis gene. A nucleotide sequence alignment of the pea gene with other plant RP genes revealed that a sequence, -TTAGGGTTTT-, was commonly found in the forward and/or the inverted orientation at or near the TATA boxes of the promoters of these genes and may have a role in regulating the coordinate production of the RP genes in plants.

Amino Acid Sequence↗

PubMed2003

Cloning and prokaryotic expression of a cDNA encoding a putative mitochondrial malate dehydrogenase in Oryza sativa.

Malate dehydrogenase (MDH) has been characterized as a key player in oxaloacetate (OAA) biosynthesis mechanism in citrate acid cycle that generates reducing powers for further assimilation in the whole cell. Here we present the cloning, characterization and prokaryotic expression of a putative Mdh (OsmMDH) in Oryza sativa. Sequence alignment shows that there is a high homology between the deduced amino acid sequence of OsmMDH and MDH portein in Eucalyptus gunnii (80%), as well as between the deduced amino acid sequence of OsmMDH and other MDHs. Moreover, pI and the mitochondrial location of OsmMDH are predicted. The tissue-specific expression pattern of OsmMDH reveals that it is abundant in young panicle and immature seed, while its expression level is mush lower in leaf and root. Its expression in E. coli BL21 as a fusion gene is studied further.

Amino Acid Sequence↗

PubMed2004

MyHits: a new interactive resource for protein annotation and domain identification.

The MyHits web server (http://myhits.isb-sib.ch) is a new integrated service dedicated to the annotation of protein sequences and to the analysis of their domains and signatures. Guest users can use the system anonymously, with full access to (i) standard bioinformatics programs (e.g. PSI-BLAST, ClustalW, T-Coffee, Jalview); (ii) a large number of protein sequence databases, including standard (Swiss-Prot, TrEMBL) and locally developed databases (splice variants); (iii) databases of protein motifs (Prosite, Interpro); (iv) a precomputed list of matches ('hits') between the sequence and motif databases. All databases are updated on a weekly basis and the hit list is kept up to date incrementally. The MyHits server also includes a new collection of tools to generate graphical representations of pairwise and multiple sequence alignments including their annotated features. Free registration enables users to upload their own sequences and motifs to private databases. These are then made available through the same web interface and the same set of analytical tools. Registered users can manage their own sequences and annotations using only web tools and freeze their data in their private database for publication purposes.

Computer Graphics↗

PubMed1999

Molecular modeling of the amyloid-beta-peptide using the homology to a fragment of triosephosphate isomerase that forms amyloid in vitro.

The main component of the amyloid senile plaques found in Alzheimer's brain is the amyloid-beta-peptide (A beta), a proteolytic product of a membrane precursor protein. Previous structural studies have found different conformations for the A beta peptide depending on the solvent and pH used. In general, they have suggested an alpha-helix conformation at the N-terminal domain and a beta-sheet conformation for the C-terminal domain. The structure of the complete A beta peptide (residues 1-40) solved by NMR has revealed that only helical structure is present in A beta. However, this result cannot explain the large beta-sheet A beta aggregates known to form amyloid under physiological conditions. Therefore, we investigated the structure of A beta by molecular modeling based on extensive homology using the Smith and Waterman algorithm implemented in the MPsrch program (Blitz server). The results showed a mean value of 23% identity with selected sequences. Since these values do not allow a clear homology to be established with a reference structure in order to perform molecular modeling studies, we searched for detailed homology. A 28% identity with an alpha/beta segment of a triosephosphate isomerase (TIM) from Culex tarralis with an unsolved three-dimensional structure was obtained. Then, multiple sequence alignment was performed considering A beta, TIM from C.tarralis and another five TIM sequences with known three-dimensional structures. We found a TIM segment with secondary structure elements in agreement with previous experimental data for A beta. Moreover, when a synthetic peptide from this TIM segment was studied in vitro, it was able to aggregate and to form amyloid fibrils, as established by Congo red binding and electron microscopy. The A beta model obtained was optimized by molecular dynamics considering ionizable side chains in order to simulate A beta in a neutral pH environment. We report here the structural implications of this study.

Algorithms↗

PubMed2005

Genetic variability of the heme uptake system among different strains of the fish pathogen Vibrio anguillarum: identification of a new heme receptor.

The ability to utilize heme compounds as iron sources was investigated in Vibrio anguillarum strains belonging to serotypes O1 to O10. All strains, regardless of their serotype or isolation origin could utilize hemin and hemoglobin as sole iron sources. Similarly, all of the isolates could bind hemin and Congo red, and this binding was mediated by cell envelope proteins. PCR and Southern hybridization were used to assay the occurrence of heme transport genes huvABCD, which have been previously described in serotype O1. Of 23 strains studied, two serotype O3 isolates proved negative for all huvABCD genes, whereas nine strains included in serotypes O2, O3, O4, O6, O7, and O10 tested negative for the outer membrane heme receptor gene huvA. A gene coding for a novel outer membrane heme receptor was cloned and characterized in a V. anguillarum serotype O3 strain lacking huvA. The new heme receptor, named HuvS, showed significant similarity to other outer membrane heme receptors described in Vibrionaceae, but little homology (39%) to HuvA. This heme receptor was present in 9 out of 11 of the V. anguillarum strains that tested negative for HuvA. Furthermore, complementation experiments demonstrated that HuvS could substitute for the HuvA function in Escherichia coli and V. anguillarum mutants. The huvS and huvA sequences alignment, as well as the analysis of their respective upstream and downstream DNA sequences, suggest that horizontal transfer and recombination might be responsible for generating this genetic diversity.

Amino Acid Sequence↗

PubMedDate not supplied

Synovial T cell receptor heterogeneity in early arthritis.

Rheumatoid arthritis (RA) has been postulated to result from a synovial immune response to an unidentified antigen(s), which should be mirrored by the T cell response. Here we investigate the T cell receptor (TCR) repertoire in the synovial tissue of patients with arthritis of early to moderate duration. We developed a nested polymerase chain reaction (PCR) technique to examine the TCR repertoire of small biopsy specimens, and show that the method is highly sensitive. We apply this technique to synovial biopsies obtained from the knee joints of patients with early to moderate duration arthritis (average duration of arthritis 1 year, range 0. 02-2.75 years). We examined biopsies from 5 normal individuals, 32 RA patients, 7 patients with seronegative spondyloarthropathy (Sp), and 12 patients with undifferentiated arthritis (UA). TCR message was detectable in 4/5 normals, 15/32 RA, 5/7 Sp, and 8/11 UA biopsies, with sampling error likely accounting for most negative biopsies. The average numbers of TCR Vbetas detected per TCR-positive biopsy were 5.0 +/- 3.7 for normals, 12.7 +/- 8.4 for RA, 18.0 +/- 7.4 for Sp patients, and 14.4 +/- 10.2 for UA. Examination of TCR messages by single-stranded conformational polymorphism analysis showed similar proportions of dominant clones in the normals compared with the patients with inflammatory arthritis. Sequence analysis was performed on 33 dominant clones from 16 patients. Sequence alignment of the third hypervariable regions showed some evidence of disease-specific sequence clustering for Sp, while some RA sequences showed similarity to previously described motifs. These data indicate greater TCR heterogeneity in early Sp and UA compared with normal synovium. Disease-specific TCR sequences may occur in early RA and Sp.

Adult↗

PubMed2005

Phylogenetic analysis of family 6 glycoside hydrolases.

Multiple sequence alignment separates members of glycoside hydrolase Family 6 into eight subfamilies: one of mainly actinobacterial endoglucanases (EGs), one of ascomycotal EGs, one of chytridiomycotal EGs and cellobiohydrolases (CBHs), one of actinobacterial and proteobacterial CBHs, one of chytridiomycotal CBHs, two of ascomycotal CBHs, and one of basidiomycotal CBHs. Each also has some proteins of unknown function. Multiple sequence alignment also extends to all of Family 6 the observation that lengths of loops that form the active-site tunnel in CBHs vary among subfamilies, and along with loop conformations, determine enzyme function.

Amino Acid Sequence↗