Computational toolbox for the analysis of protein-glycan interactions

Protein-glycan interactions play pivotal roles in numerous biological processes, ranging from cellular recognition to immune response modulation. Understanding the intricate details of these interactions is crucial for deciphering the molecular mechanisms underlying various physiological and pathological conditions. Computational techniques have emerged as powerful tools that can help in drawing, building and visualising complex biomolecules and provide insights into their dynamic behaviour at atomic and molecular levels. This review provides an overview of the main computational tools useful for studying biomolecular systems, particularly glycans, both in free state and in complex with proteins, also with reference to the principles, methodologies, and applications of all-atom molecular dynamics simulations. Herein, we focused on the programs that are generally employed for preparing protein and glycan input files to execute molecular dynamics simulations and analyse the corresponding results. The presented computational toolbox represents a valuable resource for researchers studying protein-glycan interactions and incorporates advanced computational methods for building, visualising and predicting protein/glycan structures, modelling protein-ligand complexes, and analyse MD outcomes. Moreover, selected case studies have been reported to highlight the importance of computational tools in studying protein-glycan systems, revealing the capability of these tools to provide valuable insights into the binding kinetics, energetics, and structural determinants that govern specific molecular interactions.

Systematic Preparation of a 66-IgG Library with Symmetric and Asymmetric Homogeneous Glycans and Their Functional Evaluation

Immunoglobulin G (IgG) antibodies possess a conserved N-glycosylation site in the Fc domain. In FcγRIIIa affinity column chromatography, unglycosylated, hemiglycosylated, and fully glycosylated IgG retention times differ considerably. Using retention-time differences, 66 different trastuzumab antibodies with symmetric and asymmetric homogeneous glycans were prepared systematically, substantially expanding the scope of IgGs with homogeneous glycans. Using the prepared trastuzumab with homogeneous glycans, thermal stability and antibody-dependent cellular cytotoxicity were investigated. In some glycan series, a directly proportional relationship was observed between the thermal unfolding temperature (Tm) and the calorimetric unfolding heat (ΔHcal). Antibody function could be deduced from the combination of a pair of glycans in an intact form. Controlling glycan structure through the combination of a pair of glycans permits the precise tuning of stability and effector functions of IgG. Overall, our technology can be used to investigate the effects of glycans on antibody functions.

Interaction of Asn297-Linked Glycan Ligands with the Fc Fragment of the Immunoglobulin Class G1: A Molecular Dynamics Simulation Study

The allosteric modulation of the homodimeric H10-03-6 protein to glycan ligands L1 and L2, and the STAB19 protein to glycan ligands L3 and L4, respectively, has been studied by molecular dynamics simulations and free energy calculations. The results revealed that the STAB19 protein has a significantly higher affinity for L3 (-11.38 ± 2.32 kcal/mol) than that for L4 (-5.51 ± 1.92 kcal/mol). However, the combination of the H10-03-6 protein with glycan L2 (1.23 ± 6.19 kcal/mol) is energetically unfavorable compared with that of L1 (-13.96 ± 0.35 kcal/mol). Further, the binding of glycan ligands L3 and L4 to STAB19 would result in the significant closure of the two CH2 domains of the STAB19 conformation with the decrease of the centroid distances between the two CH2 domains compared with the H10-03-6/L1/L2 complex. The CH2 domain closure of STAB19 relates directly to the formation of new hydrogen bonds and hydrophobic interactions between the residues Ser239, Val240, Asp265, Glu293, Asn297, Thr299, Ser337, Asp376, Thr393, Pro395, and Pro396 in STAB19 and glycan ligands L3 and L4, which suggests that these key residues would contribute to the specific regulation of STAB19 to L3 and L4. In addition, the distance analysis revealed that the EF loop in the H10-03-6/L1/L2 model presents a high flexibility and partial disorder compared with the stabilized STAB19/L3/L4 complex. These results will be helpful in understanding the specific regulation through the asymmetric structural characteristics in the CH2 and CH3 domains of the H10-03-6 and STAB19 proteins.

Impact of BA.1, BA.2, and BA.4/BA.5 Omicron mutations on therapeutic monoclonal antibodies

The emergence of Omicron SARS-CoV-2 subvariants (BA.1, BA.2, BA.4, and BA.5), with an unprecedented number of mutations in their receptor-binding domain (RBD) of the spike-protein, has fueled a resurgence of COVID-19 infections, posing a major challenge to the efficacy of existing vaccines and monoclonal antibody (mAb) therapeutics. We conducted a systematic molecular dynamics (MD) simulation to investigate how the RBD mutations of these subvariants affect the interactions with broad mAbs including AstraZeneca (COV2-2196 and COV2-2130), Brii Biosciences (BRII-196), Celltrion (CT-P59), Eli Lilly (LY-CoV555 and LY-CoV016), Regeneron (REGN10933 and REGN10987), Vir Biotechnology (S309), and S2X259. Our results show a complete loss of binding for COV2-2196, BRII-196, CT-P59, and LY-CoV555 with all Omicron RBDs. Additionally, REGN10987 totally loses its binding with BA.1, but retains a partial binding with BA.2 and BA.4/5. The binding reduction is significant for LY-CoV016 and REGN10933 but moderate for COV2-2130. S309 and S2X259 retain their binding with BA.1 but exhibit decreased binding with other subvariants. We introduce a mutational escape map for each mAb to identify the key RBD sites and the corresponding critical mutations. Overall, our findings suggest that the majority of therapeutic mAbs have diminished or missing activity against Omicron subvariants, indicating the urgent need for a new therapeutic mAb with a better design.

Effect of Fc core fucosylation and light chain isotype on IgG1 flexibility

N-glycosylation plays a key role in modulating the bioactivity of monoclonal antibodies (mAbs), as well as the light chain (LC) isotype can influence their physicochemical properties. However, investigating the impact of such features on mAbs conformational behavior is a big challenge, due to the very high flexibility of these biomolecules. In this work we investigate, by accelerated molecular dynamics (aMD), the conformational behavior of two commercial immunoglobulins G1 (IgG1), representative of κ and λ LCs antibodies, in both their fucosylated and afucosylated forms. Our results show, through the identification of a stable conformation, how the combination of fucosylation and LC isotype modulates the hinge behavior, the Fc conformation and the position of the glycan chains, all factors potentially affecting the binding to the FcγRs. This work also represents a technological enhancement in the conformational exploration of mAbs, making aMD a suitable approach to clarify experimental results.

Extensive substrate recognition by the streptococcal antibody-degrading enzymes IdeS and EndoS

Enzymatic cleavage of IgG antibodies is a common strategy used by pathogenic bacteria to ablate immune effector function. The Streptococcus pyogenes bacterium secretes the protease IdeS and the glycosidase EndoS, which specifically catalyse cleavage and deglycosylation of human IgG, respectively. IdeS has received clinical approval for kidney transplantation in hypersensitised individuals, while EndoS has found application in engineering antibody glycosylation. We present crystal structures of both enzymes in complex with their IgG1 Fc substrate, which was achieved using Fc engineering to disfavour preferential Fc crystallisation. The IdeS protease displays extensive Fc recognition and encases the antibody hinge. Conversely, the glycan hydrolase domain in EndoS traps the Fc glycan in a "flipped-out" conformation, while additional recognition of the Fc peptide is driven by the so-called carbohydrate binding module. In this work, we reveal the molecular basis of antibody recognition by bacterial enzymes, providing a template for the development of next-generation enzymes.

Immobilization of Biantennary N-Glycans Leads to Branch Specific Epitope Recognition by LSECtin

The molecular recognition features of LSECtin toward asymmetric N-glycans have been scrutinized by NMR and compared to those occurring in glycan microarrays. A pair of positional glycan isomers (LDN3 and LDN6), a nonelongated GlcNAc4Man3 N-glycan (G0), and the minimum binding epitope (the GlcNAcβ1-2Man disaccharide) have been used to shed light on the preferred binding modes under both experimental conditions. Strikingly, both asymmetric LDN3 and LDN6 N-glycans are recognized by LSECtin with similar affinities in solution, in sharp contrast to the results obtained when those glycans are presented on microarrays, where only LDN6 was efficiently recognized by the lectin. Thus, different results can be obtained using different experimental approaches, pointing out the tremendous difficulty of translating in vitro results to the in vivo environment.

Exploration, Representation, and Rationalization of the Conformational Phase Space of N-Glycans

Despite their fundamental biological relevance, structure-property relationships in N-glycans are fundamentally lacking, and their highly multidimensional compositional and conformational phase spaces remain largely unexplored. The torsional flexibility of the glycosidic linkages and the ring dynamics result in wide, rugged free-energy landscapes that are difficult to sample in molecular dynamics simulations. We show that a novel enhanced-sampling scheme combining replica exchange with solute and collective-variable tempering, enabling transitions over all relevant energy barriers, delivers converged distributions of solvated N-glycan conformers. Several dimensionality-reduction algorithms are compared and employed to generate conformational free-energy maps in two dimensions. Together with an originally developed conformation-based nomenclature scheme that uniquely identifies glycan conformers, our modeling procedure is applied to reveal the effect of chemical substitutions on the conformational ensemble of selected high-mannose-type and complex glycans. Moreover, the structure-prediction capabilities of two commonly used glycan force fields are assessed via the theoretical prediction of experimentally available nuclear magnetic resonance J-coupling constants. The results especially confirm the key role of ω and ψ torsion angles in discriminating between different conformational states and suggest an intriguing correlation between the torsional and ring-puckering degrees of freedom that may be biologically relevant.

Molecular simulations of complex carbohydrates and glycoconjugates

Complex carbohydrates (glycans) are the most abundant and versatile biopolymers in nature. The broad diversity of biochemical functions that carbohydrates cover is a direct consequence of the variety of 3D architectures they can adopt, displaying branched or linear arrangements, widely ranging in sizes, and with the highest diversity of building blocks of any other natural biopolymer. Despite this unparalleled complexity, a common denominator can be found in the glycans' inherent flexibility, which hinders experimental characterization, but that can be addressed by high-performance computing (HPC)-based molecular simulations. In this short review, I present and discuss the state-of-the-art of molecular simulations of complex carbohydrates and glycoconjugates, highlighting methodological strengths and weaknesses, important insights through emblematic case studies, and suggesting perspectives for future developments.

Successive remodeling of IgG glycans using a solid-phase enzymatic platform

The success of glycoprotein-based drugs in various disease treatments has become widespread. Frequently, therapeutic glycoproteins exhibit a heterogeneous array of glycans that are intended to mimic human glycopatterns. While immunogenic responses to biologic drugs are uncommon, enabling exquisite control of glycosylation with minimized microheterogeneity would improve their safety, efficacy and bioavailability. Therefore, close attention has been drawn to the development of glycoengineering strategies to control the glycan structures. With the accumulation of knowledge about the glycan biosynthesis enzymes, enzymatic glycan remodeling provides a potential strategy to construct highly ordered glycans with improved efficiency and biocompatibility. In this study, we quantitatively evaluate more than 30 enzymes for glycoengineering immobilized immunoglobulin G, an impactful glycoprotein class in the pharmaceutical field. We demonstrate successive glycan remodeling in a solid-phase platform, which enabled IgG glycan harmonization into a series of complex-type N-glycoforms with high yield and efficiency while retaining native IgG binding affinity.

A solid-phase glycan remodeling (SPGR) platform is presented. Over thirty enzymes were evaluated for successive glycoengineering of immobilized antibodies with outstanding performance in several SPGR reactions.

Conformational preferences of triantennary and tetraantennary hybrid N-glycans in aqueous solution: Insights from 20 μs long atomistic molecular dynamic simulations

In the current study, we have investigated the conformational dynamics of a triantennary (N-glycan1) and tetraantennary (N-glycan2) hybrid N-glycans found on the surface of the HIV glycoprotein using 20 μs long all-atom molecular dynamics (MD) simulations. The main objective of the present study is to elucidate the influence of adding a complex branch on the overall glycan structural dynamics. Our investigation suggests that the average RMSD value increases when a complex branch is added to N-glycan1. However, the RMSD distribution is relatively wider in the case of N-glycan1 compared to N-glycan2, which indicates that multiple complex branches restrict the conformational variability of glycans. A similar observation is obtained from the principal component analysis of both glycans. All the puckering states (⁴C₁ to ¹C₄) of each monosaccharide except mannose are sampled in our simulations, although the ⁴C₁ chair form is energetically more favorable than ¹C₄. In N-glycan1, the 1-6 linkage in the mannose branch [Man(9)-α(1-6)-Man(5)] stays in the gauche-gauche cluster, whereas it moves towards trans-gauche in N-glycan2. For both glycans, mannose branches are more flexible than the complex branches, and adding a complex branch does not influence the dynamics of the mannose branches. We have noticed that the end-to-end distance of the complex branch shortens by ∼ 10 Å in the presence of another complex branch. This suggests that in the presence of an additional complex branch, the other complex branch adopts a close folded structure. All these conformational changes involve the selective formation of inter-residue and water-mediated hydrogen-bond networks.

Fine-tuning the spike: role of the nature and topology of the glycan shield in the structure and dynamics of the SARS-CoV-2 S

The dense glycan shield is an essential feature of the SARS-CoV-2 spike (S) architecture, key to immune evasion and to the activation of the prefusion conformation. Recent studies indicate that the occupancy and structures of the SARS-CoV-2 S glycans depend not only on the nature of the host cell, but also on the structural stability of the trimer; a point that raises important questions about the relative competence of different glycoforms. Moreover, the functional role of the glycan shield in the SARS-CoV-2 pathogenesis suggests that the evolution of the sites of glycosylation is potentially intertwined with the evolution of the protein sequence to affect optimal activity. Our results from multi-microsecond molecular dynamics simulations indicate that the type of glycosylation at N234, N165 and N343 greatly affects the stability of the receptor binding domain (RBD) open conformation, and thus its exposure and accessibility. Furthermore, our results suggest that the loss of glycosylation at N370, a newly acquired modification in the SARS-CoV-2 S glycan shield's topology, may have contributed to increase the SARS-CoV-2 infectivity as we find that N-glycosylation at N370 stabilizes the closed RBD conformation by binding a specific cleft on the RBD surface. We discuss how the absence of the N370 glycan in the SARS-CoV-2 S frees the RBD glycan binding cleft, which becomes available to bind cell-surface glycans, and potentially increases host cell surface localization.

IgG1 conformational behavior: elucidation of the N-glycosylation role via molecular dynamics

Currently, monoclonal antibodies (mAbs) are the most used biopharmaceuticals for human therapy. One of the key aspects in their development is the control of effector functions mediated by the interaction between fragment crystallizable (Fc) and Fcγ receptors, which is a secondary mechanism of the action of biotherapeutics. N-glycosylation at the Fc portion can regulate these mechanisms, and much experimental evidence suggests that modifications of glycosidic chains can affect antibody binding to FcγRIIIa, consequently impacting the immune response. In this work, we try to elucidate via in silico procedures the structural role exhibited by glycans, particularly fucose, in mAb conformational freedom that can potentially affect the receptor recognition. By using adalimumab, a marketed IgG1, as a general template, after rebuilding its three-dimensional (3D) structure through homology modeling approaches, we carried out molecular dynamics simulations of three differently glycosylated species: aglycosylated, afucosylated, and fucosylated antibody. Trajectory analysis showed different dynamical behaviors and pointed out that sugars can influence the overall 3D structure of the antibody. As a result, we propose a putative structural mechanism by which the presence of fucose introduces conformational constraints in the whole antibody and not only in the Fc domain, preventing a conformation suitable for the interaction with the receptor. As secondary evidence, we observed a high flexibility of the antibodies that is translated into an asymmetric behavior of Fab portions shown by all the simulated biopolymers, making the dynamical asymmetry a new, to our knowledge, molecular aspect that may be further investigated. In conclusion, these findings can help understand the contribution of sugars on the structural architecture of mAbs, paving the way to novel strategies of pharmaceutical development.

The Role of Fucose-Containing Glycan Motifs Across Taxonomic Kingdoms

The extraordinary diversity of glycans leads to large differences in the glycomes of different kingdoms of life. Yet, while most monosaccharides are solely found in certain taxonomic groups, there is a small set of monosaccharides with widespread distribution across nearly all domains of life. These general monosaccharides are particularly relevant for glycan motifs, as they can readily be used by commensals and pathogens to mimic host glycans or hijack existing glycan recognition systems. Among these, the monosaccharide fucose is especially interesting, as it frequently presents itself as a terminal monosaccharide, primed for interaction with proteins. Here, we analyze fucose-containing glycan motifs across all taxonomic kingdoms. Using a hereby presented large species-specific glycan dataset and a plethora of methods for glycan-focused bioinformatics and machine learning, we identify characteristic as well as shared fucose-containing glycan motifs for various taxonomic groups, demonstrating clear differences in fucose usage. Even within domains, fucose is used differentially based on an organism’s physiology and habitat. We particularly highlight differences in fucose-containing motifs between vertebrates and invertebrates. With the example of pathogenic and non-pathogenic Escherichia coli strains, we also demonstrate the importance of fucose-containing motifs in molecular mimicry and thereby pathogenic potential. We envision that this study will shed light on an important class of glycan motifs, with potential new insights into the role of fucosylated glycans in symbiosis, pathogenicity, and immunity.

Novel diagnostic options for endometriosis - Based on the glycome and microbiome

Background

Endometriosis is a chronic gynaecological disease whose aetiology is still unknown. Despite its prevalence among women of reproductive age, the pathology of the disease has not yet been elucidated and only symptomatic treatment is available. Endometriosis has high latency and diagnostic methods are both limited and invasive.

Aim of review

The aim of this review is to summarise minimally invasive or non-invasive diagnostic methods for endometriosis and their diagnostic efficiencies. Furthermore, we discuss the identification and diagnostic potential of novel disease biomarkers of microbial or glycan origin.

Key scientific concepts of review

Great efforts have been made to develop minimally invasive or non-invasive diagnostic methods in endometriosis. The problem with most potential biomarker candidates is that they have high accuracy only in cases of severe disease. Therefore, it is necessary to examine other potential biomarkers more closely. Associations between gastrointestinal and genital tract microbial health and endometriosis have been identified. For instance, irritable bowel syndrome is more common in women with endometriosis, and hormonal imbalance has a negative impact on the microbiome of both the genital tract and the gastrointestinal system. Further interrogation of these associations may have potential diagnostic significance and may identify novel therapeutic avenues. Glycomics may also be a potent source of biomarkers of endometriosis, with a number of glyco-biomarkers already approved by the FDA. Endometriosis-associated microbial and glycomic profiles may represent viable targets for development of innovative diagnostics in this debilitating disease.

Fc Galactosylation Promotes Hexamerization of Human IgG1, Leading to Enhanced Classical Complement Activation

Human IgG contains one evolutionarily conserved N-linked glycan in its Fc region at position 297. This glycan is crucial for Fc-mediated functions, including its induction of the classical complement cascade. This is induced after target recognition through the IgG-Fab regions, allowing neighboring IgG-Fc tails to associate through Fc:Fc interaction, ultimately leading to hexamer formation. This hexamerization seems crucial for IgG to enable efficient interaction with the globular heads of the first complement component C1q and subsequent complement activation. In this study, we show that galactose incorporated in the IgG1-Fc enhances C1q binding, C4, C3 deposition, and complement-dependent cellular cytotoxicity in human erythrocytes and Raji cells. IgG1-Fc sialylation slightly enhanced binding of C1q, but had little effect on downstream complement activation. Using various mutations that decrease or increase hexamerization capacity of IgG1, we show that IgG1-Fc galactosylation has no intrinsic effect on C1q binding to IgG1, but enhances IgG1 hexamerization potential and, thereby, complement activation. These data suggest that the therapeutic potential of Abs can be amplified without introducing immunogenic mutations, by relatively simple glycoengineering.

SARS-CoV-2 RBD antibodies that maximize breadth and resistance to escape

An ideal therapeutic anti-SARS-CoV-2 antibody would resist viral escape1-3, have activity against diverse sarbecoviruses4-7, and be highly protective through viral neutralization8-11 and effector functions12,13. Understanding how these properties relate to each other and vary across epitopes would aid the development of therapeutic antibodies and guide vaccine design. Here we comprehensively characterize escape, breadth and potency across a panel of SARS-CoV-2 antibodies targeting the receptor-binding domain (RBD). Despite a trade-off between in vitro neutralization potency and breadth of sarbecovirus binding, we identify neutralizing antibodies with exceptional sarbecovirus breadth and a corresponding resistance to SARS-CoV-2 escape. One of these antibodies, S2H97, binds with high affinity across all sarbecovirus clades to a cryptic epitope and prophylactically protects hamsters from viral challenge. Antibodies that target the angiotensin-converting enzyme 2 (ACE2) receptor-binding motif (RBM) typically have poor breadth and are readily escaped by mutations despite high neutralization potency. Nevertheless, we also characterize a potent RBM antibody (S2E128) with breadth across sarbecoviruses related to SARS-CoV-2 and a high barrier to viral escape. These data highlight principles underlying variation in escape, breadth and potency among antibodies that target the RBD, and identify epitopes and features to prioritize for therapeutic development against the current and potential future pandemics.

An atomistic perspective on antibody-dependent cellular cytotoxicity quenching by core-fucosylation of IgG1 Fc N-glycans from enhanced sampling molecular dynamics

The immunoglobulin type G (IgG) Fc N-glycans are known to modulate the interaction with membrane-bound Fc γ receptors (FcγRs), fine-tuning the antibody's effector function in a sequence-dependent manner. Particularly interesting in this respect are the roles of galactosylation, which levels are linked to autoimmune conditions and aging, of core fucosylation, which is known to reduce significantly the antibody-dependent cellular cytotoxicity (ADCC), and of sialylation, which also reduces antibody-dependent cellular cytotoxicity (ADCC) but only in the context of core-fucosylation. In this article, we provide an atomistic level perspective through enhanced sampling computer simulations, based on replica exchange molecular dynamics (REMD), to understand the molecular determinants linking the Fc N-glycans sequence to the observed IgG1 function. Our results indicate that the two symmetrically opposed N-glycans interact extensively through their core trimannose residues. At room temperature, the terminal galactose on the α (1-6) arm is restrained to the protein through a network of interactions that keep the arm outstretched; meanwhile, the α (1-3) arm extends toward the solvent where a terminal sialic acid remains fully accessible. We also find that the presence of core fucose interferes with the extended sialylated α (1-3) arm, altering its conformational propensity and as a consequence of steric hindrance, significantly enhancing the Fc dynamics. Furthermore, structural analysis shows that the core-fucose position within the Fc core obstructs the access of N162 glycosylated FcγRs very much like a "door-stop," potentially decreasing the IgG/FcγR binding free energy. These results provide an atomistic level-of-detail framework for the design of high potency IgG1 Fc N-glycoforms.

Antibodies to the SARS-CoV-2 receptor-binding domain that maximize breadth and resistance to viral escape

An ideal anti-SARS-CoV-2 antibody would resist viral escape 1-3 , have activity against diverse SARS-related coronaviruses 4-7 , and be highly protective through viral neutralization 8-11 and effector functions 12,13 . Understanding how these properties relate to each other and vary across epitopes would aid development of antibody therapeutics and guide vaccine design. Here, we comprehensively characterize escape, breadth, and potency across a panel of SARS-CoV-2 antibodies targeting the receptor-binding domain (RBD), including S309 4 , the parental antibody of the late-stage clinical antibody VIR-7831. We observe a tradeoff between SARS-CoV-2 in vitro neutralization potency and breadth of binding across SARS-related coronaviruses. Nevertheless, we identify several neutralizing antibodies with exceptional breadth and resistance to escape, including a new antibody (S2H97) that binds with high affinity to all SARS-related coronavirus clades via a unique RBD epitope centered on residue E516. S2H97 and other escape-resistant antibodies have high binding affinity and target functionally constrained RBD residues. We find that antibodies targeting the ACE2 receptor binding motif (RBM) typically have poor breadth and are readily escaped by mutations despite high neutralization potency, but we identify one potent RBM antibody (S2E12) with breadth across sarbecoviruses closely related to SARS-CoV-2 and with a high barrier to viral escape. These data highlight functional diversity among antibodies targeting the RBD and identify epitopes and features to prioritize for antibody and vaccine development against the current and potential future pandemics.

Oligomannose N-Glycans 3D Architecture and Its Response to the FcγRIIIa Structural Landscape

Oligomannoses are evolutionarily the oldest class of N-glycans, where the arms of the common pentasaccharide unit, i.e., Manα(1-6)-[Manα(1-3)]-Manβ(1-4)-GlcNAcβ(1-4)-GlcNAcβ1-Asn, are functionalized exclusively with branched arrangements of mannose (Man) monosaccharide units. In mammalian species oligomannose N-glycans can have up to 9 Man; meanwhile structures can grow to over 200 units in yeast mannan. The highly dynamic nature, branching complexity, and 3D structure of oligomannoses have been recently highlighted for their roles in immune escape and infectivity of enveloped viruses, such as HIV-1 and SARS-CoV2. The architectural features that allow these N-glycans to perform their functions are yet unclear, due to their intrinsically disordered nature that hinders their structural characterization. In this work we will discuss the results of over 54 μs of cumulative sampling by molecular dynamics (MD) simulations of differently processed, free (not protein-linked) oligomannose N-glycans common in vertebrates. We then discuss the effects of a protein surface on their structural equilibria based on over 4 μs cumulative MD sampling of the fully glycosylated CD16a Fc γ receptor (FcγRIIIa), where the type of glycosylation is known to modulate its binding affinity for IgG1s, regulating the antibody-dependent cellular cytotoxicity (ADCC). Our results show that the protein's structural constraints shift the oligomannoses conformational ensemble to promote conformers that satisfy the steric requirements and hydrogen bonding networks demanded by the protein's surface landscape. More importantly, we find that the protein does not actively distort the N-glycans into structures not populated in the unlinked forms in solution. Ultimately, the highly populated conformations of the Man5 linked glycans support experimental evidence of high levels of hybrid complex forms at N45 and show a specific presentation of the arms at N162, which may be involved in mediating binding affinity to the IgG1 Fc.

Circulating SARS-CoV-2 spike N439K variants maintain fitness while evading antibody-mediated immunity

SARS-CoV-2 can mutate and evade immunity, with consequences for efficacy of emerging vaccines and antibody therapeutics. Here, we demonstrate that the immunodominant SARS-CoV-2 spike (S) receptor binding motif (RBM) is a highly variable region of S and provide epidemiological, clinical, and molecular characterization of a prevalent, sentinel RBM mutation, N439K. We demonstrate N439K S protein has enhanced binding affinity to the hACE2 receptor, and N439K viruses have similar in vitro replication fitness and cause infections with similar clinical outcomes as compared to wild type. We show the N439K mutation confers resistance against several neutralizing monoclonal antibodies, including one authorized for emergency use by the US Food and Drug Administration (FDA), and reduces the activity of some polyclonal sera from persons recovered from infection. Immune evasion mutations that maintain virulence and fitness such as N439K can emerge within SARS-CoV-2 S, highlighting the need for ongoing molecular surveillance to guide development and usage of vaccines and therapeutics.

Multiscale Simulations Examining Glycan Shield Effects on Drug Binding to Influenza Neuraminidase

Influenza neuraminidase is an important drug target. Glycans are present on neuraminidase and are generally considered to inhibit antibody binding via their glycan shield. In this work, we studied the effect of glycans on the binding kinetics of antiviral drugs to the influenza neuraminidase. We created all-atom in silico systems of influenza neuraminidase with experimentally derived glycoprofiles consisting of four systems with different glycan conformations and one system without glycans. Using Brownian dynamics simulations, we observe a two- to eightfold decrease in the rate of ligand binding to the primary binding site of neuraminidase due to the presence of glycans. These glycans are capable of covering much of the surface area of neuraminidase, and the ligand binding inhibition is derived from glycans sterically occluding the primary binding site on a neighboring monomer. Our work also indicates that drugs preferentially bind to the primary binding site (i.e., the active site) over the secondary binding site, and we propose a binding mechanism illustrating this. These results help illuminate the complex interplay between glycans and ligand binding on the influenza membrane protein neuraminidase.

Beyond Shielding: The Roles of Glycans in the SARS-CoV-2 Spike Protein

The ongoing COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in more than 28,000,000 infections and 900,000 deaths worldwide to date. Antibody development efforts mainly revolve around the extensively glycosylated SARS-CoV-2 spike (S) protein, which mediates host cell entry by binding to the angiotensin-converting enzyme 2 (ACE2). Similar to many other viral fusion proteins, the SARS-CoV-2 spike utilizes a glycan shield to thwart the host immune response. Here, we built a full-length model of the glycosylated SARS-CoV-2 S protein, both in the open and closed states, augmenting the available structural and biological data. Multiple microsecond-long, all-atom molecular dynamics simulations were used to provide an atomistic perspective on the roles of glycans and on the protein structure and dynamics. We reveal an essential structural role of N-glycans at sites N165 and N234 in modulating the conformational dynamics of the spike's receptor binding domain (RBD), which is responsible for ACE2 recognition. This finding is corroborated by biolayer interferometry experiments, which show that deletion of these glycans through N165A and N234A mutations significantly reduces binding to ACE2 as a result of the RBD conformational shift toward the "down" state. Additionally, end-to-end accessibility analyses outline a complete overview of the vulnerabilities of the glycan shield of the SARS-CoV-2 S protein, which may be exploited in the therapeutic efforts targeting this molecular machine. Overall, this work presents hitherto unseen functional and structural insights into the SARS-CoV-2 S protein and its glycan coat, providing a strategy to control the conformational plasticity of the RBD that could be harnessed for vaccine development.

Beyond Shielding: The Roles of Glycans in SARS-CoV-2 Spike Protein

The ongoing COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in more than 15,000,000 infections and 600,000 deaths worldwide to date. Antibody development efforts mainly revolve around the extensively glycosylated SARS-CoV-2 spike (S) protein, which mediates the host cell entry by binding to the angiotensin-converting enzyme 2 (ACE2). Similar to many other viruses, the SARS-CoV-2 spike utilizes a glycan shield to thwart the host immune response. Here, we built a full-length model of glycosylated SARS-CoV-2 S protein, both in the open and closed states, augmenting the available structural and biological data. Multiple microsecond-long, all-atom molecular dynamics simulations were used to provide an atomistic perspective on the roles of glycans, and the protein structure and dynamics. We reveal an essential structural role of N-glycans at sites N165 and N234 in modulating the conformational dynamics of the spike's receptor binding domain (RBD), which is responsible for ACE2 recognition. This finding is corroborated by biolayer interferometry experiments, which show that deletion of these glycans through N165A and N234A mutations significantly reduces binding to ACE2 as a result of the RBD conformational shift towards the "down" state. Additionally, end-to-end accessibility analyses outline a complete overview of the vulnerabilities of the glycan shield of SARS-CoV-2 S protein, which may be exploited by therapeutic efforts targeting this molecular machine. Overall, this work presents hitherto unseen functional and structural insights into the SARS-CoV-2 S protein and its glycan coat, providing a strategy to control the conformational plasticity of the RBD that could be harnessed for vaccine development.

How and why plants and human N-glycans are different: Insight from molecular dynamics into the "glycoblocks" architecture of complex carbohydrates

The N-glycosylation is one of the most abundant and diverse post-translational modifications of proteins, implicated in protein folding and structural stability, and mediating interactions with receptors and with the environment. All N-glycans share a common core from which linear or branched arms stem from, with functionalization specific to different species and to the cells’ health and disease state. This diversity generates a rich collection of structures, all diversely able to trigger molecular cascades and to activate pathways, which also include adverse immunogenic responses. These events are inherently linked to the N-glycans’ 3D architecture and dynamics, which remain for the large part unresolved and undetected because of their intrinsic structural disorder. In this work we use molecular dynamics (MD) simulations to provide insight into N-glycans’ 3D structure by analysing the effects of a set of very specific modifications found in plants and invertebrate N-glycans, which are immunogenic in humans. We also compare these structural motifs and combine them with mammalian N-glycan motifs to devise strategies for the control of the N-glycan 3D structure through sequence. Our results suggest that the N-glycans’ architecture can be described in terms of the local spatial environment of groups of monosaccharides. We define these “glycoblocks” as self-contained 3D units, uniquely identified by the nature of the residues they comprise, their linkages and structural/dynamic features. This alternative description of glycans’ 3D architecture can potentially lead to an easier prediction of sequence-to-structure relationships in complex carbohydrates, with important implications in glycoengineering design.

Effects of sialic acid linkage on antibody-fragment crystallizable receptor binding and antibody dependent cytotoxicity depend on levels of fucosylation/bisecting

Aim: Fragment crystallizable (Fc) glycosylation of immunoglobulin G-type monoclonal antibodies applied to therapeutic applications is regarded a critical quality attribute and can influence bioactivity, pharmacokinetics and/or immunogenicity/safety. Investigating the impact of certain Fc N-glycans is therefore of importance to assess its criticality for a therapeutic product. This has been done for N-glycan types like fucosylation, galactosylation or sialylation. There were contradictory results reported for functionality especially with regard to sialylation. Material & methods: We elucidated the effect of terminal sialic acid residues on Fcγ receptor binding and antibody dependent cytotoxicity activity of two immunoglobulin G1 antibodies with different levels of fucosylation/bi-secting. Conclusion: We found the impact to be specific to the sialylation linkage type, in other words, α2,3- versus α2,6-linked sialic acid attached to the terminal galactose residues.

Identification of the epitope of 10-7G glycan antibody to recognize cancer-associated haptoglobin

We previously identified fucosylated haptoglobin (Fuc-Hpt) as a clinical serum biomarker of pancreatic cancer and established the novel glycan monoclonal antibody (mAb) 10-7G. This antibody recognizes cancer-associated haptoglobin including Fuc-Hpt and the precursor of haptoglobin. Interestingly, Western blot analysis showed that the 10-7G mAb reacts with the haptoglobin α chain, which has no N-glycan potential sites; haptoglobin β chain has four N-glycan sites. In this study, we identified the epitope for the 10-7G mAb using haptoglobin deletion mutants, as well as inhibition ELISA with recombinant peptides. We illustrated molecular graphics to show a relationship between the epitope and the β chain. Furthermore, we hypothesized that the 10-7G mAb minimally recognizes normal haptoglobin, but aberrant glycosylation on the β chain causes conformational changes, enabling the 10-7G mAb to easily access the epitope within the α chain. Because 10-7G values, an enzyme-linked immunosorbent assay-immobilized 10-7G mAb, in patients with pancreatic cancer varied by haptoglobin phenotype, the amount of aberrant glycosylation needed to affect haptoglobin conformation probably depends on haptoglobin phenotype. Taken together, the 10-7G mAb recognized characteristic peptides on the haptoglobin α chain as a result of conformational changes and is a promising tool for diagnosing pancreatic cancer by haptoglobin phenotype.

Structural basis for the specific cleavage of core-fucosylated N-glycans by endo-β-N-acetylglucosaminidase from the fungus Cordyceps militaris

N-Linked glycans play important roles in various cellular and immunological events. Endo-β-N-acetylglucosaminidase (ENGase) can release or transglycosylate N-glycans and is a promising tool for the chemoenzymatic synthesis of glycoproteins with homogeneously modified glycans. The ability of ENGases to act on core-fucosylated glycans is a key factor determining their therapeutic utility because mammalian N-glycans are frequently α-1,6-fucosylated. Although the biochemistries and structures of various ENGases have been studied extensively, the structural basis for the recognition of the core fucose and the asparagine-linked GlcNAc is unclear. Herein, we determined the crystal structures of a core fucose-specific ENGase from the caterpillar fungus Cordyceps militaris (Endo-CoM), which belongs to glycoside hydrolase family 18. Structures complexed with fucose-containing ligands were determined at 1.75-2.35 Å resolutions. The fucose moiety linked to GlcNAc is extensively recognized by protein residues in a round-shaped pocket, whereas the asparagine moiety linked to the GlcNAc is exposed to the solvent. The N-glycan-binding cleft of Endo-CoM is Y-shaped, and several lysine and arginine residues are present at its terminal regions. These structural features were consistent with the activity of Endo-CoM on fucose-containing glycans on rituximab (IgG) and its preference for a sialobiantennary substrate. Comparisons with other ENGases provided structural insights into their core fucose tolerance and specificity. In particular, Endo-F3, a known core fucose-specific ENGase, has a similar fucose-binding pocket, but the surrounding residues are not shared with Endo-CoM. Our study provides a foothold for protein engineering to develop enzymatic tools for the preparation of more effective therapeutic antibodies.

The Impact of Immunoglobulin G1 Fc Sialylation on Backbone Amide H/D Exchange

The usefulness of higher-order structural information provided by hydrogen/deuterium exchange-mass spectrometry (H/DX-MS) for the structural impact analyses of chemical and post-translational antibody modifications has been demonstrated in various studies. However, the structure-function assessment for protein drugs in biopharmaceutical research and development is often impeded by the relatively low-abundance (below 5%) of critical quality attributes or by overlapping effects of modifications, such as glycosylation, with chemical amino acid modifications; e.g., oxidation or deamidation. We present results demonstrating the applicability of the H/DX-MS technique to monitor conformational changes of specific Fc glycosylation variants produced by in vitro glyco-engineering technology. A trend towards less H/DX in Fc Cγ2 domain segments correlating with larger glycan structures could be confirmed. Furthermore, significant deuterium uptake differences and corresponding binding properties to Fc receptors (as monitored by SPR) between α-2,3- and α-2,6-sialylated Fc glycosylation variants were verified at sensitive levels.

Effects of terminal galactose residues in mannose α1-6 arm of Fc-glycan on the effector functions of therapeutic monoclonal antibodies

Typical crystallizable fragment (Fc) glycans attached to the CH2 domain in therapeutic monoclonal antibodies (mAbs) are core-fucosylated and asialo-biantennary complex-type glycans, e.g., G2F (full galactosylation), G1aF (terminal galactosylation on the Man α1-6 arm), G1bF (terminal galactosylation on the Man α1-3 arm), and G0F (non-galactosylation). Terminal galactose (Gal) residues of Fc-glycans are known to influence effector functions such as antibody-dependent cell-mediated cytotoxicity and complement-dependent cytotoxicity (CDC), but the impact of the G1F isomers (G1aF and G1bF) on the effector functions has not been reported. Here, we prepared four types of glycoengineered anti-CD20 mAbs bearing homogeneous G2F, G1aF, G1bF, or G0F (G2F mAb, G1aF mAb, G1bF mAb, or G0F mAb, respectively), and evaluated their biological activities. Interestingly, G1aF mAb showed higher C1q- and FcγR-binding activities, CDC activity, and FcγR-activation property than G1bF mAb. The activities of G1aF mAb and G1bF mAb were at the same level as G2F mAb and G0F mAb, respectively. Hydrogen–deuterium exchange/mass spectrometry analysis of dynamic structures of mAbs revealed the greater involvement of the terminal Gal residue on the Man α1-6 arm in the structural stability of the CH2 domain. Considering that mAbs interact with FcγR and C1q via their hinge proximal region in the CH2 domain, the structural stabilization of the CH2 domain by the terminal Gal residue on the Man α1-6 arm of Fc-glycans may be important for the effector functions of mAbs. To our knowledge, this is the first report showing the impact of G1F isomers on the effector functions and dynamic structure of mAbs.

Abbreviations: ABC, ammonium bicarbonate solution; ACN, acetonitrile; ADCC, antibody-dependent cell-mediated cytotoxicity; C1q, complement component 1q; CDC, complement-dependent cytotoxicity; CQA, critical quality attribute; Endo, endo-β-N-acetylglucosaminidase; FA, formic acid; Fc, crystallizable fragment; FcγR, Fcγ receptors; Fuc, fucose; Gal, galactose; GlcNAc, N-acetylglucosamine; GST, glutathione S-transferase; HER2, human epidermal growth factor receptor 2; HDX, hydrogen–deuterium exchange; HILIC, hydrophilic interaction liquid chromatography; HLB-SPE, hydrophilic-lipophilic balance–solid-phase extraction; HPLC, high-performance liquid chromatography; mAb, monoclonal antibody; Man, mannose; MS, mass spectrometry; PBS, phosphate-buffered saline; SGP, hen egg yolk sialylglycopeptides.

Molecular Basis of Broad Spectrum N-Glycan Specificity and Processing of Therapeutic IgG Monoclonal Antibodies by Endoglycosidase S2

Immunoglobulin G (IgG) glycosylation critically modulates antibody effector functions. Streptococcus pyogenes secretes a unique endo-β-N-acetylglucosaminidase, EndoS2, which deglycosylates the conserved N-linked glycan at Asn297 on IgG Fc to eliminate its effector functions and evade the immune system. EndoS2 and specific point mutants have been used to chemoenzymatically synthesize antibodies with customizable glycosylation for gain of functions. EndoS2 is useful in these schemes because it accommodates a broad range of N-glycans, including high-mannose, complex, and hybrid types; however, its mechanism of substrate recognition is poorly understood. We present crystal structures of EndoS2 alone and bound to complex and high-mannose glycans; the broad N-glycan specificity is governed by critical loops that shape the binding site of EndoS2. Furthermore, hydrolytic experiments, domain-swap chimeras, and hydrogen-deuterium exchange mass spectrometry reveal the importance of the carbohydrate-binding module in the mechanism of IgG recognition by EndoS2, providing insights into engineering enzymes to catalyze customizable glycosylation reactions.