I have simulated a protein-lipid system as cg with gromacs. Before running , 10 proteins were placed on top of the membrane with the gmx insert command. I checked PDBs. All were right. The simulation was run for 1.5 µs. At the end, PBC was removed in different ways. For example, the following steps were run.
gmx trjconv -f production.xtc -s production.tpr -o production_whole.xtc -pbc whole
gmx trjconv -f production_whole.xtc -s production.tpr -o production_whole_nojump.xtc -pbc nojump
gmx trjconv -f production_whole_nojump.xtc -s production.tpr -o production_whole_nojump_center.xtc -pbc mol -center -n index.ndx
Finally, the structure result from last frame was checked, The result showed that some peptids are placed in inner and some in the outer part of membrane. I think it's because of the rotation. But with the following command, I could not remove rotaion
gmx trjconv -f production_whole.xtc -s production.tpr -o production_whole_nojump.xtc -fit rotatexy
So any one have suggestion to remove rotation? It is appreciated for any help.
The 3D representation of lipid-protein from last frame have been attached.