Multiscale Modelling for investigating single molecule effects on the mechanics of actin filaments

Abstract

This work presents a preliminary multiscale computational investigation of the effects of nucleotides and cations on the mechanics of actin filaments (F-actin). At the molecular level, Molecular Dynamics (MD) simulations are employed to characterize the rearrangements of the actin monomers (G-actin) in terms of secondary structures evolution in physiological conditions. At the mesoscale level, a coarse grain (CG) procedure is adopted where each monomer is represented by means of Elastic Network Modeling (ENM) technique. At the macroscale level, actin filaments up to hundreds of nanometers are assumed as isotropic and elastic beams and characterized via Rotation Translation Block (RTB) analysis. F-actin bound to adenosine triphosphate (ATP) shows a persistence length around 5 μm, while actin filaments bound to adenosine diphosphate (ADP) have a persistence length of about 3 μm. With magnesium bound to the high affinity binding site of G-actin, the persistence length of F-actin decreases to about 2 μm only in the ADP-bound form of the filament, while the same ion has no effects, in terms of stiffness variation, on the ATP-bound form of F-actin. The molecular mechanisms behind these changes in flexibility are herein elucidated. Thus, this study allows to analyze how the local binding of cations and nucleotides on G-actin induce molecular rearrangements that transmit to the overall F-actin, characterizing shifts of mechanical properties, that can be related with physiological and pathological cellular phenomena, as cell migration and spreading. Further, this study provides the basis for upcoming investigating of network and cellular remodelling at higher length scales.