Team:Imperial College/Motility Control

We modelled the motility of motile B.subtilis using a simple mechanical model as shown above. Using experimental data from cell tracking, we were able to extract parameters which describe the model. The model and its parameters are defined here. In summary, parameter A is the ratio of the flagellum force to the viscosity of the medium, B is the initial velocity of the cell and alpha is the ratio of the vicosity of the medium to the cell's mass.

The following figure shows the results of our model fitting. We have introduced a change in flagellar force at certain points of the cell trajectory so as to achieve a better fit. A maximum of two runs were allowed for each cell trajectory.

The MATLAB Distribution Fitting Tool was used to model the distribution of parameter A, which is a representation of the flagellum force, assuming that the medium is homogenous and its viscosity is constant throughout the medium. Parameter A was the found to be exponentially distributed. The following figures describe the probability density function (PDF) and cumulative distribution function (CDF).

We fitted the data according to the simple mechanical model we have developed. From our model fitting process, we can see that flagellum force is exponentially distributed. Our mechanical model though simple, fits the cell trajectory data extremely well as shown in the figure above. Further work which can be done would be to utilise a movable stage to track the movement of B. subtilis over its entire run so as to obtain a distribution of other motility parameters associated running and tumbling events.