Underground isoleucine biosynthesis pathways in E. coli
Abstract
The promiscuous activities of enzymes provide fertile ground for the evolution of new metabolic pathways. Here, we systematically explore the ability of E. coli to harness underground metabolism to compensate for the deletion of an essential biosynthetic pathway. By deleting all threonine deaminases, we generated a strain in which isoleucine biosynthesis was interrupted at the level of 2-ketobutyrate. Incubation of this strain under aerobic conditions resulted in the emergence of a novel 2-ketobutyrate biosynthesis pathway based upon the promiscuous cleavage of O-succinyl-L-homoserine by cystathionine γ-synthase (MetB). Under anaerobic conditions, pyruvate formate-lyase enabled 2-ketobutyrate biosynthesis from propionyl-CoA and formate. Surprisingly, we found this anaerobic route to provide a substantial fraction of isoleucine in a wild-type strain when propionate is available in the medium. This study demonstrates the selective advantage underground metabolism offers, providing metabolic redundancy and flexibility which allow for the best use of environmental carbon sources.
Introduction
The patchwork model suggests that novel metabolic pathways emerge from the promiscuous activities of enzymes participating in diverse metabolic processes (Lazcano and Miller, 1999; Jensen, 1976; Khersonsky and Tawfik, 2010; Noda-Garcia et al., 2018). According to this, underground metabolism – the network of metabolic conversions which are catalyzed as side reactions of enzymes that have evolved to support other activities (D'Ari and Casadesús, 1998) – provides fertile ground for the evolution of new pathways. A computational analysis suggested that about half of all underground reactions generate metabolites that already exist in the endogenous metabolic network, thus enabling the emergence of metabolic bypasses for the production of key cellular building blocks (Notebaart et al., 2014). For example, promiscuous activities of the arginine biosynthesis enzymes enabled proline production in an E. coli strain deleted in the canonical proline biosynthesis route (Itikawa et al., 1968); several pathways, each based on the promiscuous activities of different enzymes, have been shown to relieve pyridoxal phosphate auxotrophy in E. coli (Kim et al., 2010; Oberhardt et al., 2016); and adaptive evolution of E. coli harnessed enzyme promiscuity to enable growth on the non-natural feedstock 1,2-propanediol (Lee and Palsson, 2010).
The biosynthesis of isoleucine provides multiple examples of enzyme promiscuity and structural similarity to other pathways, suggesting underground metabolism as a likely origin (Jensen, 1976). Cellular production of isoleucine and valine are catalyzed by the same enzymes; in valine biosynthesis, pyruvate is self-condensed whereas in isoleucine biosynthesis, pyruvate is condensed with 2-ketobutyrate (2 KB). In most organisms, 2 KB is produced from threonine cleavage. Despite this, different biosynthetic routes, which mirror other metabolic pathways, are known to support the production of 2 KB in specific lineages (Figure 1). For example, the citramalate route mirrors the TCA cycle, where pyruvate, rather than oxaloacetate, is condensed with acetyl-CoA, and 2 KB, rather than 2-ketoglutarate, is the pathway product (Charon et al., 1974). Similarly, in several anaerobic microorganisms, propionate is converted to 2 KB in a reductive pathway that mirrors acetate conversion to pyruvate via ligation with CoA and ferredoxin-dependent carboxylation (Buchanan, 1969; Monticello et al., 1984; Eikmanns et al., 1983). In plants, methionine γ-lyase, catalyzing a reaction similar to that of cystathionine γ-lyase, is a complementary source of cellular 2 KB (Joshi and Jander, 2009; Joshi et al., 2010).
Here, we tested whether the underground metabolism of E. coli could compensate for the deletion of the canonical 2 KB production route in isoleucine biosynthesis. Under aerobic conditions, we show that a mutation in a cysteine biosynthesis enzyme decreased its activity and lowered the steady-state concentration of this amino acid; this freed an intermediate of methionine biosynthesis, which usually reacts with cysteine, to be converted to 2 KB. Under anaerobic conditions, E. coli was able to use the enzyme pyruvate formate-lyase (PFL) to generate 2 KB from propionyl-CoA and formate. Surprisingly, we found this pathway to provide a substantial fraction of the cellular 2 KB also in a wild-type (WT) strain. This study thus demonstrates the inherent ability of microorganisms to effectively exploit underground metabolism to generate patchwork pathways.
Results
A latent aerobic isoleucine biosynthesis pathway
Biosynthesis of isoleucine in E. coli starts with the deamination of threonine to give 2 KB, which is then condensed with pyruvate to produce 2-aceto-2-hydroxybutanoate. This intermediate is subsequently isomerized, reduced, dehydrated, and aminated to generate isoleucine. To generate an isoleucine auxotrophy, we constructed a strain in which the two genes encoding for threonine deaminase were deleted (ilvA and tdcB). When provided with isoleucine, this strain grew identically to the WT strain (dark red line in Figure 2). Without the addition of isoleucine, we observed no growth in the first 70 hr, but within 70–120 hr, different replicates started growing with a growth rate somewhat lower than that of the WT strain (purple lines in Figure 2).
One explanation of the observed growth without isoleucine is an increased activity of serine deaminases, which are known to accept threonine as a substrate (Cicchillo et al., 2004). We, therefore, deleted all the genes encoding for serine deaminase (sdaA, sdaB, and tdcG). We termed the resulting strain Δ5 (ΔilvA ΔtdcB ΔsdaA ΔsdaB ΔtdcG). However, strain Δ5 displayed similar behavior to the ΔilvA ΔtdcB strain, that is, while growth without isoleucine was not observed in the first 70 hr, after 70 hr, the replicates started growing (green lines in Figure 2). This strongly suggests the emergence of a latent threonine-independent isoleucine biosynthesis pathway.
To check whether this route still depends on the generation of 2 KB or rather bypasses this metabolic intermediate altogether, we constructed a strain deleted in the gene ilvC. This gene encodes a ketol-acid reductoisomerase that operates downstream of 2 KB in the isoleucine biosynthesis pathway. While cultivation of the ΔilvC strain with valine and isoleucine enabled growth similar to that of the WT strain (pink line in Figure 2; valine is required as IlvC participates also in valine biosynthesis), growth without isoleucine was not observed even after 150 hr (orange line in Figure 2). As a further confirmation, we observed that addition of 2 KB to the cultivation medium rescued the growth of the ΔilvA ΔtdcB strain and the Δ5 strain but not the ΔilvC strain (Figure 2). This confirms that 2 KB still serves as a metabolic intermediate in the underground isoleucine biosynthesis pathway.
In some microorganisms, 2 KB is produced via the citramalate pathway (Xu et al., 2004; Drevland et al., 2007; Hochuli et al., 1999; Risso et al., 2008). This route is a catalytic parallel to the first half of the TCA cycle – from citrate synthase to isocitrate dehydrogenase – where pyruvate replaces oxaloacetate, reacting with acetyl-CoA to give citramalate, which is subsequently metabolized to 2 KB. We wondered whether the enzymes of the TCA cycle, or their isozymes, catalyze the reactions of the citramalate pathway in E. coli. To check if this were the case, we performed a carbon labeling assay, feeding the WT strain, the ΔilvA ΔtdcB strain, and the Δ5 strain with either glucose-1-13C or glucose-3-13C. The expected labeling pattern of isoleucine should completely change depending on the biosynthesis route of 2 KB (Figure 3A). In the WT strain, isoleucine is expected to be roughly half labeled upon feeding with glucose-3-13C and follow a 1:2:1 pattern (unlabeled: once labeled: double labeled) with glucose-1-13C. By contrast, if 2 KB is produced via the citramalate pathway, then isoleucine would be completely unlabeled when feeding with glucose-3-13C and would roughly follow a 1:2:2:1 pattern (unlabeled: once labeled: double labeled: triple labeled) with glucose-1-13C. The labeling pattern of isoleucine in the ΔilvA ΔtdcB and the Δ5 strains is almost identical to that of the WT strain (Figure 3B, the small deviations from the expected labeling can be attributed to the ambient abundance of 13C and the shuffling of labeled carbon by the activity of the pentose phosphate pathway and the TCA cycle). This effectively rules out the activity of the citramalate pathway and indicates that the biosynthesis of 2 KB in both gene deletion strains is closely related to the natural production pathway.
2-Ketobutyrate biosynthesis from O-succinyl-L-homoserine
In several organisms, the biosynthesis of isoleucine is at least partially tied to that of the sulfur-containing amino acids cysteine and methionine. For example, in organisms in which cysteine biosynthesis depends on γ-elimination of cystathionine, 2 KB serves as a byproduct (Steegborn et al., 1999). Also, in plants, 2 KB can be produced from γ-elimination of methionine (Joshi and Jander, 2009; Joshi et al., 2010). To test whether the production of 2 KB in the Δ5 strain is related to methionine biosynthesis, we deleted metA, encoding for homoserine O-succinyltransferase, which catalyzes the first committed step of the methionine biosynthesis pathway. The Δ5 ΔmetA strain could grow with the addition of methionine and isoleucine (Figure 4, pink line), while no growth was observed without isoleucine (Figure 4, orange line). This indicates that the production of 2 KB in the Δ5 strain is indeed dependent on the methionine biosynthesis route.
At first glance, the biosynthesis of sulfur-containing amino acids in E. coli would not be expected to give rise to 2 KB. Specifically, in this bacterium, cystathionine does not undergo γ-elimination but rather β-elimination which generates homocysteine and pyruvate (Figge, 2006). Also, E. coli lacks a characterized annotated methionine γ-lyase. While an E. coli gene was previously suggested to encode for such an enzyme (Manukhov et al., 2005), the inability of the Δ5 ΔmetA strain to grow when supplemented with methionine but not isoleucine confirms that E. coli lacks methionine γ-lyase activity.
We, therefore, searched the literature for reactions that could convert intermediates of the methionine biosynthesis pathway to 2 KB. We found multiple such candidate reactions from a range of organisms (Figure 5). First, rather than catalyzing the condensation of cysteine with O-succinyl-L-homoserine, cystathionine γ-synthase (MetB) can cleave the latter intermediate to succinate and 2 KB (reaction one in Figure 5; kcat = 7.7 s−1 and KM(O-succinyl-L-homoserine)=0.33 mM for E. coli enzyme [Figge, 2006; Holbrook et al., 1990]). MetB can also condense O-succinyl-L-homoserine and homocysteine to give homolanthionine (reaction two in Figure 5; KM(homocysteine)=0.54 mM for C. glutamicum enzyme [Krömer et al., 2006]), which can be subsequently cleaved by cystathionine-β-lyase (MetC) to regenerate homocysteine and produce 2 KB (reaction three in Figure 5; e.g. kcat ≈ 180 s−1 and KM(homolanthionine)=4.5 mM for E. coli enzyme [Figge, 2006; Krömer et al., 2006; Alting et al., 1995; Dwivedi et al., 1982]). In addition, MetC from various organisms can catalyze both the β-elimination and the γ-elimination of cystathionine, giving rise to both pyruvate and 2 KB (reaction four in Figure 5; e.g. kcat(γ-elimination) ≈ 0.9 s−1 and KM(cystathionine)=0.20 mM for Streptomyces phaeochromogenes enzyme where β-elimination proceeds at ≈ 1/7 rate of γ-elimination [Alting et al., 1995; Nagasawa et al., 1984]). As MetC can also cleave cysteine to pyruvate (kcat ≈ 1.8 s−1 for E. coli enzyme [Dwivedi et al., 1982; Awano et al., 2003; Flint et al., 1996]), it is possible that it could also cleave homocysteine to produce 2 KB (reaction five in Figure 5). Finally, some MetC variants can act as methionine γ-lyase, releasing 2 KB by directly cleaving methionine (reaction six in Figure 5; e.g. kcat ≈ 0.01 s−1 for Lactococcus lactis subsp. cremoris B78 enzyme [Alting et al., 1995; Amarita et al., 2004; Irmler et al., 2008]).
To identify which of the possible biosynthesis routes is responsible for the production of 2 KB in the Δ5 strain, we performed several enzymatic experiments (Materials and methods and Figure 5—figure supplement 1). We found that MetC from E. coli did not catalyze the degradation of cystathionine, homocysteine, or methionine to 2 KB (ruling out reactions 4, 5, and 6 in Figure 5). To further confirm that MetC is not involved in 2 KB production, we deleted its encoding gene. The Δ5 ΔmetC strain could grow even without the addition of methionine (which is in line with the ΔmetC strain in the Keio collection [Baba et al., 2006]), presumably due to the existence of multiple PLP-dependent enzymes that can catalyze the MetC reaction, for example MalY, Alr, and FimE (Zdych et al., 1995; Patrick et al., 2007). More importantly, this strain could directly grow without isoleucine and without the need for adaptation (green line in Figure 4).
Supporting previous studies (Figge, 2006; Holbrook et al., 1990), we found MetB to catalyze the cleavage of O-succinyl-L-homoserine to 2 KB (reaction one in Figure 5) with kcat = 9.3 ± 0.4 s−1, KM(O-succinyl-L-homoserine)=0.60 ± 0.08 mM, and thus kcat/KM(O-succinyl-L-homoserine)=16 ± 2 mM−1 s−1 (Figure 5—figure supplement 2A). Yet, it was previously reported that the formation of 2 KB is suppressed in the presence of cysteine (Holbrook et al., 1990). We, therefore, characterized the cleavage reaction in the presence of either cysteine or homocysteine (as an alternative substrate), at physiological concentrations of ~0.3 mM (Wheldrake, 1967; Bennett et al., 2009; Guo et al., 2013) and at artificially high concentrations of 3–6 mM. We found that formation of 2 KB was indeed suppressed at physiological concentrations of cysteine, while homocysteine suppressed 2 KB formation only at high concentrations (Figure 5—figure supplement 2B–E). Nonetheless, as the deletion of metB in the Δ5 strain completely abolishes growth without isoleucine (brown line in Figure 4), it seems very likely that MetB is involved in the production of 2 KB.
Disruption of MetC or a mutation in serine acetyltransferase enable steady 2-ketobutyrate production from O-succinyl-L-homoserine
Next, we aimed to understand the genetic basis underlying 2 KB biosynthesis in the ΔilvA ΔtdcB and Δ5 strains. We sequenced the genomes of several strains isolated from the cultures growing without isoleucine. We found that in most of the sequenced strains (11 strains out of the 16 sequenced) metC was either deleted or mutated (Supplementary file 1). This is in line with the findings reported above, that is, Δ5 ΔmetC strain is able to directly grow without isoleucine and without lag time (green line in Figure 4). It therefore seems that a lower metabolic flux toward methionine biosynthesis, as expected by the deletion of metC, enhances the side reactivity of the pathway enzymes and results in a higher conversion rate of O-succinyl-L-homoserine to 2 KB. Indeed, we found that, as compared to the WT strain and the Δ5 strain, in the Δ5 ΔmetC strain the concentration of methionine was ≈3-fold lower (Figure 6A), while the concentration of O-succinyl-L-homoserine was ≈3-fold higher (Figure 6A). The side reactivity of MetB toward O-succinyl-L-homoserine cleavage therefore seems to be enhanced by the high concentration of this metabolite in the Δ5 ΔmetC strain.
Another mutated strain harbored a single mutation in the gene coding for serine acetyltransferase (CysE): Ala33Thr (Supplementary file 1). We used Multiplex Automated Genomic Engineering (MAGE [Wang et al., 2009]) to introduce this mutation into the ∆5 strain. The resulting strain (Δ5 cysE*) grew immediately (with no lag time) on a medium without isoleucine and with a growth rate identical to the WT strain (purple lines in Figure 4; note that the growth of the Δ5 cysE* strain was considerably faster than that of the Δ5 ΔmetC strain). Therefore, it seems that the mutation in cysE enabled a steady, efficient production of 2 KB.
CysE catalyzes the first committed step in the biosynthesis of cysteine (Denk and Bock, 1987). To understand the effect of the Ala33Thr mutation on enzyme activity we measured the kinetics of the purified enzyme (Materials and methods and Figure 5—figure supplements 1 and 3). Despite the mutation occurring far from the active site, we found that the apparent kcat of the enzyme decreased 2-fold, from 350 ± 30 s−1 to 170 ± 30 s−1, while the KM for acetyl-CoA increased by more than 8-fold, from 0.6 ± 0.2 mM to 5.0 ± 2.0 mM (interestingly, KM for serine changed only slightly from 0.8 ± 0.3 mM in the WT to 0.5 ± 0.2 mM in the mutant). As the concentration of acetyl-CoA in E. coli lies in the range of 0.6–0.75 mM (Bennett et al., 2009), the increase in KM for acetyl-CoA directly affects the reaction rate. Overall, the Ala33Thr mutation is expected to decrease the rate of the CysE reaction by more than 17-fold under physiological conditions.
The decreased activity of CysE halved the intracellular cysteine concentration relative to that of the WT strain and the Δ5 strain (Figure 7). This deprived MetB of its co-substrate, thus lowering the methionine biosynthesis flux and resulting in more than a 2-fold reduction in methionine concentration (Figure 6A). The lower availability of methionine led to overexpression of genes participating in its biosynthesis, increasing the transcript level of metA and metB more than 3-fold relative to the WT strain and the Δ5 strain (Figure 6B; interestingly, such increase in transcript level was not observed in the Δ5 ΔmetC strain). This higher expression, as well as the partial lifting of MetA inhibition by methionine (Born and Blanchard, 1999), led to a dramatic 30- to 60-fold increase in the concentration of O-succinyl-L-homoserine (Figure 6A). Overall, the higher concentrations of both MetB and O-succinyl-L-homoserine, coupled with the lower concentration of cysteine, resulted in the diversion of more O-succinyl-L-homoserine toward cleavage and 2 KB production. Interpreted so, the mutation in CysE – an enzyme that is not directly involved in 2 KB biosynthesis – enhances a previously negligible underground reaction for 2 KB production, thus awakening a latent isoleucine biosynthesis route (Figure 7).
Anaerobic 2-ketobutyrate biosynthesis from a reversible 2-ketobutyrate formate-lyase activity
Next, we aimed to explore underground isoleucine biosynthesis routes under anaerobic conditions (Figure 8). The evolved ΔilvA ΔtdcB cysE* strain and the constructed Δ5 cysE* strain grew anaerobically without isoleucine (Figure 8—figure supplement 1). However, we did not observe the reemergence of the metC deletion, the cysE* mutation, or any other mutation: the Δ5 strain failed to grow without the addition of isoleucine, even after 120 hr (black lines in Figure 8A). We wondered whether we could enable 2 KB production by the addition of small metabolites that E. coli might encounter in its natural habitat. We focused on propionate, a short-chain fatty acid abundant in the mammalian intestine (McNeil et al., 1978), which is known to be naturally activated to propionyl-CoA (Hesslinger et al., 1998; Liu et al., 2014). As we had previously found that PFL (encoded by pflB) supported the in vivo condensation of acetyl-CoA and formate to produce pyruvate (Zelcbuch et al., 2016), we hypothesized that this enzyme could also catalyze a 2 KB formate-lyase (KBFL) reaction, that is, accepting propionyl-CoA to produce 2 KB. We indeed found that a purified PFL catalyzes the condensation of propionyl-CoA and formate to give 2 KB with a specific activity of 0.6–1.1 μmol/min/mg, kcat of 0.9–1.7 s−1, KM(propionyl-CoA)=0.83 ± 0.6 mM and KM(formate)=69 ± 11 mM (Table 1, Materials and methods, Figure 5—figure supplement 1, and Figure 8—figure supplements 2 and 3).
We found that the Δ5 strain could grow anaerobically when propionate and formate were added to the medium instead of isoleucine (red line in Figure 8A). As formate is naturally produced by the PFL-dependent cleavage of pyruvate, we reasoned that the addition of this compound to the medium might be redundant. Indeed, we found that propionate alone, when added to the medium, supported the growth of the Δ5 strain, albeit with longer lag time and lower growth rate than with the addition of formate (purple line in Figure 8A, all replicates displayed an identical growth phenotype, thus ruling out genetic mutations). To confirm that propionate is assimilated via the activity of PFL, rather than via a reductive carboxylation route (Buchanan, 1969; Monticello et al., 1984; Eikmanns et al., 1983), we deleted pflB. When provided with isoleucine and acetate, the Δ5 ΔpflB strain was able to grow (green line in Figure 8B, acetate was added as the deletion of pflB disrupted the endogenous acetyl-CoA biosynthesis route under anaerobic conditions [Hasona et al., 2004]). However, replacement of isoleucine with propionate, with or without further addition of formate, resulted in no growth (red and purple lines in Figure 8B). This confirms that PFL is indeed responsible for 2 KB biosynthesis.
To provide a further validation that 2 KB is produced from propionate, we followed the 13C labeling of isoleucine upon addition of propionate-1-13C to the growth medium under anaerobic conditions. Surprisingly, we found that even within a WT strain, the addition of labeled propionate resulted in 30% labeled isoleucine (Figure 9), suggesting that the KBFL route operates naturally if propionate is present in the medium. When formate was also added to the medium, the fraction of labeled isoleucine increased to almost 60% (Figure 9), making the KBFL-dependent pathway the dominant isoleucine biosynthesis route within a WT strain. In the Δ5 strain, in which the canonical isoleucine biosynthesis pathways are deleted, isoleucine was almost 80% labeled upon addition of formate and labeled propionate (Figure 9). Deletion of the metA increased this labeling to 90% (Figure 9), indicating that the MetB-dependent route operates anaerobically although with a low flux. As 10% of the isoleucine was not labeled even in the Δ5 ΔmetA strain, it seems that there might be another undiscovered yet marginal route for isoleucine biosynthesis.
We wondered if the KBFL reaction is supported solely by PFL or whether TdcE, which was previously shown to catalyze the 2 KB formate-lyase reaction (Sawers et al., 1998), also contributed to the production of 2 KB. We found that purified TdcE is catalytically superior to PFL, generating 2 KB with a specific activity of 1.8–2.1 μmol/min/mg, kcat of 2.7–3.2 s−1, KM(propionyl-CoA)=0.52 ± 0.08 mM, and KM(formate)=20 ± 3 mM (Table 1, Materials and methods, Figure 5—figure supplement 1, and Figure 8—figure supplements 2 and 3). Yet, while the deletion of pflB (in which case acetate was further added to the medium) abolished isoleucine labeling, deletion of tdcE resulted in identical labeling to that of the WT strain (Figure 9). This indicates that PFL is solely responsible for the production of 2 KB from propionate and formate. Still, it could be that TdcE does not contribute to 2 KB biosynthesis only because it is not expressed under the tested conditions. Hence, we explored the consequences of overexpressing tdcE on a plasmid. We found that overexpression of tdcE in a Δ5 ΔpflB strain enabled fast growth with propionate (purple line in Figure 8C), indicating that TdcE can replace PFL both in catalyzing pyruvate cleavage to generate formate (as noted previously [Sawers et al., 1998]) and in catalyzing 2 KB production from propionyl-CoA and formate. Moreover, overexpression of tdcE in a WT strain fed with propionate-1-13C resulted in a high proportion of labeled isoleucine: more than 70% without externally added formate and ~90% with externally added formate, both are considerably higher than without tdcE overexpression (Figure 9). Overexpression of tdcE in the Δ5 ΔmetA strain resulted in 95% labeled isoleucine, the highest labeling we observed (Figure 9). Overall, these results confirm that TdcE can support in vivo 2 KB synthesis from propionyl-CoA in a comparable manner to that of PFL.
Discussion
This study serves to demonstrate the selective advantage underground metabolism may offer in nature, providing metabolic redundancy and flexibility which allows for the best use of environmental carbon sources (Notebaart et al., 2014). We demonstrated two metabolic pathways for the biosynthesis of isoleucine which are based on the activity of promiscuous enzymes. The first pathway, uncovered under aerobic conditions, is based on rerouting the metabolism of O-succinyl-L-homoserine from methionine biosynthesis toward the production of 2 KB. This represents a rather classic case of underground metabolism: a promiscuous reaction, the activity of which normally has a negligible physiological contribution, is enhanced by the rewiring of cellular metabolism under conditions which couple it to growth. Interestingly, a change in the steady-state concentration of a seemingly unrelated metabolite – cysteine – greatly enhances this promiscuous activity and enables a new biosynthetic route.
Under anaerobic conditions, we uncovered a pathway that relies on the highly expressed PFL (Neidhardt et al., 1996) to condense propionyl-CoA and formate to 2 KB. This route, unlike the aerobic cleavage of O-succinyl-L-homoserine, was highly active also in a WT strain, providing a substantial fraction of cellular isoleucine when propionate was present in the medium. As short-chain fatty acids, such as propionate, are abundant in the mammalian intestine (McNeil et al., 1978), the native habitat of E. coli, the KBFL pathway is likely to play a key role in isoleucine biosynthesis under physiological conditions. Hence, this pathway does not actually represent underground metabolism, but rather should be regarded an auxiliary biosynthesis route based on promiscuous enzyme activities. The use of the KBFL pathway is favorable over the canonical threonine cleavage route as it better uses the carbon sources available in the environment. As PFL supports the KBFL reaction under physiological conditions, we suggest renaming it ‘pyruvate formate-lyase/2-ketobutyrate formate-lyase’, in a similar manner to TdcE, which was named as ‘2-ketobutyrate formate-lyase/pyruvate formate-lyase’.
Several previous studies have suggested that, in various microbial lineages, 2 KB is derived from propionate via the activity of 2 KB:ferredoxin oxidoreductase (Buchanan, 1969; Monticello et al., 1984; Eikmanns et al., 1983; Sauer et al., 1975). The findings presented here indicate that, in at least some of these microorganisms, propionate might be converted to 2 KB via the KBFL route, where formate is produced internally via pyruvate cleavage, or even from the reduction of CO2 via formate dehydrogenase (Maia et al., 2017). Even the microorganism in which the reductive carboxylation of propionyl-CoA was first demonstrated, that is, Clostridium pasteurianum, harbors a genomic PFL (Kanehisa, 2000) and thus may also assimilate propionate into isoleucine via the KBFL pathway. Moreover, in several microorganisms, propionyl-CoA is produced from succinyl-CoA as part of anaerobic fermentation (Swick and Wood, 1960; Stams et al., 1984; Schink et al., 1987; Tholozan et al., 1994); these organisms might be able to use the KBFL route without the need for an external supply of propionate. The prevalence of both the genes and the substrates involved in the KBFL pathway indicates that this route may be widespread.
Many cellular building blocks have alternative biosynthesis pathways. Nevertheless, isoleucine seems unusually rich in the sheer number of alternative routes for its bioproduction. There appear to be at least nine routes for 2 KB production: threonine cleavage (canonical pathway), the citramalate pathway (Charon et al., 1974), the glutamate mutase pathway (Abramsky et al., 1962), cystathionine cleavage (Steegborn et al., 1999), methionine cleavage (Joshi and Jander, 2009; Joshi et al., 2010), homolanthionine cleavage (Krömer et al., 2006), O-succinyl-L-homoserine cleavage (this study), propionyl-CoA reductive carboxylation (Buchanan, 1969; Monticello et al., 1984; Eikmanns et al., 1983), and propionyl-CoA condensation with formate (this study). The large repertoire of alternative pathways makes isoleucine biosynthesis the epitome of metabolic flexibility.
Materials and methods
Strains and plasmids
Request a detailed protocolAll strains used in this study are listed in Table 2. The E. coli SIJ488 strain based upon K-12 MG1655 (Jensen et al., 2016) was used for the generation of deletion strains. SIJ488 is engineered to carry the gene deletion machinery in its genome (inducible recombinase and flippase). All gene deletions were carried out by successive rounds of λ-Red recombineering using kanamycin cassettes (FRT-PGK-gb2-neo-FRT [KAN], Gene Bridges, Heidelberg, Germany) or chloramphenicol cassettes (pKD3 [Datsenko and Wanner, 2000]) as described by Baba et al., 2006. Homologous extensions (50 bp) for the deletion cassettes were generated by PCR using oligonucleotides listed in Supplementary file 2.
For the overexpression of the 2-ketobutyrate lyase (tdcE), endogenous genes were amplified from E. coli genomic DNA using a two-step PCR (to remove cloning system relevant restriction sites [Wenk et al., 2018] – in this case a single site) using in the first round of PCR: primer pairs tdcE_HIS_A + tdcE_B and tdcE_C + tdcE_CAP_D. The corresponding PCR products were combined in a second PCR using primer pairs tdcE_HIS_A and tdcE_CAP_D. The tdcE gene was subsequently cloned into cloning vector pNivC (Wenk et al., 2018) using restriction enzymes Mph1103I and XhoI, generating pNivC-tdcE. The tdcE gene was subsequently cloned from pNiv-tdcE into expression vector pZ-ASS using enzymes EcoRI and PstI, resulting in pZASS-tdcE (Wenk et al., 2018).
Plasmids for the overexpression of MetB, MetC, and CysE (pCA24N-metB, pCA24N-metC, and pCA24N-cysE) were isolated from the ASKA collection (Kitagawa et al., 2005). The CysE Ala33Thr point mutation was generated by linear amplification of the pCA24N-cysE plasmid by PCR using following mismatch primers: CysE_A33T_fw and CysE_A33T_rv (Supplementary file 2). The template plasmid was digested by DpnI (Thermo Fisher Scientific, Darmstadt, Germany) at 37°C and E. coli DH5α was transformed with the reaction mixture to amplify the pCA24N-cysE_A33T plasmid. The mutation was confirmed by sequencing.
Cultivation conditions
Request a detailed protocolFor strain maintenance, generation of deletion strains, and for growth during cloning LB medium (1% NaCl, 1% tryptone, 0.5% yeast extract) was used. Antibiotics for selection were used at the following concentrations: chloramphenicol, 30 μg/mL; kanamycin, 50 μg/mL; ampicillin, 100 μg/mL; streptomycin, 100 μg/mL. M9 minimal media was used for growth experiments (50 mM Na2HPO4, 20 mM KH2PO4, 1 mM NaCl, 20 mM NH4Cl, 2 mM MgSO4,100 μM CaCl2, 134 μM EDTA, 13 μM FeCl3·6H2O, 6.2 μM ZnCl2, 0.76 μM CuCl2·2H2O, 0.42 μM CoCl2·2H2O, 1.62 μM H3BO3, 0.081 μM, MnCl2·4H2O) with carbon sources added to the concentrations specified in the text and figures.
For growth experiments, overnight cultures were incubated in 4 mL LB medium with no antibiotics for deletion strains, with strains having being previously plated on antibiotic media to ensure genotype. In the event a strain carried a plasmid, the relevant antibiotic was also added to the preculture medium. Before inoculation of the experiment, cultures were harvested and washed four times in M9 medium without carbon source by centrifugation (6000 rpm, 3 min) to remove residual carbon sources from the cells. Plate reader experiments were inoculated with a starting OD600 of 0.01 using the washed culture. Plate reader experiments were carried out in 96-well microtiter plates (Nunclon Delta Surface, Thermo Fisher Scientific). Each well contained 150 μL of cell culture covered with 50 μL mineral oil (Sigma-Aldrich, Taufkirchen, Germany), to avoid evaporation. An Infinite M200 Pro plate reader (Tecan) was used for incubation (37°C), shaking, and OD600 measurements. Three cycles of four shaking phases, each of 1 min were used (1. linear shaking at an amplitude of 3 mm, 2. orbital shaking at an amplitude of 3 mm, 3. linear shaking at an amplitude of 2 mm, and 4. orbital shaking at an amplitude of 2 mm). Optical density (OD 600 nm) was measured after each round of shaking (~12.5 min). Plate reader OD measurements were converted to cuvette values according to the formula ODcuvette = ODplate/0.23. Growth curves were processed in MATLAB and represent averages of (technical) triplicate measurements where the variability between triplicate measurements was less than 5% (unless explicitly stated, as in Figure 2). Within a given plate experiment, technical replicates averaged are replicates of a single genetic strain; biological replicates were generated for the Δ5 ΔmetC strain and Δ5 cysE* strains, these were engineered in parallel and were shown to behave in an identical manner.
Anaerobic growth experiments were performed in an Infinite M200 Pro plate reader (Tecan) inside a vinyl anaerobic chamber (N2 with 10% CO2, 2.5% H2, model B, Coy Laboratory Products, Grass Lake, MI). M9 mediu for growth experiments were placed for at least 24 hr in the anaerobic chamber to allow for the exchange of dissolved oxygen from the media.
Whole-genome sequencing
Request a detailed protocolGenomic DNA was extracted using the NucleoSpin Microbial DNA kit (MACHERY-NAGEL, Düren, Germany) following manufacturer’s instructions. Library construction and genome sequencing were performed by Novogene (Cambridge, United Kingdom) using the paired-end Illumina sequencing platform. Analysis of the sequencing data was performed using the breseq pipeline (Deatherage and Barrick, 2014) and the SIJ488 reference genome sequence, which is derived from the MG1655 (NC_000913, GenBank). Supplementary file 1 shows the identified mutations, as compared to the parent strains.
Multiplex automated genome engineering (MAGE)
Request a detailed protocolIntroducing genomic point mutations was achieved by using multiplex automated genome engineering (MAGE; Wang et al., 2009). A single colony of desired strain(s) transformed with pORTMAGE (Nyerges et al., 2016; Addgene catalog no. 72680) was incubated in LB medium supplemented with 100 mg L−1 of ampicillin at 30°C in a shaking incubator. To start the MAGE cycle, overnight cultures were diluted by 100 times in the same medium and cultivated to an optical density of 0.4–0.5 at 600 nm. 1 mL of each culture was transferred to sterile microcentrifuge tubes, and then transferred to 42°C thermomixer (Thermomixer C, Eppendorf, Hamburg, Germany) to express λ-Red genes by heat shock for 15 min at 1000 rpm. After induction, cells were quickly chilled on ice for at least 15 min, and then made electrocompetent by washing three times with ice-cold ddH2O. 40 uL of electrocompetent cell was mixed with 2 uL of 50 uM of oligomer stock solution and the final volume of the suspension was adjusted to 50 uL. The oligomer used to introduce the mutation was CysE_A33T_MAGE (Supplementary file 2). Electroporation was done on Gene Pulser XCell (Bio-Rad) set to 1.8 kV, 25 μF capacitance, and 200 Ω resistance for 1 mm gap cuvette. Immediately after electroporation, 1 mL of LB was added to cuvette and the electroporation mixtures in LB were transferred to sterile culture tubes and cultured with shaking at 30°C, 240 rpm for 1 hr to allow for recovery. After recovery, 2 mL of LB medium supplemented with ampicillin was added and then further incubated in the same condition. When the culture reached an OD600 of 0.4–0.5, cells were either subjected to additional MAGE cycles or analyzed for genotype via PCR and sequencing.
Isotopic-labeling of proteinogenic amino acids
Request a detailed protocolStationary isotope tracing experiments using 13C-1-labeled glucose, 13C-3-labeled glucose, and labeled sodium propionate (Sigma-Aldrich) were performed in order to understand the metabolic activities responsible for isoleucine biosynthesis. All experiments were performed in duplicate. Strains were grown in M9 with 10 mM labeled glucose as sole carbon source in the case of labeled glucose experiments, and 10 mM unlabeled glucose with 5 mM sodium propionate-1-13C in the case of the anaerobic experiments (extra carbon sources added where relevant at concentrations indicated). Anaerobic experiments were performed also with an anaerobic pre-culture in M9 medium. The equivalent of 1 mL of culture at OD600one was harvested by centrifugation and washed twice in water. Cellular biomass was hydrolyzed by boiling in acid (95°C in 6 M HCl for 24 hr [You et al., 2012]). After drying the samples at 95°C, the samples were resuspended in ultra-pure water and the masses of amino acids analyzed with UPLC–ESI–MS as previously described (Giavalisco et al., 2011). Chromatography was performed with a Waters Acquity UPLC system (Waters), using an HSS T3 C18 reversed phase column (100 × 2.1 mm2, 1.8 μm; Waters). 0.1% formic acid in H2O (A) and 0.1% formic acid in acetonitrile (B) were the mobile phases. Flow rate was 0.4 mL/min and the gradient was: 0 to 1 min – 99% A; 1 to 5 min – linear gradient from 99% A to 82%; 5 to 6 min – linear gradient from 82% A to 1% A; 6 to 8 min – kept at 1% A; 8 to 8.5 min – linear gradient to 99% A; 8.5 to 11 min – re-equilibrate. Mass spectra were acquired using an Exactive mass spectrometer (Thermo Fisher Scientific) in positive ionization mode, with a scanning range of 50.0 to 300.0 m/z. Spectra were recorded during the first 5 min of the LC gradients. Retention times for amino acids under these conditions were determined by analyzing amino-acid standards (Sigma-Aldrich) under the same conditions. Data analysis was performed using Xcalibur (Thermo Fisher Scientific).
Enzyme assays for MetB, MetC, and CysE
Request a detailed protocolE. coli BL21 (DE3) was transformed with the overexpression plasmids pCA24N-metB, pCA24N-metC, pCA24N-cysE, or pCA24N-cysE_A33T, and grown on LB agar containing 34 µg/mL chloramphenicol. An expression culture in TB containing 34 µg/mL chloramphenicol was inoculated from the plate and grown at 37°C while shaking at 110 rpm until an OD600 of roughly 0.9 was reached. The culture was shaken for another 4 hr at 37°C or overnight at 20°C and 110 rpm for gene expression. The culture was harvested by centrifugation at 4°C and 5000 g. Alternatively, the pellet was stored at −20°C. The pellet was resuspended in Buffer A (50 mM Tris-HCl pH 7.9, 500 mM NaCl) containing 10 µg/mL DNaseI and 5 mM MgCl2. The cells were lysed by ultrasonication and the lysate was cleared by ultracentrifugation at 50,000 g and 4°C for 45 min followed by filtration through a 0.45 µm syringe filter. The lysate was loaded onto a 1 mL HisTrap FF column (GE Healthcare, Freiburg, Germany) equilibrated in Buffer A. Unspecifically bound protein was washed with Buffer A containing 50 mM imidazole. The protein was eluted with Buffer A containing 500 mM imidazole and subsequently desalted over a HiTrap 5 mL desalting column (GE Healthcare) using 20 mM Tris-HCl pH 7.9, 50 mM NaCl. The protein was concentrated using Amicon Ultra-4 centrifugal filters (Merck Millipore, Darmstadt, Germany). The protein was stored in 50% (v/v) glycerol at −20°C.
The MetB elimination reaction (red arrow, reaction one in Figure 5) was kinetically characterized using a spectrophotometric assay. The assays were performed at 30°C in 10 mm quartz cuvettes (Hellma Analytics, Germany) on a Cary-60 UV/Vis spectrometer (Agilent Technologies Inc, Santa Clara, CA). In the elimination reaction, MetB cleaves O-succinyl-L-homoserine to succinate and 2-ketobutyrate. The latter can be reduced by the coupling enzyme lactate dehydrogenase using NADH as reducing agent. The NADH consumption was followed at 340 nm (εNADH = 6.22 mM−1 cm−1). The reaction was started by adding O-succinyl-L-homoserine to 22 nM MetB in 200 mM Tris-HCl pH 8.0 containing ~1 U of lactate dehydrogenase (Sigma-Aldrich Chemie GmbH, Germany) and 300 µM NADH. Kinetic parameters were determined from a Michaelis-Menten fit of 18 data points.
In the presence of O-succinyl-L-homoserine and cysteine or homocysteine, MetB can catalyze the replacement reaction yielding succinate and cystathionine or homolanthionine, respectively (MetB black arrow and reaction two in Figure 5). However, the presence of cysteine or homocysteine does not necessarily exclude the elimination reaction activity of the MetB enzyme (reaction one in Figure 5). To determine the reaction route, MetB was incubated with O-succinyl-L-homoserine and cysteine or homocysteine in separate experiments. To investigate the MetB reaction route at physiological cysteine or homocysteine concentrations, we ran an assay with 100 nM of MetB in 200 mM Tris-HCl pH 8.0 and 50 µM pyridoxal-5’-phosphate in the presence of 6 mM O-succinyl-L-homoserine and 300 µM of either cysteine or homocysteine. The assay mixture was incubated at 30°C and samples were quenched with 5% formic acid at specific time points. To investigate the MetB reaction route at high concentrations of cysteine and homocysteine, 90 nM of MetB was assayed in 200 mM Tris-HCl pH 8.0 in the presence of 6 mM O-succinyl-L-homoserine and either 3 mM L-cysteine or 6 mM L-homocysteine. The assay mixture was incubated at 30°C and samples were quenched with 10% formic acid at specific time points. In both experiments, the formation of 2 KB and succinate was quantitatively determined by HPLC-MS/MS analysis. The levels of succinate shown in Figure 5—figure supplement 1B–E are the levels of succinate generated from the replacement reaction of MetB exclusively; the succinate generated by the elimination reaction has already been taken into account by subtracting the concentration of 2 KB in the sample (elimination reaction produces an equimolar amount of 2 KB and succinate).
To test whether MetC can produce 2 KB, we incubated 24 µM of MetC and 10 mM of homocysteine or 87 µM MetC and 33 mM of methionine in 50 mM Tris-HCl pH 8.5 at 30°C (reaction 5 or 6 in Figure 5). No 2 KB generation could be detected when following the NADH consumption of the coupling enzyme lactate dehydrogenase (~1 U of lactate dehydrogenase and 300 µM NADH). Some MetC variants from other organisms have also been reported to catalyze both the β-elimination (canonical E. coli) and the γ-elimination of cystathionine resulting in pyruvate and 2 KB, respectively. In this case, as both these products are accepted by the coupling enzyme lactate dehydrogenase (Kim and Whitesides, 1988), the reaction route could not be determined spectrophotometrically. Instead, 20 nM of MetC was incubated with 3 mM cystathionine in 50 mM Tris-HCl pH 8.5. The assay mixture was incubated at 30°C and samples were quenched with 10% formic acid at specific time points. Formation of 2 KB and pyruvate was quantitatively determined using HPLC-MS/MS analysis. Only pyruvate could be detected in these samples indicating no γ-elimination of cystathionine by the E. coli MetC enzyme.
The activity of CysE and its variant CysE Ala33Thr were assayed by photospectrometrically following the reaction of released CoA with 5,5’-dithiobis-(2-nitrobenzoic acid) (DTNB) at 412 nm using an extinction coefficient of 14.4 mM−1 cm−1. The assays were performed at 30°C in 10 mm quartz cuvettes (Hellma Analytics) on a Cary-60 UV/Vis spectrometer (Agilent Technologies). The reaction was started by adding the CysE variant to 0.2 mM DTNB in 100 mM potassium phosphate buffer pH 8. The dependence of the turnover number on acetyl-CoA concentration (synthesized as previously described [Peter et al., 2016]) was measured with 0.3 nM CysE or 1.3 nM CysE A33T at a saturating L-serine concentration of 20 mM. The dependence of the turnover number on L-serine concentration was measured with 1.0 nM CysE or 2.0 nM CysE A33T at a non-saturating acetyl-CoA concentration of 0.6 mM.
Succinate, pyruvate, and 2-ketobutyrate in the enzymatic assays were quantitatively determined using LC-MS/MS based on external calibration curves. The calibration curves were generated with at least seven concentrations within the linear range of the analyte in the matrix of the assay samples. The compounds were separated over a Kinetex EVO C18 chromatography column (50 × 2.1 mm2; 3.5 µm, 100 Å, Phenomenex) equipped with a 20 × 2.1 mm2 guard column of similar specificity on an Agilent Infinity II 1290 HPLC system. The column was heated to 55°C and 5 µL of sample or standard were injected. A constant flow rate of 0.2 mL/min was applied. Mobile phase A (0.1% formic acid in water [Honeywell, Morristown, New Jersey]) and mobile phase B (0.1% formic acid in acetonitrile [Honeywell]) were mixed according to the following mobile phase profile: 0–1.6 min constant at 0% B; 1.6–5 min from 0% to 95% B; 5–6 min constant at 95% B; 6–6.5 min from 95% to 0% B; 6.5–10 min constant at 0% B. An Agilent 6495 ion funnel mass spectrometer equipped with an electrospray ionization source in negative mode was used for analysis under the following conditions: ESI spray voltage 1500 V, sheath gas 200° C at 12 L/min, nebulizer pressure 30 psig, and drying gas 180° C at 11 L/min. Mass transition and retention time of the compounds were compared to standards for identification. The MassHunter software (Agilent, Santa Clara, CA, USA) was used for chromatogram integration.
Anaerobic production of 2-ketobutyrate from PFL (i.e. PflB) and TdcE
Request a detailed protocolPropionyl-CoA was synthesized and purified as described previously (Peter et al., 2016). The dried powder was kept at −20°C and dissolved in 50 mM acetate (pH 4.5) before use. The concentration was determined by two independent methods. First, a NanoDrop 2000 Spectrophotometer (Thermo Fisher Scientific) was used with the extinction coefficient of saturated CoA-esters at 260 nm (ε260nm = 16.4 mM−1 cm−1). Second, the depletion of the compound when incubated with the enzyme PduP (CoA-acylating propionaldehyde dehydrogenase; Zarzycki et al., 2017), following the oxidation of NADH at 340 nm (ε340nm = 6.22 mM−1 cm−1). Both methods gave the same results within a 5% error.
Overexpression plasmids for PFL-AE (PFL-activating enzyme), PFL, and TdcE were obtained from the ASKA collection. PFL-AE, PFL, and TdcE were produced in E. coli BL21(DE3). 1 L terrific broth (TB) containing 34 µg/mL chloramphenicol was inoculated with freshly transformed cells and incubated at 37°C. After reaching an OD600 of 0.8, the expression was induced by adding IPTG to a final concentration of 0.5 mM and the incubation temperature was lowered to 25°C. Cells were harvested after 16 hr by centrifugation (4500 g, 10 min) and resuspended in buffer A (50 mM HEPES-KOH pH 7.8, 500 mM KCl) containing 10 µg/mL DNaseI and 5 mM MgCl2. If not used immediately, cell pellets were flash-frozen in liquid nitrogen and stored at −20°C. The cell lysate obtained by sonication was clarified by centrifugation 75,000 g at 4°C for 45 min. The supernatant was filtered through a 0.4 µm syringe tip filter (Sarstedt). Ni-affinity purification was performed with an Äkta FPLC system from GE Healthcare. The filtered soluble lysate was loaded onto a 1 mL Ni-Sepharose Fast Flow column (HisTrap FF, GE Healthcare) that had been equilibrated with 10 mL buffer A. After washing with 20 mL 85% buffer A, 15% buffer B (50 mM HEPES-KOH pH 7.8, 500 mM KCl, 500 mM imidazole), the protein was eluted with 100% buffer B. Fractions containing purified protein were pooled and the buffer was exchanged to storage buffer (50 mM HEPES-KOH pH 7.8, 150 mM KCl) using a desalting column (HiTrap, GE Healthcare). Proteins were concentrated by ultrafiltration (Amicon Ultra). Protein concentration was determined on a NanoDrop 2000 Spectrophotometer (Thermo Fisher Scientific) using the extinction coefficient at 280 nm, as calculated by ProtParam (Gasteiger, 2005). Enzyme purity was confirmed using SDS-PAGE. The purified proteins were flash-frozen in liquid nitrogen and stored at −80°C.
PFL-AE was expressed and purified as PFL and TdcE. However, after the HisTrap purification step, the protein was transferred to an anaerobic glovebox. The reconstitution of the iron-sulfur cluster was performed following the protocol of Broderick et al. with a few adaptations (Byer et al., 2018). Using a PD-10 desalting column (GE Healthcare), the buffer was exchanged to an anaerobic reconstitution buffer (50 mM Tris-SO4 pH 7.5, 150 mM KCl, 5 mM dithiothreitol [DTT]). The DTT was added fresh to the buffer right before use. After buffer exchange the protein was diluted with reconstitution buffer to a concentration of ~3 mg/mL (~100 µM) and transferred to a glass vial with a magnetic stirring bar. The vial was placed on ice on a magnetic stirrer. A 10 mM stock solution of FeCl3 was added in 10 aliquots over the course of 1 hr to a final concentration of 600 µM. Then, a 50 mM stock of Na2S·(H2O)9 was added in 10 aliquots over the course of 1 hr to a final concentration of 600 µM. This mixture was incubated on ice with stirring for 6 hr. Then the protein was desalted with a PD-10 column (equilibrated with reconstitution buffer), and concentrated on an Amicon to 30 mg/mL. UV-Vis spectra were recorded before and after the reconstitution with a Cary 4000 UV-Vis spectrometer (Agilent Technologies) using a FiberMate2 Fiber Optic Coupler (Harrick Scientific Products) with a path length of 1 cm (Figure 8—figure supplement 2). The reconstituted PFL-AE was transferred in aliquots to rubber stoppered HPLC vials, flash frozen in liquid nitrogen and stored at −80°C.
All anaerobic steps were performed in an anaerobic glovebox (Coy Laboratories) under an N2 atmosphere containing 3% to 3.5% H2. Residual O2 was removed by palladium catalysts. O2 concentration was monitored and maintained below 5 ppm at all times. All stocks were prepared under anaerobic conditions or prepared outside of the glovebox and equilibrated for at least 3 hr inside the glovebox. Dithionite and DTT stocks were freshly prepared under anaerobic conditions each day.
To activate PFL, 100 mM Tris-HCl pH 7.5, 50 mM KCl, 10 mM DTT, 50 mM sodium oxamate, 125 µM SAM, and 25 µM PFL-AE were incubated for 10 min at 30°C. Then, 25 µM PFL and 62.5 µM sodium dithionite were added. After 30 min at 30°C, the activation mix was placed on ice.
To activate TdcE, 100 mM Tris-HCl pH 7.5, 50 mM KCl, 10 mM DTT, 50 mM sodium oxamate, 50 µM SAM, and 10 µM PFL-AE were incubated for 10 min at 30°C. Then, 10 µM TdcE and 25 µM sodium dithionite were added. After 30 min at 30°C, the activation mix was placed on ice.
Glycyl radical formation was confirmed by exposing the activation mix to oxygen followed by SDS-PAGE analysis. The glycyl radical rapidly reacts with oxygen, cleaving off a ~3.5 kDa fragment from the C-terminus. PFL and TdcE were activated to a similar degree, as demonstrated in Figure 8—figure supplement 2.
Activation and stability of the activation mix was tested by the pyruvate cleavage assay (Figure 8—figure supplement 2). For this, 100 mM Tris-HCl pH 7.5, 50 mM KCl, 10 mM malate, 1 mM NAD, 20 mM pyruvate, 250 µM CoA, 4.5 U/mL citrate synthase from the porcine heart (Sigma-Aldrich), 27.5 U/mL malic dehydrogenase from the porcine heart (Sigma-Aldrich) were mixed and the reaction initiated by adding the activation mix in 100-fold dilution. The reduction of NAD was followed at 340 nm.
2-Ketobutyrate formation kinetics were measured in a coupled enzyme assay: L-lactate dehydrogenase reduces 2-ketobutyrate to 2-hydroxybutyrate with the consumption of NADH. Absorption at 340 nm was recorded with a Cary 4000 UV-Vis spectrometer using a FiberMate2 Fiber Optic Coupler with a path length of 1 cm. Each substrate concentration was measured in triplicate. One activation mix was used up to 1.5 hr storage on ice, then a fresh activation mix was prepared for kinetic measurements. PFL and TdcE kinetic parameters take the half-of-the-sites activity (Unkrig et al., 1989) into account, that is, the dimer mass was used to calculate the protein concentration. Because the exact active site concentration (i.e. the glycyl radical content after activation) was not determined, Vmax is reported as specific activity. For Michaelis-Menten parameter calculation and fitting, GraphPad Prism 8 was used. For PFL with propionyl-CoA as substrate the data was fit to the substrate inhibition equation.
Propionyl-CoA kinetics: 100 mM Tris-HCl pH 7.5, 50 mM KCl, 500 mM sodium formate, 200 µM NADH, 25 U/mL L-lactate dehydrogenase from bovine heart (LDH; Sigma-Aldrich). The reaction was started by adding an activation mix in 100-fold dilution. Thus, the final concentration was 250 nM PFL and 100 nM TdcE, respectively. For PFL, the propionyl-CoA concentrations were 100 µM, 250 µM, 500 µM, 1 mM, 2.5 mM and for TdcE, 150 µM, 350 µM, 550 µM, 800 µM, 1.6 mM, 1.9 mM. In both cases, there was no activity in the absence of propionyl-CoA.
Formate kinetics: 100 mM Tris-HCl pH 7.5, 50 mM KCl, 200 µM NADH, 25 U/mL LDH. Propionyl-CoA was added to a final concentration of 800 µM (PFL) or 1.9 mM (TdcE). The reaction was started by adding an activation mix in 100-fold dilution. Thus, the final concentration was 250 nM PFL and 100 nM TdcE, respectively. For PFL, the formate concentrations were 12 mM, 40 mM, 100 mM, 250 mM, 500 mM, and for TdcE, 5 mM, 12 mM, 40 mM, 100 mM, 500 mM. In both cases, there was no activity in the absence of formate.
Transcript level analysis by reverse transcriptase quantitative PCR
Request a detailed protocolTo determine mRNA levels, total RNA was extracted from three biological replicates from cells in exponential phase growing on M9 minimal medium with 10 mM glucose (supplemented with 2 mM isoleucine in case of the Δ5 strain). Total RNA was purified using the RNeasy Mini Kit (Qiagen, Hilden, Germany) as instructed by the manufacturer. In brief,~2.5×108 cells (0.5 mL of OD600 0.5) were mixed with 2 volumes of RNAprotect Bacteria Reagent (Qiagen) and pelleted, followed by enzymatic lysis, on-column removal of genomic DNA with RNase-free DNase (Qiagen) and spin-column-based purification of RNA. Concentration and integrity of the isolated RNA were determined by NanoDrop and gel electrophoresis. cDNA was synthesized via reverse transcription of 500 ng RNA with the qScript cDNA Synthesis Kit (QuantaBio, Beverly, MA). Quantitative real-time PCR was performed two times in technical triplicates using the Maxima SYBR Green/ROX qPCR Master Mix (Thermo Fisher Scientific). An input corresponding to 25 pg total RNA/cDNA was used per reaction. Non-specific amplification products were excluded by melting curve analysis. The gene encoding 16S rRNA (rrsA) was chosen as a well-established reference transcript for expression normalization (Zhou et al., 2011). Two primer pairs for amplification of metA and metB, respectively, were used for qPCR (Supplementary file 2). Equal amplification efficiencies between the primers for the genes of interest and the reference gene were assumed. Differences in expression levels were calculated according to the 2-∆∆Ct method (Livak and Schmittgen, 2001). Reported data represents the average of 2-∆∆Ct values that were calculated for each sample individually relative to the average of all biological WT replicate ∆Ct(Ct(GOI)-Ct(rrsA)) values. Negative control assays with the direct input of RNA (without previous reverse transcription) confirmed that the observed fold change of metA and metB in the Δ5 cysE* sample was not due to an increase in genomic DNA contamination, that is, it accounted for less than ~10% of the signal (∆Ct between +RT/–RT samples >3).
Cellular concentrations of methionine and succinyl-homoserine
Request a detailed protocolTo assess the relative intracellular concentration changes of succinyl-homoserine and methionine among different strains, we cultivated the strains in four independent replicates of 20 mL media in 100 mL Erlenmeyer flasks (n = 4) to mid-exponential phase (OD600 approx. 0.4–0.8). After measuring OD600, 12 mL culture were immediately harvested and quenched with 36 mL precooled 60% methanol (v/v). Cells were spun down by centrifugation (Beckman Allegra 25R refrigerated centrifuge): 5000 g, −10°C, 10 min, and stored at −80°C before extraction. On ice, the cells were added to precooled extraction fluid (10 mM Tris-HCl, 1 mM EDTA, 50% MeOH (v/v), pH 7.0) and chloroform 0.5 mL of each, briefly vortexed and shaken at 4°C for 2 hr for extraction. After centrifugation (5000 g, −10°C, 10 min), the upper phase of the extracts were filtered (Fisher syringe filter, PTFE, 0.2 µM, 13 mm diameter), and stored at −80°C until further processing.
Quantitative determination was performed using an HRES-LC-MS. The chromatographic separation was performed on a Vanquish HPLC system using a ZicHILIC SeQuant column (150 × 2.1 mm2, 3.5 μm particle size, 100 Å pore size) connected to a ZicHILIC guard column (20 × 2.1 mm2, 5 μm particle size) (Merck KgAA) with a constant flow rate of 0.3 mL/min with mobile phase A being 0.1% formic acid in 99:1 water:acetonitrile (Honeywell) and phase B being 0.1% formic acid 99:1 water:acetonitrile (Honeywell) at 45° C.
The injection volume was 2 µL. The mobile phase profile consisted of the following steps and linear gradients: 0–1 min constant at 90% B; 1–8 min from 90% to 20% B; 8–9 min constant at 20% B; 9–9.1 min from 20% to 90% B; 9.1–11 min constant at 90% B. An ID-X Orbitrap mass spectrometer equipped with a HESI electrospray ion source (Thermo Fisher Scientific) was used at the following conditions: ESI spray voltage 3500 V, sheath gas at 50 AU, auxiliary gas at 10 AU and sweep gas at 1 AU, with the vaporizor temperature being 350°C and the ion transfer tube temperature being 325°C. Compounds were identified based on their accurate mass and retention time compared to standards. Chromatograms were integrated using Xcalibur software (Thermo Fisher Scientific). Absolute concentrations were calculated based on an external calibration curve.
Relative cellular concentrations of cysteine and homocysteine
Request a detailed protocolThe concentration of cysteine could not be determined using the above method due to its instability during the measurement process. Therefore, we used another method to assess the relative intercellular concentration of these amino acids.
To assess the relative intracellular concentration changes of cysteine, homocysteine, and 2 KB among different strains, we cultivated the strains in four independent replicates of 4 mL media in glass test tubes (n = 4) to mid-exponential phase (OD600 approx. 0.4–0.8). About 109 cells (equivalent of 1 mL of OD600 = 1 culture) were collected on 0.45 μm Durapore membrane filters (Merck Millipore, Ireland). After washing the filters with 1 mL of fresh medium in under 10 s, the filters with cells were placed in 5 mL of −20°C 40:40:20 (v/v/v) acetonitrile/methanol/water overnight. 4 mL of debris-free extract per sample was dried by vacuum centrifugation at 35°C overnight. The dried metabolites were stored at −20°C until further processing.
Metabolites were methoxyaminated and trimethylsilylated before gas chromatography-atmospheric pressure chemical ionization-quadrupole time of flight mass spectrometry (GC/APCI-qTOF-MS) as described earlier (Kopka et al., 2017). GC/APCI-qTOF-MS analyses were performed with an Agilent 7890B gas chromatograph (Agilent Technologies Deutschland GmbH, Germany) hyphenated to a Bruker Impact II mass spectrometer (Bruker Daltonik GmbH, Bremen, Germany). The machine parameters and analysis procedures were exactly as detailed in Kopka et al., 2017. Data mining was performed with Profile Analysis Version 2.2 software (Bruker Daltonik GmbH).
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files. Source data files have been provided for Figures 2 and 7 as well as for the metabolomic analysis.
References
-
The formation of isoleucine from beta-methylaspartic acid in Escherichia coli WThe Journal of Biological Chemistry 240:PC265–PC266.
-
Identification and functional analysis of the gene encoding methionine-gamma-lyase in Brevibacterium linensApplied and Environmental Microbiology 70:7348–7354.https://doi.org/10.1128/AEM.70.12.7348-7354.2004
-
Effect of cysteine desulfhydrase gene disruption on L-cysteine overproduction in Escherichia coliApplied Microbiology and Biotechnology 62:239–243.https://doi.org/10.1007/s00253-003-1262-2
-
Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the keio collectionMolecular Systems Biology 2:2006–2008.https://doi.org/10.1038/msb4100050
-
Absolute metabolite concentrations and implied enzyme active site occupancy in Escherichia coliNature Chemical Biology 5:593–599.https://doi.org/10.1038/nchembio.186
-
Role of ferredoxin in the synthesis of alpha-ketobutyrate from propionyl coenzyme A and carbon dioxide by enzymes from photosynthetic and nonphotosynthetic BacteriaThe Journal of Biological Chemistry 244:4218–4223.
-
Escherichia coli L-serine deaminase requires a [4Fe-4S] cluster in catalysisJournal of Biological Chemistry 279:32418–32425.https://doi.org/10.1074/jbc.M404381200
-
Identification of mutations in laboratory-evolved microbes from next-generation sequencing data using breseqMethods in Molecular Biology 1151:165–188.https://doi.org/10.1007/978-1-4939-0554-6_12
-
Enzymology and evolution of the pyruvate pathway to 2-oxobutyrate in Methanocaldococcus jannaschiiJournal of Bacteriology 189:4391–4400.https://doi.org/10.1128/JB.00166-07
-
Propionate assimilation by methanogenic BacteriaArchives of Microbiology 136:106–110.https://doi.org/10.1007/BF00404782
-
BookMethionine Biosynthesis in Escherichia coli and Corynebacterium glutamicumIn: Wendisch VF, editors. Amino Acid Biosynthesis - Pathways, Regulation and Metabolic Engineering. Springer. pp. 163–193.https://doi.org/10.1007/978-3-540-48596-4
-
Studies on the synthesis of the Fe-S cluster of dihydroxy-acid dehydratase in Escherichia coli crude extract isolation of O-acetylserine sulfhydrylases A and B and beta-cystathionase based on their ability to mobilize sulfur from cysteine and to participate in Fe-S cluster synthesisThe Journal of Biological Chemistry 271:16053–16067.https://doi.org/10.1074/jbc.271.27.16053
-
BookThe Proteomics Protocols HandbookIn: Walker J. M, editors. Protein Identification and Analysis Tools on the ExPASy Server. Humana Press. pp. 571–607.https://doi.org/10.1385/1-59259-890-0:571
-
ECMDB: the E. coli Metabolome DatabaseNucleic Acids Research 41:D625–D630.https://doi.org/10.1093/nar/gks992
-
Amino acid biosynthesis in the halophilic archaeon Haloarcula hispanicaJournal of Bacteriology 181:3226–3237.https://doi.org/10.1128/JB.181.10.3226-3237.1999
-
Cloning and characterization of two Lactobacillus casei genes encoding a cystathionine lyaseApplied and Environmental Microbiology 74:99–106.https://doi.org/10.1128/AEM.00745-07
-
Enzymic basis for a genetic suppression: accumulation and deacylation of N-acetylglutamic γ-semialdehyde in enterobacterial mutantsBiochimica Et Biophysica Acta (BBA) - Enzymology 159:547–550.https://doi.org/10.1016/0005-2744(68)90142-3
-
Enzyme recruitment in evolution of new functionAnnual Review of Microbiology 30:409–425.https://doi.org/10.1146/annurev.mi.30.100176.002205
-
KEGG: kyoto encyclopedia of genes and genomesNucleic Acids Research 28:27–30.https://doi.org/10.1093/nar/28.1.27
-
Enzyme promiscuity: a mechanistic and evolutionary perspectiveAnnual Review of Biochemistry 79:471–505.https://doi.org/10.1146/annurev-biochem-030409-143718
-
L-Lactate dehydrogenase: substrate specificity and use as a catalyst in the synthesis of homochiral 2-hydroxy acidsJournal of the American Chemical Society 110:2959–2964.https://doi.org/10.1021/ja00217a044
-
On the origin of metabolic pathwaysJournal of Molecular Evolution 49:424–431.https://doi.org/10.1007/PL00006565
-
Adaptive evolution of Escherichia coli K-12 MG1655 during growth on a nonnative carbon source, L-1,2-propanediolApplied and Environmental Microbiology 76:4158–4168.https://doi.org/10.1128/AEM.00373-10
-
Acs is essential for propionate utilization in Escherichia coliBiochemical and Biophysical Research Communications 449:272–277.https://doi.org/10.1016/j.bbrc.2014.05.015
-
Isoleucine synthesis by Clostridium sporogenes from propionate or alpha-methylbutyrateJournal of General Microbiology 130:309–318.https://doi.org/10.1099/00221287-130-2-309
-
Cystathionine gamma-lyase of Streptomyces phaeochromogenes the occurrence of cystathionine gamma-lyase in filamentous Bacteria and its purification and characterizationThe Journal of Biological Chemistry 259:10393–10403.
-
Metabolite-Enzyme coevolution: from single enzymes to metabolic pathways and networksAnnual Review of Biochemistry 87:187–216.https://doi.org/10.1146/annurev-biochem-062917-012023
-
Multicopy suppression underpins metabolic evolvabilityMolecular Biology and Evolution 24:2716–2722.https://doi.org/10.1093/molbev/msm204
-
Elucidation of an alternate isoleucine biosynthesis pathway in Geobacter sulfurreducensJournal of Bacteriology 190:2266–2274.https://doi.org/10.1128/JB.01841-07
-
Amino acid biosynthesis in mixed rumen culturesBiochemical Journal 150:357–372.https://doi.org/10.1042/bj1500357
-
The glycyl radical enzyme TdcE can replace pyruvate formate-lyase in glucose fermentationJournal of Bacteriology 180:3509–3516.https://doi.org/10.1128/JB.180.14.3509-3516.1998
-
Pathway of propionate formation from ethanol in Pelobacter propionicusArchives of Microbiology 147:321–327.https://doi.org/10.1007/BF00406127
-
Pathway of propionate formation in Desulfobulbus propionicusArchives of Microbiology 139:167–173.https://doi.org/10.1007/BF00401994
-
Kinetics and inhibition of recombinant human cystathionine gamma-lyase toward the rational control of transsulfurationThe Journal of Biological Chemistry 274:12675–12684.https://doi.org/10.1074/jbc.274.18.12675
-
Metabolic pathway to propionate of Pectinatus frisingensis, a strictly anaerobic beer-spoilage bacteriumArchives of Microbiology 162:401–408.https://doi.org/10.1007/BF00282104
-
An engineering approach for rewiring microbial metabolismMethods in Enzymology 608:329–367.https://doi.org/10.1016/bs.mie.2018.04.026
-
Metabolic pathway confirmation and discovery through (13)C-labeling of proteinogenic amino acidsJournal of Visualized Experiments 26:e3583.https://doi.org/10.3791/3583
-
MalY of Escherichia coli is an enzyme with the activity of a beta C-S lyase (cystathionase)Journal of Bacteriology 177:5035–5039.https://doi.org/10.1128/JB.177.17.5035-5039.1995
Article and author information
Author details
Funding
Max Planck Society
- Charles AR Cotton
- Iria Bernhardsgrütter
- Hai He
- Simon Burgener
- Luca Schulz
- Nicole Paczia
- Beau Dronsella
- Alexander Erban
- Stepan Toman
- Marian Dempfle
- Alberto De Maria
- Joachim Kopka
- Steffen N Lindner
- Tobias J Erb
- Arren Bar-Even
Federal Ministry of Education and Research (031B0194 (FormatPlant))
- Charles AR Cotton
- Simon Burgener
- Luca Schulz
- Stepan Toman
- Steffen N Lindner
- Tobias J Erb
- Arren Bar-Even
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This study was funded by the Max Planck Society and by the German Federal Ministry of Education amd Research, Grant no. 031B0194 (project FormatPlant) .
Copyright
© 2020, Cotton et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 3,992
- views
-
- 496
- downloads
-
- 25
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
Dynamic conformational and structural changes in proteins and protein complexes play a central and ubiquitous role in the regulation of protein function, yet it is very challenging to study these changes, especially for large protein complexes, under physiological conditions. Here, we introduce a novel isobaric crosslinker, Qlinker, for studying conformational and structural changes in proteins and protein complexes using quantitative crosslinking mass spectrometry. Qlinkers are small and simple, amine-reactive molecules with an optimal extended distance of ~10 Å, which use MS2 reporter ions for relative quantification of Qlinker-modified peptides derived from different samples. We synthesized the 2-plex Q2linker and showed that the Q2linker can provide quantitative crosslinking data that pinpoints key conformational and structural changes in biosensors, binary and ternary complexes composed of the general transcription factors TBP, TFIIA, and TFIIB, and RNA polymerase II complexes.
-
- Biochemistry and Chemical Biology
- Stem Cells and Regenerative Medicine
Human induced pluripotent stem cells (hiPSCs) have great potential to be used as alternatives to embryonic stem cells (hESCs) in regenerative medicine and disease modelling. In this study, we characterise the proteomes of multiple hiPSC and hESC lines derived from independent donors and find that while they express a near-identical set of proteins, they show consistent quantitative differences in the abundance of a subset of proteins. hiPSCs have increased total protein content, while maintaining a comparable cell cycle profile to hESCs, with increased abundance of cytoplasmic and mitochondrial proteins required to sustain high growth rates, including nutrient transporters and metabolic proteins. Prominent changes detected in proteins involved in mitochondrial metabolism correlated with enhanced mitochondrial potential, shown using high-resolution respirometry. hiPSCs also produced higher levels of secreted proteins, including growth factors and proteins involved in the inhibition of the immune system. The data indicate that reprogramming of fibroblasts to hiPSCs produces important differences in cytoplasmic and mitochondrial proteins compared to hESCs, with consequences affecting growth and metabolism. This study improves our understanding of the molecular differences between hiPSCs and hESCs, with implications for potential risks and benefits for their use in future disease modelling and therapeutic applications.