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.

More Javad Zamani's questions See All
Similar questions and discussions