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.     

Development of non-peptide ACE inhibitors as novel and potent cardiovascular therapeutics: An in silico modelling approach

Angiotensin-converting enzyme (ACE) inhibitors have been acknowledged as first-line agents for the treatment of hypertension and a variety of cardiovascular disorders. In this context, quantitative structure–activity relationship (QSAR) models for a series of non-peptide compounds as ACE inhibitors are developed based on Simplified Molecular Input-Line Entry System (SMILES) notation and local graph invariants. Three random splits into the training and test sets are used. The Monte Carlo method is applied for model development. Molecular docking studies are used for the final assessment of the developed QSAR model and the design of novel inhibitors. The statistical quality of the developed model is good. Molecular fragments responsible for the increase/decrease of the studied activity are calculated. The computer-aided design of new compounds, as potential ACE inhibitors, is presented. The predictive potential of the applied approach is tested, and the robustness of the model is proven using different methods. The results obtained from molecular docking studies are in excellent correlation with the results from QSAR studies. The presented study may be useful in the search for novel cardiovascular therapeutics based on ACE inhibition.     

Integration of genetic virtual screening patterns and latent multivariate modeling techniques for QSAR optimization based on combinations and/or interactions between peptides and proteins

Both the concept and the model of snug quantitative structure-activity relationship (QSAR) were pro-posed and developed for molecular design through constructing QSAR based on some known mode of receptor/ligand interactions. Many disadvantages of traditional models can be avoided by using the proposed method because the traditional models only determined upon molecular structural features in sample sets themselves. A genetic virtual screening of peptide/protein combinations (GVSPPC) is proposed for the first time by utilizing this idea to examine peptide/protein affinity activities. A genetic algorithm (GA) was developed for screening combinative targets with an interaction mode for virtual receptors. GVSPPC succeeds in disposing difficulties in rational QSAR,in order to search for the ligand/receptor interactions on conditions of unknown structures. Some bioactive oligo-/poly-peptide systems covering 58 angiotensin converting enzyme (ACE) inhibitors and 18 double site mutation residues in camel antibody protein cAb-Lys3 were investigated by GVSPPC with satisfactory results (R 2 cu>0.91,Q 2 cv > 0.86,ERMS=0.19-0.95),respectively,which demonstrates that GVSPPC is more inter-pretable in the ligand-receptor interaction than the traditional QSAR method.     

A combined molecular modeling study on gelatinases and their potent inhibitors

Zinc‐dependent matrix metalloproteinase (MMP) family is considered to be an attractive target because of its important role in many physiological and pathological processes. In the present work, a molecular modeling study combining protein‐, ligand‐ and complex‐based computational methods was performed to analyze a new series of β‐N‐biaryl ether sulfonamide hydroxamates as potent inhibitors of gelatinase A (MMP‐2) and gelatinase B (MMP‐9). Firstly, the similarities and differences between the binding sites of MMP‐2 and MMP‐9 were analyzed through sequence alignment and structural superimposition. Secondly, in order to extract structural features influencing the activities of these inhibitors, quantitative structure‐activity relationship (QSAR) models using genetic algorithm‐multiple linear regression (GA‐MLR), comparative molecular field (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were developed. The proposed QSAR models could give good predictive ability for the studied inhibitors. Thirdly, docking study was employed to further explore the binding mode between the ligand and protein. The results from all the above analyses could provide the information about the similarities and differences of the binding mode between the MMP‐2, MMP‐9 and their potent inhibitors. The obtained results can provide very useful information for the design of new potential inhibitors. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

酪氨酸激酶抑制剂的三维定量构效关系研究

利用柔性原子受体模型(FLARM)方法对一系列的异黄酮和喹诺酮衍生物表皮生长因子受体酪氨酸激酶抑制剂进行了三维定量构效关系研究，得到了合理的构效关系模型.FLARM方法的计算结果还给出了虚拟的受体模型，该模型说明了抑制剂与受体之间可能的相互作用.由该虚拟受体模型得到的受体-配体相互作用与Novartis药效团模型比较类似.     

Nonfitting protein–ligand interaction scoring function based on first‐principles theoretical chemistry methods: Development and application on kinase inhibitors

Targeted therapy is currently a hot topic in the fields of cancer research and drug design. An important requirement for this approach is the development of potent and selective inhibitors for the identified target protein. However, current ways to estimate inhibitor efficacy rely on empirical protein–ligand interaction scoring functions which, suffering from their heavy parameterizations, often lead to a low accuracy. In this work, we develop a nonfitting scoring function, which consists of three terms: (1) gas‐phase protein‐ligand binding enthalpy obtained by the eXtended ONIOM hybrid method based on an integration of density functional theory (DFT) methods (XYG3 and ωB97X‐D) and the semiempirical PM6 method, (2) solvation free energy based on DFT‐SMD solvation model, and (3) entropy effect estimated by using DFT frequency analysis. The new scoring function is tested on a cyclin‐dependent kinase 2 (CDK2) inhibitor database including 76 CDK2 protein inhibitors and a p21‐activated kinase 1 (PAK1) inhibitor database including 20 organometallic PAK1 protein inhibitors. From the results, good correlations are found between the calculated scores and the experimental inhibitor efficacies with the square of correlation coefficient R² of 0.76–0.88. This suggests a good predictive power of this scoring function. To the best of our knowledge, this is the first high level theory‐based nonfitting scoring function with such a good level of performance. This scoring function is recommended to be used in the final screening of lead structure derivatives. © 2013 Wiley Periodicals, Inc.     

An alignment‐free methodology for modelling field‐based 3D‐structure activity relationships using inductive logic programming

Traditional 3D‐quantitative structure–activity relationship (QSAR)/structure–activity relationship (SAR) methodologies are sensitive to the quality of an alignment step which is required to make molecular structures comparable. Even though many methods have been proposed to solve this problem, they often result in a loss of model interpretability. The requirement of alignment is a restriction imposed by traditional regression methods due to their failure to represent relations between data objects directly. Inductive logic programming (ILP) is a class of machine‐learning methods able to describe relational data directly. We propose a new methodology which is aimed at using the richness in molecular interaction fields (MIFs) without being restricted by any alignment procedure. A set of MIFs is computed and further compressed by finding their minima corresponding to the sites of strongest interaction between a molecule and the applied test probe. ILP uses these minima to build easily interpretable rules about activity expressed as pharmacophore rules in the powerful language of first‐order logic. We use a set of previously published inhibitors of factor Xa of the benzamidine family to discuss the problems, requirements and advantages of the new methodology. Copyright © 2007 John Wiley & Sons, Ltd.     

Searching Inhibitors of Adenosine Kinase by Simulation Methods

Searching new inhibitors of adenosine kinase （AK） is still drawing attention of experimental scientists. A better and solid model is here proposed by means of simulation methods from different ways, the direct analysis of receptor itself, the conventional 3D-QSAR methods and the integration of docking method and the conventional QSAR analysis.     

Modeling ligand-receptor interaction for some MHC class II HLA-DR4 peptide mimetic inhibitors using several molecular docking and 3D QSAR techniques

The ligand-receptor interaction between some peptidomimetic inhibitors and a class II MHC peptide presenting molecule, the HLA-DR4 receptor, was modeled using some three-dimensional (3D) quantitative structure-activity relationship (QSAR) methods such as the Comparative Molecular Field Analysis (CoMFA), Comparative Molecular Similarity Indices Analysis (CoMSIA), and a pharmacophore building method, the Catalyst program. The structures of these peptidomimetic inhibitors were generated theoretically, and the conformations used in the 3D QSAR studies were defined by docking them into the known structure of HLA-DR4 receptor through the GOLD, GLIDE Rigidly, GLIDE Flexible, and Xscore programs. Some of the parameters used in these docking programs were selected by docking an X-ray ligand into the receptor and comparing the root-means-square difference (RMSD) computed between the coordinates of the X-ray and docked structure. However, the goodness of a docking result for docking a series of peptidomimetic inhibitors into the HLA-DR4 receptor was judged by comparing the Spearman's rank correlation coefficient computed between each docking result and the activity data taken from the literature. The best CoMFA and CoMSIA models were constructed using the aligned structures of the best docking result. The CoMSIA was conducted in a stepwise manner to identify some important molecular features that were further employed in a pharmacophore building process by the Catalyst program. It was found that most inhibitors of the training set were accurately predicted by the best pharmacophore model, the Hypo1 hypothesis constructed. The deviation or conflict found between the actual and predicted activities of some inhibitors of both the training and the test sets were also investigated by mapping the Hypo1 hypothesis onto the corresponding structures of the inhibitors.