Dosing optimization of rituximab for primary membranous nephropathy by population pharmacokinetic and pharmacodynamic study

Primary membranous nephropathy (PMN) is the most common cause for adult nephrotic syndrome. Rituximab has demonstrated promising clinical efficacy by random controlled trials and the off-label use is widely adopted in PMN. However, the standard dosage is borrowed from B cell lymphoma treatment with far more antigens and is oversaturated for PMN treatment, accompanied with additional safety risk and unnecessary medical cost. More than 15% serious adverse events were observed under standard dosage and low dose therapies were explored recently. Dose optimization by clinical trials is extremely time- and cost-consuming and can be significantly accelerated with the aid of model-informed drug development. Here, we aim to establish the first population pharmacokinetic and pharmacodynamic (PPK/PD) model for rituximab in PMN to guide its dosage optimization. Rituximab pharmacokinetic and pharmacodynamic data from 41 PMN patients in a retrospective study under a newly proposed monthly mini-dose were used to construct quantitative dose-exposure-response relationship via mechanistic target-mediated drug disposition (TMDD) model followed by regression between the reduction of anti-PLA2R titer and time after the treatment. The final model, validated by goodness-of-fit plots, visual predictive checks and bootstrap, was used to recommend the optimized dosing regimen by simulations. The model was well validated for PK/PD prediction. The systemic clearance and half-life are 0.54 L/h and 14.7 days, respectively. Simulation of a novel regimen (6 monthly doses of 100 mg) indicated the comparable ability and superior duration time of CD20⁺ B cell depletion compared with standard dosage, while the cumulative dosage and safety risk was significantly decreased. We established the first PPK/PD model and provide evidence to support the dosage optimization based on monthly mini-dose. Our study can also efficiently accelerate dosage optimization of novel anti-CD20 antibodies in PMN and other indications.

1 Introduction

Membranous nephropathy (MN) is the one of the most predominate causes of adult nephrotic syndrome (NS). About 80% of MN patients are primary MN (PMN) with no underlying cause, while one-third PMN patients progress to advanced kidney failure (Ronco et al., 2021). Podocyte phospholipase A2 receptor (PLA2R) was identified as the major antigen in PMN responsible for the formation of immune complex deposition on basement membranes that induces podocyte injury and proteinuria (Beck et al., 2009). The pathophysiological role of PLA2R provides a clear rational that the decrease of anti-PLA2R titer precedes the remission of proteinuria (Ruggenenti et al., 2015). Thus, PMN remission can be predicted by anti-PLA2R antibody reduction whereas relapse occurred with antibody re-emergence.

As an autoimmune disease, immunosuppressive therapies are widely used to treat PMN. However, traditional immunosuppressive treatments have long been suffered from high incidence of adverse effect and relapse rate. Rituximab (RTX), a chimeric monoclonal anti-CD20 antibody, eliminates CD20⁺ B cells through antibody-mediated cytotoxicity (ADCC) and complement-mediated cytotoxicity (CDC) effect (Boross and Leusen, 2012). RTX was primarily approved for the treatment of B cell lymphoma and then further extended to autoimmune disease. RTX offers a selective B cell targeting approach that prevent antibody production, and thus was introduced into PMN treatment (Remuzzi et al., 2002). Random controlled trials demonstrated that RTX was noninferior to cyclosporine in terms of 12-month proteinuria remission and was superior in the duration of proteinuria remission up to 24 months (Dahan et al., 2017; Fervenza et al., 2019). Accordingly, RTX has been suggested as first-line treatment for PMN by the latest KIDGO guideline (Rovin et al., 2021).

Although the off-label use of RTX was widely adopted for PMN treatment, the optimal dosage is still in debate. The dosage was borrowed from B cell lymphoma treatment and recommended as four weekly infusions of 375 mg/m² or two 1000 mg infusions with 2 weeks apart (Rovin et al., 2021). However, the B cell counts in PMN patients are only 1/300 of that in lymphoma patients., Thus, the standard dose may be oversaturated for B cell depletion in PMN patients even with the urinary loss of rituximab and cause unnecessary medical cost (Bensalem et al., 2022), though the safety risk is relatively low when compared with other immunosuppressive drugs (Fervenza et al., 2019). It is still debatable for the relationship between larger doses and improved clinical outcomes. For example, studies by Gabriella Moroni et al. and Roberta Fenoglio et al. indicated that the efficacy is similar under low dose or high dose of rituximab (Moroni et al., 2017; Fenoglio et al., 2021). Although some studies demonstrated that high dose therapy yield higher remission rate (Seitz-Polski et al., 2019), further analysis indicated that the remission is associated with residual rituximab levels at month-3 (Boyer-Suavet et al., 2019), which can serve as an effective predictor for remission (Teisseyre et al., 2021). Accordingly, it is more effective to maintain the residual rituximab level for a longer time than to have high peak concentration. A retrospective study was carried out by our group and revealed that monthly mini-dose (100 mg) has better short-term (6 months) clinical efficacy compared to the standard dose regimens while the cumulative dose and the rate of adverse event is much lower (Wang et al., 2023). Inspired by the promising preliminary results, further studies are urgently warranted to compare with standard dosage and better optimize the dosage.

It is very time- and cost-consuming to explore the optimal dosage by conducting clinical trials, especially for PMN that requires long-term (up to 24 months) observations. Model-informed drug development strategy is a powerful tool to aid and accelerate dose optimization by constructing the quantitative dose-exposure-response relationship and accurately predicting the clinical efficacy of different regimens in advance (Cardone et al., 2011). However, the pharmacokinetic (PK) characteristics were rarely described for RTX in PMN treatment and the pharmacokinetic and pharmacodynamic (PK/PD) relationship of RTX has not been established yet (Bensalem et al., 2022). Although CD20⁺ B cell depletion is not the clinical endpoint for the treatment of PMN, it can serve as an efficient surrogate that bridge the PK of RTX and clinical efficacy. Here, we constructed the first population PK/PD (PPK/PD) model that revealed the quantitative relationship between dosage, RTX concentration, CD20⁺ B cell count and anti-PLA2R titer, and enabled us to effectively optimize dosage according to simulation results. Our study provides concrete evidence to find the optimal dose with significantly reduced SAE occurrence and medical burden for PMN patients on the basis of monthly mini-dose.

2 Materials and methods

2.1 Patient population and study design

In our retrospective study, a total of 41 PMN patients who received RTX treatment from March 2019 to December 2021 at Department of Nephrology, Peking University Third Hospital were included. Most patients received monthly mini-dose of 100 mg RTX, whereas some of them received 200–500 mg RTX for a few months with the intention to better control disease progress. Due to COVID-19 epidemic, part of patient population experienced a gap up to 10 months and received RTX infusion again afterwards. The trial was consistent with the Good Clinical Practice (GCP) standards and was registered at Chinese clinical trial registry (ChiCTR2200057381).

2.2 PK and PD data

The trough and peak RTX blood samples were collected before and after RTX infusion, respectively. RTX concentrations were detected by ELISA kit purchased from Abcam (Cambridge, MA, United States). RTX PK observations below the limit of quantification (3 ng/mL) were excluded from the analysis. CD20⁺ B cell count was measured before every RTX intravenous (IV) infusion. CD20⁺ B cell depletion was defined as an absolute count < 5/μL in peripheral blood (Cravedi et al., 2007). Serum anti-PLA2R titer was assessed by an ELISA kit purchased from Euroimmune, Lubeck, Germany.

2.3 PPK/PD model development

PPK/PD models were constructed by nonlinear mixed-effects modeling methods to describe the PK characteristics and variability for RTX concentration and CD20⁺ B cell count (Supplementary Figure S1). The observed RTX concentrations were converted to molar units. The CD20⁺ B cell counts were transformed to CD20 concentration in molar units by multiplying 94000 CD20 molecules per cell and dividing by Avogadro’s constant (Ginaldi et al., 1998).

A total of four structural PPK models were developed as follows. At first, a basic two-compartment model was constructed as reported by many previous studies that is characterized by linear elimination (Ng et al., 2005; Regazzi et al., 2005; Li et al., 2007; Blasco et al., 2009). Based on the basic two-compartment model, two additional models were constructed with nonlinear elimination approximation described by Michaelis-Menten elimination (Tout et al., 2017) and time-varying elimination (Li et al., 2012), respectively. Finally, a mechanistic target-mediated drug disposition (TMDD) model with quasi-steady-state (QSS) approximation was constructed. In this model, the free RTX in central compartment is directly eliminated by a first-order elimination rate constant (k_el). In addition, central compartment RTX distributes to and from peripheral compartment by first-order distribution rate constant k_pt and k_tp, respectively. In plasma, RTX binds to CD20 antigen (R) with a second-order rate constant (k_on) to form drug-receptor complex (RC), whereas RC dissociates at a first-order rate (k_off) and internalizes for degradation (k_int). CD20 is reported to be a non-internalizing receptor, and the elimination of CD20 and the subsequent depletion of CD20⁺ B cells is thought to be via ADCC and CDC effect (Boross and Leusen, 2012). Thus, the internalization rate constant k_int was replaced by the RTX-CD20 complex target-mediated elimination rate constant (k_tmd), which is assumed to be the maximal rate of elimination observed for mAbs in clinical use (Glassman and Balthasar, 2017). On the other hand, CD20⁺ B cell is synthesized at a zero-order rate (k_syn) and degraded at a first-order rate (k_deg) in the absence of RTX binding. Due to its non-internalization character, the turnover rate of CD20 was assumed to be similar to the disappearance rate of B cells from circulation (∼3.9% per day) (Macallan et al., 2005). The RTX-CD20 complex elimination rate constant (k_tmd) is much greater than the dissociation rate constant k_off (Bondza et al., 2017). Thus, QSS approximation of TMDD model was used, where QSS constant k_ss is defined as the sum of dissociation constant k_D and the quotient of k_tmd and dissociation rate constant (k_off) (Gibiansky et al., 2008). Models were evaluated with objective function value (OFV), while the model with the lowest OFV and better fitting of the diagnostic plots was selected as structural model.

Between-subject variability (BSV) was modeled in an exponential form (Eq. 1):

$P_i=TVP\times e^{{\eta }_i}(1)$

where $P_i$ is the ith individual value of the parameter; TVP is the typical value of the population parameter; η_i is the empirical Bayes estimates of BSV for the ith individual that is normally distributed with mean zero and variance ω².

The residual variabilities were assessed as an error model in a proportional form (Eq. 2):

$C_{ij}=\stackrel{\wedge }{C_{ij}}\times \left(1+{\epsilon }_{pro}\right)(2)$

where $C_{ij}$ and $\stackrel{\wedge }{C_{ij}}$ is the ith observed and predicted concentration for the jth individual, respectively; ${\epsilon }_{pro}$ is the proportional error component that is normally distributed with mean zero and variance σ².

2.4 Covariate model development

Covariates with high correlation (correlation coefficient > 0.3) were taken into consideration for further the covariate analysis, including age, sex, eGFR, etc. Continuous and categorical covariates were established by a power function and a linear function (Eqs 3, 4), respectively:

$P_i={\theta }_1\times {\left(\frac{Co{v}_{con}}{Co{v}_{median}}\right)}^{{\theta }_2}(3)$

$P_i={\theta }_1\times \left(1+{\theta }_2\times Co{v}_{cat}\right)(4)$

Where P_i is the parameter for the ith individual; $Co{v}_{con}$ and $Co{v}_{cat}$ is the value of continuous and categorical covariate for the ith individual, respectively; ${\theta }_1$ is the typical value for the population parameter; ${\theta }_2$ is the estimated value for the covariate effect.

If two or more covariates were highly correlated, only the most clinically and practically relevant covariate was tested on the PK parameter in the analysis. These covariates were further screened using the SCM implemented in PsN.

2.5 Model evaluation

Goodness-of-fit (GOF) plots were used to evaluate the final model, which defined as the combination of four plots including the observed concentration (DV) versus the population predicted concentration (PRED) or individual predicted concentrations (IPRED) and conditional weighted residuals (CWRES) versus PRED or time. The predictive performance of the final model was evaluated using visual predictive checks (VPCs), which were repeated 1000 times to compare the DV data and simulated value with the 95% confidence intervals for the 5^th,50th and 95th percentiles.

The results of bootstraps with 500 replicates were finally used to assess the final model. The RSEs of the parameter estimates and the condition number of the model were also considered.

2.6 Simulations

Using the established model, we carried out simulations for both standard dosage and monthly mini-dose to predicted RTX concentrations and CD20⁺ B cell counts under different regimen groups. In addition, mini-dose with larger infusion interval was also simulated.

2.7 Exponentially decrease function of anti-PLA2R titer

To explore the relationship between CD20⁺ B cell counts and anti-PLA2R titer, the individual anti-PLA2R titer data at descending stage was extracted, while serology relapse data that was defined as two consecutive ascending anti-PLA2R titer data was excluded. Anti-PLA2R titers detected with interval longer than 3 months or after 12 months since the first RTX infusion were excluded as well. The anti-PLA2R titer data was plotted versus time and described using an exponentially decrease function (Eq. 5). For patients experienced relapse and second round of anti-PLA2R reduction, two exponentially decrease functions using respective data were used.

$anti-PLA2Rtiter=A\times e^{-k_{e,PLA2R}\times t}(5)$

2.8 Non-compartmental analysis

Non-compartmental analysis was conducted for the simulated results using typical PK parameters. The area under the concentration-time curve (AUC), the elimination half-life (t_1/2), the apparent clearance (CL/F) and apparent volume of distribution during the terminal phase Vz (Vz/F) were determined.

2.9 Software

All the modeling and simulation were carried out with the first order conditional estimation method with interaction that is implemented within NONMEM 7.2 and aided by PsN. The dataset was treated and processed by R package. Pirana was used for documentation of the development process. Non-compartmental analysis was conducted using Phoenix WinNonlin V8.3.3 (Certara, Princeton, NJ, United States).

3 Results

3.1 Demographic characteristics and clinical response

A summary of subject demographic characteristics is shown (Table 1). A total of 31 male (75.6%) and 10 female (24.4%) patients were enrolled. The ages of the patients ranged from 19 to 76 years, with weight ranging from 50 to 101 kg (75.3 ± 11.3 kg). The albumin concentration and urine protein at baseline is 25.8 ± 5.7 g/L and 8.0 ± 3.6 g/day, respectively. The anti-PLA2R antibody titer ranged from 5.4–2695 U/mL (260.3 ± 453.2 U/mL), while 17 patients had a high anti-PLA2R antibody titer (>150 U/mL). The ratio of CD19⁺ B cell depletion at month 3 was 80% (33/41). The mean follow-up period was 15.9 (6–44) months. The cumulative dose of rituximab at month 6, month 12 and last follow-up was 578 ± 350, 923 ± 544 and 1060 ± 705 mg, respectively. The urinary protein levels decreased and serum albumin levels increased gradually after rituximab treatment, while the eGFR remained relatively stable. The remission rate was 21.2% (7/33) with 0% complete remission (CR) at 3 months, 50% (19/38) with 5.3% CR (2/38) at 6 months, 72.7% (24/33) with 18.2% CR (6/33) at 12 months, and 81.6% (30/38) with 28.9% (11/38) CR at last follow-up, respectively. Rate of relapse was 7.9% (3/38) at last follow-up.

Table 1

Table 1. Demographic characteristics and clinical response of the study population at baseline and after treatment.

3.2 PK and PD data

A total of 171 RTX concentration, 220 CD20⁺ B cell count and 276 PLA2R titer data from 41 PMN patients were collected. Most patient receive multiple doses of 100 mg RTX infusion with at least 1 month apart, while part of them received 200–500 mg RTX infusion.

3.3 PPK/PD model construction

Four PPK models were constructed to describe RTX PK profile, including two-compartment model with linear, time-varying and Michaelis-Menten. An additional PPK/PD model described by TMDD model with QSS approximation was established. Among four models, TMDD model best captured the RTX concentration and CD20⁺ B cell count profile with lowest OFV (Supplementary Table S1), and was thus selected as structural model. The model has included both specific elimination (k_tmd) and non-specific elimination (k_el). The contribution of both proteolysis and urinary loss have been considered in our model, although it is difficult to exactly evaluate their respective contribution due to limited data.

3.4 Covariate model and final model

More than 30 baseline physiological parameters were used for covariates screening. Notably, covariates have been identified for RTX in many other diseases, such as body weight, sex, age on clearance and body surface area on central volume (Wang et al., 2020). However, no covariates were identified in our study. Thus, the final model was replaced by the base structural model and the results of estimated model parameters were presented (Table 2).

Table 2

Table 2. Estimated PPK/PD final model parameters.

3.5 Model evaluation

The goodness-of-fit diagnostic plots of the final model were shown (Figures 1, 2). There was no obvious bias in the model fitting. The VPC and bootstrap results suggested that the model adequately describe the both the RTX concentrations and CD20⁺ B cell counts (Supplementary Figure S2). The general trend and the observed variability were well captured.

Figure 1

Figure 1. Goodness-of-fit plots for RTX of the final model. The red line is generated by linear method and shows the trend of the data, and the blue circles represent the observed data.

Figure 2

Figure 2. Goodness-of-fit plots for CD20 of the final model. The red line is generated by linear method and shows the trend of the data, and the blue circles represent the observed data.

3.6 Exponentially decrease function of anti-PLA2R titer

Mechanically, anti-PLA2R antibody is excreted by plasma cells, which is in turn transformed from activated B cells that express CD20, and eliminated through catabolism. Thus, once CD20⁺ B cells are depleted by RTX, the plasma B cell count slowly decreases and the production of anti-PLA2R is halted, leading to the subsequent decrease of anti-PLA2R titer. However, we failed to establish the quantitatively relationship between CD20⁺ B cell counts and anti-PLA2R titers. The major reason is that monthly (or even longer) monitor of CD20⁺ B cell counts is incapable to capture the descending process of CD20⁺ B cells as they are rapidly depleted in peripheral blood within 72 h after receiving RTX. Thus, the majority of CD20⁺ B cell count data was obtained either before RTX administration or after depletion (mostly 0/μL) (Supplementary Figure S3A). In addition, the relationship between the recovery of CD20⁺ B cell count and anti-PLA2R titer is controversial. The reemergence of CD20⁺ B cell may be followed either by the increase of anti-PLA2R titer (Supplementary Figure S3B) or not (Supplementary Figures S3C, D). Accordingly, these data are also helpless for modeling.

The anti-PLA2R titers of 36 in all 41 patients were analyzed without the integration of CD20⁺ B cell count. Notably, we found that the anti-PLA2R titer data can be well described using an exponentially decrease function whose R² > 0.8 for all patient except one (R² = 0.67). The function is parameterized with initial anti-PLA2R titer and individual elimination constant k_e,PLA2R, which is independent of CD20⁺ B cell counts. The mean k_e,PLA2R is 0.033 ± 0.017, corresponding to the mean half-life 21 days, similar with typical IgG4 antibody half-life. Similar phenomena were also observed in previous studies treated with standard dosage, where the k_e,PLA2R is 0.025 with a R² of 0.98 using titers in the first 3 months (Ruggenenti et al., 2015). We conducted Pearson correlation analysis to explore the physiological factor that influence k_e,PLA2R. Unfortunately, no significant relationship (Pearson correlation coefficient > 0.7) was identified.

3.7 Simulation of different dosage

The simulations were carried out to compare the efficacy of standard dosage 1 (four weekly infusions of 375 mg/m2), standard dosage 2 (two 1000 mg infusion with 2 weeks apart) and mini-dose (six monthly 100 mg). RTX rapidly achieve peak concentration at 215, 332, and 39 μg/mL, whereas the trough concentration collected 1 week after first infusion is 52, 118 and 7.8 μg/mL for standard dosage 1, standard dosage 2 and mini-dose, respectively (Figure 3). The apparent clearance and half-life are 0.54 L/day and 14.7 days, respectively, as calculated by non-compartmental analysis for monthly mini-dose.

Figure 3

Figure 3. Simulation of RTX concentration after infusion of (A) four weekly 375 mg/m² RTX; (B) two 1000 mg RTX with two weeks apart; (C) six monthly 100 mg RTX; (D) three 100 mg RTX with two months apart. The solid line represents the mean RTX concentration versus time, while the shadow represents the 90% prediction intervals.

CD20⁺ B cell depletion was used as a surrogate for clinical efficacy, since the decrease of anti-PLA2R titer occur shortly after the depletion and is independent of B cell recovery as demonstrated above. One the other hand, B cell recovery can be regarded as a marker for safety since it is related with reduced infection incident, which is the predominate adverse effect of RTX in PMN treatment. Though the concentration is much lower for mini-dose, the ability to deplete CD20⁺ B cell is similar (Figure 4). Under all three dosages, the depletions of CD20⁺ B cells occur within 24 h, and slowly recover after last dose. Notably, the duration of depletion is longer for mini-dose than standard dose. CD20⁺ B cell count increases back to more than 5/μL in ∼ 5.5 months for both standard dosages after first infusion, while the depletion can last for more than 7 months for mini-dosage. On the other hand, it takes shorter time for mini-dosage to recover back to normal (90% of baseline CD20⁺ B cell count) than standard dosages (10 months vs 12 months after last dose). Thus, the mini-dose outperforms standard dose with lower total doses (600 mg vs 2000–2400 mg), longer depletion duration (7 months vs 5.5 months) and faster recovery (10 months vs 12 months).

Figure 4

Figure 4. Simulation of CD20⁺ B cell counts after infusion of (A) four weekly 375 mg/m² RTX; (B) two 1000 mg RTX with two weeks apart; (C) six monthly 100 mg RTX; (D) three 100 mg RTX with two months apart. The solid line represents the mean CD20⁺ B cell counts versus time, while the shadow represents the 90% prediction intervals. Dashed line indicated CD20⁺ B cell depletion criterion (5/μL).

To explore the possibility to further reduce the dosage with larger interval, we simulated the RTX concentration and CD20⁺ B cell count receiving 3 doses of 100 mg every 2 months. Although the RTX concentrations and time for CD20⁺ B cell depletion remain the same with monthly mini-dose in the first week, the duration of CD20⁺ B cell depletion is largely reduced. The median CD20⁺ B cell count increases up to more than 5/μL in less than 2 months. Thus, 100 mg every 2 months dosage is insufficient for long time CD20⁺ B cell depletion and thus should not be recommended for PMN treatment.

4 Discussion

PPK model have been established to describe the PK profile of RTX under various disease conditions, including non-Hodgkin’s lymphoma (NHL) (Ternant et al., 2019), diffuse large B cell lymphoma (DLCBL) (Blasco et al., 2009), chronic lymphocytic leukemia (CLL) (Li et al., 2007), follicular lymphoma (FL) (Regazzi et al., 2005), rheumatoid arthritis (RA) (Ng et al., 2005) and anti-neutrophil cytoplasmic autoantibody-associated vasculitis (AAV) (Bensalem et al., 2020). However, PK profile of RTX in PMN patients may be different due to distinct CD20⁺ B cell count (in comparison with B cell lymphoma) or additional elimination of antibody through proteinuria. Indeed, previous study indicated that the clearance of RTX in PMN is ∼2 folds than in follicular lymphoma and autoimmune disorders (Fogueri et al., 2019). Similar PK profile was also observed in our study. The systematic clearance in our study is 0.54 L/h, same as previously reported RTX clearance in PMN and higher than that in B cell lymphoma and other autoimmune disease (Bensalem et al., 2022). The central and peripheral volumes of distribution is 2.48 L and 4.68 L, respectively. The central volume of distribution in our study is slightly lower than that in other disease but 1.9-fold higher than that in previous PMN study. The half-life of a drug means the time it costs for the body to eliminate half of the remaining drug. Shorter half-life is associated with the faster decrease of RTX concentration and may lead to reduced clinical efficacy. The RTX half-life in our study (14.7 days) is comparable to that in previous kidney disease study, such as PMN (11.5 days) (Fogueri et al., 2019), pediatric primary glomerulonephritis (14.6 days) (Zachwieja et al., 2012) and adult minimal change nephrotic syndrome (10–15 days) (Iwabuchi et al., 2018) and is much shorter than that in follicular lymphoma and autoimmune disorders (Regazzi et al., 2005). These evidences indicated that the PK and PD profiles are significantly affected by the urinary loss of RTX, which is consistent with previous studies (Del Vecchio et al., 2021; Teisseyre et al., 2022).

Two-compartment model is frequently used to describe RTX PK in most studies, a simple approach but ignores the nonlinear elimination due to antigen binding. It was recently demonstrated that TMDD model can better capture the PK profile of rituximab (Bensalem et al., 2022). In our study, we constructed four two-compartment models with linear or nonlinear elimination described by time-varying, Michaelis-Menten or TMDD approaches. Consistent with previous study, TMDD model outperform other models with lowest OFV and better fitting. In addition to RTX concentration, TMDD is capable to capture the profile of CD20⁺ B cell count, a PD marker, providing the basis to construct the PK/PD relationship to guide dosage optimization.

CD20⁺ B cell plays a critical role in PMN progress and prognosis. A prospective study was conducted to compare risk and benefit profile of standard four weekly 375 mg/m² RTX dose treatment or B cell-driven treatment (patients receive a second infusion of RTX only if they had CD20⁺ B cell count more than 5/μL). B cell-driven treatments achieved similar reduction of proteinuria with less adverse events and much lower medical costs (Cravedi et al., 2007). A recent study revealed that serum RTX concentration at month-3 (3 month after RTX injection) can well predict the remission rate at month-6 and month-12. Patients with RTX concentrations less than 2 μg/mL at month-3 exhibited weaker CD20⁺ B cell depletion ability, lower clinical remission rate and need longer time to achieve clinical remission (Teisseyre et al., 2021). Collectively, these studies indicated the rational and reliability to use the duration of CD20⁺ B cell depletion as a marker to indicate the therapeutic effect and to balance risk and benefit.

Although the PK/PD model was not established for anti-PLA2R titer due to the sparse sampling of CD20⁺ B cells, an exponentially decrease function was used to predict the reduction of anti-PLA2R. The mean elimination constant k_e,PLA2R is 0.033 ± 0.017, corresponding to the mean half-life of 21 days. Six monthly doses of 100 mg RTX infusion is able to maintain the depletion of CD20⁺ B cells for more than 7 months, which is sufficient for the fully elimination of anti-PLA2R antibody (∼10 half-lives). However, a large variability was observed for k_e,PLA2R that ranges from 0.01 to 0.079. For some patients, 6 months infusion may be insufficient or unnecessary. Although the factors that influence the individualized elimination constant k_e,PLA2R remains elusive, monitoring of anti-PLA2R titer can be used to fit the parameter. After that, the reduction of anti-PLA2R is predictable and dosage can be adjusted accordingly. Therefore, a personalized treatment is feasible, which will be further explored in the upcoming study.

There are still many defects in our study. First, the sample size is still limited, especially when sparse sampling was adopted. Second, we failed to identify any covariates that influence both RTX concentration clearance and anti-PLA2R antibody elimination. Previous studies have identified body weight, sex and body surface area as covariate for clearance under other disease conditions (Wang et al., 2020). For PMN, however, all these covariates have little influence on RTX PK profile. No covariate was reported by other PPK study for RTX usage in PMN as well (Fogueri et al., 2019). More clinical studies will be conducted in the near future to further optimize our model and dosage.

5 Conclusion

In this study, we established the quantitative relationship between RTX dosage, PK characteristics, CD20⁺ B cell counts and anti-PLA2R titers by population PK/PD model. The simulation results supported the use of monthly mini-dose for the treatment of PMN patient with comparable efficacy and much lower cost. Our study can also provide insightful view for the clinical development and precision medicine of other anti-CD20 antibodies.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by Ethics Institute of Peking University Third Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

HoL, DZ, DL designed the study. ZD performed the experiment and curated the data. SN, WK, and YL analyzed the data. HoL and WK wrote the manuscript. SW, HiL and YW critically reviewed the manuscript. All authors contributed to the article and approved the submitted version.

Funding

The study is supported by Bill & Melinda Gates Foundation (INV-007625), the Clinical Cohort Construction Program of Peking University Third Hospital (BYSYDL2021017) and Capital’s Funds for Health Improvement and Research (CFH 2022-2Z-40917).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2024.1197651/full#supplementary-material

References

Beck, L. H., Bonegio, R. G., Lambeau, G., Beck, D. M., Powell, D. W., Cummins, T. D., et al. (2009). M-type phospholipase A2 receptor as target antigen in idiopathic membranous nephropathy. N. Engl. J. Med. 361 (1), 11–21. doi:10.1056/NEJMoa0810457

PubMed Abstract | CrossRef Full Text | Google Scholar

Bensalem, A., Cartron, G., Specks, U., Mulleman, D., Gyan, E., Cornec, D., et al. (2022). The influence of underlying disease on rituximab pharmacokinetics may be explained by target-mediated drug disposition. Clin. Pharmacokinet. 61 (3), 423–437. doi:10.1007/s40262-021-01081-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Bensalem, A., Mulleman, D., Paintaud, G., Azzopardi, N., Gouilleux-Gruart, V., Cornec, D., et al. (2020). Non-linear rituximab pharmacokinetics and complex relationship between rituximab concentrations and anti-neutrophil cytoplasmic antibodies (ANCA) in ANCA-associated vasculitis: the RAVE trial revisited. Clin. Pharmacokinet. 59 (4), 519–530. doi:10.1007/s40262-019-00826-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Blasco, H., Chatelut, E., de Bretagne, I. B., Congy-Jolivet, N., and Le Guellec, C. (2009). Pharmacokinetics of rituximab associated with CHOP chemotherapy in B-cell non-Hodgkin lymphoma. Fundam. Clin. Pharmacol. 23 (5), 601–608. doi:10.1111/j.1472-8206.2009.00714.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bondza, S., Foy, E., Brooks, J., Andersson, K., Robinson, J., Richalet, P., et al. (2017). Real-time characterization of antibody binding to receptors on living immune cells. Front. Immunol. 8, 455. doi:10.3389/fimmu.2017.00455

PubMed Abstract | CrossRef Full Text | Google Scholar

Boross, P., and Leusen, J. H. (2012). Mechanisms of action of CD20 antibodies. Am. J. Cancer Res. 2 (6), 676–690.

PubMed Abstract | Google Scholar

Boyer-Suavet, S., Andreani, M., Cremoni, M., Brglez, V., Benzaken, S., Bernard, G., et al. (2019). Rituximab bioavailability in primary membranous nephropathy. Nephrol. Dial. Transpl. 34 (8), 1423–1425. doi:10.1093/ndt/gfz041

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardone, K. E., Lodise, T. P., Patel, N., Hoy, C. D., Meola, S., Manley, H. J., et al. (2011). Pharmacokinetics and pharmacodynamics of intravenous daptomycin during continuous ambulatory peritoneal dialysis. Clin. J. Am. Soc. Nephrol. 6 (5), 1081–1088. doi:10.2215/CJN.08510910

PubMed Abstract | CrossRef Full Text | Google Scholar

Cravedi, P., Ruggenenti, P., Sghirlanzoni, M. C., and Remuzzi, G. (2007). Titrating rituximab to circulating B cells to optimize lymphocytolytic therapy in idiopathic membranous nephropathy. Clin. J. Am. Soc. Nephrol. 2 (5), 932–937. doi:10.2215/CJN.01180307

PubMed Abstract | CrossRef Full Text | Google Scholar

Dahan, K., Debiec, H., Plaisier, E., Cachanado, M., Rousseau, A., Wakselman, L., et al. (2017). Rituximab for severe membranous nephropathy: a 6-month trial with extended follow-up. J. Am. Soc. Nephrol. 28 (1), 348–358. doi:10.1681/ASN.2016040449

PubMed Abstract | CrossRef Full Text | Google Scholar

Del Vecchio, L., Allinovi, M., Rocco, P., and Brando, B. (2021). Rituximab therapy for adults with nephrotic syndromes: standard schedules or B cell-targeted therapy? J. Clin. Med. 10 (24), 5847. doi:10.3390/jcm10245847

PubMed Abstract | CrossRef Full Text | Google Scholar

Fenoglio, R., Baldovino, S., Sciascia, S., De Simone, E., Del Vecchio, G., Ferro, M., et al. (2021). Efficacy of low or standard rituximab-based protocols and comparison to Ponticelli's regimen in membranous nephropathy. J. Nephrol. 34 (2), 565–571. doi:10.1007/s40620-020-00781-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Fervenza, F. C., Appel, G. B., Barbour, S. J., Rovin, B. H., Lafayette, R. A., Aslam, N., et al. (2019). Rituximab or cyclosporine in the treatment of membranous nephropathy. N. Engl. J. Med. 381 (1), 36–46. doi:10.1056/NEJMoa1814427

PubMed Abstract | CrossRef Full Text | Google Scholar

Fogueri, U., Cheungapasitporn, W., Bourne, D., Fervenza, F. C., and Joy, M. S. (2019). Rituximab exhibits altered pharmacokinetics in patients with membranous nephropathy. Ann. Pharmacother. 53 (4), 357–363. doi:10.1177/1060028018803587

PubMed Abstract | CrossRef Full Text | Google Scholar

Gibiansky, L., Gibiansky, E., Kakkar, T., and Ma, P. (2008). Approximations of the target-mediated drug disposition model and identifiability of model parameters. J. Pharmacokinet. Pharmacodyn. 35 (5), 573–591. doi:10.1007/s10928-008-9102-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Ginaldi, L., De Martinis, M., Matutes, E., Farahat, N., Morilla, R., and Catovsky, D. (1998). Levels of expression of CD19 and CD20 in chronic B cell leukaemias. J. Clin. Pathol. 51 (5), 364–369. doi:10.1136/jcp.51.5.364

PubMed Abstract | CrossRef Full Text | Google Scholar

Glassman, P. M., and Balthasar, J. P. (2017). Physiologically-based modeling to predict the clinical behavior of monoclonal antibodies directed against lymphocyte antigens. MAbs 9 (2), 297–306. doi:10.1080/19420862.2016.1261775

PubMed Abstract | CrossRef Full Text | Google Scholar

Iwabuchi, Y., Moriyama, T., Itabashi, M., Takei, T., and Nitta, K. (2018). Rituximab as a therapeutic option for steroid-sensitive minimal change nephrotic syndrome in adults. Contrib. Nephrol. 195, 12–19. doi:10.1159/000486930

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Levi, M., Charoin, J.-E., Frey, N., Kheoh, T., Ren, S., et al. (2007). Rituximab exhibits a long half-life based on a population pharmacokinetic analysis in non-hodgkin's lymphoma (NHL) patients. Blood 110 (11), 2371. doi:10.1182/blood.v110.11.2371.2371

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Zhi, J., Wenger, M., Valente, N., Dmoszynska, A., Robak, T., et al. (2012). Population pharmacokinetics of rituximab in patients with chronic lymphocytic leukemia. J. Clin. Pharmacol. 52 (12), 1918–1926. doi:10.1177/0091270011430506

PubMed Abstract | CrossRef Full Text | Google Scholar

Macallan, D. C., Wallace, D. L., Zhang, Y., Ghattas, H., Asquith, B., de Lara, C., et al. (2005). B-cell kinetics in humans: rapid turnover of peripheral blood memory cells. Blood 105 (9), 3633–3640. doi:10.1182/blood-2004-09-3740

PubMed Abstract | CrossRef Full Text | Google Scholar

Moroni, G., Depetri, F., Del Vecchio, L., Gallelli, B., Raffiotta, F., Giglio, E., et al. (2017). Low-dose rituximab is poorly effective in patients with primary membranous nephropathy. Nephrol. Dial. Transpl. 32 (10), 1691–1696. doi:10.1093/ndt/gfw251

PubMed Abstract | CrossRef Full Text | Google Scholar

Ng, C. M., Bruno, R., Combs, D., and Davies, B. (2005). Population pharmacokinetics of rituximab (anti-CD20 monoclonal antibody) in rheumatoid arthritis patients during a phase II clinical trial. J. Clin. Pharmacol. 45 (7), 792–801. doi:10.1177/0091270005277075

PubMed Abstract | CrossRef Full Text | Google Scholar

Regazzi, M. B., Iacona, I., Avanzini, M. A., Arcaini, L., Merlini, G., Perfetti, V., et al. (2005). Pharmacokinetic behavior of rituximab: a study of different schedules of administration for heterogeneous clinical settings. Ther. Drug Monit. 27 (6), 785–792. doi:10.1097/01.ftd.0000184162.60197.c1

PubMed Abstract | CrossRef Full Text | Google Scholar

Remuzzi, G., Chiurchiu, C., Abbate, M., Brusegan, V., Bontempelli, M., and Ruggenenti, P. (2002). Rituximab for idiopathic membranous nephropathy. Lancet 360 (9337), 923–924. doi:10.1016/S0140-6736(02)11042-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Ronco, P., Beck, L., Debiec, H., Fervenza, F. C., Hou, F. F., Jha, V., et al. (2021). Membranous nephropathy. Nat. Rev. Dis. Prim. 7 (1), 69. doi:10.1038/s41572-021-00303-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Rovin, B. H., Adler, S. G., Barratt, J., Bridoux, F., Burdge, K. A., Chan, T. M., et al. (2021). Executive summary of the KDIGO 2021 guideline for the management of glomerular diseases. Kidney Int. 100 (4), 753–779. doi:10.1016/j.kint.2021.05.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruggenenti, P., Debiec, H., Ruggiero, B., Chianca, A., Pelle, T., Gaspari, F., et al. (2015). Anti-phospholipase A2 receptor antibody titer predicts post-rituximab outcome of membranous nephropathy. J. Am. Soc. Nephrol. 26 (10), 2545–2558. doi:10.1681/ASN.2014070640

PubMed Abstract | CrossRef Full Text | Google Scholar

Seitz-Polski, B., Dahan, K., Debiec, H., Rousseau, A., Andreani, M., Zaghrini, C., et al. (2019). High-dose rituximab and early remission in pla2r1-related membranous nephropathy. Clin. J. Am. Soc. Nephrol. 14 (8), 1173–1182. doi:10.2215/CJN.11791018

PubMed Abstract | CrossRef Full Text | Google Scholar

Teisseyre, M., Cremoni, M., Boyer-Suavet, S., Crepin, T., Benzaken, S., Zorzi, K., et al. (2021). Rituximab immunomonitoring predicts remission in membranous nephropathy. Front. Immunol. 12, 738788. doi:10.3389/fimmu.2021.738788

PubMed Abstract | CrossRef Full Text | Google Scholar

Teisseyre, M., Cremoni, M., Boyer-Suavet, S., Ruetsch, C., Graca, D., Esnault, V. L. M., et al. (2022). Advances in the management of primary membranous nephropathy and rituximab-refractory membranous nephropathy. Front. Immunol. 13, 859419. doi:10.3389/fimmu.2022.859419

PubMed Abstract | CrossRef Full Text | Google Scholar

Ternant, D., Monjanel, H., Venel, Y., Prunier-Aesch, C., Arbion, F., Colombat, P., et al. (2019). Nonlinear pharmacokinetics of rituximab in non-Hodgkin lymphomas: a pilot study. Br. J. Clin. Pharmacol. 85 (9), 2002–2010. doi:10.1111/bcp.13991

PubMed Abstract | CrossRef Full Text | Google Scholar

Tout, M., Gagez, A. L., Lepretre, S., Gouilleux-Gruart, V., Azzopardi, N., Delmer, A., et al. (2017). Influence of fcgr3a-158V/F genotype and baseline CD20 antigen count on target-mediated elimination of rituximab in patients with chronic lymphocytic leukemia: a study of filo group. Clin. Pharmacokinet. 56 (6), 635–647. doi:10.1007/s40262-016-0470-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Deng, Z., Wang, Y., Bao, W., Zhou, S., Cui, Z., et al. (2023). Monthly mini-dose rituximab for primary anti-PLA2R-positive membranous nephropathy: a personalized approach. BMC Nephrol. 24 (1), 146. doi:10.1186/s12882-023-03206-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Du, W., Zhang, X., and Li, P. (2020). The influence of different disease states on rituximab pharmacokinetics. Curr. Drug Metab. 21, 938–946. doi:10.2174/1389200221666200719004035

PubMed Abstract | CrossRef Full Text | Google Scholar

Zachwieja, J., Silska, M., Ostalska-Nowicka, D., Soltysiak, J., Lipkowska, K., Blumczynski, A., et al. (2012). Efficacy and safety of rituximab treatment in children with primary glomerulonephritis. J. Nephrol. 25 (6), 1060–1066. doi:10.5301/jn.5000096

PubMed Abstract | CrossRef Full Text | Google Scholar