Chimeric virus-like particles of human norovirus constructed by structure-guided epitope grafting elicit cross-reactive immunity against both GI.1 and GII.4 genotypes.

IMPORTANCE: Human norovirus (HuNoV) is highly infectious and can result in severe illnesses in the elderly and children. So far, there is no effective antiviral drug to treat HuNoV infection, and thus, the development of HuNoV vaccines is urgent. However, NoV evolves rapidly, and currently, at least 10 genogroups with numerous genotypes have been found. The genetic diversity of NoV and the lack of cross-protection between different genotypes pose challenges to the development of broadly protective vaccines. In this study, guided by structural alignment between GI.1 and GII.4 HuNoV VP1 proteins, several chimeric-type virus-like particles (VLPs) were designed through surface-exposed loop grafting. Mouse immunization studies show that two of the designed chimeric VLPs induced cross-immunity against both GI.1 and GII.4 HuNoVs. To our knowledge, this is the first designed chimeric VLPs that can induce cross-immune activities across different genogroups of HuNoV, which provides valuable strategies for the development of cross-reactive HuNoV vaccines.

Evolution of the interactions between GII.4 noroviruses and histo-blood group antigens: Insights from experimental and computational studies.

Norovirus (NoV) is the major pathogen causing the outbreaks of the viral gastroenteritis across the world. Among the various genotypes of NoV, GII.4 is the most predominant over the past decades. GII.4 NoVs interact with the histo-blood group antigens (HBGAs) to invade the host cell, and it is believed that the receptor HBGAs may play important roles in selecting the predominate variants by the nature during the evolution of GII.4 NoVs. However, the evolution-induced changes in the HBGA-binding affinity for the GII.4 NoV variants and the mechanism behind the evolution of the NoV-HBGA interactions remain elusive. In the present work, the virus-like particles (VLPs) of the representative GII.4 NoV stains epidemic in the past decades were expressed by using the Hansenula polymorpha yeast expression platform constructed by our laboratory, and then the enzyme linked immunosorbent assay (ELISA)-based HBGA-binding assays as well as the molecular dynamics (MD) simulations combined with the molecular mechanics/generalized born surface area (MMGBSA) calculations were performed to investigate the interactions between various GII.4 strains and different types of HBGAs. The HBGA-binding assays show that for all the studied types of HBGAs, the evolution of GII.4 NoVs results in the increased NoV-HBGA binding affinities, where the early epidemic strains have the lower binding activity and the newly epidemic strains exhibit relative stronger binding intensity. Based on the MD simulation and MMGBSA calculation results, a physical mechanism that accounts for the increased HBGA-binding affinity was proposed. The evolution-involved residue mutations cause the conformational rearrangements of loop-2 (residues 390-396), which result in the narrowing of the receptor-binding pocket and thus tighten the binding of the receptor HBGAs. Our experimental and computational studies are helpful for better understanding the mechanism behind the evolution-induced increasing of HBGA-binding affinity, which may provide useful information for the drug and vaccine designs against GII.4 NoVs.

The functional motions and related key residues behind the uncoating of coxsackievirus A16.

The coxsackievirus A16 (CVA16) is a highly contagious virus that causes the hand, foot, and mouth disease, which seriously threatens the health of children. At present, there are still no available antiviral drugs or effective treatments against the infection of CVA16, and thus it is of great significance to develop anti-CVA16 vaccines. However, the intrinsic uncoating property of the capsid may destroy the neutralizing epitopes and influence its immunogenicity, which hinders the vaccine developments. In the present work, the functional-quantity-based elastic network model analysis method developed by our group was extended to combine with group theory to investigate the uncoating motions of the CVA16 capsid, and then the functionally key residues controlling the uncoating motions were identified by our functional-quantity-based perturbation method. Several motion modes encoded in the topological structure of the capsid were revealed to be responsible for the uncoating of CVA16 particle. These modes predominantly contribute to the fluctuation of the gyration radius of the capsid. Then, by using the perturbation method, four clusters of key sites involved in the uncoating motions were identified, whose perturbations induce significant changes in the fluctuation of the gyration radius. These key residues are mainly located at the 2-fold channels, the quasi 3-fold channels, the bottom of the canyons, and the inter-subunit interfaces around the 3-fold axes. Our studies are helpful for better understanding the uncoating mechanism of the CVA16 capsid and provide potential target sites to prevent the uncoating motions, which is valuable for the vaccine design against CVA16.

Identification of key sites controlling protein functional motions by using elastic network model combined with internal coordinates.

The elastic network model (ENM) is an effective method to extract the intrinsic dynamical properties encoded in protein tertiary structures. We have proposed a new ENM-based analysis method to reveal the motion modes directly responsible for a specific protein function, in which an internal coordinate related to the specific function was introduced to construct the internal/Cartesian hybrid coordinate space. In the present work, the function-related internal coordinates combined with a linear perturbation method were applied to identify the key sites controlling specific protein functional motions. The change in the fluctuations of the internal coordinate in response to residue perturbation was calculated in the hybrid coordinate space by using the linear response theory. The residues with the large fluctuation changes were identified to be the key sites that allosterically control the specific protein function. Two proteins, i.e., human DNA polymerase ß and the chaperonin from Methanococcus maripaludis, were investigated as case studies, in which several collective and local internal coordinates were applied to identify the functionally key residues of these two studied proteins. The calculation results are consistent with the experimental observations. It is found that different collective internal coordinates lead to similar results, where the predicted functionally key sites are located at similar positions in the protein structure. While for the local internal coordinates, the predicted key sites tend to be situated at the region near to the coordinate-involving residues. Our studies provide a starting point for further exploring other function-related internal coordinates for other interesting proteins.

Enhancement of protein mechanical stability: Correlated deformations are handcuffed by ligand binding.

As revealed by previous experiments, protein mechanical stability can be effectively regulated by ligand binding with the binding site distant from the force-bearing region. However, the mechanism for such long-range allosteric control of protein mechanics is still largely unknown. In this work, we use protein topology-based elastic network model (ENM) and all-atomic steered molecular dynamics (SMD) simulations to study the impact of ligand binding on protein mechanical stability in two systems, i.e., GB1 and CheY-binding P2-domain of CheA (CBDCheA). Both ENM and SMD results show that the ligand binding has considerable and negligible effects on the mechanical stability of these two proteins, respectively. These results are consistent with the experimental observations. A physical mechanism for the enhancement of protein mechanical stability was then proposed: the correlated deformations of the force-bearing region and the binding site are handcuffed by the binding of ligand. The handcuff effect suppresses the propagation of internal force in the force-bearing region, thus improving the resistance to the loading force. Our study indicates that ENM method can effectively identify the structure motifs allosterically related to the deformation in the force bearing region, as well as the force propagation pathway within the structure of the studied proteins. Hence, it should be helpful to understand the molecular origin of the different mechanical properties in response to ligand binding for GB1 and CBDCheA.

Allosteric transitions of ATP-binding cassette transporter MsbA studied by the adaptive anisotropic network model.

The transporter MsbA is a kind of multidrug resistance ATP-binding cassette transporter that can transport lipid A, lipopolysaccharides, and some amphipathic drugs from the cytoplasmic to the periplasmic side of the inner membrane. In this work, we explored the allosteric pathway of MsbA from the inward- to outward-facing states during the substrate transport process with the adaptive anisotropic network model. The results suggest that the allosteric transitions proceed in a coupled way. The large-scale closing motions of the nucleotide-binding domains occur first, accompanied with a twisting motion at the same time, which becomes more obvious in middle and later stages, especially for the later. This twisting motion plays an important role for the rearrangement of transmembrane helices and the opening of transmembrane domains on the periplasmic side that mainly take place in middle and later stages respectively. The topological structure plays an important role in the motion correlations above. The conformational changes of nucleotide-binding domains are propagated to the transmembrane domains via the intracellular helices IH1 and IH2. Additionally, the movement of the transmembrane domains proceeds in a nonrigid body, and the two monomers move in a symmetrical way, which is consistent with the symmetrical structure of MsbA. These results are helpful for understanding the transport mechanism of the ATP-binding cassette exporters.

Molecular simulation assisted identification of Ca(2+) binding residues in TMEM16A.

Calcium-activated chloride channels (CaCCs) play vital roles in a variety of physiological processes. Transmembrane protein 16A (TMEM16A) has been confirmed as the molecular counterpart of CaCCs which greatly pushes the molecular insights of CaCCs forward. However, the detailed mechanism of Ca(2+) binding and activating the channel is still obscure. Here, we utilized a combination of computational and electrophysiological approaches to discern the molecular mechanism by which Ca(2+) regulates the gating of TMEM16A channels. The simulation results show that the first intracellular loop serves as a Ca(2+) binding site including D439, E444 and E447. The experimental results indicate that a novel residue, E447, plays key role in Ca(2+) binding. Compared with WT TMEM16A, E447Y produces a 30-fold increase in EC50 of Ca(2+) activation and leads to a 100-fold increase in Ca(2+) concentrations that is needed to fully activate the channel. The following steered molecular dynamic (SMD) simulation data suggests that the mutations at 447 reduce the Ca(2+) dissociation energy. Our results indicated that both the electrical property and the size of the side-chain at residue 447 have significant effects on Ca(2+) dependent gating of TMEM16A.

Domain Motions and Functionally-Key Residues of L-Alanine Dehydrogenase Revealed by an Elastic Network Model.

Mycobacterium tuberculosis L-alanine dehydrogenase (L-MtAlaDH) plays an important role in catalyzing L-alanine to ammonia and pyruvate, which has been considered to be a potential target for tuberculosis treatment. In the present work, the functional domain motions encoded in the structure of L-MtAlaDH were investigated by using the Gaussian network model (GNM) and the anisotropy network model (ANM). The slowest modes for the open-apo and closed-holo structures of the enzyme show that the domain motions have a common hinge axis centered in residues Met133 and Met301. Accompanying the conformational transition, both the 1,4-dihydronicotinamide adenine dinucleotide (NAD)-binding domain (NBD) and the substrate-binding domain (SBD) move in a highly coupled way. The first three slowest modes of ANM exhibit the open-closed, rotation and twist motions of L-MtAlaDH, respectively. The calculation of the fast modes reveals the residues responsible for the stability of the protein, and some of them are involved in the interaction with the ligand. Then, the functionally-important residues relevant to the binding of the ligand were identified by using a thermodynamic method. Our computational results are consistent with the experimental data, which will help us to understand the physical mechanism for the function of L-MtAlaDH.

The Intrinsic Dynamics and Unfolding Process of an Antibody Fab Fragment Revealed by Elastic Network Model.

Antibodies have been increasingly used as pharmaceuticals in clinical treatment. Thermal stability and unfolding process are important properties that must be considered in antibody design. In this paper, the structure-encoded dynamical properties and the unfolding process of the Fab fragment of the phosphocholine-binding antibody McPC603 are investigated by use of the normal mode analysis of Gaussian network model (GNM). Firstly, the temperature factors for the residues of the protein were calculated with GNM and then compared with the experimental measurements. A good result was obtained, which provides the validity for the use of GNM to study the dynamical properties of the protein. Then, with this approach, the mean-square fluctuation (MSF) of the residues, as well as the MSF in the internal distance (MSFID) between all pairwise residues, was calculated to investigate the mobility and flexibility of the protein, respectively. It is found that the mobility and flexibility of the constant regions are higher than those of the variable regions, and the six complementarity-determining regions (CDRs) in the variable regions also exhibit relative large mobility and flexibility. The large amplitude motions of the CDRs are considered to be associated with the immune function of the antibody. In addition, the unfolding process of the protein was simulated by iterative use of the GNM. In our method, only the topology of protein native structure is taken into account, and the protein unfolding process is simulated through breaking the native contacts one by one according to the MSFID values between the residues. It is found that the flexible regions tend to unfold earlier. The sequence of the unfolding events obtained by our method is consistent with the hydrogen-deuterium exchange experimental results. Our studies imply that the unfolding behavior of the Fab fragment of antibody McPc603 is largely determined by the intrinsic dynamics of the protein.

Conformational Motions and Functionally Key Residues for Vitamin B12 Transporter BtuCD-BtuF Revealed by Elastic Network Model with a Function-Related Internal Coordinate.

BtuCD-BtuF from Escherichia coli is a binding protein-dependent adenosine triphosphate (ATP)-binding cassette (ABC) transporter system that uses the energy of ATP hydrolysis to transmit vitamin B12 across cellular membranes. Experimental studies have showed that during the transport cycle, the transporter undergoes conformational transitions between the "inward-facing" and "outward-facing" states, which results in the open-closed motions of the cytoplasmic gate of the transport channel. The opening-closing of the channel gate play critical roles for the function of the transporter, which enables the substrate vitamin B12 to be translocated into the cell. In the present work, the extent of opening of the cytoplasmic gate was chosen as a function-related internal coordinate. Then the mean-square fluctuation of the internal coordinate, as well as the cross-correlation between the displacement of the internal coordinate and the movement of each residue in the protein, were calculated based on the normal mode analysis of the elastic network model to analyze the function-related motions encoded in the structure of the system. In addition, the key residues important for the functional motions of the transporter were predicted by using a perturbation method. In order to facilitate the calculations, the internal coordinate was introduced as one of the axes of the coordinate space and the conventional Cartesian coordinate space was transformed into the internal/Cartesian space with linear approximation. All the calculations were carried out in this internal/Cartesian space. Our method can successfully identify the functional motions and key residues for the transporter BtuCD-BtuF, which are well consistent with the experimental observations.

Allosteric transitions of the maltose transporter studied by an elastic network model.

The maltose transporter from Escherichia coli is one of the ATP-binding cassette (ABC) transporters that utilize the energy from ATP hydrolysis to translocate substrates across cellular membranes. Until 2011, three crystal structures have been determined for maltose transporter at different states in the process of transportation. Here, based on these crystal structures, the allosteric pathway from the resting state (inward-facing) to the catalytic intermediate state (outward-facing) is studied by applying an adaptive anisotropic network model. The results suggest that the allosteric transitions proceed in a coupled way. The closing of the nucleotide-binding domains occurs first, and subsequently this conformational change is propagated to the transmembrane domains (TMD) via the EAA and EAS loops, and then to the maltose-binding protein, which facilitates the translocation of the maltose. It is also found that there exist nonrigid-body and asymmetric movements in the TMD. The cytoplasmic gate may only play the role of allosteric propagation during the transition from the pretranslocation to outward-facing states. In addition, the results show that the movment of the helical subdomain towards the RecA-like subdomain mainly occurs in the earlier stages of the transition. These results can provide some insights into the understanding of the mechanism of ABC transporters.

Study, by use of coarse-grained models, of the functionally crucial residues and allosteric pathway of anesthetic regulation of the Gloeobacter violaceus ligand-gated ion channel.

Although pentameric ligand-gated ion channels (pLGICs) have been found to be the targets of general anesthetics, the mechanism of the effects of anesthetics on pLGICs remains elusive. pLGICs from Gloeobacter violaceus (GLIC) can be inhibited by the anesthetic ketamine. X-ray crystallography has shown that the ketamine binding site is distant from the channel gate of the GLIC. It is still not clear how ketamine controls the function of the GLIC by long-range allosteric regulation. In this work, the functionally crucial residues and allosteric pathway of anesthetic regulation of the GLIC were identified by use of a coarse-grained thermodynamic method developed by our group. In our method, the functionally crucial sites were identified as the residues thermodynamically coupled with binding of ketamine. The results from calculation were highly consistent with experimental data. Our study aids understanding of the mechanism of the anesthetic action of ketamine on the GLIC by long-range allosteric modulation.

Design of hepadnavirus core protein-based chimeric virus-like particles carrying epitopes from respiratory syncytial virus.

Respiratory syncytial virus (RSV) is one of the most important pathogens causing respiratory tract infection in humans, especially in infants and the elderly. The identification and structural resolution of the potent neutralizing epitopes on RSV fusion (F) protein enable an "epitope-focused" vaccine design. However, the display of RSV F epitope II on the surface of the widely-used human hepatitis B virus core antigen (HBcAg) has failed to induce neutralizing antibody response in mice. Here, we used the hepadnavirus core protein (HcAg) from different mammalian hosts as scaffolds to construct chimeric virus-like particles (VLPs) presenting the RSV F epitope II. Mouse immunization showed that different HcAg-based chimeric VLPs elicited significantly different neutralizing antibody responses, among which the HcAg derived from roundleaf bat (RBHcAg) is the most immunogenic. Furthermore, RBHcAg was used as the scaffold platform to present multiple RSV F epitopes, and the immunogenicity was further improved in comparison to that displaying a single epitope II. The designed RBHcAg-based multiple-epitope-presenting VLP formulated with MF59-like adjuvant elicited a potent and balanced Th1/Th2 immune response, and offered substantial protection in mice against the challenge of live RSV A2 virus. The designed chimeric VLPs may serve as the potential starting point for developing epitope-focused vaccines against RSV. Our study also demonstrated that RBHcAg is an effective VLP carrier for presenting foreign epitopes, providing a promising platform for epitope-focused vaccine design.

Detection of persistent organic pollutants binding modes with androgen receptor ligand binding domain by docking and molecular dynamics.

BACKGROUND: Persistent organic pollutants (POPs) are persistent in the environment after release from industrial compounds, combustion productions or pesticides. The exposure of POPs has been related to various reproductive disturbances, such as reduced semen quality, testicular cancer, and imbalanced sex ratio. Among POPs, dichlorodiphenyldichloroethylene (4,4'-DDE) and polychlorinated biphenyls (PCBs) are the most widespread and well-studied compounds. Recent studies have revealed that 4,4'-DDE is an antagonist of androgen receptor (AR). However, the mechanism of the inhibition remains elusive. CB-153 is the most common congener of PCBs, while the action of CB-153 on AR is still under debate. RESULTS: Molecular docking and molecular dynamics (MD) approaches have been employed to study binding modes and inhibition mechanism of 4,4'-DDE and CB-153 against AR ligand binding domain (LBD). Several potential binding sites have been detected and analyzed. One possible binding site is the same binding site of AR natural ligand androgen 5α-dihydrotestosterone (DHT). Another one is on the ligand-dependent transcriptional activation function (AF2) region, which is crucial for the co-activators recruitment. Besides, a novel possible binding site was observed for POPs with low binding free energy with the receptor. Detailed interactions between ligands and the receptor have been represented. The disrupting mechanism of POPs against AR has also been discussed. CONCLUSIONS: POPs disrupt the function of AR through binding to three possible biding sites on AR/LBD. One of them shares the same binding site of natural ligand of AR. Another one is on AF2 region. The third one is in a cleft near N-terminal of the receptor. Significantly, values of binding free energy of POPs with AR/LBD are comparable to that of natural ligand androgen DHT.

Identification of the Intrinsic Motions and Related Key Residues Responsible for the Twofold Channel Opening of Poliovirus Capsid by Using an Elastic Network Model Combined with an Internal Coordinate.

Poliovirus (PV) is an infectious virus that causes poliomyelitis, which seriously threatens the health of children. The release of viral RNA is a key step of PV in host cell infection, and multiple lines of evidence have demonstrated that RNA release is initiated by the opening of the twofold channels of the PV capsid. However, the mechanism that controls the twofold channel opening is still not well understood. In addition, the channel opening motion of the recombinant PV capsid leads to the destruction of predominant neutralizing epitopes and thus hinders the capsid as a vaccine immunogen. Therefore, it is important to identify the intrinsic motions and the related key residues controlling the twofold channel opening for understanding the virus infection mechanism and developing capsid-based vaccines. In the present work, the width of the channel was selected as an internal coordinate directly related to the channel opening, and then the elastic network model (ENM) combined with the group theory were employed to extract the intrinsic motion modes that mostly contribute to the opening of the twofold channels. Our results show that the channel opening predominately induced by the breathing motion and the overall rotation of each protomer in the capsid. Then, an internal coordinate-based perturbation method was used to identify the key residues regulating the twofold channel opening of PV. The calculation results showed that the predicted key residues are mainly located at the twofold axes, the bottom of the canyons and the quasi threefold axes. Our study is helpful for better understanding the twofold channel opening mechanism and provides a potential target for preventing the opening of the channels, which is of great significance for PV vaccine design. The source code of this study is available at https://github.com/SJGLAB/CapsidKeyRes.git.

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.

Safety and immunogenicity of a mosaic vaccine booster against Omicron and other SARS-CoV-2 variants: a randomized phase 2 trial.

An ongoing randomized, double-blind, controlled phase 2 trial was conducted to evaluate the safety and immunogenicity of a mosaic-type recombinant vaccine candidate, named NVSI-06-09, as a booster dose in subjects aged 18 years and older from the United Arab Emirates (UAE), who had administered two or three doses of inactivated vaccine BBIBP-CorV at least 6 months prior to enrollment. The participants were randomly assigned with 1:1 to receive a booster dose of NVSI-06-09 or BBIBP-CorV. The primary outcomes were immunogenicity and safety against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant, and the exploratory outcome was cross-immunogenicity against other circulating strains. Between May 25 and 30, 2022, 516 adults received booster vaccination with 260 in NVSI-06-09 group and 256 in BBIBP-CorV group. Interim results showed a similar safety profile between two booster groups, with low incidence of adverse reactions of grade 1 or 2. For immunogenicity, by day 14 post-booster, the fold rises in neutralizing antibody geometric mean titers (GMTs) from baseline elicited by NVSI-06-09 were remarkably higher than those by BBIBP-CorV against the prototype strain (19.67 vs 4.47-fold), Omicron BA.1.1 (42.35 vs 3.78-fold), BA.2 (25.09 vs 2.91-fold), BA.4 (22.42 vs 2.69-fold), and BA.5 variants (27.06 vs 4.73-fold). Similarly, the neutralizing GMTs boosted by NVSI-06-09 against Beta and Delta variants were also 6.60-fold and 7.17-fold higher than those by BBIBP-CorV. Our findings indicated that a booster dose of NVSI-06-09 was well-tolerated and elicited broad-spectrum neutralizing responses against divergent SARS-CoV-2 variants, including Omicron and its sub-lineages.

A new residue-nucleotide propensity potential with structural information considered for discriminating protein-RNA docking decoys.

Understanding the key factors that influence the preferences of residue-nucleotide interactions in specific protein-RNA interactions has remained a research focus. We propose an effective approach to derive residue-nucleotide propensity potentials through considering both the types of residues and nucleotides, and secondary structure information of proteins and RNAs from the currently largest nonredundant and nonribosomal protein-RNA interaction database. To test the validity of the potentials, we used them to select near-native structures from protein-RNA docking poses. The results show that considering secondary structure information, especially for RNAs, greatly improves the predictive power of pair potentials. The success rate is raised from 50.7 to 65.5% for the top 2000 structures, and the number of cases in which a near-native structure is ranked in top 50 is increased from 7 to 13 out of 17 cases. Furthermore, the exclusion of ribosomes from the database contributes 8.3% to the success rate. In addition, some very interesting findings follow: (i) the protein secondary structure element π-helix is strongly associated with RNA-binding sites; (ii) the nucleotide uracil occurs frequently in the most preferred pairs in which the unpaired and non-Watson-Crick paired uracils are predominant, which is probably significant in evolution. The new residue-nucleotide potentials can be helpful for the progress of protein-RNA docking methods, and for understanding the mechanisms of protein-RNA interactions.

Immunogenicity and safety of NVSI-06-07 as a heterologous booster after priming with BBIBP-CorV: a phase 2 trial.

The increased coronavirus disease 2019 (COVID-19) breakthrough cases pose the need of booster vaccination. We conducted a randomised, double-blinded, controlled, phase 2 trial to assess the immunogenicity and safety of the heterologous prime-boost vaccination with an inactivated COVID-19 vaccine (BBIBP-CorV) followed by a recombinant protein-based vaccine (NVSI-06-07), using homologous boost with BBIBP-CorV as control. Three groups of healthy adults (600 individuals per group) who had completed two-dose BBIBP-CorV vaccinations 1-3 months, 4-6 months and ≥6 months earlier, respectively, were randomly assigned in a 1:1 ratio to receive either NVSI-06-07 or BBIBP-CorV boost. Immunogenicity assays showed that in NVSI-06-07 groups, neutralizing antibody geometric mean titers (GMTs) against the prototype SARS-CoV-2 increased by 21.01-63.85 folds on day 28 after vaccination, whereas only 4.20-16.78 folds of increases were observed in control groups. For Omicron variant, the neutralizing antibody GMT elicited by homologous boost was 37.91 on day 14, however, a significantly higher neutralizing GMT of 292.53 was induced by heterologous booster. Similar results were obtained for other SARS-CoV-2 variants of concerns (VOCs), including Alpha, Beta and Delta. Both heterologous and homologous boosters have a good safety profile. Local and systemic adverse reactions were absent, mild or moderate in most participants, and the overall safety was quite similar between two booster schemes. Our findings indicated that NVSI-06-07 is safe and immunogenic as a heterologous booster in BBIBP-CorV recipients and was immunogenically superior to the homologous booster against not only SARS-CoV-2 prototype strain but also VOCs, including Omicron.

Design of a mutation-integrated trimeric RBD with broad protection against SARS-CoV-2.

The continuous emergence of SARS-CoV-2 variants highlights the need of developing vaccines with broad protection. Here, according to the immune-escape capability and evolutionary convergence, the representative SARS-CoV-2 strains carrying the hotspot mutations were selected. Then, guided by structural and computational analyses, we present a mutation-integrated trimeric form of spike receptor-binding domain (mutI-tri-RBD) as a broadly protective vaccine candidate, which combined heterologous RBDs from different representative strains into a hybrid immunogen and integrated immune-escape hotspots into a single antigen. When compared with a homo-tri-RBD vaccine candidate in the stage of phase II trial, of which all three RBDs are derived from the SARS-CoV-2 prototype strain, mutI-tri-RBD induced significantly higher neutralizing antibody titers against the Delta and Beta variants, and maintained a similar immune response against the prototype strain. Pseudo-virus neutralization assay demonstrated that mutI-tri-RBD also induced broadly strong neutralizing activities against all tested 23 SARS-CoV-2 variants. The in vivo protective capability of mutI-tri-RBD was further validated in hACE2-transgenic mice challenged by the live virus, and the results showed that mutI-tri-RBD provided potent protection not only against the SARS-CoV-2 prototype strain but also against the Delta and Beta variants.