Skip to main content

MINI REVIEW article

Front. Mol. Biosci., 25 January 2021
Sec. Molecular Recognition
This article is part of the Research Topic Interactions Between Proteins and Biomacromolecules: Tools and Applications View all 11 articles

Elucidating the Interactions Between Heparin/Heparan Sulfate and SARS-CoV-2-Related Proteins—An Important Strategy for Developing Novel Therapeutics for the COVID-19 Pandemic

\nMingjia Yu&#x;Mingjia Yu1Tianji Zhang&#x;Tianji Zhang2Wei ZhangWei Zhang2Qianyun SunQianyun Sun3Hongmei Li
Hongmei Li2*Jin-ping Li,
Jin-ping Li1,4*
  • 1Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, Beijing, China
  • 2Division of Chemistry and Analytical Science, National Institute of Metrology, Beijing, China
  • 3Division of Chemistry, Shandong Institute of Metrology, Jinan, China
  • 4Department of Medical Biochemistry and Microbiology, University of Uppsala, Uppsala, Sweden

Owing to the high mortality and the spread rate, the infectious disease caused by SARS-CoV-2 has become a major threat to public health and social economy, leading to over 70 million infections and 1. 6 million deaths to date. Since there are currently no effective therapeutic or widely available vaccines, it is of urgent need to look for new strategies for the treatment of SARS-CoV-2 infection diseases. Binding of a viral protein onto cell surface heparan sulfate (HS) is generally the first step in a cascade of interaction that is required for viral entry and the initiation of infection. Meanwhile, interactions of selectins and cytokines (e.g., IL-6 and TNF-α) with HS expressed on endothelial cells are crucial in controlling the recruitment of immune cells during inflammation. Thus, structurally defined heparin/HS and their mimetics might serve as potential drugs by competing with cell surface HS for the prevention of viral adhesion and modulation of inflammatory reaction. In this review, we will elaborate coronavirus invasion mechanisms and summarize the latest advances in HS–protein interactions, especially proteins relevant to the process of coronavirus infection and subsequent inflammation. Experimental and computational techniques involved will be emphasized.

Introduction

Currently, the whole world is facing the deadly coronavirus disease 2019 (COVID-19) outbreak caused by the coronavirus (CoV) SARS-CoV-2, which has been far beyond the outbreaks caused by the other two major coronaviruses (SARS and MERS) in the past 20 years (Drosten et al., 2003; Zaki et al., 2012). So far, there are no specific therapeutic and effective drugs available. Vaccines, although achieving success worldwide, are still far from being widely accessible. In this regard, multidimensional antiviral strategies are strongly needed in preventing the spread of COVID-19 and treating infected individuals.

Glycosaminoglycans (GAGs) are a group of anionic polysaccharides composed of repeating disaccharide building blocks, including heparin/heparan sulfate (HS) (-4GlcAβ/IdoAα1-4GlcNxα1-, x = Ac, SO3H or H), chondrotin/dermatan sulfate (-4GlcAβ/IdoAα1-3GalNxβ1-), keratan sulfate (-3GalAβ1-4GlcNAcβ1-), and hyaluronic acid (HA) (-4GlcAβ1-3GlcNAcβ1-) (Figure 1). Sulfation at various positions of the sugar residues could occur except for the HA, making their structures heterogeneous and extremely difficult to characterize (Tianji Zhang et al., 2019). Among them, heparin and HS exhibit the most diverse biological activities, most of which are mediated by their interactions with proteins (Li and Kusche-Gullberg, 2016). Recent work identified HS on the cell surface as a co-receptor for the SARS-CoV-2 spike protein (S protein) (Clausen et al., 2020), making the HS–S protein interaction an extremely appealing target for manipulating SARS-CoV-2 infection. Designing competing HS mimetics requires the elucidation of the mode of interaction, particularly sequence specificities of the HS.

FIGURE 1
www.frontiersin.org

Figure 1. Structure of GAGs. Structures of heparin/heparan sulfate (A), chondrotin/dermatan sulfate (B), keratan sulfate (C), and hyaluronic acid (D) were illustrated. Reprinted from ref Tianji Zhang et al. (2019) with permission.

In this review, we will summarize the classification and invasion mechanisms of the major CoVs and elaborate possible antiviral strategies based on the interactions between heparin/HS and proteins. Relevant technologies involved in elucidating heparin/HS–protein interactions are crucial for developing (sequence) specific antiviral molecules and will thus be underlined.

Classification and Invasion Mechanisms of the Coronavirus

CoVs are highly diverse, enveloped, and positive-sense single-stranded (up to 30,000 bp) RNA viruses (Coutard et al., 2020) belonging to the Nidovirales order in the subfamily of Othocoronavirinae (Wang et al., 2020). Infection by the viruses can cause severe diseases affecting upper respiratory, gastrointestinal, and central nervous systems in humans and other animals (Gallagher and Buchmeier, 2001). Based on systematic analysis of viral nucleic acid sequence, CoVs can be classified into four genera: alpha, beta, gamma, and delta according to the 10th Report on Virus Taxonomy from the International Committee on Taxonomy of Viruses (ICTV) (Fehr and Perlman, 2015). Among them, alpha- and beta-coronaviruses can infect mammals, gamma-coronaviruses can infect avian species, while delta-coronaviruses can infect both (Li, 2016). Currently, there are seven representative strains of human coronaviruses (HCoVs) including four low-pathogenic coronaviruses [HCoV-229E, HCoV-NL63 (alpha-coronaviruses), HCoV-OC43, and HCoV-HKU1 (beta-coronaviruses)], which cause mild respiratory diseases in humans (Su et al., 2016), and three high pathogenic strains including HCoVs {severe acute respiratory syndrome coronavirus (SARS-CoV) (Drosten et al., 2003), Middle East respiratory syndrome coronavirus (MERS-CoV) (Elfiky et al., 2017), and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (beta-coronaviruses) (Hui et al., 2020)}. The three highly pathogenic strains have caused deadly pneumonia in humans since the beginning of the twenty-first century. Unfortunately, so far, there are still no specific therapeutics approved against these human-infecting coronaviruses, mainly due to lacking sufficient knowledge in the pathological process of viral infection. Thus, in-depth understanding of the infection mechanisms will facilitate the development of effective interventions against these highly pathogenic coronaviruses and are of high urgency for the control and treatment of COVID-19.

CoVs share similar genome identities. Two-thirds of the genome at the 5′-terminus contain two large overlapping open reading frames (ORFs), ORF 1a and ORF 1b, which encode polyproteins 1a (pp1a) and pp1b/1ab, respectively. The polyproteins can be further cleaved into 15–16 non-structural proteins (nsp2-nsp16 or nsp1-nsp16). One-third of the genome at the 3′-terminus encodes four common structural proteins in the order of Spike (S) that characterizes all coronaviruses, Envelope (E), Membrane (M), and Nucleocapsid (N) (Wang et al., 2020) (Figure 2A). The S protein is a trimeric class I fusion protein that protrudes from the virion surface and mediates receptor recognition, membrane fusion, virus entry, and antibody neutralization (Gallagher and Buchmeier, 2001). Considering its significant functions during viral infection (Liu et al., 2004), the S protein serves as a main target for the development of antiviral drugs (Du et al., 2017), antibodies (He et al., 2006), entry inhibitors (Lu et al., 2014), and vaccines (Du et al., 2009). Each monomer of the trimeric S protein is ~180 kDa containing two subunits—a receptor-binding subunit (S1) and a membrane-fusion subunit (S2), which are linked through a fusion peptide. The S1 subunit contains two independent domains—the N-terminal domain (NTD) and the C-terminal domain (C-domain) (Ou et al., 2020) (Figure 2B), either of which can serve as the receptor-binding domain (RBD) depending on the virus strains (Kubo et al., 1994; Ou et al., 2017). The S2 subunit consists of four main domains—the heptad repeat 1 (HR1) domain, heptad repeat 2 (HR2) domain, transmembrane domain (TM), and cytoplasm domain (CP) (Xia et al., 2020). During viral infection, a two-step sequential protease cleavage process triggers the activation of S proteins (Belouzard et al., 2009; Millet and Whittaker, 2014), which is modulated by host range and cell tropism. The first cleavage occurs between the S1 and S2 subunits, leading to the release of the S1 subunit and its transition to the post-fusion conformation (Su et al., 2016). Then, as the RBD of the S1 subunit binds to a host cell receptor [CoVs can recognize both angiotensin-converting enzyme 2 (ACE2) and dipeptidyl peptidase 4 (DPP4, also known as CD26)] (Kuhn et al., 2004; Raj et al., 2013), another cleavage site on S2 is exposed and cleaved by host proteases at the S2' site located upstream of the fusion peptide (FP). The HR1 and HR2 domains of S2 form a six-helix bundle (6-HB) fusion core in order to bring viral and cellular membranes into close proximity for subsequent fusion and infection (Figure 3) (Bosch et al., 2004). The host proteases such as furin (Millet and Whittaker, 2014), cathepsins (Bertram et al., 2013), human airway trypsin-like protease (Berman et al., 2000; Bertram et al., 2011), and transmembrane protease serine protease-2 (TMPRSS-2) (Gierer et al., 2013) are widely expressed in many important organs, which is a critical reason for the systematic infection, serious pathogenicity, and high mortality of the CoVs. Therefore, the RBD and 6-HB fusion core of CoVs and the proteases on infected cells have become potential targets for the development of virus attachment/fusion inhibitors, neutralizing antibodies, and vaccines.

FIGURE 2
www.frontiersin.org

Figure 2. Schematic structure of CoVs and the S protein. (A) Schematic structure of virion of CoVs and their four structural proteins, including spike (S), envelope (E), membrane (M), and nucleocapsid (N) proteins. (B) Schematic representation of S protein of CoVs, including receptor-binding domain (RBD), N-terminal domain (NTD), or C-terminal domain (C-domain), fusion peptide (FP), heptad repeat 1 (HR1), heptad repeat 2 (HR2), transmembrane domain (TM), and cytoplasm domain (CP).

FIGURE 3
www.frontiersin.org

Figure 3. Putative antiviral mechanism of CoVs during viral entry, including the cleavage of S protein, activation of S2 domain, and the fusion of viral and host membranes.

Invasion of CoVs occurs in two steps, initial binding to the receptor on the cell surface and fusion of S protein with the host cell membrane to deliver their nucleocapsid to the target cell. It has been known that CoVs often initially interact with cell surface molecules to promote their binding to specific receptors. HS proteoglycans (HSPGs) are abundantly present in almost all mammalian cells and serve as a co-receptor for a number of viruses (Gomes and Dietrich, 1982). HSPGs could initially bind to the surface proteins of CoVs, promote subsequent recognition of a secondary Receptor (ACE2/DPP4), and facilitate the attachment and entry of virus by increasing their local concentrations. Studies also suggested that different compositions in HS could impact the tropism of viruses (Wickramasinghe et al., 2011) and HS co-receptors on host cell surface leads to conformational changes of the CoVs' S proteins (Lang et al., 2011; Milewska et al., 2014; Mycroft-West et al., 2020a), possibly through the formation of a ternary complex (Clausen et al., 2020) (Figure 3). These findings suggested that the HS–S protein interactions might serve as a potential target to attenuate virus infection.

Heparin/HS and their Interactions With Protein

Structures of Heparin/HS

Heparin is a significant anti-coagulant that has been used in clinic over decades. The heparin polysaccharide chains are linear and polyanionic, with repeating disaccharide units of α-L-iduronic acid (IdoA) or β-D-glucuronic acid (GlcA) residue linked to glucosamine (GlcN) residue by a 1-4 glycosidic bond. The sugar units are sulfated at N-, 6- and 3-O on the GlcN residues as well as 2-O on the hexuronic acid by site-specific sulfotransferases. The 3-O-sulfation is rare but critical for heparin to form a specific pentasaccharide domain that specifically bind to anti-thrombin with high affinity, which is essential for its anti-coagulant activity (Lindahl et al., 1980). In addition to the anti-coagulant activity, heparin and its derivatives have also been studied for their anti-inflammatory, antiviral, anti-angiogenesis, anti-neoplastic, and anti-metastatic effects (Hao et al., 2019). HS shares high structural similarity with heparin (Linhardt and Toida, 2004), but generally with lower level of sulfation and epimerization, therefore displaying distinct domain structures (Figure 4). The functionalities of heparin and HS are mediated by their interaction with various proteins including proteases, protease inhibitors, chemokines, cytokines, growth factors, and their respective receptors (Xu and Esko, 2014; Seffer et al., 2020), with variable specificities.

FIGURE 4
www.frontiersin.org

Figure 4. Synthesis and structures of heparin/HS.

An Overview of Heparin/HS and Protein Interactions

Due to the highly anionic nature of heparin/HS, the interactions between heparin/HS and proteins are primarily through the interaction between negatively charged sulfate and carboxyl groups on heparin/HS and positively charged lysine and arginine residues on the proteins. The role of electrostatics in heparin/HS–protein interactions was elucidated in several studies (Olson et al., 1991; Thompson et al., 1994; Friedrich et al., 2001). Meanwhile, nonionic interactions such as hydrogen bonding and van der Waals packing also contribute to the free energy for the binding reactions (Thompson et al., 1994).

Some heparin/HS binding proteins can be identified by amino acid sequences known as Cardin-Weintraub motifs corresponding to “XBBXBX” and “XBBBXXBX”, where X is a hydropathic residue and B is a basic residue, such as arginine and lysine, responsible for interacting with the sulfate groups present on heparin/HS (Cardin and Weintraub, 1989; Hileman et al., 1998). On the other hand, well-characterized heparin/HS–protein interactions revealed specific requirement of the carbohydrate sequence. The most prominent example is the binding of antithrombin with the unique pentasaccharide sequence, -GlcNS/Ac6S-GlcA-GlcNS3S6S-IdoA2S-GlcNS6S- in heparin, where the 3-O-sulfation is critical (Richard et al., 2009). Unlike the extremely rigorous sequence requirement as of antithrombin, or the purely non-specific interaction as in the case of heparin and protamine (Hubbard and Jennings, 1985), majority of the heparin/HS–protein interactions are selective. For instance, the interaction between HS and FGF2, a member of the fibroblast growth factor family, prefers the disaccharide unit of IdoA2S and GlcNS on heparin/HS (Turnbull et al., 1992; Jemth et al., 2002). More evidence is emerging, indicating that binding of HS and proteins is somewhat between purely specific and generally non-specific (Forsten-Williams et al., 2008; Nugent et al., 2013). Non-specific bindings solely depend on the high negative charge density of the carbohydrate chain and positively charged residues of the proteins, while modifications or domains on heparin/HS determine specificity levels of the interactions (Xu and Esko, 2014).

Experimental Technologies Involved in Studying Heparin/HS–Protein Interactions

The major challenge for elucidating heparin/HS–protein interactions is to decipher the carbohydrate sequence that is commonly of high heterogeneity. Multidimensional technologies that facilitate understanding sequence specificities have been comprehensively summarized in a recent review (Yang and Chi, 2017), including X-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy, and mass spectrometry (MS).

Heparin/HS oligosaccharide microarrays are valuable tools that can be used to probe the interactions between structurally defined oligosaccharides and proteins with relatively small amounts of samples. The bottleneck of the microarray assay is to synthesize oligosaccharide libraries of intensive diversity. Zong et al. prepared a tetrasaccharide library consisting of 47 unique structures, which is one the most comprehensive HS microarrays covering a large portion of possible structural variabilities (Zong et al., 2017). In a recent study, chemoenzymatic strategies have been successfully applied to construct microarrays composed of tetrasaccharide to 18-mer containing various N-, 6-O-, 2-O-, and 3-O-sulfation modifications (Horton et al., 2020). Cell-based microarrays have also been developed, aiming at demonstrating the functionality of specific heparin/HS saccharides in real cell signaling (Puvirajesinghe et al., 2012; Sterner et al., 2013).

The surface plasmon resonance (SPR) sensor (Thompson et al., 1994) is one of the most convenient tools for detecting heparin/HS–protein interactions through changes of the refractive index signals. One of the major advantages of SPR is the capacity of probing biomolecular interactions at the thermodynamic level, offering real-time and label-free measurement of reaction rate constants (kon, koff) and resultant equilibrium constants (KA, KD) (Homola, 2008). In a recent study, interactions between heparin/HS and various cytokines were characterized by coupling surface plasmon resonance imaging for thermodynamic analysis method and Matrix-Assisted Laser Desorption/Ionization Time of Flight Mass Spectrometry (MALDI-TOF-MS) for structural determination (Przybylski et al., 2020). A self-assembled monolayer of short polyethylene oxide chains was used for grafting cytokines. Captured carbohydrates were carried out directly on the biochip surface using MALDI-TOF-MS, while MS identification was enhanced by on-chip digestion of the cytokine-bound GAGs by heparinase treatment.

Computational Techniques

Despite the advances in the experimental techniques, there are limitations (e.g., failure in the acquisition of co-crystal structures) in obtaining the information regarding molecular interactions. Thus, computational techniques are indispensable tools for comprehensively understanding the heparin/HS–protein interactions. On the basis of theoretical models, the computational techniques are especially helpful in designing novel drug, performing wide scale analysis against large database (e.g., the PDB database), and for understanding the interaction dynamics. The interactions between protein and heparin/HS can be weak or strong, transient or stable, non-permanent or permanent, which can be evaluated by basic parameters of binding affinities such as the equilibrium dissociation constant, KD (Ma et al., 2018); Gibbs free energy of binding, ΔG (Steinbrecher and Labahn, 2010); inhibition constant, Ki (Pekkarinen et al., 2007); half maximal inhibitory concentration, IC50 (Sebaugh, 2011); and electrostatic potential energy (Bitencourt-Ferreira et al., 2019). Several techniques are available for a wide application in assessment of protein-heparin/HS interaction.

Homology Modeling

The homology modeling method uses 3D structures deposited in the PDB database to predict protein structures of sequential similarities. Homology modeling can give spatial structures with the highest accuracy (Werner et al., 2012) and thus has been widely applied for rational analysis of interactions between small organic molecule (ligand) and target protein during the docking and virtual screening for drug discovery (Cheng et al., 2012). Homology modeling can be built by four methods, including rigid body assembly [by tools like SWISS-MODEL (Arnold et al., 2006)], segment matching [by tools like SEGMOD/ENCAD (Levitt, 1992)], spatial restraint [by tools like MODELER (Sali and Blundell, 1993)], and artificial evolution [by tools like NEST (Petrey et al., 2003)].

Molecular Docking

Molecular docking is a computational procedure extensively used in novel drug discovery, through which a small molecule (ligand) is docked into a macromolecule (target protein) at the binding sites to predict the binding conformations and affinity. The conformation of a ligand binding with the receptor depends on its state variables (including position: x-, y-, and z-translations; orientation: euler angles, axis angle, and quaternion; and conformation: torsion angles for each rotatable bond), which decides the extent of the multidimensional search space within the protein–ligand interaction. All docking methods require a scoring function to rank various candidate protein–ligand binding modes and a search method to explore their state variables. Scoring functions are computational approximations to predict protein–ligand binding affinity based on empiricism, force field, and knowledge, while search methods are classified into local and global ones by the extent of search space. Local search methods [such as solis and wets (Solis and Wets, 1981) and the pattern search (Lewis and Torczon, 2002)] tend to find the nearest or local minimum energy to the current conformation, whereas global methods [such as Monte Carlo (MC) simulated annealing (SA) (Kirkpatrick et al., 1983) and the genetic algorithm (GA) (Goldberg, 1988)] search for the best or global minimum energy within the defined search space. Hybrid global-local search methods have been shown to perform even better due to the higher efficiency in finding lower energies among the different candidate protein–ligand binding modes (Morris and Lim-Wilby, 1999). Molecular docking can be performed by various docking programs such as AutoDock (Goodsell et al., 1996), AutoDock Vina (Trott and Olson, 2010), FlexX (Rarey et al., 1996), GOLD (Jones et al., 1997), and Molegro Virtual Docker (MVD) (De Azevedo, 2010).

Electrostatic Potential Energy

Since electrostatic interactions are part of scoring functions that widely influence the binding affinity of protein–ligand complexes, electrostatic potential energies are calculated in computational models to compare protein–ligand binding affinities. The electrostatic force is conservative as it only depends on the initial and final positions and most protein–ligand complexes show only partial charges. Partial Equalization of Orbital Electronegativity (PEOE) is the most widely used method (Gasteiger and Marsili, 1980) provided by the software AutoDockTools4 (Morris et al., 2009) that estimates partial charges of target protein and ligands in order to calculate electrostatic potential energy.

Molecular Dynamics

Molecular dynamics (MD) is a computational simulation technique that can obtain not only multiple conformations of target proteins and ligands, but also a wealth of energetic status about the interactions in a time-dependent manner. MD simulations combined with binding free energy calculations can improve the accuracy of binding prediction and are thus suitable for studying the motions of target protein on ligand binding (Radkiewicz and Brooks, 2000; Salsbury et al., 2001). The Newtonian equation of motions is applied for each atom in the MD simulations for approximations (Schlick, 2010), which requires the information of initial coordinates (obtained from experimental structures, models, or combination of the two), potential (obtained from different force fields along with the coordinates) (MacKerell et al., 1998), and algorithms. Given the diverse complexity of the protein–ligand structures, different force field models [such as CHARMM (Miller et al., 2008), AMBER (Guvench and MacKerell, 2008), and GROMACS (Van Der Spoel et al., 2005)] are flexibly used during simulation, which are associated with modeling suites of CHARMM (Brooks et al., 2009), AMBER (Case et al., 2005), GROMACS (Hess et al., 2008), and NAMD (Phillips et al., 2005). Owing to the advances in computers and algorithms, the complexes of biomacromolecules can be simulated in nanoseconds with whole atoms, generating numerous conformations. The characterization of each conformation is accomplished by sophisticated methods that could be divided up to four types, including gross measures of protein and simulation stability, clustering analysis, quasiharmonic and principal component analysis, and correlation function analysis. The gross measures of protein and simulation stability is the most widely used approach for checking simulation integrity and estimating equilibration timescale of the simulation. Parameters such as root-mean-square deviation (RMSD), structural clustering, free energy of binding and native contacts, and average temperature and pressure are generally calculated for obtaining their fluctuations.

Docking an HS fragment of proper size (≥4 monosaccharide units) to a protein is challenging due to the flexibility brought by the glycosidic rings, linkages, and the high density of negative charges. Sapay et al. proposed a two-step method based on molecular docking and MD simulation to explore the binding modes of HS to cellular growth factors (FGF2 and CXCL12α) (Sapay et al., 2011). The method provided dynamical modeling of the protein–ligand complex by building the docking models of HS fragment on protein surface and refining the contacts between HS fragment and the protein.

Computational Study of SARS-CoV-2 Infection

The powerful tools of computational technology has made significant contribution to the studies on viruses including SARS-CoV-2. The method of homology modeling has made initially important contribution. Based on the rich genomic information and bioinformatics analysis of the proteins encoded by the novel coronavirus genes, Wu et al. built 19 structures of SARS-CoV-2 by homology modeling through the Fold and Function Assignment System server, including viral papain like protease (PLpro), main protease (3CLpro, also named 3-chymotrypsin-like protease), RNA-dependent RNA polymerase (RdRp), helicase, and S protein (Wu et al., 2020). As the most efficient way to find anti-SARS-CoV-2 drugs is to screen those that are commonly used in clinic, small-molecule compounds from several resources including the U.S Food and Drug Administration (FDA)-approved drug database (ZINC drug database, ZDD), traditional Chinese medicine/natural products database, and the antiviral drugs database were docked into these computational models by ICM 3.7.3 modeling software (MolSoft LLC) to virtually screen potential druggable targets. Successfully predicted targets and potential drug compounds can be further tested in vitro and in vivo for treating SARS-CoV-2 infections.

The calculations of electrostatic potential energy were performed to estimate protein-heparin/HS binding affinities combined with the docking technique. Clausen et al. calculated the electrostatic potential map of both SARS-CoV-1 RBD (PDB ID: 3BGF) (Pak et al., 2009) and SARS-CoV-2 RBD (PDB ID: 6M17) (Yan et al., 2020) by the Molecular Operating Environment (MOE) software (Clausen et al., 2020). Combined with docking studies of oligosaccharide fragments derived from heparin with RBD, it revealed an extended electropositive surface of RBD composed of positively charged residues including R346, R355, K444, R466, and possibly R509 that could coordinate the electronegative oligosaccharides through hydrogen bonds and hydrophobic interactions. This study demonstrated that the SARS-CoV-2 S protein may mediate an enhanced interaction with HS analogs and heparinoid derivatives compared to SARS-CoV-1 by the evolution of Lys444 and Glu354 (SARS-CoV-1) to Thr444 and Asn354 (SARS-CoV-2), respectively.

MD simulations have also been proven as a convenient method to describe the motions and binding affinities of ligand into target proteins. Han et al. designed and simulated several potential peptide inhibitors (including virus-binding domain α-helices extracted from the protease domain of ACE2) against the SARS-CoV-2 coronavirus (Han and Kral, 2020). Classical MD simulations were performed by the modeling suites NAMD (Phillips et al., 2005) and CHARMM36 protein force field (MacKerell et al., 1998), which screened the most suitable peptide inhibitor with good binding affinity yet low RMSD for critical amino acids, indicative of relatively high binding energies. The novel designed peptide inhibitors have provided insights for researchers to develop therapeutic antiviral inhibitors by offering the α1 helix of ACE2 a sulfated ligand. Other molecules of similar structures, the heparin/HS for instance, could also attach to positively charged residues at the bottom of the RBD.

The timescale of the MD simulations is also a determinant for the convergence of structural clustering, free energy of binding, and native contacts between the GAGs and target proteins. Bojarski et al. analyzed the structure of fibroblast growth factor 1 (FGF1) complexed with heparin [PDB ID: 2AXM (DiGabriele et al., 1998)] through microsecond-scale simulations by the force field of AMBER16 (Bojarski et al., 2019). The analysis revealed a conformational selection mechanism of GAGs binding and determined the structural specificity in the FGF1–heparin complex. Their findings could potentially contribute to the development of novel biomaterial resembling GAGs in the field of regenerative medicine.

Utilizing Heparin/HS–Protein Interactions to Explore Novel Strategies for Treatment of SARS-CoV-2 Infections

Interactions Between Heparin/HS and the S Protein

Based on the knowledge of virus–heparin/HS interaction, it is assumed that exogenous added heparin/HS may interfere with viral infection. Several excellent studies have focused on the interactions between heparin/HS and the SARS-CoV-2 S protein, especially structure specificity of the carbohydrate chains. Nevertheless, it is worth noting that results from different research groups exhibited inconsistency to some extent, even if the similar analytical methods were conducted. This could be possibly attributed to experimental parameter setting, and the complexity and heterogeneity of the heparin/HS structures.

Using the SPR technique, Mycroft-West et al. first reported the SARS-CoV-2 S1 RBD binding to unfractionated heparin (Mycroft-West et al., 2020a). Through circular dichroism (Martino et al., 2020) spectroscopy, the authors further indicated that the RBD underwent conformational change in the presence of heparin, including helix and beta-sheet content alterations. The changes demonstrated that the RBD interacted with heparin in aqueous solution of physiological significance, whereby the major changes induced by heparin were those associated with antiparallel and helix content. In a subsequent study (Mycroft-West et al., 2020b), the authors found that the addition of heparin to Vero cells between 6.25 and 200 μg ml−1 inhibited invasion of SARS-CoV-2 by 44–80%. Additionally, SPR data revealed that 2-O, 6-O, and completely desulfated heparin had no inhibitory activity on heparin–RBD binding, proving the significance of 2-O and 6-O-sulfation on heparin/HS-spike interactions. On the other hand, persulfated, N-desulfated/N-re-acetylated, and strikingly 2-O/6-O doubly desulfated heparin possessed inhibitory activity. The authors attributed this phenomenon to a preference of RBD for a particular spatial arrangement of charged groups.

In another study, Kim et al. found that both monomeric and trimeric SARS-CoV-2 S proteins bound to immobilized heparin (KD = 40 and 73 pM, respectively) more tightly than the SARS-CoV and MERS-CoV S proteins (500 and 1 nM, respectively) (Kim et al., 2020). Heparin-derived oligosaccharides (dp 4 to 18), N-desulfated, 2-O-desulfated, and 6-O-desulfated heparin failed to compete with immobilized heparin for binding to the S protein, suggesting that chain length and all the sulfate groups within heparin were critical in the interaction. On the other hand, in responses to heparin, tri-sulfated non-anti-coagulant HS, and non-anti-coagulant low-molecular-weight heparin (LMWH), the binding of SARS-CoV-2 S protein to the surface-immobilized heparin decreased in a concentration-dependent fashion. The IC50 were determined to be 0.056, 0.12, and 26.4 μM, respectively. Additionally, unbiased computational ligand docking predicted putative heparin/HS-binding motifs on the S protein: 453–459 (YRLFRKS), 681–686 (PRRARS), and 810–816 (SKPSKRS), among which the 681–686 (PRRARS) site between the S1 and S2 subunits was a novel insertion not present in the SARS and MERS S proteins.

Liu et al. performed microarray binding experiments using an extensive HS oligosaccharide library (Liu et al., 2020). Their data suggested that the SARS-CoV-2 S protein can bind HS in a length- and sequence-dependent manner, while hexa- and octasaccharides composed of IdoA2S-GlcNS6S repeating units were identified as optimal ligands. Notably, 3-O-sulfation on the GlcN residue was proven not essential for efficient binding. In support of the microarray data, SPR experiments showed that the SARS-CoV-2 S protein bound with higher affinity to heparin (KD = 55 nM) compared to the RBD (KD = 1 μM) alone. The previously determined octasaccharide composed of IdoA2S-GlcNS6S repeating subunits could inhibit heparin–S protein interaction with an IC50 of 38 nM. Their data supported a model in which the RBD of the SARS-CoV-2 S protein conferred sequence specificity and agreed with Kim et al. (2020) where an additional HS binding site in the S1/S2 proteolytic cleavage site could enhance the avidity of binding.

To obtain insights into heparin/HS–S protein binding and virus infection in a safer circumstance, Tandon et al. pseudotyped SARS-CoV-2 S protein on a third-generation lentiviral (pLV) vector for testing the impact of various sulfated polysaccharides on transduction efficiency in mammalian cells (Tandon et al., 2020). The pLV vector pseudotyped the S protein efficiently and produced high titers on HEK293T cells. Both unfractionated heparin (UFH) and enoxaparin (a low-molecular-weight heparin drug) exhibited high apparent inhibitory activity. However, in contrast with previous SPR results, the authors found that selective desulfation at the 6-O-position of the GlcN residue did not significantly reduce inhibitory activity of either UFH or enoxaparin. Concentration–response curves showed that pLV-S particles were efficiently neutralized by a range of concentrations of UFH, enoxaparin, 6-O-desulfated UFH, and 6-O-desulfated enoxaparin with an IC50 of 5.99 μg/L, 1.08 mg/L, 1.77 μg/L, and 5.86 mg/L, respectively. This study also enabled SPR analysis using pseudotyped lentiviral virions instead of isolated S protein, which was of more biological relevance. In the binding competition experiments, soluble heparin, non-anti-coagulant heparin, and a non-anti-coagulant low-molecular-weight heparin (NACH) showed IC50 values of 125 nM, 500 nM, and 25 μM, respectively.

Tiwari et al. used a model of cellular cell-to-cell fusion assay to show that the SARS-CoV-2 S protein-mediated cell-to-cell fusion could arise even in the absence of ACE2 (Tiwari et al., 2020). Further, they demonstrated that the S protein differentially recognized the 3-O-sulfated HS structures generated by the two different isoforms, 3OST-3B and 3OST-5, and that the S2 subunit was critical for cell-to-cell fusion mediated by the S protein−3-O-sulfated HSPG pathway. SPGG, a synthetic, small, and highly sulfated non-sugar compound, was capable of serving as an effective inhibitor of cell-to-cell fusion.

In an elegant study, Clausen et al. demonstrated the dependence of HS on SARS-CoV-2 infection (Clausen et al., 2020). Molecular modeling identified the putative binding surface for oligosaccharides that resided in the RBD of the S protein and were adjacent to, but separate from the ACE2 binding site. Interactions between RBD/ectodomain and HS were proved by affinity-based approaches. A ternary complex of heparin, ACE2, and the S protein was demonstrated by binding of S protein to immobilized heparin-BSA and titrating with biotinylated ACE2, in which case the binding of ACE2 increased in proportion to the amount of S protein bound to the heparin-BSA. Through flow cytometry, the authors proved that HS was essential for the spike ectodomains binding to ACE2 and several human cell lines, while heparin lyases treatment dramatically reduced binding. Similarly, targeting heparin/HS synthesis enzymes including NDST1, HS6ST1, HS6ST2, and B4GALT7 (required for GAG assembly) significantly reduced binding. Consistent with Liu et al. (2020), the authors claimed that the interaction between heparin and the S protein was independent of the anti-coagulant activity. Furthermore, infection of pseudotyped vesicular stomatitis virus (VSV) expressing the full-length S protein and SARS-CoV-2 virus was proven to be dependent on cellular HS. Interestingly, Hep3B cells with inactivated HS6ST1/2 responded differently to VSV and SARS-CoV-2.

Relevant to the previously established dependency of HS on SARS-CoV-2 infection, Martino et al. showed that commensal host bacterial communities capable of modifying HS changed with host age and sex in adult COVID-19 patients. The prevalence of those bacteria and the expression of key microbial glycosidases, which were capable of blocking SARS-CoV-2 S protein binding to human lung adenocarcinoma cells in vitro, was lower in bronchoalveolar lavage fluid (BALF) compared to healthy controls (Martino et al., 2020). Zhang et al. performed a drug repurposing screen and identified Mitoxantrone, an FDA-approved cancer treatment drug that also directly targets HSPGs and inhibits pseudo-coronavirus infection (Zhang et al., 2020a). Several other drugs, Sunitinib BNTX and Latrunculin, which disrupt actin dynamics on the cell surface, were also proven to inhibit SARS-CoV-2 cell entry. The fact that structurally unrelated actin inhibitors all blocked coronavirus entry strongly suggested that the endocytosis of coronavirus required the actin cytoskeleton in addition to ACE2 and HS.

Interactions Between Heparin/HS and Cytokines

Severity of SARS, MERS, and COVID-2019 are associated with the presence of lymphopenia and inflammatory cytokine storm (de Wit et al., 2016; Tan et al., 2020; Zhou et al., 2020). The process of inflammation storm is divided into three steps accompanied by a series of inflammatory responses and recruitment of leukocytes on the infected areas. (1) After initial invasion of virus, macrophages and mast cells immediately release macrophage inflammatory protein 1-α (TNF-α) and interleukin-1 (IL-1) at the site of pathogen adhesiveness in order to activate leukocyte extravasation. (2) Selectins (E-, L-, and P-selectins) on leukocytes interact with endothelial surface-associated HS, which allows leukocyte tethering and rolling along vessel wall (Wang et al., 2002). (3) An array of inflammatory chemokines and cytokines are activated by HS presented on the endothelial surface, which triggers integrin adhesion molecules binding onto leukocytes and subsequent leukocyte extravasation out of the blood vessel (Butcher, 1991; Norgard-Sumnicht and Varki, 1995; Tanaka et al., 1996; Luo et al., 2001). Rich evidence has shown that various inflammatory chemokines and cytokines including microphage inflammatory protein (MIP-1α), RANTES, granulocyte colony-stimulating factor, interferon-γ-inducible protein 10, monocyte chemoattractant protein 1, tumor necrosis factor-α (Kuschert et al., 1999; Huang et al., 2020), as well as IL-2 (Najjam et al., 1998), IL-7 (Borghesi et al., 1999), IL-8 (Spillmann et al., 1998), and IL-10 (Salek-Ardakani et al., 2000) selectively bind to distinct domains of the heterogeneous HS with various affinities and sequence specificities. Gao et al. reported that periodate-oxidized, borohydride-reduced heparin (RO-heparin) could inhibit thioglycollate-induced peritoneal inflammation by preventing neutrophil recruitment dependent on the release of L- and P-selectin (Gao et al., 2005). This is an indication that RO-heparin could attenuate L- and P-selectin-mediated acute inflammation.

Current/Potential Clinical Applications of Heparin/HS in COVID-19

With the evidence described above, it is most likely that disrupting HS–protein interactions by exogenous and competitive heparin/HS mimetics could interfere with virus infection and/or suppress the inflammatory responses. In fact, COVID-19 patients commonly suffer from hyper-coagulopathy and are routinely treated with heparin/LMWH. Significant differences in 28-days mortality were observed in the subgroup of patients with a concentration of D-dimer (>3 μg/ml) higher than sixfold of the normal upper limit, or who had a sepsis-induced coagulopathy (SIC) score ≥4 (40.0 vs 64.2%, P = 0.029) (Shi et al., 2020a). Recent studies have shown that hospitalized patients with severe COVID-19 treated with LMWH or fondaparinux (an ultra-low-molecular-weight heparin) had better prognosis in relation to mortality (Lin et al., 2020; Russo et al., 2020; Tang et al., 2020). It needs to be noted that side effects like heparin-induced thrombocytopenia (HIT) have also been reported in heparin-treated COVID-19 patients (Daviet et al., 2020; Lozano and Franco, 2020). Correct dosage and real-time monitoring of the anti-Xa activity are crucial in heparin/LMWH treatment of COVID-2019 (Duranteau et al., 2018).

Apart from its anti-coagulant effects, a retrospective cohort study found that IL-6 levels were significantly reduced while the percentage of lymphocytes was remarkably increased in the hospitalized COVID-19 treated with LMWH in comparison to the non-LMWH-treated group (Shi et al., 2020b), demonstrating the anti-inflammatory activity of the drug. In addition to the systemic administration, local application of heparin/LMWH through intranasal or inhalation route have also been reported for the treatment of lung diseases and inhalation injury (Yildiz-Pekoz and Ozsoy, 2017; Zielinski et al., 2019). Considering the antiviral activities of heparin/LMWH, along with data suggesting that the nasal epithelium is a portal for initial infection and transmission, Tandon et al. suggested that intranasal administration of UFH may be an effective and safe prophylactic treatment SARS-CoV-2 transmission. Due to the low bioavailability of intranasally administered heparin (Bendstrup et al., 2002), this approach might avoid dangerous side effects or complications with anti-coagulation treatments while potentially still providing a prophylactic or therapeutic benefit (Tandon et al., 2020).

Conclusive Remarks

Despite the well-established anti-coagulant activity and the observed anti-inflammatory effects, the potential anti-SARS-CoV-2 activity of heparin/HS was only recently proposed. It is still controversial regarding the structure specificities of the heparin/HS chains for its interaction with the S protein; however, in vitro experiments and some clinical data have provided promising evidence of heparin/HS (including their mimetics) and heparin/HS-interacting molecules as anti-SARS-CoV-2 drugs. Further elucidation of the heparin/HS–S protein interaction will facilitate the construction of structurally defined oligosaccharide sequences that can be prepared through several methods reported (Roy et al., 2014; Hansen et al., 2015; Baytas and Linhardt, 2020; Zhang et al., 2020b). Non-anti-coagulant heparin, the anti-coagulant activity of which is selectively eliminated, may also be an option to be explored for the treatment of COVID-19 patients without the risk of bleeding complications (Cassinelli et al., 2020; Lindahl and Li, 2020)

Author Contributions

MY and TZ contributed equally to this work by writing the first draft. J-pL contributed to editing the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants from the National Natural Science Foundation of China (31961133004), the Swedish Research Council (2018-02503, 2018-06016), and Beijing Advanced Innovation Center for Soft Matter Science and Engineering.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The handling editor declared a past collaboration with several of the authors TZ, HL, and J-pL.

Acknowledgments

Figures 2, 3 were created with tools obtained from BioRender.com.

References

Arnold, K., Bordoli, L., Kopp, J., and Schwede, T. (2006). The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics 22, 195–201. doi: 10.1093/bioinformatics/bti770

PubMed Abstract | CrossRef Full Text | Google Scholar

Baytas, S. N., and Linhardt, R. J. (2020). Advances in the preparation and synthesis of heparin and related products. Drug Discov Today 25, 2095–2109. doi: 10.1016/j.drudis.2020.09.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Belouzard, S., Chu, V. C., and Whittaker, G. R. (2009). Activation of the SARS coronavirus spike protein via sequential proteolytic cleavage at two distinct sites. Proc. Natl. Acad. Sci. U. S. A. 106, 5871–5876. doi: 10.1073/pnas.0809524106

PubMed Abstract | CrossRef Full Text | Google Scholar

Bendstrup, K. E., Gram, J., and Jensen, J. I. (2002). Effect of inhaled heparin on lung function and coagulation in healthy volunteers. Eur. Respir. J. 19, 606–610. doi: 10.1183/09031936.02.00105202

PubMed Abstract | CrossRef Full Text | Google Scholar

Berman, H. M., Westbrook, J., Feng, Z., Gilliland, G., Bhat, T. N., Weissig, H., et al. (2000). The protein data bank. Nucl. Acids Res. 28, 235–242. doi: 10.1093/nar/28.1.235

CrossRef Full Text | Google Scholar

Bertram, S., Dijkman, R., Habjan, M., Heurich, A., Gierer, S., Glowacka, I., et al. (2013). TMPRSS2 activates the human coronavirus 229E for cathepsin-independent host cell entry and is expressed in viral target cells in the respiratory epithelium. J. Virol. 87, 6150–6160. doi: 10.1128/JVI.03372-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Bertram, S., Glowacka, I., Muller, M. A., Lavender, H., Gnirss, K., Nehlmeier, I., et al. (2011). Cleavage and activation of the severe acute respiratory syndrome coronavirus spike protein by human airway trypsin-like protease. J. Virol. 85, 13363–13372. doi: 10.1128/JVI.05300-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Bitencourt-Ferreira, G., Veit-Acosta, M., and de Azevedo, W. F. (2019). Electrostatic energy in protein-ligand complexes. Methods Mol. Biol. 2053, 67–77. doi: 10.1007/978-1-4939-9752-7_5

PubMed Abstract | CrossRef Full Text | Google Scholar

Bojarski, K. K., Sieradzan, A. K., and Samsonov, S. A. (2019). Molecular dynamics insights into protein-glycosaminoglycan systems from microsecond-scale simulations. Biopolymers 110:e23252. doi: 10.1002/bip.23252

PubMed Abstract | CrossRef Full Text | Google Scholar

Borghesi, L. A., Yamashita, Y., and Kincade, P. W. (1999). Heparan sulfate proteoglycans mediate interleukin-7-dependent B lymphopoiesis. Blood 93, 140–148. doi: 10.1182/blood.V93.1.140

PubMed Abstract | CrossRef Full Text | Google Scholar

Bosch, B. J., Martina, B. E. E., Van Der Zee, R., Lepault, J., Haijema, B. J., Versluis, C., et al. (2004). Severe acute respiratory syndrome coronavirus (SARS-CoV) infection inhibition using spike protein heptad repeat-derived peptides. Proc. Natl. Acad. Sci. U. S. A. 101, 8455–8460. doi: 10.1073/pnas.0400576101

PubMed Abstract | CrossRef Full Text | Google Scholar

Brooks, B. R., Brooks, C. L., Mackerell, A. D. Jr., Nilsson, L., Petrella, R. J., Roux, B., Won, Y., et al. (2009). CHARMM: the biomolecular simulation program. J. Comput. Chem. 30, 1545–1614. doi: 10.1002/jcc.21287

PubMed Abstract | CrossRef Full Text | Google Scholar

Butcher, E. C. (1991). Leukocyte-endothelial cell recognition: three (or more) steps to specificity and diversity. Cell 67, 1033–1036. doi: 10.1016/0092-8674(91)90279-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardin, A. D., and Weintraub, H. J. (1989). Molecular modeling of protein-glycosaminoglycan interactions. Arteriosclerosis 9, 21–32. doi: 10.1161/01.ATV.9.1.21

PubMed Abstract | CrossRef Full Text | Google Scholar

Case, D. A., Cheatham, T. E. 3rd, Darden, T., Gohlke, H., Luo, R., Merz, K. M. Jr., Onufriev, A., et al. (2005). The Amber biomolecular simulation programs. J. Comput. Chem. 26, 1668–1688. doi: 10.1002/jcc.20290

PubMed Abstract | CrossRef Full Text | Google Scholar

Cassinelli, G., Torri, G., and Naggi, A. (2020). Non-anticoagulant heparins as heparanase inhibitors. Adv. Exp. Med. Biol. 1221, 493–522. doi: 10.1007/978-3-030-34521-1_20

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, T., Li, Q., Zhou, Z., Wang, Y., and Bryant, S. H. (2012). Structure-based virtual screening for drug discovery: a problem-centric review. AAPS J. 14, 133–141. doi: 10.1208/s12248-012-9322-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Clausen, T. M., Sandoval, D. R., Spliid, C. B., Pihl, J., Perrett, H. R., Painter, C. D., et al. (2020). SARS-CoV-2 infection depends on cellular heparan sulfate and ACE2. Cell 2020:201616. doi: 10.1101/2020.07.14.201616

PubMed Abstract | CrossRef Full Text | Google Scholar

Coutard, B., Valle, C., de Lamballerie, X., Canard, B., Seidah, N. G., and Decroly, E. (2020). The spike glycoprotein of the new coronavirus 2019-nCoV contains a furin-like cleavage site absent in CoV of the same clade. Antiviral Res. 176:104742. doi: 10.1016/j.antiviral.2020.104742

PubMed Abstract | CrossRef Full Text | Google Scholar

Daviet, F., Guervilly, C., Baldesi, O., Bernard-Guervilly, F., Pilarczyk, E., Genin, A., et al. (2020). Heparin induced thrombocytopenia in severe COVID-19 patients. Circulation 2020:49015. doi: 10.1161/CIRCULATIONAHA.120.049015

PubMed Abstract | CrossRef Full Text | Google Scholar

De Azevedo, W. F. (2010). MolDock applied to structure-based virtual screening. Curr. Drug Targets. 11, 327–334. doi: 10.2174/138945010790711941

PubMed Abstract | CrossRef Full Text | Google Scholar

de Wit, E., van Doremalen, N., Falzarano, D., and Munster, V. J. (2016). SARS and MERS: recent insights into emerging coronaviruses. Nat. Rev. Microbiol. 14, 523–534. doi: 10.1038/nrmicro.2016.81

PubMed Abstract | CrossRef Full Text | Google Scholar

DiGabriele, A. D., Lax, I., Chen, D. I., Svahn, C. M., Jaye, M., Schlessinger, J., et al. (1998). Structure of a heparin-linked biologically active dimer of fibroblast growth factor. Nature 393, 812–817. doi: 10.1038/31741

PubMed Abstract | CrossRef Full Text | Google Scholar

Drosten, C., Gunther, S., Preiser, W., van der Werf, S., Brodt, H. R., Becker, S., et al. (2003). Identification of a novel coronavirus in patients with severe acute respiratory syndrome. N. Engl. J. Med. 348, 1967–1976. doi: 10.1056/NEJMoa030747

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, L., He, Y., Zhou, Y., Liu, S., Zheng, B. J., and Jiang, S. (2009). The spike protein of SARS-CoV–a target for vaccine and therapeutic development. Nat. Rev. Microbiol. 7, 226–236. doi: 10.1038/nrmicro2090

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, L., Yang, Y., Zhou, Y., Lu, L., Li, F., and Jiang, S. (2017). MERS-CoV spike protein: a key target for antivirals. Expert. Opin. Ther. Targets. 21, 131–143. doi: 10.1080/14728222.2017.1271415

CrossRef Full Text | Google Scholar

Duranteau, J., Taccone, F. S., Verhamme, P., Ageno, W., and Force, E. V. G. T. (2018). European guidelines on perioperative venous thromboembolism prophylaxis: intensive care. Eur. J. Anaesthesiol. 35, 142–146. doi: 10.1097/EJA.0000000000000707

CrossRef Full Text | Google Scholar

Elfiky, A. A., Mahdy, S. M., and Elshemey, W. M. (2017). Quantitative structure-activity relationship and molecular docking revealed a potency of anti-hepatitis C virus drugs against human corona viruses. J. Med. Virol. 89, 1040–1047. doi: 10.1002/jmv.24736

CrossRef Full Text | Google Scholar

Fehr, A. R., and Perlman, S. (2015). Coronaviruses: an overview of their replication and pathogenesis. Methods Mol. Biol. 1282, 1–23. doi: 10.1007/978-1-4939-2438-7_1

CrossRef Full Text | Google Scholar

Forsten-Williams, K., Chu, C. L., Fannon, M., Buczek-Thomas, J. A., and Nugent, M. A. (2008). Control of growth factor networks by heparan sulfate proteoglycans. Ann. Biomed. Eng. 36, 2134–2148. doi: 10.1007/s10439-008-9575-z

CrossRef Full Text | Google Scholar

Friedrich, U., Blom, A. M., Dahlback, B., and Villoutreix, B. O. (2001). Structural and energetic characteristics of the heparin-binding site in antithrombotic protein C. J. Biol. Chem. 276, 24122–24128. doi: 10.1074/jbc.M011567200

CrossRef Full Text | Google Scholar

Gallagher, T. M., and Buchmeier, M. J. (2001). Coronavirus spike proteins in viral entry and pathogenesis. Virology 279, 371–374. doi: 10.1006/viro.2000.0757

CrossRef Full Text | Google Scholar

Gao, Y., Li, N., Fei, R., Chen, Z., Zheng, S., and Zeng, X. (2005). P-Selectin-mediated acute inflammation can be blocked by chemically modified heparin, RO-heparin. Mol. Cells 19, 350–355. Available online at: http://www.molcells.org/journal/view.html?year=2005&volume=19&number=3&spage=350

PubMed Abstract | Google Scholar

Gasteiger, J., and Marsili, M. (1980). Iterative partial equalization of orbital electronegativity—a rapid access to atomic charges. Tetrahedron 36, 3219–3228. doi: 10.1016/0040-4020(80)80168-2

CrossRef Full Text | Google Scholar

Gierer, S., Bertram, S., Kaup, F., Wrensch, F., Heurich, A., Kramer-Kuhl, A., et al. (2013). The spike protein of the emerging betacoronavirus EMC uses a novel coronavirus receptor for entry, can be activated by TMPRSS2, and is targeted by neutralizing antibodies. J. Virol. 87, 5502–5511. doi: 10.1128/JVI.00128-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Goldberg, D. (1988). Genetic Algorithms in Search Optimization and Machine Learning. New York, NY: Addison-Wesley Publishing Company, Inc.

Google Scholar

Gomes, P. B., and Dietrich, C. P. (1982). Distribution of heparin and other sulfated glycosaminoglycans in vertebrates. Comp. Biochem. Physiol. B. 73, 857–863. doi: 10.1016/0305-0491(82)90329-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Goodsell, D. S., Morris, G. M., and Olson, A. J. (1996). Automated docking of flexible ligands: applications of AutoDock. J. Mol. Recognit. 9, 1–5. doi: 10.1002/(SICI)1099-1352(199601)9:1<1::AID-JMR241>3.0.CO;2-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Guvench, O., and MacKerell, A. D. (2008). Comparison of protein force fields for molecular dynamics simulations. Methods Mol. Biol. 443, 63–88. doi: 10.1007/978-1-59745-177-2_4

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, Y., and Kral, P. (2020). Computational design of ACE2-based peptide inhibitors of SARS-CoV-2. ACS Nano. 14, 5143–5147. doi: 10.1021/acsnano.0c02857

PubMed Abstract | CrossRef Full Text | Google Scholar

Hansen, S. U., Miller, G. J., Cliff, M. J., Jayson, G. C., and Gardiner, J. M. (2015). Making the longest sugars: a chemical synthesis of heparin-related [4] n oligosaccharides from 16-mer to 40-mer. Chem. Sci. 6, 6158–6164. doi: 10.1039/C5SC02091C

PubMed Abstract | CrossRef Full Text | Google Scholar

Hao, C., Xu, H., Yu, L., and Zhang, L. (2019). Heparin: an essential drug for modern medicine. Prog. Mol. Biol. Transl. Sci. 163, 1–19. doi: 10.1016/bs.pmbts.2019.02.002

PubMed Abstract | CrossRef Full Text | Google Scholar

He, Y., Li, J., Heck, S., Lustigman, S., and Jiang, S. (2006). Antigenic and immunogenic characterization of recombinant baculovirus-expressed severe acute respiratory syndrome coronavirus spike protein: implication for vaccine design. J. Virol. 80, 5757–5767. doi: 10.1128/JVI.00083-06

PubMed Abstract | CrossRef Full Text | Google Scholar

Hess, B., Kutzner, C., van der Spoel, D., and Lindahl, E. (2008). GROMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J. Chem. Theory Comput. 4, 435–447. doi: 10.1021/ct700301q

PubMed Abstract | CrossRef Full Text | Google Scholar

Hileman, R. E., Fromm, J. R., Weiler, J. M., and Linhardt, R. J. (1998). Glycosaminoglycan-protein interactions: definition of consensus sites in glycosaminoglycan binding proteins. Bioessays 20, 156–167. doi: 10.1002/(SICI)1521-1878(199802)20:2<156::AID-BIES8>3.0.CO;2-R

PubMed Abstract | CrossRef Full Text | Google Scholar

Homola, J. (2008). Surface plasmon resonance sensors for detection of chemical and biological species. Chem. Rev. 108, 462–493. doi: 10.1021/cr068107d

PubMed Abstract | CrossRef Full Text | Google Scholar

Horton, M., Su, G., Yi, L., Wang, Z., Xu, Y., Pagadala, V., et al. (2020). Construction of heparan sulfate microarray for investigating the binding of specific saccharide sequences to proteins. Glycobiology 2020:cwaa068. doi: 10.1093/glycob/cwaa068

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, C., Wang, Y., Li, X., Ren, L., Zhao, J., Hu, Y., et al. (2020). Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet 395, 497–506. doi: 10.1016/S0140-6736(20)30183-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Hubbard, A. R., and Jennings, C. A. (1985). Neutralisation of heparan sulphate and low molecular weight heparin by protamine. Thromb Haemost. 53, 86–89. doi: 10.1055/s-0038-1661242

PubMed Abstract | CrossRef Full Text | Google Scholar

Hui, D. S., E, I. A., Madani, T. A., Ntoumi, F., Kock, R., Dar, O., et al. (2020). The continuing 2019-nCoV epidemic threat of novel coronaviruses to global health—the latest 2019 novel coronavirus outbreak in Wuhan, China. Int. J. Infect. Dis. 91, 264–266. doi: 10.1016/j.ijid.2020.01.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Jemth, P., Kreuger, J., Kusche-Gullberg, M., Sturiale, L., Gimenez-Gallego, G., and Lindahl, U. (2002). Biosynthetic oligosaccharide libraries for identification of protein-binding heparan sulfate motifs. Exploring the structural diversity by screening for fibroblast growth factor (FGF)1 and FGF2 binding. J. Biol. Chem. 277, 30567–73. doi: 10.1074/jbc.M203404200

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, G., Willett, P., Glen, R. C., Leach, A. R., and Taylor, R. (1997). Development and validation of a genetic algorithm for flexible docking. J. Mol. Biol. 267, 727–748. doi: 10.1006/jmbi.1996.0897

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, S. Y., Jin, W., Sood, A., Montgomery, D. W., Grant, O. C., Fuster, M. M., et al. (2020). Characterization of heparin and severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) spike glycoprotein binding interactions. Antiviral. Res. 181:104873. doi: 10.1016/j.antiviral.2020.104873

PubMed Abstract | CrossRef Full Text | Google Scholar

Kirkpatrick, S., Gelatt, C. D. Jr., and Vecchi, M. P. (1983). Optimization by simulated annealing. Science 220, 671–680. doi: 10.1126/science.220.4598.671

PubMed Abstract | CrossRef Full Text | Google Scholar

Kubo, H., Yamada, Y. K., and Taguchi, F. (1994). Localization of neutralizing epitopes and the receptor-binding site within the amino-terminal 330 amino acids of the murine coronavirus spike protein. J. Virol. 68, 5403–5410. doi: 10.1128/JVI.68.9.5403-5410.1994

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuhn, J. H., Li, W., Choe, H., and Farzan, M. (2004). Angiotensin-converting enzyme 2: a functional receptor for SARS coronavirus. Cell Mol. Life Sci. 61, 2738–2743. doi: 10.1007/s00018-004-4242-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuschert, G. S., Coulin, F., Power, C. A., Proudfoot, A. E., Hubbard, R. E., Hoogewerf, A. J., et al. (1999). Glycosaminoglycans interact selectively with chemokines and modulate receptor binding and cellular responses. Biochemistry 38, 12959–12968. doi: 10.1021/bi990711d

PubMed Abstract | CrossRef Full Text | Google Scholar

Lang, J., Yang, N., Deng, J., Liu, K., Yang, P., Zhang, G., et al. (2011). Inhibition of SARS pseudovirus cell entry by lactoferrin binding to heparan sulfate proteoglycans. PLoS ONE 6:e23710. doi: 10.1371/journal.pone.0023710

PubMed Abstract | CrossRef Full Text | Google Scholar

Levitt, M. (1992). Accurate modeling of protein conformation by automatic segment matching. J. Mol. Biol. 226, 507–533. doi: 10.1016/0022-2836(92)90964-L

PubMed Abstract | CrossRef Full Text | Google Scholar

Lewis, R. M., and Torczon, V. (2002). A globally convergent augmented lagrangian pattern search algorithm for optimization with general constraints and simple bounds. SIAM J. Optimizat. 12, 1075–1089. doi: 10.1137/S1052623498339727

CrossRef Full Text | Google Scholar

Li, F. (2016). Structure, function, and evolution of coronavirus spike proteins. Annu. Rev. Virol. 3, 237–261. doi: 10.1146/annurev-virology-110615-042301

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J. P., and Kusche-Gullberg, M. (2016). Heparan sulfate: biosynthesis, structure, and function. Int. Rev. Cell Mol. Biol. 325, 215–273. doi: 10.1016/bs.ircmb.2016.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, L., Lu, L., Cao, W., and Li, T. (2020). Hypothesis for potential pathogenesis of SARS-CoV-2 infection-a review of immune changes in patients with viral pneumonia. Emerg. Microbes Infect. 9, 727–732. doi: 10.1080/22221751.2020.1746199

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindahl, U., Backstrom, G., Thunberg, L., and Leder, I. G. (1980). Evidence for a 3-O-sulfated D-glucosamine residue in the antithrombin-binding sequence of heparin. Proc. Natl. Acad. Sci. U. S. A. 77, 6551–6555. doi: 10.1073/pnas.77.11.6551

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindahl, U., and Li, J. P. (2020). Heparin - an old drug with multiple potential targets in Covid-19 therapy. J. Thromb. Haemost. 2020:14898. doi: 10.1111/jth.14898

PubMed Abstract | CrossRef Full Text | Google Scholar

Linhardt, R. J., and Toida, T. (2004). Role of glycosaminoglycans in cellular communication. Acc. Chem. Res. 37, 431–438. doi: 10.1021/ar030138x

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, L., Chopra, P., Li, X., Wolfert, M. A., Tompkins, S. M., and Boons, G. J. (2020). SARS-CoV-2 spike protein binds heparan sulfate in a length- and sequence-dependent manner. bioRxiv. doi: 10.1101/2020.05.10.087288

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, S., Xiao, G., Chen, Y., He, Y., Niu, J., Escalante, C. R., et al. (2004). Interaction between heptad repeat 1 and 2 regions in spike protein of SARS-associated coronavirus: implications for virus fusogenic mechanism and identification of fusion inhibitors. Lancet 363, 938–947. doi: 10.1016/S0140-6736(04)15788-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Lozano, R., and Franco, M. E. (2020). Incidence of heparin-induced thrombocytopenia in patients with 2019 coronavirus disease. Med. Clin. 2020:23. doi: 10.1016/j.medcle.2020.05.023

CrossRef Full Text | Google Scholar

Lu, L., Liu, Q., Zhu, Y., Chan, K. H., Qin, L., Li, Y., et al. (2014). Structure-based discovery of Middle East respiratory syndrome coronavirus fusion inhibitor. Nat. Commun. 5:3067. doi: 10.1038/ncomms4067

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, J., Kato, M., Wang, H., Bernfield, M., and Bischoff, J. (2001). Heparan sulfate and chondroitin sulfate proteoglycans inhibit E-selectin binding to endothelial cells. J. Cell Biochem. 80, 522–531. doi: 10.1002/1097-4644(20010315)80:4<522::AID-JCB1006>3.0.CO;2-H

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, W., Yang, L., and He, L. (2018). Overview of the detection methods for equilibrium dissociation constant KD of drug-receptor interaction. J. Pharm. Anal. 8, 147–152. doi: 10.1016/j.jpha.2018.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

MacKerell, A. D., Bashford, D., Bellott, M., Dunbrack, R. L., Evanseck, J. D., Field, M. J., et al. (1998). All-atom empirical potential for molecular modeling and dynamics studies of proteins. J. Phys. Chem. B 102, 3586–3616. doi: 10.1021/jp973084f

PubMed Abstract | CrossRef Full Text | Google Scholar

Martino, C., Kellman, B. P., Sandoval, D. R., Clausen, T. M., Marotz, C. A., Song, S. J., et al. (2020). Bacterial modification of the host glycosaminoglycan heparan sulfate modulates SARS-CoV-2 infectivity. bioRxiv 2020:238444. doi: 10.1101/2020.08.17.238444

CrossRef Full Text | Google Scholar

Milewska, A., Zarebski, M., Nowak, P., Stozek, K., Potempa, J., and Pyrc, K. (2014). Human coronavirus NL63 utilizes heparan sulfate proteoglycans for attachment to target cells. J. Virol. 88, 13221–13230. doi: 10.1128/JVI.02078-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, B. T., Singh, R. P., Klauda, J. B., Hodoscek, M., Brooks, B. R., and Woodcock, H. L. (2008). CHARMMing: a new, flexible web portal for CHARMM. J. Chem. Inf. Model. 48, 1920–1929. doi: 10.1021/ci800133b

PubMed Abstract | CrossRef Full Text | Google Scholar

Millet, J. K., and Whittaker, G. R. (2014). Host cell entry of Middle East respiratory syndrome coronavirus after two-step, furin-mediated activation of the spike protein. Proc. Natl. Acad. Sci. U. S. A. 111, 15214–15219. doi: 10.1073/pnas.1407087111

PubMed Abstract | CrossRef Full Text | Google Scholar

Morris, G. M., Huey, R., Lindstrom, W., Sanner, M. F., Belew, R. K., Goodsell, D. S., et al. (2009). AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility. J. Comput. Chem. 30, 2785–2791. doi: 10.1002/jcc.21256

PubMed Abstract | CrossRef Full Text | Google Scholar

Morris, G. M., and Lim-Wilby, M. (1999). Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function. J. Comput. Chem. 443, 365–382.

Google Scholar

Mycroft-West, C., Su, D., Elli, S., Li, Y., Guimond, S., Miller, G., et al. (2020a). The 2019 coronavirus (SARS-CoV-2) surface protein (Spike) S1 Receptor Binding Domain undergoes conformational change upon heparin binding. bioRxiv 2020:971093. doi: 10.1101/2020.02.29.971093

CrossRef Full Text | Google Scholar

Mycroft-West, C., Su, D., Pagani, I., Rudd, T., Elli, S., Guimond, S., et al. (2020b). Heparin inhibits cellular invasion by SARS-CoV-2: structural dependence of the interaction of the surface protein (spike) S1 receptor binding domain with heparin. bioRxiv 2020:66761. doi: 10.1101/2020.04.28.066761

CrossRef Full Text | Google Scholar

Najjam, S., Mulloy, B., Theze, J., Gordon, M., Gibbs, R., and Rider, C. C. (1998). Further characterization of the binding of human recombinant interleukin 2 to heparin and identification of putative binding sites. Glycobiology 8, 509–516. doi: 10.1093/glycob/8.5.509

PubMed Abstract | CrossRef Full Text | Google Scholar

Norgard-Sumnicht, K., and Varki, A. (1995). Endothelial heparan sulfate proteoglycans that bind to L-selectin have glucosamine residues with unsubstituted amino groups. J. Biol. Chem. 270, 12012–12024. doi: 10.1074/jbc.270.20.12012

PubMed Abstract | CrossRef Full Text | Google Scholar

Nugent, M. A., Zaia, J., and Spencer, J. L. (2013). Heparan sulfate-protein binding specificity. Biochemistry 78, 726–735. doi: 10.1134/S0006297913070055

PubMed Abstract | CrossRef Full Text | Google Scholar

Olson, S. T., Halvorson, H. R., and Bjork, I. (1991). Quantitative characterization of the thrombin-heparin interaction. Discrimination between specific and nonspecific binding models. J. Biol. Chem. 266, 6342–6352.

PubMed Abstract | Google Scholar

Ou, X., Guan, H., Qin, B., Mu, Z., Wojdyla, J. A., Wang, M., et al. (2017). Crystal structure of the receptor binding domain of the spike glycoprotein of human betacoronavirus HKU1. Nat. Commun. 8:15216. doi: 10.1038/ncomms15216

PubMed Abstract | CrossRef Full Text | Google Scholar

Ou, X., Liu, Y., Lei, X., Li, P., Mi, D., Ren, L., et al. (2020). Characterization of spike glycoprotein of SARS-CoV-2 on virus entry and its immune cross-reactivity with SARS-CoV. Nat. Commun. 11:1620. doi: 10.1038/s41467-020-15562-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Pak, J. E., Sharon, C., Satkunarajah, M., Auperin, T. C., Cameron, C. M., Kelvin, D. J., et al. (2009). Structural insights into immune recognition of the severe acute respiratory syndrome coronavirus S protein receptor binding domain. J. Mol. Biol. 388, 815–823. doi: 10.1016/j.jmb.2009.03.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Pekkarinen, A. I., Longstaff, C., and Jones, B. L. (2007). Kinetics of the inhibition of fusarium serine proteinases by barley (Hordeum vulgare L.) inhibitors. J. Agric. Food Chem. 55, 2736–2742. doi: 10.1021/jf0631777

PubMed Abstract | CrossRef Full Text | Google Scholar

Petrey, D., Xiang, Z., Tang, C. L., Xie, L., Gimpelev, M., Mitros, T., et al. (2003). Using multiple structure alignments, fast model building, and energetic analysis in fold recognition and homology modeling. Proteins 53(Suppl.6), 430–435. doi: 10.1002/prot.10550

PubMed Abstract | CrossRef Full Text | Google Scholar

Phillips, J. C., Braun, R., Wang, W., Gumbart, J., Tajkhorshid, E., Villa, E., et al. (2005). Scalable molecular dynamics with NAMD. J. Comput. Chem. 26, 1781–1802. doi: 10.1002/jcc.20289

PubMed Abstract | CrossRef Full Text | Google Scholar

Przybylski, C., Gonnet, F., Saesen, E., Lortat-Jacob, H., and Daniel, R. (2020). Surface plasmon resonance imaging coupled to on-chip mass spectrometry: a new tool to probe protein-GAG interactions. Anal. Bioanal. Chem. 412, 507–519. doi: 10.1007/s00216-019-02267-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Puvirajesinghe, T. M., Ahmed, Y. A., Powell, A. K., Fernig, D. G., Guimond, S. E., and Turnbull, J. E. (2012). Array-based functional screening of heparin glycans. Chem. Biol. 19, 553–558. doi: 10.1016/j.chembiol.2012.03.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Radkiewicz, J. L., and Brooks, C. L. (2000). Protein dynamics in enzymatic catalysis: exploration of dihydrofolate reductase. J. Am. Chem. Soc. 122, 225–231. doi: 10.1021/ja9913838

CrossRef Full Text | Google Scholar

Raj, V. S., Mou, H., Smits, S. L., Dekkers, D. H., Muller, M. A., Dijkman, R., et al. (2013). Dipeptidyl peptidase 4 is a functional receptor for the emerging human coronavirus-EMC. Nature 495, 251–254. doi: 10.1038/nature12005

PubMed Abstract | CrossRef Full Text | Google Scholar

Rarey, M., Kramer, B., Lengauer, T., and Klebe, G. (1996). A fast flexible docking method using an incremental construction algorithm. J. Mol. Biol. 261, 470–489. doi: 10.1006/jmbi.1996.0477

PubMed Abstract | CrossRef Full Text | Google Scholar

Richard, B., Swanson, R., and Olson, S. T. (2009). The signature 3-O-sulfo group of the anticoagulant heparin sequence is critical for heparin binding to antithrombin but is not required for allosteric activation. J. Biol. Chem. 284, 27054–27064. doi: 10.1074/jbc.M109.029892

PubMed Abstract | CrossRef Full Text | Google Scholar

Roy, S., El Hadri, A., Richard, S., Denis, F., Holte, K., Duffner, J., et al. (2014). Synthesis and biological evaluation of a unique heparin mimetic hexasaccharide for structure-activity relationship studies. J. Med. Chem. 57, 4511–4520. doi: 10.1021/jm4016069

PubMed Abstract | CrossRef Full Text | Google Scholar

Russo, V., Cardillo, G., Viggiano, G. V., Mangiacapra, S., Cavalli, A., Fontanella, A., et al. (2020). Fondaparinux use in patients with COVID-19: a preliminary multicenter real-world experience. J. Cardiovasc. Pharmacol. 76, 369–371. doi: 10.1097/FJC.0000000000000893

PubMed Abstract | CrossRef Full Text | Google Scholar

Salek-Ardakani, S., Arrand, J. R., Shaw, D., and Mackett, M. (2000). Heparin and heparan sulfate bind interleukin-10 and modulate its activity. Blood 96, 1879–1888. doi: 10.1182/blood.V96.5.1879

PubMed Abstract | CrossRef Full Text | Google Scholar

Sali, A., and Blundell, T. L. (1993). Comparative protein modelling by satisfaction of spatial restraints. J. Mol. Biol. 234, 779–815. doi: 10.1006/jmbi.1993.1626

PubMed Abstract | CrossRef Full Text | Google Scholar

Salsbury, F. R. Jr., Crowley, M. F., and Brooks, C. L. (2001). Modeling of the metallo-beta-lactamase from B. fragilis: structural and dynamic effects of inhibitor binding. Proteins 44, 448–459. doi: 10.1002/prot.1110

PubMed Abstract | CrossRef Full Text | Google Scholar

Sapay, N., Cabannes, E., Petitou, M., and Imberty, A. (2011). Molecular modeling of the interaction between heparan sulfate and cellular growth factors: bringing pieces together. Glycobiology 21, 1181–1193. doi: 10.1093/glycob/cwr052

PubMed Abstract | CrossRef Full Text | Google Scholar

Schlick, T. (2010). “Molecular Modeling and Simulation: An Interdisciplinary Guide,” in Molecular Dynamics: Basics (New York, NY: Springer), 429–435.

Google Scholar

Sebaugh, J. L. (2011). Guidelines for accurate EC50/IC50 estimation. Pharm. Stat. 10, 128–134. doi: 10.1002/pst.426

PubMed Abstract | CrossRef Full Text | Google Scholar

Seffer, M. T., Cottam, D., Forni, L. G., and Kielstein, J. T. (2020). Heparin 2.0: a new approach to the infection crisis. Blood Purif. 47, 1–7. doi: 10.1159/000508647

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, C., Tingting, W., Li, J.-P., Sullivan, M. A., Wang, C., Wang, H., et al. (2020a). Comprehensive landscape of heparin therapy for COVID-19. Carbohydrate Polym. 254:117232. doi: 10.1016/j.carbpol.2020.117232

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, C., Wang, C., Wang, H., Yang, C., Cai, F., Zeng, F., et al. (2020b). The potential of low molecular weight heparin to mitigate cytokine storm in severe COVID-19 patients: a retrospective cohort study. Clin. Transl. Sci. 2020:12880. doi: 10.1111/cts.12880

PubMed Abstract | CrossRef Full Text | Google Scholar

Solis, F. J., and Wets, R. J. B. (1981). Minimization by random search techniques. Math. Operat. Res. 6, 19–30. doi: 10.1287/moor.6.1.19

CrossRef Full Text | Google Scholar

Spillmann, D., Witt, D., and Lindahl, U. (1998). Defining the interleukin-8-binding domain of heparan sulfate. J. Biol. Chem. 273, 15487–15493. doi: 10.1074/jbc.273.25.15487

PubMed Abstract | CrossRef Full Text | Google Scholar

Steinbrecher, T., and Labahn, A. (2010). Towards accurate free energy calculations in ligand protein-binding studies. Curr. Med. Chem. 17, 767–785. doi: 10.2174/092986710790514453

PubMed Abstract | CrossRef Full Text | Google Scholar

Sterner, E., Meli, L., Kwon, S. J., Dordick, J. S., and Linhardt, R. J. (2013). FGF-FGFR signaling mediated through glycosaminoglycans in microtiter plate and cell-based microarray platforms. Biochemistry 52, 9009–9019. doi: 10.1021/bi401284r

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, S., Wong, G., Shi, W., Liu, J., Lai, A. C. K., Zhou, J., et al. (2016). Epidemiology, genetic recombination, and pathogenesis of coronaviruses. Trends Microbiol. 24, 490–502. doi: 10.1016/j.tim.2016.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, L., Wang, Q., Zhang, D., Ding, J., Huang, Q., Tang, Y. Q., et al. (2020). Lymphopenia predicts disease severity of COVID-19: a descriptive and predictive study. Signal. Transduct. Target Ther. 5:33. doi: 10.1038/s41392-020-0148-4

CrossRef Full Text | Google Scholar

Tanaka, Y., Kimata, K., Wake, A., Mine, S., Morimoto, I., Yamakawa, N., et al. (1996). Heparan sulfate proteoglycan on leukemic cells is primarily involved in integrin triggering and its mediated adhesion to endothelial cells. J. Exp. Med. 184, 1987–1997. doi: 10.1084/jem.184.5.1987

PubMed Abstract | CrossRef Full Text | Google Scholar

Tandon, R., Sharp, J. S., Zhang, F., Pomin, V. H., Ashpole, N. M., Mitra, D., et al. (2020). Effective inhibition of SARS-CoV-2 entry by heparin and enoxaparin derivatives. bioRxiv 2020:140236. doi: 10.1101/2020.06.08.140236

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, N., Bai, H., Chen, X., Gong, J., Li, D., and Sun, Z. (2020). Anticoagulant treatment is associated with decreased mortality in severe coronavirus disease 2019 patients with coagulopathy. J. Thromb. Haemost. 18, 1094–1099. doi: 10.1111/jth.14817

CrossRef Full Text | Google Scholar

Thompson, L. D., Pantoliano, M. W., and Springer, B. A. (1994). Energetic characterization of the basic fibroblast growth factor-heparin interaction: identification of the heparin binding domain. Biochemistry 33, 3831–3840. doi: 10.1021/bi00179a006

PubMed Abstract | CrossRef Full Text | Google Scholar

Tianji Zhang, R. Z., Yongqin, L.v., Meng, W., Hongmei, L., Tianwei, T., and Jin-Ping, L. (2019). Glycosaminoglycans in biological samples – towards identification of novel biomarkers. TrAC 2019:115732. doi: 10.1016/j.trac.2019.115732

CrossRef Full Text | Google Scholar

Tiwari, V., Tandon, R., Sankaranarayanan, N. V., Beer, J. C., Kohlmeir, E. K., Swanson-Mungerson, M., et al. (2020). Preferential recognition and antagonism of SARS-CoV-2 spike glycoprotein binding to 3- O -sulfated heparan sulfate. bioRxiv 2020:331751. doi: 10.1101/2020.10.08.331751

CrossRef Full Text | Google Scholar

Trott, O., and Olson, A. J. (2010). AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 31, 455–461. doi: 10.1002/jcc.21334

PubMed Abstract | CrossRef Full Text | Google Scholar

Turnbull, J. E., Fernig, D. G., Ke, Y., Wilkinson, M. C., and Gallagher, J. T. (1992). Identification of the basic fibroblast growth factor binding sequence in fibroblast heparan sulfate. J. Biol. Chem. 267, 10337–10341.

PubMed Abstract | Google Scholar

Van Der Spoel, D., Lindahl, E., Hess, B., Groenhof, G., Mark, A. E., and Berendsen, H. J. (2005). GROMACS: fast, flexible, and free. J. Comput. Chem. 26, 1701–1718. doi: 10.1002/jcc.20291

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Brown, J. R., Varki, A., and Esko, J. D. (2002). Heparin's anti-inflammatory effects require glucosamine 6-O-sulfation and are mediated by blockade of L- and P-selectins. J. Clin. Invest. 110, 127–136. doi: 10.1172/JCI0214996

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, N., Shang, J., Jiang, S., and Du, L. (2020). Subunit vaccines against emerging pathogenic human coronaviruses. Front. Microbiol. 11:298. doi: 10.3389/fmicb.2020.00298

PubMed Abstract | CrossRef Full Text | Google Scholar

Werner, T., Morris, M. B., Dastmalchi, S., and Church, W. B. (2012). Structural modelling and dynamics of proteins for insights into drug interactions. Adv. Drug Deliv. Rev. 64, 323–343. doi: 10.1016/j.addr.2011.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Wickramasinghe, I. N., de Vries, R. P., Grone, A., de Haan, C. A., and Verheije, M. H. (2011). Binding of avian coronavirus spike proteins to host factors reflects virus tropism and pathogenicity. J. Virol. 85, 8903–8912. doi: 10.1128/JVI.05112-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, C., Liu, Y., Yang, Y., Zhang, P., Zhong, W., Wang, Y., et al. (2020). Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods. Acta Pharm. Sin. B 10, 766–788. doi: 10.1016/j.apsb.2020.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, S., Zhu, Y., Liu, M., Lan, Q., Xu, W., Wu, Y., et al. (2020). Fusion mechanism of 2019-nCoV and fusion inhibitors targeting HR1 domain in spike protein. Cell Mol. Immunol. 17, 765–767. doi: 10.1038/s41423-020-0374-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, D., and Esko, J. D. (2014). Demystifying heparan sulfate-protein interactions. Annu. Rev. Biochem. 83, 129–157. doi: 10.1146/annurev-biochem-060713-035314

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, R., Zhang, Y., Li, Y., Xia, L., Guo, Y., and Zhou, Q. (2020). Structural basis for the recognition of SARS-CoV-2 by full-length human ACE2. Science 367, 1444–1448. doi: 10.1126/science.abb2762

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., and Chi, L. (2017). Characterization of structural motifs for interactions between glycosaminoglycans and proteins. Carbohydr. Res. 452, 54–63. doi: 10.1016/j.carres.2017.10.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Yildiz-Pekoz, A., and Ozsoy, Y. (2017). Inhaled heparin: therapeutic efficacy and recent formulations. J. Aerosol. Med. Pulm. Drug Deliv. 30, 143–156. doi: 10.1089/jamp.2015.1273

PubMed Abstract | CrossRef Full Text | Google Scholar

Zaki, A. M., van Boheemen, S., Bestebroer, T. M., Osterhaus, A. D., and Fouchier, R. A. (2012). Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia. N. Engl. J. Med. 367, 1814–1820. doi: 10.1056/NEJMoa1211721

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Q., Chen, C. Z., Swaroop, M., Xu, M., Wang, L., Lee, J., et al. (2020a). Targeting heparan sulfate proteoglycan-assisted endocytosis as a COVID-19 therapeutic option. bioRxiv 2020:2202549. doi: 10.1101/2020.07.14.202549

CrossRef Full Text | Google Scholar

Zhang, X., Lin, L., Huang, H., and Linhardt, R. J. (2020b). Chemoenzymatic synthesis of glycosaminoglycans. Acc. Chem. Res. 53, 335–346. doi: 10.1021/acs.accounts.9b00420

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, W., Liu, Y., Tian, D., Wang, C., Wang, S., Cheng, J., et al. (2020). Potential benefits of precise corticosteroids therapy for severe 2019-nCoV pneumonia. Signal Transduct. Target Ther. 5:18. doi: 10.1038/s41392-020-0127-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Zielinski, M., Wroblewski, P., and Kozielski, J. (2019). Is inhaled heparin a viable therapeutic option in inhalation injury? Adv. Respir. Med. 87, 184–188. doi: 10.5603/ARM.2019.0029

PubMed Abstract | CrossRef Full Text | Google Scholar

Zong, C., Venot, A., Li, X., Lu, W., Xiao, W., Wilkes, J. L., et al. (2017). Heparan sulfate microarray reveals that heparan sulfate-protein binding exhibits different ligand requirements. J. Am. Chem. Soc. 139, 9534–9543. doi: 10.1021/jacs.7b01399

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: heparin, heparan sulfate, COVID-19, coronavirus, interaction

Citation: Yu M, Zhang T, Zhang W, Sun Q, Li H and Li J-p (2021) Elucidating the Interactions Between Heparin/Heparan Sulfate and SARS-CoV-2-Related Proteins—An Important Strategy for Developing Novel Therapeutics for the COVID-19 Pandemic. Front. Mol. Biosci. 7:628551. doi: 10.3389/fmolb.2020.628551

Received: 12 November 2020; Accepted: 16 December 2020;
Published: 25 January 2021.

Edited by:

Lianli Chi, Shandong University, China

Reviewed by:

Chao Cai, Ocean University of China, China
Xudong Qu, Shanghai Jiao Tong University, China

Copyright © 2021 Yu, Zhang, Zhang, Sun, Li and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jin-ping Li, amluLXBpbmcubGkmI3gwMDA0MDtpbWJpbS51dS5zZQ==; Hongmei Li, bGlobSYjeDAwMDQwO25pbS5hYy5jbg==

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.