[gmx-users] restraints failure problem

[Sadaf Rani]

Dear Gromacs users

I am doing a test calculation in MD simulation between ligand and a protein
residue for bond angle and dihedral restraints. I have set them in my
topology file as below:-
; distance restraints
[ bonds ]
;    i        j     type     r0A     r1A     r2A    fcA    r0B     r1B
r2B    fcB
    71    29    10     0.403   0.403   10.0   0.0    0.403   0.403   10.0
4184.000

[ angle_restraints ]
;   ai    aj    ak    al  type    thA      fcA    multA  thB      fcB
 multB
    29    71    69    71   1  104.92   0.0    1     104.92   41.840   1
    71    29    31    29   1  147.02   0.0    1     147.02   41.840   1

[ dihedral_restraints ]
;   ai    aj    ak    al   type    phiA     dphiA  fcA    phiB      dphiB
 fcB
    69    71    29    31   1    146.72     0.0    0.0   146.72     0.0
41.840
    72    71 29    31   1     17.55      0.0    0.0    17.55     0.0
41.840
    72    71 29    25   1     173.10     0.0    0.0   173.10     0.0
41.840

However, after energy minimization, I observed that my restraints are not
working and I find the following change:-

bonds changed from 0.403 *0.413*
angles:-
104.92   changed to *109.24*
147.02  changed to *147.65*
dihedrals:-
 146.72  changed to *133.39*
 17.55  changed to *12.94*
173.10  changed to *-166.16*
How should I work with these so that they may not change much? Dihedrals
seem to be changing more than bonds and angles. I have to apply these in
the free energy calculation which already crashed because of restraint
failure.
I would really be thankful for your kind suggestions, please.

Thanks.
Sadaf