A Tutorial on Target-Mediated Drug Disposition (TMDD) Models

Structure-Guided Elaboration of a Fragment-Like Hit into an Orally Efficacious Leukotriene A4 Hydrolase Inhibitor

Leukotriene A4 hydrolase (LTA4H) is the final and rate-limiting enzyme in the biosynthesis of pro-inflammatory leukotriene B4 (LTB4). Preclinical studies have provided strong evidence that LTA4H is an attractive drug target for the treatment of chronic inflammatory diseases. Here, we describe the transformation of compound 2, a fragment-like hit, into the potent inhibitor of LTA4H 3. Our strategy involved two key steps. First, we aimed to increase the polarity of fragment 2 to improve its drug-likeness, particularly its solubility, while preserving both its promising potency and low molecular weight. Second, we utilized structural information and incorporated a basic amino function, which allowed for the formation of an essential hydrogen bond with Q136 of LTA4H and consequently enhanced the potency. Compound 3 exhibited exceptional selectivity and showed oral efficacy in a KRN passive serum-induced arthritis model in mice. The anticipated human dose to achieve 90% target engagement at the trough concentration was determined to be 40 mg administered once daily.

Antibody Drug Clearance: An Underexplored Marker of Outcomes with Checkpoint Inhibitors

Immune-checkpoint inhibitor (ICI) therapy has dramatically changed the clinical landscape for several cancers, and ICI use continues to expand across many cancer types. Low baseline clearance (CL) and/or a large reduction of CL during treatment correlates with better clinical response and longer survival. Similar phenomena have also been reported with other monoclonal antibodies (mAb) in cancer and other diseases, highlighting a characteristic of mAb clinical pharmacology that is potentially shared among various mAbs and diseases. Though tempting to attribute poor outcomes to low drug exposure and arguably low target engagement due to high CL, such speculation is not supported by the relatively flat exposure-response relationship of most ICIs, where a higher dose or exposure is not likely to provide additional benefit. Instead, an elevated and/or increasing CL could be a surrogate marker of the inherent resistant phenotype that cannot be reversed by maximizing drug exposure. The mechanisms connecting ICI clearance, therapeutic efficacy, and resistance are unclear and likely to be multifactorial. Therefore, to explore the potential of ICI CL as an early marker for efficacy, this review highlights the similarities and differences of CL characteristics and CL-response relationships for all FDA-approved ICIs, and we compare and contrast these to selected non-ICI mAbs. We also discuss underlying mechanisms that potentially link mAb CL with efficacy and highlight existing knowledge gaps and future directions where more clinical and preclinical investigations are warranted to clearly understand the value of baseline and/or time-varying CL in predicting response to ICI-based therapeutics.

A mechanism-based pathway toward administering highly active N-phage cocktails

Bacteriophage (phage) therapy is being explored as a possible response to the antimicrobial resistance public health emergency. Administering a mixture of different phage types as a cocktail is one proposed strategy for therapeutic applications, but the optimal method for formulating phage cocktails remains a major challenge. Each phage strain has complex pharmacokinetic/pharmacodynamic (PK/PD) properties which depend on the nano-scale size, target-mediated, self-dosing nature of each phage strain, and rapid selection of resistant subpopulations. The objective of this study was to explore the pharmacodynamics (PD) of three unique and clinically relevant anti-Pseudomonas phages after simulation of dynamic dosing strategies. The Hollow Fiber Infection Model (HFIM) is an in vitro system that mimics in vivo pharmacokinetics (PK) with high fidelity, providing an opportunity to quantify phage and bacteria concentration profiles over clinical time scales with rich sampling. Exogenous monotherapy-bolus (producing max concentrations of Cmax = 7 log10 PFU/mL) regimens of phages LUZ19, PYO2, and E215 produced Pseudomonas aeruginosa nadirs of 0, 2.14, or 2.99 log10 CFU/mL after 6 h of treatment, respectively. Exogenous combination therapy bolus regimens (LUZ19 + PYO2 or LUZ19 + E215) resulted in bacterial reduction to <2 log10 CFU/mL. In contrast, monotherapy as a continuous infusion (producing a steady-state concentration of Css,avg = 2 log10PFU/mL) was less effective at reducing bacterial densities. Specifically, PYO2 failed to reduce bacterial density. Next, a mechanism-based mathematical model was developed to describe phage pharmacodynamics, phage-phage competition, and phage-dependent adaptive phage resistance. Monte Carlo simulations supported bolus dose regimens, predicting lower bacterial counts with bolus dosing as compared to prolonged phage infusions. Together, in vitro and in silico evaluation of the time course of phage pharmacodynamics will better guide optimal patterns of administration of individual phages as a cocktail.

Sabatolimab (MBG453) model-informed drug development for dose selection in patients with myelodysplastic syndrome/acute myeloid leukemia and solid tumors

Sabatolimab is a novel immunotherapy with immuno-myeloid activity that targets T-cell immunoglobulin domain and mucin domain-3 (TIM-3) on immune cells and leukemic blasts. It is being evaluated for the treatment of myeloid malignancies in the STIMULUS clinical trial program. The objective of this analysis was to support the sabatolimab dose-regimen selection in hematologic malignancies. A population pharmacokinetic (PopPK) model was fit to patients with solid tumors and hematologic malignancies, which included acute myeloid leukemia, myelodysplastic syndrome (including intermediate-, high-, and very high-risk per Revised International Prognostic Scoring System), and chronic myelomonocytic leukemia. The PopPK model, together with a predictive model of sabatolimab distribution to the bone marrow and binding to TIM-3 was used to predict membrane-bound TIM-3 bone marrow occupancy. In addition, the total soluble TIM-3 (sTIM-3) kinetics and the pharmacokinetic (PK) exposure-response relationship in patients with hematologic malignancies were examined. At intravenous doses above 240 mg Q2w and 800 mg Q4w, we observed linear PK, a plateau in the accumulation of total sTIM-3, and a flat exposure-response relationship for both safety and efficacy. In addition, the model predicted membrane-bound TIM-3 occupancy in the bone marrow was above 95% in over 95% of patients. Therefore, these results support the selection of the 400 mg Q2w and 800 mg Q4w dosing regimens for the STIMULUS clinical trial program.

Preclinical development of 1B7/CD3, a novel anti-TSLPR bispecific antibody that targets CRLF2-rearranged Ph-like B-ALL

Patients harboring CRLF2-rearranged B-lineage acute lymphocytic leukemia (B-ALL) face a 5-year survival rate as low as 20%. While significant gains have been made to position targeted therapies for B-ALL treatment, continued efforts are needed to develop therapeutic options with improved duration of response. Here, first we have demonstrated that patients with CRLF2-rearranged Ph-like ALL harbor elevated thymic stromal lymphopoietin receptor (TSLPR) expression, which is comparable with CD19. Then we present and evaluate the anti-tumor characteristics of 1B7/CD3, a novel CD3-redirecting bispecific antibody (BsAb) that co-targets TSLPR. In vitro, 1B7/CD3 exhibits optimal binding to both human and cynomolgus CD3 and TSLPR. Further, 1B7/CD3 was shown to induce potent T cell activation and tumor lytic activity in both cell lines and primary B-ALL patient samples. Using humanized cell- or patient-derived xenograft models, 1B7/CD3 treatment was shown to trigger dose-dependent tumor remission or growth inhibition across donors as well as induce T cell activation and expansion. Pharmacokinetic studies in murine models revealed 1B7/CD3 to exhibit a prolonged half-life. Finally, toxicology studies using cynomolgus monkeys found that the maximum tolerated dose of 1B7/CD3 was ≤1 mg/kg. Overall, our preclinical data provide the framework for the clinical evaluation of 1B7/CD3 in patients with CRLF2-rearranged B-ALL.

Clinical and research updates on the VISTA immune checkpoint: immuno-oncology themes and highlights

Immune checkpoints limit the activation of the immune system and serve an important homeostatic function but can also restrict immune responses against tumors. Inhibition of specific immune checkpoint proteins such as the B7:CD28 family members programmed cell death protein-1 (PD-1) and cytotoxic T-lymphocyte antigen-4 (CTLA-4) has transformed the treatment of various cancers by promoting the anti-tumor activation of immune cells. In contrast to these effects, the V-domain immunoglobulin suppressor of T-cell activation (VISTA) regulates the steady state of the resting immune system and promotes homeostasis by mechanisms distinct from PD-1 and CTLA-4. The effects of VISTA blockade have been shown to include a decrease in myeloid suppression coupled with proinflammatory changes by mechanisms that are separate and distinct from other immune checkpoint proteins; in some preclinical studies these immune effects appear synergistic. Given the potential benefits of VISTA blockade in the context of cancer therapy, the second Annual VISTA Symposium was convened virtually on September 23, 2022, to review new research from investigators and immuno-oncology experts. Discussions in the meeting extended the knowledge of VISTA biology and the effects of VISTA inhibition, particularly on cells of the myeloid lineage and resting T cells, as three candidate anti-VISTA antibodies are in, or nearing, clinical development.

Relationship Between Cetuximab Target-Mediated Pharmacokinetics and Progression-Free Survival in Metastatic Colorectal Cancer Patients

BACKGROUND AND OBJECTIVE

Cetuximab, an anti-epidermal growth factor receptor (EGFR) monoclonal immunoglobulin (Ig)G1 antibody, has been approved for the treatment of metastatic colorectal cancer (mCRC). The influence of target-antigen on cetuximab pharmacokinetics has never been investigated using target-mediated drug disposition (TMDD) modelling. This study aimed to investigate the relationship between cetuximab concentrations, target kinetics and progression-free survival (PFS).

METHODS

In this ancillary study (NCT00559741), 91 patients with mCRC treated with cetuximab were assessed. Influence of target levels on cetuximab pharmacokinetics was described using TMDD modelling. The relationship between cetuximab concentrations, target kinetics and time-to-progression (TTP) was described using a joint pharmacokinetic-TTP model, where unbound target levels were assumed to influence hazard of progression by an Emax model. Mitigation strategies of concentration-response relationship, i.e., time-varying endogenous clearance and mutual influences of clearance and time-to-progression were investigated.

RESULTS

Cetuximab concentration-time data were satisfactorily described using the TMDD model with quasi-steady-state approximation and time-varying endogenous clearance. Estimated target parameters were baseline target levels (R0 = 43 nM), and complex elimination rate constant (kint = 0.95 day-1). Estimated time-varying clearance parameters were time-invariant component of CL (CL0= 0.38 L/day-1), time-variant component of CL (CL1= 0.058 L/day-1) and first-order rate of CL1 decreasing over time (kdes = 0.049 day-1). Part of concentration-TTP was TTP-driven, where clearance and TTP were inversely correlated. In addition, increased target occupancy was associated with increased TTP.

CONCLUSION

This is the first study describing the complex relationship between cetuximab target-mediated pharmacokinetics and PFS in mCRC patients using a joint PK-time-to-progression model. Further studies are needed to provide a more in-depth description of this relationship.

Prediction of relative change in free nerve growth factor following subcutaneous administration of tanezumab, a novel monoclonal antibody to nerve growth factor

Tanezumab is a monoclonal antibody against nerve growth factor (NGF). We investigated tanezumab pharmacokinetic (PK)-NGF relationships and predicted the extent of systemic free NGF suppression with target-mediated drug disposition (TMDD) modeling using data from three pivotal phase III interventional studies (NCT02697773, NCT02709486, and NCT02528188) in patients with osteoarthritis. Patients received tanezumab 2.5 mg or 5 mg every 8 weeks (q8w) subcutaneously. A TMDD model using a previously established population PK model was used to describe plasma tanezumab and serum total NGF concentration data, and simulations were performed to predict "unobserved" free NGF versus time profiles and dose-response relationships for free NGF. A total of 2992 patients had available data for plasma tanezumab or serum total NGF concentrations and were included in the analysis; 706 of these had data for both tanezumab and total NGF concentrations. The model generally performed well to predict observed total NGF concentrations up to ~24 weeks after each dose. Simulations suggested free NGF concentration would be suppressed by ~75% (median) near the peak of tanezumab concentration and by less than 5% (median) around the trough tanezumab concentration with a tanezumab 2.5 mg q8w regimen. Free NGF concentration was predicted to return to baseline level at ~8 weeks (95% prediction interval: 5-16 weeks) after the last tanezumab dose. This model adequately described plasma tanezumab and serum total NGF concentrations following s.c. administration of tanezumab 2.5 or 5 mg q8w, allowed prediction of relative change in systemic free NGF following s.c. administration of tanezumab.

A 20-Year Research Overview: Quantitative Prediction of Hepatic Clearance Using the In Vitro-In Vivo Extrapolation Approach Based on Physiologically Based Pharmacokinetic Modeling and Extended Clearance Concept

Understanding the extended clearance concept and establishing a physiologically based pharmacokinetic (PBPK) model are crucial for investigating the impact of changes in transporter and metabolizing enzyme abundance/functions on drug pharmacokinetics in blood and tissues. This mini-review provides an overview of the extended clearance concept and a PBPK model that includes transporter-mediated uptake processes in the liver. In general, complete in vitro and in vivo extrapolation (IVIVE) poses challenges due to missing factors that bridge the gap between in vitro and in vivo systems. By considering key in vitro parameters, we can capture in vivo pharmacokinetics, a strategy known as the top-down or middle-out approach. We present the latest progress, theory, and practice of the Cluster Gauss-Newton method, which is used for middle-out analyses. As examples of poor IVIVE, we discuss "albumin-mediated hepatic uptake" and "time-dependent inhibition" of OATP1Bs. The hepatic uptake of highly plasma-bound drugs is more efficient than what can be accounted for by their unbound concentration alone. This phenomenon is referred to as "albumin-mediated" hepatic uptake. IVIVE was improved by measuring hepatic uptake clearance in vitro in the presence of physiologic albumin concentrations. Lastly, we demonstrate the application of Cluster Gauss-Newton method-based analysis to the target-mediated drug disposition of bosentan. Incorporating saturable target binding and OATP1B-mediated hepatic uptake into the PBPK model enables the consideration of nonlinear kinetics across a wide dose range and the prediction of receptor occupancy over time. SIGNIFICANCE STATEMENT: There have been multiple instances where researchers' endeavors to unravel the underlying mechanism of poor in vitro-in vivo extrapolation have led to the discovery of previously undisclosed truths. These include 1) albumin-mediated hepatic uptake, 2) the target-mediated drug disposition in small molecules, and 3) the existence of a trans-inhibition mechanism by inhibitors for OATP1B-mediated hepatic uptake of drugs. Consequently, poor in vitro-in vivo extrapolation and the subsequent inquisitiveness of scientists may serve as a pivotal gateway to uncover hidden mechanisms.

Improving the efficacy of plant-made anti-HIV monoclonal antibodies for clinical use

Introduction

Broadly neutralising antibodies are promising candidates for preventing and treating Human Immunodeficiency Virus/Acquired Immunodeficiency Syndrome (HIV/AIDS), as an alternative to or in combination with antiretroviral therapy (ART). These mAbs bind to sites on the virus essential for virus attachment and entry, thereby inhibiting entry into the host cell. However, the cost and availability of monoclonal antibodies, especially combinations of antibodies, hampers implementation of anti-HIV bNAb therapies in low- to middle- income countries (LMICs) where HIV-1 prevalence is highest.

Methods

We have produced three HIV broadly neutralizing antibodies (bNAbs), 10-1074, VRC01 and 3BNC117 in the Nicotiana benthamiana transient expression system. The impact of specific modifications to enhance potency and efficacy were assessed. To prolong half-life and increase bioavailability, a M252Y/S254T/T256E (YTE) or M428L/N434S (LS) mutation was introduced. To increase antibody dependent cellular cytotoxicity (ADCC), we expressed an afucosylated version of each antibody using a glycoengineered plant line.

Results

The majority of bNAbs and their variants could be expressed at yields of up to 47 mg/kg. Neither the expression system nor the modifications impacted the neutralization potential of the bNAbs. Afucosylated bNAbs exhibit enhanced ability to bind to FcγRIIIa and trigger ADCC, regardless of the presence of Fc amino acid mutations. Lastly, we demonstrated that Fc-modified variants expressed in plants show enhanced binding to FcRn, which results in a favourable in vivo pharmacokinetic profile compared to their unmodified counterparts.

Conclusion

Tobacco plants are suitable expression hosts for anti-HIV bNAbs with increased efficacy and an improved pharmacokinetic profile.

Pharmacokinetics and tolerability of prefilled syringe and auto-injector presentations of MSB11456: results of a randomized, single-dose study in healthy adults

BACKGROUND

Tocilizumab is a monoclonal immunoglobulin G interleukin-6 receptor antagonist. MSB11456 is a proposed tocilizumab biosimilar.

OBJECTIVE

To determine the pharmacokinetic equivalence of a single subcutaneous injection of MSB11456, when delivered via autoinjector (AI) and prefilled syringe (PFS), in healthy adult subjects.

RESEARCH DESIGN AND METHODS

In this randomized, open-label, single fixed-dose, crossover study, 91 subjects received subcutaneous administration of tocilizumab 162 mg via AI and PFS presentations. The primary endpoint pharmacokinetic parameters were analyzed using analysis of variance. Safety data were summarized descriptively.

RESULTS

There were no differences in pharmacokinetic parameters between presentations, and safety parameters were comparable. The 90% confidence intervals for the geometric least squares mean ratios of all primary pharmacokinetic parameters were contained within the predefined 80.00% to 125.00% bioequivalence limits, indicating pharmacokinetic equivalence between the AI and PFS.

CONCLUSIONS

MSB11456 administration via AI was bioequivalent to administration via PFS. MSB11456 can be administered by AI or PFS, increasing the available range of self-injection devices.

TRIAL REGISTRATION

The trial is registered at EudraCT, number 2020-003419-86.

Pharmacokinetics of a proposed tocilizumab biosimilar (MSB11456) versus US-licensed tocilizumab: results of a randomized, double-blind, single-intravenous dose study in healthy adults

BACKGROUND

Tocilizumab, a recombinant monoclonal immunoglobulin G, targets the interleukin-6 receptor. MSB11456 is a proposed tocilizumab biosimilar.

OBJECTIVES

To assess pharmacokinetic equivalence of intravenous MSB11456 to US-licensed tocilizumab.

RESEARCH DESIGN AND METHODS

In this double-blind, parallel-group, single-dose study, 128 healthy adults were randomized to a single one-hour 8 mg/kg IV infusion of either MSB11456 or US-licensed tocilizumab. Blood samples were collected pre-dose and at regular intervals up to day 48 post-dose. The primary endpoint pharmacokinetic parameter was analyzed using analysis of variance (ANOVA) model on the natural logarithm of the endpoint (AUC_0-last), with treatment as a fixed effect. Immunogenicity and safety data were summarized descriptively.

RESULTS

Subjects received either MSB11456 (N = 62) or US-licensed tocilizumab (N = 66). Pharmacokinetic bioequivalence, defined as 90% confidence intervals for the geometric least squares mean ratio entirely contained within the 80.00% to 125.00% equivalence limits, was demonstrated between MSB11456 and US-licensed tocilizumab for the primary and secondary pharmacokinetic endpoints. Anti-drug antibody responses, frequency of neutralizing antibodies against tocilizumab, and safety profiles showed no notable between-treatment differences. Safety was comparable between treatments.

CONCLUSIONS

Pharmacokinetic similarity of MSB11456 and US-licensed tocilizumab was demonstrated, with comparable immunogenicity and safety profiles, supporting MSB11455 as a biosimilar to US-licensed tocilizumab. The trial is registered at EudraCT, number 2019-003484-22.

Best Practices in mAb and Soluble Target Assay Selection for Quantitative Modelling and Qualitative Interpretation

Biologics, especially monoclonal antibodies (mAbs), are an increasingly important part of the drug discovery and development portfolio across the pharmaceutical industry. To enable robust demonstration of pillars 1 and 2 [1] for mAbs, specialised assays are required to measure the complex interactions between mAb and target. This is especially important for the interpretation of soluble target interactions. In some instances, multiple assays with overlapping purposes (e.g., developing both complex and total assays) have been developed. In retrospect, these efforts may have led to excessive time and resources spent in assay development and the generation of data that is contradictory or misleading. Our recommendation is to invest resources early into the development of total assays for both mAb and target. Free target assay data may be inaccurate and report higher levels of free target than are present in the sample at collection due to re-equilibrium during measurement. Total assay formats are inherently less sensitive to the effects of sample preparation, assay conditions, and re-equilibration than free or complex assays. It is acknowledged that pathology/pharmacology is ultimately driven by the free target and knowledge of its dynamics are critical. However, generation of appropriate total target data and using model-based estimation of free target concentrations is a more robust approach than utilisation of direct assay derived estimates. Where free data are utilised, the potential biases should be prospectively considered when developing the assay and utilising the data for quantitative analyses.

A mechanistic PK/PD model to enable dose selection of the potent anti-tryptase antibody (MTPS9579A) in patients with moderate-to-severe asthma

Tryptase, a protease implicated in asthma pathology, is secreted from mast cells upon activation during an inflammatory allergic response. MTPS9579A is a novel monoclonal antibody that inhibits tryptase activity by irreversibly dissociating the active tetramer into inactive monomers. This study assessed the relationship between MTPS9579A concentrations in healthy subjects and tryptase levels in serum and nasal mucosal lining fluid from healthy subjects and patients with moderate-to-severe asthma. These data were used to develop a mechanistic pharmacokinetic/pharmacodynamic (PK/PD) model that quantitatively inter-relates MTPS9579A exposure and inhibition of active tryptase in the airway of patients with asthma. From initial estimates of airway tryptase levels and drug partitioning, the PK/PD model predicted almost complete neutralization of active tryptase in the airway of patients with asthma with MTPS9579A doses of 900 mg and greater, administered intravenously (i.v.) once every 4 weeks (q4w). Suppression of active tryptase during an asthma exacerbation event was also evaluated using the model by simulating the administration of MTPS9579A during a 100-fold increase in tryptase secretion in the local tissue. The PK/PD model predicted that 1800 mg MTPS9579A i.v. q4w results in 95.7% suppression of active tryptase at the steady-state trough concentration. Understanding how the exposure-response relationship of MTPS9579A in healthy subjects translates to patients with asthma is critical for future clinical studies assessing tryptase inhibition in the airway of patients with moderate-to-severe asthma.

Review of the Existing Translational Pharmacokinetics Modeling Approaches Specific to Monoclonal Antibodies (mAbs) to Support the First-In-Human (FIH) Dose Selection

The interest in therapeutic monoclonal antibodies (mAbs) has continuously growing in several diseases. However, their pharmacokinetics (PK) is complex due to their target-mediated drug disposition (TMDD) profiles which can induce a non-linear PK. This point is particularly challenging during the pre-clinical and translational development of a new mAb. This article reviews and describes the existing PK modeling approaches used to translate the mAbs PK from animal to human for intravenous (IV) and subcutaneous (SC) administration routes. Several approaches are presented, from the most empirical models to full physiologically based pharmacokinetic (PBPK) models, with a focus on the population PK methods (compartmental and minimal PBPK models). They include the translational approaches for the linear part of the PK and the TMDD mechanism of mAbs. The objective of this article is to provide an up-to-date overview and future perspectives of the translational PK approaches for mAbs during a model-informed drug development (MIDD), since the field of PK modeling has gained recently significant interest for guiding mAbs drug development.

Application of physiologically-based pharmacokinetic models for therapeutic proteins and other novel modalities

The past two decades have seen diversification of drug development pipelines and approvals from traditional small molecule therapies to alternative modalities including monoclonal antibodies, engineered proteins, antibody drug conjugates (ADCs), oligonucleotides and gene therapies. At the same time, physiologically-based pharmacokinetic (PBPK) models for small molecules have seen increased industry and regulatory acceptance.This review focusses on the current status of the application of PBPK models to these newer modalities and give a perspective on the successes, challenges and future directions of this field.There is greatest experience in the development of PBPK models for therapeutic proteins, and PBPK models for ADCs benefit from prior experience for both therapeutic proteins and small molecules. For other modalities, the application of PBPK models is in its infancy.Challenges are discussed and a common theme is lack of availability of physiological and experimental data to characterise systems and drug parameters to enable a priori prediction of pharmacokinetics. Furthermore, sufficient clinical data are required to build confidence in developed models.The PBPK modelling approach provides a quantitative framework for integrating knowledge and data from multiple sources and can be built on as more data becomes available.

Early Feasibility Assessment: A Method for Accurately Predicting Biotherapeutic Dosing to Inform Early Drug Discovery Decisions

The application of model-informed drug discovery and development (MID3) approaches in the early stages of drug discovery can help determine feasibility of drugging a target, prioritize between targets, or define optimal drug properties for a target product profile (TPP). However, applying MID3 in early discovery can be challenging due to the lack of pharmacokinetic (PK) and pharmacodynamic (PD) data at this stage. Early Feasibility Assessment (EFA) is the application of mechanistic PKPD models, built from first principles, and parameterized by data that is readily available early in drug discovery to make effective dose predictions. This manuscript demonstrates the ability of EFA to make accurate predictions of clinical effective doses for nine approved biotherapeutics and outlines the potential of extending this approach to novel therapeutics to impact early drug discovery decisions.

Assessment of Strategies for Safe Drug Discontinuation and Transition of Denosumab Treatment in PMO—Insights From a Mechanistic PK/PD Model of Bone Turnover

Denosumab (Dmab) treatment against postmenopausal osteoporosis (PMO) has proven very efficient in increasing bone mineral density (BMD) and reducing the risk of bone fractures. However, concerns have been recently raised regarding safety when drug treatment is discontinued. Mechanistic pharmacokinetic-pharmacodynamic (PK-PD) models are the most sophisticated tools to develop patient specific drug treatments of PMO to restore bone mass. However, only a few PK-PD models have addressed the effect of Dmab drug holidays on changes in BMD. We showed that using a standard bone cell population model (BCPM) of bone remodelling it is not possible to account for the spike in osteoclast numbers observed after Dmab discontinuation. We show that inclusion of a variable osteoclast precursor pool in BCPMs is essential to predict the experimentally observed rapid rise in osteoclast numbers and the associated increases in bone resorption. This new model also showed that Dmab withdrawal leads to a rapid increase of damage in the bone matrix, which in turn decreases the local safety factor for fatigue failure. Our simulation results show that changes in BMD strongly depend on Dmab concentration in the central compartment. Consequently, bone weight (BW) might play an important factor in calculating effective Dmab doses. The currently clinically prescribed constant Dmab dose of 60 mg injected every 6 months is less effective in increasing BMD for patients with high BW (2.5% for 80 kg in contrast to 8% for 60 kg after 6 years of treatment). However, bone loss observed 24 months after Dmab withdrawal is less pronounced in patients with high BW (3.5% for 80kg and 8.5% for 60 kg). Finally, we studied how to safely discontinue Dmab treatment by exploring several transitional and combined drug treatment strategies. Our simulation results indicate that using transitional reduced Dmab doses are not effective in reducing rapid bone loss. However, we identify that use of a bisphosphonate (BP) is highly effective in avoiding rapid bone loss and increase in bone tissue damage compared to abrupt withdrawal of Dmab. Furthermore, the final values of BMD and damage were not sensitive to the time of administration of the BP.

Recent advances in the translation of drug metabolism and pharmacokinetics science for drug discovery and development

Drug metabolism and pharmacokinetics (DMPK) is an important branch of pharmaceutical sciences. The nature of ADME (absorption, distribution, metabolism, excretion) and PK (pharmacokinetics) inquiries during drug discovery and development has evolved in recent years from being largely descriptive to seeking a more quantitative and mechanistic understanding of the fate of drug candidates in biological systems. Tremendous progress has been made in the past decade, not only in the characterization of physiochemical properties of drugs that influence their ADME, target organ exposure, and toxicity, but also in the identification of design principles that can minimize drug-drug interaction (DDI) potentials and reduce the attritions. The importance of membrane transporters in drug disposition, efficacy, and safety, as well as the interplay with metabolic processes, has been increasingly recognized. Dramatic increases in investments on new modalities beyond traditional small and large molecule drugs, such as peptides, oligonucleotides, and antibody-drug conjugates, necessitated further innovations in bioanalytical and experimental tools for the characterization of their ADME properties. In this review, we highlight some of the most notable advances in the last decade, and provide future perspectives on potential major breakthroughs and innovations in the translation of DMPK science in various stages of drug discovery and development.

Nonlinear Population Pharmacokinetics of Anifrolumab in Healthy Volunteers and Patients With Systemic Lupus Erythematosus

We characterized the population pharmacokinetics of anifrolumab, a type I interferon receptor–blocking antibody. Pharmacokinetic data were analyzed from the anifrolumab (intravenous [IV], every 4 weeks) arms from 5 clinical trials in patients with systemic lupus erythematosus (SLE) (n = 664) and healthy volunteers (n = 6). Population pharmacokinetic modeling was performed using a 2‐compartment model with parallel linear and nonlinear elimination pathways. The impact of covariates (demographics, interferon gene signature [IFNGS, high/low], disease characteristics, renal/hepatic function, SLE medications, and antidrug antibodies) on pharmacokinetics was evaluated. Time‐varying clearance (CL) was characterized using an empirical sigmoidal time‐dependent function. Anifrolumab exposure increased more than dose‐proportionally from 100 to 1000 mg IV every 4 weeks. Based on population pharmacokinetics modeling, the baseline median linear CL was 0.193 L/day in IFNGS‐high patients and 0.153 L/day in IFNGS‐low/healthy volunteers. After a year, median anifrolumab linear CL decreased by 8.4% from baseline. Body weight and IFNGS were significant pharmacokinetic covariates, whereas age, sex, race, disease activity, SLE medications, and presence of antidrug antibodies had no significant effect on anifrolumab pharmacokinetics. Anifrolumab at a concentration of 300 mg IV every 4 weeks was predicted to be below the lower limit of quantitation in 95% of patients ≈10 weeks after a single dose and ≈16 weeks after stopping dosing at steady state. To conclude, anifrolumab exhibited nonlinear pharmacokinetics and time‐varying linear CL; doses ≥300 mg IV every 4 weeks provided sustained anifrolumab concentrations. This study provides further evidence to support the use of anifrolumab 300 mg IV every 4 weeks in patients with moderate to severe SLE.

Bayesian Nonlinear Models for Repeated Measurement Data: An Overview, Implementation, and Applications

Nonlinear mixed effects models have become a standard platform for analysis when data is in the form of continuous and repeated measurements of subjects from a population of interest, while temporal profiles of subjects commonly follow a nonlinear tendency. While frequentist analysis of nonlinear mixed effects models has a long history, Bayesian analysis of the models has received comparatively little attention until the late 1980s, primarily due to the time-consuming nature of Bayesian computation. Since the early 1990s, Bayesian approaches for the models began to emerge to leverage rapid developments in computing power, and have recently received significant attention due to (1) superiority to quantify the uncertainty of parameter estimation; (2) utility to incorporate prior knowledge into the models; and (3) flexibility to match exactly the increasing complexity of scientific research arising from diverse industrial and academic fields. This review article presents an overview of modeling strategies to implement Bayesian approaches for the nonlinear mixed effects models, ranging from designing a scientific question out of real-life problems to practical computations.

Berkeley Madonna Version 10-A simulation package for solving mathematical models

Berkeley Madonna is a software program that provides an easy and intuitive environment for graphically building and numerically solving mathematical equations. Our users range from college undergraduates with little or no mathematical experience to academic researchers and professionals building and simulating sophisticated mathematical models that represent complex systems in the biological, chemical, and engineering fields. Here we briefly describe our recent advances including a new Java‐based user interface introduced in Version 9 and our transition from a 32‐ to 64‐bit architecture with the release of Version 10. We take the reader through an example tutorial that illustrates how to construct a mathematical model in Berkeley Madonna while highlighting some of the recent changes to the software. Specifically, we construct a standard pharmacokinetic model of the antifungal medication amphotericin B taken from the literature and discuss aspects related to model building, key numerical considerations, data fitting, and graphical visualization. We end by discussing planned functionality and features intended for future releases.

Characterizing the pharmacokinetics and biodistribution of therapeutic proteins: an industry white paper

Characterization of the pharmacokinetics (PK) and biodistribution of therapeutic proteins (TPs) is a hot topic within the pharmaceutical industry, particularly with an ever-increasing catalog of novel modality TPs. Here, we review the current practices, and provide a summary of extensive cross-company discussions as well as a survey completed by International Consortium for Innovation and Quality (IQ consortium) members on this theme. A wide variety of in vitro, in vivo and in silico techniques are currently used to assess PK and biodistribution of TPs, and we discuss the relevance of these from an industry perspective, focusing on PK/PD understanding at the preclinical stage of development, and translation to human. We consider that the 'traditional in vivo biodistribution study' is becoming insufficient as a standalone tool, and thorough characterization of the interaction of the TP with its target(s), target biology, and off-target interactions at a microscopic scale are key to understand the overall biodistribution at a full-body scale. Our summary of the current challenges and our recommendations to address these issues could provide insight into the implementation of best practices in this area of drug development, and continued cross-company collaboration will be of tremendous value. Significance Statement The Innovation & Quality Consortium (IQ) Translational and ADME Sciences Leadership Group (TALG) working group for the ADME of therapeutic proteins evaluates the current practices, recent advances, and challenges in characterizing the PK and biodistribution of therapeutic proteins during drug development, and proposes recommendations to address these issues. Incorporating the in vitro, in vivo and in silico approaches discussed herein may provide a pragmatic framework to increase early understanding of PK/PD relationships, and aid translational modelling for first-in-human dose predictions.

An Extended Model Including Target Turnover, Ligand-Target Complex Kinetics, and Binding Properties to Describe Drug-Receptor Interactions

Since the beginning of this century, target-mediated drug disposition has become a central concept in modeling drug action in drug development. It combines a range of processes, such as turnover, protein binding, internalization, and non-specific elimination, and often serves as a nucleus of more complex pharmacokinetic models. It is simple enough to comprehend but complex enough to be able to describe a wide range of phenomena and data sets. However, the complexity comes at a price: many parameters. In this chapter, we present an overview of the temporal development of the compounds involved after different types of drug doses and offer convenient handles for dissecting data sets in a sophisticated manner in order to estimate the values of these parameters, such as rate constants and pertinent concentrations.

Interpretation of Clinical Efficacy of Biologics in Chronic Rhinosinusitis With Nasal Polyps via Understanding the Local and Systemic Pathomechanisms

Chronic rhinosinusitis with nasal polyps (CRSwNP) is a heterogeneous disease treated with medication or surgery. For recalcitrant type 2 CRSwNP, biological agents have been effectively used to improve nasal polyp score, nasal congestion score, daily symptoms related to CRSwNP, and time to systemic corticosteroid use or revision surgery. Although general guidelines for using biologics to treat CRSwNP were proposed by the European Position Paper on Rhinosinusitis and Nasal Polyps in 2020 and various studies have tested their efficacy, there is much more to learn about biologics—specific indication and choice of biologics based on the endotypes, for instance. Understanding the vascular distribution of monoclonal antibodies and the differences in the vascularity of the non-polyp mucosa and nasal polyp tissue will not only aid understanding of each biologic’s clinical effect but also provide insights to establishing a more personalized approach to treating recalcitrant CRSwNP with biologics.

Population Pharmacokinetic-B Cell Modeling for Ofatumumab in Patients with Relapsing Multiple Sclerosis

BACKGROUND

Ofatumumab, a fully human anti-CD20 monoclonal antibody indicated for the treatment of relapsing forms of multiple sclerosis (RMS), binds to a unique conformational epitope, thereby depleting B cells very efficiently and allowing subcutaneous administration at lower doses.

OBJECTIVES

The aims were to characterize the relationship between ofatumumab concentration and B cell levels, including the effect of covariates such as body weight, age, or baseline B cell count, and use simulations to confirm the chosen therapeutic dose.

METHODS

Graphical and regression analyses previously performed based on data from a dose-range finding study provided the B cell depletion target used in the present work. All available adult phase 2/3 data for ofatumumab in RMS patients were pooled to develop a population pharmacokinetics (PK)-B cell count model, using nonlinear mixed-effects modeling. The population PK-B cell model was used to simulate B cell depletion and repletion times and the effect of covariates on PK and B cell metrics, as well as the dose response across a range of subcutaneous ofatumumab monthly doses.

RESULTS

The final PK-B cell model was developed using data from 1486 patients. The predetermined B cell target was best achieved and sustained with the 20-mg dose regimen, with median B cell count reaching 8 cells/µL in 11 days and negligible repletion between doses. Only weight had a significant effect on PK, which did not translate into any clinically relevant effect on B cell levels.

CONCLUSION

The PK-B cell modeling confirms the dose chosen for the licensed ofatumumab regimen and demonstrates no requirement for dose adjustment based on adult patient characteristics.

Infliximab Efficacy May Be Linked to Full TNF-α Blockade in Peripheral Compartment-A Double Central-Peripheral Target-Mediated Drug Disposition (TMDD) Model

Infliximab is an anti-TNF-α monoclonal antibody approved in chronic inflammatory bowel diseases (IBD). This study aimed at providing an in-depth description of infliximab target-mediated pharmacokinetics in 133 IBD patients treated with 5 mg/kg infliximab at weeks 0, 2, 14, and 22. A two-compartment model with double target-mediated drug disposition (TMDD) in both central and peripheral compartments was developed, using a rich database of 26 ankylosing spondylitis patients as a reference for linear elimination kinetics. Population approach and quasi-steady-state (QSS) approximation were used. Concentration-time data were satisfactorily described using the double-TMDD model. Target-mediated parameters of central and peripheral compartments were respectively baseline TNF concentrations (R^C₀ = 3.3 nM and R^P₀ = 0.46 nM), steady-stated dissociation rates (K^C_SS = 15.4 nM and K^P_SS = 0.49 nM), and first-order elimination rates of complexes (k^C_int = 0.17 day⁻¹ and k^P_int = 0.0079 day⁻¹). This model showed slower turnover of targets and infliximab-TNF complex elimination rate in peripheral compartment than in central compartment. This study allowed a better understanding of the multi-scale target-mediated pharmacokinetics of infliximab. This model could be useful to improve model-based therapeutic drug monitoring of infliximab in IBD patients.

Target-mediated drug disposition modeling of an anti-TFPI antibody (MG1113) in cynomolgus monkeys to predict human pharmacokinetics and pharmacodynamics

BACKGROUND

MG1113 is a human monoclonal antibody of tissue factor pathway inhibitor (TFPI) under development for prophylaxis for hemophilia patients with or without inhibitors against factor VIII products, which have been used for the treatment of hemophilia. Because TFPI is a negative regulator in the extrinsic coagulation pathway, neutralization of TFPI function by MG1113 can potentially increase coagulation activity by bypassing the intrinsic coagulation pathway, which factor VIII activates.

OBJECTIVES

This study aims to determine the correlation between pharmacokinetics (PK) and pharmacodynamics (PD) after administering MG1113 to monkeys and to predict the PK and PD of MG1113 in humans by the Target-Mediated Drug Disposition (TMDD) model using the results from monkeys.

METHODS

The PK profile of MG1113 and the PD effect on the free TFPI level were evaluated after intravenous (IV) and subcutaneous (SC) administrations of MG1113 (2.5, 5, and 10 mg/kg) to male cynomolgus monkeys. After setting up the PK/PD model on monkeys, PK parameters on humans were calculated using allometric scaling, and then clinically effective doses were predicted applying the TMDD model.

RESULTS AND CONCLUSIONS

MG1113 showed nonlinear PK after both IV and SC administrations at the dosing range from 2.5 to10 mg/kg. The concentrations of MG1113 versus TFPI could be characterized a dose-response relationship using a TMDD model. The TMDD modeling and simulation built in this study were used to simulate various dosage regimens of MG1113 to apply to the first-in-human study design, and moreover expected to be referred to establish the dose for further clinical trials.

Modeling Pharmacokinetics and Pharmacodynamics of Therapeutic Antibodies: Progress, Challenges, and Future Directions

With more than 90 approved drugs by 2020, therapeutic antibodies have played a central role in shifting the treatment landscape of many diseases, including autoimmune disorders and cancers. While showing many therapeutic advantages such as long half-life and highly selective actions, therapeutic antibodies still face many outstanding issues associated with their pharmacokinetics (PK) and pharmacodynamics (PD), including high variabilities, low tissue distributions, poorly-defined PK/PD characteristics for novel antibody formats, and high rates of treatment resistance. We have witnessed many successful cases applying PK/PD modeling to answer critical questions in therapeutic antibodies’ development and regulations. These models have yielded substantial insights into antibody PK/PD properties. This review summarized the progress, challenges, and future directions in modeling antibody PK/PD and highlighted the potential of applying mechanistic models addressing the development questions.

Safety and pharmacokinetics of a new biosimilar trastuzumab (HL02): a Phase I bioequivalence study in healthy Chinese men

Objectives: The study aimed to explore the bioequivalence of a proposed biosimilar HL02 vs. its reference products (US-trastuzumab) among healthy Chinese men.Methods: In this study, nine healthy male subjects received single ascending doses of trastuzumab biosimilar (HL02, 2-8 mg/kg), and then a randomized, double-blind, two-arm, parallel study was conducted to investigate the PK similarity of HL02 (6 mg/kg) with that of US-trastuzumab as a reference drug.Results: PK properties exhibited by HL02 (N = 55) were similar to those of US-trastuzumab (N = 52). The comparison of biosimilarity with US-trastuzumab showed that the 90% confidence intervals (CIs) of the ratios for C_max, AUC_0-t, and AUC_0-∞ were within 80-125%. Nineteen subjects were positive for ADA and negative for NAb in both HL02 and US-trastuzumab groups. In total, 81.67% of subjects in HL02 and 78.95% in US-trastuzumab groups showed treatment-related mild or moderate adverse events, mild elevation of transaminase level being the most common adverse events (AE) recorded.Conclusions: The PK characteristics and immunogenicity exhibited by HL02 were similar to that of the reference product, US-trastuzumab. The safety profiles were similar in both the treatment groups with mild-moderate adverse effects.

Real world use of biologic drug levels and anti-drug antibodies in patients with psoriasis - does therapeutic drug monitoring have a place in routine clinical practice?

Background: Psoriasis is a chronic disorder with increasing new treatments targeting the T-helper cell (Th)-1/Th17 axis. There remains a subset of patients who experience a primary or secondary failure to biologic treatments.Methods: We present ten patients with psoriasis who failed biologic therapy with measurement of serum drug levels and anti-drug antibody levels (ADAs) with review of the current literature. Our objective was to identify demographic factors, disease status, drug level and ADAs which might correlate with primary and secondary failure.Results: There are a number of factors affecting drug levels in patients with psoriasis on biologics including the presence of ADAs, patient adherence to treatment regimes, pharmacogenetics and the pharmacokinetic properties of the drug following subcutaneous injection. Our results demonstrate that biologic failure is related to low serum drug levels subtherapeutic in 80% of our cohort. Primary failure may correlate with the presence of ADAs but not with serum drug levels. All patients were ANA negative and there remains considerable debate on the utility of routine ANA testing.Conclusions: The role of therapeutic drug monitoring in dermatology remains uncertain and requires further study. We aim to promote debate in the dermatology community as to the utility of therapeutic drug monitoring in routine practice.

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.

Physiologically based pharmacokinetic (PBPK) modeling of RNAi therapeutics: Opportunities and challenges

Physiologically based pharmacokinetic (PBPK) modeling is a powerful tool with many demonstrated applications in various phases of drug development and regulatory review. RNA interference (RNAi)-based therapeutics are a class of drugs that have unique pharmacokinetic properties and mechanisms of action. With an increasing number of RNAi therapeutics in the pipeline and reaching the market, there is a considerable amount of active research in this area requiring a multidisciplinary approach. The application of PBPK models for RNAi therapeutics is in its infancy and its utility to facilitate the development of this new class of drugs is yet to be fully evaluated. From this perspective, we briefly discuss some of the current computational modeling approaches used in support of efficient development and approval of RNAi therapeutics. Considerations for PBPK model development are highlighted both in a relative context between small molecules and large molecules such as monoclonal antibodies and as it applies to RNAi therapeutics. In addition, the prospects for drawing upon other recognized avenues of PBPK modeling and some of the foreseeable challenges in PBPK model development for these chemical modalities are briefly discussed. Finally, an exploration of the potential application of PBPK model development for RNAi therapeutics is provided. We hope these preliminary thoughts will help initiate a dialogue between scientists in the relevant sectors to examine the value of PBPK modeling for RNAi therapeutics. Such evaluations could help standardize the practice in the future and support appropriate guidance development for strengthening the RNAi therapeutics development program.

A Phase I Clinical Study Comparing the Tolerance, Immunogenicity, and Pharmacokinetics of Proposed Biosimilar BAT1806 and Reference Tocilizumab in Healthy Chinese Men

Objective: The study aimed to explore the bioequivalence of a proposed biosimilar BAT1806 to its reference products marketed in the EU and US (RoActemra-EU and Actemra-US) among healthy Chinese men. The tolerance, immunogenicity, and pharmacokinetics (PK) of the three drugs were also investigated.

Methods: In this randomized, double-blind, single-dose, three-arm, parallel study, a single-dose of 4 mg/kg of the reference products, or the biosimilar was administered to the participants. The participants were followed up for 57 days, and PK, immunogenicity, and tolerance evaluations were completed during this period.

Results: The PK parameters were similar in all three groups: BAT1806 (n = 45), RoActemra-EU (n = 42), and Actemra-US (n = 42). The 90% confidence intervals (CIs) for the ratios of C_max, AUC_0–t and AUC_0–∞ were 86.90–104.41% for BAT1806 vs. RoActemra-EU, 91.70–106.15% for BAT1806 vs Actemra-US, and 90.04–105.53% for Actemra-US vs RoActemra-EU. For all comparisons, the 90% CIs for the C_max, AUC_0–t, and AUC_0–∞ were within the predefined bioequivalence limit of 80.00–125.00%. The intersubject variability ranged from 14.5% to 21.5%, which was considerably low. Among the participants, 19 (42.2%), 10 (23.8%), and 12 (28.6%) from the BAT1806, RoActemra-EU, and Actemra-US groups were, respectively, found to be positive for anti-drug antibodies, while 14 (31.1%), nine (21.4%), and 12 (28.6%) were positive for neutralizing antibodies. Nevertheless, these antibodies did not affect the drug concentrations, and the outcomes in the bioequivalence tests were similar after sensitivity analysis. Treatment-related and treatment-emergent adverse events (TEAEs) were recorded in 27, 34, and 32 participants in the BAT1806, RoActemra-EU, and Actemra-US groups, respectively. The most common treatment-related adverse events observed were a decrease in neutrophil, and white blood cell counts.

Conclusion: The PK characteristics of BAT1806 were similar to those of the reference products, RoActemra-EU and Actemra-US. Both BAT1806 and the reference products exhibited low intersubject variability and similar safety profiles.

Clinical trial registration number:http://www.chinadrugtrials.org.cn/index.html, CTR20180039; https://clinicaltrials.gov/NCT03606876

Development of a Pharmacokinetic Model Describing Neonatal Fc Receptor-Mediated Recycling of HL2351, a Novel Hybrid Fc-Fused Interleukin-1 Receptor Antagonist, to Optimize Dosage Regimen

HL2351 (hIL‐1Ra‐hyFc) is a novel recombinant protein formed by the fusion of two human interleukin‐1 receptor antagonist components into one antibody‐derived fragment crystallizable portion. Although HL2351 has a pharmacological mechanism of action similar to that of anakinra as a commercialized biopharmaceutical drug, HL2351 has been desired to reduce the dose frequency and improve therapeutic efficacy due to its long circulation half‐life. In this study, we aimed to develop a population pharmacokinetic (PK) model for HL2351 using a neonatal Fc receptor (FcRn)–mediated recycling model based on a quasi‐steady‐state approximation of target‐mediated drug disposition (TMDD) for the description of interactions between the drug and FcRn. FcRn recycling was expected in the case of HL2351 because of PK related to the antibody portion. A TMDD model was also applied to describe interactions of IL1R with HL2351 or anakinra. PK data were collected from a phase I study conducted in six groups (1, 2, 4, 8, 12 mg/kg HL2351 and 100 mg anakinra single subcutaneous administration; n = 8 per group). In consequence, the PK of anakinra and HL2351 following administration of multiple doses at different dosages were simulated. Optimized doses were considered based on average concentrations of IL1R bound to anakinra and HL2351. HL2351 at doses of 326 mg or 4.267, 4.982, 5.288, 5.458, or 5.748 mg/kg once weekly or HL2351 at 1726 mg or 21.92, 26.86, 29.10, 30.36, or 32.53 mg/kg once biweekly would have similar therapeutic effects with anakinra at a dose of 100 mg or 1, 2, 3, 4, or 8 mg/kg administered once daily, respectively.

Population pharmacokinetics and pharmacodynamics of a novel vascular adhesion protein-1 inhibitor using a multiple-target mediated drug disposition model

ASP8232 is a novel inhibitor of vascular adhesion protein-1 that was under evaluation for reducing residual albuminuria in patients with diabetic kidney disease. To characterize the pharmacokinetics (PK) of ASP8232 and its effect on vascular adhesion protein 1 (VAP-1) plasma activity and VAP-1 concentrations (pharmacodynamics, PD) in an integrated and quantitative manner, a target mediated drug disposition model was developed based on pooled data from four completed clinical trials with ASP8232 in healthy volunteers, and in patients with diabetic kidney disease and diabetic macular edema, respectively. In this model, the binding of ASP8232 to its soluble and membrane-bound target in the central and peripheral compartments were included. The model was able to adequately describe the non-linear PK and PD of ASP8232. The observed difference in PK between healthy volunteers and renally impaired patients could be explained by an effect of baseline estimated glomerular filtration rate on ASP8232 clearance and relative bioavailability. The relationship between ASP8232 concentration and VAP-1 inhibition was successfully established and can be applied to simulate drug exposure and degree of VAP-1 inhibition for any given dose of ASP8232 across the spectrum of renal function.

Population Pharmacokinetics of Sarilumab in Patients with Rheumatoid Arthritis

Background and Objective

Sarilumab binds to the interleukin-6 receptor with high affinity, inhibiting cis and trans signaling by interleukin-6. Sarilumab has demonstrated efficacy and safety in patients with rheumatoid arthritis. The objective of this study was to develop a population-pharmacokinetic model using data from 1770 patients with rheumatoid arthritis across phase I–III studies.

Methods

Potential covariates were identified using a stepwise forward-addition and backward-deletion strategy, and the final model was evaluated by visual predictive check and bootstrap methods.

Results

Sarilumab pharmacokinetics is described by a two-compartment model with first-order absorption and parallel linear and nonlinear Michaelis–Menten elimination. A subcutaneous dose of sarilumab 200 mg every 2 weeks resulted in more pronounced saturation of the nonlinear clearance pathway over the dosing interval than 150 mg every 2 weeks. Steady-state exposure (area under the plasma concentration–time curve from day 0 to day 14) increased twofold with dose escalation from 150 to 200 mg every 2 weeks. Body weight, anti-drug antibody status, sarilumab drug product, sex, creatinine clearance, albumin, and baseline C-reactive protein levels were identified as significant covariates according to the predefined statistical significance criteria in stepwise covariate searches. The main intrinsic source of pharmacokinetic variability in exposure was body weight. Compared with a typical 71-kg patient, the area under the plasma concentration–time curve from day 0 to day 14 was 20–23% lower for an 83-kg patient and 20–25% higher for a 62-kg patient.

Conclusions

These findings, combined with the safety and efficacy data, indicated limited clinical relevance of body-weight effect on sarilumab exposure. No adjustment in sarilumab dose is required for body weight or any other demographics assessed.

Electronic supplementary material

The online version of this article (10.1007/s40262-019-00765-1) contains supplementary material, which is available to authorized users.

Hemopexin as an Inhibitor of Hemolysis-Induced Complement Activation

Hemopexin is the main plasmatic scavenger of cell-free heme, released in the context of intravascular hemolysis or major cell injury. Heme is indispensable for the oxygen transport by hemoglobin but when released outside of the erythrocytes it becomes a danger-associated molecular pattern, contributing to tissue injury. One of the mechanisms of pro-inflammatory action of heme is to activate the innate immune complement cascade. Therefore, we hypothesized that injection of hemopexin will prevent hemolysis-induced complement activation. Human plasma-derived hemopexin is compatible with the heme clearance machinery of the mice. 100 or 500 mg/kg of hemopexin was injected in C57Bl/6 mice before treatment with phenylhydrazine (inducer of erythrocytes lysis) or with PBS as a control. Blood was taken at different timepoints to determine the pharmacokinetic of injected hemopexin in presence and absence of hemolysis. Complement activation was determined in plasma, by the C3 cleavage (western blot) and in the kidneys (immunofluorescence). Kidney injury was evaluated by urea and creatinine in plasma and renal NGAL and HO-1 gene expression were measured. The pharmacokinetic properties of hemopexin (mass spectrometry) in the hemolytic mice were affected by the target-mediated drug disposition phenomenon due to the high affinity of binding of hemopexin to heme. Hemolysis induced complement overactivation and signs of mild renal dysfunction at 6 h, which were prevented by hemopexin, except for the NGAL upregulation. The heme-degrading capacity of the kidney, measured by the HO-1 expression, was not affected by the treatment. These results encourage further studies of hemopexin as a therapeutic agent in models of diseases with heme overload.

Pharmacokinetics and Pharmacodynamics of Garetosmab (Anti-Activin A): Results From a First-in-Human Phase 1 Study

We describe outcomes from the first‐in‐human study of garetosmab (a fully human monoclonal antibody that inhibits activin A) under development for the treatment of fibrodysplasia ossificans progressiva (FOP). In a double‐blind, placebo‐controlled phase 1 study, 40 healthy women of nonchildbearing potential were randomized to receive a single dose of intravenous garetosmab 0.3, 1, 3, or 10 mg/kg; subcutaneous garetosmab 300 mg; or placebo. Serum concentrations of functional garetosmab (with ≥1 arm free to bind to target), total activin A, and antidrug antibodies were measured predose and up to 113 days post–first dose. Garetosmab demonstrated an acceptable safety profile with no dose‐limiting toxicities. Garetosmab displayed nonlinear pharmacokinetics with target‐mediated elimination. With increasing doses of intravenous garetosmab, mean peak concentration increased in a dose‐proportional manner; mean steady‐state estimates ranged from 41.4 to 47.8 mL/kg. A greater than dose‐proportional increase in mean area under the concentration‐time curve from time zero extrapolated to infinity (range, 72.2‐7520 mg*day/L) was observed, consistent with decreasing mean clearance (range, 4.35‐1.34 mL/day/kg). Following administration of intravenous garetosmab, mean concentrations of total activin A increased in a dose‐dependent manner. At 10 mg/kg, total activin A levels reached a state of little or no change between weeks 4 and 12, suggesting saturation of the target‐mediated pathway. No safety signals were seen in this study to preclude investigation in patients. Following intravenous administration, garetosmab concentrations decreased quickly, then decreased over time (reflecting linear elimination), and finally decreased in a nonlinear phase, reflecting target‐mediated elimination. Results here support further investigation. Garetosmab 10 mg/kg every 4 weeks intravenously is being evaluated in patients with FOP (NCT03188666).

Pharmacokinetics, Pharmacodynamics, Safety, and Tolerability of Dupilumab in Healthy Adult Subjects

Dupilumab is a fully human monoclonal antibody directed against the interleukin (IL)‐4 receptor α subunit (IL‐4Rα) of IL‐4 heterodimeric type I and type II receptors that mediate IL‐4/IL‐13 signaling through this pathway. Blockade of these receptors broadly suppresses type 2 inflammation associated with atopic/allergic diseases, including atopic dermatitis and asthma. Six phase 1 studies investigated the pharmacokinetics, pharmacodynamics, safety, and tolerability of dupilumab in healthy subjects. Two randomized, double‐blind, placebo‐controlled, sequential studies assessed safety and tolerability of single escalating dupilumab doses administered intravenously or subcutaneously (one included various racial groups, and one included exclusively Japanese subjects); 3 randomized, parallel‐group, single‐dose studies compared the pharmacokinetic profiles of different dupilumab products and formulations after single subcutaneous doses; and one study assessed dupilumab administered as fast versus slow subcutaneous injections. Dupilumab concentrations in serum were measured in all studies, and total immunoglobulin E (IgE) and thymus‐ and activation‐regulated chemokine (TARC) concentrations were measured in 2 studies as pharmacodynamic markers. Across the phase 1 studies, dupilumab exhibited target‐mediated pharmacokinetics consisting of parallel linear and nonlinear elimination, with the target‐mediated phase highly dominated by nonlinearity at lower drug concentrations. Systemic exposure and tolerability of dupilumab were consistent irrespective of differences in product, formulation, or racial background. Dupilumab reduced circulating concentrations of total IgE and TARC, indicating blockade of IL‐4Rα–mediated signaling. Dupilumab had a favorable safety profile across the wide range of doses administered. Together, these findings support the continued development and use of dupilumab in treatment of type 2 diseases.

Safety, pharmacokinetics and pharmacodynamics of branebrutinib (BMS-986195), a covalent, irreversible inhibitor of Bruton's tyrosine kinase: Randomised phase I, placebo-controlled trial in healthy participants

Aims

Branebrutinib (BMS‐986195) is a potent, highly selective, oral, small‐molecule, covalent inhibitor of Bruton's tyrosine kinase (BTK). This study evaluated safety, pharmacokinetics and pharmacodynamics of branebrutinib in healthy participants.

Methods

This double‐blind, placebo‐controlled, single‐ and multiple‐ascending dose (SAD; MAD) Phase I study (NCT02705989) enrolled participants into 3 parts: SAD, MAD and JMAD (MAD in first‐generation Japanese participants). In each part, participants were randomised 3:1 to receive branebrutinib (SAD: 0.3–30 mg; [J]MAD: 0.3–10 mg) or placebo. Participants in the MAD parts received branebrutinib daily for 14 days and were followed for 14 days postdosing. Safety was assessed by monitoring, laboratory and physical examinations, vital signs, and recording adverse events (AEs). Pharmacodynamics were assessed with a mass spectrometry assay that measured drug‐occupied and free BTK.

Results

The SAD, MAD and JMAD parts of the study included 40, 32 and 24 participants. Branebrutinib was well tolerated and AEs were mild/moderate, except for 1 serious AE that led to discontinuation. Branebrutinib was rapidly absorbed, with maximum plasma concentration occurring within 1 hour and a half‐life of 1.2—1.7 hours, dropping to undetectable levels within 24 hours. BTK occupancy was rapid, with 100% occupancy reached after a single 10‐mg dose. BTK occupancy decayed predictably over time (mean half‐life in MAD panels: 115–154 hours), such that pharmacodynamic effects were maintained after branebrutinib plasma levels fell below the lower limit of quantification.

Conclusion

Rapid and high occupancy of BTK and the lack of notable safety findings support further clinical development of branebrutinib.

Assessment of romosozumab efficacy in the treatment of postmenopausal osteoporosis: Results from a mechanistic PK-PD mechanostat model of bone remodeling

This paper introduces a theoretical framework for the study of the efficacy of romosozumab, a humanized monoclonal antibody targeting sclerostin for the treatment of osteoporosis. We developed a comprehensive mechanistic pharmacokinetic-pharmacodynamic (PK-PD) model of the effect of drug treatment on bone remodeling in postmenopausal osteoporosis (PMO). We utilized a one-compartment PK model to represent subcutaneous injections of romosozumab and subsequent absorption into serum. The PD model is based on a recently-developed bone cell population model describing the bone remodeling process at the tissue scale. The latter accounts for mechanical feedback by incorporating nitric oxide (NO) and sclerostin (Scl) as biochemical feedback molecules. Utilizing a competitive binding model, where Wnt and Scl compete for binding to LRP5/6, allows to regulate anabolic bone remodeling responses. Here, we extended this model with respect to romosozumab binding to sclerostin. For the currently approved monthly injections of 210 mg, the model predicted a 6.59%, 10.38% and 15.25% increase in BMD at the lumbar spine after 6, 12 and 24 months, respectively. These results are in good agreement with the data reported in the literature. Our model is also able to distinguish the bone-site specific drug effects. For instance, at the femoral neck, our model predicts a BMD increase of 3.85% after 12 months of 210 mg injections, which is consistent with literature observations. Finally, our simulations indicate rapid bone loss after treatment discontinuation, indicating that some additional interventions such as use of bisphosphonates are required to maintain bone mass.

Simultaneous Target-Mediated Drug Disposition Model for Two Small-Molecule Compounds Competing for Their Pharmacological Target: Soluble Epoxide Hydrolase

1-(1-propanoylpiperidin-4-yl)-3-[4-(trifluoromethoxy)phenyl]urea (TPPU) and 1-(4-trifluoro-methoxy-phenyl)-3-(1-cyclopropanecarbonyl-piperidin-4-yl)-urea (TCPU) are potent inhibitors of soluble epoxide hydrolase (sEH) that have much better efficacy in relieving nociceptive response than the Food and Drug Administration-approved drug gabapentin in a rodent model of diabetic neuropathy. Experiments conducted in sEH knockout mice or with coadministration of a potent sEH displacer demonstrated that the pharmacokinetics of TPPU and TCPU were influenced by the specific binding to their pharmacologic target sEH, a phenomenon known as target-mediated drug disposition (TMDD). To quantitatively characterize the complex pharmacokinetics of TPPU and TCPU and gain better understanding on their target occupancy, population pharmacokinetics analysis using a nonlinear mixed-effect modeling approach was performed in the current study. The final model was a novel simultaneous TMDD interaction model, in which TPPU and TCPU compete for sEH, with TCPU binding to an additional unknown target pool with larger capacity that we refer to as a refractory pool. The total amount of sEH enzyme in mice was predicted to be 16.2 nmol, which is consistent with the experimental value of 10 nmol. The dissociate rate constants of TPPU and TCPU were predicted to be 2.24 and 2.67 hours^-1, respectively, which is close to the values obtained from in vitro experiments. Our simulation result predicted that 90% of the sEH will be occupied shortly after a low dose of 0.3 mg/kg TPPU administration, with ≥40% of sEH remaining to be bound with TPPU for at least 7 days. Further efficacy experiments are warranted to confirm the predicted target occupancy. SIGNIFICANCE STATEMENT: Although target-mediated drug disposition (TMDD) models have been well documented, most of them were established in a single compound scenario. Our novel model represents the first TMDD interaction model for two small-molecule compounds competing for the same pharmacological target.

Clinical Pharmacokinetic and Pharmacodynamic Considerations in the Treatment of Ulcerative Colitis

Ulcerative colitis (UC) is an inflammatory bowel disease (IBD) of unknown etiology, probably caused by a combination of genetic and environmental factors. The treatment of patients with active UC depends on the severity, localization and history of IBD medication. According to the classic step-up approach, treatment with 5-aminosalicylic acid compounds is the first step in the treatment of mild to moderately active UC. Corticosteroids, such as prednisolone are used in UC patients with moderate to severe disease activity, but only for remission induction therapy because of side effects associated with long-term use. Thiopurines are the next step in the treatment of active UC but monotherapy during induction therapy in UC patients is not preferred because of their slow onset. Therapeutic drug monitoring (TDM) of the pharmacologically active metabolites of thiopurines, 6-thioguanine nucleotide (6-TGN), has proven to be beneficial. Thiopurine S-methyltransferase (TMPT) plays a role in the metabolic conversion pathway of thiopurines and exhibits genetic polymorphism; however, the clinical benefit and relevance of TPMT genotyping is not well established. In patients with severely active UC refractory to corticosteroids, calcineurin inhibitors such as ciclosporin A (CsA) and tacrolimus are potential therapeutic options. These agents usually have a rather rapid onset of action. Monoclonal antibodies (anti-tumor necrosis factor [TNF] agents, vedolizumab) are the last pharmacotherapeutic option for UC patients before surgery becomes inevitable. Body weight, albumin status and antidrug antibodies contribute to the variability in the pharmacokinetics of anti-TNF agents. Additionally, the use of concomitant immunomodulators (thiopurines/methotrexate) lowers the rate of immunogenicity, and therefore the concomitant use of anti-TNF therapy with an immunomodulator may confer some advantage compared with monotherapy in certain patients. TDM of anti-TNF agents could be beneficial in patients with primary nonresponse and secondary loss of response. The potential benefit of applying TDM during vedolizumab treatment has yet to be determined.

Modeling and Simulation of the Pharmacokinetics and Target Engagement of an Antagonist Monoclonal Antibody to Interferon-γ-Induced Protein 10, BMS-986184, in Healthy Participants to Guide Therapeutic Dosing

BMS‐986184 is a human, second‐generation, anti–interferon‐γ–induced protein 10 (IP‐10) monoclonal antibody. In this study the pharmacokinetics and target engagement (TE) of BMS‐986184 in healthy participants were characterized using population‐based target‐mediated drug disposition (TMDD) modeling and data from a first‐in‐human study (NCT02864264). The results of the first‐in‐human study and the model generated were used to conduct stochastic simulations of a virtual population of healthy participants to predict pharmacokinetic exposures and TE responses for different dosage regimens. A 2‐compartment, 2‐target, TMDD structural model, assuming quasi‐steady‐state and stimulated production on treatment, was developed by simultaneous fitting of the total drug, serum‐free IP‐10, and serum total IP‐10 concentration data, with the second unobservable target contribution to drug elimination described by the Michaelis‐Menten elimination term. Model evaluation confirmed agreement between model predictions and observed data. Simulation of a virtual population of healthy individuals demonstrated that steady state was reached at the eighth dosing interval, and that around 150 mg subcutaneously every other week could be a suitable target dosage regimen for future clinical trials. Integrated modeling strategies such as this can be used to help guide rational clinical trial development of drugs with TMDD, leading to improved dose selection and greater patient benefits.

Pharmacokinetics and pharmacodynamic effect of crenezumab on plasma and cerebrospinal fluid beta-amyloid in patients with mild-to-moderate Alzheimer's disease

Background

Crenezumab, a fully humanized anti-beta-amyloid (Aβ) immunoglobulin G4 (IgG4) monoclonal antibody, binds to both monomeric and aggregated forms of Aβ. We assessed the pharmacokinetics (PK)/pharmacodynamics (PD) of crenezumab and its interaction with monomeric Aβ(1–40) and Aβ(1–42) peptides in serum/plasma and cerebrospinal fluid (CSF) samples from the phase II ABBY and BLAZE studies and the phase Ib GN29632 study.

Methods

In ABBY, BLAZE, and GN29632 studies, patients with mild-to-moderate AD were treated with either placebo or crenezumab (300 mg subcutaneously every 2 weeks [q2w], or 15 mg/kg, 30 mg/kg, 45 mg/kg, 60 mg/kg, or 120 mg/kg intravenously q4w). Serum/plasma PK/PD analyses included samples from 131 patients who received crenezumab in all three studies. CSF PK/PD analyses included samples from 76 patients who received crenezumab in ABBY or BLAZE. The impact of baseline patient factors on Aβ profiles was also evaluated.

Results

The serum concentration of crenezumab increased in a dose-proportional manner between 15 and 120 mg/kg q4w. Total monomeric plasma Aβ(1–40) and Aβ(1–42) levels significantly increased after crenezumab administration. The mean crenezumab CSF to serum ratio was ~ 0.3% and was similar across dosing cohorts/routes of administration. No clear correlation was observed between crenezumab concentration and Aβ(1–42) increase in CSF at week 69. The target-mediated drug disposition (TMDD) model described the observed plasma concentration–time profiles of crenezumab and Aβ well. Elimination clearance (CL_el) and central volume of distribution (V_cent) of crenezumab were estimated at 0.159 L/day and 2.89 L, respectively, corresponding to a half-life of ~ 20 days. Subcutaneous bioavailability was estimated at 66.2%.

Conclusions

Crenezumab PK was dose proportional up to 120 mg/kg, with a half-life consistent with IgG monoclonal antibodies. Our findings provide evidence for peripheral target engagement in patients with mild-to-moderate AD. The study also showed that a model-based approach is useful in making inference on PK/PD relationship with unmeasured species such as free plasma Aβ levels.

Trial registrations

ABBY: ClinicalTrials.gov, NCT01343966. Registered April 28, 2011. BLAZE: ClinicalTrials.gov, NCT01397578. Registered July 19, 2011. GN29632: ClinicalTrials.gov, NCT02353598. Registered February 3, 2015.

Dose Correction for a Michaelis-Menten Approximation of a Target-Mediated Drug Disposition Model with a Multiple Intravenous Dosing Regimens

This study aimed to develop a method for implementing dose correction in a Michaelis-Menten (M-M) approximation of a target-mediated drug disposition (TMDD) model with multiple intravenous (IV) bolus administrations. We derived the formula of a correction factor (Fcorr) for each dose in a multiple IV bolus dosing regimens for M-M model. Fcorr depends on the residual free drug amount prior IV bolus dosing event and dose amount. We conducted a stochastic simulation and estimation (SSE) exercise based on therapeutic antibody PK parameters to evaluate the effect of Fcorr on parameter estimation. Previously published clinical PK data of recombinant human erythropoietin (rHuEPO) from four clinical trials in healthy subjects receiving multiple IV bolus doses were analyzed by both M-M model with and without dose correction (MMC and MMNC) as well as the rapid-binding/quasi-steady-state (RB/QSS) TMDD models. Our results showed that MMNC introduced bias to fixed-effect parameter estimates and overestimated random-effect variables. Compared with MMC, MMNC was not able to adequately characterize the nonlinearity in the PK data of antibody and rHuEPO. The MMC-based simulation demonstrated that thricely weekly 10 IU/kg rHuEPO dosing regimen yielded Fcorr = 0.5. This result suggested that the lower-than-expected exposure for rHuEPO at low dose is due to target binding. An M-M approximation of the TMDD model should include a dose correction to avoid model misfitting and potential bias in the estimated PK parameters.

Concept of Pharmacologic Target-Mediated Drug Disposition in Large-Molecule and Small-Molecule Compounds

Target-mediated drug disposition (TMDD) is a term to describe a nonlinear pharmacokinetic (PK) phenomenon that is caused by high-affinity binding of a compound to its pharmacologic targets. As the interaction between a drug and its pharmacologic target belongs to the process of pharmacodynamics (PD), TMDD can be viewed as a consequence of "PD affecting PK." Although there are numerous TMDD-related articles in the literature, most of them focus on characterizing TMDD using various mathematical models, which may not be suitable for those readers who have little interest in mathematical modeling and only want to have an understanding of the basic concept. The goal of this review is to serve as a "primer" on TMDD. This review explains (1) how TMDD happens; (2) why large-molecule and small-molecule compounds exhibiting TMDD demonstrate substantially different nonlinear PK behaviors; (3) what nonlinear PK profiles look like in large-molecule and small-molecule compounds exhibiting TMDD, using pegfilgrastim, erythropoietin, ABT-384, and linagliptin as case examples; and (4) how to identify whether the nonlinear PK of a compound is because of TMDD.