Bioinformatics analysis of the differences in the binding profile of the wild-type and mutants of the SARS-CoV-2 spike protein variants with the ACE2 receptor

Impacts of neonicotinoid compounds on the structure and function of Apis mellifera OBP14: Insights from SPR, ITC, multispectroscopy, and molecular modeling

Honeybees are vital for biodiversity and agricultural productivity, yet their populations are declining globally, partly due to exposure to neonicotinoid pesticides. Odorant-binding protein 14 (OBP14) plays an important role in honeybee chemosensation, but its involvement in neonicotinoid toxicity remains underexplored due to limitations in traditional fluorescence spectroscopy techniques. This gap hampers our understanding of neonicotinoid risks to honeybee health. Here, we explored the molecular interactions between OBP14 from Apis mellifera and three widely used neonicotinoids (imidacloprid, thiamethoxam, and clothianidin) using molecular modeling, surface plasmon resonance (SPR), isothermal titration calorimetry (ITC), and multispectroscopy. SPR and ITC characterized the binding affinity, specificity, and thermodynamic parameters of AmelOBP14 interacting with three neonicotinoid compounds, revealing that the binding process is spontaneous and primarily driven by hydrophobic and electrostatic interactions. Molecular modeling highlighted that phenylalanine residue Phe54, near the binding site, plays a critical role in these interactions. UV-vis absorption spectroscopy and synchronous fluorescence spectroscopy (SFS) support slight changes in the microenvironment around the aromatic amino acids of OBP14. Fourier Transform Infrared Spectroscopy (FTIR) and circular dichroism spectroscopy (CD) indicate a decrease in the α-helix content of OBP14, suggesting a change in its secondary structure, while three-dimensional (3D) fluorescence spectroscopy confirms the non-fluorescent nature of the OBP14 polypeptide backbone. The study results revealed its potential as a biomarker for pesticide risk assessment, providing important insights into the molecular mechanisms by which neonicotinoids may impair bee chemosensory function, and offering guidance for the design of safer pesticides to minimize harm to these important pollinators.

Analysis of the bonding affinities between human PSGL-1 and Vpu derived from the different HIV-1 groups - in silico insights

Human P-selectin glycoprotein ligand 1 (PSGL-1) and HIV-1 viral protein U (Vpu) play major roles in limiting and increasing the ability of HIV-1 to infect cells, respectively. There is currently no published data reporting on the specific interactions between PSGL-1 and Vpu, and possible outcomes and consequences of these interactions. To date, it has only been established that Vpu binds human PSGL-1 to degrade PSGL-1 and therefore promote HIV replication. There are, however, four different types of HIV-1, and it would be helpful to know how VpuM, VpuN, VpuO, and VpuP can bind to and possibly inhibit PSGL-1 expression. Bioinformatics methods were used to find out how strongly each type of Vpu found in the different HIV-1 groups bonds with human PSGL-1. Thus, we used molecular docking (MD) and molecular dynamics simulations (MDS) to figure out how PSGL-1 and VpuM, VpuN, VpuO, and VpuP interact with each other. To ensure the reliability of the predicted outcomes, the binding energy of each model was calculated using the MM/GBSA technique. Our findings show that PSGL-1-VpuP (4 H bonds, 2 salt bridges) and PSGL-1-VpuM (3 H bonds, 2 salt bridges) have stronger bonding affinities than PSGL-1-VpuN (4 H bonds, no salt bridges) and PSGL-1-VpuO (2 H bonds, 1 salt bridge). The MDS test also shows that PSGL-1-VpuM and PSGL-1-VpuP protein complexes are more stable and compact, with lower residual fluctuations compared to PSGL-1-VpuN and PSGL-1-VpuO protein complexes. Binding free energies of -82.27 ± 1.35 kcal/mol, -82.17 ± 0.84 kcal/mol, -67.84 ± 0.63 kcal/mol, and -131.86 ± 1.08 kcal/mol were recorded for each of PSGL1-VpuM, PSGL1-VpuN, PSGL1-VpuO, and PSGL1-VpuP, respectively, which further supports our results. Our research shows that Vpu from the M and P HIV-1 groups may be better at blocking human PSGL-1 than VpuO and VpuN groups. These results are novel in this specific realm of HIV research, and as such, further investigations in more robust experimental studies are warranted.

Medicinal Phytocompounds as Potential Inhibitors of p300-HIF1α Interaction: A Structure-Based Screening and Molecular Dynamics Simulation Study

Background: Hypoxia plays a key role in cancer progression, mainly by stabilizing and activating hypoxia-inducible factor-1 (HIF-1). For HIF-1 to function under low oxygen conditions, it must interact with the transcriptional coactivator p300, a critical step for promoting cancer cell survival and adaptation in hypoxic environments. Methods: Consequently, we used drug design and molecular simulation techniques to screen phytochemical databases, including traditional Chinese and African medicine sources, for compounds that could disrupt the p300/HIF-1 interaction. Results: In this study, we identified potential compounds with high docking scores such as EA-176920 (-8.719), EA-46881231 (-8.642), SA-31161 (-9.580), SA-5280863 (-8.179), NE-5280362 (-10.287), NE-72276 (-9.017), NA-11210533 (-10.366), NA-11336960 (-7.818), TCM-5281792 (-12.648), and TCM-6441280 (-9.470 kcal/mol) as lead compounds. Furthermore, the compound with the highest docking score from each database (EA-176920, SA-31161, NE-5280362, NA-11210533, and TCM-5281792) was subjected to further analysis. The stable binding affinity of these compounds with p300 was confirmed by Post-simulation binding free energy (-22.0020 kcal/mol, -25.4499 kcal/mol, -32.4530 kcal/mol, -33.9918 kcal/mol, and -57.7755 kcal/mol, respectively) and KD analysis. Moreover, the selected compounds followed the Lipinski rules with favorable ADMET properties like efficient intestinal absorption, high water solubility, and no toxicity. Conclusions: Our findings highlight the potential of natural compounds to target key protein-protein interactions in cancer and lay the groundwork for future in vitro and in vivo studies to explore their therapeutic potential. Specifically, disrupting the p300/HIF-1 interaction could interfere with hypoxia-driven pathways that promote tumor growth, angiogenesis, and metastasis, offering a promising strategy to suppress cancer progression at the molecular level.

Screening of medicinal phytocompounds with structure-based approaches to target key hotspot residues in tyrosyl-DNA phosphodiesterase 1: augmenting sensitivity of cancer cells to topoisomerase I inhibitors

One of cancer's well-known hallmarks is DNA damage, yet it's intriguing that DNA damage has been explored as a therapeutic strategy against cancer. Tyrosyl-DNA phosphodiesterase 1, involved in DNA repair from topoisomerase I inhibitors, a chemotherapy class for cancer treatment. Inhibiting TDP1 can increase unresolved Top1 cleavage complexes in cancer cells, inducing DNA damage and cell death. TDP1's catalytic activity depends on His263 and His493 residues. Using molecular simulation, structure-based drug design, and free energy calculation, we identified potential drugs against TDP1. A multi-step screening of medicinal plant compound databases (North Africa, East Africa, Northeast Africa, and South Africa) identified the top four candidates. Docking scores for top hits 1-4 were -7.76, -7.37, -7.35, and -7.24 kcal/mol. Top hit 3 exhibited the highest potency, forming a strong bonding network with both His263 and His493 residues. All-atoms simulations showed consistent dynamics for top hits 1-4, indicating stability and potential for efficient interaction with interface residues. Minimal fluctuations in residue flexibility suggest these compounds can stabilize internal flexibility upon binding. The binding free energies of -35.11, -36.70, -31.38, and -23.85 kcal/mol were calculated for the top hit 1-4 complexes. Furthermore, the chosen compounds demonstrate outstanding ADMET characteristics, such as excellent water solubility, effective gastrointestinal absorption, and the absence of hepatotoxicity. Cytotoxicity analysis revealed top hit 2 higher probability of activity against 24 cancer cell lines. Our findings suggest that these compounds (top hits 1-4) hold promise for innovative drug therapies, suitable for both in vivo and in vitro experiments.

Exploring global natural product databases for NLRP3 inhibition: Unveiling novel combinatorial therapeutic strategy for hidradenitis suppurativa

BACKGROUND

Hidradenitis suppurativa (HS) is a chronic inflammatory skin condition of the terminal hair follicle, which can present in sporadic, familial, or syndromic forms. The exact pathogenesis of HS remains elusive, posing a challenge for the development of effective treatments. Among the various immunological mechanisms, the NLRP3 inflammasome is thought to contribute to the pathogenesis of HS, releasing cytokines such as IL-1β and IL-18 which initiates and exacerbates inflammation. Consequently, targeting NLRP3 offers a potential strategy for mitigating inflammation in HS-affected skin.

METHODS

In this study we used the docking, molecular dynamics simulation and binding free energy approaches to identify the potent inhibitor of NLRP3 by screening the African phytocompounds and traditional Chinese medicine databases.

RESULTS

Our virtual drug screening analysis identified two lead compounds from each database, characterized by high docking scores such as SA-21676268 (-8.135 kcal/mol), SA-167673 (-10.251 kcal/mol), EA-45360194 (-10.376 kcal/mol), EA-46881231 (-10.011 kcal/mol), NEA-44258150 (-9.856 kcal/mol), NEA-135926572 (-7.662 kcal/mol), NA-163089376 (-9.237 kcal/mol), NA-440735 (-8.826 kcal/mol), TCM-392442 (-10.438 kcal/mol), and TCM-10043097 (-9.046 kcal/mol) which highlighted the strong binding affinity as compared to the control NP3-146 drug (-5.09 kcal/mol). Moreover, the values of dissociation constant further validated the strong binding affinity between the identified lead compounds and NLRP3. The dynamic stability and strong bonding energies of the lead compounds-NLRP3 complexes were confirmed by the molecular dynamic simulation and binding free energy calculation. The analysis of ADMET properties for all compounds indicated high intestinal absorption, water solubility, absence of hepatotoxicity, and skin sensitivity.

CONCLUSION

In conclusion, our molecular simulations and binding free energy calculations confirmed the strong affinity of these lead compounds for NLRP3 as compared to the control drug, highlighting their potential as part of a combinatorial therapeutic strategy for HS to effectively reduce disease-related inflammation.

The apoptotic and anti-proliferative effects of Neosetophomone B in T-cell acute lymphoblastic leukaemia via PI3K/AKT/mTOR pathway inhibition

The phosphatidylinositol 3-kinase/Protein Kinase B/mammalian target of rapamycin (PI3K/AKT/mTOR) signalling pathway is pivotal in various cancers, including T-cell acute lymphoblastic leukaemia (T-ALL), a particularly aggressive type of leukaemia. This study investigates the effects of Neosetophomone B (NSP-B), a meroterpenoid fungal metabolite, on T-ALL cell lines, focusing on its anti-cancer mechanisms and therapeutic potential. NSP-B significantly inhibited the proliferation of T-ALL cells by inducing G0/G1 cell cycle arrest and promoting caspase-dependent apoptosis. Additionally, NSP-B led to the dephosphorylation and subsequent inactivation of the PI3K/AKT/mTOR signalling pathway, a critical pathway in cell survival and growth. Molecular docking studies revealed a strong binding affinity of NSP-B to the active site of AKT, primarily involving key residues crucial for its activity. Interestingly, NSP-B treatment also induced apoptosis and significantly reduced proliferation in phytohemagglutinin-activated primary human CD3+ T cells, accompanied by a G0/G1 cell cycle arrest. Importantly, NSP-B did not affect normal primary T cells, indicating a degree of selectivity in its action, targeting only T-ALL cells and activated T cells. In conclusion, our findings highlight the potential of NSP-B as a novel therapeutic agent for T-ALL, specifically targeting the aberrantly activated PI3K/AKT/mTOR pathway and being selective in action. These results provide a strong basis for further investigation into NSP-B's anti-cancer properties and potential application in T-ALL clinical therapies.

Immunoinformatic based designing of highly immunogenic multi-epitope subunit vaccines to stimulate an adaptive immune response against Junin virus

The Junin virus causes Argentine hemorrhagic fever, leading to severe complications such as high fever, malaise, muscle pain, and bleeding disorders, including hemorrhages in the skin and mucous membranes. Neurological issues like confusion, seizures, and coma can also occur. Without prompt and effective treatment, the disease can be fatal, with mortality rates reaching up to 30%. Taking serious measures is essential to mitigate the spread of the disease. Vaccination is the most effective choice to neutralize the Junin virus in the current situation. Consequently, to design the highly immunogenic and non-allergenic multi-epitope subunit vaccine against the Junin virus, we employed the immunoinformatic approach to screen the glycoprotein, nucleoprotein, and RDRP protein for potential immunogenic CTL (Cytotoxic T Lymphocyte), HTL (Helper T Lymphocyte) and B (B Lymphocyte) cell epitopes. Afterward, the predicted epitopes were subjected to 3D modeling and validation. The strong binding affinity of the constructed vaccines with the human TLR3 was confirmed through molecular docking, with scores of - 333 kcal/mol for glycoprotein, - 297 kcal/mol for nucleoprotein, - 308 kcal/mol for RDRP, and - 305 kcal/mol for combined vaccines. Additionally, the binding free energies recorded were - 63.54 kcal/mol, - 64.16 kcal/mol, - 56.81 kcal/mol, and - 51.52 kcal/mol, respectively. Furthermore, the dynamic stability, residual fluctuation, and compactness of vaccine-TLR-3 complexes were confirmed by the molecular dynamic simulation. The codon adaptation index (CAI) values and high GC content confirmed the stable expression of constructed vaccines in the pET-28a ( +) expression vector. The immune simulation analysis demonstrated that administering booster doses of the developed vaccines resulted in a notable increase in IgG, IgM, interleukins, and cytokines levels, indicating effective antigen clearance over time. In conclusion, our study provides preclinical evidence for designing a highly effective Junin virus vaccine, necessitating further in-vitro and in-vivo experiments.

Natural Compounds Targeting Thymic Stromal Lymphopoietin (TSLP): A Promising Therapeutic Strategy for Atopic Dermatitis

Atopic dermatitis (AD) is a chronic inflammatory skin disease with rising prevalence, marked by eczematous lesions, itching, and a weakened skin barrier often tied to filaggrin gene mutations. This breakdown allows allergen and microbe entry, with thymic stromal lymphopoietin (TSLP) playing a crucial role by activating immune pathways that amplify the allergic response. TSLP's central role in AD pathogenesis makes it a promising therapeutic target. Consequently, in this study, we used the virtual drug screening, molecular dynamics simulation, and binding free energies calculation approaches to explore the African Natural Product Database against the TSLP protein. The molecular screening identified four compounds with high docking scores, namely SA_0090 (-7.37), EA_0131 (-7.10), NA_0018 (-7.03), and WA_0006 (-6.99 kcal/mol). Furthermore, the KD analysis showed a strong binding affinity of these compounds with TSLP, with values of -5.36, -5.36, -5.34, and -5.32 kcal/mol, respectively. Moreover, the strong binding affinity of these compounds was further validated by molecular dynamic simulation analysis, which revealed that the WA_0006-TSLP is the most stable complex with the lowest average RMSD. However, the total binding free energies were -40.5602, -41.0967, -27.3293, and -51.3496 kcal/mol, respectively, showing the strong interaction between the selected compounds and TSLP. Likewise, these compounds showed excellent pharmacokinetics characteristics. In conclusion, this integrative approach provides a foundation for the development of safe and effective treatments for AD, potentially offering relief to millions of patients worldwide.

Discovery of novel MLK4 inhibitors against colorectal cancer through computational approaches

Colorectal cancer (CRC) is a significant health issue globally, affecting approximately 10 % of the world's population. The prevalence of CRC highlights the need for effective treatments and prevention strategies. The current therapeutic option, such as chemotherapy, has significant side effects. Thus, this study investigated the anticancer properties of Sanguinarine derivatives, an alkaloid found in traditional herbs via chemoinformatic approaches. Six Sanguinarine derivatives were discovered through virtual screening and molecular docking to determine their binding affinities against the mixed lineage kinase (MLK4) protein which is responsible for CRC. All the compounds were found to be more effective than standard drug used for colorectal cancer treatment, with Sanguinarine derivative 11 showing the highest affinity. The stability of the drug was confirmed through molecular dynamics simulations at 500 ns. This suggests that compound 11 has a higher chance of replacing 5-Fluorouracil, which is currently a widely used chemotherapy drug. Before molecular dynamics simulations, the pharmacokinetic and chemical properties of Sanguinarine derivatives were determined using pkCSM server and DFT method, respectively. The results support that compound 11 is a good drug candidate, as evidenced by Lipinski's Rule of Five. Therefore, compound 11 is recommended for further analysis via in vivo and in vitro studies to confirm its efficacy and safety.

Examining sialic acid derivatives as potential inhibitors of SARS-CoV-2 spike protein receptor binding domain

Severe acute respiratory syndrome coronavirus 2 (SARS CoV-2) has been the primary reason behind the COVID-19 global pandemic which has affected millions of lives worldwide. The fundamental cause of the infection is the molecular binding of the viral spike protein receptor binding domain (SP-RBD) with the human cell angiotensin-converting enzyme 2 (ACE2) receptor. The infection can be prevented if the binding of RBD-ACE2 is resisted by utilizing certain inhibitors or drugs that demonstrate strong binding affinity towards the SP RBD. Sialic acid based glycans found widely in human cells and tissues have notable propensity of binding to viral proteins of the coronaviridae family. Recent experimental literature have used N-acetyl neuraminic acid (Sialic acid) to create diagnostic sensors for SARS-CoV-2, but a detailed interrogation of the underlying molecular mechanisms is warranted. Here, we perform all atom molecular dynamics (MD) simulations for the complexes of certain Sialic acid-based molecules with that of SP RBD of SARS CoV-2. Our results indicate that Sialic acid not only reproduces a binding affinity comparable to the RBD-ACE2 interactions, it also assumes the longest time to dissociate completely from the protein binding pocket of SP RBD. Our predictions corroborate that a combination of electrostatic and van der Waals energies as well the polar hydrogen bond interactions between the RBD residues and the inhibitors influence free energy of binding.Communicated by Ramaswamy H. Sarma.

Molecular screening of phytocompounds targeting the interface between influenza A NS1 and TRIM25 to enhance host immune responses

BACKGROUND

Influenza A virus causes severe respiratory illnesses, especially in developing nations where most child deaths under 5 occur due to lower respiratory tract infections. The RIG-I protein acts as a sensor for viral dsRNA, triggering interferon production through K63-linked poly-ubiquitin chains synthesized by TRIM25. However, the influenza A virus's NS1 protein hinders this process by binding to TRIM25, disrupting its association with RIG-I and preventing downstream interferon signalling, contributing to the virus's evasion of the immune response.

METHODS

In our study we used structural-based drug designing, molecular simulation, and binding free energy approaches to identify the potent phytocompounds from various natural product databases (>100,000 compounds) able to inhibit the binding of NS1 with the TRIM25.

RESULTS

The molecular screening identified EA-8411902 and EA-19951545 from East African Natural Products Database, NA-390261 and NA-71 from North African Natural Products Database, SA-65230 and SA- 4477104 from South African Natural Compounds Database, NEA- 361 and NEA- 4524784 from North-East African Natural Products Database, TCM-4444713 and TCM-6056 from Traditional Chinese Medicines Database as top hits. The molecular docking and binding free energies results revealed that these compounds have high affinity with the specific active site residues (Leu95, Ser99, and Tyr89) involved in the interaction with TRIM25. Additionally, analysis of structural dynamics, binding free energy, and dissociation constants demonstrates a notably stronger binding affinity of these compounds with the NS1 protein. Moreover, all selected compounds exhibit exceptional ADMET properties, including high water solubility, gastrointestinal absorption, and an absence of hepatotoxicity, while adhering to Lipinski's rule.

CONCLUSION

Our molecular simulation findings highlight that the identified compounds demonstrate high affinity for specific active site residues involved in the NS1-TRIM25 interaction, exhibit exceptional ADMET properties, and adhere to drug-likeness criteria, thus presenting promising candidates for further development as antiviral agents against influenza A virus infections.

Abrogation of ORF8-IRF3 binding interface with Carbon nanotube derivatives to rescue the host immune system against SARS-CoV-2 by using molecular screening and simulation approaches

The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has led to over six million deaths worldwide. In human immune system, the type 1 interferon (IFN) pathway plays a crucial role in fighting viral infections. However, the ORF8 protein of the virus evade the immune system by interacting with IRF3, hindering its nuclear translocation and consequently downregulate the type I IFN signaling pathway. To block the binding of ORF8-IRF3 and inhibit viral pathogenesis a quick discovery of an inhibitor molecule is needed. Therefore, in the present study, the interface between the ORF8 and IRF3 was targeted on a high-affinity carbon nanotube by using computational tools. After analysis of 62 carbon nanotubes by multiple docking with the induced fit model, the top five compounds with high docking scores of - 7.94 kcal/mol, - 7.92 kcal/mol, - 7.28 kcal/mol, - 7.19 kcal/mol and - 7.09 kcal/mol (top hit1-5) were found to have inhibitory activity against the ORF8-IRF3 complex. Molecular dynamics analysis of the complexes revealed the high compactness of residues, stable binding, and strong hydrogen binding network among the ORF8-nanotubes complexes. Moreover, the total binding free energy for top hit1-5 was calculated to be - 43.21 ± 0.90 kcal/mol, - 41.17 ± 0.99 kcal/mol, - 48.85 ± 0.62 kcal/mol, - 43.49 ± 0.77 kcal/mol, and - 31.18 ± 0.78 kcal/mol respectively. These results strongly suggest that the identified top five nanotubes (hit1-5) possess significant potential for advancing and exploring innovative drug therapies. This underscores their suitability for subsequent in vivo and in vitro experiments, marking them as promising candidates worthy of further investigation.

Tackling suppressive cancer microenvironment by NARF-derived immune modulatory vaccine and its validation using simulation strategies

Introduction: Targeting tumor microenvironment is beneficial and present an ideal setting for the development of futuristic immunotherapy. Here, we make use of Nuclear prelamin A recognition factor (NARF), a protein linked to the coactivation of transcriptional regulators in human breast cancer stem cells (CSC) in our investigation.

Methods: In this study, we initially computed the epitope regions possessing the ability to stimulate both T and B cells within the NARF protein. These identified epitope areas were fused with an adjuvant such as RpfB and RpfE as well as linkers like AAY, GPGPG, KK, and EAAAK. The constructed vaccine was further characterized by assessing its physicochemical properties and population coverage. The potential interactions of the designed vaccine with different toll-like receptors were examined by a sequence of computational studies. Of note, docking study were employed to understand its mechanism of action. Molecular dynamics and immune simulation studies were conducted to comprehend more into their structural stability and immune responses. The resultant vaccine was back-translated, codon-optimised and introduced into pET-28 (+) vector.

Results and discussion: We hypothesize from the results that the designed NARF protein-based vaccine in our analysis could effectively provoke the immune responses in the target organism through TLR-7 binding and promotes MHC class-II mediated antigen presentation. Indeed, comprehensive evaluations conducted in both in vitro and in vivo settings are imperative to substantiate the safety and efficacy of the developed vaccine.

Elucidating the binding mechanism of SARS-CoV-2 NSP6-TBK1 and structure-based designing of phytocompounds inhibitors for instigating the host immune response

SARS-CoV-2, also referred to as severe acute respiratory syndrome coronavirus 2, is the virus responsible for causing COVID-19, an infectious disease that emerged in Wuhan, China, in December 2019. Among its crucial functions, NSP6 plays a vital role in evading the human immune system by directly interacting with a receptor called TANK-binding kinase (TBK1), leading to the suppression of IFNβ production. Consequently, in the present study we used the structural and biophysical approaches to analyze the effect of newly emerged mutations on the binding of NSP6 and TBK1. Among the identified mutations, four (F35G, L37F, L125F, and I162T) were found to significantly destabilize the structure of NSP6. Furthermore, the molecular docking analysis highlighted that the mutant NSP6 displayed its highest binding affinity with TBK1, exhibiting docking scores of -1436.2 for the wildtype and -1723.2, -1788.6, -1510.2, and -1551.7 for the F35G, L37F, L125F, and I162T mutants, respectively. This suggests the potential for an enhanced immune system evasion capability of NSP6. Particularly, the F35G mutation exhibited the strongest binding affinity, supported by a calculated binding free energy of -172.19 kcal/mol. To disrupt the binding between NSP6 and TBK1, we conducted virtual drug screening to develop a novel inhibitor derived from natural products. From this screening, we identified the top 5 hit compounds as the most promising candidates with a docking score of -6.59 kcal/mol, -6.52 kcal/mol, -6.32 kcal/mol, -6.22 kcal/mol, and -6.21 kcal/mol. The molecular dynamic simulation of top 3 hits further verified the dynamic stability of drugs-NSP6 complexes. In conclusion, this study provides valuable insight into the higher infectivity of the SARS-CoV-2 new variants and a strong rationale for the development of novel drugs against NSP6.

Exploring the natural products chemical space to abrogate the F3L-dsRNA interface of monkeypox virus to enhance the immune responses using molecular screening and free energy calculations

Amid the ongoing monkeypox outbreak, there is an urgent need for the rapid development of effective therapeutic interventions capable of countering the immune evasion mechanisms employed by the monkeypox virus (MPXV). The evasion strategy involves the binding of the F3L protein to dsRNA, resulting in diminished interferon (IFN) production. Consequently, our current research focuses on utilizing virtual drug screening techniques to target the RNA binding domain of the F3L protein. Out of the 954 compounds within the South African natural compound database, only four demonstrated notable docking scores: -6.55, -6.47, -6.37, and -6.35 kcal/mol. The dissociation constant (KD) analysis revealed a stronger binding affinity of the top hits 1-4 (-5.34, -5.32, -5.29, and -5.36 kcal/mol) with the F3L in the MPXV. All-atom simulations of the top-ranked hits 1 to 4 consistently exhibited stable dynamics, suggesting their potential to interact effectively with interface residues. This was further substantiated through analyses of parameters such as radius of gyration (Rg), Root Mean Square Fluctuation, and hydrogen bonding. Cumulative assessments of binding free energy confirmed the top-performing candidates among all the compounds, with values of -35.90, -52.74, -28.17, and -32.11 kcal/mol for top hits 1-4, respectively. These results indicate that compounds top hit 1-4 could hold significant promise for advancing innovative drug therapies, suggesting their suitability for both in vivo and in vitro experiments.

Multi-omic analysis characterizes molecular susceptibility of receptors to SARS-CoV-2 spike protein

In the post COVID-19 era, new SARS-CoV-2 variant strains may continue emerging and long COVID is poised to be another public health challenge. Deciphering the molecular susceptibility of receptors to SARS-CoV-2 spike protein is critical for understanding the immune responses in COVID-19 and the rationale of multi-organ injuries. Currently, such systematic exploration remains limited. Here, we conduct multi-omic analysis of protein binding affinities, transcriptomic expressions, and single-cell atlases to characterize the molecular susceptibility of receptors to SARS-CoV-2 spike protein. Initial affinity analysis explains the domination of delta and omicron variants and demonstrates the strongest affinities between BSG (CD147) receptor and most variants. Further transcriptomic data analysis on 4100 experimental samples and single-cell atlases of 1.4 million cells suggest the potential involvement of BSG in multi-organ injuries and long COVID, and explain the high prevalence of COVID-19 in elders as well as the different risks for patients with underlying diseases. Correlation analysis validated moderate associations between BSG and viral RNA abundance in multiple cell types. Moreover, similar patterns were observed in primates and validated in proteomic expressions. Overall, our findings implicate important therapeutic targets for the development of receptor-specific vaccines and drugs for COVID-19.

Cupressus sempervirens L. flavonoids as potent inhibitors to xanthine oxidase: in vitro, molecular docking, ADMET and PASS studies

Excessive intake of purine-rich foods such as seafood and red meat leads to excess xanthine oxidase activity and provokes gout attacks. The aim of this paper is to evaluate in vitro and in silico, the inhibition effect of Cupressus sempervirens plant extracts (flavonoids (Cae) and alkaloids (CaK)) and its six derivative compounds on bovine xanthine oxidase (BXO). The in silico study consists of molecular docking with GOLD v4.0 based on the best PLPchem score (PLP) and prediction of biological activity with the PASS server tool. The inhibitors used were lignan (cp1), Amentoflavone (cp2), Cupressuflavone (cp3), Isocryptomerin (cp4), Hinokiflavone (cp5), and Neolignan (cp6). The in vitro results showed that CaK gives an IC50 of 3.52 ± 0.04 μg/ml. Similarly, Cae saved an IC50 of 8.46 ± 1.98 μg/ml compared with the control (2.82 ± 0.10 μg/ml). The in silico results show that cp1 was the best inhibitor model (PLP of 88.09) with approved pharmacokinetics. These findings suggest that cp1 and cp2 may offer good alternatives for the treatment of hyperuricemia; cp3 was moderate, while the others (cp4 to cp6) were considered weak inhibitors according to their PLP.Communicated by Ramaswamy H. Sarma.

How helpful were molecular dynamics simulations in shaping our understanding of SARS-CoV-2 spike protein dynamics?

The SARS-CoV-2 spike protein (S) represents an important viral component that is required for successful viral infection in humans owing to its essential role in recognition of and entry to host cells. The spike is also an appealing target for drug designers who develop vaccines and antivirals. This article is important as it summarizes how molecular simulations successfully shaped our understanding of spike conformational behavior and its role in viral infection. MD simulations found that the higher affinity of SARS-CoV-2-S to ACE2 is linked to its unique residues that add extra electrostatic and van der Waal interactions in comparison to the SARS-CoV S. This illustrates the spread potential of the pandemic SARS-CoV-2 relative to the epidemic SARS-CoV. Different mutations at the S-ACE2 interface, which is believed to increase the transmission of the new variants, affected the behavior and binding interactions in different simulations. The contributions of glycans to the opening of S were revealed via simulations. The immune evasion of S was linked to the spatial distribution of glycans. This help the virus to escape the immune system recognition. This article is important as it summarizes how molecular simulations successfully shaped our understanding of spike conformational behavior and its role in viral infection. This will pave the way to us preparing for the next pandemic as the computational tools are tailored to help fight new challenges.

The XBB.1.5 slightly increase the binding affinity for host receptor ACE2 and exhibit strongest immune escaping features: molecular modeling and free energy calculation

Introduction: The current XBB variant of SARS-CoV-2 with the strongest immune escaping properties is currently the most dominant variant circulating around the world. With the emergence of XBB global morbidities and mortalities have raised again. In the current scenario, it was highly required to delineate the binding capabilities of NTD of XBB subvariant towards human neutralizing antibodies and to dig out the binding affinity of RBD of XBB subvariant with ACE2 receptor. Materials and Methods: The current study uses molecular interaction and simulation-based approaches to decipher the binding mechanism of RBD with ACE2 and mAb interaction with NTD of the spike protein. Results: Molecular docking of the Wild type NTD with mAb revealed a docking score of -113.2 ± 0.7 kcal/mol while XBB NTD docking with mAb reported -76.2 ± 2.3 kcal/mol. On the other hand, wild-type RBD and XBB RBD with ACE2 receptor demonstrated docking scores of -115.0 ± 1.5 kcal/mol and -120.8 ± 3.4 kcal/mol respectively. Moreover, the interaction network analysis also revealed significant variations in the number of hydrogen bonds, salt-bridges, and non-bonded contacts. These findings were further validated by computing the dissociation constant (KD). Molecular simulation analysis such as RMSD, RMSF, Rg and hydrogen bonding analysis revealed variation in the dynamics features of the RBD and NTD complexes due to the acquired mutations. Furthermore, the total binding energy for the wild-type RBD in complex with ACE2 reported -50.10 kcal/mol while XBB-RBD coupled with ACE2 reported -52.66 kcal/mol respectively. This shows though the binding of XBB is slightly increased but due to the variation in the bonding network and other factors makes the XBB variant to enter into the host cell efficiently than the wild type. On the other hand, the total binding free energy for the wildtype NTD-mAb was calculated to be -65.94 kcal/mol while for XBB NTD-mAb was reported to be -35.06 kcal/mol respectively. The significant difference in the total binding energy factors explains that the XBB variant possess stronger immune evasion properties than the others variants and wild type. Conclusions: The current study provides structural features for the XBB variant binding and immune evasion which can be used to design novel therapeutics.

PTRH2 Gene Variants: Recent Review of the Phenotypic Features and Their Bioinformatics Analysis

Peptidyl-tRNA hydrolase 2 (PTRH2) is an evolutionarily highly conserved mitochondrial protein. The biallelic mutations in the PTRH2 gene have been suggested to cause a rare autosomal recessive disorder characterized by an infantile-onset multisystem neurologic endocrine and pancreatic disease (IMNEPD). Patients with IMNEPD present varying clinical manifestations, including global developmental delay associated with microcephaly, growth retardation, progressive ataxia, distal muscle weakness with ankle contractures, demyelinating sensorimotor neuropathy, sensorineural hearing loss, and abnormalities of thyroid, pancreas, and liver. In the current study, we conducted an extensive literature review with an emphasis on the variable clinical spectrum and genotypes in patients. Additionally, we reported on a new case with a previously documented mutation. A bioinformatics analysis of the various PTRH2 gene variants was also carried out from a structural perspective. It appears that the most common clinical characteristics among all patients include motor delay (92%), neuropathy (90%), distal weakness (86.4%), intellectual disability (84%), hearing impairment (80%), ataxia (79%), and deformity of head and face (~70%). The less common characteristics include hand deformity (64%), cerebellar atrophy/hypoplasia (47%), and pancreatic abnormality (35%), while the least common appear to be diabetes mellitus (~30%), liver abnormality (~22%), and hypothyroidism (16%). Three missense mutations were revealed in the PTRH2 gene, the most common one being Q85P, which was shared by four different Arab communities and was presented in our new case. Moreover, four different nonsense mutations in the PTRH2 gene were detected. It may be concluded that disease severity depends on the PTRH2 gene variant, as most of the clinical features are manifested by nonsense mutations, while only the common features are presented by missense mutations. A bioinformatics analysis of the various PTRH2 gene variants also suggested the mutations to be deleterious, as they seem to disrupt the structural confirmation of the enzyme, leading to loss of stability and functionality.

Structural insights into the effect of mutations in the spike protein of SARS-CoV-2 on the binding with human furin protein

The SARS COV-2 and its variants are spreading around the world at an alarming speed, due to its higher transmissibility and the conformational changes caused by mutations. The resulting COVID-19 pandemic has imposed severe health consequences on human health. Several countries of the world including Pakistan have studied its genome extensively and provided productive findings. In the current study, the mCSM, DynaMut2, and I-Mutant servers were used to analyze the effect of identified mutations on the structural stability of spike protein however, the molecular docking and simulations approaches were used to evaluate the dynamics of the bonding network between the wild-type and mutant spike proteins with furin. We addressed the mutational modifications that have occurred in the spike protein of SARS-COV-2 that were found in 215 Pakistani's isolates of COVID-19 patients to study the influence of mutations on the stability of the protein and its interaction with the host cell. We found 7 single amino acid substitute mutations in various domains that reside in spike protein. The H49Y, N74K, G181V, and G446V were found in the S1 domain while the D614A, V622F, and Q677H mutations were found in the central helices of the spike protein. Based on the observation, G181V, G446V, D614A, and V622F mutants were found highly destabilizing and responsible for structural perturbation. Protein-protein docking and molecular simulation analysis with that of furin have predicted that all the mutants enhanced the binding efficiency however, the V622F mutant has greatly altered the binding capacity which is further verified by the K_D value (7.1 E^-14) and therefore may enhance the spike protein cleavage by Furin and increase the rate of infectivity by SARS-CoV-2. On the other hand, the total binding energy for each complex was calculated which revealed -50.57 kcal/mol for the wild type, for G181V -52.69 kcal/mol, for G446V -56.44 kcal/mol, for D614A -59.78 kcal/mol while for V622F the TBE was calculated to be -85.84 kcal/mol. Overall, the current finding shows that these mutations have increased the binding of Furin for spike protein and shows that D614A and V622F have significant effects on the binding and infectivity.

Genomic and structural mechanistic insight to reveal the differential infectivity of omicron and other variants of concern

BACKGROUND

The genome of SARS-CoV-2, is mutating rapidly and continuously challenging the management and preventive measures adopted and recommended by healthcare agencies. The spike protein is the main antigenic site that binds to the host receptor hACE-2 and is recognised by antibodies. Hence, the mutations in this site were analysed to assess their role in differential infectivity of lineages having these mutations, rendering the characterisation of these lineages as variants of concern (VOC) and variants of interest (VOI).

METHODS

In this work, we examined the genome sequence of SARS-CoV-2 VOCs and their phylogenetic relationships with the other PANGOLIN lineages. The mutational landscape of WHO characterized variants was determined and mutational diversity was compared amongst the different severity groups. We then computationally studied the structural impact of the mutations in receptor binding domain of the VOCs. The binding affinity was quantitatively determined by molecular dynamics simulations and free energy calculations.

RESULTS

The mutational frequency, as well as phylogenetic distance, was maximum in the case of omicron followed by the delta variant. The maximum binding affinity was for delta variant followed by the Omicron variant. The increased binding affinity of delta strain followed by omicron as compared to other variants and wild type advocates high transmissibility and quick spread of these two variants and high severity of delta variant.

CONCLUSION

This study delivers a foundation for discovering the improved binding knacks and structural features of SARS-CoV-2 variants to plan novel therapeutics and vaccine candidates against the virus.

Single-Nanoparticle-Based Digital SERS Sensing Platform for the Accurate Quantitative Detection of SARS-CoV-2

To prevent the ongoing spread of the highly infectious severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), accurate and early detection based on a rapid, ultrasensitive, and highly reliable sensing method is crucially important. Here, we present a bumpy core-shell surface-enhanced Raman spectroscopy (SERS) nanoprobe-based sensing platform with single-nanoparticle (SNP)-based digital SERS analysis. The tailorable bumpy core-shell SERS nanoprobe with an internal self-assembled monolayer of 4-nitrobenzenethiol Raman reporters, synthesized using HEPES biological buffer, generates a strong, uniform, and reproducible SERS signal with an SNP-level sensitive and narrowly distributed enhancement factor (2.1 × 108 to 2.2 × 109). We also propose an SNP-based digital SERS analysis method that provides direct visualization of SNP detection at ultralow concentrations and reliable quantification over a wide range of concentrations. The bumpy core-shell SERS nanoprobe-based sensing platform with SNP-based digital SERS analysis achieves the ultrasensitive and quantitative detection of the SARS-CoV-2 spike protein with a limit of detection of 7.1 × 10-16 M over a wide dynamic range from 3.7 × 10-15 to 3.7 × 10-8 M, far outperforming the conventional enzyme-linked immunosorbent assay method for the target protein. Furthermore, it can detect mutated spike proteins from the SARS-CoV-2 variants, representing the key mutations of Alpha, Beta, Gamma, Delta, and Omicron variants. Therefore, this sensing platform can be effectively and efficiently used for the accurate and early detection of SARS-CoV-2 and be adapted for the ultrasensitive and reliable detection of other highly infectious diseases.

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.

An integrated understanding of the evolutionary and structural features of the SARS-CoV-2 spike receptor binding domain (RBD)

Conventional drug development strategies typically use pocket in protein structures as drug-target sites. They overlook the plausible effects of protein evolvability and resistant mutations on protein structure which in turn may impair protein-drug interaction. In this study, we used an integrated evolution and structure guided strategy to develop potential evolutionary-escape resistant therapeutics using receptor binding domain (RBD) of SARS-CoV-2 spike-protein/S-protein as a model. Deploying an ensemble of sequence space exploratory tools including co-evolutionary analysis and deep mutational scans we provide a quantitative insight into the evolutionarily constrained subspace of the RBD sequence-space. Guided by molecular simulation and structure network analysis we highlight regions inside the RBD, which are critical for providing structural integrity and conformational flexibility. Using fuzzy C-means clustering we combined evolutionary and structural features of RBD and identified a critical region. Subsequently, we used computational drug screening using a library of 1615 small molecules and identified one lead molecule, which is expected to target the identified region, critical for evolvability and structural stability of RBD. This integrated evolution-structure guided strategy to develop evolutionary-escape resistant lead molecules have potential general applications beyond SARS-CoV-2.

Investigation of the binding and dynamic features of A.30 variant revealed higher binding of RBD for hACE2 and escapes the neutralizing antibody: A molecular simulation approach

With the emergence of Delta and Omicron variants, many other important variants of SARS-CoV-2, which cause Coronavirus disease-2019, including A.30, are reported to increase the concern created by the global pandemic. The A.30 variant, reported in Tanzania and other countries, harbors spike gene mutations that help this strain to bind more robustly and to escape neutralizing antibodies. The present study uses molecular modelling and simulation-based approaches to investigate the key features of this strain that result in greater infectivity. The protein-protein docking results for the spike protein demonstrated that additional interactions, particularly two salt-bridges formed by the mutated residue Lys484, increase binding affinity, while the loss of key residues at the N terminal domain (NTD) result in a change to binding conformation with monoclonal antibodies, thus escaping their neutralizing effects. Moreover, we deeply studied the atomic features of these binding complexes through molecular simulation, which revealed differential dynamics when compared to wild type. Analysis of the binding free energy using MM/GBSA revealed that the total binding free energy (TBE) for the wild type receptor-binding domain (RBD) complex was -58.25 kcal/mol in contrast to the A.30 RBD complex, which reported -65.59 kcal/mol. The higher TBE for the A.30 RBD complex signifies a more robust interaction between A.30 variant RBD with ACE2 than the wild type, allowing the variant to bind and spread more promptly. The BFE for the wild type NTD complex was calculated to be -65.76 kcal/mol, while the A.30 NTD complex was estimated to be -49.35 kcal/mol. This shows the impact of the reported substitutions and deletions in the NTD of A.30 variant, which consequently reduce the binding of mAb, allowing it to evade the immune response of the host. The reported results will aid the development of cross-protective drugs against SARS-CoV-2 and its variants.

MutCov: A pipeline for evaluating the effect of mutations in spike protein on infectivity and antigenicity of SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), causing an outbreak of coronavirus disease 2019 (COVID-19), is a major threat to public health worldwide. Previous studies have shown that the spike protein of SARS-CoV-2 determines viral infectivity and major antigenicity. However, the spike protein has been undergoing various mutations, which bring a great challenge to the prevention and treatment of COVID-19. Here we present the MutCov, a pipeline for evaluating the effect of mutations in spike protein on infectivity and antigenicity of SARS-CoV-2 by calculating the binding free energy between spike protein and angiotensin-converting enzyme 2 (ACE2) or neutralizing monoclonal antibody (mAb). The predicted infectivity and antigenicity were highly consistent with biologically experimental results, and demonstrated that the MutCov achieved good prediction performance. In conclusion, the MutCov is of high importance for systematically evaluating the effect of novel mutations and improving the prevention and treatment of COVID-19. The source code and installation instruction of MutCov are freely available at http://jianglab.org.cn/MutCov.

The Omicron (B.1.1.529) variant of SARS-CoV-2 binds to the hACE2 receptor more strongly and escapes the antibody response: Insights from structural and simulation data

As SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) continues to inflict chaos globally, a new variant officially known as B.1.1.529 was reported in South Africa and was found to harbor 30 mutations in the spike protein. It is too early to speculate on transmission and hospitalizations. Hence, more analyses are required, particularly to connect the genomic patterns to the phenotypic attributes to reveal the binding differences and antibody response for this variant, which can then be used for therapeutic interventions. Given the urgency of the required analysis and data on the B.1.1.529 variant, we have performed a detailed investigation to provide an understanding of the impact of these novel mutations on the structure, function, and binding of RBD to hACE2 and mAb to the NTD of the spike protein. The differences in the binding pattern between the wild type and B.1.1.529 variant complexes revealed that the key substitutions Asn417, Ser446, Arg493, and Arg498 in the B.1.1.529 RBD caused additional interactions with hACE2 and the loss of key residues in the B.1.1.529 NTD resulted in decreased interactions with three CDR regions (1-3) in the mAb. Further investigation revealed that B.1.1.529 displayed a stable dynamic that follows a global stability trend. In addition, the dissociation constant (KD), hydrogen bonding analysis, and binding free energy calculations further validated the findings. Hydrogen bonding analysis demonstrated that significant hydrogen bonding reprogramming took place, which revealed key differences in the binding. The total binding free energy using MM/GBSA and MM/PBSA further validated the docking results and demonstrated significant variations in the binding. This study is the first to provide a basis for the higher infectivity of the new SARS-CoV-2 variants and provides a strong impetus for the development of novel drugs against them.

Computational prediction of the effect of mutations in the receptor-binding domain on the interaction between SARS-CoV-2 and human ACE2

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causing COVID-19 continues to mutate. Numerous studies have indicated that this viral mutation, particularly in the receptor-binding domain area, may increase the viral affinity for human angiotensin-converting enzyme 2 (hACE2), the receptor for viral entry into host cells, thereby increasing viral virulence and transmission. In this study, we investigated the binding affinity of SARS-CoV-2 variants (Delta plus, Iota, Kappa, Mu, Lambda, and C.1.2) on hACE2 using computational modeling with a protein-protein docking approach. The simulation results indicated that there were differences in the interactions between the RBD and hACE2, including hydrogen bonding, salt bridge interactions, non-bonded interactions, and binding free energy differences among these variants. Molecular dynamics simulations revealed that mutations in the RBD increase the stability of the hACE2-spike protein complex relative to the wild type, following the global stability trend and increasing the binding affinity. The value of binding-free energy calculated using molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) indicated that all mutations in the spike protein increased the contagiousness of SARS-CoV-2 variants. The findings of this study provide a foundation for developing effective interventions against these variants. Computational modeling elucidates that the spike protein of SARS-CoV-2 variants binds considerably stronger than the wild-type to hACE2.

Viral and Host Genetic and Epigenetic Biomarkers Related to SARS-CoV-2 Cell Entry, Infection Rate, and Disease Severity

The rapid spread of COVID-19 outbreak lead to a global pandemic declared in March 2020. The common features of corona virus family helped to resolve structural characteristics and entry mechanism of SARS-CoV-2. However, rapid mutagenesis leads to the emergence of new strains that may have different reproduction rates or infectivity and may impact the course and severity of the disease. Host related factors may also play a role in the susceptibility for infection as well as the severity and outcomes of the COVID-19. We have performed a literature and database search to summarize potential viral and host-related genomic and epigenomic biomarkers, such as genetic variability, miRNA, and DNA methylation in the molecular pathway of SARS-CoV-2 entry into the host cell, that may be related to COVID-19 susceptibility and severity. Bioinformatics tools may help to predict the effect of mutations in the spike protein on the binding to the ACE2 receptor and the infectivity of the strain. SARS-CoV-2 may also target several transcription factors and tumour suppressor genes, thus influencing the expression of different host genes and affecting cell signalling. In addition, the virus may interfere with RNA expression in host cells by exploiting endogenous miRNA and its viral RNA. Our analysis showed that numerous human miRNA may form duplexes with different coding and non-coding regions of viral RNA. Polymorphisms in human genes responsible for viral entry and replication, as well as in molecular damage response and inflammatory pathways may also contribute to disease prognosis and outcome. Gene ontology analysis shows that proteins encoded by such polymorphic genes are highly interconnected in regulation of defense response. Thus, virus and host related genetic and epigenetic biomarkers may help to predict the course of the disease and the response to treatment.

Omicron SARS-CoV-2 Variant Spike Protein Shows an Increased Affinity to the Human ACE2 Receptor: An In Silico Analysis

The rise of SARS-CoV-2 variants, with changes that could be related to an increased virus pathogenicity, have received the interest of the scientific and medical community. In this study, we evaluated the changes that occurred in the viral spike of the SARS-CoV-2 Omicron variant and whether these changes modulate the interactions with the angiotensin-converting enzyme 2 (ACE2) host receptor. The mutations associated with the Omicron variant were retrieved from the GISAID and covariants.org databases, and a structural model was built using the SWISS-Model server. The interaction between the spike and the human ACE2 was evaluated using two different docking software, Zdock and Haddock. We found that the binding free energy was lower for the Omicron variant as compared to the WT spike. In addition, the Omicron spike protein showed an increased number of electrostatic interactions with ACE2 than the WT spike, especially the interactions related to charged residues. This study contributes to a better understanding of the changes in the interaction between the Omicron spike and the human host ACE2 receptor.