A novel mathematical model describing adaptive cellular drug metabolism and toxicity in the chemoimmune system

Cells cope with the threat of xenobiotic stress by activating a complex molecular network that recognizes and eliminates chemically diverse toxic compounds. This "chemoimmune system" consists of cellular Phase I and Phase II metabolic enzymes, Phase 0 and Phase III ATP Binding Cassette (ABC) membrane transporters, and nuclear receptors regulating these components. In order to provide a systems biology characterization of the chemoimmune network, we designed a reaction kinetic model based on differential equations describing Phase 0-III participants and regulatory elements, and characterized cellular fitness to evaluate toxicity. In spite of the simplifications, the model recapitulates changes associated with acquired drug resistance and allows toxicity predictions under variable protein expression and xenobiotic exposure conditions. Our simulations suggest that multidrug ABC transporters at Phase 0 significantly facilitate the defense function of successive network members by lowering intracellular drug concentrations. The model was extended with a novel toxicity framework which opened the possibility of performing in silico cytotoxicity assays. The alterations of the in silico cytotoxicity curves show good agreement with in vitro cell killing experiments. The behavior of the simplified kinetic model suggests that it can serve as a basis for more complex models to efficiently predict xenobiotic and drug metabolism for human medical applications.