Theoretical analysis of binding specificity of influenza viral hemagglutinin to avian and human receptors based on the fragment molecular orbital method

The hemagglutinin (HA) protein of the influenza virus binds to the host cell receptor in the early stage of viral infection. A change in binding specificity from avian 2-3 to human 2-6 receptor is essential for optimal human-to-human transmission and pandemics. Therefore, it is important to reveal the key factors governing the binding affinity of HA-receptor complex at the molecular level for the understanding and prediction of influenza pandemics. In this work, on the basis of ab initio fragment molecular orbital (FMO) method, we have carried out the interaction energy analysis of HA-receptor complexes to quantitatively elucidate the binding specificity of HAs to avian and human receptors. To discuss the binding property of influenza HA comprehensively, a number of HAs from human H1, swine H1, avian H3 and avian H5 viruses were analyzed. We performed detailed investigations about the interaction patterns of complexes of various HAs and receptor analogues, and revealed that intra-molecular interactions between conserved residues in HA play an important role for HA-receptor binding. These results may provide a hint to understand the role of conserved acidic residues at the receptor binding site which are destabilized by the electrostatic repulsion with sialic acid. The calculated binding energies and interaction patterns between receptor and HAs are consistent with the binding specificities of each HA and thus explain the receptor binding mechanism. The calculated results in the present analysis have provided a number of viewpoints regarding the models for the HA-receptor binding specificity associated with mutated residues. Examples include the role of Glu190 and Gln226 for the binding specificity of H5 HA. Since H5 HA has not yet been adapted to human receptor and the mechanism of the specificity change is unknown, this result is helpful for the prediction of the change in receptor specificity associated with forthcoming possible pandemics.     

One‐Pot Synthesis of N‐Acetyl‐ and N‐Glycolylneuraminic Acid Capped Trisaccharides and Evaluation of Their Influenza A(H1 N1) Inhibition

Human lung epithelial cells natively offer terminal N‐acetylneuraminic acid (Neu5Ac) α(2→6)‐linked to galactose (Gal) as binding sites for influenza virus hemagglutinin. N‐Glycolylneuraminic acid (Neu5Gc) in place of Neu5Ac is known to affect hemagglutinin binding in other species. Not normally generated by humans, Neu5Gc may find its way to human cells from dietary sources. To compare their influence in influenza virus infection, six trisaccharides with Neu5Ac or Neu5Gc α(2→6) linked to Gal and with different reducing end sugar units were prepared using one‐pot assembly and divergent transformation. The sugar assembly made use of an N‐phthaloyl‐protected sialyl imidate for chemoselective activation and α‐stereoselective coupling with a thiogalactoside. Assessment of cytopathic effect showed that the Neu5Gc‐capped trisaccharides inhibited the viral infection better than their Neu5Ac counterparts.     

Monomeric inhibitors of influenza neuraminidase enhance the hemagglutination inhibition activities of polyacrylamides presenting multiple C-sialoside groups

Background: Influenza viruses use hemagglutinin (HA) arrays to bind to sialic acid moieties on the surface of cells; crosslinking of erythrocytes by this mechanism leads to hemagglutination. A number of synthetic polymers containing multiple sialic acid (Neu5Ac) groups as side chains are potent inhibitors of this process. Inhibition may be due to two mechanisms: polyvalent binding of the inhibitor's multiple Neu5Ac side chains to multiple HA sites on the viral surface, or steric stabilization of the viral particle by a layer of the adsorbed, water-swollen polymer, which prevents adhesion to the erythrocyte. The balance between these two effects is not yet known.Results: Polyacrylamides with multiple C-sialosides (PA(Neu5Ac)) were 2–20 fold more effective as inhibitors of virally mediated hemagglutination when assayed in the presence of Neu2en-NH₂, a potent monomeric inhibitor of influenza neuraminidase (NA). The ability of monomeric inhibitors of NA to enhance the inhibition of hemagglutination in this assay correlated with the affinity of the monomer for NA.Conclusions: We propose that inhibitors of NA act by competing with the C-sialosides of PA(Neu5Ac) for binding to the active sites of the NA. Competitive displacement of Neu5Ac causes an expansion of the layer of polymeric gel adsorbed to the virus, enhancing its inhibitory effect. This study provides an example of synergy between two ligands directed toward the active sites of two different proteins, and reinforces the conclusion that steric stabilization is important for the activity of polyvalent inhibitors.     

Intramolecular oximercuration-demercuration reaction: a new stereocontrolled approach to sialic acid containing disaccharides

Methyl 2,3,4-tri-O-benzyl-α,-D-glucopyranoside has been converted into Neu5Ac α(2-6)Glu and Neu5Ac β (2-6)Glu derivatives, using an intramolecular oximercuration-demercuration reaction.     

禽流感病毒HA蛋白与人受体的分子对接

血凝素(hemagglutinin,HA)是位于禽流感病毒表面的糖蛋白。在病毒感染过程中,HA与禽类宿主细胞表面受体结合,介导病毒膜与宿主核内体膜的融合,在传染过程中发挥关键作用。自然界中的禽流感病毒处于不断演化之中,其HA的禽受体结合位点常常发生氨基酸变异。因此,当HA变异体与人受体结合能力较强时,禽流感病毒往往会发生跨种传播而感染人。为预防禽流感的跨种传播,人们迫切需要发展大规模快速检测或预测HA变异体与人受体结合亲和力的方法,以评估各种新发禽流感病毒的跨种传播能力,提前筛选出有潜在危险的病毒株。针对此问题,本研究以H7N9亚型的HA蛋白H7为研究对象,发展了一种运用分子对接的计算方法,预测HA变异体与人受体的结合亲和力。该方法的计算结果表明,H7与人受体的结合亲和力普遍弱于有较强传染人能力的H1,说明H7N9亚型病毒的跨种传播能力普遍较弱;但是,计算分析也揭示,部分新发的H7N9毒株的HA有强的人受体结合亲和力,提示在自然演化过程中,H7N9病毒有可能演化出具有较强的感染人能力的新毒株,这与2013年禽流感疫情的实际发生情况相一致。因此,本文所发展的计算方法可用于快速预测新发禽流感病毒HA与人受体的结合亲和力,为新发禽流感病毒的跨种传播风险评估提供理论依据。     

Synthesis of α- and β-Methyl NEU5Ac α-(2→8) NEU5Ac Disaccharides

ABSTRACT

Acid hydrolysis of colominic acid, an α-(2→8)-linked oligomer of sialic acid, yielded Neu5Ac α-(2→8) Neu5Ac (di-Neu5Ac) 2 as one of the products. Starting from this disaccharide, it was possible to prepare two potential di-Neu5Ac donors, 5 and 8, as their corresponding 2-chloro derivatives. Subsequent reaction of the donor 8 with methanol as a simple acceptor led to the α- and β-methyl Neu5Ac α-(2→8) Neu5Ac glycosides.     

Theoretical investigation on the binding specificity of fluorinated sialyldisaccharides Neu5Acα(2–3)Gal and Neu5Acα(2–6)Gal with influenza hemagglutinin H1 – A Molecular Dynamics Study

In the present study, the hydroxyl groups at the C4 and C7 positions of sialic acid and C6 position of galactose in Neu5Acα(2–3)Gal (N23G) and the hydroxyl groups at the C8 position of sialic acid and C3 and C4 positions of galactose in Neu5Acα(2–6)Gal (N26G) were substituted with fluorine atoms, respectively. Molecular dynamics simulations of 100 ns duration were carried out to investigate the structural and dynamical behavior of H1 bound with the tri-fluorinated N23G and N26G (FN23G and FN26G). Based on energy analysis, it was concluded that FN26G should be a better binder for hemagglutinin (H1) than FN23G and it might act as an inhibitor for influenza.     

Transition state analysis of Vibrio cholerae sialidase-catalyzed hydrolyses of natural substrate analogues

A series of isotopically labeled natural substrate analogues (phenyl 5-N-acetyl-α-d-neuraminyl-(2→3)-β-d-galactopyranosyl-(1→4)-1-thio-β-d-glucopyranoside; Neu5Acα2,3LacβSPh, and the corresponding 2→6 isomer) were prepared chemoenzymatically in order to characterize, by use of multiple kinetic isotope effect (KIE) measurements, the glycosylation transition states for Vibrio cholerae sialidase-catalyzed hydrolysis reactions. The derived KIEs for Neu5Acα2,3LacβSPh for the ring oxygen ((18)V/K), leaving group oxygen ((18)V/K), C3-S deuterium ((D)V/K(S)) and C3-R deuterium ((D)V/K(R)) are 1.029 ± 0.002, 0.983 ± 0.001, 1.034 ± 0.002, and 1.043 ± 0.002, respectively. In addition, the KIEs for Neu5Acα2,6βSPh for C3-S deuterium ((D)V/K(S)) and C3-R deuterium ((D)V/K(R)) are 1.021 ± 0.001 and 1.049 ± 0.001, respectively. The glycosylation transition state structures for both Neu5Acα2,3LacβSPh and Neu5Acα2,6LacβSPh were modeled computationally using the experimental KIE values as goodness of fit criteria. Both transition states are late with largely cleaved glycosidic bonds coupled to pyranosyl ring flattening ((4)H(5) half-chair conformation) with little or no nucleophilic involvement of the enzymatic tyrosine residue. Notably, the transition state for the catalyzed hydrolysis of Neu5Acα2,6βSPh appears to incorporate a lesser degree of general-acid catalysis, relative to the 2,3-isomer.     

Synthesis of spacer armed Kdn(2→6') and (2→3')-lactosamines for immunochemical research

Glycosylation of lactosamine acceptors with Kdn thioglycoside donors in the presence of NIS/TfOH as a promoter affords products with both α- and β-ketosidic linkage (2–6' or 2–3') between the Kdn and Gal residues. After deprotection, the synthesized trisaccharides and glycans containing Neu5Ac were printed to a chip and their comparative interaction with human serum antibodies was explored.     

检测禽流感H5亚型病毒的阻抗型免疫研究

研制了一种可用于H5亚型禽流感病毒快速检测的阻抗型免疫传感器。通过蛋白A将H5N1表面抗原血凝素(HA)的单克隆抗体固定于金叉指阵列微电极表面,并与待测溶液中的目标抗原H5N1进行免疫反应。在[Fe(CN)6]3"/4"溶液中进行电化学阻抗谱扫描,表征电极的表面修饰及抗原捕获过程。当H5N1病毒浓度在21～26 HA unit/50μL范围时,其浓度的对数值与叉指阵列微电极的电子传递阻抗的变化值呈线性关系,相关系数为0.9885;检出限为20 HA unit/50μL,检测时间为1 h。此传感器特异性好,灵敏度高,可以重复使用,在病原微生物快速检测领域具有良好的应用前景。     

Analysis of sialyl oligosaccharides by high-performance liquid chromatography-electrospray ionisation-mass spectrometry with differentiation of alpha2-3 and alpha2-6 sialyl linkages

A method for analysing sialyl oligosaccharides from bovine colostrum using high-performance liquid chromatography-electrospray ionisation-mass spectrometry (HPLC-ESI-MS) is described. Under positive ionisation mode, mass spectra of alpha2-3 and alpha2-6 linkages were different, and the former produced a prominent B2 (or B3 in disialyl lactose) mass fragment. This fragment was absent from mass spectra with alpha2-6 linkages. Two sialyl oligosaccharides, which have not been reported previously, were tentatively identified. One comprises a N-acetyl neuraminic acid (Neu5Ac), two hexoses (Hex), and one N-acetyl hexosamine (HexNAc) residue ((Neu5Ac)1 (Hex)2 (HexNAc)1), and the other comprises one Neu5Ac and one Hex residue ((Neu5Ac)1(Hex)1).     

Elucidating the Anomalous Binding Enhancement of Isoquinoline Boronic Acid for Sialic Acid Under Acidic Conditions: Expanding Biorecognition Beyond Vicinal Diols

A zwitterionic heterocyclic boronic acid based on 4-isoquinolineboronic acid (IQBA) exhibits the highest reported binding affinity for sialic acid or N-acetylneuraminic acid (Neu5Ac, K=5390±190 m ⁻¹) through the formation of a cyclic boronate ester complex under acidic conditions (pH 3). This anomalous pH-dependent binding enhancement does not occur with common neutral saccharides (e.g., glucose, fructose, sorbitiol), because it is mediated via selective complexation to a α-hydroxycarboxylate moiety forming a stable ion pair and ternary complex with Neu5Ac in phosphate buffer. IQBA expands biorecognition beyond classical vicinal diols under neutral or alkaline buffer conditions, which enables the direct analysis of Neu5Ac by native fluorescence with sub-micromolar detection limits.     

Preferential Lectin Binding of Cancer Cells upon Sialic Acid Treatment Under Nutrient Deprivation

The terminal monosaccharide of glycoconjugates on a eukaryotic cell surface is typically a sialic acid (Neu5Ac). Increased sialylation usually indicates progression and poor prognosis of most carcinomas. Here, we utilize two human mammary epithelial cell lines, HB4A (breast normal cells) and T47D (breast cancer cells), as a model system to demonstrate differential surface glycans when treated with sialic acid under nutrient deprivation. Under a starved condition, sialic acid treatment of both cells resulted in increased activities of α2→3/6 sialyltransferases as demonstrated by solid phase assay using lectin binding. However, a very strong Maackia amurensis agglutinin I (MAL-I) staining on the membrane of sialic acid-treated T47D cells was observed, indicating an increase of Neu5Acα2→3Gal on the cell surface. To our knowledge, this is a first report showing the utility of lectins, particularly MAL-I, as a means to discriminate between normal and cancer cells after sialic acid treatment under nutrient deprivation. This method is sensitive and allows selective detection of glycan sialylation on a cancer cell surface.     

3 beta-Hydroxysialic acid glycosides. I. Calcium-binding ability and chemical and enzymatic stabilities.

Methyl alpha- and beta-glycosides of N-acetylneuraminic acid (Neu5Ac) and N-acetyl-3 beta-hydroxyneuraminic acid (Neu5Ac beta 3OH) (1-4) were prepared to evaluate their calcium-binding ability. (Methyl alpha-glucopyranosidonyl) alpha- and beta-, and 4-methylumbelliferyl alpha-glycosides of Neu5Ac and Neu5Ac beta 3OH (5-10) were also synthesized for the comparison of chemical and enzymatic stabilities, respectively. Methyl beta-glycosides of Neu5Ac and Neu5Ac beta 3OH, 3 and 4, respectively, showed intense calcium-binding abilities, while no such ability was observed in the corresponding alpha-glycosides, 1 and 2. The Neu5Ac beta 3OH glycosides, 6, 8, and 10, showed much stronger resistance to acidic hydrolysis and sialidase digestion than the corresponding Neu5Ac glycosides, 5, 7, and 9.     

Differential receptor binding affinities of influenza hemagglutinins on glycan arrays

A library of 27 sialosides, including seventeen 2,3-linked and ten 2,6-linked glycans, has been prepared to construct a glycan array and used to profile the binding specificity of different influenza hemagglutinins (HA) subtypes, especially from the 2009 swine-originated H1N1 and seasonal influenza viruses. It was found that the HAs from the 2009 H1N1 and the seasonal Brisbane strain share similar binding profiles yet different binding affinities toward various α2,6 sialosides. Analysis of the binding profiles of different HA subtypes indicate that a minimum set of 5 oligosaccharides can be used to differentiate influenza H1, H3, H5, H7, and H9 subtypes. In addition, the glycan array was used to profile the binding pattern of different influenza viruses. It was found that most binding patterns of viruses and HA proteins are similar and that glycosylation at Asn27 is essential for receptor binding.     

Characterization of derivatization of sialic acid with 2-aminoacridone and determination of sialic acid content in glycoproteins by capillary electrophoresis and high performance liquid chromatography--ion trap mass spectrometry.

A simple and highly sensitive capillary electrophoresis (CE) method for determining the content of N-acetylneuraminic acid (Neu5Ac) in glycoproteins was developed. Neu5Ac was derivatized with 2-aminoacridone (AMAC) by reductive amination, and the AMAC-Neu5Ac adduct could be readily separated from the other 11 AMAC-derivatized neutral and acidic monosaccharides usually present in glycoproteins by CE in a 0.3 mol/L borate buffer, pH 10.5, and detected at 260 nm. The derivatization of Neu5Ac was achieved at 55 degrees C for 4 h. AMAC-Neu5Ac was stable at 20 degrees C in the dark for at least 12 h while at room temperature it spontaneously converted into another substance with a lower electrophoretic mobility, which was identified as decarboxylated AMAC-Neu5Ac by high performance liquid chromatography - ion trap mass spectrometry (HPLC-ITMS). Concentration and mass of Neu5Ac as low as 1 micromol/L and 35 fmol could be detected. The linear correlation coefficient between the ratio of peak area to migration time of AMAC-Neu5Ac and the concentration of Neu5Ac ranging from 10 to 120 micromol/L was 0.9978 (n=8). This method was successfully applied to the analysis of sialic acid in human urinary trypsin inhibitor (hu-UTI), bovine alpha1-acid glycoprotein (alpha1-AGP) and recombinant human erythropoietin (rhu-EPO). By combination of CE and HPLC-ITMS we found that N-glycolylneuraminic acid (Neu5Gc) was present in bovine alpha1-AGP in addition to Neu5Ac, with a quantity comparable to that of the latter.     

Binding interaction analysis of the active site and its inhibitors for neuraminidase (N1 subtype) of human influenza virus by the integration of molecular docking, FMO calculation and 3D-QSAR CoMFA modeling

Recently, the worldwide spread of A/H5N1 avian influenza with high virulence has highlighted the potential threat of human influenza pandemic. Tamiflu and Relenza are currently the only two anti-influenza drugs targeting the neuraminidase (NA) enzyme of human influenza virus. Reports of the emergence of drug resistance further make the development of new potent anti-influenza inhibitors a priority. The X-ray crystallographic study of A/H5N1 avian influenza NA subtypes (Russell, R. J. Nature 2006, 443, 45-49) has demonstrated that there exist two genetically distinct groups, group-1 (N1, N4, N5 and N8) and group-2 (N2, N3, N6, N7 and N9), whose conformations are substantially different. The detailed comparison of their active sites has established, heretofore, the most accurate and solid molecular basis of structure and mechanism for the development of new anti-influenza drugs. In the present study, a three-dimensional structure of N1 subtype of human influenza type A virus (N1hA) has been generated by homology modeling using the X-ray crystallographic structure of N1 subtype of avian influenza virus (N1aA) as the template. Binding interaction analysis between the active site and its inhibitors has been performed by combining ab initio fragment molecular orbital (FMO) calculations and three-dimensional quantitative structure-activity relationship with comparative molecular field analysis (3D-QSAR CoMFA) modeling. Integrated with docking-based 3D-QSAR CoMFA modeling, molecular surface property (electrostatic and steric) mapping and FMO pair interaction analysis, a set of new receptor-ligand binding models and bioaffinity predictive models for rational design and virtual screening of more potent inhibitors of N1hA are established. In addition, the flexibility of the loop-150 of N1hA and N1aA has been examined by a series of molecular dynamics simulations.     

Synthesis of selective inhibitors against V. cholerae sialidase and human cytosolic sialidase NEU2

Sialidases or neuraminidases catalyze the hydrolysis of terminal sialic acid residues from sialyl oligosaccharides and glycoconjugates. Despite successes in developing potent inhibitors specifically against influenza virus neuraminidases, the progress in designing and synthesizing selective inhibitors against bacterial and human sialidases has been slow. Guided by sialidase substrate specificity studies and sialidase crystal structural analysis, a number of 2-deoxy-2,3-dehydro-N-acetylneuraminic acid (DANA or Neu5Ac2en) analogues with modifications at C9 or at both C5 and C9 were synthesized. Inhibition studies of various bacterial sialidases and human cytosolic sialidase NEU2 revealed that Neu5Gc9N(3)2en and Neu5AcN(3)9N(3)2en are selective inhibitors against V. cholerae sialidase and human NEU2, respectively.     

A Water Soluble Pd2L4 Cage for Selective Binding of Neu5Ac

The sialic acid N-acetylneuraminic acid (Neu5Ac) and its derivatives are involved in many biological processes including cell-cell recognition and infection by influenza. Molecules that can recognize Neu5Ac might thus be exploited to intervene in or monitor such events. A key obstacle in this development is the sparse availability of easily prepared molecules that bind to this carbohydrate in its natural solvent; water. Here, we report that the carbohydrate binding pocket of an organic soluble [Pd₂L₄]⁴⁺ cage could be equipped with guanidinium-terminating dendrons to give the water soluble [Pd₂L₄][NO₃]₁₆ cage 7 . It was shown by means of NMR spectroscopy that 7 binds selectively to anionic monosaccharides and strongest to Neu5Ac with K_a=24 M⁻¹. The cage had low to no affinity for the thirteen neutral saccharides studied. Aided by molecular modeling, the selectivity for anionic carbohydrates such as Neu5Ac could be rationalized by the presence of charge assisted hydrogen bonds and/or the presence of a salt bridge with a guanidinium solubilizing arm of 7 . Establishing that a simple coordination cage such as 7 can already selectively bind to Neu5Ac in water paves the way to improve the stability, affinity and/or selectivity properties of M₂L₄ cages for carbohydrates and other small molecules.     

Anomalous binding profile of phenylboronic acid with N-acetylneuraminic acid (Neu5Ac) in aqueous solution with varying pH

Borates are known to interact with carbohydrate moieties expressed on the surface of biological membranes of a variety of cells, viruses, bacteria, and fungi. This study revealed the anomalous binding profile of borate in aqueous solution with N-acetylneuraminic acid (Neu5Ac, sialic acid) as a potential receptor site on the surfaces of biological membranes using (11)B, (1)H, (13)C, and (15)N nuclear magnetic resonance spectroscopies. 3-(Propionamido)phenylboronic acid (PAPBA) was chosen as the model borate compound. The equilibrium constant (K) for Neu5Ac binding to PAPBA was compared with those for glucose, mannose, and galactose, which are the major carbohydrate constituents of glycoproteins and glycolipids expressed on biological membranes. In the Neu5Ac/PAPBA system, the unusual pH dependency of the K values, a decrease in K with increasing pH, was observed, suggesting the formation of a trigonal-formed complex stabilized by the coordination of an amide group of Neu5Ac at the C-5 position to the boron atom, forming intramolecular B-N or B-O bonding. Furthermore, the anomalously high complexing ability at physiological pH 7.4 was confirmed for this system, with the K value 37.6 which is approximately 7 times higher than that for glucose. This exceptionally high value of K at physiological pH, compared to those of other sugars, strongly suggests that the boronic acid selectively recognizes the Neu5Ac residues of the glycosylated components including glycoproteins and gangliosides existing on the surface of the biological membranes.