Probing Polyoxometalate-Protein Interactions Using Molecular Dynamics Simulations.

The molecular interactions between the Ce -substituted Keggin anion [PW O Ce(OH ) ] (CeK) and hen egg-white lysozyme (HEWL) were investigated by molecular dynamics simulations. The analysis of CeK was compared with the Ce -substituted Keggin dimer [(PW O ) Ce] (CeK ) and the Zr -substituted Lindqvist anion [W O Zr(OH )(OH)] (ZrL) to understand how POM features such as shape, size, charge, or type of incorporated metal ion influence the POM⋅⋅⋅protein interactions. Simulations revealed two regions of the protein in which the CeK anion interacts strongly: cationic sites formed by Arg21 and by Arg45 and Arg68. The POMs chiefly interact with the side chains of the positively charged (arginines, lysines) and the polar uncharged residues (tyrosines, serines, aspargines) via electrostatic attraction and hydrogen bonding with the oxygen atoms of the POM framework. The CeK anion shows higher protein affinity than the CeK and ZrL anions, because it is less hydrophilic and it has the right size and shape for establishing interactions with several residues simultaneously. The larger, more negatively charged CeK anion has a high solvent-accessible surface, which is sub-optimal for the interaction, while the smaller ZrL anion is highly hydrophilic and cannot efficiently interact with several residues simultaneously.