Abstract
Recent withdrawal of several drugs from the market due to elevated levels of N-nitrosamine impurities underscores the need for computational approaches to assess the carcinogenicity risk of nitrosamines. However, current approaches are limited because robust animal carcinogenicity data are only available for a few simple nitrosamines, which do not represent the structural diversity of the many possible nitrosamine drug substance related impurities (NDSRIs). In this paper, we present a novel method that uses data on CYP-mediated metabolic hydroxylation of CH(2) groups in non-nitrosamine xenobiotics to identify structural features that may also help in predicting the likelihood of metabolic α-carbon hydroxylation in N-nitrosamines. Our approach offers a new avenue for tapping into potentially large experimental data sets on xenobiotic metabolism to improve risk assessment of nitrosamines. As α-carbon hydroxylation is the crucial rate-limiting step in nitrosamine metabolic activation, identifying and quantifying the influence of various structural features on this step can provide valuable insights into their carcinogenic potential. This is especially important considering the scarce information available on factors that affect NDSRI metabolic activation. We have identified hundreds of structural features and calculated their impact on hydroxylation, a significant advancement compared to the limited findings from the small nitrosamine carcinogenicity data set. While relying solely on α-carbon hydroxylation prediction is insufficient for forecasting carcinogenic potency, the identified features can help in the selection of relevant structural analogues in read across studies and assist experts who, after considering other factors such as the reactivity of the resulting electrophilic diazonium species, can establish the acceptable intake (AI) limits for nitrosamine impurities.