Folding of Trp-cage Mini Protein Using Temperature and Biasing Potential Replica—Exchange Molecular Dynamics Simulations
The folding process of the 20 residue Trp-cage mini-protein was investigated using standard temperature replica exchange molecular dynamics (T-RexMD) simulation and a biasing potential RexMD (BP-RexMD) method. In contrast to several conventional molecular dynamics simulations, both RexMD methods sam...
Main Authors: | , |
---|---|
Format: | Article |
Language: | English |
Published: |
MDPI AG
2009-03-01
|
Series: | International Journal of Molecular Sciences |
Subjects: | |
Online Access: | http://www.mdpi.com/1422-0067/10/3/1121/ |
id |
doaj-e497407fbed9490391d7beda4589fbc5 |
---|---|
record_format |
Article |
spelling |
doaj-e497407fbed9490391d7beda4589fbc52020-11-25T00:24:17ZengMDPI AGInternational Journal of Molecular Sciences1422-00672009-03-011031121113710.3390/ijms10031121Folding of Trp-cage Mini Protein Using Temperature and Biasing Potential Replica—Exchange Molecular Dynamics SimulationsSrinivasaraghavan KannanMartin ZachariasThe folding process of the 20 residue Trp-cage mini-protein was investigated using standard temperature replica exchange molecular dynamics (T-RexMD) simulation and a biasing potential RexMD (BP-RexMD) method. In contrast to several conventional molecular dynamics simulations, both RexMD methods sampled conformations close to the native structure after 10-20 ns simulation time as the dominant conformational states. In contrast, to T-RexMD involving 16 replicas the BP-RexMD method achieved very similar sampling results with only five replicas. The result indicates that the BP-RexMD method is well suited to study folding processes of proteins at a significantly smaller computational cost, compared to T-RexMD. Both RexMD methods sampled not only similar final states but also agreed on the sampling of intermediate conformations during Trp-cage folding. The analysis of the sampled potential energy contributions indicated that Trp-cage folding is favored by both van der Waals and to a lesser degree electrostatic contributions. Folding does not introduce any significant sterical strain as reflected by similar energy distributions of bonded energy terms (bond length, bond angle and dihedral angle) of folded and unfolded Trp-cage structures. http://www.mdpi.com/1422-0067/10/3/1121/Conformational samplingmolecular dynamics simulationprotein foldingReplica Exchange. |
collection |
DOAJ |
language |
English |
format |
Article |
sources |
DOAJ |
author |
Srinivasaraghavan Kannan Martin Zacharias |
spellingShingle |
Srinivasaraghavan Kannan Martin Zacharias Folding of Trp-cage Mini Protein Using Temperature and Biasing Potential Replica—Exchange Molecular Dynamics Simulations International Journal of Molecular Sciences Conformational sampling molecular dynamics simulation protein folding Replica Exchange. |
author_facet |
Srinivasaraghavan Kannan Martin Zacharias |
author_sort |
Srinivasaraghavan Kannan |
title |
Folding of Trp-cage Mini Protein Using Temperature and Biasing Potential Replica—Exchange Molecular Dynamics Simulations |
title_short |
Folding of Trp-cage Mini Protein Using Temperature and Biasing Potential Replica—Exchange Molecular Dynamics Simulations |
title_full |
Folding of Trp-cage Mini Protein Using Temperature and Biasing Potential Replica—Exchange Molecular Dynamics Simulations |
title_fullStr |
Folding of Trp-cage Mini Protein Using Temperature and Biasing Potential Replica—Exchange Molecular Dynamics Simulations |
title_full_unstemmed |
Folding of Trp-cage Mini Protein Using Temperature and Biasing Potential Replica—Exchange Molecular Dynamics Simulations |
title_sort |
folding of trp-cage mini protein using temperature and biasing potential replica—exchange molecular dynamics simulations |
publisher |
MDPI AG |
series |
International Journal of Molecular Sciences |
issn |
1422-0067 |
publishDate |
2009-03-01 |
description |
The folding process of the 20 residue Trp-cage mini-protein was investigated using standard temperature replica exchange molecular dynamics (T-RexMD) simulation and a biasing potential RexMD (BP-RexMD) method. In contrast to several conventional molecular dynamics simulations, both RexMD methods sampled conformations close to the native structure after 10-20 ns simulation time as the dominant conformational states. In contrast, to T-RexMD involving 16 replicas the BP-RexMD method achieved very similar sampling results with only five replicas. The result indicates that the BP-RexMD method is well suited to study folding processes of proteins at a significantly smaller computational cost, compared to T-RexMD. Both RexMD methods sampled not only similar final states but also agreed on the sampling of intermediate conformations during Trp-cage folding. The analysis of the sampled potential energy contributions indicated that Trp-cage folding is favored by both van der Waals and to a lesser degree electrostatic contributions. Folding does not introduce any significant sterical strain as reflected by similar energy distributions of bonded energy terms (bond length, bond angle and dihedral angle) of folded and unfolded Trp-cage structures. |
topic |
Conformational sampling molecular dynamics simulation protein folding Replica Exchange. |
url |
http://www.mdpi.com/1422-0067/10/3/1121/ |
work_keys_str_mv |
AT srinivasaraghavankannan foldingoftrpcageminiproteinusingtemperatureandbiasingpotentialreplicaexchangemoleculardynamicssimulations AT martinzacharias foldingoftrpcageminiproteinusingtemperatureandbiasingpotentialreplicaexchangemoleculardynamicssimulations |
_version_ |
1725352875940380672 |