Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Multiple sequence alignment”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28Linked to original sources

Isolation and characterization of pcp, a gene encoding a pyrrolidone carboxyl peptidase in Staphylococcus aureus.

The pcp gene, encoding a pyrrolidone carboxyl peptidase (PYRase), was cloned from a lambda GT11 genomic library prepared from Staphylococcus aureus FDA 574 and sequenced. The pcp gene is located 740 bp downstream from cna, a gene that encodes a collagen-binding adhesin in S. aureus. S. aureus pcp encodes a 212-amino-acid (aa) polypeptide. The pcp gene was overexpressed in Escherichia coli and the PYRase purified to homogeneity. The recombinant enzyme exhibited biological activity, as determined using the chromogenic substrate L-pyroglutamyl-beta-napthylamide. Biochemical analysis of the PYRase using thiol-blocking chemicals suggested that the enzyme belongs to the cysteine peptidase family. Moreover, multiple sequence alignment revealed a high degree of similarity to previously described bacterial PYRases. This family of peptidases has been used to selectively remove the N-terminal pyrrolidone carboxylic acid residue found on certain blocked proteins and peptides prior to aa sequencing. However, the exact biological role of PYRases has yet to be elucidated.

Amino Acid Sequence↗

Cloning and characterisation of a gene from Plasmodium vivax and P. knowlesi: homology with valine-tRNA synthetase.

We have previously described a lambdagt11 clone detected by immune screening with a monoclonal antibody (mAb) A12. This mAb is capable of completely blocking Plasmodium vivax transmission in the mosquito vector. An epitope recognised by A12 was mapped to six amino acids (aa) within the translated sequence of this clone. Here, we describe the complete sequence of the gene within which we mapped this epitope. Surprisingly, the translated sequence of the full-length open reading frame shows homology with that of valine-tRNA synthetases (Val-tRS) from other organisms. DNA cross-hybridisation with several of these species was observed by Southern blot. In addition, the corresponding gene has been obtained from the closely related simian malaria parasite, P. knowlesi. The two aa sequences show 66% identity and yet are very divergent from other Val-tRS sequences, apart from conserved blocks related to functional activity. Multiple sequence alignments reflect this dichotomy, as do predicted differences in antigenicity.

Amino Acid Sequence↗

Structure prediction of proteins--where are we now?

Although the 'structure from sequence' prediction problem remains fundamentally unsolved, new and promising methods in one, two and three dimensions have reopened the field. Significantly improved one-dimensional prediction of secondary structure from multiple sequence alignments is now in routine use. In the two-dimensional approach, inter-residue contacts can be detected by analysis of correlated mutations, albeit with low accuracy. Finally, three-dimensional methods, in which pseudopotentials or information values are derived from the databases, are proving their value for distinguishing between correct and incorrect models.

Amino Acid Sequence↗

Characterization and modeling of membrane proteins using sequence analysis.

The current libraries of amino acid sequences of membrane proteins are a valuable resource for the analysis of elements common to these proteins. Multiple-sequence alignment techniques and the identification of conserved features of transmembrane segments have improved the prediction of membrane protein topology. Molecular modeling in combination with structural studies or site-directed mutagenesis is proving to be a powerful link between theory and experiment. Unfortunately, the number of high-resolution structures of intrinsic membrane proteins, although increased recently, presents a restricted and perhaps biased view of membrane protein structure.

Amino Acid Sequence↗

Progress in protein structure prediction?

Prediction of protein secondary structure is an old problem and progress has been slow. Recently, spectacular success has been claimed in the blind prediction of the catalytic subunit of the cAMP-dependent protein kinase. When predictions in this and other test cases are assessed critically, some claims of prediction success turn out to be exaggerated, but a kernel of real progress remains: protein structure prediction can be improved substantially when a family of related sequences is available. Enough so that molecular biologists equipped with a new amino acid sequence and a multiple sequence alignment in hand may be tempted to test the new prediction methods.

Amino Acid Sequence↗

Sequence analysis of cytochrome bd oxidase suggests a revised topology for subunit I.

Numerous sequences of the cytochrome bd quinol oxidase (cytochrome bd) have recently become available for analysis. The analysis has revealed a small number of conserved residues, a new topology for subunit I and a phylogenetic tree involving extensive horizontal gene transfer. There are 20 conserved residues in subunit I and two in subunit II. Algorithms utilizing multiple sequence alignments predicted a revised topology for cytochrome bd, adding two transmembrane helices to subunit I to the seven that were previously indicated by the analysis of the sequence of the oxidase from E. coli. This revised topology has the effect of relocating the N-terminus and C-terminus to the periplasmic and cytoplasmic sides of the membrane, respectively. The new topology repositions I-H19, the putative ligand for heme b595, close to the periplasmic edge of the membrane, which suggests that the heme b595/heme d active site of the oxidase is located near the outer (periplasmic) surface of the membrane. The most highly conserved region of the sequence of subunit I contains the sequence GRQPW and is located in a predicted periplasmic loop connecting the eighth and ninth transmembrane helices. The potential importance of this region of the protein was previously unsuspected, and it may participate in the binding of either quinol or heme d. There are two very highly conserved glutamates in subunit I, E99 and E107, within the third transmembrane helix (E. coli cytochrome bd-I numbering). It is speculated that these glutamates may be part of a proton channel leading from the cytoplasmic side of the membrane to the heme d oxygen-reactive site, now placed near the periplasmic surface. The revised topology and newly revealed conserved residues provide a clear basis for further experimental tests of these hypotheses. Phylogenetic analysis of the new sequences of cytochrome bd reveals considerable deviation from the 16sRNA tree, suggesting that a large amount of horizontal gene transfer has occurred in the evolution of cytochrome bd.

Amino Acid Sequence↗

A new multidrug resistance protein at the blood-brain barrier.

Porcine brain capillary endothelial cells (PBCEC) cultured in serum-free and hydrocortisone supplemented medium are characterised by high transendothelial electrical resistances and low cell monolayer permeabilities for small solutes very similar to the blood-brain barrier (BBB) in vivo. Differential screening of a subtracted cDNA library disclosed a 2.1-kb mRNA that is overexpressed in hydrocortisone treated PBCEC relative to untreated cells. The mRNA encodes a 656-aa member of the ATP-binding cassette (ABC) superfamily of transporters that we named brain multidrug resistance protein (BMDP). Phylogenetic analysis and multiple sequence alignment showed that porcine BMDP is most related to the human and mouse breast cancer resistance protein (BCRP). Northern blot analysis revealed that BMDP is expressed in brain tissue in vivo and was predominantly localised within the endothelial cells isolated from brain capillaries. Thus, we identified a new transport protein at the BBB that might play an important role in the exclusion of xenobiotics from the brain.

ATP Binding Cassette Transporter, Subfamily B↗

Roles of conserved methionine residues in tobacco acetolactate synthase.

Acetolactate synthase (ALS) catalyzes the first common step in the biosynthesis of valine, leucine, and isoleucine. ALS is the target of several classes of herbicides, including the sulfonylureas, the imidazolinones, and the triazolopyrimidines. The conserved methionine residues of ALS from plants were identified by multiple sequence alignment using ClustalW. The alignment of 17 ALS sequences from plants revealed 149 identical residues, seven of which were methionine residues. The roles of three well-conserved methionine residues (M350, M512, and M569) in tobacco ALS were determined using site-directed mutagenesis. The mutation of M350V, M512V, and M569V inactivated the enzyme and abolished the binding affinity for cofactor FAD. Nevertheless, the secondary structure of each of the mutants determined by CD spectrum was not affected significantly by the mutation. Both M350C and M569C mutants were strongly resistant to three classes of herbicides, Londax (a sulfonylurea), Cadre (an imidazolinone), and TP (a triazolopyrimidine), while M512C mutant did not show a significant resistance to the herbicides. The mutant M350C was more sensitive to pH change, while the mutant M569C showed a profile for pH dependence activity similar to that of wild type. These results suggest that M512 residue is likely located at or near the active site, and that M350 and M569 residues are probably located at the overlapping region between the active site and a common herbicide binding site.

Acetolactate Synthase↗

Evolutionary relationship between K(+) channels and symporters.

The hypothesis is presented that at least four families of putative K(+) symporter proteins, Trk and KtrAB from prokaryotes, Trk1,2 from fungi, and HKT1 from wheat, evolved from bacterial K(+) channel proteins. Details of this hypothesis are organized around the recently determined crystal structure of a bacterial K(+) channel: i. e., KcsA from Streptomyces lividans. Each of the four identical subunits of this channel has two fully transmembrane helices (designated M1 and M2), plus an intervening hairpin segment that determines the ion selectivity (designated P). The symporter sequences appear to contain four sequential M1-P-M2 motifs (MPM), which are likely to have arisen from gene duplication and fusion of the single MPM motif of a bacterial K(+) channel subunit. The homology of MPM motifs is supported by a statistical comparison of the numerical profiles derived from multiple sequence alignments formed for each protein family. Furthermore, these quantitative results indicate that the KtrAB family of symporters has remained closest to the single-MPM ancestor protein. Strong sequence evidence is also found for homology between the cytoplasmic C-terminus of numerous bacterial K(+) channels and the cytoplasm-resident TrkA and KtrA subunits of the Trk and KtrAB symporters, which in turn are homologous to known dinucleotide-binding domains of other proteins. The case for homology between bacterial K(+) channels and the four families of K(+) symporters is further supported by the accompanying manuscript, in which the patterns of residue conservation are demonstrated to be similar to each other and consistent with the known 3D structure of the KcsA K(+) channel.

Amino Acid Sequence↗

Pseudomonas fluorescens mannitol 2-dehydrogenase and the family of polyol-specific long-chain dehydrogenases/reductases: sequence-based classification and analysis of structure-function relationships.

Multiple sequence alignment and analysis of evolutionary relationships have been used to characterize a family of polyol-specific long-chain dehydrogenases/reductases (PSLDRs). At the present time, 66 known and putative NAD(P)H-dependent oxidoreductases of mainly prokaryotic origin and between 357 and 544 amino acids in length constitute this family. The family is shown to include D-mannitol 2-dehydrogenase, D-mannonate 5-oxidoreductase, D-altronate 5-oxidoreductase, D-arabinitol 4-dehydrogenase, and D-mannitol-1-phosphate 5-dehydrogenase which form individual sub-families (defined by internal sequence identity of >/=30%) having distant origin and divergent substrate specificity but clearly displaying entire-chain relationship. When all forms are aligned, only three residues, Gly-33, Asp-230, and Lys-295 (in the numbering of Pseudomonas fluorescens D-mannitol 2-dehydrogenase (PsM2DH)) are strictly conserved. By combining sequence alignment with the known structure of PsM2DH and results from site-directed mutagenesis, we have developed a structure/function analysis for the family. Gly-33 is in the N-terminal coenzyme-binding domain and part of a nucleotide fingerprint region for the family, and Asp-230 and Lys-295 are at an interdomain segment contributing to the active site in which the lysine likely functions as the catalytic general acid/base. PSLDRs do not require a metal cofactor for activity and are specific for transferring the 4-pro-S hydrogen from NAD(P)H. Comparisons reveal that the core part of the two-domain fold has been conserved throughout all family members, perhaps reflecting the recruitment of a stable oxidoreductase structure and extensive trimming thereof to acquire functional properties specific to each sub-family. They also identify interactions that define the chemical mechanism of oxidoreduction and likely contribute to substrate and co-substrate specificities and are thus relevant for protein engineering.

Amino Acid Sequence↗

The molecular peculiarities of catalase-peroxidases.

In developing ideas of how protein structure modifies haem reactivity, the activity of Class I of the plant peroxidase superfamily (including cytochrome c peroxidase, ascorbate peroxidase and catalase-peroxidases (KatGs)) is an exciting field of research. Despite striking sequence homologies, there are dramatic differences in catalytic activity and substrate specificity with KatGs being the only member with substantial catalase activity. Based on multiple sequence alignment performed for Class I peroxidases, we present a hypothesis for the pronounced catalase activity of KatGs. In their catalytic domains KatGs are shown to possess three large insertions, two of them are typical for KatGs showing highly conserved sequence patterns. Besides an extra C-terminal copy of the ancestral hydroperoxidase gene resulting from gene duplication, these two large loops are likely to control the orientation of both the haem group and of essential residues in the active site. They seem to modulate the access of substrates to the prosthetic group at the distal side as well as the flexibility and character of the bond between the proximal histidine and the ferric iron. The hypothesis presented opens new possibilities in the rational engineering of peroxidases.

Amino Acid Sequence↗

A refined structure of human aquaporin-1.

A refined structure of the human water channel aquaporin-1 is presented. The model rests on the high resolution X-ray structure of the homologous bacterial glycerol transporter GlpF, electron crystallographic data at 3.8 A resolution and a multiple sequence alignment of the aquaporin superfamily. The crystallographic R and free R values (36.7% and 37.8%) for the refined structure are significantly lower than for previous models. Improved geometry and enhanced stability in molecular dynamics simulations demonstrate a significant improvement of the aquaporin-1 structure. Comparison with previous aquaporin-1 models shows significant differences, not only in the loop regions, but also in the core of the water channel.

Aquaporin 1↗

BTK-2, a new inhibitor of the Kv1.1 potassium channel purified from Indian scorpion Buthus tamulus.

A novel inhibitor of voltage-gated potassium channel was isolated and purified to homogeneity from the venom of the red scorpion Buthus tamulus. The primary sequence of this toxin, named BTK-2, as determined by peptide sequencing shows that it has 32 amino acid residues with six conserved cysteines. The molecular weight of the toxin was found to be 3452 Da. It was found to block the human potassium channel hKv1.1 (IC(50)=4.6 microM). BTK-2 shows 40-70% sequence similarity to the family of the short-chain toxins that specifically block potassium channels. Multiple sequence alignment helps to categorize the toxin in the ninth subfamily of the K+ channel blockers. The modeled structure of BTK-2 shows an alpha/beta scaffold similar to those of the other short scorpion toxins. Comparative analysis of the structure with those of the other toxins helps to identify the possible structure-function relationship that leads to the difference in the specificity of BTK-2 from that of the other scorpion toxins. The toxin can also be used to study the assembly of the hKv1.1 channel.

Amino Acid Sequence↗

Diversification and evolution of L-myo-inositol 1-phosphate synthase.

L-myo-Inositol 1-phosphate synthase (MIPS, EC 5.5.1.4), the key enzyme in the inositol and phosphoinositide biosynthetic pathway, is present throughout evolutionarily diverse organisms and is considered an ancient protein/gene. Analysis by multiple sequence alignment, phylogenetic tree generation and comparison of newly determined crystal structures provides new insight into the origin and evolutionary relationships among the various MIPS proteins/genes. The evolution of the MIPS protein/gene among the prokaryotes seems more diverse and complex than amongst the eukaryotes. However, conservation of a 'core catalytic structure' among the MIPS proteins implies an essential function of the enzyme in cellular metabolism throughout the biological kingdom.

Amino Acid Sequence↗

Evolution of dnaQ, the gene encoding the editing 3' to 5' exonuclease subunit of DNA polymerase III holoenzyme in Gram-negative bacteria.

The nucleotide sequences of the dnaQ genes from Salmonella typhimurium and Buchnera aphidicola, encoding the epsilon-subunit of the DNA polymerase III holoenzyme, have been determined. The Salmonella typhimurium dnaQ protein consists of 243 amino acid residues with a calculated molecular weight of 27224. The Buchnera aphidicola dnaQ protein contains 233 amino acid residues with a calculated molecular weight of 27170. A multiple sequence alignment of the amino acid sequences of the dnaQ proteins and those of DNA polymerase IIIs from Gram-positive bacteria produced six homologous segments. These homologous segments contain highly conserved amino acid sequence motifs involved in catalytically important metal ion bindings (ligands 1, 2 and 3). However, metal ligand 4 is found to be altered in the 3'-5' exonuclease domain of the family C DNA polymerases and dnaQ proteins in Gram-negative bacteria. From these results, we propose that the last common ancestor of the dnaQ gene of Gram-negative bacteria and the DNA polymerase III gene (pol C gene) of Gram-positive bacteria was a single gene containing both 3'-5' exonuclease and DNA polymerase domains and then the dnaQ gene separated from the polymerase gene in Gram-negative bacteria.

Amino Acid Sequence↗

The arrangement of the transmembrane helices in the secretin receptor family of G-protein-coupled receptors.

The members of the secretin receptor family of G-protein-coupled receptors share no significant sequence similarity to the more familiar rhodopsin-like family. However, multiple sequence alignment analysis reveals seven hydrophobic regions with significant alpha-helical periodicity. Residues that are likely to be buried on the interior of the helical bundle and others that are likely to contact the lipid bilayer are identified. A predicted arrangement of the helical bundle is described in which, by comparison with the arrangement in the rhodopsin family, helices 2 and 7 are more buried within the bundle while helix 3 is more exposed to the lipid bilayer.

Amino Acid Sequence↗

Evidence for a new class of scorpion toxins active against K+ channels.

cDNAs encoding novel long-chain scorpion toxins (64 amino acid residues, including only six cysteines) were isolated from cDNA libraries produced from the venom glands of the scorpions Androctonus australis from Old World and Tityus serrulatus from New World. The encoded peptides were very similar to a recently identified toxin from T. serrulatus, which is active against the voltage-sensitive 'delayed-rectifier' potassium channel, but they were completely different from the long-chain and short-chain scorpion toxins already characterised. However, there was some sequence similarity (42%) between these new toxins, Aa TX Kbeta and Ts TX Kbeta, and scorpion defensins purified from the hemolymph of Buthidae scorpions Leiurus quinquestriatus and A. australis. Thus, according to a multiple sequence alignment using CLUSTAL, these new toxins seem to be related to the scorpion defensins.

Amino Acid Sequence↗

The FxRxHrS motif: a conserved region essential for DNA binding of the VirR response regulator from Clostridium perfringens.

The VirSR two-component signal transduction pathway regulates virulence and toxin production in Clostridium perfringens, the causative agent of gas gangrene. The response regulator, VirR, binds to repeat sequences located upstream of the promoter and is directly responsible for the transcriptional activation of pfoA, the structural gene for the cholesterol-dependent cytolysin, perfringolysin O. Comparative sequence analysis of the 236 amino acid residue VirR protein revealed a two-domain structure: a typical N-terminal response regulator domain and an uncharacterised C-terminal domain. Database searching revealed that over 40 other proteins, many of which appeared to be response regulators or transcriptional activators, had homology with the VirR C-terminal domain (VirRc). Multiple sequence alignment of this VirRc family revealed a highly conserved region that was designated the FxRxHrS motif. By deletion analysis this motif was shown to be essential for the functional integrity of the VirR protein. Alanine scanning mutagenesis and subsequent phenotypic analysis indicated that conserved residues located within the motif were required for activity. These residues extended from L179 to N194. More detailed site-directed mutagenesis showed that amino acid residues R186, H188 and S190 were essential for activity since even conservative substitutions in these positions resulted in non-functional proteins. Three of the mutant proteins, R186K, S190A and S190C, were purified and shown by in vitro gel shift analysis to be unable to bind to the specific target DNA with the same efficiency as the wild-type protein. These data reveal for the first time that VirRc functions as a DNA binding domain in which the highly conserved FxRxHrS motif has a functional role. These studies have important implications for this new family of transcriptional factors since they imply that the conserved FxRxHrS motif may be involved in DNA binding in all of these proteins, irrespective of their biological role.

Amino Acid Motifs↗