AM1/d-CB1: A semiempirical model for QM/MM simulations of chemical glycobiology systems

A semiempirical method based on the AM1/d Hamiltonian is introduced to model chemical glycobiological systems. We included in the parameter training set glycans and the chemical environment often found about them in glycoenzymes. Starting with RM1 and AM1/d-PhoT models we optimized H, C, N, O, and P atomic parameters targeting the best performing molecular properties that contribute to enzyme catalyzed glycan reaction mechanisms. The training set comprising glycans, amino acids, phosphates and small organic model systems was used to derive parameters that reproduce experimental data or high-level density functional results for carbohydrate, phosphate and amino acid heats of formation, amino acid proton affinities, amino acid and monosaccharide dipole moments, amino acid ionization potentials, water-phosphate interaction energies, and carbohydrate ring pucker relaxation times. The result is the AM1/d-Chemical Biology 1 or AM1/d-CB1 model that is considerably more accurate than existing NDDO methods modeling carbohydrates and the amino acids often present in the catalytic domains of glycoenzymes as well as the binding sites of lectins. Moreover, AM1/d-CB1 computed proton affinities, dipole moments, ionization potentials and heats of formation for transition state puckered carbohydrate ring conformations, observed along glycoenzyme catalyzed reaction paths, are close to values computed using DFT M06-2X. AM1/d-CB1 provides a platform from which to accurately model reactions important in chemical glycobiology.