PubMed Health⌕ Search

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 991 records · Page 55Linked to original sources

Prediction of the structural similarity between spermidine/putrescine-binding protein and maltose-binding protein.

A structural similarity between spermidine/putrescine-binding protein and maltose-binding protein has been predicted by a sequence-structure compatibility method. The sequence alignment obtained by this method revealed a consensus sequence motif located on the surface loop between the first alpha-helix and the second beta-strand, and the further analysis identified a similar motif in iron-binding protein. The conservation of this motif among certain bacterial periplasmic binding proteins suggests a common functional role for this region as well as an evolutionary relationship between them.

ATP-Binding Cassette Transporters↗

Structure and hydride transfer mechanism of a moderate thermophilic dihydrofolate reductase from Bacillus stearothermophilus and comparison to its mesophilic and hyperthermophilic homologues.

Dihydrofolate reductase (DHFR) from a moderate thermophilic organism, Bacillus stearothermophilus, has been cloned and expressed. Physical characterization of the protein (BsDHFR) indicates that it is a monomeric protein with a molecular mass of 18,694.6 Da (0.8), coincident with the mass of 18 694.67 Da calculated from the primary sequence. Determination of the X-ray structure of BsDHFR provides the first structure for a monomeric DHFR from a thermophilic organism, indicating a high degree of conservation of structure in relation to all chromosomal DHFRs. Structurally based sequence alignment of DHFRs indicates the following levels of sequence identity and similarity for BsDHFR: 38 and 58% with Escherichia coli, 35 and 56% with Lactobacillus casei, and 23 and 40% with Thermotoga maritima, respectively. Steady state kinetic isotope effect studies indicate an ordered kinetic mechanism at elevated temperatures, with NADPH binding first to the enzyme. This converts to a more random mechanism at reduced temperatures, reflected in a greatly reduced K(m) for dihydrofolate at 20 degrees C in relation to that at 60 degrees C. A reduction in either temperature or pH reduces the degree to which the hydride transfer step is rate-determining for the second-order reaction of DHF with the enzyme-NADPH binary complex. Transient state kinetics have been used to study the temperature dependence of the isotope effect on hydride transfer at pH 9 between 10 and 50 degrees C. The data support rate-limiting hydride transfer with a moderate enthalpy of activation (E(a) = 5.5 kcal/mol) and a somewhat greater temperature dependence for the kinetic isotope effect than predicted from classical behavior [A(H)/A(D) = 0.57 (0.15)]. Comparison of kinetic parameters for BsDHFR to published data for DHFR from E. coli and T. maritima shows a decreasing trend in efficiency of hydride transfer with increasing thermophilicity of the protein. These results are discussed in the context of the capacity of each enzyme to optimize H-tunneling from donor (NADPH) to acceptor (DHF) substrates.

Amino Acid Sequence↗

Identification of evolutionary conserved regulatory sequences in the 5' untranscribed region of the neural-specific ubiquitin C-terminal hydrolase (PGP9.5) gene.

The structure at the 5' end of the gene encoding neural-specific protein gene product 9.5 (PGP9.5) has been compared between two evolutionary distant species: the human and Monodelphis domestica. In contrast to the highly conserved coding sequences of the gene, only a 48% identity was found across a 1-kb stretch of 5' untranslated and untranscribed DNA. Promoter function studies performed on the human sequence identified a 233-bp CpG-rich minimal promoter. Truncation mutagenesis revealed the presence of essential positive cis-acting regulatory sequences within the region -182 to -123 relative to the transcription initiation site. Sequence alignment analysis of the human and Monodelphis promoter sequences showed 76% identity in this 59-bp region of the gene. A perfectly conserved 12-bp sequence (PSN) located within this region acts as a non-cell-specific activator of transcription in a heterologous reporter gene (pBLCAT2). PGP9.5 gene expression can be readily detected in human neuroblastoma cell lines but is absent in nonneuronal cell lines such as HeLa. Studies on the cell type specificity of the human PGP9.5 promoter demonstrated that in contrast to the endogenous gene, the promoter is active in HeLa cells. However, the promoter displays higher levels of activity in human neuroblastoma cell lines. A conserved 16-bp sequence located at -356 (motif 5) was able to reduce the activity of a heterologous minimal promoter specifically in HeLa cells. In conclusion, we have shown that expression of the PGP9.5 gene is regulated by evolutionary conserved positive and negative cis-acting sequences located in the untranscribed region of the gene.

Animals↗

Rapidly evolving aphid gall effector proteins exhibit saposin-like folds.

Many insects manipulate plants by injecting effector proteins. In one extreme example of this molecular "hijacking," Hormaphis cornu aphids inject bicycle proteins into Hamamelis virginiana, contributing to the development of novel organs called galls. Bicycle proteins share no amino acid sequence similarity with proteins of known function. Here, we report the crystal structures of two divergent bicycle proteins. Both proteins contain saposin-like folds: one with multiple disulfide bonds exhibits a swapped domain topology; the other has no disulfide bonds and possesses two distinct, tandem domains. To explore the structural evolution of bicycle proteins, we attempted to predict bicycle protein structures with Alphafold2 (AF2) and other deep learning programs. While AF2 did not recover the two experimental structures using existing databases, it succeeded when provided with multiple sequence alignments (MSAs) of protein sequences from newly sequenced closely related species. Using this approach, we generated 2,400 high-confidence bicycle protein predictions from seven aphid species. While all aphid bicycle proteins contain predicted saposin-like folds, they display a vast diversity of structural and physicochemical properties. While this diversity thwarts prediction of conserved functions encoded in structure, it suggests that bicycle proteins have evolved to target diverse plant processes and/or to evade plant immune surveillance. Our extension of AF2 with custom MSAs of proteins from closely related species provides a generalizable, powerful approach for predicting structures of rapidly evolving protein families.

Animals↗

Molecular cloning of nuclear receptor coactivator-1 gene in pig.

Nuclear receptor coactivator-1 (NCOA-1) as a member of Steroid receptor coactivator (SRC) family can interact with some DNA-bound nuclear receptors and enhance their transcriptional activation function, which is physiologically and pathologically important for animals and human. Much work has been done on the gene in human and mouse, but few in pig. Here we cloned and sequenced full length cDNA of NCOA-1 gene in pig. The putative protein includes 1440 amino acids. Sequences alignment showed that pig NCOA-1 gene cDNA sequence has identity of 93.73% with human and mouse, and putative pig NCOA-1 protein has identity of 94.99% with human, mouse and rat. The motifs, LXXLL, which have been demonstrated to be important for NCOA-1 mediating ligand-dependent and direct interaction with nuclear receptor, were conserved in pig.

Amino Acid Sequence↗

Homology modeling of the central catalytic domain of insertion sequence ISLC3 isolated from Lactobacillus casei ATCC 393.

The tertiary structure of the central catalytic domain of insertion sequence ISLC3 isolated from Lactobacillus casei ATCC 393 was predicted using the homology modeling approach. The novel insertion sequence was isolated by us from the template bacteriophage phiA3 of L.casei ATCC 393. The number of amino acid residues of the ISLC3 central catalytic domain was 116 and was treated as the query sequence. There were five Web-available threading methods used to find some primary structure templates for the query sequence. These primary templates were further screened using the SWISS-MODEL Protein Modeling Server and the default parameter settings therein to give six final structure templates. All of these final structure templates were the integrase (IN) protein of retroviruses. Multiple sequence alignment using these IN sequences against the query one revealed the signature DDE motif. Based on the structures of these final templates, the structure of the query sequence was constructed using the InsightII/Discover/Homology programs. A metal ion, Mg(2+), was inserted into the center of the putative catalytic pocket formed by the DDE residues of the predicted structure in the final rounds of refinement by molecular dynamics (MD) simulations. The structure with a metal ion included was designated with Mg and that without a metal ion was designated free Mg. The average exposed surface area of some hydrophobic residues of both the predicted free Mg and with Mg structures were computed and compared with those computed for the six structure templates. Whereas the predicted with Mg structure was slightly more exposed than the predicted free Mg structure, the former appeared to be more stable than the latter, as revealed by the lower conformation energy recorded for the former during the structure refinement by MD simulations. To verify further the predicted structures, the coordinates of both predicted structures were fed into the ERRAT Protein Verification Server. It was found that the quality of the predicted with Mg structure was much better than that of the free Mg structure. The validation results also indicated that regions of the predicted with Mg structure that can be rejected at the 95% confidence level were approximately 20% whereas those which can be rejected at the same level for the six structure templates were approximately 10%. The predicted with Mg structure was also docked into a short oligonucleotide representing the substrate of the ISLC3 transposase using the DOCK_4.0.2 program. It was found that both Glu140 and Asp68 residues of the DDE motif of the predicted with Mg structure were able to form hydrogen bonds with the DNA substrate, which was similar to what was observed in a docking study using the retrovirus IN 1asu and its DNA substrate.

Amino Acid Sequence↗

Ty1-copia group retrotransposons in persimmon (Diospyros kaki Thunb.).

We cloned and characterized Ty1-copia group retrotransposons in persimmon (Diospyros kaki Thunb.). Genomic DNA or methyl jasmonate (MJA)-treated cDNA were used as templates to amplify the reverse transcriptase region of Ty1-copia group retrotransposons. About 280 bp fragments were amplified and cloned, and 97 clones were sequenced. Forty-nine clones included frameshift or the stop codon, or both. Among 48 clones containing complete reading frames, 42 clones had unique nucleotide sequences. Alignment and phylogenetic analysis of putative amino acid sequences in the 42 clones indicated that these clones (named Tdk; retroTransposon in Diospyros kaki) fell into seven subgroups and six ungrouped sequences, indicating high sequence heterogeneity in Tdk clones. Phylogenetic analysis comparing unrelated plant species shows that some Tdk clones are more closely related to Ty1-copia group retrotransposons in the orders Solanales and Sapindales rather than to other Tdk clones. Southern blot analysis using Tdk2B, Tdk4c, Tdk6Ac, Tdk12K and Tdk13G clones as probes showed that persimmon and its related species, D. lotus, D. lotus var. glabba, D. oleifera, D. rhombifolia and D. virginiana, contained multiple Tdk-like sequences, indicating that homologous elements exist in other Diospyros species.

Acetates↗

ProtoMatch: a tool for analyzing high-density, sequential eye gaze and cursor protocols.

ProtoMatch is a software tool for integrating and analyzing fixed-location and movement eye gaze and cursor data. It provides a comprehensive collection of protocol analysis tools that support sequential data analyses for eye fixations and scanpaths as well as for cursor "fixations" (dwells at one location) and "cursorpaths" (movements between locations). ProtoMatch is modularized software that integrates both eye gaze and cursor protocols into a unified stream of data and provides an assortment of filters and analyses. ProtoMatch subsumes basic analyses (i.e., fixation duration, number of fixations, etc.) and introduces a method of objectively computing the similarity between scanpaths or cursorpaths using sequence alignment. The combination of filters, basic analyses, and sequence alignment in ProtoMatch provides researchers with a versatile system for performing both confirmatory and exploratory sequential data analyses (Sanderson & Fisher, 1994).

Eye Movements↗

Fast and reliable prediction of noncoding RNAs.

We report an efficient method for detecting functional RNAs. The approach, which combines comparative sequence analysis and structure prediction, already has yielded excellent results for a small number of aligned sequences and is suitable for large-scale genomic screens. It consists of two basic components: (i) a measure for RNA secondary structure conservation based on computing a consensus secondary structure, and (ii) a measure for thermodynamic stability, which, in the spirit of a z score, is normalized with respect to both sequence length and base composition but can be calculated without sampling from shuffled sequences. Functional RNA secondary structures can be identified in multiple sequence alignments with high sensitivity and high specificity. We demonstrate that this approach is not only much more accurate than previous methods but also significantly faster. The method is implemented in the program rnaz, which can be downloaded from www.tbi.univie.ac.at/~wash/RNAz. We screened all alignments of length n > or = 50 in the Comparative Regulatory Genomics database, which compiles conserved noncoding elements in upstream regions of orthologous genes from human, mouse, rat, Fugu, and zebrafish. We recovered all of the known noncoding RNAs and cis-acting elements with high significance and found compelling evidence for many other conserved RNA secondary structures not described so far to our knowledge.

Algorithms↗

Parameterization studies for the SAM and HMMER methods of hidden Markov model generation.

Multiple sequence alignment of distantly related viral proteins remains a challenge to all currently available alignment methods. The hidden Markov model approach offers a new, flexible method for the generation of multiple sequence alignments. The results of studies attempting to infer appropriate parameter constraints for the generation of de novo HMMs for globin, kinase, aspartic acid protease, and ribonuclease H sequences by both the SAM and HMMER methods are described.

Aspartic Acid Endopeptidases↗

Can three-dimensional contacts in protein structures be predicted by analysis of correlated mutations?

A method has been developed to detect pairs of positions with correlated mutations in protein multiple sequence alignments. The method is based on reconstruction of the phylogenetic tree for a set of sequences and statistical analysis of the distribution of mutations in the branches of the tree. The database of homology-derived protein structures (HSSP) is used as the source of multiple sequence alignments for proteins of known three-dimensional structure. We analyse pairs of positions with correlated mutations in 67 protein families and show quantitatively that the presence of such positions is a typical feature of protein families. A significant but weak tendency is observed for correlated residue pairs to be close in the three-dimensional structure. With further improvements, methods of this type may be useful for the prediction of residue--residue contacts and subsequent prediction of protein structure using distance geometry algorithms. In conclusion, we suggest a new experimental approach to protein structure determination in which selection of functional mutants after random mutagenesis and analysis of correlated mutations provide sufficient proximity constraints for calculation of the protein fold.

Amino Acid Sequence↗

Three-dimensional structure of kynureninase from Pseudomonas fluorescens.

Kynureninase [E.C. 3.7.1.3] is a pyridoxal-5'-phosphate (PLP)-dependent enzyme that catalyzes the hydrolytic cleavage of l-kynurenine to anthranilic acid and l-alanine. Sequence alignment with other PLP-dependent enzymes indicated that kynureninase is in subgroup IVa of the aminotransferases, along with nifS, CsdB, and serine-pyruvate aminotransferase, which suggests that kynureninase has an aminotransferase fold. Crystals of Pseudomonas fluorescens kynureninase were obtained, and the structure was solved by molecular replacement using the CsdB coordinates combined with multiple isomorphous heavy atom replacement. The coordinates were deposited in the PDB (ID code 1QZ9). The structure, refined to an R factor of 15.5% to 1.85 A resolution, is dimeric and has the aminotransferase fold. The structure also confirms the prediction from sequence alignment that Lys-227 is the PLP-binding residue in P. fluorescens kynureninase. The conserved Asp-201, expected for an aminotransferase fold, is located near the PLP nitrogen, but Asp-132 is also strictly conserved and at a similar distance from the pyridinium nitrogen. Mutagenesis of both conserved aspartic acids shows that both contribute equally to PLP binding, but Asp-201 has a greater role in catalysis. The structure shows that Tyr-226 donates a hydrogen bond to the phosphate of PLP. Unusual among PLP-dependent enzymes, Trp-256, which is also strictly conserved in kynureninases from bacteria to humans, donates a hydrogen bond to the phosphate through the indole N1-hydrogen.

Alanine↗

Homology modeling and active-site residues probing of the thermophilic Alicyclobacillus acidocaldarius esterase 2.

The moderate thermophilic eubacterium Alicyclobacillus (formerly Bacillus) acidocaldarius expresses a thermostable carboxylesterase (esterase 2) belonging to the hormone-sensitive lipase (HSL)-like group of the esterase/lipase family. Based on secondary structures predictions and a secondary structure-driven multiple sequence alignment with remote homologous protein of known three-dimensional (3D) structure, we previously hypothesized for this enzyme the alpha/beta-hydrolase fold typical of several lipases and esterases and identified Ser155, Asp252, and His282 as the putative members of the catalytic triad. In this paper we report the construction of a 3D model for this enzyme based on the structure of mouse acetylcholinesterase complexed with fasciculin. The model reveals the topological organization of the fold corroborating our predictions. As regarding the active-site residues, Ser155, Asp252, and His282 are located close to each other at hydrogen bond distances. Their catalytic role was here probed by biochemical and mutagenic studies. Moreover, on the basis of the secondary structure-driven multiple sequence alignment and the 3D structural model, a residue supposed important for catalysis, Gly84, was mutated to Ser. The activity of the mutated enzyme was drastically reduced. We propose that Gly84 is part of a putative "oxyanion hole" involved in the stabilization of the transition state similar to the C group of the esterase/lipase family.

Amino Acid Sequence↗

Sequence logos: a new way to display consensus sequences.

A graphical method is presented for displaying the patterns in a set of aligned sequences. The characters representing the sequence are stacked on top of each other for each position in the aligned sequences. The height of each letter is made proportional to its frequency, and the letters are sorted so the most common one is on top. The height of the entire stack is then adjusted to signify the information content of the sequences at that position. From these 'sequence logos', one can determine not only the consensus sequence but also the relative frequency of bases and the information content (measured in bits) at every position in a site or sequence. The logo displays both significant residues and subtle sequence patterns.

Amino Acid Sequence↗

Comparing sequences without using alignments: application to HIV/SIV subtyping.

BACKGROUND: In general, the construction of trees is based on sequence alignments. This procedure, however, leads to loss of informationwhen parts of sequence alignments (for instance ambiguous regions) are deleted before tree building. To overcome this difficulty, one of us previously introduced a new and rapid algorithm that calculates dissimilarity matrices between sequences without preliminary alignment. RESULTS: In this paper, HIV (Human Immunodeficiency Virus) and SIV (Simian Immunodeficiency Virus) sequence data are used to evaluate this method. The program produces tree topologies that are identical to those obtained by a combination of standard methods detailed in the HIV Sequence Compendium. Manual alignment editing is not necessary at any stage. Furthermore, only one user-specified parameter is needed for constructing trees. CONCLUSION: The extensive tests on HIV/SIV subtyping showed that the virus classifications produced by our method are in good agreement with our best taxonomic knowledge, even in non-coding LTR (Long Terminal Repeat) regions that are not tractable by regular alignment methods due to frequent duplications/insertions/deletions. Our method, however, is not limited to the HIV/SIV subtyping. It provides an alternative tree construction without a time-consuming aligning procedure.

Animals↗

Searching databases of conserved sequence regions by aligning protein multiple-alignments.

A general searching method for comparing multiple sequence alignments was developed to detect sequence relationships between conserved protein regions. Multiple alignments are treated as sequences of amino acid distributions and aligned by comparing pairs of such distributions. Four different comparison measures were tested and the Pearson correlation coefficient chosen. The method is sensitive, detecting weak sequence relationships between protein families. Relationships are detected beyond the range of conventional sequence database searches, illustrating the potential usefulness of the method. The previously undetected relation between flavoprotein subunits of two oxidoreductase families points to the potential active site in one of the families. The similarity between the bacterial RecA, DnaA and Rad51 protein families reveals a region in DnaA and Rad51 proteins likely to bind and unstack single-stranded DNA. Helix--turn--helix DNA binding domains from diverse proteins are readily detected and shown to be similar to each other. Glycosylasparaginase and gamma-glutamyltransferase enzymes are found to be similar in their proteolytic cleavage sites. The method has been fully implemented on the World Wide Web at URL: http://blocks.fhcrc.org/blocks-bin/LAMAvsearch.

Amino Acid Sequence↗

The hemocyanin from a living fossil, the cephalopod Nautilus pompilius: protein structure, gene organization, and evolution.

By electron microscopic and immunobiochemical analyses we have confirmed earlier evidence that Nautilus pompilius hemocyanin (NpH) is a ring-like decamer (M(r) = approximately 3.5 million), assembled from 10 identical copies of an approximately 350-kDa polypeptide. This subunit in turn is substructured into seven sequential covalently linked functional units of approximately 50 kDa each (FUs a-g). We have cloned and sequenced the cDNA encoding the complete polypeptide; it comprises 9198 bp and is subdivided into a 5' UTR of 58 bp, a 3' UTR of 365 bp, and an open reading frame for a signal peptide of 21 amino acids plus a polypeptide of 2903 amino acids (M(r) = 335,881). According to sequence alignments, the seven FUs of Nautilus hemocyanin directly correspond to the seven FU types of the previously sequenced hemocyanin "OdH" from the cephalopod Octopus dofleini. Thirteen potential N-glycosylation sites are distributed among the seven Nautilus hemocyanin FUs; the structural consequences of putatively attached glycans are discussed on the basis of the published X-ray structure for an Octopus dofleini and a Rapana thomasiana FU. Moreover, the complete gene structure of Nautilus hemocyanin was analyzed; it resembles that of Octopus hemocyanin with respect to linker introns but shows two internal introns that differ in position from the three internal introns of the Octopus hemocyanin gene. Multiple sequence alignments allowed calculation of a rather robust phylogenetic tree and a statistically firm molecular clock. This reveals that the last common ancestor of Nautilus and Octopus lived 415 +/- 24 million years ago, in close agreement with fossil records from the early Devonian.

Amino Acid Sequence↗

Exploring substrate binding and discrimination in fructose1, 6-bisphosphate and tagatose 1,6-bisphosphate aldolases.

Fructose 1,6-bisphosphate aldolase catalyses the reversible condensation of glycerone-P and glyceraldehyde 3-phosphate into fructose 1,6-bisphosphate. A recent structure of the Escherichia coli Class II fructose 1,6-bisphosphate aldolase [Hall, D.R., Leonard, G.A., Reed, C.D., Watt, C.I., Berry, A. & Hunter, W.N. (1999) J. Mol. Biol. 287, 383-394] in the presence of the transition state analogue phosphoglycolohydroxamate delineated the roles of individual amino acids in binding glycerone-P and in the initial proton abstraction steps of the mechanism. The X-ray structure has now been used, together with sequence alignments, site-directed mutagenesis and steady-state enzyme kinetics to extend these studies to map important residues in the binding of glyceraldehyde 3-phosphate. From these studies three residues (Asn35, Ser61 and Lys325) have been identified as important in catalysis. We show that mutation of Ser61 to alanine increases the Km value for fructose 1, 6-bisphosphate 16-fold and product inhibition studies indicate that this effect is manifested most strongly in the glyceraldehyde 3-phosphate binding pocket of the active site, demonstrating that Ser61 is involved in binding glyceraldehyde 3-phosphate. In contrast a S61T mutant had no effect on catalysis emphasizing the importance of an hydroxyl group for this role. Mutation of Asn35 (N35A) resulted in an enzyme with only 1.5% of the activity of the wild-type enzyme and different partial reactions indicate that this residue effects the binding of both triose substrates. Finally, mutation of Lys325 has a greater effect on catalysis than on binding, however, given the magnitude of the effects it is likely that it plays an indirect role in maintaining other critical residues in a catalytically competent conformation. Interestingly, despite its proximity to the active site and high sequence conservation, replacement of a fourth residue, Gln59 (Q59A) had no significant effect on the function of the enzyme. In a separate study to characterize the molecular basis of aldolase specificity, the agaY-encoded tagatose 1,6-bisphosphate aldolase of E. coli was cloned, expressed and kinetically characterized. Our studies showed that the two aldolases are highly discriminating between the diastereoisomers fructose bisphosphate and tagatose bisphosphate, each enzyme preferring its cognate substrate by a factor of 300-1500-fold. This produces an overall discrimination factor of almost 5 x 105 between the two enzymes. Using the X-ray structure of the fructose 1,6-bisphosphate aldolase and multiple sequence alignments, several residues were identified, which are highly conserved and are in the vicinity of the active site. These residues might potentially be important in substrate recognition. As a consequence, nine mutations were made in attempts to switch the specificity of the fructose 1,6-bisphosphate aldolase to that of the tagatose 1,6-bisphosphate aldolase and the effect on substrate discrimination was evaluated. Surprisingly, despite making multiple changes in the active site, many of which abolished fructose 1, 6-bisphosphate aldolase activity, no switch in specificity was observed. This highlights the complexity of enzyme catalysis in this family of enzymes, and points to the need for further structural studies before we fully understand the subtleties of the shaping of the active site for complementarity to the cognate substrate.

Aldehyde-Lyases↗