PBPK modeling of intestinal and liver enzymes and transporters in drug absorption and sequential metabolism

Experimental strategies have long been applied for in vitro or in vivo evaluation of the effect of transporters and/or enzymes on the bioavailability. However, the lack of specific inhibitors or inducers of transporters and enzymes and the multiplicity of nuclear receptors in gene regulation and cross talk have led to compromised assessment of these effects in vivo. These and other causes have resulted in confusion and controversy in transporter-enzyme interplay. In this review, physiologically-based pharmacokinetic (PBPK) intestinal and liver models are utilized to predict the contributions of enzymes and transporters on intestinal availability (F(I)) and hepatic availability (F(H)), with the aim to fully understand the impact of these variables on bioavailability (F(sys)) in vivo. We emphasize the often overlooked impact of influx and efflux clearances, and apply the PBPK models and their solutions to examine individual organ clearances of the intestine and the liver. In order to accurately predict oral bioavailability, these organ models are incorporated into the whole body PBPK model, and additional complicated scenarios such as segmental differences and zonal heterogeneity of transporters and enzymes in the intestine and liver and segregated blood flow patterns of the intestine are further discussed. The sequential metabolism of a drug to form primary and secondary metabolites in the first-pass organs is considered in PBPK modeling, revealing that the segregated flow model (SFM) of the intestine is more appropriate than the traditional PBPK intestinal model (TM). Examples are included to highlight the potential application of these PBPK models on the quantitative prediction of bioavailability.     

Physiologically-Based PK/PD Modelling of Therapeutic Macromolecules

Therapeutic proteins are a diverse class of drugs consisting of naturally occurring or modified proteins, and due to their size and physico-chemical properties, they can pose challenges for the pharmacokinetic and pharmacodynamic studies. Physiologically-based pharmacokinetics (PBPK) modelling has been effective for early in silico prediction of pharmacokinetic properties of new drugs. The aim of the present workshop was to discuss the feasibility of PBPK modelling of macromolecules. The classical PBPK approach was discussed with a presentation of the successful example of PBPK modelling of cyclosporine A. PBPK model was performed with transport of the cyclosporine across cell membranes, affinity to plasma proteins and active membrane transporters included to describe drug transport between physiological compartments. For macromolecules, complex PBPK modelling or permeability-limited and/or target-mediated distribution was discussed. It was generally agreed that PBPK modelling was feasible and desirable. The role of the lymphatic system should be considered when absorption after extravascular administration is modelled. Target-mediated drug disposition was regarded as an important feature for generation of PK models. Complex PK-models may not be necessary when a limited number of organs are affected. More mechanistic PK/PD models will be relevant when adverse events/toxicity are included in the PK/PD modelling.     

An integrated approach to model hepatic drug clearance.

It has been well accepted that hepatic drug extraction depends on the blood flow, vascular binding, transmembrane barriers, transporters, enzymes and cosubstrate and their zonal heterogeneity. Models of hepatic drug clearances have been appraised with respect to their utility in predicting drug removal by the liver. Among these models, the "well-stirred" model is the simplest since it assumes venous equilibration, with drug emerging from the outflow being in equilibrium with drug within the liver, and the concentration is the same throughout. The "parallel tube" and dispersion models, and distributed model of Goresky and co-workers have been used to account for the observed sinusoidal concentration gradient from the inlet and outlet. Departure from these models exists to include heterogeneity in flow, enzymes, and transporters. This article utilized the physiologically based pharmacokinetic (PBPK) liver model and its extension that include heterogeneity in enzymes and transporters to illustrate how in vitro uptake and metabolic data from zonal hepatocytes on transport and enzymes may be used to predict the kinetics of removal in the intact liver; binding data were also necessary. In doing so, an integrative platform was provided to examine determinants of hepatic drug clearance. We used enalapril and digoxin as examples, and described a simple liver PBPK model that included transmembrane transport and metabolism occurring behind the membrane, and a zonal model in which the PBPK model was expanded three sets of sub-compartments that are arranged sequentially to represent zones 1, 2, and 3 along the flow path. The latter model readily accommodated the heterogeneous distribution of hepatic enzymes and transporters. Transport and metabolic data, piecewise information that served as initial estimates, allowed for the unknown efflux and other intrinsic clearances to be estimated. The simple or zonal PBPK model provides predictive views on the hepatic removal of drugs and metabolites.     

Impact of genetic polymorphisms of transporters on the pharmacokinetic, pharmacodynamic and toxicological properties of anionic drugs

As the importance of drug transporters in the clinical pharmacokinetics of drugs is recognized, genetic polymorphisms of drug transporters have emerged as one of the determinant factors to produce the inter-individual variability of pharmacokinetics. Many clinical studies have shown the influence of genetic polymorphisms of drug transporters on the pharmacokinetics and subsequent pharmacological and toxicological effects of drugs. The functional change in a transporter in clearance organs such as liver and kidney affects the drug concentration in the blood circulation, while that in the pharmacological or toxicological target can alter the local concentration at the target sites without changing its plasma concentration. As for the transporters for organic anions, some single nucleotide polymorphisms (SNPs) or haplotypes occurring with high frequency in organic anion transporting polypeptide (OATP) 1B1, multidrug resistance 1 (MDR1), and breast cancer resistance protein (BCRP) have been extensively investigated in both human clinical studies and in vitro functional assays. We introduce some examples showing the relationship between haplotypes in transporters and pharmacokinetics and pharmacological effects of drugs. We also discuss how to predict the effect of functional changes in drug transporters caused by genetic polymorphisms on the pharmacokinetics of drugs from in vitro data.     

慢性肾功能不全对药动学影响研究进展

慢性肾功能不全是临床常见疾病之一,主要经肾消除及主要经非肾途径消除的药物,药代动力学及药效学均受慢性肾功能不全疾病的影响。肾功能不全尤其是终末期肾病不仅直接减少药物的排泄,还能影响药物吸收、药物转运及药物代谢,从而影响药物的非肾消除。肾功能不全可直接抑制药物代谢酶及转运体的活性并减少酶及转运体在胃肠道、肝脏等多个器官的mRNA及蛋白的表达。引起这种抑制作用的原因尚不明确,血液透析可以部分抵消或逆转这种抑制作用并存在透析消除。     

Expression of thirty-six drug transporter genes in human intestine, liver, kidney, and organotypic cell lines.

This study was designed to quantitatively assess the mRNA expression of 36 important drug transporters in human jejunum, colon, liver, and kidney. Expression of these transporters in human organs was compared with expression in commonly used cell lines (Caco-2, HepG2, and Caki-1) originating from these organs to assess their value as in vitro transporter system models, and was also compared with data obtained from the literature on expression in rat tissues to assess species differences. Transporters that were highly expressed in the intestine included HPT1, PEPT1, BCRP, MRP2, and MDR1, whereas, in the liver, OCT1, MRP2, OATP-C, NTCP and BSEP were the main transporters. In the kidney, OAT1 was expressed at the highest levels, followed by OAT3, OAT4, MCT5, MDR1, MRP2, OCT2, and OCTN2. The best agreement between human tissue and the representative cell line was observed for human jejunum and Caco-2 cells. Expression in liver and kidney ortholog cell lines was not correlated with that in the associated tissue. Comparisons with rat transporter gene expression revealed significant species differences. Our results allowed a comprehensive quantitative comparison of drug transporter expression in human intestine, liver, and kidney. We suggest that it would be beneficial for predictive pharmacokinetic research to focus on the most highly expressed transporters. We hope that our comparison of rat and human tissue will help to explain the observed species differences in in vivo models, increase understanding of the impact of active transport processes on pharmacokinetics and distribution, and improve the quality of predictions from animal studies to humans.     

A Simplified PBPK Modeling Approach for Prediction of Pharmacokinetics of Four Primarily Renally Excreted and CYP3A Metabolized Compounds During Pregnancy

During pregnancy, a drug’s pharmacokinetics may be altered and hence anticipation of potential systemic exposure changes is highly desirable. Physiologically based pharmacokinetics (PBPK) models have recently been used to influence clinical trial design or to facilitate regulatory interactions. Ideally, whole-body PBPK models can be used to predict a drug’s systemic exposure in pregnant women based on major physiological changes which can impact drug clearance (i.e., in the kidney and liver) and distribution (i.e., adipose and fetoplacental unit). We described a simple and readily implementable multitissue/organ whole-body PBPK model with key pregnancy-related physiological parameters to characterize the PK of reference drugs (metformin, digoxin, midazolam, and emtricitabine) in pregnant women compared with the PK in nonpregnant or postpartum (PP) women. Physiological data related to changes in maternal body weight, tissue volume, cardiac output, renal function, blood flows, and cytochrome P450 activity were collected from the literature and incorporated into the structural PBPK model that describes HV or PP women PK data. Subsequently, the changes in exposure (area under the curve (AUC) and maximum concentration (C _max)) in pregnant women were simulated. Model-simulated PK profiles were overall in agreement with observed data. The prediction fold error for C _max and AUC ratio (pregnant vs. nonpregnant) was less than 1.3-fold, indicating that the pregnant PBPK model is useful. The utilization of this simplified model in drug development may aid in designing clinical studies to identify potential exposure changes in pregnant women a priori for compounds which are mainly eliminated renally or metabolized by CYP3A4.     

Pharmacokinetic interplay of phase II metabolism and transport: a theoretical study

Understanding of the interdependence of cytochrome P450 enzymes and P-glycoprotein in disposition of drugs (also termed "transport-metabolism interplay") has been significantly advanced in recent years. However, whether such "interplay" exists between phase II metabolic enzymes and efflux transporters remains largely unknown. The objective of this article is to explore the role of efflux transporters (acting on the phase II metabolites) in disposition of the parent drug in Caco-2 cells, liver, and intestine via simulations utilizing a catenary model (for Caco-2 system) and physiologically based pharmacokinetic (PBPK) models (for the liver and intestine). In all three models, "transport-metabolism interplay" (i.e., inhibition of metabolite efflux decreases the metabolism) can be observed only when futile recycling (or deconjugation) occurred. Futile recycling appeared to bridge the two processes (i.e., metabolite formation and excretion) and enable the interplay thereof. Without futile recycling, metabolite formation was independent on its downstream process excretion, thus impact of metabolite excretion on its formation was impossible. Moreover, in liver PBPK model with futile recycling, impact of biliary metabolite excretion on the exposure of parent drug [(systemic (reservoir) area under the concentration-time curve (AUC(R1))] was limited; a complete inhibition of efflux resulted in AUC(R1) increases of less than 1-fold only. In intestine PBPK model with futile recycling, even though a complete inhibition of efflux could result in large elevations (e.g., 3.5-6.0-fold) in AUC(R1), an incomplete inhibition of efflux (e.g., with a residual activity of ≥ 20% metabolic clearance) saw negligible increases (<0.9-fold) in AUC(R1). In conclusion, this study presented mechanistic observations of pharmacokinetic interplay between phase II enzymes and efflux transporters. Those studying such "interplay" are encouraged to adequately consider potential consequences of inhibition of efflux transporters in humans.     

In vitro ADME phenotyping in drug discovery: current challenges and future solutions

Drug-metabolizing enzymes and drug transporters are key regulators of drug disposition and pharmacodynamics, which are closely linked to drug efficacy and safety. In this article, current challenges and future solutions to predicting their influence on pharmacokinetics and inter-organ distribution in humans, from data generated during the drug discovery decision-making process, are presented. In vitro phenotyping strategies for drug metabolizing enzymes (eg, CYP3A4, UGT1A1) and transporters (eg, OATP1B1) are offered, including perspectives on a selection of in vitro systems, novel in vitro phenotyping reagents and remaining technology gaps, challenges in extrapolating in vitro data to the in vivo situation, in silico models for the prediction of whether compounds are enzyme or transporter substrates, and the impact of pharmacogenomics.     

肾脏转运体和/或代谢酶介导药物排泄及药物相互作用的研究进展

肾脏是人体最重要的排泄器官。肾单元近端小管细胞具有多种药物转运体和代谢酶,在药物及其代谢物处置中发挥关键作用。近端小管细胞中主要转运体包括有机阴离子转运体、有机阳离子转运体、有机阳离子/肉毒碱转运体、多药及毒素外排转运蛋白、P-糖蛋白、乳腺癌耐药蛋白和多药耐药相关蛋白;主要代谢酶包括细胞色素P450酶,UDP-葡萄糖醛酸基转移酶、磺酸基转移酶、谷胱甘肽S-转移酶。肾脏转运体和/或代谢酶介导药物相互作用(DDIs)是临床关注的重要问题。肾脏转运体和代谢酶存在密切协作关系,在肾脏也存在多种相互作用现象(包括转运-转运相互作用,代谢-代谢相互作用和转运-代谢相互作用),其显著影响药物肾脏处置、临床疗效和肾毒性。本文系统阐述了这些相互作用对药物及其代谢物的肾脏排泄、药动学、DDIs和肾毒性的影响。今后需要进一步阐明肾脏转运-代谢相互作用机制,将有助于研究体内药物肾脏处置和DDIs,促进临床合理用药。     

The use of PBPK modeling across the pediatric age range using propofol as a case

The project SAFEPEDRUG aims to provide guidelines for drug research in children, based on bottom-up and top-down approaches. Propofol, one of the studied model compounds, was selected because it is extensively metabolized in liver and kidney, with an important role for the glucuronidation pathway. Besides, being a lipophilic molecule, it is distributed into fat tissues, from where it redistributes into the systemic circulation. In the past, both bottom-up (Physiologically based pharmacokinetic, PBPK) and top-down approaches (population pharmacokinetic, popPK) were applied to describe its pharmacokinetics (PK). In this work, a combination of the two was used to check their performance to describe PK in children and neonates (both term and preterm) using propofol as a case compound. First, in vitro data was generated in human liver microsomes and recombinant enzymes and used to develop an adult PBPK model in Simcyp^®. Activity adjustment factors (AAFs) were calculated to account for differences between in vitro and in vivo enzyme activity. Clinical data were analyzed using a 3-compartment model in NONMEM. These data were used to construct a retrograde PBPK model and for qualification of the PBPK models. Once an accurate in vivo clearance was obtained accounting for the contribution of the different metabolic pathways, the resulting PBPK models were challenged with new data for qualification. After that, the constructed adult PPBK model for propofol was extrapolated to the pediatric population. Both the default built-in and in vivo derived ontogeny functions were used to do so. The models were qualified by comparing their predicted PK parameters to published values, and by comparison of predicted concentration–time profiles to available clinical data. Clearance values were predicted well, especially when compared with values obtained from trials where long-term sampling was applied, whereas volume of distribution was lower compared to the most common popPK model predictions. Concentration–time profiles were predicted well up until and including the preterm neonatal population. In this work, it was thus shown that PBPK can be used to predict the PK up to and including the preterm neonatal population without the use of pediatric in vivo data. This work adds weight to the need for further development of PBPK models, especially regarding distribution modeling and the use of in vivo derived ontogeny functions.     

中西药相互作用的药代动力学机制研究进展

近年来,中西药联合应用日益增多,引起诸多药物相互作用和不良反应。药物代谢酶和转运体在药物的体内代谢和处置过程中发挥重要的作用,中西药联用不仅会改变药物的理化性质和人体的生理特性,而且可以调节药物代谢酶和转运体的表达和功能,进而产生药动学相互作用,引起药物疗效的改变。本文对中西药相互作用的药动学基础及其评价方法进行总结,旨在为中西药的联合应用和研究提供参考。     

Physiologically based pharmacokinetic (PBPK) modelling tools: how to fit with our needs?

In 2005, a survey compared a number of commercial PBPK software available at the time, mainly focusing on 'ready to use' modelling tools. Since then, these tools and software have been further developed and improved to allow modellers to perform WB-PBPK modelling including ADME processes at a high level of sophistication. This review presents a comparison of the features, values and limitations of both the 'ready to use' software and of the traditional user customizable software that are frequently used for the building and use of PBPK models, as well as the challenges associated with the various modelling approaches regarding their current and future use. PBPK models continue to be used more and more frequently during the drug development process since they represent a quantitative, physiologically realistic platform with which to simulate and predict the impact of various potential scenarios on the pharmacokinetics and pharmacodynamics of drugs. The 'ready to use' PBPK software has been a major factor in the increasing use of PBPK modelling in the pharmaceutical industry, opening up the PBPK approach to a broader range of users. The challenge is now to educate and to train scientists and modellers to ensure their appropriate understanding of the assumptions and the limitations linked both to the physiological framework of the 'virtual body' and to the scaling methodology from in vitro to in vivo (IVIVE).     

Pharmacometric Models for Characterizing the Pharmacokinetics of Orally Inhaled Drugs

During the last decades, the importance of modeling and simulation in clinical drug development, with the goal to qualitatively and quantitatively assess and understand mechanisms of pharmacokinetic processes, has strongly increased. However, this increase could not equally be observed for orally inhaled drugs. The objectives of this review are to understand the reasons for this gap and to demonstrate the opportunities that mathematical modeling of pharmacokinetics of orally inhaled drugs offers. To achieve these objectives, this review (i) discusses pulmonary physiological processes and their impact on the pharmacokinetics after drug inhalation, (ii) provides a comprehensive overview of published pharmacokinetic models, (iii) categorizes these models into physiologically based pharmacokinetic (PBPK) and (clinical data-derived) empirical models, (iv) explores both their (mechanistic) plausibility, and (v) addresses critical aspects of different pharmacometric approaches pertinent for drug inhalation. In summary, pulmonary deposition, dissolution, and absorption are highly complex processes and may represent the major challenge for modeling and simulation of PK after oral drug inhalation. Challenges in relating systemic pharmacokinetics with pulmonary efficacy may be another factor contributing to the limited number of existing pharmacokinetic models for orally inhaled drugs. Investigations comprising in vitro experiments, clinical studies, and more sophisticated mathematical approaches are considered to be necessary for elucidating these highly complex pulmonary processes. With this additional knowledge, the PBPK approach might gain additional attractiveness. Currently, (semi-)mechanistic modeling offers an alternative to generate and investigate hypotheses and to more mechanistically understand the pulmonary and systemic pharmacokinetics after oral drug inhalation including the impact of pulmonary diseases.KEY WORDS: inhalation, mathematical models, physiologically based pharmacokinetic models, population pharmacokinetics, pulmonary absorption     

Pharmacogenomics of the triazole antifungal agent voriconazole

Genetic polymorphisms in drug-metabolizing enzymes are frequently responsible for high variability in the pharmacokinetics of certain drugs leading to large variations in drug efficacy and adverse drug effects, or large ranges of the doses required for optimal drug efficacy. Voriconazole is a triazole antifungal agent which has been available for several years and has potent in vitro and in vivo activity against a broad spectrum of medically important pathogens, including Aspergillus, Cryptococcus and Candida. Voriconazole is extensively metabolized by the cytochrome P450 system with CYP2C19 being the major route for elimination. Thus, polymorphisms in the CYP2C19 gene have substantial impact on the pharmacokinetics of voriconazole and its interactions with other drugs. This article summarizes the current knowledge regarding CYP2C19 and discusses the influences of other drug-metabolizing enzymes and drug transporters on voriconazole disposition.     

In vivo pharmacokinetics, activation of MAPK signaling and induction of phase II/III drug metabolizing enzymes/transporters by cancer chemopreventive compound BHA in the mice

Phenolic antioxidant butylated hydroxyanisole (BHA) is a commonly used food preservative with broad biological activities, including protection against chemical-induced carcinogenesis, acute toxicity of chemicals, modulation of macromolecule synthesis and immune response, induction of phase II detoxifying enzymes, as well as its undesirable potential tumor-promoting activities. Understanding the molecular basis underlying these diverse biological actions of BHA is thus of great importance. Here we studied the pharmacokinetics, activation of signaling kinases and induction of phase II/III drug metabolizing enzymes/transporter gene expression by BHA in the mice. The peak plasma concentration of BHA achieved in our current study after oral administration of 200 mg/kg BHA was around 10 microM. This in vivo concentration might offer some insights for the many in vitro cell culture studies on signal transduction and induction of phase II genes using similar concentrations. The oral bioavailability (F) of BHA was about 43% in the mice. In the mouse liver, BHA induced the expression of phase II genes including NQO-1, HO-1, gamma-GCS, GST-pi and UGT 1A6, as well as some of the phase III transporter genes, such as MRP1 and Slcolb2. In addition, BHA activated distinct mitogen-activated protein kinases (MAPKs), c-Jun N-terminal kinase (JNK), extracellular signal-regulated protein kinase (ERK), as well as p38, suggesting that the MAPK pathways may play an important role in early signaling events leading to the regulation of gene expression including phase II drug metabolizing and some phase III drug transporter genes. This is the first study to demonstrate the in vivo pharmacokinetics of BHA, the in vivo activation of MAPK signaling proteins, as well as the in vivo induction of Phase II/III drug metabolizing enzymes/transporters in the mouse livers.     

Simulation of Monoclonal Antibody Pharmacokinetics in HumansUsing a Minimal Physiologically Based Model

Compared to small chemical molecules, monoclonal antibodies and Fc-containing derivatives (mAbs) have unique pharmacokinetic behaviour characterised by relatively poor cellular permeability, minimal renal filtration, binding to FcRn, target-mediated drug disposition, and disposition via lymph. A minimal physiologically based pharmacokinetic (PBPK) model to describe the pharmacokinetics of mAbs in humans was developed. Within the model, the body is divided into three physiological compartments; plasma, a single tissue compartment and lymph. The tissue compartment is further sub-divided into vascular, endothelial and interstitial spaces. The model simultaneously describes the levels of endogenous IgG and exogenous mAbs in each compartment and sub-compartment and, in particular, considers the competition of these two species for FcRn binding in the endothelial space. A Monte-Carlo sampling approach is used to simulate the concentrations of endogenous IgG and mAb in a human population. Existing targeted-mediated drug disposition (TMDD) models are coupled with the minimal PBPK model to provide a general platform for simulating the pharmacokinetics of therapeutic antibodies using primarily pre-clinical data inputs. The feasibility of utilising pre-clinical data to parameterise the model and to simulate the pharmacokinetics of adalimumab and an anti-ALK1 antibody (PF-03446962) in a population of individuals was investigated and results were compared to published clinical data.     

生理药代动力学模型在药物基因组学中的应用

生理药代动力学模型作为药代动力学的新兴工具,可以有效用于表征各种情况下由基因组学多态性导致的体内药动学变化.本文介绍生理药代动力学模型在药物基因组学研究中的基本原理并讨论其适用性,从I相代谢酶、II相代谢酶和转运体的基因多态性角度总结了模型的应用情况,以期了解药物基因组学的新型研究手段,为促进药物临床应用和提高药物研发...