Computational model, method, and system for kinetically-tailoring multi-drug chemotherapy for individuals

• US 7,286,970 B2
• Filed: 01/07/2003
• Issued: 10/23/2007
• Est. Priority Date: 01/07/2002
• Status: Expired due to Fees

• Current Assignee: Lawrence Livermore National Security LLC (Government of the United States of America)
• Original Assignee: Regents of the University of California (University of California)
• Inventors: Gardner, Shea Nicole
• Primary Examiner(s): Moran; Marjorie A.
• Application Number:
  US10/338,559
• Publication Number:
  US 20030144798A1
• Time in Patent Office: 1,750 Days
• Field of Search: None
• US Class Current: 703/2
• CPC Class Codes: G16H 20/10 relating to drugs or medications, e.g. for ensuring correct administration to patients
• CPC Class Codes: G16H 20/40 relating to mechanical, radiation or invasive therapies, e.g. surgery, laser therapy, dialysis or acupuncture
• CPC Class Codes: G16H 50/50 for simulation or modelling of medical disorders
• CPC Class Codes: G16Z 99/00 Subject matter not provided for in other main groups of this subclass

Abstract

A method and system for tailoring treatment regimens to individual patients with diseased cells exhibiting evolution of resistance to such treatments. A mathematical model is provided which models rates of population change of proliferating and quiescent diseased cells using cell kinetics and evolution of resistance of the diseased cells, and pharmacokinetic and pharmacodynamic models. Cell kinetic parameters are obtained from an individual patient and applied to the mathematical model to solve for a plurality of treatment regimens, each having a quantitative efficacy value associated therewith. A treatment regimen may then be selected from the plurlaity of treatment options based on the efficacy value.

Claims

1. A method for tailoring multidrug chemotherapy treatment regimens to individual patients with diseased cells exhibiting evolution of resistance to drugs used in such treatments, said method comprising the steps of:

providing a set of differential equations representing a mathematical model of rates of population change of proliferating and quiescent diseased cells, said mathematical model including parameters corresponding to cell kinetics and evolution of drug resistance of the diseased cells to cell-cycle phase-specific cytotoxic drugs, cell-cycle phase non-specific cytotoxic drugs, and cytostatic drugs;

obtaining cell kinetic parameter values from an individual patient including proliferative index, apoptotic index, S-phase fraction, cell cycle time, and level of drug resistance;

using the cell kinetic parameter values obtained from the patient to solve the set of differential equations of the mathematical model to determine a plurality of treatment regimens for the patient, each of said treatment regimens having a quantitative efficacy value associated therewith representing the efficacy thereof in reducing the diseased cell population of the patient; and

selecting one of the treatment regimens having an efficacy value that is desirable for treating the patient, said selected treatment regimen thereby being tailored to the patient.

2. The method of claim 1,

wherein the set of differential equations representing the mathematical model further incorporates a Gompertzian curve for modeling tumor growth.

3. The method of claim 1,

wherein the set of differential equations representing the mathematical model additionally includes parameters corresponding to evolution of cell kinetics.

4. The method of claim 1,

wherein the set of differential equations representing the mathematical model additionally includes parameters corresponding to a transition rate from proliferating to quiescent compartments, and a transition rate from quiescent to proliferating compartments, with the proliferative and apoptotic indices obtained from the patient used to determine transition rate from proliferating to quiescent compartments.

5. The method of claim 1,

wherein the set of differential equations representing the mathematical model additionally includes a parameter corresponding to a time-dependent drug concentration.

6. The method of claim 1,

wherein the set of differential equations representing the mathematical model additionally includes a parameter corresponding to relative average dose intensity, for modeling toxicity.

7. The method of claim 1,

wherein the efficacy value is the probability of survival.

8. A computer program product comprising:

a computer-readable medium having a computer-readable set of instructions embodied in said medium for enabling a computer system to provide multi-drug chemotherapy cancer treatment alternatives to individual patients based on tumor cell kinetics of the patient, the set of instructions comprising;

a system model code for causing the computer system to provide a set of differential equations representing a mathematical model of rates of population change of proliferating and quiescent diseased cells, said mathematical model including parameters corresponding to cell kinetics and evolution of drug resistance of the diseased cells to cell-cycle phase-specific cytotoxic drugs, cell-cycle phase non-specific cytotoxic drugs, and cytostatic drugs;

an input code for causing the computer system to receive tumor cell kinetic parameters obtained from the individual patient, the cell kinetic data including proliferative index, apoptotic index, S-phase fraction, cell cycle time, and level of drug resistance;

a processing code for causing the computer system to solve the set of differential equations of the mathematical model using the tumor cell kinetic parameters obtained from the patient to produce a plurality of alternative treatment regimens for the patient, each of said treatment regimens having a quantitative efficacy value associated therewith representing the efficacy thereof in reducing the tumor cell population of the patient; and

an output code for causing the computer system to present the plurality of treatment regimens so that one of the treatment regimens may be selected based on an efficacy value that is desirable for the patient.

9. A system for providing multi-drug chemotherapy cancer treatment alternatives to individual patients based on tumor cell kinetics of the patient, said system comprising:

an input module for receiving cell kinetic data obtained from an individual patient, the cell kinetic data including proliferative index, apoptotic index, S-phase fraction, cell cycle time, and level of drug resistance;

a system processor for using the cell kinetic data obtained from the patient in the mathematical model to solve a set of differential equations representing a mathematical model of rates of population change of proliferating and quiescent tumor cells, said mathematical model including parameters corresponding to cell kinetics and evolution of drug resistance of the tumor cells to cell-cycle phase-specific cytotoxic drugs, cell-cycle phase non-specific cytotoxic drugs, and cytostatic drugs, and to determine for a plurality of treatment regimens for the patient, each of said treatment regimens having a quantitative efficacy value associated therewith representing the efficacy thereof in reducing the tumor cell population of the patient; and

an output module for presenting the plurality of treatment regimens, so that one of the treatment regimens may be selected based on an efficacy value that is desirable for the patient.

10. The system of claim 9,

wherein the set of differential equations representing the mathematical model additionally includes parameters corresponding to evolution of cell kinetics.

11. The system of claim 9,

wherein the set of differential equations representing the mathematical model additionally includes parameters corresponding to a transition rate from proliferating to quiescent compartments, and a transition rate from quiescent to proliferating compartments, with the proliferative and apoptotic indices obtained from the patient used to determine transition rate from proliferating to quiescent compartments.

12. The system of claim 9,

wherein the set of differential equations representing the mathematical model additionally includes a parameter corresponding to a time-dependent drug concentration.

13. The system of claim 9,

wherein the set of differential equations representing the mathematical model additionally includes a parameter corresponding to relative average dose intensity, for modeling toxicity.

14. The system of claim 9,

wherein the efficacy value is the probability of survival.