Evaluation of extensional and torsional stiffness of single actin filaments by molecular dynamics analysis

It is essential to investigate the mechanical behaviour of cytoskeletal actin filaments in order to understand their critical role as mechanical components in various cellular functional activities. These actin filaments consisting of monomeric molecules function in the thermal fluctuations. Hence, it is important to understand their mechanical behaviour on the microscopic scale by comparing the stiffness based on thermal fluctuations with the one experimentally measured on the macroscopic scale. In this study, we perform a large-scale molecular dynamics (MD) simulation for a half-turn structure of an actin filament. We analyse its longitudinal and twisting Brownian motions in equilibrium and evaluated its apparent extensional and torsional stiffness on the nanosecond scale. Upon increasing the sampling-window durations for analysis, the apparent stiffness gradually decreases and exhibits a trend to converge to a value that is close to the experimental value. This suggests that by extrapolating the data obtained in the MD analysis, we can estimate the experimentally determined stiffness on the microsecond to millisecond scales. For shorter temporal scales, the apparent stiffness is larger than experimental values, indicating that fast, local motions of the molecular structure are dominant. To quantify the local structural changes within the filament on the nanosecond scale and investigate the molecular mechanisms, such as the binding of the actin-regulatory proteins to the filaments, it is preferable to analyse the mechanical behaviour on the nanometre and nanosecond scales using MD simulation.