Influenza poses a major health issue globally. Neutralizing antibodies targeting the highly conserved stem region of hemagglutinin (HA) of the influenza virus provide considerable protection against the infection. Using an array of advanced simulation technologies, we developed a high-resolution structural model of full-length, Fab-bound HA in a native viral membrane to characterize direct membrane interactions that govern the efficacy of the antibody. We reveal functionally important residues beyond the antibody's complementary-determining regions that contribute to its membrane binding. Mutagenesis experiments and infectivity assays confirm that deactivating the membrane-binding residues of the antibody decreases its neutralization activity. Therefore, we propose that the association with the viral membrane plays a key role in the neutralization activity of these antibodies. Given the rapid evolution of the influenza virus, the developed model provides a structural framework for the rational design and development of more effective therapeutic antibodies.

Glycosaminoglycans (GAGs) play a key role in a variety of biological processes in the extracellular matrix (ECM) via interactions with their protein targets. Due to their high flexibility, periodicity and electrostatics-driven interactions, GAG-containing complexes are very challenging to characterize both experimentally and in silico. In this study, we, for the first time, systematically analyzed the interactions of endostatin, a proteolytic fragment of collagen XVIII known to be anti-angiogenic and anti-tumoral, with heparin (HP) and representative heparan sulfate (HS) oligosaccharides of various lengths, sequences and sulfation patterns. We first used conventional molecular docking and a docking approach based on a repulsive scaling–replica exchange molecular dynamics technique, as well as unbiased molecular dynamic simulations, to obtain dynamically stable GAG binding poses. Then, the corresponding free energies of binding were calculated and the amino acid residues that contribute the most to GAG binding were identified. We also investigated the potential influence of Zn2+ on endostatin–HP complexes using computational approaches. These data provide new atomistic details of the molecular mechanism of HP’s binding to endostatin, which will contribute to a better understanding of its interplay with proteoglycans at the cell surface and in the extracellular matrix.

ABSTRACT Ebola virus glycoprotein (EBOV GP) is one of the most heavily O-glycosylated viral glycoproteins, yet we still lack a fundamental understanding of the structure of its large O-glycosylated mucin-like domain and to what degree the host O-glycosylation capacity influences EBOV replication. Using tandem mass spectrometry, we identified 47 O-glycosites on EBOV GP and found similar glycosylation signatures on virus-like particle- and cell lysate-derived GP. Furthermore, we performed quantitative differential O-glycoproteomics on proteins produced in wild-type HEK293 cells and cell lines ablated for the three key initiators of O-linked glycosylation, GalNAc-T1, -T2, and -T3. The data show that 12 out of the 47 O-glycosylated sites were regulated, predominantly by GalNAc-T1. Using the glycoengineered cell lines for authentic EBOV propagation, we demonstrate the importance of O-linked glycan initiation and elongation for the production of viral particles and the titers of progeny virus. The mapped O-glycan positions and structures allowed to generate molecular dynamics simulations probing the largely unknown spatial arrangements of the mucin-like domain. The data highlight targeting GALNT1 or C1GALT1C1 as a possible way to modulate O-glycan density on EBOV GP for novel vaccine designs and tailored intervention approaches. IMPORTANCE Ebola virus glycoprotein acquires its extensive glycan shield in the host cell, where it is decorated with N-linked glycans and mucin-type O-linked glycans. The latter is initiated by a family of polypeptide GalNAc-transferases that have different preferences for optimal peptide substrates resulting in a spectrum of both very selective and redundant substrates for each isoform. In this work, we map the exact locations of O-glycans on Ebola virus glycoprotein and identify subsets of sites preferentially initiated by one of the three key isoforms of GalNAc-Ts, demonstrating that each enzyme contributes to the glycan shield integrity. We further show that altering host O-glycosylation capacity has detrimental effects on Ebola virus replication, with both isoform-specific initiation and elongation playing a role. The combined structural and functional data highlight glycoengineered cell lines as useful tools for investigating molecular mechanisms imposed by specific glycans and for steering the immune responses in future vaccine designs.

Siglec-7 regulates immune cell activity and is a promising target for immunomodulation. Here, we report the discovery of novel sialic acid derivatives binding to Siglec-7. Synthesis and affinity measurements are complemented by high-quality models of sialoside-Siglec-7 complexes based on molecular dynamics (MD) simulations on the microsecond time scale. We provide details for the predicted binding modes for the new ligands, e.g., that an extension of the carbon backbone leads to a different molecular interaction pattern with the receptor and the nearby water structure than found for known Siglec-7 ligands. Further on, we uncover some shortcomings of the GLYCAM06 and GAFF2 force fields when used for the simulation of sialoside-based glycomimetics. Our results open new opportunities for the rational design of Siglec-7 inhibitors. In addition, we provide strategies on how to use and visualize MD simulations to describe and investigate sialoside-Siglec complexes in general.

Influenza poses a major health issue globally. Characterized by annual epidemics with mild to severe symptoms and even a significant number of deaths, the virus has historically resulted in substantial economic and societal problems in the world. Neutralizing antibodies that bind to the membrane-proximal region of hemagglutinin (HA) provide considerable protection against the infection by blocking the fusion of virus to the host cell. Here, using an array of advanced modeling and simulation techniques employed to a very large molecular assembly, for the first time, we develop a high-resolution structural model of full-length, antibody-bound HA in a native viral membrane and describe key molecular interactions that determine the binding affinity of the antibody to the viral membrane.

As small and stable high-affinity antigen binders, VHHs boast attractive characteristics both for therapeutic use in various disease indications, and as versatile reagents in research and diagnostics. To further increase the versatility of VHHs, we explored the VHH scaffold in a structure-guided approach to select regions where the introduction of an N-glycosylation N-X-T sequon and its associated glycan should not interfere with protein folding or epitope recognition. We expressed variants of such glycoengineered VHHs in the Pichia pastoris GlycoSwitchM5 strain, allowing us to pinpoint preferred sites at which Man 5 GlcNAc 2 -glycans can be introduced at high site occupancy without affecting antigen binding. A VHH carrying predominantly a Man 5 GlcNAc 2 N-glycan at one of these preferred sites showed highly efficient, glycan-dependent uptake by Mf4/4 macrophages in vitro and by alveolar lung macrophages in vivo, illustrating one potential application of glyco-engineered VHHs: a glycan-based targeting approach for lung macrophage endolysosomal system delivery. The set of optimal artificial VHH N-glycosylation sites identified in this study can serve as a blueprint for targeted glyco-engineering of other VHHs, enabling site-specific functionalization through the rapidly expanding toolbox of synthetic glycobiology.

Nonulosonic acids or non-2-ulosonic acids (NulOs) are an ancient family of 2-ketoaldonic acids (α-ketoaldonic acids) with a 9-carbon backbone. In nature, these monosaccharides occur either in a 3-deoxy form (referred to as “sialic acids”) or in a 3,9-dideoxy “sialic-acid-like” form. The former sialic acids are most common in the deuterostome lineage, including vertebrates, and mimicked by some of their pathogens. The latter sialic-acid-like molecules are found in bacteria and archaea. NulOs are often prominently positioned at the outermost tips of cell surface glycans, and have many key roles in evolution, biology and disease. The diversity of stereochemistry and structural modifications among the NulOs contributes to more than 90 sialic acid forms and 50 sialic-acid-like variants described thus far in nature. This paper reports the curation of these diverse naturally occurring NulOs at the NCBI sialic acid page (https://www.ncbi.nlm.nih.gov/glycans/sialic.html) as part of the NCBI-Glycans initiative. This includes external links to relevant Carbohydrate Structure Databases. As the amino and hydroxyl groups of these monosaccharides are extensively derivatized by various substituents in nature, the Symbol Nomenclature For Glycans (SNFG) rules have been expanded to represent this natural diversity. These developments help illustrate the natural diversity of sialic acids and related NulOs, and enable their systematic representation in publications and online resources.

Coronavirus spike proteins mediate receptor binding and membrane fusion, making them prime targets for neutralizing antibodies. In the cases of severe acute respiratory syndrome coronavirus, severe acute respiratory syndrome coronavirus 2 and Middle East respiratory syndrome coronavirus, spike proteins transition freely between open and closed conformations to balance host cell attachment and immune evasion 1 – 5 . Spike opening exposes domain S1 B , allowing it to bind to proteinaceous receptors 6 , 7 , and is also thought to enable protein refolding during membrane fusion 4 , 5 . However, with a single exception, the pre-fusion spike proteins of all other coronaviruses studied so far have been observed exclusively in the closed state. This raises the possibility of regulation, with spike proteins more commonly transitioning to open states in response to specific cues, rather than spontaneously. Here, using cryogenic electron microscopy and molecular dynamics simulations, we show that the spike protein of the common cold human coronavirus HKU1 undergoes local and long-range conformational changes after binding a sialoglycan-based primary receptor to domain S1 A . This binding triggers the transition of S1 B domains to the open state through allosteric interdomain crosstalk. Our findings provide detailed insight into coronavirus attachment, with possibilities of dual receptor usage and priming of entry as a means of immune escape.

Glycans of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein are speculated to play functional roles in the infection processes as they extensively cover the protein surface and are highly conserved across the variants. The spike protein has been the principal target for vaccine and therapeutic development while the exact effects of its glycosylation remain elusive. Analytical reports have described the glycan heterogeneity of the spike protein. Subsequent molecular simulation studies provided a knowledge basis of the glycan functions. However, experimental data on the role of discrete glycoforms on the spike protein pathobiology remains scarce. Building an understanding of their roles in SARS-CoV-2 is important as we continue to develop effective medicines and vaccines to combat the disease. Herein, we used designed combinations of glycoengineering enzymes to simplify and control the glycosylation profile of the spike protein receptor-binding domain (RBD). Measurements of the receptor-binding affinity revealed opposite regulatory effects of the RBD glycans with and without sialylation, which presents a potential strategy for modulating the spike protein behaviors through glycoengineering. Moreover, we found that the reported anti-SARS-CoV-(2) antibody, S309, neutralizes the impact of different RBD glycoforms on the receptor-binding affinity. In combination with molecular dynamics simulation, this work reports the regulatory roles that glycosylation plays in the interaction between the viral spike protein and host receptor, providing new insights into the nature of SARS-CoV-2. Beyond this study, enzymatic glycan remodeling offers the opportunity to understand the fundamental role of specific glycoforms on glycoconjugates across molecular biology.

We have used NMR experiments to explore binding of selected glycans and glycomimetics to the SARS CoV‑2 spike glycoprotein (S‑protein) and to its receptor binding domain (RBD). STD NMR experiments confirm binding of sialoglycans to the S‑protein of the prototypic Wuhan strain virus and yield dissociation constants in the mM range. The absence of STD effects for sialoglycans in the presence of the Omicron/BA.1 S‑protein reflects a loss of binding as a result of S‐protein evolution. Likewise, no STD effects are observed for the deletion mutant Δ 143‐145 of the Wuhan S‐protein, supporting localization of the binding site in the N‐terminal domain (NTD). The glycomimetics Oseltamivir and Zanamivir bind weakly to the S‑protein of both virus strains. Binding of blood group antigens to the Wuhan S‑protein cannot be confirmed by STD NMR. Using 1 H, 15 N TROSY HSQC based chemical shift perturbation (CSP) experiments we exclude binding of any of the ligands studied to the RBD of the Wuhan S‑protein. Our results put reported data on glycan binding into perspective and shed new light on the potential role of glycan‐binding to the S‐protein.

SummaryEbola virus glycoprotein is one of the most heavily O-glycosylated viral envelope glycoproteins. Using glycoengineered cell lines we demonstrate that O-linked glycan truncation and perturbed initiation retarded the production of viral particles and decreased infectivity of progeny virus. Next, using TMT isobaric labelling, we performed quantitative differential O-glycoproteomics on proteins produced in wild type HEK293 cells and cell lines ablated for three key GalNAc-transferases, GalNAc-T1, -T2, and -T3, as well as compared it to patterns on wild type virus-like particles. We demonstrate selective initiation of a subset of O-glycosites by each enzyme, with GalNAc-T1 having the largest contribution. This work represents a comprehensive site-specific analysis of EBOV GP, with 47 O-glycosites identified, and sheds light on differential regulation of EBOV GP glycosylation initiated by host GalNAc-Ts. Together with the effect on viral propagation it opens prospective avenues for tailored intervention approaches and means for modulating immunogen O-glycan density.

We have used chemical shift perturbation (CSP) and saturation transfer difference (STD) NMR experiments to identify and characterize the binding of selected ligands to the receptor-binding domain (RBD) of the spike glycoprotein (S-protein) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We also subjected full-length S-protein to STD NMR experiments, allowing correlations with RBD-based results. CSPs reveal the binding sites for heparin and fondaparinux, and affinities were measured using CSP titrations. We then show that α-2,3-sialyllactose binds to the S-protein but not to the RBD. Finally, combined CSP and STD NMR experiments show that lifitegrast, a compound used for the treatment of dry eye, binds to the linoleic acid (LA) binding pocket with a dissociation constant in the μM range. This is an interesting finding, as lifitegrast lends itself well as a blueprint for medicinal chemistry, eventually furnishing novel entry inhibitors targeting the highly conserved LA binding site.

As coronavirus disease 2019 (COVID-19) persists, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) emerge, accumulating spike (S) glycoprotein mutations. S receptor binding domain (RBD) comprises a free fatty acid (FFA)–binding pocket. FFA binding stabilizes a locked S conformation, interfering with virus infectivity. We provide evidence that the pocket is conserved in pathogenic β-coronaviruses (β-CoVs) infecting humans. SARS-CoV, MERS-CoV, SARS-CoV-2, and VOCs bind the essential FFA linoleic acid (LA), while binding is abolished by one mutation in common cold–causing HCoV-HKU1. In the SARS-CoV S structure, LA stabilizes the locked conformation, while the open, infectious conformation is devoid of LA. Electron tomography of SARS-CoV-2–infected cells reveals that LA treatment inhibits viral replication, resulting in fewer deformed virions. Our results establish FFA binding as a hallmark of pathogenic β-CoV infection and replication, setting the stage for FFA-based antiviral strategies to overcome COVID-19.

Significant interest in the development of high-affinity ligands for Siglecs exists due to the various therapeutically relevant functions of these proteins. Here, we report a new strategy to develop and design Siglec ligands as disialyl-oligosaccharide mimetics exemplified on Siglec-2 (CD22). We report insights into development of dimeric ligands with high affinity and avidity to cell surface-expressed CD22, assay development, tool compounds, structure activity relationships, and biological data on calcium flux regulation in B-cells. The binding modes of selected ligands have been modeled based on state-of-the-art molecular dynamics simulations on the microsecond timescale, providing detailed views on ligand binding and opening a new perspective on drug design efforts for Siglecs. High-avidity dimeric ligands containing a linker opening the way towards bispecifics are presented as well.

The myelomonocytic receptor CD33 (Siglec-3) inhibits innate immune reactivity by extracellular V-set domain recognition of sialic acid (Sia)-containing “self-associated molecular patterns” (SAMPs). We earlier showed that V-set domain-deficient CD33-variant allele, protective against late-onset Alzheimer’s Disease (LOAD), is derived and specific to the hominin lineage. We now report multiple hominin-specific CD33 V-set domain mutations. Due to hominin-specific, fixed loss-of-function mutation in the CMAH gene, humans lack N-glycolylneuraminic acid (Neu5Gc), the preferred Sia-ligand of ancestral CD33. Mutational analysis and molecular dynamics (MD)-simulations indicate that fixed change in amino acid 21 of hominin V-set domain and conformational changes related to His45 corrected for Neu5Gc-loss by switching to N-acetylneuraminic acid (Neu5Ac)-recognition. We show that human-specific pathogens Neisseria gonorrhoeae and Group B Streptococcus selectively bind human CD33 (huCD33) as part of immune-evasive molecular mimicry of host SAMPs and that this binding is significantly impacted by amino acid 21 modification. In addition to LOAD-protective CD33 alleles, humans harbor derived, population-universal, cognition-protective variants at several other loci. Interestingly, 11 of 13 SNPs in these human genes (including CD33) are not shared by genomes of archaic hominins: Neanderthals and Denisovans. We present a plausible evolutionary scenario to compile, correlate, and comprehend existing knowledge about huCD33-evolution and suggest that grandmothering emerged in humans.

Population analysis in terms of glycosidic torsion angles is frequently used to reveal preferred conformers of glycans. However, due to high structural diversity and flexibility of carbohydrates, conformational characterization of complex glycans can be a challenging task. Herein, we present a conformation module of oligosaccharide fragments occurring in natural glycan structures developed on the platform of the Carbohydrate Structure Database. Currently, this module deposits free energy surface and conformer abundance maps plotted as a function of glycosidic torsions for 194 "inter"residue bonds. Data are automatically and continuously derived from explicit-solvent molecular dynamics (MD) simulations. The module was also supplemented with high-temperature MD data of saccharides (2,403 maps) provided by GlycoMapsDB (hosted by GLYCOSCIENCES.de project). Conformational data defined by up to 4 torsional degrees of freedom can be freely explored using a web interface of the module available at http://csdb.glycoscience.ru/database/core/search_conf.html.

The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) presents a public health crisis, and the vaccines that can induce highly potent neutralizing antibodies are essential for ending the pandemic. The spike (S) protein on the viral envelope mediates human angiotensin-converting enzyme 2 binding and thus is the target of a variety of neutralizing antibodies. In this work, we built various S trimer-antibody complex structures on the basis of the fully glycosylated S protein models described in our previous work and performed all-atom molecular dynamics simulations to gain insight into the structural dynamics and interactions between S protein and antibodies. Investigation of the residues critical for S-antibody binding allows us to predict the potential influence of mutations in SARS-CoV-2 variants. Comparison of the glycan conformations between S-only and S-antibody systems reveals the roles of glycans in S-antibody binding. In addition, we explored the antibody binding modes and the influences of antibody on the motion of S protein receptor binding domains. Overall, our analyses provide a better understanding of S-antibody interactions, and the simulation-based S-antibody interaction maps could be used to predict the influences of S mutation on S-antibody interactions, which will be useful for the development of vaccine and antibody-based therapy.

Angiotensin-converting enzyme 2 (ACE2) is the main entry point in airway epithelial cells for SARS-CoV-2. ACE2 binding to the SARS-CoV-2 protein spike triggers viral fusion with the cell plasma membrane, resulting in viral RNA genome delivery into the host. Despite ACE2's critical role in SARS-CoV-2 infection, full understanding of ACE2 expression, including in response to viral infection, remains unclear. ACE2 was thought to encode five transcripts and one protein of 805 amino acids. In the present study, we identify a novel short isoform of ACE2 expressed in the airway epithelium, the main site of SARS-CoV-2 infection. Short ACE2 is substantially upregulated in response to interferon stimulation and rhinovirus infection, but not SARS-CoV-2 infection. This short isoform lacks SARS-CoV-2 spike high-affinity binding sites and, altogether, our data are consistent with a model where short ACE2 is unlikely to directly contribute to host susceptibility to SARS-CoV-2 infection.

Enveloped viruses hijack not only the host translation processes, but also its glycosylation machinery, and to a variable extent cover viral surface proteins with tolerogenic host-like structures. SARS-CoV-2 surface protein S presents as a trimer on the viral surface and is covered by a dense shield of N-linked glycans, and a few O-glycosites have been reported. The location of O-glycans is controlled by a large family of initiating enzymes with variable expression in cells and tissues and hence is difficult to predict. Here, we used our well-established O-glycoproteomic workflows to map the precise positions of O-linked glycosylation sites on three different entities of protein S—insect cell or human cell-produced ectodomains, or insect cell derived receptor binding domain (RBD). In total 25 O-glycosites were identified, with similar patterns in the two ectodomains of different cell origin, and a distinct pattern of the monomeric RBD. Strikingly, 16 out of 25 O-glycosites were located within three amino acids from known N-glycosites. However, O-glycosylation was primarily found on peptides that were unoccupied by N-glycans, and otherwise had low overall occupancy. This suggests possible complementary functions of O-glycans in immune shielding and negligible effects of O-glycosylation on subunit vaccine design for SARS-CoV-2.

C-Mannosylation is a relatively rare form of protein glycosylation involving the attachment of an α-mannopyranosyl residue to C-2 of the indole moiety of the amino acid tryptophan. This type of linkage was initially discovered in RNase 2 from human urine but later confirmed to be present in many other important proteins. Based on NMR experiments and extensive molecular dynamics simulations on the hundred microsecond timescale we demonstrate that, for isolated glycopeptides and denatured RNase 2, the C-linked mannopyranosyl residue exists as an ensemble of conformations, among which 1C4 is the most abundant. However, for native RNase 2, molecular dynamics and NMR studies revealed that the mannopyranosyl residue favors a specific conformation, which optimally stabilizes the protein fold through a network of hydrogen bonds and which leads to a significant reduction of the protein dynamics on the microsecond timescale. Our findings contribute to the understanding of the biological role of C-mannosylation.