Understanding the molecular interaction of SARS-CoV-2 spike mutants with ACE2 (angiotensin converting enzyme 2)

Comparative docking and molecular dynamics studies of molnupiravir (EIDD-2801): implications for novel mechanisms of action on influenza and SARS-CoV-2 protein targets

Molnupiravir (EIDD-2801) (MLN) is an oral antiviral drug for COVID-19 treatment, being integrated into viral RNA through RNA-dependent RNA polymerase (RdRp). Upon ingestion, MLN is transformed into two active metabolites: β-d-N4-hydroxycytidine (NHC) (EIDD-1931) in the host plasma, and EIDD-1931-triphosphate (MTP) within the host cells. However, recent studies provide increasing evidence of MLN's interactions with off-target proteins beyond the viral genome, suggesting that the complete mechanisms of action of MLN remain unclear. The aim of this study was therefore to investigate the molecular interactions of MLN in the form of NHC and MTP with the non-RNA structural components of avian influenza (hemagglutinin, neuraminidase) and SARS-CoV-2 (spike glycoprotein, Mpro, and RdRp) viruses and to elucidate whether these two metabolites possess the ability to form stable complexes with these major viral components. Molecular docking of NHC and MTP was performed using AutoDock 4.2.6 and the obtained protein-drug complexes were submitted to 200-ns molecular dynamics simulations in triplicate with subsequent free energy calculations using GROMACS. Docking scores, molecular dynamics and MM/GBSA results showed that MTP was tightly bound within the active site of SARS-CoV-2 RdRp and remained highly stable throughout the 200-ns simulations. Besides, it was also shown that NHC and MTP formed moderately-to-highly stable molecular complexes with off-target receptors hemagglutinin, neuraminidase and Mpro, but rather weak interactions with spike glycoprotein. Our computational findings suggest that NHC and MTP may directly inhibit these receptors, and propose that additional studies on the off-target effects of MLN, i.e. real-time protein binding assays, should be performed.Communicated by Ramaswamy H. Sarma.

Mutations in the Receptor Binding Domain of Severe Acute Respiratory Coronavirus-2 Omicron Variant Spike Protein Significantly Stabilizes Its Conformation

The Omicron variant and its sub-lineages are the only current circulating SARS-CoV-2 viruses worldwide. In this study, the conformational stability of the isolated Receptor Binding Domain (RBD) of Omicron's spike protein is examined in detail. The parent Omicron lineage has over ten mutations in the ACE2 binding region of the RBD that are specifically associated with its β hairpin loop domain. It is demonstrated through biophysical molecular computations that the mutations in the β hairpin loop domain significantly increase the intra-protein interaction energies of intra-loop and loop-RBD interactions. The interaction energy increases include the formation of new hydrogen bonds in the β hairpin loop domain that help stabilize this critical ACE2 binding region. Our results also agree with recent experiments on the stability of Omicron's core β barrel domain, outside of its loop domain, and help demonstrate the overall conformational stability of the Omicron RBD. It is further shown here through dynamic simulations that the unbound state of the Omicron RBD remains closely aligned with the bound state configuration, which was not observed for the wild-type RBD. Overall, these studies demonstrate the significantly increased conformational stability of Omicron over its wild-type configuration and raise a number of questions on whether conformational stability could be a positive selection feature of SARS-CoV-2 viral mutational changes.

Inhibitory effects of SARS-CoV-2 spike protein and BNT162b2 vaccine on erythropoietin-induced globin gene expression in erythroid precursor cells from patients with β-thalassemia

During the recent coronavirus disease 2019 (COVID-19) pandemic several patients with β-thalassemia have been infected by severe acute respiratory syndrome coronavirus (SARS-CoV-2), and most patients were vaccinated against SARS-CoV-2. Recent studies demonstrate an impact of SARS-CoV-2 infection on the hematopoietic system. The main objective of this study was to verify the effects of exposure of erythroid precursor cells (ErPCs) from patients with β-thalassemia to SARS-CoV-2 spike protein (S-protein) and the BNT162b2 vaccine. Erythropoietin (EPO)-cultured ErPCs have been either untreated or treated with S-protein or BNT162b2 vaccine. The employed ErPCs were from a β-thalassemia cellular Biobank developed before the COVID-19 pandemic. The genotypes were β+-IVSI-110/β+-IVSI-110 (one patient),  β039/β+-IVSI-110 (3 patients), and β039/ β039 (2 patients). After treatment with S-protein or BNT162b2 for 5 days, lysates were analyzed by high performance liquid chromatography (HPLC), for hemoglobin production, and isolated RNA was assayed by RT-qPCR, for detection of globin gene expression. The main conclusions of the results obtained are that SARS-CoV-2 S-protein and BNT162b2 vaccine (a) inhibit fetal hemoglobin (HbF) production by β-thalassemic ErPCs and (b) inhibit γ-globin mRNA accumulation. In addition, we have performed in silico studies suggesting a high affinity of S-protein to HbF. Remarkably, the binding interaction energy of fetal hemoglobin to S-protein was comparable with that of angiotensin-converting enzyme 2 (ACE2). Our results are consistent with the hypothesis of a relevant impact of SARS-CoV-2 infection and COVID-19 vaccination on the hematopoietic system.

Chronicling the 3-year evolution of the COVID-19 pandemic: analysis of disease management, characteristics of major variants, and impacts on pathogenicity

Announced on December 31, 2019, the novel coronavirus arising in Wuhan City, Hubei Province resulted in millions of cases and lives lost. Following intense tracking, coronavirus disease 2019 (COVID-19) was declared a pandemic by the World Health Organization (WHO) in 2020. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was identified as the cause of COVID-19 and the continuous evolution of the virus has given rise to several variants. In this review, a comprehensive analysis of the response to the pandemic over the first three-year period is provided, focusing on disease management, development of vaccines and therapeutics, and identification of variants. The transmissibility and pathogenicity of SARS-CoV-2 variants including Alpha, Beta, Gamma, Delta, and Omicron are compared. The binding characteristics of the SARS-CoV-2 spike protein to the angiotensin-converting enzyme 2 (ACE2) receptor and reproduction numbers are evaluated. The effects of major variants on disease severity, hospitalisation, and case-fatality rates are outlined. In addition to the spike protein, open reading frames mutations are investigated. We also compare the pathogenicity of SARS-CoV-2 with SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV). Overall, this study highlights the strengths and weaknesses of the global response to the pandemic, as well as the importance of prevention and preparedness. Monitoring the evolution of SARS-CoV-2 is critical in identifying and potentially predicting the health outcomes of concerning variants as they emerge. The ultimate goal would be a position in which existing vaccines and therapeutics could be adapted to suit new variants in as close to real-time as possible.

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.

The SARS-CoV-2 Hydra, a tiny monster from the 21st century: Thermodynamics of the BA.5.2 and BF.7 variants

SARS-CoV-2 resembles the ancient mythical creature Hydra. Just like with the Hydra, when one head is cut, it is followed by appearance of two more heads, suppression of one SARS-CoV-2 variant causes appearance of newer variants. Unlike Hydra that grows identical heads, newer SARS-CoV-2 variants are usually more infective, which can be observed as time evolution of the virus at hand, which occurs through acquisition of mutations during time. The appearance of new variants is followed by appearance of new COVID-19 pandemic waves. With the appearance of new pandemic waves and determining of sequences, in the scientific community and general public the question is always raised of whether the new variant will be more virulent and more pathogenic. The two variants characterized in this paper, BA.5.2 and BF.7, have caused a pandemic wave during the late 2022. This paper gives full chemical and thermodynamic characterization of the BA.5.2 and BF.7 variants of SARS-CoV-2. Having in mind that Gibbs energy of binding and biosynthesis represent the driving forces for the viral life cycle, based on the calculated thermodynamic properties we can conclude that the newer variants are more infective than earlier ones, but that their pathogenicity has not changed.

Strain wars 5: Gibbs energies of binding of BA.1 through BA.4 variants of SARS-CoV-2

This paper reports, for the first time, standard Gibbs energies of binding of the BA.1, BA.2, BA.3, BA.2.13, BA.2.12.1 and BA.4 Omicron variants of SARS-CoV-2, to the Human ACE2 receptor. Variants BA.1 through BA.3 exhibit a trend of decreasing standard Gibbs energy of binding and hence increased infectivity. The BA.4 variant exhibits a less negative standard Gibbs energy of binding, but also more efficient evasion of the immune response. Therefore, it was concluded that all the analyzed strains evolve in accordance with expectations of the theory of evolution, albeit using different strategies.

Omicron BA.2.75 Sublineage (Centaurus) Follows the Expectations of the Evolution Theory: Less Negative Gibbs Energy of Biosynthesis Indicates Decreased Pathogenicity

SARS-CoV-2 belongs to the group of RNA viruses with a pronounced tendency to mutate. Omicron BA.2.75 is a subvariant believed to be able to suppress the currently dominant BA.5 and cause a new winter wave of the COVID-19 pandemic. Omicron BA.2.75 is characterized by a greater infectivity compared to earlier Omicron variants. However, the Gibbs energy of the biosynthesis of virus particles is slightly less negative compared to those of other variants. Thus, the multiplication rate of Omicron BA.2.75 is lower than that of other SARS-CoV-2 variants. This leads to slower accumulation of newly formed virions and less damage to host cells, indicating evolution of SARS-CoV-2 toward decreasing pathogenicity.

Omicron BA.2.75 Subvariant of SARS-CoV-2 Is Expected to Have the Greatest Infectivity Compared with the Competing BA.2 and BA.5, Due to Most Negative Gibbs Energy of Binding

Omicron BA.2.75 may become the next globally dominant strain of COVID-19 in 2022. The BA.2.75 sub-variant has acquired more mutations (9) in spike protein and other genes of SARS-CoV-2 than any other variant. Thus, its chemical composition and thermodynamic properties have changed compared with earlier variants. In this paper, the Gibbs energy of the binding and antigen-receptor binding rate was reported for the BA.2.75 variant. Gibbs energy of the binding of the Omicron BA.2.75 variant is more negative than that of the competing variants BA.2 and BA.5.

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.

Applications of Molecular Dynamics Simulation in Protein Study

Molecular Dynamics (MD) Simulations is increasingly used as a powerful tool to study protein structure-related questions. Starting from the early simulation study on the photoisomerization in rhodopsin in 1976, MD Simulations has been used to study protein function, protein stability, protein-protein interaction, enzymatic reactions and drug-protein interactions, and membrane proteins. In this review, we provide a brief review for the history of MD Simulations application and the current status of MD Simulations applications in protein studies.

Strain Wars: Competitive interactions between SARS-CoV-2 strains are explained by Gibbs energy of antigen-receptor binding

Since the beginning of the COVID-19 pandemic, SARS-CoV-2 has mutated several times into new strains, with an increased infectivity. Infectivity of SARS-CoV-2 strains depends on binding affinity of the virus to its host cell receptor. In this paper, we quantified the binding affinity using Gibbs energy of binding and analyzed the competition between SARS-CoV-2 strains as an interference phenomenon. Gibbs energies of binding were calculated for several SARS-SoV-2 strains, including Hu-1 (wild type), B.1.1.7 (alpha), B.1.351 (beta), P.1 (Gamma), B.1.36 and B.1.617 (Delta). The least negative Gibbs energy of binding is that of Hu-1 strain, -37.97 kJ/mol. On the other hand, the most negative Gibbs energy of binding is that of the Delta strain, -49.50 kJ/mol. We used the more negative Gibbs energy of binding to explain the increased infectivity of newer SARS-CoV-2 strains compared to the wild type. Gibbs energies of binding was found to decrease chronologically, with appearance of new strains. The ratio of Gibbs energies of binding of mutated strains and wild type was used to define a susceptibility coefficient, which is an indicator of viral interference, where a virus can prevent or partially inhibit infection with another virus.

Decreased Interfacial Dynamics Caused by the N501Y Mutation in the SARS-CoV-2 S1 Spike:ACE2 Complex

Coronavirus Disease of 2019 (COVID-19) caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has resulted in a massive health crisis across the globe, with some genetic variants gaining enhanced infectivity and competitive fitness, and thus significantly aggravating the global health concern. In this regard, the recent SARS-CoV-2 alpha, beta, and gamma variants (B.1.1.7, B.1.351, and P.1 lineages, respectively) are of great significance in that they contain several mutations that increase their transmission rates as evident from clinical reports. By the end of March 2021, these variants were accounting for about two-thirds of SARS-CoV-2 variants circulating worldwide. Specifically, the N501Y mutation in the S1 spike receptor binding domain (S1-RBD) of these variants have been reported to increase its affinity for ACE2, although the basis for this is not entirely clear yet. Here, we dissect the mechanism underlying the increased binding affinity of the N501Y mutant for ACE2 using molecular dynamics (MD) simulations of the available ACE2-S1-RBD complex structure (6M0J) and show a prolonged and stable interfacial interaction of the N501Y mutant S1-RBD with ACE2 compared to the wild type S1-RBD. Additionally, we find that the N501Y mutant S1-RBD displays altered dynamics that likely aids in its enhanced interaction with ACE2. By elucidating a mechanistic basis for the increased affinity of the N501Y mutant S1-RBD for ACE2, we believe that the results presented here will aid in developing therapeutic strategies against SARS-CoV-2 including designing of therapeutic agents targeting the ACE2-S1-RBD interaction.

Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2

Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.

Physicochemical effect of the N501Y, E484K/Q, K417N/T, L452R and T478K mutations on the SARS-CoV-2 spike protein RBD and its influence on agent fitness and on attributes developed by emerging variants of concern

The spike protein comprises one of the main structural components of SARS-CoV-2 because it is directly involved in the infection process and viral transmission, and also because of its immunogenic properties, as an inducer of the protective antibodies production and as a vaccine component. The occurrence of mutations in this region or in other the virus genome regions, comprises a natural phenomenon in its evolution. However, they also occur due to the selective immune pressure, to which the agent is continuously subjected, especially in the spike protein immunodominant regions, such as the RBD. Mutations in the spike protein can change the virus' fitness, increasing its affinity for target cells, its transmissibility and its virulence. In addition, these mutations can giving it the potential ability to evade the protective antibodies action obtained from convalescent sera or vaccine origin, as well as those used in therapy, which may favor the virus expansion and compromise the infection control. Five mutations N501Y, E484K/Q, K417N/T, L452R and T478K, located in the spike protein RBD, have had a greater impact because they are associated with new attributes developed by the virus, which characterize the emerging variants of concern (VOCs) of SARS-Cov-2 identified so far. The occurrence of these mutations induces complex physicochemical effects that can alter the spike protein's structure and its function, which in turn, lead to changes in the agents' fitness. This manuscript discusses the attributes of VOCs associated with the physicochemical effects caused by the aforementioned mutations.

COVID-19 Genetic Variants and Their Potential Impact in Vaccine Development

In the two years since the SARS-CoV-2 pandemic started, it has caused over 5 million deaths and 400 million infected cases, and the world continues to be on high alert for COVID-19. Among the variants of interest and concern of SARS-CoV-2, the current Omicron (B.1.1.529) and stealth Omicron (BA.2) raised serious concerns due to rapid rates of infection caused by numerous mutations in the spike protein, which could escape from the antibody-mediated neutralization and increase the risk of reinfections. Hence, this work aims to describe the most relevant mutations in the SARS-CoV-2 spike protein, discuss vaccine against variant of concerns, describe rare adverse events after COVID-19 vaccination, introduce the most available promising COVID-19 vaccine candidates, and provide few perspectives of the future 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.