[AMBER] How to define COM for groups in ASMD

From: Matthew Guberman-Pfeffer <matthew.guberman-pfeffer.uconn.edu>
Date: Sat, 6 Mar 2021 03:12:38 -0500

Dear Amber Community,

Similar to a recent post, I want to pull apart two proteins by applying a distance restraint to their centers of mass. I know that a distance restraint of the form

&rst
iat=-1,-1
r2=15.,
r2a=20.
rk2 = 5000
igr1=...
igr2=…
/

would apply to the coordinate-averaged positions of the groups. Would a work around to instead use COMs be (I’m not sure what to do with the iat variable):
&rst
restraint = "distance( COM(...) COM(...) )”
/

Or, is there no way to implement this?

Best regards,
Matthew

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sat Mar 06 2021 - 00:30:03 PST
Custom Search