A Pro to Gly mutation in the hinge of the arabinose-binding protein enhances binding and alters specificity. Sugar-binding and crystallographic studies

The L-arabinose-binding protein (ABP) of Escherichia coli consists structurally of two distinct globular domains connected by a hinge of three separate peptide segments. Arabinose is bound and completely sequestered within the deep cleft between the two domains. With reduced affinity, ABP also binds D-galactose (approximately 2-fold reduction) and D-fucose (approximately 40-fold reduction). Experiments have been conducted to explore the role in sugar binding of the hinge connecting the two domains of ABP. To increase the flexibility of the hinge region, a glycine was substituted for a proline at position 254 by site-directed mutagenesis. Unexpectedly, this mutation resulted in the dramatic enhancement of galactose binding over that of arabinose. The affinity of the mutant ABP for galactose increased by over 20-fold, while that for arabinose and fucose remained relatively unchanged. We have measured association and dissociation rates of the Gly-254 ABP with L-arabinose, D-galactose, and D-fucose and have determined the crystallographic structure of the protein complexed with each of the three sugars. Both the ligand-binding kinetic measurements and structure analysis indicate that the altered specificity is due to an effective increase in the rigidity of the hinge in the closed conformation which is induced upon galactose binding. Stabilizing contacts are formed between the strands of the hinge in the Gly-254 ABP when galactose is bound which are not found in complexes with the other sugars or the liganded wild-type protein.     

Computational analysis of the binding affinities of the natural-product cyclopentapeptides argifin and argadin to chitinase B from Serratia marcescens

Molecular dynamics (MD) simulations and the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method were applied to study the interaction of the natural-product cyclopentapeptide chitinase inhibitors argifin and argadin with chitinase B (ChiB) from Serratia marcescens. Argadin inhibited ChiB with an inhibition constant (K(i)) value of 20 nM, which was three orders of magnitude greater than that of argifin (K(i)=33,000 nM). The MM-PBSA free-energy analysis provided absolute binding free energies of -6.98 and -11.16 kcal/mol for the argifin and argadin complexes, respectively. These estimates were in good agreement with the free energies derived from the experimental K(i) values (-6.36 and -10.92 kcal/mol for the argifin and argadin complexes, respectively). The energetic analysis revealed that the van der Waals and nonpolar solvation energies drove the binding of both argifin and argadin. We found that the binding of argadin gained approximately 12 kcal/mol more van der Waals energy than that of argifin, which was mainly responsible for the difference in binding free energy between argifin and argadin. In particular, W220 and W403 of ChiB were found to contribute to the more favorable van der Waals interaction with argadin. We also designed argifin derivatives with better binding affinity, in which a constituent amino-acid residue of argifin was mutated to one with a bulky side chain. The derivative in which D-Ala of argifin was replaced with D-Trp appeared to possess a binding affinity that was equally potent to that of argadin.     

Interactions of wheat-germ agglutinin with GlcNAc beta 1,6Gal sequence

The interactions of wheat-germ agglutinin (WGA) with the GlcNAc beta 1,6Gal sequence, a characteristic component of branched poly-N-acetyllactosaminoglycans, were investigated using isothermal titration calorimetry and X-ray crystallography. GlcNAc beta 1,6Gal exhibited an affinity greater than GlcNAc beta 1,4GlcNAc to all WGA isolectins, whereas Gal beta 1,6GlcNAc showed much less affinity than GlcNAc beta 1,4GlcNAc. X-ray structural analyses of the glutaraldehyde-crosslinked WGA isolectin 3 crystals in complex with GlcNAc beta 1,6Gal, GlcNAc beta 1,4GlcNAc and GlcNAc beta 1,6Gal beta 1,4Glc were performed at 2.4, 2.2 and 2.2 A resolution, respectively. In spite of different glycosidic linkages, GlcNAc beta 1,6Gal and GlcNAc beta 1,4GlcNAc exhibited basically similar binding modes to each other, in contact with side chains of two aromatic residues, Tyr64 and His66. However, the conformations of the ligands in the two primary binding sites were not always identical. GlcNAc beta 1,6Gal showed more extensive variation in the parameters defining the glycosidic linkage structure compared to GlcNAc beta 1,4GlcNAc, demonstrating large conformational flexibility of the former ligand in the interaction with WGA. The difference in the ligand binding conformation was accompanied by alterations of the side chain conformation of the amino acid residues involved in the interactions. The hydrogen bond between Ser62 and the non-reducing end GlcNAc was always observed regardless of the ligand type, indicating the key role of this interaction. In addition to the hydrogen bonding and van der Waals interactions, CH--pi interactions involving Tyr64, His66 and Tyr73 are suggested to play an essential role in determining the ligand binding conformation in all complexes. One of the GlcNAc beta 1,6Gal ligands had no crystal packing contact with another WGA molecule, therefore the conformation might be more relevant to the interaction mode in solution.     

Statistical analysis of atomic contacts at RNA-protein interfaces

Forty-five crystals of complexes between proteins and RNA molecules from the Protein Data Bank have been statistically surveyed for the number of contacts between RNA components (phosphate, ribose and the four bases) and amino acid side chains. Three groups of complexes were defined: the tRNA synthetases; the ribosomal complexes; and a third group containing a variety of complexes. The types of atomic contacts were a priori classified into ionic, neutral H-bond, C-H...O H-bond, or van der Waals interaction. All the contacts were organized into a relational database which allows for statistical analysis. The main conclusions are the following: (i) in all three groups of complexes, the most preferred amino acids (Arg, Asn, Ser, Lys) and the less preferred ones (Ala, Ile, Leu, Val) are the same; Trp and Cys are rarely observed (respectively 15 and 5 amino acids in the ensemble of interfaces); (ii) of the total number of amino acids located at the interfaces 22% are hydrophobic, 40% charged (positive 32%, negative 8%), 30% polar and 8% are Gly; (iii) in ribosomal complexes, phosphate is preferred over ribose, which is preferred over the bases, but there is no significant preference in the other two groups; (iv) there is no significant prevalence of a base type at protein-RNA interfaces, but specifically Arg and Lys display a preference for phosphate over ribose and bases; Pro and Asn prefer bases over ribose and phosphate; Met, Phe and Tyr prefer ribose over phosphate and bases. Further, Ile, Pro, Ser prefer A over the others; Leu prefers C; Asp and Gly prefer G; and Asn prefers U. Considering the contact types, the following conclusions could be drawn: (i) 23% of the contacts are via potential H-bonds (including CH...O H-bonds and ionic interactions), 72% belong to van der Waals interactions and 5% are considered as short contacts; (ii) of all potential H-bonds, 54% are standard, 33% are of the C-H...O type and 13% are ionic; (iii) the Watson-Crick sites of G, O6(G) and principally N2(G) and the hydroxyl group O2' is more often involved in H-bonds than expected; the protein main chain is involved in 32% and the side chains in 68% of the H-bonds; considering the neutral and ionic H-bonds, the following couples are more frequent than expected-base A-Ser, base G-Asp/Glu, base U-Asn. The RNA CH groups interact preferentially with oxygen atoms (62% on the main chain and 19% on the side chains); (iv) the bases are involved in 38% of all H-bonds and more than 26% of the H-bonds have the H donor group on the RNA; (v) the atom O2' is involved in 21% of all H-bonds, a number greater than expected; (vi) amino acids less frequently in direct contact with RNA components interact frequently via their main chain atoms through water molecules with RNA atoms; in contrast, those frequently observed in direct contact, except Ser, use instead their side chain atoms for water bridging interactions.     

Protein-ligand energetics assessed using deoxy and fluorodeoxy sugars in equilibrium binding and high resolution crystallographic studies.

Hydrogen bonding interactions are one of the most important single factors in protein-ligand interactions and molecular recognition. To probe the energetics of the interactions, we have analyzed the binding of 1-deoxy-, 2-deoxy- and 6-fluoro-6-deoxy- analogues of D-galactose (Gal) to a primary high-affinity periplasmic receptor for monosaccharide active transport. Kd values and atomic structures refined at 1.81 to 1.45 A resolution of the complexes have been determined and compared with those of Gal binding. With binding site residues and the bound modified sugars in nearly identical positions as found in the complex with Gal, the binding of 1-deoxy-Gal or 2-deoxy-Gal reflects the overall contribution of 1.8 kcal mol-1 per hydrogen bond (neutral-charge type) to the affinity of Gal. Neglected in these estimates is the contribution of van der Waals' forces that accompany the formation of hydrogen bonds with each sugar hydroxyl. Contrary to expectations, the 6-fluoro-6-deoxy analogue proved to be an inadequate probe of Gal OH6 as a hydrogen bond donor due to the binding of a new water molecule and structural changes arising from the electronegative fluoro group. This study sheds new light on the energetics of protein-ligand interactions and the use of engineered ligands in assessing these interactions.     

Protein-RNA interactions: structural analysis and functional classes

A data set of 89 protein-RNA complexes has been extracted from the Protein Data Bank, and the nucleic acid recognition sites characterized through direct contacts, accessible surface area, and secondary structure motifs. The differences between RNA recognition sites that bind to RNAs in functional classes has also been analyzed. Analysis of the complete data set revealed that van der Waals interactions are more numerous than hydrogen bonds and the contacts made to the nucleic acid backbone occur more frequently than specific contacts to nucleotide bases. Of the base-specific contacts that were observed, contacts to guanine and adenine occurred most frequently. The most favored amino acid-nucleotide pairings observed were lysine-phosphate, tyrosine-uracil, arginine-phosphate, phenylalanine-adenine and tryptophan-guanine. The amino acid propensities showed that positively charged and polar residues were favored as expected, but also so were tryptophan and glycine. The propensities calculated for the functional classes showed trends similar to those observed for the complete data set. However, the analysis of hydrogen bond and van der Waal contacts showed that in general proteins complexed with messenger RNA, transfer RNA and viral RNA have more base specific contacts and less backbone contacts than expected, while proteins complexed with ribosomal RNA have less base-specific contacts than the expected. Hence, whilst the types of amino acids involved in the interfaces are similar, the distribution of specific contacts is dependent upon the functional class of the RNA bound.     

Topography of the combining region of a Thomsen-Friedenreich-antigen-specific lectin jacalin (Artocarpus integrifolia agglutinin). A thermodynamic and circular-dichroism spectroscopic study.

Thermodynamic analysis of carbohydrate binding by Artocarpus integrifolia (jackfruit) agglutinin (jacalin) shows that, among monosaccharides, Me alpha GalNAc (methyl-alpha-N-acetylgalactosamine) is the strongest binding ligand. Despite its strong affinity for Me alpha GalNAc and Me alpha Gal, the lectin binds very poorly when Gal and GalNAc are in alpha-linkage with other sugars such as in A- and B-blood-group trisaccharides, Gal alpha 1-3Gal and Gal alpha 1-4Gal. These binding properties are explained by considering the thermodynamic parameters in conjunction with the minimum energy conformations of these sugars. It binds to Gal beta 1-3GalNAc alpha Me with 2800-fold stronger affinity over Gal beta 1-3GalNAc beta Me. It does not bind to asialo-GM1 (monosialoganglioside) oligosaccharide. Moreover, it binds to Gal beta 1-3GalNAc alpha Ser, the authentic T (Thomsen-Friedenreich)-antigen, with about 2.5-fold greater affinity as compared with Gal beta 1-3GalNAc. Asialoglycophorin A was found to be about 169,333 times stronger an inhibitor than Gal beta 1-3GalNAc. The present study thus reveals the exquisite specificity of A. integrifolia lectin for the T-antigen. Appreciable binding of disaccharides Glc beta 1-3GalNAc and GlcNAc beta 1-3Gal and the very poor binding of beta-linked disaccharides, which instead of Gal and GalNAc contain other sugars at the reducing end, underscore the important contribution made by Gal and GalNAc at the reducing end for recognition by the lectin. The ligand-structure-dependent alterations of the c.d. spectrum in the tertiary structural region of the protein allows the placement of various sugar units in the combining region of the lectin. These studies suggest that the primary subsite (subsite A) can accommodate only Gal or GalNAc or alpha-linked Gal or GalNAc, whereas the secondary subsite (subsite B) can associate either with GalNAc beta Me or Gal beta Me. Considering these factors a likely arrangement for various disaccharides in the binding site of the lectin is proposed. Its exquisite specificity for the authentic T-antigen, Gal beta 1-3GalNAc alpha Ser, together with its virtual non-binding to A- and B-blood-group antigens, Gal beta 1-3GalNAc beta Me and asialo-GM1 should make A. integrifolia lectin a valuable probe for monitoring the expression of T-antigen on cell surfaces.     

Analysis of the interactions of sulfur‐containing amino acids in membrane proteins

The interactions of Met and Cys with other amino acid side chains have received little attention, in contrast to aromatic–aromatic, aromatic–aliphatic or/and aliphatic–aliphatic interactions. Precisely, these are the only amino acids that contain a sulfur atom, which is highly polarizable and, thus, likely to participate in strong Van der Waals interactions. Analysis of the interactions present in membrane protein crystal structures, together with the characterization of their strength in small‐molecule model systems at the ab‐initio level, predicts that Met–Met interactions are stronger than Met–Cys ≈ Met–Phe ≈ Cys–Phe interactions, stronger than Phe–Phe ≈ Phe–Leu interactions, stronger than the Met–Leu interaction, and stronger than Leu–Leu ≈ Cys–Leu interactions. These results show that sulfur‐containing amino acids form stronger interactions than aromatic or aliphatic amino acids. Thus, these amino acids may provide additional driving forces for maintaining the 3D structure of membrane proteins and may provide functional specificity.     

Molecular docking studies on quinazoline antifolate derivatives as human thymidylate synthase inhibitors

We have performed molecular docking on quinazoline antifolates complexed with human thymidylate synthase to gain insight into the structural preferences of these inhibitors. The study was conducted on a selected set of one hundred six compounds with variation in structure and activity. The structural analyses indicate that the coordinate bond interactions, the hydrogen bond interactions, the van der Waals interactions as well as the hydrophobic interactions between ligand and receptor are responsible simultaneously for the preference of inhibition and potency. In this study, fast flexible docking simulations were performed on quinazoline antifolates derivatives as human thymidylate synthase inhibitors. The results indicated that the quinazoline ring of the inhibitors forms hydrophobic contacts with Leu192, Leu221 and Tyr258 and stacking interaction is conserved in complex with the inhibitor and cofactor.     

<Emphasis Type="Italic">Ab initio</Emphasis> fragment molecular orbital studies of influenza virus hemagglutinin–sialosaccharide complexes toward chemical clarification about the virus host range determination

If we predict the host range of new or mutant influenza virus in advance, we are able to measure against pandemic human influenza immediately after the new virus emerges somewhere. Influenza viral hemagglutinin(HA)–sialoside receptor interaction is a target event for in silico chemical prediction studies about the virus host range determination. We theoretically studied avian and human influenza A virus HA H3 subtype complexed with avian or human type receptor Neu5Acα(2-3 or 2-6)Gal analogues by ab initio fragment molecular orbital (FMO) method at the second order Møller–Plesset (MP2)/6–31G level, which can evaluate correctly not only electrostatic interactions but also lipophilic interactions based on van der Waals dispersion force. Avian H3 bound to avian α2-3 11.4 kcal/mol stronger than to human α2-6 in the model complexes with taking account of intermolecular lipophilic interaction. A substitution at the position 226 between Gln(avian) and Leu(human) on influenza H3 HA1 has altered its virus host range between avian and human. In the ab initio FMO studies, binding energy of avian Gln226Leu H3–human α2-6 was quite similar to that in the human H3–human α2-6 complex with amino acid sequence differences at nine positions in the models. This similarity indicates that avian Gln226Leu H3 virus can infect human with the same level as human H3 virus. Opposite mutation Leu226Gln in the human H3 gave the moderate binding energies to avian α2-3 with similarity to avian H3–α2-3 complex that supported our previous virus-sialoside binding assay. Ab initio FMO studies have revealed the relationship between influenza H3 virus host range and H3–α(2-3 or 2-6) receptors binding. Our theoretical approach may predict the infectious level of new viruses and point out some unknown dangerous mutation positions on HA in advance.     

Monoclonal antibody CC3C195, which detects cancer-associated antigens in serum, binds to the human Lea blood group antigen and to its sialylated derivative

Mouse monoclonal antibody CC3C195, which detects elevated levels of its antigen in sera from many patients with colon and pancreatic cancer, binds with high affinity to the sialylated human Lea blood group antigen NeuAc alpha 2-3Gal beta 1-3 [Fuc alpha 1-4]GlcNac . . . and with lower affinity to the Lea blood group antigen itself.     

Structural basis for the selective inhibition of JNK1 by the scaffolding protein JIP1 and SP600125

The c-jun N-terminal kinase (JNK) signaling pathway is regulated by JNK-interacting protein-1 (JIP1), which is a scaffolding protein assembling the components of the JNK cascade. Overexpression of JIP1 deactivates the JNK pathway selectively by cytoplasmic retention of JNK and thereby inhibits gene expression mediated by JNK, which occurs in the nucleus. Here, we report the crystal structure of human JNK1 complexed with pepJIP1, the peptide fragment of JIP1, revealing its selectivity for JNK1 over other MAPKs and the allosteric inhibition mechanism. The van der Waals contacts by the three residues (Pro157, Leu160, and Leu162) of pepJIP1 and the hydrogen bonding between Glu329 of JNK1 and Arg156 of pepJIP1 are critical for the selective binding. Binding of the peptide also induces a hinge motion between the N- and C-terminal domains of JNK1 and distorts the ATP-binding cleft, reducing the affinity of the kinase for ATP. In addition, we also determined the ternary complex structure of pepJIP1-bound JNK1 complexed with SP600125, an ATP-competitive inhibitor of JNK, providing the basis for the JNK specificity of the compound.     

CH/π interactions in the crystal structure of class I MHC antigens and their complexes with peptides

The crystal structure of class I major histocompatibility complex antigens (MHC) bound to their specific ligand peptides were analyzed, in the context of the CH/π interaction, with use of a computer program CHPI. A number of short CH/C_sp² distances have been shown at the boundary of the heavy chain and β2 microglobulin. These interactions are conserved between species, human versus murine. A number of contacts shorter than the conventional van der Waals distance have been disclosed between CH hydrogens and aromatic side-chain groups in the MHC/peptide complexes. The CH/π interaction has been suggested to contribute to the specificity in the complex formation of class I MHC.     

Direct measurements of forces between phosphatidylcholine and phosphatidylethanolamine bilayers in aqueous electrolyte solutions

We report direct measurements of the full interbilayer force laws (force vs. distance) between bilayers of various phosphatidylcholines and phosphatidylethanolamine in aqueous solutions. Bilayers were first deposited on molecularly smooth (mica) surfaces and the interbilayer forces then measured at a resolution of 1 A. Three types of forces were identified: attractive van der Waals forces, repulsive electrostatic (double-layer) forces, and (at short range) repulsive steric hydration forces. Double-layer forces, which arise from ion binding, were insignificant in monovalent salt solutions, e.g., NaCl up to 1 M, but were already present in solutions containing millimolar levels of CaCl2 and MgCl2, giving rise to forces in excellent agreement with theory. Ca2+ binds more strongly than Mg2+, and both bind less to lecithin bilayers in the fluid state (T greater than Tc). The plane of charge coincides with the location of the negative phosphate groups, while the effective plane of origin of the van der Waals force is 4-5 A farther out. In water, the adhesion energies are in the range 0.10-0.15 erg/cm2 for lecithins and approximately 0.8 erg/cm2 for phosphatidylethanolamine. The adhesion energies vary on addition of salt due to changes in the repulsive double-layer and hydration forces rather than to a change in the attractive van der Waals force. The short-range repulsive forces which balance the van der Waals force at separations of 10-30 A are due to a combination of hydration and steric repulsions, the latter arising from thermal motions of head groups and thickness fluctuations of fluid bilayers (above Tc). It is also concluded that bilayer fusion is not simply related to the interbilayer force law.     

Structural modification of indomethacin toward selective inhibition of COX-2 with a significant increase in van der Waals contributions

We have synthesized a fluorinated analogue of indomethacin bearing a 3,3,3-trifluoroprop-1-enyl group at its 2-position and evaluated its inhibitory activity towards the COX-1 and COX-2 enzymes in vitro. The results revealed that this fluorinated analogue exhibited much greater inhibitory activity and selectivity towards COX-2 than indomethacin. The increased affinity between the fluorinated analogue and COX-2 was attributed to a significant increase in van der Waals contacts (i.e. van der Waals contributions in ΔG were ?13.80?kcal/mol for COX-1 and ?18.46?kcal/mol for COX-2), explaining an effect of the fluorine substituent in enzyme selectivity. This newly synthesized fluorinated analogue therefore represents a potent and selective COX-2 inhibitor.     

Carbonic anhydrase inhibitors: X-ray crystallographic studies for the binding of 5-amino-1,3,4-thiadiazole-2-sulfonamide and 5-(4-amino-3-chloro-5-fluorophenylsulfonamido)-1,3,4-thiadiazole-2-sulfonamide to human isoform II

The X-ray crystal structures of 5-amino-1,3,4-thiadiazole-2-sulfonamide (the acetazolamide precursor) and 5-(4-amino-3-chloro-5-fluorophenylsulfonamido)-1,3,4-thiadiazole-2-sulfonamide in complex with the human isozyme II of carbonic anhydrase (CA, EC 4.2.1.1) are reported. The thiadiazole-sulfonamide moiety of the two compounds binds in the canonic manner to the zinc ion and interacts with Thr199, Glu106, and Thr200. The substituted phenyl tail of the second inhibitor was positioned in the hydrophobic part of the binding pocket, at van der Waals distance from Phe131, Val 135, Val141, Leu198, Pro202, and Leu204. These structures may help in the design of better inhibitors of these widespread zinc-containing enzymes.     

Caenorhabditis elegans proteins captured by immobilized Galβ1-4Fuc disaccharide units: assignment of 3 annexins

Galβ1-4Fuc is a key structural motif in Caenorhabditis elegans glycans and is responsible for interaction with C. elegans galectins. In animals of the clade Protostomia, this unit seems to have important roles in glycan–protein interactions and corresponds to the Galβ1-4GlcNAc unit in vertebrates. Therefore, we prepared an affinity adsorbent having immobilized Galβ1-4Fuc in order to capture carbohydrate-binding proteins of C. elegans, which interact with this disaccharide unit. Adsorbed C. elegans proteins were eluted with ethylenediaminetetraacetic acid (EDTA) and followed by lactose (Galβ1-4Glc), digested with trypsin, and were then subjected to proteomic analysis using LC–MS/MS. Three annexins, namely NEX-1, -2, and -3, were assigned in the EDTA-eluted fraction. Whereas, galectins, namely LEC-1, -2, -4, -6, -9, -10, and DC2.3a, were assigned in the lactose-eluted fraction. The affinity of annexins for Galβ1-4Fuc was further confirmed by adsorption of recombinant NEX-1, -2, and -3 proteins to the Galβ1-4Fuc column in the presence of Ca²⁺. Furthermore, frontal affinity chromatography analysis using an immobilized NEX-1 column showed that NEX-1 has an affinity for Galβ1-4Fuc, but no affinity toward Galβ1-3Fuc and Galβ1-4GlcNAc. We would hypothesize that the recognition of the Galβ1-4Fuc disaccharide unit is involved in some biological processes in C. elegans and other species of the Protostomia clade.     

Studies on the structure of a lithium-treated soybean pectin: characteristics of the fragments and determination of the carbohydrate substituents of galacturonic acid

Two galacturonic-acid-containing polysaccharide fractions (ChSS and P) were isolated from soybean meal and subjected to lithium treatment. The fragments obtained were analyzed by using monosaccharide and methylation analyses, and NMR spectroscopy. Lithium degradation of ChSS, followed by sodium borodeuteride reduction, hydrolysis, sodium borohydride reduction, and acetylation afforded alditol acetates, of which the labeled ones reflected residues linked to GalA. As followed from quantifications of the labeled and non-labeled alditols from each constituent monosaccharide by GLC-EIMS, 6 mol% of Ara, 22 mol% of Fuc, 13 mol% of Gal, 53 mol% of Rha, and 57 mol% of Xyl are glycosidically linked to GalA. Analysis of the lithium-treated polymer revealed that it contains arabinogalactan side chains linked to Rha O-4, which consist of a beta-(1 --> 4)-linked galactan substituted with highly branched arabinan chains. On average, an arabinogalactan chain contains up to 29 Gal and 25 Ara residues. Surface plasmon resonance was used to determine conditions for affinity chromatography. Furthermore, this technique confirmed the presence of terminal alpha-Fuc residues in ChSS. Polysaccharide P turned out to be relatively resistant to lithium degradation.     

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

Influenza A viruses attach to alpha-sialosides on the target cell surface by their hemagglutinins, which strictly recognize the difference in sialic acid-galactose linkage. Why does avian virus H3 subtype bind to avian receptor Neu5Ac(alpha2-3)Gal stronger than to human receptor Neu5Ac(alpha2-6)Gal? Why does avian H3 mutated Gln226 to Leu preferentially bind to human receptor? In this paper, we theoretically answer the questions by molecular mechanics and ab initio fragment molecular orbital (FMO) calculations. The binding energy between avian H3 and avian receptor is 8.2kcal/mol larger than that of the avian H3-human receptor complex estimated at the FMO-HF/STO-3G level, which is a reason that avian H3 binds to avian receptor stronger than to human receptor. Avian Leu226 H3 clashes to Gal unit on the avian receptor to quite decrease its binding affinity. In contrast, Gal unit on the human receptor forms intermolecular hydrophobic interaction with avian Leu226 H3 to afford moderate binding affinity.     

Halogenated Benzenes Bound within a Non-polar Cavity in T4 Lysozyme Provide Examples of I?S and I?Se Halogen-bonding

We showed earlier that the mutation of Leu99 to alanine in bacteriophage T4 lysozyme creates an internal cavity of volume ∼ 150 Å³ that binds benzene and a variety of other ligands. As such, this cavity provides an excellent target to study protein-ligand interaction. Here, we use low-temperature crystallography and related techniques to analyze the binding of halogen-incorporated benzenes typified by C₆F₅X, where X = H, F, Cl, Br or I, and C₆H₅X, where X = H or I was also studied. Because of the increased electron density of fluorine relative to hydrogen, the geometry of binding of the fluoro compounds can often be determined more precisely than their hydrogen-containing analogs. All of the ligands bind in essentially the same plane but the center of the phenyl ring can translate by up to 1.2 Å. In no case does the ligand rotate freely within the cavity. The walls of the cavity consist predominantly of hydrocarbon atoms, and in several cases it appears that van der Waals interactions define the geometry of binding. In comparing the smallest with the largest ligand, the cavity volume increases from 181 Å³ to 245 Å³. This shows that the protein is flexible and adapts to the size and shape of the ligand. There is a remarkably close contact of 3.0 Å between the iodine atom on C₆F₅I and the sulfur or selenium atom of Met or SeMet102. This interaction is 1.0 Å less than the sum of the van der Waals radii and is a clear example of a so-called halogen bond. Notwithstanding this close approach, the increase in binding energy for the halogen bond relative to a van der Waals contact is estimated to be only about 0.5-0.7 kcal/mol.