Structure of the novel monomeric glyoxalase I from Zea mays

Linking the key physiological functions of essential micronutrients to their deficiency symptoms in plants

In this review, we untangle the physiological key functions of the essential micronutrients and link them to the deficiency responses in plants. Knowledge of these responses at the mechanistic level, and the resulting deficiency symptoms, have improved over the last decade and it appears timely to review recent insights for each of them. A proper understanding of the links between function and symptom is indispensable for an accurate and timely identification of nutritional disorders, thereby informing the design and development of sustainable fertilization strategies. Similarly, improved knowledge of the molecular and physiological functions of micronutrients will be important for breeding programmes aiming to develop new crop genotypes with improved nutrient-use efficiency and resilience in the face of changing soil and climate conditions.

Metal-Binding Protein TaGlo1 Improves Fungal Resistance to Arsenite (AsIII) and Methylarsenite (MAsIII) in Paddy Soil

Trivalent arsenicals such as arsenite (AsIII) and methylarsenite (MAsIII) are thought to be ubiquitous in flooded paddy soils and have higher toxicity than pentavalent forms. Fungi are widely prevalent in the rice rhizosphere, and the latter is considered a hotspot for As uptake. However, few studies have focused on alleviating As toxicity in paddy soils using fungi. In this study, we investigated the mechanism by which the protein TaGlo1, derived from the As-resistant fungal strain Trichoderma asperellum SM-12F1, mitigates AsIII and MAsIII toxicity in paddy soils. Taglo1 gene expression in Escherichia coli BL21 conferred strong resistance to AsIII and MAsIII, while purified TaGlo1 showed a high affinity for AsIII and MAsIII. Three cysteine residues (Cys13, Cys18, and Cys71) play crucial roles in binding with AsIII, while only two (Cys13 and Cys18) play crucial roles for MAsIII binding. TaGlo1 had a stronger binding strength for MAsIII than AsIII. Importantly, up to 90.2% of the homologous TaGlo1 proteins originate from fungi by GenBank searching. In the rhizospheres of 14 Chinese paddy soils, Taglo1 was widely distributed and its gene abundance increased with porewater As. This study highlights the potential of fungi to mitigate As toxicity and availability in the soil-rice continuum and suggests future microbial strategies for bioremediation.

Viridiplantae-specific GLXI and GLXII isoforms co-evolved and detoxify glucosone in planta

Reactive carbonyl species (RCS) such as methylglyoxal (MGO) and glyoxal (GO) are highly reactive, unwanted side-products of cellular metabolism maintained at harmless intracellular levels by specific scavenging mechanisms.MGO and GO are metabolized through the glyoxalase (GLX) system, which consists of two enzymes acting in sequence, GLXI and GLXII. While plant genomes encode a number of different GLX isoforms, their specific functions and how they arose during evolution are unclear. Here, we used Arabidopsis (Arabidopsis thaliana) as a model species to investigate the evolutionary history of GLXI and GLXII in plants and whether the GLX system can protect plant cells from the toxicity of RCS other than MGO and GO. We show that plants possess two GLX systems of different evolutionary origins and with distinct structural and functional properties. The first system is shared by all eukaryotes, scavenges MGO and GO, especially during seedling establishment, and features Zn2+-type GLXI proteins with a metal cofactor preference that were present in the last eukaryotic common ancestor. GLXI and GLXII of the second system, featuring Ni2+-type GLXI, were acquired by the last common ancestor of Viridiplantae through horizontal gene transfer from proteobacteria and can together metabolize keto-D-glucose (KDG, glucosone), a glucose-derived RCS, to D-gluconate. When plants displaying loss-of-function of a Viridiplantae-specific GLXI were grown in KDG, D-gluconate levels were reduced to 10%-15% of those in the wild type, while KDG levels showed an increase of 48%-67%. In contrast to bacterial GLXI homologs, which are active as dimers, plant Ni2+-type GLXI proteins contain a domain duplication, are active as monomers, and have a modified second active site. The acquisition and neofunctionalization of a structurally, biochemically, and functionally distinct GLX system indicates that Viridiplantae are under strong selection to detoxify diverse RCS.

The Human Glyoxalase Gene Family in Health and Disease

The glyoxalase gene family consists of six structurally and functionally diverse enzymes with broad roles in metabolism. The common feature that defines this family is based on structural motifs that coordinate divalent cations which are required for activity. These family members have been implicated in a variety of physiological processes, including amino-acid metabolism (4-hydroxyphenylpyruvate dioxygenase; HPD), primary metabolism (methylmalonyl-CoA epimerase; MCEE), and aldehyde detoxication (glyoxalase 1; GLO1) and therefore have significant associations with disease. A central function of this family is the detoxification of reactive dicarbonyls (e.g., methylglyoxal), which react with cellular nucleophiles, resulting in the modification of lipids, proteins, and DNA. These damaging modifications activate canonical stress responses such as heat shock, unfolded protein, antioxidant, and DNA damage responses. Thus, glyoxalases serve an important role in homeostasis, preventing the pathogenesis of metabolic disease states, including obesity, diabetes, cardiovascular disease, renal failure, and aging. This review presents a thorough overview of the literature surrounding this diverse enzyme class. Although extensive literature exists for some members of this family (e.g., GLO1), little is known about the physiological role of glyoxalase domain-containing protein 4 (GLOD4) and 5 (GLOD5), paving the way for exciting avenues for future research.

Genome-Wide Expression Analysis of Glyoxalase I Genes Under Hyperosmotic Stress and Existence of a Stress-Responsive Mitochondrial Glyoxalase I Activity in Durum Wheat (Triticum durum Desf.)

Glyoxalase I (GLYI) catalyzes the rate-limiting step of the glyoxalase pathway that, in the presence of GSH, detoxifies the cytotoxic molecule methylglyoxal (MG) into the non-toxic D-lactate. In plants, MG levels rise under various abiotic stresses, so GLYI may play a crucial role in providing stress tolerance. In this study, a comprehensive genome database analysis was performed in durum wheat (Triticum durum Desf.), identifying 27 candidate GLYI genes (TdGLYI). However, further analyses of phylogenetic relationships and conserved GLYI binding sites indicated that only nine genes encode for putative functionally active TdGLYI enzymes, whose distribution was predicted in three different subcellular compartments, namely cytoplasm, plastids and mitochondria. Expression profile by qRT-PCR analysis revealed that most of the putative active TdGLYI genes were up-regulated by salt and osmotic stress in roots and shoots from 4-day-old seedlings, although a different behavior was observed between the two types of stress and tissue. Accordingly, in the same tissues, hyperosmotic stress induced an increase (up to about 40%) of both GLYI activity and MG content as well as a decrease of GSH (up to about -60%) and an increase of GSSG content (up to about 7-fold) with a consequent strong decrease of the GSH/GSSG ratio (up to about -95%). Interestingly, in this study, we reported the first demonstration of the existence of GLYI activity in highly purified mitochondrial fraction. In particular, GLYI activity was measured in mitochondria from durum wheat (DWM), showing hyperbolic kinetics with Km and Vmax values equal to 92 ± 0.2 μM and 0.519 ± 0.004 μmol min-1 mg-1 of proteins, respectively. DWM-GLYI resulted inhibited in a competitive manner by GSH (Ki = 6.5 ± 0.7 mM), activated by Zn2+ and increased, up to about 35 and 55%, under salt and osmotic stress, respectively. In the whole, this study provides basis about the physiological significance of GLYI in durum wheat, by highlighting the role of this enzyme in the early response of seedlings to hyperosmotic stress. Finally, our results strongly suggest the existence of a complete mitochondrial GLYI pathway in durum wheat actively involved in MG detoxification under hyperosmotic stress.

Two-Substrate Glyoxalase I Mechanism: A Quantum Mechanics/Molecular Mechanics Study

Glyoxalase I (GlxI) is an important enzyme that catalyzes the detoxification of methylglyoxal (MG) with the help of glutathione (H-SG). It is currently unclear whether MG and H-SG are substrates of GlxI or whether the enzyme processes hemithioacetal (HTA), which is nonenzymatically formed from MG and H-SG. Most previous studies have concentrated on the latter mechanism. Here, we study the two-substrate reaction mechanism of GlxI from humans (HuGlxI) and corn (ZmGlxI), which are Zn(II)-active and -inactive, respectively. Hybrid quantum mechanics/molecular mechanics calculations were used to obtain geometrical structures of the stationary points along reaction paths, and big quantum mechanical systems with more than 1000 atoms and free-energy perturbations were used to improve the quality of the calculated energies. We studied, on an equal footing, all reasonable reaction paths to the S- and R-enantiomers of HTA from MG and H-SG (the latter was considered in two different binding modes). The results indicate that the MG and H-SG reaction in both enzymes can follow the same path to reach S-HTA. However, the respective overall barriers and reaction energies are different for the two enzymes (6.1 and -9.8 kcal/mol for HuGlxI and 15.7 and -2.2 kcal/mol for ZmGlxI). The first reaction step to produce S-HTA is facilitated by a crystal water molecule that forms hydrogen bonds with a Glu and a Thr residue in the active site. The two enzymes also follow similar paths to R-HTA. However, the reactions reach a deprotonated and protonated R-HTA in the human and corn enzymes, respectively. The production of deprotonated R-HTA in HuGlxI is consistent with other theoretical and experimental works. However, our calculations show a different behavior for ZmGlxI (both S- and R-HTA can be formed in the enzyme with the alcoholic proton on HTA). This implies that Glu-144 of corn GlxI is not basic enough to keep the alcoholic proton. In HuGlxI, the two binding modes of H-SG that lead to S- and R-HTA are degenerate, but the barrier leading to R-HTA is lower than the barrier to S-HTA. On the other hand, ZmGlxI prefers the binding mode, which produces S-HTA; this observation is consistent with experiments. Based on the results, we present a modification for a previously proposed two-substrate reaction mechanism for ZmGlxI.

Quantum Mechanics/Molecular Mechanics Study of the Reaction Mechanism of Glyoxalase I

Glyoxalase I (GlxI) is a member of the glyoxalase system, which is important in cell detoxification and converts hemithioacetals of methylglyoxal (a cytotoxic byproduct of sugar metabolism that may react with DNA or proteins and introduce nucleic acid strand breaks, elevated mutation frequencies, and structural or functional changes of the proteins) and glutathione into d-lactate. GlxI accepts both the S and R enantiomers of hemithioacetal, but converts them to only the S-d enantiomer of lactoylglutathione. Interestingly, the enzyme shows this unusual specificity with a rather symmetric active site (a Zn ion coordinated to two glutamate residues; Glu-99 and Glu-172), making the investigation of its reaction mechanism challenging. Herein, we have performed a series of combined quantum mechanics and molecular mechanics calculations to study the reaction mechanism of GlxI. The substrate can bind to the enzyme in two different modes, depending on the direction of its alcoholic proton (H2; toward Glu-99 or Glu-172). Our results show that the S substrate can react only if H2 is directed toward Glu-99 and the R substrate only if H2 is directed toward Glu-172. In both cases, the reactions lead to the experimentally observed S-d enantiomer of the product. In addition, the results do not show any low-energy paths to the wrong enantiomer of the product from neither the S nor the R substrate. Previous studies have presented several opposing mechanisms for the conversion of R and S enantiomers of the substrate to the correct enantiomer of the product. Our results confirm one of them for the S substrate, but propose a new one for the R substrate.

Alr2321, a multiple stress inducible glyoxalase I of Anabaena sp. PCC7120 detoxifies methylglyoxal and reactive species oxygen

Abiotic stresses enhance the cellular level of reactive oxygen species (ROS) which consequently leads to toxic methylglyoxal (MG) production. Glyoxalases (GlyI & GlyII) catalyze the conversion of toxic MG into non-toxic lactic acid but their properties and functions have been overlooked in cyanobacteria. This is the first attempt to conduct a genome-wide analysis of GlyI protein (PF00903) from Anabaena sp. PCC7120. Out of total nine GlyI domain possessing proteins, only three (Alr2321, Alr4469, All1022) harbour conserve His/Glu/His/Glu metal binding site at their homologous position and are deficient in conserved region specific for Zn²⁺ dependent members. Their biochemical, structural and functional characterization revealed that only Alr2321 is a homodimeric Ni²⁺ dependent active GlyI with catalytic efficiency 11.7 × 10⁶ M^-1 s^-1. It has also been found that Alr2321 is activated by various divalent metal ions and has maximum GlyI activity with Ni²⁺ followed by Co²⁺ > Mn²⁺ > Cu²⁺ and no activity with Zn^2+. Moreover, the expression of alr2321 was found to be maximally up-regulated under heat (19 fold) followed by cadmium, desiccation, arsenic, salinity and UV-B stresses. BL21/pGEX-5X2-alr2321 showed improved growth under various abiotic stresses as compared to BL21/pGEX-5X2 by increased scavenging of intracellular MG and ROS levels. Taken together, these results suggest noteworthy links between intracellular MG and ROS, its detoxification by Alr2321, a member of GlyI family of Anabaena sp. PCC7120, in relation to abiotic stress.

Preliminary Characterization of a Ni2+-Activated and Mycothiol-Dependent Glyoxalase I Enzyme from Streptomyces coelicolor

The glyoxalase system consists of two enzymes, glyoxalase I (Glo1) and glyoxalase II (Glo2), and converts a hemithioacetal substrate formed between a cytotoxic alpha-ketoaldehyde, such as methylglyoxal (MG), and an intracellular thiol, such as glutathione, to a non-toxic alpha-hydroxy acid, such as d-lactate, and the regenerated thiol. Two classes of Glo1 have been identified. The first is a Zn2+-activated class and is exemplified by the Homo sapiens Glo1. The second class is a Ni2+-activated enzyme and is exemplified by the Escherichia coli Glo1. Glutathione is the intracellular thiol employed by Glo1 from both these sources. However, many organisms employ other intracellular thiols. These include trypanothione, bacillithiol, and mycothiol. The trypanothione-dependent Glo1 from Leishmania major has been shown to be Ni2+-activated. Genetic studies on Bacillus subtilis and Corynebacterium glutamicum focused on MG resistance have indicated the likely existence of Glo1 enzymes employing bacillithiol or mycothiol respectively, although no protein characterizations have been reported. The current investigation provides a preliminary characterization of an isolated mycothiol-dependent Glo1 from Streptomyces coelicolor. The enzyme has been determined to display a Ni2+-activation profile and indicates that Ni2+-activated Glo1 are indeed widespread in nature regardless of the intracellular thiol employed by an organism.

Deciphering the number and location of active sites in the monomeric glyoxalase I of Zea mays

Detoxification of methylglyoxal, a toxic by-product of central sugar metabolism, is a major issue for all forms of life. The glyoxalase pathway evolved to effectively convert methylglyoxal into d-lactate via a glutathione hemithioacetal intermediate. Recently, we have shown that the monomeric glyoxalase I from maize exhibits a symmetric fold with two cavities, potentially harboring two active sites, in analogy with homodimeric enzyme surrogates. Here we confirm that only one of the two cavities exhibits glyoxalase I activity and show that it adopts a tunnel-shaped structure upon substrate binding. Such conformational change gives rise to independent binding sites for glutathione and methylglyoxal in the same active site, with important implications for the molecular reaction mechanism, which has been a matter of debate for several decades. DATABASE: Structural data are available in The Protein Data Bank database under the accession numbers 6BNN, 6BNX, and 6BNZ.

Comparative analysis of glyoxalase pathway genes in Erianthus arundinaceus and commercial sugarcane hybrid under salinity and drought conditions

Background

Glyoxalase pathway is a reactive carbonyl species (RCS) scavenging mechanism involved in the detoxification of methylglyoxal (MG), which is a reactive α-ketoaldehyde. In plants under abiotic stress, the cellular toxicity is reduced through glyoxalase pathway genes, i.e. Glyoxalase I (Gly I), Glyoxalase II (Gly II) and Glyoxalase III (Gly III). Salinity and water deficit stresses produce higher amounts of endogenous MG resulting in severe tissue damage. Thus, characterizing glyoxalase pathway genes that govern the MG metabolism should provide new insights on abiotic stress tolerance in Erianthus arundinaceus, a wild relative of sugarcane and commercial sugarcane hybrid (Co 86032).

Results

In this study, three glyoxalase genes (Glyoxalase I, II and III) from E. arundinaceus (a wild relative of sugarcane) and commercial sugarcane hybrid (Co 86032) were characterized. Comparative gene expression profiles (qRT-PCR) of Glyoxalase I, II and III under salinity and water deficit stress conditions revealed differential transcript expression with higher levels of Glyoxalase III in both the stress conditions. Significantly, E. arundinaceus had a higher expression level of glyoxalase genes compared to commercial sugarcane hybrid. On the other hand, gas exchange parameters like stomatal conductance and transpiration rate were declined to very low levels under both salt and drought induced stresses in commercial sugarcane hybrid when compared to E. arundinaceus. E. arundinaceus maintained better net photosynthetic rate compared to commercial sugarcane hybrid. The phylogenetic analysis of glyoxalase proteins showed its close evolutionary relationship with Sorghum bicolor and Zea mays. Glyoxalase I and II were predicted to possess 9 and 7 isoforms respectively whereas, Glyoxalase III couldn’t be identified as it comes under uncharacterized protein identified in recent past. Chromosomal mapping is also carried out for glyoxalase pathway genes and its isoforms. Docking studies revealed the binding affinities of glyoxalase proteins in both E. arundinaceus and commercial sugarcane hybrid with their substrate molecules.

Conclusions

This study emphasizes the role of Glyoxalase pathway genes in stress defensive mechanism which route to benefit in progressive plant adaptations and serves as potential candidates for development of salt and drought tolerant crops.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5349-7) contains supplementary material, which is available to authorized users.

Dissecting the Physiological Function of Plant Glyoxalase I and Glyoxalase I-Like Proteins

The Arabidopsis genome annotation include 11 glyoxalase I (GLXI) genes, all encoding for protein members of the vicinal oxygen chelate (VOC) superfamily. The biochemical properties and physiological importance of three Arabidopsis GLXI proteins in the detoxification of reactive carbonyl species has been recently described. Analyses of phylogenetic relationships and conserved GLXI binding sites indicate that the other eight GLXI genes (GLXI-like) do not encode for proteins with GLXI activity. In this perspective article we analyse the structural features of GLXI and GLXI-like proteins, and explore splice forms and transcript abundance under abiotic stress conditions. Finally, we discuss future directions of research on this topic with respect to the substrate identification of GLXI and GLXI-like proteins and the need of reliable quantitative measurements of reactive carbonyl species in plant tissues.

Defense against Reactive Carbonyl Species Involves at Least Three Subcellular Compartments Where Individual Components of the System Respond to Cellular Sugar Status

Methylglyoxal (MGO) and glyoxal (GO) are toxic reactive carbonyl species generated as by-products of glycolysis. The pre-emption pathway for detoxification of these products, the glyoxalase (GLX) system, involves two consecutive reactions catalyzed by GLXI and GLXII. In Arabidopsis thaliana, the GLX system is encoded by three homologs of GLXI and three homologs of GLXII, from which several predicted GLXI and GLXII isoforms can be derived through alternative splicing. We identified the physiologically relevant splice forms using sequencing data and demonstrated that the resulting isoforms have different subcellular localizations. All three GLXI homologs are functional in vivo, as they complemented a yeast GLXI loss-of-function mutant. Efficient MGO and GO detoxification can be controlled by a switch in metal cofactor usage. MGO formation is closely connected to the flux through glycolysis and through the Calvin Benson cycle; accordingly, expression analysis indicated that GLXI is transcriptionally regulated by endogenous sugar levels. Analyses of Arabidopsis loss-of-function lines revealed that the elimination of toxic reactive carbonyl species during germination and seedling establishment depends on the activity of the cytosolic GLXI;3 isoform. The Arabidopsis GLX system involves the cytosol, chloroplasts, and mitochondria, which harbor individual components that might be used at specific developmental stages and respond differentially to cellular sugar status.

Characteristic Variations and Similarities in Biochemical, Molecular, and Functional Properties of Glyoxalases across Prokaryotes and Eukaryotes

The glyoxalase system is the ubiquitous pathway for the detoxification of methylglyoxal (MG) in the biological systems. It comprises two enzymes, glyoxalase I (GLYI) and glyoxalase II (GLYII), which act sequentially to convert MG into d-lactate, thereby helping living systems get rid of this otherwise cytotoxic byproduct of metabolism. In addition, a glutathione-independent GLYIII enzyme activity also exists in the biological systems that can directly convert MG to d-lactate. Humans and Escherichia coli possess a single copy of GLYI (encoding either the Ni- or Zn-dependent form) and GLYII genes, which through MG detoxification provide protection against various pathological and disease conditions. By contrast, the plant genome possesses multiple GLYI and GLYII genes with a role in abiotic stress tolerance. Plants possess both Ni²⁺- and Zn²⁺-dependent forms of GLYI, and studies on plant glyoxalases reveal the various unique features of these enzymes distinguishing them from prokaryotic and other eukaryotic glyoxalases. Through this review, we provide an overview of the plant glyoxalase family along with a comparative analysis of glyoxalases across various species, highlighting similarities as well as differences in the biochemical, molecular, and physiological properties of these enzymes. We believe that the evolution of multiple glyoxalases isoforms in plants is an important component of their robust defense strategies.

Genome-wide analysis and expression profiling of glyoxalase gene families in soybean (Glycine max) indicate their development and abiotic stress specific response

BACKGROUND

Glyoxalase pathway consists of two enzymes, glyoxalase I (GLYI) and glyoxalase II (GLYII) which detoxifies a highly cytotoxic metabolite methylglyoxal (MG) to its non-toxic form. MG may form advanced glycation end products with various cellular macro-molecules such as proteins, DNA and RNA; that ultimately lead to their inactivation. Role of glyoxalase enzymes has been extensively investigated in various plant species which showed their crucial role in salinity, drought and heavy metal stress tolerance. Previously genome-wide analysis of glyoxalase genes has been conducted in model plants Arabidopsis and rice, but no such study was performed in any legume species.

RESULTS

In the present study, a comprehensive genome database analysis of soybean was performed and identified a total of putative 41 GLYI and 23 GLYII proteins encoded by 24 and 12 genes, respectively. Detailed analysis of these identified members was conducted including their nomenclature and classification, chromosomal distribution and duplication, exon-intron organization, and protein domain(s) and motifs identification. Expression profiling of these genes has been performed in different tissues and developmental stages as well as under salinity and drought stresses using publicly available RNAseq and microarray data. The study revealed that GmGLYI-7 and GmGLYII-8 have been expressed intensively in all the developmental stages and tissues; while GmGLYI-6, GmGLYI-9, GmGLYI-20, GmGLYII-5 and GmGLYII-10 were highly abiotic stress responsive members.

CONCLUSIONS

The present study identifies the largest family of glyoxalase proteins to date with 41 GmGLYI and 23 GmGLYII members in soybean. Detailed analysis of GmGLYI and GmGLYII genes strongly indicates the genome-wide segmental and tandem duplication of the glyoxalase members. Moreover, this study provides a strong basis about the biological role and function of GmGLYI and GmGLYII members in soybean growth, development and stress physiology.