Why does avian influenza A virus hemagglutinin bind to avian receptor stronger than to human receptor? Ab initio fragment molecular orbital studies

Molecular recognition of SARS-CoV-2 spike glycoprotein: quantum chemical hot spot and epitope analyses

Due to the COVID-19 pandemic, researchers have attempted to identify complex structures of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike glycoprotein (S-protein) with angiotensin-converting enzyme 2 (ACE2) or a blocking antibody. However, the molecular recognition mechanism-critical information for drug and antibody design-has not been fully clarified at the amino acid residue level. Elucidating such a microscopic mechanism in detail requires a more accurate molecular interpretation that includes quantum mechanics to quantitatively evaluate hydrogen bonds, XH/π interactions (X = N, O, and C), and salt bridges. In this study, we applied the fragment molecular orbital (FMO) method to characterize the SARS-CoV-2 S-protein binding interactions with not only ACE2 but also the B38 Fab antibody involved in ACE2-inhibitory binding. By analyzing FMO-based interaction energies along a wide range of binding interfaces carefully, we identified amino acid residues critical for molecular recognition between S-protein and ACE2 or B38 Fab antibody. Importantly, hydrophobic residues that are involved in weak interactions such as CH-O hydrogen bond and XH/π interactions, as well as polar residues that construct conspicuous hydrogen bonds, play important roles in molecular recognition and binding ability. Moreover, through these FMO-based analyses, we also clarified novel hot spots and epitopes that had been overlooked in previous studies by structural and molecular mechanical approaches. Altogether, these hot spots/epitopes identified between S-protein and ACE2/B38 Fab antibody may provide useful information for future antibody design, evaluation of the binding property of the SARS-CoV-2 variants including its N501Y, and small or medium drug design against the SARS-CoV-2.

Enhancing erythrocyte-influenza virus specificity by glycan-conjugated gold nanoparticle: Validation of hemagglutination by aptamer and neuraminidases

This study demonstrated the terminated sialo-sugar chains (Neu5Acα2,6Gal and Neu5Acα2,3Gal)-mediated specificity enhancement of influenza virus and chicken red blood cell (RBC) by hemagglutination assay. These glycan chains were immobilized on the gold nanoparticle (GNP) to withhold the higher numbers. With the preliminary optimization, a clear button formation with 0.5% RBC was visualized. On the other hand, intact B/Tokio/53/99 with 750 nM hemagglutinin (HA) displayed a nice hemagglutination. The interference on the specificity of RBC and influenza virus was observed by anti-influenza aptamer at the concentration 31 nM; however, there is no hemagglutination prevention was noticed in the presence of complementary aptamer sequences. Spiking GNP-conjugated Neu5Acα2,6Gal or Neu5Acα2,3Gal or a mixture of these two to the reaction promoted the hemagglutination to 63-folds higher with 12 nM virus, whereas under the same condition the heat-inactivated viruses were lost the hemagglutination. Neuraminidases from Clostridium perfringens and Arthrobacter ureafaciens at 0.0025 neuraminidase units are able to abolish the hemagglutination. Other enzymes, Glycopeptidase F (Elizabethkingia meningoseptica) and Endoglycosidase H (Streptomyces plicatus) did not show the changes with agglutination. Obviously, sialyl-Gal-terminated glycan-conjugated GNP amendment has enhanced the specificity of erythrocyte-influenza virus and able to be controlled by aptamer or neuraminidases.

Energy decomposition analysis approaches and their evaluation on prototypical protein–drug interaction patterns

The partitioning of the interaction energy into chemical components such as electrostatics, polarization, and charge transfer is possible with energy decomposition analysis approaches. We review and evaluate these for biomolecular applications.

Rapid estimation of binding activity of influenza virus hemagglutinin to human and avian receptors

A critical step for avian influenza viruses to infect human hosts and cause epidemics or pandemics is acquisition of the ability of the viral hemagglutinin (HA) to bind to human receptors. However, current global influenza surveillance does not monitor HA binding specificity due to a lack of rapid and reliable assays. Here we report a computational method that uses an effective scoring function to quantify HA-receptor binding activities with high accuracy and speed. Application of this method reveals receptor specificity changes and its temporal relationship with antigenicity changes during the evolution of human H3N2 viruses. The method predicts that two amino acid differences at 222 and 225 between HAs of A/Fujian/411/02 and A/Panama/2007/99 viruses account for their differences in binding to both avian and human receptors; this prediction was verified experimentally. The new computational method could provide an urgently needed tool for rapid and large-scale analysis of HA receptor specificities for global influenza surveillance.

Extending the Power of Quantum Chemistry to Large Systems with the Fragment Molecular Orbital Method

Following the brief review of the modern fragment-based methods and other approaches to perform quantum-mechanical calculations of large systems, the theoretical development of the fragment molecular orbital method (FMO) is covered in detail, with the emphasis on the physical properties, which can be computed with FMO. The FMO-based polarizable continuum model (PCM) for treating the solvent effects in large systems and the pair interaction energy decomposition analysis (PIEDA) are described in some detail, and a range of applications of FMO to biological studies is introduced. The factors determining the relative stability of polypeptide conformers (alpha-helix, beta-turn, and extended form) are elucidated using FMO/PCM and PIEDA, and the interactions in the Trp-cage miniprotein construct (PDB: 1L2Y) are analyzed using PIEDA.

The Fragment Molecular Orbital Method for Geometry Optimizations of Polypeptides and Proteins

The fragment molecular orbital method (FMO) has been used with a large number of wave functions for single-point calculations, and its high accuracy in comparison to ab initio methods has been well established. We have developed the analytic derivative of the electrostatic interaction between far separated fragments and performed a number of restricted Hartree-Fock (RHF) geometry optimizations using FMO and ab initio methods. In particular, the alpha-helix, beta-turn, and extended conformers of a 10-residue polyalanine were studied and the good FMO accuracy was established (the rms deviations for the former two forms were about 0.2 A and for the latter structure about 0.001 A). Met-enkephalin dimer was used as a model for the polypeptide binding and computed at the 3-21G and 6-31G* levels with a similar accuracy achieved; the error in the binding energy predictions (FMO vs ab initio) was 1-3 kcal/mol. Chignolin (PDB: 1uao) and an agonist polypeptide of the erythropoietin receptor protein (emp1) were optimized at the 3-21(+)G level, with the rms deviation from ab initio of about 0.2 A, or 0.5 degrees in terms of bond angles. The effect of solvation on the structure optimization was studied in chignolin and the Trp-cage miniprotein construct (PDB:1l2y), by describing water with TIP3P. The computed structures in gas phase and solution are compared to each other and experiment.

Influenza viral hemagglutinin complicated shape is advantageous to its binding affinity for sialosaccharide receptor

Do the complexity and the bulkiness of a protein affect the affinity between protein and ligand? We attempted to investigate this problem by using ab initio fragment molecular orbital (FMO) method to calculate the binding energy between human influenza viral hemagglutinin (HA) and human oligo-saccharide receptor. We compared the binding energies of 4 different sizes of human A virus HA H3 subtype complexed with human receptor Neu5Ac(alpha2-6)Gal as a model. The full shape receptor binding domain complexed with Neu5Ac(alpha2-6)Gal had the highest binding energy 170.3kcal/mol at the FMO-HF/STO-3G level, which was 52.3kcal/mol higher than that of the smallest domain-receptor complex. These data provide the consideration of the backyard bulkiness beyond the binding site of protein to the protein-ligand stability.