Problem with orientational restraints

GROMACS version: 2020
GROMACS modification: No
Hi everyone, I am facing with one pretty annoying problem while setting up a protein simulation with orientational restraints (RDC data from NMR). I’ve searched everywhere for an answer and it look like on one person had similar problem, but never posted a solution.
I am setting up simulations and for a test I am trying to incorporate some constraints to protein .itp file right after all bonded interactions (I also tried to relocate it to the corresponding place in the topol.top file – same result):
[ orientation_restraints ]
;ai aj type exp. label alpha const. obs. weight
32 33 1 1 3 3 6.083 2.9235 1.00
52 53 1 1 4 3 6.083 9.8683 1.00
55 56 1 1 5 3 6.083 -13.2275 1.00
66 67 1 1 6 3 6.083 -8.1949 1.00
102 103 1 1 7 3 6.083 2.9235 1.00
166 167 1 1 8 3 6.083 3.9235 1.00

after that I add this to mdp file (here is an example for final md run, but I also tried it on any step nvt, npt – same problem):
orire = yes
orire-fc = 1 ; Orientation restraints force constant
orire-tau = 0 ; Turn off time averaging
orire-fitgrp = backbone

I run : gmx mdrun -v -dlb yes -update cpu -deffnm md_0_1

and get this error:
Fatal error:
Found 44630 copies of a molecule with orientation restrains while the current
code only supports a single copy. If you want to ensemble average, run
multiple copies of the system using the multi-sim feature of mdrun.

As 44630 is exact number of my SOL molecules I assume that for some reason the constraints are applied to them, not the protein molecule. Maybe someone can help? I can upload all my file, but those are pretty standard and if I remove “orire= yes” everything works. Thank you for looking.

Please attach all input files for grompp so I can check what the issue is.

Thank you. As a new user it doesn’t allow me to upload a file. One think I discover, it looks like it runs in gromacs 2016.4

2016 does not have the check that gave the error in 2020, so we don’t know if your results are correct or not.
I tried a similar system with protein+water in 2020 and I did not get the error.
You should at least be able to upload the top and itp files. Could you dod that? If that doesn’t work, could you upload the input files for grompp somewhere else and post a link?

please take a look at the files (this is the same system, but has more restraints): https://www.dropbox.com/sh/9ed0ey31inmr7ob/AACvthhmSnDdg7Bdrl-MzwAma?dl=0

I don’t know if you have the same problem with the system on dropbox as with the system you first posted about. To me this seems to be a different issue, but maybe not? The issue here is that there seems to be a check missing in Gromacs and I get an assertion failure about orires type indices not being consecutive. The code is written for a single molecule only, but there is no check for this. You could solve this by merging the two chains into one moleculetype. I can also see if I can fix this issue, but I’m currently very busy.