Cytotoxicity and 2D-QSAR study of some heterocyclic compounds

Herein we have studied the cytotoxicity and quantitative structure–activity relationship (QSAR) of heterocyclic compounds containing cyclic urea and thiourea nuclei. A set of 22 hydantoin and thiohydantoin related heterocyclic compounds were investigated with respect to their LC₅₀ values (Log of LC₅₀) against brine shrimp lethality bioassay in order to derive the 2D-QSAR models using MLR, PLS and ANN methods. The best predictive models by MLR, PLS and ANN methods gave highly significant square correlation coefficient (R²) values of 0.83, 0.81 and 0.91 respectively. The model also exhibited good predictive power confirmed by the high value of cross validated correlation coefficient Q² (0.74).     

Docking-based CoMFA and CoMSIA study of azaindole carboxylic acid derivatives as promising HIV-1 integrase inhibitors

Three-dimensional quantitative structure–activity relationship (3D-QSAR) studies were performed based on a series of azaindole carboxylic acid derivatives that had previously been reported as promising HIV-1 integrase inhibitors. Docking studies to explore the binding mode were performed based on the highly active molecule 36. The best docked conformation of molecule 36 was used as template for alignment. The comparative molecular field analysis (CoMFA) model (including steric and electrostatic fields) yielded the cross validation q ² = 0.655, non-cross validation r ² = 0.989 and predictive r ² _pred = 0.979. The best comparative molecular similarity indices analysis (CoMSIA) model (including steric, electrostatic, hydrophobic and hydrogen-bond acceptor fields) yielded the cross validation q ² = 0.719, non-cross validation r ² = 0.992 and predictive r ² _pred = 0.953. A series of new azaindole carboxylic acid derivatives were designed and the HIV-1 integrase inhibitory activities of these designed compounds were predicted based on the CoMFA and CoMSIA models.     

DFT‐based QSAR study and molecular design of AHMA derivatives as potent anticancer agents

A quantitative structure–activity relationship (QSAR) of 3‐(9‐acridinylamino)‐5‐hydroxymethylaniline (AHMA) derivatives and their alkylcarbamates as potent anticancer agents has been studied using density functional theory (DFT), molecular mechanics (MM+), and statistical methods. In the best established QSAR equation, the energy (E_NL) of the next lowest unoccupied molecular orbital (NLUMO) and the net charges (Q_FR) of the first atom of the substituent R, as well as the steric parameter (MR₂) of subsituent R₂ are the main independent factors contributing to the anticancer activity of the compounds. A new scheme determining outliers by “leave‐one‐out” (LOO) cross‐validation coefficient (q) was suggested and successfully used. The fitting correlation coefficient (R²) and the “LOO” cross‐validation coefficient (q²) values for the training set of 25 compounds are 0.881 and 0.829, respectively. The predicted activities of 5 compounds in the test set using this QSAR model are in good agreement with their experimental values, indicating that this model has excellent predictive ability. Based on the established QSAR equation, 10 new compounds with rather high anticancer activity much greater than that of 34 compounds have been designed and await experimental verification. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

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.     

3D－QSAR Study on Diindolylmethane and Its Analogues with Comparative Molecular Field Analysis(CoMFA)

Comparative molecular field analysis (CoMFA),a three dimensional quantitative structure-activity relationship (3D-QSAR) method was applied to a series of diindolylmethane(DIM) analogs to study the relationship between their structure and their induction of CYP 1A1-associated ethoxyresorufin-O-deethylase(EROD) activity.A DISCO model of pharmacophore was derved to guide the superposition of the compounds.The coefficient of cross-validation (q^2) and non cross-validation(r^2) for the model established by the study are 0.827 and 0.988 respectively,the value of variance ratio (F) is 103.53 and standard error estimate (SEE)is 0.044.These values indicate that the CoMFA model derived is significant and might have a good prediction for the catalytic activity of DIM compounds.As a consequence,the predicted activity values of new designed compounds were all higher than that of the reported value.     

Molecular modeling studies of human immunodeficiency virus type 1 protease inhibitors using three‐dimensional quantitative structure‐activity relationship,virtual screening,and docking simulations

In this paper, two 3‐dimensional quantitative structure‐activity relationship models for 60 human immunodeficiency virus (HIV)‐1 protease inhibitors were established using random sampling analysis on molecular surface and translocation comparative molecular field vector analysis (Topomer CoMFA). The non–cross‐validation (r²), cross‐validation (q²), correlation coefficient of external validation (Q²_ext), and F of 2 models were 0.94, 0.80, 0.79, and 198.84 and 0.94, 0.72, 0.75, and 208.53, respectively. The results indicated that 2 models were reasonable and had good prediction ability. Topomer Search was used to search R groups in the ZINC database, 20 new compounds were designed, and the Topomer CoMFA model was used to predicate the biological activity. The results showed that 18 new compounds were more active than the template molecule. So the Topomer Search is effective in screening and can guide the design of new HIV/AIDS drugs. The mechanism of action was studied by molecular docking, and it showed that the protease inhibitors and Ile50, Asp25, and Arg8 sites of HIV‐1 protease have interactions. These results have provided an insight for the design of new potent inhibitors of HIV‐1 protease.     

Using the QSPR Approach for Estimating the Density of Azole‐based Energetic Compounds

The relationship between density of energetic azole‐based compounds and their molecular structure is investigated through quantitative structure‐property relationship (QSPR) approach. The methodology of this work introduces a new model, which related density of azole‐based energetic compounds to the optimum elemental composition, the degree of unsaturation (DoU) of the compounds, presence of nitroimino group in the structural formula, as well as several non‐additive structural parameters. The presence of nitroimino functional group and also increasing the value of n_O/n_N in the formula of these compounds can enhance their density. The correlation is derived on the basis of experimental density values of 100 azole‐based energetic compounds with different molecular structure as training set. The determination coefficient of the new correlation is 0.923. Also, it has the root mean square deviation (RMSD) and the average absolute deviation (AAD) of 0.038 and 0.030 g · cm^–3, respectively. In addition, the correlation gives good predictions for further 25 azole‐based energetic compounds as test set (Q²_EXT = 0.901). The predictive ability of the correlation is checked using a cross validation method (Q²_LMO = 0.918). The proposed method can also apply for designing novel azole‐based energetic compounds.     

Exploring the molecular basis for selective cytotoxicity of lamellarins against human hormone-dependent T47D and hormone-independent MDA-MB-231 breast cancer cells

Abstract  

The common structural requirements for cytotoxicity of lamellarins against two human breast cancer cell lines were determined using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) techniques. Twenty lamellarins were selected to serve as the training set, whereas another group of six compounds were used as the test set. The best CoMFA and CoMSIA models for both cell lines yielded satisfactory predictive ability with r _cv² values in the range of 0.659–0.728. Additionally, the contour maps obtained from both the CoMFA and CoMSIA models agreed well with the experimental results and may be used in the design of more potent cytotoxic compounds for human breast cancers. Both analyses not only suggested structural requirements of various substituents around the lamellarin skeleton for their cytotoxic activity against both human breast cancer cell lines but also revealed the molecular basis for the differences between the saturated and unsaturated D-rings of the lamellarins.