Up State of the SARS-COV-2 Spike Homotrimer Favors an Increased Virulence for New Variants

Bacterial Membrane Vesicles as a Novel Vaccine Platform against SARS-CoV-2

Throughout history, infectious diseases have plagued humanity, with outbreaks occurring regularly worldwide. Not every outbreak affects people globally; however, in the case of Coronavirus Disease 2019 (COVID-19), caused by a novel coronavirus (SARS-CoV-2), it reached a pandemic level within a remarkably short period. Fortunately, advancements in medicine and biotechnology have facilitated swift responses to the disease, resulting in the development of therapeutics and vaccines. Nevertheless, the persistent spread of the virus and the emergence of new variants underscore the necessity for protective interventions, leading researchers to seek more effective vaccines. Despite the presence of various types of vaccines, including mRNA and inactivated vaccines against SARS-CoV-2, new platforms have been investigated since the pandemic, and research on bacterial membrane vesicles (BMVs) has demonstrated their potential as a novel COVID-19 vaccine platform. Researchers have explored different strategies for BMV-based COVID-19 vaccines, such as mixing the vesicles with antigenic components of the virus due to their adjuvant capacity or decorating the vesicles with the viral antigens to create adjuvanted delivery systems. These approaches have presented promising results in inducing robust immune responses, but obstacles such as reproducibility in obtaining and homogeneous characterization of BMVs remain in developing vesicle-based vaccines. Overall, the development of BMV-based vaccines represents a novel and promising strategy in the fight against COVID-19. Additional research and clinical trials are needed to further evaluate the potential of these vaccines to offer long-lasting protection against SARS-CoV-2 and its evolving variants.

Analytical Ultracentrifugation Detects Quaternary Rearrangements and Antibody-Induced Conformational Selection of the SARS-CoV-2 Spike Trimer

Analytical ultracentrifugation (AUC) analysis shows that the SARS-CoV-2 trimeric Spike (S) protein adopts different quaternary conformations in solution. The relative abundance of the "open" and "close" conformations is temperature-dependent, and samples with different storage temperature history have different open/close distributions. Neutralizing antibodies (NAbs) targeting the S receptor binding domain (RBD) do not alter the conformer populations; by contrast, a NAb targeting a cryptic conformational epitope skews the Spike trimer toward an open conformation. The results highlight AUC, which is typically applied for molecular mass determination of biomolecules as a powerful tool for detecting functionally relevant quaternary protein conformations.

Differences between Omicron SARS-CoV-2 RBD and other variants in their ability to interact with cell receptors and monoclonal antibodies

SARS-CoV-2 remains a health threat with the continuous emergence of new variants. This work aims to expand the knowledge about the SARS-CoV-2 receptor-binding domain (RBD) interactions with cell receptors and monoclonal antibodies (mAbs). By using constant-pH Monte Carlo simulations, the free energy of interactions between the RBD from different variants and several partners (Angiotensin-Converting Enzyme-2 (ACE2) polymorphisms and various mAbs) were predicted. Computed RBD-ACE2-binding affinities were higher for two ACE2 polymorphisms (rs142984500 and rs4646116) typically found in Europeans which indicates a genetic susceptibility. This is amplified for Omicron (BA.1) and its sublineages BA.2 and BA.3. The antibody landscape was computationally investigated with the largest set of mAbs so far in the literature. From the 32 studied binders, groups of mAbs were identified from weak to strong binding affinities (e.g. S2K146). These mAbs with strong binding capacity and especially their combination are amenable to experimentation and clinical trials because of their high predicted binding affinities and possible neutralization potential for current known virus mutations and a universal coronavirus.Communicated by Ramaswamy H. Sarma.

Novel bispecific human antibody platform specifically targeting a fully open spike conformation potently neutralizes multiple SARS-CoV-2 variants

Rapid emergence of new variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has prompted an urgent need for the development of broadly applicable and potently neutralizing antibody platform against the SARS-CoV-2, which can be used for combatting the coronavirus disease 2019 (COVID-19). In this study, based on a noncompeting pair of phage display-derived human monoclonal antibodies (mAbs) specific to the receptor-binding domain (RBD) of SARS-CoV-2 isolated from human synthetic antibody library, we generated K202.B, a novel engineered bispecific antibody with an immunoglobulin G4-single-chain variable fragment design, with sub- or low nanomolar antigen-binding avidity. Compared with the parental mAbs or mAb cocktail, the K202.B antibody showed superior neutralizing potential against a variety of SARS-CoV-2 variants in vitro. Furthermore, structural analysis of bispecific antibody-antigen complexes using cryo-electron microscopy revealed the mode of action of K202.B complexed with a fully open three-RBD-up conformation of SARS-CoV-2 trimeric spike proteins by simultaneously interconnecting two independent epitopes of the SARS-CoV-2 RBD via inter-protomer interactions. Intravenous monotherapy using K202.B exhibited potent neutralizing activity in SARS-CoV-2 wild-type- and B.1.617.2 variant-infected mouse models, without significant toxicity in vivo. The results indicate that this novel approach of development of immunoglobulin G4-based bispecific antibody from an established human recombinant antibody library is likely to be an effective strategy for the rapid development of bispecific antibodies, and timely management against fast-evolving SARS-CoV-2 variants.

NS1 from Two Zika Virus Strains Differently Interact with a Membrane: Insights to Understand Their Differential Virulence

Zika virus (ZIKV) from Uganda (UG) expresses a phenotype related to fetal loss, whereas the variant from Brazil (BR) induces microcephaly in neonates. The differential virulence has a direct relation to biomolecular mechanisms that make one strain more aggressive than the other. The nonstructural protein 1 (NS1) is a key viral toxin to comprehend these viral discrepancies because of its versatility in many processes of the virus life cycle. Here, we aim to examine through coarse-grained models and molecular dynamics simulations the protein-membrane interactions for both NS1_ZIKV-UG and NS1_ZIKV-BR dimers. A first evaluation allowed us to establish that the NS1 proteins, in the membrane presence, explore new conformational spaces when compared to systems simulated without a lipid bilayer. These events derive from both differential coupling patterns and discrepant binding affinities to the membrane. The N-terminal domain, intertwined loop, and greasy finger proposed previously as binding membrane regions were also computationally confirmed by us. The anchoring sites have aromatic and ionizable residues that manage the assembly of NS1 toward the membrane, especially for the Ugandan variant. Furthermore, in the presence of the membrane, the difference in the dynamic cross-correlation of residues between the two strains is particularly pronounced in the putative epitope regions. This suggests that the protein-membrane interaction induces changes in the distal part related to putative epitopes. Taken together, these results open up new strategies for the treatment of flaviviruses that would specifically target these dynamic differences.

S Protein, ACE2 and Host Cell Proteases in SARS-CoV-2 Cell Entry and Infectivity; Is Soluble ACE2 a Two Blade Sword? A Narrative Review

Since the spread of the deadly virus SARS-CoV-2 in late 2019, researchers have restlessly sought to unravel how the virus enters the host cells. Some proteins on each side of the interaction between the virus and the host cells are involved as the major contributors to this process: (1) the nano-machine spike protein on behalf of the virus, (2) angiotensin converting enzyme II, the mono-carboxypeptidase and the key component of renin angiotensin system on behalf of the host cell, (3) some host proteases and proteins exploited by SARS-CoV-2. In this review, the complex process of SARS-CoV-2 entrance into the host cells with the contribution of the involved host proteins as well as the sequential conformational changes in the spike protein tending to increase the probability of complexification of the latter with angiotensin converting enzyme II, the receptor of the virus on the host cells, are discussed. Moreover, the release of the catalytic ectodomain of angiotensin converting enzyme II as its soluble form in the extracellular space and its positive or negative impact on the infectivity of the virus are considered.

Towards an optimal monoclonal antibody with higher binding affinity to the receptor-binding domain of SARS-CoV-2 spike proteins from different variants

A highly efficient and robust multiple scales in silico protocol, consisting of atomistic Molecular Dynamics (MD), coarse-grain (CG) MD, and constant-pH CG Monte Carlo (MC), has been developed and used to study the binding affinities of selected antigen-binding fragments of the monoclonal antibody (mAbs) CR3022 and several of its here optimized versions against 11 SARS-CoV-2 variants including the wild type. Totally 235,000 mAbs structures were initially generated using the RosettaAntibodyDesign software, resulting in top 10 scored CR3022-like-RBD complexes with critical mutations and compared to the native one, all having the potential to block virus-host cell interaction. Of these 10 finalists, two candidates were further identified in the CG simulations to be the best against all SARS-CoV-2 variants. Surprisingly, all 10 candidates and the native CR3022 exhibited a higher affinity for the Omicron variant despite its highest number of mutations. The multiscale protocol gives us a powerful rational tool to design efficient mAbs. The electrostatic interactions play a crucial role and appear to be controlling the affinity and complex building. Studied mAbs carrying a more negative total net charge show a higher affinity. Structural determinants could be identified in atomistic simulations and their roles are discussed in detail to further hint at a strategy for designing the best RBD binder. Although the SARS-CoV-2 was specifically targeted in this work, our approach is generally suitable for many diseases and viral and bacterial pathogens, leukemia, cancer, multiple sclerosis, rheumatoid, arthritis, lupus, and more.

Characterization and phylogenetic analysis of Iranian SARS-CoV-2 genomes: A phylogenomic study

Background and Aim

Characterization of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) based on analyzing the evolution and mutations of viruses is crucial for tracking viral infections, potential mutants, and other pathogens. The purpose was to study the complete sequences of SARS-CoV-2 to reveal genetic distance and mutation rate among different provinces of Iran.

Methods

As of March 2020-April 2021, a total of 131 SARS-CoV-2 whole genome sequences submitted from Tehran and 133 SARS-CoV-2 full-length sequences from 24 cities with high coverage submitted to EpiCoV GISAID database were analyzed to infer clades and mutation annotation compared with the wild-type variant Wuhan-Hu-1.

Results

The results of variant annotation were revealed 11,204 and 9468 distinct genomes were identified among the samples from different cities and Tehran, respectively. The phylogenetic analysis of genomic sequences showed the presence of eight GISAID clades, namely GH, GR, O, GRY, G, GK, L, and GV, and six Nextstrain clades; that is, 19A, 20A, 20B, 20I (alpha, V1), 20H (Beta, V2), and 21I (Delta) in Iran. The GH (GISAID clade), 20A (Nextstrain clade), and B.1 (Pango lineage) were predominant in Iran. Notably, analysis of the spike protein revealed D614G mutation (S_D614G) in 56% of the sequences. Also, the delta variant of the coronavirus, the super-infectious strain that was first identified among the sequences submitted from the southern cities of the country such as Zahedan, Yazd and Bushehr, and most likely from these places to other cities of Iran as well has expanded.

Conclusions

Our results indicate that most of the circulated viruses in Iran in the early time of the pandemic had collected in eight GISAID clades. Therefore, a continuous and extensive genome sequence analysis would be necessary to understand the genomic epidemiology of SARS-CoV-2 in Iran.

A bias of Asparagine to Lysine mutations in SARS-CoV-2 outside the receptor binding domain affects protein flexibility

Introduction

COVID-19 pandemic has been threatening public health and economic development worldwide for over two years. Compared with the original SARS-CoV-2 strain reported in 2019, the Omicron variant (B.1.1.529.1) is more transmissible. This variant has 34 mutations in its Spike protein, 15 of which are present in the Receptor Binding Domain (RBD), facilitating viral internalization via binding to the angiotensin-converting enzyme 2 (ACE2) receptor on endothelial cells as well as promoting increased immune evasion capacity.

Methods

Herein we compared SARS-CoV-2 proteins (including ORF3a, ORF7, ORF8, Nucleoprotein (N), membrane protein (M) and Spike (S) proteins) from multiple ancestral strains. We included the currently designated original Variant of Concern (VOC) Omicron, its subsequent emerged variants BA.1, BA2, BA3, BA.4, BA.5, the two currently emerging variants BQ.1 and BBX.1, and compared these with the previously circulating VOCs Alpha, Beta, Gamma, and Delta, to better understand the nature and potential impact of Omicron specific mutations.

Results

Only in Omicron and its subvariants, a bias toward an Asparagine to Lysine (N to K) mutation was evident within the Spike protein, including regions outside the RBD domain, while none of the regions outside the Spike protein domain were characterized by this mutational bias. Computational structural analysis revealed that three of these specific mutations located in the central core region, contribute to a preference for the alteration of conformations of the Spike protein. Several mutations in the RBD which have circulated across most Omicron subvariants were also analysed, and these showed more potential for immune escape.

Conclusion

This study emphasizes the importance of understanding how specific N to K mutations outside of the RBD region affect SARS-CoV-2 conformational changes and the need for neutralizing antibodies for Omicron to target a subset of conformationally dependent B cell epitopes.

Energetics of Spike Protein Opening of SARS-CoV-1 and SARS-CoV-2 and Its Variants of Concern: Implications in Host Receptor Scanning and Transmission

The receptor binding domain(s) (RBD) of spike (S) proteins of SARS-CoV-1 and SARS-CoV-2 (severe acute respiratory syndrome coronavirus) undergoes closed to open transition to engage with host ACE2 receptors. In this study, using multi atomistic (equilibrium) and targeted (non-equilibrium) molecular dynamics simulations, we have compared energetics of RBD opening pathways in full-length (modeled from cryo-EM structures) S proteins of SARS-CoV-1 and SARS-CoV-2. Our data indicate that amino acid variations at the RBD interaction interface can culminate into distinct free energy landscapes of RBD opening in these S proteins. We further report that mutations in the S protein of SARS-CoV-2 variants of concern can reduce the protein-protein interaction affinity of RBD(s) with its neighboring domains and could favor its opening to access ACE2 receptors. The findings can also aid in predicting the impact of future mutations on the rate of S protein opening for rapid host receptor scanning.

Electrostatic Features for the Receptor Binding Domain of SARS-COV-2 Wildtype and Its Variants. Compass to the Severity of the Future Variants with the Charge-Rule

Electrostatic intermolecular interactions are important in many aspects of biology. We have studied the main electrostatic features involved in the interaction of the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein with the human receptor Angiotensin-converting enzyme 2 (ACE2). As the principal computational tool, we have used the FORTE approach, capable to model proton fluctuations and computing free energies for a very large number of protein-protein systems under different physical-chemical conditions, here focusing on the RBD-ACE2 interactions. Both the wild-type and all critical variants are included in this study. From our large ensemble of extensive simulations, we obtain, as a function of pH, the binding affinities, charges of the proteins, their charge regulation capacities, and their dipole moments. In addition, we have calculated the pK_as for all ionizable residues and mapped the electrostatic coupling between them. We are able to present a simple predictor for the RBD-ACE2 binding based on the data obtained for Alpha, Beta, Gamma, Delta, and Omicron variants, as a linear correlation between the total charge of the RBD and the corresponding binding affinity. This "RBD charge rule" should work as a quick test of the degree of severity of the coming SARS-CoV-2 variants in the future.

Molecular dynamics studies reveal structural and functional features of the SARS-CoV-2 spike protein

The SARS-CoV-2 virus is responsible for the COVID-19 pandemic the world experience since 2019. The protein responsible for the first steps of cell invasion, the spike protein, has probably received the most attention in light of its central role during infection. Computational approaches are among the tools employed by the scientific community in the enormous effort to study this new affliction. One of these methods, namely molecular dynamics (MD), has been used to characterize the function of the spike protein at the atomic level and unveil its structural features from a dynamic perspective. In this review, we focus on these main findings, including spike protein flexibility, rare S protein conformational changes, cryptic epitopes, the role of glycans, drug repurposing, and the effect of spike protein variants.

The Local Topological Free Energy of the SARS-CoV-2 Spike Protein

The novel coronavirus SARS-CoV-2 infects human cells using a mechanism that involves binding and structural rearrangement of its Spike protein. Understanding protein rearrangement and identifying specific amino acids where mutations affect protein rearrangement has attracted much attention for drug development. In this manuscript, we use a mathematical method to characterize the local topology/geometry of the SARS-CoV-2 Spike protein backbone. Our results show that local conformational changes in the FP, HR1, and CH domains are associated with global conformational changes in the RBD domain. The SARS-CoV-2 variants analyzed in this manuscript (alpha, beta, gamma, delta Mink, G614, N501) show differences in the local conformations of the FP, HR1, and CH domains as well. Finally, most mutations of concern are either in or in the vicinity of high local topological free energy conformations, suggesting that high local topological free energy conformations could be targets for mutations with significant impact of protein function. Namely, the residues 484, 570, 614, 796, and 969, which are present in variants of concern and are targeted as important in protein function, are predicted as such from our model.

Self-association features of NS1 proteins from different flaviviruses

Flaviviruses comprise a large group of arboviral species that are distributed in several countries of the tropics, neotropics, and some temperate zones. Since they can produce neurological pathologies or vascular damage, there has been intense research seeking better diagnosis and treatments for their infections in the last decades. The flavivirus NS1 protein is a relevant clinical target because it is involved in viral replication, immune evasion, and virulence. Being a key factor in endothelial and tissue-specific modulation, NS1 has been largely studied to understand the molecular mechanisms exploited by the virus to reprogram host cells. A central part of the viral maturation processes is the NS1 oligomerization because many stages rely on these protein-protein assemblies. In the present study, the self-associations of NS1 proteins from Zika, Dengue, and West Nile viruses are examined through constant-pH coarse-grained biophysical simulations. Free energies of interactions were estimated for different oligomeric states and pH conditions. Our results show that these proteins can form both dimers and tetramers under conditions near physiological pH even without the presence of lipids. Moreover, pH plays an important role mainly controlling the regimes where van der Waals interactions govern their association. Finally, despite the similarity at the sequence level, we found that each flavivirus has a well-characteristic protein-protein interaction profile. These specific features can provide new hints for the development of binders both for better diagnostic tools and the formulation of new therapeutic drugs.

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.

VLP-Based COVID-19 Vaccines: An Adaptable Technology against the Threat of New Variants

Virus-like particles (VLPs) are a versatile, safe, and highly immunogenic vaccine platform. Recently, there are developmental vaccines targeting SARS-CoV-2, the causative agent of COVID-19. The COVID-19 pandemic affected humanity worldwide, bringing out incomputable human and financial losses. The race for better, more efficacious vaccines is happening almost simultaneously as the virus increasingly produces variants of concern (VOCs). The VOCs Alpha, Beta, Gamma, and Delta share common mutations mainly in the spike receptor-binding domain (RBD), demonstrating convergent evolution, associated with increased transmissibility and immune evasion. Thus, the identification and understanding of these mutations is crucial for the production of new, optimized vaccines. The use of a very flexible vaccine platform in COVID-19 vaccine development is an important feature that cannot be ignored. Incorporating the spike protein and its variations into VLP vaccines is a desirable strategy as the morphology and size of VLPs allows for better presentation of several different antigens. Furthermore, VLPs elicit robust humoral and cellular immune responses, which are safe, and have been studied not only against SARS-CoV-2 but against other coronaviruses as well. Here, we describe the recent advances and improvements in vaccine development using VLP technology.