4D-QSAR Analysis of a Set of Ecdysteroids and a Comparison to CoMFA Modeling

Abstract
The ecdysteroid-responsive Drosophila melanogaster BII cell line is a prototypical homologous inducible gene expression system. A training set of 71 ecdysteroids, for which the −log(EC50) potencies in the ecdysteroid-responsive BII cell line were measured, was used to construct 4D-QSAR models. Four nearly equivalent optimum 4D-QSAR models, for two modestly different alignments, were identified (Q2 = 0.76−0.80). These four models, together with two CoMFA models, were used in consensus modeling to arrive at a three-dimensional pharmacophore. The C-2 and C-22 hydroxyls are identified as hydrogen-bond acceptor sites which enhance activity. A hydrophobic site near C-12 is consistent with increasing activity. The side-chain substituents at C-17 are predicted to adopt semiextended “active” conformations which could fit into a cylinder-shaped binding pocket lined largely with nonpolar residues for enhanced activity. A test set of 20 ecdysteroids was used to evaluate the QSAR models. Two 4D-QSAR models for one alignment were identified to be superior to the others based on having the smallest average residuals of prediction for the prediction set (0.69 and 1.13 −log[EC50] units). The correlation coefficients of the optimum 4D-QSAR models (R2 = 0.87 and 0.88) are nearly the same as those of the best CoMFA model (R2 = 0.92) determined for the same training set. However, the cross-validation correlation coefficient of the CoMFA model is less significant (Q2 = 0.59) than those of the 4D-QSAR models (Q2 = 0.80 and 0.80).

This publication has 34 references indexed in Scilit: