[AMBER] error in leap -- Glycam residue building

From: lara lara <lara.4884.gmail.com>
Date: Tue, 10 Sep 2013 15:29:47 -0400

Dear AMber !

I have an leap error when i generated the glycam residues

I would like to join two sugar units with 1O of first to the second that is
done by the letter code
but I am not able to add sulfate on the 1O free sugars.

like I want to add sulfate in position 2 of first residue.

I used the command as below , How to fix this .

> m = sequence { 1uA}
> edit m
> set m tail m.1.O1
> m = sequence { m 2uA}
> edit m
> remove m m.1.H2O
> set m.1.O2 charge -0.4630
> set m tail m.1.O2
> m = sequence { m SUL}
> edit m
> set m.2.O2 charge -0.4630
> set m tail m.2.O2
> m = sequence { m SUL}
> edit m
> saveamberparm m l.prmtop l.inpcrd
Checking Unit.
WARNING: The unperturbed charge of the unit: -4.000000 is not zero.

 -- ignoring the warning.

Building topology.
Building atom parameters.
Building bond parameters.
Could not find bond parameter for: Oh - S
Building angle parameters.
Could not find angle parameter: Cg - Oh - S
Could not find angle parameter: Oh - S - O2
Could not find angle parameter: Oh - S - O2
Could not find angle parameter: Oh - S - O2
Building proper torsion parameters.
 ** No torsion terms for Cg-Cg-Oh-S
 ** No torsion terms for Cg-Cg-Oh-S
 ** No torsion terms for Cg-Oh-S-O2
 ** No torsion terms for Cg-Oh-S-O2
 ** No torsion terms for Cg-Oh-S-O2
 ** No torsion terms for H1-Cg-Oh-S
Building improper torsion parameters.
old PREP-specified impropers:
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Sep 10 2013 - 13:00:02 PDT
Custom Search