Abstract
Density functional theory (DFT) has provided detailed information on the molecular structure and spin⁻spin coupling constants of heparin tetrasaccharide (GlcNS,6S-IdoA2S-GlcNS,6S-IdoA2S-OMe) representing the predominant heparin repeating-sequence. The fully optimised molecular structures of two tetrasaccharide conformations (differing from each other in the conformational form of the sulphated iduronic acid residue⁻one ¹C₄ and the other ²S₀) were obtained using the B3LYP/6-311+G(d,p) level of theory and applying explicit water molecules to simulate the presence of a solvent. The theoretical data provided insight into variations of the bond lengths, bond angles and torsion angles, formations of intra- and intermolecular hydrogen bonds and ionic interactions. Optimised molecular structures indicated the formation of a complex hydrogen bond network, including interresidue and intraresidue bonds. The ionic interactions strongly influence the first hydration shell and, together with hydrogen bonds, play an important role in shaping the 3D tetrasaccharide structure. DFT-derived indirect three⁻bond proton⁻proton coupling constants (³J(H-C-C-H)) showed that the best agreement with experiment was obtained with a weighted average of 67:33 (¹C₄:²S₀) of the IdoA2S forms. Detailed analysis of Fermi-contact contributions to ³J(H-C-C-H) showed that important contributions arise from the oxygen lone pairs of neighbouring oxygen atoms. The analysis also showed that the magnitude of diamagnetic spin⁻orbit contributions are sufficiently large to determine the magnitude of some proton⁻proton coupling constants. The data highlight the need to use appropriate quantum-chemical calculations for a detailed understanding of the solution properties of heparin oligosaccharides.