Glycosidic α-linked mannopyranose disaccharides: an NMR spectroscopy and molecular dynamics simulation study employing additive and Drude polarizable force fields

NMR investigations of glycan conformation, dynamics, and interactions

Glycans are ubiquitous in nature, decorating our cells and serving as the initial points of contact with any visiting entities. These glycan interactions are fundamental to host-pathogen recognition and are related to various diseases, including inflammation and cancer. Therefore, understanding the conformations and dynamics of glycans, as well as the key features that regulate their interactions with proteins, is crucial for designing new therapeutics. Due to the intrinsic flexibility of glycans, NMR is an essential tool for unravelling these properties. In this review, we describe the key NMR parameters that can be extracted from the different experiments, and which allow us to deduce the necessary geometry and molecular motion information, with a special emphasis on assessing the internal motions of the glycosidic linkages. We specifically address the NMR peculiarities of various natural glycans, from histo-blood group antigens to glycosaminoglycans, and also consider the special characteristics of their synthetic analogues (glycomimetics). Finally, we discuss the application of NMR protocols to study glycan-related molecular recognition events, both from the carbohydrate and receptor perspectives, including the use of stable isotopes and paramagnetic NMR methods to overcome the inherent degeneracy of glycan chemical shifts.

Refinement of the Drude Polarizable Force Field for Hexose Monosaccharides: Capturing Ring Conformational Dynamics with Enhanced Accuracy

We present a revised version of the Drude polarizable carbohydrate force field (FF), focusing on refining the ring and exocyclic torsional parameters for hexopyranose monosaccharides. This refinement addresses the previously observed discrepancies between calculated and experimental NMR 3J coupling values, particularly in describing ring dynamics and exocyclic rotamer populations within major hexose monosaccharides and their anomers. Specifically, α-MAN, β-MAN, α-GLC, β-GLC, α-GAL, β-GAL, α-ALT, β-ALT, α-IDO, and β-IDO were targeted for optimization. The optimization process involved potential energy scans (PES) of the ring and exocyclic dihedral angles computed using quantum mechanical (QM) methods. The target data for the reoptimization included PES of the inner ring dihedrals (C1-C2-C3-C4, C2-C3-C4-C5, C5-O5-C1-C2, C4-C5-O5-C1, O5-C1-C2-C3, C3-C4-C5-O5) and the exocyclic torsions, other than the pseudo ring dihedrals (O1-C1-O5-C5, O2-C2-C1-O5, and O4-C4-C5-O5) and hydroxyl torsions used in the previous parametrization efforts. These parameters, in conjunction with previously developed Drude parameters for hexopyranose monosaccharides, were validated against experimental observations, including NMR data and conformational energetics, in aqueous environments. The resulting polarizable model is shown to be in good agreement with a range of QM data, experimental NMR data, and conformational energetics of monosaccharides in aqueous solutions. This offers a significant improvement of the Drude carbohydrate force field, wherein the refinement enhances the accuracy of accessing the conformational dynamics of carbohydrates in biomolecular simulations.

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

NMR Coupling Constants, Karplus Equations, and Adjusted MD Statistics: Detecting Diagnostic Torsion Angles for the Solution Geometry of 6-[α-d-Mannopyranosyl]-d-Mannopyranose (Mannobiose)

The quantitative solution conformations of 2-(hydroxymethyl)-tetrahydropyran, α-methyl-d-mannopyranoside, and 6-[α-d-mannopyranosyl]-d-mannopyranose (mannobiose) are described. Parametrized Karplus equations for redundant spin pairs across the terminal ω-torsion and the glycosidic ω-torsion for mannobiose are developed, including ω/θ-hypersurfaces for the terminal hydroxymethylene group. Experimental NMR data, algorithmic spectral simulation (clustered Hamiltonian method), molecular dynamics (MD) simulations (GLYCAM06), energy minimizations by DFT, and adjusted torsion angle populations weighted over the Karplus-type equations are used. We demonstrate that spectral simulation is a powerful tool in the refinement of initial J values obtained from static GAIO DFT calculations. We also show that only as few as one of multiple redundant torsions can be diagnostic for conformational analysis of the disaccharide.

Targeting the PANoptosis signaling pathway for myocardial protection: therapeutic potential of Xian Ling Gu Bao capsule

Introduction: Myocardial infarction (MI), the most prevalent ischemic heart disease, constitutes a primary cause of global cardiovascular disease with incidence and mortality. The pathogenesis of MI is exceedingly intricate, with PANoptosis playing a pivotal role in its pathological process. Xian Ling Gu Bao capsule (XLGB) contains various active components, including flavonoids, terpenes, and phenylpropanoids, and exhibits a wide range of pharmacological activities. However, it remains unclear whether XLGB can protect the myocardium from damage after MI. This study aimed to investigate the impact of XLGB on isoprenaline (ISO)-induced MI in mice and its potential mechanisms. Methods: This study assessed the protective effects of XLGB against ISO-induced MI through techniques such as echocardiography, HE staining, Masson staining, and enzyme-linked immunosorbent assay (ELISA). Furthermore, the potential mechanisms of XLGB's protective effects on MI were explored using bioinformatics, molecular docking, and molecular dynamics simulations. These mechanisms were further validated through immunofluorescence staining and Western blotting. Results: The results demonstrated that various doses of XLGB exhibited a significant reduction in myocardial injury induced by myocardial infarction. Intriguingly, higher dosages of XLGB displayed superior therapeutic efficacy compared to the positive control metoprolol. This protective effect is primarily achieved through the inhibition of oxidative stress and the inflammatory processes. Furthermore, we have elucidated that XLGB protected the myocardium from MI-induced damage by suppressing PANoptosis, with a critical role played by the NLRP3/Caspase3/RIP1 signaling pathway. Of particular note, the primary compounds of XLGB were found to directly interact with NLRP3/Caspase3/RIP1, a discovery further validated through molecular docking and molecular dynamics simulations. This suggests that NLRP3/Caspase3/RIP1 may be a therapeutic target for XLGB-induced myocardial protection. Conclusion: In summary, our findings reveal a novel property of XLGB: reverses myocardial damage following MI by inhibiting the NLRP3/Caspase3/RIP1-mediated PANoptosis pathway.

Insights into the transcriptomic mechanism and characterization of endoglucanases from Aspergillus terreus in cellulose degradation

Filamentous fungi are the main industrial source of cellulases which are important in the process of converting cellulose to fermentable sugars. In this study, transcriptome analysis was conducted on Aspergillus terreus NEAU-7 cultivated using corn stover and glucose as carbon sources. Four putative endoglucanases (EG5A, EG7A, EG12A, and EG12C) from A. terreus NEAU-7 were efficiently expressed in Pichia pastoris. Among them, EG7A exhibited the highest enzyme activity (75.17 U/mg) with an optimal temperature of 40 °C and pH 5.0. EG5A and EG12A displayed specific activities of 19.92 U/mg and 14.62 U/mg, respectively, at 50 °C. EG12C showed acidophilic characteristics with an optimal pH of 3.0 and a specific activity of 12.21 U/mg at 40 °C. With CMC-Na as the substrate, the Km value of EG5A, EG7A, EG12A or, EG12C was, 11.08 ± 0.87 mg/mL, 6.82 ± 0.74 mg/mL, 7.26 ± 0.64 mg/mL, and 9.88 ± 0.86 mg/mL, with Vmax values of 1258.23 ± 51.62 μmol∙min-1∙mg-1, 842.65 ± 41.53 μmol∙min-1∙mg-1, 499.38 ± 20.42 μmol∙min-1∙mg-1, and 681.41 ± 30.08 μmol∙min-1∙mg-1, respectively. The co-treatment of EG7A with the commercial cellulase increased the yield of reducing sugar by 155.77 % (filter paper) and 130.49 % (corn stover). Molecular docking assay showed the interaction energy of EG7A with cellotetraose at -10.50 kcal/mol, surpassing EG12A (-10.43 kcal/mol), EG12C (-10.28 kcal/mol), and EG5A (-9.00 kcal/mol). Root Mean Square Deviation (RMSD) and Solvent Accessible Surface Area (SASA) values revealed that the presence of cellotetraose stabilized the molecular dynamics simulation of the cellotetraose-protein complex over a 100 ns time scale. This study provides valuable insights for developing recombinant enzymes and biomass degradation technologies.

Does size matter? - Comparing pyranoses with septanoses as ligands of the bacterial lectin FimH

The pharmacological modulation of disease-relevant carbohydrate-protein interactions represents an underexplored area of medicinal chemistry. One particular challenge in the design of glycomimetic compounds is the inherent instability of the glycosidic bond toward enzymatic cleavage. This problem has traditionally been approached by employing S-, N-, or C-glycosides with reduced susceptibility toward glycosidases. The application of ring-extended glycomimetics is an innovative approach to circumvent this issue. On the example of the bacterial adhesin FimH, it was explored how design principles from pyranose glycomimetics transfer to analogous septanose structures. A series of ring-extended FimH antagonists exhibiting the well-proven pharmacophore necessary for targeting the tyrosine-gate of FimH was synthesized. The resulting septanoses were evaluated for their affinity to the conformationally rigid isolated lectin domain of FimH (FimHLD), as well as a structurally flexible full-length FimH (FimHFL) construct. Some elements of potent mannoside-based FimH antagonists could be successfully transferred to septanose-based ligands, ultimately resulting in a 32-fold increase in binding affinity. Interestingly, the canonical ca. 100-fold loss of binding affinity between FimHLD and FimHFL is partly mitigated by the more flexible septanose antagonists, hinting at potentially differing interaction features of the flexible glycomimetics with intermediately populated states during the conformational transition of FimHFL.

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.

Glycan Shape, Motions, and Interactions Explored by NMR Spectroscopy

Glycans in the form of oligosaccharides, polysaccharides, and glycoconjugates are ubiquitous in nature, and their structures range from linear assemblies to highly branched and decorated constructs. Solution state NMR spectroscopy facilitates elucidation of preferred conformations and shapes of the saccharides, motions, and dynamic aspects related to processes over time as well as the study of transient interactions with proteins. Identification of intermolecular networks at the atomic level of detail in recognition events by carbohydrate-binding proteins known as lectins, unraveling interactions with antibodies, and revealing substrate scope and action of glycosyl transferases employed for synthesis of oligo- and polysaccharides may efficiently be analyzed by NMR spectroscopy. By utilizing NMR active nuclei present in glycans and derivatives thereof, including isotopically enriched compounds, highly detailed information can be obtained by the experiments. Subsequent analysis may be aided by quantum chemical calculations of NMR parameters, machine learning-based methodologies and artificial intelligence. Interpretation of the results from NMR experiments can be complemented by extensive molecular dynamics simulations to obtain three-dimensional dynamic models, thereby clarifying molecular recognition processes involving the glycans.

Chitosan Polysaccharides from a Polarizable Multiscale Approach

We report simulations of chitosan polysaccharides in the aqueous phase, at infinite dilute conditions and zero ionic strength. Those simulations are performed by means of a polarizable multiscale modeling scheme that relies on a polarizable all atom force field to model solutes and on a polarizable solvent coarse grained approach. Force field parameters are assigned only from quantum chemistry ab initio data. We simulate chitosan monomer units, dimers and 50-long chains. Regarding the 50-long chains we simulate three sets of ten randomly built chain replica at three different pH conditions (corresponding to different chain protonation states, the chain degree of deacetylation is 85%). Our simulations show the persistence length of 50-long chitosan chains at strong acidic conditions (pH <5) to be 24 ± 2 nm (at weak/negligible ionic strength conditions), and to be 1 order of magnitude shorter at usual pH conditions. Our simulation data support the most recent simulation and experimental studies devoted to chitosan polysaccharides in the aqueous phase.