Binding mode of SARS-CoV-2 fusion peptide to human cellular membrane

The critical role of aromatic residues in the binding of the SARS-CoV-2 fusion peptide to phospholipid bilayer membranes

Based on the SARS-CoV-2 fusion peptide (FP) structure determined from the NMR experiment, we created six FP models under different environmental conditions to explore the effects of salt and cholesterol on FP-membrane binding. The all-atom molecular dynamics (MD) simulation results indicated that ionic environments notably impact the FP structure as well as the stability of the helical elements within the peptide. Our findings highlighted the unpredictable influence of ions on the secondary structures and dynamics of the FP, emphasizing the complexity and sensitivity of the peptide's conformations to ionic conditions. When exploring the peptide's interaction with a cholesterol-free phospholipid bilayer membrane, we found that the helical elements of the FP remain stable irrespective of the salt type (Na+ or Ca2+). This result emphasizes the crucial role of phospholipid bilayer membranes in supporting the secondary structures of the FP. The MD simulation results showed that Ca2+ ions facilitated deeper membrane penetration than Na+ ions, highlighting the critical role of calcium ions in the FP-membrane binding. Our study indicates the essential role of the aromatic residues (such as Phe833 and Tyr837) in the FP-membrane binding process. Finally, we investigated the FP-membrane binding patterns in the presence of cholesterol. The MD simulation results demonstrated that the coupling of Ca2+ ions and cholesterol would also benefit the FP-membrane binding. Furthermore, our findings reveal that while the type of ion and cholesterol content exert varied and unpredictable influences on FP-membrane binding patterns, aromatic residues like tyrosine (Tyr) and phenylalanine (Phe) play an essential role in FP-membrane binding. In particular, deep mutational scanning (DMS) experiments have confirmed that mutating phenylalanine in the FP significantly decreases viral mutational fitness, emphasizing the pivotal role of phenylalanine residues in membrane fusion. This knowledge can aid in developing more effective therapeutic strategies targeting the viral fusion peptide and its key amino acids, ultimately contributing to developing treatments and vaccines against the virus.

Positive residues of the SARS-CoV-2 fusion domain are key contributors to the initiation of membrane fusion

SARS-CoV-2 is one of the most infectious viruses ever recorded. Despite a plethora of research over the last several years, the viral life cycle is still not well understood, particularly membrane fusion. This process is initiated by the fusion domain (FD), a highly conserved stretch of amino acids consisting of a fusion peptide (FP) and fusion loop (FL), which in synergy perturbs the target cells' lipid membrane to lower the energetic cost necessary for fusion. In this study, through a mutagenesis-based approach, we have investigated the basic residues within the FD (K825, K835, R847, K854) utilizing an in vitro fusion assay and 19F NMR, validated by traditional 13C 15N techniques. Alanine and charge-conserving mutants revealed every basic residue plays a highly specific role within the mechanism of initiating fusion. Intriguingly, K825A led to increased fusogenecity which was found to be correlated to the number of amino acids within helix one, further implicating the role of this specific helix within the FD's fusion mechanism. This work has found basic residues to be important within the FDs fusion mechanism and highlights K825A, a specific mutation made within the FD of the SARS-CoV-2 spike protein, as requiring further investigation due to its potential to contribute to a more virulent strain of SARS-CoV-2.

Atomistic characterization of β2-glycoprotein I domain V interaction with anionic membranes

BACKGROUND

Interaction of β2-glycoprotein I (β2GPI) with anionic membranes is crucial in antiphospholipid syndrome (APS), implicating the role of its membrane-binding domain, domain V (DV). The mechanism of DV binding to anionic lipids is not fully understood.

OBJECTIVES

This study aimed to elucidate the molecular details of β2GPI DV binding to anionic membranes.

METHODS

We utilized molecular dynamics simulations to investigate the structural basis of anionic lipid recognition by DV. To corroborate the membrane-binding mode identified in the highly mobile membrane mimetic simulations, we conducted additional simulations using a full membrane model.

RESULTS

The study identified critical regions in DV, namely the lysine-rich loop and the hydrophobic loop, which are essential for membrane association via electrostatic and hydrophobic interactions, respectively. A novel lysine pair contributing to membrane binding was also discovered, providing new insights into β2GPI's membrane interaction. Simulations revealed 2 distinct binding modes of DV to the membrane, with mode 1 characterized by the insertion of the hydrophobic loop into the lipid bilayer, suggesting a dominant mechanism for membrane association. This interaction is pivotal for the pathogenesis of APS, as it facilitates the recognition of β2GPI by antiphospholipid antibodies.

CONCLUSION

The study advances our understanding of the molecular interactions between β2GPI's DV and anionic membranes, which are crucial for APS pathogenesis. It highlights the importance of specific regions in DV for membrane binding and reveals a predominant binding mode. These findings have significant implications for APS diagnostics and therapeutics, offering a deeper insight into the molecular basis of the syndrome.

Ebola Virus Glycoprotein Strongly Binds to Membranes in the Absence of Receptor Engagement

Ebola virus (EBOV) is an enveloped virus that must fuse with the host cell membrane in order to release its genome and initiate infection. This process requires the action of the EBOV envelope glycoprotein (GP), encoded by the virus, which resides in the viral envelope and consists of a receptor binding subunit, GP1, and a membrane fusion subunit, GP2. Despite extensive research, a mechanistic understanding of the viral fusion process is incomplete. To investigate GP-membrane association, a key step in the fusion process, we used two approaches: high-throughput measurements of single-particle diffusion and single-molecule measurements with optical tweezers. Using these methods, we show that the presence of the endosomal Niemann-Pick C1 (NPC1) receptor is not required for primed GP-membrane binding. In addition, we demonstrate this binding is very strong, likely attributed to the interaction between the GP fusion loop and the membrane's hydrophobic core. Our results also align with previously reported findings, emphasizing the significance of acidic pH in the protein-membrane interaction. Beyond Ebola virus research, our approach provides a powerful toolkit for studying other protein-membrane interactions, opening new avenues for a better understanding of protein-mediated membrane fusion events.

The carboxy terminus causes interfacial assembly of oleate hydratase on a membrane bilayer

The soluble flavoprotein oleate hydratase (OhyA) hydrates the 9-cis double bond of unsaturated fatty acids. OhyA substrates are embedded in membrane bilayers; OhyA must remove the fatty acid from the bilayer and enclose it in the active site. Here, we show that the positively charged helix-turn-helix motif in the carboxy terminus (CTD) is responsible for interacting with the negatively charged phosphatidylglycerol (PG) bilayer. Super-resolution microscopy of Staphylococcus aureus cells expressing green fluorescent protein fused to OhyA or the CTD sequence shows subcellular localization along the cellular boundary, indicating OhyA is membrane-associated and the CTD sequence is sufficient for membrane recruitment. Using cryo-electron microscopy, we solved the OhyA dimer structure and conducted 3D variability analysis of the reconstructions to assess CTD flexibility. Our surface plasmon resonance experiments corroborated that OhyA binds the PG bilayer with nanomolar affinity and we found the CTD sequence has intrinsic PG binding properties. We determined that the nuclear magnetic resonance structure of a peptide containing the CTD sequence resembles the OhyA crystal structure. We observed intermolecular NOE from PG liposome protons next to the phosphate group to the CTD peptide. The addition of paramagnetic MnCl2 indicated the CTD peptide binds the PG surface but does not insert into the bilayer. Molecular dynamics simulations, supported by site-directed mutagenesis experiments, identify key residues in the helix-turn-helix that drive membrane association. The data show that the OhyA CTD binds the phosphate layer of the PG surface to obtain bilayer-embedded unsaturated fatty acids.

SARS-CoV-2 fusion peptide sculpting of a membrane with insertion of charged and polar groups

The fusion peptide of SARS-CoV-2 spike is essential for infection. How this charged and hydrophobic domain occupies and affects membranes needs clarification. Its depth in zwitterionic, bilayered micelles at pH 5 (resembling late endosomes) was measured by paramagnetic NMR relaxation enhancements used to bias molecular dynamics simulations. Asp830 inserted deeply, along with Lys825 or Lys835. Protonation of Asp830 appeared to enhance agreement of simulated and NMR-measured depths. While the fusion peptide occupied a leaflet of the DMPC bilayer, the opposite leaflet invaginated with influx of water and choline head groups in around Asp830 and bilayer-inserted polar side chains. NMR-detected hydrogen exchange found corroborating hydration of the backbone of Thr827-Phe833 inserted deeply in bicelles. Pinching of the membrane at the inserted charge and the intramembrane hydration of polar groups agree with theory. Formation of corridors of hydrated, inward-turned head groups was accompanied by flip-flop of head groups. Potential roles of the defects are discussed.

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.

Cholesterol and Ceramide Facilitate Membrane Fusion Mediated by the Fusion Peptide of the SARS-CoV-2 Spike Protein

SARS-CoV-2 entry into host cells is mediated by the Spike (S) protein of the viral envelope. The S protein is composed of two subunits: S1 that induces binding to the host cell via its interaction with the ACE2 receptor of the cell surface and S2 that triggers fusion between viral and cellular membranes. Fusion by S2 depends on its heptad repeat domains that bring membranes close together and its fusion peptide (FP) that interacts with and perturbs the membrane structure to trigger fusion. Recent studies have suggested that cholesterol and ceramide lipids from the cell surface may facilitate SARS-CoV-2 entry into host cells, but their exact mode of action remains unknown. We have used a combination of in vitro liposome-liposome and in situ cell-cell fusion assays to study the lipid determinants of S-mediated membrane fusion. Our findings reveal that both cholesterol and ceramide lipids facilitate fusion, suggesting that targeting these lipids could be effective against SARS-CoV-2. As a proof of concept, we examined the effect of chlorpromazine (CPZ), an antipsychotic drug known to perturb membrane structure. Our results show that CPZ effectively inhibits S-mediated membrane fusion, thereby potentially impeding SARS-CoV-2 entry into the host cell.

Host Cell Membrane Capture by the SARS-CoV-2 Spike Protein Fusion Intermediate

Cell entry by SARS-CoV-2 is accomplished by the S2 subunit of the spike S protein on the virion surface by capture of the host cell membrane and fusion with the viral envelope. Capture and fusion require the prefusion S2 to transit to its potent fusogenic form, the fusion intermediate (FI). However, the FI structure is unknown, detailed computational models of the FI are unavailable, and the mechanisms and timing of membrane capture and fusion are not established. Here, we constructed a full-length model of the SARS-CoV-2 FI by extrapolating from known SARS-CoV-2 pre- and postfusion structures. In atomistic and coarse-grained molecular dynamics simulations the FI was remarkably flexible and executed giant bending and extensional fluctuations due to three hinges in the C-terminal base. The simulated configurations and their giant fluctuations are quantitatively consistent with SARS-CoV-2 FI configurations measured recently using cryo-electron tomography. Simulations suggested a host cell membrane capture time of ∼2 ms. Isolated fusion peptide simulations identified an N-terminal helix that directed and maintained binding to the membrane but grossly underestimated the binding time, showing that the fusion peptide environment is radically altered when attached to its host fusion protein. The large configurational fluctuations of the FI generated a substantial exploration volume that aided capture of the target membrane, and may set the waiting time for fluctuation-triggered refolding of the FI that draws the viral envelope and host cell membrane together for fusion. These results describe the FI as machinery that uses massive configurational fluctuations for efficient membrane capture and suggest novel potential drug targets.

Blood pH Analysis in Combination with Molecular Medical Tools in Relation to COVID-19 Symptoms

The global outbreak of SARS-CoV-2/COVID-19 provided the stage to accumulate an enormous biomedical data set and an opportunity as well as a challenge to test new concepts and strategies to combat the pandemic. New research and molecular medical protocols may be deployed in different scientific fields, e.g., glycobiology, nanopharmacology, or nanomedicine. We correlated clinical biomedical data derived from patients in intensive care units with structural biology and biophysical data from NMR and/or CAMM (computer-aided molecular modeling). Consequently, new diagnostic and therapeutic approaches against SARS-CoV-2 were evaluated. Specifically, we tested the suitability of incretin mimetics with one or two pH-sensitive amino acid residues as potential drugs to prevent or cure long-COVID symptoms. Blood pH values in correlation with temperature alterations in patient bodies were of clinical importance. The effects of biophysical parameters such as temperature and pH value variation in relation to physical-chemical membrane properties (e.g., glycosylation state, affinity of certain amino acid sequences to sialic acids as well as other carbohydrate residues and lipid structures) provided helpful hints in identifying a potential Achilles heel against long COVID. In silico CAMM methods and in vitro NMR experiments (including ³¹P NMR measurements) were applied to analyze the structural behavior of incretin mimetics and SARS-CoV fusion peptides interacting with dodecylphosphocholine (DPC) micelles. These supramolecular complexes were analyzed under physiological conditions by ¹H and ³¹P NMR techniques. We were able to observe characteristic interaction states of incretin mimetics, SARS-CoV fusion peptides and DPC membranes. Novel interaction profiles (indicated, e.g., by ³¹P NMR signal splitting) were detected. Furthermore, we evaluated GM1 gangliosides and sialic acid-coated silica nanoparticles in complex with DPC micelles in order to create a simple virus host cell membrane model. This is a first step in exploring the structure-function relationship between the SARS-CoV-2 spike protein and incretin mimetics with conserved pH-sensitive histidine residues in their carbohydrate recognition domains as found in galectins. The applied methods were effective in identifying peptide sequences as well as certain carbohydrate moieties with the potential to protect the blood-brain barrier (BBB). These clinically relevant observations on low blood pH values in fatal COVID-19 cases open routes for new therapeutic approaches, especially against long-COVID symptoms.

High-throughput identification of prefusion-stabilizing mutations in SARS-CoV-2 spike

Designing prefusion-stabilized SARS-CoV-2 spike is critical for the effectiveness of COVID-19 vaccines. All COVID-19 vaccines in the US encode spike with K986P/V987P mutations to stabilize its prefusion conformation. However, contemporary methods on engineering prefusion-stabilized spike immunogens involve tedious experimental work and heavily rely on structural information. Here, we establish a systematic and unbiased method of identifying mutations that concomitantly improve expression and stabilize the prefusion conformation of the SARS-CoV-2 spike. Our method integrates a fluorescence-based fusion assay, mammalian cell display technology, and deep mutational scanning. As a proof-of-concept, we apply this method to a region in the S2 domain that includes the first heptad repeat and central helix. Our results reveal that besides K986P and V987P, several mutations simultaneously improve expression and significantly lower the fusogenicity of the spike. As prefusion stabilization is a common challenge for viral immunogen design, this work will help accelerate vaccine development against different viruses.

SARS-CoV-2 Protein S Fusion Peptide Is Capable of Wrapping Negatively-Charged Phospholipids

COVID-19, caused by SARS-CoV-2, which is a positive-sense, single-stranded RNA enveloped virus, emerged in late 2019 and was declared a worldwide pandemic in early 2020 causing more than 600 million infections so far and more than 6 million deaths in the world. Although new vaccines have been implemented, the pandemic continues to impact world health dramatically. Membrane fusion, critical for the viral entry into the host cell, is one of the main targets for the development of novel antiviral therapies to combat COVID-19. The S2 subunit of the viral S protein, a class I membrane fusion protein, contains the fusion domain which is directly implicated in the fusion mechanism. The knowledge of the membrane fusion mechanism at the molecular level will undoubtedly result in the development of effective antiviral strategies. We have used all-atom molecular dynamics to analyse the binding of the SARS-CoV-2 fusion peptide to specific phospholipids in model membranes composed of only one phospholipid plus cholesterol in the presence of either Na⁺ or Ca²⁺. Our results show that the fusion peptide is capable of binding to the membrane, that its secondary structure does not change significantly upon binding, that it tends to preferentially bind electronegatively charged phospholipids, and that it does not bind cholesterol at all. Understanding the intricacies of the membrane fusion mechanism and the molecular interactions involved will lead us to the development of antiviral molecules that will allow a more efficient battle against these viruses.

Towards de novo design of transmembrane α-helical assemblies using structural modelling and molecular dynamics simulation

Computational de novo protein design involves iterative processes consisting of amino acid sequence design, structural modelling and scoring, and design validation by synthesis and experimental characterisation. Recent advances in protein structure prediction and modelling methods have enabled the highly efficient and accurate design of water-soluble proteins. However, the design of membrane proteins remains a major challenge. To advance membrane protein design, considering the higher complexity of membrane protein folding, stability, and dynamic interactions between water, ions, lipids, and proteins is an important task. For introducing explicit solvents and membranes to these design methods, all-atom molecular dynamics (MD) simulations of designed proteins provide useful information that cannot be obtained experimentally. In this review, we first describe two major approaches to designing transmembrane α-helical assemblies, consensus and de novo design. We further illustrate recent MD studies of membrane protein folding related to protein design, as well as advanced treatments in molecular models and conformational sampling techniques in the simulations. Finally, we discuss the possibility to introduce MD simulations after the existing static modelling and screening of design decoys as an additional step for refinement of the design, which considers membrane protein folding dynamics and interactions with explicit membranes.

The Spread of SARS-CoV-2 Omicron Variant in CALABRIA: A Spatio-Temporal Report of Viral Genome Evolution

We investigated the evolution of SARS-CoV-2 spread in Calabria, Southern Italy, in 2022. A total of 272 RNA isolates from nasopharyngeal swabs of individuals infected with SARS-CoV-2 were sequenced by whole genome sequencing (N = 172) and/or Sanger sequencing (N = 100). Analysis of diffusion of Omicron variants in Calabria revealed the prevalence of 10 different sub-lineages (recombinant BA.1/BA.2, BA.1, BA.1.1, BA.2, BA.2.9, BA.2.10, BA.2.12.1, BA.4, BA.5, BE.1). We observed that Omicron spread in Calabria presented a similar trend as in Italy, with some notable exceptions: BA.1 disappeared in April in Calabria but not in the rest of Italy; recombinant BA.1/BA.2 showed higher frequency in Calabria (13%) than in the rest of Italy (0.02%); BA.2.9, BA.4 and BA.5 emerged in Calabria later than in other Italian regions. In addition, Calabria Omicron presented 16 non-canonical mutations in the S protein and 151 non-canonical mutations in non-structural proteins. Most non-canonical mutations in the S protein occurred mainly in BA.5 whereas non-canonical mutations in non-structural or accessory proteins (ORF1ab, ORF3a, ORF8 and N) were identified in BA.2 and BA.5 sub-lineages. In conclusion, the data reported here underscore the importance of monitoring the entire SARS-CoV-2 genome.

Effect of the Graphene Nanosheet on Functions of the Spike Protein in Open and Closed States: Comparison between SARS-CoV-2 Wild Type and the Omicron Variant

The spread of coronavirus disease 2019 caused by SARS-CoV-2 and its variants has become a global health crisis. Although there were many attempts to use nanomaterials-based devices to fight against SARS-CoV-2, it still remains elusive as to how the nanomaterials interact with SARS-CoV-2 and affect its biofunctions. Here, taking the graphene nanosheet (GN) as the model nanomaterial, we investigate its interaction with the spike protein in both WT and Omicron by molecular simulations. In the closed state, the GN can insert into the region between the receptor binding domain (RBD) and the N-terminal domain (NTD) in both wild type (WT) and Omicron, which keeps the RBD in the down conformation. In the open state, the GN can hamper the binding of up RBD to ACE2 in WT, but it has little impact on up RBD and, even worse, stimulates the down-to-up transition of down RBDs in Omicron. Moreover, the GN can insert in the vicinity of the fusion peptide in both WT and Omicron and prevents the detachment of S1 from the whole spike protein. The present study reveals the effect of the SARS-CoV-2 variant on the nanomaterial-spike protein interaction, which informs prospective efforts to design functional nanomaterials against SARS-CoV-2.

Effect of Disulfide Bridge on the Binding of SARS-CoV-2 Fusion Peptide to Cell Membrane: A Coarse-Grained Study

In this paper, we present the parameterization of the CAVS coarse-grained (CG) force field for 20 amino acids, and our CG simulations show that the CAVS force field could accurately predict the amino acid tendency of the secondary structure. Then, we used the CAVS force field to investigate the binding of a severe acute respiratory syndrome-associated coronavirus fusion peptide (SARS-CoV-2 FP) to a phospholipid bilayer: a long FP (FP-L) containing 40 amino acids and a short FP (FP-S) containing 26 amino acids. Our CAVS CG simulations displayed that the binding affinity of the FP-L to the bilayer is higher than that of the FP-S. We found that the FP-L interacted more strongly with membrane cholesterol than the FP-S, which should be attributed to the stable helical structure of the FP-L at the C-terminus. In addition, we discovered that the FP-S had one major and two minor membrane-bound states, in agreement with previous all-atom molecular dynamics (MD) studies. However, we found that both the C-terminal and N-terminal amino acid residues of the FP-L can strongly interact with the bilayer membrane. Furthermore, we found that the disulfide bond formed between Cys840 and Cys851 stabilized the helices of the FP-L at the C-terminus, enhancing the interaction between the FP-L and the bilayer membrane. Our work indicates that the stable helical structure is crucial for binding the SARS-CoV-2 FP to cell membranes. In particular, the helical stability of FP should have a significant influence on the FP-membrane binding.

High-throughput identification of prefusion-stabilizing mutations in SARS-CoV-2 spike

Designing prefusion-stabilized SARS-CoV-2 spike is critical for the effectiveness of COVID-19 vaccines. All COVID-19 vaccines in the US encode spike with K986P/V987P mutations to stabilize its prefusion conformation. However, contemporary methods on engineering prefusion-stabilized spike immunogens involve tedious experimental work and heavily rely on structural information. Here, we established a systematic and unbiased method of identifying mutations that concomitantly improve expression and stabilize the prefusion conformation of the SARS-CoV-2 spike. Our method integrated a fluorescence-based fusion assay, mammalian cell display technology, and deep mutational scanning. As a proof-of-concept, this method was applied to a region in the S2 domain that includes the first heptad repeat and central helix. Our results revealed that besides K986P and V987P, several mutations simultaneously improved expression and significantly lowered the fusogenicity of the spike. As prefusion stabilization is a common challenge for viral immunogen design, this work will help accelerate vaccine development against different viruses.

S2 Subunit of SARS-CoV-2 Spike Protein Induces Domain Fusion in Natural Pulmonary Surfactant Monolayers

Pulmonary surfactant has been attempted as a supportive therapy to treat COVID-19. Although it is mechanistically accepted that the fusion peptide in the S2 subunit of the S protein plays a predominant role in mediating viral fusion with the host cell membrane, it is still unknown how the S2 subunit interacts with the natural surfactant film. Using combined bio-physicochemical assays and atomic force microscopy imaging, it was found that the S2 subunit inhibited the biophysical properties of the surfactant and induced microdomain fusion in the surfactant monolayer. The surfactant inhibition has been attributed to membrane fluidization caused by insertion of the S2 subunit mediated by its fusion peptide. These findings may provide novel insight into the understanding of bio-physicochemical mechanisms responsible for surfactant interactions with SARS-CoV-2 and may have translational implications in the further development of surfactant replacement therapy for COVID-19 patients.

Exploring the pH dependence of the SARS‐CoV ‐2 complete fusion domain and the role of its unique structural features

SARS‐CoV‐2 may enter target cells through the process of membrane fusion at either the plasma (~pH 7.4–7.0) or endosomal (~pH 6.5–5.0) membrane in order to deliver its genetic information. The fusion domain (FD) of the spike glycoprotein is responsible for initiating fusion and is thus integral to the viral life cycle. The FD of SARS‐CoV‐2 is unique in that it consists of two structurally distinctive regions referred to as the fusion peptide (FP) and the fusion loop (FL); yet the molecular mechanisms behind how this FD perturbs the membrane to initiate fusion remains unclear. In this study via solution NMR, we witnessed only a slight conformational change in the FD between pH 7.4 and pH 5.0, resulting in a minor elongation of helix 1. However, we found that the FD's ability to mediate membrane fusion has a large and significant pH dependence, with fusion events being more readily induced at low pH. Interestingly, a biphasic relationship between the environmental pH and fusogenicity was discovered, suggesting a preference for the FD to initiate fusion at the late endosomal membrane. Furthermore, the conserved disulfide bond and hydrophobic motif “LLF” were found to be critical for the function of the complete FD, with minimal activity witnessed when either was perturbed. In conclusion, these findings indicate that the SARS‐CoV‐2 FD preferably initiates fusion at a pH similar to the late endosome through a mechanism that heavily relies on the internal disulfide bond of the FL and hydrophobic LLF motif within the FP.

Different Binding Modes of SARS-CoV-1 and SARS-CoV-2 Fusion Peptides to Cell Membranes: The Influence of Peptide Helix Length

Although the amino acid sequences of SARS-CoV-1 and SARS-CoV-2 fusion peptides (FPs) are highly conserved, the cryo-electron microscopy structures of the SARS-CoV-1 and SARS-CoV-2 spike proteins show that the helix length of SARS-CoV-1 FP is longer than that of SARS-CoV-2 FP. In this work, we simulated the membrane-binding models of SARS-CoV-1 and SARS-CoV-2 FPs and compared the binding modes of the FPs with the POPC/POPE/cholesterol bilayer membrane. Our simulation results show that the SARS-CoV-2 FP binds to the bilayer membrane more effectively than the SARS-CoV-1 FP. It is seen that the short N-terminal helix of SARS-CoV-2 FP is more favorable to insert into the target membrane than the long N-terminal helix of SARS-CoV-1 FP. Meanwhile, the potential of mean force calculations showed that the SARS-CoV-2 FP would prefer only one binding mode (N-terminal binding), whereas the SARS-CoV-1 FP has two favorable membrane-binding modes (C-terminal and N-terminal binding modes). Moreover, in the case of SARS-CoV-1 FP binding to the target membrane, the transition between the two binding modes is relatively fast. Finally, we discovered that the membrane-binding mode would influence the helix length of SARS-CoV-1 FP, while the helix length of SARS-CoV-2 FP could be stably maintained in the membrane-bound configurations. This work suggests that the short helix might endow the FP with high membrane-anchoring strength. In particular, the membrane-penetrating residues (Phe, Ile, and Leu) of short α-helix interact with the cell membrane more strongly than those of long α-helix.

COVID-19 and Neurodegenerative Diseases: Prion-Like Spread and Long-Term Consequences

COVID-19 emerged as a global pandemic starting from Wuhan in China and spread at a lightning speed to the rest of the world. One of the potential long-term outcomes that we speculate is the development of neurodegenerative diseases as a long-term consequence of SARS-CoV-2 especially in people that have developed severe neurological symptoms. Severe inflammatory reactions and aging are two very strong common links between neurodegenerative diseases and COVID-19. Thus, patients that have very high viral load may be at high risk of developing long-term adverse neurological consequences such as dementia. We hypothesize that people with neurodegenerative diseases such as Alzheimer’s disease, Parkinson’s disease, and aged people are at higher risk of getting the COVID-19 than normal adults. The basis of this hypothesis is the fact that SARS-CoV-2 uses as a receptor angiotensin-converting enzyme 2 to enter the host cell and that this interaction is calcium-dependent. This could then suggest a direct relationship between neurodegenerative diseases, ACE-2 expression, and the susceptibility to COVID-19. The analysis of the available literature showed that COVID-19 virus is neurotropic and was found in the brains of patients infected with this virus. Furthermore, that the risk of having the infection increases with dementia and that infected people with severe symptoms could develop dementia as a long-term consequence. Dementia could be developed following the acceleration of the spread of prion-like proteins. In the present review we discuss current reports concerning the prevalence of COVID-19 in dementia patients, the individuals that are at high risk of suffering from dementia and the potential acceleration of prion-like proteins spread following SARS-CoV-2 infection.

Computational methods to study enveloped viral entry

The protein-membrane interactions that mediate viral infection occur via loosely ordered, transient assemblies, creating challenges for high-resolution structure determination. Computational methods and in particular molecular dynamics simulation have thus become important adjuncts for integrating experimental data, developing mechanistic models, and suggesting testable hypotheses regarding viral function. However, the large molecular scales of virus-host interaction also create challenges for detailed molecular simulation. For this reason, continuum membrane models have played a large historical role, although they have become less favored for high-resolution models of protein assemblies and lipid organization. Here, we review recent progress in the field, with an emphasis on the insight that has been gained using a mixture of coarse-grained and atomic-resolution molecular dynamics simulations. Based on successes and challenges to date, we suggest a multiresolution strategy that should yield the best mixture of computational efficiency and physical fidelity. This strategy may facilitate further simulations of viral entry by a broader range of viruses, helping illuminate the diversity of viral entry strategies and the essential common elements that can be targeted for antiviral therapies.

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

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

Identifying Distinct Structural Features of the SARS-CoV-2 Spike Protein Fusion Domain Essential for Membrane Interaction

The SARS-CoV-2 spike protein is the primary antigenic determinant of the virus and has been studied extensively, yet the process of membrane fusion remains poorly understood. The fusion domain (FD) of viral glycoproteins is well established as facilitating the initiation of membrane fusion. An improved understanding of the structural plasticity associated with these highly conserved regions aids in our knowledge of the molecular mechanisms that drive viral fusion. Within the spike protein, the FD of SARS-CoV-2 exists immediately following S2' cleavage at the N-terminus of the S2 domain. Here we have shown that following the introduction of a membrane at pH 7.4, the FD undergoes a transition from a random coil to a more structurally well-defined postfusion state. Furthermore, we have classified the domain into two distinct regions, a fusion peptide (FP, S816-G838) and a fusion loop (FL, D839-F855). The FP forms a helix-turn-helix motif upon association with a membrane, and the favorable entropy gained during this transition from a random coil is likely the driving force behind membrane insertion. Membrane depth experiments then revealed the FP is found inserted within the membrane below the lipid headgroups, while the interaction of the FL with the membrane is shallower in nature. Thus, we propose a structural model relevant to fusion at the plasma membrane in which the FP inserts itself just below the phospholipid headgroups and the FL lays upon the lipid membrane surface.

Sterically confined rearrangements of SARS-CoV-2 Spike protein control cell invasion

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is highly contagious, and transmission involves a series of processes that may be targeted by vaccines and therapeutics. During transmission, host cell invasion is controlled by a large-scale (200–300 Å) conformational change of the Spike protein. This conformational rearrangement leads to membrane fusion, which creates transmembrane pores through which the viral genome is passed to the host. During Spike-protein-mediated fusion, the fusion peptides must be released from the core of the protein and associate with the host membrane. While infection relies on this transition between the prefusion and postfusion conformations, there has yet to be a biophysical characterization reported for this rearrangement. That is, structures are available for the endpoints, though the intermediate conformational processes have not been described. Interestingly, the Spike protein possesses many post-translational modifications, in the form of branched glycans that flank the surface of the assembly. With the current lack of data on the pre-to-post transition, the precise role of glycans during cell invasion has also remained unclear. To provide an initial mechanistic description of the pre-to-post rearrangement, an all-atom model with simplified energetics was used to perform thousands of simulations in which the protein transitions between the prefusion and postfusion conformations. These simulations indicate that the steric composition of the glycans can induce a pause during the Spike protein conformational change. We additionally show that this glycan-induced delay provides a critical opportunity for the fusion peptides to capture the host cell. In contrast, in the absence of glycans, the viral particle would likely fail to enter the host. This analysis reveals how the glycosylation state can regulate infectivity, while providing a much-needed structural framework for studying the dynamics of this pervasive pathogen.

Binding of SARS-CoV-2 Fusion Peptide to Host Endosome and Plasma Membrane

During infection the SARS-CoV-2 virus fuses its viral envelope with cellular membranes of its human host. The viral spike (S) protein mediates both the initial contact with the host cell and the subsequent membrane fusion. Proteolytic cleavage of S at the S2′ site exposes its fusion peptide (FP) as the new N-terminus. By binding to the host membrane, the FP anchors the virus to the host cell. The reorganization of S2 between virus and host then pulls the two membranes together. Here we use molecular dynamics (MD) simulations to study the two core functions of the SARS-CoV-2 FP: to attach quickly to cellular membranes and to form an anchor strong enough to withstand the mechanical force during membrane fusion. In eight 10 μs long MD simulations of FP in proximity to endosomal and plasma membranes, we find that FP binds spontaneously to the membranes and that binding proceeds predominantly by insertion of two short amphipathic helices into the membrane interface. Connected via a flexible linker, the two helices can bind the membrane independently, yet binding of one promotes the binding of the other by tethering it close to the target membrane. By simulating mechanical pulling forces acting on the C-terminus of the FP, we then show that the bound FP can bear forces up to 250 pN before detaching from the membrane. This detachment force is more than 10-fold higher than an estimate of the force required to pull host and viral membranes together for fusion. We identify a fully conserved disulfide bridge in the FP as a major factor for the high mechanical stability of the FP membrane anchor. We conclude, first, that the sequential binding of two short amphipathic helices allows the SARS-CoV-2 FP to insert quickly into the target membrane, before the virion is swept away after shedding the S1 domain connecting it to the host cell receptor. Second, we conclude that the double attachment and the conserved disulfide bridge establish the strong anchoring required for subsequent membrane fusion. Multiple distinct membrane-anchoring elements ensure high avidity and high mechanical strength of FP–membrane binding.

Ca2+-dependent mechanism of membrane insertion and destabilization by the SARS-CoV-2 fusion peptide

Cell penetration after recognition of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus by the ACE2 receptor and the fusion of its viral envelope membrane with cellular membranes are the early steps of infectivity. A region of the Spike protein of the virus, identified as the "fusion peptide" (FP), is liberated at its N-terminal site by a specific cleavage occurring in concert with the interaction of the receptor-binding domain of the Spike. Studies have shown that penetration is enhanced by the required binding of Ca²⁺ ions to the FPs of coronaviruses, but the mechanisms of membrane insertion and destabilization remain unclear. We have predicted the preferred positions of Ca²⁺ binding to the SARS-CoV-2-FP, the role of Ca²⁺ ions in mediating peptide-membrane interactions, the preferred mode of insertion of the Ca²⁺-bound SARS-CoV-2-FP, and consequent effects on the lipid bilayer from extensive atomistic molecular dynamics simulations and trajectory analyses. In a systematic sampling of the interactions of the Ca²⁺-bound peptide models with lipid membranes, SARS-CoV-2-FP penetrated the bilayer and disrupted its organization only in two modes involving different structural domains. In one, the hydrophobic residues F833/I834 from the middle region of the peptide are inserted. In the other, more prevalent mode, the penetration involves residues L822/F823 from the LLF motif, which is conserved in CoV-2-like viruses, and is achieved by the binding of Ca²⁺ ions to the D830/D839 and E819/D820 residue pairs. FP penetration is shown to modify the molecular organization in specific areas of the bilayer, and the extent of membrane binding of the SARS-CoV-2 FP is significantly reduced in the absence of Ca²⁺ ions. These findings provide novel mechanistic insights regarding the role of Ca²⁺ in mediating SARS-CoV-2 fusion and provide a detailed structural platform to aid the ongoing efforts in rational design of compounds to inhibit SARS-CoV-2 cell entry.