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

Mechanism of proton-powered c-ring rotation in a mitochondrial ATP synthase

Proton-powered c-ring rotation in mitochondrial ATP synthase is crucial to convert the transmembrane protonmotive force into torque to drive the synthesis of adenosine triphosphate (ATP). Capitalizing on recent cryo-EM structures, we aim at a structural and energetic understanding of how functional directional rotation is achieved. We performed multi-microsecond atomistic simulations to determine the free energy profiles along the c-ring rotation angle before and after the arrival of a new proton. Our results reveal that rotation proceeds by dynamic sliding of the ring over the a-subunit surface, during which interactions with conserved polar residues stabilize distinct intermediates. Ordered water chains line up for a Grotthuss-type proton transfer in one of these intermediates. After proton transfer, a high barrier prevents backward rotation and an overall drop in free energy favors forward rotation, ensuring the directionality of c-ring rotation required for the thermodynamically disfavored ATP synthesis. The essential arginine of the a-subunit stabilizes the rotated configuration through a salt bridge with the c-ring. Overall, we describe a complete mechanism for the rotation step of the ATP synthase rotor, thereby illuminating a process critical to all life at atomic resolution.

Revealing a Hidden Intermediate of Rotatory Catalysis with X-ray Crystallography and Molecular Simulations

The mechanism of rotatory catalysis in ATP-hydrolyzing molecular motors remains an unresolved puzzle in biological energy transfer. Notwithstanding the wealth of available biochemical and structural information inferred from years of experiments, knowledge on how the coupling between the chemical and mechanical steps within motors enforces directional rotatory movements remains fragmentary. Even more contentious is to pinpoint the rate-limiting step of a multistep rotation process. Here, using vacuolar or V₁-type hexameric ATPase as an exemplary rotational motor, we present a model of the complete 4-step conformational cycle involved in rotatory catalysis. First, using X-ray crystallography, a new intermediate or "dwell" is identified, which enables the release of an inorganic phosphate (or P_i) after ATP hydrolysis. Using molecular dynamics simulations, this new dwell is placed in a sequence with three other crystal structures to derive a putative cyclic rotation path. Free-energy simulations are employed to estimate the rate of the hexameric protein transformations and delineate allosteric effects that allow new reactant ATP entry only after hydrolysis product exit. An analysis of transfer entropy brings to light how the side-chain-level interactions transcend into larger-scale reorganizations, highlighting the role of the ubiquitous arginine-finger residues in coupling chemical and mechanical information. An inspection of all known rates encompassing the 4-step rotation mechanism implicates the overcoming of the ADP interactions with V₁-ATPase to be the rate-limiting step of motor action.

Simulating a chemically fueled molecular motor with nonequilibrium molecular dynamics

Most computer simulations of molecular dynamics take place under equilibrium conditions-in a closed, isolated system, or perhaps one held at constant temperature or pressure. Sometimes, extra tensions, shears, or temperature gradients are introduced to those simulations to probe one type of nonequilibrium response to external forces. Catalysts and molecular motors, however, function based on the nonequilibrium dynamics induced by a chemical reaction's thermodynamic driving force. In this scenario, simulations require chemostats capable of preserving the chemical concentrations of the nonequilibrium steady state. We develop such a dynamic scheme and use it to observe cycles of a particle-based classical model of a catenane-like molecular motor. Molecular motors are frequently modeled with detailed-balance-breaking Markov models, and we explicitly construct such a picture by coarse graining the microscopic dynamics of our simulations in order to extract rates. This work identifies inter-particle interactions that tune those rates to create a functional motor, thereby yielding a computational playground to investigate the interplay between directional bias, current generation, and coupling strength in molecular information ratchets.

Bending-torsional elasticity and energetics of the plus-end microtubule tip

SignificanceThe mechanochemical basis of microtubule growth, which is essential for the normal function and division of eukaryotic cells, has remained elusive and controversial, despite extensive work. In particular, recent findings have created the paradox that the microtubule plus-end tips look very similar during both growing and shrinking phases, thereby challenging the traditional textbook picture. Our large-scale atomistic simulations resolve this paradox and explain microtubule growth and shrinkage dynamics as a process governed by energy barriers between protofilament conformations, the heights of which are in turn fine-tuned by different nucleotide states, thus implementing an information-driven Brownian ratchet.

SUMOylation of spastin promotes the internalization of GluA1 and regulates dendritic spine morphology by targeting microtubule dynamics

Dendritic spines are specialized structures involved in neuronal processes on which excitatory synaptic contact occurs. The microtubule cytoskeleton is vital for maintaining spine morphology and mature synapses. Spastin is related to microtubule-severing proteases and is involved in synaptic bouton formation. However, it is not yet known if spastin can be modified by Small Ubiquitin-like Modifier (SUMO) or how this modification regulates dendritic spines. Spastin was shown to be SUMOylated at K427, and its deSUMOylation promoted microtubule stability. In addition, SUMOylation of spastin was shown to affect signalling pathways associated with long term synaptic depression. SUMOylated spastin promoted the development of dendrites and dendritic spines. Moreover, SUMOylated spastin regulated endocytosis and affected the transport of the AMPA receptor, GluA1. Our findings suggest that SUMOylation of spastin promotes GluA1 internalization and regulates dendritic spine morphology through targeting of microtubule dynamics.

Theory of Nonequilibrium Free Energy Transduction by Molecular Machines

Biomolecular machines are protein complexes that convert between different forms of free energy. They are utilized in nature to accomplish many cellular tasks. As isothermal nonequilibrium stochastic objects at low Reynolds number, they face a distinct set of challenges compared with more familiar human-engineered macroscopic machines. Here we review central questions in their performance as free energy transducers, outline theoretical and modeling approaches to understand these questions, identify both physical limits on their operational characteristics and design principles for improving performance, and discuss emerging areas of research.

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.

Evidence for a Partially Stalled γ Rotor in F1-ATPase from Hydrogen-Deuterium Exchange Experiments and Molecular Dynamics Simulations

F₁-ATPase uses ATP hydrolysis to drive rotation of the γ subunit. The γ C-terminal helix constitutes the rotor tip that is seated in an apical bearing formed by α₃β₃. It remains uncertain to what extent the γ conformation during rotation differs from that seen in rigid crystal structures. Existing models assume that the entire γ subunit participates in every rotation. Here we interrogated E. coli F₁-ATPase by hydrogen-deuterium exchange (HDX) mass spectrometry. Rotation of γ caused greatly enhanced deuteration in the γ C-terminal helix. The HDX kinetics implied that most F₁ complexes operate with an intact rotor at any given time, but that the rotor tip is prone to occasional unfolding. A molecular dynamics (MD) strategy was developed to model the off-axis forces acting on γ. MD runs showed stalling of the rotor tip and unfolding of the γ C-terminal helix. MD-predicted H-bond opening events coincided with experimental HDX patterns. Our data suggest that in vitro operation of F₁-ATPase is associated with significant rotational resistance in the apical bearing. These conditions cause the γ C-terminal helix to get "stuck" (and unfold) sporadically while the remainder of γ continues to rotate. This scenario contrasts the traditional "greasy bearing" model that envisions smooth rotation of the γ C-terminal helix. The fragility of the apical rotor tip in F₁-ATPase is attributed to the absence of a c₁₀ ring that stabilizes the rotation axis in intact F_oF₁. Overall, the MD/HDX strategy introduced here appears well suited for interrogating the inner workings of molecular motors.

The regulatory subunit ε in Escherichia coli FOF1-ATP synthase

F-type ATP synthases are extraordinary multisubunit proteins that operate as nanomotors. The Escherichia coli (E. coli) enzyme uses the proton motive force (pmf) across the bacterial plasma membrane to drive rotation of the central rotor subunits within a stator subunit complex. Through this mechanical rotation, the rotor coordinates three nucleotide binding sites that sequentially catalyze the synthesis of ATP. Moreover, the enzyme can hydrolyze ATP to turn the rotor in the opposite direction and generate pmf. The direction of net catalysis, i.e. synthesis or hydrolysis of ATP, depends on the cell's bioenergetic conditions. Different control mechanisms have been found for ATP synthases in mitochondria, chloroplasts and bacteria. This review discusses the auto-inhibitory behavior of subunit ε found in F_OF₁-ATP synthases of many bacteria. We focus on E. coli F_OF₁-ATP synthase, with insights into the regulatory mechanism of subunit ε arising from structural and biochemical studies complemented by single-molecule microscopy experiments.

Free-energy simulations reveal molecular mechanism for functional switch of a DNA helicase

Helicases play key roles in genome maintenance, yet it remains elusive how these enzymes change conformations and how transitions between different conformational states regulate nucleic acid reshaping. Here, we developed a computational technique combining structural bioinformatics approaches and atomic-level free-energy simulations to characterize how the Escherichia coli DNA repair enzyme UvrD changes its conformation at the fork junction to switch its function from unwinding to rezipping DNA. The lowest free-energy path shows that UvrD opens the interface between two domains, allowing the bound ssDNA to escape. The simulation results predict a key metastable 'tilted' state during ssDNA strand switching. By simulating FRET distributions with fluorophores attached to UvrD, we show that the new state is supported quantitatively by single-molecule measurements. The present study deciphers key elements for the 'hyper-helicase' behavior of a mutant and provides an effective framework to characterize directly structure-function relationships in molecular machines.

Two zinc-binding domains in the transporter AdcA from facilitate high-affinity binding and fast transport of zinc

Zinc is an essential metal in bacteria. One important bacterial zinc transporter is AdcA, and most bacteria possess AdcA homologs that are single-domain small proteins due to better efficiency of protein biogenesis. However, a double-domain AdcA with two zinc-binding sites is significantly overrepresented in Streptococcus species, many of which are major human pathogens. Using molecular simulation and experimental validations of AdcA from Streptococcus pyogenes, we found here that the two AdcA domains sequentially stabilize the structure upon zinc binding, indicating an organization required for both increased zinc affinity and transfer speed. This structural organization appears to endow Streptococcus species with distinct advantages in zinc-depleted environments, which would not be achieved by each single AdcA domain alone. This enhanced zinc transport mechanism sheds light on the significance of the evolution of the AdcA domain fusion, provides new insights into double-domain transporter proteins with two binding sites for the same ion, and indicates a potential target of antimicrobial drugs against pathogenic Streptococcus species.

Evolution and molecular mechanism of PitAs in iron transport of Streptococcus species

Iron is an essential element for almost all bacteria. The iron ATP-binding cassette (ABC) transporters located on the cell membrane affects bacterial virulence and infection. Although a variety of Fe³⁺-transporters have been found in bacteria, their evolutionary processes are rarely studied. Pneumococcal iron ABC transporter (PitA), a highly conserved Fe³⁺-transporter in most pathogenic bacteria, influences the capsule formation and virulence of bacteria. However, multiple sequence alignment revealed that PitA is expressed in four different variants in bacteria, and the structural complexity of these variants increases progressively. To more efficiently import Fe³⁺ ions into bacterial cells, bacteria have evolved a fused PitA from two separately expressed PitA-1 (SPD_0227) and PitA-2 (SPD_0226) proteins. Further biochemical characterization indicated that both PitA-1 and PitA-2 have weaker Fe³⁺-binding ability than their protein complex. More importantly, Glutathione S-Transferase (GST) pull-down and isothermal titration calorimetry (ITC) detection showed that PitA-1 and PitA-2 interact with each other via Tyr111-Leu37, Asn112-Gln38, Asn103-Leu33, and Asn103-Thr34. Further molecular dynamics (MD) simulations demonstrated that this interaction in full-length PitA is stronger than that in the two individual proteins. Deletion of PitA family genes could lead to decrease in the ability of iron acquisition and of adhesion and invasion of S. pneumoniae. Our study revealed the evolving state and molecular mechanism of Fe³⁺-transporter PitAs in bacteria and provided important information for understanding the iron transportation mechanism in bacteria and designing new antibacterial drugs.