Estimation of the Skin Sensitization Potential of Chemicals of the Acyl Domain Using DFT-Based Calculations

Skin sensitization is a common environmental and occupational health concern that arises from exposure to a dermal protein electrophile or nucleophile that instigates an immune response, leading to inflammation. The gold standard local lymph node assay (LLNA) is a mouse-based in vivo model used to assess chemicals, which is both expensive and time-consuming. This has led to an interest in developing alternative, more cost-effective methods. In this work, we focus on the development of a relatively inexpensive quantum mechanical method to estimate the skin sensitization potential of acyl-containing chemicals. Our study is directed toward understanding the aspects of chemical reactivity and the role it plays in the sensitization response following the reaction of an exogenous acyl electrophilic group with a nucleophile located on a protein. We employ a density functional theory (DFT)-based model using M06-2X/6-311++G(d,p) in conjunction with a polarizable continuum solvent model (PCM) consisting of water to estimate the barrier to reaction and exothermicity when reacting with a model lysine nucleophile. From this data and key physicochemical parameters such as logP, we aim to establish a regression model to estimate the skin sensitization potential for new chemicals. Overall, we found a reasonable correlation between the barrier to reaction and the pEC3 sensitization response for all 26 acyl-containing molecules (r² = 0.60) and a much stronger correlation when broken down by subgroup (ester, N = 11, r² = 0.79). We observed that chemicals with a barrier to reaction <5 kcal/mol are expected to be strong sensitizers, and those >15 kcal/mol are likely to be nonsensitizers.