Abstract:The conformational features of the glycosidic linkage are the most important variable to consider when studying di-, oligo-, and polysaccharide molecules using molecular dynamics (MD) simulations. The accuracy of the theoretical model describing this degree of freedom influences the quality of the results obtained from MD calculations based on this model. This article focuses on the following two issues related to the conformation of the glycosidic linkage. First, we describe the results of a comparative analysis of the predictions of three carbohydrate-dedicated classical force fields for MD simulations, namely, CHARMM, GLYCAM, and GROMOS, in the context of different parameters of structural and energetic nature related to the conformation of selected types of glycosidic linkages, α(1 → 4), β(1 → 3), and β(1 → 4), connecting glucopyranose units. This analysis revealed several differences, mainly concerning the energy levels of the secondary and tertiary conformers and the linkage flexibility within the dominant exo-syn conformation for α(1 → 4) and β(1 → 3) linkages. Some aspects of the comparative analysis also included the newly developed, carbohydrate-dedicated Martini 3 coarse-grained force field. Second, to overcome the time-scale problem associated with sampling slow degrees of freedom in polysaccharide chains during MD simulations, we developed a coarse-grained (CG) model based on the data from MD simulations and designed for Monte Carlo modeling. This model (CG MC) is based on information from simulations of short saccharide chains, effectively sampled in atomistic MD simulations, and is capable of extrapolating local conformational properties to the case of polysaccharides of arbitrary length. The CG MC model has the potential to estimate the conformations of very long polysaccharide chains, taking into account the influence of secondary and tertiary conformations of glycosidic linkages. With respect to the comparative analysis of force fields, the application of CG MC modeling showed that relatively small differences in the predictions of individual force fields with respect to a single glycosidic linkage accumulate when considering their effect on the structure of longer chains, leading to drastically different predictions with respect to parameters describing the polymer conformation, such as the persistence length.

Molecular Modeling of Glycosylated Catalytic Domain of Human Carbonic Anhydrase IX

Network of Conformational Transitions Revealed by Molecular Dynamics Simulations of the Carbonic Anhydrase II Apo-Enzyme

Ligand based pharmacophore modelling and integrated computational approaches in the quest for small molecule inhibitors against hCA IX

Comparative study of stability and activity of wild-type and mutant human carbonic anhydrase II enzymes using molecular dynamics and docking simulations

Control of Substrate Conformation by Hydrogen Bonding in a Retaining β-Endoglycosidase

Modelling and Molecular Dynamics Predict the Structure and Interactions of the Glycine Receptor Intracellular Domain

Coarse-Grained Model of Glycosaminoglycans for Biomolecular Simulations

Identification of novel CA IX inhibitor: Pharmacophore modeling, docking, DFT, and dynamic simulation

Computational modeling of novel inhibitory peptides targeting proteoglycan like region of carbonic anhydrase IX and in vitro validation in HeLa cells

Markov-state Model for CO2 Binding with Carbonic Anhydrase under Confinement

A combination strategy of structure-based virtual screening, MM-GBSA, cross docking, molecular dynamics and metadynamics simulations used to investigate natural compounds as potent and specific inhibitors of tumor linked human carbonic anhydrase IX

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

Privileged Scaffold Hybridization in the Design of Carbonic Anhydrase Inhibitors

Discovery of non-sulfonamide carbonic anhydrase IX inhibitors through structure-based virtual screening

Large Conformational Changes of Insertion 3 in Human Glycyl-tRNA Synthetase (hglyrs) During Catalysis*

A computational approach for deciphering the organization of glycosaminoglycans

1H, 15N and 13C resonance assignments of the S2A and H64A double mutant of human carbonic anhydrase II

Comparative Bindings of Glycosaminoglycans with CXCL8 Monomer and Dimer: Insights from Conformational Dynamics and Kinetics of Hydrogen Bonds

A Case Study of the Glycoside Hydrolase Enzyme Mechanism Using an Automated QM-Cluster Model Building Toolkit

Detachment of HCO3– from the Active Site of Carbonic Anhydrase: Molecular Dynamics Simulation and Machine Learning

The Conformation of Glycosidic Linkages According to Various Force Fields: Monte Carlo Modeling of Polysaccharides Based on Extrapolation of Short-Chain Properties