The investigation of the secondary structure propensities and free-energy landscapes of peptide ligands by replica exchange molecular dynamics simulations

Demir K., Kilic N., DUDAK F. C., BOYACI İ. H., YAŞAR F.

MOLECULAR SIMULATION, vol.40, no.13, pp.1015-1025, 2014 (SCI-Expanded) identifier identifier

  • Publication Type: Article / Article
  • Volume: 40 Issue: 13
  • Publication Date: 2014
  • Doi Number: 10.1080/08927022.2013.835486
  • Journal Indexes: Science Citation Index Expanded (SCI-EXPANDED), Scopus
  • Page Numbers: pp.1015-1025
  • Hacettepe University Affiliated: Yes


The conformational states of two peptide sequences that bind to staphylococcal enterotoxin B are sampled by replica exchange molecular dynamic (REMD) simulations in explicit water. REMD simulations were treated with 52 replicas in the range of 280-501K for both peptides. The conformational ensembles of both peptides are dominated by random coil, bend and turn structures with a small amount of helical structures for each temperature. In addition, while an insignificant presence of beta-bridge structures were observed for both peptides, the beta-sheet structure was observed only for peptide 3. The results obtained from simulations at 300K are consistent with the experimental results obtained from circular dichroism spectroscopy. From the analysis of REMD results, we also calculated hydrophobic and hydrophilic solvent accessible surface areas for both peptides, and it was observed that the hydrophobic segments of the peptides tend to form bend or turn structures. Moreover, the free-energy landscapes of both peptides were obtained by principal component analysis to understand how the secondary structural properties change according to their complex space. From the free-energy analysis, we have found several minima for both peptides at decreased temperature. For these obvious minima of both peptides, it was observed that the random coil, bend and turn structures are still dominant and the helix, beta-bridge or beta-sheet structures can appear or disappear with respect to minima. On the other hand, when we compare the results of REMD with conventional MD simulations for these peptides, the configurations of peptide 3 might be trapped in energy minima during the conventional MD simulations. Hence, it can be said that the REMD simulations have provided a sufficiently high sampling efficiency.