Distant Residues Modulate Conformational Opening in SARS-CoV-2 Spike Protein

Infection by SARS-CoV-2 involves the attachment of the receptor binding domain (RBD) of its spike proteins to the ACE2 receptors on the peripheral membrane of host cells. Binding is initiated by a down-to-up conformational change in the spike protein, the change that presents the RBD to the receptor. To date, computational and experimental studies that search for therapeutics have concentrated, for good reason, on the RBD. However, the RBD region is highly prone to mutations, and is therefore a hotspot for drug resistance. In contrast, we here focus on the correlations between the RBD and residues distant to it in the spike protein. This allows for a deeper understanding of the underlying molecular recognition events and prediction of the highest-effect key mutations in distant, allosteric sites, with implications for therapeutics. Also, these sites can appear in emerging mutants with possibly higher transmissibility and virulence, and pre-identifying them can give clues for designing pancoronavirus vaccines against future outbreaks. Our model, based on time-lagged independent component analysis (tICA) and protein graph connectivity network, is able to identify multiple residues that exhibit long-distance coupling with the RBD opening. Residues involved in the most ubiquitous D614G mutation and the A570D mutation of the highly contagious UK SARS-CoV-2 variant are predicted ab-initio from our model. Conversely, broad spectrum therapeutics like drugs and monoclonal antibodies can target these key distant-but-conserved regions of the spike protein. Significance StatementThe novel coronavirus (SARS-CoV-2) pandemic resulted in the largest public health crisis in recent times. Significant drug design effort against SARS-CoV-2 is focused on the receptor binding domain (RBD) of the spike protein, although this region is highly prone to mutations causing therapeutic resistance. We applied deep data analysis methods on all-atom molecular dynamics simulations to identify key non-RBD residues that play a crucial role in spike-receptor binding and infection. textcol-orredBecause the non-RBD residues are typically conserved across multiple coronaviruses, they can be targeted by broad spectrum antibodies and drugs to treat infections from new strains that might appear during future epidemics.