Tao Peng, Parquette Jon R, Hadad Christopher M
Department of Chemistry, The Ohio State University , 100 West 18th Avenue, Columbus, Ohio 43210, United States.
J Chem Theory Comput. 2012 Dec 11;8(12):5137-49. doi: 10.1021/ct2009335. Epub 2012 Sep 18.
Some unnatural polymers/oligomers have been designed to adopt a well-defined, compact, three-dimensional folding capability. Azobenzene units are common linkages in these oligomer designs. Two alternating pyridinedicarboxamide/m-(phenylazo)azobenzene oligomers that can fold into both right- and left-handed helices were studied computationally in order to understand their dynamical properties. Helical structures were shown to be the global minima among the many different conformations generated from the Monte Carlo simulations, and extended conformations have higher potential energies than compact ones. To understand the interconversion process between right- and left-handed helices, replica-exchange molecular dynamic (REMD) simulations were performed on both oligomers, and with this method, both right- and left-handed helices were successfully sampled during the simulations. REMD trajectories revealed twisted conformations as intermediate structures in the interconversion pathway between the two helical forms of these azobenzene oligomers. This mechanism was observed in both oligomers in current study and occurred locally in the larger oligomer. This discovery indicates that the interconversion between helical structures with different handedness goes through a compact and partially folded structure instead of globally unfold and extended structure. This is also verified by the nudged elastic band (NEB) calculations. The temperature weighted histogram analysis method (T-WHAM) was applied on the REMD results to generate contour maps of the potential of mean force (PMF). Analysis showed that right- and left-handed helices are equally sampled in these REMD simulations. In large oligomers, both right- and left-handed helices can be adopted by different parts of the molecule simultaneously. The interconversion between two helical forms can occur in the middle of the helical structure and not necessarily at the termini of the oligomer.
一些非天然聚合物/低聚物已被设计成具有明确、紧凑的三维折叠能力。偶氮苯单元是这些低聚物设计中常见的连接基团。为了了解两种交替的吡啶二甲酰胺/m-(苯基偶氮)偶氮苯低聚物的动力学性质,对其进行了计算研究,这两种低聚物都能折叠成右手螺旋和左手螺旋。在蒙特卡罗模拟产生的许多不同构象中,螺旋结构被证明是全局最小值,伸展构象的势能比紧凑构象更高。为了理解右手螺旋和左手螺旋之间的相互转化过程,对这两种低聚物都进行了副本交换分子动力学(REMD)模拟,通过这种方法,在模拟过程中成功地采样到了右手螺旋和左手螺旋。REMD轨迹显示,扭曲构象是这些偶氮苯低聚物两种螺旋形式相互转化途径中的中间结构。在本研究的两种低聚物中都观察到了这种机制,并且在较大的低聚物中局部发生。这一发现表明,不同手性螺旋结构之间的相互转化是通过紧凑的部分折叠结构,而不是全局展开和伸展的结构。这也通过推挤弹性带(NEB)计算得到了验证。将温度加权直方图分析方法(T-WHAM)应用于REMD结果,以生成平均力势能(PMF)的等高线图。分析表明,在这些REMD模拟中,右手螺旋和左手螺旋的采样是相等的。在大的低聚物中,分子的不同部分可以同时采用右手螺旋和左手螺旋。两种螺旋形式之间的相互转化可以发生在螺旋结构的中间,而不一定在低聚物的末端。