Single point mutations can potentially enhance infectivity of SARS-CoV-2 revealed by in silico affinity maturation and SPR assay

Binding affinity between coronavirus spike protein and human ACE2 receptor

Coronaviruses (CoVs) pose a major risk to global public health due to their ability to infect diverse animal species and potential for emergence in humans. The CoV spike protein mediates viral entry into the cell and plays a crucial role in determining the binding affinity to host cell receptors. With particular emphasis on α- and β-coronaviruses that infect humans and domestic animals, current research on CoV receptor use suggests that the exploitation of the angiotensin-converting enzyme 2 (ACE2) receptor poses a significant threat for viral emergence with pandemic potential. This review summarizes the approaches used to study binding interactions between CoV spike proteins and the human ACE2 (hACE2) receptor. Solid-phase enzyme immunoassays and cell binding assays allow qualitative assessment of binding but lack quantitative evaluation of affinity. Surface plasmon resonance, Bio-layer interferometry, and Microscale Thermophoresis on the other hand, provide accurate affinity measurement through equilibrium dissociation constants (KD). In silico modeling predicts affinity through binding structure modeling, protein-protein docking simulations, and binding energy calculations but reveals inconsistent results due to the lack of a standardized approach. Machine learning and deep learning models utilize simulated and experimental protein-protein interaction data to elucidate the critical residues associated with CoV binding affinity to hACE2. Further optimization and standardization of existing approaches for studying binding affinity could aid pandemic preparedness. Specifically, prioritizing surveillance of CoVs that can bind to human receptors stands to mitigate the risk of zoonotic spillover.

Unraveling the Dynamics of SARS-CoV-2 Mutations: Insights from Surface Plasmon Resonance Biosensor Kinetics

Surface Plasmon Resonance (SPR) technology is known to be a powerful tool for studying biomolecular interactions because it offers real-time and label-free multiparameter analysis with high sensitivity. This article summarizes the results that have been obtained from the use of SPR technology in studying the dynamics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) mutations. This paper will begin by introducing the working principle of SPR and the kinetic parameters of the sensorgram, which include the association rate constant (ka), dissociation rate constant (kd), equilibrium association constant (KA), and equilibrium dissociation constant (KD). At the end of the paper, we will summarize the kinetic data on the interaction between angiotensin-converting enzyme 2 (ACE2) and SARS-CoV-2 obtained from the results of SPR signal analysis. ACE2 is a material that mediates virus entry. Therefore, understanding the kinetic changes between ACE2 and SARS-CoV-2 caused by the mutation will provide beneficial information for drug discovery, vaccine development, and other therapeutic purposes.

Revealing the SARS-CoV-2 Spike Protein and Specific Antibody Immune Complex Formation Mechanism for Precise Evaluation of Antibody Affinity

The profound understanding and detailed evaluation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (SCoV2-S) protein and specific antibody interaction mechanism is of high importance in the development of immunosensors for COVID-19. In the present work, we studied a model system of immobilized SCoV2-S protein and specific monoclonal antibodies by molecular dynamics of immune complex formation in real time. We simultaneously applied spectroscopic ellipsometry and quartz crystal microbalance with dissipation to reveal the features and steps of the immune complex formation. We showed direct experimental evidence based on acoustic and optical measurements that the immune complex between covalently immobilized SCoV2-S and specific monoclonal antibodies is formed in two stages. Based on these findings it was demonstrated that applying a two-step binding mathematical model for kinetics analysis leads to a more precise determination of interaction rate constants than that determined by the 1:1 Langmuir binding model. Our investigation showed that the equilibrium dissociation constants (KD) determined by a two-step binding model and the 1:1 Langmuir model could differ significantly. The reported findings can facilitate a deeper understanding of antigen-antibody immune complex formation steps and can open a new way for the evaluation of antibody affinity towards corresponding antigens.

Free Energy Perturbation Calculations of Mutation Effects on SARS-CoV-2 RBD::ACE2 Binding Affinity

The strength of binding between human angiotensin converting enzyme 2 (ACE2) and the receptor binding domain (RBD) of viral spike protein plays a role in the transmissibility of the SARS-CoV-2 virus. In this study we focus on a subset of RBD mutations that have been frequently observed in infected individuals and probe binding affinity changes to ACE2 using surface plasmon resonance (SPR) measurements and free energy perturbation (FEP) calculations. Our SPR results are largely in accord with previous studies but discrepancies do arise due to differences in experimental methods and to protocol differences even when a single method is used. Overall, we find that FEP performance is superior to that of other computational approaches examined as determined by agreement with experiment and, in particular, by its ability to identify stabilizing mutations. Moreover, the calculations successfully predict the observed cooperative stabilization of binding by the Q498R N501Y double mutant present in Omicron variants and offer a physical explanation for the underlying mechanism. Overall, our results suggest that despite the significant computational cost, FEP calculations may offer an effective strategy to understand the effects of interfacial mutations on protein-protein binding affinities and, hence, in a variety of practical applications such as the optimization of neutralizing antibodies.

Functional evolution of SARS-CoV-2 spike protein: Maintaining wide host spectrum and enhancing infectivity via surface charge of spike protein

The SARS-CoV-2 virus, which causes the COVID-19, is rapidly accumulating mutations to adapt to the hosts. We collected SARS-CoV-2 sequence data from the end of 2019 to January 2023 to analyze for their evolutionary features during the pandemic. We found that most of the SARS-CoV-2 genes are undergoing negative purifying selection, while the spike protein gene (S-gene) is undergoing rapid positive selection. From the original strain to the alpha, delta and omicron variant types, the Ka/Ks of the S-gene increases, while the Ka/Ks within one variant type decreases over time. During the evolution, the codon usage did not evolve towards optimal translation and protein expression. In contrast, only S-gene mutations showed a remarkable trend on accumulating more positive charges. This facilitates the infection via binding human ACE2 for cell entry and binding furin for cleavage. Such a functional evolution emphasizes the survival strategy of SARS-CoV-2, and indicated new druggable target to contain the viral infection. The nearly fully positively-charged interaction surfaces indicated that the infectivity of SARS-CoV-2 virus may approach a limit.

Structural heterogeneity and precision of implications drawn from cryo-electron microscopy structures: SARS-CoV-2 spike-protein mutations as a test case

Protein structures may be used to draw functional implications at the residue level, but how sensitive are these implications to the exact structure used? Calculation of the effects of SARS-CoV-2 S-protein mutations based on experimental cryo-electron microscopy structures have been abundant during the pandemic. To understand the precision of such estimates, we studied three distinct methods to estimate stability changes for all possible mutations in 23 different S-protein structures (3.69 million ΔΔG values in total) and explored how random and systematic errors can be remedied by structure-averaged mutation group comparisons. We show that computational estimates have low precision, due to method and structure heterogeneity making results for single mutations uninformative. However, structure-averaged differences in mean effects for groups of substitutions can yield significant results. Illustrating this protocol, functionally important natural mutations, despite individual variations, average to a smaller stability impact compared to other possible mutations, independent of conformational state (open, closed). In summary, we document substantial issues with precision in structure-based protein modeling and recommend sensitivity tests to quantify these effects, but also suggest partial solutions to the problem in the form of structure-averaged “ensemble” estimates for groups of residues when multiple structures are available.

Spike protein of SARS-CoV-2 variants: a brief review and practical implications

The scientific community has been alarmed by the possible immunological evasion, higher infectivity, and severity of disease caused by the newest variants of SARS-CoV-2. The spike protein has an important role in the cellular invasion of viruses and is the target of several vaccines and therapeutic resources, such as monoclonal antibodies. In addition, some of the most relevant mutations in the different variants are on the spike (S) protein gene sequence that leads to structural alterations in the predicted protein, thus causing concern about the protection mediated by vaccines against these new strains. The present review highlights the most recent knowledge about COVID-19 and vaccines, emphasizing the different spike protein structures of SARS-CoV-2 and updating the reader about the emerging viral variants and their classifications, the more common viral mutations described and their distribution in Brazil. It also compiles a table with the most recent knowledge about all of the Omicron spike mutations.

Structural bioinformatics analysis of SARS-CoV-2 variants reveals higher hACE2 receptor binding affinity for Omicron B.1.1.529 spike RBD compared to wild type reference

To date, more than 263 million people have been infected with SARS-CoV-2 during the COVID-19 pandemic. In many countries, the global spread occurred in multiple pandemic waves characterized by the emergence of new SARS-CoV-2 variants. Here we report a sequence and structural-bioinformatics analysis to estimate the effects of amino acid substitutions on the affinity of the SARS-CoV-2 spike receptor binding domain (RBD) to the human receptor hACE2. This is done through qualitative electrostatics and hydrophobicity analysis as well as molecular dynamics simulations used to develop a high-precision empirical scoring function (ESF) closely related to the linear interaction energy method and calibrated on a large set of experimental binding energies. For the latest variant of concern (VOC), B.1.1.529 Omicron, our Halo difference point cloud studies reveal the largest impact on the RBD binding interface compared to all other VOC. Moreover, according to our ESF model, Omicron achieves a much higher ACE2 binding affinity than the wild type and, in particular, the highest among all VOCs except Alpha and thus requires special attention and monitoring.

Design and simulation of a millifluidic device for differential detection of SARS-CoV-2 and H1N1 based on triboelectricity

Differential diagnosis of pathogenic diseases, presently coronavirus disease 2019 (COVID-19) and influenza, is crucial with due attention to their superspreading events, presumably long incubation period, particular complications, and treatments. In this paper, a label-free, self-powered, and ultrafast immunosensor device working based on triboelectric effect is proposed. Equilibrium constants of specific antibody-antigen reactions are accompanied by IEP-relevant electric charges of antigens to recognize SARS-CoV-2 and H1N1. Simulation attributes including fluid flow and geometrical parameters are optimized so that the maximum capture efficiency of 85.63% is achieved. Accordingly, antibody-antigen complexes form electric double layers (EDLs) across the channel interfaces. The resultant built-in electric field affects the following external electric field derived from triboelectricity, leading to the variation of open-circuit voltage as a sensing metric. The device is flexible to operate in conductor-to-dielectric single-electrode and contact-separation modes simultaneously. While the detection limit is reduced utilizing the single-electrode mode compared to the latter one, surface treatment of the triboelectric pair contributes to the sensitivity enhancement. A threshold value equal to −4.113 V is featured to discriminate these two viruses in a vast detectable region; however, further surface engineering can allow the on-site detection of any electrically-charged pathogen applying the emerging triboelectric immunosensor enjoying a lower detection limit.

Epistasis at the SARS-CoV-2 Receptor-Binding Domain Interface and the Propitiously Boring Implications for Vaccine Escape

At the time of this writing, December 2021, potential emergence of vaccine escape variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a grave global concern. The interface between the receptor-binding domain (RBD) of SARS-CoV-2 spike (S) protein and the host receptor (ACE2) overlaps the binding site of principal neutralizing antibodies (NAb), limiting the repertoire of viable mutations. Nonetheless, variants with multiple RBD mutations have risen to dominance. Nonadditive, epistatic relationships among RBD mutations are apparent, and assessing the impact of such epistasis on the mutational landscape, particularly the risk of vaccine escape, is crucial. We employed protein structure modeling using Rosetta to compare the effects of all single mutants at the RBD-NAb and RBD-ACE2 interfaces for the wild type and Delta, Gamma, and Omicron variants. Overall, epistasis at the RBD interface appears to be limited, and the effects of most multiple mutations are additive. Epistasis at the Delta variant interface weakly stabilizes NAb interaction relative to ACE2 interaction, whereas in Gamma, epistasis more substantially destabilizes NAb interaction. Despite bearing many more RBD mutations, the epistatic landscape of Omicron closely resembles that of Gamma. Thus, although Omicron poses new risks not observed with Delta, structural constraints on the RBD appear to hamper continued evolution toward more complete vaccine escape. The modest ensemble of mutations relative to the wild type that are currently known to reduce vaccine efficacy is likely to contain the majority of all possible escape mutations for future variants, predicting the continued efficacy of the existing vaccines. IMPORTANCE Emergence of vaccine escape variants of SARS-CoV-2 is arguably the most pressing problem during the COVID-19 pandemic as vaccines are distributed worldwide. We employed a computational approach to assess the risk of antibody escape resulting from mutations in the receptor-binding domain of the spike protein of the wild-type SARS-CoV-2 virus as well as the Delta, Gamma, and Omicron variants. The efficacy of the existing vaccines against Omicron could be substantially reduced relative to the wild type, and the potential for vaccine escape is of grave concern. Our results suggest that although Omicron poses new evolutionary risks not observed for Delta, structural constraints on the RBD make continued evolution toward more complete vaccine escape from either Delta or Omicron unlikely. The modest set of escape-enhancing mutations already identified for the wild type likely include the majority of all possible mutations with this effect.

Epistasis at the SARS-CoV-2 RBD Interface and the Propitiously Boring Implications for Vaccine Escape

At the time of this writing, December 2021, potential emergence of vaccine escape variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a grave global concern. The interface between the receptor-binding domain (RBD) of SARS-CoV-2 spike (S) protein and the host receptor (ACE2) overlap with the binding site of principal neutralizing antibodies (NAb), limiting the repertoire of viable mutations. Nonetheless, variants with multiple mutations in the RBD have rose to dominance. Non-additive, epistatic relationships among RBD mutations are apparent, and assessing the impact of such epistasis on the mutational landscape is crucial. Epistasis can substantially increase the risk of vaccine escape and cannot be completely characterized through the study of the wild type (WT) alone. We employed protein structure modeling using Rosetta to compare the effects of all single mutants at the RBD-NAb and RBD-ACE2 interfaces for the WT, Delta, Gamma, and Omicron variants. Overall, epistasis at the RBD interface appears to be limited and the effects of most multiple mutations are additive. Epistasis at the Delta variant interface weakly stabilizes NAb interaction relative to ACE2 interaction, whereas in the Gamma variant, epistasis more substantially destabilizes NAb interaction. Although a small, systematic trend towards NAb destabilization not observed for Delta or Gamma was detected for Omicron, and despite bearing significantly more RBD mutations, the epistatic landscape of the Omicron variant closely resembles that of Gamma. These results suggest that, although Omicron poses new risks not observed with Delta, structural constraints on the RBD hamper continued evolution towards more complete vaccine escape. The modest ensemble of mutations relative to the WT that are currently known to reduce vaccine efficacy is likely to comprise the majority of all possible escape mutations for future variants, predicting continued efficacy of the existing vaccines.

Molecular Level Dissection of Critical Spike Mutations in SARS-CoV-2 Variants of Concern (VOCs): A Simplified Review

SARS-CoV-2 virus during its spread in the last one and half year has picked up critical changes in its genetic code i.e. mutations, which have leads to deleterious epidemiological characteristics. Due to critical role of spike protein in cell entry and pathogenesis, mutations in spike regions have been reported to enhance transmissibility, disease severity, possible escape from vaccine-induced immune response and reduced diagnostic sensitivity/specificity. Considering the structure-function impact of mutations, understanding the molecular details of these key mutations of newly emerged variants/lineages is of urgent concern. In this review, we have explored the literature on key spike mutations harbored by alpha, beta, gamma and delta 'variants of concern' (VOCs) and discussed their molecular consequences in the context of resultant virus biology. Commonly all these VOCs i.e. B.1.1.7, B.1.351, P.1 and B.1.617.2 lineages have decisive mutation in Receptor Binding Motif (RBM) region and/or region around Furin cleavage site (FCS) of spike protein. In general, mutation induced disruption of intra-molecular interaction enhances molecular flexibility leading to exposure of spike protein surface in these lineages to make it accessible for inter-molecular interaction with hACE2. A disruption of spike antigen-antibody inter-molecular interactions in epitope region due to the chemical nature of substituting amino acid hampers the neutralization efficacy. Simplified surveillance of mutation induced changes and their consequences at molecular level can contribute in rationalizing mutation's impact on virus biology. It is believed that molecular level dissection of these key spike mutation will assist the future investigations for a more resilient outcome against severity of COVID-19.