Molecular Dynamics Simulation Study of the Interaction between Human Angiotensin Converting Enzyme 2 and Spike Protein Receptor Binding Domain of the SARS-CoV-2 B.1.617 Variant

Intrinsically Disordered Regions Function as a Cervical Collar to Remotely Regulate the Nodding Dynamics of SARS-CoV-2 Prefusion Spike Heads

The SARS-CoV-2 prefusion spike heads (receptor binding domains, RBDs) frequently nod down and up to interact with host cell receptors. As the spike protein is a trimeric unit of significant size, to understand its large-scale structural dynamics associated with the nodding mechanism and the mutational impact on the same, we develop a topological symmetry-information-loaded coarse-grained structure-based model of a spike trimer using recent cryo-EM structural data. Our study reveals the control of two distant intrinsically disordered regions (IDRs), namely, 630 and FPPR loops, over the nodding dynamics of spike heads. We find that the order-disorder transition of IDRs becomes more evident in the variants of concern (VOCs) that are associated with the characteristic mutation, D614G, in the proximity of these IDRs. In some VOCs, the two other mutations A570D and S982A also show an integral effect. The driver mutation D614G instigates a salt-bridge disruption, altering the order-disorder dynamics of both 630 and FPPR loops and their interaction with the C-terminal domains (CTD1/CTD2). This altered connectivity in these mutants allows the two IDRs to act collectively as a "cervical collar" for the RBD, supporting various spike head postures, consistent with cryo-EM results available for specific cases. The IDRs' control over the spike structure and dynamics presents an exciting opportunity where they can be targeted as remote operational switches to artificially maneuver the nod for effective therapeutic interventions.

Analytical Ultracentrifugation Detects Quaternary Rearrangements and Antibody-Induced Conformational Selection of the SARS-CoV-2 Spike Trimer

Analytical ultracentrifugation (AUC) analysis shows that the SARS-CoV-2 trimeric Spike (S) protein adopts different quaternary conformations in solution. The relative abundance of the "open" and "close" conformations is temperature-dependent, and samples with different storage temperature history have different open/close distributions. Neutralizing antibodies (NAbs) targeting the S receptor binding domain (RBD) do not alter the conformer populations; by contrast, a NAb targeting a cryptic conformational epitope skews the Spike trimer toward an open conformation. The results highlight AUC, which is typically applied for molecular mass determination of biomolecules as a powerful tool for detecting functionally relevant quaternary protein conformations.

Delta variant (B.1.617.2) of SARS-CoV-2: current understanding of infection, transmission, immune escape, and mutational landscape

The Delta variant is one of the alarming variants of the SARS-CoV-2 virus that have been immensely detrimental and a significant cause of the prolonged pandemic (B.1.617.2). During the SARS-CoV-2 pandemic from December 2020 to October 2021, the Delta variant showed global dominance, and afterwards, the Omicron variant showed global dominance. Delta shows high infectivity rate which accounted for nearly 70% of the cases after December 2020. This review discusses the additional attributes that make the Delta variant so infectious and transmissible. The study also focuses on the significant mutations, namely the L452R and T478K present on the receptor-binding domain of spike (S)-glycoprotein, which confers specific alterations to the Delta variant. Considerably, we have also highlighted other notable factors such as the immune escape, infectivity and re-infectivity, vaccine escape, Ro number, S-glycoprotein stability, cleavage pattern, and its binding affinity with the host cell receptor protein. We have also emphasized clinical manifestations, symptomatology, morbidity, and mortality for the Delta variant compared with other significant SARS-CoV-2 variants. This review will help the researchers to get an elucidative view of the Delta variant to adopt some practical strategies to minimize the escalating spread of the SARS-CoV-2 Delta variant.

Entropic Overcompensation of the N501Y Mutation on SARS-CoV-2 S Binding to ACE2

Recent experimental work has shown that the N501Y mutation in the SARS-CoV-2 S glycoprotein's receptor binding domain (RBD) increases binding affinity to the angiotensin-converting enzyme 2 (ACE2), primarily by overcompensating for a less favorable enthalpy of binding by greatly reducing the entropic penalty for complex formation, but the basis for this entropic overcompensation is not clear [Prévost et al. J. Biol. Chem.2021, 297, 101151]. We use all-atom molecular dynamics simulations and free-energy calculations to qualitatively assess the impact of the N501Y mutation on the enthalpy and entropy of binding of RBD to ACE2. Our calculations correctly predict that N501Y causes a less favorable enthalpy of binding to ACE2 relative to the original strain. Furthermore, we show that this is overcompensated for by a more entropically favorable increase in large-scale quaternary flexibility and intraprotein root mean square fluctuations of residue positions upon binding in both RBD and ACE2. The enhanced quaternary flexibility stems from N501Y's ability to remodel the inter-residue interactions between the two proteins away from interactions central to the epitope and toward more peripheral interactions. These findings suggest that an important factor in determining protein-protein binding affinity is the degree to which fluctuations are distributed throughout the complex and that residue mutations that may seem to result in weaker interactions than their wild-type counterparts may yet result in increased binding affinity thanks to their ability to suppress unfavorable entropy changes upon binding.

SARS-CoV-2 Delta Variant: Interplay between Individual Mutations and Their Allosteric Synergy

Since its first appearance in April 2021, B.1.617.2, also termed variant Delta, catalyzed one major worldwide wave dominating the second year of coronavirus disease 2019 (COVID-19) pandemic. Despite its quick disappearance worldwide, the strong virulence caused by a few point mutations remains an unsolved problem largely. Along with the other two sublineages, the Delta variant harbors an accumulation of Spike protein mutations, including the previously identified L452R, E484Q, and the newly emerged T478K on its receptor binding domain (RBD). We used molecular dynamics (MD) simulations, in combination with free energy perturbation (FEP) calculations, to examine the effects of two combinative mutation sets, L452R + E484Q and L452R + T478K. Our dynamic trajectories reveal an enhancement in binding affinity between mutated RBD and the common receptor protein angiotensin converting enzyme 2 (ACE2) through a net increase in the buried molecular surface area of the binary complex. This enhanced binding, mediated through Gln493, sets the same stage for all three sublineages due to the presence of L452R mutation. The other mutation component, E484Q or T478K, was found to impact the RBD-ACE2 binding and help the variant to evade several monoclonal antibodies (mAbs) in a distinct manner. Especially for L452R + T478K, synergies between mutations are mediated through a complex residual and water interaction network and further enhance its binding to ACE2. Taking together, this study demonstrates that new variants of SARS-CoV-2 accomplish both "attack" (infection) and "defense" (antibody neutralization escape) with the same "polished sword" (mutated Spike RBD).

Molecular dynamics study on the strengthening behavior of Delta and Omicron SARS-CoV-2 spike RBD improved receptor-binding affinity

The COVID-19 pandemic caused by a virus that can be transmitted from human to human via air droplets has changed the quality of life and economic systems all over the world. The viral DNA has mutated naturally over time leading to the diversity of coronavirus victims which has posed a serious threat to human security on a massive scale. The current variants have developed in a dominant way and are considered "Variants of Concern" by the World Health Organization (WHO). In this work, Kappa (B.1.617.1), Delta (B.1.617.2), and Omicron (B.1.1.529) variants were obtained to evaluate whether naturally occurring mutations have strengthened viral infectivity. We apply reliable in silico structural dynamics and energetic frameworks of the mutated S-RBD protein for ACE2-binding to analyze and compare the structural information related to the wild-type. In particular, the hotspot residues at Q493, Q498, and N501 on the S-RBD protein were determined as contributing factors to the employment stability of the relevant binding interface. The L452R mutation induces an increment of the hydrogen bonds formed by changing the Q493 environment for ACE2 binding. Moreover, the Q493K exchange in Omicron enables the formation of two additional salt bridges, leading to a strong binding affinity by increased electrostatic interaction energy. These results could be used in proposing concrete informative data for a structure-based design engaged in finding better therapeutics against novel variants.

RBD spatial orientation of the spike protein and its binding to ACE2: insight into the high infectivity of the SARS-CoV-2 Delta variant from MD simulations

The spike glycoprotein on the surface of the SARS-CoV-2 envelope plays an important role in its invasion into host cells. The binding of the spike glycoprotein RBD to the angiotensin-converting enzyme 2 (ACE2) receptor as a critical step in the spread of the virus has been explored intensively since the outbreak of COVID-19, but the high transmissibility of the virus such as the Delta variant is still not fully understood. Here, molecular simulations on the binding interactions of the wild-type spike protein and its four variants (Beta, Kappa, Delta, and Mu) with ACE2 and the antibody were performed, and the present results reveal that the residue mutations will not strengthen the binding affinity of the variant for ACE2, but remarkably influences the spatial orientation of the spike protein. Only the up-right conformational receptor binding domain (RBD) can bind ACE2, which is stabilized by the nearby RBDs in the down state, revealing that the RBD bears dual functional characteristics. The present results provide new insights into plausible mechanisms for high infectivity of the virus variants and their immune escape.

Evolution of Stronger SARS-CoV-2 Variants as Revealed Through the Lens of Molecular Dynamics Simulations

Using molecular dynamics simulations, the protein–protein interactions of the receptor-binding domain of the wild-type and seven variants of the severe acute respiratory syndrome coronavirus 2 spike protein and the peptidase domain of human angiotensin-converting enzyme 2 were investigated. These variants are alpha, beta, gamma, delta, eta, kappa, and omicron. Using 100 ns simulation data, the residue interaction networks at the protein–protein interface were identified. Also, the impact of mutations on essential protein dynamics, backbone flexibility, and interaction energy of the simulated protein–protein complexes were studied. The protein–protein interface for the wild-type, delta, and omicron variants contained several stronger interactions, while the alpha, beta, gamma, eta, and kappa variants exhibited an opposite scenario as evident from the analysis of the inter-residue interaction distances and pair-wise interaction energies. The study reveals that two distinct residue networks at the central and right contact regions forge stronger binding affinity between the protein partners. The study provides a molecular-level insight into how enhanced transmissibility and infectivity by delta and omicron variants are most likely tied to a handful of interacting residues at the binding interface, which could potentially be utilized for future antibody constructs and structure-based antiviral drug design.

Dynamics of the interaction between the receptor-binding domain of SARS-CoV-2 Omicron (B.1.1.529) variant and human angiotensin-converting enzyme 2

Background

The COVID-19 pandemic is still a global public health issue. Omicron, a SARS-CoV-2 B.1.1.529 variant, has raised concerns about transmission and vaccine effectiveness. Omicron currently has the greatest number of variantions.

Methods

To gain a better understanding of the significance of these variations and the dynamics of the interaction between the Omicron spike (S) protein and its human host factor angiotensin-converting enzyme 2 (ACE2), triplicate 500 ns molecular dynamics simulations were run using the structure of the S protein's receptor-binding domain (RBD) in complex with ACE2. The interaction and binding energy, determined using the molecular mechanics-generalized Born surface area approach, were compared to the original SARS-CoV-2 and the B.1.617 variant.

Results

Though mutations K417N and G496S in the S protein RBD disrupt interactions found in the original SARS-CoV-2 complex, mutations Q493R and N501Y introduce interactions not found in the original complex. Interaction at a key viral hotspot and hydrophobic contacts at ACE2's N-terminus were preserved, but intermolecular hydrogen bonds and polar contacts in the S-ACE2 interface were lower than in the original SARS-CoV-2 interface.

Deciphering the Impact of Mutations on the Binding Efficacy of SARS-CoV-2 Omicron and Delta Variants With Human ACE2 Receptor

The pandemic of COVID-19, caused by SARS-CoV-2, has globally affected the human health and economy. Since the emergence of the novel coronavirus SARS-CoV-2, the life-threatening virus continues to mutate and evolve. Irrespective of acquired natural immunity and vaccine-induced immunity, the emerging multiple variants are growing exponentially, crossing the territorial barriers of the modern world. The rapid emergence of SARS-CoV-2 multiple variants challenges global researchers regarding the efficacy of available vaccines and variant transmissibility. SARS-CoV-2 surface-anchored S-protein recognizes and interacts with the host-cell ACE2, facilitating viral adherence and entrance into the cell. Understanding the interfacial interactions between the spike protein of SARS-CoV-2 variants and human ACE2 receptor is important for the design and development of antiviral therapeutics against SARS-CoV-2 emerging variants. Despite extensive research, the crucial determinants related to the molecular interactions between the spike protein of SARS-CoV-2 variants and host receptors are poorly understood. Thus, in this study, we explore the comparative interfacial binding pattern of SARS-CoV-2 spike RBD of wild type, Delta, and Omicron with the human ACE2 receptor to determine the crucial determinants at the atomistic level, using MD simulation and MM/GBSA energy calculations. Based on our findings, the substitution of Q493R, G496S, Q498R, and Y505H induced internal conformational changes in Omicron spike RBD, which leads to higher binding affinity than Delta spike RBD with the human ACE2 receptor, eventually contributing to higher transmission and infectivity. Taken together, these results could be used for the structure-based design of effective antiviral therapeutics against SARS-CoV-2 variants.

Insights into the structure and dynamics of SARS-CoV-2 spike glycoprotein double mutant L452R-E484Q

The Receptor Binding Domain (RBD) of SARS-CoV-2, located on the S1 subunit, plays a vital role in the virus binding and its entry into the host cell through angiotensin-converting enzyme 2 (ACE2) receptor. Therefore, understanding the dynamic effects of mutants on the SARS-CoV-2 RBD is essential for discovering drugs to inhibit the virus binding and disrupt its entry into the host cells. A recent study reported a double mutant of SARS-CoV-2, L452R-E484Q, located in the RBD region. Thus, this study employed various computational algorithms and methods to understand the structural impact of both individual variants L452R, E484Q, and the double mutant L452R-E484Q on the native RBD of spike glycoprotein. The effects of the mutations on native RBD structure were predicted by in silico algorithms, which predicted changes in the protein structure and function upon the mutations. Subsequently, molecular dynamics (MD) simulations were employed to understand the conformational stability and functional changes on the RBD upon the mutations. The comparative results of MD simulation parameters displayed that the double mutant induces significant conformational changes in the spike glycoprotein RBD, which may alter its biological functions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03151-0.

Generalized Methodology for the Quick Prediction of Variant SARS-CoV-2 Spike Protein Binding Affinities with Human Angiotensin-Converting Enzyme II

Variants of the SARS-CoV-2 virus continue to remain a threat 2 years from the beginning of the pandemic. As more variants arise, and the B.1.1.529 (Omicron) variant threatens to create another wave of infections, a method is needed to predict the binding affinity of the spike protein quickly and accurately with human angiotensin-converting enzyme II (ACE2). We present an accurate and convenient energy minimization/molecular mechanics Poisson-Boltzmann surface area methodology previously used with engineered ACE2 therapeutics to predict the binding affinity of the Omicron variant. Without any additional data from the variants discovered after the publication of our first model, the methodology can accurately predict the binding of the spike/ACE2 variant complexes. From this methodology, we predicted that the Omicron variant spike has a Kd of ∼22.69 nM (which is very close to the experimental Kd of 20.63 nM published during the review process of the current report) and that spike protein of the new "Stealth" Omicron variant (BA.2) will display a Kd of ∼12.9 nM with the wild-type ACE2 protein. This methodology can be used with as-yet discovered variants, allowing for quick determinations regarding the variant's infectivity versus either the wild-type virus or its variants.

Mutations in the receptor-binding domain of human SARS CoV-2 spike protein increases its affinity to bind human ACE-2 receptor

The severe acute respiratory syndrome virus-2 (SARS CoV-2) infection has resulted in the current global pandemic. The binding of SARS CoV-2 spike protein receptor-binding domain (RBD) to the human angiotensin converting enzyme-2 (ACE-2) receptor causes the host infection. The spike protein has undergone several mutations with reference to the initial strain isolated during December 2019 from Wuhan, China. A number of these mutant strains have been reported as variants of concern and as variants being monitored. Some of these mutants are known to be responsible for increased transmissibility of the virus. The reason for the increased transmissibility caused by the point mutations can be understood by studying the structural implications and inter-molecular interactions in the binding of viral spike protein RBD and human ACE-2. Here, we use the crystal structure of the RBD in complex with ACE-2 available in the public domain and analyse the 250 ns molecular dynamics (MD) simulations of wild-type and mutants; K417N, K417T, N440K, N501Y, L452R, T478K, E484K and S494P. The ionic, hydrophobic and hydrogen bond interactions, amino acid residue flexibility, binding energies and structural variations are characterized. The MD simulations provide clues to the molecular mechanisms of ACE-2 receptor binding in wild-type and mutant complexes. The mutant spike proteins RBD were associated with greater binding affinity with ACE-2 receptor. Communicated by Ramaswamy H. Sarma.