Abstract
It has become commonplace (270+ article citations to date) to measure the fraction unbound (FrUn) of drugs in tissue homogenates and diluted plasma and then use a Correction Factor Equation (CFE) to extrapolate to the undiluted state. The CFE is based on assumptions of nonspecific binding with experimental use of very low drug concentrations. There are several possible determinants of apparent nonspecific binding as measured by methods such as equilibrium dialysis: true macromolecule binding and lipid partitioning along with receptor, enzyme, and transporter interactions. Theoretical calculations based on nonlinear protein binding indicate that the CFE will be most reliable to obtain FrUn when added drug concentration is small, binding constants are weak, protein concentrations are relatively high, and tissue dilution is minimal. When lipid partitioning is the sole factor determining apparent tissue binding, the CFE should be perfectly accurate. Use of very low drug concentrations, however, makes it more likely that specific binding to receptors and other targets may occur, and thus FrUn may reflect some binding to such components. Inclusion of trapped blood can clearly cause minor to marked discrepancies from purely tissue binding alone, which can be corrected. Furthermore, assessment of the occurrence of ionization/pH shifts, drug instability, and tissue metabolism may be necessary. Caution is needed in the use and interpretation of results from tissue dilution studies and other assessments of nonspecific binding, particularly for very strongly bound drugs with very small FrUn values and in tissues with metabolic enzymes, receptors, and trapped blood.
SIGNIFICANCE STATEMENT The use of tissue, plasma, and cell preparations to help obtain fraction unbound and tissue-to-plasma partition coefficients in pharmacokinetics has grown commonplace, especially for brain. This report examines theoretical, physiological, and experimental issues that need consideration before trusting such measurements and calculations.
Introduction
The “free drug (hormone) hypothesis” (Mendel, 1989; Lin, 2006) is commonly invoked to explain many aspects of drug distribution in pharmacokinetics and pharmacodynamics. Thus, the assessment of unbound drug concentrations in plasma and tissues is often sought. Although such measurements are relatively easy in serum or plasma by ultrafiltration or equilibrium dialysis (except for very highly bound drugs), it is challenging to do this in tissues owing to the need for excision and processing of the samples, which particularly necessitate homogenization. The latter results in a mixture of tissue and blood components and dilution of both the drug and tissue components and, if used for assessing fraction unbound (FrUn), requires extrapolation to the original undiluted values [corrected fraction unbound (cfu)]. The employment of diluted tissue homogenates to assess tissue-to-plasma partition coefficients was first proposed by the Hanano laboratory, which used several dilutions for extrapolation to whole tissues (Lin et al., 1982, Harashima et al., 1984). Many tissues and drugs showed good agreement with in vivo values, but some did not.
Kalvass and Maurer (2002) have proposed a Correction Factor Equation (CFE) for adjusting for tissue dilutions in circumstances wherein tissue binding is assumed to be nonspecific, namely linear with drug concentration:(1)in which measured fraction of drug unbound (mfu) is the measured FrUn, and Dil is the single X-fold homogenate dilution factor. This reference has been cited at least 274 times (Web of Science, accessed 08/17/2020), particularly with respect to assessing brain-tissue binding of various CNS drugs and providing implications regarding pharmacology. The methodology is also offered by at least one commercial vendor (Cyprotex, Watertown, MA).
As reviewed (Jusko and Gretch, 1976; Fichtl et al., 1991), numerous tissue components, including albumin, in interstitial fluids can contribute to apparent tissue binding of drugs. Although these will vary with type of tissue, consider the composition of the brain in particular. The brain consists of about 75% water, 5 w/v% protein, and about 5 w/v% lipid. Myelin in normal human brain is largely composed of diverse lipids at 78%–81% (by dry weight) (O’Brien and Sampson, 1965). Brain lipids in rats are composed of fractions of total mass of 0.039 neutral lipids, 0.0015 phospholipids, and 0.036 acid phospholipids along with 0.75 water, which are values that differ with species (Rodgers et al., 2012). Furthermore, the brain and other tissues contain numerous transporters, enzymes, and receptors along with blood, which may contribute to apparent tissue binding of drugs as well, particularly at very low drug concentrations (Ballard et al., 1974; Harashima et al., 1984; Fichtl et al., 1991; Fridén et al., 2010; Stieger and Gao, 2015).
The use of the CFE eq. 1 with tissue homogenates and sometimes in diluted plasma (Kalvass et al., 2018) was originally proposed and has been used by many without full assessment of the conditions in which nonspecific binding will occur. The accuracy of the CFE for various experimental conditions and contributions of possible binding and partitioning components of tissues needs further consideration. This report examines the accuracy of the CFE in the context of possible nonlinear protein binding and simple lipid partitioning, considers the potential impacts of trapped blood and drug metabolism, and addresses the possible contributions from specific binding components, such as receptors, enzymes, and transporters.
Theoretical
Ligand Binding.
A standard relationship (Goldstein, 1949) that describes the nonlinear binding of drug to a single type of macromolecule is:(2)in which Db is concentration of bound drug, Df is concentration of free drug, n is the number of binding sites, Pt is the protein concentration, and KD is the equilibrium dissociation constant of drug and target (proteins, receptors). The concentration of total drug (Dt) is:
(3)
The fraction of drug unbound based on a specific measurement or calculation (fu) is typically measured experimentally as Df/Dt and follows from eqs. 2 and 3 as:(4)
For calculation purposes, it is necessary to generate values of fu from known concentration Dt as the positive root of the quadratic:(5)in which: a = Dt, b = [KD + n·Pt – Dt], and c = -KD.
Equation 4 includes the four indicated determinants of actual protein binding over the full range of possible drug and protein concentrations for one type of binding macromolecule (Pt). The common definition of nonspecific binding is the condition in which Df in eq. 4 is very small and thus negligible in relation to both KD and n·Pt, yielding:(6)
This equation leads to the CFE relationship (eq. 1) in which Dil accounts for differences in n·Pt for the diluted homogenate (Kalvass and Maurer, 2002).
Lipid Partitioning.
The model shown in Fig. 1 depicts an equilibrium dialysis cell in which the tissue homogenate consists of lipid (fat) suspended in an aqueous liquid on one side of a semipermeable membrane with an aqueous medium on the other side. The lipid/aqueous partition coefficient is:(7)
Lipid partitioning: model of determinants of tissue drug distribution and apparent binding for simple lipid [concentration of drug in fat (Cfat)] to aqueous (Cw) partitioning. The partition coefficient Kp = Cfat/Cw. For equilibrium dialysis, the total tissue homogenate concentration CT = Cw.Fw + Cfat.Ffat, in which Fw and Ffat are the fractional aqueous and fat contents.
The measured total drug concentration (CT) in the homogenate side at equilibrium is:(8)
in which Fw is the fraction of tissue or homogenate comprising water, Ffat is the fraction of tissue or homogenate comprising fat, and Fw = 1 − Ffat. This leads to:(9)
Since in this situation, fu = Cw/CT, then(10)
With dilution (Dil) of the homogenate changing values of Ffat, the relationship between the measured homogenate (mfu) and corrected fraction unbound (cfu) in whole tissue is:(11)
This equation is functionally equivalent to CFE eq. 1.
Correction for Trapped Blood.
The inclusion of trapped blood or plasma in a tissue homogenate is unavoidable except by carrying out either whole-body or organ perfusion prior to or just after sacrifice of an animal. This problem is well recognized in analyzing tissue samples in physiologically based pharmacokinetics (PBPK) studies (Khor and Mayersohn, 1991). To correct tissue homogenate–binding measurements for this problem, it is necessary to first obtain FrUn in undiluted plasma or blood [fraction unbound of drug in whole blood or plasma (fuB)]. In turn, the expected FrUn attributed to true tissue in the corrected tissue homogenate (cfu) can be calculated.
Consider the diagram in Fig. 2 to represent the composition of the whole tissue [volume of tissue (VT)], with FB as the fraction of plasma/blood in tissue or homogenate and FT as the fraction of actual tissue separate from blood. The following mass balance would hold:(12)in which Cf is concentration of free drug, CbB is concentration of bound drug in blood, and CbT is concentration of bound drug in tissue. The total tissue-homogenate concentration (CT) is:
(13)
Model for relationship of trapped blood in tissues and role of binding. Concentration of free drug (Cf) is assumed to equilibrate with concentration of bound drug in blood (CbB) and concentration of bound drug in tissue (CbT). FB and fraction of actual tissue separate from blood (FT) are noted. VT, volume of tissue.
The cfu is:(14)
Division by Cf and recognizing the FrUn relationships for blood and remaining tissue:(15a)and
(15b)
yields:(16)
Rearrangement results in:(17)
Note that these calculations need consistency in using only either plasma binding and fraction of plasma in tissues or corresponding whole-blood binding and fractions in these equations. They assume that there will be no cross-interference for binding between plasma/blood and tissue components in the homogenate. Because of the negative sign in eq. 17, this equation will only operate properly when cfu falls between fuB and fuT. Note that eq. 17 will apply after the CFE correction of the tissue homogenate is made.
The composition of the vascular space in brain has been assessed using several methods as reviewed (Friden et al., 2010). Residual blood in rat brain comprises water 10.3, plasma proteins 7.99, and erythrocytes 2.13 μl/g brain, namely about 2% of brain mass. They introduced a correction model for calculating the brain/plasma-unbound partition coefficient.
Receptor Binding.
Drug bound in presence of enzymes, receptors, transporters, or other targets with strong binding in tissues occurring along with nonspecific binding can be described as:(18)
in which Bmax is total receptor concentration, KD is the drug-receptor dissociation constant, and Kns is the nonspecific binding constant of drug (Hazra et al., 2007).
Methods
Previous experimental data for dexamethasone (DEX) binding to receptors in rat liver cytosol were obtained (Hazra et al., 2007) using a developed radioligand-binding assay (Boudinot et al., 1986) with some modifications. Binding assessments of methylprednisolone (MPL) and DEX in rat tissue homogenates were recently carried out (Ayyar et al., 2019b).
Results
Ligand Binding.
Simulations were performed using eqs. 4 and 5 to assess the expected values of FrUn in relation to varying drug and protein concentrations for hypothetical drugs with weak KD = 10−4 M, moderate KD = 10−5 M, and strong KD = 10−7 M equilibrium dissociation constant values. The standard practice in applying the CFE for brain tissue is to employ an added drug concentration of 1 or 5 µM (Kalvass et al., 2007; Friden et al., 2011). To mimic the presence of a very low or trace drug concentration (eq. 6) and apparent nonspecific binding, a drug concentration of 0.01 µM was compared with an assumed added concentration of 10 μM. Ultimately, the ability of the CFE (eq. 1) to recover the “true” FrUn was assessed. For these and subsequent simulations, we assumed one class of binding proteins as described by eq. 2, with the protein having the molecular weight of albumin (69,000).
Figure 3 shows values of FrUn in relation to a wide range of protein concentrations for drugs with three values of KD in which the tracer and experimental drug concentrations were used in operation of eqs. 4 and 5. As expected, the lowest KD value produced strongest binding and the lowest values of FrUn, ranging from about 0.5 at very low-protein concentrations to fu = 0.0001 at the Pt of 10−3 M. The latter is approximately the total protein concentration in brain of 6.9%. The FrUn is similar for the two drug concentrations at the higher protein concentrations and more weakly bound drugs. However, there is increasing divergence of FrUn values for the two drug concentrations when protein concentrations are low, especially for the more strongly bound compound. It can be noted that when starting with a set Dt concentration, there will be generation of lower FrUn values because of reductions in Df as binding increases.
Ligand binding: relationship of fu vs. Pt for compounds with three different KD values according to eq. 2 for the indicated values of Df.
Figure 4 provides simulations in which FrUn was calculated for a wide range of total drug concentrations but when values were determined for relatively low 10−6 M and high 10−4 M protein concentrations. The FrUn values are of course lowest when KD was strongest (10−7 M) and protein concentrations were higher. The values of FrUn were nonlinear with drug concentration as expected, but all of the profiles show constant FrUn values for a wide range of very low drug concentrations. The determination of possible consistency of FrUn thus requires consideration of all three contributing factors of Dt, Pt, and KD.
Ligand binding: relationship of fu vs. Df for two Pt values according to eq. 2.
Figure 5 depicts the differences expected in FrUn over a range of protein concentrations when added drug concentrations are 10 versus 0.01 µM. The latter concentration will produce least saturation of binding and best reflect strongest (lowest FrUn) binding conditions. The smallest differences are found at relatively high protein concentrations (Pt = 10−4 to 10−3 M), and substantial differences are found when protein concentrations are very low. Interestingly, the drug with moderate binding (KD = 10−5 M) exhibited a wider range of FrUn ratios closest to a 1% difference compared with the others.
Ligand binding: differences in fu in relation to Pt for compounds with three KD values when drug concentrations are 10 and 0.01 µM according to eq. 2.
Figure 6 demonstrates the ability of the CFE to recover the true FrUn for the three binding constants over the range of protein concentrations when it was assumed that the drug concentration added was 10 µM, and the tissue was diluted 3-fold as is commonly done. In this case, the true FrUn was calculated using eqs. 4 and 5 for a given drug and protein concentration, and a “measured” FrUn was obtained with the Pt reduced to one-third. In turn, the calculated or experimental FrUn was generated using eq. 1. The graph provides the ratio of the cfu/true fu. The CFE was found to function well for drugs with all KD values when protein concentrations were high but diverged appreciably for the stronger binding compound as protein concentrations decreased.
Ligand binding: ratio of experimental fu to true fu in relation to Pt for compounds with three KD values. The experimental fu was calculated based on the CFE (eq. 1) assuming a Dil of 3.
These simulations indicate that for simple binding alone, the reliability of the CFE will depend strongly on utilizing the lowest possible drug concentration to seek linearity and greatest degree of binding but will be most dependable for more weakly binding drugs and when the binding protein is present in relatively higher concentrations. These conditions are helped by diluting the tissue homogenate as little as possible.
Lipid Partitioning.
Simulations of apparent FrUn were performed based on the assumption that equilibrium dialysis would be performed (Fig. 1) using diluted tissue homogenates with a classic partition coefficient [lipid/water partition coefficient or in vivo tissue/plasma ratio (Kp)] and lipid partitioning solely responsible for drug concentrations in the aqueous (Cw) and homogenate (CT) chambers. Figure 7 shows the relationship of apparent FrUn versus the fraction of fat in the homogenate for differing values of Kp. Although the Ffat axis is scaled to 1.0, practical values for undiluted tissue would range from about 0.05 for brain to 0.8 for adipose tissue (Rodgers et al., 2012). As expected, higher values of Kp produce smaller values of FrUn. It can be noted that the largest Kp represents a logP of 6, and the range includes a wide array of drugs. These simulations assume that the drug is not ionized where the situation becomes more complicated.
Lipid partitioning: relationship between fu and Ffat for compounds with the listed Kp.
Figure 8 was constructed in a similar manner to Fig. 6 in providing an assessment of how well the CFE would recover the true FU in a situation in which 10-µM concentrations of drug were added to a tissue homogenate diluted 3-fold. A “measured” FrUn was generated using eq. 10 for the dilution, and then eq. 1 was employed to obtain the calculated cfu. The ratios of cfu to true fu are graphed in relation to Ffat in Fig. 8. For all values of Kp and Ffat, the CFE produced exact predictions of the true FrUn. Since lipid partitioning was assumed to be a linear process with constant Kp values, this would be true for any added drug concentration.
Presence of Trapped Blood/Plasma.
Figure 9 shows simulations using eq. 17 of possible effects of trapped blood/plasma in tissue homogenates used to assess tissue binding. The calculations assume that the true fuT = 0.05 and consider the effects of various degrees of blood/plasma-binding fuB values as denoted for the fractions of FB ranging from 0 to 0.5 as shown in the graph. It can be seen that deviations from true fuT values can vary in both the positive and negative directions depending on fuB in relation to fuT. When blood/plasma binding is less than tissue binding, there is a dilution effect from the trapped blood that skews the deviation in the positive direction. When plasma binding is stronger than tissue binding, the deviation is skewed in the negative direction, sometimes considerably when blood FrUn is very strong. When they are similar, there is no correction needed, as indicated by the horizontal line in the graph. This simulation shows only one set of possibilities since the variable degrees of binding and tissue compositions will all contribute to a wide range of possible deviations from true values. It can be noted that some drugs have an unbound fraction in plasma (fup) as low as 0.000012 (venetoclax), which requires special approaches for measurement (Kalvass et al., 2018).
Deviations from true tissue-binding fuT with varying FB and various degrees of blood/plasma binding (as indicated by fuB values) present in the tissue homogenate.
Specific Target Binding.
An example of the combination of specific receptor binding and nonspecific binding that can occur for many ligands and tissues is shown in Fig. 10. This graph depicts an experimental assessment of dexamethasone binding to components of rat liver cytosol (Hazra et al., 2007). To quantitate joint glucocorticoid receptor binding and nonspecific binding, the radiolabeled drug (H3) is added along with a range of Dt values, bound drug-receptor complex (Db) is separated by ultracentrifugation, and the relationship between Db and Df is assessed. The contribution of Kns is determined by saturating receptor binding with added unlabeled drug in very high concentrations (100× the values on the y-axis) in companion samples. The operative binding equation in this case is eq. 18.
Specific target binding. Relationship of Db vs. Df in an experiment where H3-dexamethasone binding to glucocorticoid receptors in rat liver cytosol was assessed for calculation of Bmax and affinity (KD). The free-drug concentration scale for nonspecific binding was about 100-fold higher than that labeled. Joint assessment of total and Kns was necessary to resolve all binding parameters. In this study, values of Bmax = 1079 fmol, KD = 844 fmol, and Kns = 0.021 fmol were calculated. From Hazra et al. (2007).
Joint fitting of the data when there is a range of low and very high drug concentrations allows identification of Bmax, KD, and Kns. This example demonstrates that the common practice of screening tissue homogenates for calculation of FrUn using a single very low drug concentration does not necessarily reflect nonspecific binding but may include measurement of binding to possible receptors or other biologic targets.
Presence of Metabolic Enzymes or Drug Instability.
Figure 11 depicts the metabolic and binding assessments of the corticosteroid MPL in diluted homogenates freshly prepared from male rat livers. The metabolic loss rate of MPL was more rapid in the less diluted homogenates but linear with time and dilution. These data were used to assess hepatic clearance of MPL (Ayyar et al., 2019b). To quantitate the FrUn, the tissue homogenates were diluted (3- to 10-fold), and drug was added at a concentration of 10 µg/ml. Preliminary validation experiments performed at 1 µg/ml yielded a similar Db/Df value as compared with 10 µg/ml, suggesting concentration-independent binding within this range. These concentrations correspond well with observed MPL liver exposures in our in vivo PK-pharmacodynamic assessments in rats (Ayyar et al., 2019a). Then samples were loaded into ultrafiltration devices and incubated at 37°C for a predetermined time period (to achieve an equilibrium in binding), free and bound drug was separated by ultrafiltration at 37°C, and the relationship between values of the binding ratio Db/Df and the Dil were assessed at each dilution. The observed linear relationship can be derived from eq. 2 when Df is very low:(19)in which the slope = n·Pt/KD.
(Top) Time course of in vitro stability of MPL in liver homogenates prepared at 3- (red), 4- (blue), 6- (orange), and 10-fold (green) dilutions from freshly harvested male livers. Symbols are the mean ± S.D. (n = 3 per time point). (Bottom) Binding of MPL in homogenates prepared from male rat liver at initial concentrations of 10 µg/ml. Symbols depict the mean ± S.D. of the binding ratios (Db/Df) across four dilutions of tissue homogenate. The dashed line represents the best-fit line (eq. 19) extrapolated to an undiluted state of tissue (Dil = 1). From Ayyar et al. (2019b).
With back-extrapolation to unity (i.e., an undiluted tissue state), fuT was computed using the relationship (Lin et al., 1982; Ayyar et al., 2019b):(20)
As shown in Fig. 11, MPL binding was proportional to the total protein concentration in liver homogenates across 3- to 10-fold dilutions. Employing the method in eq. 2, the fuT for MPL in liver was 0.162 ± 0.01. Since most screening studies employ a single dilution, the measured FrUn was also corrected to the undiluted state using eq. 1. Using only the 3-fold dilution values, the CFE method yielded a value of 0.150 ± 0.01. The use of several dilutions of homogenates is advantageous in providing more robust data for extrapolation and confirming linearity of binding with varying tissue dilutions. The plasma binding of MPL at 37°C was linear with moderate binding of 61% (fup = 0.39). Our studies of corticosteroids are not complicated by changes or differences in tissue pH since these are neutral compounds. It can be noted that the published calculation methods (Poulin and Theil, 2002; Berezhkovskiy, 2004; Rodgers and Rowland, 2006) for estimating tissue:plasma ratios yielded reasonable estimates of Kp in muscle and lung as compared with our directly measured values, but estimates in liver differed appreciably (Ayyar et al., 2019b). Although we did not, at the time, recognize the importance of trapped blood, recalculation using the CFE and eqs. 16 and 17 (blood/plasma partition coefficient equal to 1) revealed that tissue FrUn values for MPL are essentially the same in muscle but change from 0.17 to 0.14 in lung, 0.16 to 0.13 in male liver, and 0.063 to 0.055 in female rat liver. Corresponding values for DEX also changed from 0.12 to 0.11 in lung and 0.066 to 0.059 in liver (Ayyar et al., 2019b).
Discussion
Binding to plasma and tissues is a major determinant of drug distribution in the body and thus plays an important role in pharmacokinetics. The ratio of fup and fuT generally defines the Kp, an important parameter in PBPK that determines total drug concentrations within tissues and the volume of distribution at steady-state (Vss). Since Kp is often governed by reversible binding to macromolecules in blood and tissues, in vitro assessments of drug binding in these media are used to estimate tissue Kp values in vivo. However, in vivo drug uptake into tissues may be more complicated because the “permissive” availability from both plasma-binding proteins and transporter-mediated tissue influx can occur (Herve et al., 1994). Transporter-mediated efflux is also commonplace, particularly in brain and liver.
Numerous tissue components can bind drugs, including immunophilins, phospholipids, microsomes, mitochondria, DNA, and others (Jusko and Gretch, 1976, Fichtl et al., 1991). Organs like muscle and heart contain contractile proteins that bind cardiac glycosides. A PBPK model for doxorubicin included nonlinear binding to DNA and cardiolipin (Gustafson et al., 2002). Many compounds exhibit nonlinear target-mediated binding and disposition (Mager and Jusko, 2001), including vildagliptin that binds to dipeptidyl peptidase IV both in plasma and tissues (Landersdorfer et al., 2012). Albumin and globulins are both found in interstitial fluids at concentrations that are usually about one-half of that in plasma (Jusko and Gretch, 1976). Interestingly, the methods of estimating tissue Kp values attribute involvement of the three types of tissue lipids and thus invoke simple partitioning (Rodgers and Rowland, 2006). Uncertainties exist if the tissue dilution and low drug concentrations used in measurements are adequate to allow for extrapolations of FrUn to whole-tissue values. Notably, our simulations—on one hand—argue for use of low drug concentrations and minimal dilutions (higher protein concentrations); our demonstration—on the other hand—of receptor binding (Fig. 10) indicates that drug concentrations considerably higher than KD values (usually nanomolars) are needed to see linearity.
Our simulations that assess the CFE method for drug binding to tissue components with possibly nonlinear binding demonstrate generally good but not perfect agreement of the methodology with theory (Figs. 3⇑⇑–6). Similar simulations assessing the role of protein concentration affecting either FrUn or fraction bound were published (Jusko and Gretch, 1976; Wan et al., 2007). All show that there is a range of protein concentrations (and thus dilutions) that will produce similar FrUn values, but caution is needed that such range may not be foreseen without experimentation. This was done in some studies in tissues (Lin et al., 1982; Wan et al., 2007; Ayyar et al., 2019b) and plasma (Kalvass et al., 2018).
If fat distribution with linear partition coefficients accounts for tissue “binding,” then the CFE method works perfectly in theory (Fig. 8). Body tissues of rats contain various percentages of neutral lipids, with the fat content ranging from about 4% in brain to 86% in adipose tissues (Rodgers et al., 2012). The contribution of the Ffat to values of FrUn in tissues varies with the partition coefficient and fat content (Fig. 7). Drugs exhibit a wide range of logP and logD values, and this factor is considered in computational methods that predict tissue partition coefficients. Our simulations extend to tissue distribution expectations for compounds with very high logP values (6) in which the in vivo Kp is determined by the relative concentrations of lipid in plasma/blood versus tissue. This appears to produce an upper limit in in vivo Kp values (Haddad et al., 2000). However, the CFE appears to work reasonably well in brain for many neutral compounds owing to its low-protein, high-lipid, and minimal–trapped blood contents.
Trapped blood offers a major possible artifact when examining tissue samples from in vivo studies to obtain Kp values (Khor and Mayersohn, 1991). It is commonplace in PBPK studies to correct for trapped blood using values for fractional tissue blood space (Bernareggi and Rowland, 1991). This is clearly a possible artifact for using tissue homogenates and slices to assess binding as well. The presence of albumin (and perhaps other plasma proteins) at a concentration of about 0.8–1.0 µM in brain homogenates is considered to reflect contamination from trapped blood (Longhi et al., 2011; Loryan et al., 2016). Mouse-brain homogenates were found to contain about 1% trapped blood using chromium-labeled red blood cells (Garg and Balthasar, 2009). Except for one report (Friden et al., 2010), this problem has largely gone unrecognized for tissue homogenates. Although this artifact seemingly can be obviated by perfusion of the animal organ or whole body with saline (PBS) before processing the tissue (Taves et al., 2010), this process produces a dilution factor if the perfusion fluid remains “trapped” in the tissue (similar to 0 FrUn binding shown in Fig. 9). Studies with 36 compounds were carried out in rats in which whole-body perfusions with PBS were carried out before collection of 14 sets of tissues (Berry et al., 2010). Diluted tissue homogenates were analyzed, and the CFE was applied to obtain fuT and ultracentrifugation to measure fup. Values of Vss factoring in the blood/plasma partition coefficient (B/P) ratio were generated by:(21)
These summed Vss values were within 2-fold of PK-calculated values for 77% of neutral compounds and for 61% for all compounds (neutral, acids, bases, zwitterions). There is clearly room for improvement in using tissue homogenates for predicting tissue binding and Vss values.
Our assessments pertain to simple factors that control tissue binding and accuracy of the CFE method and do not account for additional determinants, such as ionization for acids and bases, relevant pH gradients, subcompartments, and other differences among tissues. Ionization can be accounted for when using computational methods when the pKa of the drug is employed (Poulin and Theil, 2000; Berezhkovskiy 2004; Rodgers and Rowland, 2006; Poulin, 2015). Brain-slice measurements of tissue partition coefficients that preserve tissue architecture are considered more reliable than tissue homogenates owing to the better integrity of drug-pH partitioning (Friden et al., 2007, 2011). These authors also provide useful insights into tissue subcompartments that may attain far differing drug concentrations than those indicated by in vitro measurements that destroy tissue integrity. Weak acids will be primarily found in interstitial fluids owing to their extensive ionization that reduces cellular access. Weak bases that are more highly ionized at lower pH values will concentrate in lysosomes with low pH (4.5–5.5) compared with the cytosol pH of 7.0–7.5 (Yokogawa et al., 2002; Trapp et al., 2008).
The role of trapped blood in skewing use of homogenates for tissue binding will, of course, vary with the type of tissue or organ and its content of trapped blood (Fig. 9). Although the brain and muscle contain the smallest fractions of trapped blood, the compounds that are most effective in the CNS tend to be lipophilic weak bases with high degrees of plasma protein and lipid binding. Large deviations from the true tissue FrUn are expected even with small percentages of trapped blood when plasma or blood FrUn is much larger than tissue FrUn (Fig. 9). Most of the 42 CNS-active compounds assessed in mouse plasma and brain homogenates generally showed fuT ≥ fup values (Wan et al., 2007). Furthermore, 108 diverse compounds demonstrated a strong correlation (r = −0.78) of fuT with cLogP values.
Although the CFE method is most commonly applied to brain homogenates, the brain has ATP-binding cassette transporters (Stieger and Gao, 2015; Gomez-Zepeda et al., 2019). Several studies have compared Kp values obtained from brain homogenates with in vivo measured values using a variety of compounds. Good concordance for nine compounds was found when passive CNS distribution was expected, but appreciable overestimates of fu-plasma to fu-brain ratios occurred for compounds with an active efflux mechanism (Kalvass and Maurer, 2002). Brain homogenate FrUn values were assessed for 56 compounds, of which 13 were neutral compounds not subject to pH considerations (Friden et al., 2011). The ratios of observed to predicted FrUn values were within 20% for eight neutral compounds and off by varying degrees for the five others. Assessments of brain distribution of CNS drugs for 33 compounds showed that 23 exhibited unbound brain-to-plasma ratios that were within 3-fold (Kalvass et al., 2007). Most were weak bases, but the acidic and neutral compounds showed closer agreement. Thus, there is considerable uncertainty in use of tissue-binding measurements to estimate in vivo Kp values.
Drug loss in tissue homogenates prepared from metabolizing organs can significantly impact measurements of FrUn in tissue when times of processing and equilibrium dialysis are factors. For example, brain tissue in human and rat can convert alprazolam to active hydroxy metabolites (Agarwal et al., 2008). This effect was apparent during our assessments of MPL binding, particularly with use of liver homogenates (Fig. 11). To reduce the impact of drug metabolism or degradation in various media, some homogenate-binding studies are conducted at 4°C (Lin et al., 1982; Berry et al., 2010). However, the unbound fractions for several drugs tend to increase at lower temperatures (Ballard, 1974). In addition to metabolic instability, chemical instability may need confirmation. Sirolimus is unstable, with a half-life of about 15 hours in rat whole blood, which likely will also occur in tissues and homogenates (Ferron and Jusko, 1998).
Methods that are frequently used to assess in vitro drug metabolism employ microsomes, hepatocytes, and occasionally liver slices. It is a common practice to measure nonspecific binding to calculate unbound intrinsic clearance in microsomes (Obach, 1999; Austin et al., 2002). These studies typically employ equilibrium dialysis (over 5 hours) with a single low drug concentration (e.g., 1.0 µM) and varying (0.3–10 mg/ml) microsomal protein. Constant nonspecific binding and drug stability are assumed, but tissue uptake of drug results in a dilution effect during dialysis. When using hepatocytes, the incubation medium may contain no serum proteins, 10% fetal calf serum, or whole human serum (Riley et al., 2005). Drug binding in the cell cytosol is only sometimes measured (Naritomi et al., 2003). Again, single drug concentrations (1 μM) are used, and constant, nonspecific binding is assumed. Equations using these data for prediction of in vivo clearances employ both the in vitro intrinsic clearances and FrUn values. The ability for making these predictions fare no better than yielding average differences of 4-fold (Riley et al., 2005; Wood et al., 2017). More rigorous attention to obtaining appropriate in vitro FrUn values has been called for (Wood et al., 2017).
Summary
Several conditions must be met for suitable application of diluted tissue homogenates along with the CFE (eq. 1) in measuring, predicting, and interpreting the tissue unbound fractions of drugs: 1) nonspecific and linear binding to macromolecules in plasma and tissues, along with lipid partitioning, dominate the tissue retention process; 2) drug distributes and binds relatively uniformly within each organ; 3) trapped blood drug contributions are either minor (e.g., brain) or are corrected for; 4) nonlinear processes, such as receptor binding, are noncontributory; 5) active influx or efflux transport do not contribute significantly; and 6) metabolic (elimination) processes and chemical instability are minimal. In vitro assessment of tissue binding or assessment of nonspecific binding using at least three dilutions or concentrations can help confirm linearity for extrapolation purposes. Some of these factors may contribute to the difficulty observed in using tissue homogenates to predict in vivo tissue partition coefficients and in assessing FrUn in microsomal and hepatocyte media. Although confirmation of tissue partition coefficients with in vivo studies is always advisable, the combined use of tissue preparations, in vivo measurements, and computational methods can provide unique insights into determinants of drug distribution (Ayyar et al., 2019b).
Authorship Contributions
Participated in research design: Jusko, Molins, Ayyar.
Conducted experiments: Molins, Ayyar.
Performed data analysis: Molins, Ayyar.
Wrote or contributed to the writing of the manuscript: Jusko, Ayyar
Footnotes
- Received May 14, 2020.
- Accepted July 23, 2020.
This work was supported by National Institutes of Health [Grant R35-131800].
Abbreviations
- CFE
- Correction Factor Equation
- cfu
- corrected fraction unbound
- CNS
- central nervous system
- CT
- total drug concentration
- Cw
- concentration of drug in water
- Db
- concentration of bound drug
- DEX
- dexamethasone
- Df
- concentration of free drug
- Dil
- homogenate dilution factor
- Dt
- concentration of total drug
- FB
- fraction of plasma/blood in tissue or homogenate
- Ffat
- fraction of tissue or homogenate comprising fat
- FrUn
- fraction unbound
- fu
- fraction of drug unbound based on a specific measurement or calculation
- fuB
- fraction unbound of drug in whole blood or plasma
- fup
- unbound fraction in plasma
- fuT
- fraction unbound of drug in whole tissue(s)
- Fw
- fraction of tissue or homogenate comprising water
- KD
- equilibrium dissociation constant of drug and target (proteins, receptors)
- Kns
- nonspecific binding constant of drug
- Kp
- lipid/water partition coefficient or in vivo tissue/plasma ratio
- MPL
- methylprednisolone
- PBPK
- physiologically based PK PD pharmacodynamic
- PK
- pharmacokinetics
- Pt
- protein concentration
- Vss
- volume of distribution at steady state
- VT
- volume of tissue
- Copyright © 2020 by The American Society for Pharmacology and Experimental Therapeutics