[gmx-users] Change restraints during md

[David van der Spoel]

Michal Kolinski wrote:
> Hello all.
> 
> I wonder is there a way for turning off the distance restraints during 
> the md simulation.
> I’m currently investigation different binding modes of a small
> ligand inside the active site of the receptor protein. 
> I use restraints to force a proper orientation of ligand
> inside the active site during the first stages of the simulation.
> Then I would like to slowly turn off the restraints between protein and 
> small
> molecule (without a rapid change in distance restraints force) to
> investigate the complex stability. 
> Maybe someone knows how to implement
> it in Gromacs code?  Please give me some advice.

You want to modify the force constant from the initial value to zero, 
much like simulated annealing.

This is not implemented but could be done in md.c

the variable you want to change is inputrec->dr_fc


>  Thank you for help.
> Best wishes.
> Michal
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> gmx-users mailing list    gmx-users at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please don't post (un)subscribe requests to the list. Use the 
> www interface or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/mailing_lists/users.php


-- 
David.
________________________________________________________________________
David van der Spoel, PhD, Assoc. Prof., Molecular Biophysics group,
Dept. of Cell and Molecular Biology, Uppsala University.
Husargatan 3, Box 596,  	75124 Uppsala, Sweden
phone:	46 18 471 4205		fax: 46 18 511 755
spoel at xray.bmc.uu.se	spoel at gromacs.org   http://folding.bmc.uu.se
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++