Insights on the Structural Variations of the Furin-Like Cleavage Site Found Among the December 2019-July 2020 SARS-CoV-2 Spike Glycoprotein: A Computational Study Linking Viral Evolution and Infection

Capturing a Crucial ‘Disorder-to-Order Transition’ at the Heart of the Coronavirus Molecular Pathology—Triggered by Highly Persistent, Interchangeable Salt-Bridges

The COVID-19 origin debate has greatly been influenced by genome comparison studies of late, revealing the emergence of the Furin-like cleavage site at the S1/S2 junction of the SARS-CoV-2 Spike (FLCSSpike) containing its 681PRRAR685 motif, absent in other related respiratory viruses. Being the rate-limiting (i.e., the slowest) step, the host Furin cleavage is instrumental in the abrupt increase in transmissibility in COVID-19, compared to earlier onsets of respiratory viral diseases. In such a context, the current paper entraps a ‘disorder-to-order transition’ of the FLCSSpike (concomitant to an entropy arrest) upon binding to Furin. The interaction clearly seems to be optimized for a more efficient proteolytic cleavage in SARS-CoV-2. The study further shows the formation of dynamically interchangeable and persistent networks of salt-bridges at the Spike–Furin interface in SARS-CoV-2 involving the three arginines (R682, R683, R685) of the FLCSSpike with several anionic residues (E230, E236, D259, D264, D306) coming from Furin, strategically distributed around its catalytic triad. Multiplicity and structural degeneracy of plausible salt-bridge network archetypes seem to be the other key characteristic features of the Spike–Furin binding in SARS-CoV-2, allowing the system to breathe—a trademark of protein disorder transitions. Interestingly, with respect to the homologous interaction in SARS-CoV (2002/2003) taken as a baseline, the Spike–Furin binding events, generally, in the coronavirus lineage, seems to have preference for ionic bond formation, even with a lesser number of cationic residues at their potentially polybasic FLCSSpike patches. The interaction energies are suggestive of characteristic metastabilities attributed to Spike–Furin interactions, generally to the coronavirus lineage, which appears to be favorable for proteolytic cleavages targeted at flexible protein loops. The current findings not only offer novel mechanistic insights into the coronavirus molecular pathology and evolution, but also add substantially to the existing theories of proteolytic cleavages.

Structural Insights on the SARS-CoV-2 Variants of Concern Spike Glycoprotein: A Computational Study With Possible Clinical Implications

Coronavirus disease 2019 (COVID-19) pandemic has been attributed to SARS-CoV-2 (SARS2) and, consequently, SARS2 has evolved into multiple SARS2 variants driving subsequent waves of infections. In particular, variants of concern (VOC) were identified to have both increased transmissibility and virulence ascribable to mutational changes occurring within the spike protein resulting to modifications in the protein structural orientation which in-turn may affect viral pathogenesis. However, this was never fully elucidated. Here, we generated spike models of endemic HCoVs (HCoV 229E, HCoV OC43, HCoV NL63, HCoV HKU1, SARS CoV, MERS CoV), original SARS2, and VOC (alpha, beta, gamma, delta). Model quality check, structural superimposition, and structural comparison based on RMSD values, TM scores, and contact mapping were all performed. We found that: 1) structural comparison between the original SARS2 and VOC whole spike protein model have minor structural differences (TM > 0.98); 2) the whole VOC spike models putatively have higher structural similarity (TM > 0.70) to spike models from endemic HCoVs coming from the same phylogenetic cluster; 3) original SARS2 S1-CTD and S1-NTD models are structurally comparable to VOC S1-CTD (TM = 1.0) and S1-NTD (TM > 0.96); and 4) endemic HCoV S1-CTD and S1-NTD models are structurally comparable to VOC S1-CTD (TM > 0.70) and S1-NTD (TM > 0.70) models belonging to the same phylogenetic cluster. Overall, we propose that structural similarities (possibly ascribable to similar conformational epitopes) may help determine immune cross-reactivity, whereas, structural differences (possibly associated with varying conformational epitopes) may lead to viral infection (either reinfection or breakthrough infection).

Variant Prediction by Analyzing RdRp/S Gene Double or Low Amplification Pattern in Allplex SARS-CoV-2 Assay

The spread of delta variants (B.1.671.2) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a severe global threat. Multiplex real-time PCR is a common method for confirming SARS-CoV-2 infection, however, additional tests, such as whole genomic sequencing, are required to reveal the presence or type of viral mutation. Moreover, applying whole genomic sequencing to all SARS-CoV-2 positive samples is challenging due to time and cost constraints. Here, we report that the double or low amplification curve observed during RNA-dependent RNA polymerase (RdRp) gene/S gene amplification in the Allplex SARS-CoV-2 Assay is related to delta/alpha variants. We analyzed the RdRp/S gene amplification curve using 94 samples confirmed as SARS-CoV-2 infection by the Allplex SARS-CoV-2 Assay from January to August, 2021. These positive samples identified variant types using the Novaplex SARS-CoV-2 Variants I and IV Assays. Overall, 17 samples showing a double curve and 11 samples showing a low amplification pattern were associated with alpha-/delta-type strains with variants in the P681 region. The double or low curve shown in the RdRp gene amplification curve had 100% sensitivity and 100% specificity for diagnosing delta/alpha variants. During the SARS-CoV-2 virus diagnostic RT-PCR test using the Allplex SARS-CoV-2 Assay, we could consider the presence of delta/alpha variants in the samples with double or low amplification curve of the RdRp/S gene channel. This PCR amplification curve abnormality enables rapid and cost-effective variant type prediction during SARS-CoV-2 diagnostic testing in clinical laboratories.