We report the results of protein folding （219M, C34, N36, 2KES, 2KHK） by the method of accelerated molecular dynamics （aMD） at room temperature with the implicit solvent model. Starting from the linear structures,...We report the results of protein folding （219M, C34, N36, 2KES, 2KHK） by the method of accelerated molecular dynamics （aMD） at room temperature with the implicit solvent model. Starting from the linear structures, these proteins successfully fold to the native structure in a lO0-ns aMD simulation. In contrast, they are failed under the traditional MD simulation in the same simulation time. Then we find that the lowest root mean square deviations of helix structures from the native structures are 0.36 A, 0.63 A, 0.52 A, 1.1 A and 0.78 A. What is more, native contacts, cluster and free energy analyses show that the results of the aMD method are in accordance with the experiment very well. All analyses show that the aMD can accelerate the simulation process, thus we may apply it to the field of computer aided drug designs.展开更多
All-atom molecular dynamics （MD） simulations and chemical shifts were used to study interactions and structures in the glycine-water system. Radial distribution functions and the hydrogen-bond network were applied i...All-atom molecular dynamics （MD） simulations and chemical shifts were used to study interactions and structures in the glycine-water system. Radial distribution functions and the hydrogen-bond network were applied in MD simulations. Aggregates in the aqueous glycine solution could be classified into different regions by analysis of the hydrogen-bonding network. Temperature-dependent NMR spectra and the viscosity of glycine in aqueous solutions were measured to compare with the results of MD simulations. The variation tendencies of the hydrogen atom chemical shifts and viscosity with concentration of glycine agree with the statistical results of hydrogen bonds in the MD simulations.展开更多
The hybrid atomistic structure-based model has been validated to be effective in investigation of G-quadruplex folding.In this study,we performed large-scale conventional all-atom simulations to complement the folding...The hybrid atomistic structure-based model has been validated to be effective in investigation of G-quadruplex folding.In this study,we performed large-scale conventional all-atom simulations to complement the folding mechanism of human telomeric sequence Htel24 revealed by a multi-basin hybrid atomistic structure-based model.Firstly,the real time-scale of folding rate,which cannot be obtained from the structure-based simulations,was estimated directly by constructing a Markov state model.The results show that Htel24 may fold as fast as on the order of milliseconds when only considering the competition between the hybrid-1 and hybrid-2 G-quadruplex conformations.Secondly,in comparison with the results of structure-based simulations,more metastable states were identified to participate in the formation of hybrid-1 and hybrid-2 conformations.These findings suggest that coupling the hybrid atomistic structure-based model and the conventional all-atom model can provide more insights into the folding dynamics of DNA G-quadruplex.As a result,the multiscale computational framework adopted in this study may be useful to study complex processes of biomolecules involving large conformational changes.展开更多
Ionic liquids (ILs) have been widely used in separation, catalysis, electrochemistry, etc., and one of the most outstanding characteristics is that ILs can be tailored and tuned for specific tasks. In order to design ...Ionic liquids (ILs) have been widely used in separation, catalysis, electrochemistry, etc., and one of the most outstanding characteristics is that ILs can be tailored and tuned for specific tasks. In order to design and make better use of ionic liquids, the structures and properties relationship is indispensable. Both molecular dynamics and Monte Carlo simulations have been proved useful to understand the behavior of molecules at the microscale and the properties of the system. However, the quality of such simulations depends on force field parameters describing the interactions between atoms. All-atom (AA) or the united-atom (UA) force fields will be chosen because of the demand for more exact results or the lower computational cost, respectively. In order to make a systematic comparison of the two force fields, molecular simulations for four kinds of acyclic guanidinium-based ionic liquids (cations: (R2N)2C=N+<, anion: nitric or perchloric acid) were performed based on the AA and the UA force fields in this work. AA force field parameters were derived from our previous work (Fluid Phase Equilib., 2008, 272: 1-7), and the UA parameters were proposed in this work. Molecular dynamics simulation results for the AA and UA force fields were compared. Simulation densities are very similar to each other. Center of mass radial distribution functions (RDFs), site to site RDFs and spatial distribution functions (SDFs) were also investigated to depict the microscopic structures of the ILs.展开更多
为了准确预测全氮材料的生成焓,基于原子化反应,采用B3PW91、B3P86、B3LYP、X3LYP、O3LYP、M052X、M062X、M06HF、B2PLYP 9种密度泛函分别对52种多氮化合物的气相生成焓进行了计算。通过与实验数据对比,双杂化泛函B2PLYP的平均绝对偏差...为了准确预测全氮材料的生成焓,基于原子化反应,采用B3PW91、B3P86、B3LYP、X3LYP、O3LYP、M052X、M062X、M06HF、B2PLYP 9种密度泛函分别对52种多氮化合物的气相生成焓进行了计算。通过与实验数据对比,双杂化泛函B2PLYP的平均绝对偏差最小,为30.1 k J·mol^(-1),且优于G4方法,选择该泛函计算了N_4(T_d)、_N6(D_(3h))、N_8(O_h)、N_(10)(D_(5h))及N_(12)(D_(6h))5种全氮分子的气相生成焓,计算结果依次为756.4,1338.2,1878.5,2144.3,2787.0 k J·mol^(-1)。展开更多
基金Supported by the National Natural Science Foundation of China under Grant Nos 31200545,11274206 and 11574184
文摘We report the results of protein folding （219M, C34, N36, 2KES, 2KHK） by the method of accelerated molecular dynamics （aMD） at room temperature with the implicit solvent model. Starting from the linear structures, these proteins successfully fold to the native structure in a lO0-ns aMD simulation. In contrast, they are failed under the traditional MD simulation in the same simulation time. Then we find that the lowest root mean square deviations of helix structures from the native structures are 0.36 A, 0.63 A, 0.52 A, 1.1 A and 0.78 A. What is more, native contacts, cluster and free energy analyses show that the results of the aMD method are in accordance with the experiment very well. All analyses show that the aMD can accelerate the simulation process, thus we may apply it to the field of computer aided drug designs.
基金Supported by the National Natural Science Foundation of China(No.20903026)the Talents Introduction Foundation for Universities of Guangdong Province(2011)Scientific Research Foundation of the Natural Science Foundation of Guangdong Province(No.S2011010002483)
文摘All-atom molecular dynamics （MD） simulations and chemical shifts were used to study interactions and structures in the glycine-water system. Radial distribution functions and the hydrogen-bond network were applied in MD simulations. Aggregates in the aqueous glycine solution could be classified into different regions by analysis of the hydrogen-bonding network. Temperature-dependent NMR spectra and the viscosity of glycine in aqueous solutions were measured to compare with the results of MD simulations. The variation tendencies of the hydrogen atom chemical shifts and viscosity with concentration of glycine agree with the statistical results of hydrogen bonds in the MD simulations.
基金the National Natural Science Foundation of China(Grant Nos.11504043,61671107,31670727,and 61771093)the Science Foundation of Shandong Province of China(Grant No.ZR2016JL027)+1 种基金the Taishan Young Scholars Program of Shandong Province of China(Grant No.tsqn20161049)the Youth Science and Technology Innovation Plan of Universities in Shandong,China(Grant No.2019KJE007)。
文摘The hybrid atomistic structure-based model has been validated to be effective in investigation of G-quadruplex folding.In this study,we performed large-scale conventional all-atom simulations to complement the folding mechanism of human telomeric sequence Htel24 revealed by a multi-basin hybrid atomistic structure-based model.Firstly,the real time-scale of folding rate,which cannot be obtained from the structure-based simulations,was estimated directly by constructing a Markov state model.The results show that Htel24 may fold as fast as on the order of milliseconds when only considering the competition between the hybrid-1 and hybrid-2 G-quadruplex conformations.Secondly,in comparison with the results of structure-based simulations,more metastable states were identified to participate in the formation of hybrid-1 and hybrid-2 conformations.These findings suggest that coupling the hybrid atomistic structure-based model and the conventional all-atom model can provide more insights into the folding dynamics of DNA G-quadruplex.As a result,the multiscale computational framework adopted in this study may be useful to study complex processes of biomolecules involving large conformational changes.
基金This work was supported by the National Natural Science Foundation of China （No.20903026）, the Doctoral Scientific Research Foundation of the Natural Science Foundation of Guangdong Province （No.7301567）, and the Research Foundation of Guangdong Pharmaceutical University （No.2006YKX05）.
基金This work was supported by the Doctoral Scientific Research Foundation of the Natural Science Foundation of Guangdong Province, China （No.7301567） and the Research Foundation of Guangdong Pharmaceutical University, China （No.2006YKX05）.
基金supported by the General Program Youth of National Natural Science Foundation of China (20903098, 21073194, 21106146)State Key Laboratory of Multiphase Complex Systems (MPCS-2011-D-05)
文摘Ionic liquids (ILs) have been widely used in separation, catalysis, electrochemistry, etc., and one of the most outstanding characteristics is that ILs can be tailored and tuned for specific tasks. In order to design and make better use of ionic liquids, the structures and properties relationship is indispensable. Both molecular dynamics and Monte Carlo simulations have been proved useful to understand the behavior of molecules at the microscale and the properties of the system. However, the quality of such simulations depends on force field parameters describing the interactions between atoms. All-atom (AA) or the united-atom (UA) force fields will be chosen because of the demand for more exact results or the lower computational cost, respectively. In order to make a systematic comparison of the two force fields, molecular simulations for four kinds of acyclic guanidinium-based ionic liquids (cations: (R2N)2C=N+<, anion: nitric or perchloric acid) were performed based on the AA and the UA force fields in this work. AA force field parameters were derived from our previous work (Fluid Phase Equilib., 2008, 272: 1-7), and the UA parameters were proposed in this work. Molecular dynamics simulation results for the AA and UA force fields were compared. Simulation densities are very similar to each other. Center of mass radial distribution functions (RDFs), site to site RDFs and spatial distribution functions (SDFs) were also investigated to depict the microscopic structures of the ILs.
文摘为了准确预测全氮材料的生成焓,基于原子化反应,采用B3PW91、B3P86、B3LYP、X3LYP、O3LYP、M052X、M062X、M06HF、B2PLYP 9种密度泛函分别对52种多氮化合物的气相生成焓进行了计算。通过与实验数据对比,双杂化泛函B2PLYP的平均绝对偏差最小,为30.1 k J·mol^(-1),且优于G4方法,选择该泛函计算了N_4(T_d)、_N6(D_(3h))、N_8(O_h)、N_(10)(D_(5h))及N_(12)(D_(6h))5种全氮分子的气相生成焓,计算结果依次为756.4,1338.2,1878.5,2144.3,2787.0 k J·mol^(-1)。
