Synthesis,In vitro Coagulation Activities and Molecular Docking Studies on Three L-Histidine Amide Derivatives

Three novel L-histidine amide derivatives were synthesized and the corresponding chemical structures were characterized by means of melting point analysis, IR, MS, ¹H NMR as well as ¹³C NMR. The coagulation acti- vities of the compounds were evaluated by an MOE(molecular operating environment) docking technique and coagulation test. The results obtained from molecular docking show that the interactions between the compounds and thrombin exhibit procoagulant activity in combination with an improved combinatory effect. Moreover, the results of in vitro coagulation tests show that the L-histidine amide derivatives feature coagulant activities in common coagulation pathways. Compared with the blank control group, the optimal shortening rates of compounds 1―3 were 39.08%(0.5 mmol/L), 22.94%(1.0 mmol/L) and 15.38%(0.0625 mmol/L), respectively.     

分子对接方法的应用与发展

分子对接方法加快了药物开发周期,具有快速、准确度高等优点.本文详述了分子对接方法的基本原理,及分子对接空间和能量的匹配要求和优化时的各种方法.综述了该方法在药物设计、药理分析和探测生命体系等方面的应用.     

Molecular Docking of 3-Methylindole-containing Drugs Binding into CYP3A4

Drugs SPD-304(6,7-dimethyl-3-{[methyl-(2-{methyl-[1-(3-trifluoromethyl-phenyl)-1H-indol-3-ylmethyl]- amino}-ethyl)-amino]-methyl}-chromen-4-one) and zafirlukast contain a common structural element of 3-substituted indole moiety which closely relates to a dehydrogenated reaction catalyzed by cytochrome P450s(CYPs). It was reported that the dehydrogenation can produce a reactive electrophilic intermediate which cause toxicities and inactivate CYPs. Drug L-745,870(3-{[4-(4-chlorophenyl)piperazin-1-yl]-methyl}-1H-pyrrolo- 2,3-β-pyridine) might have similar effect since it contains the same structural element. We used molecular docking approach combined with molecular dynamics(MD) simulation to model three-dimensional(3D) complex structures of SPD-304, zafirlukast and L-745,870 into CYP3A4, respectively. The results show that these three drugs can stably bind into the active site and the 3-methylene carbons of the drugs keep a reasonable reactive distance from the heme iron. The complex structure of SPD-304-CYP3A4 is in agreement with experimental data. For zafirlukast, the calculation results indicate that 3-methylene carbon might be the dehydrogenation reaction site. Docking model of L-745,870-CYP3A4 shows a potential possibility of L-745,870 dehydrogenated by CYP3A4 at 3-methylene carbon which is in agreement with experiment in vivo. In addition, residues in the phenylalanine cluster as well as S119 and R212 play a critical role in the ligands binding based on our calculations. The docking models could provide some clues to understand the metabolic mechanism of the drugs by CYP3A4.     

Molecular Modeling and Docking of Mannose-binding Lectin from Lycoris radiata

Lycoris radiata mannose-binding lectin（LRL） is a protein which binds mannose residues specifically. The maturation peptide and three mannose-binding domains（residues 49-57, 80-88 and 113--121） of LRL were identified by sequence analysis. The 3D structure of LRL constructed by homology modeling shaped a flstular triangular prism. Three flanks of the prism are mainly composed of β-sheets and each flank has a mannose-binding domain. According to the docking and dynamics simulation, the bindings of residues 49--57 and 80--88 with mannose are more stable than that of residues 113--121 with it. The key residues for binding mannose are Gin80, Asp82, Ash84 and Tyr88. The study preliminarily analyzed the interaction sites and mechanism of LRL with mannoses, which could be useful for the study on insect-resistance and related drug discovery of LRL.     

Protein structure predictions by parallel simulated annealing molecular dynamics using genetic crossover

We propose a conformational search method to find a global minimum energy structure for protein systems. The simulated annealing is a powerful method for local conformational search. On the other hand, the genetic crossover can search the global conformational space. Our method incorporates these attractive features of the simulated annealing and genetic crossover. In the previous works, we have been using the Monte Carlo algorithm for simulated annealing. In the present work, we use the molecular dynamics algorithm instead. To examine the effectiveness of our method, we compared our results with those of the normal simulated annealing molecular dynamics simulations by using an α-helical miniprotein. We used genetic two-point crossover here. The conformations, which have lower energy than those obtained from the conventional simulated annealing, were obtained.     

3C-like蛋白酶抑制剂的构效关系、分子对接和分子动力学

3C-like蛋白酶是中东呼吸综合征冠状病毒（MERS-CoV）等其它冠状病毒的繁殖过程中极为重要的蛋白酶。它已成为人类在抗冠状病毒领域中的研究热点。本文基于计算生物学方法对与MERS-CoV同属的蝙蝠冠状病毒HKU4（HKU4-CoV）的43个肽类3C-like蛋白酶抑制剂分子,建立三维定量构效关系（3D-QSAR）模型。在基于配体叠合的基础上,发现比较分子相似性指数分析法（CoMSIA）中的四个场组合（位阻场、静电场、氢键供体场与氢键受体场）为最优的模型（Q²=0.522,R_ncv²=0.996,R_pre²=0.904;Q²：交叉验证相关系数,R_ncv²：非交叉验证相关系数,R_pre²：验证集分子的预测值相关系数）,并借助该模型通过分子对接（docking）与分子动力学（MD）方法阐明了配受体结合作用。实验结果表明：（1）基于最优的CoMSIA模型基础上的三维等势图形象地说明了分子基团的位阻作用、静电作用、氢键供体与氢键受体作用对分子生物活性的影响;（2）分子对接研究结果显示了疏水性以及结晶水、氨基酸His166和Glu169在配体和受体结合过程中产生重要作用;（3）分子动力学模拟进一步验证了分子对接模型的可靠性,并发现了两个新的关键氨基酸Ser24与Gln192,它们与配体产生了两个较强的氢键。此外,根据这些结果,一些新的具有潜在抑制活性的肽类化合物作为3C-like蛋白酶抑制剂被获得。以上结果能够帮助深入了解3C-like蛋白酶与肽类抑制剂的作用机理,并且能够为今后的抗MERS-CoV药物设计提供有价值的参考。     

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

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

禽流感病毒HA蛋白与人受体的分子对接

血凝素(hemagglutinin,HA)是位于禽流感病毒表面的糖蛋白。在病毒感染过程中,HA与禽类宿主细胞表面受体结合,介导病毒膜与宿主核内体膜的融合,在传染过程中发挥关键作用。自然界中的禽流感病毒处于不断演化之中,其HA的禽受体结合位点常常发生氨基酸变异。因此,当HA变异体与人受体结合能力较强时,禽流感病毒往往会发生跨种传播而感染人。为预防禽流感的跨种传播,人们迫切需要发展大规模快速检测或预测HA变异体与人受体结合亲和力的方法,以评估各种新发禽流感病毒的跨种传播能力,提前筛选出有潜在危险的病毒株。针对此问题,本研究以H7N9亚型的HA蛋白H7为研究对象,发展了一种运用分子对接的计算方法,预测HA变异体与人受体的结合亲和力。该方法的计算结果表明,H7与人受体的结合亲和力普遍弱于有较强传染人能力的H1,说明H7N9亚型病毒的跨种传播能力普遍较弱;但是,计算分析也揭示,部分新发的H7N9毒株的HA有强的人受体结合亲和力,提示在自然演化过程中,H7N9病毒有可能演化出具有较强的感染人能力的新毒株,这与2013年禽流感疫情的实际发生情况相一致。因此,本文所发展的计算方法可用于快速预测新发禽流感病毒HA与人受体的结合亲和力,为新发禽流感病毒的跨种传播风险评估提供理论依据。     

Synthesis,Biological Evaluation and Molecular Docking of Novel Phenylpyrimidine Derivatives as Potential Anticancer Agents

Based on our previous researches, a novel phenylpyrimidine pharmacophore model was proposed and fifteen derivatives were synthesized and characterized by means of spectroscopy methods. The inhibitory effects of them were screened against HeLa cell line by virtue of MTT assay in vitro. The results indicate some of the phenylpyrimidine derivatives exhibit potent biological activities. Among them, compounds 6g and 6h exhibit the best activity at half maximal inhibitory concentrations of 1.5 and 2.8 μmol/L, respectively. These compounds also exhibit good activities against HepG2 cell line and MCF-7 cell line. FLT-3 kinase was screened as the most potent molecular target. Computational docking between compound 6g and FLT-3 was carried out to interpret the binding mode. The results show phenylpyrimidine derivatives have effective antitumor activities, which provides a base for further research of them as antitumor agents.     

分子对接和动力学模拟提高嗜热蛋白酶PhpI的活力

通过分子对接和动力学模拟对嗜热蛋白酶的分子进行改造, 确定蛋白酶PH1704(PhpI)定点突变残基, 并通过分子生物学实验进行验证. 突变体K43C的蛋白酶活力提高了5.8倍. 分子动力学模拟结果表明, 经过8 ns的动力学模拟后, K43C突变体二级结构由野生型的S2片层(F11-E12-D13)变成环状结构. E12和K43均是活性位点的重要残基, 这种变化将导致活性位点的柔性增强, 有利于催化反应的发生.     

多光谱法与分子对接模型研究西维来司钠与弹性蛋白酶的相互作用

采用多种光谱法及分子对接技术对西维来司钠（ONO-5046）与弹性蛋白酶的相互作用进行研究,利用荧光光谱法判断出ONO-5046的加入对弹性蛋白酶产生荧光猝灭作用,与弹性蛋白酶相互作用的猝灭类型为静态猝灭,ONO-5046与弹性蛋白酶以结合比为1∶1的比例构成复合物。经计算,ONO-5046与弹性蛋白酶体系的ΔH<0和ΔS<0,判断出氢键与范德华力为其主要结合作用力,由ΔG<0可知该反应可以自发进行。同步荧光光谱法、紫外-可见分光光谱法和圆二色光谱法探讨了ONO-5046的加入使弹性蛋白酶中氨基酸残基周围环境的疏水性降低,二级结构发生改变,α-螺旋结构的比例减少,无规则卷曲的比例增多;分子对接模拟结果表明：在与弹性蛋白酶相互作用过程中,ONO-5046分子骨架上的多个酰基结构、2个苯环以及较长的共轭结构能提供氢键与范德华力的作用位点,同时,端碳上的3个甲基结构有利于疏水作用,ONO-5046上的氧负离子与弹性蛋白酶的精氨酸存在盐桥作用。以上结果表明,ONO-5046与弹性蛋白酶的结合存在多种作用力,能形成稳定的复合物,从而抑制了弹性蛋白酶的活性。     

对羟基杏仁酸合成酶三维结构模建及其与底物的分子对接研究

以对羟基苯丙酮酸双氧化酶(HPPD)的晶体结构为模板, 利用同源模建方法构建了与其高度同源、底物相同但催化功能存在明显差别的对羟基杏仁酸合成酶(HMS)的三维结构, 并对模建结构的合理性进行了分析. 在模建结果的基础上, 对HPPD和HMS分别与底物羟苯基丙酮酸(HPP)进行分子对接计算, 比较了二者结合模式的异同, 为两种同源酶在催化方面差异性的合理阐释提供了一些有益的信息.     

Thermus thermophilus木糖异构酶与木糖醇的分子对接及模型分析

在1BXB结构基础上, 通过分子对接方法构建木糖异构酶与抑制剂木糖醇的复合物模型, 为合理设计解除木糖醇对木糖异构酶的抑制及进一步揭示木糖醇对该酶抑制机理提供参考.     

金属核酸酶及寡聚酰胺与双链DNA分子对接模式的理论研究

采用分子对接软件AutoDock分别研究了金属核酸酶、寡聚酰胺与DNA对接模式. 研究结果表明, 金属核酸酶、寡聚酰胺与DNA对接结果与实验数据吻合得较好. 在金属核酸酶与DNA的对接中, 最优构型都是金属核酸酶对接在DNA的小沟内. 随着中心金属上电荷的增加, 最大对接能呈现上升趋势. 尤其是配体体积较小的金属核酸酶, 最佳对接位点受中心金属上电荷影响较大, 即随着中心金属上电荷的从少到多, 最佳对接位点呈现多样化的趋势. 在柔性分子寡聚酰胺与双链DNA分子的对接中, 分子随DNA小沟的柔性变形是对接结合的关键. 其最佳的结合物位点也都是在DNA的小沟内, 其对接结果与晶体结构一致.     

双金属存在下整合酶和抑制剂5CITEP的分子对接研究

在HIV-1整合酶(IN)和5CITEP复合物晶体结构的基础上, 用分子对接程序(Affinity)将含有单Mg2+和双Mg2+ 的HIV-1 IN核心区与抑制剂5CITEP进行对接, 获得了能形成复合物结构的理论模型. 通过配体与受体之间的相互作用能和结构分析给出此种抑制剂的结合模式, 并与晶体结构进行比较, 揭示出引入的第二个Mg2+原子在整合过程中所起的重要作用. 前后相互作用能的变化趋势很明显, 配体和受体的作用模式比单Mg2+体系更加清晰. 由单Mg2+体系的4种作用方式改变到双Mg2+体系的两种作用方式, 相互作用能提高了将近40 kJ/mol. 为基于整合酶结构的药物设计提供了参考信息.     

Spectroscopic and Molecular Docking Approaches for Investigation of Interaction of Phellopterin with Human Serum Albumin

The mechanism of interaction between human serum albumin (HSA) and natural product phellopterin (PL) from Angelica dahurica was investigated by spectroscopic techniques with molecular docking under simulated physiological conditions. The experimental results showed that the fluorescence of HSA was regularly quenched by PL, and the quenching constants (K_SV) decreased with increasing temperature, which indicated that the quenching mechanism was a static quenching procedure. The binding constants (K_A) were larger than 10^?5 M^?1 and the number of binding sites (n) was approximate to 1 at different temperatures, which indicated that the binding affinity was hige and there was just one main binding site in HSA for PL. According to thermodynamic parameters from Van't Hoff equation, the binding process of PL with HSA was spontaneous and exothermic process due to ΔG < 0, and the electrostatic force played major role in the binding between PL and HSA according to ΔH < 0 and ΔS > 0. The binding distance (r) was calculated to be about 3.35 nm, which implied that the energy transfer from HSA to PL occurred with high possibility according to the theory of Förster's non-radiation energy transfer. The microenvironment and conformation of HSA changed with the addition of PL based on the results of synchronous and three-dimensional fluorescence methods. The molecular docking analysis revealed the binding locus of PL to HSA in subdomain IIIA (Sudlow's site II).     

FWAVina: A novel optimization algorithm for protein-ligand docking based on the fireworks algorithm

Protein-ligand docking is an essential process that has accelerated drug discovery. How to accurately and effectively optimize the predominant position and orientation of ligands in the binding pocket of a target protein is a major challenge. This paper proposed a novel ligand binding pose search method called FWAVina based on the fireworks algorithm, which combined the fireworks algorithm with the efficient Broyden-Fletcher-Goldfarb-Shannon local search method adopted in AutoDock Vina to address the pose search problem in docking. The FWA was used as a global optimizer to rapidly search promising poses, and the Broyden-Fletcher-Goldfarb-Shannon method was incorporated into FWAVina to perform an exact local search. FWAVina was developed and tested on the PDBbind and DUD-E datasets. The docking performance of FWAVina was compared with the original Vina program. The results showed that FWAVina achieves a remarkable execution time reduction of more than 50 % than Vina without compromising the prediction accuracies in the docking and virtual screening experiments. In addition, the increase in the number of ligand rotatable bonds has almost no effect on the efficiency of FWAVina. The higher accuracy, faster convergence and improved stability make the FWAVina method a better choice of docking tool for computer-aided drug design. The source code is available at https://github.com/eddyblue/FWAVina/.     

Combined QSAR/QSPR,Molecular Docking,and Molecular Dynamics Study of Environmentally Friendly PBDEs with Improved Insulating Properties

To improve the insulating properties of polybrominated diphenyl ethers(PBDEs), we studied the molecular structures and energy gap(Eg) values of 209 PBDEs using a three-dimensional quantitative structure-activity relationship(3D-QSAR) model, molecular docking, and molecular dynamics. We also analyzed the interaction mechanisms of PBDEs using a 2D-QSAR model, molecular substitution characteristics, and molecular docking. The 3D-QSAR model showed that the 2-, 4-, 5-, and 6-positions significantly influenced the PBDE insulating properties. Using BDE-34 as a template molecule, we designed six derivatives with 0.47%-28.44% higher insulation tlian BDE-34. Compared with BDE-34, the stability and flame retardancy of the above six derivatives were not adversely affected. These derivatives, except for 2,6-cyanomethyl-BDE, 2-cyanomethyl-BDE, and 2-aniinomethyl-BDE, were more toxic and biodegradable than BDE-34, but showed weaker bioaccumulation and migration abilities than BDE-34. Mechanism analysis showed that the highest occupied orbital energy, the most negative charge, and the dipole moment were the main quantitative parameters that aflected the PBDE insulating properties. PBDE insulation gradually decreased as the number of Br atoms increased. The level of similarity between the substitution patterns on the two benzene rings was significantly correlated with PBDE insulation, with hydrophobic groups having a more significant efiect on PBDE insulation.     

A comparison of heuristic search algorithms for molecular docking

This paper describes the implementation and comparison of four heuristic search algorithms (genetic algorithm, evolutionary programming, simulated annealing and tabu search) and a random search procedure for flexible molecular docking. To our knowledge, this is the first application of the tabu search algorithm in this area. The algorithms are compared using a recently described fast molecular recognition potential function and a diverse set of five protein–ligand systems. Statistical analysis of the results indicates that overall the genetic algorithm performs best in terms of the median energy of the solutions located. However, tabu search shows a better performance in terms of locating solutions close to the crystallographic ligand conformation. These results suggest that a hybrid search algorithm may give superior results to any of the algorithms alone.     

Synthesis,Molecular Docking and Potential Antioxidant Activity of Di/Trisubstituted Pyridazinone Derivatives

A number of pyridazinone derivatives bearing substituted benzylidene and heterocyclic/aromatic rings at 4^th and 6^th positions, respectively were synthesized in good to moderate yields and screened for antioxidant activity. Antioxidant activity of pyridazinone derivatives was evaluated by using several in vitro radical scavenging methods such as 1,1‐diphenylpicrylhydrazyl (DPPH), hydrogen peroxide (H₂O₂), nitric oxide (NO), reducing power, and metal chelating assay etc. Molegro virtual docker software was used to study the binding affinity of the title compounds with the xanthine oxidoreductase enzyme. Amongst the tested compounds, 5a, 5d, 5g & 5j were found to exhibit excellent antioxidant activity at par with the positive control, ascorbic acid. The molecular docking studies of these compounds demonstrated a good selectivity profile with xanthine oxidoreductase receptors. A preliminary study of the structural‐activity relationship showed that the presence of electron withdrawing group and heterocyclic ring on pyridazinone nucleus are associated with the best potency and selectivity profile. It could be proposed that xanthine oxidoreductase receptor may be involved in observed antioxidant activity of pyridazinone derivatives bearing aromatic ring and benzylidene substituents and thus the synthesized compounds are worthy of further exploration.