Skip to main content

Pyruvate kinase, a metabolic sensor powering glycolysis, drives the metabolic control of DNA replication



In all living organisms, DNA replication is exquisitely regulated in a wide range of growth conditions to achieve timely and accurate genome duplication prior to cell division. Failures in this regulation cause DNA damage with potentially disastrous consequences for cell viability and human health, including cancer. To cope with these threats, cells tightly control replication initiation using well-known mechanisms. They also couple DNA synthesis to nutrient richness and growth rate through a poorly understood process thought to involve central carbon metabolism. One such process may involve the cross-species conserved pyruvate kinase (PykA) which catalyzes the last reaction of glycolysis. Here we have investigated the role of PykA in regulating DNA replication in the model system Bacillus subtilis.


On analysing mutants of the catalytic (Cat) and C-terminal (PEPut) domains of B. subtilis PykA we found replication phenotypes in conditions where PykA is dispensable for growth. These phenotypes are independent from the effect of mutations on PykA catalytic activity and are not associated with significant changes in the metabolome. PEPut operates as a nutrient-dependent inhibitor of initiation while Cat acts as a stimulator of replication fork speed. Disruption of either PEPut or Cat replication function dramatically impacted the cell cycle and replication timing even in cells fully proficient in known replication control functions. In vitro, PykA modulates activities of enzymes essential for replication initiation and elongation via functional interactions. Additional experiments showed that PEPut regulates PykA activity and that Cat and PEPut determinants important for PykA catalytic activity regulation are also important for PykA-driven replication functions.


We infer from our findings that PykA typifies a new family of cross-species replication control regulators that drive the metabolic control of replication through a mechanism involving regulatory determinants of PykA catalytic activity. As disruption of PykA replication functions causes dramatic replication defects, we suggest that dysfunctions in this new family of universal replication regulators may pave the path to genetic instability and carcinogenesis.


Propagation of genomes between generations requires timely and accurate duplication of DNA by well-known (hereafter termed “classical”) regulatory mechanisms that control the activity of replication initiators and origins [1, 2]. However, for coordinating DNA synthesis with nutrient richness and growth rate, replication must also be under a metabolic control. This additional level of control operates by modulating the initiation and elongation phases of replication. In bacteria, the net result is a precise and reproducible timing of DNA synthesis in the cell cycle across a wide range of nutritional conditions and growth rates [3,4,5,6]. In eukaryotes, the metabolic control is thought to confine DNA synthesis to the reduction phase of a redox metabolic cycle, reiterated several times per cell cycle [7,8,9,10,11,12,13]. The mechanism of the metabolic control of replication remains largely unknown as well as its interface with “classical” replication control functions and importance for cell survival.

In bacteria, long-standing hypotheses postulating that the metabolic control of replication operates by modulating the concentration of replication precursors or the concentration of the active form of the replication initiator DnaA (DnaA-ATP), have been challenged [14,15,16,17,18]. It is also unlikely that this control operates by modulating the activity of DnaA and oriC regulators, as replication still responds to metabolism in cells lacking such regulators [15, 19, 20]. Hence, several groups have argued that the metabolic control of replication is a multifactorial process which senses nutrient richness and communicates this information to the replication machinery (see for instance [7, 9, 13, 21,22,23,24]. The guanosine tetra- and penta-phosphate that signals the cellular metabolic state and impacts replication initiation and elongation at basal concentrations may be one of these regulators [25,26,27] (see however [15, 28]).

Nutrient richness is exquisitely sensed by a metabolic area termed central carbon metabolism (CCM) [29, 30]. This metabolic area encompasses about 30 highly conserved enzymes and reactions that extract from nutrients the precursors and energy needed for macromolecular synthesis and biomass production. By directly sensing the supply and the demand in biosynthetic reactions, CCM may be at a strategic position for orchestrating the metabolic control of replication. This hypothesis is supported by an ever-increasing core of data uncovering CCM-replication links from bacteria to higher eukaryotes.

In Escherichia coli, several CCM-replication links connect the replication initiator DnaA to CCM metabolites. They may thus couple initiation to CCM activity. First, cyclic AMP, a metabolite that regulates CCM, interacts with DnaA to stimulate its binding to the replication origin and facilitate DnaA rejuvenation from the inactive DnaA-ADP to the active DnaA-ATP form [31]. Second, changes in the pyruvate-acetate node suppress initiation defects of a DnaA mutant (dnaA46) [32,33,34]. Third, two metabolites of the pyruvate-acetate node (acetyl-CoA and acetyl-phosphate) drive DnaA acetylation to prevent DnaA from binding ATP and oriC [35]. Other E. coli studies uncovered links that may couple elongation to metabolism [32, 33, 36]. In another proteobacteria (Caulobacter crescentus), the citrate synthase CitA which ensures the first reaction of the Krebs (tricarboxylic acid) cycle promotes initiation by negatively regulating the master driver of the cell cycle, CtrA [37].

CCM-replication links have also been discovered in the model Gram-positive bacterium Bacillus subtilis. First, subunits of the pyruvate dehydrogenase (PdhC) and related enzymes bind the origin region, the replication enzymes DnaC (helicase) and DnaG (primase) and/or inhibit replication initiation [38,39,40,41]. Second, CCM mutations suppress initiation and/or elongation defects in mutants of DnaC, DnaG, and the lagging-strand polymerase DnaE. They also disturb replication initiation and elongation in a medium-dependent manner [42,43,44]. Third, the metabolic control of replication is disturbed in cells mutated for PdhB, the glyceraldehyde 3-phosphate dehydrogenase (GapA), or the pyruvate kinase (PykA) [15]. Collectively, these results suggest that replication coupling to CCM activity in B. subtilis may depend on metabolic signals originating from the CCM area converting glyceraldehyde 3-phosphate to acetate (thick bars Fig. 1A) and on CCM-driven modulation of activities of replication enzymes involved in initiation and elongation (DnaC, DnaG, and DnaE) [43, 45,46,47,48,49]. Although the underlying mechanism remains elusive, it was proposed that in the rich LB medium, CCM modulates initiation by altering the functional recruitment of DnaC, DnaG, and DnaE to oriC and affects elongation by altering the activity of the lagging-strand polymerase DnaE in replication forks [44].

Fig. 1
figure 1

Replication analysis in wild-type (WT) and pykA null cells growing exponentially in MC medium. A Schematic representation of CCM. Glyco: Glycolysis; Gluco: gluconeogenesis; PPP: pentose phosphate pathway; TCA: tricarboxylic acid cycle; O: overflow pathway; G3P: glyceraldehyde 3-phosphate; PEP: phosphoenolpyruvate; Gray arrows: carbon flux; Thick bars: CCM area genetically linked to replication genes [42, 44]. B Growth rate and pyruvate kinase activity. Cells were grown exponentially for more than 20 generations using successive dilutions. Growth was assessed by spectrophotometry (OD650nm, a typical experiment is shown) and pyruvate kinase activity was determined in crude extracts (mean from six independent experiments). C Ori/ter ratio: Ori/ter ratios were determined by qPCR using as template total DNA extracted from growing cells. Numbers in brackets stand for the number of independent measurements. The Mann-Whitney U test showed that the ratios in wild-type and ΔpykA cells are significantly different at p < .05 (Table S1). D DNA elongation: Parameters of DNA elongation were determined using a marker frequency analysis by qPCR. The nearly monotonous decrease of marker frequencies from the origin to the terminus showed that there is no pause site along the chromosome (a typical experiment is shown). Numbers in brackets refer to C period (mean and SD from at least 3 experiments) and mean fork speed. E Number of origins per cell: To determine the number of origins/cell, chloramphenicol was added to exponentially growing cells. The drug inhibits replication initiation and cell division but allows completion of ongoing rounds of replication. After 4 h of drug treatment (runout experiment), cells were analyzed by flow cytometry after DNA staining. Panels: typical runout DNA histograms with the % of cells containing 4 and 8 chromosomes. Numbers in brackets stand for the number of independent reiterations of the experiment. The number of origins/cell (mean and SD) is given below the strain name. See “Methods” for details

CCM-replication links were also found in eukaryotes ( [50,51,52,53,54] and see below). At first, these findings are surprising as metabolism and replication occur in different cellular compartments (cytoplasm and nucleus, respectively). However, data accumulating since the late 1950s show that this dogma has exceptions, with numerous CCM enzymes shown to be present in the nucleus for ensuring non-metabolic (moonlighting) functions (see for instance [55,56,57,58,59,60,61,62]. Hence, CCM-replication links in eukaryotes may involve nuclear CCM determinants. Several data support this notion. For instance, the timing of origin firing in eukaryotic cells depends on an increase in acetyl-CoA that promotes histone acetylation in the nucleus. This increase is geared by a metabolic cycle in yeast and by nuclear forms of the ATP-citrate lyase and Pdh complexes in mammalian cells [63,64,65]. In addition, the nuclear form of GAPDH and the lactate dehydrogenase are an integral part of a cofactor which, in association with the transcription factor Oct-1, stimulates S phase progression by inducing expression of histone H2B [66, 67]. In contrast, repression of histone genes by exogenous pyruvate delays S phase entry [68]. It was also found that a nuclear form of the phosphoglycerate kinase (PGK) interacts with the protein kinase CDC7 to positively regulate replication initiation by keeping the stimulatory effect of CDC7 on the MCM helicase [69]. Finally, some nuclear CCM enzymes (PGK, GAPDH and lactate dehydrogenase) modulate the activity of eukaryotic replicative polymerases (Polα, Polε, and Polδ) in vitro [70,71,72]. Collectively, these data suggest that CCM may be of prime importance for the metabolic control of replication from bacteria to human cells.

Here, we report on the relationship between PykA and DNA replication in B. subtilis. PykA is an ancestral protein that is highly conserved and abundant throughout the evolution tree. It catalyzes the irreversible transfer of phosphate from phosphoenolpyruvate (PEP) to ADP resulting in the production of pyruvic acid and ATP in the last reaction of glycolysis. The 3D structure of the functional homotetramer has been solved for many organisms [73]. It highlights the conservation of their global architecture and of their active site and enabled the identification of residues important for substrates (PEP and ADP) binding and for phosphoryl transfer in B. subtilis PykA (Fig. S1A). Almost all pyruvate kinases are homotropically activated by the substrate PEP. They are also allosterically regulated by heterotropic effectors. Crystal structures of the four PykA states (apo, substrate-complex, effector-complex, and substrate-effector-complex) show that the mechanism of PykA regulation involves conformational changes within each subunit and between neighboring subunits promoted by the binding of effectors. This suggests that cells harbor different conformers at levels depending on effector concentrations and typifying nutritional conditions, with the active R-state conformer at a concentration sufficient to fulfill cellular biosynthetic needs. The heterotropic effectors of B. subtilis PykA are AMP and ribose 5-phosphate [74]. They are predicted to stabilize the active conformer in a process involving an extra C-terminal domain (see below) [75]. Interestingly, B. subtilis PykA is constitutively expressed and produced in rather large amounts across a wide range of nutritional conditions, including gluconeogenic media where pykA is dispensable for growth [76, 77]. This suggests that PykA may have more than one function.

The PykA proteins of B. subtilis and related species have an extra C-terminal sequence (residues 477-585) [78, 79]. It interacts with the catalytic domain (termed thereafter Cat) through a hydrogen bond between E209 and L536 of neighboring subunits [74]. Interestingly, the sequence and 3D structure of the extra C-terminus are homologous to one of the PEP utilizer (PEPut) domains found in the EI component of the phosphoenolpyruvate:carbohydrate phosphotransferase system, the pyruvate phosphate dikinase (PPDK), and the PEP synthase (PEPS) (Fig. S1B) [74, 78]. Interestingly, a TSH motif mapping immediately downstream of the Cat-PEPut interacting L536 residue is conserved in PEPuts (Fig. S1B) [78]. In EI, the histidine residue of the motif is phosphorylated at the expense of PEP and the phosphoryl group is transferred to a protein (Hpr) during the process of sugar transport [80, 81]. In PPDK and PEPS, this H residue is a key part of the catalytic site: it is essential for the reversible conversion of pyruvate into PEP and the reaction depends on the transfer of a phosphoryl group from ATP or PEP to pyruvate or AMP through transient phosphorylation of the conserved histidine [82, 83]. The T residue of the TSH motif is phosphorylated for inhibiting the catalytic activity of H in PEPS and PPDK and its phosphorylation/dephosphorylation is catalyzed by a serine/threonine kinases of the DUF299 family using ADP as donor [84,85,86,87]. The function of the S residue of the TSH motif is unknown [84]. This residue is however phosphorylated in B. subtilis and Arabidopsis thaliana [88,89,90]. As the Cat-PEPut interaction stabilized by a hydrogen bond between E209 and L536 is predicted to assist heterotrophic effectors in regulating PykA activity [75] and as L536 is in immediate proximity with the TSH motif, posttranslational modifications of the motif may contribute to PykA regulation by modulating the Cat-PEPut interaction. This hypothesis is consistent with studies showing that posttranslational modifications, which are prevalent from bacteria to eukaryotes, affect CCM enzyme activities through different mechanisms and help the cell to respond to nutrient availability and fluctuations [91, 92].

Here, we demonstrate that PykA is a medium-dependent effector of replication with PEPut operating as an inhibitor of initiation and Cat, as a stimulator of elongation. These functions depend neither on PykA catalytic activity nor on changes in the metabolome measurable by LC/MS. In contrast, they may depend on functional interactions of PykA with the helicase DnaC, primase DnaG, and lagging-strand polymerase DnaE. Surprisingly, the PEPut and Cat replication functions are of prime importance for replication timing even in the presence of “classical” replication control functions. Finally, we found that PEPut regulates PykA catalytic activity and that PykA allosteric regulation and PykA replication functions share requirements in substrate binding to Cat, the Cat-PEPut interaction, and the TSH motif and its phosphorylation status. We propose from these findings that cells use the metabolic sensor capacity of PykA for integrating metabolism with DNA replication and that PykA typifies a new family of regulators gearing the metabolic control of replication. An intertwining regulatory model is discussed.


Replication defects in pykA null cells occur in a medium where PykA is dispensable for growth

Replication phenotypes in ΔpykA cells were found in glycolytic regimens where PykA depletion causes a dramatic reduction in cellular metabolism and growth rate [15, 42, 44]. Here, we showed that a mere decrease in metabolism and/or growth rate is not mandatory for replication phenotypes to occur as replication defects in ΔpykA cells still happen in a gluconeogenic medium where PykA is dispensable for growth (MC medium, Fig. 1B). In this medium, runout flow cytometry studies showed a lower number of origin sequences per cell in the mutant in comparison to the wild-type strain, while qPCR studies showed that the mutant has a higher ori/ter ratio (the relative copy number of origin sequences versus terminus sequences), an extended C period (the time required for replicating the entire chromosome), and a lower speed of replication forks (Fig. 1C–E and Table S1). It is inferred from this that a moderate metabolic change ensuing depletion of PykA catalytic activity (this activity is reduced 30-fold in the mutant compared to the wild-type strain (Fig. 1B)) and/or the absence of the PykA protein, suffices for replication defects to occur in ΔpykA cells.

Cat and PEPut affect replication independently of their metabolic functions

Ori/ter ratio analysis in Cat mutants

To further investigate the relationship between PykA and replication, DNA synthesis was analyzed in PykA catalytic mutants grown in MC. The mutants used were either deleted for the catalytic domain (pykAΔcat) or contained point mutations impeding either PEP binding, ADP binding (pykAR32A, pykAR73A, pykAGD245-6AA, and pykAT278A), or the phosphoryl transfer during catalysis (pykAK220A) (Fig. S1A) [73, 93]. We also used a strain containing a 27-amino-acid deletion (residues 208-234; pykAJP) encompassing residues involved in the stabilization of the Cat-PEPut interaction (E209) and in the phosphoryl transfer (K220). This deletion is a spontaneous suppressor of thermosensitive DnaE mutations [42]. Replication was assessed using the ori/ter ratio, a parameter sensitive to initiation and elongation defects and convenient for screening a rather large number of strains.

Four classes of mutants characterized by a low, wild-type, high and notably high ratio were identified using the Mann-Whitney statistical U test run at a significance p < 0.01 (Fig. 2A, Table S1). Control experiments showed that the tested mutations dramatically inhibit (or abolish) PykA metabolic activity (Fig. 2A) and do not affect growth in MC (Fig. S2, left panel). Hence, replication defects in Cat mutants occur in the absence of a mere decrease in growth. Residues important for a proper ori/ter ratio are R32, R73, T278, GD245/6, and those deleted in PykAJP (208-234). Since mutants display different ori/ter ratios and similar residual pyruvate kinase activities, no obvious relationship links PykA catalytic activity to replication.

Fig. 2
figure 2

Ori/ter ratios and PykA catalytic activities in Cat and PEPut mutants. Ori/ter ratios were grouped (color code) according to the Mann-Whitney U test at a significance of p < 0.01 (see Table S1 for details). Numbers in brackets refer to the number of independent iterations. The horizontal gray bars highlight the wild-type ratio area. Pyruvate kinase activities (bolded numbers expressed in % of wild-type activity) in crude extracts were determined from at least three independent experiments (SD/means < 10%). A Cat mutant analysis. B PEPut mutant analysis. See Fig. 1 for details

Ori/ter ratio analysis in PEPut mutants

Next, we analyzed the ori/ter ratio in a collection of PEPut mutants comprising a strain encoding a PykA protein deleted for PEPut (pykAΔPEPut) and mutants of the TSH motif in which individual amino acids, or the whole motif, were replaced by A or D, the latter residue mimicking phosphorylation. We also tested mutants of the Cat-PEPut interaction (pykAE209A and pykAL536A).

Three classes of PEPut mutants characterized by a very low, low, and wild-type ratio were identified by the Mann-Whitney U test run at a significance p < 0.01 (Fig. 2B, Table S1). Residues important for a proper ratio are the S residue of the TSH motif and the amino acids stabilizing the Cat-PEPut interaction (E209 and L536), as corresponding A mutants have a low ratio. The T and H residues of the TSH motif may also be important for replication in conditions favoring their phosphorylation, as the corresponding D mutants are associated with a very low and low ratio, respectively (Fig. 2B). Interestingly, PEPut mutations notably affect (positively and negatively) PykA catalytic activity (Fig. 2B) and do not alter growth in MC (Fig. S2, right panel). Moreover, PykA activity and ori/ter ratio do not covary (Fig. 2B).

We inferred from these findings that Cat and PEPut are endowed with non-metabolic functions allowing PykA to modulate DNA replication independently of its catalytic activity, and that the engagement of PEPut in replication may depend on posttranslational modifications of the TSH motif. Our results also show that PEPut is a new, potent regulator of PykA catalytic activity.

Cat and PEPut mutations do not affect the metabolome in MC

To better characterize the effect of PykA mutations on cellular metabolism in MC, the metabolome of cells hosting a PykA catalytic activity ranging from 3 to 160% (the wild-type strain and mutants pykAGD245/6AA, pykAT278A, pykAT>D, pykAL536A, and pykAT>A which have a PykA activity 3, 4, 10, 53, and 160% of wild-type cells, respectively) was analyzed. For each strain, metabolomes were prepared from 4 independent cultures and analyzed by LC/MS. Untargeted analysis was carried out using XCMS and yielded 1960 and 14335 features in the positive and negative ionization modes, respectively, with automatic compound annotation (Kegg database; Table S2). Principal components analysis (PCA) was conducted on the full data matrix. The PCA score plots (Fig. 3A, B) showed that Dimension 1 capture most of the variance in both modes. In the positive mode, this axis clearly separates the wild-type and pykAT>D samples while the remaining data are mixed into each other. In the negative mode, pykAGD245/6AA, pykAT278A, pykAL536A, and pykAT>A are better segregated along Dimension 2, while wild-type and pykAT>D remain well separated in Dimension 1, with pykAT>D samples clearly away from the others. A closer inspection of PCA loadings plots identified 14 variables having the largest effect on each component. Of these, six have a biased signal due to poor correction of retention times and three are not significantly affected (fold change <3) (Stuani et al. 2014) as compared to the wild-type strain. Remaining variables pointed to molecules that do not correspond to any compound of the Kegg database or to molecules that are not expected to participate in B. subtilis metabolism. Altogether, this analysis did not detect significant changes in metabolite (including detectable nucleotides) concentrations and uncovered no relationship between score plots and PykA activity.

Fig. 3
figure 3

Metabolome analysis of wild-type and pykA mutants in MC. A, B PCA scores plots (first two components) in the positive (A) and negative (B) ionization mode. Ellipses are the 95% confidence regions. Indicated are the relative activities of PykA. Data correspond to 4 independent extractions (liquid cultures). C Comparison of PykA substrate and product contents in pykAT>D and wild-type (WT) cells. D Comparison of the contents of compounds annotated as pseudaminic acid (Pse) and CMP-pseudaminic acid (CMP-Pse) in pykAT>D and wild-type cells. Data in C and D correspond to 3 independent extractions (solid cultures).*, p > 0.05 ; **, p < 0.01 (Welch’s T-test). Values in bold panels C and D indicate the fold change for each metabolite (wild-type versus pykAT>D)

To further characterized the effect of PykA mutation on the metabolome, a study based on a procedure allowing detection of metabolites with short lifetime was carried out on the wild-type and pykAT>D strains (see Experimental Procedures) [94]. This analysis confirmed the inference of Fig. 3A,B and the lack of significant changes in NTP pools (Fig. S3). A targeted analysis on PykA substrates and products revealed no significant change in ATP, pyruvate, and ADP concentrations and a 4.9 lower PEP content in the mutant (Fig. 3C). The data also uncovered an ion that accumulated in the pykAT>D mutant, annotated as pseudaminic acid (Pse) in the Kegg database (Fig. 3D). Consistently, another accumulating ion annotated as a derivative of Pse (CMP-Pse) was found in the mutant. Putative Pse was observed mainly in the negative mode at m/z 333.1296 but was also detected in the positive mode at m/z 335.1443 (Fig. S4A, B). These two m/z values were anticipated to be the deprotonated [M−H] and protonated [M+H]+ forms of the molecule, respectively. Fragmentation pattern and hydrogen/deuterium exchanges (HDX) were also fully consistent with the Pse structure. Pse harbors 6 mobile protons. We determined the mobile proton number of m/z 333 in positive and negative ionization modes using in-solution HDX approach. Mobile protons (7 and 5) were observed in the positive and negative mode, respectively; i.e., 6 mobile protons for the neutral form (Fig. S4C-F). However, data were also in agreement with a previous study which deciphered the biosynthetic pathway of a Pse isomer: legionaminic acid (Leg) [95]. As this isomer should display the same fragmentation patterns under collision-induced dissociation (CID) as Pse, we considered m/z 333 in mutants as Leg. The signal was indeed absent in cells deleted for the Leg synthase (ΔspsE), confirming m/z 333 as Leg. Moreover, Leg accumulated in cells inactivated for the enzyme converting Leg into CMP-Leg (CMP-Leg synthase; ΔspsF) and the CMP-Leg signal was detected in wild-type and pykAT>D cells but not in mutants of the Leg pathway (Fig. S4G, H). Finally, PEP depletion in pykAT>D likely result from Leg production as PEP is a precursor of Leg [95] and as metabolome data suggested that PEP concentration inversely correlates with Leg concentration (Fig. 3C, D).

Changes in Leg, CMP-Leg and PEP concentration prompted us to analyze replication in Leg pathway mutants (ΔspsE, ΔspsF, and ΔspsABCDEF) in wild-type and pykAT>D contexts. Results showed that neither the depletion of Leg, CMP-Leg, and/or PEP nor the accumulation of Leg affect the ori/ter ratio in both contexts showing that the Leg pathway is not involved in replication (ratios in single and double mutants vary by less than 5% compared to the wild-type and pykAT>D parental strains). Collectively, our findings show that PykA mutations do not notably disturb the cellular metabolome (including nucleotide pools) in MC. Some mutants may however significantly (>3-fold) impact PEP, Leg, and CMP-Leg concentration but these changes are not responsible for the replication phenotypes.

PEPut and Cat affect initiation and elongation, respectively

Changes in ori/ter ratio may result from changes in initiation and/or elongation. To determine which replication step is primarily defective in PykA mutants, cells replicating the chromosome from a plasmid replicon rather than from the cellular DnaA/oriC initiation complex were used. As argued previously [15, 44], if replication defects in pykA mutants were to result from changes in initiation, cells replicating their genome from the plasmid replicon would not suffer from metabolic mutation and would thus exhibit the ori/ter ratio of PykA+ oriN-dependent cells. In contrast, if replication defects were to result from changes in elongation, plasmid replicon-dependent cells would still suffer from the pykA mutation and would thus have an ori/ter ratio different from PykA+ oriN-dependent cells. The analysis was carried out in the ΔoriC context rather than in the ΔdnaA context to avoid any interference of DnaA depletion on genome expression [96].

Results showed a ratio typical to PykA+ oriN-dependent cells in PEPut and Cat-PEPut interaction mutants and a clearly reduced (pykAJP, pykAGD256/6AA) or increased (pykAT278A) ratio in Cat mutants (Fig. 4). Hence, PEPut and the Cat-PEPut interaction mutants are primarily affected in initiation while Cat mutants are primarily affected in elongation. This suggests that the commitment of PykA in replication involves two separate non-metabolic functions, one geared by PEPut and its interaction with Cat which affects initiation, and the other geared by Cat which affects elongation.

Fig. 4
figure 4

Primarily replication defect in pykA mutants. Ori/ter ratios were measured in cells replicating the chromosome from a plasmid replicon instead of the natural chromosomal replication origin. See Fig. 1 for details

PEPut is an inhibitor of initiation and Cat, a stimulator of elongation

Microscopy studies showed that nucleoids are efficiently segregated in pykA mutants grown in MC (Fig. S5, top raw). Therefore, a cell cycle analysis was carried out in wild-type cells and pykA mutants using flow cytometry (runout experiments) and qPCR (marker frequency analysis), as described previously [97]. Representative results are given Fig. S5 (middle and bottom rows), and the complete set of data is summarized in Table S3. Results showed that PykA mutations can affect replication and cell cycle to remarkably high degrees and confirmed the prime importance of PEPut and Cat in initiation and elongation, respectively. More specifically and by comparison to the wild-type strain, cells lacking Cat (pykAΔCat) and most (5/7) of the remaining Cat mutants replicate DNA remarkably slowly (600–740 bp/s instead of 830 bp/s in the wild-type strain) and have a marginal initiation defect (Table S3). In contrast, cells lacking PEPut (pykAΔPEPut) and 3 other PEPut mutants (pykAT>D, pykATSH>DDD, and pykATSH>AAA) suffer from a dramatic initiation phenotype: they initiate replication in mother cells (i.e., in cells containing two origins) rather than in grandmother cells (i.e., in cells containing four origins) and their cell cycle lasts 2 instead of 3 generations (Table S3). Additionally, two of these mutants (pykAT>D and pykATSH>DDD) replicate DNA extremely fast (about 1000 bp/s). The dramatic effects the pykAT>D (initiation) and pykAGD245/6AA (elongation) mutations have on the cell cycle and replication pattern are highlighted Fig. 5.

Fig. 5
figure 5

Cell cycle and replication pattern of wild-type, pykAGD245/6AA and pykAT>D cells. Vertical red lines stand for three successive generations: grandmother: 0τ to 1τ; mother: 1τ to 2τ; daughter: 2τ to 3τ; 28.4 min each). The C and D periods of cell cycles ending at time 1τ, 2τ, and 3τ are indicated as lines and boxes colored in black, blue and yellow, respectively. The cell cycle in green will end at time 4τ (not shown). Gray areas refer to replication periods spanning every cell cycle and numbers correspond to the proportion of cells in each period. Replication patterns typifying each period are given above the gray areas using the cell cycle color code. Ai: age of replication initiation; At: age of replication termination. Activity stand for pyruvate kinase activity

Changes in nucleotide pools affect replication. In pykA mutants, such changes may result from an indirect effect of PykA mutations on CCM reactions producing nucleotide precursors or from a direct effect on ATP concentration (PykA produces ATP during the PEP to pyruvate conversion) and/or (d)NTPs production (the PykA catalytic site is promiscuous, using PEP as phosphoryl donor to phosphorylate (d)NDPs with various efficiencies [98]). If so, cells with similar defects in PykA activity would have similar pools and hence, similar replication phenotypes. Results do not support this assumption. Initiation is dramatically stimulated in four PEPut mutants having a low to moderate PykA activity (10–40%), while it is nearly unchanged in 13 mutants having a PykA activity ranging from 3 to 160% (Table S3). Moreover, very low, wild-type, and very fast fork speeds were found in strains with similar residual PykA activities (3–10%) (Table S3). Additionally, the metabolome analysis did not reveal significant changes in nucleotide pools (Fig. 3 and Table S2).

Collectively, our findings showed that PykA has two distinct replication functions with Cat acting as a stimulator of replication fork speed and PEPut, as an inhibitor of initiation. As the inhibitory activity of PEPut is massively abrogated in pykATSH>AAA, pykATSH>DDD, and pykAT>D cells, our results strengthen the notion that the TSH motif and its phosphorylation at T are of prime importance for initiation. Finally, given that replication phenotypes were found in cells fully proficient in “classical” replication control functions, this work showed that PykA is a new master regulator of replication.

PEPut affects replication initiation and PykA catalytic activity in a medium-dependent manner

Results presented above showed that PykA is endowed with two replication functions in MC, a rich gluconeogenic medium. To investigate whether these functions are expressed in other nutritional conditions, replication parameters (ori/ter ratio, C period and number of ori/cell) of wild-type, pykAT>A, and pykAT>D cells were assessed in another gluconeogenic broth (minimal medium supplemented with malate, M) and in two glycolytic media (minimal medium supplemented with glucose (G) or glucose + casa (GC)). Results showed wild-type parameters under the three tested conditions except in pykAT>D cells grown in M (Fig. 6A, B). In this condition, a low ori/ter ratio, a low number of ori/cell, and an enlarged C period were found. However, while the three strains have similar growth rates in GC and G, the pykAT>D strain grows poorly in M compared to wild-type and pykAT>A cells (generation times 140 versus 70 min) (Fig. 6B and S5, left panel). This phenotype is a direct consequence of the pykAT>D mutation, as catalytic mutants (pykAGD245/6AA and pykAK220A) and pykA knockout cells grow as fast as the wild-type and pykAT>A strains in M (Fig. S6, left panel). In order to further characterize the effect of the pykAT>D mutation on replication in M, replication parameters of wild-type cells growing in proline at the similar rate as pykAT>D in M were determined. Results (Fig. S6, right panel) showed that the mutant has a lower ori/ter ratio (1.3 versus 1.5), an enlarged C period (77 versus 58 min) and a lower number of ori/cell (compare runout DNA histograms) than wild-type cells in proline. In addition, they showed that the mutant likely suffers from an early initiation phenotype, as runout DNA histograms uncovered a significant proportion of non-initiated cells with 1 chromosome in pykAT>D grown in M while the vast majority of non-initiated cells in proline cultures of the wild-type strain has 2 chromosomes. We concluded from these findings that PEPut inhibits initiation in gluconeogenic media (MC and M) through a mechanism abrogated by the pykAT>D mutation. The lack of phenotype in pykAT>D cells grown in GC and G either suggests that PykA does not impact initiation in glycolytic conditions or does so through a mechanism insensitive to the pykAT>D mutation. Hence, the replication function of PEPut is medium dependent.

Fig. 6
figure 6

Replication and PykA catalytic activity analysis in wild-type, pykAT>A, and pykAT>D cells grown in different media. A Ori/ter ratio and C period. B Runout DNA histograms. Numbers in brackets stand for generation times (min). C PykA catalytic activity. Numbers refer to ratios as follows: top row: pykAT>A/wild-type; bottom row: wild-type/pykAT>D. Data in A and C are from 3–6 biological repeats (mean and SD). Representative runout DNA histograms are shown in (B)

Next, we analyzed the role of PEPut in PykA catalytic activity in the strains and growth conditions used above. Results showed that PykA activity is increased in pykAT>A and reduced in pykAT>D contexts by comparison to wild-type cells in M, GC, and G (Fig. 6C) as it is in MC (Fig. 2) suggesting that T phosphorylation negatively regulates PykA catalytic activity in a large range of nutritional conditions. Interestingly, the level of T phosphorylation may be higher in glycolytic than in gluconeogenic media as the activity produced by wild-type cells is closer to that of pykAT>D cells in glycolytic media and closer to that of pykAT>A cells in gluconeogenic growth conditions (compare the ratios pykAT>A versus wild-type and wild-type versus pykAT>D Fig. 6C). It is inferred from this that the PEPut-driven regulation of PykA catalytic activity is sensitive to the polarity of the carbon flux travelling CCM.

Overall, our findings suggested that the level of T phosphorylation in the TSH motif varies with growth conditions and has two distinct functions: it negatively regulates PykA catalytic activity and alleviates the inhibitory activity of PEPut on replication initiation. However, while the former activity occurs at mid-levels of T phosphorylation, the latter activity may operate at high T phosphorylation levels.

The purified PykA protein modulates replication activities of DnaE and DnaC

The effect of PykA on replication may involve the recruitment of the glycolytic enzyme to the initiation and/or elongation replication machineries. It was previously shown that PykA is produced at a rather high concentration [77] and is diffused throughout the cytoplasm [99]. Hence, for localizing PykA, we constructed strains encoding PykA-mCherryBSU or mCherryBSU-PEPut fusions (wild-type or mutated) from an inducible promoter (Physpank) at hundreds to thousands proteins per cell [46]. The constructions were inserted at an ectopic position (the amyE locus) in a context lacking or containing the natural pykA gene. Microscopy studies of living cells showed a diffused mCherryBSU signal at all inducer concentrations (Fig. S7, top panels). However, upon cell fixation and at low inducer concentration (0-5 μM), irregular spotty signals were observed (Fig. S7, bottom panels). Unfortunately, we were unable to improve the fixative procedure and the quality of the images, preventing further analyses on PykA localization.

In order to investigate whether PykA affects replication by directly modulating the activity of replication enzymes, an in vitro approach was used. The PykA protein was heterologously expressed and purified, and control experiments showed that the purified PykA was metabolically active and formed the expected stable functional tetramer (Fig. S8). The effect of PykA was then tested on DnaC, DnaG, and DnaE, as these replication enzymes are genetically linked to PykA [42]. The DnaC helicase melts the duplex DNA at oriC during initiation and separates the DNA strands in replication forks during elongation. The DnaG primase synthesizes RNA primers at the origin and in replication forks. These primers are extended by DnaE to start leading-strand synthesis (which is mainly carried out by PolC) and to ensure partial or complete synthesis of the lagging strand. Previous studies showed that DnaC, DnaG, and DnaE form a ternary complex ensuring important roles during replication initiation and elongation [43, 45, 47, 49]. The three replication enzymes were purified and their activities were tested using previously described assays [43, 49] in the presence or absence of PykA. It is of interest to note here that these assays do not allow PykA activity (they lack the PykA substrates PEP and ADP) and probe key features of DNA elongation: chain polymerization and double-strand separation at forks. In contrast, they do not probe initiation as they give no information on origin melting and protein recruitment at melted origins.

Replication assays with DnaE were carried out at a polymerase concentration (10 nM) that produces a low amount of replication products in order to facilitate the detection of stimulatory effects. In reaction mixtures containing M13 ssDNA annealed to a 20-mer DNA oligonucleotide and DnaE in combination or not with equimolar concentrations of PykA, a substantial stimulation of DnaE polymerase activity by PykA was found in terms of both sizes and amounts of nascent DNA synthesized (Fig. 7A, left panel). A similar stimulation was observed with a 60-mer oligonucleotide annealed onto M13 ssDNA and a 15-mer oligonucleotide annealed onto a 110-mer oligonucleotide (Fig. S9A and data not shown). The stimulation was specific to PykA as it was not observed with equimolar amounts of BSA and M13 ssDNA primed with a 20-mer (Fig. 7A, right panel) or a 60-mer (Fig. S9B). The lack of DnaE stimulation by BSA was further confirmed at a 50-fold excess concentration over DnaE with the 20-mer primed M13 ssDNA (Fig. S9C) (note that the marginal stimulation observed at very high (500-fold) BSA excess is artifactual, acting likely as a blocking agent preventing adhesion of DnaE to the plastic reaction tubes). Titration experiments and gel shift assays showed that the stimulation of DNA synthesis by PykA was not due to a stimulation of DnaE binding to a primed template (Fig. S10).

Fig. 7
figure 7

PykA stimulates the DNA polymerase activity of DnaE. A Representative alkaline agarose gels showing DnaE primer extension time courses (30, 60, 90, 120, and 150 s) with DnaE (10 nM) alone and in the presence or absence of PykA (40 nM tetramer) or BSA (40 nM), as described in “Methods”. Lanes M and C represent DNA size markers and the control radioactive substrate in the absence of any proteins, respectively. B Representative alkaline agarose gels (from three independent experiments) showing DnaE primer extension time courses (30, 60, 90, 120, and 150 s) with DnaE (10 nM) and in the presence or absence of PykA (40 nM tetramer), with molar ratios of DnaE monomer: DnaN dimer: HolB monomer:YqeN monomer:DnaX trimer:PykA tetramer set to 1:1:1:1:1:1, considering the oligomeric states of these proteins. C As in B with 2 nM DnaE and molar ratios of 1:1:1:1:1:1

Within the replisome, the polymerase activity of DnaE is stimulated by several proteins [43, 46, 49, 100], with DnaN being a potent stimulator. This ring-shaped protein is loaded at the 3’ end of primed sites by the clamp loader, encircles the DNA, and binds DnaE to form a stable complex that slides along the DNA template. In order to determine whether the polymerase activity of the DnaE-DnaN complex can be further stimulated by PykA, we carried out primer extension assays with DnaE, DnaN, proteins of the clamp loader (HolB, YqeN and DnaX), and PykA (Fig. 7B). As previously observed [43], we found that the polymerase activity of DnaE (10 nM) is strongly stimulated by DnaN and the clamp loader (compare the left panels of Fig. 7A, B). In this condition of high activity, the effect of PykA is unclear (Fig. 7B). However, this glycolytic enzyme may confer an additional stimulation to the DnaE activity as its presence leads at the last timepoints to a greater accumulation of large fragments and a clear reduction in labelled primer (Fig. 7B, compare with and without PykA). To confirm this, we carried out similar primer extension assays in the presence of lower, suboptimal concentrations of DnaE (2 nM), DnaN, clamp loader proteins, and PykA (in both sets of assay conditions, the molar ratios of proteins were kept identical). At this suboptimal DnaE concentration, no nascent DNA was detectable in the absence of PykA and significant amounts of nascent DNA fragments were synthesized by DnaE in the presence of PykA (Fig. 7C). This suggests that PykA stimulates the activity of DnaE even in conditions where its activity is strongly stimulated by DnaN. We concluded from these data that PykA stimulates the DnaE polymerase activity when the polymerase is slow (i.e., alone) and fast (i.e., in a complex with DnaN) probably via a direct interaction between PykA and DnaE. Since the purified PEPut domain does not affect DnaE activity (Fig. S11-12), the stimulation probably depends on a direct interaction between DnaE and the Cat domain or the interaction interface involves structural features of the PykA tetramer that are not preserved in the purified PEPut. The stimulation occurs with short (20-mer) and long (60-mer) primers suggesting that PykA may stimulate DnaE polymerase activity during extension of RNA primers generated by DnaG and during lagging-strand synthesis.

The helicase activity of DnaC was assayed by monitoring the displacement of a labelled 104-mer oligonucleotide annealed onto M13 ssDNA and forming a double-fork substrate with poly(dT) tails as previously [49]. To assemble a functional DnaC hexamer onto the DNA substrate, the helicase loader DnaI was added to reaction mixtures at equimolar concentrations. Results showed that DnaC is marginally inhibited in the presence of PykA (Fig. 8A). As we previously showed that DnaC activity is stimulated by DnaG [49], DnaG was added to the abovementioned reaction mixtures at equimolar concentrations. This analysis confirmed the stimulation of helicase activity by DnaG and showed that this stimulation is cancelled by PykA (Fig. 8B). Thus, in some contexts, PykA can significantly inhibit the helicase activity of DnaC. Collectively, these in vitro results suggest that PykA can modulate replication initiation and elongation through direct functional interactions with replication enzymes.

Fig. 8
figure 8

PykA directly inhibits the helicase activity of the helicase DnaC and via the primase DnaG. A Time courses (5, 10, 15, 25, and 30 min) showing the helicase activity of DnaC/DnaI in the presence or absence of PykA, as indicated. B Time courses (5, 10, 15, 25, and 30 min) showing the helicase activity of DnaC/DnaI in the presence or absence of PykA and/or DnaG, as indicated. The reactions were carried out as described in “Methods”. Representative native PAGE gels (from two independent experiments) are shown with lanes a and b representing control annealed and fully displaced (boiled) control DNA substrates. Data were plotted as a percentage of displaced primer versus time using GrapPad Prism 4 software. Error bars show the standard deviation from two independent repeat experiments


PykA is one of the most conserved and abundant proteins of the evolution tree and its ability to convert PEP and ADP into pyruvate and ATP during glycolysis is known in exquisite detail [73, 75]. Its activity is allosterically regulated through conformational changes within each subunit and between neighboring subunits promoted by the binding of effectors. Interestingly, an ever-increasing core of data showed that pyruvate kinases from bacteria to higher eukaryotes, akin other CCM enzymes, often ensure non-metabolic (moonlighting) functions in processes as diverse as transcription, viral replication, angiogenesis, pathogenicity, and tumorigenesis [60,61,62, 101, 102]. Here, we provide new insights on how the catalytic activity of B. subtilis PykA is regulated and demonstrate that this enzyme ensures non-metabolic functions important for replication initiation, elongation, and timing.

Regulation of PykA catalytic activity

The B. subtilis PykA protein comprises the canonical catalytic domain (Cat) and an extra C-terminal peptide 110-amino-acid-long. This peptide, termed PEPut, is homologous to domains embedded in other metabolic enzymes and known to ensure catalytic and regulatory functions depending on phosphorylation of a conserved TSH motif (see  the "Background" section). Crystallographic studies suggested that Cat and PEPut interact through a hydrogen bond between E209 and L536 [74], and this interaction was proposed to assist heterotrophic effectors in stabilizing the active R-state conformation [75]. Here, we demonstrate that PEPut negatively regulates PykA catalytic activity up to 16-fold and that this process depends on the Cat-PEPut interaction, the TSH motif and its phosphorylation status (Fig. 2). A key determinant of this regulation may be the phosphorylation level of T in the TSH motif since the most dramatic changes in PykA activity occur in conditions mimicking the absence (T>A context, no inhibition) or massive (T>D context, strong inhibition) T phosphorylation (Fig. 2). Interestingly, T phosphorylation may vary with the polarity of the carbon flux travelling CCM, being higher in glycolytic than in gluconeogenic growth conditions (Fig. 6C). Our findings thus enforce the stabilization model of the active R-state conformation by Cat-PEPut interaction and further suggest that this stabilization is negatively regulated by T phosphorylation in the TSH motif. Hence, and assuming that knowledge in TSH motif phosphorylation prevails in B. subtilis PykA (see above), PykA regulation may depend on the concentration of phosphoryl donors (potentially PEP and ADP) and the activity of kinases/phosphatases for TSH phosphorylation, in addition to the concentration of PykA effectors (PEP, AMP, and ribose 5-phosphate). PykA can thus be viewed as a metabolic sensor mediating multiple interactions with allosteric regulators and exhibiting different TSH phosphorylation profiles depending on signaling metabolite concentration and CCM activity. This suggests that the thousands of PykA protein present in the B. subtilis cytoplasm exist in multiple forms and proportions that may vary with cellular metabolism for providing the active R-state PykA conformer at concentrations fulfilling biosynthetic needs and for conveying information on the nutritional environment as illustrated Fig. 9A. Additional posttranslational modifications (acetylation, succinylation…) prevalent in pyruvate kinases from bacteria to eukaryotes [91, 92] may further increase the diversity of PykA forms within cells enabling a higher level of PykA regulation and a more comprehensive information on the nutritional environment. This sophisticated regulation of PykA catalytic activity may provide an elegant solution to the problem posed by the constitutive and abundant production of PykA in a wide range of nutritional conditions [76, 77].

Fig. 9
figure 9

Intertwining model of PykA catalytic activity regulation and PykA-driven metabolic control of DNA replication. A Regulation of PykA catalytic activity. Several metabolites (PEP, ATP, ADP, AMP, R5P…) are used as signaling molecules for PykA regulation. They gear Cat-metabolite interactions, Cat-PEPut interaction, and TSH motif phosphorylation in PEPut. In response to these signals, PykA adopts different conformations with different types and levels of TSH motif phosphorylation. These multiple forms of PykA typify the cellular metabolic state and drive allosteric regulation of PykA catalytic activity. This regulation impacts the metabolome and the concentration of signaling metabolites. This mechanism can be implemented by posttranslational modifications (PTMs) of PykA. B Metabolic control of DNA replication by PykA. The ability of PykA to sense signaling metabolites (effectors and phosphoryl donors) and adopt multiple forms typifying the nutritional environment is used by cells to convey a metabolic signal to the replication machinery and regulate DNA synthesis in a large range of nutritional conditions. Top panel: metabolic control of initiation. (i) In response to the concentration of a phosphoryl donor and the activity of a kinase/phosphatase system, the phosphorylation level of the T residue in the TSH motif of PEPut varies. Initiation is inhibited at low- and activated at high phosphorylation levels. High phosphorylation level may be powered by a metabolic cycle causing periodic accumulation of the phosphoryl donor at the age of initiation. Bottom panel: metabolic control of elongation. The balance between the multiple forms of PykA is sensed by replication enzymes for modulating elongation. Previous data suggest that receptors of metabolic signals conveyed by PykA are the helicase DnaC, the primase DnaG, and the lagging strand polymerase DnaE which are essential for both replication initiation and elongation (see the “Background” section). The metabolic control of initiation and elongation help cells to properly time replication in the cell cycle in a large range of nutritional conditions

PykA encodes initiation and elongation replication functions geared by PEPut and Cat, respectively

Replication phenotypes in B. subtilis cells lacking PykA were initially discovered in nutritional conditions where this glycolytic enzyme is important for growth and metabolism. Under these conditions, groups reported that pykA deletion allows (i) suppression of initiation and elongation defects due to thermosensitive mutations in the replication enzymes DnaC, DnaG, and DnaE (but not DnaI, DnaD, PolC, DnaX, and DnaN) [42, 43], (ii) stimulation of initiation [44], and (iii) alteration of the metabolic control of replication [15]. Here we show that pykA knockout cells cultivated in a medium (MC) where PykA is dispensable for metabolism and growth, also suffer from replication initiation and elongation defects (Fig. 1 and Table S3).

To gain insights into the PykA replication relationship, cell cycle parameters of Cat and PEPut mutants were monitored in MC. Results uncovered significant to strong replication phenotypes in most of the 18 tested mutants (Fig. 2 and Table S3). These phenotypes do not result from changes in PykA catalytic activity as this activity and replication parameters do not covary (Fig. 2). Additionally, they are not due to mere changes in cellular metabolism as the metabolome of B. subtilis in MC is resilient to changes in PykA catalytic activity (Fig. 3 and Table S2). In contrast, these phenotypes result from inactivation of either of two non-metabolic functions encoded by PykA: the one carried by PEPut which inhibits initiation, and the one carried by Cat which stimulates elongation. Interestingly, these moonlighting functions may operate in a medium-dependent manner (Fig. 6A,B). In vitro studies further showed that PykA functionally interacts with three enzymes essential for replication initiation and elongation. Indeed, the purified PykA protein stimulates the activity of the lagging-strand polymerase DnaE in conditions mimicking primer extension and lagging-strand synthesis, likely via physical contacts between Cat and DnaE (Fig. 7, S8-9 and S11). Moreover, the purified metabolic enzyme inhibits the helicase activity of DnaC and the stimulatory effect of DnaG on DnaC (Fig. 8).

Collectively, our report shows that PykA is endowed with a replication initiation and elongation function geared by PEPut and Cat, respectively, and that these functions may involve functional interactions with replication enzymes at oriC and replication forks. Although the strains used in this study are fully proficient in “classical” replication functions, most of the mutants lacking either PykA replication function exhibit remarkable cell cycle perturbations (Figs. 5 and 6 and Table S3). This surprising result suggests that the control geared by “classical” replication functions is far from sufficient to properly time DNA synthesis in the cell cycle and that an additional level of regulation by PykA is needed to meet this goal (at least in the gluconeogenic media used here). Therefore, we propose that PykA typifies a new family of replication control factor that operates in a medium-dependent manner along “classical” replication functions for gearing a metabolic control of DNA replication and for properly timing replication in the cell cycle.

Basic principles of the metabolic control of replication by PykA

A targeted mutagenesis showed that the inhibitory activity of PEPut on initiation depends on the Cat-PEPut interaction, the TSH motif, and its phosphorylation status (Fig. 2 and Table S3). The importance of phosphorylation is highlighted by data suggesting that high T phosphorylation abrogates the initiation inhibitory activity of PEPut in gluconeogenic (but not glycolytic) growth conditions. In these conditions, pykAT>D cells have an early initiation phenotype while initiation is properly regulated in the wild-type and pykAT>A contexts (Fig. 6A, B). Assuming that key metabolites transiently accumulate through cell cycle progression in B. subtilis as in other systems [103,104,105], one can propose a model for initiation regulation by PEPut. In this model, initiation is inhibited for most of the cell cycle by a weakly T phosphorylated form of PEPut (equivalent to the wild-type and pykAT>A contexts). However, late in the pre-initiation period, the initiation inhibitory activity of PEPut is alleviated by a massive and transient T phosphorylation (equivalent to pykAT>D context). We propose that this transient phosphorylation occurs in response to the combined action of a cell cycle-dependent accumulation of a phosphoryl donor and of a strong kinase activity of a protein of the DUF299 family [85, 86] (a sequence analysis suggested that the poorly characterized protein YqfL is the B. subtilis DUF299 like protein. Its involvement in replication is under investigation). Intriguingly, the involvement of PEPut in initiation and PykA catalytic activity have similar requirements: they both depend on the Cat-PEPut interaction, the TSH motif, and its phosphorylation status. This suggests that the mechanism of the PEPut-driven metabolic control of initiation involves determinants of the mechanism of the PEPut-driven regulation of PykA catalytic activity, as illustrated Fig. 9. In this line of thinking, we propose that the system used for massive and transient T phosphorylation late in the pre-initiation period regulates PykA catalytic activity most of the cell cycle by ensuring a nutritional-dependent basal level of T phosphorylation.

The targeted mutagenesis also showed that residues R32, GD245/6, and amino acids 208–234 are of prime importance for Cat-driven stimulation of DNA elongation (Fig. 5 and Table S3). R32 and GD245/6 are involved in ATP/ADP and PEP binding [73] while the 208–234 region resides at the heart of the catalytic site and is genetically linked to the lagging-strand polymerase DnaE [42]. In vitro data further suggested that Cat interacts with the lagging-strand DnaE polymerase stimulating its activity (Fig. 7 and S11) and studies on different systems showed that lagging-strand synthesis gears fork speed [106,107,108,109]. Collectively, these results suggest that Cat stimulates elongation through a mechanism involving Cat-metabolites and Cat-DnaE interactions. As Cat-metabolite interactions induce complex conformational changes that cells use to allosterically regulate PykA catalytic activity in a wide range of nutritional conditions [73, 75], we propose that the mechanism of the Cat-driven metabolic control of elongation involves determinants of the mechanism of the Cat-driven regulation of PykA catalytic activity as illustrated Fig. 9. Again, these functions can be implemented by PykA posttranslational modifications [91, 92]. We conclude that cells use the ability of PykA to sense the metabolic status and adopt multiple forms typifying this status to convey metabolic information to the replication machinery and drive a metabolic control to DNA synthesis. We propose that the numerous CCM-replication links described in the literature (see  the "Background" section) highlight potential determinants of the metabolic control of replication in other systems.

The findings reported here show that side chain modifications disrupting PykA replication functions can dramatically disturb replication timing. As replication timing alterations cause accumulation of chromosomal lesions (double-strand breaks) with a high mutagenic potential [110], determinants of the metabolic control of replication may form a new family of genetic stability factors. Two hallmarks of cancer cells are aerobic glycolysis and a high rate of mutagenesis. Nowadays, cancer is thought to originate from mutations in oncogenes, tumor suppressors, and/or mutator genes, and other hallmarks, including their altered metabolism, are thought to be a result of these mutations rather than their cause [111, 112]. Extrapolating the results presented here from B. subtilis to humans, we speculate that CCM mutations (or changes in metabolite pools) that disrupt the metabolic control of replication may pave the way to genetic instability and may thus be an additional root cause of cancer initiation.


Genome duplication is carried out once per cell cycle by a complex machinery tightly regulated by well-known, overlapping mechanisms. Failure in this regulation increases the risk of genetic instability and favors the emergence of diseases such as cancer. Since the 60s, genome duplication is also known to be under a metabolic control. The importance of this control for cell survival is unknown as well as its mechanism and interplay with classical replication control functions. Interestingly, an increasing core of data suggests that this mechanism may involve central carbon metabolism from bacteria to human. Accordingly, links between the enzyme ensuring the last reaction of glycolysis (the pyruvate kinase PykA) and replication in the bacterium Bacillus subtilis were found. Here, we show that the catalytic domain of B. subtilis PykA (Cat) stimulates DNA elongation while its C-terminal domain (PEPut) inhibits initiation and that these activities may involve functional interactions between Cat and enzymes involved in replication initiation and elongation. Very surprisingly, dramatic replication timing defects can occur in cells inactivated for either PykA replication activity and these activities depend on factors involved in regulation of PykA catalytic activity. We therefore propose that PykA typifies a novel family of replication regulator that acts in a nutrient-dependent manner to orchestrate the metabolic control of replication and replication timing. Since PykA’s replication functions and PykA regulation of catalytic activity share basic requirements, we propose that the mechanism allowing PykA to sense the cellular metabolic status and regulate its own activity is used by cells to control replication. We extrapolate from these results that the metabolic control of replication involves CCM enzymes endowed with a metabolic sensor activity from bacteria to human and that mutations that disrupt this control may disrupt replication timing and thus pave the way to genetic instability and diseases such as cancer.


Strains and plasmids

Strains and plasmids are listed in Supplementary Table S4. E. coli StellarTM competent cells (Clontech, Saint-Germain-en-Laye, France) were used as host for plasmid constructions. Most of the B. subtilis strains were constructed by transforming competent cells with genomic DNA (see Table S4 for details). Strains constructed using PCR products and plasmid DNA as transforming material are described below. The DGRM1173 strain containing the markerless spsABCDEF deletion was constructed as previously [95]. Genotypes of constructed plasmids and strains were checked by phenotypic analyses, endonuclease restriction, PCR analysis, and/or Sanger sequencing (Eurofins Genomics, Germany).

Construction of the pykA-tet strain (DGRM295)

In order to facilitate the construction of pykA mutants, we first inserted a tet gene immediately downstream from the transcription terminator of pykA. We first amplified with the Q5 high-fidelity DNA polymerase (New England Biolabs, Evry, Fr) the tet gene of plasmid pTB19 [113] and two chromosomal sequences flanking the site of insertion and containing the 3′ end of the pykA gene on the one hand and the ytzA plus the 3′ end of ytvI on the other hand. Then, we separated the PCR products from the parental genomic DNA by gel electrophoresis and purified the PCR fragments using the Monarch DNA Gel Extraction Kit (New England Biolabs, Evry, Fr). The purified fragments were then mixed and fused together using assembly PCR carried out with the Q5 high-fidelity DNA polymerase and the two primers flanking the fused fragment. This reaction depends on the addition of 20-bp sequences homologous to the 5′ and 3′ tails of the tet fragment to the internal side of the chromosomal PCR fragments. Competent cells of a wild-type strain cured of prophages (TF8A) were then transformed with the assembly PCR product and double crossover events integrating the tet gene downstream of pykA were selected on Tet containing plates (before plating, the tet gene was induced by incubating cells 1 h at 37 °C in the presence of 1.5 μg/mL Tet). A representative transformant, named DGRM295, was selected by PCR and DNA sequencing.

Construction of Cat and PEPut mutants

To generate Cat and PEPut mutants, pairs of PCR reactions were carried out using as template the genomic DNA of DGRM295 (pykA-tet). In each reaction, one external (i.e., in pykA or ytvI) primer and one mutagenic primer mapping in Cat or PEPut were used to generate PCR products with the desired pykA mutation at one end. PCR fragments were then assembled and the assembly products used to transformed TF8A competent cells, as described above (inactive pyruvate kinase mutants were selected on LB + Tet + Malate 0,4% w/v). Three representative transformants were selected by sequencing for nearly all constructions.

Construction of pykA mutants replicating the chromosome from oriN

To construct pykA mutants replicating the chromosome from oriN, competent cells of a TF8A derivative carrying oriN and the cat gene at the spoIIIJ locus and deleted for oriC (DGRM589, [44]) were transformed with genomic DNA of pykA mutants. Transformants were then selected on plates supplemented with appropriate antibiotic (and 0.4% w/v malate when the mutation inactivated PykA activity). Three representative transformants were generally selected. The presence of the pykA mutation was confirmed by sequencing while the presence of the oriN-cat structure and the oriC deletion were checked by measuring the size and sensitivity to EcoRI restriction of the corresponding PCR products.

Construction of strains encoding mCherryBSU fusions

To generate strains encoding mCherryBSU fusions, plasmids were first constructed (their sequence is available upon request). To construct pLJH225-Cm, a reverse PCR on pIC610 was carried out with the Q5 high-fidelity DNA polymerase to produce a linear fragment with the Physpank promoter (oriented outward) at one end and the PEPut ORF (oriented inward) at the other end. Using the same DNA polymerase and pBS1C3 as template, we then produced a PCR product containing the mCherryBSU ORF fused to a linker (GGTGGAGGTGGATCT 3n) at the C-terminus. The fragment was engineered to carry at its extremities 20-bp-long sequences homologous to the ends of the pIC610 PCR product. The two fragments were then gel purified, ligated to each other using the NEBuilder Hifi DNA Assembly Master Mix (NEB, Evry, France) and the ligation product was used to transform StellarTM E. coli competent cells to AmpR. Transformants encoding the mCherryBSU-linker-PEPut fusion from Physpank were screened by colony PCR. The plasmid DNA was extracted from three representative positive clones and checked by DNA sequencing. A validated construct was named pLJH225-Cm and used to transform B. subtilis competent cells deleted for the natural pykA gene (DGRM25). A strain containing the fusion cassette inserted by a double crossover at the amyE locus was selected by PCR and named DGRM1110.

The construction of pLJH227 which contains the PykA-linker-mCherryBSU fusion was carried as follows. First, we restricted the pDR111 DNA with NheI and HindIII and generated with the Q5 enzyme two PCR products, one containing the pykA ORF and its RBS region, the other the mCherryBSU ORF fused in 5′ to the linker. The DNA fragments (containing appropriate 20 bp sequence homologies at their ends) were then gel purified and ligated using the NEBuilder Hifi DNA Assembly Master Mix. Upon transformation of StellarTM competent cells with the ligation mixture, selection on Amp containing plates, colony PCR screening, plasmid extraction, and DNA sequencing, we obtained pLJH227. To construct pLJH226, a similar strategy was followed using a PCR fragment containing the pykAΔPEPut ORF (and its RBS region) instead of the entire pykA ORF (and its RBS region). For generating pLJH226-Cm, the BglII fragment of pLJH226 containing the SpR marker and pUC sequences was replaced by the corresponding restriction fragment of pIC610 which contains a CmR marker in place of the SpR marker. The pLJH226-Cm and pLJH227 DNA were then used to transformed DGRM25-Km (ΔpykA) competent cells and representative transformants containing the fusion cassette inserted by a double crossover at the amyE locus were selected by PCR and named DGRM1109 (construction encoding the PykA-linker-mCherryBSU fusion) or DGRM1111 (construction encoding the PykAΔPEPut-linker-mCherryBSU fusion).

Growth conditions

Routinely, B. subtilis and E. coli cells were grown at 37 °C in LB with or without antibiotics at the following concentrations: spectinomycin (Sp, 60 μg/mL); tetracycline (Tet, 7.5 μg/ml); chloramphenicol (Cm, 5 μg/mL); phleomycin (Phl, 7.5 μg/mL); kanamycin (Kan, 5 μg/mL); ampicillin (Amp, 50 μg/mL). Malate (0.4% w/v) was added to LB for routine cultivation of pykA mutants. Other studies were carried out with cells grown at 37 °C in a minimal medium (K2HPO4: 80 mM; KH2PO4: 44 mM; (NH4)2SO4: 15 mM; C6H5Na3O7 2H20: 3, 4 mM; CaCl2: 50 mM; MgSO4: 2 mM; FeIII citrate: 11 μg/mL; MnCl2: 10 μM; FeSO4: 1 μM; FeCl3: 4 μg/mL; tryptophan: 0.01% w/v) supplemented with various carbon sources as follows (in w/v): MC medium: malate (0.4%) + enzymatic casein hydrolysate (casa, 0.02%); GC medium: glucose (0.4%) + enzymatic casein hydrolysate (casa, 0.02%); M medium: malate (0.4%); G medium: glucose (0.4%); C medium: casa (0.02%).

Quantitative PCR

For monitoring ori/ter ratios, 6–14 cultures inoculated from generally three independent constructs were first grown overnight at 30 °C in MC supplemented with antibiotic. Upon saturation, cultures were diluted 1000-fold in the same medium without antibiotic and growth at 37 °C was carefully monitored using spectrophotometry. Samples for qPCR analysis were collected at low cell concentration (OD600 nm = 0.06 to 0.15) to ensure that cell cycle parameters are determined in steady-state cells and are not affected by the approach to the stationary phase or by changes in medium composition. The genomic DNA was extracted as described previously [44] or using the PureLink Genomic DNA mini kit (Invitrogen by Thermo Fisher Scientific, Courtaboeuf, Fr). Every qPCR reaction was carried out using two technical repeats of 4 serial dilutions. A non-replicating control DNA (stage II sporlets, [114]) was analyzed simultaneously with the samples in about 1/4 of the qPCR plates. Reactions and amplifications were carried out as previously described in 1× SYBR qPCR Premix Ex Taq (Tli RNaseH Plus) (Ozyme, St Quentin en Yvelines, France) mix and on a Mastercycler®ep realplex (Eppendorf, Le Pecq, Fr) [44]. Ratios were normalized using the mean of about 125 measures of the non-replicating control DNA (0.5885 ± 0.006) and compared using the nonparametric Mann-Whitney U test ( run at a significance level of p < 0.01 and a two-tailed hypothesis.

C periods were determined from three independent cultures and using 10 pairs of primers arranged regularly (from oriC to terC) along the right arm of the chromosome, as previously described [44]. Mean fork velocity was calculated using the C period (min) and the actual size of the TF8A genome (4,000,631 bp; this genome is deleted for prophages SPβ, PBSX, and skin).

Flow cytometry analysis

Strains were grown as indicated in the previous section, and at OD600nm = 0.1–0.15, chloramphenicol (200 μg/mL) was added to the cultures to impede replication initiation and cell division and allow completion of ongoing rounds of replication [115]. After 4 h of drug treatment, 1.5 mL of cells was fixed in filtered ethanol 70% v/v and stored at 4 °C. Stored cells were then washed twice in 1 mL of filtered Tris-buffered saline buffer (TBS 150) (20 mM Tris-HCl pH 7.5, 150 mM NaCl) and stained with Hoechst 33258 (1.5 μg/mL) for at least 30 min, as described elsewhere [97]. Flow cytometry analysis was carried out using a MoFlow Astrios cell sorter (Beckman Coulter, Life Sciences) equipped with a 355-nm krypton laser and a 448/59-nm bandpass filter used to collect Hoechst 33258 fluorescence data. Data were analyzed with the Kaluza software (Beckman Coulter, Life Sciences). We counted 100,000 events. Cell cycle parameters were determined as previously described [97].

Pyruvate kinase activity measurement in crude extracts

Cells (25 mL; OD600nm = 0.3) growing exponentially in MC were carefully collected by centrifugation (7 min; 4,922 RFC; room temperature) and resuspended in 75 μL of lysis buffer (Na2HPO4 2H2O: 60 mM; NaH2PO4 4H2O: 4 mM; KCl: 10 mM; MgSO4 7H2O: 1 mM; DTT: 1mM; Lyzozyme: 0.1 mg/mL; DNase I: 40 U/mL). They were then incubated 20 min on ice, 5 min at 37 °C, and 15 min at room temperature. Crude extracts were then collected by centrifugation (10 min; 14,000 rpm; 4 °C). The PykA activity was determined using the colorimetric/fluorometric assay kit K709-100 (CliniScience, Nanterre, Fr) and fluorescence (Ex/Em = 535/587 nm) was assessed using a ClarioStar apparatus (BMG Labtech, Champigny-sur-Marne, Fr). Protein concentration was monitored using the standard Bradford assay.


Cells were grown in MC medium as for qPCR analysis. For DNA segregation analysis, 1.5 mL of cells at OD600nm = 0.15 were briefly (2 min) incubated in the presence of 1 μg/mL 4’,6’-diamidino-2-phenylindole (DAPI, Sigma-Aldrich, Saint-Quentin-Falavier, France) and/or 2 μg/mL FM4-94 (Invitrogen, Villebon-sur-Yvette, France) for nucleoid and membrane staining, respectively. They were then collected by gentle centrifugation at room temperature (11,000 rpm for 1 min), resuspended in 10 μL MC medium and immediately mounted on a glass slide covered with a 1.0% w/v agarose pad in MC. A 0.17-mm glass coverslip was then placed on top of the agarose pad. For mCherryBSU fusion localization, cells prepared as above were collected by gentle centrifugation at room temperature (11,000 rpm for 1 min) and resuspended in MC or in a PBS solution containing 1% v/v paraformaldehyde (Alfa Aesar by Thermo Fisher Scientific, Kandel, Germany) (750 μL). Cells treated with paraformaldehyde were incubated 10 min at 4 °C. Samples treated or not with paraformaldehyde were then stained with DAPI, collected, and resuspended in 10 μL MC medium as above. Microscopy was carried out using an epifluorescence microscope (Zeiss, Axio Observer.Z1) with a × 100 magnification oil-immersion objective (Plan-APOCHROMAT Ph3) and a CMOS camera (Orca Flash 4.0 LT Hamamatsu). Digital images were acquired and analyzed using the Zen 2.6 (blue edition) software.

Mass spectrometry

Metabolome preparation

Metabolite extraction was carried out from cells grown at 37 °C and an OD600 = 0.4 using a protocol adapted from the Metabolomics Service Protocols from the University of Glasgow ( and modified as previously described [116]. For comparative analysis, pykAT>D and wild-type cells were grown on solid medium and metabolomes were prepared as previously [94].

Chromatographic conditions

Analyses were conducted using a Dionex Ultimate 3000 Rapid Separation LC (Thermo Fisher Scientific). The method for metabolome analysis was previously reported [94]. For hydrogen/deuterium exchange (HDX) experiments, chromatographic conditions were also described [116].

HRMS analyses

High-resolution measurements were obtained with a VelosPro Orbitrap Elite mass spectrometer (Thermo Fisher Scientific) fitted with a heated electrospray ionization source operating in positive and negative ionization modes. The mass spectrometer settings were as previously reported [116]. Collision-induced dissociation spectra (CID) in resonant excitation conditions were acquired using data-dependent scanning function for identification purpose. Non-targeted CID experiments using data-dependent scan were conducted at 27% Normalized Collision Energy (NCE) in positive and negative modes. CID experiments to confirm the presence of Pse (Leg) in metabolome were conducted at 22% NCE in positive and negative modes.

Metabolomic data processing

LC/MS raw data were converted to mzXML using MSConvert 3.0.20338 [117] and processed with R package XCMS [118, 119] with an “in-house” modified KEGGREST (Tenenbaum and Maintainer, 2021) R function for automatic compound annotation. Multicore processing was made possible through R package BiocParallel (Morgan et al., 2021). Multivariate analyses were conducted using muma for data pretreatment (pareto scaling) and FactominR for PCA visualization [120, 121]. Metabolite identification was aided by MS2 interpretation and retention time matching to commercial reference standards. Metabolomics raw data are available in MetaboLights ( [122].

Protein biochemistry

Replication enzymes

Replication enzymes DnaC, DnaG, DnaE, DnaN, HolB, YqeN, and DnaX were purified and tested as described previously [43, 49].

Helicase assays

Briefly, DnaC helicase assays were carried out with 633 nM (referring to monomer) each DnaC and DnaI, in 50 mM NaCl, 50 mM Tris-HCl pH 7.5, 10 mM MgCl2, 1 mM DTT, 2.5 mM ATP, and 2 nM DNA substrate, in the presence and absence of DnaG (313 nM) and/or PykA (633 nM monomer). The DNA substrate was constructed by annealing a 5′-32P radioactively labelled 100mer oligonucleotide (5′-CACACACACACACACACACACACACACACACACACACACACACACACACACACACACACCCCTTTAAAAAAAAAAAAAAAAGCCAAAAGCAGTGCCAAGCTTGCATGCC-3′) onto M13 ssDNA. The helicase activity was assayed by monitoring the displacement of the radioactive oligonucleotide from the M13 ssDNA through non-denaturing PAGE using 10% w/v polyacrylamide gels. Data were quantified using a Personal Molecular Imager with Quantity One 1-D analysis software (Bio-Rad) and analyzed with GraphPad Prism 4 software.

Polymerase assays

Time course polymerase assays were carried out by monitoring the DnaE primer extension activity using a 5′-32P radioactively labelled 20mer (5′-CAGTGCCAAGCTTGCATGCC-3′) or 60mer (5′-CAGTGCCAAGCTTGCATGCCTGCAGGTCGACTCTAGAGGATCCCCGGGTACCGAGCTCGA-3′) annealed to M13 ssDNA substrates (2 nM), in 50 mM Tris-HCl pH 7.5, 50 mM NaCl, 10 mM MgCl2, 1 mM DTT, 1 mM dNTPs, and 10 nM DnaE in the presence or absence of 40 nM PykA tetramer or 40 nM BSA. In some reactions with a lower, suboptimal 2 nM DnaE, 8 nM PykA tetramer was used. Nascent DNA was resolved through alkaline gel electrophoresis, as described before [43, 49]. Visualization and quantification were carried out using a Personal Molecular Imager with Quantity One 1-D analysis software (Bio-Rad) and data were analyzed with GraphPad Prism 4 software.

The effect of purified PEPut on the DnaE activity was investigated using a short DNA substrate comprising a 5′-32P radioactively 15mer oligonucleotide annealed onto a 110mer oligonucleotide as explained in Fig. S11.

Cloning, expression, and purification of the B. subtilis PykA

Cloning of pykA

A DNA fragment of 1755 bp carrying the B. subtilis pykA gene was amplified from genomic DNA using the PyKAF (5′-TACTTCCAATCCAATGCAAGAAAAACTAAAATTGTTTGTACCATCG-3′) and PykAR (5′-TTATCCACTTCCAATGTTATTAAAGAACGCTCGCACG-3′) forward and reverse primers, respectively in colony PCR reactions using Q5 high-fidelity DNA polymerase. Typically, a B. subtilis single colony was suspended in 20 mL LB and grown at 37 °C until the optical density reached 0.4–0.8. Thereafter, colony PCR reactions were carried out with genomic DNA (10 μL) acting as the template at 10-fold dilution. The PCR reactions were carried out in a volume of 50 μL with 1 unit Q5 high-fidelity DNA polymerase, 0.5 μM PykAF and PykAR, and 0.25 mM dNTPs in 1XQ5 polymerase buffer. PCR products were cleaned up with the Clean-up kit (Thermo Scientific), resolved by agarose electrophoresis, gel extracted using the GeneJET Gel Extraction kit (Thermo Scientific) and cloned into the p2CT plasmid (gift by James Berger) using ligation independent cloning to construct the p2CT-BsuPykA expression vector. This vector codes for an N-terminally His-tagged/MBP PykA protein with the His-MBP tag removable by proteolysis using the TEV protease.

Expression of PykA

For heterologous expression of the B. subtilis PykA, the p2CT-BsuPykA expression vector was transformed into Rosetta(DE3) E. coli. Single colonies were used to inoculate two 600 mL 2xYT cultures, supplemented with 60 μL carbenicillin (50 mg/mL) in 2-Lt conical flasks. The flasks were incubated at 37 °C, with shaking (180 rpm) until the optical density reached 0.6–0.8. Expression of PykA was induced by the addition of 0.5 mM IPTG and a further 3 h of growth the cells were harvested by centrifugation at 3000g for 15 min. Cells were suspended in 30 mL of buffer A (500 mM NaCl, 50 mM Tris-HCl pH 7.5, 20 mM imidazole) supplemented with 1 mM PMSF and 50 μL protease inhibitor cocktail (Fischer), sonicated at 15 amplitude microns for 1 min, 4 times with 1 min intervals on ice in between. Then, benzonase (20 μL) was added to the cell lysate which was further clarified at 40,000g for 40 min. The soluble crude extract was clarified and filtered through a 0.22-μm filter.

Purification of PykA

The PykA protein was purified from the filtered crude extract using a combination of IMAC (Immobilized Metal Affinity Chromatography) and gel filtration. First, the filtered crude extract was loaded onto a 5-mL HisTrap HP column (GE Healthcare) equilibrated in buffer A. The column was washed thoroughly with buffer A and the PykA protein was eluted using gradient elution with buffer B (500 mM NaCl, 50 mM Tris-HCl pH 7.5, 1 M imidazole). The eluted PykA protein was collected and quantified spectrophotometrically (extinction coefficient 76,780 M−1 cm−1). TEV protease was added at 1:20 molar ratio while dialyzing the protein solution overnight in dialysis buffer (500 mM NaCl, 50 mM Tris-HCl pH 7.5) at 4 °C in order to remove the His-MBP tag. The untagged PykA protein was then loaded back onto a 5 mL HisTrap HP column equilibrated in buffer A and the flow-through containing the untagged PykA was collected. Finally, the PykA protein solution was spun concentrated to 5–7 mL using a vivaspin 10 kDa cut-off filter. EDTA was added to 1 mM and the PykA was then loaded onto a HiLoad 26/60 Superdex 200 Prep Grade gel filtration column (GE Healthcare) equilibrated in buffer C (500 mM NaCl, 50 mM Tris-HCl pH 7.5, and 1 mM EDTA). Fractions containing the PykA protein were pooled and the protein was quantified spectrophotometrically (extinction coefficient 8940 M−1 cm−1), aliquoted and stored in – 80 °C.

PykA activity assay

The activity of purified PykA was assayed coupling the PykA catalyzed reaction (conversion of phosphoenolpyruvate to pyruvate) to the conversion of pyruvate into lactate catalyzed by LDH (Lactate Dehydrogenase) in the presence of NADH at 25 °C. The oxidation of NADH to NAD was followed spectrophotometrically at 340 nm and this is directly proportional to the activity of PykA. The LDH-catalyzed reaction was first optimized to ensure that it does not become a limiting factor when measuring the activity of PykA. The PykA catalyzed reactions were carried out in 96-well plates using the reaction master mix (10 mM Tri-HCl pH 7.5, 10 mM MgCl2, 50 mM KCl, 0.5 mM NADH, 2 mM ADP, 9.375 × 10−4 mg/mL LDH, and 5.7 μg/mL PykA) at 25 °C, at increasing PEP concentrations (0, 0.001, 0.05, 0.1, 0.25, 0.5, 0.75, 1, 2, 4, 8, and 10 mM). Data were analyzed using GraphPad Prism 4 software to plot a Hill plot and a Michaelis-Menten plot from where Vmax, Km, and n (the Hill coefficient) were calculated using GraphPad Prism 4 software.

Characterization of the oligomeric state of B. subtilis PykA

PykA assembles into a functional tetramer. The oligomeric state of the B. subtilis PykA was assessed by native mass spectrometry (MS) and gel filtration. Purified PykA protein was buffer exchanged into 200 mM ammonium acetate (pH 6.8) with a 75-μL Zeba (7k MWCO) desalting column. Protein (4 μM) was loaded onto a borosilicate emitter tip (prepared in-house) and back-fitted with a platinum wire electrode. Native mass spectra were acquired by nano-electrospray ionization on a modified Q-Tof instrument (Micromass) with 4.8 mbar backing pressure, capillary voltage of 1.3 kV, cone 30 V, source temperature 50 °C, collision 20 V, mass range 1500–15000 m/z, and a 32k quad profile for transmitting high m/z ions. The instrument was externally calibrated with NaCI (2 ng/μL). Mass spectra were processed with MassLynx v4.1 with minimal smoothing and background subtraction. Deconvoluted zero-charge spectra were generated using UniDec v4.2.2 with settings for low-resolution Q-ToF spectra.

The native MS spectrum revealed a very clean and extremely stable tetramer with miniscule amounts of dimer and monomer (Fig. S8E-F). Equally, comparative analytical gel filtration was consistent with a PykA tetramer as it eluted before the γ-globulin (158 kDa) and not as a monomer (62,314 Da) (Fig. S8G-H).

Cloning, expression, and purification of the PEPut domain of B. subtilis PykA

Cloning of the PEPut domain

The DNA fragment coding for the PEPut domain, with the preceding 10 amino acids was isolated by PCR using genomic DNA and the pepF (5′-TACTTCCAATCCAATGCAGCACAAAATGCAAAAGAAGCT-3′) and pepR (5′-TTATCCACTTCCAATGTTATTAAAGAACGCTCGCACG-3′) primers and cloned into the p2CT plasmid, as described above for PykA. The resulting p2CT-PEPut expression construct produced an N-terminally His-MBP tagged PEPut protein. The His-MBP tag was removable by proteolysis using the TEV protease.

Expression and purification of PEPut

Expression and purification of the PEPut were carried out as described for the full length PykA protein but the last gel filtration step during purification was omitted. Quantification of the final untagged PEPut (MW 9582.8 Da) was carried out spectrophotometrically using the extinction coefficients 69,330 M−1 cm−1 (for the His-MBP tagged PEPut) and 1490 M−1 cm−1 (for the untagged PEPut after TEV protease treatment). Purified PEPut is shown in Fig. S11).

Availability of data and materials

All data generated and analyzed in this study are available, either in this published article or the supplementary material or in the MetaboLights repository ( [122]. Strains and plasmids are available upon request from the authors as well as sequence details.


  1. Siddiqui K, On KF, Diffley JFX. Regulating DNA replication in eukarya. Cold Spring Harb Perspect Biol. Cold Spring Harbor Lab. 2013;5:a012930.

    Article  CAS  Google Scholar 

  2. Jameson K, Wilkinson AJ. Control of initiation of DNA replication in Bacillus subtilis and Escherichia coli. Genes. 2017;8:22–32.

    Article  CAS  PubMed Central  Google Scholar 

  3. Bipatnath M, Dennis PP, Bremer H. Initiation and velocity of chromosome replication in Escherichia coli B/r and K-12. J Bacteriol. 1998;180(2):265–73.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. Helmstetter CE. Timing of synthetic activities in the cell cycle. In: Neidhart FC, Curtis RI, Ingraham EC, Lin KB, editors. (null) LECC Escherichia coli and Salmonella. Washington DC: ASM Press; 1996. p. 1591–605.

    Google Scholar 

  5. Schaechter M, Maaloe O, O KN. Dependency on medium and temperature of cell size and chemical composition during balanced growth of Salmonella typhimurium. J Gen Microbiol. 1958;19:592–606.

    CAS  Article  PubMed  Google Scholar 

  6. Sharpe ME, Hauser PM, Sharpe RG, Errington J. Bacillus subtilis cell cycle as studied by fluorescence microscopy: constancy of cell length at initiation of DNA replication and evidence for active nucleoid partitioning. J Bacteriol. 1998;180(3):547–55.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. Buchakjian MR, Kornbluth S. The engine driving the ship: metabolic steering of cell proliferation and death. Nat Rev Microbiol. 2010;11:715–27.

    CAS  Google Scholar 

  8. Burnetti AJ, Aydin M, Buchler NE. Cell cycle Start is coupled to entry into the yeast metabolic cycle across diverse strains and growth rates. Mol Biol Cell. 2015;27:64–74.

    Article  CAS  PubMed  Google Scholar 

  9. Ewald JC. How yeast coordinates metabolism, growth and division. Curr Opin Microbiol. 2018;45:1–7.

    CAS  Article  PubMed  Google Scholar 

  10. Klevecz RR, Bolen J, Forrest G, Murray DB. A genomewide oscillation in transcription gates DNA replication and cell cycle. Proc Natl Aca Sci USA. 2004;101(5):1200–5.

    CAS  Article  Google Scholar 

  11. Papagiannakis A, Niebel B, Wit EC, Heinemann M. Autonomous metabolic oscillations robustly gate the early and late cell cycle. Mol Cell. 2017;65:285–95.

    CAS  Article  PubMed  Google Scholar 

  12. Tu BP, Kudlicki A, Rowicka M, McKnight SL. Logic of the yeast metabolic cycle: temporal compartmentalization of cellular processes. Science. 2005;310:1152–8.

    CAS  Article  PubMed  Google Scholar 

  13. Yu F-X, Dai R-P, Goh S-R, Zheng L, Luo Y. Logic of a mammalian metabolic cycle: an oscillated NAD+/NADH redox signaling regulates coordinated histone expression and S-phase progression. Cell Cycle. 2009;8(5):773–9.

    CAS  Article  PubMed  Google Scholar 

  14. Flåtten I, Fossum-Raunehaug S, Taipale R, Martinsen S, Skarstad K. The DnaA protein is not the limiting factor for initiation of replication in Escherichia coli. PLoS Genet. 2015;11(6):e1005276.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Murray H, Koh A. Multiple regulatory systems coordinate DNA replication with cell growth in Bacillus subtilis. PLoS Genet. 2014;10(10):e1004731.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Mathews CK. Deoxyribonucleotide metabolism, mutagenesis and cancer. Nat Rev Cancer. 2015;15:528–39.

    CAS  Article  PubMed  Google Scholar 

  17. Maya-Mendoza A, Moudry P, Merchut-Maya JM, Lee M, Strauss R, Bartek J. High speed of fork progression induces DNA replication stress and genomic instability. Nature. 2018;559:279–84.

    CAS  Article  PubMed  Google Scholar 

  18. Hu C-M, Tien S-C, Hsieh P-K, Jeng Y-M, Chang M-C, Chang Y-T, et al. High glucose triggers nucleotide imbalance through O-GlcNAcylation of key enzymes and induces KRAS mutation in pancreatic cells. Cell Metab. 2019;29:1334–49.

    CAS  Article  PubMed  Google Scholar 

  19. Lu M, Campbell JL, Boye E, Kleckner N. SeqA: a negative modulator of replication nitiation in E. coli. Cell. 1994;77:413–26.

    CAS  Article  PubMed  Google Scholar 

  20. Ishida T, Akimitsu N, Kashioka T, Hatano M, Kubota T, Ogata Y, et al. DiaA, a novel DnaA-binding protein, ensures the timely initiation of Escherichia coli chromosome replication. J Biol Chem. 2004;279(44):45546–55.

    CAS  Article  PubMed  Google Scholar 

  21. Baranska S, Glinkowska M, Herman-Antosiewicz A, Maciag-Dorszynska M, Nowicki D, Szalewska-Palasz A, et al. Replicating DNA by cell factories: roles of central carbon metabolism and transcription in the control of DNA replication in microbes, and implications for understanding this process in human cells. Microb Cell Fact. 2013;12:55.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. Boye E, Nordström K. Coupling the cell cycle to cell growth. EMBO Rep. 2003;4(8):757–60.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. Wang JD, Levin PA. Metabolism, cell growth and the bacterial cell cycle. Nat Rev Microbiol. 2009;7(11):822–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  24. Du Y-CN, Stillman B. Yph1p, an ORC-interacting protein: potential links between cell proliferation control, DNA replication, and ribosome biogenesis. Cell. 2002;109:835–48.

    CAS  Article  PubMed  Google Scholar 

  25. Beaufay F, Coppine J, Hallez R. When the metabolism meets the cell cycle in bacteria. Curr Opin Microbiol. 2021;60:104–13.

    CAS  Article  PubMed  Google Scholar 

  26. Fernández-Coll L, Maciag-Dorszynska M, Tailor K, Vadia S, Levin PA, Szalewska-Palasz A, et al. The absence of (p)ppGpp renders initiation of Escherichia coli chromosomal DNA synthesis independent of growth rates. mBio. 2020;11(2):e03223–19.

    Article  PubMed  PubMed Central  Google Scholar 

  27. DeNapoli J, Tehranchi AK, Wang JD. Dose-dependent reduction of replication elongation rate by (p)ppGpp in Escherichia coli and Bacillus subtilis. Mol Microbiol. 2013;88(11):93–104.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. Hernandez JV, Bremer H. Characterization of RNA and DNA synthesis in Escherichia coli strains devoid of ppGpp. J Biol Chem. 1993;268(15):10851–62.

    CAS  Article  PubMed  Google Scholar 

  29. Chubukov V, Gerosa L, Kochanowski K, Sauer U. Coordination of microbial metabolism. Nat Rev Microbiol. 2014;12:327–40.

    CAS  Article  PubMed  Google Scholar 

  30. Liu GY, Sabatini DM. mTOR at the nexus of nutrition, growth, ageing and disease. Nat Rev Mol Cell Biol. 2021;21:183–203.

    Article  CAS  Google Scholar 

  31. Hughes P, Landoulsi A, Kohiyama M. A novel role for cAMP in the control of the activity of the E. coli chromosome replication initiator protein, DnaA. Cell. 1988;55:343–50.

    CAS  Article  PubMed  Google Scholar 

  32. Maciąg M, Nowicki D, Jannière L, Szalewska-Pałasz A, Wegrzyn G. Genetic response to metabolic fluctuations: correlation between central carbon metabolism and DNA replication in Escherichia coli. Microb Cell Fact. 2011;10:19.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Maciag-Dorszynska M, Ignatowska M, Jannière L, Wegrzyn G, Szalewska-Pałasz A. Mutations in central carbon metabolism genes suppress defects in nucleoid position and cell division of replication mutants in Escherichia coli. Gene. 2012;503:31–5.

    CAS  Article  PubMed  Google Scholar 

  34. Tymecka-Mulik J, Boss L, Maciag-Dorszynska M, Matias Rodrigues JF, Gaffke L, Wosinski A, et al. Suppression of the Escherichia coli dnaA46 mutation by changes in the activities of the pyruvate-acetate node links DNA replication regulation to central carbon metabolism. PLoS One. 2017;12(4):e0176050.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Zhang Q, Zhou A, Li S, Ni J, Tao J, Lu J, et al. Reversible lysine acetylation is involved in DNA replication initiation by regulating activities of initiator DnaA in Escherichia coli. Sci Rep. 2016;6:30837.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. Krause K, Maciag-Dorszynska M, Wosinski A, Gaffke L, Morcinek-Orłowska J, Rintz E, et al. The role of metabolites in the link between DNA replication and central carbon metabolism in Escherichia coli. Genes. 2020;11:447.

    CAS  Article  PubMed Central  Google Scholar 

  37. Bergé M, Pezzatti J, Gonzalez-Ruiz V, Degeorges L, Mottet-Osman G, Rudaz S, et al. Bacterial cell cycle control by citrate synthase independent of enzymatic activity. eLife. 2020;9:e52272.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Laffan J, Firshein W. Membrane protein binding to the origin region of Bacillus subtilis. J Bacteriol. 1987;169(9):4135–40.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. Laffan J, Firshein W. Origin-specific DNA-binding membrane-associated protein may be involved in repression of initiation in Bacillus subtilis. Proc Natl Aca Sci USA. 1988;85:7452–6.

    CAS  Article  Google Scholar 

  40. Stein A, Firshein W. Probable identification of a membrane-associated repressor of Bacillus subtilis DNA replication as the E2 subunit of the pyruvate dehydrogenase complex. J Bacteriol. 2000;182(8):2119–24.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  41. Noirot-Gros M-F, Dervyn E, Wu LJ, Mervelet P, Errington J, Ehrlich SD, et al. An expanded view of bacterial DNA replication. Proc Natl Aca Sci USA. 2002;99(12):8342–7.

    CAS  Article  Google Scholar 

  42. Jannière L, Canceill D, Suski C, Kanga S, Dalmais B, Lestini R, et al. Genetic evidence for a link between glycolysis and DNA replication. PLoS One. 2007;5:e447.

    Article  CAS  Google Scholar 

  43. Paschalis V, Le Chatelier E, Green M, Képès F, Soultanas P, Jannière L. Interactions of the Bacillus subtilis DnaE polymerase with replisomal proteins modulate its activity and fidelity. Open Biol. 2017;7:170146.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Nouri H, Monnier A-F, Fossum-Raunehaug S, Maciag-Dorszynska M, Cabin-Flaman A, Képès F, et al. Multiple links connect central carbon metabolism to DNA replication initiation and elongation in Bacillus subtilis. DNA Res. 2018;25(6):641–53.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  45. Dervyn E, Suski C, Daniel R, Bruand C, Chapuis J, Errington J, et al. Two essential DNA polymerases at the bacterial replication fork. Science. 2001;294:1716–9.

    CAS  Article  PubMed  Google Scholar 

  46. Le Chatelier E, Becherel OJ, d'Alencon E, Canceill D, Ehrlich SD, Fuchs RPP, et al. Involvement of DnaE, the second replicative DNA polymerase from Bacillus subtilis, in DNA mutagenesis. J Biol Chem. 2004;279(3):1757–67.

    Article  CAS  PubMed  Google Scholar 

  47. Sanders GM, Dallmann HG, McHenry CS. Reconstitution of the B. subtilis replisome with 13 proteins including two distinct replicases. Mol Cell. 2010;37:273–81.

    CAS  Article  PubMed  Google Scholar 

  48. Soultanas P. Loading mechanisms of ring helicases at replication origins. Mol Microbiol. 2012;84(1):6–16.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  49. Rannou O, Le Chatelier E, Larson MA, Nouri H, Dalmais B, Laughton C, et al. Functional interplay of DnaE polymerase, DnaG primase and DnaC helicase within a ternary complex, and primase to polymerase hand-off during lagging strand DNA replication in Bacillus subtilis. Nucleic Acids Res. 2013;41(10):5303–20.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  50. Dickinson JR, Williams AS. The cdc30 mutation in Saccharomyces cerevisiae results in a temperature-sensitive isoenzyme of phosphoglucose isomerase. J Gen Microbiol. 1987;133:135–40.

    CAS  PubMed  Google Scholar 

  51. Fornalewicz K, Wieczorek A, Wegrzyn G, Łyżeń R. Silencing of the pentose phosphate pathway genes influences DNA replication in human fibroblasts. Gene. 2017;635:33–8.

    CAS  Article  PubMed  Google Scholar 

  52. Konieczna A, Szczepańska A, Sawiuk K, Wegrzyn G, Łyżeń R. Effects of partial silencing of genes coding for enzymes involved in glycolysis and tricarboxylic acid cycle on the enterance of human fibroblasts to the S phase. BMC Cell Biol. 2015;16:16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Sprague GFJ. Isolation and characterization of a Saccharomyces cerevisiae mutant deficient in pyruvate kinase activity. J Bacteriol. 1977;130(1):232–41.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  54. Wieczorek A, Fornalewicz K, Mocarski Ł, Łyżeń R, Wegrzyn G. Double silencing of relevant genes suggests the existence of the direct link between DNA replication/repair and central carbon metabolism in human fibroblasts. Gene. 2018;650:1–6.

    CAS  Article  PubMed  Google Scholar 

  55. Ronai Z. Glycolytic enzymes as DNA binding proteins. Int J Biochem Cell Biol. 1993;25(7):1073–6.

    CAS  Google Scholar 

  56. Sirover MA. New insights into an old protein: The functional diversity of mammalian glyceraldehyde-3-phosphate dehydrogenase. Biochim Biophys Acta. 1999;1432:159–84.

    CAS  Article  PubMed  Google Scholar 

  57. Sirover MA. On the functional diversity of glyceraldehyde-3-phosphate dehydrogenase: Biochemical mechanisms and regulatory control. Biochim Biophys Acta. 2011;1810:741–51.

    CAS  Article  PubMed  Google Scholar 

  58. Kim J-W, Dang CV. Multifaceted roles of glycolytic enzymes. Trends Biochem Sci. 2005;30(3):142–50.

    CAS  Article  PubMed  Google Scholar 

  59. Konieczna A, Szczepańska A, Sawiuk K, Łyżeń R, Wegrzyn G. Enzymes of the central carbon metabolism: are they linkers between transcription, DNA replication, and carcinogenesis? Med Hypotheses. 2015;84:58–67.

    CAS  Article  PubMed  Google Scholar 

  60. Boukouris AE, Zervopoulos SD, Michelakis ED. Metabolic enzymes moonlighting in the nucleus: metabolic regulation of gene transcription. Trends Biochem Sci. 2016;41(8):712–30.

    CAS  Article  PubMed  Google Scholar 

  61. Lu Z, Hunter T. Metabolic kinases moonlighting as protein kinases. Trends Biochem Sci. 2018;43(4):301–10.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  62. Snaebjornsson MT, Schulze A. Non-canonical functions of enzymes facilitate cross-talk between cell metabolic and regulatory pathways. Exp Mol Med. 2018;50:34.

    Article  CAS  PubMed Central  Google Scholar 

  63. Cai L, Sutter BM, Li B, Tu BP. Acetyl-CoA induces cell growth and proliferation by promoting the acetylation of histones at growth genes. Mol Cell. 2011;42:426–37.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  64. Sutendra G, Kinnaird A, Dromparis P, Paulin R, Stenson TH, Haromy A, et al. A nuclear pyruvate dehydrogenase complex is important for the generation of Acetyl-CoA and histone acetylation. Cell. 2014;158:84–97.

    CAS  Article  PubMed  Google Scholar 

  65. Wellen HG, Sachdeva UM, Bui TV, Cross JR, Thompson CB. ATP-citrate lyase links cellular metabolism to histone acetylation. Science. 2009;324:1076–80.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  66. Zheng L, Roeder RG, Luo Y. S phase activation of the histone H2B promoter by OCA-S, a coactivator complex that contains GAPDH as a key component. Cell. 2003;114:255–66.

    CAS  Article  PubMed  Google Scholar 

  67. Dai RP, Yu FX, Goh SR, Chng HW, Tan YL, Fu JL, et al. Histone 2B (H2B) expression is confined to a proper NAD+/NADH redox status. J Biol Chem. 2008;283(4):26894–901.

    CAS  Article  PubMed  Google Scholar 

  68. Ma R, Wu Y, Zhai Y, Hu B, Ma W, Yang W, et al. Exogenous pyruvate represses histone gene expression and inhibits cancer cell proliferation via the NAMPT–NAD+–SIRT1 pathway. Nucleic Acids Res. 2019;47(21):11132–50.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  69. Li X, Qian X, Jiang H, Xia Y, Zheng Y, Li J, et al. Nuclear PGK1 alleviates ADP-dependent inhibition of CDC7 to promote DNA replication. Mol Cell. 2018;72:650–60.

    CAS  Article  PubMed  Google Scholar 

  70. Grosse F, Nasheuer H-P, Scholtissek S, Schomburg U. Lactate dehydrogenase and glyceraldehyde-phosphate dehydrogenase are single-stranded DNA-binding proteins that affect the DNA-polymerase-α–primase complex. Eur J Biochem. 1986;160:459–67.

    CAS  Article  PubMed  Google Scholar 

  71. Popanda O, Fox G, Thielmann HW. Modulation of DNA polymerases alpha, delta and epsilon by lactate dehydrogenase and 3-phosphoglycerat kinase. Biochim Biophys Acta. 1998;1397:102–17.

    CAS  Article  PubMed  Google Scholar 

  72. Jindal HK, Vishwanatha JK. Functional identity of a primer recognition protein as phosphoglycerate kinase. J Biol Chem. 1990;265(12):6540–3.

    CAS  Article  PubMed  Google Scholar 

  73. Schormann N, Hayden KL, Lee P, Banerjee S, Chattopadhyay D. An overview of structure, function, and regulation of pyruvate kinases. Protein Sci. 2019;28(1771):1784.

    Google Scholar 

  74. Suzuki K, Ito S, Shimizu-Ibuka A, Sakai H. Crystal structure of pyruvate kinase from Geobacillus stearothermophilus. J Biochem. 2008;144(3):305–12.

    CAS  Article  PubMed  Google Scholar 

  75. Morgan HP, Zhong W, McNae IW, Michels PAM, Fothergill-Gilmore LA, Walkinshaw MD. Structures of pyruvate kinases display evolutionarily divergent allosteric strategies. R Soc Open Sci. 2014;1:140120.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Nicolas P, Mäder U, Dervyn E, Rochat T, Leduc A, Pigeonneau N, et al. Condition-dependent transcriptome reveals high-level regulatory architecture in Bacillus subtilis. Science. 2012;335:1103–6.

    CAS  Article  PubMed  Google Scholar 

  77. Muntel J, Fromion V, Goelzer A, Maab S, Mäder U, Büttner K, et al. Comprehensive absolute quantification of the cytosolic proteome of Bacillus subtilis by data independent, parallel fragmentation in liquid chromatographis/mass spectrometry. Mol Cell Proteomics. 2014;13(4):1008–19.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  78. Sakai H. Possible structure and function of the extra C-terminal sequence of pyruvate kinase from Bacillus stearothermophilus. J Biochem. 2004;136:471–6.

    CAS  Article  PubMed  Google Scholar 

  79. Nguyen CC, Saier MH Jr. Phylogenetic analysis of the putative phosphorylation domain in the pyruvate kinase of Bacillus stearothermophilus. Res Microbiol. 1995;146:713–9.

    CAS  Article  PubMed  Google Scholar 

  80. Alpert CA, Frank R, Stüber K, Deutscher J, Hengstenberg W. Phosphoenolpyruvate-dependent protein kinase enzyme I of Streptococcus faecalis: purification and properties of the enzyme and characterization of its active center. Biochemistry. 1985;24:959–64.

    CAS  Article  PubMed  Google Scholar 

  81. Teplyakov A, Lim K, Zhu PP, Kapadia G, Chen CCH, Schwartz J, et al. Structure of phosphorylated enzyme I, the phosphoenolpyruvate:sugar phosphotransferase system sugar translocation signal protein. Proc Natl Aca Sci USA. 2006;103(44):16218–23.

    CAS  Article  Google Scholar 

  82. Goss NH, Evans CT, Wood HG. Pyruvate phosphate dikinase: sequence of the histidyl peptide, the pyrophosphoryl and phosphoryl carrier. Biochemistry. 1980;19:5805–9.

    CAS  Article  PubMed  Google Scholar 

  83. Herzberg O, Chen CC, Kapadia G, McGuire M, Carroll LJ, Noh SJ, et al. Swiveling-domain mechanism for enzymatic phosphotransfer between remote reaction sites. Proc Natl Aca Sci USA. 1996;93:2652–7.

    CAS  Article  Google Scholar 

  84. Tolentino R, Chastain C, Burnell J. Identification of the amino acid involved in the regulation of bacterial pyruvate, orthophosphate dikinase and phosphoenolpyruvate synthetase. Adv Biol Chem. 2013;03:12–21.

    Article  Google Scholar 

  85. Burnell JN, Chastain CJ. Cloning and expression of maize-leaf pyruvate, Pi dikinase regulatory protein gene. Biochem Biophys Res Commun. 2006;345:675–80.

    CAS  Article  PubMed  Google Scholar 

  86. Burnell JN. Cloning and characterization of Escherichia coli DUF299: a bifunctional ADP-dependent kinase - Pi-dependent pyrophosphorylase from bacteria. BMC Biochem. 2010;11:1–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. Burnell JN, Hatch MD. Regulation of C4 photosynthesis: Identification of a catalycally important histidine residue and its role in the regulation of pyruvate. Pi dikinase Arch Biochem Biophys. 1984;231(1):175–82.

    CAS  Article  PubMed  Google Scholar 

  88. Eymann C, Dreisbach A, Albrecht D, Bernhardt J, Becher D, Gentner S, et al. A comprehensive proteome map of growing Bacillus subtilis cells. Proteomics. 2004;4:2849–76.

    CAS  Article  PubMed  Google Scholar 

  89. Mäder U, Schmeisky AG, Flórez LA, Stülke J. SubtiWiki--a comprehensive community resource for the model organism Bacillus subtilis. Nucleic Acids Res. 2012;40(D1):D1278–87.

    Article  CAS  PubMed  Google Scholar 

  90. Reiland S, Messerli G, Baerenfaller K, Gerrits B, Endler A, Grossmann J, et al. Large-scale Arabidopsis phosphoproteome profiling reveals novel chloroplast kinase substrates and phosphorylation networks. Plant Physiol. 2009;150(2):889–903.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  91. Pisithkul T, Patel NM, Amador-Noguez D. Post-translational modifications as key regulators of bacterial metabolic fluxes. Curr Opin Microbiol. 2015;24:29–37.

    CAS  Article  PubMed  Google Scholar 

  92. Brunk E, Chang RL, Xia J, Hefzi H, Yurkovich JT, Kim D, et al. Characterizing posttranslational modifications in prokaryotic metabolism using a multiscale workflow. Proc Natl Aca Sci USA. 2018;115(43):11096–101.

    CAS  Article  Google Scholar 

  93. Bollenbach TJ, Mesecar AD, Nowak T. Role of lysine 240 in the mechanism of yeast pyruvate kinase catalysis. Biochemistry. 1999;38:9137–45.

    CAS  Article  PubMed  Google Scholar 

  94. Stuani L, Lechaplais C, Salminen AV, Sérugens B, Durot M, Castelli V, et al. Novel metabolic features in Acinetobacter baylyi ADP1 revealed by a multiomics approach. Metabolomics. 2014;10:1223–38.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  95. Dubois T, Krzewinski F, Yamakawa N, Lemy C, Hamiot A, Brunet L, et al. Genes encode an original legionaminic acid pathway required for crust assembly in Bacillus subtilis. mBio. 2020;11:e01153–20.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  96. Washington TA, Smith JL, Grossman AD. Genetic networks controlled by the bacterial replication initiator and transcription factor DnaA in Bacillus subtilis. Mol Microbiol. 2017;106(1):109–28.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  97. Morigen OI, Skarstad K. Growth rate dependent numbers of SeqA structures organize the multiple replication forks in rapidly growing Escherichia coli. Genes Cells. 2009;14:643–57.

    CAS  Article  PubMed  Google Scholar 

  98. Gao S, Bao J, Gu X, Xin X, Chen C, Ryu DDY. Substrate promiscuity of pyruvate kinase on various deoxynucleoside diphosphates for synthesis of deoxynucleoside triphosphates. Enzyme Microb Technol. 2008;43:455–9.

    CAS  Article  Google Scholar 

  99. Meile J-C, Wu LJ, Ehrlich SD, Errington J, Noirot P. Systematic localisation of proteins fused to the green fluorescent protein inBacillus subtilis: Identification of new proteins at the DNA replication factory. Proteomics. 2006;6(7):2135–46.

    CAS  Article  PubMed  Google Scholar 

  100. Bruck I, O'Donnell M. The DNA replication machine of a gram-positive organism. J Biol Chem. 2000;275(37):28971–83.

    CAS  Article  PubMed  Google Scholar 

  101. Chuang C, Prasanth KR, Nagy PD. The glycolytic pyruvate kinase is recruited directly into the viral replicase complex to generate ATP for RNA synthesis. Cell Host Microbe. 2017;22:639–52.

    CAS  Article  PubMed  Google Scholar 

  102. Pancholi V, Chhatwal GS. Housekeeping enzymes as virulence factors for pathogens. Int J Med Microbiol. 2003;293:391–401.

    CAS  Article  PubMed  Google Scholar 

  103. Tu BP, Mohler RE, Liu JC, Dombek KM, Young ET, Synovec RE, et al. Cyclic changes in metabolic state during the life of a yeast cell. Proc Natl Aca Sci USA. 2007;104(43):16886–91.

    CAS  Article  Google Scholar 

  104. Ewald JC, Kuehne A, Zamboni N, Skotheim JM. The yeast cyclin-dependent kinase routes carbon fluxes to fuel cell cycle progression. Mol Cell. 2016;62:532–45.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  105. Hartl J, Kiefer P, Kaczmarczyk A, Mittelviefhaus M, Meyer F, Vonderach T, et al. Untargeted metabolomics links glutathione to bacterial cell cycle progression. Nat Metab. 2020;2:153–66.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  106. Lee J-B, Hite RK, Hamdan SM, Xie SX, Richardson CC, van Oijen AM. DNA primase acts as a molecular break in DNA replication. Nature. 2006;439(7076):621–4.

    CAS  Article  PubMed  Google Scholar 

  107. Tanner NA, Hamdan SM, Jergic S, Schaeffer PM, Dixon NE, van Oijen AM. Single-molecule studies of fork dynamics in Escherichia coli DNA replication. Nat Struct Mol Biol. 2008;15(2):170–6.

    CAS  Article  PubMed  Google Scholar 

  108. Yao NY, Georgescu RE, Finkelstein J, O'Donnell ME. Single-molecule analysis reveals that the lagging strand increases replisome processivity but slows replication fork progression. Proc Natl Aca Sci USA. 2009;106(32):13236–41.

    CAS  Article  Google Scholar 

  109. Georgescu RE, Yao N, Indiani C, Yurieva O, O'Donnell ME. Replisome mechanics: lagging strand events that influence speed and processivity. Nucleic Acids Res. 2014;42(10):6497–510.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  110. Rodgers K, McVey M. Error-Prone repair of DNA double-strand breaks. J Cell Physiol. 2016;231:15–24.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  111. Loeb LA. Human cancers express mutator phenotypes: origin, consequences and targeting. Nat Rev Cancer. 2011;11(6):450–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  112. Macheret M, Halazonetis TD. DNA replication stress as a hallmark of cancer. Annu Rev Pathol Mech Dis. 2015;10:425–48.

    CAS  Article  Google Scholar 

  113. Oskam L, Hillenga DJ, Venema G, Bron S. The large Bacillus plasmid pTB19 contains two integrated rolling-circle plasmids carrying mobilization functions. Plasmid. 1991;26:30–9.

    CAS  Article  PubMed  Google Scholar 

  114. Magill NG, Setlow P. Properties of purified sporlets produced by spoII mutants of Bacillus subtilis. J Bacteriol Am Soc Microbiol. 1992;174(24):8148–51.

    CAS  Article  Google Scholar 

  115. Séror SJ, Casarégola S, Vannier F, Zouari N, Dahl M, Boye E. A mutant cysteinyl-tRNA synthetase affecting timing of chromosomal replication initiation in B. subtilis and conferring resistance to a protein kinase C inhibitor. EMBO J. 1994;13(10):2472–80.

    Article  PubMed  PubMed Central  Google Scholar 

  116. Thomas M, Stuani L, Darii E, Lechaplais C, Pateau E, Tabet J-C, et al. De novo structure determination of 3-((3-aminopropyl)amino)- 4-hydroxybenzoic acid, a novel and abundant metabolite in Acinetobacter baylyi ADP1. Metabolomics. 2019;15:45.

    Article  CAS  PubMed  Google Scholar 

  117. Chambers MC, Maclean B, Burke R, Amode D, Ruderman DL, Neumann S, et al. A cross-platform toolkit for mass spectrometry and proteomics. Nat Biotechnol. 2012;30(10):918–20.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  118. Smith CA, Want EJ, OMaille G, Abagyan R, Siuzdak G. XCMS: Processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Anal Chem. 2006;78(3):779–87.

    CAS  Article  PubMed  Google Scholar 

  119. Tautenhahn R, Böttcher C, Neumann S. Highly sensitive feature detection for high resolution LC/MS. BMC Bioinformatics. 2008;9:504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  120. Gaude E, Chignola F, Spiliotopoulos D, Spitaleri A, Ghitti M, Garcia-Manteiga J, et al. muma, an R package for metabolomics univariate and multivariate statistical analysis. Curr Metabolomics. 2013;1:180–9.

    CAS  Article  Google Scholar 

  121. Lê S, Josse J, Husson F. FactoMineR: An R package for multivariate analysis. J Stat Softw. 2008;25(1):1–18.

    Article  Google Scholar 

  122. Haug. MetaboLights : a resource evolving in response to the needs of its scientific community. Nucleic Acids Res. 2020;48(D1):D440–4.

    CAS  PubMed  Google Scholar 

Download references


This work has benefited from the facilities and expertise of the France-BioImaging infrastructure supported by the Agence Nationale de la Recherche (ANR-10-INSB-04, call “investissements d’Avenir”). Strains BTS8 and LAS26 were generous gifts from Lyle A. Simmons. Plasmid p2CT, pIC610, and pMADΔspsABCDEF were gifts from James Berger, Dominique Le Coq and Thomas Dubois, respectively. We thank Ioana Popescu for assistance with the ClarioStar (BMG LABTECH, Fr), Olivier Espeli for help in cell fixation, and Dr Jedd Bellamy-Carter for helping us with the mass spectrometry experiments to assess the stoichiometry and stability of the PykA protein. We are grateful to Dominique Le Coq, Nathalie Declercq, and Thomas Dubois for stimulating discussions.


This work was supported by the Biotechnology Biological Sciences Research Council (BBSRC) grant (BB/R013357/1) to P.S., a BBSRC sub-contract (RIS 1165589) to L.J. and by recurring research funds from CNRS (Centre National de la Recherche Scientifique) and UEVE (Université d’Evry Val d’Essonne) to L.J.. S.H. was funded by a PhD fellowship from the MESR (Ministère de l’Enseignement Supérieur et de la Recherche; ED SDSV, Université Paris-Saclay, Université d’Evry Val d’Essonne). M.P. is a PhD student partially funded by a Vice Chancellor’s Excellence Award at the University of Nottingham. L.J. is on the CNRS staff. Funders had no role in the design of the study, the analysis and interpretation of data, and in the manuscript writing.

Author information




L.J. and P.S. designed the study; S.H. and L. J. conducted and analyzed in vivo experiments; M.P. and A.H. and P.S. conducted and analyzed in vitro experiments; E.P. conducted metabolic experiments; E.D. conducted mass spectrometry experiments for metabolic structure confirmation; C.L. conducted metabolic data pretreatment and statistical data analysis; E.P., E.D., C.L. and A.P. analyzed metabolic data; A.P. supervised metabolic experiments. L.J. and P.S. wrote the paper with the assistance of A.P. for the metabolomic sections. All the authors read and approved the final version of the manuscript.

Author’s information

Twitter handle: @ChemistryUoN (Panos Soultanas).

Corresponding authors

Correspondence to Panos Soultanas or Laurent Janniere.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors and funders have given their consent for publication.

Competing interests

The authors declare they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Fig. S1

. Key amino-acids of the Cat and PEPut domains of PykA. A. Cat domain analysis. Clustalw and Chimera analysis of the pyruvate kinase of B. subtilis (PykA), human cells (PKM2) and Mycobacterium tuberculosis (PYK) identified key amino acids of the catalytic site of the B. subtilis protein. B. PEPut domain analysis. Alignment of the PEPut domain of PykA to related domains of various metabolic enzymes. The red arrow highlights the conserved LTSH motif (coordinates 536-539). Fig. S2. Effect of Cat and PEPut mutations on growth in MC. Wild-type and pykA mutants were first grown over-night in MC supplemented with antibiotic when appropriate. Upon saturation, cultures were diluted 1000-fold in the same medium without antibiotic and growth was monitored spectrophometrically. Left panel: Analysis of catalytic mutants (pykAΔcat, pykAR32A, pykAR73A, pykAK220A, pykAGD245/6AA, pykAT278A, pykAJP). Right panel: Analysis of PEPut and Cat-PEPut interaction mutants (pykAΔPEP, pykAT>A, pykAS>A, pykAH>A, pykATSH>AAA, pykAT>D, pykAS>D, pykAH>D, pykATSH>DDD, pykAE209A, pykAL536A). Controls: TF8A (wild-type) and ΔpykA. Fig. S3. Analysis of NTP in the metabolome of wild-type and pykAT>D cells. ATP, GTP and CTP were detected in the positive ionization mode. UTP was detected in the negative ionization mode. Note that TTP signals were too low for quantifications. Data correspond to 3 independent extractions (solid cultures).*, p > 0.05 ; **, p < 0.05 (Welch's T-test). Values in bold indicate the fold change for each metabolite (WT vs pykAT>D). Fig. S4. LC/MS analysis of legionaminic acid in the metabolome. A. Extracted ion chromatogram (EIC) corresponds to the deprotonated molecule [M-H]- at m/z 333.1303 (5 ppm accuracy). B. Zoom on the mass spectrum of legionaminic acid in the negative mode. C. Collision Induced dissociation (CID) spectrum of legionaminic acid in the negative mode at 22% Normalized Collision Energy (NCE). D. CID spectrum of legionaminic acid in the positive mode at 22% NCE. E and F. Zoom on the mass spectrum of the deuterated forms of legionaminic acid in the negative and positive ionization mode, respectively. G and H. Comparison of legionaminic acid (G) and CMP-legionaminic (H) acid contents in wild-type (WT), ΔspsE and ΔspsF cells, respectively. Data correspond to 3 independent extractions (liquid cultures). ***, p < 0.001; **, p < 0.01 (Welch's T-test). Fig. S5. Representative cell cycle results in Cat and PEPut mutants. Top raw: Microscopy of exponentially growing cells stained with FM4-64 (membrane staining, red) and DAPI (nucleoid staining, blue). Middle raw: Representative runout DNA histograms (experiments were reiterated 3-12 times). Bottom raw: Representative marker frequency analysis along the right arm of the chromosome (experiments were reiterated at least three times). Fig. S6. Cell cycle parameters of wild-type and pykAT>D cells grown in proline and malate, respectively. Left panel: Growth in malate of wild-type (WT), ΔpykA (Δ), pykAT>D (T>D) and other pykA mutants (pykAT>A, pykAGD245/6AA and pykAK220A, blue lines). Right panel: Runout DNA histograms and cell cycle parameters. Fig. S7. PykA-mCherryBSU localization. Strains deleted for the natural pykA gene and encoding the PykA-mCherryBSU fusion from an inducible promoter (Physpank) were grown in MC and microscopy analysis was carried out at OD600nm = 0.1 to 0.2. Top raw: analysis of the mCherryBSU signal produced at different IPTG concentrations. Bottom raw: analysis of cells grown in the absence of IPTG, fixed in a 1x PBS solution supplemented with 1% paraformaldehyde and stained with DAPI. Scale bar: 4 μm. Similar results were obtained with fusions mutated in the Cat or PEPut domain of PykA. Fig. S8. PykA purification and characterization of its function and oligomeric state. A. SDS-PAGE (15% polyacrylamide gel) showing over-expression of the 6His-MBP-PykA in Rosetta (DE3) E. coli. The soluble expressed tagged PykA protein is shown in a red rectangular in lane 4, whereas lanes M, 1, 2 and 3 show molecular weight standards, the control uninduced insoluble fraction, the control uninduced soluble fraction and the IPTG-induced insoluble fraction, respectively. B.1. SDS-PAGE (15% polyacrylamide gel) showing fractions from the first IMAC purification step of PykA. From left to right, lanes represent molecular weight standards (M), the flow-through (1), the eluted tagged PykA (2), the overnight TEV treated tagged PykA (3), the flow-through fractions containing untagged PykA from the second IMAC step after TEV proteolysis (4-10). B.2. SDS-PAGE (15% polyacrylamide gel) showing the final gel filtration column (HiLoad 26/60 Superdex 200 Prep Grade Gel Filtration Column). From left to right, lanes represent molecular weight standards (M) and fractions of the size exclusion chromatography (1-9). C. The graph shows a Hill plot for the activity of PykA at 25°C. The Rate/(Vmax-Rate) (Y-axis) was plotted against the PEP substrate concentration (X-axis) using GraphPad Prism 4 software and the Vmax (19.3 μmol/min), Km (2.7 mM) and the Hill coefficient n (0.8111) values are shown below the graph. The n value is <1 indicating negative cooperative binding of PykA to its PEP substrate. D. The graph shows a Michaelis-Menten plot for the activity of PykA at 25°C. The initial rate of the reaction (Y-axis) was plotted against the PEP substrate concentration (X-axis) using GraphPad Prism 4 software and the Vmax (16.3 μmol/min) and Km (1.7 mM) values are shown below the graph. E. A native mass spectrum showing the PykA tetramer and miniscule amounts of the dimer and monomer. The theoretical mass of the PykA monomer (62,314.9 Da), dimer (124,629.8 Da) and tetramer (249,259.6 Da). Native mass spectrometry showed that PykA was found to be predominantly tetrameric (250,092 ± 72 Da), with very low abundance dimer (124,811 ± 38 Da) and monomer (62,437 ± 8 Da) peaks. F. Collision induced dissociation of the 33+ charge state of the tetramer shows it to be very stable in the gas-phase, with no apparent dissociation to lower-order oligomers. G. Comparative analytical gel filtration of the PykA tetramer against molecular weight standards (Thyroglobulin 670 kDa, g-globulin 158 kDa, ovalbumin 44 kDa, myoglobin 17 kDa and vitamin B12 1.3 kDa) through a Superdex 200 10/300 GL prepacked Tricorn gel filtration column (GE Healthcare). H. Selectivity trendline constructed from the molecular weight standards (shown in graph G) for the estimation of the PykA MW. The x axis is in logarithmic scale. Graphpad was used for plotting the data points. The theoretical value of our PykA (249,259.6 Da) is close to the estimated (285,000 Da) which along with the MS data verifies the tetramer in solution. Kd is the equilibrium distribution coefficient. The numbers (1-5) on the data points correspond to the proteins shown in graph G. Fig. S9. Stimulation of DnaE activity by PykA but not by BSA. A. Primer extension assays monitoring the extension of a 5′-32P-radioactively labelled 60mer DNA primer annealed onto M13 ssDNA over time by the B. subtilis DnaE. The activity of DnaE polymerase (10 nM) was monitored in the presence and absence of PykA (10 nM, tetramer) through a time course (30-150 sec). Lanes in the gels from left to right indicate: (M): DNA-ladder and then the time course (0, 30, 60, 90, 120 and 150 sec) depicted by the rectangular triangle. B. Primer extension assays as above with or without 10 nM (monomer) BSA instead of PykA. C. DnaE (1 nM) polymerase activity at increasing BSA concentrations (0, 5, 50, 500 nM), as indicated by the rectangular triangle, monitored by alkaline agarose electrophoresis. The DNA substrate is a labelled 20mer (5′-CAGTGCCAAGCTTGCATGCC-3′) primer annealed onto ssM13 ssDNA (2nM). The primer extension reaction was carried out for a longer time than above (5 min instead of 30-150 sec) and the film was over-exposed to compensate for the lower DnaE concentration. The assay was carried out at 37 °C in 50 mM Tris-HCl 7.5, 50 mM NaCl, 10 mM MgCl2 mM DTT, 1 mM dNTPs. No stimulation of the DnaE polymerase activity was observed in the presence of 5 and 50 nM BSA. The marginal stimulation observed at 500 nM BSA excess is likely because at this high concentration, BSA acts as a blocking agent preventing adhesion of DnaE to the plastic reaction tubes. Fig. S10. Stimulation of DnaE activity by PykA does not result from stimulation of DnaE binding to primed templates. EMSA investigation of the effect of PykA on the DNA binding of DnaE polymerase. The DNA substrate was constructed by annealing a 5′-32P-radioactively labelled 15mer (5′-AAGGGGGTGTGTGTG-3′) primer annealed onto a 30mer (5′-ACACACACACACACACACACACACCCCCTT-3′) oligonucleotide. Binding reactions were carried out with 1 nM DNA substrate, DnaE (500nM) and increasing concentrations (0, 12.5, 125 and 1,250 nM tetramer) of PykA, as indicated by the rectangular triangle for 10 min at 37°C in 50 mM NaCl, 10 mM MgCl2, 50 mM Tris-HCl pH 7.5. Lanes C and PykA represent the radioactive substrate in the absence of any proteins and in the presence of PykA (1,250 nM tetramer), respectively, showing that PykA does not bind to the DNA substrate. No stimulation of DnaE binding to DNA was observed in the presence of increasing concentrations of PykA indicating that PykA does not enhance the DNA binding activity of DnaE. Fig. S11. PEPut purification. A. SDS-PAGE showing overexpression of the His-MBP tagged PEPut in Rosetta (DE3) E. coli. From left to right, lanes show protein MW markers (M), the insoluble uninduced (1), soluble uninduced (2), insoluble induced (3) and soluble induced (4) fractions. The expressed soluble His-MBP tagged PEPut is shown by a red rectangular. B. SDS-PAGE showing the final purified untagged PEPut after removal of the His-MBP tag with TEV proteolysis. Lanes from left to right show protein MW markers (M) and fractions from the flow through the HisTrap column containing the pure untagged PEPut (lanes 1,2 and 3). Fig. S12. PEPut does not stimulate DnaE activity. Primer extension time course (30, 60, 90, 120 and 150 sec) assays using a primed DNA substrate (133 pM) constructed by annealing a radioactively labelled 5′-32P 15mer primer (5′-AAGGGGGTGTGTGTG-3′) onto a 110mer oligonucleotide (5′-CACACACACACACACACACACACACACACACACACACACACACACACACACACACACACCCCTTTAAAAAAAAAAAAAAAAGCCAAAAGCAGTGCCAAGCTTGCATGCC-3′), at suboptimal 25 pM DnaE concentration (left), in the presence of 25 pM PykA tetramer (middle) and 25 pM PEPut domain monomer (right). At this suboptimal DnaE concentration, there is no detectable DnaE primer extension activity in the absence of PykA but clear activity is visible in the presence of PykA. By comparison, no DnaE activity is detectable in the presence of the purified PEPut domain. These data show that full length PykA stimulates the DnaE activity while the PEPut domain alone does not.

Additional file 2: Tables S1-S4


Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Horemans, S., Pitoulias, M., Holland, A. et al. Pyruvate kinase, a metabolic sensor powering glycolysis, drives the metabolic control of DNA replication. BMC Biol 20, 87 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Replication timing
  • Central carbon metabolism
  • Glycolytic enzymes
  • Replication enzymes
  • Cell cycle
  • Allosteric regulation
  • Protein phosphorylation
  • PEPut domain
  • Signaling