The CH-π Interaction in Protein-Carbohydrate Binding: Bioinformatics and In Vitro Quantification

A Computational and Chemical Design Strategy for Manipulating Glycan-Protein Recognition

Glycans are complex biomolecules that encode rich information and regulate various biological processes, such as fertilization, host-pathogen binding, and immune recognition, through interactions with glycan-binding proteins. A key driving force for glycan-protein recognition is the interaction between the π electron density of aromatic amino acid side chains and polarized C─H groups of the pyranose (termed the CH-π interaction). However, the relatively weak binding affinity between glycans and proteins has hindered the application of glycan detection and imaging. Here, computational modeling and molecular dynamics simulations are employed to design a chemical strategy that enhances the CH-π interaction between glycans and proteins by genetically incorporating electron-rich tryptophan derivatives into a lectin PhoSL, which specifically recognizes core fucosylated N-linked glycans. This significantly enhances the binding affinity of PhoSL with the core fucose ligand and enables sensitive detection and imaging of core fucosylated glycans in vitro and in xenograft tumors in mice. Further, the study showed that this strategy is applicable to improve the binding affinity of GafD lectin for N-acetylglucosamine-containing glycans. The approach thus provides a general and effective way to manipulate glycan-protein recognition for glycoscience applications.

GPER binding site detection and description: A flavonoid-based docking and molecular dynamics simulations study

Flavonoids, a phenolic compounds class widely distributed in the plant kingdom, have attracted much interest for their implications on several health and disease processes. Usually, the consumption of this type of compounds is approximately 1 g/d, primarily obtained from cereals, chocolate, and dry legumes ensuring its beneficial role in maintaining the homeostasis of the human body. In this context, in cancer disease prominent data points to the role of flavonoids as adjuvant treatment aimed at the regression of the disease. GPER, an estrogen receptor on the cell surface, has been postulated as a probable orchestrator of the beneficial effects of several flavonoids through modulation/inhibition of various mechanisms that lead to cancer progression. Therefore, applying pocket and cavity protein detection and docking and molecular dynamics simulations (MD), we generate, from a cluster composed of 39 flavonoids, crucial insights into the potential role as GPER ligands, of Puerarin, Isoquercetin, Kaempferol 3-O-glucoside and Petunidin 3-O-glucoside, aglycones whose sugar moiety delimits a new described sugar-acceptor sub-cavity into the cavity binding site on the receptor, as well as of the probable activation mechanism of the receptor and the pivotal residues involved in it. Altogether, our results shed light on the potential use of the aforementioned flavonoids as GPER ligands and for further evaluations in in vitro and in vivo assays to elucidate their probable anti-cancer activity.

Crystal structure of mango α1,3/α1,4-fucosyltransferase elucidates unique elements that regulate Lewis A-dominant oligosaccharide assembly

In various organisms, α1,3/α1,4-fucosyltransferases (CAZy GT10 family enzymes) mediate the assembly of type I (Galβ1,3GlcNAc) and/or type II (Galβ1,4GlcNAc)-based Lewis structures that are widely distributed in glycoconjugates. Unlike enzymes of other species, plant orthologues show little fucosyltransferase activity for type II-based glycans and predominantly catalyze the assembly of the Lewis A structure [Galβ1,3(Fucα1,4)GlcNAc] on the type I disaccharide unit of their substrates. However, the structural basis underlying this unique substrate selectivity remains elusive. In this study, we investigated the structure-function relationship of MiFUT13A, a mango α1,3/α1,4-fucosyltransferase. The prepared MiFUT13A displayed distinct α1,4-fucosyltransferase activity. Consistent with the enzymatic properties of this molecule, X-ray crystallography revealed that this enzyme has a typical GT-B fold-type structure containing a set of residues that are responsible for its SN2-like catalysis. Site-directed mutagenesis and molecular docking analyses proposed a rational binding mechanism for type I oligosaccharides. Within the catalytic cleft, the pocket surrounding Trp121 serves as a binding site, anchoring the non-reducing terminal β1,3-galactose that belongs to the type I disaccharide unit. Furthermore, Glu177 was postulated to function as a general base catalyst through its interaction with the 4-hydroxy group of the acceptor N-acetylglucosamine residue. Adjacent residues, specifically Thr120, Thr157 and Asp175 were speculated to assist in binding of the reducing terminal residues. Intriguingly, these structural elements were not fully conserved in mammalian orthologue which also shows predominant α1,4-fucosyltransferase activity. In conclusion, we have proposed that MiFUT13A generates the Lewis A structure on type I glycans through a distinct mechanism, divergent from that of mammalian enzymes.

C-H···π interactions disrupt electrostatic interactions between non-aqueous electrolytes to increase solubility

Grid-scale energy storage applications, such as redox flow batteries, rely on the solubility of redox-active organic molecules. Although redox-active pyridiniums exhibit exceptional persistence in multiple redox states at low potentials (desirable properties for energy storage applications), their solubility in non-aqueous media remains low, and few practical molecular design strategies exist to improve solubility. Here we convey the extent to which discrete, attractive interactions between C-H groups and π electrons of an aromatic ring (C-H···π interactions) can describe the solubility of N-substituted pyridinium salts in a non-aqueous solvent. We find a direct correlation between the number of C-H···π interactions for each pyridinium salt and its solubility in acetonitrile. The correlation presented in this work highlights a consequence of disrupting strong electrostatic interactions with weak dispersion interactions, showing how minimal structural change can dramatically impact pyridinium solubility.

Interacting Quantum Atoms and Multipolar Electrostatic Study of XH···π Interactions

The interaction energies of nine XH···π (X = C, N, and O) benzene-containing van der Waals complexes were analyzed, at the atomic and fragment levels, using QTAIM multipolar electrostatics and the energy partitioning method interacting quantum atoms/fragment (IQA/IQF). These descriptors were paired with the relative energy gradient method, which solidifies the connection between quantum mechanical properties and chemical interpretation. This combination provides a precise understanding, both qualitative and quantitative, of the nature of these interactions, which are ubiquitous in biochemical systems. The formation of the OH···π and NH···π systems is electrostatically driven, with the Qzz component of the quadrupole moment of the benzene carbons interacting with the charges of X and H in XH. There is the unexpectedly intramonomeric role of X-H (X = O, N) where its electrostatic energy helps the formation of the complex and its covalent energy thwarts it. However, the CH···π interaction is governed by exchange-correlation energies, thereby establishing a covalent character, as opposed to the literature's designation as a noncovalent interaction. Moreover, dispersion energy is relevant, statically and in absolute terms, but less relevant compared to other energy components in terms of the formation of the complex. Multipolar electrostatics are similar across all systems.

Mechanisms of interactions in pattern-recognition of common glycostructures across pectin-derived heteropolysaccharides by Toll-like receptor 4

Complex pectin, originating from terrestrial plant cell walls has been attracting research attention as a promising source of a new innate immune modulator. Numerous bioactive polysaccharides associated with pectin are newly reported every year, but the general mechanism of their immunological action remains unclear owing to the complexity and heterogeneity of pectin. Herein, we systematically investigated the interactions in pattern-recognition for common glycostructures of pectic heteropolysaccharides (HPSs) by Toll-like receptors (TLRs). The compositional similarity of glycosyl residues derived from pectic HPS was confirmed by conducting systematic reviews, leading to molecular modeling of representative pectic segments. Via structural investigation, the inner concavity of leucine-rich repeats of TLR4 was predicted to act as a binding motif for carbohydrate recognition, and subsequent simulations predicted the binding modes and conformations. We experimentally demonstrated that pectic HPS exhibits the non-canonical and multivalent binding aspects for TLR4 resulting in receptor activation. Furthermore, we showed that pectic HPSs were selectively clustered with TLR4 during endocytosis, inducing downstream signals to cause phenotypic activation of macrophages. Overall, we have presented a better explanation for the pattern recognition of pectic HPS and further proposed an approach to understand the interaction between complex carbohydrates and proteins.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Martini 3 Coarse-Grained Force Field for Carbohydrates

The Martini 3 force field is a full reparametrization of the Martini coarse-grained model for biomolecular simulations. Due to the improved interaction balance, it allows for a more accurate description of condensed phase systems. In the present work, we develop a consistent strategy to parametrize carbohydrate molecules accurately within the framework of Martini 3. In particular, we develop a canonical mapping scheme which decomposes arbitrarily large carbohydrates into a limited number of fragments. Bead types for these fragments have been assigned by matching physicochemical properties of mono- and disaccharides. In addition, guidelines for assigning bonds, angles, and dihedrals were developed. These guidelines enable a more accurate description of carbohydrate conformations than in the Martini 2 force field. We show that models obtained with this approach are able to accurately reproduce osmotic pressures of carbohydrate water solutions. Furthermore, we provide evidence that the model differentiates correctly the solubility of the polyglucoses dextran (water-soluble) and cellulose (water insoluble but soluble in ionic liquids). Finally, we demonstrate that the new building blocks can be applied to glycolipids. We show they are able to reproduce membrane properties and induce binding of peripheral membrane proteins. These test cases demonstrate the validity and transferability of our approach.

Cycloalkyl Groups as Building Blocks of Artificial Carbohydrate Receptors: Studies with Macrocycles Bearing Flexible Side-Arms

The cyclopentyl group was expected to act as a building block for artificial carbohydrate receptors and to participate in van der Waals contacts with the carbohydrate substrate in a similar way as observed for the pyrrolidine ring of proline in the crystal structures of protein-carbohydrate complexes. Systematic binding studies with a series of 1,3,5-trisubstituted 2,4,6-triethylbenzenes bearing various cycloalkyl groups as recognition units provided indications of the involvement of these groups in the complexation process and showed the influence of the ring size on the receptor efficiency. Representatives of compounds that exhibit a macrocyclic backbone and flexible side arms were now chosen as further model systems to investigate whether the previously observed effects represent a general trend. Binding studies with these macrocycles towards β-D-glucopyranoside, an all-equatorial substituted carbohydrate substrate, included 1H NMR spectroscopic titrations and microcalorimetric investigations. The performed studies confirmed the previously observed tendency and showed that the compound bearing cyclohexyl groups displays the best binding properties.

The evolutionary advantage of an aromatic clamp in plant family 3 glycoside exo-hydrolases

In the barley β-D-glucan glucohydrolase, a glycoside hydrolase family 3 (GH3) enzyme, the Trp286/Trp434 clamp ensures β-D-glucosides binding, which is fundamental for substrate hydrolysis during plant growth and development. We employ mutagenesis, high-resolution X-ray crystallography, and multi-scale molecular modelling methods to examine the binding and conformational behaviour of isomeric β-D-glucosides during substrate-product assisted processive catalysis that operates in GH3 hydrolases. Enzyme kinetics reveals that the W434H mutant retains broad specificity, while W434A behaves as a strict (1,3)-β-D-glucosidase. Investigations of reactant movements on the nanoscale reveal that processivity is sensitive to mutation-specific alterations of the tryptophan clamp. While wild-type and W434H utilise a lateral cavity for glucose displacement and sliding of (1,3)-linked hydrolytic products through the catalytic site without dissociation, consistent with their high hydrolytic rates, W434A does not adopt processive catalysis. Phylogenomic analyses of GH3 hydrolases disclose the evolutionary advantage of the tryptophan clamp that confers broad specificity, high catalytic efficiency, and processivity.

Novel Synthesis Methods of New Imidazole-Containing Coordination Compounds Tc(IV, V, VII)-Reaction Mechanism, Xrd and Hirshfeld Surface Analysis

In this work, we have proposed two new methods for the synthesis of [TcO₂L₄]⁺ (where L = imidazole (Im), methylimidazole (MeIm)) complexes using thiourea (Tu) and Sn(II) as the reducing agents. The main and by-products of the reactions were determined, and possible reaction mechanisms were proposed. We have shown that the reduction of Tc(VII) with thiourea is accompanied by the formation of the Tc(III) intermediate and further oxidation to Tc(V). The reaction conditions’ changing can lead to the formation of Tc(VII) and Tc(IV) salts. Seven new crystal structures are described in this work: Tc(V) complexes, salts with Tc(VII) and Tc(IV) anions. For the halide salts of Tu the cell parameters were determined. In all of the obtained compounds, except for [TcO₂(MeIm)₄]TcO₄, there are π–stacking interactions between the aromatic rings. An increase in the anion size lead to weakening of the intermolecular interactions. The halogen bonds and anion-π interactions were also found in the hexahalide-containing compounds. The Hirshfeld surface analysis showed that the main contribution to the crystal packing is created by the van der Waals interactions of the H···H type (42.5–55.1%), H···C/C···H (17.7–21.3%) and hydrogen bonds, which contribute 15.7–25.3% in total.

The Role of Tryptophan in π Interactions in Proteins: An Experimental Approach

In proteins, the amino acids Phe, Tyr, and especially Trp are frequently involved in π interactions such as π–π, cation−π, and CH−π bonds. These interactions are often crucial for protein structure and protein–ligand binding. A powerful means to study these interactions is progressive fluorination of these aromatic residues to modulate the electrostatic component of the interaction. However, to date no protein expression platform is available to produce milligram amounts of proteins labeled with such fluorinated amino acids. Here, we present a Lactococcus lactis Trp auxotroph-based expression system for efficient incorporation (≥95%) of mono-, di-, tri-, and tetrafluorinated, as well as a methylated Trp analog. As a model protein we have chosen LmrR, a dimeric multidrug transcriptional repressor protein from L. lactis. LmrR binds aromatic drugs, like daunomycin and riboflavin, between Trp96 and Trp96′ in the dimer interface. Progressive fluorination of Trp96 decreased the affinity for the drugs 6- to 70-fold, clearly establishing the importance of electrostatic π–π interactions for drug binding. Presteady state kinetic data of the LmrR–drug interaction support the enthalpic nature of the interaction, while high resolution crystal structures of the labeled protein–drug complexes provide for the first time a structural view of the progressive fluorination approach. The L. lactis expression system was also used to study the role of Trp68 in the binding of riboflavin by the membrane-bound riboflavin transport protein RibU from L. lactis. Progressive fluorination of Trp68 revealed a strong electrostatic component that contributed 15–20% to the total riboflavin-RibU binding energy.

Transcriptomic and Drug Discovery Analyses Reveal Natural Compounds Targeting the KDM4 Subfamily as Promising Adjuvant Treatments in Cancer

KDM4 proteins are a subfamily of histone demethylases that target the trimethylation of lysines 9 and 36 of histone H3, which are associated with transcriptional repression and elongation respectively. Their deregulation in cancer may lead to chromatin structure alteration and transcriptional defects that could promote malignancy. Despite that KDM4 proteins are promising drug targets in cancer therapy, only a few drugs have been described as inhibitors of these enzymes, while studies on natural compounds as possible inhibitors are still needed. Natural compounds are a major source of biologically active substances and many are known to target epigenetic processes such as DNA methylation and histone deacetylation, making them a rich source for the discovery of new histone demethylase inhibitors. Here, using transcriptomic analyses we determined that the KDM4 family is deregulated and associated with a poor prognosis in multiple neoplastic tissues. Also, by molecular docking and molecular dynamics approaches, we screened the COCONUT database to search for inhibitors of natural origin compared to FDA-approved drugs and DrugBank databases. We found that molecules from natural products presented the best scores in the FRED docking analysis. Molecules with sugars, aromatic rings, and the presence of OH or O- groups favor the interaction with the active site of KDM4 subfamily proteins. Finally, we integrated a protein-protein interaction network to correlate data from transcriptomic analysis and docking screenings to propose FDA-approved drugs that could be used as multitarget therapies or in combination with the potential natural inhibitors of KDM4 enzymes. This study highlights the relevance of the KDM4 family in cancer and proposes natural compounds that could be used as potential therapies.

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.

Prebiotic Peptides Based on the Glycocodon Theory Analyzed with FRET

In modern protein–carbohydrate interactions, carbohydrate–aromatic contact with CH–π interactions are used. Currently, they are considered driving forces of this complexation. In these contacts, tryptophan, tyrosine, and histidine are preferred. In this study, we focus on primary prebiotic chemistry when only glycine, alanine, aspartic acid, and valine are available in polypeptides. In this situation, when the aromatic acids are not available, hydrogen-bonding aspartic acid must be used for monosaccharide complexation. It is shown here that (DAA)n polypeptides play important roles in primary “protein”–glucose recognition, that (DGG)n plays an important role in “protein”–ribose recognition, and that (DGA)n plays an important role in “protein”–galactose recognition. Glucose oxidase from Aspergillus niger, which still has some ancient prebiotic sequences, is chosen here as an example for discussion.

Visualization of hydrogen atoms in a perdeuterated lectin-fucose complex reveals key details of protein-carbohydrate interactions

Carbohydrate-binding proteins from pathogenic bacteria and fungi have been shown to be implicated in various pathological processes, where they interact with glycans present on the surface of the host cells. These interactions are part of the initial processes of infection of the host and are very important to study at the atomic level. Here, we report the room temperature neutron structures of PLL lectin from Photorhabdus laumondii in its apo form and in complex with deuterated L-fucose, which is, to our knowledge, the first neutron structure of a carbohydrate-binding protein in complex with a fully deuterated carbohydrate ligand. A detailed structural analysis of the lectin-carbohydrate interactions provides information on the hydrogen bond network, the role of water molecules, and the extent of the CH-π stacking interactions between fucose and the aromatic amino acids in the binding site.

Anatomy of noncovalent interactions between the nucleobases or ribose and π-containing amino acids in RNA-protein complexes

A set of >300 nonredundant high-resolution RNA–protein complexes were rigorously searched for π-contacts between an amino acid side chain (W, H, F, Y, R, E and D) and an RNA nucleobase (denoted π–π interaction) or ribose moiety (denoted sugar–π). The resulting dataset of >1500 RNA–protein π-contacts were visually inspected and classified based on the interaction type, and amino acids and RNA components involved. More than 80% of structures searched contained at least one RNA–protein π-interaction, with π–π contacts making up 59% of the identified interactions. RNA–protein π–π and sugar–π contacts exhibit a range in the RNA and protein components involved, relative monomer orientations and quantum mechanically predicted binding energies. Interestingly, π–π and sugar–π interactions occur more frequently with RNA (4.8 contacts/structure) than DNA (2.6). Moreover, the maximum stability is greater for RNA–protein contacts than DNA–protein interactions. In addition to highlighting distinct differences between RNA and DNA–protein binding, this work has generated the largest dataset of RNA–protein π-interactions to date, thereby underscoring that RNA–protein π-contacts are ubiquitous in nature, and key to the stability and function of RNA–protein complexes.

Binding modes of methyl α-d-glucopyranoside to an artificial receptor in crystalline complexes

Complexes formed between methyl α-d-glucopyranoside and an artificial receptor represent a valuable source of information about the basic molecular features of carbohydrate recognition.