Restricted N-glycan conformational space in the PDB and its implication in glycan structure modeling

Molecular Modeling of Glycosylated Catalytic Domain of Human Carbonic Anhydrase IX

Glycans exhibit significant structural diversity due to the flexibility of glycosidic bonds linking their constituent monosaccharides and the formation of numerous hydrogen bonds. The present work searches a simulated ensemble of glycan chain conformations attached to the catalytic domain of N-glycosylated human carbonic anhydrase IX (HCA IX-c) to identify conformations pointed away or back-folded toward the protein surface guided by different amino acid residues. A series of classical molecular dynamics (MD) simulation studies for a total of 30 μs followed by accelerated MD simulations for a total of 2 μs have been performed using two different force fields to capture varying degrees of fluctuations of both glycan chain and HCA IX. From the underlying free energy profile and kinetics derived using hidden Markov state model, several stable glycan orientations are identified that extend away from the protein surface and convert among each other with rate constants of the order 107-1010 S-1. Most importantly, we have identified a rare glycan conformation which reaches close to a catalytically important amino acid residue, Glu-106. We further enlist the protein residues that couple such less frequent event of the glycan chain back-folding toward the surface of the protein.

Generalized biomolecular modeling and design with RoseTTAFold All-Atom

Deep-learning methods have revolutionized protein structure prediction and design but are presently limited to protein-only systems. We describe RoseTTAFold All-Atom (RFAA), which combines a residue-based representation of amino acids and DNA bases with an atomic representation of all other groups to model assemblies that contain proteins, nucleic acids, small molecules, metals, and covalent modifications, given their sequences and chemical structures. By fine-tuning on denoising tasks, we developed RFdiffusion All-Atom (RFdiffusionAA), which builds protein structures around small molecules. Starting from random distributions of amino acid residues surrounding target small molecules, we designed and experimentally validated, through crystallography and binding measurements, proteins that bind the cardiac disease therapeutic digoxigenin, the enzymatic cofactor heme, and the light-harvesting molecule bilin.

Conformational dynamics of α-1 acid glycoprotein (AGP) in cancer: A comparative study of glycosylated and unglycosylated AGP

α-1 acid glycoprotein (AGP) is one of the most abundant plasma proteins. It fulfills two important functions: immunomodulation, and binding to various drugs and receptors. These different functions are closely associated and modulated via changes in glycosylation and cancer missense mutations. From a structural point of view, glycans alter the local biophysical properties of the protein leading to a diverse ligand-binding spectrum. However, glycans can typically not be observed in the resolved X-ray crystallography structure of AGP due to their high flexibility and microheterogeneity, so limiting our understanding of AGP's conformational dynamics 70 years after its discovery. We here investigate how mutations and glycosylation interfere with AGP's conformational dynamics changing its biophysical behavior, by using molecular dynamics (MD) simulations and sequence-based dynamics predictions. The MD trajectories show that glycosylation decreases the local backbone flexibility of AGP and increases the flexibility of distant regions through allosteric effects. We observe that mutations near the glycosylation site affect glycan's conformational preferences. Thus, we conclude that mutations control glycan dynamics which modulates the protein's backbone flexibility directly affecting its accessibility. These findings may assist in the drug design targeting AGP's glycosylation and mutations in cancer.

Advances in understanding N-glycosylation structure, function, and regulation in health and disease

N-linked glycosylation is a post-translational modification crucial for membrane protein folding, stability and other cellular functions. Alteration of membrane protein N-glycans is implicated in wide range of pathological conditions including cancer metastasis, chronic inflammatory diseases, and viral pathogenesis. Even though the roles of N-glycans have been studied extensively, our knowledge of their mechanisms remains unclear due to the lack of detailed structural analysis of the N-glycome. Mapping the N-glycome landscape will open new avenues to explore disease mechanisms and identify novel therapeutic targets. This review discusses the diverse structure of N-linked glycans, the function and regulation of N-glycosylation in health and disease, and ends with a focus on recent approaches to target N-glycans in rheumatoid arthritis and cancer metastasis.

Challenges and limitations in the studies of glycoproteins: A computational chemist's perspective

Experimenters face challenges and limitations while analyzing glycoproteins due to their high flexibility, stereochemistry, anisotropic effects, and hydration phenomena. Computational studies complement experiments and have been used in characterization of the structural properties of glycoproteins. However, recent investigations revealed that computational studies face significant challenges as well. Here, we introduce and discuss some of these challenges and weaknesses in the investigations of glycoproteins. We also present requirements of future developments in computational biochemistry and computational biology areas that could be necessary for providing more accurate structural property analyses of glycoproteins using computational tools. Further theoretical strategies that need to be and can be developed are discussed herein.

Glycans in Virus-Host Interactions: A Structural Perspective

Many interactions between microbes and their hosts are driven or influenced by glycans, whose heterogeneous and difficult to characterize structures have led to an underappreciation of their role in these interactions compared to protein-based interactions. Glycans decorate microbe glycoproteins to enhance attachment and fusion to host cells, provide stability, and evade the host immune system. Yet, the host immune system may also target these glycans as glycoepitopes. In this review, we provide a structural perspective on the role of glycans in host-microbe interactions, focusing primarily on viral glycoproteins and their interactions with host adaptive immunity. In particular, we discuss a class of topological glycoepitopes and their interactions with topological mAbs, using the anti-HIV mAb 2G12 as the archetypical example. We further offer our view that structure-based glycan targeting strategies are ready for application to viruses beyond HIV, and present our perspective on future development in this area.

Dissecting the conformation of glycans and their interactions with proteins

The use of in silico strategies to develop the structural basis for a rational optimization of glycan-protein interactions remains a great challenge. This problem derives, in part, from the lack of technologies to quantitatively and qualitatively assess the complex assembling between a glycan and the targeted protein molecule. Since there is an unmet need for developing new sugar-targeted therapeutics, many investigators are searching for technology platforms to elucidate various types of molecular interactions within glycan-protein complexes and aid in the development of glycan-targeted therapies. Here we discuss three important technology platforms commonly used in the assessment of the complex assembly of glycosylated biomolecules, such as glycoproteins or glycosphingolipids: Biacore analysis, molecular docking, and molecular dynamics simulations. We will also discuss the structural investigation of glycosylated biomolecules, including conformational changes of glycans and their impact on molecular interactions within the glycan-protein complex. For glycoproteins, secreted protein acidic and rich in cysteine (SPARC), which is associated with various lung disorders, such as chronic obstructive pulmonary disease (COPD) and lung cancer, will be taken as an example showing that the core fucosylation of N-glycan in SPARC regulates protein-binding affinity with extracellular matrix collagen. For glycosphingolipids (GSLs), Globo H ceramide, an important tumor-associated GSL which is being actively investigated as a target for new cancer immunotherapies, will be used to demonstrate how glycan structure plays a significant role in enhancing angiogenesis in tumor microenvironments.

IMRPS: Inserted and Modified Residues in Protein Structures. A database

Modified residues present in proteins are the result of post-translational modifications (PTMs). These PTMs increase the functional diversity of the proteome and influence various biological processes and diseased conditions. Therefore, identification and understanding of PTMs in various protein structures is of great significance. In view of this, an online database, Inserted and Modified Residues in Protein Structures (IMRPS), has been developed. IMRPS is a derived database that furnishes information on the residues modified and inserted in the protein structures available in the Protein Data Bank (PDB). The database is equipped with a graphical user interface and has an option to view the data for non-redundant protein structures (25 and 90%) as well. A quality criteria cutoff has been incorporated to assist in displaying the specific set of PDB codes. The entire protein structure along with the inserted or modified residues can be visualized in JSmol. This database will be updated regularly (presently, every three months) and can be accessed through the URL http://cluster.physics.iisc.ac.in/imrps/.

The effect of lipidation and glycosylation on short cationic antimicrobial peptides

The global health threat surrounding bacterial resistance has resulted in antibiotic researchers shifting their focus away from 'traditional' antibiotics and concentrating on other antimicrobial agents, including antimicrobial peptides. These low molecular weight "mini-proteins" exhibit broad-spectrum activity against bacteria, including multi-drug resistant strains, viruses, fungi and protozoa and constitute a major element of the innate-immune system of many multicellular organisms. Some naturally occurring antimicrobial peptides are lipidated and/or glycosylated and almost all antimicrobial peptides in clinical use are either lipopeptides (Daptomycin and Polymyxin E and B) or glycopeptides (Vancomycin). Lipidation, glycosylation and PEGylation are an option for improving stability and activity in serum and for reducing the rapid clearing via the kidneys and liver. Two broad-spectrum antimicrobial peptides NH₂-RIRIRWIIR-CONH₂ (A1) and NH₂-KRRVRWIIW-CONH₂ (B1) were conjugated via a linker, producing A2 and B2, to individual fatty acids of C₈, C₁₀, C₁₂ and C₁₄ and in addition, A2 was conjugated to either glucose, N-acetyl glucosamine, galactose, mannose, lactose or polyethylene glycol (PEG). Antimicrobial activity against two Gram-positive strains (methicillin resistant Staphylococcus aureus (MRSA) and vancomycin resistant Enterococcus faecalis (VRE)) and three Gram-negative strains (Salmonella typhimurium, E. coli and Pseudomonas aeruginosa) were determined. Activity patterns for the lipidated versions are very complex, dependent on sequence, bacteria and fatty acid. Two reciprocal effects were measured; compared to the parental peptides, some combinations led to a 16-fold improvement whereas other combinations let to a 32-fold reduction in antimicrobial activity. Glycosylation decreased antimicrobial activity by 2 to 16-fold in comparison to A1, respectively on the sugar-peptide combination. PEGylation rendered the peptide inactive. Antimicrobial activity in the presence of 25% human serum of A1 and B1 was reduced 32-fold and 8-fold, respectively. The longer chain fatty acids almost completely restored this activity; however, these fatty acids increased hemolytic activity. B1 modified with C8 increased the therapeutic index by 2-fold for four bacterial strains. Our results suggest that finding the right lipid-peptide combination can lead to improved activity in the presence of serum and potentially more effective drug candidates for animal studies. Glycosylation with the optimal sugar and numbers of sugars at the right peptide position could be an alternative route or could be used in addition to lipidation to counteract solubility and toxicity issues.

CHARMM-GUI Glycan Modeler for modeling and simulation of carbohydrates and glycoconjugates

Characterizing glycans and glycoconjugates in the context of three-dimensional structures is important in understanding their biological roles and developing efficient therapeutic agents. Computational modeling and molecular simulation have become an essential tool complementary to experimental methods. Here, we present a computational tool, Glycan Modeler for in silico N-/O-glycosylation of the target protein and generation of carbohydrate-only systems. In our previous study, we developed Glycan Reader, a web-based tool for detecting carbohydrate molecules from a PDB structure and generation of simulation system and input files. As integrated into Glycan Reader in CHARMM-GUI, Glycan Modeler (Glycan Reader & Modeler) enables to generate the structures of glycans and glycoconjugates for given glycan sequences and glycosylation sites using PDB glycan template structures from Glycan Fragment Database (http://glycanstructure.org/fragment-db). Our benchmark tests demonstrate the universal applicability of Glycan Reader & Modeler to various glycan sequences and target proteins. We also investigated the structural properties of modeled glycan structures by running 2-μs molecular dynamics simulations of HIV envelope protein. The simulations show that the modeled glycan structures built by Glycan Reader & Modeler have the similar structural features compared to the ones solved by X-ray crystallography. We also describe the representative examples of glycoconjugate modeling with video demos to illustrate the practical applications of Glycan Reader & Modeler. Glycan Reader & Modeler is freely available at http://charmm-gui.org/input/glycan.

Conformational Heterogeneity of the HIV Envelope Glycan Shield

To better understand the conformational properties of the glycan shield covering the surface of the HIV gp120/gp41 envelope (Env) trimer, and how the glycan shield impacts the accessibility of the underlying protein surface, we performed enhanced sampling molecular dynamics (MD) simulations of a model glycosylated HIV Env protein and related systems. Our simulation studies revealed a conformationally heterogeneous glycan shield with a network of glycan-glycan interactions more extensive than those observed to date. We found that partial preorganization of the glycans potentially favors binding by established broadly neutralizing antibodies; omission of several specific glycans could increase the accessibility of other glycans or regions of the protein surface to antibody or CD4 receptor binding; the number of glycans that can potentially interact with known antibodies is larger than that observed in experimental studies; and specific glycan conformations can maximize or minimize interactions with individual antibodies. More broadly, the enhanced sampling MD simulations described here provide a valuable tool to guide the engineering of specific Env glycoforms for HIV vaccine design.

Identification and Characterization of Complex Glycosylated Peptides Presented by the MHC Class II Processing Pathway in Melanoma

The MHC class II (MHCII) processing pathway presents peptides derived from exogenous or membrane-bound proteins to CD4+ T cells. Several studies have shown that glycopeptides are necessary to modulate CD4+ T cell recognition, though glycopeptide structures in these cases are generally unknown. Here, we present a total of 93 glycopeptides from three melanoma cell lines and one matched EBV-transformed line with most found only in the melanoma cell lines. The glycosylation we detected was diverse and comprised 17 different glycoforms. We then used molecular modeling to demonstrate that complex glycopeptides are capable of binding the MHC and may interact with complementarity determining regions. Finally, we present the first evidence of disulfide-bonded peptides presented by MHCII. This is the first large scale study to sequence glyco- and disulfide bonded MHCII peptides from the surface of cancer cells and could represent a novel avenue of tumor activation and/or immunoevasion.

Delineating the conformational flexibility of trisaccharides from NMR spectroscopy experiments and computer simulations

The conformation of saccharides in solution is challenging to characterize in the context of a single well-defined three-dimensional structure. Instead, they are better represented by an ensemble of conformations associated with their structural diversity and flexibility. In this study, we delineate the conformational heterogeneity of five trisaccharides via a combination of experimental and computational techniques. Experimental NMR measurements target conformationally sensitive parameters, including J couplings and effective distances around the glycosidic linkages, while the computational simulations apply the well-calibrated additive CHARMM carbohydrate force field in combination with efficient enhanced sampling molecular dynamics simulation methods. Analysis of conformational heterogeneity is performed based on sampling of discreet states as defined by dihedral angles, on root-mean-square differences of Cartesian coordinates and on the extent of volume sampled. Conformational clustering, based on the glycosidic linkage dihedral angles, shows that accounting for the full range of sampled conformations is required to reproduce the experimental data, emphasizing the utility of the molecular simulations in obtaining an atomic detailed description of the conformational properties of the saccharides. Results show the presence of differential conformational preferences as a function of primary sequence and glycosidic linkage types. Significant differences in conformational ensembles associated with the anomeric configuration of a single glycosidic linkage reinforce the impact of such changes on the conformational properties of carbohydrates. The present structural insights of the studied trisaccharides represent a foundation for understanding the range of conformations adopted in larger oligosaccharides and how these molecules encode their conformational heterogeneity into the monosaccharide sequence.

Enhanced conformational sampling using replica exchange with concurrent solute scaling and hamiltonian biasing realized in one dimension

Replica exchange (REX) is a powerful computational tool for overcoming the quasi-ergodic sampling problem of complex molecular systems. Recently, several multidimensional extensions of this method have been developed to realize exchanges in both temperature and biasing potential space or the use of multiple biasing potentials to improve sampling efficiency. However, increased computational cost due to the multidimensionality of exchanges becomes challenging for use on complex systems under explicit solvent conditions. In this study, we develop a one-dimensional (1D) REX algorithm to concurrently combine the advantages of overall enhanced sampling from Hamiltonian solute scaling and the specific enhancement of collective variables using Hamiltonian biasing potentials. In the present Hamiltonian replica exchange method, termed HREST-BP, Hamiltonian solute scaling is applied to the solute subsystem, and its interactions with the environment to enhance overall conformational transitions and biasing potentials are added along selected collective variables associated with specific conformational transitions, thereby balancing the sampling of different hierarchical degrees of freedom. The two enhanced sampling approaches are implemented concurrently allowing for the use of a small number of replicas (e.g., 6 to 8) in 1D, thus greatly reducing the computational cost in complex system simulations. The present method is applied to conformational sampling of two nitrogen-linked glycans (N-glycans) found on the HIV gp120 envelope protein. Considering the general importance of the conformational sampling problem, HREST-BP represents an efficient procedure for the study of complex saccharides, and, more generally, the method is anticipated to be of general utility for the conformational sampling in a wide range of macromolecular systems.

Preferred conformations of N-glycan core pentasaccharide in solution and in glycoproteins

N-linked glycans are on protein surfaces and have direct and water/ion-mediated interactions with surrounding amino acids. Such contacts could restrict their conformational freedom compared to the same glycans free in solution. In this work, we have examined the conformational freedom of the N-glycan core pentasaccharide moiety in solution using standard molecular dynamics (MD) simulations as well as temperature replica-exchange MD simulations. Both simulations yield the comparable conformational variability of the pentasaccharide in solution, indicating the convergence of both simulations. The glycoprotein crystal structures are analyzed to compare the conformational freedom of the N-glycan on the protein surface with the simulation result. Surprisingly, the pentasaccharide free in solution shows more restricted conformational variability than the N-glycan on the protein surface. The interactions between the carbohydrate and the protein side chain appear to be responsible for the increased conformational diversity of the N-glycan on the protein surface. Finally, the transfer entropy analysis of the simulation trajectory also reveals an unexpected causality relationship between intramolecular hydrogen bonds and the conformational states in that the hydrogen bonds play a role in maintaining the conformational states rather than driving the change in glycosidic torsional states.

Revealing the Mechanisms of Protein Disorder and N-Glycosylation in CD44-Hyaluronan Binding Using Molecular Simulation

The extracellular N-terminal hyaluronan binding domain (HABD) of CD44 is a small globular domain that confers hyaluronan (HA) binding functionality to this large transmembrane glycoprotein. When recombinantly expressed by itself, HABD exists as a globular water-soluble protein that retains the capacity to bind HA. This has enabled atomic-resolution structural biology experiments that have revealed the structure of HABD and its binding mode with oligomeric HA. Such experiments have also pointed to an order-to-disorder transition in HABD that is associated with HA binding. However, it had remained unclear how this structural transition was involved in binding since it occurs in a region of HABD distant from the HA-binding site. Furthermore, HABD is known to be N-glycosylated, and such glycosylation can diminish HA binding when the associated N-glycans are capped with sialic acid residues. The intrinsic flexibility of disordered proteins and of N-glycans makes it difficult to apply experimental structural biology approaches to probe the molecular mechanisms of how the order-to-disorder transition and N-glycosylation can modulate HA binding by HABD. We review recent results from molecular dynamics simulations that provide atomic-resolution mechanistic understanding of such modulation to help bridge gaps between existing experimental binding and structural biology data. Findings from these simulations include: Tyr42 may function as a molecular switch that converts the HA-binding site from a low affinity to a high affinity state; in the partially disordered form of HABD, basic amino acids in the C-terminal region can gain sufficient mobility to form direct contacts with bound HA to further stabilize binding; and terminal sialic acids on covalently attached N-glycans can form charge-paired hydrogen bonding interactions with basic amino acids that could otherwise bind to HA, thereby blocking HA binding to glycosylated CD44 HABD.

GS-align for glycan structure alignment and similarity measurement

MOTIVATION

Glycans play critical roles in many biological processes, and their structural diversity is key for specific protein-glycan recognition. Comparative structural studies of biological molecules provide useful insight into their biological relationships. However, most computational tools are designed for protein structure, and despite their importance, there is no currently available tool for comparing glycan structures in a sequence order- and size-independent manner.

RESULTS

A novel method, GS-align, is developed for glycan structure alignment and similarity measurement. GS-align generates possible alignments between two glycan structures through iterative maximum clique search and fragment superposition. The optimal alignment is then determined by the maximum structural similarity score, GS-score, which is size-independent. Benchmark tests against the Protein Data Bank (PDB) N-linked glycan library and PDB homologous/non-homologous N-glycoprotein sets indicate that GS-align is a robust computational tool to align glycan structures and quantify their structural similarity. GS-align is also applied to template-based glycan structure prediction and monosaccharide substitution matrix generation to illustrate its utility.

AVAILABILITY AND IMPLEMENTATION

http://www.glycanstructure.org/gsalign.

CONTACT

wonpil@ku.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Carbohydrate Microarrays

Carbohydrates, like nucleic acids and proteins, are essential biological molecules. Owing to their intrinsic physicochemical properties, carbohydrates are capable of generating structural diversity in a multitude of ways and are prominently displayed on the surfaces of cell membranes or on the exposed regions of macromolecules. Recent studies highlight that carbohydrate moieties are critical for molecular recognition, cell-cell interactions, and cell signaling in many physiological and pathological processes, and for biocommunication between microbes and host species. Modern carbohydrate microarrays emerged in 2002 and brought in new high-throughput tools for “glyco code” exploration. In this section, some basic concepts of sugar chain diversity, glyco-epitope recognition, and the evolving area of glyco-epitomics and biomarker discovery are discussed. Two complementary technologies, carbohydrate antigen arrays and photogenerated glyco-chips, serve as models to illustrate how to apply carbohydrate microarrays to address biomedical questions.

Terminal sialic acids on CD44 N-glycans can block hyaluronan binding by forming competing intramolecular contacts with arginine sidechains

Specific sugar residues and their linkages form the basis of molecular recognition for interactions of glycoproteins with other biomolecules. Seemingly small changes, like the addition of a single monosaccharide in the covalently attached glycan component of glycoproteins, can greatly affect these interactions. For instance, the sialic acid capping of glycans affects protein-ligand binding involved in cell-cell and cell-matrix interactions. CD44 is a single-pass transmembrane glycoprotein whose binding with its carbohydrate ligand hyaluronan (HA), an extracellular matrix component, mediates processes such as leukocyte homing, cell adhesion, and tumor metastasis. This binding is highly regulated by glycosylation of the N-terminal extracellular hyaluronan-binding domain (HABD); specifically, sialic acid capped N-glycans of HABD inhibit ligand binding. However, the molecular mechanism behind this sialic acid mediated regulation has remained unknown. Two of the five N-glycosyation sites of HABD have been previously identified as having the greatest inhibitory effect on HA binding, but only if the glycans contain terminal sialic acid residues. These two sites, Asn25 and Asn120, were chosen for in silico glycosylation in this study. Here, from extensive standard molecular dynamics simulations and biased simulations, we propose a molecular mechanism for this behavior based on spontaneously-formed charge-paired hydrogen bonding interactions between the negatively-charged sialic acid residues and positively-charged Arg sidechains known to be critically important for binding to HA, which itself is negatively charged. Such intramolecular hydrogen bonds would preclude associations critical to hyaluronan binding. This observation suggests how CD44 and related glycoprotein binding is regulated by sialylation as cellular environments fluctuate.

Effects of galactosylation in immunoglobulin G from all-atom molecular dynamics simulations

Explicit water atomistic molecular dynamics simulations (200 ns, ∼330,000 atoms) were performed to study the effects of galactosylation in the Fc domain of immunoglobulin G1. Two glycoforms were simulated to observe changes in protein-carbohydrate interactions and carbohydrate structure. A high degree of flexibility was observed in the small hinge region of the protein, while large domains remained stable. The hinge region flexibility allowed both translation and rotation of the domains relative to each other, resulting in a large number of possible conformations available. The distributions of rotational orientations between the Fab1 and Fab2 domains showed that while these domains are able to orient themselves rather freely pointing in space they rotated in unison to remain rotationally oriented at specific angles. Additionally, removing specific terminal galactose residues increased the mobility of the carbohydrate, resulting in different protein-carbohydrate interactions. Glycosylation has been suggested as a route to improve the aggregation resistance of monoclonal antibodies for therapeutic treatments to aid the immune system. The results presented here may provide insight into the search for IgG molecules with increased aggregation resistance to be used as monoclonal antibodies.

PTM-SD: a database of structurally resolved and annotated posttranslational modifications in proteins

Posttranslational modifications (PTMs) define covalent and chemical modifications of protein residues. They play important roles in modulating various biological functions. Current PTM databases contain important sequence annotations but do not provide informative 3D structural resource about these modifications. Posttranslational modification structural database (PTM-SD) provides access to structurally solved modified residues, which are experimentally annotated as PTMs. It combines different PTM information and annotation gathered from other databases, e.g. Protein DataBank for the protein structures and dbPTM and PTMCuration for fine sequence annotation. PTM-SD gives an accurate detection of PTMs in structural data. PTM-SD can be browsed by PDB id, UniProt accession number, organism and classic PTM annotation. Advanced queries can also be performed, i.e. detailed PTM annotations, amino acid type, secondary structure, SCOP class classification, PDB chain length and number of PTMs by chain. Statistics and analyses can be computed on a selected dataset of PTMs. Each PTM entry is detailed in a dedicated page with information on the protein sequence, local conformation with secondary structure and Protein Blocks. PTM-SD gives valuable information on observed PTMs in protein 3D structure, which is of great interest for studying sequence-structure- function relationships at the light of PTMs, and could provide insights for comparative modeling and PTM predictions protocols. Database URL: PTM-SD can be accessed at http://www.dsimb.inserm.fr/dsimb_tools/PTM-SD/.

'Naked' and hydrated conformers of the conserved core pentasaccharide of N-linked glycoproteins and its building blocks

N-glycosylation of eukaryotic proteins is widespread and vital to survival. The pentasaccharide unit -Man3GlcNAc2- lies at the protein-junction core of all oligosaccharides attached to asparagine side chains during this process. Although its absolute conservation implies an indispensable role, associated perhaps with its structure, its unbiased conformation and the potential modulating role of solvation are unknown; both have now been explored through a combination of synthesis, laser spectroscopy, and computation. The proximal -GlcNAc-GlcNAc- unit acts as a rigid rod, while the central, and unusual, -Man-β-1,4-GlcNAc- linkage is more flexible and is modulated by the distal Man-α-1,3- and Man-α-1,6- branching units. Solvation stiffens the 'rod' but leaves the distal residues flexible, through a β-Man pivot, ensuring anchored projection from the protein shell while allowing flexible interaction of the distal portion of N-glycosylation with bulk water and biomolecular assemblies.

Molecular mechanisms of inhibition of influenza by surfactant protein D revealed by large-scale molecular dynamics simulation

Surfactant protein D (SP-D), a mammalian C-type lectin, is the primary innate inhibitor of influenza A virus (IAV) in the lung. Interactions of SP-D with highly branched viral N-linked glycans on hemagglutinin (HA), an abundant IAV envelope protein and critical virulence factor, promote viral aggregation and neutralization through as yet unknown molecular mechanisms. Two truncated human SP-D forms, wild-type (WT) and double mutant D325A+R343V, representing neck and carbohydrate recognition domains are compared in this study. Whereas both WT and D325A+R343V bind to isolated glycosylated HA, WT does not inhibit IAV in neutralization assays; in contrast, D325A+R343V neutralization compares well with that of full-length native SP-D. To elucidate the mechanism for these biochemical observations, we have determined crystal structures of D325A+R343V in the presence and absence of a viral nonamannoside (Man9). On the basis of the D325A+R343V-Man9 structure and other crystallographic data, models of complexes between HA and WT or D325A+R343V were produced and subjected to molecular dynamics. Simulations reveal that whereas WT and D325A+R343V both block the sialic acid receptor site of HA, the D325A+R343V complex is more stable, with stronger binding caused by additional hydrogen bonds and hydrophobic interactions with HA residues. Furthermore, the blocking mechanism of HA differs for WT and D325A+R343V because of alternate glycan binding modes. The combined results suggest a mechanism through which the mode of SP-D-HA interaction could significantly influence viral aggregation and neutralization. These studies provide the first atomic-level molecular view of an innate host defense lectin inhibiting its viral glycoprotein target.