
MBW:Senescence Can Explain Microbial PersistenceFrom MathBioArticle: Senescence can explain microbial persistence Authors: I. Klapper, P. Gilbert, B. P. Ayati, J. Dockery and P. S. Stewart Reviewed by: Ryan Brown ContentsExecutive SummaryMany antimicrobial medications like penicillin have been observed to significantly reduce bacterial populations, though they are unable to completely eliminate them. Regardless of the exposure time, some bacterial cells, called persisters, manage to survive. After exposure to the antimicrobial has subsided the persisters are able to reproduce intolerant daughter cells. Persister cells possess many interesting qualities and the mechanism of persister cell formation is of interest. Most theories suggest that the bacteria are able to switch between levels of persister states depending on environmental conditions. In their article Klapper et al. introduce a mechanism where cells senesce with time and this increase in senescence results in increased tolerance to the drug. Using this idea as a foundation they develop a mathematical model to simulate populations of persisters vs time in a culture. IntroductionAfter a population of bacteria is exposed to some antimicrobial any bacteria that survive are called persisters. It is important to discern the difference between resistance and persistence. The persister cells do not have an immunology to the antimicrobial medications but rather possess the capability to tolerate exposure by not dying. These persister cells have many interesting characteristics including:
The basis of the model relies on previous work observations of microbial senescence. Senescence can be thought of as maturity or old age and it is believed that with increased senescence, mother cells display different behaviors than their daughters. Classical views of cell division describe a parent cell dividing into two identical daughters with no age associated with them. Certain yeasts like Saccharomyces cerevisiae have assymetric cell division (physically larger mother cell) and have been the focus of studies of microbial aging ^{[2]}. Stewart et al. were able to track aging in the morphologically symmetric dividers E. Coli and concluded that it is a result of a physiological asymmetry in the cells^{[2]}. They tracked roughly 35,000 E. coli to show that the cells had an increase in senescence in the form of slower growth rate. Their model assumed that the cells contained poles of different ages and this age was increased with each division. See Figure A^{[2]} The purpose Klapper et al.'s article is to create a mathematical model to describe the behaviors of persisters using the concept of senescence. They make several assumptions to develop their model. Below they are listed straight from the original article^{[3]}.
Assumptions 13 are essential to the model, while 47 are for completeness of the model and used for simulation purposes. Older cells are regarded as the persisters and grow more slowly than the younger, non persister cells. New cells are produced from older cells as shown in Figure 1^{[3]}. .
TheoryMicrobial ratesThe Table in the section below describes all important variables and parameters that will be used in the following discussion. Populations are given in units of c.f.u. They start by defining the substrate usage per c.f.u at time t and cell age a, , as where is the growth media concentration, is the substrate usage factor, is the minimum substrate usage factor, is what they refer to as the senescence time parameter. They then define a birth rate as where is the cell yeild coefficient. Notice that with increasing age the cell will use less substrate and therefore yield less offspring as asserted by assumption 4. Finally, by assumption 5 the death rate is with being the death rate constant. Age StructureThe differential equation used to describe the bacterial population density, , is This equation is not valid for b(a=0,t) i.e. for the cells that are born from the existing population at a given time t. To account for this the authors show that Similarly, as these new bacteria are created more substrate is used so there is a changing concentration given by the following At this point we can introduce the senescence factor, s(a) and the senescence weighted total bacteria population B(t) Now equations (24) are rewritten as The initial conditions chosen by the authors as follows
AntimicrobialTo simulate the exposure to an antimicrobial Klapper et al. modify equation (5) with a killing rate, , that only applies to young and therefore non persistent cells. If is a killing rate coefficient and the antimicrobial concentration then equation (5) becomes with The tolerance coefficient parameter, , is a value that corresponds to how old a cell has to be before it is considered persistent. For small values of a cell becomes persistent at a relatively young age. This parameter is not experimentally found but is estimated by the authors. They discuss their choice its value in the Methods section of the original article. Notice also that they chose tolerance age to have a linear relationship with antimicrobial concentration. This choice is consistent with the linear relation that senescence has with time and the authors explain that this is for definiteness and that some other relationship can be used as long as senescence and tolerance age are consistent with each other. Equation 8 is numerically solved in the simulation with the initial conditions and parameter values given in the table below. Table of VariablesNote: This is a copy of Table 1 from the original article. Deciding values of each parameter are discussed in the Methods section of the article. The important constraints are:
ResultsUsing numerical methods^{[4]} equations (6)(8) are solved for population vs time. Figure 2 shows the result. The antimicrobial is applied to the batch during different phases of the bacteria growth. After some time of exposure it is removed, the sample is recultured (simulated) and the cells are allowed to repopulate. In Figure (a) the antimicrobial is applied during a stationary growth phase while in (b) it is applied during the exponential growth phase. The plots are loglog plots hence the linear relation in (b) for t<5h. Both graphs tell similar stories. The overall population of cells increases rapidly in the initial stages before the introduction of the antimicrobial. Notice that while the overall population increases dramatically the amount of persister cells does not follow the same sharp increase. When the antimicrobial is applied there is a drastic decline in the population of all cells (solid line), to the point where all cells present are persister cells (dashed line). The persisters then create several daughter cells and this causes the population to shoot back up. The persister cell population does not grow significantly until the stationary phase when the total cell population begins to level out. This is a consequence of the time delay (dependant on the value of ) when the daughter cells age to persisters. It appears at large time all cells will become persisters, but will likely die to limited resources and/or natural causes. Extensions
ConclusionsKnowing the mechanism that allows for the presence of persisters after exposure to an antimicrobial has obvious implications to medicine. Understanding the mechanism that causes persistence could lead to the development of more powerful antibiotics. The idea that cells develop persistence as a consequence of senescence might be a step toward a better understanding. With the senescence model, Klapper et al. have been able to mathematically simulate some of the characteristics observed by persisters.
References
