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.

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.

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.

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 mosaic-type trimeric RBD-based COVID-19 vaccine candidate induces potent neutralization against Omicron and other SARS-CoV-2 variants.

Large-scale populations in the world have been vaccinated with COVID-19 vaccines, however, breakthrough infections of SARS-CoV-2 are still growing rapidly due to the emergence of immune-evasive variants, especially Omicron. It is urgent to develop effective broad-spectrum vaccines to better control the pandemic of these variants. Here, we present a mosaic-type trimeric form of spike receptor-binding domain (mos-tri-RBD) as a broad-spectrum vaccine candidate, which carries the key mutations from Omicron and other circulating variants. Tests in rats showed that the designed mos-tri-RBD, whether used alone or as a booster shot, elicited potent cross-neutralizing antibodies against not only Omicron but also other immune-evasive variants. Neutralizing antibody ID50 titers induced by mos-tri-RBD were substantially higher than those elicited by homo-tri-RBD (containing homologous RBDs from prototype strain) or the BIBP inactivated COVID-19 vaccine (BBIBP-CorV). Our study indicates that mos-tri-RBD is highly immunogenic, which may serve as a broad-spectrum vaccine candidate in combating SARS-CoV-2 variants including Omicron.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic continues to pose a serious threat to public health and has so far resulted in over six million deaths worldwide. Mass vaccination programs have reduced the risk of serious illness and death in many people, but the virus continues to persist and circulate in communities across the globe. Furthermore, the current vaccines may be less effective against the new variants of the virus, such as Omicron and Delta, which are continually emerging and evolving. Therefore, it is urgent to develop effective vaccines that can provide broad protection against existing and future forms of SARS-CoV-2. There are several different types of SARS-CoV-2 vaccine, but they all work in a similar way. They contain molecules that induce immune responses in individuals to help the body recognize and more effectively fight SARS-CoV-2 if they happen to encounter it in the future. These immune responses may be so specific that new variants of a virus may not be recognized by them. Therefore, a commonly used strategy for producing vaccines with broad protection is to make multiple vaccines that each targets different variants and then mix them together before administering to patients. Here, Zhang et al. took a different approach by designing a new vaccine candidate against SARS-CoV2 that contained three different versions of part of a SARS-CoV2 protein ­ the so-called spike protein ­ all linked together as one molecule. The different versions of the spike protein fragment were designed to include key features of the fragments found in Omicron and several other SARS-CoV-2 variants. The experiments found that this candidate vaccine elicited a much higher immune response against Omicron and other SARS-CoV-2 variants in rats than an existing SARS-CoV-2 vaccine. It was also effective as a booster shot after a first vaccination with the existing SARS-CoV-2 vaccine. These findings demonstrate that the molecule developed by Zhang et al. induces potent and broad immune responses against different variants of SARS-CoV-2 including Omicron in rats. The next steps following on from this work are to evaluate the safety and immunogenicity of this vaccine candidate in clinical trials. In the future, it may be possible to use a similar approach to develop new broad-spectrum vaccines against other viruses.

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.

Safety and immunogenicity of a hybrid-type vaccine booster in BBIBP-CorV recipients in a randomized phase 2 trial.

NVSI-06-08 is a potential broad-spectrum recombinant COVID-19 vaccine that integrates the antigens from multiple SARS-CoV-2 strains into a single immunogen. Here, we evaluate the safety and immunogenicity of NVSI-06-08 as a heterologous booster dose in BBIBP-CorV recipients in a randomized, double-blind, controlled, phase 2 trial conducted in the United Arab Emirates (NCT05069129). Three groups of healthy adults over 18 years of age (600 participants per group) who have administered two doses of BBIBP-CorV 4-6-month, 7-9-month and >9-month earlier, respectively, are randomized 1:1 to receive either a homologous booster of BBIBP-CorV or a heterologous booster of NVSI-06-08. The incidence of adverse reactions is low, and the overall safety profile is quite similar between two booster regimens. Both Neutralizing and IgG antibodies elicited by NVSI-06-08 booster are significantly higher than those by BBIBP-CorV booster against not only SARS-CoV-2 prototype strain but also multiple variants of concerns (VOCs). Especially, the neutralizing antibody GMT against Omicron variant induced by heterologous NVSI-06-08 booster reaches 367.67, which is substantially greater than that boosted by BBIBP-CorV (GMT: 45.03). In summary, NVSI-06-08 is safe and immunogenic as a booster dose following two doses of BBIBP-CorV, which is immunogenically superior to the homologous boost with another dose of BBIBP-CorV.

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.

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.

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.

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.

Effects of ligand binding on the mechanical stability of protein GB1 studied by steered molecular dynamics simulation.

Regulation of the mechanical properties of proteins plays an important role in many biological processes, and sheds light on the design of biomaterials comprised of protein. At present, strategies to regulate protein mechanical stability focus mainly on direct modulation of the force-bearing region of the protein. Interestingly, the mechanical stability of GB1 can be significantly enhanced by the binding of Fc fragments of human IgG antibody, where the binding site is distant from the force-bearing region of the protein. The mechanism of this long-range allosteric control of protein mechanics is still elusive. In this work, the impact of ligand binding on the mechanical stability of GB1 was investigated using steered molecular dynamics simulation, and a mechanism underlying the enhanced protein mechanical stability is proposed. We found that the external force causes deformation of both force-bearing region and ligand binding site. In other words, there is a long-range coupling between these two regions. The binding of ligand restricts the distortion of the binding site and reduces the deformation of the force-bearing region through a long-range allosteric communication, which thus improves the overall mechanical stability of the protein. The simulation results are very consistent with previous experimental observations. Our studies thus provide atomic-level insights into the mechanical unfolding process of GB1, and explain the impact of ligand binding on the mechanical properties of the protein through long-range allosteric regulation, which should facilitate effective modulation of protein mechanical properties.

Impacts of the charged residues mutation S48E/N62H on the thermostability and unfolding behavior of cold shock protein: insights from molecular dynamics simulation with Go model.

The cold shock protein from the hyperthermophile Thermotoga maritima (Tm-Csp) exhibits significantly higher thermostability than its homologue from the thermophile Bacillus caldolyticus (Bc-Csp). Experimental studies have shown that the electrostatic interactions unique to Tm-Csp are responsible for improving its thermostability. In the present work, the favorable charged residues in Tm-Csp were grafted into Bc-Csp by a double point mutation of S48E/N62H, and the impacts of the mutation on the thermostability and unfolding/folding behavior of Bc-Csp were then investigated by using a modified Go model, in which the electrostatic interactions between charged residues were considered in the model. Our simulation results show that this Tm-Csp-like charged residue mutation can effectively improve the thermostability of Bc-Csp without changing its two-state folding mechanism. Besides that, we also studied the unfolding kinetics and unfolding/folding pathway of the wild-type Bc-Csp and its mutant. It is found that this charged residue mutation obviously enhanced the stability of the C-terminal region of Bc-Csp, which decreases the unfolding rate and changes the unfolding/folding pathway of the protein. Our studies indicate that the thermostability, unfolding kinetics and unfolding/folding pathway of Bc-Csp can be artificially changed by introducing Tm-Csp-like favorable electrostatic interactions into Bc-Csp.

Analysis of conformational motions and related key residue interactions responsible for a specific function of proteins with elastic network model.

Protein collective motions play a critical role in many biochemical processes. How to predict the functional motions and the related key residue interactions in proteins is important for our understanding in the mechanism of the biochemical processes. Normal mode analysis (NMA) of the elastic network model (ENM) is one of the effective approaches to investigate the structure-encoded motions in proteins. However, the motion modes revealed by the conventional NMA approach do not necessarily correspond to a specific function of protein. In the present work, a new analysis method was proposed to identify the motion modes responsible for a specific function of proteins and then predict the key residue interactions involved in the functional motions by using a perturbation approach. In our method, an internal coordinate that accounts for the specific function was introduced, and the Cartesian coordinate space was transformed into the internal/Cartesian space by using linear approximation, where the introduced internal coordinate serves as one of the axes of the coordinate space. NMA of ENM in this internal/Cartesian space was performed and the function-relevant motion modes were identified according to their contributions to the specific function of proteins. Then the key residue interactions important for the functional motions of the protein were predicted as the interactions whose perturbation largely influences the fluctuation along the internal coordinate. Using our proposed methods, the maltose transporter (MalFGK2) from E. Coli was studied. The functional motions and the key residue interactions that are related to the channel-gating function of this protein were successfully identified.

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.