skip to main content

Title: SARS-CoV-2 spike opening dynamics and energetics reveal the individual roles of glycans and their collective impact
Abstract The trimeric spike (S) glycoprotein, which protrudes from the SARS-CoV-2 viral envelope, binds to human ACE2, initiated by at least one protomer’s receptor binding domain (RBD) switching from a "down” (closed) to an "up” (open) state. Here, we used large-scale molecular dynamics simulations and two-dimensional replica exchange umbrella sampling calculations with more than a thousand windows and an aggregate total of 160 μ s of simulation to investigate this transition with and without glycans. We find that the glycosylated spike has a higher barrier to opening and also energetically favors the down state over the up state. Analysis of the S-protein opening pathway reveals that glycans at N165 and N122 interfere with hydrogen bonds between the RBD and the N-terminal domain in the up state, while glycans at N165 and N343 can stabilize both the down and up states. Finally, we estimate how epitope exposure for several known antibodies changes along the opening path. We find that the BD-368-2 antibody’s epitope is continuously exposed, explaining its high efficacy.  more » « less
Award ID(s):
1828187 1806833
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Communications Biology
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The glycosylation on the spike (S) protein of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, modulates the viral infection by altering conformational dynamics, receptor interaction and host immune responses. Several variants of concern (VOCs) of SARS-CoV-2 have evolved during the pandemic, and crucial mutations on the S protein of the virus have led to increased transmissibility and immune escape. In this study, we compare the site-specific glycosylation and overall glycomic profiles of the wild type Wuhan-Hu-1 strain (WT) S protein and five VOCs of SARS-CoV-2: Alpha, Beta, Gamma, Delta and Omicron. Interestingly, both N- and O-glycosylation sites on the S protein are highly conserved among the spike mutant variants, particularly at the sites on the receptor-binding domain (RBD). The conservation of glycosylation sites is noteworthy, as over 2 million SARS-CoV-2 S protein sequences have been reported with various amino acid mutations. Our detailed profiling of the glycosylation at each of the individual sites of the S protein across the variants revealed intriguing possible association of glycosylation pattern on the variants and their previously reported infectivity. While the sites are conserved, we observed changes in the N- and O-glycosylation profile across the variants. The newly emerged variants, which showed higher resistance to neutralizing antibodies and vaccines, displayed a decrease in the overall abundance of complex-type glycans with both fucosylation and sialylation and an increase in the oligomannose-type glycans across the sites. Among the variants, the glycosylation sites with significant changes in glycan profile were observed at both theN-terminal domain and RBD of S protein, with Omicron showing the highest deviation. The increase in oligomannose-type happens sequentially from Alpha through Delta. Interestingly, Omicron does not contain more oligomannose-type glycans compared to Delta but does contain more compared to the WT and other VOCs. O-glycosylation at the RBD showed lower occupancy in the VOCs in comparison to the WT. Our study on the sites and pattern of glycosylation on the SARS-CoV-2 S proteins across the VOCs may help to understand how the virus evolved to trick the host immune system. Our study also highlights how the SARS-CoV-2 virus has conserved bothN- andO- glycosylation sites on the S protein of the most successful variants even after undergoing extensive mutations, suggesting a correlation between infectivity/ transmissibility and glycosylation.

    more » « less
  2. Infection by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) involves the attachment of the receptor-binding domain (RBD) of its spike proteins to the ACE2 receptors on the peripheral membrane of host cells. Binding is initiated by a down-to-up conformational change in the spike protein, the change that presents the RBD to the receptor. To date, computational and experimental studies that search for therapeutics have concentrated, for good reason, on the RBD. However, the RBD region is highly prone to mutations, and is therefore a hotspot for drug resistance. In contrast, we here focus on the correlations between the RBD and residues distant to it in the spike protein. This allows for a deeper understanding of the underlying molecular recognition events and prediction of the highest-effect key mutations in distant, allosteric sites, with implications for therapeutics. Also, these sites can appear in emerging mutants with possibly higher transmissibility and virulence, and preidentifying them can give clues for designing pan-coronavirus vaccines against future outbreaks. Our model, based on time-lagged independent component analysis (tICA) and protein graph connectivity network, is able to identify multiple residues that exhibit long-distance coupling with the RBD opening. Residues involved in the most ubiquitous D614G mutation and the A570D mutation of the highly contagious UK SARS-CoV-2 variant are predicted ab initio from our model. Conversely, broad-spectrum therapeutics like drugs and monoclonal antibodies can target these key distant-but-conserved regions of the spike protein. 
    more » « less
  3. With the increased prevalence of new SARS-CoV-2 variants of concern, such as Delta and Omicron, the COVID-19 pandemic has become an ongoing human health disaster, killing millions worldwide. SARS-CoV-2 invades its host through the interaction of its spike (S) protein with a host cell receptor, angiotensin-converting enzyme 2 (ACE2). In addition, heparan sulfate (HS) on the surface of host cells plays an important role as a co-receptor for this viral pathogen–host cell interaction. Our previous studies demonstrated that many sulfated glycans, such as heparin, fucoidans, and rhamnan sulfate have anti-SARS-CoV-2 activities. In the current study, a small library of sulfated glycans and highly negatively charged compounds, including pentosan polysulfate (PPS), mucopolysaccharide polysulfate (MPS), sulfated lactobionic acid, sulodexide, and defibrotide, was assembled and evaluated for binding to the S-proteins and inhibition of viral infectivity in vitro. These compounds inhibited the interaction of the S-protein receptor-binding domain (RBD) (wild type and different variants) with immobilized heparin, a highly sulfated HS, as determined using surface plasmon resonance (SPR). PPS and MPS showed the strongest inhibition of interaction of heparin and S-protein RBD. The competitive binding studies showed that the IC50 of PPS and MPS against the S-protein RBD binding to immobilized heparin was ~35 nM and ~9 nM, respectively, much lower than the IC50 for soluble heparin (IC50 = 56 nM). Both PPS and MPS showed stronger inhibition than heparin on the S-protein RBD or spike pseudotyped lentiviral particles binding to immobilized heparin. Finally, in an in vitro cell-based assay, PPS and MPS exhibited strong antiviral activities against pseudotyped viral particles of SARS-CoV-2 containing wild-type or Delta S-proteins. 
    more » « less
  4. The COVID-19 pandemic is a major human health concern. The pathogen responsible for COVID-19, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), invades its host through the interaction of its spike (S) protein with a host cell receptor, angiotensin-converting enzyme 2 (ACE2). In addition to ACE2, heparan sulfate (HS) on the surface of host cells also plays a significant role as a co-receptor. Our previous studies demonstrated that sulfated glycans, such as heparin and fucoidans, show anti-COVID-19 activities. In the current study, rhamnan sulfate (RS), a polysaccharide with a rhamnose backbone from a green seaweed, Monostroma nitidum, was evaluated for binding to the S-protein from SARS-CoV-2 and inhibition of viral infectivity in vitro. The structural characteristics of RS were investigated by determining its monosaccharide composition and performing two-dimensional nuclear magnetic resonance. RS inhibition of the interaction of heparin, a highly sulfated HS, with the SARS-CoV-2 spike protein (from wild type and different mutant variants) was studied using surface plasmon resonance (SPR). In competitive binding studies, the IC50 of RS against the S-protein receptor binding domain (RBD) binding to immobilized heparin was 1.6 ng/mL, which is much lower than the IC50 for heparin (~750 ng/mL). RS showed stronger inhibition than heparin on the S-protein RBD or pseudoviral particles binding to immobilized heparin. Finally, in an in vitro cell-based assay, RS showed strong antiviral activities against wild type SARS-CoV-2 and the delta variant. 
    more » « less
  5. Abstract

    The continued emergence of new SARS‐CoV‐2 variants has accentuated the growing need for fast and reliable methods for the design of potentially neutralizing antibodies (Abs) to counter immune evasion by the virus. Here, we report on the de novo computational design of high‐affinity Ab variable regions (Fv) through the recombination of VDJ genes targeting the most solvent‐exposed hACE2‐binding residues of the SARS‐CoV‐2 spike receptor binding domain (RBD) protein using the software toolOptMAVEn‐2.0. Subsequently, we carried out computational affinity maturation of the designed variable regions through amino acid substitutions for improved binding with the target epitope. Immunogenicity of designs was restricted by preferring designs that match sequences from a 9‐mer library of “human Abs” based on a human string content score. We generated 106 different antibody designs and reported in detail on the top five that trade‐off the greatest computational binding affinity for the RBD with human string content scores. We further describe computational evaluation of the top five designs produced byOptMAVEn‐2.0using a Rosetta‐based approach. We used RosettaSnugDockfor local docking of the designs to evaluate their potential to bind the spike RBD and performed “forward folding” withDeepAbto assess their potential to fold into the designed structures. Ultimately, our results identified one designed Ab variable region, P1.D1, as a particularly promising candidate for experimental testing. This effort puts forth a computational workflow for the de novo design and evaluation of Abs that can quickly be adapted to target spike epitopes of emerging SARS‐CoV‐2 variants or other antigenic targets.

    more » « less