Classical force field parameters for the heme prosthetic group of cytochrome c.
Accurate force fields are essential for describing biological systems in a molecular dynamics simulation. To analyze the docking of the small redox protein cytochrome c (cyt c) requires simulation parameters for the heme in both the reduced and oxidized states. This work presents parameters for the partial charges and geometries for the heme in both redox states with ligands appropriate to cyt c. The parameters are based on both protein X-ray structures and ab initio density functional theory (DFT) geometry optimizations at the B3LYP/6-31G* level. The simulations with the new parameter set reproduce the geometries of the X-ray structures and the interaction energies between water and heme prosthetic group obtained from B3LYP/6-31G* calculations. The parameter set developed here will provide new insights into docking processes of heme containing redox proteins.
Autenrieth
Felix
F
School of Chemical Sciences, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA.
Tajkhorshid
Emad
E
Baudry
Jerome
J
Luthey-Schulten
Zaida
Z
D000465:Algorithms; D003198:Computer Simulation; D018360:Crystallography, X-Ray; D045304:Cytochromes c / Q000737:chemistry*; D006418:Heme / Q000737:chemistry*; D008958:Models, Molecular; D008968:Molecular Conformation; D015394:Molecular Structure; D010084:Oxidation-Reduction; D055672:Static Electricity; D013816:Thermodynamics
D006418:Heme; D045304:Cytochromes c