Excited states of fluorescent proteins were studied using symmetry-adapted cluster-configuration interaction (SAC-CI) method. Protein-environmental effect on the excitation and fluorescence energies was investigated. In green fluorescent protein (GFP), the overall protein-environmental effect on the first excitation energy is not significant. However, glutamine (Glu) 94 and arginine (Arg96) have the red-shift contribution as reported in a previous study (Laino et al., Chem Phys 2004, 298, 17). The excited states of GFP active site (GFP-W22-Ser205-Glu222-Ser65) were also calculated. Such large-scale SAC-CI calculations were performed with an improved code containing a new algorithm for the perturbation selection. The SAC-CI results indicate that a charge-transfer state locates at 4.19 eV, which could be related to the channel of the photochemistry as indicated in a previous experimental study. We also studied the excitation and fluorescence energies of blue fluorescent protein, cyan fluorescent protein, and Y66F. The SAC-CI results are very close to the experimental ones. The protonation state of blue fluorescent protein was determined. Conformation of cyan fluorescent protein indicated by the present calculation agrees to the experimentally observed structure.

units of the environment give smaller corrections, of ca 100 cm -1 , for the ligand-field excitations energies, which is explained by the larger ligands of this cluster. However, the correction for the energy of the lowest LMCT transition is found to be ca 600 cm -1 for the CASSCF and ca 1300 cm -1 for the CASPT2 calculation.3

In this article, we present the multicore (mc) QM/MM method, a QM/MM method that can optimize the structure of chromophore aggregate in protein. A QM region is composed of the sum of the QM subregions that are small enough to apply practical electronic structure calculations. QM/MM energy gradient calculations are performed for each QM subregion. Several benchmark examinations were carried out to figure out availabilities and limitations. In the interregion distances of more than 3.5-4.0 A, the mcQM/MM energy gradient is very close to that obtained by the ordinary QM/MM method in which all the QM subregions were treated together as a single QM region. In van der Waals complex, the error exponentially drops with the distance, while the error decreases slowly in a hydrogen bonding complex. On the other hand, the optimized structures were reproduced with reasonable accuracy in both cases. The computational efficiency is the best advantage in the mcQM/MM approach, especially when the QM region is significantly large and the QM method used is computationally demanding. With this approach, we could optimize the structures of a bacterial photosynthetic reaction center protein in the ground and excited states, which consists of more than 14,000 atoms.