Correlation between square of electron tunneling matrix element and donor-acceptor distance in fluctuating protein media.

Correlation between fluctuations of the square of electron tunneling matrix element TDA2 and the donor-acceptor distance RDA in the electron transfer (ET) reaction from bacteriopheophytin anion to the primary quinone of the reaction center in the photosynthetic bacteria Rhodobacter sphaeroides is investigated by a combined study of molecular dynamics simulations of the protein conformation fluctuation and quantum chemical calculations. We adopted two kinds of RDA ; edge-to-edge distance REE and center-to-center distance RCC . The value of TDA2 distributed over more than 5 orders of magnitude and the fluctuation of the value of RDA distributed over more than 1.8 Å for the 106 instantaneous conformations of 1 ns simulation. We made analysis of the time-averaged correlation step by step as follows. We divide the 106 simulation data into 1000/t parts of small data set to obtain the averaged data points of <TDA2> t and <REE > t or <RCC > t . Plotting the 1000/t sets of log10 <TDA2> t as a function of <REE > t or <RCC > t , we made a principal coordinate analysis for these distributions. The slopes <βE > t and <βC > t of the primary axis are very large at small value of t and they are decreased considerably as t becomes large. The ellipticity for the distribution of <TDA2> tvs <REE > t which can be a measure for the degree of correlation became very small when t is large, while it does not hold for the distribution of <TDA2> tvs <RCC > t . These results indicate that only the correlation between <TDA2> t and <REE > t for large t satisfies the well-known linear relation ("Dutton law"), although the slope is larger than the original value 1.4 Å-1. Based on the present result, we examined the analysis of the dynamic disorder by means of the single-molecule spectroscopy by Xie and co-workers with use of the "Dutton law".

The electron transfer (ET) in biological systems takes place by means of the electron tunneling between donor (D) and acceptor (A) in protein media. Usually the distance between the donor and acceptor in such ET is separated by 5–20 Å. In order to make such a long-range ET possible, the superexchange mechanism in which the electronic state of the protein media is virtually incorporated into the electronic states of donor and acceptor to facilitate the electron tunneling works1. The electron tunneling proceeds by passing the bonds with lower energy barriers and the spaces with shorter distance under the restriction that the overall pathway length should be short as much as possible2. As a result, the chosen tunneling pathways and the ET rate are expected to be sensitive to the kind of protein conformation3.

Experimentally it was found that the ET rate k in protein media has much correlation with the donor-acceptor distance R as follows: where β is a parameter with a value 1.4 Å−1 [4]. Here, the distance was obtained by calculating the edge-to-edge distance of π-orbitals of donor and acceptor molecules with the closest distance. This relation was obtained by combining the various data of photosynthetic ET systems and others in stationary state. We call this empirical relation the “Dutton law” hereafter. More sophisticated form of the “Dutton law” was also given with a concept of the packing density of atoms5.

Theoretically, the ET rate k in the elastic tunneling mechanism6 is written as follows1: where < > represents the thermal average and (FC) denote the thermally averaged Franck-Condon factor or called the nuclear factor. Marcus showed that (FC) is expressed by the Gaussian function of the free energy difference ΔG between the initial and final states, which is called the Marcus energy gap law1.

The R-dependence of (FC) enters through the reorganization energy and the energy gap. However, those R-dependences are weak so far as the R is larger than about 5 Å7,8, which applies to many redox proteins. The R-dependence in equation (1) mostly comes from the R-dependence of . Therefore, we write where A is a constant and we replaced R with the average value because the x-ray crystallographic value of R is already averaged. Equation (3) is the relation obtained in the ordinary thermal equilibrium state.

On the other hand, recent theoretical calculations using the MD simulations and quantum chemistry showed that instantaneous values of the electron tunneling matrix element T are greatly varied under the fluctuation of the protein conformation9. The protein conformational fluctuation involves the fluctuation of conformation of amino acids as well as the fluctuation of donor-acceptor distance R. Recently experimental studies using the single-molecule spectroscopy were made for extracting dynamical information of the function and the mechanism of proteins10–13. In such case, protein conformation fluctuation plays a significant role. In some case, the conformation fluctuation with time scale longer than the biochemical reaction is present. This slow fluctuation was called the dynamic disorder14. Under such situation, Xie and co-workers exploited a new experimental technique to derive this dynamic disorder using the fluorescence quenching phenomenon due to the ET reaction15,16. When the distance between the fluorescer and the quencher is small, the ET reaction rate is large and the fluorescence quenching time is reduced and vice versa. By observing the fluctuation of the fluorescence quenching time, the fluctuation of the ET rate was obtained, and the fluctuation of the fluorescer-quencher distance was calculated by using the “Dutton law” of equation (1). Applying this new method of single-molecule spectroscopy combined with the photon-by-photon technique15,16 to flavin reductase, Xie and co-workers obtained a detailed information of the dynamic disorder of protein in the time range of milliseconds to minutes17,18. From this result, they derived power law of the distance autocorrelation function. These results together with their earlier study13 prompted many theoretical studies19–23.

Here, there remains a question whether the similar relation to the “Dutton law” holds between the fluctuating ET rate and the fluctuating distance R in single protein. So far, no one examined this problem in protein environment, although some calculations were made for the fluctuating ET in solutions24,25.

In the present study, we theoretically examine what kind of relationship holds between the fluctuating ET rate and fluctuating R. Based on this analysis, we examine the validity of using the “Dutton law” in the analysis of the dynamic disorder. The constitution of this paper is as follows: In the Method section, the calculation method is described. We applied the method to the ET from the bacteriopheophytin anion (BPhe−) to the primary quinone (QA) in the bacterial photosynthetic reaction center (BPRC) of Rhodobacter sphaeroides. In the Results section the results are presented. The Discussion section is devoted to the discussion. In the Conclusion section conclusion is made.

Method

The procedures for the MD simulations are the same as before6,9,33. Here, we briefly describe it. The initial configuration of the reaction center (RC) was obtained from the Protein Data Bank, entry code 1AIJ 26. The MD simulations are performed for the whole RC including 828 amino acids, all cofactors, 262 water molecules, and 4 detergent molecules of lauryldimethylamineoxide (LDAO) which are attached on the surface of the RC. The phytil chain of the primary quinine is included up to C31. We restrict the position of the hydrophobic phytil chain in the region C21 to C31 since the carbon atoms of this region are located outside the protein surface and no LDAO molecule is attached to them33. The MD program PRESTO27 was used with the AMBER force field28. The SHAKE algorithm29 was used for bond stretchings of hydrogen atoms. We first obtained the energy minimum conformation of protein. Then, the temperature T of the system was gradually elevated to 300 K. Integration time step of 1 fs was employed. After 600 ps run for equilibration, we generated a trajectory for 1 ns. We collected the conformations of 106 at every 1 fs, which were used to calculate T and R.

The procedures for the quantum chemistry (QC) calculations of T are the same as before6,9. Here, we briefly describe it. The chemical structures of donor and acceptor which we adopted in calculating the molecular orbitals of donor and acceptor are shown in Figure 1. In these calculations, long hydrocarbon chains of native bacteriopheophytin a and the native ubiquinone 10 are truncated to hydrogen atoms. It has been shown that the similar result of the time variation of T is obtained in both cases when 60 amino acids surrounding donor and acceptor are considered and when 3 amino acids between donor and acceptor are considered9. Then we consider the three amino acids TrpM252, MetM218, and HisM219 in the present calculations. The coordination of the three amino acids is shown in Figure 1. The electronic structure of donor and acceptor were calculated by the PM3 method30 in the Gaussian package31. To save the computation time, we fix the atomic orbital coefficients in the molecular orbitals which are solved in the isolated state of BPhe and QA. The electronic structures of the pruned protein are solved at the extended Hückel level. We referred to the FORTICON8 program32. We take into account the variable hyperconjugation effect of the methyl group with the π-conjugated part of quinone in the course of rotation of the methyl group33. The tunneling matrix element T is calculated for each protein conformation by the same method as before6,9. The tunneling energy is chosen at −9.5 eV.

Figure 1.

General view of the coordination of donor (bacterio-pheophytin), acceptor (quinone), and the three amino acids (TrpM252, MetM218, and HisM219) which we investigate in this paper. The molecular structure of donor and acceptor is represented by a ball and stick model. The molecular structure of the three amino acids is represented by a stick model. Blue balls, red balls, deep blue balls, and white balls indicate carbon atoms, oxygen atoms, nitrogen atoms, and hydrogen atoms, respectively. The gold stick in the MetM218 represents the sulfur atom.

Results

We first define two kinds of the distance between donor and acceptor. The first is the shortest edge-to-edge distance of π-conjugated region, R, which measures between oxygen atom bonded to C131 of BPhe and C6 of QA. This kind of edge-to-edge distance was used by the Dutton group4. The second is the center-to-center distance R which measures the distance between the center of the polycyclic ring of BPhe and the center of the aromatic ring of QA. This kind of distance was used by Luo et al.23. In Figure 2 we represents these two kinds of distance by the straight line (R) and the dotted line (R). In the same Figure, the average tunneling routes (Trp route and Met route) which were recently obtained are drawn by the red pipes (the pipe width is proportional to the magnitude of the mean-square of the tunneling current corrected by the quantum interference effect)33. The distance R is closely along the Trp route and roughly along the Met route.

Figure 2.

Visualization of the edge-to-edge distance (black straight line) and the center-to-center distance (black broken line) drawn on the calculated map of the mean-square tunneling currents corrected by the quantum-interference (taken from Ref. 33). Protein conformation together with that of donor and acceptor are drawn by the ball and stick model. The orange and green balls indicate the atoms constituting the conjugated part and non-conjugated part of donor and acceptor, respectively. The red pipes represent the electron tunneling pathways whose widths are proportional to the magnitude of . We observed two main tunneling routes: Trp route passing through Trp residue and Met route passing through a part of Met.

In Figure 3 we show a typical example of the time variation of T2. We find that T2 changes very quickly (in tens fs) and the very sharp peaks appear randomly. Such quick fluctuations of T2 also occurs in the other redox proteins34–37.

Figure 3.

Example of rapid variations of T2 as a function of time t.

In Figure 4 we show the calculated correlation diagram of log T2 and R. In this diagram, the 10,000 data points were plotted at every 100 fs from the simulation data of 1 ns. We find that the values of T2 distribute over more than 5 orders of magnitude for the distribution of about 2 Å of R. We see that very little correlation exists between the values of log T2 and R calculated for the instantaneous conformations of fluctuating protein.

Figure 4.

Plot of T2 as a function of R for the 104 points at every 100 fs in the 1 ns simulation data. The averaged point is shown by the yellow star.

In Figure 5 we plotted the correlation diagram of log T2 and R calculated by same way as Figure 4. We find that the values of T2 distribute over more than 5 orders of magnitude for the distribution of about 1.8 Å of R. In this case also we find very little correlation between the values of log T2 and R calculated for the instantaneous conformations of fluctuating protein.

Figure 5.

Plot of T2 as a function of R for the 104 points at every 100 fs in the 1 ns simulation data. The averaged point is shown by the yellow star.

Next, we investigate the dynamical property of the correlation between T2 and R or R. We calculate the mutual correlation function r(τ) as follows where < > denotes the average taken for the time range T, R is R or R, and τ is the delay time between T2 and R.

We calculated r(τ) using the 106 sets of data for T2 and R in the time range T=1 ns. In Figure 6, we plotted the calculated r(τ) as a function of τ by the red curve in case of R= R and the green curve in case of R= R. We find that r(τ) of the red curve has a rather sharp peak while r(τ) of the green curve has broad peak. We also find that these r(τ)’s are almost symmetrical with respect to τ=0. The value of r(0) in the case R= R or R= R is −0.23 or −0.15, respectively. The value of |r(0)| represents how much synchronously the variation of T2 takes place with that of R. The maximum value of |r(0)| is 1. The sign of r(0) is plus or minus depending on the variation of T2 is in phase or out of phase with that of R. The above value of r(0) for R denotes that the variations of T2 and R are out of phase with a certain amount of synchronism to one another. The smaller value of |r(0)| for R denotes that the degree of synchronism is smaller than that for R. We define the time where the value of r(τ) becomes 1/e time of r(0) as the mutual correlation time τ. We read 90 fs for the red curve and 240 fs for the green curve in Figure 6. We consider that the dynamical correlation between the fluctuations of T2 and R continues on time scale of τ. Therefore, we can say that the dynamical correlation between T2 and R is maintained for larger time than that between T2 and R. In Figure 7, we plotted the traces of the data points of log10 T2 and R during the time τ starting from six typical times in the 1 ns simulation. We see that the traces of the data points abruptly change. We tried to fit the data points to the following linear function The straight lines in Figure 7 show those fitted lines. We observe that the slope β changes abruptly form traces to traces. This fact indicates that almost no dynamical correlation exists between ln T2 and R in the very short time range of tens fs. The similar result was obtained for ln T2 and R.

Figure 6.

Plot of the mutual correlation function r(τ) for T2 and R or R.

Figure 7.

Traces of the data points of log10 T2 vs R for 90 fs starting from the time at 90, 180, 270, 360, 450, and 540 ps in the 1 ns simulation data. The straight line is the best fitted line to equation (5).

Next we investigate the correlation between T2 and R by taking time average step by step. We divide the 1 ns simulation data (106 data points) into 1000/t parts of time length t to obtain the averages and or . When t is very large, the number of the averaged data points of , and is very small. In such case, we shift the starting time for dividing the 1 ns simulation data in plural ways until we can obtain more than 20 averaged data points. For example, when t=200 ps, we adopt the starting times for dividing the simulation data at 40i ps (i=0, ..., 3) to give 21 averaged data points as a total. In Figure 8, we showed the calculated diagram between log10 and for three kinds of t; 2, 20, and 200 ps represented by the symbols of green plus, blue closed square and red closed circles, respectively. We also marked the totally averaged (for 1 ns) values of log10 and by the yellow star. We see that the distribution of the averaged data points appears to converge to a straight line as t becomes large. Then, we try to fit the averaged data points to the following function: For this purpose, we made the principal coordinate analysis for the scattered data points. Here, we treat the coordinate log10 and the coordinate equivalently in the analysis. The primary axis corresponds to the linear function of equation (6). Then, we calculated the standard deviation <σ1> along the primary axis, the standard deviation <σ2> along the secondary axis, and the ratio <σ2>/<σ1> (≡ <ε>) which is the ellipticity when the distribution of the averaged data points for each t is approximated as ellipsoid as a whole. In Table 1 we listed the calculated values of the slope <β> and the ellipticity <ε> for 8 kinds of t. We see that <β> is very large at small value of t and it is decreased smoothly as t is increased (<β> = 8.452 for t = 1 ps; <β> = 1.817 for t=333 ps). We have drawn these primary axes for t=2, 20, and 200 ps by the green, blue, and red lines, respectively, in Figure 8. We clearly observe that the slope of the primary axis decreases very much with time. In Table 1, we see that <ε> remains at a constant level of about 0.5 until about 10 ps and then it rapidly decreases down to 0.155 as t is increased from 10 ps to 333 ps. When the ellipticity is small, we can say that the correlation between log10 and is very strong. Then, we conclude that the correlation between log10 and becomes stronger and converges to a straight line more and more as t is increased.

Figure 8.

Plot of log10 as a function of for three kinds of t (2, 20, and 200 ps). The primary axis for the distribution of each time average is drawn. The totally averaged point is represented by the yellow star.

Table 1.

Calculated values of the slope <β> and ellipticity <ε> for some typical average times in the case of edge-to-edge distance

t	<βE>t (Å−1)	<ε>t
1 ps	8.452	0.486
2 ps	7.327	0.484
5 ps	4.206	0.519
10 ps	3.085	0.481
20 ps	2.594	0.399
100 ps	2.166	0.262
200 ps	1.916	0.251
333 ps	1.817	0.155

We also calculated the correlation diagram of log10 and and plotted in Figure 9. We fit the data points to the following linear function: Then, we made the principal coordinate analysis for the scattered data points for each time. In Table 2 we listed the calculated values of the slope <β> and the ellipticity <ε> for 8 kinds of t. We see that <β> is large at small value of t and it is decreased as t is increased (<β> = 8.944 for t=1 ps; <β> = 2.456 for t=333 ps). We also see that <ε> remains almost at the level of 0.5 to 0.6 until t=333 ps. This fact indicates that the correlation between log10 and remains weak for all the value of t until 333 ps. We have drawn the primary axes for t=2, 20, and 200 ps by the green, blue, and red lines in the Figure 9. We clearly see that the slope <β> decreases very much with time.

Figure 9.

Plot of log10 as a function of for three kinds of t (2, 20, and 200 ps). The primary axis for the distribution of each time average is drawn. The totally averaged point is represented by the yellow star.

Table 2.

Calculated values of the slope <β> and ellipticity <ε> for some typical average times in the case of center-to-center distance

t	<βC>t (Å−1)	<ε>t
1 ps	8.944	0.477
2 ps	9.958	0.469
5 ps	7.058	0.520
10 ps	5.290	0.541
20 ps	4.213	0.477
100 ps	3.812	0.532
200 ps	2.684	0.563
333 ps	2.456	0.615

In order to see the time variation of <ε> and <β> or <β> in more detail, we plotted in Figure 10 the diagram for the ellipticity <ε> and the slope <β> or <β> using the data of 90 kinds of t for 1–10 ps, the data points of 18 kinds of t for 10–100 ps and the data points of 8 kinds of t for 100–333 ps. We find that the red line which passes the red circles sharply declines when t is larger than 10 ps and approaches the stationary state value represented by the red square. The evaluation of this point is discussed in the next section. In the case of <β>, <ε> does not decrease until 333 ps.

Figure 10.

Plot of <ε> vs <β> or <β> for the various values of t. The red pluses represent the data points for the time in the range 1–10 ps, the red open circles represent the data points for the time in the range 10–100 ps and the red closed circles represent the data points for the time in the range 100–333 ps for the edge-to-edge distance. The red solid line is extrapolation to the larger time. The red square denoted by S represents the data point of <ε>S and <β>S in the stationary state. The green symbols are the data points for the center-to-center distance. The notation is the same as that of the red symbols.

Using the 106 data points for R and R of 1 ns simulation, we calculate a mean force potential (MFP) as a function of R or R as follows: where P(R) and P(R) are the frequencies of R and R which appear in each bin of 0.01 Å, respectively, and N is the total number of the data points. In Figure 11, we plotted the calculated MFP. We found that the simulation points of F(R) in the upper graph arrange almost in the parabolic form in the region 8.2–10.2 Å. We see that the mean-force potential F(R) is fitted to the quadratic function very well, although some scatterer is seen in the region 10.2–10.7 Å. The green line is the best fitted curve of the quadratic function of R. Similarly, we found that the simulation points of F(R) of the lower graph arrange almost in the parabolic form in the region 13.2–14.9 Å. We see that the MFP is fitted to the quadratic function (green line) very well. The curvature of F(R) is 1.58 time as large as that of F(R). This fact indicates that the change of the center-to-center distance is harder than that of the edge-to-edge distance.

Figure 11.

Mean force potentials as a function of R (top) and as a function of R (bottom). The calculated data points are represented by the red pluses. The fitted parabolic potentials are drawn by green lines.

Discussion

In the first, we evaluate the lower limit of <β> in Figure 10 which corresponds to the value in the stationary state. As we see in Figure 11A, the mean force potential (MFP) as a function of R fits very well with the parabolic function in the region 8.5 Å < R < 10.0 Å. Therefore, we consider that the fluctuation of the distance in this region is well in the stationary state. We make the histogram for the distribution of R using the data of Figure 4. We divide the R coordinate into some parts with a bin of 0.02 Å. We count the number of the data points belonging to each bin and select those bins in which the number of the data points involved is larger than 1000. Next we calculate the average value of T2 in each selected bin. In such treatment, the time correlation among T2(t)’s and R (t)’s is completely neglected. So, we can consider that the average in each bin corresponds to the stationary state. We denote the averaged value of T2 in each bin as S.

In Figure 12, we plot the value of log10 S as a function of R. Then, we make the principal coordinate analysis. The obtained primary axis is expressed as Here, the value 1.736 Å−1 is the slope <β>S in the stationary state. We also evaluate the ellipticity <ε>S as 0.058. This data point is plotted in Figure 10 by the red square named S. In Figure 10 we find that the value of <ε> decreases rapidly and ceases to stop at the stationary state as t is increased from 100 ps. The value of <β> decreases slowly from about 2.1 Å−1 to the lower limit as t is increased from 100 ps.

Figure 12.

Plot of log10 S vs R. S is calculated by averaging T2 which belongs to the bin of R. The bin length is 0.02 Å. The green line represents the primary axis.

Let us consider the mechanism of the property that the values of <β> and <β> are decreased as t is increased and converge to the stationary state value. In the present paper, we have adopted the principal coordinate analysis. In such analysis, the primary axis is chosen along the axis where the dispersion of the data points becomes the largest. When the average time t is small (such as t=2 ps in Figure 8), the dispersion along the log10 axis is much larger than that along the or the axis. Then, the slope <β> of the primary axis is large. When t is increased, the dispersion along the log10 axis is decreased considerably but that along the or the axis is decreased only slightly. Due to such different behaviors, the slope <β> or <β> of the primary coordinate is decreased when t is increased. The decrease of the slope <β> or <β> ceases when t is very large. The mechanism by which the dispersion along the log10 axis is decreased more rapidly than that along the or axis with increase of t is considered as follows: The variation of T2 with time arises from the variation of the amino acids conformation as well as the variation of or . Fast variation of the most amino acids conformation is averaged more rapidly than the slow variation of or . By such mechanism, the dispersion along the log10 axis decreases more rapidly than that along the or axis. When t is very large, the fast variation of the amino acids conformation is almost averaged out and the averaged conformation will change in accordance with the variation of or . Under such situation, the distribution of log10 vs or is virtually determined by the variation of or . By such mechanism, <β> or <β> converges to a certain value when t is very large.

In Figure 10 we have also shown that the ellipticity <ε> for the edge-to-edge case decreases very much when t is very large, while <ε> for the center-to-center case does not decrease from 0.5 when t is very large. This fact will indicate the followings: As we have shown in Figure 2, almost corresponds to the ET route while R is not directly correlated with the ET route. Therefore, the correlation between log10 and is expected to be a limited one. We can generally say that the correlation between log10 and is strong when <ε> is small. Therefore, we conclude that the edge-to-edge distance R is a better measure of the donor-acceptor distance for the ET than the center-to-center distance. However, we should remind of the fact that the above conclusion will hold true only when the edge-to-edge line is close to the ET route.

So far we have used the average time t as a parameter. This average time has a significant meaning as follows: The probability for the ET is accumulated for each infinitesimal time (1 fs in the present simulation). When the probability is attained to a value of 1/e at time t, the ET takes place on the average. Unless the ET rate is very fast, it is reasonable to estimate the value of t by the ET time (inverse of the ET rate in the thermal equilibrium state). The ET time in the present system was experimentally obtained as 200 ps38. If we choose t=200 ps, <β> = 1.916 Å−1. Therefore, the β value in the newly found “Dutton law” which holds between the fluctuating distance and fluctuating ET rate in one protein is considerably larger than that of β (1.4 Å−1) corresponding to the ordinary “Dutton law” which holds between the average distance and average ET rate among many redox proteins.

Now, we refer to the analysis of the single-molecule spectroscopy by Yang et al. for flavin reductase17. The ET time in this system is around 200 ps17. Therefore, if the property of the fluctuation of this system is assumed to be similar to ours, the variation of the due to variation of R in one state remains as represented by the red points in Figure 8. The value of <β> at t=200 ps is 1.916 Å−1. Using the data of Figure 8, we evaluated the standard deviation for the fluctuation of as 0.086 Å and for the fluctuation of as 1.06 time on the average. These values are relatively small although the used data points (21 points) are much limited. Therefore, we can say that some effect of the rapid fluctuations of R and in one state (local energy minimum) may be reflected in the average life time γ−1 in the single-molecule experiment, but its effect will be small as compared with that of the fluctuation of the possible dynamic disorder due to the transition among many local energy minima.

Using the population of the average life time γ−1, Yang et al. calculated the mean force potential (MFP) as a function of the R which fluctuates in the time regions of milliseconds to tenseconds17. We call the fluctuation in the region of milliseconds to tenseconds the slow fluctuation, hereafter. In such calculation they made use of the ordinary “Dutton law”. The calculated MFP is simulated by the parabolic form with curvature 3.2 kcal/mol/Å2, even though this MFP retains certain skewness. On the other hand, Luo et al. calculated the MFP for flavin reductase in the picoseconds region and its curvature was about 8.6 kcal/mol/Å2 when we read from Figure 2A of Ref. 23. This value is not so much different from the curvature 11.1 kcal/mol/Å2 of the MFP obtained by the present study in Figure 11B. (It should be reminded that the distance calculated by Luo et al. is for the center-to-center distance23). From these results we may say that the stochastic property of the distance fluctuation in the picoseconds region is rather similar between the reaction center in our system and the flavin reductase. We call the fluctuations in the region of picoseconds the fast fluctuation, hereafter.

Luo et al. also calculated the MFP by applying nine kinds of external potentials to extend the region of the distance fluctuation23. By means of the umbrella sampling technique, they calculated the MFP over the wide range of the distance from 6.25 Å to 9.4 Å23. It should be noted that this MFP is made of a convolution of the MFP due to the fast fluctuation and the MFP due to the slow fluctuation whose effect was supposed to be incorporated by the umbrella sampling. (Similarly it should be noted that the dispersions of atomic fluctuations in a protein in crystal are made of the dispersions due to the convolution of fast fluctuation and slow fluctuation when one compares the experimentally obtained fluctuation data with the crystallographic data17). In order to obtain the MFP due to only the slow fluctuation, it is necessary to deconvolute the MFP due to the fast fluctuation. The deconvolution in the MFP is approximately performed by subtracting the width of the MFP due to the fast fluctuation from the width of the total MFP. This width of the MFP is tentatively measured at the 1 kcal/mol height from the bottom. The calculated result is listed in Table 3. The calculated widths of the MFP of the present system and that of the single-molecule spectroscopy are also listed in the same table. The width (1.0 Å) of the MFP due to the fast fluctuation by Luo et al. is similar to the width (0.82 Å) of the MFP of the present system measured by the center-to-center distance although the former is a little larger than the latter. (This result is consistent with the above discussion about the curvature of the MFP). This fact indicates that the property of the fast fluctuation of the present system is rather similar to that of flavin reductase. The width (0.5 Å) of the MFP due to the slow fluctuation of Luo et al. measured by the center-to-center distance will be enlarged to 0.64 Å if it is rescaled by the edge-to-edge distance by the same ratio as Figure 11 in the present system. The width (1.7 Å) of the MFP due to the slow fluctuation obtained by the single-molecule spectroscopy by Yang et al. is too large as compared with the above corrected value 0.64 Å. However, if we recalculate the MFP of Yang et al. with use of β=2.0 Å−1 instead of 1.4 Å−1, the width becomes 0.83 Å, which is close to the above corrected value 0.64 Å of the MD simulation due to only the slow fluctuation. These results indicate that the β value of the “Dutton law” which holds for the slow fluctuation should be also around 2.0 Å−1.

Table 3.

Width of the mean force potential at 1 kcal/mol height from the bottom. All the units are in Å

present ET system		flavin reductase
MD simulation		single-moleculea	MD simulationb
edge-to-edge	center-to-center	edge-to-edge	center-to-center
fast	fast	slow	fast and slow	fast	slow
1.05	0.82	1.7 (0.83)c	1.5	1.0	0.5 (0.64)d

Yang et al.17

Luo et al.23

Corrected using β =2.0 Å−1.

Rescaled by the edge-to-edge distance.

Summarizing the above, we found that the width of the MFP due to the fast fluctuation is similar between the present system and flavin reductase. If one assumes that the β value in the slow fluctuation is 2.0 Å−1 or larger, the width of the MFP obtained by the single molecule spectroscopy becomes rather close to the MFP calculated by the MD simulation of umbrella sampling which was supposed to incorporate the slow fluctuation. This fact would indicate that the β value of the “Dutton law” which holds in the fluctuating state is about 2.0 Å−1 for the fast and slow fluctuations.

The origin of the above slow fluctuation will be the slow movement of the large portion of the protein. The amplitude of this movement of protein conformation must be of considerably large scale. Accompanied with this slow and large scale of movement of the protein, a slow and small scale of passive movement of the donor-acceptor distance will be induced. At present it is not certain whether the magnitude of the small scale of movement of the donor-acceptor distance is proportional to that of the large scale of movement or not. It is true that only in the case when their movements are directly correlated to each other, the slow and small scale of the movement of donor-acceptor distance can be used as a measure of the slow and large scale of movement of the protein conformation which deserves the name of dynamic disorder. Its check should be made hereafter.

In the above discussion, we have shown that the linear relationship holds true between the fluctuating distance and the fluctuating ln for the ET in the RC and also for the ET in flavin reductase. The former ET takes place inside the membrane protein and the latter ET takes place near the surface of the globular protein. Then, we deduce that the above “Dutton law” for the fluctuating state will hold true for many kinds of redox proteins.

Conclusion

In the present study we investigated the correlation between fluctuations of the square of electron tunneling matrix element and the donor-acceptor distance for the ET in the photosynthetic reaction center. At first we found that the distribution in the correlation diagram of log10 T2 vs R or R for the instantaneous conformations is very broad. We also found that the dynamic correlation between log10 T2 and R is very poor even within the correlation time. The distribution of the time-averaged diagram of ln vs or was much reduced as the average time t is increased. Applying the principal coordinate analysis, we determined the primary axis from which we obtain the slope <β> or <β>. We defined the ellipticity <ε> as the ratio between the standard deviation along the secondary axis and that along the primary axis for each t. The time averaged ellipticity <ε> of the edge-to-edge case remains large value at small value of t until about 10 ps but it becomes rapidly small when t is increased from 10 ps. This fact indicates that the linear relationship between the fluctuation of ln and the fluctuation of holds true very well. We found that the value of the slope <β> or <β> decreases very much when t is increased and converges to a certain value in the stationary state. The lower bound of <β> expressed by <β>S is 1.74 Å−1, which is larger than the β-value 1.4 Å−1 in the ordinary “Dutton law”. On the other hand, <ε> of the center-to-center case remains large value even for large values of t. Then, the correlation between ln and is worse than that between ln and . Combining the present results with the MD simulation results of Luo et al. and the single-molecule spectroscopic study by Yang et al. for flavin reductase, we deduced that the linear relationship between the fluctuation of ln and the fluctuation of with larger slope than that of the ordinary “Dutton law” will holds true in many kinds of redox proteins.