《生物分子的结构及其相互作用的计算化学研究》由会员分享，可在线阅读，更多相关《生物分子的结构及其相互作用的计算化学研究（34页珍藏版）》请在金锄头文库上搜索。
1、有机化学专业优秀论文有机化学专业优秀论文 生物分子的结构及其相互作用的计算化学生物分子的结构及其相互作用的计算化学研究研究关键词：可卡因关键词：可卡因 尼古丁尼古丁 磷酸二酯酶磷酸二酯酶 分子力学方法分子力学方法 分子动力学模拟分子动力学模拟 计算化学计算化学 生物分子生物分子摘要：本文运用一系列的理论化学方法研究了一些生物小分子(尼古丁、可卡因、 神经递质等)和生物大分子(磷酸二酯酶)的结构及其相互作用。全文由三章构成： 第一章简单介绍了溶液中从头算理论的最新进展以及磷酸二酯酶体系活性中心 的结构；第二章采用几种溶剂化模型计算了一些生物小分子在溶液中的离解常 数，通过溶液中频率的计算分析了组
2、胺在溶液中的主要构象；第三章运用量子 力学/分子力学和分子动力学模拟方法研究磷酸二酯酶-5(PDE5)活性中心的结构， 利用分子对接和 3D-QSAR 方法研究了环鸟嘌呤衍生物同磷酸二酯酶-5(PDE5)之 间的相互作用。 第二章的第一部分运用第一原理电子结构方法对 24 种胺类 化合物，包括可卡因、尼古丁、10 种神经递质和 12 种苯胺分子进行 pKlt;,agt;值的计算，采用了四种不同的自洽反应场溶剂化模型： 表面极化和体极化作用(SVPE)的方法、极化连续介质模型(PCM)、积分连续介质 模型(IEFPCM)、极化导体模型(COSMO，亦称为 CPCM)。由 SVPE 方法计算的绝对
3、 pKlt;,agt;值同实验值最接近，计算值和实验值之间的均方根差 (RMSD)只有 1.18，而用 PCM、IEFPCM 和 COSMO 方法计算的 pKlt;,agt;值和实验值之间的 RMSD 分别为 3.21、2.72 和 3.08。当 所有溶剂化模型所计算的 pKlt;,agt;值同实验值进行一元线性拟合 后，RMSD 都变小了，在 0.510.83 之间。最小 RMSD 的值(0.51)也是对应着表 面极化和体极化作用(SVPE)的方法。所有的计算结果都表明运用 SVPE 模型的第 一原理电子结构计算是预测胺类化合物比较合理的方法。第二部分利用第一原 理电子结构方法计算了一价和二
4、价组胺阳离子在气相和溶液中各种可能的构象， 并且用 IEFPCM 模型对组胺一价阳离子进行了频率计算，根据校正的计算频率以 及同位素的移动同实验的红外和拉曼光谱进行比较。最后的结果表明， glt;,3gt;H 构象在溶液中的自由能最低，溶液中约占 57的比率； 根据计算频率和实验光谱的比较，glt;,3gt;H 构象的计算频率和实 验光谱吻合的较好，平均偏差约为 10 cmlt;#39;-1gt;。所有 结果都表明：在溶液中侧链质子化的一价组胺阳离子在溶液中主要是 glt;,3gt;H 构象，而不是 tlt;,3gt;H 构象。这里所得 到的溶液中侧链质子化的一价组胺离子详细的结构信息对将来研
5、究组胺分子同 各种生物大分子体系的相互作用具有重要的指导意义。 第三章的第一部分通 过量子力学/分子力学(QM/MM)和分子动力学模拟计算得到了两种 PDE5 酶活性位 点的结构。为了验证：PDE5 活性中心第二个桥配体(BL2)的性质：到底是氢氧 根 HOlt;#39;-gt;还是 Hlt;,2gt;O 水分子?分别从 X-射线的晶体结构和分子动力学模拟平衡后结构出发进行 QM/MM 优化计算。计 算结果清楚的表明只有当 BL2 是氢氧根离子时，它才能真正的充当桥配体连接 两个带正电的金属离子 Znlt;#39;2+gt;和 Mglt;#39;2+gt;；当 BL2 是水分子的时候，它只能同
6、 MR2+离 子配位而离开 Znlt;#39;2+gt;离子。而且计算结果表明 QM/MM优化的几何构型受到溶剂水分子、蛋白质的动力学环境以及 embed MM 部分电荷 的影响，而且这几个因素是相互依赖的。当完全忽略或完全考虑这三种因素的 QM/MM 计算结果是一致的，表明这三种因素对几何构型优化的影响可以相互抵 消。然而，如果仅仅考虑其中的几个，就会得到完全不同的结果。因此要么全 部考虑，要么全部忽略。第二部分根据得到的 PDE5 活性中心的结构，运用分子 对接和 3D-QSAR 分析来研究 PDE5 和 PDE6 酶同一系列环鸟嘌呤衍生物的相互作 用。通过分子对接得到化合物的初始构象，运
7、用 CoMFA 和 CoMSIA 分析得到定量 构效关系式和定量结构选择性关系式模型(具有比较高的交叉验证系数 qlt;#39;2gt;和传统的相关系数 rlt;#39;2gt;的值)，这些关系式能够预测 PDE5 抑制剂的活性 和 PDE5 同 PDE6 之间选择性的抑制剂。这些比较高的验证系数 qlt;#39;2gt;和 rlt;#39;2gt;的值以及后面 所做的进一步的测试都表明所得的 3D-QSAR 和 3D-QSSR 模型是合理的，这些关 系式在提高环鸟嘌呤衍生物对靶标蛋白质的抑制性和选择性方面具有很高的应 用价值。正文内容正文内容本文运用一系列的理论化学方法研究了一些生物小分子(
8、尼古丁、可卡因、 神经递质等)和生物大分子(磷酸二酯酶)的结构及其相互作用。全文由三章构成： 第一章简单介绍了溶液中从头算理论的最新进展以及磷酸二酯酶体系活性中心 的结构；第二章采用几种溶剂化模型计算了一些生物小分子在溶液中的离解常 数，通过溶液中频率的计算分析了组胺在溶液中的主要构象；第三章运用量子 力学/分子力学和分子动力学模拟方法研究磷酸二酯酶-5(PDE5)活性中心的结构， 利用分子对接和 3D-QSAR 方法研究了环鸟嘌呤衍生物同磷酸二酯酶-5(PDE5)之 间的相互作用。 第二章的第一部分运用第一原理电子结构方法对 24 种胺类 化合物，包括可卡因、尼古丁、10 种神经递质和 12
9、 种苯胺分子进行 pKlt;,agt;值的计算，采用了四种不同的自洽反应场溶剂化模型： 表面极化和体极化作用(SVPE)的方法、极化连续介质模型(PCM)、积分连续介质 模型(IEFPCM)、极化导体模型(COSMO，亦称为 CPCM)。由 SVPE 方法计算的绝对 pKlt;,agt;值同实验值最接近，计算值和实验值之间的均方根差 (RMSD)只有 1.18，而用 PCM、IEFPCM 和 COSMO 方法计算的 pKlt;,agt;值和实验值之间的 RMSD 分别为 3.21、2.72 和 3.08。当 所有溶剂化模型所计算的 pKlt;,agt;值同实验值进行一元线性拟合 后，RMSD
10、都变小了，在 0.510.83 之间。最小 RMSD 的值(0.51)也是对应着表 面极化和体极化作用(SVPE)的方法。所有的计算结果都表明运用 SVPE 模型的第 一原理电子结构计算是预测胺类化合物比较合理的方法。第二部分利用第一原 理电子结构方法计算了一价和二价组胺阳离子在气相和溶液中各种可能的构象， 并且用 IEFPCM 模型对组胺一价阳离子进行了频率计算，根据校正的计算频率以 及同位素的移动同实验的红外和拉曼光谱进行比较。最后的结果表明， glt;,3gt;H 构象在溶液中的自由能最低，溶液中约占 57的比率； 根据计算频率和实验光谱的比较，glt;,3gt;H 构象的计算频率和实
11、验光谱吻合的较好，平均偏差约为 10 cmlt;#39;-1gt;。所有 结果都表明：在溶液中侧链质子化的一价组胺阳离子在溶液中主要是 glt;,3gt;H 构象，而不是 tlt;,3gt;H 构象。这里所得 到的溶液中侧链质子化的一价组胺离子详细的结构信息对将来研究组胺分子同 各种生物大分子体系的相互作用具有重要的指导意义。 第三章的第一部分通 过量子力学/分子力学(QM/MM)和分子动力学模拟计算得到了两种 PDE5 酶活性位 点的结构。为了验证：PDE5 活性中心第二个桥配体(BL2)的性质：到底是氢氧 根 HOlt;#39;-gt;还是 Hlt;,2gt;O 水分子?分别从 X-射线的
12、晶体结构和分子动力学模拟平衡后结构出发进行 QM/MM 优化计算。计 算结果清楚的表明只有当 BL2 是氢氧根离子时，它才能真正的充当桥配体连接 两个带正电的金属离子 Znlt;#39;2+gt;和 Mglt;#39;2+gt;；当 BL2 是水分子的时候，它只能同 MR2+离 子配位而离开 Znlt;#39;2+gt;离子。而且计算结果表明 QM/MM 优化的几何构型受到溶剂水分子、蛋白质的动力学环境以及 embed MM 部分电荷 的影响，而且这几个因素是相互依赖的。当完全忽略或完全考虑这三种因素的 QM/MM 计算结果是一致的，表明这三种因素对几何构型优化的影响可以相互抵 消。然而，如果
13、仅仅考虑其中的几个，就会得到完全不同的结果。因此要么全 部考虑，要么全部忽略。第二部分根据得到的 PDE5 活性中心的结构，运用分子对接和 3D-QSAR 分析来研究 PDE5 和 PDE6 酶同一系列环鸟嘌呤衍生物的相互作 用。通过分子对接得到化合物的初始构象，运用 CoMFA 和 CoMSIA 分析得到定量 构效关系式和定量结构选择性关系式模型(具有比较高的交叉验证系数 qlt;#39;2gt;和传统的相关系数 rlt;#39;2gt;的值)，这些关系式能够预测 PDE5 抑制剂的活性 和 PDE5 同 PDE6 之间选择性的抑制剂。这些比较高的验证系数 qlt;#39;2gt;和 rlt
14、;#39;2gt;的值以及后面 所做的进一步的测试都表明所得的 3D-QSAR 和 3D-QSSR 模型是合理的，这些关 系式在提高环鸟嘌呤衍生物对靶标蛋白质的抑制性和选择性方面具有很高的应 用价值。 本文运用一系列的理论化学方法研究了一些生物小分子(尼古丁、可卡因、神经 递质等)和生物大分子(磷酸二酯酶)的结构及其相互作用。全文由三章构成：第 一章简单介绍了溶液中从头算理论的最新进展以及磷酸二酯酶体系活性中心的 结构；第二章采用几种溶剂化模型计算了一些生物小分子在溶液中的离解常数， 通过溶液中频率的计算分析了组胺在溶液中的主要构象；第三章运用量子力学/ 分子力学和分子动力学模拟方法研究磷酸二
15、酯酶-5(PDE5)活性中心的结构，利 用分子对接和 3D-QSAR 方法研究了环鸟嘌呤衍生物同磷酸二酯酶-5(PDE5)之间 的相互作用。 第二章的第一部分运用第一原理电子结构方法对 24 种胺类化 合物，包括可卡因、尼古丁、10 种神经递质和 12 种苯胺分子进行 pKlt;,agt;值的计算，采用了四种不同的自洽反应场溶剂化模型： 表面极化和体极化作用(SVPE)的方法、极化连续介质模型(PCM)、积分连续介质 模型(IEFPCM)、极化导体模型(COSMO，亦称为 CPCM)。由 SVPE 方法计算的绝对 pKlt;,agt;值同实验值最接近，计算值和实验值之间的均方根差 (RMSD)
16、只有 1.18，而用 PCM、IEFPCM 和 COSMO 方法计算的 pKlt;,agt;值和实验值之间的 RMSD 分别为 3.21、2.72 和 3.08。当 所有溶剂化模型所计算的 pKlt;,agt;值同实验值进行一元线性拟合 后，RMSD 都变小了，在 0.510.83 之间。最小 RMSD 的值(0.51)也是对应着表 面极化和体极化作用(SVPE)的方法。所有的计算结果都表明运用 SVPE 模型的第 一原理电子结构计算是预测胺类化合物比较合理的方法。第二部分利用第一原 理电子结构方法计算了一价和二价组胺阳离子在气相和溶液中各种可能的构象， 并且用 IEFPCM 模型对组胺一价阳离子进行了频率计算，根据校正的计算频率以 及同位素的移动同实验的红外和拉曼光谱进行比较。最后的结果表明， glt;,3gt;H 构象在溶液中的自由能最低，溶液中约占 57的比率； 根据计算频率和实验光谱的比较，glt;,3gt;H 构象的计算频率和实 验光谱吻合的较好，平均偏差约为 10 cmlt;#39;-1gt;。所有 结果都表明：在溶液中侧
