Physiologically-Based Pharmacokinetic Modelling to Investigate Baricitinib and Tofacitinib Dosing Recommendations for COVID-19 in Geriatrics.

Janus kinase (JAK) inhibitors baricitinib and tofacitinib are recommended by the US National Institutes of Health as immunomodulatory drugs for coronavirus disease 2019 (COVID-19) treatment. In addition, baricitinib has recently received Emergency Use Authorization from the US Food and Drug Administration, although the instruction provided dosing information only for adults. Geriatrics with organ dysfunction are one of the most vulnerable cohorts when combating the pandemic. The aim of the present work was to evaluate current dosing strategies of baricitinib and tofacitinib for potential COVID-19 treatment for White and Chinese geriatric patients with chronic renal impairment. An established physiologically-based pharmacokinetic (PBPK) modeling framework for age-dependent simulations was utilized. PBPK drug models adopted from literature were first verified. Several population models representing different age groups, ethnicities, and stages of renal impairment were used for prospective simulations. Notwithstanding the increase in systemic exposure of both drugs resulting from renal dysfunction was more pronounced for geriatrics than general White populations, our simulations confirmed their current dosage adjustments based on renal functions are broadly adequate. The exception being White older subjects with mild renal impairment where current recommendation of 4 mg baricitinib yielded a 2.31-fold increase in systemic exposure, and reduction to 2 mg could mitigate the potential risk to an acceptable 1.15-fold. Comparable relationships between systemic exposure and renal dysfunction were observed for both drugs in the Chinese population. In summary, PBPK modeling of both JAK inhibitors supports the rational and prudent dose adjustments of these COVID-19 therapeutics among adult patients of different age groups and renal functions.

WHAT IS THE CURRENT KNOWLEDGE ON THE TOPIC?

☑ Janus kinase inhibitors baricitinib and tofacitinib are recommended by the National Institutes of Health (NIH) for coronavirus disease 2019 (COVID‐19) treatment. Current Emergency Use Authorization instruction issued by the US Food and Drug Administration (FDA) for baricitinib provides information only for adults, whereas geriatric patients with organ dysfunction are under‐represented when combating the pandemic.

WHAT QUESTION DID THIS STUDY ADDRESS?

☑ The research question is whether the current dosing regimens of baricitinib and tofacitinib have to be adjusted for potential COVID‐19 treatment among geriatric patients, especially those with renal impairment?

WHAT DOES THIS STUDY ADD TO OUR KNOWLEDGE?

☑ Current dosage adjustments based on renal functions are broadly adequate for geriatric patients. The exception being older White subjects with mild renal impairment where current recommendation of 4 mg baricitinib yielded a 2.31‐fold increase in systemic exposure, and reduction to 2 mg could mitigate the potential risk.

HOW MIGHT THIS CHANGE CLINICAL PHARMACOLOGY OR TRANSLATIONAL SCIENCE?

☑ Physiologically‐based pharmacokinetic modeling provided prospective evidence for the rational dosing of baricitinib in a vulnerable cohort and supports dose adjustments of baricitinib and tofacitinib among adult patients with COVID‐19 of different age groups and renal functions.

The pathogenesis of coronavirus disease 2019 (COVID‐19) could be driven by a dysregulated immune/inflammatory response to severe acute respiratory syndrome‐coronavirus 2, which often appears at the later clinical course culminating in tissue damage. , Consequently, anti‐inflammatory/immunosuppressive therapies have been included in the US National Institutes of Health (NIH) COVID‐19 Treatment Guideline for hospitalized adults requiring oxygen. The Janus kinase‐Signal Transducer and Activator of Transcription (JAK‐STAT) signaling pathway has been shown to be related to the proinflammatory cytokines and chemokines in the pathology of COVID‐19. Treatment with corticosteroid dexamethasone has conferred a survival benefit of patients with COVID‐19, and for those on dexamethasone who also have rapidly increasing oxygen needs and systemic inflammation, JAK inhibitor baricitinib is recommended by the NIH as an immunomodulatory drug. The US Food and Drug Administration (FDA) has also issued an Emergency Use Authorization (EUA) in November 2020 to permit the emergency use of baricitinib for the treatment of COVID‐19 in hospitalized adults and pediatric patients 2 years of age or older. Another JAK inhibitor, tofacitinib, has been listed in the NIH guideline as an alternative if baricitinib is unavailable or not feasible to be prescribed. Both of them are currently included in the Treatment Guidelines for COVID‐19 issued by Singapore National Centre for Infectious Diseases. Approximately 70% and 30% of baricitinib and tofacitinib undergo renal clearance in humans, respectively.

The risk of severe COVID‐19 if an individual becomes infected is known to be higher in older individuals and the share of the population at increased risk was highest in countries with older populations. Prescribing information of both baricitinib and tofacitinib suggested caution should be used for rheumatoid arthritis treatment among older adults, due to decreased renal function and increased incidence of adverse effects. , Moreover, current EUA for baricitinib is provided only for hospitalized adults and pediatric patients. As advocated by a recently published white paper, a comprehensive plan is needed to ensure adequate evaluation of the safety and effectiveness of drugs in the geriatric population as aging can modify the pharmacokinetics (PKs), pharmacodynamics, and likelihood of adverse effects of a drug. Hence, a question arose whether the current dosing regimens in adults for baricitinib and tofacitinib have to be adjusted for potential COVID‐19 treatment for the White and Chinese older patients, especially for those with chronic and notable renal impairment.

The physiologically‐based pharmacokinetic (PBPK) model is a mathematical framework by which the absorption, distribution, metabolism, and excretion (ADME of a drug can be simulated. It has been recommended by the health authorities to investigate clinical scenarios, which are difficult to study in clinical practice. The aim of the present work was to evaluate current dosing strategies and investigate potential dose adjustment of baricitinib and tofacitinib therapy of COVID‐19 in White and Chinese older adults, using our previously established PBPK modeling framework for age‐dependent simulations. ,

METHODS

Model development and verification

All PBPK modeling was performed using a population‐based absorption, distribution, metabolism, and excretion (simulator (Simcyp, version 19; Sheffield, UK)). Drug models were first developed based on two published PBPK models of baricitinib and tofacitinib. Key drug‐dependent parameters of baricitinib, tofacitinib, and probenecid necessary for PK simulation are reported in the original articles. Midazolam and simvastatin models were adopted from the Simcyp compound library and are presented in Table  as well.

Plasma concentration‐time profiles were simulated, and PK parameters were subsequently generated to verify the performance of the PBPK models. The participants, dose, and regimen selected for the simulations were matched to the corresponding clinical study designs. “Healthy Volunteers” and “Chinese Healthy Volunteers” provided in Simcyp were adopted as populations for the simulations for respective healthy subjects. A 10 × 10 trial was simulated each time. Model accuracy was assessed via fold error of the mean value of each simulated PK parameter to its corresponding observed clinical value.

Prospective simulations

Doses of baricitinib were selected based on its COVID‐19 EUA, which is 4 mg once daily (q.d.) for subjects with normal renal function. Dosage of 5 mg twice daily (b.i.d.) for rheumatoid arthritis from FDA prescribing information of tofacitinib were utilized for simulations. Several population models provided in Simcyp were utilized for prospective simulations of different clinical scenarios. “NEurCaucasian” (20–85 years) and “Chinese” (20–70 years) were used to represent the respective general populations. “Geriatric NEC” was used to simulate the White geriatric population (65–85 years), whereas the Chinese geriatric population model developed by Cui et al. was utilized for Chinese geriatric patients. The methods that characterized key parameters of each population model are summarized in Table  . For patients with chronic renal failure, the population models were further modified to recapitulate the different stages of renal dysfunction based on previous studies. , , Specifically, the default value of 60 million proximal tubular cells per gram kidney for a healthy population (creatinine clearance (CrCL) > 80 mL/min) was scaled down proportionally according to glomerular filtration rate (GFR) to 28.6, 17.6, and 8.8 million cells, for mild (CrCL 50–79 mL/min), moderate (CrCL 30–49 mL/min) and severe renal impairment (CrCL < 29 mL/min), respectively. Additionally, we applied a relative activity factor/relative expression factor of 0.73 and 0.41 for the baricitinib drug model in the moderate and severe renal failure populations to recapitulate the disproportionately greater deterioration in human organic anion transporter 3 (hOAT3)‐mediated secretion in comparison to GFR. , Such modification was not incorporated in the tofacitinib model due to the relatively minor contribution of renal transporters to the excretion of tofacitinib. Simulated area under the curve (AUC; mean value with 95% confidence interval (CI)) of each scenario was compared with that in healthy subjects to derive the AUC ratio for evaluating the potential overexposure of baricitinib and tofacitinib. Specifically, an AUC ratio of 1.84, obtained from the DDI simulation between baricitinib and probenecid in healthy subjects, was predefined as the upper limit of acceptable in vivo exposure of baricitinib.

RESULTS

Model verification and validation

As outlined in Table  , all simulated PK parameters fell within 1.6‐fold error and were comparable to the clinically observed values. , , , , The mechanistic PBPK models accurately recapitulated clinically observed PK profiles of both baricitinib and tofacitinib across different ethnicities and in the presence/absence of the hOAT3 inhibitor. Specifically, simulations of tofacitinib utilizing subjects with declined renal functions corroborated clinical observations of increase in tofacitinib systemic exposure. The Chinese geriatric population model was judiciously examined and verified via two probe substrates, simvastatin and midazolam.

Table 1

Observed and predicted pharmacokinetic parameters (mean ± SD) for baricitinib, probenecid, tofacitinib, midazolam, and simvastatin

Verification purpose	Substrate	Population	Dose	Parameter	Observed	Predicted	Fold error of mean	References
Baricitinib	Baricitinib	Healthy population	4 mg single	AUC0–∞ (h*ng/mL)	236.0 ± 51.9	223.8 ± 85.4	0.95	14
				C max (ng/mL)	36.2 ± 8.0	39.6 ± 12.9	1.09
				CL/F (L/hour)	16.9 ± 3.7	21.3 ± 10.0	1.26
				CLR (L/hour)	11.0 ± 2.4	12.3 ± 5.6	1.12
			4 mg single with probenecid	AUC0–∞ (hour*ng/mL)	480.0 ± 67.2	410.8 ± 121.7	0.86	14
				C max (ng/mL)	37.3 ± 7.5	46.3 ± 14.1	1.24
				CL/F (L/hour)	8.3 ± 1.2	10.8 ± 4.0	1.30
				CLR (L/hour)	3.4 ± 0.7	5.1 ± 1.5	1.49
		Chinese healthy population	4 mg single	AUC0–∞ (hour*ng/mL)	274.0 ± 45.2	245.3 ± 80.4	0.90	18
				C max (ng/mL)	51.3 ± 20.4	42.8 ± 13.8	0.83
				CL/F (L/hour)	15.0 ± 2.5	18.2 ± 6.5	1.22
	Probenecid	Healthy population	1000 mg b.i.d.	C max (µg/mL)	145.0	165.9	1.14	14
				AUCss (µg/mL)	115.0	122.4	1.06
Tofacitinib	Tofacitinib	Healthy population	10 mg single	AUC0–∞ (hour*ng/mL)	268.0 ± 71.5	261.5 ± 89.3	0.98	19
				C max (ng/mL)	94.2 ± 25.3	71.1 ± 24.8	0.75
		Subjects with mild RI	10 mg single	AUC0–∞ (hour*ng/mL)	370.0 ± 109.0	327.5 ± 105.3	0.89	19
				C max (ng/mL)	87.3 ± 23.2	74.5 ± 24.8	0.85
		Subjects with moderate RI	10 mg single	AUC0–∞ (hour*ng/mL)	396.0 ± 154.0	359.2 ± 142.9	0.91	19
				C max (ng/mL)	104.0 ± 47.5	75.4 ± 26.1	0.72
		Subjects with severe RI	10 mg single	AUC0–∞ (hour*ng/mL)	615.0 ± 214.0	512.8 ± 173.0	0.83	18
				C max (ng/mL)	111.0 ± 28.6	86.8 ± 31.8	0.78
		Chinese healthy population	10 mg single	AUC0–∞ (hour*ng/mL)	277.1 ± 37.0	354.3 ± 118.6	1.28	20
				C max (ng/mL)	105.3 ± 42.0	90.3 ± 28.3	0.86
Chinese geriatric population	Midazolam	Chinese elderly population (66‐75 years)	15 mg single	AUC0–∞ (hour*ng/mL)	229.0 ± 67.3	285.5 ± 243.5	1.25	16
				C max (ng/mL)	98.0 ± 41.5	89.2 ± 83.6	0.91
				CL/F (L/hour)	74.6 ± 25.7	103.7 ± 101.1	1.39
		Chinese elderly population (>76 years)	15 mg single	AUC0–∞ (hour*ng/mL)	254.0 ± 87.3	319.5 ± 269.0	1.26	16
				C max (ng/mL)	96.8 ± 57.5	97.3 ± 75.8	1.01
				CL/F (L/hour)	66.0 ± 23.2	95.5 ± 88.1	1.45
	Simvastatin	Chinese elderly population (60–91 years)	20 mg QD	AUC0–∞ (hour*ng/mL)	17.2	22.1	1.29	16
				CL/F (L/hour)	1020.0	1630.8	1.60

AUC0 – ∞, area under the curve from zero to infinity; AUCss, area under the curve at steady‐state; C max, peak plasma concentration; CL/F, apparent clearance; CLR, renal clearance; RI, renal impairment.

Simulations of baricitinib among White and Chinese populations with different stages of chronic kidney disease are illustrated in Figure  a. Most of the AUC ratios were lower than the upper limit of 1.84, corroborating that current dosage adjustments based on GFR and in the presence of a strong hOAT3 inhibitor were adequate for both general and older populations. The AUC elevation resulting from renal dysfunction was more pronounced for geriatrics than general populations. Notably, for older White subjects with mild renal impairment, current recommendation of 4 mg dose of baricitinib could lead to a 2.31‐fold AUC increase (95% CI 2.10–2.55), and reduction to 2 mg could normalize the ratio to an acceptable 1.15‐fold (95% CI 1.05–1.27). Although a comparable trend of AUC fold‐change being correlated with renal dysfunction was observed for the Chinese population, the overall AUC increase was less pronounced than in the White population for both general and older populations. In the Chinese geriatric population with mild renal impairment receiving 4 mg dose, the AUC ratio of baricitinib was elevated to an acceptable 1.68‐fold (95% CI 1.52–1.85) and no dose adjustment simulation was performed.

Figure 1

Simulated area under the curve (AUC) ratio of baricitinib (a) and tofacitinib (b) normalized to that in healthy subjects across different ethnicities, age groups, renal functions, and in the presence/absence of probenecid. Data are presented as mean value with 95% confidence interval (CI). Healthy subjects: Simcyp “Healthy Volunteers” or “Chinese Healthy Volunteers” population; General population: Simcyp “NEurCaucasian” (20–85 years) or “Chinese” (20–70 years) population; Geriatric population: Simcyp “Geriatric NEC” population or Chinese geriatric population based on Cui et al.’s work. b.i.d., twice daily; GFR, glomerular filtration rate; q.d., once daily; RI, renal impairment.

Similar trends of increased systemic exposure with declined renal function were observed for tofacitinib, among both White and Chinese populations (Figure  b). As expected, older subjects yielded a higher AUC ratio of tofacitinib when compared with the general population. AUC ratios were 1.52 (95% CI 1.37–1.69) and 1.27 (95% CI 1.13–1.42) for older and general White adults with moderate renal impairment, respectively. Recommended dose reduction from 5 mg b.i.d. to 5 mg q.d. for older White subjects with severe renal impairment normalized AUC ratio to a comparable level observed with normal renal function (AUC ratio of 0.97, 95% CI 0.87–1.08).

DISCUSSION

Both baricitinib and tofacitinib are partially eliminated by the kidneys. Baricitinib is predominantly eliminated unchanged in urine (70% of dose) where its renal tubular secretion is mainly mediated by basolaterally expressed OAT3 and several apically expressed transporters. , Probenecid, a strong hOAT3 inhibitor, increased the AUC of baricitinib by 2‐fold and decreased its renal clearance (CLR) to 69% of control in healthy subjects. Our PBPK modeling reproduced the reduced CLR of baricitinib when it was co‐administered with probenecid. According to the dosage adjustment stipulated by EUA of baricitinib, the dose should be reduced from 4 mg q.d. to 2 mg or 1 mg q.d. according to GFR and be further halved when co‐administered with strong hOAT3 inhibitors. Our simulations among general White and Chinese adult populations rationalized these dosing recommendations. Current prescribing information of baricitinib for rheumatoid arthritis treatment recommends 2 mg q.d. as the baseline dose and does not recommend for its use in patients with severe renal impairment, which are different from the EUA instruction. Our simulation results revealed that a dose reduction strategy could be applied to this patient cohort as well for adequate treatment of rheumatoid arthritis.

Renal excretion accounted for 30% of the elimination of tofacitinib. Both the FDA and European Medicines Agency (EMA) recommend dose adjustment from 5 mg b.i.d. to 5 mg q.d. of tofacitinib for subjects with severe renal impairment and the FDA recommends the same dosage reduction for moderate renal impairment. , Active renal secretion may contribute to renal elimination of tofacitinib, whereas the potential transporter has not been identified yet. A PBPK simulation study by its manufacturer suggested limited contribution of active CLR secretion to overall clearance of tofacitinib. The PBPK model was then adopted for our investigation and recapitulated the clinical‐observed increase in AUC of tofacitinib when renal function progressed from normal to severe. The overall increase in AUC of tofacitinib due to renal impairment was relatively less pronounced than that observed for baricitinib. This is not unexpected due to the relatively minor contribution of renal elimination to the total clearance of tofacitinib.

The COVID‐19 pandemic underscored the challenge of under‐representation of older adults in clinical trials, especially those living in long‐term care facilities. A recent analysis of drug trials for COVID‐19 concluded that 23% excluded older adults based on a chronologic age restriction, and an additional 53% had indirect age‐related exclusions for comorbidities, functional impairments, lack of access to the internet or information technology, or other broad, poorly defined, or supported exclusions. Addressing this challenge by the FDA, a workshop hosted in 2020, entitled “Roadmap to 2030 for New Drug Evaluation in Older Adults,” thoroughly discussed the deficits in the evaluation of drugs in clinical trials for older adults who are often accompanied by multiple chronic health conditions, polypharmacy, or frailty. The under‐representation of this cohort sets similar barriers for making adequate decisions on optimal drug therapy of patients with COVID‐19.

Physiologically‐based pharmacokinetic modeling has previously demonstrated its predictive power by simulating the impact of aging in defining the PK of drugs in White and Chinese geriatric populations. , A combination of measured in vitro and clinically observed in vivo data could be applied to predict complex scenarios in young and older adults for dose adjustment recommendations. Using PBPK modeling, our laboratory reported that rivaroxaban, a substrate of hOAT3, could have a clinically significant increase in systemic exposure when it is co‐administered with amiodarone among older patients with moderate renal impairment, and proposed further dose reduction of rivaroxaban to ensure its efficacy and safety. , Our subsequent sensitivity analysis found that age and hOAT3 activity jointly contributed to the extent of this DDI and elevated systemic exposure of rivaroxaban.

Here, we applied validated PBPK models to simulate the PKs of baricitinib and tofacitinib in White and Chinese older subjects with different severities of renal impairment. Our results recapitulated the augmented AUC ratios of both JAK inhibitors in the presence of renal dysfunction where the older subjects are shown to be more susceptible. Extra caution is therefore needed during the treatment of COVID‐19 for geriatric patients with renal impairment using these drugs.

Although single doses up to 40 mg and multiple doses of up to 20 mg daily for 10 days of baricitinib have been administered in clinical trials without dose‐limiting toxicity, concomitant probenecid would require dose modification of baricitinib due to its 2‐fold increase in systemic exposure. Our simulation results observed similar risk of overexposure when older adults with mild renal impairment received baricitinib (AUC ratio of 2.31). A 2 mg q.d. dose of baricitinib, which is not currently listed as a recommended dosage adjustment, yielded an acceptable AUC ratio of 1.15. For tofacitinib, caution has been suggested for its geriatric use due to a higher incidence of infections in the older population. Our simulation illuminated another potential risk where there was a simulated 1.52‐fold increase in the AUC of tofacitinib when it was prescribed at 5 mg b.i.d. to older White subjects with moderate renal impairment. The increase in AUC was expectedly higher than those observed for general adults, and substantiated the dosage reduction of tofacitinib recommended by the FDA.

Exposure of baricitinib and tofacitinib in Chinese subjects were less affected by renal dysfunction as compared with their White counterparts, especially among older adults. AUC only increased by 1.68‐fold (95% CI 1.52–1.85) in older subjects with mild renal impairment receiving 4 mg q.d. baricitinib, and 1.36‐fold (95% CI 1.23–1.51) in those with moderate renal impairment receiving 5 mg b.i.d. tofacitinib. Notably, the age‐related distribution of physiological parameters was different between the Chinese and White geriatric populations (Table  ). Although geriatric populations were set as 65–85 years in our simulations, mean age distribution for the White population was 75.9 years (95% CI 66.7–84.1 years), whereas a shift was observed for Chinese population with less variability (mean 69.3 years, 95% CI 65.3–76.3 years). Such mean age distribution disparity possibly yielded differences in key renal function parameters, such as kidney weight and serum creatinine. Another possibility is that there is currently no validated PBPK model for the Chinese geriatric population with chronic renal disease and parameters were derived from White population models. Taken together, further optimization and verification of the renal impairment population model for Chinese geriatric patients are warranted.

The geriatric population with chronic organ dysfunction is one of the most vulnerable cohorts when combating the COVID‐19 pandemic, for whom treatment indications might not be well‐represented in proportion to the prevalence of infection. Our proposed dose reduction of baricitinib in White older patients with mild renal impairment provided prospective evidence for the rational dosing of this cohort. In conclusion, PBPK modeling of both JAK inhibitors supports the rational and prudent dose adjustments of these COVID‐19 therapeutics among adult patients of different age groups and renal functions.

FUNDING

Z.W. received postdoctoral fellowship from the National University of Singapore. E.C.Y.C. received research grant funding from the Joseph Lim Boon Tiong Urology Cancer Research Initiative (Grant: A‐0002678‐01‐00).

CONFLICT OF INTEREST

The authors declared no competing interests for this work.

AUTHOR CONTRIBUTIONS

Z.W. and E.C.Y.C. designed the research and wrote the manuscript. Z.W. performed the research and analyzed the data.

Table S1

Click here for additional data file.

Table S2

Click here for additional data file.