Modeling pyranose ring pucker in carbohydrates using machine learning and semi-empirical quantum chemical methods

A neural network potential based on pairwise resolved atomic forces and energies

Molecular simulations have become a key tool in molecular and materials design. Machine learning (ML)-based potential energy functions offer the prospect of simulating complex molecular systems efficiently at quantum chemical accuracy. In previous work, we have introduced the ML-based PairF-Net approach to neural network potentials, that adopts a pairwise interatomic scheme to predicting forces within a molecular system. Here, we further develop the PairF-Net model to intrinsically incorporate energy conservation and couple the model to a molecular mechanical (MM) environment within the OpenMM package. The updated PairF-Net model yields energy and force predictions and dynamical distributions in good agreement with the rMD17 dataset of ten small organic molecules in the gas-phase. We further show that these in vacuo ML models of small molecules can be applied to force predictions in aqueous solution via hybrid ML/MM simulations. We present a new benchmark dataset for these ten molecules in solution, obtained from QM/MM simulations, which we denote as rMD17-aq (https://zenodo.org/records/10048644); and assess the ability of PairF-Net to reproduce the molecular energy, atomic forces and dynamical distributions of these solution conformations via ML/MM simulations.

Reinvestigation of the internal glycan rearrangement of Lewis a and blood group type H1 epitopes

Protonated ions of fucose-containing oligosaccharides are prone to undergo internal glycan rearrangement which results in chimeric fragments that obfuscate mass-spectrometric analysis. Lack of accessible tools that would facilitate systematic analysis of glycans in the gas phase limits our understanding of this phenomenon. In this work, we use density functional theory modeling to interpret cryogenic IR spectra of Lewis a and blood group type H1 trisaccharides and to establish whether these trisaccharides undergo the rearrangement during gas-phase analysis. Structurally unconstrained search reveals that none of the parent ions constitute a thermodynamic global minimum. In contrast, predicted collision cross sections and anharmonic IR spectra provide a good match to available experimental data which allowed us to conclude that fucose migration does not occur in these antigens. By comparing the predicted structures with those obtained for Lewis x and blood group type H2 epitopes, we demonstrate that the availability of the mobile proton and a large difference in the relative stability of the parent ions and rearrangement products constitute the prerequisites for the rearrangement reaction.

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.