How to strike a conformational balance in protein force fields for molecular dynamics simulations?

Molecular dynamics (MD) simulation is a powerful tool for exploring the conformational energy landscape of proteins, and the reliability of MD results is crucially dependent on the underlying force field (FF). An accurate FF capable of producing balanced distributions of diverse conformations at multiple levels has been a long-sought goal. Towards this, several decades of joint efforts have been made to address FF deficiencies, manifested by conformational biases at different levels (local conformations, secondary structures, and global extendedness of polypeptide chain). We first present the major FF biases, then review the strategies to address them separately. Specifically, both nonresidue-specific and residue-specific strategies for torsional parameter optimization have been applied to achieve local conformation and secondary structure balances. Significant improvements can be gained with residue-specific torsional parameters especially when explicit dihedral couplings are considered. Further, the additional balance between protein–protein and protein–water interactions has been optimized via multiple ways to reproduce the global extendedness of polypeptide chains, especially for unfolded or disordered proteins. This review aims to summarize the most valuable experience and lessons gained from the past, which, we hope, can facilitate further improvements of both classical FFs and more sophisticated models such as polarizable FFs. This article is categorized under: Molecular and Statistical Mechanics > Molecular Dynamics and Monte-Carlo Methods Molecular and Statistical Mechanics > Molecular Mechanics.