Motor protein with nonequilibrium potential: Its thermodynamics and efficiency

Physical Review E - Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics

Volumn 69, Issue 1, 2004, Pages 4-

  Qian, Hong ^a  

^a UNIVERSITY OF WASHINGTON   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

ADENOSINETRIPHOSPHATE; BOUNDARY CONDITIONS; BROWNIAN MOVEMENT; CELLS; COMPUTATION THEORY; CONFORMATIONS; HYDROLYSIS; PROTEINS; QUANTUM THEORY;

LANGEVIN EQUATIONS; NONEQUILIBRIUM STEADY STATE (NESS);

MOLECULAR PHYSICS;

ADENOSINE DIPHOSPHATE; ADENOSINE TRIPHOSPHATE; MOLECULAR MOTOR; PROTEIN;

BIOLOGICAL MODEL; CHEMICAL MODEL; CHEMISTRY; COMPUTER SIMULATION; ENERGY TRANSFER; METABOLISM; MOTION; PHYSIOLOGY; THERMODYNAMICS;

ADENOSINE DIPHOSPHATE; ADENOSINE TRIPHOSPHATE; COMPUTER SIMULATION; ENERGY TRANSFER; MODELS, BIOLOGICAL; MODELS, CHEMICAL; MOLECULAR MOTOR PROTEINS; MOTION; PROTEINS; THERMODYNAMICS;

EID: 85036373693     PISSN: 1063651X     EISSN: None     Source Type: Journal    
DOI: 10.1103/PhysRevE.69.012901     Document Type: Article

References (41)

1. P.G. Bergmann and J.L. Lebowitz, Phys. Rev. 99, 578 (1955). The paper sets up the foundation for using stochastic dynamics to study nonequilibrium processes in both closed and open systems. Its Sec. IV provided a mathematical justification for the existence of NESS in open systems.
2. Section III of Ref. 1 also showed that in a closed system Helmholtz potential decreases monotonically in the course of time and is at its minimum when equilibrium canonical distribution is reached. This result has been rediscovered independently in the mesoscopic theory of single macromolecules: H. Qian, Phys. Rev. E 63, 042103 (2001);
3. Phys. Rev. EH. Qian65, 016102 (2002).
4. H. Qian, Proc. R. Soc. London, Ser. A 457, 1645 (2001);
5. H. QianJ. Phys. Chem. 106, 2065 (2002);
6. H. Qian, M. Qian, and X. Tang, J. Stat. Phys. 107, 1129 (2002);
7. K.H. Kim and H. Qian, e-print physics/0303016.
8. R. Graham and H. Haken, Phys. A. 243, 289 (1971);
9. R. Graham, Springer Tracts Mod. Phys. 66, 1 (1973);G. Nicolis and I. Prigogine, Self-organization in Nonequilibrium Systems (Wiley, New York, 1977).
10. F. Jülicher, A. Aidari, and J. Prost, Rev. Mod. Phys. 69, 1269 (1997);
11. C. Bustamante, D. Keller, and G. Oster, Acc. Chem. Res. 34, 412 (2001);
12. P. Reimann, Phys. Rep. 361, 57 (2002);
13. M.R. Parrondo and B.J. de Cisneros, Appl. Phys. A: Mater. Sci. Process. 75, 2, 179 (2002).
14. M.E. Fisher and A.B. Kolomeisky, Proc. Natl. Acad. Sci. U.S.A. 96, 6597 (1999);
15. M.E. FisherA.B. KolomeiskyPhysica A 274, 241 (1999);
16. H. Qian, Biophys. Chem. 67, 263 (1997);
17. Biophys. Chem.H. Qian83, 35 (2000).
18. H. Wang and G. Oster, Nature (London) 396, 279 (1998).
19. M.O. Magnasco, Phys. Rev. Lett. 72, 2656 (1994);
20. D. Keller and C. Bustamante, Biophys. J. 78, 541 (2000). The important difference between this class of models and that of rectified Brownian motion 9 and fluctuating potential models 10 is that the driving force is introduced into the x direction a priori. Hence, the energy transduction is built into the model explicitly rather than being an emerging property.
21. R.F. Fox, Phys. Rev. E 57, 2177 (1998);
22. Phys. Rev. ER.F. Fox and M.H. Choi, 63, 051901 (2001).
23. H. Qian, Phys. Rev. Lett. 81, 3063 (1998);
24. H. QianJ. Math. Chem. 27, 219 (2000);
25. H. Qiane-print cond-mat/0106302.
26. J. Villá and A. Warshel, J. Phys. Chem. 105, 7887 (2001);
27. A. Warshel, Acc. Chem. Res. 35, 385 (2002);
28. W. Yang, Proc. Natl. Acad. Sci. U.S.A. 100, 874 (2003).
29. If in the exchange process the lifetime of the motor protein without nucleotide is significantly long, and the dynamics of the unligated motor protein is complex, a second molecular dynamic simulation should be introduced for the bare protein. This further complicates the model but introduces no conceptual difficulties. Hence we assume that the exchange process is fast and simple in the present work.
30. (Formula presented) ’s are the macroscopic rate constants for the ADP-ATP exchange and (Formula presented) characterizes the location of the transition state of the reaction. See Eq. (6) for more discussions.
31. C.-X. Wu and M.-Z. Guo, Acta Math. (China) 1, 195 (1981).
32. H.-C. Kaiser, H. Neidhardt, and J. Rehberg, J. Math. Phys. 43, 5325 (2002).
33. With the presence of (Formula presented) the boundary conditions are (Formula presented) and (Formula presented)
34. A. Parmeggiani, F. Jülicher, A. Ajdari, and J. Prost, Phys. Rev. E 60, 2127 (1999).
35. I. Derényi, M. Bier, and R.D. Astumian, Phys. Rev. Lett. 83, 903 (1999).
36. H. Wang and G. Oster, Europhys. Lett. 57, 134 (2002).
37. The frictional coefficient of motor protein is (Formula presented) So when a frictional object with (Formula presented) is attached, assuming the force generated by the motor is F, then the velocity of the complex is (Formula presented) where f is the tension between the motor and the object. Hence, the total heat dissipation rate (Formula presented) and the (Formula presented) associated with dragging the object is (Formula presented) Hence (Formula presented)
38. The phase space of our model is called cylinder in dynamical systems theory. Stochasticity aside, there are two classes of periodic solutions to a dynamical system on a cylinder: zero and nonzero rotation numbers. The former is a futile cycle in our theory. The tight coupling is related to periodic solution with rotation number 1. Higher rotation numbers correspond to multiple steps with single hydrolysis.
39. B. Derrida, Phys. Rep. 301, 65 (1998);
40. J.L. Lebowitz and H. Spohn, J. Stat. Phys. 95, 333 (1999).
41. Computing the potential of mean force (Formula presented) which involves chemical bond breaking requires a hybrid quantum mechanics/molecular mechanics computation 11. After that, the standard state chemical potentials of the discrete states can be computed as follows: (Formula presented)and similarly for (Formula presented)