Ustekinumab Dosing Individualization in Crohn’s Disease Guided by a Population Pharmacokinetic–Pharmacodynamic Model

Ustekinumab is a monoclonal antibody used in Crohn’s disease (CD). Dose optimization in case of non-response and the role of pharmacokinetic–pharmacodynamic (PK-PD) monitoring remain unresolved dilemmas in clinical practice. We aimed to develop a population PK-PD model for ustekinumab in CD and simulate efficacy of alternative dosing regimens. We included 57 patients and recorded their characteristics during 32 weeks after starting with ustekinumab therapy. Serum ustekinumab concentration was prospectively measured and fecal calprotectin (FC) concentration was used to monitor the disease activity. Ustekinumab PK-PD was described by a two-compartment target-mediated drug disposition model linked to an indirect response model. Lower fat-free mass, higher serum albumin, previous non-exposure to biologics, FCGR3A-158 V/V variant and lower C-reactive protein were associated with higher ustekinumab exposure. Model-based simulation suggested that 41.9% of patients receiving standard dosing achieve biochemical remission at week 32. In patients not achieving remission with standard dosing at week 16, transition to 4-weekly subcutaneous maintenance dosing with or without intravenous reinduction resulted in comparably higher remission rates at week 32 (51.1% vs. 49.2%, respectively). Our findings could be used to guide stratified ustekinumab treatment in CD, particularly in patients with unfavorable characteristics, who might benefit from early transition to 4-weekly maintenance dosing.


Introduction
Crohn's disease (CD) is a debilitating, relapsing-remitting, incurable inflammatory disease of the digestive tract [1]. Treatment goals have shifted from symptomatic improvement to a combination of endoscopic and clinical remission [2]. Biomarker remission was identified as an adjunct treatment goal with fecal calprotectin (FC) being sensitive for endoscopic disease activity [3,4].
Ustekinumab is an IgG1κ monoclonal antibody that binds with high affinity to the p40 subunit shared by interleukin 12 (IL12) and interleukin 23 (IL23) [5]. Standard dosing with weight-based intravenous induction followed by 90 mg subcutaneous injections every 8 weeks is not effective for all patients, with real-world studies reporting that only up to one half of the patients achieve clinical remission, and up to 25% of patients achieve endoscopic remission [6][7][8][9][10][11]. The exposure-response relationship between peak ustekinumab serum concentration after the induction dose and endoscopic remission at 24 weeks suggests that a subset of patients not responding to conventional dosing may benefit from dose escalation [11]. Recent real-world studies assessing off-label intensified maintenance therapy with 90 mg administered every four weeks in patients with insufficient response have shown favorable rates of clinical, biochemical and endoscopic remission [12][13][14][15].
It remains unclear which patients are best suited for dose optimization, what is the most suitable dosing strategy (intravenous re-induction or dosing interval shortening for subcutaneous therapy) and how pharmacokinetic-pharmacodynamic (PK-PD) monitoring could improve the effectiveness of ustekinumab. Moreover, clinical PK studies of ustekinumab in CD are scarce and data on how demographic, pathophysiological, and genetic factors impact the drug's PK and PD, are limited.
Our objective was to develop a population PK-PD model linking ustekinumab dosing regimen to the time course of FC concentrations. Through model-based simulations, which integrate both PK and PD data, we aimed to explore different dosing strategies, which are increasingly used in daily clinical practice.

Patients and Study Design
This was a prospective observational study in patients with CD starting treatment with ustekinumab at a single tertiary referral center. The study design and clinical results have partially been reported elsewhere [11]. Briefly, consecutive patients aged ≥18 years with CD who started treatment with ustekinumab between October 2017 and June 2019 were screened for eligibility and included in the current study, regardless of baseline endoscopic disease activity. They received a weight-based intravenous induction dose (≤55 kg: 260 mg; 55-85 kg: 390 mg; >85 kg: 520 mg) infused over one hour at baseline, followed by a fixed subcutaneous maintenance dose of 90 mg every 8 weeks. The duration of follow-up was 32 weeks.

Pharmacokinetic and Pharmacodynamic Data
Serum samples for the PK analysis were prospectively collected at baseline, 1 h after the end of the intravenous infusion (peak), and subsequently at weeks 2, 4, 8, 9, 10, 12, 16, 20, 24, and 32. Fecal samples of the first morning bowel movement were collected for the PD analysis at baseline, and subsequently at weeks 8, 16, 24, and 32. Ileocolonoscopies were performed within 2 months before starting ustekinumab to determine the presence of mucosal ulcerations to define endoscopically active luminal disease at baseline. Blood samples for genotyping analyses were collected at baseline.
Genomic DNA was isolated from whole blood collected in EDTA tubes by using a FlexiGene DNA kit (Qiagen, Hilden, Germany). The concentration and purity of DNA were measured on NanoDrop™ One/OneC Microvolume UV-Vis Spectrophotometer (Thermo Fischer Scientific, Waltham, MA, USA). Genotyping of single nucleotide polymorphisms (SNPs) was performed by TaqMan ® Pre-Designed SNP Genotyping Assays (Thermo Fischer Scientific, Waltham, MA, USA) on a LightCycler 480II real-time polymerase chain reaction instrument (Roche, Basel, Switzerland) according to manufacturer's recommendations. SNPs were selected based on previously identified associations with ustekinumab treatment outcomes in psoriasis or a possible role in the mediation of IgG clearance [16][17][18]. The analyzed SNPs and the specific assays used were: IL12B (rs3212227, assay ID: C_2084293_10; rs3213094, assay ID: C_29927086_10; and rs6887695, assay ID: C_1994992_10), FCGR2A (rs1801274, assay ID: C_9077561_20) and FCGR3A (rs396991, assay ID: C_25815666_10). To validate our results, 20% of the samples were re-genotyped for each SNP, and the results were found to be reproducible with no discrepancies noted.
Demographic characteristics (age, weight, height, sex), clinical data (co-morbidities, disease history and disease location, concomitant and previous treatment) and biochemical markers (C-reactive protein (CRP), serum albumin) were recorded. Fat-free mass (FFM) was determined using the Janmahasatian model [19].

Statistical Analysis and Pharmacokinetic-Pharmacodynamic Modeling
Descriptive statistics were used to present the data as non-normally distributed with medians and inter-quartile ranges (IQR). The nonlinear mixed-effects methodology and NONMEM ® software (version 7.3, Icon Development Solutions, Ellicott City, MD, USA) was used to analyze the PK and PD data and develop a population PK-PD model, linking ustekinumab concentrations, patients' characteristics and FC at all time points. A population PK model was developed using the ustekinumab concentrations converted to nanomolar concentration, assuming its molecular weight of 149 kDa. A linear oneand two-compartment models, and approximations of the target mediated drug disposition (TMDD) model were tested [20][21][22][23]. First-order absorption of ustekinumab was assumed after the subcutaneous dosing. Considering the concentration of ustekinumab target (p40 subunit of IL12 and IL23) was not measured, it was modeled as a latent variable in TMDD models. The PK studies in monkeys report a bi-exponential decline of the cytokines containing the p40 subunit (IL12 and IL23) [24,25]. Therefore, the PK data were extracted from these studies and PK parameters for the IL12 and IL23 were estimated and allometrically scaled to human. These estimates were used as initial estimates of the PK parameters of the target disposition in the TMDD models. One-and two-compartment models were tested to describe the distribution of the target.
Additive, proportional and combination (additive + proportional) error models were tested for residual variability. The logit transformation was used to describe the interindividual variability (IIV) of fraction of absorbed ustekinumab and the exponential model was used to describe the IIV of other parameters. The model with the lowest Akaike information criterion (AIC) was used as a base PK model (Table S1), for the subsequent covariate model building.
To explain the estimated IIV in the PK parameters, the candidate parameter-covariate relationships were selected based on scientific plausibility, previously reported relationships, and trends in correlation plots between individual PK parameters and covariates.
Stepwise covariate procedure (p < 0.05 in the forward inclusion, p < 0.01 in the backward elimination) was used to test the significance of parameter-covariate relationships. To estimate the effect of the body size, FFM was chosen over total body weight, considering the distribution of monoclonal antibodies is predominantly limited to extracellular fluids. The evaluated continuous covariates were: disease duration, FFM, baseline serum CRP concentration and serum albumin concentration. The evaluated categorical covariates were: previous biological therapy (bio-naïve), smoking, and SNPs in IL12B (rs3212227, rs3213094, rs6887695), FCGR2A (rs1801274) and FCGR3A (rs396991). Linear and power model were tested for continuous covariates and dominant and recessive grouping combinations were tested for SNPs.
The individual PK parameters obtained from the final PK model were used for the subsequent PD analysis [26], to describe the relationship between ustekinumab PK, target disposition and FC concentration. The bio-phase distribution model, indirect response Pharmaceutics 2021, 13, 1587 4 of 16 model and signal transduction (transit compartment indirect response) model were tested to describe the delay between ustekinumab PK and PD [27]. Exponential model, as well as Box-Cox transformation were tested to describe the IIV of baseline FC concentration (FC 0 ).
Laplacian estimation with interaction was used for parameter estimation. The M3 method was used for data below the lower or above the upper LOQ [28]. AIC was used for model comparison. The model was internally validated with a visual predictive check (VPC, n = 2000) and parameter uncertainty was assessed with bootstrap with replacement method (n = 2000).

Model Based Simulations
The final PK-PD model and NONMEM ® were used for simulation of various treatment regimens and scenarios of PK and PD monitoring. The proportion of patients achieving biochemical remission (FC < 100 mg/kg) at week 8, week 16, week 24 and week 32 after the first dose was estimated for every treatment regimen, to compare the efficacy. This cut-off was chosen based on test characteristics identified by meta-analyses [29], with emphasis on studies using the same assay as our center, considering the large inter-assay variability [30].
First, a virtual population of 10,000 CD patients was created, which resembled the observed cohort regarding the distributions and correlations of patients' characteristics. Subsequently, several clinically relevant scenarios were simulated, as follows: (a). All patients received standard ustekinumab treatment, with weight-based induction dose at baseline, followed by fixed 90 mg maintenance doses every eight weeks (standard treatment); (b). All patients received weight-based induction dose at baseline followed by fixed 90 mg maintenance doses every four weeks; (c). All patients received weight-based induction doses every eight weeks; (d). Patients receiving standard treatment who were not in remission at week 16, switched to maintenance doses every four weeks from week 20; (e). Patients receiving standard treatment who were not in remission at week 16 received a weight-based reinduction dose at week 16, and continued with maintenance doses every eight weeks; (f). Patients receiving standard treatment who were not in remission at week 16 received a weight-based reinduction dose at week 16 and switched to maintenance doses every four weeks from week 20.

Pharmacokinetic and Pharmacodynamic Data
A total of 574 serum samples was available for ustekinumab measurement. Five samples had an ustekinumab concentration below the lower LOQ, and none of the analyzed samples had measurable antibodies to ustekinumab. SNPs were analyzed for all patients, and all genotype frequencies were in Hardy-Weinberg equilibrium ( Table 1). A total of 224 samples was available for FC measurement. There were 15 samples with FC concentration below the lower LOQ and 11 samples above the upper LOQ.

Pharmacokinetic-Pharmacodynamic Model
A quasi-equilibrium approximation of the TMDD model [23] best described the time profile of the ustekinumab concentration (Table S1). The model ( Figure 1) assumes rapid binding of ustekinumab to the target antigen and was extended with distribution of the unbound ustekinumab and unbound target into a peripheral compartment since it resulted in an improved fit (Table S1). The binding was assumed in the central compartment. Estimated parameters comprised linear ustekinumab PK, linear target PK and parameters related to the binding of ustekinumab to the target ( Table 2). The terminal half-life of ustekinumab after the induction dose in a typical patient was estimated at 17 days and the volume of distribution was low, with volumes of central and peripheral compartment of ustekinumab estimated at 3.57 L and 3.30 L, respectively. There was a high IIV of the rate constant of target synthesis (K syn ), with the coefficient of variation estimated at 99.5%.
). e. For variability terms % shrinkage is presented in brackets.   Covariates associated with PK parameters were as follows: ustekinumab clearance was higher with increasing FFM, decreasing serum albumin, and in patients with previous exposure to biologics (Table 2). A power model described the increase of volumes of the central and peripheral compartment of ustekinumab with increasing FFM ( Table 2). The fraction of absorbed ustekinumab was higher in valine homozygous (V/V) patients of rs396991 FCGR3A polymorphism, compared to phenylalanine homozygotes (F/F) and heterozygotes (V/F) combined (Table 2, Figure S1). Additionally, K syn increased with increasing baseline CRP ( Table 2).
The concentration of the unbound target was linked to the FC concentration via the indirect response model, which best described the stimulating effect of the unbound target on the FC production (Table 3). Thus, the model explains the ustekinumab effect on the decrease in FC concentration via the intermediate target layer (Figure 1). The initial concentration of the target and FC were modelled as K syn /K deg and FC 0 , respectively. In the absence of ustekinumab, steady concentration of the target and FC is assumed. In contrast, when ustekinumab is administered and bound to the target, the unbound target concentration decreases, which inhibits the FC production. Exponential model best described the IIV on FC 0 and patients who started ustekinumab due to endoscopically active disease, had higher baseline FC concentration, which was accounted for in the PD model as a binary covariate on FC 0 (Table 3). There was a good agreement between model predictions and the observed data ( Figure 2 and Figure S2).

Simulations
The simulation of patients receiving standard ustekinumab regimen demonstrates that patients with weight ≤55 kg have lower ustekinumab concentration in the induction phase, and higher ustekinumab concentration in the maintenance phase, compared to patients with weight >55 kg. In addition, patients in the lowest weight group had higher unbound p40 and FC concentrations during the induction phase ( Figure S3).
After the first ustekinumab administration, the unbound target rapidly decreases. A subsequent increase in the target concentration is observed, due to the constant target synthesis. A faster increase in the target concentration in patients with weight >55 kg is observed, compared to patients with weight ≤55 kg, because of the fixed 90 mg subcutaneous dose in the maintenance phase and thus higher exposure to ustekinumab of patients with weight ≤55 kg ( Figure S3). In addition, initial target concentration influences the peak ustekinumab concentration. A patient with the 10-fold higher K syn compared to the   The dynamics of the FC are similar across the weight groups, with faster increase of the FC concentration in higher weight groups ( Figure S3). Patients with the endoscopically active disease at baseline had higher FC concentration, and only these patients were considered for the rest of simulations.
The proportions of patients in remission at week 32 receiving standard treatment with maintenance dosing every 8 weeks, treatment with initial induction and subsequent maintenance dosing every four weeks from week 4 or treatment with induction dosing every eight weeks are 41.9%, 52.2% and 56.0%, respectively (Figure 3). In patients not achieving remission with standard treatment at week 16, transition to subcutaneous maintenance dosing every 4 weeks with or without intravenous reinduction resulted in similar biochemical remission at week 32 (51.1% vs. 49.2%, respectively), while reinduction and continuation with maintenance dosing every 8 weeks was comparable to standard dosing (Figure 4). maintenance dosing every four weeks from week 4 or treatment with induction dosing every eight weeks are 41.9%, 52.2% and 56.0%, respectively (Figure 3). In patients not achieving remission with standard treatment at week 16, transition to subcutaneous maintenance dosing every 4 weeks with or without intravenous reinduction resulted in similar biochemical remission at week 32 (51.1% vs. 49.2%, respectively), while reinduction and continuation with maintenance dosing every 8 weeks was comparable to standard dosing (Figure 4).

Discussion
Clinical PK studies of ustekinumab are relatively scarce, and there is currently no consensus on the optimal management of inadequate response to treatment. To our knowledge, this is the first study to present a semi-mechanistic PK-PD model linking ustekinumab PK, CD patients' clinical characteristics, binding to p40 subunit and FC

Discussion
Clinical PK studies of ustekinumab are relatively scarce, and there is currently no consensus on the optimal management of inadequate response to treatment. To our knowledge, this is the first study to present a semi-mechanistic PK-PD model linking ustekinumab PK, CD patients' clinical characteristics, binding to p40 subunit and FC dynamics. Additionally, we are the first to evaluate the influence of genetic factors on ustekinumab PK-PD in CD. Model-based simulations suggest the possibility of achieving better outcomes with intensified dosing and provide further insight into managing patients with insufficient response at week 16.
Building on previous work identifying an exposure-response relationship for ustekinumab for biochemical and endoscopic outcomes [8,10,11,31], our PK-PD model suggests that rates of biochemical remission could be improved by using intensified dosing regimens to increase ustekinumab exposure (Figure 4). A subgroup of patients at risk of insufficient exposure to ustekinumab (albumin < 43 g/L, previously exposed to biological treatment, V/F or F/F genotype of rs396991, CRP > 5 mg/L) were predicted to benefit the most from dosing escalation. These patients could also be considered for potential proactive dosing escalation in future studies ( Figures S5 and S6). The proportion of patients achieving biochemical remission decreased with an increasing number of unfavorable covariates, and this trend was attenuated by dosing escalation ( Figure S6).
Different approaches have been used in patients not responding or losing response to standard dosing of ustekinumab: dosing interval shortening [12][13][14][15], intravenous reinduction [32] or a combination of both [32,33]. These modifications were made at week 16 of treatment or later. Clinical remission, defined as a Harvey-Bradshaw Index ≤ 4, was reported in 28-40% of patients after interval shortening to every 4 weeks [12][13][14], CRP normalization occurred in 22% of patients [12]. The interpretation of endoscopic remission rates around 35% (4/11, 14/39) is confounded by low numbers and potential selection bias where patients with biochemically active disease did not undergo endoscopic assessment [12,14]. None of the studies conducted so far were designed and powered to compare the efficacy of different strategies. Simulations based on our model showed that switching to maintenance dosing every 4 weeks from week 16 leads to similar rates of biochemical remission at week 32 regardless of an additional intravenous reinduction dose. A single additional reinduction dose followed by maintenance every 8 weeks is expected to provide only transient benefit ( Figure 4). This is broadly concordant with the findings of a retrospective real-world series, where 15/18 patients responded clinically to an additional intravenous infusion and 10/15 required escalated subcutaneous dosing in ongoing maintenance [32]. Even patients already on maximized subcutaneous dosing may benefit from an additional intravenous infusion, although none of these studies included PK assessment [33].
Our findings should be compared with the recently reported results of the STARDUST trial (NCT03107793) evaluating the efficacy of a treat-to-target (T2T) approach to guide treatment escalation decisions compared to standard of care (SOC) [34]. Clinical responders at week 16 underwent randomization, dosing in the T2T arm was first determined based on the change in endoscopic score from baseline and subsequently on combined clinical and biomarker (CRP, FC) monitoring. At week 48, endoscopic response (T2T: 37.7%; SOC: 29.9%) and remission (T2T: 11.4%; SOC: 14.5%) were not significantly different between the two groups. Notably, only 17% of patients in the T2T arm were on four-weekly dosing at week 48, although this percentage is also affected by the 20% discontinuation rate in this arm as patients not achieving treatment targets within four weeks of starting four-weekly dosing discontinued the study. It also remains to be seen what proportion of discontinuations was driven by clinical disease activity despite biochemical remission, as both were required to meet the target. Despite these results, dose escalation may still hold added benefit, particularly if applied early after induction to a selected subpopulation of patients unlikely to respond to standard dosing due to low drug exposure.
The estimated terminal half-life in our study is in accordance with PK studies in psoriasis, psoriatic arthritis, and IBD [35][36][37][38], reporting half-life between 19 and 22 days. The distribution of ustekinumab into tissues was limited, which is consistent with its high molecular weight and hydrophilic properties. Ustekinumab linear clearance was in range of reported clearance values in other studies in IBD [37][38][39], and the TMDD PK model is consistent with the recent findings of the TNF-mediated PK of infliximab in ulcerative colitis [40]. Notably, the elimination of ustekinumab in our model was captured in two catabolic pathways, distinct for IgG antibodies [41]: a linear, non-specific clearance, and a nonlinear, specific clearance, mediated by the binding of the antibody to its target. Therefore, a higher concentration of the target reflecting increased inflammatory burden may lead to accelerated ustekinumab clearance and lower ustekinumab concentration. In CD, the concentration of the p40 subunit is related to the disease activity, although the reported serum concentrations of IL12 and IL23 are highly variable [42,43].
The PK model additionally includes a bi-exponential decline of the concentration of the target, which reflects the possible distribution of the target to other tissues and is in line with the in vivo studies of IL12 and IL23 [24,25]. Thus, the initial binding of ustekinumab to the target in the serum may explain the influence of the target on the peak ustekinumab concentration and volume of distribution ( Figure S4). The subsequent redistribution of the target from other tissues into the serum, as well as the constant target synthesis, may explain the target-mediated ustekinumab clearance.
We have identified similar factors influencing ustekinumab PK compared to other studies [35][36][37][38]. Ustekinumab exposure was found to increase with decreasing FFM. However, dosing regimens seemed similarly efficient in patients stratified by the weight according to the induction dose ( Figure S6). This indicates that the weight is not related with ustekinumab efficacy and is in line with the recent post hoc analysis of IM-UNITI study, assessing the impact of body mass index on clinical efficacy [44]. The exposure increased with increasing serum albumin, as well. A more severe disease may result in a leakage of the albumin and drug into the bowel. The correlation may additionally be explained by the role of FcRn in homeostatic regulation of both, IgG and albumin [45]. Moreover, ustekinumab exposure was higher in bio-naïve patients. Presumably, these patients are less likely to have a more severe type of disease compared to those who have previously failed therapy with other biologics. Therefore, bio-naïve patients may express lower target concentrations, which leads to lower ustekinumab clearance and higher exposure.
We also observed higher fraction of absorbed ustekinumab, and higher exposure related with V/V variant of the rs396991 polymorphism, otherwise known as FCGR3A-158 V/F ( Figure S1, Table 2). This is a functionally significant polymorphism in FCGR3A gene, which encodes the FcγRIIIa receptor. Our finding is highly interesting, since it mirrors the studies reporting association between the V/V genotype and better outcomes of CD treatment with infliximab [46,47]. The V variant is presumably related with higher affinity between FcγRIIIa and the IgG antibodies, compared to the F variant, which may result in altered binding at the absorption site and lead to a variable response to treatment in IBD [48]. However, further studies are necessary to clarify the clinical relevance and the mechanistic involvement of FCGR3A-158 polymorphism in ustekinumab PK-PD.
The indirect response PD model for FC has been previously used in a PK-PD studies of infliximab and ustekinumab in CD [38,49]. Notably, the rates of FC degradation are similar, which may additionally account for the different FC measurement assays. In our study, all patients regardless of the endoscopic disease activity were included, and the presence of mucosal ulcerations at baseline was associated with higher baseline FC and CRP, indicating a higher inflammatory burden in patients with confirmed active luminal disease. Antibodies to ustekinumab were not detected in any of the tested samples, which is consistent with the low observed rates of immunogenicity in registration studies [50][51][52].
Specifically, our PK-PD model is in line with the recent study by Wang et al. [38], who reported a PK-PD model linking ustekinumab, FC and endoscopic outcomes, based on the sparse data from an exposure-response study by Verstockt et al. [10]. Notably, the estimated ustekinumab linear CL is similar (Wang et  Additionally, both studies used indirect response PD models with comparable K out estimate for FC (Wang et al.: 0.0416 day −1 ; vs. our study: 0.0581 day −1 ), although there was a large observed difference in the average FC concentration due to the different assay [53], and in our model the FC concentration was linked to the latent target instead of the ustekinumab. The residual unexplained PD variability was relatively high in both studies, at more than 55%. Nevertheless, both studies suggested the potential benefit of the shorter interval of the maintenance dosing.
Our findings should be interpreted in the context of the study's limitations. The model was developed in a small cohort of patients and the target concentration was not measured. Despite the plausibility of covariates explaining the PK of ustekinumab, which overlap with findings from other immune-mediated inflammatory diseases and models for other monoclonal antibodies in IBD, the model should be further validated in an independent cohort. Notwithstanding the correlation between endoscopic outcomes and FC concentrations, normalization of the latter is neither an independent treatment goal nor does its normalization guarantee endoscopic remission as observed in the STARDUST trial. Finally, outcomes of simulations based on our model should be regarded as hypothesisgenerating for future study design and not used to inform clinical practice here and now.
In conclusion, this article presents a semi-mechanistic population PK-PD model for ustekinumab in CD, and the findings could be used to guide future attempts to personalize treatment with ustekinumab in CD. Especially patients with low serum albumin, previous failure of biological therapy, FCGR3A-158 V/F or F/F variant and high CRP might benefit from early optimization of ustekinumab therapy to 4-weekly maintenance dosing.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/pharmaceutics13101587/s1, Figure S1: Fraction of absorbed ustekinumab after subcutaneous administration according to FCGR3A rs396991 polymorphism, Figure S2: Diagnostic plots of the final pharmacokinetic and pharmacodynamic model, Figure S3: A population simulation using the final pharmacokinetic-pharmacodynamic model and representative virtual patient population, Figure  S4: Serum unbound ustekinumab concentration time-profile in a patient with high, typical, or low target synthesis rate constant, Figure S5: Simulated ustekinumab and fecal calprotectin concentration time-profile of a typical patient according to ustekinumab exposure and target synthesis, Figure S6: Proportion of patients who achieved biochemical remission at week 32 stratified by each covariate, Table S1: Comparison of the tested base pharmacokinetic models, Supplementary information S1: NONMEM control stream for the population PK and PD model.