pH-dependent 11° F1FO ATP synthase sub-steps reveal insight into the FO torque generating mechanism
Abstract
Most cellular ATP is made by rotary F1FO ATP synthases using proton translocation-generated clockwise torque on the FO c-ring rotor, while F1-ATP hydrolysis can force counterclockwise rotation and proton pumping. The FO torque-generating mechanism remains elusive even though the FO interface of stator subunit-a, which contains the transmembrane proton half-channels, and the c-ring is known from recent F1FO structures. Here, single-molecule F1FO rotation studies determined that the pKa values of the half-channels differ, show that mutations of residues in these channels change the pKa values of both half-channels, and reveal the ability of FO to undergo single c-subunit rotational stepping. These experiments provide evidence to support the hypothesis that proton translocation through FO operates via a Grotthuss mechanism involving a column of single water molecules in each half-channel linked by proton translocation-dependent c-ring rotation. We also observed pH-dependent 11° ATP synthase-direction sub-steps of the Escherichia coli c10-ring of F1FO against the torque of F1-ATPase-dependent rotation that result from H+ transfer events from FO subunit-a groups with a low pKa to one c-subunit in the c-ring, and from an adjacent c-subunit to stator groups with a high pKa. These results support a mechanism in which alternating proton translocation-dependent 11° and 25° synthase-direction rotational sub-steps of the c10-ring occur to sustain F1FO ATP synthesis.
Editor's evaluation
This manuscript describes single molecule spectroscopic work to probe the mechanisms by which protonation/deprotonation steps produce torque between the a-subunit and the C10 ring, which is subsequently conveyed to F1 to couple to ATP synthesis. This is an important bioenergetics question and the approach yields some tantalizing clues as to which protonation steps are involved. In principle, this knowledge could lead to direct experimental or computational tests to resolve the overall mechanism. The previous issues in the manuscript were well addressed by the authors in the new revision. The reviewers found a few small additional issues, which the authors can address in the final revisions.
https://doi.org/10.7554/eLife.70016.sa0Introduction
The F1FO ATP synthase (Figure 1) that is found in all animals, plants, and eubacteria is comprised of two molecular motors that are attached by their rotors and by their stators (Kühlbrandt, 2019; Spetzler et al., 2012). The FO motor, which is embedded in bioenergetic membranes, uses a non-equilibrium transmembrane chemiosmotic proton gradient also known as a proton-motive force (pmf) to power clockwise (CW) rotation of its ring of c-subunits relative to the subunit-a and -b stator proteins as viewed from the Escherichia coli periplasm. These subunits contribute to the peripheral stalk bound to one side of the F1 (αβ)3-subunit ring where each αβ-heterodimer comprises a catalytic site that synthesizes ATP from ADP and Pi. Subunit-γ, which docks to the c-ring along with subunit-ε, extends into the core of the (αβ)3-ring (Figure 1B). The rotary position of subunit-γ within the (αβ)3-ring causes the conformations of the three catalytic sites to differ such that one site contains ADP and Pi, a second site contains ATP, and the third site is empty. Pmf-powered CW rotation of subunit-γ forces conformational changes to all catalytic sites in the (αβ)3-ring, which releases ATP from one catalytic site with each 120° rotational step (Kühlbrandt, 2019; Spetzler et al., 2012). In this manner, F1FO converts the energy from the pmf (ΔμH+) into energy in the form of a non-equilibrium chemical gradient (ΔμATP) where the ATP/ADP•Pi concentration ratio is far in excess of that found at equilibrium.
When ΔGATP is significantly higher than n*pmf, the F1-ATPase motor can overpower the FO motor and catalyze net ATP hydrolysis (Steigmiller et al., 2008; Fischer et al., 2000). This results in ATPase-dependent power strokes that rotate continuously CCW for 120° at saturating ATP concentrations (Yasuda et al., 2001; Spetzler et al., 2006), and pump protons across the membrane to the periplasm in E. coli (Spetzler et al., 2012). Power strokes are separated by catalytic dwells that last a few ms, during which ATP is hydrolyzed (Yasuda et al., 2001; Spetzler et al., 2006; Spetzler et al., 2009). At rate-limiting ATP concentrations, an ATP-binding dwell can interrupt the E. coli F1 power stroke ~34° after the catalytic dwell (Yasuda et al., 2001; Martin et al., 2014). At this same rotary position, high ADP concentrations can compete with ATP to bind to the empty catalytic site resulting in ADP inhibition (Martin et al., 2014). All ATP synthases can catalyze ATP hydrolysis to some extent, although many have evolved regulatory mechanisms to minimize this energy wasteful process (Kühlbrandt, 2019). Under some circumstances, E. coli employs F1FO as an ATPase-driven H+ pump to maintain a pmf as an energy source for other metabolic processes (Spetzler et al., 2012).
The means by which H+ translocation from the periplasm generates CW rotational torque on the c-ring is poorly understood. Molecular motors are believed to operate by either power stroke or by a Brownian ratchet mechanism, which has been postulated for FO (Hwang and Karplus, 2019; Oster et al., 2000). Although evidence clearly supports a power stroke mechanism for F1-ATPase-dependent rotation (Martin et al., 2014; Martin et al., 2018; Pu and Karplus, 2008), there is little direct evidence in support of either mechanism for FO-driven rotation in the ATP synthase direction. Protons enter and exit FO via half-channels in stator subunit-a. Two c-subunits in the ring contact subunit-a at a time where the leading cD61 during CW rotation (synthase direction) accepts a proton from the input channel, while the lagging cD61 donates its proton to the output channel.
The half-channels are separated by the highly conserved subunit-a arginine (aR210 in E. coli), which has long been thought to be responsible for displacing the H+ from cD61 into the output channel during ATP synthesis (Lightowlers et al., 1987; Cain and Simoni, 1989; Cain and Simoni, 1988; Vik et al., 1988). As the result of mutations of E. coli subunit-a residues aN214, aE219, aH245, aQ252, and aE196 to other groups including leucine that decreased ATP synthase activity, ATPase-dependent H+ pumping, and altered the ATP hydrolysis activity, these residues were proposed to translocate protons directly along the half-channels (Lightowlers et al., 1988; Howitt et al., 1990; Eya et al., 1991; Hartzog and Cain, 1994; Hatch et al., 1995; Hahn et al., 2018). Cryo-EM F1FO structures that reveal details of subunit-a confirmed that these residues are positioned along possible half-channels that are separated by aR210 (Hahn et al., 2018; Zhou et al., 2015; Pinke et al., 2020; Sobti et al., 2020; Martin et al., 2015).
Alternatively, H+ translocation through FO has also been postulated to occur via a Grotthuss mechanism (Cukierman, 2006) where a column of single water molecules that are hydrogen-bonded to specific protein groups behave in a coherent manner to transfer protonic charge over long distances via rapid exchange of H+ between H3O and H2O (Cukierman, 2006; Wraight, 2006). A recent F1FO structure from bovine mitochondria was of sufficient resolution to observe density near the input channel residues consistent with Grotthuss-type water molecules in this half-channel (Spikes et al., 2020). A consequence of this coherent behavior of a Grotthuss water column is that the rate of H+ transfer can be much faster than the rate via free diffusion (Cukierman, 2006). The possibility that FO operates via a Grotthuss mechanism was first suggested from the observation of an astounding H+ translocation rate of 6240 H+ s–1 from a driving force of 100 mV across Rhodobacter capsulatus vesicles containing FO that lacked F1 (Feniouk et al., 2004). The R. capsulatus FO rates of H+ transfer exceed the rate of delivery of protons by free diffusion from the bulk aqueous solution at a concentration of 10–8 M (pH 8) such that the ability to supply protons to a Grotthuss water column should be rate-limiting (Wraight, 2006). To achieve this rate of H+ translocation, the existence of a H+ antenna at the entrance to the FO input channel has been postulated (Wraight, 2006), which in R. capsulatus, was calculated to consist of a hemispherical Coulomb cage with a H+ capture radius of ~40 Å surrounding the entrance to the input channel.
Single-molecule studies of F1FO molecules embedded in lipid bilayer nanodiscs (Figure 1C) revealed that the 120° CCW ATPase power strokes can be interrupted by transient dwells (TDs) at ~36° intervals with a duration of ~150 μs (Martin et al., 2015; Ishmukhametov et al., 2010; Yanagisawa and Frasch, 2017). In more than 70% of TDs, the FO motor not only halted F1-ATPase CCW rotation, but the c-ring rotated CW in synthase-direction steps (Martin et al., 2015; Yanagisawa and Frasch, 2017). Complete assembly of F1FO nanodisc complexes from the membrane scaffold protein (MSP), lipids, and detergent-solubilized F1FO was verified by 2D electrophoresis where the first nondenaturing gel dimension contained a single band, which after the second denaturing dimension, contained bands corresponding to MSP and all F1FO subunits (Ishmukhametov et al., 2010). The ATPase activity of nanodisc preparations is DCCD-sensitive, remains unchanged after 8 hr at room temperature, and is 1.5-fold higher than the initial activity of detergent-solubilized F1FO, which loses all activity in <4 hr. The observation that a subunit-a insertion mutation, which disrupts the interface between subunit-a and the c-ring, retained ATPase activity but lost the ability to form TDs indicates that these dwells result from an interaction between subunit-a and successive c-subunits in the c-ring (Ishmukhametov et al., 2010).
The occurrence of TDs was found to increase at pH 8 when viscous drag on the nanorod is sufficient to slow the angular velocity of the F1 ATPase-driven power stroke (Martin et al., 2015; Ishmukhametov et al., 2010). These results showed that there is a kinetic component that affects the probability that the interaction between subunit-a and the c-ring will occur relative to the F1-ATPase power stroke. A kinetic dependence for F1FO-catalyzed ATP synthesis versus hydrolysis has been theorized based on energetic calculations (Gao et al., 2005). Occurrence of TDs, including those with synthase-direction steps, is also known to increase inversely with pH between pH 5 and pH 7 (Yanagisawa and Frasch, 2017). This suggests that synthase-direction steps depend on H+ transfer from the protonated groups with a low pKa from the subunit-a input channel to the c-ring, and from the c-ring to unprotonated groups with a high pKa in subunit-a output channel.
Maximal ATP synthase rates in membrane vesicles where E. coli F1FO is oriented with F1 on the outer surface are typically achieved with inner and outer pH values of 5.0 and 8.5, respectively (Fischer et al., 2000; Fischer et al., 1994). The pmf can be derived from non-equilibrium differences in pH (ΔpH) or membrane potential (Δψ). Although the energy from ΔpH and Δψ to drive ATP synthesis are interconvertible, the latter is the dominant energy source for E. coli in vivo.
Each of the three successive 120° F1-ATPase power strokes required for a full revolution of the F1FO rotor is unique because the rotary positions of subunit-γ differ relative to the peripheral stalk that includes subunit-a (Figure 1A–C). These power strokes also differ because the 36° stepping of the c10-ring and the 120° power strokes can only be aligned during one of the three catalytic dwells. As a result, the three power strokes require the translocation of 4 H+, 3 H+, and 3 H+ that result in net c-ring rotations of 144°, 108°, and 108°, respectively. Consequently, the c-ring and subunit-γ become misaligned by +14° and −14° during two of the catalytic dwells. The elasticities of the peripheral stalk, subunit-δ, and to some extent, subunit-γ accommodate these rotary differences (Sobti et al., 2020; Murphy et al., 2019).
The positive and negative torsion on the c10-ring from the elastic energy needed to accommodate the +14° and −14° misalignments during rotation affects the ability to form TDs and their associated synthase-direction steps (Yanagisawa and Frasch, 2017), which along with their pH dependence and occurrence every 36°, indicate that they correspond to single c-subunit stepping relative to subunit-a. To determine the percentage of power strokes in which TDs were observed, data sets were collected from one of the three 120° power strokes from each single-molecule of nanodisc-embedded E. coli F1FO examined. For each data set that comprised ~300 power strokes, the percent of power strokes that contained TDs was determined. For each molecule examined there was an equal chance that the c-ring and catalytic dwell was aligned, or subject to the positive or negative torsion from misalignment. The distribution of data sets collected from many molecules versus the percent of power strokes containing TDs fit to the sum of three Gaussians that corresponded to low, medium, and high probabilities of TD formation. The high and low TD percentages were consistent with the torsion from misalignment that provides additional energy to promote TD formation, or to inhibit it when torsion is in the opposite direction (Yanagisawa and Frasch, 2017). Similar effects of the c10-ring and catalytic dwell misalignments have also been observed in other single-molecule F1FO measurements (Sielaff et al., 2019).
We have now examined mutations of residues in the putative subunit-a half-channels that have been implicated to participate in H+ transfer events in both the input and output channels with the goal of distinguishing whether FO rotation occurs via a Brownian ratchet versus a power stroke mechanism, and whether the half-channel residues transfer protons individually or act together to support water column that transfers protons via a Grotthuss mechanism.
The ability to form ATP synthase-direction steps as a function of pH during ATP hydrolysis-driven CCW power strokes was characterized in single-molecule studies of F1FO molecules embedded in lipid bilayer nanodiscs. Formation of synthase-direction steps during a TD was maximal at the pH value when unprotonated and protonated forms of the input and output channels were optimal such that synthase-direction steps required H+ transfer both from the input channel to the leading c-subunit and from the lagging c-subunit to the output channel. Mutation of a residue from either the input or output channel altered both the low and high pKa values of TD formation indicating that input and output channels communicate. This is consistent with a Grotthuss mechanism where a water column in each of the two channels is connected by rotation-dependent H+ transfer to and from the leading and lagging c-subunits in the c-ring, which is also supported by features from a variety of F1FO structures.
The extent of rotation during ATP synthase-direction steps was unexpectedly found to rotate 11° in the WT and all mutants. Cryo-EM structures of sub-states with subunit-a:c-ring differences of 11° that position the lagging c-subunit cD61 adjacent to output residues aS199 and aE196 (Sobti et al., 2020) are consistent with the synthase-direction steps observed here. When combined with structural information, the results do not support the hypothesis that the role of aR210 is to displace the H+ from cD61, but are consistent with a Grotthuss H+ translocation mechanism involving both half-channels for sustained ATP synthase-direction c-ring rotation that results from successive alternating 11° and 25° synthase-direction sub-steps for each c-subunit in the c10-ring. Direct evidence of a mixed model was observed in which some synthase-direction steps show characteristics of a power stroke while others exhibit oscillations consistent with a Brownian ratchet mechanism.
Results
Contributions of subunit-a residues putatively involved in the ATP synthase H+ half-channels were assessed by the effects on TD formation caused by mutations that converted charged or polar groups in subunit-a to hydrophobic leucine. Changes in rotational position were measured by a 35 × 75 nm gold nanorod (AuNR) bound to the biotinylated c-ring of individual E. coli FOF1 molecules embedded in lipid bilayer nanodiscs (Ishmukhametov et al., 2010), hereafter F1FO (Figure 1C). Changes in rotational position during F1-ATPase power strokes in the presence of saturating 1 mM MgATP were monitored by the intensity of polarized red light scattered from the AuNR (Spetzler et al., 2006; Hornung et al., 2011). Prior to data collection, the polarizer was adjusted so that the scattered red light intensity was at a minimum during one of the three F1 catalytic dwells (Figures 1D and 2A). The subsequent power stroke caused an increase in light intensity to a maximum when the AuNR had rotated 90° (Ragunathan et al., 2017). Rotational data sets of each F1FO molecule examined were collected for 5 s, which included ~300 of these power strokes (Yanagisawa and Frasch, 2017). Ten data sets were collected for each molecule. The number of F1FO molecules examined at each pH for WT and mutants is indicated in Figure 2—figure supplement 3. Using WT at pH 5.0 as an example where data from 103 F1FO molecules were collected, this was equivalent to 1030 data sets, and ~309,000 power strokes examined. For each molecule examined, rotational position versus time was calculated from scattered light intensity versus time using an arcsine1/2 function from which the number of TDs observed during the first 90° of rotation were determined (Ragunathan et al., 2017).
Example power strokes from WT and mutant F1FO molecules at pH 5.0 where TDs were present (black dots) and absent (blue dots) are shown in Figure 2A and Figure 2—figure supplement 1, respectively. When present, TDs either stopped F1-ATPase CCW rotation momentarily (green dots) or exhibited CW rotation in the ATP synthase direction, hereafter synthase-direction steps (red dots). None of the mutations examined eliminated the ability of F1FO to form TDs. Power strokes typically contained two to three TDs, when present. These were separated by an average of ~36°, consistent with an interaction between subunit-a and successive c-subunits in the c10-ring of E. coli F1FO.
A power stroke mechanism has been defined as the generation of a large free energy gradient over a distance comparable to the step size of the molecular motion so that transition to the forward position occurs nearly irreversibly (Hwang and Karplus, 2019). By contrast, in a Brownian ratchet mechanism the motor is thought to visit previous and forward positions through thermal motion, where stabilization in the forward position results by conformational changes triggered by the fuel processing event. While some synthase-direction steps shown of Figure 2A and Figure 2—figure supplement 1 rotated CW in a concerted, and apparently irreversible manner characteristic of a power stroke, others indicated by brackets were observed to oscillate back and forth during the TD. These oscillations most commonly occurred late in the F1 power stroke (~70–80°) and were more pronounced in all mutations examined except aN214L (Figure 2—figure supplement 1). Such oscillations are direct evidence of a Brownian ratchet mechanism and are likely the result of a close balance between the energy that powers the F1-ATPases power stroke with the energy that powers synthase-direction rotation, which suggests that these mutations cause a decrease in the energy to power synthase-direction rotation.
Subunit-a mutations alter pKas of TD formation
We postulated that mutation of subunit-a residues involved in H+ translocation related to c-ring rotation would alter the pKa of the half-channel in which the is located if each residue contributes independently to the H+ translocation process. Consequently, we precisely determined the pKa values of groups that contribute to TD formation (Figure 2B) using equations applied to the pH dependence of enzyme inhibition kinetics (Cook and Cleland, 2007).
TDs occur when subunit-a binds to the c-ring to stop F1 ATPase-driven rotation for a period of time. Thus, a TD represents an extent that FO inhibited the F1ATPase motor, which occur as often as 3.6 times per F1 power stroke. Kinetically, the ATPase power stroke duration without TDs is ~200 μs, while the average duration of each TD is ~150 μs (Martin et al., 2015; Ishmukhametov et al., 2010; Yanagisawa and Frasch, 2017). In data sets where TDs occur in 100% of the power strokes, for example, aN214L at pH 6.0, all of the ~300 power strokes in that data set will contain TDs (i.e. they look like the power strokes in Figure 2A where TDs are denoted by green or red dots) such that the efficiency of TD formation is 100%. Such a data set represents a 64% inhibition of the F1ATPase power stroke kinetics. When all ~300 power strokes in a data set look like those in Figure 2A represented by blue dots, the efficiency of TD formation in that data set is 0%.
Precise pKa determination depended upon the fits of the efficiency at each pH examined to the pH dependence over the range of pH values in Figure 2B. The number of molecules examined for WT and each mutant is indicated in Figure 2—figure supplement 3. For example, a total of 553 F1FO molecules (~92 million power strokes) were examined to establish the pH dependence of WT. A maximum average of 47.5% of WT power strokes from all three efficiency groups occurred at pH 5.0, which decreased with increasing pH until it plateaued at a minimum of ~22% at pH values > 7.5 (Figure 2B). The pH dependences for WT and mutants were fit to Equation 1 where T is the total average TD occurrence, Tmin is the minimum TD occurrence, and K1 and K2 are the inhibition constants that define the increase and maximum TD occurrence versus pH as the result of either a residue that is protonated with pKa1, or unprotonated with pKa2, respectively. It is noteworthy that K1 is similar to a dissociation constant because a smaller K1 increases the ability of subunit-a to bind to, and stop, c-ring rotation with decreasing pH (Figure 2—figure supplement 2). Conversely, a smaller K2 value decreases TD formation with decreasing pH because it is the unprotonated form of that residue that binds and inhibits.
The fit of the data to Equation 1 defines the slope of the curve as well as the high and low plateau values. Because these are log-log plots, the pKa values (Figure 2B, dotted lines) are determined by the intercept of the slope with the high and low plateau values (solid lines). None of the mutations changed Tmin significantly. Using parameters derived by the fits of the data to Equation 1 for WT and mutants (Table 1), the WT group(s) that must be protonated to induce a TD had pKa1 and K1 of 5.6, and 6.4, respectively, while the group(s) that must be unprotonated to induce a TD had pKa2 and K2 values of 7.5 and 6.75, respectively.
The aN214L mutation, which had the greatest effect on the pH dependence of TD formation, increased the maximum percent of TDs formed at low pH to 61% (1.3-fold) and shifted the pH dependence in the alkaline direction from WT. These changes were due to increases in K1 and K2 to 6.4 and 6.75, respectively, that increased pKa1 and pKa2 by 0.9 and 0.7 pH units. The differential increases in K1 and K2 by 0.6 and 0.75 units led to the aN214L-dependent increase in maximum TD formation at low pH because an equal shift of these values in the same direction causes the curve to shift to higher pH values without affecting the maximum occurrence of TDs formed (Figure 2—figure supplement 2). Similar but smaller effects were observed with aQ252L and aE219L (Figure 2B) where K1 increased by 0.6 and 0.5 units, respectively, resulting in a pKa1 increase of almost 1 pH unit from that of WT. However, aQ252L and aE219L decreased K2 by 0.35 and 0.40 units from WT such that the increase in pKa2 was proportionally smaller than that observed for aN214L. Consequently, while both mutants shifted the pH dependence in the alkaline direction from that of WT, only aQ252L showed an increase in the maximum TD occurrence (52%).
Mutations aH245L and aE196L caused the smallest changes in the pH dependence of TD formation. The former increased K1 and K2 by 0.1 and 0.12 units, which had no effect on pKa1, and increased pKa2 by 0.2 units. The latter was the only mutation to decrease the values of both K1 and K2, which decreased pKa2 by 0.1 pH units from that of WT. It is noteworthy that aE196 is a component of the H+ output channel.
In summary, mutation of residues examined in either the input or the output channel changed both the high and low pKa values of TD formation. These results support the conclusion that they participate in H+ translocation as first proposed from ensemble studies (Lightowlers et al., 1987; Cain and Simoni, 1989; Cain and Simoni, 1988; Vik et al., 1988; Lightowlers et al., 1988; Howitt et al., 1990; Eya et al., 1991; Hartzog and Cain, 1994; Hatch et al., 1995), but also show that the two channels communicate via the c-ring, which would occur if the residues in both channels support a Grotthuss-type water column connected by c-ring rotation.
Subunit-a mutations affect TD formation efficiency
The percent of TDs observed per data set was fit to three Gaussian distributions with low (blue), medium (orange), and high (green) efficiencies as shown at pH 6.0 (Figure 2C), and at all pH values examined (Figure 2—figure supplement 3). These efficiency differences were proposed to result from elastic energy resulting from the 14° rotational mismatch between two of the three catalytic dwells and the c10-ring that supplements or subtracts from the binding energy required for subunit-a to stop F1-ATPase-driven rotation momentarily, resulting in a TD. If TDs result from H+ translocation-dependent interactions between subunit-a and the c-ring, mutations that impact H+ translocation should alter the TD formation efficiency.
Subunit-a mutations affected the percent of TDs formed per data set during power strokes in each of these efficiencies, which correlate to the three rotary positions of the central stalk relative to the peripheral stalk (Yanagisawa and Frasch, 2017; Sielaff et al., 2019). The proportional differences of efficiencies of TD formation are shown relative to the average low efficiency for WT (Figure 3). Medium and high efficiency distributions of TDs in WT increased 1.5-fold and 2.2-fold, respectively, relative to low efficiency. The aN214L mutation increased the percent of TDs per data set for high, medium, and low efficiencies by 3-fold, 2-fold, and 1.2-fold, respectively, from the WT low efficiency. Mutations aQ252L and aE219L also increased TDs per data set for the high (2.7-fold and 2.5-fold) and medium (1.7-fold and 1.6-fold), but not the low efficiency distributions. Mutations aH245L and aE196L either did not increase the efficiency or slightly decreased the efficiency of the distributions of TD formation per data set.
Synthase-direction steps rotate CW an average of ~11°
The proportion of TDs with and without a synthase-direction step for WT and mutants are shown in Figure 4A at the pH values when the proportion of synthase-direction steps was minimum (black bars) and maximum (red bars), and at all pH values examined in Figure 4—figure supplement 1. The minimum proportion of synthase-direction steps was observed at pH 5.5 for WT and all mutants except aN214L that occurred at pH 6.0. Even at these low pH values, synthase-direction steps accounted for 62–68% of all TDs. In WT, a maximum of ~80% of TDs contained synthase-direction steps at pH 7.0, which was an increase of 13% from the minimum. These plots also show the distributions of the extent of CW rotation during a synthase-direction step, for which the 11° and 9° average and median values of CW rotation, respectively, were not changed significantly by the mutations (Figure 4B).
After subtracting the occurrence of the extent of synthase-direction step CW rotation at the pH when it was at a minimum (black bars) from that observed at other pH values (Figure 4—figure supplement 1) including that at its maximum (red bars), a Gaussian distribution of the increase in the extent of synthase-direction step CW rotation was observed (Figure 4C). During a synthase-direction step, the mean and standard deviations in the extent of CW rotation (Figure 4—figure supplement 2) was 12° ± 3° for WT, with little variation resulting from the mutations including: 11° ± 3° (aN214L), 11° ± 4° (aQ252L), 11° ± 3° (aH245L), 10° ± 3° (aE219L), and 11° ± 3° (aE196L) . In all cases, the distributions were truncated with minimum CW rotational steps of 6°. At their maxima, the extents of CW c-ring rotation during synthase events rotated 25° and 36° about 1% and 0.1% of the time, respectively.
Subunit-a mutations affect the proportion of TDs with synthase-direction steps
We tested the hypothesis that synthase-direction steps result when both the input and output channels are in the correct protonation state to enable H+ transfer to the carboxyl groups of the leading, and from the lagging c-subunits, respectively. If correct, the pH dependence of synthase-direction steps should follow the sum of the proportions of the protonated input and unprotonated output channels even when these pKa values change as the result of mutations. Alternatively, it was possible that the 11° rotations that we attribute to synthase-direction steps resulted instead from twisting of subunit-a and the c-ring as a single unit in response to the mismatch of the c-ring and catalytic dwell positions, which would not be subject to changes by subunit-a mutations.
The subset of TDs that forced the c-ring to rotate CW (synthase-direction steps) against the CCW force of F1-ATPase rotation was pH dependent (Figure 4D). A maximum of 80% of TDs contained synthase-direction steps in WT at ~pH 7.3, and a minimum of 67% at pH 5.5. At pH values > 7.5, the proportion of synthase-direction steps decreased to 71% at pH 9.0.
Because a TD either contains (TS) or lacks (TN) a synthase-direction step, the pH dependence of TDs with a synthase-direction step (Figure 4D) was the inverse of that without a synthase-direction step (Figure 4E) per Equation 2.
For WT, the minimum TN of 20% at pH 7.5 increased 1.7-fold and 1.5-fold at pH 5.5 and at pH 9.0, respectively. At these extremes of low and high pH values, TD formation was dominated by groups where either pKa1 is protonated or by unprotonated groups with pK2. This conclusion is supported by the good fits of the pH dependencies of TDs without synthase-direction steps for WT and subunit-a mutants (Figure 4E) to Equation 3, where the probability of forming a TD without a synthase-direction step (TN) is the sum of the probability (P1) of the protonated group(s) with pKa1 (X1), and the probability (P2) of unprotonated group(s) with pKa2 (Y2). Thus, these results support the conclusion that a TD without a synthase-direction step can result from a H+ transfer event from the protonated group with pKa1 or from a H+ transfer event to the unprotonated group with pKa2.
Fits of the pH dependence of TDs without synthase-direction steps from Equation 3 (black line) were based on the pKa values (Figure 4E), and probabilities summarized in Table 2. The WT data fit to probabilities of 38% and 33% for protonated groups (pKa1 6.5) and unprotonated groups (pKa2 7.7), respectively, such that the difference between the pKa values was 1.2 pH units. Consequently, TN showed a minimum at ~pH 7.3, and maxima at high and low pH values when only the group(s) with either pKa1 or pKa2 were protonated and unprotonated, respectively.
As a result of the subunit-a mutations, P1 values changed to a smaller extent (32–42%) than did P2 values (5–35%). Except for aE219L, all mutations decreased P2, including a >6-fold decrease with aN214L. The difference between pKa values observed with the mutants was from 0.3 to 0.5 pH units compared to the 1.2 pH unit difference of WT. Both pKa1 and pKa2 of aN214L increased by 1.5 and 0.7 pH units such that the minimum TN of ~18% at pH 8.0 represented an increase of 0.7 pH units from that of WT. At pH 5.5, TNs comprised 38% of all TDs in aN214L. A similar but smaller shift of the minimum TN occurrence to pH 7.5 was also observed for aH245L, which primarily resulted from an increase in pKa1 by 0.8 pH units from WT. A striking effect of mutations aQ252L, aE219L, and aE196L was that they suppressed the pH dependence of synthase-direction step formation. Of these, aE219L was most suppressed where TS varied between 66% and 71% of TDs over the pH range examined.
In all cases, the occurrence of synthase-direction steps reached a maximum at the crossover point between the fractions of protonated groups with pKa1 and unprotonated groups with pKa2. This is the point at which the largest fractions of both groups were in the correct protonation state where H+ transfer events could occur from the pKa1 groups to the c-ring, and from the c-ring to the pKa2 groups. These results eliminate the alternative hypothesis that synthase-direction steps result from the twisting of FO relative to F1 as the result of elastic energy from c10-ring:catalytic dwell mismatches because the elastic energy resulting from the mismatch in rotary positions would not be affected by these mutations.
Discussion
The results presented here provide new insight into the mechanism by which the FO motor uses the energy from H+ translocation to generate CW rotational torque on the c-ring to catalyze ATP synthesis. These studies support the hypothesis that residues associated with the half-channels work together to support a water column that transfers protons across the membrane in a coherent manner coupled to c-ring rotation in lieu of transferring protons directly and independently. These single-molecule investigations also tested the hypothesis that c-rotation operates via a Brownian ratchet versus a power stroke mechanism, and the results provide the first direct evidence that synthase-direction steps can occur by both mechanisms. Finally, the results presented here show that the proton translocation-dependent synthase-direction rotation occurs in 11° steps. These results do not support the hypothesis that the function of the essential aR210 is to deprotonate cD61 because recent FO structures show that the unprotonated lagging cD61 carboxyl is still 7.3 Å away from aR210 after an 11° c-ring rotation. Alternatively, an alternating two-step mechanism is proposed below to resolve this discrepancy.
FO uses a Grotthuss mechanism to translocate protons through both half-channels
The results presented here support a Grotthuss mechanism in FO where water columns in each half-channel communicate via rotation-dependent H+ transfer to and from the leading and lagging c-ring cD61 carboxyls. The coherent behavior of the water columns enables the release of a H+ to the cytoplasm concurrent with each H+ that enters the subunit-a input channel from the periplasm. This conclusion is supported by observations that: (i) ATP synthase-direction steps were maximal in WT and mutants when the fractions of protonated groups and unprotonated groups with low and high pKa values, respectively, were optimal for H+ transfer both from the lagging cD61 to the output channel and from the leading cD61 to the input channel; (ii) mutation of a residue from either the input or output channel altered both the low and high pKa values of TDs indicating that the channels communicate; (iii) all mutations changed the ability to form TDs, indicating all the groups examined participate; and (iv) none of the mutations completely eliminated the ability to form TDs.
This conclusion is also consistent with the fact that participating residues aS199, aN214, and aQ252 are polar but not ionizable, and distances between channel residues that are too far apart for direct H+ transfer but are positioned at distances able to support a water column. Although aQ252 is highly conserved (Figure 5—figure supplement 1), glycine or hydrophobic groups are naturally substituted: (i) for aN214 (Toxoplasma gondii, Thermus thermophilus); (ii) for aH245 (Mycobacterium phlei, Tetrahymena thermophila, Acetobacter woodii, Toxoplasma gondii, Ilyobacter tartaricus, Fusobacterium nucleatum, Thermus thermophilus); and (iii) for aE219 (Tetrahymena thermophila, Euglena gracilis, Bacillus pseudofirmus OF4, Pichia angusta, Saccharomyces cerevisiae, Arthrospira platensis PCC9438). Such mutations can be tolerated if the primary role of these groups was to support a water column that transferred protons.
A recent F1FO structure from bovine mitochondria was of sufficient resolution to observe density near the input channel residues consistent with Grotthuss-type water molecules in this half-channel (Spikes et al., 2020). Unidentified electron densities near subunit-a input channel residues in F1FO structures from E. coli (Sobti et al., 2020) and from Polytomella (Murphy et al., 2019) also suggest the presence of bound waters. The observation of a water column in both half-channels of the VO complex (Roh et al., 2020) suggests that Grotthuss-based H+ translocation is a commonly shared trait among the greater family of rotary ATPases.
Additional structural evidence that supports the existence of a Grotthuss H+ translocation mechanism (Figure 5) is the presence of an ~30 Å diameter funnel that is lined with carboxylate and imidazole residues as the funnel narrows (Figure 5). The aE219-carboxyl examined here, which we propose to be the start of the Grotthuss column is positioned at the apex of this funnel. A Grotthuss mechanism was first proposed to explain extremely high rates of FO-dependent H+ translocation across R. capsulatus membranes (Feniouk et al., 2004). The rates were so fast that an ~40 Å diameter Coulomb cage of charged and polar groups was proposed to be required to serve as a H+ ‘antenna’ to increase the delivery rate of protons from the aqueous solution to the entrance of the input channel water column (Wraight, 2006).
The results presented here indicate that multiple residues contribute to the pKa values that enable synthase-direction steps. The funnel extends the input channel from the residues investigated here that are located near the middle of the membrane to the periplasm, which provides the ultimate supply of protons to the c-ring. The charged and polar residues in the funnel are likely to contribute to the pKa value, and may also facilitate the protonation of the input channel relative to the output channel even though both half-channels are exposed to the same buffer as the result of the lipid bilayer nanodiscs employed. More work is necessary to address this issue.
FO undergoes H+ translocation-dependent 11° c-ring synthase-direction rotation steps
The extent of rotation during ATP synthase-direction steps was unexpectedly found to rotate CW by 11° in the WT and all mutants. Evidence presented here supports the conclusion that synthase-direction steps result from protonation of the leading cD61 from the input channel and from deprotonation of the lagging cD61 to the output channel to rotate the c-ring relative to subunit-a. These results include that: (i) synthase-direction steps depend on a group of residues with a low pKa that must be protonated, and a second group with a high pKa that must be unprotonated; (ii) formation of synthase-direction steps reached a maximum at the pH when the fractions of protonated groups with the low pKa and unprotonated groups with the high pKa were optimal; and (iii) mutating subunit-a residues in either the input or output half-channels altered both high and low pKa values, and altered the extent and pH dependence of synthase-direction step formation. The effects of the subunit-a mutants on the synthase-direction steps rule out the possibility that these steps result from twisting the entire FO relative to F1.
During a CCW F1-ATPase power stroke, TDs occur every 36°, which is equivalent to an interaction between subunit-a and each successive c-subunit in the E. coli c10-ring. Since synthase-direction steps rotate by 11°, rotation by an additional 25° is required to advance the c-ring by one full c-subunit, which we observed in only 0.1% of the synthase-direction steps. Rotational sub-state structures (pdb-IDs 6OQR and 6OQS) of E. coli F1FO that differ by a 25° rotation of the c-ring relative to subunit-a were obtained by cryo-EM (Sobti et al., 2020). Since advancing the c-ring by one c-subunit involves rotation by 36°, the difference between these sub-state structures also reveals information relevant to the 11° sub-step reported here.
The E. coli F1FO rotational sub-state structures that differ by the 25° rotation of the c-ring relative to subunit-a were obtained when the complex was inhibited by ADP (Sobti et al., 2020). Similar 11° and 25° rotational sub-states have also been observed with ADP-inhibited F1FO from B. taurus (Zhou et al., 2015) and from M. smegmatis (Guo et al., 2021). In M. smegmatis F1FO, the binding of bedaquiline stabilizes a rotational sub-state that is either 25° CW or 8° CCW from the equivalent rotational state in the absence of the drug (Guo et al., 2021). The rotational position of the c-ring in the cryo-EM structure of S. cerevisiae F1FO is also changed by ~9° when the inhibitor oligomycin is bound to FO (Srivastava et al., 2018).
Several structural features of E. coli FO (Sobti et al., 2020) are relevant to its ability to undergo synthase-direction steps, and in combination with the results presented here, they provide insight into the mechanism of sustained CW rotation to power ATP synthesis. A transmembrane view of subunit-a (Figure 5) shows that aE196, aS199, aR201, aN214, and aQ252 are aligned along the plane of cD61 rotation. This plane is surrounded by hydrophobic residues that form a vestibule. Between aS199 and aR210, polar groups line the vestibule above (aS202 and aS206) and below (aK203 and aY263) the cD61 rotation plane. Although these polar groups do not directly participate in H+ translocation (Eya et al., 1991), they enable water to access the vestibule (Angevine and Fillingame, 2003; Angevine et al., 2003) to make it less hydrophobic than the lipid phase of the membrane. Residues that provide a possible path for the output channel from aE196 to the cytoplasm include aQ181, aE177, and the subunit-a C-terminal carboxyl, which span this distance at ~4 Å intervals (Sobti et al., 2020), consistent with that needed to stabilize a Grotthuss water channel. More work is required to characterize this channel, especially since aE196 and aS199 are the only output channel residues conserved among other species (Figure 5—figure supplement 1).
A mechanism where FO uses alternating 11° (Figure 6A,B) and 25° (Figure 6B,C) sub-steps to power c-ring rotation that drives ATP synthesis is consistent with the data presented here, and with E. coli F1FO structures 6OQR and 6OQS (Sobti et al., 2020). In 6OQS (Figure 6A), the lagging cD61 (orange) is 3.5 Å from aS199, which enables H+ transfer to aS199 and aE196 via bound water. The leading cD61 (pink) is 3.8 Å from the aR210-guanidinium, consistent with intervening water. This cD61 is also proximal to aN214 and aQ252, which positions it for protonation from the input channel via bound water. In our model, the pH-dependent 11° sub-step (Figure 6A) occurs upon H+ transfer from water bound to aN214 and aQ252 to the leading cD61, and H+ transfer from the lagging cD61 to aS199 and aE196.
The low, medium, and high efficiencies of TD formation reported here (Figure 2B) were attributed to torsional strain resulting from the asymmetry between 36° c10-ring stepping, and the 120 ° F1 power strokes (Yanagisawa and Frasch, 2017; Sielaff et al., 2019). Based on this asymmetry observed in ADP-inhibited F1FO structures (Sobti et al., 2016), high efficiency TD formation was proposed to occur (Sielaff et al., 2019) in the rotary state comparable to that in which rotary sub-state structures PDB-IDs 6OQR and 6OQS were subsequently observed at 3.1 Å resolution (Sobti et al., 2020). Sobti et al., 2020, concurred that torsional strain contributed to their ability to resolve the 6OQR and 6OQS sub-state structures. However, in results presented here, catalytically active F1Fo in lipid bilayer nanodiscs show successive 11° ATP synthase-direction steps every 36° including at the rotary position of the ATPase power stroke where ADP inhibits rotation (Figure 2A, dashed line). Because ATP synthase-direction steps can also occur with low efficiency when torsional strain decreases the probability of forming a synthase-direction step, it is clear that torsional strain is not the primary contributing factor to the ability of FO to undergo 11° ATP synthase-direction steps.
After the 11° sub-step (Figure 6B), the distance of the now unprotonated lagging cD61 to the aR210-guanidinium decreased from ~11.5 to ~7.5 Å. These distances are inconsistent with the long-held belief that the role of aR210 is to displace the proton from the lagging cD61. Instead, we postulate that the electrostatic attraction between the negatively charged lagging cD61 and aR210 is sufficient to induce the 25° sub-step. As the result of this sub-step, the distance between these groups decreases from 7.5 to 3.8 Å (Figure 6C). As the 11° sub-step repeats, the subsequent loss of negative charge when the lagging cD61 is protonated by aN214 and aQ252 then allows this c-subunit to rotate away from aR210 into the lipid bilayer.
The probability that a TD occurs may appear to be stochastic. However, its occurrence depends on the kinetics and the energetics of the system. Slowing the angular velocity of the F1 ATPase-driven power stroke increases TD occurrence at pH 8.0 (shown here to be suboptimal) indicating that the ability to form a TD depends on the rate that an interaction can form between subunit-a and each c-subunit relative to the angular velocity of F1-ATPase-driven rotation (Ishmukhametov et al., 2010).
Evidence supports the hypothesis that the energy for F1-ATPase power strokes is derived from ATP binding-dependent closure of the β-subunit lever domain upon subunit-γ, which is initiated at ~36° after the catalytic dwell in E. coli F1 (Martin et al., 2014). Based on the KD of ATP at 36° measured in Geobacillus stearothermophilus F1, the energy available for the power stroke from ATP binding is ~13.5 kBT (Adachi et al., 2012).
The results here suggest that the energy required to power the 11° synthase-direction step is close to that of the F1-ATPase power stroke including: (i) that some synthase-direction steps oscillate consistent with a Brownian ratchet, especially those steps that occur late in the F1 power stroke when the affinity for ATP is the highest; and (ii) that TD formation efficiency is increased or decreased (high and low efficiencies) by the 0.4 kBT of torsional energy from the ±14° rotary mismatch between F1 and FO calculated from Equation 4, where θ is the rotational displacement in radians using the spring constant, κ, of 12.6 kBT radian–2 measured for E. coli F1FO(Sielaff et al., 2008).
Consequently, the FO motor must have at least 13.5 kBT available to cause a TD. Possible sources of energy for TDs in addition to the 0.4 kBT of torsional energy include: (i) as much as 4.4 kBT from the difference of pKa values between the input and output channels; (ii) 5.9 kBT from the exclusion of the lagging charged cD61 carboxyl from the lipid bilayer into the aqueous vestibule, based on its measured desolvation energy (White and Wimley, 1999). The energy penalty of 0.8 kBT to insert the leading protonated cD61 carboxyl into the lipid bilayer is avoided by its conversion to the closed and locked position in the c-ring (Pogoryelov et al., 2010); and (iii) as much as 38.1 kBT from the electrostatic attraction of aR210 to unprotonated lagging cD61 at a distance of 11.5 Å when it is exposed directly to the lipid bilayer. The energy of this attractive force, which is highly dependent on the hydrophobicity of its environment, is calculated by the modified Coulomb equation (Equation 5), where qi and qj are elementary charges, rij is the interatomic distance (in Ångstroms), and e is the dielectric constant, which is a measure of the hydrophobicity of the environment that ranges from 2 (lipid bilayer) to 80 (aqueous solvent). Although we do not yet know how wet the vestibule is, a dielectric constant of 13 and an 11.5 Å aR210-cD61 distance results in 3.8 kBT, which when summed with the other energy sources totals 14.1 kBT without input of torsional energy. Since F1-ATPase rotation from the catalytic dwell to the point that ATP binds is powered by no more than 4 kBT (Martin et al., 2018), this explains why synthase-direction rotation at these rotational positions typically has power stroke characteristics.
After the 25° rotation step when the unprotonated cD61-aR210 distance is 3.8 Å, the electrostatic force is 11.4 kBT or 73.8 kBT when the dielectric constant is 13 or 2, respectively. Thus, the electrostatic interaction in a hydrophobic environment would be far too strong for any rotation to occur. More work is required to quantify the energetics of these sub-steps in the ATP synthesis mechanism because, when understood in combination with the steady-state pmf values and the dissociation constants of ATP, ADP, and Pi versus rotary position, these energy contributions will determine the non-equilibrium ATP/ADP•Pi concentration ratio that can be maintained by F1FO at steady-state in vivo.
Materials and methods
Mutagenesis and purification of nanodisc FOF1
Request a detailed protocolThe FO subunit-a mutant plasmids were constructed from the pNY1-Ase plasmid containing the entire unc operon with a 6-His tag on the N-terminus of subunit-β and a cysteine inserted at the second position of subunit c (c2∇Cys), described previously by Ishmukhametov et al., 2010. The aN214L, aQ252L, aH245L, aE219L, and aE196L mutations were made on the plasmid using QuikChangesII-XL Site-Directed Mutagenesis Kit (Agilent). E. coli strain DK8, which lacks the unc operon, was transformed with the mutant plasmid. Cells were grown in 8 L of LB medium containing 50 μg/mL of ampicillin while shaking at 170 rpm at 37°C. About 40 g wet weight of cells was pelleted by centrifugation at 7700× g for 10 min at 4°C and stored at –80°C.
All subsequent steps were carried out at 4°C. The cell pellet was thawed, resuspended in 40 mL of French Press Buffer containing 200 mM Tris-HCl (pH 8.0), 100 mM KCl, 5 mM MgCl2, 0.1 mM EDTA, and 2.5% (v/v) glycerol. The cells were lysed by running through a French Press at 12,000 psi, twice. Unbroken cells and cell debris were pelleted down at 7700× g for 15 min. The supernatant was transferred to an ultracentrifuge tube and the cell membrane was pelleted down by centrifugation at 180,000× g for 3 hr. To detergent-solubilize the FOF1, the membrane pellet was resuspended in Extraction buffer containing 6% (v/v) glycerol, 50 mM Tris-HCl (pH 8.0), 100 mM NaCl, 40 mM ε-aminocaproic acid, 15 mM p-aminobenzamidine, 1% octyl glucopyranoside, 0.5% sodium deoxycholate, 0.5% sodium cholate, 0.03% phosphatidylcholine, 30 mM imidazole, and 5 mM MgCl2. The membrane suspension was incubated while gently shaking for 90 min and ultracentrifuged at 180,000× g for 2 hr. The supernatant was loaded into Ni-NTA column containing 1.5 mL Ni-NTA slurry that was equilibrated with Extraction buffer. The column was washed with 20 mL of Extraction buffer, and protein was eluted with 3 mL of Extraction buffer containing 200 mM imidazole. The protein concentration of the Ni-NTA column elution was determined by BCA assay and FOF1 was incorporated into lipid bilayer nanodiscs.
To form nanodiscs, membrane scaffold protein MSP-1E3D1 was used. The MSP was prepared by removing the His-tag with overnight TEV protease digestion at room temperature (at 25:1 ratio, w/w). Cleaved MSP was purified by passing through a Ni-NTA column. To incorporate the FOF1 into nanodiscs, a 1:4 molar ratio of FOF1:MSP was mixed with the addition of 10% stock sodium cholate in Extraction buffer to make the final sodium cholate concentration of 1%. To biotinylate the cysteine residue inserted at the N-terminus of subunit c, a 10-fold molar excess of biotin maleimide was added to FOF1. The mixture was gently shaken for 15 min to form biotinylated FOF1 nanodisc. The sample was desalted by running through a Sephadex G-50 column equilibrated with Buffer A containing 6% (v/v) glycerol, 50 mM Tris-HCl (pH 8.0), 100 mM NaCl, 4 mM p-aminobenzamidine, and 5 mM MgCl2. The sample was aliquoted into 50 μL and quick frozen with liquid nitrogen until use. The presence of all F1FO subunits and the MSP in each preparation was confirmed by SDS-PAGE.
Gold-nanorod single-molecule experiments
Request a detailed protocolRotation of individual nanodisc FOF1 molecules were observed by single-molecule rotation assay. Sample slides were prepared with modifications of previously described methods (Yanagisawa and Frasch, 2017). Briefly, purified nanodisc FOF1 were immobilized on a microscope slide by the His-tag on subunit-β, unbound enzymes were washed off the slide with wash buffer (30 mM Tris, 30 mM PIPES, 10 mM KCl, at the appropriate pH), 80 × 40 nm AuNR coated with avidin was bound to the biotinylated c-ring of E. coli nanodisc FOF1, excess AuNRs were washed off with the wash buffer, and rotation buffer (1 mM Mg2+ ATP, 30 mM Tris, 30 mM PIPES, 10 mM KCl, at the pH indicated) was added to the slide. The rotation of individual molecules was observed by measuring the change in intensity of polarized red light scattered from the AuNR using a single-photon detector. In each molecule observed, the rotation of the nanorod attached to an active nanodisc FOF1 complex was confirmed by the change in scattered light intensity as a function of the rotational position of the polarizing filter as described previously (Spetzler et al., 2006; Hornung et al., 2011). To make the measurement of nanodisc FOF1 undergoing power strokes, the orientation of the polarizing filter was adjusted to align with the minimum light intensity position that that corresponded to one of the three catalytic dwells. The sinusoidal change of polarized red light intensity was measured as the AuNR rotated from 0° to 90° relative to the catalytic dwell position. Measurements were taken in the form of 5 s data set at frame rate of 100 kHz. The occurrence of TDs in each subunit-a mutant was analyzed at varying pH from 5.0 to 8.0. TDs that occurred during the power strokes in the recorded data sets were analyzed by determining the arcsin1/2 of the intensity at each time point (Martin et al., 2018; Martin et al., 2015 Sielaff et al., 2016).
Data availability
All data generated or analyzed during this study are included in the manuscript and supporting files. Data source files for all figures have been uploaded to DRYAD and can be located at: https://doi.org/10.5061/dryad.9cnp5hqhw.
-
Dryad Digital RepositoryData from: pH-dependent 11° F1FO ATP synthase sub-steps reveal insight into the FO torque generating mechanism.https://doi.org/10.5061/dryad.9cnp5hqhw
References
-
Aqueous access channels in subunit a of rotary ATP synthaseThe Journal of Biological Chemistry 278:6066–6074.https://doi.org/10.1074/jbc.M210199200
-
Interaction between Glu-219 and His-245 within the a subunit of F1F0-ATPase in Escherichia coliJournal of Biological Chemistry 263:6606–6612.https://doi.org/10.1016/S0021-9258(18)68684-3
-
Proton translocation by the F1F0ATPase of Escherichia coli Mutagenic analysis of the a subunitThe Journal of Biological Chemistry 264:3292–3300.https://doi.org/10.1016/S0021-9258(18)94065-2
-
BookFundamentals of Protein Structure and FunctionIn: Buxbaum E, editors. Enzyme Kinetics and Mechanism. London: Garland Science. pp. 325–364.https://doi.org/10.4324/9780203833575
-
Et tu, Grotthuss! and other unfinished storiesBiochimica et Biophysica Acta 1757:876–885.https://doi.org/10.1016/j.bbabio.2005.12.001
-
Role of the carboxyl terminal region of H(+)-ATPase (F0F1) a subunit from Escherichia coliArchives of Biochemistry and Biophysics 284:71–77.https://doi.org/10.1016/0003-9861(91)90265-k
-
ATP synthesis catalyzed by the ATP synthase of Escherichia coli reconstituted into liposomesEuropean Journal of Biochemistry 225:167–172.https://doi.org/10.1111/j.1432-1033.1994.00167.x
-
The activity of the ATP synthase from Escherichia coli is regulated by the transmembrane proton motive forceThe Journal of Biological Chemistry 275:30157–30162.https://doi.org/10.1074/jbc.275.39.30157
-
Second-site suppressor mutations at glycine 218 and histidine 245 in the alpha subunit of F1F0 ATP synthase in Escherichia coliThe Journal of Biological Chemistry 269:32313–32317.https://doi.org/10.1016/S0021-9258(18)31637-5
-
Microsecond resolution of single-molecule rotation catalyzed by molecular motorsMethods in Molecular Biology 778:273–289.https://doi.org/10.1007/978-1-61779-261-8_18
-
Mutational analysis of the function of the a-subunit of the F0F1-APPase of Escherichia coliBiochimica et Biophysica Acta 1015:264–268.https://doi.org/10.1016/0005-2728(90)90030-8
-
Structure and Mechanisms of F-Type ATP SynthasesAnnual Review of Biochemistry 88:515–549.https://doi.org/10.1146/annurev-biochem-013118-110903
-
The proton pore in the Escherichia coli F0F1-ATPase: a requirement for arginine at position 210 of the a-subunitBiochimica et Biophysica Acta 894:399–406.https://doi.org/10.1016/0005-2728(87)90118-6
-
Fo-driven Rotation in the ATP Synthase Direction against the Force of F1 ATPase in the FoF1 ATP SynthaseThe Journal of Biological Chemistry 290:10717–10728.https://doi.org/10.1074/jbc.M115.646430
-
How Fo-ATPase generates rotary torquePhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 355:523–528.https://doi.org/10.1098/rstb.2000.0593
-
Cryo-EM structure of the entire mammalian F-type ATP synthaseNature Structural & Molecular Biology 27:1077–1085.https://doi.org/10.1038/s41594-020-0503-8
-
Microscopic rotary mechanism of ion translocation in the F(o) complex of ATP synthasesNature Chemical Biology 6:891–899.https://doi.org/10.1038/nchembio.457
-
The uniqueness of subunit α of mycobacterial F-ATP synthases: An evolutionary variant for niche adaptationThe Journal of Biological Chemistry 292:11262–11279.https://doi.org/10.1074/jbc.M117.784959
-
Power Stroke Angular Velocity Profiles of Archaeal A-ATP Synthase Versus Thermophilic and Mesophilic F-ATP Synthase Molecular MotorsThe Journal of Biological Chemistry 291:25351–25363.https://doi.org/10.1074/jbc.M116.745240
-
BookEnergy transduction by the two molecular motors of the F1FOATP synthase.In: Spetzler D, editors. Photosynthesis: Plastid Biology, Energy Conversion and Carbon Assimilation, Advances in Photosynthesis and Respiration. Dordrecht: Springer. pp. 561–583.https://doi.org/10.1093/aob/mcs281
-
Mutagenesis of the alpha subunit of the F1Fo-ATPase from Escherichia coli Mutations at Glu-196, Pro-190, and Ser-199The Journal of Biological Chemistry 263:6599–6605.https://doi.org/10.1016/S0021-9258(18)68683-1
-
Membrane protein folding and stability: physical principlesAnnual Review of Biophysics and Biomolecular Structure 28:319–365.https://doi.org/10.1146/annurev.biophys.28.1.319
-
Chance and design--proton transfer in water, channels and bioenergetic proteinsBiochimica et Biophysica Acta 1757:886–912.https://doi.org/10.1016/j.bbabio.2006.06.017
-
Protonation-dependent stepped rotation of the F-type ATP synthase c-ring observed by single-molecule measurementsThe Journal of Biological Chemistry 292:17093–17100.https://doi.org/10.1074/jbc.M117.799940
Article and author information
Author details
Funding
National Institute of General Medical Sciences (R01GM097510)
- Wayne D Frasch
National Science Foundation (2119963)
- Wayne D Frasch
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was funded by NSF-BII 2119963, and from NIH R01GM097510 to WDF.
Copyright
© 2021, Yanagisawa and Frasch
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 981
- views
-
- 154
- downloads
-
- 17
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Structural Biology and Molecular Biophysics
Peroxisome proliferator-activated receptor gamma (PPARγ) is a nuclear receptor transcription factor that regulates gene expression programs in response to ligand binding. Endogenous and synthetic ligands, including covalent antagonist inhibitors GW9662 and T0070907, are thought to compete for the orthosteric pocket in the ligand-binding domain (LBD). However, we previously showed that synthetic PPARγ ligands can cooperatively cobind with and reposition a bound endogenous orthosteric ligand to an alternate site, synergistically regulating PPARγ structure and function (Shang et al., 2018). Here, we reveal the structural mechanism of cobinding between a synthetic covalent antagonist inhibitor with other synthetic ligands. Biochemical and NMR data show that covalent inhibitors weaken—but do not prevent—the binding of other ligands via an allosteric mechanism, rather than direct ligand clashing, by shifting the LBD ensemble toward a transcriptionally repressive conformation, which structurally clashes with orthosteric ligand binding. Crystal structures reveal different cobinding mechanisms including alternate site binding to unexpectedly adopting an orthosteric binding mode by altering the covalent inhibitor binding pose. Our findings highlight the significant flexibility of the PPARγ orthosteric pocket, its ability to accommodate multiple ligands, and demonstrate that GW9662 and T0070907 should not be used as chemical tools to inhibit ligand binding to PPARγ.
-
- Structural Biology and Molecular Biophysics
SCARF1 (scavenger receptor class F member 1, SREC-1 or SR-F1) is a type I transmembrane protein that recognizes multiple endogenous and exogenous ligands such as modified low-density lipoproteins (LDLs) and is important for maintaining homeostasis and immunity. But the structural information and the mechanisms of ligand recognition of SCARF1 are largely unavailable. Here, we solve the crystal structures of the N-terminal fragments of human SCARF1, which show that SCARF1 forms homodimers and its epidermal growth factor (EGF)-like domains adopt a long-curved conformation. Then, we examine the interactions of SCARF1 with lipoproteins and are able to identify a region on SCARF1 for recognizing modified LDLs. The mutagenesis data show that the positively charged residues in the region are crucial for the interaction of SCARF1 with modified LDLs, which is confirmed by making chimeric molecules of SCARF1 and SCARF2. In addition, teichoic acids, a cell wall polymer expressed on the surface of gram-positive bacteria, are able to inhibit the interactions of modified LDLs with SCARF1, suggesting the ligand binding sites of SCARF1 might be shared for some of its scavenging targets. Overall, these results provide mechanistic insights into SCARF1 and its interactions with the ligands, which are important for understanding its physiological roles in homeostasis and the related diseases.