Nguyen Phuong H, Gorbunov Roman D, Stock Gerhard
Institute of Physical and Theoretical Chemistry, J. W. Goethe University, Frankfurt, Germany.
Biophys J. 2006 Aug 15;91(4):1224-34. doi: 10.1529/biophysj.106.084996. Epub 2006 May 26.
Employing nonequilibrium molecular dynamics simulations, a comprehensive computational study of the photoinduced conformational dynamics of a photoswitchable bicyclic azobenzene octapeptide is presented. The calculation of time-dependent probability distributions along various global and local reaction coordinates reveals that the conformational rearrangement of the peptide is rather complex and occurs on at least four timescales: 1) After photoexcitation, the azobenzene unit of the molecule undergoes nonadiabatic photoisomerization within 0.2 ps. 2) On the picosecond timescale, the cooling (13 ps) and the stretching (14 ps) of the photoexcited peptide is observed. 3) Most reaction coordinates exhibit a 50-100 ps component reflecting a fast conformational rearrangement. 4) The 500-1000 ps component observed in the simulation accounts for the slow diffusion-controlled conformational equilibration of the system. The simulation of the photoinduced molecular processes is in remarkable agreement with time-resolved optical and infrared experiments, although the calculated cooling as well as the initial conformational rearrangements of the peptide appear to be somewhat too slow. Based on an ab initio parameterized vibrational Hamiltonian, the time-dependent amide I frequency shift is calculated. Both intramolecular and solvent-induced contributions to the frequency shift were found to change by < or = 2 cm(-1), in reasonable agreement with experiment. The potential of transient infrared spectra to characterize the conformational dynamics of peptides is discussed in some detail.
利用非平衡分子动力学模拟,对一种可光开关的双环偶氮苯八肽的光诱导构象动力学进行了全面的计算研究。沿各种全局和局部反应坐标计算时间相关概率分布表明,该肽的构象重排相当复杂,且至少在四个时间尺度上发生:1)光激发后,分子中的偶氮苯单元在0.2皮秒内经历非绝热光异构化。2)在皮秒时间尺度上,观察到光激发肽的冷却(13皮秒)和拉伸(14皮秒)。3)大多数反应坐标显示出一个50 - 100皮秒的分量,反映了快速的构象重排。4)模拟中观察到的500 - 1000皮秒分量解释了系统缓慢的扩散控制构象平衡。光诱导分子过程的模拟与时间分辨光学和红外实验结果显著一致,尽管计算得到的肽的冷却以及初始构象重排似乎有点过慢。基于从头算参数化的振动哈密顿量,计算了随时间变化的酰胺I频率位移。发现分子内和溶剂诱导对频率位移的贡献变化均≤2厘米⁻¹,与实验结果合理相符。详细讨论了瞬态红外光谱表征肽构象动力学的潜力。