AMBER Archive (2006)

Subject: AMBER: bad connectivity - more detail

From: Seth Lilavivat (sethl_at_gatech.edu)
Date: Tue Dec 26 2006 - 12:02:27 CST


Dear Amber Community,

I am a novice AMBER user, creating modified nucleic acid structures. The
current model that I am working on is a bicyclo RNA structure. Here is a brief
summary of my current protocol:

1. Start with closest residue, in this case RA5, RA, RA3, RU5, RU, and RU3
2. Modifiy them in xleap to have the correct structure.
3. Create a residue library file in one of the following two ways:
   1. Generate a PDB of the new residue, then use antechamber to parameterize
   2. Save the modified residue as a prep file from xleap. Then switch out all
of the partial charges by hand, check all of the atom types, and try and modify
necessary force field parameters.

One problem that I am running into is:

Connectivity - when use antechamber/gaff generated prep residues, I get a link
from C4 (if it is an A) to O1P or O2 to O1P (if it is a U). The link should be
from the O2* to the P (since this is a 5' - 2' structure).

Attached are some residue files and a pdb. I hope this helps.

Thanks,
SETH















-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu