Are all measures of liver Kpuu a function of FH, as determined following oral dosing, or have we made a critical error in defining hepatic drug clearance?

Here we present, utilizing universally accepted relationships for hepatic clearance at steady state, that for all models of hepatic elimination the ratio of unbound liver drug concentration to unbound systemic blood concentration, Kpuu, is a function of or related to the hepatic bioavailability for that drug, FH. According to the derivation for the well-stirred model, Kpuu can never exceed unity, can frequently be a function of hepatic blood flow, and is equivalent to the value of FH as determined following oral dosing. For the parallel tube model, Kpuu will not equal FH but will be a function of FH and will also never be a value greater than 1. When hepatic clearance is rate limited by basolateral transporters, Kpuu will be less than 1, and less than FH. We believe that such outcomes are highly unlikely, and that the error arises from a basic assumption concerning hepatic clearance that leads to the mechanistic models of hepatic elimination, the well-stirred, parallel tube and dispersion models. That basic assumption is that the steady-state systemic concentration multiplied by the hepatic systemic clearance is equal to the product of the average unbound liver steady-state concentration and the intrinsic hepatic clearance (Css · CL = CH,u · CLint). Calculations of Kpuu and FH based on present methods of analysis provide a strong argument as to why this universally accepted relationship is not correct. Alternatively, we have shown in recent publications that hepatic clearance may be adequately determined based on Kirchhoff's Laws where no assumption of the above equality concerning hepatic intrinsic clearance is required, and where Kpuu is independent of hepatic extraction ratio and FH.

A Complete Extension of Classical Hepatic Clearance Models Using Fractional Distribution Parameter fd in Physiologically Based Pharmacokinetics

The classical organ clearance models have been proposed to relate the plasma clearance CLp to probable mechanism(s) of hepatic clearance. However, the classical models assume the intrinsic capability of drug elimination (CLu,int) that is physically segregated from the vascular blood but directly acts upon the unbound drug concentration in the blood (fubCavg), and do not handle the transit-time delay between the inlet/outlet concentrations in their closed-form clearance equations. Therefore, we propose unified model structures that can address the internal blood concentration patterns of clearance organs in a more mechanistic/physiological manner, based on the fractional distribution parameter fd operative in PBPK. The basic partial/ordinary differential equations for four classical models are revisited/modified to yield a more complete set of extended clearance models, i.e., the Rattle, Sieve, Tube, and Jar models, which are the counterparts of the dispersion, series-compartment, parallel-tube, and well-stirred models. We demonstrate the feasibility of applying the resulting extended models to isolated perfused rat liver data for 11 compounds and an example dataset for in vitro-in vivo extrapolation of the intrinsic to the systemic clearances. Based on their feasibilities to handle such real data, these models may serve as an improved basis for applying clearance models in the future.

In Silico and In Vivo Pharmacokinetic Evaluation of 84-B10, a Novel Drug Candidate against Acute Kidney Injury and Chronic Kidney Disease

Acute kidney injury (AKI) and chronic kidney disease (CKD) have become public health problems due to high morbidity and mortality. Currently, drugs recommended for patients with AKI or CKD are extremely limited, and candidates based on a new mechanism need to be explored. 84-B10 is a novel 3-phenylglutaric acid derivative that can activate the mitochondrial protease, Lon protease 1 (LONP1), and may protect against cisplatin-induced AKI and unilateral ureteral obstruction- or 5/6 nephrectomy [5/6Nx]-induced CKD model. Preclinical studies have shown that 84-B10 has a good therapeutic effect, low toxicity, and is a good prospect for further development. In the present study, the UHPLC-MS/MS method was first validated then applied to the pharmacokinetic study and tissue distribution of 84-B10 in rats. Physicochemical properties of 84-B10 were then acquired in silico. Based on these physicochemical and integral physiological parameters, a physiological based pharmacokinetic (PBPK) model was developed using the PK-Sim platform. The fitting accuracy was estimated with the obtained experimental data. Subsequently, the validated model was employed to predict the pharmacokinetic profiles in healthy and chronic kidney injury patients to evaluate potential clinical outcomes. Cmax in CKD patients was about 3250 ng/mL after a single dose of 84-B10 (0.41 mg/kg), and Cmax,ss was 1360 ng/mL after multiple doses. This study may serve in clinical dosage setting in the future.

Kirchhoff's Laws and Hepatic Clearance, Well-Stirred Model - Is There Common Ground?

Clearance concepts are extensively applied in drug development and drug therapy. The well-stirred model (WSM) of hepatic elimination is the most widely adopted physiologic model in pharmacokinetics owing to its simplicity. A common feature of this organ model is its use to relate hepatic clearance of a compound to the physiologic variables: organ blood flow rate, binding within blood, and hepatocellular metabolic and excretory activities. Recently, Kirchhoff's laws of electrical network have been applied to organ clearance (Pachter et al., 2022; Benet and Sodhi, 2023) with the claim that they yield the same equation for hepatic clearance as the WSM, and that the equation is independent of a mechanistic model. This commentary analyzes this claim and shows that implicit in the application of Kirchhoff's approaches are the same assumptions as those of the WSM. Concern is also expressed in the interpretation of permeability or transport parameters and related equations, as well as the inappropriateness of the corresponding equation defining hepatic clearance. There is no value, and some dangers, in applying Kirchhoff's electrical laws to organ clearance. SIGNIFICANCE STATEMENT: This commentary refutes this claim by Pachter et al. (2022), and Benet and Sodhi, (2023), who suggest that the well-stirred model (WSM) of hepatic elimination, the most widely applied physiologic model of hepatic clearance, provides the same equation as Kirchhoff's laws of electrical network that is independent of a physiologic model. A careful review shows that the claim is groundless and fraught with errors. We conclude that there is no place for the application of Kirchhoff's laws to organ clearance models.

Utility of Minimal Physiologically Based Pharmacokinetic Models for Assessing Fractional Distribution, Oral Absorption, and Series-Compartment Models of Hepatic Clearance

Minimal physiologically based pharmacokinetic (mPBPK) models are physiologically relevant, require less information than full PBPK models, and offer flexibility in pharmacokinetics (PK). The well-stirred hepatic model (WSM) is commonly used in PBPK, whereas the more plausible dispersion model (DM) poses computational complexities. The series-compartment model (SCM) mimics the DM but is easier to operate. This work implements the SCM and mPBPK models for assessing fractional tissue distribution, oral absorption, and hepatic clearance using literature-reported blood and liver concentration-time data in rats for compounds mainly cleared by the liver. Further handled were various complexities, including nonlinear hepatic binding and metabolism, differing absorption kinetics, and sites of administration. The SCM containing one to five (n) liver subcompartments yields similar fittings and provides comparable estimates for hepatic extraction ratio (ER), prehepatic availability (Fg ), and first-order absorption rate constants (ka ). However, they produce decreased intrinsic clearances (CLint ) and liver-to-plasma partition coefficients (Kph ) with increasing n as expected. Model simulations demonstrated changes in intravenous and oral PK profiles with alterations in Kph and ka and with hepatic metabolic zonation. The permeability (PAMPA P) of the various compounds well explained the fitted fractional distribution (fd ) parameters. The SCM and mPBPK models offer advantages in distinguishing systemic, extrahepatic, and hepatic clearances. The SCM allows for incorporation of liver zonation and is useful in assessing changes in internal concentration gradients potentially masked by similar blood PK profiles. Improved assessment of intraorgan drug concentrations may offer insights into active moieties driving metabolism, biliary excretion, pharmacodynamics, and hepatic toxicity. SIGNIFICANCE STATEMENT: The minimal physiologically based pharmacokinetic model and the series-compartment model are useful in assessing oral absorption and hepatic clearance. They add flexibility in accounting for various drug- or system-specific complexities, including fractional distribution, nonlinear binding and saturable hepatic metabolism, and hepatic zonation. These models can offer improved insights into the intraorgan concentrations that reflect physiologically active moieties often driving disposition, pharmacodynamics, and toxicity.

Modeled Rat Hepatic and Plasma Concentrations of Chemicals after Virtual Administrations Using Two Sets of in Silico Liver-to-Plasma Partition Coefficients

The hepatic elimination of chemical substances in pharmacokinetic models requires hepatic intrinsic clearance (CLh,int) parameters for unbound drug in the liver, and these are regulated by the liver-to-plasma partition coefficients (Kp,h). Both Poulin and Theil and Rodgers and Rowland have proposed in silico expressions for Kp,h for a variety of chemicals. In this study, two sets of in silico Kp,h values for 14 model substances were assessed using experimentally reported in vivo steady-state Kp,h data and time-dependent virtual internal exposures in the liver and plasma modeled by forward dosimetry in rats. The Kp,h values for 14 chemicals independently calculated using the primary Poulin and Theil method in this study were significantly correlated with those obtained using the updated Rodgers and Rowland method and with reported in vivo steady-state Kp,h data in rats. When pharmacokinetic parameters were derived based on individual in vivo time-dependent data for diazepam, phenytoin, and nicotine in rats, the modeled liver and plasma concentrations after intravenous administration of the selected substrates in rats using two sets of in silico Kp,h values were mostly similar to the reported time-dependent in vivo internal exposures. Similar results for modeled liver and plasma concentrations were observed with input parameters estimated by machine-learning systems for hexobarbital, fingolimod, and pentazocine, with no reference to experimental pharmacokinetic data. These results suggest that the output values from rat pharmacokinetic models based on in silico Kp,h values derived from the primary Poulin and Theil model would be applicable for estimating toxicokinetics or internal exposure to substances.

Exploring the Pharmacokinetic Mysteries of the Liver: Application of Series Compartment Models of Hepatic Elimination

Among the basic hepatic clearance models, the dispersion model (DM) is the most physiologically sound compared with the well-stirred model and the parallel tube model. However, its application in physiologically-based pharmacokinetic (PBPK) modeling has been limited due to computational complexities. The series compartment models (SCM) of hepatic elimination that treats the liver as a cascade of well-stirred compartments connected by hepatic blood flow exhibits some mathematical similarities to the DM but is easier to operate. This work assesses the quantitative correlation between the SCM and DM and demonstrates the operation of the SCM in PBPK with the published single-dose blood and liver concentration-time data of six flow-limited compounds. The predicted liver concentrations and the estimated intrinsic clearance (CLint ) and PBPK-operative tissue-to-plasma partition coefficient (Kp ) values were shown to depend on the number of liver sub-compartments (n) and hepatic enzyme zonation in the SCM. The CLint and Kp decreased with increasing n, with more remarkable differences for drugs with higher hepatic extraction ratios. Given the same total CLint , the SCM yields a higher Kp when the liver perivenous region exhibits a lower CLint as compared with a high CLint at this region. Overall, the SCM nicely approximates the DM in characterizing hepatic elimination and offers an alternative flexible approach as well as providing some insights regarding sequential drug concentrations in the liver. SIGNIFICANCE STATEMENT: The SCM nicely approximates the DM when applied in PBPK for characterizing hepatic elimination. The number of liver sub-compartments and hepatic enzyme zonation are influencing factors for the SCM resulting in model-dependent predictions of total/internal liver concentrations and estimates of CLint and the PBPK-operative Kp . Such model-dependency may have an impact when the SCM is used for in vitro-to-in vivo extrapolation (IVIVE) and may also be relevant for PK/PD/toxicological effects when it is the driving force for such responses.