2D and 3D‐QSAR studies on antiproliferative thiazolidine analogs

Two‐dimensional (2D) and three‐dimensional (3D) quantitative structure–activity relationships (QSARs) of 22 thiazolidine analogs with antiproliferative activity expressed as pIC₅₀, which is defined as the negative value of the logarithm of necessary molar concentration of these compounds to cause 50% growth inhibition against melanoma cell lines WM‐164, have been studied by using a combined method of the DFT, MM2 and statistics for 2D, as well as the comparative molecular field analysis (CoMFA) method for 3D. The established 2D‐QSAR model in training set comprised of random 18 compounds shows not only significant statistical quality, but also predictive ability, with the square of adjusted correlation coefficient (R = 0.832) and the square of the cross‐validation coefficient (q² = 0.803). The same model was further applied to predict pIC₅₀ values of the four compounds in the test set, and the resulting R reaching 0.784, further confirms that this 2D‐QSAR model has high predictive ability. The 3D‐QSAR model also shows good correlative and predictive capabilities in terms of R² (0.956) and q² (0.615) obtained from CoMFA model. Further, the robustness of the CoMFA model was verified by bootstrapping analysis (100 runs) with R (0.979) and SD_bs (0.056). It is very interesting to find that the results from 2D‐ and 3D‐QSAR analyses accord with each other, and they all show that the steric interaction plays a crucial role in determining the cytotoxicities of the compounds, and that selecting a moderate‐size or appropriate‐hydrophobicity substituent R as well as increasing the negative charges of C₄ on phenyl ring at the same time are advantageous to improving the cytotoxicity. Such results can offer some useful theoretical references for directing the molecular design and understanding the action mechanism of this kind of compound with antiproliferative activity. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Grid potential analysis,virtual screening studies and ADME/T profiling on N‐arylsulfonylindoles as anti‐HIV‐1 agents

A grid potential analysis employing a novel approach of 3D quantitative structure–activity relationships (QSAR) as AutoGPA module in MOE2009.10 was performed on a dataset of 42 compounds of N‐arylsulfonylindoles as anti‐HIV‐1 agents. The uniqueness of AutoGPA module is that it automatically builds the 3D‐QSAR model on the pharmacophore‐based molecular alignment. The AutoGPA‐based 3D‐QSAR model obtained in the present study gave the cross‐validated Q² value of 0.588, r²_pred value of 0.701, r²_m statistics of 0.732 and Fisher value of 94.264. The results of 3D‐QSAR analysis indicated that hydrophobic groups at R₁ and R₂ positions and electron releasing groups at R₃ position are favourable for good activity. To find similar analogues, virtual screening on ZINC database was carried out using generated AutoGPA‐based 3D‐QSAR model and showed good prediction. In addition to those mentioned earlier, in‐silico ADME absorption, distribution, metabolism and excretion profiling and toxicity risk assessment test was performed, and results showed that majority of compounds from current dataset and newly virtually screened hits generated were within their standard limit. Copyright © 2013 John Wiley & Sons, Ltd.     

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     

Pharmacophore mapping and electronic feature analysis for a series of nitroaromatic compounds with antitubercular activity

A five point pharmacophore was generated using PHASE for a series of nitroaromatic compounds and their congeners as antitubercular agents. The generated pharmacophore yielded significant 3D‐QSAR model with r² of 0.890 for a training set of 92 molecules. The model also showed excellent predictive power with correlation coefficient Q² of 0.857 for a test set of 31 compounds. The pharmacophore indicated that presence of a nitro group, a piperazine moiety, one aromatic ring feature and two acceptor features are necessary for potent antitubercular activity. The pharmacophore was supported by electronic property analysis using density functional theory (DFT) at B3LYP/3–21*G level. Molecular electrostatic profile of the compounds was consistent with the generated pharmacophore model, particularly appearance of localized negative potential regions near both the oxygen atoms of nitro group extending laterally to the isoxazole ring system/amide bond in the most active compounds. Calculated data further revealed that all active compounds have smaller LUMO energies located over the nitro group, furan ring, and isoxazole ring/amide bond attached to it. Higher negative values of LUMO energies concentrated over the nitro group are indicative of the electron acceptor capacity of the compounds, suggesting that these compounds are prodrugs and must be activated by TB‐nitroreductase. The results obtained from this study should aid in efficient design and development of nitroaromatic compounds as antitubercular agents. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Theoretical studies on QSAR and mechanism of 2‐indolinone derivatives as tubulin inhibitors

The theoretical studies on three‐dimensional quantitative structure activity relationship (3D‐QSAR) and action mechanism of a series of 2‐indolinone derivatives as tubulin inhibitors against human breast cancer cell line MDA‐MB‐231 have been carried out. The established 3D‐QSAR model from the comparative molecular field analysis (CoMFA) shows not only significant statistical quality but also predictive ability, with high correlation coefficient (R² = 0.986) and cross‐validation coefficient (q² = 0.683). In particular, the appropriate binding orientations and conformations of these 2‐indolinone derivatives interacting with tubulin are located by docking study, and it is very interesting to find that the plot of the energy scores of these compounds in DOCK versus the corresponding experimental pIC₅₀ values exhibits a considerable linear correlation. Therefore, the inhibition mechanism that 2‐indolinone derivatives were regarded as tubulin inhibitors can be theoretically confirmed. Based on such an inhibition mechanism along with 3D‐QSAR results, some important factors improving the activities of these compounds were discussed in detail. These factors can be summarized as follows: the H atom adopted as substituent R₁, the substituent R₂ with higher electropositivity and smaller bulk, the substituents R₄–R₆ (on the phenyl ring) with higher electropositivity and larger bulk, and so on. These results can offer useful theoretical references for understanding the action mechanism, designing more potent inhibitors, and predicting their activities prior to synthesis. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Three-dimensional Quantitative Structure-Activity Relationship Analyses of a Series of Butenolide ETA Antagonists

Two kinds of Three-dimensional Quantitative Structure-activity Relationship(3D-QSAR) methods,comparative molecular filed analysis(CoMFA) and comparative molecular similarity indices analysis (CoMSIA) ,were applied to analyze the structure-activity relationship of a series of 63 butenolide ETA selective antagonists with respect to their inhibition against human ETA receptor,The CoMFA and CoMSIA models were developed for the conceivable alignment of the molecules based on a template structure from the crystallized data.The statistical results from the initial orientation of the aligned molecules show that the 3D-QSAR model from CoMFA(q^2=0.543) is obviously superior to that from the conventional CoMSIA(q^2=0.407).In order to refine the model,all-space search (ASS) was applied to minimize the field sampling process.By rotating and translating the molecular aggregate within the grid systematically,all the possible samplings of the molecular fields were tested and subsequently the one with the highest q^2 was picked out .The comparison of the sensitivity of CoMFA and CoMSIA to different space orientation shows that the CoMFA q^2 values are more sensitive to the translations and rotations of the aligned molecules with respect to the lattice than those of CoMSIA.The best CoMFA model from ASS was further refined by the region focused technique.The high quality of the best model is indicated by the high corss-validated correlation and the prediction on the external test set.The CoMFA coefficient contour plots identify several key features that explain the wide range of activities,which may help us to design new effective ETA selective antagonists.     

How dependent are molecular and atomic properties on the electronic structure method? Comparison of Hartree‐Fock,DFT, and MP2 on a biologically relevant set of molecules

This article compares molecular properties and atomic properties defined by the quantum theory of atoms in molecules (QTAIM) obtained from three underlying levels of theory: MP2(full), density functional theory (DFT) (B3LYP), and Hartree‐Fock (H‐F). The same basis set (6‐311++G(d,p)) has been used throughout the study. The calculations and comparisons were applied to a set of 30 small molecules representing common fragments of biological molecules. The molecular properties investigated are the energies and the electrostatic moments (up to and including the quadrupoles), and the atomic properties include electron populations (and atomic charge), atomic dipolar and quadrupolar polarizations, atomic volumes, and corrected and raw atomic energies. The Cartesian distance between dipole vectors and the Frobenius distance between the quadrupole tensors calculated at the three levels of theory provide a measure of their correlation (or lack thereof). With the exception of energies (atomic and molecular), it is found that both DFT and H‐F are in excellent agreement with MP2, especially with regards to the electrostatic mutipoles up to the quadrupoles, but DFT and MP2 agree better in almost all studied properties (with the exception of molecular geometries). QTAIM properties whether obtained from H‐F, DFT(B3LYP), or MP2 calculations when used in the construction of empirical correlations with experiment such as quantitative structure‐activity‐(or property)‐relationships (QSAR/QSPR) are equivalent (because the properties calculated at the three levels are very highly correlated among themselves with r² typically >0.95, and therefore preserving trends). These results suggest that the massive volume of results that were published in the older literature at the H‐F level is valid especially when used to study trends or in QSAR or QSPR studies, and, as long as our test set of molecules is representative, there is no pressing need to re‐evaluate them at other levels of theory except when inadequate basis sets were used by today's standards. Extensive tabulation of molecular and atomic properties at the three theoretical levels is available in the Supporting Information, including optimized geometries, molecular energies, virial ratios, molecular electrostatic moments up to and including hexadecapoles, atomic populations, atomic volumes, atomic electrostatic moments up to and including the quadrupoles, and atomic energies. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Theoretical study of silicon–sulfur clusters (SiS2)n− (n=1–6)

The possible geometrical structures and relative stability of silicon–sulfur clusters (SiS₂) (n=1–6) are explored by means of density functional theory (DFT) quantum chemical calculations. We also compare DFT with second‐order Møller–Plesset (MP2) and Hartree–Fock (HF) methods. The effects of polarization functions, diffuse functions, and electron correlation are included in MP2 and B3LYP quantum chemical calculations, and B3LYP is effective in larger cluster structure optimization, so we can conclude that the DFT approach is useful in establishing trends. The electronic structures and vibrational spectra of the most stable geometrical structures of (SiS₂)_n⁻ are analyzed by B3LYP. As a result, the regularity of the (SiS₂)_n⁻ cluster growing is obtained, and the calculation may predict the formation mechanism of the (SiS₂)_n⁻ cluster. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 280–290, 2001     

DFT Study for a Series of Less-Symmetrical Crown Ethers and Their Complexes with Alkali Metal Cations

A series of ring‐contracted (14‐crown‐5, 17‐crown‐6) and ring‐enlarged (16‐crown‐5, 17‐crown‐5, 19‐crown‐6, 20‐crown‐6) crown ethers and their complexes with alkali‐metal cations Na⁺ and K⁺ had been explored using density functional theory (DFT) at B3LYP/6‐31G* level in order to reveal the effects of the methylene‐chain length in a crown ether. The nucleophilicity of all crown ethers had been investigated by the Fukui functions. The quantum chemistry parameters, such as the energy gap (ΔE), the highest occupied molecular orbital energy (E_HOMO) and the lowest unoccupied molecular orbital energy (E_LUMO) for less‐symmetrical crown ethers and symmetrical frameworks (15‐crown‐5, 18‐crown‐6) had been calculated. In addition, the thermodynamic energies of complexation reactions had also been studied. The results of the DFT calculations show that the methylene‐chain length plays an important role in determining the structure characters of the crown ethers and also strongly influences the properties of the ethers. Some of the calculated results are in a good agreement with the experimental values.     

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.     

Comparative molecular field analysis and comparative molecular similarity index analysis studies on 1H NMR chemical shift of NH group of diaryl triazene derivatives

Comparative molecular field analysis (CoMFA), comparative molecular field analysis region focusing (CoMFA‐RF) for optimizing the region for the final partial least square analysis, and comparative molecular similarity indices analysis (CoMSIA) methods were employed to develop three‐dimensional quantitative structure–activity relationship (3D‐QSAR) models of ¹H NMR chemical shift of NH proton of diaryl triazene derivatives. The best orientation was searched by all‐orientation search (AOS) strategy to minimize the effect of the initial orientation of the structures. The predictive abilities of CoMFA‐RF and CoMSIA models were determined using a test set of ten compounds affording predictive correlation coefficients of 0.721 and 0.754, respectively, indicating good predictive power. For further model validation, cross validation (leave one out), progressive scrambling, and bootstrapping were also applied. The accuracy and speed of obtained 3D‐QSAR models for the prediction of ¹H NMR chemical shifts of NH group of diaryl triazene derivatives were greater compared to some computational well‐known procedures. Copyright © 2013 John Wiley & Sons, Ltd.     

The Weakly Coordinating Tris(trichlorosilyl)silyl Anion

Closely following the procedure for the preparation of the base‐stabilized dichlorosilylene complex NHC^Dipp⋅SiCl₂ reported by Roesky, Stalke, and co‐workers (Angew. Chem. Int. Ed . 2009 , 48 , 5683–5686), a few crystals of the salt [NHC^Dipp−H⋅⋅⋅Cl⋅⋅⋅H−NHC^Dipp]Si(SiCl₃)₃ were isolated, aside from the reported byproduct [NHC^Dipp−H⁺⋅⋅⋅Cl⁻], and characterized by X‐ray crystallography (NHC^Dipp=N,N‐di(2,6‐diisopropylphenyl)imidazo‐2‐ylidene). They contain the weakly coordinating anion Si(SiCl₃)₃⁻, which was also obtained in high yields upon deprotonation of the conjugate Brønsted acid HSi(SiCl₃)₃ with NHC^Dipp or PMP (PMP=1,2,2,6,6‐pentamethylpiperidine). The acidity of HSi(SiCl₃)₃ was estimated by DFT calculations to be substantially higher than those of other H‐silanes. Further DFT studies on the electronic structure of Si(SiCl₃)₃⁻, including the electrostatic potential and the electron localizability, confirmed its low basicity and nucleophilicity compared with other silyl anions.     

Study on the Three-dimensional Quantitative Structure-activity Relationship of Pyridazinonyl-substituted 1,3,4-Thiadiazoles

The three-dimensional quantitative structure-activity relationships of a series of 5-[ 1-aryl-1,4-dihydro-6-methylpyridazin-4-one-3-yl]-2-arylamino-1,3,4-thiadiazoles, related to the fungicidal activity, were studied using the comparative molecular field analysis （CoMFA）. The results show that the contributions of steric and electrostatic fields to the activity are 0.505 and 0.495, respectively. The cross-validated q^2 and the correlation coefficient r^2 for the model established by the study are 0.769 and 0.938, respectively, with the F value of 60.996, and the standard deviation s of 0.074. These values indicate that the model is significant and has good predictability. The analysis results are in good agreement well with the study of 2D-QSAR, and offered important structural insights into designing highly active compounds prior to synthesis.     

Binding conformations,QSAR, and molecular design of Alkene‐3‐quinolinecarbonitriles as Src inhibitors

A theoretical study on binding orientations and quantitative structure–activity relationship (QSAR) of a novel series of alkene‐3‐quinolinecarbonitriles acting as Src inhibitors has been carried out by using the docking study and three‐dimensional QSAR (3D‐QSAR) analyses. The appropriate binding orientations and conformations of these compounds interacting with Src kinase were revealed by the docking studies, and the established 3D‐QSAR models show significant statistical quality and satisfactory predictive ability, with high R² values and q² values: comparative molecular field analysis (CoMFA) model (q² = 0.748, R² = 0.972), comparative molecular similarity indices analysis (CoMSIA) model (q² = 0.731, R² = 0.987). The systemic external validation indicated that both CoMFA and CoMSIA models possessed high predictive powers with $ R{^2}_{\!\!\!\rm pred} $ values of 0.818 and 0.892, $ {r^2}_{\!\!\!\rm m} $ values of 0.879 and 0.886, $ {r^2}_{\!\!\!\rm m(LOO)} $ values of 0.874 and 0.874, $ r^2_{\rm m(overall)} $ values of 0.879 and 0.885, respectively. Several key structural features of the compounds responsible for inhibitory activity were discussed in detail. Based on these structural factors, eight new compounds with quite higher predicted Src‐inhibitory activities have been designed and presented. We hope these theoretical results can offer some valuable references for the pharmaceutical molecular design as well as the action mechanism analysis. © 2012 Wiley Periodicals, Inc.     

Controlled Synthesis of 3D Flower‐like Ni2P Composed of Mesoporous Nanoplates for Overall Water Splitting

Developing efficient non‐noble metal and earth‐abundant electrocatalysts with tunable microstructures for overall water splitting is critical to promote clean energy technologies for a hydrogen economy. Herein, novel three‐dimensional (3D) flower‐like Ni₂P composed of mesoporous nanoplates with controllable morphology and high surface area was prepared by a hydrothermal method and low‐temperature phosphidation as efficient electrocatalysts for the hydrogen evolution reaction (HER) and oxygen evolution reaction (OER). Compared with the urchin‐like Ni_x P_y , the 3D flower‐like Ni₂P with a diameter of 5 μm presented an efficient and stable catalytic performance in 0.5 m H₂SO₄, with a small Tafel slope of 79 mV dec⁻¹ and an overpotential of about 240 mV at a current density of 10 mA cm⁻² with a mass loading density of 0.283 mg cm⁻². In addition, the catalyst also exhibited a remarkable performance for the OER in 1.0 m KOH electrolyte, with an overpotential of 320 mV to reach a current density of 10 mA cm⁻² and a small Tafel slope of 72 mV dec⁻¹. The excellent catalytic performance of the as‐prepared Ni₂P may be ascribed to its novel 3D morphology with unique mesoporous structure.     

Thermochemistry and Kinetics of the Thermal Decomposition of 1‐Chlorohexane

A theoretical kinetic study of the thermal decomposition of 1‐chlorohexane in gas phase between 600 and 1000 K was performed. Transition‐state theory and unimolecular reaction rate theory were combined with molecular information provided by quantum chemical calculations. Particularly, the B3LYP, BMK, M05–2X, and M06–2X formulations of the density functional theory (DFT) and the high‐level ab initio methods G3B3 and G4 were employed. The possible reaction channels for the thermal decomposition of 1‐chlorohexane were investigated, and the reaction takes place through the elimination of HCl with the formation of 1‐hexene. The derived high‐pressure limit rate coefficients are k _∞(600–1000 K) = (8 ± 5) × 10¹³ exp[‐((56.7 ± 0.4) kcal mol⁻¹/RT )] s⁻¹. The pressure effect over the reaction was analyzed from the calculation of the low‐pressure limit rate coefficients and the falloff curves. In addition, the standard enthalpies of formation at 298 K of −46.9 ± 1.5 kcal mol⁻¹ for 1‐chlorohexane and 5.8 ± 1.5 kcal mol⁻¹ for C₆H₁₃ radical were derived from isodesmic and isogiric reactions at high levels of theory.