Next Article in Journal
T-Cell Metabolic Reprogramming in Atherosclerosis
Previous Article in Journal
Characterization of the Secretome from Spheroids of Adipose-Derived Stem Cells (SASCs) and Its Potential for Tissue Regeneration
Previous Article in Special Issue
Population Pharmacokinetic Model of Adalimumab Based on Prior Information Using Real World Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Upholding or Breaking the Law of Superposition in Pharmacokinetics

1
Faculty of Pharmacy and Pharmaceutical Sciences, University of Alberta, Edmonton, AB T6G 2E1, Canada
2
Escuela de Posgrado, Universidad Internacional Iberoamericana, Campeche 24560, Mexico
*
Author to whom correspondence should be addressed.
Biomedicines 2024, 12(8), 1843; https://doi.org/10.3390/biomedicines12081843
Submission received: 25 June 2024 / Revised: 23 July 2024 / Accepted: 7 August 2024 / Published: 13 August 2024
(This article belongs to the Special Issue Pharmacokinetics and Pharmacodynamics of Therapeutic Biologics)

Abstract

:
The law of superposition underpins first-order linear pharmacokinetic relationships. Most drugs, therefore, after a single dose can be described by first-order or linear processes, which can be superposed to understand multiple-dose regimen behavior. However, there are a number of situations where drugs could display behaviors after multiple dosing that leads to capacity-limited or saturation non-linear kinetics and the law of superposition is overruled. This review presents a practical guide to understand the equations and calculations for single and multiple-dosing regimens after intravenous and oral administration. It also provides the pharmaceutical basis for saturation in ADME processes and the consequent changes in the area under the concentration–time curve, which represents drug exposure that can lead to the modulation of efficacy and/or toxic effects. The pharmacokineticist must implicitly understand the principles of superposition, which are a central tenet of drug behavior and disposition during drug development.

1. Introduction

Although there are some drugs that can be taken acutely as a single dose, such as anti-migraine drugs [1,2,3,4] or hypnotics [5,6,7,8] for a limited period of time, the vast majority of drugs, such as antibiotics [9,10,11,12], anti-diabetics [13,14,15], and blood pressure medication [16,17,18], are designed to be repeatedly administered orally. Chronic disease states often require a multiple-dose regimen to enable the accumulation of the drug and attain steady-state concentrations [19,20,21]. The law of superposition or sometimes referred to as the principle of superposition enables the understanding of the pharmacokinetics of a drug after a single dose and for the accurate prediction of the steady-state conditions to be made for repeated dosing or, potentially, vice versa [22,23]. While some scientists, particularly in the pharmaceutical industry, may have training during pharmacy education on pharmacokinetics, others have a varied biomedical background and are less versed in the underlying pharmacokinetic theorems. In 21st century drug development, there are a myriad of pharmacokinetic technology programs and simulation software that will enable investigators to take and input single-dose pharmacokinetic concentration–time data and predict steady-state concentrations after a specific dose regimen. Modeling programs offer the advantages of large data handling capacities and execution speed; however, the conceptual and mathematical limitations to this modeling and the interpretation of the results must be clearly understood and delineated.
The law of superposition is a fundamental aspect of pharmacokinetics, and it is often dealt with in a cursory manner in textbooks. Few articles have tried to deal with a conceptual understanding of superposition in a practical pharmacokinetic manner to demonstrate the application of the principle [19,22,23,24,25,26,27,28,29,30,31,32]. In this article, we first present a review of the principles of superposition and then discuss its application to multiple dosing. Finally, we apply this principle through simulations to gain an understanding of where the law of superposition does not hold and illustrate the non-linearity in the processes that can, in fact, overturn this law.

Law of Superposition

Pharmacokinetically, the body can be considered to be a construct function that in mathematical terms converts the input (dose rate) into an output or response (plasma concentration of drug), as a function of time [32,33,34]. A linear system is characterized by input–output relationships that are completely defined by one or more linear equations, exhibiting properties of homogeneity and additivity. Linear systems applied to the field of pharmacokinetics were explored by Cutler [24,25], where the differential equations involved define a function. The most important feature of a linear system is that it is additive and obeys the law of superposition, where the output of a linear system equals the sum of the outputs of separately applied inputs [23]. The application of a linear model is contingent upon the condition that the concentrations involved are within a range where saturation effects are absent, indicating a limited scope of inputs [34,35]. For a specific input range, the linearity of the system can be experimentally tested to either reject or disprove the null hypothesis.
Simply put, the law of superposition is a mathematical concept that enables pharmaceutical scientists to predict concentration–time data [36,37]. It has, however, an underlying assumption that the preceding dose of a drug and its pharmacokinetics do not alter or affect in any way subsequent doses of a drug and the expressed pharmacokinetic profile [38,39,40]. In practice, the plasma concentration–time data after the first, second, or nth dose will be superimposed by the plasma concentration reached after the (n − 1)th dose [39]. For practical purposes, the area under the plasma concentration–time curve after a single dose is equal to the area under the plasma concentration–time curve after multiple doses when steady-state concentrations are attained [39] (Figure 1).
The law of superposition allows the pharmacokineticist to project the plasma drug concentration–time curve of a drug after multiple consecutive doses, based on the plasma drug concentration–time curve obtained after a single dose [39,40]. It is important in a clinical setting to have the patient exposed to the target drug concentration without excessive fluctuation and drug accumulation outside of the therapeutic window [37]. Following a single dose, the extent of overall drug exposure over time can be quantified by estimating the area under the concentration–time curve from zero to infinity (AUC0−∞). During multiple dosing, the area under the concentration–time curve over a steady-state dosing interval (AUC0−τ) serves as a measure of overall drug exposure. When the clearance (CL), the aggregate of all the elimination processes, remains constant (i.e., it does not change with increasing drug concentrations), the AUC0−τ at the steady state will be equivalent to the AUC0−∞ following a single-dose administration [39,41].
In the analysis of AUC data in pharmacokinetics, the law of superposition provides a simple test for detecting the presence of non-linear processes that do obey the law, i.e., processes that are first-order processes [34]. Practically, a greater than or less than response, in terms of the area under the drug concentration–time curve (AUC) with repeated doses is the expressed result of the involvement of zero-order processes [39]. AUC equivalence presumes that the dose administered as a single dose matches the dose given in a multiple dose. This is useful for determining pharmacokinetic parameters [23,27,42]. For example, after a single IV bolus dose, the total body clearance (CL) can be computed using Equation (1).
C L = D o s e A U C 0
AUC equivalence enables the estimation of clearance under steady-state conditions using the same equation, but using AUC0−τ instead of AUC0−∞. The latter clearance estimate is often referred to as the steady-state clearance (CLss). In the case of a single oral dose, the AUC is calculated taking into consideration the bioavailability (F) of the drug when using CLtotal (Equation (2)) or taking into consideration the fraction unbound in blood (fu) when using the intrinsic clearance (CLint) [22] (Equation (3)).
A U C = F × D o s e C L t o t a l
A U C = F × D o s e f u × C L i n t
It becomes readily apparent in Figure 2 that in multiple dosing the concentrations oscillate at steady state between a minimum concentration at a steady state (CSSmin) and a maximum concentration at a steady state (CSSmax) and that an average steady-state concentration Cav,SS can be calculated [19].
The drug clearance can be determined by dividing the dose by either of the aforementioned two AUC values, or by dividing the dosing rate by the average steady-state concentration (Cav,ss). The Cav,ss is calculated by dividing either of the two AUC values by the dosing interval (Tau, T) [27]. Thus, the average steady-state plasma level (Cp,ss) can by calculated by Equation (4).
C p , s s = A U C 0 T
Superposition is based on mathematical equations describing linear processes [34,35]. Both experimental data and drug development results can uphold or refute the mathematical relationship [43]. Pharmaceutical development that relies solely on equations or theoretical superposition can be misleading or incorrect. There are observational limits to when superposition occurs and when it may deviate from linearity. Ab initio superposition involves multiple untested assumptions that might not hold for a particular drug in a particular disease state, as discussed in further detail below. Modeling approaches required for predictive analytics, including superposition, can utilize simple inductive reasoning and are only sometimes correlative with observational findings [30]. Further, a causality-based understanding of drugs is achieved through deductive reasoning and the application of a scientific method.
The law of superposition would predict that for a drug, the AUC after a single dose is equal to the AUC at the steady state. Therefore, Equation (5) would hold true.
A U C S i n g l e D o s e A U C S S = 1
This is independent of a compartmental model and is simply a function of linearity. Furthermore, the CL (after IV administration) or CL/F (after oral administration) following a single dose, will equal the clearance values at the steady state. In the case of oral administration, the bioavailability remains constant between a single dose and the steady state. The AUC at the steady state (AUCSS) can be defined as in Equation (6).
A U C S S = lim n t n t n + 1 C n t d t
Non-linearity can cause the AUCSS to differ from the AUCSingle Dose [32,43], this will be discussed later in the section on breaking the law of superposition. Most drugs will exhibit some deviation from the ideal accumulation, as described [37]. Nevertheless, comprehending the principle of superposition permits reasonable predictions of the pharmacokinetic behavior of drugs under repeat-dose conditions for a broad range of compounds. This understanding proves especially beneficial when transitioning from single-dose studies to repeat-dose studies during the initial stages of drug development. Additionally, this knowledge facilitates the design of repeat-dose regimens that achieve the desired drug concentrations efficiently and reliably within a clinically acceptable timeframe. When determining multiple-dose regimens, it is crucial to assess whether successive doses will impact the effects of previous doses. The principle of superposition assumes that earlier doses do not influence the pharmacokinetics of later doses. As a result, the blood levels following the second, third, or nth dose will overlay those achieved after the (n − 1)th dose [39]. In addition, the AUC0−∞ following single-dose administration is equal to the steady-state area between doses (Figure 1). The disposition processes after an oral dose are absorption, distribution, metabolism, and excretion (ADME). Therefore, the superposition principle is applicable when the ADME processes exhibit linear or first-order kinetics. With this approach, the concentrations after multiple doses can be determined by summing up the concentrations from each dose [24,25,42]. Furthermore, doubling the dose will result in proportional doubling of the concentrations for each time. However, this principle does not apply when any of the disposition processes are non-linear [24,25,34].
Mathematically a single-dose concentration–time curve C(t) is denoted as a function over time f(t) [42]. Multiple dosing occurs when the administered doses remain constant over time, where the drug is administered at fixed time intervals (tau) (Figure 3). The concentration at time t after n doses is equal to the sum of the concentration at t after n − 1 doses and the single-dose concentration at t − tn (from the last dose) [22]. Therefore, the contribution of the last dose can be superimposed to that of the first n − 1 doses [22]. The concentration–time curve Cn(t) is described in Equations (7) and (8).
C n t = C n 1 t + f ( t t n )
C n t = i = 1 n f ( t t i )

2. Modeling a Hypothetical Drug “Canadamycin”

2.1. Pharmacokinetics after a Single Intravenous Dose

For practical purposes, we named our hypothetical drug Canadamycin, which was to be administered via the intravenous (IV) route at a dose of 100 mg. It has a volume of distribution of 5 L, a half-life of 4 h, and dosing interval (tau) of 8 h. Therefore, a typical three times daily dosing (TID) is required. The C0 was estimated to be 20 mg/L (µg/mL), after dividing the dose by the volume of distribution (100 mg/5 L).
A direct intravenous bolus injection into a vein results in the near instantaneous achievement of C0, followed by a first-order reduction in the concentration. This route is frequently employed for drugs administered in hospital clinics or emergency settings, where a quick onset of action is critical. Pharmacokinetics textbooks typically utilize this approach to elucidate key concepts associated with drug accumulation [37,39]. Regardless of the administration route, as the number of doses progressively increases during a therapeutic drug regimen, drug concentrations within the dosing intervals (tau) will eventually stabilize, forming a plateau [42]. This phenomenon consistently occurs under the law of superposition, when conditions exhibit constant total clearance. At this plateau, the rate of drug administration matches the rate of drug elimination. When this equilibrium is achieved, plasma concentrations are considered to be at a steady state [42].
After a single IV bolus dose administration, the concentration of the drug will begin to decline and will be removed (cleared) and eliminated by the body. The concentration will decrease over time in terms of the half-life (the time required for the concentration to decline by 50%). Canadamycin has a half-life of 4 h, meaning that 8 h post-dose we would observe a concentration of 5 mg/L, a concentration of 1.25 mg/L at 16 h post-dose, and a concentration of 0.3124 mg/L at 24 h post-dose, etc. (Table 1). If no further doses are administered, the concentration will continue to decline by an additional 50% every 4 h. In this scenario, each subsequent dose is given before the previous dose is fully cleared from the body. As a result, the concentration at the beginning of the next dose is the sum of the drug amount remaining in the body (Cmin, denoted as Cτ above), plus the Cmax from a single IV bolus dose (or C0). This cumulative effect of drug concentrations from previous doses and additional doses is the principle of superposition, which implies that each dose will additively influence the overall concentration.
Drug accumulation has a limit because, as the plasma concentration increases, the amount of drug eliminated (1 − r) during the dosing interval also grows. This is due to the fact that the elimination rate is directly proportional to the drug amount in the body, multiplied by the rate constant for first-order elimination [21,42]. As seen above, accumulation will occur if the dosing interval is less than 6 half-lives. Accumulation is inversely proportional to the fraction of the dose lost during an interval. The accumulation factor is equal to 1 over the fraction of the dose lost during a dosing interval = 1/(1 − fraction remaining) [21,42]. If a drug is given once every half-life, in the case of Canadamycin every 4 h, the accumulation of the dose is 2 [1/(1 − 0.5)]. If Canadamycin is administered every 8 h (or 2 half-lives), the accumulation factor is 1.33 [1/(1 − 0.25)], since there would be 25% of the drug remaining.
When a drug is administered at a fixed dose and a regular interval, as in multiple-dose regimens, the amount of drug in the body will increase and then plateau at a mean plasma level higher than the peak concentration (C0). If the second dose is administered before the previous dose is entirely eliminated, drug accumulation will occur. This accumulation is quantified by the R index and is influenced by the elimination constant and the dosing interval, but it is independent of the dose itself [21,42]. For a drug administered in repeated oral doses, the time required to reach the steady state depends on the elimination half-life of the drug. This time is independent of the size of the dose, the length of the dosing interval, and the number of doses [39]. For instance, if the dose or dosing interval of a drug is modified, the time needed to reach the steady state remains unchanged, but the final steady-state plasma level adjusts proportionately. Consequently, the blood concentration after the second, third, or nth dose will superimpose the blood concentrations attained after the (n − 1)th dose.
The principle of superposition allows for the prediction of the plasma drug concentration–time profile after multiple consecutive doses, by utilizing the plasma concentration–time data from a single dose. This method is based on two key assumptions: that the drug follows first-order elimination kinetics, and that its pharmacokinetics remain unchanged with repeated dosing. As a result, the plasma concentration after multiple doses can be forecasted from those observed after an initial dose. The first-order pharmacokinetics can apply to multiple dosing; consider Equation (9) to describe the plasma concentration at any time after a single IV bolus drug administration dose that follows linear pharmacokinetics.
C t = 1 n A i e t / τ j = 1 m A i e a j t
where Ai is the concentration coefficient, τj is the time constant, and aj is the rate constant (aj = 1/τj).
The concentration at the end of the first dosing interval (tau) is described in Equation (10). Equation (11) describes the concentration at the start of the second interval, while Equation (12) describes the concentration at the end of the second dose interval.
C p 1 τ = C p 1 0 e k e l t
C p 2 0 = C p 1 τ + C p 1 0 = C p 1 0 e k e l t + C p 1 0
C p 2 τ = C p 1 0 e k e l t + C p 1 0 e k e l t
As described in Table 1, there is a fraction of the initial plasma concentration remaining at the end of the dosing interval. Intuitively there is both a fraction lost (R) and a fraction remaining in the body (1 − R) for every dosing interval. The remaining fraction in the body before the next dose depends only on the first-order rate constant of elimination and the time after the dose, as described in Equation (13).
R = e k e l t
Thus, replacing R in Equation (12) to describe the concentration at the end of the second dose interval is shown in Equation (14). Equation (15) describes the concentration at the start of the third dose interval.
C p 2 τ = C p 1 0 R + C p 1 0 R 2
C p 2 0 = C p 1 0 + C p 1 0 R + C p 1 0 R 2
These fractions of the drug that remain in the body before the next dose mathematically form a geometric series of remaining concentrations, with each term R times the preceding term. This leads to drug accumulation. In the case of multiple dosing, it appears that the concentration profiles of the individual doses can be totaled or summed up over time (Figure 4). The steady state can then be calculated by the addition of areas, without the need to use any pharmacokinetic model, provided that first-order linearity holds [34].

2.2. Pharmacokinetics after Multiple Intravenous Dosing

In linear pharmacokinetic systems, the drug concentration achieved by a second dose is comparable to that produced by the first dose. For instance, when a second 100 mg IV dose of Canadamycin is administered at the 8-h mark, 5 mg/L of the drug remains from the first dose, as shown in Table 1. Consequently, the plasma concentration (Cp) immediately after the second dose would be 25 mg/L, consisting of the 5 mg/L remaining from the initial dose and an additional 20 mg/L from the new dose. This results in a combined concentration of 25 mg/L. In linear systems, each subsequent 100 mg dose will generate a concentration similar to that of the first dose. However, the resulting concentration will be the sum of the residual concentrations from all prior doses and the concentration from the most recent dose. For example, administering a third dose at the 24-h mark will also produce an initial concentration of 20 mg/L. The concentration remaining from the second dose would be 15 mg/L, and the concentration remaining from the first dose would be 2.5 mg/L. When these concentrations are combined, the observed total concentration will be approximately 26.25 mg/L, calculated as 20 mg/L (from the third dose) + 3.75 mg/L (remaining from the second dose) + 1.25 mg/L (remaining from the first dose). This additive process of combining concentrations from multiple doses to determine the observed concentration exemplifies the principle of superposition [21,42]. The described pattern of superposition assumes that each dose behaves similarly, despite an increasing concentration. While simple additive superposition is generally applicable to most drugs, it is important to note that the superposition of exposures can become more complex if the pharmacokinetic behavior of each dose changes with an increasing concentration. This complexity arises in drugs with saturable clearance, such as phenytoin, where the pharmacokinetics do not remain consistent as the concentration increases [44]. For the remainder of this review, the simple additive superposition scenario will be used.
Consecutive multiple doses of a drug will lead to an increasing concentration in the body until a plateau, known as the steady state, is reached. At a steady state, the amount of drug administered with each dose is balanced by an equivalent amount of drug eliminated from the body between doses, resulting in a constant rate of drug input and output (rate in = rate out) [39]. For the majority of drugs, reaching a steady state typically takes about five half-lives. The time required to attain a steady state in a repeat-dose regimen is governed exclusively by the half-life of the drug and the dosing interval, as accumulation and the remaining fraction are influenced by specific pharmacokinetic equations. Although changing the dose or dosing interval can influence the concentration at the steady state, it does not affect the time needed to reach the steady state. Some drugs have prolonged half-lives, ranging from days to weeks or even longer. In such cases, waiting for five half-lives to achieve the desired steady-state concentration may not be necessary. When time is critical, as in the administration of antibiotics for critically ill patients, a loading dose can be used to reach the steady state more quickly [39].
For our hypothetical drug, Canadamycin, the Cmax and Cmin values using Equations (16) and (17), respectively, can be estimated taking into consideration the following calculations:
C m a x = C p 1 0 ( 1 R ) = 20 ( 1 0.257 ) = 26.92   m g / L
C m i n = C p 1 R 0 ( 1 R ) = C max R = 26.92 × 0.257 = 6.92   m g / L
Dose = 100 mg
Vd = 5 L
t½ = 4 h
τ = 8 h
C0 = 20 mg/L
kel = 0.693/t½ = 0.693/4 h = 0.17 h−1
R = e k e l t = e 0.17 × 8 = 0.257
Therefore, the plasma concentration will fluctuate between 26.92 mg/L and 6.92 mg/L during each dosing interval when the steady state is reached.

2.3. Pharmacokinetics after Multiple Oral Dosing

For oral dosing, the equations are more complex, but it is possible to model and predict the concentration with repeated doses and to estimate the plasma concentration and dosing regimens using superposition. The Cmax value corresponds to the concentration at the peak time (tpeak) after reaching the steady state. In this case, Cmin is estimated based on the assumption that the plasma concentration has peaked and e−ka t is close to zero and that the next dose is administered once the absorption phase is complete, which is described in Equation (18) and a simplified version is represented by Equation (19).
C m i n = F × D o s e × k a V × ( k a k e l ) e k e l τ 1 e k e l τ
C m i n = A R 1 R
The relationship between the loading dose and maintenance dose and, thus, drug accumulation during multiple-dose administration can be examined by analyzing the ratio between the minimum concentration at the steady state and the concentration during one dosing interval (τ) after the first dose, as demonstrated in Equation (20).
C m i n C p 1 τ = 1 1 e k e l τ = 1 ( 1 R )
Equations have been developed that facilitate the calculation of the plasma concentration achieved following multiple oral administrations. To start, the plasma concentration achieved following a single oral dose that does not have flip-flop kinetics [45] can be described by Equation (21).
C p = F × D o s e × k a V × ( k a k e l ) e k e l τ e k a τ
An intriguing outcome of this equation is that the average plasma concentration remains the same whether the dose is given as a single dose every dosing interval (τ) or divided into shorter dosing intervals.
For multiple oral dosing, we will use our hypothetical drug, Canadamycin, with an oral dose of 300 mg every 12 h, a bioavailability (F) of 100%, a volume of distribution of 30 L, and a half-life of 6 h.
Dose = 300 mg
Vd = 30 L
t½ = 6 h
τ = 12 h
kel = 0.693/t½ = 0.693/6 h = 0.116 h−1
We will estimate the required dose (maintenance dose) to be administered every 12 h to achieve an average plasma concentration of 15 mg/L using Equations (22) and (23).
C p = F × D o s e V × k e l × τ
D o s e = C p × V × k e l × τ F = 15 × 30 × 0.116 × 12 1 = 624   m g
The loading dose can be calculated using Equation (24), with the use of the R value.
R = e k e l t = e 0.116 × 12 = 0.25
L o a d i n g   D o s e = M a i n t e n a n c e   D o s e 1 R = 624 1 0.25 = 832   m g
To estimate the fluctuations in the plasma concentration, the Cmin can be calculated using Equation (25).
C m i n = F × D o s e V e k e l τ 1 e k e l τ = 1 × 624 30 0.25 1 0.25 = 6.93   m g / m L
Therefore, the plasma concentration would likely fluctuate between approximately 7 and 23 mg/L, with an average concentration of around 15 mg/L. In this case, the average concentration can be estimated using an approximation of Cmax = [Cav + (Cav − Cmin)] or [23 = 15 + (15 − 7)].
Alternatively, half the maintenance dose (312 mg) could be administered more frequently (half the original tau), every 6 h, as calculated using Equation (26).
C m i n = F × D o s e V e k e l τ 1 e k e l τ = 1 × 312 30 0.5 1 0.5 = 10.4   m g / m L
Thus, the plasma concentration would fluctuate between about 10.4 to 20 with an average of 15 mg/L.

3. Breaking the Law of Superposition

The law of superposition applies when all the disposition processes are linear or first order, this includes absorption, distribution, metabolism, and elimination (ADME) [23,32,34,42,43]. Thus, the concentration after multiple doses can be calculated by adding together the concentrations from each dose. The assumptions underlying the law of superposition are as follows:
  • The dosing is in a range where the pharmacokinetics are linear and dose proportional;
  • The rate and extent of absorption and average clearance are the same for each dosing interval;
  • Each dose acts independently from one another, and the sum of all these dosing events provides the total concentration of the drug in circulation.
There are many situations, however, in which the superposition principle does not apply. In these cases, the pharmacokinetics of the drug change after multiple dosing due to various factors, including the changing pathophysiology of the patient, the saturation of the drug carrier system, enzyme induction, and enzyme inhibition [43]. Drugs that follow non-linear pharmacokinetics generally do not have predictable plasma drug concentrations after multiple doses, according to the superposition principle. When the law of superposition is broken, the prediction of the pharmacokinetics, efficacy, and toxicity can become more complicated [23,32,34,42,43]. As described above, the underlying theorem is that the AUC for a single dose will equal the AUC at the steady state. However, the AUC after an oral dose can be affected not only by absorption, but also by protein binding and elimination [43]. Therefore, the AUC at the steady state (AUCSS) can either be higher or lower than the AUC after a single dose (AUCSingle Dose). Figure 5 demonstrates the linearity of pharmacokinetics and the validity of superposition, as well as in instances where non-linearity in pharmacokinetics causes the breakdown of the superposition principle. An overview of possible causes related to various ADME processes, along with examples, is provided in Table 2.

3.1. Less than Proportional Increase in AUC (AUCSS < AUCSingle Dose)

We first consider saturable processes such as absorption, distribution, and elimination. There are at least three processes that can affect drug absorption, namely the saturate–dissolution rate, the GI transit time, and the ability or inability of drugs to cross intestinal barriers [43]. While most drugs are absorbed through passive diffusion, some utilize specific transport processes in the small intestine. When passive diffusion plays a minor role, a less than proportional increase in the area under the curve (AUC) may occur due to the saturation of the active transport process. Drug absorption can be described by a combination of Michaelis–Menten kinetics and first-order kinetics [23,32,34,42,43]. This situation does not always hold. It remains that the absorption of most drugs follows linear kinetics, and the pharmacokinetic parameters describing the absorption of a drug do not change over the therapeutic dose range. However, there are cases where transporters are involved in absorption and, when high doses are utilized (or in overdose situations), the linear range may be exceeded [43]. Sorafenib is an oral tyrosine kinase inhibitor that has been approved for treating advanced renal cell carcinoma and hepatocellular carcinoma. It has been found that the administration of sorafenib at doses exceeding 800 mg/day would result in a less than proportional increase in the area under the curve (AUC) of sorafenib with respect to repeated doses. This phenomenon was attributed to the saturation of the absorption of sorafenib at higher doses. This study has systematically considered and excluded all plausible contributing factors. It has also been demonstrated that partitioning of the daily dose into three discrete administrations would ameliorate the absorption saturation and lead to higher drug exposure at higher drug doses, compared with once or twice daily doses at which the absorption saturation was more notable [57].
Since the concentration of a drug after oral dosing depends on both absorption and elimination, the area under the curve (AUC) represents the combined outcome of these processes. Elimination (clearance) influences this relationship, as a less than proportional increase in AUC at the steady state may result from both increased elimination and decreased absorption. It is generally accepted that only the unbound drug can diffuse across membranes, enabling distribution between the vascular and tissue compartments. Accordingly, changes in drug–protein binding in plasma and tissues can impact the distribution of drugs within the body [32,34,43]. Drug distribution is often measured by the volume of distribution, usually defined as the ratio of the amount of drug in the body to the plasma concentration. The most basic quantitative expression relating to the volume Vd of distribution to binding in plasma and tissue is provided by Equation (27).
V d = V p + V t f u f t
where Vd is the volume of distribution in plasma, Vp is the plasma volume, Vt is the tissue volume, and fu and ft are the fraction of the unbound drug in the plasma and tissue, respectively.
It is clear that the volume of distribution depends on plasma protein binding, but also on tissue binding. For drugs with low clearance, protein binding significantly influences drug elimination, and total body clearance is directly proportional to the unbound fraction in plasma [23,34,39,43]. A drug can exhibit concentration-dependent protein binding, with the unbound fraction changing [43]. A less than proportional increase in the AUC due to the saturation of protein-binding sites in plasma may result in an increase in the free fraction, which will result in enhanced clearance. The AUCSS would be increased less than proportionally, and the bioavailability reduced; thus, the AUCSS would be less than the AUCSingle Dose. Assuming a drug binds to n discrete sites on a protein molecule, with all sites having the same affinity for the drug and acting independently, the fraction of the unbound drug, fu, can be described by Equation (28).
f u = K d + D f n p t + K d + D f
where Kd is the dissociation constant at the equilibrium, Df and pt are the concentration of the unbound drug and total protein (bound and unbound), respectively, and n is the number of binding sites on the protein molecule. The binding of a drug to a protein is determined by four factors: the drug concentration, the protein concentration, the number of binding sites per protein molecule, and the dissociation constant [43]. Irrespective of the affinity, the maximum possible binding capacity is the product of the molar concentration of the protein and the number of binding sites per protein molecule. Mathematically, it can be observed that with multiple dosing, an elevated unbound fraction in plasma, fu, becomes apparent in plasma. This results in an increased volume of distribution, lowered plasma concentrations, and a less than proportional increase in the AUC. Vismodegib (Erivedge™, GDC-0449; Genentech, South San Francisco, CA, USA) is an approved drug in the United States for the treatment of advanced basal cell carcinoma (BCC). An evaluation of the intravenous pharmacokinetic (PK) parameters, following both single and multiple oral dosing regimens, showed increased clearance and an increased volume of distribution, with values 81% and 63% higher, respectively, following repeated dosing. Moreover, a significant reduction in bioavailability (77% lower) was observed with multiple doses. Intriguingly, continuous daily dosing of vismodegib led to a 2.4-fold increase in the unbound fraction compared to a single dose [58].
In tissues, drugs can bind to non-specific binding proteins, enzymes, and receptors [43]. Enzymes and receptors are particularly important pharmacologically, as they contribute to the therapeutic effect. However, distinguishing these from non-specific binding sites is challenging because the binding capacity of enzymes and receptors are often low and overshadowed by non-specific binding. Similar to plasma protein binding, tissue protein binding can become saturated as drug concentrations increase over time, resulting in changes in the volume of distribution. However, the saturation of tissue binding and of plasma protein binding can occur simultaneously [43]. The volume of distribution is directly related to plasma protein binding and inversely related to tissue binding, so changes in one type of binding can counterbalance the other. This complexity is further compounded by the fact that different tissues may not become saturated to the same extent. Consequently, the volume of distribution of a drug may not be significantly affected by changes in tissue binding, depending on the degree of non-linear tissue binding. Thus, it is essential to evaluate the non-linearity of tissue binding by examining the ratio of the drug concentration in a given tissue to that in plasma. At equilibrium, the tissue/plasma concentration ratio should match the ratio of the unbound fraction in plasma to the unbound fraction in a non-eliminating tissue [43]. In eliminating tissues, it is more complex; the ratio may be influenced by blood flow to the tissue and the eliminating process. Moreover, elimination could also be impacted by changes in intrinsic clearance (CLint’). Enzyme autoinduction is a time-dependent phenomenon, where the elimination of a drug increases with multiple doses [43]. Thus, an AUCSS less than the AUCSingle Dose may be observed. One such example is the anticonvulsant carbamazepine (CBZ), that is known to auto-induce its own metabolism. Throughout the course of multiple 200 mg doses of CBZ, the daily trough CBZ levels consistently declined, reaching their lowest point after two weeks. Notably, the observed CBZ concentration was only 53.9% of what had been predicted under the assumption of no autoinduction. Furthermore, the total CBZ clearance following multiple doses was higher than after a single dose. Digging deeper into CBZ pharmacokinetics, considerable increases in the AUCs of CBZ metabolites during multiple-dose administration were noted. These changes were accompanied by increased urinary excretion of these metabolites compared to CBZ, all due to the phenomenon of autoinduction [59].
Phenytoin serves as another illustrative case, akin to carbamazepine. Like carbamazepine, phenytoin, an antiepileptic medication, also exhibits violation of the law of superposition due to enzyme autoinduction. In the case of phenytoin, after the third dose, the calculated AUC was found to be statistically significantly smaller than that observed after the initial dose. This phenomenon was attributed to the ability of phenytoin to induce its own metabolism. Consequently, this autoinduction led to an increased rate of elimination, which in turn resulted in a reduction in the AUC following multiple doses [60].

3.2. More than Proportional Increase in AUC (AUCSS > AUCSingle Dose)

A pharmacokinetic situation in which an increase in the AUC at the steady state is greater than after a single dose is also possible. In this case, a lack of proportionality is more likely due to a decrease in elimination than to an increase in absorption [43]. However, in certain instances, absorption can increase with higher doses, due to specific mechanisms.
The more than proportional increase in the AUC is likely due to non-linear clearance. This disproportionate increase in the AUC is most likely attributable to decreased elimination clearance after multiple doses [32,43]. For drugs with a high level of extraction in the liver, a significant first-pass effect occurs. If this first-pass effect diminishes with multiple dosing, it can lead to a disproportionate increase in the AUC. This increase is due to both a decrease in elimination and an increase in bioavailability, the latter resulting from the saturation of hepatic metabolism during the initial passage of the drug through the liver [43]. The disproportionate increase in the AUC with multiple oral doses can be linked to enhanced bioavailability resulting from a reduced first-pass effect. The effects of capacity-limited metabolism are more apparent under steady-state conditions than after a single dose. This is exemplified by the relationship between the plasma phenytoin concentration at the steady state and the administered dose. At the steady state, the plasma concentration of phenytoin increases more than proportionally in terms of the rate of the input [29,44]. Verapamil, a drug used to treat cardiovascular conditions, exhibits low oral bioavailability, primarily due to extensive pre-systemic (first pass) liver elimination. When the pharmacokinetics of verapamil were assessed, it was observed that the AUC at the steady state (1999 ± 435 ng/mL h) exceeded that recorded after the initial dose (788 ± 224 ng/mL h). That increase in AUC was attributed to the saturation of the first-pass effect that occurred with multiple dosing. As the first-pass effect became saturated, a greater proportion of the drug managed to evade pre-systemic elimination. That, in turn, led to the accumulation of verapamil and an increase in the AUC with multiple doses [61,62].
Drugs are eliminated from the body through various processes, including biotransformation and excretion. The efficiency with which the body eliminates a drug is commonly expressed as the total clearance, which is the sum of the individual clearances of different organs, occurring simultaneously. CLint is the intrinsic clearance, a measure of the intracellular removal of a drug, described by the Michaelis–Menten equation, where V is the capacity of the enzyme or carrier systems, K is the Michaelis constant, and C is the unbound drug concentration [39]. It is evident that the CLint remains constant when the concentration is lower than K, with a limited value of V/K. However, non-linearity arises when the drug concentration approaches the value of K. Since drug elimination typically involves multiple processes, determining whether elimination is linear or non-linear becomes increasingly complex as the number of these processes increases [43]. Depending on the relative contribution of saturable pathways to overall elimination, non-linearity may not be easily discernible.
When saturation occurs, the accumulation of the drug will be greater at the steady state than predicted from a single dose [21,42,43]. Drug elimination often exhibits zero-order kinetics at high concentrations and first-order kinetics at low concentrations, a phenomenon known as concentration-dependent kinetics. At multiple doses, which lead to a higher plasma concentration, zero-order kinetics are observed. Conversely, at lower doses, the kinetics are linear, or first order. This pattern is particularly common with drugs that are extensively metabolized [43]. A typical characteristic of enzymatic reactions and active transport is their limited capacity. The liver has a finite amount of enzymes, establishing a maximum rate at which metabolism can occur. Additionally, the rate of metabolism can be further constrained by the limited availability of co-substances or co-factors necessary for the enzymatic process [43].
Most of our understanding of enzyme kinetics comes from in vitro studies, where the concentrations of the substrates, enzymes, and co-factors are meticulously controlled. Numerous factors are involved in in vivo studies, making it difficult to isolate each one and assess it in detail. Nevertheless, the fundamental principles of enzyme kinetics are applicable to pharmacokinetics. At single dose concentrations, where Km > Cp, Km + Cp is approximately equal to Km, as described in Equation (29).
C p t = V m K m C p
where Vm/Km is a first-order elimination rate constant and the whole equation now looks like that for first-order elimination.
Therefore, at a low plasma concentration following a single dose, first-order kinetics would be expected. This is a typical situation for most drugs, as the Km is usually larger than the plasma concentration achieved. However, a high concentration at multiple dosing is achieved, that is Cp > Km, then Km + Cp is approximately equal to Cp leading to a zero-order elimination process, as described in Equation (30).
C p t = V m × C p C p = V m
At high steady-state plasma concentrations, zero-order or concentration-independent kinetics occur. The presence of saturation kinetics can be significant when multiple doses of certain drugs are given, or in cases of overdose. At the steady state, the effective elimination rate constant is reduced, leading to excessive drug accumulation if the saturation kinetics are not properly understood.
The renal excretion of drugs usually involves three processes: glomerular filtration, renal tubular secretion, and reabsorption from the renal tubular lumen [43]. Glomerular filtration is a passive process that depends on the unbound concentration of a drug in plasma. The relationship between renal clearance CLR and these processes is described in Equation (31).
CL R = f u × G F R + C L
where GFR, CL, and CLR are the glomerular filtration rate, secretion clearance, and reabsorption clearance, respectively, and fu is the unbound fraction of the drug in plasma. When the CLR and GFR of a drug are much greater than one, it indicates that renal tubular secretion of the drug is occurring [43]. Conversely, when the ratio is much less than one, it suggests reabsorption of the drug from the tubular lumen. Renal tubular secretion is a specialized and saturable process, whereas tubular reabsorption can occur through passive diffusion or active transport, while active reabsorption is also saturable [43]. In addition, there are processes whose rates are affected by pharmacological actions of the drug itself. Elimination could also be affected through changes in intrinsic clearance (CLint’). Enzyme inhibition can be a time-dependent phenomenon, in which the elimination clearance of a drug decreases following multiple doses and the decrease in clearance is significant compared to a single dose [43]. Thus, the AUCSS would be higher than the AUCSingle Dose. Propranolol and itraconazole provide examples of how changes in intrinsic clearance (CLint’) can affect elimination. When administered in multiple doses, these drugs have been found to inhibit their own metabolizing enzymes. That inhibition resulted in a decrease in their hepatic clearance and led to an increase in the AUC at the steady state [63,64].

3.3. Other Factors Affecting Adherence to the Law of Superposition

The reasons discussed in the previous two sections pertain to situations where the law of superposition may not accurately characterize drug behavior due to drug-related factors. However, in other instances, a variety of factors can play a role. Patient-related causes that can lead to deviations from the law of superposition include non-compliance, irregular dosing schedules, and inconsistent dosing amounts [19,20,65]. If the patient does not take the prescribed dosing amount or does not adhere to the dosing schedule, the drug concentration–time curve can exhibit an increase in fluctuation. Similarly, variations in the dosage amount taken by the patient may result in varying drug levels in the body, which can impact the pharmacokinetic behavior of some drugs. Inconsistent dosing sizes may result in varying drug levels in the body, which can lead to the appearance of the superposition principle not being upheld. In all the described cases, the variability in drug administration may cause deviations from the expected linear pharmacokinetic behavior described by the law of superposition.
When patients take medications, some drugs may interact with specific components in the food they ingest or those released during the digestive process, leading to a reduction in the anticipated absorption [66,67]. When patients deviate from prescribed dosage instructions, these interactions can impact the expected pharmacokinetic behavior of drugs, especially when multiple doses are involved. The impact of meals on drug absorption is multifaceted. In particular, the concurrent administration of certain drugs with food has been shown to enhance their absorption. This phenomenon is particularly relevant for drugs that rely on surfactants to facilitate their absorption [66]. Also, lipophilic drugs when taken with a fatty meal can follow another absorption pathway, bypassing the traditional portal vein route and instead entering the systemic circulation via the lymphatic system, alongside absorbed dietary fats [68,69].
Drug–drug interactions stand as another influential factor that can disrupt the adherence to the law of superposition for certain drugs. When certain drugs are co-administered with others, several intricate interactions can come into play, ultimately impacting their pharmacokinetic behavior. One common scenario involves competitive interactions, where co-administered drugs may compete for the same transporters responsible for drug absorption or elimination [70]. This competitive struggle for transporters can lead to deviations from the expected linear pharmacokinetics of the involved drugs. Furthermore, drug interactions can also affect metabolizing enzymes [71]. Some drugs may induce or inhibit these enzymes, which play critical roles in processes like pre-systemic first-pass metabolism and overall liver clearance [72,73]. These enzyme-related changes can significantly influence the pharmacokinetic profile of drugs, causing it to deviate from the anticipated linear behavior described by the law of superposition. Another intricate effect arises from the interaction between co-administered drugs with high plasma protein-binding affinity [74,75]. When two drugs that both extensively bind to plasma proteins are taken together, they can engage in what is known as protein-binding substitution. This process can influence the distribution and clearance of both drugs, as they essentially compete for available protein-binding sites.

4. Conclusions

The pharmacokinetics and area under the concentration–time curve of the majority of drugs after administration of a single dose can be described by first-order or linear processes and can be used to predict a similar steady-state area under the concentration–time curve exposure and, thus, follow the law of superposition. However, there are a number of scenarios and situations where drugs could display behaviors after multiple dosing that leads to capacity-limited or saturation non-linear kinetics and the law of superposition is overruled and not followed. As a consequence, when ADME processes are saturated, the AUC at the steady state may result in substantial changes to the plasma concentration, possibly leading to changes in the pharmacodynamics effect or a toxic effect. Thus, a sound understanding of the principles and influence of dosing (single or multiple dosing) on pharmacokinetics is important in the analysis of pharmacokinetic data by analysts during drug development.

Author Contributions

Conceptualization, M.Y., J.A.Y., R.L. and N.M.D.; methodology M.Y., J.A.Y. and N.M.D.; formal analysis M.Y., J.A.Y., R.L. and N.M.D., investigation M.Y., J.A.Y., R.L. and N.M.D.; resources R.L. and N.M.D.; data curation M.Y. and J.A.Y.; writing—original draft preparation M.Y., J.A.Y., R.L. and N.M.D.; writing—review and editing M.Y., J.A.Y., R.L. and N.M.D.; visualization M.Y., J.A.Y. and N.M.D.; supervision R.L. and N.M.D.; project administration R.L. and N.M.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the MITACS Accelerate Internship (IT40698) for M.Y.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

References

  1. Garza, I.; Swanson, J.W. Prophylaxis of migraine. Neuropsychiatr. Dis. Treat. 2006, 2, 281–291. [Google Scholar] [CrossRef] [PubMed]
  2. Ong, J.J.Y.; De Felice, M. Migraine Treatment: Current Acute Medications and Their Potential Mechanisms of Action. Neurotherapeutics 2018, 15, 274–290. [Google Scholar] [CrossRef] [PubMed]
  3. Diener, H.-C.; Holle-Lee, D.; Nägel, S.; Dresler, T.; Gaul, C.; Göbel, H.; Heinze-Kuhn, K.; Jürgens, T.; Kropp, P.; Meyer, B.; et al. Treatment of migraine attacks and prevention of migraine: Guidelines by the German Migraine and Headache Society and the German Society of Neurology. Clin. Transl. Neurosci. 2019, 3, 3. [Google Scholar] [CrossRef]
  4. Do, T.P.; Guo, S.; Ashina, M. Therapeutic novelties in migraine: New drugs, new hope? J. Headache Pain 2019, 20, 37. [Google Scholar] [CrossRef] [PubMed]
  5. Borbély, A.A.; Mattmann, P.; Loepfe, M.; Fellmann, I.; Gerne, M.; Strauch, I.; Lehmann, D. A single dose of benzodiazepine hypnotics alters the sleep EEG in the subsequent drug-free night. Eur. J. Pharmacol. 1983, 89, 157–161. [Google Scholar] [CrossRef]
  6. Johnson, L.C.; Chernik, D.A. Sedative-hypnotics and human performance. Psychopharmacology 1982, 76, 101–113. [Google Scholar] [CrossRef] [PubMed]
  7. Malpas, A.; Rowan, A.J.; Boyce, C.R.; Scott, D.F. Persistent behavioural and electroencephalographic changes after single doses of nitrazepam and amylobarbitone sodium. Br. Med. J. 1970, 2, 762–764. [Google Scholar] [CrossRef]
  8. Kornetsky, C.; Vates, T.S.; Kessler, E.K. A comparison of hypnotic and residual psychological effects of single doses of chlorpromazine and secobarbital in man. J. Pharmacol. Exp. Ther. 1959, 127, 51–54. [Google Scholar] [CrossRef]
  9. Eckburg, P.B.; Jain, A.; Walpole, S.; Moore, G.; Utley, L.; Manyak, E.; Dane, A.; Melnick, D. Safety, Pharmacokinetics, and Food Effect of Tebipenem Pivoxil Hydrobromide after Single and Multiple Ascending Oral Doses in Healthy Adult Subjects. Antimicrob. Agents Chemother. 2019, 63, e00618-19. [Google Scholar] [CrossRef]
  10. Chua, H.C.; Tse, A.; Smith, N.M.; Mergenhagen, K.A.; Cha, R.; Tsuji, B.T. Combatting the Rising Tide of Antimicrobial Resistance: Pharmacokinetic/Pharmacodynamic Dosing Strategies for Maximal Precision. Int. J. Antimicrob. Agents 2021, 57, 106269. [Google Scholar] [CrossRef]
  11. Xie, J.; Roberts, J.A.; Lipman, J.; Cai, Y.; Wang, H.; Zhao, N.; Xu, X.; Yang, S.; Li, Y.; Zhang, K. Pharmacokinetic/pharmacodynamic adequacy of polymyxin B against extensively drug-resistant Gram-negative bacteria in critically ill, general ward and cystic fibrosis patient populations. Int. J. Antimicrob. Agents 2020, 55, 105943. [Google Scholar] [CrossRef] [PubMed]
  12. Bunke, C.M.; Aronoff, G.R.; Luft, F.C. Pharmacokinetics of common antibiotics used in continuous ambulatory peritoneal dialysis. Am. J. Kidney Dis. 1983, 3, 114–117. [Google Scholar] [CrossRef] [PubMed]
  13. Kang, W.Y.; Lee, H.W.; Gwon, M.R.; Cho, S.; Shim, W.S.; Lee, K.T.; Yang, D.H.; Seong, S.J.; Yoon, Y.R. A Pharmacokinetic Drug Interaction Between Fimasartan and Linagliptin in Healthy Volunteers. Drug Des. Dev. Ther. 2020, 14, 2101–2111. [Google Scholar] [CrossRef] [PubMed]
  14. Moon, S.J.; Yu, K.S.; Kim, M.G. An Assessment of Pharmacokinetic Interaction Between Lobeglitazone and Sitagliptin After Multiple Oral Administrations in Healthy Men. Clin. Ther. 2020, 42, 1047–1057. [Google Scholar] [CrossRef] [PubMed]
  15. Ghim, J.L.; Phuong, N.T.T.; Kim, M.J.; Kim, E.J.; Song, G.S.; Ahn, S.; Shin, J.G.; Kim, E.Y. Pharmacokinetics of fixed-dose combination of atorvastatin and metformin compared with individual tablets. Drug Des. Dev. Ther. 2019, 13, 1623–1632. [Google Scholar] [CrossRef] [PubMed]
  16. Sheng, L.; Cao, W.; Lin, P.; Chen, W.; Xu, H.; Zhong, C.; Yuan, F.; Chen, H.; Li, H.; Liu, C.; et al. Safety, Tolerability and Pharmacokinetics of Single and Multiple Ascending Doses of Benfotiamine in Healthy Subjects. Drug Des. Dev. Ther. 2021, 15, 1101–1110. [Google Scholar] [CrossRef]
  17. Kim, J.R.; Kim, S.; Huh, W.; Ko, J.W. No pharmacokinetic interactions between candesartan and amlodipine following multiple oral administrations in healthy subjects. Drug Des. Dev. Ther. 2018, 12, 2475–2483. [Google Scholar] [CrossRef]
  18. Vaidyanathan, S.; Jarugula, V.; Dieterich, H.A.; Howard, D.; Dole, W.P. Clinical pharmacokinetics and pharmacodynamics of aliskiren. Clin. Pharmacokinet. 2008, 47, 515–531. [Google Scholar] [CrossRef]
  19. Wang, W.; Husan, F.; Chow, S.C. The impact of patient compliance on drug concentration profile in multiple doses. Stat. Med. 1996, 15, 659–669. [Google Scholar] [CrossRef]
  20. Wang, W. Patient Compliance and its Impact on Steady State Pharmacokinetics. In Applied Statistics in the Pharmaceutical Industry: with Case Studies Using S-Plus, 1st ed.; Millard, S.P., Krause, A., Eds.; Springer: New York, NY, USA, 2001; pp. 217–236. [Google Scholar] [CrossRef]
  21. van Rossum, J.M. Pharmacokinetics of accumulation. J. Pharm. Sci. 1968, 57, 2162–2165. [Google Scholar] [CrossRef]
  22. Wang, W.; Ouyang, S.P. The formulation of the principle of superposition in the presence of non-compliance and its applications in multiple dose pharmacokinetics. J. Pharmacokinet. Biopharm. 1998, 26, 457–469. [Google Scholar] [CrossRef] [PubMed]
  23. Lin, J.H. Dose-dependent pharmacokinetics: Experimental observations and theoretical considerations. Biopharm. Drug Dispos. 1994, 15, 1–31. [Google Scholar] [CrossRef] [PubMed]
  24. Cutler, D.J. Numerical deconvolution by least squares: Use of prescribed input functions. J. Pharmacokinet. Biopharm. 1978, 6, 227–241. [Google Scholar] [CrossRef] [PubMed]
  25. Cutler, D.J. Numerical deconvolution by least squares: Use of polynomials to represent the input function. J. Pharmacokinet. Biopharm. 1978, 6, 243–263. [Google Scholar] [CrossRef] [PubMed]
  26. Fernández-Campos, F.; Ferrero, C.; Colom, H.; Jiménez-Castellanos, M.R. In vivo absorption behaviour of theophylline from starch-methyl methacrylate matrix tablets in beagle dogs. Int. J. Pharm. 2015, 478, 684–692. [Google Scholar] [CrossRef] [PubMed]
  27. Chiou, W.L. Rapid compartment- and model-independent estimation of times required to attain various fractions of steady-state plasma level during multiple dosing of drugs obeying superposition principle and having various absorption or infusion kinetics. J. Pharm. Sci. 1979, 68, 1546–1547. [Google Scholar] [CrossRef] [PubMed]
  28. Gupta, P.; Hutmacher, M.M.; Frame, B.; Miller, R. An alternative method for population pharmacokinetic data analysis under noncompliance. J. Pharmacokinet. Pharmacodyn. 2008, 35, 219–233. [Google Scholar] [CrossRef] [PubMed]
  29. Ma, L. Analysis of Nonlinear Pharmacokinetic Systems and the Nonlinear Disposition of Phenylbutazone in Equine (Horses). Ph.D. Thesis, Oregon State University, Corvallis, OR, USA, 13 June 2013. [Google Scholar]
  30. McCoy, A.T.; Bartels, M.J.; Rick, D.L.; Saghir, S.A. TK Modeler version 1.0, a Microsoft® Excel®-based modeling software for the prediction of diurnal blood/plasma concentration for toxicokinetic use. Regul. Toxicol. Pharmacol. 2012, 63, 333–343. [Google Scholar] [CrossRef]
  31. Mengozzi, G.; Intorre, L.; Bertini, S.; Giorgi, M.; Soldani, G. Comparative bioavailability of two sustained-release theophylline formulations in the dog. Pharmacol. Res. 1998, 38, 481–485. [Google Scholar] [CrossRef]
  32. Thron, C.D. Linearity and superposition in pharmacokinetics. Pharmacol. Rev. 1974, 26, 3–31. [Google Scholar]
  33. Tuntland, T.; Ethell, B.; Kosaka, T.; Blasco, F.; Zang, R.X.; Jain, M.; Gould, T.; Hoffmaster, K. Implementation of pharmacokinetic and pharmacodynamic strategies in early research phases of drug discovery and development at Novartis Institute of Biomedical Research. Front. Pharmacol. 2014, 5, 174. [Google Scholar] [CrossRef]
  34. van Rossum, J.M.; de Bie, J.E. Systems dynamics in clinical pharmacokinetics. An introduction. Clin. Pharmacokinet. 1989, 17, 27–44. [Google Scholar] [CrossRef]
  35. Peletier, L.A.; de Winter, W. Impact of saturable distribution in compartmental PK models: Dynamics and practical use. J. Pharmacokinet. Pharmacodyn. 2017, 44, 1–16. [Google Scholar] [CrossRef] [PubMed]
  36. Gabrielsson, J.; Meibohm, B.; Weiner, D. Pattern Recognition in Pharmacokinetic Data Analysis. AAPS J. 2016, 18, 47–63. [Google Scholar] [CrossRef] [PubMed]
  37. Mehrotra, N.; Gupta, M.; Kovar, A.; Meibohm, B. The role of pharmacokinetics and pharmacodynamics in phosphodiesterase-5 inhibitor therapy. Int. J. Impot. Res. 2007, 19, 253–264. [Google Scholar] [CrossRef] [PubMed]
  38. Shen, J.; Boeckmann, A.; Vick, A. Implementation of dose superimposition to introduce multiple doses for a mathematical absorption model (transit compartment model). J. Pharmacokinet. Pharmacodyn. 2012, 39, 251–262. [Google Scholar] [CrossRef] [PubMed]
  39. Siwale, R.C.; Sani, S.N. Multiple-Dosage Regimens. In Applied Biopharmaceutics & Pharmacokinetics, 8th ed.; Ducharme, M.P., Shargel, L., Eds.; McGraw Hill Education: New York, NY, USA, 2022. [Google Scholar]
  40. Mehvar, R. Dependence of time to reach steady-state on the length of dosage interval. Ann. Pharmacother. 2008, 42, 1518–1519. [Google Scholar] [CrossRef]
  41. Hieb, B.R.; Shrewsbury, B. Consecutive intravenous infusions: Simulation of two compartment pharmacokinetic drugs. Comput. Programs Biomed. 1980, 12, 96–104. [Google Scholar] [CrossRef] [PubMed]
  42. Brocks, D.R.; Mehvar, R. Rate and extent of drug accumulation after multiple dosing revisited. Clin. Pharmacokinet. 2010, 49, 421–438. [Google Scholar] [CrossRef]
  43. van Rossum, J.M.; van Lingen, G.; Burgers, J.P.T. Dose-dependent pharmacokinetics. Pharmacol. Ther. 1983, 21, 77–99. [Google Scholar] [CrossRef]
  44. Nation, R.L.; Evans, A.M.; Milne, R.W. Pharmacokinetic drug interactions with phenytoin (Part I). Clin. Pharmacokinet. 1990, 18, 37–60. [Google Scholar] [CrossRef]
  45. Yáñez, J.A.; Remsberg, C.M.; Sayre, C.L.; Forrest, M.L.; Davies, N.M. Flip-flop pharmacokinetics–delivering a reversal of disposition: Challenges and opportunities during drug development. Ther. Deliv. 2011, 2, 643–672. [Google Scholar] [CrossRef]
  46. Stewart, B.H.; Kugler, A.R.; Thompson, P.R.; Bockbrader, H.N. A saturable transport mechanism in the intestinal absorption of gabapentin is the underlying cause of the lack of proportionality between increasing dose and drug levels in plasma. Pharm. Res. 1993, 10, 276–281. [Google Scholar] [CrossRef]
  47. Proctor, W.R.; Bourdet, D.L.; Thakker, D.R. Mechanisms underlying saturable intestinal absorption of metformin. Drug Metab. Dispos. 2008, 36, 1650–1658. [Google Scholar] [CrossRef]
  48. Hanada, K.; Ikemi, Y.; Kukita, K.; Mihara, K.; Ogata, H. Stereoselective first-pass metabolism of verapamil in the small intestine and liver in rats. Drug Metab. Dispos. 2008, 36, 2037–2042. [Google Scholar] [CrossRef]
  49. Kalam, M.N.; Rasool, M.F.; Rehman, A.U.; Ahmed, N. Clinical pharmacokinetics of propranolol hydrochloride: A review. Curr. Drug Metab. 2020, 21, 89–105. [Google Scholar] [CrossRef]
  50. Karlsson, M.O.; Molnar, V.; Freijs, A.; Nygren, P.; Bergh, J.; Larsson, R. Pharmacokinetic models for the saturable distribution of paclitaxel. Drug Metab. Dispos. 1999, 27, 1220–1223. [Google Scholar]
  51. Eling, T.E.; Pickett, R.D.; Orton, T.C.; Anderson, M.W. A study of the dynamics of imipramine accumulation in the isolated perfused rabbit lung. Drug Metab. Dispos. 1975, 3, 389–399. [Google Scholar]
  52. Ha, H.R.; Chen, J.; Freiburghaus, A.U.; Follath, F. Metabolism of theophylline by cDNA-expressed human cytochromes P-450. Br. J. Clin. Pharmacol. 1995, 39, 321–326. [Google Scholar] [CrossRef]
  53. Bernus, I.; Dickinson, R.G.; Hooper, W.D.; Eadie, M.J. Dose-dependent metabolism of carbamazepine in humans. Epilepsy Res. 1996, 24, 163–172. [Google Scholar] [CrossRef]
  54. Patel, B.A.; Chu, C.K.; Boudinot, F.D. Pharmacokinetics and saturable renal tubular secretion of zidovudine in rats. J. Pharm. Sci. 1989, 78, 530–534. [Google Scholar] [CrossRef] [PubMed]
  55. Hendel, J.; Nyfors, A. Nonlinear renal elimination kinetics of methotrexate due to saturation of renal tubular reabsorption. Eur. J. Clin. Pharmacol. 1984, 26, 121–124. [Google Scholar] [CrossRef] [PubMed]
  56. Babu, E.; Takeda, M.; Narikawa, S.; Kobayashi, Y.; Yamamoto, T.; Cha, S.H.; Sekine, T.; Sakthisekaran, D.; Endou, H. Human organic anion transporters mediate the transport of tetracycline. Jpn. J. Pharmacol. 2002, 88, 69–76. [Google Scholar] [CrossRef] [PubMed]
  57. Hornecker, M.; Blanchet, B.; Billemont, B.; Sassi, H.; Ropert, S.; Taieb, F.; Mir, O.; Abbas, H.; Harcouet, L.; Coriat, R.; et al. Saturable absorption of sorafenib in patients with solid tumors: A population model. Investig. New Drug. 2012, 30, 1991–2000. [Google Scholar] [CrossRef] [PubMed]
  58. Graham, R.A.; Hop, C.E.C.A.; Borin, M.T.; Lum, B.L.; Colburn, D.; Chang, I.; Shin, Y.G.; Malhi, V.; Low, J.A.; Dresser, M.J. Single and multiple dose intravenous and oral pharmacokinetics of the hedgehog pathway inhibitor vismodegib in healthy female subjects. Br. J. Clin. Pharmacol. 2012, 74, 788–796. [Google Scholar] [CrossRef]
  59. Yoon, Y.-R.; Shin, J.-G.; Cha, I.-J.; Kim, K.-A.; Shim, J.-C.; Kim, Y.-H.; Shin, J.-B. Pharmacokinetic analysis on autoinduction of carbamazepine metabolism. J. Korean Soc. Clin. Pharmacol. Ther. 1996, 4, 139–147. [Google Scholar] [CrossRef]
  60. Chetty, M.; Miller, R.; Seymour, M.A. Phenytoin auto-induction. Ther. Drug Monit. 1998, 20, 60–62. [Google Scholar] [CrossRef]
  61. Shand, D.G.; Hammill, S.C.; Aanonsen, L.; Pritchett, E.L.C. Reduced verapamil clearance during long-term oral administration. Clin. Pharmacol. Ther. 1981, 30, 701–703. [Google Scholar] [CrossRef]
  62. Abernethy, D.R.; Wainer, I.W.; Anacleto, A.I. Verapamil metabolite exposure in older and younger men during steady-state oral verapamil administration. Drug Met. Dispos. 2000, 28, 760–765. [Google Scholar]
  63. Lalonde, R.L.; Pieper, J.A.; Straka, R.J.; Bottorff, M.B.; Mirvis, D.M. Propranolol pharmacokinetics and pharmacodynamics after single doses and at steady-state. Eur. J. Clin. Pharmacol. 1987, 33, 315–318. [Google Scholar] [CrossRef]
  64. Barone, J.A.; Koh, J.G.; Bierman, R.H.; Colaizzi, J.L.; Swanson, K.A.; Gaffar, M.C.; Moskovitz, B.L.; Mechlinski, W.; Van de Velde, V. Food interaction and steady-state pharmacokinetics of itraconazole capsules in healthy male volunteers. Antimicrob. Agents Chemother. 1993, 37, 778–784. [Google Scholar] [CrossRef]
  65. Vrijens, B.; Goetghebeur, E. The impact of compliance in pharmacokinetic studies. Stat. Methods Med. Res. 1999, 8, 247–262. [Google Scholar] [CrossRef]
  66. Martinez, M.N.; Amidon, G.L. A mechanistic approach to understanding the factors affecting drug absorption: A review of fundamentals. J. Clin. Pharmacol. 2002, 42, 620–643. [Google Scholar] [CrossRef]
  67. Singh, B.N. Effects of food on clinical pharmacokinetics. Clin. Pharmacokinet. 1999, 37, 213–255. [Google Scholar] [CrossRef]
  68. Gershkovich, P.; Hoffman, A. Effect of a high-fat meal on absorption and disposition of lipophilic compounds: The importance of degree of association with triglyceride-rich lipoproteins. Eur. J. Pharm. Sci. 2007, 32, 24–32. [Google Scholar] [CrossRef]
  69. Yousef, M.; Silva, D.; Bou-Chacra, N.; Davies, N.M.; Löbenberg, R. The lymphatic system: A sometimes forgotten compartment in pharmaceutical sciences. J. Pharm. Pharm. Sci. 2021, 24, 533–547. [Google Scholar] [CrossRef]
  70. König, J.; Müller, F.; Fromm, M.F. Transporters and drug-drug interactions: Important determinants of drug disposition and effects. Pharmacol. Rev. 2013, 65, 944–966. [Google Scholar] [CrossRef]
  71. Tornio, A.; Filppula, A.M.; Niemi, M.; Backman, J.T. Clinical studies on drug–drug interactions involving metabolism and transport: Methodology, pitfalls, and interpretation. Clin. Pharmacol. Ther. 2019, 105, 1345–1361. [Google Scholar] [CrossRef]
  72. Hall, S.D.; Thummel, K.E.; Watkins, P.B.; Lown, K.S.; Benet, L.Z.; Paine, M.F.; Mayo, R.R.; Turgeon, D.K.; Bailey, D.G.; Fontana, R.J.; et al. Molecular and physical mechanisms of first-pass extraction. Drug Met. Dispos. 1999, 27, 161–166. [Google Scholar]
  73. Bachmann, K. Drug–drug interactions with an emphasis on drug metabolism and transport. In Pharmacology Principles and Practice, Hacker, M., Messer, W., Bachmann, K., Eds.; Academic Press: Cambridge, MA, USA, 2009; pp. 303–325. [Google Scholar] [CrossRef]
  74. DeVane, C.L. Clinical significance of drug binding, protein binding, and binding displacement drug interactions. Psychopharmacol. Bull. 2002, 36, 5–21. [Google Scholar]
  75. Heuberger, J.; Schmidt, S.; Derendorf, H. When is protein binding important? J. Pharm. Sci. 2013, 102, 458–3467. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Concentration time profile after a single dose and multiple dosing when reaching the steady state.
Figure 1. Concentration time profile after a single dose and multiple dosing when reaching the steady state.
Biomedicines 12 01843 g001
Figure 2. Concentration–time profile of an ideal scenario after multiple-dose regimen 132, administered at the same time interval and reaching a steady state. Modified from Wang et al. [19].
Figure 2. Concentration–time profile of an ideal scenario after multiple-dose regimen 132, administered at the same time interval and reaching a steady state. Modified from Wang et al. [19].
Biomedicines 12 01843 g002
Figure 3. Representation of a concentration–time profile after six sequential doses administered at the same time interval until the steady state. Modified from Wang et al. [19].
Figure 3. Representation of a concentration–time profile after six sequential doses administered at the same time interval until the steady state. Modified from Wang et al. [19].
Biomedicines 12 01843 g003
Figure 4. The law of superposition used to generate the concentration profile of multiple dosing from the profile of a single dose. The area of the first dose during the interval of the second dose adds to the area of the second dose, and so on. Modified from Van Rossum and de Bie [34].
Figure 4. The law of superposition used to generate the concentration profile of multiple dosing from the profile of a single dose. The area of the first dose during the interval of the second dose adds to the area of the second dose, and so on. Modified from Van Rossum and de Bie [34].
Biomedicines 12 01843 g004
Figure 5. (a) Linearity of pharmacokinetics and the validity of the superposition principle, along with (b) scenarios where non-linearity in pharmacokinetics leads to the breakdown of this principle.
Figure 5. (a) Linearity of pharmacokinetics and the validity of the superposition principle, along with (b) scenarios where non-linearity in pharmacokinetics leads to the breakdown of this principle.
Biomedicines 12 01843 g005
Table 1. Concentration of Canadamycin (a hypothetical drug) in mg/L before each subsequent multiple dose is administered.
Table 1. Concentration of Canadamycin (a hypothetical drug) in mg/L before each subsequent multiple dose is administered.
DoseCanadamycin (mg/L)Start
Multiple
Trough End (r)Time (h)Concentration Lost during Dosage Interval
120.0020.005.0815.00
25.0025.006.251618.75
31.2526.256.562419.69
40.312526.566.643219.92
50.07812526.646.664019.98
60.0195312526.666.674819.99
70.00488281226.676.675420.00
120.0020.005.0815.00
Table 2. Primary potential causes for deviations from the law of superposition in various ADME processes, with examples.
Table 2. Primary potential causes for deviations from the law of superposition in various ADME processes, with examples.
ADME ProcessMain Related Causes of Non-LinearityExamples
AbsorptionCarrier-mediated transportGabapentin [46], Metformin [47]
Saturable pre-systemic lossGut metabolismVerapamil [48]
First-pass metabolismPropranolol [49]
DistributionSaturable plasma protein bindingPaclitaxel [50]
Saturable tissue bindingImipramine [51]
EliminationMetabolismSaturable metabolismTheophylline [52]
Enzyme inductionCarbamazepine [53]
ExcretionSaturable renal excretionSaturable active secretionZidovudine [54]
Saturable reabsorptionMethotrexate [55]
Saturable biliary excretionTetracycline [56]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yousef, M.; Yáñez, J.A.; Löbenberg, R.; Davies, N.M. Upholding or Breaking the Law of Superposition in Pharmacokinetics. Biomedicines 2024, 12, 1843. https://doi.org/10.3390/biomedicines12081843

AMA Style

Yousef M, Yáñez JA, Löbenberg R, Davies NM. Upholding or Breaking the Law of Superposition in Pharmacokinetics. Biomedicines. 2024; 12(8):1843. https://doi.org/10.3390/biomedicines12081843

Chicago/Turabian Style

Yousef, Malaz, Jaime A. Yáñez, Raimar Löbenberg, and Neal M. Davies. 2024. "Upholding or Breaking the Law of Superposition in Pharmacokinetics" Biomedicines 12, no. 8: 1843. https://doi.org/10.3390/biomedicines12081843

APA Style

Yousef, M., Yáñez, J. A., Löbenberg, R., & Davies, N. M. (2024). Upholding or Breaking the Law of Superposition in Pharmacokinetics. Biomedicines, 12(8), 1843. https://doi.org/10.3390/biomedicines12081843

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop