Oxidation Potential of 2,6-Dimethyl-1,4-dihydropyridine Derivatives Estimated by Structure Descriptors

Abstract

Linear relationships, expressing the electrochemical properties of molecules as functions of structure, give insight into the associated electrochemical process and are a tool for prediction. Many biological activities rely on water-based dissociation, making electrochemical properties a bridge between structure and activity. Motivated by a previous study, a replica is made here on a different dataset in order to validate/invalidate the previously reported results. There are several methods for obtaining structure-based descriptors. Some of the methods have been devised to account for molecular topology, some to account for molecular geometry, and others to account for both. Two methods are involved here to derive structure-based descriptors and further obtain structure–property relationships (FMPI and ChPE). In order to express structure descriptors, both FMPI and ChPE express first the topology of the molecule, using the heavy atoms identity matrix and the heavy atoms adjacency matrix, both square symmetric matrices in the belief that symmetry is one major factor of molecular stability. A set of 2,6-dimethyl-1,4-dihydropyridine derivatives with oxidation peak potentials and coulometrically determined number of electrons experimental data is subjected to the search for structure–activity relationships. Even if the 2,6-dimethyl-1,4-dihydropyridine is a symmetric compound (of Cs point group), their derivatives are generally not symmetric (9 out of 24 are asymmetric). The dataset is subjected to descriptive and inferential statistics in order to filter out the most relevant structure–activity relationship. The geometry is built using three levels of theory (one from molecular mechanics and two others from density functionals, of which one accounts for the interaction with water as solvent). One challenge of picking one out of two reported measured values is dealt with by calculating the likelihood associated with the two choices. Relevant structure–activity models are extracted and discussed. The use of in vivo (in water, SM8 model) models in geometry optimization (from MMFF94 and B3LYP and to M06 + Water SM8) results in a precision gain, but this is, in most of the cases, not statistically significant, and this can be considered a negative result.