I want to come back to the question of the NMA in the Gromacs :)<br><br>I've found in manual possible algorithms of this analysis- I must calculate Hessian matrix via Md-run and then calculate modes with g_nmens program<br>
<br><br>1- I've performed CG minimization of my initianl structure<br><br>2- I've loadet to the gromp .tpr file of my minimized structure<br><br>3- I've used mdrun and obtain error <br>Fatal error:<br>Constraints present with Normal Mode Analysis, this combination is not supported<br>
<br>So I suppose that the problem with my .mdp file wich consist of the parametres for the NMA of Lyzosyme in water :)<br><br>integrator = nm ; Normal Mode Analysis<br><br>; Parameters describing how to find the neighbors of each atom and how to calculate the interactions<br>
nstlist = 1 ; Frequency to update the neighbor list and long range forces<br>ns_type = grid ; Method to determine neighbor list (simple, grid)<br>rlist = 1.0 ; Cut-off for making neighbor list (short range forces)<br>
coulombtype = PME ; Treatment of long range electrostatic interactions<br>rcoulomb = 1.0 ; Short-range electrostatic cut-off<br>rvdw = 1.0 ; Short-range Van der Waals cut-off<br>pbc = xyz ; Periodic Boundary Conditions (yes/no)<br>
<br><br>How I should edit this configure file for my NMA? <br>If I edit the above parametries for rvdw etc in accardance to the <a href="http://www.gromacs.org/Documentation/How-tos/Normal_Mode_Analysis?highlight=nm">http://www.gromacs.org/Documentation/How-tos/Normal_Mode_Analysis?highlight=nm</a> I've obtained error that<br>
<br>ERROR 1 [file nma.mdp]:<br> With coulombtype = PME, rcoulomb must be equal to rlist<br> If you want optimal energy conservation or exact integration use PME-Switch<br><br><br>ERROR 2 [file nma.mdp]:<br> With vdwtype = Cut-off, rvdw must be >= rlist<br>
<br>What are the most suitable coulombtype for the NMA and combination of the other options?<br><br><br>James<br><br> <br>