Team:ETH Zurich/Modeling/Microfluidics

Overview

In our inital model, we chose a 5 cm long cylindrical channel with a radius of 1 mm (or equivalently 2 mm diameter). This channel is connected to a reservoir filled with medium and acetaldehyde. For the purpose of modeling the diffusion and degradation of acetaldehyde in the channel, we assumed that the reservoir has constant acetaldehyde concentration (or equivalently has infinite size).

Diffusion

First of all, we have to model diffusion of Acetaldehyde from the reservoir to the closed end of the microfluidics channel. The relevant partial differential equation is the one for isotropic diffusion, i.e. diffusion that is uniform in all directions:

For the relevant diffusion parameter, D_AcAl, see the microfludics model section on the parameters page.

Degradation

Cell Density

In order to able to determine how much Acetaldehyde our smoColi cells are able to degrade, we first need to determine the cell density of the immobilized cells in agarose. For this calculation, we assume we have a cell culture with OD₆₀₀ with which we fill the channel. Parameter-wise, to determine the relative cell density in the channel voume, we need to know:

• the mean cell density in [cells/volume] at OD₆₀₀=1, cd_{OD600, undiluted}, and
• the mean cell volume of a single cell, V_eColi.

A parameter we can choose suitably is the dilution of the cells in agarose, c_Dilution. Then we can calculate the relative cell density in the smoColi-agarose-filled microfluidics channel:

Enzymatic Degradation

Now that we know what the cell concentration in the channel is, we can model the degradation of the toxic molecule by the cell culture in agarose. This process, together with diffusion of more toxic molecules from the reservoir, will finally generate the toxic molecule gradient.

In our enzymatic degradation model, we assume diffusion into the cells is fast compared to the degradation process. The degradation process then is modeled using Michaelis-Menten enzyme reaction kinetics:

The constants are the usual ones for Michaelis-Menten-type reactions:

• v_max,AcAl is the maximum reaction rate, i.e. the rate the toxic molecule degradation reaction occurs with at saturation point, and
• K_M,AcAl is the Michaelis Menten constant, i.e. the Acetaldehyde concentration at which the reaction rate is half of v_max,AcAl

Estimation of v_max,AcAl

--- Degradation of toxic molecule ---

Simulation

To see if we can indeed create an acetaldehyde gradient, we first solved for the gradient steady state in the channel. We chose an acetaldehyde concentration of 100 mg/l for the reservior in these simulations.

Steady State

We can see that the diameter does not influence the mean dynamics of the acetaldehyde gradient. We note that there may be stochastic behavior for very thin channels, since there will be only very few cells. For the channel diameter we chose a mean field approximation should be sufficient though as the amount of cells involved in the degradation is sufficiently high.

Dynamics

We were also interested in finding out the time scale which is needed to create the gradient. Therefore we also analyzed the dynamics in a spatiotemporal reaction (degradation) - diffusion simulation.