Insights into the role of the aromatic residue in galactose-binding sites: MP2/6-311G++** study on galactose- and glucose-aromatic residue analogue complexes

The presence of an aromatic residue (Trp, Phe, Tyr) facing the nonpolar face of galactose is a common feature of galactose-specific lectins. The interactions such as those between the C-H groups of galactose and the pi-electron cloud of aromatic residues have been characterized as weak hydrogen bonds between soft acids and soft bases, largely governed by dispersive and charge transfer interactions. An analysis of the binding sites of several galactose-specific lectins revealed that the spatial position-orientation of galactose relative to the binding site aromatic residue varies substantially. The effect of variations in position-orientations of galactose on the interaction energies of galactose-aromatic residue complexes has not been determined so far. In view of this, MP2/6-311G++** calculations were performed on galactose- and glucose-aromatic residue analogue complexes in eight position-orientations. The results show that the strength of the C-H...pi interactions in galactose-aromatic residue complexes is comparable to that of a hydrogen bond. Rather than the type of aromatic residue, the position-orientation of the saccharide appears to be more critical in determining the strength of their interactions. Earlier studies have found the binding site aromatic residue to be critical, but its role was not clear. This study shows that the aromatic residue is important for discriminating galactose from glucose, in addition to its contribution to binding energy.     

Role of CH/pi interactions in substrate binding by Escherichia coli beta-galactosidase

Interactions between carbohydrates and aromatic amino-acid residues are often observed in structures of carbohydrate-protein complexes. They are characterized by an orientation of the pyranose or furanose ring parallel with the aromatic ring of amino-acid residues. An important role in the formation of these complexes is supposed to be played by CH/pi interactions. This paper presents an ab initio quantum chemistry study of CH/pi interactions between beta-galactosidase from E. coli and its substrates and products. The energy stabilizing the interaction between Trp999 residue and substrate bound in the shallow binding mode was calculated at the MP2/6-31+G(d) level as 5.2kcalmol(-1) for the glucose moiety of allolactose, 2.4kcalmol(-1) for the galactose moiety of allolactose and 5.0kcalmol(-1) for the glucose moiety of lactose. The energy stabilizing the interaction between Trp568 residue and galactose in the deep binding mode was calculated as 2.7kcalmol(-1). Interaction energies at the HF/6-31+G(d) and B3LYP/6-31+G(d) levels were small or repulsive; therefore, highly correlated ab initio methods were necessary to study these interactions. These unexpectedly strong interactions give a rationale for allolactose formation and illustrate the role of the Trp999 residue. In addition, this illustrates the importance of CH/pi interactions for the function of carbohydrate-binding proteins and carbohydrate-processing enzymes.     

Escherichia coli Heat-Labile Enterotoxin Binds to Glycosylated Proteins with Lactose by Amino Carbonyl Reaction

The binding of Escherichia coli heat-labile enterotoxin (LT) type I to glycosylated proteins with lactose (Galβ1-4Glc) by amino carbonyl reaction was studied by the Western blot assay and by the microtiter well binding assay. LT bound to a lactose-α-lactalbumin amino carbonyl product (Lac-LA), whereas cholera toxin did not. The binding ability of Lac-LA was abolished by β-galactosidase treatment, indicating that the terminal galactose is essential for the binding of LT. The binding of LT to Lac-LA was inhibited by galactose and lactose, and most effectively inhibited by lactulose (Galβ1-4Fru), which is a structural analog of the Amadori rearrangement product of the amino carbonyl reaction between lactose and an ε-amino group of a lysine residue (lactuloselysine). The results suggest that LT recognizes the portion of lactuloselysine in Lac-LA. LT also bound to a melibiose (Galα1-6Glc)-α-lactalbumin amino carbonyl product (Mel-LA), but the binding ability of Mel-LA was weaker than that of Lac-LA, suggesting that the β1-4 linked terminal galactose is dispensable but preferable for the binding. Furthermore, LT bound to the amino carbonyl products of lactose with β-lactoglobulin, caseins, bovine serum albumin, and ovalbumin. These results indicate that LT binds to the amino carbonyl products between proteins and sugars containing the terminal galactose, such as lactose.     

Separation and characteristics of galactose-1-phosphate and glucose-1-phosphate uridyltransferase from fruit peduncles of cucumber

Galactose-1-phosphate uridyltransferase (EC 2.7.7.10), responsible for the conversion of galactose-1-phosphate (Gal-1-P) to uridine diphosphate galactose (UDPgal) was examined in fruit peduncles of Cucumis sativus L. Two uridyltransferases (pyrophosphorylases), from I and II, were partially purified and resolved on a diethylamino-ethyl-cellulose column. Form I can utilize glucose-1-phosphate (Glc-1-P), while form II can utilize either Gal-1-P or Glc-1-P, with a preference for Gal-1-P. Form I was more heat stable than form II. Both Glc-1-P and Gal-1-P activities of form II were inactivated at the same rate by heating. The finding of a uridyltransferase with preference for Gal-1-P indicates that cucumber may have a Gal-1-P uridyltransferase (pyrophosphorylase) pathway for the catabolism of stachyose in the peduncles. The absence of the enzyme UDP-glucose-hexose-1-phosphate uridyltransferase (EC 2.7.7.12) in this tissue rules out catabolism by the classical Leloir pathway. The incorporation of carbon from UDPglc into Glc-1-P as opposed to sucrose may be regulated by the activities of the uridyltransferases. Pyrophosphate, in the same concentration range, inhibits UDP-gal formation (K_i=0.58±0.10 mM) and stimulates Glc-1-P formation. The ratio of units of pyrophosphatase to units of Gal-1-P uridyltransferase was higher in peduncles from growing fruit than from unpollinated fruit. Modulation of carbon partitioning through a uridyltransferase pathway may be a factor controlling growth of the cucumber fruit.Abbreviations Gal-1-P Galactose-1-phosphate - Glc-1-P glucose-1-phosphate - UDPgal uridine diphosphate galactose - UDPglc uridine diphosphate glucose Paper No. 6908 of the Journal Series of the North Carolina Agricultural Research Service, Raleigh. The use of trade names in this publication does not imply endorsement by the North Carolina Agricultural Research Service of products named, nor criticism of similar ones not mentioned     

Architecture of the sugar binding sites in carbohydrate binding proteins-a computer modeling study

Different sugars, Gal, GalNAc and Man were docked at the monosaccharide binding sites of Erythrina corallodenron (EcorL), peanut lectin (PNA), Lathyrus ochrus (LOLI), and pea lectin (PSL). To study the lectin-carbohydrate interactions, in the complexes, the hydroxymethyl group in Man and Gal favors, gg and gt conformations respectively, and is the dominant recognition determination. The monosaccharide binding site in lectins that are specific to Gal/GalNAc is wider due to the additional amino acid residues in loop D as compared to that in lectins specific to Man/Glc, and affects the hydrogen bonds of the sugar involving residues from loop D, but not its orientation in the binding site. The invariant amino acid residues Asp from loop A, and Asn and an aromatic residue (Phe or Tyr) in loop C provides the basic architecture to recognize the common features in C4 epimers. The invariant Gly in loop B together with one or two residues in the variable region of loop D/A holds the sugar tightly at both ends. Loss of any one of these hydrogen bonds leads to weak interaction. While the subtle variations in the sequence and conformation of peptide fragment that resulted due to the size and location of gaps present in amino acid sequence in the neighborhood of the sugar binding site of loop D/A seems to discriminate the binding of sugars which differ at C4 atom (galacto and gluco configurations). The variations at loop B are important in discriminating Gal and GalNAc binding. The present study thus provides a structural basis for the observed specificities of legume lectins which uses the same four invariant residues for binding. These studies also bring out the information that is important for the design/engineering of proteins with the desired carbohydrate specificity.     

Cloning and expression analysis of a UDP-galactose/glucose pyrophosphorylase from melon fruit provides evidence for the major metabolic pathway of galactose metabolism in raffinose oligosaccharide metabolizing plants

The Cucurbitaceae translocate a significant portion of their photosynthate as raffinose and stachyose, which are galactosyl derivatives of sucrose. These are initially hydrolyzed by alpha-galactosidase to yield free galactose (Gal) and, accordingly, Gal metabolism is an important pathway in Cucurbitaceae sink tissue. We report here on a novel plant-specific enzyme responsible for the nucleotide activation of phosphorylated Gal and the subsequent entry of Gal into sink metabolism. The enzyme was antibody purified, sequenced, and the gene cloned and functionally expressed in Escherichia coli. The heterologous protein showed the characteristics of a dual substrate UDP-hexose pyrophosphorylase (PPase) with activity toward both Gal-1-P and glucose (Glc)-1-P in the uridinylation direction and their respective UDP-sugars in the reverse direction. The two other enzymes involved in Glc-P and Gal-P uridinylation are UDP-Glc PPase and uridyltransferase, and these were also cloned, heterologously expressed, and characterized. The gene expression and enzyme activities of all three enzymes in melon (Cucumis melo) fruit were measured. The UDP-Glc PPase was expressed in melon fruit to a similar extent as the novel enzyme, but the expressed protein was specific for Glc-1-P in the UDP-Glc synthesis direction and did not catalyze the nucleotide activation of Gal-1-P. The uridyltransferase gene was only weakly expressed in melon fruit, and activity was not observed in crude extracts. The results indicate that this novel enzyme carries out both the synthesis of UDP-Gal from Gal-1-P as well as the subsequent synthesis of Glc-1-P from the epimerase product, UDP-Glc, and thus plays a key role in melon fruit sink metabolism.     

Structural basis of carbohydrate recognition by lectin II from Ulex europaeus, a protein with a promiscuous carbohydrate-binding site

Protein-carbohydrate interactions are the language of choice for inter- cellular communication. The legume lectins form a large family of homologous proteins that exhibit a wide variety of carbohydrate specificities. The legume lectin family is therefore highly suitable as a model system to study the structural principles of protein-carbohydrate recognition. Until now, structural data are only available for two specificity families: Man/Glc and Gal/GalNAc. No structural data are available for any of the fucose or chitobiose specific lectins.The crystal structure of Ulex europaeus (UEA-II) is the first of a legume lectin belonging to the chitobiose specificity group. The complexes with N-acetylglucosamine, galactose and fucosylgalactose show a promiscuous primary binding site capable of accommodating both N-acetylglucos amine or galactose in the primary binding site. The hydrogen bonding network in these complexes can be considered suboptimal, in agreement with the low affinities of these sugars. In the complexes with chitobiose, lactose and fucosyllactose this suboptimal hydrogen bonding network is compensated by extensive hydrophobic interactions in a Glc/GlcNAc binding subsite. UEA-II thus forms the first example of a legume lectin with a promiscuous binding site and illustrates the importance of hydrophobic interactions in protein-carbohydrate complexes. Together with other known legume lectin crystal structures, it shows how different specificities can be grafted upon a conserved structural framework.     

Structural basis for langerin recognition of diverse pathogen and mammalian glycans through a single binding site

Langerin mediates the carbohydrate-dependent uptake of pathogens by Langerhans cells in the first step of antigen presentation to the adaptive immune system. Langerin binds to an unusually diverse number of endogenous and pathogenic cell surface carbohydrates, including mannose-containing O-specific polysaccharides derived from bacterial lipopolysaccharides identified here by probing a microarray of bacterial polysaccharides. Crystal structures of the carbohydrate-recognition domain from human langerin bound to a series of oligomannose compounds, the blood group B antigen, and a fragment of β-glucan reveal binding to mannose, fucose, and glucose residues by Ca²⁺ coordination of vicinal hydroxyl groups with similar stereochemistry. Oligomannose compounds bind through a single mannose residue, with no other mannose residues contacting the protein directly. There is no evidence for a second Ca²⁺-independent binding site. Likewise, a β-glucan fragment, Glcβ1-3Glcβ1-3Glc, binds to langerin through the interaction of a single glucose residue with the Ca²⁺ site. The fucose moiety of the blood group B trisaccharide Galα1-3(Fucα1-2)Gal also binds to the Ca²⁺ site, and selective binding to this glycan compared to other fucose-containing oligosaccharides results from additional favorable interactions of the nonreducing terminal galactose, as well as of the fucose residue. Surprisingly, the equatorial 3-OH group and the axial 4-OH group of the galactose residue in 6SO₄-Galβ1-4GlcNAc also coordinate Ca²⁺, a heretofore unobserved mode of galactose binding in a C-type carbohydrate-recognition domain bearing the Glu-Pro-Asn signature motif characteristic of mannose binding sites. Salt bridges between the sulfate group and two lysine residues appear to compensate for the nonoptimal binding of galactose at this site.     

Computational approach for the design of AP1867 analogs: aiming at new synthetic routes for potential immunosuppressant agents

Molecular modelling and synthetic arguments are valuable tools for the design of potential immunosuppressant agents. In this paper, eight proline-based compounds related to the AP1867 structure are studied and at least one of them is found to be a structurally good candidate for the inhibition of FKBP protein. Theoretical calculations were carried out to locate the most energetically favorable chemical substituent group relative to a core skeleton group on interaction with the FKBP binding cavity. Connolly accessible surface calculations have complemented the molecular mechanics and dynamics approaches. Calculated results were also analyzed on the basis of hydrogen bond interactions, relative energies of interaction, root-mean square deviations of amino acid residues of the crystallized protein, and orientation of the substituent groups within the active site. The results show a significant reduction in the relative interaction energies and very good shape complementarities between our final analog compound and the FKBP binding pocket.     

Density functional theory calculations on entire proteins for free energies of binding: Application to a model polar binding site

In drug optimization calculations, the molecular mechanics Poisson‐Boltzmann surface area (MM‐PBSA) method can be used to compute free energies of binding of ligands to proteins. The method involves the evaluation of the energy of configurations in an implicit solvent model. One source of errors is the force field used, which can potentially lead to large errors due to the restrictions in accuracy imposed by its empirical nature. To assess the effect of the force field on the calculation of binding energies, in this article we use large‐scale density functional theory (DFT) calculations as an alternative method to evaluate the energies of the configurations in a “QM‐PBSA” approach. Our DFT calculations are performed with a near‐complete basis set and a minimal parameter implicit solvent model, within the self‐consistent calculation, using the ONETEP program on protein–ligand complexes containing more than 2600 atoms. We apply this approach to the T4‐lysozyme double mutant L99A/M102Q protein, which is a well‐studied model of a polar binding site, using a set of eight small aromatic ligands. We observe that there is very good correlation between the MM and QM binding energies in vacuum but less so in the solvent. The relative binding free energies from DFT are more accurate than the ones from the MM calculations, and give markedly better agreement with experiment for six of the eight ligands. Furthermore, in contrast to MM‐PBSA, QM‐PBSA is able to correctly predict a nonbinder. Proteins 2014; 82:3335–3346. © 2014 Wiley Periodicals, Inc.     

Evidence for involvement of tryptophan residue in the low-affinity saccharide binding site of ricin D

The nature of the saccharide-binding site of ricin D, which is a galactose- and N-acetylgalactosamine-specific lectin, was studied by chemical modification and spectroscopy. With excitation at 290 nm, ricin D displayed a fluorescence spectrum with a maximum at 335 nm. Upon binding of the specific saccharides, the spectrum shifted to shorter wavelength by 3 nm. However, binding of galactosamine and N-acetylgalactosamine failed to induce such a change in the fluorescence spectrum. The interaction of ricin D with its specific saccharides was analyzed in terms of the variation of the intensity at 320 nm as a function of saccharide concentration. The results indicate that the change in the fluorescence spectrum induced by saccharide binding is attributable to the binding of saccharide to the low-affinity (LA-) binding site of ricin D. The cytoagglutinating activity of ricin D decreased to 2% upon modification of two tryptophan residues/mol with N-bromosuccinimide at pH 4.0, but in the presence of galactose or lactose one tryptophan residue/mol remained unmodified, and a fairly high cytoagglutinating activity was retained. Galactosamine and N-acetylgalactosamine did not show such a protective effect. Spectroscopic analyses indicate that the decrease in the cytoagglutinating activity of ricin D upon tryptophan modification is principally due to the loss of the saccharide binding activity of the LA-binding site. The results suggest that one tryptophan residue is essential for saccharide binding at the LA-binding site, which can bind galactose and lactose but lacks the ability to bind N-acetylgalactosamine and galactosamine.     

Do electrostatic interactions destabilize protein–nucleic acid binding?

The negatively charged phosphates of nucleic acids are often paired with positively charged residues upon binding proteins. It was thus counter-intuitive when previous Poisson-Boltzmann (PB) calculations gave positive energies from electrostatic interactions, meaning that they destabilize protein-nucleic acid binding. Our own PB calculations on protein-protein binding have shown that the sign and the magnitude of the electrostatic component are sensitive to the specification of the dielectric boundary in PB calculations. A popular choice for the boundary between the solute low dielectric and the solvent high dielectric is the molecular surface; an alternative is the van der Waals (vdW) surface. In line with results for protein-protein binding, in this article, we found that PB calculations with the molecular surface gave positive electrostatic interaction energies for two protein-RNA complexes, but the signs are reversed when the vdW surface was used. Therefore, whether destabilizing or stabilizing effects are predicted depends on the choice of the dielectric boundary. The two calculation protocols, however, yielded similar salt effects on the binding affinity. Effects of charge mutations differentiated the two calculation protocols; PB calculations with the vdW surface had smaller deviations overall from experimental data.     

Effect of water coordination on competition between π and non-π cation binding sites in aromatic amino acids: l-phenylalanine, l-tyrosine, and l-tryptophan Li+, Na+, and K+ complexes

Quantum chemistry methods have been applied to charged complexes of the alkali metals Li(+), Na(+), and K(+) with the aromatic amino acids (AAAs) phenylalanine (Phe), tyrosine (Tyr), and tryptophan (Trp). The geometries of 72 different complexes (Phe·M, Tyr·M, Trp·M, M is Li(+), Na(+), or K(+)) were completely optimized at the B3LYP/6-311+G(d,p) level of density functional theory. The solvent effect on the geometry and stability of individual complexes was studied by making use of a microsolvation model. The interaction enthalpies, entropies, and Gibbs energies of nine different complexes of the systems Phe·M, Tyr·M, and Trp·M (M is Li(+), Na(+), or K(+)) were also determined at the B3LYP density functional level of theory. The calculated Gibbs binding energies of the M(+)-AAA complexes follow the order Phe < Tyr < Trp for all three metal cations studied. Among the three AAAs studied, the indole ring of Trp is the best π donor for alkali metal cations. Our calculations demonstrated the existence of strong cation-π interactions between the alkali metals and the aromatic side chains of the three AAAs. These AAAs comprise about 8% of all known protein sequences. Thus, besides the potential for hydrogen-bond interaction, aromatic residues of Phe, Tyr, and Trp show great potential for π-donor interactions. The existence of cation-π interaction in proteins has also been demonstrated experimentally. However, more complex experimental studies of metal cation-π interaction in diverse biological systems will no doubt lead to more exact validation of these investigations.     

Investigation of the nature of the methionine-pi interaction in beta-hairpin peptide model systems

There are frequent contacts between aromatic rings and sulfur atoms in proteins. However, it is unclear to what degree this putative interaction is stabilizing and what the nature of the interaction is. We have investigated the aryl-sulfur interaction by placing a methionine residue diagonal to an aromatic ring on the same face of a beta-hairpin, which places the methionine side chain in close proximity to the aryl side chain. The methionine (Met)-aryl interaction was compared with an equivalent hydrophobic and cation-pi interaction in the context of the beta-hairpin. The interaction between phenylalanine (Phe), tryptophan (Trp), or cyclohexylalanine (Cha) and Met stabilized the beta-hairpin by -0.3 to -0.5 kcal mole(-1), as determined by double-mutant cycles. The peptides were subjected to thermal denaturations that suggest a hydrophobic driving force for the interactions between Met and Trp or Cha. The observed interaction of Met or norleucine (Nle) with Trp or Cha are quite similar, implying a hydrophobic driving force for the Met-pi interaction. However, the thermodynamic data suggest that there may be some differences between the interaction of Met with Trp and Phe and that there may be a small thermodynamic component to the Met...Phe interaction.     

Protein–ligand interactions: Probing the energetics of a putative cation–π interaction

In order to probe the energetics associated with a putative cation–π interaction, thermodynamic parameters are determined for complex formation between the Grb2 SH2 domain and tripeptide derivatives of RCO–pTyr–Ac6c–Asn wherein the R group is varied to include different alkyl, cycloalkyl, and aryl groups. Although an indole ring is reputed to have the strongest interaction with a guanidinium ion, binding free energies, ΔG°, for derivatives of RCO–pTyr–Ac6c–Asn bearing cyclohexyl and phenyl groups were slightly more favorable than their indolyl analog. Crystallographic analysis of two complexes reveals that test ligands bind in similar poses with the notable exception of the relative orientation and proximity of the phenyl and indolyl rings relative to an arginine residue of the domain. These spatial orientations are consistent with those observed in other cation–π interactions, but there is no net energetic benefit to such an interaction in this biological system. Accordingly, although cation–π interactions are well documented as important noncovalent forces in molecular recognition, the energetics of such interactions may be mitigated by other nonbonded interactions and solvation effects in protein–ligand associations.     

Effect of ligand-induced conformational changes on the reactivity of specific sulfhydryl residues in rat brain hexokinase

Rat brain hexokinase (ATP:D-hexose 6-phosphotransferase, EC 2.7.1.1) contains 21 cysteine residues. On the basis of the amino acid sequence of the enzyme, these are predicted to be distributed among 14 peptides produced by tryptic digestion. Ten of these peptides, containing cysteine residues derivatized by reaction with the specific sulfhydryl reagent 2-bromoacetamido-4-nitrophenol have been identified in HPLC peptide maps; the four missing peptides are predicted to be relatively large and hydrophobic in character, properties that may have prevented their detection under the present conditions. The sequences encompassed by the 10 identified peptides include 12 of the 21 cysteine residues in the enzyme. The relative reactivity of these sulfhydryl groups with 2-bromoacetamido-4-nitrophenol has been assessed, and is in general accord with what might be predicted on the basis of their accessibility in the previously proposed structure for this enzyme. The effect of various ligands on reactivity of identified sulfhydryl groups has been determined; unique patterns of altered reactivity, resulting from ligand-induced conformational changes, have been observed. Biphasic effects were observed with increasing concentrations of either glucose 6-phosphate (Glc-6-P) or Pi. In both cases, decreased reactivity of sulfhydryls in the N-terminal half of the molecule was observed at low concentrations of the ligand, while further increase in ligand concentration resulted in decreased reactivity of sulfhydryl groups in the C-terminal half. In contrast, sulfhydryls in both N- and C-terminal halves were protected concomitantly by increasing concentrations of Glc. These results are consistent with previous studies that indicated (a) the existence of two sites for binding of Glc-6-P or Pi, a high affinity site in the N-terminal half and a site with lower affinity in the C-terminal half of the brain hexokinase molecule, and (b) binding of Glc to a single site located in the C-terminal half but evoking conformational effects throughout the molecule; the glucose analog, N-acetylglucosamine, previously shown to have more limited effects on conformation, affected reactivity of sulfhydryl groups only in the C-terminal half of the molecule. As reflected by effects on reactivity of sulfhydryl groups, conformational changes induced by binding of nucleotides depends markedly on the specific nature of the purine or pyrimidine base as well as the length and chelation status of the polyphosphate side chain. These results focus attention on specific regions of the molecule (the immediate environment of the sulfhydryl groups) that are affected by the binding of these ligands.     

Studies on pituitary follitropin. III. Functional role of the amino groups in subunit and receptor interactions of the ovine hormone.

The effectiveness of Glc, mannose, 2-deoxyglucose, fructose, galactose, arabinose, and N-acetylglucosamine at protecting rat brain hexokinase (ATP: d-hexose 6-phosphotransferase, EC 2.7.1.1) from inactivation by chymotrypsin, glutaraldehyde, heat, and Ellman's reagent have been compared. The relative effectiveness at protecting against these inactivating agents decreases in the order Glc > mannose > 2-deoxyglucose > fructose, galactose, and arabinose, the last three providing no significant protection at all. The nonphosphorylatable substrate analog, N-acetylglucosamine, provides substantial protection against heat inactivation, but is ineffective against inactivation by the other agents. Similar inactivation studies were conducted using several hexose 6-phosphates. Glc-6-P and 1,5-anhydroglucitol-6-P provided substantial protection while 2-deoxyglucose-6-P, fructose-6-P, mannose-6-P, and galactose-6-P were all relatively ineffective at protecting hexokinase activity. The protective effect of these ligands is taken as an indication of ligand-induced conformational changes in brain hexokinase. The results are interpreted in terms of, and considered to support, a recently proposed model (J. E. Wilson, 1978, Arch. Biochem. Biophys.185, 88–99) in which the suitability of a carbohydrate as a substrate depends directly on its ability to induce specific conformational changes prerequisite for subsequent catalytic events while the inhibitory effectiveness of a hexose 6-phosphate is likewise related to its ability to evoke appropriate conformational change in the enzyme. Synergistic interactions between hexose and hexose-6-P binding sites, first reported for Glc and Glc-6-P by Ellison et al. (1975, J. Biol. Chem.250, 1864–1871), have been confirmed. Thus, although fructose and galactose were themselves quite ineffective at providing protection against inactivation of hexokinase by chymotrypsin or glutaraldehyde, they greatly increased the protection afforded by suboptimal (with respect to degree of protection in the absence of added hexose) levels of Glc-6-P. Conversely, the 6-phosphates of fructose, galactose, mannose, and 2-deoxyglucose, which were themselves ineffective at protecting the enzyme activity, markedly enhanced the protection provided by suboptimal levels of Glc or mannose. Based on the relationship between this enhancement of Glc-dependent protection and the hexose-6-P concentration, the dissociation constants for the complexes of hexokinase with 2-deoxyglucose-6-P and mannose-6-P were estimated to be ?0.5 mm.     

Sugar-binding sites of the HA1 subcomponent of Clostridium botulinum type C progenitor toxin

Clostridium botulinum type C 16S progenitor toxin contains a hemagglutinin (HA) subcomponent, designated HA1, which appears to play an important role in the effective internalization of the toxin in gastrointestinal epithelial cells and in creating a broad specificity for the oligosaccharide structure that corresponds to various targets. In this study, using the recombinant protein fused to glutathione S-transferase, we investigated the binding specificity of the HA1 subcomponent to sugars and estimated the binding sites of HA1 based on X-ray crystallography and soaking experiments using various sugars. N-Acetylneuraminic acid, N-acetylgalactosamine, and galactose effectively inhibited the binding that occurs between glutathione S-transferase-HA1 and mucins, whereas N-acetylglucosamine and glucose did not inhibit it. The crystal structures of HA1 complex with N-acetylneuraminic acid, N-acetylgalactosamine, and galactose were also determined. There are two sugar-binding sites, sites I and II. Site I corresponds to the electron densities noted for all sugars and is located at the C-terminal β-trefoil domain, while site II corresponds to the electron densities noted only for galactose. An aromatic amino acid residue, Trp176, at site I has a stacking interaction with the hexose ring of the sugars. On the other hand, there is no aromatic residue at site II; thus, the interaction with galactose seems to be poor. The double mutant W176A at site I and D271F at site II has no avidity for N-acetylneuraminic acid but has avidity for galactose. In this report, the binding specificity of botulinum C16S toxin HA1 to various sugars is demonstrated based on its structural features.     

Binding free energies and free energy components from molecular dynamics and Poisson-Boltzmann calculations. Application to amino acid recognition by aspartyl-tRNA synthetase

Specific amino acid binding by aminoacyl-tRNA synthetases (aaRS) is necessary for correct translation of the genetic code. Engineering a modified specificity into aminoacyl-tRNA synthetases has been proposed as a means to incorporate artificial amino acid residues into proteins in vivo. In a previous paper, the binding to aspartyl-tRNA synthetase of the substrate Asp and the analogue Asn were compared by molecular dynamics free energy simulations. Molecular dynamics combined with Poisson-Boltzmann free energy calculations represent a less expensive approach, suitable for examining multiple active site mutations in an engineering effort. Here, Poisson-Boltzmann free energy calculations for aspartyl-tRNA synthetase are first validated by their ability to reproduce selected molecular dynamics binding free energy differences, then used to examine the possibility of Asn binding to native and mutant aspartyl-tRNA synthetase. A component analysis of the Poisson-Boltzmann free energies is employed to identify specific interactions that determine the binding affinities. The combined use of molecular dynamics free energy simulations to study one binding process thoroughly, followed by molecular dynamics and Poisson-Boltzmann free energy calculations to study a series of related ligands or mutations is proposed as a paradigm for protein or ligand design.The binding of Asn in an alternate, "head-to-tail" orientation observed in the homologous asparagine synthetase is analyzed, and found to be more stable than the "Asp-like" orientation studied earlier. The new orientation is probably unsuitable for catalysis. A conserved active site lysine (Lys198 in Escherichia coli) that recognizes the Asp side-chain is changed to a leucine residue, found at the corresponding position in asparaginyl-tRNA synthetase. It is interesting that the binding of Asp is calculated to increase slightly (rather than to decrease), while that of Asn is calculated, as expected, to increase strongly, to the same level as Asp binding. Insight into the origin of these changes is provided by the component analyses. The double mutation (K198L,D233E) has a similar effect, while the triple mutation (K198L,Q199E,D233E) reduces Asp binding strongly. No binding measurements are available, but the three mutants are known to have no ability to adenylate Asn, despite the "Asp-like" binding affinities calculated here. In molecular dynamics simulations of all three mutants, the Asn ligand backbone shifts by 1-2 A compared to the experimental Asp:AspRS complex, and significant side-chain rearrangements occur around the pocket. These could reduce the ATP binding constant and/or the adenylation reaction rate, explaining the lack of catalytic activity in these complexes. Finally, Asn binding to AspRS with neutral K198 or charged H449 is considered, and shown to be less favorable than with the charged K198 and neutral H449 used in the analysis.     

Ricin D-saccharide interaction as studied by ultraviolet difference spectroscopy

The interaction of ricin D with specific saccharides was investigated by ultraviolet difference spectroscopy. Upon binding to saccharides, ricin D displayed ultraviolet difference spectra with maxima at 280 nm and 288 nm. Such difference spectra suggest that the environment of a tyrosine residue(s) located at or near the saccharide-binding site is changed by the binding of saccharide. In addition to the two positive peaks, a small trough was observed around 300 nm in the complexes with galactose-containing saccharides but not in the complex with N-acetylgalactosamine or galactosamine, suggesting the participation of tryptophan in the binding with galactose-containing saccharides. The magnitude of the difference maxima increased with increasing concentration of saccharides until the binding site was saturated. From the variation of the maximum at 288 nm as a function of saccharide concentration, the association constants were obtained for the binding of saccharides to ricin D at various temperatures and pH's. The saccharide binding of ricin D decreased with increasing temperature and with decreasing pH below pH 6.0. It was suggested that difference maximum at 288 nm observed in the ricin D-saccharide interaction reflects the binding of saccharides to the high-affinity saccharide-binding site of ricin D.