Hi, everyone,
I did a short run (<1ns) for a system of a small peptide in TIP3P water box using both pmemd and sander7 and found that the restart structures from equilibration (minimization followed by 75ps md) by pmemd and sander7 were quite different. The difference was maintained throughout production runs. Pmemd did not disrupt the molecule as much as sander7 when RMSDs are considered. But when examining individual residues, I found that for some residues, (Gly-Arg, in a beta-turn,) pmemd set disallowed (in Ramachandran plot) phi-psi values while phi-psi’s from sander7 are more allowed. By the end of the simulation, those disallowed phi-psi’s were not changed yet.
My questions:
Is it normal that some torsion angles are disallowed throughout simulation?
What is the possible cause of the inconsistency between pmemd and sander7 for my system?
The system was created with ff99 (Amber7) and temp=300 in production runs. The input files are identical. The starting structure (before minimization) was resolved from NMR spectra, without disallowed torsion angles.
Thanks for your help.
Chen
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sat Aug 23 2003 - 06:53:01 PDT
