Institute of Biochemical Physics, Russian Academy of Sciences, Moscow, Russia.
Photochemical properties of the green fluorescent protein (GFP) [1,2]which is the final
emitter in the bioluminescence reaction of the jellyfish Aequorea Victoria have motivated a
number of experimental and theoretical studies. Furthermore, other recently discovered
fluorescent proteins (FPs) homological to GFP reveal even more surprising properties of
important practical applications. Due to that, detailed theoretical studies of FPs structure and
properties as well as their chromophores in various environments (protein matix, water
solution, gas phase and etc) form a wide scope of modern research. Gas phase GFP
chromophore has been under investigation in quite a few theoretical studies in the passed
decade, methods ranging from semiempirical to highly accurate ab initio multiconfigurational
theories with perturbation theory corrections. One of the important aspects of such research is
investigation of cis to trans isomerisation of the model chromophore.
As a model chromophore for gas or condenced phase studies and modelling one can use 4'-
hydroxybenzylidene-2,3-dimethylimidazolinone (HBDI) which closely resembles electronic
properties of denaturated GFP chromophore. He et al[3] obtained values for cis to trans
isomerisation barrier along with energy difference of the isomers of HBDI in water solution by
using NMR technique. The height of the barrier was found to be incredibly low compared to
other systems having structural similarity. To our best knowledge no theoretical modelling
results of this reaction consistent with the experiment are reported, although understanding of
such processes is crucial for better comprehension of behavior of several GFP-like proteins
(asFP595, Dronpa).
We apply modern quantum chemistry methods to modelling of HBDI both in gas phase
and solution. Our focus is on modelling cis to trans isomerisation and properties of the model
chromophore consistent with experimental data using modern ab initio methods. The study
reveals problematic behavior of density functional theory (DFT) due to multiconfigurational
character of the wave function at twisted geometry whereas complete active space self
consistent field (CASSCF) yields reasonable results and explains DFT failure. Solvent effects
are found to be important and can reduce activation energy by more than 10 kcal/mol. PC
GAMESS[4] and GAMESS (US) [5]quantum chemistry software were used to perform the
calculations. We also used original hybrid QM/MM flexible fragments technique which is
based on PC GAMESS and Tinker[6] programms. We report results of CASSCF(12/11)/cc-
pVDZ approximation in cluster of 49 TIP3P water molecules (first solvation shell) to be in
good agreement with experimental data[3].
Authors thank SKIF-Siberia and SKIF-MSU for providing computational resources
References
1. M. Zimmer. Chem. Rev. 102, 759 (2002).
2. R.Y. Tsien. Annu. Rev. Biochem. 67, 509 (1998).
3. X. He, A.F. Bell, P.J. Tonge. FEBS Letters 549, 35 (2003).
4. A.V. Nemukhin, B.L. Grigorenko, A.A. Granovsky. Moscow State Univ. Res. Bull. Khimia
45, 75 (2004).
5. M.W. Schmidt, K. K. Baldridge, J. A. Boatz, S. T. Elbert, M. S. Gordon, J. H. Jensen, S.
Koseki, N. Matsunaga, S. Su. Nguyen, T. L. Windus, M. Dupuis, J. A. Mont-gomery. J.
Comput. Chem. 14, 1347 (1993).