Three-Dimensional Potential Energy Surface of Selected Carbohydrates' CH/[pi] Dispersion Interactions Calculated by High-Level Quantum Mechanical Methods

• Published in: Chemistry : a European journal 2011-05, Vol.17 (20), p.5680
• Main Authors: Kozmon, Stanislav, Matuska, Radek, Spiwok, Vojtch, Koca, Jaroslav
• Format: Article
• Language: English
• Subjects: Carbohydrates; Chemistry; Hydrocarbons
• ISSN: 0947-6539; 1521-3765
• DOI: 10.1002/chem.201002876

In this study we present the first systematic computational three-dimensional scan of carbohydrate hydrophobic patches for the ability to interact through CH/π dispersion interactions. The carbohydrates [beta]-d-glucopyranose, [beta]-d-mannopyranose and [alpha]-l-fucopyranose were studied in a complex with a benzene molecule, which served as a model of the CH/π interaction in carbohydrate/protein complexes. The 3D relaxed scans were performed at the SCC-DFTB-D level with 3757 grid points for both carbohydrate hydrophobic sides. The interaction energy of all grid points was recalculated at the DFT-D BP/def2-TZVPP level. The results obtained clearly show highly delimited and separated areas around each CHgroup, with an interaction energy up to -5.40kcalmol-1. The results also show that with increasing Hπ distance these delimited areas merge and form one larger region, which covers all hydrogen atoms on that specific carbohydrate side. Simultaneously, the interaction becomes weaker with an energy of -2.5kcalmol-1. All local energy minima were optimized at the DFT-D BP/def2-TZVPP level and the interaction energies of these complexes were refined by use of the high-level ab initio computation at the CCSD(T)/CBS level. Results obtained from the optimization suggest that the CH group hydrogen atoms are not equivalent and the interaction energy at the CCSD(T)/CBS level range from -3.54 to -5.40kcalmol-1. These results also reveal that the optimal Hπ distance for the CH/π dispersion interaction is approximately (2.310±0.030)Å, and the angle defined as carbon-hydrogen-benzene geometrical centre is (180±30)°. These results reveal that whereas the dispersion interactions with the lowest interaction energies are quite strictly located in space, the slightly higher interaction energy regions adopt a much larger space.