Antigenicity of the Mu (B.1.621) and A.2.5 SARS-CoV-2 Spikes

The rapid emergence of SARS-CoV-2 variants is fueling the recent waves of the COVID-19 pandemic. Here, we assessed ACE2 binding and antigenicity of Mu (B.1.621) and A.2.5 Spikes. Both these variants carry some mutations shared by other emerging variants. Some of the pivotal mutations such as N501Y and E484K in the receptor-binding domain (RBD) detected in B.1.1.7 (Alpha), B.1.351 (Beta) and P.1 (Gamma) are now present within the Mu variant. Similarly, the L452R mutation of B.1.617.2 (Delta) variant is present in A.2.5. In this study, we observed that these Spike variants bound better to the ACE2 receptor in a temperature-dependent manner. Pseudoviral particles bearing the Spike of Mu were similarly neutralized by plasma from vaccinated individuals than those carrying the Beta (B.1.351) and Delta (B.1.617.2) Spikes. Altogether, our results indicate the importance of measuring critical parameters such as ACE2 interaction, plasma recognition and neutralization ability of each emerging variant.

Covalent coupling of Spike’s receptor binding domain to a multimeric carrier produces a high immune response against SARS-CoV-2

The receptor binding domain (RBD) of the Spike protein from SARS-CoV-2 is a promising candidate to develop effective COVID-19 vaccines since it can induce potent neutralizing antibodies. We have previously reported the highly efficient production of RBD in Pichia pastoris, which is structurally similar to the same protein produced in mammalian HEK-293T cells. In this work we designed an RBD multimer with the purpose of increasing its immunogenicity. We produced multimeric particles by a transpeptidation reaction between RBD expressed in P. pastoris and Lumazine Synthase from Brucella abortus (BLS), which is a highly immunogenic and very stable decameric 170 kDa protein. Such particles were used to vaccinate mice with two doses 30 days apart. When the particles ratio of RBD to BLS units was high (6–7 RBD molecules per BLS decamer in average), the humoral immune response was significantly higher than that elicited by RBD alone or by RBD-BLS particles with a lower RBD to BLS ratio (1–2 RBD molecules per BLS decamer). Remarkably, multimeric particles with a high number of RBD copies elicited a high titer of neutralizing IgGs. These results indicate that multimeric particles composed of RBD covalent coupled to BLS possess an advantageous architecture for antigen presentation to the immune system, and therefore enhancing RBD immunogenicity. Thus, multimeric RBD-BLS particles are promising candidates for a protein-based vaccine.

Molecular Aspects of Spike–ACE2 Interaction

A new betacoronavirus (CoV-2) is responsible for the pandemic of severe acute respiratory syndrome (SARS) that began in China at the end of 2019, today known as COronaVIrus Disease 2019 (COVID-19). Subsequent studies confirmed the human angiotensin-converting enzyme 2 (hACE2) as the main cell receptor of spike trimeric glycoprotein, located on the viral envelope, mediating the CoV-2 invasion into the host cells through the receptor-binding domain (RBD) of the spike. Computational analysis of the known experimental 3D structures of spike–ACE2 complexes evidenced distinguishing features in the molecular interactions at the RBD-cell receptor binding interface between CoV-2 and previous CoV-1. The spike represents a key target for drug design as well as an optimal antigen for RNA/viral vector vaccines and monoclonal antibodies in order to maximize prevention and therapy of COVID-19.

Molecular Docking for Ligand-Receptor Binding Process Based on Heterogeneous Computing

Molecular docking aims to predict possible drug candidates for many diseases, and it is computationally intensive. Particularly, in simulating the ligand-receptor binding process, the binding pocket of the receptor is divided into subcubes, and when the ligand is docked into all cubes, there are many molecular docking tasks, which are extremely time-consuming. In this study, we propose a heterogeneous parallel scheme of molecular docking for the binding process of ligand to receptor to accelerate simulating. The parallel scheme includes two layers of parallelism, a coarse-grained layer of parallelism implemented in the message-passing interface (MPI) and a fine-grained layer of parallelism focused on the graphics processing unit (GPU). At the coarse-grain layer of parallelism, a docking task inside one lattice is assigned to one unique MPI process, and a grouped master-slave mode is used to allocate and schedule the tasks. Meanwhile, at the fine-gained layer of parallelism, GPU accelerators undertake the computationally intensive computing of scoring functions and related conformation spatial transformations in a single docking task. The results of the experiments for the ligand-receptor binding process show that on a multicore server with GPUs the parallel program has achieved a speedup ratio as high as 45 times in flexible docking and as high as 54.5 times in semiflexible docking, and on a distributed memory system, the docking time for flexible docking and that for semiflexible docking gradually decrease as the number of nodes used in the parallel program gradually increases. The scalability of the parallel program is also verified in multiple nodes on a distributed memory system and is approximately linear.

Online Hydrophilic Interaction Chromatography (HILIC) Enhanced Top-Down Mass Spectrometry Characterization of SARS-Cov-2 Spike Receptor Binding Domain

SARS-CoV-2 cellular infection is mediated by the heavily glycosylated spike protein. Recombinant versions of the spike protein and the receptor binding domain (RBD) are necessary for seropositivity assays and can potentially serve as vaccines against viral infection. RBD plays key roles in the spike protein’s structure and function, and thus comprehensive characterization of recombinant RBD is critically important for biopharmaceutical applications. Liquid-chromatography coupled to mass spectrometry (LCMS) has been widely used to characterize post-translational modifications in proteins including glycosylation. Most studies of RBDs were performed at the proteolytic peptide (bottom-up proteomics) or released glycan level because of the technical challenges in resolving highly heterogenous glycans at the intact protein level. Herein, we evaluated several online separation techniques: 1. C2 reverse-phase liquid chromatography (RPLC), 2. capillary zone electrophoresis (CZE), and 3. acrylamide-based monolithic hydrophilic interaction chromatography (HILIC) to separate intact recombinant RBDs with varying combinations of glycosylations (glycoforms) for top-down MS. Within the conditions we explored, the HILIC method was superior to RPLC and CZE at separating RBD glycoforms, which differ significantly in neutral glycan groups. In addition, our top-down analysis readily captured unexpected modifications (e.g., cysteinylation, N-terminal sequence variation) and low abundance, heavily glycosylated proteoforms that may be missed by using glycopeptide data alone. The HILIC top-down MS platform holds great potential in resolving heterogenous glycoproteins for facile comparison of biosimilars in quality control applications.

Comparative Immunogenicity of the Recombinant Receptor-Binding Domain of Protein S SARS-CoV-2 Obtained in Prokaryotic and Mammalian Expression Systems

The receptor-binding domain (RBD) of the protein S SARS-CoV-2 is considered to be one of the appealing targets for developing a vaccine against COVID-19. The choice of an expression system is essential when developing subunit vaccines, as it ensures the effective synthesis of the correctly folded target protein, and maintains its antigenic and immunogenic properties. Here, we describe the production of a recombinant RBD protein using prokaryotic (pRBD) and mammalian (mRBD) expression systems, and compare the immunogenicity of prokaryotic and mammalian-expressed RBD using a BALB/c mice model. An analysis of the sera from mice immunized with both variants of the protein revealed that the mRBD expressed in CHO cells provides a significantly stronger humoral immune response compared with the RBD expressed in E.coli cells. A specific antibody titer of sera from mice immunized with mRBD was ten-fold higher than the sera from the mice that received pRBD in ELISA, and about 100-fold higher in a neutralization test. The data obtained suggests that mRBD is capable of inducing neutralizing antibodies against SARS-CoV-2.

Mutational landscape and in silico structure models of SARS-CoV-2 spike receptor binding domain reveal key molecular determinants for virus-host interaction

Background SARS-CoV-2, the causative agent of COVID-19 pandemic is a RNA virus prone to mutations. Formation of a stable binding interface between the Receptor Binding Domain (RBD) of SARS-CoV-2 Spike (S) protein and Angiotensin-Converting Enzyme 2 (ACE2) of host is pivotal for viral entry. RBD has been shown to mutate frequently during pandemic. Although, a few mutations in RBD exhibit enhanced transmission rates leading to rise of new variants of concern, most RBD mutations show sustained ACE2 binding and virus infectivity. Yet, how all these mutations make the binding interface constantly favourable for virus remain enigmatic. This study aims to delineate molecular rearrangements in the binding interface of SARS-CoV-2 RBD mutants. Results Here, we have generated a mutational and structural landscape of SARS-CoV-2 RBD in first six months of the pandemic. We analyzed 31,403 SARS-CoV-2 genomes randomly across the globe, and identified 444 non-synonymous mutations in RBD that cause 49 distinct amino acid substitutions in contact and non-contact amino acid residues. Molecular phylogenetic analysis suggested independent emergence of RBD mutants. Structural mapping of these mutations on the SARS-CoV-2 Wuhan reference strain RBD and structural comparison with RBDs from bat-CoV, SARS-CoV, and pangolin-CoV, all bound to human or mouse ACE2, revealed several changes in the interfacial interactions in all three binding clusters. Interestingly, interactions mediated via N487 residue in cluster-I and Y449, G496, T500, G502 residues in cluster-III remained largely unchanged in all RBD mutants. Further analysis showed that these interactions are evolutionarily conserved in sarbecoviruses which use ACE2 for entry. Importantly, despite extensive changes in the interface, RBD-ACE2 stability and binding affinities were maintained in all the analyzed mutants. Taken together, these findings reveal how SARS-CoV-2 uses its RBD residues to constantly remodel the binding interface. Conclusion Our study broadly signifies understanding virus-host binding interfaces and their alterations during pandemic. Our findings propose a possible interface remodelling mechanism used by SARS-CoV-2 to escape deleterious mutations. Future investigations will focus on functional validation of in-silico findings and on investigating interface remodelling mechanisms across sarbecoviruses. Thus, in long run, this study may provide novel clues to therapeutically target RBD-ACE2 interface for pan-sarbecovirus infections.

Wild Type and Omicron SARS-CoV-2 Spike Receptor Binding Domains Bind Similarly to the Human ACE2 Receptor: An MM-GBSA Study

SARS-CoV-2 is a coronavirus that has created a global pandemic. The virus contains a spike protein which has been shown to bind to the ACE2 receptor on the surface of human cells. Vaccines have been developed that recognize elements of the SARS-CoV-2 spike protein and they have been successful in preventing infection. Recently, the omicron variant of the SARS-CoV-2 virus was reported and quickly became a variant of concern due to its transmissibility. This variant contained an unusually large number (32) of point mutations, of which 15 of those mutations are in the receptor binding domain of the spike protein. In order to assess the differential binding ability of the wild type and omicron variant of the RBD spike protein to human ACE2 receptors, we conducted 2 μs of molecular dynamics simulation to estimate the binding affinities and behaviors. Based upon MM-GBSA binding affinity, center of mass distance measurements, ensemble clustering, pairwise residue decomposition and hydrogen bonding analysis, we can conclude that the 15 point mutations in the receptor binding domain do not increase the affinity of the spike protein for the human ACE2 receptor. The MM-GBSA binding estimations over a 2 μs trajectory, suggest that the wild type binds to ACE2 with a value of -29.69 kcal/mol while the omicron mutant binds with an energy value of -26.67 kcal/mol. These values are within the error estimates of the MM-GBSA method. While some mutations increase binding, more mutations diminish binding, leading to an overall similar picture of binding.

Towards an optimal monoclonal antibody with higher binding affinity to the receptor-binding domain of SARS-CoV-2 spike proteins from different variants

A highly efficient and robust multiple scales in silico protocol, consisting of atomistic constant charge Molecular Dynamics (MD), constant-charge coarse-grain (CG) MD and constant-pH CG Monte Carlo (MC), has been used to study the binding affinities, the free energy of complexation of selected antigen-binding fragments of the monoclonal antibody (mAbs) CR3022 (originally derived from SARS-CoV-1 patients almost two decades ago) and 11 SARS-CoV-2 variants including the wild type. CR3022 binds strongly to the receptor-binding domain (RBD) of SARS-CoV-2 spike protein, but chooses a different site rather than the receptor-binding motif (RBM) of RBD, allowing its combined use with other mAbs against new emerging virus variants. Totally 235,000 mAbs structures were generated using the RosettaAntibodyDesign software, resulting in top 10 scored CR3022-RBD complexes with critical mutations and compared to the native one, all having the potential to block virus-host cell interaction. Of these 10 finalists, two candidates were further identified in the CG simulations to be clearly best against all virus variants, and surprisingly, all 10 candidates and the native CR3022 did exhibit a higher affinity for the Omicron variant with its highest number of mutations (15) of them all considered in this study. The multiscale protocol gives us a powerful rational tool to design efficient mAbs. The electrostatic interactions play a crucial role and appear to be controlling the affinity and complex building. Clearly, mAbs carrying a lower net charge show a higher affinity. Structural determinants could be identified in atomistic simulations and their roles are discussed in detail to further hint at a strategy towards designing the best RBD binder. Although the SARS-CoV-2 was specifically targeted in this work, our approach is generally suitable for many diseases and viral and bacterial pathogens, leukemia, cancer, multiple sclerosis, rheumatoid, arthritis, lupus, and more.