A Bayesian Method for Population-wide Cardiotoxicity Hazard and Risk Characterization Using an In Vitro Human Model

Incorporating new approach methods (NAMs) data in dose-response assessments: The future is now!

Regulatory dose-response assessments traditionally rely on in vivo data and default assumptions. New Approach Methods (NAMs) present considerable opportunities to both augment traditional dose-response assessments and accelerate the evaluation of new/data-poor chemicals. This review aimed to determine the potential utilization of NAMs through a unified conceptual framework that compartmentalizes derivation of toxicity values into five sequential Key Dose-response Modules (KDMs): (1) point-of-departure (POD) determination, (2) test system-to-human (e.g. inter-species) toxicokinetics and (3) toxicodynamics, (4) human population (intra-species) variability in toxicodynamics, and (5) toxicokinetics. After using several "traditional" dose-response assessments to illustrate this framework, a review is presented where existing NAMs, including in silico, in vitro, and in vivo approaches, might be applied across KDMs. Further, the false dichotomy between "traditional" and NAMs-derived data sources is broken down by organizing dose-response assessments into a matrix where each KDM has Tiers of increasing precision and confidence: Tier 0: Default/generic values, Tier 1: Computational predictions, Tier 2: Surrogate measurements, and Tier 3: Direct measurements. These findings demonstrated that although many publications promote the use of NAMs in KDMs (1) for POD determination and (5) for human population toxicokinetics, the proposed matrix of KDMs and Tiers reveals additional immediate opportunities for NAMs to be integrated across other KDMs. Further, critical needs were identified for developing NAMs to improve in vitro dosimetry and quantify test system and human population toxicodynamics. Overall, broadening the integration of NAMs across the steps of dose-response assessment promises to yield higher throughput, less animal-dependent, and more science-based toxicity values for protecting human health.

Characterizing PFAS hazards and risks: a human population-based in vitro cardiotoxicity assessment strategy

Per- and poly-fluoroalkyl substances (PFAS) are emerging contaminants of concern because of their wide use, persistence, and potential to be hazardous to both humans and the environment. Several PFAS have been designated as substances of concern; however, most PFAS in commerce lack toxicology and exposure data to evaluate their potential hazards and risks. Cardiotoxicity has been identified as a likely human health concern, and cell-based assays are the most sensible approach for screening and prioritization of PFAS. Human-induced pluripotent stem cell (iPSC)-derived cardiomyocytes are a widely used method to test for cardiotoxicity, and recent studies showed that many PFAS affect these cells. Because iPSC-derived cardiomyocytes are available from different donors, they also can be used to quantify human variability in responses to PFAS. The primary objective of this study was to characterize potential human cardiotoxic hazard, risk, and inter-individual variability in responses to PFAS. A total of 56 PFAS from different subclasses were tested in concentration-response using human iPSC-derived cardiomyocytes from 16 donors without known heart disease. Kinetic calcium flux and high-content imaging were used to evaluate biologically-relevant phenotypes such as beat frequency, repolarization, and cytotoxicity. Of the tested PFAS, 46 showed concentration-response effects in at least one phenotype and donor; however, a wide range of sensitivities were observed across donors. Inter-individual variability in the effects could be quantified for 19 PFAS, and risk characterization could be performed for 20 PFAS based on available exposure information. For most tested PFAS, toxicodynamic variability was within a factor of 10 and the margins of exposure were above 100. This study identified PFAS that may pose cardiotoxicity risk and have high inter-individual variability. It also demonstrated the feasibility of using a population-based human in vitro method to quantify population variability and identify cardiotoxicity risks of emerging contaminants.

In vitro to in vivo extrapolation from 3D hiPSC-derived cardiac microtissues and physiologically based pharmacokinetic modeling to inform next-generation arrhythmia risk assessment

Proarrhythmic cardiotoxicity remains a substantial barrier to drug development as well as a major global health challenge. In vitro human pluripotent stem cell-based new approach methodologies have been increasingly proposed and employed as alternatives to existing in vitro and in vivo models that do not accurately recapitulate human cardiac electrophysiology or cardiotoxicity risk. In this study, we expanded the capacity of our previously established 3D human cardiac microtissue model to perform quantitative risk assessment by combining it with a physiologically based pharmacokinetic model, allowing a direct comparison of potentially harmful concentrations predicted in vitro to in vivo therapeutic levels. This approach enabled the measurement of concentration responses and margins of exposure for 2 physiologically relevant metrics of proarrhythmic risk (i.e. action potential duration and triangulation assessed by optical mapping) across concentrations spanning 3 orders of magnitude. The combination of both metrics enabled accurate proarrhythmic risk assessment of 4 compounds with a range of known proarrhythmic risk profiles (i.e. quinidine, cisapride, ranolazine, and verapamil) and demonstrated close agreement with their known clinical effects. Action potential triangulation was found to be a more sensitive metric for predicting proarrhythmic risk associated with the primary mechanism of concern for pharmaceutical-induced fatal ventricular arrhythmias, delayed cardiac repolarization due to inhibition of the rapid delayed rectifier potassium channel, or hERG channel. This study advances human-induced pluripotent stem cell-based 3D cardiac tissue models as new approach methodologies that enable in vitro proarrhythmic risk assessment with high precision of quantitative metrics for understanding clinically relevant cardiotoxicity.

Informing Hazard Identification and Risk Characterization of Environmental Chemicals by Combining Transcriptomic and Functional Data from Human-Induced Pluripotent Stem-Cell-Derived Cardiomyocytes

Environmental chemicals may contribute to the global burden of cardiovascular disease, but experimental data are lacking to determine which substances pose the greatest risk. Human-induced pluripotent stem cell (iPSC)-derived cardiomyocytes are a high-throughput cardiotoxicity model that is widely used to test drugs and chemicals; however, most studies focus on exploring electro-physiological readouts. Gene expression data may provide additional molecular insights to be used for both mechanistic interpretation and dose-response analyses. Therefore, we hypothesized that both transcriptomic and functional data in human iPSC-derived cardiomyocytes may be used as a comprehensive screening tool to identify potential cardiotoxicity hazards and risks of the chemicals. To test this hypothesis, we performed concentration-response analysis of 464 chemicals from 12 classes, including both pharmaceuticals and nonpharmaceutical substances. Functional effects (beat frequency, QT prolongation, and asystole), cytotoxicity, and whole transcriptome response were evaluated. Points of departure were derived from phenotypic and transcriptomic data, and risk characterization was performed. Overall, 244 (53%) substances were active in at least one phenotype; as expected, pharmaceuticals with known cardiac liabilities were the most active. Positive chronotropy was the functional phenotype activated by the largest number of tested chemicals. No chemical class was particularly prone to pose a potential hazard to cardiomyocytes; a varying proportion (10-44%) of substances in each class had effects on cardiomyocytes. Transcriptomic data showed that 69 (15%) substances elicited significant gene expression changes; most perturbed pathways were highly relevant to known key characteristics of human cardiotoxicants. The bioactivity-to-exposure ratios showed that phenotypic- and transcriptomic-based POD led to similar results for risk characterization. Overall, our findings demonstrate how the integrative use of in vitro transcriptomic and phenotypic data from iPSC-derived cardiomyocytes not only offers a complementary approach for hazard and risk prioritization, but also enables mechanistic interpretation of the in vitro test results to increase confidence in decision-making.

Reducing uncertainty in dose-response assessments by incorporating Bayesian benchmark dose modeling and in vitro data on population variability

There are two primary sources of uncertainty in the interpretability of toxicity values, like the reference dose (RfD): estimates of the point of departure (POD) and the absence of chemical-specific human variability data. We hypothesize two solutions-employing Bayesian benchmark dose (BBMD) modeling to refine POD determination and combining high-throughput toxicokinetic modeling with population-based toxicodynamic in vitro data to characterize chemical-specific variability. These hypotheses were tested by deriving refined probabilistic estimates for human doses corresponding to a specific effect size (M) in the Ith population percentile (HDM I) across 19 Superfund priority chemicals. HDM I values were further converted to biomonitoring equivalents in blood and urine for benchmarking against human data. Compared to deterministic default-based RfDs, HDM I values were generally more protective, particularly influenced by chemical-specific data on interindividual variability. Incorporating chemical-specific in vitro data improved precision in probabilistic RfDs, with a median 1.4-fold reduction in uncertainty variance. Comparison with US Environmental Protection Agency's Exposure Forecasting exposure predictions and biomonitoring data from the National Health and Nutrition Examination Survey identified chemicals with margins of exposure nearing or below one. Overall, to mitigate uncertainty in regulatory toxicity values and guide chemical risk management, BBMD modeling and chemical-specific population-based human in vitro data are essential.

Evaluating scientific confidence in the concordance of in vitro and in vivo protective points of departure

To fulfil the promise of reducing reliance on mammalian in vivo laboratory animal studies, new approach methods (NAMs) need to provide a confident basis for regulatory decision-making. However, previous attempts to develop in vitro NAMs-based points of departure (PODs) have yielded mixed results, with PODs from U.S. EPA's ToxCast, for instance, appearing more conservative (protective) but poorly correlated with traditional in vivo studies. Here, we aimed to address this discordance by reducing the heterogeneity of in vivo PODs, accounting for species differences, and enhancing the biological relevance of in vitro PODs. However, we only found improved in vitro-to-in vivo concordance when combining the use of Bayesian model averaging-based benchmark dose modeling for in vivo PODs, allometric scaling for interspecies adjustments, and human-relevant in vitro assays with multiple induced pluripotent stem cell-derived models. Moreover, the available sample size was only 15 chemicals, and the resulting level of concordance was only fair, with correlation coefficients <0.5 and prediction intervals spanning several orders of magnitude. Overall, while this study suggests several ways to enhance concordance and thereby increase scientific confidence in vitro NAMs-based PODs, it also highlights challenges in their predictive accuracy and precision for use in regulatory decision making.

Mapping Interindividual Variability of Toxicodynamics Using High-Throughput Transcriptomics and Primary Human Hepatocytes from Fifty Donors

BACKGROUND

Understanding the variability across the human population with respect to toxicodynamic responses after exposure to chemicals, such as environmental toxicants or drugs, is essential to define safety factors for risk assessment to protect the entire population. Activation of cellular stress response pathways are early adverse outcome pathway (AOP) key events of chemical-induced toxicity and would elucidate the estimation of population variability of toxicodynamic responses.

OBJECTIVES

We aimed to map the variability in cellular stress response activation in a large panel of primary human hepatocyte (PHH) donors to aid in the quantification of toxicodynamic interindividual variability to derive safety uncertainty factors.

METHODS

High-throughput transcriptomics of over 8,000 samples in total was performed covering a panel of 50 individual PHH donors upon 8 to 24 h exposure to broad concentration ranges of four different toxicological relevant stimuli: tunicamycin for the unfolded protein response (UPR), diethyl maleate for the oxidative stress response (OSR), cisplatin for the DNA damage response (DDR), and tumor necrosis factor alpha (TNF α ) for NF- κ B signaling. Using a population mixed-effect framework, the distribution of benchmark concentrations (BMCs) and maximum fold change were modeled to evaluate the influence of PHH donor panel size on the correct estimation of interindividual variability for the various stimuli.

RESULTS

Transcriptome mapping allowed the investigation of the interindividual variability in concentration-dependent stress response activation, where the average of BMCs had a maximum difference of 864-, 13-, 13-, and 259-fold between different PHHs for UPR, OSR, DDR, and NF- κ B signaling-related genes, respectively. Population modeling revealed that small PHH panel sizes systematically underestimated the variance and gave low probabilities in estimating the correct human population variance. Estimated toxicodynamic variability factors of stress response activation in PHHs based on this dataset ranged between 1.6 and 6.3.

DISCUSSION

Overall, by combining high-throughput transcriptomics and population modeling, improved understanding of interindividual variability in chemical-induced activation of toxicity relevant stress pathways across the human population using a large panel of plated cryopreserved PHHs was established, thereby contributing toward increasing the confidence of in vitro-based prediction of adverse responses, in particular hepatotoxicity. https://doi.org/10.1289/EHP11891.

Hazard and risk characterization of 56 structurally diverse PFAS using a targeted battery of broad coverage assays using six human cell types

Per- and poly-fluoroalkyl substances (PFAS) are extensively used in commerce leading to their prevalence in the environment. Due to their chemical stability, PFAS are considered to be persistent and bioaccumulative; they are frequently detected in both the environment and humans. Because of this, PFAS as a class (composed of hundreds to thousands of chemicals) are contaminants of very high concern. Little information is available for the vast majority of PFAS, and regulatory agencies lack safety data to determine whether exposure limits or restrictions are needed. Cell-based assays are a pragmatic approach to inform decision-makers on potential health hazards; therefore, we hypothesized that a targeted battery of human in vitro assays can be used to determine whether there are structure-bioactivity relationships for PFAS, and to characterize potential risks by comparing bioactivity (points of departure) to exposure estimates. We tested 56 PFAS from 8 structure-based subclasses in concentration response (0.1-100 μM) using six human cell types selected from target organs with suggested adverse effects of PFAS - human induced pluripotent stem cell (iPSC)-derived hepatocytes, neurons, and cardiomyocytes, primary human hepatocytes, endothelial and HepG2 cells. While many compounds were without effect; certain PFAS demonstrated cell-specific activity highlighting the necessity of using a compendium of in vitro models to identify potential hazards. No class-specific groupings were evident except for some chain length- and structure-related trends. In addition, margins of exposure (MOE) were derived using empirical and predicted exposure data. Conservative MOE calculations showed that most tested PFAS had a MOE in the 1-100 range; ∼20% of PFAS had MOE<1, providing tiered priorities for further studies. Overall, we show that a compendium of human cell-based models can be used to derive bioactivity estimates for a range of PFAS, enabling comparisons with human biomonitoring data. Furthermore, we emphasize that establishing structure-bioactivity relationships may be challenging for the tested PFAS.

Proof of concept for quantitative adverse outcome pathway modeling of chronic toxicity in repeated exposure

Adverse Outcome Pathway (AOP) is a useful tool to glean mode of action (MOE) of a chemical. However, in order to use it for the purpose of risk assessment, an AOP needs to be quantified using in vitro or in vivo data. Majority of quantitative AOPs developed so far, were for single exposure to progressively higher doses. Limited attempts were made to include time in the modeling. Here as a proof-of concept, we developed a hypothetical AOP, and quantified it using a virtual dataset for six repeated exposures using a Bayesian Network Analysis (BN) framework. The virtual data was generated using realistic assumptions. Effects of each exposure were analyzed separately using a static BN model and analyzed in combination using a dynamic BN (DBN) model. Our work shows that the DBN model can be used to calculate the probability of adverse outcome when other upstream KEs were observed earlier. These probabilities can help in identification of early indicators of AO. In addition, we also developed a data driven AOP pruning technique using a lasso-based subset selection, and show that the causal structure of AOP is itself dynamic and changes over time. This proof-of-concept study revealed the possibility for expanding the applicability of the AOP framework to incorporate biological dynamism in toxicity appearance by repeated insults.

Risk-based prioritization of PFAS using phenotypic and transcriptomic data from human induced pluripotent stem cell-derived hepatocytes and cardiomyocytes

Per- and polyfluoroalkyl substances (PFAS) are chemicals with important applications; they are persistent in the environment and may pose human health hazards. Regulatory agencies are con­sidering restrictions and bans of PFAS; however, little data exists for informed decisions. Several prioritization strategies were proposed for evaluation of potential hazards of PFAS. Structure-based grouping could expedite the selection of PFAS for testing; still, the hypothesis that structure-effect relationships exist for PFAS requires confirmation. We tested 26 structurally diverse PFAS from 8 groups using human induced pluripotent stem cell-derived hepatocytes and cardiomyocytes, and tested concentration-response effects on cell function and gene expression. Few phenotypic effects were observed in hepatocytes, but negative chronotropy was observed in cardiomyocytes for 8 PFAS. Substance- and cell type-dependent transcriptomic changes were more prominent but lacked substantial group-specific effects. In hepatocytes, we found upregulation of stress-related and extracellular matrix organization pathways, and down-regulation of fat metabolism. In car­diomyocytes, contractility-related pathways were most affected. We derived phenotypic and transcriptomic points of departure and compared them to predicted PFAS exposures. Conservative estimates for bioactivity and exposure were used to derive a bioactivity-to-exposure ratio (BER) for each PFAS; 23 of 26 PFAS had BER > 1. Overall, these data suggest that structure-based PFAS grouping may not be sufficient to predict their biological effects. Testing of individual PFAS may be needed for scientifically-supported decision-making. Our proposed strategy of using two human cell types and considering phenotypic and transcriptomic effects, combined with dose-response analysis and calculation of BER, may be used for PFAS prioritization.

Assessing proarrhythmic potential of environmental chemicals using a high throughput in vitro-in silico model with human induced pluripotent stem cell-derived cardiomyocytes

QT prolongation and the potentially fatal arrhythmia Torsades de Pointes are common causes for withdrawing or restricting drugs; however, little is known about similar liabilities of environmental chemicals. Current in vitro-in silico models for testing proarrhythmic liabilities, using human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CM), provide an opportunity to address this data gap. These methods are still low- to medium-throughput and not suitable for testing the tens of thousands of chemicals in commerce. We hypothesized that combining high-throughput population- based in vitro testing in hiPSC-CMs with a fully in silico data analysis workflow can offer sensitive and specific predictions of proarrhythmic potential. We calibrated the model with a published hiPSC-CM dataset of drugs known to be positive or negative for proarrhythmia and tested its performance using internal cross-validation and external validation. Additionally, we used computational down-sampling to examine three study designs for hiPSC-CM data: one replicate of one donor, five replicates of one donor, and one replicate of a population of five donors. We found that the population of five donors had the best performance for predicting proarrhythmic potential. The resulting model was then applied to predict the proarrhythmic potential of environmental chemicals, additionally characterizing risk through margin of exposure (MOE) calculations. Out of over 900 environmental chemicals tested, over 150 were predicted to have proarrhythmic potential, but only seven chemicals had a MOE < 1. We conclude that a high-throughput in vitro-in silico approach using population-based hiPSC-CM testing provides a reasonable strategy to screen environmental chemicals for proarrhythmic potential.

Advancing probabilistic risk assessment by integrating human biomonitoring, new approach methods, and Bayesian modeling: A case study with the mycotoxin deoxynivalenol

Deoxynivalenol (DON) is a mycotoxin frequently observed in cereals and cereal-based foods, with reported toxicological effects including reduced body weight, immunotoxicity and reproductive defects. The European Food Safety Authority used traditional risk assessment approaches to derive a deterministic Tolerable Daily Intake (TDI) of 1 μg/kg-day, however data from human biomarkers studies indicate widespread and variable exposure worldwide, necessitating more sophisticated and advanced methods to quantify population risk. The World Health Organization/International Programme on Chemical Safety (WHO/IPCS) has previously used DON as a case example in replacing the TDI with a probabilistic toxicity value, using default uncertainty and variability distributions to derive the Human Dose corresponding to an effect size M in the Ith percentile of the population (HDMI) for M = 5 % decrease in body weight and I = 1 %. In this study, we extend this case study by incorporating (1) Bayesian modeling approaches, (2) using both in vivo data and in vitro population new approach methods to replace default distributions for interspecies toxicokinetic (TK) differences and intraspecies TK and toxicodynamic (TD) variability, and (3) integrating biomonitoring data and probabilistic dose-response functions to characterize population risk distributions. We first derive an HDMI of 5.5 [1.4-24] μg/kg-day, also using TK modeling to converted the HDMI to Biomonitoring Equivalents, BEMI for comparison with biomonitoring data, with a blood BEMI of 0.53 [0.17-1.6] μg/L and a urinary excretion BEMI of 3.9 [1.0-16] μg/kg-day. We then illustrate how this integrative approach can advance quantitative risk characterization using two human biomonitoring datasets, estimating both the fraction of population with an effect size M ≥ 5 % as well as the distribution of effect sizes. Overall, we demonstrate that integration of Bayesian modeling, human biomonitoring data, and in vitro population-based TD data within the WHO/IPCS probabilistic framework yields more accurate, precise, and comprehensive risk characterization.

Cumulative Risk Meets Inter-Individual Variability: Probabilistic Concentration Addition of Complex Mixture Exposures in a Population-Based Human In Vitro Model

Although humans are continuously exposed to complex chemical mixtures in the environment, it has been extremely challenging to investigate the resulting cumulative risks and impacts. Recent studies proposed the use of “new approach methods,” in particular in vitro assays, for hazard and dose−response evaluation of mixtures. We previously found, using five human cell-based assays, that concentration addition (CA), the usual default approach to calculate cumulative risk, is mostly accurate to within an order of magnitude. Here, we extend these findings to further investigate how cell-based data can be used to quantify inter-individual variability in CA. Utilizing data from testing 42 Superfund priority chemicals separately and in 8 defined mixtures in a human cell-based population-wide in vitro model, we applied CA to predict effective concentrations for cytotoxicity for each individual, for “typical” (median) and “sensitive” (first percentile) members of the population, and for the median-to-sensitive individual ratio (defined as the toxicodynamic variability factor, TDVF). We quantified the accuracy of CA with the Loewe Additivity Index (LAI). We found that LAI varies more between different mixtures than between different individuals, and that predictions of the population median are generally more accurate than predictions for the “sensitive” individual or the TDVF. Moreover, LAI values were generally <1, indicating that the mixtures were more potent than predicted by CA. Together with our previous studies, we posit that new approach methods data from human cell-based in vitro assays, including multiple phenotypes in diverse cell types and studies in a population-wide model, can fill critical data gaps in cumulative risk assessment, but more sophisticated models of in vitro mixture additivity and bioavailability may be needed. In the meantime, because simple CA models may underestimate potency by an order of magnitude or more, either whole-mixture testing in vitro or, alternatively, more stringent benchmarks of cumulative risk indices (e.g., lower hazard index) may be needed to ensure public health protection.

A Population-Based Human In Vitro Approach to Quantify Inter-Individual Variability in Responses to Chemical Mixtures

Human cell-based population-wide in vitro models have been proposed as a strategy to derive chemical-specific estimates of inter-individual variability; however, the utility of this approach has not yet been tested for cumulative exposures in mixtures. This study aimed to test defined mixtures and their individual components and determine whether adverse effects of the mixtures were likely to be more variable in a population than those of the individual chemicals. The in vitro model comprised 146 human lymphoblastoid cell lines from four diverse subpopulations of European and African descent. Cells were exposed, in concentration−response, to 42 chemicals from diverse classes of environmental pollutants; in addition, eight defined mixtures were prepared from these chemicals using several exposure- or hazard-based scenarios. Points of departure for cytotoxicity were derived using Bayesian concentration−response modeling and population variability was quantified in the form of a toxicodynamic variability factor (TDVF). We found that 28 chemicals and all mixtures exhibited concentration−response cytotoxicity, enabling calculation of the TDVF. The median TDVF across test substances, for both individual chemicals or defined mixtures, ranged from a default assumption (101/2) of toxicodynamic variability in human population to >10. The data also provide a proof of principle for single-variant genome-wide association mapping for toxicity of the chemicals and mixtures, although replication would be necessary due to statistical power limitations with the current sample size. This study demonstrates the feasibility of using a set of human lymphoblastoid cell lines as an in vitro model to quantify the extent of inter-individual variability in hazardous properties of both individual chemicals and mixtures. The data show that population variability of the mixtures is unlikely to exceed that of the most variable component, and that similarity in genome-wide associations among components may be used to accrue additional evidence for grouping of constituents in a mixture for cumulative assessments.

Beyond pharmaceuticals: Fit-for-purpose new approach methodologies for environmental cardiotoxicity testing

Environmental factors play a substantial role in determining cardiovascular health, but data informing the risks presented by environmental toxicants is insufficient. In vitro new approach methodologies (NAMs) offer a promising approach with which to address the limitations of traditional in vivo and in vitro assays for assessing cardiotoxicity. Driven largely by the needs of pharmaceutical toxicity testing, considerable progress in developing NAMs for cardiotoxicity analysis has already been made. As the scientific and regulatory interest in NAMs for environmental chemicals continues to grow, a thorough understanding of the unique features of environmental cardiotoxicants and their associated cardiotoxicities is needed. Here, we review the key characteristics of as well as important regulatory and biological considerations for fit-for-purpose NAMs for environmental cardiotoxicity. By emphasizing the challenges and opportunities presented by NAMs for environmental cardiotoxicity we hope to accelerate their development, acceptance, and application.

In vitro-in silico-based prediction of inter-individual and inter-ethnic variations in the dose-dependent cardiotoxicity of R- and S-methadone in humans

New approach methodologies predicting human cardiotoxicity are of interest to support or even replace in vivo-based drug safety testing. The present study presents an in vitro–in silico approach to predict the effect of inter-individual and inter-ethnic kinetic variations in the cardiotoxicity of R- and S-methadone in the Caucasian and the Chinese population. In vitro cardiotoxicity data, and metabolic data obtained from two approaches, using either individual human liver microsomes or recombinant cytochrome P450 enzymes (rCYPs), were integrated with physiologically based kinetic (PBK) models and Monte Carlo simulations to predict inter-individual and inter-ethnic variations in methadone-induced cardiotoxicity. Chemical specific adjustment factors were defined and used to derive dose–response curves for the sensitive individuals. Our simulations indicated that Chinese are more sensitive towards methadone-induced cardiotoxicity with Margin of Safety values being generally two-fold lower than those for Caucasians for both methadone enantiomers. Individual PBK models using microsomes and PBK models using rCYPs combined with Monte Carlo simulations predicted similar inter-individual and inter-ethnic variations in methadone-induced cardiotoxicity. The present study illustrates how inter-individual and inter-ethnic variations in cardiotoxicity can be predicted by combining in vitro toxicity and metabolic data, PBK modelling and Monte Carlo simulations. The novel methodology can be used to enhance cardiac safety evaluations and risk assessment of chemicals.

A tiered approach to population-based in vitro testing for cardiotoxicity: Balancing estimates of potency and variability

Population-wide in vitro studies for characterization of cardiotoxicity hazard, risk, and population variability show that human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CM) are a powerful and high-throughput testing platform for drugs and environmental chemicals alike. However, studies in multiple donor-derived hiPSC-CMs, across large libraries of chemicals tested in concentration-response are technically complex, and study design optimization is needed to determine sufficient and fit-for-purpose population size considerations. Therefore, we tested a hypothesis that a computational down-sampling analysis based on the data from hiPSC-CM screening of 136 diverse compounds in a population of 43 non-diseased donors, including multiple replicates of the "standard" donor hiPSC-CMs, will inform optimal study designs depending on the decision context (hazard, risk and/or inter-individual variability in cardiotoxicity). Through 50 independent random subsamples of 5, 10, or 20 donors, we estimated accuracy and precision for quantifying potency, inter-individual variability, and QT prolongation risk; the results were compared to the full 43-donor cohort. We found that for potency and clinical risk of QT prolongation, a cohort of 5 randomly-selected unique donors provides accurate and precise estimates. Larger cohort sizes afforded marginal improvements, and 5 replicates of a single donor performed worse. For estimating inter-individual variability, cohorts of at least 20 donors are needed, with smaller populations on average showing bias towards underestimation in population variance. Collectively, this study shows that a variable-size hiPSC-CM-based population-wide in vitro model can be used in a number of decision scenarios for identifying cardiotoxic hazards of drugs and environmental chemicals in the population context.

Potential Human Health Hazard of Post-Hurricane Harvey Sediments in Galveston Bay and Houston Ship Channel: A Case Study of Using In Vitro Bioactivity Data to Inform Risk Management Decisions

Natural and anthropogenic disasters may be associated with redistribution of chemical contaminants in the environment; however, current methods for assessing hazards and risks of complex mixtures are not suitable for disaster response. This study investigated the suitability of in vitro toxicity testing methods as a rapid means of identifying areas of potential human health concern. We used sediment samples (n = 46) from Galveston Bay and the Houston Ship Channel (GB/HSC) areas after hurricane Harvey, a disaster event that led to broad redistribution of chemically-contaminated sediments, including deposition of the sediment on shore due to flooding. Samples were extracted with cyclohexane and dimethyl sulfoxide and screened in a compendium of human primary or induced pluripotent stem cell (iPSC)-derived cell lines from different tissues (hepatocytes, neuronal, cardiomyocytes, and endothelial) to test for concentration-dependent effects on various functional and cytotoxicity phenotypes (n = 34). Bioactivity data were used to map areas of potential concern and the results compared to the data on concentrations of polycyclic aromatic hydrocarbons (PAHs) in the same samples. We found that setting remediation goals based on reducing bioactivity is protective of both "known" risks associated with PAHs and "unknown" risks associated with bioactivity, but the converse was not true for remediation based on PAH risks alone. Overall, we found that in vitro bioactivity can be used as a comprehensive indicator of potential hazards and is an example of a new approach method (NAM) to inform risk management decisions on site cleanup.

A new approach method for characterizing inter-species toxicodynamic variability

Inter-species differences in toxicodynamics are often a critical source of uncertainty in safety evaluations and typically dealt with using default adjustment factors. In vitro studies that use cells from different species demonstrated some success for estimating the relationships between life span and/or body weight and sensitivity to cytotoxicity; however, no apparent investigation evaluated the utility of these models for risk assessment. It was hypothesized that an in vitro model using dermal fibroblasts derived from diverse species and individuals might be utilized to inform the extent of inter-species and inter-individual variability in toxicodynamics. To test this hypothesis and characterize both inter-species and inter-individual variability in cytotoxicity, concentration-response cytotoxicity screening of 40 chemicals in primary dermal fibroblasts from 68 individuals of 54 diverse species was conducted. Chemicals examined included drugs, environmental pollutants, and food/flavor/fragrance agents; most of these were previously assessed either in vivo or in vitro for inter-species or inter-individual variation. Species included humans, the typical preclinical species and representatives from other orders of mammals and birds. Data demonstrated that both inter-species and inter-individual components of variability contribute to the observed differences in sensitivity to cell death. Further, it was found that the magnitude of the observed inter-species and inter-individual differences was chemical-dependent. This study contributes to the paradigm shift in risk assessment from reliance on in vivo toxicity testing to higher-throughput in vitro or alternative approaches, extending the strategy to replace use of default adjustment factors with experimental characterization of toxicodynamic inter-individual variability and to also address toxicodynamic inter-species variability.

Cardiotoxicity Hazard and Risk Characterization of ToxCast Chemicals Using Human Induced Pluripotent Stem Cell-Derived Cardiomyocytes from Multiple Donors

Heart disease remains a significant human health burden worldwide with a significant fraction of morbidity attributable to environmental exposures. However, the extent to which the thousands of chemicals in commerce and the environment may contribute to heart disease morbidity is largely unknown, because in contrast to pharmaceuticals, environmental chemicals are seldom tested for potential cardiotoxicity. Human induced pluripotent stem cell (iPSC)-derived cardiomyocytes have become an informative in vitro model for cardiotoxicity testing of drugs with the availability of cells from multiple individuals allowing in vitro testing of population variability. In this study, we hypothesized that a panel of iPSC-derived cardiomyocytes from healthy human donors can be used to screen for the potential cardiotoxicity hazard and risk of environmental chemicals. We conducted concentration-response testing of 1029 chemicals (drugs, pesticides, flame retardants, polycyclic aromatic hydrocarbons (PAHs), plasticizers, industrial chemicals, food/flavor/fragrance agents, etc.) in iPSC-derived cardiomyocytes from 5 donors. We used kinetic calcium flux and high-content imaging to derive quantitative measures as inputs into Bayesian population concentration-response modeling of the effects of each chemical. We found that many environmental chemicals pose a hazard to human cardiomyocytes in vitro with more than half of all chemicals eliciting positive or negative chronotropic or arrhythmogenic effects. However, most of the tested environmental chemicals for which human exposure and high-throughput toxicokinetics data were available had wide margins of exposure and, thus, do not appear to pose a significant human health risk in a general population. Still, relatively narrow margins of exposure (<100) were estimated for some perfuoroalkyl substances and phthalates, raising concerns that cumulative exposures may pose a cardiotoxicity risk. Collectively, this study demonstrated the value of using a population-based human in vitro model for rapid, high-throughput hazard and risk characterization of chemicals for which little to no cardiotoxicity data are available from guideline studies in animals.

Human induced pluripotent stem cell (iPSC)-derived cardiomyocytes as an in vitro model in toxicology: strengths and weaknesses for hazard identification and risk characterization

INTRODUCTION

Human induced pluripotent stem cell (iPSC)-derived cardiomyocytes is one of the most widely used cell-based models that resulted from the discovery of how non-embryonic stem cells can be differentiated into multiple cell types. In just one decade, iPSC-derived cardiomyocytes went from a research lab to widespread use in biomedical research and preclinical safety evaluation for drugs and other chemicals.

AREAS COVERED

This manuscript reviews data on toxicology applications of human iPSC-derived cardiomyocytes. We detail the outcome of a systematic literature search on their use (i) in hazard assessment for cardiotoxicity liabilities, (ii) for risk characterization, (iii) as models for population variability, and (iv) in studies of personalized medicine and disease.

EXPERT OPINION

iPSC-derived cardiomyocytes are useful to increase the accuracy, precision, and efficiency of cardiotoxicity hazard identification for both drugs and non-pharmaceuticals, with recent efforts beginning to demonstrate their utility for risk characterization. Notable limitations include the needs to improve the maturation of cells in culture, to better understand their potential use identifying structural cardiotoxicity, and for additional case studies involving population-wide and disease-specific risk characterization. Ultimately, the greatest future benefits are likely for non-pharmaceutical chemicals, filling a critical gap where no routine testing for cardiotoxicity is currently performed.

Relationships between constituents of energy drinks and beating parameters in human induced pluripotent stem cell (iPSC)-Derived cardiomyocytes

Consumption of energy drinks has been associated with adverse cardiovascular effects; however, little is known about the ingredients that may contribute to these effects. We therefore characterized the chemical profiles and in vitro effects of energy drinks and their ingredients on human induced pluripotent stem cell (iPSC)-derived cardiomyocytes, and identified the putative active ingredients using a multivariate prediction model. Energy drinks from 17 widely-available over-the-counter brands were evaluated in this study. The concentrations of six common ingredients (caffeine, taurine, riboflavin, pantothenic acid, adenine, and L-methionine) were quantified by coupling liquid chromatography with a triple quadrupole mass spectrometer for the acquisition of LC-MS/MS spectra. In addition, untargeted analyses for each beverage were performed with a platform combining LC, ion mobility spectrometry and mass spectrometry (LC-IMS-MS) measurements. Approximately 300 features were observed across samples in the untargeted studies, and of these ~100 were identified. In vitro effects of energy drinks and some of their ingredients were then tested in iPSC-derived cardiomyocytes. Data on the beat rate (positive and negative chronotropy), ion channel function (QT prolongation), and cytotoxicity were collected in a dilution series. We found that some of the energy drinks elicited adverse effects on the cardiomyocytes with the most common being an increase in the beat rate, while QT prolongation was also observed at the lowest concentrations. Finally, concentration addition modeling using quantitative data from the 6 common ingredients and multivariate prediction modeling was used to determine potential ingredients responsible for the adverse effects on the cardiomyocytes. These analyses suggested theophylline, adenine, and azelate as possibly contributing to the in vitro effects of energy drinks on QT prolongation in cardiomyocytes.