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

PubMed2025

Inference of Cytochrome P450 Evolutionary History Using Structural and Physicochemical Metrics.

Cytochrome P450s are a superfamily of heme-binding monooxygenases involved with the detoxification of intrinsic and extrinsic toxins. They are near ubiquitous within biological domains and are found in all domains. Members of families within the superfamily are defined based on amino acid identity thresholds, with thresholds as low as 40% in some families. Relationships among Cytochrome P450 families have proven elusive due to sub-Twilight Zone interfamily identities (<30%) that result in poor multiple sequence alignment quality and thus low levels of support for downstream phylogenetic reconstructions. Despite the low identities, Cytochrome P450 structures are remarkably well conserved both within and among families. In such cases, structural phylogenetics has the potential to unveil elusive relationships because the selectively favored physicochemical properties giving rise to the structure and function of the proteins persist despite sequence-level divergence. Recently, in two separate publications, we demonstrated that by utilizing physicochemical vectors, dynamic time warping, and hierarchical clustering (PCDTW), large swaths of protein domain families and betacoronavirus receptor-binding domain clades were congruent with validated functional/structural relationships. These were important findings because anomalous sequence alignment-based maximum likelihood phylogenetic findings, which were not congruent with the known functional relationships, were resolved. That also validated the use of physicochemical vectors in making inferences about structural/functional homology. Additionally, it illuminated that the same methods might be applied to other protein families with relationships that are difficult to resolve from sequence data alone. Herein, we used Molecular Weight and Hydrophobicity Physicochemical Dynamic Time Warping (MWHP PCDTW) along with structural and sequence alignment-based phylogenetic methodologies to analyze all of the Cytochrome P450s found both in the high-fidelity Structural Classificaction of Proteins (SCOP) database and the reviewed sequences with both experimentally resolved and de novo predicted structures in the Protein Data Bank and the AlphaFold (AF) Protein Structure Database, respectively. We compared the resulting phylogenetic topologies and found that in some cases, structure-based methods may be less able to resolve random/convergent similarity than physicochemical and sequence-based methodologies. This finding agrees with previous findings that demonstrate the usefulness of physicochemical properties in resolving both random structural similarity and potentially convergent relationships.

Cytochrome P-450 Enzyme System↗

PubMed2004

DPANN: improved sequence to structure alignments following fold recognition.

In fold recognition (FR) a protein sequence of unknown structure is assigned to the closest known three-dimensional (3D) fold. Although FR programs can often identify among all possible folds the one a sequence adopts, they frequently fail to align the sequence to the equivalent residue positions in that fold. Such failures frustrate the next step in structure prediction, protein model building. Hence it is desirable to improve the quality of the alignments between the sequence and the identified structure. We have used artificial neural networks (ANN) to derive a substitution matrix to create alignments between a protein sequence and a protein structure through dynamic programming (DPANN: Dynamic Programming meets Artificial Neural Networks). The matrix is based on the amino acid type and the secondary structure state of each residue. In a database of protein pairs that have the same fold but lack sequences-similarity, DPANN aligns over 30% of all sequences to the paired structure, resembling closely the structural superposition of the pair. In over half of these cases the DPANN alignment is close to the structural superposition, although the initial alignment from the step of fold recognition is not close. Conversely, the alignment created during fold recognition outperforms DPANN in only 10% of all cases. Thus application of DPANN after fold recognition leads to substantial improvements in alignment accuracy, which in turn provides more useful templates for the modeling of protein structures. In the artificial case of using actual instead of predicted secondary structures for the probe protein, over 50% of the alignments are successful.

Amino Acid Sequence↗

PubMed2004

Protein structure alignment by deterministic annealing.

MOTIVATION: Protein structure alignment is one of the most important computational problems in molecular biology and plays a key role in protein structure prediction, fold family classification, motif finding, phylogenetic tree reconstruction and so on. From the viewpoint of computational complexity, a pairwise structure alignment is also a NP-hard problem, in contrast to the polynomial time algorithm for a pairwise sequence alignment. RESULTS: We propose a method for solving the structure alignment problem in an accurate manner at the amino acid level, based on a mean field annealing technique. We define the structure alignment as a mixed integer-programming (MIP) problem. By avoiding complicated combinatorial computation and exploiting the special structure of the continuous partial problem, we transform the MIP into a reduced non-linear continuous optimization problem (NCOP) with a much simpler form. To optimize the reduced NCOP, a mean field annealing procedure is adopted with a modified Potts model, whose solution is generally identical to that of the MIP. There is no 'soft constraint' in our mean field model and all constraints are automatically satisfied throughout the annealing process, thereby not only making the optimization more efficient but also eliminating many unnecessary parameters that depend on problems and usually require careful tuning. A number of benchmark examples are tested by the proposed method with comparisons to several existing approaches.

Algorithms↗

PubMed1997

Aligning a DNA sequence with a protein sequence.

We develop several algorithms for the problem of aligning DNA sequence with a protein sequence. Our methods account for frameshift errors, but not for introns in the DNA sequence. Thus, they are particularly appropriate for comparing a cDNA sequence that suffers from sequencing errors with an amino acid sequence or a protein sequence database. We describe algorithms for computing optimal alignments for several definitions of DNA-protein alignment, verify sufficient conditions for equivalence of certain definitions, describe techniques for efficient implementation, and discuss experience with these ideas in a new release of the FASTA suite of database-searching programs.

Algorithms↗

PubMed1998

Self-organizing tree-growing network for the classification of protein sequences.

The self-organizing tree algorithm (SOTA) was recently introduced to construct phylogenetic trees from biological sequences, based on the principles of Kohonen's self-organizing maps and on Fritzke's growing cell structures. SOTA is designed in such a way that the generation of new nodes can be stopped when the sequences assigned to a node are already above a certain similarity threshold. In this way a phylogenetic tree resolved at a high taxonomic level can be obtained. This capability is especially useful to classify sets of diversified sequences. SOTA was originally designed to analyze pre-aligned sequences. It is now adapted to be able to analyze patterns associated to the frequency of residues along a sequence, such as protein dipeptide composition and other n-gram compositions. In this work we show that the algorithm applied to these data is able to not only successfully construct phylogenetic trees of protein families, such as cytochrome c, triosephophate isomerase, and hemoglobin alpha chains, but also classify very diversified sequence data sets, such as a mixture of interleukins and their receptors.

Algorithms↗

PubMed2002

Cloning and expression of a lombricine kinase from an echiuroid worm: insights into structural correlates of substrate specificity.

Phosphagen kinases constitute a large family of enzymes catalyzing the reversible phosphorylation of guanidino acceptor compounds. These guanidino substrates differ substantially in size and chemical properties. In spite of the appearance of X-ray crystal structures for two members of this family, creatine kinase (CK) and arginine kinase (AK), the structural correlates of substrate specificity remain to be fully elucidated. We have determined the cDNA and deduced amino acid sequences for lombricine (guanidinethylphosphoserine) kinase (LK) from the echiuroid worm Urechis caupo and expressed the cDNA in Escherichia coli. The recombinant protein was purified by affinity chromatography and showed high capacity for phosphorylation of lombricine. Phosphagen kinases consist of a small, N-terminal domain and a much larger domain connected by a linker sequence. A key event in catalysis in CK and AK, and certainly all other phosphagen kinases, is a large conformational change involving involving a rotation of the two domains and the movement of two highly conserved flexible loops (one located in the small domain; the other located in the large domain of these enzymes) which clamp down on the substrates. Multiple sequence alignments of Urechis LK with the only other LK sequence available and CK, AK and glycocyamine kinase sequences, confirm the importance of the small flexible loop located in the N-terminal domain of phosphagen kinases as one component of the structural determinants of guanidine specificity. The role of the other flexible loop in the large domain in terms of substrate specificity remains questionable.

Amino Acid Sequence↗

PubMed1994

Analysis of the single-stranded DNA bacteriophage phi X174, refined at a resolution of 3.0 A.

The structure of the bacteriophage phi X174 was examined in a 2.7 A resolution map and refined, using 6.0 A to 3.0 A resolution data with F > or = 5 sigma (F). The final R-factor was 20.9% and the root-mean-square deviation from idealized bond lengths was 0.021 A. The Hendrickson-Konnert refinement was restrained by the phases derived from the molecular replacement icosahedral averaging procedure. The mature phage capsid consists of 60 copies of the F protein with 426 amino acids, the G protein with 175 amino acids and the J protein with 37 amino acids, as well as 12 copies of the H protein with 328 amino acids. The entire polypeptide chain of the F and G protein, all but the first N-terminal residue of the J protein, and 178 solvent molecules were included in the refinement calculations. The secondary structural features of the F, G and J proteins and their interactions with each other are described. The majority of the protein-protein interactions are between the icosahedral 5-fold related interfaces of the F and of the G proteins. These pentameric units of the F and G proteins form the 9S and 6S assembly intermediates, respectively. The J protein lacks any secondary structure and acts as a linking arm between the icosahedral 5-fold related F proteins. Water molecules were introduced only after phase extension to 2.7 A resolution had been completed. The F protein is associated with lower "thermal" parameters and exhibits greater water order in its environment than the G and J proteins. The largest thermal parameters occur in residues on the viral surface. The solvent contributes to the interactions between the proteins. There is an interface of solvent molecules between the F and the G pentamers which stabilizes the pentameric G protein spikes in a crater centered at each of the icosahedral 5-fold vertices of the F protein capsid. Sequence alignments of the F, G and J amino acid sequences for the homologous bacteriophages G4, alpha 3, phi K and phi X174 with respect to the phi X174 structure demonstrated the conservation of functionally important residues on the viral surface.

Amino Acid Sequence↗

PubMed2001

Sequence and structure of the mouse connexin45 gene.

The connexin45 (Cx45) gene was cloned from a mouse genomic Bacterial Artificial Chromosome library. Approximately 8.4 kb of the genomic DNA was sequenced, and the structure of the Cx45 gene was determined. The mouse Cx45 gene is composed of 3 exons, with the entire coding sequence contained within exon III (EMBL Accession Number AJ300716). This structure is unique for the Cx45 gene, since all other members of the connexin family have only two exons. In addition, computer analysis reveals a potential TATA box and two putative AP-1 binding sites in the 5' region of the gene. Sequence alignment with connexin43 indicates substantial homology in the intronic sequences upstream of the 3' exons of the two genes, suggesting that the Cx45 gene is inherently similar to the rest of the connexin family, and that it probably evolved from an ancestor common to the other connexins.

Animals↗

PubMed1997

A prediction of the amino acids and structures involved in DNA recognition by type I DNA restriction and modification enzymes.

The S subunits of type I DNA restriction/modification enzymes are responsible for recognising the DNA target sequence for the enzyme. They contain two domains of approximately 150 amino acids, each of which is responsible for recognising one half of the bipartite asymmetric target. In the absence of any known tertiary structure for type I enzymes or recognisable DNA recognition motifs in the highly variable amino acid sequences of the S subunits, it has previously not been possible to predict which amino acids are responsible for sequence recognition. Using a combination of sequence alignment and secondary structure prediction methods to analyse the sequences of S subunits, we predict that all of the 51 known target recognition domains (TRDs) have the same tertiary structure. Furthermore, this structure is similar to the structure of the TRD of the C5-cytosine methyltransferase, Hha I, which recognises its DNA target via interactions with two short polypeptide loops and a beta strand. Our results predict the location of these sequence recognition structures within the TRDs of all type I S subunits.

Amino Acid Sequence↗

PubMed2004

The PredictProtein server.

PredictProtein (http://www.predictprotein.org) is an Internet service for sequence analysis and the prediction of protein structure and function. Users submit protein sequences or alignments; PredictProtein returns multiple sequence alignments, PROSITE sequence motifs, low-complexity regions (SEG), nuclear localization signals, regions lacking regular structure (NORS) and predictions of secondary structure, solvent accessibility, globular regions, transmembrane helices, coiled-coil regions, structural switch regions, disulfide-bonds, sub-cellular localization and functional annotations. Upon request fold recognition by prediction-based threading, CHOP domain assignments, predictions of transmembrane strands and inter-residue contacts are also available. For all services, users can submit their query either by electronic mail or interactively via the World Wide Web.

Internet↗

PubMed2002

Clustering patterns of cytotoxic T-lymphocyte epitopes in human immunodeficiency virus type 1 (HIV-1) proteins reveal imprints of immune evasion on HIV-1 global variation.

The human cytotoxic T-lymphocyte (CTL) response to human immunodeficiency virus type 1 (HIV-1) has been intensely studied, and hundreds of CTL epitopes have been experimentally defined, published, and compiled in the HIV Molecular Immunology Database. Maps of CTL epitopes on HIV-1 protein sequences reveal that defined epitopes tend to cluster. Here we integrate the global sequence and immunology databases to systematically explore the relationship between HIV-1 amino acid sequences and CTL epitope distributions. CTL responses to five HIV-1 proteins, Gag p17, Gag p24, reverse transcriptase (RT), Env, and Nef, have been particularly well characterized in the literature to date. Through comparing CTL epitope distributions in these five proteins to global protein sequence alignments, we identified distinct characteristics of HIV amino acid sequences that correlate with CTL epitope localization. First, experimentally defined HIV CTL epitopes are concentrated in relatively conserved regions. Second, the highly variable regions that lack epitopes bear cumulative evidence of past immune escape that may make them relatively refractive to CTLs: a paucity of predicted proteasome processing sites and an enrichment for amino acids that do not serve as C-terminal anchor residues. Finally, CTL epitopes are more highly concentrated in alpha-helical regions of proteins. Based on amino acid sequence characteristics, in a blinded fashion, we predicted regions in HIV regulatory and accessory proteins that would be likely to contain CTL epitopes; these predictions were then validated by comparison to new sets of experimentally defined epitopes in HIV-1 Rev, Tat, Vif, and Vpr.

Amino Acid Sequence↗

PubMed2001

Nuclear factor 1 is a negative regulator of gadd153 gene expression in vascular smooth muscle cells.

Growth arrest and DNA damage inducible gene 153 (gadd153) is expressed at very low levels in growing cells but is markedly induced in response to cellular stresses, including glucose deprivation, exposure to genotoxic agents, and other growth-arresting situations. Forced expression of GADD153 can induce cell cycle arrest and/or apoptosis in many types of cells. Recently, we reported that GADD153 was induced in vascular smooth muscle cells (VSMCs) in neointimal lesions of balloon-injured carotid arteries. To investigate the underlying molecular mechanisms of gadd153 gene expression in VSMCs, we isolated and characterized a promoter region of the rat gadd153 gene. Sequence alignments of this region revealed 1 TATA-like sequence and several well-known cis elements. The 5'-deletion analysis for this region showed that a domain spanning -447 through -368 drastically reduced the promoter activity to almost equal levels of promoterless control. Because this domain contained a consensus sequence for the nuclear factor 1 family of proteins (NF1), DNA-binding studies were performed by use of 2 types of NF1 consensus probes. Both probes were specifically shifted by nuclear extracts from proliferating VSMCs and were supershifted by antiserum against CCAAT transcription factor/NF1. In addition, promoter activity of a mutant luciferase vector, which was generated by a point mutation at the NF1 binding motif of the gadd153 gene, was 14-fold higher than that of a wild-type one. These results suggest that gadd153 gene expression in VSMCs is negatively regulated by an NF1-binding motif, and NF1 may act as an antiapoptotic factor by continuously suppressing gadd153 gene expression in growing VSMCs.

Animals↗

PubMed2006

Multiple mapping method: a novel approach to the sequence-to-structure alignment problem in comparative protein structure modeling.

A major bottleneck in comparative protein structure modeling is the quality of input alignment between the target sequence and the template structure. A number of alignment methods are available, but none of these techniques produce consistently good solutions for all cases. Alignments produced by alternative methods may be superior in certain segments but inferior in others when compared to each other; therefore, an accurate solution often requires an optimal combination of them. To address this problem, we have developed a new approach, Multiple Mapping Method (MMM). The algorithm first identifies the alternatively aligned regions from a set of input alignments. These alternatively aligned segments are scored using a composite scoring function, which determines their fitness within the structural environment of the template. The best scoring regions from a set of alternative segments are combined with the core part of the alignments to produce the final MMM alignment. The algorithm was tested on a dataset of 1400 protein pairs using 11 combinations of two to four alignment methods. In all cases MMM showed statistically significant improvement by reducing alignment errors in the range of 3 to 17%. MMM also compared favorably over two alignment meta-servers. The algorithm is computationally efficient; therefore, it is a suitable tool for genome scale modeling studies.

Algorithms↗

PubMed1990

Conservation of functional residues between yeast and E. coli inorganic pyrophosphatases.

The alignments of the amino acid sequences of inorganic pyrophosphatase (PPase) from Saccharomyces cerevisiae (Y1-PPase, 286 amino acids) and Escherichia coli (E-PPase, 175 amino acids) are examined in the light of crystallographic and chemical modification results placing specific amino acid residues at the active site of the yeast enzyme. The major results are: (1) the full E-PPase sequence aligns within residues 28-225 of Y1-PPase, raising the possibility that the N-terminal and C-terminal portions of Y1-PPase may not be essential for activity, and (2) that whereas the overall identity between the two sequences is only modest (22-27% depending on the choice of alignment parameters), of some 17 putative active site residues, 14-16 are identical between Y-PPase and E-PPase. PPase thus appears to be an example of enzymes from widely divergent species that conserve common functional elements within the context of substantial overall sequence variation.

Amino Acid Sequence↗

PubMed2006

An amino acid substitution matrix for protein conformation identification.

Amino acid substitution matrices play an essential role in protein sequence alignment, a fundamental task in bioinformatics. Most widely used matrices, such as PAM matrices derived from homologous sequences and BLOSUM matrices derived from aligned segments of PROSITE, did not integrate conformation information in their construction. There are a few structure-based matrices, which are derived from limited data of structure alignment. Using databases PDB_SELECT and DSSP, we create a database of sequence-conformation blocks which explicitly represent sequence-structure relationship. Members in a block are identical in conformation and are highly similar in sequence. From this block database, we derive a conformation-specific amino acid substitution matrix CBSM60. The matrix shows an improved performance in conformational segment search and homolog detection.

Amino Acid Sequence↗

PubMed2025

Columba: fast approximate pattern matching with optimized search schemes.

MOTIVATION: Aligning sequencing reads to reference genomes is a fundamental task in bioinformatics. Aligners can be classified as lossy or lossless: lossy aligners prioritize speed by reporting only one or a few high-scoring alignments, whereas lossless aligners output all optimal alignments, ensuring completeness and sensitivity. RESULTS: This paper introduces Columba, a high-performance lossless aligner tailored for Illumina sequencing data. Columba processes single or paired-end reads in FASTQ format and outputs alignments in SAM format. By utilizing advanced search schemes and bit-parallel alignment techniques, Columba achieves exceptional speed. Columba is available in two variants. The first, based on the bidirectional FM-index, prioritizes speed. The second, Columba RLC, uses run-length compression using a bidirectional move structure, significantly reducing memory usage for large, repetitive datasets like pan-genomes. Benchmarks on the human genome, as well as bacterial and human pan-genome datasets, demonstrate that Columba is much faster than existing lossless aligners and even competitive with lossy tools. We integrated Columba into the OptiType HLA genotyping pipeline, where it substantially reduced computational time while maintaining accuracy. These results position Columba as a versatile, state-of-the-art tool for high-sensitivity genomic analyses. AVAILABILITY AND IMPLEMENTATION: The source code of Columba is available at https://github.com/biointec/columba under AGPL license. Scripts to reproduce the benchmarks and analyses are available at https://doi.org/10.5281/zenodo.15849246.

Software↗

PubMed2004

Structural exon database, SEDB, mapping exon boundaries on multiple protein structures.

UNLABELLED: Comparative analysis of exon/intron organization of genes and their resulting protein structures is important for understanding evolutionary relationships between species, rules of protein organization and protein functionality. We present Structural Exon Database (SEDB), with a Web interface, an application that allows users to retrieve the exon/intron organization of genes and map the location of the exon boundaries and the intron phase onto a multiple structural alignment. SEDB is linked with Friend, an integrated analytical multiple sequence/structure viewer, which allows simultaneous visualization of exon boundaries on structure and sequence alignments. With SEDB researchers can study the correlations of gene structure with the properties of the encoded three-dimensional protein structures across eukaryotic organisms. AVAILABILITY: SEDB is publicly available at http://glinka.bio.neu.edu/SEDB/SEDB.html SUPPLEMENTARY INFORMATION: On the SEDB Web site.

Databases, Genetic↗

PubMed1998

Crystal structure of the lysozyme from bacteriophage lambda and its relationship with V and C-type lysozymes.

Like other lysozymes, the bacteriophage lambda lysozyme is involved in the digestion of bacterial walls. This enzyme is remarkable in that its mechanism of action is different from the classical lysozyme's mechanism. From the point of view of protein evolution, it shows features of lysozymes from different classes. The crystal structure of the enzyme in which all tryptophan residues have been replaced by aza-tryptophan has been solved by X-ray crystallography at 2.3 A using a combination of multiple isomorphous replacement, non-crystallographic symmetry averaging and density modification techniques. There are three molecules in the asymmetric unit. The characteristic structural elements of lysozymes are conserved: each molecule is organized in two domains connected by a helix and the essential catalytic residue (Glu19) is located in the depth of a cleft between the two domains. This cleft shows an open conformation in two of the independent molecules, while access to the cavity is much more restricted in the last one. A structural alignment with T4 lysozyme and hen egg white lysozyme allows us to superpose about 60 C alpha atoms with a rms distance close to 2 A. The best alignments concern the helix preceding the catalytic residue, some parts of the beta sheets and the helix joining the two domains. The results of sequence alignments with the V and C lysozymes, in which weak local similarities had been detected, are compared with the structural results.

Amino Acid Sequence↗

Refine your search to explore more results.