Mechanochemical simulations of actomyosin networks are traditionally based on one-dimensional models of actin filaments having zero width. Here, and in the follow up paper (arXiv, DOI 10.48550/arXiv.2203.01284), approaches are presented for more efficient modeling that incorporates stretching, shearing, and twisting of actin filaments. Our modeling of a semiflexible filament with a small but finite width is based on the Cosserat theory of elastic rods, which allows for six degrees of freedom at every point on the filament’s backbone. In the variational models presented in this paper, a small and discrete set of parameters is used to describe a smooth filament shape having all degrees of freedom allowed in the Cosserat theory. Two main approaches are introduced: one where polynomial spline functions describe the filament’s configuration, and one in which geodesic curves in the space of the configurational degrees of freedom are used. We find that in the latter representation the strain energy function can be calculated without resorting to a small-angle expansion, so it can describe arbitrarily large filament deformations without systematic error. These approaches are validated by a dynamical model of a Cosserat filament, which can be further extended by using multiresolution methods to allow more detailed monomer-based resolution in certain parts of the actin filament, as introduced in the follow up paper. The presented framework is illustrated by showing how torsional compliance in a finite-width filament can induce broken chiral symmetry in the structure of a cross-linked bundle.

We present a multi-resolution methodology for modelling F-actin filaments. It provides detailed microscopic information at the level of individual monomers at a lower computational cost by replacing the monomer-based model in parts of the simulated filament by a rod-based macroscopic model. In the monomer-based description, G-actin is represented by ellipsoids bound at the surface in a double helical configuration to form F-actin. The rod-based model is coarser, in which Factin is described using a Cosserat model, as seen in the preceding paper. 1 The multi-resolution methodology is illustrated using three case studies, designed to test the properties of F-actin under stretching, bending, shearing and twisting. The methodology is especially suited for situations where filaments are subject to bending deformations. We investigate the limitations of using the standard Cosserat model to capture the complete torsional behaviour of F-actin, presenting its extensions which account for curvature dependent rigidities and a twist-stretch coupling to improve accuracy of the overall multi-resolution scheme.

This work investigates multi-resolution methodologies for simulating dimer models. The solvent particles which make up the heat bath interact with the monomers of the dimer either through direct collisions (short-range) or through harmonic springs (long-range). Two types of multi-resolution methodologies are considered in detail: (a) describing parts of the solvent far away from the dimer by a coarser approach; (b) describing each monomer of the dimer by using a model with different level of resolution. These methodologies are then used to investigate the effect of a shared heat bath versus two uncoupled heat baths, one for each monomer. Furthermore, the validity of the multi-resolution methods is discussed by comparison to dynamics of macroscopic Langevin equations.

Mechanochemical simulations of actomyosin networks are traditionally based on onedimensional models of actin filaments having zero width. Here, and in the follow up paper, 1 approaches are presented for more efficient modelling which incorporates stretching, bending, shearing and twisting of actin filaments. Our modelling of a semiflexible filament with a small but finite width is based on the Cosserat theory of elastic rods, which allows for six degrees of freedom at every point on the filament's backbone. In the variational models presented in this paper, a small and discrete set of parameters is used to describe a smooth filament shape having all degrees of freedom allowed in the Cosserat theory. Two main approaches are introduced: one where polynomial spline functions describe the filament's configuration, and one in which geodesic curves in the space of the configurational degrees of freedom are used. We find that in the latter representation the strain energy function can be calculated without resorting to a small-angle expansion, so it can describe arbitrarily large filament deformations without systematic error. These approaches are validated by a dynamical model of a Cosserat filament, which can be further extended by using multi-resolution methods to allow more detailed monomer-based resolution in certain parts of the actin filament, as introduced in the follow up paper. 1 The presented framework is illustrated by showing how torsional compliance in a finite-width filament can induce broken chiral symmetry in the structure of a cross-linked bundle.

No abstract