A Comparative Study of Two Different Force Fields on Structural and Thermodynamics Character of H1 Peptide via Molecular Dynamics Simulations

The dynamics and thermodynamics character of H1 peptide in aqueous solution has been investigated through temperature replica exchange molecular dynamics (T-REMD) simulations using two different force fields (OPLS-AA and GROMOS 43A1). The two independent T-REMD simulations were completed starting front initial conformations alpha-helix and beta-sheet, respectively. Each replica was run for 300 ns. The performance of each force field was assessed from the parameters such as the distributions of backbone dihedral angles, the number of native hydrogen bond, root mean square deviations (RMSD) of C-alpha atoms and all heavy atoms, formation of beta-turn, the stability of folded beta-hairpin structure and the favorite conformations of different force fields. The simulation using GROMOS 43A1 force field starting from alpha-helix structure sampled the conformation cluster which C-alpha RMSD was 0.05 nm from beta-sheet structure and the cluster contains 39% of all conformations. The simulation using OPLS-AA force field produced more sampling in P-II region than in GROMOS 43A1 force field. The both force field simulations produced some sampling in the alpha region, but the probabilities of the conformations including any helical content were only 1-2%. Under the both force fields, the beta-turn structures exhibited higher stability than alpha-helix structures and the folded beta-hairpin structures. In the GROMOS 43A1 force field, the free energy change from the unfolded state to the hairpin state was in good agreement with the results of several experiments about some beta-peptides (not the H1 peptide) and the other molecular dynamics Simulations of H1 peptide. However, the folded-hairpin structure was more destabilized in the OPLS-AA force field than in the GROMOS 43A1 force field and experiments.