Hi everybody
Thank you very much Tsjerk. I'm sorry for asking all these questions
but some points are still unclear to me, and I didn't find much
information on the mailing list archives or in the manual that could
help me.
Grompp is complaining about the angles because in the .top file there
are no entries for the angle type, because pdb2gmx just creates the
bond, but doesn't specify the type like you told me.
I've edited the ff...bon.itp file as follows but it seems useless.
; bond-, angle- and dihedraltypes for specbonds:
[ bondtypes ]
S S 2 gb_33
NR FE 2 gb_32
; cystine - heme link (is CR1-S, use CH2-S):
S CR1 2 gb_30
; cystine - fs link:
FE S 2 gb_48
[ angletypes ]
CH1 CH2 S 2 ga_15
CH2 S S 2 ga_5
CR1 NR FE 2 ga_33
NR FE NR 2 ga_16
; cystine - heme link (is CH2-S-CR1, use CH2-S-CH2):
CH2 S CR1 2 ga_3
; cystine - heme link (is S-CR1-C/CH2, use CHn-CH2-S):
S CR1 C 2 ga_15
S CR1 CH2 2 ga_15
; cystine - fs link
S FE S 2 ga_47
FE S FE 2 ga_48
etc...
where I specified ga_47, ga_48, gb_48...
Is there a way that grompp could take these types or do I have to edit
the topology file manually to specify the type?
It seems strange to me because there are also no entries for the bond
types, and grompp doesn't give me any warnings for that, why?
Thank you again
Nicolas
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
