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.

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.

E484K mutation in SARS-CoV-2 RBD enhances binding affinity with hACE2 but reduces interactions with neutralizing antibodies and nanobodies: Binding free energy calculation studies.

The pandemic of the COVID-19 disease caused by SARS-CoV-2 has led to more than 200 million infections and over 4 million deaths worldwide. The progress in the developments of effective vaccines and neutralizing antibody therapeutics brings hopes to eliminate the threat of COVID-19. However, SARS-CoV-2 continues to mutate, and several new variants have been emerged. Among the various naturally-occurring mutations, the E484K mutation shared by many variants attracted serious concerns, which may potentially enhance the receptor binding affinity and reduce the immune response. In the present study, the molecular mechanism behind the impacts of E484K mutation on the binding affinity of the receptor-binding domain (RBD) with the receptor human angiotensin-converting enzyme 2 (hACE2) was investigated by using the molecular dynamics (MD) simulations combined with the molecular mechanics-generalized Born surface area (MMGBSA) method. Our results indicate that the E484K mutation results in more favorable electrostatic interactions compensating the burial of the charged and polar groups upon the binding of RBD with hACE2, which significantly improves the RBD-hACE2 binding affinity. Besides that, the E484K mutation also causes the conformational rearrangements of the loop region containing the mutant residue, which leads to tighter binding interface of RBD with hACE2 and formation of some new hydrogen bonds. The tighter binding interface and the new hydrogen bonds formation also contribute to the improved binding affinity of RBD to the receptor hACE2. In addition, six neutralizing antibodies and nanobodies complexed with RBD were selected to explore the effects of E484K mutation on the recognition of these antibodies to RBD. The simulation results show that the E484K mutation significantly reduces the binding affinities to RBD for most of the studied neutralizing antibodies/nanobodies, and the decrease in the binding affinities is mainly owing to the unfavorable electrostatic interactions caused by the mutation. Our studies revealed that the E484K mutation may improve the binding affinity between RBD and the receptor hACE2, implying more transmissibility of the E484K-containing variants, and weaken the binding affinities between RBD and the studied neutralizing antibodies/nanobodies, indicating reduced effectiveness of these antibodies/nanobodies. Our results provide valuable information for the effective vaccine development and antibody/nanobody drug design.

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.

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.

Energy transport pathway in proteins: Insights from non-equilibrium molecular dynamics with elastic network model.

Intra-molecular energy transport between distant functional sites plays important roles in allosterically regulating the biochemical activity of proteins. How to identify the specific intra-molecular signaling pathway from protein tertiary structure remains a challenging problem. In the present work, a non-equilibrium dynamics method based on the elastic network model (ENM) was proposed to simulate the energy propagation process and identify the specific signaling pathways within proteins. In this method, a given residue was perturbed and the propagation of energy was simulated by non-equilibrium dynamics in the normal modes space of ENM. After that, the simulation results were transformed from the normal modes space to the Cartesian coordinate space to identify the intra-protein energy transduction pathways. The proposed method was applied to myosin and the third PDZ domain (PDZ3) of PSD-95 as case studies. For myosin, two signaling pathways were identified, which mediate the energy transductions form the nucleotide binding site to the 50 kDa cleft and the converter subdomain, respectively. For PDZ3, one specific signaling pathway was identified, through which the intra-protein energy was transduced from ligand binding site to the distant opposite side of the protein. It is also found that comparing with the commonly used cross-correlation analysis method, the proposed method can identify the anisotropic energy transduction pathways more effectively.