Semiempirical PM5 molecular orbital study on chlorophylls and bacteriochlorophylls: comparison of semiempirical, ab initio, and density functional results.

The semiempirical PM5 method has been used to calculate fully optimized structures of magnesium-bacteriochlorin, magnesium-chlorin, magnesium-porphin, mesochlorophyll a, chlorophylls a, b, c(1), c(2), c(3), and d, and bacteriochlorophylls a, b, c, d, e, f, g, and h with all homologous structures. Hartree-Fock/6-31G* ab initio and density functional B3LYP/6-31G* methods were used to optimize structures of methyl chlorophyllide a, chlorophyll c(1), and methyl bacteriochlorophyllides a and c for comparison. Spectroscopic transition energies of the chromophores and their 1:1 or 1:2 solvent complexes were calculated with the Zindo/S CIS method. The self-consistent reaction field model was used to estimate solvent shifts. The PM5 calculations predict planar structure of the porphyrin ring and central position of the four coordinated magnesium atoms in all pigments studied, in accord with the experimental, ab initio, and density functional results, a significant improvement as compared to the older semiempirical PM3 approach. Only small differences in PM5 and B3LYP/6-31G* or Hartree-Fock/6-31G* minimum energy geometries of the reference molecules were observed. Calculations show that in 1:1 solvent complexes, where the magnesium atom is five coordinated, the magnesium atom is shifted out of the plane of the porphyrin ring towards the solvent molecule, while the hexa coordinated 1:2 complexes are again planar. The PM5 method gives atomic charges that are comparable with those obtained from the Hartree-Fock/6-31G* and B3LYP/6-31G* calculations. The single point ZINDO/S CIS calculations with PM5 minimum energy structure gave excellent correlations between calculated and experimental transition energies of the chlorophylls and bacteriochlorophylls studied. Such correlations may be used for prediction of transition energies of the chromophores in protein binding sites. Calculations also predict existence of dark electronic states below the main Soret absorption band in all chromophores studied. The results suggest that the semiempirical PM5 method is a fairly reliable and computationally efficient method in predicting molecular parameters of porphyrin-like molecules. Copyright 2003 Wiley Periodicals, Inc. J Comput Chem 1: 123-137, 2004