Molecular dynamics simulations were performed on both apo and copper forms of the human copper chaperone, Hah1. Wild-type Hah1 and a methionine (M10) to serine mutant were investigated. We have evidenced the central role of residue M10 in stabilizing the hydrophobic core of Hah1 as well as the internal structure of the metal-binding site. When copper(I) is bound, the mobility of Hah1 is reduced whereas mutation of M10 implies a drastic increase of the mobility of apoHah1, stressing the importance of this highly conserved hydrophobic residue for copper sequestration by the apoprotein.

We have developed new force field and parameters for copper(I) and mercury(II) to be used in molecular dynamics simulations of metalloproteins. Parameters have been derived from fitting of ab initio interaction potentials calculated at the MP2 level of theory, and results compared to experimental data when available. Nonbonded parameters for the metals have been calculated from ab initio interaction potentials with TIP3P water. Due to high charge transfer between Cu(I) or Hg(II) and their ligands, the model is restricted to a linear coordination of the metal bonded to two sulfur atoms. The experimentally observed asymmetric distribution of metal ligand bond lengths (r) is accounted for by the addition of an anharmonic (r3) term in the potential. Finally, the new parameters and potential, introduced into the CHARMM force field, are tested in short molecular dynamics simulations of two metal thiolates fragments in water. (Brooks BR et al. J Comput Chem 1983, 4, 1987.1).