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

Mechanistic study of the transmission pattern of the SARS-CoV-2 omicron variant

The omicron variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) characterized by 30 mutations in its spike protein, has rapidly spread worldwide since November 2021, significantly exacerbating the ongoing COVID-19 pandemic. In order to investigate the relationship between these mutations and the variant's high transmissibility, we conducted a systematic analysis of the mutational effect on spike-angiotensin-converting enzyme-2 (ACE2) interactions and explored the structural/energy correlation of key mutations, utilizing a reliable coarse-grained model. Our study extended beyond the receptor-binding domain (RBD) of spike trimer through comprehensive modeling of the full-length spike trimer rather than just the RBD. Our free-energy calculation revealed that the enhanced binding affinity between the spike protein and the ACE2 receptor is correlated with the increased structural stability of the isolated spike protein, thus explaining the omicron variant's heightened transmissibility. The conclusion was supported by our experimental analyses involving the expression and purification of the full-length spike trimer. Furthermore, the energy decomposition analysis established those electrostatic interactions make major contributions to this effect. We categorized the mutations into four groups and established an analytical framework that can be employed in studying future mutations. Additionally, our calculations rationalized the reduced affinity of the omicron variant towards most available therapeutic neutralizing antibodies, when compared with the wild type. By providing concrete experimental data and offering a solid explanation, this study contributes to a better understanding of the relationship between theories and observations and lays the foundation for future investigations.

Molecular modeling and simulation approaches to characterize potential molecular targets for burdock inulin to instigate protection against autoimmune diseases

In the current study, we utilized molecular modeling and simulation approaches to define putative potential molecular targets for Burdock Inulin, including inflammatory proteins such as iNOS, COX-2, TNF-alpha, IL-6, and IL-1β. Molecular docking results revealed potential interactions and good binding affinity for these targets; however, IL-1β, COX-2, and iNOS were identified as the best targets for Inulin. Molecular simulation-based stability assessment demonstrated that inulin could primarily target iNOS and may also supplementarily target COX-2 and IL-1β during DSS-induced colitis to reduce the role of these inflammatory mechanisms. Furthermore, residual flexibility, hydrogen bonding, and structural packing were reported with uniform trajectories, showing no significant perturbation throughout the simulation. The protein motions within the simulation trajectories were clustered using principal component analysis (PCA). The IL-1β-Inulin complex, approximately 70% of the total motion was attributed to the first three eigenvectors, while the remaining motion was contributed by the remaining eigenvectors. In contrast, for the COX2-Inulin complex, 75% of the total motion was attributed to the eigenvectors. Furthermore, in the iNOS-Inulin complex, the first three eigenvectors contributed to 60% of the total motion. Furthermore, the iNOS-Inulin complex contributed 60% to the total motion through the first three eigenvectors. To explore thermodynamically favorable changes upon mutation, motion mode analysis was carried out. The Free Energy Landscape (FEL) results demonstrated that the IL-1β-Inulin achieved a single conformation with the lowest energy, while COX2-Inulin and iNOS-Inulin exhibited two lowest-energy conformations each. IL-1β-Inulin and COX2-Inulin displayed total binding free energies of - 27.76 kcal/mol and - 37.78 kcal/mol, respectively, while iNOS-Inulin demonstrated the best binding free energy results at - 45.89 kcal/mol. This indicates a stronger pharmacological potential of iNOS than the other two complexes. Thus, further experiments are needed to use inulin to target iNOS and reduce DSS-induced colitis and other autoimmune diseases.

SARS-CoV-2 omicron RBD forms a weaker binding affinity to hACE2 compared to Delta RBD in in-silico studies

The COVID-19 pandemic sparked an unprecedented race in biotechnology in a search for effective therapies and a preventive vaccine. The continued appearance of SARS-CoV-2 variants of concern (VoCs) further swept the world. The entry of SARS-CoV-2 into cells is mediated by binding the receptor-binding domain (RBD) of the S protein to the cell-surface receptor, human angiotensin-converting enzyme 2 (hACE2). In this study, using a coarse-grained force field to parameterize the system, we employed steered-molecular dynamics (SMD) simulations to reveal the binding of SARS-CoV-2 Delta/Omicron RBD to hACE2. Our benchmarked results demonstrate a good correlation between computed rupture force and experimental binding free energy for known protein-protein systems. Moreover, our findings show that the Omicron RBD has a weaker binding affinity to hACE2, consistent with the respective experimental results. This indicates that our method can effectively be applied to other emerging SARS-CoV-2 strains.Communicated by Ramaswamy H. Sarma.

Biosensing circulating MicroRNAs in autoinflammatory skin diseases: Focus on Hidradenitis suppurativa

MicroRNAs (miRNAs) play a crucial role in the early diagnosis of autoinflammatory diseases, with Hidradenitis Suppurativa (HS) being a notable example. HS, an autoinflammatory skin disease affecting the pilosebaceous unit, profoundly impacts patients' quality of life. Its hidden nature, with insidious initial symptoms and patient reluctance to seek medical consultation, often leads to a diagnostic delay of up to 7 years. Recognizing the urgency for early diagnostic tools, recent research identified significant differences in circulating miRNA expression, including miR-24-1-5p, miR-146a-5p, miR26a-5p, miR-206, miR338-3p, and miR-338-5p, between HS patients and healthy controls. These miRNAs serve as potential biomarkers for earlier disease detection. Traditional molecular biology techniques, like reverse transcription quantitative-polymerase chain reaction (RT-qPCR), are employed for their detection using specific primers and probes. Alternatively, short peptides offer a versatile and effective means for capturing miRNAs, providing specificity, ease of synthesis, stability, and multiplexing potential. In this context, we present a computational simulation pipeline designed for crafting peptide sequences that can capture circulating miRNAs in the blood of patients with autoinflammatory skin diseases, including HS. This innovative approach aims to expedite early diagnosis and enhance therapeutic follow-up, addressing the critical need for timely intervention in HS and similar conditions.

Structure-guided engineering and molecular simulations to design a potent monoclonal antibody to target aP2 antigen for adaptive immune response instigation against type 2 diabetes

Introduction

Diabetes mellitus (DM) is recognized as one of the oldest chronic diseases and has become a significant public health issue, necessitating innovative therapeutic strategies to enhance patient outcomes. Traditional treatments have provided limited success, highlighting the need for novel approaches in managing this complex disease.

Methods

In our study, we employed graph signature-based methodologies in conjunction with molecular simulation and free energy calculations. The objective was to engineer the CA33 monoclonal antibody for effective targeting of the aP2 antigen, aiming to elicit a potent immune response. This approach involved screening a mutational landscape comprising 57 mutants to identify modifications that yield significant enhancements in binding efficacy and stability.

Results

Analysis of the mutational landscape revealed that only five substitutions resulted in noteworthy improvements. Among these, mutations T94M, A96E, A96Q, and T94W were identified through molecular docking experiments to exhibit higher docking scores compared to the wild-type. Further validation was provided by calculating the dissociation constant (KD), which showed a similar trend in favor of these mutations. Molecular simulation analyses highlighted T94M as the most stable complex, with reduced internal fluctuations upon binding. Principal components analysis (PCA) indicated that both the wild-type and T94M mutant displayed similar patterns of constrained and restricted motion across principal components. The free energy landscape analysis underscored a single metastable state for all complexes, indicating limited structural variability and potential for high therapeutic efficacy against aP2. Total binding free energy (TBE) calculations further supported the superior performance of the T94M mutation, with TBE values demonstrating the enhanced binding affinity of selected mutants over the wild-type.

Discussion

Our findings suggest that the T94M substitution, along with other identified mutations, significantly enhances the therapeutic potential of the CA33 antibody against DM by improving its binding affinity and stability. These results not only contribute to a deeper understanding of antibody-antigen interactions in the context of DM but also provide a valuable framework for the rational design of antibodies aimed at targeting this disease more effectively.

Imaging metabolic mechanisms and the binding behavior of nutrients/transporters of edible Matricaria flowers VOCs

To promote the consumption of flowers and to utilize the nutritional value of proteins, the efficacy of the beneficial components of flowers has been intensively studied. Anthemis nobilis was used as the study object, and all its volatile components (VOCs) were fingerprinted using headspace solid-phase micro-extraction gas-mass spectrometry (HS-SPME/GC-MS). GC-MS fingerprints of five parts of Anthemis nobilis were established using three proteins, bovine lactoferrin (BLF), bovine lactoglobulin (β-Lg), and human serum albumin (HSA), as nutrient transporters. The interactions between the volatile components from different parts of the mother chrysanthemum plant and the nutrient/transport proteins were investigated. The results of fingerprinting showed that the flavor components were dominated by alkenes. In addition, this study revealed that among the three nutrient transporters, the strongest binding to the adsorbed volatile components was HSA, followed by BLF, and β-Lg was second. In addition, a characteristic molecule, camphene, was screened. Integrated molecular simulation using fluorescence spectroscopy was used to validate the results of the interaction of the nutrient/transport proteins systems with characteristic molecule. The properties of the characteristic molecules such as absorption, distribution, metabolism, excretion and toxicity in vivo were analyzed using ADMET to provide a theoretical basis for the preparation of flower-flavored dairy products.

Structural basis for the mechanism of interaction of SARS-CoV-2 B.1.640.2 variant RBD with the host receptors hACE2 and GRP78

The inflicted chaos instigated by the SARS-CoV-2 (Severe acute respiratory syndrome coronavirus 2) globally continues with the emergence of novel variants. The current global outbreak is aggravated by the manifestation of novel variants, which affect the effectiveness of the vaccine, attachment with hACE2 (human Angiotensin-converting enzyme 2) and immune evasion. Recently, a new variant named University Hospital Institute (IHU) (B.1.640.2) was reported in France in November 2021 and is spreading globally affecting public healthcare. The B.1.640.2 SARS-CoV-2 strain revealed 14 mutations and 9 deletions in spike protein. Thus, it is important to understand how these variations in the spike protein impact the communication with the host. A protein coupling approach along with molecular simulation protocols was used to interpret the variation in the binding of the wild type (WT) and B.1.640.2 variant with hACE2 and Glucose-regulating protein 78 (GRP78) receptors. The initial docking scores revealed a stronger binding of the B.1.640.2-RBD with both the hACE2 and GRP78. To further understand the crucial dynamic changes, we looked at the structural and dynamic characteristics and also explored the variations in the bonding networks between the WT and B.1.640.2-RBD (receptor-binding domain) in association with hACE2 and GRP78, respectively. Our findings revealed that the variant complex demonstrated distinct dynamic properties in contrast to the wild type due to the acquired mutations. Finally, to provide conclusive evidence on the higher binding by the B.1.640.2 variant the TBE was computed for each complex. For the WT with hACE2 the TBE was quantified to be-61.38 ± 0.96 kcal/mol and for B.1.640.2 variant the TBE was estimated to be -70.47 ± 1.00 kcal/mol. For the WT-RBD-GRP78 the TBE -was computed to be 32.32 ± 0.56 kcal/mol and for the B.1.640.2-RBD a TBE of -50.39 ± 0.88 kcal/mol was reported. This show that these mutations are the basis for higher binding and infectivity produced by B.1.640.2 variant and can be targeted for drug designing against it.Communicated by Ramaswamy H. Sarma.

Structural insights into the mechanism of resistance to bicalutamide by the clinical mutations in androgen receptor in chemo-treatment resistant prostate cancer

Despite advanced diagnosis and detection technologies, prostate cancer (PCa) is the most prevalent neoplasms in males. Dysregulation of the androgen receptor (AR) is centrally involved in the tumorigenesis of PCa cells. Acquisition of drug resistance due to modifications in AR leads to therapeutic failure and relapse in PCa. An overhaul of comprehensive catalogues of cancer-causing mutations and their juxta positioning on 3D protein can help in guiding the exploration of small drug molecules. Among several well-studied PCa-specific mutations, T877A, T877S and H874Y are the most common substitutions in the ligand-binding domain (LBD) of the AR. In this study, we combined structure as well as dynamics-based in silico approaches to infer the mechanistic effect of amino acid substitutions on the structural stability of LBD. Molecular dynamics simulations allowed us to unveil a possible drug resistance mechanism that acts through structural alteration and changes in the molecular motions of LBD. Our findings suggest that the resistance to bicalutamide is partially due to increased flexibility in the H12 helix, which disturbs the compactness, thereby reducing the affinity for bicalutamide. In conclusion, the current study helps in understanding the structural changes caused by mutations and could assist in the drug development process.Communicated by Ramaswamy H. Sarma.

Novel Omicron Variants Enhance Anchored Recognition of TMEM106B: A New Pathway for SARS-CoV-2 Cellular Invasion

The recent discovery that TMEM106B serves as a receptor mediating ACE2-independent SARS-CoV-2 entry into cells deserves attention, especially in the background of the frequent emergence of mutant strains. Here, the structure-dynamic features of this novel pathway are dissected deeply. Our investigation revealed that the large loop (RBD@471-491) could anchor TMEM106B, which was then firmly locked by another loop (RBD@444-451). The novel and widely disseminated Omicron variants (BA.2.86/EG.5.1) affect the anchoring recognition of proteins, with BA.2.86 being more likely to impact cells with limited or undetectable ACE2 expression. The large loop of the EG.5.1 variant captures TMEM106B poorly due to impaired electrostatic complementarity. Furthermore, we emphasize that antibody design against these two loops could enhance the protection of ACE2 low-expressing cells according to the alanine scanning mutagenesis of multiple antibodies. We hope this study will provide a novel perspective for the prevention and treatment against this new viral invasion pathway.

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.

Mutational analysis of SARS-CoV-2 ORF6-KPNA2 binding interface and identification of potent small molecule inhibitors to recuse the host immune system

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) surfaced on 31 December, 2019, and was identified as the causative agent of the global COVID-19 pandemic, leading to a pneumonia-like disease. One of its accessory proteins, ORF6, has been found to play a critical role in immune evasion by interacting with KPNA2 to antagonize IFN signaling and production pathways, resulting in the inhibition of IRF3 and STAT1 nuclear translocation. Since various mutations have been observed in ORF6, therefore, a comparative binding, biophysical, and structural analysis was used to reveal how these mutations affect the virus's ability to evade the human immune system. Among the identified mutations, the V9F, V24A, W27L, and I33T, were found to have a highly destabilizing effect on the protein structure of ORF6. Additionally, the molecular docking analysis of wildtype and mutant ORF6 and KPNA2 revealed the docking score of - 53.72 kcal/mol for wildtype while, -267.90 kcal/mol, -258.41kcal/mol, -254.51 kcal/mol and -268.79 kcal/mol for V9F, V24A, W27L, and I33T respectively. As compared to the wildtype the V9F showed a stronger binding affinity with KPNA2 which is further verified by the binding free energy (-42.28 kcal/mol) calculation. Furthermore, to halt the binding interface of the ORF6-KPNA2 complex, we used a computational molecular search of potential natural products. A multi-step virtual screening of the African natural database identified the top 5 compounds with best docking scores of -6.40 kcal/mol, -6.10 kcal/mol, -6.09 kcal/mol, -6.06 kcal/mol, and -6.03 kcal/mol for tophit1-5 respectively. Subsequent all-atoms simulations of these top hits revealed consistent dynamics, indicating their stability and their potential to interact effectively with the interface residues. In conclusion, our study represents the first attempt to establish a foundation for understanding the heightened infectivity of new SARS-CoV-2 variants and provides a strong impetus for the development of novel drugs against them.

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.

Structural and molecular investigation of the impact of S30L and D88N substitutions in G9R protein on coupling with E4R from Monkeypox virus (MPXV)

Understanding the pathogenesis mechanism of the Monkeypox virus (MPXV) is essential to guide therapeutic development against the Monkeypox virus. In the current study, we investigated the impact of the only two reported substitutions, S30L, D88N, and S30L-D88N on the G9R of the replication complex in 2022 with E4R using structural modeling, simulation, and free energy calculation methods. From the molecular docking and dissociation constant (KD) results, it was observed that the binding affinity did not increase in the mutants, but the interaction paradigm was altered by these substitutions. Molecular simulation data revealed that these mutations are responsible for destabilization, changes in protein packing, and internal residue fluctuations, which can cause functional variance. Additionally, hydrogen bonding analysis revealed that the estimated number of hydrogen bonds are almost equal among the wild-type G9R and each mutant. The total binding free energy for the wild-type G9R with E4R was -85.00 kcal/mol while for the mutants the TBE was -42.75 kcal/mol, -43.68 kcal/mol, and -48.65 kcal/mol respectively. This shows that there is no direct impact of these two reported mutations on the binding with E4R, or it may affect the whole replication complex or any other mechanism involved in pathogenesis. To explore these variations further, we conducted PCA and FEL analyses. Based on our findings, we speculate that within the context of interaction with E4R, the mutations in the G9R protein might be benign, potentially leading to functional diversity associated with other proteins.Communicated by Ramaswamy H. Sarma.

A novel fully human anti-NT-ANGPTL3 antibody from phage display library exhibits potent ApoB, TG, and LDL-C lowering activities in hyperlipidemia mice

Dyslipidemia is characterized by elevated plasma levels of low-density lipoprotein cholesterol (LDL-C), triglycerides (TG), and TG-rich lipoprotein (TGRLs) in circulation, and is closely associated with the incidence and development of cardiovascular disease. Angiopoietin-like protein 3 (ANGPTL3) deficiency has been identified as a cause of familial combined hypolipidemia in humans, which allows it to be an important therapeutic target for reducing plasma lipids. Here, we report the discovery and characterization of a novel fully human antibody F1519-D95aA against N-terminal ANGPTL3 (NT-ANGPTL3), which potently inhibits NT-ANGPTL3 with a KD as low as 9.21 nM. In hyperlipidemic mice, F1519-D95aA shows higher apolipoprotein B (ApoB) and TG-lowering, and similar LDL-C reducing activity as compared to positive control Evinacumab (56.50% vs 26.01% decrease in serum ApoB levels, 30.84% vs 25.28% decrease in serum TG levels, 23.32% vs 22.52% decrease in serum LDLC levels, relative to vehicle group). Molecular docking and binding energy calculations reveal that the F1519-D95aA-ANGPTL3 complex (10 hydrogen bonds, -65.51 kcal/mol) is more stable than the Evinacumab-ANGPTL3 complex (4 hydrogen bonds, -63.76 kcal/mol). Importantly, F1519-D95aA binds to ANGPTL3 with different residues in ANGPTL3 from Evinacumab, suggesting that F1519-D95aA may be useful for the treatment of patients resistant to Evinacumab. In conclusion, F1519-D95aA is a novel fully human anti-NT-ANGPTL3 antibody with potent plasma ApoB, TG, and LDL-C lowering activities, which can potentially serve as a therapeutic agent for hyperlipidemia and relevant cardiovascular diseases.

Structural plasticity of omicron BA.5 and BA.2.75 for enhanced ACE-dependent entry into cells

The current study investigated the binding variations among the wilt type, Omicron sub-variants BA.2.75 and BA.5, using protein-protein docking, protein structural graphs (P SG), and molecular simulation methods. HADDOCK predicted docking scores and dissociation constant (KD) revealed tighter binding of these sub-variants in contrast to the WT. Further investigation revealed variations in the hub residues, protein sub-networks, and GlobalMetapath in these variants as compared to the WT. A very unusual dynamic for BA.2.75 and BA.5 was observed, and secondary structure transition can also be witnessed in the loops (44-505). The results show that the flexibility of these three loops is increased by the mutations as an allosteric effect and thus enhances the chances of bonding with the nearby residues to connect and form a stable connection. Furthermore, the additional hydrogen bonding contacts steer the robust binding of these variants in contrast to the wild type. The total binding free energy for the wild type was calculated to be -61.38 kcal/mol, while for BA.2.75 and BA.5 variants the T BE was calculated to be -70.42 kcal/mol and 69.78 kcal/mol, respectively. We observed that the binding of BA.2.75 is steered by the electrostatic interactions while the BA.5 additional contacts are due to the vdW (Van der Waal) energy. From these findings, it can be observed the Spike (S) protein is undergoing structural adjustments to bind efficiently to the hACE2 (human angiotensin-converting enzyme 2) receptor and, in turn, increase entry to the host cells. The current study will aid the development of structure-based drugs against these variants.Communicated by Ramaswamy H. Sarma.

Cooperative and structural relationships of the trimeric Spike with infectivity and antibody escape of the strains Delta (B.1.617.2) and Omicron (BA.2, BA.5, and BQ.1)

Herein, we conducted simulations of trimeric Spike from several SARS-CoV-2 variants of concern (Delta and Omicron sub-variants BA.2, BA.5, and BQ.1) and investigated the mechanisms by which specific mutations confer resistance to neutralizing antibodies. We observed that the mutations primarily affect the cooperation between protein domains within and between protomers. The substitutions K417N and L452R expand hydrogen bonding interactions, reducing their interaction with neutralizing antibodies. By interacting with nearby residues, the K444T and N460K mutations in the SpikeBQ.1 variant potentially reduces solvent exposure, thereby promoting resistance to antibodies. We also examined the impact of D614G, P681R, and P681H substitutions on Spike protein structure that may be related to infectivity. The D614G substitution influences communication between a glycine residue and neighboring domains, affecting the transition between up- and -down RBD states. The P681R mutation, found in the Delta variant, enhances correlations between protein subunits, while the P681H mutation in Omicron sub-variants weakens long-range interactions that may be associated with reduced fusogenicity. Using a multiple linear regression model, we established a connection between inter-protomer communication and loss of sensitivity to neutralizing antibodies. Our findings underscore the importance of structural communication between protein domains and provide insights into potential mechanisms of immune evasion by SARS-CoV-2. Overall, this study deepens our understanding of how specific mutations impact SARS-CoV-2 infectivity and shed light on how the virus evades the immune system.

WGX50 mitigates doxorubicin-induced cardiotoxicity through inhibition of mitochondrial ROS and ferroptosis

BACKGROUND

Doxorubicin (DOX)-induced cardiotoxicity (DIC) is a major impediment to its clinical application. It is indispensable to explore alternative treatment molecules or drugs for mitigating DIC. WGX50, an organic extract derived from Zanthoxylum bungeanum Maxim, has anti-inflammatory and antioxidant biological activity, however, its function and mechanism in DIC remain unclear.

METHODS

We established DOX-induced cardiotoxicity models both in vitro and in vivo. Echocardiography and histological analyses were used to determine the severity of cardiac injury in mice. The myocardial damage markers cTnT, CK-MB, ANP, BNP, and ferroptosis associated indicators Fe2+, MDA, and GPX4 were measured using ELISA, RT-qPCR, and western blot assays. The morphology of mitochondria was investigated with a transmission electron microscope. The levels of mitochondrial membrane potential, mitochondrial ROS, and lipid ROS were detected using JC-1, MitoSOX™, and C11-BODIPY 581/591 probes.

RESULTS

Our findings demonstrate that WGX50 protects DOX-induced cardiotoxicity via restraining mitochondrial ROS and ferroptosis. In vivo, WGX50 effectively relieves doxorubicin-induced cardiac dysfunction, cardiac injury, fibrosis, mitochondrial damage, and redox imbalance. In vitro, WGX50 preserves mitochondrial function by reducing the level of mitochondrial membrane potential and increasing mitochondrial ATP production. Furthermore, WGX50 reduces iron accumulation and mitochondrial ROS, increases GPX4 expression, and regulates lipid metabolism to inhibit DOX-induced ferroptosis.

CONCLUSION

Taken together, WGX50 protects DOX-induced cardiotoxicity via mitochondrial ROS and the ferroptosis pathway, which provides novel insights for WGX50 as a promising drug candidate for cardioprotection.

Computational screening of natural products to identify potential inhibitors for human neuropilin-1 (NRP1) receptor to abrogate the binding of SARS-CoV-2 and host cell

Recently, a new variant B.1.1.529 or Omicron variant and its sub-variants (BA2.75, BA.5) of SARS-CoV-2 (Severe acute respiratory virus 2) have been reported with a larger number of mutations in the spike protein and particularly in the RBD (receptor-binding domain). The omicron (B.1.1.529) variant has aggravated the pandemic situation further and needs more analysis for therapeutic development. Keeping in view the urgency of the required data, the current study used molecular modeling and simulation-based methods to target the NRP1 (Neuropilin 1) protein to halt the entry into the host cell. Employing a molecular screening approach to screen the North-East African natural compounds database (NEANCDB) revealed Subereamine B with a docking score of -8.44 kcal/mol, Zinolol with the docking score of -8.05 while Subereamine A with a docking score of -7.88 kcal/mol as the best hits against NRP1. Molecular simulation-based further validation revealed stable dynamics, good structural packing, and dynamic residues flexibility index. Moreover, hydrogen bonding fraction analysis demonstrated the interactions remained sustained during the simulation. Furthermore, the total binding free energy for Subereamine B was -44.24 ±0.91 kcal/mol, for Zinolol -34.32 ±0.40 kcal/mol while for Subereamine A the TBE was calculated to be -41.78 ± 0.36 kcal/mol respectively. This shows that the two arginine-based alkaloids, i.e. Subereamine B and Subereamine A could inhibit the NRP1 more strongly than Zinolol. In conclusion, this study provides a basis for the development of novel drugs against SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

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

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

Molecular dynamics simulations of the spike trimeric ectodomain of the SARS-CoV-2 Omicron variant: structural relationships with infectivity, evasion to immune system and transmissibility

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron is currently the most prevalent SARS-CoV-2 variant worldwide. Herein, we calculated molecular dynamics simulations of the trimeric spikeWT and SpikeBA.1 for 300 ns. Our results show that SpikeBA.1 has more conformational flexibility than SpikeWT. Our principal component analysis (PCA) allowed us to observe a broader spectrum of different conformations for SpikeBA.1, mainly at N-terminal domain (NTD) and receptor-binding domain (RBD). Such increased flexibility could contribute to decreased neutralizing antibody recognition of this variant. Our molecular dynamics data show that the RBDBA.1 easily visits an up-conformational state and the prevalent D614G mutation is pivotal to explain molecular dynamics results for this variant because to lost hydrogen bonding interactions between the residue pairs K854SC/D614SC, Y837MC/D614MC, K835SC/D614SC, T859SC/D614SC. In addition, SpikeBA.1 residues near the furin cleavage site are more flexible than in SpikeWT, probably due to P681H and D614G substitutions. Finally, dynamical cross-correlation matrix (DCCM) analysis reveals that D614G and P681H may allosterically affect the cleavage site S1/S2. Conversely, S2' site may be influenced by residues located between NTD and RBD of a neighboring protomer of the SpikeWT. Such communication may be lost in SpikeBA.1, explaining the changes of the cell tropism in the viral infection. In addition, the movements of the NTDWT and NTDBA.1 may modulate the RBD conformation through allosteric effects. Taken together, our results explain how the structural aspects may explain the observed gains in infectivity, immune system evasion and transmissibility of the Omicron variant.Communicated by Ramaswamy H. Sarma.

A multi-tier computational screening framework to effectively search the mutational space of SARS-CoV-2 receptor binding motif to identify mutants with enhanced ACE2 binding abilities

SARS-CoV-2 gained crucial mutations at the receptor binding domain (RBD) that often changed the course of the pandemic leading to new waves with increased case fatality. Variants are observed with enhanced transmission and immune invasion abilities. Thus, predicting future variants with enhanced transmission ability is a problem of utmost research interest. Here, we have developed a multi-tier exhaustive SARS-CoV-2 mutation screening platform combining MM/GBSA, extensive molecular dynamics simulations, and steered molecular dynamics to identify RBD mutants with enhanced ACE2 binding capability. We have identified four RBM mutations (F490K, S494K, G504F, and the P499L) with significantly higher ACE2 binding abilities than wild-type RBD. Compared to wild-type RBD, they all form stable complexes with more hydrogen bonds and salt-bridge interactions with ACE2. Our simulation data suggest that these mutations allosterically alter the packing of the RBM interface of the RBD-ACE2 complex. As a result, the rupture force required to break the RBD-ACE2 contacts is significantly higher for these mutants.

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.

Targeting SARS-CoV-2 Macrodomain-1 to Restore the Innate Immune Response Using In Silico Screening of Medicinal Compounds and Free Energy Calculation Approaches

Among the different drug targets of SARS-CoV-2, a multi-domain protein known as NSP3 is a critical element of the translational and replication machinery. The macrodomain-I, in particular, has been reported to have an essential role in the viral attack on the innate immune response. In this study, we explore natural medicinal compounds and identify potential inhibitors to target the SARS-CoV-2-NSP3 macrodomain-I. Computational modeling and simulation tools were utilized to investigate the structural-dynamic properties using triplicates of 100 ns MD simulations. In addition, the MM/GBSA method was used to calculate the total binding free energy of each inhibitor bound to macrodomain-I. Two significant hits were identified: 3,5,7,4'-tetrahydroxyflavanone 3'-(4-hydroxybenzoic acid) and 2-hydroxy-3-O-beta-glucopyranosyl-benzoic acid. The structural-dynamic investigation of both compounds with macrodomain-I revealed stable dynamics and compact behavior. In addition, the total binding free energy for each complex demonstrated a robust binding affinity, of ΔG -61.98 ± 0.9 kcal/mol for Compound A, while for Compound B, the ΔG was -45.125 ± 2.8 kcal/mol, indicating the inhibitory potential of these compounds. In silico bioactivity and dissociation constant (KD) determination for both complexes further validated the inhibitory potency of each compound. In conclusion, the aforementioned natural products have the potential to inhibit NSP3, to directly rescue the host immune response. The current study provides the basis for novel drug development against SARS-CoV-2 and its variants.

Structural communication fingerprinting and dynamic investigation of RBD-hACE2 complex from BA.1 × AY.4 recombinant variant (Deltacron) of SARS-CoV-2 to decipher the structural basis for enhanced transmission

The BA.1 × AY.4 recombinant variant (Deltacron) continues to inflict chaos globally due to its rapid transmission and infectivity. To decipher the mechanism of pathogenesis by the BA.1 × AY.4 recombinant variant (Deltacron), a protein coupling, protein structural graphs (PSG), residue communication and all atoms simulation protocols were used. We observed that the bonding network is altered by this variant; engaging new residues that helps to robustly bind. HADDOCK docking score for the wild type has been previously reported to be -111.8 ± 1.5 kcal/mol while the docking score for the Deltacron variant was calculated to be -128.3 ± 2.5 kcal/mol. The protein structural graphs revealed variations in the hub residues, number of nodes, inter and intra residues communities, and path communication perturbation caused by the acquired mutations in the Deltacron-RBD thus alter the binding approach and infectivity. Moreover, the dynamic behaviour reported a highly flexible structure with enhanced residues flexibility particularly by the loops required for interaction with ACE2. It was observed that these mutations have altered the secondary structure of the RBD mostly transited to the loops thus acquired higher flexible dynamics than the native structure during the simulation. The total binding free energy for each of these complexes, that is, WT-RBD and Deltacron-RBD were reported to be -61.38 kcal/mol and -70.47 kcal/mol. Protein's motion revealed a high trace value in the Deltacron variant that clearly depict more structural flexibility. The broad range of phase space covered by the Deltacron variant along PC1 and PC2 suggests that these mutations are important in contributing conformational heterogeneity or flexibility that consequently help the variant to bind more efficiently than the wild type. The current study provides a basis for structure-based drug designing against SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Exploring the natural products chemical space through a molecular search to discover potential inhibitors that target the hypoxia-inducible factor (HIF) prolyl hydroxylase domain (PHD)

Introduction: Hypoxia-inducible factor (HIF) prolyl hydroxylase domain (PHD) enzymes are major therapeutic targets of anemia and ischemic/hypoxia diseases. To overcome safety issues, liver failure, and problems associated with on-/off-targets, natural products due to their novel and unique structures offer promising alternatives as drug targets. Methods: In the current study, the Marine Natural Products, North African, South African, East African, and North-East African chemical space was explored for HIF-PHD inhibitors discovery through molecular search, and the final hits were validated using molecular simulation and free energy calculation approaches. Results: Our results revealed that CMNPD13808 with a docking score of -8.690 kcal/mol, CID15081178 with a docking score of -8.027 kcal/mol, CID71496944 with a docking score of -8.48 kcal/mol and CID11821407 with a docking score of -7.78 kcal/mol possess stronger activity than the control N-[(4-hydroxy-8-iodoisoquinolin-3-yl)carbonyl]glycine, 4HG (-6.87 kcal/mol). Interaction analysis revealed that the target compounds interact with Gln239, Tyr310, Tyr329, Arg383 and Trp389 residues, and chelate the active site iron in a bidentate manner in PHD2. Molecular simulation revealed that these target hits robustly block the PHD2 active site by demonstrating stable dynamics. Furthermore, the half-life of the Arg383 hydrogen bond with the target ligands, which is an important factor for PHD2 inhibition, remained almost constant in all the complexes during the simulation. Finally, the total binding free energy of each complex was calculated as CMNPD13808-PHD2 -72.91 kcal/mol, CID15081178-PHD2 -65.55 kcal/mol, CID71496944-PHD2 -68.47 kcal/mol, and CID11821407-PHD2 -62.06 kcal/mol, respectively. Conclusion: The results show the compounds possess good activity in contrast to the control drug (4HG) and need further in vitro and in vivo validation for possible usage as potential drugs against HIF-PHD2-associated diseases.

Beta, Delta, and Omicron, Deadliest Among SARS-CoV-2 Variants: A Computational Repurposing Approach

SARS-CoV-2 has been highly susceptible to mutations since its emergence in Wuhan, China, and its subsequent propagation due to containing an RNA as its genome. The emergence of variants with improved transmissibility still poses a grave threat to global health. The spike protein mutation is mainly responsible for higher transmissibility and risk severity. This study retrieved SARS-CoV-2 variants structural and nonstructural proteins (NSPs) sequences from several geographic locations, including Africa, Asia, Europe, Oceania, and North and South America. First, multiple sequence alignments with BioEdit and protein homology modeling were performed using the SWISS Model. Then the structure visualization and structural analysis were performed by superimposing against the Wuhan sequence by Pymol to retrieve the RMSD values. Sequence alignment revealed familiar, uncommon regional among variants and, interestingly, a few unique mutations in Beta, Delta, and Omicron. Structural analysis of such unique mutations revealed that they caused structural deviations in Beta, Delta, and Omicron spike proteins. In addition, these variants were more severe in terms of hospitalization, sickness, and higher mortality, which have a substantial relationship with the structural deviations because of those unique mutations. Such evidence provides insight into the SARS-CoV-2 spike protein vulnerability toward mutation and their structural and functional deviations, particularly in Beta, Delta, and Omicron, which might be the cause of their broader coverage. This knowledge can help us with regional vaccine strain selection, virus pathogenicity testing, diagnosis, and treatment with more specific vaccines.

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.

Effects of vaccines on clinical characteristics of convalescent adult patients infected with SARS-CoV-2 Omicron variant: A retrospective study

Introduction

The protective effect of SARS-CoV-2 vaccines has become a global focus due to Omicron variant pandemic. The effects of various SARS-CoV-2 vaccines are diverse. However, studies on the effect of domestic vaccines on clinical characteristics in convalescent adult patients infected with the Omicron variant are lacking.

Methods

In this retrospective, single-center cohort study, the effect of three domestic vaccines on clinical characteristics of convalescent adult patients infected with the Omicron variant was investigated in the initial largest outbreak of the Omicron variant infection between January and February 2022 in Tianjin, China. The primary endpoint was COVID-19 severity and the secondary endpoints were re-positive results on nucleic acid tests, liver and kidney function, and inflammation levels during recovery.

Results

A total of 320 adult patients infected with the Omicron variant were enrolled, including 296 post-vaccination and 24 unvaccinated patients. The median age of the unvaccinated patients was higher than that of vaccinated patients, but no significant difference was detected in the sex composition ratio between the different groups. Binary logistic regression results suggested that Sinopharm and Sinovac vaccine was an independent protective factor for relieving the severity of the Omicron variant infection. Regrettably, the vaccines did not showed any protective effect on the liver and kidney function of convalescent adult patients. Three domestic vaccines significantly relieved inflammation and increased the SARS-CoV-2-specific antibody levels. Furthermore, Sinovac and CanSino vaccines had a better immune stimulation effect on increasing T lymphocytes levels in convalescent adult patients. In addition, three domestic vaccines have protective effects on preventing re-detectable positive (RP) result in convalescent adult patients.

Conclusion

Although the three domestic vaccines cannot prevent the infection of the Omicron variant, it has a significant protective effect in adult patients. This study supports the policy of accelerating to vaccination worldwide combat the evolving and mutating SARS-CoV-2.

Discussion

Omicron spreads faster and might escape antibodies more readily than previous variants, increasing the cases of reinfection and breakthrough infections in vaccinated people. Although vaccinated people are likely to have a much lower risk of severe disease from Omicron infection, many issues still need to be considered. Concerns about lower vaccine efficacy because of new variants might have changed our understanding of the COVID-19 endgame, disabusing the world of the notion that global vaccination is by itself adequate for controlling SARS-CoV-2 infection. The current data showed that vaccination with three domestic SARS-CoV-2 vaccines alleviates the disease severity of adult patients with COVID-19, reduces the inflammation level and the RP rate of convalescent adult patients, and enhances body’s defense against the virus in convalescent adult patients. Moreover, our study has highlighted that a combination prevention approach of vaccination and public health measures would be an effective strategy.

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

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

Deciphering the mechanism of resistance by novel double mutations in pncA in Mycobacterium tuberculosis using protein structural graphs (PSG) and structural bioinformatic approaches

The evolution of MDR and XDR-TB is a growing concern and public health safety threat around the world. Gene mutations are the prime cause of drug resistance in tuberculosis, however the reports of double mutations further aggravated the situation. Despite the large-scale genomic sequencing and identification of novel mutations, structure investigation of the protein is still required to structurally and functionally characterize these novel mutations to design novel drugs for improved clinical outcome. Hence, we used structural bioinformatics approaches i.e. molecular modeling, residues communication and molecular simulation to understand the impact of novel double S59Y-L85P, D86G-V180F and S104G-V130 M mutation on the structure, function of pncA encoded Pyrazinamidase (PZase) and resistance of Pyrazinamide (PZA). Our results revealed that these mutations alter the binding paradigm and destabilize the protein to release the drug. Protein commination network (PCN) revealed variations in the hub residues and sub-networks which consequently alter the internal communication and signaling. The region 1-75 demonstrated higher flexibility in the mutant structures and minimal by the wild type which destabilize of the internally arranged beta-sheets which consequently reduce the binding of PZA and potentially Fe ion in the mutants. Hydrogen bonding analysis further validated the findings. The total binding free energy (ΔG) for each complex i.e. wild type -7.46 kcal/mol, S59Y-L85P -5.21 kcal/mol, S104G-V130 M -5.33 kcal/mol while for the D86G-V180F mutant the TBE was calculated to be -6.26 kcal/mol. This further confirms that these mutations reduce the binding energy of PZA for PZase and causes resistance in the effective therapy for TB. The trajectories motion was also observed to be affected by these mutations. In conclusion, these mutations use destabilizing approach to reduce the binding of PZA and causes resistance. These features can be used to design novel structure-based drugs against Tuberculosis.

Baculoviral COVID-19 Delta DNA vaccine cross-protects against SARS-CoV2 variants in K18-ACE2 transgenic mice

After severe acute respiratory syndrome coronavirus-2 (SARS-CoV2) made the world tremble with a global pandemic, SARS-CoV2 vaccines were developed. However, due to the coronavirus's intrinsic nature, new variants emerged, such as Delta and Omicron, refractory to the vaccines derived using the original Wuhan strain. We developed an HERV-enveloped recombinant baculoviral DNA vaccine against SARS-CoV2 (AcHERV-COVID19S). A non-replicating recombinant baculovirus that delivers the SARS-CoV2 spike gene showed a protective effect against the homologous challenge in a K18-hACE2 Tg mice model; however, it offered only a 50 % survival rate against the SARS-CoV2 Delta variant. Therefore, we further developed the AcHERV-COVID19 Delta vaccine (AcHERV-COVID19D). The AcHERV-COVID19D induced higher neutralizing antibodies against the Delta variant than the prototype or Omicron variant. On the other hand, cellular immunity was similarly high for all three SARS-CoV2 viruses. Cross-protection experiments revealed that mice vaccinated with the AcHERV-COVID19D showed 100 % survival upon challenge with Delta and Omicron variants and 71.4 % survival against prototype SARS-CoV2. These results support the potential of the viral vector vaccine, AcHERV-COVID19D, in preventing the spread of coronavirus variants such as Omicron and SARS-CoV2 variants.

Discovery of Isojacareubin as a covalent inhibitor of SARS-CoV-2 main protease using structural and experimental approaches

The ongoing pandemic with the emergence of immune evasion potential and, particularly, the current omicron subvariants intensified the situation further. Although vaccines are available, the immune evasion capabilities of the recent variants demand further efficient therapeutic choices to control the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic. Hence, considering the necessity of the small molecule inhibitor, we target the main protease (3CLpro), which is an appealing target for the development of antiviral drugs against SARS-CoV-2. High-throughput molecular in silico screening of South African natural compounds database reported Isojacareubin and Glabranin as the potential inhibitors for the main protease. The calculated docking scores were reported to be -8.47 and -8.03 kcal/mol, respectively. Moreover, the structural dynamic assessment reported that Isojacareubin in complex with 3CLpro exhibit a more stable dynamic behavior than Glabranin. Inhibition assay indicated that Isojacareubin could inhibit SARS-CoV-2 3CLpro in a time- and dose-dependent manner, with half maximal inhibitory concentration values of 16.00 ± 1.35 μM (60 min incubation). Next, the covalent binding sites of Isojacareubin on SARS-CoV-2 3CLpro was identified by biomass spectrometry, which reported that Isojacareubin can covalently bind to thiols or Cysteine through Michael addition. To evaluate the inactivation potency of Isojacareubin, the inactivation kinetics was further investigated. The inactivation kinetic curves were plotted according to various concentrations with gradient-ascending incubation times. The KI value of Isojacareubin was determined as 30.71 μM, whereas the Kinact value was calculated as 0.054 min-1 . These results suggest that Isojacareubin is a covalent inhibitor of SARS-CoV-2 3CLpro .

Network pharmacology- and molecular simulation-based exploration of therapeutic targets and mechanisms of heparin for the treatment of sepsis/COVID-19

Critically infected patients with COVID-19 (coronavirus disease 2019) are prone to develop sepsis-related coagulopathy as a result of a robust immune response. The mechanism underlying the relationship between sepsis and COVID-19 is largely unknown. LMWH (low molecular weight heparin) exhibits both anti-inflammatory and anti-coagulating properties that result in a better prognosis of severely ill patients with COVID-19 co-associated with sepsis-induced coagulopathy or with a higher D-dimer value. Heparin-associated molecular targets and their mechanism of action in sepsis/COVID-19 are not well understood. In this work, we characterize the pharmacological targets, biological functions and therapeutic actions of heparin in sepsis/COVID-19 from the perspective of network pharmacology. A total of 38 potential targets for heparin action against sepsis/COVID-19 and 8 core pharmacological targets were identified, including IL6, KNG1, CXCL8, ALB, VEGFA, F2, IL10 and TNF. Moreover, enrichment analysis showed that heparin could help in treating sepsis/COVID-19 through immunomodulation, inhibition of the inflammatory response, regulation of angiogenesis and antiviral activity. The pharmacological effects of heparin against these targets were further confirmed by molecular docking and simulation analysis, suggesting that heparin exerts effective binding capacity by targeting the essential residues in sepsis/COVID-19. Prospective clinical practice evaluations may consider the use of these key prognostic indicators for the treatment of sepsis/COVID-19.Communicated by Ramaswamy H. Sarma.

Recent Developments in Vaccines for Viral Diseases

The world is continuously facing the threat of emerging infectious diseases [...].

Global emerging Omicron variant of SARS-CoV-2: Impacts, challenges and strategies

Newly emerging variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are continuously posing high global public health concerns and panic resulting in waves of coronavirus disease 2019 (COVID-19) pandemic. Depending on the extent of genomic variations, mutations and adaptation, few of the variants gain the ability to spread quickly across many countries, acquire higher virulency and ability to cause severe disease, morbidity and mortality. These variants have been implicated in lessening the efficacy of the current COVID-19 vaccines and immunotherapies resulting in break-through viral infections in vaccinated individuals and recovered patients. Altogether, these could hinder the protective herd immunity to be achieved through the ongoing progressive COVID-19 vaccination. Currently, the only variant of interest of SARS-CoV-2 is Omicron that was first identified in South Africa. In this review, we present the overview on the emerging SARS-CoV-2 variants with a special focus on the Omicron variant, its lineages and hybrid variants. We discuss the hypotheses of the origin, genetic change and underlying molecular mechanism behind higher transmissibility and immune escape of Omicron variant. Major concerns related to Omicron including the efficacy of the current available immunotherapeutics and vaccines, transmissibility, disease severity, and mortality are discussed. In the last part, challenges and strategies to counter Omicron variant, its lineages and hybrid variants amid the ongoing COVID-19 pandemic are presented.

Mutational Insight into Allosteric Regulation of Kir Channel Activity

Potassium (K+) channels are regulated in part by allosteric communication between the helical bundle crossing, or inner gate, and the selectivity filter, or outer gate. This network is triggered by gating stimuli. In concert, there is an allosteric network which is a conjugated set of interactions which correlate long-range structural rearrangements necessary for channel function. Inward-rectifier K+ (Kir) channels favor inward K+ conductance, are ligand-gated, and help establish resting membrane potentials. KirBac1.1 is a bacterial Kir (KirBac) channel homologous to human Kir (hKir) channels. Additionally, KirBac1.1 is gated by the anionic phospholipid ligand phosphatidylglycerol (PG). In this study, we use site-directed mutagenesis to investigate residues involved in the KirBac1.1 gating mechanism and allosteric network we previously proposed using detailed solid-state NMR (SSNMR) measurements. Using fluorescence-based K+ and sodium (Na+) flux assays, we identified channel mutants with impaired function that do not alter selectivity of the channel. In tandem, we performed coarse grain molecular dynamics simulations, observing changes in PG-KirBac1.1 interactions correlated with mutant channel activity and contacts between the two transmembrane helices and pore helix tied to this behavior. Lipid affinity is closely tied to the proximity of two tryptophan residues on neighboring subunits which lure anionic lipids to a cationic pocket formed by a cluster of arginine residues. Thus, these simulations establish a structural and functional basis for the role of each mutated site in the proposed allosteric network. The experimental and simulated data provide insight into key functional residues involved in gating and lipid allostery of K+ channels. Our findings also have direct implications on the physiology of hKir channels due to conservation of many of the residues identified in this work from KirBac1.1.

Insights from computational studies on the potential of natural compounds as inhibitors against SARS-CoV-2 spike omicron variant

Coronavirus disease 2019 (COVID-19) is a major global health emergency, with more than six million deaths worldwide. It is becoming increasingly challenging to treat COVID-19 due to the emergence of novel variants. The omicron variant is capable to evade defences and spread quickly. Among many validated COVID-19 targets, the spike (S) protein plays an important role in receptor recognition (via the S1 subunit) and membrane fusion (via the S2 subunit). The S protein is one of the vital targets for the development of drugs to combat this illness. In this research, we applied various computational methods such as molecular docking, molecular dynamics, MM-GBSA calculations, and ADMET prediction to identify potential natural products from Saudi medicinal plants against the spike omicron variant. As a result, three compounds (LTS0002490, LTS0117007, and LTS0217912) were identified with better binding affinity to the spike omicron variant compared to the reference compound (VE607). In addition, these compounds showed stable interactions with the target during molecular dynamics simulations for 140 ns. Last, these compounds have optimal ADMET properties. We suggest that these compounds may be considered promising hits to treat COVID-19 if experimentally validated.

The role of SARS-CoV-2 accessory proteins in immune evasion

Many questions on the SARS-CoV-2 pathogenesis remain to answer. The SARS-CoV-2 genome encodes some accessory proteins that are essential for infection. Notably, accessory proteins of SARS-CoV-2 play significant roles in affecting immune escape and viral pathogenesis. Therefore SARS-CoV-2 accessory proteins could be considered putative drug targets. IFN-I and IFN-III responses are the primary mechanisms of innate antiviral immunity in infection clearance. Previous research has shown that SARS-CoV-2 suppresses IFN-β by infecting host cells via ORF3a, ORF3b, ORF6, ORF7a, ORF7b, ORF8, and ORF9b. Furthermore, ORF3a, ORF7a, and ORF7b have a role in blocking IFNα signaling, and ORF8 represses IFNβ signaling. The ORF3a, ORF7a, and ORF7b disrupt the STAT1/2 phosphorylation. ORF3a, ORF6, ORF7a, and ORF7b could prevent the ISRE promoter activity. The main SARS-CoV-2 accessory proteins involved in immune evasion are discussed here for comprehensive learning on viral entry, replication, and transmission in vaccines and antiviral development.

Insight into SARS-CoV-2 Omicron variants in Saudi Arabian genomic isolates

OBJECTIVES

To investigate the emergent mutations involved in the evolutionary stages of the virus for better management of pandemic.

METHODS

This cross-sectional genomic investigation was performed on February 28, 2022, at the Biology Department, Faculty of Science, Tabuk University. Numerous mutations were searched in genomic isolates of Omicron variant prevalent in the Kingdom of Saudi Arabian. Whole-genome sequences were retrieved from genomic databases and were subjected to the Global Initiative on Sharing Avian Influenza Data (GISAID) CoVsurver for the Omicron variant detection and mutations.

RESULTS

Approximately 8.755 million SARS-CoV-2 genomes were reported to GISAID on February 28, 2022, of which 1270 have been reported from the Kingdom of Saudi Arabia. Among the 1270 genomes, 30 were Omicron variants. Among the Saudi Arabian genomes, 30 were detected as Omicron variants. Twenty-four unique mutations have been detected in membrane, envelope, spike and non-structural proteins (NSP) 12, NSP3, and NSP2. Ten of these unique mutations have been detected in spike protein.

CONCLUSION

The current study provides useful information for further experimental investigation of mutation's effects on virus transmission, severity, and vaccine efficacy.

Evolutionary trajectory of receptor binding specificity and promiscuity of the spike protein of SARS-CoV-2

SARS-CoV-2 infects cells by attachment to its receptor-the angiotensin converting enzyme 2 (ACE2). Regardless of the wealth of structural data, little is known about the physicochemical mechanism of interactions of the viral spike (S) protein with ACE2 and how this mechanism has evolved during the pandemic. Here, we applied experimental and computational approaches to characterize the molecular interaction of S proteins from SARS-CoV-2 variants of concern (VOC). Data on kinetics, activation-, and equilibrium thermodynamics of binding of the receptor binding domain (RBD) from VOC with ACE2 as well as data from computational protein electrostatics revealed a profound remodeling of the physicochemical characteristics of the interaction during the evolution. Thus, as compared to RBDs from Wuhan strain and other VOC, Omicron RBD presented as a unique protein in terms of conformational dynamics and types of non-covalent forces driving the complex formation with ACE2. Viral evolution resulted in a restriction of the RBD structural dynamics, and a shift to a major role of polar forces for ACE2 binding. Further, we investigated how the reshaping of the physicochemical characteristics of interaction affects the binding specificity of S proteins. Data from various binding assays revealed that SARS-CoV-2 Wuhan and Omicron RBDs manifest capacity for promiscuous recognition of unrelated human proteins, but they harbor distinct reactivity patterns. These findings might contribute for mechanistic understanding of the viral tropism and capacity to evade immune responses during evolution.

Ultra-Large-Scale Screening of Natural Compounds and Free Energy Calculations Revealed Potential Inhibitors for the Receptor-Binding Domain (RBD) of SARS-CoV-2

The emergence of immune-evading variants of SARS-CoV-2 further aggravated the ongoing pandemic. Despite the deployments of various vaccines, the acquired mutations are capable of escaping both natural and vaccine-induced immune responses. Therefore, further investigation is needed to design a decisive pharmacological treatment that could efficiently block the entry of this virus into cells. Hence, the current study used structure-based methods to target the RBD of the recombinant variant (Deltacron) of SARS-CoV-2, which was used as a model variant. From the virtual drug screenings of various databases, a total of four hits were identified as potential lead molecules. Key residues were blocked by these molecules with favorable structural dynamic features. The binding free energies further validated the potentials of these molecules. The TBE for MNP was calculated to be -32.86 ± 0.10 kcal/mol, for SANC00222 the TBE was -23.41 ± 0.15 kcal/mol, for Liriodenine the TBE was -34.29 ± 0.07 kcal/mol, while for Carviolin the TBE was calculated to be -27.67 ± 0.12 kcal/mol. Moreover, each complex demonstrated distinct internal motion and a free energy profile, indicating a different strategy for the interaction with and inhibition of the RBD. In conclusion, the current study demands further in vivo and in vitro validation for the possible usage of these compounds as potential drugs against SARS-CoV-2 and its variants.

Prediction of viral-host interactions of COVID-19 by computational methods

Experimental approaches are currently used to determine viral-host interactions, but these approaches are both time-consuming and costly. For these reasons, computational-based approaches are recommended. In this study, using computational-based approaches, viral-host interactions of SARS-CoV-2 virus and human proteins were predicted. The study consists of four different stages; in the first stage viral and host protein sequences were obtained. In the second stage, protein sequences were converted into numerical expressions by various protein mapping methods. These methods are entropy-based, AVL-tree, FIBHASH, binary encoding, CPNR, PAM250, BLOSUM62, Atchley factors, Meiler parameters, EIIP, AESNN1, Miyazawa energies, Micheletti potentials, Z-scale, and hydrophobicity. In the third stage, a deep learning model was designed and BiLSTM was used for this. In the last stage, the protein sequences were classified, and the viral-host interactions were predicted. The performances of protein mapping methods were determined by accuracy, F1-score, specificity, sensitivity, and AUC scores. According to the classification results, the best classification process was obtained by the entropy-based method. With this method, 94.74% accuracy, and 0.95 AUC score were calculated. Then, the most successful classification process was performed with the Z-scale and 91.23% accuracy, and 0.96 AUC score were obtained. Although other protein mapping methods are not as efficient as Z-scale and entropy-based methods, they have achieved successful classification. AVL-tree, FIBHASH, binary encoding, CPNR, PAM250, BLOSUM62, Atchley factors, Meiler parameters and AESNN1 methods showed over 80% accuracy, F1-score, and AUC score. Accuracy scores of EIIP, Miyazawa energies, Micheletti potentials and hydrophobicity methods remained below 80%. When the results were examined in general, it was observed that the computational approaches were successful in predicting viral-host interactions between SARS-CoV-2 virus and human proteins.

SARS-CoV-2 VOCs, Mutational diversity and clinical outcome: Are they modulating drug efficacy by altered binding strength?

The global COVID-19 pandemic continues due to emerging Severe Acute Respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOC). Here, we performed comprehensive analysis of in-house sequenced SARS-CoV-2 genome mutations dynamics in the patients infected with the VOCs - Delta and Omicron, within Recovered and Mortality patients. Statistical analysis highlighted significant mutations - T4685A, N4992N, and G5063S in RdRp; T19R in NTD spike; K444N and N532H in RBD spike, associated with Delta mortality. Mutations, T19I in NTD spike, Q493R and N440K in the RBD spike were significantly associated with Omicron mortality. We performed molecular docking for possible effect of significant mutations on the binding of Remdesivir. We found that Remdesivir showed less binding efficacy with the mutant Spike protein of both Delta and Omicron mortality compared to recovered patients. This indicates that mortality associated mutations could have a modulatory effect on drug binding which could be associated with disease outcome.

Virulence and Biofilm Inhibition of 3-Methoxycinnamic Acid against Agrobacterium tumefaciens

AIMS

In the current study the anti-virulence and anti-biofilm activities of the cinnamic acid derivative, 3-methoxycinnamic acid, was investigated against Agrobacterium tumefaciens.

METHODS AND RESULTS

Based on the disc diffusion test and β-galactosidase activity assay, 3-methoxycinnamic acid was shown to interfere with the quorum sensing (QS) system of A. tumefaciens. Crystal violet staining assay, phenol-sulfuric acid method, Bradford protein assay and confocal laser scanning microscopy (CLSM) revealed that the biofilm formation of A. tumefaciens was inhibited after the treatment of 3-methoxycinnamic acid. Employing high performance liquid chromatography (HPLC) analysis of culture supernatant revealed that the production of 3-oxo-octanoylhomoserine lactone (3-oxo-C8-HSL) decreased concentration-dependently after treatment with 3-methoxycinnamic acid. Swimming and chemotaxis assays also indicated that 3-methoxycinnamic acid had a good effect on reducing the motility and chemotaxis of A. tumefaciens. In addition, the RT-qPCR, molecular docking and simulations further demonstrated that 3-methoxycinnamic acid could competitively inhibit the binding of 3-oxo-C8-HSL to TraR and down-regulate virulence-related genes.

CONCLUSIONS

3-Methoxycinnamic acid is proved to have good anti-virulence and anti-biofilm activities against A. tumefaciens.

SIGNIFICANCE AND IMPACT OF THE STUDY

This is the first study that investigates the anti-virulence and anti-biofilm activities of 3-methoxycinnamic acid against A. tumefaciens. With its potential QS-related virulence and biofilm inhibitory activities, 3-methoxycinnamic acid is expected to be developed as a potent pesticide or adjuvant for the prevention and treatment of crown gall caused by A. tumefaciens.

Density functional theory computation of the binding free energies between various mutations of SARS-CoV-2 RBD and human ACE2: molecular level roots of the contagiousness

The receptor-binding domain (RBD) of SARS-CoV-2 attaches to the human ACE2 to initiate binding of SARS-CoV-2 to human cell and leads to the infection process afterwards. In this study, various mutations of SARS-CoV-2 spike RBD and human ACE2 complexes are investigated via density functional theory (DFT) computations to obtain binding free energies. The DFT computations are performed without fragmenting the interfaces to involve longer-range quantum mechanical interactions for improving accuracy. The vibrational free energies, van der Waals dispersion forces and basis set superposition error corrections are also included in the calculations. The results show that the absolute value of the binding energy of B.1.1.7 mutated spike RBD-ACE2 complex is more than five times higher than that of the original strain. The results of this study are expected to be useful for a deeper understanding of the relation of the binding free energies and the level of contagiousness.

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.

Amomum tsao-ko essential oil, a novel anti-COVID-19 Omicron spike protein natural products: A computational study

Since the outbreak of COVID-19, this virus has been constantly mutating. The latest mutant Omicron has been identified as VOC by WHO. The main reason for its concern is the mutation of 46 amino acids in spike protein, which has brought the global epidemic prevention into another difficulty. Herbal aromatic plant Amomum tsao-ko was excavated from formula 1 and 2 for the treatment of COVID-19 in China, and its active components were extracted and identified. Molecular dynamics simulation and Fpocket were applied to find the key sites on RBDOmicron, and molecular docking was also used to reveal the interaction between A. tsao-ko essential oil (AEO) and RBDOmicron. The AEO components were analyzed and identified by GC/Q-TOF MS. There were 20 kinds of AEO with a relative area percentage of more than 1%, and they were related to the three active centres of RBDOmicron. In this study, virtual screening was used to mine the essential oil components of medicinal plants, and it was found that the components could interact with the spike protein RBD in aerosol to block the interaction of RBD and hACE2, thus cutting off the transmission route and protecting the host. This study has certain guiding significance in the modernization of Traditional Chinese medicine and the prevention of COVID-19.