Receptor binding, immune escape, and protein stability direct the natural selection of SARS-CoV-2 variants

Skin patch delivery of a SARS-CoV-2 spike DNA vaccine produces broad neutralising antibody responses

The ongoing SARS-CoV-2 pandemic continues to be a major health burden globally, especially in resource-limited areas. Continued research into more effective and accessible vaccines is required to reduce the burden of disease. Here, we use an emerging vaccine delivery system, the high-density microarray patch (HD-MAP) to deliver a plasmid DNA vaccine (Delta 6P) encoding for the SARS-CoV-2 spike protein. HD-MAP delivery of this vaccine resulted in robust IgG responses in mice against multiple domains of the spike protein. The cellular response to vaccination was also measured, and comparative analysis showed that relative to intramuscular vaccination, HD-MAP vaccination elicited spike-specific CD4+ T and CD8+ T cell responses that were largely comparable, but the number of polyfunctional CD4+ T cells was higher in the HD-MAP group. Collectively, this work suggests that HD-MAP delivery of the Delta 6P vaccine is effective against SARS-CoV-2, warranting further investigation.

A comprehensive review of current insights into the virulence factors of SARS-CoV-2

The evolution of SARS-CoV-2 pathogenicity has been a major focus of attention. However, the determinants of pathogenicity are still unclear. Various hypotheses have attempted to elucidate the mechanisms underlying the evolution of viral pathogenicity, but a definitive conclusion has yet to be reached. Here, we review the potential impact of all proteins in SARS-CoV-2 on the viral pathogenic process and analyze the effects of their mutations on pathogenicity evolution. We aim to summarize which virus-encoded proteins are crucial in influencing viral pathogenicity, defined as disease severity following infection. Mutations in these key proteins, which are the virulence factors in SARS-CoV-2, may be the driving forces behind the evolution of viral pathogenicity. Mutations in the S protein can impact viral entry and fusogenicity. Mutations in proteins such as NSP2, NSP5, NSP14, and ORF7a can alter the virus's ability to suppress host protein synthesis and innate immunity. Mutations in NSP3, NSP4, NSP6, N protein, NSP5, and NSP12 may alter viral replication efficiency. The combined effects of mutations in the S protein and NSP6 can significantly reduce viral replication. In addition, various viral proteins, including ORF3a, ORF8, NSP4, Spike protein, N protein, and E protein, directly participate in the inflammatory process. Mutations in these proteins can modulate the levels of inflammation following infection. Collectively, these viral protein mutations can influence SARS-CoV-2 pathogenicity by impacting viral immune evasion, replication capacity, and the level of inflammation mediated by infection. In conclusion, the evolution of SARS-CoV-2 pathogenicity is likely determined by multiple virulence factors.

Longitudinal analysis of SARS-CoV-2 IgG antibody durability in Puerto Rico

Understanding the dynamics of antibody responses following vaccination and SARS-CoV-2 infection is important for informing effective vaccination strategies and other public health interventions. This study investigates SARS-CoV-2 antibody dynamics in a Puerto Rican cohort, analyzing how IgG levels vary by vaccination status and previous infection. We assess waning immunity and the distribution of hybrid immunity with the aim to inform public health strategies and vaccination programs in Puerto Rico and similar settings. We conducted a prospective, longitudinal cohort study to identify SARS-CoV-2 infections and related outcomes in Ponce, Puerto Rico, from June 2020-August 2022. Participants provided self-collected nasal swabs every week and serum every six months for RT-PCR and IgG testing, respectively. IgG reactivity against nucleocapsid (N) antigens, which generally indicate previous infection, and spike (S1) and receptor-binding domain (RBD) antigens, which indicate history of either infection or vaccination, was assessed using the Luminex Corporation xMAP® SARS-CoV-2 Multi-Antigen IgG Assay. Prior infection was defined by positive RT-PCRs, categorized by the predominant circulating SARS-CoV-2 variant at the event time. Demographic information, medical history, and COVID-19 vaccination history were collected through standardized questionnaires. Of 882 participants included in our analysis, 34.0% experienced at least one SARS-CoV-2 infection, with most (78.7%) occurring during the Omicron wave (December 2021 onwards). SARS-CoV-2 antibody prevalence increased over time, reaching 98.4% by the final serum collection, 67.0% attributable to vaccination alone, 1.6% from infection alone, and 31.4% from both. Regardless of prior infection status, RBD and S1 IgG levels gradually declined following two vaccine doses. A third dose boosted these antibody levels and showed a slower decline over time. N-antibody levels peaked during the Omicron surge and waned over time. Vaccination in individuals with prior SARS-CoV-2 infection elicited the highest and most durable antibody responses. N or S1 seropositivity was associated with lower odds of a subsequent positive PCR test during the Omicron period, with N antibodies showing a stronger association. By elucidating the differential decay of RBD and S1 antibodies following vaccination and the complexities of N-antibody response following infection, this study in a Puerto Rican cohort strengthens the foundation for developing targeted interventions and public health strategies.

Rheostatic contributions to protein stability can obscure a position's functional role

Rheostat positions, which can be substituted with various amino acids to tune protein function across a range of outcomes, are a developing area for advancing personalized medicine and bioengineering. Current methods cannot accurately predict which proteins contain rheostat positions or their substitution outcomes. To compare the prevalence of rheostat positions in homologs, we previously investigated their occurrence in two pyruvate kinase (PYK) isozymes. Human liver PYK contained numerous rheostat positions that tuned the apparent affinity for the substrate phosphoenolpyruvate (Kapp-PEP) across a wide range. In contrast, no functional rheostat positions were identified in Zymomonas mobilis PYK (ZmPYK). Further, the set of ZmPYK substitutions included an unusually large number that lacked measurable activity. We hypothesized that the inactive substitution variants had reduced protein stability, precluding detection of Kapp-PEP tuning. Using modified buffers, robust enzymatic activity was obtained for 19 previously-inactive ZmPYK substitution variants at three positions. Surprisingly, both previously-inactive and previously-active substitution variants all had Kapp-PEP values close to wild-type. Thus, none of the three positions were functional rheostat positions, and, unlike human liver PYK, ZmPYK's Kapp-PEP remained poorly tunable by single substitutions. To directly assess effects on stability, we performed thermal denaturation experiments for all ZmPYK substitution variants. Many diminished stability, two enhanced stability, and the three positions showed different thermal sensitivity to substitution, with one position acting as a "stability rheostat." The differences between the two PYK homologs raises interesting questions about the underlying mechanism(s) that permit functional tuning by single substitutions in some proteins but not in others.

Biophysical principles predict fitness of SARS-CoV-2 variants

SARS-CoV-2 employs its spike protein's receptor binding domain (RBD) to enter host cells. The RBD is constantly subjected to immune responses, while requiring efficient binding to host cell receptors for successful infection. However, our understanding of how RBD's biophysical properties contribute to SARS-CoV-2's epidemiological fitness remains largely incomplete. Through a comprehensive approach, comprising large-scale sequence analysis of SARS-CoV-2 variants and the identification of a fitness function based on binding thermodynamics, we unravel the relationship between the biophysical properties of RBD variants and their contribution to viral fitness. We developed a biophysical model that uses statistical mechanics to map the molecular phenotype space, characterized by dissociation constants of RBD to ACE2, LY-CoV016, LY-CoV555, REGN10987, and S309, onto an epistatic fitness landscape. We validate our findings through experimentally measured and machine learning (ML) estimated binding affinities, coupled with infectivity data derived from population-level sequencing. Our analysis reveals that this model effectively predicts the fitness of novel RBD variants and can account for the epistatic interactions among mutations, including explaining the later reversal of Q493R. Our study sheds light on the impact of specific mutations on viral fitness and delivers a tool for predicting the future epidemiological trajectory of previously unseen or emerging low-frequency variants. These insights offer not only greater understanding of viral evolution but also potentially aid in guiding public health decisions in the battle against COVID-19 and future pandemics.

Biophysical principles predict fitness of SARS-CoV-2 variants

SARS-CoV-2 employs its spike protein's receptor binding domain (RBD) to enter host cells. The RBD is constantly subjected to immune responses, while requiring efficient binding to host cell receptors for successful infection. However, our understanding of how RBD's biophysical properties contribute to SARS-CoV-2's epidemiological fitness remains largely incomplete. Through a comprehensive approach, comprising large-scale sequence analysis of SARS-CoV-2 variants and the discovery of a fitness function based on binding thermodynamics, we unravel the relationship between the biophysical properties of RBD variants and their contribution to viral fitness. We developed a biophysical model that uses statistical mechanics to map the molecular phenotype space, characterized by binding constants of RBD to ACE2, LY-CoV016, LY-CoV555, REGN10987, and S309, onto a epistatic fitness landscape. We validate our findings through experimentally measured and machine learning (ML) estimated binding affinities, coupled with infectivity data derived from population-level sequencing. Our analysis reveals that this model effectively predicts the fitness of novel RBD variants and can account for the epistatic interactions among mutations, including explaining the later reversal of Q493R. Our study sheds light on the impact of specific mutations on viral fitness and delivers a tool for predicting the future epidemiological trajectory of previously unseen or emerging low frequency variants. These insights offer not only greater understanding of viral evolution but also potentially aid in guiding public health decisions in the battle against COVID-19 and future pandemics.

Molecular Dynamics Simulations on Spike Protein Mutants Binding with Human β Defensin Type 2

Human β defensin type 2 (hBD-2), a cationic cysteine-rich peptide secreted from the human innate immune system, can bind Spike-RBD at the same site as receptor ACE2, thus blocking viral entry into ACE2-expressing cells. In order to find out the impact of CoV-2 mutations on hBD-2's antiviral activity, it is important to investigate the binding and interaction of hBD-2 with RBD mutants. All-atom molecular dynamics simulations were conducted on typical RBD mutants, including N501Y, E484K, P479S, T478I, S477N, N439K, K417N, and N501Y-E484K-K417N, binding with hBD-2. Starting from the stable binding structure of hBD-2 and wt-RBD and ClusPro and HADDOCK docking-predicted initial structures, the RBD variants bound with hBD-2 simulations were set up, and NAMD simulations were conducted. Based on the structure and dynamics analysis, it was found that most RBD variants can still form a similar number of hydrogen bonds with hBD-2, in addition to having a similar-sized buried surface area (BSA) and a similar binding interface to the RBD wildtype. However, the RBD triple mutant formed a less stable binding structure with hBD-2 than other variants. Additionally, the free energy perturbation (FEP) method was applied to calculate the contribution of key mutant residues to the binding and the free energy change caused by the mutations. The result shows that N439K, K417N, and the trimutation increase the binding free energy of RBD with hBD-2; thus, RBD should bind less stably with hBD-2. E484K decreases the binding free energy, thus it should bind more stably with hBD-2, while N501Y, S477N, T478I, and P479S almost do not change the binding free energy with hBD-2. The MM-GBSA method predicted the binding interaction energy which shows that the trimutant should be able to escape the binding with hBD-2 but N501Y should not. The result can provide insight into understanding the functional mechanism of hBD-2 combating SARS-CoV-2 mutants.

Cryptic-site-specific antibodies to the SARS-CoV-2 receptor binding domain can retain functional binding affinity to spike variants

IMPORTANCE

Multiple SARS-CoV-2 variants of concern have emerged and caused a significant number of infections and deaths worldwide. These variants of concern contain mutations that might significantly affect antigen-targeting by antibodies. It is therefore important to further understand how antibody binding and neutralization are affected by the mutations in SARS-CoV-2 variants. We highlighted how antibody epitope specificity can influence antibody binding to SARS-CoV-2 spike protein variants and neutralization of SARS-CoV-2 variants. We showed that weakened spike binding and neutralization of Beta (B.1.351) and Omicron (BA.1) variants compared to wildtype are not universal among the panel of antibodies and identified antibodies of a specific binding footprint exhibiting consistent enhancement of spike binding and retained neutralization to Beta variant. These data and analysis can inform how antigen-targeting by antibodies might evolve during a pandemic and prepare for potential future sarbecovirus outbreaks.

Biophysical evolution of the receptor-binding domains of SARS-CoVs

With hundreds of coronaviruses (CoVs) identified in bats that can infect humans, it is essential to understand how CoVs that affected the human population have evolved. Seven known CoVs have infected humans, of which three CoVs caused severe disease with high mortalities: severe acute respiratory syndrome (SARS)-CoV emerged in 2002, Middle East respiratory syndrome-CoV in 2012, and SARS-CoV-2 in 2019. SARS-CoV and SARS-CoV-2 belong to the same family, follow the same receptor pathway, and use their receptor-binding domain (RBD) of spike protein to bind to the angiotensin-converting enzyme 2 (ACE2) receptor on the human epithelial cell surface. The sequence of the two RBDs is divergent, especially in the receptor-binding motif that directly interacts with ACE2. We probed the biophysical differences between the two RBDs in terms of their structure, stability, aggregation, and function. Since RBD is being explored as an antigen in protein subunit vaccines against CoVs, determining these biophysical properties will also aid in developing stable protein subunit vaccines. Our results show that, despite RBDs having a similar three-dimensional structure, they differ in their thermodynamic stability. RBD of SARS-CoV-2 is significantly less stable than that of SARS-CoV. Correspondingly, SARS-CoV-2 RBD shows a higher aggregation propensity. Regarding binding to ACE2, less stable SARS-CoV-2 RBD binds with a higher affinity than more stable SARS-CoV RBD. In addition, SARS-CoV-2 RBD is more homogenous in terms of its binding stoichiometry toward ACE2 compared to SARS-CoV RBD. These results indicate that SARS-CoV-2 RBD differs from SARS-CoV RBD in terms of its stability, aggregation, and function, possibly originating from the diverse receptor-binding motifs. Higher aggregation propensity and decreased stability of SARS-CoV-2 RBD warrant further optimization of protein subunit vaccines that use RBD as an antigen by inserting stabilizing mutations or formulation screening.

Structure, dynamics and free energy studies on the effect of point mutations on SARS-CoV-2 spike protein binding with ACE2 receptor

The ongoing COVID-19 pandemic continues to infect people worldwide, and the virus continues to evolve in significant ways which can pose challenges to the efficiency of available vaccines and therapeutic drugs and cause future pandemic. Therefore, it is important to investigate the binding and interaction of ACE2 with different RBD variants. A comparative study using all-atom MD simulations was conducted on ACE2 binding with 8 different RBD variants, including N501Y, E484K, P479S, T478I, S477N, N439K, K417N and N501Y-E484K-K417N on RBD. Based on the RMSD, RMSF, and DSSP results, overall the binding of RBD variants with ACE2 is stable, and the secondary structure of RBD and ACE2 are consistent after the point mutation. Besides that, a similar buried surface area, a consistent binding interface and a similar amount of hydrogen bonds formed between RBD and ACE2 although the exact residue pairs on the binding interface were modified. The change of binding free energy from point mutation was predicted using the free energy perturbation (FEP) method. It is found that N501Y, N439K, and K417N can strengthen the binding of RBD with ACE2, while E484K and P479S weaken the binding, and S477N and T478I have negligible effect on the binding. Point mutations modified the dynamic correlation of residues in RBD based on the dihedral angle covariance matrix calculation. Doing dynamic network analysis, a common intrinsic network community extending from the tail of RBD to central, then to the binding interface region was found, which could communicate the dynamics in the binding interface region to the tail thus to the other sections of S protein. The result can supply unique methodology and molecular insight on studying the molecular structure and dynamics of possible future pandemics and design novel drugs.

Multiplex Assays Enable Simultaneous Detection and Identification of SARS-CoV-2 Variants of Concern in Clinical and Wastewater Samples

The targeted screening and sequencing approaches for COVID-19 surveillance need to be adjusted to fit the evolving surveillance objectives which necessarily change over time. We present the development of variant screening assays that can be applied to new targets in a timely manner and enable multiplexing of targets for efficient implementation in the laboratory. By targeting the HV69/70 deletion for Alpha, K417N for Beta, K417T for Gamma, and HV69/70 deletion plus K417N for sub-variants BA.1, BA.3, BA.4, and BA.5 of Omicron, we achieved simultaneous detection and differentiation of Alpha, Beta, Gamma, and Omicron in a single assay. Targeting both T478K and P681R mutations enabled specific detection of the Delta variant. The multiplex assays used in combination, targeting K417N and T478K, specifically detected the Omicron sub-variant BA.2. The limits of detection for the five variants of concern were 4-16 copies of the viral RNA per reaction. Both assays achieved 100% clinical sensitivity and 100% specificity. Analyses of 377 clinical samples and 24 wastewater samples revealed the Delta variant in 100 clinical samples (nasopharyngeal and throat swab) collected in November 2021. Omicron BA.1 was detected in 79 nasopharyngeal swab samples collected in January 2022. Alpha, Beta, and Gamma variants were detected in 24 wastewater samples collected in May-June 2021 from two major cities of Alberta (Canada), and the results were consistent with the clinical cases of multiple variants reported in the community.

Viral Capsid, Antibody, and Receptor Interactions: Experimental Analysis of the Antibody Escape Evolution of Canine Parvovirus

Canine parvovirus (CPV) is a small nonenveloped single-stranded DNA virus that causes serious diseases in dogs worldwide. The original strain of the virus (CPV-2) emerged in dogs during the late 1970s due to a host range switch of a virus similar to the feline panleukopenia virus that infected another host. The virus that emerged in dogs had altered capsid receptor and antibody binding sites, with some changes affecting both functions. Further receptor and antibody binding changes arose when the virus became better adapted to dogs or to other hosts. Here, we used in vitro selection and deep sequencing to reveal how two antibodies with known interactions select for escape mutations in CPV. The antibodies bound two distinct epitopes, and one largely overlapped the host receptor binding site. We also generated mutated antibody variants with altered binding structures. Viruses were passaged with wild-type (WT) or mutated antibodies, and their genomes were deep sequenced during the selective process. A small number of mutations were detected only within the capsid protein gene during the first few passages of selection, and most sites remained polymorphic or were slow to go to fixation. Mutations arose both within and outside the antibody binding footprints on the capsids, and all avoided the transferrin receptor type 1 binding footprint. Many selected mutations matched those that have arisen in the natural evolution of the virus. The patterns observed reveal the mechanisms by which these variants have been selected in nature and provide a better understanding of the interactions between antibody and receptor selections. IMPORTANCE Antibodies protect animals against infection by many different viruses and other pathogens, and we are gaining new information about the epitopes that induce antibody responses against viruses and the structures of the bound antibodies. However, less is known about the processes of antibody selection and antigenic escape and the constraints that apply in this system. Here, we used an in vitro model system and deep genome sequencing to reveal the mutations that arose in the virus genome during selection by each of two monoclonal antibodies or their mutated variants. High-resolution structures of each of the Fab:capsid complexes revealed their binding interactions. The wild-type antibodies or their mutated variants allowed us to examine how changes in antibody structure influence the mutational selection patterns seen in the virus. The results shed light on the processes of antibody binding, neutralization escape, and receptor binding, and they likely have parallels for many other viruses.

Adsorption of SARS-CoV-2 Spike (N501Y) RBD to Human Angiotensin-Converting Enzyme 2 at a Lipid/Water Interface

The receptor binding domain (RBD) of spike proteins plays a crucial role in the process of severe acute respiratory syndrome corona virus 2 (SARS-CoV-2) attachment to the human angiotensin-converting enzyme 2 (ACE2). The N501Y mutation and later mutations introduced extra positive charges on the spike RBD and resulted in higher transmissibility, likely due to stronger binding with the highly negatively charged ACE2. Consequently, many studies have been devoted to understanding the molecular mechanism of spike protein binding with the ACE2 receptor. Most of the theoretical studies, however, have been done on isolated proteins. ACE2 is a transmembrane protein; thus, it is important to understand the interaction of spike proteins with ACE2 in a lipid matrix. In this study, the adsorption of ACE2 and spike (N501Y) RBD at a lipid/water interface was studied using the heterodyne-detected vibrational sum frequency generation (HD-VSFG) technique. The technique is a non-linear optical spectroscopy which measures vibrational spectra of molecules at an interface and provides information on their structure and orientation. It is found that ACE2 is effectively adsorbed at the positively charged 1,2-dipalmitoyl-3-trimethylammonium-propane (DPTAP) lipid monolayer via electrostatic interactions. The adsorption of ACE2 at the DPTAP monolayer causes a reorganization of interfacial water (D2O) from the D-down to the D-up orientation, indicating that the originally positively charged DPTAP interface becomes negatively charged due to ACE2 adsorption. The negatively charged interface (DPTAP/ACE2) allows further adsorption of positively charged spike RBD. HD-VSFG spectra in the amide I region show differences for spike (N501Y) RBD adsorbed at D2O, DPTAP, and DPTAP/ACE2 interfaces. A red shift observed for the spectra of spike RBD/DPTAP suggests that spike RBD oligomers are formed upon contact with DPTAP lipids.

Viral capsid, antibody, and receptor interactions: experimental analysis of the antibody escape evolution of canine parvovirus

Canine parvovirus (CPV) is a small non-enveloped single-stranded DNA virus that causes serious diseases in dogs worldwide. The original strain of the virus (CPV-2) emerged in dogs during the late-1970s due to a host range switch of a virus similar to the feline panleukopenia virus (FPV) that infected another host. The virus that emerged in dogs had altered capsid receptor- and antibody-binding sites, with some changes affecting both functions. Further receptor and antibody binding changes arose when the virus became better adapted to dogs or to other hosts. Here, we use in vitro selection and deep sequencing to reveal how two antibodies with known interactions select for escape mutations in CPV. The antibodies bind two distinct epitopes, and one largely overlaps the host receptor binding site. We also engineered antibody variants with altered binding structures. Viruses were passaged with the wild type or mutated antibodies, and their genomes deep sequenced during the selective process. A small number of mutations were detected only within the capsid protein gene during the first few passages of selection, and most sites remained polymorphic or were slow to go to fixation. Mutations arose both within and outside the antibody binding footprints on the capsids, and all avoided the TfR-binding footprint. Many selected mutations matched those that have arisen in the natural evolution of the virus. The patterns observed reveal the mechanisms by which these variants have been selected in nature and provide a better understanding of the interactions between antibody and receptor selections.

IMPORTANCE

Antibodies protect animals against infection by many different viruses and other pathogens, and we are gaining new information about the epitopes that induce antibody responses against viruses and the structures of the bound antibodies. However, less is known about the processes of antibody selection and antigenic escape and the constraints that apply in this system. Here, we use an in vitro model system and deep genome sequencing to reveal the mutations that arise in the virus genome during selection by each of two monoclonal antibodies or their engineered variants. High-resolution structures of each of the Fab: capsid complexes revealed their binding interactions. The engineered forms of the wild-type antibodies or mutant forms allowed us to examine how changes in antibody structure influence the mutational selection patterns seen in the virus. The results shed light on the processes of antibody binding, neutralization escape, and receptor binding, and likely have parallels for many other viruses.

An in-silico pharmacophore-based molecular docking study to evaluate the inhibitory potentials of novel fungal triterpenoid Astrakurkurone analogues against a hypothetical mutated main protease of SARS-CoV-2 virus

BACKGROUND

The main protease is an important structural protein of SARS-CoV-2, essential for its survivability inside a human host. Considering current vaccines' limitations and the absence of approved therapeutic targets, Mpro may be regarded as the potential candidate drug target. Novel fungal phytocompound Astrakurkurone may be studied as the potential Mpro inhibitor, considering its medicinal properties reported elsewhere.

METHODS

In silico molecular docking was performed with Astrakurkurone and its twenty pharmacophore-based analogues against the native Mpro protein. A hypothetical Mpro was also constructed with seven mutations and targeted by Astrakurkurone and its analogues. Furthermore, multiple parameters such as statistical analysis (Principal Component Analysis), pharmacophore alignment, and drug likeness evaluation were performed to understand the mechanism of protein-ligand molecular interaction. Finally, molecular dynamic simulation was done for the top-ranking ligands to validate the result.

RESULT

We identified twenty Astrakurkurone analogues through pharmacophore screening methodology. Among these twenty compounds, two analogues namely, ZINC89341287 and ZINC12128321 showed the highest inhibitory potentials against native and our hypothetical mutant Mpro, respectively (-7.7 and -7.3 kcal mol-1) when compared with the control drug Telaprevir (-5.9 and -6.0 kcal mol-1). Finally, we observed that functional groups of ligands namely two aromatic and one acceptor groups were responsible for the residual interaction with the target proteins. The molecular dynamic simulation further revealed that these compounds could make a stable complex with their respective protein targets in the near-native physiological condition.

CONCLUSION

To conclude, Astrakurkurone analogues ZINC89341287 and ZINC12128321 can be potential therapeutic agents against the highly infectious SARS-CoV-2 virus.

Exhaustive mutational analysis of severe acute respiratory syndrome coronavirus 2 ORF3a: An essential component in the pathogen's infectivity cycle

Detailed knowledge of a protein's key residues may assist in understanding its function and designing inhibitors against it. Consequently, such knowledge of one of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)'s proteins is advantageous since the virus is the etiological agent behind one of the biggest health crises of recent times. To that end, we constructed an exhaustive library of bacteria differing from each other by the mutated version of the virus's ORF3a viroporin they harbor. Since the protein is harmful to bacterial growth due to its channel activity, genetic selection followed by deep sequencing could readily identify mutations that abolish the protein's function. Our results have yielded numerous mutations dispersed throughout the sequence that counteract ORF3a's ability to slow bacterial growth. Comparing these data with the conservation pattern of ORF3a within the coronavirinae provided interesting insights: Deleterious mutations obtained in our study corresponded to conserved residues in the protein. However, despite the comprehensive nature of our mutagenesis coverage (108 average mutations per site), we could not reveal all of the protein's conserved residues. Therefore, it is tempting to speculate that our study unearthed positions in the protein pertinent to channel activity, while other conserved residues may correspond to different functionalities of ORF3a. In conclusion, our study provides important information on a key component of SARS-CoV-2 and establishes a procedure to analyze other viroporins comprehensively.

Thermodynamic analysis of the interactions between human ACE2 and spike RBD of Betacoronaviruses (SARS-CoV-1 and SARS-CoV-2)

There are many scientific reports on the interaction of the SARS-CoV-2 virus S protein (and its RBD) with the human ACE2 receptor protein. However, there are no reliable data on how this interaction differs from the interaction of the receptor binding domain of SARS-CoV-1 with ACE2, in terms of binding strength and changes in reaction enthalpy and entropy. Our studies have revealed these differences and the impact of zinc ions on this interaction. Intriguingly, the binding affinity of both RBDs (of SARS-CoV-1 and of SARS-CoV-2) to the ACE2 receptor protein is almost identical; however, there are some differences in the entropic and enthalpic contributions to these interactions.

A comprehensive mini-review on amyloidogenesis of different SARS-CoV-2 proteins and its effect on amyloid formation in various host proteins

Amyloidogenesis is the inherent ability of proteins to change their conformation from native state to cross β-sheet rich fibrillar structures called amyloids which result in a wide range of diseases like Parkinson's disease, Alzheimer's disease, Finnish familial amyloidosis, ATTR amyloidosis, British and Danish dementia, etc. COVID-19, on the other hand is seen to have many similarities in symptoms with other amyloidogenic diseases and the overlap of these morbidities and symptoms led to the proposition whether SARS-CoV-2 proteins are undergoing amyloidogenesis and whether it is resulting in or aggravating amyloidogenesis of any human host protein. Thus the SARS-CoV-2 proteins in infected cells, i.e., Spike (S) protein, Nucleocapsid (N) protein, and Envelope (E) protein were tested via different machinery and amyloidogenesis in them were proven. In this review, we will analyze the pathway of amyloid formation in S-protein, N-protein, E-protein along with the effect that SARS-CoV-2 is creating on various host proteins leading to the unexpected onset of many morbidities like COVID-induced Acute Respiratory Distress Syndrome (ARDS), Parkinsonism in young COVID patients, formation of fibrin microthrombi in heart, etc., and their future implications.

High-Throughput Molecular Dynamics-Based Alchemical Free Energy Calculations for Predicting the Binding Free Energy Change Associated with the Selected Omicron Mutations in the Spike Receptor-Binding Domain of SARS-CoV-2

The ongoing pandemic caused by SARS-CoV-2 has gone through various phases. Since the initial outbreak, the virus has mutated several times, with some lineages showing even stronger infectivity and faster spread than the original virus. Among all the variants, omicron is currently classified as a variant of concern (VOC) by the World Health Organization, as the previously circulating variants have been replaced by it. In this work, we have focused on the mutations observed in omicron sub lineages BA.1, BA.2, BA.4 and BA.5, particularly at the receptor-binding domain (RBD) of the spike protein that is responsible for the interactions with the host ACE2 receptor and binding of antibodies. Studying such mutations is particularly important for understanding the viral infectivity, spread of the disease and for tracking the escape routes of this virus from antibodies. Molecular dynamics (MD) based alchemical free energy calculations have been shown to be very accurate in predicting the free energy change, due to a mutation that could have a deleterious or a stabilizing effect on either the protein itself or its binding affinity to another protein. Here, we investigated the significance of five spike RBD mutations on the stability of the spike protein binding to ACE2 by free energy calculations using high throughput MD simulations. For comparison, we also used conventional MD simulations combined with a Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) based approach, and compared our results with the available experimental data. Overall, the alchemical free energy calculations performed far better than the MM-GBSA approach in predicting the individual impact of the mutations. When considering the experimental variation, the alchemical free energy method was able to produce a relatively accurate prediction for N501Y, the mutant that has previously been reported to increase the binding affinity to hACE2. On the other hand, the other individual mutations seem not to have a significant effect on the spike RBD binding affinity towards hACE2.

SARS-CoV-2 variants of concern and spike protein mutational dynamics in a Swedish cohort during 2021, studied by Nanopore sequencing

BACKGROUND

Since the beginning of the COVID-19 pandemic, new variants of significance to public health have emerged. Consequently, early detection of new mutations and variants through whole-genome sequencing remains crucial to assist health officials in employing appropriate public health measures.

METHODS

We utilized the ARTIC Network SARS-CoV-2 tiled amplicon approach and Nanopore sequencing to sequence 4,674 COVID-19 positive patient samples from Uppsala County, Sweden, between week 15 and 52 in 2021. Using this data, we mapped the circulating variants of concern (VOC) in the county over time and analysed the Spike (S) protein mutational dynamics in the Delta variant throughout 2021.

RESULTS

The distribution of the SARS-CoV-2 VOC matched the national VOC distribution in Sweden, in 2021. In the S protein of the Delta variant, we detected mutations attributable to variants under monitoring and variants of interest (e.g., E484Q, Q613H, Q677H, A222V and Y145H) and future VOC (e.g., T95I and Y144 deletion, which are signature mutations in the Omicron variant). We also frequently detected some less well-described S protein mutations in our Delta sequences, that might play a role in shaping future emerging variants. These include A262S, Q675K, I850L, Q1201H, V1228L and M1237I. Lastly, we observed that some of the Delta variant's signature mutations were underrepresented in our study due to artifacts of the used bioinformatics tools, approach and sequencing method. We therefore discuss some pitfalls and considerations when sequencing SARS-CoV-2 genomes.

CONCLUSION

Our results suggest that genomic surveillance in a small, representative cohort can be used to make predictions about the circulating variants nationally. Moreover, we show that detection of transient mutations in currently circulating variants can give valuable clues to signature mutations of future VOC. Here we suggest six such mutations, that we detected frequently in the Delta variant during 2021. Lastly, we report multiple systematic errors that occurred when following the ARTIC Network SARS-CoV-2 tiled amplicon approach using the V3 primers and Nanopore sequencing, which led to the masking of some of the important signature mutations in the Delta sequences.

Evolutionary progression of collective mutations in Omicron sub-lineages towards efficient RBD-hACE2: Allosteric communications between and within viral and human proteins

The interaction between the Spike (S) protein of SARS-CoV-2 and the human angiotensin converting enzyme 2 (hACE2) is essential for infection, and is a target for neutralizing antibodies. Consequently, selection of mutations in the S protein is expected to be driven by the impact on the interaction with hACE2 and antibody escape. Here, for the first time, we systematically characterized the collective effects of mutations in each of the Omicron sub-lineages (BA.1, BA.2, BA.3 and BA.4) on both the viral S protein receptor binding domain (RBD) and the hACE2 protein using post molecular dynamics studies and dynamic residue network (DRN) analysis. Our analysis suggested that Omicron sub-lineage mutations result in altered physicochemical properties that change conformational flexibility compared to the reference structure, and may contribute to antibody escape. We also observed changes in the hACE2 substrate binding groove in some sub-lineages. Notably, we identified unique allosteric communication paths in the reference protein complex formed by the DRN metrics betweenness centrality and eigencentrality hubs, originating from the RBD core traversing the receptor binding motif of the S protein and the N-terminal domain of the hACE2 to the active site. We showed allosteric changes in residue network paths in both the RBD and hACE2 proteins due to Omicron sub-lineage mutations. Taken together, these data suggest progressive evolution of the Omicron S protein RBD in sub-lineages towards a more efficient interaction with the hACE2 receptor which may account for the increased transmissibility of Omicron variants.

Structural insights into the binding of SARS-CoV-2, SARS-CoV, and hCoV-NL63 spike receptor-binding domain to horse ACE2

Severe acute respiratory syndrome coronavirus (SARS-CoV), SARS-CoV-2, and human coronavirus (hCoV)-NL63 utilize ACE2 as the functional receptor for cell entry, which leads to zoonotic infection. Horses (Equus caballus) attracted our attention because the spike protein receptor-binding domains (RBDs) of SARS-CoV-2 and SARS-CoV-2-related coronaviruses bind equine ACE2 (eACE2) with high affinity. Here we show that eACE2 binds the RBDs of these three coronaviruses and also SARS-CoV-2 variants but with lower affinities compared with human ACE2 (hACE2). Structural analysis and mutation assays indicated that eACE2-H41 accounts for the lower binding affinity of eACE2 to the RBDs of SARS-CoV-2 variants (Alpha, Beta, and Gamma), SARS-CoV, and hCoV-NL63. Pseudovirus infection assays showed that the SARS-CoV-2 Delta strain (B.1.617.2) displayed a significantly increased infection efficiency in eACE2-expressing HeLa cells. Our results reveal the molecular basis of eACE2 binding to the RBDs of SARS-CoV, SARS-CoV-2, and hCoV-NL63, which provides insights into the potential animal transmission of these ACE2-dependent coronaviruses.

Biophysical Fitness Landscape of the SARS-CoV-2 Delta Variant Receptor Binding Domain

Among the five known SARS-CoV-2 variants of concern, Delta is the most virulent leading to severe symptoms and increased mortality among infected people. Our study seeks to examine how the biophysical parameters of the Delta variant correlate to the clinical observations. Receptor binding domain (RBD) is the first point of contact with the human host cells and is the immunodominant form of the spike protein. Delta variant RBD contains two novel mutations L452R and T478K. We examined the effect of single as well as the double mutations on RBD expression in human Expi293 cells, RBD stability using urea and thermal denaturation, and RBD binding to angiotensin converting enzyme 2 (ACE2) receptor and to neutralizing antibodies using isothermal titration calorimetry. Delta variant RBD showed significantly higher expression compared to the wild-type RBD, and the increased expression is due to L452R mutation. Despite their non-conservative nature, none of the mutations significantly affected RBD structure and stability. All mutants showed similar binding affinity to ACE2 and to Class 1 antibodies (CC12.1 and LY-CoV016) as that of the wild-type. Delta double mutant L452R/T478K showed no binding to Class 2 antibodies (P2B-2F6 and LY-CoV555) and a hundred-fold weaker binding to a Class 3 antibody (REGN10987), and the decreased antibody binding is determined by the L452R mutation. These results indicate that the immune escape from neutralizing antibodies, rather than increased receptor binding, is the main biophysical parameter that determined the fitness landscape of the Delta variant RBD.

A Paradigm Shift in the Combination Changes of SARS-CoV-2 Variants and Increased Spread of Delta Variant (B.1.617.2) across the World

Since September 2020, the SARS-CoV-2 variants have gained their dominance worldwide, especially in Kenya, Italy, France, the UK, Turkey, Indonesia, India, Finland, Ireland, Singapore, Denmark, Germany, and Portugal. In this study, we developed a model on the frequency of delta variants across 28 countries (R²= 0.1497), displaying the inheritance of mutations during the generation of the delta variants with 123,526 haplotypes. The country-wise haplotype network showed the distribution of haplotypes in USA (10,174), Denmark (5,637), India (4,089), Germany (2,350), Netherlands (1,899), Sweden (1,791), Italy (1,720), France (1,293), Ireland (1,257), Belgium (1,207), Singapore (1,193), Portugal (1,184) and Spain (1,133). Our analysis shows the highest haplotype in Europe with 84% and the lowest in Australia with 0.00001%. A model of scatter plot was generated with a regression line which provided the estimated rate of mutation, including 24.048 substitutions yearly. Our study concluded that the high global prevalence of the delta variants is due to a high frequency of infectivity, supporting the paradigm shift of the viral variants.

Amyloidogenesis of SARS-CoV-2 Spike Protein

SARS-CoV-2 infection is associated with a surprising number of morbidities. Uncanny similarities with amyloid-disease associated blood coagulation and fibrinolytic disturbances together with neurologic and cardiac problems led us to investigate the amyloidogenicity of the SARS-CoV-2 spike protein (S-protein). Amyloid fibril assays of peptide library mixtures and theoretical predictions identified seven amyloidogenic sequences within the S-protein. All seven peptides in isolation formed aggregates during incubation at 37 °C. Three 20-amino acid long synthetic spike peptides (sequence 192-211, 601-620, 1166-1185) fulfilled three amyloid fibril criteria: nucleation dependent polymerization kinetics by ThT, Congo red positivity, and ultrastructural fibrillar morphology. Full-length folded S-protein did not form amyloid fibrils, but amyloid-like fibrils with evident branching were formed during 24 h of S-protein coincubation with the protease neutrophil elastase (NE) in vitro. NE efficiently cleaved S-protein, rendering exposure of amyloidogenic segments and accumulation of the amyloidogenic peptide 194-203, part of the most amyloidogenic synthetic spike peptide. NE is overexpressed at inflamed sites of viral infection. Our data propose a molecular mechanism for potential amyloidogenesis of SARS-CoV-2 S-protein in humans facilitated by endoproteolysis. The prospective of S-protein amyloidogenesis in COVID-19 disease associated pathogenesis can be important in understanding the disease and long COVID-19.

Biochemical Characterization of SARS-CoV-2 Spike RBD Mutations and Their Impact on ACE2 Receptor Binding

Infection of mammalian cells by SARS-CoV-2 coronavirus requires primary interaction between the receptor binding domain (RBD) of the viral spike protein and the host cell surface receptor angiotensin-converting enzyme 2 (ACE2) glycoprotein. Several mutations in the RBD of SARS-CoV-2 spike protein have been reported for several variants and resulted in wide spread of the COVID pandemic. For instance, the double mutations L452R and E484Q present in the Indian B.1.617 variant have been suggested to cause evasion of the host immune response. The common RBD mutations N501Y and E484K were found to enhance the interaction with the ACE2 receptor. In the current study, we analyzed the biosynthesis and secretion of the RBD double mutants L452R and E484Q in comparison to the wild-type RBD and the individual mutations N501 and E484K in mammalian cells. Moreover, we evaluated the interaction of these variants with ACE2 by means of expression of the S protein and co-immunoprecipitation with ACE2. Our results revealed that the double RBD mutations L452R and E484Q resulted in a higher expression level and secretion of spike S1 protein than other mutations. In addition, an increased interaction of these mutant forms with ACE2 in Calu3 cells was observed. Altogether, our findings highlight the impact of continuous S1 mutations on the pathogenicity of SARS-CoV-2 and provide further biochemical evidence for the dominance and high transmissibility of the double Indian mutations.

Convergent Evolution of Multiple Mutations Improves the Viral Fitness of SARS-CoV-2 Variants by Balancing Positive and Negative Selection

Multiple mutations have been seen to undergo convergent evolution in SARS-CoV-2 variants of concern. One such evolution occurs in Beta, Gamma, and Omicron variants at three amino acid positions K417, E484, and N501 in the receptor binding domain of the spike protein. We examined the physical mechanisms underlying the convergent evolution of three mutations K417T/E484K/N501Y by delineating the individual and collective effects of mutations on binding to angiotensin converting enzyme 2 receptor, immune escape from neutralizing antibodies, protein stability, and expression. Our results show that each mutation serves a distinct function that improves virus fitness supporting its positive selection, even though individual mutations have deleterious effects that make them prone to negative selection. Compared to the wild-type, K417T escapes Class 1 antibodies and has increased stability and expression; however, it has decreased receptor binding. E484K escapes Class 2 antibodies; however, it has decreased receptor binding, stability, and expression. N501Y increases receptor binding; however, it has decreased stability and expression. When these mutations come together, the deleterious effects are mitigated due to the presence of compensatory effects. Triple mutant K417T/E484K/N501Y has increased receptor binding, escapes both Class 1 and Class 2 antibodies, and has similar stability and expression as that of the wild-type. These results show that the convergent evolution of multiple mutations enhances viral fitness on different fronts by balancing both positive and negative selection and improves the chances of selection of mutations together.

Mutations and Evolution of the SARS-CoV-2 Spike Protein

The SARS-CoV-2 spike protein mediates target recognition, cellular entry, and ultimately the viral infection that leads to various levels of COVID-19 severities. Positive evolutionary selection of mutations within the spike protein has led to the genesis of new SARS-CoV-2 variants with greatly enhanced overall fitness. Given the trend of variants with increased fitness arising from spike protein alterations, it is critical that the scientific community understand the mechanisms by which these mutations alter viral functions. As of March 2022, five SARS-CoV-2 strains were labeled "variants of concern" by the World Health Organization: the Alpha, Beta, Gamma, Delta, and Omicron variants. This review summarizes the potential mechanisms by which the common mutations on the spike protein that occur within these strains enhance the overall fitness of their respective variants. In addressing these mutations within the context of the SARS-CoV-2 spike protein structure, spike/receptor binding interface, spike/antibody binding, and virus neutralization, we summarize the general paradigms that can be used to estimate the effects of future mutations along SARS-CoV-2 evolution.

COVID-19 Genetic Variants and Their Potential Impact in Vaccine Development

In the two years since the SARS-CoV-2 pandemic started, it has caused over 5 million deaths and 400 million infected cases, and the world continues to be on high alert for COVID-19. Among the variants of interest and concern of SARS-CoV-2, the current Omicron (B.1.1.529) and stealth Omicron (BA.2) raised serious concerns due to rapid rates of infection caused by numerous mutations in the spike protein, which could escape from the antibody-mediated neutralization and increase the risk of reinfections. Hence, this work aims to describe the most relevant mutations in the SARS-CoV-2 spike protein, discuss vaccine against variant of concerns, describe rare adverse events after COVID-19 vaccination, introduce the most available promising COVID-19 vaccine candidates, and provide few perspectives of the future variants.

A Detailed Overview of Immune Escape, Antibody Escape, Partial Vaccine Escape of SARS-CoV-2 and Their Emerging Variants With Escape Mutations

The infective SARS-CoV-2 is more prone to immune escape. Presently, the significant variants of SARS-CoV-2 are emerging in due course of time with substantial mutations, having the immune escape property. Simultaneously, the vaccination drive against this virus is in progress worldwide. However, vaccine evasion has been noted by some of the newly emerging variants. Our review provides an overview of the emerging variants' immune escape and vaccine escape ability. We have illustrated a broad view related to viral evolution, variants, and immune escape ability. Subsequently, different immune escape approaches of SARS-CoV-2 have been discussed. Different innate immune escape strategies adopted by the SARS-CoV-2 has been discussed like, IFN-I production dysregulation, cytokines related immune escape, immune escape associated with dendritic cell function and macrophages, natural killer cells and neutrophils related immune escape, PRRs associated immune evasion, and NLRP3 inflammasome associated immune evasion. Simultaneously we have discussed the significant mutations related to emerging variants and immune escape, such as mutations in the RBD region (N439K, L452R, E484K, N501Y, K444R) and other parts (D614G, P681R) of the S-glycoprotein. Mutations in other locations such as NSP1, NSP3, NSP6, ORF3, and ORF8 have also been discussed. Finally, we have illustrated the emerging variants' partial vaccine (BioNTech/Pfizer mRNA/Oxford-AstraZeneca/BBIBP-CorV/ZF2001/Moderna mRNA/Johnson & Johnson vaccine) escape ability. This review will help gain in-depth knowledge related to immune escape, antibody escape, and partial vaccine escape ability of the virus and assist in controlling the current pandemic and prepare for the next.

E484K and N501Y SARS-CoV 2 spike mutants Increase ACE2 recognition but reduce affinity for neutralizing antibody

SARS-CoV2 mutants B.1.1.7, B.1.351, and P.1 contain a key mutation N501Y. B.1.135 and P.1 lineages have another mutation, E484K. Here, we decode the effect of these two mutations on the host receptor, ACE2, and neutralizing antibody (B38) recognition. The N501Y RBD mutant binds to ACE2 with higher affinity due to improved π-π stacking and π-cation interactions. The higher binding affinity of the E484K mutant is caused due to the formation of additional hydrogen bond and salt-bridge interactions with ACE2. Both the mutants bind to the B38 antibody with reduced affinity due to the loss of several hydrogen-bonding interactions. The insights obtained from the study are crucial to interpret the increased transmissibility and reduced neutralization efficacy of rapidly emerging SARS-CoV2 VOCs.

The Nuts and Bolts of SARS-CoV-2 Spike Receptor-Binding Domain Heterologous Expression

COVID-19 is a highly infectious disease caused by a newly emerged coronavirus (SARS-CoV-2) that has rapidly progressed into a pandemic. This unprecedent emergency has stressed the significance of developing effective therapeutics to fight the current and future outbreaks. The receptor-binding domain (RBD) of the SARS-CoV-2 surface Spike protein is the main target for vaccines and represents a helpful "tool" to produce neutralizing antibodies or diagnostic kits. In this work, we provide a detailed characterization of the native RBD produced in three major model systems: Escherichia coli, insect and HEK-293 cells. Circular dichroism, gel filtration chromatography and thermal denaturation experiments indicated that recombinant SARS-CoV-2 RBD proteins are stable and correctly folded. In addition, their functionality and receptor-binding ability were further evaluated through ELISA, flow cytometry assays and bio-layer interferometry.