Hi Kaasi,<br>
<br>
The error message indicates that some residues are incomplete. You have
to first complete these residues and also check for residues which are
missing entirely, which you should add, as was done for the work you're
trying to reproduce. <br>
<br>
By the way, I noticed that you tried to use the gmx force field, which
is a bad choice. Didn't the authors state which force field they used?
Gromos perhaps?<br>
<br>
Good luck,<br>
<br>
Tsjerk<br><br><div><span class="gmail_quote">On 12/25/05, <b class="gmail_sendername">Kaasi Sridhara</b> <<a href="mailto:research.gromacs@gmail.com">research.gromacs@gmail.com</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi Everyone,<br>
This is the first time I am working on proteins in lipid bilayer.<br>
I downloaded 1bl8.pdb( KcsA potassium channel) and tried to do pdb2gmx on this pdb file.<br>
I got the following error:<br>
Fatal Error : Atom CZ not found in residue ARG 95 while adding hydrogens.<br>
I used Gromacs Forcefield (4 option from 0-5).<br>
I tried to look into the archives, I could find only one thread, but which did not finally solve the problem.<br>
<br>
---------<br>
About the pdb file, the residue numbering starts from 23 .<br>
<br>
My aim is to replicate the results obtained from the following journal<br>
Simulations of Ion Permeation Through a Potassium Channel: Molecular
Dynamics of KcsA in a Phospholipid Bilayer by Shrivastava and Sansom.<br>
<br>
-------<br>
They described a bit about the procedure, but I am not favourable with what they were talking.<br>
<br>
----I am just adding directly the lines from the journal----<br>
""The protein structure was that in PDB file 1bl8, with the modification that<br>
atoms for all side chains were included. The "missing" residues were added<br>
by building in stereochemically preferred conformers. Although not<br>
present in the x-ray structure, an inter-subunit salt bridge (D80 to R89) was<br>
formed during the early stages of the simulations. Note that residues 1–22<br>
and 120–160 are missing from this model, which therefore represents the<br>
core channel-forming domain of the protein.""<br>
___________<br>
Did anyone tried to work on this problem earlier, or can anyone suggest
how to deal with this problem of adding missing residues.<br>
<br>
My problem now is to form a force-field file for this protein.<br>
<br>
Thanks a lot, and Wish you all a very happy festive mood......<br><span class="sg">
-Kaasi.<br>
</span><br>_______________________________________________<br>gmx-users mailing list<br><a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br><a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">
http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>Please don't post (un)subscribe requests to the list. Use the<br>www interface or send it to <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users-request@gromacs.org">
gmx-users-request@gromacs.org</a>.<br><br></blockquote></div><br><br clear="all"><br>-- <br><br>Tsjerk A. Wassenaar, M.Sc.<br>Groningen Biomolecular Sciences and Biotechnology Institute (GBB)<br>Dept. of Biophysical Chemistry
<br>University of Groningen<br>Nijenborgh 4<br>9747AG Groningen, The Netherlands<br>+31 50 363 4336<br>