Growth-mediated negative feedback shapes quantitative antibiotic response

Dose–response relationships are a general concept for quantitatively describing biological systems across multiple scales, from the molecular to the whole‐cell level. A clinically relevant example is the bacterial growth response to antibiotics, which is routinely characterized by dose–response curves. The shape of the dose–response curve varies drastically between antibiotics and plays a key role in treatment, drug interactions, and resistance evolution. However, the mechanisms shaping the dose–response curve remain largely unclear. Here, we show in Escherichia coli that the distinctively shallow dose–response curve of the antibiotic trimethoprim is caused by a negative growth‐mediated feedback loop: Trimethoprim slows growth, which in turn weakens the effect of this antibiotic. At the molecular level, this feedback is caused by the upregulation of the drug target dihydrofolate reductase (FolA/DHFR). We show that this upregulation is not a specific response to trimethoprim but follows a universal trend line that depends primarily on the growth rate, irrespective of its cause. Rewiring the feedback loop alters the dose–response curve in a predictable manner, which we corroborate using a mathematical model of cellular resource allocation and growth. Our results indicate that growth‐mediated feedback loops may shape drug responses more generally and could be exploited to design evolutionary traps that enable selection against drug resistance.

Quantifying the Determinants of Evolutionary Dynamics Leading to Drug Resistance

The emergence of drug resistant pathogens is a serious public health problem. It is a long-standing goal to predict rates of resistance evolution and design optimal treatment strategies accordingly. To this end, it is crucial to reveal the underlying causes of drug-specific differences in the evolutionary dynamics leading to resistance. However, it remains largely unknown why the rates of resistance evolution via spontaneous mutations and the diversity of mutational paths vary substantially between drugs. Here we comprehensively quantify the distribution of fitness effects (DFE) of mutations, a key determinant of evolutionary dynamics, in the presence of eight antibiotics representing the main modes of action. Using precise high-throughput fitness measurements for genome-wide Escherichia coli gene deletion strains, we find that the width of the DFE varies dramatically between antibiotics and, contrary to conventional wisdom, for some drugs the DFE width is lower than in the absence of stress. We show that this previously underappreciated divergence in DFE width among antibiotics is largely caused by their distinct drug-specific dose-response characteristics. Unlike the DFE, the magnitude of the changes in tolerated drug concentration resulting from genome-wide mutations is similar for most drugs but exceptionally small for the antibiotic nitrofurantoin, i.e., mutations generally have considerably smaller resistance effects for nitrofurantoin than for other drugs. A population genetics model predicts that resistance evolution for drugs with this property is severely limited and confined to reproducible mutational paths. We tested this prediction in laboratory evolution experiments using the "morbidostat", a device for evolving bacteria in well-controlled drug environments. Nitrofurantoin resistance indeed evolved extremely slowly via reproducible mutations-an almost paradoxical behavior since this drug causes DNA damage and increases the mutation rate. Overall, we identified novel quantitative characteristics of the evolutionary landscape that provide the conceptual foundation for predicting the dynamics of drug resistance evolution.

Systematic discovery of drug interaction mechanisms

Drug combinations are increasingly important in disease treatments, for combating drug resistance, and for elucidating fundamental relationships in cell physiology. When drugs are combined, their individual effects on cells may be amplified or weakened. Such drug interactions are crucial for treatment efficacy, but their underlying mechanisms remain largely unknown. To uncover the causes of drug interactions, we developed a systematic approach based on precise quantification of the individual and joint effects of antibiotics on growth of genome-wide Escherichia coli gene deletion strains. We found that drug interactions between antibiotics representing the main modes of action are highly robust to genetic perturbation. This robustness is encapsulated in a general principle of bacterial growth, which enables the quantitative prediction of mutant growth rates under drug combinations. Rare violations of this principle exposed recurring cellular functions controlling drug interactions. In particular, we found that polysaccharide and ATP synthesis control multiple drug interactions with previously unexplained mechanisms, and small molecule adjuvants targeting these functions synthetically reshape drug interactions in predictable ways. These results provide a new conceptual framework for the design of multidrug combinations and suggest that there are universal mechanisms at the heart of most drug interactions.

Functional coupling between HIV-1 integrase and the SWI/SNF chromatin remodeling complex for efficient in vitro integration into stable nucleosomes

Establishment of stable HIV-1 infection requires the efficient integration of the retroviral genome into the host DNA. The molecular mechanism underlying the control of this process by the chromatin structure has not yet been elucidated. We show here that stably associated nucleosomes strongly inhibit in vitro two viral-end integration by decreasing the accessibility of DNA to integrase. Remodeling of the chromatinized template by the SWI/SNF complex, whose INI1 major component interacts with IN, restores and redirects the full-site integration into the stable nucleosome region. These effects are not observed after remodeling by other human remodeling factors such as SNF2H or BRG1 lacking the integrase binding protein INI1. This suggests that the restoration process depends on the direct interaction between IN and the whole SWI/SNF complex, supporting a functional coupling between the remodeling and integration complexes. Furthermore, in silico comparison between more than 40,000 non-redundant cellular integration sites selected from literature and nucleosome occupancy predictions also supports that HIV-1 integration is promoted in the genomic region of weaker intrinsic nucleosome density in the infected cell. Our data indicate that some chromatin structures can be refractory for integration and that coupling between nucleosome remodeling and HIV-1 integration is required to overcome this natural barrier.

Thermodynamics of intragenic nucleosome ordering

The nucleosome ordering observed in vivo along yeast genes is described by a thermodynamical model of nonuniform fluid of 1D hard rods confined by two excluding energy barriers at gene extremities. For interbarrier distances L less than or approximately equal to 1.5 kbp, nucleosomes equilibrate into a crystal-like configuration with a nucleosome repeat length (NRL) L/n approximately 165 bp, where n is the number of regularly positioned nucleosomes. We also observe "bistable" genes with a fuzzy chromatin resulting from a statistical mixing of two crystal states, one with an expanded chromatin (NRL approximately L/n) and the other with a compact one (NRL approximately L/(n+1)). By means of single nucleosome switching, bistable genes may drastically alter their expression level as suggested by their higher transcriptional plasticity. These results enlighten the role of the intragenic chromatin on gene expression regulation.

A novel strategy of transcription regulation by intragenic nucleosome ordering

Numerous studies of chromatin structure showed that nucleosome free regions (NFRs) located at 5' gene ends contribute to transcription initiation regulation. Here, we determine the role of intragenic chromatin structure on gene expression regulation. We show that, along Saccharomyces cerevisiae genes, nucleosomes are highly organized following two types of architecture that depend only on the distance between the NFRs located at the 5' and 3' gene ends. In the first type, this distance constrains in vivo the positioning of n nucleosomes regularly organized in a "crystal-like" array. In the second type, this distance is such that the corresponding genes can accommodate either n or (n + 1) nucleosomes, thereby displaying two possible crystal-like arrays of n weakly compacted or n + 1 highly compacted nucleosomes. This adaptability confers "bi-stable" properties to chromatin and is a key to its dynamics. Compared to crystal-like genes, bi-stable genes present higher transcriptional plasticity, higher sensitivity to chromatin regulators, higher H3 turnover rate, and lower H2A.Z enrichment. The results strongly suggest that transcription elongation is facilitated by higher chromatin compaction. The data allow us to propose a new paradigm of transcriptional control mediated by the stability and the level of compaction of the intragenic chromatin architecture and open new ways for investigating eukaryotic gene expression regulation.

Open chromatin encoded in DNA sequence is the signature of 'master' replication origins in human cells

For years, progress in elucidating the mechanisms underlying replication initiation and its coupling to transcriptional activities and to local chromatin structure has been hampered by the small number (approximately 30) of well-established origins in the human genome and more generally in mammalian genomes. Recent in silico studies of compositional strand asymmetries revealed a high level of organization of human genes around 1000 putative replication origins. Here, by comparing with recently experimentally identified replication origins, we provide further support that these putative origins are active in vivo. We show that regions approximately 300-kb wide surrounding most of these putative replication origins that replicate early in the S phase are hypersensitive to DNase I cleavage, hypomethylated and present a significant enrichment in genomic energy barriers that impair nucleosome formation (nucleosome-free regions). This suggests that these putative replication origins are specified by an open chromatin structure favored by the DNA sequence. We discuss how this distinctive attribute makes these origins, further qualified as 'master' replication origins, priviledged loci for future research to decipher the human spatio-temporal replication program. Finally, we argue that these 'master' origins are likely to play a key role in genome dynamics during evolution and in pathological situations.