Rigorous 3-dimensional spectral data activity relationship approach modeling strategy for ToxCast estrogen receptor data classification, validation, and feature extraction

Determination of structural factors affecting binding to mu, kappa and delta opioid receptors

Addiction is a complex behavioral phenomenon in which naturally occurring or synthetic chemicals modulate the response of the reward system through their binding to a variety of neuroreceptors, resulting in compulsive substance-seeking and use despite harmful consequences to the individual. Among these, the opioid receptor (OR) family and more specifically, the mu-opioid receptor (MOR) subtype plays a critical role in the addiction to powerful prescription and illicit drugs such as hydrocodone, oxycodone, fentanyl, cocaine, and methamphetamine (Contet et al. in Curr Opin Neurobiol 14(3):370-378, 2004). Conversely, agonists binding to kappa (KOR) and antagonists binding to delta opioid receptors (DOR) have been reported to induce negative reinforcing effects. As more than 700 new psychoactive substances were illegally sold between 2009 and 2016 (DEA-DCT-DIR-032-18), most of them lacking basic toxicological and pharmacological profiles, molecular modeling approaches that could quickly and reliably fill the gaps in our knowledge would be highly desirable tools for determining the effects of these synthetics. Here, we report accurate 3D-spectrometric data-activity relationship classification models for large and diverse datasets of MOR, KOR and DOR binders with areas under the receiver operating characteristic curve for the "blind" prediction sets exceeding 0.88. Structural features associated with (selective) binding to MOR, KOR and/or DOR were identified. These models could assist regulatory agencies in evaluating the health risks associated with the use of unprofiled substances as well as to help the pharmaceutical industry in its search for new drugs to combat addiction.

Computational identification of structural factors affecting the mutagenic potential of aromatic amines: study design and experimental validation

A grid-based, alignment-independent 3D-SDAR (three-dimensional spectral data-activity relationship) approach based on simulated ¹³C and ¹⁵N NMR chemical shifts augmented with through-space interatomic distances was used to model the mutagenicity of 554 primary and 419 secondary aromatic amines. A robust modeling strategy supported by extensive validation including randomized training/hold-out test set pairs, validation sets, "blind" external test sets as well as experimental validation was applied to avoid over-parameterization and build Organization for Economic Cooperation and Development (OECD 2004) compliant models. Based on an experimental validation set of 23 chemicals tested in a two-strain Salmonella typhimurium Ames assay, 3D-SDAR was able to achieve performance comparable to 5-strain (Ames) predictions by Lhasa Limited's Derek and Sarah Nexus for the same set. Furthermore, mapping of the most frequently occurring bins on the primary and secondary aromatic amine structures allowed the identification of molecular features that were associated either positively or negatively with mutagenicity. Prominent structural features found to enhance the mutagenic potential included: nitrobenzene moieties, conjugated π-systems, nitrothiophene groups, and aromatic hydroxylamine moieties. 3D-SDAR was also able to capture "true" negative contributions that are particularly difficult to detect through alternative methods. These include sulphonamide, acetamide, and other functional groups, which not only lack contributions to the overall mutagenic potential, but are known to actively lower it, if present in the chemical structures of what otherwise would be potential mutagens.