On Wed, 2014-04-30 at 09:42 +0000, Valentina Romano wrote:
> Dear Amber users
>
> I want to run a minimization of a protein binds to an aromatic planar molecule (2 aromatic rings fused).
>
> I would like to keep the molecule planar during the minimization applying restraints on it, but i have not experience in it.
> Do you have any suggestions?
After having a conversation with one of my colleagues, it turns out that
purine planarity is a known problem in nucleic acid simulations. If you
make the improper torsions you were modifying even stronger (by
increasing the force constant), does it solve your planarity problem?
It should still be slightly flexible (so don't make the force constant
*too* stiff), but apparently not as flexible as you are seeing.
Personally I think if you are going to force planarity with a restraint,
you might as well use an improper torsion to do it instead. That said,
you can use the NMR restraints (chapter 23 of the Amber 14 manual, and
see the test case in $AMBERHOME/test/nmropt/pme/nmropt_1_torsion for an
example).
HTH,
Jason
--
Jason M. Swails
BioMaPS,
Rutgers University
Postdoctoral Researcher
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Apr 30 2014 - 05:00:03 PDT
