Article
A Bacterial Chromosome Structuring Protein Binds Overtwisted DNA to Stimulate Type II Topoisomerases and Enable DNA Replication Graphical Abstract
Authors rapid fork progression
overtwist / (+) supercoiled DNA
DNA replication
Monica S. Guo, Diane L. Haakonsen, Wenjie Zeng, Maria A. Schumacher, Michael T. Laub
+
generates
+ +
Correspondence
[email protected] (M.A.S.),
[email protected] (M.T.L.)
In Brief topoisomerases remove (+) supercoils
GapR clamps overtwisted DNA GapR stimulates topoisomerases
topoisomerase activity
GapR
Highlights d
GapR is a DNA structuring protein conserved throughout a-proteobacteria
d
GapR promotes DNA replication initiation and elongation
d
GapR encircles overtwisted DNA as a dimer-of-dimers
d
GapR stimulates the relaxation of positive supercoils by gyrase and topo IV
Guo et al., 2018, Cell 175, 583–597 October 4, 2018 ª 2018 Elsevier Inc. https://doi.org/10.1016/j.cell.2018.08.029
A bacterial protein that specifically recognizes and encircles overtwisted DNA is required to stimulate the activity of type II topoisomerases and enable DNA replication.
Article A Bacterial Chromosome Structuring Protein Binds Overtwisted DNA to Stimulate Type II Topoisomerases and Enable DNA Replication Monica S. Guo,1,4 Diane L. Haakonsen,1,4,5 Wenjie Zeng,3 Maria A. Schumacher,3,* and Michael T. Laub1,2,6,* 1Department
of Biology, Massachusetts Institute of Technology, Cambridge, MA 02139, USA Hughes Medical Institute, Massachusetts Institute of Technology, Cambridge, MA 02139, USA 3Department of Biochemistry, Duke University School of Medicine, Durham, NC 27710, USA 4These authors contributed equally 5Present address: Department of Molecular and Cell Biology, University of California, Berkeley, Berkeley, CA 94720, USA 6Lead Contact *Correspondence:
[email protected] (M.A.S.),
[email protected] (M.T.L.) https://doi.org/10.1016/j.cell.2018.08.029 2Howard
SUMMARY
When DNA is unwound during replication, it becomes overtwisted and forms positive supercoils in front of the translocating DNA polymerase. Unless removed or dissipated, this superhelical tension can impede replication elongation. Topoisomerases, including gyrase and topoisomerase IV in bacteria, are required to relax positive supercoils ahead of DNA polymerase but may not be sufficient for replication. Here, we find that GapR, a chromosome structuring protein in Caulobacter crescentus, is required to complete DNA replication. GapR associates in vivo with positively supercoiled chromosomal DNA, and our biochemical and structural studies demonstrate that GapR forms a dimer-of-dimers that fully encircles overtwisted DNA. Further, we show that GapR stimulates gyrase and topo IV to relax positive supercoils, thereby enabling DNA replication. Analogous chromosome structuring proteins that locate to the overtwisted DNA in front of replication forks may be present in other organisms, similarly helping to recruit and stimulate topoisomerases during DNA replication. INTRODUCTION In all organisms, chromosomes must be compacted 1,000-fold to fit within the confines of a cell. However, DNA cannot be haphazardly packed and instead must be organized to facilitate a range of cellular processes, including DNA replication, transcription, and chromosome segregation. Chromosome compaction is achieved in large part by the supercoiling of DNA, forming a series of plectonemic loops in bacteria (Higgins, 2016). DNA can be either positively (+) or negatively () supercoiled, depending on whether the DNA duplex winds around a superhelical axis in a left- or right-handed fashion, respectively. The supercoiling of DNA is tightly controlled by topoisomerases,
enzymes that can break and rejoin DNA segments to introduce or relieve superhelical strain (Vos et al., 2011). Despite major advances in our understanding of topoisomerases, the molecular mechanisms that regulate when, where, and how these enzymes act in vivo remain poorly understood. Supercoils are also introduced by cellular processes that unwind the double helix, such as transcription and DNA replication. The ‘‘twin supercoiled domain’’ model posits that, as RNA or DNA polymerase translocates along and unwinds the DNA, it produces overtwisted and (+) supercoiled DNA ahead of it, with underwound, () supercoiled DNA arising in its wake (Liu and Wang, 1987; Wu et al., 1988). These supercoils must be dissipated by topoisomerases to prevent buildup of torsional stress that would inhibit unwinding of the duplex and prevent translocation of RNA or DNA polymerase along the DNA (Postow et al., 2001). The action of topoisomerases is especially critical during DNA replication as 50–100 (+) supercoils per second are generated ahead of the translocating replisome (Postow et al., 2001). Despite the importance of topoisomerases to replication, how the activities of these enzymes are controlled to promote replication elongation remains largely unknown. Topoisomerases fall into two classes, type I and type II, which respectively cleave either one or two strands of DNA (Vos et al., 2011). Bacteria, which maintain their chromosomes with net-negative supercoil, typically have a type I topoisomerase, topo I, that can relieve excess () supercoils, and two type II topoisomerases, DNA gyrase and topo IV, that can relax (+) supercoils. Gyrase also introduces () supercoils into relaxed DNA, while topo IV helps decatenate replicated chromosomes. Because gyrase, and to a lesser extent topo IV, are found in nearly all bacteria but are absent from eukaryotes, these enzymes are targeted by a number of clinically relevant antibiotics, such as the coumerins and quinolones (Vos et al., 2011). Bacteria treated with such drugs cannot replicate their DNA, underscoring the importance of topoisomerases to DNA replication (Khodursky et al., 2000). Although gyrase, perhaps with topo IV, is often assumed to be sufficient to manage the (+) supercoils that arise during replication (Khodursky et al., 2000), a role for other factors cannot be ruled out. In particular, bacteria, which do not encode histones, Cell 175, 583–597, October 4, 2018 ª 2018 Elsevier Inc. 583
harbor a bevy of so-called nucleoid-associated proteins (NAPs) that may influence topoisomerases. NAPs are a functionally heterogeneous set of proteins that are not conserved at the sequence level, but share many similarities such as small size, high basic amino-acid content, and DNA-binding ability (Badrinarayanan et al., 2015). NAPs typically bind DNA sequence non-specifically by recognizing DNA shape, e.g., narrow minor grooves, and many structure the chromosome by bending, wrapping, or bridging DNA, activities that can introduce or constrain supercoils. As NAPs are often highly abundant, their loss can lead to global changes in supercoiling (Lal et al., 2016; Weitao et al., 2000). Several NAPs in Escherichia coli, including SeqA and the larger chromosome structuring protein MukB, can bind to and stimulate topo IV to relax supercoils or decatenate interlinked circular DNA (Hayama and Marians, 2010; Kang et al., 2003; Li et al., 2010). However, to the best of our knowledge, no NAPs have been shown to stimulate gyrase to relax (+) supercoils, and a direct role for NAPs in replication elongation has not been established. Here, we demonstrate that a highly conserved NAP from the a-proteobacterium Caulobacter crescentus called GapR is essential for DNA replication in fast growth conditions by promoting the ability of gyrase and topo IV to relieve the (+) supercoils that accumulate ahead of the replication fork. GapR was identified as an essential DNA-binding protein that binds ATrich DNA and has pleiotropic effects when depleted (Arias-Cartin et al., 2017; Ricci et al., 2016; Taylor et al., 2017), but its molecular function and precise roles in the bacterial cell cycle have not been elucidated. Our studies show that a loss of GapR leads to global changes in superhelicity and an inability to complete DNA replication. Intriguingly, we find that, in addition to AT-rich DNA, GapR globally associates with the 30 ends of highly expressed genes, but only when they are actively transcribed. These results suggest that GapR recognizes and helps dissipate the (+) supercoils that arise during unwinding of the DNA duplex. Consistent with this model, we demonstrate that GapR preferentially binds overtwisted DNA and can stimulate gyrase and topo IV to remove (+) supercoils. Crystal structures of GapR reveal a unique fold with dimeric GapR that transitions, upon DNA binding, to form a dimer-of-dimers that completely encircles DNA. Strikingly, the structure shows that the GapR clamp specifically recognizes the shape of overtwisted DNA, and cannot accommodate B-form DNA, supporting a model in which GapR localizes to the overtwisted DNA ahead of replication forks where it can stimulate type II topoisomerases. RESULTS Identification of the Nucleoid-Associated Protein, GapR Despite their importance in chromosome organization, most NAPs are not broadly conserved. To identify new NAPs in Caulobacter, we isolated intact nucleoids using sucrose gradient centrifugation and performed mass spectrometry on the proteins that remained tightly associated (Figure S1A). After sedimentation, fractions were DAPI stained to identify the nucleoid fraction and the slower migrating cytoplasmic fraction, which contains RNA and highly sheared DNA (Figure S1A). We verified by immunoblotting that a known NAP, integration host factor (IHF), was
584 Cell 175, 583–597, October 4, 2018
stably associated with the nucleoid, while DNA binding proteins with higher dissociation rates, such as the transcription factor CtrA, were found in the cytoplasmic fraction (Figure 1A). Mass spectrometry identified multiple, independent peptides of known NAPs, including HU and IHF, and other tight DNA-binding proteins, including subunits of RNA polymerase (RNAP) and single-stranded DNA binding protein; no known canonical transcription factors were recovered (Table S1). To identify NAPs, we examined our mass spectrometry data for proteins with no annotated function but with properties common to many NAPs: (1) small (<25 kDa), (2) a basic region to promote tight DNA binding, (3) high expression during exponential growth, and (4) conservation across a-proteobacteria. One candidate was an 89 amino-acid protein encoded by the highly conserved gene gapR (CCNA_03428, Figures 1B and S1B). We confirmed that GapR was nucleoid-associated by repeating the nucleoid purification on cells expressing GapR-3xFLAG from its native locus. Epitope-tagged GapR was highly enriched in the nucleoid fraction and not detectable in the cytoplasmic fraction, in contrast to the a-subunit of RNAP that was present in both fractions (Figure 1A), suggesting that GapR binds the nucleoid stably and with high affinity. GapR was also identified in searches for essential Caulobacter DNA-binding proteins (Arias-Cartin et al., 2017; Ricci et al., 2016) and origin-binding proteins (Taylor et al., 2017). These studies suggested that GapR is an essential and highly expressed protein that has pleiotropic effects when depleted, but the molecular function of GapR was not determined. Depleting GapR Leads to Cellular Filamentation and DNA Replication Defects To understand the essential function of GapR, we constructed a depletion strain in which gapR is driven by a glucose-repressible, xylose-inducible promoter on a low-copy plasmid (DgapR PxylgapR). Even in the presence of xylose, this strain has reduced levels (15-fold) of GapR relative to the wild-type. The strain remains viable in xylose, but has a significant growth defect (Figures 1C, 1D, S1C, and S1D). Further inhibiting gapR expression by shifting cells to a medium with glucose led to undetectable levels of GapR protein and a five-log defect in plating viability. The GapR depletion strain was not able to grow in glucose and exhibited cell-division defects, becoming highly filamentous, often with multiple, incomplete constriction sites after a shift to a rich medium containing glucose (Figure 1E). Flow cytometry analysis of chromosome content indicated that most GapRdepleted cells did not accumulate chromosomes despite becoming filamentous (Figure 1E). Instead, the culture was enriched for G1 cells suggesting that DNA replication is disrupted at the level of initiation and/or elongation. We used RNA sequencing (RNA-seq) to analyze changes in gene expression after 6 hr of GapR depletion in glucose relative to wild-type cells grown in glucose. We identified 1,245 genes that changed expression at least 2-fold (32% of all genes; 996 genes upregulated and 249 genes downregulated; Figure 1F) indicating that GapR has broad, pervasive effects on gene expression (Arias-Cartin et al., 2017; Ricci et al., 2016), including a general stress response, indicated by a correlation with the expression changes seen in cells treated with ethanol (Pearson’s
A
B
chromosome
D
cytoplasm
sheared
intact
WT
∆gapR Pxyl-gapR
fraction # 1
2
3
4
5
6
7
8
+xyl
9 10 11
+glu
+xyl
+glu
*
IHF
C
CtrA
WT
∆gapR Pxyl-gapR
GapR
+ xyl GapR-3xFLAG
+glu RpoA
E
F
∆gapR Pxyl-gapR
WT +xyl
+glu 6 hr
+glu 12 hr
# cells
GapR depletion + glu 6 hr log2(RPKM)
+xyl
15
12
9
6
3
0 0 1N
2N
1N
2N
1N
2N
1N
G
2N
1
2
9
12
15
H
recA 0
6
WT + glu log2(RPKM)
log2(expression)
LexA regulated genes
3
supercoiling-sensitive genes
3
r
p-value
∆gapR Pxyl-gapR +xyl
∆gapR Pxyl-gapR +xyl
0.53 10-21
∆gapR Pxyl-gapR +glu
∆gapR Pxyl-gapR +glu
0.61 10-30
+ novobiocin
+ novobiocin
1.00
--
DNA damage
DNA damage
0.29
10-5
SocB induced
SocB induced
0.33
10-8
+ ethanol
+ ethanol
0.17
10-3
SciP depletion
SciP depletion
~0.0
101
-2
log2(expression)
2
-2
log2(expression)
2
Figure 1. Identification of GapR and Phenotypic Characterization of GapR-Depleted Cells (A) Immunoblots of fractions from a nucleoid isolation experiment with antibodies for IHF, CtrA, FLAG, and RpoA. (B) Schematic of predicted GapR domains. (C) Growth of the DgapR Pxyl-gapR depletion strain assessed by 10-fold serial dilutions on plates containing xylose (gapR-induced) or glucose (gapR-repressed). (D) Immunoblots of GapR in wild-type and depletion strains grown in xylose or in glucose for 6 hr. *Non-specific band used as a loading control. (E) Phase-contrast images and flow cytometry analysis of chromosome content of wild-type and DgapR Pxyl-gapR cells in xylose and at the indicated times in glucose. Arrows, sites of incomplete constriction. (F) Gene expression in GapR-depleted (+glu 6 hr) versus wild-type cells in glucose (+glu). Dashes indicate 2-fold changes in expression. (G) GapR-depleted cells induce a DNA damage response. Expression changes of LexA regulated genes (left), including recA (right) in GapR-reduced (+xyl) and GapR-depleted (+glu) cells compared to cells treated with novobiocin, DNA damaging agents, the toxin SocB, ethanol stress, or cell-cycle arrested by SciP depletion (left). (H) Supercoiling-sensitive genes change in GapR-depleted cells. As in (G), except for genes affected >4-fold by novobiocin treatment. Pearson correlation coefficient (r) and p value for comparison to novobiocin treatment. See also Figure S1 and Tables S1 and S2.
r = 0.72, p = 1.2 3 1067; Figure S1E). We observed similar, but more modest, changes when GapR levels were reduced in xylose (Figure S1E; Table S2). The GapR-depletion strain in
glucose also exhibited a DNA damage response (Arias-Cartin et al., 2017), with expression changes similar to cells treated with DNA damaging agents or expressing the toxin SocB, which
Cell 175, 583–597, October 4, 2018 585
C
4
bits
GapR-3xFLAG ChIP-seq
rpm
3
1 0
2
2 bits
1 0
control ChIP-seq
1 0
0.5
1.0
1.5
2.0
2.5
1 0
3.0
3.5
4.0
A ATA TAA AT
T
C A
1 2 3 4 5 6
2
ATAT
5 4
GapR ChIP (rpm)
B
5
3 2 1 0
1 2 3 4 5
A
0
10
20
30
chromosome position (Mbp)
50
60
70
80
E GapR ChIP log2(rpm) at ORF
% AT content
D
rpm
GapR GapR + rif
2
0
-2
-4
-6
chromosome position (bp)
GapR + rif 2175000
2176000
2177000
2178000
2179000
2180000
-2
0
2
4
6
8
2181000
3' ends of transcription units
5' ends of transcription units
change in GapR occupancy (no rif) - (+ rif)
change in GapR occupancy (no rif) - (+ rif)
normalized count
RNAP
ChIP (rpm)
GapR
ChIP (rpm)
% AT content
G
RNAP
-4
log2(GapR-depleted fold change)
F
GapR
40 % AT
2182000
chromosome position (bp)
low expression (< 50 RPKM) transcription units high expression (> 150 RPKM) transcription units
Figure 2. GapR Binds AT-Rich Genomic Regions and the 30 Ends of Highly Expressed Genes (A) ChIP-seq profiles of GapR-3xFLAG (orange) and wild-type cells (black) using an anti-FLAG antibody. (B) Motifs enriched in GapR-bound regions. (C) GapR ChIP-seq versus AT content across the genome calculated using a 100-bp window centered at each position. Mean enrichment at a given % AT (red dots). (D) AT content (top) and ChIP-seq (bottom) of untreated (orange) and rifampicin-treated (pink) GapR-3xFLAG cells. AT content below the genomic average (32%) is plotted in reverse. (E) Average GapR ChIP-seq at a given ORF versus change in gene expression in GapR-depleted cells. Shaded area indicates <2-fold change upon GapR depletion. (F) GapR accumulates at the end of highly expressed operons. GapR-3xFLAG (orange), GapR-3xFLAG post rifampicin treatment (pink), and RNAP (black) ChIPseq profiles are shown for a region containing two highly expressed tRNAs. AT content (top) and positions of annotated genes (bottom, tRNA genes in solid black). (G) GapR accumulates at the 30 end of highly expressed operons. Normalized change in GapR ChIP signal was computed from the difference in enrichment (untreated minus rifampicin-treated cells in a 1-kb window) before and after the 30 ends (left) or 50 ends (right) of transcription units. See also Figure S2.
leads to replication fork collapse (Pearson’s r = 0.57, p = 1.6 3 106, Figure 1G). A DNA damage response was not observed when cells were challenged by ethanol or subjected to cell-cycle disruption by depleting SciP. GapR Binding Is Not Correlated with GapR-Dependent Changes in Gene Expression To gain more insight into GapR function, we analyzed its global DNA binding profile by chromatin immunoprecipitation sequencing (ChIP-seq) using a strain expressing GapR-3xFLAG as the sole copy of GapR from its native locus (Figure S1F). This profile revealed widespread GapR-3xFLAG binding across the
586 Cell 175, 583–597, October 4, 2018
genome (Figure 2A). Using the 34 highest peaks, we identified two highly enriched motifs, each consisting almost exclusively of A and T, whereas the Caulobacter genome overall is 68% GC, suggesting that GapR preferentially binds AT-rich DNA (Figures 2B–2D). These data agree with previous studies indicating that AT-content of the DNA is one determinant of GapR binding (Arias-Cartin et al., 2017; Ricci et al., 2016). However, there are important exceptions to this trend not previously noted that are discussed below. To assess whether GapR directly affects gene expression, we computed the average GapR enrichment within each open reading frame (Figure 2E) or at each promoter (Figure S2A) in
the genome and plotted these values against the expression changes in GapR-depleted cells. There was no correlation, suggesting that GapR does not directly affect gene expression. GapR May Associate with (+) Supercoils at the 30 End of Highly Expressed Genes and Operons How, then, does GapR affect global patterns of transcription? Because some NAPs influence the supercoiling state of cells and consequently alter the expression of supercoiling-sensitive genes, we compared the expression profile of cells depleted of GapR in glucose to that of cells treated with novobiocin, an inhibitor of gyrase and, to a reduced extent, topo IV (Table S2). For the 289 genes that changed at least 4-fold in novobiocin-treated cells, there was a strong correlation with the changes observed in GapR-depleted cells (Pearson’s r = 0.61, p = 5.0 3 1031), although the magnitudes of the changes due to GapR-depletion were typically smaller (Figure 1H). The expression profiles of cells expressing SocB toxin, treated with ethanol, or cell-cycle arrested by SciP depletion, did not similarly correlate with the novobiocin response, suggesting that the changes in supercoiling-sensitive genes following GapR depletion are not just a result of general stress or a disruption of the cell cycle (Figure 1H). The change in supercoiling-sensitive genes suggested to us that GapR-depleted cells may have reduced () supercoiling. Consistent with this conclusion, we found that the net () superhelicity of chromosomes from GapR-depleted cells was reduced relative to that of wild-type cells, but higher than that of novobiocin-treated cells (Figure S1G). Interestingly, we noticed that GapR accumulated in genomic regions that were not AT-rich, but instead corresponded to the 30 ends of highly expressed transcription units such as tRNAs (Figures 2F and S2B). At these locations, the binding of GapR was broad (spanning several kb) and lower in magnitude compared to the AT-rich locations noted above (Figures 2B– 2D). GapR binding in these regions was dependent on active transcription, as it was reduced by the addition of the transcriptional inhibitor rifampicin (Figures 2F and S2B). This transcription-dependent accumulation of GapR was strongly dependent on the expression level of the transcription unit (two-sided t test: p = 1.8 3 1033) and was specific to the 30 ends of transcription units (Figures 2G, S2D, and S2E). In contrast, the binding of GapR at AT-rich genomic regions was not affected by rifampicin (Figures 2D and S2C). Although it is possible that GapR helps terminate transcription, like E. coli H-NS (Ray-Soni et al., 2016), our RNA-seq analysis of GapR-depleted cells did not reveal any transcriptional read-through (Figure S2B). Why would GapR associate with the 30 ends of highly expressed operons? Unwinding of the DNA duplex during transcription leads to the overtwisting of DNA and the accumulation of (+) supercoils ahead of RNAP. Thus, GapR may associate preferentially with the DNA structures that arise at the ends of highly expressed transcription units. Further, we hypothesized that GapR may help eliminate or dissipate (+) supercoils as GapR-depleted cells exhibited decreased () superhelicity (Figures 1H and S1G). If GapR indeed helps relieve (+) supercoil stress, then depleting GapR should disrupt DNA replication, as the unwinding of DNA during replication leads to a significant accumulation of (+) supercoils ahead of replication forks.
GapR Is Required for Proper DNA Replication Initiation and Elongation To assess the role of GapR in DNA replication, we used flow cytometry to measure DNA content in synchronized, GapRdepleted cells over time (Figure 3A). The GapR depletion strain was grown in glucose for 2 hr to deplete GapR, synchronized to isolate cells in G1, and released into fresh media with glucose to continue repressing gapR expression. After release, DNA content was monitored by flow cytometry (Figure 3B). For the wild-type, nearly all cells initiated replication by 40 min post-synchronization, with cell division occurring within 100 min, indicating that replication of the entire chromosome took <60 min (Figures 3B and 3C). In contrast, GapR-depleted cells showed significant defects in both initiation and elongation. Only 20% of cells had initiated within 40 min postsynchronization and only 60% had initiated after 120 min (Figure 3C). The rate of replication elongation was also slowed significantly (Figure 3D), with few cells achieving 2N content (Figure 3B). We also examined the GapR depletion strain under xylose-inducing conditions, which showed more modest but significant defects in replication initiation and elongation rate (Figures S3A–S3D). We then used DNA sequencing to more precisely assess the replication defects of synchronized wild-type and GapRdepleted cells. Replicated regions of the chromosome are easily detected and tracked over time as they have 2-fold more reads than unreplicated regions (Figure 3E). Using this assay, we found that nearly all wild-type cells (90%) initiated replication 20 min post-synchronization (Figure 3F), with replication progressing steadily toward the terminus at a rate of 32–34 kb/min (Figure 3G), consistent with prior estimates (Jensen, 2006). In contrast, for GapR-depleted cells, <20% of cells had initiated after 20 min, and only 35% had initiated by 60 min post-release (Figure 3F). Additionally, the replication forks of GapR-depleted cells showed a progressive slowdown over time, with a rate of only 12 kb/min between 60 and 80 min, and with no cells having finished replication by the 80 min time point when all wild-type cells had completed. We also performed this experiment on the GapR depletion strain shifted to glucose only upon synchronization rather than for 2 hr before synchronization, again observing delays in replication initiation and a progressive decrease in elongation rate (Figures 3F, 3G, and S3E). Together, our flow cytometry and DNA-sequencing data suggested that GapR-depleted cells are not competent to finish replication. We directly tested this model by first monitoring replication by flow cytometry over an extended (4 hr) period. Whereas wild-type cells replicated twice during this period, GapR-depleted cells accumulated intermediate (<2N) chromosome content (Figure S4A). Second, if any GapR-depleted cells finish replication, they should constrict and divide; while wildtype cells divided by 120 min, GapR-depleted cells did not divide and became progressively more filamentous over time (Figure S4B). Third, the origin-terminus ratio in GapR-depleted cells remained constant after 60 min (Figure S4C), indicating that replication initiated but did not complete. Taken together, our results demonstrate that GapR is required for the timely and successful completion of DNA replication.
Cell 175, 583–597, October 4, 2018 587
B
A
WT ∆gapR Pxyl-gapR + glu (2 hr)
isolate G1 cells (1N) release
deplete in glucose (2 hr)
120 min 100 min 80 min
fresh media
60 min 1N
2N
G1
S
40 min
G2/M
sample every 20 min
20 min
stain DNA, measure DNA content with flow cytometry or DNA sequencing
D replication rate (% of WT)
fraction initiated
WT
0.6 0.4 ∆gapR Pxyl-gapR +glu (2 hr)
0.2 0.0 0
20
40
60
80
100
120
time post release (min)
F
1.6
60
1.2
WT
+ xyl->glu
+ glu (2 hr)
1.4
2.0 1.8 1.6 1.4 1.2 1.0
20 40 60 time post release (min)
20 min ter
ori
1.8
0
35 WT
∆gapR Pxyl-gapR xyl->glu
20
ori
0.8 2.0
20
25
time post release
1.0
40
30
ter WT ∆gapR Pxyl-gapR + glu (2 hr)
1.8
80
G replication rate (kb/min)
relative DNA content (ori)
100
0 min
ori
2.0
WT ∆gapR Pxyl-gapR +glu (2 hr)
∆gapR Pxyl-gapR
1N 2N
E
1.0 0.8
2N
1.6
relative DNA content
C
1N
40 min
1.4 1.2 1.0
ter
0.8 2.0
ori
1.8 1.6
60 min
1.4 1.2 1.0
15 10
∆gapR Pxyl-gapR glu (2 hr)
5
ter
0.8 2.0
ori
1.8 1.6
0
80 min
1.4
20-40
40-60
60-80
time range (min post release)
1.2 1.0
ter
0.8
3043
0
1000
2000
3000
chromosome position (kb)
Figure 3. Replication Initiation and Elongation Are Strongly Impaired in GapR-Depleted Cells (A) Assay to assess DNA replication. Wild-type and DgapR Pxyl-gapR cells were grown for 2 hr in glucose, synchronized in G1, and released into medium containing glucose. DNA content was measured by SYTOX staining and flow cytometry or DNA sequencing. (B) Flow cytometry time courses for wild-type and GapR-depleted cells. Blue lines represent 1N and 2N DNA content. (C) Fraction of cells from (B) that initiated replication, with DNA content >1N used as a proxy for initiation. Data are mean ± SEM, n R 2. (D) Replication rates for cells in (B), calculated from a linear fit of the DNA content of cells, normalized to the wild-type. Data are mean ± SEM, n R 2. (E) DNA sequencing of wild-type cells grown in glucose and DgapR Pxyl-gapR cells grown for 2 hr in glucose. DNA content, normalized to synchronized G1 cells, is shown as a function of chromosome position using LOWESS smoothing. Blue lines, replication origin and terminus. Cartoon of wild-type fork progression is shown (right). (F) Replication initiation, based on origin DNA content over time, using DNA sequencing data in (E). (G) Rates of replication calculated from a linear fit of fork position (from E) in the indicated time ranges. See also Figures S3 and S4.
GapR-Depleted Cells Are Sensitized to Inhibitors of Type II Topoisomerases Because our ChIP experiments suggested that GapR associates with (+) supercoiled DNA, we hypothesized that GapR promotes DNA replication by helping to alleviate the (+) supercoils that accumulate ahead of replication forks. This model predicts that cells depleted of GapR should be sensitive to sub-lethal doses of antibiotics that inhibit gyrase and topo IV, the enzymes that remove (+) supercoils. To test this, we measured the plating
588 Cell 175, 583–597, October 4, 2018
viability of the GapR depletion strain grown in xylose to reduce GapR levels (Figure 1C) after treatment with sub-lethal levels of novobiocin, an antibiotic that inhibits gyrase and topo IV. Notably, the GapR depletion strain exhibited nearly 5-logs lower plating viability than wild-type in novobiocin (Figure 4A). The GapR depletion strain was also sensitized to rifampicin and hydroxyurea (Figure S4D), which also inhibit DNA replication. However, the GapR depletion strain was not more sensitive than the wild-type to low doses of mitomycin C, which directly
A +novobiocin
+mitomycin C
+chloramphenicol
WT ∆gapR Pxyl-gapR ∆recA
B WT
WT
∆gapR Pxyl-gapR time post release (min) 0 +novo
∆gapR Pxyl-gapR
+novo
*
# of cells
*
20 * 40 60 80 100 120
*novobiocin added 1N
2N
1N
1N
2N
2N
1N
2N
DNA content
C DNA content (a.u.)
15 13
*
11 9 7 0
20
40
60
80
100
120
replication rate (% of WT)
D WT WT + novo ∆gapR Pxyl-gapR ∆gapR Pxyl-gapR + novo
17
100 80 60 40 20 0
time post release (min)
-
+ WT
-
+ novobiocin
∆gapR Pxyl-gapR
Figure 4. GapR-Depleted Cells Are Sensitized to Inhibition of Type II Topoisomerase Activity (A) 10-fold serial dilutions of wild-type, DgapR Pxyl-gapR + xylose, and DrecA cells on plates containing xylose and 0.1 mg/mL novobiocin, mitomycin C, or chloramphenicol. (B) Wild-type and DgapR Pxyl-gapR cells were grown in xylose, synchronized in G1, and released into medium containing xylose. After 20 min, 2.5 mg/mL novobiocin was added (red *). DNA content was measured with flow cytometry. Blue lines represent 1N and 2N DNA content. (C) DNA content of cells from (B). *Time of novobiocin addition. The fraction of cells with DNA content >1N was used as a proxy for initiation. Data are mean ± SEM, n R 2. (D) Relative replication rate of cells from (B), calculated from a linear fit of the DNA content of cells, normalized to the wild-type. Data are mean ± SEM, n R 2. See also Figure S4.
damages DNA, or chloramphenicol, which inhibits translation (Figure 4A). The antibiotic sensitivities of the GapR depletion strain contrast with that of DrecA cells, which had nearly identical plating viability as the wild-type on novobiocin and chloramphenicol, but was highly sensitized to mitomycin C. Thus, GapR is likely not required to repair DNA damage or for a general stress response. Instead, these results support a model in which GapR promotes DNA replication by relieving (+) supercoils in front of the elongating replisome. To test whether the synergistic effects of reduced GapR and novobiocin treatment on viability result from defects in DNA replication elongation, we used flow cytometry to monitor DNA content in synchronous populations of wild-type and GapRdepleted cells in the presence of novobiocin. Whereas wildtype cells continued to replicate at a slightly reduced rate in
novobiocin, cells with reduced levels of GapR showed very little accumulation of DNA and never reached a 2N state when treated with novobiocin (Figures 4B and 4C), having a rate of replication elongation 20% that of untreated wild-type cells (Figure 4D). These data confirm that without sufficient GapR, the replication fork becomes highly sensitized to topoisomerase inhibition, suggesting that GapR may affect topoisomerase activity. GapR Binds DNA and Can Alter DNA Topology Collectively, our in vivo results suggested that GapR promotes DNA replication by attenuating (+) supercoils, thereby relieving the topological stress associated with fork progression. To directly test whether GapR binds and affects DNA topology, we expressed and purified untagged GapR. Gel filtration of GapR produced a single peak that eluted faster than expected
Cell 175, 583–597, October 4, 2018 589
210 bp
150 bp
210/150 bp mix
0 0.01 0.05 0.1 0.25 0.5 1.0 2.5 5.0 10.0 20.0
[GapR] (μM)
0 0.05 0.1 0.25 0.5 1.0 2.5 5.0 10.0 20.0
0 0.001 0.025 0.05 0.1 0.25 0.5 1.0 2.5 5.0 10.0 20.0
[GapR] (μM)
C
16bp
0 0.01 0.05 0.1 0.25 0.5 1.0 2.5 5.0 10.0 20.0
B
300bp
0 0.01 0.05 0.1 0.25 0.5 1.0 2.5 5.0 10.0 20.0
A
D
[GapR] (μM)
E
relaxed plasmid
nicked plasmid GapR
+
0 0.25 0.5 1.0 2.5 5
30°C, 30 min
F
+
+
+
+
+
+
+
[GapR] 10 20 (μM) + + topo I
GapR
[GapR]
30°C, 30 min
+ topo I
+ ligase
+ Prot K
+ Prot K
0 0.25 0.5 1.0 2.5 5
10 20 (μM)
+
+
+
+
+
+
+
+ ligase
(2) + 4.5 μg/ml chloroquine
L
+
-
N
(+)
R
(-)
(-)
(+)
L
(1)
[GapR] (μM) 0
+
(-)
0.25
(+)
0.5
1.0
2.5
5
10
20
* *
(+)
Figure 5. GapR Binds DNA and Alters DNA Topology (A–C) EMSAs of GapR binding to linear dsDNA. Purified GapR was incubated with a 300 (A) or 16 bp probe (B), or with 210 and 150 bp probes individually (left and middle, C) and together (right, C). (D) GapR alters DNA topology. GapR was incubated with relaxed plasmid and then treated with calf thymus topoisomerase I (far left) before reactions were quenched, deproteinized with Proteinase K, electrophoresed, and imaged. Agarose gel analysis of plasmid topology (right), with supercoiled and relaxed standards (left). (E) GapR alters DNA topology. GapR was incubated with nicked plasmid, treated with T4 DNA ligase (far left) before reactions were quenched and analyzed as in (D) (right). Supercoiled, relaxed, and nicked standards are shown (left). (F) GapR constrains positive supercoils. Ligation products from (E) were analyzed by 2D-chloroquine electrophoresis (bottom, see STAR Methods). Migration of different plasmid forms are diagrammed (left) with standards shown (top): N, nicked; R, relaxed; L, linear; (), negatively supercoiled; (+), positively supercoiled. See also Figure S5.
for a 10 kDa protein (Figure S5A). Size-exclusion chromatography coupled to multi-angle light scattering (SEC-MALS) revealed a single, monodisperse (Mw/Mn = 1.000) peak with an observed molecular weight of 26.4 kDa (Figure S5A). To investigate how GapR binds DNA, we performed electrophoretic mobility shift assays (EMSAs) in the presence of 300 bp, linear double-stranded DNA derived from a GapR ChIP peak (Figure 5A). We observed the formation of GapR-DNA complexes, with multiple species observed at sub-saturating levels of GapR, with an apparent KD of 220 nM (Figure S5C). As we progressively reduced the probe length, the number of shifted species also decreased, with two species for a 20 bp probe and a single shift for a 16 bp probe (Figures 5B and S5B), suggesting that GapR binds to an 10 bp stretch of DNA. For 590 Cell 175, 583–597, October 4, 2018
probes >16 bp, the multiple bands seen could reflect multiple binding events or GapR-mediated bridging of different probes. To test the latter possibility, we performed an EMSA on a mixture of probes of two different lengths (Figure 5C). The resulting EMSA pattern was a union of the individual EMSAs, indicating that GapR does not bridge DNA in vitro. The apparent KD of GapR on DNA was only slightly (3-fold) lower on a 20 bp probe compared to a 300 bp probe or a 3 kb plasmid (Figure S5C), suggesting that GapR does not spread along DNA. To test whether GapR affects DNA topology, we incubated GapR with relaxed plasmid DNA. If GapR introduces or stabilizes any (+) supercoils that arise, then compensatory () supercoils will accumulate elsewhere in the plasmid to maintain a net relaxed state. Addition of calf-thymus topoisomerase I (topo I)
will then remove the () supercoils; subsequent Proteinase K treatment will remove topo I and GapR, yielding plasmid DNA harboring any supercoils introduced or stabilized by GapR. Changes in supercoiling were assessed by agarose electrophoresis, with each band representing a unique topoisomer. Notably, when both GapR and topo I were added, we observed conversion of relaxed plasmids into supercoiled forms (Figure 5D), but only when the concentration of GapR was >2.5 mM. At the highest concentration of GapR (20 mM), we saw reduced changes in DNA topology as GapR may coat the plasmid, preventing topo I from accessing the DNA (Figure S5D). Indeed, we found that such concentrations of GapR also protected plasmid DNA from DNase digestion (Figure S5F). We also tested whether GapR could affect DNA topology using a similar assay in which a circular, nicked plasmid was first incubated with GapR and then treated with T4 DNA ligase to trap any supercoils introduced or constrained by GapR; any compensatory supercoils induced by GapR binding would dissipate because of the nick present before ligase is added (Figure 5E). As above, bound proteins were then degraded and plasmid topology examined by agarose electrophoresis. Ligating nicked plasmid in the absence of GapR yielded relaxed plasmid, as expected. However, in the presence of GapR, we observed a clear increase in plasmid migration rates, again at concentrations of GapR >2.5 mM. To determine if these GapRconstrained topoisomers were (+) or () supercoils, we performed two-dimensional chloroquine gel electrophoresis. In the presence of >2.5 mM GapR, the ligated DNA migrated as expected for a (+) supercoiled plasmid, indicating that GapR can constrain or stabilize (+) supercoils (Figure 5F). GapR-mediated changes in (+) supercoiling could arise from a preference for binding DNA with (+) writhe or for overtwisted DNA, which can convert into (+) writhe. We favored the latter because the increases in (+) writhe measured by gel electrophoresis (Figures 5D–5F) required GapR concentrations >2.5 mM. As GapR binds DNA with a KD of 200–300 nM, many GapR binding events are likely required to generate the abrupt change in plasmid writhe observed (Figures 5D–5F and S5C). The simplest explanation is that each GapR binding event introduces or stabilizes a small amount of twist in the DNA, with twist from multiple binding events inducing a change in writhe. Importantly, however, GapR probably does not directly introduce substantial numbers of supercoils in vivo as immunoblots estimated only 3,000 GapR molecules per cell (Figure S1F). Instead, our findings collectively suggest that GapR preferentially recognizes and associates with overtwisted, (+) supercoiled DNA, as occurs in front of active replication forks and at the 30 ends of highly expressed genes. GapR Forms a Dimer-of-Dimer Clamp that Encircles Overly Twisted DNA GapR shows no homology to any structurally characterized DNA binding protein and, in fact, contains no known DNA binding motif. To determine how GapR binds DNA and recognizes overtwisted DNA, we solved a crystal structure of C. crescentus GapR in complex with DNA. We obtained this structure by single wavelength anomalous diffraction (SAD) using crystals obtained with selenomethionine substituted GapR in complex with an 11-mer, AT-rich site and refined the structure to final Rwork/Rfree
values of 24.1%/26.8% to 2.3 A˚ resolution (Figures 6A and 6B; Table S4). Strikingly, the structure shows that GapR oligomerizes to form a dimer-of-dimers, which completely encircles the 11-mer DNA inside a large, central cavity (Figure 6A; Video S1), consistent with our EMSA indicating a binding site of 10 bp (Figure 5B). Each GapR subunit consists of three helices. An N-terminal helix (residues 14–51), a1, forms an antiparallel coiled coil dimer with a1 of a second GapR subunit. The a1-a10 dimer (where 0 indicates the other subunit in a dimer) buries 1,340 A˚2 of protein surface from solvent. a2 (residues 55–66) and a3 (residues 69–86) are positioned at opposite ends of the a1-a10 antiparallel coiled coil dimer. The a2-a3 helical pairs from each GapR subunit interdigitate with a a2-a3 helical pair from a second GapR dimer to form the dimer-of-dimers. These cross a2-a3 zipper interactions seal the ends of the dimer-of-dimers, creating the large central DNA binding cavity (Figure 6A). An electrostatic surface representation of the GapR dimer-of-dimers shows that it forms a DNA-binding clamp with a remarkably electropositive inner DNA binding surface (Figure 6C; Video S1). DALI searches revealed no structures with similarity to the GapR monomer, dimer, or dimer-of-dimers. The electron density for the GapR-bound DNA is clearly resolved (Figure 6B) and forms a pseudocontinuous ladder in the crystal, resulting in multiple dimers-of-dimers along the DNA (Figure 6D). There are no significant contacts between any of the dimers-of-dimers, supporting the notion that GapR does not cooperatively coat or bridge DNA. Consistent with the idea that GapR does not bind DNA sequence specifically, GapR makes no base contacts and only interacts with the DNA phosphate backbone. Residues from both a1 and a2 helices make key contacts with the DNA. Specifically, Lys34, Lys42, and Lys49 from each a1 helix and Lys56, Lys59, Arg63, and Lys66 from the basic regions of the four a2 helices provide a total of 28 contacts to the phosphates of DNA (Figure S6A; Video S1). The a1 helices (Figure 6A) and these basic residues track along the phosphate backbone suggesting that they ‘‘read’’ a specific DNA conformation. Indeed, analysis of GapR-bound DNA revealed a significantly narrowed minor groove of 9.8 A˚ (compared to 12 A˚ for B-DNA) and widened major groove of 21 A˚ (compared to 17.2 A˚ for B-DNA) (Figure 6E). However, most strikingly, GapR-bound DNA is significantly overtwisted in the structure, with the helical DNA twist being 44.5 compared to 36.0 for B-DNA (Figure 6E). Importantly, attempts to fit B-DNA, which has a greater diameter than overtwisted DNA, within the GapR DNA binding cavity resulted in clash. Thus, the GapR DNA binding clamp appears to recognize an overtwisted DNA conformation that features narrow minor grooves and widened major grooves, not a specific DNA sequence. Notably, the oligonucleotide used in the crystal structure had high AT-content, which is known to more readily adopt an overtwisted state with narrow minor grooves. GapR does not bind as well to oligonucleotides with reduced AT-content (Figure S6D). To investigate whether apo GapR adopts a different conformation, we sought the apo GapR structure. Although apo C. crescentus GapR failed to produce data quality crystals, we obtained the structure of apo Bosea sp. Root381 GapR (herein termed Bosea GapR), which shares 60% sequence identity with the C. crescentus protein (Figure S6B). We confirmed that
Cell 175, 583–597, October 4, 2018 591
A
1 14
N
51 55 66 69
α1
α2
B
86 89
α3
C
dimer 2 dimer 1
C
D
dimer-ofdimer 1
E
GapR-bound DNA
major groove (21 Å)
minor groove (9.8 Å)
twist 44.5°
B-DNA
F
GapR-DNA dimer
dimer-ofdimer 2
dimer-ofdimer 3
apo GapR dimer
major groove (17.2 Å)
minor groove (12 Å)
twist 36°
Figure 6. Crystal Structures of apo GapR and GapR-DNA Complex (A) Structure of the C. crescentus GapR-DNA complex. Diagram of GapR helices (top left). GapR subunits are colored red and cyan (dimer 1), green and yellow (dimer 2). DNA is shown as a gray spheres model. Views are down the axis of the DNA (bottom left) and from the side (right). (B) Cut away view of the GapR-DNA complex with a sigma A weighted electron density 2Fo-Fc map (blue mesh) contoured at 1.1 s around DNA. (C) Electrostatic surface representation of the GapR dimer-of-dimers, with the centrally bound DNA shown as a cartoon. Blue and red represent electropositive and electronegative regions, respectively. (D) Packing of the GapR dimer-of-dimers in the crystal showing how multiple GapR oligomers could bind a long overtwisted region of DNA. (E) Comparison of GapR-bound DNA with B-DNA. The widths of the major and minor grooves and the twist are indicated on GapR-bound DNA and B-DNA. (F) Comparison of DNA-bound and apo GapR dimers. Structures of the DNA-bound GapR and apo GapR a1-a10 dimers with one subunit colored cyan and the other red (top). Also labeled is Asp68, located between a2 and a3. Close up of the Asp68 side chain, solvent-exposed as part of a continuous helix in the apo structure, while forming N-cap interactions with the amide nitrogen of a3 when GapR binds DNA (bottom). See also Figure S6 and Table S4.
592 Cell 175, 583–597, October 4, 2018
the Bosea GapR ortholog also adopts a dimer-of-dimers conformation with DNA in the central cavity (see STAR Methods, Table S4, and Figure S6C). The apo Bosea GapR structure was solved by SAD phasing and refined to final Rwork/Rfree values of 22.2%/ 25.9% to 2.4 A˚ resolution (Table S4). Interestingly, the apo Bosea GapR forms a dimer rather than a dimer-of-dimers (Figure 6F). The apo structure harbors the same a1-a10 antiparallel dimer interaction observed in the DNA-bound Caulobacter GapR structure, but the a2 and a3 helices form a continuous helix (Figure 6F). Although we did not obtain the apo Caulobacter GapR structure, our SEC-MALS analysis was consistent with Caulobacter GapR forming an elongated dimer in solution (Figure S5A). The finding that GapR is dimeric in the absence of DNA indicates that DNA is not threaded through the dimer-of-dimer clamp, which is also implausible, but rather it adopts the dimer-of-dimers state upon binding target DNA. GapR Stimulates the Activities of Topo IV and Gyrase Collectively, our results indicate that GapR recognizes overtwisted DNA such as that arising in front of replication forks. At these locations GapR may stimulate gyrase and/or topo IV to relieve superhelical stress and promote replication elongation. To test whether GapR promotes the activity of either gyrase or topo IV, we purified the two subunits of both Caulobacter gyrase and topo IV, and reconstituted each enzyme. We first tested whether GapR can stimulate topo IV to remove (+) supercoils, deliberately selecting reaction conditions that yield relatively low topoisomerase activity to maximize our ability to observe any stimulation by GapR. We found that GapR at 0.4 mM stimulated the rate of (+) supercoil relaxation by topo IV 4-fold (Figure 7A). After 60 min, a reaction lacking GapR still had 60% of the initial supercoiled substrate remaining, whereas the reaction containing GapR had almost no starting substrate left. Notably, the concentrations of GapR that stimulated topo IV activity (Figure S7A) are insufficient to introduce (+) supercoils alone (Figures 5D, 5E, and S5E). We also tested whether GapR stimulates topo IV to relax () supercoils and observed a modest, but reproducible, 1.3-fold increase in the rate of this reaction (Figure S7B). The stimulatory effects of GapR on topo IV were specific, as other proteins, including glucose-6-phosphate dehydrogenase (G6PD) and the DNAbinding protein GcrA, did not stimulate topo IV activity (Figure S7C). We also tested whether GapR can stimulate other bacterial type II topoisomerases, by repeating relaxation assays with E. coli topo IV. We found that 0.4 mM GapR was unable to stimulate (+) supercoil relaxation by E. coli topo IV, but produced a modest increase in () supercoil relaxation (Figure S7D). We then tested whether GapR stimulates Caulobacter gyrase. As with topo IV, we observed a clear and reproducible, GapRdependent increase of 2-fold for the rates at which gyrase relaxes (+) supercoiled DNA and introduces () supercoils into relaxed DNA (Figures 7B, S7E, and S7F). The stimulatory effects on gyrase were specific as neither G6PD nor GcrA stimulated gyrase (Figure S7G). GapR had only a modest, but reproducible, stimulatory effect on (+) supercoil relaxation and () supercoil introduction by E. coli gyrase (Figure S7H). Collectively, these results demonstrate that GapR promotes the activities of both type II topoisomerases from Caulobacter, gyrase and topo IV,
including their removal of (+) supercoils. Taken together with our in vivo studies, we conclude that the association of GapR with overtwisted DNA ahead of replication forks helps stimulate the type II topoisomerases, promoting replication fork progression, the timely removal of (+) supercoils ahead of replisomes, and, consequently, the completion of DNA replication in Caulobacter. DISCUSSION The successful completion of DNA replication in any organism requires the resolution of complex topological structures that build up around the fork. Although the importance of topoisomerases in this process is well established, how topoisomerases are recruited to, and how their activity is controlled at, the replication forks is still poorly understood. Here, we identified GapR as a crucial regulator of topoisomerases and replisome progression in a-proteobacteria. Loss of GapR sensitizes cells to the inhibition of type II topoisomerases and leads to chromosomal relaxation, significantly slowed replication, and eventual cell death. We propose that the essential function of GapR is to promote the activities of type II topoisomerases, both to maintain global supercoiling levels and, importantly, to attenuate the (+) supercoils that would otherwise impede replication forks (Figure 7C). Type II topoisomerases have long been known to enable DNA replication by alleviating the torsional stress that accumulates near replication forks. Indeed, the efficacy of many front-line antibiotics (e.g., coumerins and quinolones) is based on their ability to inhibit gyrase and thereby block bacterial DNA replication (Vos et al., 2011). Although type II topoisomerases such as gyrase or topo IV are required for DNA replication in vivo, whether these enzymes are sufficient has not been established. Our results indicate that they may not be sufficient, and instead may require activators or co-factors, such as GapR. Recognition of Highly Twisted DNA by GapR GapR specifically binds to overtwisted DNA (Figures 5E, 5F, and 6E), which is precisely the form of DNA that initially arises in front of replication forks. Overtwisted DNA also arises in front of translocating RNA polymerase and our ChIP-seq results demonstrated that GapR binds such regions in a sequenceindependent, but transcription-dependent manner. Consistent with these results, the structure of GapR revealed that the protein makes no base-specific contacts with DNA, but instead recognizes the shape of the DNA. When bound to DNA, GapR adopts a structure with an internal cavity lined with basic residues that contact the negatively charged backbone of the encircled DNA. Importantly, our results indicate that GapR is only able to recognize overtwisted DNA, as the larger diameter of B-DNA cannot be accommodated within the DNA-binding cavity of GapR. Thus, the structure of GapR in complex with DNA reveals a simple, but compelling mechanism for its localization to regions of the genome that are more highly twisted by the action of translocating DNA or RNA polymerase. Our ChIP-seq analyses and prior studies of GapR revealed that GapR shows a preference for binding AT-rich DNA (AriasCartin et al., 2017; Ricci et al., 2016). Notably, such DNA often contain A-tracts that can adopt intrinsically bent structures,
Cell 175, 583–597, October 4, 2018 593
B
A positive (+)
relaxed
positive (+)
+topo IV
topo IV
relaxed
+gyrase
gyrase
(+)
0
7.5
15
30
100
-GapR
80
*
60 40
+GapR
20 0
0
10
20 30 40 time (min)
50
+GapR 60
*
C
(+)
60 time (min)
replication
% initial supercoiled DNA
% initial supercoiled DNA
(+)
(+)
0
6
12 25 38 55 72 time (min)
100 80
-GapR
*
60 40 20 0
+GapR 0
+GapR
10 20 30 40 50 60 70 time (min)
*
rapid fork progression
GapR stimulates gyrase (or topo IV) relaxes (+)sc
fork progression leads to overtwist
+ gyrase
overtwisted DNA
stimulates
+ GapR
GapR localizes to overtwisted DNA
overtwist converts to (+) writhe
+ +
GapR
+
+
(+) supercoiling blocks fork progression
+ +
Figure 7. GapR Stimulates Gyrase and Topo IV Activity (A) GapR stimulates relaxation of (+) supercoils by topo IV. Topo IV was incubated with GapR (0.4 mM) before adding (+) supercoiled pUC19 and ATP. Reactions were quenched at the times indicated and analyzed by gel electrophoresis (right). Slower migration in the gel indicates relaxation of (+) supercoils. Relaxed and (+) supercoiled standards are shown (left of the time course). Plot quantifies the % initial supercoiled plasmid (asterisk) at each time point (left). Data indicate mean ± SEM, n R 2. (B) GapR stimulates relaxation of (+) supercoils by gyrase. Assay performed as in (A) except using gyrase. Plot quantifies % initial supercoiled plasmid (asterisk) at each time point (left). (C) Model of GapR and topoisomerases in replication. Replication fork progression generates overtwisted DNA, which can convert to (+) writhe and block fork movement. GapR binds to overtwisted DNA and stimulates gyrase and/or topo IV to relieve (+) writhe, allowing replication to resume. See also Figure S7.
with narrow minor grooves and increased superhelicity (Haran and Mohanty, 2009). Thus, AT-rich regions and the DNA ahead of replisomes or RNA polymerase may be the primary binding sites for GapR by virtue of a common, overtwisted shape. Our biochemical studies showed that GapR can also associate with relaxed or () supercoiled DNA, likely via local regions of overtwisted DNA within these substrates, arising from sequences with intrinsically narrow minor grooves or from thermodynamic fluctuations in twist. GapR forms a snug clamp around overtwisted DNA raising the question of how it is loaded. The apo structure revealed that
594 Cell 175, 583–597, October 4, 2018
GapR is a dimer in the absence of DNA, with a2 and a3 forming a continuous helix (Figure 6F), suggesting that the region between a2 and a3 is kinked upon DNA binding. A specific residue, Asp68, appears to function as the pivot point between a2 and a3. Among GapR homologs, an aspartic acid is highly conserved at this position, although glutamic acid, glutamine, and asparagine are also found. Notably, these residues are all known to function as N-cap residues and helix breakers. Indeed, while Asp68 points into the solvent in the apo GapR structure as part of the continuous a2 helix, it forms stabilizing N-cap contacts with a3 in the DNA-bound structures (Figure 6F). Combined, these
data suggest a unique mechanism for GapR clamp formation and loading in which the protein uses its a1-a10 helices to track along the DNA as a caliper to search for sites with narrow minor and expanded major grooves, which then results in stable dimerof-dimers formation through kinking to form a2-a3 and subsequent interactions between a2-a3 pairs around the DNA. Thus, unlike other DNA binding clamps, GapR likely does not require a clamp loader to dock onto its target DNA. The preferential binding of GapR to highly twisted DNA explains why many GapR binding events were required to generate increases in plasmid writhe (Figures 5D and 5E). Because supercoiling is a combination of twist and writhe, each GapR binding event introduces a small amount of twist into DNA, such that multiple binding events are necessary before inducing a change in writhe. Assuming GapR binds a 10 bp site with half the available sites bound, the change of +3 writhe we observed implies an overtwisting of the DNA by 8 per GapR bound, similar to the 8.5 change measured in the structure. However, GapR likely does not directly introduce large numbers of supercoils in vivo as there are only 3,000 GapR molecules present per cell (Figure S1F), significantly less than NAPs such as HU, which is present at 60,000 copies per cell (Rouvie`re-Yaniv et al., 1979). Additionally, we found that GapR stimulates gyrase and topo IV at lower concentrations than is required to introduce topological changes. These observations, along with the structure of GapR bound to DNA, indicate that the primary role of GapR is not to directly alter DNA topology, but rather to associate with regions of (+) supercoiling, promoting the activity of type II topoisomerases precisely where their activities are needed the most. Role of GapR in Modulating Topoisomerase Activity How, mechanistically, might GapR promote gyrase and topo IV activity? The simplest model is that GapR directly binds gyrase and topo IV to promote their binding to DNA or catalytic activity (Figure S7I, left). By recognizing overtwisted DNA that arises in front of replication forks, GapR may recruit or stimulate gyrase and topo IV to act on the (+) supercoiled DNA that arises in front of forks and that must be dissipated for replication to progress. Notably, the structure of GapR revealed two acidic patches that are solvent-exposed in the DNA-bound, dimer-of-dimers configuration. Most of these acidic residues are highly conserved (Figure S6B) and they may mediate a direct interaction with the topoisomerases. Alternatively, GapR could stimulate gyrase and topo IV without directly binding these enzymes, instead stabilizing DNA in a way that promotes topoisomerase binding and/or processivity (Figure S7I, right). For instance, (+) writhed DNA that arises in front of replication forks may normally diffuse away from the fork. GapR binding to overtwisted DNA could ‘‘pin’’ a plectonemic region of DNA, effectively trapping DNA for relaxation by topoisomerases. Further biochemical and enzymatic studies are needed to fully dissect the detailed mechanisms by which GapR stimulates gyrase and topo IV. Regulation and Control of DNA Replication Elongation Whether the loss of viability in GapR-depleted cells arises from the defects in replication or the broad changes in transcription, or both, is unclear. The transcriptional changes are, however, likely an indirect effect of GapR depletion as the GapR ChIP pro-
file did not correlate with the global changes in gene expression. We favor the notion, supported by our biochemical and structural studies, that the primary function of GapR is to promote type II topoisomerase activity during DNA replication. Consistent with this model, GapR-venus was suggested to localize to the DNA ahead of the replisome (Arias-Cartin et al., 2017). Other NAPs have been implicated in controlling bacterial DNA replication. For instance, E. coli IHF and Fis are both critical for DNA replication initiation (Magnan and Bates, 2015). There are also several proteins in E. coli, including MukB and SeqA, that interact with topo IV to stimulate () supercoil removal and chromosome decatenation, although neither is required for replication elongation (Hayama and Marians, 2010; Kang et al., 2003; Li et al., 2010). E. coli HU was suggested to stimulate gyrase’s decatenation activity in vitro, but whether this function of HU promotes DNA replication in vivo is unknown (Marians, 1987). E. coli YejK stimulates topo IV to relax (+) supercoils in vitro, but YejK inhibits gyrase and DyejK strains have only modest, possibly indirect, effects on DNA replication and exhibit no growth defect (Lee and Marians, 2013). Thus, to our knowledge, GapR is the first NAP required for replication elongation and the first characterized activator of gyrase. In sum, we suggest that GapR plays a critical role in resolving the superhelical stress associated with DNA replication. GapR may also be required to alleviate the torsional stress that arises in regions of the genome where (+) supercoils accumulate and cannot diffuse, e.g., when the replication forks converge near the terminus or, more frequently, when the transcription and replication machineries converge. At such sites, (+) supercoils must be resolved for replication to continue. GapR is highly conserved throughout the a-proteobacteria where it likely plays a similar role in activating type II topoisomerases. GapR homologs are not found in other classes of bacteria, but a role for NAPs in replication elongation has not been well explored. There also may be as yet unidentified NAPs that share no homology to GapR, but recognize highly twisted DNA to similarly promote replication. In eukaryotes, there may also be analogous proteins. For example, high mobility group (HMG) proteins bind DNA relatively non-specifically, contributing to chromosome architecture in poorly defined ways. Intriguingly, yeast Hmo1 and the type II topoisomerase Top2 were recently shown to bind together near genes transcribed during S phase to somehow suppress chromosome fragility (Bermejo et al., 2009), and human HMGB1 interacts with and stimulates Top2a decatenation activity (Stros et al., 2007). We anticipate that interactions between stably bound, chromosome architecture proteins and topoisomerases, like that described here for GapR, will ultimately prove essential to the maintenance of genome integrity and the successful execution of DNA replication in most organisms. STAR+METHODS Detailed methods are provided in the online version of this paper and include the following: d d
KEY RESOURCES TABLE CONTACT FOR REAGENT AND RESOURCE SHARING
Cell 175, 583–597, October 4, 2018 595
d
d
d
d
EXPERIMENTAL MODEL AND SUBJECT DETAILS B Growth conditions and chemical treatments B Strain construction B Plasmid construction METHOD DETAILS B Nucleoid purification and mass spectrometry analysis B Expression profiling with microarrays and RNAsequencing B Chromatin Immunoprecipitation Sequencing (ChIP-seq) B Caulobacter synchronizations B Immunoblotting B Microscopy B Flow cytometry B DNA sequencing B qPCR analysis of origin-terminus ratio B Plating viability and drug sensitivity B Phylogeny analysis B Protein expression and purification B SEC-MALS B Preparation of DNA substrates B DNA electrophoresis and visualization B Electrophoretic mobility shift assays B DNA topology assays B DNase I protection assay B Crystallization of C. crescentus GapR-DNA complex, Bosea sp. apo and Bosea sp. GapR-DNA-bound complex B Structure determination and refinement of the C. crescentus GapR-DNA complex, Bosea sp. apo and DNAbound complex B Fluorescence polarization-based DNA binding experiments B Topo IV relaxation assays B Gyrase supercoiling assays QUANTIFICATION AND STATISTICAL ANALYSIS B Calculation of growth rates B Quantification of protein levels by immunoblotting B Analysis of RNA sequencing B ChIP analysis B Analysis of replisome progression by flow cytometry B Analysis of DNA replication from DNA sequencing B Quantification of cell division B Quantification of origin-terminus ratio by qPCR B Quantification of DNA binding and topoisomerase assays DATA AND SOFTWARE AVAILABILITY
SUPPLEMENTAL INFORMATION Supplemental Information includes seven figures, four tables, and one video and can be found with this article online at https://doi.org/10.1016/j.cell. 2018.08.029. ACKNOWLEDGMENTS We thank J. Berger for expression constructs; S. Goodman for antibodies; and M. LeRoux, M. Guzzo, and S. Bell for comments on the manuscript. This work supported by NIH (R01GM115563 to M.A.S. and R01GM082899 to M.T.L.), a
596 Cell 175, 583–597, October 4, 2018
Howard Hughes Medical Institute (HHMI) International Predoctoral Fellowship (to D.L.H.), and a Jane Coffin Childs Memorial Fellowship (to M.S.G.). M.T.L. is also an Investigator of the HHMI. AUTHOR CONTRIBUTIONS D.L.H. performed the cell fractionation, RNA sequencing, and ChIP sequencing. M.S.G. performed the flow cytometry, DNA sequencing, and biochemical experiments. W.Z. purified proteins for structural studies. M.A.S. performed the structural studies. D.L.H., M.S.G., M.A.S., and M.T.L. conceived experiments, analyzed data, and wrote the manuscript. DECLARATION OF INTERESTS The authors declare no competing interests. Received: March 9, 2018 Revised: June 13, 2018 Accepted: August 15, 2018 Published: September 13, 2018 REFERENCES Arias-Cartin, R., Dobihal, G.S., Campos, M., Surovtsev, I.V., Parry, B., and Jacobs-Wagner, C. (2017). Replication fork passage drives asymmetric dynamics of a critical nucleoid-associated protein in Caulobacter. EMBO J. 36, 301–318. Badrinarayanan, A., Le, T.B.K., and Laub, M.T. (2015). Bacterial chromosome organization and segregation. Annu. Rev. Cell Dev. Biol. 31, 171–199. Bermejo, R., Capra, T., Gonza´lez-Huici, V., Fachinetti, D., Cocito, A., Natoli, G., Katou, Y., Mori, H., Kurokawa, K., Shirahige, K., and Foiani, M. (2009). Genome-organizing factors Top2 and Hmo1 prevent chromosome fragility at sites of S phase transcription. Cell 138, 870–884. Chen, V.B., Arendall, W.B., 3rd, Headd, J.J., Keedy, D.A., Immormino, R.M., Kapral, G.J., Murray, L.W., Richardson, J.S., and Richardson, D.C. (2010). MolProbity: all-atom structure validation for macromolecular crystallography. Acta Crystallogr. D Biol. Crystallogr. 66, 12–21. Haran, T.E., and Mohanty, U. (2009). The unique structure of A-tracts and intrinsic DNA bending. Q. Rev. Biophys. 42, 41–81. Hayama, R., and Marians, K.J. (2010). Physical and functional interaction between the condensin MukB and the decatenase topoisomerase IV in Escherichia coli. Proc. Natl. Acad. Sci. USA 107, 18826–18831. Higgins, N.P. (2016). Species-specific supercoil dynamics of the bacterial nucleoid. Biophys. Rev. 8 (Suppl 1 ), 113–121. Jensen, R.B. (2006). Coordination between chromosome replication, segregation, and cell division in Caulobacter crescentus. J. Bacteriol. 188, 2244–2253. Jones, T.A., Zou, J.Y., Cowan, S.W., and Kjeldgaard, M. (1991). Improved methods for building protein models in electron density maps and the location of errors in these models. Acta Crystallogr. A 47, 110–119. Kang, S., Han, J.S., Park, J.H., Skarstad, K., and Hwang, D.S. (2003). SeqA protein stimulates the relaxing and decatenating activities of topoisomerase IV. J. Biol. Chem. 278, 48779–48785. Khodursky, A.B., Peter, B.J., Schmid, M.B., DeRisi, J., Botstein, D., Brown, P.O., and Cozzarelli, N.R. (2000). Analysis of topoisomerase function in bacterial replication fork movement: use of DNA microarrays. Proc. Natl. Acad. Sci. USA 97, 9419–9424. Lal, A., Dhar, A., Trostel, A., Kouzine, F., Seshasayee, A.S.N., and Adhya, S. (2016). Genome scale patterns of supercoiling in a bacterial chromosome. Nat. Commun. 7, 11055. Lee, C., and Marians, K.J. (2013). Characterization of the nucleoid-associated protein YejK. J. Biol. Chem. 288, 31503–31516. Li, Y., Stewart, N.K., Berger, A.J., Vos, S., Schoeffler, A.J., Berger, J.M., Chait, B.T., and Oakley, M.G. (2010). Escherichia coli condensin MukB stimulates
topoisomerase IV activity by a direct physical interaction. Proc. Natl. Acad. Sci. USA 107, 18832–18837. Liu, L.F., and Wang, J.C. (1987). Supercoiling of the DNA template during transcription. Proc. Natl. Acad. Sci. USA 84, 7024–7027. Magnan, D., and Bates, D. (2015). Regulation of DNA replication initiation by chromosome structure. J. Bacteriol. 197, 3370–3377.
Rouvie`re-Yaniv, J., Yaniv, M., and Germond, J.E. (1979). E. coli DNA binding protein HU forms nucleosomelike structure with circular double-stranded DNA. Cell 17, 265–274. Stros, M., Bacı´kova´, A., Polanska´, E., Stokrova´, J., and Strauss, F. (2007). HMGB1 interacts with human topoisomerase IIalpha and stimulates its catalytic activity. Nucleic Acids Res. 35, 5001–5013.
Marians, K.J. (1987). DNA gyrase-catalyzed decatenation of multiply linked DNA dimers. J. Biol. Chem. 262, 10362–10368.
Taylor, J.A., Panis, G., Viollier, P.H., and Marczynski, G.T. (2017). A novel nucleoid-associated protein coordinates chromosome replication and chromosome partition. Nucleic Acids Res. 45, 8916–8929.
Modell, J.W., Hopkins, A.C., and Laub, M.T. (2011). A DNA damage checkpoint in Caulobacter crescentus inhibits cell division through a direct interaction with FtsW. Genes Dev. 25, 1328–1343.
Vos, S.M., Tretter, E.M., Schmidt, B.H., and Berger, J.M. (2011). All tangled up: how cells direct, manage and exploit topoisomerase function. Nat. Rev. Mol. Cell Biol. 12, 827–841.
Postow, L., Crisona, N.J., Peter, B.J., Hardy, C.D., and Cozzarelli, N.R. (2001). Topological challenges to DNA replication: conformations at the fork. Proc. Natl. Acad. Sci. USA 98, 8219–8226.
Weitao, T., Nordstro¨m, K., and Dasgupta, S. (2000). Escherichia coli cell cycle control genes affect chromosome superhelicity. EMBO Rep. 1, 494–499.
Ray-Soni, A., Bellecourt, M.J., and Landick, R. (2016). Mechanisms of bacterial transcription termination: all good things must end. Annu. Rev. Biochem. 85, 319–347. Ricci, D.P., Melfi, M.D., Lasker, K., Dill, D.L., McAdams, H.H., and Shapiro, L. (2016). Cell cycle progression in Caulobacter requires a nucleoid-associated protein with high AT sequence recognition. Proc. Natl. Acad. Sci. USA 113, E5952–E5961.
Worcel, A., and Burgi, E. (1972). On the structure of the folded chromosome of Escherichia coli. J. Mol. Biol. 71, 127–147. Wu, H.Y., Shyy, S.H., Wang, J.C., and Liu, L.F. (1988). Transcription generates positively and negatively supercoiled domains in the template. Cell 53, 433–440. Zheng, G., Lu, X.-J., and Olson, W.K. (2009). Web 3DNA–a web server for the analysis, reconstruction, and visualization of three-dimensional nucleic-acid structures. Nucleic Acids Res. 37, W240-6.
Cell 175, 583–597, October 4, 2018 597
STAR+METHODS KEY RESOURCES TABLE
REAGENT or RESOURCE
SOURCE
IDENTIFIER
Mouse monoclonal anti-FLAG
Sigma
Cat #: F1804; RRID: AB_262044
Mouse monoclonal anti-rpoA
Neoclone
Cat #: W0003; RRID: AB_2687386
Rabbit polyclonal anti-CtrA
Lab collection
Cat #: N/A
Rabbit polyclonal anti-GapR
This work
Cat #: N/A
Rabbit polyclonal anti-IHF
Gift from Stephen Goodman
Cat #: N/A
Goat anti-Rabbit IgG (H+L) Secondary Antibody, HRP
Invitrogen
Cat #: 32460; RRID: AB1185567
Goat anti-Mouse IgG (H+L) Secondary Antibody, HRP
Invitrogen
Cat #: 32430; RRID: AB_1185566
Lab collection
ML76
Antibodies
Bacterial and Virus Strains Caulobacter crescentus synchronizable wild-type CB15N (NA1000) ML76 gapR::gapR-3xFLAG
This study
ML2794
ML76 gapR::tet Pxyl-gapR
This study
ML2795
ML76 recA::gent
Lab collection
ML2796
E. coli C41(DE3)
Sigma
Cat #: CMC0017
E. coli BL21(DE3)
NEB
Cat #: C2527I
C41(DE3) pRET3a-Rgy
This study
ML2797
BL21(DE3) pET28-GapR
This study
ML2798
BL21(DE3) pKS22b-hSUMO-GapR
This study
ML2799
BL21(DE3) pKS22b-hSUMO-GyrA
This study
ML2800
BL21(DE3) pKS22b-hSUMO-GyrB
This study
ML2801
BL21(DE3) pKS22b-hSUMO-ParC
This study
ML2802
BL21(DE3) pKS22b-hSUMO-ParE
This study
ML2803
BL21(DE3) pET28-his-TEV-GyrA
Gift from James Berger
ML2892
BL21(DE3) pET28-his-TEV-GyrB
Gift from James Berger
ML2893
BL21(DE3) pET15b-GapR(Ccres)
This study
ML2894
BL21(DE3) pET15b-GapR(L83M)
This study
ML2895
BL21(DE3) pET15b-GapR(Bosea)
This study
ML2896
BL21(DE3) pET15b-GapR(L48M-I54M-L73M)
This study
ML2897
ML76 xyl::PgapR-gapR-kan gapR::tet
This study
ML2954
2X-HiRPM hybridization buffer
Agilent
Cat #: 5188-6420
Adenosine triphosphate (ATP)
Sigma
Cat #: A2383
Benzonase
Millipore
Cat #: 70746-3
Betaine
Sigma
Cat #: B2754
Bovine Serum Albumin (ultrapure)
Ambion
Cat #: AM2616
Brij-58
ThermoFisher
Cat #: 28386
Ethidium Bromide
Sigma
Cat #: E1510
Carbenicillin
Goldbio
Cat #: C-103-100
Chloramphenicol
Sigma
Cat #: 60378
Chloroquine phosphate
Santa Cruz Biotech
Cat #: sc-205629
CTAB
Sigma
Cat #: 85,582-0
Chemicals, Peptides, and Recombinant Proteins
(Continued on next page)
e1 Cell 175, 583–597.e1–e16, October 4, 2018
Continued REAGENT or RESOURCE
SOURCE
IDENTIFIER
Cy3-dCTP
Sigma
Cat #: GEPA53021
Cy5-dCTP
Sigma
Cat #: GEPA55021
DAPI (4’,6-diamidino-2-phenylindole)
Invitrogen
Cat #: D3571
Deoxycholate
Sigma
Cat #: D6750
DNA Polymerase I
NEB
Cat #: M0209S
dUTP
Promega
Cat #: U119A
Dynabeads Protein A for Immunoprecipitation
Invitrogen
Cat #: 10002D
E. coli DNA ligase
NEB
Cat #: M0205S
Formaldehyde
Sigma
Cat #: F8775-25ML
Gentamycin
Sigma
Cat #: G3632
Glycine
Sigma
Cat #: G8898
GlycoBlue
Ambion
Cat #: AM9516
Glycogen
Invitrogen
Cat #: 10814010
Isopropyl b-D-1-thiogalactopyranoside
Goldbio
Cat #: I2481C
Kanamycin
Sigma
Cat #: K1377
Kapa HiFi HotStart DNA Polymerase
Kapa Biosystems
Cat #: KK2502
Kapa 2X qPCR Master Mix
Kapa Biosystems
Cat #: KK4609
Klenow (large) fragment
NEB
Cat #: M0210L
Klenow Fragment (30 /50 exo-)
NEB
Cat #: M0212L
LiCl
Mallinckrodt
Cat #: 5852
Lysozyme
ThermoFisher
Cat #: 89833
Mitomycin C
Sigma
Cat #: M0503
Novobiocin
Sigma
Cat #: N1628
NP-40 substitute
Fluka
Cat #: 74385
Nt.BspQI
NEB
Cat #: R0644S
Oxytetracycline
Sigma
Cat #: T3383
Paraformaldehyde
Sigma
Cat #: P6148
PEG 8000
Hampton Research
Cat #: HR2-515
Percoll
GE Healthcare
Cat #: 17-0891-01
Phenol – chloroform – isoamyl alcohol mixture (25:24:1)
Sigma
Cat #: 77617-100ML
Phusion DNA polymerase
NEB
Cat #: E0553L
PMSF
ThermoFisher
Cat #: 36978
Proteinase K
NEB
Cat #: P8107S
Quick Ligase
NEB
Cat #: M2200L
Random Primers
Invitrogen
Cat #: 48190-011
Ready-Lyse Lysozyme Solution
Epicenter
Cat #: R1802M
Rifampicin
Sigma
Cat #: R3501
RNase A
QIAGEN
Cat #: 19101
RNase H
NEB
Cat #: M0297L
Roche cOmplete EDTA-free tablets
Sigma
Cat #: 11873580001
Second Strand Buffer
ThermoFisher
Cat #: 10812014
SIGMAFAST Protease Inhibitor Tablet
Sigma
Cat #: S8830-20TAB
SmaI
NEB
Cat #: R0141S
Spectinomycin
Sigma
Cat #: C90007
Spermidine
Sigma
Cat #: 85580
Superase-In
ThermoFisher
Cat #: AM2696 (Continued on next page)
Cell 175, 583–597.e1–e16, October 4, 2018 e2
Continued REAGENT or RESOURCE
SOURCE
IDENTIFIER
Superscript II
ThermoFisher
Cat #: 18064014
Superscript III
ThermoFisher
Cat #: 18080093
SYBR Gold
Invitrogen
Cat #: S11494
SYTOX Green Nucleic Acid Stain
Invitrogen
Cat #: S7020
T4 DNA Ligase
NEB
Cat #: M0202L
T4 DNA Polymerase
NEB
Cat #: M0203L
T4 PNK
NEB
Cat #: M0201L
TCEP-HCl
ThermoFisher
Cat #: 20491
TEV protease
Sigma
Cat #: T4455
Topoisomerase I (calf thymus)
Invitrogen
Cat #: 38042024
Topoisomerase I (E. coli)
NEB
Cat #: M0301L
Topoisomerase IV (E. coli)
TopoGEN
Cat #: TG2000EC-1
Triton X-100
Sigma
Cat #: T9284
Trizol
ThermoFisher
Cat #: 15596026
Turbo DNase
Ambion
Cat #: AM2238
USER enzyme
NEB
Cat #: M5505S
Yeast tRNA
Sigma
Cat #: R4018
Hydroxyurea
Sigma
Cat #: H8627
GcrA
Lab collection
N/A
Glucose-6-Phosphate Dehydrogenase
Lab collection
N/A
Direct-zol RNA MiniPrep
Zymo
Cat #: R2050
Pierce Coomassie Plus (Bradford) Assay Reagent
ThermoFisher
Cat #: 23238
QIAGEN Gel Extraction Kit
QIAGEN
Cat #: 28706
Critical Commercial Assays
QIAGEN PCR Purification Kit
QIAGEN
Cat #: 28106
QIAGEN Plasmid Mini Kit
QIAGEN
Cat #: 27106
Ribo-Zero Kit
Illumina
Cat #: MRZB12424
RNA fragmentation reagents
Ambion
Cat #: AM8740
SuperSignal West Femto Maximum Sensitivity Substrate
ThermoFisher
Cat #: 34095
SuperSignal West Femto Maximum Sensitivity Substrate
ThermoFisher
Cat #: 34095
RNeasy Mini Kit
QIAGEN
Cat #: 74104
Custom Agilent Arrays
Agilent, Lab collection
N/A
Caulobacter nucleoid mass spectrometry
This paper
Table S1
DgapR Pxyl-gapR + xyl microarray data
This paper
Table S2
Wild-type + novobiocin microarray data
This paper
Table S2
DgapR Pxyl-gapR + glu RNA-seq data
This paper
GSE100657
Deposited Data
UV and MMC microarray data
Modell et al., 2011
N/A
SocB induction microarray data
ArrayExpress
E-MEXP-3990
SciP depletion microarray data
Gene Expression Omnibus
GSE22062
Ethanol treatment RNA-seq data
Gene Expression Omnibus
GSE90030
ChIP-seq data
This paper
GSE100657
RNA polymerase ChIP
Gene Expression Omnibus
GSE73925
DNA-seq data
This paper
GSE100657
Apo Bosea sp. Root 381 GapR
This paper
6CFY (Continued on next page)
e3 Cell 175, 583–597.e1–e16, October 4, 2018
Continued REAGENT or RESOURCE
SOURCE
IDENTIFIER
Caulobacter crescentus GapR-DNA
This paper
6CG8
Apo Bosea sp. Root 381 GapR-DNA
This paper
6CFX
This paper
N/A
Oligonucleotides For primer table, see Table S3 Recombinant DNA pNPTS138-UP-gapR-3xFLAG-DW
This paper
N/A
pMT4427-CCNA_03428
This paper
N/A
pNTPS-SPEC-CCNA_03428-TET
This paper
N/A
pUC19
NEB
Cat #: N3041L
pRET3a-Rgy
Gift from James Berger
N/A
pET28-CCNA_03428
This study
N/A
pKS22-hSUMO-GapR
This study
N/A
pKS22-hSUMO-ParC
This study
N/A
pKS22-hSUMO-ParE
This study
N/A
pKS22-hSUMO-GyrA
This study
N/A
pKS22-hSUMO-GyrB
This study
N/A
pET28-his-TEV-GyrA
Gift from James Berger
N/A
pET28-his-TEV-GyrB
Gift from James Berger
N/A
pET15b-GapR(Ccres)
This study
N/A
pET15b-GapR(L83M)
This study
N/A
pET15b-GapR(Bosea)
This study
N/A
pET15b-GapR(L48M-I54M-L73M)
This study
N/A
pMT4260-PgapR-gapR
This study
N/A
ASTRA software 5.3.4
Wyatt Technology, 2017
https://www.wyatt.com/products/ software/astra.html
Bowtie
Bowtie
http://bowtie-bio.sourceforge.net/ index.shtml
Bowtie2
Bowtie2
http://bowtie-bio.sourceforge.net/ bowtie2/index.shtml
BWA
Burrows-Wheeler Aligner
http://bio-bwa.sourceforge.net/
DREME
MEME Suite
http://meme-suite.org/doc/dreme. html?man_type=web
Fiji/ImageJ
Fiji/ImageJ
https://fiji.sc/
FlowJo (X 10.0.7r2)
FlowJo LLC, 2013-2017
https://www.flowjo.com/
JalView
JalView
http://www.jalview.org/
Kaleigraph
Synergy Software
http://www.synergy.com/wordpress_ 650164087/kaleidagraph/
MACS software package 1.4.2
Model-based Analysis for ChIP-Seq
http://liulab.dfci.harvard.edu/MACS/
Mascot
Matrix Science
http://www.matrixscience.com/
Metamorph
Molecular Devices, 2017
https://www.moleculardevices.com/ systems/metamorph-research-imaging/ metamorph-microscopy-automationand-image-analysis-software
MicrobeJ
MicrobeJ by the Brun Lab
http://www.microbej.com/
Microbes Online
Virtual Institute for Microbial Stress and Survival
http://www.microbesonline.org/ operons/gnc565050.html
Software and Algorithms
(Continued on next page)
Cell 175, 583–597.e1–e16, October 4, 2018 e4
Continued REAGENT or RESOURCE
SOURCE
IDENTIFIER
Molprobity
Chen et al., 2010
http://molprobity.biochem.duke.edu/
O
Jones et al., 1991
N/A
Phenix Autosol
The PHENIX Industrial Consortium
https://www.phenix-online.org/
Prism 7
GraphPad Software, 2017
https://www.graphpad.com/scientificsoftware/prism/
Python (2.7.6, 2.7.3)
Python Software Foundation, 2017
https://www.python.org/download/ releases/2.7/
Sprot_extra_022412
Koch Institute Biopolymers and Proteomics Facility
https://ki.mit.edu/sbc/biopolymers
STRING database
String Consortium
https://string-db.org/cgi/input.pl
w3DNA
Zheng et al., 2009
https://w3dna.rutgers.edu/
Weblogo version 2.8
Weblogo
http://weblogo.berkeley.edu/
Custom Python scripts for data analysis
This study
N/A
15% Tris-HCl gel
BioRad
Cat #: 1611103
4-20% Tris-HCl gels
BioRad
Cat #: 4561096
8% TBE acrylamide gel
Invitrogen
Cat #: EC62155BOX
Agencourt Ampure XP beads
Beckman Coulter
Cat #: A63880
Amicon Ultra-15 Centrifugal Filter Units, 50 kD MWCO
Millipore
Cat #: UFC905024
Amicon Ultra-4 Centrifugal Filter Units, 3 kD MWCO
Millipore
Cat #: UFC800324
DNA retardation gels
Invitrogen
Cat #: EC63655BOX
Ni-NTA agarose resin
QIAGEN
Cat #: 30230
Oligo aCGH/Chip-on-chip wash buffer 1
Agilent
Cat #: 5185-5221
Other
Oligo aCGH/Chip-on-chip wash buffer 2
Agilent
Cat #: 5185-5222
PD-10 Desalting Columns
GE Healthcare
Cat #: 17085101
Spin-X 0.22 mm cellulose acetate column
Costar
Cat #: 8160
Stabilization and Drying Solution
Agilent
Cat #: 5185-5979
CONTACT FOR REAGENT AND RESOURCE SHARING Further information and requests for resources and reagents should be directed to and will be fulfilled by the Lead Contact, Michael T. Laub (
[email protected]). EXPERIMENTAL MODEL AND SUBJECT DETAILS Growth conditions and chemical treatments Caulobacter strains were grown in PYE (2 g/L bactopeptone, 1g/L yeast extract, 0.3 g/L MgSO4, 0.5 mM 0.5M CaCl2) at 30 C unless noted. The Pxyl promoter was induced by supplementing media with xylose (0.3%) and repressed with glucose (0.2%). Antibiotics were used at the following concentrations unless noted (liquid/plates): oxytretracycline (1 mg mL-1 / 2 mg mL-1), spectinomycin (25 mg mL-1 / 200 mg mL-1), kanamycin (5 mg mL-1 / 25 mg mL-1), gentamycin (NA / 5 mg mL-1). E. coli strains were grown in LB (10 g/L NaCl, 10 g/L tryptone, 5 g/L yeast extract) and supplemented with antibiotics at the following concentrations (liquid/plates): carbenicillin (50 mg mL-1 / 100 mg mL-1), oxytretracycline (12 mg mL-1 / 12 mg mL-1), spectinomycin (50 mg mL-1 / 50 mg mL-1), kanamycin (30 mg mL-1 / 50 mg mL-1), gentamycin (15 mg mL-1 / 20 mg mL-1). When required, 400 mM IPTG was used to induce gene expression unless noted. Optical density was measured at 600 nm using a Genesys 10 Bio Spectrophotometer. Strain construction All Caulobacter strains were derivatives of ML76, an isolate of strain CB15N/NA1000. The insertion of a 3xFLAG at the C terminus of GapR was constructed via a two-step recombination method using sacB as a counterselection marker. Strain ML2794 was constructed by electroporating plasmid pNPTS138-UP-gapR-3xFLAG-DW into ML76 and first integrants were selected on kanamycin plates. A second recombination step was performed by growing first integrants in PYE overnight and then plating on sucrose.
e5 Cell 175, 583–597.e1–e16, October 4, 2018
Sucrose resistant clones were restreaked to verify they were kan sensitive and colonies were verified by PCR with primer pair CCNA_3428-3xflag-2-R and CCNA_03428-UP-F-less-hindiii. The GapR depletion strain, ML2795, was constructed by electroporating plasmid pMT4427-CCNA_03428 into ML76 and selecting on kanamycin plates. As a second step, plasmid pNTPS-SPEC-CCNA_03428-TET was electroporated and cells were selected on kanamycin, spectinomycin and tetracycline. First integrants were then grown overnight with kanamycin, tetracycline, and xylose and then plated on kanamycin, tetracycline, xylose + sucrose plates for sacB counterselection. Sucrose resistant plates were restreaked to verify sensitivity to spectinomycin. ML2954 was constructed by electroporating plasmid pMT4260-PgapR-gapR into ML76 and selecting on kanamycin plates. GapR was subsequently deleted by integrating pNTPS-SPEC-CCNA_03428-TET as described previously. Plasmid construction Integration plasmids pNPTS138-UP-gapR-3xFLAG-DW was constructed by amplifying a fragment containing the upstream region of gapR with primers CCNA_03428-UP-F-less-hindiii and CCNA_03428-3xflag-1-R to which the coding region for 3x-FLAG was added by a second round of PCR with primers CCNA_3428-3xflag-2-R and CCNA_03428-UP-F-less-hindiii. This upstream fragment was fused by splice-overlap-extension (SOE) PCR to a fragment containing the region downstream of gapR amplified with primers 3428 dw-fuse-3xflag-F and CCNA_03428-DW-R. The resulting PCR product was digested with NheI and HindIII restriction enzymes and ligated into pNPTS-138 digested with the same enzymes. The pNTPS-SPEC-CCNA_03428-TET plasmid was constructed by multiple SOE PCR reactions. First, a fragment containing 703 bp upstream of CCNA_03428 was amplified with the primer pair CC3319TETDELUP-R and FOR_UP3319PNTPSPSTI, the tetracycline (tet) resistance cassette was amplified from template pKO3 using primer pair OL7_tet_F and OL8_tet_R, and a fragment containing 703 bp downstream of CCNA_03428 was amplified using primers CC3319TETDELDW-F and CC3319TETDELDW_SPEI-R. Subsequent SOE PCRs were performed to fuse the three fragments together. The resulting product was digested with PstI and SpeI and ligated into pNTPS-SPEC digested with the same restriction enzymes. The pMT4260-PgapR-gapR plasmid was constructed by amplifying a fragment containing the promoter and coding region of GapR with GapR_up_KpnI and GapR_down_NheI, digesting with KpnI and NheI, and ligated pMT4260 digested with the same enzymes. Replicative plasmids pMT4427-CCNA_03428 was constructed by ligation of a PCR fragment amplified with primers CCNA_03428-F-NdeI and CCNA_03428-R-SacI and digested with NdeI and SacI. Expression plasmids pET28-CCNA_03428 was constructed by ligation of a PCR fragment amplified with primers CCNA_3428-F-NdeI and CCNA_3428_R_SacI and digested with NdeI and SacI. pKS22b-hSUMO-GapR was constructed by ligation of a PCR fragment amplified with primers GapR-F-BamHI and GapR-R-NotI and digested with BamHI and NotI. pKS22b-hSUMO-GyrA was constructed by Gibson assembly cloning with a PCR fragment amplified with primers GyraseA-FBamHI-homology and GyraseA-R-NotI-homology and the pKS22b vector digested with BamHI and NotI. pKS22b-hSUMO-GyrB was constructed by Gibson assembly cloning with a PCR fragment amplified with primers GyrB-F-BamHI-homology and GyrB-RNotI-homology and the pKS22b vector digested with BamHI and NotI. pKS22b-hSUMO-ParC was constructed by Gibson assembly cloning with a PCR fragment amplified with primers ParC-F-BamHI-homology and ParC-R-NotI-homology and the pKS22b vector digested with BamHI and NotI. pKS22b-hSUMO-ParE was constructed by Gibson assembly cloning with a PCR fragment amplified with primers ParE-F-BamHI-homology and ParE-R-NotI-homology and the pKS22b vector digested with BamHI and NotI. For pET15b-GapR(Cres) and pET15b-(Bosea) and their variants, the full length Bosea sp. Root 381 GapR protein and C. crescentus GapR residues 11-89 were purchased as E. coli codon optimized genes from Genscript with NdeI and BamHI sites and subcloned into pET15b using restriction cloning. This resulted in the placement of a cleavable hexa-histidine tag at the N terminus of the expressed proteins. C. crescentus GapR(11-89) was used for structural studies as the N-terminal 10 residues of the C. crescentus protein are not conserved among GapR homologs (and not present in several homologs including the Bosea sp. Root 381 GapR protein) and are predicted to be disordered. METHOD DETAILS Nucleoid purification and mass spectrometry analysis The nucleoid purification protocol was performed as follows. A 25 mL culture of exponentially growing ML76 was centrifuged at 8,500 rpm at 4 C for 5 min. Cell pellets were resuspended in 0.5 mL of ice-cold buffer A (10 mM Tris-HCl [pH 8.2], 100 mM NaCl, and 20% sucrose) followed by the addition of 0.1 mL of ice-cold buffer B (100 mM Tris-HCl [pH 8.2], 50 mM EDTA, 0.6 mg/mL lysozyme) and incubated for 1 minute on ice. Then, 0.5 mL of Buffer C (10 mM Tris-HCl [pH 8.2], 10 mM EDTA, 10 mM spermidine, 1% Brij-58, and 0.4% deoxycholate) was added and the mixture incubated for 4 min at 5 C. The lysed cells were loaded onto linear sucrose density gradients (15%–50%) containing 10 mM Tris-HCl (pH 8.2) and 100 mM NaCl, and centrifuged for 20 minutes at 10,000 rpm in a Beckmann SW 41 Ti rotor. Fractions (1 mL each) were collected by piercing a hole at the bottom of the ultracentrifugation tube and DNA content quantified by DAPI (Invitrogen) fluorescence using a plate reader (Spectramax). The fractions were
Cell 175, 583–597.e1–e16, October 4, 2018 e6
then incubated with 20 U of benzonase at room temperature for 30 min. Fractions were run on a 15% Tris-HCl gel (BioRad) and submitted for tandem mass spectrometry analysis at the Koch Institute Biopolymers and Proteomics core facility. The region of the gel corresponding to proteins in the range of 5-25 kDa was isolated and cut into two fragments for mass spectrometry analysis. Mass spectra were analyzed using Mascot (Matrix Science, London, UK), searching against the Sprot_extra_022412 database assuming digestion with trypsin. To compare chromosomal supercoiling, nucleoids were isolated as above and loaded on to 10%–30% linear sucrose gradients containing varying amounts of ethidium bromide. For novobiocin treated nucleoids, cells were treated with 25 mg/mL novobiocin (Sigma) for 15 min before cell harvest. Intercalation of ethidium bromide into DNA first causes relaxation of the chromosome and decreased sedimentation and then at higher concentrations of ethidium increased sedimentation due to introduction of (+) supercoils. The amount of ethidium bromide at the point of slowest migration is reflective of the concentration of negative supercoils in the chromosome (Worcel and Burgi, 1972). Gradients were centrifuged at 10,000 rpm in a Beckman SW41 for 15 min. The distance sedimented by the nucleoid from the top of the gradient was measured. Sedimentation data from multiple runs with different ethidium bromide concentrations were pooled for analysis. The concentration of ethidium required at the point of slowest migration was estimated by fitting the data with a LOWESS fit using Prism 7 software. Expression profiling with microarrays and RNA-sequencing RNA was extracted using hot trizol lysis and the Direct-zol RNA MiniPrep (Zymo). rRNAs were removed using the Ribo-Zero Kit for Gram-negative bacteria (Illumina). The rRNA-depleted RNA was then fragmented using the RNA fragmentation reagents (Ambion) at 70 C for 8 min. Fragmented RNA was recovered by ethanol precipitation and resuspended in 6 mL water. RNA-seq libraries were prepared based on the previously described dUTP protocol as follows. First-strand cDNA synthesis was performed by addition of 1 mL of 3 mg/mL random primers (Invitrogen) followed by incubation at 5 min at 65 C. Samples were then placed on ice for 1 min and the following reagents were added: 4 mL first-strand synthesis buffer, 2 mL 100 mM DTT, 1 mL 10 mM dNTPs, 1 mL Superasein (ThermoFisher) and 4 mL water. Reactions were incubated at room temperature for 2 min and 1 mL Superscript III (ThermoFisher) was added and the following thermocycler program was used: 10 min at 25 C, 1hr at 50 C and 15 min at 70 C. The reaction volumes were brought to 200 mL and a phenol:chloroform:isoamyl alcohol (25:24:1) (Sigma) extraction was performed. cDNA was resuspended in 104 mL water and second-strand synthesis reactions were set up with 30 mL of 5X second-strand synthesis buffer (ThermoFisher), 4 mL of 100 mM dNTPs with dUTP (Promega) instead of dTTP, 4 mL of 5X first-strand synthesis buffer and 2 mL 100 mM DTT. Reactions were incubated on ice for 5 min and 1 mL RNase H (NEB), 1 mL E.coli DNA ligase (NEB) and 4 mL DNA Pol I (NEB) were added, followed by an incubation at 16 C for 2.5 hr. cDNA was recovered by Ampure XP (Beckman Coulter) bead purification with 100 mL beads in 450 mL 20% PEG/NaCl solution. Sequencing libraries were built by first end repairing cDNA with 5 mL T4 DNA polymerase (NEB), 5 mL T4 PNK (NEB), and 1 mL Klenow large fragment (NEB) in 100 mL T4 DNA ligase buffer with 0.25 mM dNTPs for 30 min at room temperature. Repaired DNA was recovered by Ampure XP (Beckman Coulter) bead purification, using 100 mL beads in 300 mL 20% PEG/NaCl solution. Beads were washed twice with 80% ethanol, dried, and resuspended in 32 mL EB. The bead slurry was directly treated with 3 mL Klenow (30 /50 exo-) (NEB) in 50 mL NEB Buffer #2 with 0.2 mM ATP at 37 C for 30 min to add 30 overhangs to DNA. The reaction was cleaned up by addition of 150 mL 20% PEG/NaCl and capture of the Ampure XP beads. After two washes with 80% ethanol, the beads were resuspended in 23 mL EB and the supernatant was transferred to a clean tube. Y-shaped Illumina adapters were ligated onto the DNA in 50 mL total volume in Quick Ligase buffer with 3 mM Y-shaped adaptors and 1 mL Quick Ligase (NEB) for 15 min at room temperature. Y-shaped adaptors were prepared by annealing Illumina PE adaptor 1 and Illumina PE adaptor 2. 75 mL 20% PEG/NaCl was added, and the DNA then recovered by addition of Ampure beads followed by two ethanol washes. After resuspension in 23 mL water, the Ampure beads were captured and the DNA in the supernatant was transferred to a clean tube. Digestion of the second strand was performed by addition of 6 mL of 5x HF Phusion buffer and 1 mL USER (NEB) and incubation at 37 C for 15 min, followed by 5 min at 95 C. DNA libraries were amplified in 80 mL final volume with Phusion DNA polymerase (NEB) in Phusion High GC buffer supplemented with Betaine (0.4 M final) (Sigma). The total number of cycles was optimized for each sample such that 10-14 cycles of PCR was used for each sequencing library. Libraries were purified by addition of 240 mL 20% PEG/NaCl and capture of the Ampure XP beads. After two washes with 80% ethanol, the beads were resuspended in 20 mL EB. Elutions were then run on a 8% TBE polyacrylamide gel (ThermoFisher) for 30 min at 180 V and a 230-500 bps region was gel extracted using a Spin-X 0.22 mm cellulose acetate column (Costar). Single-end sequencing was performed on an Illumina Hi-Seq 2000 at the MIT Bio Micro Center. For DNA microarray experiments, 10 mL cultures of ML76 were grown in PYE + 0.3% xylose, M2X or M2G to mid-exponential phase and treated with 25 mg/mL novobiocin (Sigma) for 15 min or left untreated to use as a reference for the microarray. RNA was extracted using the RNeasy Mini Kit (QIAGEN). The generation of labeled cDNA was performed using 20 mg of RNA in 14 mL RNase-free water to which 1 mL of random primers (Invitrogen) was added. The samples were incubated at 65 C for 10 min, followed by 2 min incubation on ice. 9.6 mL of trimix (6 mL 5X first-strand buffer, 3 mL 0.1 M DTT, and 0.6 mL dNTP mix [25 mM A,G,T, 10 mM C]), 1 mL of Superscript II (ThermoFisher), and 3 mL of labeled dCTP (Cy5 dye for experimental and Cy3 dye for reference) (Sigma) were added and incubated covered with foil for 10 min at room temperature. Samples were then incubated at 42 C for 1 hr, after which an additional 1 mL of Superscript II was added for another hour. RNA was degraded by addition of 1.5 mL of 0.5 N NaOH and cDNA purified by addition of 1.5 mL of 0.5 N HCl followed by purification with a PCR Purification Kit (QIAGEN). After a wash with PE buffer, cDNA was eluted in 50 mL EB. 1 mL of 20 mg/mL yeast tRNA (Sigma) was added and the total volume was adjusted to 55 mL. 55 mL e7 Cell 175, 583–597.e1–e16, October 4, 2018
of 2X-HiRPM hybridization buffer (Agilent) as added to the samples and incubated for 3 min at 95 C followed by 1 min at 42 C. 100 mL was loaded onto a custom Agilent array placed in a hybridization chamber and incubated overnight rotating at 65 C. Arrays were washed for 5 min in Oligo aCGH/Chip-on-chip wash buffer 1 (Agilent), washed for 5 min in Oligo aCGH/Chip-on-chip wash buffer 2 (Agilent) at 30 C, and washed with acetonitrile for 2 min at room temperature. Arrays were dried Stabilization and Drying Solution (Agilent) and imaged using an Agilent scanner at the MIT Bio Micro Center. Chromatin Immunoprecipitation Sequencing (ChIP-seq) ChIP was performed using strain ML2794 (gapR::gapR-3xFLAG) or wild-type CB15N grown in PYE. The ChIP procedure was as follows. Briefly, cell cultures (20 mL) were grown to OD600 0.3 and fixed by the addition of 10 mM sodium phosphate [pH 7.6] and 1% formaldehyde (final concentrations) (Sigma). When appropriate, 25 mg/mL of rifampicin (Sigma) was added to cells for 20 minutes prior to fixation. Fixed cells were incubated at room temperature for 10 minutes and then quenched with 0.1 M glycine (Sigma) for 5 min at room temperature followed by 15 min on ice. Cells were washed three times with 1X PBS [pH 7.4] and resuspended in 500 mL of TES buffer (10 mM Tris-HCl [pH 7.5], 1 mM EDTA, 100 mM NaCl) to which 35,000 U of Ready-Lyse (Epicenter) was added. Following 15 min incubation at RT, 500 mL of ChIP buffer (16.7 mM Tris-HCl [pH 8.1], 167 mM NaCl, 1.1% Triton X-100, 1.2 mM EDTA) containing protease inhibitors (Roche cOmplete EDTA-free tablets) was added. After 10 min at 37 C, the lysates were sonicated on ice and cell debris cleared by centrifugation. Supernatant protein concentration was measured by Bradford assay (Thermo Scientific) and 500 mg of proteins were diluted into 1 mL of ChIP buffer + 0.01% SDS. The diluted supernatants were pre-cleared for 1 hr at 4 C on a rotator with 50 mL of Protein-A Dynabeads (Life Technologies) pre-blocked overnight in ChIP buffer + 0.01% SDS and 100 mg ultrapure BSA (Ambion). Beads were pelleted and 90 mL of the supernatant was removed as input DNA and stored at 80 C, the remaining pre-cleared supernatant was incubated rotating at 4 C overnight with 1 mL of FLAG antibody (Sigma). The immune complexes were captured for 2 hr at 4 C with 50 mL of pre-blocked Protein-A Dynabeads. Beads were then washed consecutively at 4 C for 15 min with 1 mL of the following buffers: low salt wash buffer (0.1% SDS, 1% Triton X-100, 2 mM EDTA, 20 mM Tris-HCl [pH 8.1], 150 mM NaCl), high salt wash buffer (0.1% SDS, 1% Triton X-100, 2 mM EDTA, 20 mM Tris-HCl [pH 8.1], 500 mM NaCl), LiCl wash buffer (0.25 M LiCl, 1% NP-40, 1% deoxycholate, 1 mM EDTA, 10 mM Tris-HCl [pH 8.1] and twice with TE buffer (10 mM Tris-HCl [pH 8.1], 1 mM EDTA). Complexes were then eluted twice from the beads with 250 mL of freshly prepared elution buffer (1% SDS, 0.1 M NaHCO3). To reverse crosslinking, 300 mM of NaCl and 2 mL of RNase A (0.5 mg / mL) (QIAGEN) were added to the collective eluates which were incubated at 65 C overnight. Samples were then incubated at 45 C for 2 hr with 5 mL of Proteinase K (NEB) in the presence of 40 mM EDTA (pH 8.0) and 40 mM Tris-HCl [pH 6.8]. DNA from the samples was then extracted twice with phenol:chloroform:isoamyl alcohol (25:24:1) (Sigma) and subsequently precipitated by adding sodium acetate (pH 5.2), 100 mg glycogen (ThermoFisher) and 1 volume of ice cold isopropanol, and stored at 20 C overnight. DNA was pelleted and washed with 75% ethanol and resuspended in TE buffer [pH 8.0]. Sequencing libraries were built by first end repairing cDNA with 5 mL T4 DNA polymerase (NEB), 5 mL T4 PNK (NEB), and 1 mL Klenow large fragment (NEB) in 100 mL T4 DNA ligase buffer with 0.25 mM dNTPs for 30 min at room temperature. Repaired DNA was recovered by Ampure XP (Beckman Coulter) bead purification, using 100 mL beads in 300 mL 20% PEG/NaCl solution. Beads were washed twice with 80% ethanol, dried, and resuspended in 32 mL EB. The bead slurry was directly treated with 3 mL Klenow (30 /50 exo-) (NEB) in 50 mL NEB Buffer #2 with 0.2 mM ATP at 37 C for 30 min to add 30 overhangs to DNA. The reaction was cleaned up by addition of 150 mL 20% PEG/NaCl and capture of the Ampure XP beads. After two washes with 80% ethanol, the beads were resuspended in 23 mL EB and the supernatant was transferred to a clean tube. Y-shaped Illumina adapters were ligated onto the DNA in 50 mL total volume in Quick Ligase buffer with 3 mM Y-shaped adaptors and 1 mL Quick Ligase (NEB) for 15 min at room temperature. Y-shaped adaptors were prepared by annealing Illumina PE adaptor 1 and Illumina PE adaptor 2. 75 mL 20% PEG/NaCl was added, and the DNA then recovered by addition of Ampure beads followed by two ethanol washes. After resuspension in 30 mL water, the Ampure beads were captured and the DNA in the supernatant was transferred to a clean tube. DNA libraries were amplified in 80 mL final volume with Phusion DNA polymerase (NEB) in Phusion High GC buffer supplemented with Betaine (0.4 M final) (Sigma). The total number of cycles was optimized for each sample such that 10-14 cycles of PCR was used for each sequencing library. Libraries were purified by addition of 240 mL 20% PEG/NaCl and capture of the Ampure XP beads. After two washes with 80% ethanol, the beads were resuspended in 20 mL EB. Elutions were then run on a 8% TBE polyacrylamide gel (ThermoFisher) for 30 min at 180 V and the 250-500 bps region was gel extracted using a Spin-X 0.22 mm cellulose acetate column (Costar). Single-end sequencing was performed on an Illumina Hi-Seq 2000 at the MIT Bio Micro Center. RpoC-3xFLAG ChIP profiles were from (GSE73925). Caulobacter synchronizations For synchronizations, Caulobacter strains were grown to mid-exponential phase and G1/swarmer cells were isolated using Percoll (GE Healthcare) density gradient centrifugation. Cells were pelleted at 8000 x g for 5 min and resuspended in M2 salts (0.87 g/L Na2HPO4, 0.53 g/L KH2PO4, 0.5 g/L NH4Cl) with an equal volume of Percoll. The suspension was centrifuged at 10000 x g for 20 min and G1 cells (bottom layer) were isolated. The G1 cells were then washed with 5 3 1 mL M2 salts before resuspension and release into appropriate medium.
Cell 175, 583–597.e1–e16, October 4, 2018 e8
Immunoblotting Nucleoid fractions were electrophoresed with 4%–20% Tris-HCl gels (Bio-Rad) and either stained with Coomassie or transferred onto PVDF and immunoblotted. For whole cell immunoblots, cells were harvested and their proteins precipitated with 20% TCA, washed with acetone, and normalized by OD for resuspension. Proteins were electrophoresed on 12% Bis-Tris gels in MES buffer (Invitrogen) with at least 2 replicates for each condition. For estimates of GapR copy per cell, a standard curve was generated with purified GapR protein. Antibodies were used at the following concentrations: anti-RpoA (1:7500, Neoclone), anti-FLAG (1:7500, Sigma), anti-IHF (1:1000) (generous gift of S. Goodman), anti-CtrA (1:5000), and anti-GapR (1:1000). Rabbit polyclonal anti-GapR antibody (Covance) was generated from purified His6-Thrombin-GapR and was used without purification at 1:1000. HRP-conjugated secondary antibodies (ThermoFisher) were used at the following concentrations: goat anti-rabbit, 1:4000; goat anti-mouse, 1:7500. The membranes were developed with SuperSignal West Femto Maximum Sensitivity Substrate (ThermoFisher) and visualized with a FluorChem R Imager (ProteinSimple). Microscopy Phase contrast images were taken on a Zeiss Observer Z1 microscope using a 100x/1.4 oil immersion objective and an LED-based Colibri illumination system using software Metamorph (Universal Imaging, PA). Cells were diluted to maintain OD600 < 0.4. For phase contrast images, cells were fixed with 0.5% paraformaldehyde (Sigma), pelleted, and resuspended in PBS. Fixed cells were spotted onto PBS 1.5% agarose pads and imaged. Flow cytometry Flow cytometry was performed as follows. 300 mL cells were fixed with 700 mL 70% ethanol, spun at 5000 x g for 2.5 min, resuspended in 50 mM Na2CO3 with 5 mg/mL RNase A (QIAGEN), and incubated at 50 C for at least 4 hr. Cells were diluted 1:10 into Na2CO3 buffer containing 0.5 mL/mL SYTOX Green nucleic acid stain (Invitrogen) and analyzed on a BD Accuri C6 flow cytometer. For time course experiments, synchronized swarmer cells were released into fresh media and sampled every 20 min after release. For novobiocin experiments, novobiocin (Sigma) was added to 2.5 mg / mL 20 min after release. DNA sequencing DNA sequencing was performed on 500 mL of each Caulobacter strain grown to mid-exponential phase (OD600 0.2-0.4) in PYE + 0.3% xylose. For GapR depletion, cells were grown in 0.3% xylose and either synchronized or harvested, washed 4 x with PYE, and resuspended in 0.2% glucose for 2 hr before synchronization. The wild-type control was grown in 0.2% glucose before synchronization. After synchronization, cells were released into 25 mL of PYE + glucose at OD600 0.2. At each time-point after release (0, 20, 40, 60, 80 min), 2 mL of cells were centrifuged at 14,000 rpm for 1 min, aspirated, and frozen in liquid nitrogen. DNA was isolated by CTAB purification. Pellets were resuspended in TE (10 mM Tris-HCl [pH 8.0], 1 mM EDTA) and normalized so that 1 mL of culture with OD600 0.1 was present per 10 mL resuspension. 400U of Ready-Lyse (Epicenter) at a final concentration of 0.125% was added to 50 mL of resuspension, and incubated for 5 min at room temperature. The volume was brought up to 800 mL in TE supplemented with 10 mM EDTA, 0.75% Triton X-100, and 2.5 mL of 20 mg/mL Proteinase K (NEB), and incubated at 50-65 C overnight. 100 mL 5 M NaCl and 100 mL of CTAB (Sigma) was added, and samples then incubated at 65 C for 10 min. The DNA in each sample was then extracted by phenol:chloroform:isoamyl alcohol (25:24:1) (Sigma), and precipitated with isopropanol. DNA was resuspended in 10 mM Tris-HCl [pH 8.0], and sheared to 300-800 bp fragments using a Bioruptor Plus sonicator (Diagenode). RNase A (QIAGEN) was added and the fragments were electrophoresed, and 300-600 bp fragments were size selected, purified with Gel Extraction Kit (QIAGEN), and resuspended in 50 mL EB. Sequencing libraries were built by first end repairing fragmented DNA with 4 mL T4 DNA polymerase (NEB), 4 mL T4 PNK (NEB), and 0.75 mL Klenow large fragment (NEB) in 100 mL T4 DNA ligase buffer with 0.25 mM dNTPs for 30 min at room temperature. Repaired DNA was recovered by Ampure XP (Beckman Coulter) bead purification, using 100 mL beads in 300 mL 20% PEG/NaCl solution. Beads were washed twice with 80% ethanol, dried, and resuspended in 30 mL EB. The bead slurry was directly treated with 3 mL Klenow (30 /50 exo-) (NEB) in 40 mL NEB Buffer #2 with 0.2 mM ATP at 37 C for 45 min to add 30 overhangs to DNA. The reaction was cleaned up by addition of 120 mL 20% PEG/NaCl and capture of the Ampure XP beads. After two washes with 80% ethanol, the beads were resuspended in 30 mL EB and Y-shaped Illumina adapters were ligated onto the DNA in 50 mL total volume in T4 Ligase buffer with 3 mM Y-shaped adaptors and 2 mL T4 DNA ligase (NEB) for 30 min at room temperature. Y-shaped adaptors were prepared by annealing Illumina PE adaptor 1 and Illumina PE adaptor 2. 150 mL 20% PEG/NaCl was added, and the DNA then recovered by capture of the Ampure beads followed by two ethanol washes. After resuspension in 20 mL water, the Ampure beads were captured and the DNA in the supernatant was transferred to a clean tube. DNA libraries were amplified with Kapa HiFi DNA polymerase (Kapa Biosystems) and multiplexed with Illumina PE PCR primer 1 and Illumina PE PCR primer 2 (multiplexing) in Kapa High GC buffer supplemented with DMSO. The total number of cycles was optimized for each sample such that 10-14 cycles of PCR was used for each sequencing library. DNA was paired-end sequenced on a NextSeq500 at the MIT Bio Micro Center. qPCR analysis of origin-terminus ratio Cells were grown in PYE with xylose, washed, and resuspended in PYE with glucose for 2 hr before synchronization and release into fresh glucose. 1 mL samples were harvested in 30 min increments and flash frozen. DNA was extracted by resuspending pellets in
e9 Cell 175, 583–597.e1–e16, October 4, 2018
600 mL Cell Lysis Solution (QIAGEN) and incubating at 80 C for 5 min to lyse cells. RNAs were removed by treatment with 50 mg RNase A (QIAGEN) at 37 C for 30 min. 200 mL Protein Precipitation Solution (QIAGEN) was added, the sample vortexed, and left on ice for 30 min to precipitate proteins. After spinning at 14,000 rpm for 5 min, the supernatant was transferred to a tube containing 600 mL isopropanol and mixed by inversion. DNA was harvested by spinning at 14,000 rpm for 1 min followed by a wash with 600 mL of 70% ethanol. The DNA pellet was resuspended in 20 mL H2O. For qPCR, DNAs were diluted 1:10 or 1:20 and DNAs was mixed with either Cori_2 or CCNA_01869 forward/reverse primer mix and 2X qPCR Master Mix. All experimental samples and standard curves were loaded onto a 384-well plate in triplicate for qPCR. qPCR was conducted in a LightCycler 480 system (Roche) using the following thermocycler program: 95 C/10 minutes, 95 C/15 s, 60 C/30 s, 72 C/30 s with 40 cycles of steps 2-4. Cp values were calculated from LightCycler 480 software at the second derivative maximum. Technical replicates were averaged to yield a final Cp value for each sample and standard curve point. Plating viability and drug sensitivity Strains were grown to mid-exponential phase and diluted in PYE to OD600 0.2. Strains were then 10-fold serially diluted and 2 mL of each dilution spotted onto plates containing 0.3% xylose and one of the following compounds when appropriate: novobiocin (0.1 mg/mL) (Sigma), mitomycin C (0.1 mg/mL) (Sigma), chloramphenicol (0.1 mg/mL) (Sigma), hydroxyurea (0.3 mg/mL) (Sigma), or rifampicin (0.03 mg/mL) (Sigma). Plates were incubated at 30 C for two or three days and imaged with a scanner. For experiments with minimal medium, strains were first grown in PYE to OD600 0.2, diluted into M2 + 0.2% glucose for wild-type cells or M2 + 0.3% xylose for the GapR depletion strain and grown to OD600 0.2, and then spotted onto M2 agar supplemented with 0.2% glucose or 0.3% xylose and grown for three days before imaging. Phylogeny analysis Phylogeny analysis was performed using the STRING database. Alignment logo was generated using WebLogo Version 2.8. Protein alignments were generated with JalView using Muscle with default settings. Genomes shown: Caulobacter crescentus (NA1000), Caulobacter K31, Brevundimonas subvibrioides, Caulobacter phage Cr30, Agrobacterium tumefaciens str. Cherry 2E-2-2, Rhodopseudomonas palustris, Sinorhizobium meliloti (strain SM11), Sphingomonas sp. (strain SKA58), Azospirillum phage Cd, Acetobacter pasteurianus, Wolbachia endosymbiont of Muscidifurax uniraptor, Bosea sp. Root 381, and Rickettsia monacensis. Protein expression and purification Purified GcrA and glucose-6-phosphate dehydrogenase (G6PD) were part of the lab collection and a gift from B. Wang in the Laub lab, respectively. Protein expressions was performed in BL21(DE3) cells grown in 2xYT to OD600 0.4-0.5 at 37 C and induced with 0.4 mM IPTG for 18-20 hr at 18 C unless otherwise noted. When appropriate, SUMO protease was added and the cleavage reaction adjusted to 400 mM NaCl. Cleavage proceeded for 2 hr or overnight at 4 C. His6-GapR and His7-SUMO-GapR: 1 L of cells was harvested by centrifugation and resuspended in 20 mL buffer A (50 mM sodium phosphate [pH 7.5], 0.5 M NaCl, 20 mM imidazole) supplemented with a SIGMAFAST Protease Inhibitor Tablet (Sigma). The cell resuspension was treated with 20 mg lysozyme (Sigma), 1 mM PMSF (Sigma), and 300 U benzonase (Novagen) at room temperature with gentle agitation for 20 min and then lysed by sonication. The cell debris was cleared by centrifugation for 30 min at 10,000 x g and passed over Ni-NTA agarose resin (QIAGEN) pre-equilibrated with buffer A at 4 C. Resin was washed with buffer A and proteins were eluted stepwise with 2.5 column volumes of buffer A containing increasing amounts of imidazole (50-100-150-250-500 mM imidazole). Fractions containing GapR (150-500 mM) were combined. For His7-SUMO-GapR, SUMO cleavage proceeded overnight at 4 C, supplemented with 300 U of Benzonase (Novagen) and 10 mM final concentration of MgCl2 to ensure removal of DNA during cleavage. Eluates were diluted or buffer exchanged with PD-10 columns (GE Healthcare) to have < 100 mM NaCl and run on a HiTrap Heparin HP (5mL) (GE Healthcare) column, pre-equilibrated in buffer B (50 mM sodium phosphate [pH 7.5], 10 mM NaCl). GapR and His6-GapR were eluted with a 2-step elution (buffer B + 0.5 M NaCl and buffer B + 1.0 M NaCl), each step being 5 column volumes. High salt fractions containing GapR were concentrated on an Amicon 3K Centrifugal Filter Unit (Millipore) and buffer exchanged into storage buffer (50 mM sodium phosphate pH 7.5, 200 mM NaCl, 10% (v/v) glycerol) by gel filtration using a Superdex 200 Increase 10/300 GL column (GE Healthcare). Fractions containing GapR were identified by SDS-PAGE/Coomassie staining. Purification of His7-SUMO-ParC and His7-SUMO-ParE was performed as follows. 1 L of cells was induced with 0.4 mM IPTG for 4 hr at 37 C and harvested and resuspended in 20 mL of buffer A800 (20 mM Tris-HCl [pH 8.0], 30 mM imidazole, 800 mM NaCl, 2 mM ßME, 10% (v/v) glycerol) supplemented with a SIGMAFAST Protease Inhibitor Tablet (Sigma). After lysis and sonication, the lysate was applied to Ni-NTA agarose resin (QIAGEN), washed with 10 column volumes of buffer A800, and then washed with 20 column volumes with buffer A400 (buffer A800 with 400 mM NaCl). Proteins were eluted with B400 (buffer A400 with 500 mM imidazole), and cleaved with SUMO protease. After validating SUMO cleavage by SDS-PAGE/Coomassie staining, the mixture was concentrated on an Amicon 50K Centrifugal Filter Unit (Millipore) and applied to a Superdex 200 Increase 10/300 GL column (GE Healthcare) equilibrated in 20 mM Tris-HCl [pH 8.0], 400 mM KCl, 2 mM ßME, and 10% (v/v) glycerol. Fractions containing ParC or ParE identified by SDS-PAGE/Coomassie staining were concentrated with a Amicon 50K Centrifugal Filter Unit. Glycerol was added and the buffer adjusted to make a final storage buffer of 20 mM Tris-HCl [pH 8.0], 400 mM KCl, 2 mM ßME, 30% (v/v) glycerol and protein aliquots were stored at 80 C.
Cell 175, 583–597.e1–e16, October 4, 2018 e10
His7-SUMO-GyrA and His7-SUMO-GyrB were purified as follows. For GyrA, 1 L of cells was grown in 2xYT at 37 C to OD600 0.5 and induced with 0.4 mM IPTG for 4 hr at 37 C and resuspended in 20 mL of buffer A1000 (50 mM HEPES, [pH 7.5], 20 mM imidazole, 1000 mM NaCl, 2 mM ßME, 10% (v/v) glycerol) supplemented with a SIGMAFAST Protease Inhibitor Tablets (Sigma). After lysis, sonication, and purification on Ni-NTA agarose resin (QIAGEN), His7-SUMO was cleaved off using SUMO protease. After validating SUMO cleavage by SDS-PAGE/Coomassie staining, the mixture was concentrated on an Amicon 50K Centrifugal Filter Unit (Millipore) and applied to a Superdex 200 Increase 10/300 GL column (GE Healthcare) equilibrated in 50 mM HEPES [pH 7.5], 400 mM KCl, 1 mM EDTA, 2 mM ßME, and 10% (v/v) glycerol. Fractions containing GyrA identified by SDS-PAGE/Coomassie staining were pooled and concentrated on a Amicon 50K Centrifugal Filter Unit. Glycerol was added to generate a final storage buffer of 45 mM HEPES [pH 7.5], 360 mM KCl, 0.9 mM EDTA, 1.8 mM ßME, 19% (v/v) glycerol and protein aliquots were stored at 80 C. For GyrB, 1 L of induced cells was resuspended in 20 mL of buffer A800 (20 mM HEPES [pH 7.5], 30 mM imidazole, 800 mM NaCl, 2 mM ßME, 10% (v/v) glycerol) supplemented with a SIGMAFAST Protease Inhibitor Tablets (Sigma). After lysis, sonication, and purification on Ni-NTA agarose resin (QIAGEN), His7-SUMO was cleaved off using SUMO protease. The mixture was concentrated on an Amicon 50K Centrifugal Filter Unit (Millipore) and applied to a Superdex 200 Increase 10/300 GL column (GE Healthcare) equilibrated in 50 mM HEPES [pH 7.5], 500 mM KCl, 1 mM EDTA, 2 mM ßME, and 10% (v/v) glycerol. Fractions containing GyrB identified by SDS-PAGE/Coomassie staining were pooled, concentrated on a Amicon 50K Centrifugal Filter Unit. Glycerol was added and the buffer adjusted to make a final storage buffer of 50 mM HEPES [pH 7.5], 500 mM KCl, 1 mM EDTA, 2 mM ßME, 30% (v/v) glycerol and protein aliquots were stored at 80 C. E. coli GyrA and GyrB were purified from pET28-GyrA and pET28-GyrB (generous gift from J. Berger) as follows. Proteins were grown in 2xYT at 37 C to OD600 0.5 and induced with 0.25 mM IPTG overnight. Pellets were resuspended in buffer A800 (20 mM Tris-HCl, [pH 8.0], 800 mM NaCl, 30 mM imidazole, 10% (v/v) glycerol) supplemented with a SIGMAFAST Protease Inhibitor Tablet (Sigma) and lysed, sonicated, and loaded on Ni-NTA agarose resin (QIAGEN). GyrA was washed with buffer B (buffer A1000 except with 100 mM NaCl), and eluted in 20 mM Tris-HCl, [pH 8.0], 100 mM NaCl, 500 mM imidazole, 2 mM ßME, and 10% (v/v) glycerol. Proteins were then washed back into buffer B with 2 mM ßME with a Amicon 50K Centrifugal Filter Unit (Millipore). 1:50 (w/w) TEV protease (Sigma) was added and proteins were incubated overnight at 4 C with shaking. Cleaved GyrA was recovered by passage through Ni-NTA agarose resin (QIAGEN), concentrated, and applied to a Superdex 200 Increase 10/300 GL column (GE Healthcare) equilibrated in 50 mM Tris-HCl, [pH 8.0], 500 mM KCl, 2 mM ßME, and 10% (v/v) glycerol. Fractions containing GyrA identified by SDS-PAGE/Coomassie staining were pooled, concentrated, and snap frozen. His-TEV-GyrB was eluted from Ni-NTA resin with buffer C (buffer A1000 except with 400 mM NaCl), concentrated, and directly applied to a Superdex 200 Increase 10/300 GL column (GE Healthcare) equilibrated in 50 mM HEPES [pH 7.5], 500 mM KCl, 1 mM EDTA, 2 mM ßME, and 10% (v/v) glycerol. Fractions containing His-TEV-GyrB identified by SDS-PAGE/Coomassie staining were pooled and concentrated. Glycerol was added and the buffer adjusted to make a final storage buffer of 50 mM HEPES [pH 7.5], 500 mM KCl, 1 mM EDTA, 2 mM ßME, 30% (v/v) glycerol and protein aliquots were stored at 80 C. Archaeoglobus fulgidus reverse gyrase (Rgy) was purified as follows: pSG483 (generous gift from J. Berger) was transformed into E. coli C41(DE3). Cells were grown in 2xYT to A600 = 0.5 at 37 C and induced with 0.5 mM IPTG for 18-20 h at 18 C. 1 L of cells was resuspended in 20 mL of buffer A (50 mM Tris-Cl [pH 8.0], 300 mM NaCl) supplemented with a SIGMAFAST Protease Inhibitor Tablet (Sigma). The cell resuspension was treated with 20 mg lysozyme and 1 mM PMSF at room temperature with gentle agitation for 20 min and then lysed by sonication. Lysates were spun for 30 min at 40,000 rpm in a Beckman Ti-45 rotor and the supernatant was passed over Ni-NTA agarose resin (QIAGEN) pre-equilibrated with buffer A + 10mM imidazole. The resin was washed with buffer A + 20 mM imidazole and eluted in buffer A + 150 mM imidazole. The eluate was concentrated on Amicon 50K Centrifugal filter unit (Millipore) and supplemented with MgCl2 to a final concentration of 10 mM. The concentrate was applied to a Superdex 200 Increase 10/300 GL column equilibrated in buffer B (20mM Tris-Cl [pH 8.0], 200 mM NaCl, 10 mM MgCl2). Fractions containing Rgy as determined by SDS-PAGE and Coomassie staining were combined and concentrated on Amicon 50K Centrifugal filter unit (Millipore), aliquoted, and stored at 80 C. For C. crescentus and Bosea sp. GapR structural studies, pET15b expression plasmids were transformed into C41(DE3) cells. For protein expression, cells containing the plasmids were grown to an OD600 of 0.5 at 37 C and induced with 0.5 mM IPTG overnight at 15 C. Both proteins were purified using the same protocol. Briefly, the cells containing expressed GapR protein were pelleted and then resuspended in Buffer A (50 mM Tris-Cl [pH 7.5], 300 mM NaCl, 5% glycerol) with 1.0 mM b-mercaptoethanol (BME) and 600 U benzonase. The reconstituted cells were lysed using a micro-fluidizer and then centrifuged at 16000 rpm for 40 min at 4 C. The lysate was loaded onto a cobalt-NTA column and washed overnight with Buffer A. The following day, the protein was eluted stepwise, with increasing concentrations of imidazole. Both GapR proteins eluted within the 100-500 mM imidazole range and were > 95% pure at this stage. The proteins were further purified via a size exclusion chromatography step. To produce selenomethionine-substituted proteins for phasing, a C. crescentus GapR(L83M) mutant was generated as the WT protein has only one methionine. The Bosea sp. protein has no methionines except the N-terminal methionine so a Bosea sp. GapR(L48M-I54M-L73M) mutant was constructed. The selenomethionine proteins were purified as per the WT proteins except that 5 mM BME was used in the purification buffers. After purification, the proteins were concentrated to 1 mg/mL, buffer exchanged in buffer A without imidazole and the his-tag was cleaved overnight at RT using a thrombin capture cleavage kit. The cleaved tag was removed by buffer exchanging the protein into buffer A using a 10 kDa MW cutoff centricon.
e11 Cell 175, 583–597.e1–e16, October 4, 2018
SEC-MALS SEC-MALS was performed with assistance from B. Brown on equipment owned by the T. Baker Lab at MIT. Size-exclusion chromatography was performed on a Wyatt WTC-030S5 size-exclusion column with an Agilent HPLC. Multi-angle light scattering was measured with a Wyatt DAWNHELEOS-2 instrument and protein concentration was determined using a Wyatt Optilab rEX. The system was equilibrated in GapR storage buffer (50 mM sodium phosphate [pH 7.5], 200 mM NaCl, 10% (v/v) glycerol). 100 mL of 1 mg/mL GapR was injected and run at 0.5 mL/min in GapR storage buffer at room temperature. The signals were aligned and normalized internally with the ASTRA software 5.3.4 (Wyatt Technology). Bovine serum albumin (Sigma) was used as a standard. Absorbance at 280 nm was used to calculate protein molecular weight. Preparation of DNA substrates For 50, 20, and 16 bp DNA substrates, oligos were denatured at 95 C for 2 min, and allowed to cool slowly to room temperature. Larger DNA substrates were generated by PCR and purified with PCR Purification Kits (QIAGEN). Negatively supercoiled pUC19 was harvested from E. coli using a Midiprep Kit (QIAGEN) or purchased from NEB. Relaxed pUC19 was generated with E. coli Topoisomerase I (NEB). Nicked pUC19 was generated with Nt.BspQI (NEB). Linear pUC19 was generated with SmaI-HF (NEB). Positively supercoiled pUC19 was made using reverse gyrase: a 50 mL reaction containing 16.33 mg of negatively supercoiled pUC19 and 500 nM reverse gyrase was incubated in 50 mM Tris-HCl [pH 8.0], 50 mM NaCl, 10 mM MgCl2, 1 mM ATP, and 10% PEG 8000 at 80 C for 30 min. All enzymatic reactions were stopped by adding 1% SDS and 10 mM EDTA (final concentration) and plasmid DNA recovered by phenol/chloroform extraction and isopropanol precipitation. DNA electrophoresis and visualization Short (< 300 bp) fragments were electrophoresed in 6% DNA Retardation gels (Invitrogen) at 100 V for 60-90 min at 4 C. pUC19 plasmids were electrophoresed in 1% TBE agarose gels either at 30 V for 16 hr at 4 C (binding assays) or at 130 V for 90 min (topology assays). Gels were stained in SYBR Gold (Invitrogen) and imaged with a Typhoon FLA 9500 imager (GE Lifesciences). For two-dimensional chloroquine gels, electrophoresis was performed by first running reactions on a 1% TBE agarose gel at 130 V for 90 min, then soaking the gel for 2 hr with shaking in 1X TBE supplemented with 4.5 mg/mL chloroquine phosphate (Santa Cruz Biotech). The gel was then turned and electrophoresed in the orthogonal direction at 130 V for 60 min in 1X TBE supplemented with 4.5 mg/mL chloroquine phosphate. Chloroquine is a DNA intercalator that introduces (+) supercoils. In chloroquine, (-) supercoiled plasmids will become more relaxed, and migrate more slowly, whereas (+) supercoiled DNA will be further compacted, increasing its migration speed. The gel was washed 4 3 20 min in distilled water to remove chloroquine, stained with SYBR Gold (Invitrogen), and imaged with a Typhoon FLA 9500 imager (GE Lifesciences). Electrophoretic mobility shift assays Reactions (15 mL) with indicated amounts of GapR and linear dsDNA or pUC19 DNA (40 ng) in binding buffer (40 mM sodium phosphate [pH 7.5], 100 mM NaCl, 50 mg/mL ultrapure BSA, 0.5 mM DTT) were incubated at 30 C for 30 min and then placed on ice. DNA loading buffer was added and 10 mL of the reactions were electrophoresed. For experiments with multiple DNA fragments, 40 ng total of DNA was used and the fragments were mixed at equal weight. DNA topology assays For topoisomerase I treatment, mixtures of GapR and relaxed pUC19 DNA (40 ng) in binding buffer (40 mM sodium phosphate [pH 7.5], 100 mM NaCl, 50 mg/mL ultrapure BSA, 0.5 mM DTT) were incubated at 30 C for 30 min. Calf Thymus Topoisomerase I (Invitrogen) was diluted to 1 U/mL in Topoisomerase Dilution Buffer (30 mM potassium phosphate [pH 7.0], 5 mM DTT, 0.1 mM EDTA, 0.2 mg/mL BSA, 50% (v/v) glycerol, 0.1% (w/v) Triton X-100) and 1 U was added to reactions and incubated for an additional 30 min at 30 C. Reactions were stopped by addition of 1% SDS and 10 mM EDTA (final concentration), and digested with Proteinase K (NEB) for 1 hr at 37 C. DNA loading buffer was then added and samples electrophoresed. For ligase assays, mixtures of GapR and nicked pUC19 DNA (40 ng) in 1X T4 DNA ligase buffer were incubated at 30 C for 30 min. T4 DNA ligase (NEB) was diluted to 1 U/mL in 1X T4 DNA ligase buffer and 1 U was added to reactions and incubated for an additional 2 hr at room temperature. Reactions were stopped as with topoisomerase I assays. DNA loading buffer was added and the reactions electrophoresed. GapR protein preparations additionally contained trace nuclease activity, which converted plasmids into their linear forms during the assay. The linear forms are observed as an intermediate band (see Figures 5F and S5E). DNase I protection assay Mixtures of GapR and negatively supercoiled pUC19 DNA (40 ng) in binding buffer (40 mM sodium phosphate [pH 7.5], 100 mM NaCl, 50 mg/mL ultrapure BSA, 0.5 mM DTT) were incubated at 30 C for 30 min. Turbo DNase (Invitrogen) was diluted to 1 U/mL in 1X DNase buffer and 1 mL was added to reactions and incubated for an additional 5 min at 30 C. Reactions were quenched by addition of 1 mL 10% SDS and incubation at 75 C for 10 min. Samples were subsequently analyzed by electrophoresis.
Cell 175, 583–597.e1–e16, October 4, 2018 e12
Crystallization of C. crescentus GapR-DNA complex, Bosea sp. apo and Bosea sp. GapR-DNA-bound complex To obtain crystals of the C. crescentus GapR-DNA complex, the protein was first concentrated to 20 mg/mL and 1 mM DNA (top strand 50 -TTAAAATTAAA-30 ), final concentration, was added. Excess DNA was removed by concentrating the sample with a 30 kDa MW cutoff microcon. Crystals were obtained at RT using the hanging drop vapor diffusion method by mixing the GapRDNA complex 1:1 with a crystallization solution consisting of 28% PEG 400, 0.1 M HEPES pH 7.5 and 0.2 M CaCl2 and suspending the drop over a reservoir containing the crystallization solution. The crystals grew overnight and reached their final size after 2 weeks. The crystals were cryo-preserved by direct transfer from the drop to the cryo-stream. WT GapR and selenomethionine GapR(L83M) proteins in complex with the DNA produced the same crystals. The apo Bosea sp. GapR was crystallized by hanging drop vapor diffusion at RT using protein at 40 mg/mL and mixing the protein 1:1 with a solution of 15% MPD, 10 mM MgCl2, 0.1 M HEPES pH 7.5 and 0.001 M spermine. The crystals grew overnight and could be cryo-preserved straight from the drop. The same crystals were produced with the selenomethionine-substituted Bosea sp. GapR(L48M-I54M-L73M) protein. Crystals of the WT Bosea sp. GapR bound to DNA (50 -TTAAAATTAAA-30 ) were grown via hanging drop vapor diffusion at RT by mixing the GapR protein (at 20 mg/mL) at a ratio of 1:1 (one GapR dimer-of-dimer to one DNA duplex) and combining the mixture with a crystallization solution consisting of 1.7 M ammonium sulfate and 0.1 M MES pH 6.5. The crystals reached their final size after 1-2 weeks and were cryopreserved by looping a crystal into a 2 uL drop consisting of the crystallization reagent supplemented with 15% glycerol for 1-2 s before direct placement of the crystal in the cryo-stream. Structure determination and refinement of the C. crescentus GapR-DNA complex, Bosea sp. apo and DNA-bound complex All X-ray intensity data were collected at the advanced light source (ALS) beamline 8.3.1. The selenomethionine C. crescentus GapR(L83M)-DNA complex was solved by single wavelength anomalous diffraction (SAD) using data collected from a crystal at the selenium peak. Heavy atom sites were located and density modification performed in Phenix Autosol. The density modified map from Autosol was used to construct the model in O (Jones et al., 1991). The crystal takes the space group C2 (Table S4) and there are four GapR subunits and one DNA duplex in the crystallographic asymmetric unit (ASU). Once the model was fit, it was subjected to refinement in Phenix. After several rounds of refinement, refitting and analysis with Molprobity (Chen et al., 2010), the model converged to Rwork/Rfree values of 24.1%/26.8% to 2.3 A˚ resolution. The final model contains residues 12-88 of three GapR subunits, 12-89 of the fourth GapR subunit, the entire DNA duplex and 104 water molecules. See Table S4 for final refinement statistics. Parameters of GapR-bound DNA was analyzed using the w3DNA program (Zheng et al., 2009). The apo Bosea sp. GapR(L48M-I54M-L73M structure) was solved by SAD using data collected from a crystal grown with selenomethionine protein. Autosol was used to locate the selenium sites and produce a density modified, electron density map. The crystal (P41212) contains a single GapR subunit in the ASU; the dimer is generated by crystallographic symmetry. The model was constructed in O (Jones et al., 1991). and refined using Phenix-refine. The final structure has Rwork/Rfree values of 22.2%/25.9% to 2.4 A˚ resolution and contains residues 9-78 of the GapR subunit and 3 solvent molecules. The final refinement statistics are provided in Table S4. The Bosea sp. GapR-DNA crystals were produced with the WT protein (non-selenomethionine) and the structure was solved by molecular replacement (MR). The Bosea sp. apo structure did not produce a solution in MR however, the C. crescentus GapR dimer-of-dimer with the DNA removed produced an MR solution using the program MolRep. The MR solution was first subjected to rigid body refinement and then x,y,z,b refinement using Phenix-refine to 2.0 A˚ resolution. After the first round of Phenix-refine (3 cycles) the electron density maps were used to fix the model and replace the C. crescentus GapR residues with those in the Bosea sp. GapR sequence. Molprobity was used throughout to assess the model during refinement (Chen et al., 2010). Density for the DNA was evident after initial refinement rounds however it could not be unambiguously traced. Ethidium bromide staining of the crystals also supported the presence of DNA. The final model includes 349 solvent molecules and four GapR subunits, with residues 6-78, 5-78, 3-78 and 3-80 and has final Rwork/Rfree values of 19.8%/22.2% to 2.0 A˚ resolution. Fluorescence polarization-based DNA binding experiments Fluorescence polarization (FP) experiments were performed using a PanVera Beacon 2000 FP system at 25 C. Samples were excited at 490 nM and fluorescence emission was measured at 530 nM. 50 -Fluoresceinated oligonucleotides used for DNA binding experiments were the site used for crystallization (AT-rich, top strand: 50 -F-TTAAAATTAAA-30 ), a GC-rich site (50 -FTGGGCGGGCGCCCGCCCA-30 ) and a mixed bp 12-mer (50 -F-GTGAGTACTCAC-30 ). For each FP experiment, WT C. crescentus GapR protein was titrated into 0.995 mL of reaction buffer (25 mM Tris, pH 7.5, 120 mM NaCl) containing 1 nM fluoresceinated oligonucleotide and points taken. All data were fit using Kaleigraph. The analyses were performed in triplicate (technical duplicates). Topo IV relaxation assays Caulobacter topo IV was reconstituted on ice and supercoiling assays were performed as follows. Briefly, equimolar purified ParC and ParE were mixed together on ice for 10 min to form 8.63 mM tetramer. Reconstituted topo IV was serially diluted in topo IV dilution buffer (50 mM Tris-HCl [pH 8.0], 150 mM potassium glutamate, 6 mM MgCl2, 10% (v/v) glycerol). 1.5 mL of 400 nM GapR and 1 mL of diluted topo IV tetramer were incubated together on ice for 30 min, then at 30 C for 5 min. 40 ng pUC19 DNA in topo IV reaction buffer (50 mM Tris-HCl [pH 8.0], 30 mM potassium glutamate, 6 mM MgCl2, 0.5% (v/v) glycerol, 12.5 mg/mL ultrapure BSA, 12.5 mM DTT, e13 Cell 175, 583–597.e1–e16, October 4, 2018
2 mM ATP) was added to generate 16 mL final volume (final concentration of topo IV was 17.55 nM for (-) supercoil relaxation, 2.34 nM for (+) supercoil relaxation). Topoisomerase concentrations and conditions were selected to deliberately yield relatively low topoisomerase activities. At the indicated times, reactions were quenched by addition of 1% SDS and 10 mM EDTA (final concentration), and treated with Proteinase K (NEB) for 1 hr at 37 C. Topoisomers were analyzed by electrophoresis. For reactions with G6PD or GcrA, assays proceeded with identical conditions as with GapR, using 400 nM (monomer) to maintain equimolar protein. For reactions with E. coli topo IV, topo IV (TopoGEN) was diluted in topo IV buffer and treated identically as with Caulobacter topo IV. 0.5 U total topo IV (0.03125 U per mL) was used for (-) supercoil relaxation and 0.25 U (0.0156 U per mL) was used for (+) supercoil relaxation. Gyrase supercoiling assays Caulobacter gyrase was reconstituted on ice and supercoiling assays performed as follows. Briefly, equimolar GyrA and GyrB were mixed on ice for 10 min to form 9.5 mM tetramer. Reconstituted tetramer was serially diluted in gyrase dilution buffer (25 mM Tris-HCl [pH 8.0], 150 mM potassium glutamate, 2.5 mM MgCl2, 10% (v/v) glycerol), and incubated for 5 min on ice after each dilution. 1.5 mL of 400 nM GapR and 1 mL of diluted gyrase tetramer were incubated on ice for 30 min, then at 30 C for 5 min. 40 ng of pUC19 DNA in gyrase reaction buffer (30 mM Tris-HCl [pH 8.0], 120 mM potassium glutamate, 5.5 mM MgCl2, 0.1 mg/mL ultrapure BSA, 1 mM DTT, 2 mM ATP, 10% (v/v) glycerol) was added to generate 16 mL final volume (final concentration of gyrase was 12.5 nM for introduction of (-) supercoils, 2.81 nM for (+) supercoil relaxation). Topoisomerase concentrations and conditions were selected to deliberately yield relatively low topoisomerase activities. At the indicated times, reactions were quenched as with topo IV assays. Topoisomers were analyzed by electrophoresis. For reactions with G6PD or GcrA, assays proceeded with identical conditions as with GapR, using 400 nM (monomer) to maintain equimolar protein. For reactions with E. coli gyrase, purified GyrA and His-TEV-GyrB were reconstituted and diluted as with Caulobacter gyrase. Gyrase was used at 2.3 nM for (+) supercoil relaxation and 6.8 nM for (-) supercoil introduction. QUANTIFICATION AND STATISTICAL ANALYSIS Calculation of growth rates For wild-type cells the GapR depletion strain in xylose (Figure S1C), doubling time was calculated from a linear fit of the log-transformed OD600 values. Quantification of protein levels by immunoblotting Image quantification and analysis was done with Fiji/ImageJ. GapR band intensities were normalized against an internal standard. Estimated copy number of GapR per cell was calculated from a linear regression of a standard curve generated from purified GapR. The mean and SEM of at least two replicates for each condition are reported. Analysis of RNA sequencing Data analysis was performed using custom scripts written in Python 2.7.6. Single-end sequencing reads were aligned to Caulobacter NC011916.1 with bowtie using the following parameters: bowtie -v 1 -M 1–trim3 30 (trimming the last 30 bp of each read, allowing for 1 mismatch, and if a read has more than 1 reportable alignments, report 1 of them at random). The read coverage was mapped to the Caulobacter genome by assigning each mapped base a value of 1/N where N equals the length of the read. To calculate mRNA abundance, the number of reads mapped to a gene was divided by the length of the gene to yield the mean number of reads per kilobase per million sequencing reads (RPKM). The fold change in gene expression for the GapR depletion was computed by taking the RPKM ratio for each gene from GapR-depleted cells to wild-type cells grown in glucose. We compared the gene expression of the GapR-depletion strain (GSE100657, Table S2) and cells treated with novobiocin (Table S2) to the following datasets: UV and MMC microarrays, which were merged and averaged to form a consensus ‘‘DNA damage’’ dataset (Modell et al., 2011); a microarray of the induction of the toxin SocB (ArrayExpress E-MEXP-3990); RNA-seq after treatment with 4% ethanol (GSE90030); and a microarray of a SciP depletion strain (GSE22062). We computed the Pearson correlation for GapR depletion and genes that were changed at least 2-fold in at least 2 of 3 conditions (UV, MMC, and SocB). LexA regulated genes were defined based on the presence of a consensus LexA-binding site (Modell et al., 2011). Ethanol-responsive genes were defined as those that were differentially expressed at least 4-fold in ethanol treated cells (432 genes). Supercoilingresponsive genes were defined as those differentially expressed at least 4-fold in novobiocin treated cells (289 genes total). We computed the Pearson correlation for the supercoiling regulon genes between the each dataset and novobiocin treatement data. ChIP analysis For GapR ChIP, reads were mapped to the Caulobacter crescentus NA1000 reference genome (NC011916.1) using BWA. Read extension and pile-up was performed using the MACS software package 1.4.2 (with d = 200 and keeping all tags at the same location). The data were further process using custom scripts written in Python 2.7.3. The data were smoothed by convolution with a Gaussian (m = 0, s = 50bp, x = (200 bp, +200 bp)). The data were then normalized in reads per million (rpm). ChIP of RNA polymerase subunit rpoC (rpoC::rpoC-3xFLAG) in PYE was analyzed similarly.
Cell 175, 583–597.e1–e16, October 4, 2018 e14
To identify DNA sequence motifs enriched in GapR bound sequences, we isolated regions with GapR ChIP signal above 1.5 rpm (regions less than 150bps apart were merged). A 200 bp window centered at the maximum signal intensity of each of these regions was retrieved and submitted to DREME for sequence motif analysis. In Figures 2C, 2D, and 2F, the AT content at each base pair was computed using a centered 100 bp sliding window. In Figure 2E, GapR ChIP enrichment for each gene was measured by taking the sum of the GapR signal over the coding region and normalizing by the length of the gene. In Figure S2A, promoters were defined as the 300 bp region upstream of transcription start sites identified from RNAP ChIP-seq. For Figure S2A, GapR ChIP signal at promoters was computed by taking the sum of the GapR signal over the defined promoter region divided by 300 (the length of the promoter region). To analyze the GapR ChIP signal at the ends of transcription units, we first identified transcription units in the Caulobacter genome based on the microbes online operon annotation and we included single genes when not annotated in the microbes online operon database if they were located in between two genes found on the opposite DNA strand. We then looked for a change in enrichment over the rifampicin treated signal associated with the ends of transcription units. To do this, we computed the average enrichment (untreated-rifampicin treated) over 1000 bp following and preceding the end of a transcription unit (last base of the last gene in the transcription unit), we then computed the difference of these two enrichments as a measure of GapR signal (see Figure S2D). As a control, we applied the same metric to the starts of transcription units (first base of the first gene of the transcription unit). In a separate analysis, we computed the average ChIP signal enrichment (untreated-rifampicin treated) at the ends and starts of high or low expression transcription units for each base pair over a 2000 bp window (Figure S2E). To do so, we took the difference of these two ChIP signals (untreated-rifampicin treated) at each base pair over a 2000 bp window centered at the ends or starts of transcription units and computed the average over all high or low expression transcription units. In both analysis, we limited our analysis to transcription units at least 1000 bp long and we used RNA-seq data of the GapR-3xFLAG strain to categorize transcription units as high (> 150 RPKM) or low (< 50 RPKM) expression. When a transcription unit contained more than one gene, the expression value of the last gene in the operon was used. Analysis of replisome progression by flow cytometry Flow cytometry was analyzed with FlowJo from 50,000 total events for each experiment. For each experiment, 1N was calculated from the mean of the distribution of G1 cells, and 2N was assumed to be twice the value of 1N. To estimate the fraction of cells that initiated, the fraction of cells with DNA content greater than the 95th percentile of 1N cells was calculated. To estimate replication rate, a linear regression of the mean staining of the 70-71st percentile of cells versus time was determined while accumulation of DNA content was in the linear range. The slope of this regression line was used as the replication rate. At least two independent experiments were used for each analysis, which are reported as mean ± SEM. For experiments examining replication progression after treatment with novobiocin, average DNA content was calculated using the mode of the distribution, as replisomes in GapR-depleted cells displayed heterogeneous behavior. The replication rate was calculated as a linear fit of the mode of the DNA content versus time. At least two independent experiments were used for each analysis, which are reported as mean ± SEM. Analysis of DNA replication from DNA sequencing Data analysis was performed using custom scripts written in Python 2.7.6. Paired-end sequencing reads were aligned to Caulobacter NC011916.1 with bowtie2 using the following parameters: bowtie2 –N 1 (mismatches in seed set to 1). The genome was shifted 1 MB to the right (3016947 to end, 0 to 3016947) to ensure coverage of the origin region. Duplicate reads were then removed using samtools rmdup –S. The central nucleotide of each pair of reads was then mapped to the Caulobacter genome. The Caulobacter genome was divided into 40,001 bins of 1000 bp each, and the total number of reads in each bin was determined for every sample. Samples were then normalized by dividing by the binned data from the corresponding 0 min (1N chromosome content) control, to reduce sequencing biases. The normalized data was then divided by the median of the 20 kb surrounding ter region to set the copy number of ter as 1. The normalized data was then smoothed by the nonparametric LOWESS method in Python, with a smoothing variable (frac) of 0.04 and plotted against the genomic position. LOWESS smoothed data was analyzed to measure initiation and replication rate. To measure replication initiation at each timepoint, the maximum read enrichment within the 200 kb region surrounding the origin was determined from smoothed data and divided by the mean value of a 20 kb unreplicated region near the ter. The initiated value decreases in wild-type over time as replication forks reach ter since the number of reads at ter is set to 1. Replication rate was determined as follows. The fork position was first defined as the position with DNA content halfway between the maximal and minimal DNA content found at 20 min. For each 20 min increment (i.e., 20-40 min, 40-60 min, 60-80 min), the replication rate was then calculated by estimating the distance traveled by the replication fork. This analysis was then repeated six times for positions with DNA content evenly-spaced between the maximal and minimal DNA content. These seven measurements were averaged and reported as mean ± SEM. Quantification of cell division Image analysis was performed with the MicrobeJ suite in Fiji/ImageJ. Cells were recognized with built-in algorithms with threshold offset set to 20. Non-cells were excluded based on size thresholds. Cell segmentation was performed with the built-in suite with tolerance set to 0.11. Over 250 cells for each time point was counted to determine the fraction of pre-divisional cells.
e15 Cell 175, 583–597.e1–e16, October 4, 2018
Quantification of origin-terminus ratio by qPCR On each plate, relative quantities of cDNA in a given sample were calculated by comparison to a least-square fit on a 2-fold dilution standard curve (Cp versus log-transformed standard fold dilution). Terminus region quantities were normalized to origin region quantities for each time point, and normalized by the origin abundance after synchrony (0 min). At least two independent experiments were used for this analysis. Quantification of DNA binding and topoisomerase assays Image quantification was done with Fiji/ImageJ. To calculate binding affinities and hill coefficients, the amount of bound and free oligonucleotide was determined as a fraction of the total intensity present in the lane. Curves were fitted using the Prism 7 software, using the Specific binding with Hill slope algorithm. At least two independent experiments were used for each analysis, and are reported as mean ± SEM. For topoisomerase assays, band intensities were measured and normalized based on the intensity of the input. For relaxation of (+) and (-) supercoils, the intensity of the most supercoiled band was measured and normalized by the total intensity. For relaxed plasmids, the fraction of the total intensity with < 3 topoisomers was measured and normalized by the total intensity. For Figure S7F and other negative supercoiling assays with gyrase, the fraction of the total intensity with > 4 apparent topoisomers was measured and normalized by the total intensity. At least two independent experiments were used for each analysis, which are displayed as mean ± SEM at each time point. The fold-stimulation due to GapR was calculated from the difference in the slope of the change in supercoiled fractions over time. DATA AND SOFTWARE AVAILABILITY The accession number for the RNA-seq, ChIP-seq, and DNA-seq data reported in this paper is GEO: GSE100657. The accession numbers for the coordinates and structure factor amplitudes for the C. crescentus GapR-DNA, apo Bosea sp. GapR, and Bosea sp. GapR-DNA complexes reported in this paper are PDB: 6CG8, 6CFY, and 6CFX, respectively.
Cell 175, 583–597.e1–e16, October 4, 2018 e16
Supplemental Figures
proteobacteria
δ/ε gentle lysis
α
84/90
γ
0/107
β
0/57
α-proteobacteria
DAPI staining (a.u.)
nucleoid chromosome fragments
0 -1
rre ccu
nce c
er ons
vat
2
3
4
5
6
7
8
Rhizobiales
100%
Caulobacteraceae
80%
Rhodobacterales
60%
Rhodospirillales
40%
Sphingomonadales
20%
Rickettsiales
0%
9 10 11
extract nucleoid fraction mass spectrometry of nucleoid-associated proteins
Ga
pR
-3x
FL
time (min)
WT + xyl
∆gapR Pxyl-gapR +glu
AG
+ novobiocin DNA damage SocB induced
ΔgapR Xyl::PgapR-gapR
GapR
[GapR] (pmoles)
migration into gradient (cm)
1.25
2.5
0.625
0.313
0.156
G
GapR
2259 ± 803 copies/cell
ΔgapR Pxyl-gapR + glu
+ ethanol SciP depletion -2
purified GapR
0.499 ± 0.178 pmoles
ethanol response
∆gapR Pxyl-gapR +xyl
GapR-3xFLAG
+xyl
∆gapR Pxyl-gapR
+ glu
*
WT
τ = 84 min WT
E
ΔgapR Pxyl-gapR + xyl
WT WT
3d
+glu
∆gapR Pxyl-gapR τ = 125 min
fraction #
F
+ xyl
ion
D 1
2d
+glu
OD600
o
sheared DNA + RNA
1
∆gapR Pxyl-gapR
+ xyl
+xylose
11
2
WT
0/41
separate lysate in sucrose gradient fraction 1
C
# GapR homologs
PYE (rich media)
B
M2 (minimal media)
A
log2(expression)
2
5 WT+novo
4
∆gapR Pxyl-gapR
3 2
WT
1 0 0.0
0.2
0.4
0.6
0.8
1.0
[EtBr] (μg/mL)
Figure S1. Identification of GapR, a Highly Conserved a-Proteobacterial Protein, Related to Figure 1 (A) Schematic of nucleoid purification mass spectrometry assay, with a representative DAPI staining profile. For proteins identified in mass spectrometry, see Table S1. (B) GapR is highly conserved in a-proteobacteria. The number of GapR homologs identified in and the total number of genomes searched for each class of proteobacteria (top). The occurrence (reported as a percentage of genomes searched) and conservation (reported as average of BLAST bit scores) of GapR homologs is shown for each genome in the individual clades of a-proteobacteria (bottom). (C) Growth of wild-type and DgapR Pxyl-gapR cells assessed by 10-fold serial dilutions onto PYE (rich medium) plates containing either glucose or xylose (top). Plates were imaged after the indicated number of days. Growth in liquid rich medium containing xylose (middle). Growth of wild-type and DgapR Pxyl-gapR cells on M2 minimal medium plates containing either glucose or xylose (bottom). (D) Growth of wild-type, DgapR Pxyl-gapR, and a strain where DgapR is complemented by gapR expressed from its native promoter integrated at the xyl locus (DgapR xyl::PgapR-gapR). (E) GapR-depleted cells exhibit a general stress response. Expression changes of DgapR Pxyl-gapR cells in xylose or glucose, cells treated with novobiocin, DNA damaging agents, the toxin SocB, ethanol stress, or cell-cycle arrest by SciP depletion. Genes that change > 2-fold in ethanol are shown. (F) Immunoblots of GapR in wild-type and the gapR::gapR-3xFLAG strain (top). Asterisk (*) indicates a non-specific band used as a loading control. GapR-3xFLAG levels are 2.2 ± 0.47 fold higher than wild-type GapR levels, mean ± SEM, n = 2. Immunoblots quantifying GapR protein levels in wild-type cells, using purified GapR standards (bottom). Mean ± SEM, n = 4. (G) Ethidium bromide (EtBr) nucleoid sedimentation analysis to determine relative levels of (-) supercoiling in wild-type, novobiocin-treated, and GapR-depleted cells grown for 2.5 hr in glucose. LOWESS fit of data is shown with solid lines. Approximate concentration of EtBr required to relax the chromosomes is shown with dashed lines.
RNAP
0
GapR rif
RNA-seq (rpm)
-1
GapR
-2
0
2
4
+ GapR
4
3
2
1
6
log2(GapR-depleted fold change)
RNA-seq (rpm)
-2
RNAP ChIP (rpm)
1
GapR rif-treated ChIP at AT rich peaks (rpm)
5
GapR ChIP (rpm)
2
-3 -4
C % AT content
B
GapR ChIP at promoters log (rpm) 2
A
0
- GapR
0
1
2
3
4
5
GapR ChIP at AT rich peaks (rpm)
1442000
1444000
1446000
1448000
1450000
1452000
1454000
chromosome position (bp)
2.
Calculate rif-treated GapR ChIP signal 1000 bp windows following and preceding end of TU
end of TU
GapR
3.
4.
Subtract reverse window from forward window to isolate end of gene pile-up and normalize by window size
low expression TUs:
GapR rif
end of TU
GapR GapR rif
end of TU Subtract rif-treated 1000bp windows
E
GapR GapR rif
1000
5' end 3' end
high expression TUs:
5' end 3' end
0.06
0.04
(no rif) - (+ rif)
1.
Calculate GapR ChIP signal 1000 bp windows following and preceding end of TU
change in GapR occupancy
D
0.02
0
-0.02
-0.04
-0.06
-1000
-500
0
500
1000
distance from start or end of gene (bp)
= >0
Figure S2. GapR ChIP-Seq, Related to Figure 2 (A) GapR binding and gene expression changes are not correlated. Scatterplot of gene expression changes in GapR-depleted cells versus average GapR ChIPseq enrichment at promoters. Blue shaded area encompasses genes that change < 2-fold upon GapR depletion. (B) GapR binding and gene expression changes at a ribosomal protein locus. ChIP signal (middle) and AT content (top) are plotted as in Figure 2D, RNA-seq data show reads per million (rpm), and the positions of annotated genes are indicated (bottom, with ribosomal operon genes in solid black). (C) GapR ChIP signal is unchanged at AT rich peaks in untreated versus rifampicin-treated cells. (D) Flow chart of the ChIP analysis performed to identify transcription-dependent GapR occupancy at the end of transcription units. (E) Meta-analysis of GapR occupancy at the 50 and 30 ends of genes after rifampicin treatment. Normalized change in GapR ChIP occupancy is calculated from the difference in enrichment (untreated minus rifampicin treated cells) over a 1 kb window before and after the 30 or 50 ends of transcription units (TUs). High expression TUs have RPKM > 150, low expression TUs have RPKM < 50. Only TUs > 1000 bp were analyzed.
A
B isolate G1 cells (1N)
120 min 100 min release
grow in xylose
80 min
xylose or glucose
60 min 1N
40 min
2N
G1
S
G2/M
20 min
sample every 20 min
1N 2N
1N 2N
1N 2N
+ xyl
+ xyl
+ xyl -> glu
stain DNA, measure DNA content with flow cytometry or DNA sequencing
WT
C
E
fraction initiated
∆gapR Pxyl-gapR WT + xyl
1.0
+xyl
0 min
∆gapR Pxyl-gapR WT ∆gapR Pxyl-gapR + xyl -> glu
+xyl -> glu
∆gapR Pxyl-gapR + glu (2 hr)
0.8
2.0
ori
ter
time post release ori
1.8
0.6
1.6
0.4
1.4
0.2
1.2
20 min
1.0
0.0 0
20
40
60
80
100
ter
0.8 2.0
120
time post release (min)
ori
1.8
D relative DNA content
replication rate (% of WT)
1.6
100 80 60 40 20
40 min
1.4 1.2 1.0
ter
0.8 2.0
ori
1.8 1.6
60 min
1.4 1.2
0
1.0
WT + xyl
+ xyl
+ xyl -> glu
∆gapR Pxyl-gapR
ter
0.8 2.0
ori
1.8 1.6
80 min
1.4 1.2 1.0
ter
0.8
3043
0
1000
2000
3000
chromosome position (kb)
Figure S3. Replication Initiation and Elongation Is Impaired in GapR-Reduced Cells, Related to Figure 3 (A) Schematic of assay to assess DNA replication. Wild-type and DgapR Pxyl-gapR cells were grown in xylose, synchronized in G1, and released into medium with xylose (+ xyl) or glucose (+ glu). DNA content was measured with SYTOX staining and flow cytometry or DNA sequencing at various times post release. (B) DNA replication time courses for wild-type and GapR-reduced strains. Blue lines represent 1N and 2N DNA content. (C) Fraction of cells from (B) that initiated replication over time, with DNA content > 1N used as a proxy for initiation. Data are mean ± SEM, n R 2. (D) Replication rate for cells in (B), calculated from a linear fit of the DNA content of cells over time and normalized to wild-type. Data are mean ± SEM, n R 2. (E) DNA sequencing of wild-type cells grown in glucose, DgapR Pxyl-gapR cells grown in xylose and released into glucose after synchrony as in (B), and DgapR Pxyl-gapR cells depleted for 2 hr in glucose before synchrony (from Figure 4E). DNA content, normalized to that of synchronized G1 cells, is shown as a function of chromosome position using LOWESS smoothing with 10 kb windows. Blue lines represent locations of the replication origin and terminus. Schematics indicate progression of wild-type replication forks.
A ∆gapR Pxyl-gapR
WT time post release
time post release 240 min
240 min 210 min
210 min
180 min
180 min
150 min
150 min
120 min
120 min
90 min
90 min
60 min
60 min
30 min 1N
30 min
0 min
2N
1N
B % predivisional cells
n > 250
0 min
2N
WT 90 min
WT 120 min
WT 210 min
GapR depletion 90 min
GapR depletion 120 min
GapR depletion 210 min
WT
30
20
10
∆gapR Pxyl-gapR
0 0
60
120
180
240
time post release (min)
C
WT
D
∆gapR Pxyl-gapR
drug
+glu
--
2.0
ori/ter ratio
+xyl
WT
1.5
+HU +rif
1.0
--
0.5 0 0
60
120
180
time post release (min)
240
∆gapR Pxyl-gapR
+HU +rif
Figure S4. Replication Elongation Defects of GapR-Depleted Cells, Related to Figures 3 and 4 (A) Extended DNA replication time course for wild-type and GapR-depleted cells. Cells were grown in xylose, washed into glucose for 2 hr, synchronized in G1, and released into medium with glucose. DNA content was measured with SYTOX staining and flow cytometry for 4 hr after release. (B) GapR-depleted cells do not successfully divide. Cells treated as in (A) were fixed and imaged to assess cell division by quantifying the percent of predivisional cells (n > 250 cells per time point, left). Representative images of wild-type and GapR-depleted cells at 90, 120, and 210 min post release into glucose (right). (C) Origin-terminus ratio of wild-type and GapR-depleted cells. DNA was isolated from cells treated as in (A) and the relative copy number of the origin and terminus was determined by qPCR. One representative wild-type dataset is shown, and the mean ± SEM is reported for GapR-depleted cells, n R 2. (D) Sensitivity of GapR-depleted cells to 0.3 mg/mL hydroxyurea (HU) and 0.3 mg/mL rifampicin (rif).
A
B
10
15
20
20bp
50bp
elution volume (mL)
25
C 0.7
20 0.5 15 0.3
10
0.1
5
hill coefficient 1.0
0.9
fraction bound
molecular mass (kD)
polydispersity Mw = 1.000 Mn
absorbance at 280nm
GapR 26.4 kD
30
10
11
0.5
0.0 16
0 9
12
64
256
1024
E 20.0
5.0
10.0
2.5
1.0
0.5
0.25
0
D
4096
16bp
1.46 ± 0.27
20bp
1.89 ± 0.11
50bp
3.45 ± 0.29
150bp
1.97 ± 0.12
210bp
1.50 ± 0.10
300bp
1.71 ± 0.14
pUC19
3.16 ± 0.63
16384
[GapR] (nM monomer)
elution volume (mL)
[GapR] (μM)
0 0.01 0.05 0.1 0.25 0.5 1.0 2.5 5.0 10.0 20.0
[GapR] (μM)
14.3 kD
A260nm 5
[GapR] (μM)
+
F + +
-
(-)
(+)
(+)
-
+
-
L -
(+ )
L
M
(-)
20.0 10.0 5.0 2.5 1.0 0.5 0.25 0
absorbance
A280nm
0 0.01 0.05 0.1 0.25 0.5 1.0 2.5 5.0 10.0 20.0
GapR
[GapR] (μM)
+ + + + + + + + DNase
+ 0.4 μM GapR + 30 °C 1 hr
pUC19
Figure S5. In Vitro Characterization of GapR Binding to DNA, Related to Figure 5 (A) Gel filtration trace of purified GapR (70 mM monomer) (top). Absorbance at 280 nm (blue) and at 260 nm (red) are shown. Arrow indicates approximate elution peak for lysozyme monomer (14.3 kD). SEC-MALS analysis of purified GapR (1 mg/mL = 100 mM monomer) (bottom). Apparent molar mass calculated from MALS (red) and absorbance at 280 nm (gray) are shown. GapR is monodisperse (Mw/Mn = 1.000), with an observed mass of 26.4 kDa (expected 20.7 kDa for homo-dimer). (B) EMSAs of GapR binding to linear dsDNA. Purified GapR was incubated with 20 (left) and 50 bp DNA (right). (C) GapR binding curves for various substrates (left). Apparent KD for each substrate is: 16 bp, 4605 nM GapR; 20 bp, 745 nM; 50 bp, 583 nM; 150 bp, 406 nM; 210 bp, 314 nM; 300 bp, 220 nM; and pUC19 589 nM. Data are mean ± SEM, n R 2. Hill coefficients ± SE (right). (D) EMSA of GapR binding to (-) supercoiled pUC19. (E) Trace nuclease contamination in GapR preparation. Relaxed, supercoiled, and linear standards compared to relaxed or (+) supercoiled pUC19 plasmid incubated with 0.4 mM GapR for 1 hr at 30 C, followed by reaction quenching, protein removal with Proteinase K, and electrophoresis. (F) DNase protection by GapR. GapR was incubated with pUC19 and treated with DNase. Reactions were processed as in (E).
A
B
4
bits
3 2 1
M
*
D
A
SE
D AAA
E DT T
I
T
P
D
D S D E
N T
S
A
I
A
E
N
N
S
D
G
S S
N
I
P
M
S
D
G V
G
E
A
D A E E G K V
T P
E
L Q S N
M
D
T
G N D
S V
N F
0
G E P
H T
V
A
T
D A T E
V L
S
M
Q
V R
V
L
C. cres. Caulo. K31 B. subvib. φ Phage Cr30 A. pasteur. R. palus. A. tume. S. mel. Wolbachia φ Phage Cd Sphi. sp. SKA58 Bosea sp. R. mona.
8H
S
S
T
A
H
G
I
S
G
E N A
K
**
L ER ERLE EK K V EAK G D K R RK E E L Y ALG
G RI
R
T
N
V
P V
K
N S Q T S A D
Q
L
E G N
*
AT R FVI I E K LI ADDI ED YA G FY VT VLM KVI VI RL D EDR E A I EDL L M
V ADEQ
G
AYR I
A D
T G A F T DS I
E V
*
V
L
S
H
Q
V
K
S
S Q
KA I
L L Y
TTT
V
A D
Y
KF L
18G
M A D D A I - - P H T D V L N S T A Q GQ L K S I I E R V M A D D A A - - P H A D V L N S T A Q GQ L K S I I D R V M A D D A S F D A SP D V L N SA AQ GR L R T I I ER I M T E V A S - - - - V D V I N T T A K SQ L K S I I D R V M D G ER I SGD ED A A V GG I A A D R L R S I I ER V M A T S A A - A V Q E D P A T N F A K DQ L R A I I E R I - - - - - - - M D D T S T T E T V A A A E L RQ I I E R I - - - - - - - - - - M S D A H G I A R DQ L R A F I E R I - - - - - - - - - - M ED T V K I T A E E L K GY I ER I - - - - - - - - - - M SD V G E I A A D R L K S F V ER I - - - - - - - - - - - M S E G N V A A DQ L R L F I E R I - - - - - - - - - - - - M D D P V A G DQ L K S I V E R I - - - - - - - - - - - - M S E V V V K EQ L EQ Y I S K I
A
Q
K
D
Q
I
M
I
A D AT SEQQ
S
28E
G
R
DR I
V
H
G M A
V
E
Q
V
R
K
A
G
G
N V N N
N
S
Q
Q
S
T
38M
L
M
M
AR
N
S
I FG T
L
N
T
L S
S
P
E A
V
A
M S
M
L
48A
S
A Q A SI I IK K T
TV
E
LL
S A
58L
Q
R KD Q I E
V TK QD L M M
PA A T N K R AH E M KL M V RMK D
M
K
E
E
Q
Q
Q
I
Q
S
L
M S
K
E
H
S
Q
A
T
68D
A
Q S
M Q I V A
V G
M
78I
I
E R
R
I
H
I
EA I
I
L E S
L M E V
T
88E
E R L E V E K A E I M EQ I K E V Y A E A K G N G F D V K V L K K V V R I R K Q D R A K RQ E E D A I L D L Y L S A I G E I - - - - - - - - E R L E V E K T E I A EQ I K E V Y L E A K G N G F D V K I L R K V V R L R K T D R A K RQ E E D A I L D L Y L S A I G E I - - - - - - - - E R L E E D K A A V M E DM K E V F A E A K G E G Y D V K I L R K V I R I R K Q D K A K RQ E E D A I L D L Y L S A L G E I - - - - - - - - E R L V I E K K E I Q EQ V K E V L N E A K G N G F D T K I I R K V I R I RQ Q D R A K RQ E E E S I L D L Y L T A L G E I - - - - - - - - E R L E E E R K A L A G D I K D I F T E A K S A G F D V K V I RQ I I R L R K Q E P A E V E EQ E T L L D I Y R R A L GM - - - - - - - - - E R L E E E K K T I S D D I R D V Y A E A K G N G F D V K A L R T I V RM R K Q D A N E R A EQ E T I L E T Y MQ A L GM L - - - - - - - - E R L E E E K A A I Q G D I K D V I G E A K G R G Y D T K A I R T I I R L R K K D A N E R I E E E T I L Q T Y M A A L GM E - - - - - - - - E R L E E E K K T I A D D I K D V Y G E A K SM G F D A K I L R K V I S I R K Q D A D E RM EQ E A I L D T Y L Q A L GM I P A V E A E D A A E K L EQ E K R D V Q G H I R D V Y A K A A D E GWD I K V M K Q I I R L R K M D D D D R E EQ E I L L D T Y K R A L GM S Y E E E L S E - E R L E E E K S G L Q E D I K E V Y S E A K G T G F D V K I I RQ I I R L R K M D K A D RQ EQ R A I L EM Y E E A L GM T E - - - - - - - E R L E E E K K G I G D D I K D V Y L E A K A N G Y D G K I M RQ I V R L R K MQ P H D RQ EM E A I L Q T Y L A A L GM E - - - - - - - - ER L E E EK K T I A D D I K EV Y A EA K GN GY D V K V L R K V I A I R K R D A N ER A E E EA I L D L Y LQ A V G E SA - - - - - - - E R L EQ E K A D L SQ E V K D I FQ D A S S H G F D V K A M K S I L K L K K L D K D K L A EQ D A M L E L Y R D T L G I - - - - - - - - - -
Conservation - - - - - - - - - - 00231157447 *7579699 *9 *44994496448999756 *452 *9 *6 *79959949974722643 *6759 *77 *4389 *50- - - - - - - - -
Consensus
V L N S+ A A DQ L + S I I E R I E R L E E E K K + I + E D I K D V Y A E A K G N G F D V K I L R K I I R L R K Q D A A + RQ E+ E A I L D L Y L + A L GM
α1 α1
C
α2
α2
α3
dimer dimer-of-dimers
D
Figure S6. GapR DNA Structure and Contacts, Related to Figure 6 (A) GapR-DNA contacts in stereoview. Shown are the side chains from the GapR dimer-of-dimers that interact with the DNA. (B) Sequence logo of GapR computed from select homologs (top). Conserved acidic surface residues that are not involved in structural contacts or DNA binding (*). Sequence alignment of selected GapR homologs, organized by distance from C. crescentus GapR, with residues numbered based on C. crescentus GapR. Helices of GapR are shown for both the dimer and dimer-of-dimers. (C) Superimposition of the DNA bound forms of the C. crescentus GapR (magenta) and the Bosea GapR (green), showing that they form the same dimer-of-dimers structure. (D) Fluorescence polarization assay of C. crescentus GapR binding to AT-rich DNA (blue circles), GC-rich DNA (red squares), and a DNA template without AT stretches (green triangles).
+topo IV
+topo IV
(+)
(+)
0
10
20
30
10
+topo IV
+topo IV
(+)
20
30
10
+topo IV
(+)
20
30
10
(+)
20
30
10
B
+topo IV (+)
20
30
10
20
+topo IV
(-) (-)
0
7.5
15
% initial supercoiled DNA
A
60 time (min)
30
30
* no GapR
C
0.25 μM
0.5 μM
1.0 μM
(+) (+)
0
7.5
15
30
5.0 μM GapR
+GapR
*
D
+topo IV (-)
(-)
2.5 μM
60 time (min)
0
+ E. coli topo IV 15
30
45
60
0
7.5
15
30
60
80 60
-GapR
40 20 0
+GapR 0
10
20 30 40 time (min)
50
60
E
+ E. coli topo IV
(-) (-)
time (min)
100
(+)
time (min)
(+)
+
+
+
+
+
+
+ gyrase
5 2.5 1.0 0.5 0.25 0.1 0 [GapR] μM
*
+ 0.4 μM GcrA
F
+gyrase
(-)
0
(-)
10 20 30 40 50 60 time (min)
*
+GapR
*
*
100 80
60 40 20 0
0
10
G
20 30 40 time (min)
50
60
80
10
20
30
40
50
60 min
60 40 20 0
0
10
20 30 40 time (min)
H
(-)
0
25 min
100
+gyrase
(-)
10 min
+GapR
% initial supercoiled DNA
% initial supercoiled DNA
+ 0.4 μM G6PD
*
(+)
60
+ E. coli gyrase
0
60 time (min)
50
15
30
45
+ E. coli gyrase 60
time (min)
0
15
30
45
(-)
60
time (min)
* *
*
+GapR
*
+GapR
-GapR
80 60 40 20 0
+ 0.4 μM GcrA
+GapR 0
10
20 30 40 time (min)
50
60
I
recruit
GapR
80
60 40 20 0 0
10
20 30 40 time (min)
relaxes (+)sc
+ relaxes (+)sc
* 100
% initial relaxed DNA
100
% initial supercoiled DNA
% initial relaxed DNA
+GapR + 0.4 μM G6PD
60
80
60 40 20 0 0
10
20 30 40 time (min)
50
60
+ topoisomerase
+ topoisomerase
50
100
+ GapR
GapR
Figure S7. Gyrase and Topo IV Are Specifically Stimulated by GapR, Related to Figure 7 (A) GapR stimulates topo IV. Topo IV (2.34 nM) was incubated with the indicated amounts of GapR before adding (+) supercoiled pUC19 and ATP. After incubating for the times indicated, reactions were quenched and analyzed by gel electrophoresis. Slower migration in the gel indicates relaxation of (+) supercoils. (B) GapR modestly stimulates relaxation of (-) supercoils by topo IV. Reactions are the same as in (A) except with 17.55 nM topo IV, 0.4 mM GapR, and (-) supercoiled pUC19. Relaxed and (-) supercoiled standards are shown (left). Plot quantifies the % initial supercoiled plasmid (asterisk) at each time point (right). Data indicate mean ± SEM, n R 2. (C) Relaxation of (-) supercoils by topo IV is not stimulated by non-specific proteins. Assay is the same as in (B) except with 0.4 mM E. coli glucose-6-phosphate dehydrogenase (G6PD) or C. crescentus GcrA. A representative gel from at least two experiments is shown. (D) GapR does not strongly stimulate E. coli Topo IV. Assay is the same as in (B) except with E. coli topo IV. (E) GapR stimulates gyrase. Assay is the same as in (A) except with 6.25 nM gyrase. Relaxed and (+) supercoiled controls (left). Gyrase generates slower migration as (+) supercoils are relaxed and then induces faster migration as (-) supercoils are introduced. (F) GapR stimulates introduction of (-) supercoils by gyrase. Assay is the same as in (E) except with 12.5 nM gyrase, 0.4 mM GapR, and relaxed pUC19. Faster migration in the gel indicates introduction of (-) supercoils. Plot quantifies % initial relaxed plasmid (asterisk) at each time point (bottom). Data indicate mean ± SEM, n R 2. (G) Introduction of (-) supercoils by gyrase is not stimulated by non-specific proteins. Assay is the same as in (F) except with 0.4 mM G6PD or GcrA. (H) GapR stimulates E. coli gyrase. Assay is the same as in (F) except with E. coli gyrase (2.3 nM, (+) relaxation; 6.8 nM, (-) introduction). Plot quantifies % initial supercoiled plasmid (asterisk) at each time point (bottom). Data indicate mean ± SEM, n R 2. (I) Models for stimulation of topoisomerase activity by GapR. (Left) GapR may bind overtwisted DNA in front of replication forks, directly recruiting or activating type II topoisomerases to relax positively supercoiled DNA. (Right) GapR could bind overtwisted DNA flanking plectonemic DNA, thereby trapping the substrate of type II topoisomerases.