Impact of Force Field Polarization on Correlated Motions of Proteins (2510.04897v1)

Abstract: Correlated motions of proteins underpin many physiological mechanisms, such as substrate binding, signal transduction, enzymatic activity and allostery. These motions arise from low frequency collective movements of biomolecules and have mostly been studied using molecular dynamics simulations. Here, we present the effects of two different empirical energy force fields used for molecular dynamics simulations on correlated motions -- the non-polarizable CHARMM 36m additive force field and the polarizable Drude-2019 force field. The study was conducted on two proteins, ubiquitin - a small protein with a well-described dynamic - and the nuclear receptor protein PPAR_. The ligand binding domain of PPAR_ was of particular interest since its function is to regulate transcription through ligand and coregulator protein binding. It has been previously shown that a dynamical network of correlated motions ensures the transmission of information related to PPAR___ ligand binding. We present the results of classical MD simulations where we analyze the results in terms of residue fluctuations, residue correlation maps, community network analysis and hydrophobic cluster analysis. We find that RMS fluctuations tend to be greater and correlated motions are less intense with Drude-2019 force field than with the non-polarizable all atom additive force field. Analysis of large hydrophobic clusters in the respective proteins show a greater loss of native contacts in the simulations using the Drude-2019 force field than in the simulations using the all atom force additive force field. Our results provide the first quantification of the impact of using a polarizable force field in computational studies that focus on correlated motions.