A sensitive ultra-performance liquid chromatography-tandem mass spectrometry method for the simultaneous quantification of assay and trace-level genotoxic tosylate analogs (methyl and ethyl) in empagliflozin and its tablet dosage forms.

This study performed the simultaneous quantification of assay and two alkyl sulfonate (tosylate) analogs of empagliflozin (EGZ), specifically methyl 4-methyl benzene sulfonate (MMBS) and ethyl 4-methyl benzene sulfonate (EMBS) in EGZ, and its finished dosage form using an accurate and sensitive ultra-performance liquid chromatography-mass spectrometry method. The separation was achieved on a Waters Acquity BEH Shield RP18 (100 × 2.1 mm, 1.7 µm) column in gradient elution mode with 0.1% formic acid and acetonitrile as the mobile phases and a flow rate of 0.5 mL/min. For simultaneous quantification, the multiple reaction monitoring technique was utilized. The procedure was successfully validated in accordance with the International Council for Harmonisation of Technical Requirements for Pharmaceuticals for Human Use (ICH) guidelines. The peak areas of both impurities, along with their concentrations, exhibited a good relationship with Pearson's correlation coefficient (R), which was >0.999 in the range of 0.3-6 ppm with an EGZ concentration of 2 mg/mL. The percentage recoveries from the limit of quantitation (LOQ) to 200% to the specification level were in the range of 94.82%-102.92%, whereas the percentage relative standard deviation (%RSD) was <2. Therefore, this method is rapid and accurate to quantify MMBS, EMBS, and EGZ assay simultaneously from the marketed tablet dosage forms of EGZ for commercial release and stability sample testing.

Prediction of Adverse Effects of Drug-Drug Interactions on Cardiovascular System Based on the Analysis of Structure-Activity Relationships.

Co-administration of drugs often leads to drug-drug interactions, which could be accompanied by various adverse drug reactions that pose a threat to life and health of the patient. The effect caused by adverse drug reactions on cardiovascular system is one of the most significant manifestations of drug-drug interaction. Clinical assessment of adverse drug reactions resulting from drug-drug interaction between all drug pairs used in therapeutic practice is not possible. The purpose of this work was to build models using structure-activity analysis to predict adverse effects of drugs on cardiovascular system, mediated by pairwise interactions between the drug pairs when they are taken together. Data on the adverse effects resulting from drug-drug interaction were obtained from the DrugBank database. The data on drug pairs that do not cause such effects, which are necessary for building accurate structure-activity models, were obtained from the TwoSides database, which contains the results of analysis of the spontaneous reports. Two types of descriptors were used to describe a pair of drug structures: PoSMNA descriptors and probabilistic estimates of the prediction of biological activities obtained using the PASS program. Structure-activity relationships were established using the Random Forest method. Prediction accuracy was calculated by means of five-fold cross-validation. The highest accuracy values were obtained using PASS probabilistic estimates as descriptors. The area under the ROC curve was 0.94 for bradycardia, 0.96 for tachycardia, 0.90 for arrhythmia, 0.90 for ECG QT prolongation, 0.91 for hypertension, 0.89 for hypotension.

Assessing the utility of common arguments used in expert review of in silico predictions as part of ICH M7 assessments.

Expert review of two predictions, made by complementary (quantitative) structure-activity relationship models, to an overall conclusion is a key component of using in silico tools to assess the mutagenic potential of impurities as part of the ICH M7 guideline. In lieu of a specified protocol, numerous publications have presented best practise guides, often indicating the occurrence of common prediction scenarios and the evidence required to resolve them. A semi-automated expert review tool has been implemented in Lhasa Limited's Nexus platform following collation of these common arguments and assignment to the associated prediction scenarios made by Derek Nexus and Sarah Nexus. Using datasets primarily donated by pharmaceutical companies, an automated analysis of the frequency these prediction scenarios occur, and the likelihood of the associated arguments assigning the correct resolution, could then be conducted. This article highlights that a relatively small number of common arguments may be used to accurately resolve many prediction scenarios to a single conclusion. The use of a standardised method of argumentation and assessment of evidence for a given impurity is proposed to improve the efficiency and consistency of expert review as part of an ICH M7 submission.

Migration of styrene oligomers from food contact materials: in silico prediction of possible genotoxicity.

Styrene oligomers (SO) are well-known side products formed during styrene polymerization. They consist mainly of dimers (SD) and trimers (ST) that have been shown to be still residual in polystyrene (PS) materials. In this study migration of SO from PS into sunflower oil at temperatures between 5 and 70 °C and contact times between 0.5 h and 10 days was investigated. In addition, the contents of SD and ST in the fatty foodstuffs créme fraiche and coffee cream, which are typically enwrapped in PS, were measured and the amounts detected (of up to 0.123 mg/kg food) were compared to literature data. From this comparison, it became evident, that the levels of SO migrating from PS packaging into real food call for a comprehensive risk assessment. As a first step towards this direction, possible genotoxicity has to be addressed. Due to technical and experimental limitations, however, the few existing in vitro tests available are unsuited to provide a clear picture. In order to reduce uncertainty of these in vitro tests, four different knowledge and statistics-based in silico tools were applied to such SO that are known to migrate into food. Except for SD4 all evaluated SD and ST showed no alert for genotoxicity. For SD4, either the predictions were inconclusive or the substance was assigned as being out of the chemical space (out of domain) of the respective in silico tool. Therefore, the absence of genotoxicity of SD4 requires additional experimental proof. Apart from SD4, in silico studies supported the limited in vitro data that indicated the absence of genotoxicity of SO. In conclusion, the overall migration of all SO together into food of up to 50 µg/kg does not raise any health concerns, given the currently available in silico and in vitro data.

Evaluation of the predictivity of Acute Oral Toxicity (AOT) structure-activity relationship models.

Several public efforts are aimed at discovering patterns or classifiers in the high-dimensional bioactivity space that predict tissue, organ or whole animal toxicological endpoints. The current study sought to assess and compare the predictions of the Globally Harmonized System (GHS) categories and Dangerous Goods (DG) classifications based on Lethal Dose (LD50) from several available tools (ACD/Labs, Leadscope, T.E.S.T., CATMoS, CaseUltra). External validation was done using dataset of 375 substances to demonstrate their predictive capacity. All models showed very good performance for identifying non-toxic compounds, which would be useful for DG classification, developing or triaging new chemicals, prioritizing existing chemicals for more detailed and rigorous toxicity assessments, and assessing non-active pharmaceutical intermediates. This would ultimately reduce animal use and improve risk assessments. Category-to-category prediction was not optimal, mainly due to the tendency to overpredict the outcome and the general limitations of acute oral toxicity (AOT) in vivo studies. Overprediction does not specifically pose a risk to human health, it can impact transport and material packaging requirements. Performance for compounds with LD50 ≤ 300 mg/kg (approx. 5% of the dataset) was the poorest among all groups and could be potentially improved by including expert review and read-across to similar substances.

Analytical Method Development for 19 Alkyl Halides as Potential Genotoxic Impurities by Analytical Quality by Design.

Major issues in the pharmaceutical industry involve efficient risk management and control strategies of potential genotoxic impurities (PGIs). As a result, the development of an appropriate method to control these impurities is required. An optimally sensitive and simultaneous analytical method using gas chromatography with a mass spectrometry detector (GC-MS) was developed for 19 alkyl halides determined to be PGIs. These 19 alkyl halides were selected from 144 alkyl halides through an in silico study utilizing quantitative structure-activity relationship (Q-SAR) approaches via expert knowledge rule-based software and statistical-based software. The analytical quality by design (QbD) approach was adopted for the development of a sensitive and robust analytical method for PGIs. A limited number of literature studies have reviewed the analytical QbD approach in the PGI method development using GC-MS as the analytical instrument. A GC equipped with a single quadrupole mass spectrometry detector (MSD) and VF-624 ms capillary column was used. The developed method was validated in terms of specificity, the limit of detection, quantitation, linearity, accuracy, and precision, according to the ICH Q2 guideline.

Computational Prediction of Inhibitors and Inducers of the Major Isoforms of Cytochrome P450.

Human cytochrome P450 enzymes (CYPs) are heme-containing monooxygenases. This superfamily of drug-metabolizing enzymes is responsible for the metabolism of most drugs and other xenobiotics. The inhibition of CYPs may lead to drug-drug interactions and impair the biotransformation of drugs. CYP inducers may decrease the bioavailability and increase the clearance of drugs. Based on the freely available databases ChEMBL and PubChem, we have collected over 70,000 records containing the structures of inhibitors and inducers together with the IC50 values for the inhibitors of the five major human CYPs: 1A2, 3A4, 2D6, 2C9, and 2C19. Based on the collected data, we developed (Q)SAR models for predicting inhibitors and inducers of these CYPs using GUSAR and PASS software. The developed (Q)SAR models could be applied for assessment of the interaction of novel drug-like substances with the major human CYPs. The created (Q)SAR models demonstrated reasonable accuracy of prediction. They have been implemented in the web application P450-Analyzer that is freely available via the Internet.

Assessing the impact of expert knowledge on ICH M7 (Q)SAR predictions. Is expert review still needed?

The ICH M7 (R1) guideline recommends the use of complementary (Q)SAR models to assess the mutagenic potential of drug impurities as a state-of-the-art, high-throughput alternative to empirical testing. Additionally, it includes a provision for the application of expert knowledge to increase prediction confidence and resolve conflicting calls. Expert knowledge, which considers structural analogs and mechanisms of activity, has been valuable when models return an indeterminate (equivocal) result or no prediction (out-of-domain). A retrospective analysis of 1002 impurities evaluated in drug regulatory applications between April 2017 and March 2019 assessed the impact of expert review on (Q)SAR predictions. Expert knowledge overturned the default predictions for 26% of the impurities and resolved 91% of equivocal predictions and 75% of out-of-domain calls. Of the 261 overturned default predictions, 15% were upgraded to equivocal or positive and 79% were downgraded to equivocal or negative. Chemical classes with the most overturns were primary aromatic amines (46%), aldehydes (45%), Michael-reactive acceptors (37%), and non-primary alkyl halides (33%). Additionally, low confidence predictions were the most often overturned. Collectively, the results suggest that expert knowledge continues to play an important role in an ICH M7 (Q)SAR prediction workflow and triaging predictions based on chemical class and probability can improve (Q)SAR review efficiency.

A cross-industry collaboration to assess if acute oral toxicity (Q)SAR models are fit-for-purpose for GHS classification and labelling.

This study assesses whether currently available acute oral toxicity (AOT) in silico models, provided by the widely employed Leadscope software, are fit-for-purpose for categorization and labelling of chemicals. As part of this study, a large data set of proprietary and marketed compounds from multiple companies (pharmaceutical, plant protection products, and other chemical industries) was assembled to assess the models' performance. The absolute percentage of correct or more conservative predictions, based on a comparison of experimental and predicted GHS categories, was approximately 95%, after excluding a small percentage of inconclusive (indeterminate or out of domain) predictions. Since the frequency distribution across the experimental categories is skewed towards low toxicity chemicals, a balanced assessment was also performed. Across all compounds which could be assigned to a well-defined experimental category, the average percentage of correct or more conservative predictions was around 80%. These results indicate the potential for reliable and broad application of these models across different industrial sectors. This manuscript describes the evaluation of these models, highlights the importance of an expert review, and provides guidance on the use of AOT models to fulfill testing requirements, GHS classification/labelling, and transportation needs.

Degradation of the anticancer drug flutamide by solar photo-Fenton treatment at near-neutral pH: Identification of transformation products and in silico (Q)SAR risk assessment.

Flutamide (FLUT) is a non-steroidal drug mainly used in the treatment of prostate cancer and has been detected in the aquatic environment at ng L-1 levels. The environmental fate and effects of FLUT have not yet been studied. Conventional treatment technologies fail to completely remove pharmaceuticals, so the solar photo-Fenton process (SPF) has been proposed as an alternative. In this study, the degradation of FLUT, at two different initial concentrations in ultra-pure water, was carried out by SPF. The initial SPF conditions were pH0 5, [Fe2+]0 = 5 mg L-1, and [H2O2]0 = 50 mg L-1. Preliminary elimination rates of 53.4% and 73.4%. The kinetics of FLUT degradation could be fitted by a pseudo-first order model and the kobs were 6.57 × 10-3 and 9.13 × 10-3 min-1 t30W and the half-life times were 95.62 and 73.10 min t30W were achieved for [FLUT]0 of 5 mg L-1 and 500 µg L-1, respectively. Analysis using LC-QTOF MS identified thirteen transformation products (TPs) during the FLUT degradation process. The main degradation pathways proposed were hydroxylation, hydrogen abstraction, demethylation, NO2 elimination, cleavage, and aromatic ring opening. Different in silico (quantitative) structure-activity relationship ((Q)SAR) freeware models were used to predict the toxicities and environmental fates of FLUT and the TPs. The in silico predictions indicated that these substances were not biodegradable, while some TPs were classified near the threshold point to be considered as PBT compounds. The in silico (Q)SAR predictions gave positive alerts concerning the mutagenicity and carcinogenicity endpoints. Additionally, the (Q)SAR toolbox software provided structural alerts corresponding to the positive alerts obtained with the different mutagenicity and carcinogenicity models, supporting the positive alerts with more proactive information.

Development of improved QSAR models for predicting the outcome of the in vivo micronucleus genetic toxicity assay.

All drugs entering clinical trials are expected to undergo a series of in vitro and in vivo genotoxicity tests as outlined in the International Council on Harmonization (ICH) S2 (R1) guidance. Among the standard battery of genotoxicity tests used for pharmaceuticals, the in vivo micronucleus assay, which measures the frequency of micronucleated cells mostly from blood or bone marrow, is recommended for detecting clastogens and aneugens. (Quantitative) structure-activity relationship [(Q)SAR] models may be used as early screening tools by pharmaceutical companies to assess genetic toxicity risk during drug candidate selection. Models can also provide decision support information during regulatory review as part of the weight-of-evidence when experimental data are insufficient. In the present study, two commercial (Q)SAR platforms were used to construct in vivo micronucleus models from a recently enhanced in-house database of non-proprietary study findings in mice. Cross-validated performance statistics for the new models showed sensitivity of up to 74% and negative predictivity of up to 86%. In addition, the models demonstrated cross-validated specificity of up to 77% and coverage of up to 94%. These new models will provide more reliable predictions and offer an investigational approach for drug safety assessment with regards to identifying potentially genotoxic compounds.

Evaluation of the applicability of existing (Q)SAR models for predicting the genotoxicity of pesticides and similarity analysis related with genotoxicity of pesticides for facilitating of grouping and read across: An EFSA funded project.

To facilitate the practical implementation of the guidance on the residue definition for dietary risk assessment, EFSA has organized an evaluation of applicability of existing in silico models for predicting the genotoxicity of pesticides and their metabolites, including literature survey, application of QSARs and development of Read Across methodologies. This paper summarizes the main results. For the Ames test, all (Q)SAR models generated statistically significant predictions, comparable with the experimental variability of the test. The reliability of the models for other assays/endpoints appears to be still far from optimality. Two new Read Across approaches were evaluated: Read Across was largely successful for predicting the Ames test results, but less for in vitro Chromosomal Aberrations. The worse results for non-Ames endpoints may be attributable to the several revisions of experimental protocols and evaluation criteria of results, that have made the databases qualitatively non-homogeneous and poorly suitable for modeling. Last, Parent/Metabolite structural differences (besides known Structural Alerts) that may, or may not cause changes in the Ames mutagenicity were identified and catalogued. The findings from this work are suitable for being integrated into Weight-of-Evidence and Tiered evaluation schemes. Areas needing further developments are pointed out.

Management of pharmaceutical ICH M7 (Q)SAR predictions - The impact of model updates.

Pharmaceutical applicants conduct (Q)SAR assessments on identified and theoretical impurities to predict their mutagenic potential. Two complementary models-one rule-based and one statistical-based-are used, followed by expert review. (Q)SAR models are continuously updated to improve predictions, with new versions typically released on a yearly basis. Numerous releases of (Q)SAR models will occur during the typical 6-7 years of drug development until new drug registration. Therefore, it is important to understand the impact of model updates on impurity mutagenicity predictions over time. Compounds representative of pharmaceutical impurities were analyzed with three rule- and three statistical-based models covering a 4-8 year period, with the individual time frame being dependent on when the individual models were initially made available. The largest changes in the combined outcome of two complementary models were from positive or equivocal to negative and from negative to equivocal. Importantly, the cumulative change of negative to positive predictions was small in all models (<5%) and was further reduced when complementary models were combined in a consensus fashion. We conclude that model updates of the type evaluated in this manuscript would not necessarily require re-running a (Q)SAR prediction unless there is a specific need. However, original (Q)SAR predictions should be evaluated when finalizing the commercial route of synthesis for marketing authorization.

Skin sensitization in silico protocol.

The assessment of skin sensitization has evolved over the past few years to include in vitro assessments of key events along the adverse outcome pathway and opportunistically capitalize on the strengths of in silico methods to support a weight of evidence assessment without conducting a test in animals. While in silico methods vary greatly in their purpose and format; there is a need to standardize the underlying principles on which such models are developed and to make transparent the implications for the uncertainty in the overall assessment. In this contribution, the relationship between skin sensitization relevant effects, mechanisms, and endpoints are built into a hazard assessment framework. Based on the relevance of the mechanisms and effects as well as the strengths and limitations of the experimental systems used to identify them, rules and principles are defined for deriving skin sensitization in silico assessments. Further, the assignments of reliability and confidence scores that reflect the overall strength of the assessment are discussed. This skin sensitization protocol supports the implementation and acceptance of in silico approaches for the prediction of skin sensitization.

Genetic toxicology in silico protocol.

In silico toxicology (IST) approaches to rapidly assess chemical hazard, and usage of such methods is increasing in all applications but especially for regulatory submissions, such as for assessing chemicals under REACH as well as the ICH M7 guideline for drug impurities. There are a number of obstacles to performing an IST assessment, including uncertainty in how such an assessment and associated expert review should be performed or what is fit for purpose, as well as a lack of confidence that the results will be accepted by colleagues, collaborators and regulatory authorities. To address this, a project to develop a series of IST protocols for different hazard endpoints has been initiated and this paper describes the genetic toxicity in silico (GIST) protocol. The protocol outlines a hazard assessment framework including key effects/mechanisms and their relationships to endpoints such as gene mutation and clastogenicity. IST models and data are reviewed that support the assessment of these effects/mechanisms along with defined approaches for combining the information and evaluating the confidence in the assessment. This protocol has been developed through a consortium of toxicologists, computational scientists, and regulatory scientists across several industries to support the implementation and acceptance of in silico approaches.

A defined approach for predicting skin sensitisation hazard and potency based on the guided integration of in silico, in chemico and in vitro data using exclusion criteria.

A decision tree-based defined approach (DA) has been designed using exclusion criteria based on applicability domain knowledge of in chemico/in vitro information sources covering key events 1-3 in the skin sensitisation adverse outcome pathway and an in silico tool predicting the adverse outcome (Derek Nexus). The hypothesis is that using exclusion criteria to de-prioritise less applicable assays and/or in silico outcomes produces a rational, transparent, and reliable DA for the prediction of skin sensitisation potential. Five exclusion criteria have been established: Derek Nexus reasoning level, Derek Nexus negative prediction, metabolism, lipophilicity, and lysine-reactivity. These are used to prioritise the most suitable information sources for a given chemical and results from which are used in a '2 out of 3' approach to provide a prediction of hazard. A potency category (and corresponding GHS classification) is then assigned using a k-Nearest Neighbours model containing human and LLNA data. The DA correctly identified the hazard (sensitiser/non-sensitiser) for 85% and 86% of a dataset with reference LLNA and human data. The correct potency category was identified for 59% and 68% of chemicals, and the GHS classification accurately predicted for 73% and 76% with reference LLNA and human data, respectively.

(Q)SAR Models of HIV-1 Protein Inhibition by Drug-Like Compounds.

Despite the achievements of antiretroviral therapy, discovery of new anti-HIV medicines remains an essential task because the existing drugs do not provide a complete cure for the infected patients, exhibit severe adverse effects, and lead to the appearance of resistant strains. To predict the interaction of drug-like compounds with multiple targets for HIV treatment, ligand-based drug design approach is widely applied. In this study, we evaluated the possibilities and limitations of (Q)SAR analysis aimed at the discovery of novel antiretroviral agents inhibiting the vital HIV enzymes. Local (Q)SAR models are based on the analysis of structure-activity relationships for molecules from the same chemical class, which significantly restrict their applicability domain. In contrast, global (Q)SAR models exploit data from heterogeneous sets of drug-like compounds, which allows their application to databases containing diverse structures. We compared the information for HIV-1 integrase, protease and reverse transcriptase inhibitors available in the EBI ChEMBL, NIAID HIV/OI/TB Therapeutics, and Clarivate Analytics Integrity databases as the sources for (Q)SAR training sets. Using the PASS and GUSAR software, we developed and validated a variety of (Q)SAR models, which can be further used for virtual screening of new antiretrovirals in the SAVI library. The developed models are implemented in the freely available web resource AntiHIV-Pred.

Making reliable negative predictions of human skin sensitisation using an in silico fragmentation approach.

A previously published fragmentation method for making reliable negative in silico predictions has been applied to the problem of predicting skin sensitisation in humans, making use of a dataset of over 2750 chemicals with publicly available skin sensitisation data from 18 in vivo assays. An assay hierarchy was designed to enable the classification of chemicals within this dataset as either sensitisers or non-sensitisers where data from more than one in vivo test was available. The negative prediction approach was validated internally, using a 5-fold cross-validation, and externally, against a proprietary dataset of approximately 1000 chemicals with in vivo reference data shared by members of the pharmaceutical, nutritional, and personal care industries. The negative predictivity for this proprietary dataset was high in all cases (>75%), and the model was also able to identify structural features that resulted in a lower accuracy or a higher uncertainty in the negative prediction, termed misclassified and unclassified features respectively. These features could serve as an aid for further expert assessment of the negative in silico prediction.

Construction and application of (Q)SAR models to predict chemical-induced in vitro chromosome aberrations.

In drug development, genetic toxicology studies are conducted using in vitro and in vivo assays to identify potential mutagenic and clastogenic effects, as outlined in the International Council for Harmonisation (ICH) S2 regulatory guideline. (Quantitative) structure-activity relationship ((Q)SAR) models that predict assay outcomes can be used as an early screen to prioritize pharmaceutical candidates, or later during product development to evaluate safety when experimental data are unavailable or inconclusive. In the current study, two commercial QSAR platforms were used to build models for in vitro chromosomal aberrations in Chinese hamster lung (CHL) and Chinese hamster ovary (CHO) cells. Cross-validated CHL model predictive performance showed sensitivity of 80 and 82%, and negative predictivity of 75 and 76% based on 875 training set compounds. For CHO, sensitivity of 61 and 67% and negative predictivity of 68 and 74% was achieved based on 817 training set compounds. The predictive performance of structural alerts in a commercial expert rule-based SAR software was also investigated and showed positive predictivity of 48-100% for selected alerts. Case studies examining incorrectly-predicted compounds, non-DNA-reactive clastogens, and recently-approved pharmaceuticals are presented, exploring how an investigational approach using similarity searching and expert knowledge can improve upon individual (Q)SAR predictions of the clastogenicity of drugs.

Integrating computational methods to predict mutagenicity of aromatic azo compounds.

Azo dyes have several industrial uses. However, these azo dyes and their degradation products showed mutagenicity, inducing damage in environmental and human systems. Computational methods are proposed as cheap and rapid alternatives to predict the toxicity of azo dyes. A benchmark dataset of Ames data for 354 azo dyes was employed to develop three classification strategies using knowledge-based methods and docking simulations. Results were compared and integrated with three models from the literature, developing a series of consensus strategies. The good results confirm the usefulness of in silico methods as a support for experimental methods to predict the mutagenicity of azo compounds.