Monte Carlo-based QSAR modeling of dimeric pyridinium compounds and drug design of new potent acetylcholine esterase inhibitors for potential therapy of <Emphasis Type="Italic">myasthenia gravis</Emphasis>

The Monte Carlo method was used for QSAR modeling of dimeric pyridinium compounds as acetylcholine esterase inhibitors. QSAR model was developed for a series of 39 dimeric pyridinium compounds. QSAR models were calculated with the representation of the molecular structure by the simplified molecular-input line-entry system. One split into the training and test set have been examined. The statistical quality of the developed model is very good. The calculated model for dimeric pyridinium derivatives had following statistical parameters: r ² = 0.9477 for the training set and r ² = 0.9332 the test set. Structural indicators considered as molecular fragments responsible for the increase and decrease in the inhibition activity have been defined. The computer-aided design of new dimeric pyridinium compounds potential acetylcholine esterase inhibitors with the application of defined structural alerts has been presented.     

Predicting skin permeability from complex chemical mixtures: incorporation of an expanded QSAR model

Quantitative structure–activity relationship (QSAR) models have been widely used to study the permeability of chemicals or solutes through skin. Among the various QSAR models, Abraham’s linear free-energy relationship (LFER) model is often employed. However, when the experimental conditions are complex, it is not always appropriate to use Abraham’s LFER model with a single set of regression coefficients. In this paper, we propose an expanded model in which one set of partial slopes is defined for each experimental condition, where conditions are defined according to solvent: water, synthetic oil, semi-synthetic oil, or soluble oil. This model not only accounts for experimental conditions but also improves the ability to conduct rigorous hypothesis testing. To more adequately evaluate the predictive power of the QSAR model, we modified the usual leave-one-out internal validation strategy to employ a leave-one-solute-out strategy and accordingly adjust the Q² _LOO statistic. Skin permeability was shown to have the rank order: water > synthetic > semi-synthetic > soluble oil. In addition, fitted relationships between permeability and solute characteristics differ according to solvents. We demonstrated that the expanded model (r² = 0.70) improved both the model fit and the predictive power when compared with the simple model (r² = 0.21).     

Structural investigations of anthranilimide derivatives by CoMFA and CoMSIA 3D-QSAR studies reveal novel insight into their structures toward glycogen phosphorylase inhibition

In the present work, three-dimensional quantitative structure–activity relationship (3-D QSAR) studies on a set of 70 anthranilimide compounds has been performed using docking-based as well as substructure-based molecular alignments. This resulted in the selection of more statistically relevant substructure-based alignment for further studies. Further, molecular models with good predictive power were derived using CoMFA (r ^2?=?0.997; Q ^2?=?0.578) and CoMSIA (r ^2?=?0.976; Q ^2?=?0.506), for predicting the biological activity of new compounds. The so-developed contour plots identified several key features of the compounds explaining wide activity ranges. Based on the information derived from the CoMFA contour maps, novel leads were proposed which showed better predicted activity with respect to the already reported systems. Thus, the present study not only offers a highly significant predictive QSAR model for anthranilimide derivatives as glycogen phosphorylase (GP) inhibitors which can eventually assist and complement the rational drug-design attempts, but also proposes a highly predictive pharmacophore model as a guide for further development of selective and more potent GP inhibitors as anti-diabetic agents.     

Combined 2D and 3D-QSAR,molecular modelling and docking studies of pyrazolodiazepinones as novel phosphodiesterase 2 inhibitors

Selective inhibition of phosphodiesterase 2 (PDE2) in cells where it is located elevates cyclic guanosine monophosphate (cGMP) and acts as novel analgesic with antinociceptive activity. Three-dimensional quantitative structure–activity relationship (QSAR) studies for pyrazolodiazepinone inhibitors exhibiting PDE2 inhibition were performed using comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA) and Topomer CoMFA, and two-dimensional QSAR study was performed using a Hologram QSAR (HQSAR) method. QSAR models were generated using training set of 23 compounds and were validated using test set of nine compounds. The optimum partial least squares (PLS) for CoMFA-Focusing, CoMSIA-SDH, Topomer CoMFA and HQSAR models exhibited good ‘leave-one-out’ cross validated correlation coefficient (q²) of 0.790, 0.769, 0.840 and 0.787, coefficient of determination (r²) of 0.999, 0.964, 0.979 and 0.980, and high predictive power (r²_pred) of 0.796, 0.833, 0.820 and 0.803 respectively. Docking studies revealed that those inhibitors able to bind to amino acid Gln859 by cGMP binding orientation called ‘glutamine-switch’, and also bind to the hydrophobic clamp of PDE2 isoform, could possess high selectivity for PDE2. From the results of all the studies, structure–activity relationships and structural requirements for binding to active site of PDE2 were established which provide useful guidance for the design and future synthesis of potent PDE2 inhibitors.     

Hologram QSAR Models of 4-[(Diethylamino)methyl]-phenol Inhibitors of Acetyl/Butyrylcholinesterase Enzymes as Potential Anti-Alzheimer Agents

Hologram QSAR models were developed for a series of 36 inhibitors (29 training set and seven test set compounds) of acetyl/butyrylcholinesterase (AChE/BChE) enzymes, an attractive molecular target for Alzheimer's disease (AD) treatment. The HQSAR models (N = 29) exhibited significant cross-validated (AChE, q² = 0.787; BChE, q² = 0. 904) and non-cross-validated (AChE, r² = 0.965; BChE, r² = 0.952) correlation coefficients. The models were used to predict the inhibitory potencies of the test set compounds, and agreement between the experimental and predicted values was verified, exhibiting a powerful predictive capability. Contribution maps show that structural fragments containing aromatic moieties and long side chains increase potency. Both the HQSAR models and the contribution maps should be useful for the further design of novel, structurally related cholinesterase inhibitors.     

CoMFA and CoMSIA studies on C-aryl glucoside SGLT2 inhibitors as potential anti-diabetic agents

SGLT2 has become a target of therapeutic interest in diabetes research. CoMFA and CoMSIA studies were performed on C-aryl glucoside SGLT2 inhibitors (180 analogues) as potential anti-diabetic agents. Three different alignment strategies were used for the compounds. The best CoMFA and CoMSIA models were obtained by means of Distill rigid body alignment of training and test sets, and found statistically significant with cross-validated coefficients (q ²) of 0.602 and 0.618, respectively, and conventional coefficients (r ²) of 0.905 and 0.902, respectively. Both models were validated by a test set of 36 compounds giving satisfactory predicted correlation coefficients (r ² _pred) of 0.622 and 0.584 for CoMFA and CoMSIA models, respectively. A comparison was made with earlier 3D QSAR study on SGLT2 inhibitors, which shows that our 3D QSAR models are better than earlier models to predict good inhibitory activity. CoMFA and CoMSIA models generated in this work can provide useful information to design new compounds and helped in prediction of activity prior to synthesis.     

Impact assessment of the rational selection of training and test sets on the predictive ability of QSAR models

This study performed an analysis of the influence of the training and test set rational selection on the quality and predictively of the quantitative structure–activity relationship (QSAR) model. The study was carried out on three different datasets of Influenza Neuraminidase (H1N1) inhibitors. The three datasets were divided into training and test sets using three rational selection methods: based on k-means, Kennard–Stone algorithm and Activity and the results were compared with Random selection. Then, a total of 31,490 mathematical models were developed and those models that presented a determination coefficient higher than: r²_train > 0.8, r²_loo > 0.7, r²_test > 0.5 and minimum standard deviation (SD) and minimum root-mean square error (RMS) were selected. The selected models were validated using the internal leave-one-out method and the predictive capacity was evaluated by the external test set. The results indicate that random selection could lead to erroneous results. In return, a rational selection allows for obtaining more reliable conclusions. The QSAR models with major predictive power were found using the k-means algorithm and selection by activity.     

3D-Pharmacophore mapping of thymidine-based inhibitors of TMPK as potential antituberculosis agents

Tuberculosis (TB) is the primary cause of mortality among infectious diseases. Mycobacterium tuberculosis monophosphate kinase (TMPKmt) is essential to DNA replication. Thus, this enzyme represents a promising target for developing new drugs against TB. In the present study, the receptor-independent, RI, 4D-QSAR method has been used to develop QSAR models and corresponding 3D-pharmacophores for a set of 81 thymidine analogues, and two corresponding subsets, reported as inhibitors of TMPKmt. The resulting optimized models are not only statistically significant with r ² ranging from 0.83 to 0.92 and q ² from 0.78 to 0.88, but also are robustly predictive based on test set predictions. The most and the least potent inhibitors in their respective postulated active conformations, derived from each of the models, were docked in the active site of the TMPKmt crystal structure. There is a solid consistency between the 3D-pharmacophore sites defined by the QSAR models and interactions with binding site residues. Moreover, the QSAR models provide insights regarding a probable mechanism of action of the analogues.     

Insight into the structural requirements of benzimidazole derivatives as interleukin‐2 inducible t‐cell kinase inhibitors by computational explorations

In the present work, a set of ligand‐ and receptor‐based 3D‐QSAR models were developed to explore the structure–activity relationship of 109 benzimidazole‐based interleukin‐2‐inducible T‐cell kinase (ITK) inhibitors. In order to reveal the requisite 3D structural features impacting the biological activities, a variety of in silico modeling approaches including the comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), docking, and molecular dynamics were applied. The results showed that the ligand‐based CoMFA model (Q² = 0.552, R²_ncv = 0.908, R²_pred = 0.787, SEE = 0.252, SEP = 0.558) and CoMSIA model (Q² = 0.579, R²_ncv = 0.914, R²_pred = 0.893, SEE = 0.240, SEP = 0.538) were superior to other models with greater predictive power. In addition, a combined analysis between the 3D contour maps and docking results showed that: (1) Compounds with bulky or hydrophobic substituents near ring D and electropositive or hydrogen acceptor groups around rings C and D could increase the activity. (2) The key amino acids impacting the receptor–ligand interactions in the binding pocket are Met438, Asp500, Lys391, and Glu439. The results obtained from this work may provide helpful guidelines in design of novel benzimidazole analogs as inhibitors of ITK. © 2013 Wiley Periodicals, Inc.