<meta charset="utf-8">Olá, Justin, Mark and Osmair. <div>Thanks for help.<br><br></div><meta charset="utf-8"><div>Building dipeptide has not been the problem. Actually I need to use atomic coordinates obtained </div><div>from a previous X-ray crystal structure. For this, I need to give each atom its proper name. </div>
<div>That is, I have this: </div><div>...</div>ATOM 10 C UNK 1 -2.980 6.010 1.900 1.00 0.00 <br>ATOM 11 C UNK 1 -2.380 5.770 1.130 1.00 0.00 <br>ATOM 12 H UNK 1 -4.050 7.000 1.400 1.00 0.00 <br>
ATOM 13 H UNK 1 -3.590 7.740 0.930 1.00 0.00 <br>ATOM 14 H UNK 1 -4.620 6.540 0.740 1.00 0.00 <br>ATOM 15 C UNK 1 -4.930 7.590 2.470 1.00 0.00 <div>
...</div><div><br></div><div>and I need something like this:</div><div>...<br><div>ATOM 10 CG PHE 1 -2.980 6.010 1.900 1.00 0.00</div><div>ATOM 11 CD1 PHE 1 -2.380 5.770 1.130 1.00 0.00</div>
<div>ATOM 12 HD1 PHE 1 -4.050 7.000 1.400 1.00 0.00</div><div>ATOM 13 CD2 PHE 1 -3.590 7.740 0.930 1.00 0.00</div><div>ATOM 14 HD2 PHE 1 -4.620 6.540 0.740 1.00 0.00</div>
<div>ATOM 15 CE1 PHE 1 -4.930 7.590 2.470 1.00 0.00</div><div>...</div><div><br></div><div>Of course if I try to use pdb2gmx without assigning these names, the execution fails. </div><div>I know this should be simple and that there should be several ways to revolve, but it's puzzling me.</div>
<div><br></div><div>Is there any software or script I can use to do that automatically? <br></div><div>Any hint is very welcome!</div><div><br></div><div><br></div><div>Até mais<br>eef</div><div><br></div><div>_______________________________________<br>
Eudes Eterno Fileti<br>Centro de Ciências Naturais e Humanas<br>Universidade Federal do ABC — CCNH<br>Av. dos Estados, 5001<br>Santo André - SP - Brasil<br>CEP 09210-971<br>+55.11.4996-0196<br><a href="http://fileti.ufabc.edu.br">http://fileti.ufabc.edu.br</a><br>
<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
From: Mark Abraham <<a href="mailto:mark.abraham@anu.edu.au">mark.abraham@anu.edu.au</a>><br>
Subject: Re: [gmx-users] Dipeptide generation problem [Justin]<br>
To: Discussion list for GROMACS users <<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>><br>
Message-ID: <<a href="mailto:fb1e9eeeeb59.4c772156@anu.edu.au">fb1e9eeeeb59.4c772156@anu.edu.au</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
<br>
<br>
----- Original Message -----<br>
From: Eudes Fileti <<a href="mailto:fileti@ufabc.edu.br">fileti@ufabc.edu.br</a>><br>
Date: Friday, August 27, 2010 2:15<br>
Subject: Re: [gmx-users] Dipeptide generation problem [Justin]<br>
To: <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<br>
> Olá Justin, > thank you for responding to my post. I had tried what you mentioned before. > All I get is that the some atoms are missing (see error message in the > link <a href="https://sites.google.com/site/fileti/" target="_blank">https://sites.google.com/site/fileti/</a> ). > In fact, as I just want NH3+ (from LYS) and COO- (from ASP), the other atoms are not included > in PDB file and this can give error. If I put everything into the same group ("1" for example) many > of the atoms are deleted by pdb2gmx because they are duplicates.><br>
> In fact I still do not quite understand the logic of creating a PDB from the aminoacid residues. > I've read some tips on the list but still I could not understand.<br>
<br>
Have a look some MD tutorial material (doesn't have to be GROMACS). They probably start with a well-formed PDB file. That's probably a better learning experience than trying to read the format description.<br>
<br>
Alternatively, use a molecule builder (search webpage for suggestions here) to create the coordinates, and then get pdb2gmx to do the rest of the work.<br>
<br>
The idea behind pdb2gmx is that you don't have to move mountains to get your termini right, and the corresponding topology. If you just give it two adjacent PHE (charged or not), you can tell pdb2gmx what termini you want it to have. Done right, this whole task should not require you to inspect or change a coordinate file at all.<br>
<br>
Mark<br>
<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://lists.gromacs.org/pipermail/gmx-users/attachments/20100827/543767bc/attachment-0001.html" target="_blank">http://lists.gromacs.org/pipermail/gmx-users/attachments/20100827/543767bc/attachment-0001.html</a><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Thu, 26 Aug 2010 16:36:23 +0000<br>
From: Osmair Oliveira <<a href="mailto:osmair07@hotmail.com">osmair07@hotmail.com</a>><br>
Subject: RE: [gmx-users] Dipeptide generation problem [Justin]<br>
To: <<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>><br>
Message-ID: <SNT110-W61A3EA483256A7CAF02DECC3850@phx.gbl><br>
Content-Type: text/plain; charset="windows-1252"<br>
<br>
<br>
Hi Eudes,<br>
I do not know how work CHARMM force field, but if you use the following PDB file<br>
for your dipeptide using OPLS-AA force field, you can obtain good results!.<br>
<br>
<br>
ATOM 1 N PHE 1 0.000 0.000 0.000 1.00 0.00<br>
ATOM 2 H1 PHE 1 -0.940 0.000 -0.330 1.00 0.00<br>
ATOM 3 H2 PHE 1 0.470 0.820 -0.330 1.00 0.00<br>
ATOM 4 H3 PHE 1 0.470 -0.820 -0.330 1.00 0.00<br>
ATOM 5 CA PHE 1 0.000 0.000 1.460 1.00 0.00<br>
ATOM 6 HA PHE 1 -0.500 0.920 1.790 1.00 0.00<br>
ATOM 7 CB PHE 1 -0.780 -1.200 1.990 1.00 0.00<br>
ATOM 8 HB1 PHE 1 -0.340 -2.120 1.570 1.00 0.00<br>
ATOM 9 HB2 PHE 1 -0.670 -1.250 3.080 1.00 0.00<br>
ATOM 10 CG PHE 1 -2.240 -1.170 1.650 1.00 0.00<br>
ATOM 11 CD1 PHE 1 -3.140 -0.480 2.450 1.00 0.00<br>
ATOM 12 HD1 PHE 1 -2.770 0.040 3.340 1.00 0.00<br>
ATOM 13 CD2 PHE 1 -2.730 -1.830 0.530 1.00 0.00<br>
ATOM 14 HD2 PHE 1 -2.030 -2.380 -0.120 1.00 0.00<br>
ATOM 15 CE1 PHE 1 -4.490 -0.450 2.140 1.00 0.00<br>
ATOM 16 HE1 PHE 1 -5.180 0.100 2.790 1.00 0.00<br>
ATOM 17 CE2 PHE 1 -4.070 -1.800 0.220 1.00 0.00<br>
ATOM 18 HE2 PHE 1 -4.440 -2.330 -0.670 1.00 0.00<br>
ATOM 19 CZ PHE 1 -4.950 -1.110 1.020 1.00 0.00<br>
ATOM 20 HZ PHE 1 -6.020 -1.090 0.780 1.00 0.00<br>
ATOM 21 C PHE 1 1.400 0.000 2.020 1.00 0.00<br>
ATOM 22 O PHE 1 2.160 -0.950 1.880 1.00 0.00<br>
ATOM 23 N PHE 2 1.710 1.140 2.660 1.00 0.00<br>
ATOM 24 H PHE 2 0.990 1.860 2.700 1.00 0.00<br>
ATOM 25 CA PHE 2 3.020 1.330 3.260 1.00 0.00<br>
ATOM 26 HA PHE 2 3.770 1.260 2.460 1.00 0.00<br>
ATOM 27 CB PHE 2 3.100 2.720 3.900 1.00 0.00<br>
ATOM 28 HB1 PHE 2 2.290 2.820 4.630 1.00 0.00<br>
ATOM 29 HB2 PHE 2 4.050 2.800 4.460 1.00 0.00<br>
ATOM 30 CG PHE 2 3.020 3.850 2.920 1.00 0.00<br>
ATOM 31 CD1 PHE 2 4.150 4.290 2.250 1.00 0.00<br>
ATOM 32 HD1 PHE 2 5.120 3.820 2.450 1.00 0.00<br>
ATOM 33 CD2 PHE 2 1.810 4.460 2.650 1.00 0.00<br>
ATOM 34 HD2 PHE 2 0.900 4.120 3.160 1.00 0.00<br>
ATOM 35 CE1 PHE 2 4.070 5.330 1.340 1.00 0.00<br>
ATOM 36 HE1 PHE 2 4.980 5.670 0.830 1.00 0.00<br>
ATOM 37 CE2 PHE 2 1.730 5.500 1.740 1.00 0.00<br>
ATOM 38 HE2 PHE 2 0.760 5.980 1.540 1.00 0.00<br>
ATOM 39 CZ PHE 2 2.860 5.940 1.090 1.00 0.00<br>
ATOM 40 HZ PHE 2 2.800 6.760 0.360 1.00 0.00<br>
ATOM 41 C PHE 2 3.320 0.260 4.290 1.00 0.00<br>
ATOM 42 O1 PHE 2 2.370 -0.690 4.490 1.00 0.00<br>
ATOM 43 O2 PHE 2 4.520 0.340 4.920 1.00 0.00<br>
<br>
By<br>
Osmair<br>
Federal University of Sao Carlos - Brazil<br>
<br>
Date: Thu, 26 Aug 2010 13:14:20 -0300<br>
Subject: Re: [gmx-users] Dipeptide generation problem [Justin]<br>
From: <a href="mailto:fileti@ufabc.edu.br">fileti@ufabc.edu.br</a><br>
To: <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<br>
Olá Justin, thank you for responding to my post. I had tried what you mentioned before. All I get is that the some atoms are missing (see error message in the link <a href="https://sites.google.com/site/fileti/" target="_blank">https://sites.google.com/site/fileti/</a> ).<br>
In fact, as I just want NH3+ (from LYS) and COO- (from ASP), the other atoms are not included in PDB file and this can give error. If I put everything into the same group ("1" for example) many<br>
of the atoms are deleted by pdb2gmx because they are duplicates.<br>
In fact I still do not quite understand the logic of creating a PDB from the aminoacid residues. I've read some tips on the list but still I could not understand.<br>
<br>
That's it. If you have any suggestion, it will be very useful.<br>
Até maiseef_______________________________________<br>
Eudes Eterno Fileti<br>
Centro de Ciências Naturais e Humanas<br>
<br>
Universidade Federal do ABC — CCNH<br>
Av. dos Estados, 5001<br>
Santo André - SP - Brasil<br>
CEP 09210-971<br>
+55.11.4996-0196<br>
<a href="http://fileti.ufabc.edu.br" target="_blank">http://fileti.ufabc.edu.br</a><br>
<br>
<br>
<br>
<br>
------------------------------<br>
<br>
<br>
<br>
Message: 2<br>
<br>
Date: Thu, 26 Aug 2010 09:32:06 -0300<br>
<br>
From: Eudes Fileti <<a href="mailto:fileti@ufabc.edu.br">fileti@ufabc.edu.br</a>><br>
<br>
Subject: [gmx-users] Dipeptide generation problem<br>
<br>
To: <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<br>
Message-ID:<br>
<br>
<AANLkTikWxuC1g8fQ=fqnPfiXSWaDSTsPha=<a href="mailto:9QjPTDY38@mail.gmail.com">9QjPTDY38@mail.gmail.com</a>><br>
<br>
Content-Type: text/plain; charset="windows-1252"<br>
<br>
<br>
<br>
Hello everybody,<br>
<br>
I'm trying to create a dipeptide (L-PHE-L-PHE) from the PHE CHARMM residue.<br>
<br>
This dipeptide has a positively charged site (NH3+) and a negatively charged<br>
<br>
site (COO-).<br>
<br>
My PDB file ( <a href="https://sites.google.com/site/fileti/" target="_blank">https://sites.google.com/site/fileti/</a> ) does not seem to be<br>
<br>
consistent<br>
<br>
and produces error when I execute pdb2gmx.<br>
<br>
<br>
<br>
The link above presents the PDB and the error message from pdb2gmx<br>
<br>
(which I have not found similar in the forum).<br>
<br>
<br>
<br>
Could someone give me a hand with that?<br>
<br>
<br>
<br>
Bests<br>
<br>
eef<br>
<br>
_______________________________________<br>
<br>
Eudes Eterno Fileti<br>
<br>
Centro de Ciências Naturais e Humanas<br>
<br>
Universidade Federal do ABC — CCNH<br>
<br>
Av. dos Estados, 5001<br>
<br>
Santo André - SP - Brasil<br>
<br>
CEP 09210-971<br>
<br>
+55.11.4996-0196<br>
<br>
<a href="http://fileti.ufabc.edu.br" target="_blank">http://fileti.ufabc.edu.br</a><br>
<br>
-------------- next part --------------<br>
<br>
An HTML attachment was scrubbed...<br>
<br>
URL: <a href="http://lists.gromacs.org/pipermail/gmx-users/attachments/20100826/0cfce19f/attachment-0001.html" target="_blank">http://lists.gromacs.org/pipermail/gmx-users/attachments/20100826/0cfce19f/attachment-0001.html</a><br>
<br>
<br>
<br>
<br>
------------------------------<br>
<br>
<br>
<br>
Message: 3<br>
<br>
Date: Thu, 26 Aug 2010 08:36:24 -0400<br>
<br>
From: "Justin A. Lemkul" <<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>><br>
<br>
Subject: Re: [gmx-users] Dipeptide generation problem<br>
<br>
To: Discussion list for GROMACS users <<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>><br>
<br>
Message-ID: <<a href="mailto:4C765FC8.4020306@vt.edu">4C765FC8.4020306@vt.edu</a>><br>
<br>
Content-Type: text/plain; charset=windows-1252; format=flowed<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
Eudes Fileti wrote:<br>
<br>
> Hello everybody,<br>
<br>
> I'm trying to create a dipeptide (L-PHE-L-PHE) from the PHE CHARMM residue.<br>
<br>
> This dipeptide has a positively charged site (NH3+) and a negatively<br>
<br>
> charged site (COO-).<br>
<br>
> My PDB file ( <a href="https://sites.google.com/site/fileti/" target="_blank">https://sites.google.com/site/fileti/</a> ) does not seem to<br>
<br>
> be consistent<br>
<br>
> and produces error when I execute pdb2gmx.<br>
<br>
><br>
<br>
> The link above presents the PDB and the error message from pdb2gmx<br>
<br>
> (which I have not found similar in the forum).<br>
<br>
><br>
<br>
> Could someone give me a hand with that?<br>
<br>
><br>
<br>
<br>
<br>
Your .pdb file contains various broken residues (ASP, LYS) and non-sequential<br>
<br>
numbering (i.e., those broken residues are all numbered 1). Clean up the .pdb<br>
<br>
and try again.<br>
<br>
<br>
<br>
-Justin<br>
<br>
<br>
<br>
> Bests<br>
<br>
> eef<br>
<br>
> _______________________________________<br>
<br>
> Eudes Eterno Fileti<br>
<br>
> Centro de Ciências Naturais e Humanas<br>
<br>
> Universidade Federal do ABC — CCNH<br>
<br>
> Av. dos Estados, 5001<br>
<br>
> Santo André - SP - Brasil<br>
<br>
> CEP 09210-971<br>
<br>
> +55.11.4996-0196<br>
<br>
> <a href="http://fileti.ufabc.edu.br" target="_blank">http://fileti.ufabc.edu.br</a><br>
<br>
><br>
<br>
<br>
<br>
--<br>
<br>
========================================<br>
<br>
<br>
<br>
Justin A. Lemkul<br>
<br>
Ph.D. Candidate<br>
<br>
ICTAS Doctoral Scholar<br>
<br>
MILES-IGERT Trainee<br>
<br>
Department of Biochemistry<br>
<br>
Virginia Tech<br>
<br>
Blacksburg, VA<br>
<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | (540) 231-9080<br>
<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>
========================================<br>
<br>
<br>
<br>
--</blockquote></div></div></div>