A dynamical systems model to simulate the perturbation kinetics of gene expression by antisense oligonucleotides

Antisense oligonucleotides owe their efficacy to an ability to induce RNase H-dependent suppression of RNA translation, for sufficient time to allow physiological proteolysis. The magnitude and time delay preceding the protein nadir concentration determine the extent and timing of maximum antisense oligonucleotide activity. Antisense oligonucleotide degradation underlies reversal of RNA downregulation. The kinetics of protein downregulation is therefore determined by the complex interaction of both ligand chemistry (nuclease stability, affinity and RNase H activation), and gene expression kinetics. Optimization of antisense oligonucleotide efficacy and experimental design requires understanding of these interactions. The kinetics of protein and RNA downregulation have therefore been simulated by analysing a two-compartment kinetic model incorporating RNase H-dependent transcript degradation. The system of nonlinear differential equations describing this model was solved numerically using Runge-Kutte integration. The timecourse solutions corresponding to the four state variables (RNA, protein, antisense/RNA heteroduplex and antisense oligonucleotide), were determined simultaneously. This allowed systematic in silico examination of the consequences of altering variables such as oligonucleotide concentration, affinity, and stability, or the scheduling of multiple transfections on RNA and protein perturbations. By providing a tool for examining antisense oligonucleotide action theoretically, this heuristic model should facilitate both the rational design and interpretation of antisense experiments.