The Nonclinical Disposition and Pharmacokinetic/Pharmacodynamic Properties of N-Acetylgalactosamine-Conjugated Small Interfering RNA Are Highly Predictable and Build Confidence in Translation to Human.

Conjugation of oligonucleotide therapeutics, including small interfering RNAs (siRNAs) or antisense oligonucleotides, to N-acetylgalactosamine (GalNAc) ligands has become the primary strategy for hepatocyte-targeted delivery, and with the recent approvals of GIVLAARI (givosiran) for the treatment of acute hepatic porphyria, OXLUMO (lumasiran) for the treatment of primary hyperoxaluria, and Leqvio (inclisiran) for the treatment of hypercholesterolemia, the technology has been well validated clinically. Although much knowledge has been gained over decades of development, there is a paucity of published literature on the drug metabolism and pharmacokinetic properties of GalNAc-siRNA. With this in mind, the goals of this minireview are to provide an aggregate analysis of these nonclinical absorption, distribution, metabolism, and excretion (ADME) data to build confidence on the translation of these properties to human. Upon subcutaneous administration, GalNAc-conjugated siRNAs are quickly distributed to the liver, resulting in plasma pharmacokinetic (PK) properties that reflect rapid elimination through asialoglycoprotein receptor-mediated uptake from circulation into hepatocytes. These studies confirm that liver PK, including half-life and, most importantly, siRNA levels in RNA-induced silencing complex in hepatocytes, are better predictors of pharmacodynamics (PD) than plasma PK. Several in vitro and in vivo nonclinical studies were conducted to characterize the ADME properties of GalNAc-conjugated siRNAs. These studies demonstrate that the PK/PD and ADME properties of GalNAc-conjugated siRNAs are highly conserved across species, are largely predictable, and can be accurately scaled to human, allowing us to identify efficacious and safe clinical dosing regimens in the absence of human liver PK profiles. SIGNIFICANCE STATEMENT: Several nonclinical ADME studies have been conducted in order to provide a comprehensive overview of the disposition and elimination of GalNAc-conjugated siRNAs and the pharmacokinetic/pharmacodynamic translation between species. These studies demonstrate that the ADME properties of GalNAc-conjugated siRNAs are well correlated and predictable across species, building confidence in the ability to extrapolate to human.

Reconstructing exposures from biomarkers using exposure-pharmacokinetic modeling--A case study with carbaryl.

Sources of uncertainty involved in exposure reconstruction for short half-life chemicals were characterized using computational models that link external exposures to biomarkers. Using carbaryl as an example, an exposure model, the Cumulative and Aggregate Risk Evaluation System (CARES), was used to generate time-concentration profiles for 500 virtual individuals exposed to carbaryl. These exposure profiles were used as inputs into a physiologically based pharmacokinetic (PBPK) model to predict urinary biomarker concentrations. These matching dietary intake levels and biomarker concentrations were used to (1) compare three reverse dosimetry approaches based on their ability to predict the central tendency of the intake dose distribution; and (2) identify parameters necessary for a more accurate exposure reconstruction. This study illustrates the trade-offs between using non-iterative reverse dosimetry methods that are fast, less precise and iterative methods that are slow, more precise. This study also intimates the necessity of including urine flow rate and elapsed time between last dose and urine sampling as part of the biomarker sampling collection for better interpretation of urinary biomarker data of short biological half-life chemicals. Resolution of these critical data gaps can allow exposure reconstruction methods to better predict population-level intake doses from large biomonitoring studies.

Development of a Practical Nomogram for Personalized Anemia Management in Patients Treated with Ataxia Telangiectasia and Rad3-related Inhibitor Camonsertib.

PURPOSE: Camonsertib is a highly selective and potent inhibitor of ataxia telangiectasia and Rad3-related (ATR) kinase. Dose-dependent anemia is a class-related on-target adverse event often requiring dose modifications. Individual patient risk factors for the development of significant anemia complicate the selection of a "one-size-fits-all" ATR inhibitor (ATRi) dose and schedule, possibly leading to suboptimal therapeutic doses in patients at low risk of anemia. We evaluated whether early predictors of anemia could be identified to ultimately inform a personalized dose-modification approach. PATIENTS AND METHODS: On the basis of preclinical observations and a mechanistic understanding of ATRi-related anemia, we identified several potential factors to explore in a multivariable linear regression modeling tool for predicting hemoglobin level ahead of day 22 (cycle 2) of treatment. RESULTS: In patients treated with camonsertib monotherapy (NCT04497116), we observed that hemoglobin decline is consistently preceded by reticulocytopenia, and dose- and exposure-dependent decreases in monocytes. We developed a nomogram incorporating baseline and day 8 hemoglobin and reticulocyte values that predicted the day 22 hemoglobin values of patients with clinically valuable concordance (within 7.5% of observations) 80% of the time in a cross-validation performance test of data from 60 patients. CONCLUSIONS: The prediction of future hemoglobin decrease, after a week of treatment, may enable a personalized, early dose modification to prevent development of clinically significant anemia and resulting unscheduled dose holds or transfusions.

Current Approaches for Predicting Human PK for Small Molecule Development Candidates: Findings from the IQ Human PK Prediction Working Group Survey.

Accurate prediction of human clearance (CL) and volume of distribution at steady state (Vd,ss) for small molecule drug candidates is an essential component of assessing likely efficacious dose and clinical safety margins. In 2021, the IQ Consortium Human PK Prediction Working Group undertook a survey of IQ member companies to understand the current PK prediction methods being used to estimate these parameters across the pharmaceutical industry. The survey revealed a heterogeneity in approaches being used across the industry (e.g., the use of allometric approaches, differing incorporation of binding terms, and inconsistent use of empirical correction factors for in vitro-in vivo extrapolation, IVIVE), which could lead to different PK predictions with the same input data. Member companies expressed an interest in improving human PK predictions by identifying the most appropriate compound-class specific methods, as determined by physiochemical properties and knowledge of CL pathways. Furthermore, there was consensus that increased understanding of the uncertainty inherent to the compound class-dependent prediction would be invaluable in aiding communication of human PK and dose uncertainty at the time of candidate nomination for development. The human PK Prediction Working Group is utilizing these survey findings to help interrogate clinical IV datasets from across the IQ consortium member companies to understand PK prediction accuracy and uncertainty from preclinical datasets.

Expanding RNAi therapeutics to extrahepatic tissues with lipophilic conjugates.

Therapeutics based on short interfering RNAs (siRNAs) delivered to hepatocytes have been approved, but new delivery solutions are needed to target additional organs. Here we show that conjugation of 2'-O-hexadecyl (C16) to siRNAs enables safe, potent and durable silencing in the central nervous system (CNS), eye and lung in rodents and non-human primates with broad cell type specificity. We show that intrathecally or intracerebroventricularly delivered C16-siRNAs were active across CNS regions and cell types, with sustained RNA interference (RNAi) activity for at least 3 months. Similarly, intravitreal administration to the eye or intranasal administration to the lung resulted in a potent and durable knockdown. The preclinical efficacy of an siRNA targeting the amyloid precursor protein was evaluated through intracerebroventricular dosing in a mouse model of Alzheimer's disease, resulting in amelioration of physiological and behavioral deficits. Altogether, C16 conjugation of siRNAs has the potential for safe therapeutic silencing of target genes outside the liver with infrequent dosing.

On the potential of in vitro organ-chip models to define temporal pharmacokinetic-pharmacodynamic relationships.

Functional human-on-a-chip systems hold great promise to enable quantitative translation to in vivo outcomes. Here, we explored this concept using a pumpless heart only and heart:liver system to evaluate the temporal pharmacokinetic/pharmacodynamic (PKPD) relationship for terfenadine. There was a time dependent drug-induced increase in field potential duration in the cardiac compartment in response to terfenadine and that response was modulated using a metabolically competent liver module that converted terfenadine to fexofenadine. Using this data, a mathematical model was developed to predict the effect of terfenadine in preclinical species. Developing confidence that microphysiological models could have a transformative effect on drug discovery, we also tested a previously discovered proprietary AstraZeneca small molecule and correctly determined the cardiotoxic response to its metabolite in the heart:liver system. Overall our findings serve as a guiding principle to future investigations of temporal concentration response relationships in these innovative in vitro models, especially, if validated across multiple time frames, with additional pharmacological mechanisms and molecules representing a broad chemical diversity.

Associations of Perfluoroalkyl Substances (PFAS) with Lower Birth Weight: An Evaluation of Potential Confounding by Glomerular Filtration Rate Using a Physiologically Based Pharmacokinetic Model (PBPK).

BACKGROUND: Prenatal exposure to perfluoroalkyl substances (PFAS) has been associated with lower birth weight in epidemiologic studies. This association could be attributable to glomerular filtration rate (GFR), which is related to PFAS concentration and birth weight. OBJECTIVES: We used a physiologically based pharmacokinetic (PBPK) model of pregnancy to assess how much of the PFAS-birth weight association observed in epidemiologic studies might be attributable to GFR. METHODS: We modified a PBPK model to reflect the association of GFR with birth weight (estimated from three studies of GFR and birth weight) and used it to simulate PFAS concentrations in maternal and cord plasma. The model was run 250,000 times, with variation in parameters, to simulate a population. Simulated data were analyzed to evaluate the association between PFAS levels and birth weight due to GFR. We compared simulated estimates with those from a meta-analysis of epidemiologic data. RESULTS: The reduction in birth weight for each 1-ng/mL increase in simulated cord plasma for perfluorooctane sulfonate (PFOS) was 2.72 g (95% CI: -3.40, -2.04), and for perfluorooctanoic acid (PFOA) was 7.13 g (95% CI: -8.46, -5.80); results based on maternal plasma at term were similar. Results were sensitive to variations in PFAS level distributions and the strength of the GFR-birth weight association. In comparison, our meta-analysis of epidemiologic studies suggested that each 1-ng/mL increase in prenatal PFOS and PFOA levels was associated with 5.00 g (95% CI: -21.66, -7.78) and 14.72 g (95% CI: -8.92, -1.09) reductions in birth weight, respectively. CONCLUSION: Results of our simulations suggest that a substantial proportion of the association between prenatal PFAS and birth weight may be attributable to confounding by GFR and that confounding by GFR may be more important in studies with sample collection later in pregnancy.

Lifetime PCB 153 bioaccumulation and pharmacokinetics in pilot whales: Bayesian population PBPK modeling and Markov chain Monte Carlo simulations.

Physiologically based pharmacokinetic (PBPK) models for wild animal populations such as marine mammals typically have a high degree of model uncertainty and variability due to the scarcity of information and the embryonic nature of this field. Parameters values used in marine mammals models are usually taken from other mammalian species (e.g. rats or mice) and might not be entirely suitable to properly explain the kinetics of pollutants in marine mammals. Therefore, several parameters for a PBPK model for the bioaccumulation and pharmacokinetics of PCB 153 in long-finned pilot whales were estimated in the present study using the Bayesian approach executed with Markov chain Monte Carlo (MCMC) simulations. This method uses 'prior' information of the parameters, either from the literature or from previous model runs. The advantage is that this method uses such 'prior' parameters to calculate probability distributions to determine 'posterior' values that best explain the field observations. Those field observations or datasets were PCB 153 concentrations in blubber of long-finned pilot whales from Sandy Cape and Stanley, Tasmania, Australia. The model predictions showed an overall decrease in PCB 153 levels in blubber over the lifetime of the pilot whales. All parameters from the Sandy Cape model were updated using the Stanley dataset, except for the concentration of PCB 153 in the milk. The model presented here is a promising and preliminary start to PBPK modeling in long-finned pilot whales that would provide a basis for non-invasive studies in these protected marine mammals.

Is the relationship between prenatal exposure to PCB-153 and decreased birth weight attributable to pharmacokinetics?

BACKGROUND: A recent meta-analysis based on data from > 7,000 pregnancies reported an association between prenatal polychlorinated biphenyl (PCB)-153 exposure and reduced birth weight. Gestational weight gain, which is associated negatively with PCB levels in maternal and cord blood, and positively with birth weight, could substantially confound this association. OBJECTIVE: We sought to estimate the influence of gestational weight gain on the association between PCB-153 exposure and birth weight using a pharmacokinetic model. METHODS: We modified a recently published pharmacokinetic model and ran Monte Carlo simulations accounting for variability in physiologic parameters and their correlations. We evaluated the pharmacokinetic model by comparing simulated plasma PCB-153 levels during pregnancy to serial measurements in 10 pregnant women from another study population. We estimated the association between simulated plasma PCB-153 levels and birth weight using linear regression models. RESULTS: The plasma PCB-153 level profiles generated with the pharmacokinetic model were comparable to measured levels in 10 pregnant women. We estimated a 118-g decrease in birth weight (95% CI: -129, -106 g) for each 1-µg/L increase in simulated cord plasma PCB-153, compared with the 150-g decrease estimated based on the previous meta-analysis. The estimated decrease in birth weight was reduced to -6 g (95% CI: -18, 6 g) when adjusted for simulated gestational weight gain. CONCLUSION: Our findings suggest that associations previously noted between PCB levels and birth weight may be attributable to confounding by maternal weight gain during pregnancy.

Using population physiologically based pharmacokinetic modeling to determine optimal sampling times and to interpret biological exposure markers: The example of occupational exposure to styrene.

BACKGROUND: Biomonitoring of chemicals in the workplace provides an integrated characterization of exposure that accounts for uptake through multiple pathways and physiological parameters influencing the toxicokinetics. OBJECTIVES: We used the case of styrene to (i) determine the best times to sample venous blood and end-exhaled air, (ii) characterize the inter-individual variability in biological levels following occupational exposure and (iii) propose biological limit values using a population physiologically based pharmacokinetic (PBPK) model. METHODS: We performed Monte Carlo simulations with various physiological, exposure and workload scenarios. Optimal sampling times were identified through regression analyses between levels in biological samples and 24-h area under the arterial blood concentration vs. time curve. We characterized the variability in levels of styrene in biological samples for exposures to a time weighted average (TWA) of 20ppm. RESULTS: Simulations suggest that the best times to sample venous blood are at the end of shift in poorly ventilated workplaces and 15min after the shift in highly ventilated workplaces. Exhaled air samples are most informative 15min after the shift. For a light workload, simulated styrene levels have a median (5th-95th percentiles) of 0.4mg/l (0.2-0.6) in venous blood at the end of shift and 0.5ppm (0.3-0.8) in exhaled air 15min after the end of shift. CONCLUSION: This study supports the current BEI(®) of the ACGIH of 0.2mg/l of styrene in venous blood at the end of shift and indicates a biological limit value of 0.3ppm in end-exhaled air 15min after the end of shift.

A case study addressing the reliability of polychlorinated biphenyl levels measured at the time of breast cancer diagnosis in representing early-life exposure.

BACKGROUND: To date, breast cancer epidemiologic studies have relied on blood or tissue specimens sampled at the time of diagnosis or a few years prior to assess lifetime exposure to polychlorinated biphenyls (PCB). In this study, we evaluated whether such PCB measurements are indicative of early-life levels by reconstructing lifetime toxicokinetic profiles for women included in the CECILE case-control study, using a physiologically based pharmacokinetic (PBPK) model. METHODS: We simulated lifetime toxicokinetic profiles of PCB-153 for 2,134 French women by incorporating information on body weight history, height, pregnancies, and breast-feeding in the PBPK model. Oral dose was calculated by considering measured blood PCB-153 and the temporal trend of environmental contamination. Area under the concentration versus time curve (AUC) for each decade of life and maximum blood concentration (C(max)) were compiled and compared with measured levels, using Pearson partial correlation analyses adjusting for age at diagnosis. RESULTS: When considering all individuals, simulated AUCs correlated with measured PCBs, with coefficients ranging from 0.735 to 0.981. The weakest correlations were obtained with AUCs for the first decades of life. Stratified analyses suggested that breast-feeding reduces the reliability of late-life blood levels in representing lifetime exposure. CONCLUSION: Results of this study suggest that PCB levels measured at the time of diagnosis do not fully represent early-life exposures. IMPACT: PBPK-derived estimates of early-life levels circumvent the limitations of current approaches in assessing PCB lifetime exposure and may be used to address hypothesized windows of breast vulnerability (e.g., puberty) in this population.