Insights into the quality of DnaA boxes and their cooperativity

Bacterial Replication Initiation as Precision Control by Protein Counting

Balanced biosynthesis is the hallmark of bacterial cell physiology, where the concentrations of stable proteins remain steady. However, this poses a conceptual challenge to modeling the cell-cycle and cell-size controls in bacteria, as prevailing concentration-based eukaryote models are not directly applicable. In this study, we revisit and significantly extend the initiator-titration model, proposed 30 years ago, and we explain how bacteria precisely and robustly control replication initiation based on the mechanism of protein copy-number sensing. Using a mean-field approach, we first derive an analytical expression of the cell size at initiation based on three biological mechanistic control parameters for an extended initiator-titration model. We also study the stability of our model analytically and show that initiation can become unstable in multifork replication conditions. Using simulations, we further show that the presence of the conversion between active and inactive initiator protein forms significantly represses initiation instability. Importantly, the two-step Poisson process set by the initiator titration step results in significantly improved initiation synchrony with C V ~ 1 / N scaling rather than the standard 1 / N scaling in the Poisson process, where N is the total number of initiators required for initiation. Our results answer two long-standing questions in replication initiation: (i) Why do bacteria produce almost two orders of magnitude more DnaA, the master initiator proteins, than required for initiation? (ii) Why does DnaA exist in active (DnaA-ATP) and inactive (DnaA-ADP) forms if only the active form is competent for initiation? The mechanism presented in this work provides a satisfying general solution to how the cell can achieve precision control without sensing protein concentrations, with broad implications from evolution to the design of synthetic cells.

Replication initiation in bacteria: precision control based on protein counting

Balanced biosynthesis is the hallmark of bacterial cell physiology, where the concentrations of stable proteins remain steady. However, this poses a conceptual challenge to modeling the cell-cycle and cell-size controls in bacteria, as prevailing concentration-based eukaryote models are not directly applicable. In this study, we revisit and significantly extend the initiator-titration model, proposed thirty years ago, and explain how bacteria precisely and robustly control replication initiation based on the mechanism of protein copy-number sensing. Using a mean-field approach, we first derive an analytical expression of the cell size at initiation based on three biological mechanistic control parameters for an extended initiator-titration model. We also study the stability of our model analytically and show that initiation can become unstable in multifork replication conditions. Using simulations, we further show that the presence of the conversion between active and inactive initiator protein forms significantly represses initiation instability. Importantly, the two-step Poisson process set by the initiator titration step results in significantly improved initiation synchrony with C V ~ 1 / N scaling rather than the standard 1 / N scaling in the Poisson process, where N is the total number of initiators required for initiation. Our results answer two long-standing questions in replication initiation: (1) Why do bacteria produce almost two orders of magnitude more DnaA, the master initiator proteins, than required for initiation? (2) Why does DnaA exist in active (DnaA-ATP) and inactive (DnaA-ADP) forms if only the active form is competent for initiation? The mechanism presented in this work provides a satisfying general solution to how the cell can achieve precision control without sensing protein concentrations, with broad implications from evolution to the design of synthetic cells.

Diverse Mechanisms of Helicase Loading during DNA Replication Initiation in Bacteria

Initiation of DNA replication is required for cell viability and passage of genetic information to the next generation. Studies in Escherichia coli and Bacillus subtilis have established ATPases associated with diverse cellular activities (AAA+) as essential proteins required for loading of the replicative helicase at replication origins. AAA+ ATPases DnaC in E. coli and DnaI in B. subtilis have long been considered the paradigm for helicase loading during replication in bacteria. Recently, it has become increasingly clear that most bacteria lack DnaC/DnaI homologs. Instead, most bacteria express a protein homologous to the newly described DciA (dnaC/dnaI antecedent) protein. DciA is not an ATPase, and yet it serves as a helicase operator, providing a function analogous to that of DnaC and DnaI across diverse bacterial species. The recent discovery of DciA and of other alternative mechanisms of helicase loading in bacteria has changed our understanding of DNA replication initiation. In this review, we highlight recent discoveries, detailing what is currently known about the replicative helicase loading process across bacterial species, and we discuss the critical questions that remain to be investigated.

Evidence for a chromosome origin unwinding system broadly conserved in bacteria

Genome replication is a fundamental requirement for the proliferation of all cells. Throughout the domains of life, conserved DNA replication initiation proteins assemble at specific chromosomal loci termed replication origins and direct loading of replicative helicases (1). Despite decades of study on bacterial replication, the diversity of bacterial chromosome origin architecture has confounded the search for molecular mechanisms directing the initiation process. Recently a basal system for opening a bacterial chromosome origin (oriC) was proposed (2). In the model organism Bacillus subtilis, a pair of double-stranded DNA (dsDNA) binding sites (DnaA-boxes) guide the replication initiator DnaA onto adjacent single-stranded DNA (ssDNA) binding motifs (DnaA-trios) where the protein assembles into an oligomer that stretches DNA to promote origin unwinding. We report here that these core elements are predicted to be present in the majority of bacterial chromosome origins. Moreover, we find that the principle activities of the origin unwinding system are conserved in vitro and in vivo. The results suggest that this basal mechanism for oriC unwinding is broadly functionally conserved and therefore may represent an ancestral system to open bacterial chromosome origins.

The DnaA Tale

More than 50 years have passed since the presentation of the Replicon Model which states that a positively acting initiator interacts with a specific site on a circular chromosome molecule to initiate DNA replication. Since then, the origin of chromosome replication, oriC, has been determined as a specific region that carries sequences required for binding of positively acting initiator proteins, DnaA-boxes and DnaA proteins, respectively. In this review we will give a historical overview of significant findings which have led to the very detailed knowledge we now possess about the initiation process in bacteria using Escherichia coli as the model organism, but emphasizing that virtually all bacteria have DnaA proteins that interacts with DnaA boxes to initiate chromosome replication. We will discuss the dnaA gene regulation, the special features of the dnaA gene expression, promoter strength, and translation efficiency, as well as, the DnaA protein, its concentration, its binding to DnaA-boxes, and its binding of ATP or ADP. Furthermore, we will discuss the different models for regulation of initiation which have been proposed over the years, with particular emphasis on the Initiator Titration Model.

Nucleotide-Induced Conformational Changes in Escherichia coli DnaA Protein Are Required for Bacterial ORC to Pre-RC Conversion at the Chromosomal Origin

DnaA oligomerizes when bound to origins of chromosomal replication. Structural analysis of a truncated form of DnaA from Aquifex aeolicus has provided insight into crucial conformational differences within the AAA+ domain that are specific to the ATP- versus ADP- bound form of DnaA. In this study molecular docking of ATP and ADP onto Escherichia coli DnaA, modeled on the crystal structure of Aquifex aeolicus DnaA, reveals changes in the orientation of amino acid residues within or near the vicinity of the nucleotide-binding pocket. Upon limited proteolysis with trypsin or chymotrypsin ADP-DnaA, but not ATP-DnaA generated relatively stable proteolytic fragments of various sizes. Examined sites of limited protease susceptibility that differ between ATP-DnaA and ADP-DnaA largely reside in the amino terminal half of DnaA. The concentration of adenine nucleotide needed to induce conformational changes, as detected by these protease susceptibilities of DnaA, coincides with the conversion of an inactive bacterial origin recognition complex (bORC) to a replication efficient pre-replication complex (pre-RC) at the E. coli chromosomal origin of replication (oriC).

Genome evolution in an ancient bacteria-ant symbiosis: parallel gene loss among Blochmannia spanning the origin of the ant tribe Camponotini

Stable associations between bacterial endosymbionts and insect hosts provide opportunities to explore genome evolution in the context of established mutualisms and assess the roles of selection and genetic drift across host lineages and habitats. Blochmannia, obligate endosymbionts of ants of the tribe Camponotini, have coevolved with their ant hosts for ∼40 MY. To investigate early events in Blochmannia genome evolution across this ant host tribe, we sequenced Blochmannia from two divergent host lineages, Colobopsis obliquus and Polyrhachis turneri, and compared them with four published genomes from Blochmannia of Camponotus sensu stricto. Reconstructed gene content of the last common ancestor (LCA) of these six Blochmannia genomes is reduced (690 protein coding genes), consistent with rapid gene loss soon after establishment of the symbiosis. Differential gene loss among Blochmannia lineages has affected cellular functions and metabolic pathways, including DNA replication and repair, vitamin biosynthesis and membrane proteins. Blochmannia of P. turneri (i.e., B. turneri) encodes an intact DnaA chromosomal replication initiation protein, demonstrating that loss of dnaA was not essential for establishment of the symbiosis. Based on gene content, B. obliquus and B. turneri are unable to provision hosts with riboflavin. Of the six sequenced Blochmannia, B. obliquus is the earliest diverging lineage (i.e., the sister group of other Blochmannia sampled) and encodes the fewest protein-coding genes and the most pseudogenes. We identified 55 genes involved in parallel gene loss, including glutamine synthetase, which may participate in nitrogen recycling. Pathways for biosynthesis of coenzyme A, terpenoids and riboflavin were lost in multiple lineages, suggesting relaxed selection on the pathway after inactivation of one component. Analysis of Illumina read datasets did not detect evidence of plasmids encoding missing functions, nor the presence of coresident symbionts other than Wolbachia. Although gene order is strictly conserved in four Blochmannia of Camponotus sensu stricto, comparisons with deeply divergent lineages revealed inversions in eight genomic regions, indicating ongoing recombination despite ancestral loss of recA. In sum, the addition of two Blochmannia genomes of divergent host lineages enables reconstruction of early events in evolution of this symbiosis and suggests that Blochmannia lineages may experience distinct, host-associated selective pressures. Understanding how evolutionary forces shape genome reduction in this system may help to clarify forces driving gene loss in other bacteria, including intracellular pathogens.

Two oppositely oriented arrays of low-affinity recognition sites in oriC guide progressive binding of DnaA during Escherichia coli pre-RC assembly

The onset of chromosomal DNA replication requires highly precise and reproducible interactions between initiator proteins and replication origins to assemble a pre-replicative complex (pre-RC) that unwinds the DNA duplex. In bacteria, initiator protein DnaA, bound to specific high- and low-affinity recognition sites within the unique oriC locus, comprises the pre-RC, but how complex assembly is choreographed to ensure precise initiation timing during the cell cycle is not well understood. In this study, we present evidence that higher-order DnaA structures are formed at oriC when DnaA monomers are closely positioned on the same face of the DNA helix by interaction with two oppositely oriented essential arrays of closely spaced low-affinity DnaA binding sites. As DnaA levels increase, peripheral high-affinity anchor sites begin cooperative loading of the arrays, which is extended by sequential binding of additional DnaA monomers resulting in growth of the complexes towards the centre of oriC. We suggest that this polarized assembly of unique DnaA oligomers within oriC plays an important role in mediating pre-RC activity and may be a feature found in all bacterial replication origins.

AT-rich region and repeated sequences - the essential elements of replication origins of bacterial replicons

Repeated sequences are commonly present in the sites for DNA replication initiation in bacterial, archaeal, and eukaryotic replicons. Those motifs are usually the binding places for replication initiation proteins or replication regulatory factors. In prokaryotic replication origins, the most abundant repeated sequences are DnaA boxes which are the binding sites for chromosomal replication initiation protein DnaA, iterons which bind plasmid or phage DNA replication initiators, defined motifs for site-specific DNA methylation, and 13-nucleotide-long motifs of a not too well-characterized function, which are present within a specific region of replication origin containing higher than average content of adenine and thymine residues. In this review, we specify methods allowing identification of a replication origin, basing on the localization of an AT-rich region and the arrangement of the origin's structural elements. We describe the regularity of the position and structure of the AT-rich regions in bacterial chromosomes and plasmids. The importance of 13-nucleotide-long repeats present at the AT-rich region, as well as other motifs overlapping them, was pointed out to be essential for DNA replication initiation including origin opening, helicase loading and replication complex assembly. We also summarize the role of AT-rich region repeated sequences for DNA replication regulation.

The Caulobacter crescentus chromosome replication origin evolved two classes of weak DnaA binding sites

The Caulobacter crescentus replication initiator DnaA and essential response regulator CtrA compete to control chromosome replication. The C. crescentus replication origin (Cori) contains five strong CtrA binding sites but only two apparent DnaA boxes, termed G-boxes (with a conserved second position G, TGATCCACA). Since clusters of DnaA boxes typify bacterial replication origins, this discrepancy suggested that C. crescentus DnaA recognizes different DNA sequences or compensates with novel DNA-binding proteins. We searched for novel DNA sites by scanning mutagenesis of the most conserved Cori DNA. Autonomous replication assays showed that G-boxes and novel W-boxes (TCCCCA) are essential for replication. Further analyses showed that C. crescentus DnaA binds G-boxes with moderate and W-boxes with very weak affinities significantly below DnaA's capacity for high-affinity Escherichia coli-boxes (TTATCCACA). Cori has five conserved W-boxes. Increasing W-box affinities increases or decreases autonomous replication depending on their strategic positions between the G-boxes. In vitro, CtrA binding displaces DnaA from proximal G-boxes and from distal W-boxes implying CtrA-DnaA competition and DnaA-DnaA cooperation between G-boxes and W-boxes. Similarly, during cell cycle progression, CtrA proteolysis coincides with DnaA binding to Cori. We also observe highly conserved W-boxes in other replication origins lacking E. coli-boxes. Therefore, strategically weak DnaA binding can be a general means of replication control.

Remodeling of nucleoprotein complexes is independent of the nucleotide state of a mutant AAA+ protein

DnaA protein, a member of the AAA+ (ATPase associated with various cellular activities) family, initiates DNA synthesis at the chromosomal origin of replication (oriC) and regulates the transcription of several genes, including its own. The assembly of DnaA complexes at chromosomal recognition sequences is affected by the tight binding of ATP or ADP by DnaA. DnaA with a point mutation in its membrane-binding amphipathic helix, DnaA(L366K), previously described for its ability to support growth in cells with altered phospholipid content, has biochemical characteristics similar to those of the wild-type protein. Yet DnaA(L366K) fails to initiate in vitro or in vivo replication from oriC. We found here, through in vitro dimethyl sulfate footprinting and gel mobility shift assays, that DnaA(L366K) in either nucleotide state was unable to assemble into productive prereplication complexes. In contrast, at the dnaA promoter, both the ATP and the ADP form of DnaA(L366K) generated active nucleoprotein complexes that efficiently repressed transcription in a manner similar to wild-type ATP-DnaA. Thus, it appears that unlike wild-type DnaA protein DnaA(L366K) can adopt architectures that are independent of its bound nucleotide, and instead the locus determines the functionality of the higher order DnaA(L366K)-DNA complexes.

DnaA structure, function, and dynamics in the initiation at the chromosomal origin

Escherichia coli DnaA is the initiator of chromosomal replication. Multiple ATP-DnaA molecules assemble at the oriC replication origin in a highly regulated manner, and the resultant initiation complexes promote local duplex unwinding within oriC, resulting in open complexes. DnaB helicase is loaded onto the unwound single-stranded region within oriC via interaction with the DnaA multimers. The tertiary structure of the functional domains of DnaA has been determined and several crucial residues in the initiation process, as well as their unique functions, have been identified. These include specific DNA binding, inter-DnaA interaction, specific and regulatory interactions with ATP and with the unwound single-stranded oriC DNA, and functional interaction with DnaB helicase. An overall structure of the initiation complex is also proposed. These are important for deepening our understanding of the molecular mechanisms that underlie DnaA assembly, oriC duplex unwinding, regulation of the initiation reaction, and DnaB helicase loading. In this review, we summarize recent progress on the molecular mechanisms of the functions of DnaA on oriC. In addition, some members of the AAA+ protein family related to the initiation of replication and its regulation (e.g., DnaA) are briefly discussed.

Distribution of stable DnaA-binding sites on the Bacillus subtilis genome detected using a modified ChIP-chip method

We developed a modified ChIP-chip method, designated ChAP-chip (Chromatin Affinity Precipitation coupled with tiling chip). The binding sites of Bacillus subtilis Spo0J determined using this technique were consistent with previous findings. A DNA replication initiator protein, DnaA, formed stable complexes at eight intergenic regions on the B. subtilis genome. Characterization of the binding sequences suggested that two factors—the local density of DnaA boxes and their affinities for DnaA—are critical for stable binding. We further showed that in addition to autoregulation, DnaA directly modulate the expression of sda in a positive, and ywlC and yydA in a negative manner. Examination of possible stable DnaA-binding sequences in other Bacillus species suggested that DnaA-dependent regulation of those genes is maintained in most bacteria examined, supporting their biological significance. In addition, a possible stable DnaA-binding site downstream of gcp is also suggested to be conserved. Furthermore, potential DnaA-binding sequences specific for each bacterium have been identified, generally in close proximity to oriC. These findings suggest that DnaA plays several additional roles, such as control of the level of effective initiator, ATP-DnaA, and/or stabilization of the domain structure of the genome around oriC for the proper initiation of chromosome replication.

Mechanism for the TtDnaA-Tt-oriC cooperative interaction at high temperature and duplex opening at an unusual AT-rich region in Thermoanaerobacter tengcongensis

Thermoanaerobacter tengcongensis is an anaerobic low-GC thermophilic bacterium. To further elucidate the replication initiation of chromosomal DNA at high temperature, the interaction between the replication initiator (TtDnaA) and the putative origin (Tt-oriC) in this thermophile was investigated. We found that efficient binding of TtDnaA to Tt-oriC at high temperature requires (i) at least two neighboring DnaA boxes, (ii) the specific feature of the TtDnaA Domain IV and (iii) the self-oligomerization of TtDnaA. Replacement of the TtDnaA Domain IV by the counterpart of Escherichia coli DnaA or disruption of its oligomerization by amino acid mutations (W9A/L20S) abolished the oriC-binding activity of TtDnaA at 60 degrees C, but not at 37 degrees C. Moreover, ATP-TtDnaA, but not ADP-TtDnaA or the oligomerization-deficient mutants was able to unwind the Tt-oriC duplex. The minimal oriC required for this duplex opening in vitro was demonstrated to consist of DnaA boxes 1-8 and an unusual AT-rich region. Interestingly, although no typical ATP-DnaA box was found in this AT-rich region, it was exclusively bound by ATP-TtDnaA and acted as the duplex-opening and replication-initiation site. Taken together, we propose that oligomerization of ATP-DnaA and simultaneously binding of several DnaA boxes and/or AT-rich region may be generally required in replication initiation at high temperature.

Nucleotide sequence analysis of cryptic plasmid pAM5 from Acidiphilium multivorum

Plasmid pAM5 of Acidiphilium multivorum JCM-8867 has been completely sequenced by initial cloning of HindIII-PstI fragments followed by primer walking. It has a size of 5161bp and single site for several restriction enzymes as revealed by DNA sequencing. Sequence analysis predicts five putative open reading frames. ORF1 and ORF3 show significant identity with various plasmid encoded mobilization (Mob) and replication initiation (Rep) proteins, respectively. The putative Mob protein has several characteristics of the MOB(Q) family having the motifs with conserved amino acid residues. Upstream of the Mob ORF, there exists a 34bp oriT region having a nic consensus sequence. The constructed plasmid pSK1 bearing pAM5 mob region can be mobilized to Escherichia coli in presence of conjugative plasmid pRK2013. The replication module comprises of several DnaA like boxes, several perfect direct and inverted repeats, a potential prokaryotic promoter and putative rep gene. The rep module is very similar to several theta replicating iteron family plasmids, suggesting pAM5 replication to follow the same course. Any phenotypic character determinant (e.g., metal resistance, antibiotic resistance etc.) gene is absent in pAM5, suggesting this plasmid to be cryptic in nature. However, a pAM5 derivative plasmid named pSK2, containing the putative pAM5 rep region, can replicate and be stably maintained in Acidiphilium, Acidocella, and E. coli strains; it can also carry foreign DNA fragments. Thus, pSK2 could serve as a cloning shuttle vector between these bacteria. It was observed that pAM5 Rep is essential for pSK2 to replicate in acidophiles. In its natural host, A. multivorum JCM-8867, pAM5 maintains a copy number of 50-60, and its derivative pSK2 maintains a comparatively, higher copy number in E. coli than in acidophiles.

Mechanism of origin activation by monomers of R6K-encoded pi protein

One recurring theme in plasmid duplication is the recognition of the origin of replication (ori) by specific Rep proteins that bind to DNA sequences called iterons. For plasmid R6K, this process involves a complex interplay between monomers and dimers of the Rep protein, pi, with seven tandem iterons of gamma ori. Remarkably, both pi monomers and pi dimers can bind to iterons, a new paradigm in replication control. Dimers, the predominant form in the cell, inhibit replication, while monomers facilitate open complex formation and activate the ori. Here, we investigate a mechanism by which pi monomers out-compete pi dimers for iteron binding, and in so doing activate the ori. With an in vivo plasmid incompatibility assay, we find that pi monomers bind cooperatively to two adjacent iterons. Cooperative binding is eliminated by insertion of a half-helical turn between two iterons but is diminished only slightly by insertion of a full helical turn between two iterons. These studies show also that pi bound to a consensus site promotes occupancy of an adjacent mutated site, another hallmark of cooperative interactions. pi monomer/iteron interactions were quantified using a monomer-biased pi variant in vitro with the same collection of two-iteron constructs. The cooperativity coefficients mirror the plasmid incompatibility results for each construct tested. pi dimer/iteron interactions were quantified with a dimer-biased mutant in vitro and it was found that pi dimers bind with negligible cooperativity to two tandem iterons.

Sequence characteristics required for cooperative binding and efficient in vivo titration of the replication initiator protein DnaA in E. coli

Plasmids carrying the mioC promoter region, which contains two DnaA boxes, R5 and R6 with one misfit to the consensus TT(A)/(T)TNCACA, are as efficient in in vivo titration of the DnaA protein as plasmids carrying a replication-inactivated oriC region with its eight DnaA boxes. Three additional DnaA boxes around the promoter proximal R5 DnaA box were identified and shown by mutational analysis to be necessary for the cooperative binding of DnaA required for titration. These four DnaA boxes are located in the same orientation and with a spacing of two or three base-pairs. The cooperative binding was eliminated by insertion of half a helical turn between any of the DnaA boxes. Titration strongly depends on the presence and orientation of the promoter distal R6 DnaA box located 104 bp upstream of the R5 box as well as neighbouring sequences downstream of R6. Titration depends on the integrity of a 43 bp region containing the R5 DnaA box, while repression of mioC transcription by DnaA, which is dependent on the R5 DnaA box, was independent of the two DnaA boxes downstream of R5. Repression was also independent of the spacing between the two upstream DnaA boxes and the promoter as long as a DnaA box was located less than 20 bp upstream of the -35 sequence. Thus, the architectural requirements for titration and for repression of transcription are different. A new set of rules for identifying efficiently titrating DnaA box regions was formulated and used to analyse sequences for which good titration data are available.

Pre‐replication assembly ofE. colireplisome components

The localization of SeqA, thymidylate synthase, DnaB (helicase) and the DNA polymerase components alpha and tau, has been studied by immunofluorescence microscopy. The origin has been labelled through GFP-LacI bound near oriC. SeqA was located in the cell centre for one replication factory (RF) and at 1/4 and 3/4 positions in pre-divisional cells harbouring two RFs. The transition of central to 1/4 and 3/4 positions of SeqA appeared abrupt. Labelled thymidylate synthetase was found all over the cell, thus not supporting the notion of a dNTP-synthesizing complex exclusively localized near the RF. More DnaB, alpha and tau foci were found than expected. We have hypothesized that extra foci arise at pre-replication assembly sites, where the number of sites equals the number of origins, i.e. the number of future RFs. A reasonable agreement was found between predicted and found foci. In the case of multifork replication the number of foci appeared consistent with the assumption that three RFs are grouped into a higher-order structure. The RF is probably separate from the foci containing SeqA and the hemi-methylated SeqA binding sites because these foci did not coincide significantly with DnaB as marker of the RF. Co-labelling of DnaB and oriC revealed limited colocalization, indicating that DnaB did not yet become associated with oriC at a pre-replication assembly site. DnaB and tau co-labelled in the cell centre, though not at presumed pre-replication assembly sites. By contrast, alpha and tau co-labelled consistently suggesting that they are already associated before replication starts.