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.

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.

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.

Systematic Evaluation of the Effect of Formulation Variables on In Vitro Performance of Mometasone Furoate Suspension-Metered Dose Inhalers.

The therapeutic benefits of metered dose inhalers (MDIs) in pulmonary disorders are mainly driven by aerosol performance, which depends on formulation variables (drug and excipients), device design, and patient interactions. The present study provides a comprehensive investigation to better understand the effect of formulation variables on mometasone furoate (MF) suspension-based MDI product performance. The effects of MF particle size (volume median diameter; X50) and excipient concentration (ethanol and oleic acid, cosolvent, and surfactant, respectively) on selected critical quality attributes (delivered dose (DD), fine particle dose of particles lesser than 5 µm (FPD < 5), ex-throat dose and median dissolution time (MDT)) were studied. Eight MF-MDI formulations (one per batch) were manufactured based on a reduced factorial design of experiment (DOE) approach, which included relevant formulation levels with varying X50 (1.1 and 2 µm), concentration of ethanol (0.45, 0.9, 1.8, and 3.6%w/w), and oleic acid (0.001 and 0.025%w/w). The in vitro evaluation of these MF-MDI formulations indicated the importance of drug particle's X50, oleic acid, and ethanol canister concentration as critical formulation variables governing the performance of MF suspension-based MDI products. The effect of these formulation variables on DD, FPD < 5, ex-throat dose, and MDT was subsequently utilized to develop empirical relationships linking formulation factors with effects on in vitro performance measures. The developed strategy could be useful for predicting MF-MDI product performance during MDI product development and manufacturing. The systematic DOE approach utilized in this study may provide insights into the understanding of the formulation variables governing the MF-MDI product performance.

Semi-mechanistic PK/PD model to assess pulmonary targeting of beclomethasone dipropionate and its active metabolite.

PURPOSE: The objective of this study was to describe the pulmonary targeting of beclomethasone dipropionate (BDP) and its active metabolite beclomethasone 17-monopropionate (BMP) in rats using a semi-mechanistic PK/PD model. METHODS: Rat plasma and tissue concentrations of BDP and BMP, and tissue receptor occupancies of BMP after systemic and pulmonary delivery of BDP and BMP were integrated in a newly developed semi-mechanistic PK/PD model. RESULTS: After IV administration of BDP, 95.4% of BDP was converted to BMP, while after pulmonary delivery of BDP, 46.6% of deposited BDP was absorbed as BMP. The developed semi-mechanistic PK model described plasma and tissue concentrations of BDP and BMP as well as receptor occupancies sufficiently well. The model incorporated dissolution, metabolic activation, and drug absorption processes to describe the local fate of BDP and BMP after systemic and pulmonary delivery. Dissolution rate constants of BDP and BMP were estimated to be 0.47/h and 2.01/h, respectively, and the permeabilities in central lung were estimated to be 15.0 and 2.9 × 106 cm/s for BDP and BMP, respectively. The EC50 of the binding of BMP to to the receptor was estimated to be 0.0017 ng/ml. Overall, receptor occupancies in the lung were more pronounced than those in the systemic circulation after pulmonary delivery of BDP or BMP. Simulations using the developed semi-mechanistic PK/PD model demonstrated that a slow dissolution rate and low permeability can improve pulmonary targeting. CONCLUSIONS: A semi-mechanistic model was developed to describe the fate of an inhaled glucocorticoid pro-drug and its active metabolite in lung and the systemic circulation, both after pulmonary and systemic administration , thereby facilitating the understanding of the complex interplay between drug, prodrug and pharmacodynamic properties for quantifying the degree pulmonary targeting.

Huntington's Disease Progression: A Population Modeling Approach to Characterization Using Clinical Rating Scales.

Development of effective therapeutics that slow Huntington's disease progression is a research priority that requires an understanding of natural disease progression. We applied a population-modeling approach to describe the progression of 2 routinely used rating scales - the total motor score and the total functional capacity score. Models were fitted to data from research participants aged ≥ 18 years with Huntington's disease stage I or II at study entry (total functional capacity score ≥ 7), from a controlled clinical trial (CARE-HD) and 2 observational studies (COHORT and Registry). A logistic model without shape factors was selected as the base model based on placebo data from CARE-HD and validated using data from the CARE-HD active-treatment arms. Albeit with a smaller progression rate constant than was found in CARE-HD, the proposed models provided reasonable predictions for both rating scales in the pooled data from COHORT and Registry and were considered suitable for use in clinical trial simulations. Results also showed that disease burden score (a product of age and expanded CAG length) is a significant covariate on both the progression rate constant and the baseline score in the total motor score model. These findings suggest that total motor score and total functional capacity progress fastest near their half-maximal score, implying that the efficiency of clinical trials evaluating disease-modifying therapeutics for Huntington's disease could be enhanced by enrolling patients with faster disease progression or evaluating treatment effect near their half-maximal score, provided that the evaluated therapy is expected to be efficacious at this disease stage.

Quantitative Assessment of Pulmonary Targeting of Inhaled Corticosteroids Using Ex Vivo Receptor Binding Studies.

The goal of locally acting inhaled corticosteroids is to achieve distinct pulmonary effects with reduced systemic side effects. The present work using an ex vivo receptor binding model in rats was interested in assessing pulmonary targeting for several commercially available corticosteroids by monitoring receptor occupancies in the lung and systemic organs (liver, kidney, spleen, and brain) after intravenous (IV) injection or intratracheal (IT) instillation of a dry powder administration at a dose of 100 µg/kg. Pulmonary targeting, defined as the difference in cumulative receptor occupancies (AUCE) between the lung and kidney after pulmonary delivery, differed across the investigated corticosteroids (ΔAUCE range, 33 ± 46 to 143 ± 52% *h) with the highest degree found for corticosteroids with high systemic clearance and pronounced lipophilicity (presumably allowing a long pulmonary residence time). Additionally, this study demonstrated differences in the receptor occupancies across systemic organs. Using kidney receptor occupancies as the comparator, liver receptor occupancies were reduced (ΔAUCE range: - 157 ± 43 to 178 ± 42% *h) after IV and IT administration for corticosteroids with high intrinsic clearance, while they were increased for corticosteroid prodrugs due to hepatic activation. Spleen receptor occupancies were increased after IT (ΔAUCE range: 33 ± 35 to 135 ± 28% *h), but not after IV administration. This was especially true for slowly dissolving drugs. Reduced brain uptake was also observed for ciclesonide (CIC) and des-ciclesonide (desCIC), two compounds previously not investigated. In summary, ex vivo receptor binding studies represent a powerful tool to assess the fate of ICSs.

Pharmacokinetic and pharmacodynamic modeling of gut hormone peptide YY(3-36) after pulmonary delivery.

Peptide YY(3-36) (PYY(3-36)) is an endogenous appetite suppressing peptide. The present research was to perform pharmacokinetic/pharmacodynamic (PK/PD) analysis for predicting the concentration- and response-time profiles of PYY(3-36) after systemic and pulmonary delivery in mice, with the goal of suggesting a potential pulmonary dosing regimen in humans. A PK/PD model was developed to describe PYY(3-36) plasma concentration - and relative food intake rate ratio (as % of control) - time profiles after intraperitoneal and subcutaneous administration, and inhalation in mice. The absorption of inhaled PYY(3-36) from the lungs of mice could only be described with a combined slow (absorption rate of 0.147 L/h) and fast (absorption rate of 104.4 L/h) absorption process, presumably related to absorption from the central and peripheral regions of the lungs. The estimates for IC50 and Imax were 6.8 ng/mL and 63.5%, respectively, based on inhibitory Emax model. The PK parameters, such as clearance (CL), volume of distribution at steady state (Vdss), and the absorption rates (ka), were then scaled to human's. The scaled human CL and Vdss for obese subjects were 24.8 L/h and 9.0 L, respectively. The model predicted human plasma PYY(3-36) concentrations agreed reasonably well with placebo-normalized plasma PYY(3-36) concentrations after short-term infusion and SC injection in literature. An inhalation dose of PYY(3-36) of about 100 µg was proposed for obese subjects based on simulations. This PK/PD analysis satisfactorily described PYY(3-36) concentration-time and relative food intake rate ratio- time profiles at all doses and routes. The developed model might facilitate the inhalation dose selection of PYY(3-36).

Effects of lipid formulations on clove extract spray dried powders: comparison of physicochemical properties, storage stability and in vitro intestinal permeation.

Clove is an aromatic plant spice with potent antioxidant and anti-inflammatory activity. Eugenol is the main compound which contributes to such medicinal and nutritional benefits. To date, the formulation of unstable, volatile and poorly water-soluble compounds remains a challenging task. Lipid formulations can be used to improve physicochemical and biopharmaceutical properties of poorly soluble compounds. The aim of this study is to investigate the effects of lipids, such as Gelucire and Compritol on physicochemical properties; stability and in vitro intestinal permeation of spray dried powdered formulations loaded with clove's bioactive compounds. Results showed that eugenol retention in spray-dried powders could be correlated with antioxidant activity and with mass recovery after spray drying. Adding Gelucire but not Compritol to clove extract formulations, improved solubility of spray dried powders. Stability test in high humidity environment (63.5% RH) suggested that formulations including both Gelucire and Compritol were significantly more stable compared to the formulation without any lipid at the two tested temperatures (25 °C and 40 °C). This suggests that lipid additions to clove (Syzygium aromaticum) extract formulations provide protective effects for the spray dried powders in high-humidity environments. In addition, results from in vitro intestinal permeation studies suggested that eugenol uptake, was not being hindered by transporters nor was the absorption being affected by lipid formulations.

Quantitative characterization of circadian rhythm of pulmonary function in asthmatic patients treated with inhaled corticosteroids.

The aim of this study was to characterize the circadian rhythm observed for forced expiratory volume in 1 s (FEV1) in patients with persistent asthma being treated with inhaled corticosteroids. The database included 3379 FEV1 measurements from 189 patients with mild to moderate asthma. A model using the sum of two Sine functions with periods of 12 and 24 h and a constant component of mean circadian rhythm adequately described the circadian rhythm in FEV1 measurements over time. The model adequateness was evaluated by various approaches including visual predictive check (VPC), prediction-corrected VPC, standardized VPC and normalized prediction distribution error. Covariates tested included age, body weight, height, body mass index, baseline FEV1, and sex. Age and height were found to have significant effects on the mean FEV1 level and no covariate was found to have an effect on the magnitude and timing of circadian rhythm. The model predicted that a minimum FEV1 occurred in the early morning and maximum FEV1 occurred in the early afternoon, with a population mean fluctuation of 170 mL, which is consistent with the finding that asthma symptoms usually exacerbate in the early morning for patients with persistent asthma. This developed model provides the first quantitative approach to describing FEV1 circadian rhythm with ICS background treatment and provided insight in designing future registration trials for asthma drug development.

Current Scientific and Regulatory Approaches for Development of Orally Inhaled and Nasal Drug Products: Overview of the IPAC-RS/University of Florida Orlando Inhalation Conference.

This article summarizes discussions at the March 2014 conference organized by the University of Florida (UF) and International Pharmaceutical Aerosol Consortium on Regulation and Science (IPAC-RS), entitled "Orlando Inhalation Conference: Approaches in International Regulation." The special focus of the conference was on global scientific and regulatory issues associated with the testing and demonstration of equivalence for the registration of orally inhaled drug products (OIDPs) in the United States, Europe, Brazil, China, and India. The scope included all types of OIDPs throughout their lifecycle, e.g., innovator/brand-name products, generics, modifications due to lifecycle management, device changes, etc. Details were presented for the U.S. "weight of evidence approach" for registration of generic products (which includes demonstration of in vitro and in vivo equivalence, as well as quantitative and qualitative sameness, and device similarity). The European "stepwise" approach was elucidated, and the thinking of regulatory agencies in the major emerging markets was clarified. The conference also highlighted a number of areas that would benefit from further research and discussion, especially around patient/device interface and human factor studies, statistical methods and criteria for demonstrating equivalence, the relative roles of in vivo and in vitro tests, and appropriate designs and metrics for in vivo studies of inhaled drugs.

A Systematic Analysis of the Sensitivity of Plasma Pharmacokinetics to Detect Differences in the Pulmonary Performance of Inhaled Fluticasone Propionate Products Using a Model-Based Simulation Approach.

The role of plasma pharmacokinetics (PK) for assessing bioequivalence at the target site, the lung, for orally inhaled drugs remains unclear. A validated semi-mechanistic model, considering the presence of mucociliary clearance in central lung regions, was expanded for quantifying the sensitivity of PK studies in detecting differences in the pulmonary performance (total lung deposition, central-to-peripheral lung deposition ratio, and pulmonary dissolution characteristics) between test (T) and reference (R) inhaled fluticasone propionate (FP) products. PK bioequivalence trials for inhaled FP were simulated based on this PK model for a varying number of subjects and T products. The statistical power to conclude bioequivalence when T and R products are identical was demonstrated to be 90% for approximately 50 subjects. Furthermore, the simulations demonstrated that PK metrics (area under the concentration time curve (AUC) and C max) are capable of detecting differences between T and R formulations of inhaled FP products when the products differ by more than 20%, 30%, and 25% for total lung deposition, central-to-peripheral lung deposition ratio, and pulmonary dissolution characteristics, respectively. These results were derived using a rather conservative risk assessment approach with an error rate of <10%. The simulations thus indicated that PK studies might be a viable alternative to clinical studies comparing pulmonary efficacy biomarkers for slowly dissolving inhaled drugs. PK trials for pulmonary efficacy equivalence testing should be complemented by in vitro studies to avoid false positive bioequivalence assessments that are theoretically possible for some specific scenarios. Moreover, a user-friendly web application for simulating such PK equivalence trials with inhaled FP is provided.

Application of the modified chi-square ratio statistic in a stepwise procedure for cascade impactor equivalence testing.

Equivalence testing of aerodynamic particle size distribution (APSD) through multi-stage cascade impactors (CIs) is important for establishing bioequivalence of orally inhaled drug products. Recent work demonstrated that the median of the modified chi-square ratio statistic (MmCSRS) is a promising metric for APSD equivalence testing of test (T) and reference (R) products as it can be applied to a reduced number of CI sites that are more relevant for lung deposition. This metric is also less sensitive to the increased variability often observed for low-deposition sites. A method to establish critical values for the MmCSRS is described here. This method considers the variability of the R product by employing a reference variance scaling approach that allows definition of critical values as a function of the observed variability of the R product. A stepwise CI equivalence test is proposed that integrates the MmCSRS as a method for comparing the relative shapes of CI profiles and incorporates statistical tests for assessing equivalence of single actuation content and impactor sized mass. This stepwise CI equivalence test was applied to 55 published CI profile scenarios, which were classified as equivalent or inequivalent by members of the Product Quality Research Institute working group (PQRI WG). The results of the stepwise CI equivalence test using a 25% difference in MmCSRS as an acceptance criterion provided the best matching with those of the PQRI WG as decisions of both methods agreed in 75% of the 55 CI profile scenarios.

Population pharmacokinetic modeling of the unbound levofloxacin concentrations in rat plasma and prostate tissue measured by microdialysis.

Levofloxacin is a broad-spectrum fluoroquinolone used in the treatment of both acute and chronic bacterial prostatitis. Currently, the treatment of bacterial prostatitis is still difficult, especially due to the poor distribution of many antimicrobials into the prostate, thus preventing the drug to reach effective interstitial concentrations at the infection site. Newer fluoroquinolones show a greater penetration into the prostate. In the present study, we compared the unbound levofloxacin prostate concentrations measured by microdialysis to those in plasma after a 7-mg/kg intravenous bolus dose to Wistar rats. Plasma and dialysate samples were analyzed using a validated high-pressure liquid chromatography-fluorescence method. Both noncompartmental analysis (NCA) and population-based compartmental modeling (NONMEM 6) were performed. Unbound prostate tissue concentrations represented 78% of unbound plasma levels over a period of 12 h by comparing the extent of exposure (unbound AUC0-∞) of 6.4 and 4.8 h·µg/ml in plasma and tissue, respectively. A three-compartment model with simultaneous passive diffusion and saturable distribution kinetics from the prostate to the central compartment gave the best results in terms of curve fitting, precision of parameter estimates, and model stability. The following parameter values were estimated by the population model: V1 (0.38 liter; where V1 represents the volume of the central compartment), CL (0.22 liter/h), k12 (2.27 h(-1)), k21 (1.44 h(-1)), k13 (0.69 h(-1)), Vmax (7.19 µg/h), kM (0.35 µg/ml), V3/fuprostate (0.05 liter; where fuprostate represents the fraction unbound in the prostate), and k31 (3.67 h(-1)). The interindividual variability values for V1, CL, Vmax, and kM were 21, 37, 42, and 76%, respectively. Our results suggest that levofloxacin is likely to be substrate for efflux transporters in the prostate.

Development of a placebo effect model combined with a dropout model for bipolar disorder.

The aim of this study was to develop a placebo model for bipolar disorder to help optimize clinical trial designs for studies targeting manic episodes in bipolar disorder. A bipolar disease database was built based on individual longitudinal data collected from over 3,000 patients in 11 clinical trials for 5 approved bipolar drugs. An empirical placebo effect model with an exponential decay process plus a linear progression process was developed to quantify the time course of the Young Mania Rating Scale total score based on only placebo data from the database. In order to describe the dropout pattern during the trials, a parametric survival model was developed and the Weibull distribution was identified to be the best distribution to describe the data. Based on the likelihood ratio test, it was found that patients with higher baseline score, slower disease improvement and more rapid disease progression tended to dropout earlier, and the trial features such as trial starting year and trial site were also significant covariates for dropout. A combination of the placebo effect model and the dropout model was applied to simulate new clinical trials through Monte-Carlo simulation. Both the placebo effect model and dropout model described the observed data reasonably well based on various diagnostic plots. The joint placebo response and dropout models can serve as a tool to simulate the most likely level of placebo response with the expected dropout pattern to help design a new clinical trial.

Bioequivalence of inhaled drugs: fundamentals, challenges and perspectives.

Interest in bioequivalence (BE) of inhaled drugs derives largely from the desire to offer generic substitutes to successful drug products. The complexity of aerosol dosage forms renders them difficult to mimic and raises questions regarding definitions of similarities and those properties that must be controlled to guarantee both the quality and the efficacy of the product. Despite a high level of enthusiasm to identify and control desirable properties there is no clear guidance, regulatory or scientific, for the variety of aerosol dosage forms, on practical measures of BE from which products can be developed. As more data on the pharmaceutical and clinical relevance of various techniques, as described in this review, become available, it is likely that a path to the demonstration of BE will become evident. In the meantime, debate on this topic will continue.

A sensitivity analysis of the modified chi-square ratio statistic for equivalence testing of aerodynamic particle size distribution.

Demonstration of equivalence in aerodynamic particle size distribution (APSD) is one key component for establishing bioequivalence of orally inhaled drug products. We previously proposed a modified version of the Chi-square ratio statistic (mCSRS) for APSD equivalence testing and demonstrated that the median of the distribution of the mCSRS (MmCSRS) is a robust metric when test (T) and reference (R) cascade impactor (CI) profiles are identical. Here, we systematically evaluate the behavior of the MmCSRS when T and R CI profiles differ from each other in their mean deposition and variability on a single and multiple sites. All CI profiles were generated by Monte-Carlo simulations based upon modified actual CI data. Twenty thousand sets of 30 T and 30 R CI profiles were simulated for each scenario, and the behavior of the MmCSRS was correlated to metrics that characterize the difference between T and R product in mean deposition and variability. The two key findings were, first, that the MmCSRS is more sensitive to difference between T and R CI profiles on high deposition sites, and second, that a cut-off value for APSD equivalence testing based on the MmCSRS needs to be scaled on the variability of the R product. The former is considered as beneficial for equivalence testing of CI profiles as it decreases the likelihood of failing identical CI profiles by chance, in part, due to increasing analytical variability associated with lower deposition sites. The latter is expected to be important for consistently being able to discriminate equivalent from inequivalent CI profiles.

A stability analysis of a modified version of the chi-square ratio statistic: implications for equivalence testing of aerodynamic particle size distribution.

Demonstration of equivalence in aerodynamic particle size distribution (APSD; e.g., by comparing cascade impactor (CI) profiles) constitutes one of key in vitro tests for supporting bioequivalence between test (T) and reference (R) orally inhaled drug products (OIDPs). A chi-square ratio statistic (CSRS) was previously proposed for equivalence testing of CI profiles. However, it was reported that the CSRS could not consistently discriminate between equivalent and inequivalent CI profiles. The objective of the overall project was to develop a robust and sensitive methodology for assessing equivalence of APSD profiles of T and R OIDPs. We propose here a modified version of the CSRS (mCSRS) and evaluated systematically its behavior when T and R CI profiles were identical. Different scenarios comprising CI profiles with different number of deposition sites and shapes were generated by Monte-Carlo simulation. For each scenario, the mCSRS was applied to 20,000 independent sets of 30 T and 30 R CI profiles that were identical. Different metrics (including mean and median) of the distribution of 900 mCSRSs (30 T × 30 R) were then evaluated for their suitability as a test statistic (i.e., independent of the number of sites and shape of the CI profile) for APSD equivalence testing. The median of the distribution of 900 mCSRSs (MmCSRS) was one regardless of the number of sites and shape of the CI profile. Hence, the MmCSRS is a robust metric for CI profile equivalence testing when T and R CI profiles are identical and potentially useful for APSD equivalence testing.

A pharmacokinetic simulation tool for inhaled corticosteroids.

The pharmacokinetic (PK) behavior of inhaled drugs is more complicated than that of other forms of administration. In particular, the effects of certain physiological (mucociliary clearance and differences in membrane properties in central and peripheral (C/P) areas of the lung), formulation (as it relates to drug deposition and particle dissolution rate), and patient-related factors (lung function; effects on C/P deposition ratio) affect the systemic PKs of inhaled drugs. The objectives of this project were (1) to describe a compartmental model that adequately describes the fate of inhaled corticosteroids (ICS) after administration while incorporating variability between and within subjects and (2) based upon the model, to provide a freely available tool for simulation of PK trials after ICS administration. This compartment model allows for mucociliary removal of undissolved particles from the lung, distinguishes between central and peripheral regions of the lung, and models drug entering the systemic circulation via the lung and the gastrointestinal tract. The PK simulation tool is provided as an extension package to the statistical software R ('ICSpkTS'). It allows simulation of PK trials for hypothetical ICS and of four commercially available ICS (budesonide, flunisolide, fluticasone propionate, and triamcinolone acetonide) in a parallel study design. Simulated PK data and parameters agreed well with literature data for all four ICS. The ICSpkTS package is especially suitable to explore the effect of changes in model parameters on PK behavior and can be easily adjusted for other inhaled drugs.