Physiologically Based Pharmacokinetic Modeling: The Reversible Metabolism and Tissue-Specific Partitioning of Methylprednisolone and Methylprednisone in Rats.

The pharmacokinetics (PK) of methylprednisolone (MPL) exhibited tissue-specific saturable binding and reversible conversion with its metabolite, methylprednisone (MPN). Blood and 11 tissues were collected in male rats after intravenous (i.v.) bolus doses of 50 mg/kg MPL and 20 mg/kg MPN and upon i.v. infusion of MPL and MPN at 0.3, 3, and 10 mg/h per kg. The concentrations of MPL and MPN were simultaneously measured. A comprehensive physiologically based pharmacokinetic (PBPK) model was applied to describe the plasma and tissue profiles and estimate PK parameters of the MPL/MPN interconversion system. Both dosed and formed MPL and MPN were in rapid equilibrium or achieved steady-state rapidly in plasma and tissues. MPL tissue partitioning was nonlinear, with highest capacity in liver (322.9 ng/ml) followed by kidney, heart, intestine, skin, spleen, bone, brain, muscle, and lowest in adipose (2.74 ng/ml) and displaying high penetration in lung. The tissue partition coefficient of MPN was linear but widely variable (0.15∼5.38) across most tissues, with nonlinear binding in liver and kidney. The conversion of MPL to MPN occurred in kidney, lung, and intestine with total clearance of 429 ml/h, and the back conversion occurred in liver and kidney at 1342 ml/h. The irreversible elimination clearance of MPL was 789 ml/h from liver and that of MPN was 2758 ml/h with liver accounting for 44%, lung 35%, and kidney 21%. The reversible metabolism elevated MPL exposure in rats by 13%. This highly complex PBPK model provided unique and comprehensive insights into the disposition of a major corticosteroid. SIGNIFICANCE STATEMENT: Our dual physiologically based pharmacokinetic (PBPK) study and model of methylprednisolone/methylprednisone (MPL/MPN) with multiple complexities reasonably characterized and parameterized their disposition, and provided greater insights into the interpretation of their pharmacodynamics in rats. Drug knowledge gained in this study may be translatable to higher-order species to appreciate the clinical utility of MPL. The complex model itself is instructive for advanced PBPK analysis of drugs with reversible metabolism and/or nonlinear tissue partitioning features.

Comparative Proteomic Analysis Identifies Key Metabolic Regulators of Gemcitabine Resistance in Pancreatic Cancer.

Pancreatic adenocarcinoma (PDAC) is highly refractory to treatment. Standard-of-care gemcitabine (Gem) provides only modest survival benefits, and development of Gem resistance (GemR) compromises its efficacy. Highly GemR clones of Gem-sensitive MIAPaCa-2 cells were developed to investigate the molecular mechanisms of GemR and implemented global quantitative differential proteomics analysis with a comprehensive, reproducible ion-current-based MS1 workflow to quantify ∼6000 proteins in all samples. In GemR clone MIA-GR8, cellular metabolism, proliferation, migration, and 'drug response' mechanisms were the predominant biological processes altered, consistent with cell phenotypic alterations in cell cycle and motility. S100 calcium binding protein A4 was the most downregulated protein, as were proteins associated with glycolytic and oxidative energy production. Both responses would reduce tumor proliferation. Upregulation of mesenchymal markers was prominent, and cellular invasiveness increased. Key enzymes in Gem metabolism pathways were altered such that intracellular utilization of Gem would decrease. Ribonucleoside-diphosphate reductase large subunit was the most elevated Gem metabolizing protein, supporting its critical role in GemR. Lower Ribonucleoside-diphosphate reductase large subunit expression is associated with better clinical outcomes in PDAC, and its downregulation paralleled reduced MIAPaCa-2 proliferation and migration and increased Gem sensitivity. Temporal protein-level Gem responses of MIAPaCa-2 versus GemR cell lines (intrinsically GemR PANC-1 and acquired GemR MIA-GR8) implicate adaptive changes in cellular response systems for cell proliferation and drug transport and metabolism, which reduce cytotoxic Gem metabolites, in DNA repair, and additional responses, as key contributors to the complexity of GemR in PDAC. These findings additionally suggest targetable therapeutic vulnerabilities for GemR PDAC patients.

Pharmacokinetic/Pharmacodynamic Modeling of Dexamethasone Anti-Inflammatory and Immunomodulatory Effects in LPS-Challenged Rats: A Model for Cytokine Release Syndrome.

Dexamethasone (DEX) is a potent synthetic glucocorticoid used for the treatment of variety of inflammatory and immune-mediated disorders. The RECOVERY clinical trial revealed benefits of DEX therapy in COVID-19 patients. Severe SARS-CoV-2 infection leads to an excessive inflammatory reaction commonly known as a cytokine release syndrome that is associated with activation of the toll like receptor 4 (TLR4) signaling pathway. The possible mechanism of action of DEX in the treatment of COVID-19 is related to its anti-inflammatory activity arising from inhibition of cytokine production but may be also attributed to its influence on immune cell trafficking and turnover. This study, by means of pharmacokinetic/pharmacodynamic modeling, aimed at the comprehensive quantitative assessment of DEX effects in lipopolysaccharide-challenged rats and to describe interrelations among relevant signaling molecules in this animal model of cytokine release syndrome induced by activation of TLR4 pathway. DEX was administered in a range of doses from 0.005 to 2.25 mg·kg-1 in LPS-challenged rats. Serum DEX, corticosterone (CST), tumor necrosis factor α, interleukin-6, and nitric oxide as well as lymphocyte and granulocyte counts in peripheral blood were quantified at different time points. A minimal physiologically based pharmacokinetic/pharmacodynamic (mPBPK/PD) model was proposed characterizing the time courses of plasma DEX and the investigated biomarkers. A high but not complete inhibition of production of inflammatory mediators and CST was produced in vivo by DEX. The mPBPK/PD model, upon translation to humans, may help to optimize DEX therapy in patients with diseases associated with excessive production of inflammatory mediators, such as COVID-19. SIGNIFICANCE STATEMENT: A mPBPK/PD model was developed to describe concentration-time profiles of plasma DEX, mediators of inflammation, and immune cell trafficking and turnover in LPS-challenged rats. Interrelations among DEX and relevant biomarkers were reflected in the mechanistic model structure. The mPBPK/PD model enabled quantitative assessment of in vivo potency of DEX and, upon translation to humans, may help optimize dosing regimens of DEX for the treatment of immune-related conditions associated with exaggerated immune response.

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.

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.

Transitioning from Basic toward Systems Pharmacodynamic Models: Lessons from Corticosteroids.

Technology in bioanalysis, -omics, and computation have evolved over the past half century to allow for comprehensive assessments of the molecular to whole body pharmacology of diverse corticosteroids. Such studies have advanced pharmacokinetic and pharmacodynamic (PK/PD) concepts and models that often generalize across various classes of drugs. These models encompass the "pillars" of pharmacology, namely PK and target drug exposure, the mass-law interactions of drugs with receptors/targets, and the consequent turnover and homeostatic control of genes, biomarkers, physiologic responses, and disease symptoms. Pharmacokinetic methodology utilizes noncompartmental, compartmental, reversible, physiologic [full physiologically based pharmacokinetic (PBPK) and minimal PBPK], and target-mediated drug disposition models using a growing array of pharmacometric considerations and software. Basic PK/PD models have emerged (simple direct, biophase, slow receptor binding, indirect response, irreversible, turnover with inactivation, and transduction models) that place emphasis on parsimony, are mechanistic in nature, and serve as highly useful "top-down" methods of quantitating the actions of diverse drugs. These are often components of more complex quantitative systems pharmacology (QSP) models that explain the array of responses to various drugs, including corticosteroids. Progressively deeper mechanistic appreciation of PBPK, drug-target interactions, and systems physiology from the molecular (genomic, proteomic, metabolomic) to cellular to whole body levels provides the foundation for enhanced PK/PD to comprehensive QSP models. Our research based on cell, animal, clinical, and theoretical studies with corticosteroids have provided ideas and quantitative methods that have broadly advanced the fields of PK/PD and QSP modeling and illustrates the transition toward a global, systems understanding of actions of diverse drugs. SIGNIFICANCE STATEMENT: Over the past half century, pharmacokinetics (PK) and pharmacokinetics/pharmacodynamics (PK/PD) have evolved to provide an array of mechanism-based models that help quantitate the disposition and actions of most drugs. We describe how many basic PK and PK/PD model components were identified and often applied to the diverse properties of corticosteroids (CS). The CS have complications in disposition and a wide array of simple receptor-to complex gene-mediated actions in multiple organs. Continued assessments of such complexities have offered opportunities to develop models ranging from simple PK to enhanced PK/PD to quantitative systems pharmacology (QSP) that help explain therapeutic and adverse CS effects. Concurrent development of state-of-the-art PK, PK/PD, and QSP models are described alongside experimental studies that revealed diverse CS actions.

Pharmacokinetic/Pharmacodynamic Assessment of Selective Phosphodiesterase Inhibitors in a Mouse Model of Autoimmune Hepatitis.

Autoimmune hepatitis (AIH) is a life-threatening disorder currently treated with nonspecific immunosuppressive drugs. It is postulated that phosphodiesterase (PDE) inhibitors, as agents exerting anti-inflammatory and immunomodulatory activities, may constitute a possible treatment of autoimmune disorders. This study develops a pharmacokinetic/pharmacodynamic (PK/PD) model to assess the effects of PDE-selective inhibitors, namely, cilostazol (PDE3), rolipram (PDE4), and BRL-50481 (PDE7), in a mouse model of AIH. The pharmacokinetics of the PDE inhibitors (PDEi) were assessed in male BALB/c mice after intraperitoneal administration. In pharmacodynamic studies, mice received PDEi and AIH was induced in these animals by intravenous injection of concanavalin A (ConA). Serum drug concentrations, tumor necrosis factor α (TNFα), interleukin 17 (IL-17), and aminotransferase activities were quantified. The PK/PD analysis was performed using ADAPT5 software. The PK/PD model assumes inhibition of cAMP hydrolysis in T cells by PDEi, ConA-triggered formation of TNFα and IL-17, suppression of TNFα and IL-17 production by cAMP, and stimulatory effects of TNFα and IL-17 on the hepatic release of aminotransferases. Selective blockage of PDE4 leads to the highest inhibition of cAMP degradation in T cells and amelioration of disease outcomes. However, inhibition of both PDE3 and PDE7 also contribute to this effect. The proposed PK/PD model may be used to assess and predict the activities of novel PDEi and their combinations in ConA-induced hepatitis. A balanced suppression of different types of PDE appears to be a promising treatment option for AIH; however, this hypothesis warrants testing in humans based on translation of the PK/PD model into clinical settings. SIGNIFICANCE STATEMENT: A novel PK/PD model of PDE inhibitor effects in mice with ConA-induced autoimmune hepatitis was developed involving a mechanistic component describing changes in cAMP concentrations in mouse T cells. According to model predictions, inhibition of PDE4 in T cells causes the highest cAMP elevation in T cells, but suppression of PDE3 and PDE7 also contribute to this effect. A balanced inhibition of PDE3, PDE4, and PDE7 appears to be a promising treatment strategy for AIH.

Assessing Liver-to-Plasma Partition Coefficients and In Silico Calculation Methods: When Does the Hepatic Model Matter in PBPK?

The primary models used in pharmacokinetics (PK) to assess hepatic clearance (CLh ) are the well-stirred (WSM), parallel tube (PTM), and dispersion model (DM) that differ in their internal flow patterns and assumed unbound liver concentrations. Physiologically-Based Pharmacokinetic (PBPK) models require a hepatic intrinsic clearance (CLint ) and tissue-to-plasma partition coefficient (Kp ). Given measured systemic and liver concentration-time profiles, these hepatic models perform similarly but yield model-specific CLint and Kp estimates. This work provides mathematical relationships for the three basic hepatic models and assesses their corresponding PBPK-relevant Kp values with literature-reported single-dose blood and liver concentration-time data of 14 compounds. Model fittings were performed with an open-loop approach where the CLh and extraction ratio (ER) were first estimated from fitting the blood data yielding CLint values for the three hepatic models. The pre-fitted blood data served as forcing input functions to obtain PBPK-operative Kp estimates that were compared with those obtained by the tissue/plasma area ratio (AR), Chen & Gross (C&G) and published in silico methods. The CLint and Kp values for the hepatic models increased with the ER and both showed a rank order being WSM > DM > PTM. Drugs with low ER showed no differences as expected. With model-specific CLint and Kp values, all hepatic models predict the same steady-state Kp (Kp ss ) that is comparable to those from the AR and C&G methods and reported by direct measurement. All in silico methods performed poorly for most compounds. Hepatic model selection requires cautious application and interpretation in PBPK modeling. Significance Statement The three hepatic models generate different single-dose (non-steady-state) values of CLint and Kp in PBPK models especially for drugs with high ER; however, all Kp ss values expected from constant rate infusion studies were the same. These findings are relevant when using these models for IVIVE where a model-dependent CLint is used to correct measured tissue concentrations for depletion by metabolism. This model-dependency may also have an impact when assessing the PK/pharmacodynamic relationships when effects relate to assumed hepatic concentrations.

Consideration of Fractional Distribution Parameter fd in the Chen and Gross Method for Tissue-to-Plasma Partition Coefficients: Comparison of Several Methods.

PURPOSE: The tissue-to-plasma partition coefficient (Kp) describes the extent of tissue distribution in physiologically-based pharmacokinetic (PBPK) models. Constant-rate infusion studies are common for experimental determination of the steady-state Kp,ss, while the tissue-plasma concentration ratio (CT/Cp) in the terminal phase after intravenous doses is often utilized. The Chen and Gross (C&G) method converts a terminal slope CT/Cp to Kp,ss based on assumptions of perfusion-limited distribution in tissue-plasma equilibration. However, considering blood flow (QT) and apparent tissue permeability (fupPSin) in the rate of tissue distribution, this report extends the C&G method by utilizing a fractional distribution parameter (fd). METHODS: Relevant PBPK equations for non-eliminating and eliminating organs along with lung and liver were derived for the conversion of CT/Cp values to Kp,ss. The relationships were demonstrated in rats with measured CT/Cp and Kp,ss values and the model-dependent fd for 8 compounds with a range of permeability coefficients. Several methods of assessing Kp were compared. RESULTS: Utilizing fd in an extended C&G method, our estimations of Kp,ss from CT/Cp were improved, particularly for lower permeability compounds. However, four in silico methods for estimating Kp performed poorly across tissues in comparison with measured Kp values. Mathematical relationships between Kp and Kp,ss that are generally applicable for eliminating organs with tissue permeability limitations necessitates inclusion of an extraction ratio (ER) and fd. CONCLUSION: Since many different types/sources of Kp are present in the literature and used in PBPK models, these perspectives and equations should provide better insights in measuring and interpreting Kp values in PBPK.

Pharmacodynamic model of slow reversible binding and its applications in pharmacokinetic/pharmacodynamic modeling: review and tutorial.

Therapeutic responses of most drugs are initiated by the rate and degree of binding to their receptors or targets. The law of mass action describes the rate of drug-receptor complex association (kon) and dissociation (koff) where the ratio koff/kon is the equilibrium dissociation constant (Kd). Drugs with slow reversible binding (SRB) often demonstrate delayed onset and prolonged pharmacodynamic effects. This report reviews evidence for drugs with SRB features, describes previous pharmacokinetic/pharmacodynamic (PK/PD) modeling efforts of several such drugs, provides a tutorial on the mathematics and properties of SRB models, demonstrates applications of SRB models to additional compounds, and compares PK/PD fittings of SRB with other mechanistic models. We identified and summarized 52 drugs with in vitro-confirmed SRB from a PubMed literature search. Simulations with a SRB model and observed PK/PD profiles showed delayed and prolonged responses and that increasing doses/kon or decreasing koff led to greater expected maximum effects and a longer duration of effects. Recession slopes for return of responses to baseline after single doses were nearly linear with an inflection point that approaches a limiting value at larger doses. The SRB model newly captured literature data for the antihypertensive effects of candesartan and antiallergic effects of noberastine. Their PD profiles could also be fitted with indirect response and biophase models with minimal differences. The applicability of SRB models is probably commonplace, but underappreciated, owing to the need for in vitro confirmation of binding kinetics and the similarity of PK/PD profiles to models with other mechanistic determinants.

Physiologically Based Pharmacokinetics of Lysosomotropic Chloroquine in Rat and Human.

A semimechanistic physiologically based pharmacokinetic (PBPK) model for chloroquine (CQ), a highly lysosomotropic weak base, was applied to digitized rat and human concentration versus time data. The PBPK model in rat featured plasma and red blood cell (RBC) concentrations, extensive and apparent nonlinear tissue distribution, fitted hepatic and renal intrinsic clearances, and a plasma half-life of about 1 day. Tissue-to-plasma CQ ratios at 50 hours after dosing were highest in lung, kidney, liver, and spleen (182-318) and lower in heart, muscle, brain, eye, and skin (11-66). The RBC-to-plasma ratio of 11.6 was assumed to reflect cell lipid partitioning. A lysosome-based extended model was used to calculate subcellular CQ concentrations based on tissue mass balances, fitted plasma, interstitial and free cytosol concentrations, and literature-based pH and pKa values. The CQ tissue component concentrations ranked as follows: lysosome > > acidic phospholipid > plasma = interstitial = cytosol ≥ neutral lipids. The extensive lysosome sequestration appeared to change over time and was attributed to lowering pH values caused by proton pump influx of hydrogen ions. The human-to-rat volume of distribution (Vss) ratio of 7 used to scale rat tissue partitioning to human along with estimation of hepatic clearance allowed excellent fitting of oral-dose PK (150-600 mg) of CQ with a 50-day half-life in humans. The prolonged PK of chloroquine was well characterized for rat and human with this PBPK model. The calculated intratissue concentrations and lysosomal effects have therapeutic relevance for CQ and other cationic drugs. SIGNIFICANCE STATEMENT: Sequestration in lysosomes is a major factor controlling the pharmacokinetics and pharmacology of chloroquine and other cationic drugs. This report provides comprehensive physiologic modeling of chloroquine distribution in tissues and overall disposition in rat and human that reveals expected complexities and inferences related to its subcellular association with various tissue components.

Synergistic Pharmacodynamic Effects of Gemcitabine and Fibroblast Growth Factor Receptor Inhibitors on Pancreatic Cancer Cell Cycle Kinetics and Proliferation.

Median survival of pancreatic ductal adenocarcinoma cancer (PDAC) is 6 months, with 9% 5-year survival. Standard-of-care gemcitabine (Gem) provides only modest survival benefits, and combination therapies integrating novel targeted agents could improve outcomes. Fibroblast growth factor (FGF) receptors (FGFRs) play important roles in PDAC growth and invasion. Therefore, FGFR inhibitors (FGFRi) merit further investigation. Efficacy of Gem combined with NVP-BGJ398, a pan-FGFRi, was investigated in multiple PDAC cell lines exposed to the drugs alone and combined. Cell cycle distribution and cell numbers were quantified over time. Two pharmacodynamic models were developed to investigate Gem/BGJ398 interactions quantitatively: a drug-mediated cell proliferation/death model, and a drug-perturbed cell cycle progression model. The models captured temporal changes in cell numbers, cell cycle progression, and cell death during drug exposure. Simultaneous fitting of all data provided reasonable parameter estimates. Therapeutic efficacy was then evaluated in a PDAC mouse model. Compared with Gem alone, combined Gem + FGFRi significantly downregulated ribonucleotide-diphosphate reductase large subunit 1 (RRM1), a gemcitabine resistance (GemR) biomarker, suggesting the FGFRi inhibited GemR emergence. The cell proliferation/death pharmacodynamic model estimated the drug interaction coefficient ψ death = 0.798, suggesting synergistic effects. The mechanism-based cell cycle progression model estimated drug interaction coefficient ψ cycle = 0.647, also suggesting synergy. Thus, FGFR inhibition appears to synergize with Gem in PDAC cells and tumors by sensitizing cells to Gem-mediated inhibition of proliferation and cell cycle progression. SIGNIFICANCE STATEMENT: An integrated approach of quantitative modeling and experimentation was employed to investigate the nature of fibroblast growth factor receptor inhibitor (FGFRi)/gemcitabine (Gem) interaction, and to identify mechanisms by which FGFRi exposure reverses Gem resistance in pancreatic cancer cells. The results show that FGFRi interacts synergistically with Gem to sensitize pancreatic cancer cells and tumors to Gem-mediated inhibition of proliferation and cell cycle progression. Thus, addition of FGFRi to standard-of-care Gem treatment could be a clinically deployable approach to enhance therapeutic benefit to pancreatic cancer patients.

Mathematical modeling of mammalian circadian clocks affecting drug and disease responses.

To align with daily environmental changes, most physiological processes in mammals exhibit a time-of-day rhythmicity. This circadian control of physiology is intrinsically driven by a cell-autonomous clock gene network present in almost all cells of the body that drives rhythmic expression of genes that regulate numerous molecular and cellular processes. Accordingly, many aspects of pharmacology and toxicology also oscillate in a time-of-day manner giving rise to diverse effects on pharmacokinetics and pharmacodynamics. Genome-wide studies and mathematical modeling are available tools that have significantly improved our understanding of these nonlinear aspects of physiology and therapeutics. In this manuscript current literature and our prior work on the model-based approaches that have been used to explore circadian genomic systems of mammals are reviewed. Such basic understanding and having an integrative approach may provide new strategies for chronotherapeutic drug treatments and yield new insights for the restoration of the circadian system when altered by diseases.

Population pharmacokinetic modeling of intramuscular and oral dexamethasone and betamethasone in Indian women.

Population analysis of pharmacokinetic data for five differing dosage forms and routes for dexamethasone and betamethasone in 48 healthy nonpregnant Indian women was performed that accounted for a partial and complex cross-over design. Single doses of 6 mg dexamethasone phosphate (DEX-P), betamethasone phosphate (BET-P), or 1:1 mixture of betamethasone phosphate and acetate (BET-PA) were administered orally (PO) or intramuscularly (IM). Plasma concentrations collected for two periods over 96 h were described with a two-compartment model with differing PO and IM first-order absorption inputs. Clearances and volumes were divided by the IM bioavailability [Formula: see text]. The homogeneous ages, body weights, and ethnicity of the women obviated covariate analysis. Parameter estimates were obtained by the Laplace estimation method implemented in NONMEM 7.4. Typical values for dexamethasone were clearance ([Formula: see text] of 9.29 L/h, steady-state volume ([Formula: see text] of 56.4 L, IM absorption constant [Formula: see text] of 0.460 1/h and oral absorption constant ([Formula: see text] of 0.936 1/h. Betamethasone parameters were CL/FIM of 5.95 L/h, [Formula: see text] of 72.4 L, [Formula: see text] of 0.971 1/h, and [Formula: see text] of 1.21 1/h. The PO to IM F values were close to 1.0 for both drugs. The terminal half-lives averaged about 7.5 h for DEX, 17 h for BET, and 78 h for BET from BET-PA with the latter reflecting very slow release of BET from the acetate ester. Overall, BET exhibited slower clearance, larger volume of distribution, faster absorption, and longer persistence than DEX. These data may be useful in considering exposures when substituting one form of corticosteroid for another.

Pathway-level analysis of genome-wide circadian dynamics in diverse tissues in rat and mouse.

A computational framework is developed to enable the characterization of genome-wide, multi-tissue circadian dynamics at the level of "functional groupings of genes" defined in the context of signaling, cellular/genetic processing and metabolic pathways in rat and mouse. Our aim is to identify how individual genes come together to generate orchestrated rhythmic patterns and how these may vary within and across tissues. We focus our analysis on four tissues (adipose, liver, lung, and muscle). A genome-wide pathway-centric analysis enables us to develop a comprehensive picture on how the observed circadian variation at the individual gene level, orchestrates functional responses at the pathway level. Such pathway-based "meta-data" analysis enables the rational integration and comparison across platforms and/or experimental designs evaluating emergent dynamics, as opposed to comparisons of individual elements. One of our key findings is that when considering the dynamics at the pathway level, a complex behavior emerges. Our work proposes that tissues tend to coordinate gene's circadian expression in a way that optimizes tissue-specific pathway activity, depending of each tissue's broader role in homeostasis.

Population pharmacodynamic modeling of intramuscular and oral dexamethasone and betamethasone effects on six biomarkers with circadian complexities in Indian women.

Population pharmacokinetic/pharmacodynamic (PK/PD) analysis was performed for extensive data for differing dosage forms and routes for dexamethasone (DEX) and betamethasone (BET) in 48 healthy nonpregnant Indian women in a partial and complex cross-over design. Single doses of 6 mg dexamethasone phosphate (DEX-P), betamethasone phosphate (BET-P), or 1:1 mixture of betamethasone phosphate and acetate (BET-PA) were administered orally (PO) or intramuscularly (IM) where each woman enrolled in a two-period cross-over study. Plasma concentrations collected over 96 h were described with a two-compartment model with differing PO and IM first-order absorption inputs. Overall, BET exhibited slower clearance, similar volume of distribution, faster absorption, and longer persistence than DEX with BET acetate producing extremely slow absorption but full bioavailability of BET. Six biomarkers were assessed over a 24-h baseline period with four showing circadian rhythms with complex baselines. These baselines and the strong responses seen after drug dosing were fitted with various indirect response models using the Laplace estimation methods in NONMEM 7.4. Both the PK and six biomarker responses were well-described with modest variability likely due to the homogeneous ages, weights, and ethnicities of the women. The drugs either inhibited or stimulated the influx processes with some models requiring joint inclusion of drug effects on circadian cortisol suppression. The biomarkers and order of sensitivity (lowest IC50/SC50 to highest) were: cortisol, T-helper cells, basophils, glucose, neutrophils, and T-cytotoxic cells. DEX sensitivities were generally greater than BET with corresponding mean ratios for these biomarkers of 2.86, 1.27, 1.72, 1.27, 2.69, and 1.06. Overall, the longer PK (e.g. half-life) of BET, but lesser PD activity (e.g. higher IC50), produces single-dose response profiles that appear quite similar, except for the extended effects from BET-PA. This comprehensive population modeling effort provides the first detailed comparison of the PK profiles and six biomarker responses of five commonly used dosage forms of DEX and BET in healthy women.

Across-species meta-analysis of dexamethasone pharmacokinetics utilizing allometric and scaling modeling approaches.

The pharmacokinetic (PK) parameters of dexamethasone (DEX) in 11 species were collected from the literature and clearances (CL) assessed by basic allometric methods, and concentration-time course profiles were fitted using two PK models incorporating physiological or allometric scaling. Plots of log CL vs. log body weights (BW) correlated reasonably with R2  = 0.91, with a maximum ratio of actual to fitted CL of 6 (for pig). A minimal physiologically-based pharmacokinetic (mPBPK) model containing blood and two lumped tissue compartments and integrated utilization of physiological parameters was compared to an allometric two-compartment model (a2CM). The plasma PK profiles of DEX from 11 species were analyzed jointly, with the mPBPK model having conserved partition coefficients (Kp ), physiologic blood and tissue volumes, and species-specific CL values. The DEX PK profiles were reasonably captured by the mPBPK model for 9 of 11 species in the joint analysis with three fitted parameters (besides CL) including an overall tissue-to-plasma partition coefficient of 1.07. The a2CM with distribution CL and central and peripheral volumes scaled allometrically fitted the plasma concentration profiles similarly but required a total of six parameters (besides CL). Overall, the literature reported that DEX CL values exhibit moderate variability (mean = 0.64 L/h/kg; coefficient of variation = 105%), but distribution parameters were largely conserved across most species.

Physiologically Based Pharmacokinetic Modeling Involving Nonlinear Plasma and Tissue Binding: Application to Prednisolone and Prednisone in Rats.

The pharmacokinetics (PK) of prednisolone (PNL) exhibit nonlinearity related to plasma protein binding, tissue binding, metabolic interconversion with prednisone (PN), and renal elimination. Blood and 11 tissues were collected from male Wistar rats after steady-state (SS) infusion and after subcutaneous boluses of 50 mg/kg of PNL. Concentrations of PNL and PN were measured by liquid chromatography-tandem mass spectrometry. Plasma and tissue profiles were described using a complex physiologically based pharmacokinetics (PBPK) model. Concentrations of PN and PNL were in rapid equilibrium in plasma and tissues. The tissue partition coefficients (K p ) of PNL calculated from most subcutaneously dosed tissue and plasma areas were similar to SS infusion and in silico values. The blood-to-plasma ratio of PNL was 0.71 with similar red blood cell and unbound-plasma concentrations. Plasma protein binding (60%-90%) was related to corticosteroid-binding globulin (CBG) saturation. Tissue distribution was nonlinear. The equilibrium dissociation constant (K d ) of PNL shared by all tissues was 3.01 ng/ml, with the highest binding in muscle, followed by liver, heart, intestine, and bone and the lowest binding in skin, spleen, fat, kidney, lung, and brain. Fat and bone distribution assumed access only to interstitial space. Brain PNL concentrations (K p = 0.05) were low owing to presumed P-glycoprotein-mediated efflux. Clearances of CBG-free PNL were 1789 from liver and 191.2 ml/h from kidney. The PN/PNL ratio was nonlinear for plasma, spleen, heart, intestine, bone, fat, and linear for the remaining tissues. Our PBPK model with multiple complexities well described the PK profiles of PNL and PN in blood, plasma, and diverse tissues. SIGNIFICANCE STATEMENT: Because steroids, such as prednisolone and prednisone, have similar and complex pharmacokinetics properties in various species, receptors in most tissues, and multiple therapeutic and adverse actions, this physiologically based pharmacokinetics (PBPK) model may provide greater insights into the pharmacodynamic complexities of corticosteroids. The complex properties of these compounds require innovative PBPK modeling approaches that may be instructive for other therapeutic agents.

Seeking Nonspecific Binding: Assessing the Reliability of Tissue Dilutions for Calculating Fraction Unbound.

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.

Physiologically Based Pharmacokinetics of Dexamethasone in Rats.

Blood and multitissue concentration-time profiles for dexamethasone (DEX), a synthetic corticosteroid, were measured in male rats after subcutaneous bolus and infusion dosing. A physiologically based pharmacokinetics (PBPK) model was applied for 12 measured tissues. Tissue partition coefficients (K p ) and metabolic clearance were assessed from infusion studies. Blood cell to plasma partitioning (0.664) and plasma free fraction (0.175) for DEX were found to be moderate. DEX was extensively partitioned into liver (K p = 6.76), whereas the calculated K p values of most tissues ranged between 0.1 and 1.5. Despite the moderate lipophilicity of DEX (log P = 1.8), adipose exhibited very limited distribution (K p = 0.17). Presumably due to P-glycoprotein-mediated efflux, DEX concentrations were very low in brain compared with its expected high permeability. Infusion studies yielded K p values from male and female rats at steady state that were similar. In silico K p values calculated for different tissues by using GastroPlus software were similar to in vivo values except for adipose and liver. Glucocorticoid receptors are found in diverse tissues, and these PBPK modeling results may help provide exposure profiles driving pharmacodynamic effects of DEX. SIGNIFICANCE STATEMENT: Our physiologically based pharmacokinetics model describes the experimentally determined tissue and plasma dexamethasone (DEX) pharmacokinetics (PK) profiles in rats reasonably well. This model can serve for further investigation of DEX tissue distribution in rats as the PK driving force for PD effects in different tissues. No major sex differences were found for DEX tissue distribution. Knowledge gained in this study may be translatable to higher-order species including humans.