Population Pharmacokinetics of Semaglutide for Type 2 Diabetes.

INTRODUCTION: The aim of the present analysis was to characterise the absorption, distribution and elimination of semaglutide by means of population pharmacokinetic (PK) models using data from nine clinical pharmacology trials conducted in both healthy subjects and those with type 2 diabetes.
METHODS: Data were obtained from trials with subcutaneous and intravenous administration of semaglutide that utilised frequent PK sampling and included a total of 353 subjects with 10,573 concentration values.
RESULTS: Semaglutide PK properties across trials, drug product strengths and populations were well characterised by a two-compartment model with first-order absorption and elimination. For a typical subject with type 2 diabetes, clearance was estimated to be 0.0348 L/h [95% confidence interval (CI) 0.0327-0.0369 L/h], and the central and peripheral volumes were estimated to be 3.59 L (95% CI 3.28-3.90 L) and 4.10 L (95% CI 3.78-4.42 L), respectively (i.e. a total volume of distribution of 7.7 L). Interindividual variation was low (~ 15%) for both clearance and volumes of distribution, with low residual error (< 5%). Clearance and the total volume of distribution were approximately proportional to body weight. Minor differences were identified between healthy subjects and subjects with type 2 diabetes with respect to clearance and absorption rate, and between injection sites with respect to bioavailability.
CONCLUSIONS: A novel two-compartment model was developed to provide the general characteristics of semaglutide absorption following subcutaneous administration, and of distribution and elimination across administration routes. Semaglutide PK was shown to be predictable across populations and administration routes and within subjects, and was primarily influenced by body weight. FUNDING: Novo Nordisk, Bagsværd, Denmark.

Introduction

Semaglutide is a receptor agonist of the naturally occurring peptide hormone glucagon-like peptide-1 (GLP-1). As with native GLP-1, semaglutide acts to stimulate insulin secretion and inhibit glucagon release depending on the blood glucose concentration, leading to improved blood glucose levels together with a reduced risk of hypoglycaemia [1, 2]. Semaglutide is approved for the treatment of type 2 diabetes mellitus, using the subcutaneous route of administration. With 94% amino acid sequence homology to native GLP-1, semaglutide has three structural modifications that make it less susceptible to degradation by dipeptidyl peptidase-4 (DPP-4) enzymes and prolong its half-life to approximately 1 week [1, 3]. Semaglutide is thereby appropriate for once-weekly administration.

The prolonged exposure of semaglutide as compared with native GLP-1 is mainly driven by slow systemic elimination and, to a lesser extent, by delayed absorption [4], in line with the observation that more than 99% of the semaglutide is bound to plasma albumin [5]. In a metabolism trial, semaglutide was the primary component circulating in plasma and, prior to excretion, it was metabolised by proteolytic cleavage of the peptide backbone and sequential beta-oxidation of the fatty acid side chain [4]. Semaglutide metabolites are primarily excreted in the urine and faeces [4].

Semaglutide population pharmacokinetics (PK) and exposure in terms of average concentration in the dosing interval have been previously investigated for subcutaneous dosing using phase 3 steady-state data from subjects with type 2 diabetes [6]. Of the covariates investigated, only body weight had a relevant effect on the exposure of semaglutide at either once-weekly maintenance dose of 0.5 or 1.0 mg [6].

In addition to an assessment of semaglutide exposure levels, the aim of the present analysis was to provide additional details of the PK properties of semaglutide in terms of absorption, distribution and elimination and their relationships to covariate factors, in order to develop a novel disposition model that would be generally applicable across different routes of administration. This was achieved by means of population PK models using frequently sampled PK profiles from nine clinical pharmacology trials. Inclusion of data from intravenous administration provided information about the structural two-compartment disposition model, allowing for separate estimation of clearance and bioavailability.

The analysis was applied to identify and validate the use of a one-compartment model previously applied to analyse sparsely sampled phase 3 data [6], and subsequently to investigate the consistency between PK data from phase 1 and phase 3 trials.

Methods

Clinical Data Included in the Analysis

A two-compartment PK model was estimated using data from nine clinical pharmacology trials employing single as well as multiple semaglutide doses and conducted in subjects with type 2 diabetes or obesity, or in healthy volunteers (Table 1). The included trials were those with administration of recombinant semaglutide in which PK samples, assessed via liquid chromatography and tandem mass spectrometry (LC–MS/MS) [3], were available at the time of the analysis. One trial provided data for the intravenous administration of semaglutide and also used different drug product strengths (see further descriptions in the “Methods” section of the Electronic supplementary material, ESM).

Table 1

Overview of clinical trials included in the analysis

Trial [reference]	Description	Number of subjects	Population	PK data	Number of samples/subject	Doses	Drug product strength (mg/mL)
1 [18]*	Japanese safety and PK	32	Healthy male Japanese and Caucasians	SS + SD profiles	46	0.25 mg > 0.5 mg > 1.0 mg	1.34
2 [19]*	Beta-cell function	37	Males and females with T2D	SS profiles	24	0.25 mg > 0.5 mg > 1.0 mg	1.34
3 [5]	Hepatic impairment	44	Healthy males and females and those with hepatic impairment	SD profiles	20	0.5 mg	1.34
4 [20]*	QT study	82	Healthy males and females	Trough samples + SS profiles	50	0.25 mg > 0.5 mg > 1.0 mg > 1.5 mg	1.34
5 [21]*	Hypoglycaemic counterregulation	37	Males and females with T2D	Trough samples + SS profiles	24	0.25 mg > 0.5 mg > 1.0 mg	1.34
6 [22]*	Energy intake	30	Males and females with obesity	SS profiles	16	0.25 mg > 0.5 mg > 1.0 mg	1.34
7a (data on file)	Drug product strengths (equivalence)	42	Healthy males and females	Two SD profiles	20 (s.c.) 37 (i.v.)	0.5 mg s.c. 0.25 mg i.v.	1, 3 and 10
8 [23]*	DDI study 1	23	Healthy males and females	Trough samples + SS profiles	27	0.25 mg > 0.5 mg > 1.0 mg	1.34
9 [23]*	DDI study 2	26	Healthy males and females	Trough samples + SS profiles	27	0.25 mg > 0.5 mg > 1.0 mg	1.34

All trials were used for the two-compartment model, and trials marked with an * were used for the one-compartment model

Doses were administered using the PDS290 injection pen, except for trial 7, which used NovoPen 4 for s.c. administration and a syringe for i.v. administration

ClinicalTrials.gov identifiers for the above trials: NCT02146079; NCT02212067; NCT02210871; NCT02064348; NCT02147431; NCT02079870; NCT02231684; NCT02022254; and NCT02243098

DDI drug–drug interaction, i.v. intravenous dosing, PK pharmacokinetic(s), s.c. subcutaneous dosing, SD single-dose, SS steady-state, T2D type 2 diabetes

aData on file, partly reported in the present publication

In order to compare the current analysis with the analysis conducted using phase 3 data, a one-compartment model was estimated using data from seven of the clinical pharmacology trials (Table 1) which employed multiple doses of the semaglutide drug product used in the phase 3 trials. Data from the two trials with single dose administration were not included in this comparison.

Population Pharmacokinetic Modelling

Datasets containing plasma concentration and covariate data from each trial were converted to a format compatible with the nonlinear mixed-effects modelling (NONMEM) software (ICON Development Solutions, Ellicott City, MD, USA), version 7.1.2. Out of a total of 10,573 concentration values, 43 values (i.e. < 1%) were excluded because the subjects had an inadequate dosing history. Conditional weighted residuals (CWRES) were used as a diagnostic tool to identify outliers, and an additional 24 concentration values (23 in the one-compartment model) were excluded from the final dataset as they were associated with CWRES greater than 6 in order to limit bias from influential data items. The software package R (version 3.2.3, obtained from the R Foundation: http://www.r-project.org/foundation) was used for data file processing, explorative data analysis and to plot figures.

Model Development

Model development was conducted in NONMEM as follows:

The structural disposition model, which was supported by graphical analysis, in particular of intravenous data, was a two-compartment model with first-order elimination. Different absorption models were investigated, including parallel absorption, but a simpler first-order absorption model for subcutaneous administration was found to adequately describe the data, and no significant improvement was obtained based on other absorption models.

An additive residual error model on log-transformed concentration values was assessed to adequately describe the data.

Inter-individual variability in the absorption rate (ka), clearance (CL), intercompartmental clearance (Q), central volume (Vc), and peripheral volume (Vp) was explored using a log normal distribution in a model that included covariates with clear effects: drug product strength on the absorption rate constant and body weight as a covariate factor on CL, Vc and Vp. A correlation between the central volume and CL was significant and included in the model. Due to the limited amount of intravenous data, less information was available for individual estimation of Q and for the separation of Vc and Vp. Various implementations for these parameters were investigated and selected based on objective function values; i.e. using individual (eta) values for CL to identify differences across subjects for Q, as well as for Vc and Vp.

The covariate effects of sex, race, ethnicity, age, hepatic impairment and glycaemic status were investigated for CL, Vc and ka based on a graphical analysis of individual parameter estimates versus body weight, and glycaemic status was identified as a covariate factor for CL and ka. Similarly, injection site was investigated with respect to F and ka. Covariates with 90% confidence intervals which were within the 80–125% equivalence interval were not considered important.

Covariate effects were investigated using a full model approach with inclusion of all investigated covariates in one step (Table 2). The final model included only those covariates that were significant (10% alpha level). To assess the amount of variability explained by covariate factors, a base model was investigated, including only covariates for the effects of drug product strength on the absorption rate constant.

Table 2

Covariates included in the full two-compartment model

PK Parameter	Covariates includeda
Clearance (CL)	Body weight (continuous) Sex (female, male) Race (white, black, Asian) Ethnicity (Hispanic, non-Hispanic) Glycaemic status (normoglycaemia, T2D) Age group (≤ 65, > 65 years) Hepatic impairment
Bioavailability (F)	Injection site (abdomen, thigh)
Intercompartmental clearance (Q)	Body weight
Central volume (Vc)	Body weight
Peripheral volume (Vp)	Body weight
Absorption rate constant (ka)	Drug product strength (1, 1.34, 3 and 10 mg/mL) Glycaemic status (normoglycaemia, T2D)

T2D type 2 diabetes

aAssessed at baseline

Model Description

The full model for subcutaneous administration was a model with first-order absorption, two-compartment disposition and first-order elimination, whereas intravenous data were modelled as a bolus injection into the central compartment. Inter-individual variability was estimated for ka, and CL/Q (with a common factor), and for Vc/Vp (with a common factor).

The covariate models were parameterised as follows for the ith individual: and the effect of drug product strength on the rate of absorption was described as

Covariate relationships were described as follows:

In the above equations, θ represents an estimated value for a covariate and η represents between-subject variability. The other exponents are indicator variables assigned a value of 1 for the specific covariate values and 0 otherwise. CLref, Qref, Vc,ref, Vp,ref, Fref and ka, ref are parameter values for the reference subject profile, which is a healthy, white, non-Hispanic female aged ≤ 65 years with a body weight of 85 kg, using the abdomen as the injection site and dosed with a 1.34 mg/mL semaglutide drug product. The residual error model was additive for log-transformed observations.

Steady-State Exposure Evaluation for Subcutaneous Semaglutide

A one-compartment model, similar to the one used to evaluate phase 3 data, was investigated using steady-state data from subcutaneous dosing only in clinical pharmacology trials; hence, the fraction absorbed (F) could not be estimated separately in this model. Therefore, the model was parameterised by the apparent clearance (CL/F), volume of distribution (V/F) and the absorption rate constant (ka). The inter-individual variability was estimated for each of these. The residual error model was additive for log-transformed observations.

The one-compartment model contained three covariate relationships: body weight and glycaemic status (healthy, type 2 diabetes) on CL, and body weight on volume of distribution. These were parameterised as for the two-compartment model.

Compliance with Ethics Guidelines

This article is based on previously conducted studies and does not contain any studies with human participants or animals performed by any of the authors.

Results

Population Characteristics

A total of 353 subjects with 10,573 semaglutide concentrations were included in the population PK model analysis. For the comparison to the phase 3 modelling approach, 267 subjects contributed with 8292 semaglutide concentration values. The demographic characteristics for the subjects included in each of the analyses are summarised in Table 3.

Table 3

Demographic characteristics of subjects included in the two- and one-compartment models

Category	Group	Total 2-compartment	Total 1-compartment
All	N	353	267
Population	Normoglycaemia	277	193
	T2D	76	74
Race	White	327	243
	Asian	16	16
	Other/missing	7	5
	Black/African American	3	3
Ethnicity	Not Hispanic	349	264
	Hispanic	4	3
Sex	Male	226	180
	Female	127	87
Dosing	Multiple dose	267	267
	Single dose	86	–
Dose, mg	0.25	5	5
	0.5	105	19
	1.0	166	166
	1.5	77	77
Injection site	Abdomen	286	244
	Thigh	67	23
Age, years	Mean (SD)	44.6 (11.8)	44.1 (11.5)
	Range	[19–70]	[19–70]
BMI, kg/m2	Mean (SD)	26.9 (4.3)	27.0 (4.3)
	Range	[18.7–42.8]	[20.0–42.8]
Body weight, kg	Mean (SD)	81.9 (15.1)	83.0 (14.7)
	Range	[51.9–121.2]	[54.1–121.2]
Hepatic impairment status	Mild, moderate or severe	25	25
	No impairment	19	19

BMI body mass index, N number of subjects, SD standard deviation, T2D type 2 diabetes

Semaglutide Pharmacokinetics

The final population PK model provided an accurate description of the PK time course across the nine clinical pharmacology trials. Figure 1 shows the consistency between the model and data across the eight trials which employed both single and multiple doses of marketed semaglutide with a drug product strength of 1.34 mg/mL. Figure 2 shows that the biphasic PK profile obtained from the two-compartment model is consistent with intravenous data and, furthermore, that the rate of absorption as well as the maximum concentration following subcutaneous doses of different strengths were well captured by the model.

Fig. 1

Observed and model-predicted semaglutide pharmacokinetic profiles from eight clinical pharmacology trials. All trials employed semaglutide doses of up to 1.0 mg, except trial 3, which was a single dose (0.5 mg) trial, and trial 4, which used a supratherapeutic dose of 1.5 mg. Data points with error bars are observed geometric means and 95% CIs. Model predictions (lines) are population predictions from the final two-compartment model of semaglutide PK. CI confidence interval, PK pharmacokinetics, SD single dose, SS steady state

Fig. 2a–b

Observed and model-predicted subcutaneous and intravenous semaglutide PK profiles over 7 days (a) and 30 days (b). Data are from trial 7 with semaglutide doses of up to 0.5 mg. Data points with error bars are observed geometric means and 95% CIs. Model predictions (lines) are population predictions from the final two-compartment model of semaglutide PK. In b, two values above the lower limit of quantification obtained at day 35 were not included in the figure. CI confidence interval, i.v. intravenously, PK pharmacokinetic(s), s.c subcutaneously

The model-based analysis suggested that drug product strength influences the rate of absorption but not the bioavailability (Table 4), indicating that the total exposure of semaglutide is not affected by the product strength. These model-based results are consistent with the noncompartmental results for trial 7 (see Table S1 in the ESM), in which equivalence was demonstrated for exposure (AUC0–∞) for all pairwise comparisons between the three strengths of semaglutide. For Cmax, the comparison between 1 mg/mL and 3 mg/mL met the equivalence criterion but equivalence was not shown for 1 mg/mL relative to 10 mg/mL or for 3 mg/mL relative to 10 mg/mL, as Cmax increased with increasing strength.

Table 4

Parameter estimates obtained from the final two-compartment model

Parameter	Estimate	[95% CI]	RSE (% CV)	IIV (% CV)	Shrinkage (%)
Absorption rate constant (ka), h−1	0.0253	[0.0236–0.027]	3.42	37.9	9.4157
Clearance (CL), L/h	0.0348	[0.0327–0.0369]	3.06	15.2	1.5202
Central volume (Vc), L	3.59	[3.28–3.9]	4.44	15.4	6.6167
Intercompartmental clearance (Q), L/h	0.304	[0.249–0.359]	9.19	15.2	1.5202
Peripheral volume (Vp), L	4.10	[3.78–4.42]	3.97	15.4	6.6167
Absolute bioavailability (F)	0.847	[0.798–0.896]	2.96	NA	NA
Body weight effect on CL and Q	1.01	[0.912–1.1]	4.71	NA	NA
Body weight effect on Vc and Vp	0.923	[0.833–1.01]	4.96	NA	NA
T2D effect on CL	1.12	[1.09–1.16]	1.66	NA	NA
Injection site effect on F	0.883	[0.844–0.921]	2.24	NA	NA
Drug product strength effects on ka, 1 mg	0.0346	[0.0301–0.039]	6.53	NA	NA
Drug product strength effects on ka, 3 mg	0.0526	[0.0438–0.0615]	8.58	NA	NA
Drug product strength effects on ka, 10 mg	0.139	[0.109–0.168]	11.0	NA	NA
T2D effect on ka	0.544	[0.485–0.603]	5.49	NA	NA
Residual error	0.103	NA	NA	NA	4.2914

A covariance of 0.0172 was found between CL and Vc. Data from all 9 trials were included in the model

CI confidence interval, CL clearance, CV coefficient of variation, F bioavailability, IIV interindividual variability, NA not applicable, Q intercompartmental clearance, RSE relative standard error, T2D type 2 diabetes, V central volume, V peripheral volume

All parameters were successfully estimated (relative standard error ≤ 11%) in the final model using data from the nine clinical trials (Table 4). In accordance with findings based on phase 3 semaglutide data investigating covariate effects on CL/F [6], body weight was found to be the only important covariate with respect to CL and, hence, average semaglutide concentrations under steady-state conditions with 90% CIs outside the equivalence interval. There were also clear relationships between body weight and the Vc and Vp, whereas the ka was independent of body weight (Table 4 and Fig. S1c in the ESM). Less important though statistically significant effects included an effect on semaglutide CL of glycaemic status, with a 12% higher CL (95% CI 9–16%) observed in subjects with type 2 diabetes compared to healthy subjects (Table 4 and Fig. S1a in the ESM). Moreover, ka was reduced in subjects with type 2 diabetes compared to healthy subjects (Table 4 and Supplementary Fig. S1c in the ESM). For a typical subject with type 2 diabetes and a body weight of 85 kg, clearance was estimated to be 0.0348 L/h (95% CI 0.0327–0.0369 L/h) and the total volume of distribution was estimated to be 7.7 L (Vc + Vp; Table 4). A small effect of injection site was found with respect to the bioavailability, which was 12% lower when using the thigh compared with the abdomen (95% CI 0.84–0.92). The effects of other covariates were not relevant as they were found to lie within the 80–125% equivalence interval. No effects of sex, race, ethnicity, age, or hepatic impairment on CL were found, and these covariates were excluded from the final model.

Trial 7, which provided data from both subcutaneous and intravenous administration, allowed the absolute bioavailability for subcutaneous administration of semaglutide to be estimated. Using noncompartmental analysis, the absolute bioavailability after the subcutaneous administration of semaglutide was estimated to be 89% (95% CI 0.83–0.94). This was captured well in the final model, which estimated the bioavailability to be 85% (95% CI 0.80–0.90) (Table 4).

Inclusion of covariates in the final two-compartment model explained 58%, 56% and 24% of the interindividual variance of CL, volume of distribution and ka, respectively.

Model Qualification

Parameter estimates (Table 4) were estimated with a low (≤ 11%) relative standard error in the final model; shrinkage values were all below 10%, supporting the accuracy of the individual parameter estimates.

The standard visual predictive checks indicated that the model was able to reproduce the data used for estimation in terms of both median concentrations and variability between subjects (Fig. S2 in the ESM). In addition, an illustration of variability in a dosing interval was included for the dose-normalised steady-state concentrations from the target semaglutide doses 0.5 and 1.0 mg (Fig. 3a). Subjects at other dose levels and those who did not complete the study were excluded from this analysis. Thus, the figure was limited to include data from the last PK profile in each trial to illustrate the variability at steady-state conditions.

Fig. 3a–b

Variability of semaglutide PK under steady-state conditions (a) and model fits to dose-normalised semaglutide concentration–time profiles at steady state for different body weight quantiles (b). For a, points with error bars are observed steady-state PK data (geometric mean and 90% range) for 173 completing subjects at selected time points with PK sampling across trials. Overlaid are 58 simulated replicates of data (n = 10,034) in a dosing interval after the latest treatment dose. The light blue line represents the geometric mean of the observed data and the shaded area represents the 90% prediction interval. The figure includes data from six trials with steady-state PK and subjects treated with 0.5 mg and 1.0 mg semaglutide as well as the model from nine clinical pharmacology trials of semaglutide. For b, data points with error bars are observed geometric means and 95% CIs from seven clinical pharmacology trials with semaglutide at steady state concentrations after weekly dosing of 0.5 or 1.0 mg. Model predictions are population predictions from the final two-compartment model of semaglutide PK. CI confidence interval, PK pharmacokinetic(s)

The dose-normalised time profiles across body weight in 10-kg intervals (Fig. 3b) spanned a similar range to the overall 90% range, illustrating that body weight does explain most of the variability. Furthermore, the close fit of the model to data in each of the 10-kg intervals illustrates that the influence of body weight on the different parameters was adequately described by the model.

The above findings indicate that the two-compartment PK model provided an accurate description of the PK of semaglutide across clinical pharmacology trials, including different populations, administration routes and drug product strengths, and was found suitable for predictive purposes, i.e. for the simulation of semaglutide concentration–time profiles in various populations. Moreover, the low interindividual variation of ~ 15% for both CL and Vc as well as the low residual error (< 5%; Table 4) indicate that the PK of semaglutide is highly predictable.

Estimation of a One-Compartment Model for Steady-State Evaluation of Subcutaneous Data

The population PK model was compared to a one-compartment model similar to the one used in the phase 3 analysis. The two models provided almost identical apparent exposure estimates for individual subjects (Fig. S3a in the ESM). In addition, mean simulated PK profiles during subcutaneous administration were almost identical, with similar exposure levels as well as peak-to-trough fluctuations (Fig. S3b in the ESM). The PK parameter estimates in the one-compartment model used for comparison were similar to those obtained with the two-compartment model (Table S2 in the ESM).

Comparison with the Model Estimated with Phase 3 Data

Parameter values obtained in the present analysis for a reference subject profile were compared to the corresponding parameter values obtained with phase 3 data [6] (Table 5). Due to limitations in the phase 3 data, that analysis applied a fixed value of ka. Although the total volume of distribution was lower in the present analysis than in the phase 3 analysis, profiles from these models were similar (Fig. 4), indicating that the higher volume of distribution in the phase 3 model had a relatively small impact on the overall PK profile. High consistency between the simulated profiles indicate that the PK in the phase 3 trials could be predicted using the two-compartment model that was developed using clinical pharmacology data, supporting the overall predictability of semaglutide PK. Furthermore, the estimated inter-individual variability in clearance (15.2%; Table 4) was consistent with phase 3 data (12.9%) [6].

Table 5

Comparison of pharmacokinetic parameter values from different semaglutide models

Parameter	Two-compartment final model	One-compartment model	One-compartment model from phase 3a trialsa
Absorption rate constant ka, h−1	0.0253	0.0296	0.0286 (fixed)c
CL/F for subject with T2D and body weight of 85 kg, L/h	0.0460	0.0473	0.0478
CL/F for healthy volunteer with body weight of 85 kg, L/h	0.0411	0.0415	–
Total volume/Fb, L	9.08	9.77	12.2

CL/F apparent clearance, F bioavailability, T2D type 2 diabetes, Vc central volume, Vp peripheral volume, ka absorption rate constant

aValues obtained from [6]

bTotal volume/F is (Vc + Vp)/F for the two-compartment model and V/F for the one-compartment model

cDue to sparse sampling in phase 3 trials, ka was fixed at the value obtained from a previous version of the one-compartment model based on clinical pharmacology data

Fig. 4

Simulated PK profiles under steady-state conditions for a reference subject profile based on the final two-compartment model developed from clinical pharmacology data and the one-compartment model developed from phase 3 data

Discussion

The PK of semaglutide were well characterised by a two-compartment model with first-order absorption and elimination, and were shown to be predictable across clinical pharmacology and phase 3 trials and across subjects with low interindividual variation and low residual error. This model provided a joint description of semaglutide PK across clinical pharmacology trials in both healthy subjects and those with type 2 diabetes with frequently sampled PK profiles following subcutaneous and intravenous administration in well-controlled settings.

The total volume of distribution (Vc + Vp) was approximately 8 L, similar to the volume of extracellular water in the body and supporting the idea that semaglutide distributes into plasma and peripheral tissues to the same extent as albumin [7, 8]. The terminal half-life was driven by the disposition parameters rather than the rate of absorption. The population PK model for semaglutide was similar overall to that developed for dulaglutide [9], whereas liraglutide had much faster elimination [10], and exenatide extended-release had nonlinear disposition and a substantially slower absorption [11]. Compared to semaglutide, dulaglutide has approximately 50% greater CL and a 20% lower volume of distribution, as well as approximately threefold slower absorption [9]. The interindividual variability was, however, greater for dulaglutide (40.5% coefficient of variation (CV) for CL) compared to semaglutide (CV for CL was 15.2% in the present analysis and 12.9% in the phase 3 analysis [6]). Similarly, both exenatide extended-release and liraglutide had large variabilities between subjects compared to semaglutide [10, 11].

Based on the present steady-state data, we compared the developed population PK model to a one-compartment PK model similar to the one estimated in the phase 3 analysis [6] with almost identical results. Moreover, we verified that the typical steady-state profiles obtained from phase 3 data [6] were consistent with those predicted based on clinical pharmacology data, and similar estimates for clearance CV% also verified that variability was consistent compared to phase 3 trial results.

In the phase 3 population PK analysis, only body weight had a relevant effect on the exposure of semaglutide, a finding which did not warrant dose adjustment [12]. No effects of sex, age, race, ethnicity, renal function or injection site used on semaglutide exposure were observed [6]. The present analysis confirmed the finding that only body weight was a relevant covariate for exposure. Compared to the phase 3 evaluation, the current analysis provides a more general description of semaglutide PK that is applicable across different administration routes, and informs on covariate relationships not only for exposure, but also for absorption, distribution and elimination. Additional covariate effects were identified for body weight on the Vc and Vp and glycaemic status on CL and the ka. Other investigated covariates (sex, race, ethnicity, age group, hepatic impairment) were not found to be influential and were not included in the final model. The impact of covariates on the PK of semaglutide support previous findings [6, 12] that no adjustment of the dose is needed across population subgroups.

With once-daily liraglutide, another GLP-1 receptor agonist approved for the treatment of type 2 diabetes at doses of up to 1.8 mg, body weight and sex were found to be the most important covariates for exposure which, however, did not influence liraglutide PK to any clinically relevant degree [13]. Similarly, neither age, body weight, sex, race, ethnicity nor injection site influenced once-weekly semaglutide or dulaglutide PK to any clinically relevant extent [6, 9]. For twice-daily exenatide and once-weekly exenatide extended-release, body weight was also found to have a significant relationship with the volume of distribution, and a significant relationship between renal function and elimination was found [11, 14].

The slightly higher CL of semaglutide in subjects with type 2 diabetes compared to healthy subjects has also been seen previously with liraglutide 1.8 mg [10] and with liraglutide for weight management at a higher 3.0 mg dose [15]. The mechanism behind this difference is not known.

The present analysis showed that injection site had a small effect on the bioavailability of semaglutide, which was lower following injection in the thigh compared to abdominal injection. This observation was in accordance with previous findings that the injection site has a minor and not clinically relevant effect on the CL/F of semaglutide [6]. Similarly, the bioavailability of liraglutide was found to be slightly lower after administration in the thigh versus the abdomen [16]. Differences between injection sites are generally considered to be of no clinical relevance across GLP-1 receptor agonists [9, 17]. A limitation of the present analysis is that we cannot exclude the possibility that other demographic factors that we did not evaluate were influential in the target patient population.

Conclusions

In summary, semaglutide PK was shown to be predictable across clinical pharmacology and phase 3 trials, with limited interindividual variability that was primarily influenced by body weight. A novel two-compartment model provided a general characterisation of semaglutide absorption following subcutaneous administration and of distribution and elimination across administration routes. A one-compartment model structure was found to be suitable for conducting covariate analyses of exposure levels based on sparse PK sampling in late-stage trials for subcutaneous semaglutide.

Below is the link to the electronic supplementary material.

Supplementary material 1 (DOCX 313 kb)