SARS-CoV-2 variants of concern: spike protein mutational analysis and epitope for broad neutralization

Spike Protein-Fibrinogen Interaction: A Novel Immune Evasion Strategy of SARS-CoV-2?

The host protein fibrinogen has been found to interact with the N-terminal domain (NTD) of the spike protein in SARS-CoV-2. However, the evolutionary benefit of this binding to the virus still remains unclear. Herein, we put forward with rationale and supporting evidence that the binding of fibrinogen to its more conserved NTD is an immune evasion strategy adopted by the virus to outsmart the NTD targeted neutralizing antibodies.

Post-viral lung diseases: the microbiota as a key player

Viral infections of the respiratory tract can lead to chronic lung injury through immunopathological mechanisms that remain unclear. Communities of commensal bacteria colonising the respiratory tract, known as the respiratory tract microbiota, are altered in viral infections, which can contribute to inflammation, lung epithelial damage and subsequent development of lung disease. Emerging evidence on post-viral lung injury suggests an interplay between viral infections, immune responses and airway microbiota composition in the development of viral-induced lung diseases. In this review, we present the clinical characteristics of post-viral lung injury, along with the underlying immunopathological mechanisms and host-bacteria interactions, with a focus on influenza virus, respiratory syncytial virus and coronaviruses. Additionally, considering the important role of the airway microbiota in viral-induced pulmonary sequelae, we suggest key areas for future research on respiratory microbiota involvement in the development of post-viral lung diseases.

Interferon therapy in alpha and Delta variants of SARS-CoV-2: The dichotomy between laboratory success and clinical realities

The COVID-19 pandemic has caused significant morbidity and mortality worldwide. The emergence of the Alpha and Delta variants of SARS-CoV-2 has led to a renewed interest in using interferon therapy as a potential treatment option. Interferons are a group of signaling proteins produced by host cells in response to viral infections. They play a critical role in the innate immune response to viral infections by inducing an antiviral state in infected and neighboring cells. Interferon therapy has shown promise as a potential treatment option for COVID-19. In this review paper, we review the current knowledge regarding interferon therapy in the context of the Alpha and Delta variants of SARS-CoV-2 and discuss the challenges that must be overcome to translate laboratory findings into effective clinical treatments.

Interfacial subregions of SARS-CoV-2 spike RBD to hACE2 affect intermolecular affinity by their distinct roles played in association and dissociation kinetics

SARS-CoV-2's rapid global transmission depends on spike RBD's strong affinity to hACE2. In the context of binding hot spots well defined, the work investigated how interfacial subregions of SARS-CoV-2 spike RBD to hACE2 affect intermolecular affinity and their potential distinct roles involved in association and dissociation kinetics due to their local structural characteristics. Three spatially consecutive subregions of SARS-CoV-2 RBD were structurally partitioned across RBD's receptor binding motif (RBM). Their impacts on binding affinity and kinetics were differentiated through a comprehensive SPR measurement of hACE2 binding by chimeric swap mutants of respective subdomains from SARS-CoV-2 VOCs & phylogenetically close sarbecoviruses, and further compared with those of included single mutations across RBM and around the RBD core. The data supports that the intermediate interfacial subregion of RBD involving key residue at 417 is the rate-limiting effector of association kinetics and the subregion encompassing residues at 501/498/449 is the key binding energy contributor dictating dissociation kinetics, both of which relate to SARS-CoV-2's adaptive mutational evolution and host tropism closely. The kinetic data and structural analysis of local mutations' impact on spike RBD's binding and thermal stability provide a new perspective in evaluating SARS-CoV-2 evolution and other sarbecoviruses' evolvable binding to hACE2. The inherent binding mode offers direct clues of valid epitope in designing new antibodies that the coronavirus can't elude.

Genomic Evolution of the SARS-CoV-2 Omicron Variant in Córdoba, Argentina (2021-2022): Analysis of Uncommon and Prevalent Spike Mutations

Understanding the evolutionary patterns and geographic spread of SARS-CoV-2 variants, particularly Omicron, is essential for effective public health responses. This study focused on the genomic analysis of the Omicron variant in Cordoba, Argentina from 2021 to 2022. Phylogenetic analysis revealed the dominant presence of BA.1 and BA.2 lineages, with BA.5 emerging earlier than BA.4, aligning with observations from other regions. Haplotype network analysis showed significant genetic divergence within Omicron samples, forming distinct clusters. In comparison to global datasets, we identified mutations in the Omicron genomes (A27S, Y145D, and L212I) situated within the NTD region of the Spike protein. These mutations, while not widespread globally, showed higher prevalence in our region. Of particular interest were the Y145D and L212I substitutions, previously unreported in Argentina. In silico analysis revealed that both mutations impact the binding affinity of T-cell epitopes to HLA type I and II alleles. Notably, these alleles are among the most common in the Argentinian population, with some associated with protection against and others with susceptibility to SARS-CoV-2 infection. These findings strongly suggest that these prevalent mutations likely influence the immunogenicity of the Spike protein and contribute to immune evasion mechanisms. This study provides valuable insights into the genomic dynamics of the Omicron variant in Cordoba, Argentina and highlights unique mutations with potential implications for COVID-19 vaccines.

Virion morphology and on-virus spike protein structures of diverse SARS-CoV-2 variants

The evolution of SARS-CoV-2 variants with increased fitness has been accompanied by structural changes in the spike (S) proteins, which are the major target for the adaptive immune response. Single-particle cryo-EM analysis of soluble S protein from SARS-CoV-2 variants has revealed this structural adaptation at high resolution. The analysis of S trimers in situ on intact virions has the potential to provide more functionally relevant insights into S structure and virion morphology. Here, we characterized B.1, Alpha, Beta, Gamma, Delta, Kappa, and Mu variants by cryo-electron microscopy and tomography, assessing S cleavage, virion morphology, S incorporation, "in-situ" high-resolution S structures, and the range of S conformational states. We found no evidence for adaptive changes in virion morphology, but describe multiple different positions in the S protein where amino acid changes alter local protein structure. Taken together, our data are consistent with a model where amino acid changes at multiple positions from the top to the base of the spike cause structural changes that can modulate the conformational dynamics of the S protein.

Comparison of the Neutralization Power of Sotrovimab Against SARS-CoV-2 Variants: Development of a Rapid Computational Method

BACKGROUND

The rapid evolution of SARS-CoV-2 imposed a huge challenge on disease control. Immune evasion caused by genetic variations of the SARS-CoV-2 spike protein's immunogenic epitopes affects the efficiency of monoclonal antibody-based therapy of COVID-19. Therefore, a rapid method is needed to evaluate the efficacy of the available monoclonal antibodies against the new emerging variants or potential novel variants.

OBJECTIVE

The aim of this study is to develop a rapid computational method to evaluate the neutralization power of anti-SARS-CoV-2 monoclonal antibodies against new SARS-CoV-2 variants and other potential new mutations.

METHODS

The amino acid sequence of the extracellular domain of the spike proteins of the severe acute respiratory syndrome coronavirus (GenBank accession number YP_009825051.1) and SARS-CoV-2 (GenBank accession number YP_009724390.1) were used to create computational 3D models for the native spike proteins. Specific mutations were introduced to the curated sequence to generate the different variant spike models. The neutralization potential of sotrovimab (S309) against these variants was evaluated based on its molecular interactions and Gibbs free energy in comparison to a reference model after molecular replacement of the reference receptor-binding domain with the variant's receptor-binding domain.

RESULTS

Our results show a loss in the binding affinity of the neutralizing antibody S309 with both SARS-CoV and SARS-CoV-2. The binding affinity of S309 was greater to the Alpha, Beta, Gamma, and Kappa variants than to the original Wuhan strain of SARS-CoV-2. However, S309 showed a substantially decreased binding affinity to the Delta and Omicron variants. Based on the mutational profile of Omicron subvariants, our data describe the effect of the G339H and G339D mutations and their role in escaping antibody neutralization, which is in line with published clinical reports.

CONCLUSIONS

This method is rapid, applicable, and of interest to adapt the use of therapeutic antibodies to the treatment of emerging variants. It could be applied to antibody-based treatment of other viral infections.

In Silico Design of miniACE2 Decoys with In Vitro Enhanced Neutralization Activity against SARS-CoV-2, Encompassing Omicron Subvariants

The COVID-19 pandemic has overwhelmed healthcare systems and triggered global economic downturns. While vaccines have reduced the lethality rate of SARS-CoV-2 to 0.9% as of October 2024, the continuous evolution of variants remains a significant public health challenge. Next-generation medical therapies offer hope in addressing this threat, especially for immunocompromised individuals who experience prolonged infections and severe illnesses, contributing to viral evolution. These cases increase the risk of new variants emerging. This study explores miniACE2 decoys as a novel strategy to counteract SARS-CoV-2 variants. Using in silico design and molecular dynamics, blocking proteins (BPs) were developed with stronger binding affinity for the receptor-binding domain of multiple variants than naturally soluble human ACE2. The BPs were expressed in E. coli and tested in vitro, showing promising neutralizing effects. Notably, miniACE2 BP9 exhibited an average IC50 of 4.9 µg/mL across several variants, including the Wuhan strain, Mu, Omicron BA.1, and BA.2 This low IC50 demonstrates the potent neutralizing ability of BP9, indicating its efficacy at low concentrations.Based on these findings, BP9 has emerged as a promising therapeutic candidate for combating SARS-CoV-2 and its evolving variants, thereby positioning it as a potential emergency biopharmaceutical.

Booster Vaccination with BNT162b2 Improves Cellular and Humoral Immune Response in the Pediatric Population Immunized with CoronaVac

The SARS-CoV-2 Omicron variant and its sublineages continue to cause COVID-19-associated pediatric hospitalizations, severe disease, and death globally. BNT162b2 and CoronaVac are the main vaccines used in Chile. Much less is known about the Wuhan-Hu-1 strain-based vaccines in the pediatric population compared to adults. Given the worldwide need for booster vaccinations to stimulate the immune response against new Omicron variants of SARS-CoV-2, we characterized the humoral and cellular immune response against Omicron variant BA.1 in a pediatric cohort aged 10 to 16 years who received heterologous vaccination based on two doses of CoronaVac, two doses of CoronaVac (2x) plus one booster dose of BNT162b2 [CoronaVac(2x) + BNT162b2 (1x)], two doses of CoronaVac plus two booster doses of BNT162b2 [CoronaVac(2x) + BNT162b2 (2x)], and three doses of BNT162b2. We observed that the [CoronaVac(2x) + BNT162b2 (2x)] vaccination showed higher anti-S1 and neutralizing antibody titers and CD4 and CD8 T cell immunity specific to the Omicron variant compared to immunization with two doses of CoronaVac alone. Furthermore, from all groups tested, immunity against Omicron was highest in individuals who received three doses of BNT162b2. We conclude that booster vaccination with BNT162b2, compared to two doses of CoronaVac alone, induces a greater protective immunity.

Enhanced Assessment of Cross-Reactive Antigenic Determinants within the Spike Protein

Despite successful vaccination efforts, the emergence of new SARS-CoV-2 variants poses ongoing challenges to control COVID-19. Understanding humoral responses regarding SARS-CoV-2 infections and their impact is crucial for developing future vaccines that are effective worldwide. Here, we identified 41 immunodominant linear B-cell epitopes in its spike glycoprotein with an SPOT synthesis peptide array probed with a pool of serum from hospitalized COVID-19 patients. The bioinformatics showed a restricted set of epitopes unique to SARS-CoV-2 compared to other coronavirus family members. Potential crosstalk was also detected with Dengue virus (DENV), which was confirmed by screening individuals infected with DENV before the COVID-19 pandemic in a commercial ELISA for anti-SARS-CoV-2 antibodies. A high-resolution evaluation of antibody reactivity against peptides representing epitopes in the spike protein identified ten sequences in the NTD, RBD, and S2 domains. Functionally, antibody-dependent enhancement (ADE) in SARS-CoV-2 infections of monocytes was observed in vitro with pre-pandemic Dengue-positive sera. A significant increase in viral load was measured compared to that of the controls, with no detectable neutralization or considerable cell death, suggesting its role in viral entry. Cross-reactivity against peptides from spike proteins was observed for the pre-pandemic sera. This study highlights the importance of identifying specific epitopes generated during the humoral response to a pathogenic infection to understand the potential interplay of previous and future infections on diseases and their impact on vaccinations and immunodiagnostics.

Plasmonic Imaging of Multivalent NTD-Nucleic Acid Interactions for Broad-Spectrum Antiviral Drug Analysis

The discovery and identification of broad-spectrum antiviral drugs are of great significance for blocking the spread of pathogenic viruses and corresponding variants of concern. Herein, we proposed a plasmonic imaging-based strategy for assessing the efficacy of potential broad-spectrum antiviral drugs targeting the N-terminal domain of a nucleocapsid protein (NTD) and nucleic acid (NA) interactions. With NTD and NA conjugated gold nanoparticles as core and satellite nanoprobes, respectively, we found that the multivalent binding interactions could drive the formation of core-satellite nanostructures with enhanced scattering brightness due to the plasmonic coupling effect. The core-satellite assembly can be suppressed in the presence of antiviral drugs targeting the NTD-NA interactions, allowing the drug efficacy analysis by detecting the dose-dependent changes in the scattering brightness by plasmonic imaging. By quantifying the changes in the scattering brightness of plasmonic nanoprobes, we uncovered that the constructed multivalent weak interactions displayed a 500-fold enhancement in affinity as compared with the monovalent NTD-NA interactions. We demonstrated the plasmonic imaging-based strategy for evaluating the efficacy of a potential broad-spectrum drug, PJ34, that can target the NTD-NA interactions, with the IC50 as 24.35 and 14.64 μM for SARS-CoV-2 and SARS-CoV, respectively. Moreover, we discovered that ceftazidime holds the potential as a candidate drug to inhibit the NTD-NA interactions with an IC50 of 22.08 μM from molecular docking and plasmonic imaging-based drug analysis. Finally, we validated that the potential antiviral drug, 5-benzyloxygramine, which can induce the abnormal dimerization of nucleocapsid proteins, is effective for SARS-CoV-2, but not effective against SARS-CoV. All these demonstrations indicated that the plasmonic imaging-based strategy is robust and can be used as a powerful strategy for the discovery and identification of broad-spectrum drugs targeting the evolutionarily conserved viral proteins.

A new era of antibody discovery: an in-depth review of AI-driven approaches

Given their high affinity and specificity for a range of macromolecules, antibodies are widely used in the treatment of autoimmune diseases, cancers, inflammatory diseases, and Alzheimer's disease (AD). Traditional experimental methods are time-consuming, expensive, and labor-intensive. Recent advances in artificial intelligence (AI) technologies provide complementary methods that can reduce the time and costs required for antibody design by minimizing failures and increasing the success rate of experimental tests. In this review, we scrutinize the plethora of AI-driven methodologies that have been deployed over the past 4 years for modeling antibody structures, predicting antibody-antigen interactions, optimizing antibody affinity, and generating novel antibody candidates. We also briefly address the challenges faced in integrating AI-based models with traditional antibody discovery pipelines and highlight the potential future directions in this burgeoning field.

Loss to gain: pseudogenes in microorganisms, focusing on eubacteria, and their biological significance

Pseudogenes are defined as "non-functional" copies of corresponding parent genes. The cognition of pseudogenes continues to be refreshed through accumulating and updating research findings. Previous studies have predominantly focused on mammals, but pseudogenes have received relatively less attention in the field of microbiology. Given the increasing recognition on the importance of pseudogenes, in this review, we focus on several aspects of microorganism pseudogenes, including their classification and characteristics, their generation and fate, their identification, their abundance and distribution, their impact on virulence, their ability to recombine with functional genes, the extent to which some pseudogenes are transcribed and translated, and the relationship between pseudogenes and viruses. By summarizing and organizing the latest research progress, this review will provide a comprehensive perspective and improved understanding on pseudogenes in microorganisms. KEY POINTS: • Concept, classification and characteristics, identification and databases, content, and distribution of microbial pseudogenes are presented. • How pseudogenization contribute to pathogen virulence is highlighted. • Pseudogenes with potential functions in microorganisms are discussed.

Comparison of SARS-CoV-2 variants of concern in primary human nasal cultures demonstrates Delta as most cytopathic and Omicron as fastest replicating

The SARS-CoV-2 pandemic was marked with emerging viral variants, some of which were designated as variants of concern (VOCs) due to selection and rapid circulation in the human population. Here, we elucidate functional features of each VOC linked to variations in replication rate. Patient-derived primary nasal cultures grown at air-liquid interface were used to model upper respiratory infection and compared to cell lines derived from human lung epithelia. All VOCs replicated to higher titers than the ancestral virus, suggesting a selection for replication efficiency. In primary nasal cultures, Omicron replicated to the highest titers at early time points, followed by Delta, paralleling comparative studies of population sampling. All SARS-CoV-2 viruses entered the cell primarily via a transmembrane serine protease 2 (TMPRSS2)-dependent pathway, and Omicron was more likely to use an endosomal route of entry. All VOCs activated and overcame dsRNA-induced cellular responses, including interferon (IFN) signaling, oligoadenylate ribonuclease L degradation, and protein kinase R activation. Among the VOCs, Omicron infection induced expression of the most IFN and IFN-stimulated genes. Infections in nasal cultures resulted in cellular damage, including a compromise of cell barrier integrity and loss of nasal cilia and ciliary beating function, especially during Delta infection. Overall, Omicron was optimized for replication in the upper respiratory tract and least favorable in the lower respiratory cell line, and Delta was the most cytopathic for both upper and lower respiratory cells. Our findings highlight the functional differences among VOCs at the cellular level and imply distinct mechanisms of pathogenesis in infected individuals.

IMPORTANCE

Comparative analysis of infections by SARS-CoV-2 ancestral virus and variants of concern, including Alpha, Beta, Delta, and Omicron, indicated that variants were selected for efficiency in replication. In infections of patient-derived primary nasal cultures grown at air-liquid interface to model upper respiratory infection, Omicron reached the highest titers at early time points, a finding that was confirmed by parallel population sampling studies. While all infections overcame dsRNA-mediated host responses, infections with Omicron induced the strongest interferon and interferon-stimulated gene response. In both primary nasal cultures and lower respiratory cell line, infections by Delta were most damaging to the cells as indicated by syncytia formation, loss of cell barrier integrity, and nasal ciliary function.

Altered receptor binding, antibody evasion and retention of T cell recognition by the SARS-CoV-2 XBB.1.5 spike protein

The XBB.1.5 variant of SARS-CoV-2 has rapidly achieved global dominance and exhibits a high growth advantage over previous variants. Preliminary reports suggest that the success of XBB.1.5 stems from mutations within its spike glycoprotein, causing immune evasion and enhanced receptor binding. We present receptor binding studies that demonstrate retention of binding contacts with the human ACE2 receptor and a striking decrease in binding to mouse ACE2 due to the revertant R493Q mutation. Despite extensive evasion of antibody binding, we highlight a region on the XBB.1.5 spike protein receptor binding domain (RBD) that is recognized by serum antibodies from a donor with hybrid immunity, collected prior to the emergence of the XBB.1.5 variant. T cell assays reveal high frequencies of XBB.1.5 spike-specific CD4+ and CD8+ T cells amongst donors with hybrid immunity, with the CD4+ T cells skewed towards a Th1 cell phenotype and having attenuated effector cytokine secretion as compared to ancestral spike protein-specific cells. Thus, while the XBB.1.5 variant has retained efficient human receptor binding and gained antigenic alterations, it remains susceptible to recognition by T cells induced via vaccination and previous infection.

Differential laboratory passaging of SARS-CoV-2 viral stocks impacts the in vitro assessment of neutralizing antibodies

Viral populations in natural infections can have a high degree of sequence diversity, which can directly impact immune escape. However, antibody potency is often tested in vitro with a relatively clonal viral populations, such as laboratory virus or pseudotyped virus stocks, which may not accurately represent the genetic diversity of circulating viral genotypes. This can affect the validity of viral phenotype assays, such as antibody neutralization assays. To address this issue, we tested whether recombinant virus carrying SARS-CoV-2 spike (VSV-SARS-CoV-2-S) stocks could be made more genetically diverse by passage, and if a stock passaged under selective pressure was more capable of escaping monoclonal antibody (mAb) neutralization than unpassaged stock or than viral stock passaged without selective pressures. We passaged VSV-SARS-CoV-2-S four times concurrently in three cell lines and then six times with or without polyclonal antiserum selection pressure. All three of the monoclonal antibodies tested neutralized the viral population present in the unpassaged stock. The viral inoculum derived from serial passage without antiserum selection pressure was neutralized by two of the three mAbs. However, the viral inoculum derived from serial passage under antiserum selection pressure escaped neutralization by all three mAbs. Deep sequencing revealed the rapid acquisition of multiple mutations associated with antibody escape in the VSV-SARS-CoV-2-S that had been passaged in the presence of antiserum, including key mutations present in currently circulating Omicron subvariants. These data indicate that viral stock that was generated under polyclonal antiserum selection pressure better reflects the natural environment of the circulating virus and may yield more biologically relevant outcomes in phenotypic assays. Thus, mAb assessment assays that utilize a more genetically diverse, biologically relevant, virus stock may yield data that are relevant for prediction of mAb efficacy and for enhancing biosurveillance.

Significance of understanding the genomics of host-pathogen interaction in limiting antibiotic resistance development: lessons from COVID-19 pandemic

The entire world is facing the stiff challenge of COVID-19 pandemic. To overcome the spread of this highly infectious disease, several short-sighted strategies were adopted such as the use of broad-spectrum antibiotics and antifungals. However, the misuse and/or overuse of antibiotics have accentuated the emergence of the next pandemic: antimicrobial resistance (AMR). It is believed that pathogens while transferring between humans and the environment carry virulence and antibiotic-resistant factors from varied species. It is presumed that all such genetic factors are quantifiable and predictable, a better understanding of which could be a limiting step for the progression of AMR. Herein, we have reviewed how genomics-based understanding of host-pathogen interactions during COVID-19 could reduce the non-judicial use of antibiotics and prevent the eruption of an AMR-based pandemic in future.

A multifaceted approach for identification, validation, and immunogenicity of naturally processed and in silico-predicted highly conserved SARS-CoV-2 peptides

SARS-CoV-2 remains a major global public health concern. Antibody waning and immune escape variant emergence necessitate the development of next generation vaccines that induce cross-reactive durable immune responses. T cell responses to SARS-CoV-2 demonstrate higher conservation, antigenic breadth, and longevity than antibody responses. Therefore, we sought to identify pathogen-derived T cell epitopes for a potential peptide-based vaccine. We pursued an approach leveraging: 1) liquid chromatography and tandem mass spectrometry (LC-MS/MS)-based identification of peptides from ancestral SARS-CoV-2-infected cell lines, 2) epitope prediction algorithms, and 3) overlapping peptide libraries. From this strategy, we identified 380 unique SARS-CoV-2-derived peptide sequences, including 53 antigenic HLA class I and class II peptides from multiple structural and non-structural/accessory viral proteins. These peptide sequences were highly conserved across variants of concern/interest (VoC/VoIs), and are estimated to achieve coverage of >96% of the world population. Our findings validate this discovery pipeline for peptide identification and immunogenicity testing, and are a crucial step toward the development of a next-generation multi-epitope SARS-CoV-2 peptide vaccine, and a novel vaccine platform methodology.

Evolution of the SARS-CoV-2 Omicron spike

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant of concern, first identified in November 2021, rapidly spread worldwide and diversified into several subvariants. The Omicron spike (S) protein accumulated an unprecedented number of sequence changes relative to previous variants. In this review, we discuss how Omicron S protein structural features modulate host cell receptor binding, virus entry, and immune evasion and highlight how these structural features differentiate Omicron from previous variants. We also examine how key structural properties track across the still-evolving Omicron subvariants and the importance of continuing surveillance of the S protein sequence evolution over time.

Comparison of SARS-CoV-2 variants of concern in primary human nasal cultures demonstrates Delta as most cytopathic and Omicron as fastest replicating

The SARS-CoV-2 pandemic was marked with emerging viral variants, some of which were designated as variants of concern (VOCs) due to selection and rapid circulation in the human population. Here we elucidate functional features of each VOC linked to variations in replication rate. Patient-derived primary nasal cultures grown at air-liquid-interface (ALI) were used to model upper-respiratory infection and human lung epithelial cell lines used to model lower-respiratory infection. All VOCs replicated to higher titers than the ancestral virus, suggesting a selection for replication efficiency. In primary nasal cultures, Omicron replicated to the highest titers at early time points, followed by Delta, paralleling comparative studies of population sampling. All SARS-CoV-2 viruses entered the cell primarily via a transmembrane serine protease 2 (TMPRSS2)-dependent pathway, and Omicron was more likely to use an endosomal route of entry. All VOCs activated and overcame dsRNA-induced cellular responses including interferon (IFN) signaling, oligoadenylate ribonuclease L degradation and protein kinase R activation. Among the VOCs, Omicron infection induced expression of the most IFN and IFN stimulated genes. Infections in nasal cultures resulted in cellular damage, including a compromise of cell-barrier integrity and loss of nasal cilia and ciliary beating function, especially during Delta infection. Overall, Omicron was optimized for replication in the upper-respiratory system and least-favorable in the lower-respiratory cell line; and Delta was the most cytopathic for both upper and lower respiratory cells. Our findings highlight the functional differences among VOCs at the cellular level and imply distinct mechanisms of pathogenesis in infected individuals.

Wireless, battery-free, multifunctional integrated bioelectronics for respiratory pathogens monitoring and severity evaluation

The rapid diagnosis of respiratory virus infection through breath and blow remains challenging. Here we develop a wireless, battery-free, multifunctional pathogenic infection diagnosis system (PIDS) for diagnosing SARS-CoV-2 infection and symptom severity by blow and breath within 110 s and 350 s, respectively. The accuracies reach to 100% and 92% for evaluating the infection and symptom severity of 42 participants, respectively. PIDS realizes simultaneous gaseous sample collection, biomarker identification, abnormal physical signs recording and machine learning analysis. We transform PIDS into other miniaturized wearable or portable electronic platforms that may widen the diagnostic modes at home, outdoors and public places. Collectively, we demonstrate a general-purpose technology for rapidly diagnosing respiratory pathogenic infection by breath and blow, alleviating the technical bottleneck of saliva and nasopharyngeal secretions. PIDS may serve as a complementary diagnostic tool for other point-of-care techniques and guide the symptomatic treatment of viral infections.

Chronicling the 3-year evolution of the COVID-19 pandemic: analysis of disease management, characteristics of major variants, and impacts on pathogenicity

Announced on December 31, 2019, the novel coronavirus arising in Wuhan City, Hubei Province resulted in millions of cases and lives lost. Following intense tracking, coronavirus disease 2019 (COVID-19) was declared a pandemic by the World Health Organization (WHO) in 2020. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was identified as the cause of COVID-19 and the continuous evolution of the virus has given rise to several variants. In this review, a comprehensive analysis of the response to the pandemic over the first three-year period is provided, focusing on disease management, development of vaccines and therapeutics, and identification of variants. The transmissibility and pathogenicity of SARS-CoV-2 variants including Alpha, Beta, Gamma, Delta, and Omicron are compared. The binding characteristics of the SARS-CoV-2 spike protein to the angiotensin-converting enzyme 2 (ACE2) receptor and reproduction numbers are evaluated. The effects of major variants on disease severity, hospitalisation, and case-fatality rates are outlined. In addition to the spike protein, open reading frames mutations are investigated. We also compare the pathogenicity of SARS-CoV-2 with SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV). Overall, this study highlights the strengths and weaknesses of the global response to the pandemic, as well as the importance of prevention and preparedness. Monitoring the evolution of SARS-CoV-2 is critical in identifying and potentially predicting the health outcomes of concerning variants as they emerge. The ultimate goal would be a position in which existing vaccines and therapeutics could be adapted to suit new variants in as close to real-time as possible.

Spike Protein Mutation-Induced Changes in the Kinetic and Thermodynamic Behavior of Its Receptor Binding Domains Explain Their Higher Propensity to Attain Open States in SARS-CoV-2 Variants of Concern

Spike (S) protein opening in SARS-CoV-2 controls the accessibility of its receptor binding domains (RBDs) to host receptors and immune recognition. Along the evolution of SARS-CoV-2 to its variants of concern (VOC)-alpha, beta, gamma, delta, and omicron-their S proteins showed a higher propensity to attain open states. Deciphering how mutations in S protein can shape its conformational dynamics will contribute to the understanding of viral host tropism. Here using microsecond-scale multiple molecular dynamics simulations (MDS), we provide insights into the kinetic and thermodynamic contributions of these mutations to RBD opening pathways in S proteins of SARS-CoV-2 VOCs. Mutational effects were analyzed using atomistic (i) equilibrium MDS of closed and open states of S proteins and (ii) nonequilibrium MDS for closed-to-open transitions. In MDS of closed or open states, RBDs in S proteins of VOCs showed lower thermodynamic stability with higher kinetic fluctuations, compared to S proteins of ancestral SARS-CoV-2. For closed-to-open transitions in S proteins of VOCs, we observed apparently faster RBD opening with a 1.5-2-fold decrease in the thermodynamic free-energy barrier (ΔGclosed→open). Saturation mutagenesis studies highlighted S protein mutations that may control its conformational dynamics and presentation to host receptors.

Dysregulation of intracellular redox homeostasis by the SARS-CoV-2 ORF6 protein

SARS-CoV-2 has evolved several strategies to overcome host cell defenses by inducing cell injury to favour its replication. Many viruses have been reported to modulate the intracellular redox balance, affecting the Nuclear factor erythroid 2-Related Factor 2 (NRF2) signaling pathway. Although antioxidant modulation by SARS-CoV-2 infection has already been described, the viral factors involved in modulating the NRF2 pathway are still elusive. Given the antagonistic activity of ORF6 on several cellular pathways, we investigated the role of the viral protein towards NRF2-mediated antioxidant response. The ectopic expression of the wt-ORF6 protein negatively impacts redox cell homeostasis, leading to an increase in ROS production, along with a decrease in NRF2 protein and its downstream controlled genes. Moreover, when investigating the Δ61 mutant, previously described as an inactive nucleopore proteins binding mutant, we prove that the oxidative stress induced by ORF6 is substantially related to its C-terminal domain, speculating that ORF6 mechanism of action is associated with the inhibition of nuclear mRNA export processes. In addition, activation by phosphorylation of the serine residue at position 40 of NRF2 is increased in the cytoplasm of wt-ORF6-expressing cells, supporting the presence of an altered redox state, although NRF2 nuclear translocation is hindered by the viral protein to fully antagonize the cell response. Furthermore, wt-ORF6 leads to phosphorylation of a stress-activated serine/threonine protein kinase, p38 MAPK, suggesting a role of the viral protein in regulating p38 activation. These findings strengthen the important role of oxidative stress in the pathogenesis of SARS-CoV-2 and identify ORF6 as an important viral accessory protein hypothetically involved in modulating the antioxidant response during viral infection.

Intrinsically Disordered Regions Function as a Cervical Collar to Remotely Regulate the Nodding Dynamics of SARS-CoV-2 Prefusion Spike Heads

The SARS-CoV-2 prefusion spike heads (receptor binding domains, RBDs) frequently nod down and up to interact with host cell receptors. As the spike protein is a trimeric unit of significant size, to understand its large-scale structural dynamics associated with the nodding mechanism and the mutational impact on the same, we develop a topological symmetry-information-loaded coarse-grained structure-based model of a spike trimer using recent cryo-EM structural data. Our study reveals the control of two distant intrinsically disordered regions (IDRs), namely, 630 and FPPR loops, over the nodding dynamics of spike heads. We find that the order-disorder transition of IDRs becomes more evident in the variants of concern (VOCs) that are associated with the characteristic mutation, D614G, in the proximity of these IDRs. In some VOCs, the two other mutations A570D and S982A also show an integral effect. The driver mutation D614G instigates a salt-bridge disruption, altering the order-disorder dynamics of both 630 and FPPR loops and their interaction with the C-terminal domains (CTD1/CTD2). This altered connectivity in these mutants allows the two IDRs to act collectively as a "cervical collar" for the RBD, supporting various spike head postures, consistent with cryo-EM results available for specific cases. The IDRs' control over the spike structure and dynamics presents an exciting opportunity where they can be targeted as remote operational switches to artificially maneuver the nod for effective therapeutic interventions.

Cross-protection and cross-neutralization capacity of ancestral and VOC-matched SARS-CoV-2 adenoviral vector-based vaccines

COVID-19 vaccines were originally designed based on the ancestral Spike protein, but immune escape of emergent Variants of Concern (VOC) jeopardized their efficacy, warranting variant-proof vaccines. Here, we used preclinical rodent models to establish the cross-protective and cross-neutralizing capacity of adenoviral-vectored vaccines expressing VOC-matched Spike. CoroVaxG.3-D.FR, matched to Delta Plus Spike, displayed the highest levels of nAb to the matched VOC and mismatched variants. Cross-protection against viral infection in aged K18-hACE2 mice showed dramatic differences among the different vaccines. While Delta-targeted vaccines fully protected mice from a challenge with Gamma, a Gamma-based vaccine offered only partial protection to Delta challenge. Administration of CorovaxG.3-D.FR in a prime/boost regimen showed that a booster was able to increase the neutralizing capacity of the sera against all variants and fully protect aged K18-hACE2 mice against Omicron BA.1, as a BA.1-targeted vaccine did. The neutralizing capacity of the sera diminished in all cases against Omicron BA.2 and BA.5. Altogether, the data demonstrate that a booster with a vaccine based on an antigenically distant variant, such as Delta or BA.1, has the potential to protect from a wider range of SARS-CoV-2 lineages, although careful surveillance of breakthrough infections will help to evaluate combination vaccines targeting antigenically divergent variants yet to emerge.

Interaction of SARS-CoV-2 with host cells and antibodies: experiment and simulation

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the devastating global COVID-19 pandemic announced by WHO in March 2020. Through unprecedented scientific effort, several vaccines, drugs and antibodies have been developed, saving millions of lives, but the fight against COVID-19 continues as immune escape variants of concern such as Delta and Omicron emerge. To develop more effective treatments and to elucidate the side effects caused by vaccines and therapeutic agents, a deeper understanding of the molecular interactions of SARS-CoV-2 with them and human cells is required. With special interest in computational approaches, we will focus on the structure of SARS-CoV-2 and the interaction of its spike protein with human angiotensin-converting enzyme-2 (ACE2) as a prime entry point of the virus into host cells. In addition, other possible viral receptors will be considered. The fusion of viral and human membranes and the interaction of the spike protein with antibodies and nanobodies will be discussed, as well as the effect of SARS-CoV-2 on protein synthesis in host cells.

Whole genome CRISPR screening strategy to identify genes contributing to SARS-CoV-2 spike and VSV-G mediated entry

Understanding the cellular host factors that promote and inhibit viral entry is important for identifying viral countermeasures. CRISPR whole-genome screens can be used to rapidly discover host factors that contribute to or impair viral entry. However, when using live viruses and cellular lethality for selection, these screens can identify an overwhelming number of genes without specificity for the stage of the viral infection cycle. New screening methods are needed to identify host machinery contributing to specific steps of viral infection. Here, we developed a CRISPR whole-genome screen and counter-screen strategy based on a pseudoviral platform that allowed identification of genes specific to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike and vesicular stomatitis virus glycoprotein (VSV-G) mediated entry. Screening of SARS-CoV-2 spike and VSV-G on the same lentiviral pseudovirus allowed the identification of entry-specific genes relative to genes associated with retro-transcription, integration, and reporter expression from the lentiviral pseudovirus. Second, a Cre-Gag fusion protein packaged into the pseudovirus was used to bypass retro-transcription and integration by directly activating a floxed fluorescent protein reporter upon entry  reduced the number of gene hits and increase specificity for viral entry. Our approach correctly identified SARS-CoV-2 and VSV-G receptors ACE2 and low-density lipoprotein receptors, respectively, and distinguished genes associated with retroviral reporter expression from envelope-mediated entry. Moreover, the CRE-Gag fusion/flox reporter increased the screen specificity for viral entry-associated genes. Validation of a few hits demonstrates that this approach distinguishes envelope-specific host factors from genes affecting reporter expression. Overall, this approach provides a new strategy for identifying host genes influencing viral entry without the confounding complexity of live-viral screens which produce long gene lists associated with all aspects of viral pathogenesis and replication.

Persistent topological Laplacian analysis of SARS-CoV-2 variants

Topological data analysis (TDA) is an emerging field in mathematics and data science. Its central technique, persistent homology, has had tremendous success in many science and engineering disciplines. However, persistent homology has limitations, including its inability to handle heterogeneous information, such as multiple types of geometric objects; being qualitative rather than quantitative, e.g., counting a 5-member ring the same as a 6-member ring, and a failure to describe non-topological changes, such as homotopic changes in protein-protein binding. Persistent topological Laplacians (PTLs), such as persistent Laplacian and persistent sheaf Laplacian, were proposed to overcome the limitations of persistent homology. In this work, we examine the modeling and analysis power of PTLs in the study of the protein structures of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike receptor binding domain (RBD). First, we employ PTLs to study how the RBD mutation-induced structural changes of RBD-angiotensin-converting enzyme 2 (ACE2) binding complexes are captured in the changes of spectra of the PTLs among SARS-CoV-2 variants. Additionally, we use PTLs to analyze the binding of RBD and ACE2-induced structural changes of various SARS-CoV-2 variants. Finally, we explore the impacts of computationally generated RBD structures on a topological deep learning paradigm and predictions of deep mutational scanning datasets for the SARS-CoV-2 Omicron BA.2 variant. Our results indicate that PTLs have advantages over persistent homology in analyzing protein structural changes and provide a powerful new TDA tool for data science.

Analysis of the N-glycosylation profiles of the spike proteins from the Alpha, Beta, Gamma, and Delta variants of SARS-CoV-2

N-Glycosylation plays an important role in the structure and function of membrane and secreted proteins. Viral proteins used in cell entry are often extensively glycosylated to assist in protein folding, provide stability, and shield the virus from immune recognition by its host (described as a "glycan shield"). The SARS-CoV-2 spike protein (S) is a prime example, having 22 potential sites of N-glycosylation per protein protomer, as predicted from the primary sequence. In this report, we conducted mass spectrometric analysis of the N-glycosylation profiles of recombinant spike proteins derived from four common SARS-CoV-2 variants classified as Variant of Concern, including Alpha, Beta, Gamma, and Delta along with D614G variant spike as a control. Our data reveal that the amino acid substitutions and deletions between variants impact the abundance and type of glycans on glycosylation sites of the spike protein. Some of the N-glycosylation sequons in S show differences between SARS-CoV-2 variants in the distribution of glycan forms. In comparison with our previously reported site-specific glycan analysis on the S-D614G and its ancestral protein, glycan types on later variants showed high similarity on the site-specific glycan content to S-D614G. Additionally, we applied multiple digestion methods on each sample, and confirmed the results for individual glycosylation sites from different experiment conditions to improve the identification and quantification of glycopeptides. Detailed site-specific glycan analysis of a wide variety of SARS-CoV-2 variants provides useful information toward the understanding of the role of protein glycosylation on viral protein structure and function and development of effective vaccines and therapeutics.

The SARS-CoV-2 spike protein structure: a symmetry tale on distortion trail

A preliminary step in the SARS-CoV-2 human infection process is a conformational change of the receptor binding domain (RBD) of its spike protein, characterized by a significant loss of symmetry. During this process, the residues which later on bind to the human angiotensin converting enzyme 2 (ACE2) receptor, become exposed at the surface of the protein. Symmetry analysis of a data set of 33 protein structures from experimental measurements and 32 structures from molecular dynamics simulation, show that the initial state carries clear indications on the structure of the final state, with respect to the local distortion along the sequence. This surprising finding implies that this type of analysis predicts the mechanism of change. We further show that the level of local distortion at the initial state increases with variant's transmissibility, for the wild type (WT) along with past and present variants of concern (WT ∼ alpha < beta < delta < Omicron BA.1), in accordance with the trend of their evolutionary path. In other words, the initial structure of the variant which is most infectious is also the most distorted, making its path to the final state shorter. It has been claimed that the RBD migration of the spike protein is allosterically controlled. Our analysis provides a quantitative support to a major theorem in this respect - that information about an allosteric process is encoded in the structure itself, suggesting that the path of local distortion is related to an allosteric information network.

Next-generation nanophotonic-enabled biosensors for intelligent diagnosis of SARS-CoV-2 variants

Constantly mutating SARS-CoV-2 is a global concern resulting in COVID-19 infectious waves from time to time in different regions, challenging present-day diagnostics and therapeutics. Early-stage point-of-care diagnostic (POC) biosensors are a crucial vector for the timely management of morbidity and mortalities caused due to COVID-19. The state-of-the-art SARS-CoV-2 biosensors depend upon developing a single platform for its diverse variants/biomarkers, enabling precise detection and monitoring. Nanophotonic-enabled biosensors have emerged as 'one platform' to diagnose COVID-19, addressing the concern of constant viral mutation. This review assesses the evolution of current and future variants of the SARS-CoV-2 and critically summarizes the current state of biosensor approaches for detecting SARS-CoV-2 variants/biomarkers employing nanophotonic-enabled diagnostics. It discusses the integration of modern-age technologies, including artificial intelligence, machine learning and 5G communication with nanophotonic biosensors for intelligent COVID-19 monitoring and management. It also highlights the challenges and potential opportunities for developing intelligent biosensors for diagnosing future SARS-CoV-2 variants. This review will guide future research and development on nano-enabled intelligent photonic-biosensor strategies for early-stage diagnosing of highly infectious diseases to prevent repeated outbreaks and save associated human mortalities.

Distribution and Functional Analyses of Mutations in Spike Protein and Phylogenic Diversity of SARS-CoV-2 Variants Emerged during the Year 2021 in India

Introduction

Prolonged COVID-19 pandemic accelerates the emergence and transmissibility of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) variants through the accumulation of adaptive mutations. Particularly, adaptive mutations in spike (S) protein of SARS-CoV-2 leads to increased viral infectivity, severe morbidity and mortality, and immune evasion. This study focuses on the phylodynamic distribution of SARS-CoV-2 variants during the year 2021 in India besides analyzing the functional significance of mutations in S-protein of SARS-CoV-2 variants.

Methods

Whole genome of SARS-CoV-2 sequences (n = 87957) from the various parts of India over the period of January to December 2021 was retrieved from Global Initiative on Sharing All Influenza Data. All the S-protein sequences were subjected to clade analysis, variant calling, protein stability, immune escape potential, structural divergence, Furin cleavage efficiency, and phylogenetic analysis using various in silico tools.

Results

Delta variant belonging to 21A, 21I, and 21J clades was found to be predominant throughout the year 2021 though many variants were also present. A total of 4639 amino acid mutations were found in S-protein. D614G was the most predominant mutation in the S-protein followed by P681R, L452R, T19R, T478K, and D950N. The highest number of mutations was found in the N-terminal domain of S-protein. Mutations in the crucial sites of S-protein impacting pathogenicity, immunogenicity, and fusogenicity were identified. Intralineage diversity analysis showed that certain variants of SARS-CoV-2 possess high diversification.

Conclusions

The study has disclosed the distribution of various variants including the Delta, the predominant variant, in India throughout the year 2021. The study has identified mutations in S-protein of each SARS-CoV-2 variant that can significantly impact the virulence, immune evasion, increased transmissibility, high morbidity, and mortality. In addition, it is found that mutations acquired during each viral replication cycle introduce new sub-lineages as studied by intralineage diversity analysis.

RBD and Spike DNA-Based Immunization in Rabbits Elicited IgG Avidity Maturation and High Neutralizing Antibody Responses against SARS-CoV-2

Neutralizing antibodies (nAbs) are a critical part of coronavirus disease 2019 (COVID-19) research as they are used to gain insight into the immune response to severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) infections. Among the technologies available for generating nAbs, DNA-based immunization methods are an alternative to conventional protocols. In this pilot study, we investigated whether DNA-based immunization by needle injection in rabbits was a viable approach to produce a functional antibody response. We demonstrated that three doses of DNA plasmid carrying the gene encoding the full-length spike protein (S) or the receptor binding domain (RBD) of SARS-CoV-2 induced a time-dependent increase in IgG antibody avidity maturation. Moreover, the IgG antibodies displayed high cross neutralization by live SARS-CoV-2 and pseudoviruses neutralization assays. Thus, we established a simple, low cost and feasible DNA-based immunization protocol in rabbits that elicited high IgG avidity maturation and nAbs production against SARS-CoV-2, highlighting the importance of DNA-based platforms for developing new immunization strategies against SARS-CoV-2 and future emerging epidemics.

Structural analysis of receptor engagement and antigenic drift within the BA.2 spike protein

The BA.2 sub-lineage of the Omicron (B.1.1.529) severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant rapidly supplanted the original BA.1 sub-lineage in early 2022. Both lineages threatened the efficacy of vaccine-elicited antibodies and acquired increased binding to several mammalian ACE2 receptors. Cryoelectron microscopy (cryo-EM) analysis of the BA.2 spike (S) glycoprotein in complex with mouse ACE2 (mACE2) identifies BA.1- and BA.2-mutated residues Q493R, N501Y, and Y505H as complementing non-conserved residues between human and mouse ACE2, rationalizing the enhanced S protein-mACE2 interaction for Omicron variants. Cryo-EM structures of the BA.2 S-human ACE2 complex and of the extensively mutated BA.2 amino-terminal domain (NTD) reveal a dramatic reorganization of the highly antigenic N1 loop into a β-strand, providing an explanation for decreased binding of the BA.2 S protein to antibodies isolated from BA.1-convalescent patients. Our analysis reveals structural mechanisms underlying the antigenic drift in the rapidly evolving Omicron variant landscape.

Mutational Patterns Observed in SARS-CoV-2 Genomes Sampled From Successive Epochs Delimited by Major Public Health Events in Ontario, Canada: Genomic Surveillance Study

BACKGROUND

The emergence of SARS-CoV-2 variants with mutations associated with increased transmissibility and virulence is a public health concern in Ontario, Canada. Characterizing how the mutational patterns of the SARS-CoV-2 genome have changed over time can shed light on the driving factors, including selection for increased fitness and host immune response, that may contribute to the emergence of novel variants. Moreover, the study of SARS-CoV-2 in the microcosm of Ontario, Canada can reveal how different province-specific public health policies over time may be associated with observed mutational patterns as a model system.

OBJECTIVE

This study aimed to perform a comprehensive analysis of single base substitution (SBS) types, counts, and genomic locations observed in SARS-CoV-2 genomic sequences sampled in Ontario, Canada. Comparisons of mutational patterns were conducted between sequences sampled during 4 different epochs delimited by major public health events to track the evolution of the SARS-CoV-2 mutational landscape over 2 years.

METHODS

In total, 24,244 SARS-CoV-2 genomic sequences and associated metadata sampled in Ontario, Canada from January 1, 2020, to December 31, 2021, were retrieved from the Global Initiative on Sharing All Influenza Data database. Sequences were assigned to 4 epochs delimited by major public health events based on the sampling date. SBSs from each SARS-CoV-2 sequence were identified relative to the MN996528.1 reference genome. Catalogues of SBS types and counts were generated to estimate the impact of selection in each open reading frame, and identify mutation clusters. The estimation of mutational fitness over time was performed using the Augur pipeline.

RESULTS

The biases in SBS types and proportions observed support previous reports of host antiviral defense activity involving the SARS-CoV-2 genome. There was an increase in U>C substitutions associated with adenosine deaminase acting on RNA (ADAR) activity uniquely observed during Epoch 4. The burden of novel SBSs observed in SARS-CoV-2 genomic sequences was the greatest in Epoch 2 (median 5), followed by Epoch 3 (median 4). Clusters of SBSs were observed in the spike protein open reading frame, ORF1a, and ORF3a. The high proportion of nonsynonymous SBSs and increasing dN/dS metric (ratio of nonsynonymous to synonymous mutations in a given open reading frame) to above 1 in Epoch 4 indicate positive selection of the spike protein open reading frame.

CONCLUSIONS

Quantitative analysis of the mutational patterns of the SARS-CoV-2 genome in the microcosm of Ontario, Canada within early consecutive epochs of the pandemic tracked the mutational dynamics in the context of public health events that instigate significant shifts in selection and mutagenesis. Continued genomic surveillance of emergent variants will be useful for the design of public health policies in response to the evolving COVID-19 pandemic.

Microalgae as an Efficient Vehicle for the Production and Targeted Delivery of Therapeutic Glycoproteins against SARS-CoV-2 Variants

Severe acute respiratory syndrome–Coronavirus 2 (SARS-CoV-2) can infect various human organs, including the respiratory, circulatory, nervous, and gastrointestinal ones. The virus is internalized into human cells by binding to the human angiotensin-converting enzyme 2 (ACE2) receptor through its spike protein (S-glycoprotein). As S-glycoprotein is required for the attachment and entry into the human target cells, it is the primary mediator of SARS-CoV-2 infectivity. Currently, this glycoprotein has received considerable attention as a key component for the development of antiviral vaccines or biologics against SARS-CoV-2. Moreover, since the ACE2 receptor constitutes the main entry route for the SARS-CoV-2 virus, its soluble form could be considered as a promising approach for the treatment of coronavirus disease 2019 infection (COVID-19). Both S-glycoprotein and ACE2 are highly glycosylated molecules containing 22 and 7 consensus N-glycosylation sites, respectively. The N-glycan structures attached to these specific sites are required for the folding, conformation, recycling, and biological activity of both glycoproteins. Thus far, recombinant S-glycoprotein and ACE2 have been produced primarily in mammalian cells, which is an expensive process. Therefore, benefiting from a cheaper cell-based biofactory would be a good value added to the development of cost-effective recombinant vaccines and biopharmaceuticals directed against COVID-19. To this end, efficient protein synthesis machinery and the ability to properly impose post-translational modifications make microalgae an eco-friendly platform for the production of pharmaceutical glycoproteins. Notably, several microalgae (e.g., Chlamydomonas reinhardtii, Dunaliella bardawil, and Chlorella species) are already approved by the U.S. Food and Drug Administration (FDA) as safe human food. Because microalgal cells contain a rigid cell wall that could act as a natural encapsulation to protect the recombinant proteins from the aggressive environment of the stomach, this feature could be used for the rapid production and edible targeted delivery of S-glycoprotein and soluble ACE2 for the treatment/inhibition of SARS-CoV-2. Herein, we have reviewed the pathogenesis mechanism of SARS-CoV-2 and then highlighted the potential of microalgae for the treatment/inhibition of COVID-19 infection.