Structural-Dynamics and Binding Analysis of RBD from SARS-CoV-2 Variants of Concern (VOCs) and GRP78 Receptor Revealed Basis for Higher Infectivity

Computational electrostatic engineering of nanobodies for enhanced SARS-CoV-2 receptor binding domain recognition

This study presents a novel computational approach for engineering nanobodies (Nbs) for improved interaction with receptor-binding domain (RBD) of the SARS-CoV-2 spike protein. Using Protein Structure Reliability reports, RBD (7VYR_R) was selected and refined for subsequent Nb-RBD interactions. By leveraging electrostatic complementarity (EC) analysis, we engineered and characterized five Electrostatically Complementary Nbs (ECSb1-ECSb5) based on the CeVICA library's SR6c3 Nb. Through targeted modifications in the complementarity-determining regions (CDR) and framework regions (FR), we optimized electrostatic interactions to improve binding affinity and specificity. The engineered Nbs (ECSb3, ECSb4, and ECSb5) demonstrated high binding specificity for AS3, CA1, and CA2 epitopes. Interestingly, ECSb1 and ECSb2 selectively engaged with AS3 and CA1 instead of AS1 and AS2, respectively, due to a preference for residues that conferred superior binding complementarities. Furthermore, ECSbs significantly outperformed SR6c3 Nb in MM/GBSA results, notably, ECSb4 and ECSb3 exhibited superior binding free energies of -182.58 kcal.mol-1 and -119.07 kcal.mol-1, respectively, compared to SR6c3 (-105.50 kcal.mol-1). ECSbs exhibited significantly higher thermostability (100.4-148.3 kcal·mol⁻1) compared to SR6c3 (62.6 kcal·mol⁻1). Similarly, enhanced electrostatic complementarity was also observed for ECSb4-RBD and ECSb3-RBD (0.305 and 0.390, respectively) relative to SR6c3-RBD (0.233). Surface analyses confirmed optimized electrostatic patches and reduced aggregation propensity in the engineered Nb. This integrated EC and structural engineering approach successfully developed engineered Nbs with enhanced binding specificity, increased thermostability, and reduced aggregation, laying the groundwork for novel therapeutic applications targeting the SARS-CoV-2 spike protein.

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 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.

A Suitable Membrane Distance Regulated by the RBD_ACE2 Interaction is Critical for SARS-CoV-2 Spike-Mediated Viral Invasion

The receptor-binding domain (RBD) of spike recognizing the receptor angiotensin-converting enzyme 2 (ACE2) initiates membrane fusion between severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and cell membrane. Although the structure of the RBD_ACE2 complex has been well studied, its functional mechanism in membrane fusion is still not fully understood. Here, using an in vitro cell-vesicle content-mixing assay, it is found that the cleavage at the S2' site by thrombin (Thr) protease strongly accelerates membrane fusion, compared to that of cleavage at the S1/S2 site by PreScission (3C) protease. Moreover, mutations at the RBD_ACE2 interface resulted in a positive correlation between binding affinity and fusion probability. In both the cell-vesicle and cell-cell fusion assays, by crosslinking two membranes via the neutravidin (NTV)_biotin interaction or complementary DNA strands, it is found that spike drives membrane fusion in the absence of ACE2, and a suitable distance between two membranes is critical for spike-mediated membrane fusion. Finally, unsuitable membrane crosslinkers significantly inhibited the fusion probability in the presence of ACE2. Taken together, the results suggest that the RBD_ACE2 complex may act as a crosslinker to bridge the viral and cell membranes at a suitable distance, which is critical, but also substitutable for spike-mediated SARS-CoV-2 entry.

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.

Serum glucose-regulated protein 78 (GRP78) levels in COVID-19-associated mucormycosis: results of a case-control study

BACKGROUND

In experimental models, the expression of glucose-regulated protein 78 (GRP78) in endothelial cells played a role in the pathogenesis of mucormycosis. However, the role of GRP78 in COVID-19-associated mucormycosis (CAM) has not been studied. We hypothesized that serum GRP78 levels are elevated in subjects with CAM.

OBJECTIVE

To compare the serum GRP78 levels in subjects with CAM and COVID-19 controls without mucormycosis.

DESIGN AND SETTING

We performed a hospital-based, case-control study between 1 April 2021 and 31 May 2021.

PARTICIPANTS

We enrolled 24 subjects each of CAM and COVID-19 subjects without mucormycosis. We also measured serum GRP78 levels in ten healthy controls.

EXPOSURE

The primary exposure studied was serum GRP78 concentration, estimated using a commercially available ELISA kit in stored serum samples.

RESULTS

We found the mean ± standard deviation (SD) serum GRP78 levels significantly higher (p = 0.0001) among the CAM (374.3 ± 127.3 pg/mL) than the COVID-19 (246.4 ± 67.0 pg/mL) controls. The proportion of subjects with an abnormal GRP78 level (> mean [184.8 pg/mL] plus two SD [23.2 pg/mL] of GRP78 from healthy participants) was 87.5% and 45.8% in the CAM group and COVID-19 controls, respectively. Serum GRP78 level was independently associated with CAM (odds ratio 1.011; 95% confidence interval [1.002-1.019]) after adjusting for diabetes mellitus and hypoxemia during acute COVID-19.

CONCLUSION

Serum GRP78 levels were significantly higher in CAM than in COVID-19 controls. Further studies are required to the role of GRP78 in the pathogenesis of CAM.

Strain wars 2: Binding constants, enthalpies, entropies, Gibbs energies and rates of binding of SARS-CoV-2 variants

SARS-CoV-2 virus is the cause of COVID-19 pandemic and belongs to RNA viruses, showing great tendency to mutate. Several dozens of mutations have been observed on the SARS-CoV-2 virus, during the last two years. Some of the mutated strains show a greater infectivity and are capable of suppressing the earlier strains, through interference. In this work, kinetic and thermodynamic properties were calculated for strains characterized by various numbers and locations of mutations. It was shown that mutations lead to changes in chemical composition, thermodynamic properties and infectivity. Through competition, the phenomenon of interference of various SARS-CoV-2 strains was explained, which results in suppression of the wild type by mutant strains. Standard Gibbs energy of binding and binding constant for the Omicron (B.1.1.529) strain were found to be Δ_BG⁰ = −45.96 kJ/mol and K_B = 1.13 ∙ 10⁺⁸ M⁻¹, respectively.

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.