Many Roles of Carbohydrates: A Computational Spotlight on the Coronavirus S Protein Binding

Subtle Changes at the RBD/hACE2 Interface During SARS-CoV-2 Variant Evolution: A Molecular Dynamics Study

The SARS-CoV-2 Omicron variants show different behavior compared to the previous variants, especially with respect to the Delta variant, which promotes a lower morbidity despite being much more contagious. In this perspective, we performed molecular dynamics (MD) simulations of the different spike RBD/hACE2 complexes corresponding to the WT, Delta and four Omicron variants. Carrying out a comprehensive analysis of residue interactions within and between the two partners allowed us to draw the profile of each variant by using complementary methods (PairInt, hydrophobic potential, contact PCA). PairInt calculations highlighted the residues most involved in electrostatic interactions, which make a strong contribution to the binding with highly stable interactions between spike RBD and hACE2. Apolar contacts made a substantial and complementary contribution in Omicron with the detection of two hydrophobic patches. Contact networks and cross-correlation matrices were able to detect subtle changes at point mutations as the S375F mutation occurring in all Omicron variants, which is likely to confer an advantage in binding stability. This study brings new highlights on the dynamic binding of spike RBD to hACE2, which may explain the final persistence of Omicron over Delta.

Viral entry mechanisms: the role of molecular simulation in unlocking a key step in viral infections

Viral infections are a major global health concern, affecting millions of people each year. Viral entry is one of the crucial stages in the infection process, but its details remain elusive. Enveloped viruses are enclosed by a lipid membrane that protects their genetic material and these viruses are linked to various human illnesses, including influenza, and COVID-19. Due to the advancements made in the field of molecular simulation, significant progress has been made in unraveling the dynamic processes involved in viral entry of enveloped viruses. Simulation studies have provided deep insight into the function of the proteins responsible for attaching to the host receptors and promoting membrane fusion (fusion proteins), deciphering interactions between these proteins and receptors, and shedding light on the functional significance of key regions, such as the fusion peptide. These studies have already significantly contributed to our understanding of this critical aspect of viral infection and assisted the development of effective strategies to combat viral diseases and improve global health. This review focuses on the vital role of fusion proteins in facilitating the entry process of enveloped viruses and highlights the contributions of molecular simulation studies to uncover the molecular details underlying their mechanisms of action.

Impact of glycosylation on viral vaccines

Glycosylation is the most prominent modification important for vaccines and its specific pattern depends on several factors that need to be considered when developing a new biopharmaceutical. Tailor-made glycosylation can be exploited to develop more effective and safer vaccines; for this reason, a deep understanding of both glycoengineering strategies and glycans structures and functions is required. In this review we discuss the recent advances concerning glycoprotein expression systems and the explanation of glycans immunomodulation mechanisms. Furthermore, we highlight how glycans tune the immunological properties among different vaccines platforms (whole virus, recombinant protein, nucleic acid), also comparing commercially available formulations and describing the state-of-the-art analytical technologies for glycosylation analysis. The whole review stresses the aspect of glycoprotein glycans as a potential tool to overcome nowadays medical needs in vaccine field.

Enhanced Surface Accessibility of SARS-CoV-2 Omicron Spike Protein Due to an Altered Glycosylation Profile

SARS-CoV-2 spike (S) proteins undergo extensive glycosylation, aiding in proper folding, enhancing stability, and evading host immune surveillance. In this study, we used mass spectrometric analysis to elucidate the N-glycosylation characteristics and disulfide bonding of recombinant spike proteins derived from the SARS-CoV-2 Omicron variant (B.1.1.529) in comparison with the D614G spike variant. Furthermore, we conducted microsecond-long molecular dynamics simulations on spike proteins to resolve how the different N-glycans impact spike conformational sampling in the two variants. Our findings reveal that the Omicron spike protein maintains an overall resemblance to the D614G spike variant in terms of site-specific glycan processing and disulfide bond formation. Nonetheless, alterations in glycans were observed at certain N-glycosylation sites. These changes, in synergy with mutations within the Omicron spike protein, result in increased surface accessibility of the macromolecule, including the ectodomain, receptor-binding domain, and N-terminal domain. Additionally, mutagenesis and pull-down assays reveal the role of glycosylation of a specific sequon (N149); furthermore, the correlation of MD simulation and HDX-MS identified several high-dynamic areas of the spike proteins. These insights contribute to our understanding of the interplay between structure and function, thereby advancing effective vaccination and therapeutic strategies.

Mechanistic study of the transmission pattern of the SARS-CoV-2 omicron variant

The omicron variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) characterized by 30 mutations in its spike protein, has rapidly spread worldwide since November 2021, significantly exacerbating the ongoing COVID-19 pandemic. In order to investigate the relationship between these mutations and the variant's high transmissibility, we conducted a systematic analysis of the mutational effect on spike-angiotensin-converting enzyme-2 (ACE2) interactions and explored the structural/energy correlation of key mutations, utilizing a reliable coarse-grained model. Our study extended beyond the receptor-binding domain (RBD) of spike trimer through comprehensive modeling of the full-length spike trimer rather than just the RBD. Our free-energy calculation revealed that the enhanced binding affinity between the spike protein and the ACE2 receptor is correlated with the increased structural stability of the isolated spike protein, thus explaining the omicron variant's heightened transmissibility. The conclusion was supported by our experimental analyses involving the expression and purification of the full-length spike trimer. Furthermore, the energy decomposition analysis established those electrostatic interactions make major contributions to this effect. We categorized the mutations into four groups and established an analytical framework that can be employed in studying future mutations. Additionally, our calculations rationalized the reduced affinity of the omicron variant towards most available therapeutic neutralizing antibodies, when compared with the wild type. By providing concrete experimental data and offering a solid explanation, this study contributes to a better understanding of the relationship between theories and observations and lays the foundation for future investigations.

Impact of mutations on the plant-based production of recombinant SARS-CoV-2 RBDs

Subunit vaccines based on recombinant viral antigens are valuable interventions to fight existing and evolving viruses and can be produced at large-scale in plant-based expression systems. The recombinant viral antigens are often derived from glycosylated envelope proteins of the virus and glycosylation plays an important role for the immunogenicity by shielding protein epitopes. The receptor-binding domain (RBD) of the SARS-CoV-2 spike is a principal target for vaccine development and has been produced in plants, but the yields of recombinant RBD variants were low and the role of the N-glycosylation in RBD from different SARS-CoV-2 variants of concern is less studied. Here, we investigated the expression and glycosylation of six different RBD variants transiently expressed in leaves of Nicotiana benthamiana. All of the purified RBD variants were functional in terms of receptor binding and displayed almost full N-glycan occupancy at both glycosylation sites with predominately complex N-glycans. Despite the high structural sequence conservation of the RBD variants, we detected a variation in yield which can be attributed to lower expression and differences in unintentional proteolytic processing of the C-terminal polyhistidine tag used for purification. Glycoengineering towards a human-type complex N-glycan profile with core α1,6-fucose, showed that the reactivity of the neutralizing antibody S309 differs depending on the N-glycan profile and the RBD variant.