Another approach might be to add distance restraint(s) that kick in when
the ion gets a given distance from given atom(s).
Bill
On 2/19/17 2:41 AM, Elvis Martis wrote:
> 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
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun Feb 19 2017 - 03:00:03 PST
