最新版金属离子与核酸碱基及水作用模式的理论探讨毕业设计.docx
- 文档编号:24469126
- 上传时间:2023-05-27
- 格式:DOCX
- 页数:36
- 大小:1.61MB
最新版金属离子与核酸碱基及水作用模式的理论探讨毕业设计.docx
《最新版金属离子与核酸碱基及水作用模式的理论探讨毕业设计.docx》由会员分享,可在线阅读,更多相关《最新版金属离子与核酸碱基及水作用模式的理论探讨毕业设计.docx(36页珍藏版)》请在冰豆网上搜索。
最新版金属离子与核酸碱基及水作用模式的理论探讨毕业设计
LiaoningNormalUniversity
(2014届)
本科生毕业论文(设计)
题目:
金属离子与核酸碱基及水作用模式的理论探讨
学院:
化学化工学院
专业:
化学
班级序号:
2班23号
学生姓名:
张丹
指导教师:
宫利东
2014年5月
目录
摘要(Abstract)1
关键词(Keywords)1
前言2
1理论和计算方法2
1.1量子化学方法简介2
1.2Gaussian09简介3
2计算方法4
3.结果与讨论4
3.1碱基单体的结构和电荷分布4
3.2金属离子与核酸碱基的复合物5
3.2.1结合能7
3.2.2电荷分布8
3.3金属离子与核酸碱基及水形成复合物Mn+-Base-(H2O)n的结构和结合能11
3.3.1Mn+-Base-(H2O)1体系构型11
3.3.2Mn+-Base-(H2O)2体系构型12
3.3.3Mn+-Base-(H2O)n体系结合能14
3.3.4Mn+-Base-(H2O)n体系的电荷分析16
参考文献:
17
致谢18
金属离子与核酸碱基及水作用模式的理论探讨
摘要:
金属离子与核酸碱基的相互作用对许多生理过程都会产生影响。
不同的金属离子与核酸碱基相互作用的强弱取决于金属离子所带电荷以及核酸碱基的结构和不同位点。
本文应用量子化学方法研究了Li+、Na+、K+、Be2+、Mg2+、Ca2+与核酸碱基(A、C、G、T、U)及水的相互作用,包括复合物的结合位点、体系能量、电荷分布等。
采用了B3LYP/6-311+G(d,p)、B3PW91/6-311+G(d,p)和M062X/6-311+G(d,p)三种方法对复合物进行了几何构型的优化和能量的计算,并且在HF/STO-3G水平上对几何构型进行了Mulliken电荷布局分析,并将得到的结果进行比较得出结论。
关键词:
金属离子;核酸碱基;量子化学计算
Abstract:
Theinteractionofmetalionswithnucleicacidbasesimpactsonmanyphysiologicalprocesses.Theinteractionbetweendifferentmetalionsandnucleicacid basesdependsonthemetalioncharges,thestructureofthenucleicacidbasesanddifferentregions.Inthiswork,theinteractionamongLi+/Na+/K+/Be2+/Mg2+/Ca2+,nucleicacidbasesandwaterhasbeenstudiedbyquantumchemicalcalculationmethods,includingregions,bindingenergyandchargedistributions.WeoptimizedthecomplexesandcalculatedtheenergyofthecomplexeswiththemethodsofB3LYP/6-311+G(d,p)、B3PW91/6-311+G(d,p)andM062X/6-311+G(d,p),andanalyzedthedistributionofMullikenchargesattheHF/STO-3Glevel,andcomparedtheobtainedresultstoreachaconclusion.
Keywords:
metalions;nucleicacidbases;quantumchemicalcalculation
前言
核酸是一种生物大分子化合物,由许多核苷酸聚合而成,存在范围广泛,所有微生物、植物细胞、动物、生物体体内都含有核酸,是生命所必需的基本物质之一。
由于核苷酸的排列顺序及其化学成分的不同,又组成不同的核酸。
根据其化学成分的不同,将核酸分为脱氧核糖核酸(DNA)和核糖核酸(RNA),它们都是由碱基、戊糖和磷酸组成。
核酸是生物承载和传递信息的载体,在生物体的结构、功能、甚至生物催化上都有重要的作用,而这些重要作用的发挥往往离不开金属离子。
金属离子与核酸碱基相互作用能够维持核酸碱基的结构和功能,这是由于磷酸骨架带有负电荷,金属离子能中和其负电,而起到稳定结构的作用,虽然碱金属离子通常与磷酸基团相互作用,但与碱基的结合会比与磷酸基团相互作用产生更为显著的影响[1]。
金属离子-核酸碱基的相互作用参与了许多生理过程:
如维持三倍体、四倍体DNA螺旋的稳定性以及Z-DNA的稳定性;一些金属离子可以扰乱DNA的复制过程[2]。
金属离子可以与碱基直接结合,也可以通过一个水分子而相互结合[3],由于金属离子的加入可使碱基对的稳定性得到加强从而稳定和调控多种核酸碱基。
利用金属离子对核酸的影响,可以帮助人们开发研究抗癌药物,例如人们利用顺铂与鸟嘌呤相互作用,研发出可以促发癌细胞死亡的药物。
而甲基化是蛋白质和核酸的一种重要的修饰,调节基因的表达和关闭,与癌症、衰老、老年痴呆等许多疾病密切相关,是表观遗传学的重要研究内容之一。
最常见的甲基化修饰有DNA甲基化和组蛋白甲基化。
在生物体系中,碱金属离子参与多种生理活动,在维持生命的生化过程中起到了重要的作用。
Rodgers等[8]探讨了最稳定的腺嘌呤异构体和碱金属离子(Li+、Na+)的结合情况,发现金属离子在N6、N7位置上以双齿形式配位,要比在Nl或N3位置上单齿配位得到的复合物稳定,并且这种双齿形式配位方式能够破坏碱基对A--T中的氢键。
Russo等[9]的研究结果也证实Na+、K+在N6、N7位置上以双齿形式配位得到的复合物最稳定。
为了了解金属离子在核酸生化过程中的作用,研究碱金属离子与碱基的相互作用是十分必要的。
具体来说,碱金属离子与核酸结合会产生怎样的影响取决于碱金属离子与其相互作用位置。
金属离子与碱基单体相互作用具有不同的结合位点,在不同的位点其结合能也随之不同。
对于DNA双螺旋结构,在碱基上碱金属离子的最佳成键位置为嘌呤碱基的N7位,而鸟嘌呤-碱金属复合物的稳定能大于腺嘌呤-金属复合物的稳定能[10]。
Russo等人[11,12]在2001年报道了碱金属离子和腺嘌呤及鸟嘌呤的N7、N9和O6位能够发生相互作用
由核酸碱基结构的复杂性、易变性,且易受限于实验条件或者实验手段自身的缺陷,对于某些过程中具体的相互作用细节往往无法检测,因此实验研究方法就具有一定的局限性。
随着计算机科学技术的飞速发展,可用于理论模拟的方法和程序如雨后春笋般蓬勃地发展,理论模拟已经成为人们研究生物大分子体系的一种有效手段。
在理论方法中,从头计算方法不依赖任何经验参数,最近几年已经渗透到许多研究领域,在金属离子与碱基相互作用的研究中,越来越多的运用了高精度的量子化学计算方法,如:
MP2、CCSD(T)等。
本文应用量子化学方法研究金属离子与核酸碱基及水的相互作用,计算其键长、结合能,同时分析电荷转移,为建立含金属离子体系的核酸分子力场提供参考。
1理论和计算方法
1.1量子化学方法简介
量子力学是描述微观体系运动规律的科学,量子化学从1927年Heitler和London研究氢分子的结构开始,七十多年来已发展成为一门独立的学科,同时也与化学各分支学科以及物理、生物、计算数学等互相渗透的学科,在多个领域如材料、能源、环境、化工生产以及激光技术中得到广泛应用。
从头计算(abinitio)方法,主要是近十年内发展起来的原子体系和分子体系的量子理论计算方法。
而所谓从头计算法的“头”,就是因为它只以分子轨道理论在物理模型上的三个基本近似作为出发点:
(1)采用非相对论的量子理论,即从Schrödinger方程出发;
(2)使用Born-Oppenheimer近似;(3)轨道近似,即单电子近似。
从这三点基本近似出发,采用原子轨道线性组合成分子轨道(LCAO-MO)方法,除此之外,不再引入其他近似,在计算中只利用原子序数Z、普朗克常量h、电子的质量me和电量e这四个基本物理常数,而不借助于任何经验参数。
从头计算法实质上是以分子轨道理论为基础的,选取适当的基函数,计算各种所需要的积分,利用自洽场方法求解Hartree-Fock-Roothaan方程,得到分子或其他体系的分子轨道、轨道能和波函数,进而获得体系的其他性质,如平衡几何构型、电荷密度分布、偶极矩、振动频率以及势能面等。
1965年,Pople开发出Gaussian70程序,它是量子化学领域最著名和应用最广泛的软件之一,可以应用从头计算方法、半经验计算方法等进行分子能量和结构;过渡态能量和结构;化学键及反应能量;分子轨道;偶极矩;多极矩;红外光谱和拉曼光谱,核磁共振,极化率和超极化率,热力学性质,反应路径等分子相关计算。
此后又将其发展到Gaussian98程序,目前最新版本为Gaussian09。
但是从头计算采用单组态近似,即不考虑组态相互作用,这会产生一定的误差,因为对于一对自旋相同的电子(α-α或β-β),它们的相关能均包含在库仑积分和交换积分之中;但对于自旋相反的电子,交换积分为零,因此对于相关能就考虑不够。
在考虑组态相互作用后,计算结果将得到改善,使体系的能量进一步下降。
1.2Gaussian09简介
Gaussian是一个功能强大的量子化学综合软件包,广泛应用于量子化学领域。
其可执行程序可在不同型号的大型计算机,超级计算机,工作站和个人计算机上运行,并相应有不同的版本。
Gaussian功能主要有:
分子能量和结构、过渡态能量和结构、键和反应能量、分子轨道多重矩、原子电荷和电势、振动频率、红外和拉曼光谱、核磁性质、极化率和超极化率、热力学性质、反应路径等。
Gaussian09是Gaussian系列电子结构程序的较为新的版本。
它在化学、化工、生物化学、物理化学等化学相关领域方面的功能都进行了增强。
其主要功用大体有以下几个大的方面:
(1)研究大分子的反应和光谱
(2)通过自旋-自旋耦合常数确定构像
(3)研究周期性体系
(4)预测光谱
(5)模拟在反应和分子特性中溶剂的影响
Gaussian09软件内包括的基组关键词为:
STO-3G,3-21G,6-31G,6-31+G*,6-311G,极化函数:
6-31G(d)和6-31G(d,p).Gaussian03规定,单个的第一个极化函数可以用*或**表示,同时,(d,p)和**代表同一函数,比如6-31G**和6-31G(d,p)等价;3-21G*基组代表在第二排原子上附加了极化函数。
还有些基组除了定义了多极化函数以外,还附加了弥散函数,用+和++代表之。
以6-31+G(3df,2p)基组为例来说明其中的含义,首先指明了一个附有弥散函数的6-31G基组,然后在重原子上加上三组d函数和一组f函数,并且在氢原子上再加上两组p函数,所有的这些因素加和起来就是整个的6-31+G(3df,2p)关键词代表的含义。
2计算方法
基组的选择也是获取准确信息的一个关键环节,其恰当与否将会直接影响到计算结果的精确程度.基组是对体系轨道的数学描述,对应着体系的波函数。
基组越大,所做的近似和限制越小,对轨道的描述也就越精确。
基组包括劈裂价键基组、极化基组和弥散基组。
众多的研究已经表明,对于含有氢键作用的相关体系,考虑弥散函数和极化函数是非常必要的。
劈裂价键基组是增加基组中基函数的数量,将价层电子的原子轨道用两个或两个以上基函数来表示。
弥散基组就是在劈裂价键的基础上添加了弥散函数的基组,允许轨道占据更大的空间,这样的基组可以用于非价键相互作用体系的计算。
6-311G(d,p)是在6-31G(d,p)基础上对重原子添加弥散函数。
极化基组就是在劈裂价键基组的基础上添加更高能级原子轨道所对应的基函数。
常用的极化基组对重原子添加d函数,对H原子添加p函数。
6-31G(2d,2p)是在6-31G(d,p)的基础上对重原子多添加一个d函数,对H原子多添加一个p函数。
本文分别采用B3LYP/6-311+G(d,p)、B3PW91/6-311+G(d,p)和M062X/6-311+G(d,p)三种方法对金属离子与碱基及水相互作用体系进行了几何构型优化和能量计算,在考察电荷分布时采用HF/STO-3G方法。
3.结果与讨论
3.1碱基单体的结构和电荷分布
AC
GT
U
图3.1五个碱基的结构示意图及电荷分布
由于碱基上的氮、氧原子都有孤对电子,带负电荷,即碱基上的氢原子带有正电荷,而金属离子带正电荷,所以当金属离子与碱基相互作用时,尽量避开带有正电荷的氢原子,与氮、氧原子间成键。
例如鸟嘌呤G,金属离子可与碱基中N4、O7结合成键,也可与N11、N14结合成键,共有两个位点。
3.2金属离子与核酸碱基的复合物
下列各图所示,为金属离子与不同核酸碱基在不同位点上结合所形成的复合物结构。
图3.2碱基与金属离子形成复合物结构及键长示意图
由于碱基上的氮、氧原子都有孤对电子,碱基上的氢原子带有正电荷,而金属离子带正电荷,所以当金属离子与碱基相互作用时,尽量避开带有正电荷的氢原子,形成如上图所示结构。
除胞嘧啶C外,A、G、T、U都遵守上述规律形成稳定构象。
对胞嘧啶C,由于其结构中N、O原子排列的特殊性以及甲基所带来的空间位阻,导致实验中的六种金属离子不能够全部与之有三个结合位点,因此,单电荷的金属离子只能与N10、O12共同成键,其中Li+由于核电荷数小、离子半径小等原因,有两个结合位点。
带两个电荷的金属离子不但能与N10、O12共同成键,还能与N7、N10共同结合成键,其中Be2+较为特殊,还可以与O12单独成键,原因可能是因为O12所带负电荷较多,电负性较大,而Be2+所带正电荷较大,足以形成较稳定的化学键。
由图中键长数据可知,在同一主族中,随着原子序数增加,金属离子与碱基形成的复合物中,金属离子到碱基上同一原子的距离也随之增大,同一周期中,随着金属离子电荷数增大,金属离子到碱基上同一原子的距离随之减小。
3.2.1结合能
对于金属离子与甲基取代的碱基相互作用后形成的复合物,结合能(ΔE)为金属离子的能量(EMn+)与碱基的能量(Ebase)之和再减去复合物的能量(EMn+-base),表达式为:
ΔE=EMn++Ebase-EMn+-base
其中,Ebase=EZPVE+ES;EMn+=ES;EMn+-base=EZPVE+ESP。
结合能的大小反应复合物的稳定性大小。
结合能越大,复合物越稳定,反之,稳定性小。
表1金属离子与碱基相互作用的结合能(单位:
kcal/mol)
Mn+-base
ΔE
Mn+-base
ΔE
B3LYP
B3PW91
M062X
MP2
B3LYP
B3PW91
M062X
MP2
A-Li-01
41.41
47.08
49.55
47.81
G-K-01
37.60
36.30
46.39
42.48
A-Li-02
44.98
43.45
44.75
41.06
G-K-02
13.04
12.16
15.49
14.59
A-Li-03
51.44
49.26
50.72
50.60
G-Be-01
322.88
318.97
318.95
304.83
A-Na-01
32.77
29.41
33.95
31.48
G-Be-02
240.22
237.14
236.00
224.95
A-Na-02
30.12
28.63
30.28
26.71
G-Mg-01
204.95
200.76
201.71
189.63
A-Na-03
34.08
32.02
35.50
33.01
G-Mg-02
141.67
137.20
137.99
133.37
A-K-01
19.46
21.22
24.33
23.75
G-Ca-01
154.67
152.93
152.78
140.94
A-K-02
20.13
19.40
21.24
19.57
G-Ca-02
93.37
91.79
92.03
87.75
A-K-03
22.25
21.23
24.91
24.34
T-Li-01
53.41
51.17
52.89
48.93
A-Be-01
279.03
275.05
274.92
266.42
T-Li-02
49.20
46.84
31.93
45.60
A-Be-02
250.68
247.77
246.48
235.12
T-Na-01
37.55
35.74
35.99
33.85
续表1
Mn+-base
ΔE
Mn+-base
ΔE
B3LYP
B3PW91
M062X
MP2
B3LYP
B3PW91
M062X
MP2
A-Be-03
285.74
281.54
281.38
275.41
T-Na-02
33.91
31.99
34.42
30.89
A-Mg-01
163.27
158.80
161.21
153.50
T-K-01
28.31
27.14
27.48
26.60
A-Mg-02
140.81
136.34
137.54
126.94
T-K-02
24.95
23.68
26.34
23.84
A-Mg-03
168.31
163.79
168.14
159.87
T-Be-01
248.29
242.52
245.54
237.11
A-Ca-01
112.99
111.43
111.95
106.69
T-Be-02
236.82
230.33
231.66
224.25
A-Ca-02
95.58
94.26
92.18
86.16
T-Mg-01
144.65
139.82
146.55
134.05
A-Ca-03
115.37
113.49
113.86
110.15
T-Mg-02
135.24
129.31
131.90
125.17
C-Li-02
69.97
68.03
70.22
64.99
T-Ca-01
111.71
109.63
108.33
100.92
C-Li-03
63.80
-
-
47.82
T-Ca-02
101.98
99.64
99.20
92.91
C-Na-02
52.50
50.54
53.48
47.68
U-Li-01
53.37
51.20
54.02
48.64
C-K-02
40.32
39.25
42.38
38.49
U-Li-02
46.92
44.56
47.89
43.49
C-Be-01
276.93
273.05
274.11
267.12
U-Na-01
37.81
35.98
39.25
34.07
C-Be-02
301.06
297.27
298.40
288.50
U-Na-02
32.00
30.08
33.78
29.36
C-Be-03
269.90
266.61
266.89
191.02
U-K-01
28.65
27.47
30.69
26.88
C-Mg-01
162.73
158.52
161.73
154.41
U-K-02
23.31
22.03
25.72
22.46
C-Mg-02
187.39
183.19
186.13
175.48
U-Be-01
247.15
241.40
243.80
234.11
C-Ca-01
115.36
113.87
114.89
109.36
U-Be-02
229.24
223.02
226.08
217.77
C-Ca-02
141.36
139.74
139.83
130.63
U-Mg-01
142.76
138.38
142.03
132.69
G-Li-01
74.70
71.21
74.22
68.52
U-Mg-02
128.63
123.60
128.08
120.05
G-Li-02
40.10
38.18
40.25
39.22
U-Ca-01
110.82
108.70
108.80
99.93
G-Na-01
50.50
46.84
50.38
51.93
U-Ca-02
96.93
94.53
95.59
88.41
G-Na-02
23.53
21.60
24.20
22.73
由表1中数据比较可知,同种金属离子与同种碱基在不同位点结合时,结合能是不同的。
对于复合物A-Be2+而言,A-Be2+-03结合能最大,A-Be2+-01结合能最小,而A-Be2+-02的结合能介于两者之间。
因此Be2+与腺嘌呤(A)作用最稳定的位点是A-Be2+-03。
由此法同样可以得到结论:
金属离子与胞嘧啶作用最稳定的位点是C-Mn+-02,与鸟嘌呤作用最稳定的位点是G-Mn+-01,与胸腺嘧啶作用最稳定的位点是T-Mn+-01,与尿嘧啶作用最稳定的位点是U-Mn+-01。
表2金属离子与尿嘧啶相互作用最稳定结构的结合能(单位:
kcal/mol)
B3LYP
B3PW91
M062X
MP2
B3LYP
B3PW91
M062X
MP2
U-Li-01
53.37
51.20
54.02
48.64
U-Be-01
247.15
241.40
243.80
234.11
U-Na-01
37.81
35.98
39.25
34.07
U-Mg-01
142.76
138.38
142.03
132.69
U-K-01
28.65
27.47
30.69
26.88
U-Ca-01
110.82
108.70
108.80
99.93
比较所采用的几种计算方法,可以看出由B3PW91方法计算所得到结合能最接近于MP2方法,即以此方法优化出的结构是最稳定的。
3.2.2电荷分布
表3复合物及腺嘌呤中的部分原子电荷(单位:
e)
Type
qN1
qN4
qN7
qN10
qN13
qMn+
Δq
Adenine
-0.2573
-0.2793
0.0302
-0.3189
-0.3025
A-Li-01
-0.2384
-0.2621
0.1091
-0.2592
-0.2368
0.5108
0.4892
A-Li-02
-0.2486
-0.2483
0.1355
-0.2
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 最新版 金属 离子 核酸 碱基 作用 模式 理论 探讨 毕业设计