Exploring dynamics and network analysis of spike glycoprotein of SARS-COV-2

Examining sialic acid derivatives as potential inhibitors of SARS-CoV-2 spike protein receptor binding domain

Severe acute respiratory syndrome coronavirus 2 (SARS CoV-2) has been the primary reason behind the COVID-19 global pandemic which has affected millions of lives worldwide. The fundamental cause of the infection is the molecular binding of the viral spike protein receptor binding domain (SP-RBD) with the human cell angiotensin-converting enzyme 2 (ACE2) receptor. The infection can be prevented if the binding of RBD-ACE2 is resisted by utilizing certain inhibitors or drugs that demonstrate strong binding affinity towards the SP RBD. Sialic acid based glycans found widely in human cells and tissues have notable propensity of binding to viral proteins of the coronaviridae family. Recent experimental literature have used N-acetyl neuraminic acid (Sialic acid) to create diagnostic sensors for SARS-CoV-2, but a detailed interrogation of the underlying molecular mechanisms is warranted. Here, we perform all atom molecular dynamics (MD) simulations for the complexes of certain Sialic acid-based molecules with that of SP RBD of SARS CoV-2. Our results indicate that Sialic acid not only reproduces a binding affinity comparable to the RBD-ACE2 interactions, it also assumes the longest time to dissociate completely from the protein binding pocket of SP RBD. Our predictions corroborate that a combination of electrostatic and van der Waals energies as well the polar hydrogen bond interactions between the RBD residues and the inhibitors influence free energy of binding.Communicated by Ramaswamy H. Sarma.

Rapid simulation of glycoprotein structures by grafting and steric exclusion of glycan conformer libraries

Most membrane proteins are modified by covalent addition of complex sugars through N- and O-glycosylation. Unlike proteins, glycans do not typically adopt specific secondary structures and remain very mobile, shielding potentially large fractions of protein surface. High glycan conformational freedom hinders complete structural elucidation of glycoproteins. Computer simulations may be used to model glycosylated proteins but require hundreds of thousands of computing hours on supercomputers, thus limiting routine use. Here, we describe GlycoSHIELD, a reductionist method that can be implemented on personal computers to graft realistic ensembles of glycan conformers onto static protein structures in minutes. Using molecular dynamics simulation, small-angle X-ray scattering, cryoelectron microscopy, and mass spectrometry, we show that this open-access toolkit provides enhanced models of glycoprotein structures. Focusing on N-cadherin, human coronavirus spike proteins, and gamma-aminobutyric acid receptors, we show that GlycoSHIELD can shed light on the impact of glycans on the conformation and activity of complex glycoproteins.

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.

Multiscale Molecular Dynamics Simulations of the Homodimer Accessory Protein ORF7b of SARS-CoV-2

The SARS-CoV-2 ORF7b protein has drawn attention for its potential role in viral pathogenesis, but its structural details and lateral membrane associations remain elusive. In this study, we conducted multiscale molecular dynamics simulations to provide detailed molecular insights of the protein's dimerization, which is crucial for unraveling its structural model of protein-protein interface important to regulating cellular immune response. To gain a deeper understanding of homodimer configurations, we employed a machine learning algorithm for structural-based clustering. Clusters were categorized into three distinct groups for both parallel and antiparallel orientations, highlighting the influence of the initial monomer conformation on dimer configurations. Analysis of hydrogen bonding and π-π and π-cation stacking interactions within clusters revealed variations in interactions between clusters. In parallel dimers, weak stacking interactions in the transmembrane (TM) region were observed. In contrast, antiparallel dimers exhibited strong hydrogen bonding and stacking interactions contributing to tight dimeric packing, both within and outside the TM domain. Overall, our study provides a comprehensive view of the structural dynamics of ORF7b homodimerization in both parallel and antiparallel orientations. These findings shed light on the molecular interactions involved in ORF7b dimerization, which are crucial for understanding its potential roles in SARS-CoV-2 pathogenesis. This knowledge could inform future research and therapeutic strategies targeting this viral protein.

Evolving spike-protein N-glycosylation in SARS-CoV-2 variants

Since >3 years, SARS-CoV-2 has plunged humans into a colossal pandemic. Henceforth, multiple waves of infection have swept through the human population, led by variants that were able to partially evade acquired immunity. The co-evolution of SARS-CoV-2 variants with human immunity provides an excellent opportunity to study the interaction between viral pathogens and their human hosts. The heavily N-glycosylated spike-protein of SARS-CoV-2 plays a pivotal role in initiating infection and is the target for host immune-response, both of which are impacted by host-installed N-glycans. Using highly-sensitive DeGlyPHER approach, we compared the N-glycan landscape on spikes of the SARS-CoV-2 Wuhan-Hu-1 strain to seven WHO-defined variants of concern/interest, using recombinantly expressed, soluble spike-protein trimers, sharing same stabilizing-mutations. We found that N-glycan processing is conserved at most sites. However, in multiple variants, processing of N-glycans from high mannose- to complex-type is reduced at sites N165, N343 and N616, implicated in spike-protein function.

Conformational Fluctuations in β2-Microglubulin Using Markov State Modeling and Molecular Dynamics

Conformational dynamics in proteins can give rise to aggregation prone states during folding, and these kinetically stable states could form oligomers and aggregates. In this study, we investigate the intermediate states and near-folded states of β2-microglobulin and their physico-chemical properties using molecular dynamics and Markov state modeling. Analysis of hundreds of microseconds simulation show the importance of the edge strands in the misfolded states that give rise to a high exposure of hydrophobic residues in the core of the protein that could initiate oligomerization and aggregate formation. Our study sheds light on the first step of aggregation of β2m monomers and gave a better picture of the landscape of protein misfolding and aggregation.

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.

An overview of the vaccine platforms to combat COVID-19 with a focus on the subunit vaccines

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an emerging virus that has caused the recent coronavirus disease (COVID-19) global pandemic. The current approved COVID-19 vaccines have shown considerable efficiency against hospitalization and death. However, the continuation of the pandemic for more than two years and the likelihood of new strain emergence despite the global rollout of vaccination highlight the immediate need for the development and improvement of vaccines. mRNA, viral vector, and inactivated virus vaccine platforms were the first members of the worldwide approved vaccine list. Subunit vaccines. which are vaccines based on synthetic peptides or recombinant proteins, have been used in lower numbers and limited countries. The unavoidable advantages of this platform, including safety and precise immune targeting, make it a promising vaccine with wider global use in the near future. This review article summarizes the current knowledge on different vaccine platforms, focusing on the subunit vaccines and their clinical trial advancements against COVID-19.

Network analysis uncovers the communication structure of SARS-CoV-2 spike protein identifying sites for immunogen design

The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has triggered myriad efforts to understand the structure and dynamics of this complex pathogen. The spike glycoprotein of SARS-CoV-2 is a significant target for immunogens as it is the means by which the virus enters human cells, while simultaneously sporting mutations responsible for immune escape. These functional and escape processes are regulated by complex molecular-level interactions. Our study presents quantitative insights on domain and residue contributions to allosteric communication, immune evasion, and local- and global-level control of functions through the derivation of a weighted graph representation from all-atom MD simulations. Focusing on the ancestral form and the D614G-variant, we provide evidence of the utility of our approach by guiding the selection of a mutation that alters the spike's stability. Taken together, the network approach serves as a valuable tool to evaluate communication "hot-spots" in proteins to guide design of stable immunogens.

Glycosylation and the global virome

The sugars that coat the outsides of viruses and host cells are key to successful disease transmission, but they remain understudied compared to other molecular features. Understanding the comparative zoology of glycosylation - and harnessing it for predictive science - could help close the molecular gap in zoonotic risk assessment.

Nanobody derived using a peptide epitope from the spike protein receptor-binding motif inhibits entry of SARS-CoV-2 variants

The emergence of new escape mutants of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has escalated its penetration among the human population and has reinstated its status as a global pandemic. Therefore, developing effective antiviral therapy against emerging SARS-CoV variants and other viruses in a short period becomes essential. Blocking SARS-CoV-2 entry into human host cells by disrupting the spike glycoprotein-angiotensin-converting enzyme 2 interaction has already been exploited for vaccine development and monoclonal antibody therapy. Unlike the previous reports, our study used a nine-amino acid peptide from the receptor-binding motif of the spike protein as an epitope. We report the identification of an efficacious nanobody N1.2 that blocks the entry of pseudovirus-containing SARS-CoV-2 spike as the surface glycoprotein. Moreover, using mCherry fluorescence-based reporter assay, we observe a more potent neutralizing effect against both the hCoV19 (Wuhan/WIV04/2019) and the Omicron (BA.1) pseudotyped spike virus with a bivalent version of the N1.2 nanobody. In summary, our study presents a rapid and efficient methodology to use peptide sequences from a protein-receptor interaction interface as epitopes for screening nanobodies against potential pathogenic targets. We propose that this approach can also be widely extended to target other viruses and pathogens in the future.

Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2

Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.

Principles of SARS-CoV-2 Glycosylation

The structure and post-translational processing of the SARS-CoV-2 spike glycoprotein (S) is intimately associated with the function of the virus and of sterilising vaccines. The surface of the S protein is extensively modified by glycans, and their biosynthesis is driven by both the wider cellular context, and importantly, the underlining protein structure and local glycan density. Comparison of virally derived S protein with both recombinantly derived and adenovirally induced proteins, reveal hotspots of protein-directed glycosylation that drive conserved glycosylation motifs. Molecular dynamics simulations revealed that, while the S surface is extensively shielded by N-glycans, it presents regions vulnerable to neutralising antibodies. Furthermore, glycans have been shown to influence the accessibility of the receptor binding domain and the binding to the cellular receptor. The emerging picture is one of unifying, principles of S protein glycosylation and an intimate role of glycosylation in immunogen structure and efficacy.

Vaccines for COVID-19: A Systematic Review of Immunogenicity, Current Development, and Future Prospects

The persistent coronavirus disease 2019 (COVID-19), characterized by severe respiratory syndrome, is caused by coronavirus 2 (SARS-CoV-2), and it poses a major threat to public health all over the world. Currently, optimal COVID-19 management involves effective vaccination. Vaccination is known to greatly enhance immune response against viral infections and reduce public transmission of COVID-19. However, although current vaccines offer some benefits, viral variations and other factors demand the continuous development of vaccines to eliminate this virus from host. Hence, vaccine research and development is crucial and urgent to the elimination of this pandemic. Herein, we summarized the structural and replicatory features of SARS-CoV-2, and focused on vaccine-mediated disease prevention strategies like vaccine antigen selection, vaccine research, and vaccine application. We also evaluated the latest literature on COVID-19 and extensively reviewed action mechanisms, clinical trial (CT) progresses, advantages, as well as disadvantages of various vaccine candidates against SARS-CoV-2. Lastly, we discussed the current viral treatment, prevention trends, and future prospects.

Deciphering Isoniazid Drug Resistance Mechanisms on Dimeric Mycobacterium tuberculosis KatG via Post-molecular Dynamics Analyses Including Combined Dynamic Residue Network Metrics

Resistance mutations in Mycobacterium tuberculosis (Mtb) catalase peroxidase protein (KatG), an essential enzyme in isoniazid (INH) activation, reduce the sensitivity of Mtb to first-line drugs, hence presenting challenges in tuberculosis (TB) management. Thus, understanding the mutational imposed resistance mechanisms remains of utmost importance in the quest to reduce the TB burden. Herein, effects of 11 high confidence mutations in the KatG structure and residue network communication patterns were determined using extensive computational approaches. Combined traditional post-molecular dynamics analysis and comparative essential dynamics revealed that the mutant proteins have significant loop flexibility around the heme binding pocket and enhanced asymmetric protomer behavior with respect to wild-type (WT) protein. Heme contact analysis between WT and mutant proteins identified a reduction to no contact between heme and residue His270, a covalent bond vital for the heme-enabled KatG catalytic activity. Betweenness centrality calculations showed large hub ensembles with new hubs especially around the binding cavity and expanded to the dimerization domain via interface in the mutant systems, providing possible compensatory allosteric communication paths for the active site as a result of the mutations which may destabilize the heme binding pocket and the loops in its vicinity. Additionally, an interesting observation came from Eigencentrality hubs, most of which are located in the C-terminal domain, indicating relevance of the domain in the protease functionality. Overall, our results provide insight toward the mechanisms involved in KatG-INH resistance in addition to identifying key regions in the enzyme functionality, which can be used for future drug design.

Interaction of Spike protein and lipid membrane of SARS-CoV-2 with Ursodeoxycholic acid, an in-silico analysis

Numerous repositioned drugs have been sought to decrease the severity of SARS-CoV-2 infection. It is known that among its physicochemical properties, Ursodeoxycholic Acid (UDCA) has a reduction in surface tension and cholesterol solubilization, it has also been used to treat cholesterol gallstones and viral hepatitis. In this study, molecular docking was performed with the SARS-CoV-2 Spike protein and UDCA. In order to confirm this interaction, we used Molecular Dynamics (MD) in “SARS-CoV-2 Spike protein-UDCA”. Using another system, we also simulated MD with six UDCA residues around the Spike protein at random, naming this “SARS-CoV-2 Spike protein-6UDCA”. Finally, we evaluated the possible interaction between UDCA and different types of membranes, considering the possible membrane conformation of SARS-CoV-2, this was named “SARS-CoV-2 membrane-UDCA”. In the “SARS-CoV-2 Spike protein-UDCA”, we found that UDCA exhibits affinity towards the central region of the Spike protein structure of − 386.35 kcal/mol, in a region with 3 alpha helices, which comprises residues from K986 to C1032 of each monomer. MD confirmed that UDCA remains attached and occasionally forms hydrogen bonds with residues R995 and T998. In the presence of UDCA, we observed that the distances between residues atoms OG1 and CG2 of T998 in the monomers A, B, and C in the prefusion state do not change and remain at 5.93 ± 0.62 and 7.78 ± 0.51 Å, respectively, compared to the post-fusion state. Next, in “SARS-CoV-2 Spike protein-6UDCA”, the three UDCA showed affinity towards different regions of the Spike protein, but only one of them remained bound to the region between the region's heptad repeat 1 and heptad repeat 2 (HR1 and HR2) for 375 ps of the trajectory. The RMSD of monomer C was the smallest of the three monomers with a value of 2.89 ± 0.32, likewise, the smallest RMSF was also of the monomer C (2.25 ± 056). In addition, in the simulation of “SARS-CoV-2 membrane-UDCA”, UDCA had a higher affinity toward the virion-like membrane; where three of the four residues remained attached once they were close (5 Å, to the centre of mass) to the membrane by 30 ns. However, only one of them remained attached to the plasma-like membrane and this was in a cluster of cholesterol molecules. We have shown that UDCA interacts in two distinct regions of Spike protein sequences. In addition, UDCA tends to stay bound to the membrane, which could potentially reduce the internalization of SARS-CoV-2 in the host cell.

A Linkage-specific Sialic Acid Labeling Strategy Reveals Different Site-specific Glycosylation Patterns in SARS-CoV-2 Spike Protein Produced in CHO and HEK Cell Substrates

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus utilizes the extensively glycosylated spike (S) protein protruding from the viral envelope to bind to angiotensin-converting enzyme-related carboxypeptidase (ACE2) as its primary receptor to mediate host-cell entry. Currently, the main recombinant S protein production hosts are Chinese hamster ovary (CHO) and human embryonic kidney (HEK) cells. In this study, a recombinant S protein truncated at the transmembrane domain and engineered to express a C-terminal trimerization motif was transiently produced in CHO and HEK cell suspensions. To further evaluate the sialic acid linkages presenting on S protein, a two-step amidation process, employing dimethylamine and ammonium hydroxide reactions in a solid support system, was developed to differentially modify the sialic acid linkages on the glycans and glycopeptides from the S protein. The process also adds a charge to Asp and Glu which aids in ionization. We used MALDI-TOF and LC-MS/MS with electron-transfer/higher-energy collision dissociation (EThcD) fragmentation to determine global and site-specific N-linked glycosylation patterns. We identified 21 and 19 out of the 22 predicted N-glycosites of the SARS-CoV-2 S proteins produced in CHO and HEK, respectively. It was found that the N-glycosite at 1,158 position (N1158) and at 122, 282 and 1,158 positions (N122, N282 and N1158) were absent on S from CHO and HEK cells, respectively. The structural mapping of glycans of recombinant human S proteins reveals that CHO-Spike exhibits more complex and higher sialylation (α2,3-linked) content while HEK-Spike exhibits more high-mannose and a small amount of α2,3- and α2,6-linked sialic acids. The N74 site represents the most abundant glycosite on both spike proteins. The relatively higher amount of high-mannose abundant sites (N17, N234, N343, N616, N709, N717, N801, and N1134) on HEK-Spike suggests that glycan-shielding may differ among the two constructs. HEK-Spike can also provide different host immune system interaction profiles based on known immune system active lectins. Collectively, these data underscore the importance of characterizing the site-specific glycosylation of recombinant human spike proteins from HEK and CHO cells in order to better understand the impact of the production host on this complex and important protein used in research, diagnostics and vaccines.

Accelerating COVID-19 Research Using Molecular Dynamics Simulation

The COVID-19 pandemic has emerged as a global medico-socio-economic disaster. Given the lack of effective therapeutics against SARS-CoV-2, scientists are racing to disseminate suggestions for rapidly deployable therapeutic options, including drug repurposing and repositioning strategies. Molecular dynamics (MD) simulations have provided the opportunity to make rational scientific breakthroughs in a time of crisis. Advancements in these technologies in recent years have become an indispensable tool for scientists studying protein structure, function, dynamics, interactions, and drug discovery. Integrating the structural data obtained from high-resolution methods with MD simulations has helped in comprehending the process of infection and pathogenesis, as well as the SARS-CoV-2 maturation in host cells, in a short duration of time. It has also guided us to identify and prioritize drug targets and new chemical entities, and to repurpose drugs. Here, we discuss how MD simulation has been explored by the scientific community to accelerate and guide translational research on SARS-CoV-2 in the past year. We have also considered future research directions for researchers, where MD simulations can help fill the existing gaps in COVID-19 research.

Neuropilin-1 assists SARS-CoV-2 infection by stimulating the separation of Spike protein S1 and S2

The cell surface receptor Neuropilin-1 (Nrp1) was recently identified as a host factor for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) entry. The Spike protein of SARS-CoV-2 is cleaved into two segments, the S1 (residues (res.) 1-685) and the S2 (res. 686-1273) domains by furin protease. Nrp1 predominantly binds to the C-terminal RRAR amino acid motif (res. 682-685) of the S1 domain. In this study, we firstly modeled the association of an Nrp1 protein (consisting of domains a2-b1-b2) with the Spike protein. Next, we studied the separation of S2 from the S1 domain, with and without Nrp1 bound, by utilizing molecular dynamics pulling simulations. During the separation, Nrp1 stabilizes the S1 C-terminal region (res. 640-685) and thereby assists the detachment of S2 N-terminal region (res. 686-700). Without Nrp1 bound, S1 tends to become stretched, whereas the bound Nrp1 stimulates an earlier separation of S2 from the S1 domain. The liberated S2 domain is known to mediate the fusion of virus and host membranes; thus, Nrp1 likely increases virus infectivity by facilitating the S1 and S2 separation. We further analyzed the possible topological structure of the SARS-CoV-2 Spike protein when bound with Nrp1 and angiotensin-converting enzyme 2 (ACE2). Understanding of such an Nrp1-assisted viral infection opens the gate for the generation of protein-protein inhibitors, such as antibodies, which could attenuate the infection mechanism and protect certain cells in a future Nrp1-ACE2 targeted combination therapy.