Re: AMBER: library for 8oxoguanine

From: M. L. Dodson <>
Date: Mon, 24 Jul 2006 15:07:17 -0500

On Monday 24 July 2006 14:50, Carra, Claudio (JSC-SK)[USRA] wrote:
> Dear All,
> i've been scrambling to make the correct library for the 8oxoGuanine
> in a sequence of base pairs, with a reasonably good orientation.
> i've already written to this forum for help on that but i didn't
> get any reply, so i'm simplifying my issue.
> I've followed the example proposed in the tutorial 4,
> but i got a problem with missing parameter for the force field,
> in generally fixable with parmchk, but not in this case since
> the residue is already inserted in the DNA chain.
> however the tutorial at the very and point out that
> [...]"Creating a set of parameters is a bit of an art, since one is
> often faced with unknown parameters (such as force constants for
> copper to its ligants as in this example). The purpose of this
> tutorial is simply to cover the mechanics of running AMBER and
> so I provide you all of the parameters below. Note, these are for
> tutorial purposes only, I make no claims about the validity or
> appropriateness of these parameters.[...]"
> so essentially my question is whether there is a TESTED way to
> produce a DNA chain, double and single strand, with a 8oxoguanine
> residue? the question can be extended to the general problem of
> changing one or more atom (for instance a carbossilic substituent
> with an amino group) on a given protein, and to make sure that
> is it properly bonded and describe by the force field, keeping
> a reasonably orientation from the starting guess.
> i would really appreciate any comment, even negative.
> thanks a lot
> sincerely
> claudio

Maybe I am misunderstanding you, but you do not modify the
original DNA to produce a new DNA with an 8OG residue. Instead,
you generate a new residue, 8OG, complete with all the charges,
topology, etc similar to what the regular DNA residues have. You
may need to generate a frcmod file that contains additions to the
force field if your new residue contains bonds, angles or torsions
that are not present in the regular force field description.

I suggest you use Antechamber or REDII (see the archives for many
references to the location of REDII). The Antechamber web page
gives examples of generating residue descriptor files for new
residues. REDII automates the procedure, but is somewhat oriented
toward standalone chemical compounds. Since you are generating a
polymer internal residue, you will need to become pretty familiar
with the resp software to use REDII. It may be that the new
database section for the REDII web site has tutorials to do what
you want.

Once you have the new 8OG residue description, you can build your
DNA using it for MD.

M. L. Dodson
Personal:       mldodson-at-houston-dot-rr-dot-com
Work:           bdodson-at-scms-dot-utmb-dot-edu
Work: four_zero_nine-772-2178             FAX: four_zero_nine-747-8608
The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Wed Jul 26 2006 - 06:07:08 PDT
Custom Search