# Coping with inevitable accidents in metabolism

Antoine Danchin

► To cite this version:

Antoine Danchin. Coping with inevitable accidents in metabolism. Microbial Biotechnology, 2017, 10 (1), pp.57-72. 10.1111/1751-7915.12461. hal-01496339

HAL Id: hal-01496339

https://hal.sorbonne-universite.fr/hal-01496339

Submitted on 27 Mar 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Genomics update

Coping with inevitable accidents in metabolism

Antoine Danchin*

Institute of Cardiometabolism and Nutrition, Hôpital de la Pitié-Salpêtrière, 47 Boulevard de l'Hôpital, Paris, 75013, France.

Summary

Genomic studies focus on key metabolites and pathways that, despite their obvious anthropocentric design, keep being ‘predicted’, while this is only finding again what is already known. As increasingly more genomes are sequenced, this lightpost effect may account at least in part for our failure to understand the function of a continuously growing number of genes. Core metabolism often goes astray, accidentally producing a variety of unexpected compounds. Catabolism of these forgotten metabolites makes an essential part of the functions coded in metagenomes. Here, I explore the fate of a limited number of those: compounds resulting from radical reactions and molecules derived from some reactive intermediates produced during normal metabolism. I try both to update investigators with the most recent literature and to uncover old articles that may open up new research avenues in the genome exploration of metabolism. This should allow us to foresee further developments in experimental genomics and genome annotation.

Introduction

Most of the time the cell is vibrant with chemical activity.
Dormancy – the condition of the spore or the seed – is a
state that lies in limbo, not dead but not quite alive. Life
has to go forward. It begins with the chemical transfor-
mations which make the root of metabolism. At the
ground level of life, we find dynamic transformations that

extract compounds from the environment, process them
into the cell's building blocks and energy supply, while
discarding what is toxic or cannot be used. Continuously
improved by genomic studies, a comprehensive descrip-
tive knowledge of the matching pathways is available in
textbooks and databases that provide access to all what
we would like to ask about metabolism (perhaps exclud-
ing what we failed to notice). This provides us with a
description of the basic chemical transformations that
continuously unfold in cells. They look smooth while cer-
tainly sophisticated. Yet, as in all dynamic processes,
especially those involving multiple stages, things must
sometimes go awry: even metabolites must be repaired
(Danchin et al., 2011). Despite their name, accidents are
the rule, not the exception. In a famous book, Charles
Perrow explored the inevitability of accidents, events that
are not foreseen individually, that may be separately
thought to happen extremely rarely, but that are doomed
to happen sometime, and that he named for that reason
'normal' accidents (Perrow, 1984). In the case of the
cell's metabolism, this is illustrated by the fact that mole-
cules will inevitably be produced or modified out of the
anticipated pathways. For example, metabolites become
accidentally oxidized or alkylated. Furthermore, there will
always be some shadow accompanying core metabo-
lism, with reactions going off the right track, resulting in
variations on the theme of normal metabolism, participat-
ing in what was named 'underground' metabolism (D'Ari
and Casadesus, 1998) and more recently 'paralogous'
metabolism (Chan et al., 2014). The latter term was pro-
posed to illustrate the likely involvement of enzymes that
were paralogues of those involved in normal metabolism
as a ready-made way to cope with metabolites that were
chemical variants of the normal compounds. Many pro-
cesses contribute to the way core metabolism goes
astray. Rather than follow the logic of metabolism (Dan-
chin and Sekowska, 2014), which is essentially based
on a combination of chemistry, physics and geology, I
explore here widely used genomic resources as well as
recent articles to try and identify trends in the way cells
cope with chemical accidents. Obviously, it is impossible
to explore all the chemical combinations that may have
some importance in the cell's life. Here, I propose some
tracks that I hope will be useful to genome-related
research, via identification of solutions found by living

Funding information

No funding information provided.

<footnote>*For correspondence. E-mail: antoine.danchin@normalesup.org;
Tel. +331 4217 7981; Fax +331 4217 8473.
*Microbial Biotechnology* (2017) **10**(1), 57–72
doi:10.1111/1751-7915.12461</footnote>
organisms to cope with 'normal' accidents during the course of evolution. In order to help future genome annotators to discover unsuspected functions, I split the structure of accidents into two main topics: those resulting from radicals and those resulting from other reactive intermediates (involved in ageing and senescence in particular). In the following, I illustrate the way bacteria cope with reactive molecules with a few concrete examples based on genome studies. The way I proceed may be used for further genome functional exploration, starting, for example, from the list of reactive molecules proposed by Hanson and co-workers (Lerma-Ortiz et al., 2016). I have tried, as much as possible, to go back in time and provide reference to early work that, to my knowledge, has not been followed up but might be enlightening for future genome annotation and genome-driven experiments.

## Unavoidable metabolic errors caused by free radicals

A first hint at possible accident-prone reactions is the presence of so-called spontaneous reactions (i.e. assumed to appear in the absence of an enzyme) in metabolism. The widely used reference MetaCyc knowledge base long collected such reactions. It lists at this time some 480 spontaneous reactions, basically identified by compiling articles that describe chemical pathways in the cell and require steps with apparently no associated protein. This compilation makes a fairly hap hazard, but revealing collection, based on the interests of biologists, not with much background rationale: http:// metacyc.org/META/search-query?type=REACTION&catal =spont (Caspi et al., 2016). As a consequence, a great many of those reactions correspond to quite rare com pounds. I will not explore them further here, but will gather from the catalogue some reactions that have a major importance. Among spontaneous reactions, those derived from free radicals are the most obvious reactive entities produced by the cell metabolism. A first line of radical-generating reactive species are produced in the environment, in relation with more or less random energy-rich physical processes (photons or high energy particles coming from a variety of irradiation sources such as sunlight, cosmic rays or intrinsic radioactivity). The basal level of radioactivity (in particular in the form of ⁴⁰K in all cells, as well as inescapable cosmic radiation) produces a low, but sustained level of free radicals. Also, chemical reactions involve electron transfers, and unpaired electrons are particularly reactive. While these processes have been the subject of a wealth of studies involving macromolecules and membranes, their effect on metabolites has seldom been explored. Furthermore, many reactions involve free radicals (the list provided in MetaCyc is sketchy, and far from complete), which are

so reactive that they will find a target much faster than what enzymes would help them to do, except when they are generated in situ, and this is expected to reveal a variety of accidental targets.

Free radicals will react with almost any biological compound, resulting in chemical modifications that are subsequently recovered in the metabolite complement of the cell. While natural environmental irradiation progressively went down as the Earth came of age, the emergence of dioxygen in the atmosphere considerably enhanced radical production via reactive oxygen species (ROSs), especially in environments allowing the presence of free reduced iron. The ubiquitous Fenton reaction is a general source of the hydroxyl radical that may attack almost any of the cell's components (Anglada et al., 2015). In parallel, related reactive nitrogen species (NO⁺ in particular (Pullan et al., 2008)) and reactive sulfur species (Giles and Jacob, 2002) contributed to the general free radical burden of the cell. It is important to note, though, that this should not be considered as always negative: for example, evolution recruited free radicals to be used by cells as means to compete invading organisms (Frawley and Fang, 2014). Hence, metabolism is generally replete with free radicals as well as a plethora of non-radical but highly reactive derivatives (singlet oxygen, hydrogen peroxide or hypochlorite, which easily generate radical ROSs) that will contribute to the panel of unconventional metabolites in the cell.

A first line of defence for the cell is to act before these reactive entities are generated, via destruction of some of their sources. To this aim, it uses specific enzymes (peroxiredoxins, catalase, superoxide dismutase, glutathione peroxidase, in particular). In the case of reactive nitrogen species, NO may enter the nitrogen cycle (Eady et al., 2016). In a second line of defence, cells evolved tolerance processes (Giuffre et al., 2014). Some organisms, such as Lactobacilli and Mycoplasmas, have alleviated the bulk of dioxygen-dependent radicals generation by dispensing of iron and using manganese as a ROS scavenger (Lisher and Giedroc, 2013). When annotating genome sequences, it is essential to have this countermeasure in mind as it will impact overall genome annotation (Danchin and Fang, 2016). A general shielding process (that I do not discuss here as this would ask for another full article) is for the cell to protect reactive groups such as amines or thiols using available metabolic resources. This might account for widespread methylation, acetylation or succinylation of amine, imidazole, guanidinium or thiol groups (Levitan et al., 2010; Lu et al., 2011). Some of these modifications would have subsequently be recruited for regulation (and phosphorylation appears mainly to be used as a regulatory tag) and a general exploration of protection/deprotection processes in the cell should be undertaken (see (Chan
et al., 2014)), in particular because the very protection
systems (S-adenosylmethionine in particular) may acci-
dentally modify unwanted targets. To escape unavoid-
able metabolic accidents, a general way out is for cells
to use the buffering capacity of some of its components,
as I now discuss.

Sulfur building blocks: cysteine and methionine

As an example of the way cells cope with these events,
I illustrate the situation created by the presence of the
sulfur atom, which is a particularly reactive nucleophile.
When sulfur-bearing targets have been modified, a repair
mechanism can operate, preventing accumulation of an
altered cell building block which might enter paralogous
metabolism. A typical example is provided by alkylation
of cysteines in proteins, which after proteolysis will gen-
erate S-alkyl-cysteines, analogues of proteogenic amino
acids that are therefore toxic and need to be put aside
(Chan et al., 2014). ROS will also react with cysteine
residues, which, via their chemical reactivity, manage
redox states, first by creating and disrupting disulfide
bonds. Subsequently, oxidation may go to the unstable
sulfenic (-S-OH, (Gupta and Carroll, 2014)), then sulfinic
(-S-O₂H) and finally sulfonic (-SO₃H) form of the amino
acid (van Bergen et al., 2014). The sulfenic form is fairly
easily reduced back to -SH via pathways involving glu-
tathione (or related molecules, such as N-acetyl-cysteine
or bacillithiol in Firmicutes and mycothiol in Actinobacte-
ria, and sometimes coenzyme A), thioredoxin and other
enzymes involved in thiol management. The SufL(YajL)
chaperone forms mixed disulphides with target proteins,
preferentially at sulfenylated cysteines (Gautier et al.,
2012). This sulfur oxidized form is also an entry point in
the biosynthesis of a likely ROS buffer, ergothioneine,
made by fungi and bacteria (Song et al., 2015), opening
up a new track for oxidized sulfur metabolism.

When it happens (this appears to be rare in Bacteria,
but ought to be explored in-depth), reduction in the sub-
sequent oxidation stage, sulfinic acid, requires a stronger
reducing power, and enzymes of the energy-dependent
sulfiredoxin family, acting mainly on oxidized peroxire-
doxin (Boileau et al., 2011). Because these modifications
are pervasive they have been recruited by bacteria as
sensors of ROSs and other reactive species beside their
role as buffers against their deleterious action. Finally,
when oxidation results in formation of cysteic acid, the
modification is irreversible and the protein must be
degraded. Subsequently, cysteic acid may be exported
out of the cell or used directly, or indirectly after decar-
boxylation into taurine, as a sulfur source via desulfation
(Cook et al., 1998). Coenzyme A, which is not always
protected by carrying an acyl-group, is also likely to suf-
fer similar accidents. Almost nothing is known about the

fate of such CoA derivatives except for it common disul-
fide form (Wallace et al., 2012) and the highest oxidation
state, which is likely to be catabolized as other sul-
fonates (Nakamura and Tamura, 1972).

Beside oxidation of cysteine, two diastereoisomers of methionine sulfoxide, (S)- and (R)- may be formed in oxidized proteins, depending on the radical and on the environment (Davies, 2005). Remarkably, to cope with this widespread situation, cells evolved two geometrically symmetrical catalytic sites, stemming from completely different protein descents but displaying mirror symmetry of their catalytic site, methionine sulfoxide reductase A and B (Achilli et al., 2015). There is even another class of these enzymes, belonging to molybdoenzyme families, MsrPQ (see also below), for repairing proteins found in the periplasm of diderm bacteria (Gennaris et al., 2015). Yet, despite the existence of such ROS reversible ‘sponge’-like behaviour of methionine, it is expected that some oxidized proteins will be degraded, releasing free methionine sulfoxides in the cell. Furthermore, free methionine can be oxidized directly as well. Now, methionine sulfoxide is an analogue of glutamine (Chin and Means, 1996) that could enter translation, and this needs to be dealt with. In fact, the structural similarity between oxidized methionines and glutamine may explain why glutamine tRNA synthetase did not spread to all bacterial clades. The existence of a likely ancient homeotopic modification (Danchin, 1989) of glutamate by an amidotransferase on glu~tRNAgln in many bacterial clades (e.g. Firmicutes (Strauch et al., 1988)) may be accounted for by avoidance of misacylation of tRNAgln by oxidized methionine. More generally, cells cope with this oxidation using reductases that act on free methionine sulfoxides. In Enterobacteria, for example, an enzyme that repairs oxidized biotin but not protein oxidized methionine, biotin sulfoxide reductase BisC (Denkel et al., 2013), is sufficient for growth on methionine-(S)-sulfoxide (Ezraty et al., 2005). Further, another member of the Msr family (fRMsr), MsrC(YebR) (Etienne et al., 2003; Lin et al., 2007), is present in *Escherichia coli* extracts and reduces the other diastereoisomer, free methionine-(R)-sulfoxide, to methionine. Subsequently, oxidation could go further to reach methionine sulfone (Fig. 1), a compound that is fairly stable and cannot be used as a methionine substitute in animals or a great many bacteria (Or-Rashid et al., 2003). Specific peptidases allow release of free methionine sulfone that is then excreted in the environment (Payne and Tuffnell, 1981). The very fact that the environment is not flooded by this compound demonstrates that, at least for some microorganisms, methionine sulfone could be degraded and presumably used as a sulfur source. We have to go back in time to find work that identified some of the corresponding catabolism (Tonzetich, 1976).
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Progressive oxidation of cysteine and methionine residues.</figcaption></figure>

A similar situation is met with oxidized derivatives of tryptophan or histidine, for example, but for the sake of space, I do not further explore here the way cells cope with these metabolic alterations (see (Lin et al., 2016) for cell management of 2-oxohistidine peptides), noting, however, that this is of considerable interest for microbiome studies (Kennedy et al., 2016).

## Nucleic acids: the example of purines

By contrast, to show how our present knowledge could be used in similar contexts, I dwell some more on derivatives of nucleotides, as they play a central role in energy management, mutagenesis and regulation, with many genes involved in the proofreading and maintenance processes. The list of the corresponding metabolites is far too long to be described here more than in a sketchy way, be it only because RNA nucleotides can have a considerable number of modifications (Machnicka et al., 2015). Remarkably, almost nothing is known about the fate of these nucleotides when they are released in the cell after RNA (or, to a smaller extent, DNA) degradation, which happens often concomitantly with inactivation of the corresponding molecules by reactive chemical species, radicals in particular. Now, again, the environment would be replete with these molecules if there did not exist fairly efficient degradation pathways: we certainly can expect that, among the many genes of unknown function that keep being discovered as more genome and metagenome sequences are deciphered, we will find enzymes involved in the corresponding processes.

As a case in point, here is the fate of some purine derivatives. Guanine derivatives (as related purine

derivatives) are very sensitive to ROS, in reactions similar to: GTP + hydroxyl radical → 8-oxo-GTP + H⁺ (this also happens within nucleic acids in particular when G is unpaired). When input in DNA, 8-oxoG (a tautomer of 8-hydroxy-guanine, 6,8-dihydroxy-2-amino-purine) will be mutagenic, producing transversions, and when in messenger RNA, producing missense proteins. Beside some examples of reactions involving 8-oxoG, the MetaCyc database lists a small fraction of the consequences of OH⁻ reactivity (alteration of nucleic acids via production of 8-oxo-guanine, 2-hydroxy-adenine (isoguanine) or 5-hydroxy-cytosine, but also 5-hydroxyuridine). All produce metabolites that must be exported or degraded by the cell. The pathways of nucleotide degradation should provide a rich set of novel activities. The cell protects itself first by preventing 8-oxoG assimilation into nucleic acids, via degradation of 8-oxo-(d)GTP into 8-oxo-(d)GMP [using the MutT enzyme, a member of the Nudix hydrolase family (Hamm et al., 2016)]. Enzymes of this family will cope with situations (in particular alkylation of nucleotides) when a nucleotide analogue must be prevented to cross the genome sequence quality control, by hydrolysing pyrophosphate from the triphosphate group of the nucleoside triphosphate analogue. Interestingly, quite a few proteins of this family do not yet have a firmly established function (see, e.g. in *E. coli* for modified pyrimidines, NudG(YnjG) and NudI(YfaO)).

At this point, it becomes important to depart from the current practice, that assumes that functional conservation implies structural similarity, to adopt a functional view (function first, structure second) in order to try and identify enzymes with the expected functions. We know that Firmicutes and Proteobacteria, for example, recruited quite different structures to cope with a variety
of functions associated with replication (Engelen et al., 2012). As a consequence, it should not be surprising that degradation of nucleotides differs between these clades. For example, there is no structural counterpart of *E. coli* guanine deaminase GuaD in *Bacillus subtilis*. Yet, an enzyme with guanine deaminase activity, GuaN, is present in the latter. However, this time, it belongs to the cytidine (hence a pyrimidine) deaminase superfamily (Liaw et al., 2004). Interestingly, using our Phyloprofile algorithm (Engelen et al., 2012), we can see that GuaN is co-evolving with xanthine dehydrogenase and a variety of functions related to nitrogen metabolism. Furthermore, it co-evolves with methionine sulfoxide reductases and oxygenases, substantiating indeed that it is involved in catabolism of oxidized purine derivatives (Table S1). This observation opens up another whole domain of possible structures that can cope with modified nucleotides, nucleosides and nucleic acid bases. Furthermore, of course, this does not preclude the possible existence of other families of structures with similar functions. Remarkably, GuaN co-evolves with a considerable number of proteins with no ascribed function ('Y' proteins, Table S1), suggesting that exploration of the *Bacillus* clade would surely be rewarding in terms of identification of novel catabolic activities.

In parallel, in *E. coli* YrfG preferentially hydrolyses phosphate in purine nucleotides. It could also be used to prevent other modified purines to get into nucleic acids. To be sure, a variety of modified nucleotides could be degraded by proteins of the same family, with YjjG preferring pyrimidines, whereas YieH hydrolysed both purines and pyrimidines as secondary substrates (Kuznetsova et al., 2006). The corresponding family of haloacid dehalogenase-like enzymes should be studied in priority to decipher the large number of expected activities that will be essential to bacterial communities. Indeed, the frequent production of halogenated metabolites in sea environments has been proposed to be related to metabolic accidents (Manley, 2002).

Subsequently, as an expected straightforward functional pathway, a nucleosidase (phosphorylases might also play a similar but energy saving role, as in the methionine salvage pathway (Sekowska et al., 2004)) would liberate 8-oxoguanine. In *E. coli*, deaminase GuaD may, beside guanine, accept 8-oxoguanine as substrate. However, this implies some promiscuity, which may lead to unwanted reactions so that the existence of a specific deaminase seems likely. Interestingly, the deamination of 8-oxoG would directly form uric acid, an entry point into fairly involved but ubiquitous degradation pathways (described in some details below, Fig. 2). It was therefore predicted that this specific deaminase should exist, at least in some organisms. Finding the enzyme proved difficult. It was eventually identified in

*Pseudomonas aeruginosa* (Hall et al., 2010b). Exploration of metagenome data derived from environmental DNA sequences isolated from bacteria in the Sargasso Sea under the Global Ocean Sampling Project identified further amidohydrolases of the 8-oxoguanine deaminase family, catalysing the deamination of isoxanthopterin and pterin-6-carboxylate (Hall et al., 2010a). Guanine and pterins are both metabolically and structurally related, and this is witnessed by promiscuous reactions involving both types of metabolites (Kim et al., 2009). This makes obvious that there is a wide class of enzymes with related functions that may act on pterin or guanine oxidized or nitrosated derivatives. Unknown function protein SsnA in *E. coli*, for example, belongs to the COG0402 family and has the conserved signature PxxVxTHxHxxQ present in amidohydrolases and its gene belongs to an island that may be involved in modified nucleic acid bases break down.

Following guanine deamination, xanthine oxidase (XO) oxidizes hypoxanthine. It also acts on other purines and on pterins and a variety of aldehydes. Under some conditions, it may produce superoxide rather than peroxide. The *E. coli* enzyme XdhABC is an FAD/molybdopterin-dependent dehydrogenase that uses NAD⁺ as a cofactor (EC 1.17.1.4, xanthine dehydrogenase). In diverse organisms, the enzyme can be converted to an O₂-dependent XO form (EC 1.17.3.2). Other XO-type enzymes such as *E. coli* aldehyde dehydrogenase PaoABC harbour an additional [4Fe-4S] cluster. They have a scavenging function acting on aromatic aldehydes, not purines, and act as periplasmic detoxifying enzymes involved in oxidation of aldehydes, possibly on a variety of plant-derived aromatics, but also possibly on some of the modified nucleic acid bases (Correia et al., 2016). Catabolism of aromatics is particularly significant in Pseudomonadales that, beside carbon cycles, use a variety of enzymes to degrade heterocycles. In general, their authentic substrates have not yet been deciphered (see (Belda et al., 2016), and the corresponding regularly updated knowledge base at https://www.genoscope.cns.fr/agc/microscope/image/viewer.php?O_id=2834). Future metabolomic studies should help identify some of the pathways that evolved to cope with unavoidable accumulation of modified nucleic acid degradation products.

The subsequent fate of urate is remarkably diverse. Some of it is illustrated in Fig. 2. Remembering that annotations in databases are often unreliable (Percudani et al., 2013), it is essential to try and follow the chemical constraints in the postulated pathways and reappraise current annotations. In aerobic cells, urate may be oxidized to 5-hydroxyisourate by two distinct enzymes: either a coenzyme-independent urate oxidase (uricase, EC 1.7.3.3) found in eukaryotes and bacteria like
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Degradation of guanine and 8-oxoguanine. Spontaneous reactions are as indicated, showing how specific enzymes were sometimes recruited to prevent them to occur.</figcaption></figure>

B. subtilis (PucL) or a prokaryotic flavoprotein urate
hydroxylase (HpxO) originally found in some Klebsiella
species (Michiel et al., 2013). The B. subtilis enzyme
is fused at its N-terminal end to a domain similar
to Mycobacterium tuberculosis AhpG, an alkyl

hydroxyperoxide reductase that cleaves reactive oxygen
intermediates. This activity would remove hydrogen per-
oxide formed in the uricase-catalysed oxidation of uric
acid to allantoin. Next, PucL catalyses oxidation of uric
acid to 5-hydroxyisourate (HIU) to S-allantoin and
possibly also the stereoselective decarboxylation of 2-oxo-4-hydroxy-4-carboxy-5-ureidoimidazoline (OHCU), created by hydrolase PucK (Schultz et al., 2001). Still another related spontaneous reaction has been suggested: uric acid acting as an antioxidant in the periplasm may produce 5-hydroxyisourate, so that periplasmic HiuH(YedX), that hydrolyses 5-hydroxyisourate, would belong to the arsenal protecting the cell from free radicals. In E. coli and other diderms, HiuH thus produces OHCU which spontaneously decays into a racemic mixture of R- and S-allantoin. Furthermore, nonenzymatic decomposition of both HIU and OHCU does occur and does so non-stereospecifically. On top of this, another spontaneous reaction, non-enzymatic racemization of allantoin has been observed. Genomes will tend to recruit enzymes that could cope with the outcome of such spontaneous reactions or channel them along specific pathways. To be sure, two types of allantoinase have been identified, a metal-dependent form and a metal-independent form. While these two enzymes may differ in mechanism, they both preferentially bind the (S)-enantiomer of allantoin. This suggests that HiuH may also have OHCU decarboxylase activity, although this is not straightforward. Alternatively, or in complement, there may exist an allantoin racemase, as indeed discovered in Klebsiella pneumoniae (French et al., 2011).

Subsequently, allantoin can be degraded via a variety of pathways, also with intermediates that are unstable and interconvert spontaneously, yielding potentially compounds that will be hazardous to the cell metabolism (Fig. 2). As a consequence, a number of enzyme activities were recruited in the different environments of living organisms to cope with this possibly harmful situation. Oftentimes, degradation proceeds via ureidoglycine and/or ureidoglycolate and ends up with glyoxylate, urea and/or ammonium, allowing cells to recover carbon and nitrogen. Another intermediate, oxamate, was identified as one of the end-products of allantoin degradation in Streptococcus allantoicus, and further studies demonstrated that oxamate formation is catalysed by an oxamate transcarbamylase possibly via formation of oxalurate as in B. subtilis (via ureidoglycine aminotransferase) or E. coli. To our knowledge, a direct pathway from ureidoglycine to oxamate is not known. In this context, protein YgeW has a 3D structure related to arginine, ornithine and putrescine carbamoyltransferases, but belonging to a distinct clade, which is indicative of a novel activity. It has been suggested to be an oxamate carbamoyltransferases (transcarbamylase) involved in purine catabolism. However, oxamate, allantoin and other substrates were tested but purified YgeW had no activity on any of them. Yet it might generate oxamate in a S-ureidoglycine degradation pathway: ureidoglycine could not be tested in the assays if the protein precluding identification of

this activity if it was indeed present (Li et al., 2011). Subsequently, oxamate would release oxalate and ammonium, as observed for oxamate amidohydrolase HpxW in Klebsiella michiganensis (Hicks and Ealick, 2016). To substantiate this view, I note that, using the PHYLOPROFILE software in B. subtilis shows that protein YwrD is very similar to the Klebsiella HpxW counterpart and that it evolves with most of the enzymes involved in the xanthine degradation pathway (Table S2).

Remarkably, periplasmic degradation of purines is functionally associated in E. coli with MsrPQ, another set of molybdoenzymes – consistent with their coevolution (Medigue et al., 2005), that I have previously described. The hiuHmsrPQ operon is under the control of the two-component system MsrVW(YedVW) which appears to monitor ROSs (Gennaris et al., 2015) and chlorine (Melnyk et al., 2015) species, in addition to excess copper (Urano et al., 2015). Curiously, while this protection process is acknowledged as an important feature, its explicit consequences are not displayed in-depth in most studies. As can be seen from what is described here, a large variety of metabolic pathways could cope with purine degradation. When investigating the metabolism of organisms that are far from model organisms, it is therefore important to avoid directly transposing previous knowledge, but, rather, to proceed via functional analysis, looking for necessary functions first. In general, assuming catabolic pathways as already known may be very misleading, in particular in systems biology studies (Chavarria et al., 2013).

## Other reactive intermediates

In a study quite remarkable for its implications and that should belong to the library of all those interested in metabolism, the laboratory of Andrew Hanson recently summarized the role of spontaneous chemistry in the cell's context. There are so many possibilities, however, that the authors had to make a choice among their preferred molecules. They established a list of 30 top compounds subject to accidents, prominent in metabolomic studies (Lerma-Ortiz et al., 2016). In contrast to the MetaCyc list, which is essentially based on literature, the Chemical-Damage-MINE database repository of reactive compounds http://minedatabase.mcs.anl.gov/#/top30/S1 is presented in a chemically consistent manner. As a consequence of their potential to disrupt the cell's functions, we may expect that these reactive compounds will require agents (essentially enzymes, hence genes) able to remedy the effects of their reactivity, correcting errors in metabolism while redirecting unwanted compounds back to the core pathways. This should be revealed by genome studies. A well-known example is hydrolysis of 6-phosphogluconolactone, long believed to undergo
rapid spontaneous hydrolysis, but which would react with
intracellular nucleophiles, and which is therefore metabo-
lized by a lactonase (Miclet et al., 2001). I only investi-
gated further here two families of reactive molecules
because they are likely to be important during the pro-
cess of cell ageing, a process that has only been sketch-
ily studied in bacteria: those involving RidA and RidA-
related activities and those involving glycation and
related processes.

Enamine/imine deaminases/hydrolases

Enzyme catalysis removes the obstacle of high energy
of activation in specific reactions via a variety of interme-
diary steps, involving highly reactive compounds, often
as the result of the action of a coenzyme. Pyridoxal
phosphate, in particular, will produce molecules with
reactive double bonds, which under normal circum-
stances easily react with a specific substrate, ending up
in the normal product of the reaction. Yet, in particular
during the many transitions (e.g. temperature) faced by
living organisms, the reactive intermediate may diffuse
out of its enzyme, possibly meeting neighbouring sub-
strates. In transamination, racemization or deamination
reactions, reactive amino acid derivatives such as
aminoacrylate/iminopropionate will be produced. The lab-
oratory of Diana Downs has identified the RidA protein
in *Salmonella enterica* (formerly YjgF/YabJ/UK114, often
still annotated as a ribonuclease) as the major enzyme
for scavenging these intermediates (Lambrecht *et al.*,
2013).

The *E. coli* genome contains five counterparts of RidA,
paralogs TdcF, RutC(YcdK), YoaB and YjgH. TdcF is
likely to be a direct counterpart of RidA, preferably acting
on four carbon metabolites (associated with anaerobic
threonine degradation pathway) instead of three (related
to serine via 2-aminoacrylate). It may be important in
processes that prevent the accumulation of L-*allo*-threo-
nine in the cell. RutC is involved in a catabolic pathway
of pyrimidines that produces 3-aminoacrylate (Kim *et al.*,
2010). YoaB and YjgH evolved together, as well as with
RutC (Table S3). Furthermore, the synteny with *ridA*
(yjgF) around *yjgH* is often conserved. This suggests
that the corresponding activities are perhaps overlap-
ping, dealing with deamination of enamine/imines pro-
duced at intermediary steps in a variety of metabolic
pathways. The role of these proteins may be linked to
changes in osmolarity, as co-evolving genes involved in
osmolarity control are co-evolving with them, while a link
with serine deamination is noticeable in the hydroxypyru-
vate co-evolving genes (Table S3 and Fig. 1 in de Lor-
enzo *et al.*, 2015). Hanson *et al.* also noticed that *ridA*
counterparts are often clustered with genes involved in
carbamoyl-phosphate metabolism, coding for an

unknown activity that future work will understand
(Lerma-Ortiz et al., 2016). Finally, a remarkable observa-
tion may account for diversification of the corresponding
protein sequences: RidA appears to metamorphose into
a molecular chaperone when modified by N-chlorination
(Muller et al., 2015). The role of these enzymes is still
far from completely understood. It links metabolism of
branched-chain amino acids, serine and threonine and
pyrimidines in a network that is not yet fully deciphered.
Derepressibility of the *ilvGEDA* operon in *E. coli* is a
common control point, which is still used as an easy
characterization of *relA* mutations (a *relA* mutant does
not grow on a minimal medium supplemented with either
serine, methionine, glycine or uracil, methionine, leucine,
1 mM each, unless supplemented with isoleucine (Uzan
and Danchin, 1978)). Involvement of these enzymes of
core processes of life ought to be explored in priority.

Carbonyl species and related > C=O containing
molecules

Carbon metabolism in cells produces a large number of
reactive aldehydes and ketones. These compounds are
reactive carbonyl species (RCS), critical cell-damaging
and signalling agents. To those, we may associate the
less reactive alpha-keto acids, among which we find the
precursors of amino acids and their degradation prod-
ucts. Briefly, transamination characterizes anabolism
while decarboxylation associated with formation of thioe-
sters marks catabolism. I do not discuss this further,
except to stress again that the reaction of these com-
pounds with ammonia or primary amines forms imine
derivatives, making Schiff bases and RidA-like enzymes
important to cope with their reactivity (Lerma-Ortiz et al.,
2016). In the reactive carbonyl functional group, >C=O,
the carbon atom has two remaining bonds that may be
occupied by hydrogen, alkyl or aryl substituents. In car-
boxylates, this electron distribution is shared between
the carbon and the two oxygen atoms, making the group
considerably less reactive (unless it can be associated
with a conjugated double bond or as bicarbonate in
specific environments). Some consequences of their
reactivity are now explored.

Reducing sugars. Carbohydrates that can convert into
an open chain with a terminal aldehyde group act as
reducing compounds that react with free amines
(Maillard reaction, initially discovered at high
temperature). This happens quite frequently in the cell,
depending on the environment of the amine group,
resulting in glycated proteins. The carbonyl group of the
sugar often reacts with the amino group of the amino
acid, producing N-substituted glycosylamine and water.
The initial stages of the reaction occur more rapidly with
fructose than with glucose, making fructose derivatives particularly frequent but by no means the only glycated compounds that matter. Glycated proteins are often non-functional and prone to degradation or to aggregation. They further evolve into advanced glycation end-products (AGE) proteins after Amadori reactions that impair proteolytic digestibility and alter protein conformation. While this process is ubiquitous (Mironova et al., 2001; Kram and Finkel, 2015), it has mainly been studied in human being, because of its pathogenic consequences (Ajith and Vinodkumar, 2016) and because Amadori compounds are produced during cooking, a process which is associated with a specific human ability, domestication of fire. The most frequent adducts are fructose-lysine as well as similar adducts of the N-terminal amino acid. They may be repaired in

proteins by specific deglycases (Table 1) such as previously discussed SufL(YajL) in *E. coli*, a glyoxalase that has also deglycase activity (Abdallah *et al.*, 2016). SufL converts alpha-oxoaldehydes to carboxylic acids. Its overexpression enhances cellular protection from exogenously added glyoxals and reduces the glyoxal-dependent increase in intracellular AGEs. As noticed here previously, it also acts as a covalent chaperone for sulfenylated thiol proteins involved in Fe-S construction and maintenance.

Because of widespread glycation, glycated-amino acids generated after proteolysis are ubiquitous compounds that can be used by many bacteria as carbon or nitrogen sources. For example, in *B. subtilis*, the operon *frIBOMND* codes for transport and metabolism of fructose-amino acids (Deppe *et al.*, 2011). A counterpart

<table><caption>**Table 1.** Protective metabolism against methylglyoxal and glyoxal</caption><thead><tr><th colspan="3">Escherichia coli</th><th colspan="3">Bacillus subtilis</th></tr><tr><th>Name</th><th>Activity</th><th>PubMed</th><th>Name</th><th>Activity</th><th>PubMed</th></tr></thead><tbody><tr><td>MgsA</td><td>Methylglyoxal synthase</td><td>19458924,<br>21831320</td><td>MgsA</td><td>Methylglyoxal synthase</td><td>21992469,<br>23894131</td></tr><tr><td>GloA</td><td>Glyoxalase I; hemithioacetal-glutathione<br>isomerase, lactoyl-glutathione forming</td><td>25670698</td><td>GlxA<br>YwbC</td><td>Glyoxalase I; hemithioacetal-bacillithiol<br>isomerase, lactoyl-bacillithiol forming</td><td>24330391,<br>25283443</td></tr><tr><td>GloB<br>YafR</td><td>Lactoyl-glutathione hydrolase</td><td>25670698</td><td>GlxB<br>YurT</td><td>Methylglyoxalase; lactoyl-bacillithiol<br>hydrolase</td><td>24330391,<br>25283443</td></tr><tr><td>HchA YedU</td><td>Deglycase; glyoxalase III and molecular<br>chaperone Hsp31</td><td>26774339,<br>26678554</td><td>—</td><td>—</td><td>—</td></tr><tr><td>ElbB</td><td>Oxoaldehyde oxidase; glyoxalase</td><td>9603997,<br>26678554</td><td>—</td><td>—</td><td>—</td></tr><tr><td>GlxC<br>YhbO</td><td>Glyoxalase O; oxoaldehyde oxidase;<br>deglycase</td><td>26678554,<br>26774339</td><td>SufLB YfkM</td><td>General stress protein 18;<br>deglycase</td><td>24330391,<br>25283443</td></tr><tr><td>SufL<br>YajL<br>ThiJ</td><td>Oxoaldehyde oxidase; glyoxalase;<br>chaperone, protecting proteins in<br>response to oxidative stress;<br>deglycase</td><td>26678554,<br>26774339,<br>27644758</td><td>SufL YraA</td><td>Deglycase; general stress protecting<br>enzyme; protects against<br>methylglyoxal toxicity</td><td>24330391,<br>26774339</td></tr><tr><td>DkgA<br>YqhE</td><td>Beta-keto ester reductase;<br>2,5-diketo-D-gluconate reductase;<br>MGO reductase</td><td>20676725,<br>25108218</td><td>YvgN<br>YvsB</td><td>Promiscuous GO/MGO reductase<br>active in vitro on glyceraldehyde<br>and glyceraldehyde-3-phosphate</td><td>19585557</td></tr><tr><td>YqhD</td><td>1,2-propanediol dehydrogenase; alcohol<br>dehydrogenase; aldehyde<br>dehydrogenase (NAD(P)-dependent)</td><td>20543070</td><td>YugJ YugK</td><td>Duplicated butanol dehydrogenase<br>(verified in Clostridium<br>acetobutylicum)</td><td>1809209</td></tr><tr><td>DkgB<br>YafB</td><td>2,5-diketo-D-gluconate reductase;<br>4-nitrobenzaldehyde reductase;<br>MGO reductase</td><td>22328670</td><td>YtbE</td><td>Promiscuous aldo/keto reductase<br>preferential specificity for derivatives<br>of benzaldehyde</td><td>19585557</td></tr><tr><td>YdjG</td><td>Methylglyoxal reductase<br>(NADH-dependent)</td><td>16077126,<br>22074179</td><td>AkrN<br>YhdN</td><td>NADPH-dependent aldo/keto reductase;<br>protects against methylglyoxal</td><td>24330391</td></tr><tr><td>Gpr<br>YghZ</td><td>L-glyceraldehyde 3-phosphate<br>reductase (NADPH-dependent); MGO<br>reductase; GO and glycolaldehyde<br>reductase</td><td>20015532,<br>23990306</td><td>—</td><td>—</td><td>—</td></tr><tr><td>YajO</td><td>2-carboxybenzaldehyde reductase;<br>deletion increases MGO metabolites</td><td>23990306,<br>25326299</td><td>YrpG</td><td>Aldo-keto reductase; deletion does not<br>change MGO sensitivity</td><td>24330391</td></tr><tr><td>Yqil</td><td>Fimbrial protein involved in detoxification<br>of methylglyoxal</td><td>17846588</td><td>—</td><td>—</td><td>—</td></tr><tr><td>KefG KefB</td><td>Potassium proton antiporter; activated<br>by lactoyl-glutathione</td><td>7934942,<br>21143325</td><td>—</td><td>—</td><td>—</td></tr><tr><td>KhtL YabL</td><td>Proton/potassium antiporter; activated<br>by glutathione (by inference)</td><td>7934942,<br>24603618</td><td>KhtU<br>YhaU</td><td>Proton/potassium antiporter;<br>methylglyoxal resistance;<br>activated by bacillithiol</td><td>23894131,<br>24330391</td></tr></tbody></table>
exists in *E. coli*. Fructoselysine in *E. coli* is phosphorylated by a specific kinase (FrI D), followed by the conversion of fructoselysine-6-phosphate into glucose 6-phosphate and lysine by fructoselysine-6-phosphate deglycase (FrI B). There is considerable variation on this theme, depending on the species (Gemayel *et al.*, 2007). Several bacterial operons comprise a homologue of fructoselysine-6-phosphate deglycase, an homologue of the isomerase domain of glucosamine-6-phosphate synthase, and components of a phosphotransferase system (PTS), presumably transporting the fructose-amino acids, but often with no FrI D homologue (Wiame *et al.*, 2005). In *S. enterica*, deglycases GfrE and GfrF are required for growth on glucoselysine and fructoselysine respectively. Experimental observations are consistent with a pathway in which fructoselysine and glucoselysine are phosphorylated at the C-6 position of the sugar by the GfrABCD PTS as they are transported across the membrane (Miller *et al.*, 2015). Fructosamine kinases are active in the intracellular environment, where they produce 3-deoxyglucosone, which promotes AGE formation because it is seldom metabolized before it interacts with a variety of targets. Because these modifications are frequent in cooked food it is interesting to observe that human commensals may use them as major growth supplies (Bui *et al.*, 2015).

Amadoriases, also known as fructosyl amino acid oxidases (FAODs) or fructosyl amine oxidases (FAOXs), are enzymes found in fungi and bacteria. These enzymes catalyse the de-glycosylation of fructosyl amino acid to yield a free amine, glucosone and hydrogen peroxide. They have been categorized into three groups, depending on their substrate specificity. Group I FAOX are active on alpha-fructosyl amino acids (amino acids glycated on backbone amines), and group II FAOX are active on epsilon-fructosyl amino acids (amino acids glycated on side-chain amines), while group III FAOX show similar activity on either alpha- or epsilon-fructosyl amino acids (Rigoldi *et al.*, 2016). In *E. coli*, enzymes such as SolA (acting *in vitro* on sarcosine and other secondary amino acids), which is similar to *Aspergillus fumigatus* amadoriase, may display this activity. Indeed, its co-evolution profile suggests that its function is related to a large number of transport systems and to metabolism of possibly toxic aldehydes, sulfur metabolism and degradation of oxidized molecules (Table S4).

*Methylglyoxal and other aldehydes and ketones.* One of the most prominent characters of dicarbonyls is their high electrophilicity. Among these compounds, methylglyoxal (MGO) is a major endogenous toxic metabolite. It makes covalent bonds with nucleophilic groups (e.g. found in arginine, cysteine and lysine residues in proteins) (Rabbani and Thornalley, 2012). Glyoxal (GO) is a

counterpart that displays similar toxicity but is part of an active process that may limit its production (glyoxylate shunt (Jahan *et al.*, 2016)). The toxic behaviour of MGO is the more remarkable as a great many organisms have specific enzymes (e.g. in *E. coli* MgsA(YccG) methylglyoxal synthase, EC 4.2.3.3, also present in human beings) meant to synthesize it from dihydroxyacetone phosphate (Kosmachevskaya *et al.*, 2015). Conservation of this function implies that some element in the cognate metabolism is so important that it is maintained despite the toxicity of the molecule. A common hypothesis to account for this somewhat puzzling observation is that the pathway leading to MGO is used for ATP-independent glucose catabolism. It could thus be linked to phosphate homoeostasis. In *B. subtilis*, Crh, an homologue of PTS phophocarrier protein HPr, appears to control MgsA synthesis, linking its role to carbon catabolite repression and phosphate transfer (Landmann *et al.*, 2011). A further hypothesis is that while shunting the last stages of glycolysis, the pathway would limit NADH and possibly serine production, which has harmful consequences (de Lorenzo *et al.*, 2015).

MGO and GO reactions with amino acids produce a variety of adducts resulting in advanced glycation products (Kaur *et al.*, 2016). They also react with nucleic acids (Thornalley, 2003). As a consequence, a considerable number of enzymes have evolved activities inactivating these molecules, whether as their primary target, or via moonlighting. In *E. coli*, for example, we find an arsenal of enzymes that act directly or indirectly on GO and MGO, as well as on a variety of aldehydes, also reactive by themselves (Table 1). A major, fairly widespread, detoxification pathway comprises two steps, involving GloA, glyoxalase I (with a $Ni^{2+}$ cofactor) that forms (*R*)-S-lactoylglutathione from MGO and glutathione, and GloB(YafR), glyoxalase II, which is an hydroxyacylglutathione hydrolase producing *D*-lactate and regenerating glutathione. This pathway is very general, and conserved in *B. subtilis* (*GlxA/GlxB*), despite the considerable phylogenetic distance of this species from *E. coli*. However, there it uses bacillithiol as cofactor instead of glutathione. Further, a different set of enzymes does not use reduced sulfur as cofactor: in *E. coli*, HchA(YedU) is a stationary phase glyoxalase III that catalyses a cofactor-independent conversion of MGO to *D*-lactate. It also repairs MGO- and GO-glycated proteins, releasing repaired proteins and lactate or glycolate respectively (Abdallah *et al.*, 2016). Another enzyme, ElbB, possibly involved in isoprene biosynthesis, has also moonlighting activity on glycated proteins. Finally, the conserved oxoaldehyde oxidase SufL (YraA in *B. subtilis*), again, has already been discussed. In *E. coli*, it has a counterpart, glyoxalase O (*GlxO/YhbO*), YfkM in *B. subtilis*, that converts alpha-oxoaldehydes to
carboxylic acids. It seems also to have deglycase prop-
erties. Overexpression of these enzymes enhances cel-
lular protection from exogenously added GOs and
reduces the GO-dependent increase in intracellular
AGEs. Note that the *B. subtilis* enzymes are not authen-
tic orthologs of their *E. coli* counterparts, as they depend
on bacillithiol not glutathione. In the same way, counter-
parts in Mycobacteria are expected to use mycothiol
instead of glutathione. However, all these enzymes
arose from a common descent and a relevant nomencla-
ture has to be designed to take this cofactor variation
into account.

In line with a role of MGO synthesis in the balance
between reduced and oxidized NAD cofactors, DkgA is
a NADPH-dependent beta-keto ester reductase/2,5-
diketo-D-gluconate reductase that also displays GO/MGO
reductase activity. Its paralog DkgB(YafB) is primarily a
2,5-diketo-D-gluconate reductase that reduces 4-nitro-
benzaldehyde *in vitro* and acts as a MGO reductase.
These enzymes appear to be conserved in *E. coli* and
*B. subtilis* (Table 1). YdjG (AkrN(YhdN) in *B. subtilis*)
is a further aldehyde reductase, coded in a carbohydrate
degradation operon of yet unknown specificity, that also
acts as a NADH-dependent MGO reductase. Several
other enzymes acting on aldehydes may also efficiently
detoxify GO, among which Gpr(YghZ) L-glyceraldehyde
3-phosphate reductase/MGO reductase, and YeaE are
major NADPH-dependent players (Lee *et al.*, 2013). Ulti-
mately, general purpose oxidoreductase proteins such
as AldA behave as NADH-dependent aldehyde dehydro-
genases that act on glycolaldehyde and lactaldehyde,
thus minimizing their toxicity. A great many paralogs of
this enzyme are present in most organisms. Finally, in
*E. coli*, Yqil is a fimbrial protein somehow involved in
detoxification of MGO, possibly as a buffering element
(Table 1). Overall cells display a large supply of oxidore-
ductases with considerable buffering role on aldehydes,
and in particular GO and MGO. It is, however, expected
that because these metabolites are pervasive, these
reactive metabolites will, as time elapses, contribute sig-
nificantly to protein ageing. A thorough study of the con-
tribution of these enzymes that are probably redundant
via moonlighting, should be undertaken.

An additional defence mechanism is worth emphasizing. It has been mainly deciphered in *E. coli* and *B. subtilis* (Chandrangsu *et al.*, 2014), but is likely to operate in a great many other organisms as well. MGO toxicity is highly sensitive to pH. A slight acidification of the cell's cytoplasm is enough to reduce the number of MGO targets (Ferguson *et al.*, 1993). As a consequence, cells developed a proton-based protection system where potassium is exchanged for protons during efflux. This is substantiated by the observation that the (*R*)-S-lactoylglutathione intermediate in MGO detoxification is

required for the activation of the KefGB proton/potassium
efflux pump (Ozyamak *et al.*, 2010). Finally, in *B. sub-
tilis*, the chromosome neighbourhood of the *mgsA* gene,
coding for MGO synthase, is revealing. It comprises the
*bshA* and *bshBA* genes, coding for two enzymes
required for bacillithiol synthesis, embedded within a
seven-gene operon additionally including *mgsA*, and the
essential genes *cca* and *birA*, encoding tRNA
nucleotidyltransferase (CCA transferase) and biotin-pro-
tein ligase respectively (Gaballa *et al.*, 2013). To be
sure, mutants of this system are more sensitive to MGO.

*Fumarate.* While authentic dicarbonyls are usually highly reactive and toxic, carboxylates are generally fairly innocuous because of their low spontaneous reactivity. However, the presence of a conjugated carbon double bond such as that in fumaric acid makes this molecule quite reactive. Indeed, the tricarboxylic acid cycle enzyme fumarate hydratase, that inactivates this reaction trend, has been identified as a tumour suppressor in a subset of human renal cell carcinomas. Among other possible targets, lysine in particular, cysteine residues in proteins are especially prone to react with fumarate, forming S-(2-succinyl)cysteine (Alderson *et al.*, 2006). Subsequently, upon protein degradation, this metabolite is likely to accumulate in cells and either be exported or degraded via a still unknown pathway. This reaction that appears to be commonplace in a variety of human diseases (in particular associated with mitochondria, including obesity and diabetes (Thomas *et al.*, 2012)), is likely to be ubiquitous. Remarkably, it is associated, in Parkinson's disease, with modification of a detoxification family of enzymes that I found in the previous paragraph to be involved in protection against protein glycation (Poschmann *et al.*, 2014). Taken together, these observations show that fumarate reactivity should be revisited carefully in bacterial genomes exploration. Furthermore, because protein degradation would generate unusual non-proteinogenic amino acids, that do not accumulate in the environment, unknown catabolic pathways are likely to be uncovered in the future.

*Carboxylation.* In the cell's internal medium, carbon dioxide is generally present in the hydrated form of bicarbonate. Due to the essential role of the carbon dioxide/bicarbonate buffer system in regulation of physiological pH, CO(3)(·-), it may be one of the most important ROS in biological systems, and it may contribute to formation of 8-oxoG, discussed previously (Roginskaya *et al.*, 2015). Furthermore, as shown by the process of light-dependent carbon fixation, as well as carboxylation of pyruvate in the main anaplerotic reaction in carbon metabolism, carbon dioxide is certainly not devoid of reactivity. As a matter of fact, the
addition of a carboxyl group to the epsilon-amino group of the lysine side chain is a frequent but widely overlooked post-translational modification with remarkable consequences. Lysine is positively charged at physiological pH. However, the carboxyl group inverts the electrochemical properties of the lysine side chain by complementing the roles of glutamate and aspartate in the fine tuning of enzyme activity. Providing a negative charge instead of a positive one on a long side group expands the catalytic capabilities of the amino acids in proteins beyond the 20 common amino acids. Carboxylated lysine is the structural counterpart to arginine, except that its side chain is negatively charged at the pH range useful for enzymes. Experiments with urease demonstrated that this modification occurs spontaneously under basic pH conditions involving carbon dioxide (Park and Hausinger, 1995) and without the mediation of any enzyme. Witnessing the importance of this modification, the class D beta-lactamases are characterized by the presence of a carboxylated lysine in the active site that participates in catalysis (Che et al., 2014). Despite its likely significance, the extent of lysine carboxylation in proteins is currently unknown. A lysine carboxylation software has been designed. Its outcome running on a variety of conditions suggests that about 1.3% of large proteins may contain a carboxylated lysine residue (Jimenez-Morales et al., 2014). This amply demonstrates the importance of putting some focus on this type of post-translational modifications.

Conclusion

As new genome sequences keep accumulating, the number of genes with no ascribed functions increases steadily. This is such a large amount that we cannot expect to be able to set up experimental approaches that will allow us to identify their function. Identification of gene function is somewhat similar to criminal profiling (Verde and Nurra, 2010). We have to use hypotheses, data and context to try and propose a candidate. In this quest, we have three ways to go forward: deduction, which uses a preset hypothesis imposed upon the physi- ological, biochemical and genetic observations associ- ated with the corresponding gene and deduce what it must be. Induction, that proceeds the other way around, tries to infer a function by association. Finally, abduction, a shot-in-the-dark approach (you carry a gun, while in a forest, and fire: if something cries, you go there and see what you got...), applies a harsh treatment to the data and looks for the way it behaves. The most efficient way forward is induction. Yet most genomic studies tend to unknowingly use abduction, with a considerable amount of spurious observations that plague the future of predic- tions. The design of knowledge bases constructed on

the concept of neighbourhoods (Nitschke et al., 1998) must be further developed to improve our ability to make valid inferences. I have tried here to create a picture of some relevant neighbourhoods. I hoped to show how this approach might be used in future investigations of ‘omics’ approaches meant to sort out the fate of normal accidents in metabolism.

Acknowledgements

I thank Victor de Lorenzo, Philippe Marlière and Pablo
Ivan Nikel for suggestions and critical reading of a first
version of the manuscript.

Conflict of interest

None declared.

References

Abdallah, J., Mihoub, M., Gautier, V., and Richarme, G. (2016) The DJ-1 superfamily members YhbO and YajL from *Escherichia coli* repair proteins from glycation by methylglyoxal and glyoxal. *Biochem Biophys Res Commun* **470**: 282–286.

Achilli, C., Ciana, A., and Minetti, G. (2015) The discovery of methionine sulfoxide reductase enzymes: an historical account and future perspectives. *BioFactors* **41**: 135–152.

Ajith, T.A., and Vinodkumar, P. (2016) Advanced glycation end products: association with the pathogenesis of diseases and the current therapeutic advances. *Curr Clin Pharmacol* **11**: 118–127.

Alderson, N.L., Wang, Y., Blatnik, M., Frizzell, N., Walla, M.D., Lyons, T.J., et al. (2006) S-(2-Succinyl)cysteine: a novel chemical modification of tissue proteins by a Krebs cycle intermediate. *Arch Biochem Biophys* **450**: 1–8.

Anglada, J.M., Martins-Costa, M., Francisco, J.S., and Ruiz-Lopez, M.F. (2015) Interconnection of reactive oxygen species chemistry across the interfaces of atmospheric, environmental, and biological processes. *Acc Chem Res* **48**: 575–583.

Belda, E., van Heck, R.G., Lopez-Sanchez, M.J., Cruveiller, S., Barbe, V., Fraser, C., et al. (2016) The revisited genome of *Pseudomonas putida* KT2440 enlightens its value as a robust metabolic chassis. *Environ Microbiol* **18**: 3403–3424.

van Bergen, L.A., Roos, G., and De Proft, F. (2014) From thiol to sulfonic acid: modeling the oxidation pathway of protein thiols by hydrogen peroxide. *J Phys Chem A* **118**: 6078–6084.

Boileau, C., Eme, L., Brochier-Armanet, C., Janicki, A., Zhang, C.C., and Latifi, A. (2011) A eukaryotic-like sulfiredoxin involved in oxidative stress responses and in the reduction of the sulfinic form of 2-Cys peroxiredoxin in the cyanobacterium *Anabaena* PCC 7120. *New Phytol* **191**: 1108–1118.

Bui, T.P., Ritari, J., Boeren, S., de Waard, P., Plugge, C.M., and de Vos, W.M. (2015) Production of butyrate from
lysine and the Amadori product fructoselysine by a human gut commensal. *Nat Commun* **6**: 10062.

Caspi, R., Billington, R., Ferrer, L., Foerster, H., Fulcher, C.A., Keseler, I.M., et al. (2016) The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of pathway/genome databases. *Nucleic Acids Res* **44**: D471–D480.

Chan, C.M., Danchin, A., Marliere, P., and Sekowska, A. (2014) Paralogous metabolism: S-alkyl-cysteine degradation in *Bacillus subtilis*. *Environ Microbiol* **16**: 101–117.

Chandrangsu, P., Dusi, R., Hamilton, C.J., and Helmann, J.D. (2014) Methylglyoxal resistance in *Bacillus subtilis*: contributions of bacillithiol-dependent and independent pathways. *Mol Microbiol* **91**: 706–715.

Chavarria, M., Nikel, P.I., Perez-Pantoja, D., and de Lor- enzo, V. (2013) The Entner-Doudoroff pathway empowers *Pseudomonas putida* KT2440 with a high tolerance to oxidative stress. *Environ Microbiol* **15**: 1772–1785.

Che, T., Bethel, C.R., Pusztai-Carey, M., Bonomo, R.A., and Carey, P.R. (2014) The different inhibition mechanisms of OXA-1 and OXA-24 beta-lactamases are determined by the stability of active site carboxylated lysine. *J Biol Chem* **289**: 6152–6164.

Chin, D., and Means, A.R. (1996) Methionine to glutamine substitutions in the C-terminal domain of calmodulin impair the activation of three protein kinases. *J Biol Chem* **271**: 30465–30471.

Cook, A.M., Laue, H., and Junker, F. (1998) Microbial desulfonation. *FEMS Microbiol Rev* **22**: 399–419.

Correia, M.A., Otrelo-Cardoso, A.R., Schwuchow, V., Sigfridsson Clauss, K.G., Haumann, M., Romao, M.J., et al. (2016) The *Escherichia coli* periplasmic aldehyde oxidoreductase is an exceptional member of the xanthine oxidase family of molybdoenzymes. *ACS Chem Biol* **11**: 2923–2935.

Danchin, A. (1989) Homeotopic transformation and the origin of translation. *Prog Biophys Mol Biol* **54**: 81–86.

Danchin, A., and Fang, G. (2016) Unknown unknowns: essential genes in quest for function. *Microb Biotechnol* **9**: 530–540.

Danchin, A., and Sekowska, A. (2014) The logic of metabolism and its fuzzy consequences. *Environ Microbiol* **16**: 19–28.

Danchin, A., Binder, P.M., and Noria, S. (2011) Antifragility and tinkering in biology (and in business) flexibility provides an efficient epigenetic way to manage risk. *Genes (Basel)* **2**: 998–1016.

D'Ari, R., and Casadesus, J. (1998) Underground metabolism. *BioEssays* **20**: 181–186.

Davies, M.J. (2005) The oxidative environment and protein damage. *Biochim Biophys Acta* **1703**: 93–109.

Denkel, L.A., Rhen, M., and Bange, F.C. (2013) Biotin sulfoxide reductase contributes to oxidative stress tolerance and virulence in *Salmonella enterica* serovar Typhimurium. *Microbiology* **159**: 1447–1458.

Deppe, V.M., Klatte, S., Bongaerts, J., Maurer, K.H., O'Connell, T., and Meinhardt, F. (2011) Genetic control of amadori product degradation in *Bacillus subtilis* via regulation of *frlBONMD* expression by *FrIR*. *Appl Environ Microbiol* **77**: 2839–2846.

Eady, R.R., Antonyuk, S.V., and Hasnain, S.S. (2016) Fresh insight to functioning of selected enzymes of the nitrogen cycle. *Curr Opin Chem Biol* **31**: 103–112.

Engelen, S., Vallenet, D., Medigue, C., and Danchin, A. (2012) Distinct co-evolution patterns of genes associated to DNA polymerase III DnaE and PolC. *BMC Genom* **13**: 69.

Etienne, F., Spector, D., Brot, N., and Weissbach, H. (2003) A methionine sulfoxide reductase in *Escherichia coli* that reduces the R enantiomer of methionine sulfoxide. *Biochem Biophys Res Commun* **300**: 378–382.

Ezraty, B., Bos, J., Barras, F., and Aussel, L. (2005) Methionine sulfoxide reduction and assimilation in *Escherichia coli*: new role for the biotin sulfoxide reductase BisC. *J Bacteriol* **187**: 231–237.

Ferguson, G.P., Munro, A.W., Douglas, R.M., McLaggan, D., and Booth, I.R. (1993) Activation of potassium channels during metabolite detoxification in *Escherichia coli*. *Mol Microbiol* **9**: 1297–1303.

Frawley, E.R., and Fang, F.C. (2014) The ins and outs of bacterial iron metabolism. *Mol Microbiol* **93**: 609–616.

French, J.B., Neau, D.B., and Ealick, S.E. (2011) Characterization of the structure and function of *Klebsiella pneumoniae* allantoin racemase. *J Mol Biol* **410**: 447–460.

Gaballa, A., Antelmann, H., Hamilton, C.J., and Helmann, J.D. (2013) Regulation of *Bacillus subtilis* bacillithiol biosynthesis operons by Spx. *Microbiology* **159**: 2025–2035.

Gautier, V., Le, H.T., Malki, A., Messaoudi, N., Caldas, T., Kthiri, F., et al. (2012) *YajL*, the prokaryotic homolog of the Parkinsonism-associated protein DJ-1, protects cells against protein sulfenylation. *J Mol Biol* **421**: 662–670.

Gemayel, R., Fortpied, J., Rzem, R., Vertommen, D., Veigada-Cunha, M., and Van Schaftingen, E. (2007) Many fructosamine 3-kinase homologues in bacteria are ribulosamine/erythrulosamine 3-kinases potentially involved in protein deglycation. *FEBS J* **274**: 4360–4374.

Gennaris, A., Ezraty, B., Henry, C., Agrebi, R., Vergnes, A., Oheix, E., et al. (2015) Repairing oxidized proteins in the bacterial envelope using respiratory chain electrons. *Nature* **528**: 409–412.

Giles, G.I., and Jacob, C. (2002) Reactive sulfur species: an emerging concept in oxidative stress. *Biol Chem* **383**: 375–388.

Giuffre, A., Borisov, V.B., Arese, M., Sarti, P., and Forte, E. (2014) Cytochrome bd oxidase and bacterial tolerance to oxidative and nitrosative stress. *Biochim Biophys Acta* **1837**: 1178–1187.

Gupta, V., and Carroll, K.S. (2014) Sulfenic acid chemistry, detection and cellular lifetime. *Biochim Biophys Acta* **1840**: 847–875.

Hall, R.S., Agarwal, R., Hitchcock, D., Sauder, J.M., Burley, S.K., Swaminathan, S., and Raushel, F.M. (2010a) Discovery and structure determination of the orphan enzyme isoxanthopterin deaminase. *Biochemistry* **49**: 4374–4382.

Hall, R.S., Fedorov, A.A., Marti-Arbona, R., Fedorov, E.V., Kolb, P., Sauder, J.M., et al. (2010b) The hunt for 8-oxoguanine deaminase. *J Am Chem Soc* **132**: 1762–1763.

Hamm, M.L., McFadden, E.J., Ghio, M., Lindell, M.A., Gérien, K.S., and O'Handley, S.F. (2016) Insights into the substrate specificity of the MutT pyrophosphohydrolase
using structural analogues of 8-oxo-2'-deoxyguanosine
nucleotide. *Bioorg Med Chem Lett* **26**: 2014–2017.

Hicks, K.A., and Ealick, S.E. (2016) Biochemical and struc-
tural characterization of *Klebsiella pneumoniae* oxamate
amidohydrolase in the uric acid degradation pathway.
*Acta Crystallogr D Struct Biol* **72**: 808–816.

Jahan, N., Maeda, K., Matsuoka, Y., Sugimoto, Y., and Kurata, H. (2016) Development of an accurate kinetic model for the central carbon metabolism of *Escherichia coli*. *Microb Cell Fact* **15**: 112.

Jimenez-Morales, D., Adamian, L., Shi, D., and Liang, J. (2014) Lysine carboxylation: unveiling a spontaneous post-translational modification. *Acta Crystallogr D Biol Crystallogr* **70**: 48–57.

Kaur, H., Kamalov, M. and Brimble, M.A. (2016) Chemical synthesis of peptides containing site-specific advanced glycation endproducts. *Acc Chem Res* **49**: 2199–2208.

Kennedy, P.J., Cryan, J.F., Dinan, T.G. and Clarke, G. (2016) Kynurenine pathway metabolism and the microbiota-gut-brain axis. *Neuropharmacology* (in press).

Kim, J., Park, S.I., Ahn, C., Kim, H., and Yim, J. (2009) Guanine deaminase functions as dihydropterin deaminase in the biosynthesis of aurod rosopterin, a minor red eye pigment of Drosophila. *J Biol Chem* **284**: 23426–23435.

Kim, K.S., Pelton, J.G., Inwood, W.B., Andersen, U., Kustu, S., and Wemmer, D.E. (2010) The Rut pathway for pyrimidine degradation: novel chemistry and toxicity problems. *J Bacteriol* **192**: 4089–4102.

Kosmachevskaya, O.V., Shumaev, K.B., and Topunov, A.F. (2015) Carbonyl stress in bacteria: causes and consequences. *Biochemistry (Mosc)* **80**: 1655–1671.

Kram, K.E., and Finkel, S.E. (2015) Rich medium composition affects *Escherichia coli* survival, glycation, and mutation frequency during long-term batch culture. *Appl Environ Microbiol* **81**: 4442–4450.

Kuznetsova, E., Proudfoot, M., Gonzalez, C.F., Brown, G., Omelchenko, M.V., Borozan, I., et al. (2006) Genome-wide analysis of substrate specificities of the *Escherichia coli* haloacid dehalogenase-like phosphatase family. *J Biol Chem* **281**: 36149–36161.

Lambrecht, J.A., Schmitz, G.E., and Downs, D.M. (2013) RidA proteins prevent metabolic damage inflicted by PLP-dependent dehydratases in all domains of life. *MBio* **4**: e00033–13.

Landmann, J.J., Busse, R.A., Latz, J.H., Singh, K.D., Stulke, J., and Gorke, B. (2011) Crh, the paralogue of the phosphocarrier protein HPr, controls the methylglyoxal bypass of glycolysis in *Bacillus subtilis*. *Mol Microbiol* **82**: 770–787.

Lee, C., Kim, I., and Park, C. (2013) Glyoxal detoxification in *Escherichia coli* K-12 by NADPH dependent aldo-keto reductases. *J Microbiol* **51**: 527–530.

Lerma-Ortiz, C., Jeffryes, J.G., Cooper, A.J., Niehaus, T.D., Thamm, A.M., Frelin, O., et al. (2016) ‘Nothing of chemistry disappears in biology’: the Top 30 damage-prone endogenous metabolites. *Biochem Soc Trans* **44**: 961–971.

Levitan, I., Volkov, S., and Subbaiah, P.V. (2010) Oxidized LDL: diversity, patterns of recognition, and pathophysiology. *Antioxid Redox Signal* **13**: 39–75.

Li, Y., Jin, Z., Yu, X., Allewell, N.M., Tuchman, M., and Shi, D. (2011) The *ygeW* encoded protein from *Escherichia coli* is a knotted ancestral catabolic transcarbamylase. *Proteins* **79**: 2327–2334.

Liaw, S.H., Chang, Y.J., Lai, C.T., Chang, H.C., and Chang, G.G. (2004) Crystal structure of *Bacillus subtilis* guanine deaminase: the first domain-swapped structure in the cytidine deaminase superfamily. *J Biol Chem* **279**: 35479–35485.

Lin, Z., Johnson, L.C., Weissbach, H., Brot, N., Lively, M.O., and Lowther, W.T. (2007) Free methionine-(R)-sulfoxide reductase from *Escherichia coli* reveals a new GAF domain function. *Proc Natl Acad Sci U S A* **104**: 9597–9602.

Lin, J.M., Tsai, Y.T., Liu, Y.H., Lin, Y., Tai, H.C. and Chen, C.S. (2016) Identification of 2-oxohistidine interacting proteins using *E. coli* proteome chips. *Mol Cell Proteomics* (in press).

Lisher, J.P., and Giedroc, D.P. (2013) Manganese acquisition and homeostasis at the host-pathogen interface. *Front Cell Infect Microbiol* **3**: 91.

de Lorenzo, V., Sekowska, A., and Danchin, A. (2015) Chemical reactivity drives spatiotemporal organisation of bacterial metabolism. *FEMS Microbiol Rev* **39**: 96–119.

Lu, Z., Bourdi, M., Li, J.H., Aponte, A.M., Chen, Y., Lom bard, D.B., et al. (2011) SIRT3-dependent deacetylation exacerbates acetaminophen hepatotoxicity. *EMBO Rep* **12**: 840–846.

Machnicka, M.A., Olchowik, A., Grosjean, H., and Bujnicki, J.M. (2015) Distribution and frequencies of post-transcriptional modifications in tRNAs. *RNA Biol* **11**: 1619–1629.

Manley, S. (2002) Phytogenesis of halomethanes: a product of selection or a metabolic accident? *Biogeochemistry* **60**: 163–180.

Medigue, C., Krin, E., Pascal, G., Barbe, V., Bernsel, A., Ber-
tin, P.N., et al. (2005) Coping with cold: the genome of the
versatile marine Antarctica bacterium *Pseudoalteromonas
haloplanktis* TAC125. *Genome Res* **15**: 1325–1335.

Melnyk, R.A., Youngblut, M.D., Clark, I.C., Carlson, H.K., Wetmore, K.M., Price, M.N., et al. (2015) Novel mechanism for scavenging of hypochlorite involving a periplasmic methionine-rich peptide and methionine sulfoxide reductase. *MBio* **6**: e00233–15.

Michiel, M., Perchat, N., Perret, A., Tricot, S., Papeil, A., Besnard, M., *et al.* (2013) Microbial urate catabolism: characterization of HpyO, a non-homologous isofunctional isoform of the flavoprotein urate hydroxylase HpxO. *Environ Microbiol Rep* **4**: 642–647.

Miclet, E., Stoven, V., Michels, P.A., Opperdoes, F.R., Lallemand, J.Y., and Duffieux, F. (2001) NMR spectroscopic analysis of the first two steps of the pentose-phosphate pathway elucidates the role of 6-phosphogluconolactonase. *J Biol Chem* **276**: 34840–34846.

Miller, K.A., Phillips, R.S., Kilgore, P.B., Smith, G.L., and Hoover, T.R. (2015) A mannose family phosphotransferase system permease and associated enzymes are required for utilization of fructoselysine and glucoselysine in *Salmonella enterica* serovar Typhimurium. *J Bacteriol* **197**: 2831–2839.

Mironova, R., Niwa, T., Hayashi, H., Dimitrova, R., and Iva-
nov, I. (2001) Evidence for non-enzymatic glycosylation in
*Escherichia coli*. *Mol Microbiol* **39**: 1061–1068.
Muller, A., Langklotz, S., Lupilova, N., Kuhlmann, K., Bando, J.E., and Leichert, L.I. (2015) Activation of RidA chaperone function by N-chlorination. *Nat Commun* **5**: 5804.

Nakamura, H., and Tamura, Z. (1972) Growth responses of *Bifidobacterium bifidum* to S-sulfonic acid-type pantetheine related compounds. *Jpn J Microbiol* **16**: 239-242.

Nitschke, P., Guerdoux-Jamet, P., Chiapello, H., Faroux, G., Henaut, C., Henaut, A., and Danchin, A. (1998) Indigo: a World-Wide-Web review of genomes and gene functions. *FEMS Microbiol Rev* **22**: 207-227.

Or-Rashid, M.M., Onodera, R., and Wadud, S. (2003) Studies on the utilization of methionine sulfoxide and methionine sulfone by rumen microorganisms *in vitro*. *Amino Acids* **24**: 135-139.

Ozyamak, E., Black, S.S., Walker, C.A., Maclean, M.J., Bartlett, W., Miller, S., and Booth, I.R. (2010) The critical role of S-lactoylglutathione formation during methylglyoxal detoxification in *Escherichia coli*. *Mol Microbiol* **78**: 1577-1590.

Park, I.S., and Hausinger, R.P. (1995) Requirement of carbon dioxide for *in vitro* assembly of the urease nickel metallocenter. *Science* **267**: 1156-1158.

Payne, J.W., and Tuffnell, J.M. (1981) Cleavage of peptide-bound methionine sulphoxide and methionine sulphone: possible role in enhanced proteolysis in *Escherichia coli*. *FEMS Microbiol Lett* **12**: 279-282.

Percudani, R., Carnevali, D. and Puggioni, V. (2013) Ureido-glycolate hydrolase, amidohydrolase, lyase: how errors in biological databases are incorporated in scientific papers and vice versa. *Database* (Oxford), 2013: bat071.

Perrow, C. (1984) *Normal Accidents: Living with High-Risk Technologies*. New York, NY: Basic Books.

Poschmann, G., Seyfarth, K., Besong Agbo, D., Klafki, H.W., Rozman, J., Wurst, W., et al. (2014) High-fat diet induced isoform changes of the Parkinson's disease protein DJ-1. *J Proteome Res* **13**: 2339-2351.

Pullan, S.T., Monk, C.E., Lee, L., and Poole, R.K. (2008) Microbial responses to nitric oxide and nitrosative stress: growth, "omic", and physiological methods. *Methods Enzymol* **437**: 499-519.

Rabbani, N., and Thornalley, P.J. (2012) Methylglyoxal, glyoxalase 1 and the dicarbonyl proteome. *Amino Acids* **42**: 1133-1142.

Rigoldi, F., Gautieri, A., Dalle Vedove, A., Lucarelli, A.P., Vesentini, S., and Parisini, E. (2016) Crystal structure of the deglycating enzyme Amadoriase I in its free form and substrate-bound complex. *Proteins* **84**: 744-758.

Roginskaya, M., Moore, T.J., Ampadu-Boateng, D., and Razskazovskiy, Y. (2015) Efficacy and site specificity of hydrogen abstraction from DNA 2-deoxyribose by carbonate radicals. *Free Radic Res* **49**: 1431-1437.

Schultz, A.C., Nygaard, P., and Saxild, H.H. (2001) Functional analysis of 14 genes that constitute the purine catabolic pathway in *Bacillus subtilis* and evidence for a novel regulon controlled by the PucR transcription activator. *J Bacteriol* **183**: 3293-3302.

Sekowska, A., Denervaud, V., Ashida, H., Michoud, K., Haas, D., Yokota, A., and Danchin, A. (2004) Bacterial variations on the methionine salvage pathway. *BMC Microbiol* **4**: 9.

Song, H., Hu, W., Naowarojna, N., Her, A.S., Wang, S., Desai, R., et al. (2015) Mechanistic studies of a novel C-S lyase in ergothioneine biosynthesis: the involvement of a sulfenic acid intermediate. *Sci Rep* **5**: 11870.

Strauch, M.A., Zalkin, H., and Aronson, A.I. (1988) Characterization of the glutamyl-tRNA(Gln)-to-glutaminyl-tRNA (Gln) amidotransferase reaction of *Bacillus subtilis*. *J Bacteriol* **170**: 916-920.

Thomas, S.A., Storey, K.B., Baynes, J.W., and Frizzell, N. (2012) Tissue distribution of S-(2-succino)cysteine (2SC), a biomarker of mitochondrial stress in obesity and diabetes. *Obesity* (Silver Spring) **20**: 263-269.

Thornalley, P.J. (2003) Protecting the genome: defence against nucleotide glycation and emerging role of glyoxalase I overexpression in multidrug resistance in cancer chemotherapy. *Biochem Soc Trans* **31**: 1372-1377.

Tonzetich, J. (1976) Chromatographic separation of methionine, methionine sulphoxide, methionine sulphone, and their products of oral microbial metabolism. *Anal Biochem* **73**: 290-300.

Urano, H., Umezawa, Y., Yamamoto, K., Ishihama, A., and Ogasawara, H. (2015) Cooperative regulation of the common target genes between H(2)O(2)-sensing YedVW and Cu(2)(+)-sensing CusSR in *Escherichia coli*. *Microbiology* **161**: 729-738.

Uzan, M., and Danchin, A. (1978) Correlation between the serine sensitivity and the derepressibility of the *ilv* genes in *Escherichia coli* relA- mutants. *Mol Gen Genet* **165**: 21-30.

Verde, A., and Nurra, A. (2010) Criminal profiling as a plotting activity based on abductive processes. *Int J Offender Ther Comp Criminol* **54**: 829-849.

Wallace, B.D., Edwards, J.S., Wallen, J.R., Moolman, W.J., van der Westhuyzen, R., Strauss, E., et al. (2012) Turn-over-dependent covalent inactivation of *Staphylococcus aureus* coenzyme A-disulfide reductase by coenzyme A-mimetics: mechanistic and structural insights. *Biochemistry* **51**: 7699-7711.

Wiame, E., Lamosa, P., Santos, H., and Van Schaftingen, E. (2005) Identification of glucoselysine-6-phosphate deglycase, an enzyme involved in the metabolism of the fruc-tation product glucoselysine. *Biochem J* **392**: 263-269.

## Supporting information

Additional Supporting Information may be found online in the supporting information tab for this article:

**Table S1.** Co-evolution pattern of *B. subtilis* guanine deaminase GuaN. The protein co-evolves with enzymes involved in purine degradation. Note the very large number of proteins of unknown function (213/500 "Y" proteins) suggesting that exploration of cognate pathways in Firmicutes will be rewarding.

**Table S2.** Co-evolution pattern of *B. subtilis* amidohydrolase YwrD. Note that the protein evolves with most of the enzymes involved in purine catabolism. Many of the 141 remaining "Y" proteins have a related activity that can be identified. Nomenclature should be modified accordingly.
**Table S3.** Functional co-evolution of YjgH in *E. coli*. The protein co-evolves with RutC and YoaB. Many proteins of the 136 remaining “Y” proteins have an identified activity, showing that the nomenclature needs to be modified urgently. Furthermore those associated to oxygenase activities should perhaps be studied in priority.

**Table S4.** Co-evolution pattern of *E. coli* SolA protein. SolA may be an amadoriase. It is associated to co-evolution of 163 “Y” proteins, many of them putative transporters. In addition several are putative phosphatases, that may correspond to phosphorylated carbohydrate derivatives.