Quantifying Protein-Protein Interactions in Molecular Simulations

Abstract

<p><a></a></p><p><a></a>We present simple, accurate, and efficient methods to estimate the dissociation constant K<sub>d </sub>and the second osmotic virial coefficient B<sub>2 </sub>from molecular simulations. We show that for simulations of two proteins in a box, K<sub>d </sub>is determined by B<sub>2 </sub>and the fraction of bound protein. We present two different methods to calculate B<sub>2 </sub>from Monte Carlo and molecular dynamics simulations using implicit or explicit solvent. We derive a surprisingly simple expression for B<sub>2</sub>, adding significantly to the understanding of this important quantity. Non-binding interactions of proteins and other macromolecules shape the physicochemical properties of the crowded environments inside cells and of biomolecular condensates. We show how to extract the contributions of non-binding conformations to B<sub>2 </sub>and discuss how these can be determined in analytical ultracentrifugation and SAXS experiments. We expect that our methods will prove to be instrumental in force parameterization efforts and high-throughput studies of large interactomes. </p>