Molecular basis of SARS-CoV-2 Omicron variant receptor engagement and antibody evasion and neutralization

The SARS-CoV-2 Omicron variant exhibits striking immune evasion and is spreading globally at an unprecedented speed. Understanding the underlying structural basis of the high transmissibility and greatly enhanced immune evasion of Omicron is of high importance. Here through cryo-EM analysis, we present both the closed and open states of the Omicron spike, which appear more compact than the counterparts of the G614 strain, potentially related to the Omicron substitution induced enhanced protomer-protomer and S1-S2 interactions. The closed state showing dominant population may indicate a conformational masking mechanism of immune evasion for Omicron spike. Moreover, we capture two states for the Omicron S/ACE2 complex with S binding one or two ACE2s, revealing that the substitutions on the Omicron RBM result in new salt bridges/H-bonds and more favorable electrostatic surface properties, together strengthened interaction with ACE2, in line with the higher ACE2 affinity of the Omicron relative to the G614 strain. Furthermore, we determine cryo-EM structures of the Omicron S/S3H3 Fab, an antibody able to cross-neutralize major variants of concern including Omicron, elucidating the structural basis for S3H3-mediated broad-spectrum neutralization. Our findings shed new lights on the high transmissibility and immune evasion of the Omicron variant and may also inform design of broadly effective vaccines against emerging variants.

The Structural Basis of Babesia orientalis Lactate Dehydrogenase

Glycolytic enzymes play a crucial role in the anaerobic glycolysis of apicomplexan parasites for energy generation. Consequently, they are considered as potential targets for new drug development. Previous studies revealed that lactate dehydrogenase (LDH), a glycolytic enzyme, is a potential drug target in different parasites, such as Plasmodium, Toxoplasma, Cryptosporidium, and Piroplasma. Herein, in order to investigate the structural basis of LDH in Babesia spp., we determined the crystal structure of apo Babesia orientalis (Bo) LDH at 2.67-Å resolution in the space group P1. A five-peptide insertion appears in the active pocket loop of BoLDH to create a larger catalytic pocket, like other protozoa (except for Babesia microti LDH) and unlike its mammalian counterparts, and the absence of this extra insertion inactivates BoLDH. Without ligands, the apo BoLDH takes R-state (relaxed) with the active-site loop open. This feature is obviously different from that of allosteric LDHs in T-state (tense) with the active-site loop open. Compared with allosteric LDHs, the extra salt bridges and hydrogen bonds make the subunit interfaces of BoLDH more stable, and that results in the absence of T-state. Interestingly, BoLDH differs significantly from BmLDH, as it exhibits the ability to adapt quickly to the synthetic co-factor APAD+. In addition, the enzymatic activity of BoLDH was inhibited non-competitively by polyphenolic gossypol with a Ki value of 4.25 μM, indicating that BoLDH is sensitive to the inhibition of gossypol and possibly to its new derivative compounds. The current work provides the structural basis of BoLDH for the first time and suggests further investigation on the LDH structure of other Babesia spp. That knowledge would indeed facilitate the screening and designing of new LDH inhibitors to control the intracellular proliferation of Babesia spp.

Structural insights into how vacuolar sorting receptors recognize the sorting determinants of seed storage proteins

In Arabidopsis, vacuolar sorting receptor isoform 1 (VSR1) sorts 12S globulins to the protein storage vacuoles during seed development. Vacuolar sorting is mediated by specific protein–protein interactions between VSR1 and the vacuolar sorting determinant located at the C terminus (ctVSD) on the cargo proteins. Here, we determined the crystal structure of the protease-associated domain of VSR1 (VSR1-PA) in complex with the C-terminal pentapeptide (468RVAAA472) of cruciferin 1, an isoform of 12S globulins. The 468RVA470 motif forms a parallel β-sheet with the switch III residues (127TMD129) of VSR1-PA, and the 471AA472 motif docks to a cradle formed by the cargo-binding loop (95RGDCYF100), making a hydrophobic interaction with Tyr99. The C-terminal carboxyl group of the ctVSD is recognized by forming salt bridges with Arg95. The C-terminal sequences of cruciferin 1 and vicilin-like storage protein 22 were sufficient to redirect the secretory red fluorescent protein (spRFP) to the vacuoles in Arabidopsis protoplasts. Adding a proline residue to the C terminus of the ctVSD and R95M substitution of VSR1 disrupted receptor–cargo interactions in vitro and led to increased secretion of spRFP in Arabidopsis protoplasts. How VSR1-PA recognizes ctVSDs of other storage proteins was modeled. The last three residues of ctVSD prefer hydrophobic residues because they form a hydrophobic cluster with Tyr99 of VSR1-PA. Due to charge–charge interactions, conserved acidic residues, Asp129 and Glu132, around the cargo-binding site should prefer basic residues over acidic ones in the ctVSD. The structural insights gained may be useful in targeting recombinant proteins to the protein storage vacuoles in seeds.

Novel interaction between neurotrophic factor-α1/carboxypeptidase E and serotonin receptor, 5-HTR1E, protects human neurons against oxidative/neuroexcitotoxic stress via β-arrestin/ERK signaling

Protecting neurons from death during oxidative and neuroexcitotoxic stress is key for preventing cognitive dysfunction. We uncovered a novel neuroprotective mechanism involving interaction between neurotrophic factor-α1 (NF-α1/carboxypeptidase E, CPE) and human 5-HTR1E, a G protein-coupled serotonin receptor with no previously known neurological function. Co-immunoprecipitation and pull-down assays confirmed interaction between NFα1/CPE and 5-HTR1E and 125I NF-α1/CPE-binding studies demonstrated saturable, high-affinity binding to 5-HTR1E in stably transfected HEK293 cells (Kd = 13.82 nM). Treatment of 5-HTR1E stable cells with NF-α1/CPE increased pERK 1/2 and pCREB levels which prevented a decrease in pro-survival protein, BCL2, during H2O2-induced oxidative stress. Cell survival assay in β-arrestin Knockout HEK293 cells showed that the NF-α1/CPE-5-HTR1E-mediated protection against oxidative stress was β-arrestin-dependent. Molecular dynamics studies revealed that NF-α1/CPE interacts with 5-HTR1E via 3 salt bridges, stabilized by several hydrogen bonds, independent of the serotonin pocket. Furthermore, after phosphorylating the C-terminal tail and intracellular loop 3 (ICL3) of NF-α1/CPE-5-HTR1E, it recruited β-arrestin1 by forming numerous salt bridges and hydrogen bonds to ICL2 and ICL3, leading to activation of β-arrestin1. Immunofluorescence studies showed 5-HTR1E and NF-α1/CPE are highly expressed and co-localized on cell surface of human hippocampal neurons. Importantly, knock-down of 5-HTR1E in human primary neurons diminished the NF-α1/CPE-mediated protection of these neurons against oxidative stress and glutamate neurotoxicity-induced cell death. Thus, NF-α1/CPE uniquely interacts with serotonin receptor 5-HTR1E to activate the β-arrestin/ERK/CREB/BCL2 pathway to mediate stress-induced neuroprotection.

Capturing a crucial 'disorder-to-order transition' at the heart of the coronavirus molecular pathology – triggered by highly persistent, interchangeable salt-bridges

The COVID-19 origin debate has greatly been influenced by Genome comparison studies of late, revealing the seemingly sudden emergence of the Furin-Like Cleavage Site at the S1/S2 junction of the SARS-CoV-2 Spike (FLCS_Spike) containing its 681_PRRAR_685 motif, absent in other related respiratory viruses. Being the rate-limiting (i.e., the slowest) step, the host Furin cleavage is instrumental in the abrupt increase in transmissibility in COVID-19, compared to earlier onsets of respiratory viral diseases. In such a context, the current paper entraps a disorder-to-order transition of the FLCS_Spike (concomitant to an entropy arrest) upon binding to Furin. The interaction clearly seems to be optimized for a more efficient proteolytic cleavage in SARS-CoV-2. The study further shows the formation of dynamically interchangeable and persistent networks of salt-bridges at the Spike-Furin interface in SARS-CoV-2 involving the three arginines (R682, R683, R685) of the FLCS_Spike with several anionic residues (E230, E236, D259, D264, D306) coming from Furin, strategically distributed around its catalytic triad. Multiplicity and structural degeneracy of plausible salt-bridge network archetypes seems the other key characteristic features of the Spike-Furin binding in SARS-CoV-2 allowing the system to breathe - a trademark of protein disorder transitions. Interestingly, with respect to the homologous interaction in SARS-CoV (2002/2003) taken as a baseline, the Spike-Furin binding events generally in the coronavirus lineage seems to have a preference for ionic bond formation, even with lesser number of cationic residues at their potentially polybasic FLCS_Spike patches. The interaction energies are suggestive of a characteristic metastabilities attributed to Spike-Furin interactions generally to the coronavirus lineage - which appears to be favorable for proteolytic cleavages targeted at flexible protein loops. The current findings not only offer novel mechanistic insights into the coronavirus molecular pathology and evolution but also add substantially to the existing theories of proteolytic cleavages.

Halo-Tolerant And Thermostable Mechanisms of An Endoglucanase From Marine Aspergillus Niger

The cellulase cocktail of marine Aspergillus niger exhibited halo-tolerant and thermostable properties, which is of great potential in industrial application. In order to excavate the single tolerant cellulase components from complex cellulase cocktail, constitutive homologous expression was employed for direct obtainment of the endoglucanase (AnEGL). Enzymatic property study revealed that AnEGL exhibited a property of halo-tolerance and an outstanding thermostability in high salinity environment. Significantly, its activity increased by 29% and the half-life at 65 °C increased by 26.7-fold with the presence of 4.5 M NaCl. Molecular dynamics simulation revealed that Na+ and Cl－ could form salt bridges with charged residues, and then influenced the activity of loops and the stability of substrate binding pocket, which accounted for the halo-tolerance and thermostability. Further, site-specific mutagenesis study proved that the residues Asp95 and Asp99 in the pocket were of great concern for the tolerant properties. The halo-tolerant and thermostable AnEGL was of great value in lignocellulosic utilization and the conjectural mechanisms were of referential significance for other tolerant enzymes.

Improved binding affinity of the Omicron's spike protein with hACE2 receptor is the key factor behind its increased virulence

The new variant of SARS-CoV-2, Omicron, has been quickly spreading in many countries worldwide. Compared to the original virus, Omicron is characterized by several mutations in its genomic region, including spike protein's receptor-binding domain (RBD). We have computationally investigated the interaction between RBD of both wild-type and omicron variants with hACE2 receptor using molecular dynamics and MM-GBSA based binding free energy calculations. The mode of the interaction between Omicron's RBD to the human ACE2 (hACE2) receptor is similar to the original SARS-CoV-2 RBD except for a few key differences. The bind-ing free energy difference shows that the spike protein of Omicron has increased binding affinity for the hACE-2 receptor. The mutated residues in the RBD showed strong interactions with a few amino acid residues of the hACE2. More specifically, strong electrostatic interactions (salt bridges) and hydrogen bonding were observed between R493 and R498 residues of the Omicron RBD with D30/E35 and D38 residues of the hACE2, respectively. Other mutated amino acids in the Omicron RBD, e.g. S496 and H505, also exhibited hydrogen bonding with the hACE2 receptor. The pi-stacking interaction was also observed between tyrosine residues (RBD-Tyr501: hACE2-Tyr41) in the complex, which contributes majorly to binding free energies suggesting this as one of the key interactions stabilizing the complex formation. The structural insights of RBD:hACE2 complex, their binding mode information and residue wise contributions to binding free energy provide insight on the increased transmissibility of Omicron and pave the way to design and optimize novel antiviral agents.

SARS-CoV-2 Omicron Variant: ACE2 Binding, Cryo-EM Structure of Spike Protein-ACE2 Complex and Antibody Evasion

The newly reported Omicron variant is poised to replace Delta as the most rapidly spread SARS-CoV-2 variant across the world. Cryo-EM structural analysis of the Omicron variant spike protein in complex with human ACE2 reveals new salt bridges and hydrogen bonds formed by mutated residues R493, S496 and R498 in the RBD with ACE2. These interactions appear to compensate for other Omicron mutations such as K417N known to reduce ACE2 binding affinity, explaining our finding of similar biochemical ACE2 binding affinities for Delta and Omicron variants. Neutralization assays show that pseudoviruses displaying the Omicron spike protein exhibit increased antibody evasion, with greater evasion observed in sera obtained from unvaccinated convalescent patients as compared to doubly vaccinated individuals (8- vs 3-fold). The retention of strong interactions at the ACE2 interface and the increase in antibody evasion are molecular factors that likely contribute to the increased transmissibility of the Omicron variant.

Molecular Dynamics Simulations of HLA-CW4-Β2M-KIR2DL1 Protein and Homology Modeling of a Complex Associated with Psoriasis Disease (HLA-CW6-Β 2M-KIR2DS1)

Molecular dynamics (MD) simulation was used to study the interactions of two immune proteins of HLA-Cw4-β2m-KIR2DL1 complex with small peptide QYDDAVYKL (nine amino acids) in an aqueous solution. This study aims to gain a detailed information about the conformational changes and the dynamics of the complex. The right parameters and force field for performing the MD simulations that was needed to calibrate the complex structure were determined. The non-bonded interactions (Electrostatic and van der Waals contributions), H-bond formation, and salt bridges between the ligand HLA-Cw4 and the receptor KIR2DL1 were estimated using the obtained MD trajectories. The buried surface area due to binding was calculated to get insight into the causes of specificity of receptor to ligand and explains mutations experiment. The study concluded that β2-microglobulin, one part of the complex, is not directly interacting with the peptide at the groove; therefore, it could be neglected from simulation. Our results showed that β2-microglobulin does not have any significant effect on the dynamics of the 3D-structure of the complex. This project will help in understanding to optimize candidate drug design, a small peptide that disrupts the interaction, for the optimal biological effect.

Study on the Complexation and Release Mechanism of Methylphenidate Hydrochloride Ion Exchange Resin Complex

Since the advent of ion exchange resin, it has been widely used in many fields, including drug delivery systems. The drug binds to the resin through an exchange reaction to form a drug–resin complex, which can gradually release drugs through the exchange of physiological ions in the gastrointestinal tract, to realize functions such as taste masking and regulating release. In this study, the complexes of methylphenidate hydrochloride and Amberlite IRP69 were prepared and evaluated to explore the mechanism of complexation, influencing factors and release mechanism at a molecular level. Firstly, with the properties of the selected complexes, molecular dynamics simulation was innovatively used to find that the intermolecular interaction between drug molecules and ion exchange resin molecules is mainly caused by the stacking effect of π and salt bridges. Secondly, with the drug loading status as an indicator, the factors affecting the compounding process of the drug and resin were explored. Finally, the release mechanism of the drug–resin complex was studied by mathematical model fitting. In summary, a variety of methods were used to study the mechanism of complexation and release between drug and resin, providing a theoretical basis for promoting the marketing of ion−exchange resin−mediated oral preparations.