靳玲俠，閔鎖田,王志銀,趙夢婷，任宏江

(1.陜西理工學(xué)院 化學(xué)與環(huán)境科學(xué)學(xué)院,陜西 漢中 723001; 2.西安文理學(xué)院 化學(xué)工程學(xué)院,西安 710065)

?

5-羥甲基胞嘧啶pKa值的理論研究

靳玲俠1，閔鎖田1,王志銀1,趙夢婷1，任宏江2

(1.陜西理工學(xué)院 化學(xué)與環(huán)境科學(xué)學(xué)院,陜西 漢中 723001; 2.西安文理學(xué)院 化學(xué)工程學(xué)院,西安 710065)

采用兩種不同動力學(xué)循環(huán)方案，在B3LYP/6-311++G(d,p)+ZPE和B3LYP/aug-cc-pVTZ//B3LYP/6-311++G(d,p)+ZPE水平上，同時考慮PCM和CPCM模型對5-羥甲基胞嘧啶(5-HMeCyt)中N3、N4、O2不同位點質(zhì)子化的pKa值進(jìn)行計算研究.結(jié)果表明：采用經(jīng)驗值法，5-HMeCyt中N3位置質(zhì)子化的pKa值最大；采用質(zhì)子交換法，以N3質(zhì)子化的胞嘧啶為參考酸，發(fā)現(xiàn)CPCM模型計算結(jié)果對計算方法的依賴性較大.

B3LYP；5-甲基胞嘧啶；胞嘧啶；5-羥甲基胞嘧啶

5-羥甲基胞嘧啶(5-HMeCyt)可能是一個重要的表觀遺傳學(xué)標(biāo)記，也可能是5-甲基胞嘧啶(5-MeCyt)去甲基向胞嘧啶(Cyt)轉(zhuǎn)化的中間產(chǎn)物.已有文獻(xiàn)報道，pKa值對氫原子轉(zhuǎn)移、質(zhì)子交換、電子轉(zhuǎn)移等反應(yīng)機制的闡明尤為重要[1-4].Arnott等[5]實驗發(fā)現(xiàn)在一定pH值下，胞嘧啶的N3位質(zhì)子可轉(zhuǎn)移至鳥嘌呤.同時，胞嘧啶的堿基序列對pH值會產(chǎn)生一定影響.因此，可通過pKa值控制堿基配對的最佳條件.另外，在酸性介質(zhì)中，Cyt存在3個可能的質(zhì)子化位置，分別為N3、O2和N4.綜合文獻(xiàn)報道[6-7]，Cyt在酸性條件下形成6個質(zhì)子化異構(gòu)體，用滴定方法較難準(zhǔn)確測量有機物具體位點的pKa值.因此采用理論與計算化學(xué)估算嘧啶堿基的pKa值對于理解以質(zhì)子轉(zhuǎn)移為基礎(chǔ)的化學(xué)物質(zhì)和生物進(jìn)程具有非常重要的意義.

在實驗方面，F(xiàn)loria′n等人在研究胞嘧啶質(zhì)子化異構(gòu)體紅外光譜時發(fā)現(xiàn)，胞嘧啶的N3，N4，O2位置上均可發(fā)生質(zhì)子化反應(yīng)，同時發(fā)現(xiàn)N3和O2質(zhì)子化異構(gòu)體較其它異構(gòu)體穩(wěn)定[8].Saenger等人通過實驗檢測到胞嘧啶的pKa值為4.5[9-10].由于胞嘧啶質(zhì)子化位置的復(fù)雜性且部分異構(gòu)體的不穩(wěn)定性，導(dǎo)致實驗難于準(zhǔn)確檢測不穩(wěn)定異構(gòu)體的pKa值.在理論研究方面，Moser等人[11]采用密度泛函理論方法得出胞嘧啶的pKa值為6.29～6.98，并預(yù)測甲基化對胞嘧啶pKa值的影響，而George等采用Poisson-Boltzmann方法對胞嘧啶在第三螺旋結(jié)構(gòu)的N3位置pKa值進(jìn)了計算，得出N3質(zhì)子化胞嘧啶的pKa值為4.3[12]，與實驗值極為接近(4.5)，5-HMeCyt是一種重要的表觀遺傳修飾，可能與去甲基化過程有關(guān)，但其作用機制仍不明確.鑒于此，本文將系統(tǒng)考察5-HMeCyt不同位點質(zhì)子化的pKa值(圖1)，同時探討PCM和CPCM對pKa值的影響，解決質(zhì)子化位置是否起顯著催化作用.

圖1 5-HMeCyt及其N3、O2、N4位質(zhì)子化結(jié)構(gòu)示意圖

1 計算方法

為比較計算方法對計算結(jié)果的影響，本文選用B3LYP/6-311++G(d,p)+ZPE和B3LYP/aug-cc-pVTZ//B3LYP/6-31G(d,p)+ZPE兩種計算方法[13-14]，并且使用兩種循環(huán)方案分別對5-HMeCyt不同位點質(zhì)子化的pKa值進(jìn)行計算.以氣相分子作為研究對象，應(yīng)用Gaussian03程序包[15]，采用密度泛函理論B3LYP方法對5-HMeCyt及其質(zhì)子化異構(gòu)體的初始幾何構(gòu)型進(jìn)行優(yōu)化和能量計算，并采用頻率進(jìn)行驗證，以確定其最穩(wěn)定的異構(gòu)體結(jié)構(gòu).同時考察CPM和CPCM模型對5-HMeCyt的質(zhì)子化異構(gòu)體pKa值的影響，對各異構(gòu)體的溶劑化自由能進(jìn)行了計算[16].

1.1 經(jīng)驗值法[17]

Scheme1 經(jīng)驗值法:

↑-ΔGs(C5H8N3O+)↓ΔGs(C5H7N3O2)↓ΔGs(H+)

Scheme2 質(zhì)子交換法:

↑-ΔGs(C5H8N3O+)↑-ΔGs(C4H5N3O)↓ΔGs(C5H7N3O2)↓ΔGs(C4H6N3O+)

Scheme3采用經(jīng)驗值法計算5-HMeCyt質(zhì)子化異構(gòu)體的pKa值:

(1)

pKa=ΔG(aq)/(2.303RT)

(2)

Scheme4采用質(zhì)子交換法計算5-HMeCyt質(zhì)子化異構(gòu)體的pKa值:

(3)

pKa=ΔG(aq)/(2.303RT)+pKa(C4H6N3O+)

(4)

1.2 質(zhì)子交換法[18]

表1 氣相自由能Gg、吉布斯自由能校正值GZEP、自由能在水溶液中變化值ΔGaq(a.u.)及pKa值

注:aCPCM and PCM,B3LYP/6-311++G(d,p)+ZPE；SP-CPCM and SP-PCM,B3LYP/aug-cc-pVTZ//B3LYP/6-31G(d,p)+ZPE.

2 結(jié)果與討論

本文利用上述兩種熱力學(xué)循環(huán)方法計算出質(zhì)子化5-HMeCyt在水溶液中的pKa值.

比較經(jīng)驗值法采用不同模型計算出5-HMeCyt不同位點質(zhì)子化異構(gòu)體的pKa值，發(fā)現(xiàn)不同位點質(zhì)子化異構(gòu)體的pKa值由N4、O2、N3依次增大，說明5-HMeCyt不同位點質(zhì)子化異構(gòu)體的pKa值存在顯著差別.在PCM模型下，5-HMeCyt對計算方法的依賴較小.在CPCM模型下采用B3LYP/6-311++G(d,p)+ZPE與B3LYP/aug-cc-pVTZ方法計算的結(jié)果偏差較大.

采用質(zhì)子交換法運用不同計算方法對5-HMeCyt的N3位質(zhì)子化異構(gòu)體pKa值的計算結(jié)果表明，在相同計算方法下采用CPCM和PCM模型計算結(jié)果存在較大偏差.

3 結(jié)論

本文通過密度泛函理論、溶劑化模型和兩種熱力學(xué)循環(huán)方案對5-HMeCyt不同位點的pKa值進(jìn)行了詳細(xì)計算研究，得出以下結(jié)論：(1)5-HMeCyt不同位點的pKa值具有顯著差別；(2)在PCM模型下5-HMeCyt不同位點pKa值對計算方法的依賴性較小，而CPCM模型對計算方法依賴性較大.

[1] LATASHA M S,GALINA M C.Theoretical study of gas phase tautom erization rea ctions for the ground and first excited electronic states of adenine[J].J Phys.Chem.A,2002,106：4251-4256.

[2] HAYATSU,H.Bisulfite modification of cytosine and 5-methylcytosine as used in epigenetic studies[J].Genes and Environment，2006,28：1-8.

[3] HAYATSU H.Discovery of bisulfite-mediated cytosine conversion to uracil,the key reaction for DNA methylation analysis —A personal account[J].Proceedings of the Japan Academy,Series B，2008,84(8):321-330.

[4] HUANG Y,PASTOR W A,SHEN Y,et al.The behaviour of 5-hydroxymethylcytosine in bisulfite sequencing[J] .Plosone，2010,5(1)：8888.

[5] MOTORIN Y,LYKO F,HELM M.5-Methylcytosine in RNA：detection,enzymatic formation and biological functions[J].Nucleic Acids Research，2010,38 (5)：1415-1430.

[7] PURRELLO R,MOLINA M,WANG Y，et al.Keto-iminol tautomerism of protonated cytidine monophosphate characterized by ultraviolet resonance raman spectroscopy：implications of C+ iminol tautomer for base mispairing[J].J.Am.Chem.Soc.,1993,115(2)：760-767.

[8] CHEN Z Q,ZHANG C H,KIM C K,et al.Quantum mechanics study and monte carlo simulation on the hydrolytic deamination of 5-methylcytosine glycol[J].J.Phys.Chem,2011,13：6471-6483.

[9] LABET V,MORELL C,CADET J,et al.A.Proton catalyzed hydrolytic deamination of cytosine：a computational study[J].Theor Chem Acc,2008,120(4-6):429-435.

[10]LABET V,MORELL C,CADET J,et al.Hydrolytic deamination of 5-methylcytosine in protic medium：a theoretical study[J].J.Phys.Chem.A,2009,113 (11)：2524-2533.

[11]RANGE K, LóPEZ C S, MOSER A, et al. Multilevel and density functional electronic structure calculations of proton affinities and gas-phase basicities involved in biological phosphoryl transfer[J]. J. Phys. Chem. A, 2006, 110(2): 791-797.

[12]SAENGER W.Principles of Nucleic Acid Structure[M].New York：Springer Verlag,1984.

[13]BECKE A D.Density functional thermochemistry.III.The role of exact exchange[J].J.Chem.Phys.1993,98：5648-5652.

[14]LEE C,YANG W,PARR R G.Development of the colic-salvetti correlation-energy formula into a functional of the electron density[J].Phys.Rev.B.，1988,37：785-789.

[15]FRISCH M J,SCHLEGEL H B,SCUSERIA G E,et al.Gaussian 03，Revision 02[S].Wallingford CT,Gaussian,Inc，2004.

[16]MIERTU? S,TOMASI J.Approximate evaluations of the electrostatic free energy and internal energy changes in solution processes[J].Chem.Phys, 1982,65(2)：239-245.

[17]COSSI M,BARONE V,CAMMI R,et al.Ab-initio study of solvated molecules：a new implementation of the polarizable continuum model[J].Chem.Phys.Lett,1996,255：327-335.

[18]GEORGE R P,WONG L,LAMM G.pKa of cytosine on the third strand of a triplex DNA：preliminary poisson-boltzmann calculations[J].Int.J.Quant.Chem.,1998,70：1177-1184.

[責(zé)任編輯 馬云彤]

Theoretical Research on pKa Value of the 5-Hydroxymethylcytosine

JIN Ling-xia1, MIN Suo-tian1, WANG Zhi-yin1, ZHAO Meng-ting1, REN Hong-jiang2

(1. College of Chemistry and Environmental Science, Shaanxi University of Technology,Hanzhong 723001, China; 2. School of Chemical Engineering, Xi’an University, Xi’an 710065, China)

The protonated pKa values in different positions of O2, N3, N4 for 5-hydroxymethylcytosine (5-HMeCyt) are calculated by adopting two different dynamic cycle schemes at the B3LYP/6-311++G(d, p)+ZPE and B3LYP/aug-cc-pVTZ//B3LYP/6-311++G(d, p)+ZPE level, which consider the effect of the PCM and CPCM models. The results show that the protonated pKa value is the largest in the N3 position of the 5-HMeCyt by the empirical value method. The results of CPCM model greatly depend on the calculation method when taking the N3 position of cytosine as reference by the proton exchange method.

B3LYP; the 5-methylcytocine; cytosine; the 5-hydroxymethylcytosine

1008-5564(2016)04-0053-04

2016-03-24

陜西理工學(xué)院科研基金資助項目(SLGQD14-10)

靳玲俠(1981—)，女，陜西大荔人，陜西理工學(xué)院化學(xué)與環(huán)境科學(xué)學(xué)院講師，博士，主要從事生物小分子反應(yīng)機理的理論研究.

O641

A