Ferric uptake regulator protein: Binding free energy calculations and per-residue free energy decomposition

Iron homeostasis is, in many bacterial species, mediated by the ferric uptake regulator (Fur). A regulatory site able to bind iron to activate Fur for DNA binding has been described, and a structural zinc site essential for the dimerization has also been proposed. They have been localized and named site I and site 2, respectively, from the crystal structure of a zinc-substituted Pseudomonas aeruginosa Fur (PA-Fur). Notwithstanding the studies on Fur proteins from various species, both the precise site of iron binding and the effect on DNA binding affinity are still controversial. These issues were investigated here by molecular dynamics simulations and free energy calculations. Simulations were performed for eight molecular systems represented by the three forms of Fur, that is, apo Fur, metal-substituted Fur, and Fur complexed with DNA. Because of the lack of a Fur-DNA complex crystal structure, the recently published model based on mass spectrometry experiments on Escherichia coli Fur (EC-Fur), and the crystal structure of PA-Fur, was used, after adjustment to adopt a symmetric conformation. The simulation results suggest that the formerly proposed site 2 is, in fact, the regulatory iron-sensing site. The calculations also predict that Fe2+ at site 2 is hexacoordinated having an octahedral environment with only nitrogen and oxygen atoms, which is in accordance with previous spectroscopic characterizations. Energy decomposition pinpoints H87 as an additional amino acid that defines the regulatory metal site. Finally, free energy decomposition analysis reveals a number of amino acids potentially important in dimerization and in DNA binding.