Gradient of rigidity in the lamellipodia of migrating cells revealed by atomic force microscopy

Lipid membranes supported by polydimethylsiloxane substrates with designed geometry

The membrane curvature of cells and intracellular compartments continuously adapts to enable cells to perform vital functions, from cell division to signal trafficking. Understanding how membrane geometry affects these processes in vivo is challenging because of the biochemical and geometrical complexity as well as the short time and small length scales involved in cellular processes. By contrast, in vitro model membranes with engineered curvature would provide a versatile platform for this investigation and applications to biosensing and biocomputing. Here, we present a strategy that allows fabrication of lipid membranes with designed shape by combining 3D micro-printing and replica-molding lithography with polydimethylsiloxane to create curved micrometer-sized scaffolds with virtually any geometry. The resulting supported lipid membranes are homogeneous and fluid. We demonstrate the versatility of the system by fabricating structures of interesting combinations of mean and Gaussian curvature. We study the lateral phase separation and how local curvature influences the effective diffusion coefficient. Overall, we offer a bio-compatible platform for understanding curvature-dependent cellular processes and developing programmable bio-interfaces for living cells and nanostructures.

Multidepth quantitative analysis of liver cell viscoelastic properties: Fusion of nanoindentation and finite element modeling techniques

Liver cells are the basic functional unit of the liver. However, repeated or sustained injury leads to structural disorders of liver lobules, proliferation of fibrous tissue and changes in structure, thus increasing scar tissue. Cellular fibrosis affects tissue stiffness, shear force, and other cellular mechanical forces. Mechanical force characteristics can serve as important indicators of cell damage and cirrhosis. Atomic force microscopy (AFM) has been widely used to study cell surface mechanics. However, characterization of the deep mechanical properties inside liver cells remains an underdeveloped field. In this work, cell nanoindentation was combined with finite element analysis to simulate and analyze the mechanical responses of liver cells at different depths in vitro and their internal responses and stress diffusion distributions after being subjected to normal stress. The sensitivities of the visco-hyperelastic parameters of the finite element model to the effects of the peak force and equilibrium force were compared. The force curves of alcohol-damaged liver cells at different depths were measured and compared with those of undamaged liver cells. The inverse analysis method was used to simulate the finite element model in vitro. Changes in the parameters of the cell model after injury were explored and analyzed, and their potential for characterizing hepatocellular injury and related treatments was evaluated. RESEARCH HIGHLIGHTS: This study aims to establish an in vitro hyperelastic model of liver cells and analyze the mechanical changes of cells in vitro. An analysis method combining finite element analysis model and nanoindentation was used to obtain the key parameters of the model. The multi-depth mechanical differences and internal structural changes of injured liver cells were analyzed.

Dynamically Mapping the Topography and Stiffness of the Leading Edge of Migrating Cells Using AFM in Fast-QI Mode

Cell migration profoundly influences cellular function, often resulting in adverse effects in various pathologies including cancer metastasis. Directly assessing and quantifying the nanoscale dynamics of living cell structure and mechanics has remained a challenge. At the forefront of cell movement, the flat actin modules─the lamellipodium and the lamellum─interact to propel cell migration. The lamellipodium extends from the lamellum and undergoes rapid changes within seconds, making measurement of its stiffness a persistent hurdle. In this study, we introduce the fast-quantitative imaging (fast-QI) mode, demonstrating its capability to simultaneously map both the lamellipodium and the lamellum with enhanced spatiotemporal resolution compared with the classic quantitative imaging (QI) mode. Specifically, our findings reveal nanoscale stiffness gradients in the lamellipodium at the leading edge, where it appears to be slightly thinner and significantly softer than the lamellum. Additionally, we illustrate the fast-QI mode's accuracy in generating maps of height and effective stiffness through a streamlined and efficient processing of force-distance curves. These results underscore the potential of the fast-QI mode for investigating the role of motile cell structures in mechanosensing.

Geometry-Sensitive Protrusion Growth Directs Confined Cell Migration

The migratory dynamics of cells can be influenced by the complex microenvironment through which they move. It remains unclear how the motility machinery of confined cells responds and adapts to their microenvironment. Here, we propose a biophysical mechanism for a geometry-dependent coupling between cellular protrusions and the nucleus that leads to directed migration. We apply our model to geometry-guided cell migration to obtain insights into the origin of directed migration on asymmetric adhesive micropatterns and the polarization enhancement of cells observed under strong confinement. Remarkably, for cells that can choose between channels of different size, our model predicts an intricate dependence for cellular decision making as a function of the two channel widths, which we confirm experimentally.

Lamellipodia dynamics and microrheology in endothelial cell paracellular gap closure

Vascular endothelial cells (ECs) form a semipermeable barrier separating vascular contents from the interstitium, thereby regulating the movement of water and molecular solutes across small intercellular gaps, which are continuously forming and closing. Under inflammatory conditions, however, larger EC gaps form resulting in increased vascular leakiness to circulating fluid, proteins, and cells, which results in organ edema and dysfunction responsible for key pathophysiologic findings in numerous inflammatory disorders. In this study, we extend our earlier work examining the biophysical properties of EC gap formation and now address the role of lamellipodia, thin sheet-like membrane projections from the leading edge, in modulating EC spatial-specific contractile properties and gap closure. Micropillars, fabricated by soft lithography, were utilized to form reproducible paracellular gaps in human lung ECs. Using time-lapse imaging via optical microscopy, rates of EC gap closure and motility were measured with and without EC stimulation with the barrier-enhancing sphingolipid, sphingosine-1-phosphate. Peripheral ruffle formation was ubiquitous during gap closure. Kymographs were generated to quantitatively compare the lamellipodia dynamics of sphingosine-1-phosphate-stimulated and -unstimulated ECs. Utilizing atomic force microscopy, we characterized the viscoelastic behavior of EC lamellipodia. Our results indicate decreased stiffness and increased liquid-like behavior of expanding lamellipodia compared with regions away from the cellular edge (lamella and cell body) during EC gap closure, results in sync with the rapid kinetics of protrusion/retraction motion. We hypothesize this dissipative EC behavior during gap closure is linked to actomyosin cytoskeletal rearrangement and decreased cross-linking during lamellipodia expansion. In summary, these studies of the kinetic and mechanical properties of EC lamellipodia and ruffles at gap boundaries yield insights into the mechanisms of vascular barrier restoration and potentially a model system for examining the druggability of lamellipodial protein targets to enhance vascular barrier integrity.

Inferring local molecular dynamics from the global actin network structure: A case study of 2D synthetic branching actin networks

Cells rely on their cytoskeleton for key processes including division and directed motility. Actin filaments are a primary constituent of the cytoskeleton. Although actin filaments can create a variety of network architectures linked to distinct cell functions, the microscale molecular interactions that give rise to these macroscale structures are not well understood. In this work, we investigate the microscale mechanisms that produce different branched actin network structures using an iterative classification approach. First, we employ a simple yet comprehensive agent-based model that produces synthetic actin networks with precise control over the microscale dynamics. Then we apply machine learning techniques to classify actin networks based on measurable network density and geometry, identifying key mechanistic processes that lead to particular branched actin network architectures. Extensive computational experiments reveal that the most accurate method uses a combination of supervised learning based on network density and unsupervised learning based on network symmetry. This framework can potentially serve as a powerful tool to discover the molecular interactions that produce the wide variety of actin network configurations associated with normal development as well as pathological conditions such as cancer.

Switch of cell migration modes orchestrated by changes of three-dimensional lamellipodium structure and intracellular diffusion

Cell migration plays important roles in many biological processes, but how migrating cells orchestrate intracellular molecules and subcellular structures to regulate their speed and direction is still not clear. Here, by characterizing the intracellular diffusion and the three-dimensional lamellipodium structures of fish keratocyte cells, we observe a strong positive correlation between the intracellular diffusion and cell migration speed and, more importantly, discover a switching of cell migration modes with reversible intracellular diffusion variation and lamellipodium structure deformation. Distinct from the normal fast mode, cells migrating in the newly-found slow mode have a deformed lamellipodium with swollen-up front and thinned-down rear, reduced intracellular diffusion and compartmentalized macromolecule distribution in the lamellipodium. Furthermore, in turning cells, both lamellipodium structure and intracellular diffusion dynamics are also changed, with left-right symmetry breaking. We propose a mechanism involving the front-localized actin polymerization and increased molecular crowding in the lamellipodium to explain how cells spatiotemporally coordinate the intracellular diffusion dynamics and the lamellipodium structure in regulating their migrations.

Cholesterol and Sphingomyelin Polarize at the Leading Edge of Migrating Myoblasts and Involve Their Clustering in Submicrometric Domains

Myoblast migration is crucial for myogenesis and muscular tissue homeostasis. However, its spatiotemporal control remains elusive. Here, we explored the involvement of plasma membrane cholesterol and sphingolipids in this process. In resting C2C12 mouse myoblasts, those lipids clustered in sphingomyelin/cholesterol/GM1 ganglioside (SM/chol/GM1)- and cholesterol (chol)-enriched domains, which presented a lower stiffness than the bulk membrane. Upon migration, cholesterol and sphingomyelin polarized at the front, forming cholesterol (chol)- and sphingomyelin/cholesterol (SM/chol)-enriched domains, while GM1-enriched domains polarized at the rear. A comparison of domain proportion suggested that SM/chol- and GM1-enriched domains originated from the SM/chol/GM1-coenriched domains found at resting state. Modulation of domain proportion (through cholesterol depletion, combined or not with actin polymerization inhibition, or sphingolipid synthesis inhibition) revealed that the higher the chol- and SM/chol-enriched domains, the higher the myoblast migration. At the front, chol- and SM/chol-enriched domains were found in proximity with F-actin fibers and the lateral mobility of sphingomyelin in domains was specifically restricted in a cholesterol- and cytoskeleton-dependent manner while domain abrogation impaired F-actin and focal adhesion polarization. Altogether, we showed the polarization of cholesterol and sphingomyelin and their clustering in chol- and SM/chol-enriched domains with differential properties and roles, providing a mechanism for the spatial and functional control of myoblast migration.

Prediction of Golgi Polarity in Collectively Migrating Epithelial Cells Using Graph Neural Network

In the stationary epithelium, the Golgi apparatus assumes an apical position, above the cell nucleus. However, during wound healing and morphogenesis, as the epithelial cells start migrating, it relocalizes closer to the basal plane. On this plane, the position of Golgi with respect to the cell nucleus defines the organizational polarity of a migrating epithelial cell, which is crucial for an efficient collective migration. Yet, factors influencing the Golgi polarity remain elusive. Here, we constructed a graph neural network-based deep learning model to systematically analyze the dependency of Golgi polarity on multiple geometric and physical factors. In spite of the complexity of a migrating epithelial monolayer, our simple model was able to predict the Golgi polarity with 75% accuracy. Moreover, the model predicted that Golgi polarity predominantly correlates with the orientation of maximum principal stress. Finally, we found that this correlation operates locally since progressive coarsening of the stress field over multiple cell-lengths reduced the stress polarity-Golgi polarity correlation as well as the predictive accuracy of the neural network model. Taken together, our results demonstrate that graph neural networks could be a powerful tool toward understanding how different physical factors influence collective cell migration. They also highlight a previously unknown role of physical cues in defining the intracellular organization.

Effective Force Generation During Mammalian Cell Migration Under Different Molecular and Physical Mechanisms

We have developed much understanding of actin-driven cell migration and the forces that propel cell motility. However, fewer studies focused on estimating the effective forces generated by migrating cells. Since cells in vivo are exposed to complex physical environments with various barriers, understanding the forces generated by cells will provide insights into how cells manage to navigate challenging environments. In this work, we use theoretical models to discuss actin-driven and water-driven cell migration and the effect of cell shapes on force generation. The results show that the effective force generated by actin-driven cell migration is proportional to the rate of actin polymerization and the strength of focal adhesion; the energy source comes from the actin polymerization against the actin network pressure. The effective force generated by water-driven cell migration is proportional to the rate of active solute flux and the coefficient of external hydraulic resistance; the energy sources come from active solute pumping against the solute concentration gradient. The model further predicts that the actin network distribution is mechanosensitive and the presence of globular actin helps to establish a biphasic cell velocity in the strength of focal adhesion. The cell velocity and effective force generation also depend on the cell shape through the intracellular actin flow field.

The effects of measurement parameters on the cancerous cell nucleus characterization by atomic force microscopy in vitro

Cancer is now responsible for the major leading cause of death worldwide. It is noteworthy that lung cancer has been recognised as the highest incidence (11.6%) and mortality (18.4%) for combined sexes among a variety of cancer diseases. Therefore, it is of great value to investigate the mechanical properties of lung cancerous cells for early diagnosis. This paper focus on the influence of measurement parameters on the measured central Young's moduli of single live A549 cell in vitro based on the force spectroscopy mode of atomic force microscopy (AFM). The effects of the measurement parameters on the measured central Young's moduli were analysed by fitting the force–depth curves utilising the Sneddon model. The results revealed that the Young's moduli of A549 cells increased with the larger indentation force, higher indentation speed, less retraction time, deeper Z length and lower purity percentage of serum. The Young's moduli of cells increased first and then decreased with the increasing dwell time. Hence, this research may have potential significance to provide reference for the standardised detection of a single cancerous cell in vitro using AFM methodologies.

Cancer is now responsible for the majority leading cause of death worldwide and it is noteworthy that lung cancer has been recognised as the highest incidence (11.6%) and mortality (18.4%) for combined sexes among a variety of cancer diseases. Therefore, it is of great value to investigate the mechanical properties of lung cancerous cells for early diagnosis. This paper primarily investigated the morphological properties and the influence of measurement parameters on the measured local elastic moduli of single live A549 cell in vitro using the AFM‐based force spectroscopy mode. In practice, there are many factors for incorrect or inaccurate experimental results using AFM to measure the characteristics of live cells, such as non‐homogeneous nature of cells, probe geometry and size, mechanical analysis model, substrate stiffness and different measurement parameters. The various measurement parameters have become the huge impact factor to influence the measurement result. Hence, this research may have potential significance to provide reference for the standardised detection of a single cancerous cell in vitro using AFM methodologies.

Leading edge maintenance in migrating cells is an emergent property of branched actin network growth

Animal cell migration is predominantly driven by the coordinated, yet stochastic, polymerization of thousands of nanometer-scale actin filaments across micron-scale cell leading edges. It remains unclear how such inherently noisy processes generate robust cellular behavior. We employed high-speed imaging of migrating neutrophil-like HL-60 cells to explore the fine-scale shape fluctuations that emerge and relax throughout the process of leading edge maintenance. We then developed a minimal stochastic model of the leading edge that reproduces this stable relaxation behavior. Remarkably, we find lamellipodial stability naturally emerges from the interplay between branched actin network growth and leading edge shape – with no additional feedback required – based on a synergy between membrane-proximal branching and lateral spreading of filaments. These results thus demonstrate a novel biological noise-suppression mechanism based entirely on system geometry. Furthermore, our model suggests that the Arp2/3-mediated ~70–80° branching angle optimally smooths lamellipodial shape, addressing its long-mysterious conservation from protists to mammals.

In every human cell, there are tens of millions of proteins which work together to control everything from the cell’s shape to its behavior. One of the most abundant proteins is actin, which organizes itself into filaments that mechanically support the cell and help it to move.

These filaments are very dynamic, with individual actin molecules constantly being added or removed. This allows the cell to build large structures with distinct shapes and properties. Many motile cells, for example, have a structure called a lamellipodium which protrudes at their ‘leading edge’ and pushes them forward. The lamellipodium has a very robust shape that does not vary much between different cell types, or change significantly as cells migrate. But how the tens of thousands of actin molecules inside the lamellipodium organize themselves into this large, stable structure is not fully understood.

To investigate, Garner and Theriot used high-speed video microscopy to track the shape of human cells cultured in the laboratory. As the cells crawled along a glass surface, their leading edge undulated like strings being plucked on a guitar. A computer simulation showed that these ripples can be caused by filaments randomly adding and removing actin molecules.

While these random movements could destabilize the structure of the leading edge, the simulation suggests that another aspect of actin filament growth smooths out any fluctuations in the lamellipodium’s shape. Actin networks in the lamellipodium have a branched configuration, with new strands emerging off each other at an angle like branches in a tree. Garner and Theriot found that the specific angle in which new filaments are added smooths out the lamellipodium’s shape, which may explain why this geometry has persisted throughout evolution.

These findings suggest that the way in which actin filaments join together helps to maintain the shape of large cellular structures. In the future, scientists could use this design principle to build molecular machines that can self-organize into microstructures. These engineered constructs could be used to modulate the activity of living cells that have been damaged by disease.

Viscoelasticity, Like Forces, Plays a Role in Mechanotransduction

Viscoelasticity and its alteration in time and space has turned out to act as a key element in fundamental biological processes in living systems, such as morphogenesis and motility. Based on experimental and theoretical findings it can be proposed that viscoelasticity of cells, spheroids and tissues seems to be a collective characteristic that demands macromolecular, intracellular component and intercellular interactions. A major challenge is to couple the alterations in the macroscopic structural or material characteristics of cells, spheroids and tissues, such as cell and tissue phase transitions, to the microscopic interferences of their elements. Therefore, the biophysical technologies need to be improved, advanced and connected to classical biological assays. In this review, the viscoelastic nature of cytoskeletal, extracellular and cellular networks is presented and discussed. Viscoelasticity is conceptualized as a major contributor to cell migration and invasion and it is discussed whether it can serve as a biomarker for the cells' migratory capacity in several biological contexts. It can be hypothesized that the statistical mechanics of intra- and extracellular networks may be applied in the future as a powerful tool to explore quantitatively the biomechanical foundation of viscoelasticity over a broad range of time and length scales. Finally, the importance of the cellular viscoelasticity is illustrated in identifying and characterizing multiple disorders, such as cancer, tissue injuries, acute or chronic inflammations or fibrotic diseases.

Axial Scanning Metal-Induced Energy Transfer Microscopy for Extended Range Nanometer-Sectioning Cell Imaging

Nanometer-sectioning optical microscopy has become an indispensable tool in membrane-related biomedical studies. Finally, many nanometer-sectioning imaging schemes, such as variable-angle total internal reflection fluorescence microscopy, metal-induced energy transfer (MIET) imaging, and supercritical-angle fluorescence microscopy have been introduced. However, these methods can measure a single layer of molecules, and the measurement ranges are below 100 nm, which is not large enough to cover the thickness of lamellipodium. This paper proposes an optical imaging scheme that can identify the axial locations of two layers of molecules with an extended measurement range and a nanometer-scale precision by using MIET, axial focal plane scanning, and biexponential analysis in fluorescence lifetime imaging microscopy. The feasibility of the proposed method is demonstrated by measuring an artificial sample of a known structure and the lamellipodium of a human aortic endothelial cell whose thickness ranges from 100 to 450 nm with 18.3 nm precision.

Discrete mechanical model of lamellipodial actin network implements molecular clutch mechanism and generates arcs and microspikes

Mechanical forces, actin filament turnover, and adhesion to the extracellular environment regulate lamellipodial protrusions. Computational and mathematical models at the continuum level have been used to investigate the molecular clutch mechanism, calculating the stress profile through the lamellipodium and around focal adhesions. However, the forces and deformations of individual actin filaments have not been considered while interactions between actin networks and actin bundles is not easily accounted with such methods. We develop a filament-level model of a lamellipodial actin network undergoing retrograde flow using 3D Brownian dynamics. Retrograde flow is promoted in simulations by pushing forces from the leading edge (due to actin polymerization), pulling forces (due to molecular motors), and opposed by viscous drag in cytoplasm and focal adhesions. Simulated networks have densities similar to measurements in prior electron micrographs. Connectivity between individual actin segments is maintained by permanent and dynamic crosslinkers. Remodeling of the network occurs via the addition of single actin filaments near the leading edge and via filament bond severing. We investigated how several parameters affect the stress distribution, network deformation and retrograde flow speed. The model captures the decrease in retrograde flow upon increase of focal adhesion strength. The stress profile changes from compression to extension across the leading edge, with regions of filament bending around focal adhesions. The model reproduces the observed reduction in retrograde flow speed upon exposure to cytochalasin D, which halts actin polymerization. Changes in crosslinker concentration and dynamics, as well as in the orientation pattern of newly added filaments demonstrate the model's ability to generate bundles of filaments perpendicular (actin arcs) or parallel (microspikes) to the protruding direction.

Cell Surface Binding and Lipid Interactions behind Chemotherapy-Drug-Induced Ion Pore Formation in Membranes

Chemotherapy drugs (CDs) disrupt the lipid membrane’s insulation properties by inducing stable ion pores across bilayer membranes. The underlying molecular mechanisms behind pore formation have been revealed in this study using several methods that confirm molecular interactions and detect associated energetics of drugs on the cell surface in general and in lipid bilayers in particular. Liposome adsorption and cell surface binding of CD colchicine has been demonstrated experimentally. Buffer dissolved CDs were considerably adsorbed in the incubated phospholipid liposomes, measured using the patented ‘direct detection method’. The drug adsorption process is regulated by the membrane environment, demonstrated in cholesterol-containing liposomes. We then detailed the phenomenology and energetics of the low nanoscale dimension cell surface (membrane) drug distribution, using atomic force microscopy (AFM) imaging what addresses the surface morphology and measures adhesion force (reducible to adhesive energy). Liposome adsorption and cell surface binding data helped model the cell surface drug distribution. The underlying molecular interactions behind surface binding energetics of drugs have been addressed in silico numerical computations (NCs) utilizing the screened Coulomb interactions among charges in a drug–drug/lipid cluster. Molecular dynamics (MD) simulations of the CD-lipid complexes detected primarily important CD-lipid electrostatic and van der Waals (vdW) interaction energies. From the energetics point of view, both liposome and cell surface membrane adsorption of drugs are therefore obvious findings. Colchicine treated cell surface AFM images provide a few important phenomenological conclusions, such as drugs bind generally with the cell surface, bind independently as well as in clusters of various sizes in random cell surface locations. The related adhesion energy decreases with increasing drug cluster size before saturating for larger clusters. MD simulation detected electrostatic and vdW and NC-derived charge-based interactions explain molecularly of the cause of cell surface binding of drugs. The membrane binding/association of drugs may help create drug–lipid complexes with specific energetics and statistically lead to the creation of ion channels. We reveal here crucial molecular understanding and features of the pore formation inside lipid membranes that may be applied universally for most of the pore-forming existing agents and novel candidate drugs.

Atomic Force Microscopy for Live-Cell and Hydrogel Measurement

Atomic force microscopy (AFM) has emerged as a popular method for determining the mechanical properties of cells, their components, and biomaterials. Here, we describe AFM setup and application to obtain stiffness measurements from single indentations for hydrogels and myofibroblasts.

Predictive assembling model reveals the self-adaptive elastic properties of lamellipodial actin networks for cell migration

Branched actin network supports cell migration through extracellular microenvironments. However, it is unknown how intracellular proteins adapt the elastic properties of the network to the highly varying extracellular resistance. Here we develop a three-dimensional assembling model to simulate the realistic self-assembling process of the network by encompassing intracellular proteins and their dynamic interactions. Combining this multiscale model with finite element method, we reveal that the network can not only sense the variation of extracellular resistance but also self-adapt its elastic properties through remodeling with intracellular proteins. Such resistance-adaptive elastic behaviours are versatile and essential in supporting cell migration through varying extracellular microenvironments. The bending deformation mechanism and anisotropic Poisson's ratios determine why lamellipodia persistently evolve into sheet-like structures. Our predictions are confirmed by published experiments. The revealed self-adaptive elastic properties of the networks are also applicable to the endocytosis, phagocytosis, vesicle trafficking, intracellular pathogen transport and dendritic spine formation.

Probing arsenic trioxide (ATO) treated leukemia cell elasticities using atomic force microscopy

Conventional analytical techniques allow for the diagnosis of leukemia, blood and bone marrow cancers, as well as their classification into the different subtypes. However, a better understanding of the cancer treatment through cell apoptosis staging is still required. Evaluation of the timeline and responses of acute promyelocytic leukemia (APL) cells to the arsenic trioxide (ATO) treatment is essential for determining the oral dosage in leukemia prognosis. Here, an Atomic Force Microscopy (AFM) indentation approach has been used to evaluate the mechanical responses of cellular responses of APL cells to ATO treatment, alongside well-established cell viability assays, as a novel method to determine the impact of drugs. In addition, cell morphology was quantified to monitor cellular apoptosis. Viability, morphology and elasticity changes of NB4 cells (derived from Human APL patients) were correlated to different time courses of the ATO treatments. Unveiling the relationships among structural, morphological and nanomechanical properties in response to ATO drug treatment promises to pave the way for novel diagnostic tools for drug screening and for a better understanding of the specific physical and biological effects of drugs on diseased cells.

Near Simultaneous Laser Scanning Confocal and Atomic Force Microscopy (Conpokal) on Live Cells

Techniques available for micro- and nano-scale mechanical characterization have exploded in the last few decades. From further development of the scanning and transmission electron microscope, to the invention of atomic force microscopy, and advances in fluorescent imaging, there have been substantial gains in technologies that enable the study of small materials. Conpokal is a portmanteau that combines confocal microscopy with atomic force microscopy (AFM), where a probe "pokes" the surface. Although each technique is extremely effective for the qualitative and/or quantitative image collection on their own, Conpokal provides the capability to test with blended fluorescence imaging and mechanical characterization. Designed for near simultaneous confocal imaging and atomic force probing, Conpokal facilitates experimentation on live microbiological samples. The added insight from paired instrumentation provides co-localization of measured mechanical properties (e.g., elastic modulus, adhesion, surface roughness) by AFM with subcellular components or activity observable through confocal microscopy. This work provides a step by step protocol for the operation of laser scanning confocal and atomic force microscopy, simultaneously, to achieve same cell, same region, confocal imaging, and mechanical characterization.

Potential antitumor activity of digitoxin and user-designed analog administered to human lung cancer cells

BACKGROUND

Cardiac glycosides (CGs), such as digitoxin, are traditionally used for treatment of congestive heart failure; recently they also gained attention for their anticancer properties. Previous studies showed that digitoxin and a synthetic L-sugar monosaccharide analog treatment decreases cancer cell proliferation, increases apoptosis, and pro-adhesion abilities; however, no reports are available on their potential to alter lung cancer cell cytoskeleton structure and reduce migratory ability. Herein, we investigated the anticancer effects of digitoxin and its analog, digitoxigenin-α-L-rhamnoside (D6MA), to establish whether cytoskeleton reorganization and reduced motility are drug-induced cellular outcomes.

METHODS

We treated non-small cell lung carcinoma cells (NSCLCs) with sub-therapeutic, therapeutic, and toxic concentrations of digitoxin and D6MA respectively, followed by both single point and real-time assays to evaluate changes in cellular gene and protein expression, adhesion, elasticity, and migration.

RESULTS

Digitoxin and D6MA induced a decrease in matrix metalloproteinases expression via altered focal adhesion signaling and a suppression of the phosphoinositide 3-kinases / protein kinase B pathway which lead to enhanced adhesion, altered elasticity, and reduced motility of NSCLCs. Global gene expression analysis identified dose-dependent changes to nuclear factor kappa-light-chain-enhancer, epithelial tumor, and microtubule dynamics signaling.

CONCLUSIONS

Our study demonstrates that digitoxin and D6MA can target antitumor signaling pathways to alter NSCLC cytoskeleton and migratory ability to thus potentially reduce their tumorigenicity.

SIGNIFICANCE

Discovering signaling pathways that control cancer's cell phenotype and how such pathways are affected by CG treatment will potentially allow for active usage of synthetic CG analogs as therapeutic agents in advanced lung conditions.

Direct monitoring of drug-induced mechanical response of individual cells by atomic force microscopy

Mechanical characteristics of individual cells play a vital role in many biological processes and are considered as indicators of the cells' states. Disturbances including methyl-β-cyclodextrin (MβCD) and cytochalasin D (cytoD) are known to significantly affect the state of cells, but little is known about the real-time response of single cells to these drugs in their physiological condition. Here, nanoindentation-based atomic force microscopy (AFM) was used to measure the elasticity of human embryonic kidney cells in the presence and absence of these pharmaceuticals. The results showed that depletion of cholesterol in the plasma membrane with MβCD resulted in cell stiffening whereas depolymerization of the actin cytoskeleton by cytoD resulted in cell softening. Using AFM for real-time measurements, we observed that cells mechanically responded right after these drugs were added. In more detail, the cell´s elasticity suddenly increased with increasing instability upon cholesterol extraction while it is rapidly decreased without changing cellular stability upon depolymerizing actin cytoskeleton. These results demonstrated that actin cytoskeleton and cholesterol contributed differently to the cell mechanical characteristics.

Finite element analysis reveals an important role for cell morphology in response to mechanical compression

Mechanical loading naturally controls cell phenotype, development, motility and various other biological functions; however, prolonged or substantial loading can cause cell damage and eventual death. Loading-induced mechanobiological and mechanostructural responses of different cell types affect their morphology and the internal architecture and the mechanics of the cellular components. Using single, mesenchymal stem cells, we have developed a cell-specific three-dimensional finite-element model; cell models were developed from phase-contrast microscopy images. This allowed us to evaluate the mechanostructural response of the naturally occurring variety of cell morphologies to increase sustained compressive loading. We focus on the morphology of the cytoplasm and the nucleus, as the main mechanically responsive elements, and evaluate formation of tensional strains and area changes in cells undergoing increasing uniaxial compressions. Here, we study mesenchymal stem cells as a model, due to their important role in tissue engineering and regenerative medicine; the method and findings are, however, applicable to any cell type. We observe variability in the cell responses to compression, which correlate directly with the morphology of the cells. Specifically, in cells with or without elongated protrusions (i.e., lamellipodia) tensional strains were, respectively, distributed mostly in the thin extensions or concentrated around the stiff nucleus. Thus, through cell-specific computational modeling of mechanical loading we have identified an underlying cause for stiffening (by actin recruitment) along the length of lamellipodia as well as a role for cell morphology in inducing cell-to-cell variability in mechanostructural response to loading.

Flow-enhanced priming of hESCs through H2B acetylation and chromatin decondensation

BACKGROUND

Distinct mechanical stimuli are known to manipulate the behaviors of embryonic stem cells (ESCs). Fundamental rationale of how ESCs respond to mechanical forces and the potential biological effects remain elusive. Here we conducted the mechanobiological study for hESCs upon mechanomics analysis to unravel typical mechanosensitive processes on hESC-specific fluid shear.

METHODS

hESC line H1 was subjected to systematically varied shear flow, and mechanosensitive proteins were obtained by mass spectrometry (MS) analysis. Then, function enrichment analysis was performed to identify the enriched gene sets. Under a steady shear flow of 1.1 Pa for 24 h, protein expressions were further detected using western blotting (WB), quantitative real-time PCR (qPCR), and immunofluorescence (IF) staining. Meanwhile, the cells were treated with 200 nM trichostatin (TSA) for 1 h as positive control to test chromatin decondensation. Actin, DNA, and RNA were then visualized with TRITC-labeled phalloidin, Hoechst 33342, and SYTO® RNASelect™ green fluorescent cell stain (Life Technologies), respectively. In addition, cell stiffness was determined with atomic force microscopy (AFM) and annexin V-PE was used to determine the apoptosis with a flow cytometer (FCM).

RESULTS

Typical mechanosensitive proteins were unraveled upon mechanomics analysis under fluid shear related to hESCs in vivo. Functional analyses revealed significant alterations in histone acetylation, nuclear size, and cytoskeleton for hESC under shear flow. Shear flow was able to induce H2B acetylation and nuclear spreading by CFL2/F-actin cytoskeletal reorganization. The resulting chromatin decondensation and a larger nucleus readily accommodate signaling molecules and transcription factors.

CONCLUSIONS

Shear flow regulated chromatin dynamics in hESCs via cytoskeleton and nucleus alterations and consolidated their primed state.

Measuring the Elastic Properties of Living Cells

Cell's elasticity is an integrative parameter summarizing the biophysical outcome of many known and unknown cellular processes. This includes intracellular signaling, cytoskeletal activity, changes of cell volume and morphology, and many others. Not only intracellular processes defines a cell's elasticity but also environmental factors like their biochemical and biophysical surrounding. Therefore, cell mechanics represents a comprehensive variable of life. A cell in its standard conditions shows variabilities of biochemical and biophysical processes resulting in a certain range of cell's elasticity. Changes of the standard conditions, endogenously or exogenously induced, are frequently paralleled by changes of cell elasticity. Therefore cell elasticity could serve as parameter to characterize different states of a cell. Atomic force microscopy (AFM) combines high spatial resolution with very high force sensitivity and allows investigating mechanical properties of living cells under physiological conditions. However, elastic moduli reported in the literature showed a large variability, sometimes by an order of magnitude (or even more) for the same cell type assessed in different labs. Clearly, a prerequisite for the use of cell elasticity to describe the actual cell status is a standardized procedure that allows obtaining comparable values of a cell independent from the instrument, from the lab and operator. Biologically derived variations of elasticity could not be reduced due to the nature of living cells but technically and methodologically derived variations could be minimized by a standardized procedure.This chapter provides a Standardized Nanomechanical AFM Procedure (SNAP) that reduces strongly the variability of results obtained on soft samples and living cells by a reliable method to calibrate AFM cantilevers.

Lamellipodial wrinkles in fish keratocytes as markers of imperfect coordination between extension and retraction during cell migration

Cell migration involves the precise coordination between extension at the front of the cell and retraction at the rear. This coordination is particularly evident in fast moving cells such as fish keratocytes, where it leads to highly stable gliding motion, propelled at the front by broad, 0.1-0-2 μm thick lamellipodia. Transient uncoupling between extension and retraction can occur if the rear is temporarily stuck, which might eventually lead to cell shape instabilities. We have frequently observed in fish keratocytes the presence of lamellipodial radial wrinkles, detected by confocal, scanning electron and side-view microscopy as folds in the lamellipodium up to 2 μm in height. Using a linear finite elements analysis, we simulated the displacement of cells either with perfect coordination between extension and retraction or with the rear transiently stuck while the front continues extending, and we observed that in this last condition compression stresses arise in the lamellipodium which predict the formation of the observed pattern of lamellipodial wrinkles. In support of the numerical modeling findings, we observed that the transient halting of retraction at the rear using micromanipulation induced the formation of lamellipodial wrinkles in previously flat lamellipodia. The obtained results suggest that the conspicuous lamellipodial wrinkles observed in migrating fish keratocytes are the product of transient imbalances between front and rear displacements, and are therefore useful markers of the short scale dynamics of extension and retraction coordination during cell migration.

Actin Turnover in Lamellipodial Fragments

Actin turnover is the central driving force underlying lamellipodial motility. The molecular components involved are largely known, and their properties have been studied extensively in vitro. However, a comprehensive picture of actin turnover in vivo is still missing. We focus on fragments from fish epithelial keratocytes, which are essentially stand-alone motile lamellipodia. The geometric simplicity of the fragments and the absence of additional actin structures allow us to characterize the spatiotemporal lamellipodial actin organization with unprecedented detail. We use fluorescence recovery after photobleaching, fluorescence correlation spectroscopy, and extraction experiments to show that about two-thirds of the lamellipodial actin diffuses in the cytoplasm with nearly uniform density, whereas the rest forms the treadmilling polymer network. Roughly a quarter of the diffusible actin pool is in filamentous form as diffusing oligomers, indicating that severing and debranching are important steps in the disassembly process generating oligomers as intermediates. The remaining diffusible actin concentration is orders of magnitude higher than the in vitro actin monomer concentration required to support the observed polymerization rates, implying that the majority of monomers are transiently kept in a non-polymerizable "reserve" pool. The actin network disassembles and reassembles throughout the lamellipodium within seconds, so the lamellipodial network turnover is local. The diffusible actin transport, on the other hand, is global: actin subunits typically diffuse across the entire lamellipodium before reassembling into the network. This combination of local network turnover and global transport of dissociated subunits through the cytoplasm makes actin transport robust yet rapidly adaptable and amenable to regulation.

Standardized Nanomechanical Atomic Force Microscopy Procedure (SNAP) for Measuring Soft and Biological Samples

We present a procedure that allows a reliable determination of the elastic (Young's) modulus of soft samples, including living cells, by atomic force microscopy (AFM). The standardized nanomechanical AFM procedure (SNAP) ensures the precise adjustment of the AFM optical lever system, a prerequisite for all kinds of force spectroscopy methods, to obtain reliable values independent of the instrument, laboratory and operator. Measurements of soft hydrogel samples with a well-defined elastic modulus using different AFMs revealed that the uncertainties in the determination of the deflection sensitivity and subsequently cantilever's spring constant were the main sources of error. SNAP eliminates those errors by calculating the correct deflection sensitivity based on spring constants determined with a vibrometer. The procedure was validated within a large network of European laboratories by measuring the elastic properties of gels and living cells, showing that its application reduces the variability in elastic moduli of hydrogels down to 1%, and increased the consistency of living cells elasticity measurements by a factor of two. The high reproducibility of elasticity measurements provided by SNAP could improve significantly the applicability of cell mechanics as a quantitative marker to discriminate between cell types and conditions.

Dynamical Modeling and Analysis of Viscoelastic Properties of Single Cells

A single cell can be regarded as a complex network that contains thousands of overlapping signaling pathways. The traditional methods for describing the dynamics of this network are extremely complicated. The mechanical properties of a cell reflect the cytoskeletal structure and composition and are closely related to the cellular biological functions and physiological activities. Therefore, modeling the mechanical properties of single cells provides the basis for analyzing and controlling the cellular state. In this study, we developed a dynamical model with cellular viscoelasticity properties as the system parameters to describe the stress-relaxation phenomenon of a single cell indented by an atomic force microscope (AFM). The system order and parameters were identified and analyzed. Our results demonstrated that the parameters identified using this model represent the cellular mechanical elasticity and viscosity and can be used to classify cell types.

Multi-scale computational study of the mechanical regulation of cell mitotic rounding in epithelia

Mitotic rounding during cell division is critical for preventing daughter cells from inheriting an abnormal number of chromosomes, a condition that occurs frequently in cancer cells. Cells must significantly expand their apical area and transition from a polygonal to circular apical shape to achieve robust mitotic rounding in epithelial tissues, which is where most cancers initiate. However, how cells mechanically regulate robust mitotic rounding within packed tissues is unknown. Here, we analyze mitotic rounding using a newly developed multi-scale subcellular element computational model that is calibrated using experimental data. Novel biologically relevant features of the model include separate representations of the sub-cellular components including the apical membrane and cytoplasm of the cell at the tissue scale level as well as detailed description of cell properties during mitotic rounding. Regression analysis of predictive model simulation results reveals the relative contributions of osmotic pressure, cell-cell adhesion and cortical stiffness to mitotic rounding. Mitotic area expansion is largely driven by regulation of cytoplasmic pressure. Surprisingly, mitotic shape roundness within physiological ranges is most sensitive to variation in cell-cell adhesivity and stiffness. An understanding of how perturbed mechanical properties impact mitotic rounding has important potential implications on, amongst others, how tumors progressively become more genetically unstable due to increased chromosomal aneuploidy and more aggressive.

The role of myosin II in glioma invasion: A mathematical model

Gliomas are malignant tumors that are commonly observed in primary brain cancer. Glioma cells migrate through a dense network of normal cells in microenvironment and spread long distances within brain. In this paper we present a two-dimensional multiscale model in which a glioma cell is surrounded by normal cells and its migration is controlled by cell-mechanical components in the microenvironment via the regulation of myosin II in response to chemoattractants. Our simulation results show that the myosin II plays a key role in the deformation of the cell nucleus as the glioma cell passes through the narrow intercellular space smaller than its nuclear diameter. We also demonstrate that the coordination of biochemical and mechanical components within the cell enables a glioma cell to take the mode of amoeboid migration. This study sheds lights on the understanding of glioma infiltration through the narrow intercellular spaces and may provide a potential approach for the development of anti-invasion strategies via the injection of chemoattractants for localization.

Stiffness of pancreatic cancer cells is associated with increased invasive potential

Metastasis is a fundamentally physical process in which cells are required to deform through narrow gaps as they invade surrounding tissues and transit to distant sites. In many cancers, more invasive cells are more deformable than less invasive cells, but the extent to which mechanical phenotype, or mechanotype, can predict disease aggressiveness in pancreatic ductal adenocarcinoma (PDAC) remains unclear. Here we investigate the invasive potential and mechanical properties of immortalized PDAC cell lines derived from primary tumors and a secondary metastatic site, as well as noncancerous pancreatic ductal cells. To investigate how invasive behavior is associated with cell mechanotype, we flow cells through micron-scale pores using parallel microfiltration and microfluidic deformability cytometry; these results show that the ability of PDAC cells to passively transit through pores is only weakly correlated with their invasive potential. We also measure the Young's modulus of pancreatic ductal cells using atomic force microscopy, which reveals that there is a strong association between cell stiffness and invasive potential in PDAC cells. To determine the molecular origins of the variability in mechanotype across our PDAC cell lines, we analyze RNAseq data for genes that are known to regulate cell mechanotype. Our results show that vimentin, actin, and lamin A are among the most differentially expressed mechanoregulating genes across our panel of PDAC cell lines, as well as a cohort of 38 additional PDAC cell lines. We confirm levels of these proteins across our cell panel using immunoblotting, and find that levels of lamin A increase with both invasive potential and Young's modulus. Taken together, we find that stiffer PDAC cells are more invasive than more compliant cells, which challenges the paradigm that decreased cell stiffness is a hallmark of metastatic potential.

Cancer cells become less deformable and more invasive with activation of β-adrenergic signaling

Invasion by cancer cells is a crucial step in metastasis. An oversimplified view in the literature is that cancer cells become more deformable as they become more invasive. β-adrenergic receptor (βAR) signaling drives invasion and metastasis, but the effects on cell deformability are not known. Here, we show that activation of β-adrenergic signaling by βAR agonists reduces the deformability of highly metastatic human breast cancer cells, and that these stiffer cells are more invasive in vitro We find that βAR activation also reduces the deformability of ovarian, prostate, melanoma and leukemia cells. Mechanistically, we show that βAR-mediated cell stiffening depends on the actin cytoskeleton and myosin II activity. These changes in cell deformability can be prevented by pharmacological β-blockade or genetic knockout of the β₂-adrenergic receptor. Our results identify a β₂-adrenergic-Ca²⁺-actin axis as a new regulator of cell deformability, and suggest that the relationship between cell mechanical properties and invasion might be dependent on context.

Nano-mechanical single-cell sensing of cell-matrix contacts

Extracellular protein matrices provide a rigidity interface exhibiting nano-mechanical cues that guide cell growth and proliferation. Cells sense such cues using actin-rich filopodia extensions which encourage favourable cell-matrix contacts to recruit more actin-mediated local forces into forming stable focal adhesions. A challenge remains in identifying and measuring these local cellular forces and in establishing empirical relationships between them, cell adhesion and filopodia formation. Here we investigate such relationships using a micromanipulation system designed to operate at the time scale of focal contact dynamics, with the sample frequency of a force probe being 0.1 ms, and to apply and measure forces at nano-to-micro Newton ranges for individual mammalian cells. We explore correlations between cell biomechanics, cell-matrix attachment forces and the spread areas of adhered cells as well as their relative dependence on filopodia formation using synthetic protein matrices with a proven ability to induce enhanced filopodia numbers in adherent cells. This study offers a basis for engineering exploitable cell-matrix contacts in situ at the nanoscale and single-cell levels.

A novel approach for extracting viscoelastic parameters of living cells through combination of inverse finite element simulation and Atomic Force Microscopy

Dynamic mechanical behaviour of living cells has been described by viscoelasticity. However, quantitation of the viscoelastic parameters for living cells is far from sophisticated. In this paper, combining inverse finite element (FE) simulation with Atomic Force Microscope characterization, we attempt to develop a new method to evaluate and acquire trustworthy viscoelastic index of living cells. First, influence of the experiment parameters on stress relaxation process is assessed using FE simulation. As suggested by the simulations, cell height has negligible impact on shape of the force-time curve, i.e. the characteristic relaxation time; and the effect originates from substrate can be totally eliminated when stiff substrate (Young's modulus larger than 3 GPa) is used. Then, so as to develop an effective optimization strategy for the inverse FE simulation, the parameters sensitivity evaluation is performed for Young's modulus, Poisson's ratio, and characteristic relaxation time. With the experiment data obtained through typical stress relaxation measurement, viscoelastic parameters are extracted through the inverse FE simulation by comparing the simulation results and experimental measurements. Finally, reliability of the acquired mechanical parameters is verified with different load experiments performed on the same cell.

Spatial distribution of filament elasticity determines the migratory behaviors of a cell

Any cellular response leading to morphological changes is highly tuned to balance the force generated from structural reorganization, provided by actin cytoskeleton. Actin filaments serve as the backbone of intracellular force, and transduce external mechanical signal via focal adhesion complex into the cell. During migration, cells not only undergo molecular changes but also rapid mechanical modulation. Here we focus on determining, the role of spatial distribution of mechanical changes of actin filaments in epithelial, mesenchymal, fibrotic and cancer cells with non-migration, directional migration, and non-directional migration behaviors using the atomic force microscopy. We found 1) non-migratory cells only generated one type of filament elasticity, 2) cells generating spatially distributed two types of filament elasticity showed directional migration, and 3) pathologic cells that autonomously generated two types of filament elasticity without spatial distribution were actively migrating non-directionally. The demonstration of spatial regulation of filament elasticity of different cell types at the nano-scale highlights the coupling of cytoskeletal function with physical characters at the sub-cellular level, and provides new research directions for migration related disease.

Qualitative and Quantitative Analysis of ROS-Mediated Oridonin-Induced Oesophageal Cancer KYSE-150 Cell Apoptosis by Atomic Force Microscopy

High levels of intracellular reactive oxygen species (ROS) in cells is recognized as one of the major causes of cancer cell apoptosis and has been developed into a promising therapeutic strategy for cancer therapy. However, whether apoptosis associated biophysical properties of cancer cells are related to intracellular ROS functions is still unclear. Here, for the first time, we determined the changes of biophysical properties associated with the ROS-mediated oesophageal cancer KYSE-150 cell apoptosis using high resolution atomic force microscopy (AFM). Oridonin was proved to induce ROS-mediated KYSE-150 cell apoptosis in a dose dependent manner, which could be reversed by N-acetylcysteine (NAC) pretreatment. Based on AFM imaging, the morphological damage and ultrastructural changes of KYSE-150 cells were found to be closely associated with ROS-mediated oridonin-induced KYSE-150 cell apoptosis. The changes of cell stiffness determined by AFM force measurement also demonstrated ROS-dependent changes in oridonin induced KYSE-150 cell apoptosis. Our findings not only provided new insights into the anticancer effects of oridonin, but also highlighted the use of AFM as a qualitative and quantitative nanotool to detect ROS-mediated cancer cell apoptosis based on cell biophysical properties, providing novel information of the roles of ROS in cancer cell apoptosis at nanoscale.

Carbon Nanotube Uptake Changes the Biomechanical Properties of Human Lung Epithelial Cells in a Time-dependent Manner

The toxicity of engineered nanomaterials in biological systems depends on both the nanomaterial properties and the exposure duration. Herein we used a multi-tier strategy to investigate the relationship between user-characterized multi-walled carbon nanotubes (MWCNTs) exposure duration and their induced biochemical and biomechanical effects on model human lung epithelial cells (BEAS-2B). Our results showed that exposure to MWCNTs leads to time-dependent intracellular uptake and generation of reactive oxygen species (ROS), along with time-dependent gradual changes in cellular biomechanical properties. In particular, the amount of internalized MWCNTs followed a sigmoidal curve with the majority of the MWCNTs being internalized within 6h of exposure; further, the sigmoidal uptake correlated with the changes in the oxidative levels and cellular biomechanical properties respectively. Our study provides new insights into the time-dependent induced toxicity caused by exposure to occupationally relevant doses of MWCNTs and could potentially help establish bases for early risk assessments of other nanomaterials toxicological profiles.

Mechanical coupling of cytoskeletal elasticity and force generation is crucial for understanding the migrating nature of keloid fibroblasts

One of the key features of keloid is its fibroblasts migrating beyond the original wound border. During migration, cells not only undergo molecular changes but also mechanical modulation. This process is led by actin filaments serving as the backbone of intra-cellular force and transduces external mechanical signal via focal adhesion complex into the cell. Here, we focus on determining the mechanical changes of actin filaments and the spatial distribution of forces in response to changing chemical stimulations and during cell migration. Atomic force microscopy and micropost array detector are used to determine and compare the magnitude and distribution of filament elasticity and force generation in fibroblasts and keloid fibroblasts. We found both filament elasticity and force generation show spatial distribution in a polarized and migrating cell. Such spatial distribution is disrupted when mechano-signalling is perturbed by focal adhesion kinase inhibitor and in keloid fibroblasts. The demonstration of keloid pathology at the nanoscale highlights the coupling of cytoskeletal function with physical characters at the subcellular level and provides new research directions for migration-related disease such as keloid.

Investigate correlation between mechanical property and aging biomarker in passaged human dermal fibroblasts

Skin aging is associated with changes in both the mechanical properties of the skin and extracellular matrix (ECM) components. In this study, we investigated the relationships between mechanical property and aging biomarkers in passaged human dermal fibroblasts (HDFs). The stiffness of HDFs from passages 5-20 was assessed by atomic force microscopy. The ECM components including collagen, elastin, and fibrillin-1 and that of signaling molecules (SIRTs) were determined from each passage of cells. The stiffness of HDFs increased linearly from passages 5-15 and then became saturated: the average stiffness was 0.356 N/m at passages 5 and 1.186 N/m at passages 15, respectively. Expression of all aging biomarkers, including pro-collagen I and VII, elastin, fibrillin-1, and SIRT1 and SIRT6, were down-regulated by passaging. In particular, a change in the level of procollagen Type I was significantly associated with early aging, while a change in the level of fibrillin-1 was associated with late aging. All biomarkers except elastin showed a strong correlation with the cellular stiffness of HDFs.