Abstract:The SARS-CoV-1 spike glycoprotein contains a fusion peptide (FP) segment that mediates the fusion of the viral and host cell membranes. Calcium ions are thought to position the FP optimally for membrane insertion by interacting with negatively charged residues in this segment (E801, D802, D812, E821, D825, and D830); however, which residues bind to calcium and in what combinations supportive of membrane insertion are unknown. Using biological assays and molecular dynamics studies, we have determined the functional configurations of FP-Ca2+ binding that likely promote membrane insertion. We first individually mutated the negatively charged residues in the SARS CoV-1 FP to assay their roles in cell entry and syncytia formation, finding that charge loss in the D802A or D830A mutants greatly reduced syncytia formation and pseudoparticle transduction of VeroE6 cells. Interestingly, one mutation (D812A) led to a modest increase in cell transduction, further indicating that FP function likely depends on calcium binding at specific residues and in specific combinations. To interpret these results mechanistically and identify specific modes of FP-Ca2+ binding that modulate membrane insertion, we performed molecular dynamics simulations of the SARS-CoV-1 FP and Ca2+ions. The preferred residue pairs for Ca2+ binding we identified (E801/D802, E801/D830, and D812/E821) include the two residues found to be essential for S function in our biological studies (D802 and D830). The three preferred Ca2+ binding pairs were also predicted to promote FP membrane insertion. We also identified a Ca2+ binding pair (E821/D825) predicted to inhibit FP membrane insertion. We then carried out simulations in the presence of membranes and found that binding of Ca2+ to SARS-CoV-1 FP residue pairs E801/D802 and D812/E821 facilitates membrane insertion by enabling the peptide to adopt conformations that shield the negative charges of the FP to reduce repulsion by the membrane phospholipid headgroups. This calcium binding mode also optimally positions the hydrophobic LLF region of the FP for membrane penetration. Conversely, Ca2+ binding to the FP E801/D802 and D821/D825 pairs eliminates the negative charge screening and instead creates a repulsive negative charge that hinders membrane penetration of the LLF motif. These computational results, taken together with our biological studies, provide an improved and nuanced mechanistic understanding of the dymanics of SARS-CoV-1 calcium binding and their potential effects on host cell entry.

Abstract P-14: Molecular Modeling of the Transmembrane Domain of the SARS Cov-2 S-Protein and its Interaction with the Membrane

Cryo-electron Microscopy Structures of the SARS-CoV Spike Glycoprotein Reveal a Prerequisite Conformational State for Receptor Binding

Cryo-EM Structure of the SARS Coronavirus Spike Glycoprotein in Complex with Its Host Cell Receptor ACE2

Characterization of the Heptad Repeat Regions, HR1 and HR2, and Design of a Fusion Core Structure Model of the Spike Protein from Severe Acute Respiratory Syndrome (SARS) Coronavirus.

Sequence of the SARS-CoV-2 spike transmembrane domain makes it inherently dynamic

Structure, Dynamics, Receptor Binding, and Antibody Binding of the Fully Glycosylated Full-Length SARS-CoV-2 Spike Protein in a Viral Membrane

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

Developing a Fully Glycosylated Full-Length SARS-CoV-2 Spike Protein Model in a Viral Membrane

Folding of the SARS Coronavirus Spike Glycoprotein Immunological Fragment (sars_s1b): Thermodynamic and Kinetic Investigation Correlating with Three-Dimensional Structural Modeling.

Structures and distributions of SARS-CoV-2 spike proteins on intact virions

The sequence and structural integrity of the SARS-CoV-2 Spike protein transmembrane domain is crucial for viral entry

Silico analysis of interaction between full-length SARS-CoV2 S protein with human Ace2 receptor: Modelling, docking, MD simulation

Distinct conformational states of SARS-CoV-2 spike protein

Molecular Dynamics of DHHC20 Acyltransferase Suggests Principles of Lipid and Protein Substrate Selectivity

Exploring the Regulatory Function of the N-terminal Domain of SARS-CoV-2 Spike Protein Through Molecular Dynamics Simulation

Deciphering the protein motion of S1 subunit in SARS-CoV-2 spike glycoprotein through integrated computational methods

Robotic-inspired approach to multi-domain membrane receptor conformation space: theory and SARS-CoV-2 spike protein case study

Modelling of SARS-CoV-2 spike protein structures at varying pH values

A Mechanistic Understanding of the Modes of Ca 2+ Ion Binding to the SARS-CoV-1 Fusion Peptide and Their Role in the Dynamics of Host Membrane Penetration

Exploring the Regulatory Function of the N ‐terminal Domain of SARS‐CoV‐2 Spike Protein Through Molecular Dynamics Simulation (adv. Theory Simul. 10/2021)

C-Terminal Domain of Sars-Cov Main Protease Can Form A 3d Domain-Swapped Dimer