Dear Amber users,
I want to perform MD simulations on a glycoprotein with an active site zinc (+2), water and chloride ion (-1).
The metal site (ZN, water, HD1, HD2, OE1) forcefield parameters were determined using MCPB.py. The hybridization is correctly specified in the frcmod file (HD1, HD2, OE1 sp2) but when visualizing the pdb (saved from the prmtop and rst7 of tleap) it doesn't correspond (HD1 sp3, HD2 sp2, OE1 sp3, ZN sp2, Water sp3). Can I assume that it is in fact correctly given as the aforementioned and that the conversion to pdb just doesn't capture this correctly? Or should it be set in the tleap.in file at AddAtomTypes by changing sp3?
addAtomTypes {
{ "X1" "N" "sp3" }
{ "X2" "N" "sp3" }
{ "X3" "O" "sp3" }
{ "Z1" "Zn" "sp3" }
{ "X4" "O" "sp3" }
}
Is there any way to inspect the atom hybridization in AMBER?
Since glycan addition removed the chloride ion, I copied its coordinates from the deposited pdb into the Amber-format glycoprotein pdb. The chloride atom type was set to Cl- as shown below
HETATM12492 H2O 0fA 720 -35.331 -14.332 -7.397 1.00 0.00 H
TER
HETATM12493 ZN ZN1 721 -27.766 -21.998 -36.628 1.00 0.00 Zn
TER
HETATM12494 Cl- Cl- 722 -18.015 -27.711 -36.880 1.00 0.00 Cl-
TER
ATOM 12495 O WT1 723 -27.107 -20.515 -34.545 1.00 0.00 O
ATOM 12496 H1 WT1 723 -26.150 -20.515 -34.545 1.00 0.00 H
ATOM 12497 H2 WT1 723 -27.347 -19.588 -34.545 1.00 0.00 H
TER
END
The atomic_ions.lib file lists Cl- so I assumed this was correct but when visualizing the tleap output as pdb (generated from prmtop and rst7) in UCSF Chimera and Discovery Studio, the atomtype is given as a carbon (sp3) with no formal charge. Is there a way to inspect the Cl- atom in tleap? (I've tried desc from the manual but new to AMBER and not sure how to use it)
What is the correct way to insert a chloride ion in AMBER format?
Kind regards
Lizelle Lubbe
PhD (Medical biochemistry) candidate
Department of Integrative Biomedical Sciences
University of Cape Town
Disclaimer - University of Cape Town This e-mail is subject to UCT policies and e-mail disclaimer published on our website at
http://www.uct.ac.za/about/policies/emaildisclaimer/ or obtainable from +27 21 650 9111. If this e-mail is not related to the business of UCT, it is sent by the sender in an individual capacity. Please report security incidents or abuse via csirt.uct.ac.za
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Jul 10 2017 - 05:30:03 PDT
