Probing the origins of aromatase inhibitory activity of disubstituted coumarins via QSAR and molecular docking

This study investigated the quantitative structure-activity relationship (QSAR) of imidazole derivatives of 4,7-disubstituted coumarins as inhibitors of aromatase, a potential therapeutic protein target for the treatment of breast cancer. Herein, a series of 3,7- and 4,7-disubstituted coumarin derivatives (1-34) with R1 and R2 substituents bearing aromatase inhibitory activity were modeled as a function of molecular and quantum chemical descriptors derived from low-energy conformer geometrically optimized at B3LYP/6-31G(d) level of theory. Insights on origins of aromatase inhibitory activity was afforded by the computed set of 7 descriptors comprising of F10[N-O], Inflammat-50, Psychotic-80, H-047, BELe1, B10[C-O] and MAXDP. Such significant descriptors were used for QSAR model construction and results indicated that model 4 afforded the best statistical performance. Good predictive performance were achieved as verified from the internal (comprising the training and the leave-one-out cross-validation (LOO-CV) sets) and external sets affording the following statistical parameters: R2Tr = 0.9576 and RMSETr = 0.0958 for the training set; Q2CV = 0.9239 and RMSECV = 0.1304 for the LOO-CV set as well as Q2Ext = 0.7268 and RMSEExt = 0.2927 for the external set. Significant descriptors showed correlation with functional substituents, particularly, R1 in governing high potency as aromatase inhibitor. Molecular docking calculations suggest that key residues interacting with the coumarins were predominantly lipophilic or non-polar while a few were polar and positively-charged. Findings illuminated herein serve as the impetus that can be used to rationally guide the design of new aromatase inhibitors

Molecular Docking of Aromatase Inhibitors

Aromatase is an enzyme that plays a critical role in the development of estrogen receptor positive breast cancer. As aromatase catalyzes the aromatization of androstenedione to estrone, a naturally occurring estrogen, it is a promising drug target for therapeutic management. The undesirable effects found in aromatase inhibitors (AIs) that are in clinical use necessitate the discovery of novel AIs with higher selectivity, less toxicity and improving potency. In this study, we elucidate the binding mode of all three generations of AI drugs to the crystal structure of aromatase by means of molecular docking. It was demonstrated that the docking protocol could reliably reproduce the interaction of aromatase with its substrate with an RMSD of 1.350 Å. The docking study revealed that polar (D309, T310, S478 and M374), aromatic (F134, F221 and W224) and non-polar (A306, A307, V370, L372 and L477) residues were important for interacting with the AIs. The insights gained from the study herein have great potential for the design of novel AIs