神经介素B受体与拮抗剂、激动剂的分子模拟研究

大鼠神经介素B受体(rat neuromedin B receptor, rNMBR)属于G蛋白偶联受体(G-protein coupled receptor, GPCR) A家族的成员. GPCR的结构特征和在信号传导中的重要作用决定了其可以作为很好的药物靶标. 关于rNMBR与内源性激动剂神经介素B (neuromedin B, NMB)以及与非肽类拮抗剂pd168368作用机制的研究对于合理设计受体药物分子有重要的指导意义. 在这一研究中, 我们使用同源模建, 构建受体的三维结构, 进行分子对接和分子动力学的计算. 基于受体三维结构, 通过10 ns的空载受体、激动剂-受体、拮抗剂-受体的分子动力学模拟, 探讨受体与激动剂与拮抗剂的作用机制. 研究表明rNMB-R中跨膜(transmembrane, TM)螺旋3, 5, 6, 7参与配体的结合. NMB与受体的结合, 使受体转变为活性构象, 而受体同拮抗剂pd168368恰好相反.     

毒蕈碱型M1受体的同源模建和分子对接

采用同源模建方法对M1受体的三维结构进行了模拟,将得到的模型分别与M受体完全激动剂乙酰胆碱和M1受体选择性激动剂占诺美林进行分子对接,形成非特异性激动和特异性激动的受体-配体复合物.用分子动力学模拟方法分别将未与小分子对接的M1受体、M1受体-乙酰且H碱复合物、M1受体-占诺美林复合物置于磷脂双膜中模拟10 ns.将模拟后的蛋白质结构与包含活性分子的测试库对接并将结果打分,以top5%富集因子(EF)作为评价依据,用占诺美林优化后的M1受体模型的EF为8.0,用乙酰胆碱优化后M1受体模型的EF为6.5,非复合物的EF为1.5.说明M1受体选择性激动剂复合物进行分子动力学模拟后得到的三维结构模型比较合理,可以作为化合物虚拟筛选的模型对新化合物进行虚拟筛选,为找到新的选择性M1受体激动剂奠定了基础.     

毒蕈碱型M1受体的同源模建和分子对接

采用同源模建方法对M1受体的三维结构进行了模拟, 将得到的模型分别与M受体完全激动剂乙酰胆碱和M1受体选择性激动剂占诺美林进行分子对接, 形成非特异性激动和特异性激动的受体-配体复合物. 用分子动力学模拟方法分别将未与小分子对接的M1受体、M1受体-乙酰胆碱复合物、M1受体-占诺美林复合物置于磷脂双膜中模拟10 ns. 将模拟后的蛋白质结构与包含活性分子的测试库对接并将结果打分, 以top5%富集因子(EF)作为评价依据, 用占诺美林优化后的M1受体模型的EF为8.0, 用乙酰胆碱优化后M1受体模型的EF为6.5, 非复合物的EF为1.5. 说明M1受体选择性激动剂复合物进行分子动力学模拟后得到的三维结构模型比较合理, 可以作为化合物虚拟筛选的模型对新化合物进行虚拟筛选, 为找到新的选择性M1受体激动剂奠定了基础.     

三磷酸腺苷结合盒式转运体的分子模拟

三磷酸腺苷结合盒式(ATP binding cassete,ABC)转运体是一类与生物体的生理过程和疾病密切相关的膜蛋白,它们利用水解ATP释放的能量实现底物的跨膜转运。虽然研究者已对ABC转运体进行了广泛而深入的研究,但是由于ABC转运体结构的复杂性,导致许多实验技术在ABC转运体的研究中受到较大限制。分子模拟技术可以弥补实验技术的不足,已成为不可或缺的研究工具。本文综述了近年来分子模拟技术同源模建、分子对接和分子动力学模拟在ABC转运体的三维结构模建,底物结合位点确定以及构象转换分子机理的解析等研究领域的最新进展。最后总结了分子模拟技术在ABC转运体研究中所遇到的挑战和拟解决办法。     

芬太尼类化合物与阿片μ受体相互作用的分子对接与分子动力学模拟

采用分子对接和分子动力学(MD)模拟方法研究了芬太尼类化合物与阿片μ受体的相互作用机制.先用AutoDock4.0程序将芬太尼类化合物对接到同源模建的阿片μ受体结构中,再用GROMACS程序包在水溶液体系中分别对12个芬太尼激动剂和阿片μ受体蛋白复合物进行了MD模拟研究,优化对接复合物的结构,最后利用MM-PBSA方法,在APBS程序中计算芬太尼类衍生物与阿片μ受体的结合自由能,计算出的受体配合物结合常数(Ki)与其实验值吻合较好,并预测了化合物的活性排序.结果表明,复合物蛋白结构与空载受体蛋白结构有较大差异,特别是胞内区IL2、IL3和跨膜区段TM4骨架构象变化较大,不同的化合物对受体结构影响也有差异,活性较好的化合物会增加蛋白特定区域结构的柔性.芬太尼类化合物可能是通过和受体结合后诱导阿片μ受体构象转变为活性构象,引起一系列的信号传导激活G蛋白,从而引发生理效应.     

AT1受体的中药活性成分筛选模型及其作用机理研究

传统中药对治疗心血管类疾病疗效显著,例如钩藤、黄芪、益母草等在临床应用广泛.现代药理研究表明钩藤碱可以降压;黄芪中毛蕊异黄酮能舒张血管平滑肌、保护心脑血管;益母草碱可扩张微血管,改善血液流变异常,但它们分子层面作用机制尚不明确.首先以牛视紫红质蛋白为模板,模建出心血管疾病主要靶点AT1受体的三维结构.然后将AT1受体拮抗剂和中药活性成分与受体模型结合的作用方式进行了对比研究,据此提出了中药活性成分治疗心血管疾病的作用机理,并建立了AT1受体的中药活性成分筛选模型.结果表明:黄芪毛蕊异黄酮等中药活性成分能与AT1受体活性口袋中的残基发生氢键作用,结合方式与AT1受体拮抗剂相似.每一种AT1受体拮抗剂均与His183,Lys199,His256,Gln257,Ser105,Ser109,Tyr113,Asn200中多个发生氢键作用;黄芪毛蕊异黄酮与Try113,Lys199,Gln257,Ser105发生氢键作用.本研究从分子层面上阐释了一些中药活性小分子的治病机理,为进一步挖掘中药资源,研究AT1受体相关的心脑血管类药物,合理设计和筛选AT1受体的拮抗剂提供重要依据.     

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

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

鲍曼不动杆菌LpxC的同源模建及分子对接

鲍曼不动杆菌已成为最普遍的医院致病菌,且耐药情况严峻.LpxC作为新抗菌药物靶点被大量研究,但鲍曼不动杆菌LpxC晶体尚未解析得到,基于其结构的药物设计等工作无法开展.以铜绿假单胞菌LpxC晶体结构为模板,通过同源模建方法获得鲍曼不动杆菌LpxC结构模型.较好的Ramachandran plot分布和Profile-3D结果验证了模型的合理性.用分子动力学模拟优化鲍曼不动杆菌LpxC模型,修补部分不合理构象.后续分子对接结果显示S构型的苄氧乙酰基羟肟酸类抑制剂比R构型分子能更有效地结合在F191,H237和K238组成的较浅口袋中,这可能是S构型抑制剂活性更高的主要因素,模拟结果与实验数据吻合较好.     

拮抗状态下α1A1，α1B和α1D-肾上腺素能受体的分子模拟研究

采用同源建模法对α1A1，α1B和α1D-的三维结构进行了模拟，并采用分子力学、分子动力学方法对所得同源模型进行优化，然后分别采用训练集拮抗剂对接的方法得到拮抗状态下的α1A1，α1B和α1D-三维结构模型．得到的模型再采用FRED对接软件对测试集中的18个化合物进行对接并打分，再将所得打分结果与其活性进行线性回归，其回归结果具有良好的拟合效果，由此回归方程预测的活性与化合物实验值较吻合，说明我们建立的拮抗状态下的α1A1，α1B和α1D-的三维同源模型具有一定的合理性，可作为化合物虚拟筛选模型，对新化合物进行对接虚拟筛选．     

CCK1受体的同源模拟和分子对接研究

采用同源建模法对CCK1受体的三维结构进行了模拟,并采用分子动力学方法对模型进行修正和优化,再采用与训练集激动剂和拮抗剂分子对接的方法分别得到激动状态和拮抗状态CCK1受体的三维结构模型。得到的模型使用DOCK对接软件对训练集中的分子进行对接,所得结果与其实际活性拟合度较好,说明我们建立的激动和拮抗状态下的CCK1受体的三维结构模型比较合理,可以作为化合物虚拟筛选的模型对新化合物进行虚拟筛选。     

Apicidin选择性抑制ClassⅠ HDACs分子动力学研究

采用模拟方法研究组蛋白脱乙酰酶抑制剂(Apicidin)选择性抑制组蛋白去乙酰化酶(Histone deacetylases, HDACs)中的HDAC1和HDAC8. 通过HDAC8晶体结构同源模建HDAC1三维结构模型, 将Apicidin分别与HDAC1和HDAC8对接并进行分子动力学模拟, 结果表明, HDAC1活性口袋入口处的Arg270是Apicidin-HDAC1形成稳定结构的重要因素; HDAC1中Tyr303及His178与Apicidin形成2个持续存在的氢键, 而在HDAC8中未发现, 这是Apicidin选择性抑制HDAC1高于HDAC8的另一重要原因.     

Modification of PAE-degrading Esterase(CarEW) for Higher Degradation Efficiency Through Integrated Homology Modeling,Molecular Docking,and Molecular Dynamics Simulation

Six phthalate acid esters(PAEs) priority pollutants[dimethyl phthalate(DMP), diethyl phthalate(DEP), dibutyl phthalate (DBP or DNBP), di-n-octyl phthalate(DNOP), di 2-ethyl hexyl phthalate(DEHP), and butyl benzyl phthalate(BBP)] were opted as the research object. PAE-degrading esterase CarEW(PDB ID:1C7I) isolated from Bacillus subtilis acting as a template and an iterative saturation mutation strategy was adopted to modify key amino acids to attain efficient PAE-degrading esterase substitutes with a reasonable structure constructed by homology modeling method. Present study designed a total of 285 unit-site and multi-site substitutions of PAE-degrading esterase using the homology modeling method. Among them, 207 PAE-degrading esterase substitutions, which contained the 6-site PAE-degrading esterase substitute 1C7I-6-9 with 84.21% enhancement intensity of degradation ability revealed better degradability to all the 6 PAEs after modification. Moreover, molecular dynamics simulation based on the Taguchi method reported the optimal external application environment for PAE-degrading esterase substitutes as follows:pH=6, T=35℃, the rhamnolipid concentration was 50 mg/L, the molar ratio of nitrogen to phosphorus(N:P) was 10:1, the concentration of H₂O₂ was 50 mg/L, and the voltage gradient was 1.5 V/cm. The degradation ability of PAE-degrading esterase substitutes was found to be elevated by 13.04% as compared to that of the blank control under the optimal condition. Moreover, 11 highly efficient PAE-degrading esterase substitutes with thermal stability were designed.     

Homology modeling and 3D–QSAR study of benzhydrylpiperazine δ opioid receptor agonists

The binding affinity of a series of benzhydrylpiperazine δ opioid receptor agonists were pooled and evaluated by using 3D-QSAR and homology modeling/molecular docking methods. Ligand-based CoMFA and CoMSIA 3D-QSAR analyses with 46 compounds were performed on benzhydrylpiperazine analogues by taking the most active compound BW373U86 as the template. The models were generated successfully with q² value of 0.508 and r² value of 0.964 for CoMFA, and q² value of 0.530 and r² value of 0.927 for CoMSIA. The predictive capabilities of the two models were validated on the test set with R²_pred value of 0.720 and 0.814, respectively. The CoMSIA model appeared to work better in this case. A homology model of active form of δ opioid receptor was established by Swiss-Model using a reported crystal structure of active μ opioid receptor as a template, and was further optimized using nanosecond scale molecular dynamics simulation. The most active compound BW373U86 was docked to the active site of δ opioid receptor and the lowest energy binding pose was then used to identify binding residues such as s Gln105, Lys108, Leu125, Asp128, Tyr129, Leu200, Met132, Met199, Lys214, Trp274, Ile277, Ile304 and Tyr308. The docking and 3D-QSAR results showed that hydrogen bond and hydrophobic interactions played major roles in ligand-receptor interactions. Our results highlight that an approach combining structure-based homology modeling/molecular docking and ligand-based 3D-QSAR methods could be useful in designing of new opioid receptor agonists.     

Identification of novel Plasmodium falciparum PI4KB inhibitors as potential anti-malarial drugs: Homology modeling,molecular docking and molecular dynamics simulations

The current study was set to discover selective Plasmodium falciparum phosphatidylinositol-4-OH kinase type III beta (pfPI4KB) inhibitors as potential antimalarial agents using combined structure-based and ligand-based drug discovery approach. A comparative model of pfPI4KB was first constructed and validated using molecular docking techniques. Performance of Autodock4.2 and Vina4 software in predicting the inhibitor-PI4KB binding mode and energy was assessed based on two Test Sets: Test Set I contained five ligands with resolved crystal structures with PI4KB, while Test Set II considered eleven compounds with known IC50 value towards PI4KB. The outperformance of Autodock as compared to Vina was reported, giving a correlation coefficient (R²) value of 0.87 and 0.90 for Test Set I and Test Set II, respectively. Pharmacophore-based screening was then conducted to identify drug-like molecules from ZINC database with physicochemical similarity to two potent pfPI4KB inhibitors –namely cpa and cpb. For each query inhibitor, the best 1000 hits in terms of TanimotoCombo scores were selected and subjected to molecular docking and molecular dynamics (MD) calculations. Binding energy was then estimated using molecular mechanics–generalized Born surface area (MM-GBSA) approach over 50 ns MD simulations of the inhibitor-pfPI4KB complexes. According to the calculated MM-GBSA binding energies, ZINC78988474 and ZINC20564116 were identified as potent pfPI4KB inhibitors with binding energies better than those of cpa and cpb, with ΔG_binding ≥ −34.56 kcal/mol. The inhibitor-pfPI4KB interaction and stability were examined over 50 ns MD simulation; as well the selectivity of the identified inhibitors towards pfPI4KB over PI4KB was reported.     

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

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

人类2-氨基3-羧基粘康酸6-半醛脱羧酶(ACMSD)与底物及抑制剂作用模型的理论研究

利用同源模建和分子动力学模拟方法构建了人类2-氨基3-羧基粘康酸6-半醛脱羧酶(hACMSD)的三维结构, 并利用Profile-3D和Procheck等方法评估了模型的可靠性. 在此基础上, 用分子对接程序(Affinity), 将其底物2-氨基3-羧基粘康酸6-半醛(ACMS)和抑制剂喹啉酸(QA)分别与hACMSD进行对接, 获得了复合物结构的理论模型. 通过配体与受体之间相互作用能和结构分析给出了底物和抑制剂的具体结合方式, 明确了hACMSD与底物和抑制剂结合时起重要作用的氨基酸残基.     

Comprehensive structural analysis of the open and closed conformations of Theileria annulata enolase by molecular modelling and docking

Theileria annulata is an apicomplexan parasite which is responsible for tropical theileriosis in cattle. Due to resistance of T. annulata against commonly used antitheilerial drug, new drug candidates should be identified urgently. Enolase might be a druggable protein candidate which has an important role in glycolysis, and could also be related to several cellular functions as a moonlight protein. In this study; we have described three-dimensional models of open and closed conformations of T. annulata enolase by homology modeling method for the first time with the comprehensive domain, active site and docking analyses. Our results show that the enolase has similar folding patterns within enolase superfamily with conserved catalytic loops and active site residues. We have described specific insertions, possible plasminogen binding sites, electrostatic potential surfaces and positively charged pockets as druggable regions in T. annulata enolase.     

Molecular characterization,modeling and docking of CYP107CB2 from Bacillus lehensis G1, an alkaliphile

Cytochrome P450s are a superfamily of heme monooxygenases which catalyze a wide range of biochemical reactions. The reactions involve the introduction of an oxygen atom into an inactivated carbon of a compound which is essential to produce an intermediate of a hydroxylated product. The diversity of chemical reactions catalyzed by cytochrome P450s has led to their increased demand in numerous industrial and biotechnology applications. A recent study showed that a gene sequence encoding a CYP was found in the genome of Bacillus lehensis G1, and this gene shared structural similarity with the bacterial vitamin D hydroxylase (Vdh) from Pseudonocardia autotrophica. The objectives of present study was to mine, for a novel CYP from a new isolate B. lehensis G1 alkaliphile and determine the biological properties and functionalities of CYP in this bacterium. Our study employed the usage of computational methods to search for the novel CYP from CYP structural databases to identify the conserved pattern, functional domain and sequence properties of the uncharacterized CYP from B. lehensis G1. A computational homology model of the protein’s structure was generated and a docking analysis was performed to provide useful structural knowledge on the enzyme’s possible substrate and their interaction. Sequence analysis indicated that the newly identified CYP, termed CYP107CB2, contained the fingerprint heme binding sequence motif FxxGxxxCxG at position 336-345 as well as other highly conserved motifs characteristic of cytochrome P450 proteins. Using docking studies, we identified Ser-79, Leu-81, Val-231, Val-279, Val-383, Ala-232, Thr-236 and Thr-283 as important active site residues capable of stabilizing interactions with several potential substrates, including vitamin D₃, 25-hydroxyvitamin D₃ and 1α-hydroxyvitamin D₃, in which all substrates docked proximally to the enzyme’s heme center. Biochemical analysis indicated that CYP107CB2 is a biologically active protein to produce 1α,25-dihydroxyvitamin D₃ from 1α-hydroxyvitamin D₃. Based on these results, we conclude that the novel CYP107CB2 identified from B. lehensis G1 is a putative vitamin D hydroxylase which is possibly capable of catalyzing the bioconversion of parental vitamin D₃ to calcitriol, or related metabolic products.