Skip to main content

PK/PD analysis of biapenem in patients undergoing continuous hemodiafiltration



Continuous hemodiafiltration (CHDF) is used as renal replacement therapy for critically ill patients with renal failure, and to treat hypercytokinemia. Since CHDF also clears therapeutic agents, drug pharmacokinetics (PK) should be dependent upon CHDF conditions. Although the antibiotic biapenem (BIPM) is used in patients undergoing CHDF, the optimal therapeutic regimen in such patients has not been fully clarified. In this study, we investigated the PK of BIPM in patients with various levels of renal function undergoing CHDF with polysulfone (PS) membrane, and used PK models to identify the optimal administration regimen.


BIPM (300 mg) was administered by infusion in patients undergoing CHDF (n = 7). Blood and filtrate-dialysate were collected for compartment and non-compartment analysis.


The sieving coefficient of PS membrane was 1.00 ± 0.06 (mean ± S.D., n = 7), and CHDF clearance of BIPM was found to be the sum of the dialysate flow rate (QD) and filtrate flow rate (QF). Non-CHDF clearance showed inter-individual variability (4.82 ± 2.48 L/h), depending on residual renal function and non-renal clearance. Based on the average PK parameters obtained with a compartmental model, maximal kill end point (over 40 % T > MIC4 μg/mL) was achieved with regimens of 300 mg every 6 h, 300 mg every 8 h, and 600 mg every 12 h. Monte Carlo simulation indicated that 300 mg infusion for 1 h every 6 h was optimal, and the probability of target attainment at MIC2 μg/mL was 90.2 %.


Our results establish the optimal regimen of BIPM in patients with various levels of renal function undergoing CHDF with a PS membrane.


Continuous hemodiafiltration (CHDF) is a blood purification method that provides alternative functionality for patients with kidney failure. It removes a wide range of compounds, ranging from low-molecular-weight products such as urea nitrogen and creatinine to medium-molecular-weight compounds such as β-2 microglobulin and cytokines [1, 2]. Therefore, CHDF is widely used in intensive care units to treat patients with unstable hemodynamic pathologies, including acute pancreatitis, fulminant hepatitis, and acute renal failure, as well as to remove cytokines in sepsis patients, even if they have normal renal function [36]. Since the pharmacokinetics (PK) of medications may be affected by elimination through dialysis and filtration during CHDF therapy, it is necessary to select an appropriate administration regimen (timing, dose amount, and dose rate) depending upon the conditions of CHDF [7, 8]. It has also been reported that drug clearance varies depending on the material of the dialysis membrane [9], and cytokine clearance is also dependent on the dialysis membrane material [2].

Biapenem (BIPM) is a carbapenem antibacterial agent with broad-spectrum, potent activity against Gram-positive, Gram-negative, and anaerobic bacteria [10]. It is not cleaved by dehydropeptidase-1 (DHP-1), unlike other penem antibiotics such as panipenem and imipenem, and can be used as a single agent without the need for formulation of a nephrotoxicity-reducing agent or DHP-1 inhibitor [1113]. Carbapenems are classified as time-dependent antibiotics, and the clinical outcome is closely correlated with the duration for which the drug concentration remains at or above the minimum inhibitory concentration (MIC); the time above the MIC (T > MIC) needs to be over 20 and 40 % of the dosing interval to achieve stasis and maximal kill as end points, respectively [14, 15]. For carbapenem antibiotics, sensitive MIC and intermediate MIC are defined as ≥ 4 μg/mL and ≥ 8 μg/mL, respectively [16]. To obtain optimum effects of BIPM, we need to predict T > MIC from the PK parameters of BIPM. BIPM is mainly excreted from the kidneys, and clearance of BIPM in patients with renal failure is decreased to 2.62 L/h, i.e., approximately 20 % of that of healthy adults [17]. CHDF also alters the clearance of BIPM, and Ikawa et al. have reported on PK modeling and dosage adaptation of BIPM during CHDF with a polymethyl methacrylate (PMMA) membrane in patients with renal failure [18]. The PK of single BIPM administration was calculated using multi-compartment models, and it was found that the dose amount and dosing interval were important factors determining the value of % of T > MIC. However, their models were applicable only to patients with renal failure. Since CHDF is also used for patients with various levels of renal function, models that take account of this are required. Furthermore, the effects of other types of filter membrane in CHDF should also be considered.

In the present study, we developed non-compartmental and compartmental models for repeated BIPM treatment in patients with various levels of renal function undergoing CHDF with a polysulfone (PS) membrane, and examined the optimum treatment regimen by means of Monte Carlo simulation.



This study was approved by the Clinical Research Ethics Committee of Kanazawa University Hospital (2011–052). Adult patients older than 20 years who were receiving CHDF for acute renal failure or removal of cytokines, and who were prescribed BIPM to treat infection were eligible. Prior to the start of the investigation, informed consent was obtained in writing from each participant or his/her relatives. We excluded patients taking sodium valproate, or with a history of hypersensitivity to carbapenems or fourth-generation cephems, or with contraindications for carbapenems or fourth-generation cephems, or with a history of epilepsy or central nervous system damage. Seven patients (5 males, 2 females) were enrolled. Demographic and medical characteristics of each patient are summarized in Table 1. The glomerular filtration rate (GFR) of each patient was evaluated before CHDF application.

Table 1 Patient characteristics


In the present study, a PS hemofilter with a membrane surface area of 1.3 m2 (EXCELFLO® AEF-13, Asahi Kasei Medical Co., Japan) was used for hemopurification. CHDF was performed basically with a blood flow rate (QB) of 80 mL/min, dialysate flow rate (QD) of 500 mL/h, substitution flow rate (QS) of 500 mL/h, and filtration flow rate (QF) of 1000 mL/h, and was regulated appropriately based on each patient’s status. Sublood-BS® (Fuso Pharmacy, Inc., Japan) was used as a dialysate and also served as a substitution fluid for post-dilutional infusion. During CHDF, nafamostat mesylate was serially administered at a rate of approximately 20 mg/h to avoid coagulation within the circuit.

BIPM administration and collection of blood and filtrate-dialysate samples

Three hundred milligrams of BIPM (Omegacin® 0.3 g, Meiji Seika Pharma Co., Japan) was administered by intravenous drip infusion for 0.5 or 1 h at 6, 8, or 12 h intervals (Table 1). Blood samples were taken from the blood access port in the extracorporeal circuit proximal to the filter before dosing, at 0, 0.5, 1, 2, and 4 h after infusion, and just before the next infusion. Filtrate-dialysate (FD) samples were simultaneously collected from the filtrate tube to determine the sieving coefficient (SC) and CHDF clearance. After collection, each blood sample was centrifuged immediately. The plasma and supernatants of FD samples were immediately frozen at −80 °C and stored until assay.

Method of quantitative analysis of BIPM

Quantitative analysis of BIPM in plasma and FD was performed using high-performance liquid chromatography (HPLC) [18]. Briefly, plasma or FD sample (200 μL) was mixed with 200 μL of 1 M 3-morpholinopropane-1-sulfonic acid buffer. An aliquot (300 μL) was transferred to an ultrafiltration device (Nanosep® Centrifugal Devices 10 K, Pall Life Sciences, USA), and centrifuged at 15,000 g for 15 min at 15 °C. An aliquot of the ultrafiltered solution (20 μL) was injected into an HPLC system (Prominence, Shimadzu Co., Japan) equipped with an Octadecyl silica column (Shim-pack CLC-ODS 5 μm (150 x 6 i.d.), Shimadzu Co., Japan) and a UV detector (232 nm). The mobile phase of 1.5 % acetonitrile and 98.5 % 0.1 M sodium acetate buffer (pH 4.6) was delivered isocratically at a flow rate of 1.0 mL/min. The auto-sampler temperature was set at 4 °C, and the column temperature at 40 °C. The calibration curve was linear from 0.1 to 50 μg/mL, and the lower limit of detection was 0.1 μg/mL.

Calculation of PK parameters of BIPM

Non-compartmental and compartmental analyses were performed using PK analysis software: Numeric Analysis Program for Pharmacokinetics (Napp) ver. 2.31 (Department of Pharmacy, the University of Tokyo Hospital). The area under the BIPM concentration-time curve (AUC) was calculated based on the logarithmic trapezoidal rule. The SC was determined as AUCFD/AUCplasma. The AUC after first administration was calculated from 0 h to infinity, and subsequent AUCs were estimated from 0 h to the start of the next infusion interval, since the half-life of BIPM in healthy persons is 1 h [19]. The differential equations for mass balances were as follows (Fig. 1):

Fig. 1

Multicompartment model for BIPM during CHDF. X 1, X 2 and X 3 are the amounts of drug in the central, peripheral and hemofilter compartments, respectively; V 1, V 2 and V 3 are the volumes of the central compartment, peripheral compartment and hemofilter cartridge, respectively; R inf, drug infusion rate; k 10, 30, elimination rate constants; k 12, 13, 21, rate constants connecting the compartments

$$ \frac{d{X}_1}{dt}={R}_{inf} - \left({k}_{10}+{k}_{12}+{k}_{13}\right)\cdot {X}_1 + {k}_{21}\cdot {X}_2 $$
$$ \frac{d{X}_2}{dt}={k}_{12}\cdot {X}_1-{k}_{21}\cdot {X}_2 $$
$$ \frac{d{X}_3}{dt}={k}_{13}\cdot {X}_1-{k}_{30}\cdot {X}_3 $$

where X 1, X 2 and X 3 are the amounts of drug (mg) in the central, peripheral and hemofilter compartment, respectively; R inf is the drug infusion rate (mg/hr); k 10 and k 30 are the elimination rate constant (hr−1) from the central compartment and the hemofilter compartment, respectively; k 12 and k 21 are the intercompartmental transfer rate constants (hr−1); k 13 is transfer rate constant (hr−1) from the central compartment to the hemofilter compartment.

These equations were used to obtain the following formulas for the drug concentrations (μg/mL) in the central compartment (C1) and FD (C3) at time t (hr) during and after administration in the steady state:

  1. (i)

    0 t Tinf

$$ {C}_1=\frac{R_{inf}}{V_1}\left\{\frac{\left({k}_{21}-\alpha \right)}{\alpha \left(\beta -\alpha \right)}\left({e}^{\alpha \cdot {T}_{inf}}-1\right)\left(\frac{e^{-\alpha \cdot \tau }}{1-{e}^{-\alpha \cdot \tau }}\right){e}^{-\alpha \cdotp t}+\frac{\left({k}_{21}-\beta \right)}{\beta \left(\alpha -\beta \right)}\left({e}^{\beta \cdot {T}_{inf}}-1\right)\left(\frac{e^{-\beta \cdot \tau }}{1-{e}^{-\beta \cdot \tau }}\right){e}^{-\beta \cdotp t}+\frac{\left({k}_{21}-\alpha \right)}{\alpha \left(\beta -\alpha \right)}\left(1-{e}^{-\alpha \cdotp t}\right)+\frac{\left({k}_{21}-\beta \right)}{\beta \left(\alpha -\beta \right)}\left(1-{e}^{-\beta \cdotp t}\right)\right\} $$
$$ {C}_3=\frac{R_{inf}\cdotp {k}_{13}}{V_3}\left\{\frac{\left({k}_{21}-\alpha \right)}{\alpha \left(\beta -\alpha \right)\left({k}_{30}-\alpha \right)}\left({e}^{\alpha \cdotp {T}_{inf}}-1\right)\left(\frac{e^{-\alpha \cdotp \tau }}{1-{e}^{-\alpha \cdotp \tau }}\right){e}^{-\alpha \cdotp t}+\frac{\left({k}_{21}-\beta \right)}{\beta \left(\alpha -\beta \right)\left({k}_{30}-\beta \right)}\left({e}^{\beta \cdotp {T}_{inf}}-1\right)\left(\frac{e^{-\beta \cdotp \tau }}{1-{e}^{-\beta \cdotp \tau }}\right){e}^{-\beta \cdotp t}+\frac{\left({k}_{21}-{k}_{30}\right)}{k_{30}\left(\alpha -{k}_{30}\right)\left(\beta -{k}_{30}\right)}\left({e}^{k_{30}\cdotp {T}_{inf}}-1\right)\left(\frac{e^{-{k}_{30}\cdotp \tau }}{1-{e}^{-{k}_{30}\cdotp \tau }}\right){e}^{-{k}_{30}\cdotp t}+\frac{\left({k}_{21}-\alpha \right)}{\alpha \left(\beta -\alpha \right)\left({k}_{30}-\alpha \right)}\left(1-{e}^{-\alpha \cdotp t}\right)+\frac{\left({k}_{21}-\beta \right)}{\beta \left(\alpha -\beta \right)\left({k}_{30}-\beta \right)}\left(1-{e}^{-\beta \cdotp t}\right)+\frac{\left({k}_{21}-{k}_{30}\right)}{k_{30}\left(\alpha -{k}_{30}\right)\left(\beta -{k}_{30}\right)}\left(1-{e}^{-{k}_{30}\cdotp t}\right)\right\} $$
  1. (ii)

    t > Tinf

$$ {C}_1=\frac{R_{inf}}{V_1}\left\{\frac{\left({k}_{21}-\alpha \right)}{\alpha \left(\beta -\alpha \right)}\left({e}^{\alpha \cdotp {T}_{inf}}-1\right)\left(\frac{1}{1-{e}^{-\alpha \cdotp \tau }}\right){e}^{-\alpha \cdotp t}+\frac{\left({k}_{21}-\beta \right)}{\beta \left(\alpha -\beta \right)}\left({e}^{\beta \cdotp {T}_{inf}}-1\right)\left(\frac{1}{1-{e}^{-\beta \cdotp \tau }}\right){e}^{-\beta \cdotp t}\right\} $$
$$ {C}_3=\frac{R_{inf}\cdotp {k}_{13}}{V_3}\left\{\frac{\left({k}_{21}-\alpha \right)}{\alpha \left(\beta -\alpha \right)\left({k}_{30}-\alpha \right)}\left({e}^{\alpha \cdotp {T}_{inf}}-1\right)\left(\frac{1}{1-{e}^{-\alpha \cdotp \tau }}\right){e}^{-\alpha \cdotp t}+\frac{\left({k}_{21}-\beta \right)}{\beta \left(\alpha -\beta \right)\left({k}_{30}-\beta \right)}\left({e}^{\beta \cdotp {T}_{inf}}-1\right)\left(\frac{1}{1-{e}^{-\beta \cdotp \tau }}\right){e}^{-\beta \cdotp t}+\frac{\left({k}_{21}-{k}_{30}\right)}{k_{30}\left(\alpha -{k}_{30}\right)\left(\beta -{k}_{30}\right)}\left({e}^{k_{30}\cdotp {T}_{inf}}-1\right)\left(\frac{1}{1-{e}^{-{k}_{30}\cdotp \tau }}\right){e}^{-{k}_{30}\cdotp t}\right\} $$

Where Tinf is the drug infusion time (hr); τ is the inter-dose interval (hr); V 1 is the volume of distribution of the central compartment; V 3 is the blood volume of the hemofilter cartridge; and α and β are macro rate constants (hr−1) expressed as α + β = k 12 + k 21 + k 10 + k 13, αβ = k 21 · k 10 + k 21 · k 13.

The rate constants k 13 and k 30 can be expressed as:

$$ {k}_{13}=\frac{\left({Q}_F+{Q}_D\right)\cdotp SC}{V_1} $$
$$ {k}_{30}=\frac{\left({Q}_F+{Q}_D\right)}{V_3} $$

Where QF is the filtrate flow rate (L/hr); QD is the dialysate flow rate (L/hr).

The clearance in CHDF (CLCHDF) was calculated as (QF + QD) · SC and the clearance by non-CHDF routes (CLnon-CHDF) was estimated as k 10 · V 1.

We modified the formula of Ikawa et al. [18] to include data obtained not only after the first administration of BIPM, but also those obtained during repeated administration of BIPM. The concentration-time data of BIPM in plasma and FD were concurrently fitted to the multi-compartment model described above.

Glomerular filtration rate (GFR) (mL/min/body) was calculated as follows:

$$ GFR=\frac{eGFR\cdot BSA}{1.73} $$
$$ male\ :\ eGFR=194\cdot Sc{r}^{-1.094}\cdot Ag{e}^{-0.287} $$
$$ BSA=\frac{105.29\cdot H{t}^{0.619}\cdot B{W}^{0.460}}{10,000} $$
$$ female\ : eGFR=194\cdot Sc{r}^{-1.094}\cdot Ag{e}^{-0.287}\cdot 0.739 $$
$$ BSA=\frac{82.84\cdot H{t}^{0.689}\cdot B{W}^{0.437}}{10,000} $$

Where eGFR is estimated glomerular filtration rate (mL/min/1.73 m2); BSA is body surface area (m2) [20]; Scr is serum creatinine concentration (mg/dL); Age is given in years; Ht is height (cm); BW is body weight (kg). eGFR was determined just before CHDF application.

Calculation of the PK/PD breakpoint by Monte Carlo simulation

Monte Carlo simulation was performed with the normal random number generation function of Microsoft® Excel® 2007 (Microsoft, USA). We generated a population parameter set of 10,000 cases using a mean and variance of k 21, V 1, α and β obtained by the standard two-stage method (compartmental analysis), thereby generating the BIPM plasma concentration transition of 10,000 cases. After setting the dosing interval and administration time of BIPM in the above formula (C1), the exposure time for which the BIPM plasma concentration remained at the MIC was finally calculated as the cumulative percentage over a 24 h period in the BIPM plasma concentration transition of the 10,000 cases [2123]. Then, we calculated the number of cases that showed T > MIC of 40 % or more at each PK/PD target value in each regimen. Maximum MIC of more than 80 % was set as the PK/PD breakpoint probability of target attainment (PTA) in each regimen.


In the present study, no adverse events or laboratory abnormalities were noted that were definitely attributable to BIPM.

SC of the PS membrane was calculated to be 1.00 ± 0.06 by non-compartmental analysis (Table 2). The PK parameters estimated by compartmental analysis were as follows: V 1, 13.46 ± 5.29 L; k 12, 0.75 ± 0.71 h−1; k 21, 0.95 ± 0.27 h−1; k 10, 0.39 ± 0.22 h−1; and k 13, 0.13 ± 0.04 h−1 (Table 3). Figure 2 shows typical fittings between the plasma and FD concentration time curves obtained by compartmental analysis and the real values in Pt. A and Pt. B. Time courses of BIPM concentration in plasma and FD obtained by compartmental analysis in other patients also closely fitted the observed values. CLCHDF and CLnon-CHDF were estimated to be 1.53 ± 0.10 L/h and 4.86 ± 2.50 L/h, respectively. The inter-individual variation in CLnon-CHDF was large. The sum of CLCDHFand CLnon-CHDF, calculated as (QF + QD) · SC and as k 10 · V 1, respectively [x], and CLtot obtained by non-compartmental analysis of the plasma BIPM concentration [y] showed a good correlation (y = 1.01 x − 0.02, r 2 = 1.00), supporting the validity of the model. As shown in Fig. 3, the correlation between GFR obtained just before CHDF and CLnon-CHDF was also high (r 2 = 0.97), and the slope and y-intercept were 1.86 and 1.02, respectively. In this study, the number of patients who showed T > MIC4 μg/mL values of 40 % or more was 6, and one was in the range of 20 to 40 %. Moreover, the number of patients with T > MIC8 μg/mL values of more than 40 % was one, and 6 were in the range of 20 to 40 %.

Table 2 Pharmacokinetic parameters calculated by non-compartmental analysis
Table 3 Pharmacokinetic parameters of biapenem after intravenous administration of biapenem (300 mg) during CHDF
Fig. 2

Typical fittings between the plasma and FD concentration simulation curves and the observed values. Simulation curves of plasma (open circles, solid line) and FD (closed triangle, dotted line) concentration were obtained by compartmental analysis

Fig. 3

Correlation between individual GFR and CLnon-CHDF. The solid line indicates the linear regression line

Figure 4 shows % of T > MIC values in plasma of patients undergoing CHDF based on the average parameters of the multi-compartment models on various regimens, including different values of infusion time, MIC values, sum of QF and QD, dosage and administration. As QF + QD increased, % of T > MIC slightly decreased. In addition, % of T > MIC for 1 h infusion was higher than that for 0.5 h infusion under the same conditions of QF + QD and administration dose. The most effective regimen was 300 mg every 6 h given by 1 h infusion, with a QF + QD value of 1.5 L/h. This method could achieve % of T > MIC of 82.8 % at MIC = 4 μg/mL, and 35.1 % at MIC = 8 μg/mL.

Fig. 4

Values of T > MIC of 4 μg/mL and 8 μg/mL in different regimens. Panel a represents the case of MIC = 4 μg/mL and infusion time 1 h; Panel b represents the case of MIC = 4 μg/mL and infusion time 0.5 h; Panel c represents the case of MIC = 8 μg/mL and infusion time 1 h; Panel d represents the case of MIC = 8 μg/mL and infusion time 0.5 h. The solid line represents 40 % T > MIC. The dotted line represents 20 % T > MIC

Based on Monte Carlo simulation of the case of 0.5 h infusion, the PK/PD breakpoint was 0.5 μg/mL in the regimen of dosage 300 mg every 12 h, and was 1.0 μg/mL in the regimens of dosage 300 mg every 8 h, dosage 300 mg every 6 h, and dosage 600 mg every 12 h (Fig. 5). In the case of 1 h infusion, the PK/PD breakpoint was 1.0 μg/mL in the regimen of dosage 300 mg every 12 h, and 2.0 μg/mL in the regimens of dosage 300 mg every 8 h, dosage 300 mg every 6 h, and dosage 600 mg every 12 h. The PTAs were greater for 1 h infusion than for 0.5 h infusion in all cases. Among them, PTA of the regimen of dosage 300 mg every 6 h showed the highest value (90.2 %). The Japan Society of Chemotherapy defines the clinical breakpoint as more than 80 %, though more than 90 % of the PTA has become a reference value of the PK/PD breakpoint outside of Japan.

Fig. 5

PTA of 40 % T > MIC in plasma at specific MICs for different BIPM regimens. The dotted line represents 80 % PTA. Panel a: infusion time of 1 h; Panel b: infusion time of 0.5 h


In this study, we investigated the PK of BIPM using compartmental and non-compartmental analysis in patients undergoing CHDF who had various levels of renal function, and examined the suitability of various administration regimens in terms of PK/PD breakpoint for various bacteria, by means of Monte Carlo simulation.

Ikawa et al. analyzed time courses of BIPM concentration in plasma and FD after single administration in renal failure patients undergoing CHDF [18]. Here, we aimed to obtain a model to analyze the results of repeated administration of BIPM in patients undergoing CHDF who retained various levels of renal function. Time courses of BIPM concentration in plasma and FD in patients undergoing CHDF closely fitted the observed values, suggesting that the constructed model formula successfully represents the results of repeated administration of BIPM in patients undergoing CHDF. In the present study, SC (indicating penetration from blood to dialyzer) was almost 1.0 for the PS membrane material, which is almost the same as the value for PMMA membrane [18, 24]. Thus, there appears to be essentially no difference between PS and PMMA membranes in regard to drug penetration [18]. These results seem consistent with the characteristics of BIPM, such as low molecular weight and low protein binding rate.

We found a strong correlation between the sum of CLnon-CHDF and CLCHDF obtained by compartmental analysis and total clearance obtained by non-compartmental analysis (r 2 = 1.00), indicating that the compartment model formula was appropriate. Ikawa et al. reported that CLCHDF (1.29 ± 0.08 L/h) was almost the same as the sum of QF and QD (1.4 L/h) in CHDF with a PMMA membrane [18]. Suyama et al. also reported that CLCHDF (1.28 ± 0.14 L/h) was almost the same as the sum of QF and QD (1.4 L/h) in CHDF with a PMMA membrane [24]. In this study, CLCHDF (1.5 ± 0.1 L/h) estimated by compartmental analysis was similar to the sum of QF and QD (1.5–1.7 L/h) in CHDF with a PS membrane. Considering the SC estimated by non-compartmental analysis, CLCHDF of BIPM would be determined by the sum of dialysate flow rate and filtration flow rate. On the other hand, CLnon-CHDF represents the sum of non-renal clearance and residual renal clearance. Metabolism of BIPM in kidney and other tissues would contribute to CLnon-CHDF. Nakajima et al. detected two DHP-1 metabolites of BIPM, L-cysteine and L-cystine, in urine, and reported that excretion of urinary metabolites accounted for approximately 15 % of total clearance; these metabolites were not detected in plasma, suggesting that metabolism of BIPM in kidney contributes to residual renal clearance [19]. DHP-1 is also expressed in the ascending colon and ileum as well as the kidney, and the activity of DHP-1 in ileum has been reported to be twice that in kidney, suggesting that non-renal DHP-1 function in ileum would also contribute to non-renal clearance [25]. The relationship between GFR calculated before CHDF application and CLnon-CHDF was given by y = 1.86 x + 1.02 (r 2 = 0.97) (Fig. 3), and the y-intercept (1.02 L/h) represents non-renal clearance. Non-renal clearance in healthy adults can be estimated as approximately 3.23 L/h, using reported values of total clearance of BIPM in healthy adults of 12.9 ± 1.2 L/h, and urinary excretion rates of unchanged and changed BIPM of approximately 60 and 15 %, respectively [19]. Nagashima et al. reported that non-renal clearance was 2.60 ± 1.55 L/h in patients with renal failure during dialysis with a PS membrane, while CLtot in patients was 2.62 ± 0.60 L/h without HD, so that CLtot was the same as non-renal clearance in patients with renal failure [17]. Non-renal clearance in the present study was a half to one-third of those values. After extracting patients with sepsis from Fig. 3, the y-intercept took a negative value (data not shown), indicating that non-renal clearance would be almost zero in patients with sepsis. In sepsis, hepatic clearance and renal clearance would be diminished owing to the reduced function of systemic organs and impaired blood flow [26], leading to loss of DHP-1 activity in the kidney and other organs. Thus, we consider that the reason why average non-renal clearance in this study was lower than in previous studies is the low values in patients with sepsis.

We investigated the variation in % of T > MIC based on the average values of the PK parameters obtained by compartmental analysis for different regimens, including dose, dosing interval, infusion time, and CHDF conditions applied in our hospital. The % of T > MIC4 μg/mL was more than 40 % in all regimens of 900 mg or more daily dose, regardless of the CHDF conditions and infusion time. Ikawa et al. also reported that although the regimen of 300 mg every 12 h failed to achieve T > MIC4 μg/mL of more than 40 %, the regimen of 600 mg every 12 h did do so [18]. These results indicate that the minimum dosage regimen required to achieve T > MIC4 μg/mL of more than 40 % is 300 mg every 8 h (total amount: 900 mg). On the other hand, no regimen gave % of T > MIC8 μg/mL of more than 40 %, suggesting that these regimens could not achieve the maximal kill end point. However, since the % of T > MIC values obtained using average PK parameters do not reflect variations of CLnon-CHDF associated with partial kidney function, these values may not be adequate as a clinical indicator.

Monte Carlo simulation is a computer modeling process that incorporates variability in pharmacokinetic parameters and the natural MIC distribution within a bacterial population. It can be used to develop interpretive susceptibility criteria based on PK/PD breakpoints [27]. In this study, Monte Carlo simulations with 10,000 cases were performed to examine PK/PD breakpoint using the mean and variance of the PK parameters. As regards infusion time, the PK/PD breakpoint obtained for all regimens with 1 h infusion was higher than that for 0.5 h infusion, suggesting that the 1 h infusion provides a better outcome. In this study, the highest PK/PD breakpoint (over 80 % PTA) was 2 μg/mL with 300 mg every 8 h, 300 mg every 6 h, and 600 mg every 12 h, indicating that these regimens provide a sufficient bactericidal effect in the case of bacteria with MIC2 μg/mL, but not MIC4 μg/mL. In order to obtain an antimicrobial effect of BIPM towards high MIC bacteria (more than 4 μg/mL), such as Pseudomonas aeruginosa, higher dose administration and an appropriate regimen would be needed, although the maximum permitted dose of BIPM in Japan has been set at 1,200 mg/day. Among the regimens examined in this study, 300 mg every 6 h with infusion for 1 h was the optimal regimen, and PTA at MIC2 μg/mL was 90.2 %. Although maximum MIC more than 80 % is adopted as the PK/PD breakpoint in Japan, most other countries define it as more than 90 % [27], so the regimen of 300 mg every 6 h with infusion for 1 h also meets the international standard for MIC2 μg/mL bacteria. Although the PK/PD breakpoint is expected to become a decision criterion for individualized and optimized antibacterial therapy, it is focused primarily on the antibacterial agent. In the future, it will be important to take bacterial character into account as well, for example, by means of antimicrobial susceptibility testing and identification of the causative bacteria of infectious diseases.

The present study has several limitations. Firstly, we performed the Monte Carlo simulation using the mean and variance of the PK parameters obtained with the standard two-stage method because of the small number of cases (total, 7 cases), and this might have resulted in overestimation of the inter-individual variation compared to population pharmacokinetics analysis, such as nonlinear mixed effects modelling (NONMEM®). Further, we did not include a parameter of renal function (e.g. GFR or creatinine clearance) in the Monte Carlo simulation, even though BIPM clearance is correlated with renal function. Therefore, as a next step, it would be desirable to identify the PK/PD breakpoint more precisely in patients undergoing CHDF by means of Monte Carlo simulation of BIPM using population pharmacokinetic parameters including renal function (e.g. GFR and creatinine clearance), based on larger numbers of patients with various levels of renal function. In addition, from the viewpoint of clinical applicability, we did not establish that the optimal regimen (300 mg every 6 h with infusion for 1 h) determined by Monte Carlo simulation is safe, because only one patient received the optimal regimen in this study. Therefore, it is impossible to evaluate the safety of optimal regimen. However, the simulated maximum steady-state plasma concentration with the optimal regimen in the present study was approximately 11–24 μg/mL. This concentration is lower than the maximum steady-state plasma concentration after administration of 600 mg by 1 h infusion every 12 h in healthy adults (32.4 ± 2.32 μg/mL) [19]. Therefore, it seems likely that the optimal regimen would be safe even in renal dysfunction patients undergoing CHDF. Nevertheless, it will be important to confirm the safety of the optimal regimen, and clinicians should carefully consider the appropriate regimen when administering BIPM to renal dysfunction patients undergoing CHDF, in addition to monitoring for side effects.


In the present study, we used PK modeling to establish the optimal regimen of BIPM in patients with various levels of renal function undergoing CHDF. Clearance of CHDF was determined by the flow rate of dialysate and filtration conditions, since almost BIPM flowed into the dialyzer. The results of Monte Carlo simulation indicated that the regimen of 300 mg every 6 h infusion for 1 h was optimal, compared to other regimens used in our hospital. This regimen showed effective antibacterial activity towards MIC2 μg/mL bacteria.



Continuous hemodiafiltration








Minimum inhibitory concentration

T > MIC:

Time above the MIC


Polymethyl methacrylate




Glomerular filtration rate

QB :

Blood flow rate

QD :

Dialysate flow rate

QS :

Substitution flow rate

QF :

Filtration flow rate




Sieving coefficient


High-performance liquid chromatography


Numeric Analysis Program for Pharmacokinetics


Area under the concentration-time curve


Probability of target attainment


  1. 1.

    Murray P, Hall J. Renal replacement therapy for acute renal failure. Am J Respir Crit Care Med. 2000;162:777–81.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Matsuda K, Hirasawa H, Oda S, Shiga H, Nakanishi K. Current topics on cytokine removal technologies. Ther Apher. 2001;5:306–14.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Sadahiro T, Hirasawa H, Oda S, Shiga H, Nakanishi K, Kitamura N, et al. Usefulness of plasma exchange plus continuous hemodiafiltration to reduce adverse effects associated with plasma exchange in patients with acute liver failure. Crit Care Med. 2001;29:1386–92.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Jiang HL, Xue WJ, Li DQ, Yin AP, Xin X, Li CM, et al. Influence of continuous veno-venous hemofiltration on the course of acute pancreatitis. World J Gastroenterol. 2005;11:4815–21.

    PubMed Central  Article  PubMed  Google Scholar 

  5. 5.

    Hirasawa H. Indications for blood purification in critical care. Contrib Nephol. 2010;166:21–30.

    Article  Google Scholar 

  6. 6.

    Oda S, Sadahiro T, Hirayama Y, Nakamura M, Watanabe E, Tateishi Y, et al. Non-renal indications for continuous renal replacement therapy: Current status in Japan. Contrib Nephol. 2010;166:47–53.

    Article  Google Scholar 

  7. 7.

    Trotman RL, Williamson JC, Shoemaker DM, Salzer WL. Antibiotic dosing in critically ill adult patients receiving continuous renal replacement therapy. Clin Infect Dis. 2005;41:1159–66.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Hidaka S, Goto K, Hagiwara S, Iwasaka H, Noguchi T. Doripenem pharmacokinetics in critically ill patients receiving continuous hemodiafiltration. Yakugaku Zasshi. 2010;130:87–94.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Pollard TA, Lampasona V, Akkerman S, Tom K, Hooks MA, Mullins RE, et al. Vancomycin redistribution: dosing recommendations following high-flux hemodialysis. Kidney Int. 1994;45:232–7.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Ubukata K, Hikida M, Yoshida M, Nishiki K, Furukawa Y, Tashiro K, et al. In vitro activity of LJC10,627, a new carbapenem antibiotic with high stability to dehydropeptidase I. Antimicrob Agents Chemother. 1990;34:994–1000.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  11. 11.

    Hikida M, Kawashima K, Nishiki K, Furukawa Y, Nishizawa K, Saito I, et al. Renal dehydropeptidase-I stability of LJC 10,627, a new carbapenem antibiotic. Antimicrob Agents Chemother. 1992;36:481–3.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  12. 12.

    Yang Y, Testa RT, Bhachech N, Rasmussen BA, Bush K. Biochemical characterization of novel tetrahydrofuranyl 1β-methylcarbapenems: Stability to hydrolysis by renal dehydropeptidases and bacterial β-lactamases, Binding to penicillin binding proteins, and permeability properties. Antimicrob Agents Chemother. 1999;43:2904–9.

    PubMed Central  CAS  PubMed  Google Scholar 

  13. 13.

    Fukasawa M, Sumita Y, Harabe ET, Tanio T, Nouda H, Kohzuki T, et al. Stability of meropenem and effect of 1β-methyl substitution on its stability in the presence of renal dehydropeptidase I. Antimicrob Agents Chemother. 1992;36:1577–9.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  14. 14.

    Craig WA. Pharmacokinetic/pharmacodynamics parameters: rationale for antibacterial dosing mice and men. Clin Infect Dis. 1998;26:1–12.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Drusano GL. Prevention of resistance: a goal for dose selection for antimicrobial agents. Clin Infect Dis. 2003;36:S42–50.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Clinical and laboratory standards institute. Performance standard for antimicrobial susceptibility testing; Twenty-second informational supplement. M100-S22. Wayne: Clinical and Laboratory Standards Institute; 2012. p. 1–184.

    Google Scholar 

  17. 17.

    Nagashima S, Kozawa O, Otsuka T, Kohno K, Minamoto M, Yokokawa M, et al. Pharmacokinetics of a parenteral carbapenem, biapenem, in patients with end-stage renal disease and influence of hemodialysis. J Antimicrob Chemother. 2000;46:839–42.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Ikawa K, Morikawa N, Ikeda K, Suyama H. Pharmacokinetic modeling and dosage adaptation of biapenem in Japanese patients during continuous venovenous hemodiafiltration. J Infect Chemother. 2008;14:35–9.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Nakajima M, Uematsu T, Ueno K, Watanabe K, Nagashima S, Inaba H, et al. Phase 1 study of L-627, biapenem, a new parenteral carbapenem antibiotic. Jpn Pharmacol Ther. 1994;22:255–71.

    Google Scholar 

  20. 20.

    Shitara K, Takai M, Ohta M, Wakahara T, Kanehisa H, Fukunaga T, et al. Development of an equation for predicting body surface area based on three-dimensional photonic image scanning. Jpn J Phys Fitness Sports Med. 2009;58:463–74.

    Article  Google Scholar 

  21. 21.

    Ikawa K, Morikawa N, Ikeda K, Ohge H, Sueda T. Development of breakpoints of cephems for intraabdominal infections based on pharmacokinetics and pharmacodynamics in the peritoneal fluid of patients. J Infect Chemother. 2008;14:141–6.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Kuti JL, Dandekar PK, Nightingale CH, Nicolau DP. Use of Monte Carlo simulation to design an optimized pharmacodynamic dosing strategy for meropenem. J Clin Pharmacol. 2003;43:1116–23.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Mouton JW, Dudley MN, Cars O, Derendorf H, Drusano GL. Standardization of pharmacokinetic/pharmacodynamic (PK/PD) terminology for anti-infective drugs: an update. J Antimicrob Chemother. 2005;55:601–7.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Suyama H, Ikawa K, Morikawa N, Ikeda K, Fujiue Y, Morikawa S, et al. Pharmacokinetics and pharmacodynamics of biapenem in critically ill patients under continuous venovenous hemodiafiltration. Jpn J Antibiot. 2008;61:303–13.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Murakami K. Experimental study of localization of dehydropeptidase-1 in human tissues and hydrolysis of carbapenems. J Wakayama Med Soc. 1993;44:601–15.

    CAS  Google Scholar 

  26. 26.

    De Paepe P, Belpaire FM, Buylaert WA. Pharmacokinetic and pharmacodynamic considerations when treating patients with sepsis and septic shock. Clin Pharmacokinet. 2002;41:1135–51.

    Article  PubMed  Google Scholar 

  27. 27.

    Frei CR, Wiederhold NP, Burgess DS. Antimicrobial breakpoints for gram-negative aerobic bacteria based on pharmacokinetic-pharmacodynamic models with Monte Carlo simulation. J Antimicrob Chemother. 2008;61:621–8.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

Download references


This work was supported in part by JSPS KAKENHI Grants No. 23590173 and in part by a grant from Meiji Seika Pharma Co., Ltd.

Author information



Corresponding author

Correspondence to Yoshimichi Sai.

Additional information

Competing interests

This study was supported in part by Meiji Seika Pharma Co. Ltd.

Authors’ contributions

GA collected samples, performed experiments (HPLC), analyzed data, and drafted the manuscript. YH performed experiments (HPLC) and analyzed data. KI collected samples and analyzed data. KS and KM analyzed data. TS analyzed data, and helped to draft the manuscript. TN designed and coordinated the study, and collected samples. TT designed and coordinated the study. YS designed and coordinated the study, and helped to draft the manuscript. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Akashita, G., Hosaka, Y., Noda, T. et al. PK/PD analysis of biapenem in patients undergoing continuous hemodiafiltration. J Pharm Health Care Sci 1, 31 (2015).

Download citation


  • Biapenem
  • Continuous hemodiafiltration
  • Pharmacokinetics
  • Monte Carlo simulation
  • PK/PD breakpoint