# PBC centering Issue **URL:** \ **Category:** User discussions\ **Created:** [January 29, 2026, 12:21pm UTC](https://gromacs.bioexcel.eu/t/pbc-centering-issue/13156 "2026-01-29T12:21:49Z")\ **Posts on this page:** 4\ **Page:** 1
**Author:** ![Anamika](https://avatars.discourse-cdn.com/v4/letter/a/67e7ee/32.png) [@Anamika](https://gromacs.bioexcel.eu/u/Anamika)\ **Post date:** [January 29, 2026, 12:21pm UTC](https://gromacs.bioexcel.eu/t/pbc-centering-issue/13156/1 "2026-01-29T12:21:49Z")
GROMACS version: 2022.3 GROMACS modification: Yes/No Deal all, [md.mdp](https://gromacs.bioexcel.eu/uploads/short-url/8HIqiXiUDgkejfPJpURrvoZVeyr.mdp) (2.6 KB) I am Anamika Kumari, doing PhD in IIPE, Visakhapatnam. I am working with a dimeric and tetrameric protein system containing one and two ligands, respectively, and I have performed a 100 ns molecular dynamics simulation. After completion of the simulation, I observed that the protein-ligand complex was drifting outside the simulation box during visualization, and therefore attempted to remove periodic boundary condition (PBC) artifacts. To address this issue, I applied multiple PBC correction options available in GROMACS, including molecule-based reconstruction, jump removal, and centering. Specifically, I tested combinations of `-pbc mol`, `-pbc nojump`, `-pbc whole`, and `-center` to ensure the integrity and proper positioning of the complex. I also applied PBC corrections using molecule-based and jump-removal options, followed by centering. Specifically, I used the following command: gmx trjconv -s md\_0\_100.tpr -f md\_0\_100.xtc -o md\_noPBC.xtc -pbc nojump -center After applying this correction, the overall protein–ligand complex appeared intact and centered. However, upon visualization, I observed that certain residues appeared to be missing or discontinuous, particularly within the oligomeric protein chains. I have attached md.mdp file and a part of protein where residues disappeared. ![Screenshot 2026-01-29 171804](https://europe1.discourse-cdn.com/flex017/uploads/bioexcel1/original/2X/1/1d10967fa0a9a543c7ad86a50433b9e8155f222d.png) I would appreciate guidance on the correct PBC correction workflow for multimeric protein–ligand systems, and advice on how to prevent residue discontinuities while maintaining structural integrity during post-processing and visualization. ---
**Author:** ![Karis](https://avatars.discourse-cdn.com/v4/letter/k/ed8c4c/32.png) [@Karis](https://gromacs.bioexcel.eu/u/Karis)\ **Post date:** [February 4, 2026, 10:22am UTC](https://gromacs.bioexcel.eu/t/pbc-centering-issue/13156/2 "2026-02-04T10:22:24Z")
If the protein is properly centered without any contact with the boundaries, this is likely due to an issue with the visualization software rather than a PBC issue, the residues are still there but at such an orientation that VMD isn’t visualizing it properly. You could double check by visualizing those residues with lines or CPK rather than NewCartoon. ---
**Author:** ![neilthomson](https://avatars.discourse-cdn.com/v4/letter/n/e99b99/32.png) [@neilthomson](https://gromacs.bioexcel.eu/u/neilthomson)\ **Post date:** [March 25, 2026, 10:26am UTC](https://gromacs.bioexcel.eu/t/pbc-centering-issue/13156/3 "2026-03-25T10:26:20Z")
Are you using pymol to visualise a .gro or .pdb output here? Pymol guesses bonding, so it will usually make proteins look fragmented like this. Try viewing in vmd with a .tpr that contains bond info, or load your .tpr into MDAnalysis to create a pdb file in python. PBC corrections vary depending on system sadly. But, I usually give mine a first attempt with this template - a four step trjconv protocol GROMACS trjconv post-processing pipeline (4 steps): Make molecules whole (-pbc whole) Fixes molecules broken across periodic boundaries Center protein in box (-center -pbc mol -ur compact) Centers the protein and uses compact unit cell representation Cluster & remove jumps (-pbc cluster) Clusters molecules around the centered protein, removing PBC jumps Fit to reference (-fit rot+trans) Removes rotational and translational motion by fitting to the reference structure This works for some systems, but inspect the trajectory afterwards, because sometime the solvent doesn’t get corrected fully ---
**Author:** ![Anamika](https://avatars.discourse-cdn.com/v4/letter/a/67e7ee/32.png) [@Anamika](https://gromacs.bioexcel.eu/u/Anamika)\ **Post date:** [March 25, 2026, 2:56pm UTC](https://gromacs.bioexcel.eu/t/pbc-centering-issue/13156/4 "2026-03-25T14:56:22Z")
Thank you so much for your suggestions. _ **Anamika Kumari** _ _ **UGC-JRF** _ _ **23RS9P10** _ _ **Molecular Modelling & Protein Engineering Lab** _ **_Indian Institute of Petroleum & Energy (IIPE),Visakhapatnam-530003_** ![](https://europe1.discourse-cdn.com/flex017/uploads/bioexcel1/original/2X/0/03161b9b6ee98bb652be532d585e791426fa0199.jpeg)