1
|
De Carlo A, Tosca EM, Melillo N, Magni P. A two-stages global sensitivity analysis by using the δ sensitivity index in presence of correlated inputs: application on a tumor growth inhibition model based on the dynamic energy budget theory. J Pharmacokinet Pharmacodyn 2023; 50:395-409. [PMID: 37422844 PMCID: PMC10460734 DOI: 10.1007/s10928-023-09872-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2022] [Accepted: 06/16/2023] [Indexed: 07/11/2023]
Abstract
Global sensitivity analysis (GSA) evaluates the impact of variability and/or uncertainty of the model parameters on given model outputs. GSA is useful for assessing the quality of Pharmacometric model inference. Indeed, model parameters can be affected by high (estimation) uncertainty due to the sparsity of data. Independence between model parameters is a common assumption of GSA methods. However, ignoring (known) correlations between parameters may alter model predictions and, then, GSA results. To address this issue, a novel two-stages GSA technique based on the δ index, which is well-defined also in presence of correlated parameters, is here proposed. In the first step, statistical dependencies are neglected to identify parameters exerting causal effects. Correlations are introduced in the second step to consider the real distribution of the model output and investigate also the 'indirect' effects due to the correlation structure. The proposed two-stages GSA strategy was applied, as case study, to a preclinical tumor-in-host-growth inhibition model based on the Dynamic Energy Budget theory. The aim is to evaluate the impact of the model parameter estimate uncertainty (including correlations) on key model-derived metrics: the drug threshold concentration for tumor eradication, the tumor volume doubling time and a new index evaluating the drug efficacy-toxicity trade-off. This approach allowed to rank parameters according to their impact on the output, discerning whether a parameter mainly exerts a causal or 'indirect' effect. Thus, it was possible to identify uncertainties that should be necessarily reduced to obtain robust predictions for the outputs of interest.
Collapse
Affiliation(s)
- Alessandro De Carlo
- Electrical, Computer and Biomedical Engineering, University of Pavia, Pavia, Italy
| | - Elena Maria Tosca
- Electrical, Computer and Biomedical Engineering, University of Pavia, Pavia, Italy
| | - Nicola Melillo
- Electrical, Computer and Biomedical Engineering, University of Pavia, Pavia, Italy
- Systems Forecasting UK Ltd, Lancaster, UK
| | - Paolo Magni
- Electrical, Computer and Biomedical Engineering, University of Pavia, Pavia, Italy
| |
Collapse
|
2
|
Carlo AD, Tosca EM, Melillo N, Magni P. mvLognCorrEst: an R package for sampling from multivariate lognormal distributions and estimating correlations from uncomplete correlation matrix. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2023; 235:107517. [PMID: 37040682 DOI: 10.1016/j.cmpb.2023.107517] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/04/2022] [Revised: 03/15/2023] [Accepted: 03/27/2023] [Indexed: 05/08/2023]
Abstract
BACKGROUND AND OBJECTIVE Pharmacometrics (PMX) is a quantitative discipline which supports decision-making processes in all stages of drug development. PMX leverages Modeling and Simulations (M&S), which represents a powerful tool to characterize and predict the behavior and the effect of a drug. M&S-based methods, such as Sensitivity Analysis (SA) and Global Sensitivity Analysis (GSA), are gaining interest in PMX as they allow the evaluation of model-informed inference quality. Simulations should be correctly designed to obtain reliable results. Neglecting correlations between model parameters can significantly alter the results of simulations. However, the introduction of a correlation structure between model parameters can cause some issues. Sampling from a multivariate lognormal distribution, which is the typically distribution assumed for PMX model parameters, is not straightforward when a correlation structure is introduced. Indeed, correlations need to respect some constraints which depend by the CVs (i.e., coefficients of variation) of lognormal variables. In addition, when correlation matrices have some unspecified values, they should be properly fixed preserving the positive semi-definiteness of the correlation structure. In this paper, we present mvLognCorrEst, an R package developed to address these issues. METHODS The proposed sampling strategy was based on reconducting the extraction from the multivariate lognormal distribution of interest to the underlying Normal distribution. However, with high lognormal CVs, a positive semi-definite Normal covariance matrix cannot be obtained due to the violation of some theoretical constraints. In these cases, the Normal covariance matrix was approximated to its nearest positive definite matrix using Frobenius norm as matrix distance. For the estimation of unknown correlations terms, the graph theory was used to represent the correlation structure as weighed undirected graph. Plausible value ranges for the unspecified correlations were derived considering the paths between variables. Then, their estimation was performed by solving a constrained optimization problem. RESULTS Package functions are presented and applied on a real case study, that is the GSA of a PMX model that has been recently developed to support preclinical oncological studies. CONCLUSIONS mvLognCorrEst package is an R tool to support simulation-based analysis for which sampling from multivariate lognormal distributions with correlated variables and/or estimation of partially defined correlation matrix are required.
Collapse
Affiliation(s)
- Alessandro De Carlo
- Electrical, Computer and Biomedical Engineering, University of Pavia, Pavia, Italy.
| | - Elena Maria Tosca
- Electrical, Computer and Biomedical Engineering, University of Pavia, Pavia, Italy.
| | - Nicola Melillo
- Electrical, Computer and Biomedical Engineering, University of Pavia, Pavia, Italy; Systems Forecasting UK Ltd, Lancaster, UK.
| | - Paolo Magni
- Electrical, Computer and Biomedical Engineering, University of Pavia, Pavia, Italy.
| |
Collapse
|
3
|
Tsakalozou E, Alam K, Babiskin A, Zhao L. Physiologically-Based Pharmacokinetic Modeling to Support Determination of Bioequivalence for Dermatological Drug Products: Scientific and Regulatory Considerations. Clin Pharmacol Ther 2021; 111:1036-1049. [PMID: 34231211 DOI: 10.1002/cpt.2356] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2021] [Accepted: 06/11/2021] [Indexed: 12/30/2022]
Abstract
Physiologically-based pharmacokinetic (PBPK) modeling and simulation provides mechanism-based predictions of the pharmacokinetics of an active ingredient following its administration in humans. Dermal PBPK models describe the skin permeation and disposition of the active ingredient following the application of a dermatological product on the skin of virtual healthy and diseased human subjects. These models take into account information on product quality attributes, physicochemical properties of the active ingredient and skin (patho)physiology, and their interplay with each other. Regulatory and product development decision makers can leverage these quantitative tools to identify factors impacting local and systemic exposure. In the realm of generic drug products, the number of US Food and Drug Administratioin (FDA) interactions that use dermal PBPK modeling to support alternative bioequivalence (BE) approaches is increasing. In this report, we share scientific considerations on the development, verification and validation (V&V), and application of PBPK models within the context of a virtual BE assessment for dermatological drug products. We discuss the challenges associated with model V&V for these drug products stemming from the fact that target-site active ingredient concentrations are typically not measurable. Additionally, there are no established relationships between local and systemic PK profiles, when the latter are quantifiable. To that end, we detail a multilevel model V&V approach involving validation for the model of the drug product of interest coupled with the overall assessment of the modeling platform in use while leveraging in vitro and in vivo data related to local and systemic bioavailability.
Collapse
Affiliation(s)
- Eleftheria Tsakalozou
- Division of Quantitative Methods and Modeling (DQMM), Office of Research and Standards (ORS), Office of Generic Drugs (OGD), Center for Drug Evaluation and Research (CDER), US Food and Drug Administration (FDA), Silver Spring, Maryland, USA
| | - Khondoker Alam
- Division of Quantitative Methods and Modeling (DQMM), Office of Research and Standards (ORS), Office of Generic Drugs (OGD), Center for Drug Evaluation and Research (CDER), US Food and Drug Administration (FDA), Silver Spring, Maryland, USA
| | - Andrew Babiskin
- Division of Quantitative Methods and Modeling (DQMM), Office of Research and Standards (ORS), Office of Generic Drugs (OGD), Center for Drug Evaluation and Research (CDER), US Food and Drug Administration (FDA), Silver Spring, Maryland, USA
| | - Liang Zhao
- Division of Quantitative Methods and Modeling (DQMM), Office of Research and Standards (ORS), Office of Generic Drugs (OGD), Center for Drug Evaluation and Research (CDER), US Food and Drug Administration (FDA), Silver Spring, Maryland, USA
| |
Collapse
|
4
|
In Silico Modeling and Simulation to Guide Bioequivalence Testing for Oral Drugs in a Virtual Population. Clin Pharmacokinet 2021; 60:1373-1385. [PMID: 34191255 DOI: 10.1007/s40262-021-01045-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/28/2021] [Indexed: 12/18/2022]
Abstract
Model-informed drug discovery and development (MID3) shows great advantages in facilitating drug development. A physiologically based pharmacokinetic model is one of the powerful computational approaches of MID3, and the emerging field of virtual bioequivalence is well recognized to be the future of the physiologically based pharmacokinetic model. Based on the translational link between in vitro, in silico, and in vivo, virtual bioequivalence study can evaluate the similarity and potential difference of pharmacokinetic and clinical performance between test and reference formulations. With the aid of virtual bioequivalence study, the pivotal information of clinical trials can be provided to streamline the development for both new and generic drugs. However, a regulatory framework of virtual bioequivalence study has not reached its full maturity. Therefore, this article aims to present an overview of the current status of bioequivalence study, identify the framework of virtual bioequivalence studies for oral drugs, and also discuss the future opportunities of virtual bioequivalence in supporting the waiver and optimization of in vivo clinical trials.
Collapse
|
5
|
Melillo N, Darwich AS. A latent variable approach to account for correlated inputs in global sensitivity analysis. J Pharmacokinet Pharmacodyn 2021; 48:671-686. [PMID: 34032996 PMCID: PMC8405496 DOI: 10.1007/s10928-021-09764-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2021] [Accepted: 05/06/2021] [Indexed: 12/13/2022]
Abstract
In drug development decision-making is often supported through model-based methods, such as physiologically-based pharmacokinetics (PBPK). Global sensitivity analysis (GSA) is gaining use for quality assessment of model-informed inference. However, the inclusion and interpretation of correlated factors in GSA has proven an issue. Here we developed and evaluated a latent variable approach for dealing with correlated factors in GSA. An approach was developed that describes the correlation between two model inputs through the causal relationship of three independent factors: the latent variable and the unique variances of the two correlated parameters. The latent variable approach was applied to a set of algebraic models and a case from PBPK. Then, this method was compared to Sobol’s GSA assuming no correlations, Sobol’s GSA with groups and the Kucherenko approach. For the latent variable approach, GSA was performed with Sobol’s method. By using the latent variable approach, it is possible to devise a unique and easy interpretation of the sensitivity indices while maintaining the correlation between the factors. Compared methods either consider the parameters independent, group the dependent variables into one unique factor or present difficulties in the interpretation of the sensitivity indices. In situations where GSA is called upon to support model-informed decision-making, the latent variable approach offers a practical method, in terms of ease of implementation and interpretability, for applying GSA to models with correlated inputs that does not violate the independence assumption. Prerequisites and limitations of the approach are discussed.
Collapse
Affiliation(s)
- Nicola Melillo
- Centre for Applied Pharmacokinetic Research, Division of Pharmacy & Optometry, School of Health Sciences, The University of Manchester, Manchester, UK
| | - Adam S Darwich
- Division of Health Informatics and Logistics, Department of Biomedical Engineering and Health Systems, KTH Royal Institute of Technology, Stockholm, Sweden.
| |
Collapse
|
6
|
Wu F, Cristofoletti R, Zhao L, Rostami‐Hodjegan A. Scientific considerations to move towards biowaiver for biopharmaceutical classification system class III drugs: How modeling and simulation can help. Biopharm Drug Dispos 2021; 42:118-127. [DOI: 10.1002/bdd.2274] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2021] [Revised: 02/16/2021] [Accepted: 02/21/2021] [Indexed: 12/11/2022]
Affiliation(s)
- Fang Wu
- Division of Quantitative Methods and Modeling Office of Research and Standards Office of Generic Drugs Center for Drug Evaluation and Research U.S. Food and Drug Administration Silver Spring Maryland USA
| | - Rodrigo Cristofoletti
- Department of Pharmaceutics Center for Pharmacometrics and Systems Pharmacology College of Pharmacy University of Florida Orlando Florida USA
| | - Liang Zhao
- Division of Quantitative Methods and Modeling Office of Research and Standards Office of Generic Drugs Center for Drug Evaluation and Research U.S. Food and Drug Administration Silver Spring Maryland USA
| | - Amin Rostami‐Hodjegan
- Centre for Applied Pharmacokinetic Research University of Manchester Manchester UK
- Certara UK Limited Sheffield UK
| |
Collapse
|
7
|
Sattar H, Jadoon SS, Yang N, Li S, Xu M, Han Y, Ramzan A, Li W. Role of Glucuronidation Pathway in Quetiapine Metabolism: An In vivo Drug-Drug Interaction Study between Quetiapine and Probenecid. SAUDI JOURNAL OF MEDICINE AND MEDICAL SCIENCES 2020; 8:196-200. [PMID: 32952511 PMCID: PMC7485652 DOI: 10.4103/sjmms.sjmms_441_19] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2019] [Revised: 12/24/2019] [Accepted: 02/04/2020] [Indexed: 12/03/2022] Open
Abstract
Background: Uridine 5'-diphospho-glucuronosyltransferase (UGT) enzymes play a significant role in the metabolism of quetiapine, and coadministration with a UGT inhibitor/inducer drug may change its pharmacokinetic profile. Objective: The objective of this study was to assess the impact of probenecid, a UGT enzyme inhibitor, on the pharmacokinetic profile of quetiapine. Materials and Methods: Twelve treatment-naïve, 7-week-old male Sprague–Dawley rats (weighting 161 ± 22 g) were randomly and equally divided into control, quetiapine-alone and quetiapine plus probenecid groups. The quetiapine plus probenecid group received a single oral dose of probenecid (50 mg/kg) followed by 50 mg/kg of quetiapine; the quetiapine-alone group only received 50 mg/kg of quetiapine. Blood samples (0.2 ml) were collected from all rats after 0, 0.25, 0.5, 1, 2, 4, 6, 8, 10, 12 and 24 h of the drug administration in heparinized tubes. The pre-established liquid chromatography–mass spectrometry method was utilized to ascertain the plasma concentration of quetiapine and the control group was used to prepare the controlled standard. Results: Significant pharmacokinetic differences were observed between the quetiapine-alone and quetiapine plus probenecid groups in terms of Cmax (392 ± 209 vs. 1323 ± 343 ug/L, respectively, P = 0.004), AUC0-∞ (P = 0.04) and Tmax (P = 0.004). Further, in the combined drug group, there was a decrease in drug clearance (CL/F) (from 27 ± 11 to 16 ± 3 L/h/kg; P = 0.005) and an increase in the volume of distribution (Vd) (P = 0.01), but there was no significant difference between both groups in terms of half-lives (P = 0.27). No significant within-group variability of pharmacokinetic parameters was observed (P = 0.25). Conclusion: The results of this animal study suggest that glucuronidation by UGT enzyme system may also play an important role in quetiapine metabolism, which, if proven in future human studies, would imply that the bioavailability and pharmacokinetic parameters of quetiapine may require alterations when co-administered with probenecid to avoid development of quetiapine toxicity.
Collapse
Affiliation(s)
- Haseeb Sattar
- Department of Clinical Pharmacy, Wuhan Union Hospital Affiliated to Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
| | - Sarmad Sheraz Jadoon
- Department of Pharmacology, Hubei University of Traditional Chinese Medicine, Wuhan, China
| | - Ni Yang
- Department of Clinical Pharmacy, Wuhan Union Hospital Affiliated to Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
| | - Shihong Li
- Department of Clinical Pharmacy, Wuhan Union Hospital Affiliated to Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
| | - Mingzhen Xu
- Department of Clinical Pharmacy, Wuhan Union Hospital Affiliated to Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
| | - Yong Han
- Department of Clinical Pharmacy, Wuhan Union Hospital Affiliated to Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
| | - Adil Ramzan
- Department of Internal Medicine, Pakistan institute of Medical Sciences Affiliated to Shaheed Zulfiqar Ali Bhutto Medical University, Islamabad, Pakistan
| | - Weiyong Li
- Department of Clinical Pharmacy, Wuhan Union Hospital Affiliated to Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
| |
Collapse
|
8
|
Melillo N, Grandoni S, Cesari N, Brogin G, Puccini P, Magni P. Inter-compound and Intra-compound Global Sensitivity Analysis of a Physiological Model for Pulmonary Absorption of Inhaled Compounds. AAPS J 2020; 22:116. [PMID: 32862303 PMCID: PMC7456635 DOI: 10.1208/s12248-020-00499-0] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2020] [Accepted: 08/06/2020] [Indexed: 12/25/2022] Open
Abstract
In recent years, global sensitivity analysis (GSA) has gained interest in physiologically based pharmacokinetics (PBPK) modelling and simulation from pharmaceutical industry, regulatory authorities, and academia. With the case study of an in-house PBPK model for inhaled compounds in rats, the aim of this work is to show how GSA can contribute in PBPK model development and daily use. We identified two types of GSA that differ in the aims and, thus, in the parameter variability: inter-compound and intra-compound GSA. The inter-compound GSA aims to understand which are the parameters that mostly influence the variability of the metrics of interest in the whole space of the drugs' properties, and thus, it is useful during the model development. On the other hand, the intra-compound GSA aims to highlight how much the uncertainty associated with the parameters of a given drug impacts the uncertainty in the model prediction and so, it is useful during routine PBPK use. In this work, inter-compound GSA highlighted that dissolution- and formulation-related parameters were mostly important for the prediction of the fraction absorbed, while the permeability is the most important parameter for lung AUC and MRT. Intra-compound GSA highlighted that, for all the considered compounds, the permeability was one of the most important parameters for lung AUC, MRT and plasma MRT, while the extraction ratio and the dose for the plasma AUC. GSA is a crucial instrument for the quality assessment of model-based inference; for this reason, we suggest its use during both PBPK model development and use.
Collapse
Affiliation(s)
- Nicola Melillo
- Laboratory of Bioinformatics, Mathematical Modelling and Synthetic Biology, Department of Electrical, Computer and Biomedical Engineering, Università degli Studi di Pavia, Via Ferrata 5, I-27100, Pavia, Italy
| | - Silvia Grandoni
- Laboratory of Bioinformatics, Mathematical Modelling and Synthetic Biology, Department of Electrical, Computer and Biomedical Engineering, Università degli Studi di Pavia, Via Ferrata 5, I-27100, Pavia, Italy
| | - Nicola Cesari
- Pharmacokinetics, Biochemistry and Metabolism Department, Chiesi Farmaceutici S.p.A., Parma, Italy
| | - Giandomenico Brogin
- Pharmacokinetics, Biochemistry and Metabolism Department, Chiesi Farmaceutici S.p.A., Parma, Italy
| | - Paola Puccini
- Pharmacokinetics, Biochemistry and Metabolism Department, Chiesi Farmaceutici S.p.A., Parma, Italy
| | - Paolo Magni
- Laboratory of Bioinformatics, Mathematical Modelling and Synthetic Biology, Department of Electrical, Computer and Biomedical Engineering, Università degli Studi di Pavia, Via Ferrata 5, I-27100, Pavia, Italy.
| |
Collapse
|
9
|
Scotcher D, Arya V, Yang X, Zhao P, Zhang L, Huang S, Rostami‐Hodjegan A, Galetin A. A Novel Physiologically Based Model of Creatinine Renal Disposition to Integrate Current Knowledge of Systems Parameters and Clinical Observations. CPT-PHARMACOMETRICS & SYSTEMS PHARMACOLOGY 2020; 9:310-321. [PMID: 32441889 PMCID: PMC7306622 DOI: 10.1002/psp4.12509] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/24/2019] [Accepted: 02/16/2020] [Indexed: 01/11/2023]
Abstract
Creatinine is the most common clinical biomarker of renal function. As a substrate for renal transporters, its secretion is susceptible to inhibition by drugs, resulting in transient increase in serum creatinine and false impression of damage to kidney. Novel physiologically based models for creatinine were developed here and (dis)qualified in a stepwise manner until consistency with clinical data. Data from a matrix of studies were integrated, including systems data (common to all models), proteomics-informed in vitro-in vivo extrapolation of all relevant transporter clearances, exogenous administration of creatinine (to estimate endogenous synthesis rate), and inhibition of different renal transporters (11 perpetrator drugs considered for qualification during creatinine model development and verification on independent data sets). The proteomics-informed bottom-up approach resulted in the underprediction of creatinine renal secretion. Subsequently, creatinine-trimethoprim clinical data were used to inform key model parameters in a reverse translation manner, highlighting best practices and challenges for middle-out optimization of mechanistic models.
Collapse
Affiliation(s)
- Daniel Scotcher
- Centre for Applied Pharmacokinetic ResearchUniversity of ManchesterManchesterUK
| | - Vikram Arya
- Office of Clinical PharmacologyOffice of Translational SciencesCentre for Drug Evaluation and ResearchUS Food and Drug AdministrationSilver SpringMarylandUSA
| | - Xinning Yang
- Office of Clinical PharmacologyOffice of Translational SciencesCentre for Drug Evaluation and ResearchUS Food and Drug AdministrationSilver SpringMarylandUSA
| | - Ping Zhao
- Office of Clinical PharmacologyOffice of Translational SciencesCentre for Drug Evaluation and ResearchUS Food and Drug AdministrationSilver SpringMarylandUSA
- Present address:
Bill & Melinda Gates FoundationSeattleWashingtonUSA
| | - Lei Zhang
- Office of Research and StandardsOffice of Generic DrugsCentre for Drug Evaluation and ResearchUS Food and Drug AdministrationSilver SpringMarylandUSA
| | - Shiew‐Mei Huang
- Office of Clinical PharmacologyOffice of Translational SciencesCentre for Drug Evaluation and ResearchUS Food and Drug AdministrationSilver SpringMarylandUSA
| | - Amin Rostami‐Hodjegan
- Centre for Applied Pharmacokinetic ResearchUniversity of ManchesterManchesterUK
- CertaraSheffieldUK
| | - Aleksandra Galetin
- Centre for Applied Pharmacokinetic ResearchUniversity of ManchesterManchesterUK
| |
Collapse
|
10
|
A novel mass spectrometry method for the absolute quantification of several cytochrome P450 and uridine 5'-diphospho-glucuronosyltransferase enzymes in the human liver. Anal Bioanal Chem 2020; 412:1729-1740. [PMID: 32030490 DOI: 10.1007/s00216-020-02445-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2019] [Revised: 12/22/2019] [Accepted: 01/21/2020] [Indexed: 10/25/2022]
Abstract
Cytochrome P450 (CYP450) and 5'-diphosphate glucuronosyltransferases (UGT) are the two major families of drug-metabolizing enzymes in the human liver microsome (HLM). As a result of their frequent abundance fluctuation among populations, the accurate quantification of these enzymes in different individuals is important for designing patient-specific dosage regimens in the framework of precision medicine. The preparation and quantification of internal standards is an essential step for the quantitative analysis of enzymes. However, the commonly employed stable isotope labeling-based strategy (QconCAT) suffers from requiring very expensive isotopic reagents, tedious experimental procedures, and long labeling times. Furthermore, arginine-to-proline conversion during metabolic isotopic labeling compromises the quantification accuracy. Therefore, we present a new strategy that replaces stable isotope-labeled amino acids with lanthanide labeling for the preparation and quantification of QconCAT internal standard peptides, which leads to a threefold reduction in the reagent costs and a fivefold reduction in the time consumed. The absolute amount of trypsin-digested QconCAT peptides can be obtained by lanthanide labeling and inductively coupled plasma-optical emission spectrometry (ICP-OES) analysis with a high quantification accuracy (%RE < 20%). By taking advantage of the highly selective and facile ICP-OES procedure and multiplexed large-scale absolute target protein quantification using biological mass spectrometry, this strategy was successfully used for the absolute quantification of drug-metabolizing enzymes. We obtained good linearity (correlation coefficient > 0.95) over concentrations spanning 2.5 orders of magnitude with improved sensitivity (limit of quantification = 2 fmol) in nine HLM samples, indicating the potential of this method for large-scale absolute target protein quantification in clinical samples. Graphical abstract.
Collapse
|
11
|
Yau E, Olivares-Morales A, Gertz M, Parrott N, Darwich AS, Aarons L, Ogungbenro K. Global Sensitivity Analysis of the Rodgers and Rowland Model for Prediction of Tissue: Plasma Partitioning Coefficients: Assessment of the Key Physiological and Physicochemical Factors That Determine Small-Molecule Tissue Distribution. AAPS JOURNAL 2020; 22:41. [PMID: 32016678 DOI: 10.1208/s12248-020-0418-7] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/14/2019] [Accepted: 01/07/2020] [Indexed: 12/14/2022]
Abstract
In physiologically based pharmacokinetic (PBPK) modelling, the large number of input parameters, limited amount of available data and the structural model complexity generally hinder simultaneous estimation of uncertain and/or unknown parameters. These parameters are generally subject to estimation. However, the approaches taken for parameter estimation vary widely. Global sensitivity analyses are proposed as a method to systematically determine the most influential parameters that can be subject to estimation. Herein, a global sensitivity analysis was conducted to identify the key drug and physiological parameters influencing drug disposition in PBPK models and to potentially reduce the PBPK model dimensionality. The impact of these parameters was evaluated on the tissue-to-unbound plasma partition coefficients (Kpus) predicted by the Rodgers and Rowland model using Latin hypercube sampling combined to partial rank correlation coefficients (PRCC). For most drug classes, PRCC showed that LogP and fraction unbound in plasma (fup) were generally the most influential parameters for Kpu predictions. For strong bases, blood:plasma partitioning was one of the most influential parameter. Uncertainty in tissue composition parameters had a large impact on Kpu and Vss predictions for all classes. Among tissue composition parameters, changes in Kpu outputs were especially attributed to changes in tissue acidic phospholipid concentrations and extracellular protein tissue:plasma ratio values. In conclusion, this work demonstrates that for parameter estimation involving PBPK models and dimensionality reduction purposes, less influential parameters might be assigned fixed values depending on the parameter space, while influential parameters could be subject to parameters estimation.
Collapse
Affiliation(s)
- Estelle Yau
- Centre for Applied Pharmacokinetic Research (CAPKR), The University of Manchester, Manchester, UK.,Roche Pharma and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd, Grenzacherstrasse 124, 4070, Basel, Switzerland
| | - Andrés Olivares-Morales
- Roche Pharma and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd, Grenzacherstrasse 124, 4070, Basel, Switzerland.
| | - Michael Gertz
- Roche Pharma and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd, Grenzacherstrasse 124, 4070, Basel, Switzerland
| | - Neil Parrott
- Roche Pharma and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd, Grenzacherstrasse 124, 4070, Basel, Switzerland
| | - Adam S Darwich
- Centre for Applied Pharmacokinetic Research (CAPKR), The University of Manchester, Manchester, UK.,Logistics and Informatics in Health Care, School of Engineering Sciences in Chemistry, Biotechnology and Health (CBH), KTH Royal Institute of Technology, Stockholm, Sweden
| | - Leon Aarons
- Roche Pharma and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd, Grenzacherstrasse 124, 4070, Basel, Switzerland
| | - Kayode Ogungbenro
- Roche Pharma and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd, Grenzacherstrasse 124, 4070, Basel, Switzerland
| |
Collapse
|
12
|
Couto N, Al-Majdoub ZM, Gibson S, Davies PJ, Achour B, Harwood MD, Carlson G, Barber J, Rostami-Hodjegan A, Warhurst G. Quantitative Proteomics of Clinically Relevant Drug-Metabolizing Enzymes and Drug Transporters and Their Intercorrelations in the Human Small Intestine. Drug Metab Dispos 2020; 48:245-254. [PMID: 31959703 PMCID: PMC7076527 DOI: 10.1124/dmd.119.089656] [Citation(s) in RCA: 64] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2019] [Accepted: 12/23/2019] [Indexed: 01/02/2023] Open
Abstract
The levels of drug-metabolizing enzymes (DMEs) and transporter proteins in the human intestine are pertinent to determine oral drug bioavailability. Despite the paucity of reports on such measurements, it is well recognized that these values are essential for translating in vitro data on drug metabolism and transport to predict drug disposition in gut wall. In the current study, clinically relevant DMEs [cytochrome P450 (P450) and uridine 5′-diphospho-glucuronosyltransferase (UGT)] and drug transporters were quantified in total mucosal protein preparations from the human jejunum (n = 4) and ileum (n = 12) using quantification concatemer–based targeted proteomics. In contrast to previous reports, UGT2B15 and organic anion-transporting polypeptide 1 (OATP1A2) were quantifiable in all our samples. Overall, no significant disparities in protein expression were observed between jejunum and ileum. Relative mRNA expression for drug transporters did not correlate with the abundance of their cognate protein, except for P-glycoprotein 1 (P-gp) and organic solute transporter subunit alpha (OST-α), highlighting the limitations of RNA as a surrogate for protein expression in dynamic tissues with high turnover. Intercorrelations were found within P450 [2C9-2C19 (P = 0.002, R2 = 0.63), 2C9–2J2 (P = 0.004, R2 = 0.40), 2D6-2J2 (P = 0.002, R2 = 0.50)] and UGT [1A1-2B7 (P = 0.02, R2 = 0.87)] family of enzymes. There were also correlations between P-gp and several other proteins [OST-α (P < 0.0001, R2 = 0.77), UGT1A6 (P = 0.009, R2 = 0.38), and CYP3A4 (P = 0.007, R2 = 0.30)]. Incorporating such correlations into building virtual populations is crucial for obtaining plausible characteristics of simulated individuals.
Collapse
Affiliation(s)
- Narciso Couto
- Centre for Applied Pharmacokinetic Research, University of Manchester, Manchester, United Kingdom (N.C., Z.M.A.-M., B.A., J.B., A.R.-H.); Gut Barrier Group, Inflammation and Repair, University of Manchester, Salford Royal NHS Trust, Salford, United Kingdom (S.G., P.J.D., G.C., G.W.); and Certara UK Limited (Simcyp Division), Sheffield, United Kingdom (M.D.H., A.R.-H.)
| | - Zubida M Al-Majdoub
- Centre for Applied Pharmacokinetic Research, University of Manchester, Manchester, United Kingdom (N.C., Z.M.A.-M., B.A., J.B., A.R.-H.); Gut Barrier Group, Inflammation and Repair, University of Manchester, Salford Royal NHS Trust, Salford, United Kingdom (S.G., P.J.D., G.C., G.W.); and Certara UK Limited (Simcyp Division), Sheffield, United Kingdom (M.D.H., A.R.-H.)
| | - Stephanie Gibson
- Centre for Applied Pharmacokinetic Research, University of Manchester, Manchester, United Kingdom (N.C., Z.M.A.-M., B.A., J.B., A.R.-H.); Gut Barrier Group, Inflammation and Repair, University of Manchester, Salford Royal NHS Trust, Salford, United Kingdom (S.G., P.J.D., G.C., G.W.); and Certara UK Limited (Simcyp Division), Sheffield, United Kingdom (M.D.H., A.R.-H.)
| | - Pamela J Davies
- Centre for Applied Pharmacokinetic Research, University of Manchester, Manchester, United Kingdom (N.C., Z.M.A.-M., B.A., J.B., A.R.-H.); Gut Barrier Group, Inflammation and Repair, University of Manchester, Salford Royal NHS Trust, Salford, United Kingdom (S.G., P.J.D., G.C., G.W.); and Certara UK Limited (Simcyp Division), Sheffield, United Kingdom (M.D.H., A.R.-H.)
| | - Brahim Achour
- Centre for Applied Pharmacokinetic Research, University of Manchester, Manchester, United Kingdom (N.C., Z.M.A.-M., B.A., J.B., A.R.-H.); Gut Barrier Group, Inflammation and Repair, University of Manchester, Salford Royal NHS Trust, Salford, United Kingdom (S.G., P.J.D., G.C., G.W.); and Certara UK Limited (Simcyp Division), Sheffield, United Kingdom (M.D.H., A.R.-H.)
| | - Matthew D Harwood
- Centre for Applied Pharmacokinetic Research, University of Manchester, Manchester, United Kingdom (N.C., Z.M.A.-M., B.A., J.B., A.R.-H.); Gut Barrier Group, Inflammation and Repair, University of Manchester, Salford Royal NHS Trust, Salford, United Kingdom (S.G., P.J.D., G.C., G.W.); and Certara UK Limited (Simcyp Division), Sheffield, United Kingdom (M.D.H., A.R.-H.)
| | - Gordon Carlson
- Centre for Applied Pharmacokinetic Research, University of Manchester, Manchester, United Kingdom (N.C., Z.M.A.-M., B.A., J.B., A.R.-H.); Gut Barrier Group, Inflammation and Repair, University of Manchester, Salford Royal NHS Trust, Salford, United Kingdom (S.G., P.J.D., G.C., G.W.); and Certara UK Limited (Simcyp Division), Sheffield, United Kingdom (M.D.H., A.R.-H.)
| | - Jill Barber
- Centre for Applied Pharmacokinetic Research, University of Manchester, Manchester, United Kingdom (N.C., Z.M.A.-M., B.A., J.B., A.R.-H.); Gut Barrier Group, Inflammation and Repair, University of Manchester, Salford Royal NHS Trust, Salford, United Kingdom (S.G., P.J.D., G.C., G.W.); and Certara UK Limited (Simcyp Division), Sheffield, United Kingdom (M.D.H., A.R.-H.)
| | - Amin Rostami-Hodjegan
- Centre for Applied Pharmacokinetic Research, University of Manchester, Manchester, United Kingdom (N.C., Z.M.A.-M., B.A., J.B., A.R.-H.); Gut Barrier Group, Inflammation and Repair, University of Manchester, Salford Royal NHS Trust, Salford, United Kingdom (S.G., P.J.D., G.C., G.W.); and Certara UK Limited (Simcyp Division), Sheffield, United Kingdom (M.D.H., A.R.-H.)
| | - Geoffrey Warhurst
- Centre for Applied Pharmacokinetic Research, University of Manchester, Manchester, United Kingdom (N.C., Z.M.A.-M., B.A., J.B., A.R.-H.); Gut Barrier Group, Inflammation and Repair, University of Manchester, Salford Royal NHS Trust, Salford, United Kingdom (S.G., P.J.D., G.C., G.W.); and Certara UK Limited (Simcyp Division), Sheffield, United Kingdom (M.D.H., A.R.-H.)
| |
Collapse
|
13
|
Quantitative mass spectrometry-based proteomics in the era of model-informed drug development: Applications in translational pharmacology and recommendations for best practice. Pharmacol Ther 2019; 203:107397. [DOI: 10.1016/j.pharmthera.2019.107397] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2019] [Accepted: 07/29/2019] [Indexed: 02/08/2023]
|