Dear Amber users,
I’m planning to run MM-PBSA calculations of the binding affinity between a small molecule ligand and a complex formed by a protein, DNA and a Mg2+ ion with mmpbsa.py. All the elements of the complex interact with the small molecule and therefore contribute to the its binding affinity.
When trying to include the metal ion as part of the receptor for the mmpbsa.py calculation, I read in the Amber manual that “If there is a specific metal ion that you wish to include in the calculation, you can prevent ptraj from stripping this ion by NOT specifying it in strip_mask”.
However, the MMPBSA.py publication by Miller et al., 2012, mentions that “the implicit solvent models utilized by MMPBSA.py are not parametrized to account for these species. Therefore, it is not recommended to include water or ions during the calculations”.
I’d like to know if you have previous experience in computing the binding affinity of ligands to metalloprotein + DNA complexes, and if so if you need to add anything into the calculation to adequately parameterize the Mg2+ ion so it is taken into account when calculating the binding affinity. Many thanks for your help.
Best wishes,
Diego
---------------------------------------------
Diego Alonso-Martinez
Biomolecular Mass Spectrometry and Computational Proteomics Group (BMSCPG)
Department of Chemical Engineering
Imperial College London (UK)
London SW7 2AZ
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Oct 27 2017 - 06:30:02 PDT
