[gmx-users] minimum distance between the protein and its mirror image

YanhuaOuyang 15901283893 at 163.com
Mon Aug 29 11:18:27 CEST 2016


Hi,
   I am running a REMD of a disordered protein, I visualized the trajectory in VMD and I found that the protein is very close to the box edge.
   Then I use "gmx mindist " to  check  if a protein has seen its periodic image during simulation. When I used the command "gmx mindist -f remd_01.pdb -s remd.tpr -od mindist.xvg -pi" I choose the 1 group---protein. I choose some data from  mindist.xvg, they are as follow:
        0        2.800  5.225  7.200  7.200  7.200
        2        2.793  5.136  7.200  7.200  7.200
        4        2.804  5.002  7.200  7.200  7.200
        6        2.777  5.176  7.200  7.200  7.200
        8        2.820  5.187  7.200  7.200  7.200
        10       2.871  5.043  7.200  7.200  7.200
        12       2.788  5.089  7.200  7.200  7.200
        14       2.882  4.892  7.200  7.200  7.200
       ................
        5154     4.153  3.415  7.200  7.200  7.200
        5156     4.222  3.483  7.200  7.200  7.200
        5158     4.154  3.608  7.200  7.200  7.200
        5172     3.607  4.124  7.200  7.200  7.200
        5174     3.556  4.140  7.200  7.200  7.200
        5176     3.303  4.430  7.200  7.200  7.200
        5178     3.291  4.297  7.200  7.200  7.200
        ...................
        5880     1.659  5.595  7.200  7.200  7.200
        5882     1.787  5.564  7.200  7.200  7.200
        5884     1.718  5.575  7.200  7.200  7.200
        5886     1.669  5.654  7.200  7.200  7.200
        5888     1.636  5.752  7.200  7.200  7.200
        5890     1.590  5.761  7.200  7.200  7.200
        5892     1.620  5.786  7.200  7.200  7.200
        5894     1.513  5.791  7.200  7.200  7.200
        5896     1.523  5.908  7.200  7.200  7.200
I set the short-range VDW and electrostatic cutoffs=1.0 nm, the distance between the outside of protein and the edge of box is 1.0 nm. the Minimum distance to periodic image ranges from 1.5 nm to 4.2 nm from the data above. While the deal value should be at least 2.0 nm, which is double the cutoff.
Do anyone knows the simulation is normal and if a protein has seen its periodic image during simulation?
which group should I choose when i using gmx mindist, protein, C-alpha or some else?

Best regards,
Ouyang


More information about the gromacs.org_gmx-users mailing list