Theory for rates, equilibrium constants, and Brønsted slopes in F1-ATPase single molecule imaging experiments

Towards the understanding of molecular motors and its relationship with local unfolding

Molecular motors are machines essential for life since they convert chemical energy into mechanical work. However, the precise mechanism by which nucleotide binding, catalysis, or release of products is coupled to the work performed by the molecular motor is still not entirely clear. This is due, in part, to a lack of understanding of the role of force in the mechanical-structural processes involved in enzyme catalysis. From a mechanical perspective, one promising hypothesis is the Haldane-Pauling hypothesis which considers the idea that part of the enzymatic catalysis is strain-induced. It suggests that enzymes cannot be efficient catalysts if they are fully complementary to the substrates. Instead, they must exert strain on the substrate upon binding, using enzyme-substrate energy interaction (binding energy) to accelerate the reaction rate. A novel idea suggests that during catalysis, significant strain energy is built up, which is then released by a local unfolding/refolding event known as 'cracking'. Recent evidence has also shown that in catalytic reactions involving conformational changes, part of the heat released results in a center-of-mass acceleration of the enzyme, raising the possibility that the heat released by the reaction itself could affect the enzyme's integrity. Thus, it has been suggested that this released heat could promote or be linked to the cracking seen in proteins such as adenylate kinase (AK). We propose that the energy released as a consequence of ligand binding/catalysis is associated with the local unfolding/refolding events (cracking), and that this energy is capable of driving the mechanical work.

Angle-dependent rotation velocity consistent with ADP release in bacterial F1-ATPase

A model-based method is used to extract a short-lived state in the rotation kinetics of the F1-ATPase of a bacterial species, Paracoccus denitrificans (PdF1). Imaged as a single molecule, PdF1 takes large 120ø steps during it rotation. The apparent lack of further substeps in the trajectories not only renders the rotation of PdF1 unlike that of other F-ATPases, but also hinders the establishment of its mechano-chemical kinetic scheme. We addressed these challenges using the angular velocity extracted from the single-molecule trajectories and compare it with its theoretically calculated counterpart. The theory-experiment comparison indicate the presence of a 20μs lifetime state, 40o after ATP binding. We identify a kinetic cycle in which this state is a three-nucleotide occupancy state prior to ADP release from another site. A similar state was also reported in our earlier study of the Thermophilic bacillus F1-ATPase (lifetime ∼ 10 μ s), suggesting thereby a common mechanism for removing a nucleotide release bottleneck in the rotary mechanism.

Protein Dynamics and Enzymatic Catalysis

This Perspective presents a review of our work and that of others in the highly controversial topic of the coupling of protein dynamics to reaction in enzymes. We have been involved in studying this topic for many years. Thus, this perspective will naturally present our own views, but it also is designed to present an overview of the variety of viewpoints of this topic, both experimental and theoretical. This is obviously a large and contentious topic.

F1-ATPase Rotary Mechanism: Interpreting Results of Diverse Experimental Modes With an Elastic Coupling Theory

In this chapter, we review single-molecule observations of rotary motors, focusing on the general theme that their mechanical motion proceeds in substeps with each substep described by an angle-dependent rate constant. In the molecular machine F1-ATPase, the stepping rotation is described for individual steps by forward and back reaction rate constants, some of which depend strongly on the rotation angle. The rotation of a central shaft is typically monitored by an optical probe. We review our recent work on the theory for the angle-dependent rate constants built to treat a variety of single-molecule and ensemble experiments on the F1-ATPase, and relating the free energy of activation of a step to the standard free energy of reaction for that step. This theory, an elastic molecular transfer theory, provides a framework for a multistate model and includes the probe used in single-molecule imaging and magnetic manipulation experiments. Several examples of its application are the following: (a) treatment of the angle-dependent rate constants in stalling experiments, (b) use of the model to enhance the time resolution of the single-molecule imaging apparatus and to detect short-lived states with a microsecond lifetime, states hidden by the fluctuations of the imaging probe, (c) treatment of out-of-equilibrium "controlled rotation" experiments, (d) use of the model to predict, without adjustable parameters, the angle-dependent rate constants of nucleotide binding and release, using data from other experiments, and (e) insights obtained from correlation of kinetic and cryo-EM structural data. It is also noted that in the case where the release of ADP would be a bottleneck process, the binding of ATP to another site acts to accelerate the release by 5-6 orders of magnitude. The relation of the present set of studies to previous and current theoretical work in the field is described. An overall goal is to gain mechanistic insight into the biological function in relation to structure.

Thermodynamic marking of FOF1 ATP synthase

F_OF₁ ATP synthase is a ~100% efficient molecular machine for energy conversion in biology, and holds great lessons for man-made energy technology and nanotechnology. In light of formidable biocomplexity of the F_OF₁ machinery, its modeling from pure physical principles remains difficult and rare. Here we construct a thermodynamic model of F_OF₁ from experimentally accessible quantities plus a single entropy production that generally has vanishingly small values (<1k_B). Based on the physical inputs, this model captures F_OF₁ performance observed over an exhaustively wide range of proton-motive force and nucleotide concentrations. The model predicts a distinct 1/8k_BT slope for ATP synthesis rate versus proton-motive force, which is verified by experimental data and represents a profound thermodynamic marking of this amazingly efficient machine operating near a universal limit of the 2nd law of thermodynamics. The model further predicts two symmetries of heat productions, which are testable by available experimental techniques and offer quantitative constraints on F_OF₁'s possible mechanisms behind its ~100% efficiency.

Rotary catalysis of bovine mitochondrial F1-ATPase studied by single-molecule experiments

The reaction scheme of rotary catalysis and the torque generation mechanism of bovine mitochondrial F₁ (bMF₁) were studied in single-molecule experiments. Under ATP-saturated concentrations, high-speed imaging of a single 40-nm gold bead attached to the γ subunit of bMF₁ showed 2 types of intervening pauses during the rotation that were discriminated by short dwell and long dwell. Using ATPγS as a slowly hydrolyzing ATP derivative as well as using a functional mutant βE188D with slowed ATP hydrolysis, the 2 pausing events were distinctively identified. Buffer-exchange experiments with a nonhydrolyzable analog (AMP-PNP) revealed that the long dwell corresponds to the catalytic dwell, that is, the waiting state for hydrolysis, while it remains elusive which catalytic state short pause represents. The angular position of catalytic dwell was determined to be at +80° from the ATP-binding angle, mostly consistent with other F₁s. The position of short dwell was found at 50 to 60° from catalytic dwell, that is, +10 to 20° from the ATP-binding angle. This is a distinct difference from human mitochondrial F₁, which also shows intervening dwell that probably corresponds to the short dwell of bMF₁, at +65° from the binding pause. Furthermore, we conducted "stall-and-release" experiments with magnetic tweezers to reveal how the binding affinity and hydrolysis equilibrium are modulated by the γ rotation. Similar to thermophilic F₁, bMF₁ showed a strong exponential increase in ATP affinity, while the hydrolysis equilibrium did not change significantly. This indicates that the ATP binding process generates larger torque than the hydrolysis process.

Method to extract multiple states in F1-ATPase rotation experiments from jump distributions

A method is proposed for analyzing fast (10 μs) single-molecule rotation trajectories in F₁ adenosinetriphosphatase ([Formula: see text]-ATPase). This method is based on the distribution of jumps in the rotation angle that occur in the transitions during the steps between subsequent catalytic dwells. The method is complementary to the "stalling" technique devised by H. Noji et al. [Biophys. Rev. 9, 103-118, 2017], and can reveal multiple states not directly detectable as steps. A bimodal distribution of jumps is observed at certain angles, due to the system being in either of 2 states at the same rotation angle. In this method, a multistate theory is used that takes into account a viscoelastic fluctuation of the imaging probe. Using an established sequence of 3 specific states, a theoretical profile of angular jumps is predicted, without adjustable parameters, that agrees with experiment for most of the angular range. Agreement can be achieved at all angles by assuming a fourth state with an ∼10 μs lifetime and a dwell angle about 40° after the adenosine 5'-triphosphate (ATP) binding dwell. The latter result suggests that the ATP binding in one β subunit and the adenosine 5'-diphosphate (ADP) release from another β subunit occur via a transient whose lifetime is ∼10 μs and is about 6 orders of magnitude smaller than the lifetime for ADP release from a singly occupied [Formula: see text]-ATPase. An internal consistency test is given by comparing 2 independent ways of obtaining the relaxation time of the probe. They agree and are ∼15 μs.

Insights into the origin of the high energy-conversion efficiency of F1-ATPase

Our understanding of the rotary-coupling mechanism of F₁-ATPase has been greatly enhanced in the last decade by advances in X-ray crystallography, single-molecular imaging, and theoretical models. Recently, Volkán-Kacsó and Marcus [S. Volkán-Kacsó, R. A. Marcus, Proc. Natl. Acad. Sci. U.S.A. 112, 14230 (2015)] presented an insightful thermodynamic model based on the Marcus reaction theory coupled with an elastic structural deformation term to explain the observed γ-rotation angle dependence of the adenosine triphosphate (ATP)/adenosine diphosphate (ADP) exchange rates of F₁-ATPase. Although the model is successful in correlating single-molecule data, it is not in agreement with the available theoretical results. We describe a revision of the model, which leads to consistency with the simulation results and other experimental data on the F₁-ATPase rotor compliance. Although the free energy liberated on ATP hydrolysis by F₁-ATPase is rapidly dissipated as heat and so cannot contribute directly to the rotation, we show how, nevertheless, F₁-ATPase functions near the maximum possible efficiency. This surprising result is a consequence of the differential binding of ATP and its hydrolysis products ADP and P_i along a well-defined pathway.

Elastic coupling power stroke mechanism of the F1-ATPase molecular motor

The angular velocity profile of the 120° F₁-ATPase power stroke was resolved as a function of temperature from 16.3 to 44.6 °C using a Δμ_ATP = -31.25 k_BT at a time resolution of 10 μs. Angular velocities during the first 60° of the power stroke (phase 1) varied inversely with temperature, resulting in negative activation energies with a parabolic dependence. This is direct evidence that phase 1 rotation derives from elastic energy (spring constant, κ = 50 k_BT·rad^-2). Phase 2 of the power stroke had an enthalpic component indicating that additional energy input occurred to enable the γ-subunit to overcome energy stored by the spring after rotating beyond its 34° equilibrium position. The correlation between the probability distribution of ATP binding to the empty catalytic site and the negative E_a values of the power stroke during phase 1 suggests that this additional energy is derived from the binding of ATP to the empty catalytic site. A second torsion spring (κ = 150 k_BT·rad^-2; equilibrium position, 90°) was also evident that mitigated the enthalpic cost of phase 2 rotation. The maximum ΔG^ǂ was 22.6 k_BT, and maximum efficiency was 72%. An elastic coupling mechanism is proposed that uses the coiled-coil domain of the γ-subunit rotor as a torsion spring during phase 1, and then as a crankshaft driven by ATP-binding-dependent conformational changes during phase 2 to drive the power stroke.

What can be learned about the enzyme ATPase from single-molecule studies of its subunit F1?

We summarize the different types of single molecule experiments on the F1 component of FOF1-ATP Synthase and what has been learned from them. We also describe results from our recent studies on interpreting the experiments using a chemical-mechanical theory for these biological motors.

Inseparable tandem: evolution chooses ATP and Ca2+ to control life, death and cellular signalling

From the very dawn of biological evolution, ATP was selected as a multipurpose energy-storing molecule. Metabolism of ATP required intracellular free Ca(2+) to be set at exceedingly low concentrations, which in turn provided the background for the role of Ca(2+) as a universal signalling molecule. The early-eukaryote life forms also evolved functional compartmentalization and vesicle trafficking, which used Ca(2+) as a universal signalling ion; similarly, Ca(2+) is needed for regulation of ciliary and flagellar beat, amoeboid movement, intracellular transport, as well as of numerous metabolic processes. Thus, during evolution, exploitation of atmospheric oxygen and increasingly efficient ATP production via oxidative phosphorylation by bacterial endosymbionts were a first step for the emergence of complex eukaryotic cells. Simultaneously, Ca(2+) started to be exploited for short-range signalling, despite restrictions by the preset phosphate-based energy metabolism, when both phosphates and Ca(2+) interfere with each other because of the low solubility of calcium phosphates. The need to keep cytosolic Ca(2+) low forced cells to restrict Ca(2+) signals in space and time and to develop energetically favourable Ca(2+) signalling and Ca(2+) microdomains. These steps in tandem dominated further evolution. The ATP molecule (often released by Ca(2+)-regulated exocytosis) rapidly grew to be the universal chemical messenger for intercellular communication; ATP effects are mediated by an extended family of purinoceptors often linked to Ca(2+) signalling. Similar to atmospheric oxygen, Ca(2+) must have been reverted from a deleterious agent to a most useful (intra- and extracellular) signalling molecule. Invention of intracellular trafficking further increased the role for Ca(2+) homeostasis that became critical for regulation of cell survival and cell death. Several mutually interdependent effects of Ca(2+) and ATP have been exploited in evolution, thus turning an originally unholy alliance into a fascinating success story.This article is part of the themed issue 'Evolution brings Ca(2+) and ATP together to control life and death'.

Theory of long binding events in single-molecule-controlled rotation experiments on F1-ATPase

The theory of elastic group transfer for the binding and release rate constants for nucleotides in F₁-ATPase as a function of the rotor angle is further extended in several respects. (i) A method is described for predicting the experimentally observed lifetime distribution of long binding events in the controlled rotation experiments by taking into account the hydrolysis and synthesis reactions occurring during these events. (ii) A method is also given for treating the long binding events in the experiments and obtaining the rate constants for the hydrolysis and synthesis reactions occurring during these events. (iii) The theory in the previous paper is given in a symmetric form, an extension that simplifies the application of the theory to experiments. It also includes a theory-based correction of the reported "on" and "off" rates by calculating the missed events. A near symmetry of the data about the angle of -40° and a "turnover" in the binding rate data vs. rotor angle for angles greater than [Formula: see text]40° is also discussed.

Theory of single-molecule controlled rotation experiments, predictions, tests, and comparison with stalling experiments in F1-ATPase

A recently proposed chemomechanical group transfer theory of rotary biomolecular motors is applied to treat single-molecule controlled rotation experiments. In these experiments, single-molecule fluorescence is used to measure the binding and release rate constants of nucleotides by monitoring the occupancy of binding sites. It is shown how missed events of nucleotide binding and release in these experiments can be corrected using theory, with F₁-ATP synthase as an example. The missed events are significant when the reverse rate is very fast. Using the theory the actual rate constants in the controlled rotation experiments and the corrections are predicted from independent data, including other single-molecule rotation and ensemble biochemical experiments. The effective torsional elastic constant is found to depend on the binding/releasing nucleotide, and it is smaller for ADP than for ATP. There is a good agreement, with no adjustable parameters, between the theoretical and experimental results of controlled rotation experiments and stalling experiments, for the range of angles where the data overlap. This agreement is perhaps all the more surprising because it occurs even though the binding and release of fluorescent nucleotides is monitored at single-site occupancy concentrations, whereas the stalling and free rotation experiments have multiple-site occupancy.

Biophysical comparison of ATP synthesis mechanisms shows a kinetic advantage for the rotary process

The ATP synthase (F-ATPase) is a highly complex rotary machine that synthesizes ATP, powered by a proton electrochemical gradient. Why did evolution select such an elaborate mechanism over arguably simpler alternating-access processes that can be reversed to perform ATP synthesis? We studied a systematic enumeration of alternative mechanisms, using numerical and theoretical means. When the alternative models are optimized subject to fundamental thermodynamic constraints, they fail to match the kinetic ability of the rotary mechanism over a wide range of conditions, particularly under low-energy conditions. We used a physically interpretable, closed-form solution for the steady-state rate for an arbitrary chemical cycle, which clarifies kinetic effects of complex free-energy landscapes. Our analysis also yields insights into the debated "kinetic equivalence" of ATP synthesis driven by transmembrane pH and potential difference. Overall, our study suggests that the complexity of the F-ATPase may have resulted from positive selection for its kinetic advantage.

F1-ATPase conformational cycle from simultaneous single-molecule FRET and rotation measurements

Despite extensive studies, the structural basis for the mechanochemical coupling in the rotary molecular motor F1-ATPase (F1) is still incomplete. We performed single-molecule FRET measurements to monitor conformational changes in the stator ring-α3β3, while simultaneously monitoring rotations of the central shaft-γ. In the ATP waiting dwell, two of three β-subunits simultaneously adopt low FRET nonclosed forms. By contrast, in the catalytic intermediate dwell, two β-subunits are simultaneously in a high FRET closed form. These differences allow us to assign crystal structures directly to both major dwell states, thus resolving a long-standing issue and establishing a firm connection between F1 structure and the rotation angle of the motor. Remarkably, a structure of F1 in an ε-inhibited state is consistent with the unique FRET signature of the ATP waiting dwell, while most crystal structures capture the structure in the catalytic dwell. Principal component analysis of the available crystal structures further clarifies the five-step conformational transitions of the αβ-dimer in the ATPase cycle, highlighting the two dominant modes: the opening/closing motions of β and the loosening/tightening motions at the αβ-interface. These results provide a new view of tripartite coupling among chemical reactions, stator conformations, and rotary angles in F1-ATPase.

Thermodynamics and kinetics of the FoF1-ATPase: application of the probability isotherm

We use the results of recent publications as vehicles with which to discuss the thermodynamics of the proton-driven mitochondrial F o F1-ATP synthase, focusing particularly on the possibility that there may be dissociation between rotatory steps and ATP synthesis/hydrolysis. Such stoichiometric 'slippage' has been invoked in the literature to explain observed non-ideal behaviour. Numerical solution of the Rate Isotherm (the kinetic equivalent of the more fundamental Probability Isotherm) suggests that such 'slippage' is an unlikely explanation; instead, we suggest that the experimental results may be more consistent with damage to the enzyme caused by its isolation from the biomembrane and its experimental fixation, resulting in non-physiological friction within the enzyme's rotary mechanism. We emphasize the unavoidable constraint of the Second Law as instantiated by the obligatory dissipation of Gibbs Free Energy if the synthase is to operate at anything other than thermodynamic equilibrium. We use further numerical solution of the Rate Isotherm to demonstrate that there is no necessary association of low thermodynamic efficiency with high metabolic rates in a bio-world in which the dominating mechanism of metabolic control is multifactorial enzyme activation.

A Microscopic Capacitor Model of Voltage Coupling in Membrane Proteins: Gating Charge Fluctuations in Ci-VSD

The voltage sensitivity of membrane proteins is reflected in the response of the voltage sensing domains (VSDs) to the changes in membrane potential. This response is implicated in the displacement of positively charged residues, associated with the conformational changes of VSDs. The displaced charges generate nonlinear (i.e., voltage-dependent) capacitance current called the gating current (and its corresponding gating charge), which is a key experimental quantity that characterizes voltage activation in VSMP. However, the relevant theoretical/computational approaches, aimed to correlate the structural information on VSMP to electrophysiological measurements, have been rather limited, posing a broad challenge in computer simulations of VSMP. Concomitant with the development of our coarse-graining (CG) model of voltage coupling, we apply our theoretical framework for the treatments of voltage effects in membrane proteins to modeling the VSMP activation, taking the VSDs (Ci-VSD) derived from the Ciona intestinalis voltage sensitive phosphatase (Ci-VSP) as a model system. Our CG model reproduces the observed gating charge of Ci-VSD activation in several different perspectives. In particular, a new closed-form expression of the gating charge is evaluated in both nonequilibrium and equilibrium ways, while considering the fluctuation-dissipation relation that connects a nonequilibrium measurement of the gating charge to an equilibrium measurement of charge fluctuations (i.e., the voltage-independent linear component of membrane capacitance). In turn, the expression uncovers a novel link that connects an equilibrium measurement of the voltage-independent linear capacitance to a nonequilibrium measurement of the voltage-dependent nonlinear capacitance (whose integral over voltage is equal to the gating charge). In addition, our CG model yields capacitor-like voltage dependent free energy parabolas, resulting in the free energy difference and the free energy barrier for the Ci-VSD activation at "zero" (depolarization) membrane potential. Significantly, the resultant voltage dependent energetics enables a direct evaluation of capacitance-voltage relationship (C-V curve) as well as charge-voltage relationship (Q-V curve) that is in a good agreement with the observed measurement of Ci-VSD voltage activation. Importantly, an extension of our kinetic/thermodynamic model of voltage dependent activation in VSMP allows for novel derivations of voltage-dependent rate constants, whose parameters are expressed by the intrinsic properties of VSMP. These novel closed-form expressions offer a physicochemical foundation for the semiempirical Eyring-type voltage dependent rate equations that have been the cornerstone for the phenomenological (kinetic) descriptions of gating and membrane currents in the mechanistic study of ion channels and transporters. Our extended theoretical framework developed in the present study has potential implications on the roles played by gating charge fluctuations for the spike generations in nerve cells within the framework of the Hodgkin-Huxley-type model.