The folding of many small proteins is kinetically a two-state process with one major free-energy barrier to overcome,which can be roughly regarded as the inverse process of unfolding.In this work,we first use a Gaussi...The folding of many small proteins is kinetically a two-state process with one major free-energy barrier to overcome,which can be roughly regarded as the inverse process of unfolding.In this work,we first use a Gaussian network model to predict the folding nucleus corresponding to the major free-energy barrier of protein 2 GB1,and find that the folding nucleus is located in theβ-sheet domain.High-temperature molecular dynamics simulations are then used to investigate the unfolding process of 2 GB1.We draw free-energy surface from unfolding simulations,taking RMSD and contact number as reaction coordinates,which confirms that the folding of 2 GB1 is kinetically a two-state process.The comparison of the contact maps before and after the free energy barrier indicates that the transition from native to non-native structure of the protein is kinetically caused by the destruction of theβ-sheet domain,which manifests that the folding nucleus is indeed located in theβ-sheet domain.Moreover,the constrained MD simulation further confirms that the destruction of the secondary structures does not alter the topology of the protein retained by the folding nucleus.These results provide vital information for upcoming researchers to further understand protein folding in similar systems.展开更多
基金Project supported by the Strategic Priority Research Program of Chinese Academy of Sciences(Grant No.XDA17010504)the National Natural Science Foundation of China(Grant No.11947302)。
文摘The folding of many small proteins is kinetically a two-state process with one major free-energy barrier to overcome,which can be roughly regarded as the inverse process of unfolding.In this work,we first use a Gaussian network model to predict the folding nucleus corresponding to the major free-energy barrier of protein 2 GB1,and find that the folding nucleus is located in theβ-sheet domain.High-temperature molecular dynamics simulations are then used to investigate the unfolding process of 2 GB1.We draw free-energy surface from unfolding simulations,taking RMSD and contact number as reaction coordinates,which confirms that the folding of 2 GB1 is kinetically a two-state process.The comparison of the contact maps before and after the free energy barrier indicates that the transition from native to non-native structure of the protein is kinetically caused by the destruction of theβ-sheet domain,which manifests that the folding nucleus is indeed located in theβ-sheet domain.Moreover,the constrained MD simulation further confirms that the destruction of the secondary structures does not alter the topology of the protein retained by the folding nucleus.These results provide vital information for upcoming researchers to further understand protein folding in similar systems.
