> > set UNIT.1.3 element C
I wonder if putting C in quotes: "C" would work.
Bill
"Aronica, Pietro" <pietro.aronica07.imperial.ac.uk> wrote:
> Hello everyone,
> I found a possible bug in LEaP when used with ff12SB. I haven't seen this described anywhere else, and it could be of interest to other people.
> My work involves a fair amount of custom residues, so I often have to generate new ones and give parameters to them with LEaP. The way I go about doing this is by using the set command, this way:
> > set XXX.N.M parameter object
> where XXX is the unit, N the residue, M the atom number, parameter what I need to set (element, atom type, atom name, and so on) and object what I want the parameter to be. The problem arose when I tried to add a residue after upgrading from ff99SB and setting the parameters to the correct values: if I tried to set the element to something as simple as carbon (C), with a command like this,
> > set UNIT.1.3 element C
> it would not allow me to do so, though nitrogen (N), hydrogen (H), and oxygen (O) all worked without problem; similarly, it wouldn't let me set the type of the carbonyl carbon on the backbone to its correct type (C), while all other types worked. The error is always the same: "The value must be of the type: String". This did not appear with ff99SB.
> After some prodding and poking, I figured that the problem happens because the string C is assigned to a premade residue: cytosine. LEaP misinterprets the command and fails to change the element or type to C. This occurs with any prebuilt residue name, from A to WAT, and did not happen in ff99SB because cytosine is named RC, perhaps specifically to avoid this issue.
> As I need to set things to C very often, my fix will probably be to rename cytosine back to RC in some way, but I'm open to other workarounds.
> I hope this problem can get solved.
> Kind regards
> Pietro Aronica
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Oct 26 2012 - 00:00:02 PDT
