Re: [AMBER] tleap atom type and hybridization

From: Pengfei Li <ambermailpengfei.gmail.com>
Date: Mon, 10 Jul 2017 21:27:58 -0500

Hi Lizelle,

Based on my understanding, the atom types don’t really play a role in the molecular dynamics simulation in AMBER since all the bond linkage information is specified independently.

For the ion recognition in LEaP, this webpage may be helpful: http://ambermd.org/tutorials/advanced/tutorial20/nonbonded_model.htm <http://ambermd.org/tutorials/advanced/tutorial20/nonbonded_model.htm>.

Kind regards,
Pengfei

> On Jul 10, 2017, at 7:13 AM, Lizelle Lubbe <LBBLIZ002.myuct.ac.za> wrote:
>
>
> 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

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Jul 10 2017 - 19:30:05 PDT
Custom Search