| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |



* Division of Biophysics, Faculty of Biology/Chemistry, University of Osnabrück, Osnabrück, Germany;
A.N. Belozersky Institute of Physico-Chemical Biology, Moscow State University, Moscow, Russia; and
Institute of Electrochemistry, Russian Academy of Sciences, Moscow, Russia
Correspondence: Address reprint requests to Wolfgang Junge, Abt. Biophysik, FB Biologie/Chemie, Universität Osnabrück, D-49069, Osnabrück, Germany. Tel.: 49-541-969-2872; Fax: 49-541-969-2262; E-mail: junge{at}uos.de.
| ABSTRACT |
|---|
|
|
|---|
threefold) from pH 6 to 10. The maximum conductance was
10 fS at pH 8, equivalent to 6240 H+ s1 at 100-mV driving force, which is an order-of-magnitude greater than of coupled F0F1. There was no voltage-gating of F0 even at low voltage, and proton translocation could be driven by
pH alone, without voltage. The reported voltage gating in F0F1 is thus attributable to the interaction of F0 with F1 but not to F0 proper. We simulated proton conduction by a minimal rotary model including the rotating c-ring and two relay groups mediating proton exchange between the ring and the respective membrane surface. The data fit attributed pK values of
6 and
10 to these relays, and placed them close to the membrane/electrolyte interface. | INTRODUCTION |
|---|
|
|
|---|
) and generates torque. The peripheral F1-portion synthesizes ATP at the expense of this mechanical energy. Both F0 and F1 are rotary motors/generators and steppers. They are mechanically coupled by a central rotary shaft and held together by a peripheral stator (for general reviews, see Senior et al., 2002
F0 is composed of three types of subunits in a stoichiometry of a1b2c1014. All three types of subunits are required for proton conduction (Schneider and Altendorf, 1987
; Altendorf et al., 2000
). Multiple copies of the hairpin-shaped c-subunit form a ring in the membrane to which subunits a and b are attached from the outer side (Lightowlers et al., 1988
; Angevine and Fillingame, 2003
; Fillingame and Dmitriev, 2002
; Jiang et al., 2001
; Muller et al., 2001
; Seelert et al., 2000
; Stock et al., 1999
). ATP hydrolysis in F0F1 is coupled to the rotation of the c-ring (Pänke et al., 2000
; Sambongi et al., 1999
; Tsunoda et al., 2000
, 2001
), and power is elastically and without slip transmitted between the drive in F1 and the ring of F0 (Cherepanov et al., 1999
; Cherepanov and Junge, 2001
; Pänke et al., 2001
).
The currently favored mechanism of torque generation by proton flow (Junge et al., 1997
; Elston et al., 1998
) has been based on the following assumptions: Brownian rotary motion of the c-ring relative to subunit a, two non-co-linear access channels for protons from either side, and two electrostatic constraints, namely the acid Glu residue (Asp in Escherichia coli) in the middle of the hairpin of c being deprotonated when in contact with subunit a and protonated when facing the core of the membrane (Junge et al., 1997
; Elston et al., 1998
). Modifications to this basic scheme have been formulated to account for Na+-translocation in F0 from Propionigenium modestum (Dimroth et al., 1999
) or to account for internal librational motion in the ring (Fillingame et al., 2000
). The basic principle of this rotary mechanism of proton conductionBrownian fluctuations and two electrostatic constraintshas remained (see Junge, 1999
). It is obvious that a comprehensive understanding of rotary proton conduction requires both structural and kinetic information.
Most recent studies on this enzyme have focused on the F1-portion, whereas most studies on the mechanism of proton conduction by F0 date back by more than one decade. The assessment of the conductance of F0 has been hampered by complications:
Thus it is not surprising that the reported values for the conductance at neutral pH diverged by four orders of magnitude. Some studies on the conductance of F0 have yielded rather low values in the range of 0.1 fS, which is equivalent to a rate of 62 s1 at a driving force of 100 mV (Negrin et al., 1980
; Friedl and Schairer, 1981
; Schneider and Altendorf, 1982
; Sone et al., 1981
; Cao et al., 2001
). This value would be insufficient to cope with the turnover rate of the coupled F0F1 enzyme. For chloroplast CF0 we have found a 100-fold-higher figure, 9 fS (or 5600 s1 at 100 mV) under the assumption that all exposed F0 molecules contributed to the relaxation of the transmembrane voltage (Schönknecht et al., 1986
). This conductance is compatible with the maximum proton turnover rate of the coupled enzyme,
1000 s1. Later we obtained circumstantial evidence for a majority of inactivated F0, which raised the conductance to a value of 1 pS (Lill et al., 1986
; Althoff et al., 1989
). Such high rates by far exceeded the diffusive proton supply to F0, as then noted, and it has remained questionable whether the underlying assumption was correct. This work has also provided evidence for an extreme proton specificity of the conductance (>107 over other cations), small pH-dependence between 5.5 and 8, a H/D-kinetic isotope effect of 1.7, and an activation energy of 42 kJ/mol in H2O and 47 in D2O. The latter data were not affected by the above ambiguity over the absolute magnitude of the conductance (Althoff et al., 1989
).
A conductance of similarly high magnitude was reported for CF0CF1 in a lipid bilayer formed by the dip-stick technique, which is related to patch-clamp (Wagner et al., 1989
). Gated single channel currents were observed. They peaked at 0.55 pA at 180 mV, implying a conductance of 0.4 pS, and channels were gated open with a sharp activation above 100 mV. The authors attributed this voltage-gated conductance to the proton (Wagner et al., 1989
). Whether this attribution was correct has remained an open question, in particular because proteoliposomes, which contained the purified c-subunit alone, have revealed unspecific cation channels (Schönknecht et al., 1989
; for the proton conduction at pH 2 of the purified c-subunit see Schindler and Nelson, 1982
).
The lack of information on the number of active F0 "motors" per membrane area was a major obstacle for obtaining a reliable estimate of its proton conductance. Our recent studies on isolated chromatophore vesicles from the photosynthetic bacterium Rhodobacter capsulatus (Feniouk et al., 2001
, 2002
) have paved a way to overcome the ambiguity over the proportion of active and inactive F0. The key was to prepare vesicles of such small size as to contain less than one copy of F0 on the average. Then, any rapid relaxation of the transmembrane voltage could be attributed to the subset of vesicles containing mainly a single F0 molecule.
In this work we prepared vesicles with 28-nm mean diameter as checked by transmission electron microscopy (TEM). Excitation of a suspension of such vesicles with a flash of light generated a voltage step across the membrane. The magnitude of the initial voltage jump after a saturating flash of light (70 mV) was calibrated by comparison with electrochromic absorption transients that were induced by submitting vesicles from the same batch to a salt-jump. The relaxation of the voltage was monitored by electrochromic absorption transients of intrinsic carotenoids, and the proton flow by pH-indicating dyes. The relaxation was biphasic. The fast phase (25%) was mediated by F0 as evident from its sensitivity to the F0-inhibitor oligomycin. It was mainly attributable to proton conduction in the subset of vesicles containing a single copy of F0.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
|
Chromatophore size determination by TEM
Chromatophores were suspended in a medium with 100 mM KCl, 20 mM glycylglycine, 5 mM MgCl2, pH 8.0 to yield a bacteriochlorophyll concentration of 5 µM. A drop of the suspension was applied on the carbon-coated grid, stained with 3% uranyl acetate, washed with water, dried and examined in a Zeiss 902 transmission electron microscope (Carl Zeiss, Oberkochen, Germany). F1-depleted chromatophores appeared as disks with diameters of 39.7 ± 1.0 nm (SD 7.9 nm; see Results; also see Discussion).
Flash spectrophotometry
The kinetic flash-spectrophotometer was described elsewhere (Feniouk et al., 2001
). The bandwidth of the measuring beam, 8 nm, was set by interference optical filters (Schott, Mainz, Germany). Changes in transmitted light intensity (
I) were monitored by a photomultiplier (Thorn EMI, 9801B, UK) shielded from the actinic flash with two blue filters (BG 39, 4 mm, Schott). The electrical bandwidth was 3 kHz and the digital time-per-address of the averager was 200 µs. The optical path was 1 cm, both for the exciting and for the measuring beam. The final concentration of bacteriochlorophyll in the cuvette was 812 µM.
Transient signals were generated by flashing light at a repetition rate of 0.08 Hz. Eight signals at 522 nm were averaged in recordings of the electrochromism and 32 in measurements with pH indicators. During the dark time (12.5 s) between flashes the electrochromic signal decayed to <5% of its initial value after the flash of light, even if the conductance of F0 was blocked by oligomycin (not documented). The actinic flashes were provided by a xenon flash lamp (full-width at half-maximum = 10 µs) with a red optical filter (RG 780, Schott); the total energy density on the cuvette (not attenuated) was 12 mJ cm2. The peaks of the xenon emission spectrum (e.g., 825 and 875 nm) overlap with regions of low absorption of the chromatophores (their peaks are at 800 and 850 nm). We calculated an effective optical density for the excitation by convoluting the xenon emission and chromatophore absorption. With the highest concentration (12 µM bacteriochlorophyll) in the cuvette it was 0.6 OD, one-half of the value (1.2 OD), which was calculated for the peak absorption at 850 nm. In control experiments with lower bacteriochlorophyll concentration (6 µM), the effective optical density was 0.3, which provided a rather homogenous excitation profile over the thickness of the optical cell. Neutral density filters were used to attenuate the actinic light flash when indicated (see Fig. 5 legend).
Electrochromic carotenoid bandshifts at 522 nm were used as a molecular voltmeter to monitor 
(Junge and Witt, 1968
; Clark and Jackson, 1981
; Symons et al., 1977
; Jackson and Crofts, 1971
).
Transients of the pH inside chromatophores (
pHin) were monitored by absorption changes of the amphiphilic pH indicator Neutral red at 545 nm as in Ausländer and Junge (1975)
and Mulkidjanian and Junge (1994)
. In these experiments 0.3% BSA was used as an impermeable buffer to quench pH changes in the suspending medium (Ausländer and Junge, 1975
). pH transients in the outer medium (
pHout) were monitored by absorption changes of the hydrophilic pH indicator Cresol red at 575 nm (Saphon et al., 1975
) and calibrated in pH units. For each pH indicator, control traces without the indicator were recorded and subtracted.
Calibration of the carotenoid bandshift in millivolts of transmembrane electrical potential difference
Chromatophores were suspended in a medium containing 20 mM glycylglycine (pH 8.0), 5 mM MgCl2, 2 mM NaCN, 3 µM myxothiazol, 3 µM oligomycin, and KCl. The KCl concentration was varied, and NaCl was added to keep the total concentration of KCl+NaCl constant at 250 mM. Valinomycin (1 mM stock solution in EtOH, final concentration in the cuvette 1.6 µM) was added to induce a diffusion potential of K+; the resulting absorption transients at 522 nm were recorded in the same instrument used for flash spectrophotometry.
| RESULTS |
|---|
|
|
|---|
|
If the transmembrane voltage jump was generated by the photosynthetic reaction centers (RCs) in response to a single saturating flash of light, the electrochromic absorption transient was composed of two contributions, one being the response to the delocalized and transmembrane field and the other to the local electric field originating from the oxidized primary electron donor. In the calibration procedure, only the first component had to be considered. We eliminated the delocalized transmembrane component by adding valinomycin, 1 µM, and estimated its magnitude (see Fig. 2). The residual valinomycin-insensitive absorption changes (
5% of the total flash-induced jump in our sample, see Fig. 2) relaxed in the range of 100 ms, following the reduction of the primary electron donor in the RC. Although the relative extent of this component was only 5% in the samples studied in this work, it can reach 15% if no redox mediator is present in the medium.
|
We estimated the number of RCs in these vesicles. The estimate was based on the voltage of 70 mV, on the usual figure for the specific capacitance of biomembranes,
= 1 µF cm2 (Hille, 1984)
, and on the mean vesicle radius of 14 nm, as determined by TEM. The capacitor equation reads: n x eo = A x
x
U, wherein n denotes the number of RCs per vesicle, eo, the elementary charge, A the vesicle area, and
U the transmembrane voltage. We obtained a value of n = 11 as the mean number of RCs in our preparation of F1-depleted vesicles. This number implied a surface density of 7.3 x 109 mol m2,
40% larger than the previously determined 5.3 x 109 mol m2 (Packham et al., 1978
).
Electrochromic absorption transients with and without active cytochrome bc1 complex
Fig. 2 A shows flash-induced absorption transients at 522 nm in F1-depleted chromatophores. The rapidly decaying lower traces were obtained with valinomycin. Their decay reveals the magnitude of the electrochromic response to the transmembrane and delocalized electric field. The upper traces were obtained with and without added oligomycin (Linnett and Beechey, 1979
) as an inhibitor of proton conduction through F0. The slowly rising lower traces give the difference between the upper two, corresponding to the cumulative charge flow through F0. The rise of absorption after the actinic flash was biphasic (see Fig. 2 A, trace with oligomycin added). The sharp (here not time-resolved) increase is attributable to the charge separation in the RC. Its magnitude is fairly independent of pH, redox potential, and temperature (Jackson, 1988
). We used the extent of the fast rise to normalize transients that were obtained in different samples. In Fig. 2 A the fast onset was followed by a slowly rising phase. It was attributable to the electrogenic reaction in the cytochrome bc1 complex (bc1; Crofts and Wraight, 1983
; Jackson, 1988
). This slow phase was abolished upon the addition of the bc1-inhibitor myxothiazol (see Fig. 2 B).
The decay of the electrochromic transient was rapid if F0 was conducting and slow if it was blocked. The latter proved the low leak conductance of the membrane. A clear-cut interpretation of the difference trace in Fig. 2 A was hampered by the activity of the cytochrome bc1 complex for two reasons: 1), bc1 translocates electrons and protons, and thereby generates a transmembrane pH jump in addition to a voltage jump; and 2), the turnover of bc1 is hindered by a large
. In other words, the turnover of bc1 is not independent of the presence of the F0-inhibitor oligomycin. These effects were eliminated by inactivating the bc1-complex with myxothiazol, which abolished the slow rise of the electrochromic transient as documented in Fig. 2 B. The RCs generated a voltage step of the same magnitude as when bc1 was active (compare Fig. 2, A and B; also see Mulkidjanian and Junge, 1994
). Only in this case does the difference trace (± oligomycin) reflect the relaxation of the transmembrane voltage step in those vesicles that have a single or, in a minor fraction only, several F0 complexes. Comparison of the relaxation times in Fig. 2, A and B, revealed the systematic error when assessing the voltage relaxation with active bc1 (the apparent relaxation time was 711 ms as compared with 24 ms in the presence of myxothiazol). In the first case, F0 first transferred the protons driven by the RC-generated 
and then "waited" for the protons that were ejected by the cytochrome bc1-complex. The rate of proton transfer through F0 in the presence of myxothiazol was almost independent of the presence of a penetrating pH-buffer glycylglycine (see Table 1).
|
22% of total contained a single copy of F0, and
3% two copies.
The discharge time of the vesicles' electric capacitance by F0 was then apparent as the relaxation time of the difference trace (± oligomycin) as in Fig. 2 B. The alternative interpretation of the biphasic decay of the voltage, the discharge of the membrane's electric capacitance by a voltage-gated channel, was ruled out by the observation of the same biphasic behavior after a lower starting voltage (nonsaturating energy of the exciting flash; see Fig. 3 in Feniouk et al., 2002
).
|
1 above pH 8 (Fig. 4 A). The extent of the charge flow via F0 was pH-independent within the experimental error (when care was taken to avoid the aging of the preparation). The Arrhenius activation energy of the electric relaxation was 57 ± 4 kJ mol1 at pH 8, and increased to 72 ± 3 kJ mol1 upon acidification (Fig. 4 B).
|
The Ohmic behavior of F0
The following data on the electric conductance of F0 were obtained in the presence of myxothiazol, to avoid the complications caused by the electrogenic and proton pumping activity of bc1. We changed the initial voltage jump by varying the flash energy. With a mean number of RCs per vesicle of 11 (see above) the initial extent of the voltage could be varied over one order of magnitude.
Fig. 5 shows two difference transients (± oligomycin) reflecting the charge transfer through F0. These transients were obtained at 100% and 33% saturation (judging from the initial extent of the absorption jump at 522 nm), respectively. The extent of the initial step of absorption was normalized. It was obvious that the two transients perfectly superimposed. The transient at nonsaturating flash energy was only discernible by greater noise from the one obtained at saturating flash energy (circles). The relaxation was almost perfectly monoexponential (fit curve not shown in Fig. 5) and independent of the initial voltage jump. The same relaxation time was found when the flash intensity was attenuated to yield 10% saturation (not documented).
The original traces in Fig. 5 were obtained at a bacteriochlorophyll concentration of 12 µM where the effective optical density of the sample in the excitation region was 0.6 (see Materials and Methods). It caused a drop of the flash energy from the entry to the exit of the cuvette (1-cm thickness) from 100% to 25%. This did not matter for the traces with 100% signal saturation, because the flash energy was far-oversaturating, but it did matter for the traces at 33%. In this respect, it was noteworthy that we obtained the same relaxation times both in samples with halved bacteriochlorophyll concentration (6 µM, optical density 0.3, flash energy drop across the cuvette from 100% to 50%) and in samples with 10% saturation. It justified the neglect of the saturation profile in the following interpretation of our data.
The traces at the bottom of Fig. 5 show the difference between the normalized 100% and 33% transients; the deviation between them was very small indeed. Such a voltage-independent relaxation characterizes the discharge of a capacitor (area-specific capacitance =
/AsV1 m2) by an Ohmic resistor (area-specific conductance =
/AV1 m2). It yields a voltage-independent relaxation time,
=
/
. F0 behaved as an Ohmic resistor, and it showed no voltage-gating.
The rapid relaxation of the voltage (relaxation time =
2.9 ms at pH 8; see Fig. 3) was only observed in the small subset of vesicles containing F0 (26% of total). The electric capacitance of vesicles
was calculated by their size using the usual figure for biological membranes, 1 µF cm2 (Hille, 1984)
. The conductance of F0, g, was calculated according to g = A x
with
=
/
. Herein A denotes the vesicle area,
the electric relaxation time, and
and
are the area-specific conductance and electric capacitance of the membrane, respectively.
The shortest relaxation time observed, 2.3 ms, gave rise to a calculated conductance of 10.5 fS. It implied the transport of 6500 protons/s at 100-mV steady driving force or of 13,000/s at 200 mV.
It was noteworthy that the proton-transporting properties of the Rb. capsulatus F0 were robust against the changes of the concentration of K+, Na+, Cl, of various permeating and nonpermeating buffers, and redox mediators (not documented).
The proton specificity of F0
We complemented measurements of the electric relaxation with measurements of the pH transients in both phasesthe suspending medium (hydrophilic indicator Cresol red, wavelength 575 nm) and the interior of chromatophores (amphiphilic indicator Neutral red, 545 nm, with BSA as a nonpermeating buffer). The results are documented in Fig. 6. The data were obtained with active bc1, to demonstrate the proton uptake from the bulk medium by the RC and by bc1 (Fig. 5, A and B) and the proton release into the interior of chromatophores by bc1 (Fig. 5, C and D). The alkalization of the bulk and the acidification of the interior were fully apparent when F0 was blocked by oligomycin (traces 2), and both were diminished when F0 was conducting (traces 1). The respective difference traces gave the F0-related pH-relaxation, however, with somewhat distorted kinetics because of the slow proton-pumping activity of bc1 (see above discussion in the context of Fig. 2, A and B). Ignoring the complex kinetics of these difference transients, they proved that the major ion transported by F0 was the proton at an ambient pH 7.9 against a background of 100 mM KCl. Similar results were observed with NaCl and LiCl (data not shown). In other words, the proton selectivity of F0 was at least 107.
|
pH as the only driving force for proton transport: further evidence for the absence of voltage-gating
pH (Schlodder and Witt, 1981It was obvious that F0 was not, by itself, voltage-gated. The observed electrical gating in the intact enzymes from E. coli and P. modestum was likely attributable to the interaction of F1 with F0.
A kinetic model for proton conduction by F0
We investigated whether the magnitude, the Ohmic behavior, and the weak pH-dependence of proton conduction by F0 could be accommodated within the framework of the presently assumed rotary mechanism of proton conduction (Junge et al., 1997
; Elston et al., 1998
). The proton transfer through the bare chloroplast F0 has been explained elsewhere by a model where two proton-conducting half-channels were connected by a (horizontally) rotating carrier (see Fig. 7 A for an illustration). The fit of experimental data yielded a low pK of 56 for the proton-accepting group of the input half-channel and an alkaline pK of 8 for the output group (Cherepanov et al., 1999
). The weak pH-dependence between these two pK values follows from this model. The proton-transfer rate is determined by the probability of proton binding to group A and proton release from group B. When the ambient pH inside and outside are the same, as in our experiments, the product probability is constant in this interval. This feature explained the observed weak pH-dependence of proton transfer through the bare F0. In energized membranes, on the other hand, the pK values of the relay groups tend to match the ambient pH at the p- and n-sides of the coupling membrane. This allows higher turnover numbers of F0 even against backpressure (Cherepanov et al., 1999
).
|
The main elements of the model were two noncoaxial proton-conducting half-channels both considered to rapidly equilibrate with their adjacent bulk phase. In the simplest case, each half-channel contained only one proton-binding relay group (A at the acidic side of the membrane and B at the basic side). In the absence of a transmembrane voltage, the rate of proton binding to the groups A and B was assumed to depend on the pH in the respective bulk phase in accordance with the equations
![]() |
The first term in each equation denotes the proton flux owing to the diffusion of H3O+ ions and the second term describes the hydrolysis of neutral water (see Eigen, 1963
). The factor
was introduced to account for the accelerated proton supply from surface buffering groups. The addition of penetrating pH bufferglycylglycine was tested at concentration up to 0.02 Mcaused no further acceleration of proton transfer through F0; see Table 1. This implies a rapid proton delivery along the inner chromatophore surface to F0 from many ionizable groups acting as proton antennae. (See e.g., Zhang and Unwin, 2002
, for a quantitative examination of surface proton conduction.) The estimate of
was obtained by considering the strictly Ohmic behavior of F0 with accuracy >3% (See Fig. 4). Because the overall rate of proton transfer through F0 after a saturating flash was
3.5 x 103 s1 (see above), the linear relationship implied that at neutral pH the rate constant
(which was independent of
) was >105 s1, which corresponds to the
-factor >10 (see also the discussion below). We used the latter value of
in the modeling.
The rate constants of the reverse reactions (see the kinetic scheme in Fig. 7 D) were calculated by the thermodynamic equilibria of
![]() |
The rate of proton exchange between F0 and the bulk solutions was high and not limiting to the overall reaction rate. The forward and reverse rate constants of the proton transfer in the middle of membrane (m+ and m) were interconnected by the thermodynamic constraint of
![]() |
Thus, there were only three independent kinetic parameters in the model: the pK values of buffer groups in each of the two half-channels, and the rate constant of the elementary rotary step (the rotation of the cn-ring by the angle 2
/n), which brings one proton into and another one off the ring in contact with its respective access channel.
The potential energy of a transferred proton within F0 is schematically plotted in Fig. 7 C. The transmembrane electric field affects both the protonation state of the two relay groups (A and B) in the two half-channels, and the rate of proton transfer over the center of the membrane. The magnitude of the first effect depends on the (dielectrically weighted) relative depth of groups A and B in the membrane (denoted
in Fig. 7 C, L is the full and 1 the relative thickness of the membrane). The second effect depends on the effective width of the barrier, denoted ß in Fig. 7 C. The absolute height of the barrier is not important, if it is much greater than the thermal energy kBT.
A step of the transmembrane electric potential, magnitude
(in Volts, positive at the A-side) shifts the pK values as
![]() |
![]() |
) depends on the shape of potential barrier, which for a rectangular profile is
![]() |
![]() |
![]() |
![]() |
These three differential equations have to be completed by the normalizing condition
![]() |
Pij denotes the probability of protonation (index 1) and deprotonation (index 0) of the relay groups in contact with the left (left index) and the right (right index) access channel.
The fourth differential equation describes the membrane recharging
![]() |
The solution of this system of four differential equations was obtained by numerical integration.
Herein F denotes the Faraday constant, N0 the surface density of F0 in the membrane (mol m2), and Cm the electric capacitance of membrane (F m2). Because the average diameter of chromatophores was
28 nm and they contained a mean of
0.3 copies of F0 per vesicle, the mean surface density of F0 in vesicles, which contained enzyme N0, was
6.6 x 1010 mol m2. The membrane's specific capacitance was taken as 102 F m2. On modeling we took a value of 70 mV for the average 
jump in response to a single saturating flash with the blocked bc1 complex. The fit parameters pKA, pKB, and m+ were determined from the pH-dependence of the F0 conductance by nonlinear minimization. The final set of model parameters is summarized in Table 2. The resulting theoretical pH-dependence of the proton transfer rate in H2O is plotted in Fig. 3 by the solid line. The obtained values pKA = 6.1, pKB = 10.0, and m+ = 1 x 107 s1 corresponded to the overall F0 turnover rate
.
Fig. 5 A shows the normalized transients of the transmembrane voltage at pH 8.0 after a saturating (100%) and an attenuated (33%) flash, the difference between them is plotted in Fig. 5 B by the thin noisy line. This difference, which is very small, characterizes the deviation from the Ohmic behavior. The insert to Fig. 5 B shows the distribution of the deviation between two traces, which did not reveal a systematic bias from zero. The calculated mean root-square dispersion of the difference was 3%.
The high signal/noise ratio of the traces allowed us to impose constraints on the possible values of the parameters a, b, and
. Generally, the greater was
and the smaller ß, the faster was proton transfer at a given value of the transmembrane voltage, so the relaxation dynamics was faster after saturating and slower after nonsaturating flashes (positive nonlinearity). The effect of the parameter
was opposite: because the diffusion of protons outside F0 was voltage-independent, the proton diffusion at small values of
slowed down the overall proton transfer (diffusional control of the reaction) and caused thereby a negative nonlinearity in the relaxation kinetics.
The dashed curve in Fig. 5 B shows the difference between the transmembrane voltage decay after a saturating and an attenuated flash calculated with the parameter
= 1. This curve revealed a substantial deviation from the Ohmic behavior, which exceeded the experimental error. The minimal value of
consistent with the experimental data was
10, and this value was used in the following modeling.
The dotted curve in Fig. 5 B was calculated for a set of parameters
= 0.1, ß = 0.8, and
= 10. It illustrates the positive nonlinearity in the relaxation owing to the accelerating effect by the transmembrane voltage. The deviation again exceeded the experimental error, but had the opposite sign.
We calculated the membrane discharge dynamics after saturating and attenuated flashes using the kinetic model described above and determined the range of possible values of
and ß consistent with an experimental error of <3%. Generally the relaxation rate was 10-fold more sensitive to the variation of
(the buried position of groups A and B inside the membrane) than to the variation of ß (the form of the potential barrier in the middle of the membrane). In the two-dimensional phase space (
and ß correspond to the abscissa and ordinate axes of the insert to Fig. 5 A, respectively) we determined the area where the maximal deviation between the discharge curves after saturating and attenuated flashes was being <3% consistent with the experimental error. The allowed space found by this method is shown in the insert to Fig. 5 A by the uncolored upper-left area, and the forbidden area inconsistent with the experimental data is hatched. An example of the calculated difference between two recharging traces with the maximal deviation of 3% is shown in Fig. 5 B by the thick dashed line (the parameters chosen were
= 0.04, ß = 0.64, and
= 10).
In practice, these findings implied that both proton-conducting groups A and B were placed at the membrane/water interface. The dielectrically weighted distance between them had to be >0.86 (the maximal allowed value of
is 0.07). This maximal value of
corresponded, however, to an infinitely sharp rectangular barrier in the middle of the membrane (ß = 1), physically unfeasible. A realistic value of ß is 0.6; it corresponds to an upper limit of
of
0.04. It is worth noting, however, that the positive nonlinearity of the electric field inside the membrane could be partially compensated by the negative nonlinearity of the voltage-independent diffusion outside the membrane, so we used the value
= 0.07 as an estimate of the buried position of groups A and B inside the membrane. If the dielectric permittivity was constant over the full width of the membrane, say of 4 nm, the topological distance between A and B would be >3.4 nm. The dielectric permittivity of the water-filled cavities in membrane proteins may, however, strongly differ from that in the membrane interior (see, e.g., a value of 30 in the ubiquinone binding pocket of bacterial reaction center versus the value of 4 in the middle of membrane; Cherepanov et al., 2000
). In this case, the distance between groups can be substantially smaller, on the order of 1.8 nm (45% of the geometrical membrane thickness). It was obvious though, that if either group was embedded more deeply in the membrane dielectric, one had to expect a marked deviation from an Ohmic behavior.
Our analysis reveals that the proton transfer to the conserved carboxyl residue of subunit c, which is topologically located in the middle of the membrane, is mediated by at least two residues located close to the membrane surface. It remains open how this finding relates to the recent Cys-mapping experiments of Angevine and Fillingame (2003)
, which led these authors to postulate two water-accessible half-channels in the F0F1-ATP synthase of E. coli. The relay groups A and B may be found among those residues lining these half-channels.
| DISCUSSION |
|---|
|
|
|---|
In our related previous studies the average number of F0 per vesicle was also small but larger than in this work (Feniouk et al., 2002
). Those particular preparations were, however, prepared with a sonifier of smaller output power (Branson W-15, Soest, The Netherlands). The average diameter of the chromatophores used in the previous work determined by TEM was 62.8 nm (data not shown), 1.6-fold larger than that of the F1-depleted ones used here (see Fig. 1 A). The larger chromatophores contained approximately one copy of F0F1 on the average, in good correspondence with the F0 mean surface density calculated here.
We found that the electric relaxation through F0 was almost mono-exponential and that the relaxation time was practically independent of the magnitude of the initial voltage step up to 70 mV. F0 behaved as an Ohmic conductor with a linear current-voltage relationship. There was no voltage-gating of F0. The conductance varied little (less than threefold) as function of the pH over the wide range from 6 to 10. The specificity for protons over other cationic species was very high (>107).
There was no kinetic H/D-isotope effect at pH above 8, and an isotope effect of
2 was observed at pH 6 (Fig. 4 A). The appearance of the isotope effect could be attributed to protonation of a group with the apparent pK of 78. A possible candidate for this titratable group might be Glu61 of the c-subunit. The pK of its counterpart in the E. coli enzyme (cAsp61) has the apparent pK of 8 (Valiyaveetil et al., 2002
). An alternative explanation could be secondary isotope effects: the general conformation of the enzyme might be slightly altered by pH, and the H/D-isotope effect could be inherent to the acidic conformation only. In this case a smooth pH-dependence of the H/D isotope effect like that in Fig. 4 A could be expected. The pH-dependence of both the isotope effect and the Arrhenius activation energy (Fig. 4 B) points out that the rate-limiting step of proton transport through F0 is different at pH 6 and 8. Further experiments are necessary to clarify this point.
The above properties of the bacterial F0the high specificity, the weak pH-dependence, and the kinetic isotope effectwere in a good agreement with those of the chloroplast counterpart, where data were obtained over a narrower pH interval, of 5.5 < pH < 8. (Althoff et al., 1989
). The very high proton specificity of F0 was attributable to proton binding by acid/base group(s) as selectivity filter. It excluded models with a contiguous water-wire without such groups. The same property favored a Grotthuss-type mechanism of proton transport, and it excluded the self-diffusion of the proton in the form of H3O+ or OH.
These properties of F0 were simulated based on the current rotary mechanism with the multicopy c-ring carrying one essential acid residue (Glu or Asp) in the middle of the membrane and with one access channel at each side. We assumed one proton relay group in each channel. The fit to the data attributed pK values of
6 and
10, respectively, to these relay groups because of the low dependence on the pH. The Ohmic behavior called for a location of the relays very close to their respective membrane/bulk interface. F0 conducted protons irrespective of the nature of the driving force, whether of electric (transmembrane voltage) or entropic origin (pH difference): Both this property and the Ohmic behavior under an electric driving force proved the absence of voltage gating of F0. The
-gating of the coupled enzyme in chloroplasts (Junge, 1970
, 1987
; Junge et al., 1970
; Gräber et al., 1977
; Witt et al., 1977
; Schlodder and Witt, 1980
, 1981
; Schlodder et al., 1982
) and the voltage-gating in E. coli and P. modestum (Kaim and Dimroth, 1998a
, 1999
; Dimroth et al., 2000
) is attributable to the interaction of F1 with F0, but not to F0 proper.
The straightforward method used to obtain the unitary conductance of one single copy of F0 in this work resolved the discrepancy between previous reports from our lab. In studies with spinach thylakoids we arrived at a similar magnitude for the conductance (9 fS) under the then-unproven assumption that any exposed F0 was conducting (Schönknecht et al., 1986
). Circumstantial evidence for only a minority of conducting F0, which had led us to claim a much higher conductance of up to 1 pS (Lill et al., 1986
; Althoff et al., 1989
), can now be rejected. It also appears conceivable that the similarly high conductance of 0.4 pS as attributed to whole F0F1 in one study with bilayers formed by the dip-stick technique (Wagner et al., 1989
) might have involved unspecific cation channels as observed with the purified subunit c of F0 (Schönknecht et al., 1989
).
Extrapolating the data in this work on F0 to a driving force of 200 mV (nearly sufficient for the maximum rate of ATP synthesis by the coupled enzyme), one expects a transport rate by F0 of 13,000 s1. ATP synthesis by the E. coli enzyme can reach a turnover number of 270 s1 (Etzold et al., 1997
), with coupled ATP hydrolysis in Rhodospirillum rubrum at 320 s1 (Slooten and Nuyten, 1984
), and ATP synthesis in Rb. capsulatus at 250 s1 (calculated from 900 µmol/mg BChl·h, Casadio et al., 1978
, if there was one ATPase per 1000 BChl; Packham et al., 1978
; Feniouk et al., 2002
, and with a proton/ATP stoichiometry of 3.3). Taking the proton/ATP stoichiometry into account, which is probably 3.3 for yeast FoF1 (Stock et al., 1999
) and 4 (van Walraven et al., 1996
; Turina et al., 2003
) or 4.7 (Seelert et al., 2000
) for spinacia, these figures imply the turnover of only
1000 H+/s by the coupled enzyme, F0F1. This is one order-of-magnitude lower than that of the exposed F0.
Comparison with the proton conductance of gramicidin
Gramicidin is perhaps the best studied cation channel and proton conductor. It is gated open when its two half-spanning helices join tail-to-tail to fully span the membrane. Cations are conducted along and with the internal water wire. This is well established by studies on electro-osmosis and streaming potentials. The proton, on the other hand, bypasses the single-file transport (see Finkelstein and Andersen, 1981
; Levitt et al., 1978
; for molecular dynamics simulation see (Skerra and Brickmann, 1987
; Schumaker et al., 2000
). The proton conductance of the covalently linked gramicidin A channels is proportional to the H+ concentration in the pH range from 0 to 4 (Cukierman, 2000
). Extrapolation to pH 8 yields a H+ conductance of 1 fS, 10-fold smaller than the H+ conductance of F0. In simulating the data on F0 we had to assume that the supply of protins to F0 was not rate-limiting, i.e., it was by order-of-magnitude faster than expected by bulk diffusion from a half-space to the channel mouth. It implied supply from several buffering groups at the membrane surface acting as proton donors. Recently, the lateral diffusion coefficient along the lipid/water interface has been determined to be 5.8 x 109 m2 s1 (Serowy et al., 2003
), only twofold lower than in the bulk (9.3 x 109 m2 s1), but greater than for any other ion in the bulk phase and thus compatible with a Grotthuss-mechanism involving surface water (see Gutman and Nachliel, 1997
). Proton supply in the alkaline range is not critical, because water then acts as the proton donor (see Kasianowicz et al., 1987
; Gopta et al., 1999
). These features have been considered in the above simulation.
| CONCLUSIONS |
|---|
|
|
|---|
In F0 the selectivity filter for protons over other cations is highly specific, >107, both in the purple bacterium used in this work and in thylakoids from green plants, as previously studied. The selectivity is still present but less pronounced in organisms that operate on Na+ instead of H+.
We interpreted the properties of F0 in terms of the current rotary model for proton conduction. This model is based on two proton-conducting half-channels linking the rotating ring of 1014 copies of subunit c with the respective bulk phases. The observed Ohmic conduction implies that the relay groups for protons that are involved in rate limitation sample only a small fraction of the transmembrane voltage. We