[gmx-developers] How to avoid the molecules being too close together in multiple molecules MD simulation?

[Kaifang Huang]

Dear Developers,
    I want to perform a protein-ligands system MD simulation. Unlike
conventional MD simulations, the ligand were not put into the binding
pocket, but I copied many ligands and put them around the protein to run
MD. But, during MD, I found a fatal problem, that is, many ligand Molecules
tend to clump together.
    I tried to use *[distance restraints]* to control the spacing between
small molecules，which seems to be achieved by forming a bond between a pair
of atoms. However,  this would lead to a paradox that  bond often exceeds
half of the box during the simulation,  then GMX would exit. such as "The
longest distance involved in such interactions is 4.376 nm which is above
half the box length ".
    I've turned to the GROMACS user forum
<https://gromacs.bioexcel.eu/t/how-to-avoid-the-molecules-being-too-close-together-in-multiple-molecules-md-simulation/3765>,
but no one seems to have a good idea.  I sincerely hope that you can give
some advice and always look forward to your reply.




-- 
Best,

Kaifang Huang
New York University
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-developers/attachments/20220408/bd22f747/attachment.html>