Multiple ligand simultaneous docking: Orchestrated dancing of ligands in binding sites of protein

Present docking methodologies simulate only one single ligand at a time during docking process. In reality, the molecular recognition process always involves multiple molecular species. Typical protein–ligand interactions are, for example, substrate and cofactor in catalytic cycle; metal ion coordination together with ligand(s); and ligand binding with water molecules. To simulate the real molecular binding processes, we propose a novel multiple ligand simultaneous docking (MLSD) strategy, which can deal with all the above processes, vastly improving docking sampling and binding free energy scoring. The work also compares two search strategies: Lamarckian genetic algorithm and particle swarm optimization, which have respective advantages depending on the specific systems. The methodology proves robust through systematic testing against several diverse model systems: E. coli purine nucleoside phosphorylase (PNP) complex with two substrates, SHP2NSH2 complex with two peptides and Bcl‐xL complex with ABT‐737 fragments. In all cases, the final correct docking poses and relative binding free energies were obtained. In PNP case, the simulations also capture the binding intermediates and reveal the binding dynamics during the recognition processes, which are consistent with the proposed enzymatic mechanism. In the other two cases, conventional single‐ligand docking fails due to energetic and dynamic coupling among ligands, whereas MLSD results in the correct binding modes. These three cases also represent potential applications in the areas of exploring enzymatic mechanism, interpreting noisy X‐ray crystallographic maps, and aiding fragment‐based drug design, respectively. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Protein‐protein docking by shape‐complementarity and property matching

We present a computational approach to protein‐protein docking based on surface shape complementarity (“ProBinder”). Within this docking approach, we implemented a new surface decomposition method that considers local shape features on the protein surface. This new surface shape decomposition results in a deterministic representation of curvature features on the protein surface, such as “knobs,” “holes,” and “flats” together with their point normals. For the actual docking procedure, we used geometric hashing, which allows for the rapid, translation‐, and rotation‐free comparison of point coordinates. Candidate solutions were scored based on knowledge‐based potentials and steric criteria. The potentials included electrostatic complementarity, desolvation energy, amino acid contact preferences, and a van‐der‐Waals potential. We applied ProBinder to a diverse test set of 68 bound and 30 unbound test cases compiled from the Dockground database. Sixty‐four percent of the protein‐protein test complexes were ranked with an root mean square deviation (RMSD) < 5 Å to the target solution among the top 10 predictions for the bound data set. In 82% of the unbound samples, docking poses were ranked within the top ten solutions with an RMSD < 10 Å to the target solution. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Mechanism of the decrease in catalytic activity of human cytochrome P450 2C9 polymorphic variants investigated by computational analysis

Cytochrome P450 (CYP) is deeply involved in the metabolism of chemicals including pharmaceuticals. Therefore, polymorphisms of this enzyme have been widely studied to avoid unfavorable side effects of drugs in chemotherapy. In this work, we performed computational analysis of the mechanism of the decrease in enzymatic activity for three typical polymorphisms in CYP 2C9 species: *2, *3, and *5. Based on the equilibrated structure obtained by molecular dynamics simulation, the volume of the binding pocket and the fluctuation of amino residues responsible for substrate holding were compared between the wild type and the three variants. Further docking simulation was carried out to evaluate the appropriateness of the binding pocket to accommodate substrate chemicals. Every polymorphic variant was suggested to be inferior to the wild type in enzymatic ability from the structural viewpoint. F‐G helices were obviously displaced outward in CYP2C9*2. Expansion of the binding pocket, especially the space near F′ helix, was remarkable in CYP2C9*3. Disappearance of the hydrogen bond between K helix and β4 loop was observed in CYP2C9*5. The reduction of catalytic activity of those variants can be explained from the deformation of the binding pocket and the consequent change in binding mode of substrate chemicals. The computational approach is effective for predicting the enzymatic activity of polymorphic variants of CYP. This prediction will be helpful for advanced drug design because calculations forecast unexpected change in drug efficacy for individuals. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Drug efficiency indices for improvement of molecular docking scoring functions

A dataset of protein‐drug complexes with experimental binding energy and crystal structure were analyzed and the performance of different docking engines and scoring functions (as well as components of these) for predicting the free energy of binding and several ligand efficiency indices were compared. The aim was not to evaluate the best docking method, but to determine the effect of different efficiency indices on the experimental and predicted free energy. Some ligand efficiency indices, such as ΔG/W (Wiener index), ΔG/NoC (number of carbons), and ΔG/P (partition coefficient), improve the correlation between experimental and calculated values. This effect was shown to be valid across the different scoring functions and docking programs. It also removes the common bias of scoring functions in favor of larger ligands. For all scoring functions, the efficiency indices effectively normalize the free energy derived indices, to give values closer to experiment. Compound collection filtering can be done prior or after docking, using pharmacokinetic as well as pharmacodynamic profiles. Achieving these better correlations with experiment can improve the ability of docking scoring functions to predict active molecules in virtual screening. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

含1,3,4-噁二唑的喹啉衍生物的合成、抗AChE活性及分子对接研究

为研发新型乙酰胆碱酯酶抑制剂,本文通过碘介导环化反应合成了一系列含1,3,4-噁二唑的喹啉衍生物,并通过熔点和核磁共振确定了其结构。在1 μmol/mL浓度下,对这些衍生物抗乙酰胆碱酯酶活性进行了初步的生物测定。结果表明,化合物4b和4j具有中等的抑制活性,抑制率分别为79.18%和78.46%。初步的构效关系表明,在目标化合物的1,3,4-噁二唑环上引入杂环可以其提高其活性。分子对接研究表明,化合物4b和4j和乙酰胆碱酯酶的催化活性中心部位显著结合。     

同源模建和分子对接研究HDAC1/HDAC8的选择性

组蛋白去乙酰化酶(HDACs)是近年来治疗肿瘤的重要靶标之一.由于HDACs包含多种亚型,且各亚型的生理功能存在一定的差异,其选择性抑制剂的开发已成为当前的研发热点.我们通过同源模建的HDAC1结构,与已有的HDAC8晶体结构的活性位点进行比较分析,探讨了对两者选择性有重要影响的残基,为基于受体的选择性抑制剂研究提供重要信息.同时选择了52个HDAC抑制剂,分别建立了HDAC1、HDAC8的活性值与对接打分值的线性回归模型.所建的HDAC1和HDAC8的线性构效关系模型的非交叉验证系数R2分别为0.82和0.80,表明具有一定的统计学意义.利用所建模型对已设计合成的化合物进行了预测,预测结果对HDAC1、HDAC8选择性抑制剂的优化改造提供了一定的指导意义.     

基于单孔洞近似的高纯铝部分层裂实验的数值模拟研究

基于单孔洞近似,对不同撞击速度下高纯铝的部分层裂实验进行了数值模拟研究,讨论了微孔洞长大对波传播的影响及其在自由面速度波剖面上的表现. 通过分析微孔洞周围的应力场变化,认识到实测自由面速度波剖面出现"回跳"特征并不能说明材料发生完全层裂,其直接原因是样品内部微孔洞长大所引起的局部卸载效应. 将计算得到的自由面速度波剖面和微孔洞相对体积与实验结果进行了对比分析,发现两者均符合很好,表明采用单孔洞增长来近似描述部分层裂样品中随机损伤发展及其对波传播的影响是可行的. 关键词： 层裂 孔洞增长 自由面速度波剖面 微孔洞相对体积