Hi,
You can use ntr = 1 followed by restraint_wt = ?? (high enough to hold the NA ion) and restraintmask = :residue_number_for_that_particular_ion.
And don't forget to define ref = (appropriate reference coordinate file) in your command to submit the production phase.
Best Regards
Elvis Martis
Ph.D. Student (Computational Chemistry)
at Bombay College of Pharmacy
A Kalina, Santacruz [E], Mumbai 400098, INDIA
W www.elvismartis.in
Skype. adrian_elvis12
-----Original Message-----
From: Michael Shokhen [mailto:michael.shokhen.biu.ac.il]
Sent: Sunday, February 19, 2017 3:15 PM
To: AMBER.ambermd.org
Subject: [AMBER] Na cation drifts away from metal chelating protein
Dear Amber experts,
I have on my Linux server Amber16 and AmberTools16.
I am doing molecular dynamics simulation of protein with ff14SB force field and explicit TIP3 water solvent. The initial crystal structure of my protein contains chelating by protein Na cation.
During the production phase MD simulations after 200ns the Na metal cation is moving to the water solvent far away (about 55 A) from its initial position inside protein.
What is the command line that should be added to the input file for the restrained MD simulation to prevent the Na cation drifting away from its chelating complex in protein?
Thank you,
Michael
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun Feb 19 2017 - 03:00:03 PST
