Computer simulation of leukemia therapy: combined pharmacokinetics, intracellular enzyme kinetics, and cell kinetics of the treatment of L1210 leukemia by cytosine arabinoside.
An integrated mathematic computer-based model of the pharmacokinetics, intracellular enzyme kinetics, and cell kinetics of the treatment of L1210 leukemia by cytosine arabinoside (ara-C) is described. The compartment model of Bischoff and Dedrick is extended to the intracellular level by inclusion of equations describing the phosphorylation, dephosphorylation, and deamination of ara-C with enzymatic feedback control. The activities of kinase, deaminase, and phosphatase are explicitly included in the models and are estimated from relevant data. Cell proliferation is described by a continuous-flow mathematic model in which cellular maturation and cell-to-cell variability in maturation rates are key variables. Cell proliferation is related to intracellular biochemistry through mathematic expressions which relate cell lethality and progression delay to the time course of intracellular ara-CTP. In vitro and in vivo experiments performed in a number of laboratories are compared by simulation. The most sensitive parameters in dose-response and cell-survival simulations are deoxycytidine kinase activity, ara-CTP half-life, renal clearance of ara-C, and cell-kinetic parameters for proliferation and cell killing. Progression delay is vital to the realistic simulation of divided-dose schedules. By comparative simulation we have identified areas of uncertainty which can be classified by a few additional measurements. The applications of simulations combining pharmacokinetic, biochemical, and cell-kinetic data in vitro and in vivo are discussed, exploring consistency among different measurements, and relating experimental protocols to clinical treatment.