Jul 01, 2023
High
Nature (2023)Cite this article 449 Accesses 304 Altmetric Metrics details Antibiotic binding sites are located in important domains of essential enzymes and have been extensively studied in the
Nature (2023)Cite this article
449 Accesses
304 Altmetric
Metrics details
Antibiotic binding sites are located in important domains of essential enzymes and have been extensively studied in the context of resistance mutations; however, their study is limited by positive selection. Using multiplex genome engineering1 to overcome this constraint, we generate and characterize a collection of 760 single-residue mutants encompassing the entire rifampicin binding site of Escherichia coli RNA polymerase (RNAP). By genetically mapping drug–enzyme interactions, we identify an alpha helix where mutations considerably enhance or disrupt rifampicin binding. We find mutations in this region that prolong antibiotic binding, converting rifampicin from a bacteriostatic to bactericidal drug by inducing lethal DNA breaks. The latter are replication dependent, indicating that rifampicin kills by causing detrimental transcription–replication conflicts at promoters. We also identify additional binding site mutations that greatly increase the speed of RNAP.Fast RNAP depletes the cell of nucleotides, alters cell sensitivity to different antibiotics and provides a cold growth advantage. Finally, by mapping natural rpoB sequence diversity, we discover that functional rifampicin binding site mutations that alter RNAP properties or confer drug resistance occur frequently in nature.
In addition to its importance in combating clinical resistance2,3,4, the study of antibiotic resistant mutants has also advanced our understanding of fundamental biological processes5. Because antibiotics function by targeting enzymes necessary for bacterial survival, resistance mutations have been instrumental in the characterization of essential cellular machinery such as RNA polymerase6,7,8,9,10 (RNAP) and the ribosome11,12,13,14. Resistance mutations have been shown to inadvertently alter the function of antibiotic targets, providing new insights into the various functional6,11,15, regulatory16 and physiological attributes12,17,18 of essential enzymes, identifying antibiotic binding sites as phenotypic hotspots.
The binding site of rifampicin (Rif), a specific inhibitor of bacterial RNAP, is well defined by high-resolution structures8,10,19 that confirm earlier genetic findings from resistant mutants20. High-affinity binding of Rif prevents the extension of nascent RNA beyond a length of three nucleotides and Rif binding is lost in resistant mutants. Collectively, resistance mutants manifest a fascinating breadth of phenotypes that are a consequence of changes in the properties of the mutant enzyme at each step of the transcription cycle15,16,21,22,23,24. Although Rif is notorious for high-frequency resistance4,25,26, the known spectrum of resistant mutants makes less than 5% of the non-redundant mutational space at this region of the β-subunit of RNAP. Therefore, our functional understanding of the Rif binding site is still incomplete.
Multiplex automated genome engineering (MAGE) in tandem with sequencing is free of the constraints of positive selection and low mutation rates27,28,29. This system enables the introduction and monitoring of mutations to a region of interest on the bacterial chromosome with extremely high efficiency and specificity1,30,31 and has previously been used to characterize codon optimization and sequence diversity across the essential genes infA and rpoD27,28.
In this work, we generate a collection of RNAP mutants that span all possible substitutions in each position of the Rif binding site. We provide a detailed genetic map of Rif binding interactions and identify an alpha helix spanning rpoB 521–526 where mutations can either increase or decrease Rif potency. In this region, we discover single-residue mutants that convert Rif from a bacteriostatic to a bactericidal drug by increasing its binding affinity and causing DNA damage. We show that the bactericidal effects are dependent on active replication, indicating that RNAP trapped by Rif at promoters causes detrimental transcription–replication conflicts. We delineate other single-residue mutations in the Rif binding site that produce fast RNAP and determine that high transcription speed depletes the cell of nucleotides, increasing susceptibility to antibiotics, but can also be beneficial as it provides a growth advantage at lower temperatures through poor termination efficiency. Despite high evolutionary conservation, we identify functional Rif binding site mutations that occur at a high frequency in nature. Our work demonstrates that the Rif binding site is a multifunctional domain of RNAP, and its characterization broadens our understanding of the antibacterial potential of Rif and the relationship between transcription speed and bacterial physiology.
We generated and monitored a saturated pool of single amino acid substitutions in rpoB positions 510–537 and 563–572 using MAGE in tandem with sequencing27 (Fig. 1a,b). These positions have been shown to interact with Rif either structurally or genetically8,10 and encompass resistance mutants most commonly found in clinical isolates9. Following one cycle of MAGE, we prepared sequencing libraries to estimate the efficiency of oligo recombineering and retrieved most of the expected substitutions at a frequency between 10−3 and 10−5 (Fig. 1c and Extended Data Fig. 1a). We were able to reproducibly detect greater than 99% of generated mutants before selection, which is critical in our experimental design for the identification of sensitive mutants (Fig. 1d).
a, The Rif binding site within the rpoB subunit of E. coli RNAP (in blue), magnified in the lower panel (PDB 5UAC). b, Schematic diagram of the mutant pool generation. MAGE oligos were designed to generate a saturated pool of 760 single-codon mutants (20 substitutions × 38 positions) within the Rif binding site. c, Representative heatmap of mutant abundances generated with MAGE at timepoint 0 (T0) and timepoint 1 (T1), after fitness screening. For analysis of screening data, normalization was performed between T1 and T0 by calculating the log2(T1 fitness/T0) fold change. Rows correspond with rpoB positions 510–537 and 563–572; columns correspond to amino acids. White squares denote wild-type residues. d, Schematic of drug sensitivity and resistance screen. Mutants that fail to grow in LB are designated weak mutants whereas mutants that only fail to grow in drug are considered drug sensitive. Drug-specific effects can be isolated by normalizing ‘T1 drug’ mutant abundances to ‘T1 fitness’, that is, log2(T1 drug/T1 fitness).
Source data
To distinguish drug sensitive from inviable mutants, we began by monitoring mutant fitness in regular growth conditions. We prepared samples from cells in liquid culture after two subsequent 100-fold dilutions, allowing cells to reach the same optical density (OD) before collection (Extended Data Fig. 1a). We also prepared samples after passaging the mutant pool on agar plates (Extended Data Figs. 1a and 2c,d). As expected, stop codons are quickly depleted (Fig. 1c and Extended Data Fig. 1a–d). We also observed that the distribution of fitness effects in this region of rpoB is bimodal, with a peak representing neutral substitutions at 0 and the other representing detrimental substitutions at −2 (Extended Data Fig. 2b). Other systematic mutagenesis studies in essential proteins have observed the same bimodal distribution27,28. When using 0 as the threshold for inviable mutants, 58% of substitutions are inviable, which is comparable to values of 53% and 48% for rpoD and infA, respectively27,28.
Next, we screened our collection of rpoB mutants with Rif, bicyclomycin (BCM), and 5-fluorouracil (5FU) and normalized the data to isolate drug-specific effects (Fig. 1d and Extended Data Figs. 1–3). We selected BCM, an inhibitor of Rho termination, and 5FU, a nucleotide analogue, because Rif resistant mutants have previously been shown to be deficient in Rho termination32 and highly sensitive to nucleotide analogues4. Interestingly, this relatively small region of RNAP is rich in mutations that provide differential sensitivity towards these compounds (Fig. 2a and Extended Data Fig. 1). Rif resistance occurs in known positions; however, the range of substitutions conferring resistance is much wider (Figs. 2a and 3a). We note that some mutations, such as those at position 531, are relatively featureless despite being known to cause Rif resistance10. This is because these mutants are generally slow growing, as demonstrated by fitness screening, and therefore difficult to detect as they are highly resistant in the screening format (Extended Data Figs. 1a and 2a; Methods). Resistance to BCM is observed in mutations in the region around position 567 (Fig. 2a,b and Extended Data Fig. 3a,b,e,f). 5FU resistance also appears to occur in this region (Fig. 2a,b and Extended Data Fig. 3c,d,g,h).
a, Aggregate heatmap of rpoB mutant fitness and sensitivity to Rif, BCM and 5FU. Columns correspond to each screening modality and rows represent each amino acid position. b, Protein structure of the Rif binding site with a colour map of average positional fitness or sensitivity to Rif, BCM or 5FU. Resistance or high fitness mutations are represented in red. Z-scores are calculated for each screening modality from MAGE-seq experiments conducted in n = 3 biologically independent replicates.
Source data
a, Sensitivity and resistance landscape of rpoB mutants towards Rif. Rows correspond with rpoB position and columns denote all 20 amino acid substitutions, grouped by property. White squares denote wild-type residues. Bar plot of average positional fold change is provided to the right of the heatmap. b, Volcano plot of data in a. Stop codons are in grey. c, MIC of Rif in hypersensitive mutants (P = 1.2 × 10−4, 2.1 × 10−3 from left to right). d, Meta-analysis of ChIP–seq signal from wild-type (WT) and rpoB T525D RNAP around transcription start sites (TSS) of highly transcribed genes after 1 h recovery from 1 h treatment with 80 µg ml−1 Rif. e, Plating efficiency of hypersensitive mutants following 1 h treatment of 80 µg ml−1 Rif (n = 6; P = 5.7 × 10−4, 1.8 × 10−4). f, DNA damage measured by TUNEL staining in hypersensitive mutants following 1 h treatment with 80 µg ml−1 Rif. TUNEL positive cells are the percentage of cells exceeding signal detected in more than 99% of untreated cells (n = 4; P = 1.3 × 10−2, 4.7 × 10−3). g, Plating efficiency of hypersensitive mutants in recA or recBCD knockout strains after 1 h treatment of 80 µg ml−1 Rif (n = 4, P = 5.3 × 10−2, 7.3 × 10−3, 4.2 × 10−5, 5.7 × 10−5). h, Number of DSBs per 10 kilobases (kb) at promoters and over the gene body measured using SLR-qPCR following 1 h treatment with 80 µg ml−1 Rif (P = 2.4 × 10−3, 1.7 × 10−5, 1.2 × 10−4, 3.7 × 10−3, 2.4 × 10−5, 4.2 × 10−5). i, Plating efficiency of hypersensitive mutants in the PC2 background containing a temperature-sensitive allele of the replication protein dnaC (P = 1.8 × 10−3, 7.6 × 10−5). Replication was blocked by shifting cells to 42 °C for 1.5 h before 30 min treatment with 80 µg ml−1 Rif. Bar plots represent the mean of n = 3 biologically independent experiments unless noted otherwise, with error bars denoting 95% confidence intervals. P values are calculated using the independent t-test (two-sided).
Source data
By calculating the average fitness of all substitutions in a given position (Fig. 2b), we provide a measure for mutability of the wild-type residue. As expected for a highly conserved region of an essential protein, most positions are highly immutable with substitutions to R529, G566 and G570 producing the largest average fitness defects (Fig. 2b and Extended Data Fig. 2). On the other hand, the alpha helix at positions 523–527 appears to be highly mutable (Fig. 2b and Extended Data Fig. 2). Increased average positional fitness indicates that the original wild-type residue is disadvantageous in our growth conditions in some way.
In a similar manner, average positional resistance towards each drug measures the wild-type residue’s importance in drug activity. Multiple mutations in positions 516, 526 and 572 increase average resistance towards Rif (Figs. 2b and 3a), indicating that the wild-type residue in those positions is important for Rif binding. High-level resistance in some positions is caused by a single particular substitution, such as R529K or P564K (Fig. 3a), implying that the substituted residue interferes with Rif binding. Mutations at position 567 increase average resistance towards both BCM and 5FU (Fig. 2b and Extended Data Fig. 3), indicating that the original proline increases sensitivity towards these drugs. Conversely, average positional sensitivity measures the wild-type residue’s interference with drug activity. For example, the wild-type threonine at 525 appears to interfere with Rif binding as many mutations at this position increase sensitivity (Figs. 2b and 3a). Positional sensitivity or resistance provides an indication of the functional and structural importance of each position in single amino acid resolution.
When examining the behaviour of individual mutations in response to Rif, we were intrigued to find Rif hypersensitive (RifHS) mutants in positions 514, 521 and 525 (Fig. 3a), with substitutions L521Y (21Y) and T525D (25D; throughout, mutants from our collection are abbreviated to the final two numbers of the position and the substituted amino acid) producing the strongest phenotype (Fig. 3b). These substitutions have no fitness defect, indicating that their dilution when grown with Rif is not an artefact of mutant weakness (Fig. 3b and Extended Data Fig. 4a). They are also located on an alpha helix in close proximity to mutations that confer high-level resistance (Figs. 2b and 3a).
RifHS mutants 21Y and 25D were reconstructed in a wild-type MG1655 background for further characterization. In agreement with results from the primary screen, wild-type cells resumed growth 2 h after high exposure to Rif whereas RifHS mutants took more than 8 h to recover (Extended Data Fig. 4b). These mutations also result in a decreased minimal inhibitory concentration (MIC) for Rif (Fig. 3c and Extended Data Fig. 4c).
To directly test Rif’s interaction with RifHS RNAP, we cloned mutant 21Y for recombinant expression and purification. In a modified runoff experiment in which polymerase is saturated with Rif, washed and then allowed to resume transcription, 21Y is slower to recover, indicating tighter Rif binding (Extended Data Fig. 4d). This is not due to differences in catalytic activity as the rate of 3-mer synthesis is unchanged, 3-mer being the maximum product length of a polymerase obstructed by Rif8. To test whether Rif causes RifHS RNAP to linger at promoters in vivo, we performed chromatin immunoprecipitation sequencing (ChIP–seq) on RNAP after treating cells with Rif for 1 h, washing and allowing cells to recover for an additional 1 h. After recovery, more RifHS RNAP can be found persisting at promoters when compared to wild-type cells (Fig. 3d and Extended Data Fig. 4e).
Rif is known to be bacteriostatic in Escherichia coli. To test RifHS mutants for cell killing, cells were treated with Rif, washed and serial dilutions were assayed for growth on either agar plates or liquid broth (LB)18. Notably, although Rif remains bacteriostatic in wild-type cells, it becomes bactericidal for both RifHS mutants in E. coli, decreasing viability by at least one order of magnitude (Fig. 3e and Extended Data Fig. 5a–d). This effect is specific to Rif as chloramphenicol and tetracycline remain bacteriostatic (Extended Data Fig. 5e).
Previous studies report increased sensitivity of DNA-damage repair mutants in response to Rif treatment33. To investigate the bactericidal effects of Rif, we used terminal deoxynucleotidyl transferase deoxyuridine-triphosphate (dUTP) nick-end labelling (TUNEL)34 to assess DNA fragmentation after 1 h of Rif treatment. Wild-type cells exhibit weak TUNEL staining of around 15% of positive cells after Rif treatment, which is lower than any bactericidal antibiotic tested in previous work35 (Fig. 3f). Both RifHS mutants show increases in TUNEL staining with Rif treatment of up to 35%, similar to other bactericidal drugs such as 5FU and norfloxacin (approximately 30%)35.
To determine whether DNA breaks contribute to the bactericidal effects of Rif, we generated knockouts of major double strand break (DSB) repair enzymes recA or recBCD in RifHS mutants. Notably, Rif becomes potently bactericidal in both RifHS mutants in the absence of DSB repair (Fig. 3g). In addition, wild-type cells lacking recA become sensitive to Rif killing (Fig. 3g), providing additional evidence that Rif causes low levels of DNA damage in wild-type cells (Fig. 3f). We note that these effects are Rif specific, as RifHS mutants deficient in DSB repair remain insensitive to other bacteriostatic drugs (Extended Data Fig. 5f).
As a secondary indicator of DNA damage, we measured the expression of key genes in the SOS response. On Rif treatment, we see a significant increase in the expression of genes involved in DSB repair (Extended Data Fig. 6a). We also measured protein levels of the SOS gene repressor LexA following Rif treatment and observed that LexA protein levels significantly decrease following Rif treatment, indicating that it is cleaved in response to DNA damage (Extended Data Fig. 6b,c). In agreement with previous work, the bactericidal effect of Rif does not entail increased cellular respiration and the accumulation of toxic reactive oxygen species (ROS)36,37 (Extended Data Fig. 6d–f).
To investigate the location of DNA damage in RifHS mutants, we used semi-long-range quantitative PCR (SLR-qPCR) to detect DNA breaks at specific genomic loci (Extended Data Fig. 7a)38. We arbitrarily selected two highly expressed σ70-dependent genes at distant chromosomal loci and compared the frequency of DSBs at cognate promoters to gene bodies (Extended Data Fig. 7b). Notably, we observed a significant increase in DSBs at promoters after Rif treatment in both wild type and RifHS mutants, whereas no increase occurred over the respective gene body (Fig. 3h). Furthermore, RifHS mutants exhibited a greater level of promoter-specific DSBs than wild type (Fig. 3h).
Previous work has shown that the transcription initiation complex can collide with or interrupt the replisome, resulting in interference with replication fork progression or deleterious mutations39,40. To determine whether replication conflicts at promoters are involved in DNA damage by Rif, we introduced dCas9 targeting the oriC to RifHS strains to partially block replication before Rif treatment41 (Extended Data Fig. 7c). We found that blocking replication significantly decreases the frequency of DSBs at promoters (Extended Data Fig. 7d,e). Next, we constructed RifHS mutants in cells containing a temperature-sensitive allele of the replication protein dnaC2 (PC2)39,42. We verified that shifting cells to the non-permissive temperature 42 °C inhibits replication (Extended Data Fig. 7f). We then compared the bactericidal effects of Rif with and without replication. Notably, we observe that Rif is no longer bactericidal in RifHS mutants when replication is inhibited (Fig. 3i). RifHS mutants in the PC2 background remain hypersensitive at the permissive temperature of 30 °C (Extended Data Fig. 7g).
In summary, we demonstrate that: (1) RifHS RNAP with Rif are more persistently bound at promoters; (2) the bactericidal effect of Rif in RifHS mutants is associated with double-stranded DNA breaks at promoters; (3) Rif induces the SOS response; and (4) Rif killing is much more severe in cells lacking DSB repair. Interestingly, Rif killing is not associated with a buildup of ROS, indicating that it acts in a manner distinct from other bactericidal antibiotics. Instead, we observe that bactericidal effects and DNA damage from Rif are replication dependent, implying that these effects are a result of transcription–replication conflicts at promoters.
The potential of mutations in the Rif binding pocket to affect Rif binding is obvious, but the mechanism by which they cause differences in sensitivity towards other drugs is not readily apparent. By pairwise comparison of the phenotypic screens, we noticed correlations between fitness, BCM and 5FU sensitivity whereas no such correlations exist with Rif (Fig. 4a). High fitness mutants are extremely sensitive to 5FU, especially substitutions disrupting the α-helix at positions 523–527 (Fig. 2a). In addition, rpoB mutant sensitivity towards BCM and 5FU is positively correlated, hinting at a common underlying mechanism (Fig. 4a). We selected mutants with strong phenotypes from each screening modality for reconstruction and further characterization in MG1655 (Extended Data Figs. 1–3). We found the MIC of all reconstructed rpoB mutants to closely reflect screening results. Mutants 15G, 15Y, 23P, 23W, 26V and 27W are highly sensitive to 5FU, with around a fivefold decrease in MIC (Fig. 4b); these mutants are also sensitive to BCM (Fig. 4c). Mutants 12K, 67C, 67V and 72L all display high resistance towards both BCM and 5FU (Fig. 4b,c).
a, Scatterplots of rpoB mutant response to 5FU compared to fitness or mutant response to Rif and BCM. Linear regression lines, 95% confidence intervals and Pearson R values are overlayed. Two-sided P values for 5FU correlations with fitness, Rif and BCM datasets are 2.4 × 10−43, 0.11 and 5.2 × 10−32, respectively. b,c, MIC of 5FU (b) and BCM (c) for sensitive or resistant strains reconstructed in MG1655 (n = 2). d, Transcription elongation speed of wild type and mutants. Transcription of five regions on the lacZ gene is measured after induction and plotted as a function of time and distance. The slope of the best-fit line represents transcription elongation speed. e, Mean pause frequency in wild type and 5FU sensitive mutants as measured by NETseq. Values denote the aggregated mean pause frequency per kilobase for the top 3,500 expressed genes. n = 2 biologically independent experiments for each condition (P = 1.6 × 10−3, 9.8 × 10−25, 7.2 × 10−22 from left to right). f, Pyrimidines in 5FU sensitive mutants as detected by LC–MS, normalized to labelled UTP (n = 5; P = 2.4 × 10−6, 5.6 × 10−4, 8.9 × 10−5, 6.4 × 10−5, 5.0 × 10−2, 8.4 × 10−3). g, Plating efficiency of 5FU sensitive mutants with and without thymidine supplementation. Cells from overnight cultures were serially diluted and plated on LB agar plates with the indicated concentrations of 5FU in µg ml−1 with and without 1 mg ml−1 thymidine (P = 3.1 × 10−3, 8.6 × 10−3, 4.7 × 10−3). h, Plating efficiency of 5FU sensitive mutants with the indicated concentrations of trimethoprim (P = 2.5 × 10−4, 7.6 × 10−4, 3.3 × 10−3). Bar plots represent the mean of n = 3 biologically independent experiments unless noted otherwise, with error bars denoting 95% confidence intervals. P values are calculated using the independent t-test (two-sided).
Source data
BCM acts by inhibiting the Rho transcription termination factor. Previous work has shown that polymerase speed determines the efficiency of Rho termination32. To test whether mutants sensitive to both 5FU and BCM have faster RNAP, we measured their elongation rate using a multiprobe qPCR assay along the lac operon and their pause frequency using native elongating transcript sequencing (NETseq)4,43. Indeed, these mutants have a higher elongation rate (Fig. 4d) and pause less frequently (Fig. 4e) than wild-type RNAP.
We considered the possibility that high transcription rates in fast RNAP mutants deplete nucleotide pools, resulting in 5FU sensitivity. A known target of 5FU is thyA, encoding the enzyme responsible of dTMP synthesis. dUMP, the substrate for this enzymatic reaction, is produced from pyrimidine ribonucleotides. Therefore, a strain on the pool of pyrimidine ribonucleotides would sensitize cells to 5FU44,45. Using hybrid liquid chromatography–mass spectrometry (LC–MS), we found that 5FU sensitive mutants are depleted of the two detectable pyrimidines, uridine and cytidine triphosphate, to the level of an upp mutant deficient in uracil salvage (Fig. 4f). 15Y, a mutant incredibly sensitive to 5FU, is depleted of purines as well (Extended Data Fig. 9d). Indeed, thymidine supplementation rescued 5FU sensitivity of mutants 15Y, 26V and 27W (Fig. 4g), as well as the growth defect of 15Y under normal growth conditions at 37 °C (Extended Data Fig. 9e). Of note, modest changes in the expression level of nucleotide biosynthesis genes do not account for 5FU sensitivity of fast mutants (Extended Data Figs. 8 and 9a–c).
Because 5FU inhibits bacterial growth through a complex mechanism that includes not only dTMP depletion, but also misincorporation to DNA and RNA, we sought to establish sensitivity to dTMP depletion using thymidine-specific assays. First, we tested whether fast mutants were more susceptible to ‘thymineless death’ by constructing fast mutants in a genetic background of thyA, which requires thymidine supplementation for growth. When thymidine is removed, cells are able to subsist on the existing reservoirs of dTMP during a ‘resistance phase’ before viability is lost46. In support of faster depletion of intracellular thymidine, fast mutants in a thyA background have a shorter resistance phase than wild-type RNAP (Extended Data Fig. 9f). Notably, the well-studied rpoB 513P, a slow RNAP15,47, is more resistant to 5FU and has a protracted resistance phase in the thymineless death assay (Extended Data Fig. 9g,h). Second, we tested sensitivity of fast mutants towards trimethoprim. This clinically important antibiotic inhibits dihydrofolate reductase, which is required for de novo synthesis of thymidine in addition to certain amino acids48. Because LB medium is limited for thymidine, its availability affects overall trimethoprim sensitivity. Consistent with our model for nucleotide depletion, fast RNAP mutants are also more sensitive to trimethoprim (Fig. 4h). Together, these results indicate that RNAP speed determines the availability of nucleotides within the cell.
All fast RNAP mutants we isolated appear to have a fitness advantage in the original screening format (Extended Data Figs. 1a and 2), which was done at lower temperatures suitable for MAGE recombineering strains1. When reconstructed in MG1655, fast mutants have no such growth advantage at 37 °C, but all grow faster and saturate at a higher OD than wild-type cells at a lower temperature of 25 °C (Extended Data Fig. 10a). 5FUR mutants grow slowly at 25 °C (Extended Data Fig. 10b). At 42 °C, fast mutants 15Y and 27W are temperature sensitive and grow more slowly (Extended Data Fig. 10b). Previous work has shown temperature sensitivity of naturally occurring Rif resistant mutant 526Y, which is also known to be a fast RNAP7. Because fast polymerases terminate inefficiently, we hypothesized that increased readthrough might provide a growth advantage at 25 °C. Indeed, wild-type cells with subinhibitory amounts of BCM grow faster at 25 °C (Extended Data Fig. 10c). From these data, we identify that fast RNAP provides a growth advantage at lower temperatures, at least in part through poor termination efficiency.
In summary, our characterization of fast RNAP mutants of the Rif binding site led to the observation that the speed of transcription controls the availability of nucleotides within the cell. Nucleotide pools are limiting and marginal increases in RNAP speed greatly sensitize cells to antibiotics that target nucleotide biosynthesis. In addition, fast RNAP also affects desirable growth environments for bacteria, raising the question of whether such mutations appear in nature.
We aimed to inform our study of rpoB mutants by exploring natural sequence diversity. Following alignment of around 4,000 orthologues of E. coli rpoB mined from KEGG49, we observed high evolutionary conservation at the Rif binding site by Jensen–Shannon divergence50 (Fig. 5a). We next grouped species by phylum or class and determined the number and type of substitutions that separate orthologues from E. coli MG1655 rpoB within the Rif binding site. Within this 38 amino acid window, most species available for analysis contained fewer than five substitutions, with mild substitutions appearing in closely related proteobacteria (Fig. 5b,c). As one would expect in a functionally conserved region of RNAP, the most common substitutions are between amino acids of the same functional group with the exception of S522A, a hydroxylic to aliphatic substitution found in close to 10% of species surveyed (Fig. 5b,c and Extended Data Fig. 2).
a, Evolutionary conservation of rpoB by Jensen–Shannon divergence based on alignment of 4,055 sequences. Higher values indicate higher conservation. Positions are based on MG1655 rpoB and the Rif binding site is highlighted in grey. b,c, Total number of substitutions (b) and frequency of specific residue substitutions (c) that separate the surveyed 4,055 rpoB orthologues from MG1655 grouped by phylum. Proteobacteria are further subdivided by class. Phyla with fewer than 40 representative species are combined in ‘other’. Substitutions found in fewer than 20 species are omitted for brevity. In the box plot, the lower and upper ends of the box denote Q1 and Q3, respectively. The whiskers span 1.5 × (Q3 – Q1) from each side of the box. d, Scatterplots of substitution frequency in nature compared to substitution fitness, sensitivity and resistance. Substitutions found in less than 20 species are omitted for brevity. Orange denotes substitutions with significant (P < 01) phenotypes in each screen. P values are calculated using the independent t-test (two-sided).
Source data
When comparing average positional fitness to evolutionary conservation, we see a negative correlation as expected (Extended Data Fig. 10d). Generally, highly conserved positions are less mutable, except in some instances. For example, histidine at 526 is highly conserved but is mutable in our fitness screen. Mutating this position is probably deleterious in other ways, which is why it is rarely found in evolution. For example, mutations in this position are highly sensitive to 5FU (Extended Data Fig. 10d).
Average positional fitness is limited in determining the effects of single-residue variants found in nature. To speculate how natural variance in the Rif binding site might evolve, we compared the frequency of natural variance to substitution fitness in each screening modality (Fig. 5d). Among the most common natural variants with strong phenotypes in our screen, 24L, 26N and 67Q are the most notable. 24L is found in more than 10% of species surveyed, has characteristics of a fast RNAP and appears infrequently in Proteobacteria, but more commonly in Actinobacteria and Firmicutes (Fig. 5c,d and Extended Data Figs. 2 and 3). 26N is found in 1% of species surveyed, has characteristics of a fast RNAP with high-level Rif resistance and appears in Tenericutes (Fig. 5c,d and Extended Data Figs. 2 and 3); notable genera include Mycoplasma, which have been shown to be Rif resistant9. 67Q provides resistance to both BCM and 5FU and is found in -Epsilonproteobacteria and the closely related Aquificae; notable species include Helicobacter and Campylobacter (Fig. 5c,d). Some strains of Epsilonproteobacteria have been shown to be resistant to 5FU51. Future work will need to be done to characterize the effects of the identified mutations in these organisms.
To build upon decades of work with RNAP resistant mutants while avoiding the bias of characterizing antibiotic binding sites through the lens of resistance, we present a complete collection of single-residue mutants of the Rif binding site. With our collection, we observe that the Rif binding site is surprisingly mutable and that functional mutations are not exclusive to Rif resistance and cannot be fully appreciated without saturating mutagenesis (Figs. 1 and 2 and Extended Data Figs. 1–3). We expand upon the known Rif ‘resistome’ by identifying that nearly every substitution at positions 516, 526 and 572 generates high-level resistance (Figs. 2 and 3a). Many of these resistant mutations do not occur naturally, probably because they require multiple nucleotide substitutions. Moreover, we identify a single alpha helix that has a critical role in Rif binding, where mutations either cause high-level resistance (S522 and H526) or increased sensitivity (L521 and T525) (Figs. 2 and 3a). These data provide a detailed genetic mapping of ligand–enzyme interactions that will be of value to medicinal chemists working to increase antibiotic potency.
We discover single-residue mutations of RNAP that convert Rif into a bactericidal drug in E. coli (Fig. 3). We show that these mutations act by prolonging Rif binding as the mutant polymerase is slow to recover from Rif treatment both in vivo and in vitro (Fig. 3d and Extended Data Fig. 4d), the behaviour of other bacteriostatic drugs is unchanged (Extended Data Fig. 5e,f), the mutants have no fitness defect (Extended Data Fig. 2) and we observe no gene expression reprogramming (Extended Data Figs. 8 and 9a–c). In a similar way, next-generation macrolides are bactericidal by virtue of higher affinity to the ribosome than that of the parental bacteriostatic drug14,52.
DNA-damage repair mutants have previously been observed to have an increased sensitivity to Rif33. In this work, we directly observe that Rif causes double-stranded DNA breaks at promoters, induces the SOS response and kills RifHS mutants (Fig. 3 and Extended Data Figs. 5a–d and 6a–c). We show that bactericidal effects are exacerbated in RifHS mutants deficient in DSB repair and also in wild-type recA cells, indicating that Rif damages DNA in wild-type cells as well, but not to the extent of causing cell death (Fig. 3g). In agreement with previous work, we find that the bactericidal effects of Rif are not dependent on ROS (Extended Data Fig. 6d–f), indicating a mechanism distinct from other bactericidal antibiotics36,37. Instead, we observe that Rif causes DNA damage and kills only when cells are actively replicating (Fig. 3i and Extended Data Fig. 7). RNAP arrested at the promoter is known to present an obstacle for the replisome and transcription–replication conflicts have been shown to cause detrimental mutations39,40. Our results indicate that Rif damages DNA through a similar mechanism. Future work will need to be done to test this model in other bacteria. These findings imply that improving Rif binding affinity can shorten treatment times and extend its antimicrobial spectrum of activity. Moreover, the discovery that Rif acts through a unique mechanism to cause DNA damage may inform the design of combination antimicrobial therapy.
We discover in our collection other single-residue mutations that yield fast RNAP, similar to the well-characterized Rif resistant H526Y mutant and other mutants of RNAP15,21,53. We observe that high transcription rates deplete pyrimidine nucleotides, such that these mutants become particularly vulnerable to further depletion through 5FU inhibition of thyA54 (Fig. 4). Importantly, growth media supplemented with thymidine are able to rescue 5FU sensitivity of fast RNAP (Fig. 4g). We note that thymidine supplementation does not render cells fully resistant to 5FU. Misincorporation of dUTP and its 5-fluoro analogue during replication, as well as 5FU misincorporation during transcription, may contribute to the antimicrobial effects of 5FU. Nucleotide depletion by fast RNAP also makes cells more susceptible to thymineless death (Extended Data Fig. 9f) and to the clinically important antibiotic trimethoprim (Fig. 4h). High-speed transcription also sensitizes cells to the antibiotic BCM because of poor termination efficiency; however, we discover that poor termination provides a growth advantage at lower temperatures (Extended Data Fig. 10c). Together, these results demonstrate that transcription speed determines nucleotide availability, can alter antibiotic sensitivity and can change the optimal growth temperature of bacteria.
Despite high evolutionary conservation of the Rif binding site, we identify functional RNAP mutations that exist at a high frequency in nature. 524L and 526N both share characteristics with fast RNAP mutants and are located on the alpha helix we identified to be important for Rif binding (Figs. 2 and 3). As such, 526N provides high-level Rif resistance and increases sensitivity to 5FU, but not to the same degree as 526Y (Extended Data Fig. 3). It is intriguing to see that clinical resistance, which develops relatively quickly, provides a larger fitness cost than resistance that developed over millions of years of evolution. It is also interesting to speculate why Tenericutes and other phyla would carry this mutation. Perhaps they lived in close vicinity to Rif producing bacteria. rpoB 567Q, found in 2% of species surveyed, provides high-level resistance to both 5FU and BCM. We notice that species carrying these mutations live in harsh environments, such as hydrothermal vents and animal digestive tracts. These mutations were not a focus of this Article but they are probably slow, stringent mutants of RNAP, similar to 513P16 (Extended Data Fig. 9g,h), and we speculate that they increase tolerance to living in nucleotide- or nutrient-limiting conditions.
Together, our high-resolution characterization of the Rif binding site reveals several RNAP subtypes, some of which can be found in nature. We find RNAP that is exquisitely sensitive to Rif, converting it to a bactericidal antibiotic in a gram-negative bacterium through a unique mechanism of DNA damage. We also find that transcription speed determines the availability of nucleotides within the cell and can alter antibiotic susceptibility. Our study illuminates the antibiotic binding site as a multifunctional domain of essential enzymes and inspires similar characterization of other major antibiotic targets to broaden our understanding of antibacterial mechanisms, bacterial physiology and evolution.
MAGE with the E. coli EcNR2 strain1 was used to generate a pool of single-residue Rif binding site substitutions comprising rpoB 510–537 and 563–572. 760 MAGE oligos (Integrated DNA Technologies) were designed to generate 20 possible single-residue mutants in 38 positions. To limit the detection of background mutations and further increase the specificity of our approach, we specifically designed our MAGE library to produce single codon substitutions that include multiple nucleotide alterations when possible. A single round of MAGE mutagenesis was done in four bins comprising 510–518, 519–527, 528–537 and 563–572. Cells were recovered overnight and stocked in glycerol.
To begin a screening experiment, cells were thawed, bins were combined and allowed to recover to an OD at 600 nm (OD600) of 0.4 at which point an aliquot was taken as T0. Cells were then subject to a ‘recovery format’ for Rif treatment or a ‘MIC format’ for BCM or 5FU (Sigma) treatment. In the ‘recovery format’, cells were treated with 50 µg ml−1 Rif (Sigma) for 1 h, washed with LB, diluted tenfold and allowed to recover to the OD before dilution before being collected (T1 drug). In the ‘MIC format’ cells were diluted tenfold and treated with BCM (15 or 30 µg ml−1) or 5FU (1 or 5 µg ml−1) and allowed to recover to the OD before dilution before being collected (T1 drug). Untreated cells were grown in parallel to measure and control for mutant fitness (T1 fitness). Untreated cells were also passaged on agar plates by plating cells and scraping colonies after overnight incubation. All screening was done in triplicate. The screening format we used depended on the selected drug. Because Rif is a bacteriostatic drug37, we were not expecting sensitive mutants to decrease in abundance. Fortunately, previous work has shown that Rif produces a ‘post-antibiotic effect’ in which cells show a recovery delay after Rif is washed away and cells are allowed to resume growth55,56. For this reason, we designed a recovery format in which mutant rpoB pools were first treated with Rif, washed, diluted tenfold and then allowed to recover. Sensitive mutants would be expected to recover more slowly and would then be found at a lower abundance. As both BCM and 5FU are bactericidal37,54,57, we expected sensitive mutants to die and decrease in abundance. Therefore, two concentrations were selected for each drug: one concentration at the MIC at which only resistant mutants would survive and one at a lower concentration at which sensitive mutants would be lost. We isolated drug-specific effects by normalizing mutant abundances between drug treated (T1 drug) and untreated (T1 LB).
Genomic DNA was prepared from collected cells with the Monarch Genomic DNA purification kit (NEB). Libraries were prepared in two bins comprising rpoB 510–537 and 563–572. Here 100 ng DNA was first subject to two cycles of PCR in a 50 µl reaction containing Q5 polymerase (NEB) and 1 µM primers (Supplementary Table 2) designed to attach a unique barcode of 14 degenerate nucleotides to correct for PCR duplicates during analysis. Following barcode attachment, reactions were treated with 1 µl of exonuclease I (NEB) to remove unused barcodes for 1 h at 37 °C and then DNA was purified using a 1.5× ratio of PCRClean DX beads (Aline Biosciences). Sample barcode indexes and sequencing adaptor sequences were then attached with 18 cycles of exponential PCR before 2 rounds of purification with a 1× ratio of beads. Libraries were sequenced with the Novaseq 6000 SP300 cartridge (Illumina), aiming to collect 40 million reads per library and 2 × 106 unique barcodes for each experimental sample, allowing observation of 20–20,000 instances of each mutant.
Sequencing adaptors were trimmed using CutAdapt58. Next, sequences were filtered for correct length and primer sequences, and then grouped into ‘families’ by the aforementioned nucleotide barcodes. Families representing an original genomic molecule were checked for 80% consensus by nucleotide, defaulting to the wild-type sequence in the case of no consensus and then translated to the designed amino acid sequence. Families were grouped by mutant and the number of families corresponding to a single mutant was divided by the total number of families to yield mutant frequency. Only single-residue mutants were considered for analysis. To determine mutant fitness, mutant abundances after growth (T1 fitness) were normalized to abundances before growth (T0). To determine drug-specific effects, mutant abundances after drug treatment (T1 drug) were normalized to abundances after growth in LB (T1 fitness).
For all validation experiments, wild-type E. coli MG1655 was used as the parental strain (Supplementary Table 1). To generate selected rpoB mutants for reconstruction in MG1655, we used MAGE to introduce a kanamycin cassette at the 3′ end of rpoC in the E. coli EcNR2 strain. Selected rpoB mutants were generated in the resulting E. coli EcNR2 rpoc:kn. P1 phage-mediated transduction was used to move mutants to MG1655. For knockout strains, pKD46 was used as previously described59.
Overnight cultures of the indicated strains grown in LB were diluted 10,000-fold and growth in 0.2 ml of LB was monitored at 37 °C in a Bioscreen C machine (Growth Curves USA). MIC values were defined as the minimal antibiotic concentration for more than 90% growth inhibition after 15 h. For MIC and MBC assays adhering to clinical standards60, an initial inoculum of 5 × 105 cells was incubated with increasing concentrations of Rif in Mueller–Hinton broth at 37 °C for 18 h. The MIC was determined where at least 90% of bacterial growth is inhibited based on OD600. The MBC was determined as the lowest concentration of antibacterial agent that reduces the viability of the initial bacterial inoculum by at least 99.9% after 18 h at 37 °C.
Overnight cultures of the indicated strains grown in LB were diluted 1,000-fold and grown to an OD600 of 0.2 in 2 ml of medium before the addition of the indicated concentration of Rif (Sigma), ampicillin (Sigma), chloramphenicol (Sigma), tetracycline (Sigma), 2,2′-dipyridyl (Sigma) or thiourea (Sigma). After 1 h of treatment, cells were washed twice with an equal volume of LB. Serial dilutions of the washed cells were made and either plated on LB agar or placed in a Bioscreen C machine to monitor growth. For thymidine supplementation experiments, LB agar was supplemented with the indicated concentration of thymidine (Sigma) and 5FU (Sigma). For replication block experiments, strains were constructed in the PC2 background carrying a temperature-sensitive allele of the DNA replication protein dnaC39,42. Cells were grown as indicated to an OD600 of 0.2 at 30 °C before being shifted to 42 °C for 1.5 h. Cells were then added to media with 80 µg ml–1 Rif prewarmed to 42 °C for 30 min before being washed and plated on LB agar. Viability of cells was calculated by normalizing the number of treated to untreated colonies.
The rpoB mutant was cloned onto the PVS10 vector for recombinant expression and purification of E. coli core RNAP as previously described61. Than 3 pmol of RNAP core was mixed with equimolar σ70 and template DNA with sequence: tccagatcccgaaaatttatcaaaaagagtattgacttaaagtc taacctataggatacttacagccATCGAGAGGGCCACGGCGAACAGCCAACCCAATCGAACAGGCCTGCTGGTAATCGCAGGCCTTTTTATTTGGATCCCCGGGTA (capital letters denote transcribed sequence) in 20 µl TB50 buffer (20 mM Tris-HCl pH 8; 10 mM MgCl2, 50 mM NaCl, 0.03% Igepal-60) and incubated 5 min at 37 °C. Assembled transcription complexes were then treated with 2 µl of Rif (20 µg ml–1) and incubated for 5 min at 37 °C before being transferred to 10 µl of NeutrAvidin beads (Fisher) and shaken for 5 min at 25 °C. Samples were washed four times with 1 ml of TB50. Adenosine triphosphate (ATP), guanosine triphosphate (GTP), uridine triphosphate (UTP) and 5 µM cytidine triphosphate (CTP) premixed with CTP-αP32 was added to 1 mM and incubated at 37 °C. Then 10 µl aliquots were taken after 3, 5, 10, 20, 30 or 40 min and mixed with SB (1× tris-borate-EDTA (TBE); 20 mM EDTA; 8 M urea; 0.025% xylene cyanol, 0.025% bromophenol blue), heated for 5 min at 100 °C and loaded on a 30% (20 × 20 cm2) (19:1) polyacrylamide gel with 7 M urea and TBE prerun for 6 min. The gel was run for 50 min at 50 W before being transferred to an X-ray film for overnight exposure. Gel intensity was quantified using GelQuant.NET.
Overnight cultures of rpoB wild-type rpoC-10X-His and rpoB T525D rpoC-10X-His were diluted 1,000-fold and grown to an OD600 of 0.2 in 25 ml of LB at 37 °C before the addition of 80 µg ml–1 Rif. After treatment for 1 h, cells were washed once in an equal volume of LB. Cells were allowed to recover for 1 h in LB at 37 °C before crosslinking with 1% formaldehyde at 37 °C for 20 min. A final concentration of 250 mM glycine was added to quench the reaction. Cells were collected by centrifugation at ×5,000g and washed twice with ice cold 1 phosphate buffered saline (PBS) before being resuspended in 1 ml of lysis buffer (10 mM Tris pH 8.0, 100 mM NaCl, 1 mM EDTA, 1 mM EGTA, 0.1% sodium deoxycholate, 0.5% N-lauroylsarcosine) plus 2 mg ml−1 lysozyme with protease inhibitor (Roche) and incubated at 37 °C for 30 min. DNA was sheared using an ultrasonicator Covaris M220 on a 10 s on/10 s off cycle for a total of 50 cycles. Supernatant was incubated with 5 µg ml–1 6X-His-tag antibody (Proteintech) and incubated at 4 °C overnight.
Then 50 µl per sample of Protein A/G beads were pre-equilibrated with 1 ml of ice cold 1× PBS +0.5% bovine serum albumin (BSA) and incubated with samples for 0.5 h at 4 °C. Samples were washed five times in 1 ml of RIPA buffer (50 mM HEPES pH 7.5, 250 mM LiCl, 1 mM EDTA, 1% NP40, 0.7% sodium deoxycholate) and given a final wash in 1× TE (100 mM Tris-Cl pH 8, 10 mM EDTA pH 8). Complexes were uncrosslinked in 250 µl of 1× TE + 4 µl of RNAseA for 1 h at 37 °C and then left overnight at 65 °C with 1% sodium dodecyl sulfate (SDS) and 1 mg ml−1 proteinase K. DNA was purified with ChIP Clean and Concentrate (Zymo). ChIP experiments were done in duplicate.
For sequencing, sample libraries were prepared using the NEBNext ChIP–seq library kit (NEB) according to the manufacturer’s instructions. Libraries were checked on TapeStation 2200 (Agilent) for quality control. Samples were sequenced on NextSeq 2000 (Illumina). Bowtie62 and Deeptools63 were used for alignment and analysis, respectively.
For DNA fragmentation experiments, cells were diluted 1,000-fold from overnight cultures and grown to an OD600 of 0.2 in 2 ml of medium before the addition of 80 µg ml–1 Rif or 125 ng ml−1 norfloxacin. After 1 h of treatment, cells were washed twice with an equal volume of LB before being fixed with 4% formaldehyde for terminal deoxynucleotide TUNEL. Labelling of DNA fragments was performed using the Apo-Direct Kit (BD Bioscience) following manufacturer’s instructions. Samples were collected with a FACSCalibur flow cytometer (BD Biosciences) and at least 50,000 cells were collected for each sample. The percentage of TUNEL positive cells for a given condition is the percentage of cells exceeding the signal detected in more than 99% of untreated cells.
The number of DSBs in the genomic DNA was determined using SLR-qPCR38. SLR-qPCR can be used for detection of any lesion or strand break that creates a barrier to amplification by DNA polymerase64. To quantify DSBs in genomic DNA, indicated strains were inoculated in LB. Overnight cultures of indicated strains were diluted 1,000-fold in LB and allowed to grow to an OD600 of 0.2 at 37 °C. Cells were treated either with Rif (80 µg ml–1) or the carrier solvent dimethyl sulfoxide (DMSO) for 1 h at 37 °C. Cells were harvested by centrifugation and washed twice with PBS. Genomic DNA was isolated using a Monarch Genomic DNA isolation kit (NEB).
The reaction mixture contained 1× SYBR Green mix (Applied Biosystems), 500 nM of each primer and 1 ng of template DNA in a total volume of 20 µl per well. Two sets of primers (Supplementary Table 3), yielding a short and a long amplicon provide data representing the total amount of template (an internal normalization control) and undamaged DNA, respectively. The two primer pairs had comparable efficiency of amplification and yielded a single PCR product as judged by agarose gel analysis. Data analysis was done as described38 using a modified version of the 2−ΔΔCT method. DSB density was measured using genomic DNA from untreated cells as reference. DNA damage was calculated as DSBs per 10 kilobase DNA.
For replication block experiments, a dual plasmid system expressing dCas9 and a guide RNA targeting the oriC was introduced to cells of interest as previously reported41. Cells were grown to an OD600 of 0.2 at 37 °C as described above, and dCas9 was induced with 200 ng ml−1 Tc for 30 min before adding Rif.
Cells were diluted 1,000-fold from overnight cultures and grown to an OD600 of 0.2 in 2 ml of medium before the addition of 80 µg ml–1 Rif or 10 µg ml–1 ampicillin. Cells were incubated with antibiotic for 30 min before the addition of 10 mM carboxy-H2DCFDA (Thermo), after which cells were incubated for another 30 min. Cells were then washed with an equal volume of PBS, normalized for cell number and transferred to an opaque 96-well plate, after which fluorescent intensity was measured with a microplate reader. A blank well loaded with an equal volume of PBS was used to control for background fluorescence.
Cells were diluted 1,000-fold from overnight cultures and grown to an OD600 of 0.2 in 2 ml of medium before the addition of 80 µg ml–1 Rif for 1 h. Cells were then collected, resuspended in lysis buffer (10 mM Tris-HCl pH 7.5, 4% SDS), mixed with 4× lithium dodecyl sulfate (LDS) sample buffer (Invitrogen), and boiled at 95 °C for 10 min before loading on a NuPAGE 4–12% Bis-Tris precast gel (Thermo) for electrophoresis at 200 V for 40 min. Gel was then semidry transferred to a polyvinylidene difluoride (PVDF) membrane (Millipore) and incubated with the specified antibodies at 4 °C overnight in SuperBlock T20 blocking buffer (Thermo) before development the following day. Antibodies used were E. coli RNAP beta Monoclonal Antibody clone 8RB13 mouse mAB at 1:2,000 (663905; Biolegend) and LexA Antibody (E-7) at 1:200 (sc-365999; Santa Cruz). Gel intensity was quantified using GelQuant.NET.
RNA was prepared from cells treated with 80 µg ml–1 Rif for 1 h using Trizol and Direct-zol RNA microprep (Zymo). RT–qPCR was done with 100 ng of RNA with the Luna RT–qPCR kit (NEB) and Quantstudio 7 (Applied Biosystems) using primers targeting rpsL as a control (Supplementary Table 5).
RNA-sequencing libraries were prepared from mid log cells grown in LB with the NEB Ultra II Directional RNA Library Prep Kit for Illumina (NEB, E7760S). Reads were trimmed with Cutadapt and then aligned to the E. coli genome with Bowtie2 (ref. 62). DESeq2 (ref. 65) was used to analyse differential expression.
Mid log cells grown in LB were spun down, washed with PBS and resuspended in 1 m; of dry-ice-cooled extraction buffer (50% methanol, 1% formic acid, 10 µM labelled UTP; 13C9, 15N2) optimized for UTP extraction. Cells were transferred to Beadblaster (Benchmark Scientific) tubes and homogenized with ten cycles (30 s on, 30 s off). Homogenized samples were transferred to a 2 ml glass vial and vortexed with 0.4 ml of HPLC grade water and 0.8 ml of chloroform. Samples were then incubated at 4 °C for 30 min and the supernatant was separated and dried in a SpeedVac (Thermo Scientific) before being resuspended in 40 µl of HPLC grade water. The resuspended sample was then centrifuged to remove debris and transferred into a 250 µl glass insert for LC–MS analysis at the NYU metabolomics core.
Transcription elongation rates were measured with a multiprobe (Supplementary Table 3) qPCR assay along the lac operon as previously described43. E. coli cells were grown to OD600 of 0.4 before induction by 1 mM isopropyl beta-d-thiogalactopyranoside (IPTG). Following induction, aliquots were withdrawn at 10 s intervals into a tube containing stop solution (60% EtOH, 2% phenol, 10 mM EDTA) precooled to −20 °C. RNA was extracted using Trizol (Invitrogen) and the Direct-zol RNA Microprep kit (Zymo) according to manufacturer’s instructions. Then 100 ng of total RNA was converted to complementary DNA using Superscript IV reverse transcriptase (ThermoFisher) and RNaseOut (ThermoFisher). Real-time qPCR amplification was performed with SYBR Green (Invitrogen) and Quantstudio 7 (Applied Biosystems). Data analysis was done as previously described43.
NETseq libraries were prepared and analysed as previously described4, with modifications to the pause calling algorithm. Pauses were called using 100-nucleotide sliding windows with the condition that the pause is the maximum signal within the window and the signal is 2 s.d. above the window mean.
The KEGG orthology database49 was mined for orthologues of E. coli rpoB and analysed in a similar manner to work previously done with rpoD28. In brief, orthologues were filtered for bacterial proteins that were of similar length to E. coli rpoB yielding 1,342 sequences. Next, sequences were aligned with Clustal Omega66 and positional conservation was calculated by Jensen–Shannon divergence50. Orthologue variants within the rpoB binding site were determined by comparison to E. coli MG1655 rpoB and their number and frequency were calculated.
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Sequencing data is available through the NCBI sequence read archive (SRA) under PRJNA992395. Source data are provided with this paper.
Code for analysing MAGE-seq datasets was adapted from previous work28.
Wang, H. H. et al. Programming cells by multiplex genome engineering and accelerated evolution. Nature 460, 894–898 (2009).
Article CAS PubMed PubMed Central ADS Google Scholar
Imamovic, L. & Sommer, M. O. A. Use of collateral sensitivity networks to design drug cycling protocols that avoid resistance development. Sci. Transl. Med. 5, 204ra132 (2013).
Article PubMed Google Scholar
Baym, M., Stone, L. K. & Kishony, R. Multidrug evolutionary strategies to reverse antibiotic resistance. Science 351, aad3292 (2016).
Article PubMed PubMed Central Google Scholar
Rasouly, A. et al. Analysing the fitness cost of antibiotic resistance to identify targets for combination antimicrobials. Nat Microbiol 6, 1410–1423 (2021).
Article CAS PubMed PubMed Central Google Scholar
Freedy, A. M. & Liau, B. B. Discovering new biology with drug-resistance alleles. Nat. Chem. Biol. 17, 1219–1229 (2021).
Article CAS PubMed PubMed Central Google Scholar
Ezekiel, D. H. & Hutchins, J. E. Mutations affecting RNA polymerase associated with rifampicin resistance in Escherichia coli. Nature 220, 276–277 (1968).
Article CAS PubMed ADS Google Scholar
Jin, D. J. & Gross, C. A. Characterization of the pleiotropic phenotypes of rifampin-resistant rpoB mutants of Escherichia coli. J. Bacteriol. 171, 5229–5231 (1989).
Article CAS PubMed PubMed Central Google Scholar
Campbell, E. A. et al. Structural mechanism for rifampicin inhibition of bacterial RNA polymerase. Cell 104, 901–912 (2001).
Article CAS PubMed Google Scholar
Goldstein, B. P. Resistance to rifampicin: a review. J. Antibiot. 67, 625–630 (2014).
Article CAS Google Scholar
Molodtsov, V., Scharf, N. T., Stefan, M. A., Garcia, G. A. & Murakami, K. S. Structural basis for rifamycin resistance of bacterial RNA polymerase by the three most clinically important RpoB mutations found in Mycobacterium tuberculosis. Mol. Microbiol. 103, 1034–1045 (2017).
Article CAS PubMed PubMed Central Google Scholar
Ozaki, M., Mizushima, S. & Nomura, M. Identification and functional characterization of the protein controlled by the streptomycin-resistant locus in E. coli. Nature 222, 333–339 (1969).
Article CAS PubMed ADS Google Scholar
Galas, D. J. & Branscomb, E. W. Ribosome slowed by mutation to streptomycin resistance. Nature 262, 617–619 (1976).
Article CAS PubMed ADS Google Scholar
Zhang, Y. et al. The context of the ribosome binding site in mRNAs defines specificity of action of kasugamycin, an inhibitor of translation initiation. Proc. Natl Acad. Sci. USA 119, e2118553119 (2022).
Article CAS PubMed PubMed Central Google Scholar
Svetlov, M. S., Vázquez-Laslop, N. & Mankin, A. S. Kinetics of drug-ribosome interactions defines the cidality of macrolide antibiotics. Proc. Natl Acad. Sci. USA 114, 13673–13678 (2017).
Article CAS PubMed PubMed Central ADS Google Scholar
McDowell, J. C., Roberts, J. W., Jin, D. J. & Gross, C. Determination of intrinsic transcription termination efficiency by RNA polymerase elongation rate. Science 266, 822–825 (1994).
Article CAS PubMed ADS Google Scholar
Zhou, Y. N. & Jin, D. J. The rpoB mutants destabilizing initiation complexes at stringently controlled promoters behave like ‘stringent’ RNA polymerases in Escherichia coli. Proc. Natl Acad. Sci. USA 95, 2908–2913 (1998).
Article CAS PubMed PubMed Central ADS Google Scholar
Fredriksson, A., Ballesteros, M., Dukan, S. & Nyström, T. Induction of the heat shock regulon in response to increased mistranslation requires oxidative modification of the malformed proteins. Mol. Microbiol. 59, 350–359 (2006).
Article CAS PubMed Google Scholar
Saito, K. et al. Rifamycin action on RNA polymerase in antibiotic-tolerant Mycobacterium tuberculosis results in differentially detectable populations. Proc. Natl Acad. Sci. USA 114, E4832–E4840 (2017).
Article CAS PubMed PubMed Central Google Scholar
Lin, W. et al. Structural basis of Mycobacterium tuberculosis transcription and transcription inhibition. Mol. Cell 66, 169–179.e8 (2017).
Article CAS PubMed PubMed Central Google Scholar
Jin, D. J. & Gross, C. A. Mapping and sequencing of mutations in the Escherichia coli rpoB gene that lead to rifampicin resistance. J. Mol. Biol. 202, 45–58 (1988).
Article CAS PubMed Google Scholar
Svetlov, V., Belogurov, G. A., Shabrova, E., Vassylyev, D. G. & Artsimovitch, I. Allosteric control of the RNA polymerase by the elongation factor RfaH. Nucleic Acids Res. 35, 5694–5705 (2007).
Article CAS PubMed PubMed Central Google Scholar
Jin, D. J. et al. Effects of rifampicin resistant rpoB mutations on antitermination and interaction with nusA in Escherichia coli. J. Mol. Biol. 204, 247–261 (1988).
Article CAS PubMed Google Scholar
Yanofsky, C. & Horn, V. Rifampin resistance mutations that alter the efficiency of transcription termination at the tryptophan operon attenuator. J. Bacteriol. 145, 1334–1341 (1981).
Article CAS PubMed PubMed Central Google Scholar
Shiver, A. L. et al. Chemical-genetic interrogation of RNA polymerase mutants reveals structure-function relationships and physiological tradeoffs. Mol Cell 81, 2201–2215.e9 (2021).
Article CAS PubMed PubMed Central Google Scholar
Jee, J. et al. Rates and mechanisms of bacterial mutagenesis from maximum-depth sequencing. Nature 534, 693–696 (2016).
Article CAS PubMed PubMed Central ADS Google Scholar
Adams, R. A. et al. Rifamycin antibiotics and the mechanisms of their failure. J. Antibiot. 74, 786–798 (2021).
Article CAS Google Scholar
Kelsic, E. D. et al. RNA structural determinants of optimal codons revealed by MAGE-Seq. Cell Systems 3, 563–571.e6 (2016).
Article CAS PubMed PubMed Central Google Scholar
Park, J. & Wang, H. H. Systematic dissection of σ70 sequence diversity and function in bacteria. Cell Rep 36, 109590 (2021).
Article CAS PubMed PubMed Central Google Scholar
Russ, D. et al. Escape mutations circumvent a tradeoff between resistance to a beta-lactam and resistance to a beta-lactamase inhibitor. Nat. Commun. 11, 2029 (2020).
Article CAS PubMed PubMed Central ADS Google Scholar
Ellis, H. M., Yu, D., DiTizio, T. & Court, D. L. High efficiency mutagenesis, repair, and engineering of chromosomal DNA using single-stranded oligonucleotides. Proc. Natl Acad. Sci. USA 98, 6742–6746 (2001).
Article CAS PubMed PubMed Central ADS Google Scholar
Nyerges, Á. et al. Directed evolution of multiple genomic loci allows the prediction of antibiotic resistance. Proc. Natl Acad. Sci. USA 115, E5726–E5735 (2018).
Article CAS PubMed PubMed Central Google Scholar
Jin, D. J., Burgess, R. R., Richardson, J. P. & Gross, C. A. Termination efficiency at rho-dependent terminators depends on kinetic coupling between RNA polymerase and rho. Proc. Natl Acad. Sci. USA 89, 1453–1457 (1992).
Article CAS PubMed PubMed Central ADS Google Scholar
Liu, A. et al. Antibiotic sensitivity profiles determined with an Escherichia coli gene knockout collection: generating an antibiotic bar code. Antimicrob. Agents Chemother. 54, 1393–1403 (2010).
Article CAS PubMed PubMed Central Google Scholar
Gavrieli, Y., Sherman, Y. & Ben-Sasson, S. A. Identification of programmed cell death in situ via specific labeling of nuclear DNA fragmentation. J. Cell Biol. 119, 493–501 (1992).
Article CAS PubMed Google Scholar
Dwyer, D. J., Camacho, D. M., Kohanski, M. A., Callura, J. M. & Collins, J. J. Antibiotic-induced bacterial cell death exhibits physiological and biochemical hallmarks of apoptosis. Mol Cell 46, 561–572 (2012).
Article CAS PubMed PubMed Central Google Scholar
Lobritz, M. A. et al. Antibiotic efficacy is linked to bacterial cellular respiration. Proc. Natl Acad. Sci. USA 112, 8173–8180 (2015).
Article CAS PubMed PubMed Central ADS Google Scholar
Kohanski, M. A., Dwyer, D. J., Hayete, B., Lawrence, C. A. & Collins, J. J. A common mechanism of cellular death induced by bactericidal antibiotics. Cell 130, 797–810 (2007).
Article CAS PubMed Google Scholar
Rothfuss, O., Gasser, T. & Patenge, N. Analysis of differential DNA damage in the mitochondrial genome employing a semi-long run real-time PCR approach. Nucleic Acids Res. 38, e24 (2010).
Article PubMed Google Scholar
Tehranchi, A. K. et al. The transcription factor DksA prevents conflicts between DNA replication and transcription machinery. Cell 141, 595–605 (2010).
Article CAS PubMed PubMed Central Google Scholar
Mirkin, E. V., Castro Roa, D., Nudler, E. & Mirkin, S. M. Transcription regulatory elements are punctuation marks for DNA replication. Proc. Natl Acad. Sci. USA 103, 7276–7281 (2006).
Article CAS PubMed PubMed Central ADS Google Scholar
Wiktor, J., Lesterlin, C., Sherratt, D. J. & Dekker, C. CRISPR-mediated control of the bacterial initiation of replication. Nucleic Acids Res. 44, 3801–3810 (2016).
Article CAS PubMed PubMed Central Google Scholar
Carl, P. L. Escherichia coli mutants with temperature-sensitive synthesis of DNA. Mol Gen Genet 109, 107–122 (1970).
Article CAS PubMed Google Scholar
Zhu, M., Mori, M., Hwa, T. & Dai, X. Disruption of transcription–translation coordination in Escherichia coli leads to premature transcriptional termination. Nat Microbiol 4, 2347–2356 (2019).
Article PubMed PubMed Central Google Scholar
Rosener, B. et al. Evolved bacterial resistance against fluoropyrimidines can lower chemotherapy impact in the Caenorhabditis elegans host. eLife 9, e59831 (2020).
Scott, T. A. et al. Host-microbe co-metabolism dictates cancer drug efficacy in C. elegans. Cell 169, 442–456.e18 (2017).
Article CAS PubMed PubMed Central Google Scholar
Rao, T. V. P. & Kuzminov, A. Exopolysaccharide defects cause hyper-thymineless death in Escherichia coli via massive loss of chromosomal DNA and cell lysis. Proc. Natl Acad. Sci. USA 117, 33549–33560 (2020).
Article PubMed PubMed Central ADS Google Scholar
Artsimovitch, I. & Landick, R. Pausing by bacterial RNA polymerase is mediated by mechanistically distinct classes of signals. Proc. Natl Acad. Sci. USA 97, 7090–7095 (2000).
Article CAS PubMed PubMed Central ADS Google Scholar
Giroux, X., Su, W.-L., Bredeche, M.-F. & Matic, I. Maladaptive DNA repair is the ultimate contributor to the death of trimethoprim-treated cells under aerobic and anaerobic conditions. Proc. Natl Acad. Sci. USA 114, 11512–11517 (2017).
Article CAS PubMed PubMed Central ADS Google Scholar
Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M. & Tanabe, M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 44, D457–D462 (2016).
Article CAS PubMed Google Scholar
Capra, J. A. & Singh, M. Predicting functionally important residues from sequence conservation. Bioinformatics 23, 1875–1882 (2007).
Article CAS PubMed Google Scholar
Miller, W. G. et al. The complete genome sequence and analysis of the Epsilonproteobacterium Arcobacter butzleri. PLoS ONE 2, e1358 (2007).
Article PubMed PubMed Central ADS Google Scholar
Svetlov, M. S., Cohen, S., Alsuhebany, N., Vázquez-Laslop, N. & Mankin, A. S. A long-distance rRNA base pair impacts the ability of macrolide antibiotics to kill bacteria. Proc. Natl Acad. Sci. USA 117, 1971–1975 (2020).
Article CAS PubMed PubMed Central ADS Google Scholar
Conrad, T. M. et al. RNA polymerase mutants found through adaptive evolution reprogram Escherichia coli for optimal growth in minimal media. Proc. Natl Acad. Sci. USA 107, 20500–20505 (2010).
Article CAS PubMed PubMed Central ADS Google Scholar
Cohen, S. S., Flaks, J. G., Barner, H. D., Loeb, M. R. & Lichtenstein, J. The mode of action of 5-fluorouracil and its derivatives. Proc. Natl Acad. Sci. USA 44, 1004–1012 (1958).
Article CAS PubMed PubMed Central ADS Google Scholar
Guegler, C. K. & Laub, M. T. Shutoff of host transcription triggers a toxin-antitoxin system to cleave phage RNA and abort infection. Mol Cell 81, 2361–2373.e9 (2021).
Article CAS PubMed PubMed Central Google Scholar
Chan, C. Y., Au-Yeang, C., Yew, W. W., Hui, M. & Cheng, A. F. Postantibiotic effects of antituberculosis agents alone and in combination. Antimicrob. Agents Chemother. 45, 3631–3634 (2001).
Article CAS PubMed PubMed Central Google Scholar
Malik, M. et al. Lethal synergy involving bicyclomycin: an approach for reviving old antibiotics. J. Antimicrob. Chemother. 69, 3227–3235 (2014).
Article CAS PubMed PubMed Central Google Scholar
Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal 17, 10–12 (2011).
Article Google Scholar
Datsenko, K. A. & Wanner, B. L. One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc. Natl Acad. Sci. USA 97, 6640–6645 (2000).
Article CAS PubMed PubMed Central ADS Google Scholar
Wiegand, I., Hilpert, K. & Hancock, R. E. W. Agar and broth dilution methods to determine the minimal inhibitory concentration (MIC) of antimicrobial substances. Nat. Protoc. 3, 163–175 (2008).
Article CAS PubMed Google Scholar
Svetlov, V. & Artsimovitch, I. Purification of bacterial RNA polymerase: tools and protocols. Methods Mol. Biol. 1276, 13–29 (2015).
Article CAS PubMed PubMed Central Google Scholar
Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359 (2012).
Article CAS PubMed PubMed Central Google Scholar
Ramírez, F. et al. deepTools2: a next generation web server for deep-sequencing data analysis. Nucleic Acids Res. 44, W160–W165 (2016).
Article PubMed PubMed Central Google Scholar
Sikorsky, J. A., Primerano, D. A., Fenger, T. W. & Denvir, J. DNA damage reduces Taq DNA polymerase fidelity and PCR amplification efficiency. Biochem. Biophys. Res. Commun. 355, 431–437 (2007).
Article CAS PubMed PubMed Central Google Scholar
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).
Article PubMed PubMed Central Google Scholar
Sievers, F. et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol. Syst. Biol. 7, 539 (2011).
Article PubMed PubMed Central Google Scholar
Download references
We thank P. Zappile from the NYU Genome Technology Center, which is partially supported by the Cancer Center Support grant no. P30CA016087, at the Laura and Isaac Perlmutter Cancer Center. We thank Y. Siu and D. Jones from the NYU Metabolomics Core Resource Laboratory. We thank K. Alzoubi for help with data visualization. This work was supported by the Blavatnik Family Foundation (E.N.); Howard Hughes Medical Institute (E.N.); NIH grant R01GM126891 (E.N.); NYU Medical Scientist Training Program (K.B.Y.); and Public Health Service Institutional Research Training Award no. T32 AI007180 (K.B.Y.).
Department of Biochemistry and Molecular Pharmacology, New York University Grossman School of Medicine, New York, NY, USA
Kevin B. Yang, Maria Cameranesi, Manjunath Gowder, Criseyda Martinez, Yosef Shamovsky, Vitaliy Epshtein, Zhitai Hao, Thao Nguyen, Eric Nirenstein, Ilya Shamovsky, Aviram Rasouly & Evgeny Nudler
Howard Hughes Medical Institute, New York University School of Medicine, New York, NY, USA
Aviram Rasouly & Evgeny Nudler
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
K.B.Y. and A.R. performed the MAGE experiments presented in Figs. 1–3 and Extended Data Figs. 1–3. K.B.Y. performed the data analysis and visualization of MAGE data and evolutionary analysis presented in Figs. 1–5 and Extended Data Figs. 1–3. K.B.Y., A.R. and M.C. prepared the strains and performed the validation experiments presented in Figs. 3 and 4 and Extended Data Figs. 5–9. C.M. performed the flow cytometry experiments presented in Fig. 3. M.G. performed SLR-qPCR and transcription kinetic experiments presented in Figs. 3 and 4 and Extended Data Fig. 7. Z.H. purified proteins and V.E. performed in vitro transcription assays presented in Extended Data Fig. 4. Y.S. performed the ChIP–seq experiments presented in Fig. 3. T.N. performed western blots presented in Extended Data Fig. 6. E. Nirenstein performed the MIC assays presented in Extended Data Fig. 9. I.S. performed next-generation sequencing. K.B.Y, A.R. and E. Nudler designed the project and wrote the manuscript with input from all authors. A.R. and E. Nudler supervised the project. All authors discussed the results.
Correspondence to Aviram Rasouly or Evgeny Nudler.
The authors declare no competing interests.
Nature thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
a-d, Representative heatmaps of relative mutant abundances from MAGE-seq screening under the following conditions: a, timepoints 0–2 (T0-2) and plate. T0 represents the initial MAGE-seq mutant pool. T1 and T2 are after subsequent 100 fold dilutions, allowing cells to reach the same OD before harvest. Plate represents MAGE-seq libraries prepared from mutants passaged on solid medium (LB agar). b, mutant abundance following screening with rifampicin. c, mutant abundance following screening with 15 (bcmlow) or 30 ug/mL (bcmhigh) of bicyclomycin. d, mutant abundance following screening with 1 (fulow) or 5 ug/mL (fuhigh) of 5-fluorouracil. Rows correspond with rpoB positions 510–537 and 563–572; columns correspond to all 20 amino acids along with stop codons. White squares denote wild-type residues. Each heatmap displays values from the mean of n = 3 biologically independent MAGE-seq experiments.
Source data
a, Fitness landscape of rpoB mutants comparing timepoint 1 to timepoint 0. b, Volcano plot of mutant pool fitness. Stop codons are in gray. Kernel density estimate plot of mutant fitness distribution is displayed above. c,d, Fitness landscape and volcano plot of rpoB mutants comparing those that form colonies on agar plates to timepoint 0. For (a) and (c) rows correspond with rpoB position and columns correspond with all 20 possible amino acid substitutions, grouped by property. White squares denote wild-type residues. Barplot of average positional fold change is presented right of the heatmap. All MAGE-seq experiments were conducted in n = 3 biologically independent replicates. P values are calculated using the independent t-test (two-sided). Error bars denote 95% confidence intervals.
Source data
a–d, Heatmap rows correspond with rpoB position and columns correspond with all 20 possible amino acid substitutions, grouped by property. White squares denote wild-type residues. Barplot of average positional fold change is presented right of the heatmap. e–h, Volcano plots of data in (a–d). All MAGE-seq experiments were conducted in n = 3 biologically independent replicates. P values are calculated using the independent t-test (two-sided). Error bars denote 95% confidence intervals.
Source data
a, MA plot of MAGE mutant pool sensitivity and resistance to rifampicin. Stop codons are in gray. b, Growth curves of hypersensitive mutant recovery after 1 hr of treatment with 80 ug/mL rifampicin (n = 2). c, MIC of hypersensitive strains (p < 1e-5, 1e-5 from left to right). An initial inoculum of 5 × 105 cells was incubated with increasing concentrations of rifampicin in Mueller-Hinton (MH) broth at 37C for 18 h. d, In vitro transcription recovery reaction of wild-type and rpoB L521Y RNAP following rifampicin treatment and wash. Location of the trimer and terminator (69nt) are denoted with arrows. Barplot of the intensity ratio of the terminator/trimer at each timepoint is provided below as the average of two independent experiments. e, Representative genome viewer of chromatin immunoprecipitation sequencing (ChIP-seq) signal from wild-type and rpoB T525D RNAP normalized to input. Cells were collected after 1 h recovery in LB from 1 h of treatment with 80 ug/mL rifampicin. Bar plots and growth curves represent the mean of n = 3 biologically independent experiments unless otherwise noted, with error bars denoting 95% confidence intervals. P values are calculated using the independent t-test (two-sided).
Source data
a, Viability in liquid of hypersensitive mutants following 1 h of treatment with rifampicin. Cells in exponential phase were treated with 80 ug/mL of rifampicin for 1 h at 37C before being washed, serially diluted, and measured for cell growth (n = 9, p = 3.1e-2, 3.1e-4 from left to right). b, Plating efficiency of hypersensitive mutants following 1 h of treatment with MIC normalized concentrations of rifampicin (p = 2.6e-3, 6.0e-4, 9.1e-4, 5.0e-4). Cells in exponential phase were treated with 4x or 10x their respective minimal inhibitory concentration (MIC) of rifampicin before being washed, serially diluted, and plated. c, Minimal bactericidal concentration (MBC) of hypersensitive strains (p < 1e-5, 1e-5). The MBC was determined as the lowest concentration of antibacterial agent that reduces the viability of the initial bacterial inoculum by ≥99.9% after 18 h at 37C. Rifampicin was tested up to the concentration of 80 ug/mL indicated by the asterisk. d, MBC/MIC ratio determined for rifampicin hypersensitive strains (p < 1e-5, 1e-5). e, Plating efficiency of hypersensitive mutants following 1 h of treatment with rifampicin, chloramphenicol (CM), and tetracycline (TET) (p = 6.2e-1, 3.3e-1, 5.5e-1, 7.1e-1). Cells in exponential phase were treated with 80ug/mL of rifampicin, 20 ug/mL of chloramphenicol, or 20ug/mL of tetracycline for 1 h at 37C before being washed, serially diluted, and plated on LB. f, Plating efficiency of hypersensitive mutants in recA or recBCD knockout strains after 1hr treatment with 20 ug/mL tetracycline at 37C. Cells in exponential phase were treated with 20 ug/mL of tetracycline for 1 h at 37C before being washed, serially diluted, and plated on LB. Bar plots represent the mean of n = 3 biologically independent experiments unless otherwise noted, with error bars denoting 95% confidence intervals. P values are calculated using the independent t-test (two-sided).
Source data
a, Relative expression of SOS response genes in response to 1hr treatment of 80 ug/mL rifampicin at 37C by RT-qPCR. Data was normalized to expression of rpsL and then between Rif and untreated conditions (n = 5, p = 4.7e-4, 3.7e-3, 1.8e-2, 7.5e-3 from left to right). b, LexA relative protein levels measured by western blot in exponentially growing hypersensitive mutants following 1 hr treatment of 80ug/mL rifampicin at 37C (p = 2.5e-2, 3.2e-2, 5.4e-3). Data is normalized to rpoB relative protein levels. c, Representative western blot of data in (b). d, Reactive oxygen species (ROS) of cells treated with ampicillin or rifampicin, measured by carboxy-H2DCFDA (p = 5.8e-6). Exponentially growing cells were treated for 1hr with 80 ug/mL rifampicin or 10 ug/mL ampicillin at 37C. Carboxy-H2DCFDA was added 30 min prior to sample collection. e, Plating efficiency of wild-type cells treated with 10 ug/mL ampicillin with or without the addition of 500 uM 2,2′ dipyridyl and 150 mM thiourea (n = 5; p = 3.9e-4). f, Plating efficiency of hypersensitive mutants following 1 h of treatment with 80 ug/mL rifampicin with or without the addition of 500 uM 2,2′ dipyridyl and 150 mM thiourea (p = 1.5e-4, 1.7e-4, 3.8e-5, 1.7e-4). Bar plots represent the mean of n = 3 biologically independent experiments unless otherwise noted, with error bars denoting 95% confidence intervals. P values are calculated using the independent t-test (two-sided).
Source data
a, Schematic of semi long range quantitative PCR (SLR-qPCR) used to detect double strand breaks (DSB). Two sets of primers, yielding a short and a long amplicon provide data representing the total amount of template (an internal normalization control) and undamaged DNA respectively. DSB density is measured using genomic DNA from untreated cells as reference. DNA damage is calculated as DSBs per 10 kb DNA. b, Schematic of SLR-qPCR primer design around the promoter and gene-body of selected genes. c, Schematic for replication block experiments, a dual plasmid system expressing dCas9 and a guide RNA targeting the oriC was introduced to cells of interest. dCas9 prevents binding of replication initiation proteins. dCas9 was induced with 200 ng/mL of aTc for 30 min prior to adding Rif. d,e, Number of DSB per 10 kilobases at the fepA and ilvC promoter measured using SLR-qPCR following 1 hr treatment with 80ug/mL rifampicin with or without replication (p = 1.5e-2, 7.6e-3, 9.9e-4, 3.0e-4, 2.3e-3, 8.6e-3 from left to right). Replication was blocked using dCas9 with a guide RNA targeting the oriC for 30 min prior to adding rifampicin. f, qPCR measurement of the oriC/ter ratio of hypersensitive mutants in the PC2 background before and after shifting cells to 42C, normalized to actively replicating wild-type cells at 30C. Cells with a lower oriC/ter ratio are replicating less (n = 4; p = 9.3e-4, 1.9e-4, 4.8e-4). g, Plating efficiency of hypersensitive mutants in PC2 at the replication permissive temperature 30C (p = 2.2e-4, 5.1e-5). Bar plots represent the mean of n = 3 biologically independent experiments unless otherwise noted, with error bars denoting 95% confidence intervals. P values are calculated using the independent t-test (two-sided).
Source data
a, Clustering of RNA expression counts of rpoB mutants with increased, decreased, or unchanged sensitivity towards bicyclomycin and 5-fluouracil. Cells were grown in LB and samples were collected in n = 3 biologically independent replicates. Volcano plots from RNAseq analysis is provided in Supplementary Fig. 2. b, GO enrichment analysis of differentially regulated genes in rifampicin binding site mutants as compared to wild-type E. coli. P values are calculated using the one-sided Fishers exact test with Benjamini-Hochberg correction.
Source data
a, Differential gene expression of pyrimidine de novo biosynthesis, ribonucleotide salvage, or deoxyribonucleotide salvage pathways in fast mutants rpoB 515Y and 527W compared to wild-type, as measured by RNAseq. Statistically significant differences (p < 01) are shown in blue. RNAseq was performed in n = 3 biologically independent replicates and two-sided p values were calculated using the Wald test with Benjamini-Hochberg correction. b,c, Plating efficiency of 5-fluorouracil sensitive mutants in a background of pyrD or upp. Cells from overnight cultures were serially diluted and plated on LB agar plates with the indicated concentrations of 5FU in ug/mL. d, Purines and dTDP in 5-fluorouracil sensitive rpoB mutants as detected by LCMS (n = 5; p = 4.5e-5, 5.4e-1, 6.4e-3, 7.1e-3, 6.8e-1, 7.5e-2, 1.1e-3, 3.9e-1, 8.8e-1 from left to right). e, Growth curves of rpoB 515Y compared to wild-type cells with the addition of 0, .1, or 1ug/mL of thymidine. Values represent the mean of 2 biological replicates. f, Plating efficiency of 5-fluorouracil sensitive mutants in a thyA background following thymidine starvation (p = 3.2e-3, 4.2e-3, 9.8e-3). g, Plating efficiency of rpoB 513P with the indicated concentrations of 5FU in ug/mL (p = 7.2e-6). h, Plating efficiency of rpoB 513P in a thyA background following thymidine starvation (p = 5.2e-4). Barplots and growth curves represent the mean of n = 3 biologically independent experiments unless otherwise noted, with error bars denoting 95% confidence intervals. P values are calculated using the independent t-test (two-sided).
Source data
a, Growth curves of 5-fluorouracil sensitive mutants grown at 25 °C. b, Growth curves of rpoB mutants stratified by 5-fluourouracil sensitivity at 25C and 42C. c, Growth curves of wild-type cells with subinhibitory bicyclomycin, grown at 25 °C. Growth curves represent the mean of n = 3 biologically independent experiments unless noted otherwise, with error bars denoting 95% confidence intervals. d, Scatterplots of evolutionary conservation and average positional fitness, resistance, or sensitivity for each residue position. Linear regression lines and 95% confidence intervals are overlayed.
Source data
Supplementary Figs. 1 and 2, containing the gating strategy for Fig. 3 and volcano plots for RNAseq data related to Extended Data Fig. 8, and Supplementary Tables 1–5, which list the strains and primers used in the study.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Yang, K.B., Cameranesi, M., Gowder, M. et al. High-resolution landscape of an antibiotic binding site. Nature (2023). https://doi.org/10.1038/s41586-023-06495-6
Download citation
Received: 03 August 2022
Accepted: 28 July 2023
Published: 30 August 2023
DOI: https://doi.org/10.1038/s41586-023-06495-6
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.