低聚糖是各种有机体的主要能源物质,它可以由地球上最丰富的聚合物—维素的分解转化而得到,也可以由单糖的聚合而得到。这两种转化都要靠糖苷酶的催化作用来实现。在纤维素水解过程中起催化作用的酶称为糖苷水解酶,而在单糖聚合过程中起催化作用的酶称为糖基合成酶。目前实验上已经得到了许多糖苷水解酶和糖基合成酶的晶体结构,并且通过突变实验得到了活性口袋处起催化作用的残基,另外,人们从理论上也对一些糖苷酶的催化机理进行了研究。但是总的来说,仍然有很多糖苷酶的反应机理尚不十分明确。本论文利用量子力学和分子力学相结合的方法(QM/MM)对几种糖苷酶的催化机理进行了研究。关于药物作用机制及生物大分子的分子模拟也是化学研究的一个热点领域,深入探索药物分子和受体蛋白之间的作用机制对于药物的筛选和设计具有重要的指导意义。维生素D受体、孕酮受体、人的EP3受体等均是重要的靶点,可用于治疗骨质疏松、乳腺痛、动脉硬化等一系列疾病,本论文用多种,理论化学方法(分子模拟、QM/MM和3D-QSAR)深入探索了几类新型药物分子与这些药物靶点的作用模式、生物大分子在药物作用下的构象变化以及在此研究基础上的药物设计。本论文主要工作描述如下：(1)水稻p-型葡萄糖苷酶催化机理研究。◆糖基化和去糖基化机理研究。从微生物酶进化而来的水稻p-型葡萄糖苷酶(BGlul,即3Os3bglu7)比其它植物酶的同源性更好,在水解低聚糖方面更有效。水稻BGlul葡萄糖苷酶的催化水解过程遵循的是双取代机理,主要包括糖基化和去糖基化两个步骤。为了探索这种酶的催化机理,我们用联合的量子力学／分子力学(QM/MM)方法对水稻BGlul葡萄糖苷酶催化昆布二糖水解的糖基化机理进行了研究。计算结果表明,此反应的糖基化是一个协同反应,即作为酸碱催化剂(acid/base)的E176进攻糖基氧和作为亲核试剂的E386进攻异头碳是同步的；反应形成的过渡态具有阳离子特性,反应能垒为15.7kcal/mol;在过渡中端基糖环向椅式构象转化,平面化趋势加强。本文对去糖基化过程也进行了研究。研究表明,此反应也经历一个协同过程到达最终的产物,反应能垒为21.4kcal/mol；反应形成的过渡态同样具有碳正离子特性;在过渡态中,E386与底物的C2-OH之间仍有氢桥相互作用。这种氢键相互作用有利于糖环由椅式构象变成半椅式构象,并且可以稳定过渡态的碳正离子,与实验结果相符。至此,我们已经系统地研究了水稻葡萄糖苷酶催化昆布二糖水解的糖基化和去糖基化两个过程,反应物构型与传统的量子化学方法有相似之处,但是反应的能垒更低。希望我们的研究结果可以为今后人们更好地研究糖苷酶提供理论性的指导作用。◆E386突变体催化寡糖合成机理研究。作为高效率的糖苷合成酶,β-型葡萄糖苷酶的亲核残基突变体可以用来合成低聚糖。我们采用QM/MM组合方法进一步研究了水稻BGlul葡萄糖苷酶亲核残基E386突变体(E386G、E386S和E386A)合成低聚糖的催化机理。研究结果表明,此合成机理是单步反应,即E176夺取氢质子和糖苷键的形成是一个协同过程；反应能垒受突变残基侧链空间位阻的影响比较大,突变体E386G、E386S和E386A的催化反应能垒分别为22.4、25.7和28.2kcal/mol,很好的解释了相关的突变实验活性的递变结果。(2)多形拟杆菌α-型葡萄糖苷酶水解机理研究。作为重要的革兰氏阳性菌种,多形拟杆菌(B. thetaiotaomicron)是一种主要寄生在人的内脏器官的寄生菌。作为糖苷酶GH97家族中的一员,BtGH97a是一种α-型葡萄糖苷酶。这种酶可以把碳水化合物分解成可供寄主吸收的物质。因此,我们采用QM/MM组合方法对BtGH97a催化底物pNP-Glc水解的机理进行了研究,主要分析了两种可能的分步反应路径：亲核进攻优先路径和质子进攻优先路径。在第一种路径中,首先一个水分子被亲核碱残基(E439或E508)去质子化再进攻底物pNP-Glc的异头碳完成亲核进攻过程,然后酸性氨基酸残基E532进攻糖苷氧完成质子进攻过程。在第二条反应路径中,首先是E532进攻糖苷氧完成质子进攻过程,然后去质子化的水(亲核碱为E439或E508)进攻底物的异头碳完成亲核进攻过程。计算结果表明,亲核进攻优先路径是有利的反应路径,而其亲核进攻过程是速率决定步骤(能垒为15.4kcal/mol,E508作为亲核碱)。在这个决速步中,水分子去质子化和糖苷氧的质子化过程是协同的；在质子进攻优先路径中,质子进攻糖苷氧是速率决定步骤,能垒为24.1kcal/mol。所以,通过比较我们发现此水解反应遵循的是亲核进攻优先路径。本次研究工作主要回答了哪一个氨基酸残基是最可能的辅助亲核碱以及最可能的反应路径。(3)维生素D受体构象变化和相关抑制剂的3D-QSAR研究。◆维生素D受体配体结合区的构象变化研究。维生素D受体(VDR)是一种核激素受体。VDR需要结合一定的配体(天然激素或类似物)才能实现其功能表达。我们用分子动力学(MD)模拟方法研究了人类维生素D受体(hVDR)在与1α,25(OH)2D3(天然激素配体)及其两种类似物(在文章中用b和c表示)结合时配体结合区(LBD)的构象变化情况。MD模拟结果显示,这些配体可以引起LBD的不同构象变化：1α,25(OH)2D3只引起很小的的变化,这表明VDR结合天然配体时采取典型的活性构象；而b和c则会导致LBD结构的明显重排,这种结构重排现象明显地扩大了配体结合口袋,与实验测定结果一致。本文的计算结果可以帮助人们更好地认识VDR对于不同抑制剂分子的识别,并可以对人们设计低钙活性的VDR抑制剂提供参考信息。◆维生素D326,23-内酯类抑制剂的3D-QSAR研究。三维定量构效关系(3D-QSAR)方法是药物设计领域一种非常有效的药物设计方法,基于它我们可以设计出具有更高预测活性的抑制剂分子。利用这种方法本文对82个25-脱氢-1α-羟基维生素D3-26,23-内酯类抑制剂进行了研究,并由此构建了比较分子场分析(CoMFA)(?)(?)比较分子相似性指数分析(CoMSIA)模型。CoMFA模型给出的交叉验证(q2)和非交叉验证(r2ncv)的参数值分别为0.516和0.667,而CoMSIA模型给出的q2和r2ncv值分别为0.517和0.632。这两个模型给出的预测性参数rpred2值分别为0.639和0.619。总体上来说,CoMFA模型的预测能力在这里要优于CoMSIA模型。所以,基于CoMFA模型,我们设计了一些高预测活性的维生素D3类抑制剂,并通过自由能围绕(FEP)模拟揭示了这些所设计分子的合理性。(4)孕酮受体抑制剂的3D-QSAR研究。对于54种孕酮受体(PR)抑制剂,我们构建了基于受体的CoMFA和CoMSIA模型。CoMFA模型统计出的结果为：交叉验证系数q2=0.534,非交叉验证系数r2ncv=0.947。CoMSIA模型统计结果为：q2=0.615,r2ncv=0.954。这两个模型给出的预测相关系数r2pred分别为0681(CoMFA)和0.677(CoMSIA)。基于预测性较好的CoMFA模型,我们设计了两类PR抑制剂,均表现出较高的预测活性。对于两种典型的PR抑制剂,我们用QM/MM组合方法给出了其在活性位点处可能的结合模式。这将为PR类抑制剂分子的药物设计提供有用的理论信息。(5)EP3受体抑制剂的3D-QSAR研究。本文对57个人的EP3受体抑制剂进行了3D-QSAR研究,并构建了CoMFA和CoMSIA模型。这两个模型统计出的交叉验证系数q2分别为0.505和0.492,非交叉验证系数r2ncv分别为0.937和0.921,这两个模型得到的预测相关系数分别为0.549和0.540。通过比较我们发现CoMFA模型的预测能力较好。所以,基于CoMFA模型,对相关EP3抑制剂的主要结构参数进行了确认,并设计了一些具有较高预测活性的抑制剂。这将为EP3类抑制剂的设计和开发提供有用的理论信息。本研究论文的主要创新点如下：◆用QM/MM组合方法系统地研究了水稻BGlul葡萄糖苷酶的保持型催化机理及其突变体催化低聚糖合成的机理。计算中充分考虑了蛋白质环境和活性中心残基对反应的影响,可以提供更加接近实验事实的结构和能垒数据,能够更好地揭示水稻BGlul葡萄糖苷酶催化作用的实质,并且可以得到一些单从晶体结构无法得到的信息(如对E386G、E386S (?)口E386A突变体催化低聚糖合成机理的研究,我们成功地解释三种突变体活性顺次减弱的原因),在分子水平上阐明了水稻葡萄糖苷酶的催化机制。◆用QM/MM方法系统研究了多形拟杆菌α-型葡萄糖苷酶翻转型催化机理,在研究中我们设计了两种可能的反应路径,通过对对比四种可能进行的反应情况研究,得到了最可能的反应路径和对催化反应起关键作用的残基。◆系统研究了几种重要靶蛋白如维生素D、孕酮受体和前列腺素受体与相应药物分子作用机制,应用3D-QSAR模型得到了这些药物分子的结构一活性关系并构建了具有高度预测性的构效关系模型,设计了药物活性更好的药物分子,引入QM/MM组合方法精确地给出了新设计抑制剂的结合模式,同时结合自由能围绕方法(FEP)来评价新设计抑制剂的合理性。这为今后相关的药物筛选和设计提供了有价值的理论指导信息。