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

PubMed1996

Succinyl CoA: 3-oxoacid CoA transferase (SCOT): human cDNA cloning, human chromosomal mapping to 5p13, and mutation detection in a SCOT-deficient patient.

Succinyl CoA: 3-oxoacid CoA transferase (SCOT; E.C.2.8.3.5) mediates the rate-determining step of ketolysis in extrahepatic tissues, the esterification of acetoacetate to CoA for use in energy production. Hereditary SCOT deficiency in humans causes episodes of severe ketoacidosis. We obtained human-heart SCOT cDNA clones spanning the entire 1,560-nt coding sequence. Sequence alignment of the human SCOT peptides with other known CoA transferases revealed several conserved regions of potential functional importance. A single approximately 3.2-kb SCOT mRNA is present in human tissues (heart > leukocytes >> fibroblasts), but no signal is detectable in the human hepatoma cell line HepG2. We mapped the human SCOT locus (OXCT) to the cytogenetic band 5p13 by in situ hybridization. From fibroblasts of a patient with hereditary SCOT deficiency, we amplified and cloned cDNA fragments containing the entire SCOT coding sequence. We found a homozygous C-to-G transversion at nt 848, which changes the Ser 283 codon to a stop codon. This mutation (S283X) is incompatible with normal enzyme function and represents the first documentation of a pathogenic mutation in SCOT deficiency.

Amino Acid Sequence↗

PubMed1990

Genomic divergence of an HIV-2 from a German AIDS patient probably infected in Mali.

The complete nucleotide sequence of an HIV-2 isolate derived from a German AIDS patient with predominantly neurological symptoms is reported. The HIV-2BEN sequence is highly divergent from those of previously described HIV-2 and SIV strains. Evolutionary tree analysis of eight HIV-2 sequences reveals the existence of three HIV-2 groups. HIV-2BEN belongs to a group with two isolates from Ghana and The Gambia. Based on a comparison of HIV-2BEN with six HIV-2 isolates, SIVsmm and SIVmac, the variability of the structural env and gag proteins is similar within the HIV-2/SIVsmm/mac and HIV-1 groups. In contrast, the regulatory HIV-1 proteins are more highly conserved than those from HIV-2 strains. Multiple sequence alignments reveal that some domains of the envelope and regulatory proteins are well conserved among HIV-1, HIV-2/SIVsmm/mac, SIVagm and SIVmnd. The identification of conserved domains within the external glycoprotein could help to develop broadly active vaccines.

Acquired Immunodeficiency Syndrome↗

PubMed2003

Isolates of infectious bursal disease virus from India are of very virulent phenotype.

Four Indian field isolates, a classical virulent and an attenuated vaccine strains of Infectious bursal disease virus (IBDV) have been characterized by sequence analysis of part of the VP1 gene (from nucleotide 1538-1979) comprising one of viral RNA dependent RNA polymerase motifs. Sequence alignment of these viruses with reported viruses of other countries revealed Indian IBDV field isolates to be 100% similar to very virulent Japanese (OKYM), European (UK661) and Bangladesh (BD3/99) IBD viruses at amino acid level, whereas they had 0.2-0.9% divergence at nucleotide level. Out of the total 24 nucleotide changes found in the Indian field isolates, as well as reported very virulent viruses, only one resulted in amino acid change S-P at 562 position. The Indian field isolates displayed nucleotide divergence of 10.6-11.6% and amino acid divergence of 2.8-3.5% from the classical virulent and attenuated vaccine strains. The RNA dependent RNA polymerase motif from amino acid 528-541, present in the sequence analyzed, was conserved among all the viruses, irrespective of pathotype and serotype. In the phylogenetic tree, based on nucleotide sequence, Indian field viruses were grouped with reported very virulent viruses in one lineage whereas, classical virulent, attenuated vaccine and serotype 2 strains formed part of the second lineage. But in the phylogenetic tree based on amino acid sequence alignment, the serotype 2 strain OH grouped with Indian field isolates and reported very virulent viruses in one lineage and classical virulent and attenuated vaccine strains formed the second lineage.

Amino Acid Sequence↗

PubMed2003

Pandit: a database of protein and associated nucleotide domains with inferred trees.

MOTIVATION: A large, high-quality database of homologous sequence alignments with good estimates of their corresponding phylogenetic trees will be a valuable resource to those studying phylogenetics. It will allow researchers to compare current and new models of sequence evolution across a large variety of sequences. The large quantity of data may provide inspiration for new models and methodology to study sequence evolution and may allow general statements about the relative effect of different molecular processes on evolution. RESULTS: The Pandit 7.6 database contains 4341 families of sequences derived from the seed alignments of the Pfam database of amino acid alignments of families of homologous protein domains (Bateman et al., 2002). Each family in Pandit includes an alignment of amino acid sequences that matches the corresponding Pfam family seed alignment, an alignment of DNA sequences that contain the coding sequence of the Pfam alignment when they can be recovered (overall, 82.9% of sequences taken from Pfam) and the alignment of amino acid sequences restricted to only those sequences for which a DNA sequence could be recovered. Each of the alignments has an estimate of the phylogenetic tree associated with it. The tree topologies were obtained using the neighbor joining method based on maximum likelihood estimates of the evolutionary distances, with branch lengths then calculated using a standard maximum likelihood approach.

Algorithms↗

PubMed1997

Predicting enzyme function from sequence: a systematic appraisal.

Gapped and ungapped sequence alignment were tested as possible methods to classify proteins into the functional classes defined by the International Enzyme Commission (EC). We exhaustively tested all 15,208 proteins labeled with any EC class in a recent release of the SwissProt database, evaluating all 1,327 relevant EC classes. We effectively tested all possible similarity thresholds that could be used for this assignment through the use of the ROC statistic. Approximately 60% of Enzyme Commission classes containing two or more proteins could not be perfectly discriminated by sequence similarity at any threshold. An analysis of the errors indicates that false positive matches dominate, and that various error mechanisms can be identified, including the multidomain nature of many proteins and polyproteins, convergent evolution, variation in enzyme specificity, and other factors. Many of the putatively false positives are in fact biologically relevant. This work strongly suggests that functional assignment of enzymes should attempt to delimit functionally significant subregions, or domains, before matching to EC classes.

Databases, Factual↗

PubMed1997

Sisyphus and prediction of protein structure.

The problem of predicting protein structure from the sequence remains fundamentally unsolved despite more than three decades of intensive research effort. However, new and promising methods in three-dimensional (3D), 2D and 1D prediction have reopened the field. Mean-force-potentials derived from the protein databases can distinguish between correct and incorrect models (3D). Inter-residue contacts (2D) can be detected by analysis of correlated mutations, albeit with low accuracy. Secondary structure, solvent accessibility and transmembrane helices (1D) can be predicted with significantly improved accuracy using multiple sequence alignments. Some of these new prediction methods have proven accurate and reliable enough to be useful in genome analysis, and in experimental structure determination. Moreover, the new generation of theoretical methods is increasingly influencing experiments in molecular biology.

Computers↗

PubMed2001

Molecular cloning, expression, purification, and characterization of fructose-1,6-bisphosphate aldolase from Thermus aquaticus.

Fructose-1,6-bisphosphate aldolase from the thermophilic eubacteria, Thermus aquaticus YT-1, was cloned and sequenced. Nucleotide-sequence analysis revealed an open reading frame coding for a 33-kDa protein of 305 amino acids having amino acid sequence typical of thermophilic adaptation. Multiple sequence alignment classifies the enzyme as a class II B aldolase that shares similarity with aldolases from other extremophiles: Thermotoga maritima, Aquifex aeolicus, and Helicobacter pylori (49--54% identity, 76--81% homology). Taq FBP aldolase was overexpressed under tac promoter control in Escherichia coli and purified to homogeneity using heat treatment followed by two chromatographic steps. Yields of 40--50 mg of monodisperse protein were obtained per liter of culture. The quaternary structure is that of a homotetramer stabilized by an apparent 21-amino-acid insertion sequence. The recombinant protein is thermostable for at least 45 min at 80 degrees C with little residual activity below 60 degrees C. Kinetic characterization at 70 degrees C, the optimal growth temperature for T. aquaticus, indicates extreme negative subunit cooperativity (h = 0.32) with a limiting K(m) of 305 microM. The maximal specific activity (V(max)) is 46 U/mg at 70 degrees C.

Amino Acid Sequence↗

PubMed2003

Structure-based active site profiles for genome analysis and functional family subclassification.

In previous work, structure-based functional site descriptors, fuzzy functional forms (FFFs), were developed to recognize structurally conserved active sites in proteins. These descriptors identify members of protein families according to active-site structural similarity, rather than overall sequence or structure similarity. FFFs are defined by a minimal number of highly conserved residues and their three-dimensional arrangement. This approach is advantageous for function assignment across broad families, but is limited when applied to detailed subclassification within these families. In the work described here, we developed a method of three-dimensional, or structure-based, active-site profiling that utilizes FFFs to identify residues located in the spatial environment around the active site. Three-dimensional active-site profiling reveals similarities and differences among active sites across protein families. Using this approach, active-site profiles were constructed from known structures for 193 functional families, and these profiles were verified as distinct and characteristic. To achieve this result, a scoring function was developed that discriminates between true functional sites and those that are geometrically most similar, but do not perform the same function. In a large-scale retrospective analysis of human genome sequences, this profile score was shown to identify specific functional families correctly. The method is effective at recognizing the likely subtype of structurally uncharacterized members of the diverse family of protein kinases, categorizing sequences correctly that were misclassified by global sequence alignment methods. Subfamily information provided by this three-dimensional active-site profiling method yields key information for specific and selective inhibitor design for use in the pharmaceutical industry.

Algorithms↗

PubMed2003

A story of chelatase evolution: identification and characterization of a small 13-15-kDa "ancestral" cobaltochelatase (CbiXS) in the archaea.

The cobaltochelatase required for the synthesis of vitamin B12 (cobalamin) in the archaeal kingdom has been identified as CbiX through similarity searching with the CbiX from Bacillus megaterium. However, the CbiX proteins in the archaea are much shorter than the CbiX proteins found in eubacteria, typically containing less than half the number of amino acids in their primary structure. For this reason the shorter CbiX proteins have been termed CbiXS and the longer versions CbiXL. The CbiXS proteins from Methanosarcina barkeri and Methanobacter thermoautotrophicum were overproduced in Escherichia coli as recombinant proteins and characterized. Through complementation studies of a defined chelatase-deficient strain of E. coli and by direct in vitro assays the function of CbiXS as a sirohydrochlorin cobaltochelatase has been demonstrated. On the basis of sequence alignments and conserved active site residues we suggest that CbiXS may represent a primordial chelatase, giving rise to larger chelatases such as CbiXL, SirB, CbiK, and HemH through gene duplication and subsequent variation and selection. A classification scheme for chelatases is proposed.

Amino Acid Sequence↗

PubMed1998

Sequence and structural relationships of leucocins A-, B- and C-TA33a from Leuconostoc mesenteroides TA33a.

Amino acid sequences of two of the three bacteriocins from Leuconostoc mesenteroides TA33a were determined and their sequence-structure relationships investigated. Leucocin B-TA33a consists of 31 amino acid residues, with a molecular mass of 3466 Da. Leucocin B-TA33a does not belong to the pediocin family of bacteriocins, but shares 62% homology with mesenterocin 52B. A partial sequence of 36 amino acids of leucocin C-TA33a (4598 Da) was determined. Leucocin C-TA33a belongs to the class II bacteriocins having the consensus YGNGV motif. The third bacteriocin, leucocin A-TA33a, is identical to leucocin A-UAL 187. Circular dichroism spectra of the leucocins in aqueous solution and micellar SDS indicated that they undergo a structural transition when in a membrane-mimicking environment. Theoretical predictions from circular dichroism data suggest that leucocins A-, B- and C-TA33a adopt a beta-structure (48%) in membrane-mimicking environments. Sequence alignments and secondary structure predictions for the N-terminus of leucocins A- and C-TA33a predicted that Cys-9 and Cys-14 are connected by a disulfide bridge and form two beta-strands.

Amino Acid Sequence↗

PubMed2003

Characterization and transcriptional profiles of three Spodoptera frugiperda genes encoding cysteine-rich peptides. A new class of defensin-like genes from lepidopteran insects?

The present work describes sequence and transcription of three Spodoptera frugiperda genes encoding 6-cysteine-rich peptides. Sequence alignments indicate that the predicted peptides belong to the insect defensin family, although phylogenetic analyses suggest they form a cluster distinct from that of other neopteran insect defensins. The three genes were identified in a non-immune-challenged Sf9 cells cDNA (DNA complementary to RNA) library (Landais et al., Bioinformatics, in press) and were named spodoptericin, Sf-gallerimycin and Sf-cobatoxin. Spodoptericin is a novel defensin-like gene that appears to be weakly up-regulated following injection of bacteria and fungi. Interestingly, no sequence motif clearly homologous to cis regulatory element involved in the regulation of antimicrobial genes was found. An homologue of the spodoptericin gene was identified in the SilkBase Bombyx mori cDNA library. Sf-gallerimycin is related to the Galleria mellonella gallerimycin gene and is induced after immune challenge by injection of bacteria in the larval fat body as well as in hemocytes. In silico analysis of the sequence upstream from the cDNA reveals the presence of at least one motif homologous to a nuclear factor kappaB (NF-kappaB) binding site. Finally, Sf-cobatoxin is related to the G. mellonella cobatoxin-like gene. Despite high levels of constitutive expression compared to the two previous genes, transcription of Sf-cobatoxin is increased after immune, in particular, bacterial challenge. We therefore confirm that these three genes encode potential candidate molecules involved in S. frugiperda innate humoral response.

Amino Acid Sequence↗

PubMed1996

Maize chromosomal HMGc. Two closely related structure-specific DNA-binding proteins specify a second type of plant high mobility group box protein.

The chromosomal high mobility group (HMG) proteins are small and abundant non-histone proteins common to eukaryotes. We have purified the maize HMGc protein from immature kernels and characterized it by mass spectrometry and amino acid sequence analysis. HMGc could be resolved into two similar proteins by reversed phase chromatography. Cloning and characterization of the corresponding cDNAs revealed that they encode two closely related maize HMGc proteins, now termed HMGc1 and HMGc2. Their theoretical masses of 15,316 and 15,007 Da are >300 Da lower than the masses determined for the proteins purified from maize, indicating post-translational modifications of the proteins. Despite sequence similarity to maize HMGa (and previously described homologous proteins of other species) amino acid sequence alignments reveal that HMGc is in several conserved regions distinct from these proteins. Consequently, we have identified a novel type of plant protein containing an HMG box DNA binding domain and belonging to the HMG1 protein family. HMGc1 and HMGc2 were expressed in Escherichia coli, purified to homogeneity, and analyzed for their DNA binding properties. They proved to bind to DNA structure-specifically since they formed complexes with DNA minicircles at concentrations approximately 100-fold lower than the concentrations required to form complexes with linear fragments of identical sequence. Furthermore, HMGc1 and HMGc2 can constrain negative superhelical turns in plasmid DNA.

Amino Acid Sequence↗

PubMed2003

Cloning, sequence and variability analysis of expressed immunoglobulin light chain genes from yellowtail Seriola quinqueradiata.

The cDNA clones encoding immunoglobulin (Ig) light (L) chain variable (V) region associated with constant (C) region were isolated from yellowtail (Seriola quinqueradiata) kidney by expressed sequence tag analysis (accession numbers: AB062619-AB062668, AB064322). The sequences of both VL and CL region contain well-conserved cysteine residues important for intra- and inter-domain interaction in mammals. Comparisons of the amino acid sequence of the CL domain with those of other species showed a high degree of similarity, with 88.3%, 59.8%, and 60.6% to those of wolf fish (Anarhichas lupus), rainbow trout IgL I isotype (Oncorhynchus mykiss) and channel catfish G isotype (Ictalurus punctatus), respectively. Multiple sequence alignments of the CL domain with those of higher vertebrates, however, did not readily allow it to be classified as kappa or lambda isotypes. Furthermore, the pI, hydrophobicity and variability of yellowtail VL regions were studied in 65 cDNA clones and the diversity was observed in CDR1, CDR2 and CDR3 regions.

Amino Acid Sequence↗

PubMed2004

Complete sequence of RNA1 and subgenomic RNA3 of Atlantic halibut nodavirus (AHNV).

The Nodaviridae are divided into the alphanodavirus genus, which infects insects, and the betanodavirus genus, which infects fishes. Betanodaviruses are the causative agent of viral encephalopathy and retinopathy (VER) in a number of cultivated marine fish species. The Nodaviridae are small non-enveloped RNA viruses that contain a genome consisting of 2 single-stranded positivesense RNA segments: RNA1 (3.1 kb), which encodes the viral part of the RNA-dependent RNA polymerase (RdRp); and RNA2 (1.4 kb), which encodes the capsid protein. In addition to RNA1 and RNA2, a subgenomic transcript of RNA1, RNA3, is present in infected cells. We have cloned and sequenced RNA1 from the Atlantic halibut Hippoglossus hippoglossus nodavirus (AHNV), and for the first time, the sequence of a betanodaviral subgenomic RNA3 has been determined. AHNV RNA1 was 3100 nucleotides in length and contained a main open reading frame encoding a polypeptide of 981 amino acids. Conservative motifs for RdRp were found in the deduced amino acid sequence. RNA3 was 371 nucleotides in length, and contained an open reading frame encoding a peptide of 75 amino acids corresponding to a hypothetical B2 protein, although sequence alignments with the alphanodavirus B2 proteins showed only marginal similarities. AHNV RNA replication in the fish cell-line SSN-1 (derived from striped snakehead) was analysed by Northern blot analysis, which indicated that RNA3 was synthesised in large amounts (compared to RNA1) at an early point in time post-infection.

Alternative Splicing↗

PubMed2003

Mapping structural determinants of biological activities in snake venom phospholipases A2 by sequence analysis and site directed mutagenesis.

In addition to their catalytic activity, snake venom phospholipases A2 (vPLA2) present remarkable diversity in their biological effects. Sequence alignment analyses of functionally related PLA2 are frequently used to predict the structural determinants of these effects, and the predictions are subsequently evaluated by site directed mutagenesis experiments and functional assays. In order to improve the predictive potential of computer-based analysis, a simple method for scanning amino acid variation analysis (SAVANA) has been developed and included in the analysis of the lysine 49 PLA2 myotoxins (Lys49-PLA2). The SAVANA analysis identified positions in the C-terminal loop region of the protein, which were not identified using previously available sequence analysis tools. Site directed mutagenesis experiments of bothropstoxin-I, a Lys49-PLA2 isolated from the venom of Bothrops jararacussu, reveals that these residues are exactly those involved in the determination of myotoxic and membrane damaging activities. The SAVANA method has been used to analyse presynaptic neurotoxic and anti-coagulant vPLA2s, and the predicted structural determinants of these activities are in excellent agreement with the available results of site directed mutagenesis experiments. The positions of residues involved in the myotoxic and neurotoxic determinants demonstrate significant overlap, suggesting that the multiple biological effects observed in many snake vPLA2s are a consequence of superposed structural determinants on the protein surface.

Amino Acid Sequence↗

PubMed1998

Isoschizomers of the restriction endonuclease TaqI (T/CGA) requiring different metal ion concentrations and having a range of thermal stabilities from Thermus species from different continents.

One-hundred-and-fifty-two isolates of the genus Thermus, collected from hot springs on four continents, were screened for evidence of the presence of the thermophilic Type II restriction endonuclease TaqI (T/CGA). The presence of isoschizomers of TaqI in 27 of the isolates, originating from hot springs in New Zealand, Iceland, USA, Japan, mainland Portugal and the island of São Miguel in the Azores, is reported. Six of the TaqI-containing isolates from diverse geographical locations, identified by means of DNA/DNA homology and 16S rRNA sequence alignment as belonging to the Thermus species T. aquaticus, T. filiformis, T. thermophilus, T. scotoductus and T. brockianus, were selected for comparative studies. The TaqI isoschizomers from each of the six isolates were partially purified. They differed in their magnesium ion requirements, isoelectric points, subunit molecular masses and thermal stability.

Base Sequence↗

PubMed2004

Identification of the catalytic motif of the microbial ribosome inactivating cytotoxin colicin E3.

Colicin E3 is a cytotoxic ribonuclease that specifically cleaves 16S rRNA at the ribosomal A-site to abolish protein synthesis in sensitive Escherichia coli cells. We have performed extensive mutagenesis of the 96-residue colicin E3 cytotoxic domain (E3 rRNase), assayed mutant colicins for in vivo cytotoxicity, and tested the corresponding E3 rRNase domains for their ability to inactivate ribosome function in vitro. From 21 alanine mutants, we identified five positions where mutation resulted in a colicin with no measurable cytotoxicity (Y52, D55, H58, E62, and Y64) and four positions (R40, R42, E60, and R90) where mutation caused a significant reduction in cytotoxicity. Mutations that were found to have large in vivo and in vitro effects were tested for structural integrity through circular dichroism and fluorescence spectroscopy using purified rRNase domains. Our data indicate that H58 and E62 likely act as the acid-base pair during catalysis with other residues likely involved in transition state stabilization. Both the Y52 and Y64 mutants were found to be highly destabilized and this is the likely origin of the loss of their cytotoxicity. The identification of important active site residues and sequence alignments of known rRNase homologs has allowed us to identify other proteins containing the putative rRNase active site motif. Proteins that contained this active site motif included three hemagglutinin-type adhesins and we speculate that these have evolved to deliver a cytotoxic rRNase into eukaryotic cells during pathogenesis.

Amino Acid Motifs↗

PubMed2001

Molecular cloning and characterization of Mycobacterium bovis BCG pcp gene encoding pyrrolidone carboxyl peptidase.

The Mycobacterium bovis bacilli Calmette-Guerin (BCG) pcp gene that encodes the pyrrolidone carboxyl peptidase (Pcp) was cloned from a lambdagtll genomic library and sequenced. The nucleotide sequence contains a 669 bp open reading frame coding for a protein of 222 amino acid residues with a calculated molecular mass of 23,209 Da. The deduced amino acid sequence is highly homologous to the Pcps from Bacillus amyloliquefaciens, Pseudomonas fluorescens, Bacillus subtilis, Streptococcus pyogenes, and Staphylococcus aureus. A multiple sequence alignment revealed highly conserved domains. The BCG pcp gene was overexpressed in Escherichia coli. The Pcp was purified to homogeneity. The recombinant protein was further confirmed by an enzymatic assay.

Amino Acid Sequence↗