Mechanism of the conformational change of the F1-ATPase β subunit revealed by free energy simulations

Ordered-domain unfolding of thermophilic isolated β subunit ATP synthase

The flexibility of the ATP synthase's β subunit promotes its role in the ATP synthase rotational mechanism, but its domains stability remains unknown. A reversible thermal unfolding of the isolated β subunit (Tβ) of the ATP synthase from Bacillus thermophilus PS3, tracked through circular dichroism and molecular dynamics, indicated that Tβ shape transits from an ellipsoid to a molten globule through an ordered unfolding of its domains, preserving the β-sheet residual structure at high temperature. We determined that part of the stability origin of Tβ is due to a transversal hydrophobic array that crosses the β-barrel formed at the N-terminal domain and the Rossman fold of the nucleotide-binding domain (NBD), while the helix bundle of the C-terminal domain is the less stable due to the lack of hydrophobic residues, and thus the more flexible to trigger the rotational mechanism of the ATP synthase.

Revolving hexameric ATPases as asymmetric motors to translocate double-stranded DNA genome along one strand

DsDNA translocation through nanoscale pores is generally accomplished by ATPase biomotors. The discovery of the revolving dsDNA translocation mechanism, as opposed to rotation, in bacteriophage phi29 elucidated how ATPase motors move dsDNA. Revolution-driven, hexameric dsDNA motors have been reported in herpesvirus, bacterial FtsK, Streptomyces TraB, and T7 phage. This review explores the common relationship between their structure and mechanisms. Commonalities include moving along the 5'→3' strand, inchworm sequential action leading to an asymmetrical structure, channel chirality, channel size, and 3-step channel gating for controlling motion direction. The revolving mechanism and contact with one of the dsDNA strands addresses the historic controversy of dsDNA packaging using nicked, gapped, hybrid, or chemically modified DNA. These controversies surrounding dsDNA packaging activity using modified materials can be answered by whether the modification was introduced into the 3'→5' or 5'→3' strand. Perspectives concerning solutions to the controversy of motor structure and stoichiometry are also discussed.

Strategies for elucidation of the structure and function of the large membrane protein complex, FoF1-ATP synthase, by nuclear magnetic resonance

Nuclear magnetic resonance (NMR) investigation of large membrane proteins requires well-focused questions and critical techniques. Here, research strategies for F_oF₁-ATP synthase, a membrane-embedded molecular motor, are reviewed, focusing on the β-subunit of F₁-ATPase and c-subunit ring of the enzyme. Segmental isotope-labeling provided 89% assignment of the main chain NMR signals of thermophilic Bacillus (T)F₁β-monomer. Upon nucleotide binding to Lys164, Asp252 was shown to switch its hydrogen-bonding partner from Lys164 to Thr165, inducing an open-to-closed bend motion of TF₁β-subunit. This drives the rotational catalysis. The c-ring structure determined by solid-state NMR showed that cGlu56 and cAsn23 of the active site took a hydrogen-bonded closed conformation in membranes. In 505 kDa TF_oF₁, the specifically isotope-labeled cGlu56 and cAsn23 provided well-resolved NMR signals, which revealed that 87% of the residue pairs took a deprotonated open conformation at the F_oa-c subunit interface, whereas they were in the closed conformation in the lipid-enclosed region.

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.

In Silico Repurposing of J147 for Neonatal Encephalopathy Treatment: Exploring Molecular Mechanisms of Mutant Mitochondrial ATP Synthase

BACKGROUND

Neonatal Encephalopathy (NE) is a mitochondrial ATP synthase (mATPase) disease, which results in the death of infants. The case presented here is reportedly caused by complex V deficiency as a result of mutation of Arginine to Cysteine at residue 329 in the mATPase. A recent breakthrough was the discovery of J147, which targets mATPase in the treatment of Alzheimer's disease. Based on the concepts of computational target-based drug design, this study investigated the possibility of employing J147 as a viable candidate in the treatment of NE.

OBJECTIVE/METHODS

The structural dynamic implications of this drug on the mutated enzyme are yet to be elucidated. Hence, integrative molecular dynamics simulations and thermodynamic calculations were employed to investigate the activity of J147 on the mutated enzyme in comparison to its already established inhibitory activity on the wild-type enzyme.

RESULTS

A correlated structural trend occurred between the wild-type and mutant systems whereby all the systems exhibited an overall conformational transition. Equal observations in favorable free binding energies further substantiated uniformity in the mobility, and residual fluctuation of the wild-type and mutant systems. The similarity in the binding landscape suggests that J147 could as well modulate mutant mATPase activity in addition to causing structural modifications in the wild-type enzyme.

CONCLUSION

Findings suggest that J147 can stabilize the mutant protein and restore it to a similar structural state as the wild-type which depicts functionality. These details could be employed in drug design for potential drug resistance cases due to mATPase mutations that may present in the future.

Is F1-ATPase a Rotary Motor with Nearly 100% Efficiency? Quantitative Analysis of Chemomechanical Coupling and Mechanical Slip

We present a chemomechanical network model of the rotary molecular motor F₁-ATPase which quantitatively describes not only the rotary motor dynamics driven by ATP hydrolysis but also the ATP synthesis caused by forced reverse rotations. We observe a high reversibility of F₁-ATPase, that is, the main cycle of ATP synthesis corresponds to the reversal of the main cycle in the hydrolysis-driven motor rotation. However, our quantitative analysis indicates that torque-induced mechanical slip without chemomechanical coupling occurs under high external torque and reduces the maximal efficiency of the free energy transduction to 40-80% below the optimal efficiency. Heat irreversibly dissipates not only through the viscous friction of the probe but also directly from the motor due to torque-induced mechanical slip. Such irreversible heat dissipation is a crucial limitation for achieving a 100% free-energy transduction efficiency with biological nanomachines because biomolecules are easily deformed by external torque.

Multiscale molecular dynamics simulations of rotary motor proteins

Protein functions require specific structures frequently coupled with conformational changes. The scale of the structural dynamics of proteins spans from the atomic to the molecular level. Theoretically, all-atom molecular dynamics (MD) simulation is a powerful tool to investigate protein dynamics because the MD simulation is capable of capturing conformational changes obeying the intrinsically structural features. However, to study long-timescale dynamics, efficient sampling techniques and coarse-grained (CG) approaches coupled with all-atom MD simulations, termed multiscale MD simulations, are required to overcome the timescale limitation in all-atom MD simulations. Here, we review two examples of rotary motor proteins examined using free energy landscape (FEL) analysis and CG-MD simulations. In the FEL analysis, FEL is calculated as a function of reaction coordinates, and the long-timescale dynamics corresponding to conformational changes is described as transitions on the FEL surface. Another approach is the utilization of the CG model, in which the CG parameters are tuned using the fluctuation matching methodology with all-atom MD simulations. The long-timespan dynamics is then elucidated straightforwardly by using CG-MD simulations.

Mechanochemical Energy Transduction during the Main Rotary Step in the Synthesis Cycle of F1-ATPase

F₁-ATPase is a highly efficient molecular motor that can synthesize ATP driven by a mechanical torque. Its ability to function reversibly in either direction requires tight mechanochemical coupling between the catalytic domain and the rotating central shaft, as well as temporal control of substrate binding and product release. Despite great efforts and significant progress, the molecular details of this synchronized and fine-tuned energy conversion mechanism are not fully understood. Here, we use extensive molecular dynamics simulations to reconcile recent single-molecule experiments with structural data and provide a consistent thermodynamic, kinetic and mechanistic description of the main rotary substep in the synthetic cycle of mammalian ATP synthase. The calculated free energy profiles capture a discrete pattern in the rotation of the central γ-shaft, with a metastable intermediate located-consistently with recent experimental findings-at 70° relative to the X-ray position. We identify this rotary step as the ATP-dependent substep, and find that the associated free energy input supports the mechanism involving concurrent nucleotide binding and release. During the main substep, our simulations show no significant opening of the ATP-bound β subunit; instead, we observe that mechanical energy is transmitted to its nucleotide binding site, thus lowering the affinity for ATP. Simultaneously, the empty subunit assumes a conformation that enables the enzyme to harness the free energy of ADP binding to drive ATP release. Finally, we show that ligand exchange is regulated by a checkpoint mechanism, an apparent prerequisite for high efficiency in protein nanomotors.

Chemomechanical Coupling in Hexameric Protein-Protein Interfaces Harnesses Energy within V-Type ATPases

ATP synthase is the most prominent bioenergetic macromolecular motor in all life forms, utilizing the proton gradient across the cell membrane to fuel the synthesis of ATP. Notwithstanding the wealth of available biochemical and structural information inferred from years of experiments, the precise molecular mechanism whereby vacuolar (V-type) ATP synthase fulfills its biological function remains largely fragmentary. Recently, crystallographers provided the first high-resolution view of ATP activity in Enterococcus hirae V1-ATPase. Employing a combination of transition-path sampling and high-performance free-energy methods, the sequence of conformational transitions involved in a functional cycle accompanying ATP hydrolysis has been investigated in unprecedented detail over an aggregate simulation time of 65 μs. Our simulated pathways reveal that the chemical energy produced by ATP hydrolysis is harnessed via the concerted motion of the protein-protein interfaces in the V1-ring, and is nearly entirely consumed in the rotation of the central stalk. Surprisingly, in an ATPase devoid of a central stalk, the interfaces of this ring are perfectly designed for inducing ATP hydrolysis. However, in a complete V1-ATPase, the mechanical property of the central stalk is a key determinant of the rate of ATP turnover. The simulations further unveil a sequence of events, whereby unbinding of the hydrolysis product (ADP + Pi) is followed by ATP uptake, which, in turn, leads to the torque generation step and rotation of the center stalk. Molecular trajectories also bring to light multiple intermediates, two of which have been isolated in independent crystallography experiments.

Inherent conformational flexibility of F1-ATPase α-subunit

The core of F1-ATPase consists of three catalytic (β) and three noncatalytic (α) subunits, forming a hexameric ring in alternating positions. A wealth of experimental and theoretical data has provided a detailed picture of the complex role played by catalytic subunits. Although major conformational changes have only been seen in β-subunits, it is clear that α-subunits have to respond to these changes in order to be able to transmit information during the rotary mechanism. However, the conformational behavior of α-subunits has not been explored in detail. Here, we have combined unbiased molecular dynamics (MD) simulations and calorimetrically measured thermodynamic signatures to investigate the conformational flexibility of isolated α-subunits, as a step toward deepening our understanding of its function inside the α3β3 ring. The simulations indicate that the open-to-closed conformational transition of the α-subunit is essentially barrierless, which is ideal to accompany and transmit the movement of the catalytic subunits. Calorimetric measurements of the recombinant α-subunit from Geobacillus kaustophilus indicate that the isolated subunit undergoes no significant conformational changes upon nucleotide binding. Simulations confirm that the nucleotide-free and nucleotide-bound subunits show average conformations similar to that observed in the F1 crystal structure, but they reveal an increased conformational flexibility of the isolated α-subunit upon MgATP binding, which might explain the evolutionary conserved capacity of α-subunits to recognize nucleotides with considerable strength. Furthermore, we elucidate the different dependencies that α- and β-subunits show on Mg(II) for recognizing ATP.

Molecular Mechanism of Processive 3' to 5' RNA Translocation in the Active Subunit of the RNA Exosome Complex

Recent experimental studies revealed structural details of 3' to 5' degradation of RNA molecules, performed by the exosome complex. ssRNA is channeled through its multisubunit ring-like core into the active site tunnel of its key exonuclease subunit Rrp44, which acts both as an enzyme and a motor. Even in isolation, Rrp44 can pull and sequentially cleave RNA nucleotides, one at a time, without any external energy input and release a final 3-5 nucleotide long product. Using molecular dynamics simulations, we identify the main factors that control these processes. Our free energy calculations reveal that RNA transfer from solution into the active site of Rrp44 is highly favorable, but dependent on the length of the RNA strand. While RNA strands formed by 5 nucleotides or more correspond to a decreasing free energy along the translocation coordinate toward the cleavage site, a 4-nucleotide RNA experiences a free energy barrier along the same direction, potentially leading to incomplete cleavage of ssRNA and the release of short (3-5) nucleotide products. We provide new insight into how Rrp44 catalyzes a localized enzymatic reaction and performs an action distributed over several RNA nucleotides, leading eventually to the translocation of whole RNA segments into the position suitable for cleavage.

Torque transmission mechanism via DELSEED loop of F1-ATPase

F1-ATPase (F1) is an ATP-driven rotary motor in which the three catalytic β subunits in the stator ring sequentially induce the unidirectional rotation of the rotary γ subunit. Many lines of evidence have revealed open-to-closed conformational transitions in the β subunit that swing the C-terminal domain inward. This conformational transition causes a C-terminal protruding loop with conserved sequence DELSEED to push the γ subunit. Previous work, where all residues of DELSEED were substituted with glycine to disrupt the specific interaction with γ and introduce conformational flexibility, showed that F1 still rotated, but that the torque was halved, indicating a remarkable impact on torque transmission. In this study, we conducted a stall-and-release experiment on F1 with a glycine-substituted DELSEED loop to investigate the impact of the glycine substitution on torque transmission upon ATP binding and ATP hydrolysis. The mutant F1 showed a significantly reduced angle-dependent change in ATP affinity, whereas there was no change in the equilibrium for ATP hydrolysis. These findings indicate that the DELSEED loop is predominantly responsible for torque transmission upon ATP binding but not for that upon ATP hydrolysis.

CHARMM Force-Fields with Modified Polyphosphate Parameters Allow Stable Simulation of the ATP-Bound Structure of Ca(2+)-ATPase

Adenosine triphosphate (ATP) is an indispensable energy source in cells. In a wide variety of biological phenomena like glycolysis, muscle contraction/relaxation, and active ion transport, chemical energy released from ATP hydrolysis is converted to mechanical forces to bring about large-scale conformational changes in proteins. Investigation of structure-function relationships in these proteins by molecular dynamics (MD) simulations requires modeling of ATP in solution and ATP bound to proteins with accurate force-field parameters. In this study, we derived new force-field parameters for the triphosphate moiety of ATP based on the high-precision quantum calculations of methyl triphosphate. We tested our new parameters on membrane-embedded sarcoplasmic reticulum Ca(2+)-ATPase and four soluble proteins. The ATP-bound structure of Ca(2+)-ATPase remains stable during MD simulations, contrary to the outcome in shorter simulations using original parameters. Similar results were obtained with the four ATP-bound soluble proteins. The new force-field parameters were also tested by investigating the range of conformations sampled during replica-exchange MD simulations of ATP in explicit water. Modified parameters allowed a much wider range of conformational sampling compared with the bias toward extended forms with original parameters. A diverse range of structures agrees with the broad distribution of ATP conformations in proteins deposited in the Protein Data Bank. These simulations suggest that the modified parameters will be useful in studies of ATP in solution and of the many ATP-utilizing proteins.

Mechanism of the αβ conformational change in F1-ATPase after ATP hydrolysis: free-energy simulations

One of the motive forces for F1-ATPase rotation is the conformational change of the catalytically active β subunit due to closing and opening motions caused by ATP binding and hydrolysis, respectively. The closing motion is accomplished in two steps: the hydrogen-bond network around ATP changes and then the entire structure changes via B-helix sliding, as shown in our previous study. Here, we investigated the opening motion induced by ATP hydrolysis using all-atom free-energy simulations, combining the nudged elastic band method and umbrella sampling molecular-dynamics simulations. Because hydrolysis requires residues in the α subunit, the simulations were performed with the αβ dimer. The results indicate that the large-scale opening motion is also achieved by the B-helix sliding (in the reverse direction). However, the sliding mechanism is different from that of ATP binding because sliding is triggered by separation of the hydrolysis products ADP and Pi. We also addressed several important issues: 1), the timing of the product Pi release; 2), the unresolved half-closed β structure; and 3), the ADP release mechanism. These issues are fundamental for motor function; thus, the rotational mechanism of the entire F1-ATPase is also elucidated through this αβ study. During the conformational change, conserved residues among the ATPase proteins play important roles, suggesting that the obtained mechanism may be shared with other ATPase proteins. When combined with our previous studies, these results provide a comprehensive view of the β-subunit conformational change that drives the ATPase.

Evidence for an electrostatic mechanism of force generation by the bacteriophage T4 DNA packaging motor

How viral packaging motors generate enormous forces to translocate DNA into viral capsids remains unknown. Recent structural studies of the bacteriophage T4 packaging motor have led to a proposed mechanism wherein the gp17 motor protein translocates DNA by transitioning between extended and compact states, orchestrated by electrostatic interactions between complimentarily charged residues across the interface between the N- and C-terminal subdomains. Here we show that site-directed alterations in these residues cause force dependent impairments of motor function including lower translocation velocity, lower stall force and higher frequency of pauses and slips. We further show that the measured impairments correlate with computed changes in free-energy differences between the two states. These findings support the proposed structural mechanism and further suggest an energy landscape model of motor activity that couples the free-energy profile of motor conformational states with that of the ATP hydrolysis cycle.

Rotation triggers nucleotide-independent conformational transition of the empty β subunit of F₁-ATPase

F1-ATPase (F1) is the catalytic portion of ATP synthase, a rotary motor protein that couples proton gradients to ATP synthesis. Driven by a proton flux, the F1 asymmetric γ subunit undergoes a stepwise rotation inside the α3β3 headpiece and causes the β subunits' binding sites to cycle between states of different affinity for nucleotides. These concerted transitions drive the synthesis of ATP from ADP and phosphate. Here, we study the coupling between the mechanical progression of γ and the conformations of α3β3. Using molecular dynamics simulations, we show that the nucleotide-free β subunit, initially in the open, low-affinity state, undergoes a spontaneous closing transition to the half-open state in response to the γ rotation in the synthesis direction. We estimate the kinetics of this spontaneous conformational change and analyze its mechanism and driving forces. By computing free energy profiles, we find that the isolated empty β subunit preferentially adopts the half-open conformation and that the transition to this conformation from the fully open state is accompanied by well-defined changes in the structure and interactions of the active site region. These results suggest that ADP binding to F1 occurs via conformational selection and is preceded by the transition of the active site to the half-open conformation, driven by the intrinsic elasticity of β. Our results also indicate that opening of the nucleotide-free β during hydrolysis is not spontaneous, as previously assumed. Rather, the fully open conformation observed in the F1 X-ray structure is enforced sterically by the γ subunit whose orientation is stabilized by interactions with the two other β subunits in the completely closed state. This finding supports the notion that γ acts by coupling the extreme conformational states of β subunits within the α3β3 hexamer and therefore is responsible for high efficiency of the coordinated catalysis.

Phosphate release coupled to rotary motion of F1-ATPase

F1-ATPase, the catalytic domain of ATP synthase, synthesizes most of the ATP in living organisms. Running in reverse powered by ATP hydrolysis, this hexameric ring-shaped molecular motor formed by three αβ-dimers creates torque on its central γ-subunit. This reverse operation enables detailed explorations of the mechanochemical coupling mechanisms in experiment and simulation. Here, we use molecular dynamics simulations to construct a first atomistic conformation of the intermediate state following the 40° substep of rotary motion, and to study the timing and molecular mechanism of inorganic phosphate (Pi) release coupled to the rotation. In response to torque-driven rotation of the γ-subunit in the hydrolysis direction, the nucleotide-free αβE interface forming the "empty" E site loosens and singly charged Pi readily escapes to the P loop. By contrast, the interface stays closed with doubly charged Pi. The γ-rotation tightens the ATP-bound αβTP interface, as required for hydrolysis. The calculated rate for the outward release of doubly charged Pi from the αβE interface 120° after ATP hydrolysis closely matches the ~1-ms functional timescale. Conversely, Pi release from the ADP-bound αβDP interface postulated in earlier models would occur through a kinetically infeasible inward-directed pathway. Our simulations help reconcile conflicting interpretations of single-molecule experiments and crystallographic studies by clarifying the timing of Pi exit, its pathway and kinetics, associated changes in Pi protonation, and changes of the F1-ATPase structure in the 40° substep. Important elements of the molecular mechanism of Pi release emerging from our simulations appear to be conserved in myosin despite the different functional motions.

Molecular dynamics simulations of yeast F1-ATPase before and after 16° rotation of the γ subunit

We have recently proposed the "packing exchange mechanism" for F1-ATPase, wherein the perturbation by a substrate binding/release or an ATP hydrolysis is followed by the reorganization of the asymmetric packing structure of the α3β3 complex, accompanying the γ subunit rotation. As part of a further investigation of this rotational mechanism, we performed all-atom molecular dynamics simulations for yeast F1-ATPase both before and after a 16° rotation of the γ subunit triggered by a Pi release. We analyzed the structural fluctuations, the subunit interface interactions, and the dynamics of the relative subunit arrangements before and after the rotation. We found that with the Pi release the αEβE subunit interface becomes looser, which also allosterically makes the αDPβDP subunit interface looser. This structural communication between these interfaces takes place through a tightening of the αTPβTP subunit interface. The γ subunit interacts less strongly with αDP and βDP and more strongly with αTP and βTP. After the Pi release, the tightly packed interfaces are reorganized from the interfaces around βDP to those around βTP, inducing the 16° rotation. These results, which are consistent with the packing exchange mechanism, allow us to deduce a view of the structural change during the 40° rotation.