高等学校化学学报
高等學校化學學報
고등학교화학학보
Chemical Journal of Chinese Universities
2015年
11期
2211-2219
,共9页
汪双双%刘鹏%蔡文生%邵学广
汪雙雙%劉鵬%蔡文生%邵學廣
왕쌍쌍%류붕%채문생%소학엄
轮烷%穿梭%疏水性%自由能计算%多副本自适应偏置力(MW-ABF)方法
輪烷%穿梭%疏水性%自由能計算%多副本自適應偏置力(MW-ABF)方法
륜완%천사%소수성%자유능계산%다부본자괄응편치력(MW-ABF)방법
Rotaxane%Shuttling%Hydrophobicity%Free-energy calculation%Multiple walker adaptive biasing force(MW-ABF)
研究了3种由α-环糊精(α-CD)、不同线性结构及中间连接基团组成的轮烷的穿梭行为,其中,线性结构为烷基链( ALK)或聚乙二醇链( PEG),中间连接基团为联吡啶或联苯。采用分子动力学模拟方法对室温下这些轮烷在水或二甲基亚砜( DMSO)溶液中α-CD的穿梭过程进行了研究;利用自适应偏置力( ABF)和多副本ABF( MW-ABF)方法计算了穿梭过程的自由能变化。结果表明,由PEG链组成的轮烷在2种溶剂中穿梭所需克服的能垒均比由ALK链组成的轮烷低,并且前者在DMSO中穿梭所需克服的能垒比在水中的能垒要低。由联苯基团组成的轮烷比由联吡啶基团组成的轮烷的能垒要低很多。通过对自由能曲线进行分解表明,破坏联吡啶基团的溶剂化层是穿梭过程能垒形成的主要原因,因此采用疏水的联苯基团可大大降低该能垒高度;另一方面,由于PEG链比ALK链具有更强的亲水性,也导致穿梭过程能垒降低。此外,通过对2种计算方法进行比较发现, MW-ABF能够明显改善对构象空间的采样效率,加快收敛速度。
研究瞭3種由α-環糊精(α-CD)、不同線性結構及中間連接基糰組成的輪烷的穿梭行為,其中,線性結構為烷基鏈( ALK)或聚乙二醇鏈( PEG),中間連接基糰為聯吡啶或聯苯。採用分子動力學模擬方法對室溫下這些輪烷在水或二甲基亞砜( DMSO)溶液中α-CD的穿梭過程進行瞭研究;利用自適應偏置力( ABF)和多副本ABF( MW-ABF)方法計算瞭穿梭過程的自由能變化。結果錶明,由PEG鏈組成的輪烷在2種溶劑中穿梭所需剋服的能壘均比由ALK鏈組成的輪烷低,併且前者在DMSO中穿梭所需剋服的能壘比在水中的能壘要低。由聯苯基糰組成的輪烷比由聯吡啶基糰組成的輪烷的能壘要低很多。通過對自由能麯線進行分解錶明,破壞聯吡啶基糰的溶劑化層是穿梭過程能壘形成的主要原因,因此採用疏水的聯苯基糰可大大降低該能壘高度;另一方麵,由于PEG鏈比ALK鏈具有更彊的親水性,也導緻穿梭過程能壘降低。此外,通過對2種計算方法進行比較髮現, MW-ABF能夠明顯改善對構象空間的採樣效率,加快收斂速度。
연구료3충유α-배호정(α-CD)、불동선성결구급중간련접기단조성적륜완적천사행위,기중,선성결구위완기련( ALK)혹취을이순련( PEG),중간련접기단위련필정혹련분。채용분자동역학모의방법대실온하저사륜완재수혹이갑기아풍( DMSO)용액중α-CD적천사과정진행료연구;이용자괄응편치력( ABF)화다부본ABF( MW-ABF)방법계산료천사과정적자유능변화。결과표명,유PEG련조성적륜완재2충용제중천사소수극복적능루균비유ALK련조성적륜완저,병차전자재DMSO중천사소수극복적능루비재수중적능루요저。유련분기단조성적륜완비유련필정기단조성적륜완적능루요저흔다。통과대자유능곡선진행분해표명,파배련필정기단적용제화층시천사과정능루형성적주요원인,인차채용소수적련분기단가대대강저해능루고도;령일방면,유우PEG련비ALK련구유경강적친수성,야도치천사과정능루강저。차외,통과대2충계산방법진행비교발현, MW-ABF능구명현개선대구상공간적채양효솔,가쾌수렴속도。
In order to investigate the hydrophobic effect of chain-like structures on shuttling in further, three similar rotaxanes are studied in this work. The rotaxanes are formed by anα-CD, two dodecamethylene chains ( ALK) or two poly ( ethylene glycol ) ( PEG ) for stations, one bipyridinium moiety ( PY ) or one biphenyl moiety( PH) for linkers and large end groups at both sides. The shuttling of the rotaxanes was studied by means of molecular dynamics simulations( MD) combined with free-energy calculations in water and DMSO at room temperature. Two methods, the adaptive biasing force ( ABF ) method and the multiple walker ABF ( MW-ABF) method, a modified version of ABF, are adopted for calculating the free-energy change characte-rizing the shuttling process. The potentials of mean force( PMFs) for the three rotaxanes are determined. The free-energy barriers of the PMFs for the PEG-rotaxane are lower than those for the the ALK-rotaxane both in water and DMSO. Furthermore, the barriers for the PEG-rotaxanes in DMSO are lower than that in water, which is in accordance with the ALK-rotaxanes. The barriers for the PH-rotaxanes are significantly lower than those for the PY-rotaxanes. Partitioning the PMFs into free-energy components suggests that change of the the charged group by an hydrophobic biphenyl moiety or decrease of the hydrophobicity of the chain-like structure in two stations from ALK to PEG reduces the free-energy barrier with respect to the stable states in the stations. In addition, comparison of the two free-energy calculation methods shows that the MW-ABF method can significantly improve the uniformity of sampling and hence increase the computational efficiency.