Quantitative systems toxicology modeling in pharmaceutical research and development: An industry-wide survey and selected case study examples.

Quantitative systems toxicology (QST) models are increasingly being applied for predicting and understanding toxicity liabilities in pharmaceutical research and development. A European Federation of Pharmaceutical Industries and Associations (EFPIA)-wide survey was completed by 15 companies. The results provide insights into the current use of QST models across the industry.

Endeavours made by trade associations, pharmaceutical companies and regulators in the replacement, reduction and refinement of animal experimentation in safety testing of pharmaceuticals.

Following the European Commission decision to develop a roadmap to phase out animal testing and the signing of the US Modernisation Act, there is additional pressure on regulators and the pharmaceutical industry to abandon animal experimentation in safety testing. Often, endeavours already made by governments, regulators, trade associations, and industry to replace, reduce and refine animal experimentation (3Rs) are unnoticed. Herein, we review such endeavours to promote wider application and acceptance of 3Rs.

A dynamic model of the intestinal epithelium integrates multiple sources of preclinical data and enables clinical translation of drug-induced toxicity.

We have built a quantitative systems toxicology modeling framework focused on the early prediction of oncotherapeutic-induced clinical intestinal adverse effects. The model describes stem and progenitor cell dynamics in the small intestinal epithelium and integrates heterogeneous epithelial-related processes, such as transcriptional profiles, citrulline kinetics, and probability of diarrhea. We fitted a mouse-specific version of the model to quantify doxorubicin and 5-fluorouracil (5-FU)-induced toxicity, which included pharmacokinetics and 5-FU metabolism and assumed that both drugs led to cell cycle arrest and apoptosis in stem cells and proliferative progenitors.

Prediction of inotropic effect based on calcium transients in human iPSC-derived cardiomyocytes and machine learning.

Functional changes to cardiomyocytes are undesirable during drug discovery and identifying the inotropic effects of compounds is hence necessary to decrease the risk of cardiovascular adverse effects in the clinic. Recently, approaches leveraging calcium transients in human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) have been developed to detect contractility changes, induced by a variety of mechanisms early during drug discovery projects. Although these approaches have been able to provide some predictive ability, we hypothesised that using additional waveform parameters could offer improved insights, as well as predictivity.

Deriving waveform parameters from calcium transients in human iPSC-derived cardiomyocytes to predict cardiac activity with machine learning.

Human induced pluripotent stem cell-derived cardiomyocytes have been established to detect dynamic calcium transients by fast kinetic fluorescence assays that provide insights into specific aspects of clinical cardiac activity. However, the precise derivation and use of waveform parameters to predict cardiac activity merit deeper investigation. In this study, we derived, evaluated, and applied 38 waveform parameters in a novel Python framework, including (among others) peak frequency, peak amplitude, peak widths, and a novel parameter, shoulder-tail ratio.

Applying the CiPA approach to evaluate cardiac proarrhythmia risk of some antimalarials used off-label in the first wave of COVID-19.

We applied a set of in silico and in vitro assays, compliant with the Comprehensive In Vitro Proarrhythmia Assay (CiPA) paradigm, to assess the risk of chloroquine (CLQ) or hydroxychloroquine (OH-CLQ)-mediated QT prolongation and Torsades de Pointes (TdP), alone and combined with erythromycin (ERT) and azithromycin (AZI), drugs repurposed during the first wave of coronavirus disease 2019 (COVID-19). Each drug or drug combination was tested in patch clamp assays on seven cardiac ion channels, in in silico models of human ventricular electrophysiology (Virtual Assay) using control (healthy) or high-risk cell populations, and in human-induced pluripotent stem cell (hiPSC)-derived cardiomyocytes. In each assay, concentration-response curves encompassing and exceeding therapeutic free plasma levels were generated.

Considering discrepancy when calibrating a mechanistic electrophysiology model.

Uncertainty quantification (UQ) is a vital step in using mathematical models and simulations to take decisions. The field of cardiac simulation has begun to explore and adopt UQ methods to characterize uncertainty in model inputs and how that propagates through to outputs or predictions; examples of this can be seen in the papers of this issue. In this review and perspective piece, we draw attention to an important and under-addressed source of uncertainty in our predictions-that of uncertainty in the model structure or the equations themselves.

Four Ways to Fit an Ion Channel Model.

Mathematical models of ionic currents are used to study the electrophysiology of the heart, brain, gut, and several other organs. Increasingly, these models are being used predictively in the clinic, for example, to predict the risks and results of genetic mutations, pharmacological treatments, or surgical procedures. These safety-critical applications depend on accurate characterization of the underlying ionic currents.

Rapid Characterization of hERG Channel Kinetics II: Temperature Dependence.

Ion channel behavior can depend strongly on temperature, with faster kinetics at physiological temperatures leading to considerable changes in currents relative to room temperature. These temperature-dependent changes in voltage-dependent ion channel kinetics (rates of opening, closing, inactivating, and recovery) are commonly represented with Q coefficients or an Eyring relationship. In this article, we assess the validity of these representations by characterizing channel kinetics at multiple temperatures.

Reproducible model development in the cardiac electrophysiology Web Lab.

The modelling of the electrophysiology of cardiac cells is one of the most mature areas of systems biology. This extended concentration of research effort brings with it new challenges, foremost among which is that of choosing which of these models is most suitable for addressing a particular scientific question. In a previous paper, we presented our initial work in developing an online resource for the characterisation and comparison of electrophysiological cell models in a wide range of experimental scenarios.

Sinusoidal voltage protocols for rapid characterisation of ion channel kinetics.

Key Points: Ion current kinetics are commonly represented by current-voltage relationships, time constant-voltage relationships and subsequently mathematical models fitted to these. These experiments take substantial time, which means they are rarely performed in the same cell. Rather than traditional square-wave voltage clamps, we fitted a model to the current evoked by a novel sum-of-sinusoids voltage clamp that was only 8 s long.

Corrigendum: Optimization of an Cardiac Cell Model for Proarrhythmia Risk Assessment.

[This corrects the article on p. 616 in vol. 8, PMID: 28878692.

Uncertainty Quantification Reveals the Importance of Data Variability and Experimental Design Considerations for Proarrhythmia Risk Assessment.

The Comprehensive Proarrhythmia Assay (CiPA) is a global initiative intended to improve drug proarrhythmia risk assessment using a new paradigm of mechanistic assays. Under the CiPA paradigm, the relative risk of drug-induced Torsade de Pointes (TdP) is assessed using an model of the human ventricular action potential (AP) that integrates pharmacology data from multiple ion channels. Thus, modeling predictions of cardiac risk liability will depend critically on the variability in pharmacology data, and uncertainty quantification (UQ) must comprise an essential component of the assay.

Optimization of an Cardiac Cell Model for Proarrhythmia Risk Assessment.

Drug-induced Torsade-de-Pointes (TdP) has been responsible for the withdrawal of many drugs from the market and is therefore of major concern to global regulatory agencies and the pharmaceutical industry. The Comprehensive Proarrhythmia Assay (CiPA) was proposed to improve prediction of TdP risk, using models and multi-channel pharmacology data as integral parts of this initiative. Previously, we reported that combining dynamic interactions between drugs and the rapid delayed rectifier potassium current (IKr) with multi-channel pharmacology is important for TdP risk classification, and we modified the original O'Hara Rudy ventricular cell mathematical model to include a Markov model of IKr to represent dynamic drug-IKr interactions (IKr-dynamic ORd model).

Evaluation of an in silico cardiac safety assay: using ion channel screening data to predict QT interval changes in the rabbit ventricular wedge.

Introduction: Drugs that prolong the QT interval on the electrocardiogram present a major safety concern for pharmaceutical companies and regulatory agencies. Despite a range of assays performed to assess compound effects on the QT interval, QT prolongation remains a major cause of attrition during compound development. In silico assays could alleviate such problems.