skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Whole-cell modeling of E. coli colonies enables quantification of single-cell heterogeneity in antibiotic responses
Antibiotic resistance poses mounting risks to human health, as current antibiotics are losing efficacy against increasingly resistant pathogenic bacteria. Of particular concern is the emergence of multidrug-resistant strains, which has been rapid among Gram-negative bacteria such asEscherichia coli. A large body of work has established that antibiotic resistance mechanisms depend on phenotypic heterogeneity, which may be mediated by stochastic expression of antibiotic resistance genes. The link between such molecular-level expression and the population levels that result is complex and multi-scale. Therefore, to better understand antibiotic resistance, what is needed are new mechanistic models that reflect single-cell phenotypic dynamics together with population-level heterogeneity, as an integrated whole. In this work, we sought to bridge single-cell and population-scale modeling by building upon our previous experience in “whole-cell” modeling, an approach which integrates mathematical and mechanistic descriptions of biological processes to recapitulate the experimentally observed behaviors of entire cells. To extend whole-cell modeling to the “whole-colony” scale, we embedded multiple instances of a whole-cellE.colimodel within a model of a dynamic spatial environment, allowing us to run large, parallelized simulations on the cloud that contained all the molecular detail of the previous whole-cell model and many interactive effects of a colony growing in a shared environment. The resulting simulations were used to explore the response ofE.colito two antibiotics with different mechanisms of action, tetracycline and ampicillin, enabling us to identify sub-generationally-expressed genes, such as the beta-lactamase ampC, which contributed greatly to dramatic cellular differences in steady-state periplasmic ampicillin and was a significant factor in determining cell survival.  more » « less
Award ID(s):
2019589
PAR ID:
10559565
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Editor(s):
Beard, Daniel A
Publisher / Repository:
PLOS Journals
Date Published:
Journal Name:
PLOS Computational Biology
Volume:
19
Issue:
6
ISSN:
1553-7358
Page Range / eLocation ID:
e1011232
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Antibiotic responses in bacteria are highly dynamic and heterogeneous, with sudden exposure of bacterial colonies to high drug doses resulting in the coexistence of recovered and arrested cells. The dynamics of the response is determined by regulatory circuits controlling the expression of resistance genes, which are in turn modulated by the drug’s action on cell growth and metabolism. Despite advances in understanding gene regulation at the molecular level, we still lack a framework to describe how feedback mechanisms resulting from the interdependence between expression of resistance and cell metabolism can amplify naturally occurring noise and create heterogeneity at the population level. To understand how this interplay affects cell survival upon exposure, we constructed a mathematical model of the dynamics of antibiotic responses that links metabolism and regulation of gene expression, based on the tetracycline resistancetetoperon inE. coli. We use this model to interpret measurements of growth and expression of resistance in microfluidic experiments, both in single cells and in biofilms. We also implemented a stochastic model of the drug response, to show that exposure to high drug levels results in large variations of recovery times and heterogeneity at the population level. We show that stochasticity is important to determine how nutrient quality affects cell survival during exposure to high drug concentrations. A quantitative description of how microbes respond to antibiotics in dynamical environments is crucial to understand population-level behaviors such as biofilms and pathogenesis. 
    more » « less
  2. Abstract Antibiotic treatment kills a large portion of a population, while a small, tolerant subpopulation survives. Tolerant bacteria disrupt antibiotic efficacy and increase the likelihood that a population gains antibiotic resistance, a growing health concern. We examined howE. colitranscriptional networks changed in response to lethal ampicillin concentrations. We are the first to apply transcriptional regulatory network (TRN) analysis to antibiotic tolerance by leveraging existing knowledge and our transcriptional data. TRN analysis shows that gene expression changes specific to ampicillin treatment are likely caused by specific sigma and transcription factors typically regulated by proteolysis. These results demonstrate that to survive lethal concentration of ampicillin specific regulatory proteins change activity and cause a coordinated transcriptional response that leverages multiple gene systems. 
    more » « less
  3. Abstract While antibiotic resistance poses a threat from both Gram-positive bacteria (GPB) and Gram-negative bacteria (GNB), GNB pose a more imminent public health hazard globally. GNB are a threat to growing antibiotic resistance because of the complex makeup of the membrane. The AcrAB-TolC efflux pump is a known resistance mechanism ofEscherichia coli (E. coli)cells. This study utilized molecular dynamics modeling to visualize some of the changes occurring at a molecular level when airborne bacteria are exposed to stress and antibiotics. This study was conducted to build upon previous experimental research showing that there is an increase in antibiotic resistance and efflux pump activity when exposed to aerosolization. AcrB and AcrAB-TolC proteins were simulated under standard and increased pressure to compare the effect of aerosolization on the binding to the three different antibiotics (puromycin (PUY), ampicillin (AMP) and sulfamethoxazole-trimethoprim (SXT)) to the AcrB binding site. Analysis such as root-mean-square deviation of atomic positions and root-mean-square fluctuation, the opening of TolC, and the significant molecular mechanics with generalized Born and surface area solvation (MM-GBSA) scores associated with specific ligands were recorded. Resistance in experimental data indicated a relationship between the docking scores and some ligand–protein interactions. Results showed that there was more flexibility in the proteins within simulations conducted under standard pressure for the AcrB protein and the full tripartite complex AcrAB-TolC, showing that increased pressure causes more rigidity. MM-GBSA scores, used to calculate the free energy of ligand–protein binding, did not show a significant change, but interestingly, the strongest MM-GBSA scores were for ligands that moved to another binding pocket and did not result in resistance or opening of the efflux pump. However, the ligand moved from the binding site and did not cause the opening of TolC to increase significantly, whereas PUY and AMP were bound to the binding site for the duration of all simulations. AMP ligands under increased pressure showed the largest change in opening of the TolC efflux pump and aligns with experimental data showingE. colicells had the most resistance to AMP after aerosolization. These results, in addition to other real-time changes such as OM proteins and mutations of targets within the cell, could be used to delineate and mitigate antibiotic resistance mechanisms. 
    more » « less
  4. Abstract Proliferation of multidrug-resistant (MDR) bacteria poses a threat to human health, requiring new strategies. Here we propose using fitness neutral gene expression perturbations to potentiate antibiotics. We systematically explored 270 gene knockout-antibiotic combinations inEscherichia coli, identifying 90 synergistic interactions. Identified gene targets were subsequently tested for antibiotic synergy on the transcriptomic level via multiplexed CRISPR-dCas9 and showed successful sensitization ofE. coliwithout a separate fitness cost. These fitness neutral gene perturbations worked as co-therapies in reducing aSalmonella entericaintracellular infection in HeLa. Finally, these results informed the design of four antisense peptide nucleic acid (PNA) co-therapies,csgD,fnr,recAandacrA, against four MDR, clinically isolated bacteria. PNA combined with sub-minimal inhibitory concentrations of trimethoprim against two isolates ofKlebsiella pneumoniaeandE. colishowed three cases of re-sensitization with minimal fitness impacts. Our results highlight a promising approach for extending the utility of current antibiotics. 
    more » « less
  5. Pagliara, Stefano (Ed.)
    ABSTRACT Non-genetic factors can cause significant fluctuations in gene expression levels. Regardless of growing in a stable environment, this fluctuation leads to cell-to-cell variability in an isogenic population. This phenotypic heterogeneity allows a tiny subset of bacterial cells in a population called persister cells to tolerate long-term lethal antibiotic effects by entering into a non-dividing, metabolically repressed state. We occasionally noticed a high variation in persister levels, and to explore this, we tested clonal populations starting from a single cell using a modified Luria-Delbrück fluctuation test. Although we kept the conditions same, the diversity in persistence level among clones was relatively consistent: varying from ~60- to 100- and ~40- to 70-fold for ampicillin and apramycin, respectively. Then, we divided and diluted each clone to observe whether the same clone had comparable persister levels for more than one generation. Replicates had similar persister levels even when clones were divided, diluted by 1:20, and allowed to grow for approximately five generations. This result explicitly shows a cellular memory passed on for generations and eventually lost when cells are diluted to 1:100 and regrown (>seven generations). Our result demonstrates (1) the existence of a small population prepared for stress (“primed cells”) resulting in higher persister numbers; (2) the primed memory state is reproducible and transient, passed down for generations but eventually lost; and (3) a heterogeneous persister population is a result of a transiently primed reversible cell state and not due to a pre-existing genetic mutation. IMPORTANCEAntibiotics have been highly effective in treating lethal infectious diseases for almost a century. However, the increasing threat of antibiotic resistance is again causing these diseases to become life-threatening. The longer a bacteria can survive antibiotics, the more likely it is to develop resistance. Complicating matters is that non-genetic factors can allow bacterial cells with identical DNA to gain transient resistance (also known as persistence). Here, we show that a small fraction of the bacterial population called primed cells can pass down non-genetic information (“memory”) to their offspring, enabling them to survive lethal antibiotics for a long time. However, this memory is eventually lost. These results demonstrate how bacteria can leverage differences among genetically identical cells formed through non-genetic factors to form primed cells with a selective advantage to survive antibiotics. 
    more » « less