F 1 -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°s ubstep of rotary motion, and to study the timing and molecular mechanism of inorganic phosphate (P i ) 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 P i readily escapes to the P loop. By contrast, the interface stays closed with doubly charged P i . The γ-rotation tightens the ATP-bound αβ TP interface, as required for hydrolysis. The calculated rate for the outward release of doubly charged P i from the αβ E interface 120°after ATP hydrolysis closely matches the ∼1-ms functional timescale. Conversely, P i 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 P i exit, its pathway and kinetics, associated changes in P i protonation, and changes of the F 1 -ATPase structure in the 40°substep. Important elements of the molecular mechanism of P i release emerging from our simulations appear to be conserved in myosin despite the different functional motions. F 1 -ATPase (F 1 ), the catalytic domain of F o F 1 -ATP synthase, is a rotary molecular motor that reversibly interconverts ATP hydrolysis free energy and mechanical work associated with the rotation of the central stalk (1). The minimal functional F 1 consists of a hexameric ring formed by three αβ-subunit dimers, with the rod-like γ-subunit located at its center (2). The rotation of the γ-subunit is tightly coupled to the reactions in the three catalytic sites located at the αβ interfaces and hosted mainly by the β-subunits. As a result, F 1 is a unique reversible motor that rotates γ by converting ATP hydrolysis energy at high efficiency (3-5) and, conversely, synthesizes ATP from ADP and inorganic phosphate (P i ) by forced rotation of γ in the reverse direction (6, 7). The three nucleotide-binding sites are in different phases of catalysis, reflecting the asymmetric structure of the γ-subunit. Correspondingly, the αβ-subunits hosting the catalytic interfaces are in different conformational states, empty (E), ATP-bound (TP), and ADP+P i -bound (DP), as seen in crystal structures (2). They communicate through the γ-subunit or directly within the α 3 β 3 ring (8, 9) and cooperatively drive the rotation of the γ-subunit.Single-molecule experiments have shown that the ...