A GROMOS Force Field for Furanose-Based Carbohydrates

Deciphering the conformational preference and ionization dynamics of tetrahydrofuran: Insights from advanced spectroscopic techniques

Tetrahydrofuran (THF) has garnered significant attention due to its pivotal role in biological and chemical processes. The diverse array of conformations exhibited by THF profoundly impacts its reactivity and interactions with other molecules. Understanding these conformational preferences is crucial for comprehending its molecular behavior. In this study, we utilize infrared (IR) resonant vacuum ultraviolet photoionization/mass-analyzed threshold ionization (VUV-PI/MATI) mass spectroscopies to capture distinctive vibrational spectra of individual conformers, namely, "twisted" and "bent," within THF. Our conformer-specific vibrational spectra provide valuable insights into the relative populations of these two conformers. The analysis reveals that the twisted (C2) conformer is more stable than the bent (CS) conformer by 17 ± 15 cm-1. By precisely tuning the VUV photon energy to coincide with vibrational excitation via IR absorption, we selectively ionize specific conformers, yielding two-photon IR + VUV-PI/MATI spectra corresponding to the twisted and bent conformers. This investigation conclusively affirms that both the twisted and bent conformers coexist in the neutral state, while only the twisted conformer exists in the cationic state. These findings not only bridge gaps in existing knowledge but also provide profound insights into the behavior of this pivotal molecule in the realms of biology and medicine.

Characteristic 1H NMR spectra of β-d-ribofuranosides and ribonucleosides: factors driving furanose ring conformations

A series of β-d-ribofuranosides and ribonucleosides fused with 2,3-O-isopropylidene ring was synthesized and studied in terms of their conformational preferences. Based on the 1H NMR spectra, DFT calculations, and X-ray analysis the E 0-like and E 4-like conformations adopted by these furanosides are identified. The 3 E-like and 2 E-like conformations are assigned to ribonucleosides without the 2,3-O-isopropylidene group. The studies are supported by analysis of the structural data of β-d-ribofuranosides and ribonucleosides deposited in the Cambridge Crystallographic Data Center (CCDC) database. Finally, the factors influencing the conformational preferences of the furanose ring with the β-d-ribo configuration are indicated. These are the unfavorable ecliptic orientation of the 2-OH and 3-OH groups, the 1,3-pseudodiaxial interaction of the aglycone and terminal hydroxymethyl group and the endo-anomeric effect. It is also proved that the exo-anomeric effect acts in β-d-ribofuranosides.

Extending the Martini 3 Coarse-Grained Force Field to Carbohydrates

Carbohydrates play an essential role in a large number of chemical and biochemical processes. High structural diversity and conformational heterogeneity make it problematic to link their measurable properties to molecular features. Molecular dynamics simulations carried out at the level of classical force fields are routinely applied to study the complex processes occurring in carbohydrate-containing systems, while the usefulness of such simulations relies on the accuracy of the underlying theoretical model. In this article, we present the coarse-grained force field dedicated to glucopyranose-based carbohydrates and compatible with the recent version of the Martini force field (v. 3.0). The parameterization was based on optimizing bonded and nonbonded parameters with a reference to the all-atom simulation results and the experimental data. Application of the newly developed coarse-grained carbohydrate model to oligosaccharides curdlan and cellulose displays spontaneous formation of aggregates of experimentally identified features. In contact with other biomolecules, the model is capable of recovering the protective effect of glucose monosaccharides on a lipid bilayer and correctly identifying the binding pockets in carbohydrate-binding proteins. The features of the newly proposed model make it an excellent candidate for further extensions, aimed at modeling more complex, functionalized, and biologically relevant carbohydrates.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Molecular simulations of complex carbohydrates and glycoconjugates

Complex carbohydrates (glycans) are the most abundant and versatile biopolymers in nature. The broad diversity of biochemical functions that carbohydrates cover is a direct consequence of the variety of 3D architectures they can adopt, displaying branched or linear arrangements, widely ranging in sizes, and with the highest diversity of building blocks of any other natural biopolymer. Despite this unparalleled complexity, a common denominator can be found in the glycans' inherent flexibility, which hinders experimental characterization, but that can be addressed by high-performance computing (HPC)-based molecular simulations. In this short review, I present and discuss the state-of-the-art of molecular simulations of complex carbohydrates and glycoconjugates, highlighting methodological strengths and weaknesses, important insights through emblematic case studies, and suggesting perspectives for future developments.

Anomeric Stereoauxiliary Cleavage of the C-N Bond of d-Glucosamine for the Preparation of Imidazo[1,5-a]pyridines

The targeted cleavage of the C-N bonds of alkyl primary amines in sustainable compounds of biomass according to a metal-free pathway and the conjunction of nitrogen in the synthesis of imidazo[1,5-a]pyridines are still highly challenging. Despite tremendous progress in the synthesis of imidazo[1,5-a]pyridines over the past decade, many of them can still not be efficiently prepared. Herein, we report an anomeric stereoauxiliary approach for the synthesis of a wide range of imidazo[1,5-a]pyridines after cleaving the C-N bond of d-glucosamine (α-2° amine) from biobased resources. This new approach expands the scope of readily accessible imidazo[1,5-a]pyridines relative to existing state-of-the-art methods. A key strategic advantage of this approach is that the α-anomer of d-glucosamine enables C-N bond cleavage via a seven-membered ring transition state. By using this novel method, a series of imidazo[1,5-a]pyridine derivatives (>80 examples) was synthesized from pyridine ketones (including para-dipyridine ketone) and aldehydes (including para-dialdehyde). Imidazo[1,5-a]pyridine derivatives containing diverse important deuterated C(sp2 )-H and C(sp3 )-H bonds were also efficiently achieved.

Carbohydrate Force Fields: The Role of Small Partial Atomic Charges in Preventing Conformational Collapse

Although the quality of current additive all-atom force fields for carbohydrates has been demonstrated in many applications, occasional significant differences reported for the hydrodynamic behavior of specific polysaccharides modeled with different force fields is a cause for concern. In particular, irreversible conformational collapse has been noted for some polysaccharide simulations with the GLYCAM06j force field. Here, we investigate the cause of this phenomenon through comparative simulations of a range of saccharides with both the GLYCAM06j and the CHARMM36 carbohydrate force fields. We find that conformational collapse in GLYCAM06j occurs for saccharide chains containing the deoxy sugar α-l-rhamnose after relatively long simulation intervals. Further, we explore the mechanism of conformational collapse and show that this phenomenon arises because of the anomalous low energy in GLYCAM06j (as compared to quantum mechanical calculations) of a specific orientation of α-l-Rha to α-l-Rha glycosidic linkages, which are subsequently sustained by intramolecular interactions in the saccharide chain. We identify the lack of partial charges on aliphatic hydrogens in GLYCAM as the source of this anomaly, demonstrating that addition of small partial atomic charges on the aliphatic protons in rhamnose removes the conformational collapse phenomenon. This work reveals the large cumulative impact that small partial charges may have on the dynamic behavior of polysaccharides and indicates that future reparameterization of the GLYCAM06j force field should investigate the addition of partial charges on all aliphatic hydrogens.

Low-temperature heat capacity and pseudorotation in 2 methyltetrahydrofuran

The heat capacity and phase transitions of 2-methyltetrahydrofuran in the temperature range from 7 to 350 K were measured by adiabatic calorimetry. The smoothed molar thermodynamic functions in the condensed...

Tautomers of N-acetyl-d-allosamine: an NMR and computational chemistry study

d-Allosamine is a rare sugar in Nature but its pyranoid form has been found α-linked in the core region of the lipopolysaccharide from the Gram-negative bacterium Porphyromonas gingivalis and in the chitanase inhibitor allosamidin, then β-linked and N-acetylated. In water solution the monosaccharide N-acetyl-d-allosamine (d-AllNAc) shows a significant presence of four tautomers arising from pyranoid and furanoid ring forms and anomeric configurations. The furanoid ring forms both showed ³J_H1,H2≈ 4.85 Hz and to differentiate the anomeric configurations a series of chemical shift anisotropy/dipole-dipole cross-correlated relaxation NMR experiments was performed in which the α-anomeric form showed notable different relaxation rates for its components of the H1 doublet, thereby making it possible to elucidate the anomeric configuration of each of the furanoses. The conformational preferences of the different forms of d-AllNAc were investigated by ³J_HH, ²J_CH and ³J_CH coupling constants from NMR experiments, molecular dynamics simulations and density functional theory calculations. The pyranose form resides in the ⁴C₁ conformation and the furanose ring form has the majority of its conformers located on the South-East region of the pseudorotation wheel, with a small population in the Northern hemisphere. The tautomeric equilibrium was quite sensitive to changes in temperature, where the β-anomer of the pyranoid ring form decreased upon a temperature increase while the other forms increased.

Comparison of Carbohydrate Force Fields in Molecular Dynamics Simulations of Protein-Carbohydrate Complexes

In this paper, we present the results of molecular dynamics simulations aimed at critical comparison of classical, biomolecular force fields (FFs) in the context of their capabilities to describe the structural and thermodynamic features of carbohydrate-protein interactions. We have considered the three main families of FFs (CHARMM, GROMOS, and GLYCAM/AMBER) by applying them to investigate the seven different carbohydrate-protein complexes. The results indicate that although the qualitative pattern of several structural descriptors (intermolecular hydrogen bonding, ligand dynamic location, etc.) is conserved among the compared FFs, there also exists a number of significant divergences (mainly the patterns of contacts between particular amino acid residues and bound carbohydrate). The carbohydrate-protein unbinding free energies also vary from one FF to another, displaying diversified trends in deviations from the experimental data. The magnitude of those deviations is not negligible and indicates the need for refinement in the currently existing combinations of carbohydrate- and protein-dedicated biomolecular force fields. In spite of the lack of explicit functional terms responsible for the corresponding intermolecular forces, all tested FFs are capable of adequately reproducing the CH-π interactions, crucial for carbohydrate-protein binding.

Towards a Better Understanding of Computational Models for Predicting DNA Methylation Effects at the Molecular Level

Human DNA is a very sensitive macromolecule and slight changes in the structure of DNA can have disastrous effects on the organism. When nucleotides are modified, or changed, the resulting DNA sequence can lose its information, if it is part of a gene, or it can become a problem for replication and repair. Human cells can regulate themselves by using a process known as DNA methylation. This methylation is vitally important in cell differentiation and expression of genes. When the methylation is uncontrolled, however, or does not occur in the right place, serious pathophysiological consequences may result. Excess methylation causes changes in the conformation of the DNA double helix. The secondary structure of DNA is highly dependent upon the sequence. Therefore, if the sequence changes slightly the secondary structure can change as well. These slight changes will then cause the doublestranded DNA to be more open and available in some places where large adductions can come in and react with the DNA base pairs. Computer models have been used to simulate a variety of biological processes including protein function and binding, and there is a growing body of evidence that in silico methods can shed light on DNA methylation. Understanding the anomeric effect that contributes to the structural and conformational flexibility of furanose rings through a combination of quantum mechanical and experimental studies is critical for successful molecular dynamic simulations.

Three-Dimensional Structures of Carbohydrates and Where to Find Them

Analysis and systematization of accumulated data on carbohydrate structural diversity is a subject of great interest for structural glycobiology. Despite being a challenging task, development of computational methods for efficient treatment and management of spatial (3D) structural features of carbohydrates breaks new ground in modern glycoscience. This review is dedicated to approaches of chemo- and glyco-informatics towards 3D structural data generation, deposition and processing in regard to carbohydrates and their derivatives. Databases, molecular modeling and experimental data validation services, and structure visualization facilities developed for last five years are reviewed.

Chirality Effects in Biomolecular Systems: Calculation of the Relative Free Energies by Molecular Dynamics Simulations

Chirality plays an essential role in chemical and biological sciences. At the molecular level, the effects associated with this phenomenon can be studied by using the well-established technique of molecular dynamics simulations. In this work, we present several approaches suited for the molecular dynamics-based free energy calculation in chiral systems. In particular, we have proposed and tested the following strategies relying on the application of general, enhanced sampling methods: (i) biased sampling in the two-dimensional space, along the coordinates defined by the values of the selected torsional angles; (ii) biased sampling in the one- or two-dimensional space, along the path-based coordinate(s); (iii) rational alteration of the system's Hamiltonian in order to enable the interconversion between stereoisomers and reweighting the biased distribution of configurations; (iv) using the free energy landscape generated within approaches (i) or (ii) as time-independent bias in order to further improve sampling efficiency and simultaneously account for multiple chiral centers. All approaches have been tested on a set of model compounds (fenoterol, fructofuranose, and bromochlorofluoromethane), demonstrating the good performance but also some differences in the range of their applicabilities.

ADD Force Field for Sugars and Polyols: Predicting the Additivity of Protein-Osmolyte Interaction

The protein–osmolyte interaction has been shown experimentally to follow an additive construct, where the individual osmolyte–backbone and osmolyte–side-chain interactions contribute to the overall conformational stability of proteins. Here, we computationally reconstruct this additive relation using molecular dynamics simulations, focusing on sugars and polyols, including sucrose and sorbitol, as model osmolytes. A new set of parameters (ADD) is developed for this purpose, using the individual Kirkwood–Buff integrals for sugar–backbone and sugar–side-chain interactions as target experimental data. We show that the ADD parameters can reproduce the additivity of protein–sugar interactions and correctly predict sucrose and sorbitol self-association and their interaction with water. The accurate description of the separate osmolyte–backbone and osmolyte–side-chain contributions also automatically translates into a good prediction of preferential exclusion from the surface of ribonuclease A and α-chymotrypsinogen A. The description of sugar polarity is improved compared to previous force fields, resulting in closer agreement with the experimental data and better compatibility with charged groups, such as the guanidinium moiety. The ADD parameters are developed in combination with the CHARMM36m force field for proteins, but good compatibility is also observed with the AMBER 99SB-ILDN and the OPLS-AA force fields. Overall, exploiting the additivity of protein–osmolyte interactions is a promising approach for the development of new force fields.

Conformational properties of inulin, levan and arabinan studied by molecular dynamics simulations

Inulin, levan and arabinan are the polysaccharides that consist of exclusively furanose units. To date, their conformational features studied at the molecular scale have remained largely unexplained. To tackle this issue, we have performed a series of explicit-solvent molecular dynamics simulations, carried out within the furanose-dedicated force field. None of the polysaccharides exhibits a single, dominating structure type. Instead, they create a large number of separated conformational states originating from the intensive rotation around the φ and ω glycosidic angles. 2₁-helices are the preferential conformational forms for all compounds but they appear only locally, at the length of several consecutive residues. The flexibility of all three furanose-based polysaccharides is much greater in relation to the (1-4)-linked pyranose polysaccharides and is comparable to that of (1-6)-linked pyranoses. The dynamic geometries of both furanose rings and glycosidic linkages are nearly unchanged independently if considering them at the level of mono-, di- or polysaccharides.

Atomic partial charge predictions for furanoses by random forest regression with atom type symmetry function

Furanoses that are components for many important biomolecules have complicated conformational spaces due to the flexible ring and exo-cyclic moieties. Machine learning algorithms, which require descriptors as structural inputs, can be used to efficiently compute conformational adaptive (CA) charges to capture the electrostatic potential variations caused by the conformational changes in the molecular mechanics (MM) calculations. In the present study, we introduced atom type symmetry function (ATSF) developed based on atom centered symmetry function (ACSF) for describing conformations for furanoses, in which atoms were categorized by atom types defined by their properties or connectivity in classic molecular mechanics (MM) force field parameters to generate a suitable coordinate size. Random forest regression (RFR) models with ATSF showed improvements for predicting CA charges and dipole moments for furanoses compared to those with ACSF and atom name symmetry functions where atoms were categorized by their unique atom names. The CA charges predicted by RFR models with ATSF showed more comparable reproductions of the carbohydrate-water and carbohydrate-protein interactions computed with RESP charges individually derived from QM calculations than the ensemble-averaged atomic charge sets commonly employed in molecular mechanics force fields, suggesting that the predicted CA charges were capable of including electrostatic variations in their dynamic charge values. Improvements by ATSF showed that categorizing atoms by atom types introduced chemical structural perceptions to descriptors and produced a suitable coordinate size in ATSF to capture key structural features for furanoses. This categorizing scheme also allows ATSF to be readily adopted by other biomolecules thanks to the broad implementations of MM force fields.

Efficient sampling of high-energy states by machine learning force fields

A method extending the range of applicability of machine-learning force fields is proposed. It relies on biased subsampling of the high-energy states described by the predefined coordinate(s).

Relation between the NMR data and the pseudorotational free-energy profile for oxolane

The conformation of five-membered furanose rings is a crucial issue for the structural analysis of many biologically-relevant molecules, including DNA and RNA. Oxolane can be treated as a prototypical furanose, composed only of saturated unsubstituted ring. In spite of its structural simplicity, providing the accurate quantitative description of the oxolane conformational features remains a great challenge for both the experimental and theoretical techniques. Here we show the method of recovering the free-energy profiles describing the conformational equilibrium in the oxolane ring (i.e. pseudorotation) based on the experimentally-inferred NMR data ([Formula: see text] coupling constants). The results remain in agreement with the quantum-mechanical-based molecular dynamics simulations and emphasize the large contributions of all ring conformers, even those located at the free-energy barriers. This includes the significant populations of limiting 3T2/2T3 and OE/EO shapes. Our findings provide another example of a poor applicability of the two-state model, which is routinely applied to analyze the NMR data in terms of population of different ring conformers.