Molecular basis of glutathione synthetase deficiency and a rare gene permutation event

Cis-eQTLs in seven duck tissues identify novel candidate genes for growth and carcass traits

BACKGROUND

Expression quantitative trait loci (eQTL) studies aim to understand the influence of genetic variants on gene expression. The colocalization of eQTL mapping and GWAS strategy could help identify essential candidate genes and causal DNA variants vital to complex traits in human and many farm animals. However, eQTL mapping has not been conducted in ducks. It is desirable to know whether eQTLs within GWAS signals contributed to duck economic traits.

RESULTS

In this study, we conducted an eQTL analysis using publicly available RNA sequencing data from 820 samples, focusing on liver, muscle, blood, adipose, ovary, spleen, and lung tissues. We identified 113,374 cis-eQTLs for 12,266 genes, a substantial fraction 39.1% of which were discovered in at least two tissues. The cis-eQTLs of blood were less conserved across tissues, while cis-eQTLs from any tissue exhibit a strong sharing pattern to liver tissue. Colocalization between cis-eQTLs and genome-wide association studies (GWAS) of 50 traits uncovered new associations between gene expression and potential loci influencing growth and carcass traits. SRSF4, GSS, and IGF2BP1 in liver, NDUFC2 in muscle, ELF3 in adipose, and RUNDC1 in blood could serve as the candidate genes for duck growth and carcass traits.

CONCLUSIONS

Our findings highlight substantial differences in genetic regulation of gene expression across duck primary tissues, shedding light on potential mechanisms through which candidate genes may impact growth and carcass traits. Furthermore, this availability of eQTL data offers a valuable resource for deciphering further genetic association signals that may arise from ongoing extensive endeavors aimed at enhancing duck production traits.

Targeting prostate tumor low-molecular weight tyrosine phosphatase for oxidation-sensitizing therapy

Protein tyrosine phosphatases (PTPs) play major roles in cancer and are emerging as therapeutic targets. Recent reports suggest low-molecular weight PTP (LMPTP)-encoded by the ACP1 gene-is overexpressed in prostate tumors. We found ACP1 up-regulated in human prostate tumors and ACP1 expression inversely correlated with overall survival. Using CRISPR-Cas9-generated LMPTP knockout C4-2B and MyC-CaP cells, we identified LMPTP as a critical promoter of prostate cancer (PCa) growth and bone metastasis. Through metabolomics, we found that LMPTP promotes PCa cell glutathione synthesis by dephosphorylating glutathione synthetase on inhibitory Tyr270. PCa cells lacking LMPTP showed reduced glutathione, enhanced activation of eukaryotic initiation factor 2-mediated stress response, and enhanced reactive oxygen species after exposure to taxane drugs. LMPTP inhibition slowed primary and bone metastatic prostate tumor growth in mice. These findings reveal a role for LMPTP as a critical promoter of PCa growth and metastasis and validate LMPTP inhibition as a therapeutic strategy for treating PCa through sensitization to oxidative stress.

Environmentally relevant concentrations of tris (2-chloroethyl) phosphate (TCEP) induce hepatotoxicity in zebrafish (Danio rerio): a whole life-cycle assessment

Tris (2-chloroethyl) phosphate (TCEP), a typical organophosphate flame retardant, is of increasingly great concern considering their ubiquitous presence in aquatic environments and potential ecotoxicity. The present work was aimed to investigate the potential growth inhibition and hepatic stress induced by whole life-cycle exposure to TCEP (0.8, 4, 20 and 100 μg/L) in zebrafish. The results revealed that the body length, body mass and hepatic-somatic index (HSI) of zebrafish were significantly declined after exposure to TCEP for 120 days. GPx activity and GSH content were increased in the liver of zebrafish treated with low concentrations (0.8 and 4 μg/L) of TCEP, while exposure to high concentrations (20 and 100 μg/L) of TCEP reduced antioxidative capacity and elevated lipid peroxidation (LPO) levels. Gene transcription analysis demonstrated that the mRNA levels of nrf2 were altered in a similar manner to the transcription of the downstream genes nqo1 and hmox1, suggesting that Nrf2-Keap1 pathway mediated TCEP-induced oxidative stress in zebrafish liver. In addition, TCEP exposure might alleviate inflammatory response through down-regulating transcription of inflammatory cytokines (il-1β, il-6 and inos), and induce apoptosis via activating the p53-Bax pathway. Moreover, whole life-cycle exposure to TCEP caused a series of histopathological anomalies in zebrafish liver. Overall, our results revealed that lifetime exposure to environmentally relevant concentrations of TCEP could result in growth retardation and induce significant hepatotoxicity in zebrafish.

Pharmacoinformatics approach based identification of potential Nsp15 endoribonuclease modulators for SARS-CoV-2 inhibition

In the current study, a structure-based virtual screening paradigm was used to screen a small molecular database against the Non-structural protein 15 (Nsp15) endoribonuclease of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). The SARS-CoV-2 is the causative agent of the recent outbreak of coronavirus disease 2019 (COVID-19) which left the entire world locked down inside the home. A multi-step molecular docking study was performed against antiviral specific compounds (~8722) collected from the Asinex antiviral database. The less or non-interacting molecules were wiped out sequentially in the molecular docking. Further, MM-GBSA based binding free energy was estimated for 26 compounds which shows a high affinity towards the Nsp15. The drug-likeness and pharmacokinetic parameters of all 26 compounds were explored, and five molecules were found to have an acceptable pharmacokinetic profile. Overall, the Glide-XP docking score and Prime-MM-GBSA binding free energy of the selected molecules were explained strong interaction potentiality towards the Nsp15 endoribonuclease. The dynamic behavior of each molecule with Nsp15 was assessed using conventional molecular dynamics (MD) simulation. The MD simulation information was strongly favors the Nsp15 and each identified ligand stability in dynamic condition. Finally, from the MD simulation trajectories, the binding free energy was estimated using the MM-PBSA method. Hence, the proposed final five molecules might be considered as potential Nsp15 modulators for SARS-CoV-2 inhibition.

Structural biology of plant sulfur metabolism: from sulfate to glutathione

Sulfur is an essential element for all organisms. Plants must assimilate this nutrient from the environment and convert it into metabolically useful forms for the biosynthesis of a wide range of compounds, including cysteine and glutathione. This review summarizes structural biology studies on the enzymes involved in plant sulfur assimilation [ATP sulfurylase, adenosine-5'-phosphate (APS) reductase, and sulfite reductase], cysteine biosynthesis (serine acetyltransferase and O-acetylserine sulfhydrylase), and glutathione biosynthesis (glutamate-cysteine ligase and glutathione synthetase) pathways. Overall, X-ray crystal structures of enzymes in these core pathways provide molecular-level information on the chemical events that allow plants to incorporate sulfur into essential metabolites and revealed new biochemical regulatory mechanisms, such as structural rearrangements, protein-protein interactions, and thiol-based redox switches, for controlling different steps in these pathways.

Genetic Mutations in the S-loop of Human Glutathione Synthetase: Links Between Substrate Binding, Active Site Structure and Allostery

The second step in the biosynthesis of the cellular antioxidant glutathione (GSH) is catalyzed by human glutathione synthetase (hGS), a negatively cooperative homodimer. Patients with mutations in hGS have been reported to exhibit a range of symptoms from hemolytic anemia and metabolic acidosis to neurological disorders and premature death. Several patient mutations occur in the S-loop of hGS, a series of residues near the negatively cooperative γ-GC substrate binding site. Experimental point mutations and molecular dynamic simulations show the S-loop not only binds γ-GC through a salt bridge and multiple hydrogen bonds, but the residues also modulate allosteric communication in hGS. By elucidating the role of S-loop residues in active site structure, substrate binding, and allostery, the atomic level sequence of events that leads to the detrimental effects of hGS mutations in patients are more fully understood.

Oxidative stress, cell cycle arrest, DNA damage and apoptosis in adult zebrafish (Danio rerio) induced by tris(1,3-dichloro-2-propyl) phosphate

Tris(1,3-dichloro-2-propyl)phosphate (TDCPP) is an additive flame retardant of high production volume, and frequently detected in biota and environment. However, knowledge on its potential risk and toxicological mechanism still remains limited. In this study, DNA damage, transcriptomic responses and biochemical changes in the liver of zebrafish (Danio rerio) induced by TDCPP were investigated. Zebrafish was exposed to 45.81μg/L (1/100 (96h-LC₅₀)) and 229.05μg/L (1/20 (96h-LC₅₀)) TDCPP for 7 d. The reactive oxygen species (ROS) and GSH contents, in addition to antioxidant enzyme activities in the liver changed significantly, and the mRNA levels of genes related to oxidative stress were alerted in a dose-dependent and/or sex-dependent manner after exposure to TDCPP. Significant DNA damage in zebrafish liver was found, and olive tail moment increased in a concentration-dependent manner. Moreover, exposure of TDCPP at 45.81μg/L level activated the cell cycle arrest, DNA repair system and apoptosis pathway in male zebrafish, and 229.05μg/L TDCPP exposure inhibited those pathways in both male and female zebrafish. The cell apoptosis was confirmed in TUNEL assay as higher incidence of TUNEL-positive cells were observed in zebrafish exposed to 229.05μg/L TDCPP. Our results also indicated that males were more sensitive to TDCPP exposure compared with females. Taken together, our results showed that TDCPP could induce oxidative stress, cell cycle arrest, DNA damage and apoptosis in adult zebrafish liver in sex- and concentration-dependent manners.

Drug Repurposing of the Anthelmintic Niclosamide to Treat Multidrug-Resistant Leukemia

Multidrug resistance, a major problem that leads to failure of anticancer chemotherapy, requires the development of new drugs. Repurposing of established drugs is a promising approach for overcoming this problem. An example of such drugs is niclosamide, a known anthelmintic that is now known to be cytotoxic and cytostatic against cancer cells. In this study, niclosamide showed varying activity against different cancer cell lines. It revealed better activity against hematological cancer cell lines CCRF-CEM, CEM/ADR5000, and RPMI-8226 compared to the solid tumor cell lines MDA-MB-231, A549, and HT-29. The multidrug resistant CEM/ADR5000 cells were similar sensitive as their sensitive counterpart CCRF-CEM (resistance ration: 1.24). Furthermore, niclosamide caused elevations in reactive oxygen species and glutathione (GSH) levels in leukemia cells. GSH synthetase (GS) was predicted as a target of niclosamide. Molecular docking showed that niclosamide probably binds to the ATP-binding site of GS with a binding energy of -9.40 kcal/mol. Using microscale thermophoresis, the binding affinity between niclosamide and recombinant human GS was measured (binding constant: 5.64 μM). COMPARE analyses of the NCI microarray database for 60 cell lines showed that several genes, including those involved in lipid metabolism, correlated with cellular responsiveness to niclosamide. Hierarchical cluster analysis showed five major branches with significant differences between sensitive and resistant cell lines (p = 8.66 × 10⁵). Niclosamide significantly decreased nuclear factor of activated T-cells (NFAT) activity as predicted by promoter binding motif analysis. In conclusion, niclosamide was more active against hematological malignancies compared to solid tumors. The drug was particularly active against the multidrug-resistant CEM/ADR5000 leukemia cells. Inhibition of GSH synthesis and NFAT signaling were identified as relevant mechanisms for the anticancer activity of niclosamide. Gene expression profiling predicted the sensitivity or resistance of cancer cells to niclosamide.

Effects of TBEP on the induction of oxidative stress and endocrine disruption in Tm3 Leydig cells

The flame retardant tris (2-butoxyethyl) phosphate (TBEP) is a frequently detected contaminant in the environment. In the cultured TM3 cells (originated from ATCC), effects of TBEP on the induction of oxidative stress and endocrine disruption were evaluated. It was observed that exposure to 100 μg/mL TBEP for 24 h significantly reduced the viability of TM3 cells. The mRNA levels of genes related to oxidative stress including Sod, Gpx1, Cat, and Gsta1 were changed in a dose-dependent and/or time-dependent manner after exposure to 30 and 100 μg/mL TBEP for 6, 12, or 24 h. Moreover, notable decrease in glutathione (GSH) contents and increases in oxidized glutathione (GSSG) contents as well as the antioxidant enzyme activities like superoxide dismutase, catalase, glutathione peroxidase, and glutathione S-transferase were found in the group treated with 100 μg/mL TBEP for 24 h, indicating that TBEP induced oxidative stress in TM3 Leydig cells. In addition, the expression of genes related to testosterone (T) synthesis including cytochrome P450 cholesterol side-chain cleavage enzyme (P450scc), cytochrome P450 17α-hydroxysteroid dehydrogenase (P450-17α), and 17β-hydroxysteroid dehydrogenase (17β-HSD) and T levels in medium were remarkably declined by the treatment of 100 μg/mL TBEP for 24 h. And TBEP could inhibit the expression of P450-17α and 17β-HSD and T levels up-regulated by hCG in TM3 cells. Taken together, these findings indicated that TBEP can induce oxidative stress and alter steroidogenesis in TM3 cells. © 2015 Wiley Periodicals, Inc. Environ Toxicol 31: 1276-1286, 2016.

The acetaminophen metabolite N-acetyl-p-benzoquinone imine (NAPQI) inhibits glutathione synthetase in vitro; a clue to the mechanism of 5-oxoprolinuric acidosis?

1. Metabolic acidosis due to accumulation of l-5-oxoproline is a rare, poorly understood, disorder associated with acetaminophen treatment in malnourished patients with chronic morbidity. l-5-Oxoprolinuria signals abnormal functioning of the γ-glutamyl cycle, which recycles and synthesises glutathione. Inhibition of glutathione synthetase (GS) by N-acetyl-p-benzoquinone imine (NAPQI) could contribute to 5-oxoprolinuric acidosis in such patients. We investigated the interaction of NAPQI with GS in vitro. 2. Peptide mapping of co-incubated NAPQI and GS using mass spectrometry demonstrated binding of NAPQI with cysteine-422 of GS, which is known to be essential for GS activity. Computational docking shows that NAPQI is properly positioned for covalent bonding with cysteine-422 via Michael addition and hence supports adduct formation. 3. Co-incubation of 0.77 μM of GS with NAPQI (25-400 μM) decreased enzyme activity by 16-89%. Inhibition correlated strongly with the concentration of NAPQI and was irreversible. 4. NAPQI binds covalently to GS causing irreversible enzyme inhibition in vitro. This is an important novel biochemical observation. It is the first indication that NAPQI may inhibit glutathione synthesis, which is pivotal in NAPQI detoxification. Further studies are required to investigate its biological significance and its role in 5-oxoprolinuric acidosis.

Structure of NDP-forming Acetyl-CoA synthetase ACD1 reveals a large rearrangement for phosphoryl transfer

The NDP-forming acyl-CoA synthetases (ACDs) catalyze the conversion of various CoA thioesters to the corresponding acids, conserving their chemical energy in form of ATP. The ACDs are the major energy-conserving enzymes in sugar and peptide fermentation of hyperthermophilic archaea. They are considered to be primordial enzymes of ATP synthesis in the early evolution of life. We present the first crystal structures, to our knowledge, of an ACD from the hyperthermophilic archaeon Candidatus Korachaeum cryptofilum. These structures reveal a unique arrangement of the ACD subunits alpha and beta within an α2β2-heterotetrameric complex. This arrangement significantly differs from other members of the superfamily. To transmit an activated phosphoryl moiety from the Ac-CoA binding site (within the alpha subunit) to the NDP-binding site (within the beta subunit), a distance of 51 Å has to be bridged. This transmission requires a larger rearrangement within the protein complex involving a 21-aa-long phosphohistidine-containing segment of the alpha subunit. Spatial restraints of the interaction of this segment with the beta subunit explain the necessity for a second highly conserved His residue within the beta subunit. The data support the proposed four-step reaction mechanism of ACDs, coupling acyl-CoA thioesters with ATP synthesis. Furthermore, the determined crystal structure of the complex with bound Ac-CoA allows first insight, to our knowledge, into the determinants for acyl-CoA substrate specificity. The composition and size of loops protruding into the binding pocket of acyl-CoA are determined by the individual arrangement of the characteristic subdomains.

The role of strong electrostatic interactions at the dimer interface of human glutathione synthetase

The obligate homodimer human glutathione synthetase (hGS) provides an ideal system for exploring the role of protein-protein interactions in the structural stability, activity and allostery of enzymes. The two active sites of hGS, which are 40 Å apart, display allosteric modulation by the substrate γ-glutamylcysteine (γ-GC) during the synthesis of glutathione, a key cellular antioxidant. The two subunits interact at a relatively small dimer interface dominated by electrostatic interactions between S42, R221, and D24. Alanine scans of these sites result in enzymes with decreased activity, altered γ-GC affinity, and decreased thermal stability. Molecular dynamics simulations indicate these mutations disrupt interchain bonding and impact the tertiary structure of hGS. While the ionic hydrogen bonds and salt bridges between S42, R221, and D24 do not mediate allosteric communication in hGS, these interactions have a dramatic impact on the activity and structural stability of the enzyme.

The biosynthesis of UDP-d-FucNAc-4N-(2)-oxoglutarate (UDP-Yelosamine) in Bacillus cereus ATCC 14579: Pat and Pyl, an aminotransferase and an ATP-dependent Grasp protein that ligates 2-oxoglutarate to UDP-4-amino-sugars

Surface glycan switching is often observed when micro-organisms transition between different biotic and abiotic niches, including biofilms, although the advantages of this switching to the organism are not well understood. Bacillus cereus grown in a biofilm-inducing medium has been shown to synthesize an unusual cell wall polysaccharide composed of the repeating subunit →6)Gal(α1-2)(2-R-hydroxyglutar-5-ylamido)Fuc2NAc4N(α1-6)GlcNAc(β1→, where galactose is linked to the hydroxyglutarate moiety of FucNAc-4-amido-(2)-hydroxyglutarate. The molecular mechanism involved in attaching 2-hydroxyglutarate to 4-amino-FucNAc has not been determined. Here, we show two genes in B. cereus ATCC 14579 encoding enzymes involved in the synthesis of UDP-FucNAc-4-amido-(2)-oxoglutarate (UDP-Yelosamine), a modified UDP-sugar not previously reported to exist. Using mass spectrometry and real time NMR spectroscopy, we show that Bc5273 encodes a C4″-aminotransferase (herein referred to as Pat) that, in the presence of pyridoxal phosphate, transfers the primary amino group of l-Glu to C-4″ of UDP-4-keto-6-deoxy-d-GlcNAc to form UDP-4-amino-FucNAc and 2-oxoglutarate. Pat also converts 4-keto-xylose, 4-keto-glucose, and 4-keto-2-acetamido-altrose to their corresponding UDP-4-amino-sugars. Bc5272 encodes a carboxylate-amine ligase (herein referred as Pyl) that, in the presence of ATP and Mg(II), adds 2-oxoglutarate to the 4-amino moiety of UDP-4-amino-FucNAc to form UDP-Yelosamine and ADP. Pyl is also able to ligate 2-oxoglutarate to other 4-amino-sugar derivatives to form UDP-Yelose, UDP-Solosamine, and UDP-Aravonose. Characterizing the metabolic pathways involved in the formation of modified nucleotide sugars provides a basis for understanding some of the mechanisms used by bacteria to modify or alter their cell surface polysaccharides in response to changing growth and environmental challenges.

Resilience of biochemical activity in protein domains in the face of structural divergence

Recent studies point to the prevalence of the evolutionary phenomenon of drastic structural transformation of protein domains while continuing to preserve their basic biochemical function. These transformations span a wide spectrum, including simple domains incorporated into larger structural scaffolds, changes in the structural core, major active site shifts, topological rewiring and extensive structural transmogrifications. Proteins from biological conflict systems, such as toxin-antitoxin, restriction-modification, CRISPR/Cas, polymorphic toxin and secondary metabolism systems commonly display such transformations. These include endoDNases, metal-independent RNases, deaminases, ADP ribosyltransferases, immunity proteins, kinases and E1-like enzymes. In eukaryotes such transformations are seen in domains involved in chromatin-related peptide recognition and protein/DNA-modification. Intense selective pressures from 'arms-race'-like situations in conflict and macromolecular modification systems could favor drastic structural divergence while preserving function.

Structural insights into the enzymes of the trypanothione pathway: targets for antileishmaniasis drugs

Leishmaniasis is a neglected disease that kills 60,000 people worldwide, and which is caused by the protozoa Leishmania. The enzymes of the trypanothione pathway: trypanothione synthetase-amidase, trypanothione reductase (TR) and tryparedoxin-dependent peroxidase are absent in human hosts, and are essential for parasite survival and druggable. The most promising target is trypanothione synthetase-amidase, which has been also chemically validated. However, the structural data presented in this review show that TR also should be considered as a good target. Indeed, it is strongly inhibited by silver- and gold-containing compounds, which are active against Leishmania parasites and can be used for the development of novel antileishmanial agents. Moreover, TR trypanothione-binding site is not featureless but contains a sub-pocket where inhibitors bind, potentially useful for the design of new lead compounds.

Single mutation alters the substrate specificity of L-amino acid ligase

L-Amino acid ligase (Lal) catalyzes the formation of a dipeptide from two L-amino acids in an ATP-dependent manner and belongs to the ATP-grasp superfamily. Bacillus subtilis YwfE, the first identified Lal, produces the dipeptide antibiotic bacilysin, which consists of L-Ala and L-anticapsin. Its substrate specificity is restricted to smaller amino acids such as L-Ala for the N-terminal end of the dipeptide, whereas a wide range of hydrophobic amino acids including L-Phe and L-Met are recognized for the C-terminal end in vitro. We determined the crystal structures of YwfE with bound ADP-Mg(2+)-Pi and ADP-Mg(2+)-L-Ala at 1.9 and 2.0 Å resolutions, respectively. On the basis of these structures, we generated point mutants of residues that are considered to participate in the recognition of L-Ala and measured their ATPase activity. The conserved Arg328 is suggested to be a crucial residue for L-Ala recognition and catalysis. The mutation of Trp332 to Ala caused the enzyme to hydrolyze ATP, even in the absence of l-Ala, and the structure of this mutant protein appeared to show a cavity in the N-terminal substrate-binding pocket. These results suggest that Trp332 plays a key role in restricting the substrate specificity to smaller amino acids such as L-Ala. Moreover, Trp332 mutants can alter the substrate specificity and activity depending on the size and shape of substituted amino acids. These observations provide sufficient scope for the rational design of Lal to produce desirable dipeptides. We propose that the positioning of the conserved Arg residue in Lal is important for enantioselective recognition of L-amino acids.

S-glutathionylation: relevance in diabetes and potential role as a biomarker

Glutathione is considered the main regulator of redox balance in the cellular milieu due to its capacity for detoxifying deleterious molecules. The oxidative stress induced as a result of a variety of stimuli promotes protein oxidation, usually at cysteine residues, leading to changes in their activity. Mild oxidative stress, which may take place in physiological conditions, induces the reversible oxidation of cysteines to sulfenic acid form, while pathological conditions are associated with higher rates of reactive oxygen species production, inducing the irreversible oxidation of cysteines. Among these, neurodegenerative disorders, cardiovascular diseases and diabetes have been proposed to be pathogenetically linked to this state. In diabetes-associated vascular complications, lower levels of glutathione and increased oxidative stress have been reported. S-glutathionylation has been proposed as a posttranslational modification able to protect proteins from over-oxidizing environments. S-glutathionylation has been identified in proteins involved in diabetic models both in vitro and in vivo. In all of them, S-glutathionylation represents a mechanism that regulates the response to diabetic conditions, and has been described to occur in erythrocytes and neutrophils from diabetic patients. However, additional studies are necessary to discern whether this modification represents a biomarker for the early onset of diabetic vascular complications.

Inhibition of Raf-MEK-ERK and hypoxia pathways by Phyllanthus prevents metastasis in human lung (A549) cancer cell line

BACKGROUND

Lung cancer constitutes one of the malignancies with the greatest incidence and mortality rates with 1.6 million new cases and 1.4 million deaths each year. Prognosis remains poor due to deleterious development of multidrug resistance resulting in less than 15% lung cancer patients reaching five years survival. We have previously shown that Phyllanthus induced apoptosis in conjunction with its antimetastastic action. In the current study, we aimed to determine the signaling pathways utilized by Phyllanthus to exert its antimetastatic activities.

METHODS

Cancer 10-pathway reporter array was performed to screen the pathways affected by Phyllanthus in lung carcinoma cell line (A549) to exert its antimetastatic effects. Results from this array were then confirmed with western blotting, cell cycle analysis, zymography technique, and cell based ELISA assay for human total iNOS. Two-dimensional gel electrophoresis was subsequently carried out to study the differential protein expressions in A549 after treatment with Phyllanthus.

RESULTS

Phyllanthus was observed to cause antimetastatic activities by inhibiting ERK1/2 pathway via suppression of Raf protein. Inhibition of this pathway resulted in the suppression of MMP2, MMP7, and MMP9 expression to stop A549 metastasis. Phyllanthus also inhibits hypoxia pathway via inhibition of HIF-1α that led to reduced VEGF and iNOS expressions. Proteomic analysis revealed a number of proteins downregulated by Phyllanthus that were involved in metastatic processes, including invasion and mobility proteins (cytoskeletal proteins), transcriptional proteins (proliferating cell nuclear antigen; zinc finger protein), antiapoptotic protein (Bcl2) and various glycolytic enzymes. Among the four Phyllanthus species tested, P. urinaria showed the greatest antimetastatic activity.

CONCLUSIONS

Phyllanthus inhibits A549 metastasis by suppressing ERK1/2 and hypoxia pathways that led to suppression of various critical proteins for A549 invasion and migration.

Transcriptional regulation of glutathione biosynthesis genes, γ-glutamyl-cysteine ligase and glutathione synthetase in response to cadmium and nonylphenol in Chironomus riparius

We characterized Chironomus riparius glutathione (GSH) biosynthesis genes, γ-glutamyl-cysteine ligase catalytic subunit (cr-gcl) and glutathione synthetase (cr-gs) and studied their expression after cadmium (Cd) and nonylphenol (NP) exposure. The full length cDNA of the Cr-GCL catalytic subunit was 2185 base pair (bp) in length containing an open reading frame of 1905bp, a 13bp 5' and 267bp 3' untranslated regions. The theoretical molecular mass of the deduced amino acid sequence (633) was 72.65kDa with an estimated pI of 5.42. The partial cDNA of Cr-GS was 739bp in length consisting 221 amino acids. The deduced amino acid sequence of Cr-GCL and Cr-GS cDNAs showed high conservation with homologs from other species. In phylogenetic analysis Cr-GCL and Cr-GS were grouped with equivalent genes from insects belonging to the dipteran order. The expression of cr-gcl and cr-gs was measured using quantitative real-time PCR after exposure to sub lethal concentrations of Cd (2, 10 and 20mg/L) and NP (10, 50 and 100μg/L) for 12, 24, 48 and 72h using real-time PCR methods. The mRNA expression of Cr-GCL and Cr-GS was significantly modulated after exposure to different concentrations of Cd and NP for different time periods. Total GSH levels showed a non-significant decrease after exposure to Cd for 24h. However, no change in GSH levels was observed after exposure to NP for 24h. These results suggest that Cr-GS and Cr-GCL expression is modulated by Cd and NP stress and may play an important role in detoxification of xenobiotics and antioxidant defense. We conclude that Cr-GS and Cr-GCL could be used as biomarkers of Cd and NP stress in aquatic environment for the studied species.

Probing the origins of glutathione biosynthesis through biochemical analysis of glutamate-cysteine ligase and glutathione synthetase from a model photosynthetic prokaryote

Glutathione biosynthesis catalysed by GCL (glutamate-cysteine ligase) and GS (glutathione synthetase) is essential for maintaining redox homoeostasis and protection against oxidative damage in diverse eukaroytes and bacteria. This biosynthetic pathway probably evolved in cyanobacteria with the advent of oxygenic photosynthesis, but the biochemical characteristics of progenitor GCLs and GSs in these organisms are largely unexplored. In the present study we examined SynGCL and SynGS from Synechocystis sp. PCC 6803 using steady-state kinetics. Although SynGCL shares ~15% sequence identity with the enzyme from plants and α-proteobacteria, sequence comparison suggests that these enzymes share similar active site residues. Biochemically, SynGCL lacks the redox regulation associated with the plant enzymes and functions as a monomeric protein, indicating that evolution of redox regulation occurred later in the green lineage. Site-directed mutagenesis of SynGCL establishes this enzyme as part of the plant-like GCL family and identifies a catalytically essential arginine residue, which is structurally conserved across all forms of GCLs, including those from non-plant eukaryotes and γ-proteobacteria. A reaction mechanism for the synthesis of γ-glutamylcysteine by GCLs is proposed. Biochemical and kinetic analysis of SynGS reveals that this enzyme shares properties with other prokaryotic GSs. Initial velocity and product inhibition studies used to examine the kinetic mechanism of SynGS suggest that it and other prokaryotic GSs uses a random ter-reactant mechanism for the synthesis of glutathione. The present study provides new insight on the molecular mechanisms and evolution of glutathione biosynthesis; a key process required for enhancing bioenergy production in photosynthetic organisms.

Site directed mutagenesis of Schizosaccharomyces pombe glutathione synthetase produces an enzyme with homoglutathione synthetase activity

Three different His-tagged, mutant forms of the fission yeast glutathione synthetase (GSH2) were derived by site-directed mutagenesis. The mutant and wild-type enzymes were expressed in E. coli DH5α and affinity purified in a two-step procedure. Analysis of enzyme activity showed that it was possible to shift the substrate specificity of GSH2 from Gly (k_m 0,19; wild-type) to β-Ala or Ser. One mutation (substitution of Ile471, Cy472 to Met and Val and Ala 485 and Thr486 to Leu and Pro) increased the affinity of GSH2 for β-Ala (k_m 0,07) and lowered the affinity for Gly (k_m 0,83), which is a characteristic of the enzyme homoglutathione synthetase found in plants. Substitution of Ala485 and Thr486 to Leu and Pro only, increased instead the affinity of GSH2 for Ser (k_m 0,23) as a substrate, while affinity to Gly was preserved (k_m 0,12). This provides a new biosynthetic pathway for hydroxymethyl glutathione, which is known to be synthesized from glutathione and Ser in a reaction catalysed by carboxypeptidase Y. The reported findings provide further insight into how specific amino acids positioned in the GSH2 active site facilitate the recognition of different amino acid substrates, furthermore they support the evolutionary theory that homoglutathione synthetase evolved from glutathione synthetase by a single gene duplication event.

Comparison of the functions of glutathionylspermidine synthetase/amidase from E. coli and its predicted homologues YgiC and YjfC

Protein function prediction is very important in establishing the roles of various proteins in bacteria; however, some proteins in the E. coli genome have their function assigned based on low percent sequence homology that does not provide reliable assignments. We have made an attempt to verify the prediction that E. coli genes ygiC and yjfC encode proteins with the same function as glutathionylspermidine synthetase/amidase (GspSA). GspSA is a bifunctional enzyme that catalyzes the ATP-dependent formation and hydrolysis of glutathionylspermidine (G-Sp), a conjugate of glutathione (GSH) and spermidine. YgiC and YjfC proteins show 51% identity between themselves and 28% identity to the synthetase domain of the GspSA enzyme. YgiC and YjfC proteins were expressed and purified, and the properties of GspSA, YgiC, and YjfC were compared. In contrast to GspSA, proteins YgiC and YjfC did not bind to G-Sp immobilized on the affinity matrix. We demonstrated that all three proteins (GspSA, YgiC and YjfC) catalyze the hydrolysis of ATP; however, YgiC and YjfC cannot synthesize G-Sp, GSH, or GSH intermediates. gsp, ygiC, and yjfC genes were eliminated from the E. coli genome to test the ability of mutant strains to synthesize G-Sp conjugate. E. coli cells deficient in GspSA do not produce G-Sp while synthesis of the conjugate is not affected in ΔygiC and ΔyjfC mutants. All together our results indicate that YgiC and YjfC are not glutathionylspermidine synthetases as predicted from the amino acid sequence analysis.

Transcription of genes involved in glutathione biosynthesis in the solitary tunicate Ciona intestinalis exposed to metals

Exposure to metals is known to generate oxidative stress risk in living organisms, which are able to respond with the induction of antioxidant defenses, both enzymatic and non-enzymatic. Glutathione (GSH) is considered to be an important cellular component involved in protecting cells, both as metal chelating agent and oxygen radical scavenger. In this work we used molecular techniques to analyze the nucleotide and predicted amino acid sequences of genes involved in GSH biosynthesis, γ-glutamyl-cysteine ligase catalytic subunit (ci-gclc), γ-glutamyl-cysteine ligase modifier subunit (ci-gclm) and GSH synthase (ci-gs) in the solitary tunicate Ciona intestinalis. We also studied the transcription of the above genes after in vivo exposure to Cd, Cu and Zn by semiquantitativ RT-PCR to improve our knowledge about the relationship between metal-induced oxidative stress and GSH production and locate mRNA expression by in situ hybridization (ISH). These genes exhibit a good level of sequence conservation with metazoan homologs generally, especially for residues important for the activity of the enzymes. Phylogenetic analyses indicate that the three enzymes evolved in different ways, Ci-GCLC and Ci-GS being mostly correlated with invertebrate proteins, Ci-GCLM being as sister group of vertebrate GCLMs. Our in silico analyses of the ci-gs and ci-gclc promoter regions revealed putative consensus sequences similar to mammalian metal-responsive elements (MRE) and antioxidant response elements (ARE), indicating that the transcription of these genes may directly depend on metals and/or reactive oxygen species. Results highlight a statistically significant increase in gene transcription, demonstrating that metal treatments have inducible effects on these genes. They can modulate gene transcription not only through MREs but also through AREs, as a consequence of metal-dependent ROS formation. The ISH location of Ci-GS and Ci-GCLC mRNAs shows that the cells most involved in glutathione biosynthesis are circulating hemocytes. The data presented here emphasize the importance of complex metal regulation of ci-gclc, ci-gclm and ci-gs transcription, which can create an efficient detoxification pathway allowing C. intestinalis to survive in continued elevated presence of metals in the environment.

Crystallization and preliminary X-ray diffraction analysis of Bacillus subtilis YwfE, an L-amino-acid ligase

Bacillus subtilis YwfE, an L-amino-acid ligase, catalyzes the formation of an α-dipeptide from L-amino acids in an ATP-dependent manner. In order to elucidate the substrate-recognition mode and the reaction mechanism of this ligase, native and selenomethionine-derivatized (SeMet) crystals of YwfE in the presence of ADP, MgCl(2) and the dipeptide L-Ala-L-Gln were obtained using the hanging-drop vapour-diffusion method. These crystals diffracted to 1.9 and 2.8 Å resolution, respectively. Preliminary SAD phase calculations using the data set from the SeMet crystal suggested that the crystal belonged to the hexagonal space group P6(5)22, with unit-cell parameters a = b = 90.85, c = 250.31 Å, and contained one molecule in the asymmetric unit with a solvent content of 57.3%.

The ATP-grasp enzymes

The ATP-grasp enzymes consist of a superfamily of 21 proteins that contain an atypical ATP-binding site, called the ATP-grasp fold. The ATP-grasp fold is comprised of two α+β domains that "grasp" a molecule of ATP between them and members of the family typically have an overall structural design containing three common conserved focal domains. The founding members of the family consist of biotin carboxylase, d-ala-d-ala ligase and glutathione synthetase, all of which catalyze the ATP-assisted reaction of a carboxylic acid with a nucleophile via the formation of an acylphosphate intermediate. While most members of the superfamily follow this mechanistic pathway, studies have demonstrated that two enzymes catalyze only the phosphoryl transfer step and thus are kinases instead of ligases. Members of the ATP-grasp superfamily are found in several metabolic pathways including de novo purine biosynthesis, gluconeogenesis, and fatty acid synthesis. Given the critical nature of these enzymes, researchers have actively sought the development of potent inhibitors of several members of the superfamily as antibacterial and anti-obseity agents. In this review, we will discuss the structure, function, mechanism, and inhibition of the ATP-grasp enzymes.

Improving protein crystal quality by the without-oil microbatch method: crystallization and preliminary X-ray diffraction analysis of glutathione synthetase from Pseudoalteromonas haloplanktis

Glutathione synthetases catalyze the ATP-dependent synthesis of glutathione from l-γ-glutamyl- l-cysteine and glycine. Although these enzymes have been sequenced and characterized from a variety of biological sources, their exact catalytic mechanism is not fully understood and nothing is known about their adaptation at extremophilic environments. Glutathione synthetase from the Antarctic eubacterium Pseudoalteromonas haloplanktis (PhGshB) has been expressed, purified and successfully crystallized. An overall improvement of the crystal quality has been obtained by adapting the crystal growth conditions found with vapor diffusion experiments to the without-oil microbatch method. The best crystals of PhGshB diffract to 2.34 Å resolution and belong to space group P2₁2₁2₁, with unit-cell parameters a = 83.28 Å, b = 119.88 Å, c = 159.82 Å. Refinement of the model, obtained using phases derived from the structure of the same enzyme from Escherichia coli by molecular replacement, is in progress. The structural determination will provide the first structural characterization of a psychrophilic glutathione synthetase reported to date.

Aspartate 458 of human glutathione synthetase is important for cooperativity and active site structure

Human glutathione synthetase (hGS) catalyzes the second ATP-dependent step in the biosynthesis of glutathione (GSH) and is negatively cooperative to the γ-glutamyl substrate. The hGS active site is composed of three highly conserved catalytic loops, notably the alanine rich A-loop. Experimental and computational investigations of the impact of mutation of Asp458 are reported, and thus the role of this A-loop residue on hGS structure, activity, negativity cooperativity and stability is defined. Several Asp458 hGS mutants (D458A, D458N and D458R) were constructed using site-directed mutagenesis and their activities determined (10%, 15% and 7% of wild-type hGS, respectively). The Michaelis-Menten constant (K(m)) was determined for all three substrates (glycine, GAB and ATP): glycine K(m) increased by 30-115-fold, GAB K(m) decreased by 8-17-fold, and the ATP K(m) was unchanged. All Asp458 mutants display a change in cooperativity from negative cooperativity to non-cooperative. All mutants show similar stability as compared to wild-type hGS, as determined by differential scanning calorimetry. The findings indicate that Asp458 is essential for hGS catalysis and that it impacts the allostery of hGS.

Valine 44 and valine 45 of human glutathione synthetase are key for subunit stability and negative cooperativity

It was hypothesized that residues Val44 and Val45 serve as important residues for human glutathione synthetase (hGS) function and stability given their location at the dimer interface of this enzyme. Computational studies suggest that mutation at Val45 has more impact on the structure and stability of hGS than does mutation at Val44. Experimentally, enzymes with mutations at the 44 and or 45 positions of hGS were prepared, purified and assayed for initial activity. Val45 position mutations (either to alanine or tryptophan) have a greater impact on enzyme activity than do mutations at Val44. Differential scanning calorimetry experiments reveal a loss of stability in all mutant enzymes, with V45 mutations being less stable than the corresponding Val44 mutations. The γ-GluABA substrate affinity remains unaltered in V44A and V45A mutant enzymes, but increases when tryptophan is introduced at either of these positions. Hill coefficients trend towards less negative cooperativity with the exception of V45W mutant hGS. These results imply that residues V44 and V45 are located along the allosteric pathway of this negatively cooperative dimeric enzyme, that their mutation impacts the allosteric pathway more than it does the active site of hGS, and that these residues (and by extension the dimer interface in which they are located) are integral to the stability of human glutathione synthetase.

Plant glutathione biosynthesis: diversity in biochemical regulation and reaction products

In plants, exposure to temperature extremes, heavy metal-contaminated soils, drought, air pollutants, and pathogens results in the generation of reactive oxygen species that alter the intracellular redox environment, which in turn influences signaling pathways and cell fate. As part of their response to these stresses, plants produce glutathione. Glutathione acts as an anti-oxidant by quenching reactive oxygen species, and is involved in the ascorbate-glutathione cycle that eliminates damaging peroxides. Plants also use glutathione for the detoxification of xenobiotics, herbicides, air pollutants (sulfur dioxide and ozone), and toxic heavy metals. Two enzymes catalyze glutathione synthesis: glutamate-cysteine ligase, and glutathione synthetase. Glutathione is a ubiquitous protective compound in plants, but the structural and functional details of the proteins that synthesize it, as well as the potential biochemical mechanisms of their regulation, have only begun to be explored. As discussed here, the core reactions of glutathione synthesis are conserved across various organisms, but plants have diversified both the regulatory mechanisms that control its synthesis and the range of products derived from this pathway. Understanding the molecular basis of glutathione biosynthesis and its regulation will expand our knowledge of this component in the plant stress response network.

The role of the glycine triad in human glutathione synthetase

Experimental kinetics and computational modeling of human glutathione synthetase (hGS) support the significant role of the G-loop glycine triad (G369, G370, G371) for activity of this ATP-grasp enzyme. Enzyme kinetic experiments indicate that G369V and G370V mutant hGS have little activity (<0.7 and 0.3%, respectively, versus wild-type hGS). However, G371V retains ∼13% of the activity of wild-type hGS. With respect to G-loop:A-loop interaction in hGS, mutations at Gly369 and Gly370 decrease ligand binding and prevent active site closure and protection. This research indicates that Gly369 and Gly370 have essential roles in hGS, while Gly371 has a lesser involvement. Implications for glycine-rich ensembles in other phosphate-binding enzymes are discussed.

Non-homologous isofunctional enzymes: a systematic analysis of alternative solutions in enzyme evolution

BACKGROUND

Evolutionarily unrelated proteins that catalyze the same biochemical reactions are often referred to as analogous - as opposed to homologous - enzymes. The existence of numerous alternative, non-homologous enzyme isoforms presents an interesting evolutionary problem; it also complicates genome-based reconstruction of the metabolic pathways in a variety of organisms. In 1998, a systematic search for analogous enzymes resulted in the identification of 105 Enzyme Commission (EC) numbers that included two or more proteins without detectable sequence similarity to each other, including 34 EC nodes where proteins were known (or predicted) to have distinct structural folds, indicating independent evolutionary origins. In the past 12 years, many putative non-homologous isofunctional enzymes were identified in newly sequenced genomes. In addition, efforts in structural genomics resulted in a vastly improved structural coverage of proteomes, providing for definitive assessment of (non)homologous relationships between proteins.

RESULTS

We report the results of a comprehensive search for non-homologous isofunctional enzymes (NISE) that yielded 185 EC nodes with two or more experimentally characterized - or predicted - structurally unrelated proteins. Of these NISE sets, only 74 were from the original 1998 list. Structural assignments of the NISE show over-representation of proteins with the TIM barrel fold and the nucleotide-binding Rossmann fold. From the functional perspective, the set of NISE is enriched in hydrolases, particularly carbohydrate hydrolases, and in enzymes involved in defense against oxidative stress.

CONCLUSIONS

These results indicate that at least some of the non-homologous isofunctional enzymes were recruited relatively recently from enzyme families that are active against related substrates and are sufficiently flexible to accommodate changes in substrate specificity.

Structure of Trypanosoma brucei glutathione synthetase: domain and loop alterations in the catalytic cycle of a highly conserved enzyme

Glutathione synthetase catalyses the synthesis of the low molecular mass thiol glutathione from l-gamma-glutamyl-l-cysteine and glycine. We report the crystal structure of the dimeric enzyme from Trypanosoma brucei in complex with the product glutathione. The enzyme belongs to the ATP-grasp family, a group of enzymes known to undergo conformational changes upon ligand binding. The T. brucei enzyme crystal structure presents two dimers in the asymmetric unit. The structure reveals variability in the order and position of a small domain, which forms a lid for the active site and serves to capture conformations likely to exist during the catalytic cycle. Comparisons with orthologous enzymes, in particular from Homo sapiens and Saccharomyces cerevisae, indicate a high degree of sequence and structure conservation in part of the active site. Structural differences that are observed between the orthologous enzymes are assigned to different ligand binding states since key residues are conserved. This suggests that the molecular determinants of ligand recognition and reactivity are highly conserved across species. We conclude that it would be difficult to target the parasite enzyme in preference to the host enzyme and therefore glutathione synthetase may not be a suitable target for antiparasitic drug discovery.

Structural basis for feedback and pharmacological inhibition of Saccharomyces cerevisiae glutamate cysteine ligase

Structural characterization of glutamate cysteine ligase (GCL), the enzyme that catalyzes the initial, rate-limiting step in glutathione biosynthesis, has revealed many of the molecular details of substrate recognition. To further delineate the mechanistic details of this critical enzyme, we have determined the structures of two inhibited forms of Saccharomyces cerevisiae GCL (ScGCL), which shares significant sequence identity with the human enzyme. In vivo, GCL activity is feedback regulated by glutathione. Examination of the structure of ScGCL-glutathione complex (2.5 A; R = 19.9%, R(free) = 25.1%) indicates that the inhibitor occupies both the glutamate- and the presumed cysteine-binding site and disrupts the previously observed Mg(2+) coordination in the ATP-binding site. l-Buthionine-S-sulfoximine (BSO) is a mechanism-based inhibitor of GCL and has been used extensively to deplete glutathione in cell culture and in vivo model systems. Inspection of the ScGCL-BSO structure (2.2 A; R = 18.1%, R(free) = 23.9%) confirms that BSO is phosphorylated on the sulfoximine nitrogen to generate the inhibitory species and reveals contacts that likely contribute to transition state stabilization. Overall, these structures advance our understanding of the molecular regulation of this critical enzyme and provide additional details of the catalytic mechanism of the enzyme.

Structural basis for evolution of product diversity in soybean glutathione biosynthesis

The redox active peptide glutathione is ubiquitous in nature, but some plants also synthesize glutathione analogs in response to environmental stresses. To understand the evolution of chemical diversity in the closely related enzymes homoglutathione synthetase (hGS) and glutathione synthetase (GS), we determined the structures of soybean (Glycine max) hGS in three states: apoenzyme, bound to gamma-glutamylcysteine (gammaEC), and with hGSH, ADP, and a sulfate ion bound in the active site. Domain movements and rearrangement of active site loops change the structure from an open active site form (apoenzyme and gammaEC complex) to a closed active site form (hGSH*ADP*SO(4)(2-) complex). The structure of hGS shows that two amino acid differences in an active site loop provide extra space to accommodate the longer beta-Ala moiety of hGSH in comparison to the glycinyl group of glutathione. Mutation of either Leu-487 or Pro-488 to an Ala improves catalytic efficiency using Gly, but a double mutation (L487A/P488A) is required to convert the substrate preference of hGS from beta-Ala to Gly. These structures, combined with site-directed mutagenesis, reveal the molecular changes that define the substrate preference of hGS, explain the product diversity within evolutionarily related GS-like enzymes, and reinforce the critical role of active site loops in the adaptation and diversification of enzyme function.

Comprehensive analysis of the helix-X-helix motif in soluble proteins

Alpha-helices are the most common secondary structures found in globular proteins. In this report, we analyze the stereochemical and sequence properties of helix-X-helix (HXH) motifs in which two alpha-helices are linked by a single residue, in search of characteristic structures and sequence signals. The analysis is carried out on a database of 837 nonredundant HXH motifs. The kinks are characterized by the bend angle between the axes of the N-terminal and C-terminal helices and the wobble angle corresponding to the rotation of C-terminal helix axis on the plane perpendicular to the N-terminal one. The phi-psi dihedral angles of the linker residue are clustered in six distinct areas of the Ramachandran plot: two areas are located in the additional allowed alpha region (alpha(1) and alpha(2)), two areas are in the additional allowed beta region (beta(1) and beta(2)) and two areas have positive phi values (alpha(L) and beta(M)). Each phi/psi region corresponds to characteristic bend and wobble angles and amino acid distributions. Bend angles can vary from 0 degrees to 160 degrees. Most wobble angles correspond to a counter-clockwise rotation of the C-terminal helix. Proline residues are rigorously excluded from the linker position X but have a high propensity at position X+1 of the beta(1) and beta(2) motifs (12 and 7, respectively) and at position X+3 of the alpha(1) motifs (9). Glycine linkers are located either in the alpha(L) region (20%) or in the beta(M) region (80%). This latter conformation is characterized by a marked bend angle (124 degrees +/- 18 degrees) and a clockwise wobble. Among other amino acids, Asn is remarkable for its high propensity (>3) at the linker position of the alpha(2), beta(1), and beta(2) motifs. Stabilization of HXH motifs by H-bonds between polar side chains of the linker and polar groups of the backbone is determined. A method based on position-specific scoring matrices is developed for conformational prediction. The accuracy of the predictions reaches 80% when the method is applied to proline-induced kinks or to kinks with bend angles in the 50 degrees-100 degrees range.

CPSARST: an efficient circular permutation search tool applied to the detection of novel protein structural relationships

CPSARST (Circular Permutation Search Aided by Ramachandran Sequential Transformation) is an efficient database search tool that provides a new way for rapidly detecting novel relationships among proteins.

Circular permutation of a protein can be visualized as if the original amino- and carboxyl termini were linked and new ones created elsewhere. It has been well-documented that circular permutants usually retain native structures and biological functions. Here we report CPSARST (Circular Permutation Search Aided by Ramachandran Sequential Transformation) to be an efficient database search tool. In this post-genomics era, when the amount of protein structural data is increasing exponentially, it provides a new way to rapidly detect novel relationships among proteins.

A comprehensive analysis of non-sequential alignments between all protein structures

Background

The majority of relations between proteins can be represented as a conventional sequential alignment. Nevertheless, unusual non-sequential alignments with different connectivity of the aligned fragments in compared proteins have been reported by many researchers. It is interesting to understand those non-sequential alignments; are they unique, sporadic cases or they occur frequently; do they belong to a few specific folds or spread among many different folds, as a common feature of protein structure. We present here a comprehensive large-scale study of non-sequential alignments between available protein structures in Protein Data Bank.

Results

The study has been conducted on a non-redundant set of 8,865 protein structures aligned with the aid of the TOPOFIT method. It has been estimated that between 17.4% and 35.2% of all alignments are non-sequential depending on variations in the parameters. Analysis of the data revealed that non-sequential relations between proteins do occur systematically and in large quantities. Various sizes and numbers of non-sequential fragments have been observed with all possible complexities of fragment rearrangements found for alignments consisting of up to 12 fragments. It has been found that non-sequential alignments are not limited to proteins of any particular fold and are present in more than two hundred of them. Moreover, many of them are found between proteins with different fold assignments. It has been shown that protein structure symmetry does not explain non-sequential alignments. Therefore, compelling evidences have been provided that non-sequential alignments between proteins are systematic and widespread across the protein universe.

Conclusion

The phenomenon of the widespread occurrence of non-sequential alignments between proteins might represent a missing rule of protein structure organization. More detailed study of this phenomenon will enhance our understanding of protein stability, folding, and evolution.

Reaction mechanism of glutathione synthetase from Arabidopsis thaliana: site-directed mutagenesis of active site residues

Glutathione is essential for maintaining the intracellular redox environment and is synthesized from gamma-glutamylcysteine, glycine, and ATP by glutathione synthetase (GS). To examine the reaction mechanism of a eukaryotic GS, 24 Arabidopsis thaliana GS (AtGS) mutants were kinetically characterized. Within the gamma-glutamylcysteine/glutathione-binding site, the S153A and S155A mutants displayed less than 4-fold changes in kinetic parameters with mutations of Glu-220 (E220A/E220Q), Gln-226 (Q226A/Q226N), and Arg-274 (R274A/R274K) at the distal end of the binding site resulting in 24-180-fold increases in the K(m) values for gamma-glutamylcysteine. Substitution of multiple residues interacting with ATP (K313M, K367M, and E429A/E429Q) or coordinating magnesium ions to ATP (E148A/E148Q, N150A/N150D, and E371A) yielded inactive protein because of compromised nucleotide binding, as determined by fluorescence titration. Other mutations in the ATP-binding site (E371Q, N376A, and K456M) resulted in greater than 30-fold decreases in affinity for ATP and up to 80-fold reductions in turnover rate. Mutation of Arg-132 and Arg-454, which are positioned at the interface of the two substrate-binding sites, affected the enzymatic activity differently. The R132A mutant was inactive, and the R132K mutant decreased k(cat) by 200-fold; however, both mutants bound ATP with K(d) values similar to wild-type enzyme. Minimal changes in kinetic parameters were observed with the R454K mutant, but the R454A mutant displayed a 160-fold decrease in k(cat). In addition, the R132K, R454A, and R454K mutations elevated the K(m) value for glycine up to 11-fold. Comparison of the pH profiles and the solvent deuterium isotope effects of A. thaliana GS and the Arg-132 and Arg-454 mutants also suggest distinct mechanistic roles for these residues. Based on these results, a catalytic mechanism for the eukaryotic GS is proposed.

Inborn errors in the metabolism of glutathione

Glutathione is a tripeptide composed of glutamate, cysteine and glycine. Glutathione is present in millimolar concentrations in most mammalian cells and it is involved in several fundamental biological functions, including free radical scavenging, detoxification of xenobiotics and carcinogens, redox reactions, biosynthesis of DNA, proteins and leukotrienes, as well as neurotransmission/neuromodulation. Glutathione is metabolised via the gamma-glutamyl cycle, which is catalyzed by six enzymes. In man, hereditary deficiencies have been found in five of the six enzymes. Glutathione synthetase deficiency is the most frequently recognized disorder and, in its severe form, it is associated with hemolytic anemia, metabolic acidosis, 5-oxoprolinuria, central nervous system (CNS) damage and recurrent bacterial infections. Gamma-glutamylcysteine synthetase deficiency is also associated with hemolytic anemia, and some patients with this disorder show defects of neuromuscular function and generalized aminoaciduria. Gamma-glutamyl transpeptidase deficiency has been found in patients with CNS involvement and glutathionuria. 5-Oxoprolinase deficiency is associated with 5-oxoprolinuria but without a clear association with other symptoms. Dipeptidase deficiency has been described in one patient. All disorders are very rare and inherited in an autosomal recessive manner. Most of the mutations are leaky so that many patients have residual enzyme activity. Diagnosis is made by measuring the concentration of different metabolites in the gamma-glutamyl cycle, enzyme activity and in glutathione synthetase and gamma-glutamylcysteine synthetase deficiency, also by mutation analysis. Prenatal diagnosis has been preformed in glutathione synthetase deficiency. The prognosis is difficult to predict, as few patients are known, but seems to vary significantly between different patients. The aims of the treatment of glutathione synthesis defects are to avoid hemolytic crises and to increase the defense against reactive oxygen species. No treatment has been recommended for gamma-glutamyl transpeptidase, 5-oxoprolinase and dipeptidase deficiency.

Local comparison of protein structures highlights cases of convergent evolution in analogous functional sites

Background

We performed an exhaustive search for local structural similarities in an ensemble of non-redundant protein functional sites. With the purpose of finding new examples of convergent evolution, we selected only those matching sites composed of structural regions whose residue order is inverted in the relative protein sequences.

Results

A novel case of local analogy was detected between members of the ABC transporter and of the HprK/P families in their ATP binding site. This case cannot be derived by events of circular permutation since the residues of one of the region pairs are located in reverse order in the sequence of the two protein families. One of the analogous binding sites, the one identified in HprK/P, is known to also bind pyrophosphate, which is used as preferred energy source in its kinase and phosphorylase activity.

Conclusion

The discovery of this striking molecular similarity, also associated to a functional similarity, may help in suggesting new experiments aimed at a deeper understanding of members of the ABC transporter family known to be involved in many serious human diseases.

Systematic comparison of catalytic mechanisms of hydrolysis and transfer reactions classified in the EzCatDB database

Catalytic mechanisms of 270 enzymes from 131 superfamilies, mainly hydrolases and transferases, were analyzed based on their enzyme structures. A method of systematic comparison and classification of the catalytic reactions was developed. Hydrolysis and transfer reactions closely resemble one another, displaying common mechanisms, single displacement, and double displacement. These displacement mechanisms might be further subclassified according to the type of catalytic factors and nucleophilic substitution involved. Several types of catalytic factors exist: nucleophile, acid, base, stabilizer, modulator, cofactors. Nucleophilic substitution might be categorized as S(N)1/S(N)2 (or dissociative/associative) reactions. The classification indicates that some mechanisms favor particular types of catalytic factors. In hydrolyses of amide bonds and phosphoric ester bonds, mechanisms with single displacement tend to use inorganic cofactors such as zinc and magnesium ions as important catalysts, whereas those with double displacement frequently do not use such cofactors. In contrast, hydrolyses of O-glycoside bond rarely use such cofactors, with one exception. The trypsin-like hydrolytic reaction, which is catalyzed by the classic catalytic triad comprising serine/histidine/aspartate, can be considered as a "super-reaction" because it is observed in at least three nonhomologous enzymes, whereas most reactions are singlets without any nonhomologous enzymes. By dividing complex reactions into several reactions, correlations between active site structures and catalytic functions can be suggested. This classification method is applicable to other reactions such as elimination and isomerization. Furthermore, it will facilitate annotation of enzyme functions from 3D patterns of enzyme active sites. The classification is available at http://mbs.cbrc.jp/EzCatDB/RLCP/index.html.

Catalytic loop motion in human glutathione synthetase: A molecular modeling approach

Conformational changes of three flexible loops (G, A, and S) in human glutathione synthetase (hGS) arise to accommodate the substrates inside the active site. The crystal structure of hGS, a member of the ATP-grasp superfamily, has been reported only for the product-enzyme complex. To study the function of the hGS loops, molecular dynamics simulations are performed on three different conformational models: unbound enzyme, reactant-enzyme, and product-enzyme complex of hGS. The conformational changes among the three models are analyzed and the roles of the loops during the catalytic process are described. The modeled structures of hGS show that the central portions of the G- and A-loop have a double role in the reactant complex conformation: they bind the substrates and simultaneously interact with each other through an extensive network of hydrogen bonds. The present study proposes that these favorable loop-ligand and loop-loop interactions are required for opening and closing of the active site of hGS. Additionally, this research identifies important amino acid residues and explains their function within the catalytic loops of hGS.