Dear Gromacs users,<br><br>I have to simulate a system consisting of Protein+MG+GDP where Protein is in complex with MG and GDP. I generated the GDP parameters using antechamber and then converted the topology,coordinate and parameter files to the corresponding gromacs format. <br>
Then, I tried to start with the simulation using AMBER ff03 forcefield. The system ran well till energy minimization step, but at equilibration steps it is showing lincs error at NVT equilibration step. <br><br>I used 2 coupling groups for my system with lincs constraint algorithm as follows:<br>
<b>tc-grps     = Protein_MG_GDP Water_and_ions </b><br><br>I tried turning off lincs warning too, den NVT equilibration was done finely but the system was blowing up during NPT equilibration. I tried with 3-4 steps of minimization also, but it didnt work (I have tried with lincs warning off in NPT too). I have visualized the NVT output trajectory and compared it with the initial structure too, and the trajectory was fine.<br>
<br><br>I am unable to figure out the reason for why the system is blowing up. I am following the drug-ligand simulation tutorial but still can&#39;t figure out the problem. Can anyone suggest any possible means to overcome this situation?<br>
<br><br>If anyone has simulated any similar system with Protein-ion-ligand complex or have faced similar problem, please share with what to do now to resolve this issue or to simulation such a system? <br>Any help will be highly appreciated.<br>
<br>Thanks<br><br>-----<br>Neeru Sharma,<br>Pune (India)<br><br>