1
|
Jeong YS, Jusko WJ. A Complete Extension of Classical Hepatic Clearance Models Using Fractional Distribution Parameter f d in Physiologically Based Pharmacokinetics. J Pharm Sci 2024; 113:95-117. [PMID: 37279835 PMCID: PMC10902797 DOI: 10.1016/j.xphs.2023.05.019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2023] [Revised: 05/30/2023] [Accepted: 05/30/2023] [Indexed: 06/08/2023]
Abstract
The classical organ clearance models have been proposed to relate the plasma clearance CLp to probable mechanism(s) of hepatic clearance. However, the classical models assume the intrinsic capability of drug elimination (CLu,int) that is physically segregated from the vascular blood but directly acts upon the unbound drug concentration in the blood (fubCavg), and do not handle the transit-time delay between the inlet/outlet concentrations in their closed-form clearance equations. Therefore, we propose unified model structures that can address the internal blood concentration patterns of clearance organs in a more mechanistic/physiological manner, based on the fractional distribution parameter fd operative in PBPK. The basic partial/ordinary differential equations for four classical models are revisited/modified to yield a more complete set of extended clearance models, i.e., the Rattle, Sieve, Tube, and Jar models, which are the counterparts of the dispersion, series-compartment, parallel-tube, and well-stirred models. We demonstrate the feasibility of applying the resulting extended models to isolated perfused rat liver data for 11 compounds and an example dataset for in vitro-in vivo extrapolation of the intrinsic to the systemic clearances. Based on their feasibilities to handle such real data, these models may serve as an improved basis for applying clearance models in the future.
Collapse
Affiliation(s)
- Yoo-Seong Jeong
- Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, State University of New York at Buffalo, Buffalo, NY, 14214, USA
| | - William J Jusko
- Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, State University of New York at Buffalo, Buffalo, NY, 14214, USA.
| |
Collapse
|
2
|
De Sutter PJ, Gasthuys E, Vermeulen A. Comparison of monoclonal antibody disposition predictions using different physiologically based pharmacokinetic modelling platforms. J Pharmacokinet Pharmacodyn 2023:10.1007/s10928-023-09894-4. [PMID: 37952005 DOI: 10.1007/s10928-023-09894-4] [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: 06/08/2023] [Accepted: 10/17/2023] [Indexed: 11/14/2023]
Abstract
Physiologically based pharmacokinetic (PBPK) models can be used to leverage physiological and in vitro data to predict monoclonal antibody (mAb) concentrations in serum and tissues. However, it is currently not known how consistent predictions of mAb disposition are across PBPK modelling platforms. In this work PBPK simulations of IgG, adalimumab and infliximab were compared between three platforms (Simcyp, PK-Sim, and GastroPlus). Accuracy of predicted serum and tissue concentrations was assessed using observed data collected from the literature. Physiological and mAb related input parameters were also compared and sensitivity analyses were carried out to evaluate model behavior when input values were altered. Differences in serum kinetics of IgG between platforms were minimal for a dose of 1 mg/kg, but became more noticeable at higher dosages (> 100 mg/kg) and when reference (healthy) physiological input values were altered. Predicted serum concentrations of both adalimumab and infliximab were comparable across platforms, but were noticeably higher than observed values. Tissue concentrations differed remarkably between the platforms, both for total- and interstitial fluid (ISF) concentrations. The accuracy of total tissue concentrations was within a three-fold of observed values for all tissues, except for brain tissue concentrations, which were overpredicted. Predictions of tissue ISF concentrations were less accurate and were best captured by GastroPlus. Overall, these simulations show that the different PBPK platforms generally predict similar mAb serum concentrations, but variable tissue concentrations. Caution is therefore warranted when PBPK models are used to simulate effect site tissue concentrations of mAbs without data to verify the predictions.
Collapse
Affiliation(s)
- Pieter-Jan De Sutter
- Laboratory of Medical Biochemistry and Clinical Analysis, Department of Bioanalysis, Ghent University, Ottergemsesteenweg 460, 9000, Ghent, Belgium.
| | - Elke Gasthuys
- Laboratory of Medical Biochemistry and Clinical Analysis, Department of Bioanalysis, Ghent University, Ottergemsesteenweg 460, 9000, Ghent, Belgium
| | - An Vermeulen
- Laboratory of Medical Biochemistry and Clinical Analysis, Department of Bioanalysis, Ghent University, Ottergemsesteenweg 460, 9000, Ghent, Belgium
| |
Collapse
|
3
|
Cai H, Tao X, Shim J, Bauer RN, Bremer M, Bu W, LaMar J, Basile R, Dere E, Nguyen T, Laing S, Chan P, Yi T, Koerber JT, Sperinde G, Stefanich E. Mini-PBPK-Based Population Model and Covariate Analysis to Assess the Complex Pharmacokinetics and Pharmacodynamics of RO7449135, an Anti-KLK5/KLK7 Bispecific Antibody in Cynomolgus Monkeys. AAPS J 2023; 25:64. [PMID: 37353723 DOI: 10.1208/s12248-023-00829-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2023] [Accepted: 06/05/2023] [Indexed: 06/25/2023] Open
Abstract
RO7449135, an anti-kallikrein (KLK)5/KLK7 bispecific antibody, is in development as a potential therapy against Netherton's syndrome (NS). In cynomolgus monkey studies, RO7449135 bound to KLK5 and KLK7, causing considerable accumulation of total KLKs, but with non-dose-proportional increase. To understand the complex PKPD, a population model with covariate analysis was developed accounting for target binding in skin and migration of bound targets from skin to blood. The covariate analysis suggested the animal batch as the categorical covariate impacting the different KLK5 synthesis rates between the repeat-dose study and single-dose study, and the dose as continuous covariate impacting the internalization rate of the binary and ternary complexes containing KLK7. To comprehend the mechanism underlying, we hypothesized that inhibition of KLK5 by RO7449135 prevented its cleavage of the pro-enzyme of KLK7 (pro-KLK7) and altered the proportion between pro-KLK7 and KLK7. Besides the pro-KLK7, RO7449135 can interact with other proteins like LEKTI through KLK7 connection in a dose-dependent manner. The different high-order complexes formed by RO7449135 interacting with pro-KLK7 or LEKTI-like proteins can be subject to faster internalization rate. Accounting for the dose and animal batch as covariates, the model-predicted free target suppression is well aligned with the visual target engagement check. The population PKPD model with covariate analysis provides the scientific input for the complex PKPD analysis, successfully predicts the target suppression in cynomolgus monkeys, and thereby can be used for the human dose projection of RO7449135.
Collapse
Affiliation(s)
- Hao Cai
- Preclinical and Translational Pharmacokinetics and Pharmacodynamics, Genentech Inc., 1 DNA Way, South San Francisco, California, 94080, USA
| | - Xun Tao
- Preclinical and Translational Pharmacokinetics and Pharmacodynamics, Genentech Inc., 1 DNA Way, South San Francisco, California, 94080, USA
| | - Jeongsup Shim
- BioAnalytical Sciences, Genentech Inc., 1 DNA Way, South San Francisco, California, 94080, USA
| | - Rebecca N Bauer
- OMNI Biomarker Development, Genentech Inc., 1 DNA Way, South San Francisco, California, 94080, USA
| | - Meire Bremer
- OMNI Biomarker Development, Genentech Inc., 1 DNA Way, South San Francisco, California, 94080, USA
| | - Wei Bu
- BioAnalytical Sciences, Genentech Inc., 1 DNA Way, South San Francisco, California, 94080, USA
| | - Jason LaMar
- BioAnalytical Sciences, Genentech Inc., 1 DNA Way, South San Francisco, California, 94080, USA
| | - Rachel Basile
- BioAnalytical Sciences, Genentech Inc., 1 DNA Way, South San Francisco, California, 94080, USA
| | - Edward Dere
- Safety Assessment, Genentech Inc., 1 DNA Way, South San Francisco, California, 94080, USA
| | - Tien Nguyen
- Safety Assessment, Genentech Inc., 1 DNA Way, South San Francisco, California, 94080, USA
| | - Steven Laing
- Safety Assessment, Genentech Inc., 1 DNA Way, South San Francisco, California, 94080, USA
| | - Pamela Chan
- Biochemical and Cellular Pharmacology, Genentech Inc., 1 DNA Way, South San Francisco, California, 94080, USA
| | - Tangsheng Yi
- Discovery Immunology, Genentech Inc., 1 DNA Way, South San Francisco, California, 94080, USA
| | - James T Koerber
- Antibody Engineering, Genentech Inc., 1 DNA Way, South San Francisco, California, 94080, USA
| | - Gizette Sperinde
- BioAnalytical Sciences, Genentech Inc., 1 DNA Way, South San Francisco, California, 94080, USA
| | - Eric Stefanich
- Preclinical and Translational Pharmacokinetics and Pharmacodynamics, Genentech Inc., 1 DNA Way, South San Francisco, California, 94080, USA.
| |
Collapse
|
4
|
Lu D, Yadav R, Holder P, Chiang E, Sanjabi S, Poon V, Bernett M, Varma R, Liu K, Leung I, Bogaert L, Desjarlais J, Shivva V, Hosseini I, Ramanujan S. Complex PK-PD of an engineered IL-15/IL-15Rα-Fc fusion protein in cynomolgus monkeys: QSP modeling of lymphocyte dynamics. Eur J Pharm Sci 2023; 186:106450. [PMID: 37084985 DOI: 10.1016/j.ejps.2023.106450] [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/27/2022] [Revised: 03/29/2023] [Accepted: 04/18/2023] [Indexed: 04/23/2023]
Abstract
XmAb24306 is a lymphoproliferative interleukin (IL)-15/IL-15 receptor α (IL-15Rα) Fc-fusion protein currently under clinical investigation as an immunotherapeutic agent for cancer treatment. XmAb24306 contains mutations in IL-15 that attenuate its affinity to the heterodimeric IL-15 receptor βγ (IL-15R). We observe substantially prolonged pharmacokinetics (PK) (half-life ∼ 2.5 to 4.5 days) in single- and repeat-dose cynomolgus monkey (cyno) studies compared to wild-type IL-15 (half-life ∼ 1 hour), leading to increased exposure and enhanced and durable expansion of NK cells, CD8+ T cells and CD4-CD8- (double negative [DN]) T cells. Drug clearance varied with dose level and time post-dose, and PK exposure decreased upon repeated dosing, which we attribute to increased target-mediated drug disposition (TMDD) resulting from drug-induced lymphocyte expansion (i.e., pharmacodynamic (PD)-enhanced TMDD). We developed a quantitative systems pharmacology (QSP) model to quantify the complex PKPD behaviors due to the interactions of XmAb24306 with multiple cell types (CD8+, CD4+, DN T cells, and NK cells) in the peripheral blood (PB) and lymphoid tissues. The model, which includes nonspecific drug clearance, binding to and TMDD by IL15R differentially expressed on lymphocyte subsets, and resultant lymphocyte margination/migration out of PB, expansion in lymphoid tissues, and redistribution to the blood, successfully describes the systemic PK and lymphocyte kinetics observed in the cyno studies. Results suggest that after 3 doses of every-two-week (Q2W) doses up to 70 days, the relative contributions of each elimination pathway to XmAb24306 clearance are: DN T cells > NK cells > CD8+ T cells > nonspecific clearance > CD4+ T cells. Modeling suggests that observed cellular expansion in blood results from the influx of cells expanded by the drug in lymphoid tissues. The model is used to predict lymphoid tissue expansion and to simulate PK-PD for different dose regimens. Thus, the model provides insight into the mechanisms underlying the observed PK-PD behavior of an engineered cytokine and can serve as a framework for the rapid integration and analysis of data that emerges from ongoing clinical studies in cancer patients as single-agent or given in combination.
Collapse
Affiliation(s)
- Dan Lu
- Genentech, Inc., South San Francisco, CA, USA.
| | | | | | | | | | - Victor Poon
- Genentech, Inc., South San Francisco, CA, USA
| | | | | | - Ke Liu
- Xencor, Inc. Monrovia, CA, USA
| | | | | | | | | | | | | |
Collapse
|
5
|
Dawson DE, Lau C, Pradeep P, Sayre RR, Judson RS, Tornero-Velez R, Wambaugh JF. A Machine Learning Model to Estimate Toxicokinetic Half-Lives of Per- and Polyfluoro-Alkyl Substances (PFAS) in Multiple Species. TOXICS 2023; 11:98. [PMID: 36850973 PMCID: PMC9962572 DOI: 10.3390/toxics11020098] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/12/2022] [Revised: 01/09/2023] [Accepted: 01/18/2023] [Indexed: 06/18/2023]
Abstract
Per- and polyfluoroalkyl substances (PFAS) are a diverse group of man-made chemicals that are commonly found in body tissues. The toxicokinetics of most PFAS are currently uncharacterized, but long half-lives (t½) have been observed in some cases. Knowledge of chemical-specific t½ is necessary for exposure reconstruction and extrapolation from toxicological studies. We used an ensemble machine learning method, random forest, to model the existing in vivo measured t½ across four species (human, monkey, rat, mouse) and eleven PFAS. Mechanistically motivated descriptors were examined, including two types of surrogates for renal transporters: (1) physiological descriptors, including kidney geometry, for renal transporter expression and (2) structural similarity of defluorinated PFAS to endogenous chemicals for transporter affinity. We developed a classification model for t½ (Bin 1: <12 h; Bin 2: <1 week; Bin 3: <2 months; Bin 4: >2 months). The model had an accuracy of 86.1% in contrast to 32.2% for a y-randomized null model. A total of 3890 compounds were within domain of the model, and t½ was predicted using the bin medians: 4.9 h, 2.2 days, 33 days, and 3.3 years. For human t½, 56% of PFAS were classified in Bin 4, 7% were classified in Bin 3, and 37% were classified in Bin 2. This model synthesizes the limited available data to allow tentative extrapolation and prioritization.
Collapse
Affiliation(s)
- Daniel E. Dawson
- U.S. Environmental Protection Agency, Office of Research and Development, Center for Computational Toxicology and Exposure, 109 T.W. Alexander Drive, Research Triangle Park, NC 27711, USA
| | - Christopher Lau
- U.S. Environmental Protection Agency, Office of Research and Development, Center for Public Health and Environmental Assessment, 109 T.W. Alexander Drive, Research Triangle Park, NC 277011, USA
| | - Prachi Pradeep
- U.S. Environmental Protection Agency, Office of Research and Development, Center for Computational Toxicology and Exposure, 109 T.W. Alexander Drive, Research Triangle Park, NC 27711, USA
- Oak Ridge Institutes for Science and Education, Oak Ridge, TN 37830, USA
| | - Risa R. Sayre
- U.S. Environmental Protection Agency, Office of Research and Development, Center for Computational Toxicology and Exposure, 109 T.W. Alexander Drive, Research Triangle Park, NC 27711, USA
| | - Richard S. Judson
- U.S. Environmental Protection Agency, Office of Research and Development, Center for Computational Toxicology and Exposure, 109 T.W. Alexander Drive, Research Triangle Park, NC 27711, USA
| | - Rogelio Tornero-Velez
- U.S. Environmental Protection Agency, Office of Research and Development, Center for Computational Toxicology and Exposure, 109 T.W. Alexander Drive, Research Triangle Park, NC 27711, USA
| | - John F. Wambaugh
- U.S. Environmental Protection Agency, Office of Research and Development, Center for Computational Toxicology and Exposure, 109 T.W. Alexander Drive, Research Triangle Park, NC 27711, USA
| |
Collapse
|
6
|
Hu S, Datta-Mannan A, D'Argenio DZ. Monoclonal Antibody Pharmacokinetics in Cynomolgus Monkeys Following Subcutaneous Administration: Physiologically Based Model Predictions from Physiochemical Properties. AAPS J 2022; 25:5. [PMID: 36456779 DOI: 10.1208/s12248-022-00772-4] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2022] [Accepted: 10/31/2022] [Indexed: 12/03/2022] Open
Abstract
An integrated physiologically based modeling framework is presented for predicting pharmacokinetics and bioavailability of subcutaneously administered monoclonal antibodies in cynomolgus monkeys, based on in silico structure-derived metrics characterizing antibody size, overall charge, local charge, and hydrophobicity. The model accounts for antibody-specific differences in pinocytosis, transcapillary transport, local lymphatic uptake, and pre-systemic degradation at the subcutaneous injection site and reliably predicts the pharmacokinetics of five different wild-type mAbs and their Fc variants following intravenous and subcutaneous administration. Significant associations were found between subcutaneous injection site degradation rate and the antibody's local positive charge of its complementarity-determining region (R = 0.56, p = 0.0012), antibody pinocytosis rate and its overall positive charge (R = 0.59, p = 0.00063), and antibody paracellular transport and its overall charge together with hydrophobicity (R = 0.63, p = 0.00096). Based on these results, population simulations were performed to predict the relationship between bioavailability and antibody local positive charge. In addition, model simulations were conducted to calculate the relative contribution of absorption pathways (lymphatic and blood), pre-systemic degradation pathways (interstitial and lysosomal), and the influence of injection site lymph flow on antibody bioavailability and pharmacokinetics. The proposed physiologically based modeling framework integrates fundamental mechanisms governing antibody subcutaneous absorption and disposition, with structured-based physiochemical properties, to predict antibody bioavailability and pharmacokinetics in vivo.
Collapse
Affiliation(s)
- Shihao Hu
- Department of Biomedical Engineering, University of Southern California, Los Angeles, California, 90089, USA
| | - Amita Datta-Mannan
- Department of Exploratory Medicine and Pharmacology, Lilly Research Laboratories, Lilly Corporate Center, Indianapolis, Indiana, USA
| | - David Z D'Argenio
- Department of Biomedical Engineering, University of Southern California, Los Angeles, California, 90089, USA.
| |
Collapse
|
7
|
Salimi-Moosavi H, Soto M. A Non-radiometric Approach to Determine Tissue Vascular Blood Volume in Biodistribution Studies. AAPS J 2022; 24:116. [PMID: 36376552 DOI: 10.1208/s12248-022-00770-6] [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: 08/08/2022] [Accepted: 10/28/2022] [Indexed: 11/16/2022] Open
Abstract
The aim of this research was to develop a reliable non-radiometric method to measure the residual blood in tissue without the need for perfusion or radiometric measurements in biodistribution studies. It was found that the perfusion method not only was ineffective in removing blood from tissue, but also introduced additional variability in the determination of tissue drug exposure and was not reproducible across studies. In addition, the use of hemoglobin as an endogenous protein and biomarker for tissue blood content was studied and it was found that hemoglobin measurement in tissue was not a reliable and effective approach for determination of residual blood level in tissue. To evaluate an alternative method for addressing the tissue blood level in biodistribution studies, animals were dosed with a Residual Blood Determinant Reagent (RBDR) 5 min prior to tissue harvesting. The level of RBDR, an exogenous protein, was measured in whole blood homogenate and in tissue lysate. Based on the level of the RBDR, the vascular blood volume (VBV) in tissue was calculated and then the tissue exposures were corrected based on the blood volumes. The tissue VBVs measured by the RBDR method were comparable with the literature values obtained by radiometric measurements.
Collapse
Affiliation(s)
- Hossein Salimi-Moosavi
- Translational Safety & Bioanalytical Sciences, Amgen Research, Amgen Inc, One Amgen Center Drive, Thousand Oaks, California, 91320, USA.
| | - Marcus Soto
- Pharmacokinetics & Drug Metabolism, Amgen Research, Amgen Inc, One Amgen Center Drive, Thousand Oaks, California, 91320, USA
| |
Collapse
|
8
|
Van De Vyver AJ, Walz AC, Heins MS, Abdolzade-Bavil A, Kraft TE, Waldhauer I, Otteneder MB. Investigating brain uptake of a non-targeting monoclonal antibody after intravenous and intracerebroventricular administration. Front Pharmacol 2022; 13:958543. [PMID: 36105215 PMCID: PMC9465605 DOI: 10.3389/fphar.2022.958543] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2022] [Accepted: 08/02/2022] [Indexed: 11/13/2022] Open
Abstract
Monoclonal antibodies play an important role in the treatment of various diseases. However, the development of these drugs against neurological disorders where the drug target is located in the brain is challenging and requires a good understanding of the local drug concentration in the brain. In this original research, we investigated the systemic and local pharmacokinetics in the brain of healthy rats after either intravenous (IV) or intracerebroventricular (ICV) administration of EGFRvIII-T-Cell bispecific (TCB), a bispecific monoclonal antibody. We established an experimental protocol that allows serial sampling in serum, cerebrospinal fluid (CSF) and interstitial fluid (ISF) of the prefrontal cortex in freely moving rats. For detection of drug concentration in ISF, a push-pull microdialysis technique with large pore membranes was applied. Brain uptake into CSF and ISF was characterized and quantified with a reduced brain physiologically-based pharmacokinetic model. The model allowed us to interpret the pharmacokinetic processes of brain uptake after different routes of administration. The proposed model capturing the pharmacokinetics in serum, CSF and ISF of the prefrontal cortex suggests a barrier function between the CSF and ISF that impedes free antibody transfer. This finding suggests that ICV administration may not be better suited to reach higher local drug exposure as compared to IV administration. The model enabled us to quantify the relative contribution of the blood-brain barrier (BBB) and Blood-CSF-Barrier to the uptake into the interstitial fluid of the brain. In addition, we compared the brain uptake of three monoclonal antibodies after IV dosing. In summary, the presented approach can be applied to profile compounds based on their relative uptake in the brain and provides quantitative insights into which pathways are contributing to the net exposure in the brain.
Collapse
Affiliation(s)
- Arthur J. Van De Vyver
- Roche Pharma Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, Basel, Switzerland
| | - Antje-Christine Walz
- Roche Pharma Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, Basel, Switzerland
- *Correspondence: Antje-Christine Walz,
| | | | - Afsaneh Abdolzade-Bavil
- Roche Pharma Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Munich, Penzberg, Germany
| | - Thomas E. Kraft
- Roche Pharma Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Munich, Penzberg, Germany
| | - Inja Waldhauer
- Roche Pharma Research and Early Development (pRED), Roche Innovation Center Zurich (RICZ), Schlieren, Switzerland
| | - Michael B. Otteneder
- Roche Pharma Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, Basel, Switzerland
| |
Collapse
|
9
|
Lin L, Wright MR, Hop CECA, Wong H. Physiologically-Based Pharmacokinetic Models Can be used to Predict the Unique Nonlinear Absorption Profiles of Vismodegib. Drug Metab Dispos 2022; 50:1170-1181. [PMID: 35779865 DOI: 10.1124/dmd.122.000885] [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: 03/02/2022] [Accepted: 06/23/2022] [Indexed: 11/22/2022] Open
Abstract
Predicting human pharmacokinetics (PK) during the drug discovery phase is valuable to assess doses required to reach therapeutic exposures. For orally administered compounds, however, this can be especially difficult since the absorption process is complex. Vismodegib is a compound with unique nonlinear oral PK characteristics in humans. Oral physiologically-based pharmacokinetic (PBPK) models were built using preclinical in vitro and in vivo data and successfully predicted the oral PK profiles in rats, dogs, and monkeys. Simulated drug exposures (AUC0-inf and Cmax), following oral administration were within 2-fold of observed values for the dog and monkey, and close to 2-fold for the rat, providing validation to the model structure. Adaptation of this oral PBPK model to humans, using human physiological parameters coupled with predicted human PK, resulted in underpredictions of vismodegib exposure following both single and multiple doses. When observed human PK was used to drive the oral PBPK model, oral PK profiles in humans were well predicted with fold errors in predicted vs observed drug exposures being close to 1. Importantly, the oral PBPK model captured the unique nonlinear, non-dose dependent PK of vismodegib at steady-state. The mechanism responsible for nonlinearity was consistent with oral absorption being influenced by nonsink permeation conditions. We introduce a new parameter, the permeation gradient factor, to characterize the effect of nonsink conditions on permeation. Using vismodegib as an example, we demonstrate the value of using oral PBPK models in drug discovery to predict the oral PK of compounds with nonlinear absorption characteristics in human. Significance Statement A physiologically-based pharmacokinetic model was built to demonstrate the value of these models early in the drug discovery stage for the prediction of human PK for compounds with unusual oral pharmacokinetics. In this study, our model could successfully capture the unique steady-state oral pharmacokinetics of our model compound, vismodegib. The mechanism for nonlinearity can be attributed to nonsink permeation conditions in vivo. We introduce the permeation gradient factor as a parameter to assess this effect.
Collapse
Affiliation(s)
- Louis Lin
- Faculty of Pharmaceutical Sciences, University of British Columbia, Canada
| | | | | | | |
Collapse
|
10
|
Liu S, Shah DK. Mathematical Models to Characterize the Absorption, Distribution, Metabolism, and Excretion of Protein Therapeutics. Drug Metab Dispos 2022; 50:867-878. [PMID: 35197311 PMCID: PMC11022906 DOI: 10.1124/dmd.121.000460] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2021] [Accepted: 01/31/2022] [Indexed: 11/22/2022] Open
Abstract
Therapeutic proteins (TPs) have ranked among the most important and fastest-growing classes of drugs in the clinic, yet the development of successful TPs is often limited by unsatisfactory efficacy. Understanding pharmacokinetic (PK) characteristics of TPs is key to achieving sufficient and prolonged exposure at the site of action, which is a prerequisite for eliciting desired pharmacological effects. PK modeling represents a powerful tool to investigate factors governing in vivo disposition of TPs. In this mini-review, we discuss many state-of-the-art models that recapitulate critical processes in each of the absorption, distribution, metabolism/catabolism, and excretion pathways of TPs, which can be integrated into the physiologically-based pharmacokinetic framework. Additionally, we provide our perspectives on current opportunities and challenges for evolving the PK models to accelerate the discovery and development of safe and efficacious TPs. SIGNIFICANCE STATEMENT: This minireview provides an overview of mechanistic pharmacokinetic (PK) models developed to characterize absorption, distribution, metabolism, and elimination (ADME) properties of therapeutic proteins (TPs), which can support model-informed discovery and development of TPs. As the next-generation of TPs with diverse physicochemical properties and mechanism-of-action are being developed rapidly, there is an urgent need to better understand the determinants for the ADME of TPs and evolve existing platform PK models to facilitate successful bench-to-bedside translation of these promising drug molecules.
Collapse
Affiliation(s)
- Shufang Liu
- Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, The State University of New York at Buffalo, Buffalo, New York
| | - Dhaval K Shah
- Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, The State University of New York at Buffalo, Buffalo, New York
| |
Collapse
|
11
|
Chang HP, Li Z, Shah DK. Development of a Physiologically-Based Pharmacokinetic Model for Whole-Body Disposition of MMAE Containing Antibody-Drug Conjugate in Mice. Pharm Res 2022; 39:1-24. [PMID: 35044590 DOI: 10.1007/s11095-021-03162-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2021] [Accepted: 12/21/2021] [Indexed: 12/13/2022]
Abstract
PURPOSE To quantitate and mathematically characterize the whole-body pharmacokinetics (PK) of different ADC analytes following administration of an MMAE-conjugated ADC in tumor-bearing mice. METHODS The PK of different ADC analytes (total antibody, total drug, unconjugated drug) was measured following administration of an MMAE-conjugated ADC in tumor-bearing mice. The PK of ADC analytes was compared with the whole-body PK of the antibody and drug obtained following administration of these molecules alone. An ADC PBPK model was developed by linking antibody PBPK model with small-molecule PBPK model, where the drug was assumed to deconjugate in DAR-dependent manner. RESULTS Comparison of antibody biodistribution coefficient (ABC) values for total antibody suggests that conjugation of drug did not significantly affect the PK of antibody. Comparison of tissue:plasma AUC ratio (T/P) for the conjugated drug and total antibody suggests that in certain tissues (e.g., spleen) ADC may demonstrate higher deconjugation. It was observed that the tissue distribution profile of the drug can be altered following its conjugation to antibody. For example, MMAE distribution to the liver was found to increase while its distribution to the heart was found to decrease upon conjugation to antibody. MMAE exposure in the tumor was found to increase by ~20-fold following administration as conjugate (i.e., ADC). The PBPK model was able to a priori predict the PK of all three ADC analytes in plasma, tissues, and tumor reasonably well. CONCLUSIONS The ADC PBPK model developed here serves as a platform for translational and clinical investigations of MMAE containing ADCs.
Collapse
Affiliation(s)
- Hsuan-Ping Chang
- Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, The State University of New York at Buffalo, 455 Pharmacy Building, Buffalo, New York, 14214-8033, USA
| | - Zhe Li
- Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, The State University of New York at Buffalo, 455 Pharmacy Building, Buffalo, New York, 14214-8033, USA
| | - Dhaval K Shah
- Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, The State University of New York at Buffalo, 455 Pharmacy Building, Buffalo, New York, 14214-8033, USA.
| |
Collapse
|
12
|
Norris SA, White H, Tanenbaum A, Williams EL, Cruchaga C, Tian L, Schmidt RE, Perlmutter JS. Severe acute neurotoxicity reflects absolute intra-carotid 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine dose in non-human primates. J Neurosci Methods 2022; 366:109406. [PMID: 34767855 DOI: 10.1016/j.jneumeth.2021.109406] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2021] [Revised: 10/18/2021] [Accepted: 11/02/2021] [Indexed: 11/18/2022]
Affiliation(s)
- S A Norris
- Departments of Neurology, Washington University School of Medicine, 660 S. Euclid Ave, St. Louis, MO 63110, USA; Departments of Radiology, Washington University School of Medicine, 660 S. Euclid Ave, St. Louis, MO 63110, USA.
| | - Hcb White
- Departments of Neurology, Washington University School of Medicine, 660 S. Euclid Ave, St. Louis, MO 63110, USA
| | - A Tanenbaum
- Departments of Neurology, Washington University School of Medicine, 660 S. Euclid Ave, St. Louis, MO 63110, USA
| | - E L Williams
- Departments of Neurology, Washington University School of Medicine, 660 S. Euclid Ave, St. Louis, MO 63110, USA
| | - C Cruchaga
- Departments of Psychiatry, Washington University School of Medicine, 660 S. Euclid Ave, St. Louis, MO 63110, USA
| | - L Tian
- Departments of Neurology, Washington University School of Medicine, 660 S. Euclid Ave, St. Louis, MO 63110, USA
| | - R E Schmidt
- Departments of Pathology, Washington University School of Medicine, 660 S. Euclid Ave, St. Louis, MO 63110, USA
| | - J S Perlmutter
- Departments of Neurology, Washington University School of Medicine, 660 S. Euclid Ave, St. Louis, MO 63110, USA; Departments of Radiology, Washington University School of Medicine, 660 S. Euclid Ave, St. Louis, MO 63110, USA; Departments of Neuroscience, Washington University School of Medicine, 660 S. Euclid Ave, St. Louis, MO 63110, USA; Departments of Physical, Washington University School of Medicine, 660 S. Euclid Ave, St. Louis, MO 63110, USA; Departments of Occupational Therapy, Washington University School of Medicine, 660 S. Euclid Ave, St. Louis, MO 63110, USA
| |
Collapse
|
13
|
Rafidi H, Rajan S, Urban K, Shatz-Binder W, Hui K, Ferl GZ, Kamath AV, Boswell CA. Effect of molecular size on interstitial pharmacokinetics and tissue catabolism of antibodies. MAbs 2022; 14:2085535. [PMID: 35867780 PMCID: PMC9311319 DOI: 10.1080/19420862.2022.2085535] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022] Open
Abstract
Advances in antibody engineering have enabled the construction of novel molecular formats in diverse shapes and sizes, providing new opportunities for biologic therapies and expanding the need to understand how various structural aspects affect their distribution properties. To assess the effect of antibody size on systemic pharmacokinetics (PK) and tissue distribution with or without neonatal Fc receptor (FcRn) binding, we evaluated a series of non-mouse-binding anti-glycoprotein D monoclonal antibody formats, including IgG [~150 kDa], one-armed IgG [~100 kDa], IgG-HAHQ (attenuated FcRn binding) [~150 kDa], F(ab')2 [~100 kDa], and F(ab) [~50 kDa]. Tissue-specific concentration-time profiles were corrected for blood content based on vascular volumes and normalized based on interstitial volumes to allow estimation of interstitial concentrations and interstitial:serum concentration ratios. Blood correction demonstrated that the contribution of circulating antibody on total uptake was greatest at early time points and for highly vascularized tissues. Tissue interstitial PK largely mirrored serum exposure profiles. Similar interstitial:serum ratios were obtained for the two FcRn-binding molecules, IgG and one-armed IgG, which reached pseudo-steady-state kinetics in most tissues. For non-FcRn-binding molecules, interstitial:serum ratios changed over time, suggesting that these molecules did not reach steady-state kinetics during the study. Furthermore, concentration-time profiles of both intact and catabolized molecule were measured by a dual tracer approach, enabling quantification of tissue catabolism and demonstrating that catabolism levels were highest for IgG-HAHQ. Overall, these data sets provide insight into factors affecting preclinical distribution and may be useful in estimating interstitial concentrations and/or catabolism in human tissues.
Collapse
Affiliation(s)
- Hanine Rafidi
- Departments of Preclinical and Translational Pharmacokinetics and Pharmacodynamics, Research and Early Development, Genentech, Inc, South San Francisco, CA, USA
| | - Sharmila Rajan
- Departments of Preclinical and Translational Pharmacokinetics and Pharmacodynamics, Research and Early Development, Genentech, Inc, South San Francisco, CA, USA
| | - Konnie Urban
- Safety Assessment, Research and Early Development, Genentech, Inc, South San Francisco, CA, USA
| | - Whitney Shatz-Binder
- Protein Chemistry, Research and Early Development, Genentech, Inc, South San Francisco, CA, USA
| | - Keliana Hui
- Departments of Preclinical and Translational Pharmacokinetics and Pharmacodynamics, Research and Early Development, Genentech, Inc, South San Francisco, CA, USA
| | - Gregory Z Ferl
- Departments of Preclinical and Translational Pharmacokinetics and Pharmacodynamics, Research and Early Development, Genentech, Inc, South San Francisco, CA, USA.,Biomedical Imaging, Research and Early Development, Genentech, Inc, South San Francisco, CA, USA
| | - Amrita V Kamath
- Departments of Preclinical and Translational Pharmacokinetics and Pharmacodynamics, Research and Early Development, Genentech, Inc, South San Francisco, CA, USA
| | - C Andrew Boswell
- Departments of Preclinical and Translational Pharmacokinetics and Pharmacodynamics, Research and Early Development, Genentech, Inc, South San Francisco, CA, USA.,Biomedical Imaging, Research and Early Development, Genentech, Inc, South San Francisco, CA, USA
| |
Collapse
|
14
|
Pardridge WM. Kinetics of Blood-Brain Barrier Transport of Monoclonal Antibodies Targeting the Insulin Receptor and the Transferrin Receptor. Pharmaceuticals (Basel) 2021; 15:3. [PMID: 35056060 PMCID: PMC8778919 DOI: 10.3390/ph15010003] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2021] [Revised: 12/10/2021] [Accepted: 12/18/2021] [Indexed: 12/21/2022] Open
Abstract
Biologic drugs are large molecule pharmaceuticals that do not cross the blood-brain barrier (BBB), which is formed by the brain capillary endothelium. Biologics can be re-engineered for BBB transport as IgG fusion proteins, where the IgG domain is a monoclonal antibody (MAb) that targets an endogenous BBB transporter, such as the insulin receptor (IR) or transferrin receptor (TfR). The IR and TfR at the BBB transport the receptor-specific MAb in parallel with the transport of the endogenous ligand, insulin or transferrin. The kinetics of BBB transport of insulin or transferrin, or an IRMAb or TfRMAb, can be quantified with separate mathematical models. Mathematical models to estimate the half-time of receptor endocytosis, MAb or ligand exocytosis into brain extracellular space, or receptor recycling back to the endothelial luminal membrane were fit to the brain uptake of a TfRMAb or a IRMAb fusion protein in the Rhesus monkey. Model fits to the data also allow for estimates of the rates of association of the MAb in plasma with the IR or TfR that is embedded within the endothelial luminal membrane in vivo. The parameters generated from the model fits can be used to estimate the brain concentration profile of the MAb over time, and this brain exposure is shown to be a function of the rate of clearance of the antibody fusion protein from the plasma compartment.
Collapse
|
15
|
Pardridge WM, Chou T. Mathematical Models of Blood-Brain Barrier Transport of Monoclonal Antibodies Targeting the Transferrin Receptor and the Insulin Receptor. Pharmaceuticals (Basel) 2021; 14:535. [PMID: 34205013 PMCID: PMC8226686 DOI: 10.3390/ph14060535] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2021] [Revised: 05/29/2021] [Accepted: 06/01/2021] [Indexed: 02/07/2023] Open
Abstract
We develop and analyze mathematical models for receptor-mediated transcytosis of monoclonal antibodies (MAb) targeting the transferrin receptor (TfR) or the insulin receptor (IR), which are expressed at the blood-brain barrier (BBB). The mass-action kinetic model for both the TfR and IR antibodies were solved numerically to generate predictions for the concentrations of all species in all compartments considered. Using these models, we estimated the rates of MAb endocytosis into brain capillary endothelium, which forms the BBB in vivo, the rates of MAb exocytosis from the intra-endothelial compartment into brain extracellular space, and the rates of receptor recycling from the endothelial space back to the luminal endothelial plasma membrane. Our analysis highlights the optimal rates of MAb association with the targeted receptor. An important role of the endogenous ligand, transferrin (Tf) or insulin, in receptor-mediated-transport (RMT) of the associated MAb was found and was attributed to the five order magnitude difference between plasma concentrations of Tf (25,000 nM) and insulin (0.3 nM). Our modeling shows that the very high plasma concentration of Tf leads to only 5% of the endothelial TfR expressed on the luminal endothelial membrane.
Collapse
Affiliation(s)
| | - Tom Chou
- Departments of Computational Medicine and Mathematics, UCLA, Los Angeles, CA 90095, USA;
| |
Collapse
|
16
|
Hesse F, Somai V, Kreis F, Bulat F, Wright AJ, Brindle KM. Monitoring tumor cell death in murine tumor models using deuterium magnetic resonance spectroscopy and spectroscopic imaging. Proc Natl Acad Sci U S A 2021; 118:e2014631118. [PMID: 33727417 PMCID: PMC8000230 DOI: 10.1073/pnas.2014631118] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
2H magnetic resonance spectroscopic imaging has been shown recently to be a viable technique for metabolic imaging in the clinic. We show here that 2H MR spectroscopy and spectroscopic imaging measurements of [2,3-2H2]malate production from [2,3-2H2]fumarate can be used to detect tumor cell death in vivo via the production of labeled malate. Production of [2,3-2H2]malate, following injection of [2,3-2H2]fumarate (1 g/kg) into tumor-bearing mice, was measured in a murine lymphoma (EL4) treated with etoposide, and in human breast (MDA-MB-231) and colorectal (Colo205) xenografts treated with a TRAILR2 agonist, using surface-coil localized 2H MR spectroscopy at 7 T. Malate production was also imaged in EL4 tumors using a fast 2H chemical shift imaging sequence. The malate/fumarate ratio increased from 0.016 ± 0.02 to 0.16 ± 0.14 in EL4 tumors 48 h after drug treatment (P = 0.0024, n = 3), and from 0.019 ± 0.03 to 0.25 ± 0.23 in MDA-MB-231 tumors (P = 0.0001, n = 5) and from 0.016 ± 0.04 to 0.28 ± 0.26 in Colo205 tumors (P = 0.0002, n = 5) 24 h after drug treatment. These increases were correlated with increased levels of cell death measured in excised tumor sections obtained immediately after imaging. 2H MR measurements of [2,3-2H2]malate production from [2,3-2H2]fumarate provide a potentially less expensive and more sensitive method for detecting cell death in vivo than 13C MR measurements of hyperpolarized [1,4-13C2]fumarate metabolism, which have been used previously for this purpose.
Collapse
Affiliation(s)
- Friederike Hesse
- Cancer Research UK Cambridge Institute, University of Cambridge, Cambridge CB2 0RE, United Kingdom
| | - Vencel Somai
- Cancer Research UK Cambridge Institute, University of Cambridge, Cambridge CB2 0RE, United Kingdom
- Department of Radiology, School of Clinical Medicine, University of Cambridge, Cambridge CB2 0SL, United Kingdom
| | - Felix Kreis
- Department of Information Technology and Electrical Engineering, ETH Zürich, 8092 Zürich, Switzerland
| | - Flaviu Bulat
- Cancer Research UK Cambridge Institute, University of Cambridge, Cambridge CB2 0RE, United Kingdom
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Alan J Wright
- Cancer Research UK Cambridge Institute, University of Cambridge, Cambridge CB2 0RE, United Kingdom
| | - Kevin M Brindle
- Cancer Research UK Cambridge Institute, University of Cambridge, Cambridge CB2 0RE, United Kingdom;
- Department of Biochemistry, University of Cambridge, Cambridge CB2 1GA, United Kingdom
| |
Collapse
|
17
|
AT-527, a Double Prodrug of a Guanosine Nucleotide Analog, Is a Potent Inhibitor of SARS-CoV-2 In Vitro and a Promising Oral Antiviral for Treatment of COVID-19. Antimicrob Agents Chemother 2021; 65:AAC.02479-20. [PMID: 33558299 PMCID: PMC8097421 DOI: 10.1128/aac.02479-20] [Citation(s) in RCA: 87] [Impact Index Per Article: 29.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Accepted: 02/03/2021] [Indexed: 01/17/2023] Open
Abstract
The impact of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), the causative agent of COVID-19, is global and unprecedented. Although remdesivir has recently been approved by the FDA to treat SARS-CoV-2 infection, no oral antiviral is available for outpatient treatment. AT-527, an orally administered double prodrug of a guanosine nucleotide analog, was previously shown to be highly efficacious and well tolerated in hepatitis C virus (HCV)-infected subjects. Here, we report the potent in vitro activity of AT-511, the free base of AT-527, against several coronaviruses, including SARS-CoV-2. In normal human airway epithelial cells, the concentration of AT-511 required to inhibit replication of SARS-CoV-2 by 90% (EC90) was 0.47 μM, very similar to its EC90 against human coronavirus (HCoV)-229E, HCoV-OC43, and SARS-CoV in Huh-7 cells. Little to no cytotoxicity was observed for AT-511 at concentrations up to 100 μM. Substantial levels of the active triphosphate metabolite AT-9010 were formed in normal human bronchial and nasal epithelial cells incubated with 10 μM AT-511 (698 ± 15 and 236 ± 14 μM, respectively), with a half-life of at least 38 h. Results from steady-state pharmacokinetic and tissue distribution studies of nonhuman primates administered oral doses of AT-527, as well as pharmacokinetic data from subjects given daily oral doses of AT-527, predict that twice daily oral doses of 550 mg AT-527 will produce AT-9010 trough concentrations in human lung that exceed the EC90 observed for the prodrug against SARS-CoV-2 replication. This suggests that AT-527 may be an effective treatment option for COVID-19.
Collapse
|
18
|
Chang HP, Kim SJ, Shah DK. Whole-Body Pharmacokinetics of Antibody in Mice Determined using Enzyme-Linked Immunosorbent Assay and Derivation of Tissue Interstitial Concentrations. J Pharm Sci 2020; 110:446-457. [PMID: 32502472 DOI: 10.1016/j.xphs.2020.05.025] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2020] [Revised: 05/06/2020] [Accepted: 05/20/2020] [Indexed: 12/20/2022]
Abstract
Here we have reported whole-body disposition of wild-type IgG and FcRn non-binding IgG in mice, determined using Enzyme-Linked Immunosorbent Assay (ELISA). The disposition data generated using ELISA are compared with previously published biodistribution data generated using radiolabelled IgG. In addition, we introduce a novel concept of ABCIS values, which are defined as percentage ratios of tissue interstitial and plasma AUC values. These values can help in predicting tissue interstitial concentrations of monoclonal antibodies (mAbs) based on the plasma concentrations. Tissue interstitial concentrations derived from our study are also compared with previously reported values measured using microdialysis or centrifugation method. Lastly, the new set of biodistribution data generated using ELISA are used to refine the PBPK model for mAbs.
Collapse
Affiliation(s)
- Hsuan-Ping Chang
- Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, The State University of New York at Buffalo, Buffalo, NY, USA
| | - Se Jin Kim
- Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, The State University of New York at Buffalo, Buffalo, NY, USA
| | - Dhaval K Shah
- Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, The State University of New York at Buffalo, Buffalo, NY, USA.
| |
Collapse
|