论文部分内容阅读
信息技术与生物技术的结合日益紧密。随着计算机技术的发展和成熟,计算机模拟(又称分子模拟)在蛋白质结构与功能关系的研究中扮演越来越重要的角色。本论文以实验室自主筛选的一株微生物来源的脂肪酶——LipK107为主要研究对象,采用多种不同的分子模拟方法(包括:同源建模、分子对接、分子动力学、MM/PBSA结合自由能计算)详尽分析了该脂肪酶结构与功能之间的关系。并基于上述构效关系的研究,对LipK107进行了理性设计,获得了2株催化特性得到改良的脂肪酶突变体。同源建模是目前较为准确的蛋白质三维结构预测方法。作为LipK107结构与功能研究的基础,本文采用同源建模法同时构建了LipK107活性和非活性构象的三维结构。通过比较两个模型,本文成功预测出LipK107的关键结构特性,包括:酶的“催化三联体”(S79-D232-H254)、“氧负离子洞”(L13,Q80)、和“盖子”结构(1118-D152)。通过观察LipK107在催化甘油三酯水解时表现出“界面活化”效应,我们更确证了LipK107盖子结构的存在性。将LipK107的两个同源模型叠合之后,我们发现了该脂肪酶从活性构象转变为非活性构象的过程。此外,采用分子动力学模拟技术,我们也模拟了脂肪酶在水溶液中的构象变化过程,直观并且实时地观察到LipK107通过关闭盖子从活性构象变为非活性构象的过程。分子对接被广泛用于探索配体与受体之间的结合模式及相互作用机理。基于南极假丝酵母脂肪酶(Candida antartica lipase B,简称CALB)的空间结构与催化机理,本文采用全柔性分子对接程序Affinity开发了一套适用于脂肪酶的底物虚拟筛选系统(Computer-aided substrate screening,简称CASS)。使用CASS对CALB进行底物筛选,我们发现该筛选系统的准确率高达95.7%。此外,我们也使用CASS对LipK107的底物进行了筛选。预测结果表明,10个待筛选化合物中有4个是LipK107的底物,其余6个为非底物。后续的生物催化实验结果证明上述预测结果中只有1个化合物的预测结果是错误的,这表明CASS的准确性高达90%。从作者对文献的调研情况来看,CASS是目前唯一一个酶的底物虚拟筛选系统。酶催化反应过程涉及化学键的断裂和形成,而分子对接却无法实现这个过程。所以本文采用分子动力学模拟反应过渡态的方法探索了CALB的对映体选择性,发现了一种新的慢反应对映体(slow-reacting enantiomer)与酶的有效结合模式——“M/H交换”(M/Hpermutation)。使用相同的模拟思路和方法,我们预测了LipK107转酯拆分1-苯基乙醇时的对映体选择性。预测结果表明R-型底物为快反应对映体。该预测结果的正确性随后被生物催化实验结果证实。基于LipK107盖子结构的疏水和电荷分布特性,本文使用一组分子模拟方法(按照使用顺序的先后,该方法依次包括虚拟突变、分子对接、分子动力学和MM/PBSA)研究了它们对酶活的影响,并虚拟构建了相关的突变蛋白(包括I128E/V129D, E130L/K131I, Y138V, R120P/K121P, H146S/R147T)。通过分析突变体与底物结合自由能的高低,我们预测出只有提高盖子结构疏水性的两个突变体(E130L/K131I, Y138V)才能提高酶活,其余突变体的活性都会下降。随后的定点突变实验结果表明:较之野生LipK107,所有的突变体中只有E130L/K131I、Y138V催化1-苯基乙醇转酯拆分反应的转化率和产物ee值有所提高。这样的实验结果与模拟预测的结果完全一致。这不仅证明了之前预测结果的准确性,也暗示着“虚拟突变+分子对接+分子动力学+MM/PBSA"这一分子模拟方法的组合可以作为酶理性设计的一种新方法。除上述基于盖子结构疏水和电荷分布的理性设计外,本文还采用虚拟突变技术设计了具有新催化三联体的突变蛋白(D232A/N103D)。然而,与模拟结果不一致的是,所设计的突变株在后续的真实实验中失活。使用分子动力学方法,本文解析了该突变株失活的原因:新催化三联体(S79-D103-H254)中D103与H254之间的氢键无法稳定存在。