Dear Justin,<br><br>Thanks for the suggestions.<br><br>will be back aftersome time as i am redoing it. Please check me for the accuracy in steps in near future.<br> <br><br>Ram<br><br><div class="gmail_quote">On Thu, Sep 17, 2009 at 8:20 PM, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="im"><br>
<br>
ram bio wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Dear Gromacs Users,<br>
<br>
While I was runing the Justin's tutorial that is KALP-15 in DPPC , I have some queries regarding the INFLATEGRO, Please have patience to read and answer my queries:<br>
<br>
1) does the system.gro should include the position restrained file of KALP_newbox.gro after running genrestr, as i combined the ouput of genestr to dppc128.gro to generate system.gro.<br>
<br>
</blockquote>
<br></div>
You don't include a position restraint file in a coordinate file :)<br>
<br>
I give explicit instructions on what to do in section 2 of the following page:<br>
<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/membrane_protein/03_solvate.html" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/membrane_protein/03_solvate.html</a><div class="im">
<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
2) should both the box vectors in system.gro to be deleted or just the ones under KALP structure, the other being at the end of the file.<br>
<br>
</blockquote>
<br></div>
Just the ones in the KALP structure. Box vectors only ever belong on the very last line of the .gro file.<div><div></div><div class="h5"><br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
3)then I updated the system.gro file to 17503 atoms from 138 atoms and executed this command: perl INFLATEGRO system.gro 4 DPP 14 inflated_system.gro 5 areaperlipid.dat<br>
<br>
the output was as under:<br>
Reading.....<br>
Scaling lipids....<br>
There are 128 lipids...<br>
with 50 atoms per lipid..<br>
<br>
Determining upper and lower leaflet...<br>
64 lipids in the upper...<br>
64 lipids in the lower leaflet<br>
<br>
Centering protein....<br>
Checking for overlap....<br>
...this might actually take a while....<br>
100 % done...<br>
There are 2 lipids within cut-off range...<br>
1 will be removed from the upper leaflet...<br>
1 will be removed from the lower leaflet...<br>
<br>
Writing scaled bilayer & centered protein...<br>
<br>
<br>
Calculating Area per lipid...<br>
Protein X-min/max: 26 40<br>
Protein Y-min/max: 25 39<br>
X-range: 14 A Y-range: 14 A<br>
Building 14 X 14 2D grid on protein coordinates...<br>
Calculating area occupied by protein..<br>
full TMD..<br>
upper TMD....<br>
lower TMD....<br>
Area per protein: 2 nm^2<br>
Area per lipid: 10.4716089904762 nm^2<br>
<br>
Area per protein, upper half: 1.75 nm^2<br>
Area per lipid, upper leaflet : 10.4755772444444 nm^2<br>
<br>
Area per protein, lower half: 2 nm^2<br>
Area per lipid, lower leaflet : 10.4716089904762 nm^2<br>
<br>
Writing Area per lipid...<br>
Done!<br>
<br>
then, i updated the topol_dppc_top to 126 lipids<br>
<br>
; System specifications<br>
[ system ]<br>
126-Lipid DPPC Bilayer<br>
<br>
[ molecules ]<br>
; molecule name nr.<br>
DPPC 126<br>
SOL 3655<br>
<br>
</blockquote>
<br></div></div>
Well, that's certainly wrong. There's a protein in there! As well, inflategro deletes water molecules.<br>
<br>
Follow closely the procedure I describe here:<br>
<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/membrane_protein/02_topology.html" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/membrane_protein/02_topology.html</a><br>
<br>
Because it appears that your topology is incorrect even before the inflategro step.<div><div></div><div class="h5"><br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
and also updated the number of atoms in system.gro to 17403 and added 126 DPPC and 3655 SOL to topol.top file and executed the grompp command for minimization as under:<br>
<br>
grompp -f minim.mdp -c system.gro -p topol.top -o em.tpr<br>
<br>
minim.mdp parameters are as under:<br>
; minim.mdp - used as input into grompp to generate em.tpr<br>
; Parameters describing what to do, when to stop and what to save<br>
integrator = steep ; Algorithm (steep = steepest descent minimization)<br>
emtol = 1000.0 ; Stop minimization when the maximum force < 1000.0 kJ/mol/nm<br>
emstep = 0.01 ; Energy step size<br>
nsteps = 50000 ; Maximum number of (minimization) steps to perform<br>
define = -DSTRONG_POSRES<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.2 ; Cut-off for making neighbor list (short range forces)<br>
coulombtype = PME ; Treatment of long range electrostatic interactions<br>
rcoulomb = 1.2 ; Short-range electrostatic cut-off<br>
rvdw = 1.2 ; Short-range Van der Waals cut-off<br>
pbc = xyz ; Periodic Boundary Conditions (yes/no)<br>
<br>
the ouput of grompp was :<br>
<br>
Back Off! I just backed up mdout.mdp to ./#mdout.mdp.1#<br>
checking input for internal consistency...<br>
processing topology...<br>
Opening library file /usr/local/gromacs/share/gromacs/top/ff_dum.itp<br>
Generated 813 of the 2346 non-bonded parameter combinations<br>
<br>
ERROR 1 [file topol.top, line 535]:<br>
No default G96Angle types<br>
<br>
</blockquote>
<br></div></div>
The acetyl cap is problematic; for some reason pdb2gmx does not assign a value for this angle. You can look it up in ffG53a6bon.itp<div class="im"><br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Opening library file /usr/local/gromacs/share/gromacs/top/spc.itp<br>
Opening library file /usr/local/gromacs/share/gromacs/top/ions.itp<br>
Excluding 3 bonded neighbours molecule type 'Protein'<br>
Excluding 3 bonded neighbours molecule type 'DPPC'<br>
Excluding 2 bonded neighbours molecule type 'SOL'<br>
<br>
</blockquote>
<br></div>
Now this doesn't make sense with what you've shown above. Protein is recognized in the topology, but that isn't true in the topology snippet from before. Please provide accurate information, or else you won't be able to get any reliable help!<div class="im">
<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
NOTE 1 [file topol.top, line 904]:<br>
System has non-zero total charge: 4.000011e+00<br>
<br>
<br>
processing coordinates...<br>
<br>
WARNING 1 [file topol.top, line 904]:<br>
Bad box in file system.gro<br>
<br>
</blockquote>
<br></div>
This probably comes from the fact that you've deleted the wrong set of box vectors.<div><div></div><div class="h5"><br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Generated a cubic box 6.612 x 6.639 x 6.741<br>
Warning: atom name 6439 in topol.top and system.gro does not match (OW - C1)<br>
Warning: atom name 6440 in topol.top and system.gro does not match (HW1 - C2)<br>
Warning: atom name 6441 in topol.top and system.gro does not match (HW2 - C3)<br>
Warning: atom name 6442 in topol.top and system.gro does not match (OW - N4)<br>
Warning: atom name 6443 in topol.top and system.gro does not match (HW1 - C5)<br>
Warning: atom name 6444 in topol.top and system.gro does not match (HW2 - C6)<br>
Warning: atom name 6445 in topol.top and system.gro does not match (OW - O7)<br>
Warning: atom name 6446 in topol.top and system.gro does not match (HW1 - P8)<br>
Warning: atom name 6447 in topol.top and system.gro does not match (HW2 - O9)<br>
Warning: atom name 6448 in topol.top and system.gro does not match (OW - O10)<br>
Warning: atom name 6449 in topol.top and system.gro does not match (HW1 - O11)<br>
Warning: atom name 6450 in topol.top and system.gro does not match (HW2 - C12)<br>
Warning: atom name 6451 in topol.top and system.gro does not match (OW - C13)<br>
Warning: atom name 6452 in topol.top and system.gro does not match (HW1 - O14)<br>
Warning: atom name 6453 in topol.top and system.gro does not match (HW2 - C15)<br>
Warning: atom name 6454 in topol.top and system.gro does not match (OW - O16)<br>
Warning: atom name 6455 in topol.top and system.gro does not match (HW1 - C17)<br>
Warning: atom name 6456 in topol.top and system.gro does not match (HW2 - C18)<br>
Warning: atom name 6457 in topol.top and system.gro does not match (OW - C19)<br>
Warning: atom name 6458 in topol.top and system.gro does not match (HW1 - C20)<br>
(more than 20 non-matching atom names)<br>
<br>
WARNING 2 [file topol.top, line 904]:<br>
10965 non-matching atom names<br>
atom names from topol.top will be used<br>
atom names from system.gro will be ignored<br>
<br>
</blockquote>
<br></div></div>
This is bad news. The topology is expecting water, but finding lipid. Probably due to some inconsistency that I've mentioned above. Follow the instructions carefully when setting up the topology. That's why I advise that it's one of the most important (and detailed!) processes to go through.<div class="im">
<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
double-checking input for internal consistency...<br>
renumbering atomtypes...<br>
converting bonded parameters...<br>
<br>
There was 1 note<br>
<br>
There were 2 warnings<br>
<br>
-------------------------------------------------------<br>
Program grompp, VERSION 4.0.3<br>
Source code file: grompp.c, line: 986<br>
<br>
Fatal error:<br>
There was 1 error in input file(s)<br>
-------------------------------------------------------<br>
<br>
then I executed mdrun -v -deffnm em<br>
<br>
</blockquote>
<br></div>
Don't. If grompp is throwing that many errors, everything after this is meaningless, especially if it thinks your lipids are water molecules. Never plow ahead when faced with errors; it is a waste of time.<br>
<br>
-Justin<br>
<br>
-- <br>
========================================<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<br>
_______________________________________________<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</blockquote></div><br>