- Research article
- Open Access
Pyruvate kinase, a metabolic sensor powering glycolysis, drives the metabolic control of DNA replication
BMC Biology volume 20, Article number: 87 (2022)
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 . 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 . 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 .
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) . 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 .
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 . 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 . 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 . 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 . They are predicted to stabilize the active conformer in a process involving an extra C-terminal domain (see below) . 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 . 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) . 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 . 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  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 . 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.
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.
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) . 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) . 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  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 .
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.
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 . 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.
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 ). 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.
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  and is diffused throughout the cytoplasm . 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 . 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 . 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).
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 , 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 . 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 , 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.
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 , and this interaction was proposed to assist heterotrophic effectors in stabilizing the active R-state conformation . 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].
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 , and (iii) alteration of the metabolic control of replication . 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  while the 208–234 region resides at the heart of the catalytic site and is genetically linked to the lagging-strand polymerase DnaE . 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 , 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 . 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  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, ) 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).
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%).
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  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, ) 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) . 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 (https://www.socscistatistics.com/tests/mannwhitney/default2.aspx) 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 . 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 . 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 . 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 .
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.
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 (http://www.polyomics.gla.ac.uk/assets/downloads/MSMetabolomicsPrepCells-Aug2013.pdf) and modified as previously described . For comparative analysis, pykAT>D and wild-type cells were grown on solid medium and metabolomes were prepared as previously .
Analyses were conducted using a Dionex Ultimate 3000 Rapid Separation LC (Thermo Fisher Scientific). The method for metabolome analysis was previously reported . For hydrogen/deuterium exchange (HDX) experiments, chromatographic conditions were also described .
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 . 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  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 (https://www.ebi.ac.uk/metabolights/MTBLS4415) .
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.
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 (https://www.ebi.ac.uk/metabolights/MTBLS4415) . Strains and plasmids are available upon request from the authors as well as sequence details.
Siddiqui K, On KF, Diffley JFX. Regulating DNA replication in eukarya. Cold Spring Harb Perspect Biol. Cold Spring Harbor Lab. 2013;5:a012930.
Jameson K, Wilkinson AJ. Control of initiation of DNA replication in Bacillus subtilis and Escherichia coli. Genes. 2017;8:22–32.
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.
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.
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.
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.
Buchakjian MR, Kornbluth S. The engine driving the ship: metabolic steering of cell proliferation and death. Nat Rev Microbiol. 2010;11:715–27.
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.
Ewald JC. How yeast coordinates metabolism, growth and division. Curr Opin Microbiol. 2018;45:1–7.
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.
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.
Tu BP, Kudlicki A, Rowicka M, McKnight SL. Logic of the yeast metabolic cycle: temporal compartmentalization of cellular processes. Science. 2005;310:1152–8.
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.
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.
Murray H, Koh A. Multiple regulatory systems coordinate DNA replication with cell growth in Bacillus subtilis. PLoS Genet. 2014;10(10):e1004731.
Mathews CK. Deoxyribonucleotide metabolism, mutagenesis and cancer. Nat Rev Cancer. 2015;15:528–39.
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.
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.
Lu M, Campbell JL, Boye E, Kleckner N. SeqA: a negative modulator of replication nitiation in E. coli. Cell. 1994;77:413–26.
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.
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.
Boye E, Nordström K. Coupling the cell cycle to cell growth. EMBO Rep. 2003;4(8):757–60.
Wang JD, Levin PA. Metabolism, cell growth and the bacterial cell cycle. Nat Rev Microbiol. 2009;7(11):822–7.
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.
Beaufay F, Coppine J, Hallez R. When the metabolism meets the cell cycle in bacteria. Curr Opin Microbiol. 2021;60:104–13.
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.
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.
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.
Chubukov V, Gerosa L, Kochanowski K, Sauer U. Coordination of microbial metabolism. Nat Rev Microbiol. 2014;12:327–40.
Liu GY, Sabatini DM. mTOR at the nexus of nutrition, growth, ageing and disease. Nat Rev Mol Cell Biol. 2021;21:183–203.
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.
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.
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.
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.
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.
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.
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.
Laffan J, Firshein W. Membrane protein binding to the origin region of Bacillus subtilis. J Bacteriol. 1987;169(9):4135–40.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Soultanas P. Loading mechanisms of ring helicases at replication origins. Mol Microbiol. 2012;84(1):6–16.
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.
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.
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.
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.
Sprague GFJ. Isolation and characterization of a Saccharomyces cerevisiae mutant deficient in pyruvate kinase activity. J Bacteriol. 1977;130(1):232–41.
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.
Ronai Z. Glycolytic enzymes as DNA binding proteins. Int J Biochem Cell Biol. 1993;25(7):1073–6.
Sirover MA. New insights into an old protein: The functional diversity of mammalian glyceraldehyde-3-phosphate dehydrogenase. Biochim Biophys Acta. 1999;1432:159–84.
Sirover MA. On the functional diversity of glyceraldehyde-3-phosphate dehydrogenase: Biochemical mechanisms and regulatory control. Biochim Biophys Acta. 2011;1810:741–51.
Kim J-W, Dang CV. Multifaceted roles of glycolytic enzymes. Trends Biochem Sci. 2005;30(3):142–50.
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.
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.
Lu Z, Hunter T. Metabolic kinases moonlighting as protein kinases. Trends Biochem Sci. 2018;43(4):301–10.
Snaebjornsson MT, Schulze A. Non-canonical functions of enzymes facilitate cross-talk between cell metabolic and regulatory pathways. Exp Mol Med. 2018;50:34.
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.
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.
Wellen HG, Sachdeva UM, Bui TV, Cross JR, Thompson CB. ATP-citrate lyase links cellular metabolism to histone acetylation. Science. 2009;324:1076–80.
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.
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.
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.
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.
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.
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.
Jindal HK, Vishwanatha JK. Functional identity of a primer recognition protein as phosphoglycerate kinase. J Biol Chem. 1990;265(12):6540–3.
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.
Suzuki K, Ito S, Shimizu-Ibuka A, Sakai H. Crystal structure of pyruvate kinase from Geobacillus stearothermophilus. J Biochem. 2008;144(3):305–12.
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.
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.
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.
Sakai H. Possible structure and function of the extra C-terminal sequence of pyruvate kinase from Bacillus stearothermophilus. J Biochem. 2004;136:471–6.
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.
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.
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.
Goss NH, Evans CT, Wood HG. Pyruvate phosphate dikinase: sequence of the histidyl peptide, the pyrophosphoryl and phosphoryl carrier. Biochemistry. 1980;19:5805–9.
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.
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.
Burnell JN, Chastain CJ. Cloning and expression of maize-leaf pyruvate, Pi dikinase regulatory protein gene. Biochem Biophys Res Commun. 2006;345:675–80.
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.
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.
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.
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.
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.
Pisithkul T, Patel NM, Amador-Noguez D. Post-translational modifications as key regulators of bacterial metabolic fluxes. Curr Opin Microbiol. 2015;24:29–37.
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.
Bollenbach TJ, Mesecar AD, Nowak T. Role of lysine 240 in the mechanism of yeast pyruvate kinase catalysis. Biochemistry. 1999;38:9137–45.
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.
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.
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.
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.
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.
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.
Bruck I, O'Donnell M. The DNA replication machine of a gram-positive organism. J Biol Chem. 2000;275(37):28971–83.
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.
Pancholi V, Chhatwal GS. Housekeeping enzymes as virulence factors for pathogens. Int J Med Microbiol. 2003;293:391–401.
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.
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.
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.
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.
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.
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.
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.
Rodgers K, McVey M. Error-Prone repair of DNA double-strand breaks. J Cell Physiol. 2016;231:15–24.
Loeb LA. Human cancers express mutator phenotypes: origin, consequences and targeting. Nat Rev Cancer. 2011;11(6):450–7.
Macheret M, Halazonetis TD. DNA replication stress as a hallmark of cancer. Annu Rev Pathol Mech Dis. 2015;10:425–48.
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.
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.
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.
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.
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.
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.
Tautenhahn R, Böttcher C, Neumann S. Highly sensitive feature detection for high resolution LC/MS. BMC Bioinformatics. 2008;9:504.
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.
Lê S, Josse J, Husson F. FactoMineR: An R package for multivariate analysis. J Stat Softw. 2008;25(1):1–18.
Haug. MetaboLights : a resource evolving in response to the needs of its scientific community. Nucleic Acids Res. 2020;48(D1):D440–4.
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.
Ethics approval and consent to participate
Consent for publication
All authors and funders have given their consent for publication.
The authors declare they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. 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.
About this article
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). https://doi.org/10.1186/s12915-022-01278-3
- Replication timing
- Central carbon metabolism
- Glycolytic enzymes
- Replication enzymes
- Cell cycle
- Allosteric regulation
- Protein phosphorylation
- PEPut domain