Induction of rare conformation of oligosaccharide by binding to calcium-dependent bacterial lectin: X-ray crystallography and modelling study

Comparative Assessment of Water Models in Protein-Glycan Interaction: Insights from Alchemical Free Energy Calculations and Molecular Dynamics Simulations

Accurate computational simulations of protein-glycan dynamics are crucial for a comprehensive understanding of critical biological mechanisms, including host-pathogen interactions, immune system defenses, and intercellular communication. The accuracy of these simulations, including molecular dynamics (MD) simulation and alchemical free energy calculations, critically relies on the appropriate parameters, including the water model, because of the extensive hydrogen bonding with glycan hydroxyl groups. However, a systematic evaluation of water models' accuracy in simulating protein-glycan interaction at the molecular level is still lacking. In this study, we used full atomistic MD simulations and alchemical absolute binding free energy (ABFE) calculations to investigate the performance of five distinct water models in six protein-glycan complex systems. We evaluated water models' impact on structural dynamics and binding affinity through over 5.8 μs of simulation time per system. Our results reveal that most protein-glycan complexes are stable in the overall structural dynamics regardless of the water model used, while some show obvious fluctuations with specific water models. More importantly, we discover that the stability of the binding motif's conformation is dependent on the water model chosen when its residues form weak hydrogen bonds with the glycan. The water model also influences the conformational stability of the glycan in its bound state according to density functional theory (DFT) calculations. Using alchemical ABFE calculations, we find that the OPC water model exhibits exceptional consistency with experimental binding affinity data, whereas commonly used models such as TIP3P are less accurate. The findings demonstrate how different water models affect protein-glycan interactions and the accuracy of binding affinity calculations, which is crucial in developing therapeutic strategies targeting these interactions.

Conformational Preferences at the Glycosidic Linkage of Saccharides in Solution as Deduced from NMR Experiments and MD Simulations: Comparison to Crystal Structures

Glycans are central to information content and regulation in biological systems. These carbohydrate molecules are active either as oligo- or polysaccharides, often in the form of glycoconjugates. The monosaccharide entities are joined by glycosidic linkages and stereochemical arrangements are of utmost importance in determining conformation and flexibility of saccharides. The conformational preferences and population distributions at the glycosidic torsion angles φ and ψ have been investigated for O-methyl glycosides of three disaccharides where the substitution takes place at a secondary alcohol, viz., in α-l-Fucp-(1→3)-β-d-Glcp-OMe, α-l-Fucp-(1→3)-α-d-Galp-OMe and α-d-Glcp-(1→4)-α-d-Galp-OMe, corresponding to disaccharide structural elements present in bacterial polysaccharides. Stereochemical differences at or adjacent to the glycosidic linkage were explored by solution state NMR spectroscopy using one-dimensional 1 H,1 H-NOESY NMR experiments to obtain transglycosidic proton-proton distances and one- and two-dimensional heteronuclear NMR experiments to obtain 3 JCH transglycosidic coupling constants related to torsion angles φ and ψ. Computed effective proton-proton distances from molecular dynamics (MD) simulations showed excellent agreement to experimentally derived distances for the α-(1→3)-linked disaccharides and revealed that for the bimodal distribution at the ψ torsion angle for the α-(1→4)-linked disaccharide experiment and simulation were at variance with each other, calling for further force field developments. The MD simulations disclosed a highly intricate inter-residue hydrogen bonding pattern for the α-(1→4)-linked disaccharide, including a nonconventional hydrogen bond between H5' in the glucosyl residue and O3 in the galactosyl residue, supported by a large downfield 1 H NMR chemical shift displacement compared to α-d-Glcp-OMe. Comparison of population distributions of the glycosidic torsion angles φ and ψ in the disaccharide entities to those of corresponding crystal structures highlighted the potential importance of solvation on the preferred conformation.

Mutation of two intronic nucleotides alters RNA structure and dynamics inhibiting MBNL1 and RBFOX1 regulated splicing of the Insulin Receptor

Alternative splicing (AS) of Exon 11 of the Insulin Receptor ( INSR ) is highly regulated and disrupted in several human disorders. To better understand INSR exon 11 AS regulation, splicing activity of an INSR exon 11 minigene reporter was measured across a gradient of the AS regulator muscleblind-like 1 protein (MBNL1). The RNA-binding protein Fox-1 (RBFOX1) was added to determine its impact on MBNL1-regulated splicing. The role of the RBFOX1 UGCAUG binding site within intron 11 was assessed across the MBNL1 gradient. Mutating the UGCAUG motif inhibited RBFOX1 regulation of exon 11 and had the unexpected effect of reducing MBNL1 regulation of this exon. Molecular dynamics simulations showed that exon 11 and the adjacent RNA adopts a dynamically stable conformation. Mutation of the RBFOX1 binding site altered RNA structure and dynamics, while a mutation that created an optimal MBNL1 binding site at the RBFOX1 site shifted the RNA back to wild type. An antisense oligonucleotide (ASO) was used to confirm the structure in this region of the pre-mRNA. This example of intronic mutations shifting pre-mRNA structure and dynamics to modulate splicing suggests RNA structure and dynamics should be taken into consideration for AS regulation and therapeutic interventions targeting pre-mRNA.

N-Acetyl Side-Chain Conformation in Saccharides: Solution Models Obtained from MA'AT Analysis

MA'AT analysis has been applied to model the conformational properties of N-acetyl side-chains in biologically important GlcNAc and ManNAc monosaccharides and in a βGlcNAc-(1→4)-βGlcNAc disaccharide. Density functional theory calculations were conducted to obtain parameterized equations that relate the magnitudes and signs of 10 spin-coupling constants to conformations of the C2-N2 bonds of GlcNAc and ManNAc. Six of these equations were used with experimental J-couplings, measured in H₂O/²H₂O and DMSO-d₆ solvents in selectively ¹³C-labeled compounds, to model the C1-C2-N2-C1' torsion angle (θ₁) in GlcNAc and ManNAc residues. MA'AT analysis gave mean values of θ₁ of 106° for αGlcNAc and ∼116° for βGlcNAc residues, with circular standard deviations (CSDs) of 21-22° in aqueous solution, in excellent agreement with those obtained by aqueous molecular dynamics (MD) simulation. Parameter space plots revealed unique MA'AT fits of the data, and root mean squared deviations (<0.2 Hz) were twofold smaller than those back-calculated from MD models, indicating that the MA'AT models better fit the experimental J-couplings. Context effects on both θ₁ values were found to be small in a βGlcNAc-(1→4)-βGlcNAc disaccharide. MA'AT analysis gave a mean value of θ₁ of 249° for αManNAc in H₂O/²H₂O, with a CSD of ∼19°, with both values in good agreement with MD. MA'AT models of N-acetyl side-chains are similar to those obtained previously for O-acetyl side-chains (J. Phys. Chem. B 2017, 121, 66-77). Both O- and N-acetylation conformationally constrain the C-O or C-N bonds relative to the same bonds in unsubstituted compounds. The present work confirms the ability of MA'AT analysis to reveal relatively small changes in mean molecular torsion angles in solution and provides additional evidence of the method as an experimental tool complementary to MD simulation.

Structure-function relationship of a novel fucoside-binding fruiting body lectin from Coprinopsis cinerea exhibiting nematotoxic activity

Lectins are non-immunoglobulin-type proteins that bind to specific carbohydrate epitopes and play important roles in intra- and inter-organismic interactions. Here, we describe a novel fucose-specific lectin, termed CML1, which we identified from fruiting body extracts of Coprinopsis cinerea. For further characterization, the coding sequence for CML1 was cloned and heterologously expressed in Escherichia coli. Feeding of CML1-producing bacteria inhibited larval development of the bacterivorous nematode Caenorhabditis tropicalis, but not of C. elegans. The crystal structure of the recombinant protein in its apo-form and in complex with H type I or Lewis A blood group antigens was determined by X-ray crystallography. The protein folds as a sandwich of 2 antiparallel β-sheets and forms hexamers resulting from a trimer of dimers. The hexameric arrangement was confirmed by small-angle X-ray scattering (SAXS). One carbohydrate-binding site per protomer was found at the dimer interface with both protomers contributing to ligand binding, resulting in a hexavalent lectin. In terms of lectin activity of recombinant CML1, substitution of the carbohydrate-interacting residues His54, Asn55, Trp94, and Arg114 by Ala abolished carbohydrate-binding and nematotoxicity. Although no similarities to any characterized lectin were found, sequence alignments identified many non-characterized agaricomycete proteins. These results suggest that CML1 is the founding member of a novel family of fucoside-binding lectins involved in the defense of agaricomycete fruiting bodies against predation by fungivorous nematodes.

Increasing the Affinity of an O-Antigen Polysaccharide Binding Site in Shigella flexneri Bacteriophage Sf6 Tailspike Protein

Broad and unspecific use of antibiotics accelerates spread of resistances. Sensitive and robust pathogen detection is thus important for a more targeted application. Bacteriophages contain a large repertoire of pathogen-binding proteins. These tailspike proteins (TSP) often bind surface glycans and represent a promising design platform for specific pathogen sensors. We analysed bacteriophage Sf6 TSP that recognizes the O-polysaccharide of dysentery-causing Shigella flexneri to develop variants with increased sensitivity for sensor applications. Ligand polyrhamnose backbone conformations were obtained from 2D 1 H,1 H-trNOESY NMR utilizing methine-methine and methine-methyl correlations. They agreed well with conformations obtained from molecular dynamics (MD), validating the method for further predictions. In a set of mutants, MD predicted ligand flexibilities that were in good correlation with binding strength as confirmed on immobilized S. flexneri O-polysaccharide (PS) with surface plasmon resonance. In silico approaches combined with rapid screening on PS surfaces hence provide valuable strategies for TSP-based pathogen sensor design.