P. Palenčár1, F. Vácha1,3, and M. Kutý1,2
1Institute of Physical Biology,
2Institute of Systems Biology and Ecology,
Zámek 136, 37333
3Institute of Plant Molecular Biology,
Branišovská 31, 370 05 České
e-mail: palencar@ufb.jcu.cz
First, we constructed truncated models of PS II considering protein subunits and pigment molecules from and around the RC, keeping also functionally and structurally important spatial subunits. These structures were then used as starting structural models for MD simulation runs. Force field (FF) [3] parameters (charge distribution and force constants) and topology for PS II RC pigment molecules [4] were developed for the Yamber2 FF [5, 6] by performing quantum chemical calculations and by modifying and extending Ceccarelli [7] studies concerning bacterial type a pigments. Geometry optimisation of PS II RC pigment molecules and development of new FF parameters was performed at RHF/6-31G* level of quantum chemical theory using Gaussian 98 [8]. Calculation of RESP atomic charges following AMBER FF developing scheme [9] was done using Gaussian 98, R.E.D, antechamber, and resp subroutines. Complete FF parameters and topology for the PS II RC pigment molecules were developed and introduced to YASARA [5]. Several MD simulation runs, taking into account different pigment oxidation states and various solvent properties, were performed in order to check the quality of new FF parameters for the PS II RC pigments. Subsequently, by further optimising overall geometry of new structural models using MD simulation runs with small time step of 1.0 fs, we obtained final models of PS II RC complex. Detailed analysis of new FF parameters showed realistic dynamic behaviour of pigment molecules and typical electrostatic interactions with the surrounding protein environment.
Recently, changes in excitonic interactions of PS II RC pigments upon light-induced oxidation of primary donor (P680) or reduction of primary acceptor pheophytin a (Phe a), were analysed using absorption and circular dichroism (CD) spectra [10, 11]. In contrast to the oxidation of primary donor, the light-induced change in the CD spectrum upon primary acceptor reduction was temperature-dependent. This suggests a hypothesis that at a room temperature the reduced Phe a induces conformational changes of the RC protein environment, which affects the excitonic interaction of the RC chlorophylls (Chls). Having optimised structural models of PS II RC we were able to elucidate and describe some of the details of these processes.
1.
Ferreira, K.N.,
Iverson, T.M., Maghlaoui, K., Barber, J., Iwata, S.: Architecture of the
photosynthetic oxygen-evolving center. – Science 303: 1831-1838, 2004.
2.
Biesiadka, J., Loll,
B., Kern, J., Irrgang, K.-D., Zouni, A.:
3.
MacKerell, A.D., Jr.:
Empirical force fields for biological macromolecules: overview and issues. – J.
Comput. Chem. 25: 1584-1604, 2004.
4.
Palenčár, P., Vácha,
F., Kutý, M.: Force field development on pigments of photosystem 2 reaction
center. – Photosynthetica 43: 417-420, 2005.
5.
Krieger, E., Darden,
T., Nabuurs, S.B., Finkelstein, A., Vriend, G.: Making optimal use of empirical
energy functions: force-field parameterization in crystal space. – Proteins 57:
678-683, 2004.
6.
Cornell, W.D., Cieplak,
P., Bayly, C.I., Gould, I.R., Merz, K.M.,
7.
Ceccarelli, M.,
Procacci, P., Marchi, M.: An ab initio force field for the cofactors of
bacterial photosynthesis. – J. Comput. Chem. 24: 129-142, 2003.
8.
Frisch, M.J., G.W.
Trucks, H.B. Schlegel, G.E. Scuseria, M.A. Robb, J.R. Cheeseman, V.G.
Zakrzewski, J.A. Montgomery, Jr., R.E. Stratmann, J.C. Burant, S. Dapprich,
J.M. Millam, A.D. Daniels, K.N. Kudin, M.C. Strain, O. Farkas, J. Tomasi, V.
Barone, M. Cossi, R. Cammi, B. Mennucci, C. Pomelli, C. Adamo, S. Clifford, J.
Ochterski, G.A. Petersson, P.Y. Ayala, Q. Cui, K. Morokuma, N. Rega, P.
Salvador, J.J. Dannenberg, D.K. Malick, A.D. Rabuck, K. Raghavachari, J.B.
Foresman, J. Cioslowski, J.V. Ortiz, A.G. Baboul, B.B. Stefanov, G. Liu, A.
Liashenko, P. Piskorz, I. Komaromi, R. Gomperts, R.L. Martin, D.J. Fox, T.
Keith, M.A. Al-Laham, C.Y. Peng, A. Nanayakkara, M. Challacombe, P.M.W. Gill,
B. Johnson, W. Chen, M.W. Wong, J.L. Andres, C. Gonzalez, M. Head-Gordon, E.S.
Replogle, and J.A. Pople: Gaussian 98, Revision A.11.3, Gaussian, Inc.,
Pittsburgh PA, 2002.
9.
Bayly, C.I., Cieplak,
P., Cornell, W.D., Kollman, P.A.: A well-behaved electrostatic potential based
method using charge restraints for determining atom-centered charges: The RESP
model. – J. Phys. Chem. 97: 10269-10280, 1993.
10. Vácha, F., Durchan, M., Šiffel, P.:
Excitonic interactions in the reaction centre of photosystem II studied by
using circular dichroism. – Biochim. biophys. Acta. 1554: 147-152, 2002.
11.
Vácha,
F., Pšenčík, J., Kutý, M., Durchan, M., Šiffel, P.: Evidence for localization
of accumulated chlorophyll cation on the D1-accessory chlorophyll in reaction
centre of photosystem II. – Photosynth. Res. 84: 297-302, 2005.
Acknowledgement
This
research was supported by the Ministry of Education, Youth and Sports of the