Skip to main content
Advertisement

Main menu

  • Home
  • Articles
    • Current Issue
    • Fast Forward
    • Latest Articles
    • Special Sections
    • Archive
  • Information
    • Instructions to Authors
    • Submit a Manuscript
    • FAQs
    • For Subscribers
    • Terms & Conditions of Use
    • Permissions
  • Editorial Board
  • Alerts
    • Alerts
    • RSS Feeds
  • Virtual Issues
  • Feedback
  • Submit
  • Other Publications
    • Drug Metabolism and Disposition
    • Journal of Pharmacology and Experimental Therapeutics
    • Molecular Pharmacology
    • Pharmacological Reviews
    • Pharmacology Research & Perspectives
    • ASPET

User menu

  • My alerts
  • Log in
  • My Cart

Search

  • Advanced search
Drug Metabolism & Disposition
  • Other Publications
    • Drug Metabolism and Disposition
    • Journal of Pharmacology and Experimental Therapeutics
    • Molecular Pharmacology
    • Pharmacological Reviews
    • Pharmacology Research & Perspectives
    • ASPET
  • My alerts
  • Log in
  • My Cart
Drug Metabolism & Disposition

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Fast Forward
    • Latest Articles
    • Special Sections
    • Archive
  • Information
    • Instructions to Authors
    • Submit a Manuscript
    • FAQs
    • For Subscribers
    • Terms & Conditions of Use
    • Permissions
  • Editorial Board
  • Alerts
    • Alerts
    • RSS Feeds
  • Virtual Issues
  • Feedback
  • Submit
  • Visit dmd on Facebook
  • Follow dmd on Twitter
  • Follow ASPET on LinkedIn
Research ArticleCommentary

Finding Tmax and Cmax in Multicompartmental Models

Yi Rang Han, Ping I. Lee and K. Sandy Pang
Drug Metabolism and Disposition November 2018, 46 (11) 1796-1804; DOI: https://doi.org/10.1124/dmd.118.082636
Yi Rang Han
Department of Pharmaceutical Sciences, Leslie Dan Faculty of Pharmacy, University of Toronto, Toronto, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Yi Rang Han
Ping I. Lee
Department of Pharmaceutical Sciences, Leslie Dan Faculty of Pharmacy, University of Toronto, Toronto, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ping I. Lee
K. Sandy Pang
Department of Pharmaceutical Sciences, Leslie Dan Faculty of Pharmacy, University of Toronto, Toronto, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for K. Sandy Pang
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Drug absorption data are critical in bioequivalence comparisons, and factors such as the maximum drug concentration (Cmax), time to achieve Cmax (or Tmax), as well as the area under the curve (AUC) are important metrics. It is generally accepted that the AUC is a meaningful estimate of the extent of absorption, and Tmax or Cmax may be used for assessing the rate of absorption. But estimation of the rate of absorption with Tmax or Cmax is not always feasible, as explicit solutions relating Tmax and Cmax to the absorption (ka) and elimination rate (k) constants exist only for the one and not multicompartmental oral model. Therefore, the determination of Tmax or Cmax for multicompartmental models is uncertain. Here, we propose an alternate, numerical approach that uses the point-slope method for the first and second derivative(s) of the concentration-versus-time profiles and the Newton-Raphson iteration method for the determination of Tmax and Cmax. We show that the method holds for multicompartmental oral dosing under single or steady-state conditions in the absence of known microconstants, even for flip-flop (ka < β) models. Simulations showed that the Cmax and Tmax estimates obtained with the Newton-Raphson method were more accurate than those based on the noncompartmental, observation-based method recommended by the US Food and Drug Administration. The %Bias attributable to sampling frequency and assay error were less than those determined by the noncompartmental method, showing that the Newton-Raphson method is viable for the estimation of Tmax and Cmax.

Introduction

Drug absorption data are often used in bioequivalence comparisons. The US Food and Drug Administration (FDA) recommends two measures: the area under the curve (AUC) and the maximum drug concentration (Cmax) (CDER/FDA, 2015), with AUC as the primary and Cmax as the secondary measure. The AUC is known to be independent of ka, the absorption rate constant under first-order conditions that reflects the extent of drug absorption, and is normally calculated by the trapezoidal rule, which is straightforward for both compartmental and noncompartmental models. The FDA defines Cmax as the maximum drug concentration observed in the sampled blood or plasma data, and is henceforth denoted here as Cmax,obs. This Cmax,obs partly serves to provide insight into the rate of absorption. However, an understanding of the absorption rate is more complex. As previously noted (Endrenyi and Al-Shaikh, 1995; Basson et al., 1996), the Cmax is affected by both the extent and the rate of absorption. In addition, the value of Cmax,obs is further influenced by the frequency of the sampling scheme and the magnitude of the assay errors. As the sampling frequency or assay error is increased, Cmax,obs is increased as well (Endrenyi and Al-Shaikh, 1995). Thus, the comparison between Cmax,obs values that are generated by two divergent sampling schemes or assaying methods may be inappropriate. An alternate measure of the absorption rate is the time to reach Cmax,obs or Tmax,obs, which has been suggested to be an unconfounded metric for the rate of absorption (Basson et al., 1996). However, the Tmax,obs is a categorical variable, one that can take on a limited and usually fixed number of possible values and its discriminating power depends strongly again on the sampling frequency.

Generally speaking, values for Cmax and Tmax may be obtained by fitting data to compartmental models and then computing those values on the basis of the fitted parameters. In the case of the one-compartment model, the model-guided determination of Tmax,1comp and Cmax,1comp is clear. The concentration C1comp may be expressed in terms of the absorption rate constant, ka, and elimination rate constant, k, in a biexponential expression, as follows:Embedded Image(1)where V is the central volume of the distribution. The condition at Tmax,1comp yields the Cmax,1comp.Embedded Image(1A)Explicit solution for Tmax,1comp is obtained by solving dC1comp/dt = 0 (eq. 2).Embedded Image(2)The rate constant, ka, may be estimated via curve stripping or with the Wagner-Nelson method (Wagner and Nelson, 1963) to obtain the fraction remaining to be absorbed (FRA), Embedded Image. Together with k and Embedded Image, Cmax,1comp and Tmax,1comp may be calculated (eqs. 1 and 2). The equation holds in the case of flip-flop when ka < k.

The one-compartment analytical equation for finding Tmax and Cmax can be extended to a scenario with absorption lag time. This may be derived by substituting the observed time (t + tlag) and then solving the analogous problem of the form, dC1comp,lag/dt = 0 in eq. 2. The one-compartment concentration profile is:Embedded Image(3)and the associated Tmax is:Embedded Image(4)For the two-compartmental model, a triexponential expression now describes the drug concentration with oral dosing (Gibaldi and Perrier, 1982):Embedded Image(5)An accurate estimation of the parameters is more complex for such a triexponential expression and for the two-compartment, oral case (eq. 5), several scenarios are possible. First is the scenario when ka is much faster than the hybrid decay constant β, i.e, ka > > β, and the three exponential components are easily and accurately separated upon curve stripping (Gibaldi and Perrier, 1982). With k12, first-order transfer rate from compartment 2 to 1 (k21), k10, rate constant for removal in compartment 1 and V1, the blood/plasma volume of distribution of drug in the central compartment known from intravenous dosing, ka can be obtained by the Loo-Riegelman method (Loo and Riegelman, 1968) with the usual assumption that elimination occurs in the central compartment. However, for cases when β ≥ ka, the distinguishing feature of a prominent “hat” or “nose” that occurs with rapid absorption and slow elimination now disappears, and the oral profile shrinks to one resembling that for the one-compartment model (Chan and Gibaldi, 1985). Under this circumstance, the ka so obtained by the curve stripping procedure is no longer accurate. Even if the parameters for the two-compartment model can be accurately determined, the analysis of Tmax and Cmax in drugs that exhibit multiple compartmental characteristics, however, cannot be solved explicitly. Often, the Loo-Riegelman method (Loo and Riegelman, 1968) is used first to obtain ka, then Cmax and Tmax, but the microconstants k12, k21, and k10 and the assumption that elimination occurs from the central compartment are required.

In this article, we describe the use of a numerical approach for model-guided determination of Tmax and Cmax in multicompartment models, and compare the estimates obtained to those from direct observations (Tmax,obs and Cmax,obs), a method encouraged by regulatory agencies such as the FDA. This method utilizes the second derivative of the concentration-versus-time profile and the point-slope and Newton-Raphson methods. The bias resulting from assay error and sampling frequency, as well as its performance for different ka values, are discussed.

Theory and Methods

Finding Tmax and Cmax in Multicompartment Models

The Embedded Image and Embedded Image for oral dosing for a two-compartment model may be found in a manner similar to that described for the one-compartment. The Cmax,2comp occurs when the first derivative (eq. 6) is zero.Embedded Image(6)The above problem may be solved numerically using the Newton-Raphson method (Galantai, 2000), since the second derivative, (Embedded Image) may be easily computed:Embedded Image(7)We will illustrate the solution with examples, first with the two-compartment model, then with multicompartmental models, for drugs with a fast versus a slow absorption rate constant.

Example 1: Oral Dosing for the Two-Compartment Model Drugs

The following example demonstrates use of this method to obtain Cmax,2comp and Tmax,2comp. The following parameters were selected for the scenario of fast versus slow absorption: Embedded Image (fast) or 0.1 (slow) h−1, with common values for β = 0.2 h−1, α = 0.8 h−1, k21 = 0.5 h–1, and [FsysDosepo/V1] = 100. Microconstants k10, k12, and k21 will remain unchanged when β, α, and k21 are kept constant, owing to the relations that exist between α, β, k10, k12, and k21 (α*β = k10*k21, and α + β = k12 + k10 + k21). In fact, the hybrid constants correspond to k10, k12, and k21 values of 0.32, 0.18, and 0.5 h−1, respectively. After substitution of the selected values into eq. 5, values of the drug concentration, C2comp, in the central compartment, the derivatives Embedded Image and Embedded Image, and the Tmax,2comp estimates were computed and summarized in Table 1. The log-linear concentration-time profile and the two-compartmental model are shown in Fig. 1A showing a prominent “nose,” and components of absorption (ka), distribution (α), and elimination (β) for fast absorption, or an apparent one-compartment profile for slow absorption (Fig. 2A). For fast absorption, the ka may be obtained by the method of residuals. However, for slow absorption (Fig. 2A), ka may only be estimated with the Loo and Riegelman method (Loo and Riegelman, 1968) when the microconstants k12, k21, and k10 are known. In contrast, for both fast and slow absorption, the Newton-Raphson method (Galantai, 2000) may be used to obtain Tmax,2comp without knowledge of the microconstants.

View this table:
  • View inline
  • View popup
TABLE 1

Results (error free; assay error E = 0) arising from the Newton-Raphson method for estimating of Tmax from two- or three-compartmental models

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

The log-linear blood/plasma concentration-versus-time plot for a two-compartment oral dosing case, with elimination from central compartment, for drugs exhibiting fast absorption with ka of 2 h−1 and >β. Parameter values were: α, β, and k21 are 0.8, 0.2, and 0.5 h−1, respectively. (A) The dashed line is an extrapolated log linear terminal slope, β; the method of residuals would further provide α and ka. The inset shows an expanded, zoomed-in version. (B) Plot of the first derivative versus time with the first two iterations or Testimate are highlighted in green and blue; the inset shows the zoomed-out version of the same plot. A total of five iterations were required and these were replotted on (A) and the inset (see text for details), and (C) an expanded time plot of (B). Initial estimates, when chosen from the green and not red region in (C), highlight the convergence of the Newton-Raphson method with the true or correct value. However, initial estimates chosen from the red region will not converge but approach infinity upon reaching the iteration limit.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

The log-linear blood/plasma concentration-versus-time plot for a two-compartment oral dosing case, with elimination from central compartment, for drugs exhibiting flip-flop kinetics (slow ka = 0.1 h−1 < β). Parameter values were: α, β, and k21 are 0.8, 0.2, and 0.5 h−1, respectively, and values for α, β, and k21 were identical to those for Fig. 1. (A) The colored square symbols are the iteration points of the Newton-Raphson method; a zoomed-in plot is shown in the inset. (B) Plot of the first derivative as a function of time; In (B), the green and blue symbols are points used to obtain the first two iterations of the Newton-Raphson method (see text for details). (C) Zoomed-out version of the same plot. Initial estimates made within the green but not red shaded region highlights where Newton-Raphson method will converge with the correct value.

Fast ka.

A time-expanded view of Fig. 1B showed that the dC2comp/dt profile intersected at the zero line (obtained when eq. 6 = 0), then continued negatively before rebounding upwards, rising asymptotically to parallel and approach 0. The first derivative (Fig. 1B), when set = 0 in eq. 6, provided the Tmax,2comp estimate, which equaled 0.97 h graphically from Fig. 1B (see value in Table 1). Such a fast absorption can occur in the case of immediate-release oral dosage forms of morphine and cefaclor (Barbhaiya et al., 1990; Collins et al., 1998). Even faster Tmax is possible for bolus subcutaneous routes of morphine and other drugs (Home et al., 1999; Stuart-Harris et al., 2000), for which the method presented here is also applicable.

The iterative scheme of the Newton-Raphson method to find where the function is zero is illustrated in detail in Fig. 1B. An initial time estimate (Testimate1) is required for the numerical method. For demonstration purposes, an arbitrary time of 0.25 h was selected for Testimate1. For practical applications, Tmax,obs may be used as the initial estimate for Testimate1. At Testimate1, the tangent line of the function dC2comp/dt was constructed using the “point-slope method,” where the slope was essentially Embedded Image evaluated at Testimate1. From the line equation of the tangent (y = mx + b), the x-intercept at y = 0 (i.e., the time value a1) was computed, yielding the corresponding Testimate2 on the dC2comp/dt curve (see Fig. 1B, green). The point-slope procedure was repeated as shown in the blue point in Fig. 1B to provide a2 and the third estimate, Testimate3. The continued iterations lead to Testimate7 value for which the difference between dC2comp/dt and 0 is less than a chosen tolerance level (ε). When this criterion is met, the numerical method has converged to the true time needed to reach the maximum drug concentration in the blood/plasma (central) compartment (Tmax,true).

Convergence may be achieved numerically in the critical point that may not be of interest. We illustrate this with the scenarios below in Fig. 1C, which shows the behavior of the dC2Comp/dt-versus-time curve in a more expanded time scale. When the initial Testimate1 (within the shaded green region) was appropriately selected before the minimum for the dC2comp/dt curve (at 1.9 h, or the boundary between shaded green and red regions in Fig. 1C), one would obtain the true Tmax,2comp,true at convergence. However, if the initial Testimate1 was selected from the red zone (Fig. 1C), the numerical method would not converge to Tmax,2comp,true. Hence, it is important to select an appropriate value for Testimate1. For good initial estimates between 0 and 1.9 h (where Embedded Image is 0 in Fig. 1C, with the Tmax,2comp,true at 0.97 h), the numerical method would converge and yield the true Tmax,2comp value at the minimum. Estimates larger than 1.9 h, however, may lead to the asymptotic value of 0 (Fig. 1C). When all the Testimate values are plotted on the C-versus-time plot (shown as colored points in inset of Fig. 1A), it may be seen that the estimates readily approach Tmax,2comp,true. With known Tmax,2comp,true, the corresponding Cmax,2comp,true may be calculated (eq. 5).

Slow ka.

By contrast, the log-linear concentration profile and the two-compartmental model with slow ka, shown in Fig. 2A, revealed an apparent, one-compartment characteristic. The absorption constant (ka) could be obtained only with the Loo and Riegelman (1968) method with known microconstants k12, k21, and k10 or simultaneously curve fitting with intravenous data involving these microconstants. Again, the iterative Newton-Raphson method could be readily applied here. By selecting Testimate1, finding a1, then Testimate2, the point-slope procedure was applied to provide a2 and Testimate3, then eventually Tmax,2comp,true (Fig. 2C, showing comparable green and red regions as in Fig. 1C), similar to the iteration scheme of Fig. 1 when ka is faster (2 h−1). The only difference here was that the procedure provided a “shallow well” when Embedded Image was set = 0 in eq. 7 (Fig. 2C). The first derivative (Fig. 2C), when set = 0, revealed the true Tmax,2comp estimate of 5.8 h−1. The points on C2comp-versus-time plot (Fig. 2A) showed the progression in approaching the true Cmax,2comp and Tmax,2comp when ka is slow (0.1 h−1). Again, when the initial Testimate1 (within the shaded green region, Fig. 2C) was appropriately selected before the minimum for the dC2comp/dt curve (at 12.5 h, or boundary between shaded green and red regions in Fig. 2C), then one would obtain the true Tmax,2comp,true at convergence. For good initial estimates between 0 and 12.5 h (where Embedded Image is 0 in Fig. 2C, with the Tmax,2comp,true at 5.77 h), the numerical method would converge and yield the true Tmax,2comp value at the minimum. Initial estimates larger than 12.5 h, however, may lead to the asymptotic value of 0 (Fig. 2C). Multiple examples of a long Tmax have been observed for the slow, oral absorption profiles of drugs such as sertraline (>7 h) (Allard et al., 1999) and zonisamide (5–6 h) (Mimaki, 1998).

Figures 1 and 2 both illustrate that the Newton-Raphson method works well regardless whether ka > β, or ka < β.

Example 2: Oral Dosing for Multicompartment Model Drugs

The general method is applicable not only to examine temporal data after single oral dosing for the two-compartment model cases but also for multiple dosing at steady-state. The concentration for multicompartmental models is shown below, where hybrid constants are denoted as λi:Embedded Image(8)Embedded ImageTo demonstrate this, a three-compartment, single oral dose model, with elimination from central compartment, was simulated for the case where ka > λ1 and ka < λi. The selected values are summarized in Table 1, with the method for estimation of Tmax,3comp the same as that for the two-compartment model (see Fig. 3). With known Tmax,3comp,true, the corresponding Cmax,3comp,true may be calculated.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

The log-linear blood/plasma concentration-versus-time plot for the three-compartment oral dosing, with elimination from central compartment. The right panel shows the derivative dC3comp/dt plot. The method is similar to that described for the two-compartment model. Parameter values for λ1, λ2, λ3, k21, and k31 are 1.0, 0.8, 0.2, 1.5, and 0.5 h–1, respectively, with ka of 2 (A) or 0.1 h−1 (B). The square symbols are the evaluation points used for computation of dC3comp/dt by the Newton-Raphson method; these same points also appeared in the right panels of (A) and (B). The squares are the evaluation points, and the method is similar to that described for the two-compartment model cases (see Figs. 1 and 2).

Extension of the Newton-Raphson Method for Estimating Tmax in Multiple Dosing

The method may be readily applied to multiple dosing at steady-state, as exemplified below for the two-compartment model, where u is the time elapsed after the last dose here. The domain of the numerical problem is 0 ≤ u ≤ τ, where τ is the interval of administration.Embedded Image(9)Embedded Image occurs at Embedded Image, the true Tmax,2comp or Tmax,2comp,true after the last dose can be computed by solving Embedded Image.Embedded Image(10)Here, the value of Embedded Image which satisfies the condition for eq. 11 (shown below), occurs at umax,SS, the time elapsed since last dose where maximum concentration (Cmax,2comp,SS) occurs, in the steady state.

This method may be extended to multicompartment multiple dosing systems whose steady-state concentration with dosing interval of τ is given by:Embedded Image(11)An example of the calculation is shown for the two and three compartments in Fig. 4, and detailed parameter values are summarized in Table 1. The method for estimation of Tmax,multi,true is the same as that for all other models described above.

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

The log-linear blood/plasma concentration-versus-time plot for (A) two-compartment and (B) three-compartment, multiple dosing. The dosing interval (τ) is 6 h. All other parameters for (A) are identical to those from Fig. 1. All other parameters for (B) is identical to that shown in Fig. 3A.

Simulations

Simulations were performed to assess the properties of the proposed method to compare this to noncompartmental analysis, where Cmax and Tmax use observed data as a basis. The pharmacokinetic parameters used in the previous two-compartment model were applied here— α , β, and k21 of 0.8, 0.2, and 0.5 h−1, respectively. [FsysDosepo/V1] was set to 100 [where Dosepo is the oral dose administered; Fsys is systemic oral bioavailability, product of fraction absorbed (Fabs) and availabilities of intestine, liver, and lung from the first-pass effect]. Blood/Plasma concentration data were simulated using commercial software (MATLAB; The MathWorks Inc., Natick, MA) for the two-compartment, oral model, with elimination in central compartment. The designated hybrid (α and β) and ka and k21 values were used to compute N, L, and M with eq. 5 and then used to provide values of C2comp,cal for each of the designated sampling times. We compared the noncompartmental method via visual identification and our method under different assay errors and sampling frequencies. Simulations were repeated 500 times and accuracy and precision of the two estimated metrics (X = Cmax or Tmax), each acquired through two different methods, were evaluated by computing bias and relative S.D. %Bias, which is the deviation of the mean value (x̅) from true value (xtrue), was calculated using eq. 12.Embedded Image(12)The relative S.D. (RSD) of the repeat simulations was calculated using eq. 13, where xi is the estimated metric for each of the repeats and x̅ is the mean of the estimates and N is the number of replicate simulations (N = 500).

Embedded Image(13)

Assay Error.

To understand the effect of assay error, simulations were performed using conditions similar to those proposed by Tothfalusi and Endrenyi (2003). The sampling time was set to a geometric sequence starting at 0.05 and with a geometric ratio of 1.3 for a total of 24 time points. For the two-compartment model with elimination from the central compartment, ka, was set to be 2 h−1. Then, Gaussian noise, with a mean of 0 and relative S.D. of Embedded Image (where C2comp,cal is the error-free value and E is the assay error, 1%–15%) was generated to simulate error-containing concentrations (C2comp,error) among the 500 simulations for each of the five designated assay errors (1%, 2%, 5%, 10%, and 15%). The parameters, Tmax,2comp and Cmax,2comp, were obtained directly from N, L, and M, ka, α, and β, with the set of C2comp,error (with error) data from nonlinear fitting of the two-compartment model using MATLAB, which provided a solution for the constants, N, L, M, ka, α, and β. These fitted estimates were substituted into eq. 6 and Newton-Raphson numerical method was performed to solve the resulting expression. The noncompartmental Tmax,obs was used as the initial estimate. In cases where proper convergence was not achieved, the previous time point to Tmax,obs was used as the initial estimate and the numerical analysis was repeated. The analysis yielded Tmax,2comp and Cmax,2comp for comparison with their counterparts in the absence of assay error, Tmax,2comp,true and Cmax,2comp,true. Alternately, the parameters may be used to simulate concentrations at small incremental time intervals to arrive at Cmax estimates close to the Cmax,2comp,true. The maximum concentration, Cmax,2comp,true, and its corresponding time, Tmax,2comp,true, were then identified visually.

The %Bias and relative S.D. values of Cmax,2comp (Fig. 5, A and C) and Tmax,2comp (Fig. 5, B and D) owing to assay error from our proposed model-guided estimation method were much smaller than those from the noncompartmental method, whose basis is merely visual inspection of the observed data for Cmax,2comp,obs and Tmax,2comp,obs. Essentially, there was little bias on Cmax,2comp and Tmax,2comp for our model-guided estimation method, whereas for the noncompartmental method, the %Bias for Cmax,2comp,obs was higher and proportional to the assay error (Fig. 5A). Thus, Cmax,2comp values obtained from datasets with different assay errors cannot be compared reliably when the noncompartmental method is used, even though the compartmental model-guided values allow for comparison across different studies.

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Plots of %Bias for Cmax (A) or Tmax (B) and the relative S.D. or RSD (%) for Cmax (C) or Tmax (D) owing to assay error for the noncompartmental method (blue square and line)—determined by noncompartmental analysis by direct observation and by model-guided determination (red triangles and line) for two-compartment model.

Sampling Frequency.

Next, for evaluation of the impact of sampling the frequency, 500 simulations were performed for each for the four sampling frequencies ranging from 5, 10, 20, and 40. The sampling frequency reflects the number of sampling points between 0 and Tmax,2comp,true, allowing the true Tmax,2comp to fall on one of the sampling times. Again, the hypothetical patient data were simulated with the assigned set of ka, k21, α, and β values assuming elimination from the central compartment, a specific blood sampling scheme, and a constant 5% assay error. These data were refitted by nonlinear least squares method to obtain the parameters N, L, M, ka, α, and β, such that the Cmax,2comp,true and Tmax,2comp,true could be obtained with the Newton-Raphson method, as previously described.

Even as the sampling frequency was increased, Cmax,2comp obtained by our method (Fig. 6A) showed no apparent bias, in contrast to the FDA-recommended method, which showed positive bias. Furthermore, the degree of bias increased as sampling frequency was increased, suggesting that Cmax cannot easily be compared between datasets that were generated using different sampling frequencies. The positive bias for the visual inspection method is at first counter-intuitive since in the absence of assay errors, Cmax should have negative bias as by definition all other concentration values are equal or lower. However, owing to the presence of assay errors, the measured concentrations near the true Cmax,true may be greater than Cmax,true and since the noncompartmental method whose basis is visual inspection selects the largest concentration, positive bias occurs on average.

Fig. 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 6.

Plots of %Bias for Cmax (A) or Tmax (B) and the relative S.D. or RSD (%) for Cmax (C) or Tmax (D) owing to sampling frequency for the noncompartmental method (blue square and line)—determined by noncompartmental analysis by direct observation and by model-guided determination (red triangles and line) for two-compartment model. The sampling frequency is the number of sampling points between 0 and Tmax,2comp,true. The parameters were same as that shown in Fig. 1.

In addition, both Tmax,2comp and Cmax,2comp can be obtained with more certainty (lower relative S.D.) using our method than when obtained directly from the measurements (Fig. 6, C and D). Because population variation was not added to the data, the variation measured here arose from the assay errors and the sampling scheme. Figure 6B shows that the Tmax,2comp obtained using the two methods show only minor bias. The small positive bias shown here results from the asymmetric peak shape near Tmax,true (i.e., peak is sharper in the left region than the right region). As previously noted (Basson et al., 1996), measuring Tmax,2comp is a universal method and the Tmax,2comp may be used also under flip-flop kinetics, for which separation of ka may be difficult without intravenous data. However, the determined Tmax,2comp suffers from a higher relative S.D. compared with Cmax,2comp (Fig. 6, C and D).

ka.

The influence of ka on the Tmax,2comp and Cmax,2comp estimates and the applicability of the method were appraised using simulations with five different ka values: 3, 2, 1, 0.1, and 0.05 h−1 (Table 2). The sampling times were made between 0.025 and 48 h (at 0.025, 0.05, 0.1, 0.25, 0.75, 1, 1.5, 2, 2.5, 3, 4, 6, 8, 12, 18, 24, 30, 36, 42 and 48 h) with an assay error of 5% (i.e., assay error E = 5). The method was found to be applicable to cases where ka > β and ka < β. For all ka values, the RSD was lower for the two-compartmental fitting–based method, and varied between 2.19% and 5.78% in contrast to the noncompartmental method, whose RSD varied between 3.76% and 23.3% (Table 2).

View this table:
  • View inline
  • View popup
TABLE 2

Effects of varying ka on estimates of Tmax and Cmax from simulated blood data, which contained 5% assay error

Values assigned were: α, β, and k21 = 0.8, 0.2, and 0.5 h−1, respectively.

Discussion

Generally speaking, the model-guided estimation method proposed here using the Newton-Raphson method requires additional analysis in comparison with noncompartmental analysis. The Newton-Raphson numerical method was implemented from scratch and the estimates would converge well when the Tmax,obs was used as initial estimate. Moreover, in various mathematical software applications, zero-finding methods are available as built-in functions (for example, fzero in the case of MATLAB). Not all of these use the Newton-Raphson method, but the functions used are adequate in solving the zero-finding problem numerically (i.e., f(x) = 0 problem). Despite the small amount of added work compared with the noncompartmental analysis method, the model-guided analysis provides reduced bias for the determination of Cmax (Fig. 5A; Fig. 6A; Table 2) and is associated with a higher precision in the determination of both Cmax (Fig. 5C; Fig. 6C) and Tmax (Fig. 5D; Fig. 6D). Moreover, in this analysis, we have not considered the Cmax at fast burst effect of the drug at t = 0, when the burst behaves like an intravenous bolus, owing to the accompanying complexities.

In contrast to the Loo-Riegelman method (Loo and Riegelman, 1968), the model-guided estimation method proposed here does not assume a detailed mechanistic model nor require microconstants, making it suitable even when intravenous data are not available. In addition, our method may be used to analyze data under flip-flop kinetics; this is demonstrated by the example with ka < β (Fig. 2). In fact, our method is applicable generally to all multicompartmental models that are approximated as sum of exponentials under linear conditions. While the number of exponentials is required, explicit knowledge of the detailed compartmental scheme is not required.

One disadvantage of our method is that bias may be introduced when an inappropriate number of exponentials is selected. In cases where the constants are close in values, that is, ka ≈ α or ka ≈ β, the correct number of compartments will probably be underestimated, and Tmax,true is estimated with increased bias, albeit modest. Under the condition that ka = α or when ka = β, the accurate identification of the Tmax,true is not hampered and may be estimated without bias although we may be underestimating the number of compartments (data not shown). However, when the exponents are sufficiently separated, the appropriate number of exponentials will be easily selected, and the model-guided interpolation will provide a lower bias on the Cmax estimate and a higher precision for Cmax and Tmax. Hence, the method is viable for almost all scenarios. Furthermore, this method may be used even without intravenous data or a priori knowledge about the rank order of the hybrid constants α, β, and ka. Thus, the model-guided estimation method is appropriate for comparative analyses across aggregated datasets available from various sources, with varying collection parameters and limited mechanistic knowledge about the system.

Conclusion

A proper multiexponential fit, combined with the Newton-Raphson method, can provide a better approach for Cmax and Tmax estimates than the observational, noncompartmental analysis method recommended by FDA, and is more appropriate when limited data are available.

Authorship Contributions

Participated in research design: Han, Pang.

Conducted experiments: Han.

Contributed new reagents or analytic tools: Han.

Performed data analysis: Han.

Wrote or contributed to the writing of the manuscript: Han, Lee, Pang.

Footnotes

    • Received May 30, 2018.
    • Accepted August 16, 2018.
  • This work was supported by NSERC (to P.I.L. and Y.R.H.) and the Centre for Collaborative Research (CCDR), University of Toronto (to K.S.P.)

  • https://doi.org/10.1124/dmd.118.082636.

Abbreviations

α
hybrid constant for distribution for two-compartment model
β
hybrid constant for terminal decay for two-compartment model
λi
hybrid constant for an ith compartment model
τ
dosing interval in a multiple dosing regimen
Cmax
maximum drug concentration in blood/plasma
Cmax,obs
maximum drug concentration observed among the blood/plasma sampling points
Cmax,2comp
maximum blood/plasma drug concentration in (central) compartment from two-compartment model fit
Cmax,true
true maximum drug concentration reached in the blood/plasma (central) compartment
ka
first-order absorption rate constant
k10
first order elimination rate constant from compartment 1
k12
first order transfer rate constant from compartment 1 to 2
k21
first-order transfer rate constant from compartment 2 to 1; these constants are used in both two- and three-compartment models
RSD
relative standard deviation
Tmax,obs
time observed to reach the Cmax,obs in blood/plasma (central) compartment
Tmax,true
true time needed to reach the maximum drug concentration in the blood/plasma (central) compartment
Tmax,2comp
time observed to reach the Cmax,2comp in the blood/plasma (central) compartment on the basis of the two-compartment model
u
time elapsed since last dose in multiple dosing
V
volume of central compartment for one compartment model
V1
volume of the central compartment for multicompartment models
  • Copyright © 2018 by The American Society for Pharmacology and Experimental Therapeutics

References

  1. ↵
    1. Allard S,
    2. Sainati SM, and
    3. Roth-Schechter BF
    (1999) Coadministration of short-term zolpidem with sertraline in healthy women. J Clin Pharmacol 39:184–191.
    OpenUrlPubMed
  2. ↵
    1. Barbhaiya RH,
    2. Shukla UA,
    3. Gleason CR,
    4. Shyu WC, and
    5. Pittman KA
    (1990) Comparison of the effects of food on the pharmacokinetics of cefprozil and cefaclor. Antimicrob Agents Chemother 34:1210–1213.
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Basson RP,
    2. Cerimele BJ,
    3. DeSante KA, and
    4. Howey DC
    (1996) Tmax: an unconfounded metric for rate of absorption in single dose bioequivalence studies. Pharm Res 13:324–328.
    OpenUrlPubMed
  4. ↵
    CDER/FDA (2015) Waiver of In Vivo Bioavailability and Bioequivalence Studies for Immediate Release Solid Oral Dosage Forms Based on a Biopharmaceutics Classification System. Guidance for Industry. Center for Drug Evaluation and Research, Food and Drug Administration, Silver Springs, MD, pp. 1–2.
  5. ↵
    1. Chan KK and
    2. Gibaldi M
    (1985) Assessment of drug absorption after oral administration. J Pharm Sci 74:388–393.
    OpenUrlCrossRefPubMed
  6. ↵
    1. Collins SL,
    2. Faura CC,
    3. Moore RA, and
    4. McQuay HJ
    (1998) Peak plasma concentrations after oral morphine: a systematic review. J Pain Symptom Manage 16:388–402.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Endrenyi L and
    2. Al-Shaikh P
    (1995) Sensitive and specific determination of the equivalence of absorption rates. Pharm Res 12:1856–1864.
    OpenUrlPubMed
  8. ↵
    1. Galantai A
    (2000) The theory of Newton’s method. J Comput Appl Math 124:25–44.
    OpenUrl
  9. ↵
    1. Gibaldi M and
    2. Perrier D
    , editors (1982) Pharmacokinetics, 2nd ed, CRC Press, Boca Raton, FL.
  10. ↵
    1. Home PD,
    2. Barriocanal L, and
    3. Lindholm A
    (1999) Comparative pharmacokinetics and pharmacodynamics of the novel rapid-acting insulin analogue, insulin aspart, in healthy volunteers. Eur J Clin Pharmacol 55:199–203.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Loo JCK and
    2. Riegelman S
    (1968) New method for calculating the intrinsic absorption rate of drugs. J Pharm Sci 57:918–928.
    OpenUrlCrossRefPubMed
  12. ↵
    1. Mimaki T
    (1998) Clinical pharmacology and therapeutic drug monitoring of zonisamide. Ther Drug Monit 20:593–597.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Stuart-Harris R,
    2. Joel SP,
    3. McDonald P,
    4. Currow D, and
    5. Slevin ML
    (2000) The pharmacokinetics of morphine and morphine glucuronide metabolites after subcutaneous bolus injection and subcutaneous infusion of morphine. Br J Clin Pharmacol 49:207–214.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Tothfalusi L and
    2. Endrenyi L
    (2003) Estimation of Cmax and Tmax in populations after single and multiple drug administrations. J Pharmacokinet Pharmacodyn 30:363–385.
    OpenUrlPubMed
  15. ↵
    1. Wagner JG and
    2. Nelson E
    (1963) Per cent absorbed time plots derived from blood level and/or urinary excretion data. J Pharm Sci 52:610–611.
    OpenUrlCrossRefPubMed
PreviousNext
Back to top

In this issue

Drug Metabolism and Disposition: 46 (11)
Drug Metabolism and Disposition
Vol. 46, Issue 11
1 Nov 2018
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Editorial Board (PDF)
  • Front Matter (PDF)
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for sharing this Drug Metabolism & Disposition article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Finding Tmax and Cmax in Multicompartmental Models
(Your Name) has forwarded a page to you from Drug Metabolism & Disposition
(Your Name) thought you would be interested in this article in Drug Metabolism & Disposition.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Research ArticleCommentary

Tmax and Cmax for Multicompartmental Models

Yi Rang Han, Ping I. Lee and K. Sandy Pang
Drug Metabolism and Disposition November 1, 2018, 46 (11) 1796-1804; DOI: https://doi.org/10.1124/dmd.118.082636

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero

Share
Research ArticleCommentary

Tmax and Cmax for Multicompartmental Models

Yi Rang Han, Ping I. Lee and K. Sandy Pang
Drug Metabolism and Disposition November 1, 2018, 46 (11) 1796-1804; DOI: https://doi.org/10.1124/dmd.118.082636
del.icio.us logo Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Theory and Methods
    • Discussion
    • Authorship Contributions
    • Footnotes
    • Abbreviations
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Related Articles

Cited By...

More in this TOC Section

  • Kcat assessment of UGT enzymes across different organs.
  • The Future of P450 CYP Induction
  • Clearance Models: The Case for Status Quo
Show more Commentary

Similar Articles

Advertisement
  • Home
  • Alerts
Facebook   Twitter   LinkedIn   RSS

Navigate

  • Current Issue
  • Fast Forward by date
  • Fast Forward by section
  • Latest Articles
  • Archive
  • Search for Articles
  • Feedback
  • ASPET

More Information

  • About DMD
  • Editorial Board
  • Instructions to Authors
  • Submit a Manuscript
  • Customized Alerts
  • RSS Feeds
  • Subscriptions
  • Permissions
  • Terms & Conditions of Use

ASPET's Other Journals

  • Journal of Pharmacology and Experimental Therapeutics
  • Molecular Pharmacology
  • Pharmacological Reviews
  • Pharmacology Research & Perspectives
ISSN 1521-009X (Online)

Copyright © 2022 by the American Society for Pharmacology and Experimental Therapeutics