Simulations of site-specific target-mediated pharmacokinetic models for guiding the development of bispecific antibodies.

Bispecific antibodies (BAbs) are novel constructs that are under development and show promise as new therapeutic modalities for cancer and autoimmune disorders. The aim of this study is to develop a semi-mechanistic modeling approach to elucidate the disposition of BAbs in plasma and possible sites of action in humans. Here we present two case studies that showcase the use of modeling to guide BAb development. In case one, a BAb is directed against a soluble and a membrane-bound ligand for treating systemic lupus erythematosus, and in case two, a BAb targets two soluble ligands as a potential treatment for ulcerative colitis and asthma. Model simulations revealed important differences between plasma and tissues, when evaluated for drug disposition and target suppression. Target concentrations at tissue sites and type (soluble vs membrane-bound), tissue-site binding, and binding affinity are all major determinants of BAb disposition and subsequently target suppression. For the presented case studies, higher doses and/or frequent dosing regimens are required to achieve 80 % target suppression in site specific tissue (the more relevant matrix) as compared to plasma. Site-specific target-mediated models may serve to guide the selection of first-in-human doses for new BAbs.

Ubiquity: a framework for physiological/mechanism-based pharmacokinetic/pharmacodynamic model development and deployment.

Practitioners of pharmacokinetic/pharmacodynamic modeling routinely employ various software packages that enable them to fit differential equation based mechanistic or empirical models to biological/pharmacological data. The availability and choice of different analytical tools, while enabling, can also pose a significant challenge in terms of both, implementation and transferability. A package has been developed that addresses these issues by creating a simple text-based format, which provides methods to reduce coding complexity and enables the modeler to describe the components of the model based on the underlying physiochemical processes. A Perl script builds the system for multiple formats (ADAPT, MATLAB, Berkeley Madonna, etc.), enabling analysis across several software packages and reducing the chance for transcription error. Workflows can then be built around this package, which can increase efficiency and model availability. As a proof of concept, tools are included that allow models constructed in this format to be run with MATLAB both at the scripting level and through a generic graphical application that can be compiled and run as a stand-alone application.

Exposure-Response Analyses for Belzutifan to Inform Dosing Considerations and Labeling.

Belzutifan (Welireg, Merck & Co., Inc., Rahway, NJ, USA) is an oral, potent hypoxia-inducible factor-2α inhibitor, recently approved in the United States for the treatment of von Hippel-Lindau (VHL) disease-associated renal cell carcinoma (RCC) and other VHL disease-associated neoplasms. Safety and efficacy were investigated in two clinical studies: a Phase 1 dose escalation/expansion study in solid tumors and RCC and a Phase 2 study in VHL-RCC. A population pharmacokinetic model was used to estimate belzutifan exposures to facilitate exposure-response (E-R) analyses for efficacy and safety endpoints. Relationships between exposure and efficacy (overall response rate, disease control rate, progression-free survival, best overall tumor size response, and other endpoints), safety outcomes (Grade ≥3 anemia, Grade ≥3 hypoxia, and time to first dose reduction/dose interruption), and pharmacodynamic biomarkers (erythropoietin [EPO] and hemoglobin [Hgb]) were evaluated using various regression techniques and time-to-event analyses. Efficacy E-R was generally flat with non-significant positive trends with exposure. The safety E-R analyses demonstrated a lack of relationship for Grade ≥3 hypoxia and a positive relationship for Grade ≥3 anemia, with incidences also significantly dependent on baseline Hgb. Exposure-dependent reductions in EPO and Hgb were observed. Based on the cumulative benefit-risk assessment in VHL disease-associated neoplasms using E-R, no a priori dose adjustment is recommended for any subpopulation. These analyses supported the benefit-risk profile of belzutifan 120 mg once daily dosing in patients with VHL-RCC for labeling and the overall development program.

Population pharmacokinetic analyses for belzutifan to inform dosing considerations and labeling.

Belzutifan (Welireg, Merck & Co., Inc., Rahway, NJ, USA) is an oral, potent inhibitor of hypoxia-inducible factor 2α, approved for the treatment of certain patients with von Hippel-Lindau (VHL) disease-associated renal cell carcinoma (RCC), central nervous system hemangioblastomas, and pancreatic neuroendocrine tumors. It is primarily metabolized by the polymorphic uridine 5'-diphospho-glucuronosyltransferase (UGT) 2B17 and cytochrome (CYP) 2C19. A population pharmacokinetic (PK) model was built, using NONMEM version 7.3, based on demographics/PK data from three clinical pharmacology (food effect, formulation bridging, and genotype/race effect) and two clinical studies (phase I dose escalation/expansion in patients with RCC and other solid tumors; phase II in patients with VHL). Median (range) age for the combined studies was 55 years (19-84) and body weight was 73.6 kg (42.1-165.8). Belzutifan plasma PK was well-characterized by a linear two-compartment model with first-order absorption and elimination. For patients with VHL, the predicted geometric mean (% coefficient of variation) apparent clearance was 7.3 L/h (51%), apparent total volume of distribution was 130 L (35%), and half-life was 12.39 h (42%). There were no clinically relevant differences in belzutifan PK based on the individual covariates of age, sex, ethnicity, race, body weight, mild/moderate renal impairment, or mild hepatic impairment. In this model, dual UGT2B17 and CYP2C19 poor metabolizers (PMs) were estimated to have a 3.2-fold higher area under the plasma concentration-time curve compared to UGT2B17 extensive metabolizer and CYP2C19 non-PM patients. This population PK analysis enabled an integrated assessment of PK characteristics with covariate effects in the overall population and subpopulations for belzutifan labeling.

First-in-Class Anti-immunoglobulin-like Transcript 4 Myeloid-Specific Antibody MK-4830 Abrogates a PD-1 Resistance Mechanism in Patients with Advanced Solid Tumors.

PURPOSE: In this first-in-human study (NCT03564691) in advanced solid tumors, we investigated a novel first-in-class human IgG4 monoclonal antibody targeting the immunoglobulin-like transcript 4 (ILT4) receptor, MK-4830, as monotherapy and in combination with pembrolizumab. PATIENTS AND METHODS: Patients with histologically/cytologically confirmed advanced solid tumors, measurable disease by RECIST v1.1, and evaluable baseline tumor sample received escalating doses of intravenous MK-4830 every 3 weeks as monotherapy (parts A and B) and in combination with pembrolizumab (part C). Safety and tolerability were the primary objectives. Pharmacokinetics, objective response rate per RECIST v1.1, and molecular biomarkers were also evaluated. RESULTS: Of 84 patients, 50 received monotherapy and 34 received combination therapy. No dose-limiting toxicities were observed; maximum tolerated dose was not reached. MK-4830 showed dose-related target engagement. Eleven of 34 patients in the dose-escalation phase who received combination therapy achieved objective responses; 5 previously had progressive disease on anti-PD-1/PD-L1 therapies. Exploratory evaluation of the association between response and pretreatment gene expression related to interferon-gamma signaling in tumors suggested higher sensitivity to T-cell inflammation with combination therapy than historically expected with pembrolizumab monotherapy, with greater response at more moderate levels of inflammation. CONCLUSIONS: This first-in-class MK-4830 antibody dosed as monotherapy and in combination with pembrolizumab was well tolerated with no unexpected toxicities, and demonstrated dose-related evidence of target engagement and antitumor activity. Inflammation intrinsic to the ILT4 mechanism may be facilitated by alleviating the myeloid-suppressive components of the tumor microenvironment, supporting the target of ILT4 as a potential novel immunotherapy in combination with an anti-PD-1/PD-L1 agent.

Pivotal Dose of Pembrolizumab: A Dose-Finding Strategy for Immuno-Oncology.

Despite numerous publications emphasizing the value of dose finding, drug development in oncology is dominated by the mindset that higher dose provides higher efficacy. Examples of dose finding implemented by biopharmaceutical firms can change this mindset. The purpose of this article is to outline a pragmatic dose selection strategy for immuno-oncology (IO) and other targeted monoclonal antibodies (mAbs). The approach was implemented for pembrolizumab. Selecting a recommended phase II dose (RP2D) with a novel mechanism of action is often challenging due to uncertain relationships between pharmacodynamics measurements and clinical end points. Additionally, phase I efficacy and safety data are generally inadequate for RP2D selection for IO mAbs. Here, the RP2D was estimated based on phase I (clinical study KN001 A and A2) pharmacokinetics data as the dose required for target saturation, which represents a surrogate for maximal pharmacological effect for antagonist mAbs. Due to limitations associated with collecting and analyzing tumor biopsies, characterizing intratumoral target engagement (TE) is challenging. To overcome this gap, a physiologically-based pharmacokinetic model was implemented to predict intratumoral TE. As tumors are spatially heterogeneous, TE was predicted in well-vascularized and poorly vascularized tumor regions. Additionally, impact of differences in target expression, for example, due to interindividual variability and cancer type, was simulated. Simulations showed that 200 mg every 3 weeks can achieve ≥ 90% TE in clinically relevant scenarios, resulting in the recommendation of 200 mg every 3 weeks as the RP2D. Randomized dose comparison studies (KN001 B2 and D) showing similar efficacy over a fivefold dose/exposure range confirmed the RP2D as the pivotal dose.

Type I interferon receptor is a primary regulator of target-mediated drug disposition of interferon-beta in mice.

The purpose of this study is to evaluate the primary mechanism through which interferon (IFN)-beta exhibits target-mediated drug disposition (TMDD) and whether the theoretical assumptions of TMDD models are consistent with experimental pharmacokinetic (PK) data. Recombinant murine IFN-beta was administered as an intravenous injection at two dose levels (0.5 and 1 million IU/kg) to male wild-type (WT) and type-I IFN-alpha/beta receptor subunit (IFNAR-1) knockout (KO) mice (A129S7/SvEvBrd strain). Sampling was conducted at various times (n = 3/time point), and plasma was analyzed for IFN-beta concentrations using a validated enzyme-linked immunosorbent assay. The pharmacodynamic (PD) biomarker was IP-10 mRNA that was isolated from the distal femur bone and quantified using reverse transcription-polymerase chain reaction. An integrated model that includes rapid-binding TMDD and an indirect mechanism of drug action was used to characterize the PK/PD profiles. For an experimental control, PK profiles of recombinant murine erythropoietin (muEPO), another drug that exhibits TMDD, were determined after a single intravenous dose (0.5 microg/kg) in WT and KO animals. The concentration-time profiles for IFN-beta differed substantially at initial times for the WT and KO mice at the same dose levels. These differences are characteristic of ligands exhibiting receptor-mediated disposition and were well described by a rapid-binding TMDD model. No differences in muEPO PK were observed in the control study. In summary, the intact IFNAR receptor is a primary regulator of in vivo IFN-beta exposure. An integrated PK/PD model was successfully used to assess the receptor-mediated disposition and dynamics of IFN-beta.

Interspecies scaling of receptor-mediated pharmacokinetics and pharmacodynamics of type I interferons.

PURPOSE: To develop an integrated mechanism-based modeling approach for the interspecies scaling of pharmacokinetic (PK) and pharmacodynamic (PD) properties of type I interferons (IFNs) that exhibit target-mediated drug disposition (TMDD). METHODS: PK and PD profiles of human IFN-beta1a, IFN-beta1b, and IFN-alpha2a in humans, monkeys, rats, and mice from nine studies were extracted from the literature by digitization. Concentration-time profiles from different species were fitted simultaneously using various allometric relationships to scale model-specific parameters. RESULTS: PK/PD profiles of IFN-beta1a in humans and monkeys were successfully characterized by utilizing the same rate constant parameters and scaling the volume of the central compartment to body weight using an allometric exponent of 1. Concentration and effect profiles of other IFNs were also well described by changing only the affinity of the drug to its receptor. PK profiles in rodents were simulated using an allometric exponent of -0.25 for the first-order elimination rate constant, and no receptor-binding was included given the lack of cross-reactivity. CONCLUSIONS: An integrated TMDD PK/PD model was successfully combined with classic allometric scaling techniques and showed good predictive performance. Several parameters obtained from one IFN can be effectively shared to predict the kinetic behavior of other IFN subtypes.

Mechanism-based pharmacokinetic/pharmacodynamic model of parathyroid hormone-calcium homeostasis in rats and humans.

The purpose of this study was to develop a mechanism-based pharmacokinetic/pharmacodynamic model that describes the regulation of the parathyroid hormone (PTH)-Ca(2+) system in rats and humans. Temporal concentration data for endogenous PTH and Ca(2+) were extracted from literature for rats (normal adult males) and humans. In addition, exogenous PTH was administered subcutaneously to male Sprague-Dawley rats with jugular vein catheters, and plasma concentrations were measured over time. A mathematical model was developed and fitted simultaneously to endogenous PTH, Ca(2+), and exogenous PTH concentrations in rats. Ca(2+) concentrations were described using a turnover model, with its depletion being induced by a chelating agent, and PTH concentrations were characterized using a precursor-dependent indirect response model. The same structural model was used for fitting data obtained in humans. PTH stimulation was driven by occupancy of the Ca(2+) sensing receptor, and lowering of physiological Ca(2+) concentrations increased PTH secretion, with PTH profiles being adequately described by the model. PTH stimulatory capacity was baseline-dependent in rats [S(max_rats) = 34.8 x PTH(0)] and humans [S(max_humans) = 392/PTH(0)]. Modeling results suggest that normal rats are twice as sensitive to Ca(2+)-induced PTH stimulation compared with humans. In conclusion, the developed model adequately characterizes the PTH-Ca(2+) regulation across species and may be useful in the development of therapeutic drugs targeting this system.

Optimal Affinity of a Monoclonal Antibody: Guiding Principles Using Mechanistic Modeling.

Affinity optimization of monoclonal antibodies (mAbs) is essential for developing drug candidates with the highest likelihood of clinical success; however, a quantitative approach for setting affinity requirements is often lacking. In this study, we computationally analyzed the in vivo mAb-target binding kinetics to delineate general principles for defining optimal equilibrium dissociation constant ([Formula: see text]) of mAbs against soluble and membrane-bound targets. Our analysis shows that in general [Formula: see text] to achieve 90% coverage for a soluble target is one tenth of its baseline concentration ([Formula: see text]), and is independent of the dosing interval, target turnover rate or the presence of competing ligands. For membrane-bound internalizing targets, it is equal to the ratio of internalization rate of mAb-target complex and association rate constant ([Formula: see text]). In cases where soluble and membrane-bound forms of the target co-exist, [Formula: see text] lies within a range determined by the internalization rate ([Formula: see text]) of the mAb-membrane target complex and the ratio of baseline concentrations of soluble and membrane-bound forms ([Formula: see text]). Finally, to demonstrate practical application of these general rules, we collected target expression and turnover data to project [Formula: see text] for a number of marketed mAbs against soluble (TNFα, RANKL, and VEGF) and membrane-bound targets (CD20, EGFR, and HER2).

Novel Anti-TM4SF1 Antibody-Drug Conjugates with Activity against Tumor Cells and Tumor Vasculature.

Antibody-drug conjugates (ADC) represent a promising therapeutic modality for managing cancer. Here, we report a novel humanized ADC that targets the tetraspanin-like protein TM4SF1. TM4SF1 is highly expressed on the plasma membranes of many human cancer cells and also on the endothelial cells lining tumor blood vessels. TM4SF1 is internalized upon interaction with antibodies. We hypothesized that an ADC against TM4SF1 would inhibit cancer growth directly by killing cancer cells and indirectly by attacking the tumor vasculature. We generated a humanized anti-human TM4SF1 monoclonal antibody, v1.10, and armed it with an auristatin cytotoxic agent LP2 (chemical name mc-3377). v1.10-LP2 selectively killed cultured human tumor cell lines and human endothelial cells that express TM4SF1. Acting as a single agent, v1.10-LP2 induced complete regression of several TM4SF1-expressing tumor xenografts in nude mice, including non-small cell lung cancer and pancreas, prostate, and colon cancers. As v1.10 did not react with mouse TM4SF1, it could not target the mouse tumor vasculature. Therefore, we generated a surrogate anti-mouse TM4SF1 antibody, 2A7A, and conjugated it to LP2. At 3 mpk, 2A7A-LP2 regressed several tumor xenografts without noticeable toxicity. Combination therapy with v1.10-LP2 and 2A7A-LP2 together was more effective than either ADC alone. These data provide proof-of-concept that TM4SF1-targeting ADCs have potential as anticancer agents with dual action against tumor cells and the tumor vasculature. Such agents could offer exceptional therapeutic value and warrant further investigation. Mol Cancer Ther; 14(8); 1868-76. ©2015 AACR.

Pharmacodynamic model of parathyroid hormone modulation by a negative allosteric modulator of the calcium-sensing receptor.

In this study, a pharmacodynamic model is developed, based on calcium-parathyroid hormone (PTH) homeostasis, which describes the concentration-effect relationship of a negative allosteric modulator of the calcium-sensing receptor (CaR) in rats. Plasma concentrations of drug and PTH were determined from plasma samples obtained via serial jugular vein sampling following single subcutaneous doses of 1, 5, 45, and 150 mg/kg to male Sprague-Dawley rats (n = 5/dose). Drug pharmacokinetics was described by a one-compartment model with first-order absorption and linear elimination. Concentration-time profiles of PTH were characterized using a model in which the compound allosterically modulates Ca(+2) binding to the CaR that, in turn, modulates PTH through a precursor-pool indirect response model. Additionally, negative feedback was incorporated to account for tolerance observed at higher dose levels. Model fitting and parameter estimation were conducted using the maximum likelihood algorithm. The proposed model well characterized the data and provided compound specific estimates of the K(i) and cooperativity constant (α) of 1.47 ng/mL and 0.406, respectively. In addition, the estimated model parameters for PTH turnover were comparable to that previously reported. The final generalized model is capable of characterizing both PTH-Ca(+2) homeostasis and the pharmacokinetics and pharmacodynamics associated with the negative allosteric CaR modulator. As such, the model provides a simple platform for analysis of drugs targeting the PTH-Ca(+2) system.

Mechanisms of interferon-beta effects on bone homeostasis.

Restoration of dysregulated bone homeostasis is a therapeutic goal in many diseases including osteoporosis, rheumatoid arthritis and metastatic cancer. The molecular pathways regulating bone remodeling are major therapeutic targets, and studies continue to reveal endogenous factors that may be pathologically up- or down-regulated and lead to an uncoupling of bone formation and resorption. The purpose of this commentary is to highlight new mechanisms of bone homeostatic regulation mediated through the induction of endogenous interferon-beta (IFN-beta). The receptor activator of nuclear factor-kappaB (RANK) ligand (RANKL) is an important factor in the bone resorption cascade, and the RANK-RANKL interaction has been shown to induce IFN-beta and osteoclastogenesis via induction of the c-fos gene. Subsequent binding of IFN-beta to its biological receptor initiates a signal transduction cascade through the classic JAK/STAT pathway, causing an inhibition of c-fos protein production and osteoclast proliferation and differentiation (negative feedback). Another mechanism pertinent to the anti-resorptive effect of IFN-beta is the induction of nitric oxide which has been shown to inhibit osteoclast formation. The role of IFN-beta in bone metabolism could warrant its systematic evaluation as a potential adjunct to therapeutic regimens of osteolytic diseases. Here we also provide discussion of the potential challenges to optimizing IFN-beta pharmacotherapy for such purposes.

Partial derivative-based sensitivity analysis of models describing target-mediated drug disposition.

Sensitivity analysis is commonly used to characterize the effects of parameter perturbations on model output. One use for the approach is the optimization of an experimental design enabling estimation of model parameters with improved accuracy. The primary objective of this study is to conduct a sensitivity analysis of selected target-mediated pharmacokinetic models, ascertain the effect of parameter variations on model predictions, and identify influential model parameters. One linear model (Model 1, control) and 2 target-mediated models (Models 2 and 3) were evaluated over a range of dose levels. Simulations were conducted with model parameters being perturbed at the higher and lower ends from literature mean values. Profiles of free plasma drug concentrations and their partial derivatives with respect to each parameter vs time were analyzed. Perturbations resulted in altered outputs, the extent of which reflected parameter influence. The model outputs were highly sensitive to perturbations of linear disposition parameters in all 3 models. The equilibrium dissociation constant (K(D)) was less influential in Model 2 but was influential in the terminal phase in Model 3, highlighting the role of K(D) in this region. An equation for Model 3 in support of the result for K(D) was derived. Changes in the initial receptor concentration [R(tot) (0)] paralleled the observed effects of initial plasma volume (V(c)) perturbations, with increased influence at higher values. Model 3 was also sensitive to the rates of receptor degradation and internalization. These results suggest that informed sampling may be essential to accurately estimate influential parameters of target-mediated models.