Computational and comparative investigation of hydrophobic profile of spike protein of SARS-CoV-2 and SARS-CoV

The hydrophobic force is one of the most dominant factors in protein folding. A protein becomes functional only when it achieves its three-dimensional structure and stability upon folding. For a better understanding of the hydrophobic effects and their function in protein folding, quantitative measurement of the hydrophobicity of amino acid side chains is crucial. Spike protein is the primary structural protein in SARS-CoV-2 and SARS-CoV. This study explores how protein sequences in SARS-CoV-2 and SARS-CoV spike proteins encode hydrophobic interactions. Computational tools/techniques have been utilized to investigate the protein sequences of the spike proteins of SARS-CoV-2 and SARS-CoV. Investigations provided an estimate of hydrophobic distribution and its relative strength, indicating a hydrophobic pattern. Analysis of the spike protein''s hydrophobic profile may help identify and treat the virus-caused disease; additionally, it can give an insight into the transmissibility and pathogenicity of the virus.Supplementary InformationThe online version contains supplementary material available at 10.1007/s10867-022-09615-x.     

SARS-CoV-2 spike protein dictates syncytium-mediated lymphocyte elimination

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus is highly contagious and causes lymphocytopenia, but the underlying mechanisms are poorly understood. We demonstrate here that heterotypic cell-in-cell structures with lymphocytes inside multinucleate syncytia are prevalent in the lung tissues of coronavirus disease 2019 (COVID-19) patients. These unique cellular structures are a direct result of SARS-CoV-2 infection, as the expression of the SARS-CoV-2 spike glycoprotein is sufficient to induce a rapid (~45.1 nm/s) membrane fusion to produce syncytium, which could readily internalize multiple lines of lymphocytes to form typical cell-in-cell structures, remarkably leading to the death of internalized cells. This membrane fusion is dictated by a bi-arginine motif within the polybasic S1/S2 cleavage site, which is frequently present in the surface glycoprotein of most highly contagious viruses. Moreover, candidate anti-viral drugs could efficiently inhibit spike glycoprotein processing, membrane fusion, and cell-in-cell formation. Together, we delineate a molecular and cellular rationale for SARS-CoV-2 pathogenesis and identify novel targets for COVID-19 therapy.Subject terms: Cell biology, Molecular biology     

Structural basis of biased T cell receptor recognition of an immunodominant HLA-A2 epitope of the SARS-CoV-2 spike protein

CD8⁺ T cells play an important role in vaccination and immunity against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Although numerous SARS-CoV-2 CD8⁺ T cell epitopes have been identified, the molecular basis underpinning T cell receptor (TCR) recognition of SARS-CoV-2-specific T cells remains unknown. The T cell response directed toward SARS-CoV-2 spike protein–derived S^269–277 peptide presented by the human leukocyte antigen (HLA)-A∗02:01 allomorph (hereafter the HLA-A2^S269–277 epitope) is, to date, the most immunodominant SARS-CoV-2 epitope found in individuals bearing this allele. As HLA-A2^S269–277-specific CD8⁺ T cells utilize biased TRAV12 gene usage within the TCR α-chain, we sought to understand the molecular basis underpinning this TRAV12 dominance. We expressed four TRAV12⁺ TCRs which bound the HLA-A2^S269–277 complex with low micromolar affinity and determined the crystal structure of the HLA-A2^S269–277 binary complex, and subsequently a ternary structure of the TRAV12⁺ TCR complexed to HLA-A2^S269–277. We found that the TCR made extensive contacts along the entire length of the S^269–277 peptide, suggesting that the TRAV12⁺ TCRs would be sensitive to sequence variation within this epitope. To examine this, we investigated cross-reactivity toward analogous peptides from existing SARS-CoV-2 variants and closely related coronaviruses. We show via surface plasmon resonance and tetramer studies that the TRAV12⁺ T cell repertoire cross-reacts poorly with these analogous epitopes. Overall, we defined the structural basis underpinning biased TCR recognition of CD8⁺ T cells directed at an immunodominant epitope and provide a framework for understanding TCR cross-reactivity toward viral variants within the S^269–277 peptide.     

SARS-CoV-2 spike protein interacts with and activates TLR41

Dear Editor, Accumulating clinical data suggest the main causes of death by COVID-19 include respiratory failure and the onset of sepsis.1 Importantly,sepsis ha...     

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.     

Penta-peptide ATN-161 based neutralization mechanism of SARS-CoV-2 spike protein

SARS-CoV-2 has become a big challenge for the scientific community worldwide. SARS-CoV-2 enters into the host cell by the spike protein binding with an ACE2 receptor present on the host cell. Developing safe and effective inhibitor appears an urgent need to interrupt the binding of SARS-CoV-2 spike protein with ACE2 receptor in order to reduce the SARS-CoV-2 infection. We have examined the penta-peptide ATN-161 as potential inhibitor of ACE2 and SARS-CoV-2 spike protein binding, where ATN-161 has been commercially approved for the safety and possess high affinity and specificity towards the receptor binding domain (RBD) of S1 subunit in SARS-CoV-2 spike protein. We carried out experiments and confirmed these phenomena that the virus bindings were indeed minimized. ATN-161 peptide can be used as an inhibitor of protein-protein interaction (PPI) stands as a crucial interaction in biological systems. The molecular docking finding suggests that the binding energy of the ACE2-spike protein complex is reduced in the presence of ATN-161. Protein-protein docking binding energy (-40.50 kcal/mol) of the spike glycoprotein toward the human ACE2 and binding of ATN-161 at their binding interface reduced the biding energy (-26.25 kcal/mol). The finding of this study suggests that ATN-161 peptide can mask the RBD of the spike protein and be considered as a neutralizing candidate by binding with the ACE2 receptor. Peptide-based masking of spike S1 protein (RBD) and its neutralization is a highly promising strategy to prevent virus penetration into the host cell. Thus masking of the RBD leads to the loss of receptor recognition property which can reduce the chance of infection host cells.     

Receptor-binding domain of SARS-CoV-2 spike protein efficiently inhibits SARS-CoV-2 infection and attachment to mouse lung

COVID-19, caused by a novel coronavirus, SARS-CoV-2, poses a serious global threat. It was first reported in 2019 in China and has now dramatically spread across the world. It is crucial to develop therapeutics to mitigate severe disease and viral spread. The receptor-binding domains (RBDs) in the spike protein of SARS-CoV and MERS-CoV have shown anti-viral activity in previous reports suggesting that this domain has high potential for development as therapeutics. To evaluate the potential antiviral activity of recombinant SARS-CoV-2 RBD proteins, we determined the RBD residues of SARS-CoV-2 using a homology search with RBD of SARS-CoV. For efficient expression and purification, the signal peptide of spike protein was identified and used to generate constructs expressing recombinant RBD proteins. Highly purified RBD protein fused with the Fc domain of human IgG showed potent anti-viral efficacy, which was better than that of a protein fused with a histidine tag. Intranasally pre-administrated RBD protein also inhibited the attachment of SARS-COV-2 to mouse lungs. These findings indicate that RBD protein could be used for the prevention and treatment of SARS-CoV-2 infection.     

Signatures in SARS-CoV-2 spike protein conferring escape to neutralizing antibodies

Understanding SARS-CoV-2 evolution and host immunity is critical to control COVID-19 pandemics. At the core is an arms-race between SARS-CoV-2 antibody and angiotensin-converting enzyme 2 (ACE2) recognition, a function of the viral protein spike. Mutations in spike impacting antibody and/or ACE2 binding are appearing worldwide, imposing the need to monitor SARS-CoV2 evolution and dynamics in the population. Determining signatures in SARS-CoV-2 that render the virus resistant to neutralizing antibodies is critical. We engineered 25 spike-pseudotyped lentiviruses containing individual and combined mutations in the spike protein, including all defining mutations in the variants of concern, to identify the effect of single and synergic amino acid substitutions in promoting immune escape. We confirmed that E484K evades antibody neutralization elicited by infection or vaccination, a capacity augmented when complemented by K417N and N501Y mutations. In silico analysis provided an explanation for E484K immune evasion. E484 frequently engages in interactions with antibodies but not with ACE2. Importantly, we identified a novel amino acid of concern, S494, which shares a similar pattern. Using the already circulating mutation S494P, we found that it reduces antibody neutralization of convalescent and post-immunization sera, particularly when combined with E484K and with mutations able to increase binding to ACE2, such as N501Y. Our analysis of synergic mutations provides a signature for hotspots for immune evasion and for targets of therapies, vaccines and diagnostics.     

Structural effects of spike protein D614G mutation in SARS-CoV-2

A single mutation from aspartate to glycine at position 614 has dominated all circulating variants of the severe acute respiratory syndrome coronavirus 2. D614G mutation induces structural changes in the spike (S) protein that strengthen the virus infectivity. Here, we use molecular dynamics simulations to dissect the effects of mutation and 630-loop rigidification on S-protein structure. The introduction of the mutation orders the 630-loop structure and thereby induces global structural changes toward the cryoelectron microscopy structure of the D614G S-protein. The ordered 630-loop weakens local interactions between the 614^th residue and others in contrast to disordered structures in the wild-type protein. The mutation allosterically alters global interactions between receptor-binding domains, forming an asymmetric and mobile down conformation and facilitating transitions toward up conformation. The loss of salt bridge between D614 and K854 upon the mutation generally stabilizes S-protein protomer, including the fusion peptide proximal region that mediates membrane fusion. Understanding the molecular basis of D614G mutation is crucial as it dominates in all variants of concern, including Delta and Omicron.     

Glycyrrhizic acid exerts inhibitory activity against the spike protein of SARS-CoV-2

Coronavirus causes a disease with high infectivity and pathogenicity, especially SARS in 2003, MERS in 2012, and COVID-2019 currently. The spike proteins of these coronaviruses are critical for host cell entry by receptors. Thus, searching for broad-spectrum anti-coronavirus candidates, such as spike protein inhibitors, is vital and desirable due to the mutations in the spike protein. In this study, a combination of computer-aided drug design and biological verification was used to discover active monomers from traditional Chinese medicine. Surface plasmon resonance (SPR) assays and NanoBit assays were used to verify the predicated compounds with their binding activities to spike proteins and inhibitory activities on the SARS-CoV-2 RBD/ACE2 interaction, respectively. Furthermore, an MTT assay was used to evaluate the cell toxicities of active compounds. As a result, glycyrrhizic acid (ZZY-44) was found to be the most efficient and nontoxic broad-spectrum anti-coronavirus molecule in vitro, especially, the significant effect on SARS-CoV-2, which provided a theoretical basis for the study of the pharmacodynamic material basis of traditional Chinese medicine against SARS-CoV-2 and offered a lead compound for further structural modification in order to obtain more effective candidate drugs against SARS-CoV-2.     

Deamidation drives molecular aging of the SARS-CoV-2 spike protein receptor-binding motif

The spike protein is the main protein component of the SARS-CoV-2 virion surface. The spike receptor-binding motif mediates recognition of the human angiotensin-converting enzyme 2 receptor, a critical step in infection, and is the preferential target for spike-neutralizing antibodies. Posttranslational modifications of the spike receptor-binding motif have been shown to modulate viral infectivity and host immune response, but these modifications are still being explored. Here we studied asparagine deamidation of the spike protein, a spontaneous event that leads to the appearance of aspartic and isoaspartic residues, which affect both the protein backbone and its charge. We used computational prediction and biochemical experiments to identify five deamidation hotspots in the SARS-CoV-2 spike protein. Asparagine residues 481 and 501 in the receptor-binding motif deamidate with a half-life of 16.5 and 123 days at 37 °C, respectively. Deamidation is significantly slowed at 4 °C, indicating a strong dependence of spike protein molecular aging on environmental conditions. Deamidation of the spike receptor-binding motif decreases the equilibrium constant for binding to the human angiotensin-converting enzyme 2 receptor more than 3.5-fold, yet its high conservation pattern suggests some positive effect on viral fitness. We propose a model for deamidation of the full SARS-CoV-2 virion illustrating how deamidation of the spike receptor-binding motif could lead to the accumulation on the virion surface of a nonnegligible chemically diverse spike population in a timescale of days. Our findings provide a potential mechanism for molecular aging of the spike protein with significant consequences for understanding virus infectivity and vaccine development.     

An infectivity-enhancing site on the SARS-CoV-2 spike protein targeted by antibodies

1. Download : Download high-res image (204KB)
2. Download : Download full-size image

     

Enhanced sampling protocol to elucidate fusion peptide opening of SARS-CoV-2 spike protein

Large-scale conformational transitions in the spike protein S2 domain are required during host-cell infection of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus. Although conventional molecular dynamics simulations have been extensively used to study therapeutic targets of SARS-CoV-2, it is still challenging to gain molecular insight into the key conformational changes because of the size of the spike protein and the long timescale required to capture these transitions. In this work, we have developed an efficient simulation protocol that leverages many short simulations, a dynamic selection algorithm, and Markov state models to interrogate the structural changes of the S2 domain. We discovered that the conformational flexibility of the dynamic region upstream of the fusion peptide in S2 is coupled to the proteolytic cleavage state of the spike protein. These results suggest that opening of the fusion peptide likely occurs on a submicrosecond timescale after cleavage at the S2′ site. Building on the structural and dynamical information gained to date about S2 domain dynamics, we provide proof of principle that a small molecule bound to a seam neighboring the fusion peptide can slow the opening of the fusion peptide, leading to a new inhibition strategy for experiments to confirm. In aggregate, these results will aid the development of drug cocktails to inhibit infections caused by SARS-CoV-2 and other coronaviruses.     

The spike protein of SARS-CoV-2 induces heme oxygenase-1: Pathophysiologic implications

BackgroundAcute kidney injury (AKI) is both a consequence and determinant of outcomes in COVID-19. The kidney is one of the major organs infected by the causative virus, SARS-CoV-2. Viral entry into cells requires the viral spike protein, and both the virus and its spike protein appear in the urine of COVID-19 patients with AKI. We examined the effects of transfecting the viral spike protein of SARS-CoV-2 in kidney cell lines.MethodsHEK293, HEK293-ACE2⁺ (stably overexpressing ACE2), and Vero E6 cells having endogenous ACE2 were transfected with SARS-CoV-2 spike or control plasmid. Assessment of gene and protein expression, and syncytia formation was performed, and the effects of quercetin on syncytia formation examined.FindingsSpike transfection in HEK293-ACE2⁺ cells caused syncytia formation, cellular sloughing, and focal denudation of the cell monolayer; transfection in Vero E6 cells also caused syncytia formation. Spike expression upregulated potentially nephrotoxic genes (TNF-α, MCP-1, and ICAM1). Spike upregulated the cytoprotective gene HO-1 and relevant signaling pathways (p-Akt, p-STAT3, and p-p38). Quercetin, an HO-1 inducer, reduced syncytia formation and spike protein expression.InterpretationThe major conclusions of the study are: 1) Spike protein expression in kidney cells provides a relevant model for the study of maladaptive and adaptive responses germane to AKI in COVID-19; 2) such spike protein expression upregulates HO-1; and 3) quercetin, an HO-1 inducer, may provide a clinically relevant/feasible protective strategy in AKI occurring in the setting of COVID-19.FundingR01-DK119167 (KAN), R01-AI100911 (JPG), P30-DK079337; R01-DK059600 (AA).     

Structural and antigenic variations in the spike protein of emerging SARS-CoV-2 variants

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) virus is continuously evolving, and this poses a major threat to antibody therapies and currently authorized Coronavirus Disease 2019 (COVID-19) vaccines. It is therefore of utmost importance to investigate and predict the putative mutations on the spike protein that confer immune evasion. Antibodies are key components of the human immune system’s response to SARS-CoV-2, and the spike protein is a prime target of neutralizing antibodies (nAbs) as it plays critical roles in host cell recognition, fusion, and virus entry. The potency of therapeutic antibodies and vaccines partly depends on how readily the virus can escape neutralization. Recent structural and functional studies have mapped the epitope landscape of nAbs on the spike protein, which illustrates the footprints of several nAbs and the site of escape mutations. In this review, we discuss (1) the emerging SARS-CoV-2 variants; (2) the structural basis for antibody-mediated neutralization of SARS-CoV-2 and nAb classification; and (3) identification of the RBD escape mutations for several antibodies that resist antibody binding and neutralization. These escape maps are a valuable tool to predict SARS-CoV-2 fitness, and in conjunction with the structures of the spike-nAb complex, they can be utilized to facilitate the rational design of escape-resistant antibody therapeutics and vaccines.     

New discovery of high-affinity SARS-CoV-2 spike S2 protein binding peptide selected by PhIP-Seq

Highlights
1 A peptide Spep-1 targeting S2 of SARS-CoV-2 spike protein was selected by PhIP-Seq.
2 Spep-1 showed nanomolar affinity and high specificity to spike protein.
3 S-1 based immunoassay can detect femtomolar spike antigen in spiked serum samples.
4 Spep-1 can be used in future on S2 recognition, virus tracing and drug delivery.     

Antibody escape and cryptic cross-domain stabilization in the SARS-CoV-2 Omicron spike protein

1. Download : Download high-res image (246KB)
2. Download : Download full-size image

     

A potential interaction between the SARS-CoV-2 spike protein and nicotinic acetylcholine receptors

Changeux et al. (Changeux et al. C. R. Biol. 343:33–39.) recently suggested that the SARS-CoV-2 spike protein may interact with nicotinic acetylcholine receptors (nAChRs) and that such interactions may be involved in pathology and infectivity. This hypothesis is based on the fact that the SARS-CoV-2 spike protein contains a sequence motif similar to known nAChR antagonists. Here, we use molecular simulations of validated atomically detailed structures of nAChRs and of the spike to investigate the possible binding of the Y674-R685 region of the spike to nAChRs. We examine the binding of the Y674-R685 loop to three nAChRs, namely the human α4β2 and α7 subtypes and the muscle-like αβγδ receptor from Tetronarce californica. Our results predict that Y674-R685 has affinity for nAChRs. The region of the spike responsible for binding contains a PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. The conformational behavior of the bound Y674-R685 is highly dependent on the receptor subtype; it adopts extended conformations in the α4β2 and α7 complexes but is more compact when bound to the muscle-like receptor. In the α4β2 and αβγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation, similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket in which it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1, and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of simulations of the glycosylated spike show that the Y674-R685 region is accessible for binding. We suggest a potential binding orientation of the spike protein with nAChRs, in which they are in a nonparallel arrangement to one another.     

Structural diversity of the SARS-CoV-2 Omicron spike

1. Download : Download high-res image (279KB)
2. Download : Download full-size image