Next Article in Journal
Effective Immobilization of Novel Antimicrobial Peptides via Conjugation onto Activated Silicon Catheter Surfaces
Next Article in Special Issue
Pharmacokinetic Models of Tafenoquine: Insights for Optimal Malaria Treatment Strategies
Previous Article in Journal
Drug Loading in Chitosan-Based Nanoparticles
Previous Article in Special Issue
Population Pharmacokinetic Model of Linezolid and Probability of Target Attainment in Patients with COVID-19-Associated Acute Respiratory Distress Syndrome on Veno-Venous Extracorporeal Membrane Oxygenation—A Step toward Correct Dosing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pharmacokinetic Simulation and Area under the Curve Estimation of Drugs Subject to Enterohepatic Circulation

by
Melchor Alpízar
*,
José de Jesús Reséndiz
,
Elisa García Martínez
,
Sanyog Dwivedi
and
Miguel Alejandro Trejo
*
Specialized Centre for Diabetes, Obesity and Prevention of Cardiovascular Diseases (CEDOPEC) & Bioequivalence Studies (E-BIOS-CEDOPEC), Mexico City 11650, Mexico
*
Authors to whom correspondence should be addressed.
Pharmaceutics 2024, 16(8), 1044; https://doi.org/10.3390/pharmaceutics16081044
Submission received: 6 June 2024 / Revised: 13 July 2024 / Accepted: 31 July 2024 / Published: 6 August 2024
(This article belongs to the Special Issue Optimizing Drug Safety and Efficacy: Pharmacokinetic Modeling)

Abstract

:
Enterohepatic circulation (EHC) is a complex process where drugs undergo secretion and reabsorption from the intestinal lumen multiple times, resulting in pharmacokinetic profiles with multiple peaks. The impact of EHC on area under the curve (AUC) has been a topic of extensive debate, questioning the suitability of conventional AUC estimation methods. Moreover, a universal model for accurately estimating AUC in EHC scenarios is lacking. To address this gap, we conducted a simulation study evaluating five empirical models under various sampling strategies to assess their performance in AUC estimation. Our results identify the most suitable model for EHC scenarios and underscore the critical role of meal-based sampling strategies in accurate AUC estimation. Additionally, we demonstrate that while the trapezoidal method performs comparably to other models with a large number of samples, alternative models are essential when sample numbers are limited. These findings not only illuminate how EHC influences AUC but also pave the way for the application of empirical models in real-world drug studies.

1. Introduction

Enterohepatic circulation (EHC) represents a dynamic kinetic process integral to the pharmacokinetics of certain drugs, delineating a crucial pathway within the body’s metabolic landscape [1]. This intricate cycle commences with drug absorption from the intestine, where a portion of the dissolved drug traverses through the portal system to reach the hepatic lobules. Here, within the intricate architecture of the liver, the drug navigates through branches of the hepatic artery and portal vein, eventually circulating to the central vein. Along this journey, a fraction of the drug diffuses through the sinusoidal walls, finding its way into hepatocytes [2,3]. Once within these hepatic cells, the drug may undergo active transportation, facilitated by an array of pumps, such as BCRP, MATE1, MDR3, or MRP2, towards the bile canaliculi, ultimately contributing to bile formation [2,3]. Subsequently, bile, laden with the drug, travels through the bile ducts to the gallbladder, where it awaits the signal for release prompted by food intake and the hormone cholecystokinin [4,5]. The gallbladder then empties its contents, a process spanning approximately 60 min, through the sphincter of Oddi into the duodenum [6,7,8]. Within the duodenum, the drug dissolved in bile may undergo reabsorption, initiating a new cycle of circulation. Figure 1a elucidates these intricate processes alongside parallel kinetic phenomena.
In studying EHC, researchers analyze time versus plasma drug concentration curves to uncover the dynamic interplay between absorption, distribution, metabolism, and excretion, revealing peaks and troughs corresponding to distinct phases of the enterohepatic cycle. Despite offering valuable insights, such analysis is challenged by variability in physiology, drug properties, and environmental factors, along with practical considerations like sampling frequency and timing. Understanding EHC is crucial in drug studies, as it profoundly influences pharmacokinetic parameters essential for determining pharmaceutical efficacy and safety. Consequently, accurately characterizing these plasma profile fluctuations poses a significant challenge in EHC studies, necessitating the development of precise models.
The journey to model EHC began with compartmental models [9], conceptualizing the intestine and organism as discrete compartments with drug transfer occurring continuously [10] or discontinuously [11,12]. Over time, these models evolved in complexity, incorporating additional compartments to account for the gallbladder. Dynamics governing the release of bile into the intestine were then modeled with various kinetics, including continuous [13], single [14,15,16,17,18,19], or multiple, the latter encompassing models utilizing switch-type [20,21,22,23,24,25,26], sigmoid [27], or sinusoidal [28,29] functions. This evolution facilitated a deeper comprehension of EHC and its ramifications on drug pharmacokinetics.
Physiologically based pharmacokinetic models represent the pinnacle of complexity within EHC modeling, treating organs and tissues as distinct entities with specific properties for both organism and drug [30,31,32,33]. However, despite their increasing prevalence, their translation into clinical therapeutics remains limited [34,35].
Alternatively, empirical mathematical functions, devoid of physiological or compartmental underpinnings, have been employed in EHC modeling [36]. The pioneering model in EHC employed two gamma functions to fit the plasma mycophenolic acid profile, while another model utilized the sum of two Gaussian functions to depict the plasma profile of simvastatin. Demonstrated feasibility extends these models beyond two peaks by incorporating additional terms [37,38,39]. While these models may not elucidate EHC’s underlying processes, they align with experimental data and capture the secondary peak characteristics of EHC. Consequently, they have found widespread implementation in clinical settings as valuable tools in therapeutic drug monitoring (TDM).
Understanding the magnitude of EHC for drugs assumes paramount importance, as it significantly influences key pharmacokinetic parameters such as clearance ( C l ), volume of distribution ( V d ), half-life plasma ( t 1 / 2 ), and A U C [40]. Precisely, the A U C is frequently used in TDM for dose adjustment in drugs with a narrow therapeutic range [41]. Sampling times and the number of samples are a limiting factor in TDM that requires restricting the sampling strategy. On the other hand, the A U C is an essential parameter in bioequivalence studies, and it has been argued that the conventional trapezoidal method, used to calculate the A U C , probably overestimates its value when EHC exists.
The debate over whether empirical models offer advantages over traditional methods, like the trapezoidal rule, in estimating A U C remains unresolved. Thus, this study aims to evaluate the efficacy of several empirical models in characterizing time–concentration profiles of hypothetical and real drugs sharing the trait of undergoing EHC. This evaluation, grounded in simulations of various pharmacokinetic scenarios, aims to shed light on the potential of empirical models in understanding and quantifying EHC’s impact from both conventional and limited sampling perspectives.

2. Materials and Methods

A compartmental model and a set of empirical models were developed. The compartmental model was used to perform simulations of different pharmacokinetic scenarios; the empirical models were used to estimate the A U C of these simulations.

2.1. Compartmental Pharmacokinetic Model

The compartmental model was built considering the kinetic processes involved in EHC, in addition to those that can occur in parallel and that contribute to global pharmacokinetics; see Figure 1a.
The behavior of the EHC phenomenon is schematized in Figure 1b and is given the following system of differential equations:
d D C d t = k a D G + k 21 D P k 12 D C k e h c D C k r e l D C v b i o d D P d t = k 12 D C k 21 D P d D G d t = f d i s D S + k g e m D B k a D G k g e l D G d D S d t = f d i s D S d D B d t = k e h c D C k g e m D B
where D x is the concentration of drug and whose subscript indicates the site in which it is located ( D c central compartment; D P peripheral; D G intestinal; D B biliary; D S in the solid pharmaceutical form). Except for the parameters f d i s , k g e m , and v b i o , which are set as parameters as a function of time or drug concentration, most of the parameters in this system of equations are constant during the simulation.
Fraction of dissolved dose of solid dosage form ( f d i s ) is one of the time-dependent parameters and was simulated with the Weibull function in Equation (2). This is a function without a physicochemical basis widely used in modeling drug dissolution [42], and it is expressed as follows:
f d i s t = 1 e t t l a g k d a
where t l a g corresponds to the time it takes to start the dissolution process, k d is the scale parameter that determines the fraction of drug released from the formulation, and a is the shape parameter that characterizes the curve as exponential ( a = 1 ), sigmoid ( a > 1 ), or parabolic ( a < 1 ).
Emptying of the gallbladder ( k g e m ) is a time-dependent parameter; continuous function, given by Equation (3), was used to simulate emptying of the gallbladder and allows the generation of several emptying cycles distributed over a period of 24 h; it is expressed in the following way:
k g e m t = k m a x i = 1 n e b sin π t c i 24 2
where the parameter k m a x is the maximum emptying speed of the gallbladder, the parameter b is related to the duration of each of the emptying cycles, and the parameter c i allows the establishment of the moments in which the emptying of gallbladder occurs for each of the EHC cycles according to meal times (established as hours elapsed after drug intake).
The last non-constant parameter v b i o is a function of the drug concentration in the central compartment and represents the biotransformation rate of the drug. For its calculation, the Michaelis–Menten equation, Equation (4), is used. This equation allows us to consider the effects of the competition that may exist between biotransformation and other processes that the drug may have (bile transfer, transfer to the peripheral compartment, or renal elimination), and it is expressed as follows:
v b i o = v m a x D C k m + D C
where k m is the Michaelis–Menten constant, and v m a x is the maximum biotransformation rate.
To complete the analysis, the percentage of drugs within the central compartment undergoing enterohepatic circulation was calculated. This percentage is not necessarily constant over time; Equation (5) allows the calculation of its temporal fluctuation, while Equation (6) calculates its theoretical maximum value.
E H C % = k e h c k e h c + k r e l + k 12 + v b i o × 100 %
E H C m a x % = k e h c k e h c + k r e l + k 12 × 100 %

2.2. Compartmental Model Assumptions

The compartmental model was developed with the objective of representing the EHC and some parallel pharmacokinetic processes in the most realistic way possible; however, to simplify the analysis, some reasonable assumptions are made:
  • Compartments are homogeneous, so drug concentration within each of the compartments reaches an instantaneous kinetic equilibrium;
  • Liver is located within the central compartment, so there is an instantaneous kinetic equilibrium between plasma and intrahepatic drug concentrations;
  • Elimination by routes other than renal, fecal, or biotransformation is insignificant;
  • Fraction of drug that is biotransformed does not need to be reconverted to the original drug to undergo EHC, as this usually happens under the presence of microbiota glucuronidases;
  • Elimination by biotransformation occurs through a single pathway involving a single enzyme;
  • Drug dissolved in the bile is excreted to the intestinal compartment only from the bile compartment, so direct excretion from the central compartment is non-existent;
  • Fraction of drug dissolved, transfer of drug from the gallbladder to the intestine, and rate of biotransformation can be realistically represented by equations, while the rest of the parameters follow first-order kinetics that do not consider changes in transport rates between compartments due to physiological phenomena.
Finally, during simulations, it was assumed that the gallbladder is only activated by meals three times during a 24-hour period; at other times, its emptying is null.

2.3. Pharmacokinetic Simulations

The compartmental model was implemented through the SimBiology tool of MATLAB version R2021a. Drug concentrations versus time were simulated with the sundials solvertype. The description, units, and specific values of model parameters are shown in Table 1. These values were proposed to have a starting model for subsequent scenarios.
Data shown in Table 1 correspond to those of the standard model from which the values of one parameter were varied at a time to generate 18 different scenarios. These scenarios are grouped into pairs that differ in whether the value of the modified parameter is greater or less than its value in the standard model. Table 2 shows the description, the modified parameter, and its value in each scenario.
Table 1. Parameters of the standard model.
Table 1. Parameters of the standard model.
ParameterUnitsValueDescription
k a h−10.9Intestinal absorption rate
k r e l h−10.1Renal elimination rate
k g e l h−10.1Fecal elimination rate
f d i s -See Equation (2)Dissolved drug fraction
k d h−10.5Scale parameter of the Weibull function
t l a g h0Weibull function delay time
a -1Shape parameter of the Weibull function
k g e m h−1See Equation (3)Rate of transfer from gallbladder to intestine
k m a x h−13.0Maximum transfer rate from gallbladder to intestine
b -300Biliary emptying cycle duration parameter
c i -5, 11, 17Times of maximum gallbladder emptying rate
k 12 h−10Transfer rate from central to peripheral compartment
k 21 h−10Transfer rate from peripheral to central compartment
v b i o mg/hSee Equation (4)Instantaneous biotransformation rate
v m a x mg/h0.1Maximum biotransformation rate
k m mg/LSee Table 3Michaelis–Menten constant for biotransformation
k e h c h−1See Table 3Transfer rate from the central compartment to the gallbladder
Table 2. Modified parameters in pharmacokinetic scenarios.
Table 2. Modified parameters in pharmacokinetic scenarios.
ScenariosParameterValue k m  a k e h c  aDescription
S 1 k a 1.5IARapid absorption in the intestinal compartment
S 2 0.3IASlow absorption in the intestinal compartment
S 3 k r e l 0.5IIBRapid renal elimination from central compartment
S 4 0.02IIICSlow renal elimination from central compartment
S 5 k g e l 0.5IARapid elimination from intestinal compartment
S 6 0.02IASlow elimination from intestinal compartment
S 7 t l a g 0.0IADrug release without delay time
S 8 2.0IATime-delayed drug release
S 9 a 0.5IAParabolic drug release
S 10 2.0IASigmoidal drug release
S 11 k m a x 6.0IARapid gallbladder emptying
S 12 1.5IASlow gallbladder emptying
S 13 B600IAShort gallbladder emptying duration
S 14 100IALong gallbladder emptying duration
S 15 k 12 ,   k 21 1.0, 0.1IIDFast peripheral distribution with slow central return
S 16 0.1, 1.0IESlow peripheral distribution with fast central return
S 17 v m a x 1.0IVABroad first pass elimination
S 18 0.02IAReduced first pass removal
a Values of k m and k e h c are grouped in different cases represented by Roman numerals and capital letters, respectively; see Table 3.
Table 3. Values of the constants k m and k e h c and their effect on the percentage of drugs that undergo EHC.
Table 3. Values of the constants k m and k e h c and their effect on the percentage of drugs that undergo EHC.
kmEffect on
EHC% a
kehcEHCmax%
IIIIIIIVABCDE
101100100+0.06660.33330.01330.73330.133320%
0.10.01110++0.150.750.031.650.360%
0.0010.00010.011+++0.42.00.084.40.880%
a The effect is observed as a decrease with respect to the maximum value: insignificant decrease: +; moderate decrease: ++; outstanding decrease: +++.
To increase the diversity of the data generated in each scenario, the values of the parameters k m and k e h c , which are related to the magnitude of the EHC, were varied. Table 3 shows the values of k m and k e h c for each of the cases and their effect on the EHC. Different combinations of these two parameters allowed the generation of up to nine different simulations per scenario, giving a total of one hundred and sixty-two simulations. At the end of each simulation, post-simulation calculations were performed to obtain the A U C 0 .
From each simulation, time versus plasma concentration data were extracted at specific times to represent the sampling process. The sampling schemes were divided into three categories, which are shown in Table 4.

2.4. Empirical Pharmacokinetic Models

The empirical models were fitted to the data extracted from each simulation. Five models ( M 1 to M 5 ) represented by functions were built to characterize the multiple peaks present in the EHC. They were compared with each other and with respect to the reference model of extravascular administration ( M 0 ), Equation (7). Equations (8) and (9) are modifications of previously reported models [37,38], while the rest, Equations (10)–(12), are new models built intuitively. They are described as:
C M 0 t = a e b t e c t
C M 1 t = a t c e b t + i = 1 n d i e t f i 2 g i 2
C M 2 t = e a t i = 1 n b i Γ d i 0 t c i x d i 1 e x d x
C M 3 t = t e a t i = 1 n b i 1 e t c i d i , t c i 0 0 , t c i < 0
C M 4 t = e c t i = 1 n b i e 2 c t d i 1 e 2 c t d i + 1 , t d i 0 0 , t d i < 0
C M 5 t = C M 4 t + e c t j = 1 i = 1 n e c 24 j 24 j + d i d 1 ln b n b 1 + c d 1 d n d 1 d n b i e 2 c t d i 24 j 1 e 2 c t d i 24 j + 1 , t d i 24 j 0 0 , t d i 24 j < 0
where the letters a , b , c , d , etc., are the parameters of the models; t is the time and C M t the plasma concentration estimated by each model. Except for the model given by Equation (1), the models are built with summations that describe the first absorption peak (when i = 1 ) and secondary resorption peaks caused by EHC (when i > 1 ); therefore, in these models, the number of peaks is equal to n . Another characteristic that all models have in common is that they incorporate a term or factor to the left of the sum and that characterizes the elimination process.
The model described by Equation (12) is an extension of the model in Equation (11) and is intended to predict the presence of additional peaks that may not be detected after the first 24 h. This prediction is made under the assumption that the gallbladder will empty with an intensity and at a time like the first day. In this model, the intensity and moment of gallbladder emptying are given by the parameters b i and d i , respectively. The prediction after 24 h is based on the parameters that describe the largest peak (usually the first absorption peak, b i and d i ) and the last observed peak ( b n and d n ).
The A U C 0 of all models converge to a finite value when the model is fitted to time versus concentration data; its value can be calculated exactly using explicit formulas. Equations (S1)–(S6) shown in the Supplementary Materials (Supplementary Data S1) allow us to perform this calculation exactly, and they were obtained with the help of Maple software version 2018.

2.5. Compilation of Pharmacokinetic Profiles from the Literature

A bibliographic search and compilation of pharmacokinetic profiles of drugs that undergo EHC reported in different works was carried out. The search criteria were drugs administered orally, in a single dose, and with a profile with at least two clearly visible peaks. WebPlotDigitizer version 4.5 software was used to recover time and plasma concentration data from the graph reported in the study. The software allows you to select the points in the pharmacokinetic profile to estimate the values of the time–concentration data pairs. This software has been used successfully for data collection from graphs [43].

2.6. Data Fitting, Statistical Analysis, and Evaluation of Empirical Models

Non-linear regression (NLR) of the simulated data retrieved from the literature was performed using the Curve Fitting tool of MATLAB version R2021a. When necessary, restrictions were applied to the coefficients to obtain a good parametric fit. The adjusted R 2 , S S E , and R M S E error metrics were obtained directly from said tool. The script used to execute the NLR is provided in the Supplementary Materials (Supplementary Data S2).
Evaluation of models was carried out with the corrected Akaike Criterion ( A I C c ) for small sample sizes using Equation (13):
A I C c = n ln S S E n k + 2 k + 2 k k + 1 n k 1
where k is the number of model parameters, and n is the number of observations.
As part of evaluation of the models, comparisons were made between the reference A U C 0 of the simulation, with the A U C 0 calculated by the trapezoidal method and estimated by empirical models, Equations (S1)–(S6).

3. Results

3.1. Effect of EHC Degree on AUC

Initially, we aimed to verify whether the compartmental model could accurately simulate the EHC process, allowing us to investigate its impact on A U C . Simulations following single-dose extravascular administration delineated the relationship between gallbladder emptying, gallbladder drug concentrations, and plasma profile (Figure 2). Notably, spikes in plasma profile coincided with gallbladder contractions, underscoring the dynamic nature of EHC. Additionally, drugs subject to extensive EHC (depicted by the black line) exhibited prolonged residence in the organism compared to those experiencing minimal EHC (illustrated by the yellow line).
Expanding our model to encompass biotransformation and gallbladder transport enabled us to assess the influence of their interplay. Simulations with a k m = 50 revealed a consistent EHC% value (bold lines, left panel of Figure 3), suggesting negligible competition. Conversely, simulations with k m = 0.005 yielded fluctuating EHC% values, ranging between a minimum and a theoretical maximum determined by Equation (6) (pale lines, right panel of Figure 3), indicative of significant competition. In this scenario, maximum EHC% values ranged from 88% to 33%. This competitive dynamic directly impacted A U C , where cases of negligible competition yielded higher maximum values (ranging from 20.2 to 30.9 mg × h/L) compared to scenarios with relevant competition, where A U C decreased with decreasing k m values. These findings affirm the completeness and realistic simulation capabilities of our model to simulate the EHC process.

3.2. Simulation of Pharmacokinetic Scenarios and Fitting to Empirical Models

We conducted 162 simulations across eighteen scenarios to explore their impact on A U C 0 , as summarized in Table 5. Comparative analysis revealed that scenarios S1 through S 6 and S 15 exhibited significantly different A U C 0 compared to standard model S 7 . Additional details and pharmacokinetic simulation plots for each scenario are available in the Supplementary Materials (Supplementary Data S3, Figures S1–S18).
Next, we fitted six empirical models to the obtained plasma profile data. Figure 4 illustrates the fit of these models to the data from a representative simulation. Apart from model M 0 , all models successfully captured the multiple peaks of the plasma profile, exhibiting subtle visual differences.
Criteria for selecting the best model included adjusted R 2 , AICc, and accuracy in predicting A U C . Figure 5 displays scatter plots of these criteria, while Table 6 presents their means for each model. Figure 5 allows the visual identification of the result of each regression of the empirical models in a scatter plot. Models M 4 (cyan balls) and M 5 (red balls) stand out for showing a higher degree of clustering than other models, and both models have the lowest AICc values and the highest adjusted R 2 values. The gamma distribution model (pink balls) presented poorer clustering due to the greater difficulty of performing reproducible regressions since, in this model, applying restrictions for a successful NLR was less intuitive for the authors.
Model M 5 achieved the highest values for adjusted R 2 (0.979), A I C c (−58.1), and accuracy in predicting A U C (105.9). Although M 5 exhibited the smallest spread in adjusted R 2 (±0.02), it also had the highest spread in A U C prediction accuracy (±21.4). These results indicate that M 5 demonstrated the best performance characteristics and was thus chosen for subsequent analysis.

3.3. Estimation of AUC in Different Pharmacokinetic Scenarios

Applicability of model M 5 to various pharmacokinetic scenarios was examined. Results indicated that model M 5 accurately predicted A U C for most scenarios, with close to 100% accuracy (Figure 6). However, scenarios S 3 (rapid renal elimination from central compartment) and S 15 (fast peripheral distribution with slow central return) notably overestimated the actual A U C (exceeding 110% and 150%, respectively). This overestimation persisted regardless of whether the trapezoidal method or model M 5 was employed.

3.4. Estimation of AUC under Different Sampling Schemes

Next, we investigated the impact of sampling schemes on predictions of model M 5 . Firstly, we contrasted meal-based and conventional sampling methods (Figure 7). Notably, when the percentage of drugs undergoing EHC was 40%, model-predicted A U C accuracy was nearly 100% for both sampling schemes. However, at 80% EHC, accuracy decreased to approximately 50% with conventional sampling, while remaining close to 100% with meal-based sampling. Regardless of the sampling scheme or the percentage of drugs undergoing EHC, A U C estimates by the trapezoidal method closely mirrored those by the model.
Secondly, different sampling schemes for TDM were assessed (Figure 8). It was observed that with TDM1 and TDM2 schemes, model-predicted A U C accuracy ranged between 50–75%, depending on the percentage of drugs undergoing EHC. In contrast, TDM3 and TDM4 schemes notably improved accuracy to over 75%, while TDM5 showed a modest increase. Notably, A U C predictions by the model were nearly identical with TDM1/TDM2 and TDM3/TDM4 schemes, respectively. In all cases, A U C estimates by the trapezoidal method improved gradually with increased sampling times.

3.5. Application of Model M5

Fourteen average plasma profiles of drugs undergoing EHC were gathered from the literature [44,45,46,47,48,49,50,51,52,53,54,55,56], and model M 5 was applied to fit these profiles. Additional graphs provided in the Supplementary Materials illustrate the fit of model M5 to data extracted from representative pharmacokinetic profiles (Supplementary Data S4, Figures S19–S28).
The A U C was estimated using model M 5 and compared with reported values, as shown in Figure 9, obtaining a correspondence between both values close to 100%. The largest discrepancies were observed for Ezetimibe and Meloxicam, with differences of 78.4% and 124.1%, respectively, while other estimates fell within the range of 80–120%, as depicted in Table 7.

4. Discussion

4.1. Effect of EHC Degree on AUC

Simulations of our compartmental model are in agreement with other studies in which it is stated that the EHC increases the permanence of the drug in the body by increasing t 1 / 2 apparent elimination [40,57]. We want to highlight the linear nature of our compartmental model in which an increase in dose proportionally increases A U C regardless of the degree of EHC% (Supplementary Data S5, Figure S29). This is not surprising considering that drugs such as atorvastatin and mycophenolic acid that experience EHC maintain linear pharmacokinetics with respect to A U C [58,59]. However, a less understood effect of the EHC is its impact on A U C . We take advantage of the information generated from the simulations to address a question that has been discussed in other works.
It has been discussed on several occasions how A U C is altered when the drug undergoes EHC. The prevailing position is that A U C increases when EHC exists; several papers have provided evidence to support this position [60,61,62,63]. However, in other studies, it has been proposed that the degree of EHC does not have an impact on A U C [57,64]. Recently, Ibarra et al. reported that the inclusion of a hepatic compartment, where the competition between hepatobiliary secretion and biotransformation is explicit, allows us to explain within their model why the level of A U C increases when the degree of EHC is higher [40]. Similar to the model proposed by Ibarra, in our model, this competition is considered with the difference that the elimination constant by biotransformation is not of the first order but follows Michaelis–Menten kinetics.
For drugs with a very large k m , the rate of elimination by biotransformation ( v b i o ) is negligible, but when the k m rate is very small, v b i o increases, and the competition between bile secretion and biotransformation becomes important. When this happens, v m a x is reached more quickly due to the drug’s high affinity for the enzyme. As described above, the competition between the two processes causes a fluctuation of E H C % between a theoretical minimum and a theoretical maximum. The maximum can be estimated with Equation (6), while Equation (14) allows estimating the minimum possible for E H C % , which is reached if the concentration of drug in the compartment is high enough.
E H C m i n % = k e h c k e h c + k r e l + k 12 + v m a x × 100 %
Figure 3 shows that A U C of drug is lower at E H C m a x % equal to 88% compared to those with E H C m a x % equal to 33%. Within the proposed model, this is explained by the fact that biliary transport behaves as an additional elimination pathway that results in a decrease in A U C 0 by the amount of the drug extracted from the central compartment.
According to our model, we can differentiate two ways that increase A U C by altering the degree of EHC that a drug undergoes. The first is by decreasing k m of biotransforming enzyme by adding a competitor, and the second is by suppressing drug transport into the gallbladder by decreasing the value of k e h c . We translate these two possibilities into a couple of examples.
In the first case, it could happen that a drug that is subjected to EHC and to biotransformation is administered concomitantly with another drug that is also biotransformed by the same enzyme (e.g., CYP3A4). Competition of both substrates for the same enzyme will increase k m   of the drug that EHC is experiencing. As shown in Figure 3, high values of k m   increase A U C . In fact, simvastatin, a drug that undergoes EHC, and cyclosporine are biotransformed by the same enzyme CYP3A4; when administered concomitantly, the body’s exposure to simvastatin increases enough to give rise to the toxic effects of simvastatin [65].
In the second case, a drug could be presented that undergoes EHC in the normal population, but where a mutation in the gene encoding a protein responsible for drug transport to the biliary canniculus (e.g., MRP2) is able to decrease the activity of the transporter in carriers of the mutation. According to our model, this mutation would imply a decrease in E H C m a x % and an increase in A U C with respect to the normal population. In fact, people who suffer from the so-called Dubin–Johnson syndrome have higher levels of bilirubin in plasma than a normal person because the MRP2 transporter is defective, preventing bilirubin from experiencing enterohepatic circulation to which it is normally subjected, which causes its accumulation within the body [66].
In summary, suppressing biliary transport or favoring biliary transport over biotransformation could increase A U C of drugs from two different perspectives. In one study, both perspectives were investigated and reached conclusions similar to those of our model [67], but more studies are still needed to investigate how much of the A U C increase observed in these cases can be attributed to the modification of the EHC.

4.2. Simulation of Pharmacokinetic Scenarios and Fit to Empirical Models

In addition to the biliary transfer rate ( k e h c ), we can point out from Table 5 that the parameters related to absorption ( k a ) , renal elimination ( k r e l ) , intestinal elimination ( k g e l ) , and transport between compartments ( k 12 y k 21 ) were the ones that most influenced the modification of A U C 0 with respect to the standard model. Regardless of the effect that these processes may have, all of them have the characteristic pharmacokinetic profile of “sawtooth”, a sufficient requirement for their adjustment to empirical models.
A qualitative analysis of Figure 4 reveals that except for the classical model of extravascular administration ( M 0 ), most models can capture the multi-peak phenomenon, although with some differences. For example, the model M 1 generated symmetric peaks, while the rest of the models generated asymmetric peaks more similar to those in the simulation. The models M 1 and M 2 presented a smooth transition from one peak to the next, while in the rest of these models, the transition was more abrupt. The models M 1 and M 2 tended to overestimate the magnitude of some of the peaks. The models M 4 and M 5 generated nearly identical predictions in the first 24 h, but M 5 stood out for the ability to predict subsequent peaks even if they were not identified by the sampling method.
On the other hand, a quantitative analysis of the results in Table 6 shows that increasing the number of model parameters does not necessarily imply an adjusted R 2 value closer to 1. Models that required more parameters ( M 1 , M 2 , and M 3 ) had lower R 2 values than models that required fewer parameters ( M 4 and M 5 ). However, the number of parameters increases if you increase the value of A I C c . In fact, model M 1 , which requires up to 14 parameters to fit a profile with four peaks, had the highest A I C c value of all the models evaluated. As complex as it may seem, model M 5 , described by Equation (12), requires relatively few parameters to describe various peaks, which is reflected in the smallest A I C c of all the models evaluated.
It is important to mention that the greater the number of secondary peaks that you want to model, the greater the number of parameters. Increasing the number of parameters also implies increasing the number of samples; otherwise, the NLR could not be performed. The number of samples must be greater than or equal to the number of parameters to be used to characterize one or more peaks. Table 8 shows the number of parameters (and, therefore, minimum samples required) to characterize the peaks of a plasma profile with EHC. Depending on the number of peaks observed, you might choose to use more or fewer parameters for a given model.
Regarding the accuracy of the prediction A U C , it is interesting to note that even the classical model of extravascular administration ( M 0 ) that is unable to capture secondary peaks has an accuracy in A U C estimation similar to some of the models that can capture secondary peaks. Previously, we presented an empirical model to capture the secondary peaks of the plasma profile of simvastatin [38], which is equivalent to the model M 1 of this work. Unlike our previous work, this time, the biexponential model M 0 was better at estimating A U C than the model M 1 . This suggests that model M 1 might fit well with the plasma profiles of simvastatin but not with other drugs with different pharmacokinetic characteristics.
With respect to the model first proposed by Premaud et al. [37] and that in this work, which is equivalent to the model M 2 , a couple of questions can be said. This model has undoubtedly been the most successful in practice, as it has demonstrated its superiority in the field of TDM over the fixed dose in the monitoring of mycophenonlate mofetil, a drug with a narrow therapeutic margin [68]. Adaptations to it have even been proposed for the case of profiles with three peaks, which has allowed its usefulness to be further expanded [39]. In this paper, we present the model in its general form to be able to model any number of peaks, Equation (9), and we also provide a formula for the exact calculation of A U C 0 from its parameters, Equation (S3). However, in this work, there were models that presented a better performance; this was mainly because other models required a smaller number of parameters to model the same pharmacokinetic profile. In addition, on many occasions, the execution of the NLR was less intuitive for M 2 compared to other models such as M 3 , M 4 , and M 5 .

4.3. Estimation of AUC in Different Pharmacokinetic Scenarios

Once the superiority of model M 5 had been demonstrated, we proceeded to evaluate under which conditions it could be more useful and in which it does not offer significant improvements compared to the trapezoidal method. As mentioned above, the M 5 model overestimated A U C in scenarios S 3 and S 15 . This suggests that A U C predictions from model M 5 are more likely to be inaccurate for drugs that are eliminated very quickly or those with two-compartment kinetics in which the drug remains for longer in the body. It is important to note that kinetic equilibrium is not reached instantaneously between compartments, and the addition of a peripheral compartment in scenario S 15 makes it possible to take into account the presence of highly lipophilic drugs. Considering this, model M 5 might have trouble estimating A U C of highly lipophilic drugs.
An important aspect to highlight is that A U C calculation for each scenario using Equation (S6) was not significantly different from the one calculated by the trapezoidal method. These results suggest that, if possible, the selection of a meal-based sampling scheme is more important than the selection of a complex model that fits the data well but does not offer practical improvements beyond a good fit.
In the past, it has been argued that possibly the use of the trapezoidal method in drugs undergoing EHC could lead to an overestimation of A U C . This assumption is based on the fact that the lack of resolution of the secondary peaks would lead to an overestimation of the actual area of these peaks, similar to approximating the area of a circle inscribed within a square from the area of that square. However, in this paper, we observed that, more often than not, the trapezoidal method estimates the real quite well, at least when there is an exhaustive sampling based on meals; see Figure 6.

4.4. Estimation of AUC under Different Sampling Schemes

The results suggest that it is more important to opt for a meal-based sampling method when secondary peaks are very pronounced due to a high degree of EHC than when peaks are very small due to a low degree of EHC. Even so, in practice, it would be worthwhile to opt for meal-based sampling whenever possible, as the contribution of the EHC process to the A U C of some drugs can be highly variable. For example, the variability of EHC in mycophenolic acid ranges from 10 to 60% [69].
It is important to note that in a comprehensive meal-based or conventional sampling, the A U C calculated by the trapezoidal method is very similar to that calculated by model M 5 ; see Figure 7. This result suggests that when sufficiently large numbers of samples are available, as in the case of bioequivalence studies, the differences between each of the methods for calculating A U C are negligible. In these cases, opting for the trapezoidal method would be the most practical.
When sampling was limited, we found that increasing the number of samples does not necessarily lead to an improvement in the estimation of A U C . True improvements in the estimation of A U C were observed between TDM2 to TDM3 and TDM4 to TDM5 schemes. A more in-depth analysis of Figure 8 shows that the TDM1 and TDM2 samples allow the characterization of only two peaks, TDM3 and TDM4 three peaks and TDM5 four peaks of the plasma profile. The clear conclusion is that it is more important to increase the number of properly characterized peaks than to increase the number of samples.
It is important to note that in our study, the sampling times for TDM were not chosen arbitrarily but followed the rationality of meal-based sampling. From Figure 8, it is easy to see that the sampling times coincide almost exactly with the maximums and minimums of the plasma profile. These results suggest that in the TDM of drugs subject to EHC, it is essential to establish a meal schedule for the planning of sampling and that these are taken for a sufficiently long time to characterize several peaks. In this sense, our work supports the suggestions made by Lixuan et al., who studied the effect of mealtimes on A U C of mycophenolic acid [70].

4.5. Application of Model M5

Our study has the limitation that it does not have data from individual real-world profiles; we only had access to the average plasma profiles of some drugs published in other studies. Even so, we were able to adjust model M 5 to all the profiles collected and obtained, in most cases, the estimates of A U C that were quite close to those reported in their respective articles. In the process, we obtained profiles with dynamics for which we cannot be sure that they approximate the real dynamics since several of these studies did not report or consider the mealtimes for sampling. The adjustment of the proposed model to these profiles is the first step to promote its application in data generated from real-world studies based on rational sampling.
It is clear that one of the deficiencies of model M 5 is its inability to estimate the degree of EHC; however, as seen in Table 7, this model estimates the A U C quite well. It is not easy to estimate experimentally the degree of EHC that a drug experiences, and there are few works that report this value. Due to the inability to perform an analysis of the fit of model M 5 to the profiles of the drugs studied from these data, we opted to use the elimination half-life t 1 / 2 to perform this analysis.
We previously pointed out that the M5 model had difficulties in adequately estimating the A U C in those scenarios that considered drugs with very long permanence in the body. In Table 7, we find precisely that model M 5 considerably overestimated the A U C of meloxicam and amiodarone, which have a very high t 1 / 2 (39.5 h and more than 25 days, respectively). At the other extreme, model M 5 underestimated the A U C of ezetimibe although it does not have a particularly high t 1 / 2 (between 13.8 and 15.1 h).
It should be clarified that although the empirical models were applied to data generated from simulations of the EHC process, in theory, they could be applied to any other profile that has multiple peaks caused by mechanisms other than the EHC. This would be possible since, being of an empirical nature, they do not need the determination of kinetic parameters of physiological or physicochemical basis. In this way, they could be applied to plasma profiles that show multiple peaks, such as those caused by pulsatile delivery systems [71,72] (examples of these commercially available pulsatile systems are Pulsincap®, Ritalin®, Pulsys®) or those caused by physicochemical differences between different regions of the gastrointestinal tract. It has been suggested that differences in the pH of the gastrointestinal tract that ionize mycophenolic acid and not EHC could be responsible for the double-peak phenomenon observed in this drug [73]. However, it would be inapplicable to vascular delivery drug profiles that exhibit multiple spikes regardless of the mechanism behind these spikes (e.g., fentanyl, sufentanil, propofol, and phenoperidine [74,75,76,77]).
We highlight the feasibility of applying model M 5 as a structural model in population pharmacokinetic models. Many structural models employed in population pharmacokinetic approaches have used physiological or mechanism-based concepts to describe the EHC [78]. An example of the success of this type of model contains forty-five structural parameters and simultaneously describes the extent and time course of EHC in rats, dogs, and humans for fimasartan [79]. However, empirical models with fewer parameters might present advantages as structural models. The most recent example is the non-linear mixed effect (NLME) approach based on stochastic approximation expectation–maximization for the double gamma distribution model [80].
Our model is unable to predict the effect on the pharmacokinetic profile of variables such as age, gender, body weight, pathological states, etc. In this sense, the inclusion of covariates in our model for the development of further population pharmacokinetic models would allow us to limit its application to specific clinical environments in the same way that has been done for the model proposed by Premaud et al. [37].

5. Conclusions

In this study, we conducted simulations to assess the impact of EHC on A U C in various scenarios. Our findings indicate that modulating biliary transport can effectively increase A U C . We evaluated five empirical models for characterizing EHC peaks, identifying model M 5 as the most efficient due to its simplicity, accurate A U C estimates, and ability to predict peaks beyond sampling times. While model M 5 performed well in most scenarios, it struggled with drugs exhibiting rapid renal elimination or bicompartmental kinetics.
Evaluation of different sampling schemes revealed that our proposed model did not consistently outperform the trapezoidal method, particularly with extensive sampling. However, limited sampling for TDM benefited from our model’s ability to characterize peaks, yielding superior A U C estimates compared to the trapezoidal method. In conclusion, this study underscored the importance of implementing meal-based sampling schemes and laid the groundwork for applying empirical models in real-world studies of drugs undergoing EHC.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pharmaceutics16081044/s1, Supplementary Data S1: Equations for AUC calculation; Supplementary Data S2: Source script; Supplementary Data S3: Simulations of pharmacokinetic scenarios; Supplementary Data S4: Fit of model M5 to selected plasma profiles; Supplementary Data S5: Linearity of compartmental model.

Author Contributions

Conceptualization, M.A. and M.A.T.; methodology, M.A.T.; formal analysis, M.A.T.; data curation, E.G.M. and J.d.J.R.; writing—original draft preparation, M.A.T.; writing—review and editing, S.D.; visualization, J.d.J.R.; supervision, M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Acknowledgments

The graphical abstract was created using “Server Medical Art. Servier Medical Art by Servier”, licensed under a Creative Commons Attribution 3.0 Unported License.

Conflicts of Interest

Authors were employed by the company CEDOPEC. The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Gao, Y.; Shao, J.; Jiang, Z.; Chen, J.; Gu, S.; Yu, S.; Zheng, K.; Jia, L. Drug Enterohepatic Circulation and Disposition: Constituents of Systems Pharmacokinetics. Drug Discov. Today 2014, 19, 326–340. [Google Scholar] [CrossRef]
  2. Jetter, A.; Kullak-Ublick, G.A. Drugs and Hepatic Transporters: A Review. Pharmacol. Res. 2020, 154, 104234. [Google Scholar] [CrossRef]
  3. Durník, R.; Šindlerová, L.; Babica, P.; Jurček, O. Bile Acids Transporters of Enterohepatic Circulation for Targeted Drug Delivery. Molecules 2022, 27, 2961. [Google Scholar] [CrossRef]
  4. Guiastrennec, B.; Sonne, D.; Hansen, M.; Bagger, J.; Lund, A.; Rehfeld, J.; Alskär, O.; Karlsson, M.; Vilsbøll, T.; Knop, F.; et al. Mechanism-Based Modeling of Gastric Emptying Rate and Gallbladder Emptying in Response to Caloric Intake. CPT Pharmacomet. Syst. Pharmacol. 2016, 5, 692–700. [Google Scholar] [CrossRef]
  5. Ivy, A.C. A Hormone Mechanism for Gall-Bladder Contraction & Evacuation. Am. J. Surg. 1929, 7, 455–459. [Google Scholar] [CrossRef]
  6. Shaffer, E.A.; McOrmond, P.; Duggan, H. Quantitative Cholescintigraphy: Assessment of Gallbladder Filling and Emptying and Duodenogastric Reflux. Gastroenterology 1980, 79, 899–906. [Google Scholar] [CrossRef]
  7. Portincasa, P.; Moschetta, A.; Giampaolo, M.; Palasciano, G. Diffuse Gastrointestinal Dysmotility by Ultrasonography, Manometry and Breath Tests in Colonic Inertia. Eur. Rev. Med. Pharmacol. Sci. 2000, 4, 81–87. [Google Scholar]
  8. Howard, P.J.; Murphy, G.M.; Dowling, R.H. Gall Bladder Emptying Patterns in Response to a Normal Meal in Healthy Subjects and Patients with Gall Stones: Ultrasound Study. Gut 1991, 32, 1406–1411. [Google Scholar] [CrossRef]
  9. Okour, M.; Brundage, R.C. Modeling Enterohepatic Circulation. Curr. Pharmacol. Rep. 2017, 3, 301–313. [Google Scholar] [CrossRef]
  10. Harrison, L.I.; Gibaldi, M. Influence of Cholestasis on Drug Elimination: Pharmacokinetics. J. Pharm. Sci. 1976, 65, 1346–1348. [Google Scholar] [CrossRef]
  11. Steimer, J.-L.; Plusquellec, Y.; Guillaume, A.; Boisvieux, J.-F. A Time-Lag Model for Pharmacokinetics of Drugs Subject to Enterohepatic Circulation. J. Pharm. Sci. 1982, 71, 297–302. [Google Scholar] [CrossRef]
  12. Colburn, W.A.; Hirom, P.C.; Parker, R.J.; Milburn, P. A Pharmacokinetic Model for Enterohepatic Recirculation in the Rat: Phenolphthalein, a Model Drug. Drug Metab. Dispos. 1979, 7, 100–102. [Google Scholar]
  13. Colburn, W.A. Pharmacokinetic and Biopharmaceutic Parameters During Enterohepatic Circulation of Drugs. J. Pharm. Sci. 1982, 71, 131–133. [Google Scholar] [CrossRef]
  14. Pedersen, P.V.; Miller, R. Pharmacokinetics and Bioavailability of Cimetidine in Humans. J. Pharm. Sci. 1980, 69, 394–398. [Google Scholar] [CrossRef]
  15. Pedersen, P.V.; Miller, R. Pharmacokinetics of Doxycycline Reabsorption. J. Pharm. Sci. 1980, 69, 204–207. [Google Scholar] [CrossRef]
  16. Miller, R. Pharmacokinetics and Bioavailability of Ranitidine in Humans. J. Pharm. Sci. 1984, 73, 1376–1379. [Google Scholar] [CrossRef]
  17. Ide, T.; Sasaki, T.; Maeda, K.; Higuchi, S.; Sugiyama, Y.; Ieiri, I. Quantitative Population Pharmacokinetic Analysis of Pravastatin Using an Enterohepatic Circulation Model Combined with Pharmacogenomic Information on SLCO1B1 and ABCC2 Polymorphisms. J. Clin. Pharmacol. 2009, 49, 1309–1317. [Google Scholar] [CrossRef]
  18. Funaki, T. Enterohepatic Circulation Model for Population Pharmacokinetic Analysis. J. Pharm. Pharmacol. 2010, 51, 1143–1148. [Google Scholar] [CrossRef]
  19. Strandgården, K.; Höglund, P.; Grönquist, L.; Svensson, L.; Gunnarsson, P.O. Absorption and Disposition Including Enterohepatic Circulation of (14C) Roquinimex after Oral Administration to Healthy Volunteers. Biopharm. Drug Dispos. 2000, 21, 53–67. [Google Scholar] [CrossRef]
  20. Rosner, G.; Panetta, J.; Innocenti, F.; Ratain, M. Pharmacogenetic Pathway Analysis of Irinotecan. Clin. Pharmacol. Ther. 2008, 84, 393–402. [Google Scholar] [CrossRef]
  21. Jiao, Z.; Ding, J.; Shen, J.; Liang, H.; Zhong, L.; Wang, Y.; Zhong, M.; Lu, W. Population Pharmacokinetic Modelling for Enterohepatic Circulation of Mycophenolic Acid in Healthy Chinese and the Influence of Polymorphisms in UGT1A9. Br. J. Clin. Pharmacol. 2008, 65, 893–907. [Google Scholar] [CrossRef]
  22. Sam, W.-J.; Akhlaghi, F.; Rosenbaum, S.E. Population Pharmacokinetics of Mycophenolic Acid and Its 2 Glucuronidated Metabolites in Kidney Transplant Recipients. J. Clin. Pharmacol. 2009, 49, 185–195. [Google Scholar] [CrossRef]
  23. de Winter, B.C.M.; van Gelder, T.; Sombogaard, F.; Shaw, L.M.; van Hest, R.M.; Mathot, R.A.A. Pharmacokinetic Role of Protein Binding of Mycophenolic Acid and Its Glucuronide Metabolite in Renal Transplant Recipients. J. Pharmacokinet. Pharmacodyn. 2009, 36, 541–564. [Google Scholar] [CrossRef]
  24. Yau, W.-P.; Vathsala, A.; Lou, H.-X.; Zhou, S.; Chan, E. Mechanism-Based Enterohepatic Circulation Model of Mycophenolic Acid and Its Glucuronide Metabolite: Assessment of Impact of Cyclosporine Dose in Asian Renal Transplant Patients. J. Clin. Pharmacol. 2009, 49, 684–699. [Google Scholar] [CrossRef]
  25. Berg, A.K.; Mandrekar, S.J.; Ziegler, K.L.A.; Carlson, E.C.; Szabo, E.; Ames, M.M.; Boring, D.; Limburg, P.J.; Reid, J.M. Population Pharmacokinetic Model for Cancer Chemoprevention with Sulindac in Healthy Subjects. J. Clin. Pharmacol. 2013, 53, 403–412. [Google Scholar] [CrossRef]
  26. Shepard, T.A.; Reuning, R.H.; Aarons, L.J. Estimation of Area under the Curve for Drugs Subject to Enterohepatic Cycling. J. Pharmacokinet. Biopharm. 1985, 13, 589–608. [Google Scholar] [CrossRef]
  27. Jain, L.; Woo, S.; Gardner, E.R.; Dahut, W.L.; Kohn, E.C.; Kummar, S.; Mould, D.R.; Giaccone, G.; Yarchoan, R.; Venitz, J.; et al. Population Pharmacokinetic Analysis of Sorafenib in Patients with Solid Tumours. Br. J. Clin. Pharmacol. 2011, 72, 294–305. [Google Scholar] [CrossRef]
  28. Huntjens, D.R.H.; Strougo, A.; Chain, A.; Metcalf, A.; Summerfield, S.; Spalding, D.J.M.; Danhof, M.; Della Pasqua, O. Population Pharmacokinetic Modelling of the Enterohepatic Recirculation of Diclofenac and Rofecoxib in Rats. Br. J. Pharmacol. 2008, 153, 1072–1084. [Google Scholar] [CrossRef]
  29. Wajima, T.; Yano, Y.; Oguma, T. A Pharmacokinetic Model for Analysis of Drug Disposition Profiles Undergoing Enterohepatic Circulation. J. Pharm. Pharmacol. 2010, 54, 929–934. [Google Scholar] [CrossRef]
  30. Edginton, A.N.; Zimmerman, E.I.; Vasilyeva, A.; Baker, S.D.; Panetta, J.C. Sorafenib Metabolism, Transport, and Enterohepatic Recycling: Physiologically Based Modeling and Simulation in Mice. Cancer Chemother. Pharmacol. 2016, 77, 1039–1052. [Google Scholar] [CrossRef]
  31. Abbiati, R.A.; Manca, D. Enterohepatic Circulation Effect in Physiologically Based Pharmacokinetic Models: The Sorafenib Case. Ind. Eng. Chem. Res. 2017, 56, 3156–3166. [Google Scholar] [CrossRef]
  32. Voronova, V.; Sokolov, V.; Al-Khaifi, A.; Straniero, S.; Kumar, C.; Peskov, K.; Helmlinger, G.; Rudling, M.; Angelin, B. A Physiology-Based Model of Bile Acid Distribution and Metabolism Under Healthy and Pathologic Conditions in Human Beings. Cell Mol. Gastroenterol. Hepatol. 2020, 10, 149–170. [Google Scholar] [CrossRef]
  33. Gerner, B.; Scherf-Clavel, O. Physiologically Based Pharmacokinetic Modelling of Cabozantinib to Simulate Enterohepatic Recirculation, Drug–Drug Interaction with Rifampin and Liver Impairment. Pharmaceutics 2021, 13, 778. [Google Scholar] [CrossRef]
  34. Perry, C.; Davis, G.; Conner, T.M.; Zhang, T. Utilization of Physiologically Based Pharmacokinetic Modeling in Clinical Pharmacology and Therapeutics: An Overview. Curr. Pharmacol. Rep. 2020, 6, 71–84. [Google Scholar] [CrossRef]
  35. Tan, Y.-M.; Worley, R.R.; Leonard, J.A.; Fisher, J.W. Challenges Associated with Applying Physiologically Based Pharmacokinetic Modeling for Public Health Decision-Making. Toxicol. Sci. 2018, 162, 341–348. [Google Scholar] [CrossRef]
  36. Weiss, M. Empirical Models for Fitting of Oral Concentration Time Curves with and without an Intravenous Reference. J. Pharmacokinet. Pharmacodyn. 2017, 44, 193–201. [Google Scholar] [CrossRef]
  37. Prémaud, A.; Debord, J.; Rousseau, A.; Le Meur, Y.; Toupance, O.; Lebranchu, Y.; Hoizey, G.; Le Guellec, C.; Marquet, P. A Double Absorption-Phase Model Adequately Describes Mycophenolic Acid Plasma Profiles in de Novo Renal Transplant Recipients given Oral Mycophenolate Mofetil. Clin. Pharmacokinet. 2005, 44, 837–847. [Google Scholar] [CrossRef]
  38. Alpízar-Salazar, M.; Alejandro Trejo-Rangel, M.; de Jesús Reséndiz-Rojas, J.; Frydman, T.D.; Ramos-Mundo, C. A Modified Compartmental Pharmacokinetic Model of Enterohepatic Circulation for Simvastatin in Healthy Mexican Subjects. A Pilot Study. J. Pharm. Pharmacol. Res. 2020, 4. [Google Scholar] [CrossRef]
  39. Zahr, N.; Amoura, Z.; Debord, J.; Hulot, J.-S.; Saint-Marcoux, F.; Marquet, P.; Piette, J.C.; Lechat, P. Pharmacokinetic Study of Mycophenolate Mofetil in Patients with Systemic Lupus Erythematosus and Design of Bayesian Estimator Using Limited Sampling Strategies. Clin. Pharmacokinet. 2008, 47, 277–284. [Google Scholar] [CrossRef]
  40. Ibarra, M.; Trocóniz, I.F.; Fagiolino, P. Enteric Reabsorption Processes and Their Impact on Drug Pharmacokinetics. Sci. Rep. 2021, 11, 5794. [Google Scholar] [CrossRef]
  41. van Gelder, T.; Le Meur, Y.; Shaw, L.M.; Oellerich, M.; DeNofrio, D.; Holt, C.; Holt, D.W.; Kaplan, B.; Kuypers, D.; Meiser, B.; et al. Therapeutic Drug Monitoring of Mycophenolate Mofetil in Transplantation. Ther. Drug Monit. 2006, 28, 145–154. [Google Scholar] [CrossRef]
  42. Costa, P.; Sousa Lobo, J.M. Modeling and Comparison of Dissolution Profiles. Eur. J. Pharm. Sci. 2001, 13, 123–133. [Google Scholar] [CrossRef]
  43. Drevon, D.; Fursa, S.R.; Malcolm, A.L. Intercoder Reliability and Validity of WebPlotDigitizer in Extracting Graphed Data. Behav. Modif. 2017, 41, 323–339. [Google Scholar] [CrossRef]
  44. Bae, J.-W.; Choi, C.-I.; Jang, C.-G.; Lee, S.-Y. Effects of CYP2C9*1/*13 on the Pharmacokinetics and Pharmacodynamics of Meloxicam. Br. J. Clin. Pharmacol. 2011, 71, 550–555. [Google Scholar] [CrossRef]
  45. Huang, F.; Ajavon, A.; Huang, E.; Lettieri, J.; Liu, R.; Peña, C.; Berse, M. No Effect of Levothyroxine and Levothyroxine-Induced Subclinical Thyrotoxicosis on the Pharmacokinetics of Sorafenib in Healthy Male Subjects. Thyroid 2017, 27, 1118–1127. [Google Scholar] [CrossRef]
  46. Filho, H.; Ilha, J.; Silva, L.; Borges, A.; Mendes, G.; Nucci, G. Comparative Bioavailability Study with Two Amiodarone Tablet Formulations Administered with and without Food in Healthy Subjects. Arzneimittelforschung 2011, 57, 582–590. [Google Scholar] [CrossRef]
  47. Danafar, H.; Hamidi, M. Pharmacokinetics and Bioequivalence of Methotrexate in Human Plasma Studied by Liquid Chromatography-Mass Spectrometry (LC-MS). Jundishapur J. Nat. Pharm. Prod. 2016, 11. [Google Scholar] [CrossRef]
  48. Ducharme, M.P.; Warbasse, L.H.; Edwards, D.J. Disposition of Intravenous and Oral Cyclosporine after Administration with Grapefruit Juice*. Clin. Pharmacol. Ther. 1995, 57, 485–491. [Google Scholar] [CrossRef]
  49. Madan, S.; Kumar, S.; Segal, J. Comparative Pharmacokinetic Profiles of a Novel Low-Dose Micronized-Isotretinoin 32 Mg Formulation and Lidose-Isotretinoin 40 Mg in Fed and Fasted Conditions: Two Open-Label, Randomized, Crossover Studies in Healthy Adult Participants. Acta Derm. Venereol. 2019, 100, adv00049. [Google Scholar] [CrossRef]
  50. Patel, D.S.; Sharma, N.; Patel, M.C.; Patel, B.N.; Shrivastav, P.S.; Sanyal, M. Liquid Chromatography Tandem Mass Spectrometry Method for Determination of Febuxostat in Human Plasma to Support a Bioequivalence Study. J. Adv. Pharm. Sci. Technol. 2013, 1, 37–50. [Google Scholar] [CrossRef]
  51. DeGregorio, M.W.; Wurz, G.T.; Taras, T.L.; Erkkola, R.U.; Halonen, K.H.; Huupponen, R.K. Pharmacokinetics of (Deaminohydroxy)Toremifene in Humans: A New, Selective Estrogen-Receptor Modulator. Eur. J. Clin. Pharmacol. 2000, 56, 469–475. [Google Scholar] [CrossRef] [PubMed]
  52. Patiño-Rodríguez, O.; Torres-Roque, I.; Martínez-Delgado, M.; Escobedo-Moratilla, A.; Pérez-Urizar, J. Pharmacokinetic Non-Interaction Analysis in a Fixed-Dose Formulation in Combination of Atorvastatin and Ezetimibe. Front. Pharmacol. 2014, 5, 261. [Google Scholar] [CrossRef]
  53. Palmer, J.L.; Kunhihitlu, A.; Costantini, A.; Esquivel, F.; Roush, J.; Edwards, K.; Hill, T.W.K. Pharmacokinetic Bioequivalence Crossover Study of Branded Generic and Innovator Formulations of the Cholesterol Lowering Agent Ezetimibe. Clin. Pharmacol. Drug Dev. 2014, 3, 242–248. [Google Scholar] [CrossRef]
  54. Olugemo, K.; Solorio, D.; Sheridan, C.; Young, C.L. Pharmacokinetics and Safety of Low-Dose Submicron Indomethacin 20 and 40 Mg Compared with Indomethacin 50 Mg. Postgrad. Med. 2015, 127, 223–231. [Google Scholar] [CrossRef]
  55. Limopasmanee, W.; Chansakaow, S.; Rojanasthien, N.; Manorot, M.; Sangdee, C.; Teekachunhatean, S. Effects of the Chinese Herbal Formulation (Liu Wei Di Huang Wan) on the Pharmacokinetics of Isoflavones in Postmenopausal Women. Biomed. Res. Int. 2015, 2015, 902702. [Google Scholar] [CrossRef]
  56. Yang, E.; Lee, S.; Lee, H.; Hwang, I.; Jang, I.-J.; Yu, K.-S.; Lee, S. Pharmacokinetic Comparison between Fixed-Dose Combination of Fimasartan/Amlodipine 60/10 Mg and the Corresponding Loose Combination through Partial Replicated Crossover Study in Healthy Subjects. Transl. Clin. Pharmacol. 2019, 27, 134. [Google Scholar] [CrossRef]
  57. Okour, M.; Brundage, R.C. A Gallbladder-Based Enterohepatic Circulation Model for Pharmacokinetic Studies. Eur. J. Drug Metab. Pharmacokinet. 2019, 44, 493–504. [Google Scholar] [CrossRef]
  58. Gandelman, K.; Malhotra, B.; LaBadie, R.R.; Crownover, P.; Bergstrom, T. Analytes of Interest and Choice of Dose: Two Important Considerations in the Design of Bioequivalence Studies with Atorvastatin. J. Bioequivalence Bioavailab. 2011, 3, 4. [Google Scholar] [CrossRef]
  59. Granger, D.K. Enteric-Coated Mycophenolate Sodium: Results of Two Pivotal Global Multicenter Trials. Transplant. Proc. 2001, 33, 3241–3244. [Google Scholar] [CrossRef]
  60. Ichikawa, T.; Ishida, S.; Sakiya, Y.; Sawada, Y.; Hanano, M. Biliary Excretion and Enterohepatic Cycling of Glycyrrhizin in Rats. J. Pharm. Sci. 1986, 75, 672–675. [Google Scholar] [CrossRef]
  61. Shou, M.; Lu, W.; Kari, P.H.; Xiang, C.; Liang, Y.; Lu, P.; Cui, D.; Emary, W.B.; Michel, K.B.; Adelsberger, J.K.; et al. Population Pharmacokinetic Modeling for Enterohepatic Recirculation in Rhesus Monkey. Eur. J. Pharm. Sci. 2005, 26, 151–161. [Google Scholar] [CrossRef]
  62. Zhong, Z.-Y.; Sun, B.-B.; Shu, N.; Xie, Q.-S.; Tang, X.-G.; Ling, Z.-L.; Wang, F.; Zhao, K.-J.; Xu, P.; Zhang, M.; et al. Ciprofloxacin Blocked Enterohepatic Circulation of Diclofenac and Alleviated NSAID-Induced Enteropathy in Rats Partly by Inhibiting Intestinal β-Glucuronidase Activity. Acta Pharmacol. Sin. 2016, 37, 1002–1012. [Google Scholar] [CrossRef]
  63. Meaney, C.J.; Sudchada, P.; Consiglio, J.D.; Wilding, G.E.; Cooper, L.M.; Venuto, R.C.; Tornatore, K.M. Influence of Calcineurin Inhibitor and Sex on Mycophenolic Acid Pharmacokinetics and Adverse Effects Post-Renal Transplant. J. Clin. Pharmacol. 2019, 59, 1351–1365. [Google Scholar] [CrossRef] [PubMed]
  64. Shepard, T.A.; Reuning, R.H.; Aarons, L.J. Interpretation of Area Under the Curve Measurements for Drugs Subject to Enterohepatic Cycling. J. Pharm. Sci. 1985, 74, 227–228. [Google Scholar] [CrossRef]
  65. Scarfia, R.V.; Clementi, A.; Granata, A. Rhabdomyolysis and Acute Kidney Injury Secondary to Interaction between Simvastatin and Cyclosporine. Ren. Fail. 2013, 35, 1056–1057. [Google Scholar] [CrossRef] [PubMed]
  66. Keitel, V.; Nies, A.T.; Brom, M.; Hummel-Eisenbeiss, J.; Spring, H.; Keppler, D. A Common Dubin-Johnson Syndrome Mutation Impairs Protein Maturation and Transport Activity of MRP2 (ABCC2). Am. J. Physiol. -Gastrointest. Liver Physiol. 2003, 284, G165–G174. [Google Scholar] [CrossRef]
  67. Dong, J.; Yu, X.; Wang, L.; Sun, Y.; Chen, X.; Wang, G. Effects of Cyclosporin A and Itraconazole on the Pharmacokinetics of Atorvastatin in Rats. Acta Pharmacol. Sin. 2008, 29, 1247–1252. [Google Scholar] [CrossRef]
  68. Le Meur, Y.; Büchler, M.; Thierry, A.; Caillard, S.; Villemain, F.; Lavaud, S.; Etienne, I.; Westeel, P.-F.; Hurault de Ligny, B.; Rostaing, L.; et al. Individualized Mycophenolate Mofetil Dosing Based on Drug Exposure Significantly Improves Patient Outcomes after Renal Transplantation. Am. J. Transpl. 2007, 7, 2496–2503. [Google Scholar] [CrossRef] [PubMed]
  69. Bullingham, R.E.; Nicholls, A.J.; Kamm, B.R. Clinical Pharmacokinetics of Mycophenolate Mofetil. Clin. Pharmacokinet. 1998, 34, 429–455. [Google Scholar] [CrossRef]
  70. Qian, L.; Jiao, Z.; Zhong, M. Effect of Meal Timings and Meal Content on the AUC0-12h of Mycophenolic Acid: A Simulation Study. Clin. Pharmacol. Drug Dev. 2022, 11, 1331–1340. [Google Scholar] [CrossRef]
  71. Li, J.; Chai, H.; Li, Y.; Chai, X.; Zhao, Y.; Zhao, Y.; Tao, T.; Xiang, X. A Three-Pulse Release Tablet for Amoxicillin: Preparation, Pharmacokinetic Study and Physiologically Based Pharmacokinetic Modeling. PLoS ONE 2016, 11, e0160260. [Google Scholar] [CrossRef] [PubMed]
  72. Adjei, A.L.; Chaudhary, I.; Kollins, S.H.; Padilla, A. A Pharmacokinetic Study of Methylphenidate Hydrochloride Multilayer Extended-Release Capsules (Aptensio XR®) in Preschool-Aged Children with Attention-Deficit/Hyperactivity Disorder. Paediatr. Drugs 2020, 22, 561–570. [Google Scholar] [CrossRef] [PubMed]
  73. Woillard, J.-B.; Debord, J.; Marquet, P. Comment on “Population Pharmacokinetics of Mycophenolic Acid: An Update”. Clin. Pharmacokinet. 2018, 57, 1211–1213. [Google Scholar] [CrossRef] [PubMed]
  74. Stoeckel, H.; Hengstmann, J.H.; Schüttler, J. Pharmacokinetics of Fentanyl as a Possible Explanation for Recurrence of Respiratory Depression. Br. J. Anaesth. 1979, 51, 741–745. [Google Scholar] [CrossRef]
  75. Estebe, J.P.; Le Corre, P.; Levron, J.C.; Le Moing, J.P.; Le Naoures, A.; Ecoffey, C. Pilot Study on the Effect of Tourniquet Use on Sufentanil Pharmacokinetics. J. Clin. Anesth. 2002, 14, 578–583. [Google Scholar] [CrossRef] [PubMed]
  76. Cockshott, I.D.; Briggs, L.P.; Douglas, E.J.; White, M. Pharmacokinetics of Propofol in Female Patients. Studies Using Single Bolus Injections. Br. J. Anaesth. 1987, 59, 1103–1110. [Google Scholar] [CrossRef]
  77. Fischler, M.; Levron, J.C.; Trang, H.; Vaxelaire, J.F.; Flaisler, B.; Vourc’h, G. Pharmacokinetics of Phenoperidine in Anaesthetized Patients Undergoing General Surgery. Br. J. Anaesth. 1985, 57, 872–876. [Google Scholar] [CrossRef]
  78. Sherwin, C.M.T.; Fukuda, T.; Brunner, H.I.; Goebel, J.; Vinks, A.A. The Evolution of Population Pharmacokinetic Models to Describe the Enterohepatic Recycling of Mycophenolic Acid in Solid Organ Transplantation and Autoimmune Disease. Clin. Pharmacokinet. 2011, 50, 1–24. [Google Scholar] [CrossRef]
  79. Kim, T.H.; Shin, S.; Landersdorfer, C.B.; Chi, Y.H.; Paik, S.H.; Myung, J.; Yadav, R.; Horkovics-Kovats, S.; Bulitta, J.B.; Shin, B.S. Population Pharmacokinetic Modeling of the Enterohepatic Recirculation of Fimasartan in Rats, Dogs, and Humans. AAPS J. 2015, 17, 1210–1223. [Google Scholar] [CrossRef] [PubMed]
  80. Koloskoff, K.; Benito, S.; Chambon, L.; Dayan, F.; Marquet, P.; Jacqz-Aigrain, E.; Woillard, J.-B. Limited Sampling Strategy and Population Pharmacokinetic Model of Mycophenolic Acid in Pediatric Patients with Systemic Lupus Erythematosus: Application of a Double Gamma Absorption Model with SAEM Algorithm. Eur. J. Clin. Pharmacol. 2024, 80, 83–92. [Google Scholar] [CrossRef]
Figure 1. (a) Kinetic processes present in the EHC. Additional kinetic processes may occur, such as fecal elimination of a fraction of the drug secreted in the bile or not absorbed, renal elimination, distribution to peripheral tissues, hepatic biotransformation, and release of the pharmaceutical form for drugs for extravascular administration; (b) scheme of the compartmental model implemented for the simulation of the EHC. Rectangles with rounded corners represent the model compartments. The yellow circles represent the drug transfer parameters between the compartments, the majority being first-order constants. Green ellipses represent drug transfer parameters that are defined by an equation.
Figure 1. (a) Kinetic processes present in the EHC. Additional kinetic processes may occur, such as fecal elimination of a fraction of the drug secreted in the bile or not absorbed, renal elimination, distribution to peripheral tissues, hepatic biotransformation, and release of the pharmaceutical form for drugs for extravascular administration; (b) scheme of the compartmental model implemented for the simulation of the EHC. Rectangles with rounded corners represent the model compartments. The yellow circles represent the drug transfer parameters between the compartments, the majority being first-order constants. Green ellipses represent drug transfer parameters that are defined by an equation.
Pharmaceutics 16 01044 g001
Figure 2. Pharmacokinetic simulations at different percentages of drugs undergoing enterohepatic circulation. The plasma profile in the middle of the drug concentrations in the gallbladder and below the kinetics of gallbladder emptying are shown above, expressed as variation in the kgem according to Equation (3).
Figure 2. Pharmacokinetic simulations at different percentages of drugs undergoing enterohepatic circulation. The plasma profile in the middle of the drug concentrations in the gallbladder and below the kinetics of gallbladder emptying are shown above, expressed as variation in the kgem according to Equation (3).
Pharmaceutics 16 01044 g002
Figure 3. Effect of EHC% on cumulative AUC. On the (left), the temporal dynamics of the percentage of drugs subjected to EHC% is shown, and on the (right), the temporal increase in the cumulative AUC is shown.
Figure 3. Effect of EHC% on cumulative AUC. On the (left), the temporal dynamics of the percentage of drugs subjected to EHC% is shown, and on the (right), the temporal increase in the cumulative AUC is shown.
Pharmaceutics 16 01044 g003
Figure 4. Fitting of empirical models to data from a representative simulated pharmacokinetic profile. In this figure, the adjustment was made to the data obtained through a meal-based sampling, where up to 80% of the drug undergoes EHC.
Figure 4. Fitting of empirical models to data from a representative simulated pharmacokinetic profile. In this figure, the adjustment was made to the data obtained through a meal-based sampling, where up to 80% of the drug undergoes EHC.
Pharmaceutics 16 01044 g004
Figure 5. Scatter plots of the accuracy of AUC versus AICc and R2 of each model (M0–M5).
Figure 5. Scatter plots of the accuracy of AUC versus AICc and R2 of each model (M0–M5).
Pharmaceutics 16 01044 g005
Figure 6. Accuracy of AUC predicted by model M5 and by the trapezoidal method for each pharmacokinetic scenario. Calculation was made from data obtained from a meal-based sampling.
Figure 6. Accuracy of AUC predicted by model M5 and by the trapezoidal method for each pharmacokinetic scenario. Calculation was made from data obtained from a meal-based sampling.
Pharmaceutics 16 01044 g006
Figure 7. Fit of model M5 to drug profiles submitted to 40% (top), 60% (middle), and 80% (bottom) of EHC%. Data were adjusted using meal-based (Sampling I) and conventional (Sampling II) samplings, as shown in Table 4.
Figure 7. Fit of model M5 to drug profiles submitted to 40% (top), 60% (middle), and 80% (bottom) of EHC%. Data were adjusted using meal-based (Sampling I) and conventional (Sampling II) samplings, as shown in Table 4.
Pharmaceutics 16 01044 g007
Figure 8. Fit of model M5 to drug profiles submitted to 40% (top), 60% (middle), and 80% (bottom) of EHC%. Adjustment to the data was performed by sampling for the TDM, as shown in Table 4.
Figure 8. Fit of model M5 to drug profiles submitted to 40% (top), 60% (middle), and 80% (bottom) of EHC%. Adjustment to the data was performed by sampling for the TDM, as shown in Table 4.
Pharmaceutics 16 01044 g008
Figure 9. Comparison between reported for drugs that undergo EHC and AUC predicted by model M5.
Figure 9. Comparison between reported for drugs that undergo EHC and AUC predicted by model M5.
Pharmaceutics 16 01044 g009
Table 4. Sampling schemes.
Table 4. Sampling schemes.
Sampling TypeSampling Times (h)Number of SamplesDescription
Meal-based sampling0, 0.5, 1, 2, 3, 4, 5, 6, 10, 11, 12, 16, 17, 18, 29, 53, 7717Focuses on sampling before and after food intake
Conventional sampling0, 0.25, 0.5, 0.75, 1, 1.5, 2, 2.5, 3, 4, 5, 7, 9, 12, 24, 48, 7217Prioritize sampling at the beginning to characterize t m a x   and   C m a x
Sampling for TDMTDM10, 1, 3, 4, 65Reduce the number of samples and restrict their taking to the first hours after administration
TDM20, 1, 3, 4, 6, 106
TDM30, 1, 3, 4, 6, 10, 127
TDM40, 1, 3, 4, 6, 10, 12, 168
TDM50, 1, 3, 4, 6, 10, 12, 16, 189
Table 5. Average A U C 0 per scenario.
Table 5. Average A U C 0 per scenario.
Scenario A U C 0 Scenario A U C 0 Scenario A U C 0
S 1  b26.45 ± 2.73 S 7 23.91 ± 3.24 S 13 23.77 ± 3.31
S 2  b16.51 ± 3.86 S 8 23.86 ± 3.25 S 14 24.10 ± 3.13
S 3  b4.89 ± 0.64 S 9 23.89 ± 3.24 S 15  a,b9.74 ± 4.34
S 4  a,b87.52 ± 9.84 S 10 23.89 ± 3.23 S 16 20.54 ± 4.54
S 5  b12.86 ± 3.56 S 11 24.07 ± 3.15 S 17 23.49 ± 3.38
S 6  b29.06 ± 2.09 S 12 23.38 ± 3.51 S 18 24.64 ± 3.25
a Due to the greater persistence of drugs in the organism, simulation time was longer than 96 h in these scenarios. b Scenarios with p-values less than 0.05.
Table 6. Means of the criteria for model selection.
Table 6. Means of the criteria for model selection.
Model Adjusted   R 2 A I C c Accuracy (%)  A U C 0
M00.827 ± 0.14−51.9 ± 12.388.4 ± 18.7
M10.940 ± 0.17150.4 ± 12.675.9 ± 13.0
M20.888 ± 0.2652.1 ± 31.889.6 ± 17.3
M30.931 ± 0.2063.3 ± 25.884.2 ± 14.0
M40.943 ± 0.12−51.1 ± 19.186.4 ± 16.2
M50.979 ± 0.02−58.1 ± 14.7105.9 ± 21.4
Table 7. Results of NLR and AUC estimation of representative drugs.
Table 7. Results of NLR and AUC estimation of representative drugs.
DrugNumber of
Samples
Number of
Parameters
Accuracy (%)
A U C 0
Adjusted   R 2 Half - Life   t 1 / 2
(h)
Reference
Meloxicam119124.10.801339.5[44]
Sorafenib149110.60.879125.3[45]
Amiodarone159117.60.96895.1 a[46]
Methotrexate12995.40.96845.3[47]
Cyclosporin11993.60.94948.1[48]
Isotretinoin229104.60.959123.5[49]
Febuxostat199102.10.96522.78[50]
Toremifene109104.10.991229.7[51]
Atorvastatin17998.50.94003.9[52]
Ezetimibe25978.40.852415.1[53]
Ezetimibe25984.10.858613.8[53]
Indomethacin23992.20.97366.8[54]
Genistein11995.10.95274.1[55]
Fimasartan17990.50.75525.6[56]
a The elimination half-life of amiodarone is biphasic; the half-life of the first phase of disposition is presented here since the actual half-life is between 25 and 110 days.
Table 8. Number of parameters required to characterize a plasma profile.
Table 8. Number of parameters required to characterize a plasma profile.
ModelMinimum Number of Parameters and Samples per Peak a
1234
M03333
M1581114
M2471013
M3471013
M43579
M53579
a Refers to the minimum number of parameters required by each model to characterize 1, 2, 3, and 4 peaks of the plasma profile.
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

Alpízar, M.; de Jesús Reséndiz, J.; García Martínez, E.; Dwivedi, S.; Trejo, M.A. Pharmacokinetic Simulation and Area under the Curve Estimation of Drugs Subject to Enterohepatic Circulation. Pharmaceutics 2024, 16, 1044. https://doi.org/10.3390/pharmaceutics16081044

AMA Style

Alpízar M, de Jesús Reséndiz J, García Martínez E, Dwivedi S, Trejo MA. Pharmacokinetic Simulation and Area under the Curve Estimation of Drugs Subject to Enterohepatic Circulation. Pharmaceutics. 2024; 16(8):1044. https://doi.org/10.3390/pharmaceutics16081044

Chicago/Turabian Style

Alpízar, Melchor, José de Jesús Reséndiz, Elisa García Martínez, Sanyog Dwivedi, and Miguel Alejandro Trejo. 2024. "Pharmacokinetic Simulation and Area under the Curve Estimation of Drugs Subject to Enterohepatic Circulation" Pharmaceutics 16, no. 8: 1044. https://doi.org/10.3390/pharmaceutics16081044

APA Style

Alpízar, M., de Jesús Reséndiz, J., García Martínez, E., Dwivedi, S., & Trejo, M. A. (2024). Pharmacokinetic Simulation and Area under the Curve Estimation of Drugs Subject to Enterohepatic Circulation. Pharmaceutics, 16(8), 1044. https://doi.org/10.3390/pharmaceutics16081044

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