Parameters for molecular dynamics simulations of iron‐sulfur proteins

Abstract

Iron-sulfur proteins involved in electron transfer reactions have finely tuned redox potentials, which allow them to be highly efficient and specific. Factors such as metal center solvent exposure, interaction with charged residues, or hydrogen bonds between the ligand residues and amide backbone groups have all been pointed out to cause such specific redox potentials. Here, we derived parameters compatible with the AMBER force field for the metal centers of iron-sulfur proteins and applied them in the molecular dynamics simulations of three iron-sulfur proteins. We used density-functional theory (DFT) calculations and Seminario's method for the parameterization. Parameter validation was obtained by matching structures and normal frequencies at the quantum mechanics and molecular mechanics levels of theory. Having guaranteed a correct representation of the protein coordination spheres, the amide H-bonds and the water exposure to the ligands were analyzed. Our results for the pattern of interactions with the metal centers are consistent to those obtained by nuclear magnetic resonance spectroscopy (NMR) experiments and DFT calculations, allowing the application of molecular dynamics to the study of those proteins.

Insights into the key interactions between human protein phosphatase 5 and cantharidin using molecular dynamics and site-directed mutagenesis bioassays

“…Atomic single charge was calculated using the restrained electrostatic potential (RESP) fitting protocol 31 . The Quantum Mechanics (QM) and Molecular Mechanics (MM) vibrational frequencies were compared to check the quality of the force field parameters 30 31 32 33 . We performed three individual 50 ns MD at different seeds for each system.…”

Identification of Crucial Amino Acids in Mouse Aldehyde Oxidase 3 That Determine Substrate Specificity

“…Dihedral force constants involving the Mo were set to zero, while transferable van der Waals atomic parameters were taken from the literature [48]. The parameters for the Fe 2 S 2 centers were taken from literature too [49].…”

New insights about the monomer and homodimer structures of the human AOX1

“…14 The parameters of the Fe 2 S 2 clusters (FES) were taken from the literature. 15 The THI molecule was parameterized in its neutral form, even though its pK a in aqueous solution is reported to be 9.5, 16 because, as it can be observed in the X-ray structure, it sits in a hydrophobic pocket. Additionally, it can be seen that several hydrogen bond acceptors are present around the nitrogen atom of the methyl piperidine ring and free to reorient.…”