The role of stochastic gene switching in determining the pharmacodynamics of certain drugs: basic mechanisms

In this paper we analyze the impact of the stochastic fluctuation of genes between their ON and OFF states on the pharmacodynamics of a potentially large class of drugs. We focus on basic mechanisms underlying the onset of in vitro experimental dose-response curves, by investigating two elementary molecular circuits. Both circuits consist in the transcription of a gene and in the successive translation into the corresponding protein. Whereas in the first the activation/deactivation rates of the single gene copy are constant, in the second the protein, now a transcription factor, amplifies the deactivation rate, so introducing a negative feedback. The drug is assumed to enhance the elimination of the protein, and in both cases the success of therapy is assured by keeping the level of the given protein under a threshold for a fixed time. Our numerical simulations suggests that the gene switching plays a primary role in determining the sigmoidal shape of dose-response curves. Moreover, the simulations show interesting phenomena related to the magnitude of the average gene switching time and to the drug concentration. In particular, for slow gene switching a significant fraction of cells can respond also in the absence of drug or with drug concentrations insufficient for the response in a deterministic setting. For higher drug concentrations, the non-responding fraction exhibits a maximum at intermediate values of the gene switching rates. For fast gene switching, instead, the stochastic prediction follows the prediction of the deterministic approximation, with all the cells responding or non-responding according to the drug dose.