Quantifying Protein-Protein Interactions in Molecular Simulations

23 December 2019, Version 1
This content is a preprint and has not undergone peer review at the time of posting.


We present simple, accurate, and efficient methods to estimate the dissociation constant Kd and the second osmotic virial coefficient B2 from molecular simulations. We show that for simulations of two proteins in a box, Kd is determined by B2 and the fraction of bound protein. We present two different methods to calculate B2 from Monte Carlo and molecular dynamics simulations using implicit or explicit solvent. We derive a surprisingly simple expression for B2, 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 B2 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.


dissociation constant
second osmotic virial coefficient
binding affinity
molecular simulation
molecular dynamics simulation
Monte Carlo simulation
analytical ultracentrifugation
equilibrium constant


Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.