A linear programming approach to reconstructing subcellular structures from confocal images for automated generation of representative 3D cellular models
Graphical abstract
Highlights
► We reconstruct representative cellular structural elements of single cells. ► Nucleus and F-actin network of vascular smooth muscle cells analyzed. ► Segmentation, fiber directionality and linear programming techniques utilized. ► Fully automated process creates geometries for importing into finite element models.
Introduction
While cell mechanics has been recognized as an important area of study, current computational models to interpret experimental results tend to ignore individual cellular geometries. In particular, 3D computational models could help to improve the design of experiments to characterize cell mechanical properties and interactions. This could lead to reduced times for discovery of mechanobiology principles and to faster translation of those principles from benchtop to bedside in clinically relevant devices and medications. The goal of this study is to create a fully automated algorithm capable of reconstructing the geometries of the cell membrane, nucleus, and actin stress fiber network of single cells in 3D. We seek to accomplish this by processing fluorescent confocal microscopy images of each of those cellular components in such a way that the resulting geometries are optimized for structural analysis using finite element methods. If generated, such geometries could be utilized in various types of multiscale models to bridge the gap between the nano- and macro-scale models currently in use.
The traditional primary focus of modern medical research is the investigation of molecular biology and genetic factors in disease, which sometimes leads to a tendency to ignore changes in tissue structure and mechanics that can also lead to pain and morbidity (Ingber, 2003a). However, that lack of focus on the physical basis of disease has been changing in recent years with the growing emphasis on evidence-based medicine in US hospitals (Fielding and Teutsch, 2011, Kaufman, 2010) together with the substantial growth and maturation of the field of mechanobiology over the past decade (Butler and Wang, 2011). Indeed, there has been a great deal of effort to develop geometrically accurate 3D structural models at both the tissue and molecular levels (Biswas et al., 2009, Wu et al., 2010). However, there has been much less effort focused at the single-cell level and therefore comparatively little progress has been made toward generation of equally accurate 3D representations of the structural components of single cells.
The ability to predict the behavior of cells from their sub-micron and nanoscale structures could elucidate the mechanisms behind many tissue mechanical properties (Ingber, 2003b). For as long as there have been observations of the mechanical properties of cells, there have been models put forth to attempt to describe those observations. At the most basic level, there are two categories of these models: continuum and structure-based. Continuum models, which lack internal structure, were the first type of model utilized to describe the mechanical behavior of cells and generally consider the cell to be equivalent to a simple “balloon full of molasses” (Ingber, 2003b, Li et al., 2007). These types of models therefore make predictions with minimal use of geometric variables (Cao and Chandra, 2010, Unnikrishnan et al., 2007). Despite the growing amount of evidence in support of the importance of structural elements within cells that has been published throughout the past several decades (Bathe et al., 2008, Bursac et al., 2005, Chaudhuri et al., 2007, Deng et al., 2006, Deshpande et al., 2008, Hardin and Walston, 2004, Hawkins et al., 2010, Hemmer et al., 2009, Ingber, 2003a, Ingber, 2003b, Ingber, 2003c, Kasza et al., 2007, Li, 2008, Mizuno et al., 2007, Pollard, 2003, Pullarkat et al., 2007, Stamenović, 2005, Stamenović, 2008, Stamenović et al., 2009, Suresh, 2007, Tseng et al., 2005), these types of models remained popular with bioengineers due to their relative simplicity and ease of implementation.
Structure-based models, on the other hand, are comprised of one or more networks of discrete structural elements that work in harmony to determine the mechanical responses of cells. These models tend to utilize Finite Element Analysis (FEA) to allow for analysis of complicated cellular and subcellular geometries. Many single-cell Finite Element Models (FEMs) rely on idealized geometries (Karcher et al., 2003, Peeters et al., 2005, Unnikrishnan et al., 2007), however recent efforts have incorporated geometries obtained from image segmentation. The first efforts to generate accurate 3D representations of subcellular structural components using image segmentation techniques focused primarily on nuclei (Funnell and Maysinger, 2006, Gladilin et al., 2008), and the most advanced structure-based cellular mechanics models to date utilize stacks of confocal photomicrographs of a cell to generate 3D model structures. There have been a small number of these types of models proposed in the last several years (Dailey et al., 2009, Slomka and Gefen, 2010), each of which have been important advances towards the development of a fully representative 3D model of single cell mechanics. However, none of those models has been constructed with entirely non-idealized geometries for all mechanically relevant components of a cell.
Few 3D single cell models have included any form of cytoskeletal elements inside the cells (Slomka and Gefen, 2010); yet even though these models represent a significant step towards reality, they still rely on the manual addition of a limited number of cytoskeletal components. There has not yet been a system put forth in the literature that is either fully automated or capable of reconstructing any elements of the cytoskeletal networks of cells in a representative manner. The goal of this study is to present such a fully automated cellular geometric reconstruction system based on 3D confocal microscopy images of single subconfluent cells.
Section snippets
Data acquisition: Cell culture, staining, and imaging
Primary rat aortic vascular smooth muscle cells (VSMCs) obtained from female Sprague Dawley rats are used in this study. The cells are cultured in high glucose Dulbecco’s Modified Eagle’s Medium (DMEM) (HyClone Laboratories, Logan, UT, USA) with an antibiotic solution of penicillin and streptomycin (HyClone Laboratories) added to a concentration 0.5%, and an antimycotic solution of amphotericin B (HyClone Laboratories) added to a concentration 0.5%. Cells are cultured in T75 cell culture
Cell imaging
The deconvoluted images of the cell used in the analysis are shown in Fig. 4 with separated actin and nucleus images shown on the left in each pane, and a distance between image planes of 200 nm. In addition, further examples of the segmentation technique on other cells are provided in the Supplemental material. Note that a portion the nucleus still clearly appears in the upper-most image planes whereas the actin network, though present, is much more difficult to distinguish. This is a result of
Potential applications
As mentioned above, the synthetic actin stress fiber networks generated in this study do not exactly match those of the original images they were generated from. This is intentional due to the fact that the geometries generated in this study are produced with the primary goal of utilization in finite element models for structural analysis. Therefore, it is imperative that the geometries be both an accurate representation of the original geometry yet sufficiently simple that the models can be
Summary and conclusions
This paper presents an automated method for generation of structural components of single cells based on 3D stacks of confocal microscope images for use in structural finite element analysis (Fig. 11). The major contribution of this study is the novel technique presented for generation of a representative actin stress fiber network.
Cell and nucleus boundaries are segmented using simple thresholding techniques. Generation of a representative actin stress fiber network is achieved by analyzing a
Acknowledgements
This study was funded in part by NIH K25 HL 092228, NSF EPS-0903795, and NSF CCF-0845593 grants.
References (41)
- et al.
Cytoskeletal bundle mechanics
Biophysical Journal
(2008) - et al.
Simultaneous use of class-i and class-ii force fields in CHARMM for solid–liquid multiphase simulation of protein–surface interaction
Biophysical Journal
(2009) - et al.
A bio-mechanical model for coupling cell contractility with focal adhesion formation
Journal of the Mechanics and Physics of Solids
(2008) - et al.
Models of morphogenesis: the mechanisms and mechanics of cell rearrangement
Current Opinion in Genetics & Development
(2004) - et al.
Mechanics of microtubules
Journal of Biomechanics
(2010) - et al.
A three-dimensional viscoelastic model for cell deformation with experimental verification
Biophysical Journal
(2003) - et al.
The cell as a material
Current Opinion in Cell Biology
(2007) A mechanics model of microtubule buckling in living cells
Journal of Biomechanics
(2008)- et al.
Mechanical and failure properties of single attached cells under compression
Journal of Biomechanics
(2005) - et al.
Rheological properties of the Eukaryotic cell cytoskeleton
Physics Reports
(2007)
Segmentation of fluorescence microscopy images for quantitative analysis of cell nuclear architecture
Biophysical Journal
Confocal microscopy-based three-dimensional cell-specific modeling for large deformation analyses in cellular mechanics
Journal of Biomechanics
Microtubules may harden or soften cells, depending on the extent of cell distension
Journal of Biomechanics
Cytoskeletal mechanics in airway smooth muscle cells
Respiratory Physiology & Neurobiology
Biomechanics and biophysics of cancer cells
Acta Materialia
How actin crosslinking and bundling proteins cooperate to generate an enhanced cell mechanical response
Biochemical and Biophysical Research Communications
Design of a custom angled abutment for dental implants using computer-aided design and nonlinear finite element analysis
Journal of Biomechanics
Cytoskeletal remodelling and slow dynamics in the living cell
Nature Materials
Editorial note: molecular imaging and mechanobiology
Cellular and Molecular Bioengineering
Evaluation of biological cell properties using dynamic indentation measurement
Physical Review E
Cited by (9)
A new method for reconstructing the 3D shape of single cells in fruit
2022, Food Research InternationalCitation Excerpt :According to Hertz contact theory, the radii of the curvature of two objects determines their force–displacement change law in the process of mutual contact action of the two objects (Li et al., 2017). If the estimated single-cell shape deviates significantly from the real cell shape (Wood et al., 2013), the single-cell mechanics characterized by Hertz contact theory will be inaccurate. As a result, the goal of this study is to propose a method for reconstructing the 3D shape of single cell in fruit so that a more real 3D morphology information of the fruit cells can be obtained and understood.
Extracting microtubule networks from superresolution single-molecule localization microscopy data
2017, Molecular Biology of the CellLA-iMageS: a software for elemental distribution bioimaging using LA-ICP-MS data
2016, Journal of CheminformaticsDeFiNe: An optimisation-based method for robust disentangling of filamentous networks
2015, Scientific Reports