王 瑩 易海波 李會吉 代 倩 曹治煒 路 洋

(湖南大學化學化工學院,化學生物傳感與計量學國家重點實驗,長沙410082)

鹽水溶液中離子與丙氨酸極性基團間的作用對丙氨酸締合的影響:密度泛函理論與分子動力學模擬

王 瑩 易海波*李會吉 代 倩 曹治煒 路 洋

(湖南大學化學化工學院,化學生物傳感與計量學國家重點實驗,長沙410082)

采用密度泛函理論和經(jīng)典分子動力學模擬研究了鹽水溶液中Na+、Cu2+、Zn2+、Cl-與丙氨酸分子間的相互作用對丙氨酸分子締合的影響.密度泛函理論的計算結果顯示丙氨酸分子與Na+、Cu2+、Zn2+、Cl-之間的相互作用可增強其電荷分離.經(jīng)典分子動力學模擬結果顯示在水溶液中兩性離子形式的丙氨酸存在三種締合結構.鹽水溶液中,陽離子、陰離子與丙氨酸間的相互作用均能一定程度上減弱丙氨酸分子的締合.但是陽離子與丙氨酸間的相互作用明顯受離子水合作用的影響.由于Cu2+水合作用較強,雖在氣相中Cu2+與丙氨酸分子之間相互作用明顯比Na+強,但是在水溶液中則情況剛好相反.在ZnCl2稀溶液中,Zn2+與丙氨酸間的相互作用被其第一水合殼層隔開.但這種相互作用仍能明顯影響丙氨酸分子的締合,這與Zn2+的水合殼層特征有關.另外,離子與丙氨酸之間的相互作用,不僅會削弱丙氨酸的締合,也可導致丙氨酸分子間的締合結構發(fā)生轉變.離子濃度也會影響其與丙氨酸分子間的締合形式以及丙氨酸的締合結構.

兩性離子;陽離子;水合;締合;分子間相互作用

1 引言

在所有已知的蛋白質中,幾乎有一半包含金屬輔助因子,其中生理學情況下鍵合到蛋白質上最常見的是Na+、K+、Mg2+、Ca2+、Zn2+、Cu2+、Fe3+、Co3+和Mn2+等離子.金屬離子廣泛參與生命過程,發(fā)揮各種重要的生物學功能.1-5目前臨床上難治愈的疾病大多數(shù)也與氨基酸以及一些金屬離子的異常結合有關,6-12但是一些藥理作用的機制問題仍然沒有解決.生命過程中氨基酸與金屬離子間的作用通常發(fā)生在水環(huán)境中,需考慮水合作用的影響.因此,探討水溶液中氨基酸與金屬離子之間的相互作用對進一步研究多肽、蛋白質,進而解釋生命現(xiàn)象具有重要的意義.另外,有研究也表明,氨基酸溶液中加入電解質后其溶解度明顯改變,不同陽離子和陰離子對氨基酸的分離效果有較大影響.13

中性條件下,氨基酸主要有兩性離子和中性分子兩種結構形式,在氣相和水溶液中氨基酸的結構存在明顯差異.氣相中,氨基酸主要以中性結構存在,附近電荷的存在可以穩(wěn)定氨基酸的兩性離子形式.14-17水溶液中,水分子優(yōu)先穩(wěn)定氨基酸的兩性離子形式結構.18-20因此,水溶液中電荷分離的氨基酸分子締合與水合易受到陽離子、陰離子的影響. Tomé21等以及Carta和Tola22的研究進一步表明,不同種類的陽離子和陰離子對氨基酸的影響也有所不同.鹽水溶液中陽離子、陰離子與氨基酸之間的相互作用雖然得到了較廣泛的研究,23-27但其作用機理仍沒有給出較明確的解釋.以往的工作主要側重于研究鹽水溶液中氨基酸溶解度的變化,但是很少有研究從離子水合、離子與氨基酸各極性基團之間相互作用的角度探討溶液中陽離子、陰離子對氨基酸的締合特征和溶解性的影響.氣相中,二價陽離子,如Cu2+、Zn2+與氨基酸羧基間的相互作用很強,但水溶液中二價陽離子的水合能力也很強.28,29在水溶液中,這兩種作用存在相互競爭,它們將如何影響氨基酸分子的締合呢?目前,尚無這方面的具體研究.另外,Harris30和Pye31等的研究表明在ZnCl2稀溶液中,Zn2+與Cl-間可形成較穩(wěn)定的溶劑共享離子對(SSIP)結構,說明Zn2+對處于其第二殼層的Cl-仍存在較明顯的束縛作用.在水溶液中,Zn2+的這種水合以及其與Cl-相互締合的特征是否也會影響氨基酸分子的締合?因此,有必要在考慮離子水合與締合特征的基礎上討論一些陽離子、陰離子對氨基酸分子間締合作用的影響.

本文采用量子化學計算方法研究丙氨酸(Ala)與離子間的相互作用以及離子的水合作用,并采用經(jīng)典分子動力學模擬研究不同濃度的NaCl、CuCl2和ZnCl2鹽水溶液中丙氨酸分子的締合與水合情況.基于水溶液中兩性離子形式的丙氨酸較穩(wěn)定,因而本文中所提丙氨酸分子均為兩性離子形式結構.結合量子化學計算與經(jīng)典分子動力學模擬結果,分別分析在稀、濃鹽水溶液中陽離子、陰離子與丙氨酸各極性基團之間的相互作用對丙氨酸分子締合的影響.

2 計算方法

采用密度泛函理論(DFT)研究丙氨酸與陽離子、陰離子間的相互作用和水合作用.考慮到丙氨酸兩性離子的結構形式存在電荷分離,丙氨酸與締合結構以及離子的水合結構的幾何優(yōu)化和能量的計算均采用長程校正的雜化交換相關泛函CAMB3LYP32與aug-cc-pVDZ(縮寫為aVDZ)基組.所有優(yōu)化的結構都在同一計算水平進行了振動頻率計算與分析,以確保所得結構為能量極小點的結構.相互作用能與離子水合能的計算都進行了基組重疊誤差校正.基組重疊誤差校正和自然鍵軌道(NBO)電荷布居分析也都采用與結構優(yōu)化相一致的方法和基組.所有DFT計算均采用Gaussian09軟件包.33采用經(jīng)典分子動力學模擬方法研究丙氨酸在鹽水溶液中的締合與水合行為,研究體系為加入NaCl、CuCl2和ZnCl2的三種丙氨酸鹽水溶液.模擬的周期性立方體盒子中包含20個兩性離子形式的丙氨酸分子,1600個水分子(丙氨酸的飽和溶解度為1600個水分子溶解～53個丙氨酸分子),稀、濃鹽水溶液中無機鹽的濃度分別為0.17和6.18mol·kg-1.稀鹽水濃度對應的是生理鹽水濃度.在稀溶液Zn2+與丙氨酸羧基之間無直接締合,為了研究陽離子與羧基間直接締合作用對丙氨酸分子締合的影響,考慮了較高濃度的濃鹽水溶液.所有經(jīng)典分子動力學模擬計算均采用DL_POLY_4.03軟件包34完成.模擬選擇CHARMM力場,SPCE水分子模型,非鍵截斷距離設定為1.20nm.采用Verletvelocity算法35對每一步的運動方程進行求解,積分步長設定為1fs.為確保溶液模擬的收斂,對一些樣本體系的模擬結果進行了勢能、溫度與壓力等參數(shù)波動變化的分析,并根據(jù)粒子擴散系數(shù)估算分子動力學模擬需要的平衡時間(相關體系需約0.5ns的平衡).所有模擬均先在等溫等壓(NPT)系綜中進行2ns的模擬,模擬溫度為298.15K,壓力為105Pa.通過Nosé-Hoover方法36-38控制溫度和壓力.然后在正則系綜(NVT)下進行6ns模擬,其中后4ns用于數(shù)據(jù)分析.溫度同樣也采用Nosé-Hoover方法控制.溫度和壓力控制的時間間隔分別為0.5和1ps.

圖1 丙氨酸分子的兩種締合結構Fig.1 Two associated structuresof alanines B:bidentate,D:double

表1 丙氨酸結構的能量參數(shù)aTable1 Energy parametersof alanine structuresa

3 結果與討論

3.1丙氨酸的結構參數(shù)

采用CAM-B3LYP/aVDZ方法進行幾何結構優(yōu)化,得到兩種丙氨酸分子的締合結構,進行頻率分析確保所得到結構為能量局域極小點的結構.密度泛函理論(DFT)得到兩種丙氨酸分子間的締合結構:兩個氨酸分子的氨基和羧基以單齒形式兩兩締合(用D表示)、一分子丙氨酸的氨基與另一分子丙氨酸的羧基以雙齒形式單締合(用B表示)(圖1).

由表1列出的結合能計算結果可知,Na+、Cu2+和Zn2+與丙氨酸羧基以雙齒形式締合的結構更穩(wěn)定,這一結論與喬龍光等39的結論相一致.Cl-則是與氨基的氫原子形成離子氫鍵(iHB).在氣相中,丙氨酸分子的氨基與Cl-之間締合為單齒形式(Ala-Cl--M,M代表單齒),這可能是由于氮原子與Cl-之間排斥作用的緣故.表1中數(shù)據(jù)表明,丙氨酸與陽離子、陰離子間的相互作用均強于丙氨酸分子間的締合作用,即氣相中丙氨酸分子更傾向于與其周圍的陽離子、陰離子發(fā)生相互作用,尤其是與二價金屬離子Cu2+、Zn2+間的相互作用更為顯著,這與Jensen和Gordon14的結論相一致.丙氨酸與Cl-間的作用相對弱一些,但是仍強于丙氨酸分子間的締合作用.綜上所述,離子或分子與丙氨酸相互作用的強弱順序是:X2+＞Na+＞Cl-≈Ala(X2+=Zn2+,Cu2+),這一結論與El-Dossoki40的實驗結果以及Bush等41計算的紅外光譜結果相一致.

水溶性的氨基酸在中性水溶液中的結構多為兩性離子形式,這意味著,氨基酸分子的電荷分離形式有利于氨基酸分子的水合作用與氨基酸的溶解.42NBO電荷布居分析結果顯示,丙氨酸分子間的締合使得其各基團上的電荷重新分布,氨基上的正電荷更多,羧基部位負電荷也相應增加,如表2所示.陽離子、陰離子與丙氨酸間的結合也均導致丙氨酸各基團的電荷發(fā)生相應變化.丙氨酸與Cu2+發(fā)生相互作用,使得其羧基上的負電荷量明顯減少,而氨基與甲基、亞甲基上的正電荷量則均明顯增加,金屬離子所帶的正電荷量也有所降低.亦即,二價金屬陽離子與羧基之間存在電荷轉移,尤其是Cu2+更為明顯.這一作用影響到了氨基、甲基與亞甲基等基團上的電荷分布,將會增強這些基團的水合作用.與金屬離子相比,Cl-對氨基的電荷分布影響較小,氨基上的正電荷量只有輕微的降低,而羧基上負電荷則明顯增加,這對羧基與金屬離子之間的相互作用起到一定的促進作用.

表2 丙氨酸和各結構中基團的電荷(e)aTable2 Charges(e)on somegroupsof alanine and some structuresa

圖2 CAM-B3LYP/aVDZ方法計算得到的陽離子與陰離子的水合能(ΔEb)Fig.2 Hydration energies(ΔEb)of cationsand anion obtained at CAM-B3LYP/aVDZ level

圖3 兩個丙氨酸分子中兩個α-C的徑向分布函數(shù)圖Fig.3 Radialdistribution functionsbetween twoα-C of two alaninemolecules

電荷布居分析結果顯示陽離子、陰離子與丙氨酸之間的相互作用均增強了丙氨酸兩性離子結構的電荷分離,這將促進丙氨酸各極性基團的水合作用.但是溶液中陽離子、陰離子與丙氨酸之間的相互作用情況,還需考慮離子水合作用的競爭,尤其是陽離子的水合作用.如圖2所示,CAM-B3LYP/ aVDZ方法的計算結果顯示,金屬陽離子的水合作用較強,二價陽離子尤為顯著.在稀溶液與濃溶液中,陽離子、陰離子水合作用的競爭對氨基酸的締合與水合作用影響應有所不同,因此需進一步結合經(jīng)典分子動力學方法研究鹽水溶液中丙氨酸分子的締合情況.

3.2水溶液中丙氨酸的締合結構

由圖3中α-C的徑向分布函數(shù)(RDF)結果可知,在0.69mol·kg-1丙氨酸水溶液中丙氨酸分子間存在相互締合的現(xiàn)象.RDF圖在0.39-0.64nm范圍內有兩個峰,且峰的寬度較大.結合CAM-B3LYP/aVDZ方法的計算結果可知,在水溶液中,丙氨酸分子間締合結構除了(Ala)2-B和(Ala)2-D兩種形式以外,還存在一種單齒形式的單締合結構(表示為:(Ala)2-M),如圖4所示.在氣相中,因沒考慮水合作用,陽離子與羧基(―COO-)的兩個氧原子之間可同時存在較強的相互作用,因而單齒締合結構很不穩(wěn)定.但是在鹽水溶液中由于存在極性基團的水合作用以及其與離子之間的相互作用,使得丙氨酸分子間的締合可以存在單齒形式,(Ala)2-M.由RDF中兩個α-C之間的位置關系可知,溶液中丙氨酸分子間締合結構以(Ala)2-D、(Ala)2-B形式為主,這與氣相中CAM-B3LYP/aVDZ方法計算的丙氨酸締合情況基本一致.在經(jīng)典分子動力學模擬結果中,(Ala)2-B和(Ala)2-D結構的兩個α-C之間的距離與CAMB3LYP/aVDZ方法的計算結果略有不同,這可能是水溶液中丙氨酸的水合作用導致這兩種的締合結構與氣相中的締合結構存在一定程度上的差異.

圖4 由經(jīng)典分子動力學模擬結果獲得的在水溶液中丙氨酸分子的三種典型締合結構Fig.4 Three typicalassociated structuresof alanine in aqueoussolution obtained from classicalmolecular dynam icssimulation results

表3 鹽水溶液中丙氨酸α-C原子之間徑向分布數(shù)據(jù)Table 3 Data of radialdistribution forα-C atom sof the alanine in saline solutiona

如圖3所示,丙氨酸水溶液中加入少量NaCl、CuCl2或ZnCl2時,丙氨酸分子締合程度均降低,但不同鹽水溶液中,丙氨酸分子的締合方式存在差異.結合圖4可知,溶液中加入少量NaCl、ZnCl2時,丙氨酸分子間的締合結構更傾向于(Ala)2-M形式.但少量CuCl2存在時,溶液中丙氨酸分子間的締合結構變化并不明顯,這可能是由于Cu2+與羧基間的相互作用較弱.與稀溶液中的情況不同,濃鹽水溶液中丙氨酸分子間的締合程度受到其周圍陽離子、陰離子的影響更顯著.三種濃鹽水溶液中,丙氨酸分子間締合結構的選擇性較強.在NaCl或CuCl2濃鹽水溶液中,丙氨酸分子的締合結構傾向于(Ala)2-M構型,在ZnCl2濃鹽水溶液時,丙氨酸分子的締合結構則是以(Ala)2-B形式為主.表3中列出了與圖3相對應的各體系丙氨酸α-C原子間的相互作用情況.

在稀鹽水溶液中,離子水合與丙氨酸水合之間的競爭對丙氨酸分子間締合的影響不大,但在濃鹽水溶液中則是不可忽略的因素.不同金屬陽離子對丙氨酸分子締合作用的影響程度各不同,這可能與陽離子的水合以及離子的締合強弱有關.以下部分,將具體分析陽離子、陰離子水合作用以及其與丙氨酸極性基團間的相互作用對丙氨酸分子締合的影響.

3.3 丙氨酸與陽離子間的相互作用

3.3.1 羧基與金屬陽離子([COO--Xn+])的相互作用

CAM-B3LYP/aVDZ計算結果顯示,氣相中丙氨酸分子與金屬離子均是以雙齒形式相互作用.但在鹽水溶液中,由于陽離子的水合作用以及不同陽離子的水合作用存在差異,導致其與丙氨酸分子間的相互作用結構有所不同.由圖5顯示,在稀鹽水溶液中,Na+、Cu2+與羧基氧原子發(fā)生直接相互作用,而Zn2+與羧基氧原子之間并不存在直接相互作用.羧基氧原子與Zn2+在0.4275nm處存在較明顯的峰.經(jīng)典分子動力學路徑文件的團簇分析結果顯示,溶液中存在羧基與Zn2+間的溶劑共享離子對結構(見圖6),這意味著Zn2+與處于其第二殼層的丙氨酸羧基之間仍然存在較明顯相互作用.Mohammed43和Bock44等在其量子力學/分子力學(QM/MM)的分子動力學模擬研究中也提到,Zn2+的第一水殼層比Cu2+緊湊,且在其20ps的模擬時間內,Zn2+的第一、二水殼層之間不存在水分子交換,其第二水殼層中氫鍵較規(guī)則,這都間接地證明了Zn2+對其第二水殼層的束縛仍較明顯.在稀ZnCl2水溶液中,Zn2+與處于其第二殼層的丙氨酸羧基之間的這種間接相互作用仍導致丙氨酸分子締合程度的降低(見圖3).但是在濃ZnCl2水溶液中,丙氨酸分子則是直接與Zn2+發(fā)生相互作用,如圖5所示,這可能是水溶液中丙氨酸與陽離子的水合競爭使Zn2+水合不完全導致的.團簇分析結果顯示,與丙氨酸水溶液相比,在稀ZnCl2鹽水溶液中,丙氨酸分子締合結構的數(shù)目有所減少,游離丙氨酸分子數(shù)目則明顯增多,如圖7所示.

圖5 丙氨酸羧基氧原子與金屬離子相互作用的徑向分布圖Fig.5 Radialdistribution curves for interaction between oxygen atomsof alanine carboxylandmetal ions

圖6 ZnCl2稀溶液中丙氨酸與Zn2+的相互作用結構(溶劑共享締合形式)Fig.6 Interaction structure(solventshared associated form)between alanine and Zn2+in ZnCl2dilute solution

圖7 丙氨酸純水溶液與丙氨酸-ZnCl2稀溶液中各團簇的數(shù)目Fig.7 Number of different typesof clusters in alanine purewater solution and its ZnCl2dilute solution

由圖8可知,稀鹽水溶液中,羧基碳原子與金屬離子間的[C(COO-)-Na+]徑向分布函數(shù)在0.2575和0.3275nm處出現(xiàn)兩個峰,說明Na+與羧基存在單齒、雙齒兩種相互作用形式.由于第一峰更明顯,可知雙齒作用形式更穩(wěn)定,兩種結構如圖9所示.Cu2+與羧基α-C原子間的徑向分布函數(shù)在0.2575nm處的峰較明顯,而與單齒結構相對應的峰則不明顯,由此可知主要以雙齒結構為主.濃鹽水溶液中,三種金屬離子與羧基氧原子之間的相互作用形式則基本趨于一致,兩種相互作用結構都存在,但均以雙齒結構為主.值得一提的是在丙氨酸的ZnCl2濃溶液中,Zn2+與羧基碳原子相互作用的兩個徑向分布峰較窄且區(qū)分較明顯,說明兩者之間形成的單、雙齒結構存在明顯區(qū)別.Zn2+與羧基碳原子的徑向分布圖相互作用在0.4500nm處還存在一個較明顯的峰,對應的結構為羧基與Zn2+間的溶劑共享締合結構(見圖6).如圖8所示,在稀ZnCl2與濃ZnCl2溶液中丙氨酸羧基與Zn2+間均存在這種溶劑共享締合結構.

圖8 羧基碳原子與金屬離子間相互作用的徑向分布圖Fig.8 Radialdistribution curves for interaction between the carbon atom of carboxyland metal ions

圖9 丙氨酸與Na+的兩種締合結構(單,雙齒結構)

Fig.9 Two associated structuresof alanine and Na+(monodentateand bidentate)

3.3.2 金屬離子水合作用的影響

由圖10可知,Na+、Cu2+和Zn2+三種金屬離子水合性質存在一定差異.Zn2+的第一、二水殼層徑向分布峰間隔清晰,存在較明顯的第二殼層,說明Zn2+對第二殼層水分子的束縛作用仍較強,這也與Mohammed43和Bock44等的結論相一致.Mohammed等43以及Schwenk與Rode45的研究結果顯示,Zn2+對第一水殼層束縛比Cu2+強.與Zn2+相同,Na+的水合殼層包含6個水分子,但其水合作用較Zn2+弱.44由于Jahn-Teller效應,Cu2+的水合殼層可能包含5個或6個水分子.46另外,Zn2+與羧基氧原子的徑向分布峰位于0.4275nm處(見圖5),而Zn2+與水分子中氧原子的徑向分布峰也是在同一位置處出現(xiàn)較強的峰(見圖10).這也進一步說明,當丙氨酸出現(xiàn)在Zn2+的第二殼層時,Zn2+與羧基之間仍存在較強的相互作用.與二價離子相比,Na+水合作用較弱,因而在稀溶液和濃溶液中Na+與羧基間的相互作用均較強.Na+與羧基之間這種較強的相互作用,一方面降低丙氨酸分子間締合結構的締合程度,另一方面也可使得丙氨酸分子間締合結構傾向于締合程度較低(Ala)2-M結構.

在水溶液中,離子水合和丙氨酸水合之間是相互競爭關系,金屬陽離子與丙氨酸間的相互作用會受到金屬陽離子水合或其他締合作用的制約.稀鹽水溶液中離子的水合與丙氨酸分子水合的競爭對丙氨酸分子締合的影響較小,因而影響丙氨酸分子締合的主要因素為陽離子、陰離子與丙氨酸極性基團間的相互作用.溶液中Cu2+的水合能遠高于Na+(見圖2),這不利于Cu2+與丙氨酸分子間的相互作用,因而稀溶液中Cu2+對丙氨酸分子間相互作用的影響比Na+小.另外,Xia等46和Marcus47的研究結果也表明Cu2+與Cl-間存在較強的締合作用,這也會阻礙Cu2+與丙氨酸的相互作用.與加入CuCl2相比,丙氨酸水溶液中加入NaCl后,Na+與羧基間的相互作用較明顯,丙氨酸分子締合作用也降低更顯著,如圖5和圖8中RDF圖所示.

圖10 鹽水溶液中金屬離子水合作用的徑向分布圖Fig.10 Radialdistribution curves for interaction between m etal ionsand waterm olecule in saline solution OW:oxygen atom inwatermolecules

圖11 氨基氮原子與Cl-相互作用的徑向分布圖Fig.11 Radialdistribution curves for interaction between nitrogen atom of am ino group and Cl-

濃鹽水溶液中必須考慮金屬離子水合作用對丙氨酸分子間締合以及丙氨酸與離子相互作用的影響.一方面,大量陽離子、陰離子的水合作用與丙氨酸極性基團的水合作用之間存在明顯的競爭,這將有利于丙氨酸分子的締合;另一方面,當陽離子與陰離子的濃度較高時,其與丙氨酸極性基團之間的相互作用也隨之增強,這又不利于丙氨酸分子的締合,兩方面作用相結合共同影響水溶液中丙氨酸分子締合.在濃鹽水溶液中,因為Cu2+水合受限而水合不完全,導致Cu2+與羧基間的相互作用對丙氨酸分子締合程度的影響不再弱于Na+.因而與NaCl濃溶液相比,CuCl2濃溶液中丙氨酸分子締合減弱程度更高,有別于稀鹽溶液中的情形.同樣地,在ZnCl2濃溶液中,Zn2+的水合作用也受限,Zn2+可與羧基直接相互作用(見圖8和圖10),這也就導致ZnCl2濃溶液中丙氨酸分子的締合程度大為降低,如圖3所示.

3.4 氨基與陰離子(Cl-)的相互作用

水溶液中Cl的水合作用較弱,21,46因此Cl-與氨基的相互作用不可忽略.圖11中的氨基氮原子與Cl-間的徑向分布函數(shù),在0.3125nm處存在明顯的峰,可知Cl-與氨基之間存在明顯的相互作用.圖3和表3的數(shù)據(jù)表明,與NaCl稀溶液相比,ZnCl2稀溶液中丙氨酸分子的締合程度減弱更為明顯.但是稀溶液中Na+與羧基之間存在直接相互作用,而Zn2+與羧基之間只存在間接相互作用(見圖5與圖6).這可能是因為ZnCl2溶液中Cl-的濃度是NaCl溶液的兩倍,且Cl-與氨基間的相互作用對丙氨酸分子的締合也起到了一定的阻礙作用.Cl-與氨基的相互作用影響丙氨酸分子締合的觀點與Tomé等48得出的相關結論存在一定分歧,但是其研究側重于討論Cl-與鏈端疏水甲基的相互作用,而沒有分析Cl-與氨基的相互作用.丙氨酸分子的疏水鏈較短,鏈端甲基距離兩個極性基團較近,因而影響丙氨酸分子締合程度的主要因素為具有較強親水性的羧基或氨基與離子間的作用.另外,與NaCl稀溶液相比,CuCl2稀溶液中Cl-濃度是NaCl稀溶液中的兩倍,但是在CuCl2溶液中丙氨酸分子間的締合程度減弱仍較小.由此可推知,相對于金屬陽離子與羧基的相互作用而言,Cl-與氨基間的相互作用對丙氨酸分子締合的影響較小.

4 結論

結合DFT方法與經(jīng)典分子動力學模擬,研究了水溶液中丙氨酸分子間的締合以及一些陽離子、陰離子的存在對其締合作用的影響.NBO電荷布居分析結果顯示陽離子、陰離子與丙氨酸之間的相互作用均增強了丙氨酸兩性離子結構的電荷分離,非極性基團上的電荷也發(fā)生相應的變化,這將增強丙氨酸分子的水合作用.DFT與經(jīng)典分子動力學的研究結果表明,丙氨酸在水溶液中存在單齒、雙齒形式的締合結構.鹽水溶液中的陽離子、陰離子均會影響丙氨酸分子的締合程度與締合結構形式.稀鹽水溶液中,影響丙氨酸分子締合的主要因素是陽離子、陰離子與丙氨酸極性基團之間相互作用.在丙氨酸溶液中離子濃度較低時,Cu2+與Na+與羧基均存在直接相互作用.但是Cu2+水合作用強而導致其與羧基的相互作用弱,因而與比Na+相比,Cu2+對丙氨酸分子相互作用程度影響較小.在稀ZnCl2水溶液中,Zn2+與羧基并無直接相互作用,但是Zn2+與處于其第二殼層的丙氨酸羧基之間仍有較強的相互作用,這一種間接相互作用仍能阻礙丙氨酸分子間的締合.另外,Cl-與氨基的相互作用對丙氨酸分子間締合的影響相對較小但也不可忽略.在濃鹽水溶液中,大量陽離子的水合與丙氨酸極性基團的水合之間存在明顯競爭,這將有利于丙氨酸分子的締合.但是大量陽離子、陰離子與丙氨酸極性基團相互作用,又將阻礙丙氨酸分子間的締合.因此,盡管在濃鹽水溶液中離子水合與丙氨酸水合之間的競爭不可忽略,但陽離子、陰離子與丙氨酸極性基團之間的相互作用對丙氨酸分子締合的影響可能仍起主導作用.

(1)Heaton,A.L.;Bowman,V.N.;Oomens,J.;Steill,J.D.; Armentrout,P.B.J.Phys.Chem.A 2009,113(19),5519.doi: 10.1021/jp9008064

(2)Armentrout,P.B.;Armentrout,E.I.;Clark,A.A.;Cooper,T.E.; Stennett,E.M.S.;Carl,D.R.J.Phys.Chem.B 2010,114(11), 3927.doi:10.1021/jp911219u

(3)Dunbar,R.C.;Steill,J.D.;Polfer,N.C.;Oomens,J.J.Phys. Chem.B 2009,113(31),10552.doi:10.1021/jp905060n

(4)O'Brien,J.T.;Prell,J.S.;Steill,J.D.;Oomens,J.;Williams,E. R.J.Phys.Chem.A 2008,112(43),10823.doi:10.1021/ jp805787e

(5)Prell,J.S.;Demireva,M.;Oomens,J.;Williams,E.R.J.Am. Chem.Soc.2009,131(3),1232.doi:10.1021/ja808177z

(6)Frossard,E.;Bucher,M.;M?chler,F.;Mozafar,A.;Hurrell,R. J.Sci.Food Agric.2000,80,861.

(7)Wang,C.Y.;Guo,J.S.;Tian,J.;Cui,W.L.J.Jilin Medical College 2009,30(2),99. [王春艷,郭景森,田 晶,崔萬麗.吉林醫(yī)藥學院學報,2009,30(2),99.]

(8)Wingenfeld,K.;Hellhammer,D.H.;Schmidt,I.;Wagner,D.; Meinlschm idt,G.;Heim,C.J.Psychosom.Obstet.Gynaecol. 2009,30(4),282.doi:10.3109/01674820903254732

(9)Komoroski,R.A.;Pearce,J.M.Magn.Reson.Med.2008,60, 21.doi:10.1002/m rm.v60:1

(10)Dong,X.Y.;Du,W.J.;Liu,F.F.Acta Phys.-Chim.Sin.2012, 28(11),2735. [董曉燕,都文婕,劉夫鋒.物理化學學報, 2012,28(11),2735.]doi:10.3866/PKU.WHXB201207162

(11)Zhao,Y.P.;Ai,H.Q.;Chen,J.P.;Yang,A.B.;Qi,Z.N.Acta Phys.-Chim.Sin.2010,26(12),3322.[趙永平,艾洪奇,陳金鵬,楊愛彬,齊中囡.物理化學學報,2010,26(12),3322.]doi: 10.3866/PKU.WHXB20101215

(12)Zhu,Y.C.;Wang,E.Q.;Ma,G.L.;Kang,Y.B.;Zhao,L.H.; Liu,Y.Z.Acta Phys.-Chim.Sin.2014,30(1),1.[朱云城,王二瓊,馬國林,康彥彪,趙林泓,劉揚中.物理化學學報,2014, 30(1),1.]doi:10.3866/PKU.WHXB201311263

(13)Eyal,A.M.;Bressler,E.Biotechnol.Bioeng.1993,41(3),287.

(14)Jensen,J.H.;Gordon,M.S.J.Am.Chem.Soc.1995,117(31), 8159.doi:10.1021/ja00136a013

(15)Hu,C.H.;Shen,M.Z.;Schaefer,H.F.,III.J.Am.Chem.Soc. 1993,115(7),2923.doi:10.1021/ja00060a046

(16)Qin,P.H.;Lü,W.C.;Qin,W.;Zhang,W.;Xie,H.Chem.Res. Chin.Univ.2014,30(1),125.doi:10.1007/s40242-014-3303-z

(17)Gordon,M.S.;Jensen,J.H.AccountsChem.Res.1996,29(11), 536.doi:10.1021/ar9600594

(18)Qiu,X.M.;Lei,Q.F.;Fang,W.J.;Lin,R.S.Acta Chim.Sin. 2009,67(7),607.[邱曉梅,雷群芳,方文軍,林瑞森.化學學報,2009,67(7),607.]

(19)Wyttenbach,T.;Bushnell,J.E.;Bowers,M.T.J.Am.Chem. Soc.1998,120(20),5098.doi:10.1021/ja9801238

(20)Kushwaha,P.S.;M ishra,P.C.J.Mol.Struct.-Theochem 2001, 549,229.doi:10.1016/S0166-1280(01)00423-7

(21)Tomé,L.I.N.;Pinho,S.P.;Jorge,M.;Gomes,J.R.B.; Coutinho,J.A.P.J.Phys.Chem.B 2013,117,6116.

(22)Carta,R.;Tola,G.J.Chem.Eng.Data 1996,41(3),414.doi: 10.1021/je9501853

(23)Marino,T.;Russo,N.;Toscano,M.J.Inorg.Biochem.2000,79, 179.doi:10.1016/S0162-0134(99)00242-1

(24)Pulkkinen,S.;Noguera,M.;Rodríguez-Santiago,L.;Sodupe, M.;Bertran,J.Chem.-Eur.J.2000,6(23),4393.doi:10.1002/ 1521-3765(20001201)6:23＜4393::AID-CHEM 4393＞3.0.CO;2-H

(25)Apse,M.P.;Aharon,G.S.;Snedden,W.A.;Blumwald,E.Science 1999,285,1256.doi:10.1126/science.285.5431.1256

(26)Xu,J.H.;Hu,C.W.Acta Chim.Sin.2006,64(16),1622.[徐建華,胡常偉.化學學報,2006,64(16),1622.]

(27)Khoshkbarchi,M.K.;Vera,J.H.Ind.Eng.Chem.Res.1997,36 (6),2445.doi:10.1021/ie9606395

(28)Xia,F.F.;Yi,H.B.;Zeng,D.W.J.Phys.Chem.A 2009,113 (51),14029.doi:10.1021/jp909092p

(29)Bryantsev,V.S.;Diallo,M.S.;Goddard,W.A.,III.J.Phys. Chem.B 2008,112(32),9709.doi:10.1021/jp802665d

(30)Harris,D.J.;Brodholt,J.P.;Harding,J.H.;Sherman,D.M. Mol.Phys.2001,99(10),825.doi:10.1080/ 00268970010015588

(31)Pye,C.C.;Corbeil,C.R.;Rudolph,W.W.Phys.Chem.Chem. Phys.2006,8,5428.doi:10.1039/b610084h

(32)Yanai,T.;Tew,D.P.;Handy,N.C.Chem.Phys.Lett.2004,393, 51.doi:10.1016/j.cplett.2004.06.011

(33)Frisch,M.J.;Trucks,H.B.;Schlegel,G.E.;etal.Gaussian 09, Revision A.1;Gaussian Inc.:Wallingford,CT,2009.

(34)Todorov,I.T.;Smith,W.;Trachenko,K.;Dove,M.T.J.Mater. Chem.2006,16,1911.doi:10.1039/b517931a

(35)Adcock,S.A.;M cCammon,J.A.Chem.Rev.2006,106(5), 1589.doi:10.1021/cr040426m

(36)Nosé,S.Mol.Phys.1984,52,255.doi:10.1080/ 00268978400101201

(37)Hoover,W.G.Phys.Rev.A 1985,31,1695.doi:10.1103/ PhysRevA.31.1695

(38)Martyna,G.J.;K lein,M.L.;Tuckerman,M.E.J.Chem.Phys. 1992,97(4),2635.doi:10.1063/1.463940

(39)Qiao,L.G.;Fan,J.F.;Yang,C.H.Acta Chim.Sin.2007,65 (17),1751.[喬龍光,樊建芬,楊春紅.化學學報,2007,65 (17),1751.]

(40)El-Dossoki,F.I.J.Solution Chem.2010,39,1311.doi:10.1007/ s10953-010-9580-3

(41)Bush,M.F.;Oomens,J.;Saykally,R.J.;W illiams,E.R.J.Am. Chem.Soc.2008,130(20),6463.doi:10.1021/ja711343q

(42)Meng,X.J.Acta Phys.-Chim.Sin.2006,22(1),98.[孟祥軍.物理化學學報,2006,22(1),98.]doi:10.3866/PKU. WHXB20060120

(43)Mohammed,A.M.;Loeffler,H.H.;Inada,Y.;Tanada,K.; Funahashi,S.J.Mol.Liq.2005,119,55.doi:10.1016/j. molliq.2004.10.008

(44)Bock,C.W.;Markham,G.D.;Katz,A.K.;Glusker,J.P.Theor. Chem.Acc.2006,115,100.doi:10.1007/s00214-005-0056-2

(45)Schwenk,C.F.;Rode,B.M.J.Chem.Phys.2003,119, 9523.doi:10.1063/1.1614224

(46)Xia,F.F.;Yi,H.B.;Zeng,D.W.J.Phys.Chem.A 2010,114 (32),8406.doi:10.1021/jp1000804

(47)Marcus,Y.J.Chem.Soc.Faraday Trans.1991,87(18), 2995.doi:10.1039/ft9918702995

(48)Tomé,L.I.N.;Jorge,M.;Gomes,J.R.B.;Coutinho,J.A.P. J.Phys.Chem.B 2010,114(49),16450.doi:10.1021/jp104626w

Effec ts o f In terac tions betw een Ions and A lanine Po lar Groups on A lanine Associations in Saline So lu tion:Density Func tional Theo ry and Mo lecu lar Dynam ics Sim u lation

WANG Ying YIHai-Bo*LIHui-Ji DAIQian CAO Zhi-Wei LU Yang
(State Key Laboratory ofChemo/Biosensing and Chemometrics,College ofChemistry and Chemical Engineering, Hunan University,Changsha 410082,P.R.China)

Density functional theory(DFT)and classicalmolecular dynam ics simulationswere used to study the effects of the interactions between zw itterionic alanine and some ions(Na+,Cu2+,Zn2+,and Cl-)in saline solution on the association of alanine m olecules.The DFT calculation results show that the association of alanine w ith these ions can enhance charge separation ofzw itterionic alanine.Classicalmolecular dynam ics simulation results also show that three associated structures ofzw itterionic alaninemolecules are present in alanine aqueous solution,and the associations can be weakened to a certain extentby the interactions between the cations/anions and alanine polargroups.The interaction between a cation and the carboxylgroup ofalanine can be greatly affected by hydration of the cation in dilute saline solution.The interaction between Cu2+and alanine ismuch stronger than thatbetween Na+and alanine in the gas phase,but the situation is reversed in dilute aqueous solution,because the hydration ofCu2+ismuch stronger than thatofNa+.In dilute ZnCl2aqueous solution,the interaction between Zn2+and the carboxylgroup of the alaninemolecule is less direct,because of the first hydration shell of Zn2+.However,indirect interactions between Zn2+and alanine still lead to adecreased association among alaninemolecules.In addition,the interactions of cations/anions w ith alanine notonly weaken the association between alanine molecules,butalso result in transformation between two typical conformations of associated alanine m olecules.The ion concentration affects the conformations of associated cation/anion-alanine species,and associated alaninem olecules.

Zw itterion;Cation;Hydration;Association;Intermolecular interaction

O641

icle]

10.3866/PKU.WHXB201504031www.whxb.pku.edu.cn

Received:January 22,2015;Revised:April2,2015;Published onWeb:April3,2015.

?Corresponding author.Email:hbyi@hnu.edu.cn;Tel:+86-13574869952.

The projectwas supported by the National Natural Science Foundation of China(21073056,J1210040)and Students Innovation Training Project in Hunan Province,China(2012).

國家自然科學基金(21073056,J1210040)與湖南省“大學生創(chuàng)新訓練計劃”項目(2012)資助

?Editorialofficeof Acta Physico-Chim ica Sinica