应用MM/PBSA方法研究CDK2活性口袋内溶剂水分子对CDK2-配体结合自由能的影响

应用MM/PBSA方法研究了CDK2活性口袋内溶剂水分子对CDK2-配体结合自由能的影响. 结果表明, 活性口袋内溶剂水分子对CDK2-配体相互作用自由能有一定的贡献, 其贡献的大小随配体不同而有所差异, 导致这种差异的主要原因是活性位点内溶剂水分子与蛋白残基和配体之间形成了不同的氢键相互作用网络.     

CDK2-抑制剂结合自由能计算

细胞周期蛋白依赖性激酶Ⅱ(cyclin-dependent kinase 2,CDK2)是一种重要的治疗癌症的靶标.本文中采用分子动力学取样,运用MM-PBSA/GBSA两种方法计算了CDK2-NU6102复合物的绝对结合自由能.通过能量分解的方法考察了CDK2大分子主要残基与配体NU6102之间的相互作用和识别.     

血管紧张素转换酶与抑制肽结合模式的分子动力学研究

应用分子模拟方法研究了血管紧张素转换酶(Angiotensin-converting enzyme,ACE)C端结构域(C-domain)与两种抑制肽(RIGLF/AHEPVK)的结合机制,预测了两个体系的结合模式,提出在C-domain-RIGLF中His353,Asp377,Asp453,Phe457,His513,Tyr523和Phe527为RIGLF主要结合残基,而在C-domainAHEPVK中Gln281,His353,Ser355,Glu384,Lys511,His513和Tyr523等残基起关键作用.应用结合自由能计算比较了两个体系的结合能力,结果表明,RIGLF和AHEPVK均与C-domain活性位点残基存在较强作用,且AHEPVK对C-domain的结合能力较强,与实验结果一致.     

VEGFR-2 与抑制剂Sunitinib 的分子对接及分子动力学研究

用分子对接方法研究了VEGFR-2 和抑制剂Sunitinib 的相互作用模式, 并对其复合物进行了10 ns 的分子动力学(Molecular Dynamics, MD)模拟. 结果表明, 抑制剂Sunitinib 能与VEGFR-2 中位于活性空腔的Glu885, Ile888, His1026,Asp1028, Asp1046 五个氨基酸残基形成疏水作用; 另外, VEGFR-2 中His1026, Cys1024, Asp1046 三个氨基酸残基能与Sunitinib 形成三个作用强度不同的氢键. 这些基团之间的相互作用是Sunitinib 抑制VEGFR-2 活性的关键因素. 研究结果可为VEGFR-2 抑制剂的结构改良、分子设计、合成提供理论参考, 并有助于寻找活性更高、效果更好的抗肿瘤药物.     

VEGFR2活性腔性质以及与抑制剂的结合模式研究

VEGFR2介导肿瘤诱导的血管生成作用, 是抑制肿瘤生长和转移的新靶点. 为深入探讨VEGFR2活性腔性质以及与抑制剂的结合模式, 采用多拷贝同时搜寻法(MCSS)研究VEGFR2活性腔的性质, 然后用分子对接方法对5个已上临床的VEGFR抑制剂与VEGFR2活性腔进行对接计算, 讨论它们的结合模式, 确定与配体结合相关的关键残基. 研究发现: 疏水腔I, II是配体结合的关键区域, 残基Glu915, Cys917是关键的氢键作用位点, Lys866, Glu883和Asp1044形成的极性区域对提高配体亲合力很重要, 疏水腔III和极性腔IV是额外增强配体结合力的区域, IV区的Arg1030可提供额外的氢键作用位点. 本研究可为全新VEGFR2抑制剂的合理药物设计提供理论依据, 为寻找新的抗肿瘤药物奠定基础.     

HIV-1整合酶与芳香二酮酸类抑制剂相互作用的分子模拟研究

用分子对接方法研究了一系列芳香二酮酸类抑制剂与HIV-1整合酶的识别及相互作用. 结果表明, 抑制剂结合到整合酶Asp64～Leu68, Thr115～Phe121, Gln148～Lys159和Mg2＋所构成的口袋区, 抑制机理与5CITEP相似. 采用分子动力学模拟和MM/PBSA方法计算了芳香二酮酸类抑制剂与整合酶之间的结合自由能, 计算结果与实验值相吻合, 平均绝对偏差为3.6 kJ/mol, 体系范德华相互作用和溶剂化效应的非极性项是利于形成复合物的主要因素. 相关性分析结果表明, 结合自由能值与疏水相互作用有较强的线性相关(R＝0.61), 基于此, 用多元线性回归方法给出了一个能较强预测芳香二酮酸类抑制剂与HIV-1整合酶的结合自由能预测模型, 为后续基于抑制剂结构的抗HIV-1药物分子设计提供指导.     

含磷嘧啶类CDK9抑制剂的分子对接、3D-QSAR和分子动力学模拟

选取64个具有潜力的含磷嘧啶类细胞周期依赖性蛋白激酶(CDK9)小分子抑制剂,采用分子对接方法研究了该类小分子与CDK9的结合作用,结果表明,分子构象、氢键形成、疏水性和氨基酸残基Cys106在此类抑制剂与CDK9的结合过程中具有重要作用.在配体叠合的基础上,运用比较分子力场分析(Co MFA)、比较分子相似性指数分析(Co MSIA)和Topomer Co MFA(T-COMFA)研究了分子结构与抑制活性的关系,发现由训练集立体场、静电场和疏水场组合的Co MSIA模型为最优模型,其内部交叉验证相关系数(Q2=0.557)、非交叉验证相关系数(R2=0.959)和外部预测相关系数(r2=0.863)具有统计学意义,该模型的三维等值线图直观显示了化合物的活性与其三维结构的关系.根据这些结果设计了10个具有新结构的含磷嘧啶类化合物,分子对接和分子动力学模拟结果表明,新化合物和CDK9的结合模式与原化合物64相同,自由能分析从理论上证明了新化合物64d的CDK9抑制活性优于化合物64,并且显示含磷基团与残基Asp109的静电场能在化合物与CDK9作用过程中有重要作用.     

天然产物法卡林二醇与人类GABAA受体相互作用的机制

通过分子动力学模拟、伞形采样模拟、结合自由能计算和分子对接等方法,研究了法卡林二醇在γ-氨基丁酸A型(GABA_A)受体上结合作用的模式和对GABA_A受体动态属性的影响,确定了3个有效结合位点均对此受体产生拮抗作用,为后续研究聚乙炔醇类化合物对GABA_A受体作用及相应药物开发提供了理论依据.     

分子动力学模拟研究Fedratinib-JAK2/JAK3选择性

通过动力学模拟获得JAK2高选择性抑制剂Fedratinib在JAK2和JAK3激酶中的结合构象,结合自由能的计算结果表明Fedratinib在JAK2中更稳定.将能量分解到结合位点氨基酸,分析发现当分子在JAK2中占据P-loop区的疏水口袋,并与附近Arg980和Asp994等氨基酸形成氢键时,可以增加相对于JAK2的选择性.     

多巴胺第三受体蛋白三维结构及其活性位点氨基酸残基

基于牛视紫红质模板蛋白,同源模建多巴胺第三受体（D₃R）蛋白三维结构,在1-棕榈酰-2-油酰-卵磷脂（POPC）膜-水模型环境,开展300 ns分子动力学模拟提炼优化其结构,取得稳定的D₃R蛋白三维结构（2B08-D₃R）.在该蛋白基础上,采用MP2/6-31G（d,p）方法,计算多巴胺（Dop）与氨基酸残基相互作用的结合能,确定五个残基（Asp117、Ser208、His272、Phe269和Thr276）为活性位点.五个活性位点残基分别位于D₃R蛋白跨膜螺旋区TM3、TM5和TM6,组成活性空腔结构.多巴胺分子以其苯基平面与TM2-TM7包围的圆柱体空腔平行和非共价键结合方式保留在D₃R蛋白中,与D₃R蛋白结合能E_b为-97.8 kJ·mol^-1基于3PBL D₃R突变体晶体结构,构建了另外一个含有多巴胺分子的D₃R蛋白结构（Dop-3PBL-D₃R）,确定在该蛋白结构中,多巴胺的活性位点氨基酸是Asp83、His272、Phe269、Phe268和Trp265.在该蛋白结构中,多巴胺分子同样以其苯基平面与TM2-TM7包围的圆柱体空腔平行和非共价键方式结合,与该蛋白相互作用的结合能是-80.5 kJ·mol^-1.     

An analysis of ATP binding with kinase catalytic subunit by molecular dynamics simulation of the CDK2 active kinase crystal lattice

Molecular dynamics simulation was used to analyze changes in the functionally significant structural elements of the crystal lattices of pT160-CDK2/cyclin and A/ATP-Mg²⁺/substrate complexes of the native (CDK2-G16) and mutant (CDK2-S16) active kinases at physiological temperatures (300 K). The structural rearrangement of ATP caused by changes in the kinase catalytic domain was studied. ATP was fixed by the ionic and H-bond interactions of several residues, including Lys33, Asp145, and side-chain amides of the G loop between β1 and β2. The binding of the kinases to complexes with cyclin and the phosphorylation of T160 in the active complex of the CDK2 kinase result in the ATP orientation more convenient for the transfer of the phosphate group to the substrate. An analysis of interatomic distances in the ATP active site region and Asp145, Asn132, Lys33 catalytic sites participating in the orientation of ATP phosphates revealed that the Asp 145 amino acid residue was situated noticeably closer to the ATP molecule in the native complex than in its mutant counterpart. The same is true of the arrangement of the Lys33 residue with respect to ATP.     

In silico study of porphyrin-anthraquinone hybrids as CDK2 inhibitor

Cyclin-Dependent Kinases (CDKs) are known to play crucial roles in controlling cell cycle progression of eukaryotic cell and inhibition of their activity has long been considered as potential strategy in anti-cancer drug research. In the present work, a series of porphyrin-anthraquinone hybrids bearing meso-substituents, i.e. either pyridine or pyrazole rings were designed and computationally evaluated for their Cyclin Dependent Kinase-2 (CDK2) inhibitory activity using molecular docking, molecular dynamics simulation, and binding free energy calculation. The molecular docking simulation revealed that all six porphyrin hybrids were able to bind to ATP-binding site of CDK2 and interacted with key residues constituted the active cavity of CDK2, while molecular dynamics simulation indicated that all porphyrins bound to CDK2 were stable for 6 ns. The binding free energies predicted by MM-PBSA method showed that most compounds exhibited higher affinity than that of native ligand (4-anilinoquinazoline, DTQ) and the affinity of mono-H₂PyP-AQ was about three times better than that of DTQ, indicating its potential to be advanced as a new CDK2 inhibitor.     

Virtual screening studies to design potent CDK2-cyclin A inhibitors

The cell division cycle is controlled by cyclin-dependent kinases (CDK), which consist of a catalytic subunit (CDK1-CDK8) and a regulatory subunit (cyclin A-H). Pharmacophore analysis indicates that the best inhibitor model consists of (1) two hydrogen bond acceptors, (2) one hydrogen bond donor, and (3) one hydrophobic feature. The HypoRefine pharmacophore model gave an enrichment factor of 1.31 and goodness of fit score of 0.76. Docking studies were carried out to explore the structural requirements for the CDK2-cyclin A inhibitors and to construct highly predictive models for the design of new inhibitors. Docking studies demonstrate the important role of hydrogen bond and hydrophobic interactions in determining the inhibitor-receptor binding affinity. The validated pharmacophore model is further used for retrieving the most active hits/lead from a virtual library of molecules. Subsequently, docking studies were performed on the hits, and novel series of potent leads were suggested based on the interaction energy between CDK2-cyclin A and the putative inhibitors.     

Computational Analysis for Residue-Speci c CDK2-Inhibitor Bindings

Cyclin-dependent kinase 2 (CDK2) is a key macromolecule in cell cycle regulation. In cancer cells, CDK2 is often overexpressed and its inhibition is an effective therapy of many cancers including breast carcinomas, leukemia, and lymphomas. Quantitative characterization of the interactions between CDK2 and its inhibitors at atomic level may provide a deep understanding of protein-inhibitor interactions and clues for more effective drug discovery. In this study, we have used the computational alanine scanning approach in combination with an efficient interaction entropy method to study the microscopic mechanism of binding between CDK2 and its 13 inhibitors. The total binding free energy from the method shows a correlation of 0.76?0.83 with the experimental values. The free energy component reveals two binding mode in the 13 complexes, namely van der Waals dominant, and electrostatic dominant. Decomposition of the total energy to per-residue contribution allows us to identify five hydrophobic residues as hot spots during the binding. Residues that are responsible for determining the strength of the binding were also analyzed.