酪氨酸蛋白磷酸酯酶1B抑制剂的分子对接和三维定量构效关系研究

用一种柔性分子对接方法(FlexX)将12个2-草酰胺苯甲酸类抑制剂和酪氨酸蛋白磷酸酯酶(PTP1B)活性口袋进行分子对接,对接程序预测的抑制剂和酶之间的相互作用能与抑制活性之间有很好的相关性(非线性相关系数R²达0.859),这说明对接结果可以比较准确地预测抑制剂和PTP1B之间的结合模式.然后,将33个同类抑制剂的骨架叠合在分子对接预测的结合构象上,用比较分子力场分析方法(CoMFA)对其进行三维定量活性构效关系研究,得到的CoMFA模型具有很好的统计相关性(交互验证回归系数q²为0.650),并可以准确地预测测试集6个化合物的活性(平均标准偏差为0.177).同时,由CoMFA模型得出的抑制剂改造信息与用FlexX预测的结合模式是一致的,进一步证明我们预测的结合模式是正确的.为研究这类抑制剂和PTP1B的结合模式及对抑制剂进行结构改造提供了信息.     

Evaluation of docking/scoring approaches: A comparative study based on MMP3 inhibitors

An increasing number of docking/scoring programs are available that use different sampling and scoring algorithms. A reliable scoring function is the crucial element of such approaches. Comparative studies are needed to evaluate their current capabilities. DOCK4 with force field and PMF scoring as well as FlexX were used to evaluate the predictive power of these docking/scoring approaches to identify the correct binding mode of 61 MMP-3 inhibitors in a crystal structure of stromelysin and also to rank them according to their different binding affinities. It was found that DOCK4/PMF scoring performs significantly better than FlexX and DOCK4/FF in both ranking ligands and predicting their binding modes. Most notably, DOCK4/PMF was the only scoring/docking approach that found a significant correlation between binding affinity and predicted score of the docked inhibitors. However, comparing only those cases where the correct binding mode was identified (scoring highest among sampled poses), FlexX showed the best `fine tuning' (lowest rmsd) in predicted binding modes. The results suggest that not so much the sampling procedure but rather the scoring function is the crucial element of a docking program.     

Molecular docking and competitive binding study discovered different binding modes of microsomal prostaglandin E synthase-1 inhibitors

Microsomal prostaglandin E synthase-1 (mPGES-1) is a newly recognized therapeutic target for the treatment of inflammation, pain, cancer, atherosclerosis, and stroke. Many mPGES-1 inhibitors have been discovered. However, as the structure of the binding site is not well-characterized, none of these inhibitors was designed based on the mPGES-1 structure, and their inhibition mechanism remains to be fully disclosed. Recently, we built a new structural model of mPGES-1 which was well supported by experimental data. Based on this model, molecular docking and competition experiments were used to investigate the binding modes of four representive mPGES-1 inhibitors. As the inhibitor binding sites predicted by docking overlapped with both the substrate and the cofactor binding sites, mPGES-1 inhibitors might act as dual-site inhibitors. This inhibitory mechanism was further verified by inhibitor-cofactor and inhibitor-substrate competition experiments. To investigate the potency-binding site relationships of mPGES-1 inhibitors, we also carried out molecular docking studies for another series of compounds. The docking results correlated well with the different inhibitory effects observed experimentally. Our data revealed that mPGES-1 inhibitors could bind to the substrate and the cofactor binding sites simultaneously, and this dual-site binding mode improved their potency. Future rational design and optimization of mPGES-1 inhibitors can be carried out based on this binding mechanism.     

Molecular interactions of nitrogen-containing bisphosphonates within farnesyl diphosphate synthase

Bisphosphonates, known for their effectiveness in the treatment of osteoporosis, inhibit bone resorption via mechanisms that involve binding to bone mineral and cellular effects on osteoclasts. The major molecular target of nitrogen-containing bisphosphonates (N-BPs) in osteoclasts is farnesyl diphosphate synthase (FPPS). N-BPs likely inhibit this enzyme by mimicking one or more of the natural isoprenoid lipid substrates (GPP/DMAPP and IPP) but the mode of inhibition is not established. The active site of FPPS comprises a subsite for each substrate. Kinetic studies with recombinant human FPPS indicate that both potent (risedronate) and weak (NE-58051) enzyme inhibitors compete with GPP for binding to FPPS, however, binding to this site does not completely explain the difference in potency of the two inhibitors, suggesting that a second binding site may also be a target of bisphosphonate inhibition. Using the docking software suite Autodock, we explored a dual inhibitor binding mode for recombinant human FPPS. Experimental support for dual binding is suggested by Dixon plots for the inhibitors. N-BPs may inhibit by binding to both the GPP and a second site with differences in potency at least partly arising from inhibition at the second site.     

Active site similarity between human and <Emphasis Type="Italic">Plasmodium falciparum</Emphasis> phosphodiesterases: considerations for antimalarial drug design

Abstract  

The similarity between Plasmodium falciparum phosphodiesterase enzymes (PfPDEs) and their human counterparts have been examined and human PDE9A was found to be a suitable template for the construction of homology models for each of the four PfPDE isoforms. In contrast, the architecture of the active sites of each model was most similar to human PDE1. Molecular docking was able to model cyclic guanosine monophosphate (cGMP) substrate binding in each case but a docking mode supporting cyclic adenosine monophosphate (cAMP) binding could not be found. Anticipating the potential of PfPDE inhibitors as anti-malarial drugs, a range of reported PDE inhibitors including zaprinast and sildenafil were docked into the model of PfPDEα. The results were consistent with their reported biological activities, and the potential of PDE1/9 inhibitor analogues was also supported by docking.     

HIV-1整合酶链转移抑制剂的3D-QSAR、分子对接和分子动力学模拟研究

为了获得高活性、结构新颖的整合酶链转移(INST)抑制剂,本文采用Co MFA和Co MSIA两种方法对32个萘啶类INST抑制剂进行了三维定量构效关系研究,并建立了相关模型,其交叉验证系数分别为q~2=0. 809和q~2=0. 816,拟合验证系数分别为r~2=0. 998和r~2=0. 981,表明所建立的模型是可靠的且具有一定的预测能力。利用分子对接探讨小分子化合物与INST蛋白的相互作用模式,结果表明,萘啶类化合物主要通过疏水作用和氢键作用与INSTIs蛋白结合。最后通过分子动力学模拟进一步验证对接结果发现,对接的结合模式与分子动力学模拟得到的结果是一致的。本研究获得的综合模型和推论可以为开发有效的HIV INSTIs提供重要的理论信息。     

Molecular modeling study of checkpoint kinase 1 inhibitors by multiple docking strategies and prime/MM-GBSA calculation

Developing chemicals that inhibit checkpoint kinase 1 (Chk1) is a promising adjuvant therapeutic to improve the efficacy and selectivity of DNA-targeting agents. Reliable prediction of binding-free energy and binding affinity of Chk1 inhibitors can provide a guide for rational drug design. In this study, multiple docking strategies and Prime/Molecular Mechanics Generalized Born Surface Area (Prime/MM-GBSA) calculation were applied to predict the binding mode and free energy for a series of benzoisoquinolinones as Chk1 inhibitors. Reliable docking results were obtained using induced-fit docking and quantum mechanics/molecular mechanics (QM/MM) docking, which showed superior performance on both ligand binding pose and docking score accuracy to the rigid-receptor docking. Then, the Prime/MM-GBSA method based on the docking complex was used to predict the binding-free energy. The combined use of QM/MM docking and Prime/MM-GBSA method could give a high correlation between the predicted binding-free energy and experimentally determined pIC(50) . The molecular docking combined with Prime/MM-GBSA simulation can not only be used to rapidly and accurately predict the binding-free energy of novel Chk1 inhibitors but also provide a novel strategy for lead discovery and optimization targeting Chk1.     

HDAC对接研究:苯甲酰胺类抑制剂结合方式推测

通过计算机模拟的对接过程研究，发现了MS-275— 一种苯甲酰胺类的组蛋白去乙酰酶（HDAC）抑制剂与酶的可能的全新结合方式.这种结合方式与已经阐明的组蛋白去乙酰酶类似蛋白（HDLP）与曲古柳菌素A(trichostatin A, TSA)和suberoylanilide hydroxamic acid(SAHA)形成的复合物晶体结构中配体与酶的作用方式完全不同.从对接结果看，MS-275的作用靶点在酶活性口袋的最狭窄部位，而不是直接作用于锌离子.这似乎能够解释MS-275的低毒性特点，并且为设计和筛选全新的HDAC抑制剂提供了新思路.     

Identifying the binding mode of a molecular scaffold

We describe a method for docking of a scaffold-based series and present its advantages over docking of individual ligands, for determining the binding mode of a molecular scaffold in a binding site. The method has been applied to eight different scaffolds of protein kinase inhibitors (PKI). A single analog of each of these eight scaffolds was previously crystallized with different protein kinases. We have used FlexX to dock a set of molecules that share the same scaffold, rather than docking a single molecule. The main mode of binding is determined by the mode of binding of the largest cluster among the docked molecules that share a scaffold. Clustering is based on our 'nearest single neighbor' method [J. Chem. Inf. Comput. Sci., 43 (2003) 208-217]. Additional criteria are applied in those cases in which more than one significant binding mode is found. Using the proposed method, most of the crystallographic binding modes of these scaffolds were reconstructed. Alternative modes, that have not been detected yet by experiments, could also be identified. The method was applied to predict the binding mode of an additional molecular scaffold that was not yet reported and the predicted binding mode has been found to be very similar to experimental results for a closely related scaffold. We suggest that this approach be used as a virtual screening tool for scaffold-based design processes.     

Limiting the Number of Potential Binding Modes by Introducing Symmetry into Ligands: Structure‐Based Design of Inhibitors for Trypsin‐Like Serine Proteases

In the absence of X‐ray data, the exploration of compound binding modes continues to be a challenging task. For structure‐based design, specific features of active sites in different targets play a major role in rationalizing ligand binding characteristics. For example, dibasic compounds have been reported as potent inhibitors of various trypsin‐like serine proteases, the active sites of which contain several binding pockets that can be targeted by cationic moieties. This results in several possible orientations within the active site, complicating the binding mode prediction of such compounds by docking tools. Therefore, we introduced symmetry in bi‐ and tribasic compounds to reduce conformational space in docking calculations and to simplify binding mode selection by limiting the number of possible pocket occupations. Asymmetric bisbenzamidines were used as starting points for a multistage and structure‐guided optimization. A series of 24 final compounds with either two or three benzamidine substructures was ultimately synthesized and evaluated as inhibitors of five serine proteases, leading to potent symmetric inhibitors for the pharmaceutical drug targets matriptase, matriptase‐2, thrombin and factor Xa. This study underlines the relevance of ligand symmetry for chemical biology.     

人类蛋白激酶组的多靶点分子对接系统

蛋白激酶在信号转导、基因转录和蛋白翻译等生物过程起关键性作用,因而与大量人类疾病密切相关。所以,蛋白激酶的抑制剂筛选是抗肿瘤药物开发的热点,正在向基于全激酶组的高通量多靶点筛选模式发展。为了降低大规模实验筛选的成本,提高成功率,本文构建人类蛋白激酶组的多靶点分子对接系统,对抑制剂-激酶组的相互作用进行预测。我们首先利用同源模建方法,对人类激酶组约500个激酶变异体的催化域进行结构建模;接着以催化域结构模型为受体,用已知激酶抑制剂进行分子对接,对抑制剂与各激酶变异体的结合亲和力进行了定量计算。结果显示,本文所建立的多靶点分子对接系统可以准确预测抑制剂与激酶变异体的相互作用,结合自由能的计算值与实验值有很强的相关性。所以,该分子对接系统可用于多靶点激酶抑制剂的计算筛选,为激酶抑制剂开发与抗肿瘤药物设计提供理论依据。     

Microscopic binding of butyrylcholinesterase with quinazolinimine derivatives and the structure–activity correlation

Butyrylcholinesterase (BChE) is not only an important protein for development of anti-cocaine medication but also an established drug target to develop new treatment for Alzheimer’s disease (AD). The molecular basis of interaction of a new series of quinazolinimine derivatives as BChE inhibitors has been studied by molecular docking and molecular dynamics (MD) simulations. The molecular docking and MD simulations revealed that all of these inhibitors bind with BChE in similar binding mode. Based on the similar binding mode, we have carried out three-dimensional quantitative structure–activity relationship (3D-QSAR) studies on these inhibitors using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA), to understand the structure–activity correlation of this series of inhibitors and to develop predictive models that could be used in the design of new inhibitors of BChE. The study has resulted in satisfactory 3D-QSAR models. We have also developed ligand-based 3D-QSAR models. The contour maps obtained from the 3D-QSAR models in combination with the simulated binding structures help to better interpret the structure–activity relationship and is consistent with available experimental activity data. The satisfactory 3D-QSAR models strongly suggest that the determined BChE-inhibitor binding modes are reasonable. The identified binding modes and developed 3D-QSAR models for these BChE inhibitors are expected to be valuable for rational design of new BChE inhibitors that may be valuable in the treatment of Alzheimer’s disease.     

A selectivity study of benzenesulfonamide derivatives on human carbonic anhydrase II/IX by 3D-QSAR,Molecular Docking and Molecular Dynamics Simulation

Nowadays, different approaches have been pursued with the intent to develop sulfonamide-like carbonic anhydrase inhibitors that possess better selectivity profiles toward the different human isoforms of the enzyme. Here, we used conventional 3D-QSAR methods, including comparative molecular field analysis (CoMFA), comparative molecular similarity indices analysis (CoMSIA), and Topomer CoMFA, to construct three-dimensional quantitative structure-activity relationship (3D-QSAR) models for benzenesulfonamide derivatives as human carbonic anhydrase (hCA) II/IX inhibitors. The theoretical models had good reliability (R²>0.75) and predictability (Q²>0.55), and the contour maps could graphically present the contributions of the force fields for activity and identify the structural divergence between human carbonic anhydrase II inhibitors and human carbonic anhydrase IX inhibitors. Consequently, we explored the selectivity of inhibitor for human carbonic anhydrase II and IX through molecular docking, and the difference of activity coincides with the potential binding mode well. According to the results of the predicted values and the molecule docking, we found that the inhibitors published in the literature had stronger inhibition on the hCA IX; based on the theoretical models, we designed seven new compounds with good potential activity and reasonably good ADMET profile, which could selectively inhibit hCA IX. Molecular Dynamics Simulation showed that newly-designed compound D7 had good selectivity on hCA IX. The findings from 3D-QSAR and docking studies maybe helpful in the rational drug design of isoform-selective inhibitors.     

Molecular dynamics simulation, free energy calculation and structure-based 3D-QSAR studies of B-RAF kinase inhibitors

(V600E)B-RAF kinase is the most frequent onco-genic protein kinase mutation in melanoma and is a promising target to treat malignant melanoma. In this work, a molecular modeling study combining QM-polarized ligand docking, molecular dynamics, free energy calculation, and three-dimensional quantitative structure-activity relationships (3D-QSAR) was performed on a series of pyridoimidazolone compounds as the inhibitors of (V600E)B-RAF kinase to understand the binding mode between the inhibitors and (V600E)B-RAF kinase and the structural requirement for the inhibiting activity. 3D-QSAR models, including CoMFA and CoMSIA, were developed from the conformations obtained by QM-polarized ligand docking strategy. The obtained models have a good predictive ability in both internal and external validation. Furthermore, molecular dynamics simulation and free energy calculations were employed to determine the detailed binding process and to compare the binding mode of the inhibitors with different activities. The binding free energies calculated by MM/PBSA gave a good correlation with the experimental biological activity. The decomposition of free energies by MM/GBSA indicates the van der Waals interaction is the major driving force for the interaction between the inhibitors and (V600E)B-RAF kinase. The hydrogen bond interactions between the inhibitors with Glu501 and Asp594 of the (V600E)B-RAF kinase help to stabilize the DFG-out conformation. The results from this study can provide some insights into the development of novel potent (V600E)B-RAF kinase inhibitors.     

In silico based virtual screening and mixed mode QM/MM calculation identifies caffeine scaffold for designing potential inhibitors for tyrosyl tRNA synthetase of Mycobacterium tuberculosis

Aminoacyl tRNA synthetases are novel antibacterial drug target because of their important role in protein synthesis. In this study, we performed high throughput virtual screening of 205883 compounds from Asinex ligand database to identify potential specific inhibitors for Tyrosyl tRNA synthetase of Mycobacterium tuberculosis (MtbTyrRS). Compounds are ranked based on the glide extra precision docking score. It is noted that the top ranked compounds have caffeine scaffold. The top five caffeine analogs are further evaluated for other drug‐like properties. The binding energies of caffeine analogs are estimated using mixed mode quantum mechanics/molecular mechanics calculation. The results show that these caffeine analogs have good absorption, distribution, metabolism, and excretion properties and high binding affinity to the MtbTyrRS. This suggests that caffeine could be a new scaffold for designing inhibitors against Tyrosyl tRNA synthetase of M. tuberculosis. The top five caffeine analogs are also subjected to docking calculations with human cytosolic and mitochondrial Tyrosyl tRNA synthetases to ascertain their specificities toward MtbTyrRS. The comparative docking studies indicate that the top five caffeine analogs are specific for MtbTyrRS. © 2014 Wiley Periodicals, Inc.     

人类丝氨酸消旋酶的同源模建及其与多肽类抑制剂的分子对接

利用同源模建和分子动力学模拟方法构建了人类丝氨酸消旋酶(hSR)的三维结构, 并利用profile-3D和procheck方法评估了模型的可靠性. 在此基础上用分子对接程序(affinity)将多肽类抑制剂A和B分别与hSR进行对接, 获得了其复合物结构的理论模型. 通过配体与受体之间相互作用能和结构分析给出了此类抑制剂与hSR的具体结合方式, 明确了hSR与此类抑制剂结合时起重要作用的氨基酸残基, 为基于人类丝氨酸消旋酶三维结构的药物设计提供重要的参考信息.     

Identification of the quinolinedione inhibitor binding site in Cdc25 phosphatase B through docking and molecular dynamics simulations

Cdc25 phosphatase B, a potential target for cancer therapy, is inhibited by a series of quinones. The binding site and mode of quinone inhibitors to Cdc25B remains unclear, whereas this information is important for structure-based drug design. We investigated the potential binding site of NSC663284 [DA3003-1 or 6-chloro-7-(2-morpholin-4-yl-ethylamino)-quinoline-5, 8-dione] through docking and molecular dynamics simulations. Of the two main binding sites suggested by docking, the molecular dynamics simulations only support one site for stable binding of the inhibitor. Binding sites in and near the Cdc25B catalytic site that have been suggested previously do not lead to stable binding in 50 ns molecular dynamics (MD) simulations. In contrast, a shallow pocket between the C-terminal helix and the catalytic site provides a favourable binding site that shows high stability. Two similar binding modes featuring protein-inhibitor interactions involving Tyr428, Arg482, Thr547 and Ser549 are identified by clustering analysis of all stable MD trajectories. The relatively flexible C-terminal region of Cdc25B contributes to inhibitor binding. The binding mode of NSC663284, identified through MD simulation, likely prevents the binding of protein substrates to Cdc25B. The present results provide useful information for the design of quinone inhibitors and their mechanism of inhibition.     

Structure-based design and evaluation of novel N-phenyl-1H-indol-2-amine derivatives for fat mass and obesity-associated (FTO) protein inhibition

Fat mass and obesity-associated (FTO) protein contributes to non-syndromic human obesity which refers to excessive fat accumulation in human body and results in health risk. FTO protein has become a promising target for anti-obesity medicines as there is an immense need for the rational design of potent inhibitors to treat obesity. In our study, a new scaffold N-phenyl-1H-indol-2-amine was selected as a base for FTO protein inhibitors by applying scaffold hopping approach. Using this novel scaffold, different derivatives were designed by extending scaffold structure with potential functional groups. Molecular docking simulations were carried out by using two different docking algorithm implemented in CDOCKER (flexible docking) and AutoDock programs (rigid docking). Analyzing results of rigid and flexible docking, compound MU06 was selected based on different properties and predicted binding affinities for further analysis. Molecular dynamics simulation of FTO/MU06 complex was performed to characterize structure rationale and binding stability. Certainly, Arg96 and His231 residue of FTO protein showed stable interaction with inhibitor MU06 throughout the production dynamics phase. Three residues of FTO protein (Arg96, Asp233, and His231) were found common in making H-bond interactions with MU06 during molecular dynamics simulation and CDOCKER docking.     

Theoretical studies on the conformational change of adenosine kinase induced by inhibitors

Adenosine kinase (AK) is a two‐domain protein that catalyzes the phosphorylation of adenosine to adenosine monophosphate. Inhibitors of AK could increase adenosine to levels that activate nearby adenosine receptors and produce a wide variety of therapeutically beneficial activities. To get insight into the interaction mechanism between inhibitors and AK, we chose two kinds of novel inhibitors, alkynylpyrimidine inhibitor (APy) and aryl‐nucleoside inhibitor (AN), and used docking and molecular dynamics simulation methods to study the conformational changes of human AK on binding inhibitors. The calculation results revealed that both APy and AN could induce conformational changes of AK and stabilize AK at different semiopen conformations. On binding APy, the small lid‐domain rotated 14°, and the binding pocket rearranged after MD simulation. But in AK‐AN complex, the rotation of small domain is 22°, and the sugar ring of AN is mobile in the binding pocket. Further docking calculations on APy analogues indicate that the semiopen conformation could well explain the SAR of AK inhibitors. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Synthesis and activity of a new series of (Z)-3-phenyl-2-benzoylpropenoic acid derivatives as aldose reductase inhibitors

During the course of studies directed towards the discovery of novel aldose reductase inhibitors for the treatment of diabetic complications, we synthesized a series of new (Z)-3-phenyl-2-benzoylpropenoic acid derivatives and tested their in vitro inhibitory activities on rat lens aldose reductase. Of these compounds, (Z)-3-(3,4-dihydroxyphenyl)-2-(4-methylbenzoyl)propenoicacid (3k) was identified as the most potent inhibitor, with an IC50 of 0.49 microM. The theoretical binding mode of 3k was obtained by simulation of its docking into the active site of the human aldose reductase crystal structure.