Potential role of the glyoxalase pathway as a drug target in Leishmania infantum: an exact steady-state model analysis

  • Nuno Filipe Lages1Email author,

    Affiliated with

    • Marta Sousa Silva1,

      Affiliated with

      • António EN Ferreira1,

        Affiliated with

        • Ana Maria Tomás2, 3,

          Affiliated with

          • Carlos Cordeiro1 and

            Affiliated with

            • Ana Ponces Freire1

              Affiliated with

              BMC Systems Biology20071(Suppl 1):S8

              DOI: 10.1186/1752-0509-1-S1-S8

              Published: 8 May 2007


              The glyoxalase system is the main catabolic pathway of methylglyoxal, a toxic 2-oxoaldehyde which occurs in all living cells as a by-product of glycolysis [1, 2]. Trypanosomatids, such as Leishmania infantum, are pathogenic microbial parasites which have glycolytic enzymes located in a specific cell organelle, the glycosome [3]. In these organisms, the usual glyoxalase cofactor glutathione is functionally replaced by trypanothione [4, 5]. Inhibition of the glyoxalase pathway might cause the accumulation of methylglyoxal in the glycosome, hampering glycolysis and killing the parasite.


              An ODE model of the glyoxalase pathway was developed based on a kinetic characterization of the glyoxalase system in L. infantum [6]. As depicted in Figure 1, it is assumed that methylglyoxal is formed at a constant rate, v in. Methylglyoxal reacts with one or both thiol groups of trypanothione forming mono- or bis-hemithioacetals in a reversible mass-action reaction. These hemithioacetals are the substrates of the glyoxalase system, forming D-lactate as final product and regenerating free trypanothione. v in is the input flux of methylglyoxal in the pathway; v 1 and v -1 are mass-action reaction rates; v 2 and v 3 are Michaelis-Menten enzyme-catalysed reaction rates. Steady-state concentrations and their sensitivities were derived as functions of glyoxalase I and II activities, methylglyoxal influx and total thiol concentration [6]. Necessary mathematical relations for the existence of a steady state were also calculated from the asymptotic behavior of these functions. All mathematical relations are exact and were derived symbolically using Mathematica® 4.0, from Wolfram Research, Inc.
              Figure 1

              Glyoxalase pathway of L. infantum.


              Methylglyoxal concentration achieves a steady state only when parameters remain within certain boundaries. Outside these boundaries, the typical dynamic behaviour of the system is the depletion of free trypanothione, leading to the accumulation of methylglyoxal over time. Total cofactor concentration cannot be arbitrarily low, it must exist above a defined threshold, the so-called "catalytic amount". Enzyme activities must also be greater than a threshold value, which is higher than the flux of the pathway. Conversely, for constant enzyme activities, there is a limit for methylglyoxal formation flux above which the glyoxalase pathway can not prevent the accumulation of methylglyoxal. Variations of enzyme activity have virtually no effect on methylglyoxal steady-state concentration (sensitivity of approximately 0). The non-enzymatic hemithioacetal formation has a sensitivity of -1. Methylglyoxal concentration sensitivity to input flux and total cofactor concentration are approximately 1 and -1, respectively. In the vicinity of the boundaries outside of which a steady state is not attained, the absolute values of all sensitivities rise to infinity. However, the reference parameter values, which correspond to the physiological state of L. infantum in exponential phase, are very far from the boundaries – compared to the minimum values required to attain a steady state, total cofactor is approximately 50,000 fold, glyoxalase I activity is approximately 23,000 fold, glyoxalase II activity is approximately 30,000 fold and input flux is approximately 25,000 fold.
              Figure 2

              Methylglyoxal steady-state concentration as a function of (a) total trypanothione concentration, (b) methylglyoxal formation rate, (c) glyoxalase I activity and (d) glyoxalase II activity. Parameters are normalized to reference values.


              Decreasing the activity of the glyoxalase enzymes is unsuitable to harm the parasite since no significant rise of methylglyoxal steady-state concentration is achieved. However, depletion of the cofactor trypanothione or the enhancement of methylglyoxal production leads to a significant increase of the concentration of this toxic compound. These results point to a possible strategy for treatment of leishmaniasis based on the interference with trypanothione metabolism or glycolysis. This study provides quantitative data which might be useful for future research on drug targeting this disease.

              Authors’ Affiliations

              Centro de Química e Bioquímica, Departamento de Química e Bioquímica, Faculdade de Ciências da Universidade de Lisboa
              Instituto de Ciências Biomédicas Abel Salazar, Universidade do Porto
              Instituto de Biologia Molecular e Celular, Universidade do Porto


              1. Thornalley PJ: The glyoxalase system: new developments towards functional characterization of a metabolic pathway fundamental to biological life. Biochem J 1990,269(1):1–11.PubMed
              2. Richard JP: Mechanism for the formation of methylglyoxal from triosephosphates. Biochem Soc Trans 1993,21(2):549–553.PubMed
              3. Hannaert V, Saavedra E, Duffieux F, Szikora JP, Rigden DJ, Michels PA, Opperdoes FR: Plant-like traits associated with metabolism of Trypanosoma parasites. Proc Natl Acad Sci USA 2003,100(3):1067–1071.View ArticlePubMed
              4. Muller S, Liebau E, Walter RD, Krauth-Siegel RL: Thiol-based redox metabolism of protozoan parasites. Trends Parasitol 2003,19(7):320–328.View ArticlePubMed
              5. Vickers TJ, Greig N, Fairlamb AH: A trypanothione-dependent glyoxalase I with a prokaryotic ancestry in Leishmania major. Proc Natl Acad Sci USA 2004,101(36):13186–13191.View ArticlePubMed
              6. Sousa Silva M, Ferreira AEN, Tomás AM, Cordeiro C, Ponces Freire A: Quantitative assessment of the glyoxalase pathway in Leishmania infantum as a therapeutic target by modelling and computer simulation. FEBS J 2005, 272:2388–2398.View ArticlePubMed


              © Lages et al 2007

              This article is published under license to BioMed Central Ltd.