Abstract
Iron binding to protein is common in biological processes of dioxygen transport, electron transfer as well as in stabilizing drug-protein complexes. α-Helix is the most prevalent secondary structure of proteins. In this study, Fe(2+) binding to α-helix has been studied by isothermal titration calorimetry (ITC) and explicitly solvated molecular dynamics (MD) simulation. Ferrous gluconate and α-helix-rich keratin are used for the ITC study and the results revealed followed one set of identical sites binding model. The MD simulations further revealed that only the acidic side-chain functional groups and η(2) (O,O) coordination modes are involved in the binding of Fe(2+) to α-helix. The ITC results also showed that the binding of ferrous gluconate to keratin was entropy driven and the higher the temperature, the stronger the binding free energy. The favorable entropy of Fe(2+) binding to keratin was attributed to the displacement of water molecules on the α-helix surface, and was confirmed via MD simulations. The most stable coordination states of Fe(2+) and α-helix were identified via simulation: Fe(2+) stacks between two glutamic acid side chain carboxylate groups, displacing water molecules. The binding free energies calculated using MD simulation and the theoretical values were in excellent agreement with the ITC results.