Nishioka Hirotaka, Kimura Akihiro, Yamato Takahisa, Kawatsu Tsutomu, Kakitani Toshiaki
Department of Physics, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan.
J Phys Chem B. 2005 Feb 10;109(5):1978-87. doi: 10.1021/jp046282x.
Electron tunneling routes for the electron transfer from the bacteriopheophytin anion to the primary quinone in the bacterial photosynthetic reaction center of Rhodobactor sphaeroides are investigated by a combined method of molecular dynamics simulations for the protein conformation fluctuation and quantum chemical calculations for the electronic states of the donor, acceptor, and protein medium. The analysis of the tunneling route is made by mapping interatomic electron tunneling currents for each protein conformation. We found that there are two dominant routes mainly passing through Trp(M252) (Trp route) or mainly passing through Met(M218) (Met route). Actual electron tunneling pathways alternate between the two routes, depending on the protein conformation which varies with time. When either the Trp route or the Met route dominates, the electron tunneling matrix element /T(DA)/ becomes large. When both the Trp route and the Met route dominate, /T(DA)/ becomes very small due to the destructive interference of the electron tunneling currents between the two routes. We found that a linear relationship exists between the value of /T(DA)/ and the inverse of the degree of destructive interference Q for a wide range of values (ca. 3-10(3) for Q). A similar relationship was also found previously for electron transfer in ruthenium-modified azurins, suggesting that this relationship holds true in general. From these results, we are led to the conclusion that /T(DA)/ cannot exceed a maximum value at Q = 1, even if much variation of /T(DA)/ happens due to the fluctuation of protein conformation. We also conclude that the property of the electron transfer alternates between constructive and destructive interference, due to the fluctuation of protein conformation. It is impossible to keep a system in either constructive or destructive interference because thermal fluctuation of protein conformation takes place.
通过对蛋白质构象波动进行分子动力学模拟以及对供体、受体和蛋白质介质的电子态进行量子化学计算的组合方法,研究了球形红细菌细菌光合反应中心中电子从细菌脱镁叶绿素阴离子转移到初级醌的电子隧穿途径。通过绘制每种蛋白质构象的原子间电子隧穿电流来分析隧穿途径。我们发现有两条主要途径,一条主要通过色氨酸(M252)(色氨酸途径),另一条主要通过甲硫氨酸(M218)(甲硫氨酸途径)。实际的电子隧穿途径在这两条途径之间交替,这取决于随时间变化的蛋白质构象。当色氨酸途径或甲硫氨酸途径占主导时,电子隧穿矩阵元/T(DA)/变大。当色氨酸途径和甲硫氨酸途径都占主导时,由于两条途径之间电子隧穿电流的相消干涉,/T(DA)/变得非常小。我们发现,在很宽的值范围内(Q约为3 - 10³),/T(DA)/的值与相消干涉程度Q的倒数之间存在线性关系。之前在钌修饰的天青蛋白的电子转移中也发现了类似的关系,这表明这种关系通常是成立的。从这些结果中,我们得出结论,即使由于蛋白质构象的波动/T(DA)/发生很大变化,在Q = 1时/T(DA)/也不能超过最大值。我们还得出结论,由于蛋白质构象的波动,电子转移的性质在相长干涉和相消干涉之间交替。由于蛋白质构象的热波动,不可能使系统保持在相长干涉或相消干涉状态。