[gmx-users] part2-adding metal-bonding parameters to a force field
Ahmet yıldırım
ahmedo047 at gmail.com
Tue Apr 15 12:21:02 CEST 2014
Hi,
There are two H1D residue in modified pdb file and one H1E residue. ZN atom
in H1E are bound to OW and ND1 atoms in H1E and NE2 atom in H1D. That is,
Zinc has four bond.
NE2(H1D)
ı
(H1E)ND1-ZN-OW(H1E)
ı
NE2(H1D)
specbond.dat:
Line 1:H1E ND1 1 H1E ZN 1 0.2 H1E H1E
Line 2:H1E OW 1 H1E ZN 1 0.2 H1E H1E
Line 3:H1D NE2 1 H1E ZN 1 0.2 H1D H1E
After pdb2gmx:
Linking H1D-94 NE2-754 and H1E-119 ZN-920...
Linking H1D-96 NE2-775 and H1E-119 ZN-920...
Linking H1E-119 ND1-952 and H1E-119 ZN-920...
Linking H1E-119 OW-921 and H1E-119 ZN-920..
I know the above linkings are correct.
I wonder, should I change as the following the line 3 in specbond.dat
because there are two H1D residue in structure?
H1D NE2 1 H1E ZN 2 0.2 H1D H1E
2014-04-06 11:48 GMT+03:00 Tsjerk Wassenaar <tsjerkw at gmail.com>:
> Hi Ahmet,
>
> Nothing wrong with specbond at this point. The bonds are correctly added,
> and apparently the bondtypes are read correctly too, because grompp only
> starts complaining at the point of angles. Check the atomtypes
> corresponding to the atoms involved in those angles in the .itp file that
> is generated.
>
> By the way, the message from moderator is due to the constant replying in
> the thread, including the previous mails. It's getting pretty lengthy like
> that. Sometimes it's good to keep only the parts of earlier messages that
> are directly being replied too.
>
> Cheers,
>
> Tsjerk
> --
> Gromacs Users mailing list
>
> * Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before
> posting!
>
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
>
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-request at gromacs.org.
>
--
Ahmet Yıldırım
More information about the gromacs.org_gmx-users
mailing list