Abstract
We developed force field parameters for fluorinated, aromatic amino acids enabling molecular dynamics (MD) simulations of fluorinated proteins. These parameters are tailored to the AMBER ff15ipq protein force field and enable the modeling of 4, 5, 6, and 7F-tryptophan, 3F- and 3,5F-tyrosine, and 4F- or 4-CF(3)-phenylalanine. The parameters include 181 unique atomic charges derived using the implicitly polarized charge (IPolQ) scheme in the presence of SPC/E(b) explicit water molecules and 9 unique bond, angle, or torsion terms. Our simulations of benchmark peptides and proteins maintain expected conformational propensities on the μs time scale. In addition, we have developed an open-source Python program to calculate fluorine relaxation rates from MD simulations. The extracted relaxation rates from protein simulations are in good agreement with experimental values determined by (19)F NMR. Collectively, our results illustrate the power and robustness of the IPolQ lineage of force fields for modeling the structure and dynamics of fluorine-containing proteins at the atomic level.