Mutations in the B.1.1.7 SARS-CoV-2 Spike Protein Reduce Receptor-Binding Affinity and Induce a Flexible Link to the Fusion Peptide

SARS-CoV-2 Spike Protein Induces Time-Dependent CTSL Upregulation in HeLa Cells and Alveolarspheres

Autophagy and lysosomal pathways are involved in the cell entry of SARS-CoV-2 virus. To infect the host cell, the spike protein of SARS-CoV-2 binds to the cell surface receptor angiotensin-converting enzyme 2 (ACE2). To allow the fusion of the viral envelope with the host cell membrane, the spike protein has to be cleaved. One possible mechanism is the endocytosis of the SARS-CoV-2-ACE2 complex and subsequent cleavage of the spike protein, mainly by the lysosomal protease cathepsin L. However, detailed molecular and dynamic insights into the role of cathepsin L in viral cell entry remain elusive. To address this, HeLa cells and iPSC-derived alveolarspheres were treated with recombinant SARS-CoV-2 spike protein, and the changes in mRNA and protein levels of cathepsins L, B, and D were monitored. Additionally, we studied the effect of cathepsin L deficiency on spike protein internalization and investigated the influence of the spike protein on cathepsin L promoters in vitro. Furthermore, we analyzed variants in the genes coding for cathepsin L, B, D, and ACE2 possibly associated with disease progression using data from Regeneron's COVID Results Browser and our own cohort of 173 patients with COVID-19, exhibiting a variant of ACE2 showing significant association with COVID-19 disease progression. Our in vitro studies revealed a significant increase in cathepsin L mRNA and protein levels following exposure to the SARS-CoV-2 spike protein in HeLa cells, accompanied by elevated mRNA levels of cathepsin B and D in alveolarspheres. Moreover, an increase in cathepsin L promoter activity was detected in vitro upon spike protein treatment. Notably, the knockout of cathepsin L resulted in reduced internalization of the spike protein. The study highlights the importance of cathepsin L and lysosomal proteases in the SARS-CoV-2 spike protein internalization and suggests the potential of lysosomal proteases as possible therapeutic targets against COVID-19 and other viral infections.

Some mechanistic underpinnings of molecular adaptations of SARS-COV-2 spike protein by integrating candidate adaptive polymorphisms with protein dynamics

We integrate evolutionary predictions based on the neutral theory of molecular evolution with protein dynamics to generate mechanistic insight into the molecular adaptations of the SARS-COV-2 spike (S) protein. With this approach, we first identified candidate adaptive polymorphisms (CAPs) of the SARS-CoV-2 S protein and assessed the impact of these CAPs through dynamics analysis. Not only have we found that CAPs frequently overlap with well-known functional sites, but also, using several different dynamics-based metrics, we reveal the critical allosteric interplay between SARS-CoV-2 CAPs and the S protein binding sites with the human ACE2 (hACE2) protein. CAPs interact far differently with the hACE2 binding site residues in the open conformation of the S protein compared to the closed form. In particular, the CAP sites control the dynamics of binding residues in the open state, suggesting an allosteric control of hACE2 binding. We also explored the characteristic mutations of different SARS-CoV-2 strains to find dynamic hallmarks and potential effects of future mutations. Our analyses reveal that Delta strain-specific variants have non-additive (i.e., epistatic) interactions with CAP sites, whereas the less pathogenic Omicron strains have mostly additive mutations. Finally, our dynamics-based analysis suggests that the novel mutations observed in the Omicron strain epistatically interact with the CAP sites to help escape antibody binding.

Function and mechanism of bispecific antibodies targeting SARS-CoV-2

As the dynamic evolution of SARS-CoV-2 led to reduced efficacy in monoclonal neutralizing antibodies and emergence of immune escape, the role of bispecific antibodies becomes crucial in bolstering antiviral activity and suppressing immune evasion. This review extensively assesses a spectrum of representative bispecific antibodies targeting SARS-CoV-2, delving into their characteristics, design formats, mechanisms of action, and associated advantages and limitations. The analysis encompasses factors influencing the selection of parental antibodies and strategies for incorporating added benefits in bispecific antibody design. Furthermore, how different classes of parental antibodies contribute to augmenting the broad-spectrum neutralization capability within bispecific antibodies is discussed. In summary, this review presents analyses and discussions aimed at offering valuable insights for shaping future strategies in bispecific antibody design to effectively confront the challenges posed by SARS-CoV-2 and propel advancements in antiviral therapeutic development.

Cooperative and structural relationships of the trimeric Spike with infectivity and antibody escape of the strains Delta (B.1.617.2) and Omicron (BA.2, BA.5, and BQ.1)

Herein, we conducted simulations of trimeric Spike from several SARS-CoV-2 variants of concern (Delta and Omicron sub-variants BA.2, BA.5, and BQ.1) and investigated the mechanisms by which specific mutations confer resistance to neutralizing antibodies. We observed that the mutations primarily affect the cooperation between protein domains within and between protomers. The substitutions K417N and L452R expand hydrogen bonding interactions, reducing their interaction with neutralizing antibodies. By interacting with nearby residues, the K444T and N460K mutations in the SpikeBQ.1 variant potentially reduces solvent exposure, thereby promoting resistance to antibodies. We also examined the impact of D614G, P681R, and P681H substitutions on Spike protein structure that may be related to infectivity. The D614G substitution influences communication between a glycine residue and neighboring domains, affecting the transition between up- and -down RBD states. The P681R mutation, found in the Delta variant, enhances correlations between protein subunits, while the P681H mutation in Omicron sub-variants weakens long-range interactions that may be associated with reduced fusogenicity. Using a multiple linear regression model, we established a connection between inter-protomer communication and loss of sensitivity to neutralizing antibodies. Our findings underscore the importance of structural communication between protein domains and provide insights into potential mechanisms of immune evasion by SARS-CoV-2. Overall, this study deepens our understanding of how specific mutations impact SARS-CoV-2 infectivity and shed light on how the virus evades the immune system.

Molecular dynamics simulations of the spike trimeric ectodomain of the SARS-CoV-2 Omicron variant: structural relationships with infectivity, evasion to immune system and transmissibility

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron is currently the most prevalent SARS-CoV-2 variant worldwide. Herein, we calculated molecular dynamics simulations of the trimeric spikeWT and SpikeBA.1 for 300 ns. Our results show that SpikeBA.1 has more conformational flexibility than SpikeWT. Our principal component analysis (PCA) allowed us to observe a broader spectrum of different conformations for SpikeBA.1, mainly at N-terminal domain (NTD) and receptor-binding domain (RBD). Such increased flexibility could contribute to decreased neutralizing antibody recognition of this variant. Our molecular dynamics data show that the RBDBA.1 easily visits an up-conformational state and the prevalent D614G mutation is pivotal to explain molecular dynamics results for this variant because to lost hydrogen bonding interactions between the residue pairs K854SC/D614SC, Y837MC/D614MC, K835SC/D614SC, T859SC/D614SC. In addition, SpikeBA.1 residues near the furin cleavage site are more flexible than in SpikeWT, probably due to P681H and D614G substitutions. Finally, dynamical cross-correlation matrix (DCCM) analysis reveals that D614G and P681H may allosterically affect the cleavage site S1/S2. Conversely, S2' site may be influenced by residues located between NTD and RBD of a neighboring protomer of the SpikeWT. Such communication may be lost in SpikeBA.1, explaining the changes of the cell tropism in the viral infection. In addition, the movements of the NTDWT and NTDBA.1 may modulate the RBD conformation through allosteric effects. Taken together, our results explain how the structural aspects may explain the observed gains in infectivity, immune system evasion and transmissibility of the Omicron variant.Communicated by Ramaswamy H. Sarma.

Genomic surveillance of genes encoding the SARS-CoV-2 spike protein to monitor for emerging variants on Jeju Island, Republic of Korea

Introduction

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has been fueled by new variants emerging from circulating strains. Here, we report results from a genomic surveillance study of SARS-CoV-2 on Jeju Island, Republic of Korea, from February 2021 to September 2022.

Methods

A total of 3,585 SARS-CoV-2 positive samples were analyzed by Sanger sequencing of the gene encoding the spike protein before performing phylogenetic analyses.

Results

We found that the Alpha variant (B.1.1.7) was dominant in May 2021 before being replaced by the Delta variant (B.1.617.2) in July 2021, which was dominant until December 2021 before being replaced by the Omicron variant. Mutations in the spike protein, including N440K and G446S, have been proposed to contribute to immune evasion, accelerating the spread of Omicron variants.

Discussion

Our results from Juju Island, Republic of Korea, are consistent with and contribute to global surveillance efforts crucial for identifying new variants of concern of SARS-CoV-2 and for monitoring the transmission dynamics and characteristics of known strains.

Umbelliferone and eriodictyol suppress the cellular entry of SARS-CoV-2

BACKGROUND

Artemisia argyi (A. argyi), also called Chinese mugwort, has been widely used to control pandemic diseases for thousands of years since ancient China due to its anti-microbial infection, anti-allergy, and anti-inflammation activities. Therefore, the potential of A. argyi and its constituents in reducing the infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was investigated in this study.

RESULTS

Among the phytochemicals in A. argyi, eriodictyol and umbelliferone were identified to target transmembrane serine protease 2 (TMPRSS2) and angiotensin-converting enzyme 2 (ACE2) proteins, the essential factors for the cellular entry of SARS-CoV-2, in both FRET-based enzymatic assays and molecular docking analyses. These two ingredients of A. argyi suppressed the infection of ACE2-expressed HEK-293 T cells with lentiviral-based pseudo-particles (Vpp) expressing wild-type and variants of SARS-CoV-2 spike (S) protein (SARS-CoV-2 S-Vpp) via interrupting the interaction between S protein and cellular receptor ACE2 and reducing the expressions of ACE2 and TMPRSS2. Oral administration with umbelliferone efficiently prevented the SARS-CoV-2 S-Vpp-induced inflammation in the lung tissues of BALB/c mice.

CONCLUSIONS

Eriodictyol and umbelliferone, the phytochemicals of Artemisia argyi, potentially suppress the cellular entry of SARS-CoV-2 by preventing the protein binding activity of the S protein to ACE2.

How helpful were molecular dynamics simulations in shaping our understanding of SARS-CoV-2 spike protein dynamics?

The SARS-CoV-2 spike protein (S) represents an important viral component that is required for successful viral infection in humans owing to its essential role in recognition of and entry to host cells. The spike is also an appealing target for drug designers who develop vaccines and antivirals. This article is important as it summarizes how molecular simulations successfully shaped our understanding of spike conformational behavior and its role in viral infection. MD simulations found that the higher affinity of SARS-CoV-2-S to ACE2 is linked to its unique residues that add extra electrostatic and van der Waal interactions in comparison to the SARS-CoV S. This illustrates the spread potential of the pandemic SARS-CoV-2 relative to the epidemic SARS-CoV. Different mutations at the S-ACE2 interface, which is believed to increase the transmission of the new variants, affected the behavior and binding interactions in different simulations. The contributions of glycans to the opening of S were revealed via simulations. The immune evasion of S was linked to the spatial distribution of glycans. This help the virus to escape the immune system recognition. This article is important as it summarizes how molecular simulations successfully shaped our understanding of spike conformational behavior and its role in viral infection. This will pave the way to us preparing for the next pandemic as the computational tools are tailored to help fight new challenges.

Potential differentiation of successive SARS-CoV-2 mutations by RNA: DNA hybrid analyses

The constant mutation of SARS-CoV-2 has triggered a new round of public health crises and has had a huge impact on existing vaccines and diagnostic tools. It is essential to develop a new flexible method to distinguish mutations to prevent the spread of the virus. In this work, we used the combination of density functional theory (DFT) and non-equilibrium Green's function formulation with decoherence, to theoretically study the effect of viral mutation on charge transport properties of viral nucleic acid molecules. We found that all mutation of SARS-CoV-2 on spike protein was accompanied by the change of gene sequence conductance, this is attributed to the change of nucleic acid molecular energy level caused by mutation. Among them, the mutations L18F, P26S, and T1027I caused the largest conductance change after mutation. This provides a theoretical possibility for detecting virus mutation based on the change of molecular conductance of virus nucleic acid.

Electrostatic Interactions Are the Primary Determinant of the Binding Affinity of SARS-CoV-2 Spike RBD to ACE2: A Computational Case Study of Omicron Variants

To explore the mechanistic origin that determines the binding affinity of SARS-CoV-2 spike receptor binding domain (RBD) to human angiotensin converting enzyme 2 (ACE2), we constructed the homology models of RBD-ACE2 complexes of four Omicron subvariants (BA.1, BA.2, BA.3 and BA.4/5), and compared them with wild type complex (RBD_WT-ACE2) in terms of various structural dynamic properties by molecular dynamics (MD) simulations and binding free energy (BFE) calculations. The results of MD simulations suggest that the RBDs of all the Omicron subvariants (RBD_OMIs) feature increased global structural fluctuations when compared with RBD_WT. Detailed comparison of BFE components reveals that the enhanced electrostatic attractive interactions are the main determinant of the higher ACE2-binding affinity of RBD_OMIs than RBD_WT, while the weakened electrostatic attractive interactions determine RBD of BA.4/5 subvariant (RBD_BA.4/5) lowest ACE2-binding affinity among all Omicron subvariants. The per-residue BFE decompositions and the hydrogen bond (HB) networks analyses indicate that the enhanced electrostatic attractive interactions are mainly through gain/loss of the positively/negatively charged residues, and the formation or destruction of the interfacial HBs and salt bridges can also largely affect the ACE2-binding affinity of RBD. It is worth pointing out that since Q493R plays the most important positive contribution in enhancing binding affinity, the absence of this mutation in RBD_BA.4/5 results in a significantly weaker binding affinity to ACE2 than other Omicron subvariants. Our results provide insight into the role of electrostatic interactions in determining of the binding affinity of SARS-CoV-2 RBD to human ACE2.

Combined effect of traditional Chinese herbal-based formulations Jing Si herbal tea and Jing Si nasal drop inhibits adhesion and transmission of SARS-CoV2 in diabetic SKH-1 mice

Multiple studies show increased severity of SARS-CoV2-infection in patients with comorbidities such as hypertension and diabetes. In this study, we have prepared two herbal-based formulations, a pleiotropic herbal drink (Jin Si Herbal Tea, JHT) and a nasal drop (Jin Si nasal drop, JND), to provide preventive care against SARS-CoV2 infection. The effect of JHT and JND was determined in SARS-CoV2-S-pseudotyped lentivirus-infected bronchial and colorectal cell lines and in SKH-1 mouse models. For preliminary studies, ACE2 receptor abundant bronchial (Calu-3) and colorectal cells (Caco-2) were used to determine the effect of JHT and JND on the host entry of various variants of SARS-CoV2-S-pseudotyped lentivirus. A series of experiments were performed to understand the infection rate in SKH-1 mice (6 weeks old, n = 9), find the effective dosage of JHT and JND, and determine the combination effect of JHT and JND on the entry and adhesion of various variant SARS-CoV2-S-pseudotyped lentiviruses, which included highly transmissible delta and gamma mutants. Furthermore, the effect of combined JHT and JND was determined on diabetes-induced SKH-1 mice against the comorbidity-associated intense viral entry and accumulation. In addition, the effect of combined JHT and JND administration on viral transmission from infected SKH-1 mice to uninfected cage mate mice was determined. The results showed that both JHT and JND were effective in alleviating the viral entry and accumulation in the thorax and the abdominal area. While JHT showed a dose-dependent decrease in the viral load, JND showed early inhibition of viral entry from day 1 of the infection. Combined administration of 48.66 mg of JHT and 20 µL of JND showed rapid reduction in the viral entry and reduced the viral load (97-99%) in the infected mice within 3 days of treatment. Moreover, 16.22 mg of JHT and 20 µL JND reduced the viral infection in STZ-induced diabetic SKH-1 mice. Interestingly, combined JHT and JND also inhibited viral transmission among cage mates. The results, therefore, showed that combined administration of JHT and JND is a novel and an efficient strategy to potentially prevent SARS-CoV2 infection.

Omicron (BA.1) and sub-variants (BA.1.1, BA.2, and BA.3) of SARS-CoV-2 spike infectivity and pathogenicity: A comparative sequence and structural-based computational assessment

The Omicron variant of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has now spread throughout the world. We used computational tools to assess the spike infectivity, transmission, and pathogenicity of Omicron (BA.1) and sub-variants (BA.1.1, BA.2, and BA.3) in this study. BA.1 has 39 mutations, BA.1.1 has 40 mutations, BA.2 has 31 mutations, and BA.3 has 34 mutations, with 21 shared mutations between all. We observed 11 common mutations in Omicron's receptor-binding domain (RBD) and sub-variants. In pathogenicity analysis, the Y505H, N786K, T95I, N211I, N856K, and V213R mutations in omicron and sub-variants are predicted to be deleterious. Due to the major effect of the mutations characterizing in the RBD, we found that Omicron and sub-variants had a higher positive electrostatic surface potential. This could increase interaction between RBD and negative electrostatic surface potential human angiotensin-converting enzyme 2 (hACE2). Omicron and sub-variants had a higher affinity for hACE2 and the potential for increased transmission when compared to the wild-type (WT). Negative electrostatic potential of N-terminal domain (NTD) of the spike protein value indicates that the Omicron variant binds receptors less efficiently than the WT. Given that at least one receptor is highly expressed in lung and bronchial cells, the electrostatic potential of NTD negative value could be one of the factors contributing to why the Omicron variant is thought to be less harmful to the lower respiratory tract. Among Omicron sub-lineages, BA.2 and BA.3 have a higher transmission potential than BA.1 and BA.1.1. We predicted that mutated residues in BA.1.1 (K478), BA.2 (R400, R490, and R495), and BA.3 (R397 and H499) formation of new salt bridges and hydrogen bonds. Omicron and sub-variant mutations at Receptor-binding Motif (RBM) residues such as Q493R, N501Y, Q498, T478K, and Y505H all contribute significantly to binding affinity with human ACE2. Interactions with Omicron variant mutations at residues 493, 496, 498, and 501 seem to restore ACE2 binding effectiveness lost due to other mutations like K417N.

Pseudoscience and fraudulent products for COVID-19 management

As of now, the COVID-19 pandemic has become uncontrolled and is spreading widely throughout the world. Additionally, new variants of the mutated viral variants have been found in some countries that are more dangerous than the original strain. Even vaccines cannot produce complete protective immunity against the newer strains of SARS-CoV-2. Due to such a dreadful situation, lots of fear and depression have been created among the public. People are looking for the treatment of the disease at any cost and there is a race in the market to provide treatment and make money, whether it is effective or not! In such a condition, many fraud products, remedies, and myths have come into the market, which is falsely claimed to be effective for the disease and can harm the patients. Hence, FDA has banned such products and remedies. In this review, we have compiled all such fraudulent and pseudosciences identified for COVID-19. Currently, in the pandemic time, health agencies are approving the repurposed medicines based on the small-scale clinical data for emergency uses that become ineffective (most of the cases) after large randomized clinical studies. Proper vigilance strategies need to be defined by the regulatory agencies of the nation and routine awareness programs shall be arranged for educating the people and healthcare workers on routine updates.

Heterologous immunity induced by 1st generation COVID-19 vaccines and its role in developing a pan-coronavirus vaccine

Severe acute respiratory syndrome virus-2 (SARS-CoV-2), the causative infectious agent of the COVID-19 pandemic, has led to multiple (4-6) waves of infections worldwide during the past two years. The development of vaccines against SARS-CoV-2 has led to successful mass immunizations worldwide, mitigating the worldwide mortality due the pandemic to a great extent. Yet the evolution of new variants highlights a need to develop a universal vaccine which can prevent infections from all virulent SARS-CoV-2. Most of the current first generation COVID-19 vaccines are based on the Spike protein from the original Wuhan-hu-1 virus strain. It is encouraging that they still protect from serious illnesses, hospitalizations and mortality against a number of mutated viral strains, to varying degrees. Understanding the mechanisms by which these vaccines provide heterologous protection against multiple highly mutated variants can reveal strategies to develop a universal vaccine. In addition, many unexposed individuals have been found to harbor T cells that are cross-reactive against SARS-CoV-2 antigens, with a possible protective role. In this review, we will discuss various aspects of natural or vaccine-induced heterologous (cross-reactive) adaptive immunity against SARS-CoV-2 and other coronaviruses, and their role in achieving the concept of a pan-coronavirus vaccine.

Molecular Dynamics Studies on the Structural Stability Prediction of SARS-CoV-2 Variants Including Multiple Mutants

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has caused the Coronavirus Disease (COVID-19) pandemic worldwide. The spike protein in SARS-CoV-2 fuses with and invades cells in the host respiratory system by binding to angiotensin-converting enzyme 2 (ACE2). The spike protein, however, undergoes continuous mutation from a D614G single mutant to an omicron variant, including multiple mutants. In this study, variants, including multiple mutants (double, triple mutants, B.1.620, delta, alpha, delta_E484Q, mu, and omicron) were investigated in patients. The 3D structure of the full-length spike protein was used in conformational analysis depending on the SARS-CoV-2 variants. The structural stability of the variant types was analyzed based on the distance between the receptor-binding domain (RBD) of each chain in the spike protein and the binding free energy between the spike protein and bound ACE2 in the one-, two-, and three-open-complex forms using molecular dynamics (MD) simulation. Omicron variants, the most prevalent in the recent history of the global pandemic, which consist of 32 mutations, showed higher stability in all open-complex forms compared with that of the wild type and other variants. We suggest that the conformational stability of the spike protein is the one of the important determinants for the differences in viral infectivity among variants, including multiple mutants.

Occurrence of a novel cleavage site for cathepsin G adjacent to the polybasic sequence within the proteolytically sensitive activation loop of the SARS-CoV-2 Omicron variant: The amino acid substitution N679K and P681H of the spike protein

The serine proteases neutrophil elastase (NE), proteinase 3 (PR3), cathepsin G (CatG), and neutrophil serine protease 4 (NSP4) are secreted by activated neutrophils as a part of the innate immune response against invading pathogens. However, these serine proteases might be adopted by viruses to mediate viral surface protein priming resulting in host cell entrance and productive infection. Indeed, NE and PR3 hydrolyze the scissile peptide bond within the proteolytically sensitive polybasic sequence of the activation loop of SARS-CoV-2 located at the S1/S2 interface of the Spike (S) protein; an amino acid motif which differs from SARS-CoV-1. The occurrence of novel SARS-CoV-2 variants and substitution of distinct amino acids at the polybasic sequence prompts serious concerns regarding increased transmissibility. We propose that a novel cleavage site by CatG of the Omicron variant and the increased substrate turnover of the Delta variant by furin within the polybasic sequence should be considered for increased transmission of SARS-CoV-2 variants.

Molecular dynamics simulations of the delta and omicron SARS-CoV-2 spike – ACE2 complexes reveal distinct changes between both variants

SARS-CoV-2, the virus causing the COVID-19 pandemic, changes frequently through the appearance of mutations constantly leading to new variants. However, only few variants evolve as dominating and will be considered as “Variants of Concern” (VOCs) by the world health organization (WHO). At the end of 2020 the alpha (B.1.1.7) variant appeared in the United Kingdom and dominated the pandemic situation until mid of 2021 when it was substituted by the delta variant (B.1.617.2) that first appeared in India as predominant. At the end of 2021, SARS-CoV-2 omicron (B.1.1.529) evolved as the dominating variant. Here, we use in silico modeling and molecular dynamics (MD) simulations of the receptor-binding domain of the viral spike protein and the host cell surface receptor ACE2 to analyze and compare the interaction pattern between the wild type, delta and omicron variants. We identified residue 493 in delta (glutamine) and omicron (arginine) with altered binding properties towards ACE2.

A Bacterial Cell-Based Assay To Study SARS-CoV-2 Protein-Protein Interactions

Methods for detecting and dissecting the interactions of virally encoded proteins are essential for probing basic viral biology and providing a foundation for therapeutic advances. The dearth of targeted therapeutics for the treatment of coronavirus disease 2019 (COVID-19), an ongoing global health crisis, underscores the importance of gaining a deeper understanding of the interactions of proteins encoded by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Here, we describe the use of a convenient bacterial cell-based two-hybrid (B2H) system to analyze the SARS-CoV-2 proteome. We identified 16 distinct intraviral protein-protein interactions (PPIs), involving 16 proteins. We found that many of the identified proteins interact with more than one partner. Further, our system facilitates the genetic dissection of these interactions, enabling the identification of selectively disruptive mutations. We also describe a modified B2H system that permits the detection of disulfide bond-dependent PPIs in the normally reducing Escherichia coli cytoplasm, and we used this system to detect the interaction of the SARS-CoV-2 spike protein receptor-binding domain (RBD) with its cognate cell surface receptor ACE2. We then examined how the RBD-ACE2 interaction is perturbed by several RBD amino acid substitutions found in currently circulating SARS-CoV-2 variants. Our findings illustrate the utility of a genetically tractable bacterial system for probing the interactions of viral proteins and investigating the effects of emerging mutations. In principle, the system could also facilitate the identification of potential therapeutics that disrupt specific interactions of virally encoded proteins. More generally, our findings establish the feasibility of using a B2H system to detect and dissect disulfide bond-dependent interactions of eukaryotic proteins.

Two-Component Nanoparticle Vaccine Displaying Glycosylated Spike S1 Domain Induces Neutralizing Antibody Response against SARS-CoV-2 Variants

Vaccines pave the way out of the SARS-CoV-2 pandemic. Besides mRNA and adenoviral vector vaccines, effective protein-based vaccines are needed for immunization against current and emerging variants. We have developed a virus-like particle (VLP)-based vaccine using the baculovirus-insect cell expression system, a robust production platform known for its scalability, low cost, and safety. Baculoviruses were constructed encoding SARS-CoV-2 spike proteins: full-length S, stabilized secreted S, or the S1 domain. Since subunit S only partially protected mice from SARS-CoV-2 challenge, we produced S1 for conjugation to bacteriophage AP205 VLP nanoparticles using tag/catcher technology. The S1 yield in an insect-cell bioreactor was ∼11 mg/liter, and authentic protein folding, efficient glycosylation, partial trimerization, and ACE2 receptor binding was confirmed. Prime-boost immunization of mice with 0.5 μg S1-VLPs showed potent neutralizing antibody responses against Wuhan and UK/B.1.1.7 SARS-CoV-2 variants. This two-component nanoparticle vaccine can now be further developed to help alleviate the burden of COVID-19.

Genomic Variations in the Structural Proteins of SARS-CoV-2 and Their Deleterious Impact on Pathogenesis: A Comparative Genomics Approach

A continual rise in severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection causing coronavirus disease (COVID-19) has become a global threat. The main problem comes when SARS-CoV-2 gets mutated with the rising infection and becomes more lethal for humankind than ever. Mutations in the structural proteins of SARS-CoV-2, i.e., the spike surface glycoprotein (S), envelope (E), membrane (M) and nucleocapsid (N), and replication machinery enzymes, i.e., main protease (Mpro) and RNA-dependent RNA polymerase (RdRp) creating more complexities towards pathogenesis and the available COVID-19 therapeutic strategies. This study analyzes how a minimal variation in these enzymes, especially in S protein at the genomic/proteomic level, affects pathogenesis. The structural variations are discussed in light of the failure of small molecule development in COVID-19 therapeutic strategies. We have performed in-depth sequence- and structure-based analyses of these proteins to get deeper insights into the mechanism of pathogenesis, structure-function relationships, and development of modern therapeutic approaches. Structural and functional consequences of the selected mutations on these proteins and their association with SARS-CoV-2 virulency and human health are discussed in detail in the light of our comparative genomics analysis.

The SARS-CoV-2 spike protein is vulnerable to moderate electric fields

Most of the ongoing projects aimed at the development of specific therapies and vaccines against COVID-19 use the SARS-CoV-2 spike (S) protein as the main target. The binding of the spike protein with the ACE2 receptor (ACE2) of the host cell constitutes the first and key step for virus entry. During this process, the receptor binding domain (RBD) of the S protein plays an essential role, since it contains the receptor binding motif (RBM), responsible for the docking to the receptor. So far, mostly biochemical methods are being tested in order to prevent binding of the virus to ACE2. Here we show, with the help of atomistic simulations, that external electric fields of easily achievable and moderate strengths can dramatically destabilise the S protein, inducing long-lasting structural damage. One striking field-induced conformational change occurs at the level of the recognition loop L3 of the RBD where two parallel beta sheets, believed to be responsible for a high affinity to ACE2, undergo a change into an unstructured coil, which exhibits almost no binding possibilities to the ACE2 receptor. We also show that these severe structural changes upon electric-field application also occur in the mutant RBDs corresponding to the variants of concern (VOC) B.1.1.7 (UK), B.1.351 (South Africa) and P.1 (Brazil). Remarkably, while the structural flexibility of S allows the virus to improve its probability of entering the cell, it is also the origin of the surprising vulnerability of S upon application of electric fields of strengths at least two orders of magnitude smaller than those required for damaging most proteins. Our findings suggest the existence of a clean physical method to weaken the SARS-CoV-2 virus without further biochemical processing. Moreover, the effect could be used for infection prevention purposes and also to develop technologies for in-vitro structural manipulation of S. Since the method is largely unspecific, it can be suitable for application to other mutations in S, to other proteins of SARS-CoV-2 and in general to membrane proteins of other virus types.

The SARS-CoV-2 Spike protein is essential for viral infectivity and binds to the host receptor ACE2. Here, the authors present MD simulations of the Spike protein and its variants of concern and observe that the Spike protein is destabilised by moderate static electric fields, and undergoes field-induced conformational changes that hinder binding to ACE2.

Computational decomposition reveals reshaping of the SARS-CoV-2-ACE2 interface among viral variants expressing the N501Y mutation

Variants of concern of the SARS-CoV-2 virus with an asparagine-to-tyrosine substitution at position 501 (N501Y) in the receptor-binding domain (RBD) show enhanced infectivity compared to wild-type, resulting in an altered pandemic situation in affected areas. These SARS-Cov-2 variants comprise the two Alpha variants (B.1.1.7, United Kingdom and B.1.1.7 with the additional E484K mutation), the Beta variant (B.1.351, South Africa), and the Gamma variant (P.1, Brazil). Understanding the binding modalities between these viral variants and the host cell receptor ACE2 allows to depict changes, but also common motifs of virus-host cell interaction. The trimeric spike protein expressed at the viral surface contains the RBD that forms the molecular interface with ACE2. All the above-mentioned variants carry between one and three amino acid exchanges within the interface-forming region of the RBD, thereby altering the binding interface with ACE2. Using molecular dynamics (MD) simulations and decomposition of intermolecular contacts between the RBD and ACE2, we identified phenylalanine 486, glutamine 498, threonine 500, and tyrosine 505 as important interface-forming residues across viral variants. However, especially the N501Y exchange increased contact formation for this residue and also induced some local conformational changes. Comparing here, the in silico generated B.1.1.7 RBD-ACE2 complex with the now available experimentally solved structure reveals very similar behavior during MD simulation. We demonstrate, how computational methods can help to identify differences in conformation as well as contact formation for newly emerging viral variants. Altogether, we provide extensive data on all N501Y expressing SARS-CoV-2 variants of concern with respect to their interaction with ACE2 and how this induces reshaping of the RBD-ACE2 interface.

Structural Dissection of Viral Spike-Protein Binding of SARS-CoV-2 and SARS-CoV-1 to the Human Angiotensin-Converting Enzyme 2 (ACE2) as Cellular Receptor

An outbreak by a new severe acute respiratory syndrome betacoronavirus (SARS-CoV-2) has spread CoronaVirus Disease 2019 (COVID-19) all over the world. Immediately, following studies have confirmed the human Angiotensin-Converting Enzyme 2 (ACE2) as a cellular receptor of viral Spike-Protein (Sp) that mediates the CoV-2 invasion into the pulmonary host cells. Here, we compared the molecular interactions of the viral Sp from previous SARS-CoV-1 of 2002 and SARS-CoV-2 with the host ACE2 protein by in silico analysis of the available experimental structures of Sp-ACE2 complexes. The K417 amino acid residue, located in the region of Sp Receptor-Binding Domain (RBD) of the new coronavirus SARS-CoV-2, showed to have a key role for the binding to the ACE2 N-terminal region. The R426 residue of SARS-CoV-1 Sp-RBD also plays a key role, although by interacting with the central region of the ACE2 sequence. Therefore, our study evidenced peculiarities in the interactions of the two Sp-ACE2 complexes. Our outcomes were consistent with previously reported mutagenesis studies on SARS-CoV-1 and support the idea that a new and different RBD was acquired by SARS-CoV-2. These results have interesting implications and suggest further investigations.

Allosteric Cross-Talk among Spike's Receptor-Binding Domain Mutations of the SARS-CoV-2 South African Variant Triggers an Effective Hijacking of Human Cell Receptor

The rapid and relentless emergence of novel highly transmissible SARS-CoV-2 variants, possibly decreasing vaccine efficacy, currently represents a formidable medical and societal challenge. These variants frequently hold mutations on the Spike protein's receptor-binding domain (RBD), which, binding to the angiotensin-converting enzyme 2 (ACE2) receptor, mediates viral entry into host cells. Here, all-atom molecular dynamics simulations and dynamical network theory of the wild-type and mutant RBD/ACE2 adducts disclose that while the N501Y mutation (UK variant) enhances the Spike's binding affinity toward ACE2, the concomitant N501Y, E484K, and K417N mutations (South African variant) aptly adapt to increase SARS-CoV-2 propagation via a two-pronged strategy: (i) effectively grasping ACE2 through an allosteric signaling between pivotal RBD structural elements and (ii) impairing the binding of antibodies elicited by infected or vaccinated patients. This information unlocks the molecular terms and evolutionary strategies underlying the increased virulence of emerging SARS-CoV-2 variants, setting the basis for developing the next-generation anti-COVID-19 therapeutics.