Drug interaction prediction using ontology-driven hypothetical assertion framework for pathway generation followed by numerical simulation

1471-2105-9-S6-S11 1471-2105 Research Drug interaction prediction using ontology-driven hypothetical assertion framework for pathway generation followed by numerical simulation Arikuma Takeshi atake@bio.cs.titech.ac.jp Yoshikawa Sumi sumi@bio.cs.titech.ac.jp Azuma Ryuzo azuma@bio.cs.titech.ac.jp Watanabe Kentaro kentaro.watanabe.dc@hitachi.com Matsumura Kazumi matsumura-k@admet.daiichichem.co.jp Konagaya Akihiko kona@bio.cs.titech.ac.jp

Department of Computer Science, Tokyo Institute of Technology, 2-12-1 Oookayama, Meguro, Tokyo, Japan

Genomic Sciences Center, RIKEN, 1-7-22 Suehiro, Tsurumi, Yokohama, Kanagawa, Japan

BMC Bioinformatics Symposium of Computations in Bioinformatics and Bioscience (SCBB07) Guoqing Lu, Jun Ni, Thomas L Casavant and Brian Athey Research Symposium of Computations in Bioinformatics and Bioscience (SCBB07) Iowa City, Iowa, USA

13–15 August 2007

http://www.imsccs-conference.org/imsccs07/SCBB07/index.html 1471-2105 2008 9 Suppl 6 S11 http://www.biomedcentral.com/1471-2105/9/S6/S11 18541046 10.1186/1471-2105-9-S6-S11 28 5 2008 2008 Arikuma et al; licensee BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

In accordance with the increasing amount of information concerning individual differences in drug response and molecular interaction, the role of in silico prediction of drug interaction on the pathway level is becoming more and more important. However, in view of the interferences for the identification of new drug interactions, most conventional information models of a biological pathway would have limitations. As a reflection of real world biological events triggered by a stimulus, it is important to facilitate the incorporation of known molecular events for inferring (unknown) possible pathways and hypothetic drug interactions. Here, we propose a new Ontology-Driven Hypothetic Assertion (OHA) framework including pathway generation, drug interaction detection, simulation model generation, numerical simulation, and hypothetic assertion. Potential drug interactions are detected from drug metabolic pathways dynamically generated by molecular events triggered after the administration of certain drugs. Numerical simulation enables to estimate the degree of side effects caused by the predicted drug interactions. New hypothetic assertions of the potential drug interactions and simulation are deduced from the Drug Interaction Ontology (DIO) written in Web Ontology Language (OWL).

Results

The concept of the Ontology-Driven Hypothetic Assertion (OHA) framework was demonstrated with known interactions between irinotecan (CPT-11) and ketoconazole. Four drug interactions that involved cytochrome p450 (CYP3A4) and albumin as potential drug interaction proteins were automatically detected from Drug Interaction Ontology (DIO). The effect of the two interactions involving CYP3A4 were quantitatively evaluated with numerical simulation. The co-administration of ketoconazole may increase AUC and Cmax of SN-38(active metabolite of irinotecan) to 108% and 105%, respectively. We also estimates the potential effects of genetic variations: the AUC and Cmax of SN-38 may increase to 208% and 165% respectively with the genetic variation UGT1A1*28/*28 which reduces the expression of UGT1A1 down to 30%.

Conclusion

These results demonstrate that the Ontology-Driven Hypothetic Assertion framework is a promising approach for in silico prediction of drug interactions. The following future researches for the in silico prediction of individual differences in the response to the drug and drug interactions after the administration of multiple drugs: expansion of the Drug Interaction Ontology for other drugs, and incorporation of virtual population model for genetic variation analysis, as well as refinement of the pathway generation rules, the drug interaction detection rules, and the numerical simulation models.

Background

The role of in silico prediction of drug interactions on the pathway level is becoming more and more important for solving real-world problems. Multiple-drug regimens exemplify the need for the computer-assisted prediction of drug interactions. Multiple-drug regimens are commonly prescribed for elderly patients suffering from more than one disease. However, they sometimes cause unexpected severe side effects because of the drug interactions or individual differences concerning response to the drugs 1. Therefore, the prediction of drug interactions for preventing the side effects is an important issue for these regimens.

On the other hand, information useful for in silico drug interaction prediction has increased very rapidly in recent years. Technological innovations in genomic sciences have produced an enormous amount of biomolecular information including sequences, structures, and pathways. In order to integrate the biomolecular information, ontology is attracting a lot of attention 23. In addition, pharmacokinetics modeling and simulation are emerging, promising techniques to understand the dynamic behavior of drug metabolic pathways 45.

To develop a practical in silico drug interaction prediction system by integrating the above information and techniques, the following issues must be solved.

Context dependency of drug-metabolic pathways

Drug-metabolic pathways do not exist a priori. They strongly depend on contexts and situations including the administration route, single nucleotide polymorphism (SNP) of drug-response genes, and the administration of multiple drugs and foods. Therefore, a dynamic reconstruction of drug metabolic pathways from primitive molecular events is necessary for drug interaction prediction on the pathway level. Such reconstruction requires the formal definition of molecular events and the relations among them, i.e., the Drug Interaction Ontology (DIO).

Treatment of multiscale events

Pathways triggered by drug administration consist of multiscale events: from the molecular level to the body level, ranging from nanoseconds to hours or days in terms of drug response. For example, drug administration and drug excretion are body-level events, while drug transport and enzymatic reactions are molecular-level events. A comprehensive view from the molecular level to the body level is necessary in order to understand multi-scale events.

Quantitative evaluation of interactions

Quantitative evaluation plays an essential role to estimate the degree of side effects caused by drug interactions. More than one drug interactions may occur in drug-metabolic pathways from the qualitative reasoning view point. However, not all drug interactions cause side effects because of differences in binding affinity and molecular population. Quantitative simulation models with an in silico drug interaction prediction system must be incorporated to discriminate serious drug interactions from negligible ones. It should be also noted that total drug metabolism depends on not only kinetic parameters but also physiological parameters such as organ volumes and blood flows. Incorporation of kinetic parameters and physiological parameters is necessary for a realistic simulation model to recapture experimental data.

In this paper, we propose a new drug-interaction prediction framework called "Ontology-Driven Hypothetic Assertion (OHA)" focusing mainly on automatic generation of drug-metabolic pathways, automatic detection of drug interactions on multiple-drug regimens 6, and quantitative evaluation of the drug interactions with numerical simulation. The effectiveness of the OHA framework was demonstrated in the prediction of interactions between irinotecan (CPT-11) and ketoconazole.

Design philosophy

Ontology-Driven Hypothetical Assertion

The Ontology-Driven Hypothetical Assertion framework maps inferences onto the Drug Interaction Ontology as both machine and human understandable form. The inferences include the results of qualitative reasoning and numerical simulation, i.e. generated pathways and detected drug interactions and generated differential equations and simulation results. The assertions enable to interpret the obtained results such as drug interaction candidates, simulation models, the pharmacokinetic moment values with background knowledge on pharmaceutical science and biochemistry in the ontology.

Causality-based pathway modularization

Modularization is necessary for the dynamic reconstruction of pathways that depend on dose conditions. Careful selection of primitive modules is the key to ensuring the soundness of pathway reconstruction. Molecular events, such as molecular transport and enzymatic reactions, are well-formed primitive modules for this purpose. In this paper, we refer to the primitive modules as "molecular events", and the aggregation of molecular events as "pathways".

To avoid redundant pathway branch constructions, which are non-essential for the target drug interactions, we adopt causality-based modularization in which each molecular event is defined by the unique relationship between key molecules before and after the event. The triadic relationship <trigger, situator, resultant>, proposed by Yoshikawa et al. 7, is one such causality that can be commonly found in molecular reactions. For example, in case of enzymatic reactions, substrates, enzymes and other products correspond to trigger, situator and resultant, respectively. In the case of molecular transport, extra (intra) cellular molecules, transporters and intra (extra) cellular molecules correspond to the participants of the triadic relation, respectively. The triadic relationship can be applicable to higher level events like drug dosage and drug excretion, as long as its causality is unique and clear.

Figure 1 (a)(b)(c) shows a simple example of pathway reconstruction with two primitive molecular events: an enzymatic reaction in which carboxylesterase (CE) metabolizes irinotecan into SN-38 (7-ethyl-10-hydroxycamptothecin) in the liver, and molecular transport in which SN-38 in the liver is transported to the bile by MRP2 (Multidrug resistance-associated protein 2). Two molecular events are connected at the resultant of the enzymatic reaction (TR0000019) and the trigger of the molecular transport (TR0000007) for passing SN-38 in the liver (SN-38@liver) to the bile.

Figure 1

Decomposition of a pathway with primitive molecular events

Decomposition of a pathway with primitive molecular events. a) Ordinary pathway map representation. b) Primitive molecular events including trigger (T), situator (S) and resultant (R). c) Aggregation of molecular events at SN-38@liver. Rectangles and ellipses are instances of continuants and processes, respectively, in b) and c).

Ontology-based knowledge base construction

The ontological approach in knowledge base design is adopted for resource sharing and the semantic description of molecular events and pathways. Ontology is necessary to define molecular events and pathways in a form that can be shared among computers and human beings. This enables the full use of powerful computational intelligence for dynamic pathway reconstruction in a way that human intelligence can follow and understand. Ontology is also important for establishing interoperability among web resources and thereby to make use of the latest drug reaction information published in the semantic web 28. Public biological ontologies, especially in the field of chemical biology, are now dramatically increasing, and have a great potential to develop sustainable knowledge bases for molecular reaction and pathways.

One of the unique features of our knowledge-base design, i.e. DIO, is the adoption of the triplet view: intension, extension, attribute for molecular events and pathways as seen in Figure 2. Intension defines the kinds of molecular events and pathways as controlled vocabulary of processes. Attribute defines the components of molecular events and pathways as controlled vocabulary of continuants (see Figure 3 for a part of class hierarchy of controlled vocabulary for drug interaction). These controlled vocabularies are constructed from the viewpoints of "process" and "continuant" as proposed in Basic Formal Ontology (BFO) 910. Molecular event objects are asserted in extension using the controlled vocabularies defined in intention and attributes. As for extension, prototype object modeling is adopted to represent pathway as aggregation of molecular events. This is because an infinite number of terms or classes are required to express all combinations of molecular events. We avoid this problem by treating pathways as dynamic objects deduced from prototype molecular event objects rather than treating them as instances of pathway classes.

Figure 2

Overview of the Ontology-Driven Hypothetic Assertion (OHA) framework

Overview of the Ontology-Driven Hypothetic Assertion (OHA) framework. Drug Interaction Ontology (DIO) provides triple views for events and pathways, namely, intension, extension, and attribute. White triangles in extension are molecular events such as CYP3A4 mediated metabolism of irinotecan. White rectangles and circles in attribute are molecules and organs respectively; the organs indicate the location of molecules. The gray triangles and arrowed dashed lines represent the hypothetic assertions generated by the OHA framework. The Pathway Object Constructor (POC) generates pathways by connecting molecular events. The Drug Interaction Detector (DID) detects drug interactions from the generated pathways and asserts the interactions as hypothetic assertions into DIO.

Figure 3

Class hierarchy for controlled vocabularies of DIO

Class hierarchy for controlled vocabularies of DIO. Classes prefixed with "fma:" were references for the classes defined in the FMA. The class hierarchy was stored in "dio_event.owl". According to BFO, the controlled vocabulary was divided into two parts; process (e.g. events) and continuants (e.g. things). (This figure was generated using ontoviz plug-in of Protege with manual changes.).

Dynamic pathway reconstruction and drug interaction detection

The drug metabolic pathway, due to its dynamic nature, is difficult to define a priori in the manner seen in biomolecular metabolic pathways in the Kyoto Encyclopedia of Genes and Genomics (KEGG) 11. Therefore, the OHA framework provides the Pathway Object Constructor (POC) and Drug Interaction Detector (DID). Pathway Object Constructor dynamically deduces pathways from the DIO, depending on contexts or situations. Drug Interaction Detector detects drug interaction candidates finding intersections from those generated pathways.

These inferences are mapped onto the Drug Interaction Ontology as hypothetical assertions. The generated pathways and detected drug interactions are asserted as aggregation of molecular events. The detected drug interactions can be segmentalized such as competitive inhibition, noncompetitive inhibition, and uncompetitive inhibition when binding site information is available. This is effective to select an appropriate differential equation depending on the inhibition model.

Numerical simulation

In order to incorporate a quantitative simulation system into the OHA framework, the following two aspects were considered: a methodology for automatic conversion from a generated pathway to a quantitative simulation model, and a methodology for the validation of numerical simulation models by analysing the dependences of kinetic parameters and physiological parameters. By solving these two issues, the OHA framework can can be applied to in silico prediction of individual drug interactions for multiple-drug regimens, assuming that kinetic parameters and the initial enzyme concentration are roughly estimated by individual genetic variations and health indices of bio-markers.

A simulation model is automatically generated from a pathway generated by the Pathway Object Constructor and Drug Interaction Detector. The pathway is converted to intermediate model by merging organs and molecular events, respectively, to fit a given simulation model framework such as compartment model. Differential equations for the simulation model are generated from the intermediate model by converting the merged events to mathematical expression.

We used the parameter-parameter dependency analysis system (PPD Viewer) designed by Konagaya et al. 12, with high-throughput numerical simulation engine and interactive visualization tools developed on OBIGrid 313141516. The system can predict the concentration/time profiles and those moment parameters such as area under the curve (AUC), area under the moment curve (AUMC), and mean resident time (MRT) when changing some kinetic parameters in the range of one thousandth to thousands of physiological conditions. The system visualizes dependencies among kinetic parameters and molecular concentrations in terms of moment parameters of a compound in target organs.

The generated simulation models and pharmacokinetic moment parameter values are mapped onto the Drug Interaction Ontology as hypothetical assertions. The simulation models are asserted as aggregations of objects representing terms and parameters in differential equations. Those objects have references to the components of the pathway objects from which the simulation models were generated. The moment parameter values are asserted with the drug interaction objects and the corresponding simulation model for the further analysis.

Interaction between Irinotecan and Ketoconazole

Irinotecan (CPT-11) is an anti cancer drug which is commonly used for colon and breast cancers 17. Irinotecan is a prodrug of SN-38, anti neoplastic topoisomerase I inhibitor, and is bioactivatied by carboxyl esterase (CE) 17. About 60% of irinotecan is excreated as unchanged drug from bile and kidney 18. Irinotecan is also metabolised by CYP3A4 to form APC and NPC 17. NPC is further metabolized by CE to form SN-38. SN-38 undergoes glucronidation by UGT1A1 to form the inactive glucronide, SN-38G 17. In addition, it is known that the mutations on UGT1A1, UGT1A1*28 which decreases the expression of UGT1A1 enzyme down to 30%, has strong relationship with some side effects of irinotecan 1920. Ketoconazole (KCZ) is an anti fungal drug and a well known inhibitor of CYP3A4. Ketoconazole undergoes extensive metabolism in the liver to form several metabolites 2122. About 2 to 4% of urinary radioactivity represents unchanged drug 22.

It has been reported that the inhibition of CYP3A4 by ketoconazole influences the metabolism of irinotecan, which results in 6% SN-38 increase 23.

Results and discussion

Results

Pathway generation and drug interaction detection

The pathways of intravenously administered irinotecan and orally administered ketoconazole were inferred as aggregation of molecular events by the Pathway Object Constructor. Figure 4 shows the generated irinotecan metabolic pathway object where irinotecan and its derivatives circulate through the veins, liver, bile, intestines, and portal vein, namely, the enterohepatic circulation, and are excreted through the kidneys or through the bile. These generated pathways are consistent with in vivo studies 17. The generated pathway object is sound in the sense that the object is deduced from the Drug Interaction Ontology (DIO) represented by Web Ontology Language (OWL-DL) 24.

Figure 4

Inferred pathway of irinotecan

Inferred pathway of irinotecan. Small rectangles and ellipses are instances of continuants and molecular events, respectively. Large rectangles such as small intestine are the compartments. Molecular event and trigger, situator, and resultant were connected by solid lines. Each molecular event was connected by the shared continuants. For example, TR0000020 and TR0000002 in the small intestine were connected because they share the same instance of SN-38. This connection can be interpreted as that irinotecan is metabolized into SN-38 (TR0000020) and then SN-38 is transported to the intestinal lumen. (This figure was generated by our inference system.).

Interactions between intravenously administered irinotecan and orally administered ketoconazole were detected and asserted by the Drug Interaction Detector. The detected drug interactions and the hypothetic assertion are shown in Figure 5. The assertion contains four drug interactions; two of them concern "drug binding reaction" to albumin in veins (ddi2) and arteries (ddi3), and the rest of them concern "oxidation" by CYP3A4 (ddi0 and ddi1). The detected drug interaction concerning CYP3A4 (ddi0 and ddi1) has been confirmed by the literature on in vivo studies 23.

Figure 5

Dynamically asserted instances of interactions

Dynamically asserted instances of interactions. Rectangles are owl classes corresponding to intension. Circles are owl instances corresponding to molecular events and attribute. Detected drug interactions were hypothetically asserted as owl instances (dotted circles) along with hasSubprocess properties (dotted and dashed edges).

Numerical simulation

We evaluated the effects of drug interactions concerning CYP3A4 quantitatively with numerical simulations. Three simulations were performed: sole administration of irinotecan for a patient having UGT1A1*1/*1 (wild type), co-administration of irinotecan and ketoconazole for a patient having UGT1A1*1/*1, and sole administration of irinotecan for a patient with UGT1A1*28/*28 which decreases the expression of UGT1A1 down to 30% than UGT1A1*1/*1. Intravenous drip infusion (125 mg/m², 90 min) was assumed for irinotecan, and oral administration (200 mg) was assumed for ketoconazole. Figure 6 shows the simulated concentration/time profiles of irinotecan, SN-38, APC, NPC, SN-38G in blood for the simulation of sole administration of irinotecan for a patient having UGT1A1*1/*1. The concentration/time profiles agree with experimental data reported by Slatter et al. 18.

Figure 6

Blood Concentration/time profile of irinotecan, SN-38, APC, NPC, and SN-38 in blood compartment for solo administration of irinotecan

Blood Concentration/time profile of irinotecan, SN-38, APC, NPC, and SN-38 in blood compartment for solo administration of irinotecan. The drug concentrations in blood without ketoconazole administration was estimated by numerical simulation. Intravenous drip infusion (125 mg/m², 90 min) was assumed for irinotecan administration.

Table 1 shows the comparisons of AUC and Cmax for the simulations. By the ketoconazole administration, the AUC of APC and NPC were decreased to 48.1% and 35.3%, respectively. The AUC of SN-38 were increased only to 108% by the ketoconazole administration. Similarly, the Cmax of APC and NPC were decreased to 25.6% and 20.2%, respectively, whereas the Cmax of SN-38 was increased to 105% by the ketoconazole administration. On the other hand, in case of the UGT1A1*28/*28 mutation, the AUC and Cmax of SN-38 were significantly increased: the AUC was increased to 208% and the Cmax was increased to 165%. This implies that patients with UGT1A1*28/*28 may suffer severe side effects when the doses are same as patients with UGT1A1*1/*1. These results agree with the previously published experimental papers by Kehrer et al. 23, Sai et al. 19, and Ando et al. 20. In addition, the simulation results were asserted into the ontology along with the simulation model.

Table 1

AUC and Cmax of blood concentration/time profile obtained by numerical simulations.

Irinotecan

APC

NPC

SN-38

SN-38G

AUC_KCZ-[nmol·min/ml]

672

233

14.2

29.8

137

Cmax_KCZ-[μM]

2.22

0.290

0.0113

0.0491

0.129

AUC_KCZ+[nmol·min/ml]

721

112

5.01

32.2

147

Cmax_KCZ+[μM]

2.29

0.0741

0.00228

0.0515

0.136

AUC_KCZ+/AUC_KCZ-[%]

107

48.1

35.3

108

107

Cmax_KCZ+/Cmax_KCZ-[%]

103

25.6

20.2

105

AUC_UGT[nmol·min/ml]

672

233

14.2

62.1

85.5

Cmax_UGT[μM]

2.22

0.290

0.0113

0.0812

0.0693

AUC_UGT/AUC_KCZ-[%]

100

208

62.4

Cmax_UGT/Cmax_KCZ-[%]

100

165

53.7

AUC and Cmax values for solo administration of irinotecan (AUC_KCZ-, Cmax_KCZ-), co-administration of irinotecan and ketoconazole (AUC_KCZ+, Cmax_KCZ+), and sole administration of irinotecan with patients having UGT1A1*28/*28 (AUC_UGT, Cmax_UGT) were obtained.

Discussion

Pathway reconstruction and drug interaction detection

The OHA framework automatically generated pathways for co-administration of irinotecan and ketoconazole, and detected drug interactions. The generated pathways and two of the drug interactions involving CYP3A4 were confirmed by published papers. The detected interactions with plasma albumin (ddi2 and ddi3) have not been reported in any literature as far as we know. One possible reason is that the effects of drug interactions may be negligible because of the high plasma albumin concentration compared to drug concentration. To distinguish serious drug interactions from those negligible ones, a qualitative evaluation, such as numerical simulation, is necessary.

Numerical simulation

The simulation model succeeded to reproduce Slatter's experimental data of irinotecan dosage to a large extent except for the Cmax of irinotecan blood concentration, cyclic fluctuation of irinotecan concentration and a bile cancer patient data.

The Cmax of drug concentrations in the simulation model strongly depend on renal clearance, bile clearance, and kinetic parameters of enzymes. Because these parameters were estimated from the recovery amount data in urine and feces by Slatter et al. 18 and in vitro experiments data reported as published papers 2526272829, the excess of irinotecan's Cmax in the simulation model might be explained by the false peak in the experimental data which results from time lag of sampling

Another possibility is the effect of reabsorption in intestine. The time course of the experimental data indicate slight increase of SN-38G blood concentration in the period of 12 hours. This may result from the reabsorption of SN-38 thorough small intestine. In order to include the reabsorption process into the simulation, the kinetic parameters in bacterial flora in the intestine and amount/time profile data of bile excretion are required.

Lastlym Slatter et al. reported that one patient with a bile cancer showed different bile and urine recovery pattern from other patients 18. Virtual population 30 may be helpful to reconstruct this behaviour.

Conclusion

Deduction of hypothetic drug interactions from the Drug Interaction Ontology was demonstrated for an irinotecan plus ketoconazole regimen. The prototype system detected four drug interactions. Two of them concerned cytochrome p450 (CYP3A4) and were consistent with known drug interactions. The other two concerned albumin, the effects of the interactions would be negligible as long as the drug concentration is low. We then quantitatively examined the effect of the drug interactions concerned CYP3A4 and the effect of genetic variation UGT1A1*28 using numerical simulation. The numerical simulation indicates that the drug interaction has only a limited effect on the pharmacokinetics of SN-38 for an irinotecan plus ketoconazole regimen. However, the genetic variation on UGT1A1 showed two-fold increase of SN-38's AUC.

Finally, in order to realize in silico prediction of drug interactions, the following future works are remained: expansion of the Drug Interaction Ontology for other drugs, incorporation of virtual population model for generic variation analysis, and refinements of the pathway generation rules, drug interaction detection rules, and the numerical simulation models.

Methods

The OHA framework

We implemented the following prototype system to prove the concept of the OHA framework. The system was designed to predict potential drug interactions occurring after concomitant administration of irinotecan and ketoconazole by comparing the drug metabolic pathways generated from primitive molecular events. Figure 7 shows an overview of the prototype system. The system consists of the knowledge base for DIO, the inference programs of the Pathway Object Constructor and Drug Interaction Detector, simulation model generator, and simulation engine.

Figure 7

Overview of the prototype system

Overview of the prototype system. User interface (UI), inference program, and visualization program were implemented in Prolog and JAVA. The ontology (DIO) was divided into three OWL files. 1) Class hierarchy and commonly used instances are stored in dio_event.owl. 2) Instances about irinotecan and its pathway are stored in dio_cpt-11.owl. 3) Instances about ketoconazole and its pathway are stored in dio_KCZ.owl.

DIO, shown in Figure 2, was written in OWL-DL; the controlled vocabularies of process and continuant were implemented as OWL class hierarchy, and extension and part of attributes, including molecules and organs, were implemented as OWL instances. The molecular event objects in extension were represented by OWL instances and OWL properties. The ontology referred to other taxonomies and ontologies for well-established vocabularies of biochemical terms, anatomical entities and properties. This enabled the reduction of our ontology construction cost and to concentrate our efforts on the information specific to drug interaction.

The ontology for irinotecan and ketoconazole was built with Protege 31. The information about these drugs has been collected from various articles 172232333435363738. The total number of classes and instances of the current DIO is 178 and 143, respectively. Approximately two thirds of the classes were mapped to the Unified Medical Language System (UMLS) 39 and/or Foundational Model of Anatomy (FMA) 4041; 106 classes were mapped to the UMLS; 32 classes were mapped to the FMA. The class hierarchy is shown in Figure 3. Twenty-six properties were imported from the Open Biomedical Ontology Relation Ontology (OBO Relation) 42; and five properties were newly defined to implement molecular event objects as OWL instances (Table 2). The ontology OWL files: dio_event.owl, dio_cpt-11, and dio_KCZ.owl, are available as additional files 1, 2, 3.

Additional file 1

owl file for dio_event ontology: This file includes all OWL classes and properties defined for the drug interaction ontology (DIO).

Click here for file

Additional file 2

owl file for dio_cpt-11 ontology: This file includes all OWL instances of reactions defined for irinotecan.

Click here for file

Additional file 3

owl file for dio_kcz ontology: This file includes all OWL instances of reactions defined for ketoconazole.

Click here for file

Table 2

The major properties of molecular event objects in OWL-DL.

Properties

Domain (class)

Range (class)

oboRel:relationship

Process, Continuant

-oboRel:located in

independent_entities

Independent_entities, fma:Anatomical_cavity

-oboRel:has_participant

Process

Independent_entities

-- hasResultantParticipant

Process

Independent_entities

-- hasTriggerParticipant

Process

Independent_entities

-- hasSituatedParticipant

Process

Independent_entities

hasSubProcesses

Process

occurred_in

Process

Independent_entities, fma:Anatomical_cavity

The relations among molecular events, attributes (continuants) and drug interactions shown in Figure 2 were implemented as OWL properties. We defined five properties for implementing molecular event objects in OWL-DL in addition to the properties of OBO Relation ontology. Three properties for resultant, trigger, and situated participants were defined as sub properties of oboRel:has_participant, inheriting the definition of the super property. Note that classes for Domain and Range of properties are denoted by the controlled vocabulary defined in Figure 3.

The inference programs of the Pathway Object Constructor and Drug Interaction Detector were implemented with SWI-Prolog 43 and the Semweb library. The Semweb library enables the system to access OWL and RDF files as Prolog clauses. All logical inferences, including pathway generation and drug interaction detection, were implemented in Prolog.

Systems for simulation model generation, and pathway visualization were implemented in JAVA. These systems access DIO through SWI-Prolog semweb library to get detailed information of the generated pathways. The simulation models were generated in the form of SBML, then translated in the form of Octave. Graphviz 44, a graph visualization program designed by Gansner et al., was used to visualize pathways generated by the inference programs. The PPD Viewer designed by Azuma et al. 12 was used for numerical simulation.

Numerical simulation model for co-administration of Irinotecan and Ketoconazole

In order to increase predictive performance, a simplified pathway was used for the generation of simulation models from the viewpoint of the trade off between model complexity and data availability. Reabsorption through small intestine and reactions concerning albumin were omitted due to the lack of information. Furthermore, enzymatic reactions involved metabolism of ketoconazole were integrated into a single enzymatic reaction based on the work of Chien et al. 29.

A simulation model was automatically generated from the simplified pathway. The organs and tissues were integrated into 8 compartments, i.e. blood (including rapidly equilibrating tissues: artery, heart, kidneys, lung, and veins), liver, GI (gastrointestinal consists of large intestine, small intestine, portal vein, and stomach), adipose, NET (non-eliminating tissue such as skin and muscle), GI lumen, bile lumen, and urine. Michaelis-Menten equations were used for all enzymatic reactions. We used a competitive Michaelis-Menten inhibition model for this simulation as used for midazolam and ketoconazole inhibition by Chien et al. 29. For drugs and their metabolites in blood, adipose, NET, GI, and liver, the following equations apply:

d C 1 D d t = { − C L u , r , D f B D C 1 D + ( Q 2 + Q 3 ) C 2 D K p 2 D + Q 5 C 5 D K p 5 D + Q 4 C 4 D K p 4 D − ( Q 2 + Q 3 + Q 4 + Q 5 ) C 1 D } / V 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeWabiqaaaqaaKqbaoaalaaabaGaemizaqMaem4qam0aaSbaaeaacqaIXaqmcqWGebaraeqaaaqaaiabdsgaKjabdsha0baakiabg2da9maaceaabaGaeyOeI0Iaem4qamKaemitaW0aaSbaaSqaaiabdwha1jabcYcaSiabdkhaYjabcYcaSiabdseaebqabaGccqWGMbGzdaWgaaWcbaGaemOqaiKaemiraqeabeaakiabdoeadnaaBaaaleaacqaIXaqmcqWGebaraeqaaOGaey4kaSIaeiikaGIaemyuae1aaSbaaSqaaiabikdaYaqabaGccqGHRaWkcqWGrbqudaWgaaWcbaGaeG4mamdabeaakiabcMcaPKqbaoaalaaabaGaem4qam0aaSbaaeaacqaIYaGmcqWGebaraeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeGOmaiJaemiraqeabeaaaaGccqGHRaWkcqWGrbqudaWgaaWcbaGaeGynaudabeaajuaGdaWcaaqaaiabdoeadnaaBaaabaGaeGynauJaemiraqeabeaaaeaacqWGlbWscqWGWbaCdaWgaaqaaiabiwda1iabdseaebqabaaaaaGccaGL7baaaeaadaGacaqaaiabgUcaRiabdgfarnaaBaaaleaacqaI0aanaeqaaKqbaoaalaaabaGaem4qam0aaSbaaeaacqaI0aancqWGebaraeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeGinaqJaemiraqeabeaaaaGccqGHsislcqGGOaakcqWGrbqudaWgaaWcbaGaeGOmaidabeaakiabgUcaRiabdgfarnaaBaaaleaacqaIZaWmaeqaaOGaey4kaSIaemyuae1aaSbaaSqaaiabisda0aqabaGccqGHRaWkcqWGrbqudaWgaaWcbaGaeGynaudabeaakiabcMcaPiabdoeadnaaBaaaleaacqaIXaqmcqWGebaraeqaaaGccaGL9baacqGGVaWlcqWGwbGvdaWgaaWcbaGaeGymaedabeaaaaaaaa@848F@

d C 4 D d t = ( Q 4 C 1 D − Q 4 C 4 D K p 4 D ) / V 4 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqWGKbazcqWGdbWqdaWgaaqaaiabisda0iabdseaebqabaaabaGaemizaqMaemiDaqhaaOGaeyypa0ZaaeWaaeaacqWGrbqudaWgaaWcbaGaeGinaqdabeaakiabdoeadnaaBaaaleaacqaIXaqmcqWGebaraeqaaOGaeyOeI0Iaemyuae1aaSbaaSqaaiabisda0aqabaqcfa4aaSaaaeaacqWGdbWqdaWgaaqaaiabisda0iabdseaebqabaaabaGaem4saSKaemiCaa3aaSbaaeaacqaI0aancqWGebaraeqaaaaaaOGaayjkaiaawMcaaiabc+caViabdAfawnaaBaaaleaacqaI0aanaeqaaaaa@4B47@

d C 5 D d t = ( Q 5 C 1 D − Q 5 C 5 D K p 5 D ) / V 5 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqWGKbazcqWGdbWqdaWgaaqaaiabiwda1iabdseaebqabaaabaGaemizaqMaemiDaqhaaOGaeyypa0ZaaeWaaeaacqWGrbqudaWgaaWcbaGaeGynaudabeaakiabdoeadnaaBaaaleaacqaIXaqmcqWGebaraeqaaOGaeyOeI0Iaemyuae1aaSbaaSqaaiabiwda1aqabaqcfa4aaSaaaeaacqWGdbWqdaWgaaqaaiabiwda1iabdseaebqabaaabaGaem4saSKaemiCaa3aaSbaaeaacqaI1aqncqWGebaraeqaaaaaaOGaayjkaiaawMcaaiabc+caViabdAfawnaaBaaaleaacqaI1aqnaeqaaaaa@4B53@

d C 3 I M d t = ( Q 3 C 1 I M − Q 3 C 3 I M K p 3 I M ) / V 3 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqWGKbazcqWGdbWqdaWgaaqaaiabiodaZiabdMeajjabd2eanbqabaaabaGaemizaqMaemiDaqhaaOGaeyypa0ZaaeWaaeaacqWGrbqudaWgaaWcbaGaeG4mamdabeaakiabdoeadnaaBaaaleaacqaIXaqmcqWGjbqscqWGnbqtaeqaaOGaeyOeI0Iaemyuae1aaSbaaSqaaiabiodaZaqabaqcfa4aaSaaaeaacqWGdbWqdaWgaaqaaiabiodaZiabdMeajjabd2eanbqabaaabaGaem4saSKaemiCaa3aaSbaaeaacqaIZaWmcqWGjbqscqWGnbqtaeqaaaaaaOGaayjkaiaawMcaaiabc+caViabdAfawnaaBaaaleaacqaIZaWmaeqaaaaa@4FEF@

d C 3 K C Z d t = ( k a K C Z X K C Z , G I L u m e n + Q 3 C 1 K C Z − Q 3 C 3 K C Z K p 3 K C Z ) / V 3 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqcfa4aaSaaaeaacqWGKbazcqWGdbWqdaWgaaqaaiabiodaZiabdUealjabdoeadjabdQfaAbqabaaabaGaemizaqMaemiDaqhaaOGaeyypa0ZaaeWaaeaacqWGRbWAcqWGHbqydaWgaaWcbaGaem4saSKaem4qamKaemOwaOfabeaakiabdIfaynaaBaaaleaacqWGlbWscqWGdbWqcqWGAbGwcqGGSaalcqWGhbWrcqWGjbqscqWGmbatcqWG1bqDcqWGTbqBcqWGLbqzcqWGUbGBaeqaaOGaey4kaSIaemyuae1aaSbaaSqaaiabiodaZaqabaGccqWGdbWqdaWgaaWcbaGaeGymaeJaem4saSKaem4qamKaemOwaOfabeaakiabgkHiTiabdgfarnaaBaaaleaacqaIZaWmaeqaaKqbaoaalaaabaGaem4qam0aaSbaaeaacqaIZaWmcqWGlbWscqWGdbWqcqWGAbGwaeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeG4mamJaem4saSKaem4qamKaemOwaOfabeaaaaaakiaawIcacaGLPaaacqGGVaWlcqWGwbGvdaWgaaWcbaGaeG4mamdabeaaaaa@6A6B@

d C 2 I d t = { Q 2 C 1 I + Q 3 C 3 I k p 3 I − ( Q 2 + Q 3 ) C 2 I k p 2 I − ( C L u , b , I + V m a x a a l p h a a V 2 k m a + f b I C 2 I K p 2 I + V m a x c a l p h a c V 2 k m c ( 1 + f b K C Z C 2 K C Z K p 2 K C Z K i A P C ) + f b I C 2 I K p 2 I + V m a x d a l p h a d V 2 k m d ( 1 + f b K C Z C 2 K C Z K p 2 K C Z K i N P C ) + f b I C 2 I K p 2 I ) f b I C 2 I K p 2 I } / V 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeWabmqaaaqaaKqbaoaalaaabaGaemizaqMaem4qam0aaSbaaeaacqaIYaGmcqWGjbqsaeqaaaqaaiabdsgaKjabdsha0baakiabg2da9maaceaabaGaemyuae1aaSbaaSqaaiabikdaYaqabaGccqWGdbWqdaWgaaWcbaGaeGymaeJaemysaKeabeaakiabgUcaRiabdgfarnaaBaaaleaacqaIZaWmaeqaaKqbaoaalaaabaGaem4qam0aaSbaaeaacqaIZaWmcqWGjbqsaeqaaaqaaiabdUgaRjabdchaWnaaBaaabaGaeG4mamJaemysaKeabeaaaaGccqGHsislcqGGOaakcqWGrbqudaWgaaWcbaGaeGOmaidabeaakiabgUcaRiabdgfarnaaBaaaleaacqaIZaWmaeqaaOGaeiykaKscfa4aaSaaaeaacqWGdbWqdaWgaaqaaiabikdaYiabdMeajbqabaaabaGaem4AaSMaemiCaa3aaSbaaeaacqaIYaGmcqWGjbqsaeqaaaaaaOGaay5EaaaabaGaeyOeI0YaaeqaaeaacqWGdbWqcqWGmbatdaWgaaWcbaGaemyDauNaeiilaWIaemOyaiMaeiilaWIaemysaKeabeaakiabgUcaRKqbaoaalaaabaGaemOvayLaemyBa0MaemyyaeMaemiEaG3aaSbaaeaacqWGHbqyaeqaaiabdggaHjabdYgaSjabdchaWjabdIgaOjabdggaHnaaBaaabaGaemyyaegabeaacqWGwbGvdaWgaaqaaiabikdaYaqabaaabaGaem4AaSMaemyBa02aaSbaaeaacqWGHbqyaeqaaiabgUcaRmaalaaabaGaemOzayMaemOyai2aaSbaaeaacqWGjbqsaeqaaiabdoeadnaaBaaabaGaeGOmaiJaemysaKeabeaaaeaacqWGlbWscqWGWbaCdaWgaaqaaiabikdaYiabdMeajbqabaaaaaaakiabgUcaRKqbaoaalaaabaGaemOvayLaemyBa0MaemyyaeMaemiEaG3aaSbaaeaacqWGJbWyaeqaaiabdggaHjabdYgaSjabdchaWjabdIgaOjabdggaHnaaBaaabaGaem4yamgabeaacqWGwbGvdaWgaaqaaiabikdaYaqabaaabaGaem4AaSMaemyBa02aaSbaaeaacqWGJbWyaeqaaiabcIcaOiabigdaXiabgUcaRmaalaaabaGaemOzayMaemOyai2aaSbaaeaacqWGlbWscqWGdbWqcqWGAbGwaeqaaiabdoeadnaaBaaabaGaeGOmaiJaem4saSKaem4qamKaemOwaOfabeaaaeaacqWGlbWscqWGWbaCdaWgaaqaaiabikdaYiabdUealjabdoeadjabdQfaAbqabaGaem4saSKaemyAaK2aaSbaaeaacqWGbbqqcqWGqbaucqWGdbWqaeqaaaaacqGGPaqkcqGHRaWkdaWcaaqaaiabdAgaMjabdkgaInaaBaaabaGaemysaKeabeaacqWGdbWqdaWgaaqaaiabikdaYiabdMeajbqabaaabaGaem4saSKaemiCaa3aaSbaaeaacqaIYaGmcqWGjbqsaeqaaaaaaaaakiaawIcaaaqaamaaciaabaWaaeGaaeaacqGHRaWkjuaGdaWcaaqaaiabdAfawjabd2gaTjabdggaHjabdIha4naaBaaabaGaemizaqgabeaacqWGHbqycqWGSbaBcqWGWbaCcqWGObaAcqWGHbqydaWgaaqaaiabdsgaKbqabaGaemOvay1aaSbaaeaacqaIYaGmaeqaaaqaaiabdUgaRjabd2gaTnaaBaaabaGaemizaqgabeaacqGGOaakcqaIXaqmcqGHRaWkdaWcaaqaaiabdAgaMjabdkgaInaaBaaabaGaem4saSKaem4qamKaemOwaOfabeaacqWGdbWqdaWgaaqaaiabikdaYiabdUealjabdoeadjabdQfaAbqabaaabaGaem4saSKaemiCaa3aaSbaaeaacqaIYaGmcqWGlbWscqWGdbWqcqWGAbGwaeqaaiabdUealjabdMgaPnaaBaaabaGaemOta4KaemiuaaLaem4qameabeaaaaGaeiykaKIaey4kaSYaaSaaaeaacqWGMbGzcqWGIbGydaWgaaqaaiabdMeajbqabaGaem4qam0aaSbaaeaacqaIYaGmcqWGjbqsaeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeGOmaiJaemysaKeabeaaaaaaaaGccaGLPaaajuaGdaWcaaqaaiabdAgaMjabdkgaInaaBaaabaGaemysaKeabeaacqWGdbWqdaWgaaqaaiabikdaYiabdMeajbqabaaabaGaem4saSKaemiCaa3aaSbaaeaacqaIYaGmcqWGjbqsaeqaaaaaaOGaayzFaaGaei4la8IaemOvay1aaSbaaSqaaiabikdaYaqabaaaaaaa@17F8@

d C 2 A P C d t = { Q 2 C 1 A P C + Q 3 C 3 A P C K p 3 A P C + V m a x c a l p h a c V 2 K m c ( 1 + f b K C Z C 2 K C Z K p 2 K C Z K i A P C ) + f b I C 2 I K p 2 I f b I C 2 I K p 2 I − ( Q 2 + Q 3 ) C 2 A P C K p 2 A P C − C L u , b , A P C f b A P C C 2 A P C K p 2 A P C } / V 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeWabiqaaaqaaKqbaoaalaaabaGaemizaqMaem4qam0aaSbaaeaacqaIYaGmcqWGbbqqcqWGqbaucqWGdbWqaeqaaaqaaiabdsgaKjabdsha0baakiabg2da9maaceaabaGaemyuae1aaSbaaSqaaiabikdaYaqabaGccqWGdbWqdaWgaaWcbaGaeGymaeJaemyqaeKaemiuaaLaem4qameabeaakiabgUcaRiabdgfarnaaBaaaleaacqaIZaWmaeqaaKqbaoaalaaabaGaem4qam0aaSbaaeaacqaIZaWmcqWGbbqqcqWGqbaucqWGdbWqaeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeG4mamJaemyqaeKaemiuaaLaem4qameabeaaaaGccqGHRaWkjuaGdaWcaaqaaiabdAfawjabd2gaTjabdggaHjabdIha4naaBaaabaGaem4yamgabeaacqWGHbqycqWGSbaBcqWGWbaCcqWGObaAcqWGHbqydaWgaaqaaiabdogaJbqabaGaemOvay1aaSbaaeaacqaIYaGmaeqaaaqaaiabdUealjabd2gaTnaaBaaabaGaem4yamgabeaacqGGOaakcqaIXaqmcqGHRaWkdaWcaaqaaiabdAgaMjabdkgaInaaBaaabaGaem4saSKaem4qamKaemOwaOfabeaacqWGdbWqdaWgaaqaaiabikdaYiabdUealjabdoeadjabdQfaAbqabaaabaGaem4saSKaemiCaa3aaSbaaeaacqaIYaGmcqWGlbWscqWGdbWqcqWGAbGwaeqaaiabdUealjabdMgaPnaaBaaabaGaemyqaeKaemiuaaLaem4qameabeaaaaGaeiykaKIaey4kaSYaaSaaaeaacqWGMbGzcqWGIbGydaWgaaqaaiabdMeajbqabaGaem4qam0aaSbaaeaacqaIYaGmcqWGjbqsaeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeGOmaiJaemysaKeabeaaaaaaamaalaaabaGaemOzayMaemOyai2aaSbaaeaacqWGjbqsaeqaaiabdoeadnaaBaaabaGaeGOmaiJaemysaKeabeaaaeaacqWGlbWscqWGWbaCdaWgaaqaaiabikdaYiabdMeajbqabaaaaaGccaGL7baaaeaadaGacaqaaiabgkHiTiabcIcaOiabdgfarnaaBaaaleaacqaIYaGmaeqaaOGaey4kaSIaemyuae1aaSbaaSqaaiabiodaZaqabaGccqGGPaqkjuaGdaWcaaqaaiabdoeadnaaBaaabaGaeGOmaiJaemyqaeKaemiuaaLaem4qameabeaaaeaacqWGlbWscqWGWbaCdaWgaaqaaiabikdaYiabdgeabjabdcfaqjabdoeadbqabaaaaOGaeyOeI0Iaem4qamKaemitaW0aaSbaaSqaaiabdwha1jabcYcaSiabdkgaIjabcYcaSiabdgeabjabdcfaqjabdoeadbqabaqcfa4aaSaaaeaacqWGMbGzcqWGIbGydaWgaaqaaiabdgeabjabdcfaqjabdoeadbqabaGaem4qam0aaSbaaeaacqaIYaGmcqWGbbqqcqWGqbaucqWGdbWqaeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeGOmaiJaemyqaeKaemiuaaLaem4qameabeaaaaaakiaaw2haaiabc+caViabdAfawnaaBaaaleaacqaIYaGmaeqaaaaaaaa@D350@

d C 2 N P C d t = { Q 2 C 1 N P C + Q 3 C 3 N P C K p 3 N P C + V m a x d a l p h a d V 2 K m d ( 1 + f b K C Z C 2 K C Z K p 2 K C Z K i N P C ) + f b I C 2 I K p 2 I f b I C 2 I K p 2 I − ( Q 2 + Q 3 ) C 2 N P C K p 2 N A P C − ( C L u , b , N P C V m a x b a l p h a b V 2 K m b + f b N P C C 2 N P C K p 2 N P C ) ) f b N P C C 2 N P C K p 2 N P C } / V 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeWabiqaaaqaaKqbaoaalaaabaGaemizaqMaem4qam0aaSbaaeaacqaIYaGmcqWGobGtcqWGqbaucqWGdbWqaeqaaaqaaiabdsgaKjabdsha0baakiabg2da9maaceaabaGaemyuae1aaSbaaSqaaiabikdaYaqabaGccqWGdbWqdaWgaaWcbaGaeGymaeJaemOta4KaemiuaaLaem4qameabeaakiabgUcaRiabdgfarnaaBaaaleaacqaIZaWmaeqaaKqbaoaalaaabaGaem4qam0aaSbaaeaacqaIZaWmcqWGobGtcqWGqbaucqWGdbWqaeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeG4mamJaemOta4KaemiuaaLaem4qameabeaaaaGccqGHRaWkjuaGdaWcaaqaaiabdAfawjabd2gaTjabdggaHjabdIha4naaBaaabaGaemizaqgabeaacqWGHbqycqWGSbaBcqWGWbaCcqWGObaAcqWGHbqydaWgaaqaaiabdsgaKbqabaGaemOvay1aaSbaaeaacqaIYaGmaeqaaaqaaiabdUealjabd2gaTnaaBaaabaGaemizaqgabeaacqGGOaakcqaIXaqmcqGHRaWkdaWcaaqaaiabdAgaMjabdkgaInaaBaaabaGaem4saSKaem4qamKaemOwaOfabeaacqWGdbWqdaWgaaqaaiabikdaYiabdUealjabdoeadjabdQfaAbqabaaabaGaem4saSKaemiCaa3aaSbaaeaacqaIYaGmcqWGlbWscqWGdbWqcqWGAbGwaeqaaiabdUealjabdMgaPnaaBaaabaGaemOta4KaemiuaaLaem4qameabeaaaaGaeiykaKIaey4kaSYaaSaaaeaacqWGMbGzcqWGIbGydaWgaaqaaiabdMeajbqabaGaem4qam0aaSbaaeaacqaIYaGmcqWGjbqsaeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeGOmaiJaemysaKeabeaaaaaaamaalaaabaGaemOzayMaemOyai2aaSbaaeaacqWGjbqsaeqaaiabdoeadnaaBaaabaGaeGOmaiJaemysaKeabeaaaeaacqWGlbWscqWGWbaCdaWgaaqaaiabikdaYiabdMeajbqabaaaaaGccaGL7baaaeaadaGacaqaaiabgkHiTiabcIcaOiabdgfarnaaBaaaleaacqaIYaGmaeqaaOGaey4kaSIaemyuae1aaSbaaSqaaiabiodaZaqabaGccqGGPaqkjuaGdaWcaaqaaiabdoeadnaaBaaabaGaeGOmaiJaemOta4KaemiuaaLaem4qameabeaaaeaacqWGlbWscqWGWbaCdaWgaaqaaiabikdaYiabd6eaojabdgeabjabdcfaqjabdoeadbqabaaaaOGaeyOeI0YaaeWaaeaacqWGdbWqcqWGmbatdaWgaaWcbaGaemyDauNaeiilaWIaemOyaiMaeiilaWIaemOta4KaemiuaaLaem4qameabeaajuaGdaWcaaqaaiabdAfawjabd2gaTjabdggaHjabdIha4naaBaaabaGaemOyaigabeaacqWGHbqycqWGSbaBcqWGWbaCcqWGObaAcqWGHbqydaWgaaqaaiabdkgaIbqabaGaemOvay1aaSbaaeaacqaIYaGmaeqaaaqaaiabdUealjabd2gaTnaaBaaabaGaemOyaigabeaacqGHRaWkdaWcaaqaaiabdAgaMjabdkgaInaaBaaabaGaemOta4KaemiuaaLaem4qameabeaacqWGdbWqdaWgaaqaaiabikdaYiabd6eaojabdcfaqjabdoeadbqabaaabaGaem4saSKaemiCaa3aaSbaaeaacqaIYaGmcqWGobGtcqWGqbaucqWGdbWqaeqaaaaacqGGPaqkaaaakiaawIcacaGLPaaajuaGdaWcaaqaaiabdAgaMjabdkgaInaaBaaabaGaemOta4KaemiuaaLaem4qameabeaacqWGdbWqdaWgaaqaaiabikdaYiabd6eaojabdcfaqjabdoeadbqabaaabaGaem4saSKaemiCaa3aaSbaaeaacqaIYaGmcqWGobGtcqWGqbaucqWGdbWqaeqaaaaaaOGaayzFaaGaei4la8IaemOvay1aaSbaaSqaaiabikdaYaqabaaaaaaa@013C@

d C 2 S N − 38 d t = { Q 2 C 1 S N − 38 + Q 3 C 3 S N − 38 K p 3 S N − 38 − ( Q 2 + Q 3 ) C 2 S N − 38 K p 2 S N − 38 + V m a x a a l p h a a V 2 K m a + f b I C 2 I K p 2 I f b I C 2 I K p 2 I + V m a x b a l p h a b V 2 K m b + f b N P C C 2 N P C K p 2 N P C f b N P C C 2 N P C K p 2 N P C − ( C L u , b , S N − 38 + V m a x f a l p h a f V 2 K m f + f b S N − 38 C 2 S N − 38 K p 2 S N − 38 ) f b S N − 38 C 2 S N − 38 K p 2 S N − 38 } / V 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeWabmqaaaqaaKqbaoaalaaabaGaemizaqMaem4qam0aaSbaaeaacqaIYaGmcqWGtbWucqWGobGtcqGHsislcqaIZaWmcqaI4aaoaeqaaaqaaiabdsgaKjabdsha0baakiabg2da9maaceaabaGaemyuae1aaSbaaSqaaiabikdaYaqabaGccqWGdbWqdaWgaaWcbaGaeGymaeJaem4uamLaemOta4KaeyOeI0IaeG4mamJaeGioaGdabeaakiabgUcaRiabdgfarnaaBaaaleaacqaIZaWmaeqaaKqbaoaalaaabaGaem4qam0aaSbaaeaacqaIZaWmcqWGtbWucqWGobGtcqGHsislcqaIZaWmcqaI4aaoaeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeG4mamJaem4uamLaemOta4KaeyOeI0IaeG4mamJaeGioaGdabeaaaaGccqGHsislcqGGOaakcqWGrbqudaWgaaWcbaGaeGOmaidabeaakiabgUcaRiabdgfarnaaBaaaleaacqaIZaWmaeqaaOGaeiykaKscfa4aaSaaaeaacqWGdbWqdaWgaaqaaiabikdaYiabdofatjabd6eaojabgkHiTiabiodaZiabiIda4aqabaaabaGaem4saSKaemiCaa3aaSbaaeaacqaIYaGmcqWGtbWucqWGobGtcqGHsislcqaIZaWmcqaI4aaoaeqaaaaaaOGaay5EaaaabaGaey4kaSscfa4aaSaaaeaacqWGwbGvcqWGTbqBcqWGHbqycqWG4baEdaWgaaqaaiabdggaHbqabaGaemyyaeMaemiBaWMaemiCaaNaemiAaGMaemyyae2aaSbaaeaacqWGHbqyaeqaaiabdAfawnaaBaaabaGaeGOmaidabeaaaeaacqWGlbWscqWGTbqBdaWgaaqaaiabdggaHbqabaGaey4kaSYaaSaaaeaacqWGMbGzcqWGIbGydaWgaaqaaiabdMeajbqabaGaem4qam0aaSbaaeaacqaIYaGmcqWGjbqsaeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeGOmaiJaemysaKeabeaaaaaaamaalaaabaGaemOzayMaemOyai2aaSbaaeaacqWGjbqsaeqaaiabdoeadnaaBaaabaGaeGOmaiJaemysaKeabeaaaeaacqWGlbWscqWGWbaCdaWgaaqaaiabikdaYiabdMeajbqabaaaaOGaey4kaSscfa4aaSaaaeaacqWGwbGvcqWGTbqBcqWGHbqycqWG4baEdaWgaaqaaiabdkgaIbqabaGaemyyaeMaemiBaWMaemiCaaNaemiAaGMaemyyae2aaSbaaeaacqWGIbGyaeqaaiabdAfawnaaBaaabaGaeGOmaidabeaaaeaacqWGlbWscqWGTbqBdaWgaaqaaiabdkgaIbqabaGaey4kaSYaaSaaaeaacqWGMbGzcqWGIbGydaWgaaqaaiabd6eaojabdcfaqjabdoeadbqabaGaem4qam0aaSbaaeaacqaIYaGmcqWGobGtcqWGqbaucqWGdbWqaeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeGOmaiJaemOta4KaemiuaaLaem4qameabeaaaaaaamaalaaabaGaemOzayMaemOyai2aaSbaaeaacqWGobGtcqWGqbaucqWGdbWqaeqaaiabdoeadnaaBaaabaGaeGOmaiJaemOta4KaemiuaaLaem4qameabeaaaeaacqWGlbWscqWGWbaCdaWgaaqaaiabikdaYiabd6eaojabdcfaqjabdoeadbqabaaaaaGcbaWaaiGaaeaacqGHsisldaqadaqaaiabdoeadjabdYeamnaaBaaaleaacqWG1bqDcqGGSaalcqWGIbGycqGGSaalcqWGtbWucqWGobGtcqGHsislcqaIZaWmcqaI4aaoaeqaaOGaey4kaSscfa4aaSaaaeaacqWGwbGvcqWGTbqBcqWGHbqycqWG4baEdaWgaaqaaiabdAgaMbqabaGaemyyaeMaemiBaWMaemiCaaNaemiAaGMaemyyae2aaSbaaeaacqWGMbGzaeqaaiabdAfawnaaBaaabaGaeGOmaidabeaaaeaacqWGlbWscqWGTbqBdaWgaaqaaiabdAgaMbqabaGaey4kaSYaaSaaaeaacqWGMbGzcqWGIbGydaWgaaqaaiabdofatjabd6eaojabgkHiTiabiodaZiabiIda4aqabaGaem4qam0aaSbaaeaacqaIYaGmcqWGtbWucqWGobGtcqGHsislcqaIZaWmcqaI4aaoaeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeGOmaiJaem4uamLaemOta4KaeyOeI0IaeG4mamJaeGioaGdabeaaaaaaaaGccaGLOaGaayzkaaqcfa4aaSaaaeaacqWGMbGzcqWGIbGydaWgaaqaaiabdofatjabd6eaojabgkHiTiabiodaZiabiIda4aqabaGaem4qam0aaSbaaeaacqaIYaGmcqWGtbWucqWGobGtcqGHsislcqaIZaWmcqaI4aaoaeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeGOmaiJaem4uamLaemOta4KaeyOeI0IaeG4mamJaeGioaGdabeaaaaaakiaaw2haaiabc+caViabdAfawnaaBaaaleaacqaIYaGmaeqaaaaaaaa@38BF@

d C 2 S N − 38 G d t = { Q 2 C 1 S N − 38 G + Q 3 C 3 S N − 38 K p 3 S N − 38 + V m a x f a l p h a f V 2 K m f + f b S N − 38 C 2 S N − 38 K p 2 S N − 38 f b S N − 38 C 2 S N − 38 K p 2 S N − 38 − ( Q 2 + Q 3 ) C 2 S N − 38 G K p 2 S N − 38 G − C L u , b , S N − 38 G f b S N − 38 G C 2 S N − 38 G K p 2 S N − 38 G } / V 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeWabiqaaaqaaKqbaoaalaaabaGaemizaqMaem4qam0aaSbaaeaacqaIYaGmcqWGtbWucqWGobGtcqGHsislcqaIZaWmcqaI4aaocqWGhbWraeqaaaqaaiabdsgaKjabdsha0baakiabg2da9maaceaabaGaemyuae1aaSbaaSqaaiabikdaYaqabaGccqWGdbWqdaWgaaWcbaGaeGymaeJaem4uamLaemOta4KaeyOeI0IaeG4mamJaeGioaGJaem4raCeabeaakiabgUcaRiabdgfarnaaBaaaleaacqaIZaWmaeqaaKqbaoaalaaabaGaem4qam0aaSbaaeaacqaIZaWmcqWGtbWucqWGobGtcqGHsislcqaIZaWmcqaI4aaoaeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeG4mamJaem4uamLaemOta4KaeyOeI0IaeG4mamJaeGioaGdabeaaaaGccqGHRaWkjuaGdaWcaaqaaiabdAfawjabd2gaTjabdggaHjabdIha4naaBaaabaGaemOzaygabeaacqWGHbqycqWGSbaBcqWGWbaCcqWGObaAcqWGHbqydaWgaaqaaiabdAgaMbqabaGaemOvay1aaSbaaeaacqaIYaGmaeqaaaqaaiabdUealjabd2gaTnaaBaaabaGaemOzaygabeaacqGHRaWkdaWcaaqaaiabdAgaMjabdkgaInaaBaaabaGaem4uamLaemOta4KaeyOeI0IaeG4mamJaeGioaGdabeaacqWGdbWqdaWgaaqaaiabikdaYiabdofatjabd6eaojabgkHiTiabiodaZiabiIda4aqabaaabaGaem4saSKaemiCaa3aaSbaaeaacqaIYaGmcqWGtbWucqWGobGtcqGHsislcqaIZaWmcqaI4aaoaeqaaaaaaaWaaSaaaeaacqWGMbGzcqWGIbGydaWgaaqaaiabdofatjabd6eaojabgkHiTiabiodaZiabiIda4aqabaGaem4qam0aaSbaaeaacqaIYaGmcqWGtbWucqWGobGtcqGHsislcqaIZaWmcqaI4aaoaeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeGOmaiJaem4uamLaemOta4KaeyOeI0IaeG4mamJaeGioaGdabeaaaaaakiaawUhaaaqaamaaciaabaGaeyOeI0IaeiikaGIaemyuae1aaSbaaSqaaiabikdaYaqabaGccqGHRaWkcqWGrbqudaWgaaWcbaGaeG4mamdabeaakiabcMcaPKqbaoaalaaabaGaem4qam0aaSbaaeaacqaIYaGmcqWGtbWucqWGobGtcqGHsislcqaIZaWmcqaI4aaocqWGhbWraeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeGOmaiJaem4uamLaemOta4KaeyOeI0IaeG4mamJaeGioaGJaem4raCeabeaaaaGccqGHsislcqWGdbWqcqWGmbatdaWgaaWcbaGaemyDauNaeiilaWIaemOyaiMaeiilaWIaem4uamLaemOta4KaeyOeI0IaeG4mamJaeGioaGJaem4raCeabeaajuaGdaWcaaqaaiabdAgaMjabdkgaInaaBaaabaGaem4uamLaemOta4KaeyOeI0IaeG4mamJaeGioaGJaem4raCeabeaacqWGdbWqdaWgaaqaaiabikdaYiabdofatjabd6eaojabgkHiTiabiodaZiabiIda4iabdEeahbqabaaabaGaem4saSKaemiCaa3aaSbaaeaacqaIYaGmcqWGtbWucqWGobGtcqGHsislcqaIZaWmcqaI4aaocqWGhbWraeqaaaaaaOGaayzFaaGaei4la8IaemOvay1aaSbaaSqaaiabikdaYaqabaaaaaaa@EBD7@

d C 2 K C Z d t = { Q 2 C 1 K C Z + Q 3 C 3 K C Z K p 3 K C Z − ( Q 2 + Q 3 ) C 2 K C Z K p 2 K C Z − V m a x e a l p h a e V 2 K m e + f b K C Z C 2 K C Z K p 2 K C Z f b K C Z C 2 K C Z K p 2 K C Z } / V 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeqabiqaaaqaaKqbaoaalaaabaGaemizaqMaem4qam0aaSbaaeaacqaIYaGmcqWGlbWscqWGdbWqcqWGAbGwaeqaaaqaaiabdsgaKjabdsha0baakiabg2da9maaceaabaGaemyuae1aaSbaaSqaaiabikdaYaqabaGccqWGdbWqdaWgaaWcbaGaeGymaeJaem4saSKaem4qamKaemOwaOfabeaakiabgUcaRiabdgfarnaaBaaaleaacqaIZaWmaeqaaKqbaoaalaaabaGaem4qam0aaSbaaeaacqaIZaWmcqWGlbWscqWGdbWqcqWGAbGwaeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeG4mamJaem4saSKaem4qamKaemOwaOfabeaaaaGccqGHsislcqGGOaakcqWGrbqudaWgaaWcbaGaeGOmaidabeaakiabgUcaRiabdgfarnaaBaaaleaacqaIZaWmaeqaaOGaeiykaKscfa4aaSaaaeaacqWGdbWqdaWgaaqaaiabikdaYiabdUealjabdoeadjabdQfaAbqabaaabaGaem4saSKaemiCaa3aaSbaaeaacqaIYaGmcqWGlbWscqWGdbWqcqWGAbGwaeqaaaaaaOGaay5EaaaabaWaaiGaaeaacqGHsisljuaGdaWcaaqaaiabdAfawjabd2gaTjabdggaHjabdIha4naaBaaabaGaemyzaugabeaacqWGHbqycqWGSbaBcqWGWbaCcqWGObaAcqWGHbqydaWgaaqaaiabdwgaLbqabaGaemOvay1aaSbaaeaacqaIYaGmaeqaaaqaaiabdUealjabd2gaTnaaBaaabaGaemyzaugabeaacqGHRaWkdaWcaaqaaiabdAgaMjabdkgaInaaBaaabaGaem4saSKaem4qamKaemOwaOfabeaacqWGdbWqdaWgaaqaaiabikdaYiabdUealjabdoeadjabdQfaAbqabaaabaGaem4saSKaemiCaa3aaSbaaeaacqaIYaGmcqWGlbWscqWGdbWqcqWGAbGwaeqaaaaaaaWaaSaaaeaacqWGMbGzcqWGIbGydaWgaaqaaiabdUealjabdoeadjabdQfaAbqabaGaem4qam0aaSbaaeaacqaIYaGmcqWGlbWscqWGdbWqcqWGAbGwaeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeGOmaiJaem4saSKaem4qamKaemOwaOfabeaaaaaakiaaw2haaiabc+caViabdAfawnaaBaaaleaacqaIYaGmaeqaaaaaaaa@A805@

Here Qp_n, and V_n(n = 1 to 5) are blood flows and volume of compartments, respectively. The subscripts n stands for compartments: n = 1, 2, 3, 4 and 5 represents blood, liver, GI, adipose, and NET, respectively. Kp_ns, and C_ns(n = 1 to 5, s = D, I, KCZ, IM, APC, NPC, SN-38, and SN-38G) are tissue-to-blood concentration ratios, and blood concentrations, respectively. The subscripts s stands for drugs and their metabolites: s = D, I, KCZ, IM, APC, NPC, SN-38, and SN-38G represents drugs and metabolites, irinotecan, ketoconazole, metabolites of irinotecan, APC, NPC, SN-38, and SN-38, respectively. CL_u,r,sCL_u,b,s, and f_Bs(s = D, I, KCZ, IM, APC, NPC, SN-38, and SN-38G) are renal clearance for unbound drugs, bile clearance for unbound drugs, and blood unbound fraction, respectively. ka_KCZand X_KCZ;GILumenare absorption rate constant of ketoconazole and amount of ketoconazole in the GI luminal compartment, where ketoconazole is absorbed. V max_m, km_m, and α_m(m = a, b, c, d, e, and f) are Vmax, Michaelis constant, and expression amount of enzyme, respectively. The subscripts m stands for enzymatic reactions: m = a, b, c, d, e, and f represents metabolism of irinotecan to form SN-38 by CE, metabolism of NPC to form SN-38 by CE, oxidation of irinotecan to form APC by CYP3A4, oxidation of irinotecan to form NPC by CYP3A4, glucronidation of SN-38 to form SN-38G by UGT1A1, and metabolism of ketoconazole by CYP3A4, respectively. Ki_APCand Ki_NPCare inhibition constants. During the drip infusion of irinotecan, the following equation apply for blood concentration of irinotecan.

d C 1 I d t = { k 0 I − C L u , r , I f B I C 1 I + ( Q 2 + Q 3 ) C 2 I K p 2 I + Q 5 C 5 I K p 5 I + Q 4 C 4 I K p 4 I − ( Q 2 + Q 3 + Q 4 + Q 5 ) C 1 I } / V 1 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeqabiqaaaqaaKqbaoaalaaabaGaemizaqMaem4qam0aaSbaaeaacqaIXaqmcqWGjbqsaeqaaaqaaiabdsgaKjabdsha0baakiabg2da9maaceaabaGaem4AaSMaeGimaaZaaSbaaSqaaiabdMeajbqabaGccqGHsislcqWGdbWqcqWGmbatdaWgaaWcbaGaemyDauNaeiilaWIaemOCaiNaeiilaWIaemysaKeabeaakiabdAgaMnaaBaaaleaacqWGcbGqcqWGjbqsaeqaaOGaem4qam0aaSbaaSqaaiabigdaXiabdMeajbqabaGccqGHRaWkcqGGOaakcqWGrbqudaWgaaWcbaGaeGOmaidabeaakiabgUcaRiabdgfarnaaBaaaleaacqaIZaWmaeqaaOGaeiykaKscfa4aaSaaaeaacqWGdbWqdaWgaaqaaiabikdaYiabdMeajbqabaaabaGaem4saSKaemiCaa3aaSbaaeaacqaIYaGmcqWGjbqsaeqaaaaakiabgUcaRiabdgfarnaaBaaaleaacqaI1aqnaeqaaKqbaoaalaaabaGaem4qam0aaSbaaeaacqaI1aqncqWGjbqsaeqaaaqaaiabdUealjabdchaWnaaBaaabaGaeGynauJaemysaKeabeaaaaaakiaawUhaaaqaamaaciaabaGaey4kaSIaemyuae1aaSbaaSqaaiabisda0aqabaqcfa4aaSaaaeaacqWGdbWqdaWgaaqaaiabisda0iabdMeajbqabaaabaGaem4saSKaemiCaa3aaSbaaeaacqaI0aancqWGjbqsaeqaaaaakiabgkHiTiabcIcaOiabdgfarnaaBaaaleaacqaIYaGmaeqaaOGaey4kaSIaemyuae1aaSbaaSqaaiabiodaZaqabaGccqGHRaWkcqWGrbqudaWgaaWcbaGaeGinaqdabeaakiabgUcaRiabdgfarnaaBaaaleaacqaI1aqnaeqaaOGaeiykaKIaem4qam0aaSbaaSqaaiabigdaXiabdMeajbqabaaakiaaw2haaiabc+caViabdAfawnaaBaaaleaacqaIXaqmaeqaaaaaaaa@8899@

Here k0_Iis a constant for the velocity of drip infusion.

The values of physiologic parameters, kinetic parameters, and inhibition constants and administration parameters for the simulation model are shown in Table 3, 4, 5 and 6 respectively. The values of these parameters were based on published papers except enzyme expression parameters (α) and tissue-to-blood concentration ratios (Kp). The values for α and Kp were obtained by fitting the simulation result to the experimental data published by Slatter et al. 18.

Table 3

Kinetic parameters for the protein binding, tissue distribution, and urinary and biliary excretion of ketoconazole, irinotecan, and their metabolites.

Compound

f _B

Kp _GI

Kp _Liver

Kp _Adipose

Kp _NET

CL_U,R[mL/min/kg]

CL_U,B[mL/min/kg]

Irinotecan

0.37

1.00

10.0

3.00

6.15

10.6

APC

0.37

1.00

1.50

0.06

1.47

5.45

NPC

0.37

1.00

6.00

2.00

1.49

14.5

SN-38

0.05

1.00

2.00

0.70

9.91

103

SN-38G

1.00

2.80

0.08

1.44

2.03

KCZ

0.01

1.00

15.0

1.50

130

Blood unbound fraction (f_B) values for irinotecan and SN-38 were obtained from the published paper by Oliver et al. [45]. The f_Bvalues for NPC and APC were assumed to be the same as that of Irinotecan. f_Bvalue for SN-38G was assumed to be 1.00. Urinary and biliary clearances (CL_U,R, CL_U,B) were defined for unbound blood concentration. CL_U,Rand CL_U,Bvalues for irinotecan, APC, SN-38, SN-38G were calculated using data from the publish paper [18]. Kp values were determined so that the simulation concentration/time profile fits to the experimental data from Slatter et al. [18].

Table 4

Kinetic parameters for metabolic enzymes.

Enzyme

Substrate

Product

K_m[μM]

V_max[pmol/min/mg protein]

α[mg protein/g tissue]

V_max·α/K_m[mL/min/g tissue]

Carboxylesterase

Irinotecan

SN-38

2.30

2.11

128

0.117

Carboxylesterase

NPC

SN-38

2.30

2.11

128

0.117

CYP3A4

Irinotecan

APC

18.4

26.0

73.3

0.104

CYP3A4

Irinotecan

NPC

48.2

74.1

11.7

0.0180

UGT1A1

SN-38

SN-38G

3.80

50.8

750

10.0

CYP3A4

KCZ

MOK

0.00810

12.5

44.3

68.4

Because of the lack of experimental data, K_mand V_maxvalues for metabolism of NPC to SN-38 by CE were assumed to be the same as those for metabolism of irinotecan to SN-38 by CE. K_mand V_maxvalues for other reactions were obtained from the published papers [25–29]. Values of unknown parameter α were determined so that the simulation concentration/time profile fit to the experimental data published by Slatter et al. [18]. The relatively high fitted value for α in the glucronidation of SN-38 may result from the increase in UGT1A1 expression as drug response. The relatively low value for α in the oxidation of irinotecan to NPC by CYP3A4 suggests a competition in this reaction with the oxidation of irinotecan to APC by CYP3A4.

Table 5

Physiologic parameters.

Organ

Blood flow rate [mL/min/kg]

Volume [mL/kg]

Blood

61.1

51.0

Liver

5.79

32.3

13.4

32.1

Adipose

4.45

204

NET

37.4

681

The physiologic parameters were obtained from the work by Willmann et al. [30]. These values assumes Caucasian male with body weight of 73 kg and height of 176 cm.

Competing interests

The authors have received no financial support or other benefits from commercial sources for the work reported in the manuscript, and no other financial interests that any of the authors may have could create a potential conflict of interest or the appearance of a conflict of interest with regard to the work.

Authors' contributions

Takeshi Arikuma carried out the machine inference studies, participated in Drug Interaction Ontology development, participated in numerical simulation studies and drafted the manuscript. Sumi Yoshikawa carried out the ontology studies, participated in machine inference studies, and helped to draft the manuscript. Ryuzo Azuma carried out the numerical simulation and visualization studies. Kazumi Matsumura extracted knowledge from the articles for the Drug Interaction Ontology when she worked for RIKEN GSC. Kentaro Watanabe developed the original drug-interaction prediction system when he studied in TITECH. Akihiko Konagaya conceived of the study, partcipated in its design and coordination, and helped to draft the manuscript.

Table 6

Parameters for drug interactions and drug administrations.

Parameters

Values

Ki _APC

1.40 * 10^-3

Ki _NPC

5.94 * 10^-4

ka_KCZ[min^-1]

1.83 * 10^-2

dose_Irinotecan[nmol/kg]

4.86 * 10³

t_{dose,Irinotecan}[min]

dose_KCZ[nmol/kg]

5.16 * 10³

Ki_APCand Ki_NPCwere calculated from the experimental data reported by Haaz et al. [25,26]. ka_KCZwas obtained the from published paper [29].

Acknowledgements

The authors give special thanks to Dr. Fumikazu Konishi, Mr. Shingo Ohki and Mr. Ryo Umetsu in RIKEN GSC for their contributions in providing the CONDOR pool and visualization tools on the Open Bioinformatics Grid (OBIGrid).

The computations in this work were performed using the Open Bioinformatics Grid (OBIGrid) 13.

This article has been published as part of BMC Bioinformatics Volume 9 Supplement 6, 2008: Symposium of Computations in Bioinformatics and Bioscience (SCBB07). The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/9?issue=S6.

A possible mechanism of eighteen patient deaths caused by interactions of sorivudine, a new antiviral drug, with oral 5-fluorouracil prodrugs Okuda H Ogura K Kato A Takubo H Watanabe T J Pharmacol Exp Ther 1998 287 2 791 799 9808711 Baker CJ Cheung KH Semantic Web: Revolutionizing Knowledge Discovery in the Life Sciences Springer 2007 OBIGrid: towards the Ba for sharing resources, services and knowledge for bioinformatics Konagaya A 4th International Workshop on Biomedical Computations on the Grid (BioGrid'06), Singapore 2006 A physiologically based pharmacokinetic analysis of capecitabine, a triple prodrug of 5-FU, in humans, the mechanism for tumor-selective accumulation of 5-FU Tsukamoto Y Kato Y Ura M Horii I Ishitsuka H Kusuhara H Sugiyama Y Pharm Res 2001 18 8 1190 1202 11587492 Dynamically simulating the interaction of midazolam and the CYP3A4 inhibitor itraconazole using individual coupled whole-body physiologically-based pharmacokinetic (WB-PBPK) models Vossen M Sevestre M Niederalt C Jang IJ Willmann S Edginton AN Theoretical Biology and Medical Modelling 2007 4 13 13 Ontology-Driven Hypothetic Assertion (OHA) for Drug InteractionPrediction Arikuma T Yoshikawa S Watanabe K Matsumura K Konagaya A Second International Multi-Symposiums on Computer and Computational Sciences (IMSCCS 2007) IEEE Ni J, Cremer J, Seidel E, Dongarra J, Casavant TL, Huang S, Gu G 2007 1 8 Drug interaction ontology (DIO) for inferences of possible drug-drug interactions Yoshikawa S Kenji S Konagaya A Medinfo 2004 11 454 458 15360854 The Semantic Web Berners-Lee T Hendler J Nature 2001 410 6832 1023 1024 11323639 SNAP and SPAN: Towards Dynamic Spatial Ontology Grenon P Smith B Spatial Cognition and Computation 2004 1 69 103 Biodynamic Ontology:Applying BFO in the Biomedical Domain Grenon P Smith B Goldberg L Ontologies in Medicine Amsterdam: IOS Press 2004 20 38 KEGG: Kyoto Encyclopedia of Genes and Genomes Kanehisa M Goto S Nucleic Acids Research 2000 28 27 30 102409 10592173 Discovering dynamic characteristics of biochemical pathways using geometric patterns among parameter-parameter dependencies in differential equations Azuma R Umetsu R Ohki S Konishi F Yoshikawa S Konagaya A Matsumura K New Generation Computing 2007 25 425 441 The superstructure toward Open Bioinformatics Grid Konagaya A Konishi F Hatakeyama M Satou K New Generation Computing 2004 22 167 176 Trends in life science grid: from computing grid to knowledge grid Konagaya A BMC Bioinformatics 2006 7 Suppl 5 S10 1764466 17254294 Konagaya A Satou K Eds Grid Computing in Life Science: First International Workshop on Life Science Grid, LSGRID2004 Splinger 2005 Tan TW Arzberger P Konagaya A Eds Grid Computing in Life Sciences : Proceedings of the 2nd International Life Science Grid Workshop, LSGRID 2005 World Scientific Pub 2006 Clinical Pharmacokinetics and Metabolism of Irinotecan (CPT-11) Mathijssen RHJ van Alphen RJ Verweij J Loos WJ Nooter K Stoter G Sparreboom A Clinical Cancer Research 2001 7 2182 2194 11489791 Pharmacokinetics, Metabolism, and Excretion of Irinotecan (CPT-11) Following I.V. Infusion of [14C]CPT-11 in Cancer Patients Slatter JG Schaaf LJ Sams JP Feenstra KL Johnson MG Bombardt PA Cathcart KS Verburg MT Pearson LK Compton LD Miller LL Baker DS Pesheck CV Raymond S Lord I Drug Metabolism and Disposition 2000 28 4 423 433 UGT1A1 Haplotypes Associated with Reduced Glucuronidation and Increased Serum Bilirubin in Irinotecan-administered Japanese Patients with Cancer Sai K Saeki M Saito Y Ozawa S Katori N Jinno H Hasegawa R Kaniwa N ichi Sawada J Kazuo Komamura KU Kamakura S Kitakaze M Yutaka Kitamura NK Minami H Ohtsu A Shirao K Yoshida T Saijo N Clin Pharmacol Ther 2004 75 501 515 15179405 Polymorphisms of UDP-Glucuronosyltransferase Gene and Irinotecan Toxicity: A Pharmacogenetic Analysis Ando Y Saka H Ando M Sawa T Muro K Ueoka H Yokoyama A Saitoh S Shimokata K Hasegawa Y Cancer Research 2000 60 6921 6926 11156391 Mouse hepatic metabolites of ketoconazole: Isolation and structure elucidation Whitehouse LW Information C Menzies A Dawson B Cyr TD By AW Black DB Zamecnik J J Pharm Biomed Anal 1994 12 11 1425 1441 7849138 Ketoconazole: A review of its Therapeutic Efficacy in Superficial and Systemic Funfal Infections Heel R Brogden R A Carmine PM Speight T Avery G Drugs 1982 23 1 36 6276122 Modulation of Irinotecan Metabolism by Ketoconazole Kehrer DF Mathijssen RH Verweij J de Bruijn P Sparreboom A Journal of Clinical Oncology 2002 20 14 3122 3129 12118026 OWL Web Ontology Language Reference http://www.w3.org/tr/owl-ref/ Biosynthesis of an aminopiperidino metabolite of irinotecan [7-ethyl-10-[4-(1-piperidino)-1-piperidino]carbonyloxycamptothecine] by human hepatic microsomes Haaz MC Riche C Rivory LP Robert J Drug Metabolism and Disposition 1998 26 8 769 774 Metabolism of Irinotecan (CPT-11) by Human Hepatic Microsomes: Participation of Cytochrome P-450 3A and Drug Interactions Haaz MC Rivory L Rich C Vernillet L Robert J Cancer Research 1998 58 3 468 472 9458091 Bioactivation of the Anticancer Agent CPT-11 to SN-38 by Human Hepatic Microsomal Carboxylesterases and the in Vitro Assessment of Potential Drug Interactions Slatter JG Su P Sams JP Schaaf LJ Wienkers LC Drug Metab Dispos 1997 25 10 1157 1164 9321519 Common human UGT1A polymorphisms and the altered metabolism of irinotecan active metabolite 7-ethyl-10-hydroxycamptothecin (SN-38) Gagne JF Montminy V Belanger P Journault K Gaucher G Guillemette C Molecular pharmacology 2002 62 3 608 617 12181437 Stochastic prediction of CYP3A-mediated inhibition of midazolam clearance by ketoconazole Chien JY Lucksiri A Charles S Ernest I Gorski JC Wrighton SA Hall SD Drug Metab Dispos 2006 34 7 1208 1219 16611859 Development of a physiology-based whole-body population model for assessing the influence of individual variability on the pharmacokinetics of drugs Willmann S Hohn K Edginton A Sevestre M Solodenko J Weiss W Lippert J Schmitt W J Pharmacokinet Pharmacodyn 2007 34 3 401 431 17431751 The knowledge model of Protege-2000: Combining interoperability and exibility Noy N Fergerson R Musen M Proc EKAW 2000 2000 The Development of Camptothecin Analogs in Childhood Cancers Bomgaars L Berg SL Blaney SM The Oncologist 2001 6 506 516 11743213 Clinical Pharmacokinetics of Ketoconazole Daneshmend TK Warnock DW Clinical Pharmacokinetics 1988 14 13 34 3280211 N-Deacetyl ketoconazole-induced hepatotoxicity in a primary culture system of rat hepatocytes Rodriguez RJ Acosta DJ Toxicology 1997 117 123 131 9057891 Flavin-containing Monooxygenase-mediated Metabolism of N-Deacetyl Ketoconazole by Rat Hepatic Microsomes Rodriguez RJ Proteau PJ Marquez BL Hetherington CL Buckholz CJ O'connell KL Drug Metab Dispos 1999 27 8 880 886 10421614 Pharmacogenomics of Human UDP-Glucuronosyltransferases and Irinotecan Toxicity Tukey RH Strassburg CP Mackenzie PI Molecular Pharmacology 2002 62 3 446 450 12181419 Current Perspectives on the Clinical Experience, Pharmacology, and Continued Development of the Camptothecins Garcia-Carbonero R Supko JG Clinical Cancer Research 2002 8 641 661 11895891 Role of organic anion transporter OATP1B1 (OATP-C) in hepatic uptake of irinotecan and its active metabolite, 7-ethyl-10-hydroxycamptothecin: in vitro evidence and effect of single nucleotide polymorphisms Nozawa T Minami H Sugiura S Tsuji A Tamai I Drug metabolism and disposition: the biological fate of chemicals 2004 33 3 434 439 15608127 The Unified Medical Language System (UMLS): integrating biomedical terminology Bodenreider O Nucleic Acids Research 2004 32 D267 D270 308795 14681409 A reference ontology for biomedical informatics: the Foundational Model of Anatomy Rosse C Mejino JLV Journal of Biomedical Informatics 2003 36 478 500 14759820 The Foundational Model of Anatomy in OWL: Experience and Perspectives Golbreich C Zhang S Bodenreider O Journal of Web Semantics 2006 4 181 195 2270940 18360535 Relations in biomedical ontologies Smith B Ceusters W Klagges B Kohler J Kumar A Lomax J Mungall C Neuhaus F Rector AL Rosse C Genome Biology 2005 6 R46 1175958 15892874 Wielemaker J SWI-Prolog 5.6 reference Manual, Updated for version 5.6.32. Amsterdam 2007 An open graph visualization system and its applications to software engineering Gansner ER North SC Software – Practice and Experience 2000 30 11 1203 1233 In vitro binding and partitioning of irinotecan (CPT-11) and its metabolite, SN-38, in human blood Combes O Barre J Duche JC Vernillet L Archimbaud Y Marietta MP Tillement JP Urien S Investigational New Drugs 2000 18 1 5 10830136