GROMACS version: 2020.4
GROMACS modification: No
I need a way to rotate my entire system (DNA-Protein complex and solvent) such that the axis of the DNA lies along the x-axis, but I want the periodic boundaries to be preserved. Here’s an image of my system:
When I try rotating, I end up with this:
This will cause blow-ups when I run again. How can I edit things in such a way that the entire system is rotated, but molecules stay within the box? (basically the water molecules that were rotated down and went outside the box should reappear up inside the box)
Thank you in advance!
I think the “orient” plugin in VMD is quite useful for this kind of task. There are probably many other ways to do this, but hopefully this points you in the right direction at least.
http://www.ks.uiuc.edu/Research/vmd/script_library/scripts/orient/
Hi, thank you for your reply!
I’m sorry, I should’ve specified that it’s the coordinate file on gromacs that I want modified, not the VMD visualization. I want to be able to continue my simulation on gromacs with the rotated system. Is that possible?
You can save out the new coordinates after aligning in vmd.
This will violate PBC, as once you wrap this rotated system there will be overlapping atoms. A workaround would be, align protein along x axis, and then fill box with water using gmx solvate