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

The role of antibody glycosylation in autoimmune and alloimmune kidney diseases

Immunoglobulin glycosylation is a pivotal mechanism that drives the diversification of antibody functions. The composition of the IgG glycome is influenced by environmental factors, genetic traits and inflammatory contexts. Differential IgG glycosylation has been shown to intricately modulate IgG effector functions and has a role in the initiation and progression of various diseases. Analysis of IgG glycosylation is therefore a promising tool for predicting disease severity. Several autoimmune and alloimmune disorders, including critical and potentially life-threatening conditions such as systemic lupus erythematosus, anti-neutrophil cytoplasmic antibody (ANCA)-associated vasculitis and antibody-mediated kidney graft rejection, are driven by immunoglobulin. In certain IgG-driven kidney diseases, including primary membranous nephropathy, IgA nephropathy and lupus nephritis, particular glycome characteristics can enhance in situ complement activation and the recruitment of innate immune cells, resulting in more severe kidney damage. Hypofucosylation, hypogalactosylation and hyposialylation are the most common IgG glycosylation traits identified in these diseases. Modulating IgG glycosylation could therefore be a promising therapeutic strategy for regulating the immune mechanisms that underlie IgG-driven kidney diseases and potentially reduce the burden of immunosuppressive drugs in affected patients.

Effects of glycans and hinge on dynamics in the IgG1 Fc

The crystallizable fragment (Fc) domain of immunoglobulin subclass IgG1 antibodies is engineered for a wide variety of pharmaceutical applications. Two important structural variables in Fc constructs are the hinge region connecting the Fc to the antigen binding fragments (Fab) and the glycans present in various glycoforms. These components affect receptor binding interactions that mediate immune activation. To design new antibody drugs, a robust in silico method for linking stability to structural changes is necessary. In this work, all-atom simulations were used to compare the dynamic behavior of the four structural variants arising from presence or absence of the hinge and glycans. We expressed the simplest of these constructs, the 'minimal Fc' with no hinge and no glycans, in Escherichia coli and report its crystal structure. The 'maximal Fc' that includes full hinge and G0F/G1F glycans is based on a previously reported structure, Protein Data Bank (PDB) ID: 5VGP. These, along with two intermediate structures (with only the glycans or with only the hinge) were used to independently measure the stability effects of the two structural variables using umbrella sampling simulations. Principal component analysis (PCA) was used to determine free energy effects along the Fc's dominant mode of motion. This work provides a comprehensive picture of the effects of hinge and glycans on Fc dynamics and stability.Communicated by Ramaswamy H. Sarma.

Serum disease-specific IgG Fc glycosylation as potential biomarkers for nonproliferative diabetic retinopathy using mass spectrometry

OBJECTIVE

To explore the potential of serum disease-specific immunoglobulin G (DSIgG) glycosylation as a biomarker for the diagnosis of nonproliferative diabetic retinopathy (NPDR).

METHODS

A total of 387 consecutive diabetic patients presenting in an eye clinic without proliferative diabetic retinopathy (DR) were included and divided into those with nondiabetic retinopathy (NDR) (n = 181) and NPDR (n = 206) groups. Serum was collected from all patients for DSIgG separation. The enriched glycopeptides of the tryptic digests of DSIgG were detected using matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF MS). Patients were randomly divided into discovery and validation sets (1:1). The differences in glycopeptide ratios between the groups were compared by using Student's t-test or the Mann-Whitney U test. The predictive ability of the model was assessed using the area under the receiver operating characteristic curve (AUC).

RESULTS

DSIgG1 G1FN/G0FN, G2N/G2, G2FN/G2N and DSIgG2 G1F/G0F, G1FN/G0FN, G2N/G1N, G2S/G2 were significantly different between NDR and NPDR patients (p < 0.05) in both the discovery and validation sets. The prediction model that was built comprising the seven glycopeptide ratios showed good NPDR prediction performance with an AUC of 0.85 in the discovery set and 0.87 in the validation set.

CONCLUSION

DSIgG Fc N-glycosylation ratios were associated with NPDR and can be used as potential biomarkers for the early diagnosis of DR.

A molecular dynamics investigation of N-glycosylation effects on T-cell receptor kinetics

The binding interaction between the T-cell receptor (TCR) and peptide-major histocompatibility complex (pMHC) is modulated by several factors (known and unknown), however, investigations into effects of glycosylation are limited. A fully glycosylated computational model of the TCR bound to the pMHC is developed to investigate the effects of glycosylation on dissociation kinetics from the pMHC. Here, we examine the effects of N-glycosylation on TCR-pMHC bond strength using steered molecular dynamic simulations. N-glycosylation is a post-translational modification that adds sugar moieties to molecules and can modulate the activity of several immune molecules. Using a TCR-pMHC pair found in melanoma as a case study, our study demonstrates that N-glycosylation of the TCR-pMHC alters the proteins' conformation; increases the bond lifetime; and increases the number of hydrogen bonds and Lennard-Jones Contacts involved in the TCR-pMHC bond. We find that weak glycan-protein or glycan-glycan interactions impact the equilibrated structure of the TCR and pMHC leading to an increase in the overall bond strength of the TCR-pMHC complex including the duration and energetic strength under constant load. These results indicate that N-glycosylation plays an important role in the TCR-pMHC bond and should be considered in future computational and experimental studies.Communicated by Ramaswamy H. Sarma.

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.

Phagocytosis by an HIV antibody is associated with reduced viremia irrespective of enhanced complement lysis

Increasingly, antibodies are being used to treat and prevent viral infections. In the context of HIV, efficacy is primarily attributed to dose-dependent neutralization potency and to a lesser extent Fc-mediated effector functions. It remains unclear whether augmenting effector functions of broadly neutralizing antibodies (bNAbs) may improve their clinical potential. Here, we use bNAb 10E8v4 targeting the membrane external proximal region (MPER) to examine the role of antibody-mediated effector and complement (C’) activity when administered prophylactically against SHIV challenge in rhesus macaques. With sub-protective dosing, we find a 78–88% reduction in post-acute viremia that is associated with 10E8v4-mediated phagocytosis acting at the time of challenge. Neither plasma nor tissue viremic outcomes in vivo is improved with an Fc-modified variant of 10E8v4 enhanced for C’ functions as determined in vitro. These results suggest that effector functions inherent to unmodified 10E8v4 contribute to efficacy against SHIV_SF162P3 in the absence of plasma neutralizing titers, while C’ functions are dispensable in this setting, informing design of bNAb modifications for improving protective efficacy.

While antibodies neutralize HIV via Fab recognition of viral surface antigens, antibody Fc domains mediate effector functions, including antibody-dependent cellular phagocytosis (ADCP) and cytotoxicity (ADCC), and complement (C') activity. Here, Spencer et al. modify bNAb 10E8v4 to enhance C'-mediated potency in SHIV challenged rhesus macaques to probe its function in protection, showing that in the absence of neutralization, enhancing C' activities in vitro adds no value toward reducing viremia in either blood or tissue.

Insights into substrate recognition and specificity for IgG by Endoglycosidase S2

Antibodies bind foreign antigens with high affinity and specificity leading to their neutralization and/or clearance by the immune system. The conserved N-glycan on IgG has significant impact on antibody effector function, with the endoglycosidases of Streptococcus pyogenes deglycosylating the IgG to evade the immune system, a process catalyzed by the endoglycosidase EndoS2. Studies have shown that two of the four domains of EndoS2, the carbohydrate binding module (CBM) and the glycoside hydrolase (GH) domain are critical for catalytic activity. To yield structural insights into contributions of the CBM and the GH domains as well as the overall flexibility of EndoS2 to the proteins’ catalytic activity, models of EndoS2-Fc complexes were generated through enhanced-sampling molecular-dynamics (MD) simulations and site-identification by ligand competitive saturation (SILCS) docking followed by reconstruction and multi-microsecond MD simulations. Modeling results predict that EndoS2 initially interacts with the IgG through its CBM followed by interactions with the GH yielding catalytically competent states. These may involve the CBM and GH of EndoS2 simultaneously interacting with either the same Fc CH2/CH3 domain or individually with the two Fc CH2/CH3 domains, with EndoS2 predicted to assume closed conformations in the former case and open conformations in the latter. Apo EndoS2 is predicted to sample both the open and closed states, suggesting that either complex can directly form following initial IgG-EndoS2 encounter. Interactions of the CBM and GH domains with the IgG are predicted to occur through both its glycan and protein regions. Simulations also predict that the Fc glycan can directly transfer from the CBM to the GH, facilitating formation of catalytically competent complexes and how the 734 to 751 loop on the CBM can facilitate extraction of the glycan away from the Fc CH2/CH3 domain. The predicted models are compared and consistent with Hydrogen/Deuterium Exchange data. In addition, the complex models are consistent with the high specificity of EndoS2 for the glycans on IgG supporting the validity of the predicted models.

The pathogen Streptococcus pyogenes uses the endoglycosidases S and S2 to cleave the glycans on the Fc portion of IgG antibodies, leading to a decreased cytotoxicity of the antibodies, thereby evading the host immune response. To identify potential structures of the complex of EndoS2 with IgG that could lead to the catalytic hydrolysis of the IgG glycan, molecular modeling and molecular dynamics simulations were applied. The resulting structural models predict that EndoS2 initially interacts through its carbohydrate binding module (CBM) with the IgG with subsequent interactions with the catalytic glycoside hydrolase (GH) domain yielding stable complexes. In the modeled complexes the CBM and the GH interact either simultaneously with the same Fc CH2/CH3 domain or with the two individual Fc CH2/CH3 domains separately to yield potentially catalytically competent species. In addition, apo EndoS2 is shown to assume both open and closed conformations allowing it to directly form either type of complex from which deglycosylation of either mono- or diglycosylated IgG species may occur.

Dynamic Network Modeling of Allosteric Interactions and Communication Pathways in the SARS-CoV-2 Spike Trimer Mutants: Differential Modulation of Conformational Landscapes and Signal Transmission via Cascades of Regulatory Switches

The rapidly growing body of structural and biochemical studies of the SARS-CoV-2 spike glycoprotein has revealed a variety of distinct functional states with radically different arrangements of the receptor-binding domain, highlighting a remarkable function-driven conformational plasticity and adaptability of the spike proteins. In this study, we examined molecular mechanisms underlying conformational and dynamic changes in the SARS-CoV-2 spike mutant trimers through the lens of dynamic analysis of allosteric interaction networks and atomistic modeling of signal transmission. Using an integrated approach that combined coarse-grained molecular simulations, protein stability analysis, and perturbation-based modeling of residue interaction networks, we examined how mutations in the regulatory regions of the SARS-CoV-2 spike protein can differentially affect dynamics and allosteric signaling in distinct functional states. The results of this study revealed key functional regions and regulatory centers that govern collective dynamics, allosteric interactions, and control signal transmission in the SARS-CoV-2 spike proteins. We found that the experimentally confirmed regulatory hotspots that dictate dynamic switching between conformational states of the SARS-CoV-2 spike protein correspond to the key hinge sites and global mediating centers of the allosteric interaction networks. The results of this study provide a novel insight into allosteric regulatory mechanisms of SARS-CoV-2 spike proteins showing that mutations at the key regulatory positions can differentially modulate distribution of states and determine topography of signal communication pathways operating through state-specific cascades of control switch points. This analysis provides a plausible strategy for allosteric probing of the conformational equilibrium and therapeutic intervention by targeting specific hotspots of allosteric interactions and communications in the SARS-CoV-2 spike proteins.

Immunoglobulin G Glycosylation Changes in Aging and Other Inflammatory Conditions

Among the multiple roles played by protein glycosylation, the fine regulation of biological interactions is one of the most important. The asparagine 297 (Asn₂₉₇) of IgG heavy chains is decorated by a diantennary glycan bearing a number of galactose and sialic acid residues on the branches ranging from 0 to 2. In addition, the structure can present core-linked fucose and/or a bisecting GlcNAc. In many inflammatory and autoimmune conditions, as well as in metabolic, cardiovascular, infectious, and neoplastic diseases, the IgG Asn₂₉₇-linked glycan becomes less sialylated and less galactosylated, leading to increased expression of glycans terminating with GlcNAc. These conditions alter also the presence of core-fucose and bisecting GlcNAc. Importantly, similar glycomic alterations are observed in aging. The common condition, shared by the above-mentioned pathological conditions and aging, is a low-grade, chronic, asymptomatic inflammatory state which, in the case of aging, is known as inflammaging. Glycomic alterations associated with inflammatory diseases often precede disease onset and follow remission. The aberrantly glycosylated IgG glycans associated with inflammation and aging can sustain inflammation through different mechanisms, fueling a vicious loop. These include complement activation, Fcγ receptor binding, binding to lectin receptors on antigen-presenting cells, and autoantibody reactivity. The complex molecular bases of the glycomic changes associated with inflammation and aging are still poorly understood.

Molecular Simulations and Network Modeling Reveal an Allosteric Signaling in the SARS-CoV-2 Spike Proteins

The development of computational strategies for the quantitative characterization of the functional mechanisms of SARS-CoV-2 spike proteins is of paramount importance in efforts to accelerate the discovery of novel therapeutic agents and vaccines combating the COVID-19 pandemic. Structural and biophysical studies have recently characterized the conformational landscapes of the SARS-CoV-2 spike glycoproteins in the prefusion form, revealing a spectrum of stable and more dynamic states. By employing molecular simulations and network modeling approaches, this study systematically examined functional dynamics and identified the regulatory centers of allosteric interactions for distinct functional states of the wild-type and mutant variants of the SARS-CoV-2 prefusion spike trimer. This study presents evidence that the SARS-CoV-2 spike protein can function as an allosteric regulatory engine that fluctuates between dynamically distinct functional states. Perturbation-based modeling of the interaction networks revealed a key role of the cross-talk between the effector hotspots in the receptor binding domain and the fusion peptide proximal region of the SARS-CoV-2 spike protein. The results have shown that the allosteric hotspots of the interaction networks in the SARS-CoV-2 spike protein can control the dynamic switching between functional conformational states that are associated with virus entry to the host receptor. This study offers a useful and novel perspective on the underlying mechanisms of the SARS-CoV-2 spike protein through the lens of allosteric signaling as a regulatory apparatus of virus transmission that could open up opportunities for targeted allosteric drug discovery against SARS-CoV-2 proteins and contribute to the rapid response to the current and potential future pandemic scenarios.

Analysis of the SARS-CoV-2 spike protein glycan shield reveals implications for immune recognition

Here we have generated 3D structures of glycoforms of the spike (S) glycoprotein from SARS-CoV-2, based on reported 3D structures and glycomics data for the protein produced in HEK293 cells. We also analyze structures for glycoforms representing those present in the nascent glycoproteins (prior to enzymatic modifications in the Golgi), as well as those that are commonly observed on antigens present in other viruses. These models were subjected to molecular dynamics (MD) simulation to determine the extent to which glycan microheterogeneity impacts the antigenicity of the S glycoprotein. Lastly, we have identified peptides in the S glycoprotein that are likely to be presented in human leukocyte antigen (HLA) complexes, and discuss the role of S protein glycosylation in potentially modulating the innate and adaptive immune response to the SARS-CoV-2 virus or to a related vaccine. The 3D structures show that the protein surface is extensively shielded from antibody recognition by glycans, with the notable exception of the ACE2 receptor binding domain, and also that the degree of shielding is largely insensitive to the specific glycoform. Despite the relatively modest contribution of the glycans to the total molecular weight of the S trimer (17% for the HEK293 glycoform) they shield approximately 40% of the protein surface.

Development and simulation of fully glycosylated molecular models of ACE2-Fc fusion proteins and their interaction with the SARS-CoV-2 spike protein binding domain

We develop fully glycosylated computational models of ACE2-Fc fusion proteins which are promising targets for a COVID-19 therapeutic. These models are tested in their interaction with a fragment of the receptor-binding domain (RBD) of the Spike Protein S of the SARS-CoV-2 virus, via atomistic molecular dynamics simulations. We see that some ACE2 glycans interact with the S fragments, and glycans are influencing the conformation of the ACE2 receptor. Additionally, we optimize algorithms for protein glycosylation modelling in order to expedite future model development. All models and algorithms are openly available.

TETRALEC, Artificial Tetrameric Lectins: A Tool to Screen Ligand and Pathogen Interactions

C-type lectin receptor (CLR)/carbohydrate recognition occurs through low affinity interactions. Nature compensates that weakness by multivalent display of the lectin carbohydrate recognition domain (CRD) at the cell surface. Mimicking these low affinity interactions in vitro is essential to better understand CLR/glycan interactions. Here, we present a strategy to create a generic construct with a tetrameric presentation of the CRD for any CLR, termed TETRALEC. We applied our strategy to a naturally occurring tetrameric CRD, DC-SIGNR, and compared the TETRALEC ligand binding capacity by synthetic N- and O-glycans microarray using three different DC-SIGNR constructs i) its natural tetrameric counterpart, ii) the monomeric CRD and iii) a dimeric Fc-CRD fusion. DC-SIGNR TETRALEC construct showed a similar binding profile to that of its natural tetrameric counterpart. However, differences observed in recognition of low affinity ligands underlined the importance of the CRD spatial arrangement. Moreover, we further extended the applications of DC-SIGNR TETRALEC to evaluate CLR/pathogens interactions. This construct was able to recognize heat-killed Candida albicans by flow cytometry and confocal microscopy, a so far unreported specificity of DC-SIGNR. In summary, the newly developed DC-SIGNR TETRALEC tool proved to be useful to unravel novel CLR/glycan interactions, an approach which could be applied to other CLRs.

Analysis of the SARS-CoV-2 spike protein glycan shield: implications for immune recognition

Here we have generated 3D structures of glycoforms of the spike (S) glycoprotein from SARS-CoV-2, based on reported 3D structures and glycomics data for the protein produced in HEK293 cells. We also analyze structures for glycoforms representing those present in the nascent glycoproteins (prior to enzymatic modifications in the Golgi), as well as those that are commonly observed on antigens present in other viruses. These models were subjected to molecular dynamics (MD) simulation to determine the extent to which glycan microheterogeneity impacts the antigenicity of the S glycoprotein. Lastly, we have identified peptides in the S glycoprotein that are likely to be presented in human leukocyte antigen (HLA) complexes, and discuss the role of S protein glycosylation in potentially modulating the adaptive immune response to the SARS-CoV-2 virus or to a related vaccine. The 3D structures show that the protein surface is extensively shielded from antibody recognition by glycans, with the exception of the ACE2 receptor binding domain, and also that the degree of shielding is largely insensitive to the specific glycoform. Despite the relatively modest contribution of the glycans to the total molecular weight (17% for the HEK293 glycoform) the level of surface shielding is disproportionately high at 42%.

Going Native: Synthesis of Glycoproteins and Glycopeptides via Native Linkages To Study Glycan-Specific Roles in the Immune System

Glycosylation plays a myriad of roles in the immune system: Certain glycans can interact with specific immune receptors to kickstart a pro-inflammatory response, whereas other glycans can do precisely the opposite and ameliorate the immune response. Specific glycans and glycoforms can themselves become the targets of the adaptive immune system, leading to potent antiglycan responses that can lead to the killing of altered self- or pathogenic species. This hydra-like set of roles glycans play is of particular importance in cancer immunity, where it influences the anticancer immune response, likely playing pivotal roles in tumor survival or clearance. The complexity of carbohydrate biology requires synthetic access to glycoproteins and glycopeptides that harbor homogeneous glycans allowing the probing of these systems with high precision. One particular complicating factor in this is that these synthetic structures are required to be as close to the native structures as possible, as non-native linkages can themselves elicit immune responses. In this Review, we discuss examples and current strategies for the synthesis of natively linked single glycoforms of peptides and proteins that have enabled researchers to gain new insights into glycoimmunology, with a particular focus on the application of these reagents in cancer immunology.