Structural Insights into the Cofactor Role of Heparin/Heparan Sulfate in Binding between the SARS-CoV-2 Spike Protein and Host Angiotensin-Converting Enzyme II

PRO-2000 exhibits SARS-CoV-2 antiviral activity by interfering with spike-heparin binding

Here, we report on the anti-SARS-CoV-2 activity of PRO-2000, a sulfonated polyanionic compound. In Vero cells infected with the Wuhan, alpha, beta, delta or omicron variant, PRO-2000 displayed EC50 values of 1.1 μM, 2.4 μM, 1.3 μM, 2.1 μM and 0.11 μM, respectively, and an average selectivity index (i.e. ratio of cytotoxic versus antiviral concentration) of 172. Its anti-SARS-CoV-2 activity was confirmed by virus yield assays in Vero cells, Caco2 cells and A549 cells overexpressing ACE2 and TMPRSS2 (A549-AT). Using pseudoviruses bearing the SARS-CoV-2 spike (S), PRO-2000 was shown to block the S-mediated pseudovirus entry in Vero cells and A549-AT cells, with EC50 values of 0.091 μM and 1.6 μM, respectively. This entry process is initiated by interaction of the S glycoprotein with angiotensin-converting enzyme 2 (ACE2) and heparan sulfate proteoglycans. Surface Plasmon Resonance (SPR) studies showed that PRO-2000 binds to the receptor-binding domain (RBD) of S with a KD of 1.6 nM. Similar KD values (range: 1.2 nM-2.1 nM) were obtained with the RBDs of the alpha, beta, delta and omicron variants. In an SPR neutralization assay, PRO-2000 had no effect on the interaction between the RBD and ACE2. Instead, PRO-2000 was proven to inhibit binding of the RBD to a heparin-coated sensor chip, yielding an IC50 of 1.1 nM. To conclude, PRO-2000 has the potential to inhibit a broad range of SARS-CoV-2 variants by blocking the heparin-binding site on the S protein.

Well-Defined Heparin Mimetics Can Inhibit Binding of the Trimeric Spike of SARS-CoV-2 in a Length-Dependent Manner

The emergence of new SARS-CoV-2 variants and the dangers of long-covid necessitate the development of broad-acting therapeutics that can reduce viral burden. SARS-CoV-2 employs heparan sulfate (HS) as an initial cellular attachment factor, and therefore, there is interest in developing heparin as a therapeutic for SARS-CoV-2. Its use is, however, complicated by structural heterogeneity and the risk of causing bleeding and thrombocytopenia. Here, we describe the preparation of well-defined heparin mimetics by a controlled head-to-tail assembly of HS oligosaccharides having an alkyne or azide moiety by copper-catalyzed azide-alkyne cycloaddition (CuAAC). Alkyne- and azide-containing sulfated oligosaccharides were prepared from a common precursor by modifying an anomeric linker with 4-pentynoic acid and by enzymatic extension with an N-acetyl-glucosamine having an azide moiety at C-6 (GlcNAc6N₃), respectively, followed by CuAAC. The process of enzymatic extension with GlcNAc6N₃ followed by CuAAC with the desired alkyne-containing oligosaccharides could be repeated to give compounds composed of 20 and 27 monosaccharides, respectively. The heparin mimetics could inhibit the binding of the SARS-CoV-2 spike or RBD to immobilized heparin or to Vero E6 cells. The inhibitory potency increased with increasing chain length, and a compound composed of four sulfated hexasaccharides linked by triazoles had a similar potency as unfractionated heparin. Sequence analysis and HS microarray binding studies with a wide range of RBDs of variants of concern indicate that they have maintained HS-binding capabilities and selectivities. The heparin mimetics exhibit no or reduced binding to antithrombin-III and platelet factor 4, respectively, which are associated with side effects.

SARS-CoV-2 evolved variants optimize binding to cellular glycocalyx

Viral variants of concern continue to arise for SARS-CoV-2, potentially impacting both methods for detection and mechanisms of action. Here, we investigate the effect of an evolving spike positive charge in SARS-CoV-2 variants and subsequent interactions with heparan sulfate and the angiotensin converting enzyme 2 (ACE2) in the glycocalyx. We show that the positively charged Omicron variant evolved enhanced binding rates to the negatively charged glycocalyx. Moreover, we discover that while the Omicron spike-ACE2 affinity is comparable to that of the Delta variant, the Omicron spike interactions with heparan sulfate are significantly enhanced, giving rise to a ternary complex of spike-heparan sulfate-ACE2 with a large proportion of double-bound and triple-bound ACE2. Our findings suggest that SARS-CoV-2 variants evolve to be more dependent on heparan sulfate in viral attachment and infection. This discovery enables us to engineer a second-generation lateral-flow test strip that harnesses both heparin and ACE2 to reliably detect all variants of concern, including Omicron.

Heparan Sulfate and Enoxaparin Interact at the Interface of the Spike Protein of HCoV-229E but Not with HCoV-OC43

It is known that the spike protein of human coronaviruses can bind to a secondary receptor, or coreceptor, to facilitate the virus entry. While HCoV-229E uses human aminopeptidase N (hAPN) as a receptor, HCoV-OC43 binds to 9-O-acetyl-sialic acid (9-O-Ac-Sia), which is linked in a terminal way to the oligosaccharides that decorate glycoproteins and gangliosides on the surface of the host cell. Thus, evaluating the possible inhibitory activity of heparan sulfate, a linear polysaccharide found in animal tissues, and enoxaparin sodium on these viral strains can be considered attractive. Therefore, our study also aims to evaluate these molecules’ antiviral activity as possible adsorption inhibitors against non-SARS-CoV. Once the molecules’ activity was verified in in vitro experiments, the binding was studied by molecular docking and molecular dynamic simulations confirming the interactions at the interface of the spike proteins.

Entropic Overcompensation of the N501Y Mutation on SARS-CoV-2 S Binding to ACE2

Recent experimental work has shown that the N501Y mutation in the SARS-CoV-2 S glycoprotein's receptor binding domain (RBD) increases binding affinity to the angiotensin-converting enzyme 2 (ACE2), primarily by overcompensating for a less favorable enthalpy of binding by greatly reducing the entropic penalty for complex formation, but the basis for this entropic overcompensation is not clear [Prévost et al. J. Biol. Chem.2021, 297, 101151]. We use all-atom molecular dynamics simulations and free-energy calculations to qualitatively assess the impact of the N501Y mutation on the enthalpy and entropy of binding of RBD to ACE2. Our calculations correctly predict that N501Y causes a less favorable enthalpy of binding to ACE2 relative to the original strain. Furthermore, we show that this is overcompensated for by a more entropically favorable increase in large-scale quaternary flexibility and intraprotein root mean square fluctuations of residue positions upon binding in both RBD and ACE2. The enhanced quaternary flexibility stems from N501Y's ability to remodel the inter-residue interactions between the two proteins away from interactions central to the epitope and toward more peripheral interactions. These findings suggest that an important factor in determining protein-protein binding affinity is the degree to which fluctuations are distributed throughout the complex and that residue mutations that may seem to result in weaker interactions than their wild-type counterparts may yet result in increased binding affinity thanks to their ability to suppress unfavorable entropy changes upon binding.

Computational Analysis of the SARS-CoV-2 RBD-ACE2-Binding Process Based on MD and the 3D-RISM Theory

The binding process of angiotensin-converting enzyme 2 (ACE2) to the receptor-binding domain (RBD) of the severe acute respiratory syndrome-like coronavirus 2 spike protein was investigated using molecular dynamics simulation and the three-dimensional reference interaction-site model theory. The results suggested that the protein-binding process consists of a protein-protein approaching step, followed by a local structural rearrangement step. In the approaching step, the interprotein interaction energy decreased as the proteins approached each other, whereas the solvation free energy increased. As the proteins approached, the glycan of ACE2 first established a hydrogen bond with the RBD. Thereafter, the number of interprotein hydrogen bonds increased rapidly. The solvation free energy increased because of the desolvation of the protein as it approached its partner. The spatial distribution function of the solvent revealed the presence of hydrogen bonds bridged by water molecules on the RBD-ACE2 interface. Finally, principal component analysis revealed that ACE2 showed a pronounced conformational change, whereas there was no significant change in RBD.