# Glycan-protein interactions

Glycan-Protein interactions represent a class of biomolecular interactions that occur between free or protein-bound glycans and their cognate binding partners. Intramolecular glycan-protein (protein-glycan) interactions occur between glycans and proteins that they are covalently attached to. Together with protein-protein interactions, they form a mechanistic basis for many essential cell processes, especially for cell-cell interactions and host-cell interactions. For instance, SARS-CoV-2, the causative agent of COVID-19, employs its extensively glycosylated spike (S) protein to bind to the ACE2 receptor, allowing it to enter host cells. The spike protein is a trimeric structure, with each subunit containing 22 N-glycosylation sites, making it an attractive target for vaccine search.
Glycosylation, i.e., the addition of glycans (a generic name for monosaccharides and oligosaccharides) to a protein, is one of the major post-translational modification of proteins contributing to the enormous biological complexity of life. Indeed, three different hexoses could theoretically produce from 1056 to 27,648 unique trisaccharides in contrast to only 6 peptides or oligonucleotides formed from 3 amino acids or 3 nucleotides respectively. In contrast to template-driven protein biosynthesis, the "language" of glycosylation is still unknown, making glycobiology a hot topic of current research given their prevalence in living organisms.
The study of glycan-protein interactions provides insight into the mechanisms of cell-signaling and allows to create better-diagnosing tools for many diseases, including cancer. Indeed, there are no known types of cancer that do not involve erratic patterns of protein glycosylation.

## Thermodynamics of Binding

The binding of glycan-binding proteins (GBPs) to glycans could be modeled with simple equilibrium. Denoting glycans as $G$ and proteins as $P$:
$Protein(P)+Glycan(G)\rightleftharpoons PG$
With an associated equilibrium constant of
$K_{a}={\frac {[PG]}{[P][G]}}$
Which is rearranged to give dissociation constant $K_{d}$ following biochemical conventions:
$K_{d}={\frac {[P][G]}{[PG]}}$
Given that many GBPs exhibit multivalency, this model may be expanded to account for multiple equilibria:
$P+G\rightleftharpoons PG$
$PG+G\rightleftharpoons PG_{2}$
$\dots$
$PG_{n-1}+G\rightleftharpoons PG_{n}$
Denoting cumulative equilibrium of binding with $i$ ligands as
$P+iG\rightleftharpoons PG_{i}$
With corresponding equilibrium constant:
$\beta _{i}={\frac {[PG_{i}]}{[P][G]^{i}}}$
And writing material balance for protein ($c_{P}$ denotes the total concentration of protein):
$c_{P}=[P]+[PG]+\dots +[PG_{n}]$
Expressing the terms through an equilibrium constant, a final result is found:
$c_{P}=[P](1+\beta _{1}[G]+\dots +\beta _{n}[G]^{n}$
The concentration of free protein is, thus:
$[P]={\frac {c_{P}}{1+\sum _{i=1}^{n}{\beta _{i}[G]^{i}}}}$
If $n=1$, i.e. there is only one carbohydrate receptor domain, the equation reduces to
$[P]={\frac {c_{P}}{1+\beta _{1}[G]}}$
With increasing $i$ the concentration of free protein decreases; hence, the apparent $K_{D}$ decreases too.

## Binding with aromatic rings

The chemical intuition suggests that the glycan-binding sites may be enriched in polar amino acid residues that form non-covalent interactions, such as hydrogen bonds, with polar carbohydrates. Indeed, a statistical analysis of carbohydrate-binding pockets shows that aspartic acid and asparagine residues are present twice as often as would be predicted by chance. Surprisingly, there is an even stronger preference for aromatic amino acids: tryptophan has a 9-fold increase in prevalence, tyrosine a 3-fold one, and histidine a 2-fold increase. It has been shown that the underlying force is the $CH-\pi$ interaction between the aromatic $\pi$ system and the $C-H$ in carbohydrate as shown in Figure 1. The $CH-\pi$ interaction is identified if the $\theta \leqslant 40$°, the $CH-\pi$ distance (distance from $C$ to $X$) is less than 4.5Å.

### Effects of stereochemistry

This $CH-\pi$ interaction strongly depends on the stereochemistry of the carbohydrate molecule. For example, consider the top ($\beta$) and bottom ($\alpha$) faces of $\beta$-D-Glucose and $\beta$-D-Galactose. It has been shown that a single change in the stereochemistry at C4 carbon shifts preference for aromatic residues from $\beta$ side (2.7 fold preference for glucose) to the $\alpha$ side (14 fold preference for galactose).

### Effects of electronics

The comparison of electrostatic surface potentials (ESPs) of aromatic rings in tryptophan, tyrosine, phenylalanine, and histidine suggests that electronic effects also play a role in the binding to glycans (see Figure 2). After normalizing the electron densities for surface area, the tryptophan still remains the most electron rich acceptor of $CH-\pi$ interactions, suggesting a possible reason for its 9-fold prevalence in carbohydrate binding pockets. Overall, the electrostatic potential maps follow the prevalence trend of ${\ce {Trp >> Tyr > (Phe) > His}}$.

## Carbohydrate-binding partners

There are many proteins capable of binding to glycans, including lectins, antibodies, microbial adhesins, viral agglutinins, etc.

### Lectins

Lectins is a generic name for proteins with carbohydrate-recognizing domains (CRD). Although it became almost synonymous with glycan-binding proteins, it does not include antibodies which also belong to the class.
Lectins found in plants and fungi cells have been extensively used in research as a tool to detect, purify, and analyze glycans. However, useful lectins usually have sub-optimal specificities. For instance, Ulex europaeus agglutinin-1 (UEA-1), a plant-extracted lectin capable of binding to human blood type O antigen, can also bind to unrelated glycans such as 2'-fucosyllactose, GalNAcα1-4(Fucα1-2)Galβ1-4GlcNAc, and Lewis-Y antigen.

### Antibodies

Although antibodies exhibit nanomolar affinities toward protein antigens, the specificity against glycans is very limited. In fact, available antibodies may bind only <4% of the 7000 mammalian glycan antigens; moreover, most of those antibodies have low affinity and exhibit cross-reactivity.

### Lambodies

In contrast with jawed vertebrates whose immunity is based on variable, diverse, and joining gene segments (VDJs) of immunoglobulins, the jawless invertebrates, such as lamprey and hagfish, create a receptor diversity by somatic DNA rearrangement of leucine-rich repeat (LRR) modules that are incorporate in *vlr* genes (variable leukocyte receptors). Those LRR form 3D structures resembling curved solenoids that selectively bind specific glycans.
A study from University of Maryland has shown that lamprey antibodies (lambodies) could selectively bind to tumor-associated carbohydrate antigens (such as Tn and TF$\alpha$) at nanomolar affinities. The T-nouvelle antigen (Tn) and TF$\alpha$ are present in proteins in as much as 90% of different cancer cells after post-translational modification, whereas in healthy cells those antigens are much more complex. A selection of lambodies that could bind to aGPA, a human erythrocyte membrane glycoprotein that is covered with 16 TF$\alpha$ moieties, through magnetic-activated cell sorting (MACS) and fluorescence-activated cell sorting (FACS) has yielded a leucine-rich lambody VLRB.aGPA.23. This lambody selectively stained (over healthy samples) cells from 14 different types of adenocarcinomas: bladder, esophagus, ovary, tongue, cheek, cervix, liver, nose, nasopharynx, greater omentum, colon, breast, larynx, and lung. Moreover, patients whose tissues stained positive with VLRB.aGPA.23 had a significantly smaller survival rate.
A close look at the crystal structure of VLRB.aGPA.23 reveals a tryptophan residue at position 187 right over the carbohydrate binding pocket.

### Multivalency in structure

Many glycan binding proteins (GBPs) are oligomeric and typically contain multiple sites for glycan binding (also called carbohydrate-recognition domains). The ability to form multivalent protein-ligand interactions significantly enhances the strength of binding: while $K_{D}$ values for individual CRD-glycan interactions may be in the mM range, the overall affinity of GBP towards glycans may reach nanomolar or even picomolar ranges. The overall strength of interactions is described as avidity $K_{D}$ (in contrast with an affinity $K_{D}$ which describes single equilibrium). Sometimes the avidity is also called an apparent $K_{D}$ to emphasize the non-equilibrium nature of the interaction.
Common oligomerization structures of lectins are shown below. For example, galectins are usually observed as dimers, while intelectins form trimers and pentraxins assemble into pentamers. Larger structures, like hexameric Reg proteins, may assemble into membrane penetrating pores. Collectins may form even more bizarre complexes: bouquets of trimers or even cruciform-like structures (e.g. in SP-D).

## Current Research

Given the importance of glycan-protein interactions, there is an ongoing research dedicated to the a) creation of new tools to detect glycan-protein interactions and b) using those tools to decipher the so-called sugar code.

### Glycan Arrays

One of the most widely used tools for probing glycan-protein interactions is glycan arrays. A glycan array usually is an NHS- or epoxy-activated glass slides on which various glycans were printed using robotic printing. These commercially available arrays may contain up to 600 different glycans, specificity of which has been extensively studied.
Glycan-protein interactions may be detected by testing proteins of interest (or libraries of those) that bear fluorescent tags. The structure of the glycan-binding protein may be deciphered by several analytical methods based on mass-spectrometry, including MALDI-MS, LC-MS, tandem MS-MS, and/or 2D NMR.

### Bioinformatics driven research

Computational methods have been applied to search for parameters (e.g. residue propensity, hydrophobicity, planarity) that could distinguish glycan-binding proteins from other surface patches. For example, a model trained on 19 non-homologous carbohydrate binding structures was able to predict carbohydrate-binding domains (CRDs) with an accuracy of 65% for non-enzymatic structures and 87% for enzymatic ones. Further studies have employed calculations of Van der Waals energies of protein-probe interactions and amino acid propensities to identify CRDs with 98% specificity at 73% sensitivity. More recent methods can predict CRDs even from protein sequences, by comparing the sequence with those for which structures are already known.

### Sugar code

In contrast with protein studies, where a primary protein structure is unambiguously defined by the sequence of nucleotides (the genetic code), the glycobiology still cannot explain how a certain "message" is encoded using carbohydrates or how it is "read" and "translated" by other biological entities.
An interdisciplinary effort, combining chemistry, biology, and biochemistry, studies glycan-protein interactions to see how different sequences of carbohydrates initiate different cellular responses.