Drug Combination in vitro Experiment Simulation

Lingge Li

November 2016

This web app is built for in vitro drug combination experiment simulation. In this setting, for a single drug, the dose-repsonse curve is usually modelled with a sigmoid function, known as the Emax equation. It depends on two parameters, IC50 and m (shape) where IC50 is often the target of interest.

Before simulating an experiment, we need to come up with a hypothetical drug. E0 is the number of cells we expect to see in a control well without any drug. IC50 is the median effect concentration (an implicit assumption is that the drug will kill all the cells at a very high concentraion). m the shape parameter usually ranges from 0 to 2 with higher value making the curve more sigmoidal. The default plot is what we expect to observe as raw data.

Here we focus on the simplest experiment where only one cell line is used on one plate. The experimental design involves choosing a range of drug concentrations. For example, on a 24-well plate, it is common to have 4 controls and 5 concentrations with 4 replicates. In practice, drug concentrations often follow 2-fold serial dilution. Note that concentration 0 corresponds to control.

The noise is generated with $$[E(Y)]^{\lambda}\cdot\epsilon\text{ where }\epsilon\sim N(0,\sigma).$$ Obviously, sigma is the standard deviation of Gaussian noise. When lambda is 0, the noise would be independent. When lambda is 1, the noise would be proportional to the expected cell count. In other words, there would be more variation where the number of cells is large.

The most popular data analysis method transforms the data to log scale and fits a straight line. The output shown below includes the estimated parameter values and confidence intervals.

Model A

              
Model B

              
Model M

              
L20

              
L50

              
L80