iBCS: 3. A Biopharmaceutics Classification System for Orally Inhaled Drug Products.

In this article, we specify for the first time a quantitative biopharmaceutics classification system for orally inhaled drugs. To date, orally inhaled drug product developers have lacked a biopharmaceutics classification system like the one developed to navigate the development of immediate release of oral medicines. Guideposts for respiratory drug discovery chemists and inhalation product formulators have been elusive and difficult to identify due to the complexity of pulmonary physiology, the intricacies of drug deposition and disposition in the lungs, and the influence of the inhalation delivery device used to deliver the drug as a respirable aerosol. The development of an inhalation biopharmaceutics classification system (iBCS) was an initiative supported by the Product Quality Research Institute (PQRI). The goal of the PQRI iBCS working group was to generate a qualitative biopharmaceutics classification system that can be utilized by inhalation scientists as a "rule of thumb" to identify desirable molecular properties and recognize and manage CMC product development risks based on physicochemical properties of the drug and the deposited lung dose. Herein, we define the iBCS classes quantitatively according to the dose number and permeability. The proposed iBCS was evaluated for its ability to categorize marketed inhaled drugs using data from the literature. The appropriateness of the classification of each drug was assessed based on published development, clinical and nonclinical data, and mechanistic physiologically based biopharmaceutics modeling. The inhaled drug product development challenges for each iBCS classification are discussed and illustrated for different classes of marketed inhaled drugs. Finally, it is recognized that discriminatory laboratory methods to characterize regional lung deposition, dissolution, and permeability will be key to fully realizing the benefits of an iBCS to streamline and derisk inhaled drug development.

A combined physiologically-based pharmacokinetic and quantitative systems pharmacology model for modeling amyloid aggregation in Alzheimer's disease.

Antibody-mediated removal of aggregated ß-amyloid (Aß) is the current, most clinically advanced potential disease-modifying treatment approach for Alzheimer's disease. We describe a quantitative systems pharmacology (QSP) approach of the dynamics of Aß monomers, oligomers, protofibrils, and plaque using a detailed microscopic model of Aß40 and Aß42 aggregation and clearance of aggregated Aß by activated microglia cells, which is enhanced by the interaction of antibody-bound Aß. The model allows for the prediction of Aß positron emission tomography (PET) imaging load as measured by a standardized uptake value ratio. A physiology-based pharmacokinetic model is seamlessly integrated to describe target exposure of monoclonal antibodies and simulate dynamics of cerebrospinal fluid (CSF) and plasma biomarkers, including CSF Aß42 and plasma Aß42 /Aß40 ratio biomarkers. Apolipoprotein E genotype is implemented as a difference in microglia clearance. By incorporating antibody-bound, plaque-mediated macrophage activation in the perivascular compartment, the model also predicts the incidence of amyloid-related imaging abnormalities with edema (ARIA-E). The QSP platform is calibrated with pharmacological and clinical information on aducanumab, bapineuzumab, crenezumab, gantenerumab, lecanemab, and solanezumab, predicting adequately the change in PET imaging measured amyloid load and the changes in the plasma Aß42 /Aß40 ratio while slightly overestimating the change in CSF Aß42 . ARIA-E is well predicted for all antibodies except bapineuzumab. This QSP model could support the clinical trial design of different amyloid-modulating interventions, define optimal titration and maintenance schedules, and provide a first step to understand the variability of biomarker response in clinical practice.

Current practices for QSP model assessment: an IQ consortium survey.

Quantitative Systems Pharmacology (QSP) modeling is increasingly applied in the pharmaceutical industry to influence decision making across a wide range of stages from early discovery to clinical development to post-marketing activities. Development of standards for how these models are constructed, assessed, and communicated is of active interest to the modeling community and regulators but is complicated by the wide variability in the structures and intended uses of the underlying models and the diverse expertise of QSP modelers. With this in mind, the IQ Consortium conducted a survey across the pharmaceutical/biotech industry to understand current practices for QSP modeling. This article presents the survey results and provides insights into current practices and methods used by QSP practitioners based on model type and the intended use at various stages of drug development. The survey also highlights key areas for future development including better integration with statistical methods, standardization of approaches towards virtual populations, and increased use of QSP models for late-stage clinical development and regulatory submissions.

Hallmarks of neurodegenerative disease: A systems pharmacology perspective.

Age-related central neurodegenerative diseases, such as Alzheimer's and Parkinson's disease, are a rising public health concern and have been plagued by repeated drug development failures. The complex nature and poor mechanistic understanding of the etiology of neurodegenerative diseases has hindered the discovery and development of effective disease-modifying therapeutics. Quantitative systems pharmacology models of neurodegeneration diseases may be useful tools to enhance the understanding of pharmacological intervention strategies and to reduce drug attrition rates. Due to the similarities in pathophysiological mechanisms across neurodegenerative diseases, especially at the cellular and molecular levels, we envision the possibility of structural components that are conserved across models of neurodegenerative diseases. Conserved structural submodels can be viewed as building blocks that are pieced together alongside unique disease components to construct quantitative systems pharmacology (QSP) models of neurodegenerative diseases. Model parameterization would likely be different between the different types of neurodegenerative diseases as well as individual patients. Formulating our mechanistic understanding of neurodegenerative pathophysiology as a mathematical model could aid in the identification and prioritization of drug targets and combinatorial treatment strategies, evaluate the role of patient characteristics on disease progression and therapeutic response, and serve as a central repository of knowledge. Here, we provide a background on neurodegenerative diseases, highlight hallmarks of neurodegeneration, and summarize previous QSP models of neurodegenerative diseases.

iBCS: 1. Principles and Framework of an Inhalation-Based Biopharmaceutics Classification System.

For oral drugs, the formulator and discovery chemist have a tool available to them that can be used to navigate the risks associated with the selection and development of immediate release oral drugs and drug products. This tool is the biopharmaceutics classification system (giBCS). Unfortunately, no such classification system exists for inhaled drugs. The perspective outlined in this manuscript provides the foundational principles and framework for a classification system for inhaled drugs. The proposed classification system, an inhalation-based biopharmaceutics classification system (iBCS), is based on fundamental biopharmaceutics principles adapted to an inhalation route of administration framework. It is envisioned that a classification system for orally inhaled drugs will facilitate an understanding of the technical challenges associated with the development of new chemical entities and their associated new drug products (device and drug formulation combinations). Similar to the giBCS, the iBCS will be based on key attributes describing the drug substance (solubility and permeability) and the drug product (dose and dissolution). This manuscript provides the foundational aspects of an iBCS, including the proposed scientific principles and framework upon which such a system can be developed.

iBCS: 2. Mechanistic Modeling of Pulmonary Availability of Inhaled Drugs versus Critical Product Attributes.

This work is the second in a series of publications outlining the fundamental principles and proposed design of a biopharmaceutics classifications system for inhaled drugs and drug products (the iBCS). Here, a mechanistic computer-based model has been used to explore the sensitivity of the primary biopharmaceutics functional output parameters: (i) pulmonary fraction dose absorbed (Fabs) and (ii) drug half-life in lumen (t1/2) to biopharmaceutics-relevant input attributes including dose number (Do) and effective permeability (Peff). Results show the nonlinear sensitivity of primary functional outputs to variations in these attributes. Drugs with Do < 1 and Peff > 1 × 10-6 cm/s show rapid (t1/2 < 20 min) and complete (Fabs > 85%) absorption from lung lumen into lung tissue. At Do > 1, dissolution becomes a critical drug product attribute and Fabs becomes dependent on regional lung deposition. The input attributes used here, Do and Peff, thus enabled the classification of inhaled drugs into parameter spaces with distinctly different biopharmaceutic risks. The implications of these findings with respect to the design of an inhalation-based biopharmaceutics classification system (iBCS) and to the need for experimental methodologies to classify drugs need to be further explored.

Characterizing Pharmacokinetic-Pharmacodynamic Relationships and Efficacy of PI3Kδ Inhibitors in Respiratory Models of TH2 and TH1 Inflammation.

We leveraged a clinical pharmacokinetic (PK)/pharmacodynamics (PD)/efficacy relationship established with an oral phosphatidylinositol 3-kinase (PI3K)δ inhibitor (Idelalisib) in a nasal allergen challenge study to determine whether a comparable PK/PD/efficacy relationship with PI3Kδ inhibitors was observed in preclinical respiratory models of type 2 T helper cell (TH2) and type 1 T helper cell (TH1) inflammation. Results from an in vitro rat blood basophil (CD63) activation assay were used as a PD biomarker. IC50 values for PI3Kδ inhibitors, MSD-496486311, MSD-126796721, Idelalisib, and Duvelisib, were 1.2, 4.8, 0.8, and 0.5 µM. In the ovalbumin Brown Norway TH2 pulmonary inflammation model, all PI3Kδ inhibitors produced a dose-dependent inhibition of bronchoalveolar lavage eosinophils (maximum effect between 80% and 99%). In a follow-up experiment designed to investigate PK attributes [maximum (or peak) plasma concentration (Cmax), area under the curve (AUC), time on target (ToT)] that govern PI3Kδ efficacy, MSD-496486311 [3 mg/kg every day (QD) and 100 mg/kg QD] produced 16% and 93% inhibition of eosinophils, whereas doses (20 mg/kg QD, 10 mg/kg twice per day, and 3 mg/kg three times per day) produced 54% to 66% inhibition. Our profiling suggests that impact of PI3Kδ inhibitors on eosinophils is supported by a PK target with a ToT over the course of treatment close to the PD IC50 rather than strictly driven by AUC, Cmax, or Cmin (minimum blood plasma concentration) coverage. Additional studies in an Altenaria alternata rat model, a sheep Ascaris-sensitive sheep model, and a TH1-driven rat ozone exposure model did not challenge our hypothesis, suggesting that an IC50 level of TE (target engagement) sustained for 24 hours is required to produce efficacy in these traditional models. We conclude that the PK/PD observations in our animal models appear to align with clinical results associated with a TH2 airway disease.

Long-term treatment with odanacatib maintains normal trabecular biomechanical properties in ovariectomized adult monkeys as demonstrated by micro-CT-based finite element analysis.

The cathepsin K inhibitor odanacatib (ODN) is a potent and reversible inhibitor of osteoclastic resorption activity. This drug is currently under development for the treatment of postmenopausal osteoporosis. Previously, we described data on the treatment efficacy of ODN in a preclinical estrogen-deficient model of an ovariectomized (OVX) rhesus monkey using HR-pQCT based finite element analysis (FEA) in vivo estimates of bone strength on the distal radius. To support the bone safety profile of ODN, we report ex vivo data on the apparent and hard tissue biomechanical properties of the trabecular bone of vertebrae of animals after 20 months of dosing in three treatment groups: Vehicle (VEH), ODN (2 mg/kg/day), and ALN (30 µg/kg/week). Biomechanical axial compression tests were performed on cylindrical trabecular samples cored out of the third lumbar vertebra of each animal at the end of the study. The biomechanical test results demonstrated that a normal (positive correlation between bone mineral density and bone strength) apparent material property relationship was maintained in the lumbar spine of ODN and ALN treated non-human primates (NHP). Trabecular bone hard tissue Young's modulus value was estimated using experimentally measured stiffness combined with FEA. The FEA and experimental results demonstrated that ODN treatment for 20 months maintained normal trabecular bone material hard tissue properties in the OVX-monkeys and was comparable to ALN.

Mathematical Model of Bone Remodeling Captures the Antiresorptive and Anabolic Actions of Various Therapies.

A better understanding of the molecular pathways regulating the bone remodeling process should help in the development of new antiresorptive regulators and anabolic regulators, that is, regulators of bone resorption and of bone formation. Understanding the mechanisms by which parathyroid hormone (PTH) influences bone formation and how it switches from anabolic to catabolic action is important for treating osteoporosis (Poole and Reeve in Curr Opin Pharmacol 5:612-617, 2005). In this paper we describe a mathematical model of bone remodeling that incorporates, extends, and integrates several models of particular aspects of this biochemical system (Cabal et al. in J Bone Miner Res 28(8):1830-1836, 2013; Lemaire et al. in J Theor Biol 229:293-309, 2004; Peterson and Riggs in Bone 46:49-63, 2010; Raposo et al. in J Clin Endocrinol Metab 87(9):4330-4340, 2002; Ross et al. in J Disc Cont Dyn Sys Series B 17(6):2185-2200, 2012). We plan to use this model as a bone homeostasis platform to develop anabolic and antiresorptive compounds. The model will allow us to test hypotheses about the dynamics of compounds and to test the potential benefits of combination therapies. At the core of the model is the idealized account of osteoclast and osteoblast signaling given by Lemaire et al. (J Theor Biol 229:293-309, 2004). We have relaxed some of their assumptions about the roles of osteoprotegerin, transforming growth factor [Formula: see text], and receptor activator of nuclear factor [Formula: see text]B ligand; we have devised more detailed models of the interactions of these species. We have incorporated a model of the effect of calcium sensing receptor antagonists on remodeling (Cabal et al. in J Bone Miner Res 28(8):1830-1836, 2013). We have also incorporated a basic model of the effects of vitamin D on calcium homeostasis. We have included a simple model of the mechanism proposed by Bellido et al. (2003), Ross et al. (J Disc Cont Dyn Sys Series B 17(6):2185-2200, 2012), of the influence of PTH on osteoblast apoptosis, a mechanism that accounts for the anabolic response to pulsatile PTH administration. Finally, we have devised a simple model of the administration and effects of bisphosphonates. The biomarkers in the model are procollagen type 1 amino-terminal propeptide and C-terminal telopeptide. Bone mineral density is the model's principal endpoint.

Preclinical Experimental and Mathematical Approaches for Assessing Effective Doses of Inhaled Drugs, Using Mometasone to Support Human Dose Predictions.

BACKGROUND: Understanding the relationship between dose, lung exposure, and drug efficacy continues to be a challenging aspect of inhaled drug development. An experimental inhalation platform was developed using mometasone furoate to link rodent lung exposure to its in vivo pharmacodynamic (PD) effects. METHODS: We assessed the effect of mometasone delivered directly to the lung in two different rodent PD models of lung inflammation. The data obtained were used to develop and evaluate a mathematical model to estimate drug dissolution, transport, distribution, and efficacy, following inhaled delivery in rodents and humans. RESULTS: Mometasone directly delivered to the lung, in both LPS and Alternaria alternata rat models, resulted in dose dependent inhibition of BALf cellular inflammation. The parameters for our mathematical model were calibrated to describe the observed lung and systemic exposure profiles of mometasone in humans and in animal models. We found that physicochemical properties, such as lung fluid solubility and lipophilicity, strongly influenced compound distribution and lung retention. CONCLUSIONS: Presently, we report on a novel and sophisticated mathematical model leading to improvements in a current inhaled drug development practices by providing a quantitative understanding of the relationship between PD effects and drug concentration in lungs.

Effects of odanacatib on the radius and tibia of postmenopausal women: improvements in bone geometry, microarchitecture, and estimated bone strength.

The cathepsin K inhibitor odanacatib (ODN), currently in phase 3 development for postmenopausal osteoporosis, has a novel mechanism of action that reduces bone resorption while maintaining bone formation. In phase 2 studies, odanacatib increased areal bone mineral density (aBMD) at the lumbar spine and total hip progressively over 5 years. To determine the effects of ODN on cortical and trabecular bone and estimate changes in bone strength, we conducted a randomized, double-blind, placebo-controlled trial, using both quantitative computed tomography (QCT) and high-resolution peripheral (HR-p)QCT. In previously published results, odanacatib was superior to placebo with respect to increases in trabecular volumetric BMD (vBMD) and estimated compressive strength at the spine, and integral and trabecular vBMD and estimated strength at the hip. Here, we report the results of HR-pQCT assessment. A total of 214 postmenopausal women (mean age 64.0 ± 6.8 years and baseline lumbar spine T-score -1.81 ± 0.83) were randomized to oral ODN 50 mg or placebo, weekly for 2 years. With ODN, significant increases from baseline in total vBMD occurred at the distal radius and tibia. Treatment differences from placebo were also significant (3.84% and 2.63% for radius and tibia, respectively). At both sites, significant differences from placebo were also found in trabecular vBMD, cortical vBMD, cortical thickness, cortical area, and strength (failure load) estimated using finite element analysis of HR-pQCT scans (treatment differences at radius and tibia = 2.64% and 2.66%). At the distal radius, odanacatib significantly improved trabecular thickness and bone volume/total volume (BV/TV) versus placebo. At a more proximal radial site, odanacatib attenuated the increase in cortical porosity found with placebo (treatment difference = -7.7%, p = 0.066). At the distal tibia, odanacatib significantly improved trabecular number, separation, and BV/TV versus placebo. Safety and tolerability were similar between treatment groups. In conclusion, odanacatib increased cortical and trabecular density, cortical thickness, aspects of trabecular microarchitecture, and estimated strength at the distal radius and distal tibia compared with placebo.

High-resolution peripheral quantitative computed tomography and finite element analysis of bone strength at the distal radius in ovariectomized adult rhesus monkey demonstrate efficacy of odanacatib and differentiation from alendronate.

Translational evaluation of disease progression and treatment response is critical to the development of therapies for osteoporosis. In this study, longitudinal in-vivo monitoring of odanacatib (ODN) treatment efficacy was compared to alendronate (ALN) in ovariectomized (OVX) non-human primates (NHPs) using high-resolution peripheral computed tomography (HR-pQCT). Treatment effects were evaluated using several determinants of bone strength, density and quality, including volumetric bone mineral density (vBMD), three-dimensional structure, finite element analysis (FEA) estimated peak force and biomechanical properties at the ultradistal (UD) radius at baseline, 3, 6, 9, 12, and 18 months of dosing in three treatment groups: vehicle (VEH), low ODN (2 mg/kg/day, L-ODN), and ALN (30 µg/kg/week). Biomechanical axial compression tests were performed at the end of the study. Bone strength estimates using FEA were validated by ex-vivo mechanical compression testing experiments. After 18months of dosing, L-ODN demonstrated significant increases from baseline in integral vBMD (13.5%), cortical thickness (24.4%), total bone volume fraction BV/TV (13.5%), FEA-estimated peak force (26.6%) and peak stress (17.1%), respectively. Increases from baseline for L-ODN at 18 months were significantly higher than that for ALN in DXA-based aBMD (7.6%), cortical thickness (22.9%), integral vBMD (12.2%), total BV/TV (10.1%), FEA peak force (17.7%) and FEA peak stress (11.5%), respectively. These results demonstrate a superior efficacy of ODN treatment compared to ALN at the UD radii in ovariectomized NHPs.

Effect of odanacatib on bone turnover markers, bone density and geometry of the spine and hip of ovariectomized monkeys: a head-to-head comparison with alendronate.

Odanacatib (ODN) is a selective and reversible Cathepsin K (CatK) inhibitor currently being developed as a once weekly treatment for osteoporosis. Here, effects of ODN compared to alendronate (ALN) on bone turnover, DXA-based areal bone mineral density (aBMD), QCT-based volumetric BMD (vBMD) and geometric parameters were studied in ovariectomized (OVX) rhesus monkeys. Treatment was initiated 10 days after ovariectomy and continued for 20 months. The study consisted of four groups: L-ODN (2 mg/kg, daily p.o.), H-ODN (8/4 mg/kg daily p.o.), ALN (15 µg/kg, twice weekly, s.c.), and VEH (vehicle, daily, p.o.). L-ODN and ALN doses were selected to approximate the clinical exposures of the ODN 50-mg and ALN 70-mg once-weekly, respectively. L-ODN and ALN effectively reduced bone resorption markers uNTx and sCTx compared to VEH. There was no additional efficacy with these markers achieved with H-ODN. Conversely, ODN displayed inversely dose-dependent reduction of bone formation markers, sP1NP and sBSAP, and L-ODN reduced formation to a lesser degree than ALN. At month 18 post-OVX, L-ODN showed robust increases in lumbar spine aBMD (11.4%, p<0.001), spine trabecular vBMD (13.7%, p<0.001), femoral neck (FN) integral (int) vBMD (9.0%, p<0.001) and sub-trochanteric proximal femur (SubTrPF) int vBMD, (6.4%, p<0.001) compared to baseline. L-ODN significantly increased FN cortical thickness (Ct.Th) and cortical bone mineral content (Ct.BMC) by 22.5% (p<0.001) and 21.8% (p<0.001), respectively, and SubTrPF Ct.Th and Ct.BMC by 10.9% (p<0.001) and 11.3% (p<0.001) respectively. Compared to ALN, L-ODN significantly increased FN Ct. BMC by 8.7% (p<0.05), and SubTrPF Ct.Th by 7.6% (p<0.05) and Ct.BMC by 6.2% (p<0.05). H-ODN showed no additional efficacy compared to L-ODN in OVX-monkeys in prevention mode. Taken together, the results from this study have demonstrated that administration of ODN at levels which approximate clinical exposure in OVX-monkeys had comparable efficacy to ALN in DXA-based aBMD and QCT-based vBMD. However, FN cortical mineral content clearly demonstrated superior efficacy of ODN versus ALN in this model of estrogen-deficient non-human primates.

A semimechanistic model of the time-course of release of PTH into plasma following administration of the calcilytic JTT-305/MK-5442 in humans.

JTT-305/MK-5442 is a calcium-sensing receptor (CaSR) allosteric antagonist being investigated for the treatment of osteoporosis. JTT-305/MK-5442 binds to CaSRs, thus preventing receptor activation by Ca(2+) . In the parathyroid gland, this results in the release of parathyroid hormone (PTH). Sharp spikes in PTH secretion followed by rapid returns to baseline are associated with bone formation, whereas sustained elevation in PTH is associated with bone resorption. We have developed a semimechanistic, nonpopulation model of the time-course relationship between JTT-305/MK-5442 and whole plasma PTH concentrations to describe both the secretion of PTH and the kinetics of its return to baseline levels. We obtained mean concentration data for JTT-305/MK-5442 and whole PTH from a multiple dose study in U.S. postmenopausal women at doses of 5, 10, 15, and 20 mg. We hypothesized that PTH is released from two separate sources: a reservoir that is released rapidly (within minutes) in response to reduction in Ca(2+) binding, and a second source released more slowly following hours of reduced Ca(2+) binding. We modeled the release rates of these reservoirs as maximum pharmacologic effect (Emax ) functions of JTT-305/MK-5442 concentration. Our model describes both the dose-dependence of PTH time of occurrence for maximum drug concentration (Tmax ) and maximum concentration of drug (Cmax ), and the extent and duration of the observed nonmonotonic return of PTH to baseline levels following JTT-305/MK-5442 administration.

Evaluation of high-resolution peripheral quantitative computed tomography, finite element analysis and biomechanical testing in a pre-clinical model of osteoporosis: a study with odanacatib treatment in the ovariectomized adult rhesus monkey.

This study aimed to validate finite element analysis (FEA) estimation of strength, identify high-resolution peripheral quantitative computed tomography (HR-pQCT) measures correlating with strength, and evaluate the precision of HR-pQCT measurements to longitudinally monitor effects of osteoporosis treatment in ovariectomized (OVX) non-human primates (NHPs). HR-pQCT images were acquired in three groups of NHPs: Intact (n=10), OVX-odanacatib treated (OVX-ODN 30 mg/kg, n=10) and OVX-vehicle treated (OVX-Veh, n=10) at the ultradistal (UD) and distal 1/3 radii and tibia at 12, 16 and 20 months. FEA estimates of bone strength using the Pistoia criterion were validated by ex-vivo mechanical compression (r(2)=0.95) of the UD radius. Single linear regressions of FEA-determined ultimate stress showed high correlation with HR-pQCT derived parameters: integral vBMD (r(2)=0.86), bone volume fraction (r(2)=0.84) and cortical thickness (r(2)=0.79). Precision of HR-pQCT measurements, obtained from an excised radius and tibia, showed low variation (CV=0.005%-5.6%) and helped identify possible sources of error. Comparison of OVX-Veh and Intact groups showed decreases in bone parameters demonstrating trends consistent with bone loss. Comparison of OVX-ODN and OVX-Veh groups showed a treatment effect with increases in bone parameters: integral vBMD (477±27 vs. 364±22 mgHA/cm(3)) and cortical thickness (Ct.Th) (0.90±0.07 vs. 0.64±0.04 mm) at the UD radius, Ct.Th (2.15±0.28 vs. 1.56±0.08 mm) at the distal 1/3 radius. Axial compression peak stress calculated and obtained experimentally showed the OVX-ODN group was 33% stronger than the OVX-Veh group. We conclude that HR-pQCT and FEA serve as robust techniques to longitudinally monitor bone parameters and strength in NHP's.