GROMACS version:
GROMACS modification: Yes/No
I am testing MD simulation of Tri-complex structure published in PDB. Only A and D chains are used, PDB: 9AX6
I have previously tested other samples through Protein-ligand tutorial and thought it would be feasible.
( tutorial : Protein-Ligand Complex )
However, I get the following error in NVT equilibrium.
Internal error (bug):
Step 0: The total potential energy is 5.84329e+25, which is extremely high.
The LJ and electrostatic contributions to the energy are 5.84329e+25 and
-906647, respectively. A very high potential energy can be caused by
overlapping interactions in bonded interactions or very large coordinate
values. Usually this is caused by a badly- or non-equilibrated initial
configuration, incorrect interactions or parameters in the topology.
Is there anything I need to fix in energy minimization or equilibrium?