Dynamic Interactions of Fully Glycosylated SARS-CoV-2 Spike Protein with Various Antibodies

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

Cancer cells and viruses share common glycoepitopes: exciting opportunities toward combined treatments

Aberrant glycosylation patterns of glycoproteins and glycolipids have long been recognized as one the major hallmarks of cancer cells that has led to numerous glycoconjugate vaccine attempts. These abnormal glycosylation profiles mostly originate from the lack of key glycosyltransferases activities, mutations, over expressions, or modifications of the requisite chaperone for functional folding. Due to their relative structural simplicity, O-linked glycans of the altered mucin family of glycoproteins have been particularly attractive in the design of tumor associated carbohydrate-based vaccines. Several such glycoconjugate vaccine formulations have generated potent monoclonal anti-carbohydrate antibodies useful as diagnostic and immunotherapies in the fight against cancer. Paradoxically, glycoproteins related to enveloped viruses also express analogous N- and O-linked glycosylation patterns. However, due to the fact that viruses are not equipped with the appropriate glycosyl enzyme machinery, they need to hijack that of the infected host cells. Although the resulting N-linked glycans are very similar to those of normal cells, some of their O-linked glycan patterns often share the common structural simplicity to those identified on tumor cells. Consequently, given that both cancer cells and viral glycoproteins share both common N- and O-linked glycoepitopes, glycoconjugate vaccines could be highly attractive to generate potent immune responses to target both conditions.

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

Glycosylation is one of the post-translational modifications with more than 50% of human proteins being glycosylated. The exact nature and chemical composition of glycans are inaccessible to X-ray or cryo-electron microscopy imaging techniques. Therefore, computational modeling studies and molecular dynamics must be used as a "computational microscope". The spike (S) protein of SARS-CoV-2 is heavily glycosylated, and a few glycans play a more functional role "beyond shielding". In this mini-review, we discuss computational investigations of the roles of specific S-protein and ACE2 glycans in the overall ACE2-S protein binding. We highlight different functions of specific glycans demonstrated in myriad computational models and simulations in the context of the SARS-CoV-2 virus binding to the receptor. We also discuss interactions between glycocalyx and the S protein, which may be utilized to design prophylactic polysaccharide-based therapeutics targeting the S protein. In addition, we underline the recent emergence of coronavirus variants and their impact on the S protein and its glycans.

Comparative Study of the Mutations Observed in the SARS-CoV-2 RBD Variants of Concern and Their Impact on the Interaction with the ACE2 Protein

SARS-CoV-2 strains have made an appearance across the globe, causing over 757 million cases and over 6.85 million deaths at the time of writing. The emergence of these variants shows the amplitude of genetic variation to which the wild-type strains have been subjected. The rise of the different SARS-CoV-2 variants resulting from such genetic modification has significantly affected COVD-19's major impact on proliferation, virulence, and clinics. With the emergence of the variants of concern, the spike protein has been identified as a possible therapeutic target due to its critical role in binding to human cells and pathogenesis. These mutations could be linked to functional heterogeneity and use a different infection strategy. For example, the Omicron variant's multiple mutations should be carefully examined, as they represent one of the most widely spread strains and hint to us that there may be more genetic changes in the virus. As a result, we applied a common protocol where we reconstructed SARS-CoV-2 variants of concern and performed molecular dynamics simulations to study the stability of the ACE2-RBD complex in each variant. We also carried out free energy calculations to compare the binding and biophysical properties of the different SARS-CoV-2 variants when they interact with ACE2. Therefore, we were able to obtain consistent results and uncover new crucial residues that were essential for preserving a balance between maintaining a high affinity for ACE2 and the capacity to evade RBD-targeted antibodies. Our detailed structural analysis showed that SARS-CoV-2 variants of concern show a higher affinity for ACE2 compared to the Wuhan strain. Additionally, residues K417N and E484K/A might play a crucial role in antibody evasion, whereas Q498R and N501Y are specifically mutated to strengthen RBD affinity to ACE2 and, thereby, increase the viral effect of the COVID-19 virus.

Establishing SARS-CoV-2 membrane protein-specific antibodies as a valuable serological target via high-content microscopy

The prevalence and strength of serological responses mounted toward SARS-CoV-2 proteins other than nucleocapsid (N) and spike (S), which may be of use as additional serological markers, remains underexplored. Using high-content microscopy to assess antibody responses against full-length StrepTagged SARS-CoV-2 proteins, we found that 85% (166/196) of unvaccinated individuals with RT-PCR confirmed SARS-CoV-2 infections and 74% (31/42) of individuals infected after being vaccinated developed detectable IgG against the structural protein M, which is higher than previous estimates. Compared with N antibodies, M IgG displayed a shallower time-dependent decay and greater specificity. Sensitivity for SARS-CoV-2 seroprevalence was enhanced when N and M IgG detection was combined. These findings indicate that screening for M seroconversion may be a good approach for detecting additional vaccine breakthrough infections and highlight the potential to use HCM as a rapidly deployable method to identify the most immunogenic targets of newly emergent pathogens.

CHARMM-GUI Membrane Builder: Past, Current, and Future Developments and Applications

Molecular dynamics simulations of membranes and membrane proteins serve as computational microscopes, revealing coordinated events at the membrane interface. As G protein-coupled receptors, ion channels, transporters, and membrane-bound enzymes are important drug targets, understanding their drug binding and action mechanisms in a realistic membrane becomes critical. Advances in materials science and physical chemistry further demand an atomistic understanding of lipid domains and interactions between materials and membranes. Despite a wide range of membrane simulation studies, generating a complex membrane assembly remains challenging. Here, we review the capability of CHARMM-GUI Membrane Builder in the context of emerging research demands, as well as the application examples from the CHARMM-GUI user community, including membrane biophysics, membrane protein drug-binding and dynamics, protein-lipid interactions, and nano-bio interface. We also provide our perspective on future Membrane Builder development.

Investigation of the effects of N-Acetylglucosamine on the stability of the spike protein in SARS-CoV-2 by molecular dynamics simulations

A lot of effort has been made in developing vaccine and therapeutic agents against the SARS-CoV-2, concentrating on the Spike protein that binds angiotensin-converting enzyme 2 on human cells. Nowadays, some researches study the role of the N-linked glycans as potential targets for vaccines and new agents. Due to the flexibility and diversity of the N-linked glycans, in this work, we focus on the N-Acetylglucosamine moiety, which is the precursor of nearly all eukaryotic glycans. We performed molecular dynamics simulations to study the effects of the N-Acetylglucosamine on the stability of the spike glycoprotein in SARS-CoV-2. After a 100 ns of simulation on the spike proteins without and with the N-Acetylglucosamine molecules, we found that the presence of N-Acetylglucosamine increases the local stability in their vicinity; even though their effect on the full structure is negligible. Thus; it can be inferred that the N-Acetylglucosamine moieties can potentially affect the interaction of the S protein with the ACE2 receptor. We also found that the S1 domain is more flexible than the S2 domain. We propose which of the experimentally observed glycans found on the spike may be more functional than the others. Detailed understanding of glycans is key for the development of new therapeutic strategies.

Molecular Epidemiology of the Norwegian SARS-CoV-2 Delta Lineage AY.63

Extensive genomic surveillance has given great insights into the evolution of the SARS-CoV-2 virus and emerging variants. During the summer months of 2021, Norway was dominated by the Pango lineage AY.63 which is a sub-lineage of the highly transmissible Delta variant. Strikingly, AY.63 did not spread in other countries to any significant extent. AY.63 carried a key mutation, A222V, in the spike protein, as well as the deletion of three residues in nsp1. Although these mutations are close to functionally important areas, we did not find any evidence that they induced higher fitness compared to other Delta lineages. This variant was introduced to Norway at a time when there were low levels of SARS-CoV-2 and contact-reducing measures were relaxed, which probably explains why the lineage rose so quickly. Furthermore, we found that the lack of imports of AY.63 from other countries probably led to the eventual demise of the lineage in Norway.

Thirty years of molecular dynamics simulations on posttranslational modifications of proteins

Posttranslational modifications (PTMs) are an integral component to how cells respond to perturbation. While experimental advances have enabled improved PTM identification capabilities, the same throughput for characterizing how structural changes caused by PTMs equate to altered physiological function has not been maintained. In this Perspective, we cover the history of computational modeling and molecular dynamics simulations which have characterized the structural implications of PTMs. We distinguish results from different molecular dynamics studies based upon the timescales simulated and analysis approaches used for PTM characterization. Lastly, we offer insights into how opportunities for modern research efforts on in silico PTM characterization may proceed given current state-of-the-art computing capabilities and methodological advancements.

Structural and Dynamic Insights into SARS-CoV-2 Spike-Protein-Montmorillonite Interactions

The spike (S) protein of SARS-CoV-2 has been found to play a decisive role in the cell entry mechanism of the virus and has been the prime target of most vaccine development efforts. Although numerous vaccines are already in use and more than half of the world population has been fully vaccinated, the emergence of new variants of the virus poses a challenge to the existing vaccines. Hence, developing an effective drug therapy is a crucial step in ending the pandemic. Nanoparticles can play a crucial role as a drug or a drug carrier and help tackle the pandemic effectively. Here, we performed explicit all-atom molecular dynamics simulations to probe interactions between S protein and Montmorillonite (MMT) nano clay surface. We built two systems with different counterions (Na+ and Ca2+), namely Na-MMT and Ca-MMT, to investigate the effect of different ions on S protein-MMT interaction. Structural modification of S protein was observed in the presence of MMT surface, particularly the loss of helical content of S protein. We revealed that electrostatic and hydrophobic interactions synergistically govern the S protein-MMT interactions. However, hydrophobic interactions were more pronounced in the Na-MMT system than in Ca-MMT. We also revealed residues and glycans of S protein closely interacting with the MMT surface. Interestingly, N165 and N343, which we found to be closely interacting with MMT in our simulations, also have a critical role in cell entry and in thwarting the cell's immune response in recent studies. Overall, our work provides atomistic insights into S protein-MMT interaction and enriches our understanding of the nanoparticle-S protein interaction mechanism, which will help develop advanced therapeutic techniques in the future.

The inherent flexibility of receptor binding domains in SARS-CoV-2 spike protein

Spike (S) protein is the primary antigenic target for neutralization and vaccine development for the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). It decorates the virus surface and undergoes large motions of its receptor binding domains (RBDs) to enter the host cell. Here, we observe Down, one-Up, one-Open, and two-Up-like structures in enhanced molecular dynamics simulations, and characterize the transition pathways via inter-domain interactions. Transient salt-bridges between RBDA and RBDC and the interaction with glycan at N343B support RBDA motions from Down to one-Up. Reduced interactions between RBDA and RBDB in one-Up induce RBDB motions toward two-Up. The simulations overall agree with cryo-electron microscopy structure distributions and FRET experiments and provide hidden functional structures, namely, intermediates along Down-to-one-Up transition with druggable cryptic pockets as well as one-Open with a maximum exposed RBD. The inherent flexibility of S-protein thus provides essential information for antiviral drug rational design or vaccine development.

Predicting spike protein NTD mutations of SARS-CoV-2 causing immune evasion by molecular dynamics simulations

The coronavirus disease 2019 (COVID-19) pandemic was caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Among all the potential targets studied for developing drugs and antibodies, the spike (S) protein is the most striking one, which is on the surface of the virus. In contrast with the intensively investigated immunodominant receptor-binding domain (RBD) of the protein, little is known about the neutralizing antibody binding mechanisms of the N-terminal domain (NTD), let alone the effects of NTD mutations on antibody binding and thereby the risk of immune evasion. Based on 400 ns molecular dynamics simulation for 11 NTD-antibody complexes together with other computational approaches in this study, we investigated critical residues for NTD-antibody binding and their detailed mechanisms. The results show that 36 residues on the NTD including R246, Y144, K147, Y248, L249 and P251 are critically involved in the direct interaction of the NTD with many monoclonal antibodies (mAbs), indicating that the viruses harboring these residue mutations may have a high risk of immune evasion. Binding free energy calculations and an interaction mechanism study reveal that R246I, which is present in the Beta (B.1.351/501Y.V2) variant, may have various impacts on current NTD antibodies through abolishing the hydrogen bonds and electrostatic interaction with the antibodies or affecting other interface residues. Therefore, special attention should be paid to the mutations of these key residues in future antibody and vaccine design and development.

Implications of the Immune Polymorphisms of the Host and the Genetic Variability of SARS-CoV-2 in the Development of COVID-19

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is responsible for the current pandemic affecting almost all countries in the world. SARS-CoV-2 is the agent responsible for coronavirus disease 19 (COVID-19), which has claimed millions of lives around the world. In most patients, SARS-CoV-2 infection does not cause clinical signs. However, some infected people develop symptoms, which include loss of smell or taste, fever, dry cough, headache, severe pneumonia, as well as coagulation disorders. The aim of this work is to report genetic factors of SARS-CoV-2 and host-associated to severe COVID-19, placing special emphasis on the viral entry and molecules of the immune system involved with viral infection. Besides this, we analyze SARS-CoV-2 variants and their structural characteristics related to the binding to polymorphic angiotensin-converting enzyme type 2 (ACE2). Additionally, we also review other polymorphisms as well as some epigenetic factors involved in the immunopathogenesis of COVID-19. These factors and viral variability could explain the increment of infection rate and/or in the development of severe COVID-19.