[gmx-users] DNA virtual site with Charmm force field

[Erik Marklund]

Hi Terry,

Correct. The forcefields provided with gromacs do not have a complete set of parameters for making vsites for nucleotides. You can work them out from geometry and mass moments however.

Kind regards,
Erik

> On 1 Jul 2015, at 10:23, Terry <terrencesun at gmail.com> wrote:
> 
> Hi all,
> 
> I'm trying to use virtual sites in a DNA in water simulation with Charmm27
> force field. When feeding my pdb to pdb2gmx with this cammand "pdb2gmx -f
> my_DNA.pdb -vsite hydrogens', it gives me this error
> 
> Fatal error:
> No vsite database NH2 entry for type NN1
> 
> Is this not implemented or did I miss something? Any input would be
> appreciated.
> 
> Regards
> Terry
> -- 
> 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.