A Novel Method for Assessing Drug Degradation Product Safety Using Physiologically-Based Pharmacokinetic Models and Stochastic Risk Assessment

Exploring the complexities of drug formulation selection, storage, and shelf-life for exploration spaceflight

Medications have been a part of space travel dating back to the Apollo missions. Currently, medical kits aboard the International Space Station (ISS) contain medications and supplies to treat a variety of possible medical events. As we prepare for more distant exploration missions to Mars and beyond, risk management planning for astronaut healthcare should include the assembly of a medication formulary that is comprehensive enough to prevent or treat anticipated medical events, remains safe and chemically stable, and retains sufficient potency to last for the duration of the mission. Emerging innovation and technologies in pharmaceutical development, delivery, quality maintenance, and validation offer promise for addressing these challenges. The present editorial will summarize the current state of knowledge regarding innovative formulary optimization strategies, pharmaceutical stability assessment techniques, and storage and packaging solutions that could enhance drug safety and efficacy for future exploration spaceflight missions.

Formulating a heat- and shear-labile drug in an amorphous solid dispersion: Balancing drug degradation and crystallinity

We seek to further addresss the questions posed by Moseson et al. regarding whether any residual crystal level, size, or characteristic is acceptable in an amorphous solid dispersion (ASD) such that its stability, enhanced dissolution, and increased bioavailability are not compromised. To address this highly relevant question, we study an interesting heat- and shear-labile drug in development, LY3009120. To study the effects of residual crystallinity and degradation in ASDs, we prepared three compositionally identical formulations (57–1, 59–4, and 59–5) using the KinetiSol process under various processing conditions to obtain samples with various levels of crystallinity (2.3%, 0.9%, and 0.1%, respectively) and degradation products (0.74%, 1.97%, and 3.12%, respectively). Samples with less than 1% crystallinity were placed on stability, and we observed no measurable change in the drug's crystallinity, dissolution profile or purity in the 59–4 and 59–5 formulations over four months of storage under closed conditions at 25 °C and 60% humidity. For formulations 57–1, 59–4, and 59–5, bioavailability studies in rats reveal a 44-fold, 55-fold, and 62-fold increase in mean AUC, respectively, compared to the physical mixture. This suggests that the presence of some residual crystals after processing can be acceptable and will not change the properties of the ASD over time.

Physiologically Based Pharmacokinetic Models Predicting Renal and Hepatic Concentrations of Industrial Chemicals after Virtual Oral Doses in Rats

Recently developed high-throughput in vitro assays in combination with computational models could provide alternatives to animal testing. The purpose of the present study was to model the plasma, hepatic, and renal pharmacokinetics of approximately 150 structurally varied types of drugs, food components, and industrial chemicals after virtual external oral dosing in rats and to determine the relationship between the simulated internal concentrations in tissue/plasma and their lowest-observed-effect levels. The model parameters were based on rat plasma data from the literature and empirically determined pharmacokinetics measured after oral administrations to rats carried out to evaluate hepatotoxic or nephrotic potentials. To ensure that the analyzed substances exhibited a broad diversity of chemical structures, their structure-based location in the chemical space underwent projection onto a two-dimensional plane, as reported previously, using generative topographic mapping. A high-throughput in silico one-compartment model and a physiologically based pharmacokinetic (PBPK) model consisting of chemical receptor (gut), metabolizing (liver), central (main), and excreting (kidney) compartments were developed in parallel. For 159 disparate chemicals, the maximum plasma concentrations and the areas under the concentration-time curves obtained by one-compartment models and modified simple PBPK models were closely correlated. However, there were differences between the PBPK modeled and empirically obtained hepatic/renal concentrations and plasma maximal concentrations/areas under the concentration-time curves of the 159 chemicals. For a few compounds, the lowest-observed-effect levels were available for hepatotoxicity and nephrotoxicity in the Hazard Evaluation Support System Integrated Platform in Japan. The areas under the renal or hepatic concentration-time curves estimated using PBPK modeling were inversely associated with these lowest-observed-effect levels. Using PBPK forward dosimetry could provide the plasma/tissue concentrations of drugs and chemicals after oral dosing, thereby facilitating estimates of nephrotoxic or hepatotoxic potential as a part of the risk assessment.

Characterization of an Amphetamine Interference from Gabapentin in an LC-HRMS Method

An amphetamine interference was observed during the development of an liquid chromatography-high-resolution mass spectrometry (LC-HRMS) multi-class confirmation method for the determination of 47 drugs and metabolites in urine. The interference passed all qualitative criteria for amphetamine leading to potential false-positive results. Upon investigation, it was found that the amphetamine interference was correlated with the presence of high levels of gabapentin. Gabapentin is routinely detected in patient urine specimens at levels in excess of 1 mg/mL as it is widely prescribed at high doses and does not undergo significant metabolism. The source of the interference was identified as a gabapentin in-source fragment isomeric with protonated amphetamine. Here we describe the characterization of this interference and how its effect was mitigated in the LC-HRMS method.

Physiology-based toxicokinetic modelling in the frame of the European Human Biomonitoring Initiative

Given the opportunities provided by internal dosimetry modelling in the interpretation of human biomonitoring (HBM) data, the assessment of the links between exposure to chemicals and observed HBM data can be effectively supported by PBTK modelling. This paper gives a comprehensive review of available human PBTK models for compounds selected as a priority by the European Human Biomonitoring Initiative (HBM4EU). We highlight their advantages and deficiencies and suggest steps for advanced internal dose modelling. The review of the available PBTK models highlighted the conceptual differences between older models compared to the ones developed recently, reflecting commensurate differences in research questions. Due to the lack of coordinated strategies for deriving useful biomonitoring data for toxicokinetic properties, significant problems in model parameterisation still remain; these are further increased by the lack of human toxicokinetic data due to ethics issues. Finally, questions arise as well as to the extent they are really representative of interindividual variability. QSARs for toxicokinetic properties is a complementary approach for PBTK model parameterisation, especially for data poor chemicals. This approach could be expanded to model chemico-biological interactions such as intestinal absorption and renal clearance; this could serve the development of more complex generic PBTK models that could be applied to newly derived chemicals. Another gap identified is the framework for mixture interaction terms among compounds that could eventually interact in metabolism. From the review it was concluded that efforts should be shifted toward the development of generic multi-compartmental and multi-route models, supported by targeted biomonitoring coupled with parameterisation by both QSAR approach and experimental (in-vivo and in-vitro) data for newly developed and data poor compounds.

Plasma and Hepatic Concentrations of Chemicals after Virtual Oral Administrations Extrapolated Using Rat Plasma Data and Simple Physiologically Based Pharmacokinetic Models

Only a small fraction of chemicals possesses adequate in vivo toxicokinetic data for assessing potential hazards. The aim of the present study was to model the plasma and hepatic pharmacokinetics of more than 50 disparate types of chemicals and drugs after virtual oral administrations in rats. The models were based on reported pharmacokinetics determined after oral administration to rats. An inverse relationship was observed between no-observed-effect levels after oral administration and chemical absorbance rates evaluated for cell permeability ( r = -0.98, p < 0.001, n = 17). For a varied selection of more than 30 chemicals, the plasma concentration curves and the maximum concentrations obtained using a simple one-compartment model (recently recommended as a high-throughput toxicokinetic model) and a simple physiologically based pharmacokinetic (PBPK) model (consisting of chemical receptor, metabolizing, and central compartments) were highly consistent. The hepatic and plasma concentrations and the hepatic and plasma areas under the concentration-time curves of more than 50 chemicals were roughly correlated; however, differences were evident between the PBPK-modeled values in livers and empirically obtained values in plasma. Of the compounds selected for analysis, only seven had the lowest observed effect level (LOEL) values for hepatoxicity listed in the Hazard Evaluation Support System Integrated Platform in Japan. For these seven compounds, the LOEL values and the areas under the hepatic concentration-time curves estimated using PBPK modeling were inversely correlated ( r = -0.78, p < 0.05, n = 7). This study provides important information to help simulate the high hepatic levels of potent hepatotoxic compounds. Using suitable PBPK parameters, the present models could estimate the plasma/hepatic concentrations of chemicals and drugs after oral doses using both PBPK forward and reverse dosimetry, thereby indicating the potential value of this modeling approach in predicting hepatic toxicity as a part of risk assessments of chemicals absorbed in the human body.

Particle Size Distribution Equivalency as Novel Predictors for Bioequivalence

The use of particle size distribution (PSD) similarity metrics and the development and incorporation of drug release predictions based on PSD properties into PBPK models for various drug administration routes may provide a holistic approach for evaluating the effect of PSD differences on in vitro drug release and bioavailability of disperse systems. The objectives of this study were to provide a rational approach for evaluating the utility of in vitro PSD comparators for predicting bioequivalence for subcutaneously administered test and reference drug emulsions. Two types of in vitro comparators for test and reference emulsion products were evaluated: PSD characterization comparators (overlap metrics, median, and span ratios) and release profile comparators (f₂ and various fractional time ratios). A subcutaneous-input PBPK disposition model was developed to simulate blood concentration-time profiles of reference and test emulsion products and pharmacokinetic responses (e.g., AUC, C_max, and T_max) were used to determine bioequivalence. A pool of 10,440 pairs of test and reference products was simulated using Monte Carlo experiments. The PSD and release profile comparators were correlated to pass/fail bioequivalence metrics using logistical regression. Based on the use of single in vitro comparators, the f₂ method was the best predictor of bioequivalence prediction. The use of combinations of f₂ and PSD overlap comparators (e.g., OVL or PROB) improved bioequivalence prediction to about 90%. Simulation procedures used in this study demonstrated a process for developing reliable in vitro BE predictors.