Biomechanical characterization of SARS-CoV-2 spike RBD and human ACE2 protein-protein interaction

Lectins in oncology and virology: Mechanisms of anticancer activity and SARS-CoV-2 inhibition

Lectins are proteins or glycoproteins of non-immune origin with carbohydrate-binding properties. They are found both prokaryotic and eukaryotic organisms. The most abundant source of the lectins are plants. Many lectins have anticancer effects by directly exerting cytotoxic effects on malignant cells or indirectly activating the immune system. Lectins also have antiviral activities. These proteins can recognise glycoproteins on the surface of enveloped viruses and bind to them. This creates a physical barrier between them and the corresponding receptors on the surface of the host cell, which prevents the virus from entering the cell and can thus effectively inhibit the replication of the virus. In this review, we focus on the anticancer activities of selected lectins and the underlying mechanisms. We also discuss different types of lectins with antiviral activity. We have paid special attention to lectins with inhibitory activity against SARS-CoV-2. Finally, we outline the challenges of using lectins in therapy and suggest future research directions.

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.

SARS-CoV-2 Spike Protein Enhances Carboxypeptidase Activity of Angiotensin-Converting Enzyme 2

In this study, we investigated whether severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) spike protein may modify angiotensin-converting enzyme 2 (ACE2) activity in the plasma, heart, kidney, liver, lung, and six brain regions (amygdala, brain stem, cortex, hippocampus, hypothalamus, and striatum) of diabetic and hypertensive rats. We determine ACE2 activity in the plasma and lysates of heart, kidney, liver, lung, and six brain regions. MLN-4760 inhibits ACE2 activity in the plasma and all organs. On the other hand, soluble ACE2 (sACE2) activity increased in the plasma of diabetic rats, and there was no change in the plasma of hypertensive rats. ACE2 activity was augmented in the liver, brain stem, and striatum, while it decreased in the kidney, amygdala, cortex, and hippocampus of diabetic rats. ACE2 activity increased in the kidney, liver, and lung, while it decreased in the heart, amygdala, cortex, and hypothalamus of hypertensive rats. We measured the ACE2 content via enzyme-linked immunosorbent assay and found that ACE2 protein levels increased in the heart, while it decreased in the plasma, kidney, brain stem, cortex, hippocampus, hypothalamus, and striatum of diabetic rats. ACE2 protein levels decreased in the brain stem, cortex, hippocampus, and hypothalamus of hypertensive rats. Our data showed that the spike protein enhanced ACE2 activity in the liver and lungs of diabetic rats, as well as in the heart and three of the brain regions (cortex, hypothalamus, and striatum) of hypertensive rats.

Ensemble-Based Mutational Profiling and Network Analysis of the SARS-CoV-2 Spike Omicron XBB Lineages for Interactions with the ACE2 Receptor and Antibodies: Cooperation of Binding Hotspots in Mediating Epistatic Couplings Underlies Binding Mechanism and Immune Escape

In this study, we performed a computational study of binding mechanisms for the SARS-CoV-2 spike Omicron XBB lineages with the host cell receptor ACE2 and a panel of diverse class one antibodies. The central objective of this investigation was to examine the molecular factors underlying epistatic couplings among convergent evolution hotspots that enable optimal balancing of ACE2 binding and antibody evasion for Omicron variants BA.1, BA2, BA.3, BA.4/BA.5, BQ.1.1, XBB.1, XBB.1.5, and XBB.1.5 + L455F/F456L. By combining evolutionary analysis, molecular dynamics simulations, and ensemble-based mutational scanning of spike protein residues in complexes with ACE2, we identified structural stability and binding affinity hotspots that are consistent with the results of biochemical studies. In agreement with the results of deep mutational scanning experiments, our quantitative analysis correctly reproduced strong and variant-specific epistatic effects in the XBB.1.5 and BA.2 variants. It was shown that Y453W and F456L mutations can enhance ACE2 binding when coupled with Q493 in XBB.1.5, while these mutations become destabilized when coupled with the R493 position in the BA.2 variant. The results provided a molecular rationale of the epistatic mechanism in Omicron variants, showing a central role of the Q493/R493 hotspot in modulating epistatic couplings between convergent mutational sites L455F and F456L in XBB lineages. The results of mutational scanning and binding analysis of the Omicron XBB spike variants with ACE2 receptors and a panel of class one antibodies provide a quantitative rationale for the experimental evidence that epistatic interactions of the physically proximal binding hotspots Y501, R498, Q493, L455F, and F456L can determine strong ACE2 binding, while convergent mutational sites F456L and F486P are instrumental in mediating broad antibody resistance. The study supports a mechanism in which the impact on ACE2 binding affinity is mediated through a small group of universal binding hotspots, while the effect of immune evasion could be more variant-dependent and modulated by convergent mutational sites in the conformationally adaptable spike regions.

SARS-CoV-2 Omicron Subvariants Do Not Differ Much in Binding Affinity to Human ACE2: A Molecular Dynamics Study

The emergence of the variant of concern Omicron (B.1.1.529) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) exacerbates the COVID-19 pandemic due to its high contagious ability. Studies have shown that the Omicron binds human ACE2 more strongly than the wild type. The prevalence of Omicron in new cases of COVID-19 promotes novel lineages with improved receptor binding affinity and immune evasion. To shed light on this open problem, in this work, we investigated the binding free energy of the receptor binding domain of the Omicron lineages BA.2, BA.2.3.20, BA.3, BA4/BA5, BA.2.75, BA.2.75.2, BA.4.6, XBB.1, XBB.1.5, BJ.1, BN.1, BQ.1.1, and CH.1.1 to human ACE2 using all-atom molecular dynamics simulation and the molecular mechanics Poisson-Boltzmann surface area method. The results show that these lineages have increased binding affinity compared to the BA.1 lineage, and BA.2.75 and BA.2.75.2 subvariants bind ACE2 more strongly than others. However, in general, the binding affinities of the Omicron lineages do not differ significantly from each other. The electrostatic force dominates over the van der Waals force in the interaction between Omicron lineages and human cells. Based on our results, we argue that viral evolution does not further improve the affinity of SARS-CoV-2 for ACE2 but may increase immune evasion.

Benchmarking the accuracy of structure-based binding affinity predictors on Spike-ACE2 deep mutational interaction set

Since the start of COVID-19 pandemic, a huge effort has been devoted to understanding the Spike (SARS-CoV-2)-ACE2 recognition mechanism. To this end, two deep mutational scanning studies traced the impact of all possible mutations across receptor binding domain (RBD) of Spike and catalytic domain of human ACE2. By concentrating on the interface mutations of these experimental data, we benchmarked six commonly used structure-based binding affinity predictors (FoldX, EvoEF1, MutaBind2, SSIPe, HADDOCK, and UEP). These predictors were selected based on their user-friendliness, accessibility, and speed. As a result of our benchmarking efforts, we observed that none of the methods could generate a meaningful correlation with the experimental binding data. The best correlation is achieved by FoldX (R = -0.51). When we simplified the prediction problem to a binary classification, that is, whether a mutation is enriching or depleting the binding, we showed that the highest accuracy is achieved by FoldX with a 64% success rate. Surprisingly, on this set, simple energetic scoring functions performed significantly better than the ones using extra evolutionary-based terms, as in Mutabind and SSIPe. Furthermore, we demonstrated that recent AI approaches, mmCSM-PPI and TopNetTree, yielded comparable performances to the force field-based techniques. These observations suggest plenty of room to improve the binding affinity predictors in guessing the variant-induced binding profile changes of a host-pathogen system, such as Spike-ACE2. To aid such improvements we provide our benchmarking data at https://github.com/CSB-KaracaLab/RBD-ACE2-MutBench with the option to visualize our mutant models at https://rbd-ace2-mutbench.github.io/.

Single-molecule force stability of the SARS-CoV-2-ACE2 interface in variants-of-concern

Mutations in SARS-CoV-2 have shown effective evasion of population immunity and increased affinity to the cellular receptor angiotensin-converting enzyme 2 (ACE2). However, in the dynamic environment of the respiratory tract, forces act on the binding partners, which raises the question of whether not only affinity but also force stability of the SARS-CoV-2-ACE2 interaction might be a selection factor for mutations. Using magnetic tweezers, we investigate the impact of amino acid substitutions in variants of concern (Alpha, Beta, Gamma and Delta) and on force-stability and bond kinetic of the receptor-binding domain-ACE2 interface at a single-molecule resolution. We find a higher affinity for all of the variants of concern (>fivefold) compared with the wild type. In contrast, Alpha is the only variant of concern that shows higher force stability (by 17%) compared with the wild type. Using molecular dynamics simulations, we rationalize the mechanistic molecular origins of this increase in force stability. Our study emphasizes the diversity of contributions to the transmissibility of variants and establishes force stability as one of the several factors for fitness. Understanding fitness advantages opens the possibility for the prediction of probable mutations, allowing a rapid adjustment of therapeutics, vaccines and intervention measures.

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.

Deciphering the free energy landscapes of SARS-CoV-2 wild type and Omicron variant interacting with human ACE2

The binding of the receptor binding domain (RBD) of the SARS-CoV-2 spike protein to the host cell receptor angiotensin-converting enzyme 2 (ACE2) is the first step in human viral infection. Therefore, understanding the mechanism of interaction between RBD and ACE2 at the molecular level is critical for the prevention of COVID-19, as more variants of concern, such as Omicron, appear. Recently, atomic force microscopy has been applied to characterize the free energy landscape of the RBD-ACE2 complex, including estimation of the distance between the transition state and the bound state, xu. Here, using a coarse-grained model and replica-exchange umbrella sampling, we studied the free energy landscape of both the wild type and Omicron subvariants BA.1 and XBB.1.5 interacting with ACE2. In agreement with experiment, we find that the wild type and Omicron subvariants have similar xu values, but Omicron binds ACE2 more strongly than the wild type, having a lower dissociation constant KD.

Chemically modified antiviral peptides against SARS-CoV-2

To date, the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) COVID-19 pandemic continues to be a potentially lethal disease. Although both vaccines and specific antiviral drugs have been approved, the search for more specific therapeutic approaches is still ongoing. The infection mechanism of SARS-CoV-2 consists of several stages, and each one can be selectively blocked to disrupt viral infection. Peptides are a promising class of antiviral compounds, which may be suitably modified to be more stable, more effective, and more selective towards a specific viral replication step. The latter two goals might be obtained by increasing the specificity and/or the affinity of the interaction with a specific target and often imply the stabilization of the secondary structure of the active peptide. This review is focused on modified antiviral peptides against SARS-CoV-2 acting at different stages of virus replication, including ACE2-RBD interaction, membrane fusion mechanism, and the proteolytic cleavage by different viral proteases. Therefore, the landscape presented herein provides a useful springboard for the design of new and powerful antiviral therapeutics.

Single-molecule scale quantification reveals interactions underlying protein-protein interface: from forces to non-covalent bonds

Protein-protein interactions (PPIs) between the B-cell lymphoma 2 (Bcl-2) family are considered a major driving force in cell cycle regulation and signaling. However, how this interfacial noncovalent interaction is achieved molecularly remains poorly understood. Herein, anti-apoptotic protein (Bcl-2) and pro-apoptotic protein (BAX) were used as models and their PPIs were explored for the first time using atomic force microscopy-based single-molecule force spectroscopy (SMFS) and in silico approaches. In addition, we used advanced analytical models, including multiple kinetic models, thermodynamic models, Poisson distributions, and contact angle molecular recognition to fully reveal the complexity of the BAX/Bcl-2 interaction interfaces. We propose that the binding kinetics between BAX/Bcl-2 are mainly mediated by specific (hydrogen bonding) and non-specific forces (hydrophobic interactions and electrostatic interactions) and show that the complicated multivalent binding interaction induces stable BAX/Bcl-2 complexes. This study enriches our understanding of the molecular mechanisms by which BAX interacts with Bcl-2. It provides valuable insights into the physical factors that need to be considered when designing PPI inhibitors.

Elucidation of N-/O-glycosylation and site-specific mapping of sialic acid linkage isomers of SARS-CoV-2 human receptor angiotensin-converting enzyme 2

Human angiotensin-converting enzyme 2 (hACE2) is the primary receptor for cellular entry of SARS-CoV-2 into human host cells. hACE2 is heavily glycosylated and glycans on the receptor may play a role in viral binding. Thus, comprehensive characterization of hACE2 glycosylation could aid our understanding of interactions between the receptor and SARS-CoV-2 spike (S) protein, as well as provide a basis for the development of therapeutic drugs targeting this crucial interaction. Herein, 138 N-glycan compositions were identified, most of which are complex-type N-glycans, from seven N-glycosites of hACE2. Among them, 67% contain at least one sialic acid residue. At the level of glycopeptides, the overall quantification of sialylated glycan isomers observed on the sites N322 and N546 have a higher degree of NeuAc (α2-3)Gal (over 80.3%) than that of other N-glycosites (35.6-71.0%). In terms of O-glycans, 69 glycan compositions from 12 O-glycosites were identified, and especially, the C-terminus of hACE2 is heavily O-glycosylated. The terminal sialic acid linkage type of H1N1S1 and H1N1S2 are covered highly with α2,3-sialic acid. These findings could aid the investigation of the interaction between SARS-CoV-2 and human host cells.

Biophysical investigation of interactions between SARS-CoV-2 spike protein and neuropilin-1

Recent studies have suggested that neuropilin-1 (NRP1) may serve as a potential receptor in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. However, the biophysical characteristics of interactions between NRP1 and SARS-CoV-2 remain unclear. In this study, we examined the interactions between NRP1 and various SARS-CoV-2 spike (S) fragments, including the receptor-binding domain (RBD) and the S protein trimer in a soluble form or expressed on pseudovirions, using atomic force microscopy and structural modeling. Our measurements shows that NRP1 interacts with the RBD and trimer at a higher binding frequency (BF) compared to ACE2. This NRP1-RBD interaction has also been predicted and simulated via AlphaFold2 and molecular dynamics simulations, and the results indicate that their binding patterns are very similar to RBD-ACE2 interactions. Additionally, under similar loading rates, the most probable unbinding forces between NRP1 and S trimer (both soluble form and on pseudovirions) are larger than the forces between NRP1 and RBD and between trimer and ACE2. Further analysis indicates that NRP1 has a stronger binding affinity to the SARS-CoV-2 S trimer with a dissociation rate of 0.87 s-1 , four times lower than the dissociation rate of 3.65 s-1 between NRP1 and RBD. Moreover, additional experiments show that RBD-neutralizing antibodies can significantly reduce the BF for both ACE2 and NRP1. Together, the study suggests that NRP1 can be an alternative receptor for SARS-CoV-2 attachment to human cells, and the neutralizing antibodies targeting SARS-CoV-2 RBD can reduce the binding between SARS-CoV-2 and NRP1.

Probing nanomechanical interactions of SARS-CoV-2 variants Omicron and XBB with common surfaces

The emergence of SARS-CoV-2 variants has further raised concerns about viral transmission. A fundamental understanding of the intermolecular interactions between the coronavirus and different surfaces is needed to address the transmission of SARS-CoV-2 through respiratory droplet-contaminated surfaces or fomites. The receptor-binding domain (RBD) of the spike protein is a key target for the adhesion of SARS-CoV-2 on the surface. To understand the effect of mutations on adhesion, atomic force microscopy-based single-molecule force spectroscopy (AFM-SMFS) was used to quantify the interactions between wild-type, Omicron, and XBB with several surfaces. The measurement revealed that RBD exhibits relatively higher forces on paper and gold surfaces, with the average force being 1.5 times greater compared to that on plastic surface. In addition, the force elevation on paper and gold surfaces for the variants can reach ∼28% relative to the wild type. These findings enhance our understanding of the nanomechanical interactions of the virus on common surfaces.

Interaction of SARS-CoV-2 with host cells and antibodies: experiment and simulation

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the devastating global COVID-19 pandemic announced by WHO in March 2020. Through unprecedented scientific effort, several vaccines, drugs and antibodies have been developed, saving millions of lives, but the fight against COVID-19 continues as immune escape variants of concern such as Delta and Omicron emerge. To develop more effective treatments and to elucidate the side effects caused by vaccines and therapeutic agents, a deeper understanding of the molecular interactions of SARS-CoV-2 with them and human cells is required. With special interest in computational approaches, we will focus on the structure of SARS-CoV-2 and the interaction of its spike protein with human angiotensin-converting enzyme-2 (ACE2) as a prime entry point of the virus into host cells. In addition, other possible viral receptors will be considered. The fusion of viral and human membranes and the interaction of the spike protein with antibodies and nanobodies will be discussed, as well as the effect of SARS-CoV-2 on protein synthesis in host cells.

Binding of synthetic nanobodies to the SARS-CoV-2 receptor-binding domain: the importance of salt bridges

In this study, five different SARS-CoV-2 receptor-binding domain (RBD) models were created based on the crystal structures of RBD complexes with two synthetic nanobodies (Sb16 and Sb45). Microsecond all-atom MD simulations revealed that Sb16 and Sb45 substantially stabilized the flexible RBD loop (residues GLU471-SER494) due to the salt bridges and hydrogen bonding interactions between RBD and the synthetic nanobodies. However, the calculation of binding free energy displayed that Sb45 had a higher binding affinity to RBD than Sb16, in agreement with the experimental result. This is because Sb45 has stronger electrostatic attraction to RBD as compared to Sb16. In particular, the salt bridge GLU484-ARG33 in Sb45-RBD is stronger than the GLU484-LYS32 in Sb16-RBD. Furthermore, by comparing the binding affinity of Sb16 for two RBD mutants (E484K and K417N), we found that E484K mutation substantially reduced the binding affinity to Sb16, and K417N mutation had no significant effect, qualitatively in agreement with experimental studies. According to the binding free energy calculation, the strong electrostatic repulsion between LYS32 and LYS484 caused by E484K mutation destroys the salt bridge between LYS32 and GLU484 in the RBD wild type (WT). In contrast, the binding of the K417N mutant to Sb16 effectively maintains the salt bridge between LYS32 and GLU484. Therefore, our research suggests that the salt bridges between RBD and synthetic nanobodies are crucial for binding synthetic nanobodies to RBD, and a SARS-CoV-2 variant can escape neutralization from nanobodies by creating electrostatic repulsion between them.

ACE2 in chronic disease and COVID-19: gene regulation and post-translational modification

Angiotensin-converting enzyme 2 (ACE2), a counter regulator of the renin-angiotensin system, provides protection against several chronic diseases. Besides chronic diseases, ACE2 is the host receptor for SARS-CoV or SARS-CoV-2 virus, mediating the first step of virus infection. ACE2 levels are regulated by transcriptional, post-transcriptional, and post-translational regulation or modification. ACE2 transcription is enhanced by transcription factors including Ikaros, HNFs, GATA6, STAT3 or SIRT1, whereas ACE2 transcription is reduced by the transcription factor Brg1-FoxM1 complex or ERRα. ACE2 levels are also regulated by histone modification or miRNA-induced destabilization. The protein kinase AMPK, CK1α, or MAP4K3 phosphorylates ACE2 protein and induces ACE2 protein levels by decreasing its ubiquitination. The ubiquitination of ACE2 is induced by the E3 ubiquitin ligase MDM2 or UBR4 and decreased by the deubiquitinase UCHL1 or USP50. ACE2 protein levels are also increased by the E3 ligase PIAS4-mediated SUMOylation or the methyltransferase PRMT5-mediated ACE2 methylation, whereas ACE2 protein levels are decreased by AP2-mediated lysosomal degradation. ACE2 is downregulated in several human chronic diseases like diabetes, hypertension, or lung injury. In contrast, SARS-CoV-2 upregulates ACE2 levels, enhancing host cell susceptibility to virus infection. Moreover, soluble ACE2 protein and exosomal ACE2 protein facilitate SARS-CoV-2 infection into host cells. In this review, we summarize the gene regulation and post-translational modification of ACE2 in chronic disease and COVID-19. Understanding the regulation and modification of ACE2 may help to develop prevention or treatment strategies for ACE2-mediated diseases.

S373P Mutation Stabilizes the Receptor-Binding Domain of the Spike Protein in Omicron and Promotes Binding

A cluster of several newly occurring mutations on Omicron is found at the β-core region of the spike protein's receptor-binding domain (RBD), where mutation rarely happened before. Notably, the binding of SARS-CoV-2 to human receptor ACE2 via RBD happens in a dynamic airway environment, where mechanical force caused by coughing or sneezing occurs. Thus, we used atomic force microscopy-based single-molecule force spectroscopy (AFM-SMFS) to measure the stability of RBDs and found that the mechanical stability of Omicron RBD increased by ∼20% compared with the wild type. Molecular dynamics (MD) simulations revealed that Omicron RBD showed more hydrogen bonds in the β-core region due to the closing of the α-helical motif caused primarily by the S373P mutation. In addition to a higher unfolding force, we showed a higher dissociation force between Omicron RBD and ACE2. This work reveals the mechanically stabilizing effect of the conserved mutation S373P for Omicron and the possible evolution trend of the β-core region of RBD.

Antiviral Activity against SARS-CoV-2 of Conformationally Constrained Helical Peptides Derived from Angiotensin-Converting Enzyme 2

Despite the availability of vaccines, COVID-19 continues to be aggressive, especially in immunocompromised individuals. Therefore, the development of a specific therapeutic agent with antiviral activity against SARS-CoV-2 is necessary. The infection pathway starts when the receptor binding domain of the viral spike protein interacts with the angiotensin converting enzyme 2 (ACE2), which acts as a host receptor for the RBD expressed on the host cell surface. In this scenario, ACE2 analogs binding to the RBD and preventing the cell entry can be promising antiviral agents. Most of the ACE2 residues involved in the interaction belong to the α1 helix, more specifically to the minimal fragment ACE2(24-42). In order to increase the stability of the secondary structure and thus antiviral activity, we designed different triazole-stapled analogs, changing the position and the number of bridges. The peptide called P3, which has the triazole-containing bridge in the positions 36-40, showed promising antiviral activity at micromolar concentrations assessed by plaque reduction assay. On the other hand, the double-stapled peptide P4 lost the activity, showing that excessive rigidity disfavors the interaction with the RBD.

Secondary Structures of MERS-CoV, SARS-CoV, and SARS-CoV-2 Spike Proteins Revealed by Infrared Vibrational Spectroscopy

All coronaviruses are characterized by spike glycoproteins whose S1 subunits contain the receptor binding domain (RBD). The RBD anchors the virus to the host cellular membrane to regulate the virus transmissibility and infectious process. Although the protein/receptor interaction mainly depends on the spike's conformation, particularly on its S1 unit, their secondary structures are poorly known. In this paper, the S1 conformation was investigated for MERS-CoV, SARS-CoV, and SARS-CoV-2 at serological pH by measuring their Amide I infrared absorption bands. The SARS-CoV-2 S1 secondary structure revealed a strong difference compared to those of MERS-CoV and SARS-CoV, with a significant presence of extended β-sheets. Furthermore, the conformation of the SARS-CoV-2 S1 showed a significant change by moving from serological pH to mild acidic and alkaline pH conditions. Both results suggest the capability of infrared spectroscopy to follow the secondary structure adaptation of the SARS-CoV-2 S1 to different environments.

De Novo Design of Anti-COVID Drugs Using Machine Learning-Based Equivariant Diffusion Model Targeting the Spike Protein

The drug discovery and research for an anti-COVID-19 drug has been ongoing despite repurposed drugs in the market. Over time, these drugs were discontinued due to side effects. The search for effective drugs is still under process. The role of Machine Learning (ML) is critical in the search for novel drug compounds. In the current work, using the equivariant diffusion model, we built novel compounds targeting the spike protein of SARS-CoV-2. Using the ML models, 196 de novo compounds were generated which had no hits on any major chemical databases. These novel compounds fulfilled all the criteria of ADMET properties to be lead-like and drug-like compounds. Of the 196 compounds, 15 were docked with high confidence in the target. These compounds were further subjected to molecular docking, the best compound having an IUPAC name of (4aS,4bR,8aS,8bS)-4a,8a-dimethylbiphenylene-1,4,5,8(4aH,4bH,8aH,8bH)-tetraone and a binding score of -6.930 kcal/mol. The principal compound is labeled as CoECG-M1. Density Function Theory (DFT) and Quantum optimization was carried out along with the study of ADMET properties. This suggests that the compound has potential drug-like properties. The docked complex was further subjected to MD simulations, GBSA, and metadynamics simulations to gain insights into the stability of binding. The model can be in the future modified to improve the positive docking rate.

Myricetin possesses the potency against SARS-CoV-2 infection through blocking viral-entry facilitators and suppressing inflammation in rats and mice

BACKGROUND

Myricetin (3,5,7-trihydroxy-2-(3,4,5-tri hydroxyphenyl)-4-benzopyrone) is a common flavonol extracted from many natural plants and Chinese herb medicines and has been demonstrated to have multiple pharmacological activities, such as anti-microbial, anti-thrombotic, neuroprotective, and anti-inflammatory effects. Previously, myricetin was reported to target Mpro and 3CL-Pro-enzymatic activity to SARS-CoV-2. However, the protective value of myricetin on SARS-Cov-2 infection through viral-entry facilitators has not yet been comprehensively understood.

PURPOSE

The aim of the current study was to evaluate the pharmacological efficacy and the mechanisms of action of myricetin against SARS-CoV-2 infection both in vitro and in vivo.

METHODS

The inhibitory effects of myricetin on SARS-CoV-2 infection and replication were assessed on Vero E6 cells. Molecular docking analysis and bilayer interferometry (BLI) assays, immunocytochemistry (ICC), and pseudoviruses assays were performed to evaluate the roles of myricetin in the intermolecular interaction between the receptor binding domain (RBD) of the SARS-CoV-2 spike (S) protein and angiotensin-converting enzyme 2 (ACE2). The anti-inflammatory potency and mechanisms of myricetin were examined in THP1 macrophages in vitro, as well as in carrageenan-induced paw edema, delayed-type hypersensitivity (DTH) induced auricle edema, and LPS-induced acute lung injury (ALI) animal models.

RESULTS

The results showed that myricetin was able to inhibit binding between the RBD of the SARS-CoV-2 S protein and ACE2 through molecular docking analysis and BLI assay, demonstrating its potential as a viral-entry facilitator blocker. Myricetin could also significantly inhibit SASR-CoV-2 infection and replication in Vero E6 cells (EC₅₀ 55.18 μM), which was further validated with pseudoviruses containing the RBD (wild-type, N501Y, N439K, Y453F) and an S1 glycoprotein mutant (S-D614G). Moreover, myricetin exhibited a marked suppressive action on the receptor-interacting serine/threonine protein kinase 1 (RIPK1)-driven inflammation and NF-kappa B signaling in THP1 macrophages. In animal model studies, myricetin notably ameliorated carrageenan-induced paw edema in rats, DTH induced auricle edema in mice, and LPS-induced ALI in mice.

CONCLUSION

Our findings showed that myricetin inhibited HCoV-229E and SARS-CoV-2 replication in vitro, blocked SARS-CoV-2 virus entry facilitators and relieved inflammation through the RIPK1/NF-κB pathway, suggesting that this flavonol has the potential to be developed as a therapeutic agent against COVID-19.

Rational approaches to discover SARS-CoV-2/ACE2 interaction inhibitors: Pharmacophore-based virtual screening, molecular docking, molecular dynamics and binding free energy studies

The lack of effective treatment remains a bottleneck in combating the current coronavirus family pandemic, particularly coronavirus 2019 (COVID-19), which is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The infection of host cells by SARS-CoV-2 is mediated by the binding of its receptor-binding domain (RBD) on the spike (S) glycoprotein to the host angiotensin-converting enzyme (ACE2) receptor. As all developed and available vaccines against COVID-19 do not provide long-term immunity, the creation of an effective drug for the treatment of COVID-19 is necessary and cannot be ignored. Therefore, the aim of this study is to present a computational screening method to identify potential inhibitor candidates with a high probability of blocking the binding of RBD to the ACE2 receptor. Pharmacophore mapping, molecular docking, molecular dynamics (MD) simulations, and binding free-energy analyses were performed to identify potential inhibitor candidates against ACE2/SARS-CoV-2. In conclusion, we propose the compound PubChem-84280085 as a potential inhibitor of protein-protein interactions to disrupt the binding of the SARS-CoV-2-RBD to the ACE2 receptor.

Interactions of SARS-CoV-2 and MERS-CoV fusion peptides measured using single-molecule force methods

We address the challenge of understanding how hydrophobic interactions are encoded by fusion peptide (FP) sequences within coronavirus (CoV) spike proteins. Within the FPs of severe acute respiratory syndrome CoV 2 and Middle East respiratory syndrome CoV (MERS-CoV), a largely conserved peptide sequence called FP1 (SFIEDLLFNK and SAIEDLLFDK in SARS-2 and MERS, respectively) has been proposed to play a key role in encoding hydrophobic interactions that drive viral-host cell membrane fusion. Although a non-polar triad (Leu-Leu-Phe (LLF)) is common to both FP1 sequences, and thought to dominate the encoding of hydrophobic interactions, FP1 from SARS-2 and MERS differ in two residues (Phe 2 versus Ala 2 and Asn 9 versus Asp 9, respectively). Here we explore whether single-molecule force measurements can quantify hydrophobic interactions encoded by FP1 sequences, and then ask whether sequence variations between FP1 from SARS-2 and MERS lead to significant differences in hydrophobic interactions. We find that both SARS-2 and MERS wild-type FP1 generate measurable hydrophobic interactions at the single-molecule level, but that SARS-2 FP1 encodes a substantially stronger hydrophobic interaction than its MERS counterpart (1.91 ± 0.03 nN versus 0.68 ± 0.03 nN, respectively). By performing force measurements with FP1 sequences with single amino acid substitutions, we determine that a single-residue mutation (Phe 2 versus Ala 2) causes the almost threefold difference in the hydrophobic interaction strength generated by the FP1 of SARS-2 versus MERS, despite the presence of LLF in both sequences. Infrared spectroscopy and circular dichroism measurements support the proposal that the outsized influence of Phe 2 versus Ala 2 on the hydrophobic interaction arises from variation in the secondary structure adopted by FP1. Overall, these insights reveal how single-residue diversity in viral FPs, including FP1 of SARS-CoV-2 and MERS-CoV, can lead to substantial changes in intermolecular interactions proposed to play a key role in viral fusion, and hint at strategies for regulating hydrophobic interactions of peptides in a range of contexts.

Distinct Conformations of SARS-CoV-2 Omicron Spike Protein and Its Interaction with ACE2 and Antibody

Since November 2021, Omicron has been the dominant severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant that causes the coronavirus disease 2019 (COVID-19) and has continuously impacted human health. Omicron sublineages are still increasing and cause increased transmission and infection rates. The additional 15 mutations on the receptor binding domain (RBD) of Omicron spike proteins change the protein conformation, enabling the Omicron variant to evade neutralizing antibodies. For this reason, many efforts have been made to design new antigenic variants to induce effective antibodies in SARS-CoV-2 vaccine development. However, understanding the different states of Omicron spike proteins with and without external molecules has not yet been addressed. In this review, we analyze the structures of the spike protein in the presence and absence of angiotensin-converting enzyme 2 (ACE2) and antibodies. Compared to previously determined structures for the wildtype spike protein and other variants such as alpha, beta, delta, and gamma, the Omicron spike protein adopts a partially open form. The open-form spike protein with one RBD up is dominant, followed by the open-form spike protein with two RBD up, and the closed-form spike protein with the RBD down. It is suggested that the competition between antibodies and ACE2 induces interactions between adjacent RBDs of the spike protein, which lead to a partially open form of the Omicron spike protein. The comprehensive structural information of Omicron spike proteins could be helpful for the efficient design of vaccines against the Omicron variant.

Binding of human ACE2 and RBD of Omicron enhanced by unique interaction patterns among SARS-CoV-2 variants of concern

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus causing COVID-19, has continued to mutate and spread worldwide despite global vaccination efforts. In particular, the Omicron variant, first identified in South Africa in late November 2021, has become the dominant strain worldwide. Compared to the original strain identified in Wuhan, Omicron features 50 genetic mutations, with 15 mutations in the receptor-binding domain (RBD) of the spike protein, which binds to the human angiotensin-converting enzyme 2 (ACE2) receptor for viral entry. However, it is not completely understood how these mutations alter the interaction and binding strength between the Omicron RBD and ACE2. In this study, we used a combined steered molecular dynamics (SMD) simulation and experimental microscale thermophoresis (MST) approach to quantify the interaction between Omicron RBD and ACE2. We report that the Omicron brings an enhanced RBD-ACE2 interface through N501Y, Q498R, and T478K mutations; the changes further lead to unique interaction patterns, reminiscing the features of previously dominated variants, Alpha (N501Y) and Delta (L452R and T478K). Among the Q493K and Q493R, we report that Q493R shows stronger binding to ACE2 than Q493K due to increased interactions. Our MST data confirmed that the Omicron mutations in RBD are associated with a five-fold higher binding affinity to ACE2 compared to the RBD of the original strain. In conclusion, our results could help explain the Omicron variant's prevalence in human populations, as higher interaction forces or affinity for ACE2 likely promote greater viral binding and internalization, leading to increased infectivity.

Novel perspective for protein-drug interaction analysis: atomic force microscope

Proteins are major drug targets, and drug-target interaction identification and analysis are important factors for drug discovery. Atomic force microscopy (AFM) is a powerful tool making it possible to image proteins with nanometric resolution and probe intermolecular forces under physiological conditions. We review recent studies conducted in the field of target protein drug discovery using AFM-based analysis technology, including drug-driven changes in nanomechanical properties of protein morphology and interactions. Underlying mechanisms (including thermodynamic and kinetic parameters) of the drug-target interaction and drug-modulating protein-protein interaction (PPI) on the surfaces of models or living cells are discussed. Furthermore, challenges and the outlook for the field are likewise discussed. Overall, this insight into the mechanical properties of protein-drug interactions provides an unprecedented information framework for rational drug discovery in the pharmaceutical field.

How cells wrap around virus-like particles using extracellular filamentous protein structures

Nanoparticles, such as viruses, can enter cells via endocytosis. During endocytosis, the cell surface wraps around the nanoparticle to effectively eat it. Prior focus has been on how nanoparticle size and shape impacts endocytosis. However, inspired by the noted presence of extracellular vimentin affecting viral and bacteria uptake, as well as the structure of coronaviruses, we construct a computational model in which both the cell-like construct and the virus-like construct contain filamentous protein structures protruding from their surfaces. We then study the impact of these additional degrees of freedom on viral wrapping. We find that cells with an optimal density of filamentous extracellular components (ECCs) are more likely to be infected as they uptake the virus faster and use relatively less cell surface area per individual virus. At the optimal density, the cell surface folds around the virus, and folds are faster and more efficient at wrapping the virus than crumple-like wrapping. We also find that cell surface bending rigidity helps generate folds, as bending rigidity enhances force transmission across the surface. However, changing other mechanical parameters, such as the stretching stiffness of filamentous ECCs or virus spikes, can drive crumple-like formation of the cell surface. We conclude with the implications of our study on the evolutionary pressures of virus-like particles, with a particular focus on the cellular microenvironment that may include filamentous ECCs.

How cells wrap around virus-like particles using extracellular filamentous protein structures

Nanoparticles, such as viruses, can enter cells via endocytosis. During endocytosis, the cell surface wraps around the nanoparticle to effectively eat it. Prior focus has been on how nanoparticle size and shape impacts endocytosis. However, inspired by the noted presence of extracellular vimentin affecting viral and bacteria uptake, as well as the structure of coronaviruses, we construct a computational model in which both the cell-like construct and the virus-like construct contain filamentous protein structures protruding from their surfaces. We then study the impact of these additional degrees of freedom on viral wrapping. We find that cells with an optimal density of filamentous extracellular components (ECCs) are more likely to be infected as they uptake the virus faster and use relatively less cell surface area per individual virus. At the optimal density, the cell surface folds around the virus, and folds are faster and more efficient at wrapping the virus than crumple-like wrapping. We also find that cell surface bending rigidity helps generate folds, as bending rigidity enhances force transmission across the surface. However, changing other mechanical parameters, such as the stretching stiffness of filamentous ECCs or virus spikes, can drive crumple-like formation of the cell surface. We conclude with the implications of our study on the evolutionary pressures of virus-like particles, with a particular focus on the cellular microenvironment that may include filamentous ECCs.

Bionics design of affinity peptide inhibitors for SARS-CoV-2 RBD to block SARS-CoV-2 RBD-ACE2 interactions

Coronavirus Disease 2019 (COVID-19), has already posed serious threats and impacts on the health of the population and the country's economy. Therefore, it is of great theoretical significance and practical application value to better understand the process of COVID-19 infection and develop effective therapeutic drugs. It is known that the receptor-binding structural domain (SARS-CoV-2 RBD) on the spike protein of the novel coronavirus directly mediates its interaction with the host receptor angiotensin-converting enzyme 2 (ACE2), and thus blocking SARS-CoV-2 RBD-ACE2 interaction is capable of inhibiting SARS-CoV-2 infection. Firstly, the interaction mechanism between SARS-CoV-2RBD-ACE2 was explored using molecular dynamics simulation (MD) coupled with molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) free energy calculation method. The results of energy analysis showed that the key residues R403, R408, K417, and Y505 of SARS-CoV-2 RBD and the key residues D30, E37, D38, and Y41 of ACE2 were identified. Therefore, according to the hotspot residues of ACE2 and their distribution, a short peptide library of high-affinity SARS-CoV-2 RBD was constructed. And by using molecular docking virtual screening, six short peptides including DDFEDY, DEFEDY, DEYEDY, DFVEDY, DFHEDY, and DSFEDY with high affinity for SARS-CoV-2 RBD were identified. The results of MD simulation further confirmed that DDFEDY, DEYEDY, and DFVEDY are expected to be effective inhibitors. Finally, the allergenicity, toxicity and solubility properties of the three peptide inhibitors were validated.

Force-tuned avidity of spike variant-ACE2 interactions viewed on the single-molecule level

Recent waves of COVID-19 correlate with the emergence of the Delta and the Omicron variant. We report that the Spike trimer acts as a highly dynamic molecular caliper, thereby forming up to three tight bonds through its RBDs with ACE2 expressed on the cell surface. The Spike of both Delta and Omicron (B.1.1.529) Variant enhance and markedly prolong viral attachment to the host cell receptor ACE2, as opposed to the early Wuhan-1 isolate. Delta Spike shows rapid binding of all three Spike RBDs to three different ACE2 molecules with considerably increased bond lifetime when compared to the reference strain, thereby significantly amplifying avidity. Intriguingly, Omicron (B.1.1.529) Spike displays less multivalent bindings to ACE2 molecules, yet with a ten time longer bond lifetime than Delta. Delta and Omicron (B.1.1.529) Spike variants enhance and prolong viral attachment to the host, which likely not only increases the rate of viral uptake, but also enhances the resistance of the variants against host-cell detachment by shear forces such as airflow, mucus or blood flow. We uncover distinct binding mechanisms and strategies at single-molecule resolution, employed by circulating SARS-CoV-2 variants to enhance infectivity and viral transmission.

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.

Free energy and kinetic rate calculation via non-equilibrium molecular simulation: application to biomolecules

Non-equilibrium molecular dynamics (NEMD) simulation has been recognized as a powerful tool for examining biomolecules and provides fruitful insights into not only non-equilibrium but also equilibrium processes. We review recent advances in NEMD simulation and relevant, fundamental results of non-equilibrium statistical mechanics. We first introduce Crooks fluctuation theorem and Jarzynski equality that relate free energy difference to work done on a physical system during a non-equilibrium process. The theorems are beneficial for the analysis of NEMD trajectories. We then describe rate theory, a framework to calculate molecular kinetics from a non-equilibrium process; this theoretical framework enables us to calculate a reaction time-mean-first passage time-from NEMD trajectories. We, in turn, present recent NEMD techniques that apply an external force to a system to enhance molecular dissociation and introduce their application to biomolecules. Lastly, we show the current status of an appropriate selection of reaction coordinates for NEMD simulation.

Prediction of infectivity of SARS-CoV2: Mathematical model with analysis of docking simulation for spike proteins and angiotensin-converting enzyme 2

Objectives

Variants of a coronavirus (SARS-CoV-2) have been spreading in a global pandemic. Improved understanding of the infectivity of future new variants is important so that effective countermeasures against them can be quickly undertaken. In our research reported here, we aimed to predict the infectivity of SARS-CoV-2 by using a mathematical model with molecular simulation analysis, and we used phylogenetic analysis to determine the evolutionary distance of the spike protein gene (S gene) of SARS-CoV-2.

Methods

We subjected the six variants and the wild type of spike protein and human angiotensin-converting enzyme 2 (ACE2) to molecular docking simulation analyses to understand the binding affinity of spike protein and ACE2. We then utilized regression analysis of the correlation coefficient of the mathematical model and the infectivity of SARS-CoV-2 to predict infectivity.

Results

The evolutionary distance of the S gene correlated with the infectivity of SARS-CoV-2 variants. The calculated biding affinity for the mathematical model obtained with results of molecular docking simulation also correlated with the infectivity of SARS-CoV-2 variants. These results suggest that the data from the docking simulation for the receptor binding domain of variant spike proteins and human ACE2 were valuable for prediction of SARS-CoV-2 infectivity.

Conclusion

We developed a mathematical model for prediction of SARS-CoV-2 variant infectivity by using binding affinity obtained via molecular docking and the evolutionary distance of the S gene.

Molecular Docking and Dynamics Simulation of Several Flavonoids Predict Cyanidin as an Effective Drug Candidate against SARS-CoV-2 Spike Protein

The in silico method has provided a versatile process of developing lead compounds from a large database in a short duration. Therefore, it is imperative to look for vaccinations and medications that can stop the havoc caused by SARS-CoV-2. The spike protein of SARS-CoV-2 is required for the viral entry into the host cells, hence inhibiting the virus from fusing and infecting the host. This study determined the binding interactions of 36 flavonoids along with two FDA-approved drugs against the spike protein receptor-binding domain of SARS-CoV-2 through molecular docking and molecular dynamics (MD) simulations. In addition, the molecular mechanics generalized Born surface area (MM/GBSA) approach was used to calculate the binding-free energy (BFE). Flavonoids were selected based on their in vitro assays on SARS-CoV and SARS-CoV-2. Our pharmacokinetics study revealed that cyanidin showed good drug-likeness, fulfilled Lipinski's rule of five, and conferred favorable toxicity parameters. Furthermore, MD simulations showed that cyanidin interacts with spike protein and alters the conformation and binding-free energy suited. Finally, an in vitro assay indicated that about 50% reduction in the binding of hACE2 with S1-RBD in the presence of cyanidin-containing red grapes crude extract was achieved at approximately 1.25 mg/mL. Hence, cyanidin may be a promising adjuvant medication for the SARS-CoV-2 spike protein based on in silico and in vitro research.

Computational analysis of interior mutations of SARS-CoV-2 Spike protein suggest a balance of protein stability and S2: S1 separation propensity

SARS-CoV-2 variants often include surface mutations in the Spike protein that are important for viruses to recognize host receptors and evade antibody neutralization. The Spike protein also has mutations in the interior of the protein likely to affect the Spike protein S1 - S2 subunit's separation propensity, the most important of which is the D614G mutation. Remarkably, the Omicron variant contains a large number of internal mutations at the S2: S1 interface, which have not been investigated yet. In this study, we examined the effects of such interfacial mutations on the S2: S1 and subunit domain interactions and on the subunit's dissociation process. We found that the interaction with S2 is mainly contributed by the three encapsulation domains, named INT, ED1 and ED2 of S1, which are sandwiched between the S1 RBD and N-terminal NTD domain. We found that D614 is the strongest contributor for the S2: S1 interaction which is greatly weakened by the D614G mutation. Surprisingly, we found that, mutations T547K, H655Y, N764K, N856K, N969K, L981F in the Omicron variant largely enhance the S2: ED1 interaction, partially compensating the loss of S2: ED2 interaction due to the D614G mutation. Lastly, these results, together with biological considerations, allow us to suggest that in addition to the binding strength of between the RBD and ACE2, the stability of the Spike protein and the propensity of Spike protein S2: S1 separation are critical factors which likely exist in a balance for a particular infectivity and pathogenicity of the virus.

CHARMM-GUI Enhanced Sampler for various collective variables and enhanced sampling methods

Enhanced sampling methodologies modifying underlying Hamiltonians can be used for the systems with a rugged potential energy surface that makes it hard to observe convergence using conventional unbiased molecular dynamics (MD) simulations. We present CHARMM-GUI Enhanced Sampler, a web-based tool to prepare various enhanced sampling simulations inputs with user-selected collective variables (CVs). Enhanced Sampler provides inputs for the following nine methods: accelerated MD, Gaussian accelerated MD, conformational flooding, metadynamics, adaptive biasing force, steered MD, temperature replica exchange MD, replica exchange solute tempering 2, and replica exchange umbrella sampling for the method-implemented MD packages including AMBER, CHARMM, GENESIS, GROMACS, NAMD, and OpenMM. Users only need to select a group of atoms via intuitive web-implementation in order to define commonly used nine CVs of interest: center of mass based distance, angle, dihedral, root-mean-square-distance, radius of gyration, distance projected on axis, two types of angles projected on axis, and coordination numbers. The enhanced sampling methods are tested with several biological systems to illustrate their efficiency over conventional MD. Enhanced Sampler with carefully optimized system-dependent parameters will help users to get meaningful results from their enhanced sampling simulations.

One-step asparaginyl endopeptidase (OaAEP1)-based protein immobilization for single-molecule force spectroscopy

Enzymatic protein ligation has become the most powerful and widely used method for high-precision atomic force microscopy single-molecule force spectroscopy (AFM-SMFS) study of protein mechanics. However, this methodology typically requires the functionalization of the glass surface with a corresponding peptide sequence/tag for enzymatic recognition and multiple steps are needed. Thus, it is time-consuming and a high level of experience is needed for reliable results. To solve this problem, we simplified the procedure using two strategies both based on asparaginyl endopeptidase (AEP). First, we designed a heterobifunctional peptide-based crosslinker, GL-peptide-propargylglycine, which links to an N 3-functionalized surface via the click reaction. Then, the target protein with a C-terminal NGL sequence can be immobilized via the AEP-mediated ligation. Furthermore, we took advantage of the direct ligation between primary amino in a small molecule and protein with C-terminal NGL by AEP. Thus, the target protein can be immobilized on an amino-functionalized surface via AEP in one step. Both approaches were successfully applied to the AFM-SMFS study of eGFP, showing consistent single-molecule results.

Mechanosensing view of SARS-CoV-2 infection by a DNA nano-assembly

The mechanical force between a virus and its host cell plays a critical role in viral infection. However, characterization of the virus-cell mechanical force at the whole-virus level remains a challenge. Herein, we develop a platform in which the virus is anchored with multivalence-controlled aptamers to achieve transfer of the virus-cell mechanical force to a DNA tension gauge tether (Virus-TGT). When the TGT is ruptured, the complex of binding module-virus-cell is detached from the substrate, accompanied by decreased host cell-substrate adhesion, thus revealing the mechanical force between whole-virus and cell. Using Virus-TGT, direct evidence about the biomechanical force between SARS-CoV-2 and the host cell is obtained. The relative mechanical force gap (<10 pN) at the cellular level between the wild-type virus to cell and a variant virus to cell is measured, suggesting a possible positive correlation between virus-cell mechanical force and infectivity. Overall, this strategy provides a new perspective to probe the SARS-CoV-2 mechanical force.

Interdomain Linker Effect on the Mechanical Stability of Ig Domains in Titin

Titin is the largest protein in humans, composed of more than one hundred immunoglobulin (Ig) domains, and plays a critical role in muscle’s passive elasticity. Thus, the molecular design of this giant polyprotein is responsible for its mechanical function. Interestingly, most of these Ig domains are connected directly with very few interdomain residues/linker, which suggests such a design is necessary for its mechanical stability. To understand this design, we chose six representative Ig domains in titin and added nine glycine residues (9G) as an artificial interdomain linker between these Ig domains. We measured their mechanical stabilities using atomic force microscopy-based single-molecule force spectroscopy (AFM-SMFS) and compared them to the natural sequence. The AFM results showed that the linker affected the mechanical stability of Ig domains. The linker mostly reduces its mechanical stability to a moderate extent, but the opposite situation can happen. Thus, this effect is very complex and may depend on each particular domain’s property.

All-Atom Simulations of Human ACE2-Spike Protein RBD Complexes for SARS-CoV-2 and Some of its Variants: Nature of Interactions and Free Energy Diagrams for Dissociation of the Protein Complexes

The spike protein of SARS-CoV-2 is known to interact with the human ACE2 protein via its receptor binding domain (RBD). We have investigated the molecular nature of this interprotein interaction and the associated free energy diagrams for the unbinding of the two proteins for SARS-CoV-2 and some of its known variants through all-atom simulations. The present work involves generation and analysis of 2.5 μs of unbiased and 4.2 μs of biased molecular dynamics trajectories in total for five explicitly solvated RBD-ACE2 systems at full atomic level. First, we have made a comparative analysis of the details of residue-wise specific interactions of the spike protein with ACE2 for SARS-CoV-1 and SARS-CoV-2. It is found that the average numbers of both direct interprotein and water-bridged hydrogen bonds between the RBD and ACE2 are higher for SARS-CoV-2 than SARS-CoV-1. These higher hydrogen bonded interactions are further aided by enhanced nonspecific electrostatic attractions between the two protein surfaces for SARS-CoV-2. The free energy calculations reveal that there is an increase in the free energy barrier by ∼1.5 kcal/mol for the unbinding of RBD from ACE2 for SARS-CoV-2 compared to that for SARS-CoV-1. Subsequently, we considered the RBDs of three variants of SARS-CoV-2, namely N501Y, E484Q/L452R, and N440K. The free energy barrier of protein unbinding for the N501Y variant is found to be ∼4 kcal/mol higher than the wild type SARS-CoV-2 which can be attributed to additional specific interactions involving Tyr501 of RBD and Lys353 and Tyr42 of ACE2 and also enhanced nonspecific electrostatic interaction between the protein surfaces. For the other two mutant variants of E484Q/L452R and N440K, the free energy barrier for protein unbinding increases by ∼2 and ∼1 kcal/mol, respectively, compared with the wild type SARS-CoV-2, which can be attributed to an increase in the number of interprotein hydrogen bonds for the former and also to enhanced positive electrostatic potential on the RBD surfaces for both of the variants. The successive breaking of interprotein hydrogen bonds along the free energy pathway of the unbinding process is also found out for all five systems studied here.

Impact of the temperature on the interactions between common variants of the SARS-CoV-2 receptor binding domain and the human ACE2

Several key mutations in the Spike protein receptor binding domain (RBD) have been identified to influence its affinity for the human Angiotensin-Converting Enzyme 2 (ACE2). Here, we perform a comparative study of the ACE2 binding to the wild type (Wuhan) RBD and some of its variants: Alpha B.1.1.7, Beta B.1.351, Delta B.1.617.2, Kappa B.1.617.1, B.1.1.7 + L452R and Omicron B.1.1.529. Using a coiled-coil mediated tethering approach of ACE2 in a novel surface plasmon resonance (SPR)-based assay, we measured interactions at different temperatures. Binding experiments at 10 °C enhanced the kinetic dissimilarities between the RBD variants and allowed a proper fit to a Langmuir 1:1 model with high accuracy and reproducibility, thus unraveling subtle differences within RBD mutants and ACE2 glycovariants. Our study emphasizes the importance of SPR-based assay parameters in the acquisition of biologically relevant data and offers a powerful tool to deepen our understanding of the role of the various RBD mutations in ACE2 interaction binding parameters.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

SARS-CoV-2 Omicron Variant Binds to Human Cells More Strongly than the Wild Type: Evidence from Molecular Dynamics Simulation

The emergence of the variant of concern Omicron (B.1.1.529) of the severe acute respiratory syndrome coronavirus 2 has aggravated the Covid-19 pandemic due to its very contagious ability. The high infection rate may be due to the high binding affinity of Omicron to human cells, but both experimental and computational studies have yielded conflicting results on this issue. Some studies have shown that the Omicron variant binds to human angiotensin-converting enzyme 2 (hACE2) more strongly than the wild type (WT), but other studies have reported comparable binding affinities. To shed light on this open problem, in this work, we calculated the binding free energy of the receptor binding domain (RBD) of the WT and Omicron spike protein to hACE2 using all-atom molecular dynamics simulation and the molecular mechanics Poisson-Boltzmann surface area method. We showed that Omicron binds to human cells more strongly than the WT due to increased RBD charge, which enhances electrostatic interaction with negatively charged hACE2. N440K, T478K, E484A, Q493R, and Q498R mutations in the RBD have been found to play a critical role in the stability of the RBD-hACE2 complex. The effect of homogeneous and heterogeneous models of glycans coating the viral RBD and the peptidyl domain of hACE2 was examined. Although the total binding free energy is not sensitive to the glycan model, the distribution of per-residue interaction energies depends on it. In addition, glycans have a little effect on the binding affinity of the WT RBD to hACE2.

Human interaction targets of SARS-COV-2 spike protein: A systematic review

Objectives:

The development of effective targeted therapy and drug-design approaches against the SARS-CoV-2 is a universal health priority. Therefore, it is important to assess possible therapeutic strategies against SARS-CoV-2 via its most interaction targets. The present study aimed to perform a systematic review on clinical and experimental investigations regarding SARS-COV-2 interaction targets for human cell entry.

Methods:

A systematic search using relevant MeSH terms and keywords was performed in PubMed, Scopus, Embase, and Web of Science (ISI) databases up to July 2021. Two reviewers independently assessed the eligibility of the studies, extracted the data, and evaluated the methodological quality of the included studies. Additionally, a narrative synthesis was done as a qualitative method for data gathering and synthesis of each outcome measure.

Results:

A total of 5610 studies were identified, and 128 articles were included in the systematic review. Based on the results, spike antigen was the only interaction protein from SARS-CoV-2. However, the interaction proteins from humans varied including different spike receptors and several cleavage enzymes. The most common interactions of the spike protein of SARS-CoV-2 for cell entry were ACE2 (entry receptor) and TMPRSS2 (for spike priming). A lot of published studies have mainly focused on the ACE2 receptor followed by the TMPRSS family and furin. Based on the results, ACE2 polymorphisms as well as spike RBD mutations affected the SARS-CoV-2 binding affinity.

Conclusion:

The included studies shed more light on SARS-CoV-2 cellular entry mechanisms and detailed interactions, which could enhance the understanding of SARS-CoV-2 pathogenesis and the development of new and comprehensive therapeutic approaches.

SARS-CoV-2 Delta Variant Decreases Nanobody Binding and ACE2 Blocking Effectivity

The Delta variant spreads more rapidly than previous variants of SARS-CoV-2. This variant comprises several mutations on the receptor-binding domain (RBD_Delta) of its spike glycoprotein, which binds to the peptidase domain (PD) of angiotensin-converting enzyme 2 (ACE2) receptors in host cells. The RBD-PD interaction has been targeted by antibodies and nanobodies to prevent viral infection, but their effectiveness against the Delta variant remains unclear. Here, we investigated RBD_Delta-PD interactions in the presence and absence of nanobodies H11-H4, H11-D4, and Ty1 by performing 21.8 μs of all-atom molecular dynamics simulations. Unbiased simulations revealed that Delta variant mutations strengthen RBD binding to ACE2 by increasing the hydrophobic interactions and salt bridge formation, but weaken interactions with H11-H4, H11-D4, and Ty1. Among these nanobodies H11-H4 and H11-D4 bind RBD without overlapping ACE2. They were unable to dislocate ACE2 from RBD_Delta when bound side by side with ACE2 on RBD. Steered molecular dynamics simulations at comparable loading rates to high-speed atomic force microscopy (AFM) experiments estimated lower rupture forces of the nanobodies from RBD_Delta compared to ACE2. Our results suggest that existing nanobodies are less effective to inhibit RBD_Delta-PD interactions and a new generation of nanobodies is needed to neutralize the Delta variant.

The nanocaterpillar's random walk: diffusion with ligand-receptor contacts

Particles with ligand-receptor contacts bind and unbind fluctuating "legs" to surfaces, whose fluctuations cause the particle to diffuse. Quantifying the diffusion of such "nanoscale caterpillars" is a challenge, since binding events often occur on very short time and length scales. Here we derive an analytical formula, validated by simulations, for the long time translational diffusion coefficient of an overdamped nanocaterpillar, under a range of modeling assumptions. We demonstrate that the effective diffusion coefficient, which depends on the microscopic parameters governing the legs, can be orders of magnitude smaller than the background diffusion coefficient. Furthermore it varies rapidly with temperature, and reproduces the striking variations seen in existing data and our own measurements of the diffusion of DNA-coated colloids. Our model gives insight into the mechanism of motion, and allows us to ask: when does a nanocaterpillar prefer to move by sliding, where one leg is always linked to the surface, and when does it prefer to move by hopping, which requires all legs to unbind simultaneously? We compare a range of systems (viruses, molecular motors, white blood cells, protein cargos in the nuclear pore complex, bacteria such as Escherichia coli, and DNA-coated colloids) and present guidelines to control the mode of motion for materials design.

A tethered ligand assay to probe SARS-CoV-2:ACE2 interactions

In the dynamic environment of the airways, where SARS-CoV-2 infections are initiated by binding to human host receptor ACE2, mechanical stability of the viral attachment is a crucial fitness advantage. Using single-molecule force spectroscopy techniques, we mimic the effect of coughing and sneezing, thereby testing the force stability of SARS-CoV-2 RBD:ACE2 interaction under physiological conditions. Our results reveal a higher force stability of SARS-CoV-2 binding to ACE2 compared to SARS-CoV-1, causing a possible fitness advantage. Our assay is sensitive to blocking agents preventing RBD:ACE2 bond formation. It will thus provide a powerful approach to investigate the modes of action of neutralizing antibodies and other agents designed to block RBD binding to ACE2 that are currently developed as potential COVID-19 therapeutics.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections are initiated by attachment of the receptor-binding domain (RBD) on the viral Spike protein to angiotensin-converting enzyme-2 (ACE2) on human host cells. This critical first step occurs in dynamic environments, where external forces act on the binding partners and avidity effects play an important role, creating an urgent need for assays that can quantitate SARS-CoV-2 interactions with ACE2 under mechanical load. Here, we introduce a tethered ligand assay that comprises the RBD and the ACE2 ectodomain joined by a flexible peptide linker. Using magnetic tweezers and atomic force spectroscopy as highly complementary single-molecule force spectroscopy techniques, we investigate the RBD:ACE2 interaction over the whole physiologically relevant force range. We combine the experimental results with steered molecular dynamics simulations and observe and assign fully consistent unbinding and unfolding events across the three techniques, enabling us to establish ACE2 unfolding as a molecular fingerprint. Measuring at forces of 2 to 5 pN, we quantify the force dependence and kinetics of the RBD:ACE2 bond in equilibrium. We show that the SARS-CoV-2 RBD:ACE2 interaction has higher mechanical stability, larger binding free energy, and a lower dissociation rate compared to SARS-CoV-1, which helps to rationalize the different infection patterns of the two viruses. By studying how free ACE2 outcompetes tethered ACE2, we show that our assay is sensitive to prevention of bond formation by external binders. We expect our results to provide a way to investigate the roles of viral mutations and blocking agents for targeted pharmaceutical intervention.

Prefusion spike protein conformational changes are slower in SARS-CoV-2 than in SARS-CoV-1

Within the last 2 decades, severe acute respiratory syndrome coronaviruses 1 and 2 (SARS-CoV-1 and SARS-CoV-2) have caused two major outbreaks; yet, for reasons not fully understood, the coronavirus disease 2019 pandemic caused by SARS-CoV-2 has been significantly more widespread than the 2003 SARS epidemic caused by SARS-CoV-1, despite striking similarities between these two viruses. The SARS-CoV-1 and SARS-CoV-2 spike proteins, both of which bind to host cell angiotensin-converting enzyme 2, have been implied to be a potential source of their differential transmissibility. However, the mechanistic details of prefusion spike protein binding to angiotensin-converting enzyme 2 remain elusive at the molecular level. Here, we performed an extensive set of equilibrium and nonequilibrium microsecond-level all-atom molecular dynamics simulations of SARS-CoV-1 and SARS-CoV-2 prefusion spike proteins to determine their differential dynamic behavior. Our results indicate that the active form of the SARS-CoV-2 spike protein is more stable than that of SARS-CoV-1 and the energy barrier associated with the activation is higher in SARS-CoV-2. These results suggest that not only the receptor-binding domain but also other domains such as the N-terminal domain could play a crucial role in the differential binding behavior of SARS-CoV-1 and SARS-CoV-2 spike proteins.