研究论文

β发卡多肽Trpzip4折叠的副本交换分子动力学模拟

  • 廖晨伊 ,
  • 周健
展开
  • 华南理工大学化学与化工学院 广东省绿色化学产品技术重点实验室 广州 510640

收稿日期: 2013-01-04

  网络出版日期: 2013-03-01

基金资助

项目受教育部新世纪优秀人才支持计划(No. NCET-07-0313)、国家自然科学基金(No. 20876052)及广东省自然科学基金(No. S2011010002078)资助.

Replica Exchange Molecular Dynamics Simulations on the Folding of Trpzip4 β-Hairpin

  • Liao Chenyi ,
  • Zhou Jian
Expand
  • School of Chemistry and Chemical Engineering, Guangdong Provincial Key Lab for Green Chemical Product Technology, South China University of Technology, Guangzhou 510640

Received date: 2013-01-04

  Online published: 2013-03-01

Supported by

Project supported by the Program for New Century Excellent Talents in Universities, Ministry of Education, China (No. NCET-07-0313), the National Natural Science Foundation of China (No. 20876052) and the Natural Science Foundation of Guangdong Province(No. S2011010002078).

摘要

采用副本交换分子动力学对β发卡Trpzip4重折叠进行研究, 结果表明, 构象空间存在较大能垒时, 副本交换分子动力学(REMD)表现出比分子动力学(MD)更优的抽样效率. 288 K下, Trpzip4势能、骨架均方根偏差、溶剂可及表面积(SASA)在REMD中呈逐渐降低趋势. 采样得到两种特定形态的低势能构象: β发卡和螺旋-卷曲. β发卡中, 第3组氢键Asp46容易与Thr49形成氢键; 转角Thr49的羟基(OH)倾向于与Asp46的羧基(COO-)或主链上的C=O形成氢键, 强化了Asp46与Thr49的相互作用, 导致转角形成折回弯度; 与Thr49相比, Thr51的羟基倾向与水分子作用, 甲基朝内, 因此与Asp46的侧链作用微小. 螺旋-卷曲中, 吲哚环-甲基为主要疏水形式. Trpzip4在折叠成β发卡过程中, 转角氢键影响整个β发卡构象的形成. 转角氢键具有距离优势和强侧链相互作用, 最先形成.

本文引用格式

廖晨伊 , 周健 . β发卡多肽Trpzip4折叠的副本交换分子动力学模拟[J]. 化学学报, 2013 , 71(04) : 593 -601 . DOI: 10.6023/A13010015

Abstract

β-Hairpin is an essential secondary structure unit of protein. Understanding the formation mechanism of the β hairpin and kinetic stability helps to gain insight into the formation of protein secondary structure as well as the mechanism of protein folding. Replica exchange molecular dynamics (REMD) approach is applied to investigate the folding mechanism of Trpzip4 β-hairpin. The REMD method is more efficient than conventional MD in searching for a representative set of the low energy minima in the existence of a large energy gap between the native state and any of the other possible state. The potential energy, backbone RMSD and solvent-accessible surface area (SASA) present descending trends in the process of REMD at 288 K. Methyl groups (most in threonines) show hydrophobic interactions with SASA decrements. The hydrophobic interactions and intra hydrogen bond forming drive the folding and maintain the low energy structure. Two low energy structures, β-hairpin and helix-coil conformation, are sampled by REMD. In β-hairpin structure, the hydroxyl group (OH) of Thr49 at β turn is accessible to form hydrogen bonds with carboxyl group (COO-) or C=O group of Asp46; while the hydroxyl group (OH) of Thr51 is inclined to interact with water molecule with methyl orienting inwards. As a result, the strong interactions between Asp46 and Thr49 cause a bend at β turn. In helix-coil conformation, hydrophobic interactions play between indole ring and methyl group. One of the zip-in paths of Trpzip4 indicates that, the forming of hydrogen bonds at β turn affects the whole folding process of β hairpin. The forming of hydrogen bonds at β turn takes place first, since it has advantages of distance and strong interactions in both side chains and backbone. Total solvent-accessible surface area decreases significantly as it folds to the β hairpin structure. REMD method could sample large phase space for the low energy minima, and present scenarios of potential energy distribution on phase space with related variables as hydrophobic SASA, backbone RMSD, and hydrogen bonds for specific conformation. This work sheds some lights on the understanding of β hairpin folding.

参考文献

[1] Bryngelson, J. D.; Onuchic, J. N.; Socci, N. D.; Wolynes, P. G. Proteins 1995, 21, 167.

[2] Karplus, M.; Šali, A. Curr. Opin. Struct. Biol. 1995, 5, 58.

[3] Bian, L.; Yang, X.; Liu, L. Acta Chim. Sinica 2005, 63, 1081. (边六交, 杨晓燕, 刘莉, 化学学报, 2005, 63, 1081.)

[4] Yang, F.; Liang, Y.; Yang, F. J. Acta Chim. Sinica 2003, 61, 803 (杨芳, 梁毅, 杨芳(小), 化学学报, 2003, 61, 803.)

[5] Zhou, J.; Thorpe, I. F.; Izvekov, S.; Voth, G. A. Biophys. J. 2007, 92, 4289.

[6] Thorpe, I. F.; Zhou, J.; Voth, G. A. J. Phys. Chem. B 2008, 112, 13079.

[7] Dobson, C. M. Nature 2003, 426, 884.

[8] Hansmann, U. H. E. Chem. Phys. Lett. 1997, 281, 140.

[9] Kinnear, B. S.; Jarrold, M. F.; Hansmann, U. H. E. J. Mol. Graphics Modell. 2004, 22, 397.

[10] Hansmann, U. H. E. J. Chem. Phys. 2004, 120, 417.

[11] Hansmann, U. H. E.; Wille, L. T. Phys. Rev. Lett. 2002, 88, 68105.

[12] Sugita, Y.; Okamoto, Y. Chem. Phys. Lett. 1999, 314, 141.

[13] Zhang, J.; Qin, M.; Wang, W. Proteins: Struct., Funct., Bioinf. 2006, 62, 672.

[14] Munoz, V.; Thompson, P. A.; Hofrichter, J.; Eaton, W. A. Nature 1997, 390, 196.

[15] Chen, C.; Xiao, Y. Bioinformatics 2008, 24, 659.

[16] Zhou, R. Methods Mol. Biol. 2006, 350, 205.

[17] Xie, Y.; Zhou, J.; Jiang, S. Y. J. Chem. Phys. 2010, 132, 065101.

[18] Earl, D. J.; Deem, M. W. Phys. Chem. Chem. Phys. 2005, 7, 3910.

[19] Zhou, R.; Berne, B. J.; Germain, R. Proc. Natl. Acad. Sci. U. S. A. 2001, 98, 14931.

[20] Zhou, R. Proc. Natl. Acad. Sci. U. S. A. 2003, 100, 13280.

[21] Zhou, R. J. Mol. Graphics Modell. 2004, 22, 451.

[22] Haliloglu, T.; Kolinski, A.; Skolnick, J. Biopolymers 2003, 70, 548.

[23] Chen, J.; Won, H. S.; Im, W.; Dyson, H. J.; Brooks, C. L. J. Biomol. NMR. 2005, 31, 59.

[24] La Penna, G.; Mitsutake, A.; Masuya, M.; Okamoto, Y. Chem. Phys. Lett. 2003, 380, 609.

[25] Vreede, J.; Crielaard, W.; Hellingwerf, K. J.; Bolhuis, P. G. Biophys. J. 2005, 88, 3525.

[26] Jas, G. S.; Kuczera, K. Biophys. J. 2004, 87, 3786.

[27] Gnanakaran, S.; Hochstrasser, R. M.; García, A. E. Proc. Natl. Acad. Sci. U. S. A. 2004, 101, 9229.

[28] Antoine, R.; Compagnon, I.; Rayane, D.; Broyer, M.; Dugourd, P.; Breaux, G.; Hagemeister, F. C.; Pippen, D.; Hudgins, R. R.; Jarrold, M. F. J. Am. Chem. Soc. 2002, 124, 6737.

[29] Chebaro, Y.; Mousseau, N.; Derreumaux, P. J. Phys. Chem. B 2009, 113, 7668.

[30] Cecchini, M.; Rao, F.; Seeber, M.; Caflisch, A. J. Chem. Phys. 2004, 121, 10748.

[31] Lu, Y.; Wei, G.; Derreumaux, P. J. Phys. Chem. B 2011, 115, 1282.

[32] Cochran, A. G.; Skelton, N. J.; Starovasnik, M. A. Proc. Natl. Acad. Sci. U. S. A. 2001, 98, 5578.

[33] Du, D.; Zhu, Y.; Huang, C. Y.; Gai, F. Proc. Natl. Acad. Sci. U. S. A. 2004, 101, 15915.

[34] Juraszek, J. Ph.D. Dissertation, University of Amsterdam, Amsterdam, 2008.

[35] Olsen, K. A.; Fesinmeyer, R. M.; Stewart, J. M.; Andersen, N. H. Proc. Natl. Acad. Sci. U. S. A. 2005, 102, 15483.

[36] Wei, H.; Shao, Q.; Gao, Y. Q. Phys. Chem. Chem. Phys. 2010, 12, 9292.

[37] Rathore, N.; Chopra, M.; De Pablo, J. J. J. Chem. Phys. 2005, 122, 024111.

[38] Rao, F.; Caflisch, A. J. Chem. Phys. 2003, 119, 4035.

[39] Zhang, W.; Wu, C.; Duan, Y. J. Chem. Phys. 2005, 123, 154105.

[40] Sindhikara, D.; Meng, Y.; Roitberg, A. E. J. Chem. Phys. 2008, 128, 024103.

[41] Sindhikara, D. J.; Emerson, D. J.; Roitberg, A. E. J. Chem. Theory Comput. 2010, 6, 2804.

[42] Katzgraber, H. G.; Trebst, S.; Huse, D. A.; Troyer, M. J. Stat. Mech.: Theory Exp. 2006, 2006, P03018.

[43] van der Spoel, D.; Lindahl, E.; Hess, B.; van Buuren, A. R.; Apol, E.; Meulenhoff, P. J.; Tieleman, D. P.; Sijbers, A. L. T. M.; Feenstra, K. A.; van Drunen, R.; Berendsen, H. J. C. Gromacs User Manual, version 4.5, 2010, www.gromacs.org.

[44] Darden, T.; York, D.; Pedersen, L. J. Chem. Phys. 1993, 98, 10089.

[45] Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. G. J. Chem. Phys. 1995, 103, 8577.

[46] Berendsen, H. J. C.; Postma, J. P. M.; Van Gunsteren, W. F.; DiNola, A.; Haak, J. J. Chem. Phys. 1984, 81, 3684.

[47] Shaytan, A. K.; Shaitan, K. V.; Khokhlov, A. R. Biomacromolecules 2009, 10, 1224.

[48] Marsh, J. A.; Teichmann, S. A. Structure 2011, 19, 859.

[49] Lins, L.; Thomas, A.; Brasseur, R. Protein Sci. 2003, 12, 1406.

[50] Periole, X.; Mark, A. E. J. Chem. Phys. 2007, 126, 014903.

[51] Jang, S.; Kim, E.; Pak, Y. Proteins 2007, 66, 53.

[52] Blanco, F. J.; Rivas, G.; Serrano, L. Nat. Struct. Mol. Biol. 1994, 1, 584. Zagrovic, B.; Sorin, E. J.; Pande, V. J. Mol. Biol. 2001, 313, 151.
文章导航

/