Role of Q675H Mutation in Improving SARS-CoV-2 Spike Interaction with the Furin Binding Pocket

The impact of S2 mutations on Omicron SARS-CoV-2 cell surface expression and fusogenicity

SARS-CoV-2 Omicron subvariants are still emerging and spreading worldwide. These variants contain a high number of polymorphisms in the spike (S) glycoprotein that could potentially impact their pathogenicity and transmission. We have previously shown that the S:655Y and P681H mutations enhance S protein cleavage and syncytia formation. Interestingly, these polymorphisms are present in Omicron S protein. Here, we characterized the cleavage efficiency and fusogenicity of the S protein of different Omicron sublineages. Our results showed that Omicron BA.1 subvariant is efficiently cleaved but it is poorly fusogenic compared to previous SARS-CoV-2 strains. To understand the basis of this phenotype, we generated chimeric S protein using combinations of the S1 and S2 domains from WA1, Delta and Omicron BA.1 variants. We found that the S2 domain of Omicron BA.1 hindered efficient cell-cell fusion. Interestingly, this domain only contains six unique polymorphisms never detected before in ancestral SARS-CoV-2 variants. WA1614G S proteins containing the six individuals S2 Omicron mutations were assessed for their fusogenicity and S surface expression after transfection in cells. Results showed that the S:N856K and N969K substitutions decreased syncytia formation and impacted S protein cell surface levels. However, we observed that "first-generation" Omicron sublineages that emerged subsequently, had convergently evolved to an enhanced fusogenic activity and S expression on the surface of infected cells while "second-generation" Omicron variants have highly diverged and showed lineage-specific fusogenic properties. Importantly, our findings could have potential implications in the improvement and redesign of COVID-19 vaccines.

Evolution of SARS-CoV-2 in the murine central nervous system drives viral diversification

Severe coronavirus disease 2019 and post-acute sequelae of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection are associated with neurological complications that may be linked to direct infection of the central nervous system (CNS), but the selective pressures ruling neuroinvasion are poorly defined. Here we assessed SARS-CoV-2 evolution in the lung versus CNS of infected mice. Higher levels of viral divergence were observed in the CNS than the lung after intranasal challenge with a high frequency of mutations in the spike furin cleavage site (FCS). Deletion of the FCS significantly attenuated virulence after intranasal challenge, with lower viral titres and decreased morbidity compared with the wild-type virus. Intracranial inoculation of the FCS-deleted virus, however, was sufficient to restore virulence. After intracranial inoculation, both viruses established infection in the lung, but dissemination from the CNS to the lung required the intact FCS. Cumulatively, these data suggest a critical role for the FCS in determining SARS-CoV-2 tropism and compartmentalization.

Melodia: A Python Library for Protein Structure Analysis

SUMMARY

Analysing protein structure similarities is an important step in protein engineering and drug discovery. Methodologies that are more advanced than simple RMSD are available but often require extensive mathematical or computational knowledge for implementation. Grouping and optimising such tools in an efficient open-source library increases accessibility and encourages the adoption of more advanced metrics. Melodia is a Python library with a complete set of components devised for describing, comparing and analysing the shape of protein structures using differential geometry of three-dimensional curves and knot theory. It can generate robust geometric descriptors for thousands of shapes in just a few minutes. Those descriptors are more sensitive to structural feature variation than RMSD deviation. Melodia also incorporates sequence structural annotation and three-dimensional visualisations.

AVAILABILITY AND IMPLEMENTATION

Melodia is an open-source Python library freely available on https://github.com/rwmontalvao/Melodia_py, along with interactive Jupyter Notebook tutorials.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Early Emerging SARS-CoV-2 Spike Mutants Are Diversified in Virologic Properties but Elicit Compromised Antibody Responses

Despite the effective antivirals and vaccines, COVID-19 remains a public health concern. The mutations that occurred during the early stage of the pandemic can be valuable in assessing the viral fitness and evolutionary trajectory. In this study, we analyzed a panel of 2969 spike sequences deposited in GISAID before April 2020 and characterized nine representative spike single-point mutants in detail. Compared with the WA01/2020, most (8 out of 9) mutants demonstrated an equivalent or diminished protein expression or processing, pseudovirus infectivity, and cell-cell fusion. Interestingly, most of the mutants in native form elicited minimum antibody responses in mice despite unaltered CD4+ and CD8+ T cell responses. The mutants remained sensitive to the antisera and the type I interferon. Taken together, these data suggest that the early emerging mutants are virologically divergent, and some of which showed transmission fitness. Our findings have important implications for the retrospective tracing of the early SARS-CoV-2 transmission and future pandemic preparedness.

The high mutation rate at the D614G hotspot-furin cleavage site region increases the priming efficiency of the Spike protein by furin protease: analysis of Indonesian SARS-CoV-2 G614 variants obtained during the early COVID-19 pandemic

D614G mutation plays a significant role in the transmissibility of SARS-CoV-2. Identification of other mutations related to D614G mutation within the Spike protein is pivotal as they might contribute to the pathogenicity of SARS-CoV-2. This study aims to analyze the mutation rate of furin cleavage site (FCS) region of Indonesian origin SARS-CoV-2 and to predict the effect of mutation against Spike priming efficiency by furin. A total of 375 sequences of Indonesian isolates obtained during the early pandemic were used for mutation analysis. Mutation analysis includes mutation pattern, variability, frequency of mutation, amino acid conservation, and mutation rate. The effect of mutation against Spike priming efficiency by furin protease from eight sequences with mutation in the FCS region was analyzed by protein-protein docking. We showed that mutations related to the G614 variant were increasing through time, in contrast to the D614 variant. The FCS region at the position 675-692 contained the most variable (66.67%) as well as the highest mutation frequency (85.92%) and has been observed to be the hotspot mutations linked to the D614G mutation. The D614G hotspot-FCS region (residue 600-700) had the highest amino acid change per site (20.8%) as well as the highest mutation rate as 1.34 × 10^-2 substitution per site per year (95% CI 1.79 × 10^-3-2.74 × 10^-2), compared with other Spike protein regions. Mutations in the FCS region were the most common mutation found after the D614G mutation. These mutations were predicted to increase the Spike priming efficiency by furin. Thus, this study elucidates the importance of D614G mutation to other mutations located in the FCS region and their significance to Spike priming efficiency by furin.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13337-023-00827-w.

O-Linked Sialoglycans Modulate the Proteolysis of SARS-CoV-2 Spike and Likely Contribute to the Mutational Trajectory in Variants of Concern

The emergence of a polybasic cleavage motif for the protease furin in SARS-CoV-2 spike has been established as a major factor for human viral transmission. The region N-terminal to that motif is extensively mutated in variants of concern (VOCs). Besides furin, spikes from these variants appear to rely on other proteases for maturation, including TMPRSS2. Glycans near the cleavage site have raised questions about proteolytic processing and the consequences of variant-borne mutations. Here, we identify that sialic acid-containing O-linked glycans on Thr678 of SARS-CoV-2 spike influence furin and TMPRSS2 cleavage and posit O-linked glycosylation as a likely driving force for the emergence of VOC mutations. We provide direct evidence that the glycosyltransferase GalNAc-T1 primes glycosylation at Thr678 in the living cell, an event that is suppressed by mutations in the VOCs Alpha, Delta, and Omicron. We found that the sole incorporation of N-acetylgalactosamine did not impact furin activity in synthetic O-glycopeptides, but the presence of sialic acid reduced the furin rate by up to 65%. Similarly, O-glycosylation with a sialylated trisaccharide had a negative impact on TMPRSS2 cleavage. With a chemistry-centered approach, we substantiate O-glycosylation as a major determinant of spike maturation and propose disruption of O-glycosylation as a substantial driving force for VOC evolution.

SARS-CoV-2 S Mutations: A Lesson from the Viral World to Understand How Human Furin Works

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is the etiological agent responsible for the worldwide pandemic and has now claimed millions of lives. The virus combines several unusual characteristics and an extraordinary ability to spread among humans. In particular, the dependence of the maturation of the envelope glycoprotein S from Furin enables the invasion and replication of the virus virtually within the entire body, since this cellular protease is ubiquitously expressed. Here, we analyzed the naturally occurring variation of the amino acids sequence around the cleavage site of S. We found that the virus grossly mutates preferentially at P positions, resulting in single residue replacements that associate with gain-of-function phenotypes in specific conditions. Interestingly, some combinations of amino acids are absent, despite the evidence supporting some cleavability of the respective synthetic surrogates. In any case, the polybasic signature is maintained and, as a consequence, Furin dependence is preserved. Thus, no escape variants to Furin are observed in the population. Overall, the SARS-CoV-2 system per se represents an outstanding example of the evolution of substrate-enzyme interaction, demonstrating a fast-tracked optimization of a protein stretch towards the Furin catalytic pocket. Ultimately, these data disclose important information for the development of drugs targeting Furin and Furin-dependent pathogens.

Coordinated evolution at amino acid sites of SARS-CoV-2 spike

SARS-CoV-2 has adapted in a stepwise manner, with multiple beneficial mutations accumulating in a rapid succession at origins of VOCs, and the reasons for this are unclear. Here, we searched for coordinated evolution of amino acid sites in the spike protein of SARS-CoV-2. Specifically, we searched for concordantly evolving site pairs (CSPs) for which changes at one site were rapidly followed by changes at the other site in the same lineage. We detected 46 sites which formed 45 CSP. Sites in CSP were closer to each other in the protein structure than random pairs, indicating that concordant evolution has a functional basis. Notably, site pairs carrying lineage defining mutations of the four VOCs that circulated before May 2021 are enriched in CSPs. For the Alpha VOC, the enrichment is detected even if Alpha sequences are removed from analysis, indicating that VOC origin could have been facilitated by positive epistasis. Additionally, we detected nine discordantly evolving pairs of sites where mutations at one site unexpectedly rarely occurred on the background of a specific allele at another site, for example on the background of wild-type D at site 614 (four pairs) or derived Y at site 501 (three pairs). Our findings hint that positive epistasis between accumulating mutations could have delayed the assembly of advantageous combinations of mutations comprising at least some of the VOCs.

Evolution of SARS-CoV-2 during the first year of the COVID-19 pandemic in Northwestern Argentina

Studies about the evolution of SARS-CoV-2 lineages in different backgrounds such as naive populations are still scarce, especially from South America. This work aimed to study the introduction and diversification pattern of SARS-CoV-2 during the first year of the COVID-19 pandemic in the Northwestern Argentina (NWA) region and to analyze the evolutionary dynamics of the main lineages found. In this study, we analyzed a total of 260 SARS-CoV-2 whole-genome sequences from Argentina, belonging to the Provinces of Jujuy, Salta, and Tucumán, from March 31st, 2020, to May 22nd, 2021, which covered the full first wave and the early second wave of the COVID-19 pandemic in Argentina. In the first wave, eight lineages were identified: B.1.499 (76.9%), followed by N.5 (10.2%), B.1.1.274 (3.7%), B.1.1.348 (3.7%), B.1 (2.8%), B.1.600 (0.9%), B.1.1.33 (0.9%) and N.3 (0.9%). During the early second wave, the first-wave lineages were displaced by the introduction of variants of concern (VOC) (Alpha, Gamma), or variants of interest (VOI) (Lambda, Zeta, Epsilon) and other lineages with more limited distribution. Phylodynamic analyses of the B.1.499 and N.5, the two most prevalent lineages in the NWA, revealed that the rate of evolution of lineage N.5 (7.9 × 10-4 substitutions per site per year, s/s/y) was a ∼40% faster than that of lineage B.1.499 (5.6 × 10-4 s/s/y), although both are in the same order of magnitude than other non-VOC lineages. No mutations associated with a biological characteristic of importance were observed as signatures markers of the phylogenetic groups established in Northwestern Argentina, however, single sequences in non-VOC lineages did present mutations of biological importance or associated with VOCs as sporadic events, showing that many of these mutations could emerge from circulation in the general population. This study contributed to the knowledge about the evolution of SARS-CoV-2 in a pre-vaccination and without post-exposure immunization period.

Spike Mutation Profiles Associated With SARS-CoV-2 Breakthrough Infections in Delta Emerging and Predominant Time Periods in British Columbia, Canada

Background

COVID-19 vaccination is a key public health measure in the pandemic response. The rapid evolution of SARS-CoV-2 variants introduce new groups of spike protein mutations. These new mutations are thought to aid in the evasion of vaccine-induced immunity and render vaccines less effective. However, not all spike mutations contribute equally to vaccine escape. Previous studies associate mutations with vaccine breakthrough infections (BTI), but information at the population level remains scarce. We aimed to identify spike mutations associated with SARS-CoV-2 vaccine BTI in a community setting during the emergence and predominance of the Delta-variant.

Methods

This case-control study used both genomic, and epidemiological data from a provincial COVID-19 surveillance program. Analyses were stratified into two periods approximating the emergence and predominance of the Delta-variant, and restricted to primary SARS-CoV-2 infections from either unvaccinated individuals, or those infected ≥14 days after their second vaccination dose in a community setting. Each sample's spike mutations were concatenated into a unique spike mutation profile (SMP). Penalized logistic regression was used to identify spike mutations and SMPs associated with SARS-CoV-2 vaccine BTI in both time periods.

Results and Discussion

This study reports population level relative risk estimates, between 2 and 4-folds, of spike mutation profiles associated with BTI during the emergence and predominance of the Delta-variant, which comprised 19,624 and 17,331 observations, respectively. The identified mutations cover multiple spike domains including the N-terminal domain (NTD), receptor binding domain (RBD), S1/S2 cleavage region, fusion peptide and heptad regions. Mutations in these different regions imply various mechanisms contribute to vaccine escape. Our profiling method identifies naturally occurring spike mutations associated with BTI, and can be applied to emerging SARS-CoV-2 variants with novel groups of spike mutations.

SARS-CoV-2 Infects Human ACE2-Negative Endothelial Cells through an αvβ3 Integrin-Mediated Endocytosis Even in the Presence of Vaccine-Elicited Neutralizing Antibodies

Integrins represent a gateway of entry for many viruses and the Arg-Gly-Asp (RGD) motif is the smallest sequence necessary for proteins to bind integrins. All Severe Acute Respiratory Syndrome Virus type 2 (SARS-CoV-2) lineages own an RGD motif (aa 403–405) in their receptor binding domain (RBD). We recently showed that SARS-CoV-2 gains access into primary human lung microvascular endothelial cells (HL-mECs) lacking Angiotensin-converting enzyme 2 (ACE2) expression through this conserved RGD motif. Following its entry, SARS-CoV-2 remodels cell phenotype and promotes angiogenesis in the absence of productive viral replication. Here, we highlight the αvβ3 integrin as the main molecule responsible for SARS-CoV-2 infection of HL-mECs via a clathrin-dependent endocytosis. Indeed, pretreatment of virus with αvβ3 integrin or pretreatment of cells with a monoclonal antibody against αvβ3 integrin was found to inhibit SARS-CoV-2 entry into HL-mECs. Surprisingly, the anti-Spike antibodies evoked by vaccination were neither able to impair Spike/integrin interaction nor to prevent SARS-CoV-2 entry into HL-mECs. Our data highlight the RGD motif in the Spike protein as a functional constraint aimed to maintain the interaction of the viral envelope with integrins. At the same time, our evidences call for the need of intervention strategies aimed to neutralize the SARS-CoV-2 integrin-mediated infection of ACE2-negative cells in the vaccine era.

Capturing a Crucial ‘Disorder-to-Order Transition’ at the Heart of the Coronavirus Molecular Pathology—Triggered by Highly Persistent, Interchangeable Salt-Bridges

The COVID-19 origin debate has greatly been influenced by genome comparison studies of late, revealing the emergence of the Furin-like cleavage site at the S1/S2 junction of the SARS-CoV-2 Spike (FLCSSpike) containing its 681PRRAR685 motif, absent in other related respiratory viruses. Being the rate-limiting (i.e., the slowest) step, the host Furin cleavage is instrumental in the abrupt increase in transmissibility in COVID-19, compared to earlier onsets of respiratory viral diseases. In such a context, the current paper entraps a ‘disorder-to-order transition’ of the FLCSSpike (concomitant to an entropy arrest) upon binding to Furin. The interaction clearly seems to be optimized for a more efficient proteolytic cleavage in SARS-CoV-2. The study further shows the formation of dynamically interchangeable and persistent networks of salt-bridges at the Spike–Furin interface in SARS-CoV-2 involving the three arginines (R682, R683, R685) of the FLCSSpike with several anionic residues (E230, E236, D259, D264, D306) coming from Furin, strategically distributed around its catalytic triad. Multiplicity and structural degeneracy of plausible salt-bridge network archetypes seem to be the other key characteristic features of the Spike–Furin binding in SARS-CoV-2, allowing the system to breathe—a trademark of protein disorder transitions. Interestingly, with respect to the homologous interaction in SARS-CoV (2002/2003) taken as a baseline, the Spike–Furin binding events, generally, in the coronavirus lineage, seems to have preference for ionic bond formation, even with a lesser number of cationic residues at their potentially polybasic FLCSSpike patches. The interaction energies are suggestive of characteristic metastabilities attributed to Spike–Furin interactions, generally to the coronavirus lineage, which appears to be favorable for proteolytic cleavages targeted at flexible protein loops. The current findings not only offer novel mechanistic insights into the coronavirus molecular pathology and evolution, but also add substantially to the existing theories of proteolytic cleavages.