Computational Insight Into the Mechanism of SARS-CoV-2 Membrane Fusion

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.

New Strategies for Responding to SARS-CoV-2: The Present and Future of Dual-Target Drugs

The 2019 coronavirus disease (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in millions of deaths, posing a serious threat to public health and safety. Rapid mutations of SARS-CoV-2 and complex interactions among multiple targets during infection pose a risk of expiry for small molecule inhibitors. This suggests that the traditional concept of "one bug, one drug" could be ineffective in dealing with the coronavirus. The dual-target drug strategy is expected to be the key to ending coronavirus infections. However, the lack of design method and improper combination of dual-targets poses obstacle to the discovery of new dual-target drugs. In this Perspective, we summarized the profiles concerning drug design methods, structure-activity relationships, and pharmacological parameters of dual-target drugs for the treatment of COVID-19. Importantly, we underscored how target combination and rational drug design illuminate the development of dual-target drugs for SARS-CoV-2.

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.

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.

SARS-CoV-2 pseudovirus enters the host cells through spike protein-CD147 in an Arf6-dependent manner

The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its variants is threatening public health around the world. Endocytosis functions as an important way for viral infection, and SARS-CoV-2 bears no exception. However, the specific endocytic mechanism of SARS-CoV-2 remains unknown. In this study, we used endocytic inhibitors to evaluate the role of different endocytic routes in SARS-CoV-2 pseudovirus infection and found that the viral infection was associated with caveolar/lipid raft- and cytoskeleton-mediated endocytosis, but independent of the clathrin-mediated endocytosis and macropinocytosis. Meanwhile, the knockdown of CD147 and Rab5a in Vero E6 and Huh-7 cells inhibited SARS-CoV-2 pseudovirus infection, and the co-localization of spike protein, CD147, and Rab5a was observed in pseudovirus-infected Vero E6 cells, which was weakened by CD147 silencing, illustrating that SARS-CoV-2 pseudovirus entered the host cells via CD147-mediated endocytosis. Additionally, Arf6 silencing markedly inhibited pseudovirus infection in Vero E6 and Huh-7 cells, while little change was observed in CD147 knockout-Vero E6 cells. This finding indicated Arf6-mediated CD147 trafficking plays a vital role in SARS-CoV-2 entry. Taken together, our findings provide new insights into the CD147-Arf6 axis in mediating SARS-CoV-2 pseudovirus entry into the host cells, and further suggest that blockade of this pathway seems to be a feasible approach to prevent the SARS-CoV-2 infection clinically.

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.

Estimation of virus-mediated cell fusion rate of SARS-CoV-2

Several viruses have the ability to form large multinucleated cells known as syncytia. Many properties of syncytia and the role they play in the evolution of a viral infection are not well understood. One basic question that has not yet been answered is how quickly syncytia form. We use a novel mathematical model of cell-cell fusion assays and apply it to experimental data from SARS-CoV-2 fusion assays to provide the first estimates of virus-mediated cell fusion rate. We find that for SARS-CoV2, the fusion rate is in the range of 6 × 10-4-12×10-4/h. We also use our model to compare fusion rates when the protease TMPRSS2 is overexpressed (2-4 times larger fusion rate), when the protease furin is removed (one third the original fusion rate), and when the spike protein is altered (1/10th the original fusion rate). The use of mathematical models allows us to provide additional quantitative information about syncytia formation.

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.

Broad-Spectrum Antiviral Activity of the Amphibian Antimicrobial Peptide Temporin L and Its Analogs

The COVID-19 pandemic has evidenced the urgent need for the discovery of broad-spectrum antiviral therapies that could be deployed in the case of future emergence of novel viral threats, as well as to back up current therapeutic options in the case of drug resistance development. Most current antivirals are directed to inhibit specific viruses since these therapeutic molecules are designed to act on a specific viral target with the objective of interfering with a precise step in the replication cycle. Therefore, antimicrobial peptides (AMPs) have been identified as promising antiviral agents that could help to overcome this limitation and provide compounds able to act on more than a single viral family. We evaluated the antiviral activity of an amphibian peptide known for its strong antimicrobial activity against both Gram-positive and Gram-negative bacteria, namely Temporin L (TL). Previous studies have revealed that TL is endowed with widespread antimicrobial activity and possesses marked haemolytic activity. Therefore, we analyzed TL and a previously identified TL derivative (Pro³, DLeu⁹ TL, where glutamine at position 3 is replaced with proline, and the D-Leucine enantiomer is present at position 9) as well as its analogs, for their activity against a wide panel of viruses comprising enveloped, naked, DNA and RNA viruses. We report significant inhibition activity against herpesviruses, paramyxoviruses, influenza virus and coronaviruses, including SARS-CoV-2. Moreover, we further modified our best candidate by lipidation and demonstrated a highly reduced cytotoxicity with improved antiviral effect. Our results show a potent and selective antiviral activity of TL peptides, indicating that the novel lipidated temporin-based antiviral agents could prove to be useful additions to current drugs in combatting rising drug resistance and epidemic/pandemic emergencies.

Accelerating the discovery of the beyond rule of five compounds that have high affinities toward SARS-CoV-2 spike RBD

The battle against SARS-CoV-2 coronavirus is the focal point for the global pandemic that has affected millions of lives worldwide. The need for effective and selective therapeutics for the treatment of the disease caused by SARS-CoV-2 is critical. Herein, we performed a hierarchical computational approach incorporating molecular docking studies, molecular dynamics simulations, absolute binding energy calculations, and steered molecular dynamics simulations for the discovery of potential compounds with high affinity towards SARS-CoV-2 spike RBD. By leveraging ZINC15 database, a total of 1282 in-clinical and FDA approved drugs were filtered out from nearly 0.5 million protomers of relatively large compounds (MW > 500, and LogP ≤ 5). Our results depict plausible mechanistic aspects related to the blockage of SARS-CoV-2 spike RBD by the top hits discovered. We found that the most promising candidates, namely, ZINC95628821, ZINC95617623, ZINC3979524, and ZINC261494658, strongly bind to the spike RBD and interfere with the human ACE2 receptor. These findings accelerate the rational design of selective inhibitors targeting the spike RBD protein of SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Repurposing fusion inhibitor peptide against SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is continuously evolving. Although several vaccines were approved, this pandemic is still a major threat to public life. Till date, no established therapies are available against SARS-CoV-2. Peptide inhibitors hold great promise for this viral pathogen due to their efficacy, safety, and specificity. In this study, seventeen antiviral peptides which were known to inhibit SARS-CoV-1 are collected and computationally screened against heptad repeat 1 (HR1) of the SARS-CoV-2 spike protein (S2). Out of 17 peptides, Fp13 and Fp14 showed better binding affinity toward HR1 compared to a control peptide EK1 (a modified pan-coronavirus fusion inhibitor) in molecular docking. To explore the time-dependent interactions of the fusion peptide with HR1, molecular dynamics simulation was performed incorporating lipid membrane. During 100 ns MD simulation, structural and energy parameters of Fp13-HR1 and Fp14-HR1 complexes demonstrated lower fluctuations compared to the control EK1-HR1 complex. Furthermore, principal component analysis and free energy landscape study revealed that these two peptides (Fp13 and Fp14) strongly bind to the HR1 with higher affinity than that of control EK1. Tyr917, Asn919, Gln926, lys933, and Gln949 residues in HR1 protein were found to be crucial residues for peptide interaction. Notably, Fp13, Fp14 showed reasonably better binding free energy and hydrogen bond contribution than that of EK1. Taken together, Fp13 and Fp14 peptides may be highly specific for HR1 which can potentially prevent the formation of the fusion core and could be further developed as therapeutics for treatment or prophylaxis of SARS-CoV-2 infection.

COVID-19: the CaMKII-like system of S protein drives membrane fusion and induces syncytial multinucleated giant cells

The SARS-CoV-2 S protein on the membrane of infected cells can promote receptor-dependent syncytia formation, relating to extensive tissue damage and lymphocyte elimination. In this case, it is challenging to obtain neutralizing antibodies and prevent them through antibodies effectively. Considering that, in the current study, structural domain search methods are adopted to analyze the SARS-CoV-2 S protein to find the fusion mechanism. The results show that after the EF-hand domain of S protein bound to calcium ions, S2 protein had CaMKII protein activities. Besides, the CaMKII_AD domain of S2 changed S2 conformation, facilitating the formation of HR1-HR2 six-helix bundles. Apart from that, the Ca2+-ATPase of S2 pumped calcium ions from the virus cytoplasm to help membrane fusion, while motor structures of S drove the CaATP_NAI and CaMKII_AD domains to extend to the outside and combined the viral membrane and the cell membrane, thus forming a calcium bridge. Furthermore, the phospholipid-flipping-ATPase released water, triggering lipid mixing and fusion and generating fusion pores. Then, motor structures promoted fusion pore extension, followed by the cytoplasmic contents of the virus being discharged into the cell cytoplasm. After that, the membrane of the virus slid onto the cell membrane along the flowing membrane on the gap of the three CaATP_NAI. At last, the HR1-HR2 hexamer would fall into the cytoplasm or stay on the cell membrane. Therefore, the CaMKII_like system of S protein facilitated membrane fusion for further inducing syncytial multinucleated giant cells.

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.

Computational insights into the membrane fusion mechanism of SARS-CoV-2 at the cellular level

The membrane fusion mechanism of SARS-CoV-2 offers an attractive target for the development of small molecule antiviral inhibitors. Fusion involves an initial binding of the crown-like trimeric spike glycoproteins of SARS-CoV-2 to the receptor angiotensin II-converting enzyme 2 (ACE2) on the permissive host cellular membrane and a prefusion to post-fusion conversion of the spike trimer. During this conversion, the fusion peptides of the spike trimer are inserted into the host membrane to bring together the host and viral membranes for membrane fusion in highly choreographic events. However, geometric constraints due to interactions with the membranes remain poorly understood. In this study, we build structural models of super-complexes of spike trimer/ACE2 dimers based on the molecular structures of the ACE2/neutral amino acid transporter B(0)AT heterodimer. We determine the conformational constraints due to the membrane geometry on the enzymatic activity of ACE2 and on the viral fusion process. Furthermore, we find that binding three ACE2 dimers per spike trimer is essential to open the central pore as necessary for triggering productive membrane fusion through an elongation of the central stalk. The reported findings thus provide valuable insights for targeting the membrane fusion mechanism for drug design at the molecular level.

Accelerating COVID-19 Research Using Molecular Dynamics Simulation

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