Can 3D Printed Tablets Be Bioequivalent and How to Test It: A PBPK Model Based Virtual Bioequivalence Study for Ropinirole Modified Release Tablets.

As the field of personalized dosing develops, the pharmaceutical manufacturing industry needs to offer flexibility in terms of tailoring the drug release and strength to the individual patient's needs. One of the promising tools which have such capacity is 3D printing technology. However, manufacturing small batches of drugs for each patient might lead to huge test burden, including the need to conduct bioequivalence trials of formulations to support the change of equipment or strength. In this paper we demonstrate how to use 3D printing in conjunction with virtual bioequivalence trials based on physiologically based pharmacokinetic (PBPK) modeling. For this purpose, we developed 3D printed ropinirole formulations and tested their bioequivalence with the reference product Polpix. The Simcyp simulator and previously developed ropinirole PBPK model were used for the clinical trial simulations. The Weibull-fitted dissolution profiles of test and reference formulations were used as inputs for the model. The virtual bioequivalence trials were run using parallel design. The study power of 80% was reached using 125 individuals. The study demonstrated how to use PBPK modeling in conjunction with 3D printing to test the virtual bioequivalence of newly developed formulations. This virtual experiment demonstrated the bioequivalence of one of the newly developed formulations with a reference product available on a market.

3D printing combined with biopredictive dissolution and PBPK/PD modeling optimization and personalization of pharmacotherapy: Are we there yet?

Precision medicine requires selecting the appropriate dosage regimen for a patient using the right drug, at the right time. Model-Informed Precision Dosing (MIPD) is a concept suggesting utilization of model-based prediction methods for optimizing the treatment benefit-harm balance, based on individual characteristics of the patient, disease, treatment method, and other factors. Here, we discuss a theoretical workflow comprising several elements, beginning from the physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) models, through 3D printed tablets with the model proposed dose, information range and flow, and the patient themselves. We also describe each of these elements, and the connection between them, highlighting challenges and potential obstacles.

Artificial Intelligence That Predicts Sensitizing Potential of Cosmetic Ingredients with Accuracy Comparable to Animal and In Vitro Tests-How Does the Infotechnomics Compare to Other "Omics" in the Cosmetics Safety Assessment?

The aim of the current study was to develop an in silico model to predict the sensitizing potential of cosmetic ingredients based on their physicochemical characteristics and to compare the predictions with historical animal data and results from "omics"-based in vitro studies. An in silico model was developed with the use of WEKA machine learning software fed with physicochemical and structural descriptors of haptens and trained with data from published epidemiological studies compiled into estimated odds ratio (eOR) and estimated attributable risk (eAR) indices. The outcome classification was compared to the results of animal studies and in vitro tests. Of all the models tested, the best results were obtained for the Naive Bayes classifier trained with 24 physicochemical descriptors and eAR, which yielded an accuracy of 86%, sensitivity of 80%, and specificity of 90%. This model was subsequently used to predict the sensitizing potential of 15 emerging and less-studied haptens, of which 7 were classified as sensitizers: cyclamen aldehyde, N,N-dimethylacrylamide, dimethylthiocarbamyl benzothiazole sulphide, geraniol hydroperoxide, isobornyl acrylate, neral, and prenyl caffeate. The best-performing model (NaiveBayes eAR, 24 parameters), along with an alternative model based on eOR (Random Comittee eOR, 17 parameters), are available for further tests by interested readers. In conclusion, the proposed infotechnomics approach allows for a prediction of the sensitizing potential of cosmetic ingredients (and possibly also other haptens) with accuracy comparable to historical animal tests and in vitro tests used nowadays. In silico models consume little resources, are free of ethical concerns, and can provide results for multiple chemicals almost instantly; therefore, the proposed approach seems useful in the safety assessment of cosmetics.

Data on ADME parameters of bisphenol A and its metabolites for use in physiologically based pharmacokinetic modelling.

The paper presents the collection of physicochemical parameters of bisphenol A (BPA) and its sulfate (BPAS) and glucuronide (BPAG) conjugates, accompanied by data characterizing their absorption, distribution, metabolism and excretion (ADME) behavior following oral administration of BPA. The data were collected from open literature sources and publicly available databases. Additionally, data calculated by using the MarvinSketch 18.30.0 software or predicted by relevant QSAR models built in Simcyp® Simulator were also used. All data were analysed and are fit for purpose if necessary to ensure a reliable prediction of pharmacokinetics of BPA and its conjugates. The data selection process and reasoning for fitting is provided to allow critical assessment and to ensure data transparency. Finally, the sensitivity analysis was performed to assess the influence of the selected parameters on the PBPK model predictions.

Physiologically based modelling of dermal absorption and kinetics of consumer-relevant chemicals: A case study with exposure to bisphenol A from thermal paper.

Bisphenol A (BPA) is one of the best studied industrial chemicals in terms of exposure, toxicity, and toxicokinetics. This renders it an ideal candidate to exploit the recent advancements in physiologically based pharmacokinetic (PBPK) modelling to support risk assessment of BPA specifically, and of other consumer-relevant hazardous chemicals in general. Using the exposure from thermal paper as a case scenario, this study employed the multi-phase multi-layer mechanistic dermal absorption (MPML MechDermA) model available in the Simcyp® Simulator to simulate the dermal toxicokinetics of BPA at local and systemic levels. Sensitivity analysis helped to identify physicochemical and physiological factors influencing the systemic exposure to BPA. The iterative modelling process was as follows: (i) development of compound files for BPA and its conjugates, (ii) setting-up of a PBPK model for intravenous administration, (iii) extension for oral administration, and (iv) extension for exposure via skin (i.e., hand) contact. A toxicokinetic study involving hand contact to BPA-containing paper was used for model refinement. Cumulative urinary excretion of total BPA had to be employed for dose reconstruction. PBPK model performance was verified using the observed serum BPA concentrations. The predicted distribution across the skin compartments revealed a depot of BPA in the stratum corneum (SC). These findings shed light on the role of the SC to act as temporary reservoir for lipophilic chemicals prior to systemic absorption, which inter alia is relevant for the interpretation of human biomonitoring data and for establishing the relationship between external and internal measures of exposure.

An open-access data set of pig skin anatomy and physiology for modelling purposes.

The use of animal as opposed to human skin for in vitro permeation testing (IVPT) is an alternative, which can reduce logistical and economic issues. However, this surrogate also has ethical considerations and may not provide an accurate estimation of dermal absorption in humans due to physiological differences. The current project aimed to provide a detailed repository for the anatomical and physiological parameters of porcine skin, with the aim of parametrizing the Multi-phase Multi-layer Mechanistic Dermal Absorption (MPML MechDermA) Model in the Simcyp Simulator. The MPML MechDermA Model is a physiologically based pharmacokinetic (PBPK) model that accounts for the physiology and geometry of skin in a mechanistic mathematical modelling framework. The database provided herein contains information on 14 parameters related to porcine skin anatomy and physiology, namely, skin surface pH, number of stratum corneum (SC) layers, SC thickness, corneocyte thickness, corneocyte dimensions (length and width), volume fraction of water in corneocyte (where SC is divided into four parts with different water contents), intercellular lipid thickness, viable epidermis thickness, dermis thickness, hair follicle and hair shaft diameter, hair follicle depth and hair follicle density. The collected parameters can be used to parameterize PBPK models, which could be further utilized to bridge the gap between animal and human studies with interspecies extrapolation or to predict dermatokinetic properties typically assessed in IVPT experiments. Database URL: https://data.mendeley.com/datasets/mwz9xv4cpd/1.

Utilization of mechanistic modelling and simulation to analyse fenspiride proarrhythmic potency - Role of physiological and other non-drug related parameters.

WHAT IS KNOWN AND OBJECTIVE: Fenspiride, a drug that had been used for decades for the treatment of respiratory diseases, was recently withdrawn from the market due to the potential risk of QT prolongation and proarrhythmia. This is the first such withdrawal for many years and hence poses a question whether such risk could have been predicted and to what degree non-drug-specific parameters play a role in the reported QT prolongation and cases of TdP. The study aim was to test various 'what-if' scenarios to assess the influence of age, gender, heart rate, and plasma potassium concentration on QT interval prolongation due to various doses of fenspiride with the use of mechanistic mathematical modelling. METHODS: Concentration-time profiles were simulated with the use of a PBPK model developed based on published physico-chemical data, data from in vitro ADME experiments, and in vivo PK study results. Pharmacodynamic effect, that is, drug-triggered pseudoECG signal modification was simulated using a biophysically detailed model of human cardiac myocytes. Analysis of the qNet metric was also performed to classify proarrhythmic risk related to fenspiride. RESULTS: In the simulation study, arrhythmia was not observed even in the 'what-if' scenarios with extreme exposure, age, heart rate, and plasma potassium concentration. The qNet metric value positioned fenspiride in the intermediate risk class. WHAT IS NEW AND CONCLUSION: It can be hypothesized that the clinically observed arrhythmia cases were not directly caused by fenspiride alone but a combination of multiple factors, including comedications.

Development and Performance Verification of the PBPK Model for Antazoline and Its Metabolite and Its Utilization for Pharmacological Hypotheses Formulating.

Antazoline is an antihistaminic drug that is effective in the termination of paroxysmal atrial fibrillation. Despite its long presence in the market, antazoline's ADME parameters and pharmacokinetic effects in humans are poorly characterized. The objective of this study was to fill this gap by generation of in vitro and in vivo data and the development of a physiologically based pharmacokinetic model describing antazoline and its main metabolite disposition. A set of ADME parameters for the antazoline and its hydroxy metabolite is provided based on literature data, QSAR predictions, in vitro binding and metabolic stability assays. These can be used to feed PBPK models. In our current work, the developed PBPK model simulating simultaneously the pharmacokinetic profile of antazoline and its metabolite was successfully verified against the available clinical data and the presented capability to account for the clinically observed variability. When used to feed the PD model (e.g., simulating ECG), concentration-time profiles predicted by the model enable the assessment of antazoline's effect in various clinical scenarios with the possibility to account for population differences or CP mediated drug-drug interactions.

Open-access database of literature derived drug-related Torsade de Pointes cases.

Since an introduction of an ICH guidance in 2005, no new drugs were withdrawn from the market because of the causation of Torsade de Pointes (TdP). However, the risk of TdP is still a concern for marketed drugs. TdP is a type of polymorphic ventricular tachycardia which may lead to sudden cardiac death. QT/QTc interval prolongation is considered a sensitive, but not specific biomarker. To improve the effectiveness of studies' workflow related to TdP risk prediction we created an extensive, structured, open-access database of drug-related TdP cases. PubMed, Google Scholar bibliographic databases, and the Internet, via the Google search engine, were searched to identify eligible reports. A total of 424 papers with a description of 634 case reports and observational studies were included. Each paper was manually examined and listed with up to 53 variables related to patient/population characteristics, general health parameters, used drugs, laboratory measurements, ECG results, clinical management, and its outcomes, as well as suspected drug's properties and its FDA adverse reaction reports. The presented database may be considered as an extension of the recently developed and published database of drug cardiac safety-related information, part of the tox-portal project providing resources for cardiac toxicity assessment.

Magnetic Core-Shell Molecularly Imprinted Nano-Conjugates for Extraction of Antazoline and Hydroxyantazoline from Human Plasma-Material Characterization, Theoretical Analysis and Pharmacokinetics.

The aim of this study was to develop magnetic molecularly imprinted nano-conjugate sorbent for effective dispersive solid phase extraction of antazoline (ANT) and its metabolite, hydroxyantazoline (ANT-OH) in analytical method employing liquid chromatography coupled with mass spectrometry method. The core-shell material was characterized in terms of adsorption properties, morphology and structure. The heterogeneous population of adsorption sites towards ANT-OH was characterized by two Kd and two Bmax values: Kd (1) = 0.319 µg L-1 and Bmax (1) = 0.240 µg g-1, and Kd (2) = 34.6 µg L-1 and Bmax (2) = 5.82 µg g-1. The elemental composition of magnetic sorbent was as follows: 17.55, 37.33, 9.14, 34.94 wt% for Si, C, Fe and O, respectively. The extraction protocol was optimized, and the obtained results were explained using theoretical analysis. Finally, the analytical method was validated prior to application to pharmacokinetic study in which the ANT was administrated intravenously to three healthy volunteers. The results prove that the novel sorbent could be useful in extraction of ANT and ANT-OH from human plasma and that the analytical strategy could be a versatile tool to explain a potential and pharmacological activity of ANT and ANT-OH.

How circadian variability of the heart rate and plasma electrolytes concentration influence the cardiac electrophysiology - model-based case study.

The circadian rhythm of cardiac electrophysiology is dependent on many physiological and biochemical factors. Provided, that models describing the circadian patterns of cardiac activity and/or electrophysiology which have been verified to the acceptable level, modeling and simulation can give answers to many of heart chronotherapy questions. The aim of the study was to assess the performance of the circadian models implemented in Cardiac Safety Simulator v 2.2 (Certara, Sheffield, UK) (CSS), as well as investigate the influence ofcircadian rhythms on the simulation results in terms of cardiac safety. The simulations which were run in CSS accounted for inter-individual and intra-individual variability. Firstly, the diurnal variations in QT interval length in a healthy population were simulated accounting for heart rate (HR) circadian changes alone, or with concomitant diurnal variations of plasma ion concentrations. Next, tolterodine was chosen as an exemplary drug for PKPD modelling exercise to assess the role of circadian rhythmicity in the prediction of drug effects on QT interval. The results of the simulations were in line with clinical observations, what can serve as a verification of the circadian models implemented in CSS. Moreover, the results have suggested that the circadian variability of the electrolytes balance is the main factor influencing QT circadian pattern. The fluctuation of ion concentration increases the intra-subject variability of predicted drug-triggered QT corrected for HR (QTc) prolongation effect and, in case of modest drug effect on QTc interval length, allows to capture this effect.

Characterization of In Vitro and In Vivo Metabolism of Antazoline Using Liquid Chromatography-Tandem Mass Spectrometry.

Antazoline (ANT) was recently shown to be an effective and safe antiarrhythmic drug in the termination of atrial fibrillation. However, the drug is still not listed in clinical guidelines. No data on ANT metabolism in humans is available. We used liquid chromatography coupled with tandem mass spectrometry to identify and characterize metabolites of ANT. We analyzed plasma of volunteers following a single intravenous administration of 100 mg of ANT mesylate and in in vitro cultures of human hepatocytes. We revealed that ANT was transformed into at least 15 metabolites and we investigated the role of cytochrome P450 isoforms. CYP2D6 was the main one involved in the fast metabolism of ANT. The biotransformation of ANT by CYP2C19 was much slower. The main Phase I metabolite was M1 formed by the removal of phenyl and metabolite M2 with hydroxyl in the para position of phenyl. Glucuronidation was the leading Phase II metabolism. Further study on pharmacokinetics of the metabolites would allow us to better understand the activity profile of ANT and to predict its potential clinical applications. Ultimately, further investigation of the activity profile of the new hydroxylated M2 metabolite of ANT might result in an active substance with a different pharmacological profile than the parent molecule, and potentially a new drug candidate.

Better prediction of the local concentration-effect relationship: the role of physiologically based pharmacokinetics and quantitative systems pharmacology and toxicology in the evolution of model-informed drug discovery and development.

Model-informed drug discovery and development (MID3) is an umbrella term under which sit several computational approaches: quantitative systems pharmacology (QSP), quantitative systems toxicology (QST) and physiologically based pharmacokinetics (PBPK). QSP models are built using mechanistic knowledge of the pharmacological pathway focusing on the putative mechanism of drug efficacy; whereas QST models focus on safety and toxicity issues and the molecular pathways and networks that drive these adverse effects. These can be mediated through exaggerated on-target or off-target pharmacology, immunogenicity or the physiochemical nature of the compound. PBPK models provide a mechanistic description of individual organs and tissues to allow the prediction of the intra- and extra-cellular concentration of the parent drug and metabolites under different conditions. Information on biophase concentration enables the prediction of a drug effect in different organs and assessment of the potential for drug-drug interactions. Together, these modelling approaches can inform the exposure-response relationship and hence support hypothesis generation and testing, compound selection, hazard identification and risk assessment through to clinical proof of concept (POC) and beyond to the market.

Virtual Thorough QT (TQT) Trial-Extrapolation of In Vitro Cardiac Safety Data to In Vivo Situation Using Multi-Scale Physiologically Based Ventricular Cell-wall Model Exemplified with Tolterodine and Fesoterodine.

QT interval prolongation typically assessed with dedicated clinical trials called thorough QT/QTc (TQT) studies is used as surrogate to identify the proarrhythmic risk of drugs albeit with criticism in terms of cost-effectiveness in establishing the actual risk of torsade de pointes (TdP). Quantitative systems toxicology and safety (QSTS) models have potential to quantitatively translate the in vitro cardiac safety data to clinical level including simulation of TQT trials. Virtual TQT simulations have been exemplified with use of two related drugs tolterodine and fesoterodine. The impact of bio-relevant concentration in plasma versus estimated heart tissue exposure on predictions was also assessed. Tolterodine and its therapeutically equipotent metabolite formed via CYP2D6 pathway, 5-HMT, inhibit multiple cardiac ion currents (IKr, INa, ICaL). The QSTS model was able to accurately simulate the QT prolongation at therapeutic and supra-therapeutic dose levels of tolterodine well within 95% confidence interval limits of observed data. The model was able to predict the QT prolongation difference between CYP2D6 extensive and poor metaboliser subject groups at both dose levels thus confirming the ability of the model to account for electrophysiologically active metabolite. The QSTS model was able to simulate the negligible QT prolongation observed with fesoterodine establishing that the 5-HMT does not prolong QT interval even though it is a blocker of hERG channel. With examples of TOL and FESO, we demonstrated the utility of the QSTS approaches to simulate virtual TQT trials, which in turn could complement and reduce the clinical studies or help optimise clinical trial designs.

Quantitative Assessment of the Physiological Parameters Influencing QT Interval Response to Medication: Application of Computational Intelligence Tools.

Human heart electrophysiology is complex biological phenomenon, which is indirectly assessed by the measured ECG signal. ECG trace is further analyzed to derive interpretable surrogates including QT interval, QRS complex, PR interval, and T wave morphology. QT interval and its modification are the most commonly used surrogates of the drug triggered arrhythmia, but it is known that the QT interval itself is determined by other nondrug related parameters, physiological and pathological. In the current study, we used the computational intelligence algorithms to analyze correlations between various simulated physiological parameters and QT interval. Terfenadine given concomitantly with 8 enzymatic inhibitors was used as an example. The equation developed with the use of genetic programming technique leads to general reasoning about the changes in the prolonged QT. For small changes of the QT interval, the drug-related IKr and ICa currents inhibition potentials have major impact. The physiological parameters such as body surface area, potassium, sodium, and calcium ions concentrations are negligible. The influence of the physiological variables increases gradually with the more pronounced changes in QT. As the significant QT prolongation is associated with the drugs triggered arrhythmia risk, analysis of the role of physiological parameters influencing ECG seems to be advisable.

Towards Bridging Translational Gap in Cardiotoxicity Prediction: an Application of Progressive Cardiac Risk Assessment Strategy in TdP Risk Assessment of Moxifloxacin.

Drug-induced cardiac arrhythmia, especially occurrence of torsade de pointes (TdP), has been a leading cause of attrition and post-approval re-labeling and withdrawal of many drugs. TdP is a multifactorial event, reflecting more than just drug-induced cardiac ion channel inhibition and QT interval prolongation. This presents a translational gap in extrapolating pre-clinical and clinical cardiac safety assessment to estimate TdP risk reliably, especially when the drug of interest is used in combination with other QT-prolonging drugs for treatment of diseases such as tuberculosis. A multi-scale mechanistic modeling framework consisting of physiologically based pharmacokinetics (PBPK) simulations of clinically relevant drug exposures combined with Quantitative Systems Toxicology (QST) models of cardiac electro-physiology could bridge this gap. We illustrate this PBPK-QST approach in cardiac risk assessment as exemplified by moxifloxacin, an anti-tuberculosis drug with abundant clinical cardiac safety data. PBPK simulations of moxifloxacin concentrations (systemic circulation and estimated in heart tissue) were linked with in vitro measurements of cardiac ion channel inhibition to predict the magnitude of QT prolongation in healthy individuals. Predictions closely reproduced the clinically observed QT interval prolongation, but no arrhythmia was observed, even at ×10 exposure. However, the same exposure levels in presence of physiological risk factors, e.g., hypokalemia and tachycardia, led to arrhythmic event in simulations, consistent with reported moxifloxacin-related TdP events. Application of a progressive PBPK-QST cardiac risk assessment paradigm starting in early development could guide drug development decisions and later define a clinical "safe space" for post-approval risk management to identify high-risk clinical scenarios.

Drug-physiology interaction and its influence on the QT prolongation-mechanistic modeling study.

The current study is an example of drug-disease interaction modeling where a drug induces a condition which can affect the pharmacodynamics of other concomitantly taken drugs. The electrophysiological effects of hypokaliemia and heart rate changes induced by the antiasthmatic drugs were simulated with the use of the cardiac safety simulator. Biophysically detailed model of the human cardiac physiology-ten Tusscher ventricular cardiomyocyte cell model-was employed to generate pseudo-ECG signals and QTc intervals for 44 patients from four clinical studies. Simulated and observed mean QTc values with standard deviation (SD) for each reported study point were compared and differences were analyzed with Student's t test (α = 0.05). The simulated results reflected the QTc interval changes measured in patients, as well as their clinically observed interindividual variability. The QTc interval changes were highly correlated with the change in plasma potassium both in clinical studies and in the simulations (Pearson's correlation coefficient > 0.55). The results suggest that the modeling and simulation approach could provide valuable quantitative insight into the cardiological effect of the potassium and heart rate changes caused by electrophysiologically inactive, non-cardiological drugs. This allows to simulate and predict the joint effect of several risk factors for QT prolongation, e.g., drug-dependent QT prolongation due to the ion channels inhibition and the current patient physiological conditions.

Drug interaction at hERG channel: In vitro assessment of the electrophysiological consequences of drug combinations and comparison against theoretical models.

Drugs carry a proarrhythmic risk, which gets even greater when they are used in combination. In vitro assessment of the proarrhythmic potential of drugs is limited to one compound and thus neglects the potential of drug-drug interactions, including those involving active metabolites. Here we present the results of an in vitro study of potential drug-drug interactions at the level of the hERG channel for the combination of up to three compounds: loratadine, desloratadine and ketoconazole. Experiments were performed at room temperature on an automated patch-clamp device CytoPatch 2, with the use of heterogeneously, stably transfected HEK cells. Single drugs, pairs and triplets were used. The results provided as the inhibition of the IKr current for pairs were compared against the calculated theoretical interaction. Models applied to calculate the combined effect of inhibitory actions of simultaneously given drugs include: (1) simple additive model with a maximal inhibition limit of 1 (all channels blocked in 100%); (2) Bliss independence; and (3) Loewe additivity. The observed IC50 values for loratadine, desloratadine and ketoconazole were 5.15, 1.95 and 0.74 µm respectively. For the combination of drugs tested in pairs, the effect was concentration dependent. In lower concentrations, the synergistic effect was observed, while for the highest tested concentrations it was subadditive. To triple the effect, it was subadditive regardless of concentrations. The square root of sum of squares of differences between the observed and predicted total inhibition was calculated to assess the theoretical interaction models. For most of the drugs, the allotopic model offered the best fit.

Real Patient and its Virtual Twin: Application of Quantitative Systems Toxicology Modelling in the Cardiac Safety Assessment of Citalopram.

A quantitative systems toxicology (QST) model for citalopram was established to simulate, in silico, a 'virtual twin' of a real patient to predict the occurrence of cardiotoxic events previously reported in patients under various clinical conditions. The QST model considers the effects of citalopram and its most notable electrophysiologically active primary (desmethylcitalopram) and secondary (didesmethylcitalopram) metabolites, on cardiac electrophysiology. The in vitro cardiac ion channel current inhibition data was coupled with the biophysically detailed model of human cardiac electrophysiology to investigate the impact of (i) the inhibition of multiple ion currents (IKr, IKs, ICaL); (ii) the inclusion of metabolites in the QST model; and (iii) unbound or total plasma as the operating drug concentration, in predicting clinically observed QT prolongation. The inclusion of multiple ion channel current inhibition and metabolites in the simulation with unbound plasma citalopram concentration provided the lowest prediction error. The predictive performance of the model was verified with three additional therapeutic and supra-therapeutic drug exposure clinical cases. The results indicate that considering only the hERG ion channel inhibition of only the parent drug is potentially misleading, and the inclusion of active metabolite data and the influence of other ion channel currents should be considered to improve the prediction of potential cardiac toxicity. Mechanistic modelling can help bridge the gaps existing in the quantitative translation from preclinical cardiac safety assessment to clinical toxicology. Moreover, this study shows that the QST models, in combination with appropriate drug and systems parameters, can pave the way towards personalised safety assessment.