<div dir="ltr">This does sound useful, though it would be more useful if this could be implemented into the main GROMACS rather than a separate code (since otherwise it will go away as GROMACS is further developed).<br><br>
Would this be a possibility going forward, devels?<br><br>Thanks!<br><div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, Aug 16, 2013 at 3:07 AM, Floris Buelens <span dir="ltr"><<a href="mailto:floris_buelens@yahoo.com" target="_blank">floris_buelens@yahoo.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi David,<br>
<br>
I have the same requirement as you, but I've gone about it slightly differently. I've hacked the couple-moltype code path to allow decoupling of a specific residue (identified by name) instead of a molecule. This allows the residue of interest to be part of another molecule block so you can set up distance, angle and dihedral restraints using regular bonded interactions with full perturbation support.<br>
<br>
If this is useful to you or to anyone else, let me know and I'll be happy to share.<br>
<br>
best,<br>
<br>
Floris<br>
<br>
________________________________<br>
From: David Mobley <<a href="mailto:dmobley@gmail.com">dmobley@gmail.com</a>><br>
<div class="im">To: Discussion list for GROMACS development <<a href="mailto:gmx-developers@gromacs.org">gmx-developers@gromacs.org</a>><br>
</div>Sent: Friday, 14 June 2013, 21:47<br>
Subject: [gmx-developers] free energy calculations with restraints<br>
<div class="HOEnZb"><div class="h5"><br>
<br>
<br>
Hi, devs,<br>
<br>
<br>
I'm writing with an issue relating to the interplay of new free energy features with restraints.<br>
<br>
I'm very much appreciating some of the new free energy features in gromacs, such as the 'couple-moltype' option which provides a way to set up decoupling or annihilation of a specific molecule via free energy calculations without having to edit the topology file directly. This is especially great when it comes to decoupling -- charge decoupling was not previously possible via topology file editing, and vdW decoupling took substantial manipulation of the topology file.<br>
<br>
However, for binding free energies, my work employs orientational restraints between the protein and ligand. I need to be able to impose both dihedral and angle restraints on angles between the protein and ligand. Currently, I do this using angle-restraints and dihedral-restraints. This requires that both the protein and ligand be within the same logical 'molecule', which (unfortunately) means that I can't make use of the new free energy features above, since couple-moltype has to apply to a whole molecule, not just part of a molecule.<br>
<br>
So, my I see two possible solutions to the problem, and hence have these questions:<br>
1) Can dihedral and angle restraints be applied via the pull code? If not, are there any plans to add that?<br>
2) Alternatively, what about modifying the restraints code so it uses (or at least optionally allows) absolute atom numbering, rather than numbering within a specific molecule, thus allowing restraints (angle-restraints and dihedral-restraints) to be applied between 'molecules'?<br>
<br>
Thanks!<br>
David<br>
<br>
<br>
<br>
--<br>
David Mobley<br>
<a href="mailto:dmobley@gmail.com">dmobley@gmail.com</a><br>
<a href="tel:949-385-2436" value="+19493852436">949-385-2436</a><br>
<br>
</div></div><div class="HOEnZb"><div class="h5">--<br>
gmx-developers mailing list<br>
<a href="mailto:gmx-developers@gromacs.org">gmx-developers@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-developers" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-developers</a><br>
Please don't post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-developers-request@gromacs.org">gmx-developers-request@gromacs.org</a>.<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>David Mobley<br><a href="mailto:dmobley@gmail.com" target="_blank">dmobley@gmail.com</a><br>949-385-2436<br>
</div>