From: Mert Gür (gurmert_at_gmail.com)
Date: Fri Feb 13 2009 - 13:34:58 CST
This can be done by applying an constraining force. Altough this bias must
somehow be subtracted when the energy is calculated.
Lets first generate and reference file for the constrained atoms by setting
their occupancy column to the spring constant of the force.
set all [atomselect top all]
$all set occupancy 0
set A [atomselect top "protein and backbone"]
$A set occupancy 20
$all writepdb refres_20.pdb
and then adding these to your configuration file
I hope this helps.
On Thu, Feb 12, 2009 at 7:39 PM, Osman Yogurtcu <karmatech_at_yahoo.co.uk>wrote:
> I have a 300 K water NPT simulation of a ligand (one aminoacid) bound to a
> receptor. After I do minimization and equilibration (for 1 ns), I do an MD
> run. I observe that my ligand leaves the binding pocket after 2 ns.
> My question: is there a trick that I can make sure that my ligand sits in
> the binding pocket all the time or should I define bonds artificially and
> fix the ligand?
> thank you,
This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:52:22 CST