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

TL;DR

This study develops advanced simulation methods to accurately predict the thermodynamics of protein-polyelectrolyte interactions, exemplified by heparin-analog polysulfates binding to SARS-CoV-2 RBD, aiding drug design.

Contribution

The paper introduces a novel approach combining force-spectroscopy based sampling and theoretical modeling for quantitative thermodynamic predictions of PE-protein binding.

Findings

Predicted binding free energy of LPGS matches experimental data within thermal energy.

LPGS binds more strongly to Delta-variant RBD, with an additional 2.4 kT free-energy gain.

Binding is solvent-dominated and enthalpy-driven, with significant entropy-enthalpy compensation.

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.