Validation of a Fragment-Based Profiler for Thiol Reactivity for the Prediction of Toxicity: Skin Sensitization and Tetrahymena pyriformis

Eliminating Transition State Calculations for Faster and More Accurate Reactivity Prediction in Sulfa-Michael Additions Relevant to Human Health and the Environment

Fast and accurate computational approaches to predicting reactivity in sulfa-Michael additions are required for high-throughput screening in toxicology (e.g., predicting excess aquatic toxicity and skin sensitization), chemical synthesis, covalent drug design (e.g., targeting cysteine), and data set generation for machine learning. The kinetic glutathione chemoassay is a time-consuming in chemico method used to extract kinetic data in the form of log(k _GSH) for organic electrophiles. In this work, we use density functional theory to compare the use of transition states (TSs) and enolate intermediate structures following C-S bond formation in the prediction of log(k _GSH) for a diverse group of 1,4 Michael acceptors. Despite the widespread use of transition state calculations in the literature to predict sulfa-Michael reactivity, we observe that intermediate structures show much better performance for the prediction of log(k _GSH), are faster to calculate, and easier to obtain than TSs. Furthermore, we show how linear combinations of atomic charges from the isolated Michael acceptors can further improve predictions, even when using inexpensive semiempirical quantum chemistry methods. Our models can be used widely in the chemical sciences (e.g., in the prediction of toxicity relevant to the environment and human health, synthesis planning, and the design of cysteine-targeting covalent inhibitors), and represent a low-cost, sustainable approach to reactivity assessment.

Using Transition State Modeling To Predict Mutagenicity for Michael Acceptors

The Ames mutagenicity assay is a long established in vitro test to measure the mutagenicity potential of a new chemical used in regulatory testing globally. One of the key computational approaches to modeling of the Ames assay relies on the formation of chemical categories based on the different electrophilic compounds that are able to react directly with DNA and form a covalent bond. Such approaches sometimes predict false positives, as not all Michael acceptors are found to be Ames-positive. The formation of such covalent bonds can be explored computationally using density functional theory transition state modeling. We have applied this approach to mutagenicity, allowing us to calculate the activation energy required for α,β-unsaturated carbonyls to react with a model system for the guanine nucleobase of DNA. These calculations have allowed us to identify that chemical compounds with activation energies greater than or equal to 25.7 kcal/mol are not able to bind directly to DNA. This allows us to reduce the false positive rate for computationally predicted mutagenicity assays. This methodology can be used to investigate other covalent-bond-forming reactions that can lead to toxicological outcomes and learn more about experimental results.