Parametrization of a force field for metals complexed to biomacromolecules: applications to Fe(II), Cu(II) and Pb(II)

Although techniques for the simulation of biomolecules, such as proteins and RNAs, have greatly advanced in the last decade, modeling complexes of biomolecules with metal ions remains problematic. Precise calculations can be done with quantum mechanical methods but these are prohibitive for systems the size of macromolecules. More qualitative modeling can be done with molecular mechanical potentials but the parametrization of force fields for metals is often difficult, particularly if the bonding between the metal and the groups in its coordination shell has significant covalent character. In this paper we present a method for deriving bond and bond-angle parameters for metal complexes from experimental bond and bond-angle distributions obtained from the Cambridge Structural Database. In conjunction with this method, we also introduce a non-standard energy term of gaussian form that allows us to obtain a stable description of the coordination about a metal center during a simulation. The method was evaluated on Fe(II)-porphyrin complexes, on simple Cu(II) ion complexes and a number of complexes of the Pb(II) ion.