Inhibition of SARS-CoV-2 viral entry upon blocking N- and O-glycan elaboration
Abstract
The Spike protein of SARS-CoV-2, its receptor-binding domain (RBD), and its primary receptor ACE2 are extensively glycosylated. The impact of this post-translational modification on viral entry is yet unestablished. We expressed different glycoforms of the Spike-protein and ACE2 in CRISPR-Cas9 glycoengineered cells, and developed corresponding SARS-CoV-2 pseudovirus. We observed that N- and O-glycans had only minor contribution to Spike-ACE2 binding. However, these carbohydrates played a major role in regulating viral entry. Blocking N-glycan biosynthesis at the oligomannose stage using both genetic approaches and the small molecule kifunensine dramatically reduced viral entry into ACE2 expressing HEK293T cells. Blocking O-glycan elaboration also partially blocked viral entry. Mechanistic studies suggest multiple roles for glycans during viral entry. Among them, inhibition of N-glycan biosynthesis enhanced Spike-protein proteolysis. This could reduce RBD presentation on virus, lowering binding to host ACE2 and decreasing viral entry. Overall, chemical inhibitors of glycosylation may be evaluated for COVID-19.
eLife digest
COVID-19 is an infectious disease caused by the virus SARS-CoV-2. To access the internal machinery necessary for its replication, the virus needs to latch onto and then enter host cells. Such processes rely on specific ‘glycoproteins’ that carry complex sugar molecules (or glycans), and can be found at the surface of both viruses and host cells. In particular, the viral ‘Spike’ glycoprotein can attach to human proteins called ACE2, which coat the cells that line the inside of the lungs, heart, kidney and brain. Yet the roles played by glycans in these processes remains unclear.
To investigate the role of Spike and ACE-2 glycans, Yang et al. designed a form of SARS-CoV-2 that could be handled safely in the laboratory. How these viruses infect human kidney cells that carry ACE2 was then examined, upon modifying the structures of the sugars on the viral Spike protein as well as the host ACE2 receptor. In particular, the sugar structures displayed by the virus were modified either genetically or chemically, using a small molecule that disrupts the formation of the glycans. Similar methods were also applied to modify the glycans of ACE2. Together, these experiments showed that the sugars present on the Spike protein play a minor role in helping the virus stick to human cells.However, they were critical for the virus to fuse and enter the host cells. These findings highlight the important role of Spike protein sugars in SARS-CoV-2 infection, potentially offering new paths to treat COVID-19 and other coronavirus-related illnesses. In particular, molecules designed to interfere with Spike-proteins and the viral entrance into cells could be less specific to SARS-CoV-2 compared to vaccines, allowing treatments to be efficient even if the virus changes.
Introduction
SARS-CoV-2 is the virus causing the global pandemic ‘coronavirus diseases 2019’ (COVID-19). This is a zoonotic beta-coronavirus from bats that causes severe acute respiratory syndrome (Lu et al., 2020). Its effects on human physiology extends well beyond the lung as it unleashes a cytokine storm to alter immune response (Chen et al., 2020) and cause disseminated intravascular coagulation (Lillicrap, 2020). It also exhibits multiorgan tropism that impacts kidney, liver, heart and brain function (Puelles et al., 2020). Among its structural elements, the SARS-CoV-2 Spike protein is critical for viral attachment, fusion and entry into host (Hoffmann et al., 2020a; Li et al., 2003; Lan et al., 2020). The RBD region of Spike enables binding to its primary human cellular receptor, angiotensin-converting enzyme-2 (ACE2), which is ubiquitously expressed on epithelial, endothelial and blood cells (Li et al., 2003; Hamming et al., 2004). A role for heparan sulfates in viral entry has also been proposed (Clausen et al., 2020). Once bound, viral fusion and entry depends on two proteolysis sites, the furin/RRAR-S site located between the Spike S1 and S2 subunits, and a second S2’/SKR-S site (Belouzard et al., 2009). A variety of enzymes including TMPRSS2 (transmembrane protease, serine 2) and cathepsin-L may aid viral entry (Hoffmann et al., 2020b; Matsuyama et al., 2005; Shang et al., 2020). While inhibitors of viral binding, RNA transcription and protease activity are being tested to reduce viral load, this manuscript suggests an orthogonal approach based on glycan engineering.
Both the viral Spike-protein and ACE2 receptor are extensively glycosylated, with a majority of the 22 N-glycosylation sites of Spike and 7 N-glycosylation sites of ACE2 bearing carbohydrates (Walls et al., 2020; Wrapp et al., 2020; Shajahan et al., 2020a; Figure 1). The exact distribution of the oligomannose, hybrid, complex, sialylated, and fucosylated structures in these macromolecules is likely dictated both by the protein structure and host expression system (Shajahan et al., 2020b; Watanabe et al., 2020). O-linked glycans are also reported on both proteins (Shajahan et al., 2020a; Shajahan et al., 2020b). A variety of functions have been proposed for these glycans including immunological shielding (Watanabe et al., 2020), direct regulation of Spike-ACE2 binding (Bernardi et al., 2020), and control of Spike up/down conformation (Casalino et al., 2020). Experimental data supporting these concepts is currently lacking and thus our understanding of the impact of glycosylation on SARS-CoV-2 function is incomplete.
We address the above gaps in knowledge in the current manuscript. Our results suggest that both the O- and N-linked glycans of the SARS-CoV-2 Spike protein, including sialylated glycan epitopes, have a relatively minor role in regulating direct Spike-ACE2 binding interactions. However, these glycans control the rates of viral entry into HEK293T cells expressing human ACE2. Here, blocking N-glycosylation on SARS-CoV-2 pseudovirus at the high-mannose stage using CRISPR-Cas9 and also a small molecule inhibitor (kifunensine) resulted in extensive cleavage/shedding of the viral Spike protein at the time of production due to enhanced proteolysis at the S1-S2 interface. Our data also suggest that glycans may contribute to additional aspects of Spike proteolysis during viral entry. Due to these effects, viral entry into human ACE2 expressing cells was reduced by >95% when virus was produced in the CRISPR-Cas9 MGAT-1 knockout cells lacking complex N-glycans, and by 85–90% upon production in the presence of kifunensine. These observations support the need to screen chemical inhibitors of glycosylation for the treatment of COVID-19.
Results
Modest role for ACE2 sialic acids during Spike-protein molecular recognition
In order to study the impact of N- and O-linked glycosylation on SARS-CoV-2 function, we engineered cells over-expressing full-length Spike-protein and ACE2 (Figure 2A). We also histidine-tag purified soluble, dimeric Fc-his fusion proteins for the Spike S1 region, RBD and extracellular portion of ACE2. The molecules were extensively glycosylated with their apparent molecular mass being greater than their theoretical mass based on peptide alone. Thus, RBD-Fc was ~60 kDa instead of a theoretical mass of 51 kDa, S1-Fc was ~150 kDa rather than 101 kDa, and ACE2-Fc was ~140 kDa instead of 110 kDa (Figure 2B). RBD-Fc and S1-Fc readily bound to HEK293T cells upon human ACE2 over-expression, confirming that they are functionally active (Figure 2C). ACE2-Fc also specifically bound Spike-protein expressed on 293Ts. Baseline RBD-Fc and S1-Fc binding to Spike expressing 293 T cells (i.e. 293 T/S) was even lower than that of wild-type 293T (red line, Figure 2C). This suggests expression of low amounts of Spike binding proteins on WT-293Ts. These basal Spike receptors may engage cell-surface Spike on 293 T/S cell (via cis-interactions) causing RBD-Fc and S1-Fc binding to be lower in 293 T/S cells compared to WT-293Ts.
Sialidase treatment of Spike-protein expressed on 293Ts did not affect ACE2-Fc binding (Figure 2D). Sialidase treatment of ACE2 expressed on 293Ts, however, increased RBD-Fc and S1-Fc binding by 26% and 56% respectively. Control lectin binding studies confirmed the high activity of the sialidase enzyme used in this study (Figure 2—figure supplement 1A). The relatively small effect of the sialidase in these studies was not due to the absence of sialic acid on either S1-Fc or ACE2-Fc. In this regard, independent glycoprotemics mass spectrometry (MS) data showed that ~80–100% of the N-glycans expressed at specific sites of ACE2-Fc and ~40–100% of the S1-Fc glycopeptides expressed complex-type glycans (K.C., et al., manuscript in preparation). Up to ~60% of the antennae on these complex structures on ACE2 and ~20% of the structures on S1-Fc were terminated by sialic acid.
Neither sialidase treatment of Spike nor ACE2 markedly impacted viral entry
To complement the above binding studies, we measured SARS-CoV-2 pseudovirus entry into stable 293T/ACE2 cells. This assay provides an aggregate measure of both molecular binding and viral fusion/entry in the context of the physiological Spike trimeric configuration (Tai et al., 2020). Here, the control VSVG (Vesicular stomatitis virus G-protein) pseudotyped virus displayed broad tropism both for wild-type HEK293T and stable 293T/ACE2 (Figure 2E, Figure 2—figure supplement 1B). Wild-type SARS-CoV-2 Spike-protein virus (‘Spike-WT’) only entered 293T/ACE2 cells confirming the strict dependence on cell-surface ACE2 receptor. A ‘Spike-mutant’ pseudovirus, where the ‘RRAR’ furin-site was swapped with an ‘SRAS’ sequence, also only entered 293T/ACE2 (22). This mutation results in reduced proteolysis at the S1-S2 interface, as discussed later, likely due to the presence of a single arginine site at the S1-S2 interface rather than the polybasic (‘RRAR’) sequence. To study the role of sialic acid on viral entry, we titered the Spike-WT and Spike-mutant virus using p24 ELISA (Figure 2—figure supplement 2). The viruses were then treated with a pan-sialidase from Arthrobacter ureafaciens, and the same amount of Spike-WT and Spike-mutant particles were applied to stable 293T/ACE2 cells (Figure 2F, Figure 2—figure supplement 1C). Here, Spike-mutant pseudotyped virus was ~5–10 times more effective at stimulating DsRed-reporter expression compared to Spike-WT. Thus, the efficiency of cleavage at the S1-S2 interface can fine-tune SARS-CoV-2 infectivity (Hoffmann et al., 2020b; Shang et al., 2020; Walls et al., 2020). Sialidase treatment of virus did not impact viral entry, consistent with the notion that Spike sialic acid does not regulate molecular recognition. Additionally, sialidase treatment of 293T/ACE2 did not alter either VSVG or Spike-WT pseudovirus entry (Figure 2G, Figure 2—figure supplement 1D). A partial increase in Spike-mutant viral entry was measured upon sialidase treatment (p=0.08), suggesting a minor inhibitory role for ACE2 sialic acid when viral infectivity is high. Glycoproteomics analysis of purified pseudovirus produced in HEK293T cells showed that ~45–100% of the viral N-glycans were complex-type, with up to ~55% of their antennae being terminated by sialic acid (K.C., et al. manuscript in preparation). Thus, the pseudovirus used in this study was processed by glycoenzymes that are typically found in the cellular medial and trans-Golgi compartments. Together, the binding and pseudovirus data suggest that ACE2 sialic acids may modestly shield virus/Spike binding in some biological contexts and perhaps some cell systems.
Blocking Spike-protein N-glycan biosynthesis at the high-mannose stage partially reduced ACE2 binding
To determine if other aspects of N- and O-linked glycosylation may impact SARS-CoV-2 function, we applied CRISPR-Cas9 technology to develop isogenic HEK293T clones that lack either the core-1 O-glycan forming galactosyltransferase C1GalT1 or the N-glycan branching β1,2GlcNAc-transferase MGAT1 (Figure 3A; Stolfa et al., 2016; Chugh et al., 2018). Here, knocking out C1GalT1 blocks O-glycan biosynthesis at the GalNAcα-Ser/Thr (+/- sialic acid) stage as it is not elaborated to form the core-1 structure (Galβ1,3GalNAcα). These cells are termed ‘[O]-293T’. Knocking out MGAT1 blocks N-glycan biosynthesis at the Man-5 stage, and the resulting clone is called ‘[N]-293T’. Sanger sequencing confirmed that all three C1GalT1 alleles of [O]-293T contained indels (Figure 3B). Frameshift mutation was also noted at the single MGAT1 allele in [N]-293T. As anticipated, VVA-lectin (Vicia villosa agglutinin) binding to GalNAcα was dramatically augmented in the [O]-293Ts (Figure 3C). Complex glycans were absent in [N]-293Ts as they did not engage PHA-L (lectin from Phaseolus vulgaris). Additional data with a panel of fluorescent lectins further confirm specific blockade of N-glycan complex structures and lactosamine chains in [N]-293T, and inhibition of O-glycan extension in the [O]-293 T cells (Figure 3—figure supplement 1). Together, the findings confirm the absence of O- and N-glycan elaboration in the [O]-293T/C1GalT1-KO and [N]-293T/MGAT1-KO cells, respectively.
Wild-type 293Ts and the glycoenzyme knockouts were transiently transfected to express Spike and ACE2 for binding studies. Both full-length proteins expressed well and at equal levels in the different cell systems (Figure 3—figure supplement 2A). In order to quantify the effect of glycosylation on ACE2-Spike binding, ACE2-Fc was applied to Spike bearing cells, either 293T or the O/N-glycan knockouts (Figure 3D). Similarly, either S1-Fc (Figure 3E) or RBD-Fc (Figure 3F) were applied to ACE2 bearing 293T or the N-/O-glycan knockouts. These Fc-proteins were applied at sub-saturation concentrations in the cytometry binding studies in order to quantitatively evaluate differences in receptor-ligand-binding. Here, we observed a 40% decrease in ACE2-Fc binding to cell-surface Spike expressed on [N]-293Ts, and a slight 20% increase when Spike was expressed on the [O]-293Ts (Figure 3D). Glycan engineering of ACE2 in the [O]- and [N]-293Ts did not impact either S1-Fc (Figure 3E) or RBD-Fc binding (Figure 3F). Here, RBD-Fc bound ~10 times more avidly to ACE2 compared to S1-Fc. Together, the data suggest a modest role for Spike-protein glycosylation on direct Spike-ACE2 binding. Our observation that the binding of RBD for ACE2 was substantially higher than S1-ACE2 interactions is consistent with a previous report (Shang et al., 2020). This highlights the need to carefully consider RBD presentation/conformation in the context of the full protein when quantifying molecular affinity to ACE2.
Blocking N-glycan elaboration on Spike abrogated viral entry
We evaluated the impact of ACE2 and Spike glycosylation on viral entry. To determine the impact of ACE2 glycosylation, we transiently expressed this protein on wild-type 293Ts, [O]-293Ts and [N]-293Ts. The entry of the pseudovirus (VSVG, Spike-WT and Spike-mutant) to these three cell types was then evaluated (Figure 3G, Figure 3—figure supplement 2B). Here, the pattern of viral entry was similar in all cell systems. Lower entry was measured for Spike-WT and Spike-mutant pseudovirus infection of [O]-293T/ACE2 cells, but this was also noted for the control VSVG-virus. Thus, this reduced infectivity in the O-glycan knockout is not exclusively due to ACE2 glycans. We conclude that ACE2 glycans may not regulate viral entry.
In contrast to ACE2, glycan engineering of viral Spike protein dramatically affected viral entry (Figure 4). To study this, VSVG, Spike-WT and Spike-mutant pseudotyped virus were generated in wild-type 293T, [O]-293T and [N]-293T (Figure 4A). All nine viruses were used to infect stable 293T/ACE2 cells. When produced in [O]-293T, Spike-WT and Spike-mutant pseudovirus resulted in a 70–85% loss in DsRed fluorescence (Figure 4B), and a partial reduction (35–70%) in the fraction of DsRed positive cells (Figure 4C). In stark contrast, knocking out N-glycans in both pseudoviruses abrogated viral infection into 293T/ACE2 cells by >95% (Figure 4B and C). Dose studies confirmed loss of viral function upon inhibiting N-glycan elaboration over a range of viral titers (Figure 4D). This was observed for the Spike-mutant and also the Spike-WT pseudovirus (data not shown). Western blot with anti-S2 pAb (polyclonal antibody) suggest that both the Spike-WT and Spike-mutant virus expressed comparable levels of Spike protein, regardless of the glycosylation mutation (Figure 4E). The molecular mass of Spike expressed in [N]-293T was lower than that of other virus consistent with the extensive N-glycosylation of this protein. Additionally, the data suggest more extensive cleavage of Spike-WT compared to Spike-mutant, consistent with the presence of the proteolysis-sensitive polybasic furin-site in Spike-WT. The Spike pseudovirus produced in [N]-293Ts also appeared to have lower intact Spike protein. To examine this more closely as the anti-S2 pAb may have preference for binding the S2-domain over full Spike, we measured the FLAG-epitope at the C-terminus of Spike-mutant using an anti-FLAG mAb (monoclonal antibody). This may provide more uniform recognition of both the full Spike and the S2-subunit (Figure 4F). This blot suggests that N-glycan loss may trigger precipitous (95%, based on densitometry) loss of intact Spike protein due to enhanced cleavage at the S1-S2 interface. Consistent with the viral entry assay, we also noted greater proteolysis of virus produced in [O]-293Ts (52%) compared to that produced in wild-type 293Ts (32%). Together, the data suggest that both Spike N-glycans, and possibly also O-glycans, may play a role in regulating SARS-CoV-2 Spike-protein stability.
Small molecule inhibitors of N-linked glycosylation drastically reduced viral entry
Remarkably, the same observations as seen in the viral N-glycan knockouts could be reproduced upon producing viral particles in cells cultured with the potent mannosidase-I alkaloid inhibitor kifunensine (KI ~25–125 nM, Elbein et al., 1990; Figure 5A). This water-soluble inhibitor reduced the formation of complex N-glycans on cultured cells within 24–48 hr with no change in cell viability (Figure 5—figure supplement 1A). Virus produced in culture with 15 µM kifunensine displayed ~85–90% reduction in infection, as measured using DsRed-reporter in microscopy (Figure 5B, Figure 5—figure supplement 1B) and cytometry assays (Figure 5C and D). The reduction in Spike molecular mass upon addition of this small molecule was less apparent compared to that in [N]-293Ts (Figure 5E), since blockade of Mannosidase-I preferentially gives rise to Man7-9 glycans while the MGAT1/[N]-293Ts knockouts produce Man-5. Here also, the Spike-protein proteolysis was more extensive upon kifunensine addition though it was not complete, as observed in the anti-FLAG blot.
In order to determine if the enhanced S1-S2 site proteolysis of Spike protein upon N-glycan inhibition is the exclusive reason for reduced viral entry, we produced virus with a Spike variant (‘Spike-delta’) that contained the C-terminal FLAG-epitope but that lacked the furin-site (Figure 5F). The substitution of the ‘RRAR’ sequence with a single Alanine (‘A’) resulted in viral Spike that was not extensively cleaved both in the absence and presence of kifunensine. Robust viral entry into 293T/ACE2 cells was observed using this furin resistant virus, and this too could be partially blocked by kifunensine. Thus, in addition to proteolysis at the S1-S2 site, N-glycans may have additional roles during viral entry.
Discussion
This manuscript undertook a series of studies in order to comprehensively describe the role of O- and N-linked glycans on SARS-CoV-2 Spike binding to human ACE2 (Figure 6A), and related viral entry (Figure 6B). It demonstrates only a minor role for carbohydrates in regulating Spike-ACE2 direct binding. In this regard, we observed some enhancement in Spike binding and viral entry upon treating 293T/ACE2 with a pan-sialidase. Based on the published crystal structure and molecular modeling, the critical glycans regulating this process are likely located at Asn-90 or Asn-322, proximal to the Spike binding interface (Lan et al., 2020; Figure 1). The effect was nevertheless small compared to that reported for other sialic acid dependent viruses like influenza (Stencel-Baerenwald et al., 2014), reovirus, parovirus (Löfling et al., 2013) and coronavirus like the Middle-East Respiratory Syndrome virus/MERS (Li et al., 2017).
The studies performed with CRISPR-Cas9 cell lines containing disrupted O- and N-linked glycans also support the concept that glycosylation is not a critical regulator of ACE2-Spike binding. Here, blocking O- and N-glycan biosynthesis on ACE2 did not impact either viral entry or Spike-Fc binding. Thus, the functional change observed upon enzymatically removing sialic acid on ACE2 is likely a steric effect, as it can be offset by other modifications. Blocking elaboration of O- and N-glycans on Spike, also, only had a partial effect on ACE2-Fc binding. This observation is consistent with previous cryo-EM (Walls et al., 2020; Wrapp et al., 2020) and molecular dynamics simulation results (Grant et al., 2020), which show that the active/‘up’ form of Spike RBD is largely exposed with few glycans at the binding interface. It is now well appreciated that this lack of glycan shielding, makes the exposed RBD a prime target for vaccine development.
In contrast to the modest effect of glycosylation on receptor-ligand-binding, we observed a much more profound role for glycans in regulating viral entry in part due to their impact on Spike proteolysis at the S1-S2 interface. In this regard, we observed that blocking N-glycan elaboration using both genetic approaches (i.e. MGAT1-KO/[N]-293T) and the natural product mannosidase-I inhibitor, kifunensine, dramatically reduced viral entry into 293T/ACE2 cells. This was observed for two virus types: Spike-WT and Spike-mut. In the case of Spike-mut, enzymatic and small molecule inhibition resulted in enhanced proteolysis of Spike-protein during pseudovirus production. These data suggest an inverse relation between furin cleavage and viral entry into 293T/ACE2 cells. Consistent with this, Spike-mut displayed reduced S1-S2 proteolysis compared to Spike-WT, and higher viral infectivity. Others have also noted a reduction in viral entry, upon increasing furin cleavage potential by addition of more basic residues (‘RRRKR’ mutation) at the S1-S2 interface (Hoffmann et al., 2020b). Finally, a recent study observed that the introduction of the D614G mutation in Spike reduced furin cleavage, and this resulted in enhanced viral infectivity (Korber et al., 2020). SARS-CoV-2 containing the D614G mutation is currently the dominant strain worldwide (>80% prevalence, Zhang et al., 2020), and the possibility that this is due to altered furin cleavage potential is hotly debated. Overall, our studies with MGAT1-KO, C1GalT1-KO and kifunensine suggest that both N- and O-glycans may modulate the rates of proteolysis at the S1-S2 interface, thus affecting viral entry (Figure 6C). While some S1 may be associated with Spike even after partial cleavage at the S1-S2 interface (Belouzard et al., 2009), this may potentiate a quantitative reduction in RBD presentation on viral surface, resulting in reduced ACE2 binding and lower viral entry. Indeed several N-glycans are located proximal to the furin-site: N603 and N657, and also N61 (Watanabe et al., 2020). O-glycans have also been hypothesized to exist in this region at S673, T678 and S686 (Andersen et al., 2020), and some of these have been experimentally verified (Sanda et al., 2020). How this site-specific glycosylation affects SARS-CoV-2 entry, remains to be determined.
In addition to its role in regulating RBD presentation, our data suggest that N-glycans may also have other roles in regulating viral entry. In this regard, we noted that even the robust entry of Spike-delta pseudovirus into 293T/ACE2 cells could be partially inhibited by kifunensine (Figure 6C), and this is independent of S1-S2 proteolysis. Others have noted that the effect of ‘Spike-delta’ mutation may be cell type specific in that this mutation results in lower infectivity in some cell types like Calu-3 which rely on the enzyme TMPRSS2 for viral fusion, while this mutation does not affect entry into other cell types like VeroE6 which lacks TMPRSS2 (Hoffmann et al., 2020b). Based on these observations, we speculate that additional glycosylation sites proximal to S2’ may also modulate viral entry, perhaps because some proteases like TMPRSS2 have both ligand-binding and protease activity. The N-glycans proximal to S2’ include N801, and to a lesser extent N616 and N657. Additional studies with more cell types and mutant virus are necessary in order to fully elucidate the molecular mechanism by which O- and N-glycans regulate host-cell specific viral entry response.
In addition to basic science understanding, our findings have translational potential as it suggests that both O- and N-glycan truncation may be used to reduce viral entry. In this regard, it would be attractive to test additional N-glycosylation inhibitors, like Swainsonine, which has a demonstrated safety profile in humans (Goss et al., 1997). Additional potential inhibitors that may modulate viral entry include α-mannosidase inhibitors (deoxymannojirimycin, mannostatin A), α-glucosidase inhibitors (N-butyl deoxynojirimycin, N-nonyl-deoxynojirimycin, castanospermine, celgosivir) (Clarke et al., 2020), and finally compounds like SNAP (Del Solar et al., 2020; Wang et al., 2018) and 4F-GalNAc (Marathe et al., 2010) which block various aspects of N- and O-glycan biosynthesis. These potentially represent off-the-shelf drugs and compounds that may be repurposed to reduce viral load and ameliorate SARS-CoV-2 related respiratory symptoms. Studies are currently underway in our laboratory to test these concepts, and extend the assay to authentic SARS-CoV-2 strains. Overall, while SARS-CoV-2 has developed a natural ability to enhance infectivity, we propose that glycoengineering may provide a strategy to take advantage of these evolutionary features to regulate viral entry and reduce disease transmission.
Materials and methods
Materials
Human embryonic kidney 293 T cells (‘293T’) and Lenti-X 293 T cells were purchased from Clontech/Takara Bio (Mountain View, CA). Stable HEK293T/ACE2 cells were kindly provided by Michael Farzan (Scripps Research, Jupiter, FL). All cells were mycoplasma negative. These cell types were maintained in culture using Dulbecco’s Modified Eagle Medium (DMEM) containing 10% fetal bovine serum, 1% Pen-Strep and 1% GlutaMAX supplement. All lectins were from Vector labs. These include VVA (Vicia Villosa lectin, binds GalNAcα on O-glycans), SNA (Sambucus nigra lectin, binds primarily α2,6 sialic acid), ECL (Erythrina cristagalli lectin, binds desialylated lactosamine Galβ1,4GlcNAc chains) and PHA-L (Phaseolus vulgaris Leucoagglutinin, binds Galβ1,4GlcNAcβ1,6(GlcNAcβ1,2Manα1,3)Manα1,3 on complex N-glycans). Goat anti-human ACE2 polyclonal antibody (AF933) and mouse anti-human ACE2 mAb MAB9332 were available from R and D Systems (Minneapolis, MA). Rabbit anti-SARS-CoV-2 RBD (40592-T62) and anti-S2 (40590-T62) pAbs were from Sino Biologicals (Beijing, China). Anti-FLAG clone L5 was from Biolegend (San Diego, CA). All secondary antibodies (Abs) were from Jackson ImmunoResearch (West Grove, PA). When necessary, lectins and Abs were labelled by addition of 25-fold molar excess succinimidyl ester coupled Alexa-405, Alexa-488, Alexa-555 or Alexa-647 dye (Fluoroprobes, Scottsdale, AZ) to protein suspended in phosphate buffered saline (PBS, pH 7.4) for 1 hr at room temperature (RT). Following this, the reaction was quenched with 1/10th volume 1 M Tris, and unreacted Alexa-dye was removed using 7 kDa molecular mass cutoff Zeba desalting spin columns (Thermo). Unless mentioned otherwise, all other reagents were from either Thermo-Fisher or Sigma Chemicals.
Molecular dynamics simulations
Request a detailed protocolThe ACE2-RBD complex structure (6LZG, Wang et al., 2020) was superimposed to the trimeric spike-protein state (6VYB Walls et al., 2020), in which one of the RBD domains is ‘open’ so that ACE2 is associated with the RBD domain in the open conformation. Missing residues in the structure were modeled using PyRosetta (Chaudhury et al., 2010). N-glycans were then added to the structure using Glycan Reader within CHARMM-GUI (Jo et al., 2008; Jo et al., 2011). The identity of the glycans was assigned based on a published work (Watanabe et al., 2020). The final structure contained glycan modifications at 54 spike Asn and 6 ACE2 Asn.
All simulations were performed using NAMD2 (Phillips et al., 2005) using the CHARMM36 parameters (Vanommeslaeghe et al., 2010). The complex was first placed in a water box of sufficient size to ensure a minimum 12 Å separation from the wall. Net charge in the molecule was neutralized by adding 841 K+ and 743 Cl- ions, resulting in 867,166 atoms total. Protein and glycan atoms were first fixed while water (TIP3) molecules were energy minimized for 10,000 steps. Next, the entire system was relaxed for another 10,000 steps, and the system temperature was gradually increased to 310 K in increments of 1 K with 120 fs of equilibration at each temperature. A time step of 2 fs was used. Constant temperature and pressure were maintained using the Langevin framework (Martyna et al., 1994). Long range electrostatic interactions were treated using the particle mesh Ewald method (Darden et al., 1993). The simulation was continued for 26.5 ns. Intermediate structures were saved every 10 ps.
Molecular biology
Request a detailed protocolA plasmid containing the Spike protein (2019-nCov_pcDNA3.1(+)-P2A-eGFP [v1]) was kindly provided by Dr. Haisheng Yu (Institute of Laboratory Animal Science, Peking Union Medical College). To create plasmids [v3] and [v4], the full S1 domain coding region (M1 to R682) or RBD coding region (R319 to F541) was PCR amplified and cloned into the NheI/XbaI or AgeI/XbaI sites, respectively of pCSCG-19FcHisP2AdTom (Lo et al., 2013). Human ACE2 plasmid was obtained from Addgene (#1786, Li et al., 2003). The entire ACE2 coding region was amplified in two segments to silence an EcoRI cutting site in the coding region and joined by overlap extension PCR. The resulting PCR fragment was cloned into the EcoRI/NheI site of pKLV2-CMVpuroBFP vector to obtain [v2]. Additionally, [v5] was created by PCR amplifying the region encoding ACE2 extracellular domain (Q18 to S740) and cloning into the AgeI/XbaI sites of pCSCG-19FcHisP2AdTom plasmid. Plasmids [v2]-[v5] will be distributed via Addgene.
CRISPR-Cas9 knockout cell lines
Request a detailed protocolHEK293T CRISPR-Cas9 cells were generated by transfecting wild-type 293 T cells with pX330-U6-Chimeric_BB-CBh-hSpCas9 vector (Addgene, Plasmid# 42230) carrying single-guide RNA (sgRNA) targeting either the Core-1 β3Gal-T (C1GalT1, target site: GCAGATTCTAGCCAACATAA) or β1,2 GlcNAc-transferase (MGAT1, GTGGGGCGCTATCCTCTTTGTGG). While knocking out the former gene results in cells expressing O-glycans truncated at the Tn-antigen stage (GalNAcα±sialic acid), the latter prevents the N-glycan processing beyond the oligomannose/Man5 stage (Stolfa et al., 2016). Isogenic clones of both cell lines were obtained by single-cell flow cytometry sorting using fluorescently conjugated VVA and PHA-L for selection. Knockouts were confirmed by Sanger sequencing of genomic DNA. For convenience, cells with disrupted O- and N-glycoprotein processing are termed ‘[O]-293T’ and ‘[N]-293T’, respectively.
Recombinant protein expression and purification
Request a detailed protocolAll Fc-his tag proteins (vectors [v3], [v4] and [v5]) were expressed by transient transfection of HEK 293 T cells using the calcium phosphate method (Buffone et al., 2013), in 3–4 150 mm cell culture Petri dishes per protein. Protein secreted into serum-free media were passed through a 1 mL HisTrap FF column (Cytiva, Marlborough, MA) at 1 mL/min. Following extensive washing with 20 mM sodium phosphate buffer containing 0.5 M NaCl and 10 mM imidazole (pH7.2), the bound protein was eluted by collecting 0.5 mL fractions upon step increasing imidazole concentration in the above buffer to 100 mM (8 mL) followed by 500 mM (8 mL).
Protein concentration in each fraction was determined using a flow cytometry bead assay. To this end, briefly, 5 µm Polybead carboxylate microspheres (Polysciences, Warrington, PA) were covalently coupled with goat anti-human pAb using 1-Ethyl-3-(3-dimethylaminopropyl)carbodiimide/EDC chemistry (Marathe et al., 2008). 1 μL of each fraction collected during HisTRAP elution was incubated with these beads for 10 min at RT, before washing and addition of 1:200 diluted FITC conjugated secondary Ab for an additional 10 min. Bead-bound Fc-protein was then quantified using flow cytometry. All fractions containing the Fc-fusion proteins were then pooled and spin concentrated to ~150 μL volume. Buffer was exchanged to HEPES buffer (30 mM HEPES, 110 mM NaCl, 10 mM KCl, 10 mM glucose, 2 mM MgCl2) using 7 kDa Zeba desalting columns. Final protein concentration was determined in ELISA format by using a calibrant Fc-protein of known concentration that was verified to be >95% pure based on silver stain analysis. Here, a common anti-human HRP conjugated Ab was used to quantify both the Fc- proteins described in this manuscript and the calibrant Fc in the same assay.
Identity of expressed protein and also viral Spike was determined using western blotting with anti-Fc (Jackson), anti-RBD (Sino Biologicals), anti-S2 (Sino Biologicals) and anti-ACE2 (R and D Systems) pAbs. Anti-FLAG mAb L5 was also used for western blotting of Spike-mutant virus as it carried a C-terminal FLAG tag.
Flow cytometry
Request a detailed protocolStandard flow cytometry methods were applied in some instances to measure cell-surface protein expression and lectin binding, using directly conjugated fluorescent reagents (Stolfa et al., 2016). Here the lectins/Abs were added to cells for 15 min, prior to performing a quick wash followed by cytometry analysis. These data are presented as Geometric Mean Fluorescence Intensity (‘GMFI’).
Spike-protein-ACE2 binding studies
Request a detailed protocolFull-length Spike-protein and human ACE2 were expressed in 293 T cells (wild-type, [O]-293T and [N]-293T) by transient transfection using the calcium phosphate method (Buffone et al., 2013). These cells are called ‘293 T/S’ or ‘293T/ACE2’ when expressed on wild-type HEK293T cells, and ‘[O]-293 T/S’, ‘[O]-293T/ACE2’, ‘[N]-293 T/S’ or ‘[N]-293T/ACE2’ when expressed in the glycosylation-KO mutants. Here, just prior to functional studies, cells were released from the culture substrate 48–72 hr after transfection using PBS containing 5 mM EDTA, washed and re-suspended in HEPES buffer containing 1.5 mM CaCl2. During the binding studies, 0–4 μg/mL S1-Fc, RBD-Fc or ACE2-Fc were added to 20 × 106 cells/mL at 4°C for 15 min. Following this, 1:200 diluted Alexa-647 (Ax647) conjugated goat anti-human-Fc (Fab')2 Ab was added for an additional 10 min. The cells were then washed, resuspended and read immediately using a 4-laser BD Fortessa flow cytometer. Transfected cells were gated based on EGFP expression for Fc-protein binding studies performed with 293T/S-protein cells, and BFP expression for studies using 293T/ACE2 cells. cell-surface expression of ACE2 and Spike on transiently transfected cells was also monitored on experimental day using flow cytometry. These studies used Alexa-647 conjugated anti-ACE2 (MAB9332, R and D systems) and anti-RBD (Sino Biologicals) antibodies. In some cases, binding studies were also performed with stable 293T/ACE2 cells, in which case fluorescence gating was not necessary.
In some assays, cells were treated with sialidase (200 U/mL Arthrobacter ureafaciens α2–3,6,8,9-Neuraminidase, New England BioLabs) for 1 hr at 37°C, and washed using HEPES buffer prior to the binding assay. Sialidase activity was verified based on SNA and ECL staining. In this case, control cells were treated identically, except that sialidase was withheld. SNA was directly conjugated with Alexa-555 and ECL was conjugated with Alexa-647 in order to simplify the implementation of dual color cytometry measurements.
Pseudovirus production and characterization
Request a detailed protocolSARS-CoV-2 Spike-protein virus was generated by replacing the classical Vesicular stomatitis virus (VSV-G) envelope protein of 3rd generation lentivirus with either wild-type Spike protein (‘Spike-WT’, [v1]) or a mutant Spike plasmid (‘Spike-mutant’) containing an ‘SRAS’ sequence in place of furin sensitive ‘RRAR’ (gift from Michael Farzan, Quinlan et al., 2020). In an additional variant called ‘Spike-delta’, the ‘SRAS’ site in Spike-mutant was replaced by a single ‘A’. To produce these viruses, Lenti-X 293 T cells (Takara Bio, MountainView, CA) were transiently transfected with 17.8 µg of pLKO.1 TRC-DsRed, 22 µg of psPAX2 and 5.8 µg of VSVG (or 9.9 µg of v1 or 8.6 µg of Spike-mutant/Spike-delta) plasmid in a 15 cm dish using the calcium phosphate method (Buffone et al., 2013). Six hours post-transfection, the medium was changed to virus collection medium (OPTIMEM + 10% FBS). The first batch of virus was collected 18–20 hr thereafter. Virus collection medium was then supplemented with 10 mM sodium butyrate. The second virus batch was collected 18–20 hr later. Both virus batches were pooled, filtered through 0.45 μm filter and ultra-centrifuged at 50,000 × g to concentrate the particles. The viral pellet was resuspended in OPTIMEM + 10% FBS, aliquoted and stored at −80°C until use.
Following the above protocol, 14 additional virus types were generated in this manuscript. This includes virus with envelope composed of: i. VSVG, ii. Spike-WT or iii. Spike-mutant that were produced in either: a. wild-type HEK293T, b. [O]-293 T c. [N]-293T or d. wild-type HEK293T being cultured in the presence of 15 µM kifunensine (Toronto Research Chemicals, Canada). In addition, Spike-delta was produced in the presence and absence of 15 µM kifunensine. For many of the production lots, viral titer was determined using the lentiviral p24 ELISA kit from Takara Bio (MountainView, CA) following manufacturer’s instructions. Virus concentration was also verified by western blotting with anti-S2 pAb or anti-FLAG to detect equivalent amounts of Spike protein in each of the preparations.
Pseudovirus entry assay
Request a detailed protocolIn order to transduce cells, virus was added at various titers indicated in the main manuscript to either wild-type 293Ts, stable 293T/ACE2 cell lines or various glycoengineered variants that transiently overexpressed ACE2 (at 48 hr post-transfection). Such infection was performed in the presence of 8 µg/ml polybrene as described previously (Buffone et al., 2013). In some cases, either the virus or cells were treated with sialidase for 1 hr at 37°C under conditions described above, prior to transduction. Virus was removed and fresh medium was added 8 hr post-transduction. DsRed fluorescence was monitored on a daily basis. All flow cytometry and microscopy data presented in this manuscript correspond to the 72 hr time point. Here, DsRed positive cells and the arithmetic mean fluorescence intensity (‘MFI’) of all cells in the PE-channel were also quantified using flow cytometry (BD Fortessa X-20). Microscopy images were acquired using a Zeiss AxioObserver instrument (10X/0.25 NA or 20X/0.4 NA objective).
Statistics
All data are presented as mean ± standard deviation. Dual comparisons were performed using the two-tailed Student’s t-test. ANOVA followed by the Tukey post-test was used for multiple comparisons. p<0.05 was considered to be statistically significant. Number of repeats are specified in individual panels.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files. All plasmids generated by the authors will be deposited at Addgene.
References
-
The proximal origin of SARS-CoV-2Nature Medicine 26:450–452.https://doi.org/10.1038/s41591-020-0820-9
-
Silencing α1,3-fucosyltransferases in human leukocytes reveals a role for FUT9 enzyme during E-selectin-mediated cell adhesionJournal of Biological Chemistry 288:1620–1633.https://doi.org/10.1074/jbc.M112.400929
-
Beyond shielding: The roles of glycans in the SARS-CoV-2 Spike proteinACS Central Science 6:1722–1734.https://doi.org/10.1021/acscentsci.0c01056
-
Clinical and immunological features of severe and moderate coronavirus disease 2019Journal of Clinical Investigation 130:2620–2629.https://doi.org/10.1172/JCI137244
-
Particle mesh Ewald - an n.log(N) Method for Ewald sums in large systemsJournal of Chemical Physics 98:10089–10092.https://doi.org/10.1063/1.464397
-
A potent inhibitor of the glycoprotein processing mannosidase IThe Journal of Biological Chemistry 265:15599–15605.
-
Phase IB clinical trial of the oligosaccharide processing inhibitor swainsonine in patients with advanced malignanciesClinical Cancer Res 3:1077–1086.
-
CHARMM-GUI: a web-based graphical user interface for CHARMMJournal of Computational Chemistry 29:1859–1865.https://doi.org/10.1002/jcc.20945
-
Glycan reader: automated sugar identification and simulation preparation for carbohydrates and glycoproteinsJournal of Computational Chemistry 32:3135–3141.https://doi.org/10.1002/jcc.21886
-
Disseminated intravascular coagulation in patients with 2019-nCoV pneumoniaJournal of Thrombosis and Haemostasis 18:786–787.https://doi.org/10.1111/jth.14781
-
Algorithms C-PM-D Constant pressure molecular dynamics algorithmsThe Journal of Chemical Physics 101:4177–4189.https://doi.org/10.1063/1.467468
-
Scalable molecular dynamics with NAMDJournal of Computational Chemistry 26:1781–1802.https://doi.org/10.1002/jcc.20289
-
Multiorgan and renal tropism of SARS-CoV-2New England Journal of Medicine 383:590–592.https://doi.org/10.1056/NEJMc2011400
-
The sweet spot: defining virus-sialic acid interactionsNature Reviews Microbiology 12:739–749.https://doi.org/10.1038/nrmicro3346
-
CHARMM general force field: a force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fieldsJournal of Computational Chemistry 31:671–690.https://doi.org/10.1002/jcc.21367
Article and author information
Author details
Funding
National Institutes of Health (HL103411)
- Sriram Neelamegham
National Institutes of Health (GM133195)
- Sriram Neelamegham
National Institutes of Health (GM126537)
- Sriram Neelamegham
National Institutes of Health (GM139160)
- Sheldon Park
- Sriram Neelamegham
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We gratefully acknowledge Dr. Michael Farzan (Scripps Research, Jupiter, FL) for providing HEK293T/ACE2 stable cells and Spike-mutant plasmid, and Dr. Haisheng Yu (Peking Union Medical College, China) for providing Spike-protein plasmid. This work was partially supported US National Institutes of Health grants HL103411, GM133195, GM139160 and GM126537.
Copyright
© 2020, Yang et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 9,110
- views
-
- 1,359
- downloads
-
- 181
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
The Parkinson’s disease (PD)-linked protein Leucine-Rich Repeat Kinase 2 (LRRK2) consists of seven domains, including a kinase and a Roc G domain. Despite the availability of several high-resolution structures, the dynamic regulation of its unique intramolecular domain stack is nevertheless still not well understood. By in-depth biochemical analysis, assessing the Michaelis–Menten kinetics of the Roc G domain, we have confirmed that LRRK2 has, similar to other Roco protein family members, a KM value of LRRK2 that lies within the range of the physiological GTP concentrations within the cell. Furthermore, the R1441G PD variant located within a mutational hotspot in the Roc domain showed an increased catalytic efficiency. In contrast, the most common PD variant G2019S, located in the kinase domain, showed an increased KM and reduced catalytic efficiency, suggesting a negative feedback mechanism from the kinase domain to the G domain. Autophosphorylation of the G1+2 residue (T1343) in the Roc P-loop motif is critical for this phosphoregulation of both the KM and the kcat values of the Roc-catalyzed GTP hydrolysis, most likely by changing the monomer–dimer equilibrium. The LRRK2 T1343A variant has a similar increased kinase activity in cells compared to G2019S and the double mutant T1343A/G2019S has no further increased activity, suggesting that T1343 is crucial for the negative feedback in the LRRK2 signaling cascade. Together, our data reveal a novel intramolecular feedback regulation of the LRRK2 Roc G domain by a LRRK2 kinase-dependent mechanism. Interestingly, PD mutants differently change the kinetics of the GTPase cycle, which might in part explain the difference in penetrance of these mutations in PD patients.
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
Pre-mRNA splicing is catalyzed in two steps: 5ʹ splice site (SS) cleavage and exon ligation. A number of proteins transiently associate with spliceosomes to specifically impact these steps (first and second step factors). We recently identified Fyv6 (FAM192A in humans) as a second step factor in Saccharomyces cerevisiae; however, we did not determine how widespread Fyv6’s impact is on the transcriptome. To answer this question, we have used RNA sequencing (RNA-seq) to analyze changes in splicing. These results show that loss of Fyv6 results in activation of non-consensus, branch point (BP) proximal 3ʹ SS transcriptome-wide. To identify the molecular basis of these observations, we determined a high-resolution cryo-electron microscopy (cryo-EM) structure of a yeast product complex spliceosome containing Fyv6 at 2.3 Å. The structure reveals that Fyv6 is the only second step factor that contacts the Prp22 ATPase and that Fyv6 binding is mutually exclusive with that of the first step factor Yju2. We then use this structure to dissect Fyv6 functional domains and interpret results of a genetic screen for fyv6Δ suppressor mutations. The combined transcriptomic, structural, and genetic studies allow us to propose a model in which Yju2/Fyv6 exchange facilitates exon ligation and Fyv6 promotes usage of consensus, BP distal 3ʹ SS.