Energy minimisation in membrane simulation

GROMACS version:
GROMACS modification: Yes/No
Here post your question

I am trying for a membrane simulation of a 700 residue prot in dppc bilayer. I prepared the system.gro file but when I try to inflate and minimize it, the following warnings appear. What could be wrong, Pls guide me ?

ARNING: Listed nonbonded interaction between particles 11833 and 11837
at distance 6.212 which is larger than the table limit 2.050 nm.

This is likely either a 1,4 interaction, or a listed interaction inside
a smaller molecule you are decoupling during a free energy calculation.
Since interactions at distances beyond the table cannot be computed,
they are skipped until they are inside the table limit again. You will
only see this message once, even if it occurs for several interactions.

IMPORTANT: This should not happen in a stable simulation, so there is
probably something wrong with your system. Only change the table-extension
distance in the mdp file if you are really sure that is the reason.

WARNING: Listed nonbonded interaction between particles 8280 and 8301
at distance 6.245 which is larger than the table limit 2.050 nm.

This is likely either a 1,4 interaction, or a listed interaction inside
a smaller molecule you are decoupling during a free energy calculation.
Since interactions at distances beyond the table cannot be computed,
they are skipped until they are inside the table limit again. You will
only see this message once, even if it occurs for several interactions.

IMPORTANT: This should not happen in a stable simulation, so there is
probably something wrong with your system. Only change the table-extension
distance in the mdp file if you are really sure that is the reason.

Energy minimization has stopped, but the forces have not converged to the
requested precision Fmax < 100 (which may not be possible for your system).
It stopped because the algorithm tried to make a new step whose size was too
small, or there was no change in the energy since last step. Either way, we
regard the minimization as converged to within the available machine
precision, given your starting configuration and EM parameters.

Double precision normally gives you higher accuracy, but this is often not
needed for preparing to run molecular dynamics.

writing lowest energy coordinates.

Back Off! I just backed up system_inflated_em.gro to ./#system_inflated_em.gro.3#

Steepest Descents converged to machine precision in 15 steps,
but did not reach the requested Fmax < 1000.
Potential Energy = 4.9288830e+18
Maximum force = inf on atom 8594
Norm of force = inf

Hi,
I got the same problem from one of my system and I solved this problem by changing constraint= all-bonds instead of none. After that I did energy minimize again constraint= none.
I did this a couple of time.

Regards
Zohurul

Hi,

Thanks for the reply, but I got the warnings while using minim_inflategro.mdp, a step before inflation.

gmx mdrun -deffnm system_inflated_em is the command I executed.

So, Can I proceed further, as in the tutorial (ignoring the warnings) for mem simulation, given by Justin.

Following are the contents of that mdp file.

define = -DSTRONG_POSRES ; Prevent protein from moving
integrator = steep ; Algorithm (steep = steepest descent minimization)
emtol = 1000.0 ; Stop minimization when the maximum force < 1000.0 kJ/mol/nm
emstep = 0.01 ; Energy step size
nsteps = 50000 ; Maximum number of (minimization) steps to perform

; Parameters describing how to find the neighbors of each atom and how to calculate the interactions
nstlist = 1 ; Frequency to update the neighbor list and long range forces
cutoff-scheme = Verlet ; Buffered neighbor searching
ns_type = grid ; Method to determine neighbor list (simple, grid)
rlist = 1.2 ; Cut-off for making neighbor list (short range forces)
coulombtype = cutoff ; Treatment of long range electrostatic interactions
rcoulomb = 1.2 ; Short-range electrostatic cut-off
rvdw = 1.2 ; Short-range Van der Waals cut-off
pbc = xyz ; Periodic Boundary Conditions

How can we remove the overlapping atoms is the query ? Can we use maestro for refining the system.gro or is there any method to do it ?

Hope you got my point.

Regards