Abstract
In this manuscript, we use an exactly solvable stochastic binary model for the regulation of gene expression to analyze the dynamics of response to a treatment aiming to modulate the number of transcripts of a master regulatory switching gene. The challenge is to combine multiple processes with different time scales to control the treatment response by a switching gene in an unavoidable noisy environment. To establish biologically relevant timescales for the parameters of the model, we select the RKIP gene and two non-specific drugs already known for changing RKIP levels in cancer cells. We demonstrate the usefulness of our method simulating three treatment scenarios aiming to reestablish RKIP gene expression dynamics toward a pre-cancerous state: (1) to increase the promoter's ON state duration; (2) to increase the mRNAs' synthesis rate; and (3) to increase both rates. We show that the pre-treatment kinetic rates of ON and OFF promoter switching speeds and mRNA synthesis and degradation will affect the heterogeneity and time for treatment response. Hence, we present a strategy for reaching increased average mRNA levels with diminished heterogeneity while reducing drug dosage by simultaneously targeting multiple kinetic rates that effectively represent the chemical processes underlying the regulation of gene expression. The decrease in heterogeneity of treatment response by a target gene helps to lower the chances of emergence of resistance. Our approach may be useful for inferring kinetic constants related to the expression of antimetastatic genes or oncogenes and for the design of multi-drug therapeutic strategies targeting the processes underpinning the expression of master regulatory genes.