Abstract
Through a variety of mechanisms, a healthy heart is able to regulate its structure and dynamics across multiple length scales. Disruption of these mechanisms can have a cascading effect, resulting in severe structural and/or functional changes that permeate across different length scales. Due to this hierarchical structure, there is interest in understanding how the components at the various scales coordinate and influence each other. However, much is unknown regarding how myofibril bundles are organized within a densely packed cell and the influence of the subcellular components on the architecture that is formed. To elucidate potential factors influencing cytoskeletal development, we proposed a computational model that integrated interactions at both the cellular and subcellular scale to predict the location of individual myofibril bundles that contributed to the formation of an energetically favorable cytoskeletal network. Our model was tested and validated using experimental metrics derived from analyzing single-cell cardiomyocytes. We demonstrated that our model-generated networks were capable of reproducing the variation observed in experimental cells at different length scales as a result of the stochasticity inherent in the different interactions between the various cellular components. Additionally, we showed that incorporating length-scale parameters resulted in physical constraints that directed cytoskeletal architecture toward a structurally consistent motif. Understanding the mechanisms guiding the formation and organization of the cytoskeleton in individual cardiomyocytes can aid tissue engineers toward developing functional cardiac tissue.