Substitution of Coxsackievirus A16 VP1 BC and EF Loop Altered the Protective Immune Responses in Chimera Enterovirus A71.

Hand, foot and mouth disease (HFMD) is a childhood disease caused by enterovirus A71 (EV-A71) and coxsackievirus A16 (CV-A16). Capsid loops are important epitopes for EV-A71 and CV-A16. Seven chimeric EV-A71 (ChiE71) involving VP1 BC (45.5% similarity), DE, EF, GH and HI loops, VP2 EF loop and VP3 GH loop (91.3% similarity) were substituted with corresponding CV-A16 loops. Only ChiE71-1-BC, ChiE71-1-EF, ChiE71-1-GH and ChiE71-3-GH were viable. EV-A71 and CV-A16 antiserum neutralized ChiE71-1-BC and ChiE71-1-EF. Mice immunized with inactivated ChiE71 elicited high IgG, IFN-γ, IL-2, IL-4 and IL-10. Neonatal mice receiving passive transfer of WT EV-A71, ChiE71-1-EF and ChiE71-1-BC immune sera had 100%, 80.0% and no survival, respectively, against lethal challenges with EV-A71, suggesting that the substituted CV-A16 loops disrupted EV-A71 immunogenicity. Passive transfer of CV-A16, ChiE71-1-EF and ChiE71-1-BC immune sera provided 40.0%, 20.0% and 42.9% survival, respectively, against CV-A16. One-day-old neonatal mice immunized with WT EV-A71, ChiE71-1-BC, ChiE71-1-EF and CV-A16 achieved 62.5%, 60.0%, 57.1%, and no survival, respectively, after the EV-A71 challenge. Active immunization using CV-A16 provided full protection while WT EV-A71, ChiE71-1-BC and ChiE71-1-EF immunization showed partial cross-protection in CV-A16 lethal challenge with survival rates of 50.0%, 20.0% and 40%, respectively. Disruption of a capsid loop could affect virus immunogenicity, and future vaccine design should include conservation of the enterovirus capsid loops.

Insight into free energy and dynamic cross-correlations of residue for binding affinity of antibody and receptor binding domain SARS-CoV-2.

SARS-CoV-2 virus continues to evolve and mutate causing most of the mutated variants resist to many of the therapeutic monoclonal antibodies (mAbs). Despite several mAbs retained neutralizing capability for Omicron BA.1 and BA.2, reduction in neutralization potency was reported. Hence, effort of searching for mAb that is broader in neutralization breadth without losing the neutralizing ability is continued. MW06 was reported with capability in neutralizing most of the variants of concern (VOC) and it binds to the conserved region (left flank) near epitope mAb sotrovimab (S309). In this study, binding affinity of mAb MW06 and its cocktail formulation with MW05 for receptor binding domain (RBD) SARS-CoV-2 virus was investigated under molecular dynamics simulations (MDs). Binding free energies computed by Molecular Mechanics Generalised Born Surface Area (MM-GBSA) algorithm predicted the binding affinity of MW06 for RBD BA.1 (-53 kcal/mol) as strong as RBD wildtype (-58 kcal/mol) while deterioration was observed for RBD BA.2 (-43 kcal/mol). Alike S309 and MW06, simulated cocktail mAb (MW05 and MW06)-RBD interactions suggested the neutralizing capability of the cocktail formulation for RBD BA.1 and BA.2 reduced. Meanwhile, residue pairs that favour the communication between the mAb and RBD have been identified by decomposing the free energy per pairwise residue basis. Apart from understanding the effects of mutation occurred in the RBD region on human angiotensin-converting enzyme 2 (hACE2) binding, impact of heavily mutated RBD on mAb-RBD interactions was investigated in this study as well. In addition to energetic profile obtained from MDs, plotting the dynamics cross-correlation map of the mAb-RBD complex under elastic network model (ENM) was aimed to understand the cross-correlations between residue fluctuations. It allows simple and rapid analysis on the motions or dynamics of the protein residues of mAbs and RBD in complex. Protein residues having correlated motions are normally part of the structural domains of the protein and their respective motions and protein function are related. Motion of mutated RBD residues and mAb residues was less correlated while their respective interactions energy computed to be higher. The combined techniques of MDs and ENM offered simplicity in understanding dynamics and energy contribution that explain binding affinity of mAb-RBD complexes.

Protein-Protein Interactions: Insight from Molecular Dynamics Simulations and Nanoparticle Tracking Analysis.

Protein-protein interaction plays an essential role in almost all cellular processes and biological functions. Coupling molecular dynamics (MD) simulations and nanoparticle tracking analysis (NTA) assay offered a simple, rapid, and direct approach in monitoring the protein-protein binding process and predicting the binding affinity. Our case study of designed ankyrin repeats proteins (DARPins)-AnkGAG1D4 and the single point mutated AnkGAG1D4-Y56A for HIV-1 capsid protein (CA) were investigated. As reported, AnkGAG1D4 bound with CA for inhibitory activity; however, it lost its inhibitory strength when tyrosine at residue 56 AnkGAG1D4, the most key residue was replaced by alanine (AnkGAG1D4-Y56A). Through NTA, the binding of DARPins and CA was measured by monitoring the increment of the hydrodynamic radius of the AnkGAG1D4-gold conjugated nanoparticles (AnkGAG1D4-GNP) and AnkGAG1D4-Y56A-GNP upon interaction with CA in buffer solution. The size of the AnkGAG1D4-GNP increased when it interacted with CA but not AnkGAG1D4-Y56A-GNP. In addition, a much higher binding free energy (∆GB) of AnkGAG1D4-Y56A (-31 kcal/mol) obtained from MD further suggested affinity for CA completely reduced compared to AnkGAG1D4 (-60 kcal/mol). The possible mechanism of the protein-protein binding was explored in detail by decomposing the binding free energy for crucial residues identification and hydrogen bond analysis.

Dynamics and binding interactions of peptide inhibitors of dengue virus entry.

In this study, we investigate the binding interactions of two synthetic antiviral peptides (DET2 and DET4) on type II dengue virus (DENV2) envelope protein domain III. These two antiviral peptides are designed based on the domain III of the DENV2 envelope protein, which has shown significant inhibition activity in previous studies and can be potentially modified further to be active against all dengue strains. Molecular docking was performed using AutoDock Vina and the best-ranked peptide-domain III complex was further explored using molecular dynamics simulations. Molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) was used to calculate the relative binding free energies and to locate the key residues of peptide-protein interactions. The predicted binding affinity correlated well with the previous experimental studies. DET4 outperformed DET2 and is oriented within the binding site through favorable vdW and electrostatic interactions. Pairwise residue decomposition analysis has revealed several key residues that contribute to the binding of these peptides. Residues in DET2 interact relatively lesser with the domain III compared to DET4. Dynamic cross-correlation analysis showed that both the DET2 and DET4 trigger different dynamic patterns on the domain III. Correlated motions were seen between the residue pairs of DET4 and the binding site while binding of DET2 results in anti-correlated motion on the binding site. This work showcases the use of computational study in elucidating and explaining the experiment observation on an atomic level.

Polysulfonate suramin inhibits Zika virus infection.

Zika virus (ZIKV) is an arthropod-borne flavivirus that causes newborn microcephaly and Guillian-Barré syndrome in adults. No therapeutics are available to treat ZIKV infection or other flaviviruses. In this study, we explored the inhibitory effect of glycosaminoglycans and analogues against ZIKV infection. Highly sulfated heparin, dextran sulfate and suramin significantly inhibited ZIKV infection in Vero cells. De-sulfated heparin analogues lose inhibitory effect, implying that sulfonate groups are critical for viral inhibition. Suramin, an FDA-approved anti-parasitic drug, inhibits ZIKV infection with 3-5 log10 PFU viral reduction with IC50 value of ∼2.5-5 µg/ml (1.93 µM-3.85 µM). A time-of-drug-addition study revealed that suramin remains potent even when administrated at 1-24 hpi. Suramin inhibits ZIKV infection by preventing viral adsorption, entry and replication. Molecular dynamics simulation revealed stronger interaction of suramin with ZIKV NS3 helicase than with the envelope protein. Suramin warrants further investigation as a potential antiviral candidate for ZIKV infection. Heparan sulfate (HS) is a cellular attachment receptor for multiple flaviviruses. However, no direct ZIKV-heparin interaction was observed in heparin-binding analysis, and downregulate or removal of cellular HS with sodium chlorate or heparinase I/III did not inhibit ZIKV infection. This indicates that cell surface HS is not utilized by ZIKV as an attachment receptor.

Inhibition and Larvicidal Activity of Phenylpropanoids from Piper sarmentosum on Acetylcholinesterase against Mosquito Vectors and Their Binding Mode of Interaction.

Aedes aegypti, Aedes albopictus and Culex quinquefasciatus are vectors of dengue fever and West Nile virus diseases. This study was conducted to determine the toxicity, mechanism of action and the binding interaction of three active phenylpropanoids from Piper sarmentosum (Piperaceae) toward late 3rd or early 4th larvae of above vectors. A bioassay guided-fractionation on the hexane extract from the roots of Piper sarmentosum led to the isolation and identification of three active phenylpropanoids; asaricin 1, isoasarone 2 and trans-asarone 3. The current study involved evaluation of the toxicity and acetylcholinesterase (AChE) inhibition of these compounds against Aedes aegypti, Aedes albopictus and Culex quinquefasciatus larvae. Asaricin 1 and isoasarone 2 were highly potent against Aedes aegypti, Aedes albopictus and Culex quinquefasciatus larvae causing up to 100% mortality at ≤ 15 µg/mL concentration. The ovicidal activity of asaricin 1, isoasarone 2 and trans-asarone 3 were evaluated through egg hatching. Asaricin 1 and isoasarone 2 showed potent ovicidal activity. Ovicidal activity for both compounds was up to 95% at 25µg/mL. Asaricin 1 and isoasarone 2 showed strong inhibition on acetylcholinesterase with relative IC50 values of 0.73 to 1.87 µg/mL respectively. These findings coupled with the high AChE inhibition may suggest that asaricin 1 and isoasarone 2 are neuron toxic compounds toward Aedes aegypti, Aedes albopictus and Culex quinquefasciatus. Further computational docking with Autodock Vina elaborates the possible interaction of asaricin 1 and isoasarone 2 with three possible binding sites of AChE which includes catalytic triads (CAS: S238, E367, H480), the peripheral sites (PAS: E72, W271) and anionic binding site (W83). The binding affinity of asaricin 1 and isoasarone 2 were relatively strong with asaricin 1 showed a higher binding affinity in the anionic pocket.