I am going to simulate a tRNA molecule using GROMOS96 force field, which contain some modified nucleotide.
I've build the .rtp entry for each modified residues like the this: A02-mdf.top <http://gromacs.5086.n6.nabble.com/file/n5001965/A02-mdf.top> But pdb2gmx prompt out such error message: Fatal error: atom N not found in buiding block 1GUA while combining tdb and rtp I choose gromos43a1 force field. I think pdb2gmx treat my nucleotide as protein, because it cap both terminus by NH3+ and COO-. So I redo it again, with -ter option, and choose none. pdb2gmx say: Fatal error: There is a dangling bond at at least one of the terminal ends. Select a proper terminal entry. Looks like that I should edit the terminal entry, is it correct? My .rtp file almost exactly the same as proper nucleotide topology, the differences exist only in the residue part. Here is another question: why the total charge of single nucleotide molecule is -1? Sincerely Daniel LIn -- View this message in context: http://gromacs.5086.n6.nabble.com/Nucleotide-terminal-problems-in-GORMOS96-force-field-tp5001965.html Sent from the GROMACS Users Forum mailing list archive at Nabble.com. -- gmx-users mailing list gmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users * Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! * Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
