the multiple - dose pharmacokinetics of fe  999049 have been described , accounting for endogenous follicle - stimulating hormone ( fsh ) levels.exposure to fe  999049 was influenced by body weight , which may be a potential factor for individualised dosing recommendations.endogenous fsh levels were influenced by progesterone and inhibin  b levels .
this time - varying contribution of endogenous fsh may be important to consider when characterising the pharmacokinetics of recombinant fsh products .
follicle - stimulating hormone ( fsh ) is a gonadotropin that is synthesised and secreted by the anterior pituitary gland .
the major function of fsh is to regulate the reproductive processes by stimulating the gonads . in females
, fsh stimulates follicular development in the ovaries and production of inhibin  b , progesterone , and oestrogens by the ovarian follicular granulosa cells . luteinising hormone ( lh ) , another gonadotropin from the anterior pituitary , stimulates the theca cells of the follicles to deliver androgens to the granulosa cells for conversion to oestrogens .
it is also responsible for ovulation of the dominant follicle that has reached a fully mature , pre - ovulatory stage .
the ovarian hormones promote further follicular development as well as exerting negative and positive feedback loops to the hypothalamus and pituitary affecting the gonadotropin production and secretion .
in addition , gonadotropin secretion is stimulated by gonadotropin - releasing hormone ( gnrh ) produced in the hypothalamus .
female infertility can be caused by numerous factors at any level , from the hypothalamus to the ovaries and uterus .
gonadotropin therapy with either menotropins or recombinant fsh ( rfsh ) preparations can be used for infertility treatment when the cause is not primary ovarian failure , such that the ovaries are still responsive and contain primordial follicles .
the purpose of controlled ovarian stimulation with daily administration of gonadotropins prior to assisted reproductive technologies such as in vitro fertilisation ( ivf ) or intracytoplasmic sperm injection ( icsi ) is to obtain an adequate number of oocytes per retrieval , with minimum risks for the patient .
an appropriate ovarian response leading to the availability of several embryos makes it possible to select the best one(s ) for transfer .
recently , a novel recombinant human fsh ( rhfsh ; fe 999049 , ferring pharmaceuticals a / s ) was expressed for the first time in a human cell line ( per.c6 ; crucell , leiden , the netherlands ) , while existing rfsh preparations in clinical use ( e.g. follitropin- and follitropin- ) are derived from chinese hamster ovary cell lines ( cho ) .
previously , in a population pharmacokinetic analysis of first - in - human data after single ascending doses , body weight was identified as a factor that negatively correlates with serum fe  999049 concentration . in the present work , the fe  999049 population pharmacokinetics
the influence of body weight on drug exposure was investigated in order to further explore its importance as a factor for individualised fe  999049 dosing . moreover
, endogenous fsh levels may vary over time across the target patient group of infertile women , which may be an important aspect to consider in rfsh drug development and when analysing hormone data .
the endogenous fsh contribution to the total fsh levels , and the influence of other reproductive hormones on the endogenous fsh levels , were therefore also evaluated in this analysis .
data were generated in a randomised , double - blind , active control , multiple - dose trial with the aim of investigating the safety , tolerability , immunogenicity , pharmacokinetics , and pharmacodynamics of fe  999049 in healthy women .
the trial was performed according to the helsinki declaration and good clinical practice , and was approved by regulatory authorities and local ethics committees .
the trial has been described in more detail in a recent publication comparing the pharmacokinetic and pharmacodynamic properties of fe  999049 and gonal - f ( follitropin- ; emd serono , inc . ,
briefly , 49 healthy women were administered daily subcutaneous doses of 225 iu rfsh for 7  days ; 24 of the 49 women were treated with fe  999049 and 25 women received gonal - f as an active comparator . prior to the trial ( day 28 and 14 ) ,
subjects were administered two doses of a gnrh agonist ( lupron depot , 1-month depot ) to downregulate endogenous fsh .
blood samples for fsh , inhibin b , estradiol , progesterone , and lh measurements were collected at 60 and 30  min prior to administration of fe  999049 , immediately before administration , and once daily for 15  days .
in addition , after administration of the last dose on day  6 and until day  8 , the fsh concentration was measured every 4  h. analysis of serum fsh concentrations was performed at ferring pharmaceuticals a / s using a validated immunoassay based on electrochemiluminescence ( msd sector imager 2400 ) with a lower limit of quantification ( lloq ) of 0.075  g / l .
three of the 672 fsh measurements ( 0.4  % ) were below the quantification limit and were excluded from further analysis .
the personal demographics and baseline characteristics for the included subjects are listed in table  1 .
subjects with missing hormone baseline values were given the median population baseline value . between days  6 and 8 ,
fsh was measured every 4  h , and the other hormones were measured only once daily , leaving missing hormone values in between . to fill out the extra time points , the last measured hormone values were carried forward.table  1summary of subject characteristicsmean(range)age ( years)31.1(21.538.7)height ( cm)163.5(149.0175.3)weight ( kg)71.5(46.186.6)bmi ( kg / m)26.6(20.828.9)fsh ( g / l)0.211(0.0890.376)e2 ( pg / ml)24.21(7.055.0)lh ( iu / l)0.719(0.11.4)prog ( g / l)0.652(0.081.49)inhb ( pg / ml)22.0(4.089.0)demographics and measured baseline hormone levels for the 24 subjects receiving fe 999049 
 bmi body mass index , fsh follicle - stimulating hormone , e2 estradiol , lh luteinising hormone , prog progesterone , inhb inhibin b summary of subject characteristics demographics and measured baseline hormone levels for the 24 subjects receiving fe 999049 
 bmi body mass index , fsh follicle - stimulating hormone , e2 estradiol , lh luteinising hormone , prog progesterone , inhb inhibin b the pharmacokinetic model was developed using nonlinear mixed - effects modelling , where both the population parameters , interindividual variability ( iiv ) , and residual errors are estimated . for parameters with iiv , the ith subject s individual parameter ( i ) is log - normally distributed:1\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \theta_{i } = \theta \exp \left ( { \eta_{i } } \right )
$ $ \end{document}i=expiwhere  is the typical population parameter and i is the individual random effect from an approximately normal distribution with mean zero and variance  for describing the iiv of the parameter . in the model , potential influential factors
can be tested for significance as a covariate to explain some of the iiv in a parameter .
thus , the set of individual parameters ( i ) is given as a function of the typical population parameters (  ) individual values of the covariates ( ci ) and random effects ( i):2\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \vartheta_{i } = g\left ( { \vartheta , { \mathbf{c}}_{i } , \varvec{\eta}_{i } } \right ) $ $ \end{document}i = g,ci,ithe residual errors (  ) are assumed normally distributed with mean zero and variance  , and are the unexplained deviations of model predictions from the observations:3\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ y_{ij } = f\left ( {   t_{ij } , \vartheta_{i } } \right ) + \varepsilon_{ij } $ $ \end{document}yij = ftij,i+ijwhere
yij is the observation at time tij and the subscript ij denotes the jth number for subject i. the individual model prediction f (  ) at time tij is calculated from the individual parameters ( i ) .
model development was guided by changes in the nonmem objective function value ( ofv ) , precision of parameter estimates , and graphical model goodness - of - fit assessments , including visual predictive checks ( vpc ) . in the vpc ,
the observed data are compared with model predictions based on 1000 simulated trial datasets using the final model .
it displays the observations and the 2.5th , 50th , and 97.5th percentiles of observations , and the 95  % confidence intervals ( ci ) for the corresponding model predictions are plotted .
the difference in ofv ( dofv ) between two nested models is approximately chi square - distributed , with degrees of freedom ( df ) equal to the difference in the number of parameters .
based on this , the statistical significance for inclusion / exclusion of a model parameter can be judged .
a significance level of 0.05 was used for discrimination among nested models and covariate testing , corresponding to a 3.84 change in ofv for 1  df .
measurable fsh levels before drug administration indicated that endogenous fsh was not fully suppressed in this trial . the total fsh amount at time t in the central compartment was therefore modelled as the sum of the endogenous fsh and the exogenously administered rhfsh ( fe  999049):4\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { \text{fsh}}(t ) = { \text{fsh}}_{\text{en } } \left ( t \right )    { \text { + rhfsh}}\left ( t \right ) $ $ \end{document}fsh(t)=fshent+ rhfsht the endogenous fsh was assumed to have the same elimination rate constant from the central compartment as fe  999049 ; thus , the total fsh is eliminated from the central compartment with rate constant k given by clearance ( cl ) divided by volume of distribution ( v ) . since data are obtained after subcutaneous dosing , the bioavailability ( f ) is not known .
the cl and v estimated here are therefore the apparent cl ( cl / f ) and apparent v ( v / f ) .
the change in endogenous fsh amount over time in the central compartment was described as a turnover model with a zero - order production rate kin:5\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \frac{{{\text{dfsh}}_{\text{en } } \left ( t \right)}}{{{\text{d}}t } } = k_{\text{in } } - k \cdot { \text{fsh}}_{\text{en } } \left ( t \right ) $ $ \end{document}dfshentdt = kin - kfshent at time zero , before rhfsh dosing , the measured concentration in the central compartment is purely endogenous fsh .
the endogenous fsh baseline concentrations ( fshbl ) for each subject can therefore be estimated as the initial concentration in the central compartment and is the initial condition for the differential equation for the central compartment fshblv / f , i.e. the individual endogenous fsh amount at baseline .
when assuming steady state and hence zero change in the total fsh concentration before dosing , individual values of kin can be calculated from the other estimated parameters and baseline values .
one- and two - compartment models were tested for describing the distribution of fe  999049 .
a lag - time and transit compartments were investigated for modelling a potential delay in the absorption process .
body weight , age , and hormone ( estradiol , inhibin b , progesterone ) values at baseline were tested as potential covariates , explaining some of the iiv in parameter estimates .
inhibin b was also tested as a covariate over time since it has a purely inhibitory effect on fsh .
the measured hormone levels were assumed to not affect the rhfsh concentration , but only tested as covariates on the endogenous fsh baseline parameter ( fshbl ) . besides a decrease in ofv , the significance of a covariate was also evaluated by looking at the reduction in the iiv , measured as a coefficient of variation ( cv ) for the parameter s random effect .
the models were implemented and parameters estimated in nonmem 7.2.0 ( icon development solutions , ellicott city , md , usa ) . data handling and
data were generated in a randomised , double - blind , active control , multiple - dose trial with the aim of investigating the safety , tolerability , immunogenicity , pharmacokinetics , and pharmacodynamics of fe  999049 in healthy women .
the trial was performed according to the helsinki declaration and good clinical practice , and was approved by regulatory authorities and local ethics committees .
the trial has been described in more detail in a recent publication comparing the pharmacokinetic and pharmacodynamic properties of fe  999049 and gonal - f ( follitropin- ; emd serono , inc . ,
briefly , 49 healthy women were administered daily subcutaneous doses of 225 iu rfsh for 7  days ; 24 of the 49 women were treated with fe  999049 and 25 women received gonal - f as an active comparator . prior to the trial ( day 28 and 14 ) ,
subjects were administered two doses of a gnrh agonist ( lupron depot , 1-month depot ) to downregulate endogenous fsh .
blood samples for fsh , inhibin b , estradiol , progesterone , and lh measurements were collected at 60 and 30  min prior to administration of fe  999049 , immediately before administration , and once daily for 15  days .
in addition , after administration of the last dose on day  6 and until day  8 , the fsh concentration was measured every 4  h. analysis of serum fsh concentrations was performed at ferring pharmaceuticals a / s using a validated immunoassay based on electrochemiluminescence ( msd sector imager 2400 ) with a lower limit of quantification ( lloq ) of 0.075  g / l .
three of the 672 fsh measurements ( 0.4  % ) were below the quantification limit and were excluded from further analysis .
the personal demographics and baseline characteristics for the included subjects are listed in table  1 .
subjects with missing hormone baseline values were given the median population baseline value . between days  6 and 8 ,
fsh was measured every 4  h , and the other hormones were measured only once daily , leaving missing hormone values in between . to fill out the extra time points , the last measured hormone values were carried forward.table  1summary of subject characteristicsmean(range)age ( years)31.1(21.538.7)height ( cm)163.5(149.0175.3)weight ( kg)71.5(46.186.6)bmi ( kg / m)26.6(20.828.9)fsh ( g / l)0.211(0.0890.376)e2 ( pg / ml)24.21(7.055.0)lh ( iu / l)0.719(0.11.4)prog ( g / l)0.652(0.081.49)inhb ( pg / ml)22.0(4.089.0)demographics and measured baseline hormone levels for the 24 subjects receiving fe 999049 
 bmi body mass index , fsh follicle - stimulating hormone , e2 estradiol , lh luteinising hormone , prog progesterone , inhb inhibin b summary of subject characteristics demographics and measured baseline hormone levels for the 24 subjects receiving fe 999049 
 bmi body mass index , fsh follicle - stimulating hormone , e2 estradiol , lh luteinising hormone , prog progesterone , inhb inhibin b
the pharmacokinetic model was developed using nonlinear mixed - effects modelling , where both the population parameters , interindividual variability ( iiv ) , and residual errors are estimated . for parameters with iiv ,
the ith subject s individual parameter ( i ) is log - normally distributed:1\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \theta_{i } = \theta \exp \left ( { \eta_{i } } \right ) $ $ \end{document}i=expiwhere  is the typical population parameter and i is the individual random effect from an approximately normal distribution with mean zero and variance  for describing the iiv of the parameter . in the model ,
potential influential factors can be tested for significance as a covariate to explain some of the iiv in a parameter .
thus , the set of individual parameters ( i ) is given as a function of the typical population parameters (  ) individual values of the covariates ( ci ) and random effects ( i):2\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \vartheta_{i } = g\left ( { \vartheta , { \mathbf{c}}_{i } , \varvec{\eta}_{i } } \right ) $ $ \end{document}i = g,ci,ithe residual errors (  ) are assumed normally distributed with mean zero and variance  , and are the unexplained deviations of model predictions from the observations:3\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ y_{ij } = f\left ( {   t_{ij } , \vartheta_{i } } \right ) + \varepsilon_{ij } $ $ \end{document}yij = ftij,i+ijwhere yij is the observation at time tij and the subscript ij denotes the jth number for subject i. the individual model prediction f (  ) at time tij is calculated from the individual parameters ( i ) .
model development was guided by changes in the nonmem objective function value ( ofv ) , precision of parameter estimates , and graphical model goodness - of - fit assessments , including visual predictive checks ( vpc ) . in the vpc ,
the observed data are compared with model predictions based on 1000 simulated trial datasets using the final model .
it displays the observations and the 2.5th , 50th , and 97.5th percentiles of observations , and the 95  % confidence intervals ( ci ) for the corresponding model predictions are plotted .
the difference in ofv ( dofv ) between two nested models is approximately chi square - distributed , with degrees of freedom ( df ) equal to the difference in the number of parameters .
based on this , the statistical significance for inclusion / exclusion of a model parameter can be judged .
a significance level of 0.05 was used for discrimination among nested models and covariate testing , corresponding to a 3.84 change in ofv for 1  df .
measurable fsh levels before drug administration indicated that endogenous fsh was not fully suppressed in this trial .
the total fsh amount at time t in the central compartment was therefore modelled as the sum of the endogenous fsh and the exogenously administered rhfsh ( fe  999049):4\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { \text{fsh}}(t ) = { \text{fsh}}_{\text{en } } \left ( t \right )    { \text { + rhfsh}}\left ( t \right ) $ $ \end{document}fsh(t)=fshent+ rhfsht the endogenous fsh was assumed to have the same elimination rate constant from the central compartment as fe  999049 ; thus , the total fsh is eliminated from the central compartment with rate constant k given by clearance ( cl ) divided by volume of distribution ( v ) . since data are obtained after subcutaneous dosing , the bioavailability ( f ) is not known .
the cl and v estimated here are therefore the apparent cl ( cl / f ) and apparent v ( v / f ) .
the change in endogenous fsh amount over time in the central compartment was described as a turnover model with a zero - order production rate kin:5\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \frac{{{\text{dfsh}}_{\text{en } } \left ( t \right)}}{{{\text{d}}t } } = k_{\text{in } } - k \cdot { \text{fsh}}_{\text{en } } \left ( t \right ) $ $ \end{document}dfshentdt = kin - kfshent at time zero , before rhfsh dosing , the measured concentration in the central compartment is purely endogenous fsh .
the endogenous fsh baseline concentrations ( fshbl ) for each subject can therefore be estimated as the initial concentration in the central compartment and is the initial condition for the differential equation for the central compartment fshblv / f , i.e. the individual endogenous fsh amount at baseline .
when assuming steady state and hence zero change in the total fsh concentration before dosing , individual values of kin can be calculated from the other estimated parameters and baseline values .
one- and two - compartment models were tested for describing the distribution of fe  999049 .
a lag - time and transit compartments were investigated for modelling a potential delay in the absorption process .
body weight , age , and hormone ( estradiol , inhibin b , progesterone ) values at baseline were tested as potential covariates , explaining some of the iiv in parameter estimates .
inhibin b was also tested as a covariate over time since it has a purely inhibitory effect on fsh .
the measured hormone levels were assumed to not affect the rhfsh concentration , but only tested as covariates on the endogenous fsh baseline parameter ( fshbl ) . besides a decrease in ofv , the significance of a covariate was also evaluated by looking at the reduction in the iiv , measured as a coefficient of variation ( cv ) for the parameter s random effect .
the models were implemented and parameters estimated in nonmem 7.2.0 ( icon development solutions , ellicott city , md , usa ) . data handling and
the fe  999049 absorption was found to be delayed , and a one - compartment transit model best described the process .
endogenous fsh contributed significantly to the total fsh concentration in the central compartment . observed inhibin b levels
[ inhb(t ) ] suppressed endogenous fsh production rate ( kin ) over time when introduced in the model as a time - varying covariate .
the relationship was found to be best described by an inhibitory emax model , also know as an imax function , expressed as a fractional inhibition with complete inhibition possible and a parameter inhbef as the inhibin b concentration giving 50  % inhibition [ see eq .
the full compartment model is illustrated in fig .  1 and the corresponding differential equations are as follows : 6\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				
\usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \frac{{{\text{drhfsh}}_{\text{ds } } \left ( t \right)}}{{{\text{d}}t } } = - k_{\text{tr } } \cdot { \text{rhfsh}}_{\text{ds } } \left ( t \right ) $ $ \end{document}drhfshdstdt =- ktrrhfshdst7\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \frac{{{\text{drhfsh}}_{\text{tr } } \left ( t \right)}}{{{\text{d}}t } } = k_{\text{tr } } \cdot { \text{rhfsh}}_{\text{ds } } \left ( t \right ) - k_{a } \cdot { \text{rhfsh}}_{\text{tr } } \left ( t \right ) $ $ \end{document}drhfshtrtdt = ktrrhfshdst - karhfshtrt8\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \frac{{{\text{dfsh}}\left ( t \right)}}{{{\text{d}}t } } = k_{\text{in } } \cdot \left ( { 1 - \frac{{{\text{inhb}}(t)}}{{{\text{inhb}}(t ) + { \text{inhb}}_{\text{ef } }   } } } \right ) + k_{a } \cdot { \text{rhfsh}}_{\text{tr } } \left ( t \right ) - k \cdot { \text{fsh}}\left ( t \right ) $ $ \end{document}dfshtdt = kin1-inhb(t)inhb(t)+inhbef+karhfshtrt - kfshtwhere rhfshds(t ) and rhfshtr(t ) are the fe  999049 amounts at the dosing site and in the transit compartment , respectively , and ktr is the absorption rate from the dosing site to the transit compartment
. equation  ( 8) describes the change in total fsh amount in the central compartment , with the endogenous fsh contribution as the first term .
1compartment diagram illustrating the pharmacokinetic model for fe  999049 with a contribution of endogenous fsh to the central compartment .
fsh follicle - stimulating hormone , rhfsh recombinant human fsh , fsh 
 en 
 ( t ) endogenous fsh amount , inhb(t ) inhibin b level , rhfsh 
 tr 
 ( t ) rhfsh amount in the transit compartment , fsh(t ) total fsh amount in the central compartment , k 
 in endogenous fsh production rate , k 
 tr absorption rate from the dosing site , k 
 a absorption rate from the transit compartment , k elimination rate compartment diagram illustrating the pharmacokinetic model for fe  999049 with a contribution of endogenous fsh to the central compartment .
fsh follicle - stimulating hormone , rhfsh recombinant human fsh , fsh 
 en 
 ( t ) endogenous fsh amount , inhb(t ) inhibin b level , rhfsh 
 tr 
 ( t ) rhfsh amount in the transit compartment , fsh(t ) total fsh amount in the central compartment , k 
 in endogenous fsh production rate , k 
 tr absorption rate from the dosing site , k 
 a absorption rate from the transit compartment , k elimination rate iiv was introduced for cl / f , v / f , ktr , and fshbl .
the variances of the iiv on cl / f and v / f were positively correlated .
body weight was an allometrically scaled covariate at cl / f and v / f , with the power exponents fixed to 0.75 and 1 , respectively . adding body weight as a covariate reduced cv for the unexplained iiv from 18.1 to 15.6  % cv for cl / f , and from 22.0 to 18.4  % cv for v / f .
a trend towards lower estimated fsh baseline levels was seen in individuals with higher progesterone baseline levels ( fig .  2 ) .
the progesterone baseline effect was confirmed to be a statistically significant covariate when modelled as an inhibitory power function at the fshbl parameter .
points are individual - predicted fshbl values and observed progesterone baseline values with a smooth lowess trend line ( broken line ) .
the solid line is the power function used in the model describing the typical population relationship for the effect of progesterone baseline at the parameter fshbl .
fsh follicle - stimulating hormone , fsh 
 b1 endogenous fsh baseline the relationship between endogenous fsh and progesterone at baseline .
points are individual - predicted fshbl values and observed progesterone baseline values with a smooth lowess trend line ( broken line ) .
the solid line is the power function used in the model describing the typical population relationship for the effect of progesterone baseline at the parameter fshbl .
fsh follicle - stimulating hormone , fsh 
 b1 endogenous fsh baseline a combined additive and proportional error model was used to describe the residual error .
the final model parameter estimates are displayed in table  2 . the vpc plot ( fig . 
3 ) based on 1000 simulations , shows that the model overall adequately predicts the fsh concentrations over time , with appropriate variation .
the model also satisfactorily predicted the individual pharmacokinetic profile for each subject in the trial ( fig . 
the typical population cl / f and v / f for a 65-kg woman were 0.423  l / h and 24.3  l , respectively , giving an elimination rate of 0.017  h and a terminal half - life of 40  h. when the best model was found in a forward development process , it was checked to see if all the elements were still significant or if the model could be reduced .
no reduction was possible since the dofv was a minimum of 7 and up to 78.table  2pharmacokinetic parameter estimatesparameterestimate(rse%)iiv cv%(rse%)shrinkage ( % ) cl / f ( l / h)0.423(3.9)15.6(12.7)1.05v / f ( l)24.30(4.6)18.4(14.4)2.80 
 k 
 tr ( h)0.329(17.0)83.4(24.1)24.6 
 k 
 a ( h)0.148(13.2)fshbl ( g / l)0.162(9.1)27.8(16.5)5.51progblef 
 0.246(34.7)inhbef 
 100(37.2)typical population parameter estimates obtained from modellingfor cl / f and v / f the value is the typical value for a woman weighing 65  kg 
 rse relative standard error , iiv interindividual variability , cv coefficient of variation , f bioavailability , cl / f apparent clearance , v / f apparent volume of distribution , k 
 tr absorption rate from the dosing site , k 
 a absorption rate from the transit compartment , fsh 
 bl endogenous follicle - stimulating hormone baseline , progbl 
 ef power exponent for progesterone baseline covariate effect , inhb 
 ef inhibin b time - varying covariate effectfig
.  3visual predictive check for the final model , showing the individual observed fsh concentrations ( points ) and the 2.5th , 50th , and 97.5th percentiles of observations ( lines ) .
the shaded areas are the 95  % confidence intervals for the 2.5th , 50th , and 97.5th percentiles of the simulations .
4individual model fit with observed ( points ) and predicted ( lines ) total fsh concentrations versus time .
fsh follicle - stimulating hormonetable  3change in objective function value for reduced modelsremoving 
 df 
 dofvwt at v / f 
 014.36wt at cl / f 
 012.61wt at cl / f and v / f 
 08.91cov ( cl / f , v / f)111.77progesterone effect17.01inhibin b effect178.01the resulting dofv when removing covariates or the correlation between cl / f and v / f [ cov(cl / f , v / f ) ] 
 f bioavailability , cl / f apparent clearance , v / f apparent volume of distribution , wt body weight , df degrees of freedom , dofv difference in objective function value pharmacokinetic parameter estimates typical population parameter estimates obtained from modelling for cl / f and v / f the value is the typical value for a woman weighing 65  kg 
 rse relative standard error , iiv interindividual variability , cv coefficient of variation , f bioavailability , cl / f apparent clearance , v / f apparent volume of distribution , k 
 tr absorption rate from the dosing site , k 
 a absorption rate from the transit compartment , fsh 
 bl endogenous follicle - stimulating hormone baseline , progbl 
 ef power exponent for progesterone baseline covariate effect , inhb 
 ef inhibin b time - varying covariate effect visual predictive check for the final model , showing the individual observed fsh concentrations ( points ) and the 2.5th , 50th , and 97.5th percentiles of observations ( lines ) .
the shaded areas are the 95  % confidence intervals for the 2.5th , 50th , and 97.5th percentiles of the simulations .
fsh follicle - stimulating hormone individual model fit with observed ( points ) and predicted ( lines ) total fsh concentrations versus time .
fsh follicle - stimulating hormone change in objective function value for reduced models the resulting dofv when removing covariates or the correlation between cl / f and v / f [ cov(cl / f , v / f ) ] 
 f bioavailability , cl / f apparent clearance , v / f apparent volume of distribution , wt body weight , df degrees of freedom , dofv difference in objective function value in order to see what difference accounting for endogenous fsh in the model makes , the individual concentration 
time profiles for the model - predicted endogenous fsh level were compared with when the observed endogenous fsh baseline was assumed to be constant throughout the trial ( fig . 
furthermore , accounting for inhibin b suppression of the endogenous fsh contribution over time increased the v / f estimate from 18.9 to 24.3  l. not including inhibin b as a covariate in the final model resulted in an increase from 83.4 to 164.9  % cv at ktr.fig . 
the broken blue line indicates the observed endogenous fsh baseline level when assuming it to be constant throughout the trial , and the solid blue line is the model - predicted endogenous fsh level when accounting for the suppression caused by the observed inhibin b levels ( purple line ) over time , as identified by the model .
the broken blue line indicates the observed endogenous fsh baseline level when assuming it to be constant throughout the trial , and the solid blue line is the model - predicted endogenous fsh level when accounting for the suppression caused by the observed inhibin b levels ( purple line ) over time , as identified by the model .
the fe  999049 population pharmacokinetics have previously been characterised after single - dose administration using a one - compartment model .
in accordance with those results , the current analysis showed that after multiple fe  999049 dosing , the population pharmacokinetics were also best described by a one - compartment model with first - order elimination and absorption through a transit compartment .
a correlation between cl / f and v / f was observed , and body weight was an allometrically scaled covariate for both cl / f and v / f .
hereby , it was confirmed that lower average steady - state drug exposure in females correlates with higher body weight after repeated doses , which is in agreement with single - dose fe  999049 results .
in contrast to the single - dose trial , fsh downregulation in this trial was not complete since a measurable endogenous fsh level was detected before any drug administration . the observed endogenous fsh contributes to total fsh measurements and thereby affects the evaluation of rhfsh concentrations in the trial
; thus , it had to be accounted for in the model . to adjust the model accordingly , it was necessary to add an endogenous supply of fsh to the central compartment for proper prediction of the total fsh concentration . when evaluating the pharmacokinetics of drugs that are naturally occurring substances , data are often baseline corrected in order to get values representing only the exogenously administered drug
however , when dealing with hormones , fluctuations from the endogenous baseline value over time may be observed due to various feedback control systems .
even with gonadotropin downregulation , the endogenous fsh level will most likely change over time as progesterone , inhibin b , and estradiol exert inhibitory and stimulatory feedback loops at the production and release of gonadotropins .
the potential impact of varying endogenous fsh levels at the pharmacokinetic assessment of fe  999049 was tested by incorporating different covariate relationships in the model . before any drug was administered
, there was a tendency towards a decreased endogenous fsh baseline with increasing progesterone baseline .
a power function for the relation between progesterone baseline and fshbl was shown to be significantly better than using a linear , an exponential , or an imax function . when introduced as a time - varying covariate
, inhibin b was found to suppress the endogenous fsh production rate over time , with an imax function being the most significant relationship .
when accounting for the inhibitory effect of inhibin b over time , the resulting model - predicted endogenous fsh concentration profiles indicated that the observed fsh baseline is an overestimate of the endogenous fsh level at all time points other than zero ( fig . 
, the drug concentration would be underestimated if the data were baseline corrected . when multiple covariates and correlations have been added to the model , it is possible that initial significant relations have become redundant and hence can be removed from the model .
it was therefore checked whether the model could be reduced . removing any relations resulted in an increase in ofv , and
removing body weight from both cl / f and v / f increased ofv less than removing body weight from only one of these parameters , which could be due to the fact that there is a correlation between the two parameters .
the population pharmacokinetic parameters estimated with nonlinear mixed - effect modelling for single and multiple doses of fe  999049 are similar . in the modelling of single - dose data ,
cl / f and v / f for a 65  kg woman were found to be 0.430  l / h and 28.0  l , respectively , with a resulting terminal half - life of 45  h. in this trial , cl / f was 0.423  l / h and v / f was 24.3  l. by incorporating the inhibin b dynamics into the model , and suppressing the endogenous fsh contribution instead of letting the fsh baseline be constant throughout the trial , v / f increased from 18.9 to 24.3  l , suggesting that without proper estimation of endogenous fsh levels over time a bias in v / f may be obtained .
the terminal half - life in the current analysis was estimated to be 40  h , thus slightly shorter than the 45  h reported after single - dose administration .
the phenomenon of shorter half - life after repeated subcutaneous dosing compared with single administration was also observed using nca - based methods for analysing different fe  999049 datasets , but could not be further explained . the cl / f for other rfsh products has been found to be 0.99  l / h   and 0.71  l / h ; thus , our estimated cl / f for fe  999049 is consistently lower .
furthermore , the terminal half - life of fe  999049 is longer than reported values in the literature for repeated subcutaneous dosing of rfsh products , which is approximately 24  h [ 3 , 10 ] . in the same studies , a longer terminal half - life has been reported after single subcutaneous dosing ( 34.5   and 27  h ) .
one study   found a terminal half - life of 59  h after single subcutaneous dosing of 115  iu rfsh ; however , they noted that the estimate was not reliable since the low dose only gave a concentration slightly over baseline .
this work shows that when subjects have measurable endogenous fsh levels , total fsh concentration is affected , and the endogenous fsh level varies over time as a consequence of the feedback from ovarian hormones . if not accounting for this in the model , the pharmacokinetic parameter estimates can be biased
. the inclusion of endogenous hormones in the model might be even more important in the target patient group since they are likely to have higher endogenous fsh baseline levels as they will not be pituitary downregulated or suppressed .
in addition , the variation in endogenous fsh levels over time could be very different between patients due to their different types of infertility , which might affect the reproductive hormone dynamics differently .
characterising the endogenous fsh level and considering the hormone dynamics is therefore important as it could have clinical relevance for the dose needed . in summary , this early characterisation of fe  999049 population pharmacokinetics after repeated dosing is in line with previously reported findings of rhfsh population pharmacokinetics after single - dose administration .
the results provide a basis for study design and data evaluation considerations in the future development of rfsh products .
the standard method of baseline - correcting data for rfsh products has been questioned since the model indicates that endogenous fsh significantly affects the total fsh concentration and varies over time .
the multiple - dose fe  999049 population pharmacokinetics were in agreement with results obtained after single - dose administration .
it was confirmed that after repeated drug administration , fe  999049 exposure also appears lower in females with a higher body weight .
these findings add to the evidence of body weight being an important factor in designing individualised dosing schemes .
this analysis indicates that the inclusion of endogenous fsh levels and quantifying the influence of the reproductive hormone dynamics at the endogenous fsh level over time may well be essential for proper estimation of the pharmacokinetic parameters for fsh preparations .
the standard method of baseline - correcting data does not account for a variation during the trial and could potentially cause an underestimation of serum drug concentration and hence inaccurate parameter estimation . to better explore and link the time - varying impact of inhibin b on endogenous fsh levels
, a more mechanistic modelling approach may be warranted for including indirect delayed response and hormone feedback mechanisms in a pharmacokinetic 
the relevance of endogenous hormone levels should also be related to other pharmacodynamic and clinical endpoints in order to establish its significance .
this work was supported by innovation fund denmark ( industrial phd case number 11 - 117436 ) .
trine hyer rose , daniel rshammar , lars erichsen , and lars grundemar are all current or former employees of ferring pharmaceuticals a / s .
all procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 helsinki declaration and its later amendments or comparable ethical standards .
 objectivethe aim of this study was to characterise the population pharmacokinetics of fe 999049 , a novel recombinant human follicle - stimulating hormone ( fsh ) , after multiple dosing in healthy women , taking into account endogenous fsh levels and the reproductive hormone dynamics.methodslongitudinal measurements of fsh , luteinising hormone , progesterone , estradiol , and inhibin b levels were collected after repeated subcutaneous dosing with 225  iu of fe  999049 in 24 gonadotropin downregulated healthy women . 
 the fsh data were described using nonlinear mixed - effects modelling.resultsthe measured fsh levels were modelled as a sum of endogenous fsh and fe  999049 . 
 the fe  999049 population pharmacokinetics were best described using a one - compartment model with first - order absorption and elimination , and a transit model for delayed absorption . 
 the apparent clearance and volume of distribution increased with body weight in accordance with an allometrically scaled power exponent of 0.75 and 1 , respectively . 
 endogenous fsh levels were lower in individuals with higher progesterone levels at baseline and were further suppressed over time with increasing inhibin b levels.conclusionsthis characterisation of fe  999049 population pharmacokinetics after repeated dosing is in line with previous findings after single - dose administration . 
 the results provide a basis for study design and data evaluation in the future development of recombinant fsh products , and show it can be of importance to account for endogenous fsh levels and its variation over time for accurate estimation of exogenously administered fsh pharmacokinetic parameters . 
 thus , correcting fsh concentrations by the observed endogenous fsh baseline value at all time points may be incorrect . 
Key Points
Introduction
Methods
Clinical Trial Design and Data
Pharmacokinetic Modelling
Software
Results
Discussion
Conclusions
Funding
Conflict of interest
Ethical approval
Informed consent