F-type ATP synthase (F-ATPase) and vacuolar ATP hydrolase (V-ATPase) are well-known biomolecular motors, which play significant catalytic roles in ATP synthesis and ATP hydrolysis reactions. Their rotational torques are important factors involved in their rotational behavior that can be measured experimentally but with considerable difficulty. To overcome this difficulty and thereby provide an in-depth understanding of their operation mechanism, we herein carry out simple and fast computer modelling to study the two proteins, using our torque approach that relies on interatomic forces and coordinates of unequilibrated configurations taken from brief molecular dynamics (MD) simulations. As predicted by the torque approach, F-ATPase is demonstrated to be a random rotor, but it prefers to rotate in clockwise direction (as seen from the membrane toward the protein) for ATP synthesis, owing to the predominantly negative angle-averaged torques. By contrast, V-ATPase tends to rotate only in counterclockwise direction for ATP hydrolysis, due to the almost uniform averaged positive torques generated by the unidirectional rotation near the three catalytic sites. The rotational behaviors of both proteins are also affected by the surrounding solvent which can promote or hinder the internal rotation. By combining the torque approach with classic force-field MD simulations, the torques of two biomolecular motors can be calculated economically, and are found to agree with previous experiments and theoretical calculations. This work demonstrates that our torque approach can be extended to the field of biology and can help gain a deeper insight into the mechanistic rotation of biomolecular motors with modest computation time.Communicated by Ramaswamy H. Sarma.
Keywords: F1-ATPase; MD simulation; V1-ATPase; biomolecular motor; torque approach.