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.

Inter-individual variation in chlorpyrifos toxicokinetics characterized by physiologically based kinetic (PBK) and Monte Carlo simulation comparing human liver microsome and Supersome™ cytochromes P450 (CYP)-specific kinetic data as model input.

The present study compares two approaches to evaluate the effects of inter-individual differences in the biotransformation of chlorpyrifos (CPF) on the sensitivity towards in vivo red blood cell (RBC) acetylcholinesterase (AChE) inhibition and to calculate a chemical-specific adjustment factor (CSAF) to account for inter-individual differences in kinetics (HKAF). These approaches included use of a Supersome™ cytochromes P450 (CYP)-based and a human liver microsome (HLM)-based physiologically based kinetic (PBK) model, both combined with Monte Carlo simulations. The results revealed that bioactivation of CPF exhibits biphasic kinetics caused by distinct differences in the Km of CYPs involved, which was elucidated by Supersome™ CYP rather than by HLM. Use of Supersome™ CYP-derived kinetic data was influenced by the accuracy of the intersystem extrapolation factors (ISEFs) required to scale CYP isoform activity of Supersome™ to HLMs. The predicted dose-response curves for average, 99th percentile and 1st percentile sensitive individuals were found to be similar in the two approaches when biphasic kinetics was included in the HLM-based approach, resulting in similar benchmark dose lower confidence limits for 10% inhibition (BMDL10) and HKAF values. The variation in metabolism-related kinetic parameters resulted in HKAF values at the 99th percentile that were slightly higher than the default uncertainty factor of 3.16. While HKAF values up to 6.9 were obtained when including also the variability in other influential PBK model parameters. It is concluded that the Supersome™ CYP-based approach appeared most adequate for identifying inter-individual variation in biotransformation of CPF and its resulting RBC AChE inhibition.

New approach methodologies (NAMs) for human-relevant biokinetics predictions. Meeting the paradigm shift in toxicology towards an animal-free chemical risk assessment

For almost fifteen years, the availability and regulatory acceptance of new approach methodologies (NAMs) to assess the absorption, distribution, metabolism and excretion (ADME/biokinetics) in chemical risk evaluations are a bottleneck. To enhance the field, a team of 24 experts from science, industry, and regulatory bodies, including new generation toxicologists, met at the Lorentz Centre in Leiden, The Netherlands. A range of possibilities for the use of NAMs for biokinetics in risk evaluations were formulated (for example to define species differences and human variation or to perform quantitative in vitro-in vivo extrapolations). To increase the regulatory use and acceptance of NAMs for biokinetics for these ADME considerations within risk evaluations, the development of test guidelines (protocols) and of overarching guidance documents is considered a critical step. To this end, a need for an expert group on biokinetics within the Organisation of Economic Cooperation and Development (OECD) to supervise this process was formulated. The workshop discussions revealed that method development is still required, particularly to adequately capture transporter mediated processes as well as to obtain cell models that reflect the physiology and kinetic characteristics of relevant organs. Developments in the fields of stem cells, organoids and organ-on-a-chip models provide promising tools to meet these research needs in the future.

Integrating in vitro data and physiologically based kinetic modeling-facilitated reverse dosimetry to predict human cardiotoxicity of methadone.

Development of novel testing strategies to detect adverse human health effects is of interest to replace in vivo-based drug and chemical safety testing. The aim of the present study was to investigate whether physiologically based kinetic (PBK) modeling-facilitated conversion of in vitro toxicity data is an adequate approach to predict in vivo cardiotoxicity in humans. To enable evaluation of predictions made, methadone was selected as the model compound, being a compound for which data on both kinetics and cardiotoxicity in humans are available. A PBK model for methadone in humans was developed and evaluated against available kinetic data presenting an adequate match. Use of the developed PBK model to convert concentration-response curves for the effect of methadone on human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CM) in the so-called multi electrode array (MEA) assay resulted in predictions for in vivo dose-response curves for methadone-induced cardiotoxicity that matched the available in vivo data. The results also revealed differences in protein plasma binding of methadone to be a potential factor underlying variation between individuals with respect to sensitivity towards the cardiotoxic effects of methadone. The present study provides a proof-of-principle of using PBK modeling-based reverse dosimetry of in vitro data for the prediction of cardiotoxicity in humans, providing a novel testing strategy in cardiac safety studies.

Integrating physiologically based kinetic (PBK) and Monte Carlo modelling to predict inter-individual and inter-ethnic variation in bioactivation and liver toxicity of lasiocarpine.

The aim of the present study was to predict the effect of inter-individual and inter-ethnic human kinetic variation on the sensitivity towards acute liver toxicity of lasiocarpine in the Chinese and the Caucasian population, and to derive chemical specific adjustment factors (CSAFs) by integrating variation in the in vitro kinetic constants Vmax and Km, physiologically based kinetic (PBK) modelling and Monte Carlo simulation. CSAFs were derived covering the 90th and 99th percentile of the population distribution of pyrrole glutathione adduct (7-GS-DHP) formation, reflecting bioactivation. The results revealed that in the Chinese population, as compared to the Caucasian population, the predicted 7-GS-DHP formation at the geometric mean, the 90th and the 99th percentile were 2.1-, 3.3- and 4.3-fold lower respectively. The CSAFs obtained using the 99th percentile values were 8.3, 17.0 and 19.5 in the Chinese, the Caucasian population and the two populations combined, respectively, while the CSAFs were generally 3.0-fold lower at the 90th percentile. These results indicate that when considering the formation of 7-GS-DHP the Caucasian population may be more sensitive towards acute liver toxicity of lasiocarpine, and further point out that the default safety factor of 3.16 for inter-individual human kinetic differences may not be sufficiently protective. Altogether, the results obtained demonstrate that integrating PBK modelling with Monte Carlo simulations using human in vitro data is a powerful strategy to quantify inter-individual variations in kinetics, and can be used to refine the human risk assessment of pyrrolizidine alkaloids.

Use of an in vitro-in silico testing strategy to predict inter-species and inter-ethnic human differences in liver toxicity of the pyrrolizidine alkaloids lasiocarpine and riddelliine.

Lasiocarpine and riddelliine are pyrrolizidine alkaloids (PAs) known to cause liver toxicity. The aim of this study was to predict the inter-species and inter-ethnic human differences in acute liver toxicity of lasiocarpine and riddelliine using physiologically based kinetic (PBK) modelling based reverse dosimetry of in vitro toxicity data. The concentration-response curves of in vitro cytotoxicity of lasiocarpine and riddelliine defined in pooled human hepatocytes were translated to in vivo dose-response curves by PBK models developed using kinetic data obtained from incubations with pooled tissue fractions from Chinese and Caucasian individuals, providing PBK models for the average Chinese and average Caucasian, respectively. From the predicted in vivo dose-response curves, the benchmark dose lower and upper confidence limits for 5% effect (BMDL5 and BMDU5) were derived and subsequently compared to those previously obtained in rat to evaluate inter-species differences. The inter-species differences amounted to 2.0-fold for lasiocarpine and 8.2-fold for riddelliine with humans being more sensitive than rats. The inter-ethnic human differences varied 2.0-fold for lasiocarpine and 5.0-fold for riddelliine with the average Caucasian being more sensitive than the average Chinese. In conclusion, the present study provides the proof-of-principle to predict inter-species and inter-ethnic differences in in vivo liver toxicity for PAs by an alternative testing strategy integrating in vitro cytotoxicity data with PBK modelling-based reverse dosimetry.

Development of a Combined In Vitro Physiologically Based Kinetic (PBK) and Monte Carlo Modelling Approach to Predict Interindividual Human Variation in Phenol-Induced Developmental Toxicity.

With our recently developed in vitro physiologically based kinetic (PBK) modelling approach, we could extrapolate in vitro toxicity data to human toxicity values applying PBK-based reverse dosimetry. Ideally information on kinetic differences among human individuals within a population should be considered. In the present study, we demonstrated a modelling approach that integrated in vitro toxicity data, PBK modelling and Monte Carlo simulations to obtain insight in interindividual human kinetic variation and derive chemical specific adjustment factors (CSAFs) for phenol-induced developmental toxicity. The present study revealed that UGT1A6 is the primary enzyme responsible for the glucuronidation of phenol in humans followed by UGT1A9. Monte Carlo simulations were performed taking into account interindividual variation in glucuronidation by these specific UGTs and in the oral absorption coefficient. Linking Monte Carlo simulations with PBK modelling, population variability in the maximum plasma concentration of phenol for the human population could be predicted. This approach provided a CSAF for interindividual variation of 2.0 which covers the 99th percentile of the population, which is lower than the default safety factor of 3.16 for interindividual human kinetic differences. Dividing the dose-response curve data obtained with in vitro PBK-based reverse dosimetry, with the CSAF provided a dose-response curve that reflects the consequences of the interindividual variability in phenol kinetics for the developmental toxicity of phenol. The strength of the presented approach is that it provides insight in the effect of interindividual variation in kinetics for phenol-induced developmental toxicity, based on only in vitro and in silico testing.

Integrating in vitro data and physiologically based kinetic (PBK) modelling to assess the in vivo potential developmental toxicity of a series of phenols.

Toxicity outcomes derived in vitro do not always reflect in vivo toxicity values, which was previously observed for a series of phenols tested in the embryonic stem cell test (EST). Translation of in vitro data to the in vivo situation is therefore an important, but still limiting step for the use of in vitro toxicity outcomes in the safety assessment of chemicals. The aim of the present study was to translate in vitro embryotoxicity data for a series of phenols to in vivo developmental toxic potency values for the rat by physiologically based kinetic (PBK) modelling-based reverse dosimetry. To this purpose, PBK models were developed for each of the phenols. The models were parameterised with in vitro-derived values defining metabolism and transport of the compounds across the intestinal and placental barrier and with in silico predictions and data from the literature. Using PBK-based reverse dosimetry, in vitro concentration-response curves from the EST were translated into in vivo dose-response curves from which points of departure (PoDs) were derived. The predicted PoDs differed less than 3.6-fold from PoDs derived from in vivo toxicity data for the phenols available in the literature. Moreover, the in vitro PBK-based reverse dosimetry approach could overcome the large disparity that was observed previously between the in vitro and the in vivo relative potency of the series of phenols. In conclusion, this study shows another proof-of-principle that the in vitro PBK approach is a promising strategy for non-animal-based safety assessment of chemicals.

Combining in vitro embryotoxicity data with physiologically based kinetic (PBK) modelling to define in vivo dose-response curves for developmental toxicity of phenol in rat and human.

In vitro assays are often used for the hazard characterisation of compounds, but their application for quantitative risk assessment purposes is limited. This is because in vitro assays cannot provide a complete in vivo dose-response curve from which a point of departure (PoD) for risk assessment can be derived, like the no observed adverse effect level (NOAEL) or the 95 % lower confidence limit of the benchmark dose (BMDL). To overcome this constraint, the present study combined in vitro data with a physiologically based kinetic (PBK) model applying reverse dosimetry. To this end, embryotoxicity of phenol was evaluated in vitro using the embryonic stem cell test (EST), revealing a concentration-dependent inhibition of differentiation into beating cardiomyocytes. In addition, a PBK model was developed on the basis of in vitro and in silico data and data available from the literature only. After evaluating the PBK model performance, effective concentrations (ECx) obtained with the EST served as an input for in vivo plasma concentrations in the PBK model. Applying PBK-based reverse dosimetry provided in vivo external effective dose levels (EDx) from which an in vivo dose-response curve and a PoD for risk assessment were derived. The predicted PoD lies within the variation of the NOAELs obtained from in vivo developmental toxicity data from the literature. In conclusion, the present study showed that it was possible to accurately predict a PoD for the risk assessment of phenol using in vitro toxicity data combined with reverse PBK modelling.

Relative embryotoxic potency of p-substituted phenols in the embryonic stem cell test (EST) and comparison to their toxic potency in vivo and in the whole embryo culture (WEC) assay.

The applicability of the embryonic stem cell test (EST) as an alternative for in vivo embryotoxicity testing was evaluated for a series of five p-substituted phenols. To this purpose, the potency ranking for this class of compounds derived from the inhibition of cardiomyocyte differentiation in the EST was compared to in vivo embryotoxic potency data obtained from literature and to the potency ranking defined in the in vitro whole embryo culture (WEC) assay. From the results obtained it appears that the EST was able to identify the embryotoxic potential for p-substituted phenols, providing an identical potency ranking compared to the WEC assay. However, the EST was not able to predict an accurate ranking for the phenols compared to their potency observed in vivo. Only phenol, the least potent compound within this series, was correctly ranked. Furthermore, p-mercaptophenol was correctly identified as a relative potent congener of the phenols tested, but its ranking was distorted by p-heptyloxyphenol, of which the toxicity was overestimated in the EST. It is concluded that when attempting to explain the observed disparity in potency rankings between in vitro and in vivo embryotoxicity, the in vitro models should be combined with a kinetic model describing in vivo absorption, distribution, metabolism and excretion processes of the compounds.