Hi Justin , <br> Thank you , You are right .<br>Problem get solved by change in spacing in PDB file.<br>have a nice day.<br>With Best Wishes, <br>R.David <br><br><div class="gmail_quote">On Fri, Apr 6, 2012 at 6:36 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="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im"><br>
<br>
rama david wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
HI all ,<br>
sorry for above incomplete information ,<br>
I change CH3 of ACE residue of my pdb file to CA .<br>
<br>
</blockquote>
<br></div>
The formatting of a .pdb file requires fixed spacing. By turning CH3 into CA (which is correct), you deleted a space such that instead of finding 'ACE' as the residue name, pdb2gmx is finding 'CE' since it has been shifted to the left one space. Fix the spacing and you fix your problem.<br>
<br>
-Justin<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Thank you in advance<div><div class="h5"><br>
<br>
On Fri, Apr 6, 2012 at 5:34 PM, rama david <<a href="mailto:ramadavidgroup@gmail.com" target="_blank">ramadavidgroup@gmail.com</a> <mailto:<a href="mailto:ramadavidgroup@gmail.com" target="_blank">ramadavidgroup@gmail.<u></u>com</a>>> wrote:<br>
<br>
Hi GROMACS Friends and Mark..<br>
<br>
Thank you for reply ...<br>
My command line is<br>
pdb2gmx -f ....pdb -o p.pdb -p p.top -ignh -ter<br>
I choose the G96 53a6 ff , along with spc water model.<br>
<br>
Select start terminus type for ACE-1<br>
0: NH3+<br>
1: NH2<br>
2: None<br>
2<br>
Start terminus ACE-1: None<br>
<br>
Select end terminus<br>
0: COO-<br>
1: COOH<br>
2: None<br>
0<br>
End terminus : COO-<br>
<br>
So I got the following out put...<br>
Program pdb2gmx, VERSION 4.5.4<br>
Source code file: /build/buildd/gromacs-4.5.4/<u></u>src/kernel/pdb2gmx.c,<br>
line: 655<br>
<br>
Fatal error:<br>
Atom CH3 in residue ACE 1 was not found in rtp entry ACE with 3 atoms<br>
while sorting atoms.<br>
<br>
.<br>
For more information and tips for troubleshooting, please check the<br>
GROMACS<br>
website at <a href="http://www.gromacs.org/Documentation/Errors" target="_blank">http://www.gromacs.org/<u></u>Documentation/Errors</a><br>
------------------------------<u></u>-------------------------<br>
I have the limited experience for these problem..<br>
Can you tell me what to do in such conditions????<br>
<br>
I change the Atom CH3 in pdb to CA and proceed with the above<br>
command line<br>
and option .<br>
I got the following error<br>
<br>
Program pdb2gmx, VERSION 4.5.4<br>
Source code file: /build/buildd/gromacs-4.5.4/<u></u>src/kernel/resall.c,<br>
line: 581<br>
<br>
Fatal error:<br>
Residue 'CE' not found in residue topology database<br>
<br>
For more information and tips for troubleshooting, please check the<br>
GROMACS<br>
website at <a href="http://www.gromacs.org/Documentation/Errors" target="_blank">http://www.gromacs.org/<u></u>Documentation/Errors</a><br>
------------------------------<u></u>-------------------------<br>
I have the limited experience for these problem..<br>
Can you tell me what to do in such conditions????<br>
<br>
<br>
On Fri, Apr 6, 2012 at 4:38 PM, Mark Abraham<br></div></div><div><div class="h5">
<<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a> <mailto:<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.<u></u>au</a>>> wrote:<br>
<br>
On 6/04/2012 9:05 PM, Mark Abraham wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
On 6/04/2012 8:27 PM, rama david wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Gromacs Friends and Justin ,<br>
<br>
Thank you for reply and suggestion.<br>
<br>
These is short part of my PDB .<br>
<br>
ATOM 1 1H ACE 1 0.000 0.000 0.000<br>
ATOM 2 CH3 ACE 1 0.000 1.090 0.000<br>
ATOM 3 2H ACE 1 1.028 1.453 -0.000<br>
ATOM 4 3H ACE 1 -0.514 1.453 -0.890<br>
ATOM 5 C ACE 1 -0.721 1.600 1.249<br>
ATOM 6 O ACE 1 -0.839 2.806 1.453<br>
ATOM 7 N PRO 2 -1.227 0.727 2.126<br>
ATOM 8 CA PRO 2 -1.918 1.181 3.321<br>
ATOM 9 HA PRO 2 -2.553 2.031 3.073<br>
ATOM 10 C PRO 2 -0.928 1.587 4.404<br>
ATOM 11 O PRO 2 -1.069 2.645 5.013<br>
ATOM 12 CB PRO 2 -2.782 0.075 3.851<br>
ATOM 13 1HB PRO 2 -2.469 -0.177 4.864<br>
ATOM 14 2HB PRO 2 -3.822 0.401 3.862<br>
ATOM 15 CG PRO 2 -2.645 -1.140 2.964<br>
ATOM 16 1HG PRO 2 -2.232 -1.966 3.542<br>
ATOM 17 2HG PRO 2 -3.606 -1.451 2.556<br>
ATOM 18 CD PRO 2 -1.680 -0.755 1.882<br>
ATOM 19 1HD PRO 2 -2.173 -0.830 0.913<br>
ATOM 20 2HD PRO 2 -0.816 -1.419 1.895<br>
ATOM 21 N GLY 3 0.076 0.740 4.644<br>
ATOM 22 H GLY 3 0.135 -0.114 4.109<br>
ATOM 23 CA GLY 3 1.083 1.012 5.650<br>
<br>
The gromacs96 53a6 .rtp file mention the following data...<br>
<br>
bondedtypes ]<br>
; bonds angles dihedrals impropers<br>
2 2 1 2<br>
<br>
[ ACE ]<br>
[ atoms ]<br>
CA CH3 0.000 0<br>
C C 0.450 1<br>
O O -0.450 1<br>
[ bonds ]<br>
C CA gb_27<br>
C O gb_5<br>
C +N gb_19<br>
[ angles ]<br>
CA C O ga_30 CA C +N ga_19<br>
O C +N ga_33<br>
[ impropers ]<br>
C CA +N O gi_1<br>
<br>
[ NH2 ]<br>
[ atoms ]<br>
N NT -0.83 0<br>
H1 H 0.415 0<br>
H2 H 0.415 0<br>
[ bonds ]<br>
N H1 gb_2<br>
N H2 gb_2 -C N gb_9<br>
[ angles ]<br>
-O -C N ga_33<br>
-CA -C N ga_19<br>
-C N H1 ga_23<br>
-C N H2 ga_23<br>
H1 N H2 ga_24<br>
[ dihedrals ]<br>
-CA -C N H1 gd_14<br>
[ impropers ]<br>
-C -O N -CA gi_1<br>
N H1 H2 -C gi_1<br>
As per my pdb file N atom is the part of PROLINE , not the ACE...<br>
I got the topology by using the Amber 03 FF.<br>
But my priority is to the GROMOS96 53a6 force field ...<br>
<br>
Could any suggest me the way to tackle these problem???<br>
</blockquote>
<br>
The default is to try to put termini on the protein chain.<br>
Your chain already has them, so you need to tell pdb2gmx not<br>
to try to add termini. See pdb2gmx -h.<br>
</blockquote>
<br>
... and if you'd also shown the last piece of your pdb2gmx<br>
output and given us the whole command line you used, then you'd<br>
have gotten your answer faster and wasted less of people's time<br>
guessing what the real problem was.<br>
<br>
Mark<br>
<br>
--<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div></div>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><div class="im"><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/<u></u>mailman/listinfo/gmx-users</a><br>
Please search the archive at<br>
<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@<u></u>gromacs.org</a>>.<div class="im"><br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a><br>
<br>
<br>
<br>
</div></blockquote><div class="im HOEnZb">
<br>
-- <br>
==============================<u></u>==========<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<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.<u></u>vt.edu/Pages/Personal/justin</a><br>
<br>
==============================<u></u>==========<br></div><div class="HOEnZb"><div class="h5">
-- <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/<u></u>mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/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/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a><br>
</div></div></blockquote></div><br>