Assessing CYP2C19 Ontogeny in Neonates and Infants Using Physiologically Based Pharmacokinetic Models: Impact of Enzyme Maturation Versus Inhibition.

The objective of this study was to develop pediatric physiologically based pharmacokinetic (PBPK) models for pantoprazole and esomeprazole. Pediatric PBPK models were developed by Simcyp version 15 by incorporating cytochrome P450 (CYP)2C19 maturation and auto-inhibition. The predicted-to-observed pantoprazole clearance (CL) ratio ranged from 0.96-1.35 in children 1-17 years of age and 0.43-0.70 in term infants. The predicted-to-observed esomeprazole CL ratio ranged from 1.08-1.50 for children 6-17 years of age, and 0.15-0.33 for infants. The prediction was markedly improved by assuming no auto-inhibition of esomeprazole in infants in the PBPK model. Our results suggested that the CYP2C19 auto-inhibition model was appropriate for esomeprazole in adults and older children but could not be directly extended to infants. A better understanding of the complex interplay of enzyme maturation, inhibition, and compensatory mechanisms for CYP2C19 is necessary for PBPK modeling in infants. Published 2018. This article is a U.S. Government work and is in the public domain in the USA. CPT: Pharmacometrics & Systems Pharmacology published by Wiley Periodicals, Inc. on behalf of the American Society for Clinical Pharmacology and Therapeutics.

WHAT IS THE CURRENT KNOWLEDGE ON THE TOPIC?

Understanding CYP2C19 ontogeny is important for optimal prediction of pediatric safety and effectiveness of drugs that are substrates or interact with CYP2C19, but enzyme auto‐inhibition has not been investigated in neonates and infants.

WHAT QUESTION DID THIS STUDY ADDRESS?

Pantoprazole and esomeprazole share the CYP2C19 pathway but differ in their inhibition property on CYP2C19. Our PBPK study found that CL predictions for pantoprazole were within the twofold range for pediatrics, whereas the esomeprazole PBPK model significantly underpredicted CL in neonates and infants.

WHAT DOES THIS STUDY ADD TO OUR KNOWLEDGE?

The esomeprazole PBPK model with CYP2C19 auto‐inhibition could not be applied to neonates/infant. A better prediction without CYP2C19 auto‐inhibition, suggests that the interplay of CYP maturation and inhibition in the neonates and infants might be age‐dependent.

HOW MIGHT THIS CHANGE DRUG DISCOVERY, DEVELOPMENT, AND/OR THERAPEUTICS?

Application of PBPK modeling to inform drug exposure in pediatrics requires an understanding of mechanisms that alter drug CL, such as the interplay of enzyme inhibition and maturation as well as the possible compensatory pathways.

Pediatric populations undergo major growth‐related physiological changes that are known to alter drug disposition. For example, compared to adults, infant or neonate hepatic and renal clearance systems are immature, particularly in the first few weeks of life.1 Therefore, there is the need to characterize the pharmacokinetics (PKs) of drugs used in pediatric patients, especially neonates and infants, to ensure that the benefit/risk is optimized for these vulnerable patient populations. In spite of the need for pediatric PK data to inform dosing, prospective studies are difficult to perform in younger children, resulting in a scarcity of data.2, 3 Some of the challenges in conducting PK studies in pediatric patients include intense PK sample collection schemes that are typically not feasible, parental consent rates for these studies that are often low, and ethical concerns that are enhanced by the difficulty in demonstrating the potential for direct benefit for infants enrolled in clinical studies.4 To circumvent these challenges, modeling and simulation that leverages existing knowledge can be used to fill the knowledge gap.5 Among the modeling and simulation approaches, physiologically based pharmacokinetic (PBPK) modeling is emerging as a method that is particularly attractive because it can incorporate pediatric physiology that may be undergoing changes during growth and development in conjunction with enzyme/transporter ontogeny to improve the accuracy of predicting drug exposure.3, 6, 7, 8, 9, 10

Cytochrome P450 (CYP)2C19 is an important drug metabolizing enzyme, and knowledge of the influence of CYP2C19 is critical in pediatric pharmacology for understanding drug disposition. Despite its importance in drug metabolism, CYP2C19 expression reaches no more than 15% of mature levels throughout the prenatal period and its expression increases linearly in the first 5 postnatal months.11 In addition to CYP2C19 ontogeny, genetic polymorphisms are also important predictors of drug clearance, and must be accounted for in models used to characterize disposition of CYP2C19 substrates.12

The objective of this study was to use PBPK models to help understand CYP2C19 maturation and inhibition using two probe substrates, pantoprazole and esomeprazole, in pediatric patients through a learn‐confirm‐refine strategy. Both drugs are extensively metabolized in the liver through demethylation by CYP2C19 with subsequent sulfation and minor oxidation by CYP3A4. However, esomeprazole is an inhibitor of CYP2C19, whereas pantoprazole is not.13, 14 Once these models are developed, a similar approach could be used to predict the exposure of other drugs metabolized by CYP2C19 in infants.

METHODS

The exploration of CYP2C19 maturation followed a learn‐confirm‐refine strategy for each PBPK model. The workflow for pantoprazole and esomeprazole PBPK models used for pediatrics is described in Figure 1.15, 16 All the PBPK models were developed in Simcyp version 15 (Certara).

Figure 1

Modeling flowchart for pantoprazole and esomeprazole pediatric physiologically based pharmacokinetic (PBPK) model.8, 15, 16 CYP, cytochrome P450; PK, pharmacokinetic.

Development of adult PBPK models

For pantoprazole, first, an adult model was developed and the intrinsic clearance of CYP2C19 and CYP3A4 was optimized, as shown below, by using the published plasma concentration vs. time data collected following a 20 mg i.v. infusion or a 40 mg single oral dose of pantoprazole.17 The adult model was verified by simulating the pantoprazole plasma concentration vs. time data for different CYP2C19 genotypes and comparing them to the observed data.18 Pantoprazole area under the plasma concentration‐time curve (AUC), maximal concentration (Cmax), and clearance (CL) values for each subject were estimated by noncompartmental analysis (Phoenix WinNonlin 6.4). Table 1 lists the drug‐dependent parameters for the final pantoprazole PBPK model. The drug absorption of pantoprazole was predicted using the advanced distribution, absorption, and metabolism model and apparent permeability measured in Caco‐2 cell lines.19 The formulation parameters were obtained from the literature.20 The predicted log Pvo:w (the logarithm of the olive oil: buffer distribution coefficient at pH 7.4) was 1.326. The predicted volume of distribution at steady state (obtained using tissue volumes for a population representative of healthy volunteers) was 0.095 L/kg by the Rodgers and Rowland equation.21, 22 Renal CL of pantoprazole was 0.0012 L/hour, which was calculated as the product of plasma fraction unbound and the urine flow, which were 0.02 and 1 mL/minute, respectively.14, 23 Pantoprazole is predominantly metabolized by CYP2C19, with only a small fraction metabolized by CYP3A4.14, 24 The contribution percentage of each CYP isoform has not been previously reported. Pantoprazole's adult CL following i.v. administration was reported as ~15 L/hour.14, 25 The automatic sensitivity analysis and parameter estimation26 modules were used to estimate the intrinsic CL (CLint) of CYP2C19 and CYP3A4 by fitting against the plasma concentration vs. time data following a 20 mg i.v. dose and observed pantoprazole clearance.17 The estimated CLint of CYP2C19 was 17.6 μl/minute/pmol, whereas the CLint of CYP3A4 was 0.1996 μl/minute/pmol. With these CLint values and retrograde model tool in Simcyp, CYP2C19 contributed around 90% of pantoprazole metabolism, whereas CYP3A4 contributed the rest (~10%), which is consistent with the drug product labeling.27

Table 1

Drug‐dependent parameters for the pantoprazole and esomeprazole PBPK model

Parameters	Pantoprazole	Reference	Esomeprazole	Reference
Molecular weight (g/mol)	432.4	CheMBL	345.4	Simcyp librarya
LogPo/w	2.4	PubChem	2.23	Simcyp librarya
Compound type	Monoprotic acid	PubChem	Ampholyte	Simcyp librarya
pKa	3.92	PubChem	4.4, 8.7	Simcyp librarya
Blood/plasma	0.55	47	0.59	Simcyp librarya
Fu	0.02	14	0.03	b
Absorption model:	ADAM	19	Absorption model:
pH 6.5:7.4: Caco‐2 (10−6 cm/second)	18.3		Ka (L/hour): 10	c
Peff,man (10−4 cm/second)	3.329 (Predicted)		Fg: 1	d
			Qgut (L/hour): 6	Predicted by SimCYP
			Fu,gut: 0.03	Same as fup
Distribution model:	Full PBPK model	14	Minimal PBPK: Vss (L/kg): 0.2	48
Elimination model:	Enzyme kinetics:		Enzyme kinetics:
	CLint of CYP2C19: 17.57602 (μL/minute/pmol of isoform)	Sensitivity analysis and parameter estimate by fitting to clinically observed 20 mg i.v.	CLint of CYP2C19: 24.3 (μL/minute/pmol of isoform)	e
	CLint of CYP3A4:0.1996 (μL/minute/pmol of isoform)		CLint of CYP3A4: 0.36 μL/minute/pmol of isoform)	e
			fm of CYP2C19 (%): 73% fm of CYP3A4 (%): 27%	49
			Fu,mic: 1	f
	CLR (L/hour): 0.0012	23	CLR (L/hour):0.037	a
Enzyme interaction:	N/A		For irreversible inhibition KI of CYP2C19 (μM): 0.3 Kinact of CYP2C19 (L/hour): 5	g
			For reversible inhibition: Ki CYP2C19 (μM): 7.5	Assume Ki = IC50/2, IC50 value = 15 μM30
Formulation parameters
Solid formulation: Enteric coated tablets or granules	Triggering pH = 6.8	20	Solution
Intrinsic solubility (mg/mL)	0.05	PubChem

ADAM, advanced dissolution, absorption, and metabolism; ChEMBL, European Bioinformatics Institute; CLint, intrinsic clearance of enzyme; CLR, renal clearance; CYP, cytochrome P450; fa, fraction available from dosage form; fg, the fraction of drug that escapes first pass metabolism in the gut; fm, the relative contribution (fm) of the various elimination pathways for a drug; fu,gut, unbound fraction of drug in enterocytes fu,mic, unbound fraction in microsome; IC50, half‐maximal inhibitory concentration; ka, absorption rate constant (1/hour); Ki, concentration of inhibitor that supports half maximal inhibition (μM) for reversible inhibition; KI, concentration of inhibitor that supports half maximal inhibition (μM) for irreversible inhibition; Kinact, inactivation rate of enzyme (L/hour); LogPo/w, logarithm of the n‐octanol:buffer partition coefficient; N/A, not applicable; Qgut, a nominal flow in gut model (L/hour; Peff,man, effective permeability in man); Vss, volume of distribution at steady state (L/kg).

aAssumed same as omeprazole, obtained from Simcyp compound library. bLabel of Esomeprazole obtained from Drugs@FDA, http://www.accessdata.fda.gov/scripts/cder/drugsatfda/. cParameter estimated using Phoenix WinNonlin by compartmental analysis of phase I data50. It was assumed that 100% fraction of dose can be absorbed into enterocytes from solution. dGut metabolism is considered negligible. eRetrograde calculated value based on observed CLiv (L/hour) after 20 mg single dosing of esomeprazole49. fSimcyp compound library for omeprazole and model prediction. gSensitivity analysis and value of omeprazole29.

The esomeprazole PBPK model in the Simcyp repository was used, which was published previously.28 Table 1 provides the input parameters for the final esomeprazole PBPK model. The auto‐inhibition involved in esomeprazole clearance was modeled using both reversible and irreversible inhibition of CYP2C19 (i.e., time‐dependent inhibition (TDI)).29 However, the main contribution to auto‐inhibition for esomeprazole is from irreversible inhibition.29 The kinetic parameters describing the irreversible inhibition (i.e., TDI) are the maximal inactivation rate constant (kinact), the inhibitor concentration causing half‐maximal inactivation (KI), and the apparent first‐order degradation rate constant for the enzyme in vivo (kdeg). The software default values of kdeg,CYP2C19 are 0.0267 and 0.03/hour for the liver and gut, respectively. As a result of auto‐inhibition of CYP2C19 (via TDI), time‐variant intrinsic metabolic clearance of the drug by CPY2C19 in organs (CLuint,organ,CYP2C19) value becomes time‐dependent in both the gut and the liver. Note that in the base/initial model, time‐variant intrinsic metabolic clearance of the drug (CLuint) was obtained from retrograde, which may not be the “true” estimate of CLuint. Sensitivity analyses were conducted to explore plausible combinations of CLuint,organ,CYP2C19, KI, and kinact of CYP2C19 for esomeprazole using human PK data from various sources. Specifically, the CLuint,organ,CYP2C19 at the enzyme level was fixed at three different levels (16.2, 24.3, and 32.4 μL/minute/pmol of isoform; refer to the supporting information of the publication28) and a sensitivity analysis on KI and kinact was performed at each fixed level of intrinsic CL value. The best parameter values for CLuint,organ,CYP2C19, KI, and kinact were selected for the final model for esomeprazole by comparing the simulated PK profiles and parameters to the observed ones for i.v. and oral PK data. The kinetic parameter describing the reversible inhibition is Ki. The Ki of CYP2C19 used in the PBPK model is 7.5 μM, assumed as half‐maximal inhibitory concentration (IC50)/2 (competitive/reversible inhibition), where IC50 value = 15 μM.29

Pediatric PBPK models

The adult PBPK models for pantoprazole and esomeprazole were extended by using the age‐dependent changes of anatomic and physiological parameters (Simcyp version 15 default parameters) to predict pantoprazole exposure across different age groups in a stepwise manner (starting with adolescents and then extending to children, infants, and neonates, respectively).30, 31, 32, 33 Equations (1) and (2a), (2b) describe the enzyme ontogeny of hepatic CYP2C19 and CYP3A4 with Simcyp default parameters, respectively.34

Where Adultmax (maximum adult expression) is 0.98; F birth (fractional expression at birth relative to adult) is 0.3; Age50 (time to half adult expression) is 0.29; Age is the postnatal age in years; the “n” in Age is 2.44; and the age cutoff is 5 years (adult expression for CYP2C19 is used for children > 5 years of age).

Eq. 2a describes adult expression for CYP3A4 for age groups ≤ 25 years of age, and age is in units of years, whereas Adultmax is 1.06, F birth is 0.11, is 0.64, the “n” in Age is 1.91.(For age groups ≤ 25 years of age).34 Eq. 2b, 34 describes the adult expression for CYP3A4 for age groups > 25 years of age, and age is in units of years.

CYP2C19 and CYP3A4 are also expressed in the intestine. The ontogeny equations for intestinal CYP2C19 and CYP3A4 are the same as the respective hepatic enzyme (Eqs. (1) and (2a), (2b), respectively), but the parameter values are different. For intestinal CYP3A4 and CYP2C19, the parameters are: Adultmax, 1.059; F birth, 0.42; is 2.357; the “n” in is 1; and the age cutoff is 18 years (adult expression of CYP2C19 would be used when age is > 18 years; and Eq. 2b will be used for CYP3A4 if age is > 18 years).34

To predict the in vivo whole organ hepatic clearance of a drug metabolized by CYPs, the in vitro CLint via a CYP enzyme (e.g., obtained with in vitro liver microsomes, fresh or cryopreserved hepatocytes, or recombinant enzymes) is scaled by multiplying a series of scaling factors, including milligram of microsomal protein per gram of liver (MPPGL), CYP abundance, and total liver weight.35, 36 The CYP abundance values in various pediatric age groups are obtained using adult values that are multiplied by an enzyme‐specific hepatic ontogeny fraction obtained using the above‐described equations (Eqs. (1) and (2a), (2b)). Similarly, the intestinal CL for a drug metabolized by CYPs is obtained by scaling up CLint values for CYP enzymes by multiplying scaling factors, including microsomal protein per whole intestine, and relative CYP abundance in the intestine.37, 38 In addition to the enzyme ontogeny, age‐dependent changes in physiological parameters, such as organ size or volume, age‐related plasma protein binding are also incorporated in the model.34

PK simulations in pediatric populations

Virtual population simulations used 10 trials with 50 subjects each (500 subjects in total) for pantoprazole and 10 trials with 10 subjects each (100 subjects in total) for esomeprazole for each age group specified in the figure legends by matching the demographic data of the actual clinical study data (e.g., age, female/male ratio, etc.). The cutoff for each age band, as shown in the figure legends, was based on the available observed data.

PK simulations in different CYP2C19 genotypes

The PK parameters of pantoprazole in a CYP2C19 extensive metabolizer (EM) or in a CYP2C19 poor metabolizer (PM) were simulated and compared with clinical observations from adults (6 subjects of EM and 2 subjects of PM) and pediatric populations (21 subjects of EM and 3 subjects of PM, respectively; Table S1). The effect of CYP2C19 polymorphism on the exposure of esomeprazole in the adult population was previously assessed and published by using the same model.28

Evaluation of predictive performance

The predictive performance of each PBPK model was determined by using the ratio (R) of simulated CL (CL predicted) to the observed CL (CL observed). The 95% confidence interval (CI) of the ratio of the two means was calculated, as described by Fieller39 The R value was also calculated to evaluate the predictions of other PK parameters, such as AUC and Cmax. An R value within a range of 0.5–2.0 (twofold) was considered satisfactory.40 Furthermore, we considered that the models were acceptable when the clinical observations were between the 95th percentile and 5th percentile of the simulated mean plasma concentration‐time curve.41

Results

Adult PBPK models

The workflow of the modeling was described in Figure 1. The predicted pantoprazole plasma concentration vs. time profiles for the adult PBPK model following administration of a 20 mg i.v. infusion and a 40 mg oral delayed‐release tablet were shown in Figure 2. Clinical observations were within the within 5th and 95th percentile of the mean simulated concentration, which met one of our predefined model acceptance criteria.13, 14 The mean ratios (95% CI) of the predicted‐to‐observed (R values) CL estimates were 1.13 (1.04–1.24) and 1.03 (0.74–1.67) for the i.v. and oral data, respectively (Figure 3 a). In clinical studies of healthy adults who were administered 40 mg orally,18 the AUC, Cmax, and CL predicted by the pantoprazole PBPK model for different CYP2C19 genotypes were within twofold compared with clinical observations (Table S1). Previous studies have demonstrated that an esomeprazole adult PBPK model incorporating CYP2C19 auto‐inhibition could reasonably characterize PK profiles following single and multiple doses.28

Figure 2

Simulated pantoprazole plasma concentration‐time profile after administration of (a) 20 mg i.v. in adult, (b) single oral 40 mg delayed released tablet in adult, (c) single oral 40 mg delayed released tablet in children aged 12–16 years, (d) single oral 20 mg delayed released tablet in children aged 6–11 years, (e) i.v. infusion of 1.6 mg/kg in children aged 1–5 years, (f) single oral 1.25 mg delayed released granules in neonates. (a) Simulated vs. observed plasma time‐concentration profile of pantoprazole after i.v. infusion of 20 mg delayed release tablet (subjects = 12).17 The solid square denotes mean values from the clinical studies. The thick line represents the mean value of the simulated concentration, whereas the thin dash line represents 95th percentile and 5th percentile of simulated plasma concentration. (b) Simulated vs. observed plasma time‐concentration profile of pantoprazole after single oral administration of 40 mg delayed release tablet (subjects = 12).17 The solid square denotes mean values from the clinical studies. The thick line represents the mean values of the simulated concentration, whereas the thin dash curves represent 95th percentile and 5th percentile of simulated plasma concentration, respectively. (c) Simulated vs. observed plasma time‐concentration profile of pantoprazole after single oral administration of 40 mg delayed release tablet in children aged 12–16 years (subjects = 11).30 The solid square denotes mean values from the clinical studies. The thick line represents the mean values of the simulated concentration, whereas the thin dash curves represent 95th percentile and 5th percentile of simulated plasma concentration, respectively. (d) Simulated vs. observed plasma time‐concentration profile of pantoprazole after single oral administration of 20 mg delayed release tablet in children aged 6–11 years (subjects = 10).30 The solid square denotes mean values from the clinical studies. The thick line represents the mean values of the simulated concentration, whereas the thin dash curves represent 95th percentile and 5th percentile of simulated plasma concentration, respectively. (e) Simulated vs. observed plasma time‐concentration profile of pantoprazole after i.v. infusion of 1.6 mg/kg pantoprazole in children aged 1–5 years (subjects = 5).32 The solid squares denote the clinically observed mean plasma concentration sampled at different timepoints from the clinical studies. The thick line represents the mean values of the simulated concentration, whereas the thin dash curves represent 95th percentile and 5th percentile of simulated plasma concentration, respectively. (f) Simulated vs. observed plasma time‐concentration profile of pantoprazole after single oral administration of 1.25 mg delayed release granules in neonates (subjects = 14).33 The solid square with error bar denotes mean values (with SD) from the clinical studies. The thick line represents the mean values of the simulated concentration, whereas the thin dash curves represent 95th percentile and 5th percentile of simulated plasma concentration, respectively.

Figure 3

Comparison between pantoprazole (a) and esomeprazole (b) observed and predicted value of clearance (CL) ratio in adult and different age groups of pediatric populations. Results are presented as mean ratios (solid circles) in each age group with a 95% confidence interval (horizontal lines). The ratios in X‐axis are shown in log scale. Dashed lines represent where ratio = −0.301 (Log10 0.5) and 0.301 (Log10 2), respectively. “w/o inh” indicates simulation without cytochrome P450 2C19 inhibition.

Pantoprazole

Figure 2 shows the simulated vs. observed pantoprazole plasma concentration‐time profiles in neonates, and in children between 1 and 5 years, 6 and 11 years, and 12 and 16 years of age. No observed PK profiles were available for children 1 month to 1 year of age. The clinical observations were generally well aligned with the mean simulated plasma concentration. The model slightly underestimated drug exposure in neonates, but the clinical observations were still between the 95th percentile and 5th percentile of the mean simulated concentration. The R values for CL (95% CI) were 1.29 (0.84–1.78) for children 12–16 years of age given 40 mg pantoprazole orally,30 0.96 (0.70–1.30) for children 6–11 years given 20 mg orally,30 1.35 (1.25–1.51) for children 1–5 years of age given i.v. infusion of 1.6 mg/kg,32 0.43 (0.30–0.87) for children 1 month to 1 year of age given 1.25 mg/kg orally,14 and 0.70 (0.50–0.90) for neonates given 1.25 mg orally (Figure 3 a, Tables S2 and S3).33

Pantoprazole exposure was simulated in children between 2 and 14 years of age for both CYP2C19 EM and PM phenotypes using the pediatric PBPK model at a dose of 1 mg/kg pantoprazole. The predicted PK parameters for PM and EM populations were within twofold of the observed values (Table S1) with slightly overestimation of Cmax in the PM population.25

Esomeprazole

Figure 3 b shows the R values for CL when comparing the PBPK simulated results with reported values in the literature.42 We did not overlay the predicted and observed plasma concentration‐time profile data for esomeprazole in pediatrics because the above data are not available in the public domain. The esomeprazole pediatric PBPK model with enzyme auto‐inhibition reasonably described the clearance after multiple i.v. doses of esomeprazole in patients aged 6–17 years. However, the esomeprazole model underpredicted the CL (overprediction of AUC) for neonates and infants (R values for CL ratio were 0.15 and 0.33, respectively). The ratio for AUC was shown in Table S4. An improved prediction of AUC ratios, 9.25 (95% CI: 8.16–10.31) with CYP2C19 auto‐inhibition, vs. 1.19 (95% CI: 0.592–8.004) without CYP2C19 auto‐inhibition, was found for neonates (Figure 3 b, and Table S4). Consistently, the CL ratio in neonates was 0.15 (95% CI: 0.070–0.29) with CYP2C19 auto‐inhibition and was 1.00 (95% CI: 0.575–3.215) when CYP2C19 auto‐inhibition was not included. An improved prediction for age groups of 1–12 months and 1–5 years old was also observed when CYP2C19 auto‐inhibition was not included in the model.

Discussion

We followed a “learn, confirm, and refine” PBPK modeling strategy by evaluating the PBPK prediction performance with a pediatric pantoprazole model followed by confirming and refining with a pediatric esomeprazole model.16 These two probe substrates were chosen because esomeprazole is both a substrate and inhibitor of CYP2C19, whereas pantoprazole is only a substrate of CYP2C19. These models were based on a comparison of pantoprazole PK after single and multiple doses,33 and a previous model by Wu et al.28 in which a CYP2C19 auto‐inhibition was considered in an esomeprazole PBPK adult model. Our study found that clearance predictions were within a twofold range for pantoprazole in both the adult and pediatric populations using an established PBPK platform. However, the esomeprazole PBPK model with CYP2C19 auto‐inhibition significantly underpredicted CL in the younger age group of pediatric patients, especially in neonates. The differences in PBPK model predictive performance for the two drugs suggest the difficulty in the extrapolation of PK models between drugs sharing an elimination pathway.43 The PBPK model with CYP2C19 auto‐inhibition for esomeprazole, which has been verified in adults and older children, could not be directly extrapolated to neonates/infants. A good model prediction requires a thorough understanding of the complex interplay between CYP maturation and inhibition in infants.

In auto‐inhibition (mechanism‐based enzyme inhibition), enzymes are irreversibly removed from the active enzyme pool, thus showing a prolonged inhibition upon removal of the inhibitor. The only way to recover the enzyme activity is through the de novo synthesis. CYP2C19 is minimally expressed in the fetus and neonates, and the level of CYP2C enzymes, including CYP2C19, is almost undetectable in the first 24 hours after birth.44 The amount of CYP2C19 increases quickly after birth and reaches one third of the adult level after the first month of life.44 It is possible that the quick de novo synthesis of CYP2C19 in neonates provides more active enzyme and compensates for the enzyme auto‐inhibition by esomeprazole.

Our hypothesis is supported by significantly improved prediction of esomeprazole clearance in infants by not including the CYP2C19 auto‐inhibition in the esomeprazole PBPK model (“w/o inh” in Figure 3 b and Table S4). It should be noted that there are other mechanisms that may contribute to the underprediction of esomeprazole clearance in infants, including different protein binding capacity between neonates/infants and adults,43, 45 and a compensatory pathway mediated by CYP3A4 or other enzymes.46

There are a few limitations for our study. The models developed in the study did not include preterm infants due to limited observed data and a poor understanding of human physiology and enzyme/transporter ontogeny in this subpopulation. In addition, in vitro data suggest that pantoprazole might be a substrate of P‐glycoprotein (P‐gp).19 The effect of P‐gp on absorption is not considered in our current model due to the limited data regarding transporter maturation, and the contribution of P‐gp to absorption may not be significant because pantoprazole is a high permeability compound. Furthermore, pantoprazole clearance after i.v. and oral administration is similar, suggesting that intestinal efflux transporter or enzyme metabolism is minimal compared to hepatic clearance. Experimental data on the interplay between CYP2C19 de novo synthesis and enzyme inhibition would corroborate our conclusion. However, due to the difficulty in conducting experiments in younger children, and especially neonates, the experimental data in the literatures are limited. Future in vitro or in vivo studies are warranted to further understand the complex enzyme maturation mechanisms.

In summary, this study demonstrated different predictive performance of PBPK models in the neonates and infants for pantoprazole and esomeprazole, two drugs that share the same metabolic CYP2C19 pathway. These observations suggest that there is a complex interplay among CYP2C19 maturation, inhibition, and possible compensatory pathways in neonates and infants. Models involving substrates for CYP2C19 cannot be extrapolated to other CYP2C19 substrates in neonates and infants without verification. This age‐dependent interplay warrants further experimental investigation and modeling verification through the study of other CYP2C19 substrates and/or inhibitors.

Funding.

D.G. receives support for research from the National Institute of Child Health and Human Development (K23HD083465). J.M. was supported by NIH grant T32GM008562 during the drafting of the manuscript.

Conflict of Interest

The authors declared no competing interests for this work.

Author Contributions

P.D., F.W., J.N.M., J.F., V.C., D.G., G.J.B., L.Z., and J.W. wrote the manuscript. P.D., J.W., and F.W. designed the research. P.D., F.W., and J.N.J. performed the research. P.D., F.W., and J.N.J. analyzed the data.

Disclaimer.

The opinions in this article reflect the views of the authors and should not be interpreted as the position of the US Food and Drug Administration.

Table S1. Predicted or observed pantoprazole PK parameters in adult or pediatric population with different CYP2C19 genotype.

Table S2. Observed PK parameters of pantoprazole in adult and pediatric population (Mean ± SD).

Table S3. Predicted performance of the pantoprazole PBPK model in different pediatric age groups.

Table S4. Prediction of PK across pediatric age groups using the esomeprazole PBPK model with or without auto‐inhibition6.

Click here for additional data file.