SARS-CoV-2 Omicron Variant Binds to Human Cells More Strongly than the Wild Type: Evidence from Molecular Dynamics Simulation

SARS-CoV-2 Omicron Subvariants Do Not Differ Much in Binding Affinity to Human ACE2: A Molecular Dynamics Study

The emergence of the variant of concern Omicron (B.1.1.529) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) exacerbates the COVID-19 pandemic due to its high contagious ability. Studies have shown that the Omicron binds human ACE2 more strongly than the wild type. The prevalence of Omicron in new cases of COVID-19 promotes novel lineages with improved receptor binding affinity and immune evasion. To shed light on this open problem, in this work, we investigated the binding free energy of the receptor binding domain of the Omicron lineages BA.2, BA.2.3.20, BA.3, BA4/BA5, BA.2.75, BA.2.75.2, BA.4.6, XBB.1, XBB.1.5, BJ.1, BN.1, BQ.1.1, and CH.1.1 to human ACE2 using all-atom molecular dynamics simulation and the molecular mechanics Poisson-Boltzmann surface area method. The results show that these lineages have increased binding affinity compared to the BA.1 lineage, and BA.2.75 and BA.2.75.2 subvariants bind ACE2 more strongly than others. However, in general, the binding affinities of the Omicron lineages do not differ significantly from each other. The electrostatic force dominates over the van der Waals force in the interaction between Omicron lineages and human cells. Based on our results, we argue that viral evolution does not further improve the affinity of SARS-CoV-2 for ACE2 but may increase immune evasion.

Unraveling the binding mechanisms of SARS-CoV-2 variants through molecular simulations

The emergence of SARS-CoV-2 variants like Delta (AY.29) and Omicron (EG.5) poses continued challenges for vaccines and therapeutics. Mutations in the viral spike protein are key in altering infectivity and immune evasion. This study uses computational modeling to investigate the molecular binding mechanisms between spike protein variants and the ACE2 host receptor. Using the MARTNI force field, coarse-grained molecular dynamics (CGMD) simulations and nudged elastic band (NEB) calculations explore spike-ACE2 interactions for the wild type, Delta variant, and Omicron variant. The simulations reveal Omicron has the strongest binding affinity at -128.35 ± 10.91 kcal/mol, followed by Delta and wild type. Key mutations in Delta and Omicron, like Q493R and Q498R, optimize electrostatic contacts, enhancing ACE2 interactions. The wild-type spike has the highest transition state energy barrier at 17.87 kcal/mol, while Delta has the lowest barrier at 9.21 kcal/mol. Despite slightly higher dual barriers, Omicron's increased binding energy lowers its overall barrier to rapidly bind ACE2. These findings provide residue-level insights into mutation effects on SARS-CoV-2 infectivity. The computational modeling elucidates mechanisms underlying spike-ACE2 binding kinetics, aiding the development of vaccines and therapies targeting emerging viral strains.

Deciphering the free energy landscapes of SARS-CoV-2 wild type and Omicron variant interacting with human ACE2

The binding of the receptor binding domain (RBD) of the SARS-CoV-2 spike protein to the host cell receptor angiotensin-converting enzyme 2 (ACE2) is the first step in human viral infection. Therefore, understanding the mechanism of interaction between RBD and ACE2 at the molecular level is critical for the prevention of COVID-19, as more variants of concern, such as Omicron, appear. Recently, atomic force microscopy has been applied to characterize the free energy landscape of the RBD-ACE2 complex, including estimation of the distance between the transition state and the bound state, xu. Here, using a coarse-grained model and replica-exchange umbrella sampling, we studied the free energy landscape of both the wild type and Omicron subvariants BA.1 and XBB.1.5 interacting with ACE2. In agreement with experiment, we find that the wild type and Omicron subvariants have similar xu values, but Omicron binds ACE2 more strongly than the wild type, having a lower dissociation constant KD.

Comparative Computational Analysis of Spike Protein Structural Stability in SARS-CoV-2 Omicron Subvariants

The continuous emergence of new severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants with multiple spike (S) protein mutations pose serious threats to current coronavirus disease 2019 (COVID-19) therapies. A comprehensive understanding of the structural stability of SARS-CoV-2 variants is vital for the development of effective therapeutic strategies as it can offer valuable insights into their potential impact on viral infectivity. S protein mediates a virus' attachment to host cells by binding to angiotensin-converting enzyme 2 (ACE2) through its receptor-binding domain (RBD), and mutations in this protein can affect its stability and binding affinity. We analyzed S protein structural stability in various Omicron subvariants computationally. Notably, the S protein sequences analyzed in this work were obtained directly from our own sample collection. We evaluated the binding free energy between S protein and ACE2 in several complex forms. Additionally, we measured distances between the RBD of each chain in S protein to analyze conformational changes. Unlike most of the prior studies, we analyzed full-length S protein-ACE2 complexes instead of only RBD-ACE2 complexes. Omicron subvariants including BA.1, BA.2, BA.2.12.1, BA.4/BA.5, BA.2.75, BA.2.75_K147E, BA.4.6 and BA.4.6_N658S showed enhanced stability compared to wild type, potentially due to distinct S protein mutations. Among them, BA.2.75 and BA.4.6_N658S exhibited the highest and lowest level of stability, respectively.

Identification of key mutations responsible for the enhancement of receptor-binding affinity and immune escape of SARS-CoV-2 Omicron variant

The Omicron variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has raised concerns worldwide due to its enhanced transmissibility and immune escapability. The first dominant Omicron BA.1 subvariant harbors more than 30 mutations in the spike protein from the prototype virus, of which 15 mutations are located at the receptor binding domain (RBD). These mutations in the RBD region attracted significant attention, which potentially enhance the binding of the receptor human angiotensin-converting enzyme 2 (hACE2) and decrease the potency of neutralizing antibodies/nanobodies. This study applied the molecular dynamics simulations combined with the molecular mechanics-generalized Born surface area (MMGBSA) method, to investigate the molecular mechanism behind the impact of the mutations acquired by Omicron on the binding affinity between RBD and hACE2. Our results indicate that five key mutations, i.e., N440K, T478K, E484A, Q493R, and G496S, contributed significantly to the enhancement of the binding affinity by increasing the electrostatic interactions of the RBD-hACE2 complex. Moreover, fourteen neutralizing antibodies/nanobodies complexed with RBD were used to explore the effects of the mutations in Omicron RBD on their binding affinities. The calculation results indicate that the key mutations E484A and Y505H reduce the binding affinities to RBD for most of the studied neutralizing antibodies/nanobodies, mainly attributed to the elimination of the original favorable gas-phase electrostatic and hydrophobic interactions between them, respectively. Our results provide valuable information for developing effective vaccines and antibody/nanobody drugs.

Investigating the role of glycans in Omicron sub-lineages XBB.1.5 and XBB.1.16 binding to host receptor using molecular dynamics and binding free energy calculations

Omicron derived lineages viz. BA.2, BA.3, BA.4 BA.5, BF.7 and XBBs show prominence with improved immune escape, transmissibility, infectivity, and pathogenicity in general. Sub-variants, XBB.1.5 and XBB.1.16 have shown rapid spread, with mutations embedded throughout the viral genome, including the spike protein. Changing atomic landscapes in spike contributes significantly to modulate host pathogen interactions and infections thereof. In the present work, we computationally analyzed the binding affinities of spike receptor binding domains (RBDs) of XBB.1.5 and XBB.1.16 towards human angiotensin-converting enzyme 2 (hACE2) compared to Omicron. We have employed simulations and binding energy estimation of molecular complexes of spike-hACE2 to assess the interplay of interaction pattern and effect of mutations if any in the binding mode of the RBDs of these novel mutants. We calculated the binding free energy (BFE) of the RBD of the Omicron, XBB.1.5 and XBB.1.16 spike protein to hACE2. We showed that XBB.1.5 and XBB.1.16 can bind to human cells more strongly than Omicron due to the increased charge of the RBD, which enhances the electrostatic interactions with negatively charged hACE2. The per-residue decompositions further show that the Asp339His, Asp405Asn and Asn460Lys mutations in the XBBs RBD play a crucial role in enhancing the electrostatic interactions, by acquiring positively charged residues, thereby influencing the formation/loss of interfacial bonds and thus strongly affecting the spike RBD-hACE2 binding affinity. Simulation results also indicate less interference of heterogeneous glycans of XBB.1.5 spike RBD towards binding to hACE2. Moreover, despite having less interaction at the three interfacial contacts between XBB S RBD and hACE2 compared to Omicron, variants XBB.1.5 and XBB.1.16 had higher total binding free energies (ΔGbind) than Omicron due to the contribution of non-interfacial residues to the free energy, providing insight into the increased binding affinity of XBB1.5 and XBB.1.16. Furthermore, the presence of large positively charged surface patches in the XBBs act as drivers of electrostatic interactions, thus support the possibility of a higher binding affinity to hACE2.

Molecular insights and optimization strategies for the competitive binding of engineered ACE2 proteins: a multiple replica molecular dynamics study

The COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) continues to spread globally, and rapid viral evolution and the emergence of new variants pose challenges to pandemic control. During infection, the spike protein of SARS-CoV-2 interacts with the human ACE2 protein via its receptor binding domain (RBD), and it is known that engineered forms of ACE2 can compete with wild-type (WT) ACE2 for binding to inhibit infection. Here, we conducted multiple replica molecular dynamics (MRMD) simulations to study the mechanisms of the engineered ACE2 variants 3N39 and 3N94 and provide directions for optimization. Our findings reveal that engineered ACE2 is notably more efficacious in systems that show weaker binding to WT ACE2 (i.e., WT and BA.1 RBD), but also faces immune escape as the virus evolves. Moreover, by modifying residue types near the binding interface, engineered ACE2 alters the electrostatic potential distribution and reconfigures the hydrogen bonding network, which results in modified binding to the RBD. However, this structural rearrangement does not occur in all RBD variants. In addition, we identified potentially engineerable beneficial residues and potentially engineerable detrimental residues in both ACE2 and RBD. Functional conservation can thus enable the optimization of these residues and improve the binding competitiveness of engineered ACE2, which therefore provides additional immune escape prevention. Finally, we conclude that these findings have implications for understanding the mechanisms responsible for engineered ACE2 and can help us to develop engineered ACE2 proteins that show superior performance.

Physicochemical effects of emerging exchanges on the spike protein's RBM of the SARS-CoV-2 Omicron subvariants BA.1-BA.5 and its influence on the biological properties and attributes developed by these subvariants

Emerging in South Africa, SARS-CoV-2 Omicron variant was marked by the expression of an exaggerated number of mutations throughout its genome and by the emergence of subvariants, whose attributes developed by them have been associated with amino acid exchanges that occur mainly in the RBM region of the spike protein. The RBM comprises a region within the RBD and is directly involved in the SARS-CoV-2 spike protein interaction with the host cell ACE2 receptor, during the infection mechanism and viral transmission. Defined as the region from aa 437 to aa 508, there are several residues in certain positions that interact directly with the human ACE-2 receptor during these processes. The occurrence of amino acid exchanges in these positions causes physicochemical alterations in the SARS-CoV-2 spike protein, which confer additional advantages and attributes to the agent. In addition, these exchanges serve as a basis for the characterization of new variants and subvariants of SARS-CoV-2. In this review, the amino acid exchanges that have occurred in the RBM of the subvariants BA.1 to BA.5 of SARS-CoV-2 that emerged from the Omicron are described. The physicochemical effects caused by them on spike protein are also described, as well as their influence on the biological properties and attributes developed by the subvariants BA.1, BA.2, BA.3, BA.4 and BA.5.

Interaction of SARS-CoV-2 with host cells and antibodies: experiment and simulation

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the devastating global COVID-19 pandemic announced by WHO in March 2020. Through unprecedented scientific effort, several vaccines, drugs and antibodies have been developed, saving millions of lives, but the fight against COVID-19 continues as immune escape variants of concern such as Delta and Omicron emerge. To develop more effective treatments and to elucidate the side effects caused by vaccines and therapeutic agents, a deeper understanding of the molecular interactions of SARS-CoV-2 with them and human cells is required. With special interest in computational approaches, we will focus on the structure of SARS-CoV-2 and the interaction of its spike protein with human angiotensin-converting enzyme-2 (ACE2) as a prime entry point of the virus into host cells. In addition, other possible viral receptors will be considered. The fusion of viral and human membranes and the interaction of the spike protein with antibodies and nanobodies will be discussed, as well as the effect of SARS-CoV-2 on protein synthesis in host cells.

Different aspects in explaining how mutations could affect the binding mechanism of receptor binding domain of SARS-CoV-2 spike protein in interaction with ACE2

During replication, some mutations occur in SARS-CoV-2, the causal agent of COVID-19, leading to the emergence of different variants of the virus. The mutations that accrue in different variants of the virus, influence the virus' ability to bind to human cell receptors and ability to evade the human immune system, the rate of viral transmission, and effectiveness of vaccines. Some of these mutations occur in the receptor binding domain (RBD) of the spike protein that may change the affinity of the virus for the ACE2 receptor. In this study, several in silico techniques, such as MD and SMD simulations, were used to perform comparative studies to deeply understand the effect of mutation on structural and functional details of the interaction of the spike glycoprotein of SARS-CoV-2, with the ACE2 receptor. According to our results, the mutation in the RBD associated with the Omicron variant increase binding affinity of the virus to ACE2 when compared to wild type and Delta variants. We also observed that the flexibility of the spike protein of the Omicron variant was lower than in comparison to other variants. In summary, different mutations in variants of the virus can have an effect on the binding mechanism of the receptor binding domain of the virus with ACE2.

A computationally designed ACE2 decoy has broad efficacy against SARS-CoV-2 omicron variants and related viruses in vitro and in vivo

SARS-CoV-2, especially B.1.1.529/omicron and its sublineages, continues to mutate to evade monoclonal antibodies and antibodies elicited by vaccination. Affinity-enhanced soluble ACE2 (sACE2) is an alternative strategy that works by binding the SARS-CoV-2 S protein, acting as a 'decoy' to block the interaction between the S and human ACE2. Using a computational design strategy, we designed an affinity-enhanced ACE2 decoy, FLIF, that exhibited tight binding to SARS-CoV-2 delta and omicron variants. Our computationally calculated absolute binding free energies (ABFE) between sACE2:SARS-CoV-2 S proteins and their variants showed excellent agreement to binding experiments. FLIF displayed robust therapeutic utility against a broad range of SARS-CoV-2 variants and sarbecoviruses, and neutralized omicron BA.5 in vitro and in vivo. Furthermore, we directly compared the in vivo therapeutic efficacy of wild-type ACE2 (non-affinity enhanced ACE2) against FLIF. A few wild-type sACE2 decoys have shown to be effective against early circulating variants such as Wuhan in vivo. Our data suggest that moving forward, affinity-enhanced ACE2 decoys like FLIF may be required to combat evolving SARS-CoV-2 variants. The approach described herein emphasizes how computational methods have become sufficiently accurate for the design of therapeutics against viral protein targets. Affinity-enhanced ACE2 decoys remain highly effective at neutralizing omicron subvariants.

The impact of mutation sets in receptor-binding domain of SARS-CoV-2 variants on the stability of RBD–ACE2 complex

Aim: Bioinformatic analysis of mutation sets in receptor-binding domain (RBD) of currently and previously circulating SARS-CoV-2 variants of concern (VOCs) and interest (VOIs) to assess their ability to bind the ACE2 receptor. Methods: In silico sequence and structure-oriented approaches were used to evaluate the impact of single and multiple mutations. Results: Mutations detected in VOCs and VOIs led to the reduction of binding free energy of the RBD–ACE2 complex, forming additional chemical bonds with ACE2, and to an increase of RBD–ACE2 complex stability. Conclusion: Mutation sets characteristic of SARS-CoV-2 variants have complex effects on the ACE2 receptor-binding affinity associated with amino acid interactions at mutation sites, as well as on the acquisition of other viral adaptive advantages.

Is BF.7 more infectious than other Omicron subtypes: Insights from structural and simulation studies of BF.7 spike RBD variant

Fear of a fresh infection wave and a global health issue in the ongoing COVID-19 pandemic have been rekindled by the appearance of two new novel variants BF.7 and BA.4/5 of Omicron lineages. Predictions of increased antibody evasion capabilities and transmissibility have been recognised in addition to the existing lineages (BA.1.1, BA.2, BA.2.12.1 and BA.3) as cause for worry. In comparison to Omicron, BA.4 and BF.7 share nine mutations in the spike protein, Leu371Phe, Thr376Ala, Asp405Asn, Arg408Ser, Ser446Gly, Leu452Arg, Phe486Val, Arg493Gln, Ser496Gly, whereas BF.7 contains an additional mutation, Arg346Thr, in the receptor binding domain (RBD) region. Due to the critical need for analysis and data on the BA.4 and BF.7 variants, we have computationally analyzed the interaction pattern between the Omicron, BA.4 and BF.7 RBD and angiotensin-converting enzyme 2 (ACE2) to determine the influence of these unique mutations on the structures, functions, and binding affinity of RBD towards ACE2. These analyses also allow to compare molecular models to previously reported data to evaluate the robustness of our methods for quick prediction of emerging future variants. The docking results reveal that BA.4 and BF.7 have particularly strong interactions with ACE2 when compared to Omicron, as shown by several parameters such as salt bridge, hydrogen bond, and non-bonded interactions. In addition, the estimations of binding free energy corroborated the findings further. BA.4 and BF.7 were found to bind to ACE2 with similar affinities (-72.14 and - 71.54 kcal/mol, respectively) and slightly stronger than Omicron (-70.04 kcal/mol). The differences in the binding pattern between the Omicron, BA.4 and BF.7 variant complexes indicated that the BA.4 and BF.7 RBD substitutions Asp405Asn, Ser446Gly, Leu452Arg, Phe486Val and Arg493Gln caused additional interactions with ACE2. In addition, normal mode analyses also indicate more stable conformations of BA.4 and BF.7 RBDs against human ACE2. Based on these structural and simulation analyses, we hypothesized that these changes may affect the binding affinity of BA.4 and BF.7 with ACE2.

Conformational Behavior of SARS-Cov-2 Spike Protein Variants: Evolutionary Jumps in Sequence Reverberate in Structural Dynamic Differences

SARS-CoV-2 has evolved rapidly in the first 3 years of pandemic diffusion. The initial evolution of the virus appeared to proceed through big jumps in sequence changes rather than through the stepwise accumulation of point mutations on already established variants. Here, we examine whether this nonlinear mutational process reverberates in variations of the conformational dynamics of the SARS-CoV-2 Spike protein (S-protein), the first point of contact between the virus and the human host. We run extensive microsecond-scale molecular dynamics simulations of seven distinct variants of the protein in their fully glycosylated state and set out to elucidate possible links between the mutational spectrum of the S-protein and the structural dynamics of the respective variant, at global and local levels. The results reveal that mutation-dependent structural and dynamic modulations mostly consist of increased coordinated motions in variants that acquire stability and in an increased internal flexibility in variants that are less stable. Importantly, a limited number of functionally important substructures (the receptor binding domain, in particular) share the same time of movements in all variants, indicating efficient preorganization for functional regions dedicated to host interactions. Our results support a model in which the internal dynamics of the S-proteins from different strains varies in a way that reflects the observed random and non-stepwise jumps in sequence evolution, while conserving the functionally oriented traits of conformational dynamics necessary to support productive interactions with host receptors.

Modelling SARS-CoV-2 spike-protein mutation effects on ACE2 binding

The binding affinity of the SARS-CoV-2 spike (S)-protein to the human membrane protein ACE2 is critical for virus function. Computational structure-based screening of new S-protein mutations for ACE2 binding lends promise to rationalize virus function directly from protein structure and ideally aid early detection of potentially concerning variants. We used a computational protocol based on cryo-electron microscopy structures of the S-protein to estimate the change in ACE2-affinity due to S-protein mutation (ΔΔG_bind) in good trend agreement with experimental ACE2 affinities. We then expanded predictions to all possible S-protein mutations in 21 different S-protein-ACE2 complexes (400,000 ΔΔG_bind data points in total), using mutation group comparisons to reduce systematic errors. The results suggest that mutations that have arisen in major variants as a group maintain ACE2 affinity significantly more than random mutations in the total protein, at the interface, and at evolvable sites. Omicron mutations as a group had a modest change in binding affinity compared to mutations in other major variants. The single-mutation effects seem consistent with ACE2 binding being optimized and maintained in omicron, despite increased importance of other selection pressures (antigenic drift), however, epistasis, glycosylation and in vivo conditions will modulate these effects. Computational prediction of SARS-CoV-2 evolution remains far from achieved, but the feasibility of large-scale computation is substantially aided by using many structures and mutation groups rather than single mutation effects, which are very uncertain. Our results demonstrate substantial challenges but indicate ways forward to improve the quality of computer models for assessing SARS-CoV-2 mutation effects.

SARS-CoV-2 and its impact on the cardiovascular and digestive systems - The interplay between new virus variants and human cells

Since infection with the novel coronavirus SARS-CoV-2 first emerged in Wuhan, China, in December 2019, the world has been battling the pandemic COVID-19. Patients of all ages and genders are now becoming infected with the new coronavirus variant (Omicron) worldwide, and its subvariants continue to pose a threat to health and life. This article provides a literature review of cardiovascular and gastrointestinal complications resulting from SARS-CoV-2 infection. COVID-19 primarily caused respiratory symptoms, but complications can affect many vital organs. SARS-CoV-2 binds to a human cell receptor (angiotensin-converting enzyme 2 - ACE2) that is predominantly expressed primarily in the heart and gastrointestinal tract, which is why we focused on complications in these organs. Since the high transmissibility of Omicron and its ability to evade the immune system have raised worldwide concern, we have tried to summarise the current knowledge about its development from a structural point of view and to highlight the differences in its binding to human receptors and proteases compared to previous VOC.

Atomistic insights into the binding of SARS-CoV-2 spike receptor binding domain with the human ACE2 receptor: The importance of residue 493

SARS-CoV-2 is a coronavirus that has created a global pandemic. The virus contains a spike protein which has been shown to bind to the ACE2 receptor on the surface of human cells. Vaccines have been developed that recognize elements of the SARS-CoV-2 spike protein and they have been successful in preventing infection. Recently, the Omicron variant of the SARS-CoV-2 virus was reported and quickly became a variant of concern due to its transmissibility. This variant contained an unusually large number (32) of point mutations, of which 15 of those mutations are in the receptor binding domain of the spike protein. While several computational and experimental investigations comparing the binding of the Omicron and wild type RBD to the human ACE2 receptor have been conducted, many of these report contradictory findings. In order to assess the differential binding ability, we conducted 2 μs of classical molecular dynamics (cMD) simulation to estimate the binding affinities and behaviors. Based upon MM-GBSA binding affinity, per-residue energy decomposition analysis, center of mass distance measurements, ensemble clustering, pairwise residue decomposition and hydrogen bonding analysis, our results suggest that a single point mutation is responsible for the enhanced binding of the Omicron mutant relative to the WT. While the 15-point mutations in the receptor binding domain contribute positively and negatively to the affinity of the spike protein for the human ACE2 receptor, it is the point mutation Q493R that confers enhanced binding while the Q493K mutation results in similar binding. The MM-GBSA binding estimations over a 2 μs trajectory, suggest that the wild type binds to ACE2 with a value of -29.69 kcal/mol while the Q493K and Q493R Omicron mutants bind with energy values of -26.67 and -34.56 kcal/mol, respectively. These values are significantly different, given the error estimates associated with the MM-GBSA method. In general, while some mutations increase binding, more mutations diminish binding, leading to an overall similar picture of binding for Q493K and enhanced binding for Q493R.

Conformational Dynamics of the Receptor-Binding Domain of the SARS-CoV-2 Spike Protein

Variants of SARS-CoV-2 keep emerging and causing new waves of COVID-19 around the world. Effective new approaches in drug development are based on the binding of agents, such as neutralizing monoclonal antibodies to a receptor-binding domain (RBD) of SARS-CoV-2 spike protein. However, mutations in RBD may lower the affinity of previously developed antibodies. Therefore, rapid analysis of new variants and selection of a binding partner with high affinity is of great therapeutic importance. Here, we explore a computational approach based on molecular dynamics simulations and conformational clusterization techniques for the wild-type and omicron variants of RBD. Biochemical experiments support the hypothesis of the presence of several conformational states within the RBD assembly. The development of such an approach will facilitate the selection of neutralization drugs with higher affinity based on the primary structure of the target antigen.

Structural Analysis of the Spike Protein of SARS-CoV-2 Variants and Other Betacoronaviruses Using Molecular Dynamics

A structural analysis over various spike proteins from three highly pathogenic Betacoronavirus was done to understand their structural differences. The proteins were modeled using crystal structures from SARS-CoV, MERS-CoV, and other Betacoronavirus that infect bats and pangolins. The group was split in two sets; the first set corresponds to the non-mutated spike proteins, while the second set corresponds to mutated spike variants alpha, beta, gamma, delta, omicron and mu; five of them classified as variants of concern and the last one as variant of interest. A conformational space exploration was carried out for every protein by using molecular dynamic simulations. Root mean square fluctuations, principal component and cross-correlation analysis were carried out over the dynamics to analyze the flexibility and rigidity of every protein in comparison to the wild type Spike protein from the SARS-CoV-2. The obtained results indicate that the proteins, which are not spread among humans, have smooth movements compared to those of SARS-CoV-2 and its variants. In addition, a relationship between the speed of the virulence and the movement of the protein can explain the behavior of delta and omicron variants.

Electrostatic Interactions Are the Primary Determinant of the Binding Affinity of SARS-CoV-2 Spike RBD to ACE2: A Computational Case Study of Omicron Variants

To explore the mechanistic origin that determines the binding affinity of SARS-CoV-2 spike receptor binding domain (RBD) to human angiotensin converting enzyme 2 (ACE2), we constructed the homology models of RBD-ACE2 complexes of four Omicron subvariants (BA.1, BA.2, BA.3 and BA.4/5), and compared them with wild type complex (RBD_WT-ACE2) in terms of various structural dynamic properties by molecular dynamics (MD) simulations and binding free energy (BFE) calculations. The results of MD simulations suggest that the RBDs of all the Omicron subvariants (RBD_OMIs) feature increased global structural fluctuations when compared with RBD_WT. Detailed comparison of BFE components reveals that the enhanced electrostatic attractive interactions are the main determinant of the higher ACE2-binding affinity of RBD_OMIs than RBD_WT, while the weakened electrostatic attractive interactions determine RBD of BA.4/5 subvariant (RBD_BA.4/5) lowest ACE2-binding affinity among all Omicron subvariants. The per-residue BFE decompositions and the hydrogen bond (HB) networks analyses indicate that the enhanced electrostatic attractive interactions are mainly through gain/loss of the positively/negatively charged residues, and the formation or destruction of the interfacial HBs and salt bridges can also largely affect the ACE2-binding affinity of RBD. It is worth pointing out that since Q493R plays the most important positive contribution in enhancing binding affinity, the absence of this mutation in RBD_BA.4/5 results in a significantly weaker binding affinity to ACE2 than other Omicron subvariants. Our results provide insight into the role of electrostatic interactions in determining of the binding affinity of SARS-CoV-2 RBD to human ACE2.

Genome characterization, phylogenomic assessment and spatio-temporal dynamics study of highly mutated BA variants from India

Purpose

The emergence of highly mutated and transmissible BA variants has caused an unprecedented surge in COVID-19 infections worldwide. Thorough analysis of its genome structure and phylogenomic evolutionary details will serve as scientific reference for future research.

Method

Here, we have analyzed the BA variants from India using whole-genome sequencing, spike protein mutation study, spatio-temporal surveillance, phylogenomic assessment and epitope mapping.

Results

The predominance of BA.2/BA.2-like was observed in India during COVID-19 third wave. Genome analysis and mutation study highlighted the existence of 2128 amino acid changes within BA as compared to NC_045512.2. Presence of 23 unknown mutation sites (spanning region 61–831) were observed among the Indian BA variants as compared to the global BA strains. Unassigned probable Omicron showed the highest number of mutations (370) followed by BA.1 (104), BA.2.3 (56), and BA.2 (27). Presence of mutations ‘Q493R ​+ ​Q498R ​+ ​N501Y’, and ‘K417 ​N ​+ ​E484A ​+ ​N501Y’ remained exclusive to BA.2 as well as unassigned probable Omicron. The time-tree and phylogenomic network assessed the evolutionary relationship of the BA variants. Existence of 424 segregating sites and 113 parsimony informative sites within BA genomes were observed through haplotype network analysis. Epitope mapping depicted the presence of unique antigenic sites within the receptor binding domain of the BA variants that could be exploited for robust vaccine development.

Conclusion

These findings provide important scientific insights about the nature, diversity, and evolution of Indian BA variants. The study further divulges in the avenues of therapeutic upgradation for better management and eventual eradication of COVID-19.

When the Dust Has Settled: Calculation of Binding Affinities from First Principles for SARS-CoV-2 Variants with Quantitative Accuracy

Accurate determination of binding free energy is pivotal for the study of many biological processes and has been applied in a number of theoretical investigations to compare the affinity of severe acute respiratory syndrome coronavirus 2 variants toward the host cell. Diversity of these variants challenges the development of effective general therapies, their transmissibility relying either on an increased affinity toward their dedicated human receptor, the angiotensin-converting enzyme 2 (ACE2), or on escaping the immune response. Now that robust structural data are available, we have determined with utmost accuracy the standard binding free energy of the receptor-binding domain to the most widespread variants, namely, Alpha, Beta, Delta, and Omicron BA.2, as well as the wild type (WT) in complex either with ACE2 or with antibodies, namely, S2E12 and H11-D4, using a rigorous theoretical framework that combines molecular dynamics and potential-of-mean-force calculations. Our results show that an appropriate starting structure is crucial to ensure appropriate reproduction of the binding affinity, allowing the variants to be compared. They also emphasize the necessity to apply the relevant methodology, bereft of any shortcut, to account for all the contributions to the standard binding free energy. Our estimates of the binding affinities support the view that while the Alpha and Beta variants lean on an increased affinity toward the host cell, the Delta and Omicron BA.2 variants choose immune escape. Moreover, the S2E12 antibody, already known to be active against the WT (Starr et al., 2021; Mlcochova et al., 2021), proved to be equally effective against the Delta variant. In stark contrast, H11-D4 retains a low affinity toward the WT compared to that of ACE2 for the latter. Assuming robust structural information, the methodology employed herein successfully addresses the challenging protein-protein binding problem in the context of coronavirus disease 2019 while offering promising perspectives for predictive studies of ever-emerging variants.

Computation of the Binding Energies between Human ACE2 and Spike RBDs of the Original Strain, Delta and Omicron Variants of the SARS-CoV-2: A DFT Simulation Approach

The receptor binding domain (RBD) of SARS-CoV-2 binds to human ACE2 leading to infection. In this study, the complexes that are formed by the attachment of the SARS-CoV-2 spike RBDs of the original strain, delta and omicron variants to the human ACE2 are investigated via density functional theory (DFT) simulations to obtain binding energies. The DFT computations are performed without fragmenting the interfaces to involve longer-range interactions for improved accuracy, which is one of the primary features of the approach used in this study. Basis set superposition error corrections and van der Waals dispersions are also included in the DFT simulations. The binding energies of the SARS-CoV-2 spike RBDs of the original strain, delta and omicron variants to the human ACE2 are computed as -4.76, -6.68, and -11.77 eV, respectively. These binding energy values indicate that the binding of the omicron variant to the ACE2 is much more favorable than the binding of the original strain and the delta variant, which constitute a molecular reason for the takeover of the omicron variant. The binding energies and the decomposition of these energies found in this study are expected to aid in the development of neutralizing agents.