Aromatic-carbohydrate interactions: an NMR and computational study of model systems

NMR investigations of glycan conformation, dynamics, and interactions

Glycans are ubiquitous in nature, decorating our cells and serving as the initial points of contact with any visiting entities. These glycan interactions are fundamental to host-pathogen recognition and are related to various diseases, including inflammation and cancer. Therefore, understanding the conformations and dynamics of glycans, as well as the key features that regulate their interactions with proteins, is crucial for designing new therapeutics. Due to the intrinsic flexibility of glycans, NMR is an essential tool for unravelling these properties. In this review, we describe the key NMR parameters that can be extracted from the different experiments, and which allow us to deduce the necessary geometry and molecular motion information, with a special emphasis on assessing the internal motions of the glycosidic linkages. We specifically address the NMR peculiarities of various natural glycans, from histo-blood group antigens to glycosaminoglycans, and also consider the special characteristics of their synthetic analogues (glycomimetics). Finally, we discuss the application of NMR protocols to study glycan-related molecular recognition events, both from the carbohydrate and receptor perspectives, including the use of stable isotopes and paramagnetic NMR methods to overcome the inherent degeneracy of glycan chemical shifts.

CH-π Interactions Are Required for Human Galectin-3 Function

Glycan-binding proteins, or lectins, recognize distinct structural elements of polysaccharides, to mediate myriad biological functions. Targeting glycan-binding proteins involved in human disease has been challenging due to an incomplete understanding of the molecular mechanisms that govern protein-glycan interactions. Bioinformatics and structural studies of glycan-binding proteins indicate that aromatic residues with the potential for CH-π interactions are prevalent in glycan-binding sites. However, the contributions of these CH-π interactions to glycan binding and their relevance in downstream function remain unclear. An emblematic lectin, human galectin-3, recognizes lactose and N-acetyllactosamine-containing glycans by positioning the electropositive face of a galactose residue over the tryptophan 181 (W181) indole forming a CH-π interaction. We generated a suite of galectin-3 W181 variants to assess the importance of these CH-π interactions to glycan binding and function. As determined experimentally and further validated with computational modeling, variants with smaller or less electron-rich aromatic side chains (W181Y, W181F, W181H) or sterically similar but nonaromatic residues (W181M, W181R) showed poor or undetectable binding to lactose and attenuated ability to bind mucins or agglutinate red blood cells. The latter functions depend on multivalent binding, highlighting that weakened CH-π interactions cannot be overcome by avidity. Two galectin-3 variants with disrupted hydrogen bonding interactions (H158A and E184A) showed similarly impaired lactose binding. Molecular simulations demonstrate that all variants have decreased binding orientation stability relative to native galectin-3. Thus, W181 collaborates with the endogenous hydrogen bonding network to enhance binding affinity for lactose, and abrogation of these CH-π interactions is as deleterious as eliminating key hydrogen bonding interactions. These findings underscore the critical roles of CH-π interactions in carbohydrate binding and lectin function and will aid the development of novel lectin inhibitors.

Molecular recognition and proteoglycan mimic arrangement: modulating cisplatin toxicity

We have demonstrated that cisplatin (CP), an anticancer drug, showed a preference for binding the sulfated-L-iduronic acid (S-L-IdoA) unit over the sulfated-D-glucuronic acid unit of heparan sulfate. The multivalency of S-L-IdoA, such as in the proteoglycan mimic, resulted in distinct modes of cell-surface engineering in normal and cancer cells, with these disparities having a significant impact on CP-mediated toxicity.

Antiradical Activity of Lignans from Cleistanthus sumatranus: Theoretical Insights into the Mechanism, Kinetics, and Solvent Effects

Sumatranus lignans (SL) isolated from Cleistanthus sumatranus have demonstrated bioactivities, e.g., they were shown to exhibit immunosuppressive properties in previous research. Their structure suggests potential antioxidant activity that has not attracted any attention thus far. Consistently, a comprehensive analysis of the antioxidant activity of these compounds is highly desirable with the view of prospective medical applications. In this work, the mechanism and kinetics of the antiradical properties of SL against hydroperoxyl radicals were studied by using calculations based on density functional theory (DFT). In the lipid medium, it was discovered that SL reacted with HOO• through the formal hydrogen transfer mechanism with a rate constant of 101-105 M-1 s-1, whereas in aqueous media, the activity primarily occurred through the sequential proton loss electron transfer mechanism with rate constants of 102-108 M-1 s-1. In both lipidic and aqueous environments, the antiradical activity of compounds 6 and 7 exceeds that of resveratrol, ascorbic acid, and Trolox. These substances are therefore predicted to be good radical scavengers in physiological environments.

Effect of Birch Sap as Solvent and Source of Bioactive Compounds in Casein and Gelatine Films

Birch sap consists of a natural water-based solution with valuable compounds such as minerals, sugars, organic acids and phenolic compounds that can be used advantageously in the preparation of edible films. In this study, gelatine- and casein-based films were prepared using birch sap as biopolymer solvent and source of bioactive compounds with the aim of developing new bioactive materials for food packaging. The physical, mechanical, barrier, antioxidant and iron-chelating properties of the obtained films were investigated. Birch sap enhanced the mechanical properties of the films by increasing puncture strength and flexibility, as well as their ultraviolet-visible light barrier properties. In addition, the presence of bioactive compounds endowed the birch sap films with an antioxidant capacity of almost 90% and an iron-chelating capacity of 40-50% with respect to the control films. Finally, to test these films as food packaging material, a photosensitive curcumin solution was packed and exposed to ultraviolet light. Tested films were able to protect curcumin against photodegradation, and the presence of bioactive compounds inside the birch-sap-enriched materials offered an additional 10% photoprotective effect compared to control films. Results showed the potential of birch sap as an environmentally friendly biopolymer solvent and plasticizer that can improve the mechanical and photoprotective properties of the prepared materials.

How solvent determines the molecular reactive conformation and the selectivity: Solvation spheres and energy

In solution, the solvent determines the molecular conformation and the chemical reaction viability and selectivity. When solvent-solute and solvent-solvent interactions present similar strengths, explicit salvation is the best way to describe a system. The problem to solve is how big the explicit shell should be. In this paper, we want to answer one of the fundamental questions in the implementation of explicit solvation, exactly how many solvent molecules should be added and where they should be placed. Here we determine the first solvent sphere around a molecule and describe how it controls the conformation and selectivity of a selected reaction. NMR experiments were carried out to identify the number of solvent molecules around the solute that constitutes the first solvent sphere, and the interaction between this solvent sphere and the solute was detected using DFT and QTAIM calculations. A new approach to the solvation energy is presented. Finally, we established the role of solvent molecules in the conformation of the solute and in the transition states that produce the two possible products of the reaction.

CH-π Interactions in Glycan Recognition

Carbohydrate recognition is crucial for biological processes ranging from development to immune system function to host-pathogen interactions. The proteins that bind glycans are faced with a daunting task: to coax these hydrophilic species out of water and into a binding site. Here, we examine the forces underlying glycan recognition by proteins. Our previous bioinformatic study of glycan-binding sites indicated that the most overrepresented side chains are electron-rich aromatic residues, including tyrosine and tryptophan. These findings point to the importance of CH-π interactions for glycan binding. Studies of CH-π interactions show a strong dependence on the presence of an electron-rich π system, and the data indicate binding is enhanced by complementary electronic interactions between the electron-rich aromatic ring and the partial positive charge of the carbohydrate C-H protons. This electronic dependence means that carbohydrate residues with multiple aligned highly polarized C-H bonds, such as β-galactose, form strong CH-π interactions, whereas less polarized residues such as α-mannose do not. This information can guide the design of proteins to recognize sugars and the generation of ligands for proteins, small molecules, or catalysts that bind sugars.

Detection of ligand binding to glycopolymers using saturation transfer difference NMR

We report the use of saturation transfer difference (STD) NMR spectroscopy to observe the interaction of various phenylboronic acids (PBAs) with synthetic glycopolymers presenting galactose and glucose. After optimizing experimental parameters to maximize spin diffusion within the glycopolymers, STD NMR experiments were successfully used to detect binding of PBAs to the polymers. Amplification factor build-up curves in conjunction with differential epitope mapping experiments were used to generate an epitope map for the bound boronic acids. STD NMR was also used to detect the interaction between indole and a galactosylated glycopolymer, providing an indole-based view of this CH-π interaction, a common binding motif in carbohydrate recognition.

Carbohydrate-Protein Interactions: Advances and Challenges

A carbohydrate, also called saccharide in biochemistry, is a biomolecule consisting of carbon (C), hydrogen (H) and oxygen (O) atoms. For example, sugars are low molecular-weight carbohydrates, and starches are high molecular-weight carbohydrates. Carbohydrates are the most abundant organic substances in nature and essential constituents of all living things. Protein-carbohydrate interactions play important roles in many biological processes, such as cell growth, differentiation, and aggregation. They also have broad applications in pharmaceutical drug design. In this review, we will summarize the characteristic features of protein-carbohydrate interactions and review the computational methods for structure prediction, energy calculations, and kinetic studies of protein-carbohydrate complexes. Finally, we will discuss the challenges in this field.

Three-Dimensional Structures of Carbohydrates and Where to Find Them

Analysis and systematization of accumulated data on carbohydrate structural diversity is a subject of great interest for structural glycobiology. Despite being a challenging task, development of computational methods for efficient treatment and management of spatial (3D) structural features of carbohydrates breaks new ground in modern glycoscience. This review is dedicated to approaches of chemo- and glyco-informatics towards 3D structural data generation, deposition and processing in regard to carbohydrates and their derivatives. Databases, molecular modeling and experimental data validation services, and structure visualization facilities developed for last five years are reviewed.

Separation of saccharides using fullerene-bonded silica monolithic columns via π interactions in liquid chromatography

We report on a potential method to separate sugars by using the specific interaction between fullerenes and saccharides in liquid chromatography (LC). Aromatic rings with high electron density are believed to interact strongly with saccharides due to CH–π and/or OH–π interactions. In this study, the fullerene-bonded columns were used to separate saccharides by LC under aqueous conditions. As a result, 2-aminobenzamide-labeled glucose homopolymer (Glcs) was effectively separated by both C60 and C70 columns in the range of Glc-1 to Glc-20 and high blood glucose level being retained in greater quantity. Furthermore, similar separations were identified by LC–mass spectrometry with non-labeled glucose homopolymers. Theoretical study based on molecular dynamics and DFT calculation demonstrated that a supramolecular complex of saccharide–fullerene was formed through CH–π and/or OH–π interactions, and that the interactions between saccharide and fullerene increase with the increase units of the saccharide. Additionally, the C60 column retained disaccharides containing maltose, trehalose, and sucrose. In this case, it was assumed that the retention rates were determined by the difference of the dipole moment in each saccharide. These results suggest that the dipole-induced dipole interaction was dominant, and that maltose—with the higher dipole moment—was more strongly retained compared to other disaccharides having lower dipole moment.

New Insights into the Role of Hydrogen Bonding in Furanoside Binding to Protein

Furanosides have been subjected to extensive studies owing to their inherent flexibility, which is believed to play an important role in the survival and pathogenicity of different disease-causing organisms in the human body. This study reports the binding free energy (ΔG) and specificity of arabinofuranose oligosaccharides to a protein, arabinanase (Arb43A), with the use of potential of mean force (PMF) calculations using the umbrella-sampling simulations. Long molecular dynamics simulations have been carried out to understand intermolecular interactions in the arabinofuranose-protein complex. The PMF for pulling the α-(1 → 5)-linked L-arabinohexaose (ligand) from the protein provides a large free energy of binding, -16.8 kcal/mol. The ΔG of the nonreducing arabinotriose end is found to be -12.6 kcal/mol, while the ΔG of the reducing end is calculated to be -7.7 kcal/mol. In the absence of nonreducing arabinotrioside, the ΔG of the reducing arabinotrioside is -8.5 kcal/mol. Similarly, in the absence of reducing arabinotrioside, the ΔG of the nonreducing arabinotrioside is calculated to be -9.4 kcal/mol. The main contributing factor in the protein-arabinofuranose binding is hydrogen bonding. Acidic amino acid residues, Glu and Asp, with furanosides produce the strongest hydrogen bonding. Araf-A, B, and C construct the reducing arabinotriose, while Araf-D, E, and F construct the nonreducing arabinotriose. Since most of the hydrogen-bonding occupancies belong to Araf-D and Araf-E, the nonreducing arabinotriose is bound to protein more strongly than the reducing arabinotriose. This explains why the reducing arabinotriose can detach from the protein in nature.

Experimental and theoretical study of the role of CH/π interactions in the aminolysis reaction of acetyl galactoside

Molecular recognition of saccharides is a growing field, which has many implications in cancer therapy, drug discovery, and cellular communication among others. The participation of CH/π interactions in this event is well known. Nevertheless, the intrinsic role of CH/π for modulating chemical reactions is still far from being applicable. In this experimental and computational work we have evaluated the participation of CH/π interactions in the aminolysis reaction of acetyl galactoside promoted with different 6-substituted 2(1H)-pyridones. Two features have been incorporated to the promoter molecular structure, on one end the promoting pyridone group and on the other end the recognition moiety, joined together by an alkyne spacer. The small increment in the observed pseudo-first-order rate constant values (k_obs) was related to the stability of the transition state provided by noncovalent interactions, including CH/π interactions. A longer alkyne spacer was necessary to improve the molecular recognition of the galactoside substrate. The trend of the calculated activation energy values (ΔE_RTS) was in good accordance with the experimental rate constant values.

Exploring the Role of Solvent on Carbohydrate-Aryl Interactions by Diffusion NMR-Based Studies

Carbohydrate-protein interactions play an important role in many molecular recognition processes. An exquisite combination of multiple factors favors the interaction of the receptor with one specific type of sugar, whereas others are excluded. Stacking CH-aromatic interactions within the binding site provide a relevant contribution to the stabilization of the resulting sugar-protein complex. Being experimentally difficult to detect and analyze, the key CH-π interaction features have been very often dissected using a variety of techniques and simple model systems. In the present work, diffusion NMR spectroscopy has been employed to separate the components of sugar mixtures in different solvents on the basis of their differential ability to interact through CH-π interactions with one particular aromatic cosolute in solution. The experimental data show that the properties of the solvent did also influence the diffusion behavior of the sugars present in the mixture, inhibiting or improving their separation. Overall, the results showed that, for the considered monosaccharide derivatives, their diffusion coefficient values and, consequently, their apparent molecular sizes and/or shapes depend on the balance between solute/cosolute as well as solute/solvent interactions. Thus, in certain media and in the presence of the aromatic cosolute, the studied saccharides that are more suited to display CH-π interactions exhibited a lower diffusion coefficient than the noncomplexing sugars in the mixture. However, when dissolved in another medium, the interaction with the solvent strongly competes with that of the aromatic cosolute.

Computational Chemistry Tools in Glycobiology: Modelling of Carbohydrate–Protein Interactions

Molecular modelling provides a major impact in the field of glycosciences, helping in the characterisation of the molecular basis of the recognition between lectins from pathogens and human glycoconjugates, and in the design of glycocompounds with anti-infectious properties. The conformational properties of oligosaccharides are complex, and therefore, the simulation of these properties is a challenging task. Indeed, the development of suitable force fields is required for the proper simulation of important problems in glycobiology, such as the interatomic interactions responsible for oligosaccharide and glycoprotein dynamics, including O-linkages in oligo- and polysaccharides, and N- and O-linkages in glycoproteins. The computational description of representative examples is discussed, herein, related to biologically active oligosaccharides and their interaction with lectins and other proteins, and the new routes open for the design of glycocompounds with promising biological activities.

Fluoroacetamide Moieties as NMR Spectroscopy Probes for the Molecular Recognition of GlcNAc-Containing Sugars: Modulation of the CH-π Stacking Interactions by Different Fluorination Patterns

We herein propose the use of fluoroacetamide and difluoroacetamide moieties as sensitive tags for the detection of sugar–protein interactions by simple ¹H and/or ¹⁹F NMR spectroscopy methods. In this process, we have chosen the binding of N,N′‐diacetyl chitobiose, a ubiquitous disaccharide fragment in glycoproteins, by wheat‐germ agglutinin (WGA), a model lectin. By using saturation‐transfer difference (STD)‐NMR spectroscopy, we experimentally demonstrate that, under solution conditions, the molecule that contained the CHF₂CONH‐ moiety is the stronger aromatic binder, followed by the analogue with the CH₂FCONH‐ group and the natural molecule (with the CH₃CONH‐ fragment). In contrast, the molecule with the CF₃CONH‐ isoster displayed the weakest intermolecular interaction (one order of magnitude weaker). Because sugar–aromatic CH–π interactions are at the origin of these observations, these results further contribute to the characterization and exploration of these forces and offer an opportunity to use them to unravel complex recognition processes.

Achiral Molecular Recognition of Aromatic Position Isomers by Polysaccharide-Based CSPs in Relation to Chiral Recognition

Chromatographic separation of several sets of aromatic position isomers on three cellulose- and one amylose-based chiral stationary phases was performed to evaluate the potential of a polysaccharide-based chiral stationary phase (CSP) in the separation of isomeric or closely similar molecules, and to understand the interaction mechanism of this type of CSP with analytes. Their ability of molecular recognition was quite outstanding, but the selection rule was particular to each polysaccharide derivative. In the series of analytes, cellulose tris(4-methylbenzoate) and tris(3,5-dimethylphenylcarbamate) exhibited a contrasting selection rule, and the recognition mechanism was considered based on the computer-simulation of the former polymer.

A Catalytic Role for C-H/π Interactions in Base Excision Repair by Bacillus cereus DNA Glycosylase AlkD

DNA glycosylases protect genomic integrity by locating and excising aberrant nucleobases. Substrate recognition and excision usually take place in an extrahelical conformation, which is often stabilized by π-stacking interactions between the lesion nucleobase and aromatic side chains in the glycosylase active site. Bacillus cereus AlkD is the only DNA glycosylase known to catalyze base excision without extruding the damaged nucleotide from the DNA helix. Instead of contacting the nucleobase itself, the AlkD active site interacts with the lesion deoxyribose through a series of C-H/π interactions. These interactions are ubiquitous in protein structures, but evidence for their catalytic significance in enzymology is lacking. Here, we show that the C-H/π interactions between AlkD and the lesion deoxyribose participate in catalysis of glycosidic bond cleavage. This is the first demonstration of a catalytic role for C-H/π interactions as intermolecular forces important to DNA repair.

Carbohydrate-Aromatic Interactions in Proteins

Protein-carbohydrate interactions play pivotal roles in health and disease. However, defining and manipulating these interactions has been hindered by an incomplete understanding of the underlying fundamental forces. To elucidate common and discriminating features in carbohydrate recognition, we have analyzed quantitatively X-ray crystal structures of proteins with noncovalently bound carbohydrates. Within the carbohydrate-binding pockets, aliphatic hydrophobic residues are disfavored, whereas aromatic side chains are enriched. The greatest preference is for tryptophan with an increased prevalence of 9-fold. Variations in the spatial orientation of amino acids around different monosaccharides indicate specific carbohydrate C-H bonds interact preferentially with aromatic residues. These preferences are consistent with the electronic properties of both the carbohydrate C-H bonds and the aromatic residues. Those carbohydrates that present patches of electropositive saccharide C-H bonds engage more often in CH-π interactions involving electron-rich aromatic partners. These electronic effects are also manifested when carbohydrate-aromatic interactions are monitored in solution: NMR analysis indicates that indole favorably binds to electron-poor C-H bonds of model carbohydrates, and a clear linear free energy relationships with substituted indoles supports the importance of complementary electronic effects in driving protein-carbohydrate interactions. Together, our data indicate that electrostatic and electronic complementarity between carbohydrates and aromatic residues play key roles in driving protein-carbohydrate complexation. Moreover, these weak noncovalent interactions influence which saccharide residues bind to proteins, and how they are positioned within carbohydrate-binding sites.

Carbohydrate-protein interactions: molecular modeling insights

The article reviews the significant contributions to, and the present status of, applications of computational methods for the characterization and prediction of protein-carbohydrate interactions. After a presentation of the specific features of carbohydrate modeling, along with a brief description of the experimental data and general features of carbohydrate-protein interactions, the survey provides a thorough coverage of the available computational methods and tools. At the quantum-mechanical level, the use of both molecular orbitals and density-functional theory is critically assessed. These are followed by a presentation and critical evaluation of the applications of semiempirical and empirical methods: QM/MM, molecular dynamics, free-energy calculations, metadynamics, molecular robotics, and others. The usefulness of molecular docking in structural glycobiology is evaluated by considering recent docking- validation studies on a range of protein targets. The range of applications of these theoretical methods provides insights into the structural, energetic, and mechanistic facets that occur in the course of the recognition processes. Selected examples are provided to exemplify the usefulness and the present limitations of these computational methods in their ability to assist in elucidation of the structural basis underlying the diverse function and biological roles of carbohydrates in their dialogue with proteins. These test cases cover the field of both carbohydrate biosynthesis and glycosyltransferases, as well as glycoside hydrolases. The phenomenon of (macro)molecular recognition is illustrated for the interactions of carbohydrates with such proteins as lectins, monoclonal antibodies, GAG-binding proteins, porins, and viruses.

Convergent and divergent mechanisms of sugar recognition across kingdoms

Protein modules that bind specific oligosaccharides are found across all kingdoms of life from single-celled organisms to man. Different, overlapping and evolving designations for sugar-binding domains in proteins can sometimes obscure common features that often reflect convergent solutions to the problem of distinguishing sugars with closely similar structures and binding them with sufficient affinity to achieve biologically meaningful results. Structural and functional analysis has revealed striking parallels between protein domains with widely different structures and evolutionary histories that employ common solutions to the sugar recognition problem. Recent studies also demonstrate that domains descended from common ancestors through divergent evolution appear more widely across the kingdoms of life than had previously been recognized.

Conformational dynamics of dry lamellar crystals of sugar based lipids: an atomistic simulation study

The rational design of a glycolipid application (e.g. drug delivery) with a tailored property depends on the detailed understanding of its structure and dynamics. Because of the complexity of sugar stereochemistry, we have undertaken a simulation study on the conformational dynamics of a set of synthetic glycosides with different sugar groups and chain design, namely dodecyl β-maltoside, dodecyl β-cellobioside, dodecyl β-isomaltoside and a C12C10 branched β-maltoside under anhydrous conditions. We examined the chain structure in detail, including the chain packing, gauche/trans conformations and chain tilting. In addition, we also investigated the rotational dynamics of the headgroup and alkyl chains. Monoalkylated glycosides possess a small amount of gauche conformers (∼20%) in the hydrophobic region of the lamellar crystal (LC) phase. In contrast, the branched chain glycolipid in the fluid Lα phase has a high gauche population of up to ∼40%. Rotational diffusion analysis reveals that the carbons closest to the headgroup have the highest correlation times. Furthermore, its value depends on sugar type, where the rotational dynamics of an isomaltose was found to be 11-15% and more restrained near the sugar, possibly due to the chain disorder and partial inter-digitation compared to the other monoalkylated lipids. Intriguingly, the present simulation demonstrates the chain from the branched glycolipid bilayer has the ability to enter into the hydrophilic region. This interesting feature of the anhydrous glycolipid bilayer simulation appears to arise from a combination of lipid crowding and the amphoteric nature of the sugar headgroups.

Self-assembly of proteins into a three-dimensional multilayer system: investigation of the surface of the human fungal pathogen Aspergillus fumigatus

Hydrophobins are small surface active proteins that fulfil a wide spectrum of functions in fungal growth and development. The human fungal pathogen Aspergillus fumigatus expresses RodA hydrophobins that self-assemble on the outer conidial surface into tightly organized nanorods known as rodlets. AFM investigation of the conidial surface allows us to evidence that RodA hydrophobins self-assemble into rodlets through bilayers. Within bilayers, hydrophilic domains of hydrophobins point inward, thus making a hydrophilic core, while hydrophobic domains point outward. AFM measurements reveal that several rodlet bilayers are present on the conidial surface thus showing that proteins self-assemble into a complex three-dimensional multilayer system. The self-assembly of RodA hydrophobins into rodlets results from attractive interactions between stacked β-sheets, which conduct to a final linear cross-β spine structure. A Monte Carlo simulation shows that anisotropic interactions are the main driving forces leading the hydrophobins to self-assemble into parallel rodlets, which are further structured in nanodomains. Taken together, these findings allow us to propose a mechanism, which conducts RodA hydrophobins to a highly ordered rodlet structure. The mechanism of hydrophobin assembly into rodlets offers new prospects for the development of more efficient strategies leading to disruption of rodlet formation allowing a rapid detection of the fungus by the immune system.

A carbohydrate-anion recognition system in aprotic solvents

A carbohydrate–anion recognition system in nonpolar solvents is reported, in which complexes form at the B-faces of β-d-pyranosides with H1-, H3-, and H5-cis patterns similar to carbohydrate–π interactions. The complexation effect was evaluated for a range of carbohydrate structures; it resulted in either 1:1 carbohydrate–anion complexes, or 1:2 complex formation depending on the protection pattern of the carbohydrate. The interaction was also evaluated with different anions and solvents. In both cases it resulted in significant binding differences. The results indicate that complexation originates from van der Waals interactions or weak CH⋅⋅⋅A⁻ hydrogen bonds between the binding partners and is related to electron-withdrawing groups of the carbohydrates as well as increased hydrogen-bond-accepting capability of the anions.

Protein-carbohydrate interactions studied by NMR: from molecular recognition to drug design

Diseases that result from infection are, in general, a consequence of specific interactions between a pathogenic organism and the cells. The study of host-pathogen interactions has provided insights for the design of drugs with therapeutic properties. One area that has proved to be promising for such studies is the constituted by carbohydrates which participate in biological processes of paramount importance. On the one hand, carbohydrates have shown to be information carriers with similar, if not higher, importance than traditionally considered carriers as amino acids and nucleic acids. On the other hand, the knowledge on molecular recognition of sugars by lectins and other carbohydrate-binding proteins has been employed for the development of new biomedical strategies. Biophysical techniques such as X-Ray crystallography and NMR spectroscopy lead currently the investigation on this field. In this review, a description of traditional and novel NMR methodologies employed in the study of sugar-protein interactions is briefly presented in combination with a palette of NMR-based studies related to biologically and/or pharmaceutically relevant applications.

Structural and energetic basis of carbohydrate-aromatic packing interactions in proteins

Carbohydrate-aromatic interactions mediate many biological processes. However, the structure-energy relationships underpinning direct carbohydrate-aromatic packing interactions in aqueous solution have been difficult to assess experimentally and remain elusive. Here, we determine the structures and folding energetics of chemically synthesized glycoproteins to quantify the contributions of the hydrophobic effect and CH-π interactions to carbohydrate-aromatic packing interactions in proteins. We find that the hydrophobic effect contributes significantly to protein-carbohydrate interactions. Interactions between carbohydrates and aromatic amino acid side chains, however, are supplemented by CH-π interactions. The strengths of experimentally determined carbohydrate CH-π interactions do not correlate with the electrostatic properties of the involved aromatic residues, suggesting that the electrostatic component of CH-π interactions in aqueous solution is small. Thus, tight binding of carbohydrates and aromatic residues is driven by the hydrophobic effect and CH-π interactions featuring a dominating dispersive component.

Carbohydrate-aromatic interactions

The recognition of saccharides by proteins has far reaching implications in biology, technology, and drug design. Within the past two decades, researchers have directed considerable effort toward a detailed understanding of these processes. Early crystallographic studies revealed, not surprisingly, that hydrogen-bonding interactions are usually involved in carbohydrate recognition. But less expectedly, researchers observed that despite the highly hydrophilic character of most sugars, aromatic rings of the receptor often play an important role in carbohydrate recognition. With further research, scientists now accept that noncovalent interactions mediated by aromatic rings are pivotal to sugar binding. For example, aromatic residues often stack against the faces of sugar pyranose rings in complexes between proteins and carbohydrates. Such contacts typically involve two or three CH groups of the pyranoses and the π electron density of the aromatic ring (called CH/π bonds), and these interactions can exhibit a variety of geometries, with either parallel or nonparallel arrangements of the aromatic and sugar units. In this Account, we provide an overview of the structural and thermodynamic features of protein-carbohydrate interactions, theoretical and experimental efforts to understand stacking in these complexes, and the implications of this understanding for chemical biology. The interaction energy between different aromatic rings and simple monosaccharides based on quantum mechanical calculations in the gas phase ranges from 3 to 6 kcal/mol range. Experimental values measured in water are somewhat smaller, approximately 1.5 kcal/mol for each interaction between a monosaccharide and an aromatic ring. This difference illustrates the dependence of these intermolecular interactions on their context and shows that this stacking can be modulated by entropic and solvent effects. Despite their relatively modest influence on the stability of carbohydrate/protein complexes, the aromatic platforms play a major role in determining the specificity of the molecular recognition process. The recognition of carbohydrate/aromatic interactions has prompted further analysis of the properties that influence them. Using a variety of experimental and theoretical methods, researchers have worked to quantify carbohydrate/aromatic stacking and identify the features that stabilize these complexes. Researchers have used site-directed mutagenesis, organic synthesis, or both to incorporate modifications in the receptor or ligand and then quantitatively analyzed the structural and thermodynamic features of these interactions. Researchers have also synthesized and characterized artificial receptors and simple model systems, employing a reductionistic chemistry-based strategy. Finally, using quantum mechanics calculations, researchers have examined the magnitude of each property's contribution to the interaction energy.

Caffeine and sugars interact in aqueous solutions: a simulation and NMR study

Molecular dynamics simulations were carried out on several systems of caffeine interacting with simple sugars. These included a single caffeine molecule in a 3 m solution of α-D-glucopyranose, at a caffeine concentration of 0.083 m, a single caffeine in a 3 m solution of β-D-glucopyranose, and a single caffeine molecule in a 1.08 m solution of sucrose (table sugar). Parallel nuclear magnetic resonance titration experiments were carried out on the same solutions under similar conditions. Consistent with previous thermodynamic experiments, the sugars were found to have an affinity for the caffeine molecules in both the simulations and experiments, and the binding in these complexes occurs by face-to-face stacking of the hydrophobic triad of protons of the pyranose rings against the caffeine face, rather than by hydrogen bonding. For the disaccharide, the binding occurs via stacking of the glucose ring against the caffeine, with a lesser affinity for the fructose observed. These findings are consistent with the association being driven by hydrophobic hydration and are similar to the previously observed binding of glucose rings to various other planar molecules, including indole, serotonin, and phenol.

Biochemical and structural characterization of the complex agarolytic enzyme system from the marine bacterium Zobellia galactanivorans

Zobellia galactanivorans is an emerging model bacterium for the bioconversion of algal biomass. Notably, this marine Bacteroidetes possesses a complex agarolytic system comprising four β-agarases and five β-porphyranases, all belonging to the glycoside hydrolase family 16. Although β-agarases are specific for the neutral agarobiose moieties, the recently discovered β-porphyranases degrade the sulfated polymers found in various quantities in natural agars. Here, we report the biochemical and structural comparison of five β-porphyranases and β-agarases from Z. galactanivorans. The respective degradation patterns of two β-porphyranases and three β-agarases are analyzed by their action on defined hybrid oligosaccharides. In light of the high resolution crystal structures, the biochemical results allowed a detailed mapping of substrate specificities along the active site groove of the enzymes. Although PorA displays a strict requirement for C6-sulfate in the -2- and +1-binding subsites, PorB tolerates the presence of 3-6-anhydro-l-galactose in subsite -2. Both enzymes do not accept methylation of the galactose unit in the -1 subsite. The β-agarase AgaD requires at least four consecutive agarose units (DP8) and is highly intolerant to modifications, whereas for AgaB oligosaccharides containing C6-sulfate groups at the -4, +1, and +3 positions are still degraded. Together with a transcriptional analysis of the expression of these enzymes, the structural and biochemical results allow proposition of a model scheme for the agarolytic system of Z. galactanivorans.

Multiple functions of aromatic-carbohydrate interactions in a processive cellulase examined with molecular simulation

Proteins employ aromatic residues for carbohydrate binding in a wide range of biological functions. Glycoside hydrolases, which are ubiquitous in nature, typically exhibit tunnels, clefts, or pockets lined with aromatic residues for processing carbohydrates. Mutation of these aromatic residues often results in significant activity differences on insoluble and soluble substrates. However, the thermodynamic basis and molecular level role of these aromatic residues remain unknown. Here, we calculate the relative ligand binding free energy by mutating tryptophans in the Trichoderma reesei family 6 cellulase (Cel6A) to alanine. Removal of aromatic residues near the catalytic site has little impact on the ligand binding free energy, suggesting that aromatic residues immediately upstream of the active site are not directly involved in binding, but play a role in the glucopyranose ring distortion necessary for catalysis. Removal of aromatic residues at the entrance and exit of the Cel6A tunnel, however, dramatically impacts the binding affinity, suggesting that these residues play a role in chain acquisition and product stabilization, respectively. The roles suggested from differences in binding affinity are confirmed by molecular dynamics and normal mode analysis. Surprisingly, our results illustrate that aromatic-carbohydrate interactions vary dramatically depending on the position in the enzyme tunnel. As aromatic-carbohydrate interactions are present in all carbohydrate-active enzymes, these results have implications for understanding protein structure-function relationships in carbohydrate metabolism and recognition, carbon turnover in nature, and protein engineering strategies for biomass utilization. Generally, these results suggest that nature employs aromatic-carbohydrate interactions with a wide range of binding affinities for diverse functions.