Team:Goettingen/Model

Modelling of glyphosate-dependent growth inhibition

B. subtilis as a glyphosate detector

As described in the Background section, glyphosate specifically inhibits enolpyruvyl-shikimate-3-phosphate (EPSP) synthase, an enzyme present in plants, fungi, bacteria and archaea. The EPSP synthase converts shikimate-3-phosphate (S3P) and phosphoenolpyruvate (PEP) to EPSP, an intermediate in the shikimate pathway and required for synthesis of aromatic amino acids. Interaction of glyphosate with EPSP leads to depletion of cellular levels of aromatic amino acids resulting in growth inhibition. We have observed that growth of B. subtilis is strongly inhibited by glyphosate (50% reduction of the growth rate at 1 mM glyphosate, see Results section). To generate a glyphosate detection system, we wanted to make use of the glyphosate sensitivity of B. subtilis. However, the wild type allows to detect glyphosate only in the range between 0 mM and 2 mM. Fortunately, we have discovered the first glyphosate transporters of a living organism (see Results section). So we could make use of the B. subtilis ΔgltT ΔgltP mutant strain lacking both glyphosate transporters. The growth rate of this strain is reduced by 50% at a glyphosate concentration of about 8 mM! Therefore, the transporter mutant allows us develop a system with an increased detection range between 8 mM and 30 mM. To employ B. subtilis as a detector for glyphosate, we had to quantify the glyphosate-dependent reduction of the growth rate of the different B. subtilis strains. The growth rate seems to be ideal for this application because it is easy to reproduce. The relation between the growth rate and the glyphosate concentration can be obtained by modeling our growth curves with the bacterial growth model:

The bacterial growth model

Under lab conditions, the right environmental parameters and nutrients provided, bacterial growth can be observed in form of growth curves, plotting the bacterial mass as a function of time. A convenient way to observe those growth curves is provided by measuring the optical density also referred to as the turbidity.
During the growth of bacteria three distinct phases can be differentiated (Figure 1), all of which can be associated with physiological changes in the cell cultures: (1)

  1. During lag phase the bacteria are suggested to adapt to physiological changes in the medium. Adaption can be achieved by e.g. changes in the proteome of the bacteria in order to meet the needs for the new culture requirements. Bacterial growth is strongly reduced
  2. After adaption bacterial growth enters the log phase alternatively named exponential phase. During this phase the bacteria show exponential growth, the most rapid growth possible under given conditions. The actual double-timing of a certain strain is determined by the growth rate, the increase of cell mass per time.
  3. In the stationary phase no net growth can be observed. Nutrients and carbon sources are used up and increase in cell mass is no longer possible. Dead cells can lyse and provide another source of nutrients, explaining a small increase in cell mass that sometimes still occurs.

To be scientifically correct a fourth phase can be observed in the bacterial growth curve namely the death phase. As we measured our bacterial mass in terms of optical density this phase couldn't be observed in our experiments and has no relevance for our modeling approaches.

Correlation between growth rate and glyphosate concentration

Obtaining the growth rate μ

During their growth bacteria reach a phase of exponential growth as described above. In this phase every bacteria in the culture divides at the fastest rate possible, therefore the growth of the colony is directly dependend on the amount of cells already grown. This is also true for the optical density as it's proportional to cell mass. To put this growth into a context that can be calculated, a constant of proportionality is needed. This constant can be described as the growth rate of the colony, it is specific to strains and also dependend on environmental influences like glyphosate. These dependencies can be described by the following formula:

This short differential equation describes the growth behaviour of the bacteria during the log phase. The equation can be solved by first separating the derivates:

Next the equation has to be integrated at both sides yielding to the following form:

Inserting the initial condition t=0 and ODt=0 results in:

From this solution a very easy and convenient way to obtain the growth rate for every growth curve arrives. As it is probably known by everyone working in the field of microbiology plotting the logarithm of cell mass (or OD) against time leads to a linear section in the affiliated graph. This linear section is exactly what is described by the last formula, the logarithm of the exponential phase of bacterial growth with the growth rate as its slope. To reproducibly calculate the slope/growth rate we wrote a python script, that enables us to select the linear section of the graph by hand. After that a linear fit is calculated and the resulting growth rate is safed. In Figure X usage of the program is shown.

wild type

double mutant (ΔgltT and ΔgltP)

Figure X. Usage of the python script to obtain the growth rate from the linear section of growth curves.

Describing the dependency of growth rate from glyphosate concentration

After assessing the growth rates by doing fits of the growth curves, we plotted them against the glyphosate concentration in the environment of every growth experiment respectively. The graphs can be observed in Figure X and X:

Figure 2. Dependency of growth rate to glyphosate concentration in the wild type strain.

Figure 3. Dependency of growth rate to glyphosate concentration in the double mutant strain (ΔgltT and ΔgltP).

As stated in the Results section of our project the growth rate is strongly dependend on glyphosate concentration. While low concentrations of glyphosate (relative to the concentration of 50% growth reduction for every strain) do not strongly affect the growth behavior , increasing the concentration yields to very strong reduction. Moreover, higher doses of glyphosate do not change the growth rate very strongly anymore. The graphs strongly resemble the well-known Hill equation (see below) suggesting cooperativity in either the transport of glyphosate or interaction with the EPSP synthase.

We fitted the Hill equation to the data using a python script and the leastsq method of the scipy library, the fitted data are shown in Figure X and X.

Figure X. Least square fit of the dependency between growth rate and glyphosate concentration for the wild type using the hill equation. In the right corner, the resulting function with determined parameters is shown.

Figure X. Least square fit of the dependency between growth rate and glyphosate concentration for the double mutant (ΔgltT and ΔgltP) using the hill equation. In the right corner, the resulting function with determined parameters is shown.

The hill equation can be transposed to show growth rate as dependency for glyphosate concentration:

With this transposition it is possible to determine the glyphosate concentration by evaluating the growth behaviour of the strains. For every strain these modelling parameters are different, for our double mutant and the wild type we already found the right fit (the values for the parameters can be found in the two diagrams above), enabling us to detect glyphosate from 0.5 up to 50 mM.

HIER FEHLT NOCH EIN EINSETZEN DER WACHSTUMSGESCHWINDIGKEIT IN DIE GLEICHUNG UM DAMIT EINEN GLYPHOSAT-WERT AUSZURECHNEN!

But there is one catch, that made us try to further extend our consideration. We can't be exactly sure, that there are no other environmental conditions that influence our growth experiments except for glyphosate. At last correlation must not mean causality. To ensure that really the glyphosate is the dominating factor in our experiments, we set out a competition assay between the wild type and the double mutant, gradually increasing the concentration of the herbicide in the media.

Modeling the competition for different glyphosate concentrations

Modelling the optical density

During a competition assay both strains should be influenced by glyphosate in the same way they would while growing isolated from each other. The growth rate and it's reduction by glyphosate can be assumed to follow the hill equations we already determined for both strains. The difference to a normal growing experiment can be found in the competition between both strains for nutrients. Like in an isolated growing experiment, the growth curve can be plotted as logistic growth, but the limiting threshhold for the growth is not only influenced by the cell mass of one single strain but by both strains together. On top of this theory two dependencies can be written out in the following form:

Where [gly] describes the concentration of glyphosate, μ the growth rate of the strains and S the saturation of optical density

This set of differential equations can't be solved elegantly in an analytical way, therefore we used numerical integration to model the optical density (representing the cell mass) of both of our strains in competition. This is quite conveniently done by employing the NDSolve method of the computing system mathematica.

Now we had all the necessary information to model the cell mass of our both strains in competition for different glyphosate concentrations over the course of several hours. The outcome of this simulation is best depicted as ratio of the optical densitiy of our wild type and the double mutant representing the ratio of cell mass between the two (Figure X):

Figure X. Ratio of ODWild type/ODmutant against glyphosate concentration and time as it can be modeled from our previous knowledge of the dependencies between glyphosate and growth rate.

HIER FEHLT NOCH EIN VERGLEICH MIT DEN EXPERIMENTELLEN DATEN, DAFÜR BENÖTIGEN WIR AUF JEDEN FALL EINE GRAFIK, DIE DIE GESAMT-OD BIS 2 mM DARSTELLT.

At low concentrations of glyphosate, the ratio between mutant and wild type is near 1, both strains grow at the same speed and have similar cell mass. But increasing the dose of glyphosate by only as little as 2 mM, the ratio changes in favor of the mutant strain. As it can be seen in the hill plot (Figure X), the mutant strain shows a much lower reduction in growth rate, therefore it should grow much faster than the wild type shifting the ratio in cell mass.
Overall our modeling showed, that based on our assumptions of glyphosate as the dominating factor in our experiments, the ratio between the wild type and the mutant strain should pretty drastically change in favor of the mutant by increasing the glyphosate concentration in the media.

Evaluating growth behaviour of the strains by means of fluorescence

As the OD of both strains can't be measured separated from each other, no comparison between the strains by means of OD is possible. Therefore we created mutant strains expressing the fluorophores green fluorescent protein (GFP, wild type) and mOrange (double mutant) to monitor the fluorescence of both strains over the course of the experiment. We calculated the increase in fluorescence relative to the fluorescence measured at t=0 for different glyphosate concentrations. Doing this we were able to evaluate if our proposed trend of the mutant strain being the stronger competitor at higher glyphosate concentrations is right. Although it has to be mentioned that this test is more of qualitative nature, as modelling the fluorescence signals of the strains needs further research:
To properly model the fluorescence trends of the strains a link between optical density and fluorescence has to be found. To achieve this expression, maturation and degradation rates of the proteins are needed as the model of Leveau and Lindow suggests (1).
The calculated relative fluorescence and a ratio between the strands for 0 mM and 2 mM glyphosate is shown in the following Figure (Figure X):

Figure X: N-fold increases in fluorescence relative to fluorescence at t=0 (OD = 0.1) for both strains in competition. Upper left: n-fold increase in fluorescence of the wild type strain. Upper right: n-fold increase in fluorescence of the double mutant. Bottom: Ratio of relative fluorescence increases.

Without the influence of glyphosate the fluorescence of both strains increased nearly equally (Figure X, upper left and right), at the 10 h mark a ratio of 0.8 between the strains can be found.
Adding 2 mM glyphosate to the system strongly shifts the ratio of fluorescence. The double mutant is strongly expressing fluorescent proteins, nearly doubling the amount of relative fluorescence increase compared to the competition experiment without glyphosate (Figure X, upper right). In contrast the wild type strain shows no significant increase in fluorescence over the course of the experiment, at maximum a 2-fold increase can be observed (Figure X, upper left). The ratio of fluorescence decreases to less than 0.1 after ten hours of competition, suggesting a much stronger growing double mutant.

Even though these results can't be interpreted quantitatively they still show, that our assumptions might be right. While without glyphosate both strains show a rather equal increase in fluorescence intensity, adding glyphosate shifts the ratio strongly in favour of the double mutant strain. As higher fluorescence intensity correlates with a stronger expression of fluorescent proteins, it can be assumed, that the double mutant clearly dominates in growth rate in comparison to the wild type, as our modelled already suggested.

Evaluating cell growth with a β-galactosidase assay

Another method for distinction of two strains in one culture was an enzymatic solution. The sensitive strain expresses constitutively a β-galactosidase which leads to hydrolyzation of a chromogenic substrate. This assay is performed with the substrate o-nitrophenyll-β-D-galactopyranoside (ONPG). The hydrolyzation of ONPG leads to the formation of o-nitrophenol, which absorbs light at λ=420 nm. The amount of o-nitrophenol that is formed within a certain time is directly proportional to the activity of a promoter. Since the promoter is active in any cell from the sensitive strain, it is proportional to the growth. Measurements of the enzyme activity can be seen in Figure X:

Figure X. Cell mixture (1:1) of strains BP193 and BP235 (ΔgltT and ΔgltP) were incubated with increasing glyphosate concentrations (0 - 5 mM). Strain BP193 harbors the lacZ gene behind an artificial promoter, but does not tolerate glyphosate in contrast to strain BP235 (ΔgltT and ΔgltP). ß-galactosidase activity can be measured at a wavelength of 420 nm using the substrate o-nitrophenyl-ß-D-galactopyranoside. Activity decreases with increasing glyphosate concentration. BP235 does not contain the lacZ gene, hence no activity can be observed.

High concentrations of glyphosate lead to a strongly reduced enzyme activity, suggesting slower or stagnated growth of the wild type strain. These results complement what our model of the cell mass ratio suggested: For glyphosate concentrations higher than 2 mM the wild type strain is strongly inhibited in growth, reassuring that glyphosate is the dominant factor in our competition and single experiments.

This assay can also be performed with 5-bromo-4-chloro-3-indolyl-β-D-galactopyranoside (X-gal) as chromogenic substrate. The hydrolyzation leads to the formation of an insoluble indigo-dye that is located within the cells. With increasing glyphosate concentration, the intensity of the blue color decreases and is therefore an indicator for the glyphosate concentration. This method is further prove that the sensitive wild type strain is strongly reduced in growth and therefore eliminated by the glyphosate-resistant double mutant.