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 1,189 records · Page 66Linked to original sources
PubMed2006

Molecular dynamics and circular dichroism studies of human and rat C-peptides.

Proinsulin C-peptide has been recently described as an endogenous peptide hormone, responsible for important physiological functions others than its role in proinsulin processing. Accumulating evidences that C-peptide exerts beneficial effects in the treatment of long term complications of patients with type 1 diabetes mellitus indicate that this molecule may be administered together with insulin in future therapies. Despite its clear pharmacological interest, the secondary and three-dimensional (3D) structures of human C-peptide are still points of controversy. In the present work we report molecular dynamics (MD) simulations of human, rat I and rat II C-peptides. A common experimental strategy applied to all peptides consisted of homology building followed by multinanosecond MD simulations in vacuum and water. Circular dichroism (CD) experiments of each peptide in the absence and presence of 2,2,2-trifluoroethanol (TFE) were performed to support validation of the theoretical models. A multiple sequence alignment of 23 known mammalian C-peptides was constructed to identify significant conserved sites that would be important for the maintenance of secondary and tertiary structures. The analysis of the molecular dynamics trajectories for the human, rat I and rat II molecules have shown quite different general behavior, being the human C-peptide more flexible than the two others. Human and rat C-peptides exhibit very stable turn-like structures at the middle and C-terminal regions, which have been described as potential active sites of C-peptides. Human C-peptide also presented a short alpha-helix throughout the MD, which was not found in the rat molecules. CD data is in very good agreement with the MD results and both methods were able to identify a greater structural stability and potential in rat C-peptides when compared to the human C-peptide. The simulation results are discussed and validated in the light of multiple sequence alignment, recent experimental data from the literature and our own CD experiments.

Amino Acid Sequence↗
PubMed2002

Empirical determination of effective gap penalties for sequence comparison.

MOTIVATION: No general theory guides the selection of gap penalties for local sequence alignment. We empirically determined the most effective gap penalties for protein sequence similarity searches with substitution matrices over a range of target evolutionary distances from 20 to 200 Point Accepted Mutations (PAMs). RESULTS: We embedded real and simulated homologs of protein sequences into a database and searched the database to determine the gap penalties that produced the best statistical significance for the distant homologs. The most effective penalty for the first residue in a gap (q+r) changes as a function of evolutionary distance, while the gap extension penalty for additional residues (r) does not. For these data, the optimal gap penalties for a given matrix scaled in 1/3 bit units (e.g. BLOSUM50, PAM200) are q=25-0.1 * (target PAM distance), r=5. Our results provide an empirical basis for selection of gap penalties and demonstrate how optimal gap penalties behave as a function of the target evolutionary distance of the substitution matrix. These gap penalties can improve expectation values by at least one order of magnitude when searching with short sequences, and improve the alignment of proteins containing short sequences repeated in tandem.

Amino Acid Sequence↗
PubMed2006

PROTOGENE: turning amino acid alignments into bona fide CDS nucleotide alignments.

We describe Protogene, a server that can turn a protein multiple sequence alignment into the equivalent alignment of the original gene coding DNA. Protogene relies on a pipeline where every initial protein sequence is BLASTed against RefSeq or NR. The annotation associated with potential matches is used to identify the gene sequence. This gene sequence is then aligned with the query protein using Exonerate in order to extract a coding nucleotide sequence matching the original protein. Protogene can handle protein fragments and will return every CDS coding for a given protein, even if they occur in different genomes. Protogene is available from http://www.tcoffee.org/.

Base Sequence↗
PubMed2002

Comparative modelling of cytochromes P450.

The superfamily of enzymes known as the cytochromes P450 (P450s) comprises a wide-ranging class of proteins with diverse functions. They are known, amongst other things, to be involved in the hormonal regulation of metabolism and reproduction, as well as having a major clinical significance through their association with diseases such as cancer, diabetes and hepatitis. Knowledge of the three-dimensional (3D) structure of a protein gives insight into its function. The 3D structures of P450s are therefore of considerable scientific interest. A number of high-resolution structures of P450s have been determined by X-ray crystallography and studies of these structures have provided valuable insights into the mechanism of these enzymes. Only one of these structures is mammalian and as yet there is no structural information on human P450s in the public domain. Until such a structure is solved it is necessary to employ alternative methods to gain structural insight into how human P450s perform their biological function. Here we report on the use of comparative modelling to predict the structure of human P450s based on knowledge of their amino acid sequences plus the 3D structures of other (not human) P450s. As an illustrative example of these techniques we have modelled the structure of P450 2C5 using five bacterial P450 structures as templates. We examine the importance of selecting suitable templates, obtaining a good amino acid sequence alignment, and evaluating the models generated. To improve the quality of the models an iterative cycle of sequence alignment, model building, and model evaluation is employed. The result is a model with excellent stereochemistry, good amino acid side chain environment properties, and a Calpha trace similar to the crystal structure.

Amino Acid Sequence↗
PubMed1991

An evolutionary model for maximum likelihood alignment of DNA sequences.

Most algorithms for the alignment of biological sequences are not derived from an evolutionary model. Consequently, these alignment algorithms lack a strong statistical basis. A maximum likelihood method for the alignment of two DNA sequences is presented. This method is based upon a statistical model of DNA sequence evolution for which we have obtained explicit transition probabilities. The evolutionary model can also be used as the basis of procedures that estimate the evolutionary parameters relevant to a pair of unaligned DNA sequences. A parameter-estimation approach which takes into account all possible alignments between two sequences is introduced; the danger of estimating evolutionary parameters from a single alignment is discussed.

Algorithms↗
PubMed2006

Prediction of transporter family from protein sequence by support vector machine approach.

Transporters play key roles in cellular transport and metabolic processes, and in facilitating drug delivery and excretion. These proteins are classified into families based on the transporter classification (TC) system. Determination of the TC family of transporters facilitates the study of their cellular and pharmacological functions. Methods for predicting TC family without sequence alignments or clustering are particularly useful for studying novel transporters whose function cannot be determined by sequence similarity. This work explores the use of a machine learning method, support vector machines (SVMs), for predicting the family of transporters from their sequence without the use of sequence similarity. A total of 10,636 transporters in 13 TC subclasses, 1914 transporters in eight TC families, and 168,341 nontransporter proteins are used to train and test the SVM prediction system. Testing results by using a separate set of 4351 transporters and 83,151 nontransporter proteins show that the overall accuracy for predicting members of these TC subclasses and families is 83.4% and 88.0%, respectively, and that of nonmembers is 99.3% and 96.6%, respectively. The accuracies for predicting members and nonmembers of individual TC subclasses are in the range of 70.7-96.1% and 97.6-99.9%, respectively, and those of individual TC families are in the range of 60.6-97.1% and 91.5-99.4%, respectively. A further test by using 26,139 transmembrane proteins outside each of the 13 TC subclasses shows that 90.4-99.6% of these are correctly predicted. Our study suggests that the SVM is potentially useful for facilitating functional study of transporters irrespective of sequence similarity.

Amino Acid Sequence↗
PubMed1998

Identification of two new members of the mouse prolactin gene family.

The prolactin (PRL) family consists of a collection of genes expressed in the uterus, placenta, and anterior pituitary. These cytokines/hormones participate in the control of maternal-fetal adaptations to pregnancy. In this report, we establish the presence of two new members of the mouse PRL family. Novel expressed sequence tags (ESTs) with significant homology to PRL were isolated from embryonic, ectoplacental cone, and placental cDNA libraries. The cDNAs were sequenced and compared to other members of the PRL family. The two new cDNAs were assigned to the PRL family based on sequence homology and were referred to as PRL-like protein-F (PLP-F) and PRL-like protein-G (PLP-G). PLP-F cDNA encodes for a predicted 267 amino acid protein containing a 30 amino acid signal peptide and three putative N-linked glycosylation sites. PLP-G cDNA encodes for a predicted 266 amino acid protein containing a 30 amino acid signal peptide and six putative N-linked glycosylation sites. Sequence alignments of these proteins with other members of the PRL family suggest some unique features. Both sequences contain an extra amino acid segment located between exons two and three of the prototypical PRL gene and a nine amino acid carboxy terminal extension. PLP-F contained an additional 15 amino acid region situated between exons four and five of the prototypical PRL gene. Both PLP-F and PLP-G mRNAs were expressed in the placenta but not in other tissues (uterus, brain, thymus, heart, lung, diaphragm, liver, kidney, and ovary). In summary, the two newly identified members share approximately 50% amino acid sequence identity, are specifically expressed in the placenta, and represent a new subfamily within the PRL family.

Amino Acid Sequence↗
PubMed2002

In vitro selection of high affinity HspR-binding sites within the genome of Helicobacter pylori.

The major chaperone genes of Helicobacter pylori are negatively regulated by HspR, a homologue of the repressor of the dnaK operon of Streptomyces coelicolor. Using an in vitro selection and amplification approach we identified two new chromosomal binding sites of the HspR protein. Both binding sites were characterized by footprinting analysis with purified HspR protein. Intriguingly, these HspR binding sites are located at the 3prime prime or minute ends of two genes coding for predicted proteins with functions unrelated to those of chaperones. This suggests that H. pylori HspR may regulate the expression of genes encoding proteins with diverse functions. Nucleotide sequence alignment of HspR-binding sites highlights conserved nucleotides extending outside the previously proposed consensus binding sequence with structural features predicting geometry of HspR binding as an oligomer.

Bacterial Proteins↗
PubMed2003

Serine proteases and their homologs in the Drosophila melanogaster genome: an initial analysis of sequence conservation and phylogenetic relationships.

Serine proteases (SPs) and serine protease homologs (SPHs) constitute the second largest family of genes in the Drosophila melanogaster genome. Eighty-four SPs comprise less than 300 amino acid residues, and a significant portion of them are probably digestive enzymes. Some larger SPs may contain one or more regions important for protein-protein interactions, including clip domains, low-density lipoprotein receptor class A repeats, and scavenger receptor cysteine-rich domains. We identified 37 clusters of SP or SPH genes, which probably evolved from relatively recent gene duplication and sequence divergence. A majority of the SPs may be trypsin-like and activated by cleavage after a specific arginine or lysine residue. Among the 147 SPs and 57 SPHs studied, 24 SPs and 13 SPHs contain at least one regulatory clip domain. A multiple sequence alignment of the clip domains provided further information on structural conservation of these regulatory modules. Detailed sequence comparison led to an improved classification system for SPs containing clip domains. These analyses have established a framework of information about evolutionary relationships among the Drosophila SPs and SPHs, which may facilitate research on these proteins as well as homologous molecules from other invertebrate species.

Amino Acid Sequence↗
PubMed1999

Sequence analysis of genes encoding ribosomal proteins of amitochondriate protists: L1 of Trichomonas vaginalis and L29 of Giardia lamblia.

Two genes encoding the ribosomal proteins were cloned and sequenced from amitochondriate protists, L1 (L10a in mammalian nomenclature) from Trichomonas vaginalis and L29 (L35 in mammalian nomenclature) from Giardia lamblia. The deduced amino acid sequences were analyzed by sequence alignments and phylogenetic reconstructions. Both the T. vaginalis L1 and the G. lamblia L29 displayed eukaryotic sequence features, when compared with all the homologs from the three primary kingdoms.

Amino Acid Sequence↗
PubMed1998

Direct binding of p130(Cas) to the guanine nucleotide exchange factor C3G.

p130(Cas) (Cas; crk-associated substrate) belongs to a new family of docking molecules. It contains one Src homology (SH) 3 domain in its amino-terminal region followed by a region containing binding motifs for SH2 and SH3 domains. To gain further insight into Cas signaling we used the SH3 domain of Cas in a two-hybrid screen to search a human placenta library for binding partners. The screen confirmed a previous finding of its binding to the focal adhesion kinase (FAK) but also identified C3G, a guanine nucleotide exchange factor. We found direct interaction between Cas and C3G in vitro and in vivo. A series of analysis with C3G deletion mutants revealed a proline-rich Cas-binding site (Ala0-Pro1-Pro2-Lys3-Pro4-Pro5-Leu6-Pro7) located NH2-terminal to the previously characterized Crk binding motifs in C3G. Mutagenesis studies showed that Pro1, Lys3, and Pro4 within the ligand-binding site are critical for high affinity interaction. These results, combined with sequence alignments of proline-rich binding elements from proteins known for Cas binding, define the consensus sequence XXPXKPX which is recognized by the CasSH3 domain. Cas shows structural characteristics of a docking molecule and may serve to bring C3G to specific compartments within the cell.

Amino Acid Sequence↗
PubMed1997

The proofreading domain of Escherichia coli DNA polymerase I and other DNA and/or RNA exonuclease domains.

Prior sequence analysis studies have suggested that bacterial ribonuclease (RNase) Ds comprise a complete domain that is found also in Homo sapiens polymyositis-scleroderma overlap syndrome 100 kDa autoantigen and Werner syndrome protein. This RNase D 3'-->5' exoribonuclease domain was predicted to have a structure and mechanism of action similar to the 3'-->5' exodeoxyibonuclease (proofreading) domain of DNA polymerases. Here, hidden Markov model (HMM) and phylogenetic studies have been used to identify and characterise other sequences that may possess this exonuclease domain. Results indicate that it is also present in the RNase T family; Borrelia burgdorferi P93 protein, an immunodominant antigen in Lyme disease; bacteriophage T4 dexA and Escherichia coli exonuclease I, processive 3'-->5' exodeoxyribonucleases that degrade single-stranded DNA; Bacillus subtilis dinG, a probable helicase involved in DNA repair and possibly replication, and peptide synthase 1; Saccharomyces cerevisiae Pab1p-dependent poly(A) nuclease PAN2 subunit, required for shortening mRNA poly(A) tails; Caenorhabditis elegans and Mus musculus CAF1, transcription factor CCR4-associated factor 1; Xenopus laevis XPMC2, prevention of mitotic catastrophe in fission yeast; Drosophila melanogaster egalitarian, oocyte specification and axis determination, and exuperantia, establishment of oocyte polarity; H.sapiens HEM45, expressed in tumour cell lines and uterus and regulated by oestrogen; and 31 open reading frames including one in Methanococcus jannaschii . Examination of a multiple sequence alignment and two three-dimensional structures of proofreading domains has allowed definition of the core sequence, structural and functional elements of this exonuclease domain.

Amino Acid Sequence↗
PubMed2006

A Consensus Data Mining secondary structure prediction by combining GOR V and Fragment Database Mining.

The major aim of tertiary structure prediction is to obtain protein models with the highest possible accuracy. Fold recognition, homology modeling, and de novo prediction methods typically use predicted secondary structures as input, and all of these methods may significantly benefit from more accurate secondary structure predictions. Although there are many different secondary structure prediction methods available in the literature, their cross-validated prediction accuracy is generally <80%. In order to increase the prediction accuracy, we developed a novel hybrid algorithm called Consensus Data Mining (CDM) that combines our two previous successful methods: (1) Fragment Database Mining (FDM), which exploits the Protein Data Bank structures, and (2) GOR V, which is based on information theory, Bayesian statistics, and multiple sequence alignments (MSA). In CDM, the target sequence is dissected into smaller fragments that are compared with fragments obtained from related sequences in the PDB. For fragments with a sequence identity above a certain sequence identity threshold, the FDM method is applied for the prediction. The remainder of the fragments are predicted by GOR V. The results of the CDM are provided as a function of the upper sequence identities of aligned fragments and the sequence identity threshold. We observe that the value 50% is the optimum sequence identity threshold, and that the accuracy of the CDM method measured by Q(3) ranges from 67.5% to 93.2%, depending on the availability of known structural fragments with sufficiently high sequence identity. As the Protein Data Bank grows, it is anticipated that this consensus method will improve because it will rely more upon the structural fragments.

Algorithms↗
PubMed2000

Extensive sequence homology between the mycobacterium leprae LSR (12 kDa) antigen and its Mycobacterium tuberculosis counterpart.

The Mycobacterium leprae LSR (12 kDa) protein antigen has been reported to mimic whole cell M. leprae in T cell responses across the leprosy spectrum. In addition, B cell responses to specific sequences within the LSR antigen have been shown to be associated with immunopathological responses in leprosy patients with erythema nodosum leprosum. We have in the present study applied the M. leprae LSR DNA sequence as query to search for the presence of homologous genes within the recently completed Mycobacterium tuberculosis genome database (Sanger Centre, UK). By using the BLASTN search tool, a homologous M. tuberculosis open reading frame (336 bp), encoding a protein antigen of 12.1 kDa, was identified within the cosmid MTCY07H7B.25. The gene is designated Rv3597c within the M. tuberculosis H37Rv genome. Sequence alignment revealed 93% identity between the M. leprae and M. tuberculosis antigens at the amino acid sequence level. The finding that some B and T cell epitopes were localized to regions with amino acid substitutions may account for the putative differential responsiveness to this antigen in tuberculosis and leprosy.

Amino Acid Sequence↗
PubMed2002

Development of an 18S rRNA gene-targeted PCR-based diagnostic for the blue crab parasite Hematodinium sp..

The 18S rRNA gene from Hematodinum sp., a parasitic dinoflagellate that infects blue crabs, was amplified, cloned, and sequenced. The sequence showed a high similarity (95% at the nucleotide level) to sequences obtained from other dinoflagellate species, including both free-living and symbiotic species. Sequence similarity was much lower when compared with parasites of other marine invertebrates with similar life histories and with the 18S rRNA gene from the blue crab. Based on comparison of sequence alignments between Hematodinium, other dinoflagellate species, protozoan pathogens of oysters, and blue crab 18S rRNA gene sequences, 2 sets of PCR primers that specifically amplified fragments of the Hematodinium 18S rRNA gene were developed and tested. One of these primer sets (Hemat-F-1487 and Hemat-R-1654) amplified a 187 bp fragment that could be used routinely as a diagnostic test for the presence of Hematodinium in hemolymph from blue crabs. This fragment was consistently amplified from genomic DNA extracted from hemolymph of Hematodinium infected blue crabs. Comparison between the PCR technique and standard histological examination indicated that the PCR technique was reliable and provided 1000 times more sensitivity than the histological methods. The sensitivity of the PCR diagnostic was estimated to be one parasite cell among 300,000 crab hemocytes. Preliminary studies using the PCR diagnostic technique suggest that Hematodinium sp. is absent in crabs collected from waters with low salinity (5 to 10 ppt), but common in crabs from higher salinity environments in estuarine waters from southeastern Georgia (USA).

Animals↗
PubMed2000

Fast, optimal alignment of three sequences using linear gap costs.

Alignment algorithms can be used to infer a relationship between sequences when the true relationship is unknown. Simple alignment algorithms use a cost function that gives a fixed cost to each possible point mutation-mismatch, deletion, insertion. These algorithms tend to find optimal alignments that have many small gaps. It is more biologically plausible to have fewer longer gaps rather than many small gaps in an alignment. To address this issue, linear gap cost algorithms are in common use for aligning biological sequence data. More reliable inferences are obtained by aligning more than two sequences at a time. The obvious dynamic programming algorithm for optimally aligning k sequences of length n runs in O(n(k)) time. This is impractical if k>/=3 and n is of any reasonable length. Thus, for this problem there are many heuristics for aligning k sequences, however, they are not guaranteed to find an optimal alignment. In this paper, we present a new algorithm guaranteed to find the optimal alignment for three sequences using linear gap costs. This gives the same results as the dynamic programming algorithm for three sequences, but typically does so much more quickly. It is particularly fast when the (three-way) edit distance is small. Our algorithm uses a speed-up technique based on Ukkonen's greedy algorithm (Ukkonen, 1983) which he presented for two sequences and simple costs.

Algorithms↗
PubMed2003

Consensus alignment for reliable framework prediction in homology modeling.

MOTIVATION: Even the best sequence alignment methods frequently fail to correctly identify the framework regions for which backbones can be copied from the template into the target structure. Since the underprediction and, more significantly, the overprediction of these regions reduces the quality of the final model, it is of prime importance to attain as much as possible of the true structural alignment between target and template. RESULTS: We have developed an algorithm called Consensus that consistently provides a high quality alignment for comparative modeling. The method follows from a benchmark analysis of the 3D models generated by ten alignment techniques for a set of 79 homologous protein structure pairs. For 20-to-40% of the targets, these methods yield models with at least 6 A root mean square deviation (RMSD) from the native structure. We have selected the top five performing methods, and developed a consensus algorithm to generate an improved alignment. By building on the individual strength of each method, a set of criteria was implemented to remove the alignment segments that are likely to correspond to structurally dissimilar regions. The automated algorithm was validated on a different set of 48 protein pairs, resulting in 2.2 A average RMSD for the predicted models, and only four cases in which the RMSD exceeded 3 A. The average length of the alignments was about 75% of that found by standard structural superposition methods. The performance of Consensus was consistent from 2 to 32% target-template sequence identity, and hence it can be used for accurate prediction of framework regions in homology modeling.

Algorithms↗