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

Membrane fusion, a key step in the early stages of virus propagation, allows the release of the viral genome in the host cell cytoplasm. The process is initiated by fusion peptides that are small, hydrophobic components of viral membrane-embedded glycoproteins and are typically conserved within virus families. Here, we attempted to identify the correct fusion peptide region in the Spike protein of SARS-CoV-2 by all-atom molecular dynamics simulations of dual membrane systems with varied oligomeric units of putative candidate peptides. Of all of the systems tested, only a trimeric unit of a 40-amino-acid region (residues 816-855 of SARS-CoV-2 Spike) was effective in triggering the initial stages of membrane fusion, within 200 ns of simulation time. Association of this trimeric unit with dual membranes resulted in the migration of lipids from the upper leaflet of the lower bilayer toward the lower leaflet of the upper bilayer to create a structural unit reminiscent of a fusion bridge. We submit that residues 816-855 of Spike represent the bona fide fusion peptide of SARS-CoV-2 and that computational methods represent an effective way to identify fusion peptides in viral glycoproteins.