Dichloromethane mediated in vivo selection and functional characterization of rat glutathione S-transferase theta 1-1 variants.

Methylobacterium dichloromethanicum DM4 is able to grow with dichloromethane as the sole carbon and energy source by using a dichloromethane dehalogenase/glutathione S-transferase (GST) for the conversion of dichloromethane to formaldehyde. Mammalian homologs of this bacterial enzyme are also known to catalyze this reaction. However, the dehalogenation of dichloromethane by GST T1-1 from rat was highly mutagenic and toxic to methylotrophic bacteria. Plasmid-driven expression of rat GST T1-1 in strain DM4-2cr, a mutant of strain DM4 lacking dichloromethane dehalogenase, reduced cell viability 10(5)-fold in the presence of dichloromethane. This effect was exploited to select dichloromethane-resistant transconjugants of strain DM4-2cr carrying a plasmid-encoded rGSTT1 gene. Transconjugants that still expressed the GST T1 protein after dichloromethane treatment included rGSTT1 mutants encoding protein variants with sequence changes from the wild-type ranging from single residue exchanges to large insertions and deletions. A structural model of rat GST T1-1 suggested that sequence variation was clustered around the glutathione activation site and at the protein C-terminus believed to cap the active site. The enzymatic activity of purified His-tagged GST T1-1 variants expressed in Escherichia coli was markedly reduced with both dichloromethane and the alternative substrate 1,2-epoxy-3-(4'-nitrophenoxy)propane. These results provide the first experimental evidence for the involvement of Gln102 and Arg107 in catalysis, and illustrate the potential of in vivo approaches to identify catalytic residues in GSTs whose activity leads to toxic effects.

The peopling of New Guinea: evidence from class I human leukocyte antigen.

This study utilizes newly developed direct DNA typing methods for human leukocyte antigens (HLA) to provide new information about the peopling of New Guinea. The complete polymorphism of eight Melanesian populations was examined. The groups included were highlanders, northern and southern highlands fringe populations, a Sepik population, northern and southern coastal New Guinea populations, and populations from the Bismarck Archipelago and New Caledonia. The study concluded that, based on HLA and other evidence. Melanesians are likely to have evolved largely from the same ancestral stock as Aboriginal Australians but to have since differentiated. Highlanders are likely to be descendants of earlier migrations who have been isolated for a long period of time. Northern highlands fringe and Sepik populations are likely to share a closer common ancestry but to have differentiated due to long term isolation and the relative proximity to the coast of the Sepik. Southern fringe populations are likely to have a different origin, possibly from the Gulf region, although there may be some admixture with neighboring groups. Coastal populations have a wider range of polymorphisms because of the genetic trail left by later population movement along the coast from Asia that did not reach Australia or remote Oceania. Other polymorphisms found in these populations may have been introduced by the movement of Austronesian-speaking and other more recent groups of people into the Pacific, because they share many polymorphisms with contemporary southeast Asians, Polynesians, and Micronesians that are not found in highlanders or Aboriginal Australians. There is evidence suggestive of later migration to Melanesia from Polynesia and Micronesia.

Identification of novel glutathione transferases and polymorphic variants by expressed sequence tag database analysis.

The human expressed sequence tag (EST) database can be searched by different sequence alignment strategies to identify new members of gene families and allelic variants. To illustrate the value of database analysis for gene discovery, we have focused on the glutathione S-transferase (GST) super family, an approach that has led to the identification of the Zeta class. The Zeta class GSTs catalyze the glutathione-dependent biotransformation of alpha-haloacids and the isomerization of maleylacetoacetic acid to fumarylacetoacetic acid, an essential step in the catabolism of tyrosine. Allelic variants of the GST Z1 and GST A2 genes have also been identified by EST database analysis. One GST Z1 variant (GST Z1A) has significantly higher activity with dichloroacetic acid as a substrate than other GST Z1 isoforms. This variant may be important in the clinical treatment of lactic acidosis where dichloroacetic acid is prescribed. Our experience with the application of EST database searching methods suggests that it may be productively applied to other gene families of pharmacogenetic interest.

The glutathione transferase structural family includes a nuclear chloride channel and a ryanodine receptor calcium release channel modulator.

The ubiquitous glutathione transferases (GSTs) catalyze glutathione conjugation to many compounds and have other diverse functions that continue to be discovered. We noticed sequence similarities between Omega class GSTs and a nuclear chloride channel, NCC27 (CLIC1), and show here that NCC27 belongs to the GST structural family. The structural homology prompted us to investigate whether the human Omega class glutathione transferase GSTO1-1 forms or modulates ion channels. We find that GSTO1-1 modulates ryanodine receptors (RyR), which are calcium channels in the endoplasmic reticulum of various cells. Cardiac RyR2 activity was inhibited by GSTO1-1, whereas skeletal muscle RyR1 activity was potentiated. An enzymatically active conformation of GSTO1-1 was required for inhibition of RyR2, and mutation of the active site cysteine (Cys-32 --> Ala) abolished the inhibitory activity. We propose a novel role for GSTO1-1 in protecting cells containing RyR2 from apoptosis induced by Ca(2+) mobilization from intracellular stores.

Human glutathione transferase T2-2 discloses some evolutionary strategies for optimization of the catalytic activity of glutathione transferases.

Steady state, pre-steady state kinetic experiments, and site-directed mutagenesis have been used to dissect the catalytic mechanism of human glutathione transferase T2-2 with 1-menaphthyl sulfate as co-substrate. This enzyme is close to the ancestral precursor of the more recently evolved glutathione transferases belonging to Alpha, Pi, and Mu classes. The enzyme displays a random kinetic mechanism with very low k(cat) and k(cat)/K(m)((GSH)) values and with a rate-limiting step identified as the product release. The chemical step, which is fast and causes product accumulation before the steady state catalysis, strictly depends on the deprotonation of the bound GSH. Replacement of Arg-107 with Ala dramatically affects the fast phase, indicating that this residue is crucial both in the activation and orientation of GSH in the ternary complex. All pre-steady state and steady state kinetic data were convincingly fit to a kinetic mechanism that reflects a quite primordial catalytic efficiency of this enzyme. It involves two slowly interconverting or not interconverting enzyme populations (or active sites of the dimeric enzyme) both able to bind and activate GSH and strongly inhibited by the product. Only one population or subunit is catalytically competent. The proposed mechanism accounts for the apparent half-site behavior of this enzyme and for the apparent negative cooperativity observed under steady state conditions. These findings also suggest some evolutionary strategies in the glutathione transferase family that have been adopted for the optimization of the catalytic activity, which are mainly based on an increased flexibility of critical protein segments and on an optimal orientation of the substrate.

Identification of a new mutation (Gly420Ser), distal to the active site, that leads to factor XIII deficiency.

The molecular defects of the factor XIII A subunit gene were studied in a patient with factor XIII deficiency. Mutation analysis was performed on amplified DNA from each exon of this gene by single-strand conformation polymorphism (SSCP) and DNA sequencing techniques. A substitution of guanine by adenine at nucleotide 1258 in exon 10 of the coagulation factor XIII A subunit gene has been identified in the patient. The mutation results in the replacement of Gly420 by Ser in the core domain of the enzyme. Restriction enzyme analysis of amplified exon 10 DNA confirmed that the patient was homozygous for this mutation. A family study revealed that the mutation was inherited from both parents, who were first cousins. The potential effects of the mutation were predicted by molecular modeling of the amino acid substitution within the coordinates of the crystal structure. The substitution occurred within the core domain of the enzyme at a residue completely conserved among all known members of the transglutaminase family. The model of the mutant protein suggests that although the substitution of Gly420 by Ser causes only minor readjustment of the residues and does not appear to be particularly deleterious in terms of structure, the mutation is, however, likely to decrease the molecule's ability to undergo the conformational change that is thought to be required for full transglutaminase activity. Our data strongly support the previously published information about the functional significance of the residues surrounding, but not forming, the catalytic pocket in the A subunit of factor XIII.

Database analysis and gene discovery in pharmacogenetics.

The global genome research effort has resulted in the creation of extensive DNA and protein sequence databases that are a valuable resource for the identification of new genes and polymorphic variants of enzymes of pharmacogenetic interest. Previously undescribed members of gene families with novel functions and substrate specificities can be identified by database searching and sequence alignment strategies. Since the expressed sequence tag (EST) database contains sequences from many individuals, it can be searched for evidence of polymorphisms that can significantly influence enzyme function. The different approaches to these forms of analysis are reviewed and illustrated with examples from the glutathione transferase gene family.

An analysis of the helix-to-strand transition between peptides with identical sequence.

An analysis of peptide segments with identical sequence but that differ significantly in structure was performed over non-redundant databases of protein structures. We focus on those peptides, which fold into an alpha-helix in one protein but a beta-strand in another. While the study shows that many such structurally ambivalent peptides contain amino acids with a strong helical preference collocated with amino acids with a strong strand preference, the results overwhelmingly indicate that the peptide's environment ultimately dictates its structure. Furthermore, the first naturally occurring structurally ambivalent nonapeptide from evolutionary unrelated proteins is described, highlighting the intrinsic plasticity of peptide sequences. We even find seven proteins that show structural ambivalence under different conditions. Finally, a computer algorithm has been implemented to identify regions in a given sequence where secondary structure prediction programs are likely to make serious mispredictions.

Proteolytic processing of peptides in the lumen of the endoplasmic reticulum for antigen presentation by major histocompatibility class I.

We have tested the hypothesis that MHC class I molecules are actively involved as protease in the production of natural MHC class I ligands. First, the structure of a class I molecule was analyzed for homology with catalytic sites of known proteases. While several clusters of amino acids in the restriction element resembled protease active sites, structural discrepancies and the influence of nearby residues suggest that these sites are unlikely to have protease activity. Second, we have tested the presentation of viral cytotoxic T cell determinants with affinity for the same restriction element (H-2K(d) or K(k)), when targeted as tandem peptides into the endoplasmic reticulum. Peptide transporter-defective cells were used to exclude cleavage of the tandem peptides by cytosolic proteases. Cleavage by signal peptidase of the tandem peptides was ascertained. The C-terminal peptides in the tandem arrays were almost exclusively presented, suggesting that an aminopeptidase in the endoplasmic reticulum degraded the N-terminally positioned peptides. This result is inconsistent with an MHC class I-catalyzed cleavage following binding of longer peptides in the cleft of the restriction elements. Finally, we conclusively show that an aminopeptidase in the endoplasmic reticulum is also involved in antigen presentation in cells with a functional peptide transporter.

Identification, characterization, and crystal structure of the Omega class glutathione transferases.

A new class of glutathione transferases has been discovered by analysis of the expressed sequence tag data base and sequence alignment. Glutathione S-transferases (GSTs) of the new class, named Omega, exist in several mammalian species and Caenorhabditis elegans. In humans, GSTO 1-1 is expressed in most tissues and exhibits glutathione-dependent thiol transferase and dehydroascorbate reductase activities characteristic of the glutaredoxins. The structure of GSTO 1-1 has been determined at 2.0-A resolution and has a characteristic GST fold (Protein Data Bank entry code ). The Omega class GSTs exhibit an unusual N-terminal extension that abuts the C terminus to form a novel structural unit. Unlike other mammalian GSTs, GSTO 1-1 appears to have an active site cysteine that can form a disulfide bond with glutathione.

Ab initio calculations on hidden modulators of theta class glutathione transferase activity.

The glutathione transferases decrease the pKa of glutathione, allowing its deprotonation and the formation of the more reactive thiolate anion. The thiolate is maintained in the active site through a weak conventional hydrogen bond first sphere interaction donated by a Tyr hydroxyl in the Alpha, Mu, Pi, and Sigma glutathione transferase classes that can be modified by other second sphere or indirect thiolate contacts. However, the Theta and Delta class isoforms use a Ser hydroxyl for stabilizing the GSH thiolate, and as such, have a different chemical system compared with that of the Tyr possessed by other classes. We have used high level ab initio methods to investigate this interaction by using a simple methanol methanethiol system as a model. The hydrogen bond strength of this initial first sphere interaction was calculated to be less than that of the Tyr interaction. A putative second sphere interaction exists in the Theta and Delta class structures between Cys or Ser-14 and Ser-11 in the mammalian Theta subclass 1 and 2, respectively. The effect of this interaction on the first sphere interaction has also been investigated and found to significantly increase the energy of the bond.

Mutagenic analysis of conserved arginine residues in and around the novel sulfate binding pocket of the human Theta class glutathione transferase T2-2.

The human Theta class glutathione transferase GSTT2-2 has a novel sulfatase activity that is not dependent on the presence of a conserved hydrogen bond donor in the active site. Initial homology modeling and the crystallographic studies have identified three conserved Arg residues that contribute to the formation of (Arg107 and Arg239), and entry to (Arg242), a sulfate binding pocket. These residues have been individually mutated to Ala to investigate their potential role in substrate binding and catalysis. The mutation of Arg107 had a significant detrimental effect on the sulfatase reaction, while the Arg242 mutation caused only a small reduction in sulfatase activity. Surprisingly, the Arg239 had an increased activity resulting from a reduction in stability. Thus, Arg239 appears to play a role in maintaining the architecture of the active site. Electrostatic calculations performed on the wild-type and mutant forms of the enzyme are in good agreement with the experimental results. These findings, along with docking studies, suggest that prior to conjugation, the location of 1-menaphthyl sulfate, a model substrate for the sulfatase reaction, is approximately midway between the position ultimately occupied by the naphthalene ring of 1-menaphthylglutathione and the free sulfate. It is further proposed that the Arg residues in and around the sulfate binding pocket have a role in electrostatic substrate recognition.

Peptide binding motifs and specificities for HLA-DQ molecules.

HLA-DQ molecules have been associated with susceptibility to a number of autoimmune and other diseases, possibly through the peptide repertoire that can be presented by different allelic products. It is thus of importance to understand which peptides can be bound by different HLA-DQ allelic products. Recently, a model for HLA-DQ has been described and used to derive peptide positional environments for HLA-DQ allelic products. By combining the peptide positional environments with known HLA-DQ peptide binding motifs, a set of predictions of likely anchor motifs for many of the products of HLA-DQ allelic variants are made and presented in a table referred to as a roadmap for HLA-DQ peptide binding specificities.

The Ubp6 family of deubiquitinating enzymes contains a ubiquitin-like domain: SUb.

A sequence motif that is Similar to Ubiquitin (SUb) has been identified in the Saccharomyces cerevisiae ubiquitin-specific protease Ubp6. SUb is conserved in all known Ubp6 homologues from a spectrum of eukaryotic species and is also present in a group of hypothetical proteins of unknown function (Unk1-3) present in sequence databases. An N-terminal deletion mutant of Ubp6 that lacks SUb is still capable of cleaving alpha-linked ubiquitin fusions, suggesting that SUb forms a separate domain to the catalytic core of Ubp6 and demonstrating that it is not required for in vitro cleavage activity. A homology model of the 78 N-terminal amino acids of human Ubp6, based on the known fold of ubiquitin, is presented. In human Ubp6, SUb shares only 20% sequence identity with ubiquitin. Even weaker similarity occurs between S. cerevisiae SUb and ubiquitin. The homology model supports a ubiquitin-like fold for SUb and suggests that two conserved Lys residues, corresponding to Lys48 and Lys63 of ubiquitin, are functionally important.

Identification and characterization of two missense mutations causing factor XIIIA deficiency.

In this study, two amino acid substitutions, Arg260His and Val414Phe, have been identified in the factor XIIIA subunits of factor XIII deficient patients of Syrian and Indian descent, respectively. To confirm the deleterious effects of these substitutions, both variant sequences have been engineered into cDNA clones and the mutant enzymes expressed in yeast. Determination of the transglutaminase activity and immuno detection of the mutant enzymes together with mRNA hybridization revealed that the mutations dramatically reduce both the catalytic activity and the level of enzyme expressed in yeast. The mutations Arg260His and Val414Phe occur within the 'core' domain of the enzyme. Computer modelling of the mutant enzymes reveals that the substitution of the Arg260 by His results in the loss of a conserved electrostatic interaction whereas the effect of the Val414Phe substitution is a consequence of the large increase in side-chain volume. Although both mutations do not effect the active site directly, they are predicted to reduce the stability of the enzyme. The effects of these two amino acid substitutions on enzyme expression and three-dimensional structure strongly confirm that residues which are located outside of the active site can have a significant effect on protein stability and function.

A homology model for the human theta-class glutathione transferase T1-1.

A manual threading approach is used to model the human glutathione transferase T1-1 based on the coordinates of the related Theta class enzyme T2-2. The low level of sequence identity (about 20%), found in the C-terminal extension in conjunction with a relative deletion of about five residues makes this a challenging modeling problem. The C-terminal extension contributes to the active site of the molecule and is thus of particular interest for understanding the molecular mechanism of the enzyme. Manual docking of known substrates and non-substrates has implicated potential candidates for the T1-1 catalytic residues involved in the dehalogenation and epoxide-ring opening activities. These include the conserved Theta class residues Arg 107, Trp 115, and the conserved GSTT1 subclass residue His 176. Also, the residue at position 234 is implicated in the modulation of T1-1 activity with different substrates between species.

A novel Asn344 deletion in the core domain of coagulation factor XIII A subunit: its effects on protein structure and function.

In this study a previously undescribed 3 bp deletion, AAT1030-1032, in the factor XIII A subunit gene, has been detected in a Thai patient. The inframe deletion results in the translation of a factor XIII A subunit that lacks Asn344. This is the first inframe deletion to be identified in the factor XIII A subunit gene because six previously reported deletions have all caused frameshifts. The deletion has been introduced into a factor XIII A subunit cDNA and the deleted polypeptide expressed in yeast. The mRNA encoding the mutant enzyme appears to have normal stability but the translated protein is subject to premature degradation. In addition, the mutated enzyme exhibited very little transglutaminase activity compared with the wild-type enzyme. Structural modeling of the deleted enzyme suggests that the absence of Asn344 would have a potent impact on the catalytic activity by reorienting the residues associated with the catalytic center. Thus, the Asn344 deletion strongly confirms the significance of the residues surrounding the catalytic center of the factor XIII A subunit.

A viral peptide with limited homology to a self peptide can induce clinical signs of experimental autoimmune encephalomyelitis.

Molecular mimicry has been suggested as a mode of autoreactive T cell stimulation in autoimmune diseases. Myelin basic protein (MBP) peptide 1-11 induces experimental autoimmune encephalomyelitis (EAE) in susceptible strains of mice. Here we show that a herpesvirus Saimiri (HVS) peptide, AAQRRPSRPFA, with a limited homology to MBP1-11 peptide, ASQKRPSQRHG (underlined letters showing homology), can stimulate a panel of MBP-11-specific T cell hybridomas and more importantly cause EAE in mice. We demonstrate that this is due to cross-recognition of these two peptides by TCRs. Results presented in this communication are the first demonstration that a viral peptide with homology at just 5 amino acids with a self peptide can induce clinical signs of EAE in mice. These findings have important implications in understanding the breakdown of T cell tolerance to self Ags in autoimmune diseases by means of cross-reactivity with unrelated peptides.

Anti-MUC1 antibodies react directly with MUC1 peptides presented by class I H2 and HLA molecules.

Peptides bound in the groove of MHC class I molecules and detected by CTLs are not normally accessible to Ab. We now report that MUC1 peptides that are bound within the groove of MHC class I molecules (H2 and HLA) and that can be detected by CTLs can also be detected by anti-MUC1 Abs. mAbs to the middle and C-terminal regions of the class I-associated peptides but not to the N terminus were able to react with MUC1 peptides bound to H2Kb and HLA-A*0201, and only to the mid-region for H2Db, by flow cytometry and also to block CTL activity. Molecular modeling showed that the N terminus is buried (and not accessible), whereas the midpeptide residues form a loop and the C terminus is free, making these two regions accessible to Ab. The findings demonstrate for the first time that peptides associated with class I molecules can be detected by anti-peptide Abs.

A combinatorial distance-constraint approach to predicting protein tertiary models from known secondary structure.

BACKGROUND: Distance geometry methods allow protein structures to be constructed using a large number of distance constraints, which can be elucidated by experimental techniques such as NMR. New methods for gleaning tertiary structural information from multiple sequence alignments make it possible for distance constraints to be predicted from sequence information alone. The basic distance geometry method can thus be applied using these empirically derived distance constraints. Such an approach, which incorporates a novel combinatoric procedure, is reported here. RESULTS: Given the correct sheet topology and disulfide formations, the fully automated procedure is generally able to construct native-like Calpha models for eight small beta-protein structures. When the sheet topology was unknown but disulfide connectivities were included, all sheet topologies were explored by the combinatorial procedure. Using a simple geometric evaluation scheme, models with the correct sheet topology were ranked first in four of the eight example cases, second in three examples and third in one example. If neither the sheet topology nor the disulfide connectivities were given a priori, all combinations of sheet topologies and disulfides were explored by the combinatorial procedure. The evaluation scheme ranked the correct topology within the top five folds for half the example cases. CONCLUSIONS: The combinatorial procedure is a useful technique for identifying a limited number of low-resolution candidate folds for small, disulfide-rich, beta-protein structures. Better results are obtained, however, if correct disulfide connectivities are known in advance. Combinatorial distance constraints can be applied whenever there are a sufficiently small number of finite connectivities.