The protocol in that paper is completely different than yours.
I strongly recommend you find some more recent papers for DNA
stretching, that one is pretty old.
The key sentence in their paper showing the difference with what you are
doing is in:
"The DNA molecule was stretched by increasing the distance between the
two restrained atoms by 1 Å followed by an intervening 10 ps relaxation
period."
You are pulling continuously without the relaxation step.
Also, what they are plotting is the derivative of the potential energy
vs distance to get the force, while the force column in your table is
the instantaneous force during your pulling.
In your simulation, column 4 is NOT the energy, but rather the WORK,
which is a completely different object.
Adrian
On 3/25/17 2:15 PM, Ramin Salimi wrote:
> Molecular Dynamics Simulations of Duplex Stretching Reveal the Importance of Entropy in Determining the Biomechanical Properties of DNA
--
Dr. Adrian E. Roitberg
University of Florida Research Foundation Professor.
Department of Chemistry
University of Florida
roitberg.ufl.edu
352-392-6972
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sat Mar 25 2017 - 12:00:02 PDT
