Hi Ana,
I did not have this experience, but I know a similar example mentioned in the following reference "Lin, F.; Wang, R., Systematic derivation of AMBER force field parameters applicable to zinc-containing systems. J. Chem. Theory Comput. 2010, 6, 1852-1870.” Where the PM3/MM MD simulation gives a notable distortion of a zinc site.
If you just want to simulate the vibration of the metal site (no chemical reactions), pure MM model is enough, you can try to use the zinc AMBER force field for your case: http://ambermd.org/tutorials/advanced/tutorial20/ZAFF.htm <http://ambermd.org/tutorials/advanced/tutorial20/ZAFF.htm>, or use the MCPB.py to parameterize the metal site: http://ambermd.org/tutorials/advanced/tutorial20/mcpbpy.htm <http://ambermd.org/tutorials/advanced/tutorial20/mcpbpy.htm>, http://ambermd.org/tutorials/advanced/tutorial20/mcpbpy_heme.htm <http://ambermd.org/tutorials/advanced/tutorial20/mcpbpy_heme.htm>.
Hope it helps,
Pengfei
> On Mar 21, 2019, at 12:18 PM, Dreab, Ana <adreab.odu.edu> wrote:
>
> Hi all,
>
> I am trying to run DFTB3/MM simulation (using Amber18) of CCCC zinc finger and one of the cysteine residues keeps dissociating from the Zn(II). Has anyone encounter a similar problem? Any help will be much appreciated
>
> The qm region includes the Zn(II) and SCH2 fragment of each of the four cysteine residues. The overall charge of the qm region is -2.
>
> Best Regards,
> Ana Dreab
>
> Department of Chemistry and Biochemistry
> Old Dominion University
> Norfolk, VA 23529
>
>
> _______________________________________________
> 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 Sat Mar 23 2019 - 13:00:03 PDT
