On 7/10/17 2:32 AM, 조영래 wrote:
Thank you Justin
Every time your comment helped me.

I should compare two system that are Zn_binding protein and Fe_binding
protein.

So I need force field that satisfies the following conditions.
First, force field contains two metal ion (Fe and Zn)information.
Second, force field contain metal ion binding mode amino acid information
for example cysteine (CYM) in AMBER.

However many force fields did not satisfied above conditions.
So I would like to generate *.gro *itp file for metal atoms by GAMESS.
But according to your comment this is not correct.

In my case, the metal ions (Fe and Zn) binds one cysteine and two
histidines.
I want to simulate these complexes (i.e. Fe_protein and Zn_protein complex).
Please suggest me an ideal force field that in able to be used for
sumulation of the avode mentaioned complexes.


I don't know of one. Dealing with transition metals is one of the hardest tasks in parametrization. There's a reason why after many decades of development, the parameters still don't exist... Depending on what you need to do, hybrid QM/MM is probably the best/most straightforward approach. Classical force fields do not deal well with these kinds of systems.

-Justin

--
==================================================

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 629
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==================================================
--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.

Reply via email to