Abstract
A mean residence time (MRT) is an important pharmacokinetic parameter. To the author's knowledge, however, a physiologically based structure of MRT (thereafter MRT structure) has not been published so far. Primarily this is because MRT structures cannot be identified by traditional pharmacokinetic methods used for the determination of MRT. Therefore, tools from the theory of linear dynamic systems were used for the structural identification of MRT in this study. The MRT structure identified is physiologically meaningful. Accordingly, it seems that the MRT structure identified may contribute to already established knowledge about MRT.
1. Introduction
A mean residence time is an important pharmacokinetic parameter. However, to the author’s knowledge, a physiologically based structure of (thereafter structure) has not been published so far. Primarily this is because structures cannot be identified by traditional pharmacokinetic methods, commonly used to determine . Therefore, tools from the linear dynamic system theory were used for the structural identification of in this study [1].
Drug disposition is a composite procedure arising from the complexity and diversity of interactions between the drug administered, physiological mechanisms, and various exogenous factors. Furthermore, drug disposition is controlled to a greater extent by several dynamic pharmacokinetic processes [2]. For that reason, several studies described investigations of drug disposition with the aid of dynamic systems, using the following steps: (1) a construction of an -related dynamic system, (2) a development of a mathematical model of the system constructed, for example, [3–7]. -related dynamic systems are mathematical constructs (working tools) without any pharmacokinetic relevance. The meaning of the acronym ADME is explained in many pharmacokinetic studies, for example, [8]. In the present study, the ADME-related dynamic system was simply called the system.
There are the following highly significant differences in terminology between pharmacokinetics and this study, which may trigger misunderstanding for readers: The difference between the physiological nature of information conveyed by a physiological system and the functional nature of information conveyed by the system used in this study. The difference in the use of the term “dynamic”. In pharmacokinetics, the term “dynamic” is commonly used in expressions describing drug actions. In this study, the term “dynamic” was used to indicate that the system or process changes over time.
2. Methods
The structural identification of was performed employing these steps.
(1) The introduction of the following approximate assumptions: an instantaneous mixing of the drug with blood, a uniform drug distribution [9]. A linear drug disposition [9–11]. The liver is the only eliminating organ of significance for the drug administered [12].
(2) The use of the method described previously [4, 7] and a theoretical example in which it was assumed that the drug was administered in an intravenous bolus dose to a hypothetical subject.
(3) The formalization of drug disposition using the system , created with the following configuration: the drug administration was considered as an input to the system , and the concentration-time profile of the drug in arterial blood was considered as an output of the system .
(4) The development of a circulatory model of the system , taking into account the fact that drug disposition can be regarded as a result of repetitive passes of the drug around the blood circulation, for example, [10, 13–16].
(5) The determination of the transfer function of the system , for example, [3–7].
(6) The derivation of the general equation (1) for the determination of : using the transfer function ( is the Laplace variable) and all assumptions made until now, for example, [3–7, 11].
(7) The use of the general equation (1), the circulatory model developed, and the method described previously [4] to identify the structure.
3. Results and Discussion
The developed circulatory model of the system is depicted in Figure 1. As seen, major body organs are lumped into subsystems of the system , that is, into somewhat independent parts of the system , for example, [13–17]. The model takes into account the fact that drug disposition can be regarded as a result of repetitive passes of the drug around the blood circulation, for example, [10].
The identification of the structure revealed the presence of five structural components of , which were denoted by . The structural components relate directly to the drug transport to the blood circulation through the following subsystems: the cardiopulmonary subsystem [18], the corresponding structural component is , the portal-venous subsystem [19], the corresponding structural component is , the hepatic-portal subsystem [20], the corresponding structural component is , the subsystem describing drug disposition in noneliminating tissues [21], the corresponding structural component is , the subsystem , if the drug is subject to the enterohepatic cycling , for example, [22, 23], the corresponding structural component is .
If the drug is not subject to the , the structural component can be described as where is the blood flow in the subsystem , is the mean time of the drug transport through the subsystem , and is the hepatic clearance. The structural components , , and can be described as In equations above, is the blood flow in the portal vein, is the blood flow in noneliminating tissues, is the mean time of the drug transport through the subsystem , is the mean time of the drug transport through the subsystem , and is the mean time of the drug transport through the subsystem , where is the mean time of the drug transport through a noneliminating tissue; the subscript specifies the tissue [24, 25]. If the drug is subject to the , the structural component can be described as where is the mean time of the drug transport through the subsystem , and the coefficient , determines the fraction of the drug that is the subject to the .
The resulting equation (6) describes the structure identified From the text above it evident that the right-hand side of the resulting equation (6) is the sum of (2)–(5). Equation (6) looks mathematically elegant and very simple. Nevertheless, this equation provides a mathematical description of the physiologically based structure of the mean residence time of the drug administered as the intravenous bolus to the subject (here to the hypothetical subject, as specified above).
The transfer function used in the general equation (1) is the mathematical relationship between the output and input of the system . Generally, transfer functions are characteristic functions of linear dynamic systems, providing complete descriptions of linear dynamic systems in the Laplace domain, for example, [3–5].
It is well known that, after an intravenous administration, a total amount of a drug is fully and immediately available to the blood circulation for transports to all areas of the body, a drug is not destroyed by digestive enzymes, an intravenous administration offers an advantage over other routes of administration in its accuracy. The purpose of recalling these well-known facts is to explain why the assumption of the intravenous drug administration was used in the theoretical example in this study.
The circulatory model developed possesses properties of simplicity and accuracy to describe the disposition of the drug administered as the intravenous bolus to the subject. The properties of the model make the model very flexible because the model is highly capable of accurately describing drug disposition in both situations, that is, when the drug is subject to the and also when the drug is not subject to the . This indicated that the model is very suitable for identifying the structure. The model is very general and appears applicable to several drugs. Using the model, the equations (2)–(6) were determined. Equations (2)–(6) can contribute to understanding mechanisms that control , they can be used to refine already established knowledge of , and can help to gain further insights into physiological background of . The equations (2)–(6) are easy to handle and may be sufficient in pharmacokinetic studies. Primarily this is because these equations are based on the commonly available physiological and pharmacokinetic knowledge.
The structure identified in this study has not been experimentally validated up to now. Its validity can be verified by further investigations, mainly experimental investigations; consequently a full pharmacokinetic exploitation of the structure identified lies far in the future.
For the sake of conciseness, mathematical details were restricted to a bare minimum. Differences between traditional pharmacokinetic approaches to and the approach presented in this study were left unexplained. Instead of a comparison of the approaches, this study gives rise to a reasonable expectation that the structure identified may be useful for basic research in pharmacokinetics. This is because the structure identified may reveal features of not apparent from values routinely determined by traditional methods which integrate influences of several processes in the body on into single numerical quantities, not providing any information about physiological backgrounds of .
4. Conclusion
This study presented a new view on “old” principles associated with . It attempted to contribute to the current understanding of , without an intention to criticize traditional approaches to . According to the best of the author’s knowledge, and after a Medline search, it can be stated that a physiologically based structure of has not been described in the literature as yet.
Disclosure
This work is related to the activities of the author in the COST Action BM0701 Antibiotic Transport and Efflux: New strategies to combat bacterial resistance, and in the Network of Excellence Virtual Physiological Human (FP7). The most comprehensive information about the Network of Excellence, Virtual Physiological Human can be found on the website: http://www.vph-noe.eu/vph-noe-membership.
Conflict of Interests
There is no conflict of interests.
Acknowledgments
This work was partially supported by the European Union through the Network of Excellence Biosimulation: a new tool in drug development (BioSim), Contract no. LSHB-CT-2004-005137 of FP6, and by the Slovak Academy of Sciences. The author thanks L. Dedík for his contribution to this study.