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

PubMed2006

An interactive visualization tool to explore the biophysical properties of amino acids and their contribution to substitution matrices.

BACKGROUND: Quantitative descriptions of amino acid similarity, expressed as probabilistic models of evolutionary interchangeability, are central to many mainstream bioinformatic procedures such as sequence alignment, homology searching, and protein structural prediction. Here we present a web-based, user-friendly analysis tool that allows any researcher to quickly and easily visualize relationships between these bioinformatic metrics and to explore their relationships to underlying indices of amino acid molecular descriptors. RESULTS: We demonstrate the three fundamental types of question that our software can address by taking as a specific example the connections between 49 measures of amino acid biophysical properties (e.g., size, charge and hydrophobicity), a generalized model of amino acid substitution (as represented by the PAM74-100 matrix), and the mutational distance that separates amino acids within the standard genetic code (i.e., the number of point mutations required for interconversion during protein evolution). We show that our software allows a user to recapture the insights from several key publications on these topics in just a few minutes. CONCLUSION: Our software facilitates rapid, interactive exploration of three interconnected topics: (i) the multidimensional molecular descriptors of the twenty proteinaceous amino acids, (ii) the correlation of these biophysical measurements with observed patterns of amino acid substitution, and (iii) the causal basis for differences between any two observed patterns of amino acid substitution. This software acts as an intuitive bioinformatic exploration tool that can guide more comprehensive statistical analyses relating to a diverse array of specific research questions.

Amino Acid Sequence↗

PubMed1987

A strategy for the rapid multiple alignment of protein sequences. Confidence levels from tertiary structure comparisons.

An algorithm is presented for the multiple alignment of protein sequences that is both accurate and rapid computationally. The approach is based on the conventional dynamic-programming method of pairwise alignment. Initially, two sequences are aligned, then the third sequence is aligned against the alignment of both sequences one and two. Similarly, the fourth sequence is aligned against one, two and three. This is repeated until all sequences have been aligned. Iteration is then performed to yield a final alignment. The accuracy of sequence alignment is evaluated from alignment of the secondary structures in a family of proteins. For the globins, the multiple alignment was on average 99% accurate compared to 90% for pairwise comparison of sequences. For the alignment of immunoglobulin constant and variable domains, the use of many sequences yielded an alignment of 63% average accuracy compared to 41% average for individual variable/constant alignments. The multiple alignment algorithm yields an assignment of disulphide connectivity in mammalian serotransferrin that is consistent with crystallographic data, whereas pairwise alignments give an alternative assignment.

Algorithms↗

PubMed1999

Identification of a calcium binding site in Staphylococcus hyicus lipase: generation of calcium-independent variants.

In this study we have identified the presence of a high-affinity binding site for calcium in the lipase from Staphylococcus hyicus. By means of isothermal titration calorimetry we showed that the enzyme binds one calcium per molecule of enzyme with a dissociation constant of 55 microM. The residual activity of the apoenzyme compared to the activity in the presence of calcium ions varies from 65% at 10 degreesC to nearly zero at 40 degreesC. On the basis of primary sequence alignment with other staphylococcal lipases and the lipases from Bacillus thermocatenulatus and from Pseudomonas glumae in combination with site-directed mutagenesis, aspartates 354 and 357 could be identified as calcium ligands. Kinetic measurements with the D357E variant showed that replacement of Asp357 by a glutamate decreased the affinity for calcium ions 30-fold. Introduction of a lysine, an asparagine, or an alanine at position 357 and of a lysine or an asparagine at position 354 resulted in calcium-independent variants. Isothermal titration calorimetry confirmed the loss of calcium binding. Although the D357K, D357N, and D357A variants did not bind calcium, at room temperature they were nearly as active as wild-type lipase in the presence of calcium, but at elevated temperatures these calcium-independent lipases showed a reduced activity. Over the whole temperature range the activities of the D354K and D354N variants are significantly lower than wild-type enzyme in the presence of calcium and are comparable to the activity of the wild-type apoenzyme. Our results show that binding of calcium is important for the structural stabilization of staphylococcal lipases (and possibly other lipases) and that it is possible to engineer calcium-independent variants on the basis of limited structural homology with another lipase.

Amino Acid Sequence↗

PubMed2002

Cloning, overexpression, and characterization of a bacterial Ca2+-dependent phospholipase D.

Phospholipase D (PLD), an important enzyme involved in signal transduction in mammals, is also secreted by many microorganisms. A highly conserved HKD motif has been identified in most PLD homologs in the PLD superfamily. However, the Ca(2+)-dependent PLD from Streptomyces chromofuscus exhibits little homology to other PLDs. We have cloned (using DNA isolated from the ATCC type strain), overexpressed in Escherichia coli (two expression systems, pET-23a(+) and pTYB11), and purified the S. chromofuscus PLD. Based on attempts at sequence alignment with other known Ca(2+)-independent PLD enzymes from Streptomyces species, we mutated five histidine residues (His72, His171, His187, His200, His226) that could be part of variants of an HKD motif. Only H187A and H200A showed dramatically reduced activity. However, mutation of these histidine residues to alanine also significantly altered the secondary structure of PLD. Asparagine replacements at these positions yielded enzymes with structure and activity similar to the recombinant wild-type PLD. The extent of phosphatidic acid (PA) activation of PC hydrolysis by the recombinant PLD enzymes differed in magnitude from PLD purified from S. chromofuscus culture medium (a 2-fold activation rather than 4-5-fold). One of the His mutants, H226A, showed a 12-fold enhancement by PA, suggesting this residue is involved in the kinetic activation. Another notable difference of this bacterial PLD from others is that it has a single cysteine (Cys123); other Streptomyces Ca(2+)-independent PLDs have eight Cys involved in intramolecular disulfide bonds. Both C123A and C123S, with secondary structure and stability similar to recombinant wild-type PLD, exhibited specific activity reduced by 10(-5) and 10(-4). The Cys mutants still bound Ca(2+), so that it is likely that this residue is part of the active site of the Ca(2+)-dependent PLD. This would suggest that S. chromofuscus PLD is a member of a new class of PLD enzymes.

Amino Acid Sequence↗

PubMed2006

Dirigent proteins in conifer defense: gene discovery, phylogeny, and differential wound- and insect-induced expression of a family of DIR and DIR-like genes in spruce (Picea spp.).

The outer stem tissues of conifers provide a durable constitutive and inducible defense barrier consisting of suberized or lignified periderm, sclereids, a network of terpenoid-filled resin ducts, and phenolic phloem parenchyma cells. Microarray gene expression profiling of Sitka spruce (Picea sitchensis) bark attacked by stem-boring weevils (Pissodes strobi) or through mechanical wounding demonstrated significant accumulation of transcripts resembling dirigent protein (DIR) genes. To investigate this gene family and its spatial and temporal patterns of expression in conifer defense, we isolated cDNAs representing 19 unique DIR and DIR-like genes from Sitka spruce, white spruce (P. glauca), and interior spruce (P. glauca x engelmannii). Sequence alignments also identified a large number of DIR-like proteins in other plant species, which share several conserved protein motifs with known DIR proteins. Phylogenetic analysis of 72 DIR and DIR-like proteins suggests five distinct subfamilies, DIR-a and four DIR-like subfamilies (DIR-b, DIR-c, DIR-d and DIR-e). Previously characterized members of the DIR-a subfamily direct stereoselective phenolic coupling reactions in the formation of lignans and possibly lignins. The spruce genes identified here are members of the DIR-a and DIR-b subfamilies. Using gene-specific quantitative real-time PCR we measured constitutive expression for six DIR-a genes and three DIR-like genes in different stem tissues, green shoot tips, and roots of Sitka spruce. DIR-like genes revealed ubiquitous high expression in all tissues. In contrast, the six DIR-a genes showed a gradient of transcript abundance in stem tissues with highest levels in the outer cortex and lowest levels in the inner xylem. Gene-specific transcript profiling of six DIR-a genes confirmed rapid and strong accumulation (up to 500-fold) in wound- and weevil-induced stem bark and xylem. These findings suggest a role for spruce DIR genes in constitutive and induced phenolic defense mechanisms against stem-boring insects.

Amino Acid Sequence↗

PubMed2003

Analysis of the Sarcocystis neurona microneme protein SnMIC10: protein characteristics and expression during intracellular development.

Sarcocystis neurona, an apicomplexan parasite, is the primary causative agent of equine protozoal myeloencephalitis. Like other members of the Apicomplexa, S. neurona zoites possess secretory organelles that contain proteins necessary for host cell invasion and intracellular survival. From a collection of S. neurona expressed sequence tags, we identified a sequence encoding a putative microneme protein based on similarity to Toxoplasma gondii MIC10 (TgMIC10). Pairwise sequence alignments of SnMIC10 to TgMIC10 and NcMIC10 from Neospora caninum revealed approximately 33% identity to both orthologues. The open reading frame of the S. neurona gene encodes a 255 amino acid protein with a predicted 39-residue signal peptide. Like TgMIC10 and NcMIC10, SnMIC10 is predicted to be hydrophilic, highly alpha-helical in structure, and devoid of identifiable adhesive domains. Antibodies raised against recombinant SnMIC10 recognised a protein band with an apparent molecular weight of 24 kDa in Western blots of S. neurona merozoites, consistent with the size predicted for SnMIC10. In vitro secretion assays demonstrated that this protein is secreted by extracellular merozoites in a temperature-dependent manner. Indirect immunofluorescence analysis of SnMIC10 showed a polar labelling pattern, which is consistent with the apical position of the micronemes, and immunoelectron microscopy provided definitive localisation of the protein to these secretory organelles. Further analysis of SnMIC10 in intracellular parasites revealed that expression of this protein is temporally regulated during endopolygeny, supporting the view that micronemes are only needed during host cell invasion. Collectively, the data indicate that SnMIC10 is a microneme protein that is part of the excreted/secreted antigen fraction of S. neurona. Identification and characterisation of additional S. neurona microneme antigens and comparisons to orthologues in other Apicomplexa could provide further insight into the functions that these proteins serve during invasion of host cells.

Amino Acid Sequence↗

PubMed1997

Functional analysis of point mutations in human flap endonuclease-1 active site.

Human flap endonuclease-1 (hFEN-1) is highly homologous to human XPG, Saccharomyces cerevisiae RAD2 and S.cerevisiae RTH1 and shares structural and functional similarity with viral exonucleases such as T4 RNase H, T5 exonuclease and prokaryotic DNA polymerase 5'nucleases. Sequence alignment of 18 structure-specific nucleases revealed two conserved nuclease domains with seven conserved carboxyl residues and one positively charged residue. In a previous report, we showed that removal of the side chain of each individual acidic residue results in complete loss of flap endonuclease activity. Here we report a detailed analysis of substrate cleavage and binding of these mutant enzymes as well as of an additional site-directed mutation of a conserved acidic residue (E160). We found that the active mutant (R103A) has substrate binding and cleavage activity indistinguishable from the wild type enzyme. Of the inactive mutants, one (D181A) has substrate binding properties comparable to the wild type, while three others (D34A, D86A and E160A) bind with lower apparent affinity (2-, 9- and 18-fold reduced, respectively). The other mutants (D158A, D179A and D233A) have no detectable binding activity. We interpret the structural implications of these findings using the crystal structures of related enzymes with the flap endonuclease activity and propose that there are two metal ions (Mg2+or Mn2+) in hFEN enzyme. These two metal coordinated active sites are distinguishable but interrelated. One metal site is directly involved in nucleophile attack to the substrate phosphodiester bonds while the other may stabilize the structure for the DNA substrate binding. These two sites may be relatively close since some of carboxyl residues can serve as ligands for both sites.

Amino Acid Sequence↗

PubMed1998

Biochemical role of the Cryptococcus neoformans ADE2 protein in fungal de novo purine biosynthesis.

Comparative studies of 5-aminoimidazole ribonucleotide (AIR) carboxylases from Escherichia coli and Gallus gallus have identified this central step in de novo purine biosynthesis as a case for unusual divergence in primary metabolism. Recent discoveries establish the fungal AIR carboxylase, encoded by the ADE2 gene, as essential for virulence in certain pathogenic organisms. This investigation is a biochemical analysis that links the fungal ADE2 protein to the function of the E. coli AIR carboxylase system. A cDNA clone of ADE2 from Cryptococcus neoformans was isolated by genetic complementation of a purE-deficient strain of E. coli. High-level expression of the C. neoformans ADE2 was achieved, which enabled the production and purification of AIR carboxylase. Amino acid sequence alignments, C-terminal deletion mutants, and biochemical assays indicate that the ADE2 enzyme is a two-domain, bifunctional protein. The N-terminal domain is related to E. coli PurK and a series of kinetic experiments show that the ADE2-PurK activity uses AIR, ATP, and HCO3- as substrates. The biosynthetic product of the ADE2-PurK reaction was identified as N5-carboxyaminoimidazole ribonucleotide (N5-CAIR) by 1H NMR, thus confirming that the C-terminal domain contains a catalytic activity similar to that of the E. coli PurE. By using an in situ system for substrate production, the steady-state kinetic constants for turnover of N5-CAIR by ADE2 were determined and together with stoichiometry measurements, these data indicate that ADE2 has a balance in the respective catalytic turnovers to ensure efficient flux. Distinctive features of the PurE active site were probed using 4-nitro-5-aminoimidazole ribonucleotide (NAIR), an analog of the product 4-carboxy-5-aminoimidazole ribonucleotide (CAIR). NAIR was shown to be a selective inhibitor of the ADE2-PurE activity (K1 = 2.4 microM), whereas it is a slow-binding inhibitor of the G. gallus enzyme which further distinguishes the fungal ADE2 from the G. gallus AIR carboxylase. As such, this enzyme represents a novel intracellular target for the discovery of antifungal agents.

Adenosine Triphosphate↗

PubMed2003

A 3D model of SARS_CoV 3CL proteinase and its inhibitors design by virtual screening.

AIM: To constructed a three-dimensional (3D) model for the 3C like (3CL) proteinase of SARS coronavirus (SARS-CoV), and to design inhibitors of the 3CL proteinase based on the 3D model. METHODS: Bioinformatics analyses were performed to search the homologous proteins of the SARS-CoV 3CL proteinase from the GenBank and PDB database. A 3D model of the proteinase was constructed by using homology modeling technique. Targeting to the 3D model and its X-ray crystal structure of the main proteinase (Mpro) of transmissible gastroenteritis virus (TGEV), virtual screening was performed employing molecular docking method to identify possible 3CL proteinase inhibitors from small molecular databases. RESULTS: Sequence alignment indicated that the SARS-CoV 3CL proteinase was extremely homologous to TGEV Mpro, especially the substrate-binding pocket (active site). Accordingly, a 3D model for the SARS-CoV 3CL proteinase was constructed based on the crystal structure of TGEV Mpro. The 3D model adopts a similar fold of the TGEV Mpro, its structure and binding pocket feature are almost as same as that of TGEV Mpro. The tested virtual screening indicated that 73 available proteinase inhibitors in the MDDR database might dock into both the binding pockets of the TGEV Mpro and the SARS-CoV 3CL proteinase. CONCLUSIONS: Either the 3D model of the SARS-CoV 3CL proteinase or the X-ray crystal structure of the TGEV Mpro may be used as a starting point for design anti-SARS drugs. Screening the known proteinase inhibitors may be an appreciated shortcut to discover anti-SARS drugs.

Computational Biology↗

PubMed2003

Assessing functional divergence in EF-1alpha and its paralogs in eukaryotes and archaebacteria.

A number of methods have recently been published that use phylogenetic information extracted from large multiple sequence alignments to detect sites that have changed properties in related protein families. In this study we use such methods to assess functional divergence between eukaryotic EF-1alpha (eEF-1alpha), archaebacterial EF-1alpha (aEF-1alpha) and two eukaryote-specific EF-1alpha paralogs-eukaryotic release factor 3 (eRF3) and Hsp70 subfamily B suppressor 1 (HBS1). Overall, the evolutionary modes of aEF-1alpha, HBS1 and eRF3 appear to significantly differ from that of eEF-1alpha. However, functionally divergent (FD) sites detected between aEF-1alpha and eEF-1alpha only weakly overlap with sites implicated as putative EF-1beta or aminoacyl-tRNA (aa-tRNA) binding residues in EF-1alpha, as expected based on the shared ancestral primary translational functions of these two orthologs. In contrast, FD sites detected between eEF-1alpha and its paralogs significantly overlap with the putative EF-1beta and/or aa-tRNA binding sites in EF-1alpha. In eRF3 and HBS1, these sites appear to be released from functional constraints, indicating that they bind neither eEF-1beta nor aa-tRNA. These results are consistent with experimental observations that eRF3 does not bind to aa-tRNA, but do not support the 'EF-1alpha-like' function recently proposed for HBS1. We re-assess the available genetic data for HBS1 in light of our analyses, and propose that this protein may function in stop codon-independent peptide release.

Amino Acid Sequence↗

PubMed1998

A glutamate residue contributes to the exopeptidase specificity in aminopeptidase A.

Aminopeptidase A (EC 3.4.11.7, APA) is a 130 kDa membrane-bound aminopeptidase that contains the consensus sequence HEXXH (385-389) found in the zinc metalloprotease family, the zincins. Sequence alignment of the mouse APA with other monozinc-aminopeptidases indicates the presence of a highly conserved glutamate residue (Glu352 in APA) found in the conserved motif GAMEN (349-353). In monozinc-aminopeptidases, the negative charge of the glutamate side-chain carboxylate may constitute the anionic binding site involved in the recognition of the free amino group of substrates or inhibitors. The functional role of Glu352 in APA was investigated by substituting this residue with an aspartate (Asp352), a glycine (Gly352), a glutamine (Gln352) or an arginine (Arg352) residue by site-directed mutagenesis. Kinetic studies showed that the Km values of the mutant enzymes were unaffected, whereas kcat values were decreased 10-250-fold, resulting in a 10-, 30- 260- and 400-fold reduction in cleavage efficiencies for the mutants Asp352, Gly352, Gln352 and Arg352 respectively. The inhibitory potency of two different classes of inhibitors, a thiol and a phosphonate compound, was significantly (P<0.05) decreased by 10- and 4-fold respectively in the mutated enzymes. Moreover, the inhibitory potency of angiotensin I, used as a competitor of the synthetic substrate alpha-l-glutamyl beta-naphthylamide, displayed a 4-fold reduction (P<0.01) in the mutated enzymes, whereas the Ki values of its N-acetyl derivative were unchanged. These data strongly suggest that Glu352 is involved in the catalytic process of APA and contributes to the exopeptidase activity of this enzyme through interaction with the N-terminal part of substrates or inhibitors.

Amino Acid Sequence↗

PubMed2005

Structure-activity relationship of an alpha-toxin Bs-Tx28 from scorpion (Buthus sindicus) venom suggests a new alpha-toxin subfamily.

Scorpion venoms are among the most widely known source of peptidyl neurotoxins used for callipering different ion channels, e.g., for Na(+), K(+), Ca(+) or Cl(-). An alpha-toxin (Bs-Tx28) has been purified from the venom of scorpion Buthus sindicus, a common yellow scorpion of Sindh, Pakistan. The primary structure of Bs-Tx28 was established using a combination of MALDI-TOF-MS, LC-ESI-MS, and automated Edman degradation analysis. Bs-Tx28 consists of 65 amino acid residues (7274.3+/-2Da), including eight cysteine residues, and shows very high sequence identity (82-94%) with other long-chain alpha-neurotoxins, active against receptor site-3 of mammalian (e.g., Lqq-IV and Lqh-IV from scorpions Leiurus sp.) and insect (e.g., BJalpha-IT and Od-1 from Buthotus judaicus and Odonthobuthus doriae, respectively) voltage-gated Na(+) channels. Multiple sequence alignment and phylogenetic analysis of Bs-Tx28 with other known alpha- and alpha-like toxins suggests the presence of a new and separate subfamily of scorpion alpha-toxins. Bs-Tx28 which is weakly active in both, mammals and insects (LD(50) 0.088 and 14.3microg/g, respectively), shows strong induction of the rat afferent nerve discharge in a dose-dependent fashion (EC(50)=0.01microg/mL) which was completely abolished in the presence of tetrodotoxin suggesting the binding of Bs-Tx28 to the TTX-sensitive Na(+)-channel. Three-dimensional structural features of Bs-Tx28, established by homology modeling, were compared with other known classical alpha-mammal (AaH-II), alpha-insect (Lqh-alphaIT), and alpha-like (BmK-M4) toxins and revealed subtle variations in the Nt-, Core-, and RT-CT-domains (functional domains) which constitute a "necklace-like" structure differing significantly in all alpha-toxin subfamilies. On the other hand, a high level of conservation has been observed in the conserved hydrophobic surface with the only substitution of W43 (Y43/42) and an additional hydrophobic character at position F40 (L40/A/V/G39), as compared to the other mentioned alpha-toxins. Despite major differences within the primary structure and activities of Bs-Tx28, it shares a common structural and functional motif (e.g., transRT-farCT) within the RT-CT domain which is characteristic of scorpion alpha-mammal toxins.

Action Potentials↗

PubMed1993

Mutation of a highly conserved acidic residue present in the second intracellular loop of G-protein-coupled receptors does not impair hormone binding or signal transduction of the luteinizing hormone/chorionic gonadotropin receptor.

Sequence alignment shows that there is a highly conserved acidic residue (D or E) at the boundary between the third transmembrane domain and the second intracellular loop of the superfamily of G-protein-coupled receptors. Previous mutagenesis studies demonstrated that substitution of this acidic residue in the beta 2-adrenergic, muscarinic m1, and alpha 2A-adrenergic receptors by the corresponding amide preserved high affinity agonist binding, but significantly reduced or completely abolished activation of the respective effector. To determine whether the corresponding amino acid residue (E441) played a similar role in the functions of the rat LH/CG receptor, we used site-directed mutagenesis to substitute it by D or Q. The wild-type and mutant receptors (E441D or E441Q) were then transfected into human embryonic kidney 293 cells and tested for their ability to bind hCG and respond to it with increased cAMP accumulation. As predicted, the mutant LH/CG receptors were found to bind hCG with high affinity. In contrast to the results summarized above, however, an E441Q or an E441D mutation in the LH/CG receptor results in only a slight increase in the EC50 for cAMP accumulation without decreasing the maximal response attained. The most remarkable effect of these mutations was on localization of the receptor. Thus, while most of the receptors expressed in cells transfected with the E441D mutant could be detected by measuring hormone binding to intact cells, most of the receptors expressed in cells transfected with the E441Q mutant could be detected only upon solubilization of the cells with detergent.

Amino Acid Sequence↗

PubMed1996

Crystal structure of an acidic neurotoxin from scorpion Buthus martensii Karsch at 1.85 A resolution.

The crystal structure of an acidic scorpion neurotoxin, BmK M8, purified from Chinese scorpion Buthus martensii Karsch (BmK), has been determined by the molecular replacement method. It is the first structure of an acidic alpha-scorpion neurotoxin reported so far. The crystals adopt a symmetry of space group P2(1) and contain one molecule per asymmetric unit. The structure has been refined to an R factor of 18.1% using reflection data in the range of 8 to 1.85 A resolution, with standard deviations from ideal geometry of 0.017 A and 2.43 degrees for bond length and angle, respectively. The 12 residues at the C terminus with unknown sequence were determined by crystallographic refinement. The refined model shows that the structural core, consisting of a motif beta alpha beta beta, is similar to that of toxin II from Androctonus australis Hector (AaH II) or Variant 3 from Centruroides sculpturatus Ewing (CsE V3). The three conformationally variable loops protruding from this structural core are different from that of AaH II, and especially from that of CsE V3. Compared with the most potent and basic alpha-toxin AaH II, the BmK M8 is a relatively inactive toxin (1100 times less active than AaH II) with an unusually low isoelectric point (pI 5.3). Sequence alignment of the two toxins shows a difference of 26 residues (40.6%). Among them four basic or neutral residues in AaH II, namely Val10, Lys28, Val55 and Gly59, are changed to acidic glutamate in BmK M8. The residues Glu10, Glu28 and Glu55 of BmK M8 are located on a surface (Face B), opposite the "conserved hydrophobic surface" (Face A). The latter is a functionally important area proposed by Fontecilla-Camps et al. Our observations suggest that in addition to Face A, Face B may also be involved in the biological activity of scorpion toxins. The structure of BmK M8 shows an evident conformational change of the alpha-amino group at the N terminus and a deorganization of Arg2 caused by the mutation D53A. These structural changes may also be responsible for the weak toxicity of BmK M8. In association with the information from chemical modifications, a multisite binding mode for toxin-receptor interaction and three "toxic regions" in relevance to the binding process, including Face A, Face B and Site C, are proposed. Face A, mainly consisting of Tyr5, 35, 47, the alpha-amino group, Arg2 and Asp3, may be more essential for the binding. Face B, mainly comprising conserved residues Tyr14, 21, Lys28 and Val55, may contribute to the high efficacy of the binding process and substitutions by acidic residues in this area could strongly weaken the toxic activity. Site C, formed by Lys58 and Arg62 at the C terminus and Arg41 and Tyr42 from loop 38-44, may be involved in binding site specificity.

Amino Acid Sequence↗

PubMed1990

Active site labeling of the shikimate pathway enzyme, dehydroquinase. Evidence for a common substrate binding site within dehydroquinase and dehydroquinate synthase.

Dehydroquinase, the third enzyme of the shikimate biosynthetic pathway, is inactivated by iodoacetate. Iodoacetate behaves as an affinity label for the Escherichia coli enzyme with a Ki of 30 mM and a limiting inactivation rate of 0.014 min-1 at pH 7.0 and 25 degrees C. Affinity labeling is mediated by the negative charge of the reagent since iodoacetamide does not inactivate the enzyme. 2.1-2.3 mol of carboxymethyl groups are incorporated per mol of protein monomer resulting in 90% inactivation of enzymic activity. The majority of the bound label (80%) is split equally between 2 methionine residues, Met-23 and Met-205, which were identified by sequencing radiolabelled peptide fragments isolated after proteolytic digestion. An equilibrium mixture of the substrate (dehydroquinate) and product (dehydroshikimate) substantially reduces the inactivation rate and specifically decreases the incorporation of label at both of these site, implicating them as being in or near the active site of the enzyme. Sequence alignments with other biosynthetic dehydroquinases show that of the 2 methionine residues only Met-205 is conserved. N-terminal alignments of all the available dehydroquinase sequences (both catabolic and biosynthetic classes) revealed that Met-23, although itself not conserved, resides within a cluster of conserved sequence which may constitute part of the dehydroquinate binding site. A consensus sequence was derived from these alignments and used to probe the protein sequence data banks. A related sequence was found in dehydroquinate synthase, the enzyme which precedes dehydroquinase in the shikimate pathway. These results suggest that we have identified part of the dehydroquinate binding site in both enzymes.

Affinity Labels↗

PubMed1997

The crystal structure of the flavin containing enzyme dihydroorotate dehydrogenase A from Lactococcus lactis.

BACKGROUND: . Dihydroorotate dehydrogenase (DHOD) is a flavin mononucleotide containing enzyme, which catalyzes the oxidation of (S)-dihydroorotate to orotate, the fourth step in the de novo biosynthesis of pyrimidine nucleotides. Lactococcus lactis contains two genes encoding different functional DHODs whose sequences are only 30% identical. One of these enzymes, DHODA, is a highly efficient dimer, while the other, DHODB, shows optimal activity only in the presence of an iron-sulphur cluster containing protein with which it forms a complex tetramer. Sequence alignments have identified three different families among the DHODs: the two L. lactis enzymes belong to two of the families, whereas the enzyme from E. coli is a representative of the third. As no three-dimensional structures of DHODs are currently available, we set out to determine the crystal structure of DHODA from L. lactis. The differences between the two L. lactis enzymes make them particularly interesting for studying flavoprotein redox reactions and for identifying the differences between the enzyme families. RESULTS: . The crystal structure of DHODA has been determined to 2.0 resolution. The enzyme is a dimer of two crystallographically independent molecules related by a non-crystallographic twofold axis. The protein folds into and alpha/beta barrel with the flavin molecule sitting between the top of the barrel and a subdomain formed by several barrel inserts. Above the flavin isoalloxazine ring there is a small water filled cavity, completely buried beneath the protein surface and surrounded by many conserved residues. This cavity is proposed as the substrate-binding site. CONCLUSIONS: . The crystal structure has allowed the function of many of the conserved residues in DHODs to be identified: many of these are associated with binding the flavin group. Important differences were identified in some of the active-site residues which vary across the distinct DHOD families, implying significant mechanistic differences. The substrate cavity, although buried, is located beneath a highly conserved loop which is much less ordered than the rest of the protein and may be important in giving access to the cavity. The location of the conserved residues surrounding this cavity suggests the potential orientation of the substrate.

Amino Acid Sequence↗

PubMed1995

Involvement of Ser-451 and Ser-452 in the catalysis of human gamma-glutamyl transpeptidase.

The serine residue required for catalysis of gamma-glutamyl transpeptidase was identified by site-specific mutagenesis of the conserved serine residues on the basis of sequence alignment of the light subunit of human, rat, pig and two bacterial enzymes. Recombinant human gamma-glutamyl transpeptidases with replacements of these serine residues by Ala were expressed using a baculovirus-insect cell system. Substitutions of Ala at Ser-385, -413 or -425 yielded almost fully active enzymes. However, substitutions of Ala at Ser-451 or -452 yielded enzymes that were only about 1% as active as the wild-type enzyme. Further, their double mutant is only 0.002% as active as the wild type. Kinetic analysis of transpeptidation using glycylglycine as acceptor indicates that the Vmax values of Ser-451 and -452 mutants are substantially decreased (to about 3% of the wild type); however, their Km values for L-gamma-glutamyl-p-nitroanilide as donor were only increased about 5 fold compared to that of the wild type. The double mutation of Ser-451 and -452 further decreased the Vmax value to only about 0.005% of the wild type, while this mutation produced only a minor effect (2-fold increase) on the Km value for the donor. The kinetic values for the hydrolysis reaction of L-gamma-glutamyl-p-nitroanilide in the mutants followed similar trends to those for transpeptidation. The rates of inactivation of Ser-451, -452 and their double mutant enzymes by acivicin, a potent inhibitor, were less than 1% that of the wild-type enzyme. The Ki value of the double mutant for L-serine as a competitive inhibitor of the gamma-glutamyl group is only 9 fold increased over that of the wild type, whereas the Ki for the serine-borate complex, which acts as an inhibitory transition-state analog, was more than 1,000 times higher than for the wild-type enzyme. These results suggest that both Ser-451 and -452 are located at the position able to interact with the gamma-glutamyl group and participate in catalysis, probably as nucleophiles or through stabilization of the transition state.

Amino Acid Sequence↗

PubMed2001

Influence of P450 3A4 SRS-2 residues on cooperativity and/or regioselectivity of aflatoxin B(1) oxidation.

The major human liver drug-metabolizing cytochrome P450 enzymes P450 3A4 and P450 3A5 share >85% amino acid sequence identity yet exhibit different regioselectivity toward aflatoxin B(1) (AFB(1)) biotransformation [Gillam et al. (1995) Arch. Biochem. Biophys. 317, 74-384]. P450 3A4 prefers AFB1 3alpha-hydroxylation, which detoxifies and subsequently eliminates the hepatotoxin, over AFB1 exo-8,9-oxidation. P450 3A5, on the other hand, is a relatively sluggish 3alpha-hydroxylase and converts AFB(1) predominantly to the genotoxic exo-8,9-epoxide. Using a combination of approaches (sequence alignment, homology modeling and site-directed mutagenesis), we have previously identified several divergent residues in four of the six putative substrate recognition sites (SRSs) of P450 3A4, which when replaced individually with the corresponding amino acid of P450 3A5, resulted in a significant switch of the characteristic P450 3A4 AFB(1) regioselectivity toward that of P450 3A5 [Wang et al. (1998) Biochemistry 37, 12536-12545]. In particular, residues N206 and L210 in SRS-2 were found to be critical for AFB(1) detoxification via 3alpha-hydroxylation, and the corresponding mutants N206S and L210F most closely mimicked P450 3A5, not only in its regioselectivity of AFB(1) metabolism but also in its overall functional capacity. We have now further explored the plausible reasons for such relative inactivity of the SRS-2 mutants by examining N206S and additional mutants (L210A, L211F, L211A, and N206E) and found that the dramatically lowered activities of the N206S mutant are accompanied by a loss of cooperativity of AFB(1) oxidation. Molecular dynamics analyses with an existing P450 3A4 homology model [Szklarz and Halpert (1997) J. Comput. Aided Mol. Des. 11, 265] suggested that N206 (helix F) interacts with E244 (helix G), creating a salt bridge that stabilizes the protein structure and/or defines the active site cavity. To examine this possibility, several E244 mutants (E244A, V, N, S) were tested, of which E244S was the most notable for its relatively greater impairment of P450 3A4-dependent AFB(1) 3alpha-hydroxylation. However, the results with these E244 mutants failed to validate the N206-E244 interaction predicted from these molecular dynamics analyses. Collectively, our findings to date have led us to reconsider our original interpretations and to reexamine them in the light of AFB(1) molecular modeling analyses with a newly refined P450 3A4 homology model. These analyses predicted that F304 in SRS-4 (I-helix) plays a pivotal role in AFB(1) binding at the active site in either orientation leading to 3alpha- or exo-8,9-oxidation. Consistent with this prediction, conversion of F304 to Ala abolished P450 3A4-dependent AFB(1) 3alpha-hydroxylation and exo-8,9-oxidation.

Aflatoxin B1↗