A first-principles exploration of the conformational space of sodiated pyranose assisted by neural network potentials

Sampling the conformational space of monosaccharides using the first-principles methods is important and as a database of local minima provides a solid base for interpreting experimental measurements such as infrared photo-dissociation (IRPD) spectroscopy or collision-induced dissociation (CID). IRPD emphasizes low-energy conformers and CID can distinguish conformers with distinct reaction pathways. A typical computational approach is to engage empirical or semi-empirical methods to sample the conformational space first, and only selected minima are reoptimized at first-principles levels. In this work, we propose a computational scheme to explore the configurational space of 12 types of sodiated pyranoses with the assistance of a neural network potential (NNP). We demonstrated that it is possible to train an NNP based on the density functional calculations extracted from a previous study on sodiated glucose (Glc), galactose (Gal), and mannose (Man). This NNP yields a better description of the other five types of aldohexoses than the four types of ketohexoses. We further show that such a discrepancy in the accuracy of NNP can be resolved by an active learning scheme where the NNP model is engaged in generating the data and has itself updated. Through this iterative process, we can locate more than 17 000 distinct local minima at the B3LYP/6-311+G(d,p) level and an NNP with an accuracy of 1 kJ mol^-1 was created, which can be used for further studies.

A self-adapting first-principles exploration on the dissociation mechanism in sodiated aldohexose pyranoses assisted with neural network potentials

Understanding the mechanism of collision-induced dissociation (CID) in mono-saccharides with density functional theory (DFT) is challenging because of many possible reaction paths that originate from their high structural diversity. To search for the transition state (TS) from the huge number of conformers, we propose a three-step search scheme with the assistance of neural network potential (NNP). The search starts from a cross-checking of sugars, to a global search of all possible channels, and in the end, an exhaustive exploration around the low-lying channels. The cross-checking step quickly adapts the NNP from the studied molecules to the target ones. The other two steps utilize the adapted NNP to find the available pathways via random sampling of the structures. The study of the CID reactions in all eight types of aldohexose pyranoses was applied using the search scheme. The DFT calculations on AH-0 (Glc, Gal, and Man) in the previous study were utilized to construct an NNP and provide the TS structure database for searching AH-1 (All, Alt, Gul, Ido, and Tal). In total, we identified around 5200 TSs in AH-0 and AH-1, and the final NNP covers an energy range of more than 500 kJ mol^-1 with a mean absolute error of energy less than 4 kJ mol^-1. The search scheme is useful not only for saccharides but also for highly flexible bio-molecules.

Protonated α-N-Acetyl Galactose Glycopeptide Dissociation Chemistry

We recently provided mass spectrometric, H/D labeling, and computational evidence of pyranose to furanose N-acetylated ion isomerization reactions that occurred prior to glycosidic bond cleavage in both O- and N-linked glycosylated amino acid model systems (Guan et al. Phys. Chem. Chem. Phys., 2021, 23, 23256-23266). These reactions occurred irrespective of the glycosidic linkage stereochemistry (α or β) and the N-acetylated hexose structure (GlcNAc or GalNAc). In the present article, we test the generality of the preceding findings by examining threonyl α-GalNAc-glycosylated peptides. We utilize computational chemistry to compare the various dissociation and isomerization pathways accessible with collisional activation. We then interrogate the structure(s) of the resulting charged glycan and peptide fragments with infrared "action" spectroscopy. Isomerization of the original pyranose, the protonated glycopeptide [AT(GalNAc)A+H]⁺, is predicted to be facile compared to direct dissociation, as is the glycosidic bond cleavage of the newly formed furanose form, i.e., furanose oxazolinium ion structures are predicted to predominate. IR action spectra for the m/z 204, C₈H₁₄N₁O₅⁺, glycan fragment population support this prediction. The IR action spectra of the complementary m/z 262 peptide fragment were assigned as a mixture of the lowest-energy structures of [ATA+H]⁺ consistent with the literature. If general, the change to a furanose m/z 204 product ion structure fundamentally alters the ion population available for MS³ dissociation and glycopeptide sequence identification.

Collision-induced dissociation of Na+-tagged ketohexoses: experimental and computational studies on fructose

Collision-induced dissociation tandem mass spectrometry (CID-MSⁿ) and computational investigation at the MP2/6-311+G(d,p) level of theory have been employed to study Na⁺-tagged fructose, an example of a ketohexose featuring four cyclic isomers: α-fructofuranose (αFruf), β-fructofuranose (βFruf), α-fructopyranose (αFrup), and β-fructopyranose (βFrup). The four isomers can be separated by high-performance liquid chromatography (HPLC) and they show different mass spectra, indicating that CID-MSⁿ can distinguish the different fructose forms. Based on a simulation using a micro-kinetic model, we have obtained an overview of the mechanisms for the different dissociation pathways. It has been demonstrated that the preference for the C-C cleavage over the competing isomerization of linear fructose is the main reason for the previously reported differences between the CID-MS spectra of aldohexoses and ketohexoses. In addition, the kinetic modeling helped to confirm the assignment of the different measured mass spectra to the different fructose isomers. The previously reported assignment based on the peak intensities in the HPLC chromatogram had left some open questions as the preference for the dehydration channels did not always follow trends previously observed for aldohexoses. Setting up the kinetic model further enabled us to directly compare the computational and experimental results, which indicated that the model can reproduce most trends in the differences between the dissociation pathways of the four cyclic fructose isomers.

Identification of Anomericity and Linkage of Arabinose and Ribose through Collision-Induced Dissociation

Arabinose and ribose are two common pentoses that exist in both furanose and pyranose forms in plant and bacteria oligosaccharides. In this study, each pentose isomer, namely α-furanose, β-furanose, α-pyranose, and β-pyranose, was first separated through high-performance liquid chromatography followed by an investigation of collision-induced dissociation in an ion trap mass spectrometer. The major dissociation channels, dehydration and cross-ring dissociation, were analyzed by using high-level quantum chemistry calculations and transition state theory. The branching ratio of major dissociation channels was governed by two geometrical features: one being the cis or trans configuration of O1 and O2 atoms determining dehydration preferability and the other being the number of hydroxyl groups on the same side of the ring as the O1 atom determining the preferability of cross-ring dissociation. The relative branching ratios of the major channels were used to identify anomericity and the linkages of arabinose and ribose. Arabinose in the β-configuration and ribose in the α-configuration are predicted to have larger relative dehydration branching ratios than arabinose in the α-configuration and ribose in the β-configuration, respectively. Arabinose and ribose at the reducing end of oligosaccharides with 1 → 2 (pyranose and furanose), 1 → 3 (pyranose and furanose), 1 → 4 (pyranose only), and 1 → 5 (furanose only) linkages are predicted to undergo ^0,2X, ^0,3X, ^0,2A, and ^0,2A/^0,3A cross-ring dissociation, respectively. Application of the dissociation mechanism to the disaccharide linkage determination is demonstrated.