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

PubMed2006

Two-rung model of a left-handed beta-helix for prions explains species barrier and strain variation in transmissible spongiform encephalopathies.

In this study, a new beta-helical model is proposed that explains the species barrier and strain variation in transmissible spongiform encephalopathies. The left-handed beta-helix serves as a structural model that can explain the seeded growth characteristics of beta-sheet structure in PrP(Sc) fibrils. Molecular dynamics simulations demonstrate that the left-handed beta-helix is structurally more stable than the right-handed beta-helix, with a higher beta-sheet content during the simulation and a better distributed network of inter-strand backbone-backbone hydrogen bonds between parallel beta-strands of different rungs. Multiple sequence alignments and homology modelling of prion sequences with different rungs of left-handed beta-helices illustrate that the PrP region with the highest beta-helical propensity (residues 105-143) can fold in just two rungs of a left-handed beta-helix. Even if no other flanking sequence participates in the beta-helix, the two rungs of a beta-helix can give the growing fibril enough elevation to accommodate the rest of the PrP protein in a tight packing at the periphery of a trimeric beta-helix. The folding of beta-helices is driven by backbone-backbone hydrogen bonding and stacking of side-chains in adjacent rungs. The sequence and structure of the last rung at the fibril end with unprotected beta-sheet edges selects the sequence of a complementary rung and dictates the folding of the new rung with optimal backbone hydrogen bonding and side-chain stacking. An important side-chain stack that facilitates the beta-helical folding is between methionine residues 109 and 129, which explains their importance in the species barrier of prions. Because the PrP sequence is not evolutionarily optimised to fold in a beta-helix, and because the beta-helical fold shows very little sequence preference, alternative alignments are possible that result in a different rung able to select for an alternative complementary rung. A different top rung results in a new strain with different growth characteristics. Hence, in the present model, sequence variation and alternative alignments clarify the basis of the species barrier and strain specificity in PrP-based diseases.

Amino Acid Sequence↗

PubMed2004

Accurate detection of very sparse sequence motifs.

Protein sequence alignments are more reliable the shorter the evolutionary distance. Here, we align distantly related proteins using many closely spaced intermediate sequences as stepping stones. Such transitive alignments can be generated between any two proteins in a connected set, whether they are direct or indirect sequence neighbors in the underlying library of pairwise alignments. We have implemented a greedy algorithm, MaxFlow, using a novel consistency score to estimate the relative likelihood of alternative paths of transitive alignment. In contrast to traditional profile models of amino acid preferences, MaxFlow models the probability that two positions are structurally equivalent and retains high information content across large distances in sequence space. Thus, MaxFlow is able to identify sparse and narrow active-site sequence signatures which are embedded in high-entropy sequence segments in the structure based multiple alignment of large diverse enzyme superfamilies. In a challenging benchmark based on the urease superfamily, MaxFlow yields better reliability and double coverage compared to available sequence alignment software. This promises to increase information returns from functional and structural genomics, where reliable sequence alignment is a bottleneck to transferring the functional or structural characterization of model proteins to entire protein superfamilies.

Actins↗

PubMed1995

Sequence homologies between nucleotide binding regions of CFTR and G-proteins suggest structural and functional similarities.

Sequence homology between the alpha-subunits of G-proteins and other GTP-binding proteins and certain regions within the nucleotide binding domains (NBDs) of cystic fibrosis transmembrane conductance regulator (CFTR) indicates that these protein structures may be similar. A sequence alignment of the NBDs of CFTR and NBDs from other membrane transporters, forms the basis of a structural model. This model predicts that one of the conserved sequences GGQR, within which a number of CF mutations occur, forms part of the nucleotide binding pocket and serves as an ON/OFF conformational switch as observed in GTP binding proteins. Furthermore, based on subtle sequence differences between the first and second NBDs of CFTR and from structure-activity data, we suggest that the nucleotide binding site environments of the two NBDs are different.

Amino Acid Sequence↗

PubMed2000

PHAT: a transmembrane-specific substitution matrix. Predicted hydrophobic and transmembrane.

MOTIVATION: Database searching algorithms for proteins use scoring matrices based on average protein properties, and thus are dominated by globular proteins. However, since transmembrane regions of a protein are in a distinctly different environment than globular proteins, one would expect generalized substitution matrices to be inappropriate for transmembrane regions. RESULTS: We present the PHAT (predicted hydrophobic and transmembrane) matrix, which significantly outperforms generalized matrices and a previously published transmembrane matrix in searches with transmembrane queries. We conclude that a better matrix can be constructed by using background frequencies characteristic of the twilight zone, where low-scoring true positives have scores indistinguishable from high-scoring false positives, rather than the amino acid frequencies of the database. The PHAT matrix may help improve the accuracy of sequence alignments and evolutionary trees of membrane proteins.

Algorithms↗

PubMed2005

Better prediction of protein contact number using a support vector regression analysis of amino acid sequence.

BACKGROUND: Protein tertiary structure can be partly characterized via each amino acid's contact number measuring how residues are spatially arranged. The contact number of a residue in a folded protein is a measure of its exposure to the local environment, and is defined as the number of Cbeta atoms in other residues within a sphere around the Cbeta atom of the residue of interest. Contact number is partly conserved between protein folds and thus is useful for protein fold and structure prediction. In turn, each residue's contact number can be partially predicted from primary amino acid sequence, assisting tertiary fold analysis from sequence data. In this study, we provide a more accurate contact number prediction method from protein primary sequence. RESULTS: We predict contact number from protein sequence using a novel support vector regression algorithm. Using protein local sequences with multiple sequence alignments (PSI-BLAST profiles), we demonstrate a correlation coefficient between predicted and observed contact numbers of 0.70, which outperforms previously achieved accuracies. Including additional information about sequence weight and amino acid composition further improves prediction accuracies significantly with the correlation coefficient reaching 0.73. If residues are classified as being either "contacted" or "non-contacted", the prediction accuracies are all greater than 77%, regardless of the choice of classification thresholds. CONCLUSION: The successful application of support vector regression to the prediction of protein contact number reported here, together with previous applications of this approach to the prediction of protein accessible surface area and B-factor profile, suggests that a support vector regression approach may be very useful for determining the structure-function relation between primary protein sequence and higher order consecutive protein structural and functional properties.

Algorithms↗

PubMed1994

SEQSEE: a comprehensive program suite for protein sequence analysis.

SEQSEE (SEQuence SEEker) is a multi-purpose, menu-driven suite of programs designed to provide a fully integrated, state-of-the-art package for the analysis and display of protein sequences and protein databases. It is currently configured to run on most UNIX-based machines including Sun, SGI and NeXT workstations with conversion to other architectures (e.g. Vax or Cray) being a relatively simple task. SEQSEE is capable of performing nearly all of the analytical and comparative tasks found in most comprehensive commercially available software packages. These include sequence/database searching, sequence retrieval, sequence entry and editing, statistical sequence analysis, multiple sequence alignment, flexible pattern matching, and secondary structure prediction. SEQSEE also integrates a number of unique databases which allow it to perform many additional functions such as structure-based sequence alignments and homology-based secondary structure prediction. Additional enhancements to many previously published algorithms have substantially improved the performance of SEQSEE over that found for most other commercial products. The source code, the documentation and all of the required databases for SEQSEE are freely available and may be obtained by anonymous ftp.

Algorithms↗

PubMed1999

Sequence analysis of Pasteurella multocida major outer membrane protein (OmpH) and application of synthetic peptides in vaccination of chickens against homologous strain challenge.

Pasteurella multocida major outer membrane protein (OmpH) has been previously characterized as a porin. The native OmpH from strain X-73 (serotype 1) but not recombinant protein from Escherichia coli induced homologous protection in chickens. In this study OmpH sequences from 15 P. multocida serotypes as well as the CU vaccine strain were compared by sequence alignment and revealed high homology, with major variations confined to two discrete regions which were correspondingly predicted as two largest external loops. Secondary structures of OmpHs were predicted by sequence alignment of OmpHs with well defined porins and analyses of amphiphilicity, hydrophobic moment and antigenic index plots. Several synthetic peptides derived from predicted loop 2 and loop 5 of X-73 OmpH were synthesized as vaccine candidates. Vaccination studies in chickens showed that the cyclic synthetic peptide (Cyclic-L2) mimicking the predicted loop 2 induced 70% protection in chickens against strain X-73 challenge. This is the first report that a synthetic peptide mimicking the conformational epitopes of a native protein provide practical protection in target animal against bacterial infection.

Amino Acid Sequence↗

PubMed1992

Environment-specific amino acid substitution tables: tertiary templates and prediction of protein folds.

The local environment of an amino acid in a folded protein determines the acceptability of mutations at that position. In order to characterize and quantify these structural constraints, we have made a comparative analysis of families of homologous proteins. Residues in each structure are classified according to amino acid type, secondary structure, accessibility of the side chain, and existence of hydrogen bonds from the side chains. Analysis of the pattern of observed substitutions as a function of local environment shows that there are distinct patterns, especially for buried polar residues. The substitution data tables are available on diskette with Protein Science. Given the fold of a protein, one is able to predict sequences compatible with the fold (profiles or templates) and potentially to discriminate between a correctly folded and misfolded protein. Conversely, analysis of residue variation across a family of aligned sequences in terms of substitution profiles can allow prediction of secondary structure or tertiary environment.

Amino Acid Sequence↗

PubMed2003

Subunit E of mitochondrial ATP synthase: a bioinformatic analysis reveals a phosphopeptide binding motif supporting a multifunctional regulatory role and identifies a related human brain protein with the same motif.

The mitochondrial adenosine triphosphate (ATP) synthase is located in the inner membrane and consists of at least 16 subunit types in animals, one of which is subunit e, the function of which is not clearly defined. A highly homologous protein is located in the nucleus and named progesterone receptor binding protein (RBF), to designate its role in this organelle. In addition, the expression level of subunit e in mammalian cells fluctuates greatly and is induced by certain carcinogens and elevated in liver cancers. Because these previous observations suggested to us that subunit e may play multifunctional regulatory roles, we employed a bioinformatic approach to test this view. First, from sequence alignment studies, secondary structure analyses, and basic local alignment search tool (BLAST) searches, we concluded that mitochondrial subunit e and the homologous nuclear protein RBF are most likely the same protein. Second, we examined the known sequence and structure of one of the most common multifunctional cell regulatory proteins, the 14-3-3 protein, involved in phosphopeptide binding, and deduced that it has an apparent binding motif (-KX(6)R---RY-). Third, from careful examination of the conserved residues within all subunit e sequences in the database, we discovered that this protein has a comparable binding motif (-RY---KX(6)R-). Finally, in a BLAST search for additional homologs of subunit e, we found a human brain protein, KIAA1578, the C-terminal 30 amino acids of which are identical to those of human subunit e. This protein also contains a potential phosphopeptide binding motif. In summary, these studies provide support for the view that subunit e is a multifunctional cell regulator involved in cell signaling, and implicate the involvement of the KIAA1578 protein in cell signaling as well. These studies suggest also that, while functioning as a subunit of mitochondrial ATP synthases, subunit e may help regulate these complexes by binding to phosphopeptides within one or more of the other subunit types.

14-3-3 Proteins↗

PubMed2006

Phosphomevalonate kinase: functional investigation of the recombinant human enzyme.

Phosphomevalonate kinase (PMK) catalyzes a key step in isoprenoid/sterol biosynthesis, converting mevalonate 5-phosphate and ATP to mevalonate 5-diphosphate and ADP. To expedite functional and structural study of this enzyme, an expression plasmid encoding His-tagged human PMK has been constructed and recombinant enzyme isolated in an active, stable form. PMK catalyzes a reversible reaction; kinetic constants of human PMK have been determined for both forward (formation of mevalonate 5-diphosphate) and reverse (formation of mevalonate 5-phosphate) reactions. Animal and invertebrate PMKs are not orthologous to plant, fungal, or bacterial PMKs, limiting the information available from sequence alignment analysis. A homology model for the structure of human PMK has been generated. The model conforms to a nucleoside monophosphate kinase family fold. This result, together with sequence comparisons of animal and invertebrate PMKs, suggests an N-terminal basic residue rich sequence as a possible "Walker A" ATP binding motif. The functions of four basic (K17, R18, K19, K22) residues and one acidic (D23) residue in the conserved sequence have been tested by mutagenesis and characterization of isolated mutant proteins. Substrate K(m) values for K17M, R18Q, K19M, and D23N have been measured for forward and reverse reactions; in comparison with wild-type PMK values, only modest (<12-fold) changes are observed. In contrast, R18Q exhibits a V(max) decrease of 100/300-fold (forward/reverse reaction). K22M activity is too low for measurement at nonsaturating substrate concentration; specific activity is decreased by >10000-fold in both forward/reverse reactions, suggesting an active site location and an important role in phosphoryl transfer.

Amino Acid Sequence↗

PubMed2004

An approach to large scale identification of non-obvious structural similarities between proteins.

BACKGROUND: A new sequence independent bioinformatics approach allowing genome-wide search for proteins with similar three dimensional structures has been developed. By utilizing the numerical output of the sequence threading it establishes putative non-obvious structural similarities between proteins. When applied to the testing set of proteins with known three dimensional structures the developed approach was able to recognize structurally similar proteins with high accuracy. RESULTS: The method has been developed to identify pathogenic proteins with low sequence identity and high structural similarity to host analogues. Such protein structure relationships would be hypothesized to arise through convergent evolution or through ancient horizontal gene transfer events, now undetectable using current sequence alignment techniques. The pathogen proteins, which could mimic or interfere with host activities, would represent candidate virulence factors. The developed approach utilizes the numerical outputs from the sequence-structure threading. It identifies the potential structural similarity between a pair of proteins by correlating the threading scores of the corresponding two primary sequences against the library of the standard folds. This approach allowed up to 64% sensitivity and 99.9% specificity in distinguishing protein pairs with high structural similarity. CONCLUSION: Preliminary results obtained by comparison of the genomes of Homo sapiens and several strains of Chlamydia trachomatis have demonstrated the potential usefulness of the method in the identification of bacterial proteins with known or potential roles in virulence.

Bacterial Proteins↗

PubMed2002

GWFASTA: server for FASTA search in eukaryotic and microbial genomes.

Similarity searches are a powerful method for solving important biological problems such as database scanning, evolutionary studies, gene prediction, and protein structure prediction. FASTA is a widely used sequence comparison tool for rapid database scanning. Here we describe the GWFASTA server that was developed to assist the FASTA user in similarity searches against partially and/or completely sequenced genomes. GWFASTA consists of more than 60 microbial genomes, eight eukaryote genomes, and proteomes of annotatedgenomes. Infact, it provides the maximum number of databases for similarity searching from a single platform. GWFASTA allows the submission of more than one sequence as a single query for a FASTA search. It also provides integrated post-processing of FASTA output, including compositional analysis of proteins, multiple sequences alignment, and phylogenetic analysis. Furthermore, it summarizes the search results organism-wise for prokaryotes and chromosome-wise for eukaryotes. Thus, the integration of different tools for sequence analyses makes GWFASTA a powerful toolfor biologists.

Computer Systems↗

PubMed2006

Robust E-values for gapped local alignments.

We examine a Poisson heuristic for judging the significance of local sequence alignments with gaps. Model parameters are estimated directly from the sequences to be aligned, so that laborious prior simulation studies or database comparisons for the estimation of parameters describing the connection between score and E-value are unnecessary. Simulation studies give evidence that this method gives reasonable results even when the usual assumptions like the independence of sequence positions are violated.

Amino Acid Sequence↗

PubMed2000

Structural/functional assignment of unknown bacteriophage T4 proteins by iterative database searches.

Among the total of 274 orfs within bacteriophage T4, only half have been reasonably well characterized, and the functions of the rest have remained obscure. In order to predict the molecular functions of the orfs, a position-specific iterated (PSI)-BLAST search of bacteriophage T4 against the sequence database of known 3D structures was carried out. PSI-BLAST is one of the most powerful iterative sequence search methods using multiple sequence alignment, with the ability to detect many more proteins with distant homology than standard pairwise methods. The 3D structures of proteins are considered to be better preserved than the sequences, and the detected distantly homologous proteins are likely to possess highly similar 3D structures. Thirteen orfs of phage T4, whose homologues were not detected by standard pairwise methods, were found to have significantly homologous counterparts by this method. The plausibility of the results was confirmed by checking whether important residues at substrate/ligand-binding sites were conserved. Among them, two orfs, vs.1 and e.1, which are similar to Escherichia coli lytic enzyme and MutT protein, respectively, had not been studied previously. Also, gp rIIA, a rapid lysis protein, whose gene structure had been intensively studied during the development of molecular biology in the 1950s and yet whose molecular function remains unknown, has an N-terminal domain that is significantly similar to the N-terminal region of the heat shock protein Hsp90.

Amino Acid Sequence↗

PubMed2003

Prokaryotic orthologues of mitochondrial alternative oxidase and plastid terminal oxidase.

The mitochondrial alternative oxidase (AOX) and the plastid terminal oxidase (PTOX) are two similar members of the membrane-bound diiron carboxylate group of proteins. AOX is a ubiquinol oxidase present in all higher plants, as well as some algae, fungi, and protists. It may serve to dampen reactive oxygen species generation by the respiratory electron transport chain. PTOX is a plastoquinol oxidase in plants and some algae. It is required in carotenoid biosynthesis and may represent the elusive oxidase in chlororespiration. Recently, prokaryotic orthologues of both AOX and PTOX proteins have appeared in sequence databases. These include PTOX orthologues present in four different cyanobacteria as well as an AOX orthologue in an alpha-proteobacterium. We used PCR, RT-PCR and northern analyses to confirm the presence and expression of the PTOX gene in Anabaena variabilis PCC 7120. An extensive phylogeny of newly found prokaryotic and eukaryotic AOX and PTOX proteins supports the idea that AOX and PTOX represent two distinct groups of proteins that diverged prior to the endosymbiotic events that gave rise to the eukaryotic organelles. Using multiple sequence alignment, we identified residues conserved in all AOX and PTOX proteins. We also provide a scheme to readily distinguish PTOX from AOX proteins based upon differences in amino acid sequence in motifs around the conserved iron-binding residues. Given the presence of PTOX in cyanobacteria, we suggest that this acronym now stand for plastoquinol terminal oxidase. Our results have implications for the photosynthetic and respiratory metabolism of these prokaryotes, as well as for the origin and evolution of eukaryotic AOX and PTOX proteins.

Amino Acid Sequence↗

PubMed2004

CRASP: a program for analysis of coordinated substitutions in multiple alignments of protein sequences.

Recent results suggest that during evolution certain substitutions at protein sites may occur in a coordinated manner due to interactions between amino acid residues. Information on these coordinated substitutions may be useful for analysis of protein structure and function. CRASP is an Internet-available software tool for the detection and analysis of coordinated substitutions in multiple alignments of protein sequences. The approach is based on estimation of the correlation coefficient between the values of a physicochemical parameter at a pair of positions of sequence alignment. The program enables the user to detect and analyze pairwise relationships between amino acid substitutions at protein sequence positions, estimate the contribution of the coordinated substitutions to the evolutionary invariance or variability in integral protein physicochemical characteristics such as the net charge of protein residues and hydrophobic core volume. The CRASP program is available at http://wwwmgs.bionet.nsc.ru/mgs/programs/crasp/.

Amino Acid Substitution↗

PubMed2003

QAlign: quality-based multiple alignments with dynamic phylogenetic analysis.

Integrating different alignment strategies, a layout editor and tools deriving phylogenetic trees in a 'multiple alignment environment' helps to investigate and enhance results of multiple sequence alignment by hand. QAlign combines algorithms for fast progressive and accurate simultaneous multiple alignment with a versatile editor and a dynamic phylogenetic analysis in a convenient graphical user interface.

Algorithms↗

PubMed2004

Prediction of protein function in the absence of significant sequence similarity.

Tremendous progress in DNA sequencing has yielded the genomes of a host of important organisms. The utilisation of these resources requires understanding of the function of each gene. Standard methods of functional assignment involve sequence alignment to a gene of known function; however such methods often fail to find any significant matches. Here we discuss a number of recent alternative methods that may be of use when sequence alignment fails. Function can be defined in a number of ways including E.C. number and MIPS and KEGG functional classes. Phylogenetic profiles show the pattern of presence or absence of a protein between genomes. Protein-protein interactions can be identified by searching for interacting pairs of proteins that are fused to a single protein chain in another organism. The gene neighbour method uses the observation that if the genes that encode two proteins are close on a chromosome, the proteins tend to be functionally related. More general methods use sequence properties such as amino acid composition, mean hydrophobicity, predicted secondary structure and post-translational modification sites. Data mining methods devise rules in the form of IF... THEN statements that make predictions of function using sequence based attributes, predicted secondary structure and sequence similarity. Finally, structural features can be used, after modelling the structure of a protein from its sequence or solving its structure. Protein fold class can be strongly indicative of function, while other structural features, such as secondary structure content, cleft size and 3D structural motifs are also useful.

Computer Simulation↗