Quantum-Mechanics Calculations Elucidate Skin-Sensitizing Pharmaceutical Compounds.

Skin sensitization is a critical end point in occupational toxicology that necessitates the use of fast, accurate, and affordable models to aid in establishing handling guidance for worker protection. While many in silico models have been developed, the scarcity of reliable data for active pharmaceutical ingredients (APIs) and their intermediates (together regarded as pharmaceutical compounds) brings into question the reliability of these tools, which are largely constructed using publicly available nonspecialty chemicals. Here, we present the quantum-mechanical (QM) Computer-Aided Discovery and REdesign (CADRE) model, which was developed with the bioactive and structurally complex chemical space in mind by relying on the fundamentals of chemical interactions in key events (versus structural attributes of training-set data). Validated in this study on 345 APIs and intermediates, CADRE achieved 95% accuracy, sensitivity, and specificity and a combined 79% accuracy in assigning potency categories compared to the mouse local lymph node assay data. We show how historical outcomes from CADRE testing in the pharmaceutical space, generated over the past 10 years on ca. 2500 chemicals, can be used to probe the relationships between sensitization mechanisms (or the underlying chemical classes) and the probability of eliciting a sensitization response in mice of a given potency. We believe this information to be of value to both practitioners, who can use it to quickly screen and triage their data sets, as well as to model developers to fine-tune their structure-based tools. Lastly, we leverage our experimentally validated subset of APIs and intermediates to show the importance of dermal permeability on the sensitization potential and potency. We demonstrate that common physicochemical properties used to assess permeation, such as the octanol-water partition coefficient and molecular weight, are poor proxies for the more accurate energy-pair distributions that can be computed from mixed QM and classical simulations using model representations of the stratum corneum.

Mechanisms of Nitrosamine Mutagenicity and Their Relationship to Rodent Carcinogenic Potency.

A thorough literature review was undertaken to understand how the pathways of N-nitrosamine transformation relate to mutagenic potential and carcinogenic potency in rodents. Empirical and computational evidence indicates that a common radical intermediate is created by CYP-mediated hydrogen abstraction at the α-carbon; it is responsible for both activation, leading to the formation of DNA-reactive diazonium species, and deactivation by denitrosation. There are competing sites of CYP metabolism (e.g., ß-carbon), and other reactive species can form following initial bioactivation, although these alternative pathways tend to decrease rather than enhance carcinogenic potency. The activation pathway, oxidative dealkylation, is a common reaction in drug metabolism and evidence indicates that the carbonyl byproduct, e.g., formaldehyde, does not contribute to the toxic properties of N-nitrosamines. Nitric oxide (NO), a side product of denitrosation, can similarly be discounted as an enhancer of N-nitrosamine toxicity based on carcinogenicity data for substances that act as NO-donors. However, not all N-nitrosamines are potent rodent carcinogens. In a significant number of cases, there is a potency overlap with non-N-nitrosamine carcinogens that are not in the Cohort of Concern (CoC; high-potency rodent carcinogens comprising aflatoxin-like-, N-nitroso-, and alkyl-azoxy compounds), while other N-nitrosamines are devoid of carcinogenic potential. In this context, mutagenicity is a useful surrogate for carcinogenicity, as proposed in the ICH M7 (R2) (2023) guidance. Thus, in the safety assessment and control of N-nitrosamines in medicines, it is important to understand those complementary attributes of mechanisms of mutagenicity and structure-activity relationships that translate to elevated potency versus those which are associated with a reduction in, or absence of, carcinogenic potency.

N-Nitrosamine drug substance related impurities (NDSRIs) - A proposal for the addition of subcategories to carcinogenic potency categorization approach categories 1 and 2 for NDSRIs with a molecular weight > 200 Da.

The carcinogenicity potency categorization approach (CPCA) derived and harmonized by Health Authorities was a significant milestone, as it defined molecular properties that allowed for the rapid evaluation of the chemical structures of N-nitrosamine drug substance related impurities (NDSRIs) and the assignment of associated lifetime Acceptable Intake (AI) limits to inform on appropriate impurity control strategies in certain drug products. Nonetheless, it is important to continue to refine and improve on the CPCA based upon data-derived evidence. Herein, we focus on the default CPCA AI for NDSRIs, which is largely based on the small molecule N-nitrosamines (NAs). Considering the carcinogenic potency of NAs with a molecular weight >200 Da (NDSRIs molecular weight is typically 200-600 Da), we propose that in the absence of any compound specific data, the lowest lifetime Acceptable Intake for NAs, such as NDSRIs, should be 10x less (i.e., 150 ng/day) than the ICH M7 Threshold of Toxicological Concern of 1500 ng/day, (even for NDSRIs that are considered CPCA Category 1 and 2) which would conservatively result in a theoretical cancer risk of <1 in 100,000.

Acceptable intakes (AIs) for 11 small molecule N-nitrosamines (NAs).

Low levels of N-nitrosamines (NAs) were detected in pharmaceuticals and, as a result, health authorities (HAs) have published acceptable intakes (AIs) in pharmaceuticals to limit potential carcinogenic risk. The rationales behind the AIs have not been provided to understand the process for selecting a TD50 or read-across analog. In this manuscript we evaluated the toxicity data for eleven common NAs in a comprehensive and transparent process consistent with ICH M7. This evaluation included substances which had datasets that were robust, limited but sufficient, and substances with insufficient experimental animal carcinogenicity data. In the case of robust or limited but sufficient carcinogenicity information, AIs were calculated based on published or derived TD50s from the most sensitive organ site. In the case of insufficient carcinogenicity information, available carcinogenicity data and structure activity relationships (SARs) were applied to categorical-based AIs of 1500 ng/day, 150 ng/day or 18 ng/day; however additional data (such as biological or additional computational modelling) could inform an alternative AI. This approach advances the methodology used to derive AIs for NAs.

An Evaluation of the Occupational Health Hazards of Peptide Couplers.

Peptide couplers (also known as amide bond-forming reagents or coupling reagents) are broadly used in organic chemical syntheses, especially in the pharmaceutical industry. Yet, occupational health hazards associated with this chemical class are largely unexplored, which is disconcerting given the intrinsic reactivity of these compounds. Several case studies involving occupational exposures reported adverse respiratory and dermal health effects, providing initial evidence of chemical sensitization. To address the paucity of toxicological data, a pharmaceutical cross-industry task force was formed to evaluate and assess the potential of these compounds to cause eye and dermal irritation as well as corrosivity and dermal sensitization. The goal of our work was to inform health and safety professionals as well as pharmaceutical and organic chemists of the occupational health hazards associated with this chemical class. To that end, 25 of the most commonly used peptide couplers and five hydrolysis products were selected for in vivo, in vitro, and in silico testing. Our findings confirmed that dermal sensitization is a concern for this chemical class with 21/25 peptide couplers testing positive for dermal sensitization and 15 of these being strong/extreme sensitizers. We also found that dermal corrosion and irritation (8/25) as well as eye irritation (9/25) were health hazards associated with peptide couplers and their hydrolysis products (4/5 were dermal irritants or corrosive and 4/5 were eye irritants). Resulting outcomes were synthesized to inform decision making in peptide coupler selection and enable data-driven hazard communication to workers. The latter includes harmonized hazard classifications, appropriate handling recommendations, and accurate safety data sheets, which support the industrial hygiene hierarchy of control strategies and risk assessment. Our study demonstrates the merits of an integrated, in vivo -in silico analysis, applied here to the skin sensitization endpoint using the Computer-Aided Discovery and REdesign (CADRE) and Derek Nexus programs. We show that experimental data can improve predictive models by filling existing data gaps while, concurrently, providing computational insights into key initiating events and elucidating the chemical structural features contributing to adverse health effects. This interactive, interdisciplinary approach is consistent with Green Chemistry principles that seek to improve the selection and design of less hazardous reagents in industrial processes and applications.

The utility of hERG channel inhibition data in the derivation of occupational exposure limits.

Inhibition of the human ether-à-go-go (hERG) channel may lead to QT prolongation and fatal arrhythmia. While pharmaceutical drug candidates that exhibit potent hERG channel inhibition often fail early in development, many drugs with both cardiac and non-cardiac indications proceed to market. In this study, the relationship between in vitro hERG channel inhibition and published occupational exposure limit (OEL) was evaluated. A total of 23 cardiac drugs and 44 drugs with non-cardiac indications with published hERG channel IC50 and published OELs were identified. There was an apparent relationship between hERG IC50 potency and the OEL for cardiac and non-cardiac drugs. Twenty cardiac and non-cardiac drugs were identified that had a potent hERG IC50 (≤25 µM) and a contrastingly large OEL value (≥100 µg/m3). OELs or hazard banding corresponding to ≤100 µg/m3 should be sufficiently protective of effects following occupational exposure to the majority of APIs with hERG IC50 values ≤ 100 µM. It is important to consider hERG IC50 values and possible cardiac effects when deriving OEL values for drugs, regardless of indication. These considerations may be particularly important early in the drug development process for establishing exposure control bands for drugs that do not yet have full clinical safety data.

Use of the bacterial reverse mutation assay to predict carcinogenicity of N-nitrosamines.

Under ICH M7, impurities are assessed using the bacterial reverse mutation assay (i.e., Ames test) when predicted positive using in silico methodologies followed by expert review. N-Nitrosamines (NAs) have been of recent concern as impurities in pharmaceuticals, mainly because of their potential to be highly potent mutagenic carcinogens in rodent bioassays. The purpose of this analysis was to determine the sensitivity of the Ames assay to predict the carcinogenic outcome with curated proprietary Vitic (n = 131) and Leadscope (n = 70) databases. NAs were selected if they had corresponding rodent carcinogenicity assays. Overall, the sensitivity/specificity of the Ames assay was 93-97% and 55-86%, respectively. The sensitivity of the Ames assay was not significantly impacted by plate incorporation (84-89%) versus preincubation (82-89%). Sensitivity was not significantly different between use of rat and hamster liver induced S9 (80-93% versus 77-96%). The sensitivity of the Ames is high when using DMSO as a solvent (87-88%). Based on the analysis of these databases, the Ames assay conducted under OECD 471 guidelines is highly sensitive for detecting the carcinogenic hazards of NAs.

Calculating qualified non-mutagenic impurity levels: Harmonization of approaches.

The presence of impurities in drugs is unavoidable. As impurities offer no direct benefit to the patient, it is critical that impurities do not compromise patient safety. Current guidelines on the derivation of acceptable impurity levels leave aspects of calculations open for interpretation, resulting in inconsistencies across industry and regulators. To understand current impurity qualification practices from a safety standpoint, regulatory expectations and the safety risk that impurities pose, the IQ DruSafe Impurities Working Group (WG) conducted a pharmaceutical industry-wide survey. Survey results highlighted areas that could benefit from harmonization, including nonclinical species/sex selection and the application of adjustment factors (i.e., body surface area). Recommendations for alignment on these topics is included in this publication. Additionally, the WG collated repeat-dose toxicity information for 181 starting materials and intermediates, reflective of pharmaceutical impurities, to understand the toxicological risks they generally pose in relation to the drug substance (DS) and the assumptions surrounding the calculation of qualified impurity levels. An evaluation of this dataset and the survey were used to harmonize how to calculate a safe limit for an impurity based on toxicology testing of the impurity when present within the DS.

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.

Extending (Q)SARs to incorporate proprietary knowledge for regulatory purposes: is aromatic N-oxide a structural alert for predicting DNA-reactive mutagenicity?

(Quantitative) structure-activity relationship or (Q)SAR predictions of DNA-reactive mutagenicity are important to support both the design of new chemicals and the assessment of impurities, degradants, metabolites, extractables and leachables, as well as existing chemicals. Aromatic N-oxides represent a class of compounds that are often considered alerting for mutagenicity yet the scientific rationale of this structural alert is not clear and has been questioned. Because aromatic N-oxide-containing compounds may be encountered as impurities, degradants and metabolites, it is important to accurately predict mutagenicity of this chemical class. This article analysed a series of publicly available aromatic N-oxide data in search of supporting information. The article also used a previously developed structure-activity relationship (SAR) fingerprint methodology where a series of aromatic N-oxide substructures was generated and matched against public and proprietary databases, including pharmaceutical data. An assessment of the number of mutagenic and non-mutagenic compounds matching each substructure across all sources was used to understand whether the general class or any specific subclasses appear to lead to mutagenicity. This analysis resulted in a downgrade of the general aromatic N-oxide alert. However, it was determined there were enough public and proprietary data to assign the quindioxin and related chemicals as well as benzo[c][1,2,5]oxadiazole 1-oxide subclasses as alerts. The overall results of this analysis were incorporated into Leadscope's expert-rule-based model to enhance its predictive accuracy.

Discovery of velpatasvir (GS-5816): A potent pan-genotypic HCV NS5A inhibitor in the single-tablet regimens Vosevi® and Epclusa®.

Direct-acting antiviral inhibitors have revolutionized the treatment of hepatitis C virus (HCV) infected patients. Herein is described the discovery of velpatasvir (VEL, GS-5816), a potent pan-genotypic HCV NS5A inhibitor that is a component of the only approved pan-genotypic single-tablet regimens (STRs) for the cure of HCV infection. VEL combined with sofosbuvir (SOF) is Epclusa®, an STR with 98% cure-rates for genotype 1-6 HCV infected patients. Addition of the pan-genotypic HCV NS3/4A protease inhibitor voxilaprevir to SOF/VEL is the STR Vosevi®, which affords 97% cure-rates for genotype 1-6 HCV patients who have previously failed another treatment regimen.

Principles and procedures for handling out-of-domain and indeterminate results as part of ICH M7 recommended (Q)SAR analyses.

The International Council for Harmonization (ICH) M7 guideline describes a hazard assessment process for impurities that have the potential to be present in a drug substance or drug product. In the absence of adequate experimental bacterial mutagenicity data, (Q)SAR analysis may be used as a test to predict impurities' DNA reactive (mutagenic) potential. However, in certain situations, (Q)SAR software is unable to generate a positive or negative prediction either because of conflicting information or because the impurity is outside the applicability domain of the model. Such results present challenges in generating an overall mutagenicity prediction and highlight the importance of performing a thorough expert review. The following paper reviews pharmaceutical and regulatory experiences handling such situations. The paper also presents an analysis of proprietary data to help understand the likelihood of misclassifying a mutagenic impurity as non-mutagenic based on different combinations of (Q)SAR results. This information may be taken into consideration when supporting the (Q)SAR results with an expert review, especially when out-of-domain results are generated during a (Q)SAR evaluation.

In Silico Approaches In Carcinogenicity Hazard Assessment: Current Status and Future Needs.

Historically, identifying carcinogens has relied primarily on tumor studies in rodents, which require enormous resources in both money and time. In silico models have been developed for predicting rodent carcinogens but have not yet found general regulatory acceptance, in part due to the lack of a generally accepted protocol for performing such an assessment as well as limitations in predictive performance and scope. There remains a need for additional, improved in silico carcinogenicity models, especially ones that are more human-relevant, for use in research and regulatory decision-making. As part of an international effort to develop in silico toxicological protocols, a consortium of toxicologists, computational scientists, and regulatory scientists across several industries and governmental agencies evaluated the extent to which in silico models exist for each of the recently defined 10 key characteristics (KCs) of carcinogens. This position paper summarizes the current status of in silico tools for the assessment of each KC and identifies the data gaps that need to be addressed before a comprehensive in silico carcinogenicity protocol can be developed for regulatory use.

Spiropiperidine CCR5 antagonists.

A novel series of CCR5 antagonists has been identified, utilizing leads from high-throughput screening which were further modified based on insights from competitor molecules. Lead optimization was pursued by balancing opposing trends of metabolic stability and potency. Selective and potent analogs with good pharmacokinetic properties were successfully developed.

In Vivo Mutagenicity Testing of Arylboronic Acids and Esters.

Arylboronic acids and esters (referred to collectively as arylboronic compounds) are commonly used intermediates in the synthesis of pharmaceuticals but pose a challenge for chemical syntheses because they are often positive for bacterial mutagenicity in vitro. As such, arylboronic compounds are then typically controlled to levels that are acceptable for mutagenic impurities, that is, the threshold of toxicological concern (TTC). This study used ICH M7 guidance to design and conduct a testing strategy to investigate the in vivo relevance of the in vitro positive findings of arylboronic compounds. Eight arylboronic compounds representing a variety of chemical scaffolds were tested in Sprague Dawley and/or Wistar rats in the in vivo Pig-a (peripheral blood reticulocytes and mature red blood cells) and/or comet assays (duodenum and/or liver). Five of the eight compounds were also tested in the micronucleus (peripheral blood) assay. The arylboronic compounds tested orally demonstrated high systemic exposure; thus the blood and bone marrow were adequately exposed to test article. One compound was administered intravenously due to formulation stability issues. This investigation showed that arylboronic compounds that were mutagenic in vitro were not found to be mutagenic in the corresponding in vivo assays. Therefore, arylboronic compounds similar to the scaffolds tested in this article may be considered non-mutagenic and managed in accordance with the ICH Q3A/Q3B guidelines. Environ. Mol. Mutagen. 2019. © 2019 Wiley Periodicals, Inc.