title: Performance of Protein-Ligand Force Fields for the Flavodoxin-Flavin Mononucleotide System.
abstract: The ability to accurately perform molecular dynamics and free energy perturbation calculations for protein-ligand systems is of broad interest to the biophysical and pharmaceutical sciences. In this work, several popular force fields are evaluated for reproducing experimental properties of the flavodoxin/flavin mononucleotide system. Calculated (3)J couplings from molecular dynamics simulations probing φ and χ1 dihedral angles are compared to over 1000 experimental measurements. Free energy perturbation calculations were also executed between different protein mutants for comparison with experimental data for relative free energies of binding. Newer versions of popular protein force fields reproduced (3)J backbone and side chain couplings with good accuracy, with RMSD values near or below one hertz in most cases. OPLS-AA/M paired with CM5 charges for the ligand performed particularly well, both for the (3)J couplings and FEP results, with a mean unsigned error for relative free energies of binding of 0.36 kcal/mol.
authors.affiliation: Department of Chemistry, Yale University , New Haven, Connecticut 06520-8107, United States.
authors.forename: Michael J
authors.initials: MJ
authors.lastname: Robertson
authors.affiliation: Department of Chemistry, Yale University , New Haven, Connecticut 06520-8107, United States.
authors.forename: Julian
authors.initials: J
authors.lastname: Tirado-Rives
authors.affiliation: Department of Chemistry, Yale University , New Haven, Connecticut 06520-8107, United States.
authors.forename: William L
authors.initials: WL
authors.lastname: Jorgensen
mesh_terms: D005486:Flavin Mononucleotide / Q000737:chemistry*; D005418:Flavodoxin / Q000737:chemistry*; D008024:Ligands; D008958:Models, Molecular; D011487:Protein Conformation; D011506:Proteins / Q000737:chemistry*
chemical_list: D005418:Flavodoxin; D008024:Ligands; D011506:Proteins; D005486:Flavin Mononucleotide