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

PubMed1991

Molecular cloning of three cDNAs that encode cysteine proteinases in the digestive gland of the American lobster (Homarus americanus).

Three clones were isolated from a lobster digestive gland cDNA library, using oligonucleotide probes based on the partial amino terminal sequence of a digestive cysteine proteinase. The cDNAs, LCP1, LCP2 and LCP3 encode preproenzymes of 322, 323 and 321 amino acid residues, and putative mature enzymes of 217, 216 and 215 residues, respectively. Calculated mature protein molecular masses are 23386 (LCP1), 29093 (LCP2) and 23255 (LCP3) Sequence alignments show that the lobster enzymes are more similar to L (55-62% identity) than H (42-44%) or B (22-24%) cathepsins. Southern analysis indicated as many as eleven genes related to the three cDNAs.

Animals↗

PubMed2004

Reconstruction of ancestral protein sequences and its applications.

BACKGROUND: Modern-day proteins were selected during long evolutionary history as descendants of ancient life forms. In silico reconstruction of such ancestral protein sequences facilitates our understanding of evolutionary processes, protein classification and biological function. Additionally, reconstructed ancestral protein sequences could serve to fill in sequence space thus aiding remote homology inference. RESULTS: We developed ANCESCON, a package for distance-based phylogenetic inference and reconstruction of ancestral protein sequences that takes into account the observed variation of evolutionary rates between positions that more precisely describes the evolution of protein families. To improve the accuracy of evolutionary distance estimation and ancestral sequence reconstruction, two approaches are proposed to estimate position-specific evolutionary rates. Comparisons show that at large evolutionary distances our method gives more accurate ancestral sequence reconstruction than PAML, PHYLIP and PAUP*. We apply the reconstructed ancestral sequences to homology inference and functional site prediction. We show that the usage of hypothetical ancestors together with the present day sequences improves profile-based sequence similarity searches; and that ancestral sequence reconstruction methods can be used to predict positions with functional specificity. CONCLUSIONS: As a computational tool to reconstruct ancestral protein sequences from a given multiple sequence alignment, ANCESCON shows high accuracy in tests and helps detection of remote homologs and prediction of functional sites. ANCESCON is freely available for non-commercial use. Pre-compiled versions for several platforms can be downloaded from ftp://iole.swmed.edu/pub/ANCESCON/.

Amino Acid Sequence↗

PubMed1994

The xylose isomerase-encoding gene (xylA) of Clostridium thermosaccharolyticum: cloning, sequencing and phylogeny of XylA enzymes.

The xylose isomerase (XylA)-encoding gene (xylA) of the thermophilic anaerobic bacterium, Clostridium thermosaccharolyticum NCIB 9385, was cloned as a 4.0-kb DNA fragment by complementation of the Escherichia coli xylA mutant strain, DS941. The open reading frame of 1317 bp encoded a protein of 439 amino acids (aa), with a calculated M(r) of 50,236. The gene was preceeded by a typical clostridial Shine-Dalgarno sequence, and was expressed constitutively in the cloning host. Downstream, the clone appeared to carry a xylB gene (encoding xylulokinase) in the same orientation as xylA. Comparison of the deduced aa sequence of the C. thermosaccharolyticum XylA with 18 other XylA showed that this family of proteins was separated into two clusters, one comprising proteins from organisms with G + C-rich DNA, and the other proteins from organisms with a lower G + C composition. Within the second cluster, the XylA of C. thermosaccharolyticum was most closely related to the enzymes from C. thermosulfurogenes (Thermoanaerobacterium thermosulfurigenes) and C. thermohydrosulfuricum (93 and 84% identity, respectively). Analysis of the aligned sequences indicated two signatures (VXW[GP]GREG[YSTA]E and [LIVM]EPKPX]EQ]P) which may be useful in isolation of novel XylA.

Aldose-Ketose Isomerases↗

PubMed2005

Structure of adeno-associated virus type 4.

Adeno-associated virus (AAV) is a member of the Parvoviridae, belonging to the Dependovirus genus. Currently, several distinct isolates of AAV are in development for use in human gene therapy applications due to their ability to transduce different target cells. The need to manipulate AAV capsids for specific tissue delivery has generated interest in understanding their capsid structures. The structure of AAV type 4 (AAV4), one of the most antigenically distinct serotypes, was determined to 13-A resolution by cryo-electron microscopy and image reconstruction. A pseudoatomic model was built for the AAV4 capsid by use of a structure-based sequence alignment of its major capsid protein, VP3, with that of AAV2, to which AAV4 is 58% identical and constrained by its reconstructed density envelope. The model showed variations in the surface loops that may account for the differences in receptor binding and antigenicity between AAV2 and AAV4. The AAV4 capsid surface topology also shows an unpredicted structural similarity to that of Aleutian mink disease virus and human parvovirus B19, autonomous members of the genus, despite limited sequence homology.

Amino Acid Sequence↗

PubMed1995

Structural features of a superfamily of zinc-endopeptidases: the metzincins.

A large number of zinc endopeptidases contain an HEXXHXXGXXH consensus motif in their catalytic site (single letter code; X is any amino acid residue). These enzymes can be grouped into four distinct families, the astacins, the adamalysins, the serralysins and the matrix metalloproteinases (matrixins). Despite a low degree of sequence similarity, their catalytic modules are topologically similar. A topology derived sequence alignment suggests that the four families form a superfamily, called the metzincins because of a perfectly superimposable methionine residue close to the zinc-binding active site. Topological similarity to the thermolysin-like enzymes indicates that these enzymes may have had a common ancestor.

Amino Acid Sequence↗

PubMed1998

Seeking an ancient enzyme in Methanococcus jannaschii using ORF, a program based on predicted secondary structure comparisons.

We have developed a simple procedure to identify protein homologs in genomic databases. The program, called ORF, is based on comparisons of predicted secondary structure. Protein structure is far better conserved than amino acid sequence, and structure-based methods have been effective in exploiting this fact to find homologs, even among proteins with scant sequence identity. ORF is a secondary structure-based method that operates solely on predictions from sequence and requires no experimentally determined information about the structure. The approach is illustrated by an example: Thymidylate synthase, a highly conserved enzyme essential to thymidine biosynthesis in both prokaryotes and eukaryotes, is thought to be used by Archaea, but a corresponding gene has yet to be identified. Here, a candidate thymidylate synthase is identified as a previously unassigned open reading frame from the genome of Methanococcus jannaschii, viz., MJ0757. Using primary structure information alone, the optimally aligned sequence identity between MJ0757 and Escherichia coli thymidylate synthase is 7%, well below the threshold of sensitivity for detection by sequence-based methods.

Amino Acid Sequence↗

PubMed1989

A new method that simultaneously aligns and reconstructs ancestral sequences for any number of homologous sequences, when the phylogeny is given.

Among the fundamental problems in molecular evolution and in the analysis of homologous sequences are alignment, phylogeny reconstruction, and the reconstruction of ancestral sequences. This paper presents a fast, combined solution to these problems. The new algorithm gives an approximation to the minimal history in terms of a distance function on sequences. The distance function on sequences is a minimal weighted path length constructed from substitutions and insertions-deletions of segments of any length. Substitutions are weighted with an arbitrary metric on the set of nucleotides or amino acids, and indels are weighted with a gap penalty function of the form gk = a + (bxk), where k is the length of the indel and a and b are two positive numbers. A novel feature is the introduction of the concept of sequence graphs and a generalization of the traditional dynamic sequence comparison algorithm to the comparison of sequence graphs. Sequence graphs ease several computational problems. They are used to represent large sets of sequences that can then be compared simultaneously. Furthermore, they allow the handling of multiple, equally good, alignments, where previous methods were forced to make arbitrary choices. A program written in C implemented this method; it was tested first on 22 5S RNA sequences.

Algorithms↗

PubMed1994

Intragenic rearrangements in the mitochondrial NADH dehydrogenase subunit 6 gene of vertebrates.

We have sequenced the mitochondrial-encoded NADH dehydrogenase subunit 6 gene from 19 species of birds. Comparison of the derived amino-acid sequences in 22 avian species, six mammals, and two fishes, reveals an intragenic rearrangement in mammals. The C-terminal half of the mammalian protein includes an internal insertion of 10-15 amino acids and a C-terminal deletion of 8-9 amino acids. Based on comparative sequence alignments and hydropathy profile analysis, five hydrophobic segments (designated I to V) corresponding to transmembrane regions are proposed. In this structural model of NADH dehydrogenase subunit 6, the mammalian insertion is found in a variable loop region between transmembrane segments IV and V.

Amino Acid Sequence↗

PubMed2005

Glutathione transferase-like proteins encoded in genomes of yeasts and fungi: insights into evolution of a multifunctional protein superfamily.

Most fungal glutathione transferases (GSTs) do not fit easily into any of the previously characterised classes by immunological, sequence or catalytic criteria. In contrast to the paucity of studies on GSTs cloned or isolated from fungal sources, a screen of databases revealed 67 GST-like sequences from 21 fungal species. Comparison by multiple sequence alignment generated a dendrogram revealing five clusters of GST-like proteins designated clusters 1, 2, EFIBgamma, Ure2p and MAK16, the last three of which have previously been related to the GST superfamily. Surprisingly, a relatively small number of fungal GSTs belong to mainstream classes and the previously-described fungal Gamma class is not widespread in the 21 species studied. Representative crystal structures are available for the EFIBgamma and Ure2p classes and the domain structures of representative sequences are compared with these. In addition, there are some "orphan" sequences that do not fit into any previously-described class, but show similarity to genes implicated in fungal biosynthetic gene clusters. We suggest that GST-like sequences are widespread in fungi, participating in a wide range of functions. They probably evolved by a process similar to domain "shuffling".

Amino Acid Sequence↗

PubMed2005

ConFind: a robust tool for conserved sequence identification.

SUMMARY: ConFind (conserved region finder) identifies regions of conservation in multiple sequence alignments that can serve as diagnostic targets. Designed to work with a large number of closely related, highly variable sequences, ConFind provides robust handling of alignments containing partial sequences and ambiguous characters. Conserved regions are defined in terms of minimum region length, maximum informational entropy (variability) per position, number of exceptions allowed to the maximum entropy criterion and the minimum number of sequences that must contain a non-ambiguous character at a position to be considered for inclusion in a conserved region. Comparison of the calculated entropy for an alignment of 95 influenza A hemagglutinin sequences with random deletions results in a 98% reduction in the average error in ConFind relative to the 'Find Conserved Regions' option in BioEdit. REQUIREMENTS: ConFind requires Python 2.3, but Python 2.4 or an upgrade of the optparse module to Optik 1.5 is suggested. The program is known to run under Linux and DOS.

Algorithms↗

PubMed2001

Enhancement of protein modeling by human intervention in applying the automatic programs 3D-JIGSAW and 3D-PSSM.

Fourteen models were constructed and analyzed for the comparative modeling section of Critical Assessment of Techniques for Protein Structure Prediction (CASP4). Sequence identity between each target and the best possible parent(s) ranged between 55 and 13%, and the root-mean-square deviation between model and target was from 0.8 to 17.9 A. In the fold recognition section, 10 of the 11 remote homologues were recognized. The modeling protocols are a combination of automated computer algorithms, 3D-JIGSAW (for comparative modeling) and 3D-PSSM (for fold recognition), with human intervention at certain critical stages. In particular, intervention is required to check superfamily assignment, best possible parents from which to model, sequence alignments to those parents and take-off regions for modeling variable regions. There now is a convergence of algorithms for comparative modeling and fold recognition, particularly in the region of remote homology.

Algorithms↗

PubMed1997

Mitochondrial DNA polymerases from yeast to man: a new family of polymerases.

We report the sequence of a 4.5-kb cDNA clone isolated from a human melanoma library which bears high amino acid sequence identity to the yeast mitochondrial (mt) DNA polymerase (Mip1p). This cDNA contains a 3720-bp open reading frame encoding a predicted 140-kDa polypeptide that is 43% identical to Mip1p. The N-terminal part of the sequence contains a 13 glutamine stretch encoded by a CAG trinucleotide repeat which is not found in the other DNA polymerases gamma (Pol gamma). Multiple amino acid sequence alignments with Pol gamma from Saccharomyces cerevisiae, Schizosaccharomyces pombe, Pichia pastoris, Drosophila melanogaster, Xenopus laevis and Mus musculus show that these DNA polymerases form a family strongly conserved from yeast to man and are only loosely related to the Family A DNA polymerases.

Amino Acid Sequence↗

PubMed2006

Codep: maximizing co-evolutionary interdependencies to discover interacting proteins.

Approaches for the determination of interacting partners from different protein families (such as ligands and their receptors) have made use of the property that interacting proteins follow similar patterns and relative rates of evolution. Interacting protein partners can then be predicted from the similarity of their phylogenetic trees or evolutionary distances matrices. We present a novel method called Codep, for the determination of interacting protein partners by maximizing co-evolutionary signals. The order of sequences in the multiple sequence alignments from two protein families is determined in such a manner as to maximize the similarity of substitution patterns at amino acid sites in the two alignments and, thus, phylogenetic congruency. This is achieved by maximizing the total number of interdependencies of amino acids sites between the alignments. Once ordered, the corresponding sequences in the two alignments indicate the predicted interacting partners. We demonstrate the efficacy of this approach with computer simulations and in analyses of several protein families. A program implementing our method, Codep, is freely available to academic users from our website: http://www.uhnresearch.ca/labs/tillier/.

Computer Simulation↗

PubMed2002

The unique ligand-binding pocket for the human prostacyclin receptor. Site-directed mutagenesis and molecular modeling.

The human prostacyclin receptor is a seven-transmembrane alpha-helical G-protein coupled receptor, which plays important roles in both vascular smooth muscle relaxation as well as prevention of blood coagulation. The position of the native ligand-binding pocket for prostacyclin as well as other derivatives of the 20-carbon eicosanoid, arachidonic acid, has yet to be determined. Through the use of prostanoid receptor sequence alignments, site-directed mutagenesis, and the 2.8-A x-ray crystallographic structure of bovine rhodopsin, we have developed a three-dimensional model of the agonist-binding pocket within the seven-transmembrane (TM) domains of the human prostacyclin receptor. Upon mutation to alanine, 11 of 29 candidate residues within TM domains II, III, IV, V, and VII exhibited a marked decrease in agonist binding. Of this group, four amino acids, Arg-279 (TMVII), Phe-278 (TMVII), Tyr-75 (TMII), and Phe-95 (TMIII), were identified (via receptor amino acid sequence alignment, ligand structural comparison, and computer-assisted homology modeling) as having direct molecular interactions with ligand side-chain constituents. This binding pocket is distinct from that of the biogenic amine receptors and rhodopsin where the native ligands (also composed of a carbon ring and a carbon chain) are accommodated in an opposing direction. These findings should assist in the development of novel and highly specific ligands including selective antagonists for further molecular pharmacogenetic studies of the human prostacyclin receptor.

Amino Acid Sequence↗

PubMed2001

Identification and characterization of GDP-d-mannose 4,6-dehydratase and GDP-l-fucose snthetase in a GDP-l-fucose biosynthetic gene cluster from Helicobacter pylori.

In this study two open reading frames, namely HP0044 and HP0045 from H. pylori, were cloned and overexpressed in E. coli. The two recombinant proteins were demonstrated to have GDP-d-mannose 4,6-dehydratase (GMD) and GDP-l-fucose synthetase (GFS) activities, respectively. The recombinant GMD was a tetramer and had an optimum pH of 6.5. Exogenous NADP(+) was essential for its activity. The K(m) and K(cat) for GDP-d-mannose were 117.3 microM and 0.27 s(-1), respectively. The recombinant GFS was a homodimer with an optimum pH of 8.0. The K(m) and K(cat) for GDP-4-keto-6-deoxy-d-mannose were 64.08 microM and 0.75 s(-1), respectively. It can use both NADPH and NADH, but less efficient with the latter. Amino acid sequence alignment and phylogenetic analysis showed that H. pylori GFS was highly homologous to the GFS of E. coli O111 and both of them were located on a separate phylogenetic branch from other GFS. The unique clustering and origin of the two genes were also discussed.

Amino Acid Sequence↗

PubMed1998

A new structural model for P-glycoprotein.

Multidrug resistance to anti-cancer drugs is a major medical problem. Resistance is manifested largely by the product of the human MDR1 gene, P-glycoprotein, an ABC transporter that is an integral membrane protein of 1280 amino acids arranged into two homologous halves, each comprising 6 putative transmembrane alpha-helices and an ATP binding domain. Despite the plethora of data from site-directed, scanning and domain replacement mutagenesis, epitope mapping and photoaffinity labeling, a clear structural model for P-glycoprotein remains largely elusive. In this report, we propose a new model for P-glycoprotein that is supported by the vast body of previous data. The model comprises 2 membrane-embedded 16-strand beta-barrels, attached by short loops to two 6-helix bundles beneath each barrel. Each ATP binding domain contributes 2 beta-strands and 1 alpha-helix to the structure. This model, together with an analysis of the amino acid sequence alignment of P-glycoprotein isoforms, is used to delineate drug binding and translocation sites. We show that the locations of these sites are consistent with mutational, kinetic and labeling data.

ATP Binding Cassette Transporter, Subfamily B, Mem↗

PubMed2004

Phage phi29 DNA replication organizer membrane protein p16.7 contains a coiled coil and a dimeric, homeodomain-related, functional domain.

The Bacillus subtilis phage varphi29-encoded membrane protein p16.7 is one of the few proteins known to be involved in prokaryotic membrane-associated DNA replication. Protein p16.7 contains an N-terminal transmembrane domain responsible for membrane localization. A soluble variant lacking the N-terminal membrane anchor, p16.7A, forms dimers in solution, binds to DNA, and has affinity for the varphi29 terminal protein. Here we show that the soluble N-terminal half of p16.7A can form a dimeric coiled coil. However, a second domain, located in the C-terminal half of the protein, has been characterized as being the main domain responsible for p16.7 dimerization. This 70-residue C-terminal domain, named p16.7C, also constitutes the functional part of the protein as it binds to DNA and terminal protein. Sequence alignments, secondary structure predictions, and spectroscopic analyses suggest that p16.7C is evolutionarily related to DNA binding homeodomains, present in many eukaryotic transcriptional regulator proteins. Based on the results, a structural model of p16.7 is presented.

Amino Acid Sequence↗

PubMed2006

Identification of multiple distinct Snf2 subfamilies with conserved structural motifs.

The Snf2 family of helicase-related proteins includes the catalytic subunits of ATP-dependent chromatin remodelling complexes found in all eukaryotes. These act to regulate the structure and dynamic properties of chromatin and so influence a broad range of nuclear processes. We have exploited progress in genome sequencing to assemble a comprehensive catalogue of over 1300 Snf2 family members. Multiple sequence alignment of the helicase-related regions enables 24 distinct subfamilies to be identified, a considerable expansion over earlier surveys. Where information is known, there is a good correlation between biological or biochemical function and these assignments, suggesting Snf2 family motor domains are tuned for specific tasks. Scanning of complete genomes reveals all eukaryotes contain members of multiple subfamilies, whereas they are less common and not ubiquitous in eubacteria or archaea. The large sample of Snf2 proteins enables additional distinguishing conserved sequence blocks within the helicase-like motor to be identified. The establishment of a phylogeny for Snf2 proteins provides an opportunity to make informed assignments of function, and the identification of conserved motifs provides a framework for understanding the mechanisms by which these proteins function.

Adenosine Triphosphatases↗