Assessing Adherent Cell Mechanical Behavior by 6 and 12-Strut Tensegrity Model

The cytoskeleton is the main cellular component responsible for bearing mechanical loads and supporting the cell shape in the adherent cells. Number of structures have been proposed to investigate the behavior of the cytoskeleton which most of them were based on stretch elements which only stand against tension forces. actin filaments are rope-shaped filaments that only tolerate tensile forces. However, in biological systems, microtubules have a vital role in the behavior of living cells. Microtubules, which are tube-shaped filaments, have evolved to withstand compressive forces. In this way, the need to have a structure that has tensile and compression elements together for the better approximation of cell behavior is necessary. Hence, tensegrity structures were considered as a tool for modelling and understanding the cytoskeleton behavior. Multiple experiments have proven that the tensegrity model is an appropriate structure in approximating the cytoskeleton behavior of non-suspension cells. In this study, we suggested two different tensegrity structures to examine the effect of the complexities of these structures in response to mechanical stimulus. Both 6 and 12-strut structures were simulated and subjected to theoretical atomic force microscopy loading through the software ABAQUS. Next, the relation between structure reaction forces and the amount of deflection, and also cell response to the stiffness change of the substrate were investigated. After comparing behavior of these two structures, it was concluded that it is more appropriate to apply a tensegrity model with fewer struts in stiffer cells’ behavior. The structure behavior in response to changes in substrate stiffness revealed that six-strut tensegrity is more sensitive to higher substrates stiffness.

actin filaments. Each of the cell components has its own specific role in cell behavior and balance [8]. In tensegrity structure, pretension in actin filaments can be balanced with focal adhesion, that is why tensegrity is not a proper structure to consider non-adherent cells.
Because in non-adherent cells there is no focal adhesion to transfer mechanical stress through the membrane [9]. In this study, only adherent cells have been considered. Previous studies applied the structures of different complexities for studying various aspects of the cells' mechanical behavior [1,[10][11][12]. Prendergast, et al. [13] studied tissue response to mechanical stresses, trying to model a single spread cell adhered to the substrate while influenced by shear stress. He also used 6-strut tensegrity structure for the purpose of modelling cytoskeleton and solid parts for the membrane, nucleus and cytoplasm [13]. Kardas, et al. [14] modelled an osteocyte inside the lacunae.
The cell structure has involved a different type of cell component like cytoskeleton, nucleus, and integrin. Three different models made by fiber elements have been created for presenting actin network, intermediate filaments connected to the nucleus, and microtubules which were connected to the centrosome, while a 30 struts tensegrity has been chosen for modeling the nucleus. They also examined various types of cytoskeleton fibers arrangement in the cell. Their study showed that random arrangement is more consistent with experimental observations [14]. The way of estimating the nucleus by using the tensegrity structure has been recorded by Ingber [15,16]. another research on tensegrity as cell cytoskeleton structure is the Bursa and Fuis study [17]. They presented a FEM model for the eukaryotic cell. Their model involved tensegrity (cytoskeleton), membrane, cytoplasm, and nucleus. They simulated 60 peripheral cables and 30 radial cables, which had the same joint points at the cellular level.
This model is proper to show transition extracellular force to the fiber part of the cell. De Santis, et al. [18] have used the tensegrity structure for investigating cell spreading behavior on the substrate.
They used a FEM model of 6-strut tensegrity to find out how the cell feels the surface stiffness. They applied a tensegrity structure with a solid membrane, cytoplasm, and nucleus to a spread cell attached to an elastic substrate. They finally observed that there is a one-to-one relationship between microtubules' compression and the substrate elasticity [18]. To compare the complexity of the structure, Chen, et al. [19], used 6 and 12-strut to simulate intracellular force distribution and stored energy in cell spreading.
They noticed that 12-strut tensegrity has more flexibility to spread.
They also showed that during cell spreading, the stored energy in cables increased but microtubules would limit the cell spreading by decreasing energy to zero [19]. Several usages of tensegrity have been found for estimating cellular behavior, but there is a simple question here. Could the same structure be used for different cell phenotypes with different stiffness? How should it be dealt with cell stiffness?
The objective of this study is to assess the role of tensegrity's structure complexity on adherent cell behavior cell. A vertical deflection was applied to the structure to mimic BioMEM system for measuring cell compliance [20] experiment, and the reaction forces of the structure to this displacement were assessed. Then, the reaction force of this structure vs. the substrate stiffness changes was investigated to study the effect of the substrate stiffness on cellular behavior. In this structure, there are three struts in each hypothetical plane. The struts are placed together in a triangle-like structure and the head or end of the strut are connected to the nearest adjacent vertices through four cables. The 12-strut tensegrity model is completed by adding 48 cables [21]. In the cytoskeleton, the end of each microtubule which is close to the membrane is connected to the integrin. Integrin is a transversal membrane's proteins that are contacted to the extracellular environments outside of the cell and transfer all extracellular interactions into the cell. Inside the cell, integrin sends the extracellular forces toward the cytoskeleton and makes the filaments to sense tensile stresses. This initial tensile force is subjected to microtubules compression and causes a forced balance in the whole cytoskeleton. This force balance of filamentary strands eventually leads to balance in the shape and structure of the cell. In fact, the stability of the tensegrity structure is dependent on the pretension of cables. In this simulation, pres-tress of 82kPa for the initial tension of cables and 12.4kPa for compressive prestress of struts were applied [22].

Methods
For substrate, an isotropic elastic solid disk was created. Then, in order to connect the structure to the substrate, the tensegrity structure was placed on the substrate in such a way that some of the struts ends touched the substrate. The cell model structure had was created using 6 and 12-strut tensegrity, the cell was places on the substrate so three and four cables were contiguous with the substrate, respectively (Figure 1).

Dimensions and Mechanical Properties of Cytoskeleton Components
In the 6-strut tensegrity model, the length of each strut was 14µm, and the distance between two parallel struts was 7µm with a strut cross-section area of 1340 nm 2 . Young's modulus and Poisson's ratio was defined as 1.2GPa and 0.3, respectively [23,24]. Cables

Loads and Boundary Condition
AFM is a regular and standard experimental method in cell mechanical behavior studies [25]. However, AFM has a sharp point of the contact that may result in the local strain concentration.
The alternative method of cell loading using a BioMEM system for measuring cell compliance was suggested [20] that allowed application of the load to the whole upper surface of the cell; thus a similar displacement was applied to all points representing the cell's upper surface. In order to simulate the BioMEMs load, a downward displacement was applied to the upper vertices of the structure (3 and 4 vertices for 6 and 12-strut tensegrity, respectively). All degrees of freedom for the nodes located at the bottom surface of the disk (substrate) were fixed. In order to attach the tensegrity structure to the substrate, all connected vertices with the disc (located at the end of microtubules) were pinned to the nearest node of the substrate meshes. To study the cytoskeleton behavior, the applied displacement to the upper vertices of the structure (Figure 2) was incremented from 0.1 to 1 µm, and the total reaction force at these nodes was calculated. Here, parameters ρ is the density, b is the body force and a is the linear acceleration. For static condition and avoiding the presence of body forces like weight, Eq. (1) will be modified as: To solve this differential equation with the FEM, ABAQUS (version 6.4) solver was used [26]. ABAQUS generates a system of equations and solves them by calculating the stiffness matrix and the displacements at any nodes.

Meshing
To model the cables and struts, tension only and compression only truss elements (T3D2) were chosen respectively. For meshing the substrate, the 8-node 3D structural element (C3D8R) was.
After the mesh convergence 1824 elements were obtained for the substrate. Nonlinear geometry (Nlgeom) option was checked in step modulus to encounter large deformation problem.

Effect of the Variation of Number of The Strut Elements
The color representation of the stress magnitudes carried by the struts and cables are shown in Figure 3. deformation shows that increased external displacement leads to increases in the reaction force, as one should expect ( Figure 4).
The obtained data shows that the 6-strut tensegrity model exhibits larger stiffness compared to a 12-strut. Therefore, in order to model behavior of the cell with higher stiffness, application of the 6-struts tensegrity model may be advantageous.

Effect of the Substrate Stiffness
Substrate stiffness is one of the important factors affecting the adherent cell in vitro behavior. Depending on their types, cells are sensitive to a specified range of substrate stiffness. The changes in substrate stiffness cause changes in cytoskeleton balancing forces and modifies cell movement along the substrate [18,27]. In order to evaluate the effect of substrate stiffness on the cytoskeleton, substrate stiffness was changed from 10 Pa to 100 MPa and the tensegrity structure's reaction force was calculated for 1µm downward displacement. Figure 5 demonstrated the variation of the reaction force due to the change of the substrate stiffness for δ=1 μm. Let us define the area where reaction force is sensitive to the change of the structure's stiffness as an area of cell sensitivity to the substrate stiffness. This area has a range of 1 to 100 kPa for 6-strut tensegrity. While 12-strut tensegrity does not depend on stiffness greater than 1 kPa.
Since the forces of the cytoskeleton are balanced through substrate stiffness, it can be claimed that the 6-strut tensegrity compared to 12-strut structure approximates stiffer cells. So, an increase in the number of struts makes the structure more flexible and more appropriate to estimate softer cells' modeling.
The variation of the reaction force along with the stiffness of the substrate is shown in Figure 6. As demonstrated here, the sensitivity area is independent of the amount of motility, while reaction force increases by a rise in displacement. The reaction force changes against substrate stiffness for the 12-strut structure are shown in Figure 7. In a 6-strut structure, the response force sensitivity increases with increasing displacement, while the sensitivity area is the same and does not change by increasing displacement.
By comparing Figures 5 & 6, it is concluded that in the 6-struts structures, the variation of the reaction force is greater than that of the 12-strut structure. This 6-strut tensegrity is more appropriate for modelling the stiffer cell behavior.

Discussion
Each spherical tensegrity model can be built with the help of a regular polyhedron. Each vertex of the structure is aligned with the polyhedron's corners located on the surface of the sphere. The obtained results showed that a tensegrity with fewer number of struts creates a stiffer structure. By creating tensegrity structures of the same size, but with increased number of struts, the structure configuration will look more similar to the dome while the length of the struts will decrease, and stiffness of the structure will increase.
By reducing the length of the struts, and increasing distance from the center of the sphere, struts have less contribution in tolerating foreign forces, and this makes the structure more flexible. In the area of the cell sensitivity ( Figure 5), the balance between cytoskeleton structure and mechanical environment depends on the substrate stiffness. Under the applied external displacement, the shape of the structure will change until the internal forces of the components can withstand the external stimulus. The cytoskeleton attaches to integrin within the membrane by the ends of the microtubules.
Hence, integrin could be responsible for microtubule mechanical connection to ECM in biological condition and to the substrate in vitro [28] In vitro conditions the microtubules' internal forces fluctuate due to the changes in substrate properties. Due to such changes the microtubule internal force will change and, accordingly, the cytoskeleton shape will change to reach another equilibrium configuration.

Conclusion
In this study, two types of tensegrity structure were studied as possible models to approximate the cytoskeleton behavior. Cells were loaded by a uniform displacement applied to the top surface.
It is recommended to use a structure with fewer struts to simulate stiffer cell behavior. The cell was attached to the substrate to show in vitro condition. Evaluation of the cell response to alerting substrate mechanical properties (Young's modulus) showed lower sensitivity range for the structure with more compressive elements.
This outcome is compatible with the previous studies and confirms that for simulating softer cells, a tensegrity model with more struts is more suitable.