The molecular origin of the action of the F 0 proton gradient-driven rotor presents a major puzzle despite significant structural advances. Although important conceptual models have provided guidelines of how such systems should work, it has been challenging to generate a structure-based molecular model using physical principles that will consistently lead to the unidirectional protondriven rotational motion during ATP synthesis. This work uses a coarse-grained (CG) model to simulate the energetics of the F 0 -ATPase system in the combined space defined by the rotational coordinate and the proton transport (PTR) from the periplasmic side (P) to the cytoplasmic side (N). The model establishes the molecular origin of the rotation, showing that this effect is due to asymmetry in the energetics of the proton path rather than only the asymmetry of the interaction of the Asp on the c-ring helices and Arg on the subunit-a. The simulation provides a clear conceptual background for further exploration of the electrostatic basis of proton-driven mechanochemical systems.T he F 0 F 1 -ATPase is a ubiquitous nanomotor in all living cells that generate the ATP molecules essential for maintaining large gamut of cellular functions (1). This system is comprised of two rotary motors; the mechanochemical F 1 -ATPase that synthesizes or hydrolyzes ATP, utilizing the mechanical torque generated from the rotation of the stalk, and the membranebound F 0 -ATPase that drives the mechanical rotation utilizing the ion-motive force established across the membrane. The detailed nature of the conversion and the utility of energy by the F 0 F 1 -ATPase has long been one of the fundamental questions in biology. Significant progress has been achieved in understanding the mechanochemical coupling in the F 1 -ATPase, utilizing information from several high-resolution crystal structures and a wealth of biochemical and single-molecule spectroscopic data (2-6). Recently, the molecular nature of the coupling between chemical and conformational coordinates in the F 1 -ATPase has been elucidated using coarse-grained (CG) theoretical modeling approaches (7). In contrast, the molecular basis of the conversion of the pH gradient across the membrane producing directional rotation of the F 0 -ATPase is less understood, despite significant conceptual progress (3). The problems are partially due to the limited structural information about the complete membranebound F 0 -ATPase complex and also reflects the inherent difficulty in modeling coupled long-range biological processes, such as proton transfer and mechanical rotation.The F 0 motor is consisted of a rotor part, known as the c-ring, connected with the stator subunit-a and dimer subunits-b. The c-ring is a tightly packed ring-like structure composed of several α-helical hairpins whose number varies in different species (3,8). Most of the central part of the c-ring is embedded in the membrane, except for the cytoplasmic loops and the periplasmic termini. Each c-ring helices consists of a highly conserved ...