Computational models to predict blood–brain barrier permeation and CNS activity

The blood-brain permeation of a structurally diverse set of 281 compounds was modeled using linear regression and a multivariate genetic partial least squares (G/PLS) approach. Key structural features affecting the logarithm of blood-brain partitioning (logBB) were captured through statistically significant quantitative structure-activity relationship (QSAR) models. These relationships reveal the importance of logP, polar surface area, and a variety of electrotopological indices for accurate predictions of logBB. The best models reveal an excellent correlation (r > 0.9) for a training set of 58 compounds. Likewise, the comparison of the average logBB values obtained from an ensemble of QSAR models with experimental values also verifies the statistical quality of the models (r > 0.9). The models provide good agreement (r approximately 0.7) between the predicted logBB values for 34 molecules in the external validation set and the experimental values. To further validate the models for use during the drug discovery process, a prediction set of 181 drugs with reported CNS penetration data was used. A >70% success rate is obtained by using any of the QSAR models in the qualitative prediction for CNS permeable (active) drugs. A lower success rate (approximately 60%) was obtained for the best model for CNS impermeable (inactive) drugs. Combining the predictions obtained from all the models (consensus) did not significantly improve the discrimination of CNS active and CNS inactive molecules. Finally, using the therapeutic classification as a guiding tool, the CNS penetration capability of over 2000 compounds in the Synthline database was estimated. The results were very similar to the smaller set of 181 compounds.     

A new protocol for predicting novel GSK-3β ATP competitive inhibitors

Glycogen synthase kinase 3β (GSK-3β) is a potential therapeutic target for cancer, type-2 diabetes, and Alzheimer's disease. This paper proposes a new lead identification protocol that predicts new GSK-3β ATP competitive inhibitors with topologically diverse scaffolds. First, three-dimensional quantitative structure-activity relationship (3D QSAR) models were built and validated. These models are based upon known GSK-3β inhibitors, benzofuran-3-yl-(indol-3-yl) maleimides, by means of comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Second, 28?826 maleimide derivatives were selected from the PubChem database. After filtration via Lipinski's rules, 10?429 maleimide derivatives were left. Third, the FlexX-dock program was employed to virtually screen the 10?429 compounds against GSK-3β. This resulted in 617 virtual hits. Fourth, the 3D QSAR models predicted that from the 617 virtual hits, 93 compounds would have GSK-3β inhibition values of less than 15 nM. Finally, from the 93 predicted active hits, 23 compounds were confirmed as GSK-3β inhibitors from literatures; their GSK-3β inhibition ranged from 1.3 to 480 nM. Therefore, the hits rate of our virtual screening protocol is greater than 25%. The protocol combines ligand- and structure-based approaches and therefore validates both approaches and is capable of identifying new hits with topologically diverse scaffolds.     

QSAR of conformationally flexible molecules: Comparative Molecular Field Analysis of protein-tyrosine kinase inhibitors

Summary Comparative Molecular Field Analysis (CoMFA) has been applied to a study of quantitative structureactivity relationships (QSAR) of conformationally flexible molecules. The relationship between three-dimensional structure and activity of 20 styrene derivatives which inhibit protein-tyrosine kinase was determined. A technique was developed that allows accurate prediction of the inhibitory activity of these molecules and identification in each case of the active conformation. The problem of multiple energetically acceptable conformations was approached in an iterative procedure. Use was made of the varying degrees of symmetry among the molecules. First, CoMFA QSAR models were developed using only those compounds that possess a symmetrical substituent pattern on the phenyl ring. These CoMFA models were then used to select the active conformers of the less symmetrical compounds in the set. Allowing multiple conformers for each compound in the dataset yielded higher crossvalidated r² values and better predictivity of the QSAR models. Different probe atoms (C⁺, O^–, neutral C) were explored, the O^– probe atom exhibiting the highest selectivity in the conformer selection process.     

4D-QSAR analysis of a set of ecdysteroids and a comparison to CoMFA modeling.

The ecdysteroid-responsive Drosophila melanogaster B(II) cell line is a prototypical homologous inducible gene expression system. A training set of 71 ecdysteroids, for which the -log(EC(50)) potencies in the ecdysteroid-responsive B(II) 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 (Q(2) = 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[EC(50)] units). The correlation coefficients of the optimum 4D-QSAR models (R(2) = 0.87 and 0.88) are nearly the same as those of the best CoMFA model (R(2) = 0.92) determined for the same training set. However, the cross-validation correlation coefficient of the CoMFA model is less significant (Q(2) = 0.59) than those of the 4D-QSAR models (Q(2) = 0.80 and 0.80).