Energy-based modelling of single actin filament polymerization using bond graphs

Bond graphs provide an energy-based methodology for modelling complex systems hierarchically; at the moment, the method allows biological systems with both chemical and electrical subsystems to be modelled. Herein, the bond graph approach is extended to include chemomechanical transduction thus extending the range of biological systems to be modelled. Actin filament polymerization and force generation is used as an example of chemomechanical transduction, and it is shown that the TF (transformer) bond graph component provides a practical, and conceptually simple, alternative to the Brownian ratchet approach of Peskin, Odell, Oster and Mogilner. Furthermore, it is shown that the bond graph approach leads to the same equation as the Brownian ratchet approach in the simplest case. The approach is illustrated by showing that flexibility and non-normal incidence can be modelled by simply adding additional bond graph components and that compliance leads to non-convexity of the force-velocity curve. Energy flows are fundamental to life; for this reason, the energy-based approach is utilized to investigate the power transmission by the actin filament and its corresponding efficiency. The bond graph model is fitted to experimental data by adjusting the model physical parameters.