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.

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 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.

Virtual Clinical Trial Toward Polytherapy Safety Assessment: Combination of Physiologically Based Pharmacokinetic/Pharmacodynamic-Based Modeling and Simulation Approach With Drug-Drug Interactions Involving Terfenadine as an Example.

A Quantitative Systems Pharmacology approach was utilized to predict the cardiac consequences of drug-drug interaction (DDI) at the population level. The Simcyp in vitro-in vivo correlation and physiologically based pharmacokinetic platform was used to predict the pharmacokinetic profile of terfenadine following co-administration of the drug. Electrophysiological effects were simulated using the Cardiac Safety Simulator. The modulation of ion channel activity was dependent on the inhibitory potential of drugs on the main cardiac ion channels and a simulated free heart tissue concentration. ten Tusscher's human ventricular cardiomyocyte model was used to simulate the pseudo-ECG traces and further predict the pharmacodynamic consequences of DDI. Consistent with clinical observations, predicted plasma concentration profiles of terfenadine show considerable intra-subject variability with recorded Cmax values below 5 ng/mL for most virtual subjects. The pharmacokinetic and pharmacodynamic effects of inhibitors were predicted with reasonable accuracy. In all cases, a combination of the physiologically based pharmacokinetic and physiology-based pharmacodynamic models was able to differentiate between the terfenadine alone and terfenadine + inhibitor scenario. The range of QT prolongation was comparable in the clinical and virtual studies. The results indicate that mechanistic in vitro-in vivo correlation can be applied to predict the clinical effects of DDI even without comprehensive knowledge on all mechanisms contributing to the interaction.

Top-down, Bottom-up and Middle-out Strategies for Drug Cardiac Safety Assessment via Modeling and Simulations.

Cardiac safety is an issue causing early terminations at various stages of drug development. Efforts are put into the elimination of false negatives as well as false positives resulting from the current testing paradigm. In silico approaches offer mathematical system and data description from the ion current, through cardiomyocytes level, up to incorporation of inter-individual variability at the population level. The article aims to review three main modelling and simulation approaches, i.e. "top-down" which refers to models built on the observed data, "bottom-up", which stands for a mechanistic description of human physiology, and "middle-out" which combines both strategies. Modelling and simulation is a well-established tool in the assessment of drug proarrhythmic potency with an impact on research and development as well as on regulatory decisions, and it is certainly here to stay. What is more, the shift to systems biology and physiology-based models makes the cardiac effect more predictable.

In vitro-in vivo extrapolation of drug-induced proarrhythmia predictions at the population level.

Drug cardiotoxicity is a serious issue for patients, regulators, pharmaceutical companies and health service payers because they are all affected by its consequences. Despite the wide range of data they generate, existing approaches for cardiac safety testing might not be adequate and sufficiently cost-effective, probably as a result of the complexity of the problem. For this reason, translational tools (based on biophysically detailed, mathematical models) allowing for in vitro-in vivo extrapolation are gaining increasing interest. This current review describes approaches that can be used for cardiac safety assessment at the population level, by accounting for various sources of variability including kinetics of the compound of interest.

Predictive model for L-type channel inhibition: multichannel block in QT prolongation risk assessment.

Evaluation of the proarrhythmic potential of an investigated compound is now an integral element of the safety profile required for the approval of new drugs. The human ether-à-go-go-related gene (hERG) channel blocking potency is regarded as a surrogate marker of the proarrhythmic risk at the early stages of the research and development process. However, there is no straight correlation between QT prolongation and TdP occurrence probability, and hERG inhibition potential can be an inadequate predictor of QT prolongation. The L-type calcium channel plays a pivotal role in cardiomyocytes' physiology. Thus the main aim of this study was to develop a predictive model for drug-triggered CaL channel inhibition and also the assessment of drug-multichannel interaction effects on the heart rate-corrected QT interval. The data set, consisting of 123 records describing in vitro experimental settings, measured IC50 values and calculated physico-chemical properties for 72 various chemicals, was collected. The models were tested in a modified 10-fold cross-validation procedure. The generalization ability of the best model was as follows: root mean squared error (RMSE) = 1.10, normalized root mean squared error (NRMSE) = 16.09%. Out of the 10 most important variables, 5 described conditions of the in vitro experiments thus their description and experiment's conditions standardization might be the key to the models better performance. The simulations performed with the ToxComp system showed that the hERG block alone causes concentration-dependent QT prolongation, whereas when multichannel block is regarded, the effect could be reversed. For that reason, the multichannel interaction of tested compounds should be taken into consideration, in order to make the proarrhythmic risk assessment more reliable.

Virtual population generator for human cardiomyocytes parameters: in silico drug cardiotoxicity assessment.

BACKGROUND: The anatomical and histological parameters of the human ventricle depend on many factors including age and sex. Myocyte volume and electric capacitance are significant physiological parameters of left ventricle cardiomyocyte mathematical models. They allow the assessment of inter-individual variability during in vitro-in vivo extrapolation of the drug cardiotoxic effect. OBJECTIVE: The current research was carried out to analyze the relationship between age, human left ventricle cardiomyocyte volume, and electric capacitance in a healthy population. METHODS: In order to collect data describing cardiomyocyte volume and membrane area, literature searches were performed. It was assumed that the cardiomyocyte volume (VOL) and area (AREA) distribution have non-negative support and are skewed to the right. A log-linear model with constant variance was used. A simulation study was run to assess the influence of physiological parameters on action potential duration. RESULTS: The coefficient of determination for the proposed model R(2) = 0.95, that is, 95% of the variability observed in log cardiomyocyte volume can be explained by the estimated regression equation. To allow simple calculation and model performance validation, a simple Excel file was developed (Supplementary material). CONCLUSIONS: To the best of our knowledge, there is no other model available, combining age, cardiomyocyte volume, and area. The main limitations of the proposed models result from the assumptions made at the data analysis stage. The limited amount of information available in the literature and the lack of differentiation between sexes results in one common equation. The developed model is a part of the computational system for drug cardiotoxicity assessment.

Collation, assessment and analysis of literature in vitro data on hERG receptor blocking potency for subsequent modeling of drugs' cardiotoxic properties.

The assessment of the torsadogenic potency of a new chemical entity is a crucial issue during lead optimization and the drug development process. It is required by the regulatory agencies during the registration process. In recent years, there has been a considerable interest in developing in silico models, which allow prediction of drug-hERG channel interaction at the early stage of a drug development process. The main mechanism underlying an acquired QT syndrome and a potentially fatal arrhythmia called torsades de pointes is the inhibition of potassium channel encoded by hERG (the human ether-a-go-go-related gene). The concentration producing half-maximal block of the hERG potassium current (IC(50)) is a surrogate marker for proarrhythmic properties of compounds and is considered a test for cardiac safety of drugs or drug candidates. The IC(50) values, obtained from data collected during electrophysiological studies, are highly dependent on experimental conditions (i.e. model, temperature, voltage protocol). For the in silico models' quality and performance, the data quality and consistency is a crucial issue. Therefore the main objective of our work was to collect and assess the hERG IC(50) data available in accessible scientific literature to provide a high-quality data set for further studies.