[gmx-users] Topology for hetroatom (PRODRG)
Alok
alokjain at iitk.ac.in
Thu Feb 24 19:42:24 CET 2005
hello david,
I am attaching the rtp file for my LINKER
(LIN)along with this mail.
H O
| ||
-CA-C(-N-CA-CB-N1)-C-CA-
|| |
O H1
The linker section (-NH-CH2-CH2-NH-) attached to the carbonly group
of two adjacent amino acids is as shown above.
Here CA and CB within my linker (shown in bold) atoms are of atom
type CH2.
The rtp file for this linker was developed by checking the
corresponding "atom types" from the database topology file of gromacs.
rtp file:
[ LIN ]
[ atoms ]
N N -0.28000 0
H H 0.28000 0
CA CH2 0.00000 1
CB CH2 0.00000 2
N1 N -0.28000 3
H1 H 0.28000 3
[ bonds ]
N H gb_2
N CA gb_20
CA CB gb_26
CB N1 gb_20
N1 H1 gb_2
N1 +C gb_9
[ angles ]
; ai aj ak gromos type
-C N H ga_31
-C N CA ga_30
H N CA ga_17
N CA CB ga_12
CA CB N1 ga_12
CB N1 H1 ga_17
CB N1 +C ga_30
H1 N1 +C ga_31
[ impropers ]
; ai aj ak al gromos type
N -C CA H gi_1
N1 CB +C H1 gi_1
[ dihedrals ]
; ai aj ak al gromos type
-CA -C N CA gd_4
-C N CA CB gd_19
N CA CB N1 gd_17
CA CB N1 H1 gd_19
CA CB N1 +C gd_19
But as i have mentioned before,pdb2gmx program does not read my LIN
molecule.
Please suggest something.
Thanking you,
Alok Jain
David wrote:
>On Wed, 2005-02-23 at 22:26 +0530, Alok wrote:
>
>
>>Hello David,
>> I include my linker molecule (defined as LIN )in
>>the rtf and hdb files. i also add this LIN in aminoacids.dat file
>>because i want to consider this linker as a part of protein. I also
>>increase the total count in aminoacids.dat file by one.
>> but pdb2gmx program does not reconize my linker lolecule (LIN)
>>i encounter follwing fatal error..
>>
>>
>
>rtp file.
>
>
>
>>Fatal error: Residue 'LIN' not found in residue topology database
>>
>>i also check the pdb2gmx.log file (flag -debug) in which it is not
>>reading my linker molecule.
>>
>>part of the ouput from the log file generated by pdb2gmx which shows the
>>molecules read by the program from the rtf file is as :
>>
>> 80 AR 1 0 0 0 0 0
>> 81 SO42- 5 4 6 0 0 0
>> 82 ZN2+ 1 0 0 0 0 0
>> 83 NA+ 1 0 0 0 0 0
>> 84 CL- 1 0 0 0 0 0
>> 85 CA2+ 1 0 0 0 0 0
>> 86 MG2+ 1 0 0 0 0 0
>> 87 CU1+ 1 0 0 0 0 0
>> 88 CU2+ 1 0 0 0 0 0
>> 89 H2O 3 3 0 0 0 0
>> 90 HOH 3 3 0 0 0 0
>> 91 HO4 4 3 0 0 0 0
>> 92 H2OE 3 3 0 0 0 0
>> 93 CHCL3 5 9 0 0 0 0
>> 94 DMSO 4 6 0 0 0 0
>> 95 CH3OH 3 3 0 0 0 0
>> 96 CCL4 5 9 0 0 0 0
>>
>>Here in this output i do not find my LIN molecule which is defined after
>>CCL4 in the rtf file.
>>
>>Am i going wrong somewhere?? Please suggest.
>>
>> Thanking you,
>>
>> Alok Jain
>>
>>
>>David wrote:
>>
>>
>>
>>>On Wed, 2005-02-23 at 11:52 +0530, Alok wrote:
>>>
>>>
>>>
>>>
>>>>hello David,
>>>> Thank you for your reply.So we will include our
>>>>linker section molecule in the rtf and hdb files.But we are still not
>>>>sure of the parameters we are to take for the linker molecule.Should I
>>>>directly take the parameters mentioned in the ' itp file ' generated by
>>>>the PRODRG server ??????
>>>> The linker section forms a covalent bond with
>>>>the peptide.The carbonyl group at both ends of the peptide
>>>>(XXX-NH-CH2-CH2-NH-XXX) forms a covalent bond with the amino group at
>>>>the ends of the linker section.
>>>>
>>>>
>>>>
>>>>
>>>>
>>>It depends on your force field.
>>>
>>>
>>>
>>>
>>>
>>_______________________________________________
>>gmx-users mailing list
>>gmx-users at gromacs.org
>>http://www.gromacs.org/mailman/listinfo/gmx-users
>>Please don't post (un)subscribe requests to the list. Use the
>>www interface or send it to gmx-users-request at gromacs.org.
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20050225/9a847aeb/attachment.html>
More information about the gromacs.org_gmx-users
mailing list