Population pharmacokinetics of riociguat and its metabolite in patients with chronic thromboembolic pulmonary hypertension from routine clinical practice

Abstract

Pharmacokinetic data for riociguat in patients with chronic thromboembolic pulmonary hypertension (CTEPH) have previously been reported from randomized clinical trials, which may not fully reflect the population encountered in routine practice. The aim of the current study was to characterize the pharmacokinetic of riociguat and its metabolite M1 in the patients from routine clinical practice. A population pharmacokinetic model was developed in NONMEM 7.3, based on riociguat and its metabolite plasma concentrations from 49 patients with CTEPH. One sample with riociguat and M1 concentrations was available from each patient obtained at different time points after last dose. Age, bodyweight, sex, smoking status, concomitant medications, kidney and liver function markers were tested as potential covariates of pharmacokinetic of riociguat and its metabolite. Riociguat and M1 disposition was best described with one-compartment models. Apparent volume of distribution (Vd/F) for riociguat and M1 were assumed to be the same. Total bilirubin and creatinine clearance were the most predictive covariates for apparent riociguat metabolic clearance to M1 (CL_f,M1/F) and for apparent riociguat clearance through remaining pathways (CL_e,r/F), respectively. CL_f,M1/F, CL_e,r/F, Vd/F of riociguat and M1, and clearance of M1 (CL_e,M1/F) for a typical individual with 70 mL/min creatinine clearance and 0.69 mg/dL total bilirubin were 0.665 L/h (relative standard error = 17%)), 0.66 (18%) L/h, 3.63 (15%) L and 1.47 (19%) L/h, respectively. Upon visual identification of six outlying individuals, an absorption lag-time of 2.95 (6%) h was estimated for these patients. In conclusion, the only clinical characteristics related to riociguat exposure in patients with CTEPH from routine clinical practice are total bilirubin and creatinine clearance. This confirms the findings of the previous population pharmacokinetic studies based on data from randomized clinical trials.

Keywords

desmethylriociguat NONMEM creatinine clearance total bilirubin

Chronic thromboembolic pulmonary hypertension (CTEPH) is a pulmonary vascular disease caused by the chronic thrombotic obstruction of pulmonary arteries and peripheral vascular remodeling.¹ The disease is characterized by elevation of pulmonary artery mean pressure (mPAP) ≥ 25 mmHg and increase of pulmonary vascular resistance (PVR), and ultimately results in death due to the right ventricular failure.²

Most CTEPH patients are successfully treated with surgical pulmonary endarterectomy (PEA), which is the gold standard in management of CTEPH. However, up to 40% of patients with CTEPH are considered technically inoperable due to distal lesions that are not surgically accessible or due to comorbidities, and 20–30% of patients develop persistent/recurrent pulmonary hypertension after PEA.³ Both patient groups are potential candidates for the treatment with riociguat – a novel medication that relaxes vascular smooth muscle by stimulation of soluble guanylate cyclase (sGC). Riociguat displays a dual mode of action: it stimulates sGC independently of nitric oxide (NO) and increases the sensitivity of sGC to NO, resulting in increased cyclic guanosine monophosphate (cGMP) levels.^4,5

The main biotransformation pathway for riociguat is N-demethylation by cytochrome P450 enzymes, most importantly CYP1A1.⁶ Desmethylriociguat, the major circulating active metabolite M1, which exhibits 1/10th to 1/3rd of the pharmacological activity of riociguat, is further metabolized to the pharmacologically inactive N-glucuronide.⁷ The drug is eliminated in the urine (33–45%) and feces (48–59%).⁷

Pharmacokinetic (PK) data for riociguat have previously been reported from the randomized clinical trials (RCT) CHEST-1, 2 and PATENT-1, 2 in CTEPH and pulmonary arterial hypertension (PAH) patients, respectively.^8–10 RCTs are conducted in tightly controlled settings and include patients who meet stringent inclusion and exclusion criteria and may therefore not accurately reflect the population treated in clinical practice.^11,12 Therefore, the aim of the current study was to characterize the PK of riociguat and its metabolite desmethylriociguat (M1) in patients with inoperable CTEPH or persistent/recurrent pulmonary hypertension after PEA from routine clinical practice and to identify and quantify significant covariates associated with riociguat and M1 exposure in these patients.

Methods

Research design

This observational PK study was conducted at the Second Department of Internal Medicine, General University Hospital and First Faculty of Medicine, Charles University in Prague, Czech Republic. It was conducted in accordance with the principles laid down in the 18th World Medical Assembly (Helsinki, 1964), including all subsequent amendments, and in compliance with all laws and regulations of the Czech Republic. The approval of retrospective data collection was provided by ethics committee of the General University Hospital in Prague (ID 1208/18 S-IV). Written informed consent was obtained from all participants.

Patients were included in the study if they were diagnosed with inoperable CTEPH or persistent/recurrent pulmonary hypertension after PEA and received a stable riociguat dose for at least three months before the enrolment. Inoperability status was previously assessed by an interdisciplinary CTEPH team, consisting of a pulmonary hypertension specialist, a PEA surgeon, an anesthesiologist, and a radiologist. Persistent/recurrent pulmonary hypertension was diagnosed invasively by the right heart catheterization at least six months after PEA and defined as persistent elevation of mPAP ≥ 25 mmHg and PVR > 3 Wood unit.

The following data were collected from the outpatient check-up visit: time of the last dose and sampling, demography, medical history, concomitant medications, vital signs, functional capacity, and 6-minute walking distance (6MWD). Blood samples were collected for determination of N-terminal pro b-type natriuretic peptide (NT-proBNP), laboratory biochemical parameters, and riociguat and M1 concentrations. Retrospective data from the last right heart catheterization performed before initiation of riociguat treatment were used for the description of hemodynamics.

Riociguat dosing

Riociguat was prescribed according to Adempas® Summary of Product Characteristics (SmPC) and 2015 European Respiratory Society/European Society of Cardiology (ERS/ESC) treatment guidelines including required initial dose adjustments.¹ Riociguat was adjusted from a starting dose of 1 mg three times daily according to systolic systemic arterial pressure and signs or symptoms of hypotension (final range: from 1.5 mg to 2.5 mg three times daily).

Bioanalytical assay

After collection, blood samples were allowed to clot for 30 min at room temperature, and serum was separated by centrifugation (1500 × g, 15 min, 4℃) and stored frozen at –80℃ until analysis. Riociguat and M1 serum concentrations were measured using liquid chromatography with tandem mass spectrometric detection in positive ESI mode (LC-MS/MS); penta-deuterated perampanel was used as internal standard (IS). Internal standard solution (10 µL, 5000 ng mL⁻¹ in methanol) and 300 µL acetonitrile were added to 100 µL of the serum sample in a 1.5 mL Eppendorf tube. The solution was vortexed for 30 s and centrifuged (9600 × g, 3 min) and the supernatant (100 µL) was transferred to an autosampler vial. The method was developed using Nexera X2 Shimadzu HPLC (Nakagyo-ku, Kyoto, Japan) coupled with AB Sciex QTRAP 5500 (MA, USA). Mobile phase A consisted of 0.1% formic acid in water and mobile phase B consisted of acetonitrile. The analysis was performed on a Zorbax Eclipse XDB-C18 column (1.8 µm, 50 × 4.6 mm). The initial LC conditions had a flow rate of 0.5 mL min⁻¹ at a mobile phase composition of 50:50 (A: B). At 30 s the mobile phase composition was ramped to 10:90 (A: B) within 60 s and held for 90 s and then returned to initial LC conditions. Quantitation was done using multiple reaction monitoring (MRM) mode to monitor protonated precursor → product ion transition of m/z 423.022 → 109.100 for riociguat, 409.027 → 109.000 for M1 and 355.029 → 220.000 for penta-deuterated perampanel. Method performance was evaluated for riociguat and M1 following the recommendations of the Scientific Working Group for Forensic Toxicology.¹³ The test range of the assay was 5–1000 µg/L. Coefficient of variation of intra-assay was less than 11%.

Population PK analysis

The data analysis was performed using NONMEM version 7.3.0 (ICON Development Solutions, Ellicott City, MD) and PsN v3.4.2 both running under Pirana 2.9.0. The first-order estimation algorithm with interaction (FOCE-I) was used. R 3.3.2 was used for the visualization of the data and model diagnostics.

Model development was performed in three steps:

(1) Development of the structural and statistical model. For the structural model, one and two compartment models were tested to describe the distribution of riociguat and M1. Assumptions for the structural model were necessary to ascertain mathematical identifiability of the parameter values.^14,15 The same values of volume of distribution (Vd) of riociguat and its metabolite were assumed.¹⁶ For the metabolic formation clearance of M1 (CL_f,M1), the elimination clearance of M1 (CL_e,M1) and the remaining riociguat elimination clearance trough alternative pathways (CL_e,r) (Figure 2) standard first-order processes were assumed. Since sufficiently dense data in the absorption phase were lacking, the absorption rate constant (k_a) was fixed to a value of 3 h⁻¹ obtained from the literature.¹⁷ Since only one sample per patient was available, it was not possible to separately estimate intra-individual variability (IIV) and residual unexplained variability (RUV). Therefore, proportional, additive, and combination error models were tested to represent both IIV and RUV. In the absence of observations upon intravenous administration all parameters represent apparent values.

Based on visual inspection of raw data, six patients had unexpectedly high riociguat concentrations at relatively late times after their last dose (Figure 1). These patients’ data were initially excluded from model development. After the covariate analysis, these patients were reintroduced and analyzed together with the other individuals.

Figure 1.

Serum concentrations plotted against time after the last dose. ID numbers from observations from individuals with outlying riociguat concentrations are indicated. (a) Riociguat; (b) Desmethylriociguat.

Figure 2.

Schematic representation of the pharmacokinetic model for riociguat and its metabolite in patients with chronic thromboembolic pulmonary hypertension. Pharmacokinetic parameters represent apparent values. Abbreviations: k_a/F = apparent absorption rate constant, V/F = apparent distribution volume of the designated compartment (p = parent – riociguat, m = metabolite), CL_f,M1/F = apparent riociguat metabolic formation clearance to M1, CL_e,r/F = apparent riociguat clearance remaining after accounting for M1 formation, CL_e,M1/F = apparent clearance of M1, F = bioavailability.

For model selection, a decrease in objective function of more than 6.63 points between nested models (p < 0.01) was considered statistically significant, assuming a χ²-distribution. Additional criteria for model selection were relative standard error (RSE) of the estimates of structural model parameters <50%, condition number calculated by dividing the largest and smallest eigenvalue from the model fit of <1000, physiological plausibility of the obtained parameter values, and absence of bias in goodness-of-fit (GOF) plots.

(2) Covariate analysis. In the systematic covariate analysis, the following continuous covariates were tested in linear and exponential equations: age, body weight (BW), ideal body weight (IBW), body mass index (BMI), creatinine clearance (calculated by CKD-EPI (Chronic Kidney Disease Epidemiology Collaboration) equation), total bilirubin, aspartate transaminase (AST) and alanine transaminase (ALT), alkaline phosphatase (ALP), gamma-glutamyl-transferase (GGT), and NT-proBNP levels in plasma. The following categorical covariates were tested by estimating the parameter value for one category as a fraction of the parameter value for the other category: concomitant therapy (inhibitors of proton pump, diuretics, digoxin, and angiotensin converting enzyme (ACE) inhibitors), sex, smoking status. All continuous and categorical covariates were tested on the following parameters: CL_f,M1/F, CL_e,r/F, Vd/F of riociguat and M1, and CL_e,M1/F. The criteria used for model selection were the same as those described above.

After the final covariate model was developed, the patients with outlying observations were reintroduced into the analysis. As previous reports suggested that food delays riociguat absorption for 3 h,^7,18 it was investigated whether riociguat absorption in these patients was delayed. For this, both different k_a values and a delayed onset of absorption as characterized by a lag-time (Tlag) were tested based on the criteria defined above. When it was observed that this procedure did not yield large differences in the obtained parameter values, the results from the fit which included all individuals were retained.

(3) Validation of the final model. To evaluate the robustness of the final model and the precision of the parameter estimates, a bootstrap analysis was performed on the final model. Two hundred and fifty bootstrap datasets were generated by random sampling with replacement. The final model was rerun with the new datasets and median parameter values, 2.5th and 97.5th percentiles of parameter distribution and standard error of the estimates were generated and compared to the parameters of the final model.

The predictive properties of the structural and statistical model were evaluated using normalized prediction distribution errors (NPDEs), a simulation-based diagnostics. For this, the dataset was simulated 500 times, after which the observed concentrations were compared to the range of simulated values using the NPDE package developed for R.¹⁹

Results

Study population and data

In total, 49 (24 female, 25 male) patients ((median (interquartile range) age: 74 (66–78) years, BW: 74 (66–78) kg) with CTEPH receiving long-term riociguat treatment were included in our analysis. Characteristics of the patient population are summarized in Table 1. One sample with riociguat and M1 concentrations was available from each patient obtained at different time points after last dose, ranging from 1.25 to 6.75 h. Riociguat and its metabolite levels ranged between 44 and 749 µg/L, and 17 and 314 µg/L, respectively. Figure 1 shows the riociguat and M1 concentrations plotted against time after the last dose. Six outlying riociguat concentrations for ID = 6, 10, 11, 19, 25, 28 were visually identified.

Table 1.

Clinical characteristics of the study population.

Parameter (unit)	Value^a
Age (years)	74 (66–78)
Sex, female/male, n (%)	24/25 (49/51%)
Weight (kg)	80 (67–95)
Ideal bodyweight (kg)	61 (55–74)
Body mass index (kg/m²)	27.8 (23.8–30.8)
Riociguat dose (mg/day)	7.5 (6.75–7.5)
Duration of the treatment (months)	21 (15–27)
Diagnosis
Inoperable CTEPH, n (%)	37 (74%)
Residual CTEPH, n (%)	13 (26%)
Hemodynamics
Mean PAP (mm Hg), mean ± SD	43 ± 12
RAP (mm Hg), mean ± SD	7.6 ± 3.9
Cardiac output (L/min), mean ± SD	4.5 ± 1.0
PVR (WU), mean ± SD	7.8 ± 3.2
Systolic BP (mm Hg), mean ± SD	138 ± 20
Diastolic BP (mm Hg), mean ± SD	77 ± 12
Heart Rate (BPM), mean ± SD	75 ± 11
Exercise capacity
6MWD (m), mean ± SD	387 ± 119
Functional class
NYHA I/II/III/IV, n (%)	0/21/28/0 (0/42.8/57.2/0%)
Laboratory markers
Creatinine (µmol/L)	87 (70.5–102.5)
Total bilirubin (mg/dL)	0.69 (0.53–0.98)
Creatinine clearance (mL/min)	70 (59–79)
NT-proBNP (ng/L)	493 (235–1460)
Aspartate transferase (IU/L)	0.40 (0.31–0.46)
Alanine transferase (IU/L)	0.28 (0.21–0.39)
Alkaline phosphatase (IU/L)	1.23 (0.96–1.65)
Gamma-glutamyl-transferase (µIU/L)	0.52 (0.32–0.92)
Concomitant medication
Diuretics, n (%)	35 (71%)
Inhibitors of proton pump, n (%)	35 (71%)
Digoxin, n (%)	3 (6.1%)
Warfarin, n (%)	41 (84%)
ACE inhibitors, n (%)	11 (22.4%)
ARBs, n (%)	4 (8.2%)
H₂ antihistamines, n (%)	2 (4%)
Smoking status
Smoker, n (%)	1 (2%)

Values are presented as median (inter-quartile range), unless noted otherwise.

CTEPH: chronic thromboembolic pulmonary hypertension; PAP: pulmonary artery pressure; RAP: right atrial pressure; PVR: pulmonary vascular resistance; BP: blood pressure; 6MWD: 6-minute walking distance; NYHA: New York Heart Association; NT-proBNP: N-terminal pro b-type natriuretic peptide; ACE: angiotensin converting enzyme; ARBs: angiotensin II receptor blockers.

Population PK analysis

Observed riociguat and M1 serum concentrations were best described with one-compartment models and the same distribution volume was estimated for both compounds to achieve mathematical identifiability. Proportional residual error models provided the best description of the residual variability for both riociguat and M1 concentrations.

Figure 2 depicts a schematic representation of the obtained PK model. Creatinine clearance in a linear equation was found to be the most predictive covariate for CL_e,r/F, was found to increase 0.009 L/h per unit (mL/min) creatinine clearance, as depicted in Figure 3(a). For CL_f,M1/F, the most predictive covariate relationship was total bilirubin in an exponential equation with an estimated exponent of −0.463 (20%), as shown in Figure 3(b). After the inclusion of these covariate relationships, no other statistically significant covariates could be identified. CL_f,M1/F, CL_e,r/F, Vd/F of riociguat and M1, and CL_e,M1/F for a typical individual of creatinine clearance (70 mL/min) and total bilirubin level (0.69 mg/dL) were 0.665 L/h (17%)), 0.66 (18%) L/h, 3.63 L (15%) and 1.47 (19%) L/h, respectively. Additionally, the analysis showed that the riociguat absorption in the six patients exhibiting high concentrations was delayed with a lag-time of 2.95 h (6%). The parameter values obtained in the final model fit as well as the median parameter values obtained in the bootstrap procedure are presented in Table 2.

Figure 3.

Relationships between (a) total bilirubin and apparent riociguat metabolic formation clearance to M1 and (b) creatinine clearance and apparent riociguat clearance remaining after accounting for M1 formation in patients with CTEPH.

Table 2.

Parameter estimates of the final model and their corresponding bootstrap estimates.

Parameter (unit)	Final model (% RSE)	Bootstrap (95% CI)
Fixed effects
K_a/F (h⁻¹)	3 FIX
CL_e,r/F (L/h) = CL_e,rTV * (CREACL/70)
CL_e,rTV (L/h)	0.66 (17%)	0.655 (0.198–0.959)
CL_f,M1/F (L/h) = CL_f,M1TV * (BILTOT/0.69)^θBILTOT
CL_f,M1TV (L/h)	0.665 (17%)	0.671 (0.446–1.130)
θBILTOT	−0.462 (20%)	−0.448 (−0.625 to −0.143)
V_P/F = V_M/F (L)	3.63 (15%)	3.78 (2.81–5.60)
CL_e,M1/F (L/h) = CL_e,M1TV	1.47 (19%)	1.44 (0.90–2.43)
Tag (ID = 6,10,11,19,25,28) (h)	2.95 (6%)	2.98 (2.68–3.81)
Variance of residual variability
Riociguat, proportional	0.152 (21%)	0.141 (0.086–0.217)
M1, proportional	0.268 (18%)	0.254 (0.165–0.362)

Pharmacokinetic parameters represent apparent values.

k_a/F: apparent absorption rate constant; CL: apparent clearance of the designated pathway (see Figure 2 for the explanation of the symbols); CREACL: creatinine clearance in mL/min; BILTOT: total bilirubin level in mg/dL; θBILTOT: exponent for the covariate relationships between total bilirubin levels and CL_f,M1/F; V_P/F: apparent volume of distribution of the parent compound; V_M/F: apparent volume of distribution of the M1 metabolite; TV: typical value of a parameter; Tlag: lag-time.

RSE values for the structural parameters were all below 30% in the final model, indicating good precision of the estimated parameters. All median parameter values in the bootstrap procedure were within 10% of the values obtained in the final model fit indicating that the model is robust. Figures 4 and 5 present the GOF plots for riociguat and M1. Absence of bias in these plots indicates that the final model describes the observed data accurately. Finally, there was no bias in NPDE, neither over time nor over concentration range, indicating that the predictive properties of this model are also accurate (Supplementary Figures 1S and 2S).

Figure 4.

Goodness of fits (GOF) plots for the final model for riociguat pharmacokinetics in patients with chronic thromboembolic pulmonary hypertension. (a) Population predicted concentration vs. observed concentration. (b) Individual predicted concentration vs. observed concentration. (c) conditional weighted residuals (CWRES) vs. time after the last dose. (d) Conditional weighted residuals (CWRES) vs. population predicted concentration. ID numbers from observations from individuals with outlying riociguat concentrations are indicated.

Figure 5.

Goodness of fits (GOF) plots for the final model for M1 pharmacokinetics in patients with chronic thromboembolic pulmonary hypertension. (a) Population predicted concentration vs. observed concentration. (b) Individual predicted concentration vs. observed concentration. (c) conditional weighted residuals (CWRES) vs. time after the last dose. (d) Conditional weighted residuals (CWRES) vs. population predicted concentration. ID numbers from observations from individuals with outlying riociguat concentrations are indicated.

Discussion

This study used a population modelling approach to describe the PK of riociguat and its metabolite M1 in patients with CTEPH from routine clinical practice. The results of the study showed that riociguat metabolic clearance to M1 depends on the liver function, as characterized by the total bilirubin level. Creatinine clearance, a marker for kidney function, was found to be a predictive covariate for riociguat clearance remaining after accounting for M1 formation. These results indicate that impaired renal and hepatic function leads to reduced riociguat clearance and increased riociguat exposure.

Previous PK studies analyzed data obtained during RCTs, which are constrained only to patients that meet strict inclusion and exclusion criteria and may therefore not accurately reflect the population treated in the routine practice.^11,12 The current analysis included data collected in everyday practice and may therefore provide more relevant clinical information. However, this type of data are usually not as “clean” as those from RCTs. In this study for example, six patients had unexpectedly high riociguat concentrations at the time of the observation that could not be explained (Figure 1(a)). Because of the inclusion of one sample per patient in the analysis, it was not possible to distinguish whether these patients or the particular observations were outliers. As previous reports however suggested that food might delay riociguat absorption for 3 h,¹⁸ we investigated whether the observations could be explained by delayed riociguat absorption in these patients as characterized by a Tlag, or slowed down, as characterized by a different k_a/F value. The analysis showed that data from these patients were best described with a delayed riociguat absorption with a Tlag of 2.95 h, which might be due to food impact, but the design of our study does not allow for conclusions to be drawn on this matter. An important clinical question is whether this delay in absorption would result in different riociguat exposure. If these patients indeed had a delayed riociguat absorption due to food, then changes in area under curve (AUC) would not be expected. However, if there is another reason for these concentrations to be so high, AUC could in fact be impacted. To make a final conclusion, we would need to have multiple observations per patient.

Although only one sample per patient was included in the analysis, samples were obtained within a wide range of time points after the last dose (Figure 1), which allowed us to develop a structural model and identify the predictive covariates for riociguat and M1 PK. PK parameter values could be obtained with acceptable precision, as reflected in low RSE (<30%) of the parameter estimates. In addition, extensive model validation showed that the model not only described the obtained data well (Figures 4 and 5), but also predicted the data well (Supplementary Figures 1S and 2S), meaning that the conclusions regarding parameter values and covariate effects in this model are well supported by the data.

An important advantage of the current analysis is the wide range of sampling time points after the last dose of riociguat in CTEPH patients, which covers a wide range than the previous multiple-dose population PK study, which included only trough samples from the PAH and CTEPH patients and additionally samples obtained 2–3 h after the first and second dose of drug only.⁹

Direct comparison of findings between studies is difficult due to differences in parameterization and covariate relationships. Still, it is possible to make comparisons between parameter values for typical individuals. There are two studies using a population modelling approach to describe PK of riociguat and its metabolite.^9,17 One study described the PK of a single dose riociguat,¹⁷ whereas the other addressed PK of riociguat upon multiple dosing at steady state.⁹ To allow comparison of PK parameters, we calculated parameter values for the typical individual from our study with bodyweight of 80 kg, creatinine clearance of 70 mL/min and total bilirubin of 0.69 mg/dL, using the provided equations in the respective publications (Table 3). Interestingly, estimated values for Vd/F of riociguat and M1 dramatically varied between all three studies, probably due to the heterogeneity in clinical features of the patients in the studies, such as drug-protein binding, co-medication and hemodynamic characteristics, or due to differences in assumptions regarding the distribution volume. On the other hand, values for the elimination clearance of riociguat (CL_e,r/F) and M1 (CL_e,M1/F) are similar to the previously reported values, but the CL_f,M1/F value estimated in the multiple dosing study deviated significantly from the values obtained in our analysis and in the single dose study, possibly due to differences in assumptions. We found total bilirubin levels to be a significant covariate for CL_f,M1/F, contrary to previous studies, which reported high IIV of CL_f,M1/F, independently of renal or hepatic status. Creatinine clearance was found to be a predictive covariate of CL_e,r/F, which is in accordance with the previous multiple dose study. These results indicate that impaired renal and hepatic function result in reduced riociguat clearance and increased riociguat exposure. As similar findings were reported in the previous studies,^9,17 our study confirms that kidney and liver functions are the main clinical characteristics related to riociguat exposure in patients with CTEPH. Therefore, particular care should be taken in patients with renal and hepatic impairment.

Table 3.

Overview of population pharmacokinetic studies evaluating pharmacokinetic parameters of riociguat and its metabolite.

Population	Distribution Model	V_P/F (L)^a	V_M/F (L)^a	CL_e,r/F (L/h)^a	CL_f,M1/F (L/h)^a	CL_e,M1/F (L/h)^a	SC for V_P/F	SC for V_M/F	SC for CL_e,r/F	SC for CL_f,M1/F	SC for CL_e,M1/F	Reference
49 CTEPH patients BW: 80 (2.8–3.5) kg BILTOT: 0.69 (0.53–0.98) mg/dL CREACL: 70 (59–79.5) mL/min	1-CPT	3.63	3.63	0.66	0.665	1.47	–	–	CREACL	BILTOT	–	Current analysis
Single dosing study 64 patients with renal impairment BW: 83.3 (56–108.3) kg BILTOT: 0.5 (0.2–7.2) mg/dL CREACL: 75.1 (17.4–125.5) mL/min 72 patients with hepatic impairment BW: 86.5 (54.2–117) kg BILTOT: 0.4 (0.1–6.1) mg/dL CREACL:) 72.7 (7.0–129.3) mL/min	2-CPT	17.2 + 10.6^b	8.6 + 34.2^b	0.712	1.23	1.76	BW, %PB	BW	%PB	–	CREACL	¹⁷
Multiple dosing study 260 CTEPH patients BW: 74 (36–158.3) kg BILTOT: 0.6 (0.1–5.0) mg/dL CREACL: 74.4 (15.4–233) mL/min 438 PAH patients BW: 65.3 (37.7–141 kg BILTOT: 0.5 (0.1–4.6) mg/dL CREACL:) 86.9 (15.9−264) mL/min	1-CPT	34.7	133	1.76	8.424^c	3.07	BW	BW	CREACL, BILTOT, SMOK, BOS	–	CREACL, BILTOT, SMOK, BOS	⁹

CPT: compartment; CL/F: clearance of the designated pathway (see Figure 2 for the explanation of the symbols); SC: significant covariate; CREACL: creatinine clearance; BILTOT: total bilirubin level; V_P/F: parent volume of distribution; V_M/F: metabolite volume of distribution; CTEPH: chronic thromboembolic pulmonary hypertension; PAH: pulmonary arterial hypertension; BW: body weight; %PB: protein binding; SMOK: smoking; BOS: Bosentan co-medication.

Values for Vd/F and CL/F are calculated for the typical individual in this study (CREACL = 70 mL/min, BILTOT = 0.69 mg/dL, BW = 80 kg).

Value for the central compartment + value for the peripheral compartment.

Calculated from the relationship k₂₃ = CL₂₃/V_P.

It is important to note that not all PK parameters for M1 could be estimated without making assumptions. As data obtained after intravenous administration of the M1 metabolite, or data on the recovery of M1 in urine were not available, the value of Vd/F for M1 was assumed to be the same as the parent compound. The model validation confirmed that the model can accurately describe and predict the concentrations of riociguat and M1, but as a result of the assumption the absolute values of the parameters related to the metabolite should be considered in the context of the assumptions made in the current analysis. The conclusions regarding the impact of the covariates are not impacted by the assumptions.

In conclusion, we report on the PK of riociguat and its pharmacologically active metabolite desmethylriociguat in a cohort of CTEPH patients encountered in routine clinical practice. Our study confirms the findings from previous population PK studies based on data from RCTs, that the only clinical characteristics related to riociguat exposure in patients with CTEPH are total bilirubin levels and creatinine clearance.

Supplemental Material

PUL898031 Supplemental material - Supplemental material for Population pharmacokinetics of riociguat and its metabolite in patients with chronic thromboembolic pulmonary hypertension from routine clinical practice

Supplemental material, PUL898031 Supplemental material for Population pharmacokinetics of riociguat and its metabolite in patients with chronic thromboembolic pulmonary hypertension from routine clinical practice by Danica Michaličková, Pavel Jansa, Miroslava Bursová, Tomáš Hložek, Radomír Čabala, Jan Miroslav Hartinger, David Ambrož, Michael Aschermann, Jaroslav Lindner, Aleš Linhart, Ondřej Slanař and Elke H.J. Krekels in Pulmonary Circulation

Footnotes

Acknowledgments

The authors would like to thank all the patients, nurses and physicians who were part of this study. The authors also thank Dr. Parth Upadhyay for code review.

Funding

This study was supported by the Charles University projects Progres Q25 and Q38 and by the project “International Mobility of Researchers at Charles University” CZ.02.2.69/0.0/0.0/16_027/0008495.

Availability of data and materials

The data are not available in any public repository. The model codes will be made available through the model repository of DDMoRe available through: .

Ethical approval

The approval of retrospective data collection was provided by ethics committee of the General University Hospital in Prague (ID 1208/18 S-IV). Written informed consent was obtained from all participants.

Guarantor

Not applicable.

Contributorship

D.M. analyzed the data and wrote the manuscript; P.J. conceived and designed the study, performed the clinical trial, and wrote the manuscript; M.A., J.L. and A.L. performed the clinical trial; O.S. wrote the manuscript; M.B., T.H., R.Č., and D.A. developed the analytical method and performed laboratory analyses; J.M.H. wrote the manuscript; E.H.J.K. supervised the data analysis and wrote the manuscript.

Conflict of interest

P.J. has received fees and grants from Actelion Pharmaceuticals Ltd, AOP Orphan, and MSD.

Supplemental material

Supplemental material for this article is available online.

References

Galiè

Humbert

Vachiery

, et al. ESC/ERS Guidelines for the diagnosis and treatment of pulmonary hypertension: The Joint Task Force for the Diagnosis and Treatment of Pulmonary Hypertension of the European Society of Cardiology (ESC) and the European Respiratory Society (ERS): Endorsed by: Association for European Paediatric and Congenital Cardiology (AEPC), International Society for Heart and Lung Transplantation (ISHLT). Eur Heart J 2016; 37: 67–119.

Hoeper

Bogaard

Condliffe

, et al. Definitions and diagnosis of pulmonary hypertension. J Am Coll Cardiol 2013; 62: D42–D50.

Pepke-Zaba

Delcroix

Lang

, et al. Chronic thromboembolic pulmonary hypertension (CTEPH) results from an international prospective registry. Circulation 2011; 124: 1973–1981.

McLaughlin

Jansa

Nielsen-Kudsk

, et al. Riociguat in patients with chronic thromboembolic pulmonary hypertension: results from an early access study. BMC Pulm Med 2017; 17: 216.

Schermuly

Janssen

Weissmann

, et al. Riociguat for the treatment of pulmonary hypertension. Expert Opin Inv Drug 2011; 20: 567–576.

Khaybullina

Patel

Zerilli

. Riociguat (adempas): a novel agent for the treatment of pulmonary arterial hypertension and chronic thromboembolic pulmonary hypertension. Pharm Ther 2014; 39: 749.

Frey

Becker

Saleh

, et al. Clinical pharmacokinetic and pharmacodynamic profile of riociguat. Clin Pharmacokinet 2018; 57: 647–661.

Ghofrani

H-A

D'armini

Grimminger

, et al. Riociguat for the treatment of chronic thromboembolic pulmonary hypertension. N Engl J Med 2013; 369: 319–329.

Saleh

Becker

Frey

, et al. Population pharmacokinetics and the pharmacokinetic/pharmacodynamic relationship of riociguat in patients with pulmonary arterial hypertension or chronic thromboembolic pulmonary hypertension. Pulm Circ 2016; 6: S86–S96.

10.

Ghofrani

H-A

Galiè

Grimminger

, et al. Riociguat for the treatment of pulmonary arterial hypertension. N Engl J Med 2013; 369: 330–340.

11.

Grapow

von Wattenwyl

Guller

, et al. Randomized controlled trials do not reflect reality: real-world analyses are critical for treatment guidelines!. J Thorac Cardiovasc Surg 2006; 132: 5–7.

12.

Wang

D-d

Chen

, et al. Model extrapolation to a real-world dataset: evaluation of tacrolimus population pharmacokinetics and drug interaction in pediatric liver transplantation patients. Xenobiotica 2019; 3: 1–9.

13.

Toxicology SWGfF. Scientific Working Group for Forensic Toxicology (SWGTOX) standard practices for method validation in forensic toxicology. J Anal Toxicol 2013; 37: 452–474.

14.

Krekels

van Ham

Allegaert

, et al. Developmental changes rather than repeated administration drive paracetamol glucuronidation in neonates and infants. Eur J Clin Pharmacol 2015; 71: 1075–1082.

15.

Shivva

Korell

Tucker

, et al. An approach for identifiability of population pharmacokinetic–pharmacodynamic models. CPT Pharmacometrics Syst Pharmacol 2013; 2: 1–9.

16.

Bertrand

Laffont

Mentré

, et al. Development of a complex parent-metabolite joint population pharmacokinetic model. AAPS J 2011; 13: 390–404.

17.

Saleh

Becker

Frey

, et al. Population pharmacokinetics of single-dose riociguat in patients with renal or hepatic impairment. Pulm Circ 2016; 6: S75–S85.

18.

Becker

Frey

Hesse

, et al. Absorption of riociguat (BAY 63-2521): bioavailability, food effects, and dose proportionality. Pulm Circ 2016; 6: S27–S34.

19.

Comets

Brendel

Mentré

. Computing normalised prediction distribution errors to evaluate nonlinear mixed-effect models: the npde add-on package for R. Comput Methods Programs Biomed 2008; 90: 154–166.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.14 MB

0.00 MB