Quantitative Prediction of Protein-Polyelectrolyte Binding Thermodynamics: Adsorption of Heparin-Analog Polysulfates to the SARS-CoV-2 Spike Protein RBD

Abstract

Interactions of polyelectrolytes (PEs) with proteins play a crucial role in numerous biological processes, such as the internalization of virus particles into host cells. Although docking, machine learning methods, and molecular dynamics (MD) simulations are utilized to estimate binding poses and binding free energies of small-molecule drugs to proteins, quantitative prediction of the binding thermodynamics of PE-based drugs presents a significant obstacle in computer-aided drug design. This is due to the sluggish dynamics of PEs caused by their size and strong charge-charge correlations. In this paper, we introduce advanced sampling methods based on a force-spectroscopy setup and theoretical modeling to overcome this barrier. We exemplify our method with explicit solvent all-atom MD simulations of interactions of anionic PEs that show antiviral properties, namely heparin and linear polyglycerol sulfate (LPGS), with the SARS-CoV-2 spike protein receptor binding domain (RBD). Our prediction for the binding free energy of LPGS to the wild-type RBD matches experimentally measured dissociation constants within thermal energy, kT, and correctly reproduces the experimental PE-length dependence. We find that LPGS binds to the Delta-variant RBD with an additional free-energy gain of 2.4 kT, compared to the wild-type RBD, in accord with electrostatic arguments. We show that the LPGS-RBD binding is solvent-dominated and enthalpy-driven, though with a large entropy-enthalpy compensation. Our method is applicable to general polymer adsorption phenomena and predicts precise binding free energies and re-configurational friction as needed for drug and drug-delivery design.