- Open Access
Saturable human neopterin response to interferon-α assessed by a pharmacokinetic-pharmacodynamic model
Journal of Translational Medicine volume 11, Article number: 240 (2013)
In this study, we developed a pharmacokinetic (PK)- pharmacodynamic (PD) model of a new sustained release formulation of interferon-α-2a (SR-IFN-α) using the blood concentration of IFN-α and neopterin in order to quantify the magnitude and saturation of neopterin production over time in healthy volunteers. The SR-IFN-α in this study is a solid microparticular formulation manufactured by spray drying of a feeding solution containing IFN-α, a biocompatible polymer (polyethylene glycol) and sodium hyaluronate.
The full PK and PD (neopterin concentration) datasets from 24 healthy subjects obtained after single doses of 9, 18, 27 and 36 MIU of subcutaneous SR-IFN-α were used to build the mixed-effect model using NONMEM (version 7.2) with the GFORTRAN compiler.
A one-compartment model with first-order elimination and a mixture of zero- and first-order absorption was chosen to describe the PK of SR-IFN-α. The time-concentration profile of neopterin, the PD marker, was described by a turnover model combined with a single transit compartment. The saturable pattern of the neopterin response blurring the dose–response relationship of SR-IFN-α was addressed by introducing the concept of the EC50 increasing over time.
The PK-PD model of SR-IFN-α developed in this study has presented a quantitative tool to assess the time-course of a saturable neopterin response in humans.
Interferons are produced by peripheral blood leukocytes, fibroblasts, and activated T and NK cells in response to viral infection or other inducers including double-stranded RNA, lipopolysaccharide, micro-organisms, or endotoxins [1–3]. Based on immunological and physicochemical differences, human interferons are divided into α-, β-, and γ-interferon families, with numerous subtypes within each interferon family . Because of its antiviral, antiproliferative and immunomodulating properties, recombinant interferon-α (IFN-α) has been used as a treatment for various diseases . However, in the case of chronic hepatitis C, monotherapy with IFN-α has been persistently effective in only a small percentage of patients . This low response rate is thought to be due to HCV genotype variation and/or the quite short half-life of IFN-α [7, 8]. Combination therapy with other antiviral agents, such as ribavirin, is therefore recommended . Frequent administration (3 times weekly) of IFN-α has been considered to be an additional cause for therapeutic failure of interferon due to the fact that frequent administration accelerates the formation of neutralizing antibodies and causes other adverse effects due to large variations in peak-to-trough plasma drug concentrations [10, 11]. Hence, a long-acting formulation has been developed for IFN-α which can extend its effects to weeks and months.
Unlike small molecule drugs, the poor stability of protein drugs has been a hurdle to the development of long acting formulations . Technologies for producing long-acting formulations of protein drugs can be roughly classified into chemical modification such as pegylation  or formulation changes allowing delayed release from depot sites. The new sustained release formulation of IFN-α-2a (SR-IFN-α, LG Life Sciences) used in our report is a solid microparticular formulation manufactured by the spray drying of a feeding solution containing IFN-α, a biocompatible polymer (polyethylene glycol) and sodium hyaluronate .
This report is based upon a first-in-human, single ascending dose trial of SR-IFN-α in healthy volunteers where the IFN-α concentration profile showed an extended release pattern through the doses studied. Blood neopterin concentrations were also measured as a biomarker demonstrating the activity of IFN-α in this first-in-human study. Neopterin is a soluble immune activation marker released from monocytes and macrophages by IFN-α . Looking into the relationship between the exposure to IFN-α and neopterin in the healthy subjects’ data, we found that the magnitude of neopterin concentration changes was not clearly correlated with the dose of IFN-α. Although similar phenomena have been observed in animal experiments [16, 17], this has never been reported in humans despite frequent clinical trials using neopterin as a biomarker. Because clear understanding of the relationship between exposure and response is one of the fundamental goals of early-phase exploratory clinical trials, we investigated the concentration-response relationship of IFN-α and neopterin in humans, which has never been elucidated. As results, we present a pharmacokinetic-pharmacodynamic (PK-PD) model that quantifies the peculiar time-course of neopterin responses to IFN-α in humans.
Inclusion and exclusion criteria
Healthy volunteers aged 18 to 45 years with BMI ranging from 19 to 29 kg/m2, with no clinically relevant conditions identified based on medical history, physical examination, laboratory tests or electrocardiography (ECG), were eligible for inclusion. Subjects with any history that indicated a possible alteration in IFN-α metabolism or with hypersensitivity to IFN-α were excluded. The final study enrollment was 32 subjects (Table 1).
A randomized, double-blind, active controlled, dose escalation phase I clinical study was conducted on 32 healthy subjects in the clinical pharmacology unit of the Kendle International BV, located in Utrecht, Netherlands. Subjects were randomly allocated into four groups (eight subjects per group). Within each group, six were given SR-IFN-α (test formulation) and the other two were given 3 MIU Roferon-A® (Roche, active comparator) via subcutaneous injection. The doses of SR-IFN-α allocated to groups 1, 2, 3 and 4 were 9 MIU, 18 MIU, 27 MIU, and 36 MIU, respectively.
The study was performed in compliance with the European Community rules of Good Clinical Practice (GCP), the International Conference on Harmonization (ICH) Tripartite Guidelines: Guideline for GCP, the current revision of the 'Declaration of Helsinki’ (Edinburgh, amendment October 2000). The Stichting Therapeutische Evaluatie Geneesmiddelen (STEG), an independent ethics committee, approved the protocol before execution of the trial, and all participants gave written informed consent.
For the population PK analysis, peripheral blood samples (5 mL each) were taken just prior to the injection, and 0.75, 1.5, 3, 6, 8, 10, 12, 18, 24, 30, 36, 48, 60, 72, 96, 120, 144, 168, and 192 hours after the injection. For the PD marker, the sampling times differed slightly: just prior to the injection and 3, 8, 12, 18, 24, 36, 48, 72, 96, 120, 144, 168, 192, and 264 hours after the injection. The samples were collected in light-protective tubes and stored at < -70°C.
Assay of plasma concentrations of IFN-α and neopterin
IFN-α concentrations in the serum samples were analyzed using a commercial Human IFN-α Multi-Subtype ELISA Kit (product # 41105) with a detection limit of 12.5 pg/mL manufactured by Pestka Biomedical Laboratories, Inc. (Piscataway, NJ, USA). Neopterin concentrations in the serum samples were analyzed using a commercially-available Neopterin ELISA method (REF 40-371-25012, GenWay Biotech, Inc., San Diego, CA, USA) with a detection limit of 0.7 nmol/L and the specificity of about 99.95%.
Population PK-PD model
Because the aim of this study was to develop a PK-PD model for SR-IFN-α, data from the active control group participants, who were given the immediate release IFN-α formulation (8 subjects), were not included in the analysis (individual plots for PK-PD models are shown in Additional file 1).
Mean plasma concentrations of IFN-α are shown in Figure 1, and non-compartmental analysis results of the PK of SR-IFN-α are summarized in Table 2. The population PK-PD analysis was performed using NONMEM (version 7.2, Icon Development Solutions, Ellicott City, MD, USA) with the GFORTRAN compiler.
To find the model that best described the absorption profile, which showed double peaks in many subjects, first- and zero-order absorption models and their combined form, with or without lag time, were tested. Based on first-order elimination, one- and two-compartment distribution models were tested for the three absorption (first-order, zero-order and combined) processes. The Michaelis-Menten absorption and elimination models were also tested considering the potential saturable absorption or elimination using the ADVAN subroutines (Table 3).
The coupling of the PK model to the PD model was done in a sequential manner. The population PD modeling was performed using the individual PK parameters estimated from the final PK model, which were added to the PD dataset.
Turnover models, with or without transit compartments, were compared to find the most appropriate model that explained the delayed effect of IFN-α on neopterin concentrations. The turnover model was initially selected over the effect compartment model based upon the well-known action of IFN-α stimulating the release of neopterin, and the transit compartments were tested because their usefulness was reported in a previous preclinical study .
The differential equations for the drug effect model (model structures are shown in Figure 2) were:
where Ka and Ke are the absorption and elimination rate constants for IFN-α, respectively. Kin is the production rate of neopterin, a zero-order constant, Ktr is the first-order transition rate leaving the transit compartment, Kout is the first-order rate constant for the elimination of neopterin, and E(C) is the effect as a function of the individual predicted drug concentration, C. A(4), which is the amount existing in the 4th compartment, indicates the concentration of neopterin.
The stimulatory function for the drug effect, E(C), was a sigmoid function:
where Emax is the maximum effect and EC50 is the IFN-α concentration that produces 50% of the maximum effect.
In our study, differences in the mean plasma neopterin concentrations between the groups receiving different doses of SR-IFN-α were not clearly discernible (Figure 3). To account for this phenomenon in our PD model, we incorporated the concept of time-dependent attenuation of the effect parameters , especially the increasing EC50 over time, as shown in the following equation:
where ECB is the baseline of EC50 and CA and CB are coefficients to describe the EC50 increase to a certain level (ECB × CA) over time. This curve is the cumulative distribution function (CDF) of an exponential distribution to EC50. It has advantages in explaining the concave-shaped curves in relation to time.
The first-order conditional estimation (FOCE) method with interaction was used throughout the model building. Models were selected based upon a decrease in the objective function value (OFV) of more than 3.84 (P-value 0.05 in an approximate χ 2 distribution) and improvement in the individual plots, as well as other scatterplots.
A log normal distribution was assumed for inter-individual variability (η), and PK or PD parameters of the jth subject (Pj) were described as:
where TVP represents the typical population value of PK-PD parameters, such as clearance (CL), volume of distribution (V), absorption rate constant (Ka), lag time (ALAG), the first-order elimination rate of serum neopterin (Kout), and maximum stimulation effect (Emax). The inter-individual variability eta (η) for each PK-PD parameter was assumed to follow a Gaussian distribution with a mean 0 and a variance ω2. Possible correlations between the inter-individual variability were also evaluated.
As for the residual error, the additive, proportional and combined forms were tested. An example of the combined error form is shown as follows:
where IPREDij is the individual predicted concentration, Yij is the measured concentration of the jth individual at the ith sampling time, and ϵij is residual error. Residual errors (ϵ) include intra-individual variability, assay error and model misspecification. They were also assumed to follow a Gaussian distribution with a mean 0 and variance σ2.
Age, height, weight and creatinine clearance were screened as potential covariates of the parameters using Generalized Additive Modeling (GAM) implemented by Xpose version 4.2.3.
In the forward selection of covariates, variables that decreased the OFV by more than 3.84 (P < 0.05) and improved the inter-individual variability (omega value decrease) were selected. Covariates that did not increase the OFV more than 3.84 (P < 0.05) in the backward elimination step were removed from the model.
The 95% confidence intervals (CIs) for mean population PK and PD parameters were determined using a re-sampling technique based on the bootstrap method. One thousand re-sampled datasets were collected and their parameters were estimated using our population models. The models were also evaluated by visual predictive checks (VPC) using 1,000 simulated datasets.
Population PK-PD model
The results of non-compartmental PK analysis showing the trend of PK linearity of SR-IFN-α are briefly summarized in Table 2. A one-compartment model with first-order elimination and a mixture of zero- and first-order absorption best described the PK of SR-IFN-α. The delayed pattern of the time-concentration profile of neopterin was well described by the turnover model with a single transit compartment. There was no significant covariate. The structure of the final PK-PD model is shown in Figure 2. The population PK-PD parameter estimates, with corresponding standard error (SE) values, are summarized in Table 4. Basic goodness-of-fit plots are presented in Figure 4. Predicted concentration-time profiles of IFN-α and observed data from representative individuals are shown in Figure 5.
The 95% confidence intervals (CIs) for population PK-PD parameter estimates, determined using the bootstrap re-sampling method, are shown in Table 4. VPC plots simulated concentrations of 1,000 virtual datasets (nsub=1000 in the $SIMULATION block, 24,000 virtual patients) from the final model. The results from the VPC showed that the PK-PD model gave acceptable predictive performance. Curves for the 12.5th, 50th and 87.5th percentiles of concentrations were overlaid on the observed concentrations (Figure 6).
In this study, we presented results of PK-PD modeling of the time-concentration profiles of IFN-α and neopterin after administration of SR-IFN-α in healthy subjects.
For PK characteristics, the combined absorption model was successful in describing the double peak phenomenon. A zero-order absorption model appropriately described the initial increase in concentration, as measured by relatively frequent sampling, before reaching the maximum concentration. Subsequent second peaks (observed around 100 h after injection) were, however, best described by a first-order absorption model. One possible explanation for this double-peak phenomenon is in the method of SR-IFN-α administration: SR-IFN-α should be mixed with medium chain triglycerides (MCT) right before subcutaneous injection. Micro-droplets with various sizes might be formed in this mixing step, and their absorption rates may differ by droplet diameters. Despite differences in formulation, the CL (12.2 L/h) and V (691 L) in our report were not much different from those reported in a previous study for IFN-α (7~8 L/h and 700~850 L, respectively) in healthy subjects .
As for the PD model, the addition of a transit compartment before the neopterin compartment gave better outcome when compared with a simple turnover model. The transit compartment for neopterin was used to model data from monkeys , and we found that it is also useful in a human PK-PD model in this study. As there are a few mechanistic models that tried to explain little PD differences between dose groups (saturation of responses) [16, 20, 21], we tested them for our neopterin data in the preliminary PD model development step; however, none of them were successfully converged by NONMEM. Neither the precursor turnover model  that explained tolerance with depletion of precursor molecules, nor the turnover feedback model  that explained tolerance with negative feedback via a moderator compartment provided acceptable parameter estimates, and the basic goodness of fit plots were even worse than those for the descriptive model of increasing EC50. The inhibitory feedback model in monkeys  also showed similar problems. Thus, we had to use the concept of time-dependent attenuation of EC50 that did not include mechanistic reasoning to describe our neopterin response.
There are conflicting reports on saturation of neopterin by IFN-α in patients. A report on patients with hairy cell leukemia showed that neopterin responses and clinical efficacies after low doses (0.5-0.8 MIU/day) were similar to those after a conventional dose (3 MIU/day) in a 6-month clinical trial . In another clinical trial for a controlled release formulation of IFN-α, increases in the mean AUC of neopterin were marginal among the doses tested (4.09 mM·h after 20 μg IFN-α and 6.61 mM·h after 320 μg IFN-α) . The forms of IFN-α used in those studies [22, 23] were non-pegylated forms, like the SR-IFN-α used in this report; however, in phase I clinical trials of pegylated IFN-α, neopterin responses were well correlated with the doses used [24, 25]. Such a discrepancy in neopterin responses between non-pegylated and pegylated IFN-α suggests that the polyethylene glycol tail attached to IFN-α changed its PD parameters related to neopterin production. To the best of our knowledge, although pegylated formulations have long been used, it has never been reported that their PD profile may be different from that of non-pegylated forms.
Because neopterin is a frequently-used marker of cell-mediated immunity, it can be used to monitor the degree of immune activation in various clinical conditions, including infections, autoimmune diseases, malignancies, and other conditions . Neopterin is also known to mediate the cytotoxic action of activated macrophages and dendritic cells via interactions with reactive oxygen species . Thus, saturation of the neopterin response suggests that the magnitude of cytotoxic action mediated by neopterin may be similar regardless of the doses of IFN-α.
We developed a human PK-PD model revealing the saturable neopterin response to IFN-α for the first time. Our model suggests that the magnitude of cytotoxic action mediated by neopterin may be similar regardless of the doses of non-pegylated IFN-α.
Pfeffer LM, Dinarello CA, Herberman RB, Williams BR, Borden EC, Bordens R, Walter MR, Nagabhushan TL, Trotta PP, Pestka S: Biological properties of recombinant α-interferons: 40th anniversary of the discovery of interferons. Cancer Res. 1998, 58: 2489-2499.
Stark GR, Kerr IM, Williams BR, Silverman RH, Schreiber RD: How cells respond to interferons. Annu Rev Biochem. 1998, 67: 227-264. 10.1146/annurev.biochem.67.1.227.
Bao Q, Zhao Y, Niess H, Conrad C, Schwarz B, Jauch K-W, Huss R, Nelson PJ, Bruns CJ: Mesenchymal stem cell-based tumor-targeted gene therapy in gastrointestinal cancer. Stem Cells Dev. 2012, 21: 2355-2363. 10.1089/scd.2012.0060.
Samarajiwa SA, Wilson W, Hertzog PJ: Type I interferons: genetics and structure. 2006, Characterization and Application: The Interferons
Jonasch E, Haluska FG: Interferon in oncological practice: review of interferon biology, clinical applications, and toxicities. Oncologist. 2001, 6: 34-55. 10.1634/theoncologist.6-1-34.
Hoofnagle JH, Lau D: Chronic viral hepatitis — benefits of current therapies. N Engl J Med. 1996, 334: 1470-1471. 10.1056/NEJM199605303342210.
Lutchman G, Hoofnagle JH: Viral kinetics in hepatitis C. Hepatology. 2003, 37: 1257-1259. 10.1053/jhep.2003.50238.
Liang TJ, Rehermann B, Seeff LB, Hoofnagle JH: Pathogenesis, natural history, treatment, and prevention of hepatitis C. Ann Intern Med. 2000, 132: 296-305. 10.7326/0003-4819-132-4-200002150-00008.
McHutchison JG, Gordon SC, Schiff ER, Shiffman ML, Lee WM, Rustgi VK, Goodman ZD, Ling M-H, Cort S, Albrecht JK: Interferon alfa-2b alone or in combination with ribavirin as initial treatment for chronic hepatitis C. N Engl J Med. 1998, 339: 1485-1492. 10.1056/NEJM199811193392101.
Milella M, Antoneill G, Santantonio T, Currenti M, Monno L, Mariano N, Angarano G, Dianzanl F, Pastors G: Neutralizing antibodies to recombinant alpha‒interferon and response to therapy in chronic hepatitis C virus infection. Liver. 1993, 13: 146-150.
Caliceti P: Pharmacokinetics of pegylated interferons: what is misleading?. Dig Liver Dis. 2004, 36: S334-S339.
Roberts M, Bentley M, Harris J: Chemistry for peptide and protein PEGylation. Adv Drug Deliv Rev. 2002, 54: 459-476. 10.1016/S0169-409X(02)00022-4.
Grace M, Youngster S, Gitlin G, Sydor W, Xie L, Westreich L, Jacobs S, Brassard D, Bausch J, Bordens R: Structural and biologic characterization of pegylated recombinant IFN-α2b. J Interferon Cytokine Res. 2001, 21: 1103-1115. 10.1089/107999001317205240.
Choi SY, Jeh HS: Sustained release composition of protein drug. 2006,WO Patent 2,006,088,336,
Murr C, Widner B, Wirleitner B, Fuchs D: Neopterin as a marker for immune system activation. Curr Drug Metab. 2002, 3: 175-187. 10.2174/1389200024605082.
Hu X, Olivier K, Polack E, Crossman M, Zokowski K, Gronke RS, Parker S, Li Z, Nestorov I, Baker DP: In vivo pharmacology and toxicology evaluation of polyethylene glycol-conjugated interferon β-1a. J Pharmacol Exp Ther. 2011, 338: 984-996. 10.1124/jpet.111.180661.
Pepinsky RB, LePage DJ, Gill A, Chakraborty A, Vaidyanathan S, Green M, Baker DP, Whalley E, Hochman PS, Martin P: Improved pharmacokinetic properties of a polyethylene glycol-modified form of interferon-β-1a with preserved in vitro bioactivity. J Pharmacol Exp Ther. 2001, 297: 1059-1066.
Colburn WA, Eldon MA: Simultaneous pharmacokinetic/pharmacodynamic modeling. 1994, New York: Wiley
García-García I, González-Delgado CA, Valenzuela-Silva CM, Díaz-Machado A, Cruz-Díaz M, Nodarse-Cuní H, Pérez-Pérez O, Bermúdez-Badell CH, Ferrero-Bibilonia J, Páez-Meireles R: Pharmacokinetic and pharmacodynamic comparison of two "pegylated" interferon alpha-2 formulations in healthy male volunteers: a randomized, crossover, double-blind study. BMC Pharmacol. 2010, 10: 15-
Licko V, Ekblad E: Dynamics of a metabolic system: what single-action agents reveal about acid secretion. Am J Physiol Gastrointest Liver Physiol. 1992, 262: G581-G592.
Zuideveld KP, Maas HJ, Treijtel N, Hulshof J, van der Graaf PH, Peletier LA, Danhof M: A set-point model with oscillatory behavior predicts the time course of 8-OH-DPAT-induced hypothermia. Am J Physiol-Regul, Integ Comp Physiol. 2001, 281: R2059-R2071.
Gastl G, Aulitzky W, Tilg H, Nachbaur K, Troppmair J, Flener R, Huber C: A biological approach to optimize interferon treatment in hairy cell leukemia. Immunobiology. 1986, 172: 262-268. 10.1016/S0171-2985(86)80107-3.
De Leede LG, Humphries JE, Bechet AC, Van Hoogdalem EJ, Verrijk R, Spencer DG: Novel controlled-release lemna-derived IFN-α 2b (locteron): pharmacokinetics, pharmacodynamics, and tolerability in a phase I clinical trial. J Interferon Cytokine Res. 2008, 28: 113-122. 10.1089/jir.2007.0073.
Motzer RJ, Rakhit A, Ginsberg M, Rittweger K, Vuky J, Yu R, Fettner S, Hooftman L: Phase I trial of 40-kd branched pegylated interferon alfa-2a for patients with advanced renal cell carcinoma. J Clin Oncol. 2001, 19: 1312-1319.
Glue P, Fang JW, Rouzier-Panis R, Raffanel C, Sabo R, Gupta SK, Salfi M, Jacobs S: Pegylated interferon-α2b: pharmacokinetics, pharmacodynamics, safety, and preliminary efficacy data. Clin Pharmacol Ther. 2000, 68: 556-567. 10.1067/mcp.2000.110973.
Berdowska A, Zwirska‒Korczala K: Neopterin measurement in clinical diagnosis. J Clin Pharm Ther. 2001, 26: 319-329. 10.1046/j.1365-2710.2001.00358.x.
Hoffmann G, Wirleitner B, Fuchs D: Potential role of immune system activation-associated production of neopterin derivatives in humans. Inflamm Res. 2003, 52: 313-321. 10.1007/s00011-003-1181-9.
JHJ is an employee of LG Life Sciences Co., Ltd.
SJ performed the PK-PD analysis and drafted the manuscript. JHJ participated in the planning and design of the clinical study and drafting the manuscript. SH, JL, TH, and JP participated in PK-PD analysis. DSY conceived the study, supervised the process of PK-PD analysis and edited the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.